Large-scale benchmarking of circRNA detection tools reveals large differences in sensitivity but not in precision.

The detection of circular RNA molecules (circRNAs) is typically based on short-read RNA sequencing data processed using computational tools. Numerous such tools have been developed, but a systematic comparison with orthogonal validation is missing. Here, we set up a circRNA detection tool benchmarking study, in which 16 tools detected more than 315,000 unique circRNAs in three deeply sequenced human cell types. Next, 1,516 predicted circRNAs were validated using three orthogonal methods. Generally, tool-specific precision is high and similar (median of 98.8%, 96.3% and 95.5% for qPCR, RNase R and amplicon sequencing, respectively) whereas the sensitivity and number of predicted circRNAs (ranging from 1,372 to 58,032) are the most significant differentiators. Of note, precision values are lower when evaluating low-abundance circRNAs. We also show that the tools can be used complementarily to increase detection sensitivity. Finally, we offer recommendations for future circRNA detection and validation.

Early fate bias in neuroepithelial progenitors of hippocampal neurogenesis.

Hippocampal adult neural stem cells emerge from progeny of the neuroepithelial lineage during murine brain development. Hippocampus development is increasingly well understood. However, the clonal relationships between early neuroepithelial stem cells and postnatal neurogenic cells remain unclear, especially at the single-cell level. Here we report fate bias and gene expression programs in thousands of clonally related cells in the juvenile hippocampus based on single-cell RNA-seq of barcoded clones. We find evidence for early fate restriction of neuroepithelial stem cells to either neurogenic progenitor cells of the dentate gyrus region or oligodendrogenic, non-neurogenic fate supplying cells for other hippocampal regions including gray matter areas and the Cornu ammonis region 1/3. Our study provides new insights into the phenomenon of early fate restriction guiding the development of postnatal hippocampal neurogenesis.

Loss of Y in leukocytes as a risk factor for critical COVID-19 in men.

BACKGROUND: The COVID-19 pandemic, which has a prominent social and economic impact worldwide, shows a largely unexplained male bias for the severity and mortality of the disease. Loss of chromosome Y (LOY) is a risk factor candidate in COVID-19 due to its prior association with many chronic age-related diseases, and its impact on immune gene transcription. METHODS: Publicly available scRNA-seq data of PBMC samples derived from male patients critically ill with COVID-19 were reanalyzed, and LOY status was added to the annotated cells. We further studied LOY in whole blood for 211 COVID-19 patients treated at intensive care units (ICU) from the first and second waves of the pandemic. Of these, 139 patients were subject to cell sorting for LOY analysis in granulocytes, low-density neutrophils (LDNs), monocytes, and PBMCs. RESULTS: Reanalysis of available scRNA-seq data revealed LDNs and monocytes as the cell types most affected by LOY. Subsequently, DNA analysis indicated that 46%, 32%, and 29% of critically ill patients showed LOY above 5% cut-off in LDNs, granulocytes, and monocytes, respectively. Hence, the myeloid lineage that is crucial for the development of severe COVID-19 phenotype is affected by LOY. Moreover, LOY correlated with increasing WHO score (median difference 1.59%, 95% HDI 0.46% to 2.71%, p=0.025), death during ICU treatment (median difference 1.46%, 95% HDI 0.47% to 2.43%, p=0.0036), and history of vessel disease (median difference 2.16%, 95% HDI 0.74% to 3.7%, p=0.004), among other variables. In 16 recovered patients, sampled during ICU stay and 93-143 days later, LOY decreased significantly in whole blood and PBMCs. Furthermore, the number of LDNs at the recovery stage decreased dramatically (median difference 76.4 per 10,000 cell sorting events, 95% HDI 55.5 to 104, p=6e-11). CONCLUSIONS: We present a link between LOY and an acute, life-threatening infectious disease. Furthermore, this study highlights LOY as the most prominent clonal mutation affecting the myeloid cell lineage during emergency myelopoiesis. The correlation between LOY level and COVID-19 severity might suggest that this mutation affects the functions of monocytes and neutrophils, which could have consequences for male innate immunity.

