Dual coding of siRNAs and miRNAs by plant transposable elements

Non-conserved microRNAs and their roles in plants: the case for legumes

Several classes of small RNAs function to regulate stress and development pathways in all kingdoms of life. In animals and plants, microRNAs have been widely studied as important regulators of gene expression. However, non-conserved microRNAs have proven more difficult to study, raising questions as to their functionality. Using the legume family of plants as reference, we discuss this concept and provide examples where miRNAs functions have been described, highlighting their potential role in regulating important processes in these plants, such as stress responses and communication with other organisms, including bacteria and fungi. These examples suggest that non-conserved miRNAs are likely to contribute to more gene regulation circuits than currently appreciated, and in a wider range of plant species.

Scenarios for the emergence of new microRNA genes in the plant Arabidopsis halleri

MicroRNAs (miRNAs) are central players in the regulation of gene expression in eukaryotes. The repertoires of miRNA genes vary drastically even among closely related species, indicating that they are evolutionarily labile. However, the processes by which they originate over the course of evolution and the nature of their progenitors across the genome remain poorly understood. Here, we analyzed miRNA genes in Arabidopsis halleri, a plant species where we recently documented a large number of species-specific miRNA genes, likely to represent recent events of emergence. Analysis of sequence homology across the genome indicates that a diversity of sources contributes to the emergence of new miRNA genes, including inverted duplications from protein-coding genes, rearrangements of transposable element (TE) sequences, and duplications of preexisting miRNA genes. Our observations indicate that the origin from protein-coding genes was less common than was previously considered. In contrast, we estimate that almost half of the new miRNA genes likely emerged from TEs. Miniature inverted-repeat TEs (MITEs) seem to be particularly important contributors to new miRNA genes, with the Harbinger and Mariner TE superfamilies representing disproportionate sources for their emergence. We further analyzed the recent expansion of a miRNA family derived from MuDR elements and the duplication of miRNA genes formed by two hAT transposons. Overall, our results illustrate the rapid pace at which new regulatory elements can arise from the modification of preexisting sequences in a genome and highlight the central role of certain categories of TEs in this process.

Evolution of Plant Conserved microRNAs After Whole-Genome Duplications

MicroRNAs (miRNAs) are a specialized class of small silencing RNAs that regulate gene expression in numerous biological processes in eukaryotes. While the evolutionary dynamics of protein-coding genes after plant whole-genome duplications (WGDs) has been extensively studied, the patterns of evolution for conserved MIRNAs (miRNA genes) post-WGDs are less understood. In this study, we systematically investigated miRNAs and their targets in 6 plant species with varying WGD histories. Our findings reveal that WGDs significantly contribute to the expansion of conserved miRNA families. Notably, through homologous analyses of conserved miRNA families, we discovered that beyond the loci derived from WGDs and other duplication events, some conserved miRNA families have independently gained new loci and/or lost syntenic loci in specific lineage or species through evolution. Additionally, our analyses of sequence divergence in conserved miRNAs showed that the mature sequences of miRNA duplicates gradually diverge following WGDs, with this sequence divergence being correlated with that of their adjacent protein-coding genes after recent WGDs. Furthermore, expression and functional divergence analyses of duplicated targets in different miRNA-target interaction scenarios suggest that conserved miRNAs may play crucial roles in regulating the expression of duplicated genes and related regulatory networks following WGDs. In summary, our analyses reveal universal evolutionary patterns of plant conserved miRNAs following WGDs and provide evidence that some miRNA copies in conserved families originated independently during evolution.

Deciphering the impact of microRNAs in plant biology: a review of computational insights and experimental validation

Exploring the complex world of microRNA (miRNA) biogenesis and functions in plants is essential for understanding their diverse regulatory mechanisms. This review highlights the processes involved in miRNA biogenesis and their crucial roles in growth and development of plant, stress responses, and nutrient homeostasis. miRNAs play a central role in various developmental processes, including the transition from the juvenile to adult stage, the growth of shoot apical meristem, leaf and floral morphogenesis, and the determination of flowering time. By presenting the current state of research, we focus on the vital role of computational tools and databases in deciphering the regulatory networks controlled by miRNAs, which helps us navigate the intricate world of plant biology. Furthermore, it stresses the importance of experimental validation techniques in confirming computational predictions, ensuring that miRNA research is reliable and robust. As the field continues to grow, this review emphasizes the urgent need for integrated approaches, to deepen our knowledge of plant miRNA biology and its implications. These insights will pave the way for advancements in crop improvement, stress resilience, and biotechnological innovations.

An intron-split microRNA mediates cleavage of the mRNA encoded by low phosphate root in Solanaceae

MAIN CONCLUSION

A microRNA with a non-canonical precursor structure harbours an intron in between its miRNA-5p and miRNA-3p relevant for its biogenesis, is conserved across Solanaceae, and targets the mRNA of low phosphate root. Hundreds of miRNAs have been identified in plants and great advances have been accomplished in the understanding of plant miRNA biogenesis, mechanisms and functions. Still, many miRNAs, particularly those with less conventional features, remain to be discovered. Likewise, additional layers of regulation from miRNA generation to action and turnover are still being revealed. The current study describes a microRNA not previously identified given its unusual intron-split stem-loop structure, that has been previously observed only within the monocot-specific miRNA444 family. It shows its conservation across a branch of Solanales including agriculturally relevant Solanaceae family, where its transcripts had already been predicted in several species within sequence databases. The miRNA is absent in Arabidopsis thaliana but present in Solanum lycopersicum, Nicotiana benthamiana, Petunia axillaris, and Ipomoea nil. It proves that at least two different pri-miRNA variants are produced from this miRNA gene, one spliced and the other one retaining the intron. It demonstrates the dual function of its intron in the miRNA biogenesis. On the one hand, its presence in the pri-miRNA positively influences mature miRNA accumulation, but on the other hand, it needs to be removed from the pri-miRNA for efficient mature miRNA production. Finally, it sets low phosphate root as one of its targets, a protein known to be involved in root growth regulation under phosphate starvation in other plant species.

Advantage looping: Gene regulatory circuits between microRNAs and their target transcription factors in plants

The 20 to 24 nucleotide microRNAs (miRNAs) and their target transcription factors (TF) have emerged as key regulators of diverse processes in plants, including organ development and environmental resilience. In several instances, the mature miRNAs degrade the TF-encoding transcripts, while their protein products in turn bind to the promoters of the respective miRNA-encoding genes and regulate their expression, thus forming feedback loops (FBLs) or feedforward loops (FFLs). Computational analysis suggested that such miRNA-TF loops are recurrent motifs in gene regulatory networks (GRNs) in plants as well as animals. In recent years, modeling and experimental studies have suggested that plant miRNA-TF loops in GRNs play critical roles in driving organ development and abiotic stress responses. Here, we discuss the miRNA-TF FBLs and FFLs that have been identified and studied in plants over the past decade. We then provide some insights into the possible roles of such motifs within GRNs. Lastly, we provide perspectives on future directions for dissecting the functions of miRNA-centric GRNs in plants.

