Metazoan operons accelerate recovery from growth-arrested states

Selective RNA Processing and Stabilization are Multi-Layer and Stoichiometric Regulators of Gene Expression in Escherichia coli

Selective RNA processing and stabilization (SRPS) facilitates the differential expression of multiple genes in polycistronic operons. However, how the coordinated actions of SRPS-related enzymes affect stoichiometric regulation remains unclear. In the present study, the first genome-wide targetome analysis is reported of these enzymes in Escherichia coli, at a single-nucleotide resolution. A strictly linear relationship is observed between the RNA pyrophosphohydrolase processing ratio and scores assigned to the first three nucleotides of the primary transcript. Stem-loops associated with PNPase targetomes exhibit a folding free energy that is negatively correlated with the termination ratio of PNPase at the 3' end. More than one-tenth of the RNase E processing sites in the 5'-untranslated regions（UTR） form different stem-loops that affect ribosome-binding and translation efficiency. The effectiveness of the SRPS elements is validated using a dual-fluorescence reporter system. The findings highlight a multi-layer and quantitative regulatory method for optimizing the stoichiometric expression of genes in bacteria and promoting the application of SRPS in synthetic biology.

Citrobacter sp. Y3 harbouring novel gene HBCD-hd-1 mineralizes hexabromocyclododecane via new metabolic pathways according to multi-omics characterization

Hexabromocyclododecane (HBCD) is a typical persistent organic pollutant that is widely detected in the environment. Despite the significant efforts put into its mineralisation, there is still a lack of microorganism resources that can completely mineralise HBCD. Stable isotope analysis revealed that the Citrobacter sp. Y3 can use [¹³C]HBCD as its sole carbon source and degrade or even mineralise it into ¹³CO₂, with a maximum conversion rate of 100% in approximately 14 days. Strain Y3 could completely mineralise HBCD, which it used as its only carbon source, and six debromination enzymes related to HBCD degradation were found in Y3, including haloalkane dehalogenase (DhaA), haloacid dehalogenase (HAD), etc. A functional gene named HBCD-hd-1, encoding a HAD, was found to be upregulated during HBCD degradation and heterologously expressed in Escherichia coli. Recombinant E. coli with the HBCD-hd-1 gene transformed the typical intermediate 4-bromobutyric acid to 4-hydroxybutanoic acid and showed excellent degradation performance on HBCD, accompanied by nearly 100% bromine (Br) ion generation. The expression of HBCD-hd-1 in Y3 rapidly accelerated the biodegradation of HBCD. With HBCD as its sole carbon source, strain Y3 could potentially degrade HBCD, especially in a low-nutrient environment.

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.

Neighbours matter: Effects of genomic organization on gene expression plasticity in response to environmental stresses during biological invasions

Gene expression regulation has been widely recognized as an important molecular mechanism underlying phenotypic plasticity in environmental adaptation. However, it remains largely unexplored on the effects of genomic organization on gene expression plasticity under environmental stresses during biological invasions. Here, we use an invasive model ascidian, Ciona robusta, to investigate how genomic organization affects gene expression in response to salinity stresses during range expansions. Our study showed that neighboring genes were co-expressed and approximately 30% of stress responsive genes were physically clustered on chromosomes. Such coordinated expression was substantially affected by the physical distance and orientation of genes. Interestingly, the overall expression correlation of neighboring genes was significantly decreased under high salinity stresses, illustrating that the co-expression regulation could be disrupted by salinity challenges. Furthermore, the clustering of genes was associated with their function constraints and expression patterns - operon genes enriched in gene expression machinery had the highest transcriptional activity and expression stability. Notably, our analyses showed that the tail-to-tail genes, mainly involved in biological functions related to phosphorylation, homeostatic process, and ion transport, exhibited higher intrinsic expression variability and greater response to salinity challenges. Altogether, the results obtained here provide new insights into the effects of gene organization on gene expression plasticity under environmental challenges, hence improving our knowledge on mechanisms of rapid environmental adaptation during biological invasions.

Steroid hormone pathways coordinate developmental diapause and olfactory remodeling in Pristionchus pacificus

