Trans-splicing

Overcoming cancer drug-resistance calls for novel strategies targeting abnormal alternative splicing

Abnormal gene alternative splicing (AS) events are strongly associated with cancer progression. Here, we summarize AS events that contribute to the development of drug resistance and classify them into three categories: alternative cis-splicing (ACS), alternative trans-splicing (ATS), and alternative back-splicing (ABS). The regulatory mechanisms underlying AS processes through cis-acting regulatory elements and trans-acting factors are comprehensively described, and the distinct functions of spliced variants, including linear spliced variants derived from ACS, chimeric spliced variants arising from ATS, and circRNAs generated through ABS, are discussed. The identification of dysregulated spliced variants, which contribute to drug resistance and hinder effective cancer treatment, suggests that abnormal AS processes may together serve as a precise regulatory mechanism enabling drug-resistant cancer cell survival or, alternatively, represent an evolutionary pathway for cancer cells to adapt to changes in the external environment. Moreover, this review summarizes recent advancements in treatment approaches targeting AS-associated drug resistance, focusing on cis-acting regulatory elements, trans-acting factors, and specific spliced variants. Collectively, gaining an in-depth understanding of the mechanisms underlying aberrant alternative splicing events and developing strategies to target this process hold great promise for overcoming cancer drug resistance.

U6 snRNA m6A modification is required for accurate and efficient splicing of C. elegans and human pre-mRNAs

pre-mRNA splicing is a critical feature of eukaryotic gene expression. Both cis- and trans-splicing rely on accurately recognising splice site sequences by spliceosomal U snRNAs and associated proteins. Spliceosomal snRNAs carry multiple RNA modifications with the potential to affect different stages of pre-mRNA splicing. Here, we show that the conserved U6 snRNA m6A methyltransferase METT-10 is required for accurate and efficient cis- and trans-splicing of C. elegans pre-mRNAs. The absence of METT-10 in C. elegans and METTL16 in humans primarily leads to alternative splicing at 5' splice sites with an adenosine at +4 position. In addition, METT-10 is required for splicing of weak 3' cis- and trans-splice sites. We identified a significant overlap between METT-10 and the conserved splicing factor SNRNP27K in regulating 5' splice sites with +4A. Finally, we show that editing endogenous 5' splice site +4A positions to +4U restores splicing to wild-type positions in a mett-10 mutant background, supporting a direct role for U6 snRNA m6A modification in 5' splice site recognition. We conclude that the U6 snRNA m6A modification is important for accurate and efficient pre-mRNA splicing.

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.

Unique features of conventional and nonconventional introns in Euglena gracilis

BACKGROUND

Nuclear introns in Euglenida have been understudied. This study aimed to investigate nuclear introns in Euglenida by identifying a large number of introns in Euglena gracilis (E. gracilis), including cis-spliced conventional and nonconventional introns, as well as trans-spliced outrons. We also examined the sequence characteristics of these introns.

RESULTS

A total of 28,337 introns and 11,921 outrons were identified. Conventional and nonconventional introns have distinct splice site features; the former harbour canonical GT/C-AG splice sites, whereas the latter are capable of forming structured motifs with their terminal sequences. We observed that short introns had a preference for canonical GT-AG introns. Notably, conventional introns and outrons in E. gracilis exhibited a distinct cytidine-rich polypyrimidine tract, in contrast to the thymidine-rich tracts observed in other organisms. Furthermore, the SL-RNAs in E. gracilis, as well as in other trans-splicing species, can form a recently discovered motif called the extended U6/5' ss duplex with the respective U6s. We also describe a novel type of alternative splicing pattern in E. gracilis. The tandem repeat sequences of introns in this protist were determined, and their contents were comparable to those in humans.

CONCLUSIONS

Our findings highlight the unique features of E. gracilis introns and provide insights into the splicing mechanism of these introns, as well as the genomics and evolution of Euglenida.

Co-transcriptional splicing facilitates transcription of gigantic genes

Although introns are typically tens to thousands of nucleotides, there are notable exceptions. In flies as well as humans, a small number of genes contain introns that are more than 1000 times larger than typical introns, exceeding hundreds of kilobases (kb) to megabases (Mb). It remains unknown why gigantic introns exist and how cells overcome the challenges associated with their transcription and RNA processing. The Drosophila Y chromosome contains some of the largest genes identified to date: multiple genes exceed 4Mb, with introns accounting for over 99% of the gene span. Here we demonstrate that co-transcriptional splicing of these gigantic Y-linked genes is important to ensure successful transcription: perturbation of splicing led to the attenuation of transcription, leading to a failure to produce mature mRNA. Cytologically, defective splicing of the Y-linked gigantic genes resulted in disorganization of transcripts within the nucleus suggestive of entanglement of transcripts, likely resulting from unspliced long RNAs. We propose that co-transcriptional splicing maintains the length of nascent transcripts of gigantic genes under a critical threshold, preventing their entanglement and ensuring proper gene expression. Our study reveals a novel biological significance of co-transcriptional splicing.

Characterization of trans-spliced chimeric RNAs: insights into the mechanism of trans-splicing

Trans-splicing is a post-transcriptional processing event that joins exons from separate RNAs to produce a chimeric RNA. However, the detailed mechanism of trans-splicing remains poorly understood. Here, we characterize trans-spliced genes and provide insights into the mechanism of trans-splicing in the tunicate Ciona. Tunicates are the closest invertebrates to humans, and their genes frequently undergo trans-splicing. Our analysis revealed that, in genes that give rise to both trans-spliced and non-trans-spliced messenger RNAs, trans-splice acceptor sites were preferentially located at the first functional acceptor site, and their paired donor sites were weak in both Ciona and humans. Additionally, we found that Ciona trans-spliced genes had GU- and AU-rich 5' transcribed regions. Our data and findings not only are useful for Ciona research community, but may also aid in a better understanding of the trans-splicing mechanism, potentially advancing the development of gene therapy based on trans-splicing.

On a kneading theory for gene-splicing

Two well-known facets in protein synthesis in eukaryotic cells are transcription of DNA to pre-RNA in the nucleus and the translation of messenger-RNA (mRNA) to proteins in the cytoplasm. A critical intermediate step is the removal of segments (introns) containing ∼97% of the nucleic-acid sites in pre-RNA and sequential alignment of the retained segments (exons) to form mRNA through a process referred to as splicing. Alternative forms of splicing enrich the proteome while abnormal splicing can enhance the likelihood of a cell developing cancer or other diseases. Mechanisms for splicing and origins of splicing errors are only partially deciphered. Our goal is to determine if rules on splicing can be inferred from data analytics on nucleic-acid sequences. Toward that end, we represent a nucleic-acid site as a point in a plane defined in terms of the anterior and posterior sub-sequences of the site. The "point-set" representation expands analytical approaches, including the use of statistical tools, to characterize genome sequences. It is found that point-sets for exons and introns are visually different, and that the differences can be quantified using a family of generalized moments. We design a machine-learning algorithm that can recognize individual exons or introns with 91% accuracy. Point-set distributions and generalized moments are found to differ between organisms.

Nano biomaterials based strategies for enhanced brain targeting in the treatment of neurodegenerative diseases: an up-to-date perspective

