Single-cell analyses of X Chromosome inactivation dynamics and pluripotency during differentiation

The analysis of X chromosome activity of porcine embryonic stem Cells: Study based on parthenogenetic embryonic stem cells with LCDM medium

The derivation of porcine embryonic stem cell (pESC) lines remains a major challenge in this field. To date, the porcine naïve ESCs have yet to be successfully established, and standardized criteria for their characterization and evaluation are still lacking. The regulation of X-chromosome activity integrates information from embryonic development and the dosage of sex chromosomes, which is closely associated with the pluripotent state of embryonic stem cells. In this study, we aimed to establish pESC lines in LCDM medium from porcine blastocyst-stage embryos, and analyzed the features of ESCs from the sight of X chromosome activity. We assessed molecular markers and epigenetic characteristics to confirm pluripotency and X chromosome activity in porcine parthenogenetic ESCs (named as ppLCDM) using XIST RNA-FISH, immunofluorescence staining, single-cell RNA sequencing (scRNA-seq), and other techniques. Results showed that ppLCDM cells expressed most pluripotent markers. The percentage of ppLCDM cells exhibiting H3K27me3 and XIST aggregation signals increased with passage, indicating the progressive establishment of X-chromosome inactivation (XCI). Meanwhile, the pluripotency of most ppLCDM cells gradually declined during extended passaging. However, two distinct patterns of ppLCDM cells were observed from passage 35 (type I cells, P35-I) displayed normal XCI states, while type II cells (P35-II) exhibited X-chromosome erosion-like state, characterized by the loss of aggregation signals, abnormal X-linked gene ratios. Particularly, the pluripotency of ppLCDM cells with an X-chromosome erosion-like state undergoes unusual changes compared to normal cells. These findings indicate that X chromosome activity is closely associated with the pluripotent state of porcine ESCs and that heterogeneity in X chromosome activity arises during passaging. Our research provides crucial insights into X chromosome dynamics in large-animal ESC models and contribute to ongoing efforts to establish stable naïve pESC lines.

X-Linked Gene Dosage and SOX2 Act as Key Roadblocks for Human Germ Cell Specification in Klinefelter Syndrome

Klinefelter syndrome (KS), characterized by the presence of at least one extra X-chromosome, is a common cause of male infertility. However, the mechanism underlying the failure of germline specification is not well studied. Intriguingly, the differentiation efficiency of female human pluripotent stem cells (hPSCs) is often lower than that of male. This study investigates how X-linked gene dosage affects human primordial germ cell-like cells (hPGCLCs) specification in both healthy and diseased conditions. This work reveals that X-linked genes play a multifaceted role against the fate competency to hPGCLCs, with escape genes IGSF1 and CHRDL1 inhibiting the TGF-beta/Activin A and BMP pathways, respectively. Notably, this work identifies a previously unrecognized role of SOX2, upregulated by the escape gene USP9X, elucidating a species-specific function in the mammalian germline. The USP9X-SOX2 regulatory axis profoundly influenced cellular metabolism, mitochondrial morphology, and progenitor competence in hPGCLCs specification. Furthermore, the inability to downregulate SOX2 and upregulate SOX17 in response to BMP signaling impedes downstream gene activation due to motif binding competition. These findings shed novel insights into the human germline specification by elucidating the divergent roles of SOX2 versus SOX17 in mammals, influenced by X-linked gene dosage effects. These results offer potential applications for improving the induction efficiency of hPGCLCs, facilitating disease mechanistic studies.

Aging activates escape of the silent X chromosome in the female mouse hippocampus

Women live longer than men and exhibit less cognitive aging. The X chromosome contributes to sex differences, as females harbor an inactive X (Xi) and active X (Xa), in contrast to males with only an Xa. Thus, reactivation of silent Xi genes may contribute to sex differences. We use allele-specific, single-nucleus RNA sequencing to show that aging remodels transcription of the Xi and Xa across hippocampal cell types. Aging preferentially changed gene expression on the X's relative to autosomes. Select genes on the Xi underwent activation, with new escape across cells including in the dentate gyrus, critical to learning and memory. Expression of the Xi escapee Plp1, a myelin component, was increased in the aging hippocampus of female mice and parahippocampus of women. AAV-mediated Plp1 elevation in the dentate gyrus of aging male and female mice improved cognition. Understanding how the Xi may confer female advantage could lead to novel targets that counter brain aging and disease in both sexes.

Mouse totipotent blastomere-like cells model embryogenesis from zygotic genome activation to post implantation

Embryo development begins with zygotic genome activation (ZGA), eventually generating blastocysts for implantation. However, in vitro systems modeling the pre-implantation development are still absent and challenging. Here, we used mouse totipotent blastomere-like cells (TBLCs) to develop spontaneous differentiation and blastoid formation systems, respectively. We found Wnt signaling enabled the rapid expansion of TBLCs and the optimization of their culture medium. We successfully developed a TBLC-spontaneous differentiation system in which mouse TBLCs (mTBLCs) firstly converted into two types of ZGA-like cells (ZLCs) distinguished by Zscan4 expression. Surprisingly, Zscan4-, but not Zscan4+, ZLCs further passed through intermediate 4-cell and then 8-cell/morula stages to produce epiblast, primitive endoderm, and trophectoderm lineages. Significantly, single TBLCs underwent expansion, compaction, and polarization to efficiently generate blastocyst-like structures and even post-implantation egg-cylinder-like structures. Conclusively, we established TBLC-based differentiation and embryo-like structure formation systems to model early embryonic development, offering criteria for evaluating and understanding totipotency.

RNA splicing controls organ-wide maturation of postnatal heart in mice

Postnatal cardiac development requires the orchestrated maturation of diverse cellular components for which unifying control mechanisms are still lacking. Using full-length sequencing, we examined the transcriptomic landscape of the maturating mouse heart (E18.5-P28) at single-cell and transcript isoform resolution. We identified dynamically changing intercellular networks as a molecular basis of the maturing heart and alternative splicing (AS) as a common mechanism that distinguished developmental age. Manipulation of RNA-binding proteins (RBPs) remodeled the AS landscape, cardiac cell maturation, and intercellular communication through direct binding of splice targets, which were enriched for functions related to general, as well as cell-type-specific, maturation. Overexpression of an RBP nuclear cap-binding protein subunit 2 (NCBP2) in neonatal hearts repressed cardiac maturation. Together, our data suggest AS regulation by RBPs as an organ-level control mechanism in mammalian postnatal cardiac development and provide insight into the possibility of manipulating RBPs for therapeutic purposes.

Genetic gradual reduction of OGT activity unveils the essential role of O-GlcNAc in the mouse embryo

The reversible glycosylation of nuclear and cytoplasmic proteins (O-GlcNAcylation) is catalyzed by a single enzyme, namely O-GlcNAc transferase (OGT). The mammalian Ogt gene is X-linked, and it is essential for embryonic development and for the viability of proliferating cells. We perturbed OGT's function in vivo by creating a murine allelic series of four single amino acid substitutions, reducing OGT's catalytic activity to a range of degrees. The severity of the embryonic lethality was proportional to the extent of impairment of OGT's catalysis, demonstrating that the O-GlcNAc modification itself is required for early development. We identified hypomorphic Ogt alleles that perturb O-GlcNAc homeostasis while being compatible with embryogenesis. The analysis of the transcriptomes of the mutant embryos at different developmental stages suggested a sexually-dimorphic developmental delay caused by the decrease in O-GlcNAc. Furthermore, a mild reduction of OGT's enzymatic activity was sufficient to loosen the silencing of endogenous retroviruses in vivo.

