Mapping and modeling the genomic basis of differential RNA isoform expression at single-cell resolution with LR-Split-seq

The rise in throughput and quality of long-read sequencing should allow unambiguous identification of full-length transcript isoforms. However, its application to single-cell RNA-seq has been limited by throughput and expense. Here we develop and characterize long-read Split-seq (LR-Split-seq), which uses combinatorial barcoding to sequence single cells with long reads. Applied to the C2C12 myogenic system, LR-split-seq associates isoforms to cell types with relative economy and design flexibility. We find widespread evidence of changing isoform expression during differentiation including alternative transcription start sites (TSS) and/or alternative internal exon usage. LR-Split-seq provides an affordable method for identifying cluster-specific isoforms in single cells.

Nanopore sequencing of full-length circRNAs in human and mouse brains reveals circRNA-specific exon usage and intron retention

Circular RNA (circRNA) is a poorly understood class of non-coding RNAs, some of which have been shown to be functional important for cell proliferation and development. CircRNAs mainly derive from back splicing events of coding mRNAs, making it difficult to distinguish the internal exon composition of circRNA from the linearly spliced mRNA. To examine the global exon composition of circRNAs, we performed long-read sequencing of single molecules using nanopore technology for human and mouse brain-derived RNA. By applying an optimized circRNA enrichment protocol prior to sequencing, we were able to detect 7,834 and 10,975 circRNAs in human and mouse brain, respectively, of which 2,945 and 7,052 are not currently found in circBase. Alternative splicing was more prevalent in circRNAs than in linear spliced transcripts, and notably >200 not previously annotated exons were used in circRNAs. This suggests that properties associated with circRNA- specific features, e.g. the unusual back-splicing step during biogenesis, increased stability and /or their lack of translation, alter the general exon usage at steady state. We conclude that the nanopore sequencing technology provides a fast and reliable method to map the specific exon composition of circRNA.

Single-molecule imaging of transcription at damaged chromatin

How DNA double-strand breaks (DSBs) affect ongoing transcription remains elusive due to the lack of single-molecule resolution tools directly measuring transcription dynamics upon DNA damage. Here, we established new reporter systems that allow the visualization of individual nascent RNAs with high temporal and spatial resolution upon the controlled induction of a single DSB at two distinct chromatin locations: a promoter-proximal (PROP) region downstream the transcription start site and a region within an internal exon (EX2). Induction of a DSB resulted in a rapid suppression of preexisting transcription initiation regardless of the genomic location. However, while transcription was irreversibly suppressed upon a PROP DSB, damage at the EX2 region drove the formation of promoter-like nucleosome-depleted regions and transcription recovery. Two-color labeling of transcripts at sequences flanking the EX2 lesion revealed bidirectional break-induced transcription initiation. Transcriptome analysis further showed pervasive bidirectional transcription at endogenous intragenic DSBs. Our data provide a novel framework for interpreting the reciprocal interactions between transcription and DNA damage at distinct chromatin regions.

Cycles of gene expression and genome response during mammalian tissue regeneration

BackgroundCompensatory liver hyperplasia — or regeneration — induced by two-thirds partial hepatectomy (PH) permits the study of synchronized activation of mammalian gene expression, particularly in relation to cell proliferation. Here, we measured genomic transcriptional responses and mRNA accumulation changes after PH and sham surgeries.ResultsDuring the first 10–20 hours, the PH- and sham-surgery responses were very similar, including parallel early activation of cell-division-cycle genes. After 20 hours, however, whereas post-PH livers continued with a robust and coordinate cell-division-cycle gene-expression response before returning to the resting state by one week, sham-surgery livers returned directly to a resting gene-expression state. Localization of RNA polymerase II (Pol II), and trimethylated histone H3 lysine 4 (H3K4me3) and 36 (H3K36me3) on genes dormant in the resting liver and activated during the PH response revealed a general de novo promoter Pol II recruitment and H3K4me3 increase during the early 10–20 hour phase followed by Pol II elongation and H3K36me3 accumulation in gene bodies during the later proliferation phase. H3K36me3, generally appearing at the first-internal exon, was preceded 5′ by H3K36me2; 3′ of the first-internal exon, in about half of genes H3K36me3 predominated and in the other half H3K36me2 and H3K36me3 co-existed. Further, we observed some unusual gene profiles with abundant Pol II but little evident H3K4me3 or H3K36me3 modification, indicating that these modifications are neither universal nor essential partners to Pol II transcription.ConclusionsPH and sham surgical procedures on mice reveal striking early post-operatory gene expression similarities followed by synchronized mRNA accumulation and epigenetic histone mark changes specific to PH.

Protein 4.1R Exon 16 3′ Splice Site Activation Requires Coordination among TIA1, Pcbp1, and RBM39 during Terminal Erythropoiesis

ABSTRACT Exon 16 of protein 4.1R encodes a spectrin/actin-binding peptide critical for erythrocyte membrane stability. Its expression during erythroid differentiation is regulated by alternative pre-mRNA splicing. A UUUUCCCCCC motif situated between the branch point and the 3′ splice site is crucial for inclusion. We show that the UUUU region and the last three C residues in this motif are necessary for the binding of splicing factors TIA1 and Pcbp1 and that these proteins appear to act in a collaborative manner to enhance exon 16 inclusion. This element also activates an internal exon when placed in a corresponding intronic position in a heterologous reporter. The impact of these two factors is further enhanced by high levels of RBM39, whose expression rises during erythroid differentiation as exon 16 inclusion increases. TIA1 and Pcbp1 associate in a complex containing RBM39, which interacts with U2AF65 and SF3b155 and promotes U2 snRNP recruitment to the branch point. Our results provide a mechanism for exon 16 3′ splice site activation in which a coordinated effort among TIA1, Pcbp1, and RBM39 stabilizes or increases U2 snRNP recruitment, enhances spliceosome A complex formation, and facilitates exon definition through RBM39-mediated splicing regulation.

On the origins and evolution of trans-splicing of bursicon in mosquitos

Broad transcriptomic sequencing of eukaryotes has revealed the ubiquity of splicing of nuclear genes. While the vast majority of splicing events join segments of the same RNA transcript, various studies have found a few intriguing cases of trans-splicing of introns, in which splicing events within protein coding regions join segments of different RNA transcripts. The most structurally intricate case known involves the bursicon gene in mosquitos, in which an internal exon is encoded at a distinct locus, requiring multiple trans-splicing events form the mature mRNA. This arrangement is known to be ancestral to mosquitos, however the exact timing of the origin of trans-splicing and the history of the bursicon gene within mosquitos is unknown. Taking advantage of the recent availability of genomes from various Anopheles mosquitos and from relatives of mosquitos, I determined trans versus cis encoding of bursicon across Culicomorpha. I conclude that trans-splicing emerged in the last common ancestor of mosquitos, and that trans-splicing has been retained in all 19 studied Anopheles species. The retention of trans-splicing could indicate functional importance of this arrangement, or could alternatively reflect the rarity of mutations giving rise to viable allelic alternatives.