Quantitative RNA-seq meta-analysis of alternative exon usage in C. elegans

A comprehensive analysis of 3'UTRs in Caenorhabditis elegans

3'Untranslated regions (3'UTRs) are essential portions of genes containing elements necessary for pre-mRNA 3'end processing and are involved in post-transcriptional gene regulation. Despite their importance, they remain poorly characterized in eukaryotes. Here, we have used a multi-pronged approach to extract and curate 3'UTR data from 11533 publicly available datasets, corresponding to the entire collection of Caenorhabditis elegans transcriptomes stored in the NCBI repository from 2009 to 2023. We have also performed high throughput cloning pipelines to identify and validate rare 3'UTR isoforms and incorporated and manually curated 3'UTR isoforms from previously published datasets. This updated C. elegans 3'UTRome (v3) is the most comprehensive resource in any metazoan to date, covering 97.4% of the 20362 experimentally validated protein-coding genes with refined and updated 3'UTR boundaries for 23489 3'UTR isoforms. We also used this novel dataset to identify and characterize sequence elements involved in pre-mRNA 3'end processing and update miRNA target predictions. This resource provides important insights into the 3'UTR formation, function, and regulation in eukaryotes.

A nematode-specific ribonucleoprotein complex mediates interactions between the major nematode spliced leader snRNP and its target pre-mRNAs

Spliced leader trans-splicing of pre-mRNAs is a critical step in the gene expression of many eukaryotes. How the spliced leader RNA and its target transcripts are brought together to form the trans-spliceosome remains an important unanswered question. Using immunoprecipitation followed by protein analysis via mass spectrometry and RIP-Seq, we show that the nematode-specific proteins, SNA-3 and SUT-1, form a complex with a set of enigmatic non-coding RNAs, the SmY RNAs. Our work redefines the SmY snRNP and shows for the first time that it is essential for nematode viability and is involved in spliced leader trans-splicing. SNA-3 and SUT-1 are associated with the 5' ends of most, if not all, nascent capped RNA polymerase II transcripts, and they also interact with components of the major nematode spliced leader (SL1) snRNP. We show that depletion of SNA-3 impairs the co-immunoprecipitation between one of the SL1 snRNP components, SNA-2, and several core spliceosomal proteins. We thus propose that the SmY snRNP recruits the SL1 snRNP to the 5' ends of nascent pre-mRNAs, an instrumental step in the assembly of the trans-spliceosome.

Proteasome hyperactivation rewires the proteome enhancing stress resistance, proteostasis, lipid metabolism and ERAD in C. elegans

Proteasome dysfunction is implicated in the pathogenesis of neurodegenerative diseases and age-related proteinopathies. Using a C. elegans model, we demonstrate that 20S proteasome hyperactivation, facilitated by 20S gate-opening, accelerates the targeting of intrinsically disordered proteins. This leads to increased protein synthesis, extensive rewiring of the proteome and transcriptome, enhanced oxidative stress defense, accelerated lipid metabolism, and peroxisome proliferation. It also promotes ER-associated degradation (ERAD) of aggregation-prone proteins, such as alpha-1 antitrypsin (ATZ) and various lipoproteins. Notably, our results reveal that 20S proteasome hyperactivation suggests a novel role in ERAD with broad implications for proteostasis-related disorders, simultaneously affecting lipid homeostasis and peroxisome proliferation. Furthermore, the enhanced cellular capacity to mitigate proteostasis challenges, alongside unanticipated acceleration of lipid metabolism is expected to contribute to the longevity phenotype of this mutant. Remarkably, the mechanism of longevity induced by 20S gate opening appears unique, independent of known longevity and stress-resistance pathways. These results support the therapeutic potential of 20S proteasome activation in mitigating proteostasis-related disorders broadly and provide new insights into the complex interplay between proteasome activity, cellular health, and aging.

C. elegans Afadin is required for epidermal morphogenesis and functionally interfaces with the cadherin-catenin complex and RhoGAP PAC-1/ARHGAP21

During epithelial morphogenesis, the apical junctions connecting cells must remodel as cells change shape and make new connections with their neighbors. In the C. elegans embryo, new apical junctions form when epidermal cells migrate and seal with one another to encase the embryo in skin ('ventral enclosure'), and junctions remodel when epidermal cells change shape to squeeze the embryo into a worm shape ('elongation'). The junctional cadherin-catenin complex (CCC), which links epithelial cells to each other and to cortical actomyosin, is essential for C. elegans epidermal morphogenesis. RNAi genetic enhancement screens have identified several proteins that interact with the CCC to promote epidermal morphogenesis, including the scaffolding protein Afadin (AFD-1), whose depletion alone results in only minor morphogenesis defects. Here, by creating a null mutation in afd-1 , we show that afd-1 provides a significant contribution to ventral enclosure and elongation on its own. Unexpectedly, we find that afd-1 mutant phenotypes are strongly modified by diet, revealing a previously unappreciated maternal nutritional input to morphogenesis. We identify functional interactions between AFD-1 and the CCC by demonstrating that E-cadherin is required for the polarized distribution of AFD-1 to cell contact sites in early embryos. Finally, we show that afd-1 promotes the enrichment of polarity regulator and CCC-interacting protein PAC-1/ARHGAP21 to cell contact sites, and identify genetic interactions suggesting that afd-1 and pac-1 regulate epidermal morphogenesis at least in part through parallel mechanisms. Our findings reveal that C. elegans AFD-1 makes a significant contribution to epidermal morphogenesis and functionally interfaces with core and associated CCC proteins.

