Variation, coordination, and trade-offs between needle structures and photosynthetic-related traits across five Picea species: consequences on plant growth

Background Picea species are distributed and planted world-wide due to their great ecological and economic values. It has been reported that Picea species vary widely in growth traits in a given environment, which reflects genetic and phenotypic differences among species. However, key physiological processes underlying tree growth and the influencing factors on them are still unknown. Results Here, we examined needle structures, needle chemical components, physiological characteristics and growth traits across five Picea species in a common garden in Tianshui, Gansu province in China: Picea glauca, P. mariana, P. likiangensis, P. koraiensis, and P. crassifolia, among which P. glauca and P. mariana were introduced from North America, P. likiangensis was from Lijiang, Yunan province in China, P. koraiensis was from Yichun, Heilongjiang province in China, and P. crassifolia was native to the experimental site. It was found that nearly all traits varied significantly among species. Tissue-level anatomical characteristics and leaf mass per area (LMA) were affected by needle size, but the variations of them were not associated with the variations in photosynthetic and biochemical capacity among species. Variations in area-based maximum photosynthesis (Pnmax) were affected by stomatal conductance (gs), mesophyll conductance (gm) and biochemical parameters including maximum carboxylation rate (Vcmax), and maximum electron transport rate (Jmax). The fraction of N allocated to different photosynthetic apparatus displayed contrasting values among species, which contributed to the species variations in photosynthetic nitrogen use efficiency (PNUE) and Pnmax. Additionally, all growth traits were positively correlated with Pnmax and PNUE. Conclusion Needle structures are less important than needle biochemical parameters in determining the variations in photosynthetic capacity across the five Picea species. Pnmax and PNUE are closedly associated with the fraction of N allocated to photosynthetic apparatus (Pphoto) compared with leaf N content per area (Narea). The tremendous growth differences among the five Picea species were substantially related to the interspecies variation in Pnmax and PNUE. Supplementary Information The online version contains supplementary material available at 10.1186/s12870-022-03593-x.


