Effects of phosphorus supply levels on the growth and nutrient utilization of Pinus massoniana lamb hybrid saplings

Introduction: Phosphorus is vital for plant growth but is often limited in subtropical acidic soils, restricting forest productivity. Pinus massoniana, a key timber and resin species in China, shows varied shoot-sprouting characteristics and nutrient use under different phosphorus levels. Understanding its response to fertilization is crucial for optimizing plantation management and growth.

Methods: We employed a two-factor randomized block design with three Pinus massoniana hybrid varieties (with shared male or female parents) and four phosphate fertilizer levels (0, 100, 250, 500 g·plant⁻¹). After three growing seasons, tree height, diameter at breast height (DBH), annual shoot growth characteristics (fixed growth, free growth, shoot number), and growth amount were measured. Nitrogen and phosphorus contents in coniferous leaves were also analyzed.

Results: Results showed significant differences in tree height, DBH, and fixed growth among varieties and fertilizer treatments. The 500 g·plant⁻¹ treatment performed best, with fixed growth 14.35% higher than the control. Free growth and shoot flushing frequency showed no significant differences among treatments. The interaction between variety and fertilization contributed 74.10% to free growth variation. Phosphorus fertilization enhanced growth and nutrient uptake in young Pinus massoniana saplings. Increasing fertilizer intensity raised nitrogen and phosphorus contents and the nitrogen-phosphorus ratio. Correlation analysis revealed that the nutrient environment influences genetic expression among traits, affecting their interrelationships.

Discussion: This study examines how genotype, phosphate fertilizer intensity, and their interaction influence the growth and shoot-sprouting behavior of Pinus massoniana hybrids. Phosphate fertilization is a key driver of early growth and nutrient accumulation, while shoot-sprouting is largely shaped by variety-fertilizer interactions. Findings highlight the need to focus less on the genetic basis of shoot growth and more on how nutrient availability influences variety performance, supporting precise fertilization strategies for optimized forest management.

GRAPHICAL ABSTRACT

Graphical Abstract.

Highlights

● Phosphorus fertilizer application significantly promoted the growth of Pinus massoniana hybrid varieties, increasing tree height by 6.14%, DBH by 5.49%, and fixed growth by 14.35% at the 500 g·plant^-¹ application rate.

● The interaction between phosphorus fertilizer and variety significantly affected free growth, with a genetic coefficient of variation of 11.54% and an environmental coefficient of variation of 52.99%.

● Applying phosphate fertilizer is the key factor influencing the initial growth and nutrient accumulation of Pinus massoniana hybrids. Multiple shoot-sprouting behavior is primarily driven by the interaction between variety and phosphate fertilizer, while the nutritional environment affects genetic trait correlations.

1 Introduction

Pinus massoniana is one of the most important fast-growing timber species in subtropical China and a major resin-producing species, widely distributed in low hills and acidic red soil areas (Benarchid et al., 2022). However, these areas generally suffer from soil infertility and low phosphorus (P) availability, limiting the growth potential of major afforestation species such as P. massoniana, Cunninghamia, and Eucalyptus in the same region (Qiang et al., 2021). With the improvement of plantation management systems and the increase in intensive management levels, less land is needed to provide higher-yielding forest products (Klevtsov et al., 2021). Phosphorus is one of the essential elements for plant growth, playing a crucial role in cell energy metabolism, nucleic acid synthesis, root development and shoot architecture (Hanlon et al., 2018). In plantation management, the application of phosphorus fertilizer is usually an effective way to improve soil nutrients and promote forest growth (Sellin et al., 2024; Talboys et al., 2014). Currently, high-gain varieties obtained through breeding are mainly used in production, and fertilization during the sapling stage is the most commonly used plantation management measure (Bera et al., 2018). In the future, intensive or drone fertilization may become an important direction for plantation management (Gonzalez et al., 2013). Therefore, how to improve the growth potential of plantations through scientific fertilization and mastering the interaction between varieties and soil nutrients has become an important issue in current forestry research.

Shoot growth characteristics are a crucial component of tree growth patterns, directly influencing growth rate and wood productivity (Zhao et al., 2019; Zhang et al., 2020). Variations in shoot-sprouting behavior are a common strategy for P. massoniana a sapling to respond to changes in nutrient environments, depending on the effects of environment or genotype, or the interaction between genotype and environment (Martens et al., 2019). P. massoniana has entered the third round of genetic improvement, and genetic differences often lead to significant variations in tree height and diameter at breast height (DBH) among different varieties (Zhang et al., 2020). Although the genetic background of annual shoot traits may be somewhat offset by multiple rounds of genetic selection, the number of annual shoots and shoot growth characteristics are more susceptible to environmental factors (Zhao et al., 2019). For example, Martens et al. (2019) pointed out that pine trees exhibit a periodic growth pattern during annual growth, and free shoot growth is greatly affected by changes in environmental conditions. Especially when the nutrient environment is unstable, the free growth of trees will show greater environmental adaptability. Sellin et al. (2024) indicated that free shoot growth in trees is significantly affected by environmental conditions (such as nutrient environment), especially when the fertilization intensity is high, showing strong phenotypic plasticity, which is closely related to changes in different fertilization treatments. Different strains show significant differences in growth response to phosphorus fertilizer application, which is closely related to the role of genotype in nutrient absorption, utilization, and distribution (Guo W. et al., 2023). Specifically, phosphorus availability can alter the balance between primary and secondary shoot growth, as trees allocate resources to optimize nutrient capture under varying soil conditions. However, there are few studies on the interaction between high-generation hybrid varieties of P. massoniana and phosphorus fertilizer application on shoot traits, and the effects of genotype, phosphorus fertilizer, and their interaction on shoot-sprouting behavior and nutrient distribution have not been fully confirmed.

