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

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.

New Insights into X-Chromosome Reactivation during Reprogramming to Pluripotency

Dosage compensation between the sexes results in one X chromosome being inactivated during female mammalian development. Chromosome-wide transcriptional silencing from the inactive X chromosome (Xi) in mammalian cells is erased in a process termed X-chromosome reactivation (XCR), which has emerged as a paradigm for studying the reversal of chromatin silencing. XCR is linked with germline development and induction of naive pluripotency in the epiblast, and also takes place upon reprogramming somatic cells to induced pluripotency. XCR depends on silencing of the long non-coding RNA (lncRNA) X inactive specific transcript (Xist) and is linked with the erasure of chromatin silencing. Over the past years, the advent of transcriptomics and epigenomics has provided new insights into the transcriptional and chromatin dynamics with which XCR takes place. However, multiple questions remain unanswered about how chromatin and transcription related processes enable XCR. Here, we review recent work on establishing the transcriptional and chromatin kinetics of XCR, as well as discuss a model by which transcription factors mediate XCR not only via Xist repression, but also by direct targeting of X-linked genes.

Membrane Tension Gates ERK-Mediated Regulation of Pluripotent Cell Fate

Cell fate transitions are frequently accompanied by changes in cell shape and mechanics. However, how cellular mechanics affects the instructive signaling pathways controlling cell fate is poorly understood. To probe the interplay between shape, mechanics, and fate, we use mouse embryonic stem cells (ESCs), which change shape as they undergo early differentiation. We find that shape change is regulated by a β-catenin-mediated decrease in RhoA activity and subsequent decrease in the plasma membrane tension. Strikingly, preventing a decrease in membrane tension results in early differentiation defects in ESCs and gastruloids. Decreased membrane tension facilitates the endocytosis of FGF signaling components, which activate ERK signaling and direct the exit from the ESC state. Increasing Rab5a-facilitated endocytosis rescues defective early differentiation. Thus, we show that a mechanically triggered increase in endocytosis regulates early differentiation. Our findings are of fundamental importance for understanding how cell mechanics regulates biochemical signaling and therefore cell fate.

Coverage-dependent bias creates the appearance of binary splicing in single cells

Single-cell RNA sequencing provides powerful insight into the factors that determine each cell’s unique identity. Previous studies led to the surprising observation that alternative splicing among single cells is highly variable and follows a bimodal pattern: a given cell consistently produces either one or the other isoform for a particular splicing choice, with few cells producing both isoforms. Here, we show that this pattern arises almost entirely from technical limitations. We analyze alternative splicing in human and mouse single-cell RNA-seq datasets, and model them with a probabilistic simulator. Our simulations show that low gene expression and low capture efficiency distort the observed distribution of isoforms. This gives the appearance of binary splicing outcomes, even when the underlying reality is consistent with more than one isoform per cell. We show that accounting for the true amount of information recovered can produce biologically meaningful measurements of splicing in single cells.

Single-cell sequencing of miRNAs: A modified technology

The recent development of next-generation sequencing technologies has offered valuable insights into individual cells. This technology is centered on the characterization of single cells for epigenomics, genomics, and transcriptomics. Ever since the first report appeared in 2009, the single-cell RNA-sequencing saga started to explore deeper into the mechanics intrigued within a single cell. microRNA (miRNA) has been increasingly recognized as an essential molecule triggering an additional layer for gene regulation. Therefore, single-cell sequencing of miRNAs is crucial to explore the logical riddles surrounding the epigenomics, genomics, and transcriptomics of an individual cell. Scientists from the Vienna Biocenter Campus have lately performed single-cell sequencing of miRNAs in the fly, Drosophila, and nematode, Caenorhabditis elegans. In this review, we present the latest scientific explorations supported by all-inclusive data on this novel subject matter of next-generation sequencing.

Locus-Specific Regulation of Xist Expression Using the CRISPR-Cas9-Based System

