A highly specific microRNA-mediated mechanism silences LTR retrotransposons of strawberry

MicroRNAs play regulatory roles in genomic balance

Classic genetics studies found that genomic imbalance caused by changing the dosage of part of the genome (aneuploidy) has more detrimental effects than altering the dosage of the whole genome (ploidy). Previous analysis revealed global modulation of gene expression triggered by aneuploidy across various species, including maize (Zea mays), Arabidopsis, yeast, mammals, etc. Plant microRNAs (miRNAs) are a class of 20- to 24-nt endogenous small noncoding RNAs that carry out post-transcriptional gene expression regulation. That miRNAs and their putative targets are preferentially retained as duplicates after whole-genome duplication, as are many transcription factors and signaling components, indicates miRNAs are likely to be dosage-sensitive and potentially involved in genomic balance networks. This review addresses the following questions regarding the role of miRNAs in genomic imbalance. (1) How do aneuploidy and polyploidy impact the expression of miRNAs? (2) Do miRNAs play a regulatory role in modulating the expression of their targets under genomic imbalance?

Retrotransposons: How the continuous evolutionary front shapes plant genomes for response to heat stress

Long terminal repeat retrotransposons (LTR retrotransposons) are the most abundant group of mobile genetic elements in eukaryotic genomes and are essential in organizing genomic architecture and phenotypic variations. The diverse families of retrotransposons are related to retroviruses. As retrotransposable elements are dispersed and ubiquitous, their "copy-out and paste-in" life cycle of replicative transposition leads to new genome insertions without the excision of the original element. The overall structure of retrotransposons and the domains responsible for the various phases of their replication is highly conserved in all eukaryotes. The two major superfamilies of LTR retrotransposons, Ty1/Copia and Ty3/Gypsy, are distinguished and dispersed across the chromosomes of higher plants. Members of these superfamilies can increase in copy number and are often activated by various biotic and abiotic stresses due to retrotransposition bursts. LTR retrotransposons are important drivers of species diversity and exhibit great variety in structure, size, and mechanisms of transposition, making them important putative actors in genome evolution. Additionally, LTR retrotransposons influence the gene expression patterns of adjacent genes by modulating potential small interfering RNA (siRNA) and RNA-directed DNA methylation (RdDM) pathways. Furthermore, comparative and evolutionary analysis of the most important crop genome sequences and advanced technologies have elucidated the epigenetics and structural and functional modifications driven by LTR retrotransposon during speciation. However, mechanistic insights into LTR retrotransposons remain obscure in plant development due to a lack of advancement in high throughput technologies. In this review, we focus on the key role of LTR retrotransposons response in plants during heat stress, the role of centromeric LTR retrotransposons, and the role of LTR retrotransposon markers in genome expression and evolution.

The miRNome function transitions from regulating developmental genes to transposable elements during pollen maturation

Animal and plant microRNAs (miRNAs) are essential for the spatio-temporal regulation of development. Together with this role, plant miRNAs have been proposed to target transposable elements (TEs) and stimulate the production of epigenetically active small interfering RNAs. This activity is evident in the plant male gamete containing structure, the male gametophyte or pollen grain. How the dual role of plant miRNAs, regulating both genes and TEs, is integrated during pollen development and which mRNAs are regulated by miRNAs in this cell type at a genome-wide scale are unknown. Here, we provide a detailed analysis of miRNA dynamics and activity during pollen development in Arabidopsis thaliana using small RNA and degradome parallel analysis of RNA end high-throughput sequencing. Furthermore, we uncover miRNAs loaded into the two main active Argonaute (AGO) proteins in the uninuclear and mature pollen grain, AGO1 and AGO5. Our results indicate that the developmental progression from microspore to mature pollen grain is characterized by a transition from miRNAs targeting developmental genes to miRNAs regulating TE activity.

Innate, translation-dependent silencing of an invasive transposon in Arabidopsis