Deep learning-based models for preimplantation mouse and human embryos based on single-cell RNA sequencing

The rapid growth of single-cell transcriptomic technology has produced an increasing number of datasets for both embryonic development and in vitro pluripotent stem cell-derived models. This avalanche of data surrounding pluripotency and the process of lineage specification has meant it has become increasingly difficult to define specific cell types or states in vivo, and compare these with in vitro differentiation. Here we utilize a set of deep learning tools to integrate and classify multiple datasets. This allows the definition of both mouse and human embryo cell types, lineages and states, thereby maximizing the information one can garner from these precious experimental resources. Our approaches are built on recent initiatives for large-scale human organ atlases, but here we focus on material that is difficult to obtain and process, spanning early mouse and human development. Using publicly available data for these stages, we test different deep learning approaches and develop a model to classify cell types in an unbiased fashion at the same time as defining the set of genes used by the model to identify lineages, cell types and states. We used our models trained on in vivo development to classify pluripotent stem cell models for both mouse and human development, showcasing the importance of this resource as a dynamic reference for early embryogenesis.

Dosage compensation of transposable elements in mammals

In mammals, X-linked dosage compensation involves two processes: X-chromosome inactivation (XCI) to balance X chromosome dosage between males and females, and hyperactivation of the remaining X chromosome (Xa-hyperactivation) to achieve X-autosome balance in both sexes. Studies of both processes have largely focused on coding genes and have not accounted for transposable elements (TEs) which comprise 50% of the X-chromosome, despite TEs being suspected to have numerous epigenetic functions. This oversight is due in part to the technical challenge of capturing repeat RNAs, bioinformatically aligning them, and determining allelic origin. To overcome these challenges, here we develop a new bioinformatic pipeline tailored to repetitive elements with capability for allelic discrimination. We then apply the pipeline to our recent So-Smart-Seq analysis of single embryos to comprehensively interrogate whether X-linked TEs are subject to either XCI or Xa-hyperactivation. With regards to XCI, we observe significant differences in TE silencing in parentally driven "imprinted" XCI versus zygotically driven "random" XCI. Chromosomal positioning and genetic background impact TE silencing. We also find that SINEs may influence 3D organization during XCI. In contrast, TEs do not undergo Xa-hyperactivation. Thus, while coding genes are subject to both forms of dosage compensation, TEs participate only in Xi silencing. Evolutionary and functional implications are discussed.

Lineage-specific dynamics of loss of X upregulation during inactive-X reactivation

In mammals, X chromosome dosage is balanced between sexes through the silencing of one X chromosome in females. Recent single-cell RNA sequencing analysis demonstrated that the inactivation of the X chromosome is accompanied by the upregulation of the active X chromosome (Xa) during mouse embryogenesis. Here, we have investigated if the reactivation of inactive-X (Xi) leads to the loss of Xa upregulation in different cellular or developmental contexts. We find that while Xi reactivation and loss of Xa upregulation are tightly coupled in mouse embryonic epiblast and induced pluripotent stem cells, that is not the case in germ cells. Moreover, we demonstrate that partial reactivation of Xi in mouse extra-embryonic endoderm stem cells and human B cells does not result in the loss of Xa upregulation. Finally, we have established a mathematical model for the transcriptional coordination of two X chromosomes. Together, we conclude that the reactivation of Xi is not always synchronized with the loss of Xa upregulation.

saseR: Juggling offsets unlocks RNA-seq tools for fast and Scalable differential usage, Aberrant Splicing and Expression Retrieval

RNA-seq data analysis relies on many different tools, each tailored to specific applications and coming with unique assumptions and restrictions. Indeed, tools for differential transcript usage, or diagnosing patients with rare diseases through splicing and expression outliers, either lack in performance, discard information, or do not scale to massive data compendia. Here, we show that replacing the normalisation offsets unlocks bulk RNA-seq workflows for scalable differential usage, aberrant splicing and expression analyses. Our method, saseR, is much faster than state-of-the-art methods, dramatically outperforms these to detect aberrant splicing, and provides a single workflow for various short- and long-read RNA-seq applications.

Modeling X chromosome inactivation using t5iLA naive human pluripotent stem cells

BACKGROUND

X chromosome inactivation (XCI) is a critical epigenetic event for dosage compensation of X-linked genes in female mammals, ensuring developmental stability. A robust in vitro model is required for mimicking XCI during the early stages of embryonic development. This methodology article introduces an advanced framework for the in-depth study of XCI using human pluripotent stem cells (hPSCs). By focusing on the transition between naive and primed pluripotent states, we highlight the role of long non-coding RNA X-inactive specific transcript (XIST) and epigenetic alterations in mediating XCI.

RESULTS

Our methodology enables the distinction between naive and primed hESCs based on XIST expression and the activity of X-linked reporters, facilitating the investigation of XCI initiation and maintenance. Through detailed experimental procedures, we demonstrate the utility of our hESC lines in modeling the process of human XCI, including the establishment of conditions for random XCI induction and the analysis of X chromosome reactivation.

METHODS

The study outlines a comprehensive approach for characterizing the X chromosome status in hPSCs, employing dual fluorescent reporter hESC lines. These reporter lines enable real-time tracking of XCI dynamics through differentiation processes. We detailed protocols for the induction of X chromosome reactivation and inactivation, as well as the X status characterization methods including cultivation of hESCs, flow cytometric analysis, RNA fluorescence in situ hybridization (FISH), and transcriptome sequencing, providing a step-by-step guide for researchers to investigate XCI mechanisms in vitro.

CONCLUSIONS

This article provides a detailed, reproducible methodology for studying XCI mechanisms in vitro, employing hPSCs as a model system. It presents a significant advance in our ability to investigate XCI, offering potential applications in developmental biology, disease modeling, and regenerative medicine. By facilitating the study of XCI dynamics, this methodological framework paves the way for deeper understanding and manipulation of this fundamental biological process.

A system to analyze the initiation of random X-chromosome inactivation using time-lapse imaging of single cells

In female eutherian mammal development, X-chromosome inactivation (XCI) of one of the two X chromosomes is initiated early. Understanding the relationship between the initiation of XCI and cell fate is critical for understanding early female development and requires a system that can monitor XCI in single living cells. Traditional embryonic stem cells (ESCs) used for XCI studies often lose X chromosomes spontaneously during culture and differentiation, making accurate monitoring difficult. Additionally, most XCI assessment methods necessitate cell disruption, hindering cell fate tracking. We developed the Momiji (version 2) ESC line to address these difficulties, enabling real-time monitoring of X-chromosome activity via fluorescence. We inserted green and red fluorescent reporter genes and neomycin and puromycin resistance genes into the two X chromosomes of PGK12.1 ESCs, creating a female ESC line that retains two X chromosomes more faithfully during differentiation. Momiji (version 2) ESCs exhibit a more stable XX karyotype than other ESC lines, including the parental PGK12.1 line. This new tool offers valuable insights into the relationship between XCI and cell fate, improving our understanding of early female development.

