Full-length transcript sequencing and comparative transcriptomic analysis to evaluate the contribution of osmotic and ionic stress components towards salinity tolerance in the roots of cultivated alfalfa (Medicago sativa L.)

Thriving in a salty future: morpho-anatomical, physiological and molecular adaptations to salt stress in alfalfa (Medicago sativa L.) and other crops

BACKGROUND

With soil salinity levels rising at an alarming rate, accelerated by climate change and human interventions, there is a growing need for crop varieties that can grow on saline soils. Alfalfa (Medicago sativa) is a cool-season perennial leguminous crop, commonly grown as forage, biofuel feedstock and soil conditioner. It demonstrates significant potential for agricultural circularity and sustainability, for example by fixing nitrogen, sequestering carbon and improving soil structures. Although alfalfa is traditionally regarded as a moderately salt-tolerant species, modern alfalfa varieties display specific salt-tolerance mechanisms, which could be used to pave its role as a leading crop able to grow on saline soils.

SCOPE

Alfalfa's salt tolerance underlies a large variety of cascading biochemical and physiological mechanisms. These are partly enabled by its complex genome structure and out-crossing nature, but which entail impediments for molecular and genetic studies. This review first summarizes the general effects of salinity on plants and the broad-ranging mechanisms for dealing with salt-induced osmotic stress, ion toxicity and secondary stress. Second, we address the defensive and adaptive strategies that have been described for alfalfa, such as the plasticity of alfalfa's root system, hormonal crosstalk for maintaining ion homeostasis, spatiotemporal specialized metabolite profiles and the protection of alfalfa-rhizobia associations. Finally, bottlenecks for research of the physiological and molecular salt-stress responses as well as biotechnology-driven improvements of salt tolerance are identified and discussed.

CONCLUSION

Understanding morpho-anatomical, physiological and molecular responses to salinity is essential for the improvement of alfalfa and other crops in saline land reclamation. This review identifies potential breeding targets for enhancing the stability of alfalfa performance and general crop robustness for rising salt levels as well as to promote alfalfa applications in saline land management.

Exploring the mechanism of seed shattering in Psathyrostachys juncea through histological analysis and comparative transcriptomics

BACKGROUND

Seed shattering (SS) negatively impacts seed yield in Psathyrostachys juncea. Understanding and improving the SS trait requires elucidating the regulatory mechanisms of SS and identifying the key genes involved.

RESULTS

This study presents a comprehensive analysis of the abscission zone (AZ) structures at four developmental stages in two P. juncea genotypes. High-SS P. juncea (H) exhibited a significantly higher SS rate than low-SS P. juncea (L) at all four developmental stages. Anatomical analysis revealed that the degree of lignification in the AZ cell walls is related to the integrity of the abscission structure. The degradation of the AZ in H occurred earlier and was more severe compared to L. At different developmental stages of the AZ, H exhibited higher cellulase and polygalacturonase activities and higher abscisic acid contents compared to L. Conversely, L showed higher lignin, cytokinin, auxin, and gibberellin contents than H. Transcriptomic analysis identified key metabolic pathways related to SS in P. juncea, such as phenylpropanoid biosynthesis, fructose and mannose metabolism, galactose metabolism, and pentose and glucuronate interconversions. The integration of morphological, histological, physiochemical, and metabolic data led to the identification of critical genes, including AUX1, CKX, ABF, GH3, 4CL, CCoAOMT, BGAL, Gal, and PG. The roles of these genes were involved in the regulation of plant hormones and in the synthesis and degradation of cell walls within the AZ.

CONCLUSIONS

This study provides an in-depth understanding of the regulatory mechanisms of SS in P. juncea through comparative transcriptomic analysis. The SS in P. juncea may result from the degradation of the cell wall regulated by cell wall hydrolases genes. The genes identified in this study provide a basis for the genetic improvement of SS traits and serve as a reference for research on other grass species.

Transcriptome Analysis of the Seed Shattering Mechanism in Psathyrostachys juncea Using Full-Length Transcriptome Sequencing

Seed shattering (SS) functions are a survival mechanism in plants, enabling them to withstand adverse environmental conditions and ensure reproduction. However, this trait limits seed yield. Psathyrostachys juncea, a perennial forage grass with many favorable traits, is constrained by SS, limiting its broader application. To investigate the mechanisms underlying SS, second-generation Illumina sequencing and third-generation PacBio sequencing were conducted on abscission zone tissues of P. juncea at 7, 14, 21, and 28 days after heading. GO enrichment analysis identified several significant biological processes, including the "cell wall macromolecule catabolic process", "cell wall polysaccharide catabolic process", "hemicellulose catabolic process", and "xylan catabolic process", all involved in cell wall degradation. KEGG enrichment analysis showed that differentially expressed genes were predominantly enriched in pathways related to "starch and sucrose metabolism", "fructose and mannose metabolism", "phenylpropanoid biosynthesis", "pentose and glucuronate interconversions", and "galactose metabolism", each linked to both the synthesis and degradation of the cell wall. Further analysis of the "starch and sucrose metabolism" pathway revealed genes encoding fructokinase, hexokinase, β-glucosidase, sucrose phosphate synthase, sucrose synthase, and endoglucanase, all of which affected cellulose content. Reduced cellulose content can alter cell wall structure, leading to SS. These findings provide new insights into the regulation of SS in P. juncea and offer valuable references for other species within the Poaceae family.

Accurate Inference of the Polyploid Continuum Using Forward-Time Simulations

Multiple rounds of whole-genome duplication (WGD) followed by diploidization have occurred throughout the evolutionary history of angiosperms. Much work has been done to model the genomic consequences and evolutionary significance of WGD. While researchers have historically modeled polyploids as either allopolyploids or autopolyploids, the variety of natural polyploids span a continuum of differentiation across multiple parameters, such as the extent of polysomic versus disomic inheritance, and the degree of genetic differentiation between the ancestral lineages. Here we present a forward-time polyploid genome evolution simulator called SpecKS. SpecKS models polyploid speciation as originating from a 2D continuum, whose dimensions account for both the level of genetic differentiation between the ancestral parental genomes, as well the time lag between ancestral speciation and their subsequent reunion in the derived polyploid. Using extensive simulations, we demonstrate that changes in initial conditions along either dimension of the 2D continuum deterministically affect the shape of the Ks histogram. Our findings indicate that the error in the common method of estimating WGD time from the Ks histogram peak scales with the degree of allopolyploidy, and we present an alternative, accurate estimation method that is independent of the degree of allopolyploidy. Lastly, we use SpecKS to derive tests that infer both the lag time between parental divergence and WGD time, and the diversity of the ancestral species, from an input Ks histogram. We apply the latter test to transcriptomic data from over 200 species across the plant kingdom, the results of which are concordant with the prevailing theory that the majority of angiosperm lineages are derived from diverse parental genomes and may be of allopolyploid origin.

Genome-Wide Identification and Expression Profiling Analysis of the Mitochondrial Calcium Uniporter Gene Family Under Abiotic Stresses in Medicago sativa