Curcumin supplementation increases longevity and antioxidant capacity in Caenorhabditis elegans

Curcumin is well known as a potent antioxidant and free radical scavenger and has great potential for anti-aging applications. In this study, we investigate the molecular mechanism of curcumin in prolonging the lifespan of C. elegans. Four concentrations of curcumin (10, 25, 50, and 100 µM) were administered, and the optimal treatment concentration was determined by analyzing the nematode lifespan, physiology, and biochemistry. Additionally, RNA-seq and qRT-PCR were performed to explore the antioxidant effect of curcumin and its underlying mechanism. Results revealed that curcumin could significantly improve the survival capacity of C. elegans without influencing its growth. Curcumin was observed to significantly decrease the levels of reactive oxygen species (ROS) under extreme conditions such as heat stress and paraquat stress. In addition, curcumin increased the amount of nematode mitochondrial DNA (mtDNA) replication. RNA-seq results revealed that the underlying mechanism of curcumin in C. elegans is related to the mitogen-activated protein kinase (MAPK) pathway. qRT-PCR results confirmed that the expression of oxidative stress-related genes (sod-1, sod-2, sod-3, gst-4) was increased, and the expression of MAPK signaling pathway-related genes (sek-1, pmk-1, nsy-1) was significantly downregulated. Furthermore, the administration of curcumin extended the lifespan of nematodes, potentially through the enhancement of oxidative stress resistance and the downregulation of the MAPK signaling pathway. These findings improve our understanding of both lifespan extension and the potential mechanism of curcumin in C. elegans.

Two RhoGEF isoforms with distinct localisation control furrow position during asymmetric cell division

Cytokinesis partitions cellular content between daughter cells. It relies on the formation of an acto-myosin contractile ring, whose constriction induces the ingression of the cleavage furrow between the segregated chromatids. Rho1 GTPase and its RhoGEF (Pbl) are essential for this process. However, how Rho1 is regulated to sustain furrow ingression while maintaining correct furrow position remains poorly defined. Here, we show that during asymmetric division of Drosophila neuroblasts, Rho1 is controlled by two Pbl isoforms with distinct localisation. Spindle midzone- and furrow-enriched Pbl-A focuses Rho1 at the furrow to sustain efficient ingression, while Pbl-B pan-plasma membrane localization promotes the broadening of Rho1 activity and the subsequent enrichment of myosin on the entire cortex. This enlarged zone of Rho1 activity is critical to adjust furrow position, thereby preserving correct daughter cell size asymmetry. Our work highlights how the use of isoforms with distinct localisation makes an essential process more robust.

Quantitative analysis of C. elegans transcripts by Nanopore direct-cDNA sequencing reveals terminal hairpins in non trans-spliced mRNAs

In nematodes and kinetoplastids, mRNA processing involves a trans-splicing step through which a short sequence from a snRNP replaces the original 5' end of the primary transcript. It has long been held that 70% of C. elegans mRNAs are submitted to trans-splicing. Our recent work suggested that the mechanism is more pervasive but not fully captured by mainstream transcriptome sequencing methods. Here we use Oxford Nanopore's long-read amplification-free sequencing technology to perform a comprehensive analysis of trans-splicing in worms. We demonstrate that spliced leader (SL) sequences at the 5' end of the mRNAs affect library preparation and generate sequencing artefacts due to their self-complementarity. Consistent with our previous observations, we find evidence of trans-splicing for most genes. However, a subset of genes appears to be only marginally trans-spliced. These mRNAs all share the capacity to generate a 5' terminal hairpin structure mimicking the SL structure and offering a mechanistic explanation for their non conformity. Altogether, our data provide a comprehensive quantitative analysis of SL usage in C. elegans.

The ubiquitin-like protein Hub1/UBL-5 functions in pre-mRNA splicing in Caenorhabditis elegans

The ubiquitin-like protein Hub1/UBL-5 associates with proteins non-covalently. Hub1 promotes alternative splicing and splicing of precursor mRNAs with weak introns in yeast and mammalian cells; however, its splicing function has remained elusive in multicellular organisms. Here, we demonstrate the splicing function of Hub1/UBL-5 in the free-living nematode Caenorhabditis elegans. Hub1/UBL-5 binds to the HIND-containing splicing factors Snu66/SART-1 and PRP-38 and associates with other spliceosomal proteins. C. elegans hub1/ubl-5 mutants die at the Larval 3 stage and show splicing defects for selected targets, similar to the mutants in yeast and mammalian cells. UBL-5 complemented growth and splicing defects in Schizosaccharomyces pombe hub1 mutants, confirming its functional conservation. Thus, UBL-5 is important for C. elegans development and plays a conserved pre-mRNA splicing function.