ARGONAUTE1-binding Tudor domain proteins function in small interfering RNA production for RNA-directed DNA methylation

Transposable elements (TEs) contribute to plant evolution, development, and adaptation to environmental changes, but the regulatory mechanisms are largely unknown. RNA-directed DNA methylation (RdDM) is 1 TE regulatory mechanism in plants. Here, we identified that novel ARGONAUTE 1 (AGO1)-binding Tudor domain proteins Precocious dissociation of sisters C/E (PDS5C/E) are involved in 24-nt siRNA production to establish RdDM on TEs in Arabidopsis thaliana. PDS5 family proteins are subunits of the eukaryote-conserved cohesin complex. However, the double mutant lacking angiosperm-specific subfamily PDS5C and PDS5E (pds5c/e) exhibited different developmental phenotypes and transcriptome compared with those of the double mutant lacking eukaryote-conserved subfamily PDS5A and PDS5B (pds5a/b), suggesting that the angiosperm-specific PDS5C/E subfamily has a unique function in angiosperm plants. Proteome and imaging analyses revealed that PDS5C/E interact with AGO1. The pds5c/e double mutant had defects in 24-nt siRNA accumulation and CHH DNA methylation on TEs. In addition, some lncRNAs that accumulated in the pds5c/e mutant were targeted by AGO1-loading 21-nt miRNAs and 21-nt siRNAs. These results indicate that PDS5C/E and AGO1 participate in 24-nt siRNA production for RdDM in the cytoplasm. These findings indicate that angiosperm plants evolved a new regulator, the PDS5C/E subfamily, to control the increase in TEs during angiosperm evolution.

A transposable element-derived siRNAs involve DNA hypermethylation at the promoter of OsGSTZ4 for cadmium tolerance in rice

Environmental contaminants such as cadmium (Cd) pose high risks to crop production and human health. The genetic basis for regulation of Cd stress-responsive genes for plant adaptation to adverse environments remains poorly understood. In this study, we characterized a rice Zeta family glutathione-S-transferase (OsGSTZ4) gene for Cd detoxification. Heterologous expression of OsGSTZ4 in a yeast (Saccharomyces cerevisiae) conferred cellular Cd tolerance. Transgenic rice overexpressing OsGSTZ4 improved plant growth, attenuated Cd-induced toxicity, and accumulated more Cd in roots. OsGSTZ4 transcription was rapidly induced 3 h after Cd exposure and then declined to the basal level. This was followed by (days after Cd treatment) by CHH hypermethylation (by 41.2 %) at a MITE (Miniature Inverted-repeat Transposable Element) transposable element (TE) inserted in the 5'-untranscribed region (UTR) (-1,722 ∼ -1,392 bp) of OsGSTZ4. Meanwhile, three 24-nt siRNAs derived from the TE (-1,722 ∼ -1,471 bp) were detected and was also rapidly enriched under Cd stress. To validate the possibility that Cd-induced change in OsGSTZ4 expression correlates with the siRNAs-involved CHH methylation through an RdDM (RNA-directed DNA methylation) pathway, genetic analyses were performed. We found that the CHH methylation at the promoter and transcript level of OsGSTZ4 were compromised in the osdrm2 (loss of function for CHH methylation) and osrdr2i (defective in RNA-dependent RNA polymerase 2) but did not change in other types of methyltransferases such as osmet1, ossdg714 or osros1. Promoter deletion analyses confirmed that the siRNA target sequences were essential for the proper expression of OsGSTZ4. Our studies reveal an unusual feedback mechanism by which the Cd-induced rapid OsGSTZ4 expression for Cd tolerance would interplay with the late CHH hypermethylation to silence the TE through the 24-nt siRNAs- and Osdrm2-mediated RdDM pathway, and help understand the diversity of gene regulation via an epigenetic mechanism for rice adaptation to the environmental stress.

Generation of de novo miRNAs from template switching during DNA replication

The mechanisms generating novel genes and genetic information are poorly known, even for microRNA (miRNA) genes with an extremely constrained design. All miRNA primary transcripts need to fold into a stem-loop structure to yield short gene products ([Formula: see text]22 nt) that bind and repress their mRNA targets. While a substantial number of miRNA genes are ancient and highly conserved, short secondary structures coding for entirely novel miRNA genes have been shown to emerge in a lineage-specific manner. Template switching is a DNA-replication-related mutation mechanism that can introduce complex changes and generate perfect base pairing for entire hairpin structures in a single event. Here, we show that the template-switching mutations (TSMs) have participated in the emergence of over 6,000 suitable hairpin structures in the primate lineage to yield at least 18 new human miRNA genes, that is 26% of the miRNAs inferred to have arisen since the origin of primates. While the mechanism appears random, the TSM-generated miRNAs are enriched in introns where they can be expressed with their host genes. The high frequency of TSM events provides raw material for evolution. Being orders of magnitude faster than other mechanisms proposed for de novo creation of genes, TSM-generated miRNAs enable near-instant rewiring of genetic information and rapid adaptation to changing environments.

Short Interrupted Repeat Cassette (SIRC)-Novel Type of Repetitive DNA Element Found in Arabidopsis thaliana

Short interrupted repeat cassette (SIRC)-a novel DNA element found throughout the A. thaliana nuclear genome. SIRCs are represented by short direct repeats interrupted by diverse DNA sequences. The maxima of SIRC's distribution are located within pericentromeric regions. We suggest that originally SIRC was a special case of the complex internal structure of the miniature inverted repeat transposable element (MITE), and further MITE amplification, transposition, and loss of terminal inverted repeats gave rise to SIRC as an independent DNA element. SIRC sites were significantly enriched with several histone modifications associated with constitutive heterochromatin and mobile genetic elements. The majority of DNA-binding proteins, strongly associated with SIRC, are related to histone modifications for transcription repression. A part of SIRC was found to overlap highly inducible protein-coding genes, suggesting a possible regulatory role for these elements, yet their definitive functions need further investigation.

Genome-wide identification of nitrate-responsive microRNAs by small RNA sequencing in the rice restorer cultivar Nanhui 511

Rice productivity relies heavily on nitrogen fertilization, and improving nitrogen use efficiency (NUE) is important for hybrid rice breeding. Reducing nitrogen inputs is the key to achieving sustainable rice production and reducing environmental problems. Here, we analyzed the genome-wide transcriptomic changes in microRNAs (miRNAs) in the indica rice restorer cultivar Nanhui 511 (NH511) under high (HN) and low nitrogen (LN) conditions. The results showed that NH511 is sensitive to nitrogen supplies and HN conditions promoted the growth its lateral roots at the seedling stage. Furthermore, we identified 483 known miRNAs and 128 novel miRNAs by small RNA sequencing in response to nitrogen in NH511. We also detected 100 differentially expressed genes (DEGs), including 75 upregulated and 25 downregulated DEGs, under HN conditions. Among these DEGs, 43 miRNAs that exhibited a 2-fold change in their expression were identified in response to HN conditions, including 28 upregulated and 15 downregulated genes. Additionally, some differentially expressed miRNAs were further validated by qPCR analysis, which showed that miR443, miR1861b, and miR166k-3p were upregulated, whereas miR395v and miR444b.1 were downregulated under HN conditions. Moreover, the degradomes of possible target genes for miR166k-3p and miR444b.1 and expression variations were analyzed by qPCR at different time points under HN conditions. Our findings revealed comprehensive expression profiles of miRNAs responsive to HN treatments in an indica rice restorer cultivar, which advances our understanding of the regulation of nitrogen signaling mediated by miRNAs and provides novel data for high-NUE hybrid rice cultivation.

