Next Article in Journal
Experimental Study on Fiber Extraction after Saturated Steam Softening of Pleioblastus amarus
Next Article in Special Issue
A Method for Genetic Transformation Using Embryonic Callus of Pinus koraiensis
Previous Article in Journal
Comparative Study on Vascular Bundle Morphological Characteristics of Parts of Branches, Culms, and Rhizomes of Oligostachyum sulcatum
Previous Article in Special Issue
Progeny Selection and Genetic Diversity in a Pinus taeda Clonal Seed Orchard
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Screening and Site Adaptability Evaluation of Qi-Nan Clones (Aquilaria sinensis) in Southern China

1
Research Institute of Tropical Forestry, Chinese Academy of Forestry, Guangzhou 510520, China
2
College of Landscape Architecture, Nanjing Forestry University, Nanjing 210037, China
3
State Key Laboratory of Efficient Production of Forest Resources, Beijing 100091, China
4
Tropical Forestry Experimental Centre, Chinese Academy of Forestry, Pingxiang 532600, China
5
Guangzhou Institute of Forestry and Landscape Architecture, Guangzhou 510080, China
*
Author to whom correspondence should be addressed.
Forests 2024, 15(10), 1753; https://doi.org/10.3390/f15101753
Submission received: 6 September 2024 / Revised: 29 September 2024 / Accepted: 3 October 2024 / Published: 5 October 2024
(This article belongs to the Special Issue Forest Tree Breeding, Testing, and Selection)

Abstract

:
In recent years, plantations of Aquilaria sinensis in China have been dominated by Qi-nan, yet there remains limited research on the growth evaluation and breeding of these clones. In this study, a multi-point joint variance analysis, an additive main effect and multiplicative interaction (AMMI) model, a weighted average of absolute scores (WAASB) stability index, and a genotype main effect plus a genotype-by-environment interaction (GGE) biplot were used to comprehensively analyze the yield, stability, and suitable environment of 25 3-year-old Qi-Nan clones from five sites in southern China. The results showed that all the growth traits exhibited significant differences in the clones, test sites, and interactions between the clones and test sites. The phenotypic variation coefficient (PCV) and genetic variation coefficient (GCV) of the clones’ growth traits at the different sites ranged from 16.56% to 32.09% and 5.24% to 27.06%, respectively, showing moderate variation. The medium–high repeatability (R) of tree height and ground diameter ranged from 0.50 to 0.96 and 0.69 to 0.98, respectively. Among the clones, Clones G04, G05, G10, G11 and G13 showed good growth performance and could be good candidates for breeding. Environmental effects were found to be the primary source of variation, with temperature and light primarily affecting growth, while rainfall influenced survival and preservation rates. Yangjiang (YJ) was found to be the most suitable experimental site for screening high-yield and stable clones across the different sites, whereas the tree height and ground diameter at the Chengmai (CM) site were significantly higher than at the other sites, and the Pingxiang (PX) and Zhangzhou (ZZ) sites showed poor growth performance. The findings suggest that Qi-nan clones are suitable for planting in southern China. There were also abundant genetic variations in germplasm resources for the Qi-nan clones. The five selected clones could be suitable for extensive planting. Therefore, large-scale testing is necessary for determining genetic improvements in Qi-nan clones, which will be conducive to the precise localization of their promotion areas.

1. Introduction

Agarwood is a type of resin-containing wood that forms after Aquilaria, Gonystylus, and Gyrinops species are injured [1]. It is frequently used as a precious medicine and natural spice [2], causing it to have extremely high medical and economic value. Therefore, many countries and regions have attached great importance to the cultivation and breeding of easily formed agarwood and high-yield agarwood varieties [3]. In China, researchers and farmers have been breeding for many years to develop a novel special germplasm of A. sinensis (Lour.) Spreng [4], which is cultivated by grafting from a specific wild A. sinensis population that possesses outstanding induced agarwood potential [5,6]. This special clone germplasm, called “Qi-nan” [4], has the advantage of being more likely to induce agarwood and high-quality agarwood [6,7,8]. Studies have shown that Qi-Nan can be distinguished from ordinary A. sinensis through SSR and the composition of agarwood [5,9]. After one year of induction, the alcohol extract content of Qi-Nan will exceed 40% [6,10], accompanied by a richer and more diverse range of sesquiterpenes [5].
In recent years, plantations of A. sinensis in China have been dominated by Qi-nan, which is mainly planted in Hainan, Guangdong, Guang xi, Fujian and Yunnan in the south of China [4]. Compared to ordinary A. sinensis, the biggest advantage of Qi-nan is its ability to maintain the excellent traits of mother trees [11]. However, with the increased demand for its planting, the species source of many Qi-Nan clone germplasms remains unclear, and they are widely planted without systematic breeding [4], resulting in significant differences between the clones, with coefficients of variation reaching 19% for tree height, 23% for ground diameter, and 17% for agarwood yield, as has been observed in a single test site [7]. Therefore, we carried out a trial in multiple environments to study the variation in the important traits of Qi-nan clones, which is of great significance for their successful planting.
Multi-environmental trials, as an important part of the commercial breeding process, are a prerequisite for varieties promotion [12]. In research, they are often used to evaluate the productivity, stability, and adaptability of species genotypes, as well as to explain the interaction effects between the genotype and the environment [13,14,15]. Additionally, they have also been applied to Larix olgensis [16], Pinus koraiensis [17], poplar clones [18,19,20], and European beech [21]. Through these multi-environmental trials, the growth traits and genetic variation patterns of the tested varieties (lines) have been analyzed [22], the interaction effect between genotypes and the environment has been explained [23], the causes of variation have been revealed, and improved varieties have been screened out [24]. This confirms the necessity and validity of multi-environmental trials in the breeding of improved varieties.
Generally, univariate, multivariate, mixed, and non-parametric methods are usually used to study the interactions between genotypes and the environment in multi-site identification trials [25]. Variance analysis and linear regression are not comprehensive and intuitive enough for analyzing the effects of genotype × environment interactions (GEI) [26]. The AMMI model and GGE biplot are effective statistical methods in multivariate analyses which are commonly used in the analysis of multi-environmental trials [27,28]. The AMMI model combines factor analysis with variance analysis to separate the model error and interference from the residuals of the additive model [29,30]; it draws a biplot to describe the yield and interaction effect of the variety and location, and it intuitively analyzes the stability of the varieties [31]. Compared to the AMMI model, the GGE biplot is a principal component analysis of environmental centralization, and it has a strong regional genotype evaluation function [32]. It can be used to combine the main effect of the genotype with the interaction of the GE at the same time, cluster the environment, and divide the suitable area of the varieties [18,33]. Moreover, a new variety evaluation system has been introduced: the WAASB (Weighted Average of Absolute Scores) stability index, which is a mixed model incorporating a best linear unbiased prediction (BLUP) framework that takes into account all the IPCAs for stability analysis [31,34]. Through the WAASB stability index and the AMMI model, the stability of various crops and the discrimination of the experimental environment were reasonably explained [13,15], which effectively shortened the breeding process [35,36]. Combining the AMMI model with the GGE biplot can help in better understanding the effects of GEI identifying ideal test conditions, and ensuring cultivar adaptation regions [16].
By the end of 2023, it was roughly estimated that there were more than 100 variates of Qi-nan clones in China, with a planting area of more than 35,000 hectares. However, there is a lack of excellent Qi-nan clones adapted to different regions, which hinders the development of the agarwood industry. In this study, the growth performances of 25 Qi-nan clones from five sites were investigated, the effect of the GEI was examined via variance analysis, the AMMI model, the WAASB stability index, and the GGE biplot to understand the process of tree breeding. Our results can provide a scientific basis for the selection and evaluation of Qi-nan clones in southern China.

2. Materials and Methods

2.1. Study Area and Materials

In November 2020, scions from twenty-five Qi-nan clones (G01 to G25) were collected from the first generation of grafted mother trees in Maoming, Huizhou, Shenzhen, and Shanwei, and were grafted onto seedlings of A. sinensis as rootstocks in Maoming, Guangdong Province. After six months of growth, the ground diameters and tree heights of the grafted seedlings ranged from 0.8 cm to 1.3 cm and from 1.0 m to 1.3 m. In April 2021, these seedlings were transported to various test sites, and a randomized block design with 3 blocks and 15–20 replicates was implemented across five test sites (Figure 1 and Table 1), with a row spacing of 2.0 m × 1.5 m. Consistent management was applied to each site.