Needle nutrient content is widely recognized as a crucial indicator of tree adaptation to nutrient-poor environments, directly affecting tree growth rate, leaf longevity, and other processes, which are key to long-term tree productivity (Zhang et al., 2020). The nitrogen-to-phosphorus ratio (N/P ratio) is considered an important ecological indicator of plant nutrient limitation, closely related to forest growth strategies and adaptability (Guo Q. et al., 2023). The sensitivity of shoot-sprouting behavior and mechanisms to nutrient environments varies among genotypes, and the correlation between changes in shoot-sprouting behavior and needle nutrient content has not been thoroughly investigated. Therefore, we selected three hybrid varieties and established four phosphorus fertilizer application levels to quantify the growth characteristics, annual shoot behavior, and needle nutrient dynamics of P. massoniana saplings under different phosphorus fertilizer intensities. The main objectives were: (1) to investigate the effects of genotype, phosphorus fertilizer intensity, and their interaction on the growth and shoot growth characteristics of P. massoniana saplings; (2) to reveal the regulatory mechanism of phosphorus fertilizer on needle nutrient accumulation and utilization patterns; and (3) to establish the relationship between phosphorus fertilizer levels and genotype adaptability, providing a basis for P. massoniana variety promotion and precision fertilization management.

2 Materials and methods

2.1 Study area

The experiment was conducted at the Experimental Forest Farm in Qingyuan County, Zhejiang Province, China (119°01’25″ E, 27°38’48″ N), with an altitude ranging from 480 m to 510 m. The general situation of the research area is shown in Figure 1. The experimental site has a gentle slope. The region has an average annual temperature of 17.6 °C, an annual precipitation of 1721.3 mm, and a frost-free period of 245 days, characteristic of a typical subtropical monsoon climate. The experimental site comprises acidic red soil characterized by a pH of 4.67. Key soil properties include 7.8 g/kg organic matter, 0.45 g/kg total nitrogen, 0.27 g/kg total phosphorus, and available nutrient concentrations of 17.2 mg/kg nitrogen, 37.3 mg/kg potassium, and 1.02 mg/kg phosphorus.

Figure 1

Figure 1. General situation of study area.

2.2 Materials and design

A completely randomized block design with two factors (variety and fertilization) was used. Three P. massoniana hybrid varieties were selected: F32 (♀6627×♂1003), F35 (♀5163×♂1003), and F36 (♀5163×♂6627). These hybrid families were generated through artificial controlled pollination within a second-generation breeding population, with all three varieties sharing either a common paternal or maternal parent. The selected genotypes represented distinct geographic provenances: (1) genotype 6627, derived from Jiangxi provenance (central production region), exhibiting 2–3 annual shoot flushes; (2) genotype 1003, originating from Guangdong provenance (southernmost production region), characterized by 3–5 annual shoot flushes; and (3) genotype 5163, representing Zhejiang provenance (northern production region), typically showing 1–2 annual shoot flushes. Four phosphorus fertilizer treatments were established: 0 g·plant^-¹ (T0, as control), 100 g·plant^-¹ (T1), 250 g·plant^-¹ (T2), and 500 g·plant^-¹ (T3). The applied phosphorus fertilizer was diammonium phosphate ((NH₄)₂HPO₄: containing 46% phosphorus and 18% nitrogen). In March 2020, a measurement forest was established according to the scheme shown in Figure 2. The 12 treatment combinations were replicated five times (blocks), with each block containing all treatments randomly arranged. Each treatment plot comprised 40 trees (4 rows × 10 plants, 2 m × 2 m spacing), totaling 2,400 trees (12 treatments × 5 blocks × 40 trees). Fertilization was applied in May 2020–2022 via a continuous circular ditch (15 cm from trunk, 20 cm deep, 10 cm wide) around each tree. Growth traits (fixed growth: predetermined bud flushes; free growth: post-flush shoot elongation) were assessed for all 40 trees per plot (no subsampling). Protective borders of open-pollinated seedlings surrounded each block. Annual weeding ensured minimal competition, with no additional nutrient/water inputs.

Figure 2

Figure 2. Schematic arrangement of the experimental design. Each compartment contains a total of 40 individual plants, which are planted in 4 rows with 10 plants in each row. F32, ♀6627×♂1003; F35, ♀5163×♂1003; F36: ♀5163×♂6627. T0 (CK) was the unfertilized control; T1, T2 and T3 gradients were applied with 100, 250 and 500 g·plant^-¹ of diammonium phosphate.

2.3 Survey and measurement

