Deep-sequencing transcriptome analysis of field-grown Medicago sativa L. crown buds acclimated to freezing stress

Identification of Cold Tolerance Transcriptional Regulatory Genes in Seedlings of Medicago sativa L. and Medicago falcata L

Alfalfa species Medicago sativa L. (MS) and Medicago falcata L. (MF), globally prominent perennial leguminous forages, hold substantial economic value. However, our comprehension of the molecular mechanisms governing their resistance to cold stress remains limited. To address this knowledge gap, we scrutinized and compared MS and MF cold-stress responses at the molecular level following 24 h and 120 h low-temperature exposure (4 °C). Our study revealed that MF had superior physiological resilience to cold stress compared with MS, and its morphology was healthier under cold stress, and its malondialdehyde content and superoxide dismutase activity increased, first, and then decreased, while the soluble sugar content continued to accumulate. Transcriptome analysis showed that after 120 h of exposure, there were different gene-expression patterns between MS and MF, including 1274 and 2983 genes that were continuously up-regulated, respectively, and a total of 923 genes were included, including star cold-resistant genes such as ICE1 and SIP1. Gene ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed numerous inter-species differences in sustained cold-stress responses. Notably, MS-exclusive genes included a single transcription factor (TF) gene and several genes associated with a single DNA repair-related pathway, whereas MF-exclusive genes comprised nine TF genes and genes associated with 14 pathways. Both species exhibited high-level expression of genes encoding TFs belonging to AP2-EREBP, ARR-B, and bHLH TF families, indicating their potential roles in sustaining cold resistance in alfalfa-related species. These findings provide insights into the molecular mechanisms governing cold-stress responses in MS and MF, which could inform breeding programs aimed at enhancing cold-stress resistance in alfalfa cultivars.

Exploring the Structure and Substance Metabolism of a Medicago sativa L. Stem Base

The stem base of alfalfa is a critical part for its overwintering, regeneration, and yield. To better understand the specificity and importance of the stem base, we analyzed the structure, metabolic substances, and transcriptome of the stem base using anatomical techniques, ultra-high performance liquid chromatography tandem mass spectrometry (UPLC-MS/MS), and RNA sequencing (RNA-seq), and compared it with stems and roots. The anatomical structure shows that the ratio of xylem to phloem changes at the base of the stem. A total of 801 compounds involved in 91 metabolic pathways were identified from the broadly targeted metabolome. Transcriptome analysis revealed 4974 differentially expressed genes (DEGs) at the stem base compared to the stem, and 5503 DEGs compared to the root. Comprehensive analyses of differentially accumulated compounds (DACs) and DEGs, in the stem base vs. stem, identified 10 valuable pathways, including plant hormone signal transduction, zeatin biosynthesis, α-Linolenic acid metabolism, histidine metabolism, carbon metabolism, carbon fixation in photosynthetic organisms, pentose phosphate pathway, galactose metabolism, and fructose and mannose metabolism. The pathways of plant hormone signal transduction and carbon metabolism were also identified by comparing the stem base with the roots. Taken together, the stem base of alfalfa is the transition region between the stem and root in morphology; in terms of material metabolism, its growth, development, and function are regulated through hormones and sugars.

Exogenous application of salicylic acid improves freezing stress tolerance in alfalfa

Freezing stress is one of the most detrimental environmental factors that can seriously impact the growth, development, and distribution of alfalfa (Medicago sativa L.). Exogenous salicylic acid (SA) has been revealed as a cost-effective method of improving defense against freezing stress due to its predominant role in biotic and abiotic stress resistance. However, how the molecular mechanisms of SA improve freezing stress resistance in alfalfa is still unclear. Therefore, in this study, we used leaf samples of alfalfa seedlings pretreatment with 200 μM and 0 μM SA, which were exposed to freezing stress (-10°C) for 0, 0.5, 1, and 2h and allowed to recover at normal temperature in a growth chamber for 2 days, after which we detect the changes in the phenotypical, physiological, hormone content, and performed a transcriptome analysis to explain SA influence alfalfa in freezing stress. The results demonstrated that exogenous SA could improve the accumulation of free SA in alfalfa leaves primarily through the phenylalanine ammonia-lyase pathway. Moreover, the results of transcriptome analysis revealed that the mitogen-activated protein kinase (MAPK) signaling pathway-plant play a critical role in SA alleviating freezing stress. In addition, the weighted gene co-expression network analysis (WGCNA) found that MPK3, MPK9, WRKY22 (downstream target gene of MPK3), and TGACG-binding factor 1 (TGA1) are candidate hub genes involved in freezing stress defense, all of which are involved in the SA signaling pathway. Therefore, we conclude that SA could possibly induce MPK3 to regulate WRKY22 to participate in freezing stress to induced gene expression related to SA signaling pathway (NPR1-dependent pathway and NPR1-independent pathway), including the genes of non-expresser of pathogenesis-related gene 1 (NPR1), TGA1, pathogenesis-related 1 (PR1), superoxide dismutase (SOD), peroxidase (POD), ascorbate peroxidase (APX), glutathione-S-transferase (GST), and heat shock protein (HSP). This enhanced the production of antioxidant enzymes such as SOD, POD, and APX, which increases the freezing stress tolerance of alfalfa plants.

