The MicroRNA156 system: A tool in plant biotechnology

Recent allopolyploidy alters Spartina microRNA expression in response to xenobiotic-induced stress

Environmental contamination by xenobiotics represents a major threat for natural ecosystems and public health. In response, xenobiotic detoxification is a fundamental trait of organisms for developmental plasticity and stress tolerance, but the underlying molecular mechanisms remain poorly understood in plants. To decipher this process, we explored the consequences of allopolyploidy on xenobiotic tolerance in the genus Spartina Schreb. Specifically, we focused on microRNAs (miRNAs) owing to their central function in the regulation of gene expression patterns, including responses to stress. Small RNA-Seq was conducted on the parents S. alterniflora and S. maritima, their F1 hybrid S. x townsendii and the allopolyploid S. anglica under phenanthrene-induced stress (phe), a model Polycyclic Aromatic Hydrocarbon (PAH) compound. Differentially expressed miRNAs in response to phe were specifically identified within species. In complement, the respective impacts of hybridization and genome doubling were detected, through changes in miRNA expression patterns between S. x townsendii, S. anglica and the parents. The results support the impact of allopolyploidy in miRNA-guided regulation of plant response to phe. In total, we identified 17 phe-responsive miRNAs in Spartina among up-regulated MIR156 and down-regulated MIR159. We also describe novel phe-responsive miRNAs as putative Spartina-specific gene expression regulators in response to stress. Functional validation using Arabidopsis (L.) Heynh. T-DNA lines inserted in homologous MIR genes was performed, and the divergence of phe-responsive miRNA regulatory networks between Arabidopsis and Spartina was discussed.

Genome-wide identification, phylogeny and expression analysis of the SPL gene family and its important role in salt stress in Medicago sativa L

BACKGROUND

SQUAMOSA promoter-binding protein-like (SPL) transcription factors are widely present in plants and are involved in signal transduction, the stress response and development. The SPL gene family has been characterized in several model species, such as A. thaliana and G. max. However, there is no in-depth analysis of the SPL gene family in forage, especially alfalfa (Medicago sativa L.), one of the most important forage crops worldwide.

RESULT

In total, 76 putative MsSPL genes were identified in the alfalfa genome with an uneven distribution. Based on their identity and gene structure, these MsSPLs were divided into eight phylogenetic groups. Seventy-three MsSPL gene pairs arose from segmental duplication events, and the MsSPLs on the four subgenomes of individual chromosomes displayed high collinearity with the corresponding M. truncatula genome. The prediction of the cis-elements in the promoter regions of the MsSPLs detected two copies of ABA (abscisic acid)-responsive elements (ABREs) on average, implying their potential involvement in alfalfa adaptation to adverse environments. The transcriptome sequencing of MsSPLs in roots and leaves revealed that 54 MsSPLs were expressed in both tissues. Upon salt treatment, three MsSPLs (MsSPL17, MsSPL23 and MsSPL36) were significantly regulated, and the transcription level of MsSPL36 in leaves was repressed to 46.6% of the control level.

CONCLUSION

In this study, based on sequence homology, we identified 76 SPL genes in the alfalfa. The SPLs with high identity shared similar gene structures and motifs. In total, 71.1% (54 of 76) of the MsSPLs were expressed in both roots and leaves, and the majority (74.1%) preferred underground tissues to aerial tissues. MsSPL36 in leaves was significantly repressed under salt stress. These findings provide comprehensive information regarding the SPB-box gene family for improve alfalfa tolerance to high salinity.

Alfalfa transcriptome profiling provides insight into miR156-mediated molecular mechanisms of heat stress tolerance

Heat is one of the major environmental stressors that negatively affects alfalfa production. Previously, we reported the role of microRNA156 (miR156) in heat tolerance, however, mechanisms and downstream genes involved in this process were not fully studied. To provide further insight, we compared an empty vector control and miR156 overexpressing alfalfa plants (miR156+) after exposing them to heat stress (40 °C) for 24h. We collected leaf samples for transcriptome analysis to illustrate the miR156-regualted molecular mechanisms underlying the heat stress response. A total of 3579 differentially expressed genes (DEG) were detected exclusively in miR156+ plants under heat stress using the Medicago sativa genome as reference. GO and KEGG analysis indicated that these DEGs were mainly involved in "polysaccharide metabolism", "response to chemical", "secondary metabolism", "carbon metabolism" and "cell cycle". Transcription factors predicted in miR156+ plants belonged to TCP family, MYB, ABA response element-binding factor, WRKY and heat shock transcription factor. We also identified two new SPL family gene member (SPL8a and SPL12a), putatively regulated by miR156. The present study provides comprehensive transcriptome profile of alfalfa, identifies a number of genes and pathways, and reveals a miR156-regulated network of mechanisms at the gene expression level to modulate heat responses in alfalfa.