Developmental and behavioral plasticity allow animals to prioritize alternative genetic programs during fluctuating environments. Behavioral remodeling may be acute in animals that interact with host organisms, since reproductive adults and the developmentally arrested larvae often have different ethological needs for chemical stimuli. To understand the genes that coordinate the development and host-seeking behavior, we used the entomophilic nematode Pristionchus pacificus to characterize dauer-constitutive mutants (Daf-c) that inappropriately enter developmental diapause to become dauer larvae. We found two Daf-c loci with dauer-constitutive and cuticle exsheathment phenotypes that can be rescued by the feeding of Δ7-dafachronic acid, and that are dependent on the conserved canonical steroid hormone receptor Ppa-DAF-12. Specifically at one locus, deletions in the sole hydroxysteroid dehydrogenase (HSD) in P. pacificus resulted in Daf-c phenotypes. Ppa-hsd-2 is expressed in the canal-associated neurons (CANs) and excretory cells whose homologous cells in Caenorhabditis elegans are not known to be involved in the dauer decision. While in wildtype only dauer larvae are attracted to host odors, hsd-2 mutant adults show enhanced attraction to the host beetle pheromone, along with ectopic activation of a marker for putative olfactory neurons, Ppa-odr-3. Surprisingly, this enhanced odor attraction acts independently of the Δ7-DA/DAF-12 module, suggesting that Ppa-HSD-2 may be responsible for several steroid hormone products involved in coordinating the dauer decision and host-seeking behavior in P. pacificus.

RNA polymerase II CTD S2P is dispensable for embryogenesis but mediates exit from developmental diapause in C. elegans

Serine 2 phosphorylation (S2P) within the CTD of RNA polymerase II is considered a Cdk9/Cdk12-dependent mark required for 3'-end processing. However, the relevance of CTD S2P in metazoan development is unknown. We show that cdk-12 lesions or a full-length CTD S2A substitution results in an identical phenotype in Caenorhabditis elegans Embryogenesis occurs in the complete absence of S2P, but the hatched larvae arrest development, mimicking the diapause induced when hatching occurs in the absence of food. Genome-wide analyses indicate that when CTD S2P is inhibited, only a subset of growth-related genes is not properly expressed. These genes correspond to SL2 trans-spliced mRNAs located in position 2 and over within operons. We show that CDK-12 is required for maximal occupancy of cleavage stimulatory factor necessary for SL2 trans-splicing. We propose that CTD S2P functions as a gene-specific signaling mark ensuring the nutritional control of the C. elegans developmental program.

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.

Application of Transcriptional Gene Modules to Analysis of Caenorhabditis elegans' Gene Expression Data

Identification of co-expressed sets of genes (gene modules) is used widely for grouping functionally related genes during transcriptomic data analysis. An organism-wide atlas of high-quality gene modules would provide a powerful tool for unbiased detection of biological signals from gene expression data. Here, using a method based on independent component analysis we call DEXICA, we have defined and optimized 209 modules that broadly represent transcriptional wiring of the key experimental organism C. elegans. These modules represent responses to changes in the environment (e.g., starvation, exposure to xenobiotics), genes regulated by transcriptions factors (e.g., ATFS-1, DAF-16), genes specific to tissues (e.g., neurons, muscle), genes that change during development, and other complex transcriptional responses to genetic, environmental and temporal perturbations. Interrogation of these modules reveals processes that are activated in long-lived mutants in cases where traditional analyses of differentially expressed genes fail to do so. Additionally, we show that modules can inform the strength of the association between a gene and an annotation (e.g., GO term). Analysis of “module-weighted annotations” improves on several aspects of traditional annotation-enrichment tests and can aid in functional interpretation of poorly annotated genes. We provide an online interactive resource with tutorials at http://genemodules.org/, in which users can find detailed information on each module, check genes for module-weighted annotations, and use both of these to analyze their own gene expression data (generated using any platform) or gene sets of interest.

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.

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.

Genome structure predicts modular transcriptome responses to genetic and environmental conditions

Understanding the plasticity, robustness and modularity of transcriptome expression to genetic and environmental conditions is crucial to deciphering how organisms adapt in nature. To test how genome architecture influences transcriptome profiles, we quantified expression responses for distinct temperature-adapted genotypes of the nematode Caenorhabditis briggsae when exposed to chronic temperature stresses throughout development. We found that 56% of the 8,795 differentially expressed genes show genotype-specific changes in expression in response to temperature (genotype-by-environment interactions, GxE). Most genotype-specific responses occur under heat stress, indicating that cold vs. heat stress responses involve distinct genomic architectures. The 22 co-expression modules that we identified differ in their enrichment of genes with genetic vs. environmental vs. interaction effects, as well as their genomic spatial distributions, functional attributes and rates of molecular evolution at the sequence level. Genes in modules enriched for simple effects of either genotype or temperature alone tend to evolve especially rapidly, consistent with disproportionate influence of adaptation or weaker constraint on these subsets of loci. Chromosome-scale heterogeneity in nucleotide polymorphism, however, rather than the scale of individual genes predominates as the source of genetic differences among expression profiles, and natural selection regimes are largely decoupled between coding sequences and noncoding flanking sequences that contain cis-regulatory elements. These results illustrate how the form of transcriptome modularity and genome structure contribute to predictable profiles of evolutionary change.

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.

