Small RNA pathways and diversity in model legumes: lessons from genomics

Genome-wide miRNAs profiles of pearl millet under contrasting high vapor pressure deficit reveal their functional roles in drought stress adaptations

Pearl millet (Pennisetum glaucum [L.] R. Br.) is an important crop capable of growing in harsh and marginal environments, with the highest degree of tolerance to drought and heat stresses among cereals. Diverse germplasm of pearl millet shows a significant phenotypic variation in response to abiotic stresses, making it a unique model to study the mechanisms responsible for stress mitigation. The present study focuses on identifying the physiological response of two pearl millet high-resolution cross (HRC) genotypes, ICMR 1122 and ICMR 1152, in response to low and high vapor pressure deficit (VPD). Under high VPD conditions, ICMR 1152 exhibited a lower transpiration rate (Tr), higher transpiration efficiency, and lower root sap exudation than ICMR 1122. Further, Pg-miRNAs expressed in the contrasting genotypes under low and high VPD conditions were identified by deep sequencing analysis. A total of 116 known and 61 novel Pg-miRNAs were identified from ICMR 1152, while 26 known and six novel Pg-miRNAs were identified from ICMR 1122 genotypes, respectively. While Pg-miR165, 168, 170, and 319 families exhibited significant differential expression under low and high VPD conditions in both genotypes, ICMR 1152 showed abundant expression of Pg-miR167, Pg-miR172, Pg-miR396 Pg-miR399, Pg-miR862, Pg-miR868, Pg-miR950, Pg-miR5054, and Pg-miR7527 indicating their direct and indirect role in root physiology and abiotic stress responses. Drought responsive Pg-miRNA targets showed upregulation in response to high VPD stress, further narrowing down the miRNAs involved in regulation of drought tolerance in pearl millet.

Genome-Wide Identification of Key Components of RNA Silencing in Two Phaseolus vulgaris Genotypes of Contrasting Origin and Their Expression Analyses in Response to Fungal Infection

RNA silencing serves key roles in a multitude of cellular processes, including development, stress responses, metabolism, and maintenance of genome integrity. Dicer, Argonaute (AGO), double-stranded RNA binding (DRB) proteins, RNA-dependent RNA polymerase (RDR), and DNA-dependent RNA polymerases known as Pol IV and Pol V form core components to trigger RNA silencing. Common bean (Phaseolus vulgaris) is an important staple crop worldwide. In this study, we aimed to unravel the components of the RNA-guided silencing pathway in this non-model plant, taking advantage of the availability of two genome assemblies of Andean and Meso-American origin. We identified six PvDCLs, thirteen PvAGOs, 10 PvDRBs, 5 PvRDRs, in both genotypes, suggesting no recent gene amplification or deletion after the gene pool separation. In addition, we identified one PvNRPD1 and one PvNRPE1 encoding the largest subunits of Pol IV and Pol V, respectively. These genes were categorized into subgroups based on phylogenetic analyses. Comprehensive analyses of gene structure, genomic localization, and similarity among these genes were performed. Their expression patterns were investigated by means of expression models in different organs using online data and quantitative RT-PCR after pathogen infection. Several of the candidate genes were up-regulated after infection with the fungus Colletotrichum lindemuthianum.

LuluDB-The Database Created Based on Small RNA, Transcriptome, and Degradome Sequencing Shows the Wide Landscape of Non-coding and Coding RNA in Yellow Lupine (Lupinus luteus L.) Flowers and Pods

