Full-length direct RNA sequencing reveals extensive remodeling of RNA expression, processing and modification in aging Caenorhabditis elegans.

Organismal aging is marked by decline in cellular function and anatomy, ultimately resulting in death. To inform our understanding of the mechanisms underlying this degeneration, we performed standard RNA sequencing and Nanopore direct RNA sequencing over an adult time course in Caenorhabditis elegans. Long reads allowed for identification of hundreds of novel isoforms and age-associated differential isoform accumulation, resulting from alternative splicing and terminal exon choice. Genome-wide analysis reveals a decline in RNA processing fidelity and a rise in inosine and pseudouridine editing events in transcripts from older animals. In this first map of pseudouridine modifications for C. elegans, we find that they largely reside in coding sequences and that the number of genes with this modification increases with age. Collectively, this analysis discovers transcriptomic signatures associated with age and is a valuable resource to understand the many processes that dictate altered gene expression patterns and post-transcriptional regulation in aging.

Expression, not sequence, distinguishes miR-238 from its miR-239ab sister miRNAs in promoting longevity in Caenorhabditis elegans.

MicroRNAs (miRNAs) regulate gene expression by base-pairing to target sequences in messenger RNAs (mRNAs) and recruiting factors that induce translational repression and mRNA decay. In animals, nucleotides 2-8 at the 5' end of the miRNA, called the seed region, are often necessary and sometimes sufficient for functional target interactions. MiRNAs that contain identical seed sequences are grouped into families where individual members have the potential to share targets and act redundantly. A rare exception seemed to be the miR-238/239ab family in Caenorhabditis elegans, as previous work indicated that loss of miR-238 reduced lifespan while deletion of the miR-239ab locus resulted in enhanced longevity and thermal stress resistance. Here, we re-examined these potentially opposing roles using new strains that individually disrupt each miRNA sister. We confirmed that loss of miR-238 is associated with a shortened lifespan but could detect no longevity or stress phenotypes in animals lacking miR-239a or miR-239b, individually or in combination. Additionally, dozens of genes were mis-regulated in miR-238 mutants but almost no gene expression changes were detected in either miR-239a or miR-239b mutants compared to wild type animals. We present evidence that the lack of redundancy between miR-238 and miR-239ab is independent of their sequence differences; miR-239a or miR-239b could substitute for the longevity role of miR-238 when expressed from the miR-238 locus. Altogether, these studies disqualify miR-239ab as negative regulators of aging and demonstrate that expression, not sequence, dictates the specific role of miR-238 in promoting longevity.

New Roles for MicroRNAs in Old Worms.

The use of Caenorhabditis elegans as a model organism in aging research has been integral to our understanding of genes and pathways involved in this process. Several well-conserved signaling pathways that respond to insulin signaling, diet, and assaults to proteostasis have defined roles in controlling lifespan. New evidence shows that microRNAs (miRNAs) play prominent roles in regulating these pathways. In some cases, key aging-related genes have been established as direct targets of specific miRNAs. However, the precise functions of other miRNAs and their protein cofactors in promoting or antagonizing longevity still need to be determined. Here, we highlight recently uncovered roles of miRNAs in common aging pathways, as well as new techniques for the ongoing discovery of miRNA functions in aging C. elegans.

Nuclear and cytoplasmic poly(A) binding proteins (PABPs) favor distinct transcripts and isoforms.

The poly(A)-tail appended to the 3'-end of most eukaryotic transcripts plays a key role in their stability, nuclear transport, and translation. These roles are largely mediated by Poly(A) Binding Proteins (PABPs) that coat poly(A)-tails and interact with various proteins involved in the biogenesis and function of RNA. While it is well-established that the nuclear PABP (PABPN) binds newly synthesized poly(A)-tails and is replaced by the cytoplasmic PABP (PABPC) on transcripts exported to the cytoplasm, the distribution of transcripts for different genes or isoforms of the same gene on these PABPs has not been investigated on a genome-wide scale. Here, we analyzed the identity, splicing status, poly(A)-tail size, and translation status of RNAs co-immunoprecipitated with endogenous PABPN or PABPC in human cells. At steady state, many protein-coding and non-coding RNAs exhibit strong bias for association with PABPN or PABPC. While PABPN-enriched transcripts more often were incompletely spliced and harbored longer poly(A)-tails and PABPC-enriched RNAs had longer half-lives and higher translation efficiency, there are curious outliers. Overall, our study reveals the landscape of RNAs bound by PABPN and PABPC, providing new details that support and advance the current understanding of the roles these proteins play in poly(A)-tail synthesis, maintenance, and function.