Co‐evolution between hosts’ and parasites’ genomes shapes diverse pathways of acquired immunity based on silencing small (s)RNAs. In plants, sRNAs cause heterochromatinization, sequence degeneration, and, ultimately, loss of autonomy of most transposable elements (TEs). Recognition of newly invasive plant TEs, by contrast, involves an innate antiviral‐like silencing response. To investigate this response’s activation, we studied the single‐copy element EVADÉ (EVD), one of few representatives of the large Ty1/Copia family able to proliferate in Arabidopsis when epigenetically reactivated. In Ty1/Copia elements, a short subgenomic mRNA (shGAG) provides the necessary excess of structural GAG protein over the catalytic components encoded by the full‐length genomic flGAG‐POL. We show here that the predominant cytosolic distribution of shGAG strongly favors its translation over mostly nuclear flGAG‐POL. During this process, an unusually intense ribosomal stalling event coincides with mRNA breakage yielding unconventional 5’OH RNA fragments that evade RNA quality control. The starting point of sRNA production by RNA‐DEPENDENT‐RNA‐POLYMERASE‐6 (RDR6), exclusively on shGAG, occurs precisely at this breakage point. This hitherto‐unrecognized “translation‐dependent silencing” (TdS) is independent of codon usage or GC content and is not observed on TE remnants populating the Arabidopsis genome, consistent with their poor association, if any, with polysomes. We propose that TdS forms a primal defense against EVD de novo invasions that underlies its associated sRNA pattern.

Plant transfer RNA-derived fragments: Biogenesis and functions

Processing of mature transfer RNAs (tRNAs) produces complex populations of tRNA-derived fragments (tRFs). Emerging evidence shows that tRFs have important functions in bacteria, animals, and plants. Here, we review recent advances in understanding plant tRFs, focusing on their biological and cellular functions, such as regulating stress responses, mediating plant-pathogen interactions, and modulating post-transcriptional gene silencing and translation. We also review sequencing strategies and bioinformatics resources for studying tRFs in plants. Finally, we discuss future directions for plant tRF research, which will expand our knowledge of plant non-coding RNAs.

Plant Small RNA World Growing Bigger: tRNA-Derived Fragments, Longstanding Players in Regulatory Processes

In the past 2 decades, the discovery of a new class of small RNAs, known as tRNA-derived fragments (tRFs), shed light on a new layer of regulation implicated in many biological processes. tRFs originate from mature tRNAs and are classified according to the tRNA regions that they derive from, namely 3′tRF, 5′tRF, and tRF-halves. Additionally, another tRF subgroup deriving from tRNA precursors has been reported, the 3′U tRFs. tRF length ranges from 17 to 26 nt for the 3′and 5′tRFs, and from 30 to 40 nt for tRF-halves. tRF biogenesis is still not yet elucidated, although there is strong evidence that Dicer (and DICER-LIKE) proteins, as well as other RNases such as Angiogenin in mammal and RNS proteins family in plants, are responsible for processing specific tRFs. In plants, the abundance of those molecules varies among tissues, developmental stages, and environmental conditions. More recently, several studies have contributed to elucidate the role that these intriguing molecules may play in all organisms. Among the recent discoveries, tRFs were found to be involved in distinctive regulatory layers, such as transcription and translation regulation, RNA degradation, ribosome biogenesis, stress response, regulatory signaling in plant nodulation, and genome protection against transposable elements. Although tRF biology is still poorly understood, the field has blossomed in the past few years, and this review summarizes the most recent developments in the tRF field in plants.

Phaseolus vulgaris MIR1511 genotypic variations differentially regulate plant tolerance to aluminum toxicity

The common-bean (Phaseolus vulgaris), a widely consumed legume, originated in Mesoamerica and expanded to South America, resulting in the development of two geographically distinct gene pools. Poor soil condition, including metal toxicity, are often constraints to common-bean crop production. Several P. vulgaris miRNAs, including miR1511, respond to metal toxicity. The MIR1511 gene sequence from the two P. vulgaris model sequenced genotypes revealed that, as opposed to BAT93 (Mesoamerican), the G19833 (Andean) accession displays a 58-bp deletion, comprising the mature and star miR1511 sequences. Genotyping-By-Sequencing data analysis from 87 non-admixed Phaseolus genotypes, comprising different Phaseolus species and P. vulgaris populations, revealed that all the P. vulgaris Andean genotypes and part of the Mesoamerican (MW1) genotypes analyzed displayed a truncated MIR1511 gene. The geographic origin of genotypes with a complete versus truncated MIR1511 showed a distinct distribution. The P. vulgaris ALS3 (Aluminum Sensitive Protein 3) gene, known to be important for aluminum detoxification in several plants, was experimentally validated as the miR1511 target. Roots from BAT93 plants showed decreased miR1511 and increased ALS3 transcript levels at early stages under aluminum toxicity (AlT), while G19833 plants, lacking mature miR1511, showed higher and earlier ALS3 response. Root architecture analyses evidenced higher tolerance of G19833 plants to AlT. However, G19833 plants engineered for miR1511 overexpression showed lower ALS3 transcript level and increased sensitivity to AlT. Absence of miR1511 in Andean genotypes, resulting in a diminished ALS3 transcript degradation, appears to be an evolutionary advantage to high Al levels in soils with increased drought conditions.