Neurons and their connecting axons gradually degenerate in neurodegenerative diseases (NDs), leading to dysfunctionality of the neuronal cells and eventually their death. Drug delivery for the treatment of effected nervous system is notoriously complicated because of the presence of natural barriers, i.e., the blood-brain barrier and the blood cerebrospinal fluid barrier. Palliative care is currently the standard care for many diseases. Therefore, treatment programs that target the disease's origin rather than its symptoms are recommended. Nanotechnology-based drug delivery platforms offer an innovative way to circumvent these obstacles and deliver medications directly to the central nervous system, thereby enabling treatment of several common neurological problems, i.e., Alzheimer's, Parkinson's, Huntington's, and amyotrophic lateral sclerosis. Interestingly, the combination of nanomedicine and gene therapy enables targeting of selective mutant genes responsible for the progression of NDs, which may provide a much-needed boost in the struggle against these diseases. Herein, we discussed various central nervous system delivery obstacles, followed by a detailed insight into the recently developed techniques to restore neurological function via the differentiation of neural stem cells. Moreover, a comprehensive background on the role of nanomedicine in controlling neurogenesis via differentiation of neural stem cells is explained. Additionally, numerous phytoconstituents with their neuroprotective properties and molecular targets in the identification and management of NDs are also deliberated. Furthermore, a detailed insight of the ongoing clinical trials and currently marketed products for the treatment of NDs is provided in this manuscript.

A novel chimeric RNA originating from BmCPV S4 and Bombyx mori HDAC11 transcripts regulates virus proliferation

Polymerases encoded by segmented negative-strand RNA viruses cleave 5'-m7G-capped host transcripts to prime viral mRNA synthesis ("cap-snatching") to generate chimeric RNA, and trans-splicing occurs between viral and cellular transcripts. Bombyx mori cytoplasmic polyhedrosis virus (BmCPV), an RNA virus belonging to Reoviridae, is a major pathogen of silkworm (B. mori). The genome of BmCPV consists of 10 segmented double-stranded RNAs (S1-S10) from which viral RNAs encoding a protein are transcribed. In this study, chimeric silkworm-BmCPV RNAs, in which the sequence derived from the silkworm transcript could fuse with both the 5' end and the 3' end of viral RNA, were identified in the midgut of BmCPV-infected silkworms by RNA_seq and further confirmed by RT-PCR and Sanger sequencing. A novel chimeric RNA, HDAC11-S4 RNA 4, derived from silkworm histone deacetylase 11 (HDAC11) and the BmCPV S4 transcript encoding viral structural protein 4 (VP4), was selected for validation by in situ hybridization and Northern blotting. Interestingly, our results indicated that HDAC11-S4 RNA 4 was generated in a BmCPV RNA-dependent RNA polymerase (RdRp)-independent manner and could be translated into a truncated BmCPV VP4 with a silkworm HDAC11-derived N-terminal extension. Moreover, it was confirmed that HDAC11-S4 RNA 4 inhibited BmCPV proliferation, decreased the level of H3K9me3 and increased the level of H3K9ac. These results indicated that during infection with BmCPV, a novel mechanism, different from that described in previous reports, allows the genesis of chimeric silkworm-BmCPV RNAs with biological functions.

U6 snRNA m6A modification is required for accurate and efficient cis- and trans-splicing of C. elegans mRNAs

pre-mRNA splicing is a critical feature of eukaryotic gene expression. Many eukaryotes use cis-splicing to remove intronic sequences from pre-mRNAs. In addition to cis-splicing, many organisms use trans-splicing to replace the 5' ends of mRNAs with a non-coding spliced-leader RNA. Both cis- and trans-splicing rely on accurately recognising splice site sequences by spliceosomal U snRNAs and associated proteins. Spliceosomal snRNAs carry multiple RNA modifications with the potential to affect different stages of pre-mRNA splicing. Here, we show that m6A modification of U6 snRNA A43 by the RNA methyltransferase METT-10 is required for accurate and efficient cis- and trans-splicing of C. elegans pre-mRNAs. The absence of U6 snRNA m6A modification primarily leads to alternative splicing at 5' splice sites. Furthermore, weaker 5' splice site recognition by the unmodified U6 snRNA A43 affects splicing at 3' splice sites. U6 snRNA m6A43 and the splicing factor SNRNP27K function to recognise an overlapping set of 5' splice sites with an adenosine at +4 position. Finally, we show that U6 snRNA m6A43 is required for efficient SL trans-splicing at weak 3' trans-splice sites. We conclude that the U6 snRNA m6A modification is important for accurate and efficient cis- and trans-splicing in C. elegans.

Genome assemblies of the simultaneously hermaphroditic flatworms Macrostomum cliftonense and Macrostomum hystrix

The free-living, simultaneously hermaphroditic flatworms of the genus Macrostomum are increasingly used as model systems in various contexts. In particular, Macrostomum lignano, the only species of this group with a published genome assembly, has emerged as a model for the study of regeneration, reproduction, and stem-cell function. However, challenges have emerged due to M. lignano being a hidden polyploid, having recently undergone whole-genome duplication and chromosome fusion events. This complex genome architecture presents a significant roadblock to the application of many modern genetic tools. Hence, additional genomic resources for this genus are needed. Here, we present such resources for Macrostomum cliftonense and Macrostomum hystrix, which represent the contrasting mating behaviors of reciprocal copulation and hypodermic insemination found in the genus. We use a combination of PacBio long-read sequencing and Illumina shot-gun sequencing, along with several RNA-Seq data sets, to assemble and annotate highly contiguous genomes for both species. The assemblies span ∼227 and ∼220 Mb and are represented by 399 and 42 contigs for M. cliftonense and M. hystrix, respectively. Furthermore, high BUSCO completeness (∼84-85%), low BUSCO duplication rates (8.3-6.2%), and low k-mer multiplicity indicate that these assemblies do not suffer from the same assembly ambiguities of the M. lignano genome assembly, which can be attributed to the complex karyology of this species. We also show that these resources, in combination with the prior resources from M. lignano, offer an excellent foundation for comparative genomic research in this group of organisms.

Roles of alternative splicing in infectious diseases: from hosts, pathogens to their interactions

ABSTRACT

Alternative splicing (AS) is an evolutionarily conserved mechanism that removes introns and ligates exons to generate mature messenger RNAs (mRNAs), extremely improving the richness of transcriptome and proteome. Both mammal hosts and pathogens require AS to maintain their life activities, and inherent physiological heterogeneity between mammals and pathogens makes them adopt different ways to perform AS. Mammals and fungi conduct a two-step transesterification reaction by spliceosomes to splice each individual mRNA (named cis -splicing). Parasites also use spliceosomes to splice, but this splicing can occur among different mRNAs (named trans -splicing). Bacteria and viruses directly hijack the host's splicing machinery to accomplish this process. Infection-related changes are reflected in the spliceosome behaviors and the characteristics of various splicing regulators (abundance, modification, distribution, movement speed, and conformation), which further radiate to alterations in the global splicing profiles. Genes with splicing changes are enriched in immune-, growth-, or metabolism-related pathways, highlighting approaches through which hosts crosstalk with pathogens. Based on these infection-specific regulators or AS events, several targeted agents have been developed to fight against pathogens. Here, we summarized recent findings in the field of infection-related splicing, including splicing mechanisms of pathogens and hosts, splicing regulation and aberrant AS events, as well as emerging targeted drugs. We aimed to systemically decode host-pathogen interactions from a perspective of splicing. We further discussed the current strategies of drug development, detection methods, analysis algorithms, and database construction, facilitating the annotation of infection-related splicing and the integration of AS with disease phenotype.

