Identification and analysis of flavonoid pathway genes in responsive to drought and salinity stress in Medicago truncatula.

Flavonoid compounds are widely present in various organs and tissues of different plants, playing important roles when plants are exposed to abiotic stresses. Different types of flavonoids are biosynthesized by a series of enzymes that are encoded by a range of gene families. In this study, a total of 63 flavonoid pathway genes were identified from the genome of Medicago truncatula. Gene structure analysis revealed that they all have different gene structure, with most CHS genes containing only one intron. Additionally, analysis of promoter sequences revealed that many cis-acting elements responsive to abiotic stress are located in the promoter region of flavonoid pathway genes. Furthermore, analysis on M. truncatula gene chip data revealed significant changes in expression level of most flavonoid pathway genes under the induction of salt or drought treatment. qRT-PCR further confirmed significant increase in expression level of several flavonoid pathway genes under NaCl and mannitol treatments, with CHS1, CHS9, CHS10, F3'H4 and F3'H5 genes showing significant up-regulation, indicating they are key genes in response to abiotic stress in M. truncatula. In summary, our study identified key flavonoid pathway genes that were involved in salt and drought response, which provides important insights into possible modification of flavonoid pathway genes for molecular breeding of forage grass with improved abiotic resistance.

The HD-ZIP IV transcription factor GLABRA2 acts as an activator for proanthocyanidin biosynthesis in Medicago truncatula seed coat.

Proanthocyanidins (PAs), a group of flavonoids, are found in leaves, flowers, fruits, and seed coats of many plant species. PAs are primarily composed of epicatechin units in the seed coats of the model legume species, Medicago truncatula. It can be synthesized from two separate pathways, the leucoanthocyanidin reductase (MtLAR) pathway and the anthocyanidin synthase (MtANS) pathway, which produce epicatechin through anthocyanidin reductase (MtANR). These pathways are mainly controlled by the MYB-bHLH-WD40 (MBW) ternary complex. Here, we characterize a class IV homeodomain-leucine zipper (HD-ZIP IV) transcription factor, GLABRA2 (MtGL2), which contributes to PA biosynthesis in the seed coat of M. truncatula. Null mutation of MtGL2 results in dark brown seed coat, which is accompanied by reduced PAs accumulation and increased anthocyanins content. The MtGL2 gene is predominantly expressed in the seed coat during the early stages of seed development. Genetic and molecular analyses indicate that MtGL2 positively regulates PA biosynthesis by directly activating the expression of MtANR. Additionally, our results show that MtGL2 is strongly induced by the MBW activator complexes that are involved in PA biosynthesis. Taken together, our results suggest that MtGL2 acts as a novel positive regulator in PA biosynthesis, expanding the regulatory network and providing insights for genetic engineering of PA production.

Over-expression of Medicago Acyl-CoA-binding 2 genes enhance salt and drought tolerance in Arabidopsis.

Acyl-CoA-binding proteins (ACBPs) are mainly involved in acyl-CoA ester binding and trafficking in eukaryotic cells, and they function in lipid metabolism, membrane biosynthesis, cellular signaling, stress response, disease resistance, and other biological activities in plants. However, the roles of ACBP family members in Medicago remain unclear. In this study, a total of eight ACBP genes were identified in the genome of Medicago truncatula and Medicago sativa, and they were clustered into four sub-families (Class I-IV). Many cis-acting elements related to abiotic response were identified in the promoter region of these ACBP genes, in particular light-responsive elements. These ACBP genes exhibited distinct expression pattern in various tissues, and the expression level of MtACBP1/MsACBP1 and MtACBP2/MsACBP2 gene pairs were significantly increased under NaCl treatment. Subcellular localization analysis showed that MtACBP1/MsACBP1 and MtACBP2/MsACBP2 were localized in the endoplasmic reticulum of tobacco epidermal cells. Arabidopsis seedlings over-expressing MtACBP2/MsACBP2 displayed increased root length than the wild type under short light, Cu2+, ABA, PEG, and NaCl treatments. Over-expression of MtACBP2/MsACBP2 also significantly enhanced Arabidopsis tolerance under NaCl and PEG treatments in mature plants. Collectively, our study identified salt and drought responsive ACBP genes in Medicago and verified their functions in increasing resistance against salt and drought stresses.