Caenorhabditis elegans transposable elements harbor diverse transcription factor DNA-binding sites.

Transposable elements are powerful agents of evolution that can diversify transcriptional programs by distributing transcription factor DNA-binding sites throughout genomes. To investigate the extent that transposable elements provide transcription factor-binding motifs in Caenorhabditis elegans, we determined the genomic positions of DNA-binding motifs for 201 different transcription factors. Surprisingly, we found that almost all examined transcription factors have binding motifs that reside within transposable elements, and all types of transposable elements have at least 1 instance of a transcription factor motif, demonstrating that transposable elements provide previously unappreciated numbers of transcription factor-binding motifs to the C. elegans genome. After determining the occurrence of transcription factor motifs in transposable elements relative to the rest of the genome, we identified DNA-binding motifs for 45 different transcription factors that are greater than 20-fold enriched within transposable elements compared to what would be expected by chance. Consistent with potential functional roles for these transposable element-enriched transcription factor-binding sequences, we determined that all transcription factor motif types found in transposable elements have instances of residing within accessible chromatin sites associated with transcription factor binding. The overwhelming majority of transcription factor-binding motifs located within transposable elements associate with their cognate transcription factors, suggesting extensive binding of transcription factors to sequences within transposable elements. In addition, transposable elements with accessible or transcription factor-bound motifs reside in the putative promoter regions of approximately 12% of all protein-coding genes, providing widespread possibilities for influencing gene expression. This work represents the first comprehensive analysis of transposable element-transcription factor interactions in C. elegans and demonstrates that transposable element-provided transcription factor-binding sites are prevalent in this important model organism.

Recovery from heat shock requires the microRNA pathway in Caenorhabditis elegans.

The heat shock response (HSR) is a highly conserved cellular process that promotes survival during stress. A hallmark of the HSR is the rapid induction of heat shock proteins (HSPs), such as HSP-70, by transcriptional activation. Once the stress is alleviated, HSPs return to near basal levels through incompletely understood mechanisms. Here, we show that the microRNA pathway acts during heat shock recovery in Caenorhabditis elegans. Depletion of the miRNA Argonaute, Argonaute Like Gene 1 (ALG-1), after an episode of heat shock resulted in decreased survival and perdurance of high hsp-70 levels. We present evidence that regulation of hsp-70 is dependent on miR-85 and sequences in the hsp-70 3'UTR that contain target sites for this miRNA. Regulation of hsp-70 by the miRNA pathway was found to be particularly important during recovery from HS, as animals that lacked miR-85 or its target sites in the hsp-70 3'UTR overexpressed HSP-70 and exhibited reduced viability. In summary, our findings show that down-regulation of hsp-70 by miR-85 after HS promotes survival, highlighting a previously unappreciated role for the miRNA pathway during recovery from stress.

Diversification of the Caenorhabditis heat shock response by Helitron transposable elements.

Heat Shock Factor 1 (HSF-1) is a key regulator of the heat shock response (HSR). Upon heat shock, HSF-1 binds well-conserved motifs, called Heat Shock Elements (HSEs), and drives expression of genes important for cellular protection during this stress. Remarkably, we found that substantial numbers of HSEs in multiple Caenorhabditis species reside within Helitrons, a type of DNA transposon. Consistent with Helitron-embedded HSEs being functional, upon heat shock they display increased HSF-1 and RNA polymerase II occupancy and up-regulation of nearby genes in C. elegans. Interestingly, we found that different genes appear to be incorporated into the HSR by species-specific Helitron insertions in C. elegans and C. briggsae and by strain-specific insertions among different wild isolates of C. elegans. Our studies uncover previously unidentified targets of HSF-1 and show that Helitron insertions are responsible for rewiring and diversifying the Caenorhabditis HSR.

Remodeling of the Caenorhabditis elegans non-coding RNA transcriptome by heat shock.

