Total RNA sequencing reveals nascent transcription and widespread co-transcriptional splicing in the human brain

Impact of commercial RNA extraction methods on the recovery of human RNA sequence data from archival fixed tissues

Archival fixed tissues hold key insights into the evolutionary history of RNA viruses and the associated host immune response, yet access to the RNA sequence data is limited by a lack of robust methods for RNA extraction and sequence retrieval from these tissue types. Here we compared three commercial RNA extraction techniques (bead, column, and phase-based) on five fixed human brain tissues done in triplicate, that have been stored for up to 43 years. We found that for this sample set, bead-based extractions captured longer molecules and yielded a greater proportion of unique reads when aligned to the human genome, than did column and phase-based extraction methods. Via the incorporation of multiple extraction replicates, we quantified the variability in sequencing metrics resulting from tissue sample and extraction technique heterogeneity. Additionally, we compared pre- and post-sequencing metrics and found that the former poorly predicted post-sequencing on-target success. Our findings help inform future research on the recovery of RNA from archival fixed tissues.

Laser-capture microdissection for spatial transcriptomics of immunohistochemically detected neurons

We developed a versatile 'IHC/LCM-Seq' method for spatial transcriptomics of immunohistochemically detected neurons collected with laser-capture microdissection (LCM). IHC/LCM-Seq uses aluminon and polyvinyl sulfonic acid for inventive RNA-preserving strategies to maintain RNA integrity in free-floating sections of 4% formaldehyde-fixed brains. To validate IHC/LCM-Seq, we first immunostained and harvested striatal cholinergic interneurons with LCM. RNA preparations were subjected to random primer-based cDNA library preparation and bulk sequencing on the NextSeq Illumina platform. IHC/LCM-Seq detected ∼16,000 transcripts, reaching the sensitivity of a reference 'LCM-Seq method' developed for fluorescently tagged neurons microdissected from lightly formaldehyde-fixed and slide-mounted brain sections of transgenic mice. We successfully used the new IHC/LCM-Seq approach to provide unprecedented insight into the transcriptome of immunohistochemically detected gonadotropin-releasing hormone (GnRH) neurons regulating reproduction. The ∼13,000 to 14,000 transcripts identified in GnRH neurons of adult male rats and mice encoded 28 proteins implicated previously in human infertility, 35 neuropeptides, 34 nuclear receptors, and 164 G protein-coupled receptors. Functional experiments using slice electrophysiology established that the heavy Ntsr2 expression conveys a strong excitatory action of neurotensin on GnRH neurons. As an unexpected species difference, we found that GnRH neurons exclusively expressed estrogen receptor-β in rats and against the current consensus, the alpha estrogen receptor isoform in mice. The IHC/LCM-Seq technique we are reporting is a highly sensitive and accurate bulk sequencing approach to characterize the transcriptome landscape of immunohistochemically labeled neurons, including neuroendocrine GnRH cells. This method is readily applicable to any species, opening new perspectives also for future studies of the post mortem human brain.

CircaKB: a comprehensive knowledgebase of circadian genes across multiple species

Circadian rhythms, which are the natural cycles that dictate various physiological processes over a 24-h period, have been increasingly recognized as important in the management and treatment of various human diseases. However, the lack of sufficient data and reliable analysis methods have been a major obstacle to understanding the bidirectional interaction between circadian variation and human health. We have developed CircaKB, a comprehensive knowledgebase of circadian genes across multiple species. CircaKB is the first knowledgebase that provides systematic annotations of the oscillatory patterns of gene expression at a genome-wide level for 15 representative species. Currently, CircaKB contains 226 time-course transcriptome datasets, covering a wide variety of tissues, organs, and cell lines. In addition, CircaKB integrates 12 computational models to facilitate reliable data analysis and identify oscillatory patterns and their variations in gene expression. CircaKB also offers powerful functionalities to its users, including easy search, fast browsing, strong visualization, and custom upload. We believe that CircaKB will be a valuable tool and resource for the circadian research community, contributing to the identification of new targets for disease prevention and treatment. We have made CircaKB freely accessible at https://cdsic.njau.edu.cn/CircaKB.

Decoding protein-RNA interactions using CLIP-based methodologies

Protein-RNA interactions are central to all RNA processing events, with pivotal roles in the regulation of gene expression and cellular functions. Dysregulation of these interactions has been increasingly linked to the pathogenesis of human diseases. High-throughput approaches to identify RNA-binding proteins and their binding sites on RNA - in particular, ultraviolet crosslinking followed by immunoprecipitation (CLIP) - have helped to map the RNA interactome, yielding transcriptome-wide protein-RNA atlases that have contributed to key mechanistic insights into gene expression and gene-regulatory networks. Here, we review these recent advances, explore the effects of cellular context on RNA binding, and discuss how these insights are shaping our understanding of cellular biology. We also review the potential therapeutic applications arising from new knowledge of protein-RNA interactions.

Post-transcriptional regulation supports the homeostatic expression of mature RNA

Gene expression regulation is a sophisticated, multi-stage process, and its robustness is critical to normal cell function and the survival of an organism. Previous studies indicate that differential gene expression at the RNA level is typically attenuated at the protein level through translational regulation. However, how post-transcriptional regulation (PTR) influences expression change during the RNA maturation process remains unclear. In this study, we investigated this by quantifying the magnitude of expression change in precursor RNA and mature RNA across a vast range of different biological conditions. We analyzed bulk tissue RNA sequencing data from 4689 samples, including healthy and diseased tissues from human, chimpanzee, rhesus macaque, and murine sources. We demonstrated that PTR tends to support homeostatic expression of mature RNA by amplifying normal tissue-specific expression of precursor RNA, while reducing expression change of precursor RNA in disease contexts. Our study provides insight into the general influence of PTR on gene expression homeostasis. Our analysis also suggests that intronic reads in RNA-seq studies may contain under-utilized information about disease associations. Additionally, our findings may assist in identifying new disease biomarkers and more effective ways of altering gene expression as a therapeutic strategy.

The 3D Genome in Brain Development: An Exploration of Molecular Mechanisms and Experimental Methods

The human brain contains multiple cell types that are spatially organized into functionally distinct regions. The proper development of the brain requires complex gene regulation mechanisms in both neurons and the non-neuronal cell types that support neuronal function. Studies across the last decade have discovered that the 3D nuclear organization of the genome is instrumental in the regulation of gene expression in the diverse cell types of the brain. In this review, we describe the fundamental biochemical mechanisms that regulate the 3D genome, and comprehensively describe in vitro and ex vivo studies on mouse and human brain development that have characterized the roles of the 3D genome in gene regulation. We highlight the significance of the 3D genome in linking distal enhancers to their target promoters, which provides insights on the etiology of psychiatric and neurological disorders, as the genetic variants associated with these disorders are primarily located in noncoding regulatory regions. We also describe the molecular mechanisms that regulate chromatin folding and gene expression in neurons. Furthermore, we describe studies with an evolutionary perspective, which have investigated features that are conserved from mice to human, as well as human gained 3D chromatin features. Although most of the insights on disease and molecular mechanisms have been obtained from bulk 3C based experiments, we also highlight other approaches that have been developed recently, such as single cell 3C approaches, as well as non-3C based approaches. In our future perspectives, we highlight the gaps in our current knowledge and emphasize the need for 3D genome engineering and live cell imaging approaches to elucidate mechanisms and temporal dynamics of chromatin interactions, respectively.

Dynamic interactions drive early spliceosome assembly