CPGView: A package for visualizing detailed chloroplast genome structures

Chloroplast genomes have been widely used in studying plant phylogeny and evolution. Several chloroplast genome visualization tools have been developed to display the distribution of genes on the genome. However, these tools do not draw features, such as exons, introns, repetitive elements, and variable sites, disallowing in-depth examination of the genome structures. Here, we developed and validated a software package called Chloroplast Genome Viewers (CPGView). CPGView can draw three maps showing (i) the distributions of genes, variable sites, and repetitive sequences, including microsatellites, tandem and dispersed repeats; (ii) the structure of the cis-splicing genes after adjusting the exon-intron boundary positions using a coordinate scaling algorithm, and (iii) the structure of the trans-splicing gene rps12. To test the accuracy of CPGView, we sequenced, assembled, and annotated 31 chloroplast genomes from 31 genera of 22 families. CPGView drew maps correctly for all the 31 chloroplast genomes. Lastly, we used CPGView to examine 5998 publicly released chloroplast genomes from 2513 genera of 553 families. CPGView succeeded in plotting maps for 5882 but failed to plot maps for 116 chloroplast genomes. Further examination showed that the annotations of these 116 genomes had various errors needing manual correction. The test on newly generated data and publicly available data demonstrated the ability of CPGView to identify errors in the annotations of chloroplast genomes. CPGView will become a widely used tool to study the detailed structure of chloroplast genomes. The web version of CPGView can be accessed from http://www.1kmpg.cn/cpgview.

Genome assembly of the acoel flatworm Symsagittifera roscoffensis, a model for research on body plan evolution and photosymbiosis

Symsagittifera roscoffensis is a well-known member of the order Acoela that lives in symbiosis with the algae Tetraselmis convolutae during its adult stage. Its natural habitat is the eastern coast of the Atlantic, where at specific locations thousands of individuals can be found, mostly, lying in large pools on the surface of sand at low tide. As a member of the Acoela it has been thought as a proxy for ancestral bilaterian animals; however, its phylogenetic position remains still debated. In order to understand the basic structural characteristics of the acoel genome, we sequenced and assembled the genome of aposymbiotic species S. roscoffensis. The size of this genome was measured to be in the range of 910-940 Mb. Sequencing of the genome was performed using PacBio Hi-Fi technology. Hi-C and RNA-seq data were also generated to scaffold and annotate it. The resulting assembly is 1.1 Gb large (covering 118% of the estimated genome size) and highly continuous, with N50 scaffold size of 1.04 Mb. The repetitive fraction of the genome is 61%, of which 85% (half of the genome) are LTR retrotransposons. Genome-guided transcriptome assembly identified 34,493 genes, of which 29,351 are protein coding (BUSCO score 97.6%), and 30.2% of genes are spliced leader trans-spliced. The completeness of this genome suggests that it can be used extensively to characterize gene families and conduct accurate phylogenomic reconstructions.

Modulation of Virulence Factors during Trypanosoma cruzi Differentiation

Chagas disease is a neglected tropical disease caused by Trypanosoma cruzi. This protozoan developed several mechanisms to infect, propagate, and survive in different hosts. The specific expression of proteins is responsible for morphological and metabolic changes in different parasite stages along the parasite life cycle. The virulence strategies at the cellular and molecular levels consist of molecules responsible for mediating resistance mechanisms to oxidative damage, cellular invasion, and immune evasion, performed mainly by surface proteins. Since parasite surface coat remodeling is crucial to invasion and infectivity, surface proteins are essential virulence elements. Understanding the factors involved in these processes improves the knowledge of parasite pathogenesis. Genome sequencing has opened the door to high-throughput technologies, allowing us to obtain a deeper understanding of gene reprogramming along the parasite life cycle and identify critical molecules for survival. This review therefore focuses on proteins regulated during differentiation into infective forms considered virulence factors and addresses the current known mechanisms acting in the modulation of gene expression, emphasizing mRNA signals, regulatory factors, and protein complexes.

Profile of chimeric RNAs and TMPRSS2-ERG e2e4 isoform in neuroendocrine prostate cancer

Purpose

Specific gene fusions and their fusion products (chimeric RNA and protein) have served as ideal diagnostic markers and therapeutic targets for cancer. However, few systematic studies for chimeric RNAs have been conducted in neuroendocrine prostate cancer (NEPC). In this study, we explored the landscape of chimeric RNAs in different types of prostate cancer (PCa) cell lines and aimed to identify chimeric RNAs specifically expressed in NEPC.

Methods

To do so, we employed the RNA-seq data of eight prostate related cell lines from Cancer Cell Line Encyclopedia (CCLE) for chimeric RNA identification. Multiple filtering criteria were used and the candidate chimeric RNAs were characterized at multiple levels and from various angles. We then performed experimental validation on all 80 candidates, and focused on the ones that are specific to NEPC. Lastly, we studied the clinical relevance and effect of one chimera in neuroendocrine process.

Results

Out of 80 candidates, 15 were confirmed to be expressed preferentially in NEPC lines. Among them, 13 of the 15 were found to be specifically expressed in NEPC, and four were further validated in another NEPC cell line. Importantly, in silico analysis showed that tumor malignancy may be correlated to the level of these chimeric RNAs. Clinically, the expression of TMPRSS2-ERG (e2e4) was elevated in tumor tissues and indicated poor clinical prognosis, whereas the parental wild type transcripts had no such association. Furthermore, compared to the most frequently detected TMPRSS2-ERG form (e1e4), e2e4 encodes 31 more amino acids and accelerated neuroendocrine process of prostate cancer.

Conclusions

In summary, these findings painted the landscape of chimeric RNA in NEPC and supported the idea that some chimeric RNAs may represent additional biomarkers and/or treatment targets independent of parental gene transcripts.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00893-5.

A diversified and segregated mRNA spliced-leader system in the parasitic Perkinsozoa

Spliced-leader trans-splicing (SLTS) has been described in distantly related eukaryotes and acts to mark mRNAs with a short 5′ exon, giving different mRNAs identical 5′ sequence-signatures. The function of these systems is obscure. Perkinsozoa encompasses a diversity of parasitic protists that infect bivalves, toxic-tide dinoflagellates, fish and frog tadpoles. Here, we report considerable sequence variation in the SLTS-system across the Perkinsozoa and find that multiple variant SLTS-systems are encoded in parallel in the ecologically important Perkinsozoa parasite Parvilucifera sinerae. These results demonstrate that the transcriptome of P. sinerae is segregated based on the addition of different spliced-leader (SL) exons. This segregation marks different gene categories, suggesting that SL-segregation relates to functional differentiation of the transcriptome. By contrast, both sets of gene categories are present in the single SL-transcript type sampled from Maranthos, implying that the SL-segregation of the Parvilucifera transcriptome is a recent evolutionary innovation. Furthermore, we show that the SLTS-system marks a subsection of the transcriptome with increased mRNA abundance and includes genes that encode the spliceosome system necessary for SLTS-function. Collectively, these data provide a picture of how the SLTS-systems can vary within a major evolutionary group and identify how additional transcriptional-complexity can be achieved through SL-segregation.

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.