Medicago truncatula ß-glucosidase 17 contributes to drought and salt tolerance through antioxidant flavonoid accumulation.

Flavonoids are usually present in forms of glucosides in plants, which could be catabolized by ß-glucosidase (BGLU) to form their corresponding flavonoid aglycones. In this study, we isolated three abiotic-responsive BGLU genes (MtBGLU17, MtBGLU21 and MtBGLU22) from Medicago truncatula, and found only the recombinant MtBGLU17 protein could catalyse the hydrolysis of flavonoid glycosides. The recombinant MtBGLU17 protein is active towards a variety of flavonoid glucosides, including glucosides of flavones (apigenin and luteolin), flavonols (kaempferol and quercetin), isoflavones (genistein and daidzein) and flavanone (naringenin). In particular, the recombinant MtBGLU17 protein preferentially hydrolyses flavonoid-7-O-glucosides over their corresponding 3-O-glucosides. The content of luteoin-7-O-glucoside was reduced in the MtBGLU17 overexpression plants but increased in the Tnt-1 insertional mutant lines, whereas luteoin content was increased in the MtBGLU17 overexpression plants but reduced in the Tnt-1 insertional mutant lines. Under drought and salt (NaCl) treatment, the MtBGLU17 overexpression lines showed relatively higher DPPH content, and higher CAT and SOD activity than the wild type control. These results indicated that overexpression lines of MtBGLU17 possess higher antioxidant activity and thus confer drought and salt tolerance, implying MtBGLU17 could be potentially used as a candidate gene to improve plant abiotic stress tolerance.

Overexpression of MsNIP2 improves salinity tolerance in Medicago sativa.

Alfalfa (Medicago sativa) is one of the most widely cultivated forage crops in the world. However, alfalfa yield and quality are adversely affected by salinity stress. Nodulin 26-like intrinsic proteins (NIPs) play essential roles in water and small molecules transport and response to salt stress. Here, we isolated a salt stress responsive MsNIP2 gene and demonstrated its functions by overexpression in alfalfa. The open reading frame of MsNIP2 is 816 bp in length, and it encodes 272 amino acids. It has six transmembrane domains and two NPA motifs. MsNIP2 showed high identity to other known NIP proteins, and its tertiary model was similar to the crystal structure of OsNIP2-1 (7cjs) tetramer. Subcellular localization analysis showed that MsNIP2 protein fused with green fluorescent protein (GFP) was localized to the plasma membrane. Transgenic alfalfa lines overexpressing MsNIP2 showed significantly higher height and branch number compared with the non-transgenic control. The POD and CAT activity of the transgenic alfalfa lines was significantly increased and their MDA content was notably reduced compared with the control group under the treatment of NaCl. The transgenic lines showed higher capability in scavenging oxygen radicals with lighter NBT staining than the control under salt stress. The transgenic lines showed relative lower water loss rate and electrolyte leakage, but relatively higher Na+ content than the control line under salt stress. The relative expression levels of abiotic-stress-related genes (MsHSP23, MsCOR47, MsATPase, and MsRD2) in three transgenic lines were compared with the control, among them, only the expression of MsCOR47 was up-regulated. Consequently, this study offers a novel perspective for exploring the function of MsNIP2 in improving salt tolerance of alfalfa.

An inducible CRISPR activation tool for accelerating plant regeneration.

The inducible CRISPR activation (CRISPR-a) system offers unparalleled precision and versatility for regulating endogenous genes, making it highly sought after in plant research. In this study, we developed a chemically inducible CRISPR-a tool for plants called ER-Tag by combining the LexA-VP16-ER inducible system with the SunTag CRISPR-a system. We systematically compared different induction strategies and achieved high efficiency in target gene activation. We demonstrated that guide RNAs can be multiplexed and pooled for large-scale screening of effective morphogenic genes and gene pairs involved in plant regeneration. Further experiments showed that induced activation of these morphogenic genes can accelerate regeneration and improve regeneration efficiency in both eudicot and monocot plants, including alfalfa, woodland strawberry, and sheepgrass. Our study expands the CRISPR toolset in plants and provides a powerful new strategy for studying gene function when constitutive expression is not feasible or ideal.

