Skip to main content

Genotype- and tissue-specific miRNA profiles and their targets in three alfalfa (Medicago sativa L) genotypes



Alfalfa (Medicago sativa L.) is a forage legume with significant agricultural value worldwide. MicroRNAs (miRNAs) are key components of post-transcriptional gene regulation and essentially regulate many aspects of plant growth and development. Although miRNAs were reported in alfalfa, their expression profiles in different tissues and the discovery of novel miRNAs as well as their targets have not been described in this plant species.


To identify tissue-specific miRNA profiles in whole plants, shoots and roots of three different alfalfa genotypes (Altet-4, NECS-141and NF08ALF06) were used. Small RNA libraries were generated and sequenced using a high-throughput sequencing platform. Analysis of these libraries enabled identification of100 miRNA families; 21 of them belong to the highly conserved families while the remaining 79 families are conserved at the minimum between M. sativa and the model legume and close relative, M. truncatula. The profiles of the six abundantly expressed miRNA families (miR156, miR159, miR166, miR319, miR396 and miR398) were relatively similar between the whole plants, roots and shoots of these three alfalfa genotypes. In contrast, robust differences between shoots and roots for miR160 and miR408 levels were evident, and their expression was more abundant in the shoots. Additionally, 17 novel miRNAs were identified and the relative abundance of some of these differed between tissue types. Further, the generation and analysis of degradome libraries from the three alfalfa genotypes enabled confirmation of 69 genes as targets for 31 miRNA families in alfalfa.


The miRNA profiles revealed both similarities and differences in the expression profiles between tissues within a genotype as well as between the genotypes. Among the highly conserved miRNA families, miR166 was the most abundantly expressed in almost all tissues from the three genotypes. The identification of conserved and novel miRNAs as well as their targets in different tissues of multiple genotypes increased our understanding of miRNA-mediated gene regulation in alfalfa and could provide valuable insights for practical research and plant improvement applications in alfalfa and related legume species.


Alfalfa (Medicago sativa L.) is an important forage legume species with global adaptation, high forage quality and the capacity for harvesting biomass multiple times during the growing season. Alfalfa is an autotetraploid (2n = 4x = 32), perennial outcrossing species with high levels of genetic diversity in cultivated and non-cultivated populations. Besides its use as a forage, alfalfa also has potential crop for biofuel production [1]. Alfalfa has the capacity for symbiotic nitrogen fixation and can also contribute to reduce soil erosion [2, 3].

Endogenous non-coding RNAs of approximately 21–22 nucleotides represent plant miRNAs that silence gene expression by binding to complementary sequences of its target mRNA at the post-transcriptional level. Such targeting results in mRNA cleavage and degradation or repression of translation, with the former being more prevalent in plants [4,5,6,7]. The miRNA analyses in different plant species highlight the important regulatory roles of miRNAs in multiple organs (roots, stems, leaves and flowers), differentiation and development, leaf polarity, transition from juvenile to vegetative stages and vegetative to reproductive phases, and regulation of plant responses to biotic and abiotic stresses [8,9,10].

Several investigations have shown that plant miRNAs can be classified into conserved and novel lineage- or species-specific miRNAs. Conserved miRNAs and their corresponding target genes are commonly found in all or most angiosperms, with some also being described in gymnosperms as well as primitive land plants such as ferns [11, 12]. However, miRNA analysis in several legumes including M. truncatula, soybean (Glycine max L), chickpea (Cicer arietinum L.), common bean (Phaseolus vulgaris), and Lotus japonicus indicate the presence of miRNAs that seem to be specific to certain legumes and there could have important gene regulatory roles [13,14,15,16,17,18,19]. Although recent attempts were made to report miRNAs from alfalfa (M. sativa) [20,21,22], these do not include the discovery of novel miRNAs, and most importantly, the miRNA target genes have not been confirmed in this legume species. Understanding miRNAs and their target gene regulation in various tissues can provide further insights into the miRNA target networks operating in a tissue-specific manner in alfalfa.

In order to identify conserved miRNAs as well as novel miRNAs from alfalfa, we constructed and sequenced small RNA libraries from whole clonally propagated plants, roots and shoots of three alfalfa genotypes (Altet-4, NECS-141 and NF08ALF06). The sequenced reads were mapped to known miRNAs in M. truncatula, deposited in the miRBase to identify and annotate the miRNAs in alfalfa. Degradome libraries were constructed and sequenced from these three genotypes to characterize the miRNA gene targets.

Materials and methods

Plant materials and growth conditions

