A global analysis of C. elegans trans-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.

Binding profiles for 954 Drosophila and C. elegans transcription factors reveal tissue specific regulatory relationships

A catalog of transcription factor (TF) binding sites in the genome is critical for deciphering regulatory relationships. Here we present the culmination of the modERN (model organism Encyclopedia of Regulatory Networks) consortium that systematically assayed TF binding events in vivo in two major model organisms, Drosophila melanogaster (fly) and Caenorhabditis elegans (worm). We describe key features of these datasets, comprising 604 TFs identifying 3.6M sites in the fly and 350 TFs identifying 0.9 M sites in the worm. Applying a machine learning model to these data identifies sets of TFs with a prominent role in promoting target gene expression in specific cell types. TF binding data are available through the ENCODE Data Coordinating Center and at https://epic.gs.washington.edu/modERNresource, which provides access to processed and summary data, as well as widgets to probe cell type-specific TF-target relationships. These data are a rich resource that should fuel investigations into TF function during development.

Transcript accumulation rates in the early Caenorhabditis elegans embryo

Dynamic transcriptional changes are widespread in rapidly dividing developing embryos when cell fate decisions are made quickly. The Caenorhabditis elegans embryo overcomes these constraints partly through the rapid production of high levels of transcription factor mRNAs. Transcript accumulation rates for some developmental genes are known at single-cell resolution, but genome-scale measurements are lacking. We estimate zygotic mRNA accumulation rates from single-cell RNA sequencing data calibrated with single-molecule transcript imaging. Rapid transcription is common in the early C. elegans embryo with rates highest soon after zygotic transcription begins. High-rate genes are enriched for recently duplicated cell-fate regulators and share common genomic features. We identify core promoter elements associated with high rate and measure their contributions for two early endomesodermal genes, ceh-51 and sdz-31. Individual motifs modestly affect accumulation rates, suggesting multifactorial control. These results are a step toward estimating absolute transcription kinetics and understanding how transcript dosage drives developmental decisions.

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.

Fluctuations in chromatin state at regulatory loci occur spontaneously under relaxed selection and are associated with epigenetically inherited variation in C. elegans gene expression

Some epigenetic information can be transmitted between generations without changes in the underlying DNA sequence. Changes in epigenetic regulators, termed epimutations, can occur spontaneously and be propagated in populations in a manner reminiscent of DNA mutations. Small RNA-based epimutations occur in C. elegans and persist for around 3-5 generations on average. Here, we explored whether chromatin states also undergo spontaneous change and whether this could be a potential alternative mechanism for transgenerational inheritance of gene expression changes. We compared the chromatin and gene expression profiles at matched time points from three independent lineages of C. elegans propagated at minimal population size. Spontaneous changes in chromatin occurred in around 1% of regulatory regions each generation. Some were heritable epimutations and were significantly enriched for heritable changes in expression of nearby protein-coding genes. Most chromatin-based epimutations were short-lived but a subset had longer duration. Genes subject to long-lived epimutations were enriched for multiple components of xenobiotic response pathways. This points to a possible role for epimutations in adaptation to environmental stressors.

Trans-splicing in the cestode Hymenolepis microstoma is constitutive across the life cycle and depends on gene structure and composition

Spliced leader (SL) trans-splicing is a key process during mRNA maturation of many eukaryotes, in which a short sequence (SL) is transferred from a precursor SL-RNA into the 5' region of an immature mRNA. This mechanism is present in flatworms, in which it is known to participate in the resolution of polycistronic transcripts. However, most trans-spliced transcripts are not part of operons, and it is not clear if this process may participate in additional regulatory mechanisms in this group. In this work, we present a comprehensive analysis of SL trans-splicing in the model cestode Hymenolepis microstoma. We identified four different SL-RNAs which are indiscriminately trans-spliced to 622 gene models. SL trans-splicing is enriched in constitutively expressed genes and does not appear to be regulated throughout the life cycle. Operons represented at least 20% of all detected trans-spliced gene models, showed conservation to those of the cestode Echinococcus multilocularis, and included complex loci such as an alternative operon (processed as either a single gene through cis-splicing or as two genes of a polycistron). Most insertion sites were identified in the 5' untranslated region (UTR) of monocistronic genes. These genes frequently contained introns in the 5' UTR, in which trans-splicing used the same acceptor sites as cis-splicing. These results suggest that, unlike other eukaryotes, trans-splicing is associated with internal intronic promoters in the 5' UTR, resulting in transcripts with strong splicing acceptor sites without competing cis-donor sites, pointing towards a simple mechanism driving the evolution of novel SL insertion sites.

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.

Widespread transposon co-option in the Caenorhabditis germline regulatory network