Perspective for Studying the Relationship of miRNAs with Transposable Elements

Transposable elements are important sources of miRNA, long non-coding RNAs genes, and their targets in the composition of protein-coding genes in plants and animals. Therefore, the detection of expression levels of specific non-coding RNAs in various tissues and cells in normal and pathological conditions may indicate a programmed pattern of transposable elements' activation. This reflects the species-specific composition and distribution of transposable elements in genomes, which underlie gene regulation in every cell division, including during aging. TEs' expression is also regulated by epigenetic factors (DNA methylation, histone modifications), SIRT6, cytidine deaminases APOBEC3, APOBEC1, and other catalytic proteins, such as ERCC, TREX1, RB1, HELLS, and MEGP2. In evolution, protein-coding genes and their regulatory elements are derived from transposons. As part of non-coding regions and introns of genes, they are sensors for transcriptional and post-transcriptional control of expression, using miRNAs and long non-coding RNAs, that arose from transposable elements in evolution. Methods (Orbld, ncRNAclassifier) and databases have been created for determining the occurrence of miRNAs from transposable elements in plants (PlanTE-MIR DB, PlaNC-TE), which can be used to design epigenetic gene networks in ontogenesis. Based on the data accumulated in the scientific literature, the presence of 467 transposon-derived miRNA genes in the human genome has been reliably established. It was proposed to create an updated and controlled online bioinformatics database of miRNAs derived from transposable elements in healthy individuals, as well as expression changes of these miRNAs during aging and various diseases, such as cancer and difficult-to-treat diseases. The use of the information obtained can open new horizons in the management of tissue and organ differentiation to aging slow down. In addition, the created database could become the basis for clarifying the mechanisms of pathogenesis of various diseases (imbalance in the activity of transposable elements, reflected in changes in the expression of miRNAs) and designing their targeted therapy using specific miRNAs as targets. This article provides examples of the detection of transposable elements-derived miRNAs involved in the development of specific malignant neoplasms, aging, and idiopathic pulmonary fibrosis.

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.

Repetitive DNA is Functional and Encodes Parts of the Non-Coding RNA Repertoire

This is a commentary on the article by Eviatar Nevo and Kexin Li entitled "Sympatric Speciation in Mole Rats and Wild Barley and Their Genome Repeatome Evolution: A Commentary", published recently in Advanced Genetics.

A review of strategies used to identify transposition events in plant genomes

Transposable elements (TEs) were initially considered redundant and dubbed 'junk DNA'. However, more recently they were recognized as an essential element of genome plasticity. In nature, they frequently become active upon exposition of the host to stress conditions. Even though most transposition events are neutral or even deleterious, occasionally they may happen to be beneficial, resulting in genetic novelty providing better fitness to the host. Hence, TE mobilization may promote adaptability and, in the long run, act as a significant evolutionary force. There are many examples of TE insertions resulting in increased tolerance to stresses or in novel features of crops which are appealing to the consumer. Possibly, TE-driven de novo variability could be utilized for crop improvement. However, in order to systematically study the mechanisms of TE/host interactions, it is necessary to have suitable tools to globally monitor any ongoing TE mobilization. With the development of novel potent technologies, new high-throughput strategies for studying TE dynamics are emerging. Here, we present currently available methods applied to monitor the activity of TEs in plants. We divide them on the basis of their operational principles, the position of target molecules in the process of transposition and their ability to capture real cases of actively transposing elements. Their possible theoretical and practical drawbacks are also discussed. Finally, conceivable strategies and combinations of methods resulting in an improved performance are proposed.

Evolution of lmiRNAs and their targets from MITEs for rice adaptation

Twenty-four nucleotide long microRNAs (lmiRNAs) direct DNA methylation at target genes and regulate their transcription. The evolutionary origin of lmiRNAs and the range of lmiRNA-mediated regulation remain obscure. Here, we reannotated lmiRNAs and their targets in rice by applying stringent criteria. We found that the majority of lmiRNAs are derived from Miniature Inverted-repeat Transposable Elements (MITEs) and most sites targeted by MITE-derived lmiRNAs reside within MITEs, suggesting co-evolution of lmiRNAs and their targets through MITE amplification. lmiRNAs undergo dynamically changes under stress conditions and the genes targeted by lmiRNAs show an enrichment for stress-responsive genes, suggesting that lmiRNAs are widely involved in plant responses to stresses. We constructed the evolutionary histories of lmiRNAs and their targets. Nearly half of lmiRNAs emerged before or when the AA genome was diverged, while the emergence of lmiRNA targets coincided with or followed the emergence of lmiRNAs. Furthermore, we found that the sequences of a lmiRNA target site underwent variations, coincident with the divergence of rice accessions and the distribution of rice accessions in different geographical locations and climatic conditions. Our findings highlight MITEs as an important origin of lmiRNAs and suggest that the evolution of lmiRNA-target regulatory modules may contribute to rice adaptation to environmental changes.

The Epigenetic Control of the Transposable Element Life Cycle in Plant Genomes and Beyond

Within the life cycle of a living organism, another life cycle exists for the selfish genome inhabitants, which are called transposable elements (TEs). These mobile sequences invade, duplicate, amplify, and diversify within a genome, increasing the genome's size and generating new mutations. Cells act to defend their genome, but rather than permanently destroying TEs, they use chromatin-level repression and epigenetic inheritance to silence TE activity. This level of silencing is ephemeral and reversible, leading to a dynamic equilibrium between TE suppression and reactivation within a host genome. The coexistence of the TE and host genome can also lead to the domestication of the TE to serve in host genome evolution and function. In this review, we describe the life cycle of a TE, with emphasis on how epigenetic regulation is harnessed to control TEs for host genome stability and innovation.

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

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

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.

Transposable Elements: Major Players in Shaping Genomic and Evolutionary Patterns

Transposable elements (TEs) are ubiquitous genetic elements, able to jump from one location of the genome to another, in all organisms. For this reason, on the one hand, TEs can induce deleterious mutations, causing dysfunction, disease and even lethality in individuals. On the other hand, TEs can increase genetic variability, making populations better equipped to respond adaptively to environmental change. To counteract the deleterious effects of TEs, organisms have evolved strategies to avoid their activation. However, their mobilization does occur. Usually, TEs are maintained silent through several mechanisms, but they can be reactivated during certain developmental windows. Moreover, TEs can become de-repressed because of drastic changes in the external environment. Here, we describe the ‘double life’ of TEs, being both ‘parasites’ and ‘symbionts’ of the genome. We also argue that the transposition of TEs contributes to two important evolutionary processes: the temporal dynamic of evolution and the induction of genetic variability. Finally, we discuss how the interplay between two TE-dependent phenomena, insertional mutagenesis and epigenetic plasticity, plays a role in the process of evolution.