Three alfalfa genotypes NECS-141, Altet-4 and NF08ALF06 were evaluated in this study. NECS-141 is the genotype being used to sequence the tetraploid alfalfa genome [23]. Altet-4 is an aluminum tolerant genotype used to develop a mapping population [24]. NF08ALF06 is a commercial breeding line with good agronomic performance (Forage Genetics International). The three alfalfa genotypes (NECS-141, Altet-4 and NF08ALF06) were clonally propagated and grown in tissue culture. After 13 d of growth in rooting media, these were transferred to medium at pH 7 for 96 h as previously described [25]. The rooting media contains 0.55 g/L Murashige & Skoog Basal Medium with Vitamins (PhytoTechnology #M519), 1 ml Plant Preservative Mixture, PPM (PhytoTechnology), adjust the pH to 5.8, and add 12 g/L Gelzan. The plants were placed in a Conviron growth chamber (24 °C, 18 h /6 h day/night cycle, 100 μmol light intensity) for root development and growth. An additional 20 clonally propagated plants of these genotypes were grown in a Conviron growth chamber as previously described and used to evaluate the tissue-specific expression of the miRNAs. Tissue samples were harvested and immediately flash frozen in liquid nitrogen and stored at − 80 °C.

Small RNA library construction and sequencing

Total RNA was isolated from the whole plants, roots and shoots of three alfalfa genotypes using TRIzol ® Reagent (Invitrogen), according to the manufacturer’s instructions. The quality of total RNA was monitored on 1% agarose gel and their concentrations were measured using Nanodrop spectrophotometer. Small RNA libraries were generated as described previously [26] by following the protocol described for the Illumina Truseq® Small RNA Preparation kit (Illumina, San Diego, USA). Briefly, 1 μg total RNA per sample was ligated sequentially with 5′ and 3’ RNA adaptors. The ligated products were converted into cDNAs and then amplified using PCR. The amplified products were sequenced using an Illumina Hiseq® Analyzer.

Identification of conserved and novel miRNAs

The raw sequencing reads were processed as follows: adaptor sequences were trimmed off from the raw reads to obtain small RNAs. These reads were then mapped to ribosomal RNA (rRNA), transfer RNA (tRNA), small nuclear RNAs (snRNA), and the aligned and mapped reads were not used for further analysis. The remaining reads were aligned to miRBase v 20 [27] to identify miRNAs in M. sativa. The reads with 100% sequence identity were designated as conserved miRNA homologs. To identify novel miRNAs, the presence of the miRNA-star (miRNA*) sequences coupled with the predictable hairpin-like structure for the precursor sequences were used.

Degradome library construction and analyses

Degradome libraries from the alfalfa genotypes NECS-141, Altet-4 and NF08ALF06 were constructed as previously described to identify potential target mRNAs [28]. Briefly, the cleaved 5′ monophosphate containing polyadenylated mRNA fragments were ligated to an RNA oligo-nucleotide adapter containing MmeI recognition site at its 3′ end. The ligated products were converted into cDNA using reverse transcriptase and the product was amplified using only 5 PCR cycles. The PCR product was eluted, digested with MmeI restriction enzyme and then ligated to a double-stranded DNA adapter. The ligated product was again purified and amplified using 15 cycles of PCR. The final PCR product was sequenced. The reads were processed for quality and then aligned to the transcriptome assembly of M. truncatula to identify potential miRNA targets using the SeqTar pipeline [29].

Results and discussion

The analyses of small RNA libraries

High-throughput sequencing has been used to identify miRNAs and their target mRNAs in plants [15, 30, 31]. To catalogue conserved and novel miRNAs in alfalfa, a total of eight small RNA libraries from the whole plants, roots and shoots of Altet-4, NECS-141 and NF08ALF06 genotypes were constructed and sequenced. After removal of the adapter sequences and low-quality reads, the total reads ranging between 11 to 42 million, and unique reads ranging between 1.8 to 8.5 million reads from these nine libraries were obtained (Table 1). However, the quality of the small RNA library generated from the shoots of NF08ALF06 did not meet the threshold criteria, therefore only NECS-141 and Altet-4 were used for the miRNA analyses of shoot tissues.

Table 1 The mapping of total and unique reads obtained from different small RNA libraries

Quantification of miRNA abundances between the genotypes and tissues was preceded by normalizing the expression levels of miRNA families to reads per ten million (RPTM). The normalized miRNA family read frequencies ranged between 1 to 552,267 RPTM for the whole plants, between 1 to 134,679 RPTM for the root samples, and 1 to 165,310 RPTM for the shoot samples (Table 2). The range of miRNA read frequencies varied slightly between the three genotypes. As expected, the most conserved miRNAs appeared to be the most abundantly expressed in all tissues and genotypes, with the exception of miR169, miR393, miR395 and miR172 which exhibited low abundances. Specifically, miR172 levels in roots and shoots of the three genotypes were extremely low and in most cases was below 20 RPTM (Table 2). The miRNA families with the lowest expression levels, and in some cases as low as 1 RPTM, were largely represented by the non-conserved miRNAs or miRNAs that have been reported exclusively from M. truncatula (miRBase) that include miR2601, miR2674, miR5207, miR5241, miR5243, miR5244, miR5255, miR5257, miR5269, miR5282, miR5289, miR5294, miR5296, miR5299, miR5561, miR5744, and miR7701 (Table 2). miR5207 is the only miRNA that was also reported from Gossypium raimondii (miRBase). The majority of the miRNA families identified are 21 nt long, although some cases including miR2601 and miR2603 were represented by 22 nucleotides. Further, a total of 23 miRNA families included between miR5267 to miR5299 were 24 nt long. The fact that these small RNAs were initially identified in M. truncatula (miRBase), and could be identified in several independent small RNA libraries from three different alfalfa genotypes (Table 2), suggests that these sequences and their associated processing are conserved between alfalfa and its close relative M. truncatula. However, their extremely low abundances coupled with their longer read lengths could also indicate that these may be 24-nt long siRNAs. Additional studies are needed to assess the precise nature of these small RNAs, i.e., miRNAs or siRNAs.

Table 2 Identified miRNA families and their frequencies (reads per ten million [RPTM]) in whole plants, roots and shoots of three alfalfa genotypes (miRNA-stars were marked in bold)

MicroRNA profiles in alfalfa plants, roots and shoots

A total of 100 known miRNA families were identified from the small RNA libraries of the three alfalfa genotypes (Table 2). Of these, 21 families were represented by the highly conserved miRNAs, whereas the remaining 79 families could be considered as Medicago-specific miRNA families. The identification of these 79 miRNA families in alfalfa was based on their expression in M. truncatula (miRbase), therefore, these are conserved at least between M. truncatula and alfalfa.

Among the highly conserved miRNA families, miR166 was the most highly expressed family in seven of the eight samples that were surveyed in this study. The only exception to this trend was observed in the shoots of NECS-141 in which the miR2118 family was the most abundant followed by the miR166 family. The miRNA families, miR396 and miR2118 represents the second and third most abundantly expressed in the whole plants, while miR159 and miR396 were the second and third most highly expressed miRNAs in roots. Several additional miRNA families including miR398, miR160, miR168, miR319, miR408, miR1510 and miR2643 were also highly expressed but miR169, miR171, miR393, miR397 and miR395 were expressed at relatively very low levels (Table 2). On the other hand, miR159, miR156, miR319, miR398 miR1507 and miR1510 were highly expressed but miR164, miR169, miR172, miR393, miR397, miR399 and miR482 were expressed at very low levels in roots of these genotypes. Interestingly, miR160 was not sequenced from the roots of three alfalfa genotypes.

Overall, the conserved miRNA families such as the miR156, miR159, miR166, miR168, miR319, miR396, miR398 and miR408 were more highly expressed in the plants, roots and shoots of all three alfalfa genotypes. Among the legume-specific families, miR1507, miR1510, miR2118, miR2592, miR2643, miR5213, miR5232, miR5558 and miR7696 (Table 2) were also abundant in all tissues of alfalfa genotypes. Conversely, some conserved miRNA families represented by miR169 and miR393 recorded very low abundances in all samples. Other notable differences between roots and shoots include relatively low expression levels of miR160, miR167, and miR408 in roots compared to the shoots of alfalfa genotypes (Table 2).

Several miRNA families including miR482, miR1507, miR2118, miR4416 are conserved in M. truncatula, soybean, chickpea (miRBase). These miRNA families are known to regulate NBS-LRR genes that are involved in pathogen resistance. The miRNA-guided cleavage on the NBS-LRR genes initiates the generation of phasiRNAs [16, 18, 32]. In alfalfa, miR482, miR1507 and miR2118 were detected in all three tissues (Table 2), but not miR4416. Both miR2118 and miR1507 families were more abundantly expressed in all tissues and genotypes compared with miR482 family. Remarkably, miR2118 was the top most highly expressed miRNA family in shoots of NECS-141. By contrast, miR2118 levels were very low in roots of three alfalfa genotypes. On the other hand, miR1507 family displayed approximately similar levels in three tissues of alfalfa genotypes.

The miRNA-star sequences corresponding to the 12 of the 21 highly conserved miRNA families were also recovered from almost all libraries (Table 2). Additionally, miRNA-stars for the miR1510, miR4414, miR5208, and miR7696 were also detected. Furthermore, the miRNA-star expression levels for miR156, miR166 and miR168 were very high (Table 2). Intriguingly, like miR168, miR168 star levels differed greatly between different tissue. In shoots of NECS-141, miR168 star levels were slightly more than that of miR168, while both in whole plants and roots, the star levels were approximately half of the levels of miR168.

miRNA diversity in alfalfa compared with other legumes

Several miRNA families are specifically reported from the leguminous plants such as the M. truncatula, Glycine max, Lotus japonicus, Phaseolus vulgaris, Cicer arietinum, Vigna unguiculata and Acacia auriculiformis [14, 16, 18, 19, 32, 33]. These lineage-specific miRNAs include miR1507, miR1508, miR1509, miR1510, miR1512, miR1514, miR1520, miR1521, miR2118, miR2086, miR2109, miR2199, miR4414, miR5213, miR5232, and miR5234 among others (miRBase). The majority of these were reported from M. truncatula and soybean, since these legume species have been the subject of multiple studies exploring small RNAs. Most of these legume-specific miRNAs were also identified in alfalfa and a few of them including miR1507, miR1510, miR2118, miR2592, miR2643, miR5211, miR5213, miR5214, miR5232, miR5239, miR5277, miR5558, and miR7696 were specifically highly expressed in all three genotypes (Table 2).

Identification of novel miRNAs from alfalfa

The sequencing of the small RNAs from multiple tissues of three different alfalfa genotypes would allow us to identify the novel miRNAs more confidently. Novel miRNA identification was dependent on sequencing of the miRNA complementary strand (miRNA-star) coupled with the predictable fold back structure for the primary miRNA transcript. Because a stable assembly of the tetraploid alfalfa genome sequence is not available, the small RNAs were mapped to the M. truncatula genome. Mapping of the small RNAs from the three alfalfa genotypes onto the M. truncatula genome enabled the identification of novel miRNAs more confidently because they have been sequenced from M. sativa and mapped on to the M. truncatula, suggesting their conservation between M. sativa and M. truncatula. Moreover, the novel miRNA identification in this study is more robust as it includes sequencing of these small RNAs from three different genotypes. We have identified a total of 17 novel miRNAs which have been sequenced from all of the three genotypes (Table 3 and Fig. 1). Among these, t50582913 was highly expressed followed by t50063038. In roots, t50582913 was highly expressed in NECS-141 and Altet-4 but not in NF08ALF06. In shoots, t50063038 was highly expressed followed by the t50582913 and t51235783.

Table 3 Identified novel miRNAs based on sequencing both 5′ and 3′ reads and the most abundant ones that is marked in bold denotes potential novel miRNA based on their greater abundances
Fig. 1

The predicted foldback structures using the novel miRNA precursor sequences. a The fold-back structures for six novel miRNAs. b The distribution of small RNA reads on the precursors of the novel miRNAs depicted in Fig. 1a

Identification of miRNA targets in alfalfa

Although the alfalfa is one of the important legumes agronomically, the genome sequencing and annotations are not available so far. Due to this, studies have utilized the well-studied and closely related M. truncatula genome annotations as a model for alfalfa studies. The nucleotide identity for some genes was greater than 97% between M. sativa and M. truncatula [34]). Thus, using M. truncatula transcript annotations can facilitate identification of miRNA targets in alfalfa. We used SeqTar algorithm (Zheng et al., 2012) to identify miRNA targets by allowing a maximum of 4 mismatches between miRNAs and their potential target transcripts.