Integrative analysis of transcriptome and metabolome revealed the mechanisms by which flavonoids and phytohormones regulated the adaptation of alfalfa roots to NaCl stress

INTRODUCTION

Salinity critically affects the growth and development of alfalfa (Medicago sativa), making it necessary to understand the molecular mechanism of alfalfa's adaptation to salt stress.

METHODS

In this study, alfalfa roots were subjected to salt stress and transcriptomics and metabolomics analyses were performed.

RESULTS

The results showed that flavonoid synthesis, hormone synthesis, and transduction pathways may be involved in the alfalfa salt stress adaptation reaction, and that they are related. Combined analysis of differential genes and differential metabolites found that dihydroquercetin and beta-ring hydroxylase (LUT5), ABA responsive element binding factor 2 (ABF2), protein phosphatase PP2C (PP2C) and abscisic acid (ABA) receptor PYL2 (PYL), luteolinidin was significantly correlated with PP2C and phytochrome-interacting factor 4 (PIF4) and (+)-7-isomethyl jasmonate were significantly correlated with flavonol synthase (FLS) gene. (+)-7-isomethyl jasmonate and homoeriodictyol chalcone were significantly correlated with peroxidase (POD). POD was significantly up-regulated under NaCl stress for 6 and 24 h. Moreover, flavonoids, gibberellin (GA), jasmonic acid (JA) and ABA were suggested to play an important role in alfalfa's response to salt stress. Further, GA,ABA, and JA may be involved in the regulation of flavonoids to improve alfalfa's salt tolerance, and JA may be a key signal to promote the synthesis of flavonoids.

DISCUSSION

This study revealed the possible molecular mechanism of alfalfa adaptation to salt stress, and identified a number of salt-tolerance candidate genes from the synthesis and signal transduction pathways of flavonoids and plant hormones, providing new insights into the regulatory network of alfalfa response to salt stress.

Genome-wide identification and expression analysis of DREB genes in alfalfa (Medicago sativa) in response to cold stress