Clonal relations in the mouse brain revealed by single-cell and spatial transcriptomics.

The mammalian brain contains many specialized cells that develop from a thin sheet of neuroepithelial progenitor cells. Single-cell transcriptomics revealed hundreds of molecularly diverse cell types in the nervous system, but the lineage relationships between mature cell types and progenitor cells are not well understood. Here we show in vivo barcoding of early progenitors to simultaneously profile cell phenotypes and clonal relations in the mouse brain using single-cell and spatial transcriptomics. By reconstructing thousands of clones, we discovered fate-restricted progenitor cells in the mouse hippocampal neuroepithelium and show that microglia are derived from few primitive myeloid precursors that massively expand to generate widely dispersed progeny. We combined spatial transcriptomics with clonal barcoding and disentangled migration patterns of clonally related cells in densely labeled tissue sections. Our approach enables high-throughput dense reconstruction of cell phenotypes and clonal relations at the single-cell and tissue level in individual animals and provides an integrated approach for understanding tissue architecture.

Acute sleep loss alters circulating fibroblast growth factor 21 levels in humans: A randomised crossover trial.

The hormone fibroblast growth factor 21 (FGF21) modulates tissue metabolism and circulates at higher levels in metabolic conditions associated with chronic sleep-wake disruption, such as type 2 diabetes and obesity. In the present study, we investigated whether acute sleep loss impacts circulating levels of FGF21 and tissue-specific production, and response pathways linked to FGF21. A total of 15 healthy normal-weight young men participated in a randomised crossover study with two conditions, sleep loss versus an 8.5-hr sleep window. The evening before each intervention, fasting blood was collected. Fasting, post-intervention morning skeletal muscle and adipose tissue samples underwent quantitative polymerase chain reaction and DNA methylation analyses, and serum FGF21 levels were measured before and after an oral glucose tolerance test. Serum levels of FGF21 were higher after sleep loss compared with sleep, both under fasting conditions and following glucose intake (~27%-30%, p = 0.023). Fasting circulating levels of fibroblast activation protein, a protein which can degrade circulating FGF21, were not altered by sleep loss, whereas DNA methylation in the FGF21 promoter region increased only in adipose tissue. However, even though specifically the muscle exhibited transcriptional changes indicating adverse alterations to redox and metabolic homeostasis, no tissue-based changes were observed in expression of FGF21, its receptors, or selected signalling targets, in response to sleep loss. In summary, we found that acute sleep loss resulted in increased circulating levels of FGF21 in healthy young men, which may occur independent of a tissue-based stress response in metabolic peripheral tissues. Further studies may decipher whether changes in FGF21 signalling after sleep loss modulate metabolic outcomes associated with sleep or circadian disruption.

Profiling of Extracellular Small RNAs Highlights a Strong Bias towards Non-Vesicular Secretion.

The extracellular environment consists of a plethora of molecules, including extracellular miRNA that can be secreted in association with extracellular vesicles (EVs) or soluble protein complexes (non-EVs). Yet, interest in therapeutic short RNA carriers lies mainly in EVs, the vehicles conveying the great majority of the biological activity. Here, by overexpressing miRNA and shRNA sequences in parent cells and using size exclusion liquid chromatography (SEC) to separate the secretome into EV and non-EV fractions, we saw that >98% of overexpressed miRNA was secreted within the non-EV fraction. Furthermore, small RNA sequencing studies of native miRNA transcripts revealed that although the abundance of miRNAs in EVs, non-EVs and parent cells correlated well (R2 = 0.69-0.87), quantitatively an outstanding 96.2-99.9% of total miRNA was secreted in the non-EV fraction. Nevertheless, though EVs contained only a fraction of secreted miRNAs, these molecules were stable at 37 °C in a serum-containing environment, indicating that if sufficient miRNA loading is achieved, EVs can remain delivery-competent for a prolonged period of time. This study suggests that the passive endogenous EV loading strategy might be a relatively wasteful way of loading miRNA to EVs, and active miRNA loading approaches are needed for developing advanced EV miRNA therapies in the future.