Previous studies have revealed that conserved miRNAs are strongly associated with the regulation of genes that encode transcription factors [35]. These transcription factors in turn regulate key developmental processes and pathways in plants. Degradome sequencing has been very effective in identifying plant miRNA targets. Besides identifying the conserved targets, this approach can also identify non-conserved targets for the conserved miRNAs [28, 36, 37]. Degradome sequencing was used in this study to identify the cleaved mRNA fragments corresponding to the miRNA recognition sites in all three alfalfa genotypes. Approximately 30 million degradome reads were obtained from the transcripts of each of the alfalfa genotypes (Table 4) and these reads were analysed using SeqTar program. In total, we have identified 69 targets for 31 miRNA families that included 16 highly conserved families (Table 5). With respect to the conserved miRNAs, 33 targets for 16 conserved miRNA families were identified (Table 5). The known targets for miR162, miR165/166, miR398 and miR399 families were not identified in this study. Although miR165/166 family is the most abundantly expressed as scored from their read frequencies in almost all libraries but the cleaved fragments from the HD-Zip target transcripts were not recovered from degradome libraries of alfalfa genotypes.

Table 4 Mapping of the reads obtained from the degradome libraries
Table 5 miRNA targets identified in the degradome libraries generated from three alfalfa genotypes. #Mis. is number of mismatches on the miRNA complementary site; Valid reads is Reads corresponding to the expected cleavage site; Total reads is Total reads mapped to the cDNA of the gene; Percent is Percent reads at the expected cleavage site

