Piwi-interacting RNAs (piRNAs) and small interfering RNAs (siRNAs) are distinct classes of small RNAs required for proper germline development. To identify the roles of piRNAs and siRNAs in regulating gene expression in Caenorhabditis elegans, we subjected small RNAs and mRNAs from the gonads of piRNA and siRNA defective mutants to high-throughput sequencing. We show that piRNAs and an abundant class of siRNAs known as WAGO-class 22G-RNAs are required for proper expression of spermatogenic and oogenic genes. WAGO-class 22G-RNAs are also broadly required for transposon silencing, whereas piRNAs are largely dispensable. piRNAs, however, have a critical role in controlling histone gene expression. In the absence of piRNAs, histone mRNAs are misrouted into the nuclear RNAi pathway involving the Argonaute HRDE-1, concurrent with a reduction in the expression of many histone mRNAs. We also show that high-level gene expression in the germline is correlated with high level 22G-RNA production. However, most highly expressed genes produce 22G-RNAs through a distinct pathway that presumably involves the Argonaute CSR-1. In contrast, genes targeted by the WAGO branch of the 22G-RNA pathway are typically poorly expressed and respond unpredictably to loss of 22G-RNAs. Our results point to broad roles for piRNAs and siRNAs in controlling gene expression in the C. elegans germline.

Caenorhabditis elegans contains 25 Argonautes, of which, ALG-1 and ALG-2 are known to primarily interact with miRNAs. ALG-5 belongs to the AGO subfamily of Argonautes that includes ALG-1 and ALG-2, but its role in small RNA pathways is unknown. We analyzed by high-throughput sequencing the small RNAs associated with ALG-5, ALG-1 and ALG-2, as well as changes in mRNA expression in alg-5, alg-1 and alg-2 mutants. We show that ALG-5 defines a distinct branch of the miRNA pathway affecting the expression of genes involved in immunity, defense, and development. In contrast to ALG-1 and ALG-2, which associate with most miRNAs and have general roles throughout development, ALG-5 interacts with only a small subset of miRNAs and is specifically expressed in the germline where it localizes alongside the piRNA and siRNA machinery at P granules. alg-5 is required for optimal fertility and mutations in alg-5 lead to a precocious transition from spermatogenesis to oogenesis. Our results provide a near-comprehensive analysis of miRNA-Argonaute interactions in C. elegans and reveal a new role for miRNAs in the germline.

SUMMARY The germline contains an immortal cell lineage that ensures the faithful transmission of genetic and, in some instances, epigenetic information from one generation to the next. Here, we show that in Caenorhabditis elegans, the small RNA 3′-2′-O-methyltransferase henn-1/HEN1 is required for sustained fertility across generations. In the absence of henn-1, animals become progressively less fertile, becoming sterile after ~30 generations at 25°C. Sterility in henn-1 mutants is accompanied by severe defects in germline proliferation and maintenance. The requirement for henn-1 in transgenerational fertility is likely due to its role in methylating and, thereby, stabilizing Piwi-interacting RNAs (piRNAs). However, despite being essential for piRNA stability in embryos, henn-1 is not required for piRNA stability in adults. Thus, we propose that methylation is important for the role of piRNAs in establishing proper gene silencing during early stages of development but is dispensable for their role in the proliferated germline.

No abstract

