Comparative RNA-Seq analysis reveals pervasive tissue-specific alternative polyadenylation in Caenorhabditis elegans intestine and muscles

A spatiotemporally resolved atlas of mRNA decay in the C. elegans embryo reveals differential regulation of mRNA stability across stages and cell types

During embryonic development, cells undergo dynamic changes in gene expression that are required for appropriate cell fate specification. Although both transcription and mRNA degradation contribute to gene expression dynamics, patterns of mRNA decay are less well understood. Here, we directly measure spatiotemporally resolved mRNA decay rates transcriptome-wide throughout C. elegans embryogenesis by transcription inhibition followed by bulk and single-cell RNA sequencing. This allows us to calculate mRNA half-lives within specific cell types and developmental stages, and identify differentially regulated mRNA decay throughout embryonic development. We identify transcript features that are correlated with mRNA stability and find that mRNA decay rates are associated with distinct peaks in gene expression over time. Moreover, we provide evidence that, on average, mRNA is more stable in the germline than in the soma and in later embryonic stages than in earlier stages. This work suggests that differential mRNA decay across cell states and time helps to shape developmental gene expression, and it provides a valuable resource for studies of mRNA turnover regulatory mechanisms.

A comprehensive analysis of 3'UTRs in Caenorhabditis elegans

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

Single-tissue proteomics in Caenorhabditis elegans reveals proteins resident in intestinal lysosome-related organelles

The nematode intestine is the primary site for nutrient uptake and storage as well as the synthesis of biomolecules; lysosome-related organelles known as gut granules are important for many of these functions. Aspects of intestine biology are not well understood, including the export of the nutrients it imports and the molecules it synthesizes, as well as the complete functions and protein content of the gut granules. Here, we report a mass spectrometry (MS)-based proteomic analysis of the intestine of the Caenorhabditis elegans and of its gut granules. Overall, we identified approximately 5,000 proteins each in the intestine and the gonad and showed that most of these proteins can be detected in samples extracted from a single worm, suggesting the feasibility of individual-level genetic analysis using proteomes. Comparing proteomes and published transcriptomes of the intestine and the gonad, we identified proteins that appear to be synthesized in the intestine and then transferred to the gonad. To identify gut granule proteins, we compared the proteome of individual intestines deficient in gut granules to the wild type. The identified gut granule proteome includes proteins known to be exclusively localized to the granules and additional putative gut granule proteins. We selected two of these putative gut granule proteins for validation via immunohistochemistry, and our successful confirmation of both suggests that our strategy was effective in identifying the gut granule proteome. Our results demonstrate the practicability of single-tissue MS-based proteomic analysis in small organisms and in its future utility.

Transcriptome profiling of the Caenorhabditis elegans intestine reveals that ELT-2 negatively and positively regulates intestinal gene expression within the context of a gene regulatory network

ELT-2 is the major transcription factor (TF) required for Caenorhabditis elegans intestinal development. ELT-2 expression initiates in embryos to promote development and then persists after hatching through the larval and adult stages. Though the sites of ELT-2 binding are characterized and the transcriptional changes that result from ELT-2 depletion are known, an intestine-specific transcriptome profile spanning developmental time has been missing. We generated this dataset by performing Fluorescence Activated Cell Sorting on intestine cells at distinct developmental stages. We analyzed this dataset in conjunction with previously conducted ELT-2 studies to evaluate the role of ELT-2 in directing the intestinal gene regulatory network through development. We found that only 33% of intestine-enriched genes in the embryo were direct targets of ELT-2 but that number increased to 75% by the L3 stage. This suggests additional TFs promote intestinal transcription especially in the embryo. Furthermore, only half of ELT-2's direct target genes were dependent on ELT-2 for their proper expression levels, and an equal proportion of those responded to elt-2 depletion with over-expression as with under-expression. That is, ELT-2 can either activate or repress direct target genes. Additionally, we observed that ELT-2 repressed its own promoter, implicating new models for its autoregulation. Together, our results illustrate that ELT-2 impacts roughly 20-50% of intestine-specific genes, that ELT-2 both positively and negatively controls its direct targets, and that the current model of the intestinal regulatory network is incomplete as the factors responsible for directing the expression of many intestinal genes remain unknown.

stAPAminer: Mining Spatial Patterns of Alternative Polyadenylation for Spatially Resolved Transcriptomic Studies

Alternative polyadenylation (APA) contributes to transcriptome complexity and gene expression regulation and has been implicated in various cellular processes and diseases. Single-cell RNA sequencing (scRNA-seq) has enabled the profiling of APA at the single-cell level; however, the spatial information of cells is not preserved in scRNA-seq. Alternatively, spatial transcriptomics (ST) technologies provide opportunities to decipher the spatial context of the transcriptomic landscape. Pioneering studies have revealed potential spatially variable genes and/or splice isoforms; however, the pattern of APA usage in spatial contexts remains unappreciated. In this study, we developed a toolkit called stAPAminer for mining spatial patterns of APA from spatially barcoded ST data. APA sites were identified and quantified from the ST data. In particular, an imputation model based on the k-nearest neighbors algorithm was designed to recover APA signals, and then APA genes with spatial patterns of APA usage variation were identified. By analyzing well-established ST data of the mouse olfactory bulb (MOB), we presented a detailed view of spatial APA usage across morphological layers of the MOB. We compiled a comprehensive list of genes with spatial APA dynamics and obtained several major spatial expression patterns that represent spatial APA dynamics in different morphological layers. By extending this analysis to two additional replicates of the MOB ST data, we observed that the spatial APA patterns of several genes were reproducible among replicates. stAPAminer employs the power of ST to explore the transcriptional atlas of spatial APA patterns with spatial resolution. This toolkit is available at https://github.com/BMILAB/stAPAminer and https://ngdc.cncb.ac.cn/biocode/tools/BT007320.

Cultured Mesenchymal Cells from Nasal Turbinate as a Cellular Model of the Neurodevelopmental Component of Schizophrenia Etiology