After the third growing season, tree height, diameter at breast height (DBH), and annual shoot characteristics were measured for each tree. The measurement of annual shoot characteristics was based on the whorled growth pattern of P. massoniana. After overwintering, the previous year’s winter buds began to sprout in the following spring. Shoot growth was determined by measuring the length between two spring branch whorls and the number of branch whorls, thereby calculating the annual height growth and shoot number. The growth of the initial shoot (fixed growth) was the distance from the base of the spring branch whorl to the base of the first shoot. For free whorl branches, the measurement was carried out upward in the same way to obtain the annual shoot number and growth. In addition, four individual trees were randomly selected from each treatment combination in each replicate, and the current year’s needles from the upper part of the tree crown were collected and brought back to the laboratory. The needles were deactivated at 105 °C for 30 minutes, then dried at 80°C to constant weight, and their mass was measured for the determination of nitrogen and phosphorus content. Phosphorus content was determined by the molybdenum-antimony colorimetric method, and nitrogen content was determined by the Kjeldahl method (Bremner, 1960; Zhang et al., 2025).

2.4 Data analysis

We conducted comprehensive analyses of growth traits, shoot morphology, and needle nutrient content in P. massoniana saplings across hybrid varieties and phosphorus (P) fertilization treatments. The descriptive statistics (mean, range, standard deviation) and coefficients of variation were calculated for all measured parameters. The calculation formulas are as follows:

$\begin{array}{ll}\text{Average value}:\overline{x}=\mathit{\sum }xi/n& (1)\end{array}$

$\begin{array}{ll}\text{Genetic Coefficient of Variation}:\text{GCV}=100\%\times \sqrt{{\delta }_G}/\overline{ x}& (2)\end{array}$

$\begin{array}{ll}\begin{array}{c}\text{Environmental Coefficient of Variation}:\text{ECV}\\ =100\%\times \sqrt{{\delta }_E}/\overline{ x}\end{array}& (3)\end{array}$

Where $xi$ is the measured value for each individual plant, $\overline{x}$ is the mean value for a trait, n is the number of individual plants, ${\delta }_G$ and ${\delta }_E$ represent genetic and environmental variance components respectively, estimated using restricted maximum likelihood (REML) methods in a mixed-effects model framework.

The Shapiro-Wilk test was employed to assess the normal distribution of measured data across different hybrid varieties and phosphorus fertilizer application rates. One-way analysis of variance (ANOVA) and Tukey’s multiple comparison test were utilized to analyze the effects of varying phosphorus fertilizer intensities and varieties on needle nitrogen and phosphorus content. Partial least squares (PLS) regression was employed to examine multivariate relationships between functional traits and nutrient content across P gradients, addressing potential collinearity among predictors. All analyses were implemented in R 4.3.2. Effect sizes (η² for ANOVA, r for correlations) are reported alongside p-values.

3 Results

3.1 Changes in growth characteristics

After three years of cultivation, the P. massoniana exhibited an average tree height of 2.68 m, a diameter at breast height (DBH) of 24.22 mm, 1.12 flushes per year, a fixed growth of 85.02 cm, and a free growth of 42.39 cm. The free growth trait showed the highest genetic coefficient of variation (GCV=11.54%) and environmental coefficient of variation (ECV=52.99%), indicating that it has substantial plasticity. The growth and flushing characteristics of the hybrid varieties responded differently to phosphorus fertilizer application. Significant differences were observed among varieties and phosphorus treatments in tree height, DBH, and fixed growth (Table 1, Figure 3A). The T3 treatment showed the best performance in tree height, DBH, and fixed growth, with increases of 6.14%, 5.49%, and 14.35%, respectively, compared to the T0 treatment (Figure 4). This indicates that higher phosphorus fertilizer content significantly promoted the development of primary growth morphology in P. massoniana. No significant differences were found in free growth and annual flushing frequency among varieties and phosphorus treatments, suggesting that the multiple flushing behavior of hybrid varieties was not highly sensitive to phosphorus fertilizer.

Table 1

Table 1. Trait performance and analysis of variance.

Figure 3

Figure 3. Influence of varieties, fertilization and the interaction effect of the two on the measured traits. (A) analysis of variance (ANOVA) for each trait; (B) percentage of variance components of genetics, environment, and both interactions for each trait. V, variety; P, phosphate fertilizer; V×P, variety× phosphate fertilizer; D for diameter at breast height, H for tree height, Fix for fixed growth, Free for free growth, N for nitrogen content of needles, P for phosphorus content of needles, N/P for nitrogen/phosphorus ratio of needles.

Figure 4

Figure 4. Effects of different phosphorus fertilizer intensities on functional traits and nutrient contents of P. massoniana. T0 (CK) was the unfertilized control; T1, T2 and T3 gradients were applied with 100, 250 and 500 g·plant^-¹ of diammonium phosphate, respectively. * indicates P< 0.05; ** indicates P< 0.01; *** indicates P< 0.001.

3.2 Changes in conifer nutrient content

The genetic coefficients of variation for needle nitrogen and phosphorus content and the nitrogen-to-phosphorus ratio were relatively small (1.03%-3.80%), lower than the environmental coefficients of variation (10.85%-19.75%) (Table 1, Figure 3A). No significant differences were observed in needle nutrient content among varieties, but significant differences were found among fertilizer treatments (P<0.05). With increasing fertilizer intensity, needle nitrogen and phosphorus content and the nitrogen-to-phosphorus ratio showed an increasing trend (Figure 4), indicates that increased external phosphorus supply promoted uptake of both nitrogen and phosphorus in P. massoniana needle.

