Genome-Wide Identification and Expression Analysis of the SWEET Gene Family in Annual Alfalfa (Medicago polymorpha).

SWEET (Sugars will eventually be exported transporter) proteins are a group of sugar transporters that are involved in sugar efflux, phloem loading, reproductive development, plant senescence, and stress responses. In this study, 23 SWEET transporter members were identified in the Medicago polymorpha genome, heterogeneously distributed on seven chromosomes. These MpSWEET genes were divided into four subfamilies, which showed similar gene structure and motif composition within the same subfamily. Seventeen MpSWEET genes encode seven transmembrane helices (TMHs), and all MpSWEET proteins possess conserved membrane domains and putative serine phosphorylation sites. Four and three pairs of MpSWEET genes were predicted to be segmentally and tandemly duplicated, respectively, which may have contributed to their evolution of M. polymorpha. The results of microarray and RNA-Seq data showed that some MpSWEET genes were specifically expressed in disparate developmental stages (including seedling stage, early flowering stage, and late flowering stage) or tissues such as flower and large pod. Based on protein network interaction and expression patterns of MpSWEET genes, six MpSWEET genes were selected for further quantitative real-time PCR validation in different stress treatments. qRT-PCR results showed that MpSWEET05, MpSWEET07, MpSWEET12, MpSWEET15, and MpSWEET21 were significantly upregulated for at least two of the three abiotic stress treatments. These findings provide new insights into the complex transcriptional regulation of MpSWEET genes, which facilitates future research to elucidate the function of MpSWEET genes in M. polymorpha and other legume crops.

Genome-Wide Identification and Expression Analysis of the Dof Transcription Factor in Annual Alfalfa Medicago polymorpha.

The Dof transcription factor is a plant-specific transcription gene family that plays various biological functions in plant development and stress response. However, no relevant research has been conducted on Medicago polymorpha. Here, 36 MpDof genes were identified in the M. polymorpha genome and further divided into 10 groups based on the comparative phylogenetic analysis. The essential information of MpDof genes, such as chromosomal localization, gene structure, conserved motifs, and selective pressures were systematically analyzed. All 36 MpDof genes were predicted to contain more cis-acting elements related to hormone response. MpDof24 and MpDof25 were predicted to interact with MpDof11 and MpDof26 to involve in the photoperiod blooms process. The MpDof genes showed a diverse expression pattern in different tissues. Notably, MpDof29 and MpDof31 were specifically expressed in the large pod and root, respectively, suggesting their crucial role in the pod and root development. qRT-PCR analysis indicated that the expression levels of MpDof10, MpDof25, MpDof26, and MpDof29 were obviously up-regulated under drought, salt, and cold stress. Collectively, genome-wide identification, evolutionary, and expression analysis of the Dof transcription gene family in M. polymorpha will provide new information to further understand and utilize the function of these Dof genes in Medicago plants.

Screening and identification of multiple abiotic stress responsive candidate genes based on hybrid-sequencing in Vicia sativa.

Common vetch is an important leguminous forage for both livestock fodder and green manure and has a tremendous latent capacity in a sustainable agroecosystem. In the present study, a comprehensive transcriptome analysis of the aboveground leaves and underground roots of common vetch under multiple abiotic stress treatments, including NaCl, drought, cold, and cold drought, was performed using hybrid-sequencing technology, i. e. single-molecule real-time sequencing technology (SMRT) and supplemented by next-generation sequencing (NGS) technology. A total of 485,038 reads of insert (ROIs) with a mean length of 2606 bp and 228,261 full-length nonchimeric (FLNC) reads were generated. After deduplication, 39,709 transcripts were generated. Of these transcripts, we identified 1059 alternative splicing (AS) events, 17,227 simple sequence repeats (SSRs), and 1647 putative transcription factors (TFs). Furthermore, 640 candidates long noncoding RNAs (lncRNAs) and 28,256 complete coding sequences (CDSs) were identified. In gene annotation analyses, a total of 38,826 transcripts (97.78%) were annotated in eight public databases. Finally, seven multiple abiotic stress-responsive candidate genes were obtained through gene expression, annotation information, and protein-protein interaction (PPI) networks. Our research not only enriched the structural information of FL transcripts in common vetch, but also provided useful information for exploring the molecular mechanism of multiple abiotic stress tolerance between aboveground and underground tissues in common vetch and related legumes.