Study of the neurodevelopmental molecular mechanisms of schizophrenia requires the development of adequate biological models such as patient-derived cells and their derivatives. We previously used cell lines with neural progenitor properties (CNON) derived from superior or middle turbinates of patients with schizophrenia and control groups to study gene expression specific to schizophrenia. In this study, we compared single cell-RNA seq data from two CNON cell lines, one derived from an individual with schizophrenia (SCZ) and the other from a control group, with two biopsy samples from the middle turbinate (MT), also from an individual with SCZ and a control. In addition, we compared our data with previously published data from olfactory neuroepithelium (1). Our data demonstrated that CNON originated from a single cell type which is present both in middle turbinate and olfactory neuroepithelium. CNON express multiple markers of mesenchymal cells. In order to define relatedness of CNON to the developing human brain, we also compared CNON datasets with scRNA-seq data of embryonic brain (2) and found that the expression profile of CNON very closely matched one of the cell types in the embryonic brain. Finally, we evaluated differences between SCZ and control samples to assess usability and potential benefits of using single cell RNA-seq of CNON to study etiology of schizophrenia.

Executing cell-specific cross-linking immunoprecipitation and sequencing (seCLIP) in C. elegans

The single-end enhanced cross-linking immunoprecipitation (seCLIP) method is well suited for efficient and unbiased transcriptome-wide interrogation of RNA-binding protein (RBP) interaction sites. Here, we provide a protocol for executing cell-specific seCLIP for any desired RBP in Caenorhabditis elegans. We begin with steps and recommendations for transgene construction and Cas9-mediated chromosomal integration. We provide detailed procedures for isolation of RBP-associated RNA fragments, subsequent library preparation, and sequencing. We further discuss best practices for data analysis, interpretation of results, and troubleshooting. For complete details on the use and execution of this protocol, please refer to Blazie et al. (2021).1.

An updated C. elegans nuclear body muscle transcriptome for studies in muscle formation and function

The body muscle is an important tissue used in organisms for proper viability and locomotion. Although this tissue is generally well studied and characterized, and many pathways have been elucidated throughout the years, we still lack a comprehensive understanding of its transcriptome and how it controls muscle development and function. Here, we have updated a nuclear FACS sorting-based methodology to isolate and sequence a high-quality muscle transcriptome from Caenorhabditis elegans mixed-stage animals. We have identified 2848 muscle-specific protein-coding genes, including 78 transcription factors and 206 protein-coding genes containing an RNA binding domain. We studied their interaction network, performed a detailed promoter analysis, and identified novel muscle-specific cis-acting elements. We have also identified 16 high-quality muscle-specific miRNAs, studied their function in vivo using fluorochrome-based analyses, and developed a high-quality C. elegans miRNA interactome incorporating other muscle-specific datasets produced by our lab and others.Our study expands our understanding of how muscle tissue functions in C. elegans andin turn provides results that can in the future be applied to humans to study muscular-related diseases.

A Survey on Methods for Predicting Polyadenylation Sites from DNA Sequences, Bulk RNA-seq, and Single-cell RNA-seq

Alternative polyadenylation (APA) plays important roles in modulating mRNA stability, translation, and subcellular localization, and contributes extensively to shaping eukaryotic transcriptome complexity and proteome diversity. Identification of poly(A) sites (pAs) on a genome-wide scale is a critical step toward understanding the underlying mechanism of APA-mediated gene regulation. A number of established computational tools have been proposed to predict pAs from diverse genomic data. Here we provided an exhaustive overview of computational approaches for predicting pAs from DNA sequences, bulk RNA sequencing (RNA-seq) data, and single-cell RNA sequencing (scRNA-seq) data. Particularly, we examined several representative tools using bulk RNA-seq and scRNA-seq data from peripheral blood mononuclear cells and put forward operable suggestions on how to assess the reliability of pAs predicted by different tools. We also proposed practical guidelines on choosing appropriate methods applicable to diverse scenarios. Moreover, we discussed in depth the challenges in improving the performance of pA prediction and benchmarking different methods. Additionally, we highlighted outstanding challenges and opportunities using new machine learning and integrative multi-omics techniques, and provided our perspective on how computational methodologies might evolve in the future for non-3' untranslated region, tissue-specific, cross-species, and single-cell pA prediction.

Identifying the Caenorhabditis elegans vulval transcriptome

Development of the Caenorhabditis elegans vulva is a classic model of organogenesis. This system, which starts with 6 equipotent cells, encompasses diverse types of developmental event, including developmental competence, multiple signaling events to control precise and faithful patterning of three cell fates, execution and proliferation of specific cell lineages, and a series of sophisticated morphogenetic events. Early events have been subjected to extensive mutational and genetic investigations and later events to cell biological analyses. We infer the existence of dramatically changing profiles of gene expression that accompanies the observed changes in development. Yet, except from serendipitous discovery of several transcription factors expressed in dynamic patterns in vulval lineages, our knowledge of the transcriptomic landscape during vulval development is minimal. This study describes the composition of a vulva-specific transcriptome. We used tissue-specific harvesting of mRNAs via immunoprecipitation of epitope-tagged poly(A) binding protein, PAB-1, heterologously expressed by a promoter known to express GFP in vulval cells throughout their development. The identified transcriptome was small but tightly interconnected. From this data set, we identified several genes with identified functions in development of the vulva and validated more with promoter-GFP reporters of expression. For one target, lag-1, promoter-GFP expression was limited but a fluorescent tag of the endogenous protein revealed extensive expression. Thus, we have identified a transcriptome of C. elegans vulval lineages as a launching pad for exploration of functions of these genes in organogenesis.

Endomitosis controls tissue-specific gene expression during development

Polyploid cells contain more than 2 copies of the genome and are found in many plant and animal tissues. Different types of polyploidy exist, in which the genome is confined to either 1 nucleus (mononucleation) or 2 or more nuclei (multinucleation). Despite the widespread occurrence of polyploidy, the functional significance of different types of polyploidy is largely unknown. Here, we assess the function of multinucleation in Caenorhabditis elegans intestinal cells through specific inhibition of binucleation without altering genome ploidy. Through single-worm RNA sequencing, we find that binucleation is important for tissue-specific gene expression, most prominently for genes that show a rapid up-regulation at the transition from larval development to adulthood. Regulated genes include vitellogenins, which encode yolk proteins that facilitate nutrient transport to the germline. We find that reduced expression of vitellogenins in mononucleated intestinal cells leads to progeny with developmental delays and reduced fitness. Together, our results show that binucleation facilitates rapid up-regulation of intestine-specific gene expression during development, independently of genome ploidy, underscoring the importance of spatial genome organization for polyploid cell function.