3.3 Factor contribution analysis

The distribution of all indicators varied across three dimensions: variety (V), fertilization (P), and their interaction (V×P) (Figure 3B). The effect on DBH was in the order of variety effect > interaction effect (V×P) > fertilization (P) effect, with the variety effect contributing 51.22%. The effects on tree height, fixed growth, needle nitrogen content, phosphorus content, and nitrogen-to-phosphorus ratio were in the order of fertilization (P) effect > interaction effect (V×P) > variety effect. The contribution of the fertilization (P) effect ranged from 55.43% to 93.62%, and the contribution of the interaction effect (V×P) ranged from 3.13% to 22.42%. This indicates that fertilization is the dominant factor influencing initial growth and nutrient accumulation, while also being affected by the interaction effect. The interaction effect between variety and fertilization (V×P) had the highest contribution to free growth, reaching 74.1%, suggesting that the flushing behavior of different varieties exhibits heterogeneity in response to fertilization intensity.

The fixed growth of hybrid varieties 32 and 35 increased under fertilization treatments, peaking at T1 and T3 treatments, respectively; in contrast, clone 36 showed the opposite trend. The effects of different treatments on the free growth of each variety varied considerably (Figure 5). This indicates that the increase in annual flushing length is attributed to the accumulation of fixed growth and free growth. Hybrid combinations with the same paternal or maternal parent exhibited different response patterns to phosphorus input. However, regardless of low or high phosphorus fertilizer input, the increase in annual flushing length of high-generation hybrid varieties was mainly attributed to the increase in fixed growth.

Figure 5

Figure 5. Effects of different phosphate fertilizer intensities and family types on functional traits and nutrient contents of P. massoniana. T0 (CK) was the unfertilized control group; T1, T2, and T3 gradients were applied 100, 250, and 500 g·plant^-¹ of diammonium phosphate, respectively. Only intergroup differences, i.e., differences among family types, are compared in the figure.

3.4 Correlation between phosphorus fertilizer application and plant traits

As the amount of phosphorus fertilizer application increased, more correlations were observed between tree height, DBH, and nutrient factors (Figure 6). Flushing frequency showed mostly significant negative correlations with other traits (P< 0.05), especially under T1 and T2 treatments. The N:P ratio was positively correlated with tree height and DBH, but negatively correlated with other traits. Furthermore, under the T4 treatment, with increasing fertilization, nitrogen, phosphorus, and the N:P ratio showed positive correlations with multiple growth indicators, particularly with tree height (P< 0.01). This indicates that the nutritional environment influences the genetic expression of traits, thereby affecting the interrelationships between traits.

Figure 6

Figure 6. Genetic correlations between functional traits and nutrient content. T0 (CK) was the unfertilized control; T1, T2, and T3 gradients were applied with 100, 250, and 500 g·plant^-¹ of diammonium phosphate, respectively. * indicates P< 0.05; ** indicates P< 0.01; *** indicates P< 0.001. H for tree height, D for diameter at breast height, Draws for number of draws, Fix for fixed growth, Free for free growth, N for nitrogen content of needles, P for phosphorus content of needles, N/P for nitrogen/phosphorus ratio of needles.

4 Discussion

4.1 Regulation mechanism of phosphorus fertilizer intensity on needle nutrient uptake and utilization

Phosphorus fertilizer application significantly promoted the growth of P. massoniana young forests, notably improving traits such as tree height, DBH, and fixed growth (Figure 3). As a key nutrient, phosphorus input also increased nitrogen and phosphorus contents in needles, especially phosphorus uptake (Figure 4), directly enhancing biomass accumulation and morphogenesis during early growth (Na et al., 2014; Yi et al., 2023). Additionally, phosphorus input indirectly improved nitrogen use efficiency by enhancing soil nutrients or root activity (Mulligan and Sands, 1988). Growth of P. massoniana is limited by nitrogen and phosphorus deficiency in southern acidic soils, particularly phosphorus (Xu et al., 2018). This study confirms that exogenous phosphorus alleviates this limitation, improves soil nutrients and plant uptake, and significantly boosts P. massoniana growth, especially under low soil phosphorus conditions (Xu et al., 2018). Research has shown that fertilization’s effect on needle nutrients is mainly driven by environmental factors, with fertilization treatments accounting for 52.99%–64.08% of the variation in nitrogen and phosphorus content (Table 1), while genotypic influence is relatively limited. This supports Guo W. et al. (2023), who suggested needle nutrient content reflects environmental adaptation more than genotype. Under significant environmental changes, plants exhibit high plasticity in nutrient uptake and allocation, highlighting the importance of fertilization management in forestry (Rubio and Lavado, 1999; Andersen et al., 2017). Fertilization supports optimal growth in variable environments (Freiberger et al., 2013; Niu et al., 2020; Rennenberg et al., 2009).