2.2. Measurement of Growth Characters

In June 2022, the number of surviving trees of all the sample trees in the experimental plot was investigated. In December 2023, the growth traits, such as the preserved ratio (%), ground diameter (cm), and tree height (m) of all the sample trees in each experimental plot were measured.

2.3. Statistical Analysis

In this study, data entry was performed using Excel 2021. An analysis of variance (ANOVA), Duncan’s multiple-comparison analysis, and an estimation analysis of genetic parameters were performed using SPSS 28.0 (International Business Machines—IBM Corporation, Chicago, IL, USA).
A linear mixed model for a single site analysis of variance is as follows [37]:
Y j k l   = μ + B j + G k + B G j k + ε j k l
where B j is the block effect (j = 1, 2, 3), G k is the clone effect (k = 1, 2, 3, 25), B G j k is the interaction effect between kth clone and jth block, and ε j k l is the random error.
A variance linear model for multi-site joint analysis is as follows [38]:
Y i j k l = μ + S i + B i j + G k + S G i k + B G i j k + ε i j k l
where μ is the overall mean, S i is the site effect (i = 1, 2, 3, 4, 5), B i j is the block effect, G k is the clone effect, S G i k is the interaction effect of kth clones at ith site, B G i j k is the interaction effect between kth clones and jth blocks in ith site, and ε i j k l is the random error.

2.4. Estimate the Genetic Parameters

The calculation formulas of the phenotypic and genotypic coefficient of variation (PCV and GCV) are as follows [39]:
P C V % = δ P 2 X ¯ × 100
G C V % = δ g 2 X ¯ × 100
where δ P 2 and δ g 2 are the phenotypic variance and genetic variance of a trait, and X ¯ is the phenotypic mean of the trait.
The formula of clone repeatability (R) is as follows [16,19]:
R = σ A 2 +   σ b 2 σ A 2 +   σ b 2 +   σ e 2 = 1 1 F
where σ A 2 , σ b 2 and σ e 2 are the genetic variance of clones, the interaction variance between clones and environment, and the environmental random error variance, respectively. F is the F value of variance analysis.

2.5. The AMMI Model, WAASB Stability Index, and GGE Biplot

The “Argricolae” package, “metan” package, and “GGE Biplot GUI” package in the R 4.2.1 software were used to analyze the AMMI model, WAASB stability index, and GGE biplot, respectively, as well as to conduct correlation analysis.
The AMMI model is as follows [40]:
Y i k l = μ + S i + G k + r = 1 n r Ψ k r ɓ i r + P k r + ε i k l
where μ is the overall mean, S i is the average deviation of the environment, G k is the average deviation of clones, r is the eigenvalue of the r principal component analysis, and Ψ k r and ɓ i r are the clone and environmental scores of the r principal component analysis, respectively. n is the total number of principal component factor axes in model analysis, P k r is the interactive remaining items, and ε i k l is the residual error.
The GGE biplot is as follows [41]:
Y i j μ B j = 1 r i 1 δ j 1 + 2 r i 2 δ j 2 + B G i j k + ε i j
where μ is the overall mean, and B j is the mean value of all clones in site j. 1 and 2 are the eigenvalues of the first and second principal components, respectively. r i 1 and r i 2 are the eigenvectors of the first and second principal components of the i clone, respectively. δ j 1 and δ j 2 are the eigenvectors of the first and second principal components at the j location, respectively. ε i j is the residual error.

3. Results

3.1. Analysis of Variance for Growth Traits

The effects of the clones, test sites, and interactions between the clones and test sites on the tree height and the ground diameter were extremely significant across multiple environments (p < 0.01) (Table 2). The proportion of the sum of squares to the sum of squares of the total variation of the clone effect, the test sites effect, and the interactions between the clones and test sites for tree height were 19.7%, 63.3%, and 17.0%, respectively. The proportion of each effect on the ground diameter to the sum of squares of the total variation were 35.1%, 47.8%, and 17.1%, respectively. This indicated that the environmental effect was the main source of variation, followed by the clone effect. In each test site, there were significant differences in the tree height and the ground diameter between the different clones (p < 0.01). It was possible to select a variety of excellent Qi-nan clones as there were abundant variations in the Qi-nan germplasm.

3.2. Growth Performance of Different Clones and in Different Test Sites

The average tree height, ground diameter, 1-year average survival rate, and 3-year average preservation rate were 2.44 m, 5.18 cm, 86.03%, and 72.62%, respectively (Table 3). The tree height and ground diameter of G11, which exceeded the average by 14.89% and 29.65%, respectively, were found to be the largest. Conversely, G18 was found to have the smallest averages with respect to tree height and ground diameter, which were 22.95% and 25.19% smaller than the average values, respectively. The average tree height, ground diameter, 1-year average survival rate, and 3-year average preservation rate were the highest on the CM site, with values of 3.23 m, 6.57 cm, 94.42%, and 90.82%, respectively, which are significantly higher than those of other sites.

3.3. Estimating the Genetic Parameters

The PCV of the tree height and ground diameter ranged from 16.56% to 32.09% and 20.48% to 24.58%, respectively (Table 4). The GCVs of the tree height and the ground diameter ranged from 5.24% to 27.06% and 8.51% to 20.66%, respectively. The PCVs of the clone tree height and the ground diameter at each test site were greater than the GCVs, indicating abundant variations in the growth traits of the Qi-nan clones. The PCV and GCV values of the tree height at the PX site showed the highest values, 32.09% and 27.06%, respectively. Moreover, the GCVs of the tree height and the ground diameter on the ZZ site had the lowest values, 5.24% and 8.51%, respectively. The repeatability of the tree height and the ground diameter of the clones at different test sites ranged from 0.50 to 0.96 and 0.69 to 0.98, respectively. Except for the ZZ site, the repeatability of the height and the ground diameter of the clones at the other four sites was greater than 0.8. The ground diameter repeatability of the different sites was greater than the tree height except for at the FS site, indicating that the ground diameter was less affected by the environment than the tree height.

3.4. Correlation Analysis between Growth Traits and Environmental Factors

The tree height and the ground diameter had a highly significant positive correlation with the mean annual temperature (Table 5). Additionally, they also had a significant positive correlation with the mean annual sunshine duration. The 1-year survival rate had a significant positive correlation with the mean annual precipitation, and the mean annual sunshine duration. Similarly, the 3-year preservation rate had a highly significant positive correlation with the mean annual sunshine duration and a significant positive correlation with the mean annual precipitation.

3.5. AMMI Model Analysis

The AMMI biplots effectively explained the interaction information of the tree height and ground diameter (Figure 2). The AMMI biplots based on tree height and ground diameter explained 91% (PC1: 80.2%, PC2: 10.8%) and 78.1% (PC1: 58.4%, PC2: 19.7%) of the interaction effect between the clones and test sites, respectively. In the AMMI biplot, the closer the value to the coordinate origin, the more stable the clone was, and the farther away the value from the coordinate origin, the better the ability to distinguish the environment [42]. According to the selection rate of 20%, the excellent clones G14, G10, G03, G08, and G4 performed best in tree height stability. CM and YJ demonstrated the strongest ability to distinguish the tree height. G16, G21, G08, G10, and G05 performed best with respect to the stability of the ground diameter, and CM had the strongest ability to distinguish the ground diameter. In conclusion, G08 and G10 had the best stability in terms of tree height and ground diameter, and CM had the strongest ability to distinguish tree height and ground diameter, followed by YJ.

3.6. WAASB Stability Index Analysis

Based on the WAASB stability index scores (Table 6), it is apparent that the clones G08, G15, G12, G02, and G18 exhibit exceptional stability in tree height, whereas the clones G24, G08, G19, G18, and G02 demonstrate robust ground diameter stability. Notably, the clones G08, G02, and G18 stand out as having remarkable stability in both tree height and ground diameter.

3.7. GGE Biplot Analysis

