Identifying small RNAs derived from maternal- and somatic-type rRNAs in zebrafish development

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex-determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting that Chr4R transcriptomics might differ from the rest of the genome. To test this hypothesis, we conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes in the Nadia strain and identified 4 regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brains and livers validated reduced transcripts from Region-2 in somatic cells, but without sex specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. Region-2 lacks protein-coding genes with human orthologs; has zinc finger genes expressed early in zygotic genome activation; has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and a distinct set of repetitive elements. The colocalization of (1) genes silenced in ovaries but not in testes that are (2) expressed in embryos briefly at the onset of zygotic genome activation; (3) maternal-specific genes for translation machinery; (4) maternal-specific spliceosome components; and (5) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a maternal-to-zygotic-transition gene regulatory block.

A maternal-to-zygotic-transition gene block on the zebrafish sex chromosome

Wild zebrafish (Danio rerio) have a ZZ/ZW chromosomal sex determination system with the major sex locus on the right arm of chromosome-4 (Chr4R) near the largest heterochromatic block in the genome, suggesting the hypothesis that the Chr4R transcriptome might be different from the rest of the genome. We conducted an RNA-seq analysis of adult ZW ovaries and ZZ testes and identified four regions of Chr4 with different gene expression profiles. Unique in the genome, protein-coding genes in a 41.7 Mb section (Region-2) were expressed in testis but silent in ovary. The AB lab strain, which lacks sex chromosomes, verified this result, showing that testis-biased gene expression in Region-2 depends on gonad biology, not on sex-determining mechanism. RNA-seq analyses in female and male brain and liver validated few transcripts from Region-2 in somatic cells, but without sex-specificity. Region-2 corresponds to the heterochromatic portion of Chr4R and its content of genes and repetitive elements distinguishes it from the rest of the genome. In Region-2, protein-coding genes lack human orthologs; it has zinc finger genes expressed early in zygotic genome activation; it has maternal 5S rRNA genes, maternal spliceosome genes, a concentration of tRNA genes, and an distinct set of repetitive elements. The colocalization of 1) genes silenced in ovaries but not in testes that are 2) expressed in embryos briefly at the onset of zygotic genome activation; 3) maternal-specific genes for translation machinery; 4) maternal-specific spliceosome components; and 4) adjacent genes encoding miR-430, which mediates maternal transcript degradation, suggest that this is a Maternal-to-Zygotic-Transition Gene Regulatory Block.

Testicular miRNAs and tsRNAs provide insight into gene regulation during overwintering and reproduction of Onychostoma macrolepis

The late overwintering period and breeding period are two important developmental stages of testis in Onychostoma macrolepis. Small non-coding RNAs (sncRNAs) are well-known regulators of biological processes associated with numerous biological processes. This study aimed to elucidate the roles of four sncRNA classes (microRNAs [miRNAs], Piwi-interacting RNAs [piRNAs], tRNA-derived small RNAs [tsRNAs], and rRNA-derived small RNAs [rsRNAs]) across testes in the late overwintering period (in March) and breeding period (in June) by high-throughput sequencing. The testis of O. macrolepis displayed the highest levels of piRNAs and lowest levels of rsRNAs. Compared with miRNAs and tsRNAs in June, tsRNAs in March had a higher abundance, while miRNAs in March had a much lower abundance. Bioinformatics analysis identified 1,362 and 1,340 differentially expressed miRNAs and tsRNAs, respectively. Further analysis showed that miR-200-1, miR-143-1, tRFi-Lys-CTT-1, and tRFi-Glu-CTC-1 could play critical roles during the overwintering and breeding periods. Our findings provided an unprecedented insight to reveal the epigenetic mechanism underlying the overwintering and reproduction process of male O. macrolepis.

A comprehensive overview of ovarian small non-coding RNAs in the late overwintering and breeding periods of Onychostoma macrolepis

