MicroRNAs in the moss Physcomitrella patens

Biogenesis, conservation, and function of miRNA in liverworts

MicroRNAs (miRNAs) are small non-coding endogenous RNA molecules, 18-24 nucleotides long, that control multiple gene regulatory pathways via post-transcriptional gene silencing in eukaryotes. To develop a comprehensive picture of the evolutionary history of miRNA biogenesis and action in land plants, studies on bryophyte representatives are needed. Here, we review current understanding of liverwort MIR gene structure, miRNA biogenesis, and function, focusing on the simple thalloid Pellia endiviifolia and the complex thalloid Marchantia polymorpha. We review what is known about conserved and non-conserved miRNAs, their targets, and the functional implications of miRNA action in M. polymorpha and P. endiviifolia. We note that most M. polymorpha miRNAs are encoded within protein-coding genes and provide data for 23 MIR gene structures recognized as independent transcriptional units. We identify M. polymorpha genes involved in miRNA biogenesis that are homologous to those identified in higher plants, including those encoding core microprocessor components and other auxiliary and regulatory proteins that influence the stability, folding, and processing of pri-miRNAs. We analyzed miRNA biogenesis proteins and found similar domain architecture in most cases. Our data support the hypothesis that almost all miRNA biogenesis factors in higher plants are also present in liverworts, suggesting that they emerged early during land plant evolution.

The Evolution and Functional Roles of miR408 and Its Targets in Plants

MicroRNA408 (miR408) is an ancient and highly conserved miRNA, which is involved in the regulation of plant growth, development and stress response. However, previous research results on the evolution and functional roles of miR408 and its targets are relatively scattered, and there is a lack of a systematic comparison and comprehensive summary of the detailed evolutionary pathways and regulatory mechanisms of miR408 and its targets in plants. Here, we analyzed the evolutionary pathway of miR408 in plants, and summarized the functions of miR408 and its targets in regulating plant growth and development and plant responses to various abiotic and biotic stresses. The evolutionary analysis shows that miR408 is an ancient and highly conserved microRNA, which is widely distributed in different plants. miR408 regulates the growth and development of different plants by down-regulating its targets, encoding blue copper (Cu) proteins, and by transporting Cu to plastocyanin (PC), which affects photosynthesis and ultimately promotes grain yield. In addition, miR408 improves tolerance to stress by down-regulating target genes and enhancing cellular antioxidants, thereby increasing the antioxidant capacity of plants. This review expands and promotes an in-depth understanding of the evolutionary and regulatory roles of miR408 and its targets in plants.

Viral suppressor of RNA silencing in vascular plants also interferes with the development of the bryophyte Physcomitrella patens

Plant viruses are important pathogens able to overcome plant defense mechanisms using their viral suppressors of RNA silencing (VSR). Small RNA pathways of bryophytes and vascular plants have significant similarities, but little is known about how viruses interact with mosses. This study elucidated the responses of Physcomitrella patens to two different VSRs. We transformed P. patens plants to express VSR P19 from tomato bushy stunt virus and VSR 2b from cucumber mosaic virus, respectively. RNA sequencing and quantitative PCR were used to detect the effects of VSRs on gene expression. Small RNA (sRNA) sequencing was used to estimate the influences of VSRs on the sRNA pool of P. patens. Expression of either VSR-encoding gene caused developmental disorders in P. patens. The transcripts of four different transcription factors (AP2/erf, EREB-11 and two MYBs) accumulated in the P19 lines. sRNA sequencing revealed that VSR P19 significantly changed the microRNA pool in P. patens. Our results suggest that VSR P19 is functional in P. patens and affects the abundance of specific microRNAs interfering with gene expression. The results open new opportunities for using Physcomitrella as an alternative system to study plant-virus interactions.

The regulatory activities of microRNAs in non-vascular plants: a mini review

MicroRNA-mediated gene regulation in non-vascular plants is potentially involved in several unique biological functions, including biosynthesis of several highly valuable exclusive bioactive compounds, and those small RNAs could be manipulated for the overproduction of essential bioactive compounds in the future. MicroRNAs (miRNAs) are a class of endogenous, small (20-24 nucleotides), non-coding RNA molecules that regulate gene expression through the miRNA-mediated mechanisms of either translational inhibition or messenger RNA (mRNA) cleavage. In the past years, studies have mainly focused on elucidating the roles of miRNAs in vascular plants as compared to non-vascular plants. However, non-vascular plant miRNAs have been predicted to be involved in a wide variety of specific biological mechanisms; nevertheless, some of them have been demonstrated explicitly, thus showing that the research field of this plant group owns a noteworthy potential to develop novel investigations oriented towards the functional characterization of these miRNAs. Furthermore, the insights into the roles of miRNAs in non-vascular plants might be of great importance for designing the miRNA-based genetically modified plants for valuable secondary metabolites, active compounds, and biofuels in the future. Therefore, in this current review, we provide an overview of the potential roles of miRNAs in different groups of non-vascular plants such as algae and bryophytes.