The effect of data transformation on low-dimensional integration of single-cell RNA-seq

BACKGROUND

Recent developments in single-cell RNA sequencing have opened up a multitude of possibilities to study tissues at the level of cellular populations. However, the heterogeneity in single-cell sequencing data necessitates appropriate procedures to adjust for technological limitations and various sources of noise when integrating datasets from different studies. While many analysis procedures employ various preprocessing steps, they often overlook the importance of selecting and optimizing the employed data transformation methods.

RESULTS

This work investigates data transformation approaches used in single-cell clustering analysis tools and their effects on batch integration analysis. In particular, we compare 16 transformations and their impact on the low-dimensional representations, aiming to reduce the batch effect and integrate multiple single-cell sequencing data. Our results show that data transformations strongly influence the results of single-cell clustering on low-dimensional data space, such as those generated by UMAP or PCA. Moreover, these changes in low-dimensional space significantly affect trajectory analysis using multiple datasets, as well. However, the performance of the data transformations greatly varies across datasets, and the optimal method was different for each dataset. Additionally, we explored how data transformation impacts the analysis of deep feature encodings using deep neural network-based models, including autoencoder-based models and proto-typical networks. Data transformation also strongly affects the outcome of deep neural network models.

CONCLUSIONS

Our findings suggest that the batch effect and noise in integrative analysis are highly influenced by data transformation. Low-dimensional features can integrate different batches well when proper data transformation is applied. Furthermore, we found that the batch mixing score on low-dimensional space can guide the selection of the optimal data transformation. In conclusion, data preprocessing is one of the most crucial analysis steps and needs to be cautiously considered in the integrative analysis of multiple scRNA-seq datasets.

Single-Cell Transcriptome Sequence Profiling on the Morphogenesis of Secondary Hair Follicles in Ordos Fine-Wool Sheep

The Ordos fine-wool sheep is a high-quality breed in China that produces superior natural textiles and raw materials such as wool and lamb meat. However, compared to the Australian Merino sheep, there is still a gap in terms of the wool fiber fineness and wool yield. The hair follicle is the main organ that controls the type of wool fiber, and the morphological changes in the secondary hair follicle are crucial in determining wool quality. However, the process and molecular mechanisms of hair follicle morphogenesis in Ordos fine-wool sheep are not yet clear. Therefore, analyzing the molecular mechanisms underlying the process of follicle formation is of great significance for improving the fiber diameter and wool production of Ordos fine-wool sheep. The differential expressed genes, APOD, POSTN, KRT5, and KRT15, which related to primary hair follicles and secondary hair follicles, were extracted from the dermal papillae. Based on pseudo-time analysis, the differentiation trajectories of dermal lineage cells and epidermal lineage cells in the Ordos fine-wool sheep were successfully constructed, providing a theoretical basis for breeding research in Ordos fine-wool sheep.

Manual Dissociation of Mammalian Preimplantation Embryos for Single-Cell Genomics

Single-cell genomics allow the characterization and quantification of molecular heterogeneity from a wide variety of tissues. Here, we describe the manual dissociation and collection of single cells, a method adapted for the characterization of precious small tissues like preimplantation embryos. We also describe the acquisition of mouse embryos by flushing of the oviducts. The cells can then be used in multiple sequencing protocols, for example, Smart-seq2, Smart-seq3, smallseq, and scBSseq.

Chromatin-mediated silencing on the inactive X chromosome

In mammals, the second X chromosome in females is silenced to enable dosage compensation between XX females and XY males. This essential process involves the formation of a dense chromatin state on the inactive X (Xi) chromosome. There is a wealth of information about the hallmarks of Xi chromatin and the contribution each makes to silencing, leaving the tantalising possibility of learning from this knowledge to potentially remove silencing to treat X-linked diseases in females. Here, we discuss the role of each chromatin feature in the establishment and maintenance of the silent state, which is of crucial relevance for such a goal.

Single-cell analysis of human prepuce reveals dynamic changes in gene regulation and cellular communications

BACKGROUND

The cellular and molecular dynamics of human prepuce are crucial for understanding its biological and physiological functions, as well as the prevention of related genital diseases. However, the cellular compositions and heterogeneity of human prepuce at single-cell resolution are still largely unknown. Here we systematically dissected the prepuce of children and adults based on the single-cell RNA-seq data of 90,770 qualified cells.

RESULTS

We identified 15 prepuce cell subtypes, including fibroblast, smooth muscle cells, T/natural killer cells, macrophages, vascular endothelial cells, and dendritic cells. The proportions of these cell types varied among different individuals as well as between children and adults. Moreover, we detected cell-type-specific gene regulatory networks (GRNs), which could contribute to the unique functions of related cell types. The GRNs were also highly dynamic between the prepuce cells of children and adults. Our cell-cell communication network analysis among different cell types revealed a set of child-specific (e.g., CD96, EPO, IFN-1, and WNT signaling pathways) and adult-specific (e.g., BMP10, NEGR, ncWNT, and NPR1 signaling pathways) signaling pathways. The variations of GRNs and cellular communications could be closely associated with prepuce development in children and prepuce maintenance in adults.

CONCLUSIONS

Collectively, we systematically analyzed the cellular variations and molecular changes of the human prepuce at single-cell resolution. Our results gained insights into the heterogeneity of prepuce cells and shed light on the underlying molecular mechanisms of prepuce development and maintenance.

RBBP4 is an epigenetic barrier for the induced transition of pluripotent stem cells into totipotent 2C-like cells

Cellular totipotency is critical for whole-organism generation, yet how totipotency is established remains poorly illustrated. Abundant transposable elements (TEs) are activated in totipotent cells, which is critical for embryonic totipotency. Here, we show that the histone chaperone RBBP4, but not its homolog RBBP7, is indispensable for maintaining the identity of mouse embryonic stem cells (mESCs). Auxin-induced degradation of RBBP4, but not RBBP7, reprograms mESCs to the totipotent 2C-like cells. Also, loss of RBBP4 enhances transition from mESCs to trophoblast cells. Mechanistically, RBBP4 binds to the endogenous retroviruses (ERVs) and functions as an upstream regulator by recruiting G9a to deposit H3K9me2 on ERVL elements, and recruiting KAP1 to deposit H3K9me3 on ERV1/ERVK elements, respectively. Moreover, RBBP4 facilitates the maintenance of nucleosome occupancy at the ERVK and ERVL sites within heterochromatin regions through the chromatin remodeler CHD4. RBBP4 depletion leads to the loss of the heterochromatin marks and activation of TEs and 2C genes. Together, our findings illustrate that RBBP4 is required for heterochromatin assembly and is a critical barrier for inducing cell fate transition from pluripotency to totipotency.

Comparative developmental genomics of sex-biased gene expression in early embryogenesis across mammals

BACKGROUND