Analysis of tRNA-derived RNA fragments (tRFs) in Cryptococcus spp.: RNAi-independent generation and possible compensatory effects in a RNAi-deficient genotype

Small RNAs (sRNAs) are key factors in the regulation of gene expression. Recently, a new class of regulatory sRNAs derived from tRNAs was described, the tRNA-derived RNA fragments (tRFs). Such RNAs range in length from 14 to 30 nucleotides and are produced from both mature and primary tRNA transcripts, with very specific cleavage sites along the tRNA sequence. Although several mechanisms have been proposed for how tRFs mediate regulation of gene expression, the exact mechanism of tRF biogenesis and its dependency upon the RNAi pathway remain unclear. Cryptococcus gattii and Cryptococcus neoformans are basidiomycetous yeasts and important human pathogens. While C. neoformans is RNAi proficient, C. gattii VGII has lost essential RNAi genes. Here, we sought to identify the tRF production profile in C. gattii VGII and C. neoformans in order to assess the RNAi-dependency of tRF production in these fungal species. We developed a RNA-sequencing-based tRF prediction workflow designed to improve the currently available prediction tools. Using this methodology, we were able to identify tRFs in both organisms. Despite the loss of the RNAi pathway, C. gattii VGII displayed a number of identified tRFs that did not differ significantly from those observed in C. neoformans. The analysis of predicted tRF targets revealed that a higher number of targets was found for C. gattii VGII tRFs compared to C. neoformans tRFs. These results support the idea that tRFs are at least partially independent of the canonical RNAi machinery, raising questions about possible compensatory roles of alternative regulatory RNAs in the absence of a functional RNAi pathway.

Small RNA Function in Plants: From Chromatin to the Next Generation

Small RNA molecules can target a particular virus, gene, or transposable element (TE) with a high degree of specificity. Their ability to move from cell to cell and recognize targets in trans also allows building networks capable of regulating a large number of related targets at once. In the case of epigenetic silencing, small RNA may use the widespread distribution of TEs in eukaryotic genomes to coordinate many loci across developmental and generational time. Here, we discuss the intriguing role of plant small RNA in targeting transposons and repeats in pollen and seeds. Epigenetic reprogramming in the germline and early seed development provides a mechanism to control genome dosage, imprinted gene expression, and incompatible hybridizations via the "triploid block."

Mating-Type-Specific Ribosomal Proteins Control Aspects of Sexual Reproduction in Cryptococcus neoformans

The MAT locus of Cryptococcus neoformans has a bipolar organization characterized by an unusually large structure, spanning over 100 kb. MAT genes have been characterized by functional genetics as being involved in sexual reproduction and virulence. However, classical gene replacement failed to achieve mutants for five MAT genes (RPL22, RPO41, MYO2, PRT1, and RPL39), indicating that they are likely essential. In the present study, targeted gene replacement was performed in a diploid strain for both the α and a alleles of the ribosomal genes RPL22 and RPL39 Mendelian analysis of the progeny confirmed that both RPL22 and RPL39 are essential for viability. Ectopic integration of the RPL22 allele of opposite MAT identity in the heterozygous RPL22a/rpl22αΔ or RPL22α/rpl22aΔ mutant strains failed to complement their essential phenotype. Evidence suggests that this is due to differential expression of the RPL22 genes, and an RNAi-dependent mechanism that contributes to control RPL22a expression. Furthermore, via CRISPR/Cas9 technology, the RPL22 alleles were exchanged in haploid MATα and MATa strains of C. neoformans These RPL22 exchange strains displayed morphological and genetic defects during bilateral mating. These results contribute to elucidating functions of C. neoformans essential mating type genes that may constitute a type of imprinting system to promote inheritance of nuclei of both mating types.

