Search

An increasing number of long noncoding RNAs (lncRNAs) have experimentally confirmed functions, yet little is known about their transcriptional dynamics and it is challenging to determine their regulatory effects. Here, we used allele-sensitive single-cell RNA sequencing to demonstrate that, compared to messenger RNAs, lncRNAs have twice as long duration between two transcriptional bursts. Additionally, we observed increased cell-to-cell variability in lncRNA expression due to lower frequency bursting producing larger numbers of RNA molecules. Exploiting heterogeneity in asynchronously growing cells, we identified and experimentally validated lncRNAs with cell state-specific functions involved in cell cycle progression and apoptosis. Finally, we identified cis-functioning lncRNAs and showed that knockdown of these lncRNAs modulated the nearby protein-coding gene’s transcriptional burst frequency or size. In summary, we identified distinct transcriptional regulation of lncRNAs and demonstrated a role for lncRNAs in the regulation of mRNA transcriptional bursting.

Excitatory projections from the lateral hypothalamic area (LHA) to the lateral habenula (LHb) drive aversive responses. We used Patch-seq guided multimodal classification to define the structural and functional heterogeneity of the LHA-LHb pathway. Our classification identified six glutamatergic neuron types with unique electrophysiological properties, molecular profiles, and projection patterns. We found that genetically-defined LHA-LHb neurons signal distinct aspects of emotional or naturalistic behaviors: Esr1+ LHA-LHb neurons induce aversion, whereas Npy+ LHA-LHb neurons control rearing behavior. Repeated optogenetic drive of Esr1+ LHA-LHb neurons induces a behaviorally persistent aversive state, and large-scale recordings showed a region-specific neural representation of the aversive state in the prelimbic region of the prefrontal cortex. We further found that exposure to unpredictable mild shocks induced a sex-specific sensitivity to develop a stress state in female mice, which was associated with a specific shift in the intrinsic properties of bursting-type Esr1+ LHA-LHb neurons. In summary, we describe the diversity of LHA-LHb neuron types, and provide evidence for the role of Esr1+ neurons in aversion and sexually dimorphic stress sensitivity.

Understanding the routing of neuronal information requires the functional characterization of connections. Neuronal projections recruit large postsynaptic ensembles with distinct postsynaptic response types (PRTs). PRT is probed by low-throughput whole-cell electrophysiology and is not a selection criterion for single-cell RNA-sequencing (scRNA-seq). To overcome these limitations and target neurons based on specific PRTs for soma harvesting and subsequent scRNA-seq we created Voltage-Seq. To test our methodology, we established all-optical voltage imaging and recorded the PRT of 8347 periaqueductal gray (PAG) neurons evoked by the optogenetic activation of ventromedial hypothalamic (VMH) terminals. PRTs were classified and spatially resolved in the entire VMH-PAG connectome. We built an on-site analysis named VoltView to navigate soma harvesting towards target PRTs guided by a classifier which used the VMH-PAG connectome database as a reference. We demonstrated the agility of Voltage-Seq in locating VMH-driven GABAergic neurons in the PAG, solely guided by the on-site classification in VoltView.

Throughout the SARS-CoV-2 pandemic, multiple waves of variants of concern have swept across populations, leading to a chain of new and yet more contagious lineages dominating COVID-19 cases. Here, we tracked the remarkably rapid shift from Omicron BA.1 to BA.2 sub-variant dominance in the Swedish population during January–March 2022. By analysis of 174,933 clinical nasopharyngeal swab samples using a custom variant-typing RT-PCR assay, we uncover nearly two-fold higher levels of viral RNA in cases with Omicron BA.2. Importantly, increased viral load in the upper pharynx upon BA.2 infection may provide part of the explanation why Omicron BA.2 is more transmissible and currently outcompetes the BA.1 variant across populations.

Transmission mechanisms for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are incompletely understood. In particular, aerosol transmission remains unclear, with viral detection in air and demonstration of its infection potential being actively investigated. To this end, we employed a novel electrostatic collector to sample air from rooms occupied by COVID-19 patients in a major Swedish hospital. Electrostatic air sampling in conjunction with extraction-free, reverse-transcriptase polymerase chain reaction (hid-RT-PCR) enabled detection of SARS-CoV-2 in air from patient rooms (9/22; 41%) and adjoining anterooms (10/22; 45%). Detection with hid-RT-PCR was concomitant with viral RNA presence on the surface of exhaust ventilation channels in patients and anterooms more than 2 m from the COVID-19 patient. Importantly, it was possible to detect active SARS-CoV-2 particles from room air, with a total of 496 plaque-forming units (PFUs) being isolated, establishing the presence of infectious, airborne SARS-CoV-2 in rooms occupied by COVID-19 patients. Our results support circulation of SARS-CoV-2 via aerosols and urge the revision of existing infection control frameworks to include airborne transmission. Correction in: Indoor Air, 2022, vol. 32, issue 8, e13085, doi: 10.1111/ina.13085