Why do some cells contain more than one nucleus? By comparing mononucleated and multinucleated polyploid cells in C. elegans, this study shows that having multiple nuclei is important for optimal transcriptional upregulation of developmentally controlled genes.

Implications of Poly(A) Tail Processing in Repeat Expansion Diseases

Repeat expansion diseases are a group of more than 40 disorders that affect mainly the nervous and/or muscular system and include myotonic dystrophies, Huntington’s disease, and fragile X syndrome. The mutation-driven expanded repeat tract occurs in specific genes and is composed of tri- to dodeca-nucleotide-long units. Mutant mRNA is a pathogenic factor or important contributor to the disease and has great potential as a therapeutic target. Although repeat expansion diseases are quite well known, there are limited studies concerning polyadenylation events for implicated transcripts that could have profound effects on transcript stability, localization, and translation efficiency. In this review, we briefly present polyadenylation and alternative polyadenylation (APA) mechanisms and discuss their role in the pathogenesis of selected diseases. We also discuss several methods for poly(A) tail measurement (both transcript-specific and transcriptome-wide analyses) and APA site identification—the further development and use of which may contribute to a better understanding of the correlation between APA events and repeat expansion diseases. Finally, we point out some future perspectives on the research into repeat expansion diseases, as well as APA studies.

Tissue-specific DamID protocol using nanopore sequencing

DNA adenine methylation identification (DamID) is a powerful method to determine DNA binding profiles of proteins at a genomic scale. The method leverages the fusion between a protein of interest and the Dam methyltransferase of E. coli, which methylates proximal DNA in vivo. Here, we present an optimized procedure, which was developed for tissue-specific analyses in Caenorhabditis elegans and successfully used to footprint genes actively transcribed by RNA polymerases and to map transcription factor binding in gene regulatory regions. The present protocol details C. elegans-specific steps involved in the preparation of transgenic lines and genomic DNA samples, as well as broadly applicable steps for the DamID procedure, including the isolation of methylated DNA fragments, the preparation of multiplexed libraries, Nanopore sequencing, and data analysis. Two distinctive features of the approach are (i) the use of an efficient recombination-based strategy to selectively analyze rare cell types and (ii) the use of Nanopore sequencing, which streamlines the process. The method allows researchers to go from genomic DNA samples to sequencing results in less than a week, while being sensitive enough to report reliable DNA footprints in cell types as rare as 2 cells per animal.

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

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

Transcriptome changes during the initiation and progression of Duchenne muscular dystrophy in Caenorhabditis elegans

Duchenne muscular dystrophy (DMD) is a lethal, X-linked disease characterized by progressive muscle degeneration. The condition is driven by nonsense and missense mutations in the dystrophin gene, leading to instability of the sarcolemma and skeletal muscle necrosis and atrophy. Resulting changes in muscle-specific gene expression that take place in dystrophin's absence remain largely uncharacterized, as they are potentially obscured by the chronic inflammation elicited by muscle damage in humans. Caenorhabditis elegans possess a mild inflammatory response that is not active in the muscle, and lack a satellite cell equivalent. This allows for the characterization of the transcriptome rearrangements affecting disease progression independently of inflammation and regeneration. In effort to better understand these dynamics, we have isolated and sequenced body muscle-specific transcriptomes from C. elegans lacking functional dystrophin at distinct stages of disease progression. We have identified an upregulation of genes involved in mitochondrial function early in disease progression, and an upregulation of genes related to muscle repair in later stages. Our results suggest that in C. elegans, dystrophin may have a signaling role early in development, and its absence may activate compensatory mechanisms that counteract muscle degradation caused by loss of dystrophin. We have also developed a temperature-based screening method for synthetic paralysis that can be used to rapidly identify genetic partners of dystrophin. Our results allow for the comprehensive identification of transcriptome changes that potentially serve as independent drivers of disease progression and may in turn allow for the identification of new therapeutic targets for the treatment of DMD.

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.

Muscle-derived exophers promote reproductive fitness

Organismal functionality and reproduction depend on metabolic rewiring and balanced energy resources. However, the crosstalk between organismal homeostasis and fecundity and the associated paracrine signaling mechanisms are still poorly understood. Using Caenorhabditis elegans, we discovered that large extracellular vesicles (known as exophers) previously found to remove damaged subcellular elements in neurons and cardiomyocytes are released by body wall muscles (BWM) to support embryonic growth. Exopher formation (exopheresis) by BWM is sex-specific and a non-cell autonomous process regulated by developing embryos in the uterus. Embryo-derived factors induce the production of exophers that transport yolk proteins produced in the BWM and ultimately deliver them to newly formed oocytes. Consequently, offspring of mothers with a high number of muscle-derived exophers grew faster. We propose that the primary role of muscular exopheresis is to stimulate reproductive capacity, thereby influencing the adaptation of worm populations to the current environmental conditions.

Caenorhabditis elegans as a Model System for Duchenne Muscular Dystrophy

The nematode worm Caenorhabditis elegans has been used extensively to enhance our understanding of the human neuromuscular disorder Duchenne Muscular Dystrophy (DMD). With new arising clinically relevant models, technologies and treatments, there is a need to reconcile the literature and collate the key findings associated with this model.

Post-Transcriptional Regulation in Skeletal Muscle Development, Repair, and Disease

Skeletal muscle formation is a complex process that requires tight spatiotemporal control of key myogenic factors. Emerging evidence suggests that RNA processing is crucial for the regulation of these factors, and that multiple post-transcriptional regulatory pathways work dependently and independently of one another to enable precise control of transcripts throughout muscle development and repair. Moreover, disruption of these pathways is implicated in neuromuscular disease, and the recent development of RNA-mediated therapies shows enormous promise in the treatment of these disorders. We discuss the overlapping post-transcriptional regulatory pathways that mediate muscle development, how these pathways are disrupted in neuromuscular disorders, and advances in RNA-mediated therapies that present a novel approach to the treatment of these diseases.