Mammalian gonadal sex is determined by the presence or absence of a Y chromosome and the subsequent production of sex hormones contributes to secondary sexual differentiation. However, sex chromosome-linked genes encoding dosage-sensitive transcription and epigenetic factors are expressed well before gonad formation and have the potential to establish sex-biased expression that persists beyond the appearance of gonadal hormones. Here, we apply a comparative bioinformatics analysis on a pair of published single-cell datasets from mouse and human during very early embryogenesis-from two-cell to pre-implantation stages-to characterize sex-specific signals and to assess the degree of conservation among early acting sex-specific genes and pathways.

RESULTS

Clustering and regression analyses of gene expression across samples reveal that sex initially plays a significant role in overall gene expression patterns at the earliest stages of embryogenesis which potentially may be the byproduct of signals from male and female gametes during fertilization. Although these transcriptional sex effects rapidly diminish, sex-biased genes appear to form sex-specific protein-protein interaction networks across pre-implantation stages in both mammals providing evidence that sex-biased expression of epigenetic enzymes may establish sex-specific patterns that persist beyond pre-implantation. Non-negative matrix factorization (NMF) on male and female transcriptomes generated clusters of genes with similar expression patterns across sex and developmental stages, including post-fertilization, epigenetic, and pre-implantation ontologies conserved between mouse and human. While the fraction of sex-differentially expressed genes (sexDEGs) in early embryonic stages is similar and functional ontologies are conserved, the genes involved are generally different in mouse and human.

CONCLUSIONS

This comparative study uncovers much earlier than expected sex-specific signals in mouse and human embryos that pre-date hormonal signaling from the gonads. These early signals are diverged with respect to orthologs yet conserved in terms of function with important implications in the use of genetic models for sex-specific disease.

Cross-activation of FGF, NODAL, and WNT pathways constrains BMP-signaling-mediated induction of the totipotent state in mouse embryonic stem cells

Embryonic stem cells (ESCs) are an attractive model to study the relationship between signaling and cell fates. Cultured mouse ESCs can exist in multiple states resembling distinct stages of early embryogenesis, such as totipotent, pluripotent, primed, and primitive endoderm. The signaling mechanisms regulating the totipotent state and coexistence of these states are poorly understood. Here we identify bone morphogenetic protein (BMP) signaling as an inducer of the totipotent state. However, we discover that BMP's role is constrained by the cross-activation of FGF, NODAL, and WNT pathways. We exploit this finding to enhance the proportion of totipotent cells by rationally inhibiting the cross-activated pathways. Single-cell mRNA sequencing reveals that induction of the totipotent state is accompanied by suppression of primed and primitive endoderm states. Furthermore, reprogrammed totipotent cells we generate in culture resemble totipotent cells of preimplantation embryo. Our findings reveal a BMP signaling mechanism regulating both the totipotent state and heterogeneity of ESCs.

Decoding anterior-posterior axis emergence among mouse, monkey, and human embryos

Anterior-posterior axis formation regulated by the distal visceral endoderm (DVE) and anterior visceral endoderm (AVE) is essential for peri-implantation embryogenesis. However, the principles of the origin and specialization of DVE and AVE remain elusive. Here, with single-cell transcriptome analysis and pseudotime prediction, we show that DVE and AVE independently originate from the specialized primary endoderm in mouse blastocysts. Along distinct developmental paths, these two lineages, respectively, undergo four representative states with stage-specific transcriptional patterns around implantation. Further comparative analysis shows that AVE, but not DVE, is detected in human and non-human primate embryos, defining differences in polarity formation across species. Moreover, stem cell-assembled human blastoids lack DVE or AVE precursors, implying that additional induction of stem cells with DVE/AVE potential could promote the current embryo-like models and their post-implantation growth. Our work provides insight into understanding of embryonic polarity formation and early mammalian development.

DNA methyltransferases 3A and 3B target specific sequences during mouse gastrulation

In mammalian embryos, DNA methylation is initialized to maximum levels in the epiblast by the de novo DNA methyltransferases DNMT3A and DNMT3B before gastrulation diversifies it across regulatory regions. Here we show that DNMT3A and DNMT3B are differentially regulated during endoderm and mesoderm bifurcation and study the implications in vivo and in meso-endoderm embryoid bodies. Loss of both Dnmt3a and Dnmt3b impairs exit from the epiblast state. More subtly, independent loss of Dnmt3a or Dnmt3b leads to small biases in mesoderm-endoderm bifurcation and transcriptional deregulation. Epigenetically, DNMT3A and DNMT3B drive distinct methylation kinetics in the epiblast, as can be predicted from their strand-specific sequence preferences. The enzymes compensate for each other in the epiblast, but can later facilitate lineage-specific methylation kinetics as their expression diverges. Single-cell analysis shows that differential activity of DNMT3A and DNMT3B combines with replication-linked methylation turnover to increase epigenetic plasticity in gastrulation. Together, these findings outline a dynamic model for the use of DNMT3A and DNMT3B sequence specificity during gastrulation.

INSISTC: Incorporating network structure information for single-cell type classification

Uncovering gene regulatory mechanisms in individual cells can provide insight into cell heterogeneity and function. Recent accumulated Single-Cell RNA-Seq data have made it possible to analyze gene regulation at single-cell resolution. Understanding cell-type-specific gene regulation can assist in more accurate cell type and state identification. Computational approaches utilizing such relationships are under development. Methods pioneering in integrating gene regulatory mechanism discovery with cell-type classification encounter challenges such as determine gene regulatory relationships and incorporate gene regulatory network structure. To fill this gap, we developed INSISTC, a computational method to incorporate gene regulatory network structure information for single-cell type classification. INSISTC is capable of identifying cell-type-specific gene regulatory mechanisms while performing single-cell type classification. INSISTC demonstrated its accuracy in cell type classification and its potential for providing insight into molecular mechanisms specific to individual cells. In comparison with the alternative methods, INSISTC demonstrated its complementary performance for gene regulation interpretation.

In Vivo Clonal Analysis Reveals Random Monoallelic Expression in Lymphocytes That Traces Back to Hematopoietic Stem Cells

Evaluating the epigenetic landscape in the stem cell compartment at the single-cell level is essential to assess the cells’ heterogeneity and predict their fate. Here, using a genome-wide transcriptomics approach in vivo, we evaluated the allelic expression imbalance in the progeny of single hematopoietic cells (HSCs) as a read-out of epigenetic marking. After 4 months of extensive proliferation and differentiation, we found that X-chromosome inactivation (XCI) is tightly maintained in all single-HSC derived hematopoietic cells. In contrast, the vast majority of the autosomal genes did not show clonal patterns of random monoallelic expression (RME). However, a persistent allele-specific autosomal transcription in HSCs and their progeny was found in a rare number of cases, none of which has been previously reported. These data show that: 1) XCI and RME in the autosomal chromosomes are driven by different mechanisms; 2) the previously reported high frequency of genes under RME in clones expanded in vitro (up to 15%) is not found in clones undergoing multiple differentiation steps in vivo; 3) prior to differentiation, HSCs have stable patterns of autosomal RME. We propose that most RME patterns in autosomal chromosomes are erased and established de novo during cell lineage differentiation.