The development of the ovary of Onychostoma macrolepis undergoes distinct annual cyclic changes in which small non-coding RNAs (sncRNAs) could play vital roles. In this study, four sncRNA classes in the ovary, including miRNA, piRNAs, tsRNA, and rsRNAs, were systematically profiled by high-throughput sequencing. In adult ovaries of O. macrolepis, 247 miRNAs and 235 tsRNAs were identified as differentially expressing in the late overwintering period (in March) and breeding period (in June). Some up-regulated sncRNAs in March, such as miR-125-1 and tRFi-Lys-CTT-1, could be involved in inhibiting biomolecule metabolism and enhancing stress tolerance during the overwintering period. Compared with the level expression of sncRNAs in March, some sncRNAs were up-regulated in June, such as miR-146-1 and tRFi-Gly-GCC-1, and could be involved in influencing molecular synthesis and metabolism, enhancing oocyte proliferation and maturation, accelerating ovarian development, and increasing fertilization of oocytes by regulating related target mRNAs. The results suggested that sncRNAs in the ovary of Onychostoma macrolepis not only reflect characteristics of the fish's physiology at different developmental periods, but also directly affect ovarian development and oocyte maturation during the breeding period. In conclusion, these results significantly advance our understanding of the roles of sncRNA during overwintering and reproduction periods, and provide a novel perspective for uncovering characteristics of the special overwintering ecology and reproductive physiology of an atypical cavefish.

Zrsr2 and functional U12-dependent spliceosome are necessary for follicular development

ZRSR2 is a splicing factor involved in recognition of 3′-intron splice sites that is frequently mutated in myeloid malignancies and several tumors; however, the role of mutations of Zrsr2 in other tissues has not been analyzed. To explore the biological role of ZRSR2, we generated three Zrsr2 mutant mouse lines. All Zrsr2 mutant lines exhibited blood cell anomalies, and in two lines, oogenesis was blocked at the secondary follicle stage. RNA-seq of Zrsr2^mu secondary follicles showed aberrations in gene expression and showed altered alternative splicing (AS) events involving enrichment of U12-type intron retention (IR), supporting the functional Zrsr2 action in minor spliceosomes. IR events were preferentially associated with centriole replication, protein phosphorylation, and DNA damage checkpoint. Notably, we found alterations in AS events of 50 meiotic genes. These results indicate that ZRSR2 mutations alter splicing mainly in U12-type introns, which may affect peripheral blood cells, and impede oogenesis and female fertility.

•

Zrsr2^mu mice allow us to identify functions of Zrsr2 in vivo

•

Minor splicing factor Zrsr2 is essential for oogenesis and peripheral blood cells

•

Zrsr2 impairment affects the splicing of U12-type intron-containing genes

•

Zrsr2^mu aberrant splicing causes a global alteration of gene expression

Sperm-borne sncRNAs: potential biomarkers for semen fertility?

Semen infertility or sub-fertility, whether in humans or livestock species, remains a major concern for clinicians and technicians involved in reproduction. Indeed, they can cause tragedies in human relationships or have a dramatic overall negative impact on the sustainability of livestock breeding. Understanding and predicting semen fertility issues is therefore crucial and quality control procedures as well as biomarkers have been proposed to ensure sperm fertility. However, their predictive values appeared to be too limited and additional relevant biomarkers are still required to diagnose sub-fertility efficiently. During the last decade, the study of molecular mechanisms involved in spermatogenesis and sperm maturation highlighted the regulatory role of a variety of small non-coding RNAs (sncRNAs) and led to the discovery that sperm sncRNAs comprise both remnants from spermatogenesis and post-testicular sncRNAs acquired through interactions with extracellular vesicles along epididymis. This has led to the hypothesis that sncRNAs may be a source of relevant biomarkers, associated either with sperm functionality or embryo development. This review aims at providing a synthetic overview of the current state of knowledge regarding implication of sncRNA in spermatogenesis defects and their putative roles in sperm maturation and embryo development, as well as exploring their use as fertility biomarkers.