Characterization of the nuclear and cytosolic transcriptomes in human brain tissue reveals new insights into the subcellular distribution of RNA transcripts.

Transcriptome analysis has mainly relied on analyzing RNA sequencing data from whole cells, overlooking the impact of subcellular RNA localization and its influence on our understanding of gene function, and interpretation of gene expression signatures in cells. Here, we separated cytosolic and nuclear RNA from human fetal and adult brain samples and performed a comprehensive analysis of cytosolic and nuclear transcriptomes. There are significant differences in RNA expression for protein-coding and lncRNA genes between cytosol and nucleus. We show that transcripts encoding the nuclear-encoded mitochondrial proteins are significantly enriched in the cytosol compared to the rest of protein-coding genes. Differential expression analysis between fetal and adult frontal cortex show that results obtained from the cytosolic RNA differ from results using nuclear RNA both at the level of transcript types and the number of differentially expressed genes. Our data provide a resource for the subcellular localization of thousands of RNA transcripts in the human brain and highlight differences in using the cytosolic or the nuclear transcriptomes for expression analysis.

Building blocks are synthesized on demand during the yeast cell cycle.

For cells to replicate, a sufficient supply of biosynthetic precursors is needed, necessitating the concerted action of metabolism and protein synthesis during progressive phases of cell division. A global understanding of which biosynthetic processes are involved and how they are temporally regulated during replication is, however, currently lacking. Here, quantitative multiomics analysis is used to generate a holistic view of the eukaryal cell cycle, using the budding yeast Saccharomyces cerevisiae Protein synthesis and central carbon pathways such as glycolysis and amino acid metabolism are shown to synchronize their respective abundance profiles with division, with pathway-specific changes in metabolite abundance also being reflected by a relative increase in mitochondrial volume, as shown by quantitative fluorescence microscopy. These results show biosynthetic precursor production to be temporally regulated to meet phase-specific demands of eukaryal cell division.

Expression profiling and in situ screening of circular RNAs in human tissues.

Circular RNAs (circRNAs) were recently discovered as a class of widely expressed noncoding RNA and have been implicated in regulation of gene expression. However, the function of the majority of circRNAs remains unknown. Studies of circRNAs have been hampered by a lack of essential approaches for detection, quantification and visualization. We therefore developed a target-enrichment sequencing method suitable for screening of circRNAs and their linear counterparts in large number of samples. We also applied padlock probes and in situ sequencing to visualize and determine circRNA localization in human brain tissue at subcellular levels. We measured circRNA abundance across different human samples and tissues. Our results highlight the potential of this RNA class to act as a specific diagnostic marker in blood and serum, by detection of circRNAs from genes exclusively expressed in the brain. The powerful and scalable tools we present will enable studies of circRNA function and facilitate screening of circRNA as diagnostic biomarkers.

Comprehensive RNA-Sequencing Analysis in Serum and Muscle Reveals Novel Small RNA Signatures with Biomarker Potential for DMD.

Extracellular small RNAs (sRNAs), including microRNAs (miRNAs), are promising biomarkers for diseases such as Duchenne muscular dystrophy (DMD), although their biological relevance is largely unknown. To investigate the relationship between intracellular and extracellular sRNA levels on a global scale, we performed sRNA sequencing in four muscle types and serum from wild-type, dystrophic mdx, and mdx mice in which dystrophin protein expression was restored by exon skipping. Differentially abundant sRNAs were identified in serum (mapping to miRNA, small nuclear RNA [snRNA], and PIWI-interacting RNA [piRNA] loci). One novel candidate biomarker, miR-483, was increased in both mdx serum and muscle, and also elevated in DMD patient sera. Dystrophin restoration induced global shifts in miRNA (including miR-483) and snRNA-fragment abundance toward wild-type levels. Specific serum piRNA-like sRNAs also responded to exon skipping therapy. Absolute miRNA expression in muscle was positively correlated with abundance in the circulation, although multiple highly expressed miRNAs in muscle were not elevated in mdx serum, suggesting that both passive and selective release mechanisms contribute to serum miRNA levels. In conclusion, this study has revealed new insights into the sRNA biology of dystrophin deficiency and identified novel DMD biomarkers.