Under low phosphorus conditions, a negative correlation was observed between DBH and needle phosphorus content (Figure 5); as phosphorus supply increased, this shifted to a significant positive correlation, reflecting P. massoniana nutrient allocation strategy under phosphorus stress. When phosphorus is limited, the species may prioritize allocation to growth organs (e.g., roots, stems), reducing leaf phosphorus to sustain growth (Turner, 2008). This adaptive mechanism allows plants to maintain growth under stress (DeLoose et al., 2024). By analyzing needle nitrogen and phosphorus content under different phosphorus levels and their correlations with traits like flushing growth and frequency, the study found that nutrient environment strongly influences trait expression and interactions. Without phosphorus fertilizer, traits such as tree height and DBH showed strong positive correlations with other growth traits, which weakened with increased phosphorus input. Thus, appropriate phosphorus fertilization not only alters needle nutrient levels but also indirectly affects flushing traits via the nutritional environment, deepening our understanding of P. massoniana growth mechanisms under phosphorus variation.

4.2 Effects of genotype-environment interaction on growth traits of P. massoniana

This study revealed significant differences in the growth responses of P. massoniana hybrid varieties to phosphorus fertilizer intensity, highlighting the critical role of genotype–environment interaction (V×P) in shaping young forest traits. Hybrid combination 32 performed best under moderate phosphorus treatment (T2), with significantly greater tree height and DBH, while combinations 35 and 36 showed better growth under high (T3) and low (T1) phosphorus levels, respectively (Figure 4). This demonstrates genotypic variation in phosphorus response and the influence of fertilization intensity on growth potential (Debell et al., 1986). Combination 32, derived from southernmost and central provenances, exhibited high growth potential, while combinations 35 and 36, with northern provenance genes, showed different responses—indicating that growth performance is linked to parent origin. These findings align with Stovall et al. (2011), who reported that genotype × fertilization interaction significantly influenced height and stem growth in slash pine. Other studies confirm that environmental factors, especially fertilization intensity, not only directly affect growth but also interact with genotype responses, emphasizing the importance of G×E effects in forest development (Keller et al., 2018; Zhao et al., 2016).

Tree height showed higher heritability (h² = 0.37) than free growth (h² = 0.30) (Table 1), suggesting stronger genetic control over apical growth, while free growth was more influenced by environmental factors, particularly fertilization intensity. Its high plasticity (CV = 54.08%) reflects sensitivity to phosphorus availability, with trees under low phosphorus prioritizing structural growth over flushing—evidenced by the negative DBH–needle P correlation. This resource allocation is hormonally regulated via gibberellin–auxin interactions and elevated ABA under phosphorus stress (Liu et al., 2024). These patterns align with the “source–sink redistribution” model (Suriyagoda et al., 2023), highlighting the need to consider fertilization effects on plastic traits like free growth when managing P. massoniana. Genetic variation in flushing traits (h² = 0.30) suggests differential phosphorus sensitivity among clones, shaped by provenance adaptation. Free growth was most affected by variety × fertilization interaction (74.1%), while annual flushing frequency showed low responsiveness, indicating it is less sensitive to phosphorus input. Overall, fertilization interacts with genotype to regulate growth plasticity in P. massoniana.

This study found that hybrid combination 32 performed best under moderate phosphorus levels, while combinations 35 and 36 showed better growth under high and low phosphorus, respectively. This suggests that phosphorus concentration has genotype-specific effects on P. massoniana growth, with each hybrid maximizing its potential under optimal conditions (Conroy et al., 1990). Such differences likely arise from genotypic variation in nutrient absorption, utilization, and allocation (Liu et al., 2023; McKeand et al., 2023). Therefore, selecting genotypes suited to specific phosphorus levels can enhance fertilization efficiency and forest productivity (Tie et al., 2024; Yang et al., 2025).

4.3 Practical significance of optimizing P. massoniana plantation management

Shoot-sprouting behavior can reflect the growth rate of P. massoniana, with fixed growth and free growth being closely related to genetic and environmental factors, respectively. For hybrid varieties with high-generation genetic improvement, the variation in flushing frequency may not be obvious, but mainly manifested as changes in fixed growth. This means that during the genetic improvement process, the genetic stability of varieties is enhanced, and the buffering capacity against environmental changes is improved, making the flushing frequency relatively stable. Free growth, on the other hand, is more influenced by environmental stimuli. Under different phosphorus fertilizer intensities, the variation in free growth among different clones varied greatly and showed no obvious pattern, reflecting the complex role of environmental factors on it. In forestry practice, this study provides a theoretical basis for “precision genotype-environment matching”. For example, hybrid combination 32 is preferentially planted on sites with moderate phosphorus availability to maximize its root efficiency and phosphorus absorption potential; while hybrid combination 35 is promoted in high-phosphorus soils to avoid phosphorus toxicity risks and optimize resource utilization. This is consistent with the “suitable genotype for suitable site” strategy proposed by Chen et al. (2016), which emphasizes improving plantation productivity and sustainability through precise matching of genotype and environment.

