Summary Endogenous small RNAs (endo-siRNAs) interact with Argonaute (AGO) proteins to mediate sequence-specific regulation of diverse biological processes. Here, we combine deep-sequencing and genetic approaches to explore the biogenesis and function of endo-siRNAs in C. elegans. We describe conditional alleles of the dicer-related helicase, drh-3, that abrogate both RNA interference and the biogenesis of endo-siRNAs, called 22G-RNAs. DRH-3 is a core component of RNA-dependent RNA polymerase (RdRP) complexes essential for several distinct 22G-RNA systems. We show that in the germ-line, one system is dependent on worm-specific AGOs, including WAGO-1, which localizes to germ-line nuage structures called P-granules. WAGO-1 silences certain genes, transposons, pseudogenes and cryptic loci. Finally, we demonstrate that components of the nonsense-mediated decay pathway function in at least one WAGO-mediated surveillance pathway. These findings broaden our understanding of the biogenesis and diversity of 22G-RNAs and suggest novel regulatory functions for small RNAs.

In metazoans, Piwi-related Argonaute proteins have been linked to germline maintenance, and to a class of germline-enriched small RNAs termed piRNAs. Here we show that an abundant class of 21 nucleotide small RNAs (21U-RNAs) are expressed in the C. elegans germline, interact with the C. elegans Piwi family member PRG-1, and depend on PRG-1 activity for their accumulation. The PRG-1 protein is expressed throughout development and localizes to nuage-like structures called P granules. Although 21U-RNA loci share a conserved upstream sequence motif, the mature 21U-RNAs are not conserved and, with few exceptions, fail to exhibit complementarity or evidence for direct regulation of other expressed sequences. Our findings demonstrate that 21U-RNAs are the piRNAs of C. elegans and link this class of small RNAs and their associated Piwi Argonaute to the maintenance of temperature-dependent fertility.

Cell cycle progression is dependent on two major waves of gene expression. Early cell cycle gene expression occurs during G1/S to generate factors required for DNA replication, while late cell cycle gene expression begins during G2 to prepare for mitosis. Here we demonstrate that the MuvB complex-comprised of LIN9, LIN37, LIN52, LIN54, and RBBP4-serves an essential role in three distinct transcription complexes to regulate cell cycle gene expression. The MuvB complex, together with the Rb-like protein p130, E2F4, and DP1, forms the DREAM complex during quiescence and represses expression of both early and late genes. Upon cell cycle entry, the MuvB complex dissociates from p130/DREAM, binds to B-Myb, and reassociates with the promoters of late genes during S phase. MuvB and B-Myb are required for the subsequent recruitment of FoxM1 to late gene promoters during G2. The MuvB complex remains bound to FoxM1 during peak late cell cycle gene expression, while B-Myb binding is lost when it undergoes phosphorylation-dependent, proteasome-mediated degradation during late S phase. Our results reveal a novel role for the MuvB complex in recruiting B-Myb and FoxM1 to promote late cell cycle gene expression and in regulating cell cycle gene expression from quiescence through mitosis.

Diminished ovarian reserve (DOR) is a challenging diagnosis of infertility, as there are currently no tests to predict who may become affected with this condition, or at what age. We designed the present study to compare the gene expression profile of membrana granulosa cells from young women affected with DOR with those from egg donors of similar age and to determine if distinct genetic patterns could be identified to provide insight into the etiology of DOR. Young women with DOR were identified based on FSH level in conjunction with poor follicular development during an IVF cycle (n = 13). Egg donors with normal ovarian reserve (NOR) comprised the control group (n = 13). Granulosa cells were collected following retrieval, RNA was extracted and microarray analysis was conducted to evaluate genetic differences between the groups. Confirmatory studies were undertaken with quantitative RT-PCR (qRT-PCR). Multiple significant differences in gene expression were observed between the DOR patients and egg donors. Two genes linked with ovarian function, anti-Mullerian hormone (AMH) and luteinizing hormone receptor (LHCGR), were further analyzed with qRT-PCR in all patients. The average expression of AMH was significantly higher in egg donors (adjusted P-value = 0.01), and the average expression of LHCGR was significantly higher in DOR patients (adjusted P-value = 0.005). Expression levels for four additional genes, progesterone receptor membrane component 2 (PGRMC2), prostaglandin E receptor 3 (subtype EP3) (PTGER3), steroidogenic acute regulatory protein (StAR), and StAR-related lipid transfer domain containing 4 (StarD4), were validated in a group consisting of five NOR and five DOR patients. We conclude that gene expression analysis has substantial potential to determine which young women may be affected with DOR. More importantly, our analysis suggests that DOR patients fall into two distinct subgroups based on gene expression profiles, indicating that different mechanisms may be involved during development of this pathology.