Elevated temperatures activate a heat shock response (HSR) to protect cells from the pathological effects of protein mis-folding, cellular mis-organization, organelle dysfunction and altered membrane fluidity. This response includes activation of the conserved transcription factor heat shock factor 1 (HSF-1), which binds heat shock elements (HSEs) in the promoters of genes induced by heat shock (HS). The upregulation of protein-coding genes (PCGs), such as heat shock proteins and cytoskeletal regulators, is critical for cellular survival during elevated temperatures. While the transcriptional response of PCGs to HS has been comprehensively analyzed in a variety of organisms, the effect of this stress on the expression of non-coding RNAs (ncRNAs) has not been systematically examined. Here we show that in Caenorhabditis elegans HS induces up- and downregulation of specific ncRNAs from multiple classes, including miRNA, piRNA, lincRNA, pseudogene and repeat elements. Moreover, some ncRNA genes appear to be direct targets of the HSR, as they contain HSF-1 bound HSEs in their promoters and their expression is regulated by this factor during HS. These results demonstrate that multiple ncRNA genes respond to HS, some as direct HSF-1 targets, providing new candidates that may contribute to organismal survival during this stress.

Todos Santos small RNA symposium.

Worm biologists from the United States, Canada, and the United Kingdom gathered at the Colorado State University Todos Santos Center in Baja California Sur, Mexico, April 3-5, 2019 for the Todos Santos Small RNA Symposium. Meeting participants, many of whom were still recovering from the bomb cyclone that struck a large swath of North America just days earlier, were greeted by the warmth and sunshine that is nearly ubiquitous in the sleepy seaside town of Todos Santos. With only 24 speakers, the meeting had the sort of laid-back vibe you might expect amongst the palm trees and ocean breeze of the Pacific coast of Mexico. The meeting started with tracing the laboratory lineages of participants. Not surprisingly, the most common parental lineages represented at the meeting were Dr. Craig Mello, Dr. Gary Ruvkun, and Dr. Victor Ambros, whom, together with Dr. Andy Fire and Dr. David Baulcombe, pioneered the small RNA field. In sad irony, on the closing day of the meeting, participants were met with the news of Dr. Sydney Brenner's passing. By establishing the worm, Caenorhabditis elegans, as a model system Dr. Brenner paved the way for much of the research discussed here.

RNA interference may result in unexpected phenotypes in Caenorhabditis elegans.

RNA interference (RNAi) is a valuable technique to determine gene function. In Caenorhabditis elegans, RNAi can be achieved by feeding worms bacteria carrying a plasmid expressing double-stranded RNA (dsRNA) targeting a gene of interest. The most commonly used plasmid vector for this purpose is L4440. However, it has been noticed that sequences within L4440 may elicit unspecific effects. Here, we provide a comprehensive characterization of these effects and their mechanisms and describe new unexpected phenotypes uncovered by the administration of unspecific exogenous dsRNA. An example involves dsRNA produced by the multiple cloning site (MCS) of L4440, which shares complementary sequences with some widely used reporter vectors and induces partial transgene silencing via the canonical and antiviral RNAi pathway. Going beyond transgene silencing, we found that the reduced embryonic viability of mir-35-41(gk262) mutants is partially reversed by exogenous dsRNA via a mechanism that involves canonical RNAi. These results indicate cross-regulation between different small RNA pathways in C. elegans to regulate embryonic viability. Recognition of the possible unspecific effects elicited by RNAi vectors is important for rigorous interpretation of results from RNAi-based experiments.

miRNA Targeting: Growing beyond the Seed.

miRNAs are small RNAs that guide Argonaute proteins to specific target mRNAs to repress their translation and stability. Canonically, miRNA targeting is reliant on base pairing of the seed region, nucleotides 2-7, of the miRNA to sites in mRNA 3' untranslated regions. Recently, the 3' half of the miRNA has gained attention for newly appreciated roles in regulating target specificity and regulation. In addition, the extent of pairing to the miRNA 3' end can influence the stability of the miRNA itself. These findings highlight the importance of sequences beyond the seed in controlling the function and existence of miRNAs.

Tales of Detailed Poly(A) Tails.