Integrated Analysis of miRNAs Associated With Sugarcane Responses to Low-Potassium Stress

Sugarcane is among the most important global crops and a key bioenergy source. Sugarcane production is restricted by limited levels of available soil potassium (K⁺). The ability of plants to respond to stressors can be regulated by a range of microRNAs (miRNAs). However, there have been few studies regarding the roles of miRNAs in the regulation of sugarcane responses to K⁺-deficiency. To understand how these non-coding RNAs may influence sugarcane responses to low-K⁺ stress, we conducted expression profiling of miRNAs in sugarcane roots under low-K⁺ conditions via high-throughput sequencing. This approach led to the identification of 324 and 42 known and novel miRNAs, respectively, of which 36 were found to be differentially expressed miRNAs (DEMs) under low-K⁺ conditions. These results also suggested that miR156-x/z and miR171-x are involved in these responses as potential regulators of lateral root formation and the ethylene signaling pathway, respectively. A total of 705 putative targets of these DEMs were further identified through bioinformatics predictions and degradome analyses, and GO and KEGG enrichment analyses revealed these target mRNAs to be enriched for catalytic activity, binding functions, metabolic processes, plant hormone signal transduction, and mitogen-activated protein kinase (MAPK) signaling. In summary, these data provide an overview of the roles of miRNAs in the regulation of sugarcane response to low-K⁺ conditions.

Genome-wide identification, phylogenetic analysis, and expression analysis of the SPL gene family in orchardgrass (Dactylis glomerata L.)

SPL (SQUAMOSA promoter binding protein-like) is a plant-specific transcription factor family that contains the conserved SBP domain, which plays a vital role in the vegetative-to-reproductive phase transition, flowering development and regulation, tillering/branching, and stress responses. Although the SPL family has been identified and characterized in various plant species, limited information about it has been obtained in orchardgrass, which is a critical forage crop worldwide. In this study, 17 putative DgSPL genes were identified among seven chromosomes, and seven groups that share similar gene structures and conserved motifs were determined by phylogenetic analysis. Of these, eight genes have potential target sites for miR156. cis-Element and gene ontology annotation analysis indicated DgSPLs may be involved in regulating development and abiotic stress responses. The expression patterns of eight DgSPL genes at five developmental stages, in five tissues, and under three stress conditions were determined by RNA-seq and qRT-PCR. These assays indicated DgSPLs are involved in vegetative-to-reproductive phase transition, floral development, and stress responses. The transient expression analysis in tobacco and heterologous expression assays in yeast indicated that miR156-targeted DG1G01828.1 and DG0G01071.1 are nucleus-localized proteins, that may respond to drought, salt, and heat stress. Our study represents the first systematic analysis of the SPL family in orchardgrass. This research provides a comprehensive assessment of the DgSPL family, which lays the foundation for further examination of the role of miR156/DgSPL in regulating development and stress responses in forages grasses.

Involvement of the miR156/SPL module in flooding response in Medicago sativa

The highly conserved plant microRNA, miR156, affects plant development, metabolite composition, and stress response. Our previous research revealed the role of miR156 in abiotic stress response in Medicago sativa exerted by downregulating SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE transcription factors. Here we investigated the involvement and possible mechanism of action of the miR156/SPL module in flooding tolerance in alfalfa. For that, we used miR156 overexpressing, SPL13RNAi, flood-tolerant (AAC-Trueman) and -sensitive (AC-Caribou) alfalfa cultivars exposed to flooding. We also used Arabidopsis ABA insensitive (abi1-2, abi5-8) mutants and transgenic lines with either overexpressed (KIN10-OX1, KIN10-OX2) or silenced (KIN10RNAi-1, KIN10RNAi-2) catalytic subunit of SnRK1 to investigate a possible role of ABA and SnRK1 in regulating miR156 expression under flooding. Physiological analysis, hormone profiling and global transcriptome changes revealed a role for miR156/SPL module in flooding tolerance. We also identified nine novel alfalfa SPLs (SPL1, SPL1a, SPL2a, SPL7, SPL7a, SPL8, SPL13a, SPL14, SPL16) responsive to flooding. Our results also showed a possible ABA-dependent SnRK1 upregulation to enhance miR156 expression, resulting in downregulation of SPL4, SPL7a, SPL8, SPL9, SPL13, and SPL13a. We conclude that these effects induce flooding adaptive responses in alfalfa and modulate stress physiology by affecting the transcriptome, ABA metabolites and secondary metabolism.