Splicing is a critical processing step during pre-mRNA maturation in eukaryotes. The correct selection of splice sites during the early steps of spliceosome assembly is highly important and crucial for the regulation of alternative splicing. Splice site recognition and alternative splicing depend on cis-regulatory sequence elements in the RNA and trans-acting splicing factors that recognize these elements and crosstalk with the canonical splicing machinery. Structural mechanisms involving early spliceosome complexes are governed by dynamic RNA structures, protein-RNA interactions and conformational flexibility of multidomain RNA binding proteins. Here, we highlight structural studies and integrative structural biology approaches, which provide complementary information from cryo-EM, NMR, small angle scattering, and X-ray crystallography to elucidate mechanisms in the regulation of early spliceosome assembly and quality control, highlighting the role of conformational dynamics.

NERD-seq: a novel approach of Nanopore direct RNA sequencing that expands representation of non-coding RNAs

Non-coding RNAs (ncRNAs) are frequently documented RNA modification substrates. Nanopore Technologies enables the direct sequencing of RNAs and the detection of modified nucleobases. Ordinarily, direct RNA sequencing uses polyadenylation selection, studying primarily mRNA gene expression. Here, we present NERD-seq, which enables detection of multiple non-coding RNAs, excluded by the standard approach, alongside natively polyadenylated transcripts. Using neural tissues as a proof of principle, we show that NERD-seq expands representation of frequently modified non-coding RNAs, such as snoRNAs, snRNAs, scRNAs, srpRNAs, tRNAs, and rRFs. NERD-seq represents an RNA-seq approach to simultaneously study mRNA and ncRNA epitranscriptomes in brain tissues and beyond.

The Functional Relationship Between RNA Splicing and the Chromatin Landscape

Chromatin is a highly regulated and dynamic structure that has been shown to play an essential role in transcriptional and co-transcriptional regulation. In the context of RNA splicing, early evidence suggested a loose connection between the chromatin landscape and splicing. More recently, it has been shown that splicing occurs in a co-transcriptional manner, meaning that the splicing process occurs in the context of chromatin. Experimental and computational evidence have also shown that chromatin dynamics can influence the splicing process and vice versa. However, much of this evidence provides mainly correlative relationships between chromatin and splicing with just a few concrete examples providing defined molecular mechanisms by which these two processes are functionally related. Nevertheless, it is clear that chromatin and RNA splicing are tightly interconnected to one another. In this review, we highlight the current state of knowledge of the relationship between chromatin and splicing.

A critical review of the impact of candidate copy number variants on autism spectrum disorder

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder (NDD) influenced by genetic, epigenetic, and environmental factors. Recent advancements in genomic analysis have shed light on numerous genes associated with ASD, highlighting the significant role of both common and rare genetic mutations, as well as copy number variations (CNVs), single nucleotide polymorphisms (SNPs) and unique de novo variants. These genetic variations disrupt neurodevelopmental pathways, contributing to the disorder's complexity. Notably, CNVs are present in 10 %-20 % of individuals with autism, with 3 %-7 % detectable through cytogenetic methods. While the role of submicroscopic CNVs in ASD has been recently studied, their association with genomic loci and genes has not been thoroughly explored. In this review, we focus on 47 CNV regions linked to ASD, encompassing 1632 genes, including protein-coding genes and long non-coding RNAs (lncRNAs), of which 659 show significant brain expression. Using a list of ASD-associated genes from SFARI, we detect 17 regions harboring at least one known ASD-related protein-coding gene. Of the remaining 30 regions, we identify 24 regions containing at least one protein-coding gene with brain-enriched expression and a nervous system phenotype in mouse mutants, and one lncRNA with both brain-enriched expression and upregulation in iPSC to neuron differentiation. This review not only expands our understanding of the genetic diversity associated with ASD but also underscores the potential of lncRNAs in contributing to its etiology. Additionally, the discovered CNVs will be a valuable resource for future diagnostic, therapeutic, and research endeavors aimed at prioritizing genetic variations in ASD.

Defining the relationship between cellular and extracellular vesicle (EV) content in breast cancer via an integrative multi-omic analysis

Much recent research has been dedicated to exploring the utility of extracellular vesicles (EVs) as circulating disease biomarkers. Underpinning this work is the assumption that the molecular cargo of EVs directly reflects the originating cell. Few attempts have been made, however, to empirically validate this on the -omic level. To this end, we have performed an integrative multi-omic analysis of a panel of breast cancer cell lines and corresponding EVs. Whole transcriptome analysis validated that the cellular transcriptome remained stable when cultured cells are transitioned to low serum or serum-free medium for EV collection. Transcriptomic profiling of the isolated EVs indicated a positive correlation between transcript levels in cells and EVs, including disease-associated transcripts. Analysis of the EV proteome verified that HER2 protein is present in EVs, however neither the estrogen (ER) nor progesterone (PR) receptor proteins are detected regardless of cellular expression. Using multivariate analysis, we derived an EV protein signature to infer cellular patterns of ER and HER2 expression, though the ER protein could not be directly detected. Integrative analyses affirmed that the EV proteome and transcriptome captured key phenotypic hallmarks of the originating cells, supporting the potential of EVs for non-invasive monitoring of breast cancers.

A single workflow for multi-species blood transcriptomics

BACKGROUND

Blood transcriptomic analysis is widely used to provide a detailed picture of a physiological state with potential outcomes for applications in diagnostics and monitoring of the immune response to vaccines. However, multi-species transcriptomic analysis is still a challenge from a technological point of view and a standardized workflow is urgently needed to allow interspecies comparisons.

RESULTS

Here, we propose a single and complete total RNA-Seq workflow to generate reliable transcriptomic data from blood samples from humans and from animals typically used in preclinical models. Blood samples from a maximum of six individuals and four different species (rabbit, non-human primate, mouse and human) were extracted and sequenced in triplicates. The workflow was evaluated using different wet-lab and dry-lab criteria, including RNA quality and quantity, the library molarity, the number of raw sequencing reads, the Phred-score quality, the GC content, the performance of ribosomal-RNA and globin depletion, the presence of residual DNA, the strandness, the percentage of coding genes, the number of genes expressed, and the presence of saturation plateau in rarefaction curves. We identified key criteria and their associated thresholds to be achieved for validating the transcriptomic workflow. In this study, we also generated an automated analysis of the transcriptomic data that streamlines the validation of the dataset generated.

CONCLUSIONS

Our study has developed an end-to-end workflow that should improve the standardization and the inter-species comparison in blood transcriptomics studies. In the context of vaccines and drug development, RNA sequencing data from preclinical models can be directly compared with clinical data and used to identify potential biomarkers of value to monitor safety and efficacy.

KDM3A and KDM3B Maintain Naïve Pluripotency Through the Regulation of Alternative Splicing

Histone modifying enzymes play a central role in maintaining cell identity by establishing a conducive chromatin environment for lineage specific transcription factor activity. Pluripotent embryonic stem cell (ESC) identity is characterized by a lower abundance of gene repression associated histone modifications that enables rapid response to differentiation cues. The KDM3 family of histone demethylases removes the repressive histone H3 lysine 9 dimethylation (H3K9me2). Here we uncover a surprising role for the KDM3 proteins in the maintenance of the pluripotent state through post-transcriptional regulation. We find that KDM3A and KDM3B interact with RNA processing factors such as EFTUD2 and PRMT5. Acute selective degradation of the endogenous KDM3A and KDM3B proteins resulted in altered splicing independent of H3K9me2 status or catalytic activity. These splicing changes partially resemble the splicing pattern of the more blastocyst-like ground state of pluripotency and occurred in important chromatin and transcription factors such as Dnmt3b, Tbx3 and Tcf12. Our findings reveal non-canonical roles of histone demethylating enzymes in splicing to regulate cell identity.