Genome-wide identification of MITE-derived microRNAs and their targets in bread wheat

Background

Plant miRNAs are a class of small non-coding RNAs that can repress gene expression at the post-transcriptional level by targeting RNA degradation or promoting translational repression. There is increasing evidence that some miRNAs can derive from a group of non-autonomous class II transposable elements called Miniature Inverted-repeat Transposable Elements (MITEs).

Results

We used public small RNA and degradome libraries from Triticum aestivum to screen for microRNAs production and predict their cleavage target sites. In parallel, we also created a comprehensive wheat MITE database by identifying novel elements and compiling known ones. When comparing both data sets, we found high homology between MITEs and 14% of all the miRNAs production sites detected. Furthermore, we show that MITE-derived miRNAs have preference for targeting degradation sites with MITE insertions in the 3’ UTR regions of the transcripts.

Conclusions

Our results revealed that MITE-derived miRNAs can underlay the origin of some miRNAs and potentially shape a regulatory gene network. Since MITEs are found in millions of insertions in the wheat genome and are closely linked to genic regions, this kind of regulatory network could have a significant impact on the post-transcriptional control of gene expression.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-022-08364-4).

Nucleic acids movement and its relation to genome dynamics of repetitive DNA: Is cellular and intercellular movement of DNA and RNA molecules related to the evolutionary dynamic genome components?: Is cellular and intercellular movement of DNA and RNA molecules related to the evolutionary dynamic genome components?

There is growing evidence of evolutionary genome plasticity. The evolution of repetitive DNA elements, the major components of most eukaryotic genomes, involves the amplification of various classes of mobile genetic elements, the expansion of satellite DNA, the transfer of fragments or entire organellar genomes and may have connections with viruses. In addition to various repetitive DNA elements, a plethora of large and small RNAs migrate within and between cells during individual development as well as during evolution and contribute to changes of genome structure and function. Such migration of DNA and RNA molecules often results in horizontal gene transfer, thus shaping the whole genomic network of interconnected species. Here, we propose that a high evolutionary dynamism of repetitive genome components is often related to the migration/movement of DNA or RNA molecules. We speculate that the cytoplasm is probably an ideal compartment for such evolutionary experiments.

CRISPR/Cas9 mutants of tomato MICRORNA164 genes uncover their functional specialization in development

Plant MICRORNA164 (miR164) plays diverse regulatory functions by post-transcriptional repression of certain NAM/ATAF/CUC-domain transcription factors. However, the involvement of miR164 in fleshy fruit development and ripening remains poorly understood. Here, de novo prediction of tomato (Solanum lycopersicum) MIR164 genes identified four genes (SlMIR164a-d), of which SlMIR164d has an atypically long pre-miRNA. The roles of the fruit expressed SlMIR164a, b, and d were studied by analysis of their Clustered Regularly Interspaced Short Palindromic Repeats mutants. The slmir164bCR mutant plants exhibited shoot and flower abnormalities characteristic of ectopic boundary specification, whereas the shoot and flower development of slmir164aCR and slmir164dCR mutants were indistinguishable from wild-type. Strikingly, the knockout of SlMIR164a practically eliminated sly-miR164 from the developing and ripening fruit pericarp. The sly-miR164-deficient slmir164aCR fruits were smaller than the wild-type, due to reduced pericarp cell division and expansion, and displayed intense red color and matte, instead of glossy appearance, upon ripening. We found that the fruit skin phenotypes were associated with morphologically abnormal outer epidermis and thicker cuticle. Quantitation of sly-miR164 target transcripts in slmir164aCR ripening fruits demonstrated the upregulation of SlNAM3 and SlNAM2. Specific expression of their miR164-resistant versions in the pericarp resulted in the formation of extremely small fruits with abnormal epidermis, highlighting the importance of their negative regulation by sly-miR164a. Taken together, our results demonstrate that SlMIR164a and SlMIR164b play specialized roles in development: SlMIR164b is required for shoot and flower boundary specification, and SlMIR164a is required for fruit growth including the expansion of its outer epidermis, which determines the properties of the fruit skin.

A novel Transposable element-derived microRNA participates in plant immunity to rice blast disease

MicroRNAs (miRNAs) are small non-coding RNAs that direct post-transcriptional gene silencing in plant development and stress responses through cleavage or translational repression of target mRNAs. Here, we report the identification and functional characterization of a new member of the miR812 family in rice (named as miR812w) involved in disease resistance. miR812w is present in cultivated Oryza species, both japonica and indica subspecies, and wild rice species within the Oryza genus, but not in dicotyledonous species. miR812w is a 24nt-long that requires DCL3 for its biogenesis and is loaded into AGO4 proteins. Whereas overexpression of miR812w increased resistance to infection by the rice blast fungus Magnaporthe oryzae, CRISPR/Cas9-mediated MIR812w editing enhances disease susceptibility, supporting that miR812w plays a role in blast resistance. We show that miR812w derives from the Stowaway type of rice MITEs (Miniature Inverted-Repeat Transposable Elements). Moreover, miR812w directs DNA methylation in trans at target genes that have integrated a Stowaway MITE copy into their 3' or 5' untranslated region (ACO3, CIPK10, LRR genes), as well as in cis at the MIR812w locus. The target genes of miR812 were found to be hypo-methylated around the miR812 recognition site, their expression being up-regulated in transgene-free CRISPR/Cas9-edited miR812 plants. These findings further support that, in addition to post-transcriptional regulation of gene expression, miRNAs can exert their regulatory function at the transcriptional level. This relationship between miR812w and Stowaway MITEs integrated into multiple coding genes might eventually create a network for miR812w-mediated regulation of gene expression with implications in rice immunity.

Impact of Repetitive DNA Elements on Snake Genome Biology and Evolution

The distinctive biology and unique evolutionary features of snakes make them fascinating model systems to elucidate how genomes evolve and how variation at the genomic level is interlinked with phenotypic-level evolution. Similar to other eukaryotic genomes, large proportions of snake genomes contain repetitive DNA, including transposable elements (TEs) and satellite repeats. The importance of repetitive DNA and its structural and functional role in the snake genome, remain unclear. This review highlights the major types of repeats and their proportions in snake genomes, reflecting the high diversity and composition of snake repeats. We present snakes as an emerging and important model system for the study of repetitive DNA under the impact of sex and microchromosome evolution. We assemble evidence to show that certain repetitive elements in snakes are transcriptionally active and demonstrate highly dynamic lineage-specific patterns as repeat sequences. We hypothesize that particular TEs can trigger different genomic mechanisms that might contribute to driving adaptive evolution in snakes. Finally, we review emerging approaches that may be used to study the expression of repetitive elements in complex genomes, such as snakes. The specific aspects presented here will stimulate further discussion on the role of genomic repeats in shaping snake evolution.

Genome-Wide Analysis of Coding and Non-coding RNA Reveals a Conserved miR164-NAC-mRNA Regulatory Pathway for Disease Defense in Populus