The movement of selfish DNA elements can lead to widespread genomic alterations with potential to create novel functions. We show that transposon expansions in Caenorhabditis nematodes led to extensive rewiring of germline transcriptional regulation. We find that about one-third of Caenorhabditis elegans germline-specific promoters have been co-opted from two related miniature inverted repeat transposable elements (TEs), CERP2 and CELE2. These promoters are regulated by HIM-17, a THAP domain-containing transcription factor related to a transposase. Expansion of CERP2 occurred before radiation of the Caenorhabditis genus, as did fixation of mutations in HIM-17 through positive selection, whereas CELE2 expanded only in C. elegans. Through comparative analyses in Caenorhabditis briggsae, we find not only evolutionary conservation of most CERP2 co-opted promoters but also a substantial fraction that are species-specific. Our work reveals the emergence and evolutionary conservation of a novel transcriptional network driven by TE co-option with a major impact on regulatory evolution.

CMTr mediated 2'-O-ribose methylation status of cap-adjacent nucleotides across animals

Cap methyltransferases (CMTrs) O methylate the 2' position of the ribose (cOMe) of cap-adjacent nucleotides of animal, protist, and viral mRNAs. Animals generally have two CMTrs, whereas trypanosomes have three, and many viruses encode one in their genome. In the splice leader of mRNAs in trypanosomes, the first four nucleotides contain cOMe, but little is known about the status of cOMe in animals. Here, we show that cOMe is prominently present on the first two cap-adjacent nucleotides with species- and tissue-specific variations in Caenorhabditis elegans, honeybees, zebrafish, mouse, and human cell lines. In contrast, Drosophila contains cOMe primarily on the first cap-adjacent nucleotide. De novo RoseTTA modeling of CMTrs reveals close similarities of the overall structure and near identity for the catalytic tetrad, and for cap and cofactor binding for human, Drosophila and C. elegans CMTrs. Although viral CMTrs maintain the overall structure and catalytic tetrad, they have diverged in cap and cofactor binding. Consistent with the structural similarity, both CMTrs from Drosophila and humans methylate the first cap-adjacent nucleotide of an AGU consensus start. Because the second nucleotide is also methylated upon heat stress in Drosophila, these findings argue for regulated cOMe important for gene expression regulation.

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.

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.

Ubiquitin ligases and a processive proteasome facilitate protein clearance during the oocyte-to-embryo transition in Caenorhabditis elegans

The ubiquitin-mediated degradation of oocyte translational regulatory proteins is a conserved feature of the oocyte-to-embryo transition. In the nematode Caenorhabditis elegans, multiple translational regulatory proteins, including the TRIM-NHL RNA-binding protein LIN-41/Trim71 and the Pumilio-family RNA-binding proteins PUF-3 and PUF-11, are degraded during the oocyte-to-embryo transition. Degradation of each protein requires activation of the M-phase cyclin-dependent kinase CDK-1, is largely complete by the end of the first meiotic division and does not require the anaphase-promoting complex. However, only LIN-41 degradation requires the F-box protein SEL-10/FBW7/Cdc4p, the substrate recognition subunit of an SCF-type E3 ubiquitin ligase. This finding suggests that PUF-3 and PUF-11, which localize to LIN-41-containing ribonucleoprotein particles, are independently degraded through the action of other factors and that the oocyte ribonucleoprotein particles are disassembled in a concerted fashion during the oocyte-to-embryo transition. We develop and test the hypothesis that PUF-3 and PUF-11 are targeted for degradation by the proteasome-associated HECT-type ubiquitin ligase ETC-1/UBE3C/Hul5, which is broadly expressed in C. elegans. We find that several GFP-tagged fusion proteins that are degraded during the oocyte-to-embryo transition, including fusions with PUF-3, PUF-11, LIN-41, IFY-1/Securin, and CYB-1/Cyclin B, are incompletely degraded when ETC-1 function is compromised. However, it is the fused GFP moiety that appears to be the critical determinant of this proteolysis defect. These findings are consistent with a conserved role for ETC-1 in promoting proteasome processivity and suggest that proteasomal processivity is an important element of the oocyte-to-embryo transition during which many key oocyte regulatory proteins are rapidly targeted for degradation.

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.

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.

The Aphelenchus avenae genome highlights evolutionary adaptation to desiccation

Some organisms can withstand complete body water loss (losing up to 99% of body water) and stay in ametabolic state for decades until rehydration, which is known as anhydrobiosis. Few multicellular eukaryotes on their adult stage can withstand life without water. We still have an incomplete understanding of the mechanism for metazoan survival of anhydrobiosis. Here we report the 255-Mb genome of Aphelenchus avenae, which can endure relative zero humidity for years. Gene duplications arose genome-wide and contributed to the expansion and diversification of 763 kinases, which represents the second largest metazoan kinome to date. Transcriptome analyses of ametabolic state of A. avenae indicate the elevation of ATP level for global recycling of macromolecules and enhancement of autophagy in the early stage of anhydrobiosis. We catalogue 74 species-specific intrinsically disordered proteins, which may facilitate A. avenae to survive through desiccation stress. Our findings refine a molecular basis evolving for survival in extreme water loss and open the way for discovering new anti-desiccation strategies.