Food perception without ingestion leads to metabolic changes and irreversible developmental arrest in C. elegans

BACKGROUND

Developmental physiology is very sensitive to nutrient availability. For instance, in the nematode Caenorhabditis elegans, newly hatched L1-stage larvae require food to initiate postembryonic development. In addition, larvae arrested in the dauer diapause, a non-feeding state of developmental arrest that occurs during the L3 stage, initiate recovery when exposed to food. Despite the essential role of food in C. elegans development, the contribution of food perception versus ingestion on physiology has not been delineated.

RESULTS

We used a pharmacological approach to uncouple the effects of food (bacteria) perception and ingestion in C. elegans. Perception was not sufficient to promote postembryonic development in L1-stage larvae. However, L1 larvae exposed to food without ingestion failed to develop upon return to normal culture conditions, instead displaying an irreversible arrest phenotype. Inhibition of gene expression during perception rescued subsequent development, demonstrating that the response to perception without feeding is deleterious. Perception altered DAF-16/FOXO subcellular localization, reflecting activation of insulin/IGF signaling (IIS). The insulin-like peptide daf-28 was specifically required, suggesting perception in chemosensory neurons, where it is expressed, regulates peptide synthesis and possibly secretion. However, genetic manipulation of IIS did not modify the irreversible arrest phenotype caused by food perception, revealing that wild-type function of the IIS pathway is not required to produce this phenotype and that other pathways affected by perception of food in the absence of its ingestion are likely to be involved. Gene expression and Nile red staining showed that food perception could alter lipid metabolism and storage. We found that starved larvae sense environmental polypeptides, with similar molecular and developmental effects as perception of bacteria. Environmental polypeptides also promoted recovery from dauer diapause, suggesting that perception of polypeptides plays an important role in the life history of free-living nematodes.

CONCLUSIONS

We conclude that actual ingestion of food is required to initiate postembryonic development in C. elegans. We also conclude that polypeptides are perceived as a food-associated cue in this and likely other animals, initiating a signaling and gene regulatory cascade that alters metabolism in anticipation of feeding and development, but that this response is detrimental if feeding does not occur.

Trans-spliced long non-coding RNA: an emerging regulator of pluripotency

With dual capacities for unlimited self-renewal and pluripotent differentiation, pluripotent stem cells (PSCs) give rise to many cell types in our body and PSC culture systems provide an unparalleled opportunity to study early human development and disease. Accumulating evidence indicates that the molecular mechanisms underlying pluripotency maintenance in PSCs involve many factors. Among these regulators, recent studies have shown that long non-coding RNAs (lncRNAs) can affect the pluripotency circuitry by cooperating with master pluripotency-associated factors. Additionally, trans-spliced RNAs, which are generated by combining two or more pre-mRNA transcripts to produce a chimeric RNA, have been identified as regulators of various biological processes, including human pluripotency. In this review, we summarize and discuss current knowledge about the roles of lncRNAs, including trans-spliced lncRNAs, in controlling pluripotency.

Genetic Interaction Network as an Important Determinant of Gene Order in Genome Evolution

Although it is generally accepted that eukaryotic gene order is not random, the basic principles of gene arrangement on a chromosome remain poorly understood. Here, we extended existing population genetics theories that were based on two-locus models and proposed a hypothesis that genetic interaction networks drive the evolution of eukaryotic gene order. We predicted that genes with positive epistasis would move toward each other in evolution, during which a negative correlation between epistasis and gene distance formed. We tested and confirmed our prediction with computational simulations and empirical data analyses. Importantly, we demonstrated that gene order in the budding yeast could be successfully predicted from the genetic interaction network. Taken together, our study reveals the role of the genetic interaction network in the evolution of gene order, extends our understanding of the encoding principles in genomes, and potentially offers new strategies to improve synthetic biology.

Trans-splicing enhances translational efficiency in C. elegans

