Identification of novel miRNAs and their target genes in <Emphasis Type="Italic">Eucalyptus grandis</Emphasis>

Despite Eucalyptus grandis being the most widely planted hardwood tree globally, along with the availability of a sequenced genome and easily accessible functional genetic tools, the quantities and roles of miRNA in its developmental processes remains largely unknown. In this study, we constructed small RNA libraries by high-throughput sequencing from Eucalyptus grandis samples, and 386 novel miRNAs were identified by miRDeep2. We found 179 novel miRNAs, 41 miRNA families, and 456 target genes in leaf samples, and 257 novel miRNAs, 61 miRNA families, and 483 target genes in stem samples. The function of the MIR396 family of miRNAs in Eucalyptus grandis was found to be mainly associated with the process of cell growth. By annotation analysis of miRNA targets, we found that some target genes, such as GRF, expansin-A15, and RPS2, had a close correlation in stem. Finally, the three randomly selected members of the MIR396 family were confirmed to express in Eucalyptus grandis by qRT-PCR, indicating that our reported miRNAs were existed. The identification of miRNAs and their target genes will lead to a greater understanding of the role of miRNAs in the physiology, growth, and development of Eucalyptus grandis trees.     

Comparative analysis of microRNAs and putative target genes in hybrid clone <Emphasis Type="Italic">Paulownia</Emphasis> ‘yuza 1’ under drought stress

Drought stress in plants often leads to reduced productivity and limited geographic distribution, which can affect human life and ecosystems. The responses of diploid and tetraploid Paulownia tomentosa × Paulownia fortunei to drought have been reported, but the effects of drought stress on the levels of microRNA (miRNA) expression have not been published so far. Here, we constructed four small RNA (sRNA) libraries and four corresponding degradome libraries of well-watered and severe drought-treated diploid and tetraploid plants to identify the miRNAs and their putative target genes in Paulownia ‘yuza 1’, a P. tomentosa × P. fortunei hybrid clone, by sRNA and degradome sequencing. The putative target genes of miRNAs were annotated with gene ontology terms and Kyoto Encyclopedia of Genes and Genomes pathways. Three conserved and 21 novel miRNAs responsive to drought stress were found, in which 15 were identified as the main drought responsive miRNAs that conferred higher resistance in tetraploid than in diploid of Paulownia ‘yuza 1’. Our results will lay the foundation for investigating the roles of miRNAs in Paulownia and other trees in response to drought.     

The use of miRNAs as reference genes for miRNA expression normalization during <Emphasis Type="Italic">Lilium</Emphasis> somatic embryogenesis by real-time reverse transcription PCR analysis

Expression profiling of miRNAs has the ability to reveal the essence of somatic embryogenesis (SE). qRT-PCR is one of the most commonly used techniques for dynamic miRNA detection but requires optimal reference genes for data reliability. This is the first report on reference gene validation for miRNA expression normalization in Lilium (Lilium pumilum DC. Fisch. and Lilium davidii var. unicolor). In this study, seventeen miRNAs together with two snRNAs (U4, U6), one rRNA (5S rRNA) and three protein-coding genes (FP, ACT, GAPDH) were selected as reference candidates, and their expression stability was validated by qRT-PCR among eleven developing SE cultures in two lilies. Four normalization algorithms, including geNorm, BestKeeper, NormFinder and RefFinder, were also used to evaluate the stability of the reference candidates. For Lilium pumilum DC. Fisch., lpu-miR159a was the optimal reference gene during SE, followed by lpu-miR408b, while U6 was the least stable reference candidate. For Lilium davidii var. unicolor, FP presented greater stability than did half of the miRNA candidates, but the best reference gene was lda-miR162, followed by lda-miR159a. Further analysis of the expression level of miR156 and miR529 was used to evaluate the validity of the reference genes in both lilies. In general, miRNAs are superior to common protein-coding genes and snRNAs / rRNAs as reference genes for miRNA expression normalization during Lilium SE, and the most suitable reference miRNA is different between two species in the same Lilium genus. This is a pioneer study using suitable miRNAs as reference genes in Lilium and constitutes a small but essential step for the further exploration of miRNA function in Lilium, thus offering valuable references for other plants.     

Phylogeny and Molecular Evolution of miR820 and miR396 microRNA Families in <Emphasis Type="Italic">Oryza</Emphasis> AA Genomes

