3.1 Determination results of Alfalfa physiological indexes and soil enzyme activity
The results show that the growth index and nutrition index of alfalfa in the experimental group are higher than those in the control group.The biological fertilizer made alfalfa grow stronger and leaves bigger obviously.And the nutrient content of alfalfa in the treatment group is higher than that in the control group(Figure 1a,Table 1).
Table 1
Effect of biofertilizer on growth and nutritional index of Alfalfa.
|
|
|
Leaf Area (cm2)
|
Stem diameter (mm)
|
Soluble protein (mg/g)
|
Crude Protein (%)
|
Soluble sugars (%)
|
Chlorophyll a (mg/L)
|
Chlorophyll b (mg/L)
|
Carotenoids (mg/L)
|
Peroxidase (U/Fw · min)
|
MDA content (μmol/gFW)
|
|
|
30ZS
|
1.8791
|
2.482
|
4.2
|
20.58
|
0.6961
|
13.1256
|
15.1209
|
6.4786
|
7990
|
0.0030838
|
|
30HX
|
1.3881
|
2.133
|
4.045
|
17.59
|
0.5808
|
6.5448
|
9.5031
|
2.5144
|
7230
|
0.004128
|
|
In addition, this study also measured the soil enzyme activity in the experimental site.During the whole study period, other field management methods were exactly the same except adding biological fertilizer to the experimental group. To describe the influence of biological fertilizer on soil, the soil enzyme activity was measured. Soil enzyme activity of the treatment group is obviously increased, compared with the control group(Table 2).Therefore, biological fertilizer can effectively change the pollution of chemical fertilizer to soil.
Table 2
Effects of biofertilizers on soil.
|
|
|
Catalase (0.1 mol/L KMnO4) /(h · g)
|
Urease (Ure)
|
Sucrase (XII)
|
FDA Hydrolase (μg/Ml)
|
Acid phosphatase (mg/g)
|
Neutral phosphatase (mg/g)
|
Alkaline phosphatase (mg/g)
|
|
|
30ZS
|
0.015323
|
20.56145
|
41.90987
|
330.6663
|
4.227333
|
2.892667
|
1.735333
|
|
|
30HX
|
0.014118
|
20.05166
|
41.3396
|
279.8723
|
3.088667
|
2.192667
|
1.007333
|
|
|
3.2 Transcriptome sequencing, assembly and differential gene expression analysis results.
Alfalfa transcriptome under different fertilization conditions was sequenced by Illumina 2000. After removal of low-mass readings, 42.27Gb of clean data were attained. For each sample, clean data was 6.10 Gb and the proportion of Q30 bases was 92.94%or more (Table S1). A total of 61,913 Unigene were obtained after assembly. Among them, there are 18,933 Unigene with length over 1kb.
By analyzing the transcriptome of alfalfa, we found 2612 DEG, of which 1085 were up-regulated, and 1527 were down-regulated (Figure 1). Then, PCA was utilized to measure the differences in gene expression across the groups.The treatment group (30ZS PCA1) was obviously separated from the control group(30HZ), and these findings were consistent with those from the heat map(Figure 1b-d, Table S2).
We used GO classification for the DEGs to resolve the appropriate subordinate categories for the responsive genes. The matched deg were thus apportioned into three functional groups: biological processes, cell components, and molecular function (fig. 2a). Within biological processes, the most common genes were associated with cellular and metabolic processes. For cell components, genes were commonly associated with cell parts and cells. In the molecular functional group, DEGs primarily belong to the binding and catalytic activity subgroup. To ascertain the functional biological pathway rich in DEGs, the KEGG pathway database was utilized (Figure 2b). Primary enrichment pathways are plant hormone signal transduction, phenylpropanoid biosynthesis, brassinosteroid biosynthesis, Stilbenoid, diarylheptanoid and gingerol biosynthesis, Phenylalanine Metabolism and Flavonoid biosynthesis(Figure 2c).
Through KEGG enrichment of differential genes, it was found that plant signal transduction pathways were obviously enriched, in which almost all signal pathways showed up-regulation or down-regulation of genes, especially auxin signal transduction, brassinosteroid and gibberellin signal transduction pathways, and brassinosteroid organisms and pathways were also obviously enriched in KEGG enrichment of differential genes. Although some genes in ABA and ethylene signal transduction pathway are up-regulated, they are mainly concentrated in the negative regulation genes of ABA signal transduction pathway, such as PP2C gene in ABA signal transduction pathway, ETR gene and EBF1/2 gene in ethylene signal transduction pathway. In addition, ABA and cytokinin signal transduction related genes are obviously down-regulated, such as SnRK2 and ABF genes in ABA signal transduction pathway, CRE1 and B-ARR genes in cytokinin signal transduction pathway. Compared with the control transcriptome, some genes were differentially expressed only in treatment group.The genes specifically expressed by these treatment groups can be different biological pathways.Among them, 13 genes belong to phenylalanine metabolism, which is one of the most important plant secondary metabolic pathways, where it is critical in growth, development, and environmental interaction(Table S3).In addition, 6 genes and 17 genes belong to stilbenes, diarylheptanoic acid and gingerol biosynthetic pathway and flavonoid biosynthetic pathway, respectively.Phenols and flavonoids produced by stilbenes, diarylheptanoic acid and gingerol biosynthetic pathway and flavonoid biosynthetic pathway can be used as allelochemicals in plants.
It is key to the interaction between development and plant environment. In addition, 6 genes and 17 genes belong to stilbene, diarylheptic acid and gingerol biosynthesis pathway and flavonoid biosynthesis pathway, respectively. The phenols and flavonoids produced by biosynthesis of stilbene, diarylheptic acid and gingerol and biosynthesis of flavonoids can be used as allelochemicals in plants.
3.3 Real time qPCR verification results.
To verify the accuracy of transcriptome data set, real-time qPCR was used to analyze the transcription level of 10 randomly selected genes. The relative expression level of EF-1αgene was measured and calculated. These 10 genes are primary amine oxidase activity(c76636.graph_c1-F);shikimate O-hydroxycinnamoyltransferase activity(c64912.graph_c1-F);EIN3-binding F-box protein(c85195.graph_c0-F);CAP(c69462.graph_c0-F);indole-3-acetic acid-amido synthetase(c82761.graph_c0-F);SAUR-like auxin-responsive family protein(c63578.graph_c0-F);flavanone 3-hydroxylase(c70484.graph_c0-F);shikimate hydroxycinnamoyl transferase(c83150.graph_c0-F);chalcone reductase(c85005.graph_c0-F)and chalcone and stilbene synthase family protein(c67970.graph_c0-F).The results of reverse transcription polymerase chain reaction confirmed that the transcriptional changes of these 10 genes were equivalent to the fold changes obtained by our transcriptome analysis (Figure 3, Table S3).
3.4 Metabolomic analysis results.
OPLS-DA analysis was used to measure the metabolic changes to alfalfa caused by biological fertilizer, by examining all metabolites in positive and negative ion mode. With the positive ion mode, R2 =A and Q2 = B, and with negative ion mode, R2 = C and Q2 = D, which indicates the model can describe the sample well and can be employed in the subsequent search for biomarkers. In the PLS-DA scatter plot results, the treatment group and the control group are obviously separated, which indicates that biological fertilizer can affect the normal metabolic pathway of alfalfa(Figure 4a-d).
A total of 1,358 and 539 metabolites were detected in positive and negative ion mode, respectively. Among these, 74 and 28 were detected in piglet plasma in positive and negative ion mode in treatment group and control group respectively (Table S4 and Table S5).
To determine the subordinate pathways of these differential metabolites, KEGG pathway database was used. The main enrichment ways are isoflavone compound synthesis, phenylalanine metabolism, tryptophan and tyrosine synthesis, purine synthesis (Figures 4e-f).
3.5 Comprehensive enrichment analysis of transcripts and metabolite profiles.
Analysis was conducted with KEGG co-enrichment of differential genes and metabolites.Comprehensive enrichment analysis showed that bio-fertilizer could affect isoflavone biosynthesis and metabolism and phenylalanine of alfalfa. Specifically, isoflavone biosynthesis and metabolism, and phenylalanine-related differential genes and metabolites were all up-regulated under the condition of applying biological fertilizer(Figure 5).