SmcHD1 underlies the formation of H3K9me3 blocks on the inactive X chromosome in mice

Stable silencing of the inactive X chromosome (Xi) in female mammals is crucial for the development of embryos and their postnatal health. SmcHD1 is essential for stable silencing of the Xi, and its functional deficiency results in derepression of many X-inactivated genes. Although SmcHD1 has been suggested to play an important role in the formation of higher-order chromatin structure of the Xi, the underlying mechanism is largely unknown. Here, we explore the epigenetic state of the Xi in SmcHD1-deficient epiblast stem cells and mouse embryonic fibroblasts in comparison with their wild-type counterparts. The results suggest that SmcHD1 underlies the formation of H3K9me3-enriched blocks on the Xi, which, although the importance of H3K9me3 has been largely overlooked in mice, play a crucial role in the establishment of the stably silenced state. We propose that the H3K9me3 blocks formed on the Xi facilitate robust heterochromatin formation in combination with H3K27me3, and that the substantial loss of H3K9me3 caused by SmcHD1 deficiency leads to aberrant distribution of H3K27me3 on the Xi and derepression of X-inactivated genes.

Identifying cell state-associated alternative splicing events and their coregulation

Alternative splicing shapes the transcriptome and contributes to each cell's unique identity, but single-cell RNA sequencing (scRNA-seq) has struggled to capture the impact of alternative splicing. We previously showed that low recovery of mRNAs from single cells led to erroneous conclusions about the cell-to-cell variability of alternative splicing. Here, we present a method, Psix, to confidently identify splicing that changes across a landscape of single cells, using a probabilistic model that is robust against the data limitations of scRNA-seq. Its autocorrelation-inspired approach finds patterns of alternative splicing that correspond to patterns of cell identity, such as cell type or developmental stage, without the need for explicit cell clustering, labeling, or trajectory inference. Applying Psix to data that follow the trajectory of mouse brain development, we identify exons whose alternative splicing patterns cluster into modules of coregulation. We show that the exons in these modules are enriched for binding by distinct neuronal splicing factors and that their changes in splicing correspond to changes in expression of these splicing factors. Thus, Psix reveals cell type-dependent splicing patterns and the wiring of the splicing regulatory networks that control them. Our new method will enable scRNA-seq analysis to go beyond transcription to understand the roles of post-transcriptional regulation in determining cell identity.

Identification of Five Hub Genes Based on Single-Cell RNA Sequencing Data and Network Pharmacology in Patients With Acute Myocardial Infarction

Acute myocardial infarction (AMI) has a high mortality. The single-cell RNA sequencing (scRNA-seq) method was used to analyze disease heterogeneity at the single-cell level. From the Gene Expression Omnibus (GEO) database (GSE180678), AMI scRNA-seq were downloaded and preprocessed by the Seurat package. Gene expression data came from GSE182923. Cell cluster analysis was conducted. Cell types were identified. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) analyses were performed on hub genes. Drugs were predicted by protein–protein interaction (PPI) and molecular docking. In total, 7 cell clusters were defined based on the scRNA-seq dataset, and the clusters were labeled as 5 cell types by marker genes. Hematopoietic stem cell types as a differential subgroups were higher in AMI than in healthy tissues. From available databases and PPI analysis, 52 common genets were identified. Based on 52 genes, 5 clusters were obtained using the MCODE algorithm, and genes in these 5 clusters involved in immune and inflammatory pathways were determined. Correlation analysis showed that hematopoietic stem cell types were negatively correlated with ATM, CARM1, and CASP8 but positively correlated with CASP3 and PPARG. This was reversed with immune cells. Molecular docking analysis showed that DB05490 had the lowest docking score with PPARG. We identified 5 hub genes (ATM, CARM1, CASP8, CASP3, and PPARG) involved in AMI progression. Compound DB05490 was a potential inhibitor of PPAG.

A hexa-species transcriptome atlas of mammalian embryogenesis delineates metabolic regulation across three different implantation modes

Mammalian embryogenesis relies on glycolysis and oxidative phosphorylation to balance the generation of biomass with energy production. However, the dynamics of metabolic regulation in the postimplantation embryo in vivo have remained elusive due to the inaccessibility of the implanted conceptus for biochemical studies. To address this issue, we compiled single-cell embryo profiling data in six mammalian species and determined their metabolic dynamics through glycolysis and oxidative phosphorylation associated gene expression. Strikingly, we identify a conserved switch from bivalent respiration in the late blastocyst towards a glycolytic metabolism in early gastrulation stages across species, which is independent of embryo implantation. Extraembryonic lineages followed the dynamics of the embryonic lineage, except visceral endoderm. Finally, we demonstrate that in vitro primate embryo culture substantially impacts metabolic gene regulation by comparison to in vivo samples. Our work reveals a conserved metabolic programme despite different implantation modes and highlights the need to optimise postimplantation embryo culture protocols.

Preventing erosion of X-chromosome inactivation in human embryonic stem cells

X-chromosome inactivation is a paradigm of epigenetic transcriptional regulation. Female human embryonic stem cells (hESCs) often undergo erosion of X-inactivation upon prolonged culture. Here, we investigate the sources of X-inactivation instability by deriving new primed pluripotent hESC lines. We find that culture media composition dramatically influenced the expression of XIST lncRNA, a key regulator of X-inactivation. hESCs cultured in a defined xenofree medium stably maintained XIST RNA expression and coating, whereas hESCs cultured in the widely used mTeSR1 medium lost XIST RNA expression. We pinpointed lithium chloride in mTeSR1 as a cause of XIST RNA loss. The addition of lithium chloride or inhibitors of GSK-3 proteins that are targeted by lithium to the defined hESC culture medium impeded XIST RNA expression. GSK-3 inhibition in differentiating female mouse embryonic stem cells and epiblast stem cells also resulted in a loss of XIST RNA expression. Together, these data may reconcile observed variations in X-inactivation in hESCs and inform the faithful culture of pluripotent stem cells.

Derivation of totipotent-like stem cells with blastocyst-like structure forming potential

It is challenging to derive totipotent stem cells in vitro that functionally and molecularly resemble cells from totipotent embryos. Here, we report that a chemical cocktail enables the derivation of totipotent-like stem cells, designated as totipotent potential stem (TPS) cells, from 2-cell mouse embryos and extended pluripotent stem cells, and that these TPS cells can be stably maintained long term in vitro. TPS cells shared features with 2-cell mouse embryos in terms of totipotency markers, transcriptome, chromatin accessibility and DNA methylation patterns. In vivo chimera formation assays show that these cells have embryonic and extraembryonic developmental potentials at the single-cell level. Moreover, TPS cells can be induced into blastocyst-like structures resembling preimplantation mouse blastocysts. Mechanistically, inhibition of HDAC1/2 and DOT1L activity and activation of RARγ signaling are important for inducing and maintaining totipotent features of TPS cells. Our study opens up a new path toward fully capturing totipotent stem cells in vitro.

Elastic dosage compensation by X-chromosome upregulation