Introduction
Picea A. Dietrich (spruce) is the third largest genus in Pinaceae and consists of approximately 40 species mainly distributed in temperate regions and middle and high mountainous area of northern subtropical hemisphere [1,2]. The largest center of distribution and differentiation of Picea is in Asia, with the largest number of species occurring in China [1], including 16 species and 9 variants. It plays a substantial role in the boreal and subalpine forests in the northern hemisphere and features high economical value due to its high-quality timber [3]. Spruce introductions have started in China since the 1980s, mainly including early evaluation and species or provenance selection [4][5][6]. Interestingly, in these trials, some exotic (non-native) species exhibit strong adaptability and good growth performance in some regions of China such as Gansu and Jilin province [3,[5][6][7]. Similar trends were also revealed in the report of Holst [8], in which the growth of the introduced species, P. abies, could surpass native spruce species in North America. In the light of the important role Picea species play in ecology and economy, an enhanced knowledge of the physiological mechanisms underlying the interspecific growth differences may increase our ability to develop proper management plans for spruce introduction and breeding.
Plant survival and development are based on physiological processes and many of these occur within leaves, where the ability to maintain hydraulics, conduct photosynthesis, regulate temperature, gather and harvest light etc. are dependent on leaf structures [9,10]. It is widely accepted that leaf anatomy and mesophyll properties can affect carbon assimilation and leaf thickness are strongly correlated with leaf area-based photosynthesis [11][12][13]. As the key organs of conifers [14], many studies have shown that needle structures may vary tremendously among and within species e.g., [7,9,15], and these variations may have important implications for shaping plant functioning [16,17]. Considering the multi-faceted roles needles play, it is expected that needle structures should strongly affect the physiological functions at the leaf level, among which photosynthesis is an essential process. Convincing evidence has suggested that improving photosynthesis can promote plant growth although it is controversial whether plant growth is mainly limited by photosynthesis, i.e., the 'source' , or by demand of carbohydrates from new vegetative growth or the developing seeds and fruits, i.e., the 'sink' [18,19]. The most accepted view at present is that photosynthesis is a major driver of plant growth and productivity as it is the origin of carbon skeletons for plants to build up all their structures for both maintenance and growth [19], thus making it one of the most studied processes in plant physiology. Amongst leaf structural aspects, leaf mass per unit area (LMA) has been considered a vital functional trait due to its negative effects on the mass-based photosynthesis [20,21], but how LMA affects area-based photosynthesis (P n ) remains less clear when considered globally [22]. Moreover, although LMA has been demonstrated to be negatively correlated with mesophyll conductance (g m ) [23,24], an essential limiting factor for P n , this is not always true because LMA is a highly integrated trait. Different anatomical changes may lead to the uncertain relationship between LMA and g m [24,25]. Hence, whether leaf structures impose strong constraints on leaf functions remains empirical. From a functional perspective, needle internal structures should be directly related to physiological processes [26,27] and may have important implications for foliage potential photosynthesis [28][29][30]. However, little research has been conducted on the correlations between needle traits and photosynthetic characteristics at interspecies level in genus Picea, which may be essential for understanding the variation in growth rates and adaptive capacities across species.
Nitrogen (N) is one of the most important mineral nutrients that limit the growth of plants in many parts of the world [31]. The economics of N use in photosynthesis has become a research hotspot in the past several decades and has continued to this day [32,33]. As much as 75% of leaf N is invested into photosynthetic apparatus, with an average 20% invested in ribulose 1,5-bisphosphate carboxylase (Rubisco) [31,[34][35][36][37]. Strong positive correlation relationship between light-saturated net CO 2 assimilation rate (P nmax ) and leaf nitrogen content per unit area (N area ) observed in many species previously [38,39] indicate that leaf photosynthesis may be controlled by the supply and demand of leaf N content to a large extent. Photosynthetic-nitrogen use efficiency (PNUE), defined as the ratio of P nmax to N area , has been regarded as an important leaf trait for characterizing plant leaf economics, physiology, and survival strategy [40]. Species with high PNUE tend to have higher growth rates [40] and higher competitive ability in natural ecosystems [41].
Additionally, interspecific variation of P nmax and PNUE is related to N allocation in the photosynthetic apparatus (P photo , composed of leaf N allocated to Rubisco (P C ), bioenergetics (P B ) and light-harvesting components (P L )) [42,43]. Although genus Picea has been transplanted worldwide due to its huge economic and ecological values [7], to our knowledge, little information is available on how Picea species differ in photosynthetic N allocation, P nmax and PNUE, which may be directly related to the growth discrepancies among species.
In the present study, we employed five spruce species as our tested materials, out of which two were from North America (P. glauca and P. mariana), one from Heilongjiang province in China (P. koraiensis), one from Yunnan province in China (P. likiangensis), and one native to our experimental site (P. crassifolia), and cultivated them in the National Germplasm Repository Germplasm Resources Preservation of spruce in Xiaolong Mountain of Gansu province of China. Our previous research has revealed the significant growth differences among these species [3,7]. We hypothesize that the five Picea species vary widely in needle structures, biochemical parameters and photosynthetic-related traits and there is co-variation among these traits. In addition, the interspecies growth differences are closely related to the variation in photosynthetic capacity and PNUE among the five Picea species. The objectives of the present study were to (1) reveal the variation in needle structures, photosynthetic characteristics, photosynthetic N allocation and growth traits; (2) clarify the relationship among needle structures, biochemical parameters and photosynthetic-related traits; and (3) evaluate the consequences of these relationship on plant growth.
Pearson correlation analysis were conducted based on needle structural traits with significant difference among species. NT, NSA, MA, TMA and ETA were significantly positively correlated with NL. NW and NT both showed significantly positive correlations with NSA, MA, CCA and LMA while showed negative correlations with SD. In addition, LMA also showed strong positive correlations with CCA. NSA, MA and CCA were positively correlated with each other while NSA, MA and ETA were negatively correlated with CF ( Fig. 1). Table 1 Needle structures for the five studied Picea species NL needle length, NW needle width, NT needle thickness, SD stomatal density, NSA needle section area, ETA and ETF epidermis tissue area and fraction, MA and MF mesophyll area and fraction, CCA and CCF central cylinder area and fraction, TMA and TMF total mechanical tissue area and fraction. Data are mean ± SE. Different letters indicate significant differences among species (P < 0.05). *, statistical significance at P < 0.05; **, statistical significance at P < 0.01

Needle chemical traits
The five tested species showed significant differences in needle chemical traits including N area (g m − 2 ) and NSC (g m − 2 ). P. koraiensis exhibited the largest N area (6.61 g m − 2 ) while P. glauca and P. mariana exhibited the lowest (2.30 and 2.67 g m − 2 , respectivly). P.

Photosynthetic and biochemical parameters
Picea species vary considerably in photosynthetic (P nmax and iWUE) and biochemical parameters (g s , g m , V cmax and J max ) ( Table 2). Specifically, P. glauca showed the largest P nmax (9.27 umol m − 2 s − 1 ) while P. mariana and P. crassifolia showed the lowest (5.57 and 5.07 umol m − 2 s − 1 , respectively). P. likiangensis exhibited the largest iWUE (122.78 μmol mmol − 1 ) while iWUE of the remaining four species were not statistically different. V cmax , J max and g m of P. glauca, P. likiangensis and P. koraiensis were significantly higher than those of P. mariana and P. crassifolia (Table 2). To summarize, P. glauca, P. likiangensis, and P. koraiensis exhibited greater photosynthetic capacities in the five tested species, that is, higher P nmax , V cmax , and J max values. The correlation relationship between photosynthetic and biochemical parameters showed that increased iWUE was associated with lower gs while not correlated with g m , V cmax and J max . In contrast, P nmax was significantly positively correlated with g s , g m , V cmax and J max and the correlation between P nmax and g m (R 2 = 0.671, P < 0.01) was stronger than that between P nmax and g s , V cmax as well as J max (R 2 = 0.499, P = 0.000;

Fig. 2
Comparison of N area (leaf N content per area) and NSC (non-structural carbohydrate content) among species and linear relationship between LMA and N area as well as NSC. N area , leaf N content per area; NSC, non-structural carbohydrate content; LMA, leaf mass per area. Different letters indicate significant differences among species (P < 0.05). *, statistical significance at P < 0.05; **, statistical significance at P < 0.01 Table 2 Photosynthetic characteristics for the five studied Picea species P nmax , area-based maximum net photosynthesis; g s , stomatal conductance; iWUE, intrinsic water use efficiency; g m , mesophyll conductance to CO 2 ; V cmax and J max , maximum carboxylation and maximum electron transport rate, respectively. Data are mean ± SE. Different letters indicate significant differences among species (P < 0.05); *, statistical significance at P < 0.05; **, statistical significance at P < 0.01 R 2 = 0.549, P = 0.000; R 2 = 0.575, P = 0.000) (Fig. 3). P nmax showed no correlation relationship with iWUE. We used log values to evaluate the relationship between iWUE and g m /g s as well as V cmax /g s and the results showed that both were significantly positively correlated with iWUE (R 2 = 0.699, P = 0.000; R 2 = 0.648, P = 0.000) (Fig. 4).

Photosynthetic N allocation and PNUE
The N allocation in photosynthetic apparatus exhibited different modes among the five species. In general, P C , P B , P L , P photo , P non-photo and PNUE were all significantly different among species and P C was larger than P B and P L in all species (Table 3). Among the five species, P. glauca Fig. 3 Relationships between biochemical parameters and photosynthetic capacity across species. P nmax , area-based maximum net photosynthesis; g s , stomatal conductance; iWUE, intrinsic water use efficiency; g m , mesophyll conductance to CO 2 ; V cmax and J max , maximum carboxylation and maximum electron transport rate, respectively. Points represent trait values for each individual of each species. Explained variance (R 2 ) and P values are shown. Solid lines in each panel represents significant linear regression

Fig. 4
Relationships between log (g m /g s ) as well as log (V cmax /g s ) and iWUE and the relationship between log (g m /g s ) and log (V cmax /g s ). g s , stomatal conductance; g m , mesophyll conductance to CO 2 ; V cmax and J max , maximum carboxylation and maximum electron transport rate, respectively. Points represent trait values for each individual of each species. Explained variance (R 2 ) and P values are shown. Solid lines in each panel represents significant linear regression had the largest P L , P B , P C and P photo (2.74, 3.71, 11.96 and 18.41%, respectively) while P. crassifolia had the lowest P L and P photo (0.74 and 3.23%, respectively). P B , P L , P photo and P non-photo of P. mariana, P. likiangensis, P. koraiensis were not significantly different statistically. P. glauca exhibited the largest PNUE (56.89 μmol mol − 1 s − 1 ) while P. koraiensis and P. crassifolia had the lowest among species (15.93 and 16.22 μmol mol − 1 s − 1 , respectively) ( Table 3).
The results of linear regression analysis showed that both P nmax and PNUE were significantly positively correlated with P L , P B , P C and P photo , but the correlation relationship between P L , P B , P C as well as P photo and PNUE (R 2 = 0.625, P = 0.000; R 2 = 0.792, P = 0.000; R 2 = 0.710, P = 0.000; R 2 = 0.759, P = 0.000) were stronger than that between P L , P B , P C as well as P photo and P nmax (R 2 = 0.431, P = 0.002; R 2 = 0.518, P = 0.000; R 2 = 0.497, P = 0.000; R 2 = 0.522, P = 0.000) (Fig. 5).

Species differences in growth traits
The results of ANOVA showed that tree height at age 6, 9 and 11, DBH, and the number of first lateral branches were all significantly different among species, and Table 3 Photosynthetic N allocation and PNUE for the five studied Picea species P C , P B , P L , P photo and P non-photo , the fraction of leaf N allocated to Rubisco, bioenergetics, the light-harvesting components, photosynthetic apparatus and nonphotosynthetic apparatus (%), respectively. PNUE, photosynthetic nitrogen use efficiency. P photo = P L + P B + P C ; P non-photo = 1-P photo . Data are mean ± SE. Different letters indicate significant differences among species (P < 0.05); *, statistical significance at P < 0.05; **, statistical significance at P < 0.01  5 Relationships between the fraction of N allocated to different photosynthetic apparatus and P nmax and PNUE. P C , P B , P L and P photo, the fraction of leaf N allocated to Rubisco, bioenergetics, the light-harvesting components and photosynthetic apparatus (%), respectively; P nmax , area-based maximum net photosynthesis; PNUE, photosynthetic nitrogen use efficiency. Points represent trait values for each individual of each species. Explained variance (R 2 ) and P values are shown. Solid lines in each panel represents significant linear regression tree height at age 11 was significantly different among blocks (Table 4). Moreover, the interactions of species and blocks were significant for height (11a) and DBH. P. glauca showed the largest tree height at age 6, 9 and 11, DBH and the number of first lateral branches among the five species while these traits were all the lowest for P. crassifolia. To summarize, at our experimental site, P. glauca exhibited the best growth performance among species, with P. likiangensis and P. koraiensis the middle and P. mariana as well as P. crassifolia the lowest (Fig. 6). Additionally, significantly positive linear relationship was observed between P nmax as well as PNUE and growth traits. Thus, the growth discrepancies among the five studied species were primarily driven by P nmax and PNUE. However, the correlation relationship between height (11a) as well as DBH and P nmax (R 2 = 0.794, P = 0.000; R 2 = 0.754, P = 0.000) was stronger than that between the two growth traits and PNUE (R2 = 0.548, P = 0.000; R 2 = 0.475, P = 0.001), whereas both P nmax and PNUE were tightly correlated with the number of first lateral branches (R 2 = 0.731, P = 0.000; R 2 = 0.731, P = 0.000) (Fig. 7).

Discussion
Our results clearly showed the interspecies variation for needle structures, chemical components, photosynthetic capacities, biochemical parameters, photosynthetic N allocation and growth traits as well as the coordination or trade-offs between them. To our knowledge, this work represents the first relatively comprehensive study trying  . 6 Comparison of growth traits among species. Height (6a), Height (9a) and Height (11a), tree height at age 6, 9, and 11, respectively; DBH, diameter at breast height. Different letters indicate significant differences among species (P < 0.05). *, statistical significance at P < 0.05; **, statistical significance at P < 0.01 to figure out all these relationships together among multiple species within genus Picea. The variations in these characteristics and the correlations between them may make it possible for tree breeders to select the most suitable species for a given environment, for example, our experimental site with a temperate subhumid climate. We assume that the significant differences observed in needle structural characteristics were mainly attributed to genetic factors as needle structures are relatively stable at species level [44]. In addition, the significant correlation relationship between some traits (Fig. 1) may indicate the constraints of needle structural design. Leaf size affects many functional traits and is related to plant productivity and adaptive capacity to changes in its environment [45], but this has been seldom studied in genus Picea. Considering the great differences of pine needle length among species, Wang et al. reported that NL was positively correlated with mechanical tissue fraction and area while negatively correlated with MA and MF in five southern pine species, suggesting that increasing needle stiffness in these species came at a cost of photosynthetic tissue (mesophyll) when NL increased [9], but this was on the premise that NSA did not increase with NL. Similar trade-offs based on findings on tissue fractions in three Mediterranean pines were also reported by Kuusk et al. [46]. Although MA, ETA and TMA were found to be positively correlated with NL in our results, but the variation of them along needle length across species was also significantly integrated by NW or NT, and CCA was only affected by NW and NT regardless of NL. Moreover, NW was also detected to be significantly correlated with relative indicators including MF and TMF. These results may indicate that, for the studied Picea species, variations in needle internal structures should be largely, at least partially, dependent on NW and NT despite the fact that NL do have a certain impact on some traits. Moreover, LMA was also found to be positively related to NW and NT. In other words, needle anatomical changes are related to morphological adjustments [9] and the resulting coordination and trade-offs between traits may be the consequence of the interaction of various indicators regarding needle size. Needles must be rigid and flexible enough to support their own weight as well as bend or twist without breaking when subjected to external forces [47,48], which could be reflected by increasing TMA when NL increased in this study. However, Picea species show a much narrower range of NL (0.82-1.15 cm in this study) in Pinaceae comparing with genus Pinus (e.g., 7.51-34.59 cm in [9]), thus it's unlikely for Picea species to invest obviously more to needle supporting tissues at the cost of photosynthetic tissue (mesophyll) when NL increase. Additionally, our results also showed that NL, NW and NT positively affected NSA unevenly while NSA simultaneously affected internal anatomical structures. Different needle tissue area and MF were positively while TMF and CCF were negatively correlated with NSA. Interestingly, although both MA and TMA increased upon increasing NSA, MF showed a significantly negative correlation relationship with TMF. This may suggest a trade-off between photosynthetic tissue and supporting tissue when NSA increased and, seen from this study, this trade-off was mainly driven by NW. Our results add to the previous findings on the relationship between leaf size and investment to supporting tissues [9,49,50], suggesting that, in the studied species with narrow NL range, NW may also play a vital role in affecting needle anatomical traits and the resulting tradeoff between photosynthetic tissue and supporting tissue may be related to surface exposure to incoming light. In addition, the correlation relationship among anatomical traits indicate that various needle inner structures may be developmentally interrelated [7,44].
Biologists seek to distill organismal variation down to a manageable number of key functional traits [22,51], within which leaf traits may fall along a 'fast-slow' continuum [52]. In other words, variations in needle traits, for conifers, may impose strong constraints on physiological functions and may furtherly affect plant survival and development. From a functional perspective, needle anatomical characteristics may be more directly related to physiological processes [26,27]. For example, larger mesophyll tissue or vascular bundles may lead to higher photosynthetic capacity [14] or material transport capacity [53]. Additionally, LMA, one of the key leaf traits, has been demonstrated to negatively affect photosynthesis by reducing CO 2 diffusion in many studies. These assumptions have been studied at both large and small scales [22][23][24], but seldom addressed in genus Picea, especially in species introduction trials when evaluating plant adaptability to new environments. In one of our previous studies, we found that net photosynthetic rate was positively correlated with some needle anatomical traits such as MA and CCA across 17 Picea species (including the five species in this study). However, the effects of needle internal anatomical structures as well as needle biochemical parameters on the differences of potential photosynthetic capacity among species remained unknown to a large extent [7]. Here, a more comprehensive study was conducted aiming at figuring all these relationships out. We found that the variation in the inherent photosynthetic capacity among the studied species were not related to the variation in LMA and needle anatomy (Table S1) but resulted mainly from the differences in g s , g m and biochemical parameters including V cmax and J max in this study (Fig. 3). Moreover, the results of Grey correlative degree analysis (GRA) showed that g m , V cmax and J max contributed more to the interspecies variation in P nmax (Fig. S1). g s , g m and biochemical parameters including V cmax and J max were not correlated with needle anatomy and LMA, either (Table S1). Interestingly, we found that increased NSC (positively correlated with LMA, Fig. 2) was associated with decreased g m among species, which has been seldom reported previously. NSC produced by photosynthesis is essential in plant functioning, being the substrates and energy sources for metabolic processes [54,55]. Our results indicate that the accumulation of NSC may lower mesophyll CO 2 diffusion across species, but this mechanism requires more research. In summary, two conclusions may be drawn from these results. Firstly, tissue-level anatomy is not sufficient to explain the variations in photosynthetic capacity among species. Although tissue-level anatomical traits are highly expected to be related to physiological functions as they are directly involved in physiological processes, among which photosynthetic capacity may be more related to cell-level anatomical traits such as mesophyll and chloroplast surface area exposed to intercellular air space [22,[56][57][58] and the contribution of different ultracellular and cellular components [22]. Secondly, for the studied species in genus Picea, robust leaf structure did not impose any constraints on g m as well as photosynthetic capacity. Increased LMA is expected to be related to lower g m [24,28], but there are also exceptions; for example, both LMA and g m would increase upon increased mesophyll cell surface area. This result was consistent with the findings in seven Mediterranean oaks species reported by Peguero-Pina et al. [22] and the findings in wild and domesticated cotton reported by Lei et al. [58], which may reflect the multiple determinants of LMA. Therefore, further study should be aimed at evaluating the relationship between needle finer-scale anatomy and LMA and, more importantly, better figuring out the coordination or trade-offs among needle functional traits and physiological processes underlying plant growth and adaptability in Picea breeding programs.
We found that increased g s , g m , V cmax and J max could result in higher P nmax . Additionally, iWUE scaled with both g m /g s and V cmax /g s (Fig. 3), suggesting that the two parameters can be considered as potential targets to improve iWUE [59,60]. When considered under a large sample size including herbaceous, shrub and tree species, the two ratios were poorly correlated with each other and the absence of correlation indicated that improving each of them separately would improve iWUE [60]. However, g m /g s was tightly correlated with V cmax /g s in this study (Fig. 4), which implied that an enhanced iWUE required improving both g m /g s and V cmax /g s simultaneously. The discrepancies between this work and previous studies may be attributed to the limited samples involved in this study. Generally, increased g s will simultaneously increase transpiration rate and decrease water use efficiency [24,61]. Therefore, if g s was tightly correlated with g m and V cmax , it may be challenging to increase water-use efficiency without sacrificing photosynthetic rate. Although increasing g m or V cmax alone would result in increasing photosynthesis, and potentially increasing WUE, in practice, higher iWUE seems to be only achieved when there are no parallel changes in g s [59]. In this study, g s was weakly correlated with g m and not correlated with V cmax (Fig. S2, Fig. 3). In addition, no significant correlation relationship was detected between iWUE and P nmax across species (Fig. S2). These results made it possible to identify species with both relatively high photosynthetic rate and water use efficiency-P. likiangensis showed the second highest P nmax and the highest iWUE among the five studied species. For the other four species, they exhibited contrasting P nmax but similar and low iWUE values ( Table 2). Here comes the question-why don't all species have high photosynthetic rate while maintaining high water-use efficiency? Muir et al. [24] proposed that other trade-offs in nature like water or N limitation can also have a negative effect on g m or biochemical parameters, which would furtherly affect photosynthetic capacity. For example, a habitat with much water but little N may be beneficial to leaves with high stomatal density but low biochemical investment [61,62] without altering the selection on overall leaf structure. These results suggest that the adjustments of water and N use efficiencies (PNUE and WUE) may also play an important role in determining the adaptability of plant species as the coordination or trade-offs along leaf spectrum do [24], even though no obvious coordination or trade-offs between needle structures and physiological processes were detected in this study (Table S1). P nmax , PNUE and WUE progressively increases along the phylogenetic tree of land plants [19], but when comparing different species or subjecting a species to e.g. water stress, both P nmax and PNUE decrease concomitantly with increased WUE [63,64]. This negative correlation relationship was also found in invasive vs. native species [65,66] and transgenic Btcotton vs. its conventional peer [67]. Although PNUE was positively correlated with P nmax in this study, it did not show any correlation with iWUE (Fig. S3). However, all studied species showed low PNUE but high iWUE, which is consistent with the results reported by Flexas and Carriquí [19], in which gymnosperms presented the largest WUE but relatively low PNUE through the phylogeny and they were not significantly correlated with each other.
Many studies have demonstrated that P nmax increases with N area [68], but this relationship is also dependent on species [69]. For example, an exception that P nmax was inversely correlated with N area in five southern pine species native to southeastern America was reported by Wang et al. [9]. In this study, we found that variations in P nmax across species was not related to N area (Table S1) but positively correlated with P photo (Fig. 5). These results suggest that it is the fraction of N allocated to photosynthetic apparatus rather than needle total N content that is the key adaptive mechanism to maximize photosynthesis [70,71]. In addition, both P nmax and PNUE were significantly positively correlated with not only P nmax but also P L , P B and P C . This result is consistent with the findings reported in some common tress species in subtropical China [71,72]. In this study, P photo , P L , P B and P C varied tremendously among species, with P C being the largest in them for each species (Table 3), which has been well documented in previous studies (e.g., [33]). However, P photo (3.23-18.41%) and P C (1.64-11.96%) showed much lower values than those reported in previous studies; for example, P photo and P C of Fagaceae and Leguminosae tree species ranged between 13.86 and 27% and between 9.8 and 20%, respectively [73]. The differences between this work and previous studies may be attributed to species differences--Picea species may invest more N to protective and storage chemicals [9] than other species.
In this study, the five species exhibited significant differences in growth traits at the experimental site. P. glauca exhibited the best growth performance, with the largest tree height, DBH and the number of first lateral branches among the five species at age 11. P. likiangensis and P. koraiensis also exhibited good growth performance while P. crassifolia and P. mariana showed poor growth ( Table 4, Fig. 6). The considerable growth differences among the five species were largely determined by the variations in P nmax and PNUE across species. In addition, the growth period for P. glauca was 95 days at the experimental site and was longer than that of the other four species [3], which may also be associated with the best growth of P. glauca among the five species. Interestingly, the two species introduced from North America, P. glauca and P. mariana, showed diametrically opposite growth characteristics. P. glauca and P. mariana have been introduced into China for more than 20 years. Chen et al. [74] first reported that the growth of P. glauca and P. mariana planted in northeast China could surpass the native species, P. koraiensis, at age 4, 5 and 6, but no differences were observed between the two species. Similar findings with this study that the growth of P. glauca was superior than P. mariana and P. crassifolia grown in Gansu province were also reported by An et al. [75]. P. glauca and P. mariana originated from North America, and most of their allocation ranges overlap with each other. However, P. glauca mostly grows in nutrient-rich soil, while P. mariana mostly grows in soil with relatively poor nutrition [76]. Therefore, the contrasting growth performance of the two species may be related to the long-term adaptation to the environments of their original habitats.
As a final note, we wish to point out that, our study was carried out in a common garden, so current variation in the measured traits among species might be a combination of genetic differences and the adaptive nature in response to both historical and current environmental conditions at the experimental site [77]. Additionally, photosynthetic parameters were estimated at the leaflevel, while the actual consequences of various structural designs on carbon gain and water use must ultimately be evaluated from shoots (the functional unit) to crowns and canopies [9], which requires further work.

Conclusions
To the best of our knowledge, this study is a comprehensive comparison of needle structures, biochemical parameters, photosynthetic-related traits and plant growth among multiple Picea species and we found that several relationships between them can be described at the species level. To address the hypothesis posed in the Introduction, we draw four main conclusions: (1) changes in needle anatomical structures across the five studied species are affected by needle size, especially by NW; (2) compared with tissue-level anatomy and LMA, needle biochemical parameters play a more important role in determining the interspecific variation in photosynthetic capacity; (3) P nmax and PNUE were affected by the fraction of N allocated to photosynthetic apparatus rather than needle total N content; (4) the five studied Picea species vary widely in growth traits, and this variation is substantially related to the differences in P nmax and PNUE among the species. These insights provide a better understanding of the existing relationship between different biophysical structures in genus Picea and will be beneficial for the identification and selection of Picea species for a specific site condition.

Study site and field experiment
In 2008, seeds of the five Picea species were sown in seedbeds and after cultivation for 3 years, the seedlings were planted with 1.5 m × 1.5 m spacing at the same site in a nursery at the Research Institute of Forestry of Xiaolong Mountain in Gansu Province, China (34°28′50″ N, 105°54′37″ E). Details on filed design, climate type at the experimental site, and geographical information as well as climate data for the 5 species have been reported elsewhere [3,7].

Growth trait measurements
In July 2019, plant height, diameter at breast height (DBH) and the number of first lateral branches (at age 11) of the five species were determined. Additionally, tree height at age 6 and 9 was also integrated in this study for analysis.

Needle gas exchange measurements
Gas exchange measurements were all conducted in late July 2019. Light-response curves of current-year-old needles on intact upper-crown branches were measured using Li-Cor 6400 gas-exchange system (with red/blue light source, and 6400-22 L conifer chamber; Li-Cor Biosciences, Lincoln, NE, USA). Four trees of each species were randomly selected, and one set of measurements were taken from each tree (n = 4). Photosynthetic response to photosynthetic photon flux density (PPFD) and C i were determined on the same one healthy shoot. Under 400 μmol·mol − 1 of leaf chamber CO 2 concentration, photosynthetic rates were measured at photon flux densities of 1800, 1600, 1200, 1000, 800, 600, 400, 200, 100, 75, 50, 25, 0 μmol·m − 2 ·s − 1 after each leaf sample was illuminated with a saturating level of PPFD for 10-30 min to achieve fully photosynthetic induction. Photosynthetic responses to variations in leaf internal CO 2 concentration, i.e. A-Ci curves, were measured on the same shoot of each species where lightresponse curves were determined (n = 4). For each sample, a single curve over 16 chamber CO 2 concentrations (400 、300、200、100、60、400、400、600、800、1000 、1200、1400、1600、2000、2300、2600 μmol·mol − 1 ) was measured. Maximum carboxylation rate (V cmax ), and maximum electron transport rate (J max ) were estimated from A-Ci curves using the plantecophys R package [78,79]. In addition, the data of A-Ci curves were uploaded in the LeafWeb server (http:// www. leafw eb. org/) in order to have an automated analysis for mesophyll conductance (g m ) [73,80].

Determination of LMA (leaf mass per unit area), chlorophyll and N content
Following gas exchange measurements, needles were carefully harvested and then pictured for needle area calculation using ImageJ software. Afterwards, needles were weighed and frozen in liquid N 2 . For leaf dry mass measurements, needles near the shoots for gas exchange measurements were collected, pictured and weighed as described above. After being oven-dried at 65 °C for 48 h and weighed, leaf mass per area (LMA, g m − 2 ) was calculated as the ratio of needle dry mass(g) to needle area(m 2 ).
Dried needle samples were ground into fine powder in a ball mill, needle N concentration (N mass , mg/g) was determined by a by high-temperature combustion using a LECO elemental analyzer (Leco Corporation, St-Joseph, Michigan). N mass was then converted to a projected area basis (N area ) using the LMA measurements for each sample. PNUE was determined as the ratio of P nmax to N area .
Absolute chlorophyll concentration measurements were conducted using 95% (v/v) alcohol extracts of needle tissue and a Shimadzu visible-ultraviolet spectrophotometer (UV 2250, Fukuoka, Japan) following the procedure raised by Tang et al. [73].

Calculation of nitrogen allocation in the photosynthetic apparatus
Nitrogen allocation fractions of each component in the photosynthetic apparatus were calculated according to Niinemets and Tenhunen [81].
Where 6.25 is the N conversion coefficient into Rubisco (g Rubisco g − 1 leaf N) [82] and 8.068 is the N binding coefficient for cytochrome f [83]. C C is the chlorophyll concentration (mmol g − 1 ). V cr is the specific activity of Rubisco (20.78 umol CO 2 g − 1 Rubisco s − 1 ) at 25 °C, C B is the ratio of chlorophyll to organic N in light harvesting components (2.15 mmol g − 1 ) at 25 °C and J mc is the potential rate of photosynthetic electron transport (155.65 μmol electrons μmol − 1 Cyt f s − 1 ) at 25 °C [81]. P C , P B , and P L are the fraction of leaf nitrogen allocated to Rubisco, bioenergetics, and the light-harvesting components (%), respectively. P photo indicated the leaf N allocated to the photosynthetic apparatus (P Photo ) while P non-photo indicated the leaf nitrogen allocated to nonphotosynthetic apparatus.

Determination of non-structural carbohydrate content (NSC)
The concentrations of soluble sugar and starch were determined using the anthrone method as described by Li et al. [84]. NSC was calculated as the sum of the concentration of soluble sugar plus the starch concentration.

Statistical analysis
The normality of all data and the homogeneity of variance were tested prior to statistical analysis. An ANOVA to evaluate the species effect on growth traits of each year was carried out using the following model: where y ijk is the observed value of species i in block j; μ is the mean value of the species; S i is the random effect of species i = 1,2, …,5; B j is the fixed effect of block j = 1, 2, 3, 4; SB ij is the random effect of the interaction of species i and block j. e ijk is the random error.
The effect of species on needle gas exchange and needle structures were tested using one-way analysis of variance (ANOVA): P B = J max 8.068 × J mc × N area × 100% P C = V C max 6.25 × V cr × N area × 100% P photo = P L + P B + P C P non−photo = 1 − P photo (1) y ijk = µ + S i + B j + SB ij + e ijk (2) y ij = µ + S i + e ij