Hitting the brakes on transcription to extend lifespan

In aged animals from worms to humans, transcriptional elongation rates are faster, leading to changes in transcript quality and alternative splicing. Recent work by Debès et al. shows how interventions that slow elongation rates, such as mutating RNA polymerase II (Pol II) or increasing nucleosome density to impose transcriptional 'traffic', delay senescence and promote longevity.

Transcriptional elongation control in developmental gene expression, aging, and disease

The elongation stage of transcription by RNA polymerase II (RNA Pol II) is central to the regulation of gene expression in response to developmental and environmental cues in metazoan. Dysregulated transcriptional elongation has been associated with developmental defects as well as disease and aging processes. Decades of genetic and biochemical studies have painstakingly identified and characterized an ensemble of factors that regulate RNA Pol II elongation. This review summarizes recent findings taking advantage of genetic engineering techniques that probe functions of elongation factors in vivo. We propose a revised model of elongation control in this accelerating field by reconciling contradictory results from the earlier biochemical evidence and the recent in vivo studies. We discuss how elongation factors regulate promoter-proximal RNA Pol II pause release, transcriptional elongation rate and processivity, RNA Pol II stability and RNA processing, and how perturbation of these processes is associated with developmental disorders, neurodegenerative disease, cancer, and aging.

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

Diurnal oscillations of epigenetic modifications are associated with variation in rhythmic expression of homoeologous genes in Brassica napus

BACKGROUND

Epigenetic modifications that exhibit circadian oscillations also promote circadian oscillations of gene expression. Brassica napus is a heterozygous polyploid species that has undergone distant hybridization and genome doubling events and has a young and distinct species origin. Studies incorporating circadian rhythm analysis of epigenetic modifications can offer new insights into differences in diurnal oscillation behavior among subgenomes and the regulation of diverse expressions of homologous gene rhythms in biological clocks.

RESULTS

In this study, we created a high-resolution and multioscillatory gene expression dataset, active histone modification (H3K4me3, H3K9ac), and RNAPII recruitment in Brassica napus. We also conducted the pioneering characterization of the diurnal rhythm of transcription and epigenetic modifications in an allopolyploid species. We compared the evolution of diurnal rhythms between subgenomes and observed that the Cn subgenome had higher diurnal oscillation activity in both transcription and active histone modifications than the An subgenome. Compared to the A subgenome in Brassica rapa, the An subgenome of Brassica napus displayed significant changes in diurnal oscillation characteristics of transcription. Homologous gene pairs exhibited a higher proportion of diurnal oscillation in transcription than subgenome-specific genes, attributed to higher chromatin accessibility and abundance of active epigenetic modification types. We found that the diurnal expression of homologous genes displayed diversity, and the redundancy of the circadian system resulted in extensive changes in the diurnal rhythm characteristics of clock genes after distant hybridization and genome duplication events. Epigenetic modifications influenced the differences in the diurnal rhythm of homologous gene expression, and the diurnal oscillation of homologous gene expression was affected by the combination of multiple histone modifications.

CONCLUSIONS

Herein, we presented, for the first time, a characterization of the diurnal rhythm characteristics of gene expression and its epigenetic modifications in an allopolyploid species. Our discoveries shed light on the epigenetic factors responsible for the diurnal oscillation activity imbalance between subgenomes and homologous genes' rhythmic expression differences. The comprehensive time-series dataset we generated for gene expression and epigenetic modifications provides a valuable resource for future investigations into the regulatory mechanisms of protein-coding genes in Brassica napus.

Identification of Intron Retention in the Slc16a3 Gene Transcript Encoding the Transporter MCT4 in the Brain of Aged and Alzheimer-Disease Model (APPswePS1dE9) Mice

The monocarboxylate transporter 4 (MCT4; Slc16a3) is expressed in the central nervous system, notably by astrocytes. It is implicated in lactate release and the regulation of glycolytic flux. Whether its expression varies during normal and/or pathological aging is unclear. As the presence of its mature transcript in the brain of young and old mice was determined, an unexpectedly longer RT-PCR fragment was detected in the mouse frontal cortex and hippocampus at 12 vs. 3 months of age. Cultured astrocytes expressed the expected 516 base pair (bp) fragment but treatment with IL-1β to mimic inflammation as can occur during aging led to the additional expression of a 928 bp fragment like that seen in aged mice. In contrast, cultured pericytes (a component of the blood-brain barrier) only exhibited the 516 bp fragment. Intriguingly, cultured endothelial cells constitutively expressed both fragments. When RT-PCR was performed on brain subregions of an Alzheimer mouse model (APPswePS1dE9), no fragment was detected at 3 months, while only the 928 bp fragment was present at 12 months. Sequencing of MCT4 RT-PCR products revealed the presence of a remaining intron between exon 2 and 3, giving rise to the longer fragment detected by RT-PCR. These results unravel the existence of intron retention for the MCT4 gene in the central nervous system. Such alternative splicing appears to increase with age in the brain and might be prominent in neurodegenerative diseases such as Alzheimer's disease. Hence, further studies in vitro and in vivo of intron 2 retention in the Slc16a3 gene transcript are required for adequate characterization concerning the biological roles of Slc16a3 isoforms in the context of aging and Alzheimer's disease pathology.

Seizure enhances SUMOylation and zinc-finger transcriptional repression in neuronal nuclei

A single episode of pilocarpine-induced status epilepticus can trigger the development of spontaneous recurrent seizures in a rodent model for epilepsy. The initial seizure-induced events in neuronal nuclei that lead to long-term changes in gene expression and cellular responses likely contribute toward epileptogenesis. Using a transgenic mouse model to specifically isolate excitatory neuronal nuclei, we profiled the seizure-induced nuclear proteome via tandem mass tag mass spectrometry and observed robust enrichment of nuclear proteins associated with the SUMOylation pathway. In parallel with nuclear proteome, we characterized nuclear gene expression by RNA sequencing which provided insights into seizure-driven transcriptional regulation and dynamics. Strikingly, we saw widespread downregulation of zinc-finger transcription factors, specifically proteins that harbor Krüppel-associated box (KRAB) domains. Our results provide a detailed snapshot of nuclear events induced by seizure activity and demonstrate a robust method for cell-type-specific nuclear profiling that can be applied to other cell types and models.

Fine-grained cell-type specific association studies with human bulk brain data using a large single-nucleus RNA sequencing based reference panel

Brain disorders are leading causes of disability worldwide. Gene expression studies provide promising opportunities to better understand their etiology but it is critical that expression is studied on a cell-type level. Cell-type specific association studies can be performed with bulk expression data using statistical methods that capitalize on cell-type proportions estimated with the help of a reference panel. To create a fine-grained reference panel for the human prefrontal cortex, we performed an integrated analysis of the seven largest single nucleus RNA-seq studies. Our panel included 17 cell-types that were robustly detected across all studies, subregions of the prefrontal cortex, and sex and age groups. To estimate the cell-type proportions, we used an empirical Bayes estimator that substantially outperformed three estimators recommended previously after a comprehensive evaluation of methods to estimate cell-type proportions from brain transcriptome data. This is important as being able to precisely estimate the cell-type proportions may avoid unreliable results in downstream analyses particularly for the multiple cell-types that had low abundances. Transcriptome-wide association studies performed with permuted bulk expression data showed that it is possible to perform transcriptome-wide association studies for even the rarest cell-types without an increased risk of false positives.