Trans-splicing facilitated by RNA pairing greatly expands sDscam isoform diversity but not homophilic binding specificity

The Down syndrome cell adhesion molecule 1 (Dscam1) gene can generate tens of thousands of isoforms via alternative splicing, which is essential for nervous and immune functions. Chelicerates generate approximately 50 to 100 shortened Dscam (sDscam) isoforms by alternative promoters, similar to mammalian protocadherins. Here, we reveal that trans-splicing markedly increases the repository of sDscamβ isoforms in Tetranychus urticae. Unexpectedly, every variable exon cassette engages in trans-splicing with constant exons from another cluster. Moreover, we provide evidence that competing RNA pairing not only governs alternative cis-splicing but also facilitates trans-splicing. Trans-spliced sDscam isoforms mediate cell adhesion ability but exhibit the same homophilic binding specificity as their cis-spliced counterparts. Thus, we reveal a single sDscam locus that generates diverse adhesion molecules through cis- and trans-splicing coupled with alternative promoters. These findings expand understanding of the mechanism underlying molecular diversity and have implications for the molecular control of neuronal and/or immune specificity.

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.

Genion, an accurate tool to detect gene fusion from long transcriptomics reads

BACKGROUND

The advent of next-generation sequencing technologies empowered a wide variety of transcriptomics studies. A widely studied topic is gene fusion which is observed in many cancer types and suspected of having oncogenic properties. Gene fusions are the result of structural genomic events that bring two genes closely located and result in a fused transcript. This is different from fusion transcripts created during or after the transcription process. These chimeric transcripts are also known as read-through and trans-splicing transcripts. Gene fusion discovery with short reads is a well-studied problem, and many methods have been developed. But the sensitivity of these methods is limited by the technology, especially the short read length. Advances in long-read sequencing technologies allow the generation of long transcriptomics reads at a low cost. Transcriptomic long-read sequencing presents unique opportunities to overcome the shortcomings of short-read technologies for gene fusion detection while introducing new challenges.

RESULTS

We present Genion, a sensitive and fast gene fusion detection method that can also detect read-through events. We compare Genion against a recently introduced long-read gene fusion discovery method, LongGF, both on simulated and real datasets. On simulated data, Genion accurately identifies the gene fusions and its clustering accuracy for detecting fusion reads is better than LongGF. Furthermore, our results on the breast cancer cell line MCF-7 show that Genion correctly identifies all the experimentally validated gene fusions.

CONCLUSIONS

Genion is an accurate gene fusion caller. Genion is implemented in C++ and is available at https://github.com/vpc-ccg/genion .

Computational analysis of spliced leader trans-splicing in the regenerative flatworm Macrostomum lignano reveals its prevalence in conserved and stem cell related genes

In eukaryotes, trans-splicing is a process of nuclear pre-mRNA maturation where two different RNA molecules are joined together by the spliceosomal machinery utilizing mechanisms similar to cis-splicing. In diverse taxa of lower eukaryotes, spliced leader (SL) trans-splicing is the most frequent type of trans-splicing, when the same sequence derived from short small nuclear RNA molecules, called SL RNAs, is attached to the 5’ ends of different non-processed pre-mRNAs. One of the functions of SL trans-splicing is processing polycistronic pre-mRNA molecules transcribed from operons, when several genes are transcribed as one pre-mRNA molecule. However, only a fraction of trans-spliced genes reside in operons, suggesting that SL trans-splicing must also have some other, less understood functions. Regenerative flatworms are informative model organisms which hold the keys to understand the mechanism of stem cell regulation and specialization during regeneration and homeostasis. Their ability to regenerate is fueled by the division and differentiation of the adult somatic stem cell population called neoblasts. Macrostomum lignano is a flatworm model organism where substantial technological advances have been achieved in recent years, including the development of transgenesis. Although a large fraction of genes in M. lignano were estimated to be SL trans-spliced, SL trans-splicing was not studied in detail in M. lignano before. Here, we performed the first comprehensive study of SL trans-splicing in M. lignano. By reanalyzing the existing genome and transcriptome data of M. lignano, we estimate that 30 % of its genes are SL trans-spliced, 15 % are organized in operons, and almost 40 % are both SL trans-spliced and in operons. We annotated and characterized the sequence of SL RNA and characterized conserved cis- and SL transsplicing motifs. Finally, we found that a majority of SL trans-spliced genes are evolutionarily conserved and significantly over-represented in neoblast-specific genes. Our findings suggest an important role of SL trans-splicing in the regulation and maintenance of neoblasts in M. lignano.

Monodopsis and Vischeria Genomes Shed New Light on the Biology of Eustigmatophyte Algae

Members of eustigmatophyte algae, especially Nannochloropsis and Microchloropsis, have been tapped for biofuel production owing to their exceptionally high lipid content. Although extensive genomic, transcriptomic, and synthetic biology toolkits have been made available for Nannochloropsis and Microchloropsis, very little is known about other eustigmatophytes. Here we present three near-chromosomal and gapless genome assemblies of Monodopsis strains C73 and C141 (60 Mb) and Vischeria strain C74 (106 Mb), which are the sister groups to Nannochloropsis and Microchloropsis in the order Eustigmatales. These genomes contain unusually high percentages of simple repeats, ranging from 12% to 21% of the total assembly size. Unlike Nannochloropsis and Microchloropsis, long interspersed nuclear element repeats are abundant in Monodopsis and Vischeria and might constitute the centromeric regions. We found that both mevalonate and nonmevalonate pathways for terpenoid biosynthesis are present in Monodopsis and Vischeria, which is different from Nannochloropsis and Microchloropsis that have only the latter. Our analysis further revealed extensive spliced leader trans-splicing in Monodopsis and Vischeria at 36-61% of genes. Altogether, the high-quality genomes of Monodopsis and Vischeria not only serve as the much-needed outgroups to advance Nannochloropsis and Microchloropsis research, but also shed new light on the biology and evolution of eustigmatophyte algae.

Isoform-specific roles of the Drosophila filamin-type protein Jitterbug (Jbug) during development

Filamins are highly conserved actin-crosslinking proteins that regulate organization of the actin cytoskeleton. As key components of versatile signaling scaffolds, filamins are implicated in developmental anomalies and cancer. Multiple isoforms of filamins exist, raising the possibility of distinct functions for each isoform during development and in disease. Here, we provide an initial characterization of jitterbug (jbug), which encodes one of the two filamin-type proteins in Drosophila. We generate Jbug antiserum that recognizes all of the spliced forms and reveals differential expression of different Jbug isoforms during development, and a significant maternal contribution of Jbug protein. To reveal the function of Jbug isoforms, we create new genetic tools, including a null allele that deletes all isoforms, hypomorphic alleles that affect only a subset, and UAS lines for Gal4-driven expression of the major isoforms. Using these tools, we demonstrate that Jbug is required for viability and that specific isoforms are required in the formation of actin-rich protrusions including thoracic bristles in adults and ventral denticles in the embryo. We also show that specific isoforms of Jbug show differential localization within epithelia and that maternal and zygotic loss of jbug disrupts Crumbs (Crb) localization in several epithelial cell types.