Wan et al. report the genome and transcriptome of the Aphelenchus avenae nematode that can withstand long-term extreme desiccation. This study compares gene features to eight other nematode species and identifies intrinsically disordered proteins and changes in gene expression that contribute toward anhydrobiosis adaptation.

Cestodes in the genomic era

The first cestode genomes were obtained by an international consortium led by the Wellcome Sanger Institute that included representative institutions from countries where the sequenced parasites have been studied for decades, in part because they are etiological agents of endemic diseases (Argentina, Uruguay, Mexico, Canada, UK, Germany, Switzerland, Ireland, USA, Japan, and China). After this, several complete genomes were obtained reaching 16 species to date. Cestode genomes have smaller relative size compared to other animals including free-living flatworms. Moreover, the features genome size and repeat content seem to differ in the two analyzed orders. Cyclophyllidean species have smaller genomes and with fewer repetitive content than Diphyllobothriidean species. On average, cestode genomes have 13,753 genes with 6 exons per gene and 41% GC content. More than 5,000 shared cestode proteins were accurately annotated by the integration of gene predictions and transcriptome evidence being more than 40% of these proteins of unknown function. Several gene losses and reduction of gene families were found and could be related to the extreme parasitic lifestyle of these species. The application of cutting-edge sequencing technology allowed the characterization of the terminal sequences of chromosomes that possess unique characteristics. Here, we review the current status of knowledge of complete cestode genomes and place it within a comparative genomics perspective. Multidisciplinary work together with the implementation of new technologies will provide valuable information that can certainly improve our chances to finally eradicate or at least control diseases caused by cestodes.

Eukaryotic initiation factor EIF-3.G augments mRNA translation efficiency to regulate neuronal activity

The translation initiation complex eIF3 imparts specialized functions to regulate protein expression. However, understanding of eIF3 activities in neurons remains limited despite widespread dysregulation of eIF3 subunits in neurological disorders. Here, we report a selective role of the C. elegans RNA-binding subunit EIF-3.G in shaping the neuronal protein landscape. We identify a missense mutation in the conserved Zinc-Finger (ZF) of EIF-3.G that acts in a gain-of-function manner to dampen neuronal hyperexcitation. Using neuron-type-specific seCLIP, we systematically mapped EIF-3.G-mRNA interactions and identified EIF-3.G occupancy on GC-rich 5′UTRs of a select set of mRNAs enriched in activity-dependent functions. We demonstrate that the ZF mutation in EIF-3.G alters translation in a 5′UTR-dependent manner. Our study reveals an in vivo mechanism for eIF3 in governing neuronal protein levels to control neuronal activity states and offers insights into how eIF3 dysregulation contributes to neurological disorders.

Genome-wide profiling reveals functional interplay of DNA sequence composition, transcriptional activity, and nucleosome positioning in driving DNA supercoiling and helix destabilization in C. elegans

DNA topology and alternative DNA structures are implicated in regulating diverse biological processes. Although biomechanical properties of these structures have been studied extensively in vitro, characterization in vivo, particularly in multicellular organisms, is limited. We devised new methods to map DNA supercoiling and single-stranded DNA in Caenorhabditis elegans embryos and diapause larvae. To map supercoiling, we quantified the incorporation of biotinylated psoralen into DNA using high-throughput sequencing. To map single-stranded DNA, we combined permanganate treatment with genome-wide sequencing of induced double-stranded breaks. We found high levels of negative supercoiling at transcription start sites (TSSs) in embryos. GC-rich regions flanked by a sharp GC-to-AT transition delineate boundaries of supercoil propagation. In contrast to TSSs in embryos, TSSs in diapause larvae showed dramatic reductions in negative supercoiling without concomitant attenuation of transcription, suggesting developmental-stage-specific regulation. To assess whether alternative DNA structures control chromosome architecture and gene expression, we examined DNA supercoiling in the context of X-Chromosome dosage compensation. We showed that the condensin dosage compensation complex creates negative supercoils locally at its highest-occupancy binding sites but found no evidence for large-scale supercoiling domains along X Chromosomes. In contrast to transcription-coupled negative supercoiling, single-strandedness, which is most pronounced at transcript end sites, is dependent on high AT content and symmetrically positioned nucleosomes. We propose that sharp transitions in sequence composition at functional genomic elements constitute a common regulatory code and that DNA structure and propagation of torsional stress at regulatory elements are critical parameters in shaping important developmental events.

Genomics of the Parasitic Nematode Ascaris and Its Relatives