MicroRNAs (miRNAs) contribute to plant defense responses by increasing the overall genetic diversity; however, their origins and functional importance in plant defense remain unclear. Here, we employed Illumina sequencing technology to assess how miRNA and messenger RNA (mRNA) populations vary in the Chinese white poplar (Populus tomentosa) during a leaf black spot fungus (Marssonina brunnea) infection. We sampled RNAs from infective leaves at conidia germinated stage [12 h post-inoculation (hpi)], infective vesicles stage (24 hpi), and intercellular infective hyphae stage (48 hpi), three essential stages associated with plant colonization and biotrophic growth in M. brunnea fungi. In total, 8,938 conserved miRNA-target gene pairs and 3,901 Populus-specific miRNA-target gene pairs were detected. The result showed that Populus-specific miRNAs (66%) were more involved in the regulation of the disease resistance genes. By contrast, conserved miRNAs (>80%) target more whole-genome duplication (WGD)-derived transcription factors (TFs). Among the 1,023 WGD-derived TF pairs, 44.9% TF pairs had only one paralog being targeted by a miRNA that could be due to either gain or loss of a miRNA binding site after the WGD. A conserved hierarchical regulatory network combining promoter analyses and hierarchical clustering approach uncovered a miR164–NAM, ATAF, and CUC (NAC) transcription factor–mRNA regulatory module that has potential in Marssonina defense responses. Furthermore, analyses of the locations of miRNA precursor sequences reveal that pseudogenes and transposon contributed a certain proportion (∼30%) of the miRNA origin. Together, these observations provide evolutionary insights into the origin and potential roles of miRNAs in plant defense and functional innovation.

Noncoding RNAs in Medicinal Plants and their Regulatory Roles in Bioactive Compound Production

BACKGROUND

Noncoding RNAs (ncRNAs), such as microRNAs (miRNAs), small interfering RNAs (siRNAs) and long noncoding RNAs (lncRNAs), play significant regulatory roles in plant development and secondary metabolism and are involved in plant response to biotic and abiotic stresses. They have been intensively studied in model systems and crops for approximately two decades and massive amount of information have been obtained. However, for medicinal plants, ncRNAs, particularly their regulatory roles in bioactive compound biosynthesis, are just emerging as a hot research field.

OBJECTIVE

This review aims to summarize current knowledge on herbal ncRNAs and their regulatory roles in bioactive compound production.

RESULTS

So far, scientists have identified thousands of miRNA candidates from over 50 medicinal plant species and 11794 lncRNAs from Salvia miltiorrhiza, Panax ginseng, and Digitalis purpurea. Among them, more than 30 miRNAs and five lncRNAs have been predicted to regulate bioactive compound production.

CONCLUSION

The regulation may achieve through various regulatory modules and pathways, such as the miR397-LAC module, the miR12112-PPO module, the miR156-SPL module, the miR828-MYB module, the miR858-MYB module, and other siRNA and lncRNA regulatory pathways. Further functional analysis of herbal ncRNAs will provide useful information for quality and quantity improvement of medicinal plants.

When junk DNA turns functional: transposon-derived non-coding RNAs in plants

Transposable elements (TEs) are major contributors to genome complexity in eukaryotes. TE mobilization may cause genome instability, although it can also drive genome diversity throughout evolution. TE transposition may influence the transcriptional activity of neighboring genes by modulating the epigenomic profile of the region or by altering the relative position of regulatory elements. Notably, TEs have emerged in the last few years as an important source of functional long and small non-coding RNAs. A plethora of small RNAs derived from TEs have been linked to the trans regulation of gene activity at the transcriptional and post-transcriptional levels. Furthermore, TE-derived long non-coding RNAs have been shown to modulate gene expression by interacting with protein partners, sequestering active small RNAs, and forming duplexes with DNA or other RNA molecules. In this review, we summarize our current knowledge of the functional and mechanistic paradigms of TE-derived long and small non-coding RNAs and discuss their role in plant development and evolution.

Transposable elements, contributors in the evolution of organisms (from an arms race to a source of raw materials)

There is a concept proposing that the primitive lineages of prokaryotes, eukaryotes, and viruses emerged from the primordial pool of primitive genetic elements. In this genetic pool, transposable elements (TEs) became a source of raw material for primitive genomes, tools of genetic innovation, and ancestors of modern genes (e.g. ncRNAs, tRNAs, and rRNAs). TEs contributed directly to the genome evolution of three forms of life on the earth. TEs now appear as tools that were used to giving rise to sexual dimorphism and sex determination, lineage-specific expression of genes and tissue differentiation and finally genome stability and lifespan determination.

Independent recurrent evolution of MICRORNA genes converging onto similar non-canonical organisation across green plant lineages is driven by local and segmental duplication events in species, family and lineages

The relationship between evolutionary history, organisation and transcriptional regulation of genes are intrinsically linked. These have been well studied in canonically organised protein-coding genes but not of MIRNAs. In the present study, we investigated the non-canonical arrangement of MIRNAs across taxonomic boundaries from algae to angiosperms employing a combination of genome organization, phylogeny and synteny. We retrieved the complete dataset of MIRNA from twenty-five species to identify and classify based on organisational patterns. The median size of cluster was between 2-5 kb and between 1-20 % of all MIRNAs are organized in head-to-head (with bidirectional promoter), head-to-tail (tandem), and overlapping manner. Although majority of the clusters are composed of MIRNA homologs, 25% of all clusters comprises of non-homologous genes with a potential of generating functional and regulatory complexity. A comparison of phylogeny and organizational patterns revealed that multiple independent events, some of which are species-specific, and some ancient, in different lineages, are responsible for non-canonical organization. Detailed investigation of MIR395 family across the plants revealed a complex origin of non-canonical arrangement through ancient and recent, segmental and local duplications; analysis of MIR399 family revealed major expansion occurred prior to monocot-dicot split, with few lineage-specific events. Evolution of "convergent" organization pattern of non-canonical arrangement originating from independent loci through recurrent event highlights our poor understanding of evolutionary process of MIRNA genes. The present investigation thus paves way for comparative functional genomics to understand the role of non-canonical organization on transcriptional regulation and regulatory diversity in MIRNA gene families.

Domestication of High-Copy Transposons Underlays the Wheat Small RNA Response to an Obligate Pathogen

Plant genomes have evolved several evolutionary mechanisms to tolerate and make use of transposable elements (TEs). Of these, transposon domestication into cis-regulatory and microRNA (miRNA) sequences is proposed to contribute to abiotic/biotic stress adaptation in plants. The wheat genome is derived at 85% from TEs, and contains thousands of miniature inverted-repeat transposable elements (MITEs), whose sequences are particularly prone for domestication into miRNA precursors. In this study, we investigate the contribution of TEs to the wheat small RNA immune response to the lineage-specific, obligate powdery mildew pathogen. We show that MITEs of the Mariner superfamily contribute the largest diversity of miRNAs to the wheat immune response. In particular, MITE precursors of miRNAs are wide-spread over the wheat genome, and highly conserved copies are found in the Lr34 and QPm.tut-4A mildew resistance loci. Our work suggests that transposon domestication is an important evolutionary force driving miRNA functional innovation in wheat immunity.

