MicroRNA gene evolution in Arabidopsis lyrata and Arabidopsis thaliana

Principles of miRNA/miRNA* function in plant MIRNA processing

MicroRNAs (miRNAs) are essential regulators of gene expression, defined by their unique biogenesis, which requires the precise excision of the small RNA from an imperfect fold-back precursor. Unlike their animal counterparts, plant miRNA precursors exhibit variations in sizes and shapes. Plant MIRNAs can undergo processing in a base-to-loop or loop-to-base direction, with DICER-LIKE1 (DCL1) releasing the miRNA after two cuts (two-step MIRNAs) or more (sequential MIRNAs). In this study, we demonstrate the critical role of the miRNA/miRNA* duplex region in the processing of miRNA precursors. We observed that endogenous MIRNAs frequently experience suboptimal processing in vivo due to mismatches in the miRNA/miRNA* duplex, a key region that fine-tunes miRNA levels. Enhancing the interaction energy of the miRNA/miRNA* duplex in two-step MIRNAs results in a substantial increase in miRNA levels. Conversely, sequential MIRNAs display distinct and specific requirements for the miRNA/miRNA* duplexes along their foldback structure. Our work establishes a connection between the miRNA/miRNA* structure and precursor processing mechanisms. Furthermore, we reveal a link between the biological function of miRNAs and the processing mechanism of their precursors with the evolution of plant miRNA/miRNA* duplex structures.

The floral development of the allotetraploid Coffea arabica L. correlates with a small RNA dynamic reprogramming

Noncoding and coding RNAs are key regulators of plant growth, development, and stress responses. To investigate the types of transcripts accumulated during the vegetative to reproductive transition and floral development in the Coffea arabica L., we sequenced small RNA libraries from eight developmental stages, up to anthesis. We combined these data with messenger RNA and PARE sequencing of two important development stages that marks the transition of an apparent latent to a rapid growth stage. In addition, we took advantage of multiple in silico tools to characterize genomic loci producing small RNAs such as phasiRNAs, miRNAs, and tRFs. Our differential and co-expression analysis showed that some types of small RNAs such as tRNAs, snoRNAs, snRNAs, and phasiRNAs preferentially accumulate in a stage-specific manner. Members of the miR482/miR2118 superfamily and their 21-nucleotide phasiRNAs originating from resistance genes show a robust co-expression pattern that is maintained across all the evaluated developmental stages. Finally, the majority of miRNAs accumulate in a family stage-specific manner, related to modulated hormonal responses and transcription factor expression.

Role of microRNAs and their putative mechanism in regulating potato (Solanum tuberosum L.) life cycle and response to various environmental stresses

The exponentially increasing population and the demand for food is inextricably linked. This has shifted global attention to improving crop plant traits to meet global food demands. Potato (Solanum tuberosum L.) is a major non-grain food crop that is grown all over the world. Currently, some of the major global potato research work focuses on the significance of microRNAs (miRNAs) in potato. miRNAs are a type of non-coding RNAs that regulate the gene expression of their target mRNA genes by cleavage and/or their translational inhibition. This suggests an essential role of miRNAs in a multitude of plant biological processes, including maintenance of genome integrity, plant growth, development and maturation, and initiation of responses to various stress conditions. Therefore, engineering miRNAs to generate stress-resistant varieties of potato may result in high yield and improved nutritional qualities. In this review, we discuss the potato miRNAs specifically known to play an essential role in the various stages of the potato life cycle, conferring stress-resistant characteristics, and modifying gene expression. This review highlights the significance of the miRNA machinery in plants, especially potato, encouraging further research into engineering miRNAs to boost crop yields and tolerance towards stress.

Synonymous sites for accessibility around microRNA binding sites in bacterial spot and speck disease resistance genes of tomato

The major causes of mass tomato infections in both covered and open ground are agents of bacterial spot and bacterial speck diseases. MicroRNAs (miRNAs) are 16-21 nucleotides in length, non-coding RNAs that inhibit translation and trigger mRNA degradation. MiRNAs play a significant part in plant resistance to abiotic and biotic stresses by mediating gene regulation via post-transcriptional RNA silencing. In this study, we analyzed a collection of bacterial resistance genes of tomato and their binding sites for tomato miRNAs and Pseudomonas syringe pv. tomato miRNAs. Our study found that two genes, bacterial spot disease resistance gene (Bs4) and bacterial speck disease resistance gene (Prf), have a 7mer-m8 perfect seed match with miRNAs. Bs4 was targeted by one tomato miRNA (sly-miR9470-3p) and three Pseudomonas syringe pv. tomato miRNAs (PSTJ4_3p_27246, PSTJ4_3p_27246, and PSTJ4_3p_27246). Again, Prf gene was found to be targeted by two tomato miRNAs namely, sly-miR9469-5p and sly-miR9474-3p. The accessibility of the miRNA-target site and its flanking regions and the relationship between relative synonymous codon usage and tRNAs were compared. Strong access to miRNA targeting regions and decreased rate of translations suggested that miRNAs might be efficient in binding to their particular targets. We also found the existence of rare codons, which suggests that it could enhance miRNA targeting even more. The codon usage pattern analysis of the two genes revealed that both were AT-rich (Bs4 = 63.2%; Prf = 60.8%). We found a low codon usage bias in both genes, suggesting that selective restriction might regulate them. The silencing property of miRNAs would allow researchers to discover the involvement of plant miRNAs in pathogen invasion. However, the efficient validation of direct targets of miRNAs is an urgent need that might be highly beneficial in enhancing plant resistance to multiple pathogenic diseases.

Whole-transcriptome analyses identify key differentially expressed mRNAs, lncRNAs, and miRNAs associated with male sterility in watermelon

Male sterility is a valuable trait for watermelon breeding, as watermelon hybrids exhibit obvious heterosis. However, the underlying regulatory mechanism is still largely unknown, especially regarding the related non-coding genes. In the present study, approximately 1035 differentially expressed genes (DEGs), as well as 80 DE-lncRNAs and 10 DE-miRNAs, were identified, with the overwhelming majority down-regulated in male-sterile floral buds. Enrichment analyses revealed that the general phenylpropanoid pathway as well as its related metabolisms was predicted to be altered in a mutant compared to its fertile progenitor. Meanwhile, the conserved genetic pathway DYT1-TDF1-AMS-MS188-MS1, as well as the causal gene ClAMT1 for the male-sterile mutant Se18, was substantially disrupted during male reproductive development. In addition, some targets of the key regulators AMS and MS188 in tapetum development were also down-regulated at a transcriptional level, such as ABCG26 (Cla004479), ACOS5 (Cla022956), CYP703A2 (Cla021151), PKSA (Cla021099), and TKPR1 (Cla002563). Considering lncRNAs may act as functional endogenous target mimics of miRNAs, competitive endogenous RNA networks were subsequently constructed, with the most complex one containing three DE-miRNAs, two DE-lncRNAs, and 21 DEGs. Collectively, these findings not only contribute to a better understanding of genetic regulatory networks underlying male sterility in watermelon, but also provide valuable candidates for future research.

Miniature Inverted-Repeat Transposable Elements: Small DNA Transposons That Have Contributed to Plant MICRORNA Gene Evolution

Angiosperms form the largest phylum within the Plantae kingdom and show remarkable genetic variation due to the considerable difference in the nuclear genome size of each species. Transposable elements (TEs), mobile DNA sequences that can amplify and change their chromosome position, account for much of the difference in nuclear genome size between individual angiosperm species. Considering the dramatic consequences of TE movement, including the complete loss of gene function, it is unsurprising that the angiosperms have developed elegant molecular strategies to control TE amplification and movement. Specifically, the RNA-directed DNA methylation (RdDM) pathway, directed by the repeat-associated small-interfering RNA (rasiRNA) class of small regulatory RNA, forms the primary line of defense to control TE activity in the angiosperms. However, the miniature inverted-repeat transposable element (MITE) species of TE has at times avoided the repressive effects imposed by the rasiRNA-directed RdDM pathway. MITE proliferation in angiosperm nuclear genomes is due to their preference to transpose within gene-rich regions, a pattern of transposition that has enabled MITEs to gain further transcriptional activity. The sequence-based properties of a MITE results in the synthesis of a noncoding RNA (ncRNA), which, after transcription, folds to form a structure that closely resembles those of the precursor transcripts of the microRNA (miRNA) class of small regulatory RNA. This shared folding structure results in a MITE-derived miRNA being processed from the MITE-transcribed ncRNA, and post-maturation, the MITE-derived miRNA can be used by the core protein machinery of the miRNA pathway to regulate the expression of protein-coding genes that harbor homologous MITE insertions. Here, we outline the considerable contribution that the MITE species of TE have made to expanding the miRNA repertoire of the angiosperms.

