Skeletal muscle undergoes a progressive decline in mass and function with aging, a condition that in its extreme form is known as sarcopenia. This is driven by complex cellular and molecular alterations, such as shifts in myonucleus composition, increased fibrosis, and fat or immune cell infiltration. Despite extensive research, effective therapeutic interventions for sarcopenia remain limited. Recent advances in single-cell omics technologies have begun to unravel the cellular and molecular heterogeneity of mouse and human skeletal muscle across the lifespan, identifying age-enriched cell states and dynamic transcriptional changes. However, epigenetic regulation during skeletal muscle aging is less well characterized. To help address this gap, we performed single-nucleus Assay for Transposase-Accessible Chromatin using sequencing (snATAC-seq) on skeletal muscle from young adult and aged male mice, generating chromatin accessibility profiles from over 43,000 nuclei. Among other findings, our analyses reveal an age-enriched pro-atrophy subpopulation of type IIb myonuclei marked by increased chromatin accessibility at the Ampd3 locus. Furthermore, we delineate the epigenetic mechanisms underlying the transition of healthy type IIb myonuclei into Ampd3+ myonuclei, revealing key chromatin remodeling events that drive this phenotypic shift. Moreover, by integrating with an existing single-nucleus RNA sequencing dataset of the same anatomical origin, we identified thousands of cell-type-specific cis-regulatory elements related to aging programs. Within these elements, we observed a broad depletion of binding motifs for transcription factors with roles in cellular identity and muscle regeneration, concomitant with the gain of stress-responsive transcription factors. Our work helps understand the epigenetic events underlying mammalian skeletal muscle aging.

High-throughput single-cell omics of non-human primate brain tissue provides a powerful platform to investigate the molecular basis of brain aging. Here, we present a comprehensive transcriptomic and chromatin accessibility atlas of 2,955,873 nuclei from eight brain regions of 23 female cynomolgus macaques spanning the adult lifespan, including exceptionally old individuals. Our analyses reveal dynamic, cell-subtype- and region-specific age-related changes in core brain functions, including synaptic communication and axon myelination. We identify multicellular networks in the pons and medulla as a previously unrecognized hotspot of primate brain aging, highlighting white matter vulnerability as a central feature of aging. Integration with human brain aging and neurodegeneration datasets reveals both shared and divergent molecular mechanisms. We further define transcription factors and age-related chromatin remodeling programs linked to longevity and neurodegeneration. This spatiotemporal atlas establishes a foundational framework for understanding the cellular and regulatory architecture of primate brain aging and its links to disease.

The electron transport chain (ETC) is essential for NAD+ regeneration and proliferation. While many cell types tolerate ETC inhibition when pyruvate or aspartate is supplied, pluripotent stem cells (PSCs) enter a reversible paused state even at abundant pyruvate levels. Here, we show that ETC inhibition triggers severe NADH reductive stress in mouse embryonic stem cells (mESCs), driven mainly by threonine dehydrogenase (TDH). TDH-derived NADH establishes a metabolic environment that disfavors cells with compromised mitochondrial function, maintains inhibition of pyruvate dehydrogenase (PDH), and is associated with increased genomic and epigenetic stability at the cellular population level. ETC inhibition similarly induces pausing in early mouse embryos and in human pluripotent stem cells (hPSCs). In hPSCs, combined inhibition of the one-carbon metabolism enzymes serine hydroxymethyltransferase (SHMT1/2) and methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) effectively reduced reductive stress and rescued the paused phenotype. Together, these findings support a model in which limited mitochondrial redox adaptability represents a conserved metabolic feature of pluripotent stem cells and in which NADH reductive stress is associated with genomic and epigenetic stability.