The APETALA2-MYBL2 module represses proanthocyanidin biosynthesis by affecting formation of the MBW complex in seeds of Arabidopsis thaliana.

Proanthocyanidins (PAs) are the second most abundant plant phenolic natural products. PA biosynthesis is regulated by the well-documented MYB/bHLH/WD40 (MBW) complex, but how this complex itself is regulated remains ill defined. Here, in situ hybridization and ß-glucuronidase staining show that APETALA2 (AP2), a well-defined regulator of flower and seed development, is strongly expressed in the seed coat endothelium, where PAs accumulate. AP2 negatively regulates PA content and expression levels of key PA pathway genes. AP2 activates MYBL2 transcription and interacts with MYBL2, a key suppressor of the PA pathway. AP2 exerts its function by directly binding to the AT-rich motifs near the promoter region of MYBL2. Molecular and biochemical analyses revealed that AP2 forms AP2-MYBL2-TT8/EGL3 complexes, disrupting the MBW complex and thereby repressing expression of ANR, TT12, TT19, and AHA10. Genetic analyses revealed that AP2 functions upstream of MYBL2, TT2, and TT8 in PA regulation. Our work reveals a new role of AP2 as a key regulator of PA biosynthesis in Arabidopsis. Overall, this study sheds new light on the comprehensive regulation network of PA biosynthesis as well as the dual regulatory roles of AP2 in seed development and accumulation of major secondary metabolites in Arabidopsis.

Two CYP93A enzymes play a dual role in isoflavonoid biosynthesis in Glycine max L.

Glycine max L. is rich in isoflavonoids with diverse biological activities. However, isoflavonoid biosynthetic pathway is not fully elucidated in soybean. In the present study, we investigated characteristics of all the thirteen CYP93 subfamily members, and found GmCYP93A1, GmCYP93A2, and GmCYP93A3 are closely clustered, preferentially expressed in roots, and highly inducible by elicitor. When expressed in yeast, GmCYP93A1 was active towards liquiritigenin, naringenin, and 3,9-dihydroxyptercarpan, GmCYP93A2 towards 3,9-dihydroxyptercarpan with strict substrate specificity, whereas GmCYP93A3 did not show any activity towards all the tested substrates. Both GmCYP93A1 and GmCYP93A2 could catalyze 3,9-dihydroxyptercarpan into daidzein and glycinol, with both hydroxylation and aryl migration activity. Site-directed mutagenesis assays revealed that mutation in Thr446 to Ser446 in heme-binding domain increased the enzyme activity of GmCYP93A1 towards 3,9-dihydroxyptercarpan, which highlights its key amino acid residues as shown with its molecular docking with 3,9-dihydroxyptercarpan and HEM. Overexpression of GmCYP93A1 and GmCYP93A2 in the soybean hairy roots reduced the content of daidzein, whereas knockdown of these two genes increased genistein content, indicating changes in expression level of GmCYP93A1 and GmCYP93A2 altered isoflavonoid flux in soybean. Our studies on the activity of GmCYP93A1 and GmCYP93A2 enriched diverse functions of CYP93 subfamily in soybean isoflavonoid pathway, which is valuable for further understanding and bioengineering of isoflavonoid pathway in soybean.

New dual functional CYP450 gene involves in isoflavone biosynthesis in Glycine max L.

