Comprehensive identification and characterization of abiotic stress and hormone responsive glycosyl hydrolase family 1 genes in Medicago truncatula

Integrated physiological, biochemical, and transcriptomic analyses of Bruguiera gymnorhiza leaves under long-term copper stress: Stomatal size, wax crystals and composition

Copper (Cu) is a necessary mineral nutrient for plant growth and development and is involved in several morphological, physiological, and biochemical processes; however, high concentrations of Cu can negatively impact these processes. The role of stomata in responding to various biotic and abiotic stimuli has not been studied in Bruguiera gymnorhiza, particularly in terms of their coordinated interactions at the molecular, physiological, and biochemical levels. Moreover, numerous plants employ strategies such as the presence of thick waxy cuticles on their leaf epidermis and the closing of stomata to reduce water loss. Thus, this study investigates the accumulation of Cu in B. gymnorhiza and its effect on leaf morphology and the molecular response under different Cu treatments (0, 200, and 400 mg L⁻¹, Cu0, Cu200, and Cu400, respectively) during a two years stress period. The results show that Cu stress affected accumulation and transport, increased the activities of peroxidase and ascorbate peroxidase, concentrations of soluble sugar, proline, and H2O2, and decreased the activity of catalase and content of malondialdehyde. Also, Cu-induced stress decreased the uptake of phosphorus and nitrogen and inhibited plant photosynthesis, which consequently led to reduced plant growth. Scanning electron microscopy combined with gas chromatography-mass spectrometry showed that B. gymnorhiza leaves had higher wax crystals and compositions under increased Cu stress, which forced the leaf's stomata to be closed. Also, the contents of alkanes, alcohols, primary alcohol levels (C26:0, C28:0, C30:0, and C32:0), n-Alkanes (C29 and C30), and other wax loads were significantly higher, while fatty acid (C12, C16, and C18) was lower in Cu200 and Cu400 compared to Cu0. Furthermore, the transcriptomic analyses revealed 1240 (771 up- and 469 downregulated), 1000 (723 up- and 277 down-regulated), and 1476 (808 up- and 668 downregulated) differentially expressed genes in Cu0 vs Cu200, Cu0 vs Cu400, and Cu200 vs Cu400, respectively. RNA-seq analyses showed that Cu mainly affected eight pathways, including photosynthesis, cutin, suberin, and wax biosynthesis. This study provides a reference for understanding mangrove response to heavy metal stress and developing novel management practices.

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.

Long noncoding RNA TRABA suppresses β-glucosidase-encoding BGLU24 to promote salt tolerance in cotton

Salt stress severely damages the growth and yield of crops. Recently, long noncoding RNAs (lncRNAs) were demonstrated to regulate various biological processes and responses to environmental stresses. However, the regulatory mechanisms of lncRNAs in cotton (Gossypium hirsutum) response to salt stress are still poorly understood. Here, we observed that a lncRNA, trans acting of BGLU24 by lncRNA (TRABA), was highly expressed while GhBGLU24-A was weakly expressed in a salt-tolerant cotton accession (DM37) compared to a salt-sensitive accession (TM-1). Using TRABA as an effector and proGhBGLU24-A-driven GUS as a reporter, we showed that TRABA suppressed GhBGLU24-A promoter activity in double transgenic Arabidopsis (Arabidopsis thaliana), which explained why GhBGLU24-A was weakly expressed in the salt-tolerant accession compared to the salt-sensitive accession. GhBGLU24-A encodes an endoplasmic reticulum (ER)-localized β-glucosidase that responds to salt stress. Further investigation revealed that GhBGLU24-A interacted with RING-type E3 ubiquitin ligase (GhRUBL). Virus-induced gene silencing (VIGS) and transgenic Arabidopsis studies revealed that both GhBGLU24-A and GhRUBL diminish plant tolerance to salt stress and ER stress. Based on its substantial effect on ER-related degradation (ERAD)-associated gene expression, GhBGLU24-A mediates ER stress likely through the ERAD pathway. These findings provide insights into the regulatory role of the lncRNA TRABA in modulating salt and ER stresses in cotton and have potential implications for developing more resilient crops.