In addition, in actual forestry production, rational fertilization is one of the effective ways to improve the nutrient status and promote the growth of P. massoniana (Mao et al., 2021; Guo Q. et al., 2023). More attention should be paid to the regulation of DBH growth, and the timber yield of P. massoniana should be increased by optimizing fertilization and selecting suitable genotypes (Qiao et al., 2022). In the future, genomic selection technology should be combined to screen families with high phosphorus response potential, which can accelerate the genetic improvement process of P. massoniana. Future research can combine genomics approaches (such as transcriptomics and metabolomics) to analyze the molecular mechanisms of nutrient response in P. massoniana, and monitor the feedback effects of fertilization on soil fertility in the long term, to achieve sustainable management of plantations.

5 Conclusion

This study revealed the synergistic effects of hybrid varieties and phosphorus fertilizer intensity on the growth and nutrient dynamics of Pinus massoniana. Phosphorus fertilizer promoted young forest growth and needle nutrient absorption. With the increase in fertilization intensity, needle nitrogen and phosphorus content and the nitrogen-to-phosphorus ratio showed an increasing trend. Phosphorus fertilization is the dominant factor significantly affecting the initial growth and nutrient accumulation of Pinus massoniana hybrid varieties, while the multiple shoot-sprouting behavior of hybrid varieties is mainly affected by the interaction effect of variety and phosphorus fertilizer, showing different plastic responses to phosphorus fertilizer. The contribution rate of the interaction effect between variety and fertilization to free growth reached 74.1%. This suggests that when promoting and applying Pinus massoniana varieties with generational improvement, the consideration of the genetic behavior of shoot-sprouting behavior is gradually reduced, and the interaction between varieties and nutrient environment should be given more attention.

Data availability statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author contributions

WW: Conceptualization, Visualization, Writing – original draft. JH: Data curation, Formal Analysis, Writing – original draft. ZZ: Investigation, Methodology, Writing – review & editing. ZCZ: Project administration, Resources, Writing – review & editing. RW: Investigation, Validation, Writing – review & editing. S: Investigation, Validation, Writing – review & editing. XC: Resources, Supervision, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This work was funded by National Key Research and Development Project during the 14th Five-year Plan Period ( 2022YFD2200202; 2023YFD2200601). The project funding is jointly funded by two projects.

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.

Generative AI statement

The author(s) declare that no Generative AI was used in the creation of this manuscript.

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.

References

Andersen, K. M., Mayor, J. R., and Turner, B. L. (2017). Plasticity in nitrogen uptake among plant species with contrasting nutrient acquisition strategies in a tropical forest. Ecology 98, 1388–1398. doi: 10.1002/ecy.1793

PubMed Abstract | Crossref Full Text | Google Scholar

Benarchid, K., Khatori, M., and Hilali, S. (2022). Effect of pinus halepensis reforestation on soil fertilityin the forest of benisohane (ribat al kheir region -Morocco). Ecology Environ. Conserv. 28, 78–85. doi: 10.53550/eec.2022.v28i01.011

Crossref Full Text | Google Scholar

Bera, S. K., Kamdar, J. H., Kasundra, S. V., Dash, P., Maurya, A. K., Jasani, M. D., et al. (2018). Improving oil quality by altering levels of fatty acids through marker-assisted selection of ahfad2 alleles in peanut (arachis hypogaea L.). Euphytica 214, 162. doi: 10.1007/s10681-018-2241-0

Crossref Full Text | Google Scholar

Bremner, J. M. (1960). Determination of nitrogen in soil by the kjeldahl method. J. Agric. Sci. 55, 11–33. doi: 10.1017/s0021859600021572

Crossref Full Text | Google Scholar

Chen, M., Rao, R. S. P., Zhang, Y., Zhong, C., and Thelen, J. J. (2016). Metabolite variation in hybrid corn grain from a large-scale multisite study. Crop J. 4, 177–187. doi: 10.1016/j.cj.2016.03.004

Crossref Full Text | Google Scholar

Conroy, J., Milham, P., Reed, M., and Barlow, E. (1990). Increases in phosphorus requirements for Co2-enriched pine species. Plant Physiol. 92, 977–982. doi: 10.1104/pp.92.4.977

PubMed Abstract | Crossref Full Text | Google Scholar

Debell, D., Silen, R., Radwan, M., and Mandel, N. (1986). Effect of family and nitrogen-fertilizer on growth and foliar nutrients of douglas-fir saplings. For. Science. 32, 643–652. doi: 10.1093/forestscience/32.3.643

Crossref Full Text | Google Scholar

DeLoose, M., Clua, J., Cho, H., Zheng, L., Masmoudi, K., Desnos, T., et al. (2024). Recent advances in unraveling the mystery of combined nutrient stress in plants. Plant J. 117, 1764–1780. doi: 10.1111/tpj.16511

PubMed Abstract | Crossref Full Text | Google Scholar

Freiberger, M. B., Guerrini, I. A., Galetti, G., Fernandes, D. M., and Correa, J. C. (2013). Early growth and nutrition of cedar (Cedrela fissilis Vell.) as affected by nitrogen rates. Rev. Árvore. 37, 385–392. doi: 10.1590/S0100-67622013000300001

Crossref Full Text | Google Scholar

Gonzalez, E., Hébert, F, Boucher, J, Sirois, P., and Lord, D. (2013). Lichen-spruce woodland early indicators of ecological resilience following silvicultural disturbances in québec’s closed-crown forest zone. Am. J. Plant Sci. 4, 749–765. doi: 10.4236/ajps.2013.43A094