Poly(A) tails are non-templated additions of adenosines at the 3' ends of most eukaryotic mRNAs. In the nucleus, these RNAs are co-transcriptionally cleaved at a poly(A) site and then polyadenylated before being exported to the cytoplasm. In the cytoplasm, poly(A) tails play pivotal roles in the translation and stability of the mRNA. One challenge in studying poly(A) tails is that they are difficult to sequence and accurately measure. However, recent advances in sequencing technology, computational algorithms, and other assays have enabled a more detailed look at poly(A) tail length genome-wide throughout many developmental stages and organisms. With the help of these advances, our understanding of poly(A) tail length has evolved over the past 5 years with the recognition that highly expressed genes can have short poly(A) tails and the elucidation of the seemingly contradictory roles for poly(A)-binding protein (PABP) in facilitating both protection and deadenylation.

Detection of microRNA-Target Interactions by Chimera PCR (ChimP).

MicroRNAs (miRNAs) regulate gene expression by directing Argonaute proteins to target RNAs, which usually results in destabilization and translational inhibition of the target RNA. The prediction of animal miRNA target sites has remained a challenge due to the ability of miRNAs to bind target RNAs through imperfect base pairing. Recently, several labs have established methods to produce biochemical evidence of miRNA-target interactions by generating chimeric reads where the miRNA is ligated to its target RNA. Despite the insights that can be gained from chimera producing methods, the current approaches are inefficient, labor intensive and require computational expertise. Here we describe a method, called Chimera PCR (ChimP), for the validation or testing of specific miRNA-target interactions. This method allows for focused experiments to analyze miRNA targeting in a variety of conditions.

Opposing roles of microRNA Argonautes during Caenorhabditis elegans aging.

Argonaute (AGO) proteins partner with microRNAs (miRNAs) to target specific genes for post-transcriptional regulation. During larval development in Caenorhabditis elegans, Argonaute-Like Gene 1 (ALG-1) is the primary mediator of the miRNA pathway, while the related ALG-2 protein is largely dispensable. Here we show that in adult C. elegans these AGOs are differentially expressed and, surprisingly, work in opposition to each other; alg-1 promotes longevity, whereas alg-2 restricts lifespan. Transcriptional profiling of adult animals revealed that distinct miRNAs and largely non-overlapping sets of protein-coding genes are misregulated in alg-1 and alg-2 mutants. Interestingly, many of the differentially expressed genes are downstream targets of the Insulin/ IGF-1 Signaling (IIS) pathway, which controls lifespan by regulating the activity of the DAF-16/ FOXO transcription factor. Consistent with this observation, we show that daf-16 is required for the extended lifespan of alg-2 mutants. Furthermore, the long lifespan of daf-2 insulin receptor mutants, which depends on daf-16, is strongly reduced in animals lacking alg-1 activity. This work establishes an important role for AGO-mediated gene regulation in aging C. elegans and illustrates that the activity of homologous genes can switch from complementary to antagonistic, depending on the life stage.

A rADAR defense against RNAi.

Adenosine deaminases that act on RNA (ADARs) convert adenosines (A) to inosines (I) in stretches of dsRNA. The biological purpose of these editing events for the vast majority of ADAR substrates is largely unknown. In this issue of Genes & Development, Reich and colleagues (pp. 271-282) demonstrate that in Caenorhabditis elegans, A-to-I editing in double-stranded regions of protein-coding transcripts protects these RNAs from targeting by the RNAi pathway. Disruption of this safeguard through loss of ADAR activity coupled with enhanced RNAi results in developmental abnormalities and profound changes in gene expression that suggest aberrant induction of an antiviral response. Thus, editing of cellular dsRNA by ADAR helps prevent host RNA silencing and inadvertent antiviral activity.

Short poly(A) tails are a conserved feature of highly expressed genes.

Poly(A) tails are important elements in mRNA translation and stability, although recent genome-wide studies have concluded that poly(A) tail length is generally not associated with translational efficiency in nonembryonic cells. To investigate whether poly(A) tail size might be coupled to gene expression in an intact organism, we used an adapted TAIL-seq protocol to measure poly(A) tails in Caenorhabditis elegans. Surprisingly, we found that well-expressed transcripts contain relatively short, well-defined tails. This attribute appears to be dependent on translational efficiency, as transcripts enriched for optimal codons and ribosome association had the shortest tail sizes, whereas noncoding RNAs retained long tails. Across eukaryotes, short tails were a feature of abundant and well-translated mRNAs. This seems to contradict the dogma that deadenylation induces translational inhibition and mRNA decay and suggests that well-expressed mRNAs accumulate with pruned tails that accommodate a minimal number of poly(A)-binding proteins, which may be ideal for protective and translational functions.