Applications of Virus-Induced Gene Silencing in Cotton

Virus-induced gene silencing (VIGS) is an RNA-mediated reverse genetics technique that has become an effective tool to investigate gene function in plants. Cotton is one of the most important economic crops globally. In the past decade, VIGS has been successfully applied in cotton functional genomic studies, including those examining abiotic and biotic stress responses and vegetative and reproductive development. This article summarizes the traditional vectors used in the cotton VIGS system, the visible markers used for endogenous gene silencing, the applications of VIGS in cotton functional genomics, and the limitations of VIGS and how they can be addressed in cotton.

Genome-wide identification and expression analysis of the glycosyl hydrolase family 1 genes in Medicago sativa revealed their potential roles in response to multiple abiotic stresses

Glycoside hydrolase family 1 (GH1) β-glucosidases (BGLUs), are encoded by a large number of genes, which participate in the development and stress response of plants, particularly under biotic and abiotic stresses through the activation of phytohormones. However, there are few studies systematically analyzing stress or hormone-responsive BGLU genes in alfalfa. In this study, a total of 179 BGLU genes of the glycoside hydrolase family 1 were identified in the genome of alfalfa, and then were classified into five distinct clusters. Sequence alignments revealed several conserved and unique motifs among these MsBGLU proteins. Many cis-acting elements related to abiotic stresses and phytohormones were identified in the promoter of some MsBGLUs. Moreover, RNA-seq and RT-qPCR analyses showed that these MsBGLU genes exhibited distinct expression patterns in response to different abiotic stress and hormonal treatments. In summary, this study suggests that MsBGLU genes play crucial roles in response to various abiotic stresses and hormonal responses, and provides candidate genes for stress tolerance breeding in alfalfa.

QTL mapping reveals key factors related to the isoflavone contents and agronomic traits of soybean (Glycine max)

BACKGROUND

Soybean is a valuable source of edible protein and oil, as well as secondary metabolites that can be used in food products, cosmetics, and medicines. However, because soybean isoflavone content is a quantitative trait influenced by polygenes and environmental interactions, its genetic basis remains unclear.

RESULTS

This study was conducted to identify causal quantitative trait loci (QTLs) associated with soybean isoflavone contents. A mutant-based F2 population (190 individuals) was created by crossing the Korean cultivar Hwanggeum with low isoflavone contents (1,558 µg g-1) and the soybean mutant DB-088 with high isoflavone contents (6,393 µg g-1). A linkage map (3,049 cM) with an average chromosome length of 152 cM was constructed using the 180K AXIOM® SoyaSNP array. Thirteen QTLs related to agronomic traits were mapped to chromosomes 2, 3, 11, 13, 19, and 20, whereas 29 QTLs associated with isoflavone contents were mapped to chromosomes 1, 3, 8, 11, 14, 15, and 17. Notably, the qMGLI11, qMGNI11, qADZI11, and qTI11, which located Gm11_9877690 to Gm11_9955924 interval on chromosome 11, contributed to the high isoflavone contents and explained 11.9% to 20.1% of the phenotypic variation. This QTL region included four candidate genes, encoding β-glucosidases 13, 14, 17-1, and 17-2. We observed significant differences in the expression levels of these genes at various seed developmental stages. Candidate genes within the causal QTLs were functionally characterized based on enriched GO terms and KEGG pathways, as well as the results of a co-expression network analysis. A correlation analysis indicated that certain agronomic traits (e.g., days to flowering, days to maturity, and plant height) are positively correlated with isoflavone content.

CONCLUSIONS