Acute sleep loss results in tissue-specific alterations in genome-wide DNA methylation state and metabolic fuel utilization in humans.

Curtailed sleep promotes weight gain and loss of lean mass in humans, although the underlying molecular mechanisms are poorly understood. We investigated the genomic and physiological impact of acute sleep loss in peripheral tissues by obtaining adipose tissue and skeletal muscle after one night of sleep loss and after one full night of sleep. We find that acute sleep loss alters genome-wide DNA methylation in adipose tissue, and unbiased transcriptome-, protein-, and metabolite-level analyses also reveal highly tissue-specific changes that are partially reflected by altered metabolite levels in blood. We observe transcriptomic signatures of inflammation in both tissues following acute sleep loss, but changes involving the circadian clock are evident only in skeletal muscle, and we uncover molecular signatures suggestive of muscle breakdown that contrast with an anabolic adipose tissue signature. Our findings provide insight into how disruption of sleep and circadian rhythms may promote weight gain and sarcopenia.

Heterogeneity and interplay of the extracellular vesicle small RNA transcriptome and proteome.

Extracellular vesicles (EVs) mediate cell-to-cell communication by delivering or displaying macromolecules to their recipient cells. While certain broad-spectrum EV effects reflect their protein cargo composition, others have been attributed to individual EV-loaded molecules such as specific miRNAs. In this work, we have investigated the contents of vesicular cargo using small RNA sequencing of cells and EVs from HEK293T, RD4, C2C12, Neuro2a and C17.2. The majority of RNA content in EVs (49-96%) corresponded to rRNA-, coding- and tRNA fragments, corroborating with our proteomic analysis of HEK293T and C2C12 EVs which showed an enrichment of ribosome and translation-related proteins. On the other hand, the overall proportion of vesicular small RNA was relatively low and variable (2-39%) and mostly comprised of miRNAs and sequences mapping to piRNA loci. Importantly, this is one of the few studies, which systematically links vesicular RNA and protein cargo of vesicles. Our data is particularly useful for future work in unravelling the biological mechanisms underlying vesicular RNA and protein sorting and serves as an important guide in developing EVs as carriers for RNA therapeutics.

Genome-wide identification of Grainy head targets in Drosophila reveals regulatory interactions with the POU domain transcription factor Vvl.

Grainy head (Grh) is a conserved transcription factor (TF) controlling epithelial differentiation and regeneration. To elucidate Grh functions we identified embryonic Grh targets by ChIP-seq and gene expression analysis. We show that Grh controls hundreds of target genes. Repression or activation correlates with the distance of Grh-binding sites to the transcription start sites of its targets. Analysis of 54 Grh-responsive enhancers during development and upon wounding suggests cooperation with distinct TFs in different contexts. In the airways, Grh-repressed genes encode key TFs involved in branching and cell differentiation. Reduction of the POU domain TF Ventral veins lacking (Vvl) largely ameliorates the airway morphogenesis defects of grh mutants. Vvl and Grh proteins additionally interact with each other and regulate a set of common enhancers during epithelial morphogenesis. We conclude that Grh and Vvl participate in a regulatory network controlling epithelial maturation.

Massive and parallel expression profiling using microarrayed single-cell sequencing.

Single-cell transcriptome analysis overcomes problems inherently associated with averaging gene expression measurements in bulk analysis. However, single-cell analysis is currently challenging in terms of cost, throughput and robustness. Here, we present a method enabling massive microarray-based barcoding of expression patterns in single cells, termed MASC-seq. This technology enables both imaging and high-throughput single-cell analysis, characterizing thousands of single-cell transcriptomes per day at a low cost (0.13 USD/cell), which is two orders of magnitude less than commercially available systems. Our novel approach provides data in a rapid and simple way. Therefore, MASC-seq has the potential to accelerate the study of subtle clonal dynamics and help provide critical insights into disease development and other biological processes.

Visualization and analysis of gene expression in tissue sections by spatial transcriptomics.