The GGE biplot of the tree height and the ground diameter showed that AXIS1 and AXIS2 explained 92.20% and 90.56% of the total variation (the sum of the clone effect and the interaction effect between the clones and test sites), indicating that the inference results were reliable (Figure 3).
In the GGE biplot of “Mean vs. Stability” (Figure 3a,b), a green line with an arrow shows the Average Environmental Coordinate (AEC), and a line perpendicular to AEC shows the Average Environmental Axis (AEA) [43]. The intersection of the AEC and AEA represents the average value of the clones, with a higher yield being associated with a longer vertical distance from AEA in the direction of an arrow and a better stability being represented with a shorter vertical distance from the AEC. The clones G11, G13, G05, G10, and G04 exhibited high yield and good stability in tree height. The clones G11, G05, G10, G21, and G04 demonstrated high yield and good stability in ground diameter. The YJ site had good stability and a high yield in both tree height and ground diameter. The ZZ site had a poor yield and poor stability, while the CM site had the highest yield, but poor stability.
The outermost clones in the GGE biplot of “Which won where” formed a polygon shown by green straight lines, with red vertical lines from the origin dividing it into sectors, and test sites into groups (Figure 3c,d). Based on the tree height, the five sites divided the test points into two groups, with CM, PX, FS, and YJ as one group and ZZ as another group. G11 and G13 produced good tree height adaptability in CM, PX, FS, and YJ, whereas G20, G10, G05, G02, G04, G06, and G08 produced good tree height adaptability in ZZ. All of the test sites were concentrated in one area based on base diameter, which shows that G11, G10, G05, G04, G21, and G03 had strong adaptability in these test sites (Figure 3d). Overall, clone G11 had the strongest adaptability of growth traits among all the sites.
In the GGE biplot of tree height and ground diameter, the vectors of all the test sites form acute angles with the AEC, indicating that all the test sites are representative and suitable for Qi-nan clones (Figure 4a,b). It can be seen that PX, FS, and YJ have smaller angles with the AEC and stronger representativeness, while CM and ZZ have larger angles with the AEC and weaker representativeness. The length of the vector indicates the distinguishing force of the biplot; the longer the vector is, the stronger the distinguishing force. It can be seen that the CM and YJ vectors are longer and have strong discriminative ability, while the PX and ZZ vectors are shorter and have weaker discriminative ability. YJ is a representative and distinguishing test site, and it is the most suitable test site for breeding clones with good growth and stable yield.
In the ranking GGE biplot (Figure 5), the center represents the ideal test site or clone, and the closer a test site or clone is to the center, the better its overall performance. The ranking of the ideal test sites and the ideal clones in the biplot based on tree height was CM > YJ > FS > PX > ZZ and G13 > G11 > G05 > G10 > G04, respectively (Figure 5a,b). The ranking of the ideal test sites and the ideal clones in the biplot based on the ground diameter was CM > YJ > FS > PX > ZZ and G11 > G10 > G05 > G04 > G21, respectively (Figure 5c,d). Overall, CM and YJ were found to be the ideal test sites, and G11, G13, G05, G10, and G04 were the ideal clones.

4. Discussion

4.1. Growth of Qi-Nan Clones

Excellent varieties should be high-yield, stable, and have certain adaptability [44]. Generally, the growth performance of different germplasms in different test sites is not consistent, which is not only related to the stability of the germplasm genotype itself, but also to its adaptability [45]. In this study, the growth performance of Qi-nan clones in multiple environments showed extremely significant differences (p < 0.01) among the clones, the test sites, and the interactions between the clones and the test sites (Table 2). Regarding the differences among the sites, the obvious coefficient of variation in the tree height and the ground diameter may be related to the environment conditions [46]. Similar phenomena were also reported for woody species, including Eucalyptus clones, pinus pinaster, and polar clones [18,26,47], in which the different results were considered related to the varying adaptability of the clones. The difference among the clones can be attributed to their substantial genetic variations [5], and different provenances [48]. The variance difference between the clones and the sites in the interaction regarding the growth performance of the Qi-nan clones indicated that environmental conditions comprised the most important growth-limiting factor [19]. Therefore, it is necessary to study the multi-site breeding of Qi-nan clones.
There was a significant positive correlation between the height and ground diameter of the Qi-nan clones and the mean annual temperature and mean annual sunshine duration. These correlations may be related to the traits of A. sinensis, as the original species demonstrated good growth in tropical and subtropical zones [49]. This evidence reveals why the growth performance of the CM and YJ clones showed an outyield effect in comparison to the others. Further research showed significant variations in the 1-year survival and 3-year preservation rates of the Qi-nan clones across environments. Notably, the CM, YZ, and ZZ sites with abundant mean annual precipitation have significantly higher rates, suggesting that the growth of Qi-nan clones is regulated by multiple factors, with yearly rainfall being a key influence [50]. Thus, it is recommended that Qi-nan clones be planted in tropical–subtropical regions with ample light, heat, and rain for optimal growth. Additionally, multi-environment breeding should be implemented to further enhance their adaptability.

4.2. Genetic Parameters for the Growth of Qi-Nan Clones

The ranges of the PCV and GCV values of growth traits among the clones were from 16.56% to 32.09% and from 5.24% to 27.06%, respectively, with both reaching a moderate degree of variation [18,51], suggesting that there are genetically excellent clones within the experimental group [52]. Moreover, the PCVs of the growth traits of the clones in each test site were found to be greater than the GCVs, indicating that the traits of the Qi-nan clones were not only controlled by heredity, but also affected by the environment [53]. Among the growth traits, the PCV and GCV values of the ground diameter were generally higher than those of tree height, which may be due to stronger environmental influence or larger differences in genetic control and environmental stress among the clones, resulting in increased phenotypic variations in the ground diameter [54].
The repeatability values of the tree height and ground diameter of the Qi-nan clones in the CM, YJ, PX, and FS sites were all greater than 0.80, which indicates that these test sites could effectively maintain the genetic stability of the clones [55]. Meanwhile, the repeatability values of the clones at the ZZ site were only 0.50 and 0.69, respectively, which was consistent with the low GCV and high PCV values, indicating that the clones in this area were greatly affected by the environment, which limited the normal growth of some of them [56]. In addition, it was found that the repeatability of the ground diameter was generally greater than that of the tree height, and the higher the repeatability, the more credible the test results [57]. Therefore, the selection of excellent Qi-nan clones should be based on ground diameter and supplemented by tree height.

4.3. Evaluation of Qi-Nan Clones

In this study, the AMMI model and GGE biplot were used to clarify the sums of the squares of interactions, which produced more than 78% and 90% of growth traits, respectively, indicating that it was feasible and reliable to use these two methods to evaluate the clones and adaptation areas of the Qi-nan clones in this experiment [15,58].
Regarding high yield, in the GGE biplot, clones G11, G13, G05, G10, and G04 were selected in the high-yield group for tree height, while clones G11, G05, G10, G21, and G04 were selected in the high-yield group for ground diameter. This consistency with the results of multiple comparisons (Table 3) suggests that the GGE biplot can be used to effectively evaluate the high yield potential of clones [59], and this is supported by the fact that it has already been utilized in several studies [19,60].
In terms of stability, the results of the WAASB stability index differ significantly from the other two models, possibly due to the large differences in three-year preservation rates among different clones, while the WAASB stability index is based on absolute bias, resulting in insensitivity to outliers [61]. The AMMI model and the GGE biplot were used to jointly screen the clones with good tree height stability, including G03 and G04, and the clones with good ground diameter stability, including G16, G21, and G10. However, there were differences in the identification of certain clones, such as the fact that the tree height stability of clones G01, G13, G18, and G23 in the GGE biplot were higher than that of G04, whereas their stability in the AMMI model was inferior to that of G04. Similar differences in the evaluation between the AMMI and the GGE were also found in previous studies [62,63]. The coefficient variation of tree height in different test sites further revealed that G01 (with a coefficient of variation of 18.6%), G13 (15.7%), G18 (16.9%), and G23 (15.3%) were all higher than G04 (14.7%), which was stable under both models. Therefore, the stability of clones G01, G13, G18, and G23 is worse than G04, indicating that the stability of the AMMI model is more reliable [28,64]. Clone G10, with its stable ground diameter in the two models, also had similar results when compared to G01, G11, and G18. The reason for this was because the GGE biplot evaluated stability while taking into account yield [63], while the AMMI evaluated stability more holistically because stability is not related to yield [28].
In terms of adaptability, clones G11 and G13 in the GGE biplot had better tree height adaptability in the CM, PX, FS, and YJ sites. Clones G20, G10, G05, G02, G04, G06, and G08 had better tree height adaptability in the ZZ site. Clones G11, G10, G05, G04, G21, and G03 had good ground diameter adaptability at the five test sites, indicating that these clones had specific adaptability in the corresponding test sites [65]. However, 64% and 72% of the clones in the GGE biplot of the tree height and the ground diameter, respectively, were not in the corresponding sector area of the test site, and they were considered unsuitable for the test environment [59]. Therefore, more test sites with different environmental types need to be carried out in future breeding tests.
In summary, clone G04 produced the best comprehensive performance in terms of yield, stability, and adaptability. Clones G05 and G10 also showed good growth, albeit with moderate tree height stability in certain areas; notably, despite having average stability, clone G11 excelled in both yield and adaptability, making it an excellent candidate for breeding. Meanwhile, clone G13 was grew well in all of the sites except the ZZ site. This is consistent with the ideal clone results (G11, G13, G05, G10, and G04) in the GGE biplot of “Ranking Genotypes”. Moreover, we noticed that G11 showed moderate stability while having a high ranking in the GGE biplot. We believe that the excellent yield and adaptability of G11 led to this contrast, indicating the comprehensiveness of the GGE biplot in clone evaluation [66].