Herein, we reported that the major QTL associated with isoflavone contents was located in the interval from Gm11_9877690 to Gm11_9955924 (78 kb) on chromosome 11. Four β-glucosidase genes were identified that may be involved in high isoflavone contents of soybean DB-088. Thus, the mutant alleles from soybean DB-088 may be useful for marker-assisted selection in developing soybean lines with high isoflavone contents and superior agronomic traits.

The saponin bomb: a nucleolar-localized β-glucosidase hydrolyzes triterpene saponins in Medicago truncatula

Plants often protect themselves from their own bioactive defense metabolites by storing them in less active forms. Consequently, plants also need systems allowing correct spatiotemporal reactivation of such metabolites, for instance under pathogen or herbivore attack. Via co-expression analysis with public transcriptomes, we determined that the model legume Medicago truncatula has evolved a two-component system composed of a β-glucosidase, denominated G1, and triterpene saponins, which are physically separated from each other in intact cells. G1 expression is root-specific, stress-inducible, and coregulated with that of the genes encoding the triterpene saponin biosynthetic enzymes. However, the G1 protein is stored in the nucleolus and is released and united with its typically vacuolar-stored substrates only upon tissue damage, partly mediated by the surfactant action of the saponins themselves. Subsequently, enzymatic removal of carbohydrate groups from the saponins creates a pool of metabolites with an increased broad-spectrum antimicrobial activity. The evolution of this defense system benefited from both the intrinsic condensation abilities of the enzyme and the bioactivity properties of its substrates. We dub this two-component system the saponin bomb, in analogy with the mustard oil and cyanide bombs, commonly used to describe the renowned β-glucosidase-dependent defense systems for glucosinolates and cyanogenic glucosides.

Transcriptome profile analysis of Indian mustard (Brassica juncea L.) during seed germination reveals the drought stress-induced genes associated with energy, hormone, and phenylpropanoid pathways

Indian mustard (Brassica juncea L. Czern and Coss) is an important oil and vegetable crop frequently affected by seasonal drought stress during seed germination, which retards plant growth and causes yield loss considerably. However, the gene networks regulating responses to drought stress in leafy Indian mustard remain elusive. Here, we elucidated the underlying gene networks and pathways of drought response in leafy Indian mustard using next-generation transcriptomic techniques. Phenotypic analysis showed that the drought-tolerant leafy Indian mustard cv. 'WeiLiang' (WL) had a higher germination rate, antioxidant capacity, and better growth performance than the drought-sensitive cv. 'ShuiDong' (SD). Transcriptome analysis identified differentially expressed genes (DEGs) in both cultivars under drought stress during four germination time points (i.e., 0, 12, 24, and 36 h); most of which were classified as drought-responsive, seed germination, and dormancy-related genes. In the Kyoto Encyclopedia of Genes and Genome (KEGG) analyses, three main pathways (i.e., starch and sucrose metabolism, phenylpropanoid biosynthesis, and plant hormone signal transduction) were unveiled involved in response to drought stress during seed germination. Furthermore, Weighted Gene Co-expression Network Analysis (WGCNA) identified several hub genes (novel.12726, novel.1856, BjuB027900, BjuA003402, BjuA021578, BjuA005565, BjuB006596, novel.12977, and BjuA033308) associated with seed germination and drought stress in leafy Indian mustard. Taken together, these findings deepen our understanding of the gene networks for drought responses during seed germination in leafy Indian mustard and provide potential target genes for the genetic improvement of drought tolerance in this crop.

Starch and Sucrose Metabolism and Plant Hormone Signaling Pathways Play Crucial Roles in Aquilegia Salt Stress Adaption