Nematodes of the genus Ascaris are important parasites of humans and swine, and the phylogenetically related genera (Parascaris, Toxocara, and Baylisascaris) infect mammals of veterinary interest. Over the last decade, considerable genomic resources have been established for Ascaris, including complete germline and somatic genomes, comprehensive mRNA and small RNA transcriptomes, as well as genome-wide histone and chromatin data. These datasets provide a major resource for studies on the basic biology of these parasites and the host-parasite relationship. Ascaris and its relatives undergo programmed DNA elimination, a highly regulated process where chromosomes are fragmented and portions of the genome are lost in embryonic cells destined to adopt a somatic fate, whereas the genome remains intact in germ cells. Unlike many model organisms, Ascaris transcription drives early development beginning prior to pronuclear fusion. Studies on Ascaris demonstrated a complex small RNA network even in the absence of a piRNA pathway. Comparative genomics of these ascarids has provided perspectives on nematode sex chromosome evolution, programmed DNA elimination, and host-parasite coevolution. The genomic resources enable comparison of proteins across diverse species, revealing many new potential drug targets that could be used to control these parasitic nematodes.

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.

Improving Transgenesis Efficiency and CRISPR-Associated Tools Through Codon Optimization and Native Intron Addition in Pristionchus Nematodes

A lack of appropriate molecular tools is one obstacle that prevents in-depth mechanistic studies in many organisms. Transgenesis, clustered regularly interspaced short palindromic repeats (CRISPR)-associated engineering, and related tools are fundamental in the modern life sciences, but their applications are still limited to a few model organisms. In the phylum Nematoda, transgenesis can only be performed in a handful of species other than Caenorhabditis elegans, and additionally, other species suffer from significantly lower transgenesis efficiencies. We hypothesized that this may in part be due to incompatibilities of transgenes in the recipient organisms. Therefore, we investigated the genomic features of 10 nematode species from three of the major clades representing all different lifestyles. We found that these species show drastically different codon usage bias and intron composition. With these findings, we used the species Pristionchus pacificus as a proof of concept for codon optimization and native intron addition. Indeed, we were able to significantly improve transgenesis efficiency, a principle that may be usable in other nematode species. In addition, with the improved transgenes, we developed a fluorescent co-injection marker in P. pacificus for the detection of CRISPR-edited individuals, which helps considerably to reduce associated time and costs.

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.

Genomic and transcriptomic variation defines the chromosome-scale assembly of Haemonchus contortus, a model gastrointestinal worm

Haemonchus contortus is a globally distributed and economically important gastrointestinal pathogen of small ruminants and has become a key nematode model for studying anthelmintic resistance and other parasite-specific traits among a wider group of parasites including major human pathogens. Here, we report using PacBio long-read and OpGen and 10X Genomics long-molecule methods to generate a highly contiguous 283.4 Mbp chromosome-scale genome assembly including a resolved sex chromosome for the MHco3(ISE).N1 isolate. We show a remarkable pattern of conservation of chromosome content with Caenorhabditis elegans, but almost no conservation of gene order. Short and long-read transcriptome sequencing allowed us to define coordinated transcriptional regulation throughout the parasite's life cycle and refine our understanding of cis- and trans-splicing. Finally, we provide a comprehensive picture of chromosome-wide genetic diversity both within a single isolate and globally. These data provide a high-quality comparison for understanding the evolution and genomics of Caenorhabditis and other nematodes and extend the experimental tractability of this model parasitic nematode in understanding helminth biology, drug discovery and vaccine development, as well as important adaptive traits such as drug resistance.

Regulation of ribosomal protein genes: An ordered anarchy

Ribosomal protein genes are among the most highly expressed genes in most cell types. Their products are generally essential for ribosome synthesis, which is the cornerstone for cell growth and proliferation. Many cellular resources are dedicated to producing ribosomal proteins and thus this process needs to be regulated in ways that carefully balance the supply of nascent ribosomal proteins with the demand for new ribosomes. Ribosomal protein genes have classically been viewed as a uniform interconnected regulon regulated in eukaryotic cells by target of rapamycin and protein kinase A pathway in response to changes in growth conditions and/or cellular status. However, recent literature depicts a more complex picture in which the amount of ribosomal proteins produced varies between genes in response to two overlapping regulatory circuits. The first includes the classical general ribosome‐producing program and the second is a gene‐specific feature responsible for fine‐tuning the amount of ribosomal proteins produced from each individual ribosomal gene. Unlike the general pathway that is mainly controlled at the level of transcription and translation, this specific regulation of ribosomal protein genes is largely achieved through changes in pre‐mRNA splicing efficiency and mRNA stability. By combining general and specific regulation, the cell can coordinate ribosome production, while allowing functional specialization and diversity. Here we review the many ways ribosomal protein genes are regulated, with special focus on the emerging role of posttranscriptional regulatory events in fine‐tuning the expression of ribosomal protein genes and its role in controlling the potential variation in ribosome functions.

This article is categorized under:

Translation > Ribosome Biogenesis

Translation > Ribosome Structure/Function

Translation > Translation Regulation

Tissue-Specific Transcription Footprinting Using RNA PoI DamID (RAPID) in Caenorhabditis elegans

