An atlas of soybean small RNAs identifies phased siRNAs from hundreds of coding genes

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.

Unlocking the small RNAs: local and systemic modulators for advancing agronomic enhancement

Small regulatory RNAs (sRNAs) are essential regulators of gene expression across a wide range of organisms to precisely modulate gene activity based on sequence-specific recognition. In model plants like Arabidopsis thaliana, extensive research has primarily concentrated on 21 to 24-nucleotide (nt) sRNAs, particularly microRNAs (miRNAs). Recent advancements in cell and tissue isolation techniques, coupled with advanced sequencing technologies, are revealing a diverse array of preciously uncharacterized sRNA species. These include previously novel structural RNA fragments as well as numerous cell- and tissue-specific sRNAs that are active during distinct developmental stages, thereby enhancing our understanding of the precise and dynamic regulatory roles of sRNAs in plant development regulation. Additionally, a notable feature of sRNAs is their capacity for amplification and movement between cells and tissues, which facilitates long-distance communication-an adaptation critical to plants due to their sessile nature. In this review, we will discuss the classification and mechanisms of action of sRNAs, using legumes as a primary example due to their essential engagement for the unique organ establishment of root nodules and long-distance signaling, and further illustrating the potential applications of sRNAs in modern agricultural breeding and environmentally sustainable plant protection strategies.

The biogenesis, regulation and functions of transitive siRNA in plants

Small RNA (sRNA)-mediated RNA interference (RNAi) is a sequence-specific gene silencing mechanism that modulates gene expression in eukaryotes. As core molecules of RNAi, various sRNAs are encoded in the plant genome or derived from invading RNA molecules, and their biogenesis depends on distinct genetic pathways. Transitive small interfering RNAs (siRNAs), which are sRNAs produced from double-strand RNA (dsRNA) in a process that depends on RNA-dependent RNA polymerases (RDRs), can amplify and spread silencing signals to additional transcripts, thereby enabling a phenomenon termed "transitive RNAi". Members of this class of siRNAs function in various biological processes ranging from development to stress adaptation. In Arabidopsis thaliana, two RDRs participate in the generation of transitive siRNAs, acting cooperatively with various siRNA generation-related factors, such as the RNA-induced silencing complex (RISC) and aberrant RNAs. Transitive siRNAs are produced in diverse subcellular locations and structures under the control of various mechanisms, highlighting the intricacies of their biogenesis and functions. In this review, we discuss recent advances in understanding the molecular events of transitive siRNA biogenesis and its regulation, with a particular focus on factors involved in RDR recruitment. We aim to provide a comprehensive description of the generalized mechanism governing the biogenesis of transitive siRNAs. Additionally, we present an overview of the diverse biological functions of these siRNAs and raise some pressing questions in this area for further investigation.

The role of microRNAs in NBS-LRR gene expression and its implications for plant immunity and crop development

Plants evolved, over millions of years, complex defense systems against pathogens. Once infected, the interaction between pathogen effector molecules and host receptors triggers plant immune responses, which include apoptosis, systemic immune response, among others. An important protein family responsible for pathogen effector recognition is the nucleotide binding site-leucine repeat rich (NBS-LRR) proteins. The NBS-LRR gene family is the largest disease resistance gene class in plants. These proteins are widely distributed in vascular plants and have a complex multigenic cluster distribution in plant genomes. To counteract the genetic load of such a large gene family on fitness cost, plants evolved a mechanism using post transcriptional gene silencing induced by small RNAs, particularly microRNAs. For the NBS-LRR gene family, the small RNAs involved in this silencing mechanism are mainly the microRNA482/2118 superfamily. This suppression mechanism is relieved upon pathogen infection, thus allowing increased NBS-LRR expression and triggering plant immunity. In this review, we will discuss the biogenesis of microRNAs and secondary RNAs involved in this silencing mechanism, biochemical and structural features of NBS-LRR proteins in response to pathogen effectors and the evolution of microRNA-based silencing mechanism with a focus on the miR482/2118 family. Furthermore, the biotechnological manipulation of microRNA expression, using both transgenic or genome editing approaches to improve cultivated plants will be discussed, with a focus on the miR482/2118 family in soybean.

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.

Regulation of plant immunity via small RNA-mediated control of NLR expression

Plants use different receptors to detect potential pathogens: membrane-anchored pattern recognition receptors (PRRs) activated upon perception of pathogen-associated molecular patterns (PAMPs) that elicit pattern-triggered immunity (PTI); and intracellular nucleotide-binding leucine-rich repeat proteins (NLRs) activated by detection of pathogen-derived effectors, activating effector-triggered immunity (ETI). The interconnections between PTI and ETI responses have been increasingly reported. Elevated NLR levels may cause autoimmunity, with symptoms ranging from fitness cost to developmental arrest, sometimes combined with run-away cell death, making accurate control of NLR dosage key for plant survival. Small RNA-mediated gene regulation has emerged as a major mechanism of control of NLR dosage. Twenty-two nucleotide miRNAs with the unique ability to trigger secondary siRNA production from target transcripts are particularly prevalent in NLR regulation. They enhance repression of the primary NLR target, but also bring about repression of NLRs only complementary to secondary siRNAs. We summarize current knowledge on miRNAs and siRNAs in the regulation of NLR expression with an emphasis on 22 nt miRNAs and propose that miRNA and siRNA regulation of NLR levels provides additional links between PTI and NLR defense pathways to increase plant responsiveness against a broad spectrum of pathogens and control an efficient deployment of defenses.

Non-coding RNAs-mediated environmental surveillance determines male fertility in plants

Plants are continuously exposed to environmental stresses leading to significant yield losses. With the changing climatic conditions, the intensity and duration of these stresses are expected to increase, posing a severe threat to crop productivity worldwide. Male gametogenesis is one of the most sensitive developmental stages. Exposure to environmental stresses during this stage leads to male sterility and yield loss. Elucidating the underlying molecular mechanism of environment-affected male sterility is essential to address this challenge. High-throughput RNA sequencing studies, loss-of-function phenotypes of sRNA biogenesis genes and functional genomics studies with non-coding RNAs have started to unveil the roles of small RNAs, long non-coding RNAs and the complex regulatory interactions between them in regulating male fertility under different growth regimes. Here, we discuss the current understanding of the non-coding RNA-mediated environmental stress surveillance and regulation of male fertility in plants. The candidate ncRNAs emerging from these studies can be leveraged to generate environment-sensitive male sterile lines for hybrid breeding or mitigate the impact of climate change on male fertility, as the situation demands.

The role of microRNAs in responses to drought and heat stress in peanut (Arachis hypogaea)

MicroRNAs (miRNAs) are 21-24 nt small RNAs (sRNAs) that negatively regulate protein-coding genes and/or trigger phased small-interfering RNA (phasiRNA) production. Two thousand nine hundred miRNA families, of which ∼40 are deeply conserved, have been identified in ∼80 different plant species genomes. miRNA functions in response to abiotic stresses is less understood than their roles in development. Only seven peanut MIRNA families are documented in miRBase, yet a reference genome assembly is now published and over 480 plant-like MIRNA loci were predicted in the diploid peanut progenitor Arachis duranensis genome. We explored by computational analysis of a leaf sRNA library and publicly available sRNA, degradome, and transcriptome datasets the miRNA and phasiRNA space associated with drought and heat stresses in peanut. We characterized 33 novel candidate and 33 ancient conserved families of MIRNAs and present degradome evidence for their cleavage activities on mRNA targets, including several noncanonical targets and novel phasiRNA-producing noncoding and mRNA loci with validated novel targets such as miR1509 targeting serine/threonine-protein phosphatase7 and miRc20 and ahy-miR3514 targeting penta-tricopeptide repeats (PPRs), in contradistinction to other claims of miR1509/173/7122 superfamily miRNAs indirectly targeting PPRs via TAS-like noncoding RNA loci. We characterized the inverse correlations of significantly differentially expressed drought- and heat-regulated miRNAs, assayed by sRNA blots or transcriptome datasets, with target mRNA expressions in the same datasets. Meta-analysis of an expression atlas and over representation of miRNA target genes in co-expression networks suggest that miRNAs have functions in unique aspects of peanut gynophore development. Genome-wide MIRNA annotation of the published allopolyploid peanut genome can facilitate molecular breeding of value-added traits.