Plant Non-Coding RNAs: Origin, Biogenesis, Mode of Action and Their Roles in Abiotic Stress

As sessile species, plants have to deal with the rapidly changing environment. In response to these environmental conditions, plants employ a plethora of response mechanisms that provide broad phenotypic plasticity to allow the fine-tuning of the external cues related reactions. Molecular biology has been transformed by the major breakthroughs in high-throughput transcriptome sequencing and expression analysis using next-generation sequencing (NGS) technologies. These innovations have provided substantial progress in the identification of genomic regions as well as underlying basis influencing transcriptional and post-transcriptional regulation of abiotic stress response. Non-coding RNAs (ncRNAs), particularly microRNAs (miRNAs), short interfering RNAs (siRNAs), and long non-coding RNAs (lncRNAs), have emerged as essential regulators of plants abiotic stress response. However, shared traits in the biogenesis of ncRNAs and the coordinated cross-talk among ncRNAs mechanisms contribute to the complexity of these molecules and might play an essential part in regulating stress responses. Herein, we highlight the current knowledge of plant microRNAs, siRNAs, and lncRNAs, focusing on their origin, biogenesis, modes of action, and fundamental roles in plant response to abiotic stresses.

Pan-genome: A promising resource for noncoding RNA discovery in plants

Plant genomes contain both protein-coding and noncoding sequences including transposable elements (TEs) and noncoding RNAs (ncRNAs). The ncRNAs are recognized as important elements that play fundamental roles in the structural organization and function of plant genomes. Despite various hypotheses, TEs are believed to be a major precursor of ncRNAs. Transposable elements are also prime factors that cause genomic variation among members of a species. Hence, TEs pose a major challenge in the discovery and analysis of ncRNAs. With the increase in the number of sequenced plant genomes, it is now accepted that a single reference genome is insufficient to represent the complete genomic diversity and contents of a species, and exploring the pan-genome of a species is critical. In this review, we summarize the recent progress in the field of plant pan-genomes. We also discuss TEs and their roles in ncRNA biogenesis and present our perspectives on the application of pan-genomes for the discovery of ncRNAs to fully explore and exploit their biological roles in plants.

Expansion of MIR482/2118 by a class-II transposable element in cotton

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target-gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class-II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The 'GosTE' involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3-bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co-evolution of miR482/2118 and its NBS-LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.

Identification of small RNAs during cold acclimation in Arabidopsis thaliana

BACKGROUND

Cold stress causes dynamic changes in gene expression that are partially caused by small non-coding RNAs since they regulate protein coding transcripts and act in epigenetic gene silencing pathways. Thus, a detailed analysis of transcriptional changes of small RNAs (sRNAs) belonging to all known sRNA classes such as microRNAs (miRNA) and small interfering RNA (siRNAs) in response to cold contributes to an understanding of cold-related transcriptome changes.

RESULT

We subjected A. thaliana plants to cold acclimation conditions (4 °C) and analyzed the sRNA transcriptomes after 3 h, 6 h and 2 d. We found 93 cold responsive differentially expressed miRNAs and only 14 of these were previously shown to be cold responsive. We performed miRNA target prediction for all differentially expressed miRNAs and a GO analysis revealed the overrepresentation of miRNA-targeted transcripts that code for proteins acting in transcriptional regulation. We also identified a large number of differentially expressed cis- and trans-nat-siRNAs, as well as sRNAs that are derived from long non-coding RNAs. By combining the results of sRNA and mRNA profiling with miRNA target predictions and publicly available information on transcription factors, we reconstructed a cold-specific, miRNA and transcription factor dependent gene regulatory network. We verified the validity of links in the network by testing its ability to predict target gene expression under cold acclimation.

CONCLUSION

In A. thaliana, miRNAs and sRNAs derived from cis- and trans-NAT gene pairs and sRNAs derived from lncRNAs play an important role in regulating gene expression in cold acclimation conditions. This study provides a fundamental database to deepen our knowledge and understanding of regulatory networks in cold acclimation.

Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

Background

Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant.

Results

Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes.

Conclusions

Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

Rice siR109944 suppresses plant immunity to sheath blight and impacts multiple agronomic traits by affecting auxin homeostasis

Plant small RNAs (sRNAs) play significant roles in regulating various developmental processes and hormone signalling pathways involved in plant responses to a wide range of biotic and abiotic stresses. However, the functions of sRNAs in response to rice sheath blight remain unclear. We screened rice (Oryza sativa) sRNA expression patterns against Rhizoctonia solani and found that Tourist-miniature inverted-repeat transposable element (MITE)-derived small interfering RNA (siRNA) (here referred to as siR109944) expression was clearly suppressed upon R. solani infection. One potential target of siR109944 is the F-Box domain and LRR-containing protein 55 (FBL55), which encode the transport inhibitor response 1 (TIR1)-like protein. We found that rice had significantly enhanced susceptibility when siR109944 was overexpressed, while FBL55 OE plants showed resistance to R. solani challenge. Additionally, multiple agronomic traits of rice, including root length and flag leaf inclination, were affected by siR109944 expression. Auxin metabolism-related and signalling pathway-related genes were differentially expressed in the siR109944 OE and FBL55 OE plants. Importantly, pre-treatment with auxin enhanced sheath blight resistance by affecting endogenous auxin homeostasis in rice. Furthermore, transgenic Arabidopsis overexpressing siR109944 exhibited early flowering, increased tiller numbers, and increased susceptibility to R. solani. Our results demonstrate that siR109944 has a conserved function in interfering with plant immunity, growth, and development by affecting auxin homeostasis in planta. Thus, siR109944 provides a genetic target for plant breeding in the future. Furthermore, exogenous application of indole-3-acetic acid (IAA) or auxin analogues might effectively protect field crops against diseases.

Dynamic Expansion and Functional Evolutionary Profiles of Plant Conservative Gene Family SBP-Box in Twenty Two Flowering Plants and the Origin of miR156

Conservative gene families in plants, which are closely related to innovations in flowering plants, have long and complex evolutionary histories. Here, we used the SQUAMOSA promoter-binding protein (SBP-box) gene family as an example to study conservative gene families in flowering plants. In total, 11 groups, including nine angiosperm-conservative groups and two monocot- and eudicot-specific groups, were identified. Among the nine angiosperm-conservative groups, four are conserved in all land plants and the remaining five are angiosperm-specific. The five angiosperm-specific groups exhibit structural and functional diversity and evolved together, along with the evolution of flowering plants. The expansion of SBP genes was affected by miR156, and the miR156-regulated SBP genes tend to retain more copies. Our results reflect the dynamic evolutionary process of the different groups, with the identification of two genetic lines via synteny analyses. In addition, miR156 showed a close evolutionary relationship with SBP genes, suggesting that it may originate from face-to-face tandem duplication of SBP genes. SBP genes without an miR156 binding locus are usually functionally conservative or housekeeping like, belonging to the terrestrial-conservative group. In contrast, SBP genes with miR156 binding sites are selected by angiosperms to regulate more complex physiological processes.