A novel, essential trans-splicing protein connects the nematode SL1 snRNP to the CBC-ARS2 complex

Spliced leader trans-splicing is essential for gene expression in many eukaryotes. To elucidate the molecular mechanism of this process, we characterise the molecules associated with the Caenorhabditis elegans major spliced leader snRNP (SL1 snRNP), which donates the spliced leader that replaces the 5' untranslated region of most pre-mRNAs. Using a GFP-tagged version of the SL1 snRNP protein SNA-1 created by CRISPR-mediated genome engineering, we immunoprecipitate and identify RNAs and protein components by RIP-Seq and mass spectrometry. This reveals the composition of the SL1 snRNP and identifies associations with spliceosome components PRP-8 and PRP-19. Significantly, we identify a novel, nematode-specific protein required for SL1 trans-splicing, which we designate SNA-3. SNA-3 is an essential, nuclear protein with three NADAR domains whose function is unknown. Mutation of key residues in NADAR domains inactivates the protein, indicating that domain function is required for activity. SNA-3 interacts with the CBC-ARS2 complex and other factors involved in RNA metabolism, including SUT-1 protein, through RNA or protein-mediated contacts revealed by yeast two-hybrid assays, localisation studies and immunoprecipitations. Our data are compatible with a role for SNA-3 in coordinating trans-splicing with target pre-mRNA transcription or in the processing of the Y-branch product of the trans-splicing reaction.

Genome annotation of Caenorhabditis briggsae by TEC-RED identifies new exons, paralogs, and conserved and novel operons

The nematode Caenorhabditis briggsae is routinely used in comparative and evolutionary studies involving its well-known cousin Caenorhabditis elegans. The C. briggsae genome sequence has accelerated research by facilitating the generation of new resources, tools, and functional studies of genes. While substantial progress has been made in predicting genes and start sites, experimental evidence is still lacking in many cases. Here, we report an improved annotation of the C. briggsae genome using the trans-spliced exon coupled RNA end determination technique. In addition to identifying the 5' ends of expressed genes, we have discovered operons and paralogs. In summary, our analysis yielded 10,243 unique 5' end sequence tags with matches in the C. briggsae genome. Of these, 6,395 were found to represent 4,252 unique genes along with 362 paralogs and 52 previously unknown exons. These genes included 14 that are exclusively trans-spliced in C. briggsae when compared with C. elegans orthologs. A major contribution of this study is the identification of 492 high confidence operons, of which two-thirds are fully supported by tags. In addition, 2 SL1-type operons were discovered. Interestingly, comparisons with C. elegans showed that only 40% of operons are conserved. Of the remaining operons, 73 are novel, including 12 that entirely lack orthologs in C. elegans. Further analysis revealed that 4 of the 12 novel operons are conserved in Caenorhabditis nigoni. Altogether, the work described here has significantly advanced our understanding of the C. briggsae system and serves as a rich resource to aid biological studies involving this species.

Floxed exon (Flexon): A flexibly positioned stop cassette for recombinase-mediated conditional gene expression

Tools that afford spatiotemporal control of gene expression are crucial for studying genes and processes in multicellular organisms. Stop cassettes consist of exogenous sequences that interrupt gene expression and flanking site-specific recombinase sites to allow for tissue-specific excision and restoration of function by expression of the cognate recombinase. We describe a stop cassette called a flexon, composed of an artificial exon flanked by artificial introns that can be flexibly positioned in a gene. We demonstrate its efficacy in Caenorhabditis elegans for lineage-specific control of gene expression and for tissue-specific RNA interference and discuss other potential uses. The Flexon approach should be feasible in any system amenable to site-specific recombination-based methods and applicable to diverse areas including development, neuroscience, and metabolism.

Conditional gene expression is a powerful tool for genetic analysis of biological phenomena. In the widely used “lox-stop-lox” approach, insertion of a stop cassette consisting of a series of stop codons and polyadenylation signals flanked by lox sites into the 5′ untranslated region (UTR) of a gene prevents expression until the cassette is excised by tissue-specific expression of Cre recombinase. Although lox-stop-lox and similar approaches using other site-specific recombinases have been successfully used in many experimental systems, this design has certain limitations. Here, we describe the Floxed exon (Flexon) approach, which uses a stop cassette composed of an artificial exon flanked by artificial introns, designed to cause premature termination of translation and nonsense-mediated decay of the mRNA and allowing for flexible placement into a gene. We demonstrate its efficacy in Caenorhabditis elegans by showing that, when promoters that cause weak and/or transient cell-specific expression are used to drive Cre in combination with a gfp(flexon) transgene, strong and sustained expression of green fluorescent protein (GFP) is obtained in specific lineages. We also demonstrate its efficacy in an endogenous gene context: we inserted a flexon into the Argonaute gene rde-1 to abrogate RNA interference (RNAi), and restored RNAi tissue specifically by expression of Cre. Finally, we describe several potential additional applications of the Flexon approach, including more precise control of gene expression using intersectional methods, tissue-specific protein degradation, and generation of genetic mosaics. The Flexon approach should be feasible in any system where a site-specific recombination-based method may be applied.

