Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Plant Sci., 27 January 2026

Sec. Plant Nutrition

Volume 16 - 2025 | https://doi.org/10.3389/fpls.2025.1629610

Physiological and transcriptomic analyses of exogenous calcium in boosting nitrogen use efficiency via oxidative and resistance pathways in peanuts

  • 1Research Institute of Plant Nutrition and Resources & Environments/Henan Key Laboratory of Agricultural Resources and Environment, Henan Academy of Agricultural Sciences, Zhengzhou, Henan, China
  • 2Henan Soil and Fertilizer Station, Department of Agriculture of Henan Province, Zhengzhou, Henan, China

Introduction: Peanuts (Arachis hypogaea L.) exhibit a high demand for calcium, second only to nitrogen and potassium, with calcium playing a critical role in their growth, development, and nitrogen fixation. However, the mechanisms underlying calcium-mediated regulation of peanut growth and nitrogen fixation remain poorly understood.

Methods: In this study, we employed nitrogen-efficient (Puhua 66, Huayu 20) and nitrogen-inefficient (Puhua 28, Shanhua 14) peanut varieties in a two-year field experiment using a split-plot design. The main plots comprised two treatments: standard fertilization (CK) and calcium supplementation (Ca), while the sub-plots consisted of different peanut varieties. We analyzed growth parameters, physiological responses, and transcriptomic profiles.

Results: Our results demonstrated that calcium application significantly increased malondialdehyde (MDA) content in leaves while reducing peroxidase (POD) activity, enhancing pod dry matter accumulation, and promoting earlier plant maturation. Additionally, calcium application elevated the activities of nitrate reductase (NR) and glutamine synthetase (GS) (P < 0.01), thereby improving nitrogen and calcium accumulation in pods, their allocation efficiency, and the overall utilization rates of nitrogen and calcium fertilizers. Transcriptomic analysis revealed 166 differentially expressed genes (DEGs) in nitrogen-efficient varieties and 343 DEGs in nitrogen-inefficient varieties under calcium supplementation, with 67 DEGs shared between the two groups. Functional annotation and qRT-PCR validation were performed on these DEGs.Furthermore, weighted gene co-expression network analysis (WGCNA) indicated that calcium supplementation significantly up-regulated genes associated with sucrose synthase, β-amylase, GTPase-activating proteins, light-harvesting chlorophyll-protein complexes (Lhca2, Lhca3), photosynthetic electron transport (PetF, PetJ), phosphatidylinositol phospholipase C2, inositol-3-phosphate synthase, TMV resistance protein, ABC transporters, ethylene-responsive transcription factors (EIN1, EIN2, EIN3), alkylamine oxidase, glutamate dehydrogenase, and aspartate synthase.

Conclusion: These findings suggest that calcium application modulates carbohydrate metabolism, nitrogen assimilation, plant-pathogen interactions, and photosynthetic processes through differential gene expression, ultimately enhancing leaf physiological activity, dry matter partitioning, pod yield, and early maturation in peanuts.

1 Introduction

In recent years, the large-scale and intensive development of peanut cultivation has made continuous cropping increasingly common. However, long-term continuous cropping disrupts the equilibrium between soil nutrients, microbial communities, and crop requirements (Zhao et al., 2020). As a calcicolous plant, peanut requires up to 2~2.5 kg of calcium per 1000 kg of pod production (Yang et al., 2020). Yet, calcium fertilizer replenishment is frequently overlooked in agricultural practices, exacerbating soil calcium deficits (Yu et al., 2020). Calcium deficiency weakens peanut roots, causes pod shriveling, and increases the proportion of single-kernel pods, ultimately leading to empty pods and significantly reduced yields (Wu et al., 2020). Furthermore, calcium deficiency disrupts peanut metabolic processes, including reduced chlorophyll content, photosynthesis rates, soluble sugars, starch, and leaf soluble protein levels, while impairing antioxidant defense systems (Jain et al., 2011; Singh et al., 2025; Xu et al., 2025).

The application of calcium fertilizer has been shown to significantly improve the distribution ratios of nitrogen, phosphorus, potassium, and calcium in peanut plants, facilitating enhanced nutrient transport to pods and thereby increasing the number of mature pods per plant, pod dry weight, and kernel weight, which collectively contribute to improved yield (Yang et al., 2017; Li et al., 2024). Liu et al. (2021) demonstrated that calcium fertilization enhances photosynthetic rates, prolongs the peak activity of nitrogen metabolism-related enzymes, improves carbon assimilation capacity, and increases kernel yield, albeit with a suppressive effect on vegetative growth. Increasing calcium fertilizer application boosts photosynthetic product accumulation in upland peanuts by elevating the leaf net photosynthetic rate and delaying senescence, thereby significantly raising biomass and economic yield. Additionally, it enhances the activity of nitrogen metabolism enzymes throughout the growth cycle and carbon metabolism enzymes in early stages (Wang et al., 2015). Physiological analyses by Li et al., 2023 and Zuo et al. (2025) further revealed that calcium supplementation elevates antioxidant enzyme activity in leaves, reduces malondialdehyde (MDA) accumulation, and enhances overall physiological performance. Calcium nitrate application mitigates malondialdehyde (MDA) content and electrolyte leakage in peanut leaves while augmenting antioxidant enzyme activity (Huang et al., 2022).