LYSINE KETOGLUTARATE REDUCTASE TRANS-SPLICING RELATED 1 is involved in temperature-dependent root growth in rice

Root length is an important root parameter directly related to the uptake of water and nutrients. However, the molecular mechanisms controlling root length are still not fully understood. Here, we isolated a short-root mutant of rice, dice2 (defective in cell elongation 2). The cell length and meristem size of the roots were decreased in dice2, but the root function in terms of mineral element uptake, root cell width, and root anatomy were hardly altered compared with wild-type (WT) rice. The root growth defect in dice2 could be partially rescued by high temperature. Map-based cloning combined with a complementation test revealed that the short-root phenotype was caused by a nonsense mutation in a gene which was annotated to encode Lysine Ketoglutarate Reductase Trans-Splicing related 1 (OsLKRT1). OsLKRT1, encoding a cytosol-localized protein, was expressed in all cells of the root tip and elongation region as well as the shoot. RNA-seq analysis showed that there was no difference between dice2 and the WT in the expression level of genes involved in root development identified so far. These results indicate that OsLKRT1 is involved in a novel pathway required for root cell elongation in rice, although its exact role remains to be further investigated.

Widespread polycistronic gene expression in green algae

Historically, it has been understood that for gene expression in eukaryotes, each messenger RNA encodes a single protein. With the recent development of technologies to sequence full-length transcripts en masse, we have discovered hundreds of examples in two species of green algae where two, three, or more proteins are translated from a single transcript. These “polycistronic” transcripts are found in diverse species throughout the green algal lineage, which highlights their biological importance. We have leveraged these findings to coexpress pairs of genes on polycistronic transcripts in vitro, which should facilitate efforts to engineer algae for research and industrial applications.

Polycistronic gene expression, common in prokaryotes, was thought to be extremely rare in eukaryotes. The development of long-read sequencing of full-length transcript isomers (Iso-Seq) has facilitated a reexamination of that dogma. Using Iso-Seq, we discovered hundreds of examples of polycistronic expression of nuclear genes in two divergent species of green algae: Chlamydomonas reinhardtii and Chromochloris zofingiensis. Here, we employ a range of independent approaches to validate that multiple proteins are translated from a common transcript for hundreds of loci. A chromatin immunoprecipitation analysis using trimethylation of lysine 4 on histone H3 marks confirmed that transcription begins exclusively at the upstream gene. Quantification of polyadenylated [poly(A)] tails and poly(A) signal sequences confirmed that transcription ends exclusively after the downstream gene. Coexpression analysis found nearly perfect correlation for open reading frames (ORFs) within polycistronic loci, consistent with expression in a shared transcript. For many polycistronic loci, terminal peptides from both ORFs were identified from proteomics datasets, consistent with independent translation. Synthetic polycistronic gene pairs were transcribed and translated in vitro to recapitulate the production of two distinct proteins from a common transcript. The relative abundance of these two proteins can be modified by altering the Kozak-like sequence of the upstream gene. Replacement of the ORFs with selectable markers or reporters allows production of such heterologous proteins, speaking to utility in synthetic biology approaches. Conservation of a significant number of polycistronic gene pairs between C. reinhardtii, C. zofingiensis, and five other species suggests that this mechanism may be evolutionarily ancient and biologically important in the green algal lineage.

Co-Regulated Genes and Gene Clusters

There are many co-regulated genes in eukaryotic cells. The coordinated activation or repression of such genes occurs at specific stages of differentiation, or under the influence of external stimuli. As a rule, co-regulated genes are dispersed in the genome. However, there are also gene clusters, which contain paralogous genes that encode proteins with similar functions. In this aspect, they differ significantly from bacterial operons containing functionally linked genes that are not paralogs. In this review, we discuss the reasons for the existence of gene clusters in vertebrate cells and propose that clustering is necessary to ensure the possibility of selective activation of one of several similar genes.

Small-RNA Sequencing Reveals Altered Skeletal Muscle microRNAs and snoRNAs Signatures in Weanling Male Offspring from Mouse Dams Fed a Low Protein Diet during Lactation

Maternal diet during gestation and lactation affects the development of skeletal muscles in offspring and determines muscle health in later life. In this paper, we describe the association between maternal low protein diet-induced changes in offspring skeletal muscle and the differential expression (DE) of small non-coding RNAs (sncRNAs). We used a mouse model of maternal protein restriction, where dams were fed either a normal (N, 20%) or a low protein (L, 8%) diet during gestation and newborns were cross-fostered to N or L lactating dams, resulting in the generation of NN, NL and LN offspring groups. Total body and tibialis anterior (TA) weights were decreased in weanling NL male offspring but were not different in the LN group, as compared to NN. However, histological evaluation of TA muscle revealed reduced muscle fibre size in both groups at weaning. Small RNA-sequencing demonstrated DE of multiple miRs, snoRNAs and snRNAs. Bioinformatic analyses of miRs-15a, -34a, -122 and -199a, in combination with known myomiRs, confirmed their implication in key muscle-specific biological processes. This is the first comprehensive report for the DE of sncRNAs in nutrition-associated programming of skeletal muscle development, highlighting the need for further research to unravel the detailed molecular mechanisms.

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.

Rice FLOURY SHRUNKEN ENDOSPERM 5 Encodes a Putative Plant Organelle RNA Recognition Protein that Is Required for cis-Splicing of Mitochondrial nad4 Intron 1

BACKGROUND

The sequences of several important mitochondrion-encoded genes involved in respiration in higher plants are interrupted by introns. Many nuclear-encoded factors are involved in splicing these introns, but the mechanisms underlying this splicing remain unknown.

RESULTS

We isolated and characterized a rice mutant named floury shrunken endosperm 5 (fse5). In addition to having floury shrunken endosperm, the fse5 seeds either failed to germinate or produced seedlings which grew slowly and died ultimately. Fse5 encodes a putative plant organelle RNA recognition (PORR) protein targeted to mitochondria. Mutation of Fse5 hindered the splicing of the first intron of nad4, which encodes an essential subunit of mitochondrial NADH dehydrogenase complex I. The assembly and NADH dehydrogenase activity of complex I were subsequently disrupted by this mutation, and the structure of the mitochondria was abnormal in the fse5 mutant. The FSE5 protein was shown to interact with mitochondrial intron splicing factor 68 (MISF68), which is also a splicing factor for nad4 intron 1 identified previously via yeast two-hybrid (Y2H) assays.

CONCLUSION

Fse5 which encodes a PORR domain-containing protein, is essential for the splicing of nad4 intron 1, and loss of Fse5 function affects seed development and seedling growth.

Rapid protein evolution, organellar reductions, and invasive intronic elements in the marine aerobic parasite dinoflagellate Amoebophrya spp

BACKGROUND

Dinoflagellates are aquatic protists particularly widespread in the oceans worldwide. Some are responsible for toxic blooms while others live in symbiotic relationships, either as mutualistic symbionts in corals or as parasites infecting other protists and animals. Dinoflagellates harbor atypically large genomes (~ 3 to 250 Gb), with gene organization and gene expression patterns very different from closely related apicomplexan parasites. Here we sequenced and analyzed the genomes of two early-diverging and co-occurring parasitic dinoflagellate Amoebophrya strains, to shed light on the emergence of such atypical genomic features, dinoflagellate evolution, and host specialization.

RESULTS