5´XP sRNA-seq: efficient identification of transcripts with and without 5´ phosphorylation reveals evolutionary conserved small RNA

Small RNA (sRNA) sequencing has been critical for our understanding of many cellular processes, including gene regulation. Nonetheless, the varying biochemical properties of sRNA, such as 5´ nucleotide modifications, make many sRNA subspecies incompatible with common protocols for sRNA sequencing. Here we describe 5XP-seq that outlines a novel strategy that captures a more complete picture of sRNA. By tagging 5´P sRNA during library preparation, 5XP-seq combines an open approach that includes all types of 5'-terminal modifications (5´X), with a selective approach for 5-phosphorylated sRNA (5´P). We show that 5XP-seq not only enriches phosphorylated miRNA and piRNA but successfully discriminates these sRNA from all other sRNA species. We further demonstrate the importance of this strategy by successful inter-species validation of sRNAs that would have otherwise failed, including human to insect translation of several tRNA (tRFs) and rRNA (rRFs) fragments. By combining 5´ insensitive library strategies with 5´ sensitive tagging, we have successfully tackled an intrinsic bias in modern sRNA sequencing that will help us reveal the true complexity and the evolutionary significance of the sRNA world.

A New Specific and Sensitive RT-qPCR Method Based on Splinted 5' Ligation for the Quantitative Detection of RNA Species Shorter than microRNAs

Recently, we discovered a new family of unusually short RNAs mapping to 5.8S ribosomal RNA (rRNA) and which we named dodecaRNAs (doRNAs), according to the number of core nucleotides (12 nt) their members contain. To confirm these small RNA-sequencing (RNA-Seq) data, validate the existence of the two overly abundant doRNAs-the minimal core 12-nt doRNA sequence and its + 1-nt variant bearing a 5' Cytosine, C-doRNA-and streamline their analysis, we developed a new specific and sensitive splinted 5' ligation reverse transcription (RT)-quantitative polymerase chain reaction (qPCR) method. This method is based on a splint-assisted ligation of an adapter to the 5' end of doRNAs, followed by RT-qPCR amplification and quantitation. Our optimized protocol, which may discriminate between doRNA, C-doRNA, mutated and precursor sequences, can accurately detect as low as 240 copies and is quantitatively linear over a range of 7 logs. This method provides a unique tool to expand and facilitate studies exploring the molecular and cellular biology of RNA species shorter than microRNAs.

Identification of Abundant and Functional dodecaRNAs (doRNAs) Derived from Ribosomal RNA