Insplico: effective computational tool for studying splicing order of adjacent introns genome-wide with short and long RNA-seq reads

Although splicing occurs largely co-transcriptionally, the order by which introns are removed does not necessarily follow the order in which they are transcribed. Whereas several genomic features are known to influence whether or not an intron is spliced before its downstream neighbor, multiple questions related to adjacent introns' splicing order (AISO) remain unanswered. Here, we present Insplico, the first standalone software for quantifying AISO that works with both short and long read sequencing technologies. We first demonstrate its applicability and effectiveness using simulated reads and by recapitulating previously reported AISO patterns, which unveiled overlooked biases associated with long read sequencing. We next show that AISO around individual exons is remarkably constant across cell and tissue types and even upon major spliceosomal disruption, and it is evolutionarily conserved between human and mouse brains. We also establish a set of universal features associated with AISO patterns across various animal and plant species. Finally, we used Insplico to investigate AISO in the context of tissue-specific exons, particularly focusing on SRRM4-dependent microexons. We found that the majority of such microexons have non-canonical AISO, in which the downstream intron is spliced first, and we suggest two potential modes of SRRM4 regulation of microexons related to their AISO and various splicing-related features. Insplico is available on gitlab.com/aghr/insplico.

Ageing-associated changes in transcriptional elongation influence longevity

Physiological homeostasis becomes compromised during ageing, as a result of impairment of cellular processes, including transcription and RNA splicing1-4. However, the molecular mechanisms leading to the loss of transcriptional fidelity are so far elusive, as are ways of preventing it. Here we profiled and analysed genome-wide, ageing-related changes in transcriptional processes across different organisms: nematodes, fruitflies, mice, rats and humans. The average transcriptional elongation speed (RNA polymerase II speed) increased with age in all five species. Along with these changes in elongation speed, we observed changes in splicing, including a reduction of unspliced transcripts and the formation of more circular RNAs. Two lifespan-extending interventions, dietary restriction and lowered insulin-IGF signalling, both reversed most of these ageing-related changes. Genetic variants in RNA polymerase II that reduced its speed in worms5 and flies6 increased their lifespan. Similarly, reducing the speed of RNA polymerase II by overexpressing histone components, to counter age-associated changes in nucleosome positioning, also extended lifespan in flies and the division potential of human cells. Our findings uncover fundamental molecular mechanisms underlying animal ageing and lifespan-extending interventions, and point to possible preventive measures.

Integrating transcription and splicing into cell fate: Transcription factors on the block

Transcription factors (TFs) are present in all life forms and conserved across great evolutionary distances in eukaryotes. From yeast to complex multicellular organisms, they are pivotal players of cell fate decision by orchestrating gene expression at diverse molecular layers. Notably, TFs fine-tune gene expression by coordinating RNA fate at both the expression and splicing levels. They regulate alternative splicing, an essential mechanism for cell plasticity, allowing the production of many mRNA and protein isoforms in precise cell and tissue contexts. Despite this apparent role in splicing, how TFs integrate transcription and splicing to ultimately orchestrate diverse cell functions and cell fate decisions remains puzzling. We depict substantial studies in various model organisms underlining the key role of TFs in alternative splicing for promoting tissue-specific functions and cell fate. Furthermore, we emphasize recent advances describing the molecular link between the transcriptional and splicing activities of TFs. As TFs can bind both DNA and/or RNA to regulate transcription and splicing, we further discuss their flexibility and compatibility for DNA and RNA substrates. Finally, we propose several models integrating transcription and splicing activities of TFs in the coordination and diversification of cell and tissue identities. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications RNA Processing > Splicing Mechanisms.

HIV-1 Gag colocalizes with euchromatin histone marks at the nuclear periphery

The retroviral Gag protein of human immunodeficiency virus type 1 (HIV-1) plays a central role in the selection of unspliced viral genomic RNA for packaging into new virions. Previously, we demonstrated that full-length HIV-1 Gag undergoes nuclear trafficking where it associates with unspliced viral RNA (vRNA) at transcription sites. To further explore the kinetics of HIV-1 Gag nuclear localization, we used biochemical and imaging techniques to examine the timing of HIV-1 entry into the nucleus. We also aimed to determine more precisely Gag's subnuclear distribution to test the hypothesis that Gag would be associated with euchromatin, the transcriptionally active region of the nucleus. We observed that HIV-1 Gag localized to the nucleus shortly after its synthesis in the cytoplasm, suggesting that nuclear trafficking was not strictly concentration-dependent. Furthermore, we found that HIV-1 Gag preferentially localized to the transcriptionally active euchromatin fraction compared to the heterochromatin-rich region in a latently-infected CD4+ T cell line (J-Lat 10.6) treated with latency-reversal agents. Interestingly, HIV-1 Gag was more closely associated with transcriptionally-active histone markers near the nuclear periphery, where the HIV-1 provirus was previously shown to integrate. Although the precise function of Gag's association with histones in transcriptionally-active chromatin remains uncertain, together with previous reports, this finding is consistent with a potential role for euchromatin-associated Gag molecules to select newly transcribed unspliced vRNA during the initial stage of virion assembly.

Importance

The traditional view of retroviral assembly posits that HIV-1 Gag selection of unspliced vRNA begins in the cytoplasm. However, our previous studies demonstrated that HIV-1 Gag enters the nucleus and binds to unspliced HIV-1 RNA at transcription sites, suggesting that genomic RNA selection may occur in the nucleus. In the present study, we observed nuclear entry of HIV-1 Gag and co-localization with unspliced viral RNA within 8 hours post-expression. In CD4+ T cells (J-Lat 10.6) treated with latency reversal agents, as well as a HeLa cell line stably expressing an inducible Rev-dependent provirus, we found that HIV-1 Gag preferentially localized with histone marks associated with enhancer and promoter regions of transcriptionally active euchromatin near the nuclear periphery, which favors HIV-1 proviral integration sites. These observations support the hypothesis that HIV-1 Gag hijacks euchromatin-associated histones to localize to active transcription sites, promoting capture of newly synthesized genomic RNA for packaging.

Activity-Dependent Differential Regulation of Auts2 Isoforms In Vitro and In Vivo

Autism spectrum disorder (ASD) is a neurodevelopmental disorder of unknown cause, although one hypothesis suggests a potential imbalance between excitation and inhibition that leads to changes in neuronal activity and a disturbance in the brain network. However, the mechanisms through which neuronal activity contributes to the development of ASD remain largely unexplained. In this study, we described that neuronal activity at the transcriptional and translational levels regulated the expression of Auts2 isoforms. The prolonged stimulation of cultured cortical neurons significantly reduced the auts2 transcripts, accompanied by the decrease of FL-Auts2 protein, as well as one of the short isoforms (S-Auts2 var.1). Blocking neuronal activity increased the number of auts2 transcripts but not protein levels. Furthermore, blocking the NMDA receptors during stimulation could partially restore the FL-Auts2 and S-Auts2 var.1 at protein level, but not at mRNA level. Finally, Auts2 expression in the hippocampus was reduced in mice exposed to an enriched environment, a behavior paradigm designed to increase the brain activity through abundant sensory and social stimulations. Thus, our study revealed a novel regulatory effect of neuronal activity on the transcription and translation of ASD-risk gene auts2.

PTBP1-activated co-transcriptional splicing controls epigenetic status of pluripotent stem cells