Differential gene expression across cell types underlies development and cell physiology in multicellular organisms. Caenorhabditis elegans is a powerful, extensively used model to address these biological questions. A remaining bottleneck relates to the difficulty to obtain comprehensive tissue-specific gene transcription data, since available methods are still challenging to execute and/or require large worm populations. Here, we introduce the RNA Polymerase DamID (RAPID) approach, in which the Dam methyltransferase is fused to a ubiquitous RNA polymerase subunit to create transcriptional footprints via methyl marks on the DNA of transcribed genes. To validate the method, we determined the polymerase footprints in whole animals, in sorted embryonic blastomeres and in different tissues from intact young adults by driving tissue-specific Dam fusion expression. We obtained meaningful transcriptional footprints in line with RNA-sequencing (RNA-seq) studies in whole animals or specific tissues. To challenge the sensitivity of RAPID and demonstrate its utility to determine novel tissue-specific transcriptional profiles, we determined the transcriptional footprints of the pair of XXX neuroendocrine cells, representing 0.2% of the somatic cell content of the animals. We identified 3901 candidate genes with putatively active transcription in XXX cells, including the few previously known markers for these cells. Using transcriptional reporters for a subset of new hits, we confirmed that the majority of them were expressed in XXX cells and identified novel XXX-specific markers. Taken together, our work establishes RAPID as a valid method for the determination of RNA polymerase footprints in specific tissues of C. elegans without the need for cell sorting or RNA tagging.

RNA structures in alternative splicing and back-splicing

Alternative splicing greatly expands the transcriptomic and proteomic diversities related to physiological and developmental processes in higher eukaryotes. Splicing of long noncoding RNAs, and back- and trans- splicing further expanded the regulatory repertoire of alternative splicing. RNA structures were shown to play an important role in regulating alternative splicing and back-splicing. Application of novel sequencing technologies made it possible to identify genome-wide RNA structures and interaction networks, which might provide new insights into RNA splicing regulation in vitro to in vivo. The emerging transcription-folding-splicing paradigm is changing our understanding of RNA alternative splicing regulation. Here, we review the insights into the roles and mechanisms of RNA structures in alternative splicing and back-splicing, as well as how disruption of these structures affects alternative splicing and then leads to human diseases. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Structure and Dynamics > Influence of RNA Structure in Biological Systems.

Systematic analysis of long intergenic non-coding RNAs in C. elegans germline uncovers roles in somatic growth

Long intergenic non-coding RNAs (lincRNAs) are transcripts longer than 200 nucleotides that are transcribed from non-coding loci yet undergo biosynthesis similar to coding mRNAs. The disproportional number of lincRNAs expressed in testes suggests that lincRNAs are important during gametogenesis, but experimental evidence has implicated very few lincRNAs in this process. We took advantage of the relatively limited number of lincRNAs in the genome of the nematode Caenorhabditis elegans to systematically analyse the functions of lincRNAs during meiosis. We deleted six lincRNA genes that are highly and dynamically expressed in the C. elegans gonad and tested the effects on central meiotic processes. Surprisingly, whereas the lincRNA deletions did not strongly impact fertility, germline apoptosis, crossovers, or synapsis, linc-4 was required for somatic growth. Slower growth was observed in linc-4-deletion mutants and in worms depleted of linc-4 using RNAi, indicating that linc-4 transcripts are required for this post-embryonic process. Unexpectedly, analysis of worms depleted of linc-4 in soma versus germline showed that the somatic role stems from linc-4 expression in germline cells. This unique feature suggests that some lincRNAs, like some small non-coding RNAs, are required for germ-soma interactions.

Epimutations driven by small RNAs arise frequently but most have limited duration in Caenorhabditis elegans

Epigenetic regulation involves changes in gene expression independent of DNA sequence variation that are inherited through cell division. In addition to a fundamental role in cell differentiation, some epigenetic changes can also be transmitted transgenerationally through meiosis. Epigenetic alterations (epimutations) could thus contribute to heritable variation within populations and be subject to evolutionary processes such as natural selection and drift. However, the rate at which epimutations arise and their typical persistence are unknown, making it difficult to evaluate their potential for evolutionary adaptation. Here, we perform a genome-wide study of epimutations in a metazoan organism. We use experimental evolution to characterize the rate, spectrum and stability of epimutations driven by small silencing RNAs in the model nematode Caenorhabditis elegans. We show that epimutations arise spontaneously at a rate approximately 25 times greater than DNA sequence changes and typically have short half-lives of two to three generations. Nevertheless, some epimutations last at least ten generations. Epimutations mediated by small RNAs may thus contribute to evolutionary processes over a short timescale but are unlikely to bring about long-term divergence in the absence of selection.

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.

Evolutionary Dynamics of the SKN-1 → MED → END-1,3 Regulatory Gene Cascade in Caenorhabditis Endoderm Specification