Dehydration-responsive element-binding proteins (DREBs) belong to members of the AP2/ERF transcription factor superfamily, which has been reported to involve various abiotic-stress responses and tolerance in plants. However, research on the DREB-family is still limited in alfalfa (Medicago sativa L.), a forage legume cultivated worldwide. The recent genome-sequence release of the alfalfa cultivar "XinJiangDaYe" allowed us to identify 172 DREBs by a multi-step homolog search. The phylogenetic analysis indicated that such MsDREBs could be classified into 5 groups, namely A-1 (56 members), A-2 (39), A-3 (3), A-4 (61) and 13 (A-5 (13), thus adding substantial new members to the DREB-family in alfalfa. Furthermore, a comprehensive survey in silico of conserved motif, gene structure, molecular weight, and isoelectric point (pI) as well as gene expression was conducted. The resulting data showed that, for cold-stress response, 33 differentially expressed MsDREBs were identified with a threshold of Log2-fold > 1, and most of which were transcriptionally upregulated within 48 h during a cold treatment(s). Moreover, the expression profiling of MsDREBs from two ecotypes of alfalfa subspecies i.e. M. sativa ssp. falcata (F56, from a colder region of Central Asia) and M. sativa ssp. sativa (B47, from Near East) revealed that most of the cold-stress responsive MsDREBs exhibited a significantly lower expression in F56, leading to a proposal of the existence of a distinct mechanism(s) for cold tolerance regulated by DREB-related action, which would have been evolved in alfalfa with a genotypic specificity. Additionally, by examining the transcriptome of a freezing-tolerance species (M. sativa cv. Zhaodong), eight DREBs were found to be implicated in a long-term freezing-stress adaptation with a great potential. Taken together, the current genome-wide identification in alfalfa points to the importance of some MsDREBs in the cold-stress response, providing some promising molecular targets to be functionally characterized for the improvement of cold tolerance in crops including alfalfa.

Genome-Wide Identification of bHLH Transcription Factor in Medicago sativa in Response to Cold Stress

Alfalfa represents one of the most important legume forages, and it is also applied as an organic fertilizer to improve soil quality. However, this perennial plant is native to warmer temperate regions, and its valuable cold-acclimation-related regulatory mechanisms are still less known. In higher plants, the bHLH transcription factors play pleiotropic regulatory roles in response to abiotic stresses. The recently released whole genome sequencing data of alfalfa allowed us to identify 469 MsbHLHs by multi-step homolog search. Herein, we primarily identified 65 MsbHLH genes that significantly upregulated under cold stress, and such bHLHs were classified into six clades according to their expression patterns. Interestingly, the phylogenetic analysis and conserved motif screening of the cold-induced MsbHLHs showed that the expression pattern is relatively varied in each bHLH subfamily, this result indicating that the 65 MsbHLHs may be involved in a complex cold-responsive regulatory network. Hence, we analyzed the TFBSs at promoter regions that unraveled a relatively conserved TFBS distribution with genes exhibiting similar expression patterns. Eventually, to verify the core components involved in long-term cold acclimation, we examined transcriptome data from a freezing-tolerant species (cv. Zhaodong) in the field and compared the expression of cold-sensitive/tolerant subspecies of alfalfa, giving 11 bHLH as candidates, which could be important for further cold-tolerance enhancement and molecular breeding through genetic engineering in alfalfa.

Root rot a silent alfalfa killer in China: Distribution, fungal, and oomycete pathogens, impact of climatic factors and its management

Alfalfa plays a significant role in the pasture ecosystems of China's north, northeast, and northwest regions. It is an excellent forage for livestock, improves soil structure, prevents soil erosion, and has ecological benefits. Presently root rot is a significant threat to the alfalfa productivity because of the survival of the pathogens as soil-borne and because of lack of microbial competition in the impoverished nutrient-deficient soils and resistant cultivars. Furthermore, these regions' extreme ecological and environmental conditions predispose alfalfa to root rot. Moisture and temperature, in particular, have a considerable impact on the severity of root rot. Pathogens such as Fusarium spp. and Rhizoctonia solani are predominant, frequently isolated, and of major concern. These pathogens work together as disease complexes, so finding a host genotype resistant to disease complexes is challenging. Approaches to root rot control in these regions include mostly fungicides treatments and cultural practices and very few reports on the usage of biological control agents. As seed treatment, fungicides such as carbendazim are frequently used to combat root rot; however, resistance to fungicides has arisen. However, breeding and transgenic approaches could be more efficient and sustainable long-term control strategies, especially if resistance to disease complexes may be identified. Yet, research in China is mainly limited to field investigation of root rot and disease resistance evaluation. In this review, we describe climatic conditions of pastoral regions and the role of alfalfa therein and challenges of root rot, the distribution of root rot in the world and China, and the impact of root rot pathogens on alfalfa in particular R. solani and Fusarium spp., effects of environmental factors on root rot and summarize to date disease management approach.

Low Temperature Stress Tolerance: An Insight Into the Omics Approaches for Legume Crops

The change in climatic conditions is the major cause for decline in crop production worldwide. Decreasing crop productivity will further lead to increase in global hunger rate. Climate change results in environmental stress which has negative impact on plant-like deficiencies in growth, crop yield, permanent damage, or death if the plant remains in the stress conditions for prolonged period. Cold stress is one of the main abiotic stresses which have already affected the global crop production. Cold stress adversely affects the plants leading to necrosis, chlorosis, and growth retardation. Various physiological, biochemical, and molecular responses under cold stress have revealed that the cold resistance is more complex than perceived which involves multiple pathways. Like other crops, legumes are also affected by cold stress and therefore, an effective technique to mitigate cold-mediated damage is critical for long-term legume production. Earlier, crop improvement for any stress was challenging for scientific community as conventional breeding approaches like inter-specific or inter-generic hybridization had limited success in crop improvement. The availability of genome sequence, transcriptome, and proteome data provides in-depth sight into different complex mechanisms under cold stress. Identification of QTLs, genes, and proteins responsible for cold stress tolerance will help in improving or developing stress-tolerant legume crop. Cold stress can alter gene expression which further leads to increases in stress protecting metabolites to cope up the plant against the temperature fluctuations. Moreover, genetic engineering can help in development of new cold stress-tolerant varieties of legume crop. This paper provides a general insight into the "omics" approaches for cold stress in legume crops.

Identification and Characterization of Regulatory Pathways Controlling Dormancy Under Lower Temperature in Alfalfa (Medicago sativa L.)

Alfalfa (Medicago sativa L.), a kind of high-quality perennial legume forage, is widely distributed in the northern regions of China. In recent years, low temperatures have frequently occurred and limited alfalfa productivity and survival in early spring and late fall. However, the underlying molecular mechanisms of alfalfa response to cold tolerance are not well-documented. In this study, dormancy and non-dormancy alfalfa standard varieties were characterized under low-temperature stress. Our analysis revealed that plant height of the dormancy genotype was strongly inhibited by low temperature; flavonoids content, and higher expression of flavonoids biosynthesis genes (chalcone synthase, leucoanthocyanidin dioxygenase, and flavonoid 3'-monooxygenase) may play essential roles in response to low-temperature stress in dormancy genotype alfalfa. Further analyses revealed that receptor-like kinase family genes (such as cysteine-rich RLK10, lectin protein kinase, and S-locus glycoprotein like kinase), RNA and protein synthesis genes (RNA polymerases, ribosomal protein, and protein phosphatase 2C family protein), and proteasome degradation pathway genes (such as F-box family protein, RING/U-box superfamily protein, and zinc finger family protein) also highly upregulated and contributed to cold tolerance phenotype in dormancy genotype alfalfa. This will provide new insights into future studies for cold tolerance in alfalfa and offer new target genes for further functional characterization and genetic improvement of alfalfa.

Comparative Transcriptome Analysis Reveals New Insight of Alfalfa (Medicago sativa L.) Cultivars in Response to Abrupt Freezing Stress

Freezing stress is a major limiting environmental factor that affects the productivity and distribution of alfalfa (Medicago sativa L.). There is growing evidence that enhancing freezing tolerance through resistance-related genes is one of the most efficient methods for solving this problem, whereas little is known about the complex regulatory mechanism of freezing stress. Herein, we performed transcriptome profiling of the leaves from two genotypes of alfalfa, freezing tolerance "Gannong NO.3" and freezing-sensitive "WL326GZ" exposure to -10°C to investigate which resistance-related genes could improve the freezing tolerance. Our results showed that a total of 121,366 genes were identified, and there were 7,245 differentially expressed genes (DEGs) between the control and treated leaves. In particular, the DEGs in "Gannong NO.3" were mainly enriched in the metabolic pathways and biosynthesis of secondary metabolites, and most of the DEGs in "WL326GZ" were enriched in the metabolic pathways, the biosynthesis of secondary metabolites, and plant-pathogen interactions. Moreover, the weighted gene co-expression network analysis (WGCNA) showed that ATP-binding cassette (ABC) C subfamily genes were strongly impacted by freezing stress, indicating that ABCC8 and ABCC3 are critical to develop the freezing tolerance. Moreover, our data revealed that numerous Ca2+ signal transduction and CBF/DREB1 pathway-related genes were severely impacted by the freezing resistance, which is believed to alleviate the damage caused by freezing stress. Altogether, these findings contribute the comprehensive information to understand the molecular mechanism of alfalfa adaptation to freezing stress and further provide functional candidate genes that can adapt to abiotic stress.

Identification of RNA helicases in Medicago truncatula and their expression patterns under abiotic stress

RNA helicase catalyzes the denaturation of DNA or the unwinding of double-stranded RNA. It is vital to RNA splicing, transport, editing, degradation and the initiation of protein translation. However, the function of RNA helicase in Medicago truncatula has rarely been reported. In this study, 170 putative RNA helicase genes were identified in the M. truncatula genome, and classified into three subfamilies based on the presence of either a DEAD-box (52 genes), DEAH-box (38 genes), or DExD/H-box (80 genes) in their coding regions. Additionally, conserved helicase_C domains and other functional domains (e.g., the HA2, DUF, and ZnF domains) were also present in these genes. Chromosomal mapping and synteny analyses showed that there were tandem and segment duplications of RNA helicase genes. Furthermore, transcriptome and real-time PCR analysis showed that the expression of 35 RNA helicase genes was affected by abiotic stress. To be specific, 17, 12 and 19 genes were regulated by salt, drought and cold stress, respectively. It is worth noting that MtDEAD8, MtDEAH3, MtDExD/H18 and MtDExD/H23 responded to all three types of stress. These results provide valuable information for understanding the RNA helicase genes in M. truncatula and their abiotic stress-related functions.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01087-y.

Alfalfa (Medicago sativa L.) MsCML46 gene encoding calmodulin-like protein confers tolerance to abiotic stress in tobacco

KEY MESSAGE

MsCML46 enhances tolerance to abiotic stresses through alleviating osmotic stress and oxidative damage by regulating the expression of stress-related genes to optimize osmolytes levels and antioxidant enzyme activity in transgenic tobacco. Abiotic stresses are major environmental factors that constraint crop productivity worldwide. Various stimuli regulate intracellular calcium levels and calcium-mediated signal transduction, and cellular responses. Ca2+ signals are perceived by different Ca2+ receptors. Calmodulin-like protein (CML) is one of the best-characterized Ca2+ sensors which shares sequence similarity with highly conserved calmodulin (CaM) ubiquitously expressed in plants. Currently, the molecular and physiological functions of CMLs are largely unknown. In this study, the MsCML46 was characterized in alfalfa (Medicago sativa cv. Zhaodong) under freezing stress. Results showed that MsCML46 was localized to the cytoplasm of Arabidopsis, and its expression was strongly elevated by cold, drought, salt, saline-alkali, and ABA treatments. Overexpressing MsCML46 in tobacco enhanced tolerance to freezing, drought, and salt stresses as evidenced by improved contents of osmotic regulatory solutes and antioxidant enzyme activity but decreased reactive oxygen species (ROS) accumulation. Furthermore, cold, drought, and salt stresses increased the expression of stress-related genes in transgenic tobacco. MsCML46 binds free Ca2+ to promote signal transduction and maintain higher K+/Na+ ratio. In this way, it protects intracellular homeostasis under sodium ion toxicity. These results suggest that MsCML46 plays a crucial role in resisting abiotic stresses and can be exploited in genetic engineering for crops.

Qualitative and quantitative analysis of phenolic compounds in blueberries and protective effects on hydrogen peroxide-induced cell injury

This work was conducted to optimize an accelerated solvent extraction for ultra high performance liquid chromatography-mass spectrometry/mass spectrometry analysis of blueberry phenolic compounds. The conditions for accelerated solvent extraction were verified using response surface methodology to obtain the following optimized conditions: ethanol concentration (pH = 3), 48%; temperature, 50℃, and static cycle times, 3. Further, ultra high performance liquid chromatography with quadrupole Exactive Orbitrap mass spectrometry and ultra high performance liquid chromatography with triple-quadrupole tandem mass methods for determination of the detailed phenolic composition were developed and validated. Total of 81 phenolic compounds were identified by ultra high performance liquid chromatography with quadrupole Exactive Orbitrap mass spectrometry including 23 anthocyanins, 32 flavonols, 11 proanthocyanidins, 2 other flavonoids, and 13 phenolic acids. Fifty-one of these compounds have been simultaneously quantified by ultra high performance liquid chromatography with triple-quadrupole tandem mass including 31 anthocyanins, 8 flavonols, 6 proanthocyanidins, 2 other flavonoids, and 8 phenolic acids. Malvidin-dinhexoside has, for the first time, been detected in wild. Moreover, by verifying the protection on PC12 cells against oxidative damage, it was showed that the phenolic extracts (500 µg/mL) can improve significantly the viability (9.26-24.78%) of hydrogen peroxide-induced PC12 cells, activities of superoxide dismutase (34.59-37.90 U/mg) and glutathione peroxidase (6.87-14.42 mU/mg) and decrease the content of malonic dialdehyde (13.27-24.62 nmol/mg). Correlation analysis suggested that anthocyanins might contribute most to these activities.

A Review of Unreduced Gametes and Neopolyploids in Alfalfa: How to Fill the Gap between Well-Established Meiotic Mutants and Next-Generation Genomic Resources

The gene flow mediated by unreduced gametes between diploid and tetraploid plants of the Medicagosativa-coerulea-falcata complex is pivotal for alfalfa breeding. Sexually tetraploidized hybrids could represent the best way to exploit progressive heterosis simultaneously derived from gene diversity, heterozygosity, and polyploidy. Moreover, unreduced gametes combined with parthenogenesis (i.e., apomixis) would enable the cloning of plants through seeds, providing a unique opportunity for the selection of superior genotypes with permanently fixed heterosis. This reproductive strategy has never been detected in the genus Medicago, but features of apomixis, such as restitutional apomeiosis and haploid parthenogenesis, have been reported. By means of an original case study, we demonstrated that sexually tetraploidized plants maintain apomeiosis, but this trait is developmentally independent from parthenogenesis. Alfalfa meiotic mutants producing unreduced egg cells revealed a null or very low capacity for parthenogenesis. The overall achievements reached so far are reviewed and discussed along with the efforts and strategies made for exploiting reproductive mutants that express apomictic elements in alfalfa breeding programs. Although several studies have investigated the cytological mechanisms responsible for 2n gamete formation and the inheritance of this trait, only a very small number of molecular markers and candidate genes putatively linked to unreduced gamete formation have been identified. Furthermore, this scenario has remained almost unchanged over the last two decades. Here, we propose a reverse genetics approach, by exploiting the genomic and transcriptomic resources available in alfalfa. Through a comparison with 9 proteins belonging to Arabidopsis thaliana known for their involvement in 2n gamete production, we identified 47 orthologous genes and evaluated their expression in several tissues, paving the way for novel candidate gene characterization studies. An overall view on strategies suitable to fill the gap between well-established meiotic mutants and next-generation genomic resources is presented and discussed.

Identification of the trehalose-6-phosphate synthase gene family in Medicago truncatula and expression analysis under abiotic stresses

Trehalose-6-phosphate synthase (TPS) exerts important functions related to plant desiccation tolerance and responses to environmental stimuli. However, in Medicago truncatula, the TPS family has not been reported to date. This study found 11 MtTPS genes in the genome of M. truncatula, which could be divided into two subfamilies: Class I and Class II. All TPS family members have a TPS domain (Glyco transf_20) at the N-terminus and a TPP domain (Trehalose_PPase) at the C-terminus. Interestingly, the genetic structures differ between Class I and Class II, Class I members have more introns than Class II members. Furthermore, transcriptome and real-time PCR analysis showed that five MtTPS genes could be induced by drought, salt or cold. Specifically, MtTPS2, MtTPS8, MtTPS9, MtTPS11 were up-regulated under both drought and salt treatment, particularly, MtTPS8 and MtTPS9 can also be induced by cold, while MtTPS7 only responded to salt stress. In summary, this study provides the foundation for further research on TPS genes in M. truncatula and their regulatory function in response to abiotic stresses.

Genome-Wide Association Analysis Coupled With Transcriptome Analysis Reveals Candidate Genes Related to Salt Stress in Alfalfa (Medicago sativa L.)

Salt stress is the main abiotic factor affecting alfalfa yield and quality. However, knowledge of the genetic basis of the salt stress response in alfalfa is still limited. Here, a genome-wide association study (GWAS) involving 875,023 single-nucleotide polymorphisms (SNPs) was conducted on 220 alfalfa varieties under both normal and salt-stress conditions. Phenotypic analysis showed that breeding status and geographical origin play important roles in the alfalfa salt stress response. For germination ability under salt stress, a total of 15 significant SNPs explaining 9%-14% of the phenotypic variation were identified. For tolerance to salt stress in the seedling stage, a total of 18 significant SNPs explaining 12%-23% of the phenotypic variation were identified. Transcriptome analysis revealed 2,097 and 812 differentially expressed genes (DEGs) that were upregulated and 2,445 and 928 DEGs that were downregulated in the leaves and roots, respectively, under salt stress. Among these DEGs, many encoding transcription factors (TFs) were found, including MYB-, CBF-, NAC-, and bZIP-encoding genes. Combining the results of our GWAS analysis and transcriptome analysis, we identified a total of eight candidate genes (five candidate genes for tolerance to salt stress and three candidate genes for germination ability under salt stress). Two SNPs located within the upstream region of MsAUX28, which encodes an auxin response protein, were significantly associated with tolerance to salt stress. The two significant SNPs within the upstream region of MsAUX28 existed as three different haplotypes in this panel. Hap 1 (G/G, A/A) was under selection in the alfalfa domestication and improvement process.

High-Throughput Sequencing for Deciphering the Virome of Alfalfa (Medicago sativa L.)

Alfalfa (Medicago sativa L.), also known as lucerne, is a major forage crop worldwide. In the United States, it has recently become the third most valuable field crop, with an estimated value of over $9.3 billion. Alfalfa is naturally infected by many different pathogens, including viruses, obligate parasites that reproduce only inside living host cells. Traditionally, viral infections of alfalfa have been considered by breeders, growers, producers and researchers to be diseases of limited importance, although they are widespread in all major cultivation areas. However, over the past few years, due to the rapid development of high-throughput sequencing (HTS), viral metagenomics, bioinformatics tools for interpreting massive amounts of HTS data and the increasing accessibility of public data repositories for transcriptomic discoveries, several emerging viruses of alfalfa with the potential to cause serious yield losses have been described. They include alfalfa leaf curl virus (family Geminiviridae), alfalfa dwarf virus (family Rhabdoviridae), alfalfa enamovirus 1 (family Luteoviridae), alfalfa virus S (family Alphaflexiviridae) and others. These discoveries have called into question the assumed low economic impact of viral diseases in alfalfa and further suggested their possible contribution to the severity of complex infections involving multiple pathogens. In this review, we will focus on viruses of alfalfa recently described in different laboratories on the basis of the above research methodologies.

Biotechnological Perspectives of Omics and Genetic Engineering Methods in Alfalfa

For several decades, researchers are working to develop improved major crops with better adaptability and tolerance to environmental stresses. Forage legumes have been widely spread in the world due to their great ecological and economic values. Abiotic and biotic stresses are main factors limiting legume production, however, alfalfa (Medicago sativa L.) shows relatively high level of tolerance to drought and salt stress. Efforts focused on alfalfa improvements have led to the release of cultivars with new traits of agronomic importance such as high yield, better stress tolerance or forage quality. Alfalfa has very high nutritional value due to its efficient symbiotic association with nitrogen-fixing bacteria, while deep root system can help to prevent soil water loss in dry lands. The use of modern biotechnology tools is challenging in alfalfa since full genome, unlike to its close relative barrel medic (Medicago truncatula Gaertn.), was not released yet. Identification, isolation, and improvement of genes involved in abiotic or biotic stress response significantly contributed to the progress of our understanding how crop plants cope with these environmental challenges. In this review, we provide an overview of the progress that has been made in high-throughput sequencing, characterization of genes for abiotic or biotic stress tolerance, gene editing, as well as proteomic and metabolomics techniques bearing biotechnological potential for alfalfa improvement.

Identification of emerging viral genomes in transcriptomic datasets of alfalfa (Medicago sativa L.)

BACKGROUND

Publicly available transcriptomic datasets have become a valuable tool for the discovery of new pathogens, particularly viruses. In this study, several coding-complete viral genomes previously not found or experimentally confirmed in alfalfa were identified in the plant datasets retrieved from the NCBI Sequence Read Archive.

METHODS

Publicly available Medicago spp. transcriptomic datasets were retrieved from the NCBI SRA database. The raw reads were first mapped to the reference genomes of Medicago sativa and Medigago truncatula followed by the alignment of the unmapped reads to the NCBI viral genome database and de novo assembly using the SPAdes tool. When possible, assemblies were experimentally confirmed using 5'/3' RACE and RT-PCRs.

RESULTS

Twenty three different viruses were identified in the analyzed datasets, of which several represented emerging viruses not reported in alfalfa prior to this study. Among them were two strains of cnidium vein yellowing virus, lychnis mottle virus and Cactus virus X, for which coding-complete genomic sequences were obtained by a de novo assembly.

CONCLUSIONS

The results improve our knowledge of the diversity and host range of viruses infecting alfalfa, provide essential tools for their diagnostics and characterization and demonstrate the utility of transcriptomic datasets for the discovery of new pathogens.

Identification of Genes Differentially Expressed in Response to Cold in Pisum sativum Using RNA Sequencing Analyses

Low temperature stress affects growth and development in pea (Pisum sativum L.) and decreases yield. In this study, RNA sequencing time series analyses performed on lines, Champagne frost-tolerant and Térèse frost-sensitive, during a low temperature treatment versus a control condition, led us to identify 4981 differentially expressed genes. Thanks to our experimental design and statistical analyses, we were able to classify these genes into three sets. The first one was composed of 2487 genes that could be related to the constitutive differences between the two lines and were not regulated during cold treatment. The second gathered 1403 genes that could be related to the chilling response. The third set contained 1091 genes, including genes that could be related to freezing tolerance. The identification of differentially expressed genes related to cold, oxidative stress, and dehydration responses, including some transcription factors and kinases, confirmed the soundness of our analyses. In addition, we identified about one hundred genes, whose expression has not yet been linked to cold stress. Overall, our findings showed that both lines have different characteristics for their cold response (chilling response and/or freezing tolerance), as more than 90% of differentially expressed genes were specific to each of them.

Analysis of transcripts and splice isoforms in Medicago sativa L. by single-molecule long-read sequencing

The full-length transcriptome of alfalfa was analyzed with PacBio single-molecule long-read sequencing technology. The transcriptome data provided full-length sequences and gene isoforms of transcripts in alfalfa, which will improve genome annotation and enhance our understanding of the gene structure of alfalfa. As an important forage, alfalfa (Medicago sativa L.) is world-wide planted. For its complexity of genome and unfinished whole genome sequencing, the sequences and complete structure of mRNA transcripts remain unclear in alfalfa. In this study, single-molecule long-read sequencing was applied to investigate the alfalfa transcriptome using the Pacific Biosciences platform, and a total of 113,321 transcripts were obtained from young, mature and senescent leaves. We identified 72,606 open reading frames including 46,616 full-length ORFs, 1670 transcription factors from 54 TF families and 44,040 simple sequence repeats from 30,797 sequences. A total of 7568 alternative splicing events was identified and the majority of alternative splicing events in alfalfa was intron retention. In addition, we identified 17,740 long non-coding RNAs. Our results show the feasibility of deep sequencing full-length RNA from alfalfa transcriptome on a single-molecule level.

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

BACKGROUND

Alfalfa is the most extensively cultivated forage legume. Salinity is a major environmental factor that impacts on alfalfa's productivity. However, little is known about the molecular mechanisms underlying alfalfa responses to salinity, especially the relative contribution of the two important components of osmotic and ionic stress.

RESULTS

In this study, we constructed the first full-length transcriptome database for alfalfa root tips under continuous NaCl and mannitol treatments for 1, 3, 6, 12, and 24 h (three biological replicates for each time points, including the control group) via PacBio Iso-Seq. This resulted in the identification of 52,787 full-length transcripts, with an average length of 2551 bp. Global transcriptional changes in the same 33 stressed samples were then analyzed via BGISEQ-500 RNA-Seq. Totals of 8861 NaCl-regulated and 8016 mannitol-regulated differentially expressed genes (DEGs) were identified. Metabolic analyses revealed that these DEGs overlapped or diverged in the cascades of molecular networks involved in signal perception, signal transduction, transcriptional regulation, and antioxidative defense. Notably, several well characterized signalling pathways, such as CDPK, MAPK, CIPK, and PYL-PP2C-SnRK2, were shown to be involved in osmotic stress, while the SOS core pathway was activated by ionic stress. Moreover, the physiological shifts of catalase and peroxidase activity, glutathione and proline content were in accordance with dynamic transcript profiles of the relevant genes, indicating that antioxidative defense system plays critical roles in response to salinity stress.

CONCLUSIONS

Overall, our study provides evidence that the response to salinity stress in alfalfa includes both osmotic and ionic components. The key osmotic and ionic stress-related genes are candidates for future studies as potential targets to improve resistance to salinity stress via genetic engineering.

Multiple Regulatory Networks Are Activated during Cold Stress in Medicago sativa L

Cultivated alfalfa (Medicago sativa L.) is one of the most important perennial legume forages in the world, and it has considerable potential as a valuable forage crop for livestock. However, the molecular mechanisms underlying alfalfa responses to cold stress are largely unknown. In this study, the transcriptome changes in alfalfa under cold stress at 4 °C for 2, 6, 24, and 48 h (three replicates for each time point) were analyzed using the high-throughput sequencing platform, BGISEQ-500, resulting in the identification of 50,809 annotated unigenes and 5283 differentially expressed genes (DEGs). Metabolic pathway enrichment analysis demonstrated that the DEGs were involved in carbohydrate metabolism, photosynthesis, plant hormone signal transduction, and the biosynthesis of amino acids. Moreover, the physiological changes of glutathione and proline content, catalase, and peroxidase activity were in accordance with dynamic transcript profiles of the relevant genes. Additionally, some transcription factors might play important roles in the alfalfa response to cold stress, as determined by the expression pattern of the related genes during 48 h of cold stress treatment. These findings provide valuable information for identifying and characterizing important components in the cold signaling network in alfalfa and enhancing the understanding of the molecular mechanisms underlying alfalfa responses to cold stress.

Breeding approaches and genomics technologies to increase crop yield under low-temperature stress

Improved knowledge about plant cold stress tolerance offered by modern omics technologies will greatly inform future crop improvement strategies that aim to breed cultivars yielding substantially high under low-temperature conditions. Alarmingly rising temperature extremities present a substantial impediment to the projected target of 70% more food production by 2050. Low-temperature (LT) stress severely constrains crop production worldwide, thereby demanding an urgent yet sustainable solution. Considerable research progress has been achieved on this front. Here, we review the crucial cellular and metabolic alterations in plants that follow LT stress along with the signal transduction and the regulatory network describing the plant cold tolerance. The significance of plant genetic resources to expand the genetic base of breeding programmes with regard to cold tolerance is highlighted. Also, the genetic architecture of cold tolerance trait as elucidated by conventional QTL mapping and genome-wide association mapping is described. Further, global expression profiling techniques including RNA-Seq along with diverse omics platforms are briefly discussed to better understand the underlying mechanism and prioritize the candidate gene (s) for downstream applications. These latest additions to breeders' toolbox hold immense potential to support plant breeding schemes that seek development of LT-tolerant cultivars. High-yielding cultivars endowed with greater cold tolerance are urgently required to sustain the crop yield under conditions severely challenged by low-temperature.

Exploring drought stress-regulated genes in senna (Cassia angustifolia Vahl.): a transcriptomic approach

De novo assembly of reads produced by next-generation sequencing (NGS) technologies offers a rapid approach to obtain expressed gene sequences for non-model organisms. Senna (Cassia angustifolia Vahl.) is a drought-tolerant annual undershrub of Caesalpiniaceae, a subfamily of Fabaceae. There are insufficient transcriptomic and genomic data in public databases for understanding the molecular mechanism underlying the drought tolerance of senna. Therefore, the main purpose of this study was to know the transcriptome profile of senna, with special reference to drought stress. RNA from two different stages of leaf development was extracted and sequenced separately using the Illumina technology. A total of 200 million reads were generated, and a de novo assembly of processed reads in the pooled transcriptome using Trinity yielded 43,413 transcripts which were further annotated using NCBI BLAST with "green plant database (txid 33090)," Swiss Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and Gene Ontology (GO). Out of the total transcripts, 42,280 (95.0 %) were annotated by BLASTX against the green plant database of NCBI. Senna transcriptome showed the highest similarity to Glycine max (41 %), followed by Phaseolus vulgaris (16 %), Cicer arietinum (15 %), and Medicago trancatula (5 %). The highest number of GO terms were enriched for the molecular functions category; of these "catalytic activity" (GO: 0003824) (25.10 %) and "binding activity" (GO: 0005488) (20.10 %) were most abundantly represented. We used InterProscan to see protein similarity at domain level; a total of 33,256 transcripts were annotated against the Pfam domains. The transcripts were assigned with various KEGG pathways. Coding DNA sequences (CDS) encoding various drought stress-regulated pathways such as signaling factors, protein-modifying/degrading enzymes, biosynthesis of phytohormone, phytohormone signaling, osmotically active compounds, free radical scavengers, chlorophyll metabolism, leaf cuticular wax, polyamines, and protective proteins were identified through BLASTX search. The lucine-rich repeat kinase family was the most abundantly found group of protein kinases. Orphan, bHLH, and bZIP family TFs were the most abundantly found in senna. Six genes encoding MYC2 transcription factor, 9-cis-epoxycarotenoid dioxygenase (NCED), l -ascorbate peroxidase (APX), aminocyclopropane carboxylate oxidase (ACO), abscisic acid 8'-hydroxylase (ABA), and WRKY transcription factor were confirmed through reverse transcriptase-PCR (RT-PCR) and Sanger sequencing for the first time in senna. The potential drought stress-related transcripts identified in this study provide a good start for further investigation into the drought adaptation in senna. Additionally, our transcriptome sequences are the valuable resource for accelerated genomics-assisted genetic improvement programs and facilitate manipulation of biochemical pathways for developing drought-tolerant genotypes of crop plants.