Joint single-cell measurements of gene expression and DNA regulatory element activity holds great promise as a tool to understand transcriptional regulation. Towards this goal we have developed Smart3-ATAC, a highly sensitive method which allows joint mRNA and chromatin accessibility analysis genome wide in single cells. With Smart3-ATAC, we are able to obtain the highest possible quality measurements per cell. The method combines transcriptomic profiling based on the highly sensitive Smart-seq3 protocol on cytosolic mRNA, with a novel low-loss single-cell ATAC (scATAC) protocol to measure chromatin accessibility. Compared to current droplet multiome methods, the yield of both the scATAC protocol and mRNA-seq protocol is markedly higher.

Throughout the SARS-CoV-2 pandemic, multiple waves of variants of concern have swept across populations, leading to a chain of new and yet more contagious variants dominating COVID-19 cases. Here, we tracked the remarkably rapid shift from Omicron BA.1 to BA.2 sublineage dominance in the Swedish population in early 2022 at a day-by-day basis.Using a custom SARS-CoV-2 Omicron BA.1 lineage-typing RT-PCR assay, we analyzed 174,933 clinical upper airway samples collected during January to March 2022.Our study demonstrates the feasibility and reliability of parallel lineage assignment of select variants at population scale, tracking the dominant sublineage transition from BA.1 to BA.2 at day-to-day resolution and uncovering nearly 2-fold higher levels of viral RNA in cases infected with Omicron BA.2 relative to BA.1.Our data provide unique insights into the Omicron BA.1 to BA.2 transition that occurred in Sweden during early 2022, and later, across the world. This may help to understand the increased transmissibility of the BA.2 variant.

X-chromosome inactivation (XCI) and upregulation (XCU) are the major opposing chromosome-wide modes of gene regulation that collectively achieve dosage compensation in mammals, but the regulatory link between the two remains elusive. Here, we use allele-resolved single-cell RNA-seq combined with chromatin accessibility profiling to finely dissect the separate effects of XCI and XCU on RNA levels during mouse development. We uncover that balanced X dosage is flexibly attained through expression tuning by XCU in a sex- and lineage-specific manner along varying degrees of XCI and across developmental and cellular states. Male blastomeres achieve XCU upon zygotic genome activation while females experience two distinct waves of XCU, upon imprinted- and random XCI, and ablation of Xist impedes female XCU. Contrary to widely established models of mammalian dosage compensation, naïve female embryonic cells carrying two active X chromosomes do not exhibit upregulation but express both alleles at basal level, yet collectively exceeding the RNA output of a single hyperactive allele. We show, in vivo and in vitro, that XCU is kinetically driven by X-specific modulation of transcriptional burst frequency, coinciding with increased compartmentalization of the hyperactive allele. Altogether, our data provide unprecedented insights into the dynamics of mammalian XCU, prompting a revised model of the chain in events of allelic regulation by XCU and XCI in unitedly achieving stable cellular levels of X-chromosome transcripts.

Summary: Following implantation, the epiblast (EPI) cells transit from the naive to primed pluripotency, accompanied by dynamic changes in X chromosome activity in females. To investigate the molecular attributes of this process, we performed single-cell RNA-seq analysis of 1,724 cells of E5.25, E5.5, E6.25, and E6.5 mouse embryos. We identified three cellular states in the EPI cells that capture the transition along the pluripotency continuum and the acquisition of primitive streak propensity. The transition of three EPI states was driven by inductive signaling activity emanating from the visceral endoderm (VE). In the EPI of female embryos, X chromosome reactivation (XCR) was initiated prior to the completion of imprinted X chromosome inactivation (XCI), and the ensuing random XCI was highly asynchronous. Moreover, imprinted paternal XCI proceeded faster in the VE than the extraembryonic ectoderm. Our study has provided a detailed molecular roadmap of the emergent lineage commitment before gastrulation and characterized X chromosome dynamics during early mouse development. : Cheng et al. present a molecular roadmap at single-cell and allelic resolution that highlights the developmental process of epiblast cells transiting through pluripotency states and acquiring the primitive streak propensity ahead of gastrulation. In the epiblast of female embryos, the paternal X chromosome is reactivated before the completion of imprinted inactivation. Keywords: pluripotency transition, epiblast heterogeneity, pregastrula development, X chromosome inactivation dynamics