Molecular characterization of the GH3 family in alfalfa under abiotic stress.

The phytohormone auxin plays a pivotal role in regulating plant growth, development, and abiotic stress responses by promptly controlling the expression of auxin response genes. The Gretchen Hagen3 (GH3) genes are a major early auxin response gene family that contribute to auxin homeostasis by conjugating excess auxins to amino acids. To our knowledge, a genome-wide investigation of the GH3 genes in alfalfa has never been reported. Here, we present a comprehensive bioinformatics analysis of the MsGH3 gene family, including chromosomal locations, phylogenetic relationships, gene structures, conserved motifs and Gene Ontology annotation. Interestingly, the analysis revealed 31 MsGH3 genes in the alfalfa genome. These genes were classified phylogenetically into the GH3-I, GH3-II, and GH3-III subgroups. Additionally, the data analysis showed that most MsGH3 genes are tissue specific and responsive to environmental stress-related hormones. Furthermore, the analysis of cis elements and potential biological functions revealed that the MsGH3 genes play potential roles in various stress responses. Notably, qRT-PCR results following exposure to high temperature, drought, and salt treatments demonstrated that most of the MsGH3 family genes, especially MsGH3-12, MsGH3-13, and MsGH3-15, play important roles in stress responses. These findings provide invaluable insight for future practical analyses and genetic improvement of alfalfa abiotic stress tolerance.

Full-length transcriptional analysis reveals the complex relationship of leaves and roots in responses to cold-drought combined stress in common vetch.

Plant responses to single or combined abiotic stresses between aboveground and underground parts are complex and require crosstalk signaling pathways. In this study, we explored the transcriptome data of common vetch (Vicia sativa L.) subjected to cold and drought stress between leaves and roots via meta-analysis to identify the hub abiotic stress-responsive genes. A total of 4,836 and 3,103 differentially expressed genes (DEGs) were identified in the leaves and roots, respectively. Transcriptome analysis results showed that the set of stress-responsive DEGs to concurrent stress is distinct from single stress, indicating a specialized and unique response to combined stresses in common vetch. Gene Ontology (GO) enrichment analyses identified that "Photosystem II," "Defence response," and "Sucrose synthase/metabolic activity" were the most significantly enriched categories in leaves, roots, and both tissues, respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis results indicated that "ABC transporters" are the most enriched pathway and that all of the genes were upregulated in roots. Furthermore, 29 co-induced DEGs were identified as hub genes based on the consensus expression profile module of single and co-occurrence stress analysis. In transgenic yeast, the overexpression of three cross-stress tolerance candidate genes increased yeast tolerance to cold-drought combined stress. The elucidation of the combined stress-responsive network in common vetch to better parse the complex regulation of abiotic responses in plants facilitates more adequate legume forage breeding for combined stress tolerance.

Comparative Transcriptomic Analysis of Root and Leaf Transcript Profiles Reveals the Coordinated Mechanisms in Response to Salinity Stress in Common Vetch.

Owing to its strong environmental suitability to adverse abiotic stress conditions, common vetch (Vicia sativa) is grown worldwide for both forage and green manure purposes and is an important protein source for human consumption and livestock feed. The germination of common vetch seeds and growth of seedlings are severely affected by salinity stress, and the response of common vetch to salinity stress at the molecular level is still poorly understood. In this study, we report the first comparative transcriptomic analysis of the leaves and roots of common vetch under salinity stress. A total of 6361 differentially expressed genes were identified in leaves and roots. In the roots, the stress response was dominated by genes involved in peroxidase activity. However, the genes in leaves focused mainly on Ca2+ transport. Overexpression of six salinity-inducible transcription factors in yeast further confirmed their biological functions in the salinity stress response. Our study provides the most comprehensive transcriptomic analysis of common vetch leaf and root responses to salinity stress. Our findings broaden the knowledge of the common and distinct intrinsic molecular mechanisms within the leaves and roots of common vetch and could help to develop common vetch cultivars with high salinity tolerance.

Molecular Characterization of the miR156/MsSPL Model in Regulating the Compound Leaf Development and Abiotic Stress Response in Alfalfa.