Systematic transcriptome analysis associated with physiological and chronological aging in Caenorhabditis elegans

Aging is associated with changes in a variety of biological processes at the transcriptomic level, including gene expression. Two types of aging occur during a lifetime: chronological and physiological aging. However, dissecting the difference between chronological and physiological ages at the transcriptomic level has been a challenge because of its complexity. We analyzed the transcriptomic features associated with physiological and chronological aging using Caenorhabditis elegans as a model. Many structural and functional transcript elements, such as noncoding RNAs and intron-derived transcripts, were up-regulated with chronological aging. In contrast, mRNAs with many biological functions, including RNA processing, were down-regulated with physiological aging. We also identified an age-dependent increase in the usage of distal 3' splice sites in mRNA transcripts as a biomarker of physiological aging. Our study provides crucial information for dissecting chronological and physiological aging at the transcriptomic level.

Gene expression profiling of epidermal cell types in C. elegans using Targeted DamID

The epidermis of Caenorhabditis elegans is an essential tissue for survival because it contributes to the formation of the cuticle barrier as well as facilitating developmental progression and animal growth. Most of the epidermis consists of the hyp7 hypodermal syncytium, the nuclei of which are largely generated by the seam cells, which exhibit stem cell-like behaviour during development. How seam cell progenitors differ transcriptionally from the differentiated hypodermis is poorly understood. Here, we introduce Targeted DamID (TaDa) in C. elegans as a method for identifying genes expressed within a tissue of interest without cell isolation. We show that TaDa signal enrichment profiles can be used to identify genes transcribed in the epidermis and use this method to resolve differences in gene expression between the seam cells and the hypodermis. Finally, we predict and functionally validate new transcription and chromatin factors acting in seam cell development. These findings provide insights into cell type-specific gene expression profiles likely associated with epidermal cell fate patterning.

Nanopore sequencing reveals endogenous NMD-targeted isoforms in human cells

BACKGROUND

Nonsense-mediated mRNA decay (NMD) is a eukaryotic, translation-dependent degradation pathway that targets mRNAs with premature termination codons and also regulates the expression of some mRNAs that encode full-length proteins. Although many genes express NMD-sensitive transcripts, identifying them based on short-read sequencing data remains a challenge.

RESULTS

To identify and analyze endogenous targets of NMD, we apply cDNA Nanopore sequencing and short-read sequencing to human cells with varying expression levels of NMD factors. Our approach detects full-length NMD substrates that are highly unstable and increase in levels or even only appear when NMD is inhibited. Among the many new NMD-targeted isoforms that our analysis identifies, most derive from alternative exon usage. The isoform-aware analysis reveals many genes with significant changes in splicing but no significant changes in overall expression levels upon NMD knockdown. NMD-sensitive mRNAs have more exons in the 3΄UTR and, for those mRNAs with a termination codon in the last exon, the length of the 3΄UTR per se does not correlate with NMD sensitivity. Analysis of splicing signals reveals isoforms where NMD has been co-opted in the regulation of gene expression, though the main function of NMD seems to be ridding the transcriptome of isoforms resulting from spurious splicing events.

CONCLUSIONS

Long-read sequencing enables the identification of many novel NMD-sensitive mRNAs and reveals both known and unexpected features concerning their biogenesis and their biological role. Our data provide a highly valuable resource of human NMD transcript targets for future genomic and transcriptomic applications.

Molecular topography of an entire nervous system

We have produced gene expression profiles of all 302 neurons of the C. elegans nervous system that match the single-cell resolution of its anatomy and wiring diagram. Our results suggest that individual neuron classes can be solely identified by combinatorial expression of specific gene families. For example, each neuron class expresses distinct codes of ∼23 neuropeptide genes and ∼36 neuropeptide receptors, delineating a complex and expansive "wireless" signaling network. To demonstrate the utility of this comprehensive gene expression catalog, we used computational approaches to (1) identify cis-regulatory elements for neuron-specific gene expression and (2) reveal adhesion proteins with potential roles in process placement and synaptic specificity. Our expression data are available at https://cengen.org and can be interrogated at the web application CengenApp. We expect that this neuron-specific directory of gene expression will spur investigations of underlying mechanisms that define anatomy, connectivity, and function throughout the C. elegans nervous system.