Conserved chromosomal functions of RNA interference

RNA interference (RNAi), a cellular process through which small RNAs target and regulate complementary RNA transcripts, has well-characterized roles in post-transcriptional gene regulation and transposon repression. Recent studies have revealed additional conserved roles for RNAi proteins, such as Argonaute and Dicer, in chromosome function. By guiding chromatin modification, RNAi components promote chromosome segregation during both mitosis and meiosis and regulate chromosomal and genomic dosage response. Small RNAs and the RNAi machinery also participate in the resolution of DNA damage. Interestingly, many of these lesser-studied functions seem to be more strongly conserved across eukaryotes than are well-characterized functions such as the processing of microRNAs. These findings have implications for the evolution of RNAi since the last eukaryotic common ancestor, and they provide a more complete view of the functions of RNAi.

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/.

Evolution of small RNA expression following hybridization and allopolyploidization: insights from Spartina species (Poaceae, Chloridoideae)

Differential expression of mi-RNAs targeting developmental processes and progressive downregulation of repeat-associated siRNAs following genome merger and genome duplication in the context of allopolyploid speciation in Spartina. The role of small RNAs on gene expression regulation and genome stability is arousing increased interest and is being explored in various plant systems. In spite of prominence of reticulate evolution and polyploidy that affects the evolutionary history of all plant lineages, very few studies analysed RNAi mechanisms with this respect. Here, we explored small RNAs diversity and expression in the context of recent allopolyploid speciation, using the Spartina system, which offers a unique opportunity to explore the immediate changes following hybridization and genome duplication. Small RNA-Seq analyses were conducted on hexaploid parental species (S. alterniflora and S. maritima), their F1 hybrid S. x townsendii, and the neoallododecaploid S. anglica. We identified 594 miRNAs, 2197 miRNA-target genes, and 3730 repeat-associated siRNAs (mostly targeting Class I/Copia-Ivana- Copia-SIRE and LINEs elements). For both mi- and ra-siRNAs, we detected differential expression patterns following genome merger and genome duplication. These misregulations include non-additive expression of miRNAs in the F1 hybrid and additional changes in the allopolyploid targeting developmental processes. Expression of repeat-associated siRNAs indicates a strengthen of transposable element repression during the allopolyploidization process. Altogether, these results confirm the central role small RNAs play in shaping regulatory changes in naturally formed recent allopolyploids.

Plant microRNAs: biogenesis, gene silencing, web-based analysis tools and their use as molecular markers

MicroRNAs (miRNAs) are tiny (20-24 nt bp) regulatory non-protein-coding RNA molecules that have been extensively characterized and found important for many physiological and developmental processes. The miss-expression of miRNAs leads to various defects in plants. MicroRNAs repress gene expression by directing mRNA degradation or translational arrest. Several proteins such as PP43A, HYL1, DCL, HST are indispensable role players in promoting miRNA biogenesis in plants. During miRNA biogenesis, lariat RNAs are produced as by-products of pre-mRNA splicing which have a negative role in regulation of miRNA homeostasis. By acting as a decoy and by sequestering to the dicing complex, lariat RNA can prevent the processing of miRNAs. A number of bioinformatic tools with different methodologies are available to identify and validate miRNAs and their targets. Many miRNAs have been reported in different crops for different traits; however, no reports are available on their use in plant breeding. Recently, researchers have developed trait specific miRNA-based molecular markers (miRNA-SSRs/SNP) for many quantitative traits in different plant species. In the future, these molecular markers can be used for plant breeding programs. In this review, a comprehensive up-to-date information is provided on the bioinformatic tools used for analysis of plant miRNAs and their targets, the number of miRNAs, their biogenesis, gene silencing mechanism and miRNA-based molecular markers.

Identification of conserved and novel miRNAs responsive to heat stress in flowering Chinese cabbage using high-throughput sequencing

Plant microRNAs (miRNAs) are noncoding and endogenous key regulators that play significant functions in regulating plant responses to stress, and plant growth and development. Heat stress is a critical abiotic stress that reduces the yield and quality of flowering Chinese cabbage (Brassica campestris L. ssp. chinensis var. utilis Tsen et Lee). However, limited information is available on whether miRNAs are involved in the regulation of heat stress in B. campestris. A high-throughput sequencing approach was used to identify novel and conserved heat-responsive miRNAs in four small RNA libraries of flowering Chinese cabbage using leaves collected at 0 h, 1 h, 6 h and 12 h after a 38 °C heat-stress treatment. The analysis identified 41 conserved miRNAs (belonging to 19 MIR families), of which MIR156, MIR159, MIR168, MIR171 and MIR1885 had the most abundant molecules. Prediction and evaluation of novel miRNAs using the unannotated reads resulted in 18 candidate miRNAs. Differential expression analysis showed that most of the identified miRNAs were downregulated in heat-treated groups. To better understand functional importance, bioinformatic analysis predicted 432 unique putative target miRNAs involved in cells, cell parts, catalytic activity, cellular processes and abiotic stress responses. Furthermore, the Kyoto Encyclopedia of Genes and Genomes maps of flowering Chinese cabbage identified the significant role of miRNAs in stress adaptation and stress tolerance, and in several mitogen-activated protein kinases signaling pathways including cell death. This work presents a comprehensive study of the miRNAs for understanding the regulatory mechanisms and their participation in the heat stress of flowering Chinese cabbage.

In silico determination of transposon-derived miRNAs and targets in Aegilops species

Transposable elements (TEs) are found almost in all living organism, shaping organisms' genomes. miRNAs are noncoding RNA types which are especially important in gene expression regulations. Many previously determined plant miRNAs are identical/homologous to transposons (TE-MIR). The aim of this study was computational characterization of novel TE-related miRNAs and their targets in Aegilops genome by using stringent criteria. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed by BLAST2GO. Seventeen novel TE-related miRNAs in Aegilops genome were identified for the first time. GO analyses indicated that 40 targets played different roles in biological processes, cellular components and molecular functions. Moreover, these genes were involved in 10 metabolic pathways such as purine metabolism, nitrogen metabolism, oxidative phosphorylation, etc. as a result of KEGG analyses. Identification of miRNAs and their targets are significant to understand miRNA-TEs relationships and even how TEs affect plant growth and development. Obtaining results of this study are expected to provide possible new insight into Aegilops and its related species, wheat, with respect to miRNAs evolution and domestication.Communicated by Ramaswamy H. Sarma.

Prediction of regulatory long intergenic non-coding RNAs acting in trans through base-pairing interactions

Background

Long intergenic non-coding RNAs (lincRNAs) can act as regulators of expression of protein-coding genes. Trans-natural antisense transcripts (trans-NATs) are a type of lincRNAs that contain sequence complementary to mRNA from other loci. The regulatory potential of trans-NATs has been poorly studied in eukaryotes and no example of trans-NATs regulating gene expression in plants are reported. The goal of this study was to identify lincRNAs, and particularly trans-NATs, in Arabidopsis thaliana that have a potential to regulate expression of target genes in trans at the transcriptional or translational level.