Using a modified RNA-sequencing (RNA-seq) approach, we discovered a new family of unusually short RNAs mapping to ribosomal RNA 5.8S, which we named dodecaRNAs (doRNAs), according to the number of core nucleotides (12 nt) their members contain. Using a new quantitative detection method that we developed, we confirmed our RNA-seq data and determined that the minimal core doRNA sequence and its 13-nt variant C-doRNA (doRNA with a 5' Cytosine) are the two most abundant doRNAs, which, together, may outnumber microRNAs. The C-doRNA/doRNA ratio is stable within species but differed between species. doRNA and C-doRNA are mainly cytoplasmic and interact with heterogeneous nuclear ribonucleoproteins (hnRNP) A0, A1 and A2B1, but not Argonaute 2. Reporter gene activity assays suggest that C-doRNA may function as a regulator of Annexin II receptor (AXIIR) expression. doRNAs are differentially expressed in prostate cancer cells/tissues and may control cell migration. These findings suggest that unusually short RNAs may be more abundant and important than previously thought.

MicroRNAs and long non-coding RNAs as novel regulators of ribosome biogenesis

Ribosome biogenesis is the fine-tuned, essential process that generates mature ribosomal subunits and ultimately enables all protein synthesis within a cell. Novel regulators of ribosome biogenesis continue to be discovered in higher eukaryotes. While many known regulatory factors are proteins or small nucleolar ribonucleoproteins, microRNAs (miRNAs), and long non-coding RNAs (lncRNAs) are emerging as a novel modulatory layer controlling ribosome production. Here, we summarize work uncovering non-coding RNAs (ncRNAs) as novel regulators of ribosome biogenesis and highlight their links to diseases of defective ribosome biogenesis. It is still unclear how many miRNAs or lncRNAs are involved in phenotypic or pathological disease outcomes caused by impaired ribosome production, as in the ribosomopathies, or by increased ribosome production, as in cancer. In time, we hypothesize that many more ncRNA regulators of ribosome biogenesis will be discovered, which will be followed by an effort to establish connections between disease pathologies and the molecular mechanisms of this additional layer of ribosome biogenesis control.

Single-cell RNA profiling links ncRNAs to spatiotemporal gene expression during C. elegans embryogenesis

Recent studies show that non-coding RNAs (ncRNAs) can regulate the expression of protein-coding genes and play important roles in mammalian development. Previous studies have revealed that during C. elegans (Caenorhabditis elegans) embryo development, numerous genes in each cell are spatiotemporally regulated, causing the cell to differentiate into distinct cell types and tissues. We ask whether ncRNAs participate in the spatiotemporal regulation of genes in different types of cells and tissues during the embryogenesis of C. elegans. Here, by using marker-free full-length high-depth single-cell RNA sequencing (scRNA-seq) technique, we sequence the whole transcriptomes from 1031 embryonic cells of C. elegans and detect 20,431 protein-coding genes, including 22 cell-type-specific protein-coding markers, and 9843 ncRNAs including 11 cell-type-specific ncRNA markers. We induce a ncRNAs-based clustering strategy as a complementary strategy to the protein-coding gene-based clustering strategy for single-cell classification. We identify 94 ncRNAs that have never been reported to regulate gene expressions, are co-expressed with 1208 protein-coding genes in cell type specific and/or embryo time specific manners. Our findings suggest that these ncRNAs could potentially influence the spatiotemporal expression of the corresponding genes during the embryogenesis of C. elegans.

Roles of Regulatory RNAs in Nutritional Control

Small RNAs (sRNAs), including microRNAs (miRNAs), are noncoding RNA (ncRNA) molecules involved in gene regulation. sRNAs play important roles in development; however, their significance in nutritional control and as metabolic modulators is still emerging. The mechanisms by which diet impacts metabolic genes through miRNAs remain an important area of inquiry. Recent work has established how miRNAs are transported in body fluids often within exosomes, which are small cell-derived vesicles that function in intercellular communication. The abundance of other recently identified ncRNAs and new insights regarding ncRNAs as dietary bioactive compounds could remodel our understanding about how foods impact gene expression. Although controversial, some groups have shown that dietary RNAs from plants and animals (i.e., milk) are functional in consumers. In the future, regulating sRNAs either directly through dietary delivery or indirectly by altered expression of endogenous sRNA may be part of nutritional interventions for regulating metabolism.

Emerging roles of novel small non-coding regulatory RNAs in immunity and cancer

The term small non-coding RNAs (ncRNAs) refers to all those RNAs that even without encoding for a protein, can play important functional roles. Transfer RNA and ribosomal RNA-derived fragments (tRFs and rRFs, respectively) are an emerging class of ncRNAs originally considered as simple degradation products, which though play important roles in stress responses, signalling, or gene expression. They control all levels of gene expression regulating transcription and translation and affecting RNA processing and maturation. They have been linked to pivotal cellular processes such as self-renewal, differentiation, and proliferation. For this reason, mis-regulation of this novel class of ncRNAs can lead to various pathological processes such as neurodegenerative and development diseases, metabolism and immune system disorders, and cancer. In this review, we summarise the classification, biogenesis, and functions of tRFs and rRFs with a special focus on their role in immunity and cancer.

Nuclear Argonaute Piwi Gene Mutation Affects rRNA by Inducing rRNA Fragment Accumulation, Antisense Expression, and Defective Processing in Drosophila Ovaries

Drosophila key nuclear piRNA silencing pathway protein Piwi of the Argonaute family has been classically studied as a factor controlling transposable elements and fertility. Piwi has been shown to concentrate in the nucleolus for reasons largely unknown. Ribosomal RNA is the main component of the nucleolus. In this work the effect of a piwi mutation on rRNA is described. This work led to three important conclusions: A mutation in piwi induces antisense 5S rRNA expression, a processing defect of 2S rRNA orthologous to the 3′-end of eukaryotic 5.8S rRNA, and accumulation of fragments of all five rRNAs in Drosophilamelanogaster ovaries. Hypotheses to explain these phenomena are proposed, possibly involving the interaction of the components of the piRNA pathway with the RNA surveillance machinery.

Heat stress upregulates G-T mismatch binding activities in zebrafish (Danio rerio) embryos preexposed and nonexposed to a sublethal level of cadmium (Cd)

G-T mispair frequently appears in eukaryotic DNA due to the spontaneous deamination of 5-methylcytosine paired with guanine and is therefore an important target for DNA mismatch repair (MMR). Our earlier studies showed the downregulation of G-T binding activities in cadmium (Cd)-exposed (Danio rerio) embryos. Since elevation of water temperature was reported to increase Cd toxicity in zebrafish, this study explored whether heat stress affected zebrafish mismatch binding capacity in the absence or presence of Cd. Heat stress (37 °C for 30 min) induced heat shock protein 70 mRNA expression in embryos at 10 and 24 h post fertilization (hpf). Heat stress weakly upregulated normal G-T sensing machinery and inhibited G-T recognition activity in embryos preexposed to 3 μM Cd for 9 h. Either heat shock or a 23-h Cd treatment alone caused a 1.7-fold stimulation of G-T binding capacity in 24 hpf embryos and heat stress of Cd-preexposed embryos further enhanced G-T binding activity to 2.5 fold of control. Normal and Cd-downregulated loop binding activities in 10 and 24 hpf embryos were almost unreactive to heat shock. Heat stress-upregulated G-T sensing in nonexposed, but not in Cd-preexposed, 24 hpf embryos correlated with stronger gene activities encoding MMR-linked mismatch detecting factors MutS homolog 2 and 6 plus a higher DNA binding activity of the transcription factor Sp1 that regulates msh2/msh6 expression. Our results suggested the importance of heat shock response in facilitating the correction of G-T mismatch in developing zebrafish even under Cd exposure.

Small Non-Coding RNAs Derived From Eukaryotic Ribosomal RNA

The advent of RNA-sequencing (RNA-Seq) technologies has markedly improved our knowledge and expanded the compendium of small non-coding RNAs, most of which derive from the processing of longer RNA precursors. In this review article, we will present a nonexhaustive list of referenced small non-coding RNAs (ncRNAs) derived from eukaryotic ribosomal RNA (rRNA), called rRNA fragments (rRFs). We will focus on the rRFs that are experimentally verified, and discuss their origin, length, structure, biogenesis, association with known regulatory proteins, and potential role(s) as regulator of gene expression. This relatively new class of ncRNAs remained poorly investigated and underappreciated until recently, due mainly to the a priori exclusion of rRNA sequences-because of their overabundance-from RNA-Seq datasets. The situation surrounding rRFs resembles that of microRNAs (miRNAs), which used to be readily discarded from further analyses, for more than five decades, because no one could believe that RNA of such a short length could bear biological significance. As if we had not yet learned our lesson not to restrain our investigative, scientific mind from challenging widely accepted beliefs or dogmas, and from looking for the hidden treasures in the most unexpected places.