A primer for generating and using transcriptome data and gene sets

Transcriptomic approaches have provided a growing set of powerful tools with which to study genome-wide patterns of gene expression. Rapidly evolving technologies enable analysis of transcript abundance data from particular tissues and even single cells. This Primer discusses methods that can be used to collect and profile RNAs from specific tissues or cells, process and analyze high-throughput RNA-sequencing data, and define sets of genes that accurately represent a category, such as tissue-enriched or tissue-specific gene expression.

Summary: This Primer gives an overview of the considerations required when planning and conducting transcriptome profiling experiments, from sample isolation to data analysis and display, and discusses best practice in defining gene sets for use in addressing biological questions.

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

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

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.

TEnGExA: an R package based tool for tissue enrichment and gene expression analysis

RNA-seq data analysis with rapidly advancing high-throughput sequencing technology, nowadays provides large number of transcripts or genes to perform downstream analysis including functional annotation and pathway analysis. However for the data from multiple tissues, downstream analysis with tissue-specific or tissue-enriched transcripts is highly preferable. However, there is still a need of tool for quickly performing tissue-enrichment and gene expression analysis irrespective of number of input genes or tissues at various fragments per kilobase of transcript per million fragments mapped (FPKM) thresholds. To fulfill this need, we presented a freely available R package and web-interface tool, TEnGExA, which allows tissue-enrichment analysis (TEA) for any number of genes or transcripts for any species provided only a read-count or FPKM-value matrix as input. Based on the different FPKM value and fold thresholds, TEnGExA classifies the user provided gene lists into tissue-enriched or tissue-specific transcripts along with other standard classes. By analyzing the published sample data from human, plant and microorganism, we signifies that TEnGExA can easily handle complex or large data from any species to provided tissue-enriched gene list for downstream analysis in quick time. In summary, TEnGExA is quick, easy to use and an efficient tool for TEA. The R package is freely available at https://github.com/ubagithub/TEnGExA/ and the GUI web interface is accessible at http://webtom.cabgrid.res.in/tissue_enrich/.

Gene expression dynamics are a proxy for selective pressures on alternatively polyadenylated isoforms

Alternative polyadenylation (APA) produces isoforms with distinct 3′-ends, yet their functional differences remain largely unknown. Here, we introduce the APA-seq method to detect the expression levels of APA isoforms from 3′-end RNA-Seq data by exploiting both paired-end reads for gene isoform identification and quantification. We detected the expression levels of APA isoforms in individual Caenorhabditis elegans embryos at different stages throughout embryogenesis. Examining the correlation between the temporal profiles of isoforms led us to distinguish two classes of genes: those with highly correlated isoforms (HCI) and those with lowly correlated isoforms (LCI) across time. We hypothesized that variants with similar expression profiles may be the product of biological noise, while the LCI variants may be under tighter selection and consequently their distinct 3′ UTR isoforms are more likely to have functional consequences. Supporting this notion, we found that LCI genes have significantly more miRNA binding sites, more correlated expression profiles with those of their targeting miRNAs and a relative lack of correspondence between their transcription and protein abundances. Collectively, our results suggest that a lack of coherence among the regulation of 3′ UTR isoforms is a proxy for selective pressures acting upon APA usage and consequently for their functional relevance.

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.

Wide-ranging transcriptome remodelling mediated by alternative polyadenylation in response to abiotic stresses in Sorghum

Alternative polyadenylation (APA) regulates diverse developmental and physiological processes through its effects on gene expression, mRNA stability, translatability, and transport. Sorghum is a major cereal crop in the world and, despite its importance, not much is known about the role of post-transcriptional regulation in mediating responses to abiotic stresses in Sorghum. A genome-wide APA analysis unveiled widespread occurrence of APA in Sorghum in response to drought, heat, and salt stress. Abiotic stress treatments incited changes in poly(A) site choice in a large number of genes. Interestingly, abiotic stresses led to the re-directing of transcriptional output into non-productive pathways defined by the class of poly(A) site utilized. This result revealed APA to be part of a larger global response of Sorghum to abiotic stresses that involves the re-direction of transcriptional output into non-productive transcriptional and translational pathways. Large numbers of stress-inducible poly(A) sites could not be linked with known, annotated genes, suggestive of the existence of numerous unidentified genes whose expression is strongly regulated by abiotic stresses. Furthermore, we uncovered a novel stress-specific cis-element in intronic poly(A) sites used in drought- and heat-stressed plants that might play an important role in non-canonical poly(A) site choice in response to abiotic stresses.

Loss of an H3K9me anchor rescues laminopathy-linked changes in nuclear organization and muscle function in an Emery-Dreifuss muscular dystrophy model

In this study, Harr et al. use C. elegans to investigate the consequences of a missense mutation (Y45C) in lamin A (encoded by LMNA) found in the human Emery-Dreifuss muscular dystrophy (EDMD) syndrome. Using muscle-specific emerin Dam-ID and other in vivo approaches, the authors report that they were able to counteract the dominant muscle-specific defects provoked by LMNA mutation by the ablation of a lamin-associated H3K9me anchor, suggesting a novel therapeutic pathway for treating EDMD.

Mutations in the nuclear structural protein lamin A produce rare, tissue-specific diseases called laminopathies. The introduction of a human Emery-Dreifuss muscular dystrophy (EDMD)-inducing mutation into the C. elegans lamin (LMN-Y59C), recapitulates many muscular dystrophy phenotypes, and correlates with hyper-sequestration of a heterochromatic array at the nuclear periphery in muscle cells. Using muscle-specific emerin Dam-ID in worms, we monitored the effects of the mutation on endogenous chromatin. An increased contact with the nuclear periphery along chromosome arms, and an enhanced release of chromosomal centers, coincided with the disease phenotypes of reduced locomotion and compromised sarcomere integrity. The coupling of the LMN-Y59C mutation with the ablation of CEC-4, a chromodomain protein that anchors H3K9-methylated chromatin at the nuclear envelope (NE), suppressed the muscle-associated disease phenotypes. Deletion of cec-4 also rescued LMN-Y59C-linked alterations in chromatin organization and some changes in transcription. Sequences that changed position in the LMN-Y59C mutant, are enriched for E2F (EFL-2)-binding sites, consistent with previous studies suggesting that altered Rb-E2F interaction with lamin A may contribute to muscle dysfunction. In summary, we were able to counteract the dominant muscle-specific defects provoked by LMNA mutation by the ablation of a lamin-associated H3K9me anchor, suggesting a novel therapeutic pathway for EDMD.