Overexpression of miR156 and Silencing SPL6RNAi and SPL13RNAi Genes in Medicago sativa on the Changes of Carbohydrate Physiochemical, Fermentation, and Nutritional Profiles

This study aimed to explore the comparative effects of overexpressing miR156 with individually silencing SPL6RNAi and SPL13RNAi genes on carbohydrate physiochemical, fermentation, and nutritional profiles of alfalfa (Medicago sativa). Three sub-genotypes of miR156 overexpressed (miR156 OE), SPL6RNAi, and SPL13RNAi grown with the wild type (WT) in a greenhouse were harvested 3 times at an early vegetative stage. Samples were freeze-dried, ground, and analyzed for carbohydrate nutritional profiles in terms of chemical composition, CNCPS fractions, energetic values, in vitro degradation, and fermentation characteristics. Results showed that miR156 OE had lower fiber and higher energy compared to all other genotypes. Moreover, miR156 OE had higher starch compared to SPL13RNAi and higher DM degradation compared to WT and SPL13RNAi. In conclusion, overexpression of miR156 decreased the fiber content of alfalfa but increased energy and DM degradation. SPL6RNAi was more similar to miR156 OE alfalfa in chemical composition and degradation, indicating that the SPL6RNAi gene plays an important role in the miR156 overexpression event.

The Impact of Zinc Oxide Nanoparticles on Cytotoxicity, Genotoxicity, and miRNA Expression in Barley (Hordeum vulgare L.) Seedlings

Zinc oxide nanoparticles are one of the most commonly engineered nanomaterials and necessarily enter the environment because of the large quantities produced and their widespread application. Understanding the impacts of nanoparticles on plant growth and development is crucial for the assessment of probable environmental risks to food safety and human health, because plants are a fundamental living component of the ecosystem and the most important source in the human food chain. The objective of this study was to examine the impact of different concentrations of zinc oxide nanoparticles on barley Hordeum vulgare L. seed germination, seedling morphology, root cell viability, stress level, genotoxicity, and expression of miRNAs. The results demonstrate that zinc oxide nanoparticles enhance barley seed germination, shoot/root elongation, and H₂O₂ stress level and decrease root cell viability and genomic template stability and up- and downregulated miRNAs in barley seedlings.

Characterizing the Role of the miR156-SPL Network in Plant Development and Stress Response

MicroRNA (miRNA) is a short, single-stranded, non-coding RNA found in eukaryotic cells that can regulate the expression of many genes at the post-transcriptional level. Among various plant miRNAs with diverse functions, miR156 plays a key role in biological processes, including developmental regulation, immune response, metabolic regulation, and abiotic stress. MiRNAs have become the regulatory center for plant growth and development. MicroRNA156 (miR156) is a highly conserved and emerging tool for the improvement of plant traits, including crop productivity and stress tolerance. Fine-tuning of squamosa promoter biding-like (SPL) gene expression might be a useful strategy for crop improvement. Here, we studied the regulation of the miR156 module and its interaction with SPL factors to understand the developmental transition of various plant species. Furthermore, this review provides a strong background for plant biotechnology and is an important source of information for further molecular breeding to optimize farming productivity.

The interplay between miR156/SPL13 and DFR/WD40-1 regulate drought tolerance in alfalfa

BACKGROUND

Developing Medicago sativa L. (alfalfa) cultivars tolerant to drought is critical for the crop's sustainable production. miR156 regulates various plant biological functions by silencing SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors.

RESULTS

To understand the mechanism of miR156-modulated drought stress tolerance in alfalfa we used genotypes with altered expression levels of miR156, miR156-regulated SPL13, and DIHYDROFLAVONOL-4-REDUCTASE (DFR) regulating WD40-1. Previously we reported the involvement of miR156 in drought tolerance, but the mechanism and downstream genes involved in this process were not fully studied. Here we illustrate the interplay between miR156/SPL13 and WD40-1/DFR to regulate drought stress by coordinating gene expression with metabolite and physiological strategies. Low to moderate levels of miR156 overexpression suppressed SPL13 and increased WD40-1 to fine-tune DFR expression for enhanced anthocyanin biosynthesis. This, in combination with other accumulated stress mitigating metabolites and physiological responses, improved drought tolerance. We also demonstrated that SPL13 binds in vivo to the DFR promoter to regulate its expression.