Crossref Full Text | Google Scholar

Guo, W., Jiao, P., Lopez, C. M. L., Chen, Z., Wei, L., Liu, X., et al. (2023). Nitrogen and phosphorous dynamics with stand development of pinus massoniana plantations in southeast China. Front. Plant Sci. 14. doi: 10.3389/fpls.2023.1139945

PubMed Abstract | Crossref Full Text | Google Scholar

Guo, Q., Li, H., Sun, X., An, Z., and Ding, G. (2023). Patterns of needle nutrient resorption and ecological stoichiometry homeostasis along a chronosequence of pinus massoniana plantations. Forests. 14, 607. doi: 10.3390/f14030607

Crossref Full Text | Google Scholar

Hanlon, M. T., Ray, S., Saengwilai, P., Luthe, D., Lynch, J. P., and Brown, K. M. (2018). Buffered delivery of phosphate to Arabidopsis alters responses to low phosphate. J. Exp. Bot. 69, 1207–1219. doi: 10.1093/jxb/erx454

PubMed Abstract | Crossref Full Text | Google Scholar

Keller, S. R., Chhatre, V. E., and Fitzpatrick, M. C. (2018). Influence of range position on locally adaptive gene-environment associations in populus flowering time genes. J. Heredity. 109, 47–58. doi: 10.1093/jhered/esx098

PubMed Abstract | Crossref Full Text | Google Scholar

Klevtsov, D. N., Tyukavina, O. N., and Sungurova, N. R. (2021). Comparative analysis of scots pine homogeneous stands bioenergetic productivity in European north. Forestry Bulletin. 25, 15–20. doi: 10.18698/2542-1468-2021-4-15-20

Crossref Full Text | Google Scholar

Liu, S. T., Gille, C. E., Bird, T., Ranathunge, K., Finnegan, P. M., and Lambers, H. (2023). Leaf phosphorus allocation to chemical fractions and its seasonal variation in south-western Australia is a species-dependent trait. Sci. Total Environment. 901, 166395. doi: 10.1016/j.scitotenv.2023.166395

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, Y., Lan, X., Hou, H., Ji, J., Liu, X., and Lv, Z. (2024). Multifaceted ability of organic fertilizers to improve crop productivity and abiotic stress tolerance: review and perspectives. Agronomy. 14, 1141. doi: 10.3390/agronomy14061141

Crossref Full Text | Google Scholar

Mao, L., Zha, R., Chen, S., Zhang, J., Jie, L., and Zha, X. (2021). Mixture compound fertilizer and super absorbent polymer application significantly promoted growth and increased nutrient levels in pinus massoniana seedlings and soil in seriously eroded degradation region of southern China. Front. Plant Science. 12. doi: 10.3389/fpls.2021.763175

PubMed Abstract | Crossref Full Text | Google Scholar

Martens, H. J., Sørensen, S., Burow, M., and Veierskov, B. (2019). Characterization of top leader elongation in nordmann fir (abies nordmanniana). J. Plant Growth Regulation. 38, 1354–1361. doi: 10.1007/s00344-019-09938-5

Crossref Full Text | Google Scholar

McKeand, S., Mullin, T., Byram, T., and White, T. (2023). Deployment of Genetically Improved Loblolly and Slash Pines in the South. J. For. 101, 32–37. doi: 10.1093/jof/101.3.32

Crossref Full Text | Google Scholar

Mulligan, D. and Sands, R. (1988). Dry-matter, phosphorus and nitrogen partitioning in 3 eucalyptus species grown under a nutrient deficit. New Phytologist. 109, 21–28. doi: 10.1111/j.1469-8137.1988.tb00213.x

Crossref Full Text | Google Scholar

Na, L., Qiao, Z., Hong, J., Xie, Y., and Zhang, P. (2014). Effect of soluble phosphorus microbial mixed fertilizers on phosphorus nutrient and phosphorus adsorption-desorption characteristics in calcareous cinnamon soil. Chin. J. Appl. Environ. Biol. 20, 662–668. doi: 10.3724/SP.J.1145.2014.07005

Crossref Full Text | Google Scholar

Niu, Y., Zhang, M., Bai, S. H., Xu, Z., Liu, Y., Chen, F., et al. (2020). Successive mineral nitrogen or phosphorus fertilization alone significantly altered bacterial community rather than bacterial biomass in plantation soil. Appl. Microbiol. Biotechnol. 104, 7213–7224. doi: 10.1007/s00253-020-10761-2

PubMed Abstract | Crossref Full Text | Google Scholar

Qiang, L. I., Ying-Xin, H., Dao-Wei, Z., and Shan, C. (2021). Mechanism of the trade-off between biological nitrogen fixation and phosphorus acquisition strategies of herbaceous legumes under nitrogen and phosphorus addition. Chin. J. Plant Ecol. 45, 286–297. doi: 10.17521/cjpe.2020.0241

Crossref Full Text | Google Scholar

Qiao, J., Sun, Y., Pan, L., Luo, M., Ding, Z., and Sun, Z. (2022). Variability in the climate-radial growth correlation of pinus massoniana of different diameter classes. J. Forestry Res. 33, 1781–1792. doi: 10.1007/s11676-022-01472-7