We sequenced, assembled, and annotated high-quality genomes for two Amoebophrya strains (A25 and A120), using a combination of Illumina paired-end short-read and Oxford Nanopore Technology (ONT) MinION long-read sequencing approaches. We found a small number of transposable elements, along with short introns and intergenic regions, and a limited number of gene families, together contribute to the compactness of the Amoebophrya genomes, a feature potentially linked with parasitism. While the majority of Amoebophrya proteins (63.7% of A25 and 59.3% of A120) had no functional assignment, we found many orthologs shared with Dinophyceae. Our analyses revealed a strong tendency for genes encoded by unidirectional clusters and high levels of synteny conservation between the two genomes despite low interspecific protein sequence similarity, suggesting rapid protein evolution. Most strikingly, we identified a large portion of non-canonical introns, including repeated introns, displaying a broad variability of associated splicing motifs never observed among eukaryotes. Those introner elements appear to have the capacity to spread over their respective genomes in a manner similar to transposable elements. Finally, we confirmed the reduction of organelles observed in Amoebophrya spp., i.e., loss of the plastid, potential loss of a mitochondrial genome and functions.

CONCLUSION

These results expand the range of atypical genome features found in basal dinoflagellates and raise questions regarding speciation and the evolutionary mechanisms at play while parastitism was selected for in this particular unicellular lineage.

Evaluation of various RNA-seq approaches for identification of gene outrons in the flatworm Opisthorchis felineus

The parasitic flatworm Opisthorchis felineus is one of the causative agents of opisthorchiasis in humans. Recently, we assembled the O. felineus genome, but the correct genome annotation by means of standard methods was hampered by the presence of spliced leader trans-splicing (SLTS). As a result of SLTS, the original 5’-end (outron) of the transcripts is replaced by a short spliced leader sequence donated from a specialized SL RNA. SLTS is involved in the RNA processing of more than half of O. felineus genes, making it hard to determine the structure of outrons and bona fide transcription start sites of the corresponding genes and operons, being based solely on mRNA-seq data. In the current study, we tested various experimental approaches for identifying the sequences of outrons in O. felineus using massive parallel sequencing. Two of them were developed by us for targeted sequencing of already processed branched outrons. One was based on sequence-specific reverse transcription from the SL intron toward the 5’-end of the Y-branched outron. The other used outron hybridization with an immobilized single-stranded DNA probe complementary to the SL intron. Additionally, two approaches to the sequencing of rRNA-depleted total RNA were used, allowing the identification of a wider range of transcripts compared to mRNAseq. One is based on the enzymatic elimination of overrepresented cDNAs, the other utilizes exonucleolytic degradation of uncapped RNA by Terminator enzyme. By using the outron-targeting methods, we were not able to obtain the enrichment of RNA preparations by processed outrons, which is most likely indicative of a rapid turnover of these trans-splicing intermediate products. Of the two rRNA depletion methods, a method based on the enzymatic normalization of cDNA (Zymo-Seq RiboFree) showed high efficiency. Compared to mRNA-seq, it provides an approximately twofold increase in the fraction of reads originating from outrons and introns. The results suggest that unprocessed nascent transcripts are the main source of outron sequences in the RNA pool of O. felineus.

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.

How to determine the mechanism of action of CFTR modulator compounds: A gateway to theranostics

The greatest challenge of 21st century biology is to fully understand mechanisms of disease to drive new approaches and medical innovation. Parallel to this is the huge biomedical endeavour of treating people through personalized medicine. Until now all CFTR modulator drugs that have entered clinical trials have been genotype-dependent. An emerging alternative is personalized/precision medicine in CF, i.e., to determine whether rare CFTR mutations respond to existing (or novel) CFTR modulator drugs by pre-assessing them directly on patient's tissues ex vivo, an approach also now termed theranostics. To administer the right drug to the right person it is essential to understand how drugs work, i.e., to know their mechanism of action (MoA), so as to predict their applicability, not just in certain mutations but also possibly in other diseases that share the same defect/defective pathway. Moreover, an understanding the MoA of a drug before it is tested in clinical trials is the logical path to drug discovery and can increase its chance for success and hence also approval. In conclusion, the most powerful approach to determine the MoA of a compound is to understand the underlying biology. Novel large datasets of intervenients in most biological processes, namely those emerging from the post-genomic era tools, are available and should be used to help in this task.

Loss of the RNA trimethylguanosine cap is compatible with nuclear accumulation of spliceosomal snRNAs but not pre-mRNA splicing or snRNA processing during animal development

The 2,2,7-trimethylguanosine (TMG) cap is one of the first identified modifications on eukaryotic RNAs. TMG, synthesized by the conserved Tgs1 enzyme, is abundantly present on snRNAs essential for pre-mRNA splicing. Results from ex vivo experiments in vertebrate cells suggested that TMG ensures nuclear localization of snRNAs. Functional studies of TMG using tgs1 mutations in unicellular organisms yield results inconsistent with TMG being indispensable for either nuclear import or splicing. Utilizing a hypomorphic tgs1 mutation in Drosophila, we show that TMG reduction impairs germline development by disrupting the processing, particularly of introns with smaller sizes and weaker splice sites. Unexpectedly, loss of TMG does not disrupt snRNAs localization to the nucleus, disputing an essential role of TMG in snRNA transport. Tgs1 loss also leads to defective 3' processing of snRNAs. Remarkably, stronger tgs1 mutations cause lethality without severely disrupting splicing, likely due to the preponderance of TMG-capped snRNPs. Tgs1, a predominantly nucleolar protein in Drosophila, likely carries out splicing-independent functions indispensable for animal development. Taken together, our results suggest that nuclear import is not a conserved function of TMG. As a distinctive structure on RNA, particularly non-coding RNA, we suggest that TMG prevents spurious interactions detrimental to the function of RNAs that it modifies.

SLFinder, a pipeline for the novel identification of splice-leader sequences: a good enough solution for a complex problem

BACKGROUND

Spliced Leader trans-splicing is an important mechanism for the maturation of mRNAs in several lineages of eukaryotes, including several groups of parasites of great medical and economic importance. Nevertheless, its study across the tree of life is severely hindered by the problem of identifying the SL sequences that are being trans-spliced.

RESULTS

In this paper we present SLFinder, a four-step pipeline meant to identify de novo candidate SL sequences making very few assumptions regarding the SL sequence properties. The pipeline takes transcriptomic de novo assemblies and a reference genome as input and allows the user intervention on several points to account for unexpected features of the dataset. The strategy and its implementation were tested on real RNAseq data from species with and without SL Trans-Splicing.

CONCLUSIONS

SLFinder is capable to identify SL candidates with good precision in a reasonable amount of time. It is especially suitable for species with unknown SL sequences, generating candidate sequences for further refining and experimental validation.

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.

Evidence of constraint in the 3D genome for trans-splicing in human cells