Plant-pathogen interactions: MicroRNA-mediated trans-kingdom gene regulation in fungi and their host plants

MicroRNAs (miRNAs) have been prevalently studied in plants, animals, and viruses. However, recent studies show evidences of miRNA-like RNAs (milRNAs) in fungi as well. It is known that after successful infection, pathogens hijack the host machinery and use it for their own growth and multiplication. Alternatively, resistant plants can overcome the pathogen attack by a variety of mechanisms. Based on this prior knowledge, we computationally predicted milRNAs from 13 fungi, and identified their targets in transcriptomes of the respective fungi as well as their host plants. The expressions of the milRNAs and targets were confirmed using qRT-PCR. We found that plant miRNAs targeted fungal virulence genes, while fungal milRNAs targeted plant resistance genes; corroborating miRNA-mediated trans-kingdom gene regulation and the roles of miRNAs in plant-pathogen interactions. Transgenic plants with miRNAs targeting fungal virulence genes, or anti-sense of fungal milRNAs, would be expected to be highly resistant to the fungal pathogens.

Novel and Conserved miRNAs Among Brazilian Pine and Other Gymnosperms

The knowledge about plant miRNAs has increased exponentially, with thousands of miRNAs been reported in different plant taxa using high throughput sequencing technologies and bioinformatic tools. Nevertheless, several groups of plants remain unexplored, and the gap of knowledge about conifer miRNAs is considerable. There is no sequence or functional information available on miRNAs in Araucariaceae. This group is represented in Brazil by only one species, Araucaria angustifolia, an endangered species known as Brazilian pine. In the present study, Brazilian pine has its transcriptome explored with respect to small RNAs, representing the first description in a member of the Araucariaceae family. The screening for conserved miRNAs in Brazilian pine revealed 115 sequences of 30 miRNA families. A total of 106 precursors sequences were predicted. Forty one comprised conserved miRNAs from 16 families, whereas 65 were annotated as novel miRNAs. The comparison of Brazilian pine precursors with sRNA libraries of other five conifer species indicates that 9 out 65 novel miRNAs are conserved among gymnosperms, while 56 seems to be specific for Brazilian pine or restricted to Araucariaceae family. Analysis comparing novel Brazilian pine miRNAs precursors and Araucaria cunninghamii RNA-seq data identified seven orthologs between both species. Mature miRNA identified by bioinformatics predictions were validated using stem-loop RT-qPCR assays. The expression pattern of conserved and novel miRNAs was analyzed in five different tissues of 3-month-old Araucaria seedlings. The present study provides insights about the nature and composition of miRNAs in an Araucariaceae species, with valuable information on miRNAs diversity and conservation in this taxon.

MicroRNAs in Marchantia polymorpha

Contents Summary 409 I. Introduction 409 II. RNA silencing machinery in Marchantia polymorpha 410 III. miRNA prediction by integrating omics approach 410 IV. miRNAs and their targets in Marchantia polymorpha 410 V. Mpo-miR390-mediated MpTAS3 tasiRNA biogenesis and potential tasiARF target MpARF2 414 VI. Artificial miRNA and CRISPR-CAS9 edited MIR gene in Marchantia polymorpha 414 VII. Conclusions 415 Acknowledgements 415 References 415 SUMMARY: The liverwort Marchantia polymorpha occupies an important phylogenetic position for comparative studies of land plant gene regulation. Multiple gene regulatory pathways mediated by small RNAs, including microRNAs (miRNAs), trans-acting short-interfering RNAs, and heterochromatic siRNAs often associated with RNA-dependent DNA methylation, have been characterized in flowering plants. Genes for essential components for all of these small RNA-mediated gene silencing pathways are found in M. polymorpha as well as the moss Phsycomitrella patens, indicating that these pathways existed in the ancestral land plant. However, only seven miRNAs are conserved across land plants, with both ancestral and novel targets identified in M. polymorpha. There is little or no evidence that any of these conserved miRNAs are present in algae. As with other plants investigated, most miRNAs in M. polypmorpha exhibit lineage-specific evolution. Application of artificial miRNA and CRISPR-Cas9 technologies in genetic studies of M. polymorpha provide avenues to further investigate miRNA biology.

Small but powerful: function of microRNAs in plant development