m6 A-mediated alternative splicing coupled with nonsense-mediated mRNA decay regulates SAM synthetase homeostasis

Alternative splicing of pre-mRNAs can regulate gene expression levels by coupling with nonsense-mediated mRNA decay (NMD). In order to elucidate a repertoire of mRNAs regulated by alternative splicing coupled with NMD (AS-NMD) in an organism, we performed long-read RNA sequencing of poly(A)⁺ RNAs from an NMD-deficient mutant strain of Caenorhabditis elegans, and obtained full-length sequences for mRNA isoforms from 259 high-confidence AS-NMD genes. Among them are the S-adenosyl-L-methionine (SAM) synthetase (sams) genes sams-3 and sams-4. SAM synthetase activity autoregulates sams gene expression through AS-NMD in a negative feedback loop. We furthermore find that METT-10, the orthologue of human U6 snRNA methyltransferase METTL16, is required for the splicing regulation in␣vivo, and specifically methylates the invariant AG dinucleotide at the distal 3' splice site (3'SS) in␣vitro. Direct RNA sequencing coupled with machine learning confirms m⁶ A modification of endogenous sams mRNAs. Overall, these results indicate that homeostasis of SAM synthetase in C. elegans is maintained by alternative splicing regulation through m⁶ A modification at the 3'SS of the sams genes.

SLIDR and SLOPPR: flexible identification of spliced leader trans-splicing and prediction of eukaryotic operons from RNA-Seq data

BACKGROUND

Spliced leader (SL) trans-splicing replaces the 5' end of pre-mRNAs with the spliced leader, an exon derived from a specialised non-coding RNA originating from elsewhere in the genome. This process is essential for resolving polycistronic pre-mRNAs produced by eukaryotic operons into monocistronic transcripts. SL trans-splicing and operons may have independently evolved multiple times throughout Eukarya, yet our understanding of these phenomena is limited to only a few well-characterised organisms, most notably C. elegans and trypanosomes. The primary barrier to systematic discovery and characterisation of SL trans-splicing and operons is the lack of computational tools for exploiting the surge of transcriptomic and genomic resources for a wide range of eukaryotes.

RESULTS

Here we present two novel pipelines that automate the discovery of SLs and the prediction of operons in eukaryotic genomes from RNA-Seq data. SLIDR assembles putative SLs from 5' read tails present after read alignment to a reference genome or transcriptome, which are then verified by interrogating corresponding SL RNA genes for sequence motifs expected in bona fide SL RNA molecules. SLOPPR identifies RNA-Seq reads that contain a given 5' SL sequence, quantifies genome-wide SL trans-splicing events and predicts operons via distinct patterns of SL trans-splicing events across adjacent genes. We tested both pipelines with organisms known to carry out SL trans-splicing and organise their genes into operons, and demonstrate that (1) SLIDR correctly detects expected SLs and often discovers novel SL variants; (2) SLOPPR correctly identifies functionally specialised SLs, correctly predicts known operons and detects plausible novel operons.

CONCLUSIONS

SLIDR and SLOPPR are flexible tools that will accelerate research into the evolutionary dynamics of SL trans-splicing and operons throughout Eukarya and improve gene discovery and annotation for a wide range of eukaryotic genomes. Both pipelines are implemented in Bash and R and are built upon readily available software commonly installed on most bioinformatics servers. Biological insight can be gleaned even from sparse, low-coverage datasets, implying that an untapped wealth of information can be retrieved from existing RNA-Seq datasets as well as from novel full-isoform sequencing protocols as they become more widely available.

Resolution of polycistronic RNA by SL2 trans-splicing is a widely conserved nematode trait

Spliced leader trans-splicing is essential for the processing and translation of polycistronic RNAs generated by eukaryotic operons. In C. elegans, a specialized spliced leader, SL2, provides the 5' end for uncapped pre-mRNAs derived from polycistronic RNAs. Studies of other nematodes suggested that SL2-type trans-splicing is a relatively recent innovation, confined to Rhabditina, the clade containing C. elegans and its close relatives. Here we conduct a survey of transcriptome-wide spliced leader trans-splicing in Trichinella spiralis, a distant relative of C. elegans with a particularly diverse repertoire of 15 spliced leaders. By systematically comparing the genomic context of trans-splicing events for each spliced leader, we identified a subset of T. spiralis spliced leaders that are specifically used to process polycistronic RNAs-the first examples of SL2-type spliced leaders outside of Rhabditina. These T. spiralis spliced leader RNAs possess a perfectly conserved stem-loop motif previously shown to be essential for SL2-type trans-splicing in C. elegans We show that genes trans-spliced to these SL2-type spliced leaders are organized in operonic fashion, with short intercistronic distances. A subset of T. spiralis operons show conservation of synteny with C. elegans operons. Our work substantially revises our understanding of nematode spliced leader trans-splicing, showing that SL2 trans-splicing is a major mechanism for nematode polycistronic RNA processing, which may have evolved prior to the radiation of the Nematoda. This work has important implications for the improvement of genome annotation pipelines in nematodes and other eukaryotes with operonic gene organization.