Clustering pattern and evolution characteristic of microRNAs in grass carp (Ctenopharyngodon idella)

BACKGROUND

A considerable fraction of microRNAs (miRNAs) are highly conserved, and certain miRNAs correspond to genomic clusters. The clustering of miRNAs can be advantageous, possibly by allowing coordinated expression. However, little is known about the evolutionary forces responsible for the loss and acquisition of miRNA and miRNA clusters.

RESULTS

The results demonstrated that several novel miRNAs arose throughout grass carp evolution. Duplication and de novo production were critical strategies for miRNA cluster formation. Duplicates accounted for a smaller fraction of the expansion in the grass carp miRNA than de novo creation. Clustered miRNAs are more conserved and change slower, whereas unique miRNAs usually have high evolution rates and low expression levels. The expression level of miRNA expression in clusters is strongly correlated.

CONCLUSIONS

This study examines the genomic distribution, evolutionary background, and expression regulation of grass carp miRNAs. Our findings provide novel insights into the genesis and development of miRNA clusters in teleost.

Arabidopsis AAR2, a conserved splicing factor in eukaryotes, acts in microRNA biogenesis

In yeast and humans, AAR2 is involved in pre-messenger RNA (pre-mRNA) splicing through regulating U5 snRNP assembly. This study shows that Arabidopsis AAR2 promotes microRNA (miRNA) accumulation in addition to its conserved role in pre-mRNA splicing. AAR2 is associated with the microprocessor component HYL1 and promotes its dephosphorylation to produce the active form in miRNA biogenesis. The study also reveals a previously unknown role of HYL1 in causing the degradation of the primary precursors to miRNAs (pri-miRNAs) and a role of AAR2 in protecting pri-miRNAs from HYL1-depedent degradation. Taken together, our findings provide insights into the role of a conserved splicing factor in miRNA biogenesis in plants.

MicroRNAs (miRNAs) play an essential role in plant growth and development, and as such, their biogenesis is fine-tuned via regulation of the core microprocessor components. Here, we report that Arabidopsis AAR2, a homolog of a U5 snRNP assembly factor in yeast and humans, not only acts in splicing but also promotes miRNA biogenesis. AAR2 interacts with the microprocessor component hyponastic leaves 1 (HYL1) in the cytoplasm, nucleus, and dicing bodies. In aar2 mutants, abundance of nonphosphorylated HYL1, the active form of HYL1, and the number of HYL1-labeled dicing bodies are reduced. Primary miRNA (pri-miRNA) accumulation is compromised despite normal promoter activities of MIR genes in aar2 mutants. RNA decay assays show that the aar2-1 mutation leads to faster degradation of pri-miRNAs in a HYL1-dependent manner, which reveals a previously unknown and negative role of HYL1 in miRNA biogenesis. Taken together, our findings reveal a dual role of AAR2 in miRNA biogenesis and pre-messenger RNA splicing.

Parasitic plant small RNA analyses unveil parasite-specific signatures of microRNA retention, loss, and gain

Parasitism is a successful life strategy that has evolved independently in several families of vascular plants. The genera Cuscuta and Orobanche represent examples of the two profoundly different groups of parasites: one parasitizing host shoots and the other infecting host roots. In this study, we sequenced and described the overall repertoire of small RNAs from Cuscuta campestris and Orobanche aegyptiaca. We showed that C. campestris contains a number of novel microRNAs (miRNAs) in addition to a conspicuous retention of miRNAs that are typically lacking in other Solanales, while several typically conserved miRNAs seem to have become obsolete in the parasite. One new miRNA appears to be derived from a horizontal gene transfer event. The exploratory analysis of the miRNA population (exploratory due to the absence of a full genomic sequence for reference) from the root parasitic O. aegyptiaca also revealed a loss of a number of miRNAs compared to photosynthetic species from the same order. In summary, our study shows partly similar evolutionary signatures in the RNA silencing machinery in both parasites. Our data bear proof for the dynamism of this regulatory mechanism in parasitic plants.

Deciphering the role of miRNA in reprogramming plant responses to drought stress

Drought is the most prevalent environmental stress that affects plants' growth, development, and crop productivity. However, plants have evolved adaptive mechanisms to respond to the harmful effects of drought. They reprogram their: transcriptome, proteome, and metabolome that alter their cellular and physiological processes and establish cellular homeostasis. One of the crucial regulatory processes that govern this reprogramming is post-transcriptional regulation by microRNAs (miRNAs). miRNAs are small non-coding RNAs, involved in the downregulation of the target mRNA via translation inhibition/mRNA degradation/miRNA-mediated mRNA decay/ribosome drop off/DNA methylation. Many drought-inducible miRNAs have been identified and characterized in plants. Their main targets are regulatory genes that influence growth, development, osmotic stress tolerance, antioxidant defense, phytohormone-mediated signaling, and delayed senescence during drought stress. Overexpression of drought-responsive miRNAs (Osa-miR535, miR160, miR408, Osa-miR393, Osa-miR319, and Gma-miR394) in certain plants has led to tolerance against drought stress indicating their vital role in stress mitigation. Similarly, knock down (miR166/miR398c) or deletion (miR169 and miR827) of miRNAs has also resulted in tolerance to drought stress. Likewise, engineered Arabidopsis plants with miR165, miR166 using short tandem target mimic strategy, exhibited drought tolerance. Since miRNAs regulate the expression of an array of drought-responsive genes, they can act as prospective targets for genetic manipulations to enhance drought tolerance in crops and achieve sustainable agriculture. Further investigations toward functional characterization of diverse miRNAs, and understanding stress-responses regulated by these miRNAs and their utilization in biotechnological applications is highly recommended.

Mechanism for the genomic and functional evolution of the MIR2118 family in the grass lineage

The MIR2118 family has undergone tremendous expansion in the grass lineage, in which the miRNA targets numerous noncoding PHAS loci to produce 21-nt phased small interfering RNAs (phasiRNAs) involved in male fertility. However, the evolutionary trajectory of the grass MIR2118 genes and the functions of phasiRNAs have not yet been fully elucidated. We conducted comparative genomic, molecular evolution, expression and parallel analysis of RNA ends (PARE) analyses of MIR2118 and the miR2118-mediated regulatory pathway in grasses, focusing on Oryza sativa. In total, 617 MIR2118 and eight MIR1859 novel members were identified. Phylogenetic analyses showed that grass MIR2118 genes form a distinct clade from the MIR482/2118 genes of nongrass species. We reconstructed hypothetical evolutionary histories of the grass MIR2118 clusters and its MIR1859 variants, and examined the polycistronic composition and the differential expression of the osa-MIR2118 clusters. PARE data showed that osa-miR2118 might also direct the cleavage of some protein-coding gene transcripts. Importantly, we found that PARE analysis is inherently prone to false-positive target predictions when a large number of small RNAs, such as phasiRNAs, are analysed. Our results revealed the evolution and diversification of the MIR2118 family, and provide new insights into the functions of phasiRNAs in the grasses.

MicroRNAs Are Involved in Regulating Plant Development and Stress Response through Fine-Tuning of TIR1/AFB-Dependent Auxin Signaling

Auxin, primarily indole-3-acetic acid (IAA), is a versatile signal molecule that regulates many aspects of plant growth, development, and stress response. Recently, microRNAs (miRNAs), a type of short non-coding RNA, have emerged as master regulators of the auxin response pathways by affecting auxin homeostasis and perception in plants. The combination of these miRNAs and the autoregulation of the auxin signaling pathways, as well as the interaction with other hormones, creates a regulatory network that controls the level of auxin perception and signal transduction to maintain signaling homeostasis. In this review, we will detail the miRNAs involved in auxin signaling to illustrate its in planta complex regulation.