The mitochondrial calcium uniporters (MCUs) are a family of calcium unidirectional transporters important for cytoplasmic Ca2+ signals. Though the MCU proteins in several plant species have been investigated, genome-wide analysis of MCUs in alfalfa is lacking. Here, via genome-wide analysis, a total of 5, 20, and 6 MCU genes were identified in three different alfalfa cultivars, namely Medicago truncatula Jemalong A17, Medicago sativa XinJiangDaYe, and M. sativa Zhongmu No. 1, respectively. They were further phylogenetically classified into three subfamilies. Most MCU genes have only one intron, and gene duplication events of MCU genes were observed within each alfalfa accession and between different accessions. All alfalfa MCU proteins contained a highly conserved MCU domain and 10 conserved motifs, featuring two transmembrane domains and a DI/VME motif. According to the tissue expression data of M. sativa Zhongmu No. 1, MsMCU6.2 was the most abundant transcript with the highest expression in the leaf, and MsMCU5 and MsMCU1.2 showed higher expression levels in the stem than other tissues. We analyzed the expression profiles of five MCU genes (MsMCU1.1/1.2/5/6.1/6.2) under salt, drought, and cold stresses via qRT-PCR assays. All five MCU genes were induced by drought stress, except MsMCU5, whose expression was up-regulated by salt stress, while cold stress slightly altered MsMCU expression. Nine potential interacting proteins and three miRNAs targeting MtMCUs were predicted. These results provide detailed knowledge of alfalfa MCU genes and suggest their potential functions in abiotic stress response.

Differential gene expression of salt-tolerant alfalfa in response to salinity and inoculation by Ensifer meliloti

BACKGROUND

Alfalfa (Medicago sativa L.) experiences many negative effects under salinity stress, which may be mediated by recurrent selection. Salt-tolerant alfalfa may display unique adaptations in association with rhizobium under salt stress.

RESULTS

To elucidate inoculation effects on salt-tolerant alfalfa under salt stress, this study leveraged a salt-tolerant alfalfa population selected through two cycles of recurrent selection under high salt stress. After experiencing 120-day salt stress, mRNA was extracted from 8 random genotypes either grown in 0 or 8 dS/m salt stress with or without inoculation by Ensifer meliloti. Results showed 320 and 176 differentially expressed genes (DEGs) modulated in response to salinity stress or inoculation x salinity stress, respectively. Notable results in plants under 8 dS/m stress included upregulation of a key gene involved in the Target of Rapamycin (TOR) signaling pathway with a concomitant decrease in expression of the SNrK pathway. Inoculation of salt-stressed plants stimulated increased transcription of a sulfate-uptake gene as well as upregulation of the Lysine-27-trimethyltransferase (EZH2), Histone 3 (H3), and argonaute (AGO, a component of miRISC silencing complexes) genes related to epigenetic and post-transcriptional gene control.

CONCLUSIONS

Salt-tolerant alfalfa may benefit from improved activity of TOR and decreased activity of SNrK1 in salt stress, while inoculation by rhizobiumstimulates production of sulfate uptake- and other unique genes.

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.

In silico analysis identified bZIP transcription factors genes responsive to abiotic stress in Alfalfa (Medicago sativa L.)

BACKGROUND

Alfalfa (Medicago sativa L.) is the most cultivated forage legume around the world. Under a variety of growing conditions, forage yield in alfalfa is stymied by biotic and abiotic stresses including heat, salt, drought, and disease. Given the sessile nature of plants, they use strategies including, but not limited to, differential gene expression to respond to environmental cues. Transcription factors control the expression of genes that contribute to or enable tolerance and survival during periods of stress. Basic-leucine zipper (bZIP) transcription factors have been demonstrated to play a critical role in regulating plant growth and development as well as mediate the responses to abiotic stress in several species, including Arabidopsis thaliana, Oryza sativa, Lotus japonicus and Medicago truncatula. However, there is little information about bZIP transcription factors in cultivated alfalfa.

RESULT

In the present study, 237 bZIP genes were identified in alfalfa from publicly available sequencing data. Multiple sequence alignments showed the presence of intact bZIP motifs in the identified sequences. Based on previous phylogenetic analyses in A. thaliana, alfalfa bZIPs were similarly divided and fell into 10 groups. The physico-chemical properties, motif analysis and phylogenetic study of the alfalfa bZIPs revealed high specificity within groups. The differential expression of alfalfa bZIPs in a suite of tissues indicates that bZIP genes are specifically expressed at different developmental stages in alfalfa. Similarly, expression analysis in response to ABA, cold, drought and salt stresses, indicates that a subset of bZIP genes are also differentially expressed and likely play a role in abiotic stress signaling and/or tolerance. RT-qPCR analysis on selected genes further verified these differential expression patterns.

CONCLUSIONS

Taken together, this work provides a framework for the future study of bZIPs in alfalfa and presents candidate bZIPs involved in stress-response signaling.

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.

Insights into the Impact of Trans-Zeatin Overproduction-Engineered Sinorhizobium meliloti on Alfalfa (Medicago sativa L.) Tolerance to Drought Stress

Plant growth-promoting rhizobacteria have been shown to enhance plant tolerance to drought stress through various mechanisms. However, there is limited research on improving drought resistance in alfalfa by genetically modifying PGPR to produce increased levels of cytokinins. Herein, we employed synthetic biology approaches to engineer two novel strains of Sinorhizobium meliloti capable of overproducing trans-Zeatin and investigated their potential in enhancing drought tolerance in alfalfa. Our results demonstrate that alfalfa plants inoculated with these engineered S. meliloti strains exhibited reduced wilting and yellowing while maintaining higher relative water content under drought conditions. The engineered S. meliloti-induced tZ activated the activity of antioxidant enzymes and the accumulation of osmolytes. Additionally, the increased endogenous tZ content in plants alleviated the impact of drought stress on the alfalfa photosynthetic rate. However, under nondrought conditions, inoculation with the engineered S. meliloti strains had no significant effect on alfalfa biomass and nodule formation.

Genome-wide identification of B-box zinc finger (BBX) gene family in Medicago sativa and their roles in abiotic stress responses

BACKGROUND

B-box (BBX) family is a class of zinc finger transcription factors (TFs) that play essential roles in regulating plant growth, development, as well as abiotic stress. However, no systematic analysis of BBX genes has yet been conducted in alfalfa (Medica go sativa L.), and their functions have not been elucidated up to now.

RESULTS

In this study, 28 MsBBX genes were identified from the alfalfa genome, which were clustered into 4 subfamilies according to an evolutionary tree of BBX proteins. Exon-intron structure and conserved motif analysis reflected the evolutionary conservation of MsBBXs in alfalfa. Collinearity analysis showed that segmental duplication promoted the expansion of the MsBBX family. Analysis of cis-regulatory elements suggested that the MsBBX genes possessed many growth/development-, light-, phytohormone-, and abiotic stress-related elements. MsBBX genes were differentially expressed in leaves, flowers, pre-elongated stems, elongated stems, roots and nodules, and most MsBBXs were remarkably induced by drought, salt and various plant growth regulators (ABA, JA, and SA). Further functional verification demonstrated that overexpressing of the MsBBX11 gene clearly promoted salt tolerance in transgenic Arabidopsis by regulating growth and physiological processes of seedlings.

CONCLUSIONS

This research provides insights into further functional research and regulatory mechanisms of MsBBX family genes under abiotic stress of alfalfa.

QiShenYiQi pill inhibits atherosclerosis by promoting TTC39B-LXR mediated reverse cholesterol transport in liver

BACKGROUND

Tetranucleotide repeat domain protein 39B (TTC39B) was found to combine with ubiquitin ligase E3, and promote the ubiquitination modification of liver X receptor (LXR), which led to the inhibition of reverse cholesterol transport and development of atherosclerosis. QiShenYiQi pill (QSYQ) is a modern Chinese patent drug for treating ischemic cardiovascular diseases, the underlying mechanism is found to promote the expression of LXR-α/ ATP-binding cassette transporter G5 (ABCG5) in the liver of atherosclerotic mice.