Contributions of alternative splicing to muscle type development and function

Animals possess a wide variety of muscle types that support different kinds of movements. Different muscles have distinct locations, morphologies and contractile properties, raising the question of how muscle diversity is generated during development. Normal aging processes and muscle disorders differentially affect particular muscle types, thus understanding how muscles normally develop and are maintained provides insight into alterations in disease and senescence. As muscle structure and basic developmental mechanisms are highly conserved, many important insights into disease mechanisms in humans as well as into basic principles of muscle development have come from model organisms such as Drosophila, zebrafish and mouse. While transcriptional regulation has been characterized to play an important role in myogenesis, there is a growing recognition of the contributions of alternative splicing to myogenesis and the refinement of muscle function. Here we review our current understanding of muscle type specific alternative splicing, using examples of isoforms with distinct functions from both vertebrates and Drosophila. Future exploration of the vast potential of alternative splicing to fine-tune muscle development and function will likely uncover novel mechanisms of isoform-specific regulation and a more holistic understanding of muscle development, disease and aging.

Transcriptional Profiling of C. elegans Adult Cells and Tissues with Age

Multicellular organisms are composed of distinct cells and tissues that coordinate highly orchestrated responses to environmental challenges, including those that arise with age. Since C. elegans is a premier model system used to study the molecular and cellular regulators of adult and aging phenotypes, cell type- and tissue-specific approaches are needed to characterize the genome-wide expression changes associated with these responses. Here we describe a method for the FACS-based isolation and RNA sequencing of dissociated cells from adult C. elegans. This technique is amenable to profiling the cell- and tissue-specific gene expression changes in C. elegans mutants, including aging models, such as the daf-2/insulin-like signaling (IIS) pathway, and in wild-type animals with age.

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

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

A survey on identification and quantification of alternative polyadenylation sites from RNA-seq data

Alternative polyadenylation (APA) has been implicated to play an important role in post-transcriptional regulation by regulating mRNA abundance, stability, localization and translation, which contributes considerably to transcriptome diversity and gene expression regulation. RNA-seq has become a routine approach for transcriptome profiling, generating unprecedented data that could be used to identify and quantify APA site usage. A number of computational approaches for identifying APA sites and/or dynamic APA events from RNA-seq data have emerged in the literature, which provide valuable yet preliminary results that should be refined to yield credible guidelines for the scientific community. In this review, we provided a comprehensive overview of the status of currently available computational approaches. We also conducted objective benchmarking analysis using RNA-seq data sets from different species (human, mouse and Arabidopsis) and simulated data sets to present a systematic evaluation of 11 representative methods. Our benchmarking study showed that the overall performance of all tools investigated is moderate, reflecting that there is still lot of scope to improve the prediction of APA site or dynamic APA events from RNA-seq data. Particularly, prediction results from individual tools differ considerably, and only a limited number of predicted APA sites or genes are common among different tools. Accordingly, we attempted to give some advice on how to assess the reliability of the obtained results. We also proposed practical recommendations on the appropriate method applicable to diverse scenarios and discussed implications and future directions relevant to profiling APA from RNA-seq data.

ALG-1 Influences Accurate mRNA Splicing Patterns in the Caenorhabditis elegans Intestine and Body Muscle Tissues by Modulating Splicing Factor Activities

MicroRNAs (miRNAs) are known to modulate gene expression, but their activity at the tissue-specific level remains largely uncharacterized. To study their contribution to tissue-specific gene expression, we developed novel tools to profile putative miRNA targets in the Caenorhabditis elegans intestine and body muscle. We validated many previously described interactions and identified ∼3500 novel targets. Many of the candidate miRNA targets curated are known to modulate the functions of their respective tissues. Within our data sets we observed a disparity in the use of miRNA-based gene regulation between the intestine and body muscle. The intestine contained significantly more putative miRNA targets than the body muscle highlighting its transcriptional complexity. We detected an unexpected enrichment of RNA-binding proteins targeted by miRNA in both tissues, with a notable abundance of RNA splicing factors. We developed in vivo genetic tools to validate and further study three RNA splicing factors identified as putative miRNA targets in our study (asd-2, hrp-2, and smu-2), and show that these factors indeed contain functional miRNA regulatory elements in their 3'UTRs that are able to repress their expression in the intestine. In addition, the alternative splicing pattern of their respective downstream targets (unc-60, unc-52, lin-10, and ret-1) is dysregulated when the miRNA pathway is disrupted. A reannotation of the transcriptome data in C. elegans strains that are deficient in the miRNA pathway from past studies supports and expands on our results. This study highlights an unexpected role for miRNAs in modulating tissue-specific gene isoforms, where post-transcriptional regulation of RNA splicing factors associates with tissue-specific alternative splicing.

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.

APAtrap: identification and quantification of alternative polyadenylation sites from RNA-seq data

Motivation

Alternative polyadenylation (APA) has been increasingly recognized as a crucial mechanism that contributes to transcriptome diversity and gene expression regulation. As RNA-seq has become a routine protocol for transcriptome analysis, it is of great interest to leverage such unprecedented collection of RNA-seq data by new computational methods to extract and quantify APA dynamics in these transcriptomes. However, research progress in this area has been relatively limited. Conventional methods rely on either transcript assembly to determine transcript 3' ends or annotated poly(A) sites. Moreover, they can neither identify more than two poly(A) sites in a gene nor detect dynamic APA site usage considering more than two poly(A) sites.

Results

We developed an approach called APAtrap based on the mean squared error model to identify and quantify APA sites from RNA-seq data. APAtrap is capable of identifying novel 3' UTRs and 3' UTR extensions, which contributes to locating potential poly(A) sites in previously overlooked regions and improving genome annotations. APAtrap also aims to tally all potential poly(A) sites and detect genes with differential APA site usages between conditions. Extensive comparisons of APAtrap with two other latest methods, ChangePoint and DaPars, using various RNA-seq datasets from simulation studies, human and Arabidopsis demonstrate the efficacy and flexibility of APAtrap for any organisms with an annotated genome.