A comprehensive overview of miRNA targeting drought stress resistance in plants

MicroRNAs (miRNAs) are essential nonprotein-coding genes. In a range of organisms, miRNAs has been reported to play an essential role in regulating gene expressions at post-transcriptional level. They participate in most of the stress responsive processes in plants. Drought is an ultimate abiotic stress that affects the crop production. Therefore understanding drought stress responses are essential to improve the production of agricultural crops. Throughout evolution, plants have developed their own defense systems to cope with the adversities of environmental stresses. Among defensive mechanisms include the regulations of gene expression by miRNAs. Drought stress regulates the expression of some of the functionally conserved miRNAs in different plants. The given properties of miRNAs provide an insight to genetic alterations and enhancing drought resistance in cereal crops. The current review gives a summary to regulatory mechanisms in plants as well as miRNAs response to drought stresses in cereal crops. Some possible approaches and guidelines for the exploitation of drought stress miRNA responses to improve cereal crops are also described.

Biotic stress-associated microRNA families in plants

Plants and animals utilize various regulatory mechanisms for control of gene expression during development in different tissues and cell types. About 30 years ago, a new mechanism of gene regulation, termed RNA interference (RNAi), was discovered and proved revolutionary for the mechanistic understanding of gene regulation. Noncoding RNAs, including short, 21-24 nucleotide (nt) long microRNAs (miRNAs), endogenously-generated from MIR genes, are key components of RNAi processes, by post-transcriptionally controlling transcripts with antisense complementarity through either translational repression or mRNA degradation. Since their discovery, important roles in regulation of ontogenetic development, cell differentiation, proliferation, and apoptosis in eukaryotes have been elucidated. In plants, miRNAs are known regulatory elements of basic endogenous functions and responses to the environmental stimuli. While the role of miRNAs in regulation of nutrient uptake, circadian clock and general response to abiotic stress is already well understood, a comprehensive understanding of their immune-regulatory roles in response to various biotic stress factors has not yet been achieved. This review summarizes the current understanding of the function of miRNAs and their targets in plants during interaction with microbial pathogens and symbionts. Additionally, we provide a consensus conclusion regarding the typical induction or repression response of conserved miRNA families to pathogenic and beneficial fungi, bacteria, and oomycetes, as well as an outlook of agronomic application of miRNAs in plants. Further investigation of plant miRNAs responsive to microbes, aided with novel sequencing and bioinformatics approaches for discovery and prediction in non-model organisms holds great potential for development of new forms of plant protection.

Role of SNPs in the Biogenesis of Mature miRNAs

Single nucleotide polymorphisms (SNPs) play a significant role in microRNA (miRNA) generation, processing, and function and contribute to multiple phenotypes and diseases. Therefore, whole-genome analysis of how SNPs affect miRNA maturation mechanisms is important for precision medicine. The present study established an SNP-associated pre-miRNA (SNP-pre-miRNA) database, named miRSNPBase, and constructed SNP-pre-miRNA sequences. We also identified phenotypes and disease biomarker-associated isoform miRNA (isomiR) based on miRFind, which was developed in our previous study. We identified functional SNPs and isomiRs. We analyzed the biological characteristics of functional SNPs and isomiRs and studied their distribution in different ethnic groups using whole-genome analysis. Notably, we used individuals from Great Britain (GBR) as examples and identified isomiRs and isomiR-associated SNPs (iso-SNPs). We performed sequence alignments of isomiRs and miRNA sequencing data to verify the identified isomiRs and further revealed GBR ethnographic epigenetic dominant biomarkers. The SNP-pre-miRNA database consisted of 886 pre-miRNAs and 2640 SNPs. We analyzed the effects of SNP type, SNP location, and SNP-mediated free energy change during mature miRNA biogenesis and found that these factors were closely associated to mature miRNA biogenesis. Remarkably, 158 isomiRs were verified in the miRNA sequencing data for the 18 GBR samples. Our results indicated that SNPs affected the mature miRNA processing mechanism and contributed to the production of isomiRs. This mechanism may have important significance for epigenetic changes and diseases.

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.

Conservation and turnover of miRNAs and their highly complementary targets in early branching animals

MicroRNAs (miRNAs) are crucial post-transcriptional regulators that have been extensively studied in Bilateria, a group comprising the majority of extant animals, where more than 30 conserved miRNA families have been identified. By contrast, bilaterian miRNA targets are largely not conserved. Cnidaria is the sister group to Bilateria and thus provides a unique opportunity for comparative studies. Strikingly, like their plant counterparts, cnidarian miRNAs have been shown to predominantly have highly complementary targets leading to transcript cleavage by Argonaute proteins. Here, we assess the conservation of miRNAs and their targets by small RNA sequencing followed by miRNA target prediction in eight species of Anthozoa (sea anemones and corals), the earliest-branching cnidarian class. We uncover dozens of novel miRNAs but only a few conserved ones. Further, given their high complementarity, we were able to computationally identify miRNA targets in each species. Besides evidence for conservation of specific miRNA target sites, which are maintained between sea anemones and stony corals across 500 Myr of evolution, we also find indications for convergent evolution of target regulation by different miRNAs. Our data indicate that cnidarians have only few conserved miRNAs and corresponding targets, despite their high complementarity, suggesting a high evolutionary turnover.

Contributions of Adaptive Plant Architecture to Transgressive Salinity Tolerance in Recombinant Inbred Lines of Rice: Molecular Mechanisms Based on Transcriptional Networks

Genetic novelties are important nucleators of adaptive speciation. Transgressive segregation is a major mechanism that creates genetic novelties with morphological and developmental attributes that confer adaptive advantages in certain environments. This study examined the morpho-developmental and physiological profiles of recombinant inbred lines (RILs) from the salt-sensitive IR29 and salt-tolerant Pokkali rice, representing the total range of salt tolerance including the outliers at both ends of the spectrum. Morpho-developmental and physiological profiles were integrated with a hypothesis-driven interrogation of mRNA and miRNA transcriptomes to uncover the critical genetic networks that have been rewired for novel adaptive architecture. The transgressive super-tolerant FL510 had a characteristic small tiller angle and wider, more erect, sturdier, and darker green leaves. This unique morphology resulted in lower transpiration rate, which also conferred a special ability to retain water more efficiently for osmotic avoidance. The unique ability for water retention conferred by such adaptive morphology appeared to enhance the efficacy of defenses mediated by Na+ exclusion mechanism (SalTol-effects) inherited from Pokkali. The super-tolerant FL510 and super-sensitive FL499 had the smallest proportions of differentially expressed genes with little overlaps. Genes that were steadily upregulated in FL510 comprised a putative cytokinin-regulated genetic network that appeared to maintain robust growth under salt stress through well-orchestrated cell wall biogenesis and cell expansion, likely through major regulatory (OsRR23, OsHK5) and biosynthetic (OsIPT9) genes in the cytokinin signaling pathway. Meanwhile, a constitutively expressed cluster in FL510 prominently featured two transcription factors (OsIBH1, TAC3) that control tiller angle and growth habit through the brassinosteroid signaling pathway. Both the putative cytokinin-mediated and brassinosteroid-mediated clusters appeared to function as highly coordinated network synergies in FL510. In contrast, both networks appeared to be sub-optimal and inferior in the other RILs and parents as they were disjointed and highly fragmented. Transgressively expressed miRNAs (miR169, miR397, miR827) were also identified as prominent signatures of FL510, with functional implications to mechanisms that support robust growth, homeostasis, and osmotic stress avoidance. Results of this study demonstrate how genetic recombination creates novel morphology that complements inducible defenses hence transgressive adaptive phenotypes.

Host-derived artificial miRNA-mediated silencing of ecdysone receptor gene provides enhanced resistance to Helicoverpa armigera in tomato