Identification of miRNAs involved in fruit ripening by deep sequencing of Olea europaea L. transcriptome

BACKGROUND

The ripening process of olive fruits is associated with chemical and/or enzymatic specific transformations, making them particularly attractive to animals and humans. In olive drupes, including 'Cassanese' ones, ripening is usually accompanied by progressive chromatic change, resulting in a final red-brown colourization of both epidermis and mesocarp. This event has an exception in the 'Leucocarpa', in which we observed the destabilization in the equilibrium between the chlorophyll metabolism and that of the other pigments, particularly the anthocyanins, whose switch-off during maturation promotes the white colouration of the fruits. Recently, transcription profiling of 'Leucocarpa' and 'Cassanese' olives along ripening, performed through an Illumina RNA-seq approach, has provided useful insights on genes functions involved in fruit maturation such as those related to the biosynthesis of flavonoids and anthocyanins.

METHODOLOGY

To assess expression alterations of genes involved in flavonoids and anthocyanins biosynthetic pathways during ripening, possibly caused by small nuclear RNA (snRNA) in olive drupes, snRNA libraries from 'Leucocarpa' and 'Cassanese' were constructed with RNAs extracted at 100 and 130 Days After Flowering (DAF) and sequenced by an Illumina approach. 130 conserved microRNAs (miRNA) in the Viridiplantae belonging to 14 miRNA families were identified. Regarding the 130 conserved miRNAs, approximately the 48% were identified in all libraries, 5 and 18 miRNAs were shared between the "Cassanese" (C100, C130) and "Leucocarpa" (L100, L130) libraries, respectively.

CONCLUSION

For the remaining reads not-matching with known miRNAs in the Viridiplantae, we combined secondary structure and minimum free energy to discover novel olive miRNAs. Based on these analyses, 492 sequences were considered as putative novel miRNAs. The putative target genes of identified miRNA were computationally predicted by alignment with the olive drupe transcripts obtained from the same samples. A total of 218 transcripts were predicted as targets of 130 known and 492 putative novel miRNAs. Interestingly, some identified target genes are involved in negative regulation of anthocyanin metabolic process. Quantification of the expression pattern of three miRNA and their target transcripts by qRT-PCR assay confirmed the results of Illumina sequencing.

The FvPHR1 transcription factor control phosphate homeostasis by transcriptionally regulating miR399a in woodland strawberry

Plants have evolved phosphate (Pi) starvation response to adapt the low-Pi environment. The regulation of adaptive responses to phosphorus deficiency by the PHR1-miR399-PHO2 module has been well studied in Arabidopsis thaliana but not in strawberry. Transcription factor PHR1 as the central regulator in the Pi starvation signaling has been revealed in a few plant species. However, the function of PHR1 homologues in strawberry is still unknown. In this study, a total of 13 MYB-CC genes were identified in the woodland strawberry (Fragaria vesca) genome and the FvPHR1 gene was characterized. FvPHR1 contains MYB domain and coiled-coil (CC) domain and is localized in the nucleus. FvPHR1 also exhibits trans-activation ability. Furthermore, the P content in leaves of FvPHR1-overexpressing woodland strawberries was significantly increased by 1.38-fold to 1.78-fold compared with that in the wild type. FvPHR1 was also demonstrated to directly bind to the FvMIR399a promoter and positively regulate the expression of FvmiR399a in woodland strawberry. These results showed that PHR1-miR399 module is involved in the regulation of phosphate-signaling pathway and phosphate homeostasis in woodland strawberry.

Small RNA-Omics for Plant Virus Identification, Virome Reconstruction, and Antiviral Defense Characterization