Many spliceosomal introns are excised from nascent transcripts emerging from RNA polymerase II (RNA Pol II). The extent of cell-type-specific regulation and possible functions of such co-transcriptional events remain poorly understood. We examined the role of the RNA-binding protein PTBP1 in this process using an acute depletion approach followed by the analysis of chromatin- and RNA Pol II-associated transcripts. We show that PTBP1 activates the co-transcriptional excision of hundreds of introns, a surprising effect given that this protein is known to promote intron retention. Importantly, some co-transcriptionally activated introns fail to complete their splicing without PTBP1. In a striking example, retention of a PTBP1-dependent intron triggers nonsense-mediated decay of transcripts encoding DNA methyltransferase DNMT3B. We provide evidence that this regulation facilitates the natural decline in DNMT3B levels in developing neurons and protects differentiation-specific genes from ectopic methylation. Thus, PTBP1-activated co-transcriptional splicing is a widespread phenomenon mediating epigenetic control of cellular identity.

3D genome organization around nuclear speckles drives mRNA splicing efficiency

The nucleus is highly organized such that factors involved in transcription and processing of distinct classes of RNA are organized within specific nuclear bodies. One such nuclear body is the nuclear speckle, which is defined by high concentrations of protein and non-coding RNA regulators of pre-mRNA splicing. What functional role, if any, speckles might play in the process of mRNA splicing remains unknown. Here we show that genes localized near nuclear speckles display higher spliceosome concentrations, increased spliceosome binding to their pre-mRNAs, and higher co-transcriptional splicing levels relative to genes that are located farther from nuclear speckles. We show that directed recruitment of a pre-mRNA to nuclear speckles is sufficient to drive increased mRNA splicing levels. Finally, we show that gene organization around nuclear speckles is highly dynamic with differential localization between cell types corresponding to differences in Pol II occupancy. Together, our results integrate the longstanding observations of nuclear speckles with the biochemistry of mRNA splicing and demonstrate a critical role for dynamic 3D spatial organization of genomic DNA in driving spliceosome concentrations and controlling the efficiency of mRNA splicing.

A single-nucleus transcriptomics study of alcohol use disorder in the nucleus accumbens

Gene expression studies offer promising opportunities to better understand the processes underlying alcohol use disorder (AUD). As cell types differ in their function, gene expression profiles will typically vary across cell types. When studying bulk tissue, failure to account for this cellular diversity has a detrimental impact on the ability to detect disease associations. We therefore assayed the transcriptomes of 32,531 individual nuclei extracted from the nucleus accumbens (NAc) of nine donors with AUD and nine controls (72% male). Our study identified 17 clearly delineated cell types. We detected 26 transcriptome-wide significant differentially expressed genes (DEGs) that mainly involved medium spiny neurons with both D1-type and D2-type dopamine receptors, microglia (MGL) and oligodendrocytes. A higher than expected number of DEGs replicated in an existing single nucleus gene expression study of alcohol dependence in the prefrontal cortex (enrichment ratio 1.91, p value 0.019) with two genes remaining significant after a Bonferroni correction. Our most compelling result involved CD53 in MGL that replicated in the same cell type in the prefrontal cortex and was previously implicated in studies of DNA methylation, bulk gene expression and genetic variants. Several DEGs were previously reported to be associated with AUD (e.g., PER1 and MGAT5). The DEGs for MSN.3 seemed involved in neurodegeneration, disruption of circadian rhythms, alterations in glucose metabolism and changes in synaptic plasticity. For MGL, the DEGs implicated neuroinflammation and immune-related processes and for OLI, disruptions in myelination. This identification of the specific cell-types from which the association signals originate will be key for designing proper follow-up experiments and, eventually, novel clinical interventions.

Rapid nuclear deadenylation of mammalian messenger RNA

Poly(A) tails protect RNAs from degradation and their deadenylation rates determine RNA stability. Although poly(A) tails are generated in the nucleus, deadenylation of tails has mostly been investigated within the cytoplasm. Here, we combined long-read sequencing with metabolic labeling, splicing inhibition and cell fractionation experiments to quantify, separately, the genesis and trimming of nuclear and cytoplasmic tails in vitro and in vivo. We present evidence for genome-wide, nuclear synthesis of tails longer than 200 nt, which are rapidly shortened after transcription. Our data suggests that rapid deadenylation is a nuclear process, and that different classes of transcripts and even transcript isoforms have distinct nuclear tail lengths. For example, many long-noncoding RNAs retain long poly(A) tails. Modeling deadenylation dynamics predicts nuclear deadenylation about 10 times faster than cytoplasmic deadenylation. In summary, our data suggests that nuclear deadenylation might be a key mechanism for regulating mRNA stability, abundance, and subcellular localization.

Retained introns in long RNA-seq reads are not reliably detected in sample-matched short reads

BACKGROUND

There is growing interest in retained introns in a variety of disease contexts including cancer and aging. Many software tools have been developed to detect retained introns from short RNA-seq reads, but reliable detection is complicated by overlapping genes and transcripts as well as the presence of unprocessed or partially processed RNAs.

RESULTS