PURPOSE

The aim of this study is to investigate the effect of QSYQ on TTC39B-LXR mediated reverse cholesterol transport in atherosclerotic mice.

STUDY DESIGN AND METHODS

Male apolipoprotein E gene knockout mice (7 weeks old) were fed with high-fat diet and treated with low dose of QSYQ (QSYQ-l, 0.3 g/kg·d), high dose of QSYQ (QSYQ-H, 1.2 g/kg·d) and LXR-α agonist (LXR-A, GW3965 10 mg/kg·d) for 8 weeks. C57BL/6 J mice were fed with normal diet and used as negative control. Oil red O staining, HE staining, ELISA, RNA sequencing, western blot, immunohistochemistry, RT-PCR, cell culture and RNA interference were performed to analyze the effect of QSYQ on atherosclerosis.

RESULTS

HE staining showed that QSYQ reduced the atherosclerotic lesion significantly when compared to the control group. ELISA measurement showed that QSYQ decreased serum VLDL and increased serum ApoA1. Oil Red O staining showed that QSYQ reduced the lipid content of liver and protect liver function. Comparative transcriptome RNA-sequence of liver showed that DEGs after QSYQ treatment enriched in high-density lipoprotein particle, ubiquitin ligase complex, bile secretion, etc. Immunohistochemical staining and western blot proved that QSYQ increased the protein expression of hepatic SR-B1, LXR-α, LXR-β, CYP7A1 and ABCG5. Targeted inhibiting Ttc39b gene in vitro further established that QSYQ inhibited the gene expression of Ttc39b, increased the protein expression of SR-B1, LXR-α/β, CYP7A1 and ABCG5 in rat hepatocyte.

CONCLUSION

Our results demonstrated the new anti-atherosclerotic mechanism of QSYQ by targeting TTC39B-LXR mediated reverse cholesterol transport in liver. QSYQ not only promoted reverse cholesterol transport, but also improved fatty liver and protected liver function.

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/.

Physiological, Epigenetic, and Transcriptome Analyses Provide Insights into the Responses of Wheat Seedling Leaves to Different Water Depths under Flooding Conditions

Flooding stress, including waterlogging and submergence, is one of the major abiotic stresses that seriously affects the growth and development of plants. In the present study, physiological, epigenetic, and transcriptomic analyses were performed in wheat seedling leaves under waterlogging (WL), half submergence (HS), and full submergence (FS) treatments. The results demonstrate that FS increased the leaves' hydrogen peroxide (H2O2) and malondialdehyde (MDA) contents and reduced their chlorophyll contents (SPAD), photosynthetic efficiency (Fv/Fm), and shoot dry weight more than HS and WL. In addition, FS increased catalase (CAT) and peroxidase (POD) activities more than HS and WL. However, there were no significant differences in the contents of H2O2, MDA, SPAD, and Fv/Fm, and the activities of superoxide dismutase (SOD) and POD between the HS and WL treatments. The changes in DNA methylation were related to stress types, increasing under the WL and HS treatments and decreasing under the FS treatment. Additionally, a total of 9996, 10,619, and 24,949 genes were differentially expressed under the WL, HS, and FS treatments, respectively, among which the 'photosynthesis', 'phenylpropanoid biosynthesis', and 'plant hormone signal transduction' pathways were extensively enriched under the three flooding treatments. The genes involved in these pathways showed flooding-type-specific expression. Moreover, flooding-type-specific responses were observed in the three conditions, including the enrichment of specific TFs and response pathways. These results will contribute to a better understanding of the molecular mechanisms underlying the responses of wheat seedling leaves to flooding stress and provide valuable genetic and epigenetic information for breeding flood-tolerant varieties of wheat.

Breeding Efficiency for Salt Tolerance in Alfalfa

Alfalfa (Medicago sativa L.), one of the most extensively grown forage crops, is sensitive to saline soils. We measured the breeding efficiency for increased salt tolerance in alfalfa by comparing lines selected from BC79S, CS, and SII populations with their unselected parental means for forage mass and associated changes in stem length, leaf-to-stem ratio (LSR), number of nodes per stem, crude protein (CP) content, and neutral detergent fiber (NDF) content. The overall forage mass in the non-salt-stressed test (9562 kg ha-1) was greater (p < 0.001) than under salt stress (5783 kg ha-1), with a 40% production advantage. In the non-salt-stressed test, the BC79S and CS lines averaged at a 4% lower production than their parents, while SII lines had on average a 9% greater production. Conversely, in the salt-stressed test, all lines showed a 20% overall greater seasonal production than their parents. Some selected lines produced more forage mass in both the non-stressed and salt-stressed tests than their parents. The stem length, LSR, node number, CP content, and NDF content of the selected lines varied with respect to non-stressed vs. stressed, but they tended not to differ greatly from their respective parental means under either non- or salt-stressed conditions. The selection protocol provided a universal increase in forage mass under salt-stressed field conditions of the selected lines. Furthermore, we identified lines with forage mass values greater than their parental means under non- and salt-stressed field conditions.

Genome-wide investigation and expression analysis of OSCA gene family in response to abiotic stress in alfalfa

Alfalfa is an excellent leguminous forage crop that is widely cultivated worldwide, but its yield and quality are often affected by drought and soil salinization. Hyperosmolality-gated calcium-permeable channel (OSCA) proteins are hyperosmotic calcium ion (Ca2+) receptors that play an essential role in regulating plant growth, development, and abiotic stress responses. However, no systematic analysis of the OSCA gene family has been conducted in alfalfa. In this study, a total of 14 OSCA genes were identified from the alfalfa genome and classified into three groups based on their sequence composition and phylogenetic relationships. Gene structure, conserved motifs and functional domain prediction showed that all MsOSCA genes had the same functional domain DUF221. Cis-acting element analysis showed that MsOSCA genes had many cis-regulatory elements in response to abiotic or biotic stresses and hormones. Tissue expression pattern analysis demonstrated that the MsOSCA genes had tissue-specific expression; for example, MsOSCA12 was only expressed in roots and leaves but not in stem and petiole tissues. Furthermore, RT-qPCR results indicated that the expression of MsOSCA genes was induced by abiotic stress (drought and salt) and hormones (JA, SA, and ABA). In particular, the expression levels of MsOSCA3, MsOSCA5, MsOSCA12 and MsOSCA13 were significantly increased under drought and salt stress, and MsOSCA7, MsOSCA10, MsOSCA12 and MsOSCA13 genes exhibited significant upregulation under plant hormone treatments, indicating that these genes play a positive role in drought, salt and hormone responses. Subcellular localization results showed that the MsOSCA3 protein was localized on the plasma membrane. This study provides a basis for understanding the biological information and further functional analysis of the MsOSCA gene family and provides candidate genes for stress resistance breeding in alfalfa.

The physiology of drought stress in two grapevine cultivars: Photosynthesis, antioxidant system, and osmotic regulation responses