The identified miRNA targets in all three genotypes include mainly transcription factors. Specifically, five members of the squamosa promoter-binding-like protein (SPL) targeted by the miR156 family, five members of the auxin response factors targeted by both miR160 and miR167 families, five members of the apetala2 (AP2)-domain containing transcription factors, four members of the growth-regulating factor (GRFs) family targeted by miR396, two members of the TCP family transcription factors targeted by miR319, and, a NAC domain-containing transcription factor-like protein (NAC) targeted by miR164 [35]. Additionally, transcripts encoding Argonaute targeted by miR168, laccase targeted by miR397, and three plantacyanin containing proteins targeted by miR408 were also identified. Although evidence indicates that that miR398 targets Cu/Zn superoxide dismutases and a copper chaperone for the superoxide dismutases (CCS) in plants [28, 38] these relationships were not apparent in the data from this study. On the other hand, we have identified three potentially non-conserved targets (plastocyanin, protein disulphide isomerase and a hypothetical protein) for miR398 in three alfalfa genotypes. In addition to the GRFs, our analyses revealed potential non-conserved targets for miR396 including TNP1, Ulp1 protease and hypothetical proteins (Table 5).

The analyses of legume-specific miRNAs and their targets have revealed an interesting miRNA: target networks between the miRNAs and the NBS-LRR genes [14, 16, 18, 32]. In this study, we identified NBS-LRR disease resistance genes as targets for four different miRNA families including miR482, miR1507, miR1510 and miR5213 in alfalfa (Table 5).