Plant leaf patterns and shapes are spectacularly diverse. Changing the complexity of leaflet numbers is a valuable approach to increase its nutrition and photosynthesis. Alfalfa (Medicago sativa) is the most important forage legume species and has diversified compound leaf patterns, which makes it a model species for studying compound leaf development. However, transcriptomic information from alfalfa remains limited. In this study, RNA-Seq technology was used to identify 3746 differentially expressed genes (DEGs) between multifoliate and trifoliate alfalfa. Through an analysis of annotation information and expression data, SPL, one of the key regulators in modifiable plant development and abiotic stress response, was further analyzed. Here, thirty MsSPL genes were obtained from the alfalfa genome, of which 16 had the putative miR156 binding site. A tissue expression pattern analysis showed that the miR156-targeted MsSPLs were divided into two classes, namely, either tissue-specific or widely expressed in all tissues. All miR156-targeted SPLs strongly showed diversification and positive roles under drought and salt conditions. Importantly, miR156/MsSPL08 was significantly suppressed in multifoliate alfalfa. Furthermore, in the paralogous mutant of MsSPL08 isolated from Medicago truncatula, the phenotypes of mutant plants reveal that miR156/MsSPL08 is involved not only involved the branches but also especially regulates the number of leaflets. The legume is a typical compound leaf plant; the ratio of the leaflet often affects the quality of the forage. This study sheds light on new functions of SPL genes that regulate leaflet number development.

Molecular characterization of the COPT/Ctr-type copper transporter family under heavy metal stress in alfalfa.

In nature, heavy metals significantly affect crop growth and quality. Among various heavy metals, copper (Cu) is both essential and toxic to plants depending on the concentration and complex homeostatic networks. The Cu transporter family (COPT) plays important roles in Cu homeostasis, including absorption, transportation, and growth in plants; however, this gene family is still poorly understood in alfalfa (Medicago sativa L.). In this study, a total of 12 MsCOPTs were identified and characterized. Based on the conserved motif and phylogenetic analysis, MsCOPTs could be divided into four subgroups (A1, A2, A3, and B). Gene structure, chromosomal location, and synteny analyses of MsCOPTs showed that segmental and tandem duplications likely contributed to their evolution. Tissue-specific expression analysis of MsCOPT genes indicated diverse spatiotemporal expression patterns. Most MsCOPT genes had high transcription levels in roots and nodules, indicating that these genes may play vital roles in the absorption and transport of Cu through root. The complementary heterologous expression function of yeast once again indicates that root-specific COPT can supplement the growth of defective yeast strains on YPEG medium, suggesting that these genes are Cu transporters. In summary, for the first time, our research identified COPT family genes at the whole-genome level to provide guidance for effectively improving the problem of Cu deficiency in the grass-livestock chain and provide theoretical support for the subsequent development of grass and animal husbandry.

WRKY Transcription Factors in Medicago sativa L.: Genome-Wide Identification and Expression Analysis Under Abiotic Stress.

Alfalfa (Medicago sativa L.) is the most widely cultivated leguminous herb in the world. Its agricultural development has been restricted by various adverse environmental conditions, including water deficiency, high salinity, and low temperature. WRKY transcription factors (TFs) serve important roles in the regulation of plant development and stress responses. Research on the WRKY gene family has been reported for several species, but minimal information is available for alfalfa. In the present study, a total of 107 WRKY genes were identified in alfalfa and divided into 3 main groups. The classification, evolution, conserved motifs, and tissue expression were comprehensively analyzed. Meanwhile, 27 MsWRKY candidate genes that may be involved in abiotic stress were isolated through an analysis of gene expression profiles under different stresses, including cold, abscisic acid, drought, and salt treatments. Additionally, investigation of the cis-elements and potential biological functions of these genes further revealed that MsWRKY TFs may serve important roles in multiple stress resistance in alfalfa. This study provides an important foundation for future cloning and functional studies of WRKY genes in alfalfa.

Genome-wide development of miRNA-based SSR markers in Cleistogenes songorica and analysis of their transferability to Gramineae/non-Gramineae species.