X-chromosome inactivation and X-upregulation are the fundamental modes of chromosome-wide gene regulation that collectively achieve dosage compensation in mammals, but the regulatory link between the two remains elusive and the X-upregulation dynamics are unknown. Here, we use allele-resolved single-cell RNA-seq combined with chromatin accessibility profiling and finely dissect their separate effects on RNA levels during mouse development. Surprisingly, we uncover that X-upregulation elastically tunes expression dosage in a sex- and lineage-specific manner, and moreover along varying degrees of X-inactivation progression. Male blastomeres achieve X-upregulation upon zygotic genome activation while females experience two distinct waves of upregulation, upon imprinted and random X-inactivation; and ablation of Xist impedes female X-upregulation. Female cells carrying two active X chromosomes lack upregulation, yet their collective RNA output exceeds that of a single hyperactive allele. Importantly, this conflicts the conventional dosage compensation model in which naïve female cells are initially subject to biallelic X-upregulation followed by X-inactivation of one allele to correct the X dosage. Together, our study provides key insights to the chain of events of dosage compensation, explaining how transcript copy numbers can remain remarkably stable across developmental windows wherein severe dose imbalance would otherwise be experienced by the cell.

In vitro investigation of mammalian peri-implantation embryogenesis†

The embryos attach and invade into the uterus and establish the connection with their mother in peri-implantation development. During this period, the pluripotent epiblast cells of embryo undergo symmetry breaking, cell lineage allocation, and morphogenetic remodeling, accompanying with the dramatic changes of transcriptome, epigenome, and signal pathways, to prepare a state for their differentiation and gastrulation. The progresses in mouse genetics and stem cell biology have largely advanced the knowledge of these transformations which are largely hindered by the hard accessibility of natural embryos. To gain insight into mammalian peri-implantation development, great efforts have been made in the field. Recently, the advances in the prolonged in vitro culture of blastocysts, the derivation of multiple pluripotent stem cells, as well as the construction of stem cell-based embryo-like models have opened novel avenues to investigate peri-implantation development in mammals, especially for the humans. Combining with other emerging new technologies, these new models will substantially promote the comprehension of mammalian peri-implantation development, accelerating the progress of reproductive and regenerative medicine.

Single-Cell Regulatory Network Inference and Clustering Identifies Cell-Type Specific Expression Pattern of Transcription Factors in Mouse Sciatic Nerve

Advances in single-cell RNA sequencing technologies and bioinformatics methods allow for both the identification of cell types in a complex tissue and the large-scale gene expression profiling of various cell types in a mixture. In this report, we analyzed a single-cell RNA sequencing (scRNA-seq) dataset for the intact adult mouse sciatic nerve and examined cell-type specific transcription factor expression and activity during peripheral nerve homeostasis. In total, we identified 238 transcription factors expressed in nine different cell types of intact mouse sciatic nerve. Vascular smooth muscle cells have the lowest number of transcription factors expressed with 17 transcription factors identified. Myelinating Schwann cells (mSCs) have the highest number of transcription factors expressed, with 61 transcription factors identified. We created a cell-type specific expression map for the identified 238 transcription factors. Our results not only provide valuable information about the expression pattern of transcription factors in different cell types of adult peripheral nerves but also facilitate future studies to understand the function of key transcription factors in the peripheral nerve homeostasis and disease.

Induction of Rosette-to-Lumen stage embryoids using reprogramming paradigms in ESCs

Blastocyst-derived stem cell lines were shown to self-organize into embryo-like structures in 3D cell culture environments. Here, we provide evidence that embryo-like structures can be generated solely based on transcription factor-mediated reprogramming of embryonic stem cells in a simple 3D co-culture system. Embryonic stem cells in these cultures self-organize into elongated, compartmentalized embryo-like structures reflecting aspects of the inner regions of the early post-implantation embryo. Single-cell RNA-sequencing reveals transcriptional profiles resembling epiblast, primitive-/visceral endoderm, and extraembryonic ectoderm of early murine embryos around E4.5-E5.5. In this stem cell-based embryo model, progression from rosette formation to lumenogenesis accompanied by progression from naïve- to primed pluripotency was observed within Epi-like cells. Additionally, lineage specification of primordial germ cells and distal/anterior visceral endoderm-like cells was observed in epiblast- or visceral endoderm-like compartments, respectively. The system presented in this study allows for fast and reproducible generation of embryo-like structures, providing an additional tool to study aspects of early embryogenesis.

Single-Cell Genomics: Enabling the Functional Elucidation of Infectious Diseases in Multi-Cell Genomes

Since the time when detection of gene expression in single cells by microarrays to the Next Generation Sequencing (NGS) enabled Single Cell Genomics (SCG), it has played a pivotal role to understand and elucidate the functional role of cellular heterogeneity. Along this journey to becoming a key player in the capture of the individuality of cells, SCG overcame many milestones, including scale, speed, sensitivity and sample costs (4S). There have been many important experimental and computational innovations in the efficient analysis and interpretation of SCG data. The increasing role of AI in SCG data analysis has further enhanced its applicability in building models for clinical intervention. Furthermore, SCG has been instrumental in the delineation of the role of cellular heterogeneity in specific diseases, including cancer and infectious diseases. The understanding of the role of differential immune responses in driving coronavirus disease-2019 (COVID-19) disease severity and clinical outcomes has been greatly aided by SCG. With many variants of concern (VOC) in sight, it would be of great importance to further understand the immune response specificity vis-a-vis the immune cell repertoire, the identification of novel cell types, and antibody response. Given the potential of SCG to play an integral part in the multi-omics approach to the study of the host-pathogen interaction and its outcomes, our review attempts to highlight its strengths, its implications for infectious disease biology, and its current limitations. We conclude that the application of SCG would be a critical step towards future pandemic preparedness.

X-Chromosome Inactivation and Autosomal Random Monoallelic Expression as "Faux Amis"

X-chromosome inactivation (XCI) and random monoallelic expression of autosomal genes (RMAE) are two paradigms of gene expression regulation where, at the single cell level, genes can be expressed from either the maternal or paternal alleles. X-chromosome inactivation takes place in female marsupial and placental mammals, while RMAE has been described in mammals and also other species. Although the outcome of both processes results in random monoallelic expression and mosaicism at the cellular level, there are many important differences. We provide here a brief sketch of the history behind the discovery of XCI and RMAE. Moreover, we review some of the distinctive features of these two phenomena, with respect to when in development they are established, their roles in dosage compensation and cellular phenotypic diversity, and the molecular mechanisms underlying their initiation and stability.

Single-cell landscape of nuclear configuration and gene expression during stem cell differentiation and X inactivation

BACKGROUND

Mammalian development is associated with extensive changes in gene expression, chromatin accessibility, and nuclear structure. Here, we follow such changes associated with mouse embryonic stem cell differentiation and X inactivation by integrating, for the first time, allele-specific data from these three modalities obtained by high-throughput single-cell RNA-seq, ATAC-seq, and Hi-C.

RESULTS