The ability of RNA-binding proteins to form complexes with other biomolecules underpins a broad range of structural properties and functions. Understanding the subcellular distribution of RNA-binding proteins and their interacting partners in the steady state and upon perturbation can therefore shed light on these aspects. Here, we present the compartmentalized RNA-Binding Protein (or coRBP) map, an experimental resource and analytical pipeline to study subcellular RNA-binding proteins through multimodal dataset integration and machine learning. Using this approach, we generate a dataset of 1,768 known and putative RNA-binding proteins distributed in a broad panel of subcellular compartments and delineate their intermolecular and intercompartmental relationships. We also establish a hierarchy of RNA-binding protein-containing complexes at multiple scales across the cell, which suggests additional functions for multiple RNA-binding proteins. Furthermore, we investigate changes in RNA-binding protein complex composition and subcellular distribution in response to C9ORF72-associated amyotrophic lateral sclerosis/frontotemporal dementia dipeptide repeats and DNA damage stress. The coRBP map provides a resource to study the roles of RNA-binding proteins in homeostasis and disease.

BACKGROUND:Induced pluripotent stem cells induced hepatocytes (iHeps) are widely used in modeling human liver diseases and as a potential cell source for replacement therapy. However, most iHeps are relatively immature and challenging to maintain for long-term in vitro culture. METHODS:We optimized the differentiation protocol by addition of a combination of small molecules to inhibit epithelial-mesenchymal transition (EMT) in iHeps (iHeps EMTi), and further characterized their function both in vitro and in vivo analyses. RESULTS:Inhibition of EMT extended the in vitro culture period of iHeps EMTi from day 24 to day 60. In vitro analysis revealed that, compared to control, iHeps EMTi exhibited significantly higher expression levels of hepatic functional markers and enhanced hepatocyte functions, including lipid accumulation, glycogen storage, albumin secretion, and urea acid metabolism. Moreover, the molecular profiles of iHeps EMTi are closer to those of primary human hepatocytes. In addition, the in vivo engraftment efficiency of iHeps EMTi in the chimeric mice model was also improved as compared to iHeps alone. CONCLUSIONS:We established a robust protocol to generate human iHeps with improved function and capable of long-term in vitro culturing via the suppression of EMT. Moreover, those iHeps with EMT suppression have improved engraftment in human chimeric mice.

High-throughput single-cell omics of non-human primate tissues present a remarkable opportunity to study primate brain aging. Here, we introduce a transcriptomic and chromatin accessibility landscape of 1,985,317 cells from eight brain regions of 13 cynomolgus female monkeys spanning adult lifespan including exceptionally old individuals up to 29-years old. This dataset uncovers dynamic molecular changes in critical brain functions such as synaptic communication and axon myelination, exhibiting a high degree of cell type and brain region specificity. We identify the multicellular networks of the pons and medulla as a previously unrecognized hotspot for aging. Furthermore, comparative analyses with human neurodegeneration datasets highlight both shared and distinct mechanisms contributing to aging and disease. In addition, we uncover transcription factors implicated in monkey brain aging and pinpoint aging-regulated loci linked to longevity and neurodegeneration. This spatiotemporal atlas will advance our understanding of primate brain aging and its broader implications for health and disease. ### Competing Interest Statement The authors have declared no competing interest.

The mechanisms that ensure developmental progression in the early human embryo remain largely unknown. Here, we show that the family of long interspersed nuclear element 1 (LINE1) transposons prevents the reversion of naive human embryonic stem cells (hESCs) to 8-cell-like cells (8CLCs). LINE1 RNA contributes to maintenance of H3K27me3 levels, particularly at chromosome 19 (Chr19). Chr19 is enriched for key 8C regulators, H3K27me3, and genes derepressed upon LINE1 knockdown or PRC2 inhibition. Moreover, Chr19 is strongly associated with the nucleolus in hESCs but less in 8CLCs. Direct inhibition of PRC2 activity induces the 8C program and leads to a relocalization of Chr19 away from the nucleolus. LINE1 KD or PRC2 inhibition induces nucleolar stress, and disruption of nucleolar architecture is sufficient to de-repress the 8C program. These results indicate that LINE1 RNA and PRC2 maintain H3K27me3-mediated gene repression and 3D nuclear organization to prevent developmental reversion of hESCs.

Single-cell sequencing technologies have advanced our understanding of cellular heterogeneity and biological complexity. However, existing methods face limitations in throughput, capture uniformity, cell size flexibility, and technical extensibility. We present Stereo-cell, a spatial enhanced-resolution single-cell sequencing platform based on high-density DNA nanoball (DNB)-patterned arrays, which enables scalable and unbiased cell capture at a wide input range and supports high-fidelity transcriptome profiling. Stereo-cell further allows integration with imaging-based modalities and multiomics strategies, including immunofluorescence and epitope profiling. This platform is also compatible with profiling extracellular vesicles, microstructures, and large cells, whereas its spatial resolution facilitates in situ analysis of cell-cell interactions, cellular microenvironments, and subcellular transcript localization. Together, Stereo-cell provides a flexible framework for expanding single-cell research applications.