Drought stress impedes viticultural plant growth and development by modifying various metabolic pathways. However, the regulatory network response underlying drought stress is not yet clear. In this study, the leaves and roots of "Shine Muscat" ("SM," Vitis labruscana × Vitis vinifera) and "Thompson Seedless" ("TS," V. vinifera L. cv.) were subjected to drought stress to study the regulatory network used by drought stress. Morphophysiological results showed that the malondialdehyde content after 28 days of drought stress increased more significantly in "TS" than "SM." Furthermore, the multiomics analysis studies showed that a total of 3036-6714 differentially expressed genes and 379-385 differentially abundant metabolites were identified in "SM" and "TS" grapevine cultivars under drought stress. Furthermore, the retained intron was the major form of differential alternative splicing event under drought stress. The photosynthesis pathway, antioxidant system, plant hormone signal transduction, and osmotic adjustment were the primary response systems in the two grapevine cultivars under drought stress. We have identified GRIK1, RFS2, and LKR/SDH as the hub genes in the coexpression network of drought stress. In addition, the difference in the accumulation of pheophorbide-a reveals different drought resistance mechanisms in the two grapevine cultivars. Our study explained the difference in drought response between cultivars and tissues and identified drought stress-responsive genes, which provides reference data for further understanding the regulatory network of drought tolerance in grapevine.

Genome-wide identification and expression analysis of the Auxin-Response factor (ARF) gene family in Medicago sativa under abiotic stress

BACKGROUND

Alfalfa (Medicago sativa) is the most widely planted legume forage and one of the most economically valuable crops in the world. The periodic changes in its growth and development and abiotic stress determine its yield and economic benefits. Auxin controls many aspects of alfalfa growth by regulating gene expression, including organ differentiation and stress response. Auxin response factors (ARF) are transcription factors that play an essential role in auxin signal transduction and regulate the expression of auxin-responsive genes. However, the function of ARF transcription factors is unclear in autotetraploid-cultivated alfalfa.

RESULT

A total of 81 ARF were identified in the alfalfa genome in this study. Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were analyzed, identifying that ARF genes are mainly involved in transcriptional regulation and plant hormone signal transduction pathways. Phylogenetic analysis showed that MsARF was divided into four clades: I, II, III, and IV, each containing 52, 13, 7, and 9 genes, respectively. The promoter region of the MsARF gene contained stress-related elements, such as ABRE, TC-rich repeats, MBS, LTR. Proteins encoded by 50 ARF genes were localized in the nucleus without guide peptides, signal peptides, or transmembrane structures, indicating that most MsARF genes are not secreted or transported but only function in the nucleus. Protein structure analysis revealed that the secondary and tertiary structures of the 81 MsARF genes varied. Chromosomal localization analysis showed 81 MsARF genes were unevenly distributed on 25 chromosomes, with the highest distribution on chromosome 5. Furthermore, 14 segmental duplications and two sets of tandem repeats were identified. Expression analysis indicated that the MsARF was differentially expressed in different tissues and under various abiotic stressors. The quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis showed that the expression profiles of 23 MsARF genes were specific to abiotic stresses such as drought, salt, high temperature, and low temperature, as well as tissue-specific and closely related to the duration of stress.

CONCLUSION

This study identified MsARF in the cultivated alfalfa genome based on the autotetraploid level, which GO, KEGG analysis, phylogenetic analysis, sequence characteristics, and expression pattern analysis further confirmed. Together, these findings provide clues for further investigation of MsARF functional verification and molecular breeding of alfalfa. This study provides a novel approach to systematically identify and characterize ARF transcription factors in autotetraploid cultivated alfalfa, revealing 23 MsARF genes significantly involved in response to various stresses.

Identification and evolution analysis of YUCCA genes of Medicago sativa and Medicago truncatula and their expression profiles under abiotic stress

The YUCCAs (YUC) are functionally identified flavin-containing monooxidases (FMOs) in plants that act as an important rate-limiting enzyme functioning in the auxin synthesis IPA (indole-3-pyruvic acid) pathway. In this study, 12 MsYUCs and 15 MtYUCs containing characteristic conserved motifs were identified in M. sativa (Medicago sativa L.) and M. truncatula (Medicago truncatula Gaertn.), respectively. Phylogenetic analysis revealed that YUC proteins underwent an evolutionary divergence. Both tandem and segmental duplication events were presented in MsYUC and MtYUC genes. Comparative syntenic maps of M. sativa with M. truncatula, Arabidopsis (Arabidopsis thaliana), or rice (Oryza sativa L.) were constructed to illustrate the evolution relationship of the YUC gene family. A large number of cis-acting elements related to stress response and hormone regulation were revealed in the promoter sequences of MsYUCs. Expression analysis showed that MsYUCs had a tissue-specific, genotype-differential expression and a differential abiotic stress response pattern based on transcriptome data analysis of M. sativa online. In addition, RT-qPCR confirmed that salt stress significantly induced the expression of MsYUC1/MsYUC10 but significantly inhibited MsYUC2/MsYUC3 expression and the expression of MsYUC10/MsYUC11/MsYUC12 was significantly induced by cold treatment. These results could provide valuable information for functional analysis of YUC genes via gene engineering of the auxin synthetic IPA pathway in Medicago.

Optimizing an efficient ensemble approach for high-quality de novo transcriptome assembly of Thymus daenensis

Non-erroneous and well-optimized transcriptome assembly is a crucial prerequisite for authentic downstream analyses. Each de novo assembler has its own algorithm-dependent pros and cons to handle the assembly issues and should be specifically tested for each dataset. Here, we examined efficiency of seven state-of-art assemblers on ~ 30 Gb data obtained from mRNA-sequencing of Thymus daenensis. In an ensemble workflow, combining the outputs of different assemblers associated with an additional redundancy-reducing step could generate an optimized outcome in terms of completeness, annotatability, and ORF richness. Based on the normalized scores of 16 benchmarking metrics, EvidentialGene, BinPacker, Trinity, rnaSPAdes, CAP3, IDBA-trans, and Velvet-Oases performed better, respectively. EvidentialGene, as the best assembler, totally produced 316,786 transcripts, of which 235,730 (74%) were predicted to have a unique protein hit (on uniref100), and also half of its transcripts contained an ORF. The total number of unique BLAST hits for EvidentialGene was approximately three times greater than that of the worst assembler (Velvet-Oases). EvidentialGene could even capture 17% and 7% more average BLAST hits than BinPacker and Trinity. Although BinPacker and CAP3 produced longer transcripts, the EvidentialGene showed a higher collinearity between transcript size and ORF length. Compared with the other programs, EvidentialGene yielded a higher number of optimal transcript sets, further full-length transcripts, and lower possible misassemblies. Our finding corroborates that in non-model species, relying on a single assembler may not give an entirely satisfactory result. Therefore, this study proposes an ensemble approach of accompanying EvidentialGene pipelines to acquire a superior assembly for T. daenensis.

Elucidation of Physiological, Transcriptomic and Metabolomic Salinity Response Mechanisms in Medicago sativa

Alfalfa (Medicago sativa L.) is a widely grown perennial leguminous forage crop with a number of positive attributes. However, despite its moderate ability to tolerate saline soils, which are increasing in prevalence worldwide, it suffers considerable yield declines under these growth conditions. While a general framework of the cascade of events involved in plant salinity response has been unraveled in recent years, many gaps remain in our understanding of the precise molecular mechanisms involved in this process, particularly in non-model yet economically important species such as alfalfa. Therefore, as a means of further elucidating salinity response mechanisms in this species, we carried out in-depth physiological assessments of M. sativa cv. Beaver, as well as transcriptomic and untargeted metabolomic evaluations of leaf tissues, following extended exposure to salinity (grown for 3-4 weeks under saline treatment) and control conditions. In addition to the substantial growth and photosynthetic reductions observed under salinity treatment, we identified 1233 significant differentially expressed genes between growth conditions, as well as 60 annotated differentially accumulated metabolites. Taken together, our results suggest that changes to cell membranes and walls, cuticular and/or epicuticular waxes, osmoprotectant levels, antioxidant-related metabolic pathways, and the expression of genes encoding ion transporters, protective proteins, and transcription factors are likely involved in alfalfa's salinity response process. Although some of these alterations may contribute to alfalfa's modest salinity resilience, it is feasible that several may be disadvantageous in this context and could therefore provide valuable targets for the further improvement of tolerance to this stress in the future.

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

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