Translational efficiency is subject to extensive regulation. However, the factors influencing such regulation are poorly understood. In Caenorhabditis elegans, 62% of genes are trans-spliced to a specific spliced leader (SL1), which replaces part of the native 5' untranslated region (5' UTR). Given the pivotal role the 5' UTR plays in the regulation of translational efficiency, we hypothesized that SL1 trans-splicing functions to regulate translational efficiency. With genome-wide analysis on Ribo-seq data, polysome profiling experiments, and CRISPR-Cas9-based genetic manipulation of trans-splicing sites, we found four lines of evidence in support of this hypothesis. First, SL1 trans-spliced genes have higher translational efficiencies than non-trans-spliced genes. Second, SL1 trans-spliced genes have higher translational efficiencies than non-trans-spliced orthologous genes in other nematode species. Third, an SL1 trans-spliced isoform has higher translational efficiency than the non-trans-spliced isoform of the same gene. Fourth, deletion of trans-splicing sites of endogenous genes leads to reduced translational efficiency. Importantly, we demonstrated that SL1 trans-splicing plays a key role in enhancing translational efficiencies of essential genes. We further discovered that SL1 trans-splicing likely enhances translational efficiency by shortening the native 5' UTRs, hence reducing the presence of upstream start codons (uAUG) and weakening mRNA secondary structures. Taken together, our study elucidates the global function of trans-splicing in enhancing translational efficiency in nematodes, paving the way for further understanding the genomic mechanisms of translational regulation.

Activated entomopathogenic nematode infective juveniles release lethal venom proteins

Entomopathogenic nematodes (EPNs) are unique parasites due to their symbiosis with entomopathogenic bacteria and their ability to kill insect hosts quickly after infection. It is widely believed that EPNs rely on their bacterial partners for killing hosts. Here we disproved this theory by demonstrating that the in vitro activated infective juveniles (IJs) of Steinernema carpocapsae (a well-studied EPN species) release venom proteins that are lethal to several insects including Drosophila melanogaster. We confirmed that the in vitro activation is a good approximation of the in vivo process by comparing the transcriptomes of individual in vitro and in vivo activated IJs. We further analyzed the transcriptomes of non-activated and activated IJs and revealed a dramatic shift in gene expression during IJ activation. We also analyzed the venom proteome using mass spectrometry. Among the 472 venom proteins, proteases and protease inhibitors are especially abundant, and toxin-related proteins such as Shk domain-containing proteins and fatty acid- and retinol-binding proteins are also detected, which are potential candidates for suppressing the host immune system. Many of the venom proteins have conserved orthologs in vertebrate-parasitic nematodes and are differentially expressed during IJ activation, suggesting conserved functions in nematode parasitism. In summary, our findings strongly support a new model that S. carpocapsae and likely other Steinernema EPNs have a more active role in contributing to the pathogenicity of the nematode-bacterium complex than simply relying on their symbiotic bacteria. Furthermore, we propose that EPNs are a good model system for investigating vertebrate- and human-parasitic nematodes, especially regarding the function of excretory/secretory products.

Steinernema carpocapsae belongs to a special group of insect-parasitic nematodes known as entomopathogenic nematodes (EPNs). These differ from other insect parasites in at least two ways; first they kill their hosts quickly (within 2–3 days), and second they associate with bacteria to facilitate their parasitic lifestyle. The infective stage of these parasites, the infective juvenile (IJ) stage, is the only free-living stage and these IJs are developmentally arrested and only reinitiate development once they are inside a suitable host. Little is known about the early stages of parasitism and how these parasites initiate the parasitic phase of their life cycle and reinitiate development. Here we characterized the changes that occur to the nematodes' physical morphology, gene expression, and the release of protein molecules that accompany the transition from developmentally arrested IJ to active, developing parasite. We showed that contrary to long-held assumptions, the nematodes are not merely transporting pathogenic bacteria but that the nematodes contribute to parasitism by releasing toxic proteins into the host. Many of the S. carpocapsae toxins are also found in species of human-parasitic nematodes, and S. carpocapsae may serve as a valuable model for understanding the specific function of these toxins.

A functional difference between native and horizontally acquired genes in bdelloid rotifers