An insight into the roles of regulatory ncRNAs in plants: An abiotic stress and developmental perspective

Different environmental cues lead to changes in physiology, biochemistry and molecular status of plant's growth. Till date, various genes have been accounted for their role in regulating plant development and response to abiotic stress. Excluding genes that code for a functional protein in a cell, a large chunk of the eukaryotic transcriptome consists of non-coding RNAs (ncRNAs) which lack protein coding capacity but are still functional. Recent advancements in Next Generation Sequencing (NGS) technology have led to the unearthing of different types of small and large non-coding RNAs in plants. Non-coding RNAs are broadly categorised into housekeeping ncRNAs and regulatory ncRNAs which work at transcriptional, post-transcriptional and epigenetic levels. Diverse ncRNAs play different regulatory roles in nearly all biological processes including growth, development and response to changing environments. This response can be perceived and counteracted by plants using diverse evolutionarily conserved ncRNAs like miRNAs, siRNAs and lncRNAs to participate in complex molecular regimes by activating gene-ncRNA-mRNA regulatory modules to perform the downstream function. Here, we review the current understanding with a focus on recent advancements in the functional studies of the regulatory ncRNAs at the nexus of abiotic stresses and development. Also, the potential roles of ncRNAs in imparting abiotic stress tolerance and yield improvement in crop plants are also discussed with their future prospects.

The dicing activity of DCL3 and DCL4 is negatively affected by flavonoids

The dicing activities of DCL3 and DCL4 are inhibited by accumulated metabolites in soybean leaves. Epicatechin and 7,4'-dihydroxyflavone inhibited Arabidopsis DCL3 and DCL4 in vitro. Flavonoids are major secondary metabolites in plants, and soybean (Glycine max L.) is a representative plant that accumulates flavonoids, including isoflavonoids, to high levels. Naturally-occurring RNA interference (RNAi) against the chalcone synthase (CHS) gene represses flavonoid (anthocyanin) biosynthesis in an organ-specific manner, resulting in a colorless (yellow) seed coat in many soybean cultivars. To better understand seed coat-specific naturally-occurring RNAi in soybean, we characterized soybean Dicer-like (DCL) 3 and 4, which play critical roles in RNAi. Using a previously established dicing assay, two dicing activities producing 24- and 21-nt siRNAs, corresponding to DCL3 and DCL4, respectively, were detected in soybean. Dicing activity was detected in colorless seed coats where RNAi against CHS genes was found, but no dicing activity was detected in leaves where CHS expression was prevalent. Biochemical analysis revealed that soybean leaves contained two types of inhibitors effective for Arabidopsis Dicers (AtDCL3 and AtDCL4), one of which was a heat-labile high molecular weight compound of 50 to 100 kD while another was a low molecular weight substance. We found that some flavonoids, such as epicatechin and 7,4'-dihydroxyflavone, inhibited both AtDCL3 and AtDCL4, but AtDCL4 was more sensitive to these flavonoids than AtDCL3. These results suggest that flavonoids inhibit the dicing activity of DCL4 and thereby attenuate RNAi in soybean leaves.

tRNA derived small RNAs—Small players with big roles

In the past 2 decades, small non-coding RNAs derived from tRNA (tsRNAs or tRNA derived fragments; tRFs) have emerged as new powerful players in the field of small RNA mediated regulation of gene expression, translation, and epigenetic control. tRFs have been identified from evolutionarily divergent organisms from Archaea, the higher plants, to humans. Recent studies have confirmed their roles in cancers and other metabolic disorders in humans and experimental models. They have been implicated in biotic and abiotic stress responses in plants as well. In this review, we summarize the current knowledge on tRFs including types of tRFs, their biogenesis, and mechanisms of action. The review also highlights recent studies involving differential expression profiling of tRFs and elucidation of specific functions of individual tRFs from various species. We also discuss potential considerations while designing experiments involving tRFs identification and characterization and list the available bioinformatics tools for this purpose.

Identification and characterization of miRNAs and PHAS loci related to the early development of the embryo and endosperm in Fragaria × ananassa

BACKGROUND

The strawberry fleshy fruit is actually enlarged receptacle tissue, and the successful development of the embryo and endosperm is essential for receptacle fruit set. MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play indispensable regulatory roles in plant growth and development. However, miRNAs and phasiRNAs participating in the regulation of strawberry embryo and endosperm development have yet to be explored.

RESULTS

Here, we performed genome-wide identification of miRNA and phasiRNA-producing loci (PHAS) in strawberry seeds with a focus on those involved in the development of the early embryo and endosperm. We found that embryos and endosperm have different levels of small RNAs. After bioinformatics analysis, the results showed that a total of 404 miRNAs (352 known and 52 novel) and 156 PHAS genes (81 21-nt and 75 24-nt genes) could be found in strawberry seed-related tissues, of which four and nine conserved miRNA families displayed conserved expression in the endosperm and embryo, respectively. Based on refined putative annotation of PHAS loci, some auxin signal-related genes, such as CM3, TAR2, AFB2, ASA1, NAC and TAS3, were found, which demonstrates that IAA biosynthesis is important for endosperm and embryo development during early fruit growth. Additionally, some auxin signal-related conserved (miR390-TAS3) and novel (miR156-ASA1) trigger-PHAS pairs were identified.

CONCLUSIONS

Taken together, these results expand our understanding of sRNAs in strawberry embryo and endosperm development and provide a genomic resource for early-stage fruit development.

Heightened miR6024-NLR interactions facilitate necrotrophic pathogenesis in tomato

miR6024 acts as a negative regulator of R genes, hence of Tomato plant immunity, and facilitates disease by the necrotrophic pathogen A. solani. Plant resistance genes or Nucleotide-binding leucine-rich repeat (NLR) genes, integral components of plant disease stress-signaling are targeted by variable groups of miRNAs. However, the significance of miRNA-mediated regulation of NLRs during a pathogen stress response, specifically for necrotrophic fungus, is poorly understood. A thorough examination of Tomato NLRs and miRNAs could map substantial interactions of which half the annotated NLRs were targets of Solanaceae-specific and conserved miRNAs, at the NB subdomain. The Solanaceae-specific miR6024 and its NLR targets analysed in different phytopathogenic stresses revealed differential and mutually antagonistic regulation. Interestingly, miR6024-targeted cleavage of a target NLR also triggered the generation of secondary phased siRNAs which could potentially amplify the defense signal. RNA-seq analysis of leaf tissues from miR6024 overexpressing Tomato plants evidenced a perturbation in the defense transcriptome with the transgenics showing unwarranted immune response-related genes' expression with or without infection with necrotrophic Alternaria solani, though no adverse effect could be observed in the growth and development of the transgenic plants. Transgenic plants exhibited constitutive downregulation of the target NLRs, aggravated disease phenotype with an enhanced lesion, greater ROS generation and hypersusceptibility to A. solani infection, thus establishing that miR6024 negatively impacts plant immune response during necrotrophic pathogenesis. Limited knowledge about the outcome of NLR-miRNA interaction during necrotrophic pathogenesis is a hindrance to the deployment of miRNAs in crop improvement programs. With the elucidation of the necrotrophic disease-synergistic role played by miR6024, it becomes a potent candidate for biotechnological manipulation for the rapid development of pathogen-tolerant solanaceous plants.

Role of phasiRNAs in plant-pathogen interactions: molecular perspectives and bioinformatics tools

