PmiREN2.0: from data annotation to functional exploration of plant microRNAs.

Nearly 200 plant genomes have been sequenced over the last two years, and new functions of plant microRNAs (miRNAs) have been revealed. Therefore, timely update of the plant miRNA databases by incorporating miRNAs from the newly sequenced species and functional information is required to provide useful resources for advancing plant miRNA research. Here we report the update of PmiREN2.0 (https://pmiren.com/) with an addition of 19 363 miRNA entries from 91 plants, doubling the amount of data in the original version. Meanwhile, abundant regulatory information centred on miRNAs was added, including predicted upstream transcription factors through binding motifs scanning and elaborate annotation of miRNA targets. As an example, a genome-wide regulatory network centred on miRNAs was constructed for Arabidopsis. Furthermore, phylogenetic trees of conserved miRNA families were built to expand the understanding of miRNA evolution across the plant lineages. These data are helpful to deduce the regulatory relationships concerning miRNA functions in diverse plants. Beside the new data, a suite of design tools was incorporated to facilitate experimental practice. Finally, a forum named 'PmiREN Community' was added for discussion and resource and new discovery sharing. With these upgrades, PmiREN2.0 should serve the community better and accelerate miRNA research in plants.

Miniature Inverted-repeat Transposable Elements Drive Rapid MicroRNA Diversification in Angiosperms.

MicroRNAs (miRNAs) are fast evolving endogenous small RNAs that regulate organism function and behavior in both animals and plants. Although models for de novo miRNA biogenesis have been proposed, the genomic mechanisms driving swift diversification of the miRNA repertoires in plants remain elusive. Here, by comprehensively analyzing 21 phylogenetically representative plant species, ranging from green algae to angiosperms, we systematically identified de novo miRNA events associated with 8,649 miRNA loci. We found that 399 (4.6%), 466 (5.4%), and 1,402 (16.2%) miRNAs were derived from inverted gene duplication events, long terminal repeats of retrotransposons, and miniature inverted-repeat transposable elements (MITEs), respectively. Among the miRNAs of these origins, MITEs, especially those belonging to the Mutator, Tc1/Mariner, and PIF/Harbinger superfamilies, were the predominant genomic source for de novo miRNAs in the 15 examined angiosperms but not in the six non-angiosperms. Our data further illustrated a transposition-transcription process by which MITEs are converted into new miRNAs (termed MITE-miRNAs) whereby properly sized MITEs are transcribed and therefore become potential substrates for the miRNA processing machinery by transposing into introns of active genes. By analyzing the 58,038 putative target genes for the 8,095 miRNAs, we found that the target genes of MITE-miRNAs were preferentially associated with response to environmental stimuli such as temperature, suggesting that MITE-miRNAs are pertinent to plant adaptation. Collectively, these findings demonstrate that molecular conversion of MITEs is a genomic mechanism leading to rapid and continuous changes to the miRNA repertoires in angiosperm.

Quo vadis microRNAs?

Since 2002, published miRNAs have been collected and named by the online repository miRBase. However, with 11 000 annual publications this has become challenging. Recently, four specialized miRNA databases were published, addressing particular needs for diverse scientific communities. This development provides major opportunities for the future of miRNA annotation and nomenclature.

MicroRNA annotation in plants: current status and challenges.

Last two decades, the studies on microRNAs (miRNAs) and the numbers of annotated miRNAs in plants and animals have surged. Herein, we reviewed the current progress and challenges of miRNA annotation in plants. Via the comparison of plant and animal miRNAs, we pinpointed out the difficulties on plant miRNA annotation and proposed potential solutions. In terms of recalling the history of methods and criteria in plant miRNA annotation, we detailed how the major progresses made and evolved. By collecting and categorizing bioinformatics tools for plant miRNA annotation, we surveyed their advantages and disadvantages, especially for ones with the principle of mimicking the miRNA biogenesis pathway by parsing deeply sequenced small RNA (sRNA) libraries. In addition, we summarized all available databases hosting plant miRNAs, and posted the potential optimization solutions such as how to increase the signal-to-noise ratio (SNR) in these databases. Finally, we discussed the challenges and perspectives of plant miRNA annotations, and indicated the possibilities offered by an all-in-one tool and platform according to the integration of artificial intelligence.

Insights into the convergent evolution of fructan biosynthesis in angiosperms from the highly characteristic chicory genome.

Fructans in angiosperms play essential roles in physiological functions and environmental adaptations. As a major source of industrial fructans (especially inulin-type), chicory (Cichorium intybus L.) is a model species for studying fructan biosynthesis. However, the genes underlying this process and their evolutionary history in angiosperms remain elusive. We combined multiple sequencing technologies to assemble and annotate the chicory genome and scan its (epi)genomic features, such as genomic components, DNA methylation, and three-dimensional (3D) structure. We also performed a comparative genomics analysis to uncover the associations between key traits and gene families. We achieved a nearly complete chicory genome assembly and found that continuous bursts of a few highly active retrotransposon families largely shaped the (epi)genomic characteristics. The highly methylated genome with its unique 3D structure potentially influences critical biological processes. Our comprehensive comparative genomics analysis deciphered the genetic basis for the rich sesquiterpene content in chicory and indicated that the fructan-accumulating trait resulted from convergent evolution in angiosperms due to shifts in critical sites of fructan-active enzymes. The highly characterized chicory genome provides insight into Asteraceae evolution and fructan biosynthesis in angiosperms.

PmiREN: a comprehensive encyclopedia of plant miRNAs.

MicroRNAs (miRNAs) are small non-coding RNA molecules that function as diverse endogenous gene regulators at the post-transcriptional level. In the past two decades, as research effort on miRNA identification, function and evolution has soared, so has the demand for miRNA databases. However, the current plant miRNA databases suffer from several typical drawbacks, including a lack of entries for many important species, uneven annotation standards across different species, abundant questionable entries, and limited annotation. To address these issues, we developed a knowledge-based database called Plant miRNA Encyclopedia (PmiREN, http://www.pmiren.com/), which was based on uniform processing of sequenced small RNA libraries using miRDeep-P2, followed by manual curation using newly updated plant miRNA identification criteria, and comprehensive annotation. PmiREN currently contains 16,422 high confidence novel miRNA loci in 88 plant species and 3,966 retrieved from miRBase. For every miRNA entry, information on precursor sequence, precursor secondary structure, expression pattern, clusters and synteny in the genome, potential targets supported by Parallel Analysis of RNA Ends (PARE) sequencing, and references is attached whenever possible. PmiREN is hierarchically accessible and has eight built-in search engines. We believe PmiREN is useful for plant miRNA cataloguing and data mining, therefore a resource for data-driven miRNA research in plants.

Identification of Species-Specific MicroRNAs Provides Insights into Dynamic Evolution of MicroRNAs in Plants.

MicroRNAs (miRNAs) are an important class of regulatory small RNAs that program gene expression, mainly at the post-transcriptional level. Although sporadic examples of species-specific miRNAs (termed SS-miRNAs) have been reported, a genome-scale study across a variety of distant species has not been assessed. Here, by comprehensively analyzing miRNAs in 81 plant species phylogenetically ranging from chlorophytes to angiosperms, we identified 8048 species-specific miRNAs from 5499 families, representing over 61.2% of the miRNA families in the examined species. An analysis of the conservation from different taxonomic levels supported the high turnover rate of SS-miRNAs, even over short evolutionary distances. A comparison of the intrinsic features between SS-miRNAs and NSS-miRNAs (non-species-specific miRNAs) indicated that the AU content of mature miRNAs was the most striking difference. Our data further illustrated a significant bias of the genomic coordinates towards SS-miRNAs lying close to or within genes. By analyzing the 125,267 putative target genes for the 7966 miRNAs, we found the preferentially regulated functions of SS-miRNAs related to diverse metabolic processes. Collectively, these findings underscore the dynamic evolution of miRNAs in the species-specific lineages.

miRDeep-P2: accurate and fast analysis of the microRNA transcriptome in plants.

MOTIVATION: Two major challenges arise when employing next-generation sequencing methods to comprehensively identify microRNAs (miRNAs) in plants: (i) how to minimize the false-positive inheritable to computational predictions and (ii) how to minimize the computational time required for analyzing the miRNA transcriptome in plants with complex and large genomes. RESULTS: We updated miRDeep-P to miRDeep-P2 (miRDP2) by employing a new filtering strategy and overhauling the algorithm. miRDP2 has been tested against miRNA transcriptomes in plants with increasing genome sizes that included Arabidopsis, rice, tomato, maize and wheat. Compared with miRDeep-P and several other computational tools, miRDP2 processes next-generation sequencing data with superior speed. By incorporating newly updated plant miRNA annotation criteria and developing a new scoring system, the accuracy of miRDP2 outperformed other programs. Taken together, our results demonstrate miRDP2 as a fast and accurate tool for analyzing the miRNA transcriptome in plants. AVAILABILITY AND IMPLEMENTATION: The miRDP2 are freely available from https://sourceforge.net/projects/mirdp2/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Selection of suitable reference genes for qRT-PCR expression analysis of Codonopsis pilosula under different experimental conditions.

Codonopsis pilosula is a well-known medicinal plant. Although its transcriptome sequence has been published, suitable reference genes have not been systematically identified for conducting expression analyses via quantitative real-time polymerase chain reaction (qRT-PCR). To screen appropriate genes for use with this species, we applied four different methods-GeNorm, NormFinder, BestKeeper, and RefFinder-to evaluate the stability of 13 candidates: CpiEF1Bb, CpiCACS, CpiF-Box, Cpiß-Tubulin, CpiGAPDH, CpiActin2, CpiAPT1, CpiActin7, CpiActin8, CpiRPL6, CpiHAF1, CpiTubulin6, and CpiUBQ12. Expression was examined by qRT-PCR for various tissue types, chemical treatments, and developmental stages. For all tested samples, CpiGAPDH proved to be the most stable. Comprehensive analysis indicated that the most stable internal reference genes were CpiF-Box and CpiCACS in different tissues and at different developmental stages, respectively. Under NaCl stress, CpiAPT1 was the best internal reference gene. For methyl jasmonate and abscisic acid treatments, CpiGAPDH and CpiF-Box, respectively, presented the highest degree of expression stability. Based on these findings, we chose CpiSPL9 as the target gene for validating the suitability of these selected reference genes. All of these results provide a foundation for accurate quantification of expression levels by genes of interest in C. pilosula.

Expression of a novel bi-directional Brassica napus promoter in soybean.

The expression profile of a natural bi-directional promoter, derived from the Brassica napus EPSPS-A gene, was studied in transgenic soybean (Glycine max C.V. Maverick) lines. Two constructs, pDAB100331 and pDAB100333, were assembled to test the bi-directionality of the promoter. Two reporter genes, gfp and gusA, were employed and they were interchangeably placed in both constructs, one on each end of the promoter such that both proteins expressed divergently in each construct. In the T0 generation, GUS expression was more uniform throughout the leaf of pDAB100333 transgenic plants, where the gusA gene was expressed from the downstream or EPSPS-A end of the bi-directional promoter. Comparatively, GUS expression was more localized in the midrib and veins of the leaf of pDAB100331 transgenic plants, where the gusA gene was expressed from the upstream end of the bi-directional promoter. These observations indicated a unique expression pattern from each end of the promoter and consistently higher expression in genes expressed from the downstream end (e.g., EPSPS-A end) of the promoter in the tissues examined. The GFP expression pattern followed that of GUS when placed in the same position relative to the promoter. In the T1 generation, transcript analysis also showed higher expression of both gusA and gfp when those genes were located at the downstream end of the promoter. Accordingly, the pDAB100331 events exhibited a higher gfp/gusA transcript ratio, while pDAB100333 events produced a higher gusA/gfp transcript ratio consistent with the observations in T0 plants. These results demonstrated that the EPSPS-A gene bidirectional promoter can be effectively utilized to drive expression of two transgenes for the desired traits.

Identification and characterization of a subset of microRNAs in wheat (Triticum aestivum L.).

MicroRNAs (miRNAs) represent a class of endogenous regulator for post-transcriptionally modulating gene expression. Elucidating complete miRNA repertoires for individual species is a long-desired goal in miRNA research. So far only 42 have been annotated for common wheat (Triticum aestivum) due to its large genome. Here, we employed miRDeep-P, a program developed previously for retrieving miRNAs from deep sequencing data in plants, to parse 14 sequenced small RNA libraries of wheat using expression sequence tags (ESTs) as the reference in lieu of a complete genome sequence. This effort identified 145 miRNAs along with the corresponding stem-looped precursors with many differentially expressed in development and associated with powdery mildew infection. Examination of the phylogenetic distribution of these miRNAs and their target genes revealed that many are conserved in monocots. The set of miRNAs identified in this study is thus useful to further miRNA research in wheat and other important cereal crop species.

Alternative mRNA processing increases the complexity of microRNA-based gene regulation in Arabidopsis.

MicroRNAs (miRNAs) represent an important class of sequence-specific, trans-acting endogenous small RNA molecules that modulate gene expression at the post-transcriptional level. They function by binding to partial complementary cis-regulatory sites (miRNA binding sites) in their target mRNAs. Based on two recent observations from plant genome studies, namely that alternative splicing is a common phenomenon and that miRNA regulates a significant proportion of the transcriptome, we hypothesize that there may be a mechanism for gene regulation that involves both processes. In the present study, we performed a systemic search in the model plant Arabidopsis thaliana using annotated gene models as well as publically available high-throughput RNA sequencing data with a total of 570 million reads. Of the 354 high-confidence miRNA binding sites identified in Arabidopsis, at least 44 (12.4%) were affected by alternative splicing such that mRNA isoforms of the same miRNA target gene differ in the sequences encoding the miRNA binding sites. By simulation, we found that the frequency of alternative splicing at miRNA binding sites is significantly higher than at other regions. Comparative and functional analyses further indicated that the alternative splicing events are important for target gene expression and miRNA action. Together our results show that alternative splicing of miRNA binding sites is a plausible mechanism for attenuating miRNA-mediated gene regulation.

Plant MicroRNA Identification and Annotation Using Deep Sequencing Data.

MicroRNAs (miRNAs) are endogenous non-coding small RNAs, which regulate gene expression at the post-transcriptional level. A large number of studies have revealed that they play key roles in diverse life activities, such as growth and development. In the last decade, deep sequencing technology has generated substantial small RNA sequencing (sRNA-Seq) data. Meanwhile, numerous tools have been developed to identify miRNAs from these sRNA-Seq data, resulting in a surge of miRNA annotations. Among these tools, the series of miRDeep-P and miRDeep-P2 have been widely used in plant miRNA annotation. Here, we employed miRDeep-P2 to demonstrate the plant miRNA annotation processes step by step using the deep sequencing data.

Advances in alternative splicing identification: deep learning and pantranscriptome.

In plants, alternative splicing is a crucial mechanism for regulating gene expression at the post-transcriptional level, which leads to diverse proteins by generating multiple mature mRNA isoforms and diversify the gene regulation. Due to the complexity and variability of this process, accurate identification of splicing events is a vital step in studying alternative splicing. This article presents the application of alternative splicing algorithms with or without reference genomes in plants, as well as the integration of advanced deep learning techniques for improved detection accuracy. In addition, we also discuss alternative splicing studies in the pan-genomic background and the usefulness of integrated strategies for fully profiling alternative splicing.

Comparative genomics reveals a unique nitrogen-carbon balance system in Asteraceae.

The Asteraceae (daisy family) is one of the largest families of plants. The genetic basis for its high biodiversity and excellent adaptability has not been elucidated. Here, we compare the genomes of 29 terrestrial plant species, including two de novo chromosome-scale genome assemblies for stem lettuce, a member of Asteraceae, and Scaevola taccada, a member of Goodeniaceae that is one of the closest outgroups of Asteraceae. We show that Asteraceae originated ~80 million years ago and experienced repeated paleopolyploidization. PII, the universal regulator of nitrogen-carbon (N-C) assimilation present in almost all domains of life, has conspicuously lost across Asteraceae. Meanwhile, Asteraceae has stepwise upgraded the N-C balance system via paleopolyploidization and tandem duplications of key metabolic genes, resulting in enhanced nitrogen uptake and fatty acid biosynthesis. In addition to suggesting a molecular basis for their ecological success, the unique N-C balance system reported for Asteraceae offers a potential crop improvement strategy.

LettuceGDB: The community database for lettuce genetics and omics.

As a globally popular leafy vegetable and a representative plant of the Asteraceae family, lettuce has great economic and academic significance. In the last decade, high-throughput sequencing, phenotyping, and other multi-omics data in lettuce have accumulated on a large scale, thus increasing the demand for an integrative lettuce database. Here, we report the establishment of a comprehensive lettuce database, LettuceGDB (https://www.lettucegdb.com/). As an omics data hub, the current LettuceGDB includes two reference genomes with detailed annotations; re-sequencing data from over 1000 lettuce varieties; a collection of more than 1300 worldwide germplasms and millions of accompanying phenotypic records obtained with manual and cutting-edge phenomics technologies; re-analyses of 256 RNA sequencing datasets; a complete miRNAome; extensive metabolite information for representative varieties and wild relatives; epigenetic data on the genome-wide chromatin accessibility landscape; and various lettuce research papers published in the last decade. Five hierarchically accessible functions (Genome, Genotype, Germplasm, Phenotype, and O-Omics) have been developed with a user-friendly interface to enable convenient data access. Eight built-in tools (Assembly Converter, Search Gene, BLAST, JBrowse, Primer Design, Gene Annotation, Tissue Expression, Literature, and Data) are available for data downloading and browsing, functional gene exploration, and experimental practice. A community forum is also available for information sharing, and a summary of current research progress on different aspects of lettuce is included. We believe that LettuceGDB can be a comprehensive functional database amenable to data mining and database-driven exploration, useful for both scientific research and lettuce breeding.

Genome-wide mapping of the HY5-mediated gene networks in Arabidopsis that involve both transcriptional and post-transcriptional regulation.

LONG HYPOCOTYL 5 (HY5) is a basic leucine zipper transcription factor (TF) that functions downstream of multiple families of photoreceptors. Mutations in the HY5 gene cause a myriad of aberrant phenotypes in Arabidopsis, including elongated hypocotyl, reduced accumulation of pigments, halted chloroplast development in greening hypocotyls, altered root morphology, and defective hormonal and stimulus responses. HY5 thus acts as an integrator that links various gene networks to coordinate plant development. Here we report the experimental mapping of HY5-mediated gene networks in Arabidopsis by integrating genomic loci occupied by HY5 and HY5-dependent gene expression profiles. Our results indicate that HY5 binds to over 9000 genes, detectably affecting the expression of over 1100 genes, either positively or negatively. Further, HY5 indirectly regulate many other genes through sub-networks mediated by other regulators. In particular, HY5 regulates eight miRNA genes that in turn control the transcript abundance of specific target genes. Over-expressing HY5-targeted miR408 resulted in phenotypes that are opposite to the hy5 mutants. Together, our results reveal both transcriptional and post-transcriptional components of the HY5-mediated gene networks.

Asymmetric purine-pyrimidine distribution in cellular small RNA population of papaya.

BACKGROUND: The small RNAs (sRNA) are a regulatory class of RNA mainly represented by the 21 and 24-nucleotide size classes. The cellular sRNAs are processed by RNase III family enzyme dicer (Dicer like in plant) from a self-complementary hairpin loop or other type of RNA duplexes. The papaya genome has been sequenced, but its microRNAs and other regulatory RNAs are yet to be analyzed. RESULTS: We analyzed the genomic features of the papaya sRNA population from three sRNA deep sequencing libraries made from leaves, flowers, and leaves infected with Papaya Ringspot Virus (PRSV). We also used the deep sequencing data to annotate the micro RNA (miRNA) in papaya. We identified 60 miRNAs, 24 of which were conserved in other species, and 36 of which were novel miRNAs specific to papaya. In contrast to the Chargaff's purine-pyrimidine equilibrium, cellular sRNA was significantly biased towards a purine rich population. Of the two purine bases, higher frequency of adenine was present in 23nt or longer sRNAs, while 22nt or shorter sRNAs were over represented by guanine bases. However, this bias was not observed in the annotated miRNAs in plants. The 21nt species were expressed from fewer loci but expressed at higher levels relative to the 24nt species. The highly expressed 21nt species were clustered in a few isolated locations of the genome. The PRSV infected leaves showed higher accumulation of 21 and 22nt sRNA compared to uninfected leaves. We observed higher accumulation of miRNA* of seven annotated miRNAs in virus-infected tissue, indicating the potential function of miRNA* under stressed conditions. CONCLUSIONS: We have identified 60 miRNAs in papaya. Our study revealed the asymmetric purine-pyrimidine distribution in cellular sRNA population. The 21nt species of sRNAs have higher expression levels than 24nt sRNA. The miRNA* of some miRNAs shows higher accumulation in PRSV infected tissues, suggesting that these strands are not totally functionally redundant. The findings open a new avenue for further investigation of the sRNA silencing pathway in plants.

miRDeep-P: a computational tool for analyzing the microRNA transcriptome in plants.

MOTIVATION: Ultra-deep sampling of small RNA libraries by next-generation sequencing has provided rich information on the microRNA (miRNA) transcriptome of various plant species. However, few computational tools have been developed to effectively deconvolute the complex information. RESULTS: We sought to employ the signature distribution of small RNA reads along the miRNA precursor as a model in plants to profile expression of known miRNA genes and to identify novel ones. A freely available package, miRDeep-P, was developed by modifying miRDeep, which is based on a probabilistic model of miRNA biogenesis in animals, with a plant-specific scoring system and filtering criteria. We have tested miRDeep-P on eight small RNA libraries derived from three plants. Our results demonstrate miRDeep-P as an effective and easy-to-use tool for characterizing the miRNA transcriptome in plants. AVAILABILITY: http://faculty.virginia.edu/lilab/miRDP/ CONTACT: ll4jn@virginia.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Global analysis of gene-level microRNA expression in Arabidopsis using deep sequencing data.

MicroRNAs (miRNAs) regulate gene expression at the post-transcriptional level in eukaryotes. Exclusive focus on mature miRNA in most expression profiling efforts has prevented effective measurement of the expression of individual miRNA (MIR) genes. Using three sequenced small RNA libraries, we adapted miRDeep, which employs a probabilistic model of miRNA biogenesis, to analyze the miRNA transcriptome in Arabidopsis. We determined that less than 40% annotated MIR genes are expressed in shoot, root or inflorescence. We found that within paralogous families the expression pattern of individual genes correlates with the phylogenetic distance. Combining novel candidates identified in this study, we deduced the maximal number of expressed MIR genes. We further estimated the sequencing depth necessary to reach a near-saturated detection rate by curve fitting simulation. These results demonstrate that signature distribution of small RNA reads along the miRNA precursor is an effective model to profile MIR gene expression in Arabidopsis.