Results

We identified 1001 lincRNAs using an RNAseq dataset from total polyA⁺ and polysome-associated RNA of seedlings grown under high and low phosphate, or shoots and roots treated with different phytohormones, of which 550 were differentially regulated. Approximately 30% of lincRNAs showed conservation amongst Brassicaceae and 25% harbored transposon element (TE) sequences. Gene co-expression network analysis highlighted a group of lincRNAs associated with the response of roots to low phosphate. A total of 129 trans-NATs were predicted, of which 88 were significantly differentially expressed under at least one pairwise comparison. Five trans-NATs showed a positive correlation between their expression and target mRNA steady-state levels, and three showed a negative correlation. Expression of four trans-NATs positively correlated with a change in target mRNA polysome association. The regulatory potential of these trans-NATs did not implicate miRNA mimics nor siRNAs. We also looked for lincRNAs that could regulate gene expression in trans by Watson-Crick DNA:RNA base pairing with target protein-encoding loci. We identified 100 and 81 with a positive or negative correlation, respectively, with steady-state level of their predicted target. The regulatory potential of one such candidate lincRNA harboring a SINE TE sequence was validated in a protoplast assay on three distinct genes containing homologous TE sequence in their promoters. Construction of networks highlighted other putative lincRNAs with multiple predicted target loci for which expression was positively correlated with target gene expression.

Conclusions

This study identified lincRNAs in Arabidopsis with potential in regulating target gene expression in trans by both RNA:RNA and RNA:DNA base pairing and highlights lincRNAs harboring TE sequences in such activity.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5946-0) contains supplementary material, which is available to authorized users.

Non-coding RNAs and transposable elements in plant genomes: emergence, regulatory mechanisms and roles in plant development and stress responses

This review will provide evidence for the indispensable function of these elements in regulating plant development and resistance to biotic and abiotic stresses, as well as their evolutionary role in facilitating plant adaptation. Over millions of years of evolution, plant genomes have acquired a complex constitution. Plant genomes consist not only of protein coding sequences, but also contain large proportions of non-coding sequences. These include introns of protein-coding genes, and intergenic sequences such as non-coding RNA, repeat sequences and transposable elements. These non-coding sequences help to regulate gene expression, and are increasingly being recognized as playing an important role in genome organization and function. In this review, we summarize the known molecular mechanisms by which gene expression is regulated by several species of non-coding RNAs (microRNAs, long non-coding RNAs, and circular RNAs) and by transposable elements. We further discuss how these non-coding RNAs and transposable elements evolve and emerge in the genome, and the potential influence and importance of these non-coding RNAs and transposable elements in plant development and in stress responses.

LtrDetector: A tool-suite for detecting long terminal repeat retrotransposons de-novo

BACKGROUND

Long terminal repeat retrotransposons are the most abundant transposons in plants. They play important roles in alternative splicing, recombination, gene regulation, and defense mechanisms. Large-scale sequencing projects for plant genomes are currently underway. Software tools are important for annotating long terminal repeat retrotransposons in these newly available genomes. However, the available tools are not very sensitive to known elements and perform inconsistently on different genomes. Some are hard to install or obsolete. They may struggle to process large plant genomes. None can be executed in parallel out of the box and very few have features to support visual review of new elements. To overcome these limitations, we developed LtrDetector, which uses techniques inspired by signal-processing.

RESULTS

We compared LtrDetector to LTR_Finder and LTRharvest, the two most successful predecessor tools, on six plant genomes. For each organism, we constructed a ground truth data set based on queries from a consensus sequence database. According to this evaluation, LtrDetector was the most sensitive tool, achieving 16-23% improvement in sensitivity over LTRharvest and 21% improvement over LTR_Finder. All three tools had low false positive rates, with LtrDetector achieving 98.2% precision, in between its two competitors. Overall, LtrDetector provides the best compromise between high sensitivity and low false positive rate while requiring moderate time and utilizing memory available on personal computers.

CONCLUSIONS

LtrDetector uses a novel methodology revolving around k-mer distributions, which allows it to produce high-quality results using relatively lightweight procedures. It is easy to install and use. It is not species specific, performing well using its default parameters on genomes of varying size and repeat content. It is automatically configured for parallel execution and runs efficiently on an ordinary personal computer. It includes a k-mer scores visualization tool to facilitate manual review of the identified elements. These features make LtrDetector an attractive tool for future annotation projects involving long terminal repeat retrotransposons.

Retrotransposon expression in response to in vitro inoculation with two fungal pathogens of Scots pine (Pinus sylvestris L.)

OBJECTIVE

Conifer genomes show high genetic diversity in intergenic regions that contain diverse sets of transposable elements with dominating long terminal repeat (LTR) retrotransposons (RE). Transcription of RE in response to environmental stimuli could produce various types of regulatory non-coding RNAs, but global genomic methylation changes could result in a coincidental expression of normally silent genomic regions. Expression of several RE families was evaluated in Scots pine seedlings after controlled inoculations with two fungal species that exhibit different modes of pathogenicity (necrotrophic and likely biotrophic); data compared to the overall RE distribution in genome. Recognition of regulatory non-coding RNA involved in host-pathogen interplay could be valuable in understanding defence mechanisms of perennial plants.

RESULTS

In the case of necrotrophic fungi Heterobasidion annosum (HA), short activation followed by restriction of RE expression was revealed after inoculation and during the spread of the pathogen. After inoculation with Lophodermium seditiosum (LS), an early increase in RE expression was revealed with the spread of the pathogen and subsequent transcription rise in all seedlings. Our observations indicate that in the complex plant genome multiple RE families constitutively express in response to pathogen invasion and these sequences could undergo regulation related to host response or pathogen influence.

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.

Comparative analysis of miniature inverted-repeat transposable elements (MITEs) and long terminal repeat (LTR) retrotransposons in six Citrus species

BACKGROUND

Miniature inverted-repeat transposable elements (MITEs) and long terminal repeat (LTR) retrotransposons are ubiquitous in plants genomes, and highly important in their evolution and diversity. However, their mechanisms of insertion/amplification and roles in Citrus genome's evolution/diversity are still poorly understood.

RESULTS

To address this knowledge gap, we developed different computational pipelines to analyze, annotate and classify MITEs and LTR retrotransposons in six different sequenced Citrus species. We identified 62,010 full-length MITEs from 110 distinguished families. We observed MITEs tend to insert in gene related regions and enriched in promoters. We found that DTM63 is possibly an active Mutator-like MITE family in the traceable past and may still be active in Citrus. The insertion of MITEs resulted in massive polymorphisms and played an important role in Citrus genome diversity and gene structure variations. In addition, 6630 complete LTR retrotransposons and 13,371 solo-LTRs were identified. Among them, 12 LTR lineages separated before the differentiation of mono- and dicotyledonous plants. We observed insertion and deletion of LTR retrotransposons was accomplished with a dynamic balance, and their half-life in Citrus was ~ 1.8 million years.

CONCLUSIONS

These findings provide insights into MITEs and LTR retrotransposons and their roles in genome diversity in different Citrus genomes.