Degradome analyses has also been utilized to identify potential targets for several non-conserved miRNAs or miRNAs that are present only in closely related species such as the M. truncatula. To increase the confidence in identification of targets for the non-conserved miRNAs that are usually expressed at low abundances and the cleavage frequencies on those targets are relatively low, we considered as ‘targets’ only those for which the cleavages were detected at least in two of the three alfalfa genotypes. The transcripts for Medtr6g053240.1 (F-box protein interaction domain protein) had a cleavage frequency of approximately 75% and were targeted by the miR2643 in NF08ALF06 genotype. Additionally, two other F-box protein interaction domain protein genes were also identified as targets for miR2643 in alfalfa genotypes (Table 5). These results suggest that the F-box protein interaction domain protein family are regulated by this potential legume-specific miRNA. Another notable observation is that 6 different genes identified as potential targets for miR7696, and the cleavage frequency of a particular target gene (hypothetical protein, Medtr3g112250.1) was more abundant in all three alfalfa genotypes (Table 5).

Because some of the miRNA-stars are also highly expressed, we scrutinized the degradome reads for potential cleavages on the transcripts that are complementary to the miRNA-stars. This analysis has identified potential targets for at least four conserved miRNAs. Specifically, miR156-star targets a heat shock transcription factor, miR164-star targets a protein transporter Sec61 subunit alpha-like protein, miR167-star targets a GRAS family transcription factor, and, miR482-star targets an auxin response factor 1 in alfalfa (Table 5).

The confirmed targets of conserved miRNAs are known to regulate diverse developmental processes in the lifecycle of plants. For example, the SPL transcription factors (target of miR156) which regulate the transition from juvenile to adult phase of the life cycle in land plants [39]. Auxin receptors (TIR1 proteins) and ARFs targeted by miR393 and miR160, miR167, are components of the auxin signalling pathway that regulates several aspects of plant growth and development. The roles of NAC factors (targeted by miR164) include shoot meristem initiation and later root formation in Arabidopsis [40, 41]. Similarly, TCP family transcription factors have several different roles including regulating leaf morphogenesis [42, 43]. In Arabidopsis, seven out of nine GRFs are known targets for miR396 [44], and we have identified four GRFs as targets for miR396 in alfalfa (Table 5). By interacting with its coactivators called GRF-interacting factors (GIFs), this regulatory network (miR396-GRFs-GIFs) regulate leaf size, leaf growth and senescence in Arabidopsis [44]. The known targets for miR397 include laccase, which is involved in oxidative polymerization of lignin in plants [45]. Similarly, miR408 is targeting a family of plantacyanins, which could function in shuttling electron-transfer between proteins [46, 47].

The miR398 family is known to target CSDs and a copper chaperone for superoxide dismutase (CCS) genes in plants [28, 38]. In this study, we have identified plastocyanin-domain like proteins (plastocyanin is an essential electron carrier which shuttles the electrons between cytochrome b6f and PS I) represents a novel target for miR398. Protein disulphide isomerase (PDI) is a member of a family of dithiol/disulfide oxidoreductases, the thioredoxin superfamily, which functions in the formation of disulphide linkage between the cysteine residues for proper protein folding [48]. Our degradome analyses confirms that PDI represents a novel target for miR398 in alfalfa (Table 5). The other confirmed miRNA target transcripts include Leucine rich repeat resistance (LRR) proteins (TIR-NBS-LRR and CC-NBS-LRR) that play important roles in plant pathogen recognition and activation of plant innate immune responses [14, 16, 18, 32]. Yet another interesting target include the F-box protein interaction domain proteins that are regulated by miR2643, one of the very abundantly expressed miRNA in alfalfa.