Gene regulatory networks and their evolution are important in the study of animal development. In the nematode, Caenorhabditis elegans, the endoderm (gut) is generated from a single embryonic precursor, E. Gut is specified by the maternal factor SKN-1, which activates the MED → END-1,3 → ELT-2,7 cascade of GATA transcription factors. In this work, genome sequences from over two dozen species within the Caenorhabditis genus are used to identify MED and END-1,3 orthologs. Predictions are validated by comparison of gene structure, protein conservation, and putative cis-regulatory sites. All three factors occur together, but only within the Elegans supergroup, suggesting they originated at its base. The MED factors are the most diverse and exhibit an unexpectedly extensive gene amplification. In contrast, the highly conserved END-1 orthologs are unique in nearly all species and share extended regions of conservation. The END-1,3 proteins share a region upstream of their zinc finger and an unusual amino-terminal poly-serine domain exhibiting high codon bias. Compared with END-1, the END-3 proteins are otherwise less conserved as a group and are typically found as paralogous duplicates. Hence, all three factors are under different evolutionary constraints. Promoter comparisons identify motifs that suggest the SKN-1, MED, and END factors function in a similar gut specification network across the Elegans supergroup that has been conserved for tens of millions of years. A model is proposed to account for the rapid origin of this essential kernel in the gut specification network, by the upstream intercalation of duplicate genes into a simpler ancestral network.

Trans-splicing of mRNAs links gene transcription to translational control regulated by mTOR

BACKGROUND

In phylogenetically diverse organisms, the 5' ends of a subset of mRNAs are trans-spliced with a spliced leader (SL) RNA. The functions of SL trans-splicing, however, remain largely enigmatic.

RESULTS

We quantified translation genome-wide in the marine chordate, Oikopleura dioica, under inhibition of mTOR, a central growth regulator. Translation of trans-spliced TOP mRNAs was suppressed, consistent with a role of the SL sequence in nutrient-dependent translational control of growth-related mRNAs. Under crowded, nutrient-limiting conditions, O. dioica continued to filter-feed, but arrested growth until favorable conditions returned. Upon release from unfavorable conditions, initial recovery was independent of nutrient-responsive, trans-spliced genes, suggesting animal density sensing as a first trigger for resumption of development.

CONCLUSION

Our results are consistent with a proposed role of trans-splicing in the coordinated translational down-regulation of nutrient-responsive genes under growth-limiting conditions.

Post-transcriptional Regulation by 3' UTRs Can Be Masked by Regulatory Elements in 5' UTRs

In mRNA sequences, 3' UTRs are thought to contain most elements that specifically regulate localization, turnover, and translation. Although high-throughput experiments indicate that many RNA-binding proteins (RBPs) also bind 5' UTRs, much less is known about specific post-transcriptional control exerted by 5' UTRs. GLD-1 is a conserved RBP and a translational repressor with essential roles in Caenorhabditis elegans germ cell development. Previously, we showed that GLD-1 binds highly conserved sites in both 3' and 5' UTRs. Here, by targeted single-copy insertion of transgenes, we systematically tested in vivo functionality of 5' and 3' UTR binding sites individually and in combination. Our data show that sites in 5' UTRs mediate specific and strong translational repression, independent of exact position. Intriguingly, we found that the functionality of 3' UTR sites can be masked by 5' UTR sites and vice versa. We conclude that it is important to study both UTRs simultaneously.

FLAM-seq: full-length mRNA sequencing reveals principles of poly(A) tail length control

Although messenger RNAs are key molecules for understanding life, until now, no method has existed to determine the full-length sequence of endogenous mRNAs including their poly(A) tails. Moreover, although non-A nucleotides can be incorporated in poly(A) tails, there also exists no method to accurately sequence them. Here, we present full-length poly(A) and mRNA sequencing (FLAM-seq), a rapid and simple method for high-quality sequencing of entire mRNAs. We report a complementary DNA library preparation method coupled to single-molecule sequencing to perform FLAM-seq. Using human cell lines, brain organoids and Caenorhabditis elegans we show that FLAM-seq delivers high-quality full-length mRNA sequences for thousands of different genes per sample. We find that 3' untranslated region length is correlated with poly(A) tail length, that alternative polyadenylation sites and alternative promoters for the same gene are linked to different tail lengths, and that tails contain a substantial number of cytosines.

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.

FLAM-seq: full-length mRNA sequencing reveals principles of poly(A) tail length control

Although messenger RNAs are key molecules for understanding life, until now, no method has existed to determine the full-length sequence of endogenous mRNAs including their poly(A) tails. Moreover, although non-A nucleotides can be incorporated in poly(A) tails, there also exists no method to accurately sequence them. Here, we present full-length poly(A) and mRNA sequencing (FLAM-seq), a rapid and simple method for high-quality sequencing of entire mRNAs. We report a complementary DNA library preparation method coupled to single-molecule sequencing to perform FLAM-seq. Using human cell lines, brain organoids and Caenorhabditis elegans we show that FLAM-seq delivers high-quality full-length mRNA sequences for thousands of different genes per sample. We find that 3' untranslated region length is correlated with poly(A) tail length, that alternative polyadenylation sites and alternative promoters for the same gene are linked to different tail lengths, and that tails contain a substantial number of cytosines.