4.4. Evaluations of Site Adaptability for Qi-Nan Clones

Multi-environment trials were carried out to scientifically evaluate the site adaptability of the germplasm, which was not only conducive to providing a basis for the scientific promotion of germplasm [16], effectively controlling the application risk of germplasm, but also conducive to determining the adaptation range of the germplasm [19]. In our study, the AMMI model and the GGE biplot demonstrated that the Qi-nan clones at the YJ site exhibited the best comprehensive growth performance, displaying good discrimination and representativeness. The growth performance at the CM site was outstanding, featuring good discrimination, but lacking representativeness, which could potentially be attributed to the slightly lower mean annual precipitation recorded when compared to the YJ site and poor stability. The comprehensive growth performance at the FS site was general. The stability at the PX site was remarkable, but growth was hindered by limited rainfall, resulting in poor yield. Due to its high latitude and altitude, the ZZ site had the smallest frost-free period, resulting in poor stability and a poor yield, but it also had a certain degree of differentiation, which might aid in identifying cold-tolerant clone species.
Upon comprehensive evaluation, the ranking of ideal test sites in the GGE biplot was CM > YJ > FS > PX > ZZ, which is consistent with the gradient order of the mean annual temperature and mean annual sunshine duration (Table 5). This indicates that within a certain range, the higher the temperature and the more sufficient the light, the more conducive it is to tree growth [67]. Therefore, the YJ site is the most suitable experimental site for screening high-yield and stable clones. The CM site stands out for its high-yield potential but lacks stability; however, with timely irrigation during later maintenance, it could become the best-performing experimental site. The FS site is a suitable location for the growth of Qinan clones, though their performance is generally average. The PX site requires long-term watering for good growth. The temperature of the ZZ site is relatively low, and it is necessary to breed cold-resistant clones to promote growth.

5. Conclusions

There exist significant differences in the clones and sites, and the interactions between the clones and test sites. Genetic breeding is an efficient pathway for the potential selection of Qi-nan germplasms in southern China. The phenotypic variation coefficient (PCV), the genetic variation coefficient (GCV), and the repeatability of the ground diameter are generally greater than the tree height, which may indicate that the selection of excellent Qi-nan clones should focus on ground diameter. In the comprehensive evaluation of clones and site adaptability for Qi-nan clones, the AMMI model was better than the GGE biplot in screening the stability of clones, but the GGE biplot was more comprehensive in evaluating varieties. The southern region of China is suitable for the growth of Qi-nan agarwood, and different varieties have different growth performances in different regions. According to climate conditions. Targeted management measures should also be implemented in different growing regions according to the climate conditions. This study focused on the early growth performance of Qi-nan clones; thus, long-term observations are needed in the future. Additionally, a comprehensive evaluation of growth and agarwood traits should be carried out to verify the accuracy of clone selection. Furthermore, the evaluation of different environments is needed to improve the accuracy of genetic improvements.

Author Contributions

Conceptualization, H.H. and Z.H.; methodology, X.L. (Xiaojin Li) and D.X.; software, H.H. and Z.H.; validation, X.L. (Xiaofei Li), J.H. and Z.H.; formal analysis, X.L. (Xiaofei Li) and X.F.; investigation, H.H., X.L. (Xiaofei Li), X.F. and J.H.; resources, Z.H.; data curation, J.H. and Z.C.; writing—original draft preparation, H.H.; writing—review and editing, Z.H. and D.X.; visualization, Z.H.; supervision, D.X.; project administration, X.L. (Xiaojin Liu); funding acquisition, Y.S. and D.X. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Science Innovation Projects of Guangdong Province (2023KJCX009), the Key Research and Development Project of National Forestry and Grassland Administration (GLM2021-21), the Fundamental Research Funds of Chinese Academy of Forestry (CAFYBB2021QD003-03, CAFYBB2023MB008), Guangxi Forestry Science and Technology Project (2022-02), and Guangzhou Collaborative Innovation Center on Science-Tech of Ecology and Landscape (202206010058).

Data Availability Statement

The data presented in this study are available on request to the corresponding authors, and the dataset was jointly completed by the team, so the data are not publicly available.

Acknowledgments

We thank the Tropical Forestry Experimental Center of the Chinese Academy of Forestry, Baishengyuan Agarwood Planting Base in Yangjiang City, Gaoming Forestry Institute in Foshan City, Meilangwan Agarwood Planting Base in Chengmai City, and Tianma Forest Farm in Zhangzhou City for their support and assistance.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