Physiological and transcriptomic analyses provide insight into thermotolerance in desert plant Zygophyllum xanthoxylum

BACKGROUND

Heat stress has adverse effects on the growth and reproduction of plants. Zygophyllum xanthoxylum, a typical xerophyte, is a dominant species in the desert where summer temperatures are around 40 °C. However, the mechanism underlying the thermotolerance of Z. xanthoxylum remained unclear.

RESULTS

Here, we characterized the acclimation of Z. xanthoxylum to heat using a combination of physiological measurements and transcriptional profiles under treatments at 40 °C and 45 °C, respectively. Strikingly, moderate high temperature (40 °C) led to an increase in photosynthetic capacity and superior plant performance, whereas severe high temperature (45 °C) was accompanied by reduced photosynthetic capacity and inhibited growth. Transcriptome profiling indicated that the differentially expressed genes (DEGs) were related to transcription factor activity, protein folding and photosynthesis under heat conditions. Furthermore, numerous genes encoding heat transcription shock factors (HSFs) and heat shock proteins (HSPs) were significantly up-regulated under heat treatments, which were correlated with thermotolerance of Z. xanthoxylum. Interestingly, the up-regulation of PSI and PSII genes and the down-regulation of chlorophyll catabolism genes likely contribute to improving plant performance of Z. xanthoxylum under moderate high temperature.

CONCLUSIONS

We identified key genes associated with of thermotolerance and growth in Z. xanthoxylum, which provide significant insights into the regulatory mechanisms of thermotolerance and growth regulation in Z. xanthoxylum under high temperature conditions.

Full-length transcriptome sequencing and comparative transcriptome analysis of Eriocheir sinensis in response to infection by the microsporidian Hepatospora eriocheir

As a new generation of high-throughput sequencing technology, PacBio Iso-Seq technology (Iso-Seq) provides a better alternative sequencing method for the acquisition of full-length unigenes. In this study, a total of 22.27 gigabyte (Gb) subread bases and 128,614 non-redundant unigenes (mean length: 2,324 bp) were obtained from six main tissues of Eriocheir sinensis including the heart, nerve, intestine, muscle, gills and hepatopancreas. In addition, 74,732 unigenes were mapped to at least one of the following databases: Non-Redundant Protein Sequence Database (NR), Gene Ontology (GO), Kyoto Encyclopaedia of Genes and Genomes (KEGG), KEGG Orthology (KO) and Protein family (Pfam). In addition, 6696 transcription factors (TFs), 28,458 long non-coding RNAs (lncRNAs) and 94,230 mRNA-miRNA pairs were identified. Hepatospora eriocheir is the primary pathogen of E. sinensis and can cause hepatopancreatic necrosis disease (HPND); the intestine is the main target tissue. Here, we attempted to identify the key genes related to H. eriocheir infection in the intestines of E. sinensis. By combining Iso-Seq and Illumina RNA-seq analysis, we identified a total of 12,708 differentially expressed unigenes (DEUs; 6,696 upregulated and 6,012 downregulated) in the crab intestine following infection with H. eriocheir. Based on the biological analysis of these DEUs, several key processes were identified, including energy metabolism-related pathways, cell apoptosis and innate immune-related pathways. Twelve selected genes from these DEUs were subsequently verified by quantitative real-time PCR (qRT-PCR) analysis. Our findings enhance our understanding of the E. sinensis transcriptome and the specific association between E. sinensis and H. eriocheir infection.

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.

Evolution of the CBL and CIPK gene families in Medicago: genome-wide characterization, pervasive duplication, and expression pattern under salt and drought stress

BACKGROUND

Calcineurin B-like proteins (CBLs) are ubiquitous Ca²⁺ sensors that mediate plant responses to various stress and developmental processes by interacting with CBL-interacting protein kinases (CIPKs). CBLs and CIPKs play essential roles in acclimatization of crop plants. However, evolution of these two gene families in the genus Medicago is poorly understood.

RESULTS

A total of 68 CBL and 135 CIPK genes have been identified in five genomes from Medicago. Among these genomes, the gene number of CBLs and CIPKs shows no significant difference at the haploid genome level. Phylogenetic and comprehensive characteristic analyses reveal that CBLs and CIPKs are classified into four clades respectively, which is validated by distribution of conserved motifs. The synteny analysis indicates that the whole genome duplication events (WGDs) have contributed to the expansion of both families. Expression analysis demonstrates that two MsCBLs and three MsCIPKs are specifically expressed in roots, mature leaves, developing flowers and nitrogen fixing nodules of Medicago sativa spp. sativa, the widely grown tetraploid species. In particular, the expression of these five genes was highly up-regulated in roots when exposed to salt and drought stress, indicating crucial roles in stress responses.

CONCLUSIONS

Our study leads to a comprehensive understanding of evolution of CBL and CIPK gene families in Medicago, but also provides a rich resource to further address the functions of CBL-CIPK complexes in cultivated species and their closely related wild relatives.

OMICS in Fodder Crops: Applications, Challenges, and Prospects

Biomass yield and quality are the primary targets in forage crop improvement programs worldwide. Low-quality fodder reduces the quality of dairy products and affects cattle's health. In multipurpose crops, such as maize, sorghum, cowpea, alfalfa, and oat, a plethora of morphological and biochemical/nutritional quality studies have been conducted. However, the overall growth in fodder quality improvement is not on par with cereals or major food crops. The use of advanced technologies, such as multi-omics, has increased crop improvement programs manyfold. Traits such as stay-green, the number of tillers per plant, total biomass, and tolerance to biotic and/or abiotic stresses can be targeted in fodder crop improvement programs. Omic technologies, namely genomics, transcriptomics, proteomics, metabolomics, and phenomics, provide an efficient way to develop better cultivars. There is an abundance of scope for fodder quality improvement by improving the forage nutrition quality, edible quality, and digestibility. The present review includes a brief description of the established omics technologies for five major fodder crops, i.e., sorghum, cowpea, maize, oats, and alfalfa. Additionally, current improvements and future perspectives have been highlighted.

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.

NUCLEAR TRANSPORT FACTOR 2-LIKE improves drought tolerance by modulating leaf water loss in alfalfa (Medicago sativa L.)