Salt stress is one of the main abiotic stresses that strongly affects plant growth. Clarifying the molecular regulatory mechanism in ornamental plants under salt stress is of great significance for the ecological development of saline soil areas. Aquilegia vulgaris is a perennial with a high ornamental and commercial value. To narrow down the key responsive pathways and regulatory genes, we analyzed the transcriptome of A. vulgaris under a 200 mM NaCl treatment. A total of 5600 differentially expressed genes were identified. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis pointed out that starch and sucrose metabolism and plant hormone signal transduction were significantly improved. The above pathways played crucial roles when A. vulgaris was coping with salt stress, and their protein-protein interactions (PPIs) were predicted. This research provides new insights into the molecular regulatory mechanism, which could be the theoretical basis for screening candidate genes in Aquilegia.

Genome-wide identification and characterization of active ingredients related β-Glucosidases in Dendrobium catenatum

BACKGROUND

Dendrobium catenatum/D. officinale (here after D. catenatum), a well-known economically important traditional medicinal herb, produces a variety of bioactive metabolites including polysaccharides, alkaloids, and flavonoids with excellent pharmacological and clinical values. Although many genes associated with the biosynthesis of medicinal components have been cloned and characterized, the biosynthetic pathway, especially the downstream and regulatory pathway of major medicinal components in the herb, is far from clear. β-glucosidases (BGLUs) comprise a diverse group of enzymes that widely exist in plants and play essential functions in cell wall modification, defense response, phytohormone signaling, secondary metabolism, herbivore resistance, and scent release by hydrolyzing β-D-glycosidic bond from a carbohydrate moiety. The recent release of the chromosome-level reference genome of D. catenatum enables the characterization of gene families. Although the genome-wide analysis of the BGLU gene family has been successfully conducted in various plants, no systematic analysis is available for the D. catenatum. We previously isolated DcBGLU2 in the BGLU family as a key regulator for polysaccharide biosynthesis in D. catenatum. Yet, the exact number of DcBGLUs in the D. catenatum genome and their possible roles in bioactive compound production deserve more attention.

RESULTS

To investigate the role of BGLUs in active metabolites production, 22 BGLUs (DcBGLU1-22) of the glycoside hydrolase family 1 (GH1) were identified from D. catenatum genome. Protein prediction showed that most of the DcBGLUs were acidic and phylogenetic analysis classified the family into four distinct clusters. The sequence alignments revealed several conserved motifs among the DcBGLU proteins and analyses of the putative signal peptides and N-glycosylation site revealed that the majority of DcBGLU members dually targeted to the vacuole and/or chloroplast. Organ-specific expression profiles and specific responses to MeJA and MF23 were also determined. Furthermore, four DcBGLUs were selected to test their involvement in metabolism regulation. Overexpression of DcBGLU2, 6, 8, and 13 significantly increased contents of flavonoid, reducing-polysaccharide, alkaloid and soluble-polysaccharide, respectively.

CONCLUSION

The genome-wide systematic analysis identified candidate DcBGLU genes with possible roles in medicinal metabolites production and laid a theoretical foundation for further functional characterization and molecular breeding of D. catenatum.

Integrated physiological, transcriptomic and metabolomic analysis of the response of Trifolium pratense L. to Pb toxicity

Lead (Pb) interferes with plant gene expression, alters metabolite contents and affects plant growth. However, the molecular mechanism underlying the plant response to Pb is not completely understood. In the present study, Trifolium pratense L. was exposed to Pb concentrations of 0 (Pb0), 500 (Pb500), 1000 (Pb1000), 2000 (Pb2000) and 3000 (Pb3000) mg/kg in soils. Pb stress affected the ability of T. pratense to accumulate and transport Pb, increased the activity of peroxidase (POD) and the contents of malondialdehyde (MDA) and proline, decreased the amount of photosynthetic pigments and soluble proteins, and led to changes in growth and biomass. Transcriptomic and metabolomic analyses showed that Pb mainly affected eight pathways, and LHC, flavonoids, organic acids, amino acids and carbohydrates were upregulated or downregulated. Moreover, Pb500 induced the upregulation of serA, promoted the synthesis of citric acid, maintained photosynthetic pigment levels, and ultimately promoted an increase in stem length. Pb3000 induced the upregulation of ARF, GH3 and SAUR genes, but the saccharide contents and stem length decreased in response to Pb stress. We used a variety of methods to provide a molecular perspective on the mechanism underlying the response of T. pratense to Pb stress.

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.