DNA methylation inhibitor or loss and gain of function of DNA methylation key players were widely used to investigate the regulation of X inactive-specific transcript (Xist) expression by DNA methylation, which results in global change of DNA methylation. Here, we reported a novel method for regulation of Xist using the widely used clustered regularly interspaced short palindromic repeat (CRISPR)-Cas system. First, Xist expression was increased in 5-aza-2'-deoxycytidine-treated female goat fibroblast cells. Second, three single-guide RNAs (sgRNAs) that target the Xist differential methylation region (DMR) were inserted to deactivated Cas9 (dCas9) nuclease and the catalytic domain of the DNA methyltransferase Dnmt3a coexpression plasmid. Bisulfite PCR analysis and quantitative real-time PCR revealed that the methylation level of the DMR was significantly increased, while the expression of Xist was downregulated in all three sgRNAs, compared with the mock-transfected cells. Third, the methylation activity at the sites of 37 bp from the protospacer-adjacent motif sequence showed the strong change relative to the mock-transfected cells. Furthermore, genome-wide DNA methylation and expression of the DNA methylation key players were not statistically changed in all three sgRNAs. Therefore, we confirmed that Xist expression was regulated by DNA methylation, and directed DNA methylation of Xist DMR at locus-specific solution decreased Xist expression.

A Bayesian mixture model for the analysis of allelic expression in single cells

Allele-specific expression (ASE) at single-cell resolution is a critical tool for understanding the stochastic and dynamic features of gene expression. However, low read coverage and high biological variability present challenges for analyzing ASE. We demonstrate that discarding multi-mapping reads leads to higher variability in estimates of allelic proportions, an increased frequency of sampling zeros, and can lead to spurious findings of dynamic and monoallelic gene expression. Here, we report a method for ASE analysis from single-cell RNA-Seq data that accurately classifies allelic expression states and improves estimation of allelic proportions by pooling information across cells. We further demonstrate that combining information across cells using a hierarchical mixture model reduces sampling variability without sacrificing cell-to-cell heterogeneity. We applied our approach to re-evaluate the statistical independence of allelic bursting and track changes in the allele-specific expression patterns of cells sampled over a developmental time course.

Dissecting primate early post-implantation development using long-term in vitro embryo culture

The transition from peri-implantation to gastrulation in mammals entails the specification and organization of the lineage progenitors into a body plan. Technical and ethical challenges have limited understanding of the cellular and molecular mechanisms that underlie this transition. We established a culture system that enabled the development of cynomolgus monkey embryos in vitro for up to 20 days. Cultured embryos underwent key primate developmental stages, including lineage segregation, bilaminar disc formation, amniotic and yolk sac cavitation, and primordial germ cell–like cell (PGCLC) differentiation. Single-cell RNA-sequencing analysis revealed development trajectories of primitive endoderm, trophectoderm, epiblast lineages, and PGCLCs. Analysis of single-cell chromatin accessibility identified transcription factors specifying each cell type. Our results reveal critical developmental events and complex molecular mechanisms underlying nonhuman primate embryogenesis in the early postimplantation period, with possible relevance to human development.

X-chromosome upregulation is driven by increased burst frequency

Ohno's hypothesis postulates that upregulation of X-linked genes rectifies their dosage imbalance relative to autosomal genes, which are present in two active copies per cell. Here we have dissected X-chromosome upregulation into the kinetics of transcription, inferred from allele-specific single-cell RNA sequencing data from somatic and embryonic mouse cells. We confirmed increased X-chromosome expression levels in female and male cells and found that the X chromosome achieved upregulation by elevated burst frequencies. By monitoring transcriptional kinetics in differentiating female mouse embryonic stem cells, we found that increased burst frequency was established on the active X chromosome when X inactivation took place on the other allele. Thus, our study provides mechanistic insights into X-chromosome upregulation.

X Inactivation and Escape: Epigenetic and Structural Features

X inactivation represents a complex multi-layer epigenetic mechanism that profoundly modifies chromatin composition and structure of one X chromosome in females. The heterochromatic inactive X chromosome adopts a unique 3D bipartite structure and a location close to the nuclear periphery or the nucleolus. X-linked lncRNA loci and their transcripts play important roles in the recruitment of proteins that catalyze chromatin and DNA modifications for silencing, as well as in the control of chromatin condensation and location of the inactive X chromosome. A subset of genes escapes X inactivation, raising questions about mechanisms that preserve their expression despite being embedded within heterochromatin. Escape gene expression differs between males and females, which can lead to physiological sex differences. We review recent studies that emphasize challenges in understanding the role of lncRNAs in the control of epigenetic modifications, structural features and nuclear positioning of the inactive X chromosome. Second, we highlight new findings about the distribution of genes that escape X inactivation based on single cell studies, and discuss the roles of escape genes in eliciting sex differences in health and disease.