Drought is a major environmental factor that limits the production of alfalfa (Medicago sativa). In the present study, M. sativa NUCLEAR TRANSPORT FACTOR 2-LIKE (MsNTF2L) was identified as a nucleus-, cytoplasm-, and plasma membrane-localized protein. Its transcriptional expression was highly induced by ABA and drought stress. Overexpression of MsNTF2L in Arabidopsis resulted in hypersensitivity to ABA during both the seed germination and seedling growth stages. However, transgenic Arabidopsis plants exhibited enhanced tolerance to drought stress by reducing the levels of reactive oxygen species (ROS) and increasing the expression of stress/ABA-inducible genes. Consistently, analysis of MsNTF2L overexpression (OE) and RNA interference (RNAi) alfalfa plants revealed that MsNTF2L confers drought tolerance through promoting ROS scavenging, a decrease in stomatal density, ABA-induced stomatal closure, and epicuticular wax crystal accumulation. MsNTF2L highly affected epicuticular wax deposition, as a large group of wax biosynthesis and transport genes were influenced in the alfalfa OE and RNAi lines. Furthermore, transcript profiling of drought-treated alfalfa WT, OE, and RNAi plants showed a differential drought response for genes related to stress/ABA signaling, antioxidant defense, and photosynthesis. Taken together, these results reveal that MsNTF2L confers drought tolerance in alfalfa via modulation of leaf water loss (by regulating both stomata and wax deposition), antioxidant defense, and photosynthesis.

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

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

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.

Comparative Physiological and Transcriptome Profiles Uncover Salt Tolerance Mechanisms in Alfalfa

Salinity is a major limiting factor that affects crop production. Understanding of the mechanisms of plant salt tolerance is critical for improving crop yield on saline land. Alfalfa (Medicago sativa L.) is the most important forage crop, while its salt tolerance mechanisms are largely unknown. The physiological and transcriptomic responses in two contrasting salt tolerant cultivars to salinity stress were investigated in the present study. "Magnum Salt" showed higher salt tolerance than "Adrenalin," with higher relative germination rate, survival rate, biomass and K⁺/Na⁺ ratio after salt treatment. Activities of antioxidant enzymes SOD, CAT and GR, and proline concentrations were upregulated to higher levels in roots and shoots in Magnum Salt than in Adrenalin after salinity stress, except for no difference in GR activity in shoots, and lower levels of O₂ ^⋅- and H₂O₂ were accumulated in leaves. It was interesting to find that salinity caused a decrease in total unsaturated fatty acid in Adrenalin other than Magnum Salt, C18:2 was increased significantly after salinity in Magnum Salt, while it was unaltered in Adrenalin. High quality RNA sequencing (RNA-seq) data was obtained from samples of Magnum Salt and Adrenalin at different time points (0, 2, and 26 h). Generally, "phagosome," "TCA cycle" and "oxidative phosphorylation" pathways were inhibited by salinity stress. Upregulated DEGs in Magnum Salt were specifically enriched in "fatty acid metabolism," "MAPK signaling" and "hormone signal transduction" pathways. The DEGs involved in ionic homeostasis, reactive oxygen species (ROS) scavenging and fatty acid metabolism could partially explain the difference in salt tolerance between two cultivars. It is suggested that salt tolerance in alfalfa is associated with regulation of ionic homeostasis, antioxidative enzymes and fatty acid metabolism at both transcriptional and physiological level.

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 Ca²⁺ 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.

The Isolation and Full-Length Transcriptome Sequencing of a Novel Nidovirus and Response of Its Infection in Japanese Flounder (Paralichthys olivaceus)

A novel nidovirus, CSBV Bces-Po19, was isolated from the marine fish, Japanese flounder (Paralichthys olivaceus). The viral genome was 26,597 nucleotides long and shared 98.62% nucleotide identity with CSBV WHQSR4345. PacBio Sequel and Illumina sequencing were used to perform full-length transcriptome sequencing on CSBV Bces-Po19-sensitive (S) and -resistant (R) Japanese flounder. The results of negative staining revealed bacilliform and spherical virions. There were in total 1444 different genes between CSBV Bces-Po19 S and R groups, with 935 being up-regulated and 513 being down-regulated. Metabolism-, immune-, and RNA-related pathways were significantly enriched. Furthermore, CSBV Bces-Po19 infection induced alternative splicing (AS) events in Japanese flounder; the S group had a higher numbers of AS events (12,352) than the R group (11,452). The number of long non-coding RNA (lncRNA) in the S group, on the other hand, was significantly lower than in the R group. In addition to providing valuable information that sheds more light on CSBV Bces-Po19 infection, these research findings provide further clues for CSBV Bces-Po19 prevention and treatment.

Lb1G04202, an Uncharacterized Protein from Recretohalophyte Limonium bicolor, Is Important in Salt Tolerance

With global increases in saline soil, it has become increasingly important to decipher salt-tolerance mechanisms and identify strategies to improve salt tolerance in crops. Halophytes complete their life cycles in environments containing ≥200 mM NaCl; these remarkable plants provide a potential source of genes for improving crop salt tolerance. Recretohalophytes such as Limonium bicolor have salt glands that secrete Na+ on their leaf epidermis. Here, we identified Lb1G04202, an uncharacterized gene with no conserved domains, from L. bicolor, which was highly expressed after NaCl treatment. We confirmed its expression in the salt gland by in situ hybridization, and then heterologously expressed Lb1G04202 in Arabidopsis thaliana. The transgenic lines had a higher germination rate, greater cotyledon growth percentage, and longer roots than the wild type (WT) under NaCl treatments (50, 100 and 150 mM). At the seedling stage, the transgenic lines grew better than the WT and had lower Na+ and malonyldialdehyde accumulation, and higher K+ and proline contents. This corresponded with the high expression of the key proline biosynthesis genes AtP5CS1 and AtP5CS2 under NaCl treatment. Isotonic mannitol treatment showed that Lb1G04202 overexpression significantly relieved osmotic stress. Therefore, this novel gene provides a potential target for improving salt tolerance.

The Mutational, Epigenetic, and Transcriptional Effects Between Mixed High-Energy Particle Field (CR) and 7Li-Ion Beams (LR) Radiation in Wheat M1 Seedlings

Ionizing radiation (IR) is an effective approach for mutation breeding. Understanding the mutagenesis and transcriptional profiles induced by different mutagens is of great significance for improving mutation breeding efficiency. Here, using RNA sequencing and methylation-sensitive amplification polymorphism (MSAP) approaches, we compared the genetic variations, epigenetics, and transcriptional responses induced by the mixed high-energy particle field (CR) and ⁷Li-ion beam (LR) radiation in M₁ seedlings of two wheat genotypes (Yangmai 18 and Yangmai 20). The results showed that, in both wheat genotypes, CR displayed significantly a higher mutation efficiency (1.79 × 10^-6/bp) than that by LR (1.56 × 10^-6/bp). The induced mutations were not evenly distributed across chromosomes and varied across wheat genotypes. In Y18 M₁, the highest number of mutations were detected on Chr. 6B and Chr. 6D, whilst in Y20 M₁, Chr. 7A and Chr. 3A had the highest mutations. The transcript results showed that total of 4,755 CR-regulated and 1,054 LR-regulated differentially expressed genes (DEGs) were identified in the both genotypes. Gene function enrichment analysis of DEGs showed that these DEGs overlapped or diverged in the cascades of molecular networks involved in "phenylpropanoid biosynthesis" and "starch and sucrose metabolism" pathways. Moreover, IR type specific responses were observed between CR an LR irradiation, including specific TFs and response pathways. MSAP analysis showed that DNA methylation level increased in LR treatment, while decreased at CR. The proportion of hypermethylation was higher than that of hypomethylation at LR, whereas a reverse pattern was observed at CR, indicating that DNA methylation plays critical roles in response to IR irradiation. All these results support that the response to different IRs in wheat includes both common and unique pathways, which can be served as a useful resource to better understand the mechanisms of responses to different IRs in other plants.

SMRT and Illumina RNA sequencing reveal the complexity of terpenoid biosynthesis in Zanthoxylum armatum