RT-LAMP detection of SARS-CoV-2 has been shown to be a valuable approach to scale up COVID-19 diagnostics and thus contribute to limiting the spread of the disease. Here we present the optimization of highly cost-effective in-house produced enzymes, and we benchmark their performance against commercial alternatives. We explore the compatibility between multiple DNA polymerases with high strand-displacement activity and thermostable reverse transcriptases required for RT-LAMP. We optimize reaction conditions and demonstrate their applicability using both synthetic RNA and clinical patient samples. Finally, we validate the optimized RT-LAMP assay for the detection of SARS-CoV-2 in unextracted heat-inactivated nasopharyngeal samples from 184 patients. We anticipate that optimized and affordable reagents for RT-LAMP will facilitate the expansion of SARS-CoV-2 testing globally, especially in sites and settings where the need for large scale testing cannot be met by commercial alternatives.

International audience; The CD8+ T cell response to an antigen is composed of many T cell clones with unique T cell receptors, together forming a heterogeneous repertoire of effector and memory cells. How individual T cell clones contribute to this heterogeneity throughout immune responses remains largely unknown. In this study, we longitudinally track human CD8+ T cell clones expanding in response to yellow fever virus (YFV) vaccination at the single-cell level. We observed a drop in clonal diversity in blood from the acute to memory phase, suggesting that clonal selection shapes the circulating memory repertoire. Clones in the memory phase display biased differentiation trajectories along a gradient from stem cell to terminally differentiated effector memory fates. In secondary responses, YFV- and influenza-specific CD8+ T cell clones are poised to recapitulate skewed differentiation trajectories. Collectively, we show that the sum of distinct clonal phenotypes results in the multifaceted human T cell response to acute viral infections.

AbstractRT-LAMP detection of SARS-CoV-2 has been shown as a valuable approach to scale up COVID-19 diagnostics and thus contribute to limiting the spread of the disease. Here we present the optimization of highly cost-effective in-house produced enzymes, and we benchmark their performance against commercial alternatives. We explore the compatibility between multiple DNA polymerases with high strand-displacement activity and thermostable reverse transcriptases required for RT-LAMP. We optimize reaction conditions and demonstrate their applicability using both synthetic RNA and clinical patient samples. Finally, we validated the optimized RT-LAMP assay for the detection of SARS-CoV-2 in raw nasopharyngeal samples from 184 patients. We anticipate that optimized and affordable reagents for RT-LAMP will facilitate the expansion of SARS-CoV-2 testing globally, especially in sites and settings with limited economic resources.

The NK cell activation landscape in acute SARS-CoV-2 infection is associated with COVID-19 disease severity., Activated NK cells in severe COVID-19 Natural killer (NK) cells are cytotoxic lymphocytes that provide innate immune defense against viral infections and cancer, but little is known about their involvement in the host response to COVID-19. Maucourant et al. used high-dimensional flow cytometry to characterize NK cells in patients with moderate or severe COVID-19. SARS-CoV-2 infection was associated with fewer blood NK cells but a higher activation state in circulating NK cells. Severe COVID-19 resulted in an increase in “armed” NK cells containing high levels of cytotoxic proteins such as perforin. The adaptive NK subset was markedly expanded in a subset of severe patients. These findings lay the groundwork for future studies examining the mechanisms of NK cell activation in COVID-19 and their potential roles in host protection and immunopathology., Understanding innate immune responses in coronavirus disease 2019 (COVID-19) is important to decipher mechanisms of host responses and interpret disease pathogenesis. Natural killer (NK) cells are innate effector lymphocytes that respond to acute viral infections but might also contribute to immunopathology. Using 28-color flow cytometry, we here reveal strong NK cell activation across distinct subsets in peripheral blood of COVID-19 patients. This pattern was mirrored in single-cell RNA sequencing signatures of NK cells in bronchoalveolar lavage from COVID-19 patients. Unsupervised high-dimensional analysis of peripheral blood NK cells furthermore identified distinct NK cell immunotypes that were linked to disease severity. Hallmarks of these immunotypes were high expression of perforin, NKG2C, and Ksp37, reflecting increased presence of adaptive NK cells in circulation of patients with severe disease. Last, arming of CD56bright NK cells was observed across COVID-19 disease states, driven by a defined protein-protein interaction network of inflammatory soluble factors. This study provides a detailed map of the NK cell activation landscape in COVID-19 disease.