Simple sequence repeat (SSR) markers are commonly used for many genetic applications, such as map construction, fingerprinting, and genetic diversity analyses, due to their high reproducibility, polymorphism, and abundance. Endogenous miRNAs play essential roles in plant development and gene expression under diverse biotic and abiotic stress conditions. In the present study, we predicted 110 miRNA-SSR primer pairs from 287 precursor miRNAs. Among 110 primer pairs, 85 were successfully amplified and examined for transferability to other Gramineae and non-Gramineae species. The results showed that all 82 primer pairs yielded unambiguous and strong amplification, and across the 23 studied Cleistogenes accessions, a total of 385 alleles were polymorphic. The number of alleles produced per primer varied from 3 to 11, with an average of 4.69 per locus. The expected heterozygosity (He) ranged from 0.44 to 0.88, with an average of 0.74 per locus, and the PIC (Polymorphism Information Content) values ranged from 0.34 to 0.87, with an average of 0.69 per locus. Furthermore, 1422 miRNA target genes were predicted and analyzed using the GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) databases. In conclusion, the results showed that an miRNA-based microsatellite marker system can be applicable for genetic diversity and marker-assisted breeding studies.

Coordinated mechanisms of leaves and roots in response to drought stress underlying full-length transcriptome profiling in Vicia sativa L.

BACKGROUND: Common vetch (Vicia sativa L.) is an important self-pollinating annual forage legume and is of interest for drought prone regions as a protein source to feed livestock and human consumption. However, the development and production of common vetch are negatively affected by drought stress. Plants have evolved common or distinct metabolic pathways between the aboveground and underground in response to drought stress. Little is known regarding the coordinated response of aboveground and underground tissues of common vetch to drought stress. RESULTS: Our results showed that a total of 30,427 full-length transcripts were identified in 12 samples, with an average length of 2278.89 bp. Global transcriptional profiles of the above 12 samples were then analysed via Illumina-Seq. A total of 3464 and 3062 differentially expressed genes (DEGs) were identified in the leaves and roots, respectively. Gene Ontology (GO) enrichment analyses identified that the dehydrin genes and Δ1-pyrroline-5-carboxylate synthase were induced for the biosynthesis of proline and water conservation. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis results indicated that the DEGs were significantly enriched in hormone signal transduction, starch and sucrose metabolism, and arginine and proline metabolism, and various drought response candidate genes were also identified. Abscisic acid (ABA; the AREB/ABF-SnRK2 pathway) regulates the activity of AMY3 and BAM1 to induce starch degradation in leaves and increase carbon export to roots, which may be associated with the drought stress responses in common vetch. Among the co-induced transcription factors (TFs), AREB/ABF, bHLH, MYB, WRKY, and AP2/ERF had divergent expression patterns and may be key in the crosstalk between leaves and roots during adaption to drought stress. In transgenic yeast, the overexpression of four TFs increased yeast tolerance to osmotic stresses. CONCLUSION: The multipronged approach identified in the leaves and roots broadens our understanding of the coordinated mechanisms of drought response in common vetch, and further provides targets to improve drought resistance through genetic engineering.

MYB transcription factors in alfalfa (Medicago sativa): genome-wide identification and expression analysis under abiotic stresses.

BACKGROUND: Alfalfa is the most widely cultivated forage legume and one of the most economically valuable crops in the world. Its survival and production are often hampered by environmental changes. However, there are few studies on stress-resistance genes in alfalfa because of its incomplete genomic information and rare expression profile data. The MYB proteins are characterized by a highly conserved DNA-binding domain, which is large, functionally diverse, and represented in all eukaryotes. The role of MYB proteins in plant development is essential; they function in diverse biological processes, including stress and defense responses, and seed and floral development. Studies on the MYB gene family have been reported in several species, but they have not been comprehensively analyzed in alfalfa. METHODS: To identify more comprehensive MYB transcription factor family genes, the sequences of 168 Arabidopsis thaliana, 430 Glycine max, 185 Medicago truncatula, and 130 Oryza sativa MYB proteins were downloaded from the Plant Transcription Factor Database. These sequences were used as queries in a BLAST search against the M. sativa proteome sequences provided by the Noble Research Institute. RESULTS: In the present study, a total of 265 MsMYB proteins were obtained, including 50 R1-MYB, 186 R2R3-MYB, 26 R1R2R3-MYB, and three atypical-MYB proteins. These predicted MsMYB proteins were divided into 12 subgroups by phylogenetic analysis, and gene ontology (GO) analysis indicated that most of the MsMYB genes are involved in various biological processes. The expression profiles and quantitative real-time PCR analysis indicated that some MsMYB genes might play a crucial role in the response to abiotic stresses. Additionally, a total of 170 and 914 predicted protein-protein and protein-DNA interactions were obtained, respectively. The interactions between MsMYB043 and MSAD320162, MsMYB253 and MSAD320162, and MsMYB253 and MSAD308489 were confirmed by a yeast two-hybrid system. This work provides information on the MYB family in alfalfa that was previously lacking and might promote the cultivation of stress-resistant alfalfa.