Sichuan pepper (Zanthoxylum armatum DC) is a popular spice and is often prescribed in traditional Chinese medicine to treat vomiting, diarrhea, ascariasis and eczema, among other conditions. Volatile oils from Z. armatum leaves contain active ingredients, with terpenoids being one of the main components. In the present study, the combination of sequencing data of Z. armatum from PacBio single molecule real time (SMRT) and Illumina RNA sequencing (RNA-Seq) platforms facilitated an understanding of the gene regulatory network of terpenoid biosynthesis in pepper leaves. The leaves of three developmental stages from two Z. armatum cultivars, 'Rongchangwuci' (WC) and 'Zhuye' (ZY), were selected as test materials to construct sequencing libraries. A total of 143,122 predictions of unique coding sequences, 105,465 simple sequence repeats, 20,145 transcription factors and 4719 long non-coding RNAs (lncRNAs) were identified, and 142,829 transcripts were successfully annotated. The occurrence of alternative splicing events was verified by reverse transcription PCR, and quantitative real-time PCR was used to confirm the expression pattern of six randomly selected lncRNAs. A total of 96,931 differentially expressed genes were filtered from different samples. According to functional annotation, a total of 560 candidate genes were involved in terpenoid synthesis, of which 526 were differentially expressed genes (DEGs). To identify the key genes involved in terpenoid biosynthesis, the module genes in different samples, including structural and transcription factors genes, were analyzed using the weighted gene co-expression network method, and the co-expression network of genes was constructed. Thirty-one terpenoids were identified by gas chromatography-mass spectrometry. The correlation between 18 compounds with significantly different contents and genes with high connectivity in the module was jointly analyzed in both cultivars, yielding 12 candidate DEGs presumably involved in the regulation of terpenoid biosynthesis. Our findings showed that full-length transcriptome SMRT and Illumina RNA-Seq can play an important role in studying organisms without reference genomes and elucidating the gene regulation of a biosynthetic pathway.

How salt stress-responsive proteins regulate plant adaptation to saline conditions

An overview is presented of recent advances in our knowledge of candidate proteins that regulate various physiological and biochemical processes underpinning plant adaptation to saline conditions. Salt stress is one of the environmental constraints that restrict plant distribution, growth and yield in many parts of the world. Increased world population surely elevates food demands all over the globe, which anticipates to add a great challenge to humanity. These concerns have necessitated the scientists to understand and unmask the puzzle of plant salt tolerance mechanisms in order to utilize various strategies to develop salt tolerant crop plants. Salt tolerance is a complex trait involving alterations in physiological, biochemical, and molecular processes. These alterations are a result of genomic and proteomic complement readjustments that lead to tolerance mechanisms. Proteomics is a crucial molecular tool that indicates proteins expressed by the genome, and also identifies the functions of proteins accumulated in response to salt stress. Recently, proteomic studies have shed more light on a range of promising candidate proteins that regulate various processes rendering salt tolerance to plants. These proteins have been shown to be involved in photosynthesis and energy metabolism, ion homeostasis, gene transcription and protein biosynthesis, compatible solute production, hormone modulation, cell wall structure modification, cellular detoxification, membrane stabilization, and signal transduction. These candidate salt responsive proteins can be therefore used in biotechnological approaches to improve tolerance of crop plants to salt conditions. In this review, we provided comprehensive updated information on the proteomic data of plants/genotypes contrasting in salt tolerance in response to salt stress. The roles of salt responsive proteins that are potential determinants for plant salt adaptation are discussed. The relationship between changes in proteome composition and abundance, and alterations observed in physiological and biochemical features associated with salt tolerance are also addressed.

Transcriptome Analysis of Salt Stress in Hibiscus hamabo Sieb. et Zucc Based on Pacbio Full-Length Transcriptome Sequencing

Hibiscus hamabo Sieb. et Zucc is an important semi-mangrove plant with great morphological features and strong salt resistance. In this study, by combining single molecule real time and next-generation sequencing technologies, we explored the transcriptomic changes in the roots of salt stressed H. hamabo. A total of 94,562 unigenes were obtained by clustering the same isoforms using the PacBio RSII platform, and 2269 differentially expressed genes were obtained under salt stress using the Illumina platform. There were 519 differentially expressed genes co-expressed at each treatment time point under salt stress, and these genes were found to be enriched in ion signal transduction and plant hormone signal transduction. We used Arabidopsis thaliana (L.) Heynh. transformation to confirm the function of the HhWRKY79 gene and discovered that overexpression enhanced salt tolerance. The full-length transcripts generated in this study provide a full characterization of the transcriptome of H. hamabo and may be useful in mining new salt stress-related genes specific to this species, while facilitating the understanding of the salt tolerance mechanisms.

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.

Physiological and transcriptomic analyses reveal novel insights into the cultivar-specific response to alkaline stress in alfalfa (Medicago sativa L.)

Soil alkalization severely limits plant growth and development, however, the mechanisms of alkaline response in plants remain largely unknown. In this study, we performed physiological and transcriptomic analyses using two alfalfa cultivars (Medicago sativa L.) with different sensitivities to alkaline conditions. The chlorophyll content and shoot fresh mass drastically declined in the alkaline-sensitive cultivar Algonquin (AG) following alkaline treatment (0-25 mM Na2CO3 solution), while the alkaline-tolerant cultivar Gongnong NO.1 (GN) maintained relatively stable growth and chlorophyll content. Compared with AG, GN had higher contents of Ca2+ and Mg2+; the ratios of Ca2+ and Mg2+ to Na+, proline and soluble sugar, as well as higher enzyme activities of peroxidase (POD) and catalase (CAT) under the alkaline conditions. Furthermore, transcriptomic analysis identified three categories of alkaline-responsive differentially expressed genes (DEGs) between the two cultivars: 48 genes commonly induced in both the cultivars (CAR), 574 genes from the tolerant cultivar (TAR), and 493 genes from the sensitive cultivar (SAR). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that CAR genes were mostly involved in phenylpropanoid biosynthesis, lipid metabolism, and DNA replication and repair; TAR genes were significantly enriched in metabolic pathways, such as biosynthesis of amino acids and secondary metabolites including flavonoids, and the MAPK signaling pathway; SAR genes were specifically enriched in vitamin B6 metabolism. Taken together, the results identified candidate pathways associated with genetic variation in response to alkaline stress, providing novel insights into the mechanisms underlying alkaline tolerance in alfalfa.

Transcriptomic analysis of differentially expressed genes in leaves and roots of two alfalfa (Medicago sativa L.) cultivars with different salt tolerance

BACKGROUND

Alfalfa (Medicago sativa L.) production decreases under salt stress. Identification of genes associated with salt tolerance in alfalfa is essential for the development of molecular markers used for breeding and genetic improvement.

RESULT