The form of RNA processing known as SL trans-splicing involves the transfer of a short conserved sequence, the spliced leader (SL), from a noncoding SL RNA to the 5' ends of mRNA molecules. SL trans-splicing occurs in several animal taxa, including bdelloid rotifers (Rotifera, Bdelloidea). One striking feature of these aquatic microinvertebrates is the large proportion of foreign genes, i.e. those acquired by horizontal gene transfer from other organisms, in their genomes. However, whether such foreign genes behave similarly to native genes has not been tested in bdelloids or any other animal. We therefore used a combination of experimental and computational methods to examine whether transcripts of foreign genes in bdelloids were SL trans-spliced, like their native counterparts. We found that many foreign transcripts contain SLs, use similar splice acceptor sequences to native genes, and are able to undergo alternative trans-splicing. However, a significantly lower proportion of foreign mRNAs contains SL sequences than native transcripts. This demonstrates a novel functional difference between foreign and native genes in bdelloids and suggests that SL trans-splicing is not essential for the expression of foreign genes, but is acquired during their domestication.

Evolutionary Insights into RNA trans-Splicing in Vertebrates

Pre-RNA splicing is an essential step in generating mature mRNA. RNA trans-splicing combines two separate pre-mRNA molecules to form a chimeric non-co-linear RNA, which may exert a function distinct from its original molecules. Trans-spliced RNAs may encode novel proteins or serve as noncoding or regulatory RNAs. These novel RNAs not only increase the complexity of the proteome but also provide new regulatory mechanisms for gene expression. An increasing amount of evidence indicates that trans-splicing occurs frequently in both physiological and pathological processes. In addition, mRNA reprogramming based on trans-splicing has been successfully applied in RNA-based therapies for human genetic diseases. Nevertheless, clarifying the extent and evolution of trans-splicing in vertebrates and developing detection methods for trans-splicing remain challenging. In this review, we summarize previous research, highlight recent advances in trans-splicing, and discuss possible splicing mechanisms and functions from an evolutionary viewpoint.

Trans-splicing in metazoans: A link to translational control?

The trans-splicing of a spliced-leader RNA to a subset of mRNAs is a phenomenon that occurs in many species, including Caenorhabditis elegans, and yet the driving force for its evolution in disparate groups of animals remains unclear. Polycistronic mRNA resulting from the transcription of operons is resolved via trans-splicing, but operons comprise only a sub-set of trans-spliced genes. Using the marine chordate, Oikopleura dioica, we recently tested the hypothesis that metazoan operons accelerate recovery from growth arrest. We found no supporting evidence for this in O. dioica. Instead we found a striking relationship between trans-splicing and maternal mRNA in O. dioica, C. elegans and the ascidian, Ciona intestinalis. Furthermore, in O. dioica and C. elegans, we found evidence to suggest a role for mTOR signaling in the translational control of growth-related, trans-spliced maternal mRNAs. We propose that this may be a mechanism for adjusting egg number in response to nutrient levels in these species.

Hierarchical sparse coding in the sensory system of Caenorhabditis elegans

Animals with compact sensory systems face an encoding problem where a small number of sensory neurons are required to encode information about its surrounding complex environment. Using Caenorhabditis elegans worms as a model, we ask how chemical stimuli are encoded by a small and highly connected sensory system. We first generated a comprehensive library of transgenic worms where each animal expresses a genetically encoded calcium indicator in individual sensory neurons. This library includes the vast majority of the sensory system in C. elegans. Imaging from individual sensory neurons while subjecting the worms to various stimuli allowed us to compile a comprehensive functional map of the sensory system at single neuron resolution. The functional map reveals that despite the dense wiring, chemosensory neurons represent the environment using sparse codes. Moreover, although anatomically closely connected, chemo- and mechano-sensory neurons are functionally segregated. In addition, the code is hierarchical, where few neurons participate in encoding multiple cues, whereas other sensory neurons are stimulus specific. This encoding strategy may have evolved to mitigate the constraints of a compact sensory system.

Trans-splicing and operons in metazoans: translational control in maternally regulated development and recovery from growth arrest

Polycistronic mRNAs transcribed from operons are resolved via the trans-splicing of a spliced-leader (SL) RNA. Trans-splicing also occurs at monocistronic transcripts. The phlyogenetically sporadic appearance of trans-splicing and operons has made the driving force(s) for their evolution in metazoans unclear. Previous work has proposed that germline expression drives operon organization in Caenorhabditis elegans, and a recent hypothesis proposes that operons provide an evolutionary advantage via the conservation of transcriptional machinery during recovery from growth arrested states. Using a modified cap analysis of gene expression protocol we mapped sites of SL trans-splicing genome-wide in the marine chordate Oikopleura dioica. Tiled microarrays revealed the expression dynamics of trans-spliced genes across development and during recovery from growth arrest. Operons did not facilitate recovery from growth arrest in O. dioica. Instead, we found that trans-spliced transcripts were predominantly maternal. We then analyzed data from C. elegans and Ciona intestinalis and found that an enrichment of trans-splicing and operon gene expression in maternal mRNA is shared between all three species, suggesting that this may be a driving force for operon evolution in metazoans. Furthermore, we found that the majority of known terminal oligopyrimidine (TOP) mRNAs are trans-spliced in O. dioica and that the SL contains a TOP-like motif. This suggests that the SL in O. dioica confers nutrient-dependent translational control to trans-spliced mRNAs via the TOR-signaling pathway. We hypothesize that SL-trans-splicing provides an evolutionary advantage in species that depend on translational control for regulating early embryogenesis, growth and oocyte production in response to nutrient levels.