The analyses of small RNA libraries from the whole plants, shoots and roots resulted in the identification of 100 miRNA families that included highly conserved miRNAs as well as miRNAs that are at least conserved between M. truncatula and alfalfa. The conserved miRNA profiles share some similarities and a few differences between genotypes and types of tissues (roots and shoots). The tissue-specific profiles were used to identify miRNAs that are highly abundant as well as those miRNAs that are expressed at low levels. Additionally, 17 novel miRNAs with varying levels of expression were also identified in alfalfa. The present study also reports identification of 69 targets for 31 miRNA families. In addition to the conserved targets for conserved miRNAs, a few non-conserved targets such as the PDI for miR398 were confirmed. Similarly, miR2643 is targeting three transcripts encoding F-box protein interaction domain containing proteins in alfalfa. In summary, the results from this study have increased our understanding of miRNAs and miRNA-mediated gene regulation in alfalfa that could result in potential tangible targets for practical applications in alfalfa and related legume species to increase biomass yield and address abiotic and biotic limitations to agricultural productivity.


  1. 1.

    Sheaffer CC, Martin NP, Lamb JFS, Cuomo GR, Grimsbo Jewett J, Quering SR. Leaf and stem properties of alfalfa entries. Agron J. 2000;92(4):733–9.

    Article  Google Scholar 

  2. 2.

    Doole GJ, Pannell DJ. Role and value of including lucerne (Medicago sativa L.) phases in crop rotations for the management of herbicide-resistant Lolium rigidum in Western Australia. Crop Prot. 2008;27:497–504.

    Article  Google Scholar 

  3. 3.

    Latta RA, Blacklow LJ, Cocks PS. Comparative soil water, pasture production, and crop yields in phase farming systems with lucerne and annual pasture in Western Australia. Aust J Agric Res. 2001;52:295–303.

    Article  Google Scholar 

  4. 4.

    Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.

    CAS  Article  Google Scholar 

  5. 5.

    Sunkar R, Zhu JK. MicroRNAs and short-interfering RNAs in plants. J Integr Plant Biol. 2007;49:817–26.

    CAS  Article  Google Scholar 

  6. 6.

    Rogers K, Chen X. Biogenesis, turnover, and mode of action of plant microRNAs. Plant Cell. 2013;25:2383–99.

    CAS  Article  Google Scholar 

  7. 7.

    Bologna NG, Voinnet O. The diversity, biogenesis, and activities of endogenous silencing small RNAs in Arabidopsis. Ann Rev Plant Biol. 2014;65:473–503.

    CAS  Article  Google Scholar 

  8. 8.

    Chen X. Small RNAs and their roles in plant development. Ann Rev Cell and Dev Biol. 2009;25:21–44.

    Article  Google Scholar 

  9. 9.

    Sunkar R. MicroRNAs with macro effects on plant stress responses. Semin Cell Dev Biol. 2010;21(8):805–11.

    CAS  Article  Google Scholar 

  10. 10.

    Sunkar R, Li Y-F, Jagadeeswaran G. Functions of microRNAs in plant stress responses. Trends Plant Sci. 2012;17:196–203.

    CAS  Article  Google Scholar 

  11. 11.

    Sunkar R, Jagadeeswaran G. In silico identification of conserved miRNAs in large number of diverse plant species. BMC Plant Biol. 2008;8:37.

    Article  Google Scholar 

  12. 12.

    Axtell MJ, Snyder JP, Bartel DP. Common functions for diverse small RNAs of land plants. Plant Cell. 2007;19:1750–69.

    CAS  Article  Google Scholar 

  13. 13.

    Szittya G, Moxon S, Santos DM, Jing R, Fevereiro MP, Moulton V, Dalmay T. High-throughput sequencing of Medicago truncatula short RNAs identifies eight new miRNA families. BMC Genomics. 2008;9(1):593.

    Article  Google Scholar 

  14. 14.

    Jagadeeswaran G, Zheng Y, Li YF, Shukla LI, Matts J, Hoyt P, Sunkar R. Cloning and characterization of small RNAs from Medicago truncatula reveals four novel legume-specific microRNA families. New Phytol. 2009;184(1):85–98.

    CAS  Article  Google Scholar 

  15. 15.

    Lelandais-Brière C, Naya L, Sallet E, Calenge F, Frugier F, Hartmann C, et al. Genome-wide Medicagotruncatula small RNA analysis revealed novel microRNAs and isoforms differentially regulated in roots and nodules. Plant Cell. 2009;21(9):2780–96.

    Article  Google Scholar 

  16. 16.

    Zhai J, Jeong D, DePaoli E, Park S, Rosen BD, Li Y, Gonzalez AJ, Yan Z, Kitto SL, Grusak MA, Stacey G, Cook DR, Green PJ, Sherrier DJ, Meyers BC. MicroRNAs as master regulators of the plant NB-LRR defense gene family via the production of phased, trans- acting siRNAs. Genes Dev. 2001;25:2540–53.

    Article  Google Scholar 

  17. 17.

    Chi X, Yang Q, Chen X, Wang J, Pan L, Chen M, He Y, Laiang X, Yu S. Identification and characterization of microRNAs from peanut (Arachis hypogaea L.) by high-throughput sequencing. PLoS One. 2011;6(11):e27530.

    CAS  Article  Google Scholar 

  18. 18.

    Srivastava S, Zheng Y, Kudapa H, Jagadeeswaran G, Hivrale V, Varshney RK, Sunkar R. High throughput sequencing of small RNA component of leaves and inflorescence revealed conserved and novel miRNAs as well as phasiRNA loci in chickpea. Plant Sci. 2015;235:46–57.

    CAS  Article  Google Scholar 

  19. 19.

    Zheng Y, Hivrale V, Valliyodan B, Lelandais-Brière C, Farmer AD, May GD, Crespi M, Nguyen HT, Sunkar R. Small RNA profiles in soybean primary root tips under water deficit. BMC Syst Biol. 2016;10(5):126.

    Article  Google Scholar 

  20. 20.

    Fan W, Zhang S, Du H, Sun X, Shi Y, Wang C. Genome-wide identification of different dormant Medicago sativa L. MicroRNAs in response to fall dormancy. PLoS One. 2014;9(12):e114612.

    Article  Google Scholar 

  21. 21.

    Li Y, Wan L, Bi W, Wan X, Li Z, Cao J, Tong Z, Xu H, He F, Lo X. Identification of drought-responsive microRNAs from roots and leaves of alfalfa by high-throughput sequencing. Genes. 2017;8:119.

    Article  Google Scholar 

  22. 22.

    Li Z, Xu H, Li Y, Wan X, Ma Z, Cao J, Li Z, He F, Wang Y, Wan L, Tong Z, Li X. Analysis of physiological and miRNA responses to pi deficiency in alfalfa (Medicago sativa L.). Plant Mol Biol. 2018;96:473–92.

    CAS  Article  Google Scholar 

  23. 23.

    Monteros MJ, Tang H, Ramaraj T, Devitt NP, Cameron CT, Bharti AK, Mudge J, Farmer AD, Motes CM, Town CD, Brummer EC, Udvardi M. Progress in sequencing the genome of tetraploid alfalfa (Medicago sativa L). San Diego: Plant and Animal Genome Conference. 2015. Jan 10–14; 2015.

    Google Scholar 

  24. 24.

    Khu DM, Reyno R, Brummer EC, Bouton JH, Monteros MJ. Identification of aluminum tolerance QTLs in tetraploid alfalfa (Medicago sativa L.). Crop Sci. 2013;53:148–63.

    CAS  Article  Google Scholar 

  25. 25.

    Khu DM, Reyno R, Brummer EC, Monteros MJ. Screening methods for aluminum tolerance in alfalfa (Medicago sativa L.). Crop Sci. 2012;52:161–7.

    CAS  Article  Google Scholar 

  26. 26.

    Jagadeeswaran J, Nimmakayala P, Zheng Y, Gowdu K, Reddy U, Sunkar R. Characterization of the small RNA component of leaves and fruits from four different cucurbit species. BMC Genomics. 2012;13:329.

    CAS  Article  Google Scholar 

  27. 27.

    Kozomara A, Griffiths-Jones S. miRBase: annotating high confidence microRNAs using deep sequencing data. Nucleic Acids Res. 2014;42:D68–73.

    CAS  Article  Google Scholar 

  28. 28.

    Li YF, Zheng Y, Addo-Quaye C, Zhang L, Saini A, Jagadeeswaran G, Sunkar R. Transcriptome-wide identification of microRNA targets in rice. The Plant J. 2010;62(5):742–59.

    CAS  Article  Google Scholar 

  29. 29.

    Zheng Y, Li Y-F, Sunkar R, Zhang W. SeqTar: an effective method for identifying MicroRNA guided cleavage sites from Degradome of Polyadenylated transcripts in plants. Nucleic Acids Res. 2012;40:e28.

    CAS  Article  Google Scholar 

  30. 30.

    Sunkar R, Zhou X, Zheng Y, Zhang W, Zhu JK. Identification of novel and candidate miRNAs in rice by high throughput sequencing. BMC Plant Biol. 2008;8(1):25.

    Article  Google Scholar 

  31. 31.

    Li Y, Zheng Y, Jagadeeswaran G, Sunkar R. Characterization of small RNAs and their target genes in wheat seedlings using sequencing-based approaches. Plant Sci. 2013;203:17–24.

    Article  Google Scholar 

  32. 32.

    Li F, Pignatta D, Bendix C, Brunkard JO, Cohn MM, Tung J, Sun H, Kumar P, Baker B. MicroRNA regulation of plant innate immune receptors. Proc Natl Acad Sci U S A. 2012;109:1790–5.

    CAS  Article  Google Scholar 

  33. 33.

    Kohli D, Joshi G, Deokar AA, Bhardwaj AR, Agarwal M, Katiyar-Agarwal S, Srinivasan R, Jain PK. Identification and characterization of wilt and salt stress-responsive microRNAs in chickpea through high-throughput sequencing. PLoS One. 2014;9(10):e108851.

    Article  Google Scholar 

  34. 34.

    Tesfaye M, Silverstein K, Bucciarelli B, Samac CD, Vance CP. The Affymetrix Medicago GeneChip® array is applicable for transcript analysis of alfalfa (Medicago sativa). Funct Plant Biol. 2006;33:783–8.

  35. 35.

    Jones-Rhoades MW, Bartel DP, Bartel B. MicroRNAs and their regulatory roles in plants. Annu Rev Plant Biol. 2006;57:19–53.

    CAS  Article  Google Scholar 

  36. 36.

    Addo-Quaye C, Eshoo TW, Bartel DP, Axtell MJ. Endogenous siRNA and miRNA targets identified by sequencing of the Arabidopsis degradome. Curr Biol. 2008;18(10):758–62.

    CAS  Article  Google Scholar 

  37. 37.

    German MA, Pillay M, Jeong D-HH, Hetawal A, Luo S, Janardhanan P, Kannan V, Rymarquis LA, Nobuta K, German R, et al. Global identification of microRNA-target RNA pairs by parallel analysis of RNA ends. Nat Biotechnol. 2008;26:941–6.

    CAS  Article  Google Scholar 

  38. 38.

    Sunkar R, Kapoor A, Zhu JK. Post transcriptional induction of two cu/Zn superoxide dismutase genes in Arabidopsis is mediated by down-regulation of miR398 and important for oxidative stress tolerance. Plant Cell. 2006;18:2051–65.

    CAS  Article  Google Scholar 

  39. 39.

    Poethig RS. Small RNAs and developmental timing in plants. Curr Opin Genet Dev. 2009;19:374–8.

    CAS  Article  Google Scholar 

  40. 40.

    Guo HS, Xie Q, Fei JF, Chua NH. MicroRNA directs mRNA cleavage of the transcription factor NAC1 to downregulate auxin signals for Arabidopsis lateral root development. Plant Cell. 2005;17:1376–86.

    CAS  Article  Google Scholar 

  41. 41.

    Laufs P, Peaucelle A, Morin H, Traas J. MicroRNA regulation of the CUC genes is required for boundary size control in Arabidopsis meristems. Development. 2004;131:4311–22.

    CAS  Article  Google Scholar 

  42. 42.

    Palatnik JF, Allen E, Wu X, Schommer C, Schwab R, Carrington JC, Weigel D. Control of leaf morphogenesis by microRNAs. Nature. 2003;425:257–63.

    CAS  Article  Google Scholar 

  43. 43.

    Koyama T, Sato F, Ohme-Takagi M. Roles of miR319 and TCP transcription factors in leaf development. Plant Physiol. 2017;175(2):874–85.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. 44.

    Debernardi JM, Mecchia MA, Vercruyssen L, Smaczniak C, Kaufmann K, Inze D, Rodriguez RE, Palatnik JF. Post-transcriptional control of GRF transcription factors by microRNA396 and GIF co-activators affects leaf size and longevity. Plant J. 2014;79:413–26.

    CAS  Article  Google Scholar 

  45. 45.

    Zhao Q, Nakashima J, Chen F, Yin YB, Fu CX, Yun JF, Shao H, Wang XQ, Wang ZY, Dixon RA. Laccase is necessary and nonredundant with peroxidase for lignin polymerization during vascular development in Arabidopsis. Plant Cell. 2013;25:3976–87.

    CAS  Article  Google Scholar 

  46. 46.

    Nersissian AM, Immoos C, Hill MG, Hart PJ, Williams G, Herrmann RG, Valentine JS. Uclacyanins, stellacyanins, and plantacyanins ape distinct subfamilies of phytocyanins: plant-specific mononuclear blue copper proteins. Protein Sci. 1998;7:1915–29.

    CAS  Article  Google Scholar 

  47. 47.

    Ma C, Burd S, Lers A. miR408 is involved in abiotic stress responses in Arabidopsis. Plant J. 2015;84:169–87.

    CAS  Article  Google Scholar 

  48. 48.

    Wilkinson B, Gilbert HF. Protein disulfide isomerase. Boiochem Biophys Acta. 2004;1699:35–44.

    CAS  Google Scholar 