The overuse of chemical fertilizers, reduced application of organic fertilizers, and widespread straw burning have also significantly exacerbated soil salinization and acidification. These processes have led to severe calcium depletion and markedly decreased available calcium content in soils (Abdissa et al., 2018. Adding quicklime to acidic soils not only adjusts soil pH but also enhances fertilizer use efficiency and crop yields (Gao et al., 2016; Liang et al., 2021). Desalegn et al. (2017) demonstrated that quicklime application effectively reduces soil acidity, alleviates aluminum and heavy metal toxicity, improves soil physical properties (e.g., aggregate stability, hydraulic conductivity), and increases soil calcium (Ca²+) contents. Furthermore, quicklime enhances dominant fungal abundance in continuously cropped peanut soils, with a greater impact on fungal community structure than peanut variety selection (Yang et al., 2024).

High-throughput RNA sequencing (RNA-seq) technology has been widely used to decipher the key molecular pathway and candidate genes involving Ca responses in different plant species. Through transcriptomic data indicated that carbohydrates, as important regulatory factors under drought stress, and key differential genes SlGLR 3.2 and SlGLR 3.5 affecting calcium transport in tomatoes under high vapor pressure deficit were identified (Shi and Du, 2023). In high calcium American ginseng, the core differential genes and core metabolites enriched in Glutamate metabolic pathway, compared with suitable calcium American ginseng, the gene (Pg_S3572.10) related with glutathione metabolism in high calcium American ginseng significantly up-regulated by 1.29 times (Zhang et al., 2022). Wang et al. (2022) study showed that DEGs of sweet potato responding to calcium stress were mainly enriched in plant cell wall, symplast and plasmodesma, and one candidate gene related to calcium metabolism was obtained by screening analysis. RNA-seq technology has been pod size (Lv et al., 2024), seed colours (Zhao et al., 2025), and disease (Wang et al., 2024) in peanut, but the mechanism of calcium on nitrogen remains to be systematically elucidated.

Peanut varieties exhibit substantial variation in calcium absorption and utilization efficiency, resulting in differential soil calcium requirements. Previous studies have predominantly focused on the effects of calcium application across different varieties under low-calcium stress, while research on the molecular mechanisms of calcium in nitrogen utilization among different peanut varieties remains scarce. With the completion of peanut genome sequencing, there has been gained a new research opportunity to analyze the characteristics of differentially expressed genes in different peanut varieties under calcium-sufficient and calcium-deficient conditions. This study employs RNA-seq transcriptome analysis to identify differentially expressed genes in different peanut varieties under calcium-sufficient and calcium-deficient conditions, aiming to achieve three main research objectives: (1) to investigate the effects of calcium on physiological characteristics, dry matter accumulation, and nitrogen utilization in different peanut varieties; (2) to identify calcium-responsive and nitrogen utilization-related genes through transcriptome analysis; and (3) to elucidate the new molecular mechanisms by which calcium promotes peanut growth and enhances nitrogen utilization.

2 Materials and methods

2.1 Field site

The two-year field experiment (2020-2021) was conducted in ZhengYang County, Henan Province, China (32°27′25″N, 114°20′28″E). The experimental site featured fluvo-aquic soil with sandy texture, characterized by the following properties: organic matter (11.9 g/kg), total nitrogen (112.3 mg/kg), available phosphorus (18.1 mg/kg), exchangeable calcium (80 mg/kg), available potassium (162.0 mg/kg), available molybdenum (65.0 μg/kg), available boron (120.1 μg/kg), and a pH of 5.2.

2.2 Experiment design

The experiment adopts a split-plot design, with the main plot being two treatment no Ca fertilizer (0) and increased calcium fertilization (CaO, 225 kg/hm2; Liang et al., 2021), while the subplot involves peanut varieties with different nitrogen utilization efficiency types. The low-Nitrogen-efficient peanut variety Puhua 66 and Shanhua 14, and high-Nitrogen-efficient peanut variety Puhua 28 and Huayu 20 were used for the experiments. Basal Ca application was performed 3 d prior to sowing. Using the autumn-sown wheat and summer-sown peanut cultivation mode with ridging, the ridge bottom width is 75 cm, the ridge surface width is 40 cm, and the ridge height is about 18 cm. Each plot contained four rows (row length, 10 m) of peanuts with an equal row spacing of 0.20 m. A randomized design with experimental plots at 3×10.0 m in size and 23,000 plants per ha in density was applied, with three replications of each treatment. Each plot was spaced 0.50 m from each other and used for sampling consisted of four central meters in two central rows. All plots received a banded basal application of 750 kg of compound fertilizer (N: P2O5: K2O=15:15:15) per ha, thereby maintaining the same amount of NPK fertilizer for each treatment. The application of quicklime in the field should be completed one month before peanut sowing. The specific procedure is as follows: uniformly spread the measured quicklime powder across the field, and immediately after spreading, it must be deeply incorporated using tools such as a rotary tiller to ensure thorough mixing with the 0-20 cm plow layer soil. Operators must wear masks, gloves, and protective goggles to prevent lime dust from causing burns to the skin, eyes, and respiratory tract. Standard agronomic management practices were implemented to control diseases, weeds, and insect pests.

2.3 Dry matter weight, nitrogen and calcium content in different organs of peanut

On the 45th day (flowering stage), 70th day (pod-setting stage), and 120th day (harvesting stage) after peanut sowing, five representative plant samples were collected from each plot and dissected into shoots, roots, and pods. Samples were oven-dried at 105 °C for 30 minutes followed by 75 °C until constant weight was achieved. We recorded total dry weight along with organ-specific (shoot, root, pod) dry weights. Dried samples were then ground to pass through an 80-mesh sieve for nutrient analysis. Nitrogen concentration was determined using a Kjeldahl nitrogen analyzer (UPT-K1800N) following H2SO4-H2O2 digestion, while calcium concentration was measured with a flame photometer (M410; Sherwood, UK) after HNO3-HClO4-HF digestion.

The calculation formulas for relevant indicators of N or Ca accumulation and utilization characteristics are as follows:

1. N or Ca accumulation in each organ of the plant (kg/hm2) = N or Ca content in that part of the plant × the dry matter weight of the part.

2. The proportion of N or Ca distribution in different organs (%) = the cumulative amount of N or Ca in the organ/the cumulative amount of nitrogen in the whole plant × 100%.

3. The absorption and translocation amount of N or Ca (mg/plant) = The cumulative amount of N or Ca in vegetative organs at the early podding stage - The cumulative amount of N or Ca in vegetative organs at maturity.

4. The translocation rate of N or Ca (%) = The absorption and translocation amount of N or Ca/The cumulative amount of N or Ca in vegetative organs at the early podding stage × 100%.

5. The contribution rate of N or Ca (%) = The translocation amount of N or Ca/The cumulative amount of N or Ca in pods × 100%.

6. The production efficiency of N or Ca (kg/kg) = Pod weight/Total N or Ca accumulation in the whole plant at harvest time.

7. The harvest index of N or Ca (%) = N or Ca accumulation in pods at harvest time/Total N or Ca accumulation in the whole plant × 100%.

2.4 Determination of enzymatic antioxidant activity and nitrogen assimilation activity in leaves

Two fresh samples (each collected in duplicate) of the second leaf from the main stem apex and root systems from the 0~20 cm topsoil layer of peanut plants at the pod-setting stage were collected and stored at -80 °C. For peroxidase (POD) activity analysis, 0.5 g of finely ground leaf tissue was suspended in 3.0 mL of pre-cooled 50 mM phosphate-buffered saline (PBS) and centrifuged at 12,000 r for 20 min at 4 °C. One enzyme unit (U) was defined as the amount of enzyme required to inhibit the reaction by 50% under assay conditions. POD activity was measured by adding 30 µL of leaf enzyme extract to a 3 mL reaction mixture (0.2 M PBS, 0.38% guaiacol, 0.17% H2O2), with absorbance monitored at 470 nm. Enzyme activity was expressed as Δabsorbance units per minute (U/min).

Malondialdehyde (MDA) content, a marker of lipid peroxidation, was quantified via thiobarbituric acid (TBA) colorimetry as described by Vuković et al. (2022). Frozen leaf powder (0.5 g) was homogenized in 5% (w/v) trichloroacetic acid (TCA) and centrifuged at 12,000 × g for 20 min at 4 °C. The reaction mixture, consisting of 1.5 mL tissue extract and 1.5 mL reagent (0.67% TBA in 5% TCA), was incubated for 30 min at 100 °C. After cooling on ice, absorbance of the chromogenic product was measured at 450 nm, 532 nm, and 600 nm. MDA content was calculated using the formula: MDA content = 6.45 × (A532 − A600) −0.56 × A450.

Nitrate reductase (NR) activity was determined as described by Lea et al. (2006). Leaves (0.5 g) were homogenized in 4 mL of extraction buffer containing 0.1 M HEPES–KOH (pH 7.5), 3% (w/v) polyvinylpyrrolidone (PVP), 1 mM EDTA, and 7 mM cysteine (Cys). The assay mixture comprised 50 mM HEPES–KOH (pH 7.5), 100 mM NADH, 5 mM KNO3, and either 2 mM EDTA or 6 mM MgCl2 in a total volume of 2 mL. NR activity was measured in crude extracts by quantifying nitrite (NO2-) formation after adding 1% (w/v) sulfanilamide and 0.2% (w/v) N-(1-naphthyl) ethylenediamine dihydrochloride in 3 M HCl. The activity state of NR was defined as the percentage of NR activity measured in the presence of Mg²+ (and 14-3-3 proteins) relative to that measured in the presence of EDTA, reflecting the proportion of the enzyme in its nonphosphorylated, active form. All assays were conducted at 25 °C.

Glutamine synthetase (GS) activity was determined in peanut leaves using the hydroximate biosynthetic method described by Pompelli et al. (2010). GS activity was measured with a UV–VIS spectrophotometer (UV-1650PC, Shimadzu, Kyoto, Japan) by quantifying the production of γ-glutamyl-hydroxamate at 540 nm.

2.5 Total RNA extraction, transcriptome sequencing, and analysis

No Ca fertilizer and increased calcium fertilization leaf from 4 peanut cultivar used for were the same as those used for RNA sequencing (RNA-Seq). Three independent biological replicates were taken for each cultivar for a total of 24 samples. RNA-seq library construction and subsequent analyses were performed by Biomarker Technologies (Beijing, China). Total RNA was extracted from peanut leaves during the pod-setting stage and evaluated for concentration, purity, and integrity using agarose gel electrophoresis, a Qubit 2.0 fluorometer, a NanoPhotometer spectrophotometer, and an Agilent 2100 Bioanalyzer. Samples meeting quality standards were retained for downstream experiments. Following RNA purification, library preparation was executed, and the resulting cDNA libraries were sequenced on the Illumina HiSeq 4000 platform. Raw reads generated were subjected to quality control (QC) procedures, followed by adapter trimming using SOAPfuse software to produce clean data. Clean reads were mapped to the reference genome of peanut cultivar Tifrunner (available at: www.Peanutbase.org) using HISAT2 v2.1.0. Transcript assembly was conducted with StringTie (accessed via http://ccb.jhu.edu/software/stringtie/on September 15, 2022). Gene expression levels were quantified using FPKM (Fragments Per Kilobase of transcript per Million mapped reads), calculated as:FPKM = (Number of cDNA fragments)/(Total mapped reads in millions) × (Transcript length in kilobases). Differentially expressed genes (DEGs) were identified based on a false discovery rate (FDR) < 0.001 and |log2(Fold Change)| ≥ 1. Functional annotation of DEGs was performed using KOBAS v2.1.1, incorporating Gene Ontology (GO; http://geneontology.org/) and Kyoto Encyclopedia of Genes and Genomes (KEGG; https://www.kegg.jp/) databases. Additionally, DEG sequences were aligned and annotated against the Nr, Swiss-Prot, GO, and KEGG databases using BLAST software to obtain comprehensive annotation information. Principal component analysis (PCA), Pearson correlation coefficient (PCC) calculations, and weighted gene co-expression network analysis (WGCNA) were executed on the BMKCloud platform (https://international.biocloud.net) to further characterize the transcriptome data.

2.6 Real-time quantitative PCR validation

RNA extraction was performed using TRIzol reagent (Invitrogen, Waltham, USA). Quantitative fluorescence primers were designed using Primer 5.0 software (Supplementary Table S1). Total RNA was used as the template for cDNA synthesis, which was performed using reverse transcriptase (Takara Bio, Dalian, China) according to the manufacturer’s instructions. For cDNA synthesis and quantitative real-time PCR (qRT-PCR), a previously published method was employed (Xu et al., 2022). The quality of the reverse-transcribed cDNA was verified using primers for the peanut internal control gene, actin.

2.7 Statistical analysis

The statistical analysis was performed using the DPS 7.5 software with Duncan’s test (p < 0.05). Using Microsoft Excel 2007 and Origin 8.5 software to organize data and create graphs.

3 Results and analysis

3.1 The impact of calcium fertilizer on the enzyme activity of peanut leaves

Peroxidase and malondialdehyde are widely present in plants and play regulatory roles in plant development. As shown in Figure 1, the POD activity in peanut leaves decreased after the application of additional calcium fertilizer across all varieties, with significant reductions observed in Puhua 66, Puhua 28, and Shanhua 14. Conversely, MDA levels increased in the leaves of all varieties following calcium fertilizer application, with significant increases in Puhua 66 and Shanhua 14. Moreover, the application of additional calcium fertilizer resulted in more concentrated flowering across different peanut varieties, a shortened duration of functional leaves, and an accelerated completion of the entire life cycle (Figures 1E, F). These changes suggest that increased calcium fertilizer application enhances MDA activity, reduces POD activity, accelerates plant growth and development, and promotes early maturity in peanut plants. Nitrate reductase and glutamine synthetase are key enzymes in ammonium assimilation and nitrogen metabolism. As shown in Figure 1, the activities of both NR and GS in peanut leaves increased after the application of additional calcium fertilizer. Notably, GS activity reached significant levels in all varieties. This finding is consistent with T1, indicating that calcium fertilizer application promotes GS activity in the later stages of peanut growth. The T2 treatment exhibited higher GS activity across all four growth stages compared to other treatments, suggesting that increasing calcium fertilizer application while ensuring adequate nitrogen fertilizer enhances GS activity in peanuts.

Figure 1
Panel A shows a gene network with nodes labeled with gene identifiers, represented as circles of varying sizes and colors. Panel B is a heatmap with genes on the y-axis and different conditions on the x-axis, colored from green to red, indicating varying levels of expression.

Figure 1. Effect of calcium fertilizer on enzyme activity of leaves and senescence performance of pod and maturity stage for peanut. (A) POD activity of peanut leaves; (B) MDA activity of peanut leaves; (C) NR activity of peanut leaves; (D) GS activity of peanut leaves; (E) Peanut senescence performance of pod stage; (F) Peanut senescence performance of maturity stage. Different capital letters indicate significant difference among treatments (P < 0.05). The same as blow.

3.2 The impact of calcium fertilizer on dry matter accumulation of peanuts

Calcium fertilizer treatments increased dry matter accumulation in peanut nutritional organs (shoots and roots) by 3.8% to 39.9% compared to the control (CK) without calcium fertilizer, while pod dry matter accumulation was significantly increased in all four peanut varieties (Figure 2). The accumulation trends of shoots, roots, and total dry matter under different calcium fertilization treatments were similar. During the seedling stage, calcium-treated plants exhibited lower total, shoot, and root dry matter accumulation than CK, particularly for shoot and total dry matter accumulation in Puhua 28 and Huayu 20, which reached significant levels. During the pod-setting stage, total and shoot dry matter accumulation were lower than CK, while root and pod dry matter accumulation were higher than CK. Specifically, pod dry matter accumulation increased by 26.9%, 10.6%, 67.2%, and 30.2% in Puhua 28, Puhua 66, Shanhua 14, and Huayu 20, respectively. At maturity, compared to CK, total dry matter accumulation increased by 1.1%, 6.5%, 9.4%, and 1.3% in Puhua 28, Puhua 66, Shanhua 14, and Huayu 20, respectively, while pod dry matter accumulation increased by 19.7%, 20.9%, 62.9%, and 48.3%, respectively. However, shoot dry matter accumulation decreased by 11.2%, 3.79%, 25.0%, and 39.0%, respectively. These results indicate that although calcium fertilizer reduces dry matter accumulation in the early growth stage, it promotes the translocation of dry matter from shoots to pods in the middle and late stages, ultimately increasing peanut yield (Figure 2).

Figure 2
A composite image consisting of four panels: A shows a 3D plot of principal component analysis with different markers and colors representing various samples. B is a bar graph comparing upregulated and downregulated DEG gene counts across four samples: PH66, PH28, SH14, and HY20. C features a Venn diagram illustrating the overlap of different gene sets among four samples: SH14, HY20, PH28, and PH66. D presents a bubble chart depicting the statistics of pathway enrichment, listing various metabolic and biosynthesis pathways, with bubble size indicating count and color representing q-value.

Figure 2. Effect of calcium fertilizer on different organs dry matter in different stage peanut. (A) Puhua 28 different organs dry matter; (B) Shanhua 14 different organs dry matter; (C) Puhua 66 different organs dry matter; (D) Huayu 20 different organs dry matter.

3.3 The impact on peanut nitrogen absorption and utilization

Significant differences in nitrogen accumulation across different peanut organs are attributed to variety, additional calcium fertilizer application, and their interaction (Table 1). For the same variety, nitrogen accumulation in shoot dry matter is significantly reduced under additional calcium fertilizer, while pod nitrogen accumulation significantly increases, with minimal impact on roots. Among the same treatments, high-nitrogen varieties, such as Puhua 28 and Huayu 20, exhibit significantly higher pod nitrogen accumulation than low-nitrogen varieties, such as Puhua 66 and Shanhua 14. Compared to the control, calcium fertilizer treatment increases pod nitrogen accumulation by 89.13 g/kg, 143.11 g/kg, 307.03 g/kg, and 420.53 g/kg in Puhua 28, Puhua 66, Shanhua 14, and Huayu 20, respectively.

Table 1
www.frontiersin.org

Table 1. Effects of different calcium treatments on nitrogen accumulation and distribution in peanut.

Distribution rates increase by 7.7%, 5.79%, 19.61%, and 15.84%, respectively, while nitrogen translocation rates rise by 4.73%, 5.08%, 4.91%, and 7.58%, respectively. Under the interaction of variety and calcium fertilizer application, Huayu 20 with additional calcium fertilizer achieves the highest total pod nitrogen accumulation, with a nitrogen distribution ratio of 82.82%. Calcium fertilizer application promotes nitrogen allocation to reproductive organs thereby increasing the nitrogen harvest index and nitrogen use efficiency. Compared to no calcium fertilizer application, the pod harvest index of Puhua 28, Puhua 66, Shanhua 14, and Huayu 20 significantly increases by 0.07, 0.06, 0.20, and 0.16, respectively, while nitrogen use efficiency increases by 3.58 kg/kg, 2.34 kg/kg, 6.49 kg/kg, and 3.48 kg/kg under calcium fertilizer treatments.

3.4 Absorption and utilization of calcium in peanuts

As shown in Table 2, calcium accumulation in the roots is significantly reduced under additional calcium fertilizer, while pod calcium accumulation significantly increases. The overall calcium distribution across various organs follows the order: pods > shoots > roots. Compared to no calcium fertilizer application, calcium accumulation in the pods of Puhua 28, Puhua 66, Shanhua 14, and Huayu 20 significantly were increased by 602.28 kg/hm², 351.00 kg/hm², 1186.38 kg/hm², and 1041.48 kg/hm², respectively. Total plant calcium accumulation significantly increased by 475.02 kg/hm², 415.44 kg/hm², 612.54 kg/hm², and 848.16 kg/hm², respectively, while root calcium accumulation decreased by 149.04 kg/hm², 160.56 kg/hm², 58.50 kg/hm², and 52.92 kg/hm² under calcium fertilizer treatments. Calcium fertilizer application not only was increases calcium accumulation in various organs but also significantly enhances calcium translocation rate, production efficiency, and harvest index. The calcium translocation rate of Puhua 28, Puhua 66, Shanhua 14, and Huayu 20 were increased by 6.97%, 17.64%, 31.56%, and 13.85%, respectively. Among different varieties, Huayu 20 had the highest calcium harvest index, followed by Shanhua 14, then Puhua 28, with Puhua 66 having the lowest. The improvement in calcium production efficiency ranged from 8.03 to 34.36 kg/kg, with greater increases observed in nitrogen-inefficient varieties (Puhua 66 and Shanhua 14) compared to nitrogen-efficient varieties (Puhua 28 and Huayu 20), reaching significant levels.

Table 2
www.frontiersin.org

Table 2. Effects of different calcium treatments on calcium accumulation and distribution in peanut.

3.5 Calcium fertilizer application and transcriptome analysis in peanut leaves

To identify genes potentially involved in the response to calcium fertilizer, RNA-seq was performed on leaves from the pod-setting stage of different peanut varieties treated with additional calcium fertilizer compared to the control. Principal component analysis (PCA) and Pearson correlation coefficients revealed good repeatability among biological replicates, although the separation of the four varieties was not distinct, with correlation coefficients exceeding 0.9 (Figures 3a). In total, we obtained 21.17 million genes/transcripts and 137.5 Gb of clean reads from the 24 samples, averaging 5.7 Gb per sample. The Q30 values were 91.53% and 96.82%, with 96.82% of reads having quality scores at the Q20 level. Alignment efficiency ranged from 83.15% to 93.70%, yielding a total of 87,623 genes or transcripts. Among these, 66,667 genes were expressed in all four varieties (FPKM ≥ 0, FDR ≤ 0.001; Figure 3c). Pairwise comparisons identified 4,041 (Puhua 28), 7,756 (Puhua 66), 5,845 (Shanhua 14), and 3,944 (Huayu 20) differentially expressed genes (DEGs), with 3,086 overlapping DEGs. The criteria used to judge significant differences in gene expression (fold change > 2, FPKM ≥ 0.1, FDR ≤ 0.001) were stringent (Figures 3b). There were 166 genes specifically expressed in nitrogen-efficient varieties and 343 genes specific to nitrogen-inefficient varieties, with 67 overlapping DEGs (43 up-regulated and 24 down-regulated) (Figure 3c, Supplementary Table S3). And the 67 DEGS significant enrichment in KEGG pathways related to dry matter accumulation, nitrogen metabolism, early maturity (Figure 3d).

Figure 3
Panel A displays a heat map showing gene expression levels, with red and green indicating high and low expression, respectively, across multiple samples labeled on the bottom. Panel B contains six bar graphs, each illustrating the relative expression of different genes (EXL3, APS1, At3g50280, BHLH130, MES1, At1g66480) across four conditions (HY20, PH66, PH28, SH14) using RNA and qT-PCR data, indicated by different colored bars.

Figure 3. Transcriptome data of peanut accessions with calcium fertilizer application. (A) PCA analysis; (B) number of transcripts per varieties; (C) venn diagram number of up-regulated DEGs in different varieties; (D) DEG-enriched KEGG annotations.

3.6 Identification of WGCNA modules associated with calcium fertilizer in peanut leaves

To further investigate key genes responsive to calcium fertilizer, we conducted weighted gene co-expression network analysis (WGCNA). A scale-free topological model was established using a soft threshold of 14. Using the dynamic hierarchical tree cutting algorithm, we identified 14 color modules (Supplementary Figure S3a). Among these, the purple module had the highest gene count (779 genes), while the white module had the lowest (34 genes) (Supplementary Figure S3b). Module-trait relationship analysis revealed that the black module exhibited significant correlations: a negative correlation with MDA and GS levels, and a positive correlation with POD activity (Supplementary Figure S3d). We then focused on the black module for functional enrichment analysis using GO and KEGG analyses. The results showed significant enrichment in KEGG pathways related to plant hormone signal transduction, ABC transporters, spliceosome, amino acid metabolism, starch and sucrose metabolism, and cell formation (Supplementary Figure S4).

A co-expression network was constructed to identify key genes involved in the response to calcium fertilizer. The top 10 hub genes with superior MCC scores were identified, and a visual gene interaction network diagram was generated for these genes and their associated genes (with a weight value greater than 0.20) (Figure 4). Functional annotation revealed that these 10 genes encode proteins involved in various biological processes: nucleic acid binding (AT-hook motif nuclear-localized protein 23, Remorin, C-terminal region), nitrogen metabolism (Stilbene synthase 3), growth and development (NDR1/HIN1-like protein 1, Chloroplastic import inner membrane translocase subunit), and other functions (Berberine bridge enzyme-like, NADH-ubiquinone reductase complex 1) (Supplementary Table S5). Under increased calcium fertilizer application, these core genes were down-regulated in the low-efficiency variety but up-regulated in the nitrogen-efficient variety.

Figure 4
Heatmap and bar graphs illustrating gene expression data. Panel A shows a heatmap with red and green representing different expression levels for various genes across conditions. Panel B contains multiple bar graphs comparing relative expression levels between RNA and qRT-PCR for specific genes like EXL3, APS1, and MES1 across different treatments.

Figure 4. Hub genes based WGCNA. (A) Co-expression networks of the top 10 hub genes in the black; (B) Heatmap of the expression patterns of Hub genes.

3.7 Calcium improved expression of dry matter accumulation-related genes

This study focused on the DEGs in the black module, as these genes play a crucial role in responding to calcium fertilizer application in peanut (Figure 5a). Through GO and KEGG analyses of the black module, terms related to dry matter accumulation, chlorophyll metabolism, and sugar metabolism were found to be enriched. A total of 31 DEGs associated with dry matter accumulation were identified. Under calcium fertilizer application, genes involved in calcium signaling (calmodulin, HPT1), photosynthesis (PSAF, ISPF, NTH2, PSY, PCMP-H51), and sugar synthesis (AMY1.1) were up-regulated, while genes involved in lipid metabolism (At5g55050), root phototropism (RPT2), and transport (ABCC3, GAT1, ABCG14, UGT83A1) were down-regulated. Notably, ABCC3 was up-regulated in the nitrogen-efficient variety but down-regulated in the nitrogen-inefficient variety. Quantitative PCR (qPCR) results confirmed that the expression patterns of ABCC3, UGT83A1, PSAF, AH13G38220, PSY, and At3g50280 were consistent with those obtained from RNA-Seq data (Figure 5b), validating the reliability of the RNA-Seq results.

Figure 5
Bar charts and photos show dry matter distribution in different plant types: Puhua 28 (A), Shanshua 14 (B), Puhua 66 (C), and Huayu 20 (D). Each bar is divided into root, shoot, and pod segments. Dry matter is measured in grams per plant. Columns are labeled CK and Ca indicating treatment. Photographs next to each chart depict the corresponding plants.

Figure 5. Expression differences of DEGs related to dry matter accumulation in black module based on WGCNA. (A) heatmap of the expression patterns of DEGs, FPKM values were normalized by Z-score; (B) qRT-PCR validation of key dry matter accumulation genes.

3.8 Calcium improved expression of nitrogen metabolism-related genes

The black module of WGCNA identified 45 DEGs specifically associated with nitrogen metabolism and amino acid metabolism (Figure 6a). Under calcium fertilizer application, genes involved in amino acid binding (CRK29, AS2, WCRKC2), plant-pathogen interactions (BHLH30, BHLH106, BHLH148, BHLH100, BHLH130), and heat shock responses (HSFA2C, HSP26-A) were up-regulated, while genes involved in nitrogen metabolism (NHL3), domain-specific functions (LBD40), and magnesium ion transport (MRS2-I) were down-regulated. Additionally, EXL3, CBP60B, MES1, APS1, BHLH126, and At1g66480 were up-regulated in nitrogen-efficient varieties but down-regulated in nitrogen-inefficient varieties. qRT-PCR analysis further confirmed that the expression trends of EXL3, CBP60B, MES1, APS1, BHLH130, and At1g66480 were consistent with the RNA-seq results (Figure 6b), thereby demonstrating the reliability of the RNA-seq data.

Figure 6
Four bar graphs labeled A, B, C, and D show POD, MDA, NR, and GS levels across four peanut varieties (Puhua 66, Puhua 68, Shanhua 14, Huayu 20) under CK and Ca treatments. CK is in blue, and Ca is in orange, with statistically significant differences marked by letters. Two photos labeled E and F depict field comparisons of peanut plants under CK and Ca treatments. E shows greener plants, while F displays less healthy plants.

Figure 6. Expression differences of DEGs related to nitrogen metabolism in black module based on WGCNA. (A) heatmap of the expression patterns of DEGs, FPKM values were normalized by Z-score; (B) qRT-PCR validation of key nitrogen metabolism genes.

3.9 Calcium improved expression of early maturity -related genes

The black module from WGCNA identified 43 DEGs significantly enriched in terms related to programmed cell death, plant-pathogen interactions, kinases, and ethylene transcriptional synthesis (Figure 7a). Under increased calcium fertilization, genes related to ethylene (ERF1B, ERF115, ERF110), senescence (MYB58, MYB86, MYB330), plant circadian rhythm (COL9), and peroxidase activity (CSA) were up-regulated, while genes related to senescence (ERF020, ERD7), protein kinase activity (RLP6), and ethylene (ERF020, ERF106) were down-regulated. Additionally, the BTB/POZ domain protein (At5g03250), peroxidase (PER20), ethylene-related genes (ERF1B, ERF011), plant circadian rhythm gene (At3g19360), senescence-related genes (MYB4, MYB3R5), and receptor-like protein kinase (RLP6) were up-regulated in nitrogen-efficient varieties but down-regulated in nitrogen-inefficient varieties. qRT-PCR analysis confirmed that the expression trends of FLZ6, MAKR5, ERF1B, MYB3R5, SAG20, and ASR2 were consistent with the RNA-seq results (Figure 7b), thereby demonstrating the reliability of the RNA-seq data.

Figure 7
Illustration showing the effects of calcium (Ca) on peanut plants. Red arrows connect Ca to various gene expressions and biochemical changes: ABCC3, UGT83A1, PSAF, and others. Biochemical outcomes include increased GS, decreased POD, and increased NR, promoting early plant maturation. Images below show comparative growth between control (CK) and Ca-treated plants, highlighting improved nitrogen utilization and dry matter accumulation.

Figure 7. Expression differences of DEGs related to early maturity in black module based on WGCNA. (A) heatmap of the expression patterns of DEGs, FPKM values were normalized by Z-score; (B) qRT-PCR validation of key early maturity genes.

In total, 14 genes related to plant precocity were identified, enriched in 16 pathways, including MAPK signaling pathway-plant (ko04016), plant hormone signal transduction (ko04075), plant-pathogen interaction (ko04626), biosynthesis of amino acids (ko01230), glycolysis/gluconeogenesis (ko00010), and protein processing in the endoplasmic reticulum (ko00030). Based on functional annotations, five of the DEGs encode ethylene-responsive transcription factors, while one DEG encodes the transcription factor MYB3R-5, which is involved in plant hormone signal transduction. Additionally, five other differentially expressed genes encode probable disease resistance proteins, peroxidase 20, abscisic stress-ripening protein, class I heat shock protein, and serine/threonine-protein phosphatase 7 long form homolog (Supplementary Table S6).

4 Discussion

4.1 Calcium improved dry matter accumulation and distribution of peanut pods

Dry matter accumulation is a key indicator of crop population quality. Typically, the distribution of dry matter among different plant parts determines peanut yield levels (Yang et al., 2020). After pod formation, the proportion of dry matter in leaves, stems, and petioles begins to decline, while the proportion in pods increases as they become the primary growth center. By maturity, the proportion of dry matter in pods reaches its peak (Wu et al., 2020). Soil acidification and saline-alkali stress can significantly inhibit dry matter accumulation in peanuts, while calcium deficiency severely impairs pod development. Calcium deficiency in the soil hinders the transport of photosynthetic products to pods, leading to poor pod development or abortion (Bi et al., 2024). Moreover, when photosynthetic products are retained in shoot organs, more resources are diverted to vegetative growth, potentially disrupting the balance between vegetative and reproductive growth and resulting in excessive vegetative growth. Calcium fertilizer application has no significant effect on the weight of vegetative organs but significantly (p < 0.0001) increases pod biomass and harvest index in red soil dryland (Wu et al., 2020). This study demonstrates that calcium application significantly enhances pod and whole-plant biomass. Additionally, nitrogen-efficient varieties exhibit greater increases in various indicators compared to nitrogen-inefficient varieties, with reproductive organs (pods) often showing higher increases than vegetative organs (stems, leaves, and roots). Calcium application also reduces the V/R (vegetative organs/reproductive organs) ratio, consistent with findings reported by Yang et al. (2017).

Further transcriptome sequencing analysis identified 31 genes associated with dry matter accumulation. These genes showed enrichment in pathways such as Photosynthesis (ko00195), carbon metabolism (ko01200), Starch and sucrose metabolism (ko00500), Pentose and glucuronate interconversions (ko00040), starch and sucrose metabolism (ko00500), and ABC transporters (ko02010) (Supplementary Table S6). Based on functional annotations, these differentially expressed genes encode TMV resistance proteins, ABC transporters, alkylamine oxidase, glutamate dehydrogenase, and aspartate synthase. Among them, AH09G22950 and AH19G28970 encode ethylene-responsive transcription factor 106 (ERF106), a member of the largest subfamily of ethylene-responsive factors (ERFs). Previous studies have shown that ERFs are key regulatory factors in plant development and stress responses (Xie et al., 2019). Specifically, the APETALA2/ERF subfamily is implicated in growth and developmental processes mediated by plant hormones such as gibberellins, cytokinins, and brassinosteroids in Arabidopsis. Another gene, AH17G00840, encodes an ABC transporter family protein. ABC transporters (ABCC and ABCG) have been shown to facilitate the transport of secondary metabolites, such as capsaicin and dihydrocapsaicin, into ovule vacuoles during pepper fruit ripening, thereby affecting their content and promoting fruit development (Lopez-Ortiz et al., 2019). Quantitative analysis indicated that genes associated with dry matter accumulation were up regulated under calcium fertilizer application. This suggests that these genes play a major role in regulating ABC transport and sugar metabolism pathways, thereby promoting dry matter accumulation and translocation towards the pods.

4.2 Calcium application on the accumulation and utilization of nitrogen and calcium in peanuts

Previous studies have shown that calcium nitrate application increases nitrogen accumulation and promotes its transport from vegetative organs to “sink” tissues in peanuts (Sadras et al., 2019). As calcium application increases, nitrogen and phosphorus accumulation and distribution rates in pods rise, while these rates in leaves, stems, roots, pegs, and shells significantly decrease (Bedoussac et al., 2014). Calcium application may positively affect nutrient absorption and accumulation in peanut organs under salt stress (Li et al., 2024). In this study, calcium fertilizer application significantly increased nitrogen content and distribution rates in roots and pods (underground organs) compared to the control without calcium treatment, while nitrogen content in stems and leaves (aboveground organs) and nitrogen distribution rates in leaves were significantly lower. These findings are consistent with changes in dry matter and nitrogen accumulation in peanuts, indicating that calcium application enhances nitrate reductase and glutamine synthetase activities, thereby increasing nitrogen absorption and utilization and providing a foundation for nitrogen accumulation. Nitrogen-efficient peanut varieties (Puhua 28 and Huayu 20) exhibited significantly higher nitrogen accumulation and distribution rates in pods compared to nitrogen-inefficient varieties. Transcriptome analysis revealed that 45 DEGs were enriched in three nitrogen metabolism pathways (Supplementary Table S6). Four of these genes were significantly enriched in the amino acid biosynthesis pathway (ko01230), with two genes (AH09G22950, AH19G28970) involved in converting fructose-6P to glutamate-P and the other two genes (AH01G04690, AH11G03850) involved in aspartic acid synthesis. Four genes were significantly enriched in the alanine, aspartate, and glutamate metabolism pathway (ko00250), with two genes (AH01G04690, AH11G03850) involved in L-asparagine synthesis and the third gene (AH10G29320, AH03G20960) involved in glutamate synthesis from NH3. These results suggest that calcium fertilizer application enhances NR and GS activities in leaves and directs more nitrogen to reproductive organs, preventing excessive vegetative growth and laying a foundation for improved plant population quality.

In plants, calcium primarily moves from the xylem to leaves via transpiration, which is the main driving force for calcium transport (Rahman et al., 2015). Calcium uptake in leaves under high calcium supply affects the absorption of other cations and influences stomatal morphology and photosynthesis (Johnson et al., 2022; Ma et al., 2022). Stomatal conductance, size, and density influence CO2 absorption, while CO2 diffusion and fixation depend on mesophyll tissue and chloroplast distribution (Tomeo and Rosenthal, 2017). CO2 diffusion involves passage through cell walls, plasma membranes, cytosol, chloroplast membranes, intercellular spaces, chloroplast matrix, and cell membranes in mesophyll cells (Du et al., 2019). Fruit hardness is affected by pectin and total calcium concentrations, which in turn affect cell wall structural stability (Althiab-Almasaud et al., 2021). Transcriptome analysis showed that major metabolic pathways involved in pod formation include phenylpropanoid metabolism, followed by starch, sugar, and glutathione metabolism, with lignin synthesis being a key factor. AH05G03640 encodes a peroxidase involved in plant stimulus metabolism, transport, and catabolism, and participates in the synthesis of various lignin types in the phenylpropanoid biosynthesis pathway (ko00940) (Supplementary Figure S5). NewGene_16049 is a novel gene involved in coenzyme transport and metabolism in the RNA transport pathway (ko03013) (Supplementary Figure S5). Calcium application increases calcium content in stems, roots, pegs, pods, and kernels, and raises calcium distribution rates in peanut leaves while decreasing rates in stems. These changes in photosynthetic product transport capacity, cell wall structural stability, carbohydrate and amino acid levels, and differentially expressed genes effectively improve calcium distribution and single pod weight under calcium application.

4.3 The effect of calcium application on the early maturity of peanuts

Plants rely on leaves to provide energy and exchange materials with the environment. As the growth period progresses, the active oxygen scavenging mechanism in peanut functional leaves becomes imbalanced. Calcium application enhances this mechanism, increasing chlorophyll content and the activities of SOD, POD, and catalase enzymes in peanut leaves, while decreasing MDA content. This helps eliminate or avoid damage caused by the accumulation of oxygen free radicals in leaves, which is crucial for normal plant growth (Jawad et al., 2020). MDA, a final product of lipid peroxidation, reflects the degree of reactive oxygen species (ROS) accumulation and damage to cell membranes under environmental stress (Huang et al., 2022). In recent years, the role of calcium has been widely discussed. Calcium ions (Ca²+) can directly or indirectly mediate various cellular signaling processes through calmodulin (CaM), modulating ion channel activities, gene expression, and oxidative stress tolerance. Calcium deficiency can disrupt the antioxidant protection system in peanuts, reducing CAT and POD activities, increasing MDA content, relative electrical conductivity, and superoxide radical (O2-) production rates, leading to increased cell membrane permeability and accelerated peanut senescence (Singh et al., 2025; Xu et al., 2025). Calcium application can effectively reduce peanut flowering time, promoting early flowering and facilitating the transition from vegetative to reproductive growth.

This study found that calcium application increased SOD, POD, and CAT activities and decreased MDA content in leaves of most peanut varieties, reducing cell membrane damage and cell sap exudation. This delays leaf aging, improves photosynthesis, and stabilizes cell structure and function, thereby increasing plant yields. Additionally, nitrogen-efficient peanut varieties (Huayu 20 and Puhua 66) showed more significant effects compared to nitrogen-inefficient varieties (Puhua 28 and Shanhua 14) under calcium fertilizer application. Transcriptome analysis annotated 45 DEGs with GO terms related to oxygen species metabolic processes. Five DEGs encoded ethylene-insensitive (EIN) proteins involved in pathway ko04626 (Supplementary Table S6). EIN is a major regulatory factor in the ethylene signaling pathway. For example, enhanced expression of EIN3 promotes senescence and regulates ethylene-induced chlorophyll degradation in Arabidopsis thaliana leaves (Wang et al., 2020). EIN3 also activates leaf senescence under adverse environmental conditions in Sorghum bicolor. Besides its role in leaf senescence, EIN3 family genes are increasingly recognized for their involvement in defense responses to pathogens and environmental stimuli (Xu et al., 2023). For instance, overexpression of EIN3 in Arabidopsis reduces resistance to Pseudomonas syringae, while suppressing this gene or in EIL1 mutants enhances resistance, indicating a negative regulatory role in disease resistance (Wu et al., 2016). Similarly, TaEIL1, a homologue of AtEIN3 from wheat, negatively regulates defense against wheat stripe rust (Duan et al., 2013). In this study, we observed up-regulation of EIN expression in peanut leaves, potentially leading to down-regulation of pathogen defense response transcripts. We found that calcium application in peanut fields led to preferential leaf shedding, earlier fruit development, and a shortened growth period, aligning with our hypothesis.

5 Conclusions

In summary, we sequenced the transcriptomes of peanut varieties treated with exogenous calcium. The DEGs were significantly enriched in carbohydrate metabolic processes, nitrogen metabolic processes, oxidoreduction processes, and photosynthetic processes. These DEGs were up-regulated under exogenous calcium treatment, leading to increased levels of GS, MDA, and NR in leaves. This enhancement significantly improved nitrogen use efficiency and the pod allocation ratio, as well as promoted plant maturity and yellowing, resulting in a shorter growth growth period for peanuts (Figure 8). This study investigates the effects of calcium application on the molecular and physiological aspects of peanut growth and development, providing scientific evidence that calcium fertilizer application can increase peanut yield in acidic or calcium-deficient soils.

Figure 8
Panel A shows a heatmap with genes listed on the right and conditions at the bottom. Color intensity from red to green indicates expression levels. Panel B presents bar graphs for six genes, illustrating relative expression under different conditions, with blue representing RNA and orange showing qT-PCR results.

Figure 8. A molecular regulatory model for calcium application in peanut.

Data availability statement

The data presented in the study are deposited in the NCBI, accession number PRJNA1314981.

Author contributions

FX: Formal Analysis, Writing – original draft, Writing – review & editing. FX: Formal Analysis, Funding acquisition, Writing – original draft. XZ: Formal Analysis, Funding acquisition, Investigation, Writing – original draft. LL: Funding acquisition, Date curation, Formal analysis, Project administration, Writing – original draft, Writing – review & editing. XS: Writing – original draft, Writing – review & editing, Validation, Formal analysis, Methodology. YS: Writing – original draft, Writing – review & editing, Data curation, Formal analysis, Project administration. XW: Writing – original draft, Writing – review & editing, Investigation, Resources. ZZ: Software, Investigation, Writing – original draft. QL: Formal analysis, Investigation, Writing – original draft.

Funding

The author(s) declared that financial support was received for work and/or its publication. This work was funded by the Excellent Youth Science and Technology Fund of Henan Academy of Agricultural Sciences (2024YQ18); Henan Provincial Science and Technology Project (252102110146); Henan Provincial Youth Science Fund Project (252300423617); the China Agriculture Research System (CARS-13); the Henan Agriculture Research System (HARS-22-05-G2); Innovative team Science and Technology Fund of Henan Academy of Agricultural Sciences (2024TD10).

Conflict of interest

The authors declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.2025.1629610/full#supplementary-material

Supplementary Table 1 | Primer sequences used for quantitative real-time RT-PCR (qRT-PCR).

Supplementary Table 2 | RNA-Seq reads quality control result parameters.

Supplementary Table 3 | Expressions and annotation of DEGs in four comparable comple.

Supplementary Table 4 | KEGG pathways significantly enriched by 4-coexpression DEGs.

Supplementary Table 5 | Hub genes annotation.

Supplementary Figure 1 | GO and KEGG enrichment specifically in low nitrogen efficiency varieties.

Supplementary Figure 2 | GO and KEGG enrichment specifically in high nitrogen efficiency varieties.

Supplementary Figure 3 | WGCNA process.

Supplementary Figure 4 | KEGG and GO enrichment black module by WGCNA.

Supplementary Figure 5 | The key pathways Ca for peanut grow.

References

Abdissa, B., Kibebew, K., Bobe, B., Markku, Y. H., and Tesfaye, B. (2018). Effectss of lime, vermicompost, and chemical P fertilizer on selected properties of acid soils of ebantu district, western highlands of Ethiopia. Appl. Environ. Soil Sc. 2018, 1–13. doi: 10.1155/2018/8178305

Crossref Full Text | Google Scholar

Althiab-Almasaud, R., Chen, Y., Maza, E., Djari, A., Frasse, P., Mollet, J., et al. (2021). Ethylene signaling modulates tomato pollen tube growth through modifications of cell wall remodeling and calcium gradient. Plant J. 107, 893–908. doi: 10.1111/tpj.15353

PubMed Abstract | Crossref Full Text | Google Scholar

Bedoussac, L., Etienne-Pascal, J., Hauggaard-Nielsen, H., Naudin, C., Corre-Hellou, G., Jensen, E. S., et al. (2014). Ecological principles underlying the increase of productivity achieved by cereal-grain legume intercrops in organic farming. A review. Agron. Sustain Dev. 35, 911–935. doi: 10.1007/s13593-014-0277-7

Crossref Full Text | Google Scholar

Bi, W., Yin, Y., Ding, C., Tu, X., Zhou, Z., and Wang, X. (2024). Insights into the antagonistic effects of calcium on cadmium accumulation in peanuts (Arachis hypogaea L.). J. Environ. Manage 367, 122003. doi: 10.1016/j.jenvman.2024.122003

PubMed Abstract | Crossref Full Text | Google Scholar

Desalegn, T., Alemu, G., Adella, A., Debele, T., and Julin Gonzalo, J. (2017). Effects of lime and phosphorus fertilizer on acid soils and barley (hordeum vulgare L.) performance in the central highlands of Ethiopia. Exp. Agr 53, 432–444. doi: 10.1017/S0014479716000491

Crossref Full Text | Google Scholar

Du, Q., Liu, T., Jiao, X., Song, X., Zhang, J., and Li, J. (2019). Leaf anatomical adaptations have central roles in photosynthetic acclimation to humidity. J. Exp. Bot. 70, 4949–4961. doi: 10.1093/jxb/erz238

PubMed Abstract | Crossref Full Text | Google Scholar

Duan, X., Wang, X., Fu, Y., Tang, C., Li, X., Cheng, Y., et al. (2013). TaEIL1, a wheat homologue of AtEIN3, acts as a negative regulator in the wheat-stripe rust fungus interaction. Mol. Plant Pathol. 14, 728–739. doi: 10.1111/mpp.12044

PubMed Abstract | Crossref Full Text | Google Scholar

Gao, L., Wang, S., Liu, Z., Hilman, Wang, H., and Huang, J. (2016). Two cultivars of peanut (arachis hypogaea) seedlings show different tolerance to calcium deficiency. J. Plant Nutr. 39, 1016–1025. doi: 10.1080/01904167.2015.1109107

Crossref Full Text | Google Scholar

Huang, Y., Chen, J., Sun, Y., Wang, H., Zhan, J., Huang, Y., et al. (2022). Mechanisms of calcium sulfate in alleviating cadmium toxicity and accumulation in pak choi seedlings. Sci. Total Environ. 805, 150115. doi: 10.1016/j.scitotenv.2021.150115

PubMed Abstract | Crossref Full Text | Google Scholar

Jain, R. K., Singh, G. B., Singh, A. P., Goel, R. K., Aryya, N. C., and Jha, S. K. (2011). Role of measurement of antioxidant enzymes in evaluation of antioxidant therapy in tobacco abusers with oral leukoplakia. Indian J. Otolaryngol Head Neck. Surg. 63, 336–342. doi: 10.1007/s12070-011-0266-y

PubMed Abstract | Crossref Full Text | Google Scholar

Jawad, H. M., Ali, R. M., Rehman, S. U., Ansar, M., Gitari, H., Khan, I., et al. (2020). Effect of cadmium toxicity on growth, oxidative damage, antioxidant defense system and cadmium accumulation in two sorghum cultivars. Plants 9, 1575. doi: 10.3390/plants9111575

PubMed Abstract | Crossref Full Text | Google Scholar

Johnson, R., Vishwakarma, K., Hossen, M. S., Kumar, V., Shackira, A. M., Puthur, J. T., et al. (2022). Potassium in plants: Growth regulation, signaling, and environmental stress tolerance. Plant Physiol. Bioch. 172, 56–69. doi: 10.1016/j.plaphy.2022.01.001

PubMed Abstract | Crossref Full Text | Google Scholar

Lea, U. S., Leydecker, M., Quilleré, I., Meyer, C., and Lillo, C. (2006). Posttranslational regulation of nitrate reductase strongly afects the levels of free amino acids and nitrate, whereas transcriptional regulation has only minor infuence. Plant Physiol. 140, 1085–1094. doi: 10.1104/pp.105.074633

PubMed Abstract | Crossref Full Text | Google Scholar

Li, L., Dong, M., Tang, K., Zhang, H., Zeng, N., Yang, H., et al. (2023). The effects of calcium fertilization on morphological and physio-biochemical characteristics in peanut seedlings under waterlogging stress. J. Plant Nutr. 46, 2865–2881. doi: 10.1080/01904167.2022.2160757

Crossref Full Text | Google Scholar

Li, T., Wei, Q., Cui, Y., Tan, H., Zhang, H., Zeng, F., et al. (2024). Foliar application of organic Sorbitol-Chelated calcium promotes peanut (Arachis hypogaea L.) growth and changes phyllosphere bacterial community structure. J. Soil Sci. Plant Nutr. 24, 561–571. doi: 10.1007/s42729-023-01565-0

Crossref Full Text | Google Scholar

Liang, L., Song, Y., Qian, W., Ruan, J., Ding, Z., Zhang, Q., et al. (2021). Metabolomics analysis reveals the responses of tea plants to excessive calcium. J. Sci. Food Agric. 101 (13), 5678–5687. doi: 10.1002/jsfa.11222

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, C., Li, Y., Huo, W., Li, T., Wei, Q., Huang, M., et al. (2021). Effect of sorbitol calcium chelate on yield and calcium nutrient absorption of peanut. Am. J. Biochem. Biotechnol. 17, 160–173. doi: 10.3844/ajbbsp.2021.160.173

Crossref Full Text | Google Scholar

Lopez-Ortiz, C., Dutta, S. K., Natarajan, P., Peña-Garcia, Y., Abburi, V., Saminathan, T., et al. (2019). Genome-wide identification and gene expression pattern of ABC transporter gene family in Capsicum spp. PloS One 14, e0215901. doi: 10.1371/journal.pone.0215901

PubMed Abstract | Crossref Full Text | Google Scholar

Lv, Z., Lan, G., Bai, B., Yu, P., Wang, C., Zhang, H., et al. (2024). Identification of candidate genes associated with peanut pod length by combined analysis of QTL-seq and RNA-seq. Genomics. 116, 110835. doi: 10.1016/j.ygeno.2024.110835

PubMed Abstract | Crossref Full Text | Google Scholar

Ma, J., Saleem, M. H., Ali, B., Rasheed, R., Ashraf, M. A., Aziz, H., et al. (2022). Impact of foliar application of syringic acid on tomato (Solanum lycopersicum L.) under heavy metal stress-insights into nutrient uptake, redox homeostasis, oxidative stress, and antioxidant defense. Front. Plant Sci. 13, 950120. doi: 10.3389/fpls.2022.950120

PubMed Abstract | Crossref Full Text | Google Scholar

Pompelli, M. F., Martins, S. C. V., Antunes, W. C., Agnaldo, R. M., Chaves, A. R. M., and DaMatta, F. M. (2010). Photosynthesis and photoprotection in coffee leaves is affected by nitrogen and light availabilities in winter conditions. J. Plant Physiol. 167, 1052–1060. doi: 10.1016/j.jplph.2010.03.001

PubMed Abstract | Crossref Full Text | Google Scholar

Rahman, A., Mostofa, M. G., Alam, M. M., Nahar, K., Hasanuzzaman, M., and Fujita, M. (2015). Calcium mitigates arsenic toxicity in rice seedlings by reducing arsenic uptake and modulating the antioxidant defense and glyoxalase systems and stress markers. BioMed. Res. Int. 2015, 340812. doi: 10.1155/2015/340812

PubMed Abstract | Crossref Full Text | Google Scholar

Sadras, V. O., Thomas, D., Cozzolino, D., Cossani, D., and Mariano, C. (2019). Wheat yield response to nitrogen from the perspective of intraspecific competition. Field Crop Res. 243, 107632. doi: 10.1016/j.fcr.2019.107632

Crossref Full Text | Google Scholar

Shi, J. and Du, X. (2023). Transcriptome analysis reveals the regulation of cyclic nucleotide-gated ion channels in response to exogenous abscisic acid and calcium treatment under drought stress in tomato. Front. Genet. 14, 1139087. doi: 10.3389/fgene.2023.1139087

PubMed Abstract | Crossref Full Text | Google Scholar

Singh, A., Singh, S., Reddy, K. K., Singh, A. K., Pal, K. K., Dey, R., et al. (2025). Salt stress alleviation in peanut through calcium supplementation. Acta Physiol. Plant 47, 41. doi: 10.1007/s11738-025-03793-6

Crossref Full Text | Google Scholar

Tomeo, N. J. and Rosenthal, D. M. (2017). Variable mesophyll conductance among soybean cultivars sets a tradeoff between photosynthesis and water-use-efficiency. Plant Physiol. 174, 241–257. doi: 10.1104/pp.16.01940

PubMed Abstract | Crossref Full Text | Google Scholar

Vuković, R., Čamagajevac, I.Š, Vuković, A., Šunić, K., Begovi´c, L., Mlinarić, S., et al. (2022). Physiological,biochemical and molecular response of different winter wheat varieties under drought stress at germination and seedling growth stage. Antioxidants 11, 693. doi: 10.3390/antiox11040693

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, L., Fang, C., and Wang, K. (2015). Physiological responses of Leymus chinensis to long-term salt, alkali and mixed salt-alkali stresses. J. Plant Nutr. 38 (4), 526–540. doi: 10.1080/01904167.2014.937874

Crossref Full Text | Google Scholar

Wang, L., Ko, E. E., Tran, J., and Qiao, H. (2020). TREE1-EIN3-mediated transcriptional repression inhibits shoot growth in response to ethylene. Proc. Natl. Acad. Sci. U S A. 117, 29178–29189. doi: 10.1073/pnas.2018735117

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, X., Qi, F., Sun, Z., Liu, H., Wu, Y., Wu, X., et al. (2024). Transcriptome sequencing and expression analysis in peanut reveal the potential mechanism response to Ralstonia solanacearum infection. BMC Plant Biol. 24, 207. doi: 10.1186/s12870-024-04877-0

PubMed Abstract | Crossref Full Text | Google Scholar

Wang, F., Tan, W., Song, W., Yang, S., and Qiao, S. (2022). Transcriptome analysis of sweet potato responses to potassium deficiency. BMC Genomics 23, 655. doi: 10.1186/s12864-022-08870-5

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, X., Hu, W., Luo, H., Xia, Y., Zhao, Y., Wang, L., et al. (2016). Transcriptome profiling of developmental leaf senescence in sorghum (Sorghum bicolor). Plant Mol. Biol. 92, 555–580. doi: 10.1007/s11103-016-0532-1

PubMed Abstract | Crossref Full Text | Google Scholar

Wu, D., Liu, Y., Pang, J., Yong, J., Chen, Y., Bai, C., et al. (2020). Exogenous calcium alleviates Nocturnal Chilling-Induced Feedback inhibition of photosynthesis by improving Sink demand in peanut (Arachis hypogaea). Front. Plant Sci. 11, 607029. doi: 10.3389/fpls.2020.607029

PubMed Abstract | Crossref Full Text | Google Scholar

Xie, Z., Nolan, T. M., Jiang, H., and Yin, Y. (2019). AP2/ERF Transcription factor regulatory networks in hormone and abiotic stress responses in Arabidopsis. Front. Plant Sci. 10, 228. doi: 10.3389/fpls.2019.00228

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, F., Chen, S., Zhou, S., Yue, C., Yang, X., Zhang, X., et al. (2022). Genome-wide association, RNA-seq and iTRAQ analyses identify candidate genes controlling radicle length of wheat. Front. Plant Sci. 13, 939544. doi: 10.3389/fpls.2022.939544

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, M., Li, X., Xie, W., Lin, C., Wang, Q., and Tao, Z. (2023). ETHYLENE INSENSITIVE3/EIN3-LIKE1 modulate FLOWERING LOCUS C expression via histone demethylase interaction. Plant Physiol. 192, 2290–2300. doi: 10.1093/plphys/kiad131

PubMed Abstract | Crossref Full Text | Google Scholar

Xu, Y., Zhang, Z., Guo, Q., Feng, H., Qin, F., Zhang, G., et al. (2025). Effects of nitrogen allocation and photosynthetic proteins response in peanut leaves on photosynthesis under conditions of water scarcity and nitrogen deficiency. BMC Plant Biol. 25, 1366. doi: 10.1186/s12870-025-07232-z

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, L., Wang, C., He, X., Liang, H., Wu, Q., Sun, X., et al. (2024). Multi-year crop rotation and quicklime application promote stable peanut yield and high nutrient-use efficiency by regulating soil nutrient availability and bacterial/fungal community. Front. Microbiol. 15, 1367184. doi: 10.3389/fmicb.2024.1367184

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, S., Li, L., Zhang, J., Geng, Y., Guo, F., Wang, J., et al. (2017). Transcriptome and differential expression profiling analysis of the mechanism of Ca²+ regulation in peanut (Arachis hypogaea) pod development. Front. Plant Sci. 8, 1609. doi: 10.3389/fpls.2017.01609

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, S., Wang, J., Tang, Z., Guo, F., and Li, X. (2020). Transcriptome of peanut kernel and shell reveals the mechanism of calcium on peanut pod development. Sci. Rep. 10, 15723. doi: 10.1038/s41598-020-72893-9

PubMed Abstract | Crossref Full Text | Google Scholar

Yu, Z., Chen, H., Searle, E. B., and Huang, Z. (2020). Whole soil acidification and base cation reduction across subtropical China. Geoderma 361, 114107. doi: 10.1016/j.geoderma.2019.114107

Crossref Full Text | Google Scholar

Zhang, Z., Sun, H., Shao, C., Lei, H., Qian, J., Ruan, Y., et al. (2022). Calcium affects growth and physiological indices of Panax quinquefolium L. HortScience 57, 112–117. doi: 10.21273/HORTSCI16129-21

Crossref Full Text | Google Scholar

Zhao, Y., Ma, J., Li, M., Deng, L., Li, G., Xia, H., et al. (2020). Whole-genome resequencing-based QTL-seq identified AhTc1 gene encoding a R2R3-MYB transcription factor controlling peanut purple testa colour. Plant Biotechnol. J. 18, 96–105. doi: 10.1111/pbi.13175

PubMed Abstract | Crossref Full Text | Google Scholar

Zhao, K., Zhang, J., Fan, Y., Du, X., Zhu, S., Li, Z., et al. (2025). PSC1, a basic/helix-loop-helix transcription factor controlling the purplish-red testa trait in peanut. J. Integr. Plant Biol. 67, 1364–1378. doi: 10.1111/jipb.13847

PubMed Abstract | Crossref Full Text | Google Scholar

Zuo, G., Mei, W., Feng, N., and Zheng, D. (2025). Photosynthetic performance index (PIabs) and malondialdehyde (MDA) content determine rice biomass under combined salt stress and prohexadione-calcium treatment. BMC Plant Biol. 25, 823. doi: 10.1186/s12870-025-06826-x

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: peanut, calcium, transcriptomic, oxidative, nitrogen efficiency

Citation: Xu F, Li L, Si X, Suo Y, Wang X, Zhang Z, Li Q and Zhang X (2026) Physiological and transcriptomic analyses of exogenous calcium in boosting nitrogen use efficiency via oxidative and resistance pathways in peanuts. Front. Plant Sci. 16:1629610. doi: 10.3389/fpls.2025.1629610

Received: 16 May 2025; Accepted: 08 December 2025; Revised: 08 November 2025;
Published: 27 January 2026.

Edited by:

Jianling Fan, Nanjing University of Information Science and Technology, China

Reviewed by:

Feng Xu, Yangtze University, China
Ying Liu, Guangdong Ocean University, China

Copyright © 2026 Xu, Li, Si, Suo, Wang, Zhang, Li and Zhang. 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: Xiang Zhang, enh0ZjIwM0AxNjMuY29t

Disclaimer: 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.