Yellow lupine (Lupinus luteus L.) belongs to a legume family that benefits from symbiosis with nitrogen-fixing bacteria. Its seeds are rich in protein, which makes it a valuable food source for animals and humans. Yellow lupine is also the model plant for basic research on nodulation or abscission of organs. Nevertheless, the knowledge about the molecular regulatory mechanisms of its generative development is still incomplete. The RNA-Seq technique is becoming more prominent in high-throughput identification and expression profiling of both coding and non-coding RNA sequences. However, the huge amount of data generated with this method may discourage other scientific groups from making full use of them. To overcome this inconvenience, we have created a database containing analysis-ready information about non-coding and coding L. luteus RNA sequences (LuluDB). LuluDB was created on the basis of RNA-Seq analysis of small RNA, transcriptome, and degradome libraries obtained from yellow lupine cv. Taper flowers, pod walls, and seeds in various stages of development, flower pedicels, and pods undergoing abscission or maintained on the plant. It contains sequences of miRNAs and phased siRNAs identified in L. luteus, information about their expression in individual samples, and their target sequences. LuluDB also contains identified lncRNAs and protein-coding RNA sequences with their organ expression and annotations to widely used databases like GO, KEGG, NCBI, Rfam, Pfam, etc. The database also provides sequence homology search by BLAST using, e.g., an unknown sequence as a query. To present the full capabilities offered by our database, we performed a case study concerning transcripts annotated as DCL 1–4 (DICER LIKE 1–4) homologs involved in small non-coding RNA biogenesis and identified miRNAs that most likely regulate DCL1 and DCL2 expression in yellow lupine. LuluDB is available at http://luluseqdb.umk.pl/basic/web/index.php.

Argonaute Proteins: Why Are They So Important for the Legume-Rhizobia Symbiosis?

Unlike most other land plants, legumes can fulfill their nitrogen needs through the establishment of symbioses with nitrogen-fixing soil bacteria (rhizobia). Through this symbiosis, fixed nitrogen is incorporated into the food chain. Because of this ecological relevance, the genetic mechanisms underlying the establishment of the legume-rhizobia symbiosis (LRS) have been extensively studied over the past decades. During this time, different types of regulators of this symbiosis have been discovered and characterized. A growing number of studies have demonstrated the participation of different types of small RNAs, including microRNAs, in the different stages of this symbiosis. The involvement of small RNAs also indicates that Argonaute (AGO) proteins participate in the regulation of the LRS. However, despite this obvious role, the relevance of AGO proteins in the LRS has been overlooked and understudied. Here, we discuss and hypothesize the likely participation of AGO proteins in the regulation of the different steps that enable the establishment of the LRS. We also briefly review and discuss whether rhizobial symbiosis induces DNA damages in the legume host. Understanding the different levels of LRS regulation could lead to the development of improved nitrogen fixation efficiency to enhance sustainable agriculture, thereby reducing dependence on inorganic fertilizers.

Enhancement of ABA Sensitivity Through Conditional Expression of the ARF10 Gene in Brassica juncea Reveals Fertile Plants with Tolerance Against Alternaria brassicicola

Concomitant increase of auxin-responsive factors ARF16 and ARF17, along with enhanced expression of ARF10 in resistant Sinapis alba compared with that in susceptible Brassica juncea upon challenge with Alternaria brassicicola, revealed that abscisic acid (ABA)-auxin crosstalk is a critical factor for resistance response. Here, we induced the ABA response through conditional expression of ARF10 in B. juncea using the A. brassicicola-inducible GH3.3 promoter. Induced ABA sensitivity caused by conditional expression of ARF10 in transgenic B. juncea resulted in tolerance against A. brassicicola and led to enhanced expression of several ABA-responsive genes without affecting the auxin biosynthetic gene expression. Compared with ABI3 and ABI4, ABI5 showed maximum upregulation in the most tolerant transgenic lines upon pathogen challenge. Moreover, elevated expression of ARF10 by different means revealed a direct correlation between ARF10 expression and the induction of ABI5 protein in B. juncea. Through in vitro DNA-protein experiments and chromosome immunoprecipitation using the ARF10 antibody, we demonstrated that ARF10 interacts with the auxin-responsive elements of the ABI5 promoter. This suggests that ARF10 may function as a modulator of ABI5 to induce ABA sensitivity and mediate the resistance response against A. brassicicola.

Identification and profiling of narrow-leafed lupin (Lupinus angustifolius) microRNAs during seed development

BACKGROUND