The human peripheral blood displays diverse molecular characteristics across populations, understanding the drivers and underlying mechanisms of which remains challenging. Here, we introduce the Chinese Immune Multi-Omics Atlas (CIMA), elucidating sex-, age-, and genetic-related molecular variations by analyzing multi-omics data from 428 adults with over 10 million immune cells. CIMA generated an enhancer-driven gene regulatory network, identifying 237 high-quality regulons and revealing cell type-specific regulatory mechanisms. Additionally, 11,521 lead cis-expression quantitative trait loci (eQTLs) and 46,339 chromatin accessibility QTLs (caQTLs) were identified at cell type level. CIMA also uncovered pleiotropic associations among immune-related disease risk loci, eQTLs, and caQTLs in a cell type-specific manner. Lastly, a novel cell language model, CIMA-CLM, was developed to predict chromatin accessibility and noncoding variant effects using chromatin sequences and gene expressions. This work represents a population-scale multi-omics resource of human immune cells, providing a valuable reference for future investigation of immune-related diseases. ### Competing Interest Statement The applications of this research are covered in pending patents. J.W. is founder of BGI-Shenzhen, and employee of BGI have stock holdings in BGI.

The mechanism by which mammalian liver cell responses are coordinated during tissue homeostasis and perturbation is poorly understood, representing a major obstacle in our understanding of many diseases. This knowledge gap is caused by the difficulty involved with studying multiple cell types in different states and locations, particularly when these are transient. We have combined Stereo-seq (spatiotemporal enhanced resolution omics-sequencing) with single-cell transcriptomic profiling of 473,290 cells to generate a high-definition spatiotemporal atlas of mouse liver homeostasis and regeneration at the whole-lobe scale. Our integrative study dissects in detail the molecular gradients controlling liver cell function, systematically defining how gene networks are dynamically modulated through intercellular communication to promote regeneration. Among other important regulators, we identified the transcriptional cofactor TBL1XR1 as a rheostat linking inflammation to Wnt/beta-catenin signaling for facilitating hepatocyte proliferation. Our data and analytical pipelines lay the foundation for future high-definition tissue-scale atlases of organ physiology and malfunction. A Stereo-seq and scRNA-seq atlas of mouse liver in homeostasis and regeneration after partial hepatectomy identifies zonated genes, pathways, cell-cell interactions and gene regulatory networks. Functional validation finds that cooperation between TBL1XR1 and beta-catenin activates hepatocyte proliferation.

RNA-binding proteins (RBPs) regulate totipotency, pluripotency maintenance, and induction. The intricacies of how they modulate these processes through their interaction with RNAs remain to be elucidated. Here we employed Targets of RBPs Identified By Editing (TRIBE) with single-cell resolution (scTRIBE) to profile the mRNA targets of the key pluripotency regulator LIN28A in mouse embryonic stem cells (ESCs), 2-cell embryo-like cells (2CLCs), and somatic cell reprogramming. LIN28A is known to act by controlling the maturation of the let-7 microRNA, but, in addition, it binds to multiple mRNAs and influences their stability and translation efficiency. However, the mRNA targets of LIN28A in 2CLCs and reprogramming are unclear. Through quantitative single-cell analysis of the scTRIBE dataset, we observed a marked increase in the binding of LIN28A to mRNAs of ribosome biogenesis factors and a selected group of totipotency factors in 2CLCs within ESC cultures. Our results suggest that LIN28A extends the half-life of at least some of these mRNAs, providing new insights into its role in the totipotent state. We also uncovered the distinct trajectory-specific LIN28A-mRNA networks in reprogramming, helping explain how LIN28A facilitates the mesenchymal-to-epithelial transition and pluripotency acquisition. Our study not only clarifies the multifunctional role of LIN28A in these processes but also highlights the importance of decoding RNA-protein interactions at the single-cell level.