References

  1. Li, W.; Chen, H.; Wang, H.; Mei, W.; Dai, H. Natural Products in Agarwood and Aquilaria Plants: Chemistry, Biological Activities and Biosynthesis. Nat. Prod. Rep. 2021, 38, 528–565. [Google Scholar] [CrossRef]
  2. Gao, M.; Han, X.; Sun, Y.; Chen, H.; Yang, Y.; Liu, Y.; Meng, H.; Gao, Z.; Xu, Y.; Zhang, Z.; et al. Overview of Sesquiterpenes and Chromones of Agarwood Originating from Four Main Species of the Genus Aquilaria. RSC Adv. 2019, 9, 4113–4130. [Google Scholar] [CrossRef] [PubMed]
  3. Xu, Y.; Zhang, Z.; Wang, M.; Wei, J.; Chen, H.; Gao, Z.; Sui, C.; Luo, H.; Zhang, X.; Yang, Y.; et al. Identification of Genes Related to Agarwood Formation: Transcriptome Analysis of Healthy and Wounded Tissues of Aquilaria sinensis. BMC Genom. 2013, 14, 227. [Google Scholar] [CrossRef] [PubMed]
  4. Kang, Y.; Liu, P.; Lv, F.; Zhang, Y.; Yang, Y.; Wei, J. Genetic Relationship and Source Species Identification of 58 Qi-Nan Germplasms of Aquilaria Species in China That Easily Form Agarwood. PLoS ONE 2022, 17, e0270167. [Google Scholar] [CrossRef] [PubMed]
  5. Hou, W.; Liu, P.; Liu, Y.; Kang, Y.; Yang, Y.; Zhang, Y.; Gao, Z.; Yu, M.; Feng, J.; Lv, F.; et al. Chi-Nan Agarwood Germplasms Constitute a New Chemotype of Aquilaria sinensis (Lour.) Spreng. Ind. Crops Prod. 2022, 187, 115494. [Google Scholar] [CrossRef]
  6. Yu, M.; Liu, Y.; Feng, J.; Chen, D.; Yang, Y.; Liu, P.; Yu, Z.; Wei, J. Remarkable Phytochemical Characteristics of Chi-Nan Agarwood Induced from New-Found Chi-Nan Germplasm of Aquilaria Sinensis Compared with Ordinary Agarwood. Int. J. Anal. Chem. 2021, 2021, 5593730. [Google Scholar] [CrossRef] [PubMed]
  7. Li, X.; Cui, Z.; Liu, X.; Hong, Z.; Zhang, P.; Xu, D. Comparative Morphological, Anatomical and Physiological Analyses Explain the Difference of Wounding-Induced Agarwood Formation between Ordinary Agarwood Nongrafted Plants and Five Grafted Qi-Nan Clones (Aquilaria sinensis). Forests 2022, 13, 1618. [Google Scholar] [CrossRef]
  8. Li, X.; Fang, X.; Cui, Z.; Hong, Z.; Liu, X.; Li, G.; Hu, H.; Xu, D. Anatomical, Chemical and Endophytic Fungal Diversity of a Qi-Nan Clone of Aquilaria sinensis (Lour.) Spreng with Different Induction Times. Front. Plant Sci. 2024, 15, 1320226. [Google Scholar] [CrossRef]
  9. Zhang, P.; Li, X.; Cui, Z.; Xu, D. Morphological, Physiological, Biochemical and Molecular Analyses Reveal Wounding-Induced Agarwood Formation Mechanism in Two Types of Aquilaria sinensis (Lour.) Spreng. Ind. Crops Prod. 2022, 178, 114603. [Google Scholar] [CrossRef]
  10. Fang, X.; Li, X.; Zhang, Q.; Hu, H.; Hong, Z.; Liu, X.; Cui, Z.; Xu, D. Physiological and Endophytic Fungi Changes in Grafting Seedlings of Qi-Nan Clones (Aquilaria sinensis). Forests 2024, 15, 106. [Google Scholar] [CrossRef]
  11. Cañas-Gutiérrez, G.P.; Sepulveda-Ortega, S.; López-Hernández, F.; Navas-Arboleda, A.A.; Cortés, A.J. Inheritance of Yield Components and Morphological Traits in Avocado Cv. Hass from “Criollo” “Elite Trees” via Half-Sib Seedling Rootstocks. Front. Plant Sci. 2022, 13, 843099. [Google Scholar] [CrossRef] [PubMed]
  12. Studnicki, M.; Kang, M.S.; Iwańska, M.; Oleksiak, T.; Wójcik-Gront, E.; Mądry, W. Consistency of Yield Ranking and Adaptability Patterns of Winter Wheat Cultivars between Multi-Environmental Trials and Farmer Surveys. Agronomy 2019, 9, 245. [Google Scholar] [CrossRef]
  13. Greveniotis, V.; Bouloumpasi, E.; Zotis, S.; Korkovelos, A.; Kantas, D.; Ipsilandis, C.G. Genotype-by-Environment Interaction Analysis for Quantity and Quality Traits in Faba Beans Using AMMI, GGE Models, and Stability Indices. Plants 2023, 12, 3769. [Google Scholar] [CrossRef] [PubMed]
  14. Li, Y.; Suontama, M.; Burdon, R.D.; Dungey, H.S. Genotype by Environment Interactions in Forest Tree Breeding: Review of Methodology and Perspectives on Research and Application. Tree Genet. Genomes 2017, 13, 60. [Google Scholar] [CrossRef]
  15. Ma, C.; Liu, C.; Ye, Z. Influence of Genotype × Environment Interaction on Yield Stability of Maize Hybrids with AMMI Model and GGE Biplot. Agronomy 2024, 14, 1000. [Google Scholar] [CrossRef]
  16. Zhang, H.; Zhou, X.; Gu, W.; Wang, L.; Li, W.; Gao, Y.; Wu, L.; Guo, X.; Tigabu, M.; Xia, D.; et al. Genetic Stability of Larix Olgensis Provenances Planted in Different Sites in Northeast China. For. Ecol. Manag. 2021, 485, 118988. [Google Scholar] [CrossRef]
  17. Zhang, Q.; Pei, X.; Lu, X.; Zhao, C.; Dong, G.; Shi, W.; Wang, L.; Li, Y.; Zhao, X.; Tigabu, M. Variations in Growth Traits and Wood Physicochemical Properties among Pinus Koraiensis Families in Northeast China. J. For. Res. 2022, 33, 1637–1648. [Google Scholar] [CrossRef]
  18. Jiang, L.; Pei, X.; Hu, Y.; Chiang, V.L.; Zhao, X. Effects of Environment and Genotype on Growth Traits in Poplar Clones in Northeast China. Euphytica 2021, 217, 169. [Google Scholar] [CrossRef]
  19. Liu, X.; Zhao, Q.; Yin, P.; Li, H.; Li, X.; Wu, L.; Li, Y.; Hu, Y.; Zhao, X. Variation and Stability Analysis of Growth Traits of Poplar Clones in the Seedling Stage in Northeast China. J. For. Res. 2023, 34, 1107–1116. [Google Scholar] [CrossRef]
  20. Pei, X.; Jiang, L.; Ahmed, A.K.M.; Yu, H.; Chong, R.; You, X.; Zhao, X. Growth Variations and Stability Analyses of Seven Poplar Clones at Three Sites in Northeast China. J. For. Res. 2021, 32, 1673–1680. [Google Scholar] [CrossRef]
  21. Liepe, K.J.; Van Der Maaten, E.; Van Der Maaten-Theunissen, M.; Kormann, J.M.; Wolf, H.; Liesebach, M. Ecotypic Variation in Multiple Traits of European Beech: Selection of Suitable Provenances Based on Performance and Stability. Eur. J. For. Res. 2024, 143, 831–845. [Google Scholar] [CrossRef]
  22. Dias, A.; Gaspar, M.J.; Carvalho, A.; Pires, J.; Lima-Brito, J.; Silva, M.E.; Louzada, J.L. Within- and between-Tree Variation of Wood Density Components in Pinus nigra at Six Sites in Portugal. Ann. For. Sci. 2018, 75, 58. [Google Scholar] [CrossRef]
  23. De Araujo, M.J.; De Paula, R.C.; Campoe, O.C.; Carneiro, R.L. Adaptability and Stability of Eucalypt Clones at Different Ages across Environmental Gradients in Brazil. For. Ecol. Manage. 2019, 454, 117631. [Google Scholar] [CrossRef]
  24. Lebot, V.; Ranaivoson, L. Eucalyptus Genetic Improvement in Madagascar. For. Ecol. Manage. 1994, 63, 135–152. [Google Scholar] [CrossRef]
  25. Akbarpour, O.; Dehghani, H.; Sorkhi, B.; Gauch, J. Evaluation of Genotype×environment Interaction in Barley (Hordeum vulgare L.) Based on AMMI Model Using Developed SAS Program. J. Agric. Sci. Technol. 2014, 16, 909–920. Available online: https://jast.modares.ac.ir/article-23-1523-en.html (accessed on 1 July 2024).
  26. Correia, I.; Alía, R.; Yan, W.; David, T.; Aguiar, A.; Almeida, M.H. Genotype × Environment Interactions in Pinus Pinaster at Age 10 in a Multienvironment Trial in Portugal: A Maximum Likelihood Approach. Ann. For. Sci. 2010, 67, 612. [Google Scholar] [CrossRef]
  27. Mafouasson, H.N.A.; Gracen, V.; Yeboah, M.A.; Ntsomboh-Ntsefong, G.; Tandzi, L.N.; Mutengwa, C.S. Genotype-by-Environment Interaction and Yield Stability of Maize Single Cross Hybrids Developed from Tropical Inbred Lines. Agronomy 2018, 8, 62. [Google Scholar] [CrossRef]
  28. Pupin, S.; Silva, P.H.M.; Piotto, F.A.; Miranda, A.C.; Zaruma, D.U.G.; Sebbenn, A.M.; Moraes, M.L.T. Genotype x Environment Interaction, Stability, and Adaptability in Progenies of Eucalyptus urophylla S.T. BLAKE Using the AMMI Model. Silvae Genet. 2018, 67, 51–56. [Google Scholar] [CrossRef]
  29. Crossa, J. Statistical Analyses of Multilocation Trials. In Advances in Agronomy; Brady, N.C., Ed.; Academic Press: Cambridge, MA, USA, 1990; Volume 44, pp. 55–85. [Google Scholar] [CrossRef]
  30. Rodrigues, P.C.; Malosetti, M.; Gauch, H.G., Jr.; van Eeuwijk, F.A. A Weighted AMMI Algorithm to Study Genotype-by-Environment Interaction and QTL-by-Environment Interaction. Crop Sci. 2014, 54, 1555–1570. [Google Scholar] [CrossRef]
  31. de Oliveira, R.L.; Von Pinho, R.G.; Balestre, M.; Ferreira, D.V. Evaluation of Maize Hybrids and Environmental Stratification by the Methods AMMI and GGE Biplot. Crop Breed. Appl. Biotechnol. 2010, 10, 247–253. [Google Scholar] [CrossRef]
  32. Sousa, L.B.; Hamawaki, O.T.; Nogueira, A.P.O.; Batista, R.O.; Oliveira, V.M.; Hamawaki, R.L. Evaluation of Soybean Lines and Environmental Stratification Using the AMMI, GGE Biplot, and Factor Analysis Methods. Gen. Mole Resea. 2015, 14, 12660–12674. [Google Scholar] [CrossRef] [PubMed]
  33. Alwala, S.; Kwolek, T.; McPherson, M.; Pellow, J.; Meyer, D. A Comprehensive Comparison between Eberhart and Russell Joint Regression and GGE Biplot Analyses to Identify Stable and High Yielding Maize Hybrids. Field Crops Res. 2010, 119, 225–230. [Google Scholar] [CrossRef]
  34. Olivoto, T.; Lúcio, A.D.C.; da Silva, J.A.G.; Sari, B.G.; Diel, M.I. Mean Performance and Stability in Multi-Environment Trials II: Selection Based on Multiple Traits. Agron. J. 2019, 111, 2961–2969. [Google Scholar] [CrossRef]
  35. Mohammadi, R.; Jafarzadeh, J.; Armion, M.; Hatamzadeh, H.; Roohi, E. Clustering Stability Methods towards Selecting Best Performing and Stable Durum Wheat Genotypes. Euphytica. 2023, 219, 109. [Google Scholar] [CrossRef]
  36. Nataraj, V.; Bhartiya, A.; Singh, C.P.; Devi, H.N.; Deshmukh, M.P.; Verghese, P.; Singh, K.; Mehtre, S.P.; Kumari, V.; Maranna, S.; et al. WAASB-based Stability Analysis and Simultaneous Selection for Grain Yield and Early Maturity in Soybean. Agron. J. 2021, 113, 3089–3099. [Google Scholar] [CrossRef]
  37. Liang, D.; Zhang, X.; Wang, C.; Wang, X.; Li, K.; Liu, G.; Zhao, X.; Qu, G.-Z. Evaluation of Betula Platyphylla Families Based on Growth and Wood Property Traits. For. Sci. 2018, 64, 663–670. [Google Scholar] [CrossRef]
  38. Ryan, M.G.; Stape, J.L.; Binkley, D.; Alvares, C.A. Cross-Site Patterns in the Response of Eucalyptus Plantations to Irrigation, Climate and Intra-Annual Weather Variation. For. Ecol. Manage. 2020, 475, 118444. [Google Scholar] [CrossRef]
  39. Baye, T. Genotypic and Phenotypic Variability in Vernonia Galamensis Germplasm Collected from Eastern Ethiopia. J. Agric. Sci. 2002, 139, 161–168. [Google Scholar] [CrossRef]
  40. Yue, H.; Li, H.; Xu, L.; Bu, J.; Wei, J.; Chen, S.; Peng, H.; Xie, J.; Shang, S.; Jiang, X. Analysis of Genotype-Environment Interactions of Silage Maize Cultivars under Environmental Trials. Bangladesh J. Bot. 2020, 49, 55–63. [Google Scholar] [CrossRef]
  41. Yan, W. GGEbiplot—A Windows Application for Graphical Analysis of Multienvironment Trial Data and Other Types of Two-Way Data. Agron. J. 2001, 93, 1111–1118. [Google Scholar] [CrossRef]
  42. Mehareb, E.M.; Osman, M.A.M.; Attia, A.E.; Bekheet, M.A.; Abo Elenen, F.F.M. Stability Assessment for Selection of Elite Sugarcane Clones across Multi-Environment Based on AMMI and GGE-Biplot Models. Euphytica 2022, 218, 95. [Google Scholar] [CrossRef]
  43. Zhang, W.; Hu, J.; Yang, Y.; Lin, Y. One Compound Approach Combining Factor-Analytic Model with AMMI and GGE Biplot to Improve Multi-Environment Trials Analysis. J. For. Res. 2020, 31, 123–130. [Google Scholar] [CrossRef]
  44. Gill, H.S.; Halder, J.; Zhang, J.; Brar, N.K.; Rai, T.S.; Hall, C.; Bernardo, A.; Amand, P.S.; Bai, G.; Olson, E.; et al. Multi-Trait Multi-Environment Genomic Prediction of Agronomic Traits in Advanced Breeding Lines of Winter Wheat. Front. Plant Sci. 2021, 12, 709545. [Google Scholar] [CrossRef] [PubMed]
  45. da Silva, P.H.M.; Araujo, M.J.; Lee, D.J.; Bush, D.; Baroni, G.R.; de Paula, R.C. Adaptability and Stability of Novel Eucalypt Species and Provenances across Environments in Brazil at Two Assessment. New For. 2022, 53, 779–796. [Google Scholar] [CrossRef]
  46. Nabeshima, E.; Kubo, T.; Hiura, T. Variation in Tree Diameter Growth in Response to the Weather Conditions and Tree Size in Deciduous Broad-Leaved Trees. For. Ecol. Manage. 2010, 259, 1055–1066. [Google Scholar] [CrossRef]
  47. Souza, T.D.S.; Ramalho, M.A.P.; Lima, B.M.D.; Rezende, G.D.S.P. Performance of Eucalyptus Clones According to Environmental Conditions. Sci. For. 2017, 45, 601–610. [Google Scholar] [CrossRef]
  48. Zhao, X.; Bian, X.; Li, Z.; Wang, X.; Yang, C.; Liu, G.; Jiang, J.; Kentbayev, Y.; Kentbayeva, B.; Yang, C. Genetic Stability Analysis of Introduced Betula pendula, Betula kirghisorum, and Betula pubescens Families in Saline-Alkali Soil of Northeastern China. Scand. J. For. Res. 2014, 29, 639–649. [Google Scholar] [CrossRef]
  49. Liu, C.; Zhou, G.; Liu, J. Isolation and Screening of Fungi for Enhanced Agarwood Formation in Aquilaria sinensis Trees. PLoS ONE 2024, 19, e0304946. [Google Scholar] [CrossRef]
  50. Jia, B.; Sun, H.; Shugart, H.H.; Xu, Z.; Zhang, P.; Zhou, G. Growth Variations of Dahurian Larch Plantations across Northeast China: Understanding the Effects of Temperature and Precipitation. J. Environ. Manage. 2021, 292, 112739. [Google Scholar] [CrossRef]
  51. Singh, S.; Prakash, A.; Chakraborty, N.R.; Wheeler, C.; Agarwal, P.K.; Ghosh, A. Genetic Variability, Character Association and Divergence Studies in Jatropha curcas for Improvement in Oil Yield. Trees 2016, 30, 1163–1180. [Google Scholar] [CrossRef]
  52. Stovall, J.P.; Carlson, C.A.; Seiler, J.R.; Fox, T.R.; Yanez, M.A. Growth and Stem Quality Responses to Fertilizer Application by 21 Loblolly Pine Clones in the Virginia Piedmont. For. Ecol. Manag. 2011, 261, 362–372. [Google Scholar] [CrossRef]
  53. Xiao, Y.; Ma, W.; Lu, N.; Wang, Z.; Wang, N.; Zhai, W.; Kong, L.; Qu, G.; Wang, Q.; Wang, J. Genetic Variation of Growth Traits and Genotype-by-Environment Interactions in Clones of Catalpa Bungei and Catalpa fargesii f. Duclouxii. Forests 2019, 10, 57. [Google Scholar] [CrossRef]
  54. Scheiner, S.M.; Goodnight, C.J. The Comparison of Phenotypic Plasticity and Genetic Variation in Populations of the Grass Danthonia Spicata. Evolution 1984, 38, 845–855. [Google Scholar] [CrossRef] [PubMed]
  55. Yin, S.; Xiao, Z.; Zhao, G.; Zhao, X.; Sun, X.; Zhang, Y.; Wang, F.; Li, S.; Zhao, X.; Qu, G.-Z. Variation Analyses of Growth and Wood Properties of Larix olgensis Clones in China. J. For. Res. 2017, 28, 687–697. [Google Scholar] [CrossRef]
  56. Huang, R.; Wang, R.; Wei, R.; Yan, S.; Wu, G.; Zheng, H. Selection for Both Growth and Wood Properties in Chinese Fir Breeding Parents Based on a 6-Year Grafted Clone Test. Forests 2023, 14, 1774. [Google Scholar] [CrossRef]
  57. Sotelo Montes, C.; Hernández, R.E.; Beaulieu, J.; Weber, J.C. Genetic Variation in Wood Color and Its Correlations with Tree Growth and Wood Density of Calycophyllum Spruceanum at an Early Age in the Peruvian Amazon. New For. 2008, 35, 57–73. [Google Scholar] [CrossRef]
  58. Omrani, A.; Omrani, S.; Khodarahmi, M.; Shojaei, S.H.; Illés, Á.; Bojtor, C.; Mousavi, S.M.N.; Nagy, J. Evaluation of Grain Yield Stability in Some Selected Wheat Genotypes Using AMMI and GGE Biplot Methods. Agronomy 2022, 12, 1130. [Google Scholar] [CrossRef]
  59. Braga, R.C.; Paludeto, J.G.Z.; Souza, B.M.; Aguiar, A.V.; Pollnow, M.F.M.; Carvalho, A.G.M.; Tambarussi, E.V. Genetic Parameters and Genotype × Environment Interaction in Pinus taeda Clonal Tests. For. Ecol. Manag. 2020, 474, 118342. [Google Scholar] [CrossRef]
  60. Yuan, C.; Zhang, Z.; Jin, G.; Zheng, Y.; Zhou, Z.; Sun, L.; Tong, H. Genetic Parameters and Genotype by Environment Interactions Influencing Growth and Productivity in Masson Pine in East and Central China. For. Ecol. Manag. 2021, 487, 118991. [Google Scholar] [CrossRef]
  61. Olivoto, T.; Lúcio, A.D.C.; da Silva, J.A.G.; Marchioro, V.S.; de Souza, V.Q.; Jost, E. Mean Performance and Stability in Multi-Environment Trials I: Combining Features of AMMI and BLUP Techniques. Agron. J. 2019, 111, 2949–2960. [Google Scholar] [CrossRef]
  62. Yan, W.; Kang, M.S.; Ma, B.; Woods, S.; Cornelius, P.L. GGE Biplot vs. AMMI Analysis of Genotype-by-Environment Data. Crop Sci. 2007, 47, 643–653. [Google Scholar] [CrossRef]
  63. Zhang, H.; Feng, Z.; Wang, J.; Yun, X.; Qu, F.; Sun, C.; Wang, Q. Genotype by Environment Interaction for Grain Yield in Foxtail Millet (Setarai italica) Using AMMI Model and GGE Biplot. Plant Growth Regul. 2023, 99, 101–112. [Google Scholar] [CrossRef]
  64. Saeidnia, F.; Majidi, M.M.; Dehghani, M.R.; Saeidi, G.; Mirlohi, A. Drought Tolerance and Stability of Native Iranian and Foreign Tall Fescue Genotypes: Comparison of AMMI and GGE Biplot Analyses. Agron. J. 2022, 114, 2180–2185. [Google Scholar] [CrossRef]
  65. Bertoldo, J.G.; Coimbra, J.L.M.; Nodari, R.O.; Guidolin, A.F.; Hemp, S.; Barili, L.D.; Vale, N.M.; Rozzeto, D.S. Stratification of the State of Santa Catarina in Macro-Environments for Bean Cultivation. Crop. Breed. Appl. Biotechnol. 2009, 9, 335–343. [Google Scholar] [CrossRef]
  66. Taheripourfard, Z.; Izadi-Darbandi, A.; Ghazvini, H.; Ebrahimi, M.; Mortazavian, S.M.M.; Abdipour, M. Identifying Superior Barley (Hordeum vulgare L.) Genotypes Using GGE-Biplot across Warm and Moderate Environments under Irrigated Conditions in Iran. Crop Breed. J. 2017, 7, 23–35. [Google Scholar] [CrossRef]
  67. Midmore, D.J.; Cartwright, P.M.; Fischer, R.A. Wheat in Tropical Environments. II. Crop Growth and Grain Yield. Field Crops Res. 1984, 8, 207–227. [Google Scholar] [CrossRef]