Whilst information regarding small RNAs within agricultural crops is increasing, the miRNA composition of the nutritionally valuable pulse narrow-leafed lupin (Lupinus angustifolius) remains unknown.

RESULTS

By conducting a genome- and transcriptome-wide survey we identified 7 Dicer-like and 16 Argonaute narrow-leafed lupin genes, which were highly homologous to their legume counterparts. We identified 43 conserved miRNAs belonging to 16 families, and 13 novel narrow-leafed lupin-specific miRNAs using high-throughput sequencing of small RNAs from foliar and root and five seed development stages. We observed up-regulation of members of the miRNA families miR167, miR399, miR156, miR319 and miR164 in narrow-leafed lupin seeds, and confirmed expression of miR156, miR166, miR164, miR1507 and miR396 using quantitative RT-PCR during five narrow-leafed lupin seed development stages. We identified potential targets for the conserved and novel miRNAs and were able to validate targets of miR399 and miR159 using 5' RLM-RACE. The conserved miRNAs are predicted to predominately target transcription factors and 93% of the conserved miRNAs originate from intergenic regions. In contrast, only 43% of the novel miRNAs originate from intergenic regions and their predicted targets were more functionally diverse.

CONCLUSION

This study provides important insights into the miRNA gene regulatory networks during narrow-leafed lupin seed development.

Genome-wide identification and characterization of novel microRNAs in seed development of soybean

MicroRNAs (miRNAs) are important and ubiquitous regulators of gene expression in eukaryotes. However, the information about miRNAs population and their regulatory functions involving in soybean seed development remains incomplete. Base on the Dicer-like1-mediated cleavage signals during miRNA processing could be employed for novel miRNA discovery, a genome-wide search for miRNA candidates involved in seed development was carried out. As a result, 17 novel miRNAs, 14 isoforms of miRNA (isomiRs) and 31 previously validated miRNAs were discovered. These novel miRNAs and isomiRs represented tissue-specific expression and the isomiRs showed significantly higher abundance than that of their miRNA counterparts in different tissues. After target prediction and degradome sequencing data-based validation, 13 novel miRNA-target pairs were further identified. Besides, five targets of 22-nt iso-gma-miR393h were found to be triggered to produce secondary trans-acting siRNA (ta-siRNAs). Summarily, our results could expand the repertoire of miRNAs with potentially important functions in soybean.

Legume genetic resources and transcriptome dynamics under abiotic stress conditions

Grain legumes are an important source of nutrition and income for billions of consumers and farmers around the world. However, the low productivity of new legume varieties, due to the limited genetic diversity available for legume breeding programmes and poor policymaker support, combined with an increasingly unpredictable global climate is resulting in a large gap between current yields and the increasing demand for legumes as food. Hence, there is a need for novel approaches to develop new high-yielding legume cultivars that are able to cope with a range of environmental stressors. Next-generation technologies are providing the tools that could enable the more rapid and cost-effective genomic and transcriptomic studies for most major crops, allowing the identification of key functional and regulatory genes involved in abiotic stress resistance. In this review, we provide an overview of the recent achievements regarding abiotic stress resistance in a wide range of legume crops and highlight the transcriptomic and miRNA approaches that have been used. In addition, we critically evaluate the availability and importance of legume genetic resources with desirable abiotic stress resistance traits.

CRISPR/Cas9 and TALENs generate heritable mutations for genes involved in small RNA processing of Glycine max and Medicago truncatula