Availability and implementation

Freely available for download at https://apatrap.sourceforge.io.

Contact

liqq@xmu.edu.cn or xhuister@xmu.edu.cn.

Supplementary information

Supplementary data are available at Bioinformatics online.

Conserved functions of RNA-binding proteins in muscle

Animals require different types of muscle for survival, for example for circulation, motility, reproduction and digestion. Much emphasis in the muscle field has been placed on understanding how transcriptional regulation generates diverse types of muscle during development. Recent work indicates that alternative splicing and RNA regulation are as critical to muscle development, and altered function of RNA-binding proteins causes muscle disease. Although hundreds of genes predicted to bind RNA are expressed in muscles, many fewer have been functionally characterized. We present a cross-species view summarizing what is known about RNA-binding protein function in muscle, from worms and flies to zebrafish, mice and humans. In particular, we focus on alternative splicing regulated by the CELF, MBNL and RBFOX families of proteins. We discuss the systemic nature of diseases associated with loss of RNA-binding proteins in muscle, focusing on mis-regulation of CELF and MBNL in myotonic dystrophy. These examples illustrate the conservation of RNA-binding protein function and the marked utility of genetic model systems in understanding mechanisms of RNA regulation.

Ectopic Germ Cells Can Induce Niche-like Enwrapment by Neighboring Body Wall Muscle

Niche cell enwrapment of stem cells and their differentiating progeny is common and provides a specialized signaling and protective environment. Elucidating the mechanisms underlying enwrapment behavior has important basic and clinical significance in not only understanding how niches are formed and maintained but also how they can be engineered and how they are misregulated in human pathologies, such as cancer. Previous work in C. elegans found that, when germ cells, which are enwrapped by somatic gonadal niche cells, are freed into the body cavity, they embed into other tissues. We investigated this phenomenon using live-cell imaging and discovered that ectopic germ cells preferentially induce body-wall muscle to extend cellular processes that enwrap the germ cells, the extent of which was strikingly similar to the distal tip cell (DTC)-germ stem cell niche. Enwrapment was specific for escaped germ cells, and genetic analysis revealed it did not depend on pathways that control cell death and engulfment or muscle arm extension. Instead, using a large-scale RNAi screen and GFP knockin strains, we discovered that the enwrapping behavior of muscle relied upon the same suite of cell-cell adhesion molecules that functioned in the endogenous niche: the C. elegans E-cadherin HMR-1, its intracellular associates α-catenin (HMP-1) and β-catenin (HMP-2), and the L1CAM protein SAX-7. This ectopic niche-like behavior resembles the seed-and-soil model of cancer metastasis and offers a new model to understand factors regulating ectopic niche formation.

RNA processing in skeletal muscle biology and disease

RNA processing encompasses the capping, cleavage, polyadenylation and alternative splicing of pre-mRNA. Proper muscle development relies on precise RNA processing, driven by the coordination between RNA-binding proteins. Recently, skeletal muscle biology has been intensely investigated in terms of RNA processing. High throughput studies paired with deletion of RNA-binding proteins have provided a high-level understanding of the molecular mechanisms controlling the regulation of RNA-processing in skeletal muscle. Furthermore, misregulation of RNA processing is implicated in muscle diseases. In this review, we comprehensively summarize recent studies in skeletal muscle that demonstrated: (i) the importance of RNA processing, (ii) the RNA-binding proteins that are involved, and (iii) diseases associated with defects in RNA processing.

Inverted repeat structures are associated with essential and highly expressed genes on C. elegans autosome distal arms

Complementary sequences in cellular transcripts base-pair to form double-stranded RNA (dsRNA) structures. Because transposon-derived repeats often give rise to self-complementary sequences, dsRNA structures are prevalent in eukaryotic genomes, typically occurring in gene introns and untranslated regions (UTRs). However, the regulatory impact of double-stranded structures within genes is not fully understood. We used three independent methods to define loci in Caenorhabditis elegans predicted to form dsRNA and correlated these structures with patterns of gene expression, gene essentiality, and genome organization. As previously observed, dsRNA loci are enriched on distal arms of C. elegans autosomes, where genes typically show less conservation and lower overall expression. In contrast, we find that dsRNAs are associated with essential genes on autosome arms, and dsRNA-associated genes exhibit higher-than-expected expression and histone modification patterns associated with transcriptional elongation. Genes with significant repetitive sequence content are also highly expressed, and, thus, observed gene expression trends may relate either to dsRNA structures or to repeat content. Our results raise the possibility that as-yet-undescribed mechanisms promote expression of loci that produce dsRNAs, despite their well-characterized roles in gene silencing.

Comparative transcriptomics identifies genes differentially expressed in the intestine of a new fast-growing strain of common carp with higher unsaturated fatty acid content in muscle

We have created a new, fast-growing strain of common carp with higher unsaturated fatty acid content in muscle. To better understand the impacts of gene regulation in intestinal tissue on growth and unsaturated fatty acid content, we conducted a comparative RNA-Seq transcriptome analysis between intestine samples of Selected and Control groups (and corroborated selected results by PCR). After eight weeks of cage culture, weight gain of the Selected group was 20.84% higher. In muscles of the control group, monounsaturated fatty acids (FAs) were more abundant, whereas polyunsaturated FAs were more abundant in muscles of the Selected group. In total, we found 106 differentially expressed genes (DEGs) between the two groups. Only the endocytosis pathway was significantly enriched in DEGs, with two upregulated genes: il2rb and ehd1. The latter is involved in the growth hormone/insulin-like growth factor (Gh/Igf) axis, which plays a key role in the regulation of growth in animals. tll2, which is known to be associated with intestinal regeneration, was extremely highly upregulated in both transcriptomic (infinite) and qPCR (610.70) analyses. Two of the upregulated genes are associated with the fatty acid metabolism, several genes are likely to be indicators of heightened transcription levels, several are associated with metabolic and developmental roles, several with neuronal functions (including two with vision), several with the immune system, and two downregulated genes with the development of vasculature. The higher growth rate of the Selected group is likely to be at least partially attributed to increased endocytosis efficiency and genetically-driven behavioural differences (higher aggression levels). There are some indications that this new strain might have slightly impaired immune responses, and a higher propensity for inherited diseases leading to sight impairment, as well for neurodegenerative diseases in general, but these indications still need to be confirmed.

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.