Allele-specific contact decay profiles obtained by single-cell Hi-C clearly show that the inactive X chromosome has a unique profile in differentiated cells that have undergone X inactivation. Loss of this inactive X-specific structure at mitosis is followed by its reappearance during the cell cycle, suggesting a "bookmark" mechanism. Differentiation of embryonic stem cells to follow the onset of X inactivation is associated with changes in contact decay profiles that occur in parallel on both the X chromosomes and autosomes. Single-cell RNA-seq and ATAC-seq show evidence of a delay in female versus male cells, due to the presence of two active X chromosomes at early stages of differentiation. The onset of the inactive X-specific structure in single cells occurs later than gene silencing, consistent with the idea that chromatin compaction is a late event of X inactivation. Single-cell Hi-C highlights evidence of discrete changes in nuclear structure characterized by the acquisition of very long-range contacts throughout the nucleus. Novel computational approaches allow for the effective alignment of single-cell gene expression, chromatin accessibility, and 3D chromosome structure.

CONCLUSIONS

Based on trajectory analyses, three distinct nuclear structure states are detected reflecting discrete and profound simultaneous changes not only to the structure of the X chromosomes, but also to that of autosomes during differentiation. Our study reveals that long-range structural changes to chromosomes appear as discrete events, unlike progressive changes in gene expression and chromatin accessibility.

Patient derived stem cells for discovery and validation of novel pathogenic variants in inherited retinal disease

Our understanding of inherited retinal disease has benefited immensely from molecular genetic analysis over the past several decades. New technologies that allow for increasingly detailed examination of a patient's DNA have expanded the catalog of genes and specific variants that cause retinal disease. In turn, the identification of pathogenic variants has allowed the development of gene therapies and low-cost, clinically focused genetic testing. Despite this progress, a relatively large fraction (at least 20%) of patients with clinical features suggestive of an inherited retinal disease still do not have a molecular diagnosis today. Variants that are not obviously disruptive to the codon sequence of exons can be difficult to distinguish from the background of benign human genetic variations. Some of these variants exert their pathogenic effect not by altering the primary amino acid sequence, but by modulating gene expression, isoform splicing, or other transcript-level mechanisms. While not discoverable by DNA sequencing methods alone, these variants are excellent targets for studies of the retinal transcriptome. In this review, we present an overview of the current state of pathogenic variant discovery in retinal disease and identify some of the remaining barriers. We also explore the utility of new technologies, specifically patient-derived induced pluripotent stem cell (iPSC)-based modeling, in further expanding the catalog of disease-causing variants using transcriptome-focused methods. Finally, we outline bioinformatic analysis techniques that will allow this new method of variant discovery in retinal disease. As the knowledge gleaned from previous technologies is informing targets for therapies today, we believe that integrating new technologies, such as iPSC-based modeling, into the molecular diagnosis pipeline will enable a new wave of variant discovery and expanded treatment of inherited retinal disease.

Integrated analysis of Xist upregulation and X-chromosome inactivation with single-cell and single-allele resolution

To ensure dosage compensation between the sexes, one randomly chosen X chromosome is silenced in each female cell in the process of X-chromosome inactivation (XCI). XCI is initiated during early development through upregulation of the long non-coding RNA Xist, which mediates chromosome-wide gene silencing. Cell differentiation, Xist upregulation and gene silencing are thought to be coupled at multiple levels to ensure inactivation of exactly one out of two X chromosomes. Here we perform an integrated analysis of all three processes through allele-specific single-cell RNA-sequencing. Specifically, we assess the onset of random XCI in differentiating mouse embryonic stem cells, and develop dedicated analysis approaches. By exploiting the inter-cellular heterogeneity of XCI onset, we identify putative Xist regulators. Moreover, we show that transient Xist upregulation from both X chromosomes results in biallelic gene silencing right before transitioning to the monoallelic state, confirming a prediction of the stochastic model of XCI. Finally, we show that genetic variation modulates the XCI process at multiple levels, providing a potential explanation for the long-known X-controlling element (Xce) effect, which leads to preferential inactivation of a specific X chromosome in inter-strain crosses. We thus draw a detailed picture of the different levels of regulation that govern the initiation of XCI. The experimental and computational strategies we have developed here will allow us to profile random XCI in more physiological contexts, including primary human cells in vivo.

Contiguous erosion of the inactive X in human pluripotency concludes with global DNA hypomethylation

Female human pluripotent stem cells (hPSCs) routinely undergo inactive X (Xi) erosion. This progressive loss of key repressive features follows the loss of XIST expression, the long non-coding RNA driving X inactivation, and causes reactivation of silenced genes across the eroding X (Xe). To date, the sporadic and progressive nature of erosion has obscured its scale, dynamics, and key transition events. To address this problem, we perform an integrated analysis of DNA methylation (DNAme), chromatin accessibility, and gene expression across hundreds of hPSC samples. Differential DNAme orders female hPSCs across a trajectory from initiation to terminal Xi erosion. Our results identify a cis-regulatory element crucial for XIST expression, trace contiguously growing reactivated domains to a few euchromatic origins, and indicate that the late-stage Xe impairs DNAme genome-wide. Surprisingly, from this altered regulatory landscape emerge select features of naive pluripotency, suggesting that its link to X dosage may be partially conserved in human embryonic development.

Reactivation of the silenced X in human female iPSC/ESCs compromises their utility. Bansal et al. perform an integrated genomics analysis to reveal a prevalent X erosion trajectory that they validate in long-term culture. Starting with XIST loss, this trajectory indicates that reactivation may spread contiguously from escapees to silenced genes.

Sexual Dimorphism of the Heart: Genetics, Epigenetics, and Development

The democratization of genomic technologies has revealed profound sex biases in expression patterns in every adult tissue, even in organs with no conspicuous differences, such as the heart. With the increasing awareness of the disparities in cardiac pathophysiology between males and females, there are exciting opportunities to explore how sex differences in the heart are established developmentally. Although sexual dimorphism is traditionally attributed to hormonal influence, expression and epigenetic sex biases observed in early cardiac development can only be accounted for by the difference in sex chromosome composition, i.e., XX in females and XY in males. In fact, genes linked to the X and Y chromosomes, many of which encode regulatory factors, are expressed in cardiac progenitor cells and at every subsequent developmental stage. The effect of the sex chromosome composition may explain why many congenital heart defects originating before gonad formation exhibit sex biases in presentation, mortality, and morbidity. Some transcriptional and epigenetic sex biases established soon after fertilization persist in cardiac lineages, suggesting that early epigenetic events are perpetuated beyond early embryogenesis. Importantly, when sex hormones begin to circulate, they encounter a cardiac genome that is already functionally distinct between the sexes. Although there is a wealth of knowledge on the effects of sex hormones on cardiac function, we propose that sex chromosome-linked genes and their downstream targets also contribute to the differences between male and female hearts. Moreover, identifying how hormones influence sex chromosome effects, whether antagonistically or synergistically, will enhance our understanding of how sex disparities are established. We also explore the possibility that sexual dimorphism of the developing heart predicts sex-specific responses to environmental signals and foreshadows sex-biased health-related outcomes after birth.

satuRn: Scalable analysis of differential transcript usage for bulk and single-cell RNA-sequencing applications