Processing of double-stranded RNA precursors into small RNAs is an essential regulator of gene expression in plant development and stress response. Small RNA processing requires the combined activity of a functionally diverse group of molecular components. However, in most of the plant species, there are insufficient mutant resources to functionally characterize each encoding gene. Here, mutations in loci encoding protein machinery involved in small RNA processing in soya bean and Medicago truncatula were generated using the CRISPR/Cas9 and TAL-effector nuclease (TALEN) mutagenesis platforms. An efficient CRISPR/Cas9 reagent was used to create a bi-allelic double mutant for the two soya bean paralogous Double-stranded RNA-binding2 (GmDrb2a and GmDrb2b) genes. These mutations, along with a CRISPR/Cas9-generated mutation of the M. truncatula Hua enhancer1 (MtHen1) gene, were determined to be germ-line transmissible. Furthermore, TALENs were used to generate a mutation within the soya bean Dicer-like2 gene. CRISPR/Cas9 mutagenesis of the soya bean Dicer-like3 gene and the GmHen1a gene was observed in the T0 generation, but these mutations failed to transmit to the T1 generation. The irregular transmission of induced mutations and the corresponding transgenes was investigated by whole-genome sequencing to reveal a spectrum of non-germ-line-targeted mutations and multiple transgene insertion events. Finally, a suite of combinatorial mutant plants were generated by combining the previously reported Gmdcl1a, Gmdcl1b and Gmdcl4b mutants with the Gmdrb2ab double mutant. Altogether, this study demonstrates the synergistic use of different genome engineering platforms to generate a collection of useful mutant plant lines for future study of small RNA processing in legume crops.

Functional Genomics and Seed Development in Medicago truncatula: An Overview

The study of seed development in the model species Medicago truncatula has made a significant contribution to our understanding of this process in crop legumes. Thanks to the availability of comprehensive proteomics and transcriptomics databases, coupled with exhaustive mutant collections, the roles of several regulatory genes in development and maturation are beginning to be deciphered and functionally validated. Advances in next-generation sequencing and the availability of a genomic sequence have made feasible high-density SNP genotyping, allowing the identification of markers tightly linked to traits of agronomic interest. A further major advance is to be expected from the integration of omics resources in functional network construction, which has been used recently to identify "hub" genes central to important traits.

Regulation of Small RNAs and Corresponding Targets in Nod Factor-Induced Phaseolus vulgaris Root Hair Cells

A genome-wide analysis identified the set of small RNAs (sRNAs) from the agronomical important legume Phaseolus vulgaris (common bean), including novel P. vulgaris-specific microRNAs (miRNAs) potentially important for the regulation of the rhizobia-symbiotic process. Generally, novel miRNAs are difficult to identify and study because they are very lowly expressed in a tissue- or cell-specific manner. In this work, we aimed to analyze sRNAs from common bean root hairs (RH), a single-cell model, induced with pure Rhizobium etli nodulation factors (NF), a unique type of signal molecule. The sequence analysis of samples from NF-induced and control libraries led to the identity of 132 mature miRNAs, including 63 novel miRNAs and 1984 phasiRNAs. From these, six miRNAs were significantly differentially expressed during NF induction, including one novel miRNA: miR-RH82. A parallel degradome analysis of the same samples revealed 29 targets potentially cleaved by novel miRNAs specifically in NF-induced RH samples; however, these novel miRNAs were not differentially accumulated in this tissue. This study reveals Phaseolus vulgaris-specific novel miRNA candidates and their corresponding targets that meet all criteria to be involved in the regulation of the early nodulation events, thus setting the basis for exploring miRNA-mediated improvement of the common bean–rhizobia symbiosis.

Six Medicago truncatula Dicer-like protein genes are expressed in plant cells and upregulated in nodules

Here we report the existence of six putative Dicer-like genes in the Medicago truncatula genome. They are ubiquitously expressed throughout the plant and significantly induced in root nodules. Over the past decade, small noncoding RNAs (sncRNA) have emerged as widespread and important regulatory molecules influencing both the structure and expression of plant genomes. One of the key factors involved in sncRNA biogenesis in plants is a group of RNase III-type nucleases known as Dicer-like (DCL) proteins. Based on functional analysis of DCL proteins identified in Arabidopsis thaliana, four types of DCLs were distinguished (DCL1-4). DCL1 mainly produces 21 nt miRNAs. The products generated by DCL2, DCL3, and DCL4 belong to various classes of siRNAs that are 22, 24 and 21 nt in length, respectively. M. truncatula is a model legume plant closely related to many economically important cultivable species. By screening the recent M. truncatula genome assembly, we were able to identify three new DCL genes in addition to the MtDCL1-3 genes that had been earlier characterized. The newly found genes include MtDCL4 and two MtDCL2 homologs. We showed that all six M. truncatula DCL genes are expressed in plant cells. The first of the identified MtDCL2 paralogs encodes a truncated version of the DCL2 protein, while the second undergoes substantial and specific upregulation in the root nodules. Additionally, we identified an alternative splicing variant of MtDCL1 mRNA, similar to the one found in Arabidopsis. Our results indicate that DCL genes are differently activated during Medicago symbiosis with nitrogen fixing bacteria and upon pathogen infection. In addition, we hypothesize that the alternative splicing variant of MtDCL1 mRNA may be involved in tissue-specific regulation of the DCL1 level.