A sense-able microRNA.

In this issue of Genes & Development, Drexel and colleagues (pp. 2042-2047) present a beautiful example of how microRNAs (miRNAs) can regulate tissue-specific gene expression in a biologically relevant setting. They found that miR-791 is expressed in only three types of carbon dioxide (CO2)-sensing neurons in Caenorhabditis elegans, and its primary function there seems to be repression of two target genes that interfere with the behavioral response to CO2 Interestingly, these two targets are broadly expressed across other tissues. Thus, restricted miRNA expression can lead to target repression in select tissues to promote distinct cellular physiologies.

Pairing beyond the Seed Supports MicroRNA Targeting Specificity.

To identify endogenous miRNA-target sites, we isolated AGO-bound RNAs from Caenorhabditis elegans by individual-nucleotide resolution crosslinking immunoprecipitation (iCLIP), which fortuitously also produced miRNA-target chimeric reads. Through the analysis of thousands of reproducible chimeras, pairing to the miRNA seed emerged as the predominant motif associated with functional interactions. Unexpectedly, we discovered that additional pairing to 3' sequences is prevalent in the majority of target sites and leads to specific targeting by members of miRNA families. By editing an endogenous target site, we demonstrate that 3' pairing determines targeting by specific miRNA family members and that seed pairing is not always sufficient for functional target interactions. Finally, we present a simplified method, chimera PCR (ChimP), for the detection of specific miRNA-target interactions. Overall, our analysis revealed that sequences in the 5' as well as the 3' regions of a miRNA provide the information necessary for stable and specific miRNA-target interactions in vivo.

A tale of two sequences: microRNA-target chimeric reads.

In animals, a functional interaction between a microRNA (miRNA) and its target RNA requires only partial base pairing. The limited number of base pair interactions required for miRNA targeting provides miRNAs with broad regulatory potential and also makes target prediction challenging. Computational approaches to target prediction have focused on identifying miRNA target sites based on known sequence features that are important for canonical targeting and may miss non-canonical targets. Current state-of-the-art experimental approaches, such as CLIP-seq (cross-linking immunoprecipitation with sequencing), PAR-CLIP (photoactivatable-ribonucleoside-enhanced CLIP), and iCLIP (individual-nucleotide resolution CLIP), require inference of which miRNA is bound at each site. Recently, the development of methods to ligate miRNAs to their target RNAs during the preparation of sequencing libraries has provided a new tool for the identification of miRNA target sites. The chimeric, or hybrid, miRNA-target reads that are produced by these methods unambiguously identify the miRNA bound at a specific target site. The information provided by these chimeric reads has revealed extensive non-canonical interactions between miRNAs and their target mRNAs, and identified many novel interactions between miRNAs and noncoding RNAs.

Splicing remodels the let-7 primary microRNA to facilitate Drosha processing in Caenorhabditis elegans.

MicroRNAs (miRNAs) are a class of small noncoding RNAs that use partial base-pairing to recognize and regulate the expression of messenger RNAs (mRNAs). Mature miRNAs arise from longer primary transcripts (pri-miRNAs) that are processed to a shorter hairpin precursor miRNA (pre-miRNA) by the Microprocessor complex. In Caenorhabditis elegans the primary let-7 (pri-let-7) transcript undergoes trans-splicing, where pri-let-7 is cleaved at a 3' splice site and the splice-leader-1 (SL1) sequence is appended at the 5' end. Here we investigate the role of this splicing event in the biogenesis of let-7 miRNA. We hypothesized that splicing changes the secondary structure of the pri-let-7 transcript, creating a more favorable substrate for recognition by the Microprocessor. Supporting this idea, we detected conspicuous structural differences between unspliced and SL1-spliced pri-let-7 transcripts using in vitro ribonuclease (RNase) assays. Through the generation of transgenic worm strains, we found that the RNA secondary structure produced by splicing, as opposed to the act of splicing itself, optimizes processing of pri-let-7 by the Microprocessor in vivo. We also observed that the endogenous spliced, but not the unspliced, pri-let-7 transcripts bind to the Microprocessor and accumulate upon its depletion. We conclude that splicing is a key step in generating pri-let-7 transcripts with a structure that enables downstream processing events to produce appropriate levels of mature let-7.