Analysis of the pattern of proteins or messengerRNAs (mRNAs) in histological tissue sections is a cornerstone in biomedical research and diagnostics. This typically involves the visualization of a few proteins or expressed genes at a time. We have devised a strategy, which we call "spatial transcriptomics," that allows visualization and quantitative analysis of the transcriptome with spatial resolution in individual tissue sections. By positioning histological sections on arrayed reverse transcription primers with unique positional barcodes, we demonstrate high-quality RNA-sequencing data with maintained two-dimensional positional information from the mouse brain and human breast cancer. Spatial transcriptomics provides quantitative gene expression data and visualization of the distribution of mRNAs within tissue sections and enables novel types of bioinformatics analyses, valuable in research and diagnostics.

Common and distinct DNA-binding and regulatory activities of the BEN-solo transcription factor family.

Recently, the BEN (BANP, E5R, and NAC1) domain was recognized as a new class of conserved DNA-binding domain. The fly genome encodes three proteins that bear only a single BEN domain ("BEN-solo" factors); namely, Insensitive (Insv), Bsg25A (Elba1), and CG9883 (Elba2). Insv homodimers preferentially bind CCAATTGG palindromes throughout the genome to mediate transcriptional repression, whereas Bsg25A and Elba2 heterotrimerize with their obligate adaptor, Elba3 (i.e., the ELBA complex), to recognize a CCAATAAG motif in the Fab-7 insulator. While these data suggest distinct DNA-binding properties of BEN-solo proteins, we performed reporter assays that indicate that both Bsg25A and Elba2 can individually recognize Insv consensus sites efficiently. We confirmed this by solving the structure of Bsg25A complexed to the Insv site, which showed that key aspects of the BEN:DNA recognition strategy are similar between these proteins. We next show that both Insv and ELBA proteins are competent to mediate transcriptional repression via Insv consensus sequences but that the ELBA complex appears to be selective for the ELBA site. Reciprocally, genome-wide analysis reveals that Insv exhibits significant cobinding to class I insulator elements, indicating that it may also contribute to insulator function. Indeed, we observed abundant Insv binding within the Hox complexes with substantial overlaps with class I insulators, many of which bear Insv consensus sites. Moreover, Insv coimmunoprecipitates with the class I insulator factor CP190. Finally, we observed that Insv harbors exclusive activity among fly BEN-solo factors with respect to regulation of Notch-mediated cell fate choices in the peripheral nervous system. This in vivo activity is recapitulated by BEND6, a mammalian BEN-solo factor that conserves the Notch corepressor function of Insv but not its capacity to bind Insv consensus sites. Altogether, our data define an array of common and distinct biochemical and functional properties of this new family of transcription factors.

Genome-wide analysis of drosophila circular RNAs reveals their structural and sequence properties and age-dependent neural accumulation.

Circularization was recently recognized to broadly expand transcriptome complexity. Here, we exploit massive Drosophila total RNA-sequencing data, >5 billion paired-end reads from >100 libraries covering diverse developmental stages, tissues, and cultured cells, to rigorously annotate >2,500 fruit fly circular RNAs. These mostly derive from back-splicing of protein-coding genes and lack poly(A) tails, and the circularization of hundreds of genes is conserved across multiple Drosophila species. We elucidate structural and sequence properties of Drosophila circular RNAs, which exhibit commonalities and distinctions from mammalian circles. Notably, Drosophila circular RNAs harbor >1,000 well-conserved canonical miRNA seed matches, especially within coding regions, and coding conserved miRNA sites reside preferentially within circularized exons. Finally, we analyze the developmental and tissue specificity of circular RNAs and note their preferred derivation from neural genes and enhanced accumulation in neural tissues. Interestingly, circular isoforms increase substantially relative to linear isoforms during CNS aging and constitute an aging biomarker.

Diversity of miRNAs, siRNAs, and piRNAs across 25 Drosophila cell lines.