Genome-Wide Identification and Expression Profiling of the ERF Gene Family in Medicago sativa L. Under Various Abiotic Stresses.

The AP2/ERF (APETALA2/ETHYLENE RESPONSE FACTOR) transcription factor represents one of the largest plant-specific transcriptional regulators in plants. ERF plays important roles in the regulation of various developmental processes and acts as a mediator in plant external stress responses. However, the research of the ERF gene family is still limited in alfalfa (Medicago sativa L.), one of the most important forage legume species in the world. In the present study, a total of 159 ERF genes were identified, and the phylogenetic reconstruction, classification, conserved motifs, signal peptide prediction, and expression patterns under salt, drought, and low-temperature stresses of these ERF genes were comprehensively analyzed. The ERF genes family in alfalfa could be classified into 10 groups and predicted to be strongly homologous. Based on the structure and functions relationships, the III and IV subfamilies were more likely to play functions in abiotic stresses and 18 MsERF genes were selected for further quantitative real-time PCR validation in different stresses treatment. The results showed that all these MsERF genes were upregulated under three stresses except MsERF008. This study identified the possibility of abiotic tolerance candidate genes playing various roles in stress resistance at the whole-genome level, which would provide primary understanding for exploring ERF-mediated tolerance in alfalfa.

Transcriptome-Wide Characterization and Functional Identification of the Aquaporin Gene Family During Drought Stress in Common Vetch.

Aquaporins (AQPs) are transmembrane channels that are essential for the movement of water and other small molecules between biofilms in various physiological processes in plants. In this study, based on transcriptome-wide data, we identified and described a total of 21 AQP genes in common vetch (Vicia sativa subsp. sativa), which is an economically important pasture legume worldwide. Based on phylogenetic analyses, the VsAQPs were sorted into four subfamilies, including four plasma membrane intrinsic proteins (PIPs), six tonoplast intrinsic proteins (TIPs), seven NOD26-like intrinsic proteins, and four small basic intrinsic proteins. Furthermore, chemical and physical properties of these VsAQPs, including the isoelectric point and theoretical molecular weight, were analyzed. Analyses of the AQP signature sequences and key residues indicated the substrate specificity of each VsAQP. A set of VsAQPs was selected for gene expression analysis in a number of tissues and after drought stress treatments using real-time quantitative reverse transcription/polymerase chain reaction assays. Most of the PIPs and TIPs were proposed to have critical roles in regulating the flow of water during drought stress. Heterologous expression experiments in yeast indicated that VsPIP1;2 and VsPIP2;2 are key candidate genes for improving drought stress tolerance. The results reported in this study could be a crucial resource for further practical analyses and for genetic improvement of drought stress tolerance in common vetch.

Transcriptome-wide characterization and functional analysis of MATE transporters in response to aluminum toxicity in Medicago sativa L.

Multidrug and toxic compound extrusion (MATE) transporters contribute to multidrug resistance and play major determinants of aluminum (Al) tolerance in plants. Alfalfa (Medicago sativa L.) is the most extensively cultivated forage crop in the world, yet most alfalfa cultivars are not Al tolerant. The basic knowledge of the MATE transcripts family and the characterisation of specific MATE members involved in alfalfa Al stress remain unclear. In this study, 88 alfalfa MATE (MsMATE) transporters were identified at the whole transcriptome level. Phylogenetic analysis classified them into four subfamilies comprising 11 subgroups. Generally, five kinds of motifs were found in group G1, and most were located at the N-terminus, which might confer these genes with Al detoxification functions. Furthermore, 10 putative Al detoxification-related MsMATE genes were identified and the expression of five genes was significantly increased after Al treatment, indicating that these genes might play important roles in conferring Al tolerance to alfalfa. Considering the limited functional understanding of MATE transcripts in alfalfa, our findings will be valuable for the functional investigation and application of this family in alfalfa.

Mining Late Embryogenesis Abundant (LEA) Family Genes in Cleistogenes songorica, a Xerophyte Perennial Desert Plant.

Plant growth and development depends on its ability to maintain optimal cellular homeostasis during abiotic and biotic stresses. Cleistogenes songorica, a xerophyte desert plant, is known to have novel drought stress adaptation strategies and contains rich pools of stress tolerance genes. Proteins encoded by Late Embryogenesis Abundant (LEA) family genes promote cellular activities by functioning as disordered molecules, or by limiting collisions between enzymes during stresses. To date, functions of the LEA family genes have been heavily investigated in many plant species except perennial monocotyledonous species. In this study, 44 putative LEA genes were identified in the C. songorica genome and were grouped into eight subfamilies, based on their conserved protein domains and domain organizations. Phylogenetic analyses indicated that C. songorica Dehydrin and LEA_2 subfamily proteins shared high sequence homology with stress responsive Dehydrin proteins from Arabidopsis. Additionally, promoter regions of CsLEA_2 or CsDehydrin subfamily genes were rich in G-box, drought responsive (MBS), and/or Abscisic acid responsive (ABRE) cis-regulatory elements. In addition, gene expression analyses indicated that genes from these two subfamilies were highly responsive to heat stress and ABA treatment, in both leaves and roots. In summary, the results from this study provided a comprehensive view of C. songorica LEA genes and the potential applications of these genes for the improvement of crop tolerance to abiotic stresses.

Genome-wide identification and expression analysis of the fatty acid desaturase genes in Medicago truncatula.

Fatty acid desaturases (FADs) are of great importance and play critical roles in regulating plant fatty acid (FA) compositions. But to date, no reports about characterization of the FAD genes have been reported in the model dicotyledonous grass species Medicago truncatula. In this study, using database searches, 20 full-length FAD genes were identified in M. truncatula. These FAD genes were unevenly distributed on six chromosomes except the chromosome 6 and 8. Phylogenetic analysis showed the FAD genes in M. truncatula were clustered into six subfamilies and had similar exon number and intron phase in the same subfamily. Moreover, expression analysis based on qRT-PCR indicated these FAD genes were extensively involved in cold and heat responses. This study would provide an important foundation for future cloning and functional studies of FAD genes in M. truncatula and other related legume species.

Genome-Wide Development of MicroRNA-Based SSR Markers in Medicago truncatula with Their Transferability Analysis and Utilization in Related Legume Species.

Microsatellite (simple sequence repeats, SSRs) marker is one of the most widely used markers in marker-assisted breeding. As one type of functional markers, MicroRNA-based SSR (miRNA-SSR) markers have been exploited mainly in animals, but the development and characterization of miRNA-SSR markers in plants are still limited. In the present study, miRNA-SSR markers for Medicago truncatula (M. truncatula) were developed and their cross-species transferability in six leguminous species was evaluated. A total of 169 primer pairs were successfully designed from 130 M. truncatula miRNA genes, the majority of which were mononucleotide repeats (70.41%), followed by dinucleotide repeats (14.20%), compound repeats (11.24%) and trinucleotide repeats (4.14%). Functional classification of SSR-containing miRNA genes showed that all targets could be grouped into three Gene Ontology (GO) categories: 17 in biological process, 11 in molecular function, and 14 in cellular component. The miRNA-SSR markers showed high transferability in other six leguminous species, ranged from 74.56% to 90.53%. Furthermore, 25 Mt-miRNA-SSR markers were used to evaluate polymorphisms in 20 alfalfa accessions, and the polymorphism information content (PIC) values ranged from 0.39 to 0.89 with an average of 0.71, the allele number per marker varied from 3 to 18 with an average of 7.88, indicating a high level of informativeness. The present study is the first time developed and characterized of M. truncatula miRNA-SSRs and demonstrated their utility in transferability, these novel markers will be valuable for genetic diversity analysis, marker-assisted selection and genotyping in leguminous species.