CONCLUSIONS

Taken together, our results reveal that moderate relative miR156 transcript levels are sufficient to enhance drought resilience in alfalfa by silencing SPL13 and increasing WD40-1 expression, whereas higher miR156 overexpression results in drought susceptibility.

Roles of the SPL gene family and miR156 in the salt stress responses of tamarisk (Tamarix chinensis)

BACKGROUND

Accumulating evidences show that SPLs are crucial regulators of plant abiotic stress tolerance and the highly conserved module miR156/SPL appears to balance plant growth and stress responses. The halophyte Tamarix chinensis is highly resistant to salt tress. SPLs of T. chinensis (TcSPLs) and theirs roles in salt stress responses remain elusive.

RESULTS

In this study, we conducted a systematic analysis of the TcSPLs gene family including 12 members belonging to 7 groups. The physicochemical properties and conserved motifs showed divergence among groups and similarity in each group. The microRNA response elements (MREs) are conserved in location and sequence, with the exception of first MRE within TcSPL5. The miR156-targeted SPLs are identified by dual-luciferase reporter assay of MRE-miR156 interaction. The digital expression gene profiles cluster suggested potential different functions of miR156-targeted SPLs vs non-targeted SPLs in response to salt stress. The expression patterns analysis of miR156-targeted SPLs with a reverse expression trend to TcmiR156 suggested 1 h (salt stress time) could be a critical time point of post-transcription regulation in salt stress responses.

CONCLUSIONS

Our work demonstrated the post-transcription regulation of miR156-targeted TcSPLs and transcription regulation of non-targeted TcSPLs in salt stress responses, and would be helpful to expound the miR156/SPL-mediated molecular mechanisms underlying T. chinensis salt stress tolerance.

Alfalfa response to heat stress is modulated by microRNA156

Heat stress and extreme temperatures negatively affect plant development by disrupting regular cellular and biochemical functions, ultimately leading to reduced crop production. Alfalfa (Medicago sativa) is an important forage crop grown worldwide as forage for livestock feed. Limiting the effects of abiotic stress by developing alfalfa cultivars that are stress tolerant would help mitigate losses to crop production. Members of the microRNA156 (miR156) family regulate the Squamosa Promoter-Binding Protein-Like (SPL) genes that in turn impact plant growth and development by regulating downstream genes in response to various abiotic stresses. In this study, alfalfa with miR156 overexpression and SPL13 RNAi knockdown show increased tolerance to heat stress (40°C). Transgenic plants show high water potential and increased non-enzymatic antioxidant content under heat stress. Moreover, anthocyanin content and chlorophyll abundance were increased under stress. Expression of some important transcription factors and downstream genes involved in abiotic stress response were altered in miR156-overexpressing genotypes under heat. Taken together, our results demonstrate that the miR156/SPL13 network contributes to improving heat stress tolerance in alfalfa.

Salinity-associated microRNAs and their potential roles in mediating salt tolerance in rice colonized by the endophytic root fungus Piriformospora indica

Piriformospora indica (P. indica), an endophytic root fungus, supports the growth and enhanced tolerance of plants to biotic and abiotic stresses. Several recent studies showed the significant role of small RNA (sRNA) molecules including microRNAs (miRNAs) in plant adaption to environmental stress, but little is known concerning the symbiosis-mediated salt stress tolerance regulated at miRNAs level. The overarching goal of this research is to elucidate the impact of miRNAs in regulating the P. indica-mediated salt tolerance in rice. Applying sRNA-seq analysis led to identify a set of 547 differentially abundant miRNAs in response to P. indica inoculation and salt stress. These included 206 rice-specific and 341 previously known miRNAs from other plant species. In silico analysis of miRNAs predictions of the differentially abundant miRNAs led to identifying of 193 putatively target genes, most of which were encoded either genes or transcription factors involved in nutrient uptake, sodium ion transporters, growth regulators, and auxin- responsive proteins. The rice-specific miRNAs targeted the transcription factors involved in the import of potassium ions into the root cells, the export of sodium ions, and plant growth and development. Interestingly, P. indica affected the differential abundance of miRNAs regulated genes and transcription factors linked to salt stress tolerance. Our data helps to understand the molecular basis of salt stress tolerance mediated by symbionts in plant and the potential impact of miRNAs for genetic improvement of rice varieties for tolerance to salt stress.

SPL13 regulates shoot branching and flowering time in Medicago sativa