Fusion transcripts are commonly found in eukaryotes, and many aberrant fusions are associated with severe diseases, including cancer. One class of fusion transcripts is generated by joining separate transcripts through trans-splicing. However, the mechanism of trans-splicing in mammals remains largely elusive. Here we showed evidence to support an intuitive hypothesis that attributes trans-sphcing to the spatial proximity between premature transcripts. A novel trans-splicing detection tool (TSD) was developed to reliably identify intra-chromosomal trans-splicing events (iTSEs) from RNA-seq data. TSD can maintain a remarkable balance between sensitivity and accuracy, thus distinguishing it from most state-of-the-art tools. The accuracy of TSD was experimentally demonstrated by excluding potential false discovery from mosaic genome or template switching during PCR. We showed that iTSEs identified by TSD were frequently found between genomic regulatory elements, which are known to be more prone to interact with each other. Moreover, iTSE sites may be more physically adjacent to each other than random control in the tested human lymphoblastoid cell line according to Hi-C data. Our results suggest that trans-splicing and 3D genome architecture may be coupled in mammals and that our pipeline, TSD, may facilitate investigations of trans-splicing on a systematic and accurate level previously thought impossible.

Gene fragmentation and RNA editing without borders: eccentric mitochondrial genomes of diplonemids

Diplonemids are highly abundant heterotrophic marine protists. Previous studies showed that their strikingly bloated mitochondrial genome is unique because of systematic gene fragmentation and manifold RNA editing. Here we report a comparative study of mitochondrial genome architecture, gene structure and RNA editing of six recently isolated, phylogenetically diverse diplonemid species. Mitochondrial gene fragmentation and modes of RNA editing, which include cytidine-to-uridine (C-to-U) and adenosine-to-inosine (A-to-I) substitutions and 3' uridine additions (U-appendage), are conserved across diplonemids. Yet as we show here, all these features have been pushed to their extremes in the Hemistasiidae lineage. For example, Namystynia karyoxenos has its genes fragmented into more than twice as many modules than other diplonemids, with modules as short as four nucleotides. Furthermore, we detected in this group multiple A-appendage and guanosine-to-adenosine (G-to-A) substitution editing events not observed before in diplonemids and found very rarely elsewhere. With >1,000 sites, C-to-U and A-to-I editing in Namystynia is nearly 10 times more frequent than in other diplonemids. The editing density of 12% in coding regions makes Namystynia's the most extensively edited transcriptome described so far. Diplonemid mitochondrial genome architecture, gene structure and post-transcriptional processes display such high complexity that they challenge all other currently known systems.

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.

Detection of Group II Intron-Generated Chimeric mRNAs in Bacterial Cells

Chimeric RNAs are the transcripts composed of exons from two separate genes or transcripts. Although the presence of these joined RNA molecules have mainly been documented in a variety of eukaryotes, we recently demonstrated that the Ll.LtrB group II intron, from the gram-positive bacterium Lactococcus lactis, can generate chimeric mRNAs through a novel intergenic trans-splicing pathway. Here we describe the detailed experimental procedures to detect group II intron-generated mRNA-mRNA chimeras from total RNA extracts using stringent reverse transcription conditions along with a reverse splicing-deficient group II intron as a negative control.

An improved method for circular RNA purification using RNase R that efficiently removes linear RNAs containing G-quadruplexes or structured 3' ends

Thousands of eukaryotic protein-coding genes generate circular RNAs that have covalently linked ends and are resistant to degradation by exonucleases. To prove their circularity as well as biochemically enrich these transcripts, it has become standard in the field to use the 3′-5′ exonuclease RNase R. Here, we demonstrate that standard protocols involving RNase R can fail to digest >20% of all highly expressed linear RNAs, but these shortcomings can largely be overcome. RNAs with highly structured 3′ ends, including snRNAs and histone mRNAs, are naturally resistant to RNase R, but can be efficiently degraded once a poly(A) tail has been added to their ends. In addition, RNase R stalls in the body of many polyadenylated mRNAs, especially at G-rich sequences that have been previously annotated as G-quadruplex (G4) structures. Upon replacing K⁺ (which stabilizes G4s) with Li⁺ in the reaction buffer, we find that RNase R is now able to proceed through these sequences and fully degrade the mRNAs in their entirety. In total, our results provide important improvements to the current methods used to isolate circular RNAs as well as a way to reveal RNA structures that may naturally inhibit degradation by cellular exonucleases.

Functional characterization of the translation initiation factor eIF4E of Echinococcus granulosus

The eukaryotic initiation factor 4E (eIF4E) specifically recognizes the 5' mRNA cap, a rate-limiting step in the translation initiation process. Although the 7-methylguanosine cap (MMGcap) is the most common 5' cap structure in eukaryotes, the trans-splicing process that occurs in several organism groups, including nematodes and flatworms, leads to the addition of a trimethylguanosine cap (TMGcap) to some RNA transcripts. In some helminths, eIF4E can have a dual capacity to bind both MMGcap and TMGcap. In the present work, we evaluated the distribution of eIF4E protein sequences in platyhelminths and we showed that only one gene coding for eIF4E is present in most parasitic flatworms. Based on this result, we cloned the Echinococcus granulosus cDNA sequence encoding eIF4E in Escherichia coli, expressed the recombinant eIF4E as a fusion protein to GST, and tested its ability to capture mRNAs through the 5' cap using pull-down assay and qPCR. Our results indicate that the recombinant eIF4E was able to bind preferentially 5'-capped mRNAs compared with rRNAs from total RNA preparations of E. granulosus. By qPCR, we observed an enrichment in MMG-capped mRNA compared with TMG-capped mRNAs among Eg-eIF4E-GST pull-down RNAs. Eg-eIF4E structural model using the Schistosoma mansoni eIF4E as template showed to be well preserved with only a few differences between chemically similar amino acid residues at the binding sites. These data showed that E. granulosus eIF4E can be used as a potential tool to study full-length 5'-capped mRNA, besides being a potential drug target against parasitic flatworms.

A high-throughput screen for the identification of compounds that inhibit nematode gene expression by targeting spliced leader trans-splicing

Infections with parasitic nematodes are among the most significant of the neglected tropical diseases affecting about a billion people living mainly in tropical regions with low economic activity. The most effective current strategy to control nematode infections involves large scale treatment programs with anthelmintic drugs. This strategy is at risk from the emergence of drug resistant parasites. Parasitic nematodes also affect livestock, which are treated with the same limited group of anthelmintic drugs. Livestock parasites resistant to single drugs, and even multi-drug resistant parasites, are appearing in many areas. There is therefore a pressing need for new anthelmintic drugs. Here we use the nematode Caenorhabditis elegans as a model for parasitic nematodes and demonstrate that sinefungin, a competitive inhibitor of methyltransferases, causes a delay in development and reduced fecundity, and inhibits spliced leader trans-splicing. Spliced leader trans-splicing is an essential step in gene expression that does not occur in the hosts of parasitic nematodes, and is therefore a potential target for new anthelmintic drugs. We have exploited the ability of sinefungin to inhibit spliced leader trans-splicing to adapt a green fluorescent protein based reporter gene assay that monitors spliced leader trans-splicing for high-throughput screening for new anthelmintic compounds. We have established a protocol for robust high-throughput screening, combining mechanical dispensing of living C. elegans into 384- or 1536- well plates with addition of compounds using an acoustic liquid dispenser, and the detection of the inhibition of SL trans-splicing using a microplate reader. We have tested this protocol in a first pilot screen and envisage that this assay will be a valuable tool in the search for new anthelmintic drugs.

Satellite DNA-containing gigantic introns in a unique gene expression program during Drosophila spermatogenesis