Global regulatory features of alternative splicing across tissues and within the nervous system of C. elegans

Alternative splicing plays a major role in shaping tissue-specific transcriptomes. Among the broad tissue types present in metazoans, the central nervous system contains some of the highest levels of alternative splicing. Although many documented examples of splicing differences between broad tissue types exist, there remains much to be understood about the splicing factors and the cis sequence elements controlling tissue and neuron subtype-specific splicing patterns. By using translating ribosome affinity purification coupled with deep-sequencing (TRAP-seq) in Caenorhabditis elegans, we have obtained high coverage profiles of ribosome-associated mRNA for three broad tissue classes (nervous system, muscle, and intestine) and two neuronal subtypes (dopaminergic and serotonergic neurons). We have identified hundreds of splice junctions that exhibit distinct splicing patterns between tissue types or within the nervous system. Alternative splicing events differentially regulated between tissues are more often frame-preserving, are more highly conserved across Caenorhabditis species, and are enriched in specific cis regulatory motifs, when compared with other types of exons. By using this information, we have identified a likely mechanism of splicing repression by the RNA-binding protein UNC-75/CELF via interactions with cis elements that overlap a 5′ splice site. Alternatively spliced exons also overlap more frequently with intrinsically disordered peptide regions than constitutive exons. Moreover, regulated exons are often shorter than constitutive exons but are flanked by longer intron sequences. Among these tissue-regulated exons are several highly conserved microexons <27 nt in length. Collectively, our results indicate a rich layer of tissue-specific gene regulation at the level of alternative splicing in C. elegans that parallels the evolutionary forces and constraints observed across metazoa.

mRNA Editing, Processing and Quality Control in Caenorhabditis elegans

While DNA serves as the blueprint of life, the distinct functions of each cell are determined by the dynamic expression of genes from the static genome. The amount and specific sequences of RNAs expressed in a given cell involves a number of regulated processes including RNA synthesis (transcription), processing, splicing, modification, polyadenylation, stability, translation, and degradation. As errors during mRNA production can create gene products that are deleterious to the organism, quality control mechanisms exist to survey and remove errors in mRNA expression and processing. Here, we will provide an overview of mRNA processing and quality control mechanisms that occur in Caenorhabditis elegans, with a focus on those that occur on protein-coding genes after transcription initiation. In addition, we will describe the genetic and technical approaches that have allowed studies in C. elegans to reveal important mechanistic insight into these processes.

Direct full-length RNA sequencing reveals unexpected transcriptome complexity during Caenorhabditis elegans development

Massively parallel sequencing of the polyadenylated RNAs has played a key role in delineating transcriptome complexity, including alternative use of an exon, promoter, 5′ or 3′ splice site or polyadenylation site, and RNA modification. However, reads derived from the current RNA-seq technologies are usually short and deprived of information on modification, compromising their potential in defining transcriptome complexity. Here, we applied a direct RNA sequencing method with ultralong reads using Oxford Nanopore Technologies to study the transcriptome complexity in Caenorhabditis elegans. We generated approximately six million reads using native poly(A)-tailed mRNAs from three developmental stages, with average read lengths ranging from 900 to 1100 nt. Around half of the reads represent full-length transcripts. To utilize the full-length transcripts in defining transcriptome complexity, we devised a method to classify the long reads as the same as existing transcripts or as a novel transcript using sequence mapping tracks rather than existing intron/exon structures, which allowed us to identify roughly 57,000 novel isoforms and recover at least 26,000 out of the 33,500 existing isoforms. The sets of genes with differential expression versus differential isoform usage over development are largely different, implying a fine-tuned regulation at isoform level. We also observed an unexpected increase in putative RNA modification in all bases in the coding region relative to the UTR, suggesting their possible roles in translation. The RNA reads and the method for read classification are expected to deliver new insights into RNA processing and modification and their underlying biology in the future.

The full-length transcriptome of C. elegans using direct RNA sequencing

Current transcriptome annotations have largely relied on short read lengths intrinsic to the most widely used high-throughput cDNA sequencing technologies. For example, in the annotation of the Caenorhabditis elegans transcriptome, more than half of the transcript isoforms lack full-length support and instead rely on inference from short reads that do not span the full length of the isoform. We applied nanopore-based direct RNA sequencing to characterize the developmental polyadenylated transcriptome of C. elegans Taking advantage of long reads spanning the full length of mRNA transcripts, we provide support for 23,865 splice isoforms across 14,611 genes, without the need for computational reconstruction of gene models. Of the isoforms identified, 3452 are novel splice isoforms not present in the WormBase WS265 annotation. Furthermore, we identified 16,342 isoforms in the 3' untranslated region (3' UTR), 2640 of which are novel and do not fall within 10 bp of existing 3'-UTR data sets and annotations. Combining 3' UTRs and splice isoforms, we identified 28,858 full-length transcript isoforms. We also determined that poly(A) tail lengths of transcripts vary across development, as do the strengths of previously reported correlations between poly(A) tail length and expression level, and poly(A) tail length and 3'-UTR length. Finally, we have formatted this data as a publicly accessible track hub, enabling researchers to explore this data set easily in a genome browser.