Understanding innate immune responses in COVID-19 is important for deciphering mechanisms of host responses and interpreting disease pathogenesis. Natural killer (NK) cells are innate effector lymphocytes that respond to acute viral infections, but might also contribute to immune pathology. Here, using 28-color flow cytometry, we describe a state of strong NK cell activation across distinct subsets in peripheral blood of COVID-19 patients, a pattern mirrored in scRNA-seq signatures of lung NK cells. Unsupervised high-dimensional analysis identified distinct immunophenotypes that were linked to disease severity. Hallmarks of these immunophenotypes were high expression of perforin, NKG2C, and Ksp37, reflecting a high presence of adaptive NK cell expansions in circulation of patients with severe disease. Finally, arming of CD56bright NK cells was observed in course of COVID-19 disease states, driven by a defined protein-protein interaction network of inflammatory soluble factors. This provides a detailed map of the NK cell activation-landscape in COVID-19 disease.

Mammalian X-chromosome dosage balance is regulated by X-chromosome inactivation (XCI) and X-chromosome upregulation (XCU), but the dynamics of XCU as well as the interplay between the two mechanisms remain poorly understood. Here, we mapped XCU throughout early mouse embryonic development at cellular and allelic resolution, revealing sex- and lineage-specific dynamics along key events in X-chromosome regulation. Our data show that XCU is linearly proportional to the degree of XCI, indicating that dosage compensation ensues based on mRNA levels rather than number of active X chromosomes. In line with this, we reveal that the two active X chromosomes in female naïve embryonic stem cells are not hyperactive as previously thought. In all lineages, XCU was underlain by increased transcriptional burst frequencies, providing a mechanistic basis in vivo. Together, our results demonstrate unappreciated flexibility of XCU in balancing X-chromosome expression, and we propose a general model for allelic dosage balance, applicable for wider mechanisms of transcriptional regulation.

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is commonly diagnosed by reverse transcription polymerase chain reaction (RT-PCR) to detect viral RNA in patient samples, but RNA extraction constitutes a major bottleneck in current testing. Methodological simplification could increase diagnostic availability and efficiency, benefitting patient care and infection control. Here, we describe methods circumventing RNA extraction in COVID-19 testing by performing RT-PCR directly on heat-inactivated or lysed samples. Our data, including benchmarking using 597 clinical patient samples and a standardised diagnostic system, demonstrate that direct RT-PCR is viable option to extraction-based tests. Using controlled amounts of active SARS-CoV-2, we confirm effectiveness of heat inactivation by plaque assay and evaluate various generic buffers as transport medium for direct RT-PCR. Significant savings in time and cost are achieved through RNA-extraction-free protocols that are directly compatible with established PCR-based testing pipelines. This could aid expansion of COVID-19 testing., SARS-CoV-2 infection is widely diagnosed by RT-PCR, but RNA extraction is a bottleneck for fast and cheap diagnosis. Here, the authors develop protocols to perform RT-PCR directly on heat-inactivated subject samples or samples lysed with readily available detergents and benchmark performance against 597 clinically diagnosed patient samples.

An increasing number of studies have demonstrated the regulatory importance of long non-coding RNAs (lncRNAs), yet little is known about their transcriptional dynamics and it remains challenging to determine their regulatory functions. Here, we used allele-sensitive single-cell RNA-seq (scRNA-seq) to demonstrate that lncRNAs have lower burst frequencies with twice as long duration between bursts, compared to mRNAs. Additionally, we observed an increased cell-to-cell variability in lncRNA expression that was due to more sporadic bursting (lower frequency) with larger numbers of RNA molecules being produced. Exploiting heterogeneity in asynchronously growing cells, we identified and experimentally validated lncRNAs with cell-state specific functions involved in cell cycle progression and apoptosis. Finally, utilizing allele-resolved RNA expression, we identified cis functioning lncRNAs and observed that knockdown of these lncRNAs modulated either transcriptional burst frequency or size of the nearby protein-coding gene. Collectively, our results identify distinct transcriptional regulation of lncRNAs and we demonstrate a role for lncRNAs in the regulation of transcriptional bursting of mRNAs.