Our results show SPL13 plays a crucial role in regulating vegetative and reproductive development in Medicago sativa L. (alfalfa), and that MYB112 is targeted and downregulated by SPL13 in alfalfa. We previously showed that transgenic Medicago sativa (alfalfa) plants overexpressing microRNA156 (miR156) show a bushy phenotype, reduced internodal length, delayed flowering time, and enhanced biomass yield. In alfalfa, transcripts of seven SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) transcription factors, including SPL13, are targeted for cleavage by miR156. Thus, association of each target SPL gene to a trait or set of traits is essential for developing molecular markers for alfalfa breeding. In this study, we investigated SPL13 function using SPL13 overexpression and silenced alfalfa plants. Severe growth retardation, distorted branches and up-curled leaves were observed in miR156-impervious 35S::SPL13m over-expression plants. In contrast, more lateral branches and delayed flowering time were observed in SPL13 silenced plants. SPL13 transcripts were predominantly present in the plant meristems, indicating that SPL13 is involved in regulating shoot branch development. Accordingly, the shoot branching-related CAROTENOID CLEAVAGE DIOXYGENASE 8 gene was found to be significantly downregulated in SPL13 RNAi silencing plants. A R2R3-MYB gene MYB112 was also identified as being directly silenced by SPL13 based on Next Generation Sequencing-mediated transcriptome analysis and chromatin immunoprecipitation assays, suggesting that MYB112 may be involved in regulating alfalfa vegetative growth.

MicroRNA156 improves drought stress tolerance in alfalfa (Medicago sativa) by silencing SPL13

Alfalfa (Medicago sativa) is an important forage crop that is often grown in areas that frequently experience drought and water shortage. MicroRNA156 (miR156) is an emerging tool for improving various traits in plants. We tested the role of miR156d in drought response of alfalfa, and observed a significant improvement in drought tolerance of miR156 overexpression (miR156OE) alfalfa genotypes compared to the wild type control (WT). In addition to higher survival and reduced water loss, miR156OE genotypes also maintained higher stomatal conductance compared to WT during drought stress. Furthermore, we observed an enhanced accumulation of compatible solute (proline) and increased levels of abscisic acid (ABA) and antioxidants in miR156OE genotypes. Similarly, alfalfa plants with reduced expression of miR156-targeted SPL13 showed reduced water loss and enhanced stomatal conductance, chlorophyll content and photosynthetic assimilation. Several genes known to be involved in drought tolerance were differentially expressed in leaf and root of miR156 overexpression plants. Taken together, our findings reveal that miR156 improves drought tolerance in alfalfa at least partially by silencing SPL13.

The Use of Gene Modification and Advanced Molecular Structure Analyses towards Improving Alfalfa Forage

Alfalfa is one of the most important legume forage crops in the world. In spite of its agronomic and nutritive advantages, alfalfa has some limitations in the usage of pasture forage and hay supplement. High rapid degradation of protein in alfalfa poses a risk of rumen bloat to ruminants which could cause huge economic losses for farmers. Coupled with the relatively high lignin content, which impedes the degradation of carbohydrate in rumen, alfalfa has unbalanced and asynchronous degradation ratio of nitrogen to carbohydrate (N/CHO) in rumen. Genetic engineering approaches have been used to manipulate the expression of genes involved in important metabolic pathways for the purpose of improving the nutritive value, forage yield, and the ability to resist abiotic stress. Such gene modification could bring molecular structural changes in alfalfa that are detectable by advanced structural analytical techniques. These structural analyses have been employed in assessing alfalfa forage characteristics, allowing for rapid, convenient and cost-effective analysis of alfalfa forage quality. In this article, we review two major obstacles facing alfalfa utilization, namely poor protein utilization and relatively high lignin content, and highlight genetic studies that were performed to overcome these drawbacks, as well as to introduce other improvements to alfalfa quality. We also review the use of advanced molecular structural analysis in the assessment of alfalfa forage for its potential usage in quality selection in alfalfa breeding.

An Insight into microRNA156 Role in Salinity Stress Responses of Alfalfa