Teratoma, due to its remarkable ability to differentiate into multiple cell lineages, is a valuable model for studying human embryonic development. The similarity of the gene expression and chromatin accessibility patterns in these cells to those observed in vivo further underscores its potential as a research tool. Notably, teratomas derived from human naïve (pre-implantation epiblast-like) pluripotent stem cells (PSCs) have larger embryonic cell diversity and contain extraembryonic lineages, making them more suitable to study developmental processes. However, the cell type-specific epigenetic profiles of naïve PSC teratomas have not been yet characterized. Using single-cell assay for transposase-accessible chromatin sequencing (scATAC-seq), we analyzed 66,384 cell profiles from five teratomas derived from human naïve PSCs and their post-implantation epiblast-like (primed) counterparts. We observed 17 distinct cell types from both embryonic and extraembryonic lineages, resembling the corresponding cell types in human fetal tissues. Additionally, we identified key transcription factors specific to different cell types. Our dataset provides a resource for investigating gene regulatory programs in a relevant model of human embryonic development.

In contrast to rodents, the mechanisms underlying human trophectoderm and early placenta specification are understudied due to ethical barriers and the scarcity of embryos. Recent reports have shown that human pluripotent stem cells (PSCs) can differentiate into trophectoderm (TE)-like cells (TELCs) and trophoblast stem cells (TSCs), offering a valuable in vitro model to study early placenta specification. Here, we demonstrate that the VGLL1 (vestigial-like family member 1), which is highly expressed during human and non-human primate TE specification in vivo but is negligibly expressed in mouse, is a critical regulator of cell fate determination and self-renewal in human TELCs and TSCs derived from naïve PSCs. Mechanistically, VGLL1 partners with the transcription factor TEAD4 (TEA domain transcription factor 4) to regulate chromatin accessibility at target gene loci through histone acetylation and acts in cooperation with GATA3 and TFAP2C. Our work is relevant to understand primate early embryogenesis and how it differs from other mammalian species.

The family of LINE1 transposable elements underwent a massive expansion in mammalian genomes. While traditionally viewed as a mutagenic selfish element, recent studies point to roles for LINE1 in early mouse development, T cell quiescence and neurogenesis. Here we show that human LINE1 RNA is essential for self-renewal and identity of human embryonic stem cells (hESCs). Silencing of LINE1 using either antisense oligonucleotides or CRISPR interference in naïve hESCs leads to a strong induction of 8C-like cells (8CLCs). We found that genes derepressed upon LINE1 KD are not uniformly distributed across the genome, with an enrichment for chromosome 19, which includes key markers of the 8C state such as TPRX1 . Silencing of TPRX1 , but not other putative 8C regulators p53 or H3.XY , suppresses the induction of the 8C program in LINE1 KD hESCs. We found that LINE1 RNA is preferentially localized to the lamina and periphery of the nucleolus in hESCs. Sequencing of Lamina-Associated Domains (LADs) and Nucleolus-Associated Domains (NADs) reveals a preferential association of chromosome 19 with NADs in hESCs. However, 8CLCs have a distinct nucleolar morphology and a lower association of chromosome 19 and TPRX1 loci with the nucleolus relative to naïve and primed hESCs, suggesting a role for nucleolar dynamics in the 8CLC-hESC transition. In agreement, LINE1 KD leads to disruption of nucleolar architecture with signs of nucleolar stress. Independent perturbations of the nucleolus induce the 8C program in hESCs. Genes induced by LINE1 KD are enriched for targets of Polycomb Repressive Complex (PRC2), and inhibition of PRC2 leads to a strong induction of 8C genes. Our results indicate that LINE1 coordinates nuclear compartmentalization and chromatin-mediated gene repression to prevent developmental reversion of hESCs. Highlights ### Competing Interest Statement The authors have declared no competing interest.