Effects of liquid cultivation on gene expression and phenotype of C. elegans

BACKGROUND

Liquid cultures have been commonly used in space, toxicology, and pharmacology studies of Caenorhabditis elegans. However, the knowledge about transcriptomic alterations caused by liquid cultivation remains limited. Moreover, the impact of different genotypes in rapid adaptive responses to environmental changes (e.g., liquid cultivation) is often overlooked. Here, we report the transcriptomic and phenotypic responses of laboratory N2 and the wild-isolate AB1 strains after culturing P0 worms on agar plates, F1 in liquid cultures, and F2 back on agar plates.

RESULTS

Significant variations were found in the gene expressions between the N2 and AB1 strains in response to liquid cultivation. The results demonstrated that 8-34% of the environmental change-induced transcriptional responses are transmitted to the subsequent generation. By categorizing the gene expressions for genotype, environment, and genotype-environment interactions, we identified that the genotype has a substantial impact on the adaptive responses. Functional analysis of the transcriptome showed correlation with phenotypical changes. For example, the N2 strain exhibited alterations in both phenotype and gene expressions for germline and cuticle in axenic liquid cultivation. We found transcript evidence to approximately 21% of the computationally predicted genes in C. elegans by exposing the worms to environmental changes.

CONCLUSIONS

The presented study reveals substantial differences between N2 and AB1 strains for transcriptomic and phenotypical responses to rapid environmental changes. Our data can provide standard controls for future studies for the liquid cultivation of C. elegans and enable the discovery of condition-specific genes.

Genome-wide atlas of alternative polyadenylation in the forage legume red clover

Studies on prevalence and significance of alternative polyadenylation (APA) in plants have been so far limited mostly to the model plants. Here, a genome-wide analysis of APA was carried out in different tissue types in the non-model forage legume red clover (Trifolium pratense L). A profile of poly(A) sites in different tissue types was generated using so-called 'poly(A)-tag sequencing' (PATseq) approach. Our analysis revealed tissue-wise dynamics of usage of poly(A) sites located at different genomic locations. We also identified poly(A) sites and underlying genes displaying APA in different tissues. Functional categories enriched in groups of genes manifesting APA between tissue types were determined. Analysis of spatial expression of genes encoding different poly(A) factors showed significant differential expression of genes encoding orthologs of FIP1(V) and PCFS4, suggesting that these two factors may play a role in regulating spatial APA in red clover. Our analysis also revealed a high degree of conservation in diverse plant species of APA events in mRNAs encoding two key polyadenylation factors, CPSF30 and FIP1(V). Together with our previously reported study of spatial gene expression in red clover, this study will provide a comprehensive account of transcriptome dynamics in this non-model forage legume.

Developmental dynamics of gene expression and alternative polyadenylation in the Caenorhabditis elegans germline

BACKGROUND

The 3' untranslated regions (UTRs) of mRNAs play a major role in post-transcriptional regulation of gene expression. Selection of transcript cleavage and polyadenylation sites is a dynamic process that produces multiple transcript isoforms for the same gene within and across different cell types. Using LITE-Seq, a new quantitative method to capture transcript 3' ends expressed in vivo, we have characterized sex- and cell type-specific transcriptome-wide changes in gene expression and 3'UTR diversity in Caenorhabditis elegans germline cells undergoing proliferation and differentiation.

RESULTS

We show that nearly half of germline transcripts are alternatively polyadenylated, that differential regulation of endogenous 3'UTR variants is common, and that alternative isoforms direct distinct spatiotemporal protein expression patterns in vivo. Dynamic expression profiling also reveals temporal regulation of X-linked gene expression, selective stabilization of transcripts, and strong evidence for a novel developmental program that promotes nucleolar dissolution in oocytes. We show that the RNA-binding protein NCL-1/Brat is a posttranscriptional regulator of numerous ribosome-related transcripts that acts through specific U-rich binding motifs to down-regulate mRNAs encoding ribosomal protein subunits, rRNA processing factors, and tRNA synthetases.

CONCLUSIONS

These results highlight the pervasive nature and functional potential of patterned gene and isoform expression during early animal development.

Early Developmental Exposure to dsRNA Is Critical for Initiating Efficient Nuclear RNAi in C. elegans

RNAi has enabled researchers to study the function of many genes. However, it is not understood why some RNAi experiments succeed while others do not. Here, we show in C. elegans that pharyngeal muscle is resistant to RNAi when initially exposed to double-stranded RNA (dsRNA) by feeding but sensitive to RNAi in the next generation. Investigating this observation, we find that pharyngeal muscle cells as well as vulval muscle cells require nuclear rather than cytoplasmic RNAi. Further, we find in these cell types that nuclear RNAi silencing is most efficiently triggered during early development, defining a critical period for initiating nuclear RNAi. Finally, using heat-shock-induced dsRNA expression, we show that synMuv B class mutants act in part to extend this critical window. The synMuv-B-dependent early-development-associated critical period for initiating nuclear RNAi suggests that mechanisms that restrict developmental plasticity may also restrict the initiation of nuclear RNAi.

Identification of long non-coding RNA in the horse transcriptome

Background

Efforts to resolve the transcribed sequences in the equine genome have focused on protein-coding RNA. The transcription of the intergenic regions, although detected via total RNA sequencing (RNA-seq), has yet to be characterized in the horse. The most recent equine transcriptome based on RNA-seq from several tissues was a prime opportunity to obtain a concurrent long non-coding RNA (lncRNA) database.

Results

This lncRNA database has a breadth of eight tissues and a depth of over 20 million reads for select tissues, providing the deepest and most expansive equine lncRNA database. Utilizing the intergenic reads and three categories of novel genes from a previously published equine transcriptome pipeline, we better describe these groups by annotating the lncRNA candidates. These lncRNA candidates were filtered using an approach adapted from human lncRNA annotation, which removes transcripts based on size, expression, protein-coding capability and distance to the start or stop of annotated protein-coding transcripts.