Download references


Not applicable.


This research was funded by the Noble Research Institute and Forage Genetics International, a hatch grant from NIFA-0229360 (OKL02844) to RS and MM, and the National Natural Science Foundation of China (numbers 31460295 and 31760314) to YZ. This work was also partially supported by the Neustadt-Sarkeys Distinguished Professorship to RS. Publication costs are funded by the Oklahoma Agricultural Experiment Station, Oklahoma State University, Stillwater.

Availability of data and materials

The small RNA and degradome datasets generated and analyzed in the present study are available in the National Center for Biotechnology Information Gene Expression Omnibus (NCBI GEO) under accession number GSE119460 available at:

About this supplement

This article has been published as part of BMC Genomics Volume 19 Supplement 10, 2018: Proceedings of the 29th International Conference on Genome Informatics (GIW 2018): genomics. The full contents of the supplement are available online at

Author information




RS and MM conceived the idea and designed the study. CM and TH cultured the plants used in this study; RP isolated the RNA from samples and generated the small RNA libraries; YZ, SR, QW, SA and RS analyzed the small RNA and degradome libraries; RP, SA and RS wrote the manuscript; MM edited the manuscript. All authors reviewed and approved the final manuscript.

Corresponding authors

Correspondence to Yun Zheng or Ramanjulu Sunkar.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Pokoo, R., Ren, S., Wang, Q. et al. Genotype- and tissue-specific miRNA profiles and their targets in three alfalfa (Medicago sativa L) genotypes. BMC Genomics 19, 913 (2018).

Download citation