Rapid mutation rates are typical of mitochondrial genomes (mtDNAs) in animals, but it is not clear why. The difficulty of obtaining measurements of mtDNA mutation that are not biased by natural selection has stymied efforts to distinguish between competing hypotheses about the causes of high mtDNA mutation rates. Several studies which have measured mtDNA mutations in nematodes have yielded small datasets with conflicting conclusions about the relative abundance of different substitution classes (i.e. the mutation spectrum). We therefore leveraged Duplex Sequencing, a high-fidelity DNA sequencing technique, to characterize de novo mtDNA mutations in Caenorhabditis elegans. This approach detected nearly an order of magnitude more mtDNA mutations than documented in any previous nematode mutation study. Despite an existing extreme AT bias in the C. elegans mtDNA (75.6% AT), we found that a significant majority of mutations increase genomic AT content. Compared to some prior studies in nematodes and other animals, the mutation spectrum reported here contains an abundance of CG→AT transversions, supporting the hypothesis that oxidative damage may be a driver of mtDNA mutations in nematodes. Further, we found an excess of G→T and C→T changes on the coding DNA strand relative to the template strand, consistent with increased exposure to oxidative damage. Analysis of the distribution of mutations across the mtDNA revealed significant variation among protein-coding genes and as well as among neighboring nucleotides. This high-resolution view of mitochondrial mutations in C. elegans highlights the value of this system for understanding relationships among oxidative damage, replication error, and mtDNA mutation.

Differences in tRNA expression have been implicated in a remarkable number of biological processes.There is growing evidence that tRNA genes can play dramatically different roles depending on both expression and post-transcriptional modification, yet sequencing tRNAs to measure abundance and detect modifications remains challenging. Their secondary structure and extensive post-transcriptional modifications interfere with RNA-seq library preparation methods and have limited the utility of highthroughput sequencing technologies. Here, we combine two modifications to standard RNA-seq methods by treating with the demethylating enzyme AlkB and ligating with tRNA-specific adapters in order to sequence tRNAs from four species of flowering plants, a group that has been shown to have some of the most extensive rates of post-transcriptional tRNA modifications. This protocol has the advantage of detecting full-length tRNAs and sequence variants that can be used to infer many post-transcriptional modifications. We used the resulting data to produce a modification index of almost all unique reference tRNAs in Arabidopsis thaliana, which exhibited many anciently conserved similarities with humans but also positions that appear to be "hot spots" for modifications in angiosperm tRNAs. We also found evidence based on northern blot analysis and droplet digital PCR that, even after demethylation treatment, tRNA-seq can produce highly biased estimates of absolute expression levels most likely due to biased reverse transcription. Nevertheless, the generation of full-length tRNA sequences with modification data is still promising for assessing differences in relative tRNA expression across treatments, tissues or subcellular fractions and help elucidate the functional roles of tRNA modifications.

piRNAs are known to silence transposable elements, but not all piRNAs match transposon sequences. Recent studies from Shen et al. (2018) and Zhang et al. (2018) identify rules for piRNA target recognition in Caenorhabditis elegans. Permissive pairing rules allow targeting of essentially all germline mRNAs, while protective mechanisms prevent silencing self-genes.

Rapid mutation rates are typical of mitochondrial genomes (mtDNAs) in animals, but it is not clear why. The difficulty of obtaining measurements of mtDNA mutation that are not biased by natural selection has stymied efforts to distinguish between competing hypotheses about the causes of high mtDNA mutation rates. Several studies which have measured mtDNA mutations in nematodes have yielded small datasets with conflicting conclusions about the relative abundance of different substitution classes (i.e. the mutation spectrum). We therefore leveraged Duplex Sequencing, a high-fidelity DNA sequencing technique, to characterize de novo mtDNA mutations in Caenorhabditis elegans. This approach detected nearly an order of magnitude more mtDNA mutations than documented in any previous nematode mutation study. Despite an existing extreme AT bias in the C. elegans mtDNA (75.6% AT), we found that a significant majority of mutations increase genomic AT content. Compared to some prior studies in nematodes and other animals, the mutation spectrum reported here contains an abundance of CGAT transversions, supporting the hypothesis that oxidative damage may be a driver of mtDNA mutations in nematodes. Further, we found an excess of GT and CT changes on the coding DNA strand relative to the template strand, consistent with increased exposure to oxidative damage. Analysis of the distribution of mutations across the mtDNA revealed significant variation among protein-coding genes and as well as among neighboring nucleotides. This high-resolution view of mitochondrial mutations in C. elegans highlights the value of this system for understanding relationships among oxidative damage, replication error, and mtDNA mutation.