- 1Grassland Research Institute, Chinese Academy of Agricultural Science, Hohhot, Inner Mongolia, China
- 2College of Grassland, Qingdao Agricultural University, Qingdao, Shandong, China
- 3College of Life Science, Inner Mongolia University, Hohhot, Inner Mongolia, China
Introduction: Alfalfa (Medicago sativa L.) is a highly nutritious leguminous forage that plays an essential role in animal husbandry. In the middle and high latitudes of the northern hemisphere, there are problems with its low rates of overwintering and production. The application of phosphate (P) is an important measure to improve the cold resistance and production of alfalfa, but little is known about the mechanism of P in improving the cold resistance of alfalfa.
Methods: This study integrated the transcriptome and metabolome to explain the mechanism of alfalfa in response to low-temperature stress under two applications of P (50 and 200 mg kg-1) and a control of none applied.
Results: The application of P fertilizer improved the root structure and increased the content of soluble sugar and soluble protein in the root crown. In addition, there were 49 differentially expressed genes (DEGs) with 23 upregulated and 24 metabolites with 12 upregulated when 50 mg kg-1 of P was applied. In contrast, there were 224 DEGs with 173 upregulated and 12 metabolites with 6 upregulated in the plants treated with 200 mg kg-1 of P compared with the Control Check (CK). These genes and metabolites were significantly enriched in the biosynthesis of other secondary metabolites and the metabolic pathways of carbohydrates and amino acids. The integration of the transcriptome and metabolome indicated that P affected the biosynthesis of N-acetyl-L-phenylalanine, L-serine, lactose, and isocitrate during the period of increasing cold. It could also affect the expression of related genes that regulate cold tolerance in alfalfa.
Discussion: Our findings could contribute to a deeper understanding of the mechanism that alfalfa uses to tolerate cold and lay a theoretical foundation for breeding alfalfa that is highly efficient at utilizing phosphorus.
1 Introduction
As a cool-season perennial forage crop, alfalfa (Medicago sativa L.) plays a variety of beneficial roles for livestock, soil coverage, soil fertility, and human health (Gaafar et al., 2019). It has a high biomass yield, excellent nutritional quality, and wide adaption, and it can fix nitrogen (N) (Summers and Putnam, 2008). Alfalfa is grown on approximately 30 million ha of land worldwide (Acharya et al., 2020; Diatta et al., 2021), and the global production of alfalfa hay was 210.9 million metric tons (Edde, 2022). North America is the primary producer of alfalfa in the world; the United States has the largest planting area with more than 7.23 million ha of alfalfa grown annually, with Canada and Mexico close behind (Putnam et al., 2001; Wechsler et al., 2017; Edde, 2022). The cultivated area of alfalfa in China is approximately 3.77 × 106 hm2, which is approximately 13% of the world’s total arable land used to grow alfalfa. The growth of alfalfa is primarily distributed in northeastern, northwestern, and northern China (Yang et al., 2021). In these northern regions, the inability of alfalfa to withstand low subfreezing temperatures remains one of its biggest challenges (Castonguay et al., 2006; Rocher et al., 2015), which hampers the development of the alfalfa industry and animal husbandry.
The morphological and physiological characteristics of the roots are strongly associated with the cold resistance of alfalfa. The root types of alfalfa primarily include four types: tap, branch, rhizomatous, and creeping roots (Goplen, 1987; Hong et al., 1987). Tap-rooted types have a main root and a narrow, protruding crown; the branch-rooted types of alfalfa have a moderately wide crown and a number of primary roots; the rhizomatous- and creeping-rooted types have a protected crown, develop new roots readily, and can also develop stems from the roots. The shoots can separate from the maternal parent plant and become independent plants to survive; thus, they more easily regrow after freezing damage and show strong cold resistance (Goplen, 1987; Dou, 2011). A substantial amount of research has documented that the root crown, lateral roots, and root biomass affect the winter survival and persistence of alfalfa (Castonguay et al., 2006; Liu et al., 2015b; Xu et al., 2022). Previous studies have primarily focused on the root crown. It is a transitional plant structure that is located between the shoots and the root system. It is also the uppermost dormant organ of the winter plant body, which is crucial to overwintering and regeneration in the spring (Marquez-Ortiz et al., 1999; Liu et al., 2015b). The overwintering rate of alfalfa is closely related to the size and depth of the root crown; they grow more deeply and thicker in the soil, which helps the plants to successfully survive the winter (Wang, 2021b). This is considered to be a cold-sensitivity escape mechanism to prevent the exposure of the overwintering organ to low temperature. The major alfalfa varieties with strong cold resistance have more lateral roots because they help alfalfa plants to absorb more water and nutrient elements from the soil to supply the demand of overwintering plants for nutrients (Johnson et al., 1996; Wang et al., 2023). The root biomass is related to the accumulation of organic matter, and the herringbone branching is conducive to the absorption of deeper water by alfalfa, which contributes to improve the plant’s cold resistance (Viands, 1988). Its physiology, including the contents of soluble sugar and protein, is also closely related to cold resistance. Soluble sugar acts as an osmotic regulator, a cryoprotectant, and a signaling molecule to stabilize the cell membrane and scavenge reactive oxygen species (ROS) under low-temperature stress, and its contents increase during the freezing process to protect the plant cells (Trischuk et al., 2014). Therefore, the accumulation of sugar at the root crown before winter is related to the cold resistance of alfalfa (Cunningham et al., 1998; Bertrand et al., 2017). Soluble proteins are strongly hydrophilic and can enhance the water-holding capacity of cells. The accumulation of soluble proteins can bind more water to the cells and reduce the damage caused by low temperatures (Kontunen-Soppela et al., 2002; Moieni-Korbekandi et al., 2013). The increase in the soluble protein content of alfalfa in autumn and winter helps to enhance its resistance to low temperatures.
As an important means of regulating cultivation, nutrient management not only can improve the production of alfalfa, but it is also an essential measure to improve its stress resistance (Wang, 2021a). Phosphorus is one of the three essential nutrients for plant growth and development, which is absorbed in the form of phosphate (P) (Bhosale et al., 2018). It is an essential macroelement that plays a role in an array of processes, including energy generation, nucleic acid synthesis, photosynthesis, glycolysis, respiration, membrane synthesis and stability, respiration, carbohydrate metabolism, and N fixation (Vance et al., 2003; Raghothama and Karthikeyan, 2005). In addition, it plays a vital role in enhancing the adaptability of plants to the external environment (Yan et al., 2022). P is also an important component of phospholipids and ATP, which affect the resistance of plants to low temperatures (Li et al., 2016). Phospholipids in the plant cell membranes can interact with proteins, sugars, and other substances to alleviate the cell dehydration caused by low-temperature stress and protect the stability of cell membranes (Cheng, 2008). Plants can be negatively affected by the low levels of phosphorus that affect a variety of growing environments, particularly in soils that are calcareous or alkaline (Yue et al., 2019). A deficiency in P will affect plant photosynthesis, such as by limiting the distribution and utilization of carbohydrates and the absorption and transport of P (Terry and Tao, 1991). This can also cause an imbalance in the production and clearance of ROS and an array of physiological, biochemical, and metabolic changes resulting in damage to the plants (Lin et al., 2010). The application of P can increase the contents of soluble sugar in plants, and sucrose transport in the phloem requires ATP hydrolysis to provide energy so that P can regulate the metabolism and transport of sucrose in plant leaves (Tian et al., 2017). Current studies have shown that the application of moderate P fertilizer can promote the root growth of alfalfa and increase the content of cold-resistant substances, such as starch and soluble proteins, in the root crown of alfalfa to some extent (Shen et al., 2017). However, the related molecular mechanism remains unclear.
The recently developed technologies of high-throughput sequencing, high-resolution mass spectrometry, and information processing have made systems biology (omics) research a major focus of scientific investigation (Xin et al., 2019). The integrated analyses of transcriptomic and metabolic data obtained from two biological levels, i.e., transcript and metabolite levels, respectively, can serve as a useful way to study complex biological phenomena (Agarrwal et al., 2016). In addition, these integrated analyses have been applied to various plant biological processes, including the evolutionary adaption of poplar species (Populus L.) to salinity stress (Janz et al., 2010), rice (Oryza sativa L.) insect interaction research (Agarrwal et al., 2016), and the response of wild soybean (Glycine max L.) to N starvation (Liu et al., 2020). However, to our knowledge, there are few studies that have examined the mechanisms of the effect of P fertilizer on the cold resistance of alfalfa by integrating the analyses of transcriptomic and metabolic data.
In this study, alfalfa was exposed to three P treatments for 90 days. We measured the contents of the soluble sugar and protein of root crowns and some root indices. In addition, we performed transcriptomic and metabolomic studies among the three P treatments under cold temperature. Our integrated transcriptome–metabolome analysis allowed us to gain a deeper understanding of alfalfa’s response to the application of P in terms of cold tolerance. We sought to identify the strategies used by alfalfa to respond to the application of P during the increase in cold stress, which could be used to improve the cold tolerance and yield of alfalfa and breed new alfalfa varieties.
2 Materials and methods
2.1 Field design and sampling
The variety of alfalfa used in this study was Zhongcao No. 3, which was provided by the Grassland Research Institute of the Chinese Academy of Agricultural Sciences (GRI of CAAS) (Beijing, China). On September 5, 2020, an alfalfa plant that had grown for 5 years in the Agro-pastoral Experiment Station (40°34′E; 111°45′N; 1,050 m a.s.l.) of the GRI of CAAS was transferred to the greenhouse, and the seedlings were raised by stem cuttings in seedling trays on May 10, 2021. On July 15, 2021, cloned plants were selected and transplanted into pots that were 18 cm in diameter and 19 cm in depth. Each pot had only one plant and 5 kg of gray cinnamon soil from the field as a substrate. The physiochemical properties of the soil included an organic matter content of 6.25 g kg-1, total N content of 1.09 g kg-1, available N content of 69.45 mg g-1, available P content of 20.5 mg kg-1, available potassium content of 425.0 mg kg-1, and a pH of 8.5. Three levels of P levels were established according to P2O5 contents of 0, 50, and 200 mg kg-1 (designated CK, P1, and P4, respectively). Each treatment had 10 pots, which were placed into the field. The root crowns of alfalfa were collected on October 15, 2021. Three clones served as biological replicates. The root crowns were rinsed with clean water, wiped with absorbent paper, cut into small pieces of 3–5 mm, quickly placed into frozen tubes, stored in liquid N, and then preserved in at -80°C for further analysis (Ye et al., 2018). The temperature during the sampling period is shown in Figure S1.
2.2 Analyses of morphological and physiological indices
The phenotypic characteristics of alfalfa were measured as previously described (Wang et al., 2022). The root dry weight was measured by weighing the dried roots, which were oven-dried at 105°C for 30 min and then at 65°C for 48 h. The contents of soluble sugar were determined using anthrone–sulfuric acid colorimetry (Yemm and Willis, 1954). The contents of soluble protein were determined using Coomassie brilliant blue (Wu et al., 2013). Three plants were used as biological replicates for each treatment.
2.3 RNA extraction and sequencing
Samples treated with different levels of P were fully ground in liquid N. Total RNA was extracted using the TRIzol reagent (Life Technologies, Carlsbad, CA, USA) according to the manufacturer’s instructions. The RNA was purified, reverse-transcribed, used to construct a library, and sequenced according to the manufacturer’s instructions. To prepare Illumina RNA-seq libraries, the poly(A) mRNA was isolated from purified total RNA using biotin-Oligo (dT) magnetic beads and fragmented into small pieces using an RNA fragmentation kit according to the manufacturer’s instructions (Illumina, San Diego, CA, USA) (Hu et al., 2020). Fragmented RNA was then used as a template, and random oligonucleotides were used as primers to synthesize single-stranded complementary DNA (cDNA) by reverse transcription (Du, 2021). RNase H was then used to degrade the RNA strands, and a QiaQuick PCR extraction kit (Qiagen, Venlo, The Netherlands) was used to synthesize double-stranded cDNA (Du, 2021). The purified double-stranded cDNA was repaired by terminal repair. A tail was added and connected to the sequencing junction, and the cDNA that was approximately 200 bp was screened by ampure (AM) Pure Extraction-Purification (XP) beads. PCR amplification was conducted, and the PCR product was purified by AM Pure XP beads, which provided the cDNA library used for sequencing (Bao et al., 2020). To ensure the sequencing quality, a library quality inspection kit, the DNA 1000 assay kit (Agilent Technologies, Santa Clara, CA, USA), was used to inspect the library. The amplified fragments were sequenced using an Illumina Hi-Seq 4000 platform.
The reads obtained from the sequencing machines included raw reads that contained adapters or low-quality bases that would affect the following assembly and analysis. Thus, the reads were further filtered using an ultra-fast all-in-one FASTQ preprocessor (FASTP) (version 0.18.0) to obtain high-quality clean reads (Chen et al., 2018). An index of the reference genome was built, and paired-end clean reads were mapped to the reference genome using HISAT2. 2.4 with “-RNA-strandness RF” and other parameters set as the default (Kim et al., 2015; Chen et al., 2020). By comparing the value of the transcript on each Unigene, the value was standardized to the fragments per kilobase of transcript per million fragments mapped (FPKM) to represent the level of gene expression (Zhao et al., 2021). RNA differential expression analysis was performed using DESeq2 and edgeR (Robinson et al., 2010; Love et al., 2014). Genes with a false discovery rate below 0.05 and an absolute fold change ≥2 were considered to be differentially expressed genes (DEGs). Based on the expression information, we used R1 to conduct principal component analysis (PCA) and hierarchical clustering analyses. Gene Ontology (GO) functional and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed using GOATOOLS and KoBAS software, which can determine the primary biological functions of DEGs and identify the main biological pathways in which the DEGs are involved (Zhang et al., 2021). All the RNA sequencing data reported in this study have been deposited in National Center for Biotechnology Information (NCBI) under the sequence read archive with the accession number PRJNA909902.
The same batch of samples used in transcriptome sequencing was used for quantitative PCR (qPCR) detection. The qPCR primers and genes are shown in Table S1. The AKR4C9 gene was used as a reference. Real-time quantitative reverse transcription PCR (qRT-PCR) was conducted using a qTOWER 2.2 PCR System (Jena, Germany) and SYBR Green PCR Master Mix (TaKaRa, Shiga, Japan). The amplification program was as follows: 90 s at 95°C, and 40 cycles × (95°C, 5 s, 60°C, 15 s, and 72°C, 20 s). The dissolution curve of the amplified product was analyzed at 65°C–95°C. Each reaction was performed three times. The levels of the expression of candidate genes were measured using the 2-(ΔΔCt) method.
2.4 Metabolite extraction and ultra-high-performance liquid chromatography–quadrupole time-of-flight mass spectrometry
The samples treated with different levels of P were removed from storage at -80°C, vacuum frozen-freeze dried, and then crushed using a mixer mill (MM 400, RETSCH, Haan, Germany) with a zirconia bead for 1.5 min at 30 Hz. The sample (100 mg powder) was extracted overnight at 4°C with 1.0 ml of 70% aqueous methanol that contained 0.1 mg/L lidocaine as an internal standard. The supernatant was absorbed after centrifugation at 10,000 g for 10 min (CNWBOND Carbon-GCB SPE Cartridge, 250 mg, 3 ml; ANPEL, Shanghai, China, www.anpel.com.cn/cnw) and filtered (SCAA-104, 0.22 μm pore size; ANPEL) before liquid chromatography with tandem mass spectrometry (LC-MS/Metabolites (MS)) analysis (Chen et al., 2013). After the metabolites were extracted, 10 µl of each sample were mixed to serve as quality control (QC) samples, and 60 µl were evaluated by ultra-high-performance liquid chromatography–quadrupole time-of-flight mass spectrometry (UPLC-QTOF-MS) to monitor the stability of the instrument during the whole analytical process (Zhang et al., 2021).
The metabolites of the samples were qualitatively analyzed by mass spectrometry based on the metabolite database assembled by Guangzhou Chideo Biotechnology Co., Ltd. (Guangzhou, China). The structural analysis of some metabolites refers to the secondary Metabolites (MS) information in the existing mass spectrometry public databases, such as MassBank2, HMDB3 (Wishart et al., 2013), MoToDB4, and METLIN5 (Zheng et al., 2013). The quantification of metabolites was conducted by integrating the peak area of the mass spectra of all the substances and correcting the mass spectrometry peaks of the same metabolite in different samples to ensure the accuracy of quantification. The metabolites were analyzed by a PCA, partial least squares discriminant analysis, and orthogonal partial least squares discriminant analysis. The differences in metabolites between the different treatments were analyzed by a t-test and one-way analysis of variance (ANOVA; Zhao et al., 2020). The metabolites were mapped to the KEGG metabolic pathways for pathway and enrichment analyses.
2.5 Statistical analysis
The data were analyzed statistically with Statistical analysis System (SAS) 9.20 (SAS Institute, Cary, NC, USA) using a one-way analysis of variance (ANOVA), and multiple comparisons were performed based on the results of significance tests with the least significant difference method. Those with a P-value of t-test <0.05 and VIP ≥ 1 were considered differential metabolites between the two groups, and all the treatments had three biological replicates.
3 Results
3.1 Phosphate fertilizer affects the phenotypic and physiological characteristics of alfalfa
The phenotypes of alfalfa roots varied under different P treatments, and the growth and development of alfalfa roots were improved by the application of P (Figure 1A). The application of P reduced the depth of the root crown by 12.3% and 21.4% at P1 and P4 compared with the CK, respectively (Figure 1BI). The morphological indices of alfalfa increased during the P1 and P4 treatments compared with the CK. For example, the diameter of the root crown in the P1 and P4 treatments increased by 13.5% and 32.8%, respectively, while the number of branch roots from root crowns increased by 35.3% and 47.1%, respectively (Figure 1BII, III). The root dry weight increased by 48.8% and 59.3% at P1 and P4 compared with the CK, respectively (Figure 1BIV). The physiological indices of soluble sugar and soluble protein of the alfalfa roots increased following the application of P. The contents of soluble sugar in the P1 and P4 treatments increased by 76.6% and 54.5%, respectively, while those of the soluble protein increased by 62.5% and 92.7% compared with the CK, respectively (Figure 1BV, VI).
Figure 1 The responses of alfalfa phenotypic and physiological indices to phosphate (P) application under low temperature. (A) Alfalfa roots exposed to P application showed improved root growth. (B) Changes in the morphological and physiological indices of alfalfa roots, including (I) the depth of the root crown into the soil surface, (II) the diameter of the root crown, (III) the number of branch roots, (IV) root dry weight, (V) soluble sugar, and (IV) soluble protein. Values labeled with different lowercase letters indicate a significant difference between the P treatments at P< 0.05. CK, P1, and P4 represent 0, 50, and 200 P2O5 mg kg-1.
3.2 RNA Sequencing profiles of alfalfa roots in response to phosphate at low temperature
The analysis of DEGs was performed using edge R software under different P conditions. The results are shown in Figure 2. Compared with the CK, 49 genes were differentially expressed under the P1 conditions, with 23 genes upregulated and 26 genes downregulated. Compared with the CK, there were 224 DEGs under the P4 conditions, with 173 genes upregulated and 51 genes downregulated. Compared with P1, there were 392 DEGs under the P4 conditions, with 329 upregulated and 63 downregulated.
Figure 2 The total number of differentially expressed genes (DEGs) and upregulated and downregulated DEGs under different P treatments. DEGs, differentially expressed genes; P, phosphate.
To further understand the functions and the related biological process that the DEGs participated in, GO enrichment analyses were conducted, and the results are shown in Figure 3. The DEGs were classified into cellular components, molecular functions, and biological processes. With the cellular component, the enriched GO terms of the three comparison groups were the cell, cell part, organelle, and membrane, membrane part, and organelle part. In addition, CK-vs-P4 and P1-vs-P4 were enriched in the cell junction, extracellular region, and macromolecular complex. Within the molecular function, the enriched GO terms of three comparison groups were catalytic activity and binding. In addition, P1-vs-P4 was also enriched in antioxidant activity and transporter activity. Within the biological process, the enriched GO terms of three comparison groups were metabolic process, cellular process, single-organism process, response to stimulus, biological regulation, and localization. In addition, CK-vs-P4 and P1-vs-P4 were enriched in the developmental process, multicellular organismal process, multi-organism process, reproduction, reproductive process, signaling, immune system process, and negative and positive regulation of the biological process. These results indicated that the DEGs induced by P could be enriched into multiple GO terms, and the number of GO terms and upregulated genes increased with the increase in the application of P.
Figure 3 Gene Ontology (GO) enrichment analysis of the DEGs after cooling under different P treatments. DEGs, differentially expressed genes; GO, Gene Ontology; P, phosphate.
The pathway enrichment of the DEGs was analyzed using KEGG, and results are shown in Figure 4. The DEGs of CK-vs-P1 were primarily enriched in the biosynthesis of other secondary metabolites, carbohydrate metabolism, amino acid metabolism, and metabolism of other amino acids. The DEGs of CK-vs-P4 were enriched in not only the metabolic pathways described above but also the folding, sorting and degradation, transcription, and translation of genetic information processing and signal transduction and transport and catabolism. The DEGs of P1-vs-P4 were enriched in not only the same pathways of CK-vs-P4 but also the metabolism of terpenoids and polyketides and lipid metabolism. These results indicated that the number of KEGG pathways increased with the increase in P application.
Figure 4 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of DEGs after cooling under different P treatments. DEGs, differentially expressed genes; KEGG, Kyoto Encyclopedia of Genes and Genomes; P, phosphate.
3.3 Verification by real-time quantitative reverse transcription PCR
In addition to detecting metabolites to verify the RNA-Seq results, we used qRT-PCR to analyze six key upregulated genes to verify the results of the detection of transcriptional gene expression. The results are shown in Figure S2. The qRT-PCR verification results of the selected genes obtained using qRT-PCR were consistent with the expression data obtained by the corresponding RNA-seq, indicating that the RNA-seq expression data obtained are accurate and reliable.
3.4 Metabolite profiles of alfalfa roots in response to phosphate at low temperature
To obtain an overview of metabolic changes in response to the availability of P, a non-targeted metabolic analysis was performed. As shown in Figure 5, 24 metabolites were determined as having differential levels under P1 treatment compared with the CK. Among them, 12 metabolites increased, and 12 metabolites decreased. Compared with the CK, 12 metabolites were differentially produced under P4 treatment, and six of these metabolites increased, while six metabolites decreased. Compared with P1, 23 metabolites were determined as having differential levels under the P4 treatment. Among them, 10 metabolites increased, and 13 metabolites decreased. The primary types of metabolites that changed after the application of P were as follows: amino acids and derivatives, flavonoids, lipids, nucleotides and derivatives, organic acids, sugar, phenolic acids, and terpenoids. The types of variation in the specific metabolites and the main profile of metabolites that changed under the different P treatments are shown in Table S2.
Figure 5 The total number of differential metabolites with differential levels, upregulated and downregulated, under different P treatments.
KEGG metabolic pathway enrichment analysis classified the differential metabolites identified under low and high P treatments into metabolism, genetic information processing, environmental information processing, and human diseases (Figure 6). Within the metabolism, all three comparison groups were enriched in the biosynthesis of other secondary metabolites, amino acid metabolism, and lipid metabolism. In addition to these pathways, CK-vs-P4 and P1-vs-P4 were enriched in carbohydrate metabolism, energy metabolism, and the metabolism of other amino acids. These results indicated that the metabolic progress of alfalfa was significantly affected by P fertilizer at low temperature, and the application of increased amounts of P had a greater effect on metabolism.
Figure 6 KEGG pathway enrichment analysis of the changed metabolites under different P treatments. KEGG, Kyoto Encyclopedia of Genes and Genomes.
3.5 Integrated metabolome and transcriptome analyses of alfalfa
The integration of the DEGs and metabolites obtained from different P levels of alfalfa identified several metabolic pathways of coenrichment (Figure 7). The DEGs and metabolites of the CK and P1 had three common annotated metabolic pathways, including fatty acid degradation, phenyl propionic acid biosynthesis, and the biosynthesis of secondary metabolites. The DEGs and metabolites of CK and P4 had 12 common metabolic pathways, including the biosynthesis of secondary metabolites, biosynthesis of amino acids, sulfur metabolism, cyanoamino acid metabolism, ATP-binding cassette (ABC) transporters, cysteine and methionine metabolism, galactose metabolism, 2-oxocarboxylic acid metabolism, carbon metabolism, glyoxylate and dicarboxylate metabolism, aminoacyl-tRNA biosynthesis, and phenylalanine metabolism. The DEGs and metabolites of P1 and P4 had 15 common annotated metabolic pathways, including the biosynthesis of secondary metabolites, cysteine and methionine metabolism, biosynthesis of flavonoid, biosynthesis of isoflavone, biosynthesis of amino acid, sulfur metabolism, cyanoamino acid metabolism, ABC transporter, flavonoid and flavonol biosynthesis, 2-oxocarboxylic acid metabolism, aminoacyl-tRNA biosynthesis, phenylalanine metabolism, citrate cycle, and glyoxylate and dicarboxylate metabolism. The metabolic pathways of phenyl propionic acid biosynthesis and the biosynthesis of secondary metabolites were shared by the three groups.
Figure 7 Pathways of the simultaneous annotation of differential metabolites and DEGs for different P treatments. DEGs, differentially expressed genes; P, phosphate.
The related genes and metabolites in carbohydrate metabolism, amino acid metabolism, and the carbon metabolism pathway were combined in the three comparison groups with different rates of P applied, and the correlation between metabolites and genes was analyzed (Figure 8). Among them, N-acetyl-L-phenylalanine and L-serine are the primary metabolites that are related to amino acid metabolism, and there are 21 related genes. The level of production of N-acetyl-L-phenylalanine positively correlated with the following genes: methionine gamma-lyase (MGL) and homocysteine S-methyltransferase 3 (HMT3, MS. gene006429 and MS. gene34360). The production of L-serine positively correlated with methionine gamma-lyase (MGL), homocysteine S-methyltransferase 3(HMT3, MS. gene006429 and MS. gene34360), beta-cyanoalanine synthase (CAS1, MS. gene054826, MS. gene058651 and MS. gene023918), 1-aminocyclopropane-1-carboxylate oxidase 1 (ACO1, MS. gene98564, MS. gene75722, and MS. gene79848), beta glucosidase 11 (BGLU11, MS. gene70075, MS. gene037510, and MSTRG.1088), protein PAT1 homolog (PAT), cysteine synthase/L-3-cyanoalanine synthase (CAS2), and S-adenosylmethionine synthase (SAMS2). The primary metabolites related to carbon metabolism were lactose, isocitrate, and L-serine, and there were five related genes in which the production of lactose positively correlated with the gene RFS6 (r = 0.62). The level of production of isocitrate did not significantly correlate with the genes, and the production of L-serine positively correlated with RFS1 (r = 0.76).
Figure 8 Heatmaps of amino acid metabolism and carbohydrate metabolism with the genes of three P comparison group units. (A) Amino acid metabolism and (B) carbohydrate metabolism. The darker the color (red or blue), the stronger the correlation. * means that the correlation between metabolite and gene was significant.
4 Discussion
4.1 Effect of P fertilizer on the phenotypical and physiological indices of alfalfa during the increase in cold weather
Roots are the source of all the mineral elements required for plant growth; thus, root growth and development are highly plastic and vary substantially depending upon numerous soil factors (Vance et al., 2003). The root biomass is related to the accumulation of organic matter, which can contribute to improve the cold resistance of plants (Li et al., 2021). The root crown is the transitional plant structure located between the shoots and the root system, and it can influence photosynthate translocation, water transport, winter hardiness, and spring regrowth (Marquez-Ortiz et al., 1999). Alfalfa plants with deep and large root crowns survived the winter better than the plants with shallow and small root crowns (Marquez-Ortiz et al., 1999; Liu et al., 2015b; Wang et al., 2023). In addition, alfalfa with more lateral roots have strong cold resistance, and the possession of more lateral roots facilitates the absorption of more nutrients and water from the soil to meet the nutritional needs of overwintering alfalfa plants (Liu et al., 2015b). Previous research had proven the effect of a localized supply of soil P on the root proliferation of plants (Vance et al., 2003; Gruber et al., 2013). In this study, the cutting seedings had no obvious taproot, and the branch roots grew directly from the lower end of the root crown. The application of P fertilizer promoted the development of the root crown, widened the root crown, and increased the dry root weight. However, the increase in the number of branch roots was not apparent, and the depth of root crown was reduced. Our results showed that the application of P fertilizer can promote the development of the root system, particularly the root crown to some extent, and improve the cold resistance of alfalfa, which also facilitates the growth and development of alfalfa.
Soluble sugars and soluble proteins are intracellular osmotic regulators, which can be increased by stress, thus, improving the adaptability of plants to respond to stress (Zhao et al., 2015). Under low-temperature stress, amino acids, such as proline, arginine, and methionine, and soluble sugars, such as sucrose, glucose, and raffinose, accumulate in the root crowns, roots, and leaves of alfalfa (Zhang et al., 2015). In addition to reducing water loss, the enormous accumulation of these solutes can also prevent the cells from freezing and help to stabilize membrane integrity (Duca and Maria, 2015). As the amount of low-temperature stress increased, the contents of soluble protein and soluble sugar in the alfalfa root crown increased (Tao et al., 2009; Zhu et al., 2021) and improved the cold resistance of alfalfa. When the temperature increased, the resistance disappeared, and the sugar and protein content decreased (Zhang, 2008). P is an essential element in the composition of cells, which can participate in the composition of nuclear protein, lecithin, enzymes, ATP, and ADP among others and has an important impact on the synthesis, transport, and storage of sugar (Chen et al., 2013). P can enhance N fixation by legumes, promote fat metabolism, and enhance the stress resistance of crops. These changes occur because P can increase the content of soluble sugar and phospholipids in plants. Soluble sugar can reduce the freezing point of cell cytoplasm, and phospholipid can improve the adaptability of cells to temperature changes, thus, enhancing the cold resistance of crops (Wang et al., 2016). We found that the contents of soluble sugar and soluble protein both increased during exposure to increasing amounts of cold when P was applied, which is consistent with the findings of Shen et al. (2017).
4.2 Effect of P on the transcriptome of alfalfa roots under low-temperature stress
In general, when plants are subjected to low-temperature stress, they will change at the transcriptome level to resist the damage caused by low temperatures (Qi, 2017). Plant-specific upregulated DEGs under low-temperature stress play a more critical role in determining plant cold resistance than downregulated DEGs, and a higher percentage of upregulated genes increases the resistance of plants to cold (Zhang, 2015). In this study, we found that the P fertilizer had a substantial effect on the expression of alfalfa genes at low temperatures, and the number of DEGs and upregulated DEGs increased with the increase in P application.
Several key genes were detected in alfalfa following the application of P (Table S1). Cold regulated (COR) genes are genes that respond to low temperatures and are expressed by the regulation of a specific signal transduction pathway, and the products of COR genes are cryoprotective proteins that act by reducing membrane permeability during freezing and increasing the ability of membranes to expand during thawing (Wu et al., 2012). Thus, the constitutive expression of the COR genes could improve the freezing tolerance of plants. The novel cold-regulated genes CsCOR1 and CbCOR15 enhance dehydration tolerance in tobacco (Nicotiana benthamiana) (Li et al., 2010; Wu et al., 2012). Sangwan et al. (2001) proved that the rigidification of the plasma membrane could induce COR genes and result in cold acclimation. In this study, the expression of COR was upregulated under the application of P, which enhanced the cold resistance of alfalfa. CYP71D10 (P450) is a protein-coding gene present in plants, which have an immense variety of P450s that act on different substrates. The differential activities of P450s are believed to represent one of the mechanisms that enables certain crops species to be more tolerant of abiotic stress than other crops, particularly following herbicide treatment (Siminszky et al., 1999). Our research proved that the application of P under low temperatures could enrich the expression of a P450 that plays a vital role in cold tolerance. PAT is a photonuclear aspartate aminotransferase, which has not only the activity of ASPAT but also the activity of prephenate transaminase. It is primarily located in the cytoplasm, mitochondria, and plastids and is involved in not only plant metabolism but also the abiotic stress response pathways of plants (Wei et al., 2022). PAT reduces the accumulation of ROS and protects the structural integrity of cell membranes, which improves the abiotic tolerance stress of plants (Tian, 2021). S-adenosyl methionine synthetase (SAMS) uses L-methionine as a substrate, and ATP provides the energy to catalyze the synthesis of S-adenosylmethionine, which plays a very important role in the responses of plants to stress (Roje, 2006). Zhang et al. (2020) confirmed that the overexpression of SlSAMS1 significantly reduced the accumulation of superoxide (), hydrogen peroxide (H2O2) and malondialdehyde and enhanced the contents of abscisic acid and the enzymes that scavenge ROS, including superoxide dismutase, catalase, and ascorbate peroxidase, and plays an important role in improving the drought and salt tolerance of transgenic tomato (Solanum lycopersicum L.). Guo et al. (2014) showed that the overexpression of MfSAMS1 promoted polyamine synthesis and oxidation, which, in turn, improved the induction of protection against antioxidants by H2O2. As a result, this enhanced the tolerance to freezing and chilling stress in transgenic alfalfa plants. Our study proved that the expression of SAMS in the root crown of alfalfa treated with P was enhanced, which improved the plant’s resistance to cold. HSP70 plays an important role in plant growth and development and responds to various abiotic stresses, such as heat, drought, salinization, hormones, and other environmental perturbations (Cho and Choi, 2009; Montero-Barrientos et al., 2010; Usman et al., 2017). In heat-tolerant varieties of pepper (Capsicum annuum L.), the level of expression of CaHsp70 was rapidly upregulated upon heat stress, and the thermal stability of biofilm was simultaneously enhanced (Usman et al., 2015). Lv et al. (2022) found that the expression of JcHsp70s and its corresponding miRNAs, which are highly responsive to low temperature, is generally negatively regulated, and the interaction between JcHsp70s and miRNAs may be involved in the process of improving the cold tolerance of purging nut (Jatropha curcas L.) under low temperatures. In this study, treatment with P fertilizer also increased the expression of HSP70 in the root crown of alfalfa under low-temperature stress. In addition, there are other DEGs owing to the application of P fertilizer. These DEGs were enriched in three metabolic pathways, including the biosynthesis of other secondary metabolites, carbohydrate metabolism, and amino acid metabolism. They can also be involved in scavenging excess ROS, supplying energy, and maintaining cell osmotic pressure under low-temperature stress, which are very important for alfalfa to survive over the winter.
4.3 Effects of phosphate on the metabolism of alfalfa roots under low-temperature stress
The changes in plant metabolites under stress can reflect the ability of plants to adapt to stress. Plants accumulate different sugars and amino acids under low-temperature stress, such as trehalose, glucose, fructose, inositol, galactitol, raffinose, sucrose, putrescine, ascorbate, phenylalanine, and alanine (Jin et al., 2017). Under the influence of external substances, the expression of metabolites in plants will also change after responding to low-temperature stress (Zhou et al., 2017). Flavonoids are important secondary metabolites, which can enhance plant resistance and chemical defenses, thus, improving the tolerance of plants to abiotic stress (Zhu et al., 2007; Liu et al., 2015a). When plants are injured, secondary metabolites, such as alkaloids and polyphenolic acids, will gradually accumulate and form defense structures, thus, enhancing plant resistance (Hu et al., 2012). Terpenoids can directly or indirectly participate in a series of biological processes, such as hormone synthesis, cell membrane stability, and photosynthesis in plants. Plants can improve their ability to protect and defend themselves by regulating the content of terpenoids. Sugar metabolism is a critical metabolic cycle process in plant growth and development, and providing energy is the primary function of sugar metabolism. Secondly, sugar synthesis metabolism can also play an important role in the response of plants to abiotic stress (Sambe et al., 2015). Plants with low external levels of Pi have serious problems maintaining a balanced ratio of sugars, lipids, amino acids, organic acids, terpenes, and flavonoids, but the application of P can improve this situation (Mo et al., 2019; Ding et al., 2021; Nasr Esfahani et al., 2021). In this experiment, the application of P could regulate the types and contents of metabolites. When a small amount of P fertilizer is applied, it primarily affects the biosynthesis of secondary metabolites, while the application of more P fertilizer primarily affects the biosynthesis of secondary metabolites, amino acid metabolism, carbohydrate metabolism, and others.
4.4 Integrated metabolome and transcriptome analysis revealed the mechanism used by phospate to improve cold tolerance
The mechanism of cold resistance of alfalfa is very complex and involves morphology, physiology and biochemistry, molecular, transcription, and metabolism. Hence, it is crucial to integrate a variety of combinatorial techniques to comprehensively analyze the mechanism of the cold resistance of alfalfa (Wang and Xu, 2021). During cold acclimation, root morphology and the synthesis of compatible solutes were closely related to the cold resistance of alfalfa. Less root mass and lower concentrations of root total non-structural carbohydrates can minimize the survival of plants over the winter when the plants are subjected to low temperatures (Castonguay et al., 2006). The crown depth below the soil surface is viewed as a key morphological adaptation for forage legumes to overwinter successfully. The alfalfa varieties that have more branch rooting in general were better able to resist cold than the varieties that have a taproot or fewer types of branch roots (Johnson et al., 1998). The accumulation of soluble sugars under low temperature plays an important role in the acquisition of cold tolerance in plants (Guy, 1990). Several proteins were increasingly or newly synthesized during the acclimation of alfalfa to the cold (Mohapatra et al., 1987). Simultaneously, free amino acids were shown to accumulate in the taproots and crowns (Hendershot and Volenec, 1993; Dhont et al., 2003). The increase in free proline was notable, and the concentrations of arginine and histidine were also markedly increased (Dhont et al., 2003). In this study, the application of P fertilizer significantly changed the traits of root system architecture and increased the contents of soluble sugar and protein during the cold accumulation period. In addition, there were commonly annotated pathways shared by the transcriptome and metabolome. Thus, the application of P fertilizer influences the cold tolerance of alfalfa. With the increase in the application of P, the common metabolic pathways for the enrichment of DEGs and metabolites affected by P gradually increased. The biosynthesis of amino acids is regulated by a compound metabolic network that links N assimilation with carbon metabolism (Hans et al., 2008). Carbohydrate metabolism is the most important basic metabolism in the plant, which provides the necessary carbon frame and energy to synthesize amino acids, proteins, and nucleic acids in N metabolism. Under low-temperature stress, the application of P primarily affected the contents of L-serine, N-acetyl-L-phenylalanine, isocitrate, and lactobiose. A correlation analysis revealed that the production of L-serine positively correlated with MGL, HMT3, CAS1, BGLU11, PAT, CAS2, ACO1, SAMS2, and RFS1. N-acetyl-L-phenylalanine positively correlated with MGL and HMT3, and lactobiose positively correlated with RFS6. Among them, RFS1 and RFS6 are related to the synthesis of raffinose synthase, which is a crucial enzyme in the raffinose metabolism pathway. This pathway can catalyze the reaction of inositol galactoside and sucrose to raffinose, and they can also catalyze the reaction of inositol galactoside and raffinose to synthesize stachyose. Its expression can be regulated by various stress response signals and transcription factors (Nishizawa et al., 2008; Cho et al., 2010). SAMS2 facilitates the synthesis of S-adenosyl methionine synthetase and can be induced by a variety of stress factors; its overexpression can significantly improve the resistance to plant stress (Guo et al., 2014). The overexpression of PAT increased the contents of proline and soluble sugar in transgenic Arabidopsis thaliana and improved the cold, drought, and saline–alkali tolerance of plants (Yuan et al., 2016). Other genes also play essential roles in the synthesis and metabolism of amino acids and carbohydrates. Thus, the moderate application of P fertilizer can improve the cold resistance of alfalfa by regulating the synthesis and metabolism of amino acids and sugars.
5 Conclusion
Low-temperature stress affects the growth and development of alfalfa, leading to morphological, physiological, metabolic, and molecular changes and, thus, limits the overwintering and production of alfalfa. The application of P can improve the cold resistance and forage yield of alfalfa. In this study, we analyzed the changes of alfalfa under low-temperature stress by transcriptomic and metabolomic approaches after the application of P. We found that the application of P fertilizer improved the root structure and increased the content of soluble sugar and soluble protein in the root crown, and there were 49 DEGs (23 genes upregulated) and 24 metabolites (12 upregulated) following 50 mg kg-1 of P applied and 224 DEGs (173 genes upregulated) and 12 metabolites (6 upregulated) following 200 mg kg-1 of P applied, respectively, compared with the CK. These genes and metabolites were significantly enriched in the biosynthesis of other secondary metabolites, carbohydrate metabolism, and amino acid metabolism pathways. The integration of the transcriptomic and metabolomic analyses showed that P could affect the biosynthesis of N-acetyl-L-phenylalanine, L-serine, lactose, and isocitrate and the expression of related genes to regulate the cold tolerance of alfalfa. Our findings could contribute to a deeper understanding of the mechanism that alfalfa uses to respond to cold tolerance and lay a theoretical foundation for breeding highly P-efficient alfalfa.
Data availability statement
The data presented in the study are deposited in the NCBI repository, accession number PRJNA909902. The data has been released.
Author contributions
YW, LY and JX designed the experiment. YW, ZS and QW collected the samples in the field, and YW and ZS conducted laboratory analyses. YW analyzed the data, interpreted the results, and wrote the manuscript. All authors reviewed and edited the manuscript and approved the submitted version.
Funding
This research was financially supported by the Inner Mongolia Science and Technology Program (2019GG260), Achievement Transformation Project in Science and Technology of Inner Mongolia (2019CG041), High level talent introduction project of Inner Mongolia University (10000-22311201/016), Alfalfa breeding for high quality and industry demonstration (2022JBGS0020). The funding body had no contribution to the study design, data collection, analysis, interpretation of the data, the writing of this article, or the decision to submit for publication.
Acknowledgments
The authors are grateful for the sequencing platform and/or bioinformation analysis of Gene Denovo Biotechnology Co., Ltd (Guangzhou, China). We would like to thank MogoEdit (https://www.mogoedit.com) for its English editing during the preparation of this manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2023.1100601/full#supplementary-material
Footnotes
- ^ http://www.r-project.org/
- ^ http://www.massbank.jp/
- ^ http://www.hmdb.ca/
- ^ http://www.ab.wur.nl/moto/
- ^ http://metlin.scripps.edu/index.php
References
Acharya, J. P., Lopez, Y., Gouveia, B. T., de Bem Oliverira, I., Resende, M. F. R., Munoz, P. R., et al. (2020). Breeding alfalfa (Medicago sativa l.) adapted to subtropical agroecosystems. Agronomy 10 (5), 1–13. doi: 10.3390/agronomy10050742
Agarrwal, R., Padmakumari, A. P., Bentur, J. S., Nair, S. (2016). Metabolic and transcriptomic changes induced in host during hypersensitive response mediated resistance in rice against the Asian rice gall midge. Rice 9, 5. doi: 10.1186/s12284-016-0077-6
Bao, G., Tang, W., An, Q., Liu, Y., Tian, J., Zhao, N., et al. (2020). Physiological effects of the combined stresses of freezing-thawing, acid precipitation and deicing salt on alfalfa seedlings. BMC Plant Biol. 20 (1), 204. doi: 10.1186/s12870-020-02413-4
Bertrand, A., Bipfubusa, M., Claessens, A., Rocher, S., Castonguay, Y. (2017). Effect of photoperiod prior to cold acclimation on freezing tolerance and carbohydrate metabolism in alfalfa (Medicago sativa l.). Plant Sci. 264, 122–128. doi: 10.1016/j.plantsci.2017.09.003
Bhosale, R., Giri, J., Pandey, B. K., Giehl, R. F. H., Hartmann, A., Traini, R., et al. (2018). A mechanistic framework for auxin dependent arabidopsis root hair elongation to low external phosphate. Nat. Commun. 9 (1), 1409. doi: 10.1038/s41467-018-03851-3
Castonguay, Y., Laberge, S., Brummer, E. C., Volenec, J. J. (2006). Alfalfa winter hardiness: A research retrospective and integrated perspective. Adv. Agric. 90, 203–265. doi: 10.1016/s0065-2113(06)90006-6
Chen, W., Gong, L., Guo, Z., Wang, W., Zhang, H., Liu, X., et al. (2013). A novel integrated method for large-scale detection, identification, and quantification of widely targeted metabolites: Application in the study of rice metabolomics. Mol. Plant 6 (6), 1769–1780. doi: 10.1093/mp/sst080
Chen, H., Zeng, Y., Yang, Y., Huang, L., Tang, B., Zhang, H., et al. (2020). Allele-aware chromosome-level genome assembly and efficient transgene-free genome editing for the autotetraploid cultivated alfalfa. Nat. Commun. 11 (1), 2494. doi: 10.1038/s41467-020-16338-x
Chen, S., Zhou, Y., Chen, Y. (2018). Fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34 (17), 884–890. doi: 10.1093/bioinformatics/bty560
Cheng, G. (2008). Effect of phosphorus on yield、quality、nutrient absorption and seedling cold tolerance of watermelon (Wuhan: Huazhong Agricultural University).
Cho, E. K., Choi, Y. J. (2009). A nuclear-localized HSP70 confers thermoprotective activity and drought-stress tolerance on plants. Biotechnol. Lett. 31 (4), 597–606. doi: 10.1007/s10529-008-9880-5
Cho, S., Kang, E., Kim, M., Yoo, S., Im, Y., Kim, Y., et al. (2010). Jasmonate-dependent expression of a alactinol synthase gene is involved in priming of systemic fungal resistance in arabidopsis thaliana. Botany 88, 452–461. doi: 10.1139/B10-009
Cunningham, S. M., Volenec, J. J., Teuber, L. R. (1998). Plant survival and root and bud composition of alfalfa populations selected for contrasting fall dormancy. Crop Sci. 38 (4), 962–969. doi: 10.2135/cropsci1998.0011183X003800040014x
Dhont, C., Castonguay, Y., Nadeau, P., Bélanger, G., Chalifour, F. P. (2003). Alfalfa root nitrogen reserves and regrowth potential in response to fall harvests. Crop Sci. 43 (1), 181–194. doi: 10.2135/cropsci2003.1810
Diatta, A. A., Min, D., Jagadish, K. S. (2021). Drought stress responses in non-transgenic and transgenic alfalfa–current status and future research directions. Adv. Agron. 170, 2–66. doi: 10.1016/bs.agron.2021.06.002
Ding, N., Huertas, R., Torres-Jerez, I., Liu, W., Watson, B., Scheible, W. R., et al. (2021). Transcriptional, metabolic, physiological and developmental responses of switchgrass to phosphorus limitation. Plant Cell Environ. 44 (1), 186–202. doi: 10.1111/pce.13872
Dou, Y. M. (2011). Research progress on the mechanism of cold resistance of alfalfa in China. Heilongjiang Agric. Sci. 07, 146–148.
Du, S. (2021). Metabolic mechanism of cold stress in rice(Oryza sativa l.) seedlings based on metabolome and transcriptome (Nanchang: Nanchang University).
Duca, Maria (2015). “Plant physiology || physiology of plant resistance to unfavorable environmental factors,” in Biological and medical physics, biomedical engineering. Ed. Gerstman., B. S. (New York, NY: Springer), 271–308.
Edde, P. A. (2022). “Arthropod pests of alfalfa (Medicago sativa l.),” in Field crop arthropod pests of economic importance, ed. Edde., P. A. (Altria Client Services LLC, Richmond, VA, United States: Academic Press), 684–750.
Gaafar, Y., Richert-Pöggeler, K., Maaß, C., Vetten, H. J., Ziebell, H. (2019). Characterisation of a novel nucleorhabdovirus infecting alfalfa (Medicago sativa). Virol. J. 16 (1), 55. doi: 10.1186/s12985-019-1147-3
Goplen, B. P. (1987). Growing and managing alfalfa in Canada (Ottawa K1A 0C7: Communications Branch, Agriculture Canada).
Gruber, B. D., Giehl, R. F., Friedel, S., von Wiren, N. (2013). Plasticity of the arabidopsis root system under nutrient deficiencies. Plant Physiol. 163 (1), 161–179. doi: 10.1104/pp.113.218453
Guo, Z., Tan, J., Zhuo, C., Wang, C., Xiang, B., Wang, Z. (2014). Abscisic acid, H2O2 and nitric oxide interactions mediated cold-induced s-adenosylmethionine synthetase in medicago sativa subsp. falcata that confers cold tolerance through up-regulating polyamine oxidation. Plant Biotechnol. J. 12 (5), 601–612. doi: 10.1111/pbi.12166
Guy, C. (1990). Cold acclimation and freezing stress tolerance: Role of protein metabolism. Annu. Rev. Plant Physiol. Plant Molec. Biol. 41, 187–224. doi: 10.1146/annurev.pp.41.060190.001155
Hans, J. B., Henry, N., Norman, G. L. (2008). Bioengineering and molecular biology of plant pathways (Holland: ELsevier Publishing Group press).
Hendershot, K. L., Volenec, J. J. (1993). Taproot nitrogen accumulation and use in overwintering alfalfa (Medicago sativa l.). J. Plant Physiol. 141 (1), 68–74. doi: 10.1016/s0176-1617(11)80853-9
Hong, F. Z., Wu, S. Y., Yu, K. F., Bai, Y. H., Cheng, D. (1987). Study on introducing alfalfa with root tiller type. Chin. pratacult. Sci. 5, 1–4.
Hu, W., Jiang, J., Yang, H., Liu, C., He, Y. (2012). Effect of jasmonic acid methyl ester treatment on the physiological and biochemical reactions of fresh-cut apple. Sci. Technol. Food Ind. 33 (16), 338–341+346. doi: 10.13386/j.issn1002-0306.2012.16.055
Hu, T., Sun, X. Y., Zhao, Z. J., Amombo, E., Fu, J. M. (2020). High temperature damage to fatty acids and carbohydrate metabolism in tall fescue by coupling deep transcriptome and metabolome analysis. Ecotoxicol. Environ. Saf. 203, 110943. doi: 10.1016/j.ecoenv.2020.110943
Janz, D., Behnke, K., Schnitzler, J. P., Kanawati, B., Schmitt-Kopplin, P., Polle, A. (2010). Pathway analysis of the transcriptome and metabolome of salt sensitive and tolerant poplar species reveals evolutionary adaption of stress tolerance mechanisms. BMC Plant Biol. 10, 150. doi: 10.1186/1471-2229-10-150
Jin, J., Zhang, H., Zhang, J., Liu, P., Chen, X., Li, Z., et al. (2017). Integrated transcriptomics and metabolomics analysis to characterize cold stress responses in nicotiana tabacum. BMC Genomics 18 (1), 496. doi: 10.1186/s12864-017-3871-7
Johnson, L. D., Marquez-Ortiz, J. J., Barnes, D. K., Lam, J. F. S. (1996). Inheritance of root traits in alfalfa. Crop Sci. 36, 1482–1487. doi: 10.2135/cropsci1996.0011183X003600060011x
Johnson, L. D., Marquez-Ortiz, J. J., Lamb, J., Barnes, D. (1998). Root morphology of alfalfa plant introductions and cultivars. Crop Sci. 38, 497–502. doi: 10.2135/cropsci1998.0011183X003800020037x
Kim, D., Langmead, B., Salzberg., S. L. (2015). HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 12 (4), 357. doi: 10.1038/nmeth.3317
Kontunen-Soppela, S., Lankila, J., Lähdesmäki, P., Laine, K. (2002). Response of protein and carbohydrate metabolism of scots pine seedlings to low temperature. J. Plant Physiol. 159 (2), 175–180. doi: 10.1078/0176-1617-00538
Li, X. W., Feng, Z. G., Yang, H. M., Zhu, X. P., Liu, J., Yuan, H. Y. (2010). A novel cold-regulated gene from camellia sinensis, CsCOR1, enhances salt- and dehydration-tolerance in tobacco. Biochem. Biophys. Res. Commun. 394 (2), 354–359. doi: 10.1016/j.bbrc.2010.03.011
Li, M., Shi, X., Guo, C., Lin, S. (2016). Phosphorus deficiency inhibits cell division but not growth in the dinoflagellate amphidinium carterae. Front. Microbiol. 7. doi: 10.3389/fmicb.2016.00826
Li, Z., Wan, L., Li, S., Li, X., He, F., Tong, Z. (2021). Plastic response of medicago sativa l. root system traits and cold resistance to simulated rainfall events. PeerJ 9, e11962. doi: 10.7717/peerj.11962
Lin, Z. H., Chen, R. B., Guo, S. P. (2010). Research progress on physiological adaptability of plants to phosphorus deficiency. Crops 05), 5–9. doi: 10.16035/j.issn.1001-7283.2010.05.003
Liu, D., Li, M., Liu, Y., Shi, L. (2020). Integration of the metabolome and transcriptome reveals the resistance mechanism to low nitrogen in wild soybean seedling roots. Environ. Exp. Bot. 175, 104043. doi: 10.1016/j.envexpbot.2020.104043
Liu, Z., Li, X., Wang, Z., Sun, Q. (2015a). Contrasting strategies of alfalfa stem elongation in response to fall dormancy in early growth stage: The tradeoff between internode length and internode number. PloS One 10 (8), e0135934. doi: 10.1371/journal.pone.0135934
Liu, Z. Y., Xi-Liang, L. I., Feng, L. I., Tao, Y., Liu, L., Wang, Z. L., et al. (2015b). Response of alfalfa root traits to fall dormancy and its effect on winter hardiness. Sci. Agricult. Sin. 48 (9), 1689–1701. doi: 10.3864/j.issn.0578-1752.2015.09.03
Love, M. I., Huber, W., Anders, S. (2014). Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15 (12), 550. doi: 10.1186/s13059-014-0550-8
Lv, Y., Wu, D., Kong, C., Yang, Y., Gong, M. (2022). Identification and interaction analysis of the Hsp70 gene family and the corresponding miRNAs and their possible roles in low temperature adaptation in jatropha curcas. Plant Physiol. J. 58 (7), 1–16. doi: 10.13592/j.cnki.ppj.100167
Marquez-Ortiz, J. J., Lamb, J. F. S., Johnson, L. D., Barnes, D. K., Stucker, R. E. (1999). Heritability of crown traits in alfalfa. Crop Sci. 39, 38–43. doi: 10.2135/cropsci1999.0011183X003900010006x
Mo, X., Zhang, M., Liang, C., Cai, L., Tian, J. (2019). Integration of metabolome and transcriptome analyses highlights soybean roots responding to phosphorus deficiency by modulating phosphorylated metabolite processes. Plant Physiol. Biochem. 139, 697–706. doi: 10.1016/j.plaphy.2019.04.033
Mohapatra, S. S., Poole, R. J., Dhindsa, R. S. (1987). Changes in protein patterns and translatable messenger RNA populations during cold acclimation of alfalfa. Plant Physiol. J. 84 (4), 1172–1176. doi: 10.1104/pp.84.4.1172
Moieni-Korbekandi, Z., Karimzadeh, G., Sharifi, M. (2013). Evaluation of total soluble protein and antioxidant activities in two spring cultivars of canola (Brassica napus l.) in response to low temperature. Int. J. Agric. Crop Sci. 5 (4), 401–409.
Montero-Barrientos, M., Hermosa, R., Cardoza, R. E., Gutiérrez, S., Nicolás, C., Monte, E. (2010). Transgenic expression of the trichoderma harzianum hsp70 gene increases arabidopsis resistance to heat and other abiotic stresses. J. Plant Physiol. 167 (8), 659–665. doi: 10.1016/j.jplph.2009.11.012
Nasr Esfahani, M., Inoue, K., Nguyen, K. H., Chu, H. D., Watanabe, Y., Kanatani, A., et al. (2021). Phosphate or nitrate imbalance induces stronger molecular responses than combined nutrient deprivation in roots and leaves of chickpea plants. Plant Cell Environ. 44 (2), 574–597. doi: 10.1111/pce.13935
Nishizawa, A., Yabuta, Y., Shigeoka, S. (2008). Galactinol and raffinose constitute a novel function to protect plants from oxidative damage. Plant Physiol. 147 (3), 1251–1263. doi: 10.1104/pp.108.122465
Putnam, D. H., Russelle, M. P., Orloff, S. B., Kuhn, J., Fitzhugh, L., Godfrey, L. (2001). Alfalfa, wildlife and the environment: The importance and benefit of alfalfa in the 21st century. Minnesota: California Alfalfa and Forage Assiciation.
Qi, X. (2017). The Analysis of the Differentially Expressed Genes in Alfalfa under Cold Stress at Transcriptome Level. Doctor, Chinese Academy of Agricultural Science.
Raghothama, K., Karthikeyan, A. (2005). “Phosphate acquisition,” in Root physiology: From gene to function. Eds. Lambers, H., Colmer, T. D. (New Y ork: Springer).
Robinson, M. D., McCarthy, D. J., Smyth, G. K. (2010). edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 26 (1), 139–140. doi: 10.1093/bioinformatics/btp616
Rocher, S., Jean, M., Castonguay, Y., Belzile, F. (2015). Validation of genotyping-By-Sequencing analysis in populations of tetraploid alfalfa by 454 sequencing. PloS One 10 (6), e0131918. doi: 10.1371/journal.pone.0131918
Roje, S. (2006). S-Adenosyl-L-methionine: Beyond the universal methyl group donor. Phytochemistry 67 (15), 1686–1698. doi: 10.1016/j.phytochem.2006.04.019
Sambe, M. A., He, X., Tu, Q., Guo, Z. (2015). A cold-induced myo-inositol transporter-like gene confers tolerance to multiple abiotic stresses in transgenic tobacco plants. Physiol. Plant 153 (3), 355–364. doi: 10.1111/ppl.12249
Sangwan, V., Foulds, I., Singh, J., Dhindsa, R. S. (2001). Cold-activation of brassica napus BN115 promoter is mediated by structural changes in membranes and cytoskeleton, and requires Ca2+ influx. Plant J. 27 (1), 1–12. doi: 10.1046/j.1365-313x.2001.01052.x
Shen, X., Zhang, Y., Wang, X., Li, Y., Guo, Y., Yu, H. (2017). Effects of phosphorus and potassium fertilizer on the cold resistance of alfalfa roots and necks in sandy soil. J. Inner Mongolia Univ. Natl. 32 (2), 116–121.
Siminszky, B., Corbin, F. T., Ward, E. R., Fleischmann, T. J., Dewey, R. E. (1999). Expression of a soybean cytochrome P450 monooxygenase cDNA in yeast and tobacco enhances the metabolism of phenylurea herbicides. Proc. Natl. Acad. Sci. U.S.A. 96, 1750–1755. doi: 10.1073/pnas.96.4.1750
Summers, C. G., Putnam, D. H. (2008). Irrigated alfalfa management for Mediterranean and desert zones (Irvine, CA: University of California Alfalfa & Forages).
Tao, Y., Yu, Z., Sun, Q., Zhao, S. (2009). Study on physiological adaptability of alfalfa for cold resistance. Pratacult. Sci. 26 (9), 151–155.
Terry, N., Tao, I. M. (1991). Nutrient and photosynthesis:iron and phosphorus as case studies (University Printing House, Shaftesbury Rd, Cambridge CB2 8BS: Cambridge University Press).
Tian, Z. (2021). Functional study on salt tolerance of PAT1 gene from betuia platyphylla (Haerbin: Northest Forestry University).
Tian, Z. J., Li, J. P., Yang, F. (2017). Progress in crop phosphorus utilization under abiotic stresses. Chin. J. Ecol. 36 (8), 2336–2342. doi: 10.13292/j.1000-4890.201708.002
Trischuk, R. G., Schilling, B. S., Low, N. H., Gray, G. R., Gusta, L. V. (2014). Cold acclimation, de-acclimation and re-acclimation of spring canola, winter canola and winter wheat: The role of carbohydrates, cold-induced stress proteins and vernalization. Environ. Exp. Bot. 106, 156–163. doi: 10.1016/j.envexpbot.2014.02.013
Usman, M. G., Rafii, M. Y., Ismail, M. R., Malek, M. A., Latif, M. A. (2015). Expression of target gene Hsp70 and membrane stability determine heat tolerance in chili pepper. J. Am. Soc. Hortic. Sci. Am. Soc. Hortic. Sci. 140 (2), 144–150. doi: 10.21273/JASHS.140.2.144
Usman, M. G., Rafii, M. Y., Martini, M. Y., Yusuff, O. A., Miah, G. (2017). Molecular analysis of Hsp70 mechanisms in plants and their function in response to stress. Biotechnol. Genet. Eng. Rev. 33 (1), 1–14. doi: 10.1080/02648725.2017.1340546
Vance, C. P., Uhde-Stone, C., Allan, D. L. (2003). Phosphorus acquisition and use critical adaptations by plants for securing a nonrenewable resource. New Phytol. 157, 423–447. doi: 10.1046/j.1469-8137.2003.00695.x
Viands, D. R. (1988). Variability and selection for characters associated with root regeneration capability in alfalfa. Crop Sci. 28 (2), 232–236. doi: 10.1007/BF02871603
Wang, X. (2021a). Effect of fertilization and reseeding on the productivity, forage quality and soil properties of degraded alfalfa grassland (Xining: Ningxia University).
Wang, X. L. (2021b). Identification of cold tolerance and screening of cold tolerance germplasm of alfalfa (Hohhot: Inner Mongolia Agricultural University).
Wang, R., Xu, W. (2021). Development characteristics and trend of alfalfa industry in China. J. Agric. Sci. Technol. 23 (12), 7–12. doi: 10.13304/j.nykjdb.2021.0940
Wang, X., Yang, Z., Lai, Y.-c., Li, H., Zhong, P., Xu, Y.-x., et al. (2023). Effect of root traits of medicago sativa lines with fall dormancy on overwintering. Acta Pratacult. Sin. 32 (1), 144–153. doi: 10.11686/cyxb2022016
Wang, D., Yang, J., Rongsheng, L. I., Zou, W., Yin, G. (2016). Physiological response of phoebe bournei seedlings with different phosphorus supply levels to natural low temperature. J. South China Agric. Univ. 37 (2), 101–106.
Wang, Y., Zhang, J., Yu, L., Xu, Z., Samac, D. A. (2022). Overwintering and yield responses of two late-summer seeded alfalfa cultivars to phosphate supply. Agronomy 12, 327. doi: 10.3390/agronomy12020327
Wechsler, S. J., Milkove, D., U.S.D.O. AGRICULTURE (2017). Genetically modified alfalfa production in the united states (United States: Economic Research Service).
Wei, Z., Liu, C.-Y., Chen, H.-W., Li, L., Sun, L.-Q., Han, X.-S., et al. (2022). Analysis of ASPAT gene family based on drought-stressed transcriptome sequencing in vicia faba l. Acta Agro. Sin. 12, 1–12.
Wishart, D. S., Timothy, J., Chi, G. A., Michael, W., Craig, K., Liu, Y., et al. (2013). HMDB 3.0–the human metabolome database in 2013. Nucleic Acids Res. D1), D801.
Wu, H., Zhang, J., Shi, J., Fan, Z., Aliyan, R., Zhang, P. (2013). Physiological responses of cotton seedlings under low temperature stress. Acta Botanica Boreali-Occidentalia Sin. 33 (1), 74–82.
Wu, L., Zhou, M., Shen, C., Liang, J., Lin, J. (2012). Transgenic tobacco plants over expressing cold regulated protein CbCOR15b from capsella bursa-pastoris exhibit enhanced cold tolerance. J. Plant Physiol. 169 (14), 1408–1416. doi: 10.1016/j.jplph.2012.05.016
Xin, W., Zhang, L., Zhang, W., Gao, J., Yi, J., Zhen, X., et al. (2019). An integrated analysis of the rice transcriptome and metabolome reveals root growth regulation mechanisms in response to nitrogen availability. Int. J. Mol. Sci. 20 (23), 5893. doi: 10.3390/ijms20235893
Xu, Z., York, L. M., Seethepalli, A., Bucciarelli, B., Cheng, H., Samac, D. A. (2022). Objective phenotyping of root system architecture using image augmentation and machine learning in alfalfa (Medicago sativa l.). Plant Phenomics 2022, 9879610. doi: 10.34133/2022/9879610
Yan, Z., Wang, H., Kou, X., Wu, C., Fan, G., Li, T., et al. (2022). Metabolomics analysis reveals that MeJA treatment induces postharvest blueberry resistance to botrytis cinerea. Postharvest Biol. Technol. 194, 112075. doi: 10.1016/j.postharvbio.2022.112075
Yang, C. L., Duan, R. J., Wu, X. X., Qi, C. Y., Ma, Y. H., Xiong, H. Y. (2021). Genome-wide identification, sequence variation, and expression of the GPAT gene family in Medicago truncatula. Pratacult. Sci. 38 (10), 1966–1974.
Ye, C. Y., Zheng, C. P., Ying, W. W., Weng, S. S. (2018). Up-regulation of microRNA-497 inhibits the proliferation, migration and invasion but increases the apoptosis of multiple myeloma cells through the MAPK/ERK signaling pathway by targeting raf-1. Cell Cycle 17 (24), 2666–2683. doi: 10.1080/15384101.2018.1542895
Yemm, E. W., Willis, A. J. (1954). The estimation of carbohydrates in plant extracts by anthrone. Biochem. J. 57 (3), 508–514. doi: 10.1042/bj0570508
Yuan, Y., Fang, L., Karungo, S. K., Zhang, L., Gao, Y., Li, S., et al (2016). Overexpression of VaPAT1, a GRAS transcription factor from Vitis amurensis, confers abiotic stress tolerance in Arabidopsis. Plant Cell Rep. 35 (3), 655–666. doi: 10.1007/s00299-015-1910-x
Yue, Z. H., Shen, Y. H., Chen, Y. J., Liang, A. W., Chu, C. W., Chen, C., et al. (2019). Microbiological insights into the stress-alleviating property of an endophytic bacillus altitudinis WR10 in wheat under low-phosphorus and high-salinity stresses. Microorganisms 7 (11), 508. doi: 10.3390/microorganisms7110508
Zhang, X. (2008). Physiological adaptability of different alfalfa varieties to low temperature conditions in autumn and winter (Hohhot: Inner Mongolia Agriculture University).
Zhang, H. L. (2015). Transcription analysis of summer squash in response to cold stress and the development of SSR molecular marker (Jinzhong: Shanxi Agricultural University).
Zhang, X., Bao, Z., Gong, B., Shi, Q. (2020). S-adenosylmethionine synthetase 1 confers drought and salt tolerance in transgenic tomato. Environ. Exp. Bot. 179 (17), 104226. doi: 10.1016/j.envexpbot.2020.104226
Zhang, S., Shi, Y., Cheng, N., Du, H., Fan, W., Wang, C. (2015). De novo characterization of fall dormant and nondormant alfalfa (Medicago sativa l.) leaf transcriptome and identification of candidate genes related to fall dormancy. PloS One 10 (3), e0122170. doi: 10.1371/journal.pone.0122170
Zhang, Y., Zhao, P., Yue, S., Liu, M., Qiao, Y., Xu, S., et al. (2021). New insights into physiological effects of anoxia under darkness on the iconic seagrass zostera marina based on a combined analysis of transcriptomics and metabolomics. Sci. Total Environ. 768, 144717. doi: 10.1016/j.scitotenv.2020.144717
Zhao, Y., Li, M.-C., Konaté, M. M., Chen, L., Das, B., Karlovich, C., et al. (2021). TPM, FPKM, or normalized counts? a comparative study of quantification measures for the analysis of RNA-seq data from the NCI patient-derived models repository. J. Trans. Med. 19 (1), 269. doi: 10.1186/s12967-021-02936-w
Zhao, Y., Wei, Y., Zou, Y. (2015). Physiological index variation and its cultivar differences of seedling leaves and roots of early cropping rice under lower temperature stress. J. Nucl. Agric. Sci. 29 (4), 792–798. doi: 10.11869/j.issn.100-8551.2015.04.0792
Zhao, C., Zhang, R., Niu, K., Zhu, R., Wang, Y., Ma, X., et al. (2020). Metabonomics study of qinghai wild poa pratensis in response to low temperature stress. Acta Agrestia Sinca 28 (4), 904–914.
Zheng, J., Schultz, A. W., Wang, J., Johnson, C. H., Yannone, S. M., Patti, G. J., et al. (2013). Liquid chromatography quadrupole time-of-flight mass spectrometry characterization of metabolites guided by the metlin database. Nat. Protoc. 8 (3), 451–460. doi: 10.1038/nprot.2013.004
Zhou, L., Li, Y., Wang, W., Zhong, S. (2017). Research progress in the metabolomics for plants response to temperature stress. JournalofShanxi Agricult. Sci. 45 (2), 317–320. doi: 10.3969/j.issn.1002-2481.2017.02.42
Zhu, H., Hu, H., Lu, C., Li, X. (2007). Progresses on flavonoid metabolism in plants and it′s regulation. J. Xiamen Univ. (Natural Science) 46 (1), 136–143.
Keywords: alfalfa, transcriptomics, metabolomics, cold tolerance, phosphorus
Citation: Wang Y, Sun Z, Wang Q, Xie J and Yu L (2023) Transcriptomics and metabolomics revealed that phosphate improves the cold tolerance of alfalfa. Front. Plant Sci. 14:1100601. doi: 10.3389/fpls.2023.1100601
Received: 17 November 2022; Accepted: 17 February 2023;
Published: 09 March 2023.
Edited by:
Jing Zhang, Nanjing Agricultural University, ChinaReviewed by:
Boon Chin Tan, University of Malaya, MalaysiaQiguo Sun, Jiangsu Vocational College of Agriculture and Forestry, China
Copyright © 2023 Wang, Sun, Wang, Xie and Yu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jihong Xie, xiejihong223@163.com; Linqing Yu, linqing_yu@126.com