Identification of Candidate Genes Associated With Tolerance to Apple Replant Disease by Genome-Wide Transcriptome Analysis

Apple replant disease (ARD) is a worldwide economic risk in apple cultivation for fruit tree nurseries and fruit growers. Several studies on the reaction of apple plants to ARD are documented but less is known about the genetic mechanisms behind this symptomatology. RNA-seq analysis is a powerful tool for revealing candidate genes that are involved in the molecular responses to biotic stresses in plants. The aim of our work was to find differentially expressed genes in response to ARD in Malus. For this, we compared transcriptome data of the rootstock ‘M9’ (susceptible) and the wild apple genotype M. ×robusta 5 (Mr5, tolerant) after cultivation in ARD soil and disinfected ARD soil, respectively. When comparing apple plantlets grown in ARD soil to those grown in disinfected ARD soil, 1,206 differentially expressed genes (DEGs) were identified based on a log2 fold change, (LFC) ≥ 1 for up– and ≤ −1 for downregulation (p < 0.05). Subsequent validation revealed a highly significant positive correlation (r = 0.91; p < 0.0001) between RNA-seq and RT-qPCR results indicating a high reliability of the RNA-seq data. PageMan analysis showed that transcripts of genes involved in gibberellic acid (GA) biosynthesis were significantly enriched in the DEG dataset. Most of these GA biosynthesis genes were associated with functions in cell wall stabilization. Further genes were related to detoxification processes. Genes of both groups were expressed significantly higher in Mr5, suggesting that the lower susceptibility to ARD in Mr5 is not due to a single mechanism. These findings contribute to a better insight into ARD response in susceptible and tolerant apple genotypes. However, future research is needed to identify the defense mechanisms, which are most effective for the plant to overcome ARD.

Multi-omics analyses revealed key factors involved in fluorescent carbon-dots-regulated secondary metabolism in Tetrastigma hemsleyanum

BACKGROUND

Luminescent nanomaterials (LNMs), especially newly-exploited fluorescent carbon-dots (CDs), have demonstrated promising candidates for sunlight harvesting and enhanced photosynthesis efficiency of crops. However, most of the studies focus on the design and synthesis of LNMs and primary metabolism in biomass acceleration, secondary metabolism that closely associated with the quality ingredients of plants is rarely mentioned.

RESULTS

UV-absorptive and water-soluble NIR-CDs were harvested via a facile microwave-assisted carbonization method. The effect and regulatory mechanism of NIR-CDs on the secondary metabolism and bioactive ingredients accumulation in Tetrastigma hemsleyanum were explored. A total of 191 differential secondary metabolites and 6874 differentially expressed genes were identified when the NIR-CDs were adopted for enhancing growth of T. hemsleyanum. The phenolic acids were generally improved, but the flavonoids were more likely to decrease. The pivotal differentially expressed genes were involved in biosynthesis of secondary metabolites, flavonoid biosynthesis, porphyrin and chlorophyll metabolism, etc. The gene-metabolite association was constructed and 44 hub genes highly related to quality ingredients accumulation and growth were identified, among which and the top 5 genes of the PPI network might be the key regulators.

CONCLUSION

This research provided key regulatory genes and differentially accumulating quality ingredients under NIR-CDs-treatment, which could provide a theoretical basis for expanding the applications of nanomaterial in industrial crop agriculture.

VqBGH40a isolated from Chinese wild Vitis quinquangularis degrades trans-piceid and enhances trans-resveratrol