MicroRNAs (miRNAs) are a group of endogenous noncoding small RNAs frequently 21 nucleotides long. miRNAs act as negative regulators of their target genes through sequence-specific mRNA cleavage, translational repression, or chromatin modifications. Alterations of the expression of a miRNA or its targets often result in a variety of morphological and physiological abnormalities, suggesting the strong impact of miRNAs on plant development. Here, we review the recent advances on the functional studies of plant miRNAs. We will summarize the regulatory networks of miRNAs in a series of developmental processes, including meristem development, establishment of lateral organ polarity and boundaries, vegetative and reproductive organ growth, etc. We will also conclude the conserved and species-specific roles of plant miRNAs in evolution and discuss the strategies for further elucidating the functional mechanisms of miRNAs during plant development.

Sequencing of small RNAs of the fern Pleopeltis minima (Polypodiaceae) offers insight into the evolution of the microrna repertoire in land plants

MicroRNAs (miRNAs) are short, single stranded RNA molecules that regulate the stability and translation of messenger RNAs in diverse eukaryotic groups. Several miRNA genes are of ancient origin and have been maintained in the genomes of animal and plant taxa for hundreds of millions of years, playing key roles in development and physiology. In the last decade, genome and small RNA (sRNA) sequencing of several plant species have helped unveil the evolutionary history of land plants. Among these, the fern group (monilophytes) occupies a key phylogenetic position, as it represents the closest extant cousin taxon of seed plants, i.e. gymno- and angiosperms. However, in spite of their evolutionary, economic and ecological importance, no fern genome has been sequenced yet and few genomic resources are available for this group. Here, we sequenced the small RNA fraction of an epiphytic South American fern, Pleopeltis minima (Polypodiaceae), and compared it to plant miRNA databases, allowing for the identification of miRNA families that are shared by all land plants, shared by all vascular plants (tracheophytes) or shared by euphyllophytes (ferns and seed plants) only. Using the recently described transcriptome of another fern, Lygodium japonicum, we also estimated the degree of conservation of fern miRNA targets in relation to other plant groups. Our results pinpoint the origin of several miRNA families in the land plant evolutionary tree with more precision and are a resource for future genomic and functional studies of fern miRNAs.

Deep-sequence profiling of miRNAs and their target prediction in Monotropa hypopitys

Myco-heterotroph Monotropa hypopitys is a widely spread perennial herb used to study symbiotic interactions and physiological mechanisms underlying the development of non-photosynthetic plant. Here, we performed, for the first time, transcriptome-wide characterization of M. hypopitys miRNA profile using high throughput Illumina sequencing. As a result of small RNA library sequencing and bioinformatic analysis, we identified 55 members belonging to 40 families of known miRNAs and 17 putative novel miRNAs unique for M. hypopitys. Computational screening revealed 206 potential mRNA targets for known miRNAs and 31 potential mRNA targets for novel miRNAs. The predicted target genes were described in Gene Ontology terms and were found to be involved in a broad range of metabolic and regulatory pathways. The identification of novel M. hypopitys-specific miRNAs, some with few target genes and low abundances, suggests their recent evolutionary origin and participation in highly specialized regulatory mechanisms fundamental for non-photosynthetic biology of M. hypopitys. This global analysis of miRNAs and their potential targets in M. hypopitys provides a framework for further investigation of miRNA role in the evolution and establishment of non-photosynthetic myco-heterotrophs.

Identification of miRNAs and Their Targets in the Liverwort Marchantia polymorpha by Integrating RNA-Seq and Degradome Analyses

Bryophytes (liverworts, hornworts and mosses) comprise the three earliest diverging lineages of land plants (embryophytes). Marchantia polymorpha, a complex thalloid Marchantiopsida liverwort that has been developed into a model genetic system, occupies a key phylogenetic position. Therefore, M. polymorpha is useful in studies aiming to elucidate the evolution of gene regulation mechanisms in plants. In this study, we used computational, transcriptomic, small RNA and degradome analyses to characterize microRNA (miRNA)-mediated pathways of gene regulation in M. polymorpha. The data have been integrated into the open access ContigViews-miRNA platform for further reference. In addition to core components of the miRNA pathway, 129 unique miRNA sequences, 11 of which could be classified into seven miRNA families that are conserved in embryophytes (miR166a, miR390, miR529c, miR171-3p, miR408a, miR160 and miR319a), were identified. A combination of computational and degradome analyses allowed us to identify and experimentally validate 249 targets. In some cases, the target genes are orthologous to those of other embryophytes, but in other cases, the conserved miRNAs target either paralogs or members of different gene families. In addition, the newly discovered Mpo-miR11707.1 and Mpo-miR11707.2 are generated from a common precursor and target MpARGONAUTE1 (LW1759). Two other newly discovered miRNAs, Mpo-miR11687.1 and Mpo-miR11681.1, target the MADS-box transcription factors MpMADS1 and MpMADS2, respectively. Interestingly, one of the pentatricopeptide repeat (PPR) gene family members, MpPPR_66 (LW9825), the protein products of which are generally involved in various steps of RNA metabolism, has a long stem-loop transcript that can generate Mpo-miR11692.1 to autoregulate MpPPR_66 (LW9825) mRNA. This study provides a foundation for further investigations of the RNA-mediated silencing mechanism in M. polymorpha as well as of the evolution of this gene silencing pathway in embryophytes.