The Caenorhabditis elegans Transgenic Toolbox

The power of any genetic model organism is derived, in part, from the ease with which gene expression can be manipulated. The short generation time and invariant developmental lineage have made Caenorhabditis elegans very useful for understanding, e.g., developmental programs, basic cell biology, neurobiology, and aging. Over the last decade, the C. elegans transgenic toolbox has expanded considerably, with the addition of a variety of methods to control expression and modify genes with unprecedented resolution. Here, we provide a comprehensive overview of transgenic methods in C. elegans, with an emphasis on recent advances in transposon-mediated transgenesis, CRISPR/Cas9 gene editing, conditional gene and protein inactivation, and bipartite systems for temporal and spatial control of expression.

Dietary restriction induces posttranscriptional regulation of longevity genes

Dietary restriction (DR) increases life span through adaptive changes in gene expression. To understand more about these changes, we analyzed the transcriptome and translatome of Caenorhabditis elegans subjected to DR. Transcription of muscle regulatory and structural genes increased, whereas increased expression of amino acid metabolism and neuropeptide signaling genes was controlled at the level of translation. Evaluation of posttranscriptional regulation identified putative roles for RNA-binding proteins, RNA editing, miRNA, alternative splicing, and nonsense-mediated decay in response to nutrient limitation. Using RNA interference, we discovered several differentially expressed genes that regulate life span. We also found a compensatory role for translational regulation, which offsets dampened expression of a large subset of transcriptionally down-regulated genes. Furthermore, 3' UTR editing and intron retention increase under DR and correlate with diminished translation, whereas trans-spliced genes are refractory to reduced translation efficiency compared with messages with the native 5' UTR. Finally, we find that smg-6 and smg-7, which are genes governing selection and turnover of nonsense-mediated decay targets, are required for increased life span under DR.

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.

New insights from Opisthorchis felineus genome: update on genomics of the epidemiologically important liver flukes

Background

The three epidemiologically important Opisthorchiidae liver flukes Opisthorchis felineus, O. viverrini, and Clonorchis sinensis, are believed to harbour similar potencies to provoke hepatobiliary diseases in their definitive hosts, although their populations have substantially different ecogeographical aspects including habitat, preferred hosts, population structure. Lack of O. felineus genomic data is an obstacle to the development of comparative molecular biological approaches necessary to obtain new knowledge about the biology of Opisthorchiidae trematodes, to identify essential pathways linked to parasite-host interaction, to predict genes that contribute to liver fluke pathogenesis and for the effective prevention and control of the disease.

Results

Here we present the first draft genome assembly of O. felineus and its gene repertoire accompanied by a comparative analysis with that of O. viverrini and Clonorchis sinensis. We observed both noticeably high heterozygosity of the sequenced individual and substantial genetic diversity in a pooled sample. This indicates that potency of O. felineus population for rapid adaptive response to control and preventive measures of opisthorchiasis is higher than in O. viverrini and C. sinensis. We also have found that all three species are characterized by more intensive involvement of trans-splicing in RNA processing compared to other trematodes.

Conclusion

All revealed peculiarities of structural organization of genomes are of extreme importance for a proper description of genes and their products in these parasitic species. This should be taken into account both in academic and applied research of epidemiologically important liver flukes. Further comparative genomics studies of liver flukes and non-carcinogenic flatworms allow for generation of well-grounded hypotheses on the mechanisms underlying development of cholangiocarcinoma associated with opisthorchiasis and clonorchiasis as well as species-specific mechanisms of these diseases.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5752-8) contains supplementary material, which is available to authorized users.

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.

Chromatin accessibility dynamics across C. elegans development and ageing

An essential step for understanding the transcriptional circuits that control development and physiology is the global identification and characterization of regulatory elements. Here, we present the first map of regulatory elements across the development and ageing of an animal, identifying 42,245 elements accessible in at least one Caenorhabditis elegans stage. Based on nuclear transcription profiles, we define 15,714 protein-coding promoters and 19,231 putative enhancers, and find that both types of element can drive orientation-independent transcription. Additionally, more than 1000 promoters produce transcripts antisense to protein coding genes, suggesting involvement in a widespread regulatory mechanism. We find that the accessibility of most elements changes during development and/or ageing and that patterns of accessibility change are linked to specific developmental or physiological processes. The map and characterization of regulatory elements across C. elegans life provides a platform for understanding how transcription controls development and ageing.

Caenorhabditis elegans as an emerging model system in environmental epigenetics