Glycine max L. accumulates a large amount of isoflavonoid compounds, which is beneficial for plant defense, plant-microbe symbiotic interactions, and human health. Several CYP450 subfamily genes are involved in the flavonoid biosynthetic pathway in plants. In the present study, we found 24 CYP82 subfamily genes were differentially expressed in various tissues of soybean, in Phytophthora sojae-infected soybean varieties and in soybean hairy roots treated with cell wall glucan elicitor. Six of them (GmCYP82A2, GmCYP82A3, GmCYP82A4, GmCYP82A23, GmCYP82C20 and GmCYP82D26) were co-expressed with other known isoflavonoid pathway genes in soybean. Their enzymatic activity in yeast feeding assays showed that only GmCYP82D26 was able to convert naringenin to daidzein with both aryl migration and dehydration function. When GmCYP82D26 was over-expressed in soybean hairy roots, the contents of the two major isoflavonoid aglycones in soybean (daidzein and genistein) were reduced, but total flavonoids were not affected. When GmCYP82D26 was suppressed by RNAi in the hairy roots, daidzein content was decreased but genistein content was increased, with unchanged total flavonoid content. GmCYP82D26 was found to be localized in the endoplasmic reticulum at subcellular level when transiently expressed in tobacco leaf epidermis. GmCYP82D26 gene was preferentially expressed in roots, with low expression level in other tissues in soybean. Homology modeling and molecular docking showed that GmCYP82D26 could form hydrogen bond with both HEM and naringenin at C5-OH and C4 carbonyl. All these results indicated that GmCYP82D26 possesses new and dual enzymatic activity, which bridges the two branches (daidzein and genistein branch) of isoflavonoid pathway in soybean.

The discovery of a key prenyltransferase gene assisted by a chromosome-level Epimedium pubescens genome.

Epimedium pubescens is a species of the family Berberidaceae in the basal eudicot lineage, and a main plant source for the traditional Chinese medicine "Herba Epimedii". The current study achieved a chromosome-level genome assembly of E. pubescens with the genome size of 3.34 Gb, and the genome guided discovery of a key prenyltransferase (PT) in E. pubescens. Our comparative genomic analyses confirmed the absence of Whole Genome Triplication (WGT-γ) event shared in core eudicots and further revealed the occurrence of an ancient Whole Genome Duplication (WGD) event approximately between 66 and 81 Million Years Ago (MYA). In addition, whole genome search approach was successfully applied to identify 19 potential flavonoid PT genes and an important flavonoid PT (EpPT8) was proven to be an enzyme for the biosynthesis of medicinal compounds, icaritin and its derivatives in E. pubescens. Therefore, our results not only provide a good reference genome to conduct further molecular biological studies in Epimedium genus, but also give important clues for synthetic biology and industrial production of related prenylated flavonoids in future.

Genome-wide characterization and expression analysis of the HAK gene family in response to abiotic stresses in Medicago.

The high-affinity K+ transporter (HAK) family plays a vital role in K+ uptake and transport as well as in salt and drought stress responses. In the present study, we identified 22 HAK genes in each Medicago truncatula and Medicago sativa genome. Phylogenetic analysis suggested that these HAK proteins could be divided into four clades, and the members of the same subgroup share similar gene structure and conserved motifs. Many cis-acting elements related with defense and stress were found in their promoter region. In addition, gene expression profiles analyzed with genechip and transcriptome data showed that these HAK genes exhibited distinct expression pattern in different tissues, and in response to salt and drought treatments. Furthermore, co-expression analysis showed that 6 homologous HAK hub gene pairs involved in direct network interactions. RT-qPCR verified that the expression level of six HAK gene pairs was induced by NaCl and mannitol treatment to different extents. In particular, MtHK2/7/12 from M. truncatula and MsHAK2/6/7 from M. sativa were highly induced. The expression level of MsHAK1/2/11 determined by RT-qPCR showed significantly positive correlation with transcriptome data. In conclusion, our study shows that HAK genes play a key role in response to various abiotic stresses in Medicago, and the highly inducible candidate HAK genes could be used for further functional studies and molecular breeding in Medicago.

A Novel 3-O-rhamnoside: 2″-O-xylosyltransferase Responsible for Terminal Modification of Prenylflavonol Glycosides in Epimedium pubescens Maxim.