Genome-wide identification of the Phaseolus vulgaris sRNAome using small RNA and degradome sequencing

BACKGROUND

MiRNAs and phasiRNAs are negative regulators of gene expression. These small RNAs have been extensively studied in plant model species but only 10 mature microRNAs are present in miRBase version 21, the most used miRNA database, and no phasiRNAs have been identified for the model legume Phaseolus vulgaris. Thanks to the recent availability of the first version of the common bean genome, degradome data and small RNA libraries, we are able to present here a catalog of the microRNAs and phasiRNAs for this organism and, particularly, we suggest new protagonists in the symbiotic nodulation events.

RESULTS

We identified a set of 185 mature miRNAs, including 121 previously unpublished sequences, encoded by 307 precursors and distributed in 98 families. Degradome data allowed us to identify a total of 181 targets for these miRNAs. We reveal two regulatory networks involving conserved miRNAs: those known to play crucial roles in the establishment of nodules, and novel miRNAs present only in common bean, suggesting a specific role for these sequences. In addition, we identified 125 loci that potentially produce phased small RNAs, with 47 of them having all the characteristics of being triggered by a total of 31 miRNAs, including 14 new miRNAs identified in this study.

CONCLUSIONS

We provide here a set of new small RNAs that contribute to the broader knowledge of the sRNAome of Phaseolus vulgaris. Thanks to the identification of the miRNA targets from degradome analysis and the construction of regulatory networks between the mature microRNAs, we present here the probable functional regulation associated with the sRNAome and, particularly, in N2-fixing symbiotic nodules.

High throughput sequencing of small RNA component of leaves and inflorescence revealed conserved and novel miRNAs as well as phasiRNA loci in chickpea

Among legumes, chickpea (Cicer arietinum L.) is the second most important crop after soybean. MicroRNAs (miRNAs) play important roles by regulating target gene expression important for plant development and tolerance to stress conditions. Additionally, recently discovered phased siRNAs (phasiRNAs), a new class of small RNAs, are abundantly produced in legumes. Nevertheless, little is known about these regulatory molecules in chickpea. The small RNA population was sequenced from leaves and flowers of chickpea to identify conserved and novel miRNAs as well as phasiRNAs/phasiRNA loci. Bioinformatics analysis revealed 157 miRNA loci for the 96 highly conserved and known miRNA homologs belonging to 38 miRNA families in chickpea. Furthermore, 20 novel miRNAs belonging to 17 miRNA families were identified. Sequence analysis revealed approximately 60 phasiRNA loci. Potential target genes likely to be regulated by these miRNAs were predicted and some were confirmed by modified 5' RACE assay. Predicted targets are mostly transcription factors that might be important for developmental processes, and others include superoxide dismutases, plantacyanin, laccases and F-box proteins that could participate in stress responses and protein degradation. Overall, this study provides an inventory of miRNA-target gene interactions for chickpea, useful for the comparative analysis of small RNAs among legumes.

A stress-induced small RNA modulates alpha-rhizobial cell cycle progression