Genome-wide analysis of trans-splicing in the nematode Pristionchus pacificus unravels conserved gene functions for germline and dauer development in divergent operons

Discovery of trans-splicing in multiple metazoan lineages led to the identification of operon-like gene organization in diverse organisms, including trypanosomes, tunicates, and nematodes, but the functional significance of such operons is not completely understood. To see whether the content or organization of operons serves similar roles across species, we experimentally defined operons in the nematode model Pristionchus pacificus. We performed affinity capture experiments on mRNA pools to specifically enrich for transcripts that are trans-spliced to either the SL1- or SL2-spliced leader, using spliced leader-specific probes. We obtained distinct trans-splicing patterns from the analysis of three mRNA pools (total mRNA, SL1 and SL2 fraction) by RNA-seq. This information was combined with a genome-wide analysis of gene orientation and spacing. We could confirm 2219 operons by RNA-seq data out of 6709 candidate operons, which were predicted by sequence information alone. Our gene order comparison of the Caenorhabditis elegans and P. pacificus genomes shows major changes in operon organization in the two species. Notably, only 128 out of 1288 operons in C. elegans are conserved in P. pacificus. However, analysis of gene-expression profiles identified conserved functions such as an enrichment of germline-expressed genes and higher expression levels of operonic genes during recovery from dauer arrest in both species. These results provide support for the model that a necessity for increased transcriptional efficiency in the context of certain developmental processes could be a selective constraint for operon evolution in metazoans. Our method is generally applicable to other metazoans to see if similar functional constraints regulate gene organization into operons.

Operons are a conserved feature of nematode genomes

The organization of genes into operons, clusters of genes that are co-transcribed to produce polycistronic pre-mRNAs, is a trait found in a wide range of eukaryotic groups, including multiple animal phyla. Operons are present in the class Chromadorea, one of the two main nematode classes, but their distribution in the other class, the Enoplea, is not known. We have surveyed the genomes of Trichinella spiralis, Trichuris muris, and Romanomermis culicivorax and identified the first putative operons in members of the Enoplea. Consistent with the mechanism of polycistronic RNA resolution in other nematodes, the mRNAs produced by genes downstream of the first gene in the T. spiralis and T. muris operons are trans-spliced to spliced leader RNAs, and we are able to detect polycistronic RNAs derived from these operons. Importantly, a putative intercistronic region from one of these potential enoplean operons confers polycistronic processing activity when expressed as part of a chimeric operon in Caenorhabditis elegans. We find that T. spiralis genes located in operons have an increased likelihood of having operonic C. elegans homologs. However, operon structure in terms of synteny and gene content is not tightly conserved between the two taxa, consistent with models of operon evolution. We have nevertheless identified putative operons conserved between Enoplea and Chromadorea. Our data suggest that operons and "spliced leader" (SL) trans-splicing predate the radiation of the nematode phylum, an inference which is supported by the phylogenetic profile of proteins known to be involved in nematode SL trans-splicing.

Pol II docking and pausing at growth and stress genes in C. elegans

Fluctuations in nutrient availability profoundly impact gene expression. Previous work revealed postrecruitment regulation of RNA polymerase II (Pol II) during starvation and recovery in Caenorhabditis elegans, suggesting that promoter-proximal pausing promotes rapid response to feeding. To test this hypothesis, we measured Pol II elongation genome wide by two complementary approaches and analyzed elongation in conjunction with Pol II binding and expression. We confirmed bona fide pausing during starvation and also discovered Pol II docking. Pausing occurs at active stress-response genes that become downregulated in response to feeding. In contrast, "docked" Pol II accumulates without initiating upstream of inactive growth genes that become rapidly upregulated upon feeding. Beyond differences in function and expression, these two sets of genes have different core promoter motifs, suggesting alternative transcriptional machinery. Our work suggests that growth and stress genes are both regulated postrecruitment during starvation but at initiation and elongation, respectively, coordinating gene expression with nutrient availability.