Identification and expression profile analysis of NUCLEAR FACTOR-Y families in Physcomitrella patens

NUCLEAR FACTOR Y transcription factors belong to a multimember family and consist of NF-YA/B/C subunits. Members of the NF-Y family have been reported to regulate physiological processes in plant. In this study, we identified and annotated two NF-YA, nine NF-B, and twelve NF-YC proteins in the genome of Physcomitrella patens. Analyses of conserved domains demonstrated that PpNF-YA/B/C shared the same conserved domains with their orthologous proteins in Arabidopsis, O. sativa and mouse. Expression profiles indicated that PpNF-Ys were widely expressed in different tissues and developmental stages of P. patens throughout protonema and gametophores. The majority of PpNF-Y genes were responsive to abiotic stress via either ABA-independent or -dependent pathways. Some of ABA-regulated PpNF-Y expression were mediated by ABI3. To our knowledge, this study was the first to evaluate NF-Y families in Physcomitrella patens, and provides a foundation to dissect the function of PpNF-Ys.

Biogenesis, evolution, and functions of plant microRNAs

This review focuses on the biological role of one class of plant small RNAs, ~22-nt microRNAs (miRNAs). The majority of plant miRNA targets are genes encoding the effector factors of cell signaling pathways. The regulation of their expression is necessary for both ontogenesis and rapid response of plants to biotic and abiotic stress factors. We also summarized current views on the biogenesis and evolution of plant miRNAs as well as the techniques used for their investigation.

Peculiar evolutionary history of miR390-guided TAS3-like genes in land plants

PCR-based approach was used as a phylogenetic profiling tool to probe genomic DNA samples from representatives of evolutionary distant moss taxa, namely, classes Bryopsida, Tetraphidopsida, Polytrichopsida, Andreaeopsida, and Sphagnopsida. We found relatives of all Physcomitrella patens miR390 and TAS3-like loci in these plant taxa excluding Sphagnopsida. Importantly, cloning and sequencing of Marchantia polymorpha genomic DNA showed miR390 and TAS3-like sequences which were also found among genomic reads of M. polymorpha at NCBI database. Our data suggest that the ancient plant miR390-dependent TAS molecular machinery firstly evolved to target AP2-like mRNAs in Marchantiophyta and only then both ARF- and AP2-specific mRNAs in mosses. The presented analysis shows that moss TAS3 families may undergone losses of tasiAP2 sites during evolution toward ferns and seed plants. These data confirm that miR390-guided genes coding for ARF- and AP2-specific ta-siRNAs have been gradually changed during land plant evolution.

Identification of novel microRNAs in Hevea brasiliensis and computational prediction of their targets

BACKGROUND

Plants respond to external stimuli through fine regulation of gene expression partially ensured by small RNAs. Of these, microRNAs (miRNAs) play a crucial role. They negatively regulate gene expression by targeting the cleavage or translational inhibition of target messenger RNAs (mRNAs). In Hevea brasiliensis, environmental and harvesting stresses are known to affect natural rubber production. This study set out to identify abiotic stress-related miRNAs in Hevea using next-generation sequencing and bioinformatic analysis.

RESULTS

Deep sequencing of small RNAs was carried out on plantlets subjected to severe abiotic stress using the Solexa technique. By combining the LeARN pipeline, data from the Plant microRNA database (PMRD) and Hevea EST sequences, we identified 48 conserved miRNA families already characterized in other plant species, and 10 putatively novel miRNA families. The results showed the most abundant size for miRNAs to be 24 nucleotides, except for seven families. Several MIR genes produced both 20-22 nucleotides and 23-27 nucleotides. The two miRNA class sizes were detected for both conserved and putative novel miRNA families, suggesting their functional duality. The EST databases were scanned with conserved and novel miRNA sequences. MiRNA targets were computationally predicted and analysed. The predicted targets involved in "responses to stimuli" and to "antioxidant" and "transcription activities" are presented.

CONCLUSIONS

Deep sequencing of small RNAs combined with transcriptomic data is a powerful tool for identifying conserved and novel miRNAs when the complete genome is not yet available. Our study provided additional information for evolutionary studies and revealed potentially specific regulation of the control of redox status in Hevea.