Helicoverpa armigera causes huge crop losses due to its polyphagous nature. The present study demonstrates the use of artificial microRNA (amiRNA) mediated gene silencing approach to generate insect resistant tomato plants. Ecdysone receptor (HaEcR) gene of the target pest, H. armigera, which is involved in the regulation of all developmental stages of the insect life cycle, was silenced by sequence-specific amiRNA (amiRNA-HaEcR). Continuous feeding on detached tomato leaves expressing the amiRNA-319a-HaEcR resulted in reduced target gene transcripts and affected the overall growth and survival of H. armigera. Not only the target gene was down-regulated but, the feeding also affected the expression of down-stream genes involved in the ecdysone signaling pathway. The resistant trait was also observed in T₁ generation of tomato transgenic lines. These results further established the role of EcR as a master regulator in insect development and effectiveness of amiRNA technology for efficient control of H. armigera.

The Apple microR171i-SCARECROW-LIKE PROTEINS26.1 Module Enhances Drought Stress Tolerance by Integrating Ascorbic Acid Metabolism

Drought stress severely restricts crop yield and quality. Small noncoding RNAs play critical roles in plant growth, development, and stress responses by regulating target gene expression, but their roles in drought stress tolerance in apple (Malus domestica) are poorly understood. Here, we identified various small noncoding RNAs and their targets from the wild apple species Malus sieversii via high-throughput sequencing and degradome analysis. Forty known microRNAs (miRNAs) and eight new small noncoding RNAs were differentially expressed in response to 2 or 4 h of drought stress treatment. We experimentally verified the expression patterns of five selected miRNAs and their targets. We established that one miRNA, mdm-miR171i, specifically targeted and degraded SCARECROW-LIKE PROTEINS26 1 (MsSCL26 1) transcripts. Both knockout of mdm-miR171i and overexpression of MsSCL26 1 improved drought stress tolerance in the cultivated apple line 'GL-3' by regulating the expression of antioxidant enzyme genes, especially that of MONODEHYDROASCORBATE REDUCTASE, which functions in metabolism under drought stress. Transient expression analysis demonstrated that MsSCL26.1 activates MsMDHAR transcription by positively regulating the activity of the P1 region in its promoter. Therefore, the miR171i-SCL26 1 module enhances drought stress tolerance in apple by regulating antioxidant gene expression and ascorbic acid metabolism.

Genome-wide development of miRNA-based SSR markers in Cleistogenes songorica and analysis of their transferability to Gramineae/non-Gramineae species

Simple sequence repeat (SSR) markers are commonly used for many genetic applications, such as map construction, fingerprinting, and genetic diversity analyses, due to their high reproducibility, polymorphism, and abundance. Endogenous miRNAs play essential roles in plant development and gene expression under diverse biotic and abiotic stress conditions. In the present study, we predicted 110 miRNA-SSR primer pairs from 287 precursor miRNAs. Among 110 primer pairs, 85 were successfully amplified and examined for transferability to other Gramineae and non-Gramineae species. The results showed that all 82 primer pairs yielded unambiguous and strong amplification, and across the 23 studied Cleistogenes accessions, a total of 385 alleles were polymorphic. The number of alleles produced per primer varied from 3 to 11, with an average of 4.69 per locus. The expected heterozygosity (He) ranged from 0.44 to 0.88, with an average of 0.74 per locus, and the PIC (Polymorphism Information Content) values ranged from 0.34 to 0.87, with an average of 0.69 per locus. Furthermore, 1422 miRNA target genes were predicted and analyzed using the GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases. In conclusion, the results showed that an miRNA-based microsatellite marker system can be applicable for genetic diversity and marker-assisted breeding studies.

Genome-Wide microRNA Profiling Using Oligonucleotide Microarray Reveals Regulatory Networks of microRNAs in Nicotiana benthamiana During Beet Necrotic Yellow Vein Virus Infection

Beet necrotic yellow vein virus (BNYVV) infections induce stunting and leaf curling, as well as root and floral developmental defects and leaf senescence in Nicotiana benthamiana. A microarray analysis with probes capable of detecting 1596 candidate microRNAs (miRNAs) was conducted to investigate differentially expressed miRNAs and their targets upon BNYVV infection of N. benthamiana plants. Eight species-specific miRNAs of N. benthamiana were identified. Comprehensive characterization of the N. benthamiana microRNA profile in response to the BNYVV infection revealed that 129 miRNAs were altered, including four species-specific miRNAs. The targets of the differentially expressed miRNAs were predicted accordingly. The expressions of miR164, 160, and 393 were up-regulated by BNYVV infection, and those of their target genes, NAC21/22, ARF17/18, and TIR, were down-regulated. GRF1, which is a target of miR396, was also down-regulated. Further genetic analysis of GRF1, by Tobacco rattle virus-induced gene silencing, assay confirmed the involvement of GRF1 in the symptom development during BNYVV infection. BNYVV infection also induced the up-regulation of miR168 and miR398. The miR398 was predicted to target umecyanin, and silencing of umecyanin could enhance plant resistance against viruses, suggesting the activation of primary defense response to BNYVV infection in N. benthamiana. These results provide a global profile of miRNA changes induced by BNYVV infection and enhance our understanding of the mechanisms underlying BNYVV pathogenesis.

Independent origin of MIRNA genes controlling homologous target genes by partial inverted duplication of antisense-transcribed sequences

Some microRNAs (miRNAs) are key regulators of developmental processes, mainly by controlling the accumulation of transcripts encoding transcription factors that are important for morphogenesis. MADS-box genes encode a family of transcription factors which control diverse developmental processes in flowering plants. Here we study the convergent evolution of two MIRNA (MIR) gene families, named MIR444 and MIR824, targeting members of the same clade of MIKC^C -group MADS-box genes. We show that these two MIR genes most likely originated independently in monocots (MIR444) and in Brassicales (eudicots, MIR824). We provide evidence that, in both cases, the future target gene was transcribed in antisense prior to the evolution of the MIR genes. Both MIR genes then likely originated by a partial inverted duplication of their target genes, resulting in natural antisense organization of the newly evolved MIR gene and its target gene at birth. We thus propose a model for the origin of MIR genes, MEPIDAS (MicroRNA Evolution by Partial Inverted Duplication of Antisense-transcribed Sequences). MEPIDAS is a refinement of the inverted duplication hypothesis. According to MEPIDAS, a MIR gene evolves at a genomic locus at which the future target gene is also transcribed in the antisense direction. A partial inverted duplication at this locus causes the antisense transcript to fold into a stem-loop structure that is recognized by the miRNA biogenesis machinery to produce a miRNA that regulates the gene at this locus. Our analyses exemplify how to elucidate the origin of conserved miRNAs by comparative genomics and will guide future studies. OPEN RESEARCH BADGE: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/genbank/.

Small noncoding RNA discovery and profiling with sRNAtools based on high-throughput sequencing

