Plasticity of primary and secondary growth dynamics in Eucalyptushybrids: a quantitative genetics and QTL mapping perspective
© Bartholomé et al.; licensee BioMed Central Ltd. 2013
Received: 29 April 2013
Accepted: 14 August 2013
Published: 26 August 2013
The genetic basis of growth traits has been widely studied in forest trees. Quantitative trait locus (QTL) studies have highlighted the presence of both stable and unstable genomic regions accounting for biomass production with respect to tree age and genetic background, but results remain scarce regarding the interplay between QTLs and the environment. In this study, our main objective was to dissect the genetic architecture of the growth trajectory with emphasis on genotype x environment interaction by measuring primary and secondary growth covering intervals connected with environmental variations.
Three different trials with the same family of Eucalyptus urophylla x E. grandis hybrids (with different genotypes) were planted in the Republic of Congo, corresponding to two QTL mapping experiments and one clonal test. Height and radial growths were monitored at regular intervals from the seedling stage to five years old. The correlation between growth increments and an aridity index revealed that growth before two years old (r = 0.5; 0.69) was more responsive to changes in water availability than late growth (r = 0.39; 0.42) for both height and circumference. We found a regular increase in heritability with time for cumulative growth for both height [0.06 - 0.33] and circumference [0.06 - 0.38]. Heritabilities for incremental growth were more heterogeneous over time even if ranges of variation were similar (height [0-0.31]; circumference [0.19 to 0.48]). Within the trials, QTL analysis revealed collocations between primary and secondary growth QTLs as well as between early growth increments and final growth QTLs. Between trials, few common QTLs were detected highlighting a strong environmental effect on the genetic architecture of growth, validated by significant QTL x E interactions.
These results suggest that early growth responses to water availability determine the genetic architecture of total growth at the mature stage and highlight the importance of considering growth as a composite trait (such as yields for annual plants) for a better understanding of its genetic bases.
KeywordsGrowth Growth modelling Plasticity QTL x E interaction Heritability Eucalyptus
Production of wood biomass is a mandatory trait in any forest tree breeding program regardless of the final use be it pulp and paper, energy, construction or engineered wood products. Understanding the contribution of genetic and environmental factors as well as their interaction in tree growth and adaptation is a prerequisite for accelerating tree domestication to meet the increasing demand for wood . In this context, fast growing trees such as Eucalyptus will play a major role, not only as wood supply, but also as a model system to decipher the determinism of growth. Indeed, eucalypts are cultivated worldwide on more than twenty million hectares and are the most planted hardwoods in the world . This genus comprises about 700 species  distributed naturally over a wide range of pedoclimatic conditions. Within that diversity, a few species or inter-specific hybrids combining the adaptive capacities and fast growth rates of the parental species (e.g. E. urophylla x E. grandis[4, 5]), are widely used in planted forest.
Wood biomass formation is a dynamic process resulting from a wide array of physiological processes  some of which (e.g. photosynthesis) are dependent on environmental variations. It is also a product of height and circumference growth, involving the activity of the shoot apical meristem  and the vascular cambium . The dynamics of both primary and secondary growth have been extensively described and modelled up to the rotation age in many forest tree species [9–12] with the main objectives of predicting biomass production and optimizing sylvicultural practices. The first models corresponded to yield models and with the improved understanding of plant/tree physiological functioning and its interplay with climatic factors, process or mechanistic models were subsequently developed taking into account carbon allocation, nutrient and water availability, and climate effects [13, 14]. These models were often developed for a particular species [15, 16] without explicitly taking into account the genetic diversity within species, while a genetic effect of the response to environmental variation has been described in Eucalyptus on different scales of growth analysis (day, month, year) [17, 18].
The genetic determinism of tree growth and its trend over-time have been fairly well studied in most forest tree species of commercial interest [19–24]. As regards time trends, the Franklin model , suggesting three phases in variance components for growth traits (juvenile, mature and adult phases), was tested in different species [26–28]. In short rotation species with rapid initial growth such as in Eucalyptus, a more stable trend was usually observed up to the rotation age . Similarly, the ratio between dominance and additive variances has been analysed over time, with a major share of additive variance in Pinus taeda and more balanced in Eucalyptus. Age-age genetic correlations were found to be generally high between successive measurements and tended to decrease when the time lag increased i.e. between juvenile and mature stages [27, 30, 31]. Overall, these studies reveal low to moderate heritability for growth traits, indicating a strong environmental effect on the determinism of growth. More recent studies have tried to link genetic of growth response to environmental variation [23, 32, 33]. In Eucalyptus globulus, Costa e Silva et al.  reported a genetic effect on radial growth in response to temperature, rainfall and solar radiation. A similar result was presented by a more recent study showing a genetic effect for suceptibility to drought damage in E. globulus. Nevertheless, there is still a need to improve our understanding of genetic and environmental determinants of growth, through measurements covering intervals better connected with environmental variations.
An analysis of the genetic architecture of growth traits (i.e. the number, map location and effect of quantitative trait loci, QTL) can also provide a better understanding of the genetic basis of growth. Thus, over the past two decades, many studies have been undertaken to detect QTLs controlling part of the phenotypic variation of wood biomass in forest trees such as pines , poplar  and eucalyptus . In Eucalyptus, several genomic regions have been detected harbouring QTLs with low to moderate effects on the phenotypic variation of growth traits [39–44]. These studies, although carried out with a relatively small number of growth measurement time points revealed QTL instability over time [40, 41]. The functional mapping approach taking into account the growth trajectory in QTL analysis  showed similar results with early and late QTLs in poplar . However, rather few studies connecting the genetic architecture of growth traits and different environmental conditions have been carried out [47–50]. In Eucalyptus Freeman et al. [51, 52] as well as Teixeira et al. , highlighted major QTL x Environment interactions. All these studies conducted on an inter-annual scale are consistent with the idea that growth is an integrative trait involving a large number of ecophysiological  and molecular  processes on which environmental conditions can act to determine the final genetic architecture of growth. Therefore, understanding the interplay between early growth and abiotic factors (e.g. alternation between dry and wet seasons in the tropics) and how environmental variations impact the genetics of growth is a key issue, especially as part of Eucalyptus breeding programmes.
In this context, the objective of this study was to dissect the genetic basis of the growth trajectory and its interplay with seasonal variation in water availability. To that end, one interspecific cross of Eucalyptus urophylla x E. grandis was planted in three different field trials in the Republic of Congo where contrasting environmental conditions in term of water availability (rainy vs. dry seasons), were suitable for observing the effect of drought on growth. We used two complementary types of growth trait: integrative traits (cumulative growth, growth curve parameters) and responsive traits (sub-annual growth increments) in combination with climatic data to characterize growth response to changes in water availability and reveal its genetic component. This comprehensive dissection of growth enabled us to characterize the genetic determinism of primary (height) and secondary (circumference) growth using conventional quantitative genetics and QTL detection approaches.
Plant material and field experiments
A full-sib family generated from a single controlled cross between Eucalyptus urophylla (accession 14.144) x E. grandis (accession 9.21) was used in this study in three environmental settings corresponding to three field experiments. The first two trials were planted using a single tree plot design with seedlings: in April 1993 with 201 genotypes (referred to as P93) and already used in previous QTL studies [40, 55, 56] and in April 1997 (P97) with 190 genotypes (different from P93). The third trial was planted in November 1998 (P98) and corresponded to a clonal test with 60 genotypes in common with P93 replicated by rooted-cuttings with 10 copies per genotype.
The three trials were established in the Republic of Congo, in the Pointe-Noire region (4.8° S, 12° E). P93 and P97 were planted at a density of 667 trees/ha, whereas P98 was planted at 800 trees/ha. To reduce border effects a triple border row was planted around all the plantations.
DNA extractions for the full-sib progenies were carried out using dry leaves according to Doyle and Doyle .
where p is the total monthly rainfall and t the monthly mean temperature. According to De Martonne , a threshold of IDM =15 was chosen as an upper limit to define the dry season. Climatic data (temperature and rainfall) were recorded at Pointe-Noire about 25 km from the field trials.
Where Trait n is the cumulative height or circumference at “n” months old.
Where l(t) is the cumulative growth at time t; Asym is the maximum growth (the asymptote); lrc is the logarithm of the growth constant and c0 is the theoretical age for which l(t) is zero. Growth modelling was performed with R software (Version 2.15.2; http://cran.r-project.org/) using the R package: nlme[71, 72]. Trees with a high proportion of missing data were removed, based on the duration and the phase (exponential growth phase or asymptotic phase) at which missing data occurred. Selected trees presented on average 0.3, 2.2 and 0.5% of growth curve that is missing for P93, P97 and P98 respectively. For each selected genotype, a nonlinear model (monomolecular) was fitted to obtain the above-mentioned growth curve parameters (GCP) for both height and circumference. The model fitted the data well as illustrated for circumference in Additional file 1. Moreover the pseudo R-squared calculated for all adjusted models were higher than 0.9. Predicted values were also calculated for P97 for the same time points as in P93 (14, 26, 39, 51 and 59 months old).
Broad-sense heritability of growth
Where is the genetic variance and is the phenotypic variance .
Molecular markers and linkage mapping
For each parental genotype (E. urophylla and E. grandis) and each trial (P93 and P97), a genetic map was constructed according to a two-way pseudo-test cross mapping strategy . For P93, the two genetic linkage maps were constructed using RAPD SSR, EST, and STS markers genotyped on 201 individuals . For P97, only RAPD markers selected on the basis of their map position in the P93 linkage maps were used for genetic mapping with 190 individuals. The set of common polymorphic markers genotyped in both P93 and P97, enabled the unambiguous identification of homologous linkage groups (LGs) between the two parental maps of the same species. In addition, the codominant markers used in P93, enabled the identification of orthologous regions between E. urophylla and E. grandis.
JoinMap® version 4.1 [76, 77] was used to perform the linkage analysis. Chi-squared tests were performed to test whether markers followed the expected Mendelian segregation ratios. Distorted markers (p <0.01) were excluded from linkage analysis. Markers with more than 40% of missing data were also excluded. For the four maps, markers were grouped into LGs at the independence logarithm-of-the-odds (LOD) score ≥ 4.0. Then, markers within linkage groups were ordered using the Monte Carlo maximum likelihood mapping algorithm [78, 79]. Finally, mapping quality tools were used to classify markers leading to a poorer fit as accessory markers, located near the closest framework markers. Accessory markers with significant deviation from expected Mendelian segregation (p < 0.01) were indicated by * on the genetic maps. Genetic distances in centiMorgan (cM) were calculated with the Kosambi mapping function . The LG nomenclature corresponds to that defined by Brondani et al. .
Quantitative trait locus (QTL) analyses were performed with R software using the qtl package . We carried out simple and composite interval mapping using the multiple imputation method of Sen and Churchill  with a step interval of 1 cM, a genotyping error rate of 0.001 and 512 imputations per genotypes. It has been shown that this imputation method performs better with missing genotypic data . The whole-genome significance thresholds for all traits were determined with 10,000 permutations . QTLs that met or exceeded 95th and 90th percentiles were assessed to be significant or suggestive, respectively. Confidence intervals for each QTL position were obtained using Bayes credible intervals (BCI). Genetic maps and QTL positions were drawn using MapChart .
Detection of QTL x environment interaction
where y is the individual tree measurement, μ is the general mean, Trail is the fixed effect of the trial (P93 or P97), Marker is the fixed effect of the closest marker associated with a QTL detected by CIM, Trail × Marker is the interaction between the two fixed effects and ϵ is the residual random effect .
Trend in phenotypic variation and correlation
The pattern of phenotypic variations for cumulative growth (Ht, Cir) was very similar between trials. Indeed the variance increased gradually over time and the coefficient of phenotypic variation (CVp) for Ht remained rather stable ([0.07-0.09], [0.13-0.19] and [0.12-0.24] for P93, P97 and P98, respectively, Additional file 3). The same tendency was observed for Cir. The CVp of growth increments for P97 and P98 were more variable compared to cumulative growth, as exemplified by CVp values for Ht_Inc ranging from 0.18 to 0.58 for P97 and from 0.12 to 0.82 for P98. A similar range of values was observed for Cir_Inc. Given the large number of measurements in P97, we were able to characterize the evolution of CVp over time for Ht_Inc and Cir_Inc. Two growth rate peaks CVp were observed during the second and the third dry periods (Additional file 4), corresponding to a lower growth increment on average. These results suggest that growth rate differences between genotypes were enhanced during these two dry periods.
Pearson’s coefficients of correlation (r) between growth increments and I DM
Pearson’s r (p-value)
Mean of Pearson’s r [range]
% of significant correlation
0.58 [-0.26, 0.96]
0.44 [-0.36, 0.9]
0.79 [0.02, 0.98]
0.49 [-0.08, 0.92]
Trend in the genetic control of growth traits
Genetic map statistics for E. urophylla and E. grandis in both the P93 and P97 trials
Number of LGs
Total map length (cM)
Average LG length (cM)
Number of framework markers
Mean distance between markers
Number of common framework markers between P93 and P97
Missing data (%)
Among the framework markers, 72 markers for the female map and 57 markers for the male map were common between P93 and P97 maps. These “bridge markers” enabled the identification 13 homologous regions for the E. urophylla parent covering 73% (P93) and 90% (P97) of total map length. For E. grandis, 11 homologous regions were identified, covering of 79% (P93) and 89% (P97) of the total map length. As expected, the order between markers was well conserved between P93 and P97: 89% and 93% of the common markers presented the same order for the E. urophylla and the E. grandis maps, respectively. The four genetic maps with links between bridge markers revealing homologous and orthologous regions are represented in Additional file 8.
Genetic architecture of growth traits
QTL detection in P93
A total of 19 QTLs were detected for the female parent (E. urophylla) and 9 for the male parent (E. grandis) at the genome-wide type I error rate of 10% (Additional file 9). 71% of the QTLs remained significant at the 5% genome-wide error rate. These QTLs were associated with 15 and 8 different traits for E. urophylla and E. grandis respectively; representing 8 and 6 different regions distributed over 5 different linkage groups for each parental map. For the E. urophylla parent, the variance explained by all the detected QTLs associated with one trait ranged from 5.5% (Ht39) to 21.2% (Asym_h) with an average of 11.5% (± 3.4). For the male parent, this proportion of variance ranged from 4.2% (Ht39) to 13.2% (Cir39) with an average of 7.6% (± 3.1). QTLs for E. urophylla were located on LG1, LG3, LG4, LG6 and LG7 with two hotspots (LG3 at 11-12 cM and LG6 at 22-27 cM) gathering 16% and 47% of the total number of QTLs respectively. QTLs for E. grandis were located on LG1, LG2, LG3, LG5 and LG10 with two hotspots (LG5 at 113-116 cM and LG10 at 27 cM) gathering 22% and 33% of the total number of QTLs, respectively. The distribution of QTLs between the two parents was different, of the seven LGs with QTLs only two (LG1 and LG3) were common to both parents. Two specific QTLs were found for the growth curve parameter on the E. urophylla map on LG3 (for Asym_h and lrc_c). No QTL for GCP were found on the E. grandis map.
QTL detection in P97
Unlike P93, fewer QTLs were detected for the female parent (33) than for the male parent (72) at the genome-wide type I error rate of 10%, while 74% of the QTLs remained significant at 5%. These QTLs associated with 29 and 62 different traits, represented 8 and 9 different regions distributed on 6 and 7 different LGs for E. urophylla and E. grandis, respectively (Additional file 10). For E. urophylla, the variance explained by all the detected QTLs associated with one trait ranged from 5.1% (Cir11) to 15.3% (Cir52_62) with an average of 7.6% (± 2.5). This proportion of variance ranged from 5.1% (Cir8_9) to 17% (Ht52) with an average of 10.5% (± 2.9) for the male parent (E. grandis). QTLs for E. urophylla were located on LG2, LG3, LG5, LG6, LG8 and LG10 with three hotspots (LG5 at 47 cM, LG8 at 41-44 cM and LG8 at 66-68 cM ) gathering 12%, 34% and 16% of the total number of QTLs, respectively. QTLs for E. grandis were located on LG2, LG3, LG4, LG5, LG6, LG8 and LG10 with three hotspots (LG2 at 0-10 cM, LG4 at 73 cM and LG8 at 22-42 cM) gathering 6%, 6% and 60% of the total number of QTLs, respectively. All LGs with QTLs were common to both parents except for LG4, for which no QTL was detected in E. urophylla. One specific QTL region was found for the GCP located on the E. grandis map on LG10 (Asym_c, lrc_c and c0_c). The other GCP QTLs co-located with QTLs for different traits.
QTL comparison across trials
Co-locating QTLs between P93 and P97 were found on LG2 and 5 for the E. grandis parent only, involving strong overlapping of confidence intervals. These two co-localizations involved hotspots of one of two trials but represented few of all the detected QTLs. Indeed a Jaccard Index (IJ, ) was calculated between the two trials (P93 and P97) for E. grandis (IJ = 12.5%). This difference in the genetic architecture of growth traits between P93 and P97 suggested a strong QTL x environment interaction (QTL x E) that was tested using the predicted values obtained for P97 based on growth curve modelling. Thus, using closest markers already associated with QTLs, a two way analysis of variance was conducted to test QTL x E interaction for growth traits at ages 14, 26, 39, 51 and 59 months for P93 and predicted values at the same ages for P97. GCP were also used for both trials (Additional file 11). We found that most of tested markers revealed a significant interaction with the environment (trial) for both parents, E. urophylla 60% and 80%, E. grandis 40% and 90% for P93 and P97 respectively. For both trials, QTL x E was more significant for E. grandis than for E. urophylla.
Incremental growth at the early stage is related to water availability and determines biomass production at rotation age
We combined the measurement of growth on a sub-annual scale with a monthly characterisation of climate variation using an aridity index (IDM), and first showed that seasonal changes in water availability affected secondary and primary growth at the early stage of tree development. Indeed, IDM was found to be strongly correlated with radial growth increment before two years old, while it was less related to growth after two years. Besides, the correlation at early stage was more significant for radial growth (r = 0.69, p < 0.01) than for vertical growth (r = 0.5, p < 0.01), indicating that processes underlying secondary growth are more prone to environmental variation than those involved in primary growth. Other studies have already reported a relationship between apical or radial incremental growth and rainfall on different time scales (day, month and season) in different species [87–89]. In this study, we showed that the first two years corresponded to the highest biomass productivity period, during which 60% of the final height and 65% of final circumference was reached, representing a biomass productivity of 20-27 m3.ha-1.year-1. Therefore, the early growing period constitutes a key period with a strong impact on final growth, as illustrated by the high correlations values between each measurement point and the final growth at 60 months old. In tropical eucalypt plantations, it was found that the maximum growth rate occurred around two years, followed by a decline in production after canopy closure [90–92]. Therefore, our results suggest that climatic conditions underlying the aridity index (mainly rainfall) have a strong impact on early growth response and consequently determine biomass production at rotation age.
Early growth response to water availability is genetically variable
The small sample size (60 genotypes) used to calculate H2 limited the accuracy of the estimation. However, the comparison of time trends in H2 between cumulative and incremental growths highlighted contrasting patterns over time which was significant only for Ht. Indeed, cumulative growth displayed an asymptotic increase up to rotation age, contrary to incremental growth which showed high temporal heterogeneity. For Eucalyptus grown in Congo, H2 for cumulative growth calculated using a wider genetic background were found to range from 0.10 to 0.5  and to stabilize after two or three years. The same trend in H2 for cumulative growth was found in Eucalyptus in other environments [93–95]. Interestingly, this stabilization of cumulative growth heritabilities occurred at the same time as the decline in wood biomass productivity , with a decrease in or stabilization of the G x E interaction [95–97]. It could thus be suggested that the integrative nature of cumulative growth at a mature stage aggregates past genetic responses to environmental variations.
The heterogeneity of H2 for incremental growth observed for both height and radial growth suggests a variation in the impact of environmental conditions on growth response. The slowest values of incremental growth heritability were observed during the third dry period. In Pinus Pinaster, such heterogeneous trends in H2 have been shown for yearly incremental growths  and the variability described was significantly linked to environmental effects (year effects). On a lower scale, Parisi  also reported a heterogeneous pattern for narrow and broad-sense heritabilities for incremental height growth in Pinus taeda. In addition to the contrast between cumulative and incremental growth, H2 patterns for growth increment revealed opposite patterns during the 25th and 35th months between primary and secondary growth responses. Costa and Durel  also reported a clear difference between primary and secondary growth increments and suggested more stable genetic control for secondary growth response.
Between-tree variation in this early growth response, depicted by the distribution of phenotypic correlations between growth increments and IDM during early growth, suggested a genetic effect as a major component in this response. Other studies have already suggested a genetic effect in response to water availability. Drew et al.  observed such contrasting behaviour in response to climatic variations (drought) between different Eucalyptus hybrids i.e. an E. grandis x E. urophylla susceptible to drought with intermittent growth, and a less susceptible E. grandis x E. camaldulensis maintaining a continuous growth. Similar results were found between E. alba and E. urophylla x E. grandis in Congo , associated with a significant genotype x environment interaction for traits related to height growth and carbon isotope discrimination (δ13C). These genetic differences in environmental susceptibility were also suggested by Bouvet et al.  comparing different eucalyptus clones. Interestingly, in the present study, we found a QTL for the correlation between IDM and radial growth increment, suggesting that early growth response to water availability is indeed genetically controlled.
The genetic architecture of final growth is driven by the early growth response to the environment
Given the high macro-synteny and colinearity between eucalyptus genomes, we were able to compare the position of the mapped QTLs to that detected in other studies in different genotypes of the same species or different species. This comparison enabled us to identify two LGs where growth QTLs were frequently detected [39, 41–43, 52]. As already mentioned by Thumma et al. , and supported by our result, LG5 was found to contain QTLs for height and radial growth in different Eucalyptus species [39, 41–43, 52] as well as QTLs for foliar chemicals . To a lesser extent, LG8, which was associated with both incremental and cumulative growth in our study, was also associated with growth-QTLs in three other studies [41, 43, 52], as well as QTLs for rooting ability  and QTLs for resistance to biotic stresses . In addition to these two main LGs, and in agreement with the oligogenic model, i.e. few genes with large effects detected as QTLs and many genes with small effects controlled quantitative traits , QTLs were found to be distributed across the other nine chromosomes [39, 41–44, 52].
In this study, a similarity index revealed that most of the growth QTLs were not shared between two experimental trials showing different environmental trajectories (see Figure 1). The statistical power of QTL detection and of estimated effects is directly related to sample size, so called Beavis effect . Therefore QTLs with intermediate or small effects are likely to be not detected with a small sample size (e.g. 100 genotypes). Thus one cannot rule out the fact that QTL instability between P93 and P97 is the result of a rather low sample size: 200 and 190 genotypes for P93 and P97, respectively. This lack of QTL stability could also reveal true biological effect, i.e. a genotypic sensitivity to environmental conditions, as supported by the significant QTL x E interaction found in both parents, and more particularly in E. grandis. Interestingly, the two studied species display contrasting adaptability to Congolese conditions with E. urophylla better adapted than E. grandis. The instability of growth-QTLs between trials also contrasted with the QTL stability found for other traits measured on the same offspring (wood properties and δ13C, data not shown). This difference in sensitivity to environmental variations is generally used to explain the greater PEV and larger number of wood property QTLs in comparison to growth related traits [39, 42, 106]. It has been found that growth-QTLs often presented a higher level of interaction with the environment [47, 48, 107] than wood property-QTLs [41, 52].
The temporal analysis of QTL detection showed the presence of early stage, intermediate and late QTLs. Among these three patterns of QTL expression, a large proportion of QTLs associated with cumulative growth exhibited early or intermediate patterns: 50% (P93) and 75% (P97) for E. urophylla and 75% (P93) and 33% (P97) for E. grandis. These three QTL patterns were also reported by Emebiri et al.  on three-year-old pines. This instability over time is frequently reported in forest tree literature in relation to low heritabilities associated with growth. In E. globulus, a study also reported a lack of QTL stability over time for diameter at breast height between two and six years . In the same species, Bundock et al.  found partial QTL stability for stem diameter with the same location of LOD peaks between the diameter at two and six years old. Along the same lines, Thumma et al.  also reported relatively low stability for growth QTLs (height and diameter) over time (at two different ages) in one family. In our study, partial stability was found on LG2 (E. grandis) and LG5 (E. urophylla). Moreover, stable QTLs after two or three years old were found on LG6 (E. urophylla), LG8 (E. grandis) and LG10 (E. grandis). QTL stability over the years has also been reported in gymnosperms [50, 109] but their measurements were limited to the young stage (two- to ten-year-old trees) and did not therefore reflect the genetic architecture of the whole growth trajectory.
To our knowledge, our study provides the most comprehensive view of the global instability of growth genetic architecture in a forest tree species thanks to a detailed analysis of the temporal and environmental dynamics of growth over the whole rotation period. Indeed, the fine scale characterization of both growth (multiple measurements) and the environment (aridity index) enabled a precise QTL dissection from early to late growth. In the literature, growth increments are mostly calculated on a yearly basis. Consequently, differences in genetic determinism between cumulative and incremental are less noticeable [109, 110]. In our study, the wider variety of growth responses captured by short incremental data was instrumental to characterizing the genetic architecture of growth and depicting its temporal plasticity. More remarkably, thanks to the fine characterization of growth, we found that all QTLs detected at the mature stage for cumulative growth, co-localized with QTLs for early growth increment on both parents (LG5 for E. urophylla and LG2, LG8 for E. grandis). This result, which goes far beyond a simple age-age correlation analysis, suggests a partly common genetic determinism between early growth response to environmental variation and final growth. The QTL for the correlation with IDM detected on LG2 in E. grandis, as well as the pattern of LOD scores for QTLs detected for mature growth on P93 (Additional file 12), corroborate this assumption. To sum up, our study shows that the variability of growth plasticity observed at the young stage partly determined the pattern of cumulative growth genetic architecture at rotation age.
This study showed differences at phenotypic and genetic levels (phenotypic correlations, broad sense heritability, QTLs) between cumulative and incremental growth, for both primary and secondary growth in a fast growing tree species. Growth increments (especially radial growth) were found to be highly responsive to climatic variations (seasonal changes in water availability) over the two first years, contrary to cumulative growth. However, QTL analysis showed a common genetic basis between early growth increments and final growth. Moreover, strong instability was found in growth-QTLs between the two trials, which was also supported by the analysis of QTL x E. These results suggest that early growth responses (apical and radial growth) to soil water availability partially shape the genetic determinism of late growth and thus growth trajectories. Specific genomic regions were also found for growth increments and growth curve parameters, showing the merits of a combined approach with sequential (increment) and global (curve parameters) growth traits to dissect its genetic determinism. Our results highlight the importance of considering growth as a compound trait for a better understanding of its genetic bases, and suggest to combine selection on early and cumulative growth to improve the ability of trees to produce wood biomass in an optimal way.
The authors would like to thank the technicians at CRDPI (Centre de Recherches sur la Durablilité des Plantations Industrielles located in the Republic of Congo) for trial establishment and follow-up. The authors acknowledge two anonymous reviewers for their constructive comments and suggestions on the manuscript. This study was supported by grants from FEDER (Fonds Européen de Développement Régional, ABIOGEN project, n° Presage 32973). JB received a PhD fellowship from the BIOS department of CIRAD (Centre de coopération Internationale en Recherche Agronomique pour le Développement).
- FAO: State of the world’s forests, 2009. Rome: Food and Agriculture Organization of the United Nations: 2009.Google Scholar
- FAO: State of the world’s forests: 2007. Rome: Food and Agriculture Organization of the United Nations: 2007.Google Scholar
- Coppen JJW: Eucalyptus: The Genus Eucalyptus. London: Taylor & Francis: 2002.Google Scholar
- Vigneron P: Création et amélioration de variétés hybrides d’Eucalyptus au Congo. Bois et forêts des tropiques. 1992, 234: 29-42.Google Scholar
- Denison NP, Kietzka JE: The use and importance of hybrid intensive forestry in South Africa. South African Forestry Journal. 1993, 165: 55-60. 10.1080/00382167.1993.9629390.Google Scholar
- Kozowski TT, Pallardy SG: Growth control in woody plants. New York: Academic Press: 1997.Google Scholar
- Barton MK: Twenty years on: the inner workings of the shoot apical meristem, a developmental dynamo. Dev Biol. 2010, 341 (1): 95-113. 10.1016/j.ydbio.2009.11.029.PubMedGoogle Scholar
- Lachaud S, Catesson A-M, Bonnemain J-L: Structure and functions of the vascular cambium. Comptes Rendus de l’Académie des Sciences - Series III - Sciences de la Vie. 1999, 322 (8): 633-650. 10.1016/S0764-4469(99)80103-6.Google Scholar
- Vanclay JK: Growth models for tropical forests: a synthesis of models and methods. Forest Sci. 1995, 41 (1): 7-42.Google Scholar
- Vanclay JK, Skovsgaard JP: Evaluating forest growth models. Ecol Model. 1997, 98 (1): 1-12.Google Scholar
- Zeide B: The science of forestry. J Sustain For. 2008, 27 (4): 345-473. 10.1080/10549810802339225.Google Scholar
- Landsberg JJ, Sands P: Physiological ecology of forest production: principles, processes and models. London: Elsevier Science: 2010.Google Scholar
- Mäkelä A, Landsberg J, Ek AR, Burk TE, Ter-Mikaelian M, Ågren GI, Oliver CD, Puttonen P: Process-based models for forest ecosystem management: current state of the art and challenges for practical implementation. Tree Physiol. 2000, 20 (5–6): 289-298.PubMedGoogle Scholar
- Landsberg J: Physiology in forest models: history and the future. FBMIS. 2003, 1: 49-63.Google Scholar
- Almeida AC, Landsberg JJ, Sands PJ, Ambrogi MS, Fonseca S, Barddal SM, Bertolucci FL: Needs and opportunities for using a process-based productivity model as a practical tool in Eucalyptus plantations. For Ecol Manage. 2004, 193 (1–2): 167-177.Google Scholar
- Fontes L, Landsberg J, Tome J, Tome M, Pacheco CA, Soares P, Araujo C: Calibration and testing of a generalized process-based model for use in Portuguese eucalyptus plantations. Can J For Res-Rev Can Rech For. 2006, 36 (12): 3209-3221. 10.1139/x06-186.Google Scholar
- Drew D, Downes G, Grzeskowiak V, Naidoo T: Differences in daily stem size variation and growth in two hybrid eucalypt clones. Trees - Structure and Function. 2009, 23 (3): 585-595. 10.1007/s00468-008-0303-y.Google Scholar
- Villar E: Ph.D. thesis. Plasticité phénotypique et moléculaire de deux clones d’eucalyptus sous contrainte hydrique au champ. Bordeaux: Université de Bordeaux I: 2011.Google Scholar
- Namkoong G, Kang H: Quantitative Genetics of Forest Trees. Plant Breeding Reviews. New Jersey: John Wiley & Sons, Inc:1990, 139-188.Google Scholar
- Cornelius J: Heritabilities and additive genetic coefficients of variation in forest trees. Can J Forest Res. 1994, 24 (2): 372-379. 10.1139/x94-050.Google Scholar
- Hamilton M, Potts B: Review of Eucalyptus nitens genetic parameters. N Z J For Sci. 2008, 38 (1): 102-119.Google Scholar
- Volker PW, Potts BM, Borralho NMG: Genetic parameters of intra- and inter-specific hybrids of Eucalyptus globulus and E-nitens. Tree Genetics & Genomes. 2008, 4 (3): 445-460. 10.1007/s11295-007-0122-0.Google Scholar
- Bouvet JM, Saya A, Vigneron P: Trends in additive, dominance and environmental effects with age for growth traits in Eucalyptus hybrid populations. Euphytica. 2009, 165 (1): 35-54. 10.1007/s10681-008-9746-x.Google Scholar
- Retief ECL, Stanger TK: Genetic parameters of pure and hybrid populations of Eucalyptus grandis and E-urophylla and implications for hybrid breeding strategy. South Forests. 2009, 71 (2): 133-140. 10.2989/SF.2009.71.2.8.823.Google Scholar
- Franklin EC: Model relating levels of genetic variance to stand development of 4 north-american conifers. Silvae genetica. 1979, 28 (5–6): 207-212.Google Scholar
- Cotterill PP, Dean CA: Changes in the genetic-control of growth of radiata pine to 16 years and efficiencies of early selection. Silvae genetica. 1988, 37 (3–4): 138-146.Google Scholar
- Bouvet JM, Vigneron P: Age trends in variances and heritabilities in Eucalyptus factorial mating designs. Silvae genetica. 1995, 44 (4): 206-216.Google Scholar
- Stackpole D, Vaillancourt R, De Aguigar M, Potts B: Age trends in genetic parameters for growth and wood density in Eucalyptus globulus. Tree Genetics & Genomes. 2010, 6 (2): 179-193. 10.1007/s11295-009-0239-4.Google Scholar
- Xiang B, Li B, Isik F: Time trend of genetic parameters in growth traits of Pinus taeda L. Silvae genetica. 2003, 52 (3–4): 114-121.Google Scholar
- Kremer A: Predictions of age-age correlations of total height based on serial correlations between height increments in Maritime pine (Pinus pinaster Ait.). TAG Theor Appl Genet. 1992, 85 (2): 152-158.PubMedGoogle Scholar
- Gwaze DP, Bridgwater FE, Byram TD, Woolliams JA, Williams CG: Predicting age-age genetic correlations in tree-breeding programs: a case study of Pinus taeda L. Theor Appl Genet. 2000, 100 (2): 199-206. 10.1007/s001220050027.Google Scholar
- Martinez Meier A, Sanchez G, Dalla L, Salda G, Pastorino MJM, Gautry JY, Gallo LA, Rozenberg P: Genetic control of the tree-ring response of Douglas-fir (Pseudotsuga menziesii (Mirb.) Franco) to the 2003 drought and heat-wave in France. Ann For Sci. 2008, 65 (1): 102-10.1051/forest:2007074.Google Scholar
- Marron N, Ricciotti L, Bastien C, Beritognolo I, Gaudet M, Paolucci I, Fabbrini F, Salani F, Dillen SY, Ceulemans R, et al: Plasticity of growth and biomass production of an intraspecific Populus alba family grown at three sites across Europe during three growing seasons. Can J For Res-Rev Can Rech For. 2010, 40 (10): 1887-1903. 10.1139/X10-113.Google Scholar
- Costae Silva J, Potts BM, Dutkowski GW: Genotype by environment interaction for growth of Eucalyptus globulus in Australia. Tree Genetics & Genomes. 2006, 2 (2): 61-75. 10.1007/s11295-005-0025-x.Google Scholar
- Dutkowski GW, Potts BM: Genetic variation in the susceptibility of Eucalyptus globulus to drought damage. Tree Genetics & Genomes. 2012, 8 (4): 757-773. 10.1007/s11295-011-0461-8.Google Scholar
- Ritland K, Krutovsky KV, Tsumura Y, Pelgas B, Isabel N, Bousquet J: Genetic mapping in conifers. Genetics, genomics and breeding of conifers. Edited by: Piomion C, Bousquet J, Kole C. 2011, New York: Science Publishers & CRC Press, 1121-1134.Google Scholar
- Rae AM, Street NR, Rodríguez-Acosta M: Populus trees. Forest Trees. UK: Springer: 2007, 1-28.Google Scholar
- Myburg AA, Potts BM, Marques CM, Kirst M, Gion JM, Grattapaglia D, Grima-Pettenatti J: Eucalypts. Heidelberg, Germany: Springer-Verlag GmbH: 2007.Google Scholar
- Grattapaglia D, Bertolucci F, Penchel R, Sederoff R: Genetic mapping of quantitative trait loci controlling growth and wood quality traits in Eucalyptus grandis using a maternal half-sib family and RAPD markers. Genetics. 1996, 144: 1205-1214.PubMedPubMed CentralGoogle Scholar
- Verhaegen D, Plomion C, Gion JM, Poitel M, Costa P, Kremer A: Quantitative trait dissection analysis in Eucalyptus using RAPD markers: 1. Detection of QTL in interspecific hybrid progeny, stability of QTL expression across different ages. TAG Theor Appl Genet. 1997, 95 (4): 597-608. 10.1007/s001220050601.Google Scholar
- Bundock PC, Potts BM, Vaillancourt RE: Detection and stability of quantitative trait loci (QTL) in Eucalyptus globulus. Tree Genetics & Genomes. 2008, 4 (1): 85-95.Google Scholar
- Freeman JS, Whittock SP, Potts BM, Vaillancourt RE: QTL influencing growth and wood properties in Eucalyptus globulus. Tree Genetics & Genomes. 2009, 5 (4): 713-722. 10.1007/s11295-009-0222-0.Google Scholar
- Thumma BR, Baltunis BS, Bell JC, Emebiri LC, Moran GF, Southerton SG: Quantitative trait locus (QTL) analysis of growth and vegetative propagation traits in Eucalyptus nitens full-sib families. Tree Genetics & Genomes. 2010, 6 (6): 877-889. 10.1007/s11295-010-0298-6.Google Scholar
- Kullan AR, Van Dyk M, Hefer C, Jones N, Kanzler A, Myburg A: Genetic dissection of growth, wood basic density and gene expression in interspecific backcrosses of Eucalyptus grandis and E. urophylla. BMC Genet. 2012, 13 (1): 60-PubMedPubMed CentralGoogle Scholar
- Ma C-X, Casella G, Wu R: Functional Mapping of Quantitative Trait Loci Underlying the Character Process: A Theoretical Framework. Genetics. 2002, 161 (4): 1751-1762.PubMedPubMed CentralGoogle Scholar
- Wu R, Lin M: Functional mapping: how to map and study the genetic architecture of dynamic complex traits. Nat Rev Genet. 2006, 7 (3): 229-237.PubMedGoogle Scholar
- Rae A, Pinel M, Bastien C, Sabatti M, Street N, Tucker J, Dixon C, Marron N, Dillen S, Taylor G: QTL for yield in bioenergy Populus: identifying G × E interactions from growth at three contrasting sites. Tree Genetics & Genomes. 2008, 4 (1): 97-112.Google Scholar
- Dillen SY, Storme V, Marron N, Bastien C, Neyrinck S, Steenackers M, Ceulemans R, Boerjan W: Genomic regions involved in productivity of two interspecific poplar families in Europe. 1. Stem height, circumference and volume. Tree Genetics & Genomes. 2009, 5 (1): 147-164. 10.1007/s11295-008-0175-8.Google Scholar
- Rönnberg-Wästljung AC, Glynn C, Weih M: QTL analyses of drought tolerance and growth for a Salix dasyclados × Salix viminalis hybrid in contrasting water regimes. TAG Theor Appl Genet. 2005, 110 (3): 537-549. 10.1007/s00122-004-1866-7.PubMedGoogle Scholar
- Pelgas B, Bousquet J, Meirmans P, Ritland K, Isabel N: QTL mapping in white spruce: gene maps and genomic regions underlying adaptive traits across pedigrees, years and environments. BMC Genomics. 2011, 12 (1): 145-10.1186/1471-2164-12-145.PubMedPubMed CentralGoogle Scholar
- Freeman J, Potts B, Downes G, Thavamanikumar S, Pilbeam D, Hudson C, Vaillancourt R: QTL analysis for growth and wood properties across multiple pedigrees and sites in Eucalyptus globulus. BMC Proc. 2011, 5 (Suppl 7): O8-10.1186/1753-6561-5-S7-O8.PubMed CentralGoogle Scholar
- Freeman JS, Potts BM, Downes GM, Pilbeam D, Thavamanikumar S, Vaillancourt RE: Stability of quantitative trait loci for growth and wood properties across multiple pedigrees and environments in Eucalyptus globulus. New Phytol. 2013, 198 (4): 1121-1134. 10.1111/nph.12237.PubMedGoogle Scholar
- Teixeira J, Missiaggia A, Dias D, Scarpinati E, Viana J, Paula N, Paula R, Bonine C: QTL analyses of drought tolerance in Eucalyptus under two contrasting water regimes. BMC Proc. 2011, 5 (Suppl 7): P40-10.1186/1753-6561-5-S7-P40.PubMed CentralGoogle Scholar
- Grattapaglia D, Plomion C, Kirst M, Sederoff RR: Genomics of growth traits in forest trees. Curr Opin Plant Biol. 2009, 12 (2): 148-156. 10.1016/j.pbi.2008.12.008.PubMedGoogle Scholar
- Gion J-M, Rech P, Grima-Pettenati J, Verhaegen D, Plomion C: Mapping candidate genes in Eucalyptus with emphasis on lignification genes. Mol Breed. 2000, 6 (5): 441-449. 10.1023/A:1026552515218.Google Scholar
- Gion J-M, Carouche A, Deweer S, Bedon F, Pichavant F, Charpentier J-P, Bailleres H, Rozenberg P, Carocha V, Ognouabi N, et al: Comprehensive genetic dissection of wood properties in a widely-grown tropical tree: eucalyptus. BMC Genomics. 2011, 12 (1): 301-10.1186/1471-2164-12-301.PubMedPubMed CentralGoogle Scholar
- Doyle J, Doyle J: Isolation of plant DNA from fresh tissue. Focus. 1990, 12: 13-15.Google Scholar
- Mareschal L, Nzila JDD, Turpault MP, Thongo M’Bou A, Mazoumbou JC, Bouillet JP, Ranger J, Laclau JP: Mineralogical and physico-chemical properties of Ferralic Arenosols derived from unconsolidated Plio-Pleistocenic deposits in the coastal plains of Congo. Geoderma. 2011, 162 (1–2): 159-170.Google Scholar
- Kottek M, Grieser J, Beck C, Rudolf B, Rubel F: World Map of the Koppen-Geiger climate classification updated. Meteorol Z. 2006, 15 (3): 259-263. 10.1127/0941-2948/2006/0130.Google Scholar
- De Martonne E: Une Nouvelle fonction climatologique. L’Indice d’aridité. Paris: Impr. Gauthier-Villars: 1926.Google Scholar
- Von Bertalanffy L: Quantitative Laws in Metabolism and Growth. Q Rev Biol. 1957, 32 (3): 217-231.PubMedGoogle Scholar
- Richards FJ: A Flexible Growth Function for Empirical Use. J Exp Bot. 1959, 10 (2): 290-301. 10.1093/jxb/10.2.290.Google Scholar
- Pienaar LV, Turnbull KJ: The Chapman-Richards Generalization of Von Bertalanffy’s Growth Model for Basal Area Growth and Yield in Even - Aged Stands. Forest Sci. 1973, 19 (1): 2-22.Google Scholar
- Zeide B: Analysis of Growth Equations. Forest Sci. 1993, 39 (3): 594-616.Google Scholar
- Vanclay JK: Modelling forest growth and yield: applications to mixed tropical forests. Wallingford UK: CAB International: 1994.Google Scholar
- Fekedulegn D, Mac Siurtain MP, Colbert JJ: Parameter estimation of nonlinear growth models in forestry. Silva Fennica. 1999, 33 (4): 327-336.Google Scholar
- Lei YC, Zhang SY: Features and partial derivatives of Bertalanffy-Richards growth model in forestry. Nonlinear Analysis: Modelling and Control. 2004, 9 (1): 65-73.Google Scholar
- Bouvet JM: Ph.D. Evolution de la variabilite avec l’age et correlation juvenile-adulte dans des populations d’eucalyptus. Paris: Institut National Agronomique: 1995.Google Scholar
- Amaro A, Reed D, Tome M, Themido I: Modeling dominant height growth: Eucalyptus plantations in Portugal. Forest Sci. 1998, 44 (1): 37-46.Google Scholar
- Saint-Andre L, Laclau J-P, Deleporte P, Ranger J, Gouma R, Saya A, Joffre R: A Generic Model to Describe the Dynamics of Nutrient Concentrations within Stemwood across an Age Series of a Eucalyptus Hybrid. Ann Bot. 2002, 90 (1): 65-76. 10.1093/aob/mcf146.PubMedPubMed CentralGoogle Scholar
- Pinheiro JC, Bates DM: Mixed-Effects Models in S and S-Plus. New York: Springer: 2009.Google Scholar
- Pinheiro J, Bates D, Debroy S, Sarkar D, the R Development Core Team: nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1-110 edn. 2013Google Scholar
- Hadfield JD: MCMC Methods for Multi-Response Generalized Linear Mixed Models: The MCMCglmm R Package. J Stat Softw. 2010, 33 (2): 1-22.Google Scholar
- Nyquist WE, Baker RJ: Estimation of heritability and prediction of selection response in plant populations. Crit Rev Plant Sci. 1991, 10 (3): 235-322. 10.1080/07352689109382313.Google Scholar
- Grattapaglia D, Sederoff R: Genetic Linkage Maps of Eucalyptus grandis and Eucalyptus urophylla Using a Pseudo-Testcross Mapping Strategy and RAPD Markers. Genetics. 1994, 137 (4): 1121-1137.PubMedPubMed CentralGoogle Scholar
- Stam P: Construction of integrated genetic linkage maps by means of a new computer package: join map. Plant J. 1993, 3 (5): 739-744. 10.1111/j.1365-313X.1993.00739.x.Google Scholar
- Van Ooijen JW: JoinMap 4.1, Software for the Calculation of Genetic Linkage Maps in Experimental Populations. Wageningen, Netherlands: Kyazma BV: 2011.Google Scholar
- Jansen J, De Jong AG, Van Ooijen JW: Constructing dense genetic linkage maps. TAG Theor Appl Genet. 2001, 102 (6): 1113-1122.Google Scholar
- Van Ooijen JW: Multipoint maximum likelihood mapping in a full-sib family of an outbreeding species. Genet Res. 2011, 93 (5): 343-349. 10.1017/S0016672311000279.Google Scholar
- Kosambi DD: The estimation of map distances from recombination values. Ann Hum Genet. 1943, 12 (1): 172-175. 10.1111/j.1469-1809.1943.tb02321.x.Google Scholar
- Brondani R, Williams E, Brondani C, Grattapaglia D: A microsatellite-based consensus linkage map for species of Eucalyptus and a novel set of 230 microsatellite markers for the genus. BMC Plant Biol. 2006, 6 (1): 20-10.1186/1471-2229-6-20.PubMedPubMed CentralGoogle Scholar
- Broman KW, Wu H, Sen Ś, Churchill GA: R/qtl: QTL mapping in experimental crosses. Bioinformatics. 2003, 19 (7): 889-890. 10.1093/bioinformatics/btg112.PubMedGoogle Scholar
- Sen S, Churchill GA: A Statistical Framework for Quantitative Trait Mapping. Genetics. 2001, 159 (1): 371-387.PubMedPubMed CentralGoogle Scholar
- Churchill GA, Doerge RW: Empirical Threshold Values for Quantitative Trait Mapping. Genetics. 1994, 138 (3): 963-971.PubMedPubMed CentralGoogle Scholar
- Voorrips RE: MapChart: software for the Graphical Presentation of Linkage Maps and QTLs. J Hered. 2002, 93 (1): 77-78. 10.1093/jhered/93.1.77.PubMedGoogle Scholar
- Jaccard P: Distribution de la flore alpine dans le bassin des Dranses et dans quelques régions voisines. Bulletin de la Société Vaudoise des Sciences Naturelles. 1901, 37: 241-272.Google Scholar
- Zweifel R, Zimmermann L, Zeugin F, Newbery DM: Intra-annual radial growth and water relations of trees: implications towards a growth mechanism. J Exp Bot. 2006, 57 (6): 1445-1459. 10.1093/jxb/erj125.PubMedGoogle Scholar
- Craven D, Dent D, Braden D, Ashton MS, Berlyn GP, Hall JS: Seasonal variability of photosynthetic characteristics influences growth of eight tropical tree species at two sites with contrasting precipitation in Panama. For Ecol Manage. 2011, 261 (10): 1643-1653. 10.1016/j.foreco.2010.09.017.Google Scholar
- Gutierrez E, Campelo F, Camarero JJ, Ribas M, Muntan E, Nabais C, Freitas H: Climate controls act at different scales on the seasonal pattern of Quercus ilex L. stem radial increments in NE Spain. Trees-Struct Funct. 2011, 25 (4): 637-646. 10.1007/s00468-011-0540-3.Google Scholar
- Binkley D, Stape JL, Ryan MG, Barnard HR, Fownes J: Age-related decline in forest ecosystem growth: an individual-tree, stand-structure hypothesis. Ecosystems. 2002, 5 (1): 58-67. 10.1007/s10021-001-0055-7.Google Scholar
- Ryan MG, Binkley D, Fownes JH, Giardina CP, Senock RS: An experimental test of the causes of forest growth decline with stand age. Ecol Monogr. 2004, 74 (3): 393-414. 10.1890/03-4037.Google Scholar
- Ryan MG, Binkley D, Stape JL: Why don’t our stands grow even faster? Control of production and carbon cycling in eucalypt plantations. Southern Forests: a Journal of Forest Science. 2008, 71 (1): 99-104(106).Google Scholar
- Marques OG Jr, Andrade HB, Ramalho MAP: Assessment of the early selection efficiency in Eucalyptus cloeziana F. Muell. in the northwest of Minas Gerais state (Brazil). Silvae genetica. 1996, 45 (5/6): 359-361.Google Scholar
- Wei X, Borralho NMG: Genetic control of growth traits of Eucalyptus urophylla S.T. Blake in south east China. Silvae genetica. 1998, 47 (2-3): 158-165.Google Scholar
- Osorio LF, White TL, Huber DA: Age-age and trait-trait correlations for Eucalyptus grandis Hill ex Maiden and their implications for optimal selection age and design of clonal trials. Theor Appl Genet. 2003, 106 (4): 735-743.PubMedGoogle Scholar
- Osorio LF, White TL, Huber DA: Age trends of heritabilities and genotype-by-environment interactions for growth traits and wood density from clonal trials of Eucalyptus grandis Hill ex Maiden. Silvae genetica. 2001, 50 (1): 30-37.Google Scholar
- Kien ND, Jansson G, Harwood C, Thinh HH: Genetic control of growth and form in Eucalyptus urophylla in Northern Vietnam. J Trop For Sci. 2009, 21 (1): 50-56.Google Scholar
- Costa P, Durel CE: Time trends in genetic control over height and diameter in maritime pine. Can J Forest Res. 1996, 26 (7): 1209-1217. 10.1139/x26-135.Google Scholar
- Parisi LM: Shoot elongation patterns and genetic control of second-year height growth in Pinus taeda L. using clonally replicated trials. Gainesville: University of Florida: 2006.Google Scholar
- Bouvet JM, Vigneron P, Saya A: Phenotypic plasticity of growth trajectory and ontogenic allometry in response to density for Eucalyptus hybrid clones and families. Ann Bot. 2005, 96 (5): 811-821. 10.1093/aob/mci231.PubMedPubMed CentralGoogle Scholar
- Henery ML, Moran GF, Wallis IR, Foley WJ: Identification of quantitative trait loci influencing foliar concentrations of terpenes and formylated phloroglucinol compounds in Eucalyptus nitens. New Phytol. 2007, 176 (1): 82-95. 10.1111/j.1469-8137.2007.02159.x.PubMedGoogle Scholar
- Shepherd M, Kasem S, Lee D, Henry R: Mapping species differences for adventitious rooting in a Corymbia torelliana × Corymbia citriodora subspecies variegata hybrid. Tree Genetics & Genomes. 2008, 4 (4): 715-725. 10.1007/s11295-008-0145-1.Google Scholar
- Freeman JS, Potts BM, Vaillancourt RE: Few Mendelian genes underlie the quantitative response of a forest tree, Eucalyptus globulus, to a natural fungal epidemic. Genetics. 2008, 178 (1): 563-571. 10.1534/genetics.107.081414.PubMedPubMed CentralGoogle Scholar
- Morton NE, MacLean CJ: Analysis of family resemblance. 3. Complex segregation of quantitative traits. Am J Hum Gen. 1974, 26 (4): 489-503.Google Scholar
- Beavis WD: QTL analyses: power, precision, and accuracy. Molecular dissection of complex traits. 1998, 1998: 145-162.Google Scholar
- Markussen T, Fladung M, Achere V, Favre JM, Faivre-Rampant P, Aragones A, Perez DD, Harvengt L, Espinel S, Ritter E: Identification of QTLs controlling growth, chemical and physical wood property traits in Pinus pinaster. Silvae genetica. 2003, 52 (1): 8-15.Google Scholar
- Wu R, Bradshaw HD, Stettler RF: Developmental quantitative genetics of growth in Populus. Theor Appl Genet. 1998, 97 (7): 1110-1119. 10.1007/s001220050998.Google Scholar
- Emebiri LC, Devey ME, Matheson AC, Slee MU: Age-related changes in the expression of QTLs for growth in radiata pine seedlings. TAG Theor Appl Genet. 1998, 97 (7): 1053-1061. 10.1007/s001220050991.Google Scholar
- Lerceteau E, Szmidt AE, Andersson B: Detection of quantitative trait loci in Pinus sylvestris L. across years. Euphytica. 2001, 121 (2): 117-122. 10.1023/A:1012076825293.Google Scholar
- Emebiri LC: Detection and genetic mapping of quantitative trait loci influencing stem growth efficiency in radiata pine. Australia: Australian National University: 1997.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.