The genome of an organism is regulated in concert with the organized action of various genetic regulators at different hierarchical levels. Small non-coding RNAs are one of these regulators, among which microRNAs (miRNAs), a distinguished sRNA group with decisive functions in the development, growth and stress-responsive activities of both plants as well as animals, are keenly explored over a good number of years. Recent studies in plants revealed that apart from the silencing activity exhibited by miRNAs on their targets, miRNAs of specific size and structural features can direct the phasing pattern of their target loci to form phased secondary small interfering RNAs (phasiRNAs). These trigger-miRNAs were identified to target both coding and long non-coding RNAs that act as potent phasiRNA precursors or PHAS loci. The phasiRNAs produced thereby exhibit a role in enhancing further downstream regulation either on their own precursors or on those transcripts that are distinct from their genetic source of origin. Hence, these tiny regulators can stimulate an elaborative cascade of interacting RNA networks via cis and trans-regulatory mechanisms. Our review focuses on the comprehensive understanding of phasiRNAs and their trigger miRNAs, by giving much emphasis on their role in the regulation of plant defense responses, together with a summary of the computational tools available for the prediction of the same.

microRNAs and Their Roles in Plant Development

Small RNAs are short non-coding RNAs with a length ranging between 20 and 24 nucleotides. Of these, microRNAs (miRNAs) play a distinct role in plant development. miRNAs control target gene expression at the post-transcriptional level, either through direct cleavage or inhibition of translation. miRNAs participate in nearly all the developmental processes in plants, such as juvenile-to-adult transition, shoot apical meristem development, leaf morphogenesis, floral organ formation, and flowering time determination. This review summarizes the research progress in miRNA-mediated gene regulation and its role in plant development, to provide the basis for further in-depth exploration regarding the function of miRNAs and the elucidation of the molecular mechanism underlying the interaction of miRNAs and other pathways.

Progress in soybean functional genomics over the past decade

Soybean is one of the most important oilseed and fodder crops. Benefiting from the efforts of soybean breeders and the development of breeding technology, large number of germplasm has been generated over the last 100 years. Nevertheless, soybean breeding needs to be accelerated to meet the needs of a growing world population, to promote sustainable agriculture and to address future environmental changes. The acceleration is highly reliant on the discoveries in gene functional studies. The release of the reference soybean genome in 2010 has significantly facilitated the advance in soybean functional genomics. Here, we review the research progress in soybean omics (genomics, transcriptomics, epigenomics and proteomics), germplasm development (germplasm resources and databases), gene discovery (genes that are responsible for important soybean traits including yield, flowering and maturity, seed quality, stress resistance, nodulation and domestication) and transformation technology during the past decade. At the end, we also briefly discuss current challenges and future directions.

Degradome comparison between wild and cultivated rice identifies differential targeting by miRNAs

Background