Conclusion

Our equine lncRNA database has 20,800 transcripts that demonstrate characteristics unique to lncRNA including low expression, low exon diversity and low levels of sequence conservation. These candidate lncRNA will serve as a baseline lncRNA annotation and begin to describe the RNA-seq reads assigned to the intergenic space in the horse.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3884-2) contains supplementary material, which is available to authorized users.

Cell type-specific transcriptome profiling in C. elegans using the Translating Ribosome Affinity Purification technique

Organs and specific cell types execute specialized functions in multicellular organisms, in large part through customized gene expression signatures. Thus, profiling the transcriptomes of specific cell and tissue types remains an important tool for understanding how cells become specialized. Methodological approaches to detect gene expression differences have utilized samples from whole animals, dissected tissues, and more recently single cells. Despite these advances, there is still a challenge and a need in most laboratories to implement less invasive yet powerful cell-type specific transcriptome profiling methods. Here, we describe the use of the Translating Ribosome Affinity Purification (TRAP) method for C. elegans to detect cell type-specific gene expression patterns at the level of translating mRNAs. In TRAP, a ribosomal protein is fused to a tag (GFP) and is expressed under cell type-specific promoters to mark genetically defined cell types in vivo. Affinity purification of lysates of animals expressing the tag enriches for ribosome-associated mRNAs of the targeted tissue. The purified mRNAs are used for making cDNA libraries subjected to high-throughput sequencing to obtain genome-wide profiles of transcripts from the targeted cell type. The ease of exposing C. elegans to diverse stimuli, coupled with available cell type specific promoters, makes TRAP a useful approach to enable the discovery of molecular components in response to external or genetic perturbations.

Alternative Polyadenylation Directs Tissue-Specific miRNA Targeting in Caenorhabditis elegans Somatic Tissues

mRNA expression dynamics promote and maintain the identity of somatic tissues in living organisms; however, their impact in post-transcriptional gene regulation in these processes is not fully understood. Here, we applied the PAT-Seq approach to systematically isolate, sequence, and map tissue-specific mRNA from five highly studied Caenorhabditis elegans somatic tissues: GABAergic and NMDA neurons, arcade and intestinal valve cells, seam cells, and hypodermal tissues, and studied their mRNA expression dynamics. The integration of these datasets with previously profiled transcriptomes of intestine, pharynx, and body muscle tissues, precisely assigns tissue-specific expression dynamics for 60% of all annotated C. elegans protein-coding genes, providing an important resource for the scientific community. The mapping of 15,956 unique high-quality tissue-specific polyA sites in all eight somatic tissues reveals extensive tissue-specific 3'untranslated region (3'UTR) isoform switching through alternative polyadenylation (APA) . Almost all ubiquitously transcribed genes use APA and harbor miRNA targets in their 3'UTRs, which are commonly lost in a tissue-specific manner, suggesting widespread usage of post-transcriptional gene regulation modulated through APA to fine tune tissue-specific protein expression. Within this pool, the human disease gene C. elegans orthologs rack-1 and tct-1 use APA to switch to shorter 3'UTR isoforms in order to evade miRNA regulation in the body muscle tissue, resulting in increased protein expression needed for proper body muscle function. Our results highlight a major positive regulatory role for APA, allowing genes to counteract miRNA regulation on a tissue-specific basis.

PRIMA: a gene-centered, RNA-to-protein method for mapping RNA-protein interactions

Interactions between RNA binding proteins (RBPs) and mRNAs are critical to post-transcriptional gene regulation. Eukaryotic genomes encode thousands of mRNAs and hundreds of RBPs. However, in contrast to interactions between transcription factors (TFs) and DNA, the interactome between RBPs and RNA has been explored for only a small number of proteins and RNAs. This is largely because the focus has been on using 'protein-centered' (RBP-to-RNA) interaction mapping methods that identify the RNAs with which an individual RBP interacts. While powerful, these methods cannot as of yet be applied to the entire RBPome. Moreover, it may be desirable for a researcher to identify the repertoire of RBPs that can interact with an mRNA of interest-in a 'gene-centered' manner-yet few such techniques are available. Here, we present Protein-RNA Interaction Mapping Assay (PRIMA) with which an RNA 'bait' can be tested versus multiple RBP 'preys' in a single experiment. PRIMA is a translation-based assay that examines interactions in the yeast cytoplasm, the cellular location of mRNA translation. We show that PRIMA can be used with small RNA elements, as well as with full-length Caenorhabditis elegans 3' UTRs. PRIMA faithfully recapitulated numerous well-characterized RNA-RBP interactions and also identified novel interactions, some of which were confirmed in vivo. We envision that PRIMA will provide a complementary tool to expand the depth and scale with which the RNA-RBP interactome can be explored.

Tissue-specific regulation of alternative polyadenylation represses expression of a neuronal ankyrin isoform in C. elegans epidermal development

Differential mRNA polyadenylation plays an important role in shaping the neuronal transcriptome. In C. elegans, several ankyrin isoforms are produced from the unc-44 locus through alternative polyadenylation. Here, we identify a key role for an intronic polyadenylation site (PAS) in temporal- and tissue-specific regulation of UNC-44/ankyrin isoforms. Removing an intronic PAS results in ectopic expression of the neuronal ankyrin isoform in non-neural tissues. This mis-expression underlies epidermal developmental defects in mutants of the conserved tumor suppressor death-associated protein kinase dapk-1 We have previously reported that the use of this intronic PAS depends on the nuclear polyadenylation factor SYDN-1, which inhibits the RNA polymerase II CTD phosphatase SSUP-72. Consistent with this, loss of sydn-1 blocks ectopic expression of neuronal ankyrin and suppresses epidermal morphology defects of dapk-1 These effects of sydn-1 are mediated by ssup-72 autonomously in the epidermis. We also show that a peptidyl-prolyl isomerase PINN-1 antagonizes SYDN-1 in the spatiotemporal control of neuronal ankyrin isoform. Moreover, the nuclear localization of PINN-1 is altered in dapk-1 mutants. Our data reveal that tissue and stage-specific expression of ankyrin isoforms relies on differential activity of positive and negative regulators of alternative polyadenylation.