The roundworm Caenorhabitis elegans has been an established model organism for the study of genetics and developmental biology, including studies of transcriptional regulation, since the 1970s. This model organism has continued to be used as a classical model system as the field of transcriptional regulation has expanded to include scientific advances in epigenetics and chromatin biology. In the last several decades, C. elegans has emerged as a powerful model for environmental toxicology, particularly for the study of chemical genotoxicity. Here, we outline the utility and applicability of C. elegans as a powerful model organism for mechanistic studies of environmental influences on the epigenome. Our goal in this article is to inform the field of environmental epigenetics of the strengths and limitations of the well-established C. elegans model organism as an emerging model for medium-throughput, in vivo exploration of the role of exogenous chemical stimuli in transcriptional regulation, developmental epigenetic reprogramming, and epigenetic memory and inheritance. As the field of environmental epigenetics matures, and research begins to map mechanisms underlying observed associations, new toolkits and model systems, particularly manipulable, scalable in vivo systems that accurately model human transcriptional regulatory circuits, will provide an essential experimental bridge between in vitro biochemical experiments and mammalian model systems. Environ. Mol. Mutagen. 59:560-575, 2018. © 2018 Wiley Periodicals, Inc.

Transcriptome analysis of adult Caenorhabditis elegans cells reveals tissue-specific gene and isoform expression

The biology and behavior of adults differ substantially from those of developing animals, and cell-specific information is critical for deciphering the biology of multicellular animals. Thus, adult tissue-specific transcriptomic data are critical for understanding molecular mechanisms that control their phenotypes. We used adult cell-specific isolation to identify the transcriptomes of C. elegans' four major tissues (or "tissue-ome"), identifying ubiquitously expressed and tissue-specific "enriched" genes. These data newly reveal the hypodermis' metabolic character, suggest potential worm-human tissue orthologies, and identify tissue-specific changes in the Insulin/IGF-1 signaling pathway. Tissue-specific alternative splicing analysis identified a large set of collagen isoforms. Finally, we developed a machine learning-based prediction tool for 76 sub-tissue cell types, which we used to predict cellular expression differences in IIS/FOXO signaling, stage-specific TGF-β activity, and basal vs. memory-induced CREB transcription. Together, these data provide a rich resource for understanding the biology governing multicellular adult animals.

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.

Endosome maturation factors Rabenosyn-5/VPS45 and caveolin-1 regulate ciliary membrane and polycystin-2 homeostasis

Primary cilium structure and function relies on control of ciliary membrane homeostasis, regulated by membrane trafficking processes that deliver and retrieve ciliary components at the periciliary membrane. However, the molecular mechanisms controlling ciliary membrane establishment and maintenance, especially in relation to endocytosis, remain poorly understood. Here, using Caenorhabditis elegans, we describe closely linked functions for early endosome (EE) maturation factors RABS-5 (Rabenosyn-5) and VPS-45 (VPS45) in regulating cilium length and morphology, ciliary and periciliary membrane volume, and ciliary signalling-related sensory behaviour. We demonstrate that RABS-5 and VPS-45 control periciliary vesicle number and levels of select EE/endocytic markers (WDFY-2, CAV-1) and the ciliopathy membrane receptor PKD-2 (polycystin-2). Moreover, we show that CAV-1 (caveolin-1) also controls PKD-2 ciliary levels and associated sensory behaviour. These data link RABS-5 and VPS-45 ciliary functions to the processing of periciliary-derived endocytic vesicles and regulation of ciliary membrane homeostasis. Our findings also provide insight into the regulation of PKD-2 ciliary levels via integrated endosomal sorting and CAV-1-mediated endocytosis.

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.

An Expanded Role for the RFX Transcription Factor DAF-19, with Dual Functions in Ciliated and Nonciliated Neurons

Regulatory Factor X (RFX) transcription factors (TFs) are best known for activating genes required for ciliogenesis in both vertebrates and invertebrates. In humans, eight RFX TFs have a variety of tissue-specific functions, while in the worm Caenorhabditis elegans, the sole RFX gene, daf-19, encodes a set of nested isoforms. Null alleles of daf-19 confer pleiotropic effects including altered development with a dauer constitutive phenotype, complete absence of cilia and ciliary proteins, and defects in synaptic protein maintenance. We sought to identify RFX/daf-19 target genes associated with neuronal functions other than ciliogenesis using comparative transcriptome analyses at different life stages of the worm. Subsequent characterization of gene expression patterns revealed one set of genes activated in the presence of DAF-19 in ciliated sensory neurons, whose activation requires the daf-19c isoform, also required for ciliogenesis. A second set of genes is downregulated in the presence of DAF-19, primarily in nonsensory neurons. The human orthologs of some of these neuronal genes are associated with human diseases. We report the novel finding that daf-19a is directly or indirectly responsible for downregulation of these neuronal genes in C. elegans by characterizing a new mutation affecting the daf-19a isoform (tm5562) and not associated with ciliogenesis, but which confers synaptic and behavioral defects. Thus, we have identified a new regulatory role for RFX TFs in the nervous system. The new daf-19 candidate target genes we have identified by transcriptomics will serve to uncover the molecular underpinnings of the pleiotropic effects that daf-19 exerts on nervous system function.