The CD8

Immune cell differentiation is critical for adequate tissue-specific immune responses to occur. Here, we studied differentiation of human uterine natural killer cells (uNK cells). These cells reside in a tissue undergoing constant regeneration and represent the major leukocyte population at the maternal-fetal interface. However, their physiological response during the menstrual cycle and in pregnancy remains elusive. By surface proteome and transcriptome analysis as well as using humanized mice, we identify a differentiation pathway of uNK cells in vitro and in vivo with sequential acquisition of killer cell immunoglobulin-like receptors and CD39. uNK cell differentiation occurred continuously in response to the endometrial regeneration and was driven by interleukin-15. Differentiated uNK cells displayed reduced proliferative capacity and immunomodulatory function including enhanced angiogenic capacity. By studying human uterus transplantation and monozygotic twins, we found that the uNK cell niche could be replenished from circulation and that it was under genetic control. Together, our study uncovers a continuous differentiation pathway of human NK cells in the uterus that is coupled to profound functional changes in response to local tissue regeneration and pregnancy.

CD8+ T cells play essential roles in immunity to viral and bacterial infections, and to guard against malignant cells. The CD8+ T cell response to an antigen is composed of many T cell clones with unique T cell receptors, together forming a heterogenous repertoire of phenotypically and functionally distinct effector and memory cells. How individual T cell clones contribute to this heterogeneity throughout an immune response is key to understand immunity but remains largely unknown. Here, we longitudinally tracked hundreds of CD8+ T cell clones expanding in response to yellow fever virus vaccination at the single cell level in humans. We show that only a fraction of the clones detected in the acute phase of the response are detected as circulating memory T cells later, indicative of clonal selection shaping the memory repertoire. Clones persisting in the memory phase displayed biased differentiation trajectories along a gradient from stem cell memory to terminally differentiated effector memory fates. Reactivation of single memory CD8+ T cell clones revealed that they were poised to recapitulate skewed differentiation trajectories in secondary responses, and this was generalizable across individuals for both yellow fever and influenza virus. Together, we show that the sum of distinct clonal differentiation repertoires results in the multifaceted T cell response to acute viral infections in humans.

Our previous single-cell RNA-seq data from human preimplantation embryos showed that female X-chromosome mRNA levels become partly dose compensated during the timespan between zygotic genome activation (ZGA) and implantation. At the same time, XIST RNA is expressed from, and forms clouds in proximity to, both X-chromosome copies and biallelic expression of other X-linked genes persists. We proposed that X-chromosome transcription is transiently lowered on both alleles before X-chromosome inactivation (XCI) takes place. This notion was recently challenged in a reanalysis performed by Moreira de Mello et al, claiming to provide evidence against biallelic expression dampening and that instead proper XCI was responsible for the observed dosage compensation. Here we have addressed this reanalysis and highlighted methodological issues, and we conclude a current lack of evidence against biallelic X-chromosome dampening.

Transcriptional bursts render substantial biological noise in cellular transcriptomes. Here, we investigated the theoretical extent of monoallelic expression resulting from transcriptional bursting and how it compared to the amounts of monoallelic expression of autosomal genes observed in single-cell RNA-sequencing (scRNA-seq) data. We found that transcriptional bursting can explain frequent observations of autosomal monoallelic gene expression in cells. Importantly, the burst frequency largely determined the fraction of cells with monoallelic expression, whereas both burst frequency and size contributed to allelic imbalance. Allelic observations deviate from the expected when analysed across heterogeneous groups of cells, suggesting that allelic modelling can provide an unbiased assessment of heterogeneity within cells. Finally, large numbers of cells are required for analyses of allelic imbalance to avoid confounding observations from transcriptional bursting. Altogether, our results shed light on the implications of transcriptional burst kinetics on allelic expression patterns and phenotypic variation between cells.

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.

Ohno’s hypothesis postulates that X-chromosome upregulation rectifies X-dose imbalance relative to autosomal genes, present in two active copies per cell. Here we dissected X-upregulation into kinetics of transcription, inferred from allele-specific single-cell RNA-sequencing (scRNAseq) data from somatic mouse cells. We confirmed increased X-chromosome expression, and remarkably found that the X-chromosome achieved upregulation by elevated burst frequencies. This provides mechanistic insights into X-chromosome upregulation.