Mechanisms adjusting replication initiation and cell cycle progression in response to environmental conditions are crucial for microbial survival. Functional characterization of the trans-encoded small non-coding RNA (trans-sRNA) EcpR1 in the plant-symbiotic alpha-proteobacterium Sinorhizobium meliloti revealed a role of this class of riboregulators in modulation of cell cycle regulation. EcpR1 is broadly conserved in at least five families of the Rhizobiales and is predicted to form a stable structure with two defined stem-loop domains. In S. meliloti, this trans-sRNA is encoded downstream of the divK-pleD operon. ecpR1 belongs to the stringent response regulon, and its expression was induced by various stress factors and in stationary phase. Induced EcpR1 overproduction led to cell elongation and increased DNA content, while deletion of ecpR1 resulted in reduced competitiveness. Computationally predicted EcpR1 targets were enriched with cell cycle-related mRNAs. Post-transcriptional repression of the cell cycle key regulatory genes gcrA and dnaA mediated by mRNA base-pairing with the strongly conserved loop 1 of EcpR1 was experimentally confirmed by two-plasmid differential gene expression assays and compensatory changes in sRNA and mRNA. Evidence is presented for EcpR1 promoting RNase E-dependent degradation of the dnaA mRNA. We propose that EcpR1 contributes to modulation of cell cycle regulation under detrimental conditions.

Microorganisms frequently encounter adverse conditions unfavorable for cell proliferation. They have evolved diverse mechanisms, including transcriptional control and targeted protein degradation, to adjust cell cycle progression in response to environmental cues. Non-coding RNAs are widespread regulators of various cellular processes in all domains of life. In prokaryotes, trans-encoded small non-coding RNAs (trans-sRNAs) contribute to a rapid cellular response to changing environments, but so far have not been directly related to cell cycle regulation. Here, we report the first example of a trans-sRNA (EcpR1) with two experimentally confirmed targets in the core of cell cycle regulation and demonstrate that in the plant-symbiotic alpha-proteobacterium Sinorhizobium meliloti the regulatory mechanism involves base-pairing of this sRNA with the dnaA and gcrA mRNAs. Most trans-sRNAs are restricted to closely related species, but the stress-induced EcpR1 is broadly conserved in the order of Rhizobiales suggesting an evolutionary advantage conferred by ecpR1. It broadens the functional diversity of prokaryotic sRNAs and adds a new regulatory level to the mechanisms that contribute to interlinking stress responses with the cell cycle machinery.

The small RNA diversity from Medicago truncatula roots under biotic interactions evidences the environmental plasticity of the miRNAome

BACKGROUND

Legume roots show a remarkable plasticity to adapt their architecture to biotic and abiotic constraints, including symbiotic interactions. However, global analysis of miRNA regulation in roots is limited, and a global view of the evolution of miRNA-mediated diversification in different ecotypes is lacking.

RESULTS

In the model legume Medicago truncatula, we analyze the small RNA transcriptome of roots submitted to symbiotic and pathogenic interactions. Genome mapping and a computational pipeline identify 416 miRNA candidates, including known and novel variants of 78 miRNA families present in miRBase. Stringent criteria of pre-miRNA prediction yield 52 new mtr-miRNAs, including 27 miRtrons. Analyzing miRNA precursor polymorphisms in 26 M. truncatula ecotypes identifies higher sequence polymorphism in conserved rather than Medicago-specific miRNA precursors. An average of 19 targets, mainly involved in environmental responses and signalling, is predicted per novel miRNA. We identify miRNAs responsive to bacterial and fungal pathogens or symbionts as well as their related Nod and Myc-LCO symbiotic signals. Network analyses reveal modules of new and conserved co-expressed miRNAs that regulate distinct sets of targets, highlighting potential miRNA-regulated biological pathways relevant to pathogenic and symbiotic interactions.

CONCLUSIONS

We identify 52 novel genuine miRNAs and large plasticity of the root miRNAome in response to the environment, and also in response to purified Myc/Nod signaling molecules. The new miRNAs identified and their sequence variation across M. truncatula ecotypes may be crucial to understand the adaptation of root growth to the soil environment, notably in the agriculturally important legume crops.