Prenylated flavonol glycosides in Epimedium plants, as key medicinal components, are known to have great pharmaceutical activities for human health. Among the main prenylated flavonol glycosides, the modification mechanism of different sugar moieties is still not well understood. In the current study, a novel prenylated flavonol rhamnoside xylosyltransferase gene (EpF3R2″XylT) was cloned from E. pubescens, and the enzymatic activity of its decoding proteins was examined in vitro with different prenylated flavonol rhamnoside substrates and different 3-O-monosaccharide moieties. Furthermore, the functional and structural domains of EpF3R2″XylT were analyzed by bioinformatic approaches and 3-D protein structure remodeling. In summary, EpF3R2″XylT was shown to cluster with GGT (glycosyltransferase that glycosylates sugar moieties of glycosides) through phylogenetic analysis. In enzymatic analysis, EpF3R2″XylT was proven to transfer xylose moiety from UDP-xylose to prenylated flavonol rhamnoside at the 2″-OH position of rhamnose. The analysis of enzymatic kinetics showed that EpF3R2″XylT had the highest substrate affinity toward icariin with the lowest Km value of 75.96 ± 11.91 mM. Transient expression of EpF3R2″XylT in tobacco leaf showed functional production of EpF3R2″XylT proteins in planta. EpF3R2″XylT was preferably expressed in the leaves of E. pubescens, which is consistent with the accumulation levels of major prenylflavonol 3-O-triglycoside. The discovery of EpF3R2″XylT will provide an economical and efficient alternative way to produce prenylated flavonol trisaccharides through the biosynthetic approach.

Transcriptome and functional analyses reveal ERF053 from Medicago falcata as key regulator in drought resistances.

Medicago falcata L. is an important legume forage grass with strong drought resistant, which could be utilized as an important gene pool in molecular breed of forage grass. In this study, M. falcata seedlings were treated with 400 mM mannitol to simulate drought stress, and the morphological and physiological changes were investigated, as well as the transcriptome changes of M. falcata seedlings at different treatment time points (0 h, 2 h, 6 h, 12 h, 24 h, 36 h and 48 h). Transcriptome analyses revealed four modules were closely related with drought response in M. falcata by WGCNA analysis, and four ERF transcription factor genes related with drought stress were identified (MfERF053, MfERF9, MfERF034 and MfRAP2.1). Among them, MfERF053 was highly expressed in roots, and MfERF053 protein showed transcriptional activation activity by transient expression in tobacco leaves. Overexpression of MfERF053 in Arabidopsis improved root growth, number of lateral roots and fresh weight under drought, salt stress and exogenous ABA treatments. Transgenic Arabidopsis over-expressing MfERF053 gene grew significantly better than the wild type under both drought stress and salt stress when grown in soil. Taken together, our strategy with transcriptome combined WGCNA analyses identified key transcription factor genes from M. falcata, and the selected MfERF053 gene was verified to be able to enhance drought and salt resistance when over-expressed in Arabidopsis.

Genome-wide analysis of UGT gene family identified key gene for the biosynthesis of bioactive flavonol glycosides in Epimedium pubescens Maxim.

Epimedium pubescens Maxim. is a well-known traditional Chinese medicinal herb with flavonol glycosides as the major pharmaceutically active compounds. UDP-glycosyltransferases (UGTs) are a group of enzymes responsible for the glycosylation of flavonoid glycosides. In this study, a genome-wide analysis was performed to identify UGT family genes in E. pubescens. As a result, a total of 339 putative UGT genes were identified, which represents the largest UGT gene family known thus far, implying a significant expansion of the UGT gene family in E. pubescens. All EpUGTs were unevenly distributed across six chromosomes, and they were classified into 17 major groups. The expression profiles showed that UGT genes were differentially expressed in roots, leaves, flowers, shoots and fruits. In particular, several EpUGTs were highly induced by high light intensity, which was consistent with the accumulation level of bioactive flavonoids in E. pubescens. Six UGT79 genes that were preferentially expressed in roots or leaves were successfully expressed in E. coli, and only the recombinant EpGT60 protein was found to be active toward 8-prenylkaempferol and icaritin to produce the key bioactive compounds baohuoside II and baohuoside I. The optimal temperature, pH, k m and V max were determined for the recombinant EpGT60 protein. In addition, expression of recombinant EpGT60 in E. coli cell culture led to successful production of baohuoside II when fed 8-prenylkaempferol. Our study provides a foundation for further functional characterization of UGT genes in E. pubescens and provides key candidate genes for bioengineering bioactive flavonoids in E. pubescens.