Alternative splicing produces multiple functional transcripts from a single gene. Dysregulation of splicing is known to be associated with disease and as a hallmark of cancer. Existing tools for differential transcript usage (DTU) analysis either lack in performance, cannot account for complex experimental designs or do not scale to massive single-cell transcriptome sequencing (scRNA-seq) datasets. We introduce satuRn, a fast and flexible quasi-binomial generalized linear modelling framework that is on par with the best performing DTU methods from the bulk RNA-seq realm, while providing good false discovery rate control, addressing complex experimental designs, and scaling to scRNA-seq applications.

All models are wrong, but some are useful: Establishing standards for stem cell-based embryo models

Detailed studies of the embryo allow an increasingly mechanistic understanding of development, which has proved of profound relevance to human disease. The last decade has seen in vitro cultured stem cell-based models of embryo development flourish, which provide an alternative to the embryo for accessible experimentation. However, the usefulness of any stem cell-based embryo model will be determined by how accurately it reflects in vivo embryonic development, and/or the extent to which it facilitates new discoveries. Stringent benchmarking of embryo models is thus an important consideration for this growing field. Here we provide an overview of means to evaluate both the properties of stem cells, the building blocks of most embryo models, as well as the usefulness of current and future in vitro embryo models.

satuRn: Scalable analysis of differential transcript usage for bulk and single-cell RNA-sequencing applications

Alternative splicing produces multiple functional transcripts from a single gene. Dysregulation of splicing is known to be associated with disease and as a hallmark of cancer. Existing tools for differential transcript usage (DTU) analysis either lack in performance, cannot account for complex experimental designs or do not scale to massive scRNA-seq data. We introduce satuRn, a fast and flexible quasi-binomial generalized linear modelling framework that is on par with the best performing DTU methods from the bulk RNA-seq realm, while providing good false discovery rate control, addressing complex experimental designs and scaling to scRNA-seq applications.

Formative pluripotent stem cells show features of epiblast cells poised for gastrulation

The pluripotency of mammalian early and late epiblast could be recapitulated by naïve embryonic stem cells (ESCs) and primed epiblast stem cells (EpiSCs), respectively. However, these two states of pluripotency may not be sufficient to reflect the full complexity and developmental potency of the epiblast during mammalian early development. Here we report the establishment of self-renewing formative pluripotent stem cells (fPSCs) which manifest features of epiblast cells poised for gastrulation. fPSCs can be established from different mouse ESCs, pre-/early-gastrula epiblasts and induced PSCs. Similar to pre-/early-gastrula epiblasts, fPSCs show the transcriptomic features of formative pluripotency, which are distinct from naïve ESCs and primed EpiSCs. fPSCs show the unique epigenetic states of E6.5 epiblast, including the super-bivalency of a large set of developmental genes. Just like epiblast cells immediately before gastrulation, fPSCs can efficiently differentiate into three germ layers and primordial germ cells (PGCs) in vitro. Thus, fPSCs highlight the feasibility of using PSCs to explore the development of mammalian epiblast.

Advances in bulk and single-cell multi-omics approaches for systems biology and precision medicine

Multi-omics allows the systematic understanding of the information flow across different omics layers, while single omics can mainly reflect one aspect of the biological system. The advancement of bulk and single-cell sequencing technologies and related computational methods for multi-omics largely facilitated the development of system biology and precision medicine. Single-cell approaches have the advantage of dissecting cellular dynamics and heterogeneity, whereas traditional bulk technologies are limited to individual/population-level investigation. In this review, we first summarize the technologies for producing bulk and single-cell multi-omics data. Then, we survey the computational approaches for integrative analysis of bulk and single-cell multimodal data, respectively. Moreover, the databases and data storage for multi-omics, as well as the tools for visualizing multimodal data are summarized. We also outline the integration between bulk and single-cell data, and discuss the applications of multi-omics in precision medicine. Finally, we present the challenges and perspectives for multi-omics development.

Evaluating totipotency using criteria of increasing stringency

Totipotency is the ability of a single cell to give rise to all of the differentiated cell types that build the conceptus, yet how to capture this property in vitro remains incompletely understood. Defining totipotency relies on a variety of assays of variable stringency. Here, we describe criteria to define totipotency. We explain how distinct criteria of increasing stringency can be used to judge totipotency by evaluating candidate totipotent cell types in mice, including early blastomeres and expanded or extended pluripotent stem cells. Our data challenge the notion that expanded or extended pluripotent states harbour increased totipotent potential relative to conventional embryonic stem cells under in vitro and in vivo conditions.

Evaluating totipotency using criteria of increasing stringency

Totipotency is the ability of a single cell to give rise to all of the differentiated cell types that build the conceptus, yet how to capture this property in vitro remains incompletely understood. Defining totipotency relies on a variety of assays of variable stringency. Here, we describe criteria to define totipotency. We explain how distinct criteria of increasing stringency can be used to judge totipotency by evaluating candidate totipotent cell types in mice, including early blastomeres and expanded or extended pluripotent stem cells. Our data challenge the notion that expanded or extended pluripotent states harbour increased totipotent potential relative to conventional embryonic stem cells under in vitro and in vivo conditions.

X-chromosome regulation and sex differences in brain anatomy

Humans show reproducible sex-differences in cognition and psychopathology that may be contributed to by influences of gonadal sex-steroids and/or sex-chromosomes on regional brain development. Gonadal sex-steroids are well known to play a major role in sexual differentiation of the vertebrate brain, but far less is known regarding the role of sex-chromosomes. Our review focuses on this latter issue by bridging together two literatures that have to date been largely disconnected. We first consider "bottom-up" genetic and molecular studies focused on sex-chromosome gene content and regulation. This literature nominates specific sex-chromosome genes that could drive developmental sex-differences by virtue of their sex-biased expression and their functions within the brain. We then consider the complementary "top down" view, from magnetic resonance imaging studies that map sex- and sex chromosome effects on regional brain anatomy, and link these maps to regional gene-expression within the brain. By connecting these top-down and bottom-up approaches, we emphasize the potential role of X-linked genes in driving sex-biased brain development and outline key goals for future work in this field.

Of numbers and movement - understanding transcription factor pathogenesis by advanced microscopy

Transcription factors (TFs) are life-sustaining and, therefore, the subject of intensive research. By regulating gene expression, TFs control a plethora of developmental and physiological processes, and their abnormal function commonly leads to various developmental defects and diseases in humans. Normal TF function often depends on gene dosage, which can be altered by copy-number variation or loss-of-function mutations. This explains why TF haploinsufficiency (HI) can lead to disease. Since aberrant TF numbers frequently result in pathogenic abnormalities of gene expression, quantitative analyses of TFs are a priority in the field. In vitro single-molecule methodologies have significantly aided the identification of links between TF gene dosage and transcriptional outcomes. Additionally, advances in quantitative microscopy have contributed mechanistic insights into normal and aberrant TF function. However, to understand TF biology, TF-chromatin interactions must be characterised in vivo, in a tissue-specific manner and in the context of both normal and altered TF numbers. Here, we summarise the advanced microscopy methodologies most frequently used to link TF abundance to function and dissect the molecular mechanisms underlying TF HIs. Increased application of advanced single-molecule and super-resolution microscopy modalities will improve our understanding of how TF HIs drive disease.