Resveratrol (3,5,4'-trihydroxy-stilbene) is a phytoalexin that can prevent plants from pathogen attacks. Piceid is the glycosylation product of resveratrol and the main storage form of stilbenes in grapevines. Here, we reported the function of a β-glycoside hydrolase gene, VqBGH40a, from the Chinese wild grapevine Vitis quinquangularis accession Danfeng-2 in the regulation of plant resistance to powdery mildew (Uncinula necator). VqBGH40a belonging to β-glycoside hydrolase family 1 encoded 506 amino acids and was located on the cytomembrane. Its optimal induction condition was 28 or 30℃, for 4 h, with 0.1 mM IPTG in a prokaryotic expression system. Enzyme activity detection showed that purified VqBGH40a could hydrolyze trans-piceid to form trans-resveratrol in vitro. VqBGH40a was transiently overexpressed in Danfeng-2 leaves and then artificially inoculated with powdery mildew showed that VqBGH40a protein could hydrolyze trans-piceid in vivo. Additionally, a comparative family analysis between VqBGH40a and 38 VviBGHs was performed. Overall, these results demonstrate that VqBGH40a can hydrolyze trans-piceid, enhance trans-resveratrol content, and participate in the defense mechanism of grapevine against powdery mildew.

Genome-Wide Analysis of Glycoside Hydrolase Family 35 Genes and Their Potential Roles in Cell Wall Development in Medicago truncatula

Plant β-galactosidases (BGAL) function in various cell wall biogeneses and modifications, and they belong to the glycoside hydrolase family. However, the roles of BGAL family members in Medicago truncatula cell wall remodeling remain unclear. In this study, a total of 25 MtBGAL members of the glycoside hydrolase gene family 35 were identified, and they were clustered into nine sub-families. Many cis-acting elements possibly related to MeJA and abscisic acid responses were identified in the promoter region of the MtBGAL genes. Transcript analyses showed that these MtBGAL genes exhibited distinct expression patterns in various tissues and developing stem internodes. Furthermore, a stem-specific expression module associated with cell wall metabolic pathways was identified by weighted correlation network analysis (WGCNA). In particular, MtBGAL1 and MtBGAL23 within the stem-specific expression module were highly expressed in mature stems. In addition, several genes involved in lignin, cellulose, hemicellulose and pectin pathways were co-expressed with MtBGAL1 and MtBGAL23. It was also found that MtBGAL1 and MtBGAL23 were localized to the cell wall at the subcellular level, indicating their roles in the modification of cell wall metabolites in Medicago. As a whole, these results will be useful for further functional characterization and utilization of BGAL genes in cell wall modifications aiming to improve the quality of legume forage crops.

Identification and Characterization of Abiotic Stress Responsive CBL-CIPK Family Genes in Medicago

The calcineurin B-like protein (CBL) and CBL-interacting protein kinase (CIPK) play important roles in plant signal transduction and response to abiotic stress. Plants of Medicago genus contain many important forages, and their growth is often affected by a variety of abiotic stresses. However, studies on the CBL and CIPK family member and their function are rare in Medicago. In this study, a total of 23 CBL and 58 CIPK genes were identified from the genome of Medicago sativa as an important forage crop, and Medicaog truncatula as the model plant. Phylogenetic analysis suggested that these CBL and CIPK genes could be classified into five and seven groups, respectively. Moreover, these genes/proteins showed diverse exon-intron organizations, architectures of conserved protein motifs. Many stress-related cis-acting elements were found in their promoter region. In addition, transcriptional analyses showed that these CBL and CIPK genes exhibited distinct expression patterns in various tissues, and in response to drought, salt, and abscisic acid treatments. In particular, the expression levels of MtCIPK2 (MsCIPK3), MtCIPK17 (MsCIPK11), and MtCIPK18 (MsCIPK12) were significantly increased under PEG, NaCl, and ABA treatments. Collectively, our study suggested that CBL and CIPK genes play crucial roles in response to various abiotic stresses in Medicago.