RNA interference (RNAi)-based antiviral defense generates small interfering RNAs that represent the entire genome sequences of both RNA and DNA viruses as well as viroids and viral satellites. Therefore, deep sequencing and bioinformatics analysis of small RNA population (small RNA-ome) allows not only for universal virus detection and genome reconstruction but also for complete virome reconstruction in mixed infections. Viral infections (like other stress factors) can also perturb the RNAi and gene silencing pathways regulating endogenous gene expression and repressing transposons and host genome-integrated endogenous viral elements which can potentially be released from the genome and contribute to disease. This review describes the application of small RNA-omics for virus detection, virome reconstruction and antiviral defense characterization in cultivated and non-cultivated plants. Reviewing available evidence from a large and ever growing number of studies of naturally or experimentally infected hosts revealed that all families of land plant viruses, their satellites and viroids spawn characteristic small RNAs which can be assembled into contigs of sufficient length for virus, satellite or viroid identification and for exhaustive reconstruction of complex viromes. Moreover, the small RNA size, polarity and hotspot profiles reflect virome interactions with the plant RNAi machinery and allow to distinguish between silent endogenous viral elements and their replicating episomal counterparts. Models for the biogenesis and functions of small interfering RNAs derived from all types of RNA and DNA viruses, satellites and viroids as well as endogenous viral elements are presented and discussed.

SOAPnuke: a MapReduce acceleration-supported software for integrated quality control and preprocessing of high-throughput sequencing data

Quality control (QC) and preprocessing are essential steps for sequencing data analysis to ensure the accuracy of results. However, existing tools cannot provide a satisfying solution with integrated comprehensive functions, proper architectures, and highly scalable acceleration. In this article, we demonstrate SOAPnuke as a tool with abundant functions for a “QC-Preprocess-QC” workflow and MapReduce acceleration framework. Four modules with different preprocessing functions are designed for processing datasets from genomic, small RNA, Digital Gene Expression, and metagenomic experiments, respectively. As a workflow-like tool, SOAPnuke centralizes processing functions into 1 executable and predefines their order to avoid the necessity of reformatting different files when switching tools. Furthermore, the MapReduce framework enables large scalability to distribute all the processing works to an entire compute cluster.

We conducted a benchmarking where SOAPnuke and other tools are used to preprocess a ∼30× NA12878 dataset published by GIAB. The standalone operation of SOAPnuke struck a balance between resource occupancy and performance. When accelerated on 16 working nodes with MapReduce, SOAPnuke achieved ∼5.7 times the fastest speed of other tools.

Tie-Break: Host and Retrotransposons Play tRNA

tRNA fragments (tRFs) are a class of small, regulatory RNAs with diverse functions. 3'-Derived tRFs perfectly match long terminal repeat (LTR)-retroelements which use the 3'-end of tRNAs to prime reverse transcription. Recent work has shown that tRFs target LTR-retroviruses and -transposons for the RNA interference (RNAi) pathway and also inhibit mobility by blocking reverse transcription. The highly conserved tRNA primer binding site (PBS) in LTR-retroelements is a unique target for 3'-tRFs to recognize and block abundant but diverse LTR-retrotransposons that become transcriptionally active during epigenetic reprogramming in development and disease. 3'-tRFs are processed from full-length tRNAs under so far unknown conditions and potentially protect many cell types. tRFs appear to be an ancient link between RNAi, transposons, and genome stability.

Comprehensive analysis of small RNAs expressed in developing male strobili of Cryptomeria japonica

Deep sequencing of small RNAs (sRNAs) in developing male strobili of second-generation offspring originating from a nuclear genic male sterile tree of Cryptomeria japonica were performed to characterize sRNA populations in the male strobili at early pollen developmental stages. Comparing to sequences of microRNA (miRNA) families of plant species and sRNAs expressed in the reproductive organs of representative vascular plants, 37 conserved miRNA families were detected, of which eight were ubiquitously expressed in the reproductive organs of land plant species. In contrast, miR1083 was common in male reproductive organs of gymnosperm species but absent in angiosperm species. In addition to conserved miRNAs, 199 novel miRNAs candidates were predicted. The expression patterns of the obtained sRNAs were further investigated to detect the differentially expressed (DE) sRNAs between genic male sterile and fertile individuals. A total of 969 DE sRNAs were obtained and only three known miRNA families were included among them. These results suggest that both conserved and species-specific sRNAs contribute to the development of male strobili in C. japonica.

Transposon-derived small RNAs triggered by miR845 mediate genome dosage response in Arabidopsis