Transcriptional adaptation in Caenorhabditis elegans

Transcriptional adaptation is a recently described phenomenon by which a mutation in one gene leads to the transcriptional modulation of related genes, termed adapting genes. At the molecular level, it has been proposed that the mutant mRNA, rather than the loss of protein function, activates this response. While several examples of transcriptional adaptation have been reported in zebrafish embryos and in mouse cell lines, it is not known whether this phenomenon is observed across metazoans. Here we report transcriptional adaptation in C. elegans, and find that this process requires factors involved in mutant mRNA decay, as in zebrafish and mouse. We further uncover a requirement for Argonaute proteins and Dicer, factors involved in small RNA maturation and transport into the nucleus. Altogether, these results provide evidence for transcriptional adaptation in C. elegans, a powerful model to further investigate underlying molecular mechanisms.

The C. elegans 3' UTRome v2 resource for studying mRNA cleavage and polyadenylation, 3'-UTR biology, and miRNA targeting

3′ Untranslated regions (3′ UTRs) of mRNAs emerged as central regulators of cellular function because they contain important but poorly characterized cis-regulatory elements targeted by a multitude of regulatory factors. The model nematode Caenorhabditis elegans is ideal to study these interactions because it possesses a well-defined 3′ UTRome. To improve its annotation, we have used a genome-wide bioinformatics approach to download raw transcriptome data for 1088 transcriptome data sets corresponding to the entire collection of C. elegans trancriptomes from 2015 to 2018 from the Sequence Read Archive at the NCBI. We then extracted and mapped high-quality 3′-UTR data at ultradeep coverage. Here, we describe and release to the community the updated version of the worm 3′ UTRome, which we named 3′ UTRome v2. This resource contains high-quality 3′-UTR data mapped at single-base ultraresolution for 23,084 3′-UTR isoform variants corresponding to 14,788 protein-coding genes and is updated to the latest release of WormBase. We used this data set to study and probe principles of mRNA cleavage and polyadenylation in C. elegans. The worm 3′ UTRome v2 represents the most comprehensive and high-resolution 3′-UTR data set available in C. elegans and provides a novel resource to investigate the mRNA cleavage and polyadenylation reaction, 3′-UTR biology, and miRNA targeting in a living organism.

Genome-wide alternative polyadenylation dynamics in response to biotic and abiotic stresses in rice

Alternative polyadenylation (APA) is an important way to regulate gene expression at the post-transcriptional level, and is extensively involved in plant stress responses. However, the systematic roles of APA regulation in response to abiotic and biotic stresses in rice at the genome scale remain unknown. To take advantage of available RNA-seq datasets, using a novel tool APAtrap, we identified thousands of genes with significantly differential usage of polyadenylation [poly(A)] sites in response to the abiotic stress (drought, heat shock, and cadmium) and biotic stress [bacterial blight (BB), rice blast, and rice stripe virus (RSV)]. Genes with stress-responsive APA dynamics commonly exhibited higher expression levels when their isoforms with short 3' untranslated region (3' UTR) were more abundant. The stress-responsive APA events were widely involved in crucial stress-responsive genes and pathways: e.g. APA acted as a negative regulator in heat stress tolerance; APA events were involved in DNA repair and cell wall formation under Cd stress; APA regulated chlorophyll metabolism, being associated with the pathogenesis of leaf diseases under RSV and BB challenges. Furthermore, APA events were found to be involved in glutathione metabolism and MAPK signaling pathways, mediating a crosstalk among the abiotic and biotic stress-responsive regulatory networks in rice. Analysis of large-scale datasets revealed that APA may regulate abiotic and biotic stress-responsive processes in rice. Such post-transcriptome diversities contribute to rice adaption to various environmental challenges. Our study would supply useful resource for further molecular assisted breeding of multiple stress-tolerant cultivars for rice.

CRISPR editing of sftb-1/SF3B1 in Caenorhabditis elegans allows the identification of synthetic interactions with cancer-related mutations and the chemical inhibition of splicing

SF3B1 is the most frequently mutated splicing factor in cancer. Mutations in SF3B1 likely confer clonal advantages to cancer cells but they may also confer vulnerabilities that can be therapeutically targeted. SF3B1 cancer mutations can be maintained in homozygosis in C. elegans, allowing synthetic lethal screens with a homogeneous population of animals. These mutations cause alternative splicing (AS) defects in C. elegans, as it occurs in SF3B1-mutated human cells. In a screen, we identified RNAi of U2 snRNP components that cause synthetic lethality with sftb-1/SF3B1 mutations. We also detected synthetic interactions between sftb-1 mutants and cancer-related mutations in uaf-2/U2AF1 or rsp-4/SRSF2, demonstrating that this model can identify interactions between mutations that are mutually exclusive in human tumors. Finally, we have edited an SFTB-1 domain to sensitize C. elegans to the splicing modulators pladienolide B and herboxidiene. Thus, we have established a multicellular model for SF3B1 mutations amenable for high-throughput genetic and chemical screens.