Muscle atrophy and functional decline (sarcopenia) are common manifestations of frailty and are critical contributors to morbidity and mortality in older people1. Deciphering the molecular mechanisms underlying sarcopenia has major implications for understanding human ageing2. Yet, progress has been slow, partly due to the difficulties of characterizing skeletal muscle niche heterogeneity (whereby myofibres are the most abundant) and obtaining well-characterized human samples3,4. Here we generate a single-cell/single-nucleus transcriptomic and chromatin accessibility map of human limb skeletal muscles encompassing over 387,000 cells/nuclei from individuals aged 15 to 99 years with distinct fitness and frailty levels. We describe how cell populations change during ageing, including the emergence of new populations in older people, and the cell-specific and multicellular network features (at the transcriptomic and epigenetic levels) associated with these changes. On the basis of cross-comparison with genetic data, we also identify key elements of chromatin architecture that mark susceptibility to sarcopenia. Our study provides a basis for identifying targets in the skeletal muscle that are amenable to medical, pharmacological and lifestyle interventions in late life.

Cholestatic liver injuries, characterized by regional damage around the bile ductular region, lack curative therapies and cause considerable mortality. Here we generated a high-definition spatiotemporal atlas of gene expression during cholestatic injury and repair in mice by integrating spatial enhanced resolution omics sequencing and single-cell transcriptomics. Spatiotemporal analyses revealed a key role of cholangiocyte-driven signaling correlating with the periportal damage-repair response. Cholangiocytes express genes related to recruitment and differentiation of lipid-associated macrophages, which generate feedback signals enhancing ductular reaction. Moreover, cholangiocytes express high TGF beta in association with the conversion of liver progenitor-like cells into cholangiocytes during injury and the dampened proliferation of periportal hepatocytes during recovery. Notably, Atoh8 restricts hepatocyte proliferation during 3,5-diethoxycarbonyl-1,4-dihydro-collidin damage and is quickly downregulated after injury withdrawal, allowing hepatocytes to respond to growth signals. Our findings lay a keystone for in-depth studies of cellular dynamics and molecular mechanisms of cholestatic injuries, which may further develop into therapies for cholangiopathies. A comprehensive time series characterization of a mouse model of cholestatic liver injury with spatial enhanced resolution omics sequencing and single-cell RNA sequencing identifies zonal responses to insult, such as cholangiocyte signaling recruiting lipid-associated macrophages.

The N6-methyladenosine (m6A) RNA modification plays essential roles in multiple biological processes, including stem cell fate determination. To explore the role of the m6A modification in pluripotent reprogramming, we used RNA-seq to map m6A effectors in human iPSCs, fibroblasts, and H9 ESCs, as well as in mouse ESCs and fibroblasts. By integrating the human and mouse RNA-seq data, we found that 19 m6A effectors were significantly upregulated in reprogramming. Notably, IGF2BPs, particularly IGF2BP1, were among the most upregulated genes in pluripotent cells, while YTHDF3 had high levels of expression in fibroblasts. Using quantitative PCR and Western blot, we validated the pluripotency-associated elevation of IGF2BPs. Knockdown of IGF2BP1 induced the downregulation of stemness genes and exit from pluripotency. Proteome analysis of cells collected at both the beginning and terminal states of the reprogramming process revealed that the IGF2BP1 protein was positively correlated with stemness markers SOX2 and OCT4. The eCLIP-seq target analysis showed that IGF2BP1 interacted with the coding sequence (CDS) and 3’UTR regions of the SOX2 transcripts, in agreement with the location of m6A modifications. This study identifies IGF2BP1 as a vital pluripotency-associated m6A effector, providing new insight into the interplay between m6A epigenetic modifications and pluripotent reprogramming.

A formal demonstration that mammalian pluripotent stem cells possess preimplantation embryonic cell-like (naive) pluripotency is the generation of chimeric animals through early embryo complementation with homologous cells. Whereas such naive pluripotency has been well demonstrated in rodents, poor chimerism has been achieved in other species including non-human primates due to the inability of the donor cells to match the developmental state of the host embryos. Here, we have systematically tested various culture conditions for establishing monkey naive embryonic stem cells and optimized the procedures for chimeric embryo culture. This approach generated an aborted fetus and a live chimeric monkey with high donor cell contribution. A stringent characterization pipeline demonstrated that donor cells efficiently (up to 90%) incorporated into various tissues (including the gonads and placenta) of the chimeric monkeys. Our results have major implications for the study of primate naive pluripotency and genetic engineering of non-human primates.