Chromosome dosage has substantial effects on reproductive isolation and speciation in both plants and animals, but the underlying mechanisms are largely obscure 1 . Transposable elements in animals can regulate hybridity through maternal small RNA 2 , whereas small RNAs in plants have been postulated to regulate dosage response via neighboring imprinted genes3,4. Here we show that a highly conserved microRNA in plants, miR845, targets the tRNAMet primer-binding site (PBS) of long terminal repeat (LTR) retrotransposons in Arabidopsis pollen, and triggers the accumulation of 21-22-nucleotide (nt) small RNAs in a dose-dependent fashion via RNA polymerase IV. We show that these epigenetically activated small interfering RNAs (easiRNAs) mediate hybridization barriers between diploid seed parents and tetraploid pollen parents (the 'triploid block'), and that natural variation for miR845 may account for 'endosperm balance' allowing the formation of triploid seeds. Targeting of the PBS with small RNA is a common mechanism for transposon control in mammals and plants, and provides a uniquely sensitive means to monitor chromosome dosage and imprinting in the developing seed.

A multi-omics study of the grapevine-downy mildew (Plasmopara viticola) pathosystem unveils a complex protein coding- and noncoding-based arms race during infection

Fungicides are applied intensively to prevent downy mildew infections of grapevines (Vitis vinifera) with high impact on the environment. In order to develop alternative strategies we sequenced the genome of the oomycete pathogen Plasmopara viticola causing this disease. We show that it derives from a Phytophthora-like ancestor that switched to obligate biotrophy by losing genes involved in nitrogen metabolism and γ-Aminobutyric acid catabolism. By combining multiple omics approaches we characterized the pathosystem and identified a RxLR effector that trigger an immune response in the wild species V. riparia. This effector is an ideal marker to screen novel grape resistant varieties. Our study reveals an unprecedented bidirectional noncoding RNA-based mechanism that, in one direction might be fundamental for P. viticola to proficiently infect its host, and in the other might reduce the effects of the infection on the plant.

tRNA-derived small RNAs target transposable element transcripts

tRNA-derived RNA fragments (tRFs) are 18–26 nucleotide small RNAs that are not random degradation products, but are rather specifically cleaved from mature tRNA transcripts. Abundant in stressed or viral-infected cells, the function and potential targets of tRFs are not known. We identified that in the unstressed wild-type male gamete containing pollen of flowering plants, and analogous reproductive structure in non-flowering plant species, tRFs accumulate to high levels. In the reference plant Arabidopsis thaliana, tRFs are processed by Dicer-like 1 and incorporated into Argonaute1 (AGO1), akin to a microRNA. We utilized the fact that many plant small RNAs direct cleavage of their target transcripts to demonstrate that the tRF–AGO1 complex acts to specifically target and cleave endogenous transposable element (TE) mRNAs produced from transcriptionally active TEs. The data presented here demonstrate that tRFs are bona-fide regulatory microRNA-like small RNAs involved in the regulation of genome stability through the targeting of TE transcripts.

Identification of Conserved and Novel MicroRNAs in Blueberry

MicroRNAs (miRNAs) are a class of small endogenous RNAs that play important regulatory roles in cells by negatively affecting gene expression at both transcriptional and post-transcriptional levels. There have been extensive studies aiming to identify miRNAs and to elucidate their functions in various plant species. In the present study, we employed the high-throughput sequencing technology to profile miRNAs in blueberry fruits. A total of 9,992,446 small RNA tags with sizes ranged from 18 to 30 nt were obtained, indicating that blueberry fruits have a large and diverse small RNA population. Bioinformatic analysis identified 412 conserved miRNAs belonging to 29 families, and 35 predicted novel miRNAs that are likely to be unique to blueberries. Among them, expression profiles of five conserved miRNAs were validated by stem loop qRT-PCR. Furthermore, the potential target genes of conserved and novel miRNAs were predicted and subjected to Gene Ontology (GO) annotation. Enrichment analysis of the GO-represented biological processes and molecular functions revealed that these target genes were potentially involved in a wide range of metabolic pathways and developmental processes. Particularly, anthocyanin biosynthesis has been predicted to be directly or indirectly regulated by diverse miRNA families. This study is the first report on genome-wide miRNA profile analysis in blueberry and it provides a useful resource for further elucidation of the functional roles of miRNAs during fruit development and ripening.