Figure 1. Map showing positions at five different sites in South China. The red triangle mark represents the test site. CM, FS, PX, YJ and ZZ denote the Chengmai site, Foshan site, Pingxiang site, Yangjiang site, and Zhangzhou site, respectively.
Figure 1. Map showing positions at five different sites in South China. The red triangle mark represents the test site. CM, FS, PX, YJ and ZZ denote the Chengmai site, Foshan site, Pingxiang site, Yangjiang site, and Zhangzhou site, respectively.
Forests 15 01753 g001
Figure 2. Tree height analysis based on AMMI biplot (a); ground diameter analysis based on AMMI biplot (b). The green numbers stand for clones and the blue characters stand for test sites.
Figure 2. Tree height analysis based on AMMI biplot (a); ground diameter analysis based on AMMI biplot (b). The green numbers stand for clones and the blue characters stand for test sites.
Forests 15 01753 g002
Figure 3. GGE biplot of “Mean vs. Stability” analysis based on tree height (a) and ground diameter (b); GGE biplot of “Which won where” analysis based on tree height (c); GGE biplot of “Which won where” analysis based on ground diameter (d).
Figure 3. GGE biplot of “Mean vs. Stability” analysis based on tree height (a) and ground diameter (b); GGE biplot of “Which won where” analysis based on tree height (c); GGE biplot of “Which won where” analysis based on ground diameter (d).
Forests 15 01753 g003
Figure 4. GGE biplot of “Discriminativeness vs. Representativeness” analysis based on tree height (a) and ground diameter (b).
Figure 4. GGE biplot of “Discriminativeness vs. Representativeness” analysis based on tree height (a) and ground diameter (b).
Forests 15 01753 g004
Figure 5. GGE biplot of “Ranking Environments” analysis and GGE biplot of “Ranking Genotype” analysis based on tree height (a,b); GGE biplot of “Ranking Environments” analysis and GGE biplot of “Ranking Genotypes” analysis based on ground diameter (c,d).
Figure 5. GGE biplot of “Ranking Environments” analysis and GGE biplot of “Ranking Genotype” analysis based on tree height (a,b); GGE biplot of “Ranking Environments” analysis and GGE biplot of “Ranking Genotypes” analysis based on ground diameter (c,d).
Forests 15 01753 g005aForests 15 01753 g005b
Table 1. Main geographical and environmental characteristics of test sites.
Table 1. Main geographical and environmental characteristics of test sites.
Test SiteLongitudeLatitudeAltitude (m)Frost-Free Season (days)Mean Annual Precipitation (mm)Mean Annual Temperature (°C)Mean Annual Sunshine Duration (h)
Yangjiang (YJ)112°09′24″ E21°59′23″ N16358188622.32012
Foshan (FS)112°44′58″ E22°43′24″ N85356157221.61629
Chenmai (CM)110°01′58″ E19°52′13″ N38.66365178623.82059
Pingxiang (PX)106°55′12″ E22°03′20″ N251344150021.51614
Zhangzhou (ZZ)117°21′03″ E24°21′37″ N446319173621.31764
Notes: Climate information comes from China Meteorological Data Network (https://data.cma.cn/, accessed on 6 April 2024).
Table 2. Results of ANOVA analysis of tree height and ground diameter within each site.
Table 2. Results of ANOVA analysis of tree height and ground diameter within each site.
SitesSource of VariationsTree HeightGround Diameter
DFSSMSFSSMSF
MEClone24194.008.0865.743 **1617.0267.4092.497 **
Site4622.90155.731266.871 **2203.53550.90756.299 **
Clone: Site96167.701.7514.208 **786.508.2011.248 **
Residuals4329518.310.12 2998.190.70
CMClone24205.058.5473.566 **1010.7242.1154.700 **
Residuals983111.030.11 746.910.77
FSClone2435.241.4722.432 **355.3314.8129.414 **
Residuals66543.530.07 334.740.50
PXClone2415.790.664.819 **155.206.476.702 **
Residuals57177.940.13 551.110.96
YJClone24111.964.6640.446 **859.9035.8359.674 **
Residuals1174135.410.12 704.930.60
ZZClone2415.860.664.133 **125.715.245.969 **
Residuals56289.870.15 493.200.87
Notes: ME denotes the joint multi-environments. DF, SS, MS, and F denote the source of the variations, sum of squares, mean square, and F-test, respectively. **: significant difference at the 0.01 level.
Table 3. Tree height (m), ground diameter (cm), 1-year survival rate (%), and 3-year preservation rate (%) of 25 clones and at five different test sites.
Table 3. Tree height (m), ground diameter (cm), 1-year survival rate (%), and 3-year preservation rate (%) of 25 clones and at five different test sites.
Clones/SitesTree Height (m)Ground Diameter (cm)1-Year Survival Rate (%)3-Year Preservation Rate (%)
G012.11 ± 0.42 k4.06 ± 0.83 l87.11 ± 17.77 abc68.01 ± 25.51 cde
G022.60 ± 0.56 cdef5.29 ± 1.21 efg91.57 ± 12.04 abc82.02 ± 17.55 abc
G032.47 ± 0.55 fghi5.63 ± 1.19 d90.57 ± 10.81 abc78.59 ± 21.51 abcd
G042.64 ± 0.58 bcde5.74 ± 1.36 cd89.14 ± 14.21 abc81.43 ± 15.54 abc
G052.72 ± 0.58 abc6.06 ± 1.31 b90.35 ± 12.25 abc83.13 ± 14.67 abc
G062.47 ± 0.38 fghi4.98 ± 0.93 hi85.22 ± 15.23 abc79.78 ± 20.19 abcd
G072.53 ± 0.63 defgh5.18 ± 1.04 fgh88.67 ± 16.02 abc75.67 ± 26.95 abcde
G082.55 ± 0.56 defg4.85 ± 1.12 ij94.67 ± 9.07 ab87.56 ± 13.75 ab
G092.53 ± 0.61 defgh5.5 ± 1.52 def89.22 ± 13.77 abc77.33 ± 23.42 abcd
G102.68 ± 0.56 abcd5.94 ± 1.32 bc95.44 ± 6.19 a88.56 ± 15.96 a
G112.80 ± 0.78 a6.7 ± 1.69 a88.56 ± 18.6 abc82.67 ± 21.68 abc
G122.39 ± 0.54 hij4.87 ± 1.11 ij83.89 ± 19.87 abc71.44 ± 32.51 abcde
G132.76 ± 0.80 ab5.58 ± 1.46 de87.67 ± 12.15 abc76.22 ± 15.25 abcd
G142.41 ± 0.54 ghij5.31 ± 1.22 efg83.89 ± 16.7 abc74.33 ± 21.50 abcde
G152.59 ± 0.66 cdef5.31 ± 1.56 efg81.88 ± 19.62 abcd62.03 ± 21.24 defg
G162.45 ± 0.52 fghij4.68 ± 0.98 j85.48 ± 16.14 abc70.81 ± 21.56 abcde
G172.31 ± 0.41 j4.61 ± 0.99 jk82.56 ± 23.51 abc69.44 ± 27.06 bcde
G181.88 ± 0.33 l3.86 ± 0.88 l78.31 ± 23.48 cd49.94 ± 18.55 fg
G191.93 ± 0.40 l3.95 ± 0.73 l84.39 ± 16.9 abc57.90 ± 23.7 efg
G202.50 ± 0.51 efghi5.23 ± 1.06 fgh87.89 ± 11.05 abc76.31 ± 13.66 abcd
G212.65 ± 0.76 bcde5.74 ± 1.37 cd84.67 ± 14.72 abc75.44 ± 16.57 abcde
G222.31 ± 0.52 j5.25 ± 1.09 fgh85.07 ± 22.24 abc70.38 ± 23.93 abcde
G232.06 ± 0.34 k4.36 ± 0.86 k69.75 ± 23.45 d48.22 ± 19.95 g
G242.36 ± 0.70 ij5.30 ± 1.38 efg80.86 ± 19.80 bcd62.57 ± 25.02 defg
G252.38 ± 0.49 ij4.99 ± 1.26 ghi83.81 ± 18.05 abc65.78 ± 23.86 cdef
CM3.23 ± 0.68 a6.57 ± 1.54 a94.42 ± 5.13 a90.82 ± 8.08 a
FS2.14 ± 0.34 d4.66 ± 0.99 c74.49 ± 24.20 c58.67 ± 26.96 d
PX2.10 ± 0.39 d4.6 ± 1.07 cd77.82 ± 16.41 c51.36 ± 33.43 d
YJ2.61 ± 0.46 b5.57 ± 1.14 b90.40 ± 13.25 b81.93 ± 15.27 b
ZZ2.20 ± 0.44 c4.50 ± 1.05 d90.00 ± 14.82 b71.33 ± 24.01 c
Notes: The data in the table are means ± standard deviation. Different lowercase letters indicate significant differences at the 0.05 level.
Table 4. Results of the phenotypic variation coefficient (PCV), genetic variation coefficient (GCV), and repeatability (R) of tree height and ground diameter at five different test sites.
Table 4. Results of the phenotypic variation coefficient (PCV), genetic variation coefficient (GCV), and repeatability (R) of tree height and ground diameter at five different test sites.
SitesTraitsPCV (%)GCV (%)R
CMTree height21.9418.710.96
Ground diameter24.5820.660.98
FSTree height16.5610.760.90
Ground diameter21.3914.670.85
PXTree height32.0927.060.80
Ground diameter23.2910.330.85
YJTree height18.0311.740.93
Ground diameter20.4814.990.98
ZZTree height19.975.240.50
Ground diameter23.408.510.69
Notes: PCV–phenotypic coefficient of variation; GCV–genotypic coefficient of variation. R: the repeatability of the clone.
Table 5. Results of the correlation analysis between the tree height, the ground diameter, the 1-year survival rate, and the 3-year preservation rate, and geographical and environmental factors.
Table 5. Results of the correlation analysis between the tree height, the ground diameter, the 1-year survival rate, and the 3-year preservation rate, and geographical and environmental factors.
TraitsLongitudeLatitudeAltitudeFrost-Free SeasonMean Annual PrecipitationMean Annual TemperatureMean Annual Sunshine Duration
Tree height−0.18−0.81−0.570.600.650.98 **0.9 *
Ground diameter−0.28−0.87−0.680.710.620.99 **0.89 *
1-year survival rate0.28−0.3−0.060.020.88 *0.640.89 *
3-year preservation rate0.22−0.49−0.40.360.91 *0.810.97 **
Notes: ** and * indicate the highly significant and significant difference at 0.01 or 0.05 levels, respectively.
Table 6. The WAASB stability index score and ranking of tree height and ground diameter.
Table 6. The WAASB stability index score and ranking of tree height and ground diameter.
ClonesTree HeightGround Diameter
WAASB ValueRankingWAASB ValueRanking
G010.40240.3722
G020.0840.165
G030.22150.3116
G040.1190.3115
G050.25180.3420
G060.1070.188
G070.30210.3217
G080.0610.112
G090.45250.6125
G100.23170.3218
G110.26190.2110
G120.0830.2411
G130.33220.3823
G140.21140.187
G150.0820.209
G160.30200.3219
G170.34230.3521
G180.0950.154
G190.11100.153
G200.0960.2914
G210.14110.166
G220.19130.2712
G230.16120.2713
G240.1180.091
G250.22160.5324
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hu, H.; Xu, D.; Li, X.; Fang, X.; Cui, Z.; Liu, X.; Hao, J.; Su, Y.; Hong, Z. Screening and Site Adaptability Evaluation of Qi-Nan Clones (Aquilaria sinensis) in Southern China. Forests 2024, 15, 1753. https://doi.org/10.3390/f15101753

AMA Style

Hu H, Xu D, Li X, Fang X, Cui Z, Liu X, Hao J, Su Y, Hong Z. Screening and Site Adaptability Evaluation of Qi-Nan Clones (Aquilaria sinensis) in Southern China. Forests. 2024; 15(10):1753. https://doi.org/10.3390/f15101753

Chicago/Turabian Style

Hu, Houzhen, Daping Xu, Xiaofei Li, Xiaoying Fang, Zhiyi Cui, Xiaojin Liu, Jian Hao, Yu Su, and Zhou Hong. 2024. "Screening and Site Adaptability Evaluation of Qi-Nan Clones (Aquilaria sinensis) in Southern China" Forests 15, no. 10: 1753. https://doi.org/10.3390/f15101753

APA Style

Hu, H., Xu, D., Li, X., Fang, X., Cui, Z., Liu, X., Hao, J., Su, Y., & Hong, Z. (2024). Screening and Site Adaptability Evaluation of Qi-Nan Clones (Aquilaria sinensis) in Southern China. Forests, 15(10), 1753. https://doi.org/10.3390/f15101753

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop