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

Genome-wide analysis of the R2R3-MYB gene family and identification of candidate genes that regulate isoflavone biosynthesis in red clover (Trifolium pratense)

Red clover (Trifolium pratense) is a perennial legume with high feeding and medicinal value attributed to its abundant isoflavone content. Previous studies reported that R2R3-MYB transcription factors are involved in the biosynthesis of isoflavones; however, their specific role in red clover remains poorly understood. Through comprehensive genome-wide and transcriptome analyses, a total of 138 TpR2R3-MYB genes were identified and classified into 30 distinct subgroups within a phylogenetic tree. Importantly, six of these subgroups showed associations with isoflavone biosynthesis in red clover. The majority of segmental duplication events (Ka/Ks < 1) indicated that the TpR2R3-MYB gene underwent strong purifying selection during evolution. The qRT-PCR analysis demonstrated high expression levels of TpMYB79 and TpMYB53 in Minshan red clover at full flowering stage, consistent with the trend for isoflavone content determination, suggesting that TpMYB79 and TpMYB53 might be important regulators of isoflavone biosynthesis in red clover. Additionally, we observed nucleus and vacuole membrane localization of TpMYB53 and TpMYB79, with TpMYB53 primarily exerting transcriptional activation through its C-terminal activation motifs while TpMYB79 exhibited no transcriptional activity. These findings provided a foundation for the study of the biological function of R2R3-MYB transcription factors in red clover.

MYB transcription factors, their regulation and interactions with non-coding RNAs during drought stress in Brassica juncea

BACKGROUND

Brassica juncea (L.) Czern is an important oilseed crop affected by various abiotic stresses like drought, heat, and salt. These stresses have detrimental effects on the crop's overall growth, development and yield. Various Transcription factors (TFs) are involved in regulation of plant stress response by modulating expression of stress-responsive genes. The myeloblastosis (MYB) TFs is one of the largest families of TFs associated with various developmental and biological processes such as plant growth, secondary metabolism, stress response etc. However, MYB TFs and their regulation by non-coding RNAs (ncRNAs) in response to stress have not been studied in B. juncea. Thus, we performed a detailed study on the MYB TF family and their interactions with miRNAs and Long non coding RNAs.

RESULTS

Computational investigation of genome and proteome data presented a comprehensive picture of the MYB genes and their protein architecture, including intron-exon organisation, conserved motif analysis, R2R3 MYB DNA-binding domains analysis, sub-cellular localization, protein-protein interaction and chromosomal locations. Phylogenetically, BjuMYBs were further classified into different subclades on the basis of topology and classification in Arabidopsis. A total of 751 MYBs were identified in B. juncea corresponding to 297 1R-BjuMYBs, 440 R2R3-BjuMYBs, 12 3R-BjuMYBs, and 2 4R-BjuMYBs types. We validated the transcriptional profiles of nine selected BjuMYBs under drought stress through RT-qPCR. Promoter analysis indicated the presence of drought-responsive cis-regulatory components. Additionally, the miRNA-MYB TF interactions was also studied, and most of the microRNAs (miRNAs) that target BjuMYBs were involved in abiotic stress response and developmental processes. Regulatory network analysis and expression patterns of lncRNA-miRNA-MYB indicated that selected long non-coding RNAs (lncRNAs) acted as strong endogenous target mimics (eTMs) of the miRNAs regulated expression of BjuMYBs under drought stress.

CONCLUSIONS

The present study has established preliminary groundwork of MYB TFs and their interaction with ncRNAs in B. juncea and it will help in developing drought- tolerant Brassica crops.

Genome-wide identification of JAZ gene family members in autotetraploid cultivated alfalfa (Medicago sativa subsp. sativa) and expression analysis under salt stress

BACKGROUND

Jasmonate ZIM-domain (JAZ) proteins, which act as negative regulators in the jasmonic acid (JA) signalling pathway, have significant implications for plant development and response to abiotic stress.

RESULTS

Through a comprehensive genome-wide analysis, a total of 20 members of the JAZ gene family specific to alfalfa were identified in its genome. Phylogenetic analysis divided these 20 MsJAZ genes into five subgroups. Gene structure analysis, protein motif analysis, and 3D protein structure analysis revealed that alfalfa JAZ genes in the same evolutionary branch share similar exon‒intron, motif, and 3D structure compositions. Eight segmental duplication events were identified among these 20 MsJAZ genes through collinearity analysis. Among the 32 chromosomes of the autotetraploid cultivated alfalfa, there were 20 MsJAZ genes distributed on 17 chromosomes. Extensive stress-related cis-acting elements were detected in the upstream sequences of MsJAZ genes, suggesting that their response to stress has an underlying function. Furthermore, the expression levels of MsJAZ genes were examined across various tissues and under the influence of salt stress conditions, revealing tissue-specific expression and regulation by salt stress. Through RT‒qPCR experiments, it was discovered that the relative expression levels of these six MsJAZ genes increased under salt stress.

CONCLUSIONS

In summary, our study represents the first comprehensive identification and analysis of the JAZ gene family in alfalfa. These results provide important information for exploring the mechanism of JAZ genes in alfalfa salt tolerance and identifying candidate genes for improving the salt tolerance of autotetraploid cultivated alfalfa via genetic engineering in the future.

MODMS: a multi-omics database for facilitating biological studies on alfalfa (Medicago sativa L.)

Alfalfa (Medicago sativa L.) is a globally important forage crop. It also serves as a vegetable and medicinal herb because of its excellent nutritional quality and significant economic value. Multi-omics data on alfalfa continue to accumulate owing to recent advances in high-throughput techniques, and integrating this information holds great potential for expediting genetic research and facilitating advances in alfalfa agronomic traits. Therefore, we developed a comprehensive database named MODMS (multi-omics database of M. sativa) that incorporates multiple reference genomes, annotations, comparative genomics, transcriptomes, high-quality genomic variants, proteomics, and metabolomics. This report describes our continuously evolving database, which provides researchers with several convenient tools and extensive omics data resources, facilitating the expansion of alfalfa research. Further details regarding the MODMS database are available at https://modms.lzu.edu.cn/.

Characterization of Dendrobium catenatum CBL-CIPK signaling networks and their response to abiotic stress

Dendrobium catenatum is a traditional Chinese medicine listing as rare and endangered due to environmental impacts. But little is known about its stress resistance mechanism. The CBL-CIPK signaling pathway played vital roles in various stress responses. In this study, we identified 9 calcineurin B-like (CBL) genes and 28 CBL-interacting protein kinase (CIPK) genes from D. catenatum. Phylogenetic analysis showed that DcCBL and DcCIPK families could be divided into four and six subgroups, respectively. Members in each subgroup had similar gene structures. Cis-acting element analyses showed that these genes were involved in stress responses and hormone signaling. Spatial expression profiles showed that they were tissue-specific, and expressed lower in vegetative organs than reproductive organs. Gene expression analyses revealed that these genes were involved in drought, heat, cold, and salt responses and depended on abscisic acid (ABA) and salicylic acid (SA) signaling pathways. Furthermore, we cloned 19 DcCIPK genes and 9 DcCBL genes and detected ten interacting CBL-CIPK combinations using yeast two-hybrid system. Finally, we constructed 20 CBL-CIPK signaling pathways based on their expression patterns and interaction relationships. These results established CBL-CIPK signaling pathway responding to abiotic stress and provided a molecular basis for improving D. catenatum stress resistance in the future.

Genome-wide identification and analysis of the NLR gene family in Medicago ruthenica

Medicago ruthenica, important forage in the legume family, possesses high nutritional value and carries abundant tolerance genes. This study used whole-genome data of M. ruthenica to perform a genome-wide analysis of the nucleotide-binding site-leucine-rich repeat receptor (NLR) gene family, which is the largest family of plant disease resistance genes (R genes). A total of 338 NLR genes were identified in the M. ruthenica genome, including 160 typical genes that contained 80 coiled-coil (CC)-NBS-LRR (CNL) genes, 76 toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) genes, four resistance to powdery mildew 8 (RPW8)-NBS-LRR (RNL) subclass genes, and 178 atypical NLR genes encoding proteins without at least one important domain. Among its eight chromosomes, M. ruthenica chromosomes 3 and 8 contained most of the NLR genes. More than 40% of all NLR genes were located on these two chromosomes, mainly in multigene clusters. The NLR proteins of M. ruthenica had six highly conserved motifs: P-loop, GLPL, RNBS-D, kinase-2, RNBS-C, and MHDV. Phylogenetic analysis revealed that the NLR genes of M. ruthenica formed three deeply separated clades according to the N-terminal domain of the proteins encoded by these genes. Gene duplication and syntenic analysis suggested four gene duplication types in the NLR genes of M. ruthenica, namely, tandem, proximal, dispersed, and segmental duplicates, which involved 189, 49, 59, and 41 genes, respectively. A total of 41 segmental duplication genes formed 23 NLR gene pairs located on syntenic chromosomal blocks mainly between chromosomes 6 and 7. In addition, syntenic analysis between M. truncatula and M. ruthenica revealed 193 gene pairs located on syntenic chromosomal blocks of the two species. The expression analysis of M. ruthenica NLR genes showed that 303 (89.6%) of the NLR genes were expressed in different varieties. Overall, this study described the full NLR profile of the M. ruthenica genome to provide an important resource for mining disease-resistant genes and disease-resistant breeding.

Genome-Wide Identification of ERF Transcription Factor Family and Functional Analysis of the Drought Stress-Responsive Genes in Melilotus albus

As an important forage legume with high values in feed and medicine, Melilotus albus has been widely cultivated. The AP2/ERF transcription factor has been shown to play an important regulatory role in plant drought resistance, but it has not been reported in the legume forage crop M. albus. To digger the genes of M. albus in response to drought stress, we identified and analyzed the ERF gene family of M. albus at the genome-wide level. A total of 100 MaERF genes containing a single AP2 domain sequence were identified in this study, named MaERF001 to MaERF100, and bioinformatics analysis was performed. Collinearity analysis indicated that segmental duplication may play a key role in the expansion of the M. albus ERF gene family. Cis-acting element predictions suggest that MaERF genes are involved in various hormonal responses and abiotic stresses. The expression patterns indicated that MaERFs responded to drought stress to varying degrees. Furthermore, four up-regulated ERFs (MaERF008, MaERF037, MaERF054 and MaERF058) under drought stress were overexpressed in yeast and indicated their biological functions to confer the tolerance to drought. This work will advance the understanding of the molecular mechanisms underlying the drought response in M. albus. Further study of the promising potential candidate genes identified in this study will provide a valuable resource as the next step in functional genomics studies and improve the possibility of improving drought tolerance in M. albus by transgenic approaches.

The transcription factors VaERF16 and VaMYB306 interact to enhance resistance of grapevine to Botrytis cinerea infection

Botrytis cinerea is a fungus that infects cultivated grape (Vitis vinifera); the identification and characterization of resistance mechanisms in the host is of great importance for the grape industry. Here, we report that a transcription factor in the ethylene-responsive factor (ERF) family (VaERF16) from Chinese wild grape (Vitis amurensis 'Shuang You') is expressed during B. cinerea infection and in response to treatments with the hormones ethylene and methyl jasmonate. Heterologous overexpression of VaERF16 in Arabidopsis thaliana substantially enhanced resistance to B. cinerea and the bacterium Pseudomonas syringae DC3000 via the salicylic acid and jasmonate/ethylene signalling pathways. Yeast two-hybrid, bimolecular fluorescence complementation, and co-immunoprecipitation assays indicated that VaERF16 interacts with the MYB family transcription factor VaMYB306. Overexpression of VaERF16 or VaMYB306 in grape leaves increased resistance to B. cinerea and caused an up-regulation of the defence-related gene PDF1.2, which encodes a defensin-like protein. Conversely, silencing of either gene resulted in increased susceptibility to B. cinerea. Yeast one-hybrid and dual-luciferase assays indicated that VaERF16 increased the transcript levels of VaPDF1.2 by binding directly to the GCC box in its promoter. Notably, VaMYB306 alone did not bind to the VaPDF1.2 promoter, but the VaERF16-VaMYB306 transcriptional complex resulted in higher transcript levels of VaPDF1.2, suggesting that the proteins function through their mutual interaction. Elucidation of this regulatory module may be of value in enhancing resistance of grapevine to B. cinerea infection.

Integrative Analyses of Transcriptomes and Metabolomes Reveal Associated Genes and Metabolites with Flowering Regulation in Common Vetch (Vicia sativa L.)

As an important source of protein for livestock and human consumption, Vicia sativa is cultivated worldwide, but its seed production is hampered at high altitudes because of the short frost-free period. Flowering represents the transition from a vegetative to a reproductive period, and early flowering benefits plant seed production at high altitudes. However, the molecular mechanisms of flowering regulation in V. sativa remain elusive. In the present study, two V. sativa accessions with different flowering characteristics were used: Lan3 (early-flowering) was cultivated by our laboratory, and 503 (late-flowering) was selected from 222 V. sativa accessions after three years of field experiments. The shoot samples (shoot tip length = 10 cm) of these two accessions were collected 63, 70, and 77 days after sowing, and the molecular regulatory mechanism of the flowering process was identified by integrative analyses of the transcriptomes and metabolomes. Kyoto Encyclopedia of Genes and Genomes enrichment showed that the synthesis and signal transduction of plant hormone pathways were the most enriched pathways in 4274 differentially expressed genes (DEGs) and in 259 differential metabolites between Lan3 and 503. Moreover, the contents of three metabolites related to salicylic acid biosynthesis and the transcription levels of two DEGs related to salicylic acid signal transduction in Lan3 were higher than those in 503. Further verification in various accessions indicated that salicylic acid metabolism may be involved in the flowering regulation process of V. sativa. These findings provide valuable information for understanding the flowering mechanism and for promoting breeding research in V. sativa.

Evolution and functional diversification of R2R3-MYB transcription factors in plants

R2R3-MYB genes (R2R3-MYBs) form one of the largest transcription factor gene families in the plant kingdom, with substantial structural and functional diversity. However, the evolutionary processes leading to this amazing functional diversity have not yet been clearly established. Recently developed genomic and classical molecular technologies have provided detailed insights into the evolutionary relationships and functions of plant R2R3-MYBs. Here, we review recent genome-level and functional analyses of plant R2R3-MYBs, with an emphasis on their evolution and functional diversification. In land plants, this gene family underwent a large expansion by whole genome duplications and small-scale duplications. Along with this population explosion, a series of functionally conserved or lineage-specific subfamilies/groups arose with roles in three major plant-specific biological processes: development and cell differentiation, specialized metabolism, and biotic and abiotic stresses. The rapid expansion and functional diversification of plant R2R3-MYBs are highly consistent with the increasing complexity of angiosperms. In particular, recently derived R2R3-MYBs with three highly homologous intron patterns (a, b, and c) are disproportionately related to specialized metabolism and have become the predominant subfamilies in land plant genomes. The evolution of plant R2R3-MYBs is an active area of research, and further studies are expected to improve our understanding of the evolution and functional diversification of this gene family.

Genome-Wide Identification of the Q-type C2H2 Transcription Factor Family in Alfalfa (Medicago sativa) and Expression Analysis under Different Abiotic Stresses

Q-type C2H2 zinc-finger protein (C2H2-ZFP) transcription factors are associated with many plant growth development and environmental stress responses. To date, there have been few analyses of the Q-type C2H2-ZFP gene family in alfalfa (Medicago sativa subsp. sativa). In this study, we identified 58 Q-type C2H2-ZFPs across the entire alfalfa genome, and the gene structure, motif composition, chromosomal mapping, and cis-regulatory elements were explored, as well as the expression profiles of specific tissues and the response under different abiotic stresses. According to their phylogenetic features, these 58 MsZFPs were divided into 12 subgroups. Synteny analysis showed that duplication events play a vital role in the expansion of the MsZFP gene family. The collinearity results showed that a total of 26 and 42 of the 58 MsZFP genes were homologous with Arabidopsis and M. truncatula, respectively. The expression profiles showed that C2H2-ZFP genes played various roles in different tissues and abiotic stresses. The results of subsequent quantitative real-time polymerase chain reaction (qRT-PCR) showed that the nine selected MsZFP genes were rapidly induced under different abiotic stresses, indicating that C2H2-ZFP genes are closely related to abiotic stress. This study provides results on MsZFP genes, their response to various abiotic stresses, and new information on the C2H2 family in alfalfa.

Genome-Wide Identification and Expression Analysis of MYB Transcription Factor Superfamily in Dendrobium catenatum

Dendrobium catenatum is an important traditional Chinese medicine and naturally grows on tree trunks and cliffs, where it can encounter diverse environmental stimuli. MYB transcription factors are widely involved in response to abiotic stresses. However, the MYB gene family has not yet been systematically cataloged in D. catenatum. In this study, a total of 133 MYB proteins were identified in D. catenatum, including 32 MYB-related, 99 R2R3-MYB, 1 3R-MYB, and 1 4R-MYB proteins. Phylogenetic relationships, conserved motifs, gene structures, and expression profiles in response to abiotic stresses were then analyzed. Phylogenetic analysis revealed MYB proteins in D. catenatum could be divided into 14 subgroups, which was supported by the conserved motif compositions and gene structures. Differential DcMYB gene expression and specific responses were analyzed under drought, heat, cold, and salt stresses using RNA-seq and validated by qRT-PCR. Forty-two MYB genes were differentially screened following exposure to abiotic stresses. Five, 12, 11, and 14 genes were specifically expressed in response to drought, heat, cold, and salt stress, respectively. This study identified candidate MYB genes with possible roles in abiotic tolerance and established a theoretical foundation for molecular breeding of D. catenatum.

Pan-transcriptome identifying master genes and regulation network in response to drought and salt stresses in Alfalfa (Medicago sativa L.)

Alfalfa is an important legume forage grown worldwide and its productivity is affected by environmental stresses such as drought and high salinity. In this work, three alfalfa germplasms with contrasting tolerances to drought and high salinity were used for unraveling the transcriptomic responses to drought and salt stresses. Twenty-one different RNA samples from different germplasm, stress conditions or tissue sources (leaf, stem and root) were extracted and sequenced using the PacBio (Iso-Seq) and the Illumina platforms to obtain full-length transcriptomic profiles. A total of 1,124,275 and 91,378 unique isoforms and genes were obtained, respectively. Comparative analysis of transcriptomes identified differentially expressed genes and isoforms as well as transcriptional and post-transcriptional modifications such as alternative splicing events, fusion genes and nonsense-mediated mRNA decay events and non-coding RNA such as circRNA and lncRNA. This is the first time to identify the diversity of circRNA and lncRNA in response to drought and high salinity in alfalfa. The analysis of weighted gene co-expression network allowed to identify master genes and isoforms that may play important roles on drought and salt stress tolerance in alfalfa. This work provides insight for understanding the mechanisms by which drought and salt stresses affect alfalfa growth at the whole genome level.

Genome-wide identification of the MADS-box transcription factor family in autotetraploid cultivated alfalfa (Medicago sativa L.) and expression analysis under abiotic stress

Background

Alfalfa, the “queen of forage”, is the most extensively cultivated forage legume in the world. The development and yield of alfalfa are seriously limited by abiotic stress. MADS-box transcription factors are one of the largest gene families and play a pivotal role in plant development and abiotic stress. However, little is known regarding the MADS-box transcription factors in autotetraploid cultivated alfalfa.

Results

In the present study, we identified 120 MsMADS-box genes in the alfalfa genome. Phylogenetic analysis indicated that 75 type-I MsMADS-box genes were classified into the Mα, Mβ, and Mγ subgroups, and 45 type-II MsMADS-box genes were classified into 11 subgroups. The promoter region of MsMADS-box genes containing several hormone and stress related elements. Chromosomal location analysis revealed that 117 MsMADS-box genes were unevenly distributed on 32 chromosomes, and the remaining three genes were located on unmapped scaffolds. A total of nine pairs of segmental duplications and four groups of tandem duplications were found. Expression analysis showed that MsMADS-box genes were differentially expressed in various tissues and under abiotic stresses. qRT-PCR analysis revealed that the expression profiles of eight selected MsMADS-box genes were distinct under various stresses.

Conclusions

In this study, MsMADS-box genes were identified in the cultivated alfalfa genome based on autotetraploid level, and further confirmed by Gene Ontology (GO) analysis, phylogenetic analysis, sequence features and expression analysis. Taken together, these findings will provide clues for further study of MsMADS-box functions and alfalfa molecular breeding.

Our study is the first to systematically identify and characterize the MADS-box transcription factors in autotetraploid cultivated alfalfa (Medicago sativa L.), and eight MsMADS-box genes were significantly involved in response to various stresses.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07911-9.

Genome-wide analysis of MYB transcription factors of Vaccinium corymbosum and their positive responses to drought stress

Background

Blueberry (Vaccinium corymbosum L.) is an important species with a high content of flavonoids in fruits. As a perennial shrub, blueberry is characterized by shallow-rooted property and susceptible to drought stress. MYB transcription factor was reported to be widely involved in plant response to abiotic stresses, however, the role of MYB family in blueberry responding to drought stress remains elusive.

Results

In this study, we conducted a comprehensive analysis of VcMYBs in blueberry based on the genome data under drought stress, including phylogenetic relationship, identification of differentially expressed genes (DEGs), expression profiling, conserved motifs, expression correlation and protein-protein interaction prediction, etc. The results showed that 229 non-redundant MYB sequences were identified in the blueberry genome, and divided into 23 subgroups. A total of 102 MYB DEGs with a significant response to drought stress were identified, of which 72 in leaves and 69 in roots, and 8 differential expression genes with a > 20-fold change in the level of expression. 17 DEGs had a higher expression correlation with other MYB members. The interaction partners of the key VcMYB proteins were predicted by STRING analysis and in combination with physiological and morphological observation. 10 key VcMYB genes such as VcMYB8, VcMYB102 and VcMYB228 were predicted to be probably involved in reactive oxygen species (ROS) pathway, and 7 key VcMYB genes (VcMYB41, VcMYB88 and VcMYB100, etc..) probably participated in leaf regulation under drought treatment.

Conclusions

Our studies provide a new understanding of the regulation mechanism of VcMYB family in blueberry response to drought stress, and lay fundamental support for future studies on blueberry grown in regions with limited water supply for this crop.

Supplementary Information

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

Characterization of the Liriodendron Chinense MYB Gene Family and Its Role in Abiotic Stress Response

Liriodendron chinense (Lchi) is a Magnoliaceae plant, which is a basic angiosperm left behind by the Pleistocene and mainly distributed in the south of the Yangtze River. Liriodendron hybrids has good wood properties and is widely used in furniture and in other fields. It is not clear if they can adapt to different environmental conditions, such as drought and high and low temperatures, and the molecular mechanisms for this adaptation are unknown. Among plant transcription factors (TFs), the MYB gene family is one of the largest and is often involved in stress or adversity response signaling, growth, and development. Therefore, studying the role of MYBTFs in regulating abiotic stress signaling, growth, and development in Lchi is helpful to promote afforestation in different environments. In our research, a genome-wide analysis of the LchiMYB gene family was performed, including the phylogenetic relationship tree, gene exon-intron structure, collinearity, and chromosomal position. According to the evolutionary tree, 190 LchiMYBs were divided into three main branches. LchiMYBs were evenly distributed across 19 chromosomes, with their collinearity, suggesting that segment duplication events may have contributed to LchiMYB gene expansion. Transcriptomes from eight tissues, 11 stages of somatic embryogenesis, and leaves after cold, heat, and drought stress were used to analyze the function of the MYB gene family. The results of tissue expression analysis showed that most LchiMYB genes regulated bark, leaf, bud, sepal, stigma, and stamen development, as well as the four important stages (ES3, ES4, ES9, and PL) of somatic embryogenesis. More than 60 LchiMYBs responded to heat, cold, and drought stress; some of which underwent gene duplication during evolution. LchiMYB3 was highly expressed under all three forms of stress, while LchiMYB121 was strongly induced by both cold and heat stress. Eight genes with different expression patterns were selected and verified by quantitative real-time PCR (qRT-PCR) experiments. The results suggested that these LchiMYBs may regulate Lchi growth development and resistance to abiotic stress. This study shows the cross-regulatory function of LchiMYBs in the growth and development, asexual reproduction, and abiotic resistance of Lchi. This information will prove pivotal to directing further studies on the biological function of Lchi MYBTFs in genetic improvement and abiotic stress response.

A chromosome-scale genome assembly of a diploid alfalfa, the progenitor of autotetraploid alfalfa

Alfalfa (Medicago sativa L.) is one of the most important and widely cultivated forage crops. It is commonly used as a vegetable and medicinal herb because of its excellent nutritional quality and significant economic value. Based on Illumina, Nanopore and Hi-C data, we assembled a chromosome-scale assembly of Medicago sativa spp. caerulea (voucher PI464715), the direct diploid progenitor of autotetraploid alfalfa. The assembled genome comprises 793.2 Mb of genomic sequence and 47,202 annotated protein-coding genes. The contig N50 length is 3.86 Mb. This genome is almost twofold larger and contains more annotated protein-coding genes than that of its close relative, Medicago truncatula (420 Mb and 44,623 genes). The more expanded gene families compared with those in M. truncatula and the expansion of repetitive elements rather than whole-genome duplication (i.e., the two species share the ancestral Papilionoideae whole-genome duplication event) may have contributed to the large genome size of M. sativa spp. caerulea. Comparative and evolutionary analyses revealed that M. sativa spp. caerulea diverged from M. truncatula ~5.2 million years ago, and the chromosomal fissions and fusions detected between the two genomes occurred during the divergence of the two species. In addition, we identified 489 resistance (R) genes and 82 and 85 candidate genes involved in the lignin and cellulose biosynthesis pathways, respectively. The near-complete and accurate diploid alfalfa reference genome obtained herein serves as an important complement to the recently assembled autotetraploid alfalfa genome and will provide valuable genomic resources for investigating the genomic architecture of autotetraploid alfalfa as well as for improving breeding strategies in alfalfa.

Genome-wide identification of R2R3-MYB family in wheat and functional characteristics of the abiotic stress responsive gene TaMYB344

Background

MYB superfamily is one of the most abundant families in plants, which plays important roles in plant growth, development, and productivity. However, to date, researches on MYBs in wheat (Triticum aestivum L.) are scattered mostly, not comprehensive.

Results

In this study, a total of 393 R2R3-MYBs and 12 R1R2R3-MYBs were identified and analyzed including gene structure, chromosomal distribution, synteny relationship, and evolutionary relationship. Then, 29 clusters tandem duplication and 8 clusters segmental duplication genes were discovered. The expression profile of the identified genes under abiotic and biotic stress was analyzed using RNA-seq data. Based on expression patterns analysis, we screened many candidate genes involved in plant response to abiotic and biotic stress. Among them, the functional characteristics of TaMYB344 were further studied. TaMYB344 was localized in the nucleus and functioned as a weak transcriptional activator. We demonstrated that TaMYB344-overexpressing transgenic tobacco plants had enhanced tolerance to drought, heat, and high salt stress.

Conclusions

In this study, 393 R2R3-MYBs and 12 R1R2R3-MYBs in wheat were systemically identified and analyzed. Differential expression analysis indicated that many R2R3-MYBs were involved in abiotic and biotic stress response. We identified a potential candidate gene TaMYB344, overexpression of which in tobacco plants enhanced drought, heat, and salt stress tolerance. These results will provide abundant molecular data for breeding new varieties of wheat in the future.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s12864-020-07175-9.

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

The Dof transcription factor is a plant-specific transcriptional regulator that plays important roles in plant development and acts as a mediator in plant external stress responses. However, Dofs have previously been identified in several plants but not in alfalfa (Medicago sativa L.), one of the most widely cultivated forage legumes. In the present study, a total of 40 MsDof genes were identified, and the phylogenetic reconstruction, classification, conserved motifs, and expression patterns under abscisic acid (ABA), cold, heat, drought and salt stresses of these Dof genes were comprehensively analyzed. The Dof genes family in alfalfa could be classified into eight classes. Gene ontology (GO) and tissue-specific analysis indicated that most MsDof genes may be involved in biological functions during plant growth. Moreover, the expression profiles and quantitative real-time PCR analysis indicated that eight candidate abiotic tolerance genes were induced in response to four abiotic stresses. This study identified the possibility of abiotic tolerance candidate genes playing various roles in stress resistance at the whole genome level, which would provide new information on the Dof family in alfalfa.

Genome-wide characterization of the abscisic acid-, stress- and ripening-induced (ASR) gene family in wheat (Triticum aestivum L.)

BACKGROUND

Abscisic acid-, stress-, and ripening-induced (ASR) genes are a class of plant specific transcription factors (TFs), which play important roles in plant development, growth and abiotic stress responses. The wheat ASRs have not been described in genome-wide yet.

METHODS

We predicted the transmembrane regions and subcellular localization using the TMHMM server, and Plant-mPLoc server and CELLO v2.5, respectively. Then the phylogeny tree was built by MEGA7. The exon-intron structures, conserved motifs and TFs binding sites were analyzed by GSDS, MEME program and PlantRegMap, respectively.

RESULTS

In wheat, 33ASR genes were identified through a genome-wide survey and classified into six groups. Phylogenetic analyses revealed that the TaASR proteins in the same group tightly clustered together, compared with those from other species. Duplication analysis indicated that the TaASR gene family has expanded mainly through tandem and segmental duplication events. Similar gene structures and conserved protein motifs of TaASRs in wheat were identified in the same groups. ASR genes contained various TF binding cites associated with the stress responses in the promoter region. Gene expression was generally associated with the expected group-specific expression pattern in five tissues, including grain, leaf, root, spike and stem, indicating the broad conservation of ASR genes function during wheat evolution. The qRT-PCR analysis revealed that several ASRs were up-regulated in response to NaCl and PEG stress.

CONCLUSION

We identified ASR genes in wheat and found that gene duplication events are the main driving force for ASR gene evolution in wheat. The expression of wheat ASR genes was modulated in responses to multiple abiotic stresses, including drought/osmotic and salt stress. The results provided important information for further identifications of the functions of wheat ASR genes and candidate genes for high abiotic stress tolerant wheat breeding.