A method to deduce RNA secondary structure on the basis of data from microarrays of 2'-O-methyl RNA 9-mers immobilized in agarose film on glass slides is tested with a 249 nucleotide RNA from the 3' end of the R2 retrotransposon from Bombyx mori. Various algorithms incorporating binding data and free-energy minimization calculations were compared for interpreting the data to provide possible secondary structures. Two different methods give structures with 100 and 87% of the base pairs determined by sequence comparison. In contrast, structures predicted by free-energy minimization alone by Mfold and RNAstructure contain 52 and 72% of the known base pairs, respectively. This combination of high throughput microarray techniques with algorithms using free-energy calculations has potential to allow for fast determination of RNA secondary structure. It should also facilitate the design of antisense and siRNA oligonucleotides.

Binding of small oligonucleotides to the periphery of folded RNA can provide insight into the secondary structure of complex RNA in solution. To discriminate between bound and unbound fluorescein-labeled 29-O-methyl RNA probes, we use ionically coated gold nanoparticles to selectively adsorb unbound probes and quench their fluorescence. The target is the 39 untranslated region of Bombyx mori R2 RNA. Fluorescence indicates that R2 sequences complementary to some of the probes are accessible for binding in the three-dimensional structure. Hybridization occurs under homogeneous conditions in the absence of the gold nanoparticles so that steric issues associated with chip-based assays are avoided. The assay is compatible with well plate formats, takes less than 5 min, and requires only 2 pmol or less of unlabeled target RNA per probe sequence tested.

Next-generation sequencing is increasingly employed in biomedical investigations. Strong concordance between microarray and mRNA-seq levels has been reported in high quality specimens but information is lacking on formalin-fixed, paraffin-embedded (FFPE) tissues, and particularly for microRNA (miRNA) analysis. We conducted a preliminary examination of the concordance between miRNA-seq and cDNA-mediated annealing, selection, extension, and ligation (DASL) miRNA assays. Quantitative agreement between platforms is moderate (Spearman correlation 0.514–0.596) and there is discordance of detection calls on a subset of miRNAs. Quantitative PCR (q-RT-PCR) performed for several discordant miRNAs confirmed the presence of most sequences detected by miRNA-seq but not by DASL but also that miRNA-seq did not detect some sequences, which DASL confidently detected. Our results suggest that miRNA-seq is specific, with few false positive calls, but it may not detect certain abundant miRNAs in FFPE tissue. Further work is necessary to fully address these issues that are pertinent for translational research.

Hyriopsis cumingii is an important species for freshwater pearl cultivation in China. In terms of pearl production, males have larger pearls and better glossiness than females, but there are few reports focusing on the sex of H. cumingii. In this study, six mRNA and six microRNA (miRNA) libraries were prepared from ovaries and testes. Additionally, 28,502 differentially expressed genes (DEGs) and 32 differentially expressed miRNAs (DEMs) were identified. Compared with testis, 14,360 mRNAs and 20 miRNAs were up-regulated in ovary, 14,142 mRNAs and 12 miRNAs were down-regulated. In DEGs, the known genes related to sex determinism and/or differentiation were also identified, such as DMRT1, SOX9, SF1 for males, FOXL2 for females, and other potentially significant candidate genes. Three sex-related pathways have also been identified, which are Wnt, Notch, and TGF-beta. In 32 DEMs, the three miRNAs (miR-9-5p, miR-92, miR-184) were paid more attention, they predicted 28 target genes, which may also be candidates for sex-related miRNAs and genes. Differential miRNAs target genes analysis reveals the pathway associated with oocyte meiosis and spermatogenesis. Overall, the findings of the study provide significant insights to enhance our understanding of sex differentiation and/or sex determination mechanisms for H. cumingii.