Small noncoding RNAs (sRNA/sncRNAs) are generated from different genomic loci and play important roles in biological processes, such as cell proliferation and the regulation of gene expression. Next-generation sequencing (NGS) has provided an unprecedented opportunity to discover and quantify diverse kinds of sncRNA, such as tRFs (tRNA-derived small RNA fragments), phasiRNAs (phased, secondary, small-interfering RNAs), Piwi-interacting RNA (piRNAs) and plant-specific 24-nt short interfering RNAs (siRNAs). However, currently available web-based tools do not provide approaches to comprehensively analyze all of these diverse sncRNAs. This study presents a novel integrated platform, sRNAtools (https://bioinformatics.caf.ac.cn/sRNAtools), that can be used in conjunction with high-throughput sequencing to identify and functionally annotate sncRNAs, including profiling microRNAss, piRNAs, tRNAs, small nuclear RNAs, small nucleolar RNAs and rRNAs and discovering isomiRs, tRFs, phasiRNAs and plant-specific 24-nt siRNAs for up to 21 model organisms. Different modules, including single case, batch case, group case and target case, are developed to provide users with flexible ways of studying sncRNA. In addition, sRNAtools supports different ways of uploading small RNA sequencing data in a very interactive queue system, while local versions based on the program package/Docker/virtureBox are also available. We believe that sRNAtools will greatly benefit the scientific community as an integrated tool for studying sncRNAs.

Keep calm and carry on: miRNA biogenesis under stress

MicroRNAs (miRNAs) are major post-transcriptional regulators of gene expression. Their biogenesis relies on the cleavage of longer precursors by a nuclear localized processing machinery. The evolutionary preference of plant miRNAs to silence transcription factors turned these small molecules into key actors during growth and adaptive responses. Furthermore, during their life cycle plants are subject to changes in the environmental conditions surrounding them. In order to face these changes, plants display unique adaptive capacities based on an enormous developmental plasticity, where miRNAs play central roles. Many individual miRNAs have been shown to modulate the plant response to different environmental cues and stresses. In the last few years, increasing evidence has shown that not only individual genes encoding miRNAs but also the miRNA pathway as a whole is subject to regulation in response to external stimulus. In this review, we discuss the current knowledge about the miRNA pathway. We dissect the pathway to analyze the events leading to the generation of these small RNAs and emphasize the regulation of core components of the miRNA biogenesis machinery.

Fine-Tuning of MiR528 Accumulation Modulates Flowering Time in Rice

In plants, microRNA (miRNA) functions in the post-transcriptional repression of target mRNAs have been well explored. However, the mechanisms regulating the accumulation of miRNAs remain poorly understood. Here, we report that distinct mechanisms regulate accumulation of a monocot-specific miRNA, rice (Oryza sativa) miR528. At the transcriptional level, miR528 accumulated to higher levels in older plants than in young seedlings and exhibited aging-modulated gradual accumulation and diurnal rhythms in leaves; at the post-transcriptional level, aging also modulated miR528 levels by enhancing pri-miR528 alternative splicing. We found that miR528 promotes rice flowering under long-day conditions by targeting RED AND FAR-RED INSENSITIVE 2 (OsRFI2). Moreover, natural variations in the MIR528 promoter region caused differences in miR528 expression among rice varieties, which are correlated with their different binding affinities with the transcription factor OsSPL9 that activates the expression of miR528. Taken together, our findings reveal rice plants have evolved sophisticated modes fine-tuning miR528 levels and provide insight into the mechanisms that regulate MIRNA expression in plants.

Unravelling the MicroRNA-Mediated Gene Regulation in Developing Pongamia Seeds by High-Throughput Small RNA Profiling

Pongamia (Millettia pinnata syn. Pongamia pinnata) is a multipurpose biofuel tree which can withstand a variety of abiotic stresses. Commercial applications of Pongamia trees may substantially benefit from improvements in their oil-seed productivity, which is governed by complex regulatory mechanisms underlying seed development. MicroRNAs (miRNAs) are important molecular regulators of plant development, while relatively little is known about their roles in seed development, especially for woody plants. In this study, we identified 236 conserved miRNAs within 49 families and 143 novel miRNAs via deep sequencing of Pongamia seeds sampled at three developmental phases. For these miRNAs, 1327 target genes were computationally predicted. Furthermore, 115 differentially expressed miRNAs (DEmiRs) between successive developmental phases were sorted out. The DEmiR-targeted genes were preferentially enriched in the functional categories associated with DNA damage repair and photosynthesis. The combined analyses of expression profiles for DEmiRs and functional annotations for their target genes revealed the involvements of both conserved and novel miRNA-target modules in Pongamia seed development. Quantitative Real-Time PCR validated the expression changes of 15 DEmiRs as well as the opposite expression changes of six targets. These results provide valuable miRNA candidates for further functional characterization and breeding practice in Pongamia and other oilseed plants.

Characterization on the conservation and diversification of miRNA156 gene family from lower to higher plant species based on phylogenetic analysis at the whole genomic level

miRNA156 family members (miR156s) participate in regulating the transition of plant vegetative and reproductive growth, flower development, and formation of berry skin color by negatively modulating their target gene SPLs. However, the evolution and functional diversification of miR156s in plants remain elusive. Phylogenetic analysis on 310 miR156s from 51 plant species on miRBase 21.0 showed that only miR156a could be conserved in the 51 plant species, but their sequences exhibited variation; another set of miR156s, such as miR156m/n/o/p/q/r/s/t/u/v/w/x/y/z, was identified only in certain special plant species (Glycine max and Malus); also, all base variations in the sequences of 310 miR156s occurred within one miR156 seed sequence, "TGACAGAAGAGAGTGAGCAC," and the changed base sites were mainly located at the 11th and 14th bases from the 5' end of the miR156 seed sequence, in which some base variations of miR156s resulted in a difference in miR156 targeting modes; by contrast, miR156 precursor sequences are highly divergent across diverse species. Similarly, cis-regulatory motifs on the promoter sequence of MIR156s in various plants also exhibited significant discrepancy. The intragenic MIR156 genes overlapped their target SBP genes, thereby suggesting that some microRNAs (miRNAs) originate from duplication of target genes. These traits might be the reasons of the conservation and diversification of miR156 gene family. This study identified the conserved seed sequence "TGACAGAAGAGAGTGAGCAC," and the sequence variation characterization, of miR156 family evolution, also investigated the varied traits of their promoters, precursors, and mature sequences in sequence evolutions and found some miRNAs might originate from duplication of target genes. Our findings will contribute to our understanding of the functional diversification of miRNAs and the interactions of miRNA/target pairs based on the evolutionary history of miRNA genes.

The Elaboration of miRNA Regulation and Gene Regulatory Networks in Plant⁻Microbe Interactions

Plants are exposed to diverse abiotic and biotic stimuli. These require fast and specific integrated responses. Such responses are coordinated at the protein and transcript levels and are incorporated into larger regulatory networks. Here, we focus on the evolution of transcriptional regulatory networks involved in plant–pathogen interactions. We discuss the evolution of regulatory networks and their role in fine-tuning plant defense responses. Based on the observation that many of the cornerstones of immune signaling in angiosperms are also present in streptophyte algae, it is likely that some regulatory components also predate the origin of land plants. The degree of functional conservation of many of these ancient components has not been elucidated. However, ongoing functional analyses in bryophytes show that some components are conserved. Hence, some of these regulatory components and how they are wired may also trace back to the last common ancestor of land plants or earlier. Of course, an understanding of the similarities and differences during the evolution of plant defense networks cannot ignore the lineage-specific coevolution between plants and their pathogens. In this review, we specifically focus on the small RNA regulatory networks involved in fine-tuning of the strength and timing of defense responses and highlight examples of pathogen exploitation of the host RNA silencing system. These examples illustrate well how pathogens frequently target gene regulation and thereby alter immune responses on a larger scale. That this is effective is demonstrated by the diversity of pathogens from distinct kingdoms capable of manipulating the same gene regulatory networks, such as the RNA silencing machinery.

Effects of gene duplication, epistasis, recombination and gene conversion on the fixation time of compensatory mutations

Gene duplication is one of the major mechanisms of molecular evolution. Gene duplication enables copies of a gene to accumulate mutations through functional redundancy. If a gene encodes a specific protein that interacts with other proteins, RNA, or DNA, the relaxation of selective constraints caused by gene duplication might contribute to the fixation of compensatory mutations that occur at the interacting sites. In this study, we investigate the effect of gene duplication, epistasis among the duplicated copies and gene conversion on the fixation time of compensatory mutations by extending the original model of compensatory evolution proposed by Kimura in 1985. Our simulation results reveal that the time to fixation of compensatory mutations can be decreased remarkably by gene duplication if one of the duplicated loci can completely mask the deleterious effects of a mutation that occurs at the other locus. Conversely, the fixation time can be increased by gene duplication if such functional compensation is weak. We also show that the combination of the degree of functional compensation and the rate of gene conversion between duplicate loci have contrasting effects on the time to fixation of compensatory mutations.

Genome-wide comparative analysis in Solanaceous species reveals evolution of microRNAs targeting defense genes in Capsicum spp

MicroRNAs (miRNAs) play roles in various biological processes in plants including growth, development, and disease resistance. Previous studies revealed that some plant miRNAs produce secondary small interfering RNAs (siRNAs) such as phased, secondary siRNAs (phasiRNAs), and they regulate a cascade of gene expression. We performed a genome-wide comparative analysis of miRNAs in Solanaceous species (pepper, tomato, and potato), from an evolutionary perspective. Microsynteny of miRNAs was analysed based on the genomic loci and their flanking genes and most of the well-conserved miRNA genes maintained microsynteny in Solanaceae. We identified target genes of the miRNAs via degradome analysis and found that several miRNAs target many genes encoding nucleotide-binding leucine-rich repeat (NLR) or receptor-like proteins (RLPs), which are known to be major players in defense responses. In addition, disease-resistance-associated miRNAs trigger phasiRNA production in pepper, indicating amplification of the regulation of disease-resistance gene families. Among these, miR-n033a-3p, whose target NLRs have been duplicated in pepper, targets more NLRs belonging to specific subgroup in pepper than those in potato. miRNAs targeting resistance genes might have evolved to regulate numerous targets in Solanaceae, following expansion of target resistance genes. This study provides an insight into evolutionary relationship between miRNAs and their target defense genes in plants.

Plant miRNA Conservation and Evolution

Plant microRNAs do not only perform important roles in development; they also have a fascinating evolutionary dynamics. Their genes appear to originate at quite a high rate during evolution, but most of them evolve initially in an almost neutral way and hence also get lost quite rapidly. Despite the high birth and death rate, a few miRNA-encoding genes got involved in the control of important target genes and thus have been conserved during evolution. This happened obviously at all times and taxonomic levels during land plant evolution. Consequently, the genomes of extant plant species contain a mix of miRNA-encoding genes of different ages, ranging from very young, often even species-specific loci to genes that had already been established in the stem group of extant land plants more than 400 million years ago. It could well be that the evolutionary dynamics of miRNA-encoding genes contributed substantially to the evolution of developmental plasticity in plants.

microRPM: a microRNA prediction model based only on plant small RNA sequencing data

Motivation

MicroRNAs (miRNAs) are endogenous non-coding small RNAs (of about 22 nucleotides), which play an important role in the post-transcriptional regulation of gene expression via either mRNA cleavage or translation inhibition. Several machine learning-based approaches have been developed to identify novel miRNAs from next generation sequencing (NGS) data. Typically, precursor/genomic sequences are required as references for most methods. However, the non-availability of genomic sequences is often a limitation in miRNA discovery in non-model plants. A systematic approach to determine novel miRNAs without reference sequences is thus necessary.

Results

In this study, an effective method was developed to identify miRNAs from non-model plants based only on NGS datasets. The miRNA prediction model was trained with several duplex structure-related features of mature miRNAs and their passenger strands using a support vector machine algorithm. The accuracy of the independent test reached 96.61% and 93.04% for dicots (Arabidopsis) and monocots (rice), respectively. Furthermore, true small RNA sequencing data from orchids was tested in this study. Twenty-one predicted orchid miRNAs were selected and experimentally validated. Significantly, 18 of them were confirmed in the qRT-PCR experiment. This novel approach was also compiled as a user-friendly program called microRPM (miRNA Prediction Model).

Availability and implementation

This resource is freely available at http://microRPM.itps.ncku.edu.tw.

Contact

nslin@sinica.edu.tw or sarah321@mail.ncku.edu.tw.

Supplementary information

Supplementary data are available at Bioinformatics online.

In Silico Identification of Conserved MiRNAs from Physcomitrella patens ESTs and their Target Characterization

Background:

MicroRNAs (miRNAs) are groups of small non-protein-coding endogenous single stranded RNAs with approximately 18-24 nucleotides in length. High evolutionary sequence conservation of miRNAs among plant species and availability of powerful computational tools allow identification of new orthologs and paralogs.

Methods:

New conserved miRNAs in P. patens were found by EST-based homology search approaches. All candidates were screened according to a series of miRNA filtering criteria. Unigene, DFCI Gene Index (PpspGI) databases and psRNATarget algorithm were applied to identify target transcripts using P. patens putative conserved miRNA sequences.

Results:

Nineteen conserved P. patens miRNAs were identified. The sequences were homologous to known reference plant mature miRNA from 10 miRNA families. They could be folded into the typical miRNA secondary structures. RepeatMasker algorithm demonstrated that ppt-miR2919e and pptmiR1533 had simple sequence repeats in their sequences. Target sites (49 genes) were identified for 7 out of 19 miRNAs. GO and KEGG analysis of targets indicated the involvement of some in important multiple biological and metabolic processes.

Conclusion:

The majority of the registered miRNAs in databases were predicted by computational approaches while many more have remained unknown. Due to the conserved nature of miRNAs in plant species from closely to distantly related, homology search-based approaches between plants species could lead to the identification of novel miRNAs in other plant species providing baseline information for further search about the biological functions and evolution of miRNAs.

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.

Identification of Cd-responsive RNA helicase genes and expression of a putative BnRH 24 mediated by miR158 in canola (Brassica napus)

RNA helicases play crucial roles in RNA splicing, transport, editing and degradation, protein translation initiation and siRNA-mediated gene silencing. However, knowledge about their functionality in rapeseed (Brassica napus) is rare. In the study, we identified and annotated 271 RNA helicase genes from B. napus using bioinformatics and high-throughput RNA-sequencing (RNA-seq). Three subfamilies DEAD-box, DEAH-box, or DExD/H-box have been identified. One hundred and ninety-five RNA helicases were confirmed by RNA-seq and 49 were identified to differentially respond to cadmium (Cd) stress (> 1.5 fold change, p < 0.05). As an example, we functionally specified BnaA04g26450D encoding a BnRH24 under Cd exposure. BnRH24 is a constitutive gene expressing throughout the life span. Using our previously generated degradome datasets, we found that BnRH24 can be cleaved by miR158, suggesting that BnRH24 is a target of miR158 in B. napus. The mature miR158 was induced, while BnRH24 was repressed in B. napus under Cd stress. The contrasting expression pattern of B. napus miR158 and BnRH24 under the normal and Cd would support the post-transcriptional regulation of BnRH24 by miR158. Ectopic expression of BnRH24 in Arabidopsis revealed that the transgenic lines showed more sensitivity to Cd toxicity by reducing root elongation, fresh mass production, chlorophyll accumulation and increasing oxidative products such as O₂^-., H₂O₂ and thiobarbituric acid reactive substances (TBARS), indicating that the controlling the level of BnRH24 by miR158 may be required for Cd tolerance in plants.

The evolution of posttranscriptional regulation

"DNA makes RNA makes protein." After transcription, mRNAs undergo a series of intertwining processes to be finally translated into functional proteins. The "posttranscriptional" regulation (PTR) provides cells an extended option to fine-tune their proteomes. To meet the demands of complex organism development and the appropriate response to environmental stimuli, every step in these processes needs to be finely regulated. Moreover, changes in these regulatory processes are important driving forces underlying the evolution of phenotypic differences across different species. The major PTR mechanisms discussed in this review include the regulation of splicing, polyadenylation, decay, and translation. For alternative splicing and polyadenylation, we mainly discuss their evolutionary dynamics and the genetic changes underlying the regulatory differences in cis-elements versus trans-factors. For mRNA decay and translation, which, together with transcription, determine the cellular RNA or protein abundance, we focus our discussion on how their divergence coordinates with transcriptional changes to shape the evolution of gene expression. Then to highlight the importance of PTR in the evolution of higher complexity, we focus on their roles in two major phenomena during eukaryotic evolution: the evolution of multicellularity and the division of labor between different cell types and tissues; and the emergence of diverse, often highly specialized individual phenotypes, especially those concerning behavior in eusocial insects. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution Translation > Translation Regulation RNA Processing > Splicing Regulation/Alternative Splicing.

Evolution of microRNA827 targeting in the plant kingdom

Unlike most ancient microRNAs, which conservatively target homologous genes across species, microRNA827 (miR827) targets two different types of SPX (SYG1/PHO81/XPR1)-domain-containing genes, NITROGEN LIMITATION ADAPTATION (NLA) and PHOSPHATE TRANSPORTER 5 (PHT5), in Arabidopsis thaliana and Oryza sativa to regulate phosphate (Pi) transport and storage, respectively. However, how miR827 shifted its target preference and its evolutionary history are unknown. Based on target prediction analysis, we found that in most angiosperms, miR827 conservatively targets PHT5 homologs, but in Brassicaceae and Cleomaceae it preferentially targets NLA homologs, and we provide evidence for the transition of target preference during Brassicales evolution. Intriguingly, we found a lineage-specific loss of the miR827-regulatory module in legumes. Analysis of miR827-mediated cleavage efficiency and the expression of PHT5 in A. thaliana indicated that accumulation of mutations in the target site and the exclusion of the target site by alternative transcriptional initiation eliminated PHT5 targeting by miR827. Here, we identified a transition of miR827 target preference during plant evolution and revealed the uniqueness of miR827-mediated regulation among conserved plant miRNAs. Despite the change in its target preference, upregulation of miR827 by Pi starvation and its role in regulating cellular Pi homeostasis were retained.

Genome-Wide Development of MicroRNA-Based SSR Markers in Medicago truncatula with Their Transferability Analysis and Utilization in Related Legume Species

Microsatellite (simple sequence repeats, SSRs) marker is one of the most widely used markers in marker-assisted breeding. As one type of functional markers, MicroRNA-based SSR (miRNA-SSR) markers have been exploited mainly in animals, but the development and characterization of miRNA-SSR markers in plants are still limited. In the present study, miRNA-SSR markers for Medicago truncatula (M. truncatula) were developed and their cross-species transferability in six leguminous species was evaluated. A total of 169 primer pairs were successfully designed from 130 M. truncatula miRNA genes, the majority of which were mononucleotide repeats (70.41%), followed by dinucleotide repeats (14.20%), compound repeats (11.24%) and trinucleotide repeats (4.14%). Functional classification of SSR-containing miRNA genes showed that all targets could be grouped into three Gene Ontology (GO) categories: 17 in biological process, 11 in molecular function, and 14 in cellular component. The miRNA-SSR markers showed high transferability in other six leguminous species, ranged from 74.56% to 90.53%. Furthermore, 25 Mt-miRNA-SSR markers were used to evaluate polymorphisms in 20 alfalfa accessions, and the polymorphism information content (PIC) values ranged from 0.39 to 0.89 with an average of 0.71, the allele number per marker varied from 3 to 18 with an average of 7.88, indicating a high level of informativeness. The present study is the first time developed and characterized of M. truncatula miRNA-SSRs and demonstrated their utility in transferability, these novel markers will be valuable for genetic diversity analysis, marker-assisted selection and genotyping in leguminous species.

Evolution of flower color pattern through selection on regulatory small RNAs

Small RNAs (sRNAs) regulate genes in plants and animals. Here, we show that population-wide differences in color patterns in snapdragon flowers are caused by an inverted duplication that generates sRNAs. The complexity and size of the transcripts indicate that the duplication represents an intermediate on the pathway to microRNA evolution. The sRNAs repress a pigment biosynthesis gene, creating a yellow highlight at the site of pollinator entry. The inverted duplication exhibits steep clines in allele frequency in a natural hybrid zone, showing that the allele is under selection. Thus, regulatory interactions of evolutionarily recent sRNAs can be acted upon by selection and contribute to the evolution of phenotypic diversity.

Global Analysis of Small RNA Dynamics during Seed Development of Picea glauca and Arabidopsis thaliana Populations Reveals Insights on their Evolutionary Trajectories

While DNA methylation carries genetic signals and is instrumental in the evolution of organismal complexity, small RNAs (sRNAs), ~18-24 ribonucleotide (nt) sequences, are crucial mediators of methylation as well as gene silencing. However, scant study deals with sRNA evolution via featuring their expression dynamics coupled with species of different evolutionary time. Here we report an atlas of sRNAs and microRNAs (miRNAs, single-stranded sRNAs) produced over time at seed-set of two major spermatophytes represented by populations of Picea glauca and Arabidopsis thaliana with different seed-set duration. We applied diverse profiling methods to examine sRNA and miRNA features, including size distribution, sequence conservation and reproduction-specific regulation, as well as to predict their putative targets. The top 27 most abundant miRNAs were highly overlapped between the two species (e.g., miR166,-319 and-396), but in P. glauca, they were less abundant and significantly less correlated with seed-set phases. The most abundant sRNAs in libraries were deeply conserved miRNAs in the plant kingdom for Arabidopsis but long sRNAs (24-nt) for P. glauca. We also found significant difference in normalized expression between populations for population-specific sRNAs but not for lineage-specific ones. Moreover, lineage-specific sRNAs were enriched in the 21-nt size class. This pattern is consistent in both species and alludes to a specific type of sRNAs (e.g., miRNA, tasiRNA) being selected for. In addition, we deemed 24 and 9 sRNAs in P. glauca and Arabidopsis, respectively, as sRNA candidates targeting known adaptive genes. Temperature had significant influence on selected gene and miRNA expression at seed development in both species. This study increases our integrated understanding of sRNA evolution and its potential link to genomic architecture (e.g., sRNA derivation from genome and sRNA-mediated genomic events) and organismal complexity (e.g., association between different sRNA expression and their functionality).

Comprehensive survey and evolutionary analysis of genome-wide miRNA genes from ten diploid Oryza species

BACKGROUND

MicroRNAs (miRNAs) are non-coding RNAs that play versatile roles in post-transcriptional gene regulation. Although much is known about their biogenesis, and gene regulation very little is known about their evolutionary relation among the closely related species.

RESULT

All the orthologous miRNA genes of Oryza sativa (japonica) from 10 different Oryza species were identified, and the evolutionary changes among these genes were analysed. Significant differences in the expansion of miRNA gene families were observed across the Oryza species. Analysis of the nucleotide substitution rates indicated that the mature sequences show the least substitution rates among the different regions of miRNA genes, and also show a very much less substitution rates as compared to that of all protein-coding genes across the Oryza species. Evolution of miRNA genes was also found to be contributed by transposons. A non-neutral selection was observed at 80 different miRNA loci across Oryza species which were estimated to have lost ~87% of the sequence diversity during the domestication. The phylogenetic analysis revealed that O. longistaminata diverged first among the AA-genomes, whereas O. brachyantha and O. punctata appeared as the eminent out-groups. The miR1861 family organised into nine distinct compact clusters in the studied Oryza species except O. brachyantha. Further, the expression analysis showed that 11 salt-responsive miRNAs were differentially regulated between O. coarctata and O. glaberrima.

CONCLUSION

Our study provides the evolutionary dynamics in the miRNA genes of 10 different Oryza species which will support more investigations about the structural and functional organization of miRNA genes of Oryza species.

Genome-wide identification and co-expression network analysis provide insights into the roles of auxin response factor gene family in chickpea

Auxin response factors (ARFs) are the transcription factors that regulate auxin responses in various aspects of plant growth and development. Although genome-wide analysis of ARF gene family has been done in some species, no information is available regarding ARF genes in chickpea. In this study, we identified 28 ARF genes (CaARF) in the chickpea genome. Phylogenetic analysis revealed that CaARFs can be divided into four different groups. Duplication analysis revealed that 50% of CaARF genes arose from duplication events. We analyzed expression pattern of CaARFs in various developmental stages. CaARF16.3, CaARF17.1 and CaARF17.2 showed highest expression at initial stages of flower bud development, while CaARF6.2 had higher expression at later stages of flower development. Further, CaARF4.2, CaARF9.2, CaARF16.2 and CaARF7.1 exhibited differential expression under different abiotic stress conditions, suggesting their role in abiotic stress responses. Co-expression network analysis among CaARF, CaIAA and CaGH3 genes enabled us to recognize components involved in the regulatory network associated with CaARFs. Further, we identified microRNAs that target CaARFs and TAS3 locus that trigger production of trans-acting siRNAs targeting CaARFs. The analyses presented here provide comprehensive information on ARF family members and will help in elucidating their exact function in chickpea.

Improvement in Fruit Quality by Overexpressing miR399a in Woodland Strawberry

Fruit quality is an important trait in strawberry and is determined by many factors. The soluble solid content in strawberry fruits is positively related to the phosphorus content. MicroRNA399 (miR399) is involved in the regulation of phosphate (Pi) homeostasis. However, the effect of miR399 on strawberry quality remains unknown. In this study, miR399a-overexpressing transgenic woodland strawberries (Fragaria vesca) were obtained via an Agrobacterium-mediated transformation. The phosphorus (P) content was 1.1-fold to 2.1-fold higher in the leaves and fruits of the miR399a-overexpressing plants than in the wild type (WT). However, the P content in the miR399a-overexpressing plants was decreased by 25% to 45% in the roots. The primary root length of the transgenic lines in both the high-Pi and low-Pi media was shorter than that of the WT. Interestingly, the transgenic lines in pots under Pi-sufficient conditions grew better than the WT, and the fruit quality, including the contents of fructose and glucose and soluble solid, was significantly higher in the transgenic lines than in the WT. The overexpression of miR399a in strawberry can be used to improve the parameters involved in fruit quality and provides information regarding breeding nutrient-improved strawberry.

miRNAs in the alga Chlamydomonas reinhardtii are not phylogenetically conserved and play a limited role in responses to nutrient deprivation

The unicellular alga Chlamydomonas reinhardtii contains many types of small RNAs (sRNAs) but the biological role(s) of bona fide microRNAs (miRNAs) remains unclear. To address their possible function(s) in responses to nutrient availability, we examined miRNA expression in cells cultured under different trophic conditions (mixotrophic in the presence of acetate or photoautotrophic in the presence or absence of nitrogen). We also reanalyzed miRNA expression data in Chlamydomonas subject to sulfur or phosphate deprivation. Several miRNAs were differentially expressed under the various trophic conditions. However, in transcriptome analyses, the majority of their predicted targets did not show expected changes in transcript abundance, suggesting that they are not subject to miRNA-mediated RNA degradation. Mutant strains, defective in sRNAs or in ARGONAUTE3 (a key component of sRNA-mediated gene silencing), did not display major phenotypic defects when grown under multiple nutritional regimes. Additionally, Chlamydomonas miRNAs were not conserved, even in algae of the closely related Volvocaceae family, and many showed features resembling those of recently evolved, species-specific miRNAs in the genus Arabidopsis. Our results suggest that, in C. reinhardtii, miRNAs might be subject to relatively fast evolution and have only a minor, largely modulatory role in gene regulation under diverse trophic states.

Evolution of the miR5200-FLOWERING LOCUS T flowering time regulon in the temperate grass subfamily Pooideae

Flowering time is a carefully regulated trait controlled primarily through the action of the central genetic regulator, FLOWERING LOCUS T (FT). Recently it was demonstrated that a microRNA, miR5200, targets the end of the second exon of FT under short-day photoperiods in the grass subfamily Pooideae, thus preventing FT transcripts from reaching threshold levels under non-inductive conditions. Pooideae are an interesting group in that they rapidly diversified from the tropics into the northern temperate region during a major global cooling event spanning the Eocene-Oligocene transition. We hypothesize that miR5200 photoperiod-sensitive regulation of Pooideae flowering time networks assisted their transition into northern seasonal environments. Here, we test predictions derived from this hypothesis that miR5200, originally found in bread wheat and later identified in Brachypodium distachyon, (1) was present in the genome of the Pooideae common ancestor, (2) is transcriptionally regulated by photoperiod, and (3) is negatively correlated with FT transcript abundance, indicative of miR5200 regulating FT. Our results demonstrate that miR5200 did evolve at or around the base of Pooideae, but only acquired photoperiod-regulated transcription within the Brachypodium lineage. Based on expression profiles and previous data, we posit that the progenitor of miR5200 was co-regulated with FT by an unknown mechanism.

The Emergence, Evolution, and Diversification of the miR390-TAS3-ARF Pathway in Land Plants

In plants, miR390 directs the production of tasiRNAs from TRANS-ACTING SIRNA3 (TAS3) transcripts to regulate AUXIN RESPONSIVE FACTOR (ARF) genes, critical for auxin signaling; these tasiRNAs are known as tasiARFs. To understand the evolution of this miR390-TAS3-ARF pathway, we characterized homologs of these three genes from thousands of plant species, from bryophytes to angiosperms. We found the lower-stem region of MIR390 genes, critical for accurate DICER-LIKE1 processing, is conserved in sequence in seed plants. We propose a model for the transition of functional tasiRNA sequences in TAS3 genes occurred at the emergence of vascular plants, in which the two miR390 target sites of TAS3 genes showed distinct pairing patterns. Based on the cleavability of miR390 target sites and the distance between target site and tasiARF, we inferred a potential bidirectional processing mechanism exists for some TAS3 genes. We also demonstrated a tight mutual selection between tasiARF and its target genes and that ARGONAUTE7, the partner of miR390, was specified later than other factors in the pathway. All these data illuminate the evolutionary path of the miR390-TAS3-ARF pathway in land plants and demonstrate the significant variation that occurs in this functionally important and archetypal regulatory circuit.

Long noncoding miRNA gene represses wheat β-diketone waxes

The cuticle of terrestrial plants functions as a protective barrier against many biotic and abiotic stresses. In wheat and other Triticeae, β-diketone waxes are major components of the epicuticular layer leading to the bluish-white glaucous trait in reproductive-age plants. Glaucousness in durum wheat is controlled by a metabolic gene cluster at the WAX1 (W1) locus and a dominant suppressor INHIBITOR of WAX1 (Iw1) on chromosome 2B. The wheat D subgenome from progenitor Aegilops tauschii contains W2 and Iw2 paralogs on chromosome 2D. Here we identify the Iw1 gene from durum wheat and demonstrate the unique regulatory mechanism by which Iw1 acts to suppress a carboxylesterase-like protein gene, W1-COE, within the W1 multigene locus. Iw1 is a long noncoding RNA (lncRNA) containing an inverted repeat (IR) with >80% identity to W1-COE The Iw1 transcript forms a miRNA precursor-like long hairpin producing a 21-nt predominant miRNA, miRW1, and smaller numbers of related sRNAs associated with the nonglaucous phenotype. When Iw1 was introduced into glaucous bread wheat, miRW1 accumulated, W1-COE and its paralog W2-COE were down-regulated, and the phenotype was nonglaucous and β-diketone-depleted. The IR region of Iw1 has >94% identity to an IR region on chromosome 2 in Ae. tauschii that also produces miRW1 and lies within the marker-based location of Iw2 We propose the Iw loci arose from an inverted duplication of W1-COE and/or W2-COE in ancestral wheat to form evolutionarily young miRNA genes that act to repress the glaucous trait.

Genome-wide analysis of microRNA targeting impacted by SNPs in cucumber genome

Background

microRNAs (miRNAs) are endogenous small RNAs that play important regulatory functions in plant development. Genetic variations in miRNAs sequences or their target-binding sites (microRNA-target interaction sites) can alter miRNA targets in animal and human. Whether these single nucleotide polymorphisms (SNPs) in plant are functional have not yet been determined.

Results

In this study, we constructed leaf, root, and stem-derived small libraries of cucumber (Cucumis sativus) line 9930 (cultivated China-group cucumber) and C. sativus var. hardwickii (wild India group cucumber). A total of 22 conserved miRNA families, nine less-conserved miRNA families, and 49 cucumber-specific miRNAs were identified in both line 9930 and hardwickii. We employed cucumber resequencing data to perform a genome-wide scan for SNPs in cucumber miRNA-target interaction sites, including miRNA mature sequences and miRNA-target binding sites. As a result, we identified a total of 19 SNPs in mature miRNA sequences and 113 SNPs in miRNA-target binding sites with the potential to affect miRNA-target interactions. Furthermore, we experimentally confirmed that these SNPs produced 14 9930-unique targets mRNAs and 15 hardwickii-unique targets mRNA for cucumber miRNAs. This is the first experimental validation of SNPs in miRNA-target interaction sites affecting miRNA-target binding in plants.

Conclusions

Our results indicate that SNPs can alter miRNA function and produce unique miRNA targets in cultivated and wild cucumbers. Therefore, miRNA-related SNPs may have played important in events that led to the agronomic differences between domestic and wild cucumber.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3665-y) contains supplementary material, which is available to authorized users.

The evolutionary origin of plant and animal microRNAs

microRNAs (miRNAs) are a unique class of short endogenous RNAs that became known in the last few decades as major players in gene regulation at the post-transcriptional level. Their regulatory roles make miRNAs crucial for normal development and physiology in several distinct groups of eukaryotes including plants and animals. The common notion in the field is that miRNAs have evolved independently in those distinct lineages, but recent evidence from non-bilaterian metazoans, plants, as well as various algae raise the possibility that already the last common ancestor of these lineages might have employed a miRNA pathway for post-transcriptional regulation. In this review we present the commonalities and differences of the miRNA pathways in various eukaryotes and discuss the contrasting scenarios of their possible evolutionary origin and their proposed link to organismal complexity and multicellularity.