An RNA-Seq technique was applied to identify the differentially expressed genes (DEGs) associated with salt stress in two alfalfa cultivars: salt tolerant 'Halo' and salt intolerant 'Vernal'. Leaf and root tissues were sampled for RNA extraction at 0 h, 3 h, and 27 h under 12 dS m^- 1 salt stress maintained by NaCl. The sequencing generated a total of 381 million clean sequence reads and 84.8% were mapped on to the alfalfa reference genome. A total of 237 DEGs were identified in leaves and 295 DEGs in roots of the two alfalfa cultivars. In leaf tissue, the two cultivars had a similar number of DEGs at 3 h and 27 h of salt stress, with 31 and 49 DEGs for 'Halo', 34 and 50 for 'Vernal', respectively. In root tissue, 'Halo' maintained 55 and 56 DEGs at 3 h and 27 h, respectively, while the number of DEGs decreased from 42 to 10 for 'Vernal'. This differential expression pattern highlights different genetic responses of the two cultivars to salt stress at different time points. Interestingly, 28 (leaf) and 31 (root) salt responsive candidate genes were highly expressed in 'Halo' compared to 'Vernal' under salt stress, of which 13 candidate genes were common for leaf and root tissues. About 60% of DEGs were assigned to known gene ontology (GO) categories. The genes were involved in transmembrane protein function, photosynthesis, carbohydrate metabolism, defense against oxidative damage, cell wall modification and protection against lipid peroxidation. Ion binding was found to be a key molecular activity for salt tolerance in alfalfa under salt stress.

CONCLUSION

The identified DEGs are significant for understanding the genetic basis of salt tolerance in alfalfa. The generated genomic information is useful for molecular marker development for alfalfa genetic improvement for salt tolerance.

Comparative analysis of alfalfa (Medicago sativa L.) seedling transcriptomes reveals genotype-specific drought tolerance mechanisms

Drought is one of the main abiotic factors that affect alfalfa yield. The identification of genes that control this complex trait can provide important insights for alfalfa breeding. However, little is known about how alfalfa responds and adapts to drought stress, particularly in cultivars of differing drought tolerance. In this study, the drought-tolerant cultivar Dryland 'DT' and the drought-sensitive cultivar WL343HQ 'DS' were used to characterize leaf and root physiological responses and transcriptional changes in response to water deficit. Under drought stress, Dryland roots (DTR) showed more differentially expressed genes than WL343HQ roots (DSR), whereas WL343HQ leaves (DSL) showed more differentially expressed genes than Dryland leaves (DTL). Many of these genes were involved in stress-related pathways, carbohydrate metabolism, and lignin and wax biosynthesis, which may have improved the drought tolerance of alfalfa. We also observed that several genes related to ABA metabolism, root elongation, peroxidase activity, cell membrane stability, ubiquitination, and genetic processing responded to drought stress in alfalfa. We highlighted several candidate genes, including sucrose synthase, xylan 1,4-beta-xylosidase, primary-amine oxidase, and alcohol-forming fatty acyl-CoA reductase, for future studies on drought stress resistance in alfalfa and other plant species. In summary, our results reveal the unique drought adaptation and resistance characteristics of two alfalfa genotypes. These findings, which may be valuable for drought resistance breeding, warrant further gene functional analysis to augment currently available information and to clarify the drought stress regulatory mechanisms of alfalfa and other plants.

Expression of cyanobacterial genes enhanced CO2 assimilation and biomass production in transgenic Arabidopsis thaliana

Background

Photosynthesis is a key process in plants that is compromised by the oxygenase activity of Rubisco, which leads to the production of toxic compound phosphoglycolate that is catabolized by photorespiratory pathway. Transformation of plants with photorespiratory bypasses have been shown to reduce photorespiration and enhance plant biomass. Interestingly, engineering of a single gene from such photorespiratory bypasses has also improved photosynthesis and plant productivity. Although single gene transformations may not completely reduce photorespiration, increases in plant biomass accumulation have still been observed indicating an alternative role in regulating different metabolic processes. Therefore, the current study was aimed at evaluating the underlying mechanism (s) associated with the effects of introducing a single cyanobacterial glycolate decarboxylation pathway gene on photosynthesis and plant performance.

Methods

Transgenic Arabidopsis thaliana plants (GD, HD, OX) expressing independently cyanobacterial decarboxylation pathway genes i.e., glycolate dehydrogenase, hydroxyacid dehydrogenase, and oxalate decarboxylase, respectively, were utilized. Photosynthetic, fluorescence related, and growth parameters were analyzed. Additionally, transcriptomic analysis of GD transgenic plants was also performed.

Results

The GD plants exhibited a significant increase (16%) in net photosynthesis rate while both HD and OX plants showed a non-significant (11%) increase as compared to wild type plants (WT). The stomatal conductance was significantly higher (24%) in GD and HD plants than the WT plants. The quantum efficiencies of photosystem II, carbon dioxide assimilation and the chlorophyll fluorescence-based photosynthetic electron transport rate were also higher than WT plants. The OX plants displayed significant reductions in the rate of photorespiration relative to gross photosynthesis and increase in the ratio of the photosynthetic electron flow attributable to carboxylation reactions over that attributable to oxygenation reactions. GD, HD and OX plants accumulated significantly higher biomass and seed weight. Soluble sugars were significantly increased in GD and HD plants, while the starch levels were higher in all transgenic plants. The transcriptomic analysis of GD plants revealed 650 up-regulated genes mainly related to photosynthesis, photorespiratory pathway, sucrose metabolism, chlorophyll biosynthesis and glutathione metabolism.

Conclusion

This study revealed the potential of introduced cyanobacterial pathway genes to enhance photosynthetic and growth-related parameters. The upregulation of genes related to different pathways provided evidence of the underlying mechanisms involved particularly in GD plants. However, transcriptomic profiling of HD and OX plants can further help to identify other potential mechanisms involved in improved plant productivity.

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.

Cut Flower Characteristics and Growth Traits under Salt Stress in Lily Cultivars

Salt stress is a major constraint of crop productivity because it reduces yield and limits the expansion of agriculture. This study investigated salt tolerance in 26 cultivars of cut lilies (Lilium hybrids) by examining the effect of salt stress on the growth and morphological characteristics of flowers and leaves and their physiological properties (chlorophyll a fluorescence). Salt stress significantly affected the growth and development of cut lilies. Canonical discriminant analysis indicates that the middle leaf width, number of flowers, first flower diameter, petal width, and chlorophyll a fluorescence were correlated with salt stress, whereas plant height, the middle leaf length, days to flowering, and sepal width were less affected by the stress. The cultivars examined were divided into three groups: Group 1 included the salt-sensitive cultivars, which failed to develop normal flowers; Group 2 included cultivars sensitive to salt stress but tolerant to osmotic stress; and Group 3 was the salt-tolerant group, which developed commercially valuable flowers. In conclusion, the cultivars contained a variable range of cut flower characteristics and growth traits that can be employed for lily breeding programs and as material for molecular mechanisms and signaling networks under salt stress.

Full-Length Transcriptome Sequencing and Comparative Transcriptome Analysis to Evaluate Drought and Salt Stress in Iris lactea var. chinensis

Iris lactea var. chinensis (I. lactea var. chinensis) is a perennial herb halophyte with salt and drought tolerance. In this study, full-length transcripts of I. lactea var. chinensis were sequenced using the PacBio RSII sequencing platform. Moreover, the transcriptome was investigated under NaCl or polyethylene glycol (PEG) stress. Approximately 30.89 G subreads were generated and 31,195 unigenes were obtained by clustering the same isoforms by the PacBio RSII platform. A total of 15,466 differentially expressed genes (DEGs) were obtained under the two stresses using the Illumina platform. Among them, 9266 and 8390 DEGs were obtained under high concentrations of NaCl and PEG, respectively. In total, 3897 DEGs with the same expression pattern under the two stresses were obtained. The transcriptome expression profiles of I. lactea var. chinensis under NaCl or PEG stress obtained in this study may provide a resource for the same and different response mechanisms against different types of abiotic stress. Furthermore, the stress-related genes found in this study can provide data for future molecular breeding.

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

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