To grow or not to grow: nutritional control of development during Caenorhabditis elegans L1 arrest

It is widely appreciated that larvae of the nematode Caenorhabditis elegans arrest development by forming dauer larvae in response to multiple unfavorable environmental conditions. C. elegans larvae can also reversibly arrest development earlier, during the first larval stage (L1), in response to starvation. "L1 arrest" (also known as "L1 diapause") occurs without morphological modification but is accompanied by increased stress resistance. Caloric restriction and periodic fasting can extend adult lifespan, and developmental models are critical to understanding how the animal is buffered from fluctuations in nutrient availability, impacting lifespan. L1 arrest provides an opportunity to study nutritional control of development. Given its relevance to aging, diabetes, obesity and cancer, interest in L1 arrest is increasing, and signaling pathways and gene regulatory mechanisms controlling arrest and recovery have been characterized. Insulin-like signaling is a critical regulator, and it is modified by and acts through microRNAs. DAF-18/PTEN, AMP-activated kinase and fatty acid biosynthesis are also involved. The nervous system, epidermis, and intestine contribute systemically to regulation of arrest, but cell-autonomous signaling likely contributes to regulation in the germline. A relatively small number of genes affecting starvation survival during L1 arrest are known, and many of them also affect adult lifespan, reflecting a common genetic basis ripe for exploration. mRNA expression is well characterized during arrest, recovery, and normal L1 development, providing a metazoan model for nutritional control of gene expression. In particular, post-recruitment regulation of RNA polymerase II is under nutritional control, potentially contributing to a rapid and coordinated response to feeding. The phenomenology of L1 arrest will be reviewed, as well as regulation of developmental arrest and starvation survival by various signaling pathways and gene regulatory mechanisms.

Human PSENEN and U2AF1L4 genes are concertedly regulated by a genuine bidirectional promoter

Head-to-head genes with a short distance between their transcription start sites may constitute up to 10% of all genes in the genomes of various species. It was hypothesized that this intergenic space may represent bidirectional promoters which are able to initiate transcription of both genes, but the true bidirectionality was proved only for a few of them. We present experimental evidence that, according to several criteria, a 269 bp region located between the PSENEN and U2AF1L4 human genes is a genuine bidirectional promoter regulating a concerted divergent transcription of these genes. Concerted transcription of PSENEN and U2AF1L4 can be necessary for regulation of T-cell activity.

RNA processing enables predictable programming of gene expression

Complex interactions among genetic components often result in variable systemic performance in designed multigene systems. Using the bacterial clustered regularly interspaced short palindromic repeat (CRISPR) pathway we develop a synthetic RNA-processing platform, and show that efficient and specific cleavage of precursor mRNA enables reliable and predictable regulation of multigene operons. Physical separation of linked genetic elements by CRISPR-mediated cleavage is an effective strategy to achieve assembly of promoters, ribosome binding sites, cis-regulatory elements, and riboregulators into single- and multigene operons with predictable functions in bacteria. We also demonstrate that CRISPR-based RNA cleavage is effective for regulation in bacteria, archaea and eukaryotes. Programmable RNA processing using CRISPR offers a general approach for creating context-free genetic elements and can be readily used in the bottom-up construction of increasingly complex biological systems in a plug-and-play manner.

Synthesis of N²-modified 7-methylguanosine 5'-monophosphates as nematode translation inhibitors

Preparative scale synthesis of 14 new N(2)-modified mononucleotide 5' mRNA cap analogues was achieved. The key step involved use of an S(N)Ar reaction with protected 2-fluoro inosine and various primary and secondary amines. The derivatives were tested in a parasitic nematode, Ascaris suum, cell-free system as translation inhibitors. The most effective compound with IC(50) ∼0.9μM was a N(2)-p-metoxybenzyl-7-methylguanosine-5'-monophosphate 35.

PaxDb, a database of protein abundance averages across all three domains of life