Salinity is one of the major abiotic stresses affecting alfalfa productivity. Developing salinity tolerant alfalfa genotypes could contribute to sustainable crop production. The functions of microRNA156 (miR156) have been investigated in several plant species, but so far, no studies have been published that explore the role of miR156 in alfalfa response to salinity stress. In this work, we studied the role of miR156 in modulating commercially important traits of alfalfa under salinity stress. Our results revealed that overexpression of miR156 increased biomass, number of branches and time to complete growth stages, while it reduced plant height under control and salinity stress conditions. We observed a miR156-related reduction in neutral detergent fiber under non-stress, and acid detergent fiber under mild salinity stress conditions. In addition, enhanced total Kjeldahl nitrogen content was recorded in miR156 overexpressing genotypes under severe salinity stress. Furthermore, alfalfa genotypes overexpressing miR156 exhibited an altered ion homeostasis under salinity conditions. Under severe salinity stress, miR156 downregulated SPL transcription factor family genes, modified expression of other important transcription factors, and downstream salt stress responsive genes. Taken together, our results reveal that miR156 plays a role in mediating physiological and transcriptional responses of alfalfa to salinity stress.

SQUAMOSA PROMOTER BINDING PROTEIN-LIKE 2 controls floral organ development and plant fertility by activating ASYMMETRIC LEAVES 2 in Arabidopsis thaliana

A network of genes is coordinately expressed to ensure proper development of floral organs and fruits, which are essential for generating new offspring in flowering plants. In Arabidopsis thaliana, microRNA156 (miR156) plays a role in regulating the development of flowers and siliques by targeting members of the SQUAMOSA PROMOTER BINDING PROTEIN-LIKE (SPL) gene family. Despite the important roles of the miR156/SPL network, our understanding of its downstream genes that are involved in floral organ and silique growth is still incomplete. Here, we report that the miR156/SPL2 regulatory pathway regulates pollen production, fertility rate, and the elongation of floral organs, including petals, sepals, and siliques in Arabidopsis. Transgenic plants exhibiting both overexpression of miR156 and dominant-negative alleles of SPL2 had reduced ASYMMETRIC LEAVES 2 (AS2) transcript levels in their siliques. Furthermore, their fertility phenotype was similar to that of the AS2 loss-of-function mutant. We also demonstrate that the SPL2 protein binds to the 5'UTR of the AS2 gene in vivo, indicating that AS2 is directly regulated by SPL2. Our results suggest that the miR156/SPL2 pathway affects floral organs, silique development and plant fertility, as well as directly regulates AS2 expression.

Comparative transcriptome investigation of global gene expression changes caused by miR156 overexpression in Medicago sativa

Background

Medicago sativa (alfalfa) is a low-input forage and potential bioenergy crop, and improving its yield and quality has always been a focus of the alfalfa breeding industry. Transgenic alfalfa plants overexpressing a precursor of alfalfa microRNA156 (MsmiR156) were recently generated by our group. These plants (miR156OE) showed enhanced biomass yield, reduced internodal length, increased shoot branching and trichome density, and a delay in flowering time. Transcripts of three SQUAMOSA-PROMOTER BINDING PROTEIN-LIKE (SPL) genes (MsSPL6, MsSPL12, and MsSPL13) were found to be targeted for cleavage by MsmiR156 in alfalfa.

Results

To further illustrate the molecular mechanisms underlying the effects of miR156 in alfalfa, two miR156OE genotypes (A11a and A17) were subjected to Next Generation RNA Sequencing with Illumina HiSeq. More than 1.11 billion clean reads were obtained from our available sequenced samples. A total of 160,472 transcripts were generated using Trinity de novo assembly and 4,985 significantly differentially expressed genes were detected in miR156OE plants A11a and A17 using the Medicago truncatula genome as reference. A total of 17 genes (including upregulated, downregulated, and unchanged) were selected for quantitative real-time PCR (qRT-PCR) validation, which showed that gene expression levels were largely consistent between qRT-PCR and RNA-Seq data. In addition to the established SPL genes MsSPL6, MsSPL12 and MsSPL13, four new SPLs; MsSPL2, MsSPL3, MsSPL4 and MsSPL9 were also down-regulated significantly in both miR156OE plants. These seven SPL genes belong to genes phylogeny clades VI, IV, VIII, V and VII, which have been reported to be targeted by miR156 in Arabidopsis thaliana. The gene ontology terms characterized electron transporter, starch synthase activity, sucrose transport, sucrose-phosphate synthase activity, chitin binding, sexual reproduction, flavonoid biosynthesis and lignin catabolism correlate well to the phenotypes of miR156OE alfalfa plants.

Conclusions

This is the first report of changes in global gene expression in response to miR156 overexpression in alfalfa. The discovered miR156-targeted SPL genes belonging to different clades indicate miR156 plays fundamental and multifunctional roles in regulating alfalfa plant development.

Electronic supplementary material

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