We expanded the knowledge base for Drosophila cell line transcriptomes by deeply sequencing their small RNAs. In total, we analyzed more than 1 billion raw reads from 53 libraries across 25 cell lines. We verify reproducibility of biological replicate data sets, determine common and distinct aspects of miRNA expression across cell lines, and infer the global impact of miRNAs on cell line transcriptomes. We next characterize their commonalities and differences in endo-siRNA populations. Interestingly, most cell lines exhibit enhanced TE-siRNA production relative to tissues, suggesting this as a common aspect of cell immortalization. We also broadly extend annotations of cis-NAT-siRNA loci, identifying ones with common expression across diverse cells and tissues, as well as cell-restricted loci. Finally, we characterize small RNAs in a set of ovary-derived cell lines, including somatic cells (OSS and OSC) and a mixed germline/somatic cell population (fGS/OSS) that exhibits ping-pong piRNA signatures. Collectively, the ovary data reveal new genic piRNA loci, including unusual configurations of piRNA-generating regions. Together with the companion analysis of mRNAs described in a previous study, these small RNA data provide comprehensive information on the transcriptional landscape of diverse Drosophila cell lines. These data should encourage broader usage of fly cell lines, beyond the few that are presently in common usage.

The histone demethylase activity of Rph1 is not essential for its role in the transcriptional response to nutrient signaling.

Rph1 and Gis1 are two related yeast zinc finger proteins that function as downstream effectors in the Ras/PKA, TOR and Sch9 nutrient signaling pathways. Both proteins also contain JmjC histone demethylase domains, but only Rph1 is known to be an active enzyme, demethylating lysine 36 of histone H3. We have studied to what extent the demethylase activity of Rph1 contributes to its role in nutrient signaling by performing gene expression microarray experiments on a yeast strain containing a catalytically inactive allele of RPH1. We find that the enzymatic activity of Rph1 is not essential for its role in growth phase dependent gene regulation. However, the ability of Rph1 to both activate and repress transcription is partially impaired in the active site mutant, indicating that the demethylase activity may enhance its function in vivo. Consistent with this, we find that the Rph1 mutation and a deletion of the histone H3 methylase Set2 affect the same target genes in opposite directions. Genes that are differentially expressed in the Rph1 mutant are also enriched for binding of Rpd3, a downstream effector in silencing, to their promoters. The expression of some subtelomeric genes and genes involved in sporulation and meiosis are also affected by the mutation, suggesting a role for Rph1-dependent demethylation in regulating these genes. A small set of genes are more strongly affected by the active site mutation, indicating a more pronounced role for the demethylase activity in their regulation by Rph1.

Intertwined pathways for Argonaute-mediated microRNA biogenesis in Drosophila.

Although Dicer is essential for general microRNA (miRNA) biogenesis, vertebrate mir-451 is Dicer independent. Instead, its short pre-miRNA hairpin is 'sliced' by Ago2, then 3'-resected into mature miRNAs. Here, we show that Drosophila cells and animals generate functional small RNAs from mir-451-type precursors. However, their bulk maturation arrests as Ago-cleaved pre-miRNAs, which mostly associate with the RNAi effector AGO2. Routing of pre-mir-451 hairpins to the miRNA effector AGO1 was inhibited by Dicer-1 and its partner Loqs. Loss of these miRNA factors promoted association of pre-mir-451 with AGO1, which sliced them and permitted maturation into ∼ 23-26 nt products. The difference was due to the 3' modification of single-stranded species in AGO2 by Hen1 methyltransferase, whose depletion permitted 3' trimming of Ago-cleaved pre-miRNAs in AGO2. Surprisingly, Nibbler, a 3'-5' exoribonuclease that trims 'long' mature miRNAs in AGO1, antagonized miR-451 processing. We used an in vitro reconstitution assay to identify a soluble, EDTA-sensitive activity that resects sliced pre-miRNAs in AGO1 complexes. Finally, we use deep sequencing to show that depletion of dicer-1 increases the diversity of small RNAs in AGO1, including some candidate mir-451-like loci. Altogether, we document unexpected aspects of miRNA biogenesis and Ago sorting, and provide insights into maturation of Argonaute-cleaved miRNA substrates.