Although protein expression is regulated both temporally and spatially, most proteins have an intrinsic, "typical" range of functionally effective abundance levels. These extend from a few molecules per cell for signaling proteins, to millions of molecules for structural proteins. When addressing fundamental questions related to protein evolution, translation and folding, but also in routine laboratory work, a simple rough estimate of the average wild type abundance of each detectable protein in an organism is often desirable. Here, we introduce a meta-resource dedicated to integrating information on absolute protein abundance levels; we place particular emphasis on deep coverage, consistent post-processing and comparability across different organisms. Publicly available experimental data are mapped onto a common namespace and, in the case of tandem mass spectrometry data, re-processed using a standardized spectral counting pipeline. By aggregating and averaging over the various samples, conditions and cell-types, the resulting integrated data set achieves increased coverage and a high dynamic range. We score and rank each contributing, individual data set by assessing its consistency against externally provided protein-network information, and demonstrate that our weighted integration exhibits more consistency than the data sets individually. The current PaxDb-release 2.1 (at http://pax-db.org/) presents whole-organism data as well as tissue-resolved data, and covers 85,000 proteins in 12 model organisms. All values can be seamlessly compared across organisms via pre-computed orthology relationships.

RNA-seq analysis of the C. briggsae transcriptome

Curation of a high-quality gene set is the critical first step in genome research, enabling subsequent analyses such as ortholog assignment, cis-regulatory element finding, and synteny detection. In this project, we have reannotated the genome of Caenorhabditis briggsae, the best studied sister species of the model organism Caenorhabditis elegans. First, we applied a homology-based gene predictor genBlastG to annotate the C. briggsae genome. We then validated and further improved the C. briggsae gene annotation through RNA-seq analysis of the C. briggsae transcriptome, which resulted in the first validated C. briggsae gene set (23,159 genes), among which 7347 genes (33.9% of all genes with introns) have all of their introns confirmed. Most genes (14,812, or 68.3%) have at least one intron validated, compared with only 3.9% in the most recent WormBase release (WS228). Of all introns in the revised gene set (103,083), 61,503 (60.1%) have been confirmed. Additionally, we have identified numerous trans-splicing leaders (SL1 and SL2 variants) in C. briggsae, leading to the first genome-wide annotation of operons in C. briggsae (1105 operons). The majority of the annotated operons (564, or 51.0%) are perfectly conserved in C. elegans, with an additional 345 operons (or 31.2%) somewhat divergent. Additionally, RNA-seq analysis revealed over 10 thousand small-size assembly errors in the current C. briggsae reference genome that can be readily corrected. The revised C. briggsae genome annotation represents a solid platform for comparative genomics analysis and evolutionary studies of Caenorhabditis species.

Incorporating genomics into the toolkit of nematology

The study of nematode genomes over the last three decades has relied heavily on the model organism Caenorhabditis elegans, which remains the best-assembled and annotated metazoan genome. This is now changing as a rapidly expanding number of nematodes of medical and economic importance have been sequenced in recent years. The advent of sequencing technologies to achieve the equivalent of the $1000 human genome promises that every nematode genome of interest will eventually be sequenced at a reasonable cost. As the sequencing of species spanning the nematode phylum becomes a routine part of characterizing nematodes, the comparative approach and the increasing use of ecological context will help us to further understand the evolution and functional specializations of any given species by comparing its genome to that of other closely and more distantly related nematodes. We review the current state of nematode genomics and discuss some of the highlights that these genomes have revealed and the trend and benefits of ecological genomics, emphasizing the potential for new genomes and the exciting opportunities this provides for nematological studies.

Nutritional control of mRNA isoform expression during developmental arrest and recovery in C. elegans

Nutrient availability profoundly influences gene expression. Many animal genes encode multiple transcript isoforms, yet the effect of nutrient availability on transcript isoform expression has not been studied in genome-wide fashion. When Caenorhabditis elegans larvae hatch without food, they arrest development in the first larval stage (L1 arrest). Starved larvae can survive L1 arrest for weeks, but growth and post-embryonic development are rapidly initiated in response to feeding. We used RNA-seq to characterize the transcriptome during L1 arrest and over time after feeding. Twenty-seven percent of detectable protein-coding genes were differentially expressed during recovery from L1 arrest, with the majority of changes initiating within the first hour, demonstrating widespread, acute effects of nutrient availability on gene expression. We used two independent approaches to track expression of individual exons and mRNA isoforms, and we connected changes in expression to functional consequences by mining a variety of databases. These two approaches identified an overlapping set of genes with alternative isoform expression, and they converged on common functional patterns. Genes affecting mRNA splicing and translation are regulated by alternative isoform expression, revealing post-transcriptional consequences of nutrient availability on gene regulation. We also found that phosphorylation sites are often alternatively expressed, revealing a common mode by which alternative isoform expression modifies protein function and signal transduction. Our results detail rich changes in C. elegans gene expression as larvae initiate growth and post-embryonic development, and they provide an excellent resource for ongoing investigation of transcriptional regulation and developmental physiology.