We compared introns detected by 8 tools using short RNA-seq reads with introns observed in long RNA-seq reads from the same biological specimens. We found significant disagreement among tools (Fleiss' [Formula: see text]) such that 47.7% of all detected intron retentions were not called by more than one tool. We also observed poor performance of all tools, with none achieving an F1-score greater than 0.26, and qualitatively different behaviors between general-purpose alternative splicing detection tools and tools confined to retained intron detection.

CONCLUSIONS

Short-read tools detect intron retention with poor recall and precision, calling into question the completeness and validity of a large percentage of putatively retained introns called by commonly used methods.

Relationships between genome-wide R-loop distribution and classes of recurrent DNA breaks in neural stem/progenitor cells

Recent studies revealed classes of recurrent DNA double-strand breaks (DSBs) in neural stem/progenitor cells, including transcription-associated, promoter-proximal breaks and recurrent DSB clusters in late-replicating, long neural genes that may give rise to somatic brain mosaicism. The mechanistic factors promoting these different classes of DSBs in neural stem/progenitor cells are not understood. Here, we elucidated the genome-wide landscape of RNA:DNA hybrid structures called “R-loops” in primary neural stem/progenitor cells undergoing aphidicolin-induced, mild replication stress to assess the potential contribution of R-loops to the different, recurrent classes of DNA break “hotspots”. We find that R-loops in neural stem/progenitor cells undergoing mild replication stress are present primarily in early-replicating, transcribed regions and in genes with promoter GC skew that are associated with cell lineage-specific processes. Surprisingly, most long, neural genes that form recurrent DSB clusters do not show R-loop formation under conditions of mild replication stress. Our findings are consistent with a role of R-loop-associated processes in promoter-proximal DNA break formation in highly transcribed, early replicating regions but suggest that R-loops do not drive replication stress-induced, recurrent DSB cluster formation in most long, neural genes.

Comprehensive analysis of the circadian nuclear and cytoplasmic transcriptome in mouse liver

In eukaryotes, RNA is synthesised in the nucleus, spliced, and exported to the cytoplasm where it is translated and finally degraded. Any of these steps could be subject to temporal regulation during the circadian cycle, resulting in daily fluctuations of RNA accumulation and affecting the distribution of transcripts in different subcellular compartments. Our study analysed the nuclear and cytoplasmic, poly(A) and total transcriptomes of mouse livers collected over the course of a day. These data provide a genome-wide temporal inventory of enrichment in subcellular RNA, and revealed specific signatures of splicing, nuclear export and cytoplasmic mRNA stability related to transcript and gene lengths. Combined with a mathematical model describing rhythmic RNA profiles, we could test the rhythmicity of export rates and cytoplasmic degradation rates of approximately 1400 genes. With nuclear export times usually much shorter than cytoplasmic half-lives, we found that nuclear export contributes to the modulation and generation of rhythmic profiles of 10% of the cycling nuclear mRNAs. This study contributes to a better understanding of the dynamic regulation of the transcriptome during the day-night cycle.

CDK12 regulates co-transcriptional splicing and RNA turnover in human cells

The cyclin-dependent kinase CDK12 has garnered interest as a cancer therapeutic target as DNA damage response genes are particularly suppressed by loss of CDK12 activity. In this study, we assessed the acute effects of CDK12 inhibition on transcription and RNA processing using nascent RNA Bru-seq and BruChase-seq. Acute transcriptional changes were overall small after CDK12 inhibition but over 600 genes showed intragenic premature termination, including DNA repair and cell cycle genes. Furthermore, many genes showed reduced transcriptional readthrough past the end of genes in the absence of CDK12 activity. RNA turnover was dramatically affected by CDK12 inhibition and importantly, caused increased degradation of many transcripts from DNA damage response genes. We also show that co-transcriptional splicing was suppressed by CDK12 inhibition. Taken together, these studies reveal the roles of CDK12 in regulating transcription elongation, transcription termination, co-transcriptional splicing, and RNA turnover.

•

Over 600 genes showed prematurely terminated transcription when CDK12 was inhibited

•

CDK12 promotes transcriptional readthrough past transcription end sites (TESs)

•

CDK12 promotes splicing and affects transcript stability

Modeling microRNA-driven post-transcriptional regulation using exon-intron split analysis in pigs

The contribution of microRNAs (miRNAs) to mRNA post-transcriptional regulation has often been explored by the post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes in mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Göttingen minipigs. Exon-intron split analysis of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of miRNAs upregulated in fed pigs. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, exon-intron split analysis provides a valuable view of the regulation of gene expression, complementary to canonical differential expression analyses.

Genomic sequences and RNA binding proteins predict RNA splicing efficiency in various single-cell contexts

MOTIVATION

The RNA splicing efficiency is of high interest for both understanding the regulatory machinery of gene expression and estimating the RNA velocity in single cells. However, its genomic regulation and stochasticity across contexts remain poorly understood.

RESULTS

Here, by leveraging the recent RNA velocity tool, we estimated the relative splicing efficiency across a variety of single-cell RNA-Seq data sets. We further extracted large sets of genomic features and 120 RNA binding protein features and found they are highly predictive to relative RNA splicing efficiency across multiple tissues and organs on human and mouse. This predictive power brings promise to reveal the complexity of RNA processing and to enhance the analysis of single-cell transcription activities.

AVAILABILITY AND IMPLEMENTATION

In order to ensure reproducibility, all preprocessed data sets and scripts used for the prediction and figure generation are publicly available at https://doi.org/10.5281/zenodo.6513669.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Co-transcriptional splicing efficiency is a gene-specific feature that can be regulated by TGFβ

Differential splicing efficiency of specific introns is a mechanism that dramatically increases protein diversity, based on selection of alternative exons for the final mature mRNA. However, it is unclear whether splicing efficiency of introns within the same gene is coordinated and eventually regulated as a mechanism to control mature mRNA levels. Based on nascent chromatin-associated RNA-sequencing data, we now find that co-transcriptional splicing (CTS) efficiency tends to be similar between the different introns of a gene. We establish that two well-differentiated strategies for CTS efficiency exist, at the extremes of a gradient: short genes that produce high levels of pre-mRNA undergo inefficient splicing, while long genes with relatively low levels of pre-mRNA have an efficient splicing. Notably, we observe that genes with efficient CTS display a higher level of mature mRNA relative to their pre-mRNA levels. Further, we show that the TGFβ signal transduction pathway regulates the general CTS efficiency, causing changes in mature mRNA levels. Taken together, our data indicate that CTS efficiency is a gene-specific characteristic that can be regulated to control gene expression.

Co-transcriptional splicing efficiency is a gene-specific characteristic that can be regulated by TGFβ to modulate gene expression.

Longitudinal intronic RNA-Seq analysis of Parkinson's disease patients reveals disease-specific nascent transcription

Transcriptomic studies usually focus on either gene or exon-based annotations, and only limited experiments have reported changes in reads mapping to introns. The analysis of intronic reads allows the detection of nascent transcription that is not influenced by steady-state RNA levels and provides information on actively transcribed genes. Here, we describe substantial intronic transcriptional changes in Parkinson's disease (PD) patients compared to healthy controls (CO) at two different timepoints; at the time of diagnosis (BL) and three years later (V08). We used blood RNA-Seq data from the Parkinson's Progression Markers Initiative (PPMI) cohort and identified significantly changed transcription of intronic reads only in PD patients during this follow-up period. In CO subjects, only nine transcripts demonstrated differentially expressed introns between visits. However, in PD patients, 4873 transcripts had differentially expressed introns at visit V08 compared to BL, many of them in genes previously associated with neurodegenerative diseases, such as LRRK2, C9orf72, LGALS3, KANSL1AS1, and ALS2. In addition, at the time of diagnosis (BL visit), we identified 836 transcripts (e.g. SNCA, DNAJC19, PRRG4) and at visit V08, 2184 transcripts (e.g. PINK1, GBA, ALS2, PLEKHM1) with differential intronic expression specific to PD patients. In contrast, reads mapping to exonic regions demonstrated little variation indicating highly specific changes only in intronic transcription. Our study demonstrated that PD is characterized by substantial changes in the nascent transcription, and description of these changes could help to understand the molecular pathology underpinning this disease.

Genetic variants associated mRNA stability in lung

Background

Expression quantitative trait loci (eQTLs) analyses have been widely used to identify genetic variants associated with gene expression levels to understand what molecular mechanisms underlie genetic traits. The resultant eQTLs might affect the expression of associated genes through transcriptional or post-transcriptional regulation. In this study, we attempt to distinguish these two types of regulation by identifying genetic variants associated with mRNA stability of genes (stQTLs).

Results

Here, we presented a computational framework that takes advantage of recently developed methods to infer the mRNA stability of genes based on RNA-seq data and performed association analysis to identify stQTLs. Using the Genotype-Tissue Expression (GTEx) lung RNA-Seq data, we identified a total of 142,801 stQTLs for 3942 genes and 186,132 eQTLs for 4751 genes from 15,122,700 genetic variants for 13,476 genes on the autosomes, respectively. Interestingly, our results indicated that stQTLs were enriched in the CDS and 3’UTR regions, while eQTLs are enriched in the CDS, 3’UTR, 5’UTR, and upstream regions. We also found that stQTLs are more likely than eQTLs to overlap with RNA binding protein (RBP) and microRNA (miRNA) binding sites. Our analyses demonstrate that simultaneous identification of stQTLs and eQTLs can provide more mechanistic insight on the association between genetic variants and gene expression levels.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08405-y.

Non-Canonical Splicing and Its Implications in Brain Physiology and Cancer

The advance of experimental and computational techniques has allowed us to highlight the existence of numerous different mechanisms of RNA maturation, which have been so far unknown. Besides canonical splicing, consisting of the removal of introns from pre-mRNA molecules, non-canonical splicing events may occur to further increase the regulatory and coding potential of the human genome. Among these, splicing of microexons, recursive splicing and biogenesis of circular and chimeric RNAs through back-splicing and trans-splicing processes, respectively, all contribute to expanding the repertoire of RNA transcripts with newly acquired regulatory functions. Interestingly, these non-canonical splicing events seem to occur more frequently in the central nervous system, affecting neuronal development and differentiation programs with important implications on brain physiology. Coherently, dysregulation of non-canonical RNA processing events is associated with brain disorders, including brain tumours. Herein, we summarize the current knowledge on molecular and regulatory mechanisms underlying canonical and non-canonical splicing events with particular emphasis on cis-acting elements and trans-acting factors that all together orchestrate splicing catalysis reactions and decisions. Lastly, we review the impact of non-canonical splicing on brain physiology and pathology and how unconventional splicing mechanisms may be targeted or exploited for novel therapeutic strategies in cancer.

The SWI/SNF subunit BRG1 affects alternative splicing by changing RNA binding factor interactions with nascent RNA

BRG1 and BRM are ATPase core subunits of the human SWI/SNF chromatin remodelling complexes mainly associated with transcriptional initiation. They also have a role in alternative splicing, which has been shown for BRM-containing SWI/SNF complexes at a few genes. Here, we have identified a subset of genes which harbour alternative exons that are affected by SWI/SNF ATPases by expressing the ATPases BRG1 and BRM in C33A cells, a BRG1- and BRM-deficient cell line, and analysed the effect on splicing by RNA sequencing. BRG1- and BRM-affected sub-sets of genes favouring both exon inclusion and exon skipping, with only a minor overlap between the ATPase. Some of the changes in alternative splicing induced by BRG1 and BRM expression did not require the ATPase activity. The BRG1-ATPase independent included exons displayed an exon signature of a high GC content. By investigating three genes with exons affected by the BRG-ATPase-deficient variant, we show that these exons accumulated phosphorylated RNA pol II CTD, both serine 2 and serine 5 phosphorylation, without an enrichment of the RNA polymerase II. The ATPases were recruited to the alternative exons, together with both core and signature subunits of SWI/SNF complexes, and promoted the binding of RNA binding factors to chromatin and RNA at the alternative exons. The interaction with the nascent RNP, however, did not reflect the association to chromatin. The hnRNPL, hnRNPU and SAM68 proteins associated with chromatin in cells expressing BRG1 and BRM wild type, but the binding of hnRNPU to the nascent RNP was excluded. This suggests that SWI/SNF can regulate alternative splicing by interacting with splicing-RNA binding factor and influence their binding to the nascent pre-mRNA particle.

Characterization of transcript enrichment and detection bias in single-nucleus RNA-seq for mapping of distinct human adipocyte lineages

Single-cell RNA sequencing (scRNA-seq) enables molecular characterization of complex biological tissues at high resolution. The requirement of single-cell extraction, however, makes it challenging for profiling tissues such as adipose tissue, for which collection of intact single adipocytes is complicated by their fragile nature. For such tissues, single-nucleus extraction is often much more efficient and therefore single-nucleus RNA sequencing (snRNA-seq) presents an alternative to scRNA-seq. However, nuclear transcripts represent only a fraction of the transcriptome in a single cell, with snRNA-seq marked with inherent transcript enrichment and detection biases. Therefore, snRNA-seq may be inadequate for mapping important transcriptional signatures in adipose tissue. In this study, we compare the transcriptomic landscape of single nuclei isolated from preadipocytes and mature adipocytes across human white and brown adipocyte lineages, with whole-cell transcriptome. We show that snRNA-seq is capable of identifying the broad cell types present in scRNA-seq at all states of adipogenesis. However, we also explore how and why the nuclear transcriptome is biased and limited, as well as how it can be advantageous. We robustly characterize the enrichment of nuclear-localized transcripts and adipogenic regulatory lncRNAs in snRNA-seq, while also providing a detailed understanding for the preferential detection of long genes upon using this technique. To remove such technical detection biases, we propose a normalization strategy for a more accurate comparison of nuclear and cellular data. Finally, we show successful integration of scRNA-seq and snRNA-seq data sets with existing bioinformatic tools. Overall, our results illustrate the applicability of snRNA-seq for the characterization of cellular diversity in the adipose tissue.

Transcription and splicing dynamics during early Drosophila development

Widespread cotranscriptional splicing has been demonstrated from yeast to human. However, most studies to date addressing the kinetics of splicing relative to transcription used either Saccharomyces cerevisiae or metazoan cultured cell lines. Here, we adapted native elongating transcript sequencing technology (NET-seq) to measure cotranscriptional splicing dynamics during the early developmental stages of Drosophila melanogaster embryos. Our results reveal the position of RNA polymerase II (Pol II) when both canonical and recursive splicing occur. We found heterogeneity in splicing dynamics, with some RNAs spliced immediately after intron transcription, whereas for other transcripts no splicing was observed over the first 100 nt of the downstream exon. Introns that show splicing completion before Pol II has reached the end of the downstream exon are necessarily intron-defined. We studied the splicing dynamics of both nascent pre-mRNAs transcribed in the early embryo, which have few and short introns, as well as pre-mRNAs transcribed later in embryonic development, which contain multiple long introns. As expected, we found a relationship between the proportion of spliced reads and intron size. However, intron definition was observed at all intron sizes. We further observed that genes transcribed in the early embryo tend to be isolated in the genome whereas genes transcribed later are often overlapped by a neighboring convergent gene. In isolated genes, transcription termination occurred soon after the polyadenylation site, while in overlapped genes, Pol II persisted associated with the DNA template after cleavage and polyadenylation of the nascent transcript. Taken together, our data unravel novel dynamic features of Pol II transcription and splicing in the developing Drosophila embryo.

Recursive splicing is a rare event in the mouse brain

Recursive splicing (RS) is a splicing mechanism to remove long introns from messenger RNA precursors of long genes. Compared to the hundreds of RS events identified in humans and drosophila, only ten RS events have been reported in mice. To further investigate RS in mice, we analyzed RS in the mouse brain, a tissue that is enriched in the expression of long genes. We found that nuclear total RNA sequencing is an efficient approach to investigate RS events. We analyzed 1.15 billion uniquely mapped reads from the nuclear total RNA sequencing data in the mouse cerebral cortex. Unexpectedly, we only identified 20 RS sites, suggesting that RS is a rare event in the mouse brain. We also identified that RS is constitutive between excitatory and inhibitory neurons and between sexes in the mouse cerebral cortex. In addition, we found that the primary sequence context is associated with RS splicing intermediates and distinguishes RS AGGT site from non-RS AGGT sites, indicating the importance of the primary sequence context in RS sites. Moreover, we discovered that cryptic exons may use an RS-like mechanism for splicing. Overall, we provide novel findings about RS in long genes in the mouse brain.

Profiling transcription factor activity dynamics using intronic reads in time-series transcriptome data

Activities of transcription factors (TFs) are temporally modulated to regulate dynamic cellular processes, including development, homeostasis, and disease. Recent developments of bioinformatic tools have enabled the analysis of TF activities using transcriptome data. However, because these methods typically use exon-based target expression levels, the estimated TF activities have limited temporal accuracy. To address this, we proposed a TF activity measure based on intron-level information in time-series RNA-seq data, and implemented it to decode the temporal control of TF activities during dynamic processes. We showed that TF activities inferred from intronic reads can better recapitulate instantaneous TF activities compared to the exon-based measure. By analyzing public and our own time-series transcriptome data, we found that intron-based TF activities improve the characterization of temporal phasing of cycling TFs during circadian rhythm, and facilitate the discovery of two temporally opposing TF modules during T cell activation. Collectively, we anticipate that the proposed approach would be broadly applicable for decoding global transcriptional architecture during dynamic processes.

Using Drosophila to uncover molecular and physiological functions of circRNAs

Circular RNAs (circRNAs) are a class of covalently closed RNA molecules generated by backsplicing. circRNAs are expressed in a tissue-specific manner, accumulate with age in neural tissues, and are highly stable. In many cases, circRNAs are generated at the expense of a linear transcript as back-splicing competes with linear splicing. Some circRNAs regulate gene expression in cis, and some circRNAs can be translated into protein. The advent of deep sequencing and new bioinformatic tools has allowed detection of thousands of circRNAs in eukaryotes. Studying the functions of circRNAs is done using a combination of molecular and genetic methods. The unique genetic tools that can be used in studies of Drosophila melanogaster are ideal for deciphering the functions of circRNAs in vivo. These tools include the GAL4-UAS system, which can be used to manipulate the levels of circRNAs with exquisite temporal and spatial control, and genetic interaction screening, which could be used to identify pathways regulated by circRNAs. Research performed in Drosophila has revealed circRNAs production mechanisms, details of their translation, and their physiological functions. Due to their short lifecycle and the existence of excellent neurodegeneration models, Drosophila can also be used to study the role of circRNAs in aging and age-related disorders. Here, we review molecular and genetic tools and methods for detecting, manipulating, and studying circRNAs in Drosophila.

Full-length transcript sequencing of human and mouse cerebral cortex identifies widespread isoform diversity and alternative splicing

Alternative splicing is a post-transcriptional regulatory mechanism producing distinct mRNA molecules from a single pre-mRNA with a prominent role in the development and function of the central nervous system. We used long-read isoform sequencing to generate full-length transcript sequences in the human and mouse cortex. We identify novel transcripts not present in existing genome annotations, including transcripts mapping to putative novel (unannotated) genes and fusion transcripts incorporating exons from multiple genes. Global patterns of transcript diversity are similar between human and mouse cortex, although certain genes are characterized by striking differences between species. We also identify developmental changes in alternative splicing, with differential transcript usage between human fetal and adult cortex. Our data confirm the importance of alternative splicing in the cortex, dramatically increasing transcriptional diversity and representing an important mechanism underpinning gene regulation in the brain. We provide transcript-level data for human and mouse cortex as a resource to the scientific community.

•

There is widespread transcript diversity in the cortex and many novel transcripts

•

Some genes display big differences in isoform number between human and mouse cortex

•

There is evidence of differential transcript usage between human fetal and adult cortex

•

There are many novel isoforms of genes associated with human brain disease

Conditional depletion of transcriptional kinases Ctk1 and Bur1 and effects on co-transcriptional spliceosome assembly and pre-mRNA splicing

From yeast to humans, pre-mRNA splicing occurs mainly co-transcriptionally, with splicing and transcription functionally coupled such that they influence one another. The recruitment model of co-transcriptional splicing proposes that core members of the transcription elongation machinery have the potential to influence co-transcriptional spliceosome assembly and pre-mRNA splicing. Here, we tested whether the transcription elongation kinases Bur1 and Ctk1 affect co-transcriptional spliceosome assembly and pre-mRNA splicing in the budding yeast Saccharomyces cerevisiae. In S. cerevisiae, Ctk1 is the major kinase that phosphorylates serine 2 of the carboxy-terminal domain of the largest subunit of RNA polymerase II, whilst Bur1 augments the kinase activity of Ctk1 and is the major kinase for elongation factor Spt5. We used the auxin-inducible degron system to conditionally deplete Bur1 and Ctk1 kinases, and investigated the effects on co-transcriptional spliceosome assembly and pre-mRNA splicing. Depletion of Ctk1 effectively reduced phosphorylation of serine 2 of the carboxy-terminal domain but did not impact co-transcriptional spliceosome assembly or pre-mRNA splicing. In striking contrast, depletion of Bur1 did not reduce phosphorylation of serine 2 of the carboxy-terminal domain, but reduced Spt5 phosphorylation and enhanced co-transcriptional spliceosome assembly and pre-mRNA splicing, suggesting a role for this kinase in modulating co-transcriptional splicing.

Systematic functional characterization of antisense eRNA of protocadherin α composite enhancer

Enhancers generate bidirectional noncoding enhancer RNAs (eRNAs) that may regulate gene expression. At present, the eRNA function remains enigmatic. Here, we report a 5' capped antisense eRNA PEARL (Pcdh eRNA associated with R-loop formation) that is transcribed from the protocadherin (Pcdh) α HS5-1 enhancer region. Through loss- and gain-of-function experiments with CRISPR/Cas9 DNA fragment editing, CRISPRi, and CRISPRa, as well as locked nucleic acid strategies, in conjunction with ChIRP, MeDIP, DRIP, QHR-4C, and HiChIP experiments, we found that PEARL regulates Pcdhα gene expression by forming local RNA-DNA duplexes (R-loops) in situ within the HS5-1 enhancer region to promote long-distance chromatin interactions between distal enhancers and target promoters. In particular, increased levels of eRNA PEARL via perturbing transcription elongation factor SPT6 lead to strengthened local three-dimensional chromatin organization within the Pcdh superTAD. These findings have important implications regarding molecular mechanisms by which the HS5-1 enhancer regulates stochastic Pcdhα promoter choice in single cells in the brain.

Detection of Alu Exonization Events in Human Frontal Cortex From RNA-Seq Data

Alu exonization events functionally diversify the transcriptome, creating alternative mRNA isoforms and accounting for an estimated 5% of the alternatively spliced (skipped) exons in the human genome. We developed computational methods, implemented into a software called Alubaster, for detecting incorporation of Alu sequences in mRNA transcripts from large scale RNA-seq data sets. The approach detects Alu sequences derived from both fixed and polymorphic Alu elements, including Alu insertions missing from the reference genome. We applied our methods to 117 GTEx human frontal cortex samples to build and characterize a collection of Alu-containing mRNAs. In particular, we detected and characterized Alu exonizations occurring at 870 fixed Alu loci, of which 237 were novel, as well as hundreds of putative events involving Alu elements that are polymorphic variants or rare alleles not present in the reference genome. These methods and annotations represent a unique and valuable resource that can be used to understand the characteristics of Alu-containing mRNAs and their tissue-specific expression patterns.

The DNA-to-cytoplasm ratio broadly activates zygotic gene expression in Xenopus

In multicellular animals, the first major event after fertilization is the switch from maternal to zygotic control of development. During this transition, zygotic gene transcription is broadly activated in an otherwise quiescent genome in a process known as zygotic genome activation (ZGA). In fast-developing embryos, ZGA often overlaps with the slowing of initially synchronous cell divisions at the mid-blastula transition (MBT). Initial studies of the MBT led to the nuclear-to-cytoplasmic ratio model where MBT timing is regulated by the exponentially increasing amounts of some nuclear component "N" titrated against a fixed cytoplasmic component "C." However, more recent experiments have been interpreted to suggest that ZGA is independent of the N/C ratio. To determine the role of the N/C ratio in ZGA, we generated Xenopus frog embryos with ∼3-fold differences in genomic DNA (i.e., N) by using X. tropicalis sperm to fertilize X. laevis eggs with or without their maternal genome. Resulting embryos have otherwise identical X. tropicalis genome template amounts, embryo sizes, and X. laevis maternal environments. We generated transcriptomic time series across the MBT in both conditions and used X. tropicalis paternally derived mRNA to identify a high-confidence set of exclusively zygotic transcripts. Both ZGA and the increase in cell-cycle duration are delayed in embryos with ∼3-fold less DNA per cell. Thus, DNA is an important component of the N/C ratio, which is a critical regulator of zygotic genome activation in Xenopus embryos.