Small non-coding (s)RNAs are involved in the negative regulation of gene expression, playing critical roles in genome integrity, development and metabolic pathways. Targeting of RNAs by ribonucleoprotein complexes of sRNAs bound to Argonaute (AGO) proteins results in cleaved RNAs having precise and predictable 5` ends. While tools to study sliced bits of RNAs to confirm the efficiency of sRNA-mediated regulation are available, they are sub-optimal. In this study, we provide an improvised version of a tool with better efficiency to accurately validate sRNA targets.

Results

Here, we improvised the CleaveLand tool to identify additional micro (mi)RNA targets that belong to the same family and also other targets within a specified free energy cut-off. These additional targets were otherwise excluded during the default run. We employed these tools to understand the sRNA targeting efficiency in wild and cultivated rice, sequenced degradome from two rice lines, O. nivara and O. sativa indica Pusa Basmati-1 and analyzed variations in sRNA targeting. Our results indicate the existence of multiple miRNA-mediated targeting differences between domesticated and wild species. For example, Os5NG4 was targeted only in wild rice that might be responsible for the poor secondary wall formation when compared to cultivated rice. We also identified differential mRNA targets of secondary sRNAs that were generated after miRNA-mediated cleavage of primary targets.

Conclusions

We identified many differentially targeted mRNAs between wild and domesticated rice lines. In addition to providing a step-wise guide to generate and analyze degradome datasets, we showed how domestication altered sRNA-mediated cascade silencing during the evolution of indica rice.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08288-5.

Epigenetic Mechanisms of Senescence in Plants

Senescence is a major developmental transition in plants that requires a massive reprogramming of gene expression and includes various layers of regulations. Senescence is either an age-dependent or a stress-induced process, and is under the control of complex regulatory networks that interact with each other. It has been shown that besides genetic reprogramming, which is an important aspect of plant senescence, transcription factors and higher-level mechanisms, such as epigenetic and small RNA-mediated regulators, are also key factors of senescence-related genes. Epigenetic mechanisms are an important layer of this multilevel regulatory system that change the activity of transcription factors (TFs) and play an important role in modulating the expression of senescence-related gene. They include chromatin remodeling, DNA methylation, histone modification, and the RNA-mediated control of transcription factors and genes. This review provides an overview of the known epigenetic regulation of plant senescence, which has mostly been studied in the form of leaf senescence, and it also covers what has been reported about whole-plant senescence.

Fab Advances in Fabaceae for Abiotic Stress Resilience: From 'Omics' to Artificial Intelligence

Legumes are a better source of proteins and are richer in diverse micronutrients over the nutritional profile of widely consumed cereals. However, when exposed to a diverse range of abiotic stresses, their overall productivity and quality are hugely impacted. Our limited understanding of genetic determinants and novel variants associated with the abiotic stress response in food legume crops restricts its amelioration. Therefore, it is imperative to understand different molecular approaches in food legume crops that can be utilized in crop improvement programs to minimize the economic loss. 'Omics'-based molecular breeding provides better opportunities over conventional breeding for diversifying the natural germplasm together with improving yield and quality parameters. Due to molecular advancements, the technique is now equipped with novel 'omics' approaches such as ionomics, epigenomics, fluxomics, RNomics, glycomics, glycoproteomics, phosphoproteomics, lipidomics, regulomics, and secretomics. Pan-omics-which utilizes the molecular bases of the stress response to identify genes (genomics), mRNAs (transcriptomics), proteins (proteomics), and biomolecules (metabolomics) associated with stress regulation-has been widely used for abiotic stress amelioration in food legume crops. Integration of pan-omics with novel omics approaches will fast-track legume breeding programs. Moreover, artificial intelligence (AI)-based algorithms can be utilized for simulating crop yield under changing environments, which can help in predicting the genetic gain beforehand. Application of machine learning (ML) in quantitative trait loci (QTL) mining will further help in determining the genetic determinants of abiotic stress tolerance in pulses.

Uncovering miRNA-mRNA Regulatory Modules in Developing Xylem of Pinus massoniana via Small RNA and Degradome Sequencing

Xylem is required for the growth and development of higher plants to provide water and mineral elements. The thickening of the xylem secondary cell wall (SCW) not only improves plant survival, but also provides raw materials for industrial production. Numerous studies have found that transcription factors and non-coding RNAs regulate the process of SCW thickening. Pinus massoniana is an important woody tree species in China and is widely used to produce materials for construction, furniture, and packaging. However, the target genes of microRNAs (miRNAs) in the developing xylem of P. massoniana are not known. In this study, a total of 25 conserved miRNAs and 173 novel miRNAs were identified via small RNA sequencing, and 58 differentially expressed miRNAs were identified between the developing xylem (PM_X) and protoplasts isolated from the developing xylem (PM_XP); 26 of these miRNAs were significantly up-regulated in PM_XP compared with PM_X, and 32 were significantly down-regulated. A total of 153 target genes of 20 conserved miRNAs and 712 target genes of 113 novel miRNAs were verified by degradome sequencing. There may be conserved miRNA-mRNA modules (miRNA-MYB, miRNA-ARF, and miRNA-LAC) involved in softwood and hardwood formation. The results of qRT-PCR-based parallel validation were in relatively high agreement. This study explored the potential regulatory network of miRNAs in the developing xylem of P. massoniana and provides new insights into wood formation in coniferous species.

MiR1885 Regulates Disease Tolerance Genes in Brassica rapa during Early Infection with Plasmodiophora brassicae

Clubroot caused by Plasmodiophora brassicae is a severe disease of cruciferous crops that decreases crop quality and productivity. Several clubroot resistance-related quantitative trait loci and candidate genes have been identified. However, the underlying regulatory mechanism, the interrelationships among genes, and how genes are regulated remain unexplored. MicroRNAs (miRNAs) are attracting attention as regulators of gene expression, including during biotic stress responses. The main objective of this study was to understand how miRNAs regulate clubroot resistance-related genes in P. brassicae-infected Brassica rapa. Two Brassica miRNAs, Bra-miR1885a and Bra-miR1885b, were revealed to target TIR-NBS genes. In non-infected plants, both miRNAs were expressed at low levels to maintain the balance between plant development and basal immunity. However, their expression levels increased in P. brassicae-infected plants. Both miRNAs down-regulated the expression of the TIR-NBS genes Bra019412 and Bra019410, which are located at a clubroot resistance-related quantitative trait locus. The Bra-miR1885-mediated down-regulation of both genes was detected for up to 15 days post-inoculation in the clubroot-resistant line CR Shinki and in the clubroot-susceptible line 94SK. A qRT-PCR analysis revealed Bra019412 expression was negatively regulated by miR1885. Both Bra019412 and Bra019410 were more highly expressed in CR Shinki than in 94SK; the same expression pattern was detected in multiple clubroot-resistant and clubroot-susceptible inbred lines. A 5′ rapid amplification of cDNA ends analysis confirmed the cleavage of Bra019412 by Bra-miR1885b. Thus, miR1885s potentially regulate TIR-NBS gene expression during P. brassicae infections of B. rapa.

Pre-meiotic 21-nucleotide reproductive phasiRNAs emerged in seed plants and diversified in flowering plants

Plant small RNAs are important regulatory elements that fine-tune gene expression and maintain genome integrity by silencing transposons. Reproductive organs of monocots produce abundant phased, small interfering RNAs (phasiRNAs). The 21-nt reproductive phasiRNAs triggered by miR2118 are highly enriched in pre-meiotic anthers, and have been found in multiple eudicot species, in contrast with prior reports of monocot specificity. The 24-nt reproductive phasiRNAs are triggered by miR2275, and are highly enriched during meiosis in many angiosperms. Here, we report the widespread presence of the 21-nt reproductive phasiRNA pathway in eudicots including canonical and non-canonical microRNA (miRNA) triggers of this pathway. In eudicots, these 21-nt phasiRNAs are enriched in pre-meiotic stages, a spatiotemporal distribution consistent with that of monocots and suggesting a role in anther development. Although this pathway is apparently absent in well-studied eudicot families including the Brassicaceae, Solanaceae and Fabaceae, our work in eudicots supports an earlier singular finding in spruce, a gymnosperm, indicating that the pathway of 21-nt reproductive phasiRNAs emerged in seed plants and was lost in some lineages.

Non-Coding RNAs in Legumes: Their Emerging Roles in Regulating Biotic/Abiotic Stress Responses and Plant Growth and Development

Noncoding RNAs, including microRNAs (miRNAs), small interference RNAs (siRNAs), circular RNA (circRNA), and long noncoding RNAs (lncRNAs), control gene expression at the transcription, post-transcription, and translation levels. Apart from protein-coding genes, accumulating evidence supports ncRNAs playing a critical role in shaping plant growth and development and biotic and abiotic stress responses in various species, including legume crops. Noncoding RNAs (ncRNAs) interact with DNA, RNA, and proteins, modulating their target genes. However, the regulatory mechanisms controlling these cellular processes are not well understood. Here, we discuss the features of various ncRNAs, including their emerging role in contributing to biotic/abiotic stress response and plant growth and development, in addition to the molecular mechanisms involved, focusing on legume crops. Unravelling the underlying molecular mechanisms and functional implications of ncRNAs will enhance our understanding of the coordinated regulation of plant defences against various biotic and abiotic stresses and for key growth and development processes to better design various legume crops for global food security.

TarDB: an online database for plant miRNA targets and miRNA-triggered phased siRNAs

Background

In plants, microRNAs (miRNAs) are pivotal regulators of plant development and stress responses. Different computational tools and web servers have been developed for plant miRNA target prediction; however, in silico prediction normally contains false positive results. In addition, many plant miRNA target prediction servers lack information for miRNA-triggered phased small interfering RNAs (phasiRNAs). Creating a comprehensive and relatively high-confidence plant miRNA target database is much needed.

Results

Here, we report TarDB, an online database that collects three categories of relatively high-confidence plant miRNA targets: (i) cross-species conserved miRNA targets; (ii) degradome/PARE (Parallel Analysis of RNA Ends) sequencing supported miRNA targets; (iii) miRNA-triggered phasiRNA loci. TarDB provides a user-friendly interface that enables users to easily search, browse and retrieve miRNA targets and miRNA initiated phasiRNAs in a broad variety of plants. TarDB has a comprehensive collection of reliable plant miRNA targets containing previously unreported miRNA targets and miRNA-triggered phasiRNAs even in the well-studied model species. Most of these novel miRNA targets are relevant to lineage-specific or species-specific miRNAs. TarDB data is freely available at http://www.biosequencing.cn/TarDB.

Conclusions

In summary, TarDB serves as a useful web resource for exploring relatively high-confidence miRNA targets and miRNA-triggered phasiRNAs in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07680-5.

Dynamics of miRNA mediated regulation of legume symbiosis

Symbiotic nitrogen fixation in legume nodules is important in soils with low nitrogen availability. The initiation and sustainability of symbiosis require cellular reprogramming that involves the miRNA-mediated inhibition or activation of specific nodulation genes. The high-throughput sequencing of small RNA libraries has identified miRNAs and their targets, which are the major players in the post-transcriptional gene regulation (PTGS) of the different stages of legume-rhizobia symbiosis ranging from bacterial colonization and organogenesis to symbiotic nitrogen fixation. Here, we present an overview of information obtained from the miRNA libraries from nodulating tissues that have been sequenced to date. The functional analysis of miRNAs has revealed roles in phytohormone homeostasis and spatio-temporal regulation, as well as the mobility of miRNAs and their functions in shoot to root signalling that affects diverse functions, including bacterial entry, meristem division and differentiation, nitrogen fixation and senescence. Furthermore, small RNA fragments of rhizobial origin repress complementary plant mRNAs. We also consider the roles of miRNAs in determinate or indeterminate nodules. Taken together, this overview confirms that miRNAs are master regulators of the legume-rhizobia symbiosis.

Title: Hypermethylation of miRNA Genes During Nodule Development

DNA methylation has recently emerged as a powerful regulatory mechanism controlling the expression of key regulators of various developmental processes, including nodulation. However, the functional role of DNA methylation in regulating the expression of microRNA (miRNA) genes during the formation and development of nitrogen-fixing nodules remains largely unknown. In this study, we profiled DNA methylation patterns of miRNA genes during nodule formation, development, and early senescence stages in soybean (Glycine max) through the analysis of methylC-seq data. Absolute DNA methylation levels in the CG, CHH, and CHH sequence contexts over the promoter and primary transcript regions of miRNA genes were significantly higher in the nodules compared with the corresponding root tissues at these three distinct nodule developmental stages. We identified a total of 82 differentially methylated miRNAs in the nodules compared with roots. Differential DNA methylation of these 82 miRNAs was detected only in the promoter (69), primary transcript region (3), and both in the promoter and primary transcript regions (10). The large majority of these differentially methylated miRNAs were hypermethylated in nodules compared with the corresponding root tissues and were found mainly in the CHH context and showed stage-specific methylation patterns. Differentially methylated regions in the promoters of 25 miRNAs overlapped with transposable elements, a finding that may explain the vulnerability of miRNAs to DNA methylation changes during nodule development. Gene expression analysis of a set of promoter-differentially methylated miRNAs pointed to a negative association between DNA methylation and miRNA expression. Gene Ontology and pathways analyses indicate that changes in DNA methylation of miRNA genes are reprogrammed and contribute to nodule development through indirect regulation of genes involved in cellular processes and pathways with well-established roles in nodulation.

Unidirectional movement of small RNAs from shoots to roots in interspecific heterografts

Long-distance RNA movement is important for plant growth and environmental responses; however, the extent to which RNAs move between distant tissues, their relative magnitude and functional importance remain to be elucidated on a genomic scale. Using a soybean (Glycine max)-common bean (Phaseolus vulgaris) grafting system, we identified 100 shoot-root mobile microRNAs and 32 shoot-root mobile phased secondary small interfering RNAs (phasiRNAs), which were predominantly produced in shoots and transported to roots, and, in most cases, accumulated to a level similar to that observed in shoots. Many of these microRNAs or phasiRNAs enabled cleavage of their messenger RNA targets or phasiRNA precursors in roots. In contrast, most mobile-capable mRNAs were transcribed in both shoots and roots, with only small proportions transported to recipient tissues. These findings suggest that the regulatory mechanisms for small RNA movement are different from those for mRNA movement, and that the former is more strictly regulated and, probably, more functionally important than the latter.

Soybean DICER-LIKE2 Regulates Seed Coat Color via Production of Primary 22-Nucleotide Small Interfering RNAs from Long Inverted Repeats

In plants, 22-nucleotide small RNAs trigger the production of secondary small interfering RNAs (siRNAs) and enhance silencing. DICER-LIKE2 (DCL2)-dependent 22-nucleotide siRNAs are rare in Arabidopsis (Arabidopsis thaliana) and are thought to function mainly during viral infection; by contrast, these siRNAs are abundant in many crops such as soybean (Glycine max) and maize (Zea mays). Here, we studied soybean 22-nucleotide siRNAs by applying CRISPR-Cas9 to simultaneously knock out the two copies of soybean DCL2, GmDCL2a and GmDCL2b, in the Tianlong1 cultivar. Small RNA sequencing revealed that most 22-nucleotide siRNAs are derived from long inverted repeats (LIRs) and disappeared in the Gmdcl2a/2b double mutant. De novo assembly of a Tianlong1 reference genome and transcriptome profiling identified an intronic LIR formed by the chalcone synthase (CHS) genes CHS1 and CHS3 This LIR is the source of primary 22-nucleotide siRNAs that target other CHS genes and trigger the production of secondary 21-nucleotide siRNAs. Disruption of this process in Gmdcl2a/2b mutants substantially increased CHS mRNA levels in the seed coat, thus changing the coat color from yellow to brown. Our results demonstrated that endogenous LIR-derived transcripts in soybean are predominantly processed by GmDCL2 into 22-nucleotide siRNAs and uncovered a role for DCL2 in regulating natural traits.

Identification of Differentially Methylated miRNA Genes During Compatible and Incompatible Interactions Between Soybean and Soybean Cyst Nematode

DNA methylation is a widespread epigenetic mark that affects gene expression and transposon mobility during plant development and stress responses. However, the role of DNA methylation in regulating the expression of microRNA (miRNA) genes remains largely unexplored. Here, we analyzed DNA methylation changes of miRNA genes using a pair of soybean (Glycine max) near-isogenic lines (NILs) differing in their response to soybean cyst nematode (SCN; Heterodera glycines). Differences in global DNA methylation levels over miRNA genes in response to SCN infection were observed between the isogenic lines. miRNA genes with significant changes in DNA methylation levels in the promoter and primary transcript-coding regions were detected in both lines. In the susceptible isogenic line (NIL-S), 82 differentially methylated miRNAs were identified in response to SCN infection whereas, in the resistant isogenic line (NIL-R), only 16 differentially methylated miRNAs were identified. Interestingly, gma-miR5032, gma-miR5043, gma-miR1520b, and gma-2107-ch16 showed opposite methylation patterns in the isogenic lines. In addition, the miRNA paralogs gma-miR5770a and gma-miR5770b showed hypermethylation and hypomethylation in NIL-S and NIL-R, respectively. Gene expression quantification of gma-miR5032, gma-miR5043, gma-miR1520b, and gma-miR5770a/b and their confirmed targets indicated a role of DNA methylation in regulating miRNA expression and, thus, their targets upon SCN infection. Furthermore, overexpression of these four miRNAs in NIL-S using transgenic hairy root system enhanced plant resistance to SCN to various degrees with a key role observed for miR5032. Together, our results provide new insights into the role of epigenetic mechanisms in controlling miRNA regulatory function during SCN-soybean interactions.[Formula: see text] Copyright © 2020 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Genome-wide identification of microRNAs and phased siRNAs in soybean roots under long-term salt stress

BACKGROUND

Salinity stress, as the key limiting factor for agricultural productivity, can activate a series of molecular responses and alter gene expression in plants. Endogenous regulatory small RNAs, such as microRNAs (miRNAs) and phased siRNAs (phasiRNAs), play crucial roles during stress adaptation and prevent the injury from environmental circumstances.

OBJECTIVE

To identify long-term salt stress responsive miRNAs and phasiRNAs as well as their associated genes and pathways in soybean roots.

METHODS

Small RNA and degradome sequencing strategies were applied to genome widely investigate miRNAs and phasiRNAs in soybean roots under control and long-term salt stress conditions.

RESULTS

In this study, stringent bioinformatic analysis led to the identification of 253 conserved and 38 novel miRNA candidates. Results of expression profiling, target and endogenous target mimics predictions provided valuable clues to their functional roles. Furthermore, 156 genes were identified to be capable of generating 21 nt and 24 nt phasiRNAs, in which 37 candidates were confirmed by degradome data for miRNA-directed cleavage. Approximately 90% of these phasiRNA loci were protein coding genes. And GO enrichment analysis pointed to "signal transduction" and "ADP binding" entries and reflected the functional roles of identified phasiRNA genes.

CONCLUSION

Taken together, our findings extended the knowledge of salt responsive miRNAs and phasiRNAs in soybean roots, and provided valuable information for a better understanding of the regulatory events caused by small RNAs underlying plant adaptations to long-term salt stress.

Progresses, Challenges, and Prospects of Genome Editing in Soybean (Glycine max)

Soybean is grown worldwide for oil and protein source as food, feed and industrial raw material for biofuel. Steady increase in soybean production in the past century mainly attributes to genetic mediation including hybridization, mutagenesis and transgenesis. However, genetic resource limitation and intricate social issues in use of transgenic technology impede soybean improvement to meet rapid increases in global demand for soybean products. New approaches in genomics and development of site-specific nucleases (SSNs) based genome editing technologies have expanded soybean genetic variations in its germplasm and have potential to make precise modification of genes controlling the important agronomic traits in an elite background. ZFNs, TALENS and CRISPR/Cas9 have been adapted in soybean improvement for targeted deletions, additions, replacements and corrections in the genome. The availability of reference genome assembly and genomic resources increases feasibility in using current genome editing technologies and their new development. This review summarizes the status of genome editing in soybean improvement and future directions in this field.

PhasiRNAs in Plants: Their Biogenesis, Genic Sources, and Roles in Stress Responses, Development, and Reproduction

Phased secondary small interfering RNAs (phasiRNAs) constitute a major category of small RNAs in plants, but most of their functions are still poorly defined. Some phasiRNAs, known as trans-acting siRNAs, are known to target complementary mRNAs for degradation and to function in development. However, the targets or biological roles of other phasiRNAs remain speculative. New insights into phasiRNA biogenesis, their conservation, and their variation across the flowering plants continue to emerge due to the increased availability of plant genomic sequences, deeper and more sophisticated sequencing approaches, and improvements in computational biology and biochemical/molecular/genetic analyses. In this review, we survey recent progress in phasiRNA biology, with a particular focus on two classes associated with male reproduction: 21-nucleotide (accumulate early in anther ontogeny) and 24-nucloetide (produced in somatic cells during meiosis) phasiRNAs. We describe phasiRNA biogenesis, function, and evolution and define the unanswered questions that represent topics for future research.

Epigenetics: possible applications in climate-smart crop breeding

To better adapt transiently or lastingly to stimuli from the surrounding environment, the chromatin states in plant cells vary to allow the cells to fine-tune their transcriptional profiles. Modifications of chromatin states involve a wide range of post-transcriptional histone modifications, histone variants, DNA methylation, and activity of non-coding RNAs, which can epigenetically determine specific transcriptional outputs. Recent advances in the area of '-omics' of major crops have facilitated identification of epigenetic marks and their effect on plant response to environmental stresses. As most epigenetic mechanisms are known from studies in model plants, we summarize in this review recent epigenetic studies that may be important for improvement of crop adaptation and resilience to environmental changes, ultimately leading to the generation of stable climate-smart crops. This has paved the way for exploitation of epigenetic variation in crop breeding.

Identification of tRFs and phasiRNAs in tomato (Solanum lycopersicum) and their responses to exogenous abscisic acid

BACKGROUND

The non-coding small RNA tRFs (tRNA-derived fragments) and phasiRNAs (plant-specific) exert important roles in plant growth, development and stress resistances. However, whether the tRFs and phasiRNAs respond to the plant important stress hormone abscisic acid (ABA) remain enigma.

RESULTS

Here, the RNA-sequencing was implemented to decipher the landscape of tRFs and phasiRNAs in tomato (Solanum lycopersicum) leaves and their responses when foliar spraying exogenous ABA after 24 h. In total, 733 tRFs and 137 phasiRNAs were detected. The tRFs were mainly derived from the tRNA^Ala transporting alanine, which tended to be cleaved at the 5^'terminal guanine site and D loop uracil site to produce tRF^Ala with length of 20 nt. Most of phasiRNAs originated from NBS-LRR resistance genes. Expression analysis revealed that 156 tRFs and 68 phasiRNAs expressed differentially, respectively. Generally, exogenous ABA mainly inhibited the expression of tRFs and phasiRNAs. Furthermore, integrating analysis of target gene prediction and transcriptome data presented that ABA significantly downregulated the abundance of phsaiRNAs associated with biological and abiotic resistances. Correspondingly, their target genes such as AP2/ERF, WRKY and NBS-LRR, STK and RLK, were mainly up-regulated.

CONCLUSIONS

Combined with the previous analysis of ABA-response miRNAs, it was speculated that ABA can improve the plant resistances to various stresses by regulating the expression and interaction of small RNAs (such as miRNAs, tRFs, phasiRNAs) and their target genes. This study enriches the plant tRFs and phasiRNAs, providing a vital basis for further investigating ABA response-tRFs and phasiRNAs and their functions in biotic and abiotic stresses.

Molecular and genetic bases of heat stress responses in crop plants and breeding for increased resilience and productivity

To ensure the food security of future generations and to address the challenge of the 'no hunger zone' proposed by the FAO (Food and Agriculture Organization), crop production must be doubled by 2050, but environmental stresses are counteracting this goal. Heat stress in particular is affecting agricultural crops more frequently and more severely. Since the discovery of the physiological, molecular, and genetic bases of heat stress responses, cultivated plants have become the subject of intense research on how they may avoid or tolerate heat stress by either using natural genetic variation or creating new variation with DNA technologies, mutational breeding, or genome editing. This review reports current understanding of the genetic and molecular bases of heat stress in crops together with recent approaches to creating heat-tolerant varieties. Research is close to a breakthrough of global relevance, breeding plants fitter to face the biggest challenge of our time.

Integrated annotations and analyses of small RNA-producing loci from 47 diverse plants

Plant endogenous small RNAs (sRNAs) are important regulators of gene expression. There are two broad categories of plant sRNAs: microRNAs (miRNAs) and endogenous short interfering RNAs (siRNAs). MicroRNA loci are relatively well-annotated but compose only a small minority of the total sRNA pool; siRNA locus annotations have lagged far behind. Here, we used a large data set of published and newly generated sRNA sequencing data (1333 sRNA-seq libraries containing more than 20 billion reads) and a uniform bioinformatic pipeline to produce comprehensive sRNA locus annotations of 47 diverse plants, yielding more than 2.7 million sRNA loci. The two most numerous classes of siRNA loci produced mainly 24- and 21-nucleotide (nt) siRNAs, respectively. Most often, 24-nt-dominated siRNA loci occurred in intergenic regions, especially at the 5′-flanking regions of protein-coding genes. In contrast, 21-nt-dominated siRNA loci were most often derived from double-stranded RNA precursors copied from spliced mRNAs. Genic 21-nt-dominated loci were especially common from disease resistance genes, including from a large number of monocots. Individual siRNA sequences of all types showed very little conservation across species, whereas mature miRNAs were more likely to be conserved. We developed a web server where our data and several search and analysis tools are freely accessible.

Conservation and Divergence in the Meiocyte sRNAomes of Arabidopsis, Soybean, and Cucumber

Meiosis is a critical process for sexual reproduction. During meiosis, genetic information on homologous chromosomes is shuffled through meiotic recombination to produce gametes with novel allelic combinations. Meiosis and recombination are orchestrated by several mechanisms including regulation by small RNAs (sRNAs). Our previous work in Arabidopsis (Arabidopsis thaliana) meiocytes showed that meiocyte-specific sRNAs (ms-sRNAs) have distinct characteristics, including positive association with the coding region of genes that are transcriptionally upregulated during meiosis. Here, we characterized the ms-sRNAs in two important crops, soybean (Glycine max) and cucumber (Cucumis sativus). Ms-sRNAs in soybean have the same features as those in Arabidopsis, suggesting that they may play a conserved role in eudicots. We also investigated the profiles of microRNAs (miRNAs) and phased secondary small interfering RNAs in the meiocytes of all three species. Two conserved miRNAs, miR390 and miR167, are highly abundant in the meiocytes of all three species. In addition, we identified three novel cucumber miRNAs. Intriguingly, our data show that the previously identified phased secondary small interfering RNA pathway involving soybean-specific miR4392 is more abundant in meiocytes. These results showcase the conservation and divergence of ms-sRNAs in flowering plants, and broaden our understanding of sRNA function in crop species.

Identification of novel phasiRNAs loci on long non-coding RNAs in Arabidopsis thaliana

Long non-coding RNAs (lncRNAs) are the "dark matters"involved in gene regulation with complex mechanisms. However, the functions of most lncRNAs remain to be determined. Our previous work revealed a massive number of degradome-supported cleavage signatures on Arabidopsis lncRNAs. Some of them have been confirmed associated with miRNAs-like sRNAs production, while others without long stem structure remain unexplored. A systematical search for phasiRNAs generating ability of these lncRNAs was conducted. Eight novel small RNA triggered lncRNA-phasiRNA pathways were discovered and three of them were found to be conserved in Arabidopsis, Oryza sativa, Glycine max and Gossypium hirsutum. Besides, Five novel ta-siRNAs derived from these lncRNAs were further identified to be involved in the regulation of plant development, stress responses and aromatic amino acids synthesis. These results substantially expanded the gene regulation mechanisms of lncRNAs.

Genome-Wide Analysis of Cotton miRNAs During Whitefly Infestation Offers New Insights into Plant-Herbivore Interaction

Although the regulatory function of miRNAs and their targets have been characterized in model plants, a possible underlying role in the cotton response to herbivore infestation has not been determined. To investigate this, we performed small RNA and degradome sequencing between resistant and susceptible cotton cultivar following infestation with the generalist herbivore whitefly. In total, the 260 miRNA families and 241 targets were identified. Quantitative-PCR analysis revealed that several miRNAs and their corresponding targets exhibited dynamic spatio-temporal expression patterns. Moreover, 17 miRNA precursors were generated from 29 long intergenic non-coding RNA (lincRNA) transcripts. The genome-wide analysis also led to the identification of 85 phased small interfering RNA (phasiRNA) loci. Among these, nine PHAS genes were triggered by miR167, miR390, miR482a, and two novel miRNAs, including those encoding a leucine-rich repeat (LRR) disease resistance protein, an auxin response factor (ARF) and MYB transcription factors. Through combined modeling and experimental data, we explored and expanded the miR390-tasiARF cascade during the cotton response to whitefly. Virus-induced gene silencing (VIGS) of ARF8 from miR390 target in whitefly-resistant cotton plants increased auxin and jasmonic acid (JA) accumulation, resulting in increased tolerance to whitefly infestation. These results highlight the provides a useful transcriptomic resource for plant-herbivore interaction.

miR828 and miR858 regulate VvMYB114 to promote anthocyanin and flavonol accumulation in grapes

MicroRNAs are a class of non-coding small RNAs involved in the negative regulation of gene expression, which play critical roles in developmental and metabolic pathways. Studies in several plants have identified a few microRNAs and other small RNAs that target regulators of the phenylpropanoid metabolic pathway called the MYB transcription factors. However, it is not well understood how sRNA-mediated regulation of MYBs influences the accumulation of specific secondary metabolites. Using sRNA sequencing, degradome analysis, mRNA sequencing, and proteomic analysis, we establish that grape lines with high anthocyanin content express two MYB-targeting microRNAs abundantly, resulting in the differential expression of specific MYB proteins. miR828 and miR858 target coding sequences of specific helix motifs in the mRNA sequences of MYB proteins. Targeting by miR828 caused MYB RNA decay and the production of a cascade of secondary siRNAs that depend on RNA-dependent RNA polymerase 6. MYB suppression and cascade silencing was more robust in grape lines with high anthocyanin content than in a flavonol-rich grape line. We establish that microRNA-mediated silencing targeted the repressor class of MYBs to promote anthocyanin biosynthesis in grape lines with high anthocyanins. We propose that this process regulates the expression of appropriate MYBs in grape lines to produce specific secondary metabolites.

Convergent horizontal gene transfer and cross-talk of mobile nucleic acids in parasitic plants

Horizontal gene transfer (HGT), the movement and genomic integration of DNA across species boundaries, is commonly associated with bacteria and other microorganisms, but functional HGT (fHGT) is increasingly being recognized in heterotrophic parasitic plants that obtain their nutrients and water from their host plants through direct haustorial feeding. Here, in the holoparasitic stem parasite Cuscuta, we identify 108 transcribed and probably functional HGT events in Cuscuta campestris and related species, plus 42 additional regions with host-derived transposon, pseudogene and non-coding sequences. Surprisingly, 18 Cuscuta fHGTs were acquired from the same gene families by independent HGT events in Orobanchaceae parasites, and the majority are highly expressed in the haustorial feeding structures in both lineages. Convergent retention and expression of HGT sequences suggests an adaptive role for specific additional genes in parasite biology. Between 16 and 20 of the transcribed HGT events are inferred as ancestral in Cuscuta based on transcriptome sequences from species across the phylogenetic range of the genus, implicating fHGT in the successful radiation of Cuscuta parasites. Genome sequencing of C. campestris supports transfer of genomic DNA-rather than retroprocessed RNA-as the mechanism of fHGT. Many of the C. campestris genes horizontally acquired are also frequent sources of 24-nucleotide small RNAs that are typically associated with RNA-directed DNA methylation. One HGT encoding a leucine-rich repeat protein kinase overlaps with a microRNA that has been shown to regulate host gene expression, suggesting that HGT-derived parasite small RNAs may function in the parasite-host interaction. This study enriches our understanding of HGT by describing a parasite-host system with unprecedented gene exchange that points to convergent evolution of HGT events and the functional importance of horizontally transferred coding and non-coding sequences.

Rhizobial tRNA-derived small RNAs are signal molecules regulating plant nodulation

Rhizobial infection and root nodule formation in legumes require recognition of signal molecules produced by the bacteria and their hosts. Here, we show that rhizobial transfer RNA (tRNA)-derived small RNA fragments (tRFs) are signal molecules that modulate host nodulation. Three families of rhizobial tRFs were confirmed to regulate host genes associated with nodule initiation and development through hijacking the host RNA-interference machinery that involves ARGONAUTE 1. Silencing individual tRFs with the use of short tandem target mimics or by overexpressing their targets represses root hair curling and nodule formation, whereas repressing these targets with artificial microRNAs identical to the respective tRFs or mutating these targets with CRISPR-Cas9 promotes nodulation. Our findings thus uncover a bacterial small RNA-mediated mechanism for prokaryote-eukaryote interaction and may pave the way for enhancing nodulation efficiency in legumes.

Cell Cycle-Dependent Regulation and Function of ARGONAUTE1 in Plants

Regulated gene expression is key to the orchestrated progression of the cell cycle. Many genes are expressed at specific points in the cell cycle, including important cell cycle regulators, plus factors involved in signal transduction, hormonal regulation, and metabolic control. We demonstrate that post-embryonic depletion of Arabidopsis (Arabidopsis thaliana) ARGONAUTE1 (AGO1), the main effector of plant microRNAs (miRNAs), impairs cell division in the root meristem. We utilized the highly synchronizable tobacco (Nicotiana tabacum) Bright yellow 2 (BY2) cell suspension to analyze mRNA, small RNAs, and mRNA cleavage products of synchronized BY2 cells at S, G2, M, and G1 phases of the cell cycle. This revealed that in plants, only a few miRNAs show differential accumulation during the cell cycle, and miRNA-target pairs were only identified for a small proportion of the more than 13,000 differentially expressed genes during the cell cycle. However, this unique set of miRNA-target pairs could be key to attenuate the expression of several transcription factors and disease resistance genes. We also demonstrate that AGO1 binds to a set of 19-nucleotide, tRNA-derived fragments during the cell cycle progression.

Small RNA discovery in the interaction between barley and the powdery mildew pathogen

Background

Plants encounter pathogenic and non-pathogenic microorganisms on a nearly constant basis. Small RNAs such as siRNAs and miRNAs/milRNAs influence pathogen virulence and host defense responses. We exploited the biotrophic interaction between the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh), and its diploid host plant, barley (Hordeum vulgare) to explore fungal and plant sRNAs expressed during Bgh infection of barley leaf epidermal cells.

Results

RNA was isolated from four fast-neutron immune-signaling mutants and their progenitor over a time course representing key stages of Bgh infection, including appressorium formation, penetration of epidermal cells, and development of haustorial feeding structures. The Cereal Introduction (CI) 16151 progenitor carries the resistance allele Mla6, while Bgh isolate 5874 harbors the AVR_a6 avirulence effector, resulting in an incompatible interaction. Parallel Analysis of RNA Ends (PARE) was used to verify sRNAs with likely transcript targets in both barley and Bgh. Bgh sRNAs are predicted to regulate effectors, metabolic genes, and translation-related genes. Barley sRNAs are predicted to influence the accumulation of transcripts that encode auxin response factors, NAC transcription factors, homeodomain transcription factors, and several splicing factors. We also identified phasing small interfering RNAs (phasiRNAs) in barley that overlap transcripts that encode receptor-like kinases (RLKs) and nucleotide-binding, leucine-rich domain proteins (NLRs).

Conclusions

These data suggest that Bgh sRNAs regulate gene expression in metabolism, translation-related, and pathogen effectors. PARE-validated targets of predicted Bgh milRNAs include both EKA (effectors homologous to AVR_k1 and AVR_a10) and CSEP (candidate secreted effector protein) families. We also identified barley phasiRNAs and miRNAs in response to Bgh infection. These include phasiRNA loci that overlap with a significant proportion of receptor-like kinases, suggesting an additional sRNA control mechanism may be active in barley leaves as opposed to predominant R-gene phasiRNA overlap in many eudicots. In addition, we identified conserved miRNAs, novel miRNA candidates, and barley genome mapped sRNAs that have PARE validated transcript targets in barley. The miRNA target transcripts are enriched in transcription factors, signaling-related proteins, and photosynthesis-related proteins. Together these results suggest both barley and Bgh control metabolism and infection-related responses via the specific accumulation and targeting of genes via sRNAs.

Electronic supplementary material

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

Construction and analysis of degradome-dependent microRNA regulatory networks in soybean

BACKGROUND

Usually the microRNA (miRNA)-mediated gene regulatory network (GRN) is constructed from the investigation of miRNA expression profiling and target predictions. However, the higher/lower expression level of miRNAs does not always indicate the higher/lower level of cleavages and such analysis, thus, sometimes ignores the crucial cleavage events. In the present work, the degradome sequencing data were employed to construct the complete miRNA-mediated gene regulatory network in soybean, unlike the traditional approach starting with small RNA sequencing data.

RESULTS

We constructed the root-, cotyledon-, leaf- and seed-specific miRNA regulatory networks with the degradome sequencing data and the forthcoming verification of miRNA profiling analysis. As a result, we identified 205 conserved miRNA-target interactions (MTIs) involved with 6 conserved gma-miRNA families and 365 tissue-specific MTIs containing 24 root-specific, 45 leaf-specific, 63 cotyledon-specific and 225 seed-specific MTIs. We found a total of 156 miRNAs in tissue-specific MTIs including 18 tissue-specific miRNAs, however, only 3 miRNAs have consistent tissue-specific expression. Our study showed the degradome-dependent miRNA regulatory networks (DDNs) in four soybean tissues and explored their conservations and specificities.

CONCLUSIONS

The construction of DDNs may provide the complete miRNA-Target interactions in certain plant tissues, leading to the identification of the conserved and tissue-specific MTIs and sub-networks. Our work provides a basis for further investigation of the roles and mechanisms of miRNA-mediated regulation of tissue-specific growth and development in soybean.

Lineage-Specific Evolved MicroRNAs Regulating NB-LRR Defense Genes in Triticeae

Disease resistance genes encoding proteins with nucleotide binding sites and Leucine-Rich Repeat (NB-LRR) domains include many members involved in the effector-triggered immunity pathway in plants. The transcript levels of these defense genes are negatively regulated by diverse microRNAs (miRNAs) in angiosperms and gymnosperms. In wheat, using small RNA expression datasets and degradome datasets, we identified five miRNA families targeting NB-LRR defense genes in monocots, some of which arose in the Triticeae species era. These miRNAs regulate different types of NB-LRR genes, most of them with coil-coiled domains, and trigger the generation of secondary small interfering RNAs (siRNA) as a phased pattern in the target site regions. In addition to acting in response to biotic stresses, they are also responsive to abiotic stresses such as heat, drought, salt, and light stress. Their copy number and expression variation in Triticeae suggest a rapid birth and death frequency. Altogether, non-conserved miRNAs as conserved transcriptional regulators in gymnosperms and angiosperms regulating the disease resistance genes displayed quick plasticity including the variations of sequences, gene copy number, functions, and expression level, which accompanied with NB-LRR genes may be tune-regulated to plants in natural environments with various biotic and abiotic stresses.

MicroRNA-guided regulation of heat stress response in wheat

Background

With rising global temperature, understanding plants’ adaptation to heat stress has implications in plant breeding. MicroRNAs (miRNAs) are small, non-coding, regulatory RNAs guiding gene expression at the post-transcriptional level. In this study, small RNAs and the degradome (parallel analysis of RNA ends) of leaf tissues collected from control and heat-stressed wheat plants immediately at the end of the stress period and 1 and 4 days later were analysed.

Results

Sequencing of 24 small RNA libraries produced 55.2 M reads while 404 M reads were obtained from the corresponding 24 PARE libraries. From these, 202 miRNAs were ascertained, of which mature miRNA evidence was obtained for 104 and 36 were found to be differentially expressed after heat stress. The PARE analysis identified 589 transcripts targeted by 84 of the ascertained miRNAs. PARE sequencing validated the targets of the conserved members of miRNA156, miR166 and miR393 families as squamosa promoter-binding-like, homeobox leucine-zipper and transport inhibitor responsive proteins, respectively. Heat stress responsive miRNA targeted superoxide dismutases and an array of homeobox leucine-zipper proteins, F-box proteins and protein kinases. Query of miRNA targets to interactome databases revealed a predominant association of stress responses such as signalling, antioxidant activity and ubiquitination to superoxide dismutases, F-box proteins, pentatricopeptide repeat-containing proteins and mitochondrial transcription termination factor-like proteins.

Conclusion

The interlaced data set generated in this study identified and validated heat stress regulated miRNAs and their target genes associated with thermotolerance. Such accurate identification and validation of miRNAs and their target genes are essential to develop novel regulatory gene-based breeding strategies.

Electronic supplementary material

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

Post-transcriptional adaptation of the aquatic plant Spirodela polyrhiza under stress and hormonal stimuli

The Lemnaceae family comprises aquatic plants of angiosperms gaining attention due to their utility in wastewater treatment, and rapid production of biomass that can be used as feed, fuel, or food. Moreover, it can serve as a model species for neotenous growth and environmental adaptation. The latter properties are subject to post-transcriptional regulation of gene expression, meriting investigation of how miRNAs in Spirodela polyrhiza, the most basal and most thoroughly sequenced member of the family, are expressed under different growth conditions. To further scientific understanding of its capacity to adapt to environmental cues, we measured miRNA expression and processing of their target sequences under different temperatures, and in the presence of abscisic acid, copper, kinetin, nitrate, and sucrose. Using two small RNA sequencing experiments and one degradome sequencing experiment, we provide evidence for 108 miRNAs. Sequencing cleaved mRNAs validated 42 conserved miRNAs with 83 targets and 24 novel miRNAs regulating 66 targets and created a list of 575 predicted and verified targets. These analyses revealed condition-induced changes in miRNA expression and cleavage activity, and resulted in the addition of stringently reviewed miRNAs to miRBase. This combination of small RNA and degradome sequencing provided not only high confidence predictions of conserved and novel miRNAs and targets, but also a view of the post-transcriptional regulation of adaptations. A unique aspect is the role of miR156 and miR172 expression and activity in its clonal propagation and neoteny. Additionally, low levels of 24 nt sRNAs were observed, despite the lack of recent retrotransposition.

Stress-responsive miRNAome of Glycine max (L.) Merrill: molecular insights and way forward

MAIN CONCLUSION

Analysis of stress-associated miRNAs of Glycine max (L.) Merrill reveals wider ramifications of small RNA-mediated (conserved and legume-specific miRNAs) gene regulatory foot prints in molecular adaptive responses. MicroRNAs (miRNAs) are indispensable components of gene regulatory mechanism of plants. Soybean is a crop of immense commercial potential grown worldwide for its edible oil and soy meal. Intensive research efforts, using the next generation sequencing and bioinformatics techniques, have led to the identification and characterization of numerous small RNAs, especially microRNAs (miRNAs), in soybean. Furthermore, studies have unequivocally demonstrated the significance of miRNAs during the developmental processes and various stresses in soybean. In this review, we summarize the current state of understanding of miRNA-based abiotic and biotic stress responses in soybean. In addition, the molecular insights gained from the stress-related soybean miRNAs have been compared to the miRNAs of other crops, especially legumes, and the core commonalities have been highlighted, though differences among them were not ignored. Nature of response of soybean-derived conserved miRNAs during various stresses was also analyzed to gain deeper insights regarding sRNAome-based defense responses. This review further provides way forward in legume small RNA transcriptomics based on the adaptive responses of soybean and other legume-derived miRNAs.

Characterization of the Spatial and Temporal Expression of Two Soybean miRNAs Identifies SCL6 as a Novel Regulator of Soybean Nodulation

MicroRNAs (miRNAs) control expression of endogenous target genes through transcript cleavage or translational inhibition. Legume plants can form a specialized organ, the nodule, through interaction with nitrogen fixing soil bacteria. To understand the regulatory roles of miRNAs in the nodulation process, we functionally validated gma-miR171o and gma-miR171q and their target genes in soybean. These two miRNA sequences are identical in sequence but their miRNA genes are divergent and show unique, tissue-specific expression patterns. The expression levels of the two miRNAs are negatively correlated with that of their target genes. Ectopic expression of these miRNAs in transgenic hairy roots resulted in a significant reduction in nodule formation. Both gma-miR171o and gma-miR171q target members of the GRAS transcription factor superfamily, namely GmSCL-6 and GmNSP2. Transient interaction of miRNAs and their target genes in tobacco cells further confirmed their cleavage activity. The results suggest that gma-miR171o and gma-miR171q regulate GmSCL-6 and GmNSP2, which in turn, influence expression of the Nodule inception (NIN), Early Nodulin 40 (ENOD40), and Ethylene Response Factor Required for Nodulation (ERN) genes during the Bradyrhizobium japonicum-soybean nodulation process. Collectively, our data suggest a role for two miRNAs, gma-miR171o and gma-miR171q, in regulating the spatial and temporal aspects of soybean nodulation.