Intron gigantism, where genes contain megabase-sized introns, is observed across species, yet little is known about its purpose or regulation. Here we identify a unique gene expression program utilized for the proper expression of genes with intron gigantism. We find that two Drosophila genes with intron gigantism, kl-3 and kl-5, are transcribed in a spatiotemporal manner over the course of spermatocyte differentiation, which spans ~90 hours. The introns of these genes contain megabases of simple satellite DNA repeats that comprise over 99% of the gene loci, and these satellite-DNA containing introns are transcribed. We identify two RNA-binding proteins that specifically localize to kl-3 and kl-5 transcripts and are needed for the successful transcription or processing of these genes. We propose that genes with intron gigantism require a unique gene expression program, which may serve as a platform to regulate gene expression during cellular differentiation.

Introns are non-coding elements of eukaryotic genes, often containing important regulatory sequences. Curiously, some genes contain introns so large that more than 99% of the gene locus is non-coding. One of the best-studied large genes, Dystrophin, a causative gene for Duchenne Muscular Dystrophy, spans 2.2Mb, only 11kb of which is coding. This phenomenon, ‘intron gigantism’, is observed across species, yet little is known about its purpose or regulation. Here we identify a unique gene expression program utilized for the proper expression of genes with intron gigantism using Drosophila spermatogenic genes as a model system. We show that the gigantic introns of these genes are transcribed in line with the exons, likely as a single transcript. We identify two RNA-binding proteins that specifically localize to the site of transcription and are needed for the successful transcription or processing of these genes. We propose that genes with intron gigantism require a unique gene expression program, which may serve as a platform to regulate gene expression during cellular differentiation.

Disruption of the Caenorhabditis elegans Integrator complex triggers a non-conventional transcriptional mechanism beyond snRNA genes

Gene expression is generally regulated by recruitment of transcription factors and RNA polymerase II (RNAP II) to specific sequences in the gene promoter region. The Integrator complex mediates processing of small nuclear RNAs (snRNAs) as well as the initiation and release of paused RNAP II at specific genes in response to growth factors. Here we show that in C. elegans, disruption of the Integrator complex leads to transcription of genes located downstream of the snRNA loci via a non-conventional transcription mechanism based on the lack of processing of the snRNAs. RNAP II read-through generates long chimeric RNAs containing snRNA, the intergenic region and the mature mRNA of the downstream gene located in sense. These chimeric sn-mRNAs remain as untranslated long non-coding RNAs, in the case of U1- and U2-derived sn-mRNAs, but can be translated to proteins in the case of SL-derived sn-mRNAs. The transcriptional effect caused by disruption of the Integrator complex is not restricted to genes located downstream of the snRNA loci but also affects key regulators of signal transduction such as kinases and phosphatases. Our findings highlight that these transcriptional alterations may be behind the correlation between mutations in the Integrator complex and tumor transformation.

CRISPR/Cas9-mediated genome editing induces gene knockdown by altering the pre-mRNA splicing in mice

Background

Clustered Regularly Interspaced Short Palindromic Repeats/CRISPR associated protein 9 (CRISPR/Cas9) has been wildly used to generate gene knockout models through inducing indels causing frame-shift. However, there are few studies concerning the post-transcript effects caused by CRISPR-mediated genome editing.

Results

In the present study, we showed that gene knockdown model also could be generated using CRISPR-mediated gene editing by disrupting the boundary of exon and intron in mice (C57BL/6 J). CRISPR induced indel at the boundary of exon and intron (5′ splice site) caused alternative splicing and produced multiple different mRNAs, most of these mRNAs introduced premature termination codon causing down expression of the gene.

Conclusions

These results showed that alternative splicing mutants were able to generate through CRISPR-mediated genome editing by deleting the boundary of exon and intron causing disruption of 5′ splice site. Although alternative splicing was an unexpected outcome, this finding could be developed as a technology to generate gene knockdown models or to investigate pre-mRNA splicing.

Electronic supplementary material

The online version of this article (10.1186/s12896-018-0472-8) contains supplementary material, which is available to authorized users.

Characterization of spliced leader trans-splicing in a photosynthetic rhizarian amoeba, Paulinella micropora, and its possible role in functional gene transfer

Paulinella micropora is a rhizarian thecate amoeba, belonging to a photosynthetic Paulinella species group that has a unique organelle termed chromatophore, whose cyanobacterial origin is distinct from that of plant and algal chloroplasts. Because acquisition of the chromatophore was quite a recent event compared with that of the chloroplast ancestor, the Paulinella species are thought to be model organisms for studying the early process of primary endosymbiosis. To obtain insight into how endosymbiotically transferred genes acquire expression competence in the host nucleus, here we analyzed the 5′ end sequences of the mRNAs of P. micropora MYN1 strain with the aid of a cap-trapper cDNA library. As a result, we found that mRNAs of 27 genes, including endosymbiotically transferred genes, possessed the common 5′ end sequence of 28–33 bases that were posttranscriptionally added by spliced leader (SL) trans-splicing. We also found two subtypes of SL RNA genes encoded by the P. micropora MYN1 genome. Differing from the other SL trans-splicing organisms that usually possess poly(A)-less SL RNAs, this amoeba has polyadenylated SL RNAs. In this study, we characterize the SL trans-splicing of this unique organism and discuss the putative merits of SL trans-splicing in functional gene transfer and genome evolution.

Translation and Translational Control in Dinoflagellates

Dinoflagellates are unicellular protists that feature a multitude of unusual nuclear features, including large genomes, packaging of DNA without histones, and multiple gene copies organized as tandem gene arrays. Furthermore, all dinoflagellate mRNAs experience trans-splicing with a common 22-nucleotide splice leader (SL) sequence. These features challenge some of the concepts and assumptions about the regulation of gene expression derived from work on model eukaryotes such as yeasts and mammals. Translational control in the dinoflagellates, based on extensive study of circadian bioluminescence and by more recent microarray and transcriptome analyses, is now understood to be a crucial element in regulating gene expression. A picture of the translation machinery of dinoflagellates is emerging from the recent availability of transcriptomes of multiple dinoflagellate species and the first complete genome sequences. The components comprising the translational control toolkit of dinoflagellates are beginning to take shape and are outlined here.

Landscape of the spliced leader trans-splicing mechanism in Schistosoma mansoni

Spliced leader dependent trans-splicing (SLTS) has been described as an important RNA regulatory process that occurs in different organisms, including the trematode Schistosoma mansoni. We identified more than seven thousand putative SLTS sites in the parasite, comprising genes with a wide spectrum of functional classes, which underlines the SLTS as a ubiquitous mechanism in the parasite. Also, SLTS gene expression levels span several orders of magnitude, showing that SLTS frequency is not determined by the expression level of the target gene, but by the presence of particular gene features facilitating or hindering the trans-splicing mechanism. Our in-depth investigation of SLTS events demonstrates widespread alternative trans-splicing (ATS) acceptor sites occurring in different regions along the entire gene body, highlighting another important role of SLTS generating alternative RNA isoforms in the parasite, besides the polycistron resolution. Particularly for introns where SLTS directly competes for the same acceptor substrate with cis-splicing, we identified for the first time additional and important features that might determine the type of splicing. Our study substantially extends the current knowledge of RNA processing by SLTS in S. mansoni, and provide basis for future studies on the trans-splicing mechanism in other eukaryotes.