The phylogeny and evolution of the microRNA families, miR820 and miR396, was analysed across the AA genomes of the Oryza species, the close relatives of domesticated rice. A highly dynamic evolution of the miR820 family was revealed. The number of copies of MIR820 genes, their chromosomal location and the mature microRNA sequence varied greatly with a total of 16 novel miR820 variants being identified. The phylogeny of pre-MIR820 sequences revealed that MIR820 genes of recently evolved Oryza AA genomes may have derived from sequence divergence of one or a few ancestral genes found in wild Australian perennial rice populations, Taxon B (jpn2)-MIR820 genes. Genomic scale duplication played an important role in the evolution of some of the miR396 family genes in AA genome Oryza species. miR396 family contained a MIR396 gene cluster (MIR396a and MIR396c) which was conserved across the cereal genomes. Nucleotide diversity analysis at these two MIR396 loci revealed that domesticated rice has retained less than 10% of the total diversity present in wild species. In contrast, the nucleotide sequence of four MIR396 loci remained almost conserved across domesticated and wild rices, indicating that they were under extreme functional constraint and may be involved in regulating some fundamental processes which are important both for wild and domesticated rices. Expression analysis demonstrated that miR820 variants were expressed in O. glaberrima O. barthi and O. longistaminata genome. These findings pose new challenges to explain the possible role of miR820 variants identified.     

Deep sequencing discovery and profiling of conserved and novel miRNAs in the ovule of <Emphasis Type="Italic">Ginkgo biloba</Emphasis>

Key message

High-throughput sequencing and subsequent analysis identified multiple miRNAs closely related to ovule, indicating that miRNAs are important in Ginkgo biloba ovule.

Abstract

MicroRNAs (miRNAs) are small, noncoding, regulatory RNAs that play crucial regulatory roles in the process of plant growth and development. However, limited information regarding their functions in gymnosperm reproduction is available. Here, we used high-throughput sequencing combined with computational analysis to identify and characterize miRNAs from ovules of G. biloba, and identified 34 conserved miRNA families and 99 novel miRNAs. The precursor sequences of several of the conserved and novel miRNAs were further validated by RT-PCR and sequencing. Furthermore, we found that some target genes, e.g. MYB, homeodomain-leucine zipper (HD-ZIPIII) and auxin response factor (ARF), may be involved in ovule development, and that the significantly enriched pathways of some miRNA targets were related to plant–pathogen interactions and the biosynthesis of secondary metabolites. Twenty-six conserved miRNA families were found to be expressed in both leaves and ovules, while miRNA156, miRNA164, miRNA167, miRNA169, miRNA172 and miRNA390 were up-regulated in ovules. Thus, multiple miRNAs closely related to G. biloba ovule development were identified, resulting in a greater understanding of the important regulatory functions of miRNAs in plant ovules.

     

miRDeepFinder: a miRNA analysis tool for deep sequencing of plant small RNAs

miRDeepFinder is a software package developed to identify and functionally analyze plant microRNAs (miRNAs) and their targets from small RNA datasets obtained from deep sequencing. The functions available in miRDeepFinder include pre-processing of raw data, identifying conserved miRNAs, mining and classifying novel miRNAs, miRNA expression profiling, predicting miRNA targets, and gene pathway and gene network analysis involving miRNAs. The fundamental design of miRDeepFinder is based on miRNA biogenesis, miRNA-mediated gene regulation and target recognition, such as perfect or near perfect hairpin structures, different read abundances of miRNA and miRNA*, and targeting patterns of plant miRNAs. To test the accuracy and robustness of miRDeepFinder, we analyzed a small RNA deep sequencing dataset of Arabidopsis thaliana published in the GEO database of NCBI. Our test retrieved 128 of 131 (97.7%) known miRNAs that have a more than 3 read count in Arabidopsis. Because many known miRNAs are not associated with miRNA*s in small RNA datasets, miRDeepFinder was also designed to recover miRNA candidates without the presence of miRNA*. To mine as many miRNAs as possible, miRDeepFinder allows users to compare mature miRNAs and their miRNA*s with other small RNA datasets from the same species. Cleaveland software package was also incorporated into miRDeepFinder for miRNA target identification using degradome sequencing analysis. Using this new computational tool, we identified 13 novel miRNA candidates with miRNA*s from Arabidopsis and validated 12 of them experimentally. Interestingly, of the 12 verified novel miRNAs, a miRNA named AC1 spans the exons of two genes (UTG71C4 and UGT71C3). Both the mature AC1 miRNA and its miRNA* were also found in four other small RNA datasets. We also developed a tool, ??miRNA primer designer?? to design primers for any type of miRNAs. miRDeepFinder provides a powerful tool for analyzing small RNA datasets from all species, with or without the availability of genome information. miRDeepFinder and miRNA primer designer are freely available at http://www.leonxie.com/DeepFinder.php and at http://www.leonxie.com/miRNAprimerDesigner.php, respectively. A program (called RefFinder: http://www.leonxie.com/referencegene.php) was also developed for assessing the reliable reference genes for gene expression analysis, including miRNAs.