Mammalian gene expression is inherently stochastic1,2, and results in discrete bursts of RNA molecules that are synthesized from each allele3–7. Although transcription is known to be regulated by promoters and enhancers, it is unclear how cis-regulatory sequences encode transcriptional burst kinetics. Characterization of transcriptional bursting, including the burst size and frequency, has mainly relied on live-cell4,6,8 or single-molecule RNA fluorescence in situ hybridization3,5,8,9 recordings of selected loci. Here we determine transcriptome-wide burst frequencies and sizes for endogenous mouse and human genes using allele-sensitive single-cell RNA sequencing. We show that core promoter elements affect burst size and uncover synergistic effects between TATA and initiator elements, which were masked at mean expression levels. Notably, we provide transcriptome-wide evidence that enhancers control burst frequencies, and demonstrate that cell-type-specific gene expression is primarily shaped by changes in burst frequencies. Together, our data show that burst frequency is primarily encoded in enhancers and burst size in core promoters, and that allelic single-cell RNA sequencing is a powerful model for investigating transcriptional kinetics. Allele-specific single-cell RNA sequencing provides insights into transcription kinetics, with data indicating that core promoter sequences affect burst size, whereas enhancers mainly affect burst frequency.

Cellular heterogeneity can emerge from the expression of only one parental allele. However, it has remained controversial whether, or to what degree, random monoallelic expression of autosomal genes (aRME) is mitotically inherited (clonal) or stochastic (dynamic) in somatic cells, particularly in vivo. Here we used allele-sensitive single-cell RNA-seq on clonal primary mouse fibroblasts and freshly isolated human CD8+ T cells to dissect clonal and dynamic monoallelic expression patterns. Dynamic aRME affected a considerable portion of the cells' transcriptomes, with levels dependent on the cells' transcriptional activity. Notably, clonal aRME was detected, but it was surprisingly scarce (

Following implantation, the epiblast (EPI) cells transit from the naive to the primed pluripotency state in preparation for gastrulation, a momentous developmental milestone of the formation of the three primary germ layers. To investigate the molecular attributes of this critical developmental process of lineage specification and the dynamic changes in X-chromosome activity, we performed Smart-seq2 single-cell RNA-seq analysis of 1724 mouse E5.25 to E6.5 embryonic cells. We identified three distinct cellular states in the EPI cells that mirrors the transition of postimplantation EPI along the continuum of naive to primed pluripotency and the acquisition of primitive streak propensity over time. The transition of the EPI cellular state was accompanied by enhanced interaction with the visceral endoderm mediated by inductive signalling activity. By tracing X-chromosome activity at the allelic resolution, we discovered that X-chromosome reactivation (XCR) was initiated prior to the complete silencing of imprinted X-chromosome in the EPI, and followed by asynchronous random X-chromosome inactivation (XCI). Moreover, imprinted X-chromosome inactivation proceeded with a faster kinetics in the visceral endoderm compared to the extraembryonic ectoderm. Our study has outlined in detail a molecular roadmap of the emergent process of lineage specification and X-chromosome reprogramming during mouse pregastrulation development.

Summary: Determination of haplotypes is important for modelling the phenotypic consequences of genetic variation in diploid organisms, including cis-regulatory control and compound heterozygosity. We realized that single-cell RNA-seq (scRNA-seq) data are well suited for phasing genetic variants, since both transcriptional bursts and technical bottlenecks cause pronounced allelic fluctuations in individual single cells. Here we present scphaser, an R package that phases alleles at heterozygous variants to reconstruct haplotypes within transcribed regions of the genome using scRNA-seq data. The devised method efficiently and accurately reconstructed the known haplotype for ≥93% of phasable genes in both human and mouse. It also enables phasing of rare and de novo variants and variants far apart within genes, which is hard to attain with population-based computational inference. Availability and Implementation: scphaser is implemented as an R package. Tutorial and code are available at https://github.com/edsgard/scphaser Contact: rickard.sandberg@ki.se Supplementary information: Supplementary data are available at Bioinformatics online.

Massively parallel DNA sequencing of thousands of samples in a single machine-run is now possible, but the preparation of the individual sequencing libraries is expensive and time-consuming. Tagmentation-based library construction, using the Tn5 transposase, is efficient for generating sequencing libraries but currently relies on undisclosed reagents, which severely limits development of novel applications and the execution of large-scale projects. Here, we present simple and robust procedures for Tn5 transposase production and optimized reaction conditions for tagmentation-based sequencing library construction. We further show how molecular crowding agents both modulate library lengths and enable efficient tagmentation from subpicogram amounts of cDNA. The comparison of single-cell RNA-sequencing libraries generated using produced and commercial Tn5 demonstrated equal performances in terms of gene detection and library characteristics. Finally, because naked Tn5 can be annealed to any oligonucleotide of choice, for example, molecular barcodes in single-cell assays or methylated oligonucleotides for bisulfite sequencing, custom Tn5 production and tagmentation enable innovation in sequencing-based applications.

BackgroundSingle-cell RNA sequencing (scRNA-seq) offers exciting possibilities to address biological and medical questions, but a systematic comparison of recently developed protocols is still lacking.ResultsWe generated data from 447 mouse embryonic stem cells using Drop-seq, SCRB-seq, Smart-seq (on Fluidigm C1) and Smart-seq2 and analyzed existing data from 35 mouse embryonic stem cells prepared with CEL-seq. We find that Smart-seq2 is the most sensitive method as it detects the most genes per cell and across cells. However, it shows more amplification noise than CEL-seq, Drop-seq and SCRB-seq as it cannot use unique molecular identifiers (UMIs). We use simulations to model how the observed combinations of sensitivity and amplification noise affects detection of differentially expressed genes and find that SCRB-seq reaches 80% power with the fewest number of cells. When considering cost-efficiency at different sequencing depths at 80% power, we find that Drop-seq is preferable when quantifying transcriptomes of a large numbers of cells with low sequencing depth, SCRB-seq is preferable when quantifying transcriptomes of fewer cells and Smart-seq2 is preferable when annotating and/or quantifying transcriptomes of fewer cells as long one can use in-house produced transposase.ConclusionsOur analyses allows an informed choice among five prominent scRNA-seq protocols and provides a solid framework for benchmarking future improvements in scRNA-seq methodologies.

Pluripotency, differentiation, and X Chromosome inactivation (XCI) are key aspects of embryonic development. However, the underlying relationship and mechanisms among these processes remain unclear. Here, we systematically dissected these features along developmental progression using mouse embryonic stem cells (mESCs) and single-cell RNA sequencing with allelic resolution. We found that mESCs grown in a ground state 2i condition displayed transcriptomic profiles diffused from preimplantation mouse embryonic cells, whereas EpiStem cells closely resembled the post-implantation epiblast. Sex-related gene expression varied greatly across distinct developmental states. We also identified novel markers that were highly enriched in each developmental state. Moreover, we revealed that several novel pathways, including PluriNetWork and Focal Adhesion, were responsible for the delayed progression of female EpiStem cells. Importantly, we “digitalized” XCI progression using allelic expression of active and inactive X Chromosomes and surprisingly found that XCI states exhibited profound variability in each developmental state, including the 2i condition. XCI progression was not tightly synchronized with loss of pluripotency and increase of differentiation at the single-cell level, although these processes were globally correlated. In addition, highly expressed genes, including core pluripotency factors, were in general biallelically expressed. Taken together, our study sheds light on the dynamics of XCI progression and the asynchronicity between pluripotency, differentiation, and XCI.

Long noncoding RNAs (lncRNAs) regulate gene expression by association with chromatin, but how they target chromatin remains poorly understood. We have used chromatin RNA immunoprecipitation-coupled high-throughput sequencing to identify 276 lncRNAs enriched in repressive chromatin from breast cancer cells. Using one of the chromatin-interacting lncRNAs, MEG3, we explore the mechanisms by which lncRNAs target chromatin. Here we show that MEG3 and EZH2 share common target genes, including the TGF-β pathway genes. Genome-wide mapping of MEG3 binding sites reveals that MEG3 modulates the activity of TGF-β genes by binding to distal regulatory elements. MEG3 binding sites have GA-rich sequences, which guide MEG3 to the chromatin through RNA–DNA triplex formation. We have found that RNA–DNA triplex structures are widespread and are present over the MEG3 binding sites associated with the TGF-β pathway genes. Our findings suggest that RNA–DNA triplex formation could be a general characteristic of target gene recognition by the chromatin-interacting lncRNAs., Long noncoding RNAs (lncRNAs) regulate gene expression by association with chromatin. Here, the authors show that lncRNA MEG3 regulates the TGF-β pathway by bridging the interactions between polycomb repressive complex 2 and the distal regulatory elements of the TGF-β pathway genes via formation of RNA–DNA triplexes.

Single-cell RNA sequencing (scRNA-seq) offers new possibilities to address biological and medical questions. However, systematic comparisons of the performance of diverse scRNA-seq protocols are lacking. We generated data from 583 mouse embryonic stem cells to evaluate six prominent scRNA-seq methods: CEL-seq2, Drop-seq, MARS-seq, SCRB-seq, Smart-seq, and Smart-seq2. While Smart-seq2 detected the most genes per cell and across cells, CEL-seq2, Drop-seq, MARS-seq, and SCRB-seq quantified mRNA levels with less amplification noise due to the use of unique molecular identifiers (UMIs). Power simulations at different sequencing depths showed that Drop-seq is more cost-efficient for transcriptome quantification of large numbers of cells, while MARS-seq, SCRB-seq, and Smart-seq2 are more efficient when analyzing fewer cells. Our quantitative comparison offers the basis for an informed choice among six prominent scRNA-seq methods, and it provides a framework for benchmarking further improvements of scRNA-seq protocols.

Expressing One Allele at a Time Although genetic traits are often dominant or recessive, the impact of the same heterozygous genotype can vary quite a bit between individuals. Deng et al. (p. 193 ), analyzed global gene expression in hundreds of individual mouse cells and found that a substantial fraction of the genes only expressed one of the alleles, chosen randomly, at any given point in time. Such stochasticity in transcription increases the heterogeneity among cells and likely contributes to the phenotypic variance among individuals of identical genotype.

Linkage, association and expression studies previously pointed to the human QKI, KH domain containing, RNA-binding (QKI) as a candidate gene for schizophrenia. Functional studies of the mouse orthologue Qk focused mainly on its role in oligodendrocyte development and myelination, while its function in astroglia remained unexplored. Here, we show that QKI is highly expressed in human primary astrocytes and that its splice forms encode proteins targeting different subcellular localizations. Uncovering the role of QKI in astrocytes is of interest in light of growing evidence implicating astrocyte dysfunction in the pathogenesis of several disorders of the central nervous system. We selectively silenced QKI splice variants in human primary astrocytes and used RNA sequencing to identify differential expression and splice variant composition at the genome-wide level. We found that an mRNA expression of Glial fibrillary acidic protein (GFAP), encoding a major component of astrocyte intermediate filaments, was down-regulated after QKI7 splice variant silencing. Moreover, we identified a potential QKI-binding site within the 3' untranslated region of human GFAP. This sequence was not conserved between mice and humans, raising the possibility that GFAP is a target for QKI in humans but not rodents. Haloperidol treatment of primary astrocytes resulted in coordinated increases in QKI7 and GFAP expression. Taken together, our results provide the first link between QKI and GFAP, two genes with alterations previously observed independently in schizophrenic patients. Our findings for QKI, together with its well-known role in myelination, suggest that QKI is a hub regulator of glia function in humans.

Background Empirical evaluations of sexually dimorphic expression of genes on the mammalian X-chromosome are needed to understand the evolutionary forces and the gene-regulatory mechanisms controlling this chromosome. We performed a large-scale sex-bias expression analysis of genes on the X-chromosome in six different somatic tissues from mouse. Results Our results show that the mouse X-chromosome is enriched with female-biased genes and depleted of male-biased genes. This suggests that feminisation as well as de-masculinisation of the X-chromosome has occurred in terms of gene expression in non-reproductive tissues. Several mechanisms may be responsible for the control of female-biased expression on chromosome X, and escape from X-inactivation is a main candidate. We confirmed escape in case of Tmem29 using RNA-FISH analysis. In addition, we identified novel female-biased non-coding transcripts located in the same female-biased cluster as the well-known coding X-inactivation escapee Kdm5c, likely transcribed from the transition-region between active and silenced domains. We also found that previously known escapees only partially explained the overrepresentation of female-biased X-genes, particularly for tissue-specific female-biased genes. Therefore, the gene set we have identified contains tissue-specific escapees and/or genes controlled by other sexually skewed regulatory mechanisms. Analysis of gene age showed that evolutionarily old X-genes (>100 myr, preceding the radiation of placental mammals) are more frequently female-biased than younger genes. Conclusion Altogether, our results have implications for understanding both gene regulation and gene evolution of mammalian X-chromosomes, and suggest that the final result in terms of the X-gene composition (masculinisation versus feminisation) is a compromise between different evolutionary forces acting on reproductive and somatic tissues.