The C. elegans embryonic transcriptome with tissue, time, and alternative splicing resolution

We have used RNA-seq in Caenorhabditis elegans to produce transcription profiles for seven specific embryonic cell populations from gastrulation to the onset of terminal differentiation. The expression data for these seven cell populations, covering major cell lineages and tissues in the worm, reveal the complex and dynamic changes in gene expression, both spatially and temporally. Also, within genes, start sites and exon usage can be highly differential, producing transcripts that are specific to developmental periods or cell lineages. We have also found evidence of novel exons and introns, as well as differential usage of SL1 and SL2 splice leaders. By combining this data set with the modERN ChIP-seq resource, we are able to support and predict gene regulatory relationships. The detailed information on differences and similarities between gene expression in cell lineages and tissues should be of great value to the community and provides a framework for the investigation of expression in individual cells.

Heterodera glycines utilizes promiscuous spliced leaders and demonstrates a unique preference for a species-specific spliced leader over C. elegans SL1

Spliced leader trans-splicing (SLTS) plays a part in the maturation of pre-mRNAs in select species across multiple phyla but is particularly prevalent in Nematoda. The role of spliced leaders (SL) within the cell is unclear and an accurate assessment of SL occurrence within an organism is possible only after extensive sequencing data are available, which is not currently the case for many nematode species. SL discovery is further complicated by an absence of SL sequences from high-throughput sequencing results due to incomplete sequencing of the 5'-ends of transcripts during RNA-seq library preparation, known as 5'-bias. Existing datasets and novel methodology were used to identify both conserved SLs and unique hypervariable SLs within Heterodera glycines, the soybean cyst nematode. In H. glycines, twenty-one distinct SL sequences were found on 2,532 unique H. glycines transcripts. The SL sequences identified on the H. glycines transcripts demonstrated a high level of promiscuity, meaning that some transcripts produced as many as nine different individual SL-transcript combinations. Most uniquely, transcriptome analysis revealed that H. glycines is the first nematode to demonstrate a higher SL trans-splicing rate using a species-specific SL over well-conserved Caenorhabditis elegans SL-like sequences.

Modeling and analysis of RNA-seq data: a review from a statistical perspective

BACKGROUND

Since the invention of next-generation RNA sequencing (RNA-seq) technologies, they have become a powerful tool to study the presence and quantity of RNA molecules in biological samples and have revolutionized transcriptomic studies. The analysis of RNA-seq data at four different levels (samples, genes, transcripts, and exons) involve multiple statistical and computational questions, some of which remain challenging up to date.

RESULTS

We review RNA-seq analysis tools at the sample, gene, transcript, and exon levels from a statistical perspective. We also highlight the biological and statistical questions of most practical considerations.

CONCLUSIONS

The development of statistical and computational methods for analyzing RNA-seq data has made significant advances in the past decade. However, methods developed to answer the same biological question often rely on diverse statistical models and exhibit different performance under different scenarios. This review discusses and compares multiple commonly used statistical models regarding their assumptions, in the hope of helping users select appropriate methods as needed, as well as assisting developers for future method development.

SL-quant: a fast and flexible pipeline to quantify spliced leader trans-splicing events from RNA-seq data

Background

The spliceosomal transfer of a short spliced leader (SL) RNA to an independent pre-mRNA molecule is called SL trans-splicing and is widespread in the nematode Caenorhabditis elegans. While RNA-sequencing (RNA-seq) data contain information on such events, properly documented methods to extract them are lacking.

Findings

To address this, we developed SL-quant, a fast and flexible pipeline that adapts to paired-end and single-end RNA-seq data and accurately quantifies SL trans-splicing events. It is designed to work downstream of read mapping and uses the reads left unmapped as primary input. Briefly, the SL sequences are identified with high specificity and are trimmed from the input reads, which are then remapped on the reference genome and quantified at the nucleotide position level (SL trans-splice sites) or at the gene level.

Conclusions

SL-quant completes within 10 minutes on a basic desktop computer for typical C. elegans RNA-seq datasets and can be applied to other species as well. Validating the method, the SL trans-splice sites identified display the expected consensus sequence, and the results of the gene-level quantification are predictive of the gene position within operons. We also compared SL-quant to a recently published SL-containing read identification strategy that was found to be more sensitive but less specific than SL-quant. Both methods are implemented as a bash script available under the MIT license [1]. Full instructions for its installation, usage, and adaptation to other organisms are provided.