ARF2 positively regulates flavonols and proanthocyanidins biosynthesis in Arabidopsis thaliana.

MAIN CONCLUSION: Auxin response factor 2 acts as a positive regulator to fine-tune the spatial and temporal accumulation of flavonoid compounds, mainly flavonols and proanthocyanidins in Arabidopsis. Auxin response factor (ARF) proteins are reported to involve in auxin-mediated regulation of flavonoid biosynthesis. However, the detailed regulation mechanism of ARF remains still unknown. Here, we provide genetic and molecular evidence that one of the twenty-three ARF members-ARF2-positively regulates flavonoid biosynthesis at multi-level in tissue-specific manner in Arabidopsis thaliana. Loss-of-function mutation of ARF2 led to significant reduction in flavonoid content (e.g., flavonols and proanthocyanidins) in the seedlings and seeds of the Arabidopsis arf2 mutants. Over-expression of ARF2 increased flavonols and proanthocyanidins content in Arabidopsis. Additionally, the changes of flavonoid content correlate well with the transcript abundance of several regulatory genes (e.g., MYB11, MYB12, MYB111, TT2, and GL3), and key biosynthetic genes (e.g., CHS, F3'H, FLS, ANS, ANR, TT12, TT19, and TT15), in the arf2 mutant and ARF2 over-expression lines. Transient transactivation assays with site-directed mutagenesis confirmed that ARF2 directly regulates the expression of MYB12 and FLS genes in the flavonol pathway and ANR in the proanthocyanidin pathway, and indirectly regulates MYB11 and MYB111 genes in the flavonol pathway, and ANS, TT12, TT19 and TT15 genes in the proanthocyanidin pathway. Further genetic results indicated that ARF2 acts upstream of MYB12 to regulate flavonol accumulation, and of TT2 to regulate proanthocyanidins accumulation. In particular, yeast two-hybrid assays revealed that ARF2 physically interacts with TT2, a master regulator of proanthocyanidins biosynthesis. Combined together, these results indicated that ARF2 functions as a positive regulator for the fine-tuned spatial and temporal regulation of flavonoids (mainly flavonols and proanthocyanidins) accumulation in Arabidopsis.

Analyses on Flavonoids and Transcriptome Reveals Key MYB Gene for Proanthocyanidins Regulation in Onobrychis Viciifolia.

Onobrychis viciifolia (sainfoin) is one of the most high-quality legume forages, which is rich in proanthocyanidins that is beneficial for the health and production of animals. In this study, proanthocyanidins and total flavonoids in leaves of 46 different sainfoin germplasm resources were evaluated, and it showed that soluble proanthocyanidin contents varied greatly in these sainfoin germplasm resources, but total flavonoids did not show significant difference. Transcriptome sequencing with high and low proanthocyanidins sainfoin resulted in the identification of totally 52,926 unigenes in sainfoin, and they were classed into different GOC categories. Among them, 1,608 unigenes were differentially expressed in high and low proanthocyanidins sainfoin samples, including 1,160 genes that were upregulated and 448 genes that were downregulated. Analysis on gene enrichment via KEGG annotation revealed that the differentially expressed genes were mainly enriched in the phenylpropanoid biosynthetic pathway and the secondary metabolism pathway. We also analyzed the expression levels of structural genes of the proanthocyanidin/flavonoid pathway in roots, stems, and leaves in the high proanthocyanidin sainfoin via RT-qPCR and found that these genes were differentially expressed in these tissues. Among them, the expression levels of F3'5'H and ANR were higher in leaves than in roots or stems, which is consistent with proanthocyanidins content in these tissues. Among MYB genes that were differentially expressed, the expression of OvMYBPA2 was relatively high in high proanthocyanidin sainfoin. Over-expression level of OvMYBPA2 in alfalfa hairy roots resulted in decreased anthocyanin content but increased proanthocyanidin content. Our study provided transcriptome information for further functional characterization of proanthocyanidin biosynthesis-related genes in sainfoin and candidate key MYB genes for bioengineering of proanthocyanidins in plants.

Genome-Wide Identification and Characterization of Growth Regulatory Factor Family Genes in Medicago.

Growth Regulatory Factors (GRF) are plant-specific transcription factors that play critical roles in plant growth and development as well as plant tolerance against stress. In this study, a total of 16 GRF genes were identified from the genomes of Medicago truncatula and Medicago sativa. Multiple sequence alignment analysis showed that all these members contain conserved QLQ and WRC domains. Phylogenetic analysis suggested that these GRF proteins could be classified into five clusters. The GRF genes showed similar exon-intron organizations and similar architectures in their conserved motifs. Many stress-related cis-acting elements were found in their promoter region, and most of them were related to drought and defense response. In addition, analyses on microarray and transcriptome data indicated that these GRF genes exhibited distinct expression patterns in various tissues or in response to drought and salt treatments. In particular, qPCR results showed that the expression levels of gene pairs MtGRF2-MsGRF2 and MtGRF6-MsGRF6 were significantly increased under NaCl and mannitol treatments, indicating that they are most likely involved in salt and drought stress tolerance. Collectively, our study is valuable for further investigation on the function of GRF genes in Medicago and for the exploration of GRF genes in the molecular breeding of highly resistant M. sativa.

Identification and Characterization of Abiotic Stress-Responsive NF-YB Family Genes in Medicago.

Nuclear factor YB (NF-YB) are plant-specific transcription factors that play a critical regulatory role in plant growth and development as well as in plant resistance against various stresses. In this study, a total of 49 NF-YB genes were identified from the genomes of Medicago truncatula and Medicago sativa. Multiple sequence alignment analysis showed that all of these NF-YB members contain DNA binding domain, NF-YA interaction domain and NF-YC interaction domain. Phylogenetic analysis suggested that these NF-YB proteins could be classified into five distinct clusters. We also analyzed the exon-intron organizations and conserved motifs of these NF-YB genes and their deduced proteins. We also found many stress-related cis-acting elements in their promoter region. In addition, analyses on genechip for M. truncatula and transcriptome data for M. sativa indicated that these NF-YB genes exhibited a distinct expression pattern in various tissues; many of these could be induced by drought and/or salt treatments. In particular, RT-qPCR analysis revealed that the expression levels of gene pairs MsNF-YB27/MtNF-YB15 and MsNF-YB28/MtNF-YB16 were significantly up-regulated under NaCl and mannitol treatments, indicating that they are most likely involved in salt and drought stress response. Taken together, our study on NF-YB family genes in Medicago is valuable for their functional characterization, as well as for the application of NF-YB genes in genetic breeding for high-yield and high-resistance alfalfa.

Full-Length Transcriptomics Reveals Complex Molecular Mechanism of Salt Tolerance in Bromus inermis L.

Bromus inermis L. (commonly known as smooth bromegrass) is a grass species with high nutritional value, great palatability, cold tolerance, and grazing resistance, which has been widely cultivated for pasture and sand fixation in northern and northwestern China. Salt stress is a main environmental factor limiting growth and production of smooth bromegrass. In this study, we performed PacBio Iso-Seq to construct the first full-length transcriptome database for smooth bromegrass under 300 mM NaCl treatment at different time points. Third-generation full-length transcriptome sequencing yielded 19.67 G polymerase read bases, which were assembled into 355,836 full-length transcripts with an average length of 2,542 bp. A total of 116,578 differentially expressed genes were obtained by comparing the results of third-generation sequencing and second-generation sequencing. GO and KEGG enrichment analyses revealed that multiple pathways were differently activated in leaves and roots. In particular, a number of genes participating in the molecular network of plant signal perception, signal transduction, transcription regulation, antioxidant defense, and ion regulation were affected by NaCl treatment. In particular, the CBL-CIPK, MAPK, ABA signaling network, and SOS core regulatory pathways of Ca2+ signal transduction were activated to regulate salt stress response. In addition, the expression patterns of 10 salt-responsive genes were validated by quantitative real-time PCR, which were consistent with those detected by RNA-Seq. Our results reveal the molecular regulation of smooth bromegrass in response to salt stress, which are important for further investigation of critical salt responsive genes and molecular breeding of salt-tolerant smooth bromegrass.