Crossref Full Text | Google Scholar

Rennenberg, H., Dannenmann, M., Gessler, A., Kreuzwieser, J., Simon, J., and Papen, H. (2009). Nitrogen balance in forest soils: nutritional limitation of plants under climate change stresses. Plant Biol. 11, 4–23. doi: 10.1111/j.1438-8677.2009.00241.x

PubMed Abstract | Crossref Full Text | Google Scholar

Rubio, G. and Lavado, R. S. (1999). Acquisition and allocation of resources in two waterlogging-tolerant grasses. New Phytologist. 143, 539–546. doi: 10.1046/j.1469-8137.1999.00482.x

PubMed Abstract | Crossref Full Text | Google Scholar

Sellin, A., Heinsoo, K., Kupper, P., Meier, R., Õunapuu-Pikas, E., Reinthal, T., et al. (2024). Growth responses to elevated environmental humidity vary between phenological forms of picea abies. Front. Forests Global Change 7. doi: 10.3389/ffgc.2024.1370934

Crossref Full Text | Google Scholar

Stovall, J. P., Carlson, C. A., Seiler, J. R., Fox, T. R., and Yanez, M. A. (2011). Growth and stem quality responses to fertilizer application by 21 loblolly pine clones in the virginia piedmont. For. Ecol. Management. 261, 362–372. doi: 10.1016/j.foreco.2010.10.018

Crossref Full Text | Google Scholar

Suriyagoda, L. D. B., Ryan, M. H., Gille, C. E., Dayrell, R. L. C., Finnegan, P. M., Ranathunge, K., et al. (2023). Phosphorus fractions in leaves. New Phytologist. 237, 1122–1135. doi: 10.1111/nph.18588

PubMed Abstract | Crossref Full Text | Google Scholar

Talboys, P. J., Healey, J. R., Withers, P. J. A., and Jones, D. L. (2014). Phosphate depletion modulates auxin transport in triticum aestivum leading to altered root branching. J. Exp. Botany. 65 (17), 5023–5032. doi: 10.1093/jxb/eru284

PubMed Abstract | Crossref Full Text | Google Scholar

Tie, L., Penuelas, J., Huang, C., Sardans, J., Bose, A. K., Ouyang, S., et al. (2024). Phosphorus limitation of pinus massoniana reforestation increases with stand development: Evidence from plant, leaf litter, and soil. Plant Soil. 504, 817–832. doi: 10.1007/s11104-024-06661-3

Crossref Full Text | Google Scholar

Turner, B. L. (2008). Resource partitioning for soil phosphorus: a hypothesis. J. Ecology. 96, 698–702. doi: 10.1111/j.1365-2745.2008.01384.x

Crossref Full Text | Google Scholar

Xu, C., Xiang, W., Gou, M., Chen, L., Lei, P., Fang, X., et al. (2018). Effects of forest restoration on soil carbon, nitrogen, phosphorus, and their stoichiometry in Hunan, southern China. Sustainability. 10, 1874. doi: 10.3390/su10061874

Crossref Full Text | Google Scholar

Yang, C., Jiang, P., Qian, Z., Sun, J., Wu, T., Yang, Z., et al. (2025). Evolutionary history of mixed tree species improved soil nutrient content of pinus massoniana plantation. Plant Soil 21. doi: 10.1007/s11104-024-07197-2

Crossref Full Text | Google Scholar

Yi, H., Hu, S., Zhang, Y., Wang, X., Xia, Z., Lei, Y., et al. (2023). Proper delay of phosphorus application promotes wheat growth and nutrient uptake under low phosphorus condition. Agriculture. 13, 884. doi: 10.3390/agriculture13040884

Crossref Full Text | Google Scholar

Zhang, Z., Jin, G., Feng, Z., Sun, L., and Yuan, C. (2020). Joint influence of genetic origin and climate on the growth of masson pine (pinus massoniana lamb.) in China. Sci. Rep. 10, 4653. doi: 10.1038/s41598-020-61597-9

PubMed Abstract | Crossref Full Text | Google Scholar

Zhang, J., Li, C., Jin, X., Liao, Y., Ma, X., Zhang, Y., et al. (2025). Effects of exogenous phosphorus compounds interact with nitrogen addition on soil organic phosphorus fractions in a meadow steppe. J. Soils Sediments. 27, 1234–1245. doi: 10.1016/j.jsos.2025.03.021

Crossref Full Text | Google Scholar

Zhao, X., Pan, Y., Pei, X., Wang, F., Wang, C., Shao, L., et al. (2019). Forward, backward selection and variation analysis of growth traits in half-sib larix kaempferi families. Silvae Genetica. 68, 1–8. doi: 10.2478/sg-2019-0001

Crossref Full Text | Google Scholar

Zhao, X., Xia, H., Wang, X., Wang, C., Liang, D., Li, K., et al. (2016). Variance and stability analyses of growth characters in half-sib Betula platyphylla families at three different sites in China. Euphytica. 208, 173–186. doi: 10.1007/s10681-015-1617-7

Crossref Full Text | Google Scholar