Genetic dissection of growth, wood basic density and gene expression in interspecific backcrosses of Eucalyptus grandis and E. urophylla
© Kullan et al.; licensee BioMed Central Ltd. 2012
Received: 5 January 2012
Accepted: 20 July 2012
Published: 20 July 2012
F1 hybrid clones of Eucalyptus grandis and E. urophylla are widely grown for pulp and paper production in tropical and subtropical regions. Volume growth and wood quality are priority objectives in Eucalyptus tree improvement. The molecular basis of quantitative variation and trait expression in eucalypt hybrids, however, remains largely unknown. The recent availability of a draft genome sequence (http://www.phytozome.net) and genome-wide genotyping platforms, combined with high levels of genetic variation and high linkage disequilibrium in hybrid crosses, greatly facilitate the detection of quantitative trait loci (QTLs) as well as underlying candidate genes for growth and wood property traits. In this study, we used Diversity Arrays Technology markers to assess the genetic architecture of volume growth (diameter at breast height, DBH) and wood basic density in four-year-old progeny of an interspecific backcross pedigree of E. grandis and E. urophylla. In addition, we used Illumina RNA-Seq expression profiling in the E. urophylla backcross family to identify cis- and trans-acting polymorphisms (eQTLs) affecting transcript abundance of genes underlying QTLs for wood basic density.
A total of five QTLs for DBH and 12 for wood basic density were identified in the two backcross families. Individual QTLs for DBH and wood basic density explained 3.1 to 12.2% of phenotypic variation. Candidate genes underlying QTLs for wood basic density on linkage groups 8 and 9 were found to share trans-acting eQTLs located on linkage groups 4 and 10, which in turn coincided with QTLs for wood basic density suggesting that these QTLs represent segregating components of an underlying transcriptional network.
This is the first demonstration of the use of next-generation expression profiling to quantify transcript abundance in a segregating tree population and identify candidate genes potentially affecting wood property variation. The QTLs identified in this study provide a resource for identifying candidate genes and developing molecular markers for marker-assisted breeding of volume growth and wood basic density. Our results suggest that integrated analysis of transcript and trait variation in eucalypt hybrids can be used to dissect the molecular basis of quantitative variation in wood property traits.
Eucalyptus tree species and hybrids are now the most widely planted hardwoods in tropical, subtropical and temperate regions, primarily due to their fast growth, short rotation, high productivity, adaptability to a broad range of environments and suitability for pulp and paper production [1, 2]. Eucalypt plantations currently occupy more than 20 million hectares in over 50 countries. India (3.9 million hectares), Brazil (3.7 million hectares) and China (2.6 million hectares) are the leading eucalypt growers in the world (http://www.git-forestry.com). Commonly planted eucalypts are mainly from the subgenus Symphyomyrtus, sections Latoangulatae (E. grandis E. urophylla), Maidenaria (E. globulus, E. nitens) and Exsertaria (E. camaldulensis, E. tereticornis) and include hybrids of some of these species .
Eucalyptus breeding programs mainly focus on primary growth traits, such as height and volume growth (diameter at breast height, DBH) and wood quality traits including physical properties such as wood density and chemical properties such as cellulose and lignin content. Wood basic density (oven dry mass per green volume) is an important trait for kraft pulp production  as it affects specific wood consumption. The genetics of growth and wood property traits has been studied in Eucalyptus[3–5] and the heritability of growth traits is generally reported to be lower than that of wood quality traits [6, 7] due to the large numbers of genes involved and environmental effects impacting on growth. The efficiency of selection for these traits could be enhanced by molecular breeding approaches enabled by the availability of high-throughput, genome-wide genotyping technologies [8, 9] and the recent completion of a draft reference genome sequence for Eucalyptus (E. grandis V1.0, JGI, http://www.phytozome.net).
The identification of genetic factors underlying quantitative variation for growth and wood properties is important for tree breeding as well as gene discovery efforts. QTL analysis in structured pedigrees allows the identification of genomic regions harbouring candidate genes and trait-linked molecular markers . Over the past decades, with the advancement of DNA marker technologies and genome-wide linkage mapping, major efforts have been dedicated to the genetic dissection of growth and wood quality traits in eucalypts [11–18]. The number of QTLs detected for growth traits in Eucalyptus has generally been lower (1–3) than that detected for wood quality traits (3–7), which may reflect the lower heritability associated with growth traits and the lack of statistical power to detect many small-effect growth QTLs segregating in tree breeding pedigrees. In addition, the number of QTLs have, most likely, been underestimated and the magnitude of QTL effects overestimated due to the relatively small progeny sample sizes used for QTL detection (100 to 200 individuals) . QTL intervals normally span 10–30 cM regions, which may contain thousands of genes in Eucalyptus[20, 21]. The identification of positional candidate genes by linkage mapping in trees therefore remains a difficult task.
To partly overcome the limited resolution of QTL analysis, genetical genomics approaches have been used to narrow down the list of candidate genes in QTL intervals in several plant species including eucalypts [13, 22–25]. In these approaches, the transcript levels of individual genes, generally obtained from microarray analyses, are treated as separate quantitative phenotypes and chromosomal regions affecting transcript variation (i.e. expression QTLs or eQTLs) are identified using conventional QTL detection methods. Comparison of trait QTL and eQTL positions enables the identification of candidate genes, which may potentially affect phenotypic traits through variation in transcript abundance and downstream effects on protein abundance and activity . eQTL analyses performed in Arabidopsis[22, 27–29], maize , wheat  and Populus[32, 33], have allowed elucidation of the genetic control of biochemical pathways as well as genetic correlations observed in populations. In Eucalyptus, Kirst et al. identified QTLs for transcript levels (i.e., eQTLs) of lignin-associated genes, which co-localized with QTLs for growth, suggesting that the same genomic regions may affect lignin content and tree growth.
Eucalyptus grandis is extensively used for the production of pulp in subtropical regions because of its rapid growth, however the species suffers from fungal diseases in tropical regions . Eucalyptus urophylla, a tropical eucalypt found in natural forests in Indonesia, is more disease tolerant than E. grandis. Interspecific hybrids of these two species combine the fast growth and successful vegetative propagation of E. grandis with the greater tolerance to fungal disease, excellent adaptability, higher wood density and greater coppicing capability of E. urophylla[35, 36]. The identification of genetic factors that differentiate these two species and that may underlie hybrid superiority are therefore important for Eucalyptus breeding programs and for fundamental understanding of growth and development in hybrids. Currently in Eucalyptus, advanced generation hybrid breeding and gene pyramiding are hampered by long generation times, hybrid breakdown, inbreeding depression and limited knowledge of genetic factors controlling interspecific trait variation. Knowledge of gene expression (eQTL) and trait QTLs associated with the trait variation, their main effects and interactions with other QTLs will facilitate marker-assisted selection in E. grandis or E. urophylla parents (for F1 hybrid breeding), or in hybrid progeny (for advanced generation hybrid breeding).
The aim of this study was to dissect the molecular genetic basis of diameter at breast height (DBH) and wood basic density in an interspecific backcross pedigree of E. grandis x E. urophylla. Herein we report QTLs for both traits and we demonstrate the potential use of expression QTL (eQTL) analysis for dissecting cis- and trans-acting components of genetic variation underlying major QTLs for wood basic density.
Trait distribution and correlation
Summary statistics of DBH and wood basic density measured at age four years in the E. grandis and E. urophylla backcross families
E. grandisBC family
E. urophyllaBC family
Framework maps for QTL analysis
For QTL analysis, a framework of testcross (1:1) DArT markers were selected from the previously constructed, high-density genetic maps . Using a minimum interval support of LOD 3.0, framework genetic linkage maps were produced for the E. grandis and E. urophylla BC parents and the F1 hybrid parent (Additional file 1: Table S1 and Additional file 2: Table S2). The average marker interval of the parental framework maps ranged from 6.1 cM (F1 hybrid in the E. grandis BC family) to 8.5 cM (E. urophylla BC parent). On average, 92 to 96% of loci in the four parental framework maps were within 10 cM of a neighbouring marker. The linkage group numbering and the orientation of the linkage groups were in accordance with Brondani et al..
QTL analysis of DBH and wood density
Putative QTLs for DBH at age four years identified by CIM in the E. grandis and E. urophylla backcross mapping families
Variance (R2) explained by QTL (%)
Additive effect (SD)a
E. grandis BC parent
F1 hybrid parent (E. grandis BC family)
F1 hybrid parent (E. urophylla BC family)
E. urophylla BC parent
Putative QTLs for wood density at four years identified by CIM in the E. grandis and E. urophylla backcross mapping families
Variance (R2) explained by QTL (%)
Additive effect (SD)a
F1 hybrid parent (E. grandis BC family)
F1 hybrid parent (E. urophylla BC family)
E . urophylla BC parent
Between one and eight QTLs were identified for wood density (Table 3). Only one QTL for density was identified in the E. urophylla BC parent and no QTLs were identified in the E. grandis BC parent. The majority of the density QTLs were identified in the F1 hybrid maps, with three and eight QTLs detected in the E. grandis and E. urophylla BC families, respectively (Table 3, Additional file 3: Figure S1). The percentage of phenotypic variation explained by the density QTLs ranged from 3.1 to 12.2% (Table 3). The positive effects of all three QTLs identified for density in the F1 hybrid map of the E. grandis BC family were associated with the E. urophylla allele. For the F1 hybrid map of the E. urophylla BC family, positive effects of QTLs on LG4, 9 and 10 were associated with the E. grandis allele, while those on LG2, 3, 6 and 8 were associated with the E. urophylla allele. Analysis of epistatic interactions among DBH and wood density QTLs revealed two significant interactions, one between LG2 and LG8, one between LG8 and LG10 for wood density QTLs in the F1 hybrid map of the E. urophylla BC family (Additional file 4: Figure S2).
Transcriptome profiling and expression QTL mapping of genes underlying a major QTL interval for wood density
The top ten positional candidate genes located in the wood density QTL interval on linkage group 9 of the F 1 hybrid parent in the E. urophylla backcross family (52.2 cM, 28,269,541 bp to 70.9 cM, 37,191,970 bp, E. grandis genome assembly V1.0, http://www.phytozome.net/ )
Physical position (bp)
eQTL on linkage group
eQTL position (cM)
Correlation with wood density (R2)
Cyclic nucleotide-regulated ion channel family protein
4, 6, 10
12, 101, 51
Histidine kinase 1
4, 9, 10
27, 52, 74
Auxin efflux carrier family protein
Cytochrome BC1 synthesis
Calcium-binding EF-hand family protein
F-box family protein
S-adenosyl-L-methionine-dependent methyltransferases superfamily protein
Glutaredoxin family protein
While most previous QTL mapping studies of diameter growth and wood density in Eucalyptus[12, 16–18, 40, 41] were based on F1 hybrid pedigrees, the present study focused on trait dissection in an F2 interspecific backcross pedigree. QTL mapping in the shared F1 hybrid parent and the E. grandis and E. urophylla BC parents allowed assessment of the architecture of interspecific as well as intraspecific genetic variation affecting trait variation. In this approach, fixed genetic differences between the parental species are likely to be in heterozygous state in the F1 hybrid and segregate in either or both backcross families depending on the degree of dominance. In addition, genetic factors that are heterozygous in the backcross parents (i.e. intraspecific variation) also segregate in the backcross progeny. If fixed genetic differences between the pure-species parents were purely due to additive genetic effects, the majority of QTLs in the F1 hybrid would segregate in both backcross families. However, the majority of QTLs in the two F1 hybrid maps were detected in only one of the two backcross families. Only one QTL for DBH on LG6 in the F1 hybrid parent was shared in both backcross families (Additional file 3: Figure S1). Failure to detect QTLs segregating from the F1 hybrid in both backcross families may be due to dominance effects (Additional file 8: Figure S4) playing a significant role in the expression of QTLs in alternative genetic backgrounds [42–44], or may be the result of epistatic interactions , or due to differences in QTL effects for the same alleles segregating in different genetic backgrounds (i.e. in the presence of a different set of segregating QTLs). For example, we identified a significant epistatic interaction between the wood density QTLs on LG8 and LG10, and between the wood density QTLs on LG2 and LG8 in the F1 hybrid map (E. urophylla BC family). This may explain why we detected the wood density QTLs on LG2 and LG8 in the E. urophylla BC family only (Additional file 3: Figure S1).
Overall, more QTLs were identified in the F1 hybrid parent (14) than in the backcross parents (3) for DBH (three compared to one in the BC parents) and density (eleven compared to one in the respective backcross parent, Tables 2 and 3). The majority of the positive QTL effects for DBH in the F1 hybrid were associated with the E. grandis allele and most of the positive QTL effects for density were associated with the E. urophylla allele. This is congruent with the expected interspecific and intraspecific genetic variation segregating in the backcross families (Table 1, Figure 1). The number of QTLs detected for DBH in Eucalyptus has generally been lower than that detected for wood density (Additional file 9: Table S4). The lower number of QTLs identified for DBH in this study (Tables 2 and 3) is consistent with published QTL reports, reflecting the lower heritability associated with growth traits compared to wood density in Eucalyptus and the limited statistical power to detect larger numbers of small effect QTLs. The well-described Beavis effect  certainly also applies in our study which means that some QTL effects listed in Tables 2 and 3 may be inflated and we fully expect that more QTLs of lower effect would be detected if our mapping populations were to be expanded.
DBH and density QTLs were detected in different regions of the genome (Additional file 3: Figure S1) suggesting that the two traits are affected by independent polymorphic loci in this pedigree. This is further supported by the low phenotypic correlation observed between DBH and wood density in the two backcross families (0.06 for E. grandis BC and −0.02 for E. urophylla BC). Freeman et al. also identified independent QTLs for DBH and wood density in an F2 outbred pedigree of E. globulus. However, other QTL studies in Eucalyptus have identified co-located QTLs affecting DBH and wood density in addition to independent QTLs [12, 40]. This could be explained by the occurrence of different polymorphisms affecting the two traits in each mapping pedigree resulting in different levels of correlation reported for DBH and density in previous studies [46–50]. The identification of independent QTLs for DBH and wood density in this study suggests that MAB could be used to improve growth and wood density simultaneously in this hybrid pedigree by selecting for combinations of QTL alleles with positive effects on DBH as well as wood density.
Comparative genetic mapping facilitates the identification of QTLs across different environments, ages and in different genetic backgrounds. Previous comparative genetic mapping studies in Eucalyptus suggested high levels of synteny and co-linearity among the genomes of eucalypt species [37, 38, 51, 52] enabling the comparative analysis of QTLs in different species [14, 16, 17, 40, 41, 51]. QTLs for DBH and wood density were detected on homologous linkage groups in the parental maps in this study (Additional file 3: Figure S1) and E. globulus and E. nitens linkage maps in previous studies (Additional file 10: Table S5 and Additional file 11: Table S6). QTLs identified for wood density on LG1, LG6 and LG10 of the F1 hybrid map may correspond to wood density QTLs previously identified on the same linkage groups in E. globulus, while QTLs identified for wood density on LG6 and LG8 of the F1 hybrid map may correspond to wood density QTLs identified in E. nitens. A QTL identified for wood density on LG9 of the F1 hybrid map (E. urophylla BC) may correspond to a wood density QTL previously identified in E. nitens. Similarly, QTLs identified for DBH on LG6 and LG10 (Table 2) of the same parental map may represent the same genomic regions as DBH QTLs reported in E. nitens and E. globulus, respectively. Common regions affecting trait variation across species should be the priority targets for the identification of candidate genes, the development of gene-based markers, association genetic studies and eventually MAB. We expect the resolution of comparative QTL analysis to drastically improve with the use of large numbers of trans-specific and trans-pedigree markers such as microsatellite, DArT and SNP markers linked to the E. grandis reference genome sequence.
An advantage of MAB in trees is the early selection of seedlings, reducing the time and cost normally involved in growing trees to maturity in the field before being able to identify elite trees . Experiments in crop plants have indicated that major effect QTLs and candidate genes associated with these QTLs are more reliable for MAB . Most of the QTLs detected in previous studies in Eucalyptus have likely been from the top end of the distribution of segregating QTL effects, some of which could be considered major effect QTLs [12, 14–17, 40]. However, there is a bottleneck between mapped QTLs and gene discovery mainly due to the low resolution of QTL mapping in populations of only several hundred individuals. To extend the information of QTL mapping, genetical genomic approaches  have been used to identify positional candidate genes and regulatory networks underlying phenotypic variation in several plant species [13, 22–24, 27]. In this study, the majority (70.8%, Additional file 6: Figure S3) of eQTLs identified for 294 xylem expressed genes underlying a major wood density QTL on LG9 (F1 hybrid map, Figure 2, Table 3), did not co-locate with the physical positions of the genes (i.e. were trans-acting eQTLs). The trans-eQTLs detected for these genes most likely correspond to diverse regulatory factors controlling the expression of the genes located in the QTL interval on LG9 one (or more) of which could harbor the trait altering polymorphism underlying the wood density QTL. eQTLs co-localizing with the physical genome position of the gene (cis-eQTL; 29.2%) were identified for only 86 genes in the interval, which is in agreement with the lower proportion of cis-eQTL previously reported for Eucalyptus (22%) in an interspecific backcross population of E. grandis and E. globulus and more recently for Populus (23%) using whole-genome microarray analysis in an interspecific hybrid pedigree .
Schadt et al. reported that genes whose transcript levels are correlated with trait variation could be considered potential candidate genes for the trait. In the present study the transcript levels of ten genes located in the wood density QTL interval on LG9 were found to be positively correlated (R2 > 0.4) with wood density variation (Table 4). Some of these genes encoding a nucleotide-regulated ion channel family protein (DND1), histidine kinase (HK), S-adenosyl-L-methionine-dependent methyltransferase (SAM), an auxin efflux carrier family protein and calcium-binding EF-hand family protein, have previously been reported to be involved in plant growth and development and cell wall biogenesis. The transcript abundance of three of these genes (DND1 HK1 and a gene encoding a calcium binding EF-hand family protein) was affected by trans-acting eQTLs on LG10 (Table 4). Importantly, the same genomic region on LG10 (51 cM to 74 cM) co-localized with a major wood density QTL on the same linkage group suggesting that this genomic region may contain trans-acting factors affecting wood density as well as the transcript abundance of the candidate genes underlying the wood density QTL on LG9 (Figure 2). DND1 has been shown to be involved in plant defense responses in Arabidopsis. HK was reported to act as a cytokinin receptor  involved in diverse plant growth and developmental processes [57, 58]. SAM is a key enzyme for the phenylpropanoid pathway, involved in the synthesis of lignin . Auxin, essential for plant growth and development (e.g. vascular tissue differentiation, apical development, cell elongation and tropical growth) is transported from cell to cell by auxin efflux carrier proteins [60–62]. Besseau et al. showed that a reduction in the level of hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferase (HCT), a gene involved in lignin biosynthesis, was correlated with the inhibition of auxin transport in Arabidopsis, suggesting that auxin efflux carrier family proteins might be important for cell wall deposition and lignification. Similarly, plant cells contain large amounts of calcium in their cell walls and previous studies highlighted Ca2+ playing a role in secondary cell wall biosynthesis [64, 65].
The observation that the top most positively correlated genes (at the level of transcript abundance) in the wood density QTL interval on LG9 (Table 4) prominently shared trans-eQTLs on LG4 and LG10 (F1 hybrid map, E. urophylla BC family) suggested the presence of trans-acting factors that also underlie wood density QTLs at the same loci. This, together with the detection of a significant epistatic interaction between wood density QTLs on LG8 and LG10, led us to investigate the transcript abundance of genes in the QTL interval on LG8, with the hypothesis that a similar cis-trans relationship would exist between LG8 and LG10 as was observed for LG9 and LG10. We indeed found that the top most positively correlated genes in the wood density QTL interval on LG8 also shared trans-eQTLs on LG4 and LG10 (and LG6, Additional file 3: Table S3a). The top most negatively correlated genes in the wood density QTL intervals on LG 8 and 9 (Additional file 7: Table S3b) did not exhibit such a strong pattern of shared trans-acting eQTLs, but it is formally possible that any of the positively or negatively correlated genes in these two QTL intervals affect trait abundance via a cis-acting and/or trans-acting eQTLs. Together, these findings suggest that at least some of the wood density eQTLs detected in this study may represent segregating components of a transcriptional network (Figure 3). Furthermore, our results suggest that transacting genes (e.g. transcription factors) located in the QTL intervals on LG4 and LG10, together with putative target genes located in the QTL intervals on LG8, LG9 and other identified wood density QTLs should be prioritized for further investigation. Trans-acting factors for which the parental species are differentiated would be heterozygous in the F1 hybrid and could have large effect on gene expression and trait variation in backcross progeny. Cases where transcription factors as well as their target genes segregate may give rise to detectable epistatic interactions as was putatively observed for the wood density QTLs on LG10 (trans-acting) and LG8 (cis-acting).
We have detected QTLs for DBH and wood density in an interspecific backcross pedigree of E. grandis x E urophylla, with each QTL explaining between 3.1 and 12.2% of the phenotypic variation. Furthermore, our study is the first to use Illumina RNA-Seq expression profiling in a segregating tree population to quantify transcript abundance and map eQTLs for candidate genes potentially affecting wood property variation. This approach allowed us to detect cis- and trans-eQTLs for candidate genes co-locating with four wood density QTLs generating a hypothesis for the underlying mode of action of the trait QTLs and suggesting the presence of a transcriptional network of which some components may segregate in this backcross pedigree (Figure 3). Candidate gene-based markers developed from these QTL and eQTL regions will promote MAB of hybrids of E. grandis and E. urophylla and allow more detailed molecular genetic dissection of quantitative trait variation in these trees.
Two previously described  interspecific F2 backcross families sharing an E. grandis x E. urophylla F1 hybrid parent (GUSAP1, Sappi, South Africa) were analysed to detect QTLs for DBH and wood basic density in the F1 hybrid and two pure-species, backcross parents. The F1 hybrid was originally derived from an E. grandis seed parent (GSAP1) pollinated with an E. urophylla pollen mix (parent unknown). A single F1 individual (GUSAP1) was selected for backcrossing to unrelated individuals of E. grandis (GSAP2) and E. urophylla (USAP1) to avoid possible inbreeding depression in the backcross progeny. A total of 308 seedlings of the E. urophylla BC family were planted in 2005 and 286 seedlings of the E. grandis BC family were planted in 2006 near KwaMbonambi, KwaZulu-Natal (Sappi, South Africa). The site chosen for the trial is located on flat coastal land with deep sandy soils and little spatial variation. The trees were planted in 2.7 x 2.2 m spacing and standard silvicultural operations were applied to ensure that the site was weed free at planting and remained weed free until canopy closure to ensure that the trees were not subjected to any competing vegetation.
Diameter (cm) at breast height (DBH) of the main stem was assessed in the E. urophylla and E. grandis BC families at age four in 2009 and 2010, respectively. For the assessment of wood property traits, a wood disk taken at the height of 1.35 m on the main stem was used to determine wood basic density using the water displacement method .
DNA isolation and DArT genotyping
High-throughput DNA extraction and DArT genotyping were used to obtain molecular marker genotypes for the three parental trees and all of the backcross progeny. The Eucalyptus DArT array was the same as that developed by Sansaloni et al.  and subsequently used for genotyping in other studies [37, 52]. The array comprises 7680 informative, polymorphic DArT markers selected by generating genomic representations from diverse Eucalyptus species and performing segregation analyses of more than 20,000 DArT polymorphisms in various Eucalyptus mapping populations .
Genetic mapping and QTL analysis
High-density, single-tree, genetic linkage maps were previously constructed for the E. grandis BC parent, the E. urophylla BC parent and the shared F1 hybrid parent using the two-way pseudo-testcross mapping strategy [37, 67]. For QTL analysis, particularly for Composite Interval Mapping (CIM), it is desirable to have uniform spacing of marker loci  and there is not much to gain from having markers closer than 5–10 cM . Therefore a subset of testcross DArT markers (1:1) with approximately 5–10 cM spacing were selected for the construction of framework linkage maps for the individual parental trees (Additional file 1: Table S1 and Additional file 2: Table S2). The map distances for the framework maps were recalculated using the Kosambi mapping function in Joinmap® 4 . The parameter settings used were Rec = 0.40, LOD = 3.0 and Jump = 5.
Marker data in the F1 hybrid parent map were recoded (1-E. grandis, 0-E. urophylla), based on the genotype of the E. grandis parent of the F1 hybrid  so that the QTL effect at each locus could be associated with the E. grandis or E. urophylla allele inherited from the F1 hybrid. The directions of QTL effects for the E. grandis and E. urophylla BC parents are arbitrary from one linkage group to another with respect to the grandparent alleles.
QTL analysis was conducted using CIM according to Zeng  using WinQTL Cartographer 2.5 . CIM was performed using Model 6 after scanning the genetic map and estimating the likelihood of a QTL and its corresponding effect at every 2 cM map interval, while using significant marker cofactors to adjust for the phenotypic effects associated with each target interval in the genetic map. Forward and backward stepwise regression, with a threshold of p <0.1, was used to select marker cofactors for background control with a window size of 10 cM on either side of the QTL test site. QTL peaks separated by a minimum distance of 20 cM or more on the same linkage group were considered evidence for two different QTLs . For peaks within 20 cM, the highest peak was chosen to approximate the position of the QTL. QTLs were declared significant by comparing LOD scores to an empirical genome-wide significance threshold calculated from 1,000 permutations for genome-wide α = 0.05 and α = 0.01 to control for type-I error. All QTLs with LOD support greater than 2.5 were also recorded. QTL position, LOD support, coefficients of determination (R2) and additive effects were estimated for each QTL. We used QTLNetwork (V2.0)  with a mixed-model composite interval mapping approach to test for epistatic interactions among the observed QTLs for DBH and wood density. The identified QTLs were projected onto linkage maps using MapChart 2.1 software .
Transcriptome profiling and eQTL analysis
The outer, differentiating xylem tissue from 96 individuals of the E. urophylla BC family was collected immediately after removing the bark from the zone of 1.5 to 2.0 m on the main stem. Tissue was immediately frozen in liquid nitrogen to avoid RNA degradation. The samples were stored at −80 °C until RNA extraction. Total RNA was isolated as described previously  and used for Illumina RNA-Seq analysis by BGI Americas (http://bgiamericas.com). A minimum of 20 million mapable paired-end reads (PE50) per individual was obtained. After quality filtering (per base Phred score > 20), on average 96% of the reads mapped (75% mapped as proper pairs, the rest as single-reads) to predicted gene models in the draft E. grandis genome sequence (DOE-JGI V1.0, http://www.phytozome.net) using TopHat version 1.30  and gene expression values (fragments per kilobase of coding sequence per million mapped fragments, FPKM) were calculated for each predicted gene model using Cufflinks version 1.0.3 (bias correction and quartile normalization was enabled for the FPKM calculation) . To identify putative positional candidate genes, eQTL analysis was performed for the genes underlying a major QTL interval for wood basic density on LG9 using CIM as described above for trait QTLs.
The authors are grateful to Sappi Forests Research (South Africa) for making the crosses, maintaining the mapping pedigree and phenotyping the individual trees. We thank in particular G. Galloway and W. Naidoo (Sappi) and their technical teams for assistance with phenotypic measurements. We acknowledge A. Kilian of Diversity Arrays Technology Pty Ltd (DArT P/L, Canberra, Australia) for assistance with marker analysis. This work was funded by Sappi through the Forest Molecular Genetics Program and by the Technology and Human Resources for Industry Program (THRIP), the National Research Foundation (NRF) and the Department of Science and Technology (DST) of South Africa. We thank the handling editor and two anonymous reviewers for helpful suggestions to improve the manuscript.
- Eldridge K, Davidson J, Harwood C, van Wyk G: Eucalypt domestication and breeding. 1993, Claredon Press, Oxford, 288-
- Myburg AA, Potts BM, Marques CM, Kirst M, Gion JM, Grattapaglia D, Grima-Pettenati J: Eucalyptus. Genome Mapping and Molecular Breeding in Plants. Volume 7. Edited by: Kole CR. 2007, Springer, Forest trees. New York, NY, USA, 115-160.
- Raymond CA, Apiolaza LA: Incorporating wood quality and deployment traits in Eucalyptus globules and Eucalyptus nitens. Plantation forest biotechnology for the 21st Century. Edited by: Walter C, Carson M. 2004, Forest Research New Zealand, Rotorua, New Zealand, 87-99.
- Stackpole DJ, Vaillancourt RE, Alves A, Rodrigues J, Potts BM: Genetic Variation in the chemical components of Eucalyptus globulus wood. G3. Genes, Genomes, Genetics. 2011, 1 (2): 151-
- Florestais P: Estimates of genetic parameters of wood traits for sawn timber production in Eucalyptus grandis. Genetics and Molecular Biology. 2004, 27 (4): 567-573. 10.1590/S1415-47572004000400017.View Article
- Raymond CA: Genetics of Eucalyptus wood properties. Annals of Forest Science. 2002, 59 (5–6): 525-531.View Article
- Hamilton MG, Potts BM: Review of Eucalyptus nitens genetic parameters. New Zealand Journal of Forestry Science. 2008, 38 (1): 102-119.
- Sansaloni CP, Petroli CD, Carling J, Hudson CJ, Steane DA, Myburg AA, Grattapaglia D, Vaillancourt RE, Kilian A: A high-density Diversity Arrays Technology (DArT) microarray for genome-wide genotyping in Eucalyptus. Plant Methods. 2010, 6 (1): 16-10.1186/1746-4811-6-16.PubMed CentralView ArticlePubMed
- Grattapaglia D, Resende MDV: Genomic selection in forest tree breeding. Tree Genetics & Genomes. 2011, 7 (2): 241-255. 10.1007/s11295-010-0328-4.View Article
- Grattapaglia D, Kirst M: Eucalyptus applied genomics: from gene sequences to breeding tools. New Phytologist. 2008, 179 (4): 911-929. 10.1111/j.1469-8137.2008.02503.x.View ArticlePubMed
- Byrne M, Murrell JC, Owen JV, Kriedemann P, Williams ER, Moran GF: Identification and mode of action of quantitative trait loci affecting seedling height and leaf area in Eucalyptus nitens. Theoretical and Applied Genetics. 1997, 94 (5): 674-681. 10.1007/s001220050465.View Article
- Grattapaglia D, Bertolucci FLG, Penchel R, Sederoff RR: 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 (3): 1205-PubMed CentralPubMed
- Kirst M, Myburg AA, De León JPG, Kirst ME, Scott J, Sederoff RR: Coordinated genetic regulation of growth and lignin revealed by quantitative trait locus analysis of cDNA microarray data in an interspecific backcross of Eucalyptus. Plant Physiology. 2004, 135 (4): 2368-10.1104/pp.103.037960.PubMed CentralView ArticlePubMed
- Thamarus K, Groom K, Bradley A, Raymond CA, Schimleck LR, Williams ER, Moran GF: Identification of quantitative trait loci for wood and fibre properties in two full-sib pedigrees of Eucalyptus globulus. Theoretical and Applied Genetics. 2004, 109 (4): 856-864. 10.1007/s00122-004-1699-4.View ArticlePubMed
- 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. Theoretical and Applied Genetics. 1997, 95 (4): 597-608. 10.1007/s001220050601.View Article
- 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.View Article
- Thumma BR, Southerton SG, Bell JC, Owen JV, Henery ML, Moran GF: Quantitative trait locus (QTL) analysis of wood quality traits in Eucalyptus nitens. Tree Genetics & Genomes. 2010, 6 (2): 305-317. 10.1007/s11295-009-0250-9.View Article
- Gion JM, Carouche A, Deweer S, Bedon F, Pichavant F, Charpentier JP, Bailleres H, Rozenberg P, Carocha V, Ognouabi N, Verhaegen D, Grima-Pettenati J, Vigneron P, Plomion C: Comprehensive genetic dissection of wood properties in a widely-grown tropical tree: Eucalyptus. BMC Genomics. 2011, 12: 301-10.1186/1471-2164-12-301.PubMed CentralView ArticlePubMed
- Beavis WD: QTL analyses: power, precision, and accuracy. Molecular Dissection of Complex Traits. Edited by: Paterson AHCRC, Boca Raton FL. 1998, 145-162.
- Grattapaglia D, Bradshaw HD: Nuclear DNA content of commercially important Eucalyptus species and hybrids. Canadian Journal of Forest Research. 1994, 24: 1074-1078. 10.1139/x94-142.View Article
- Grattapaglia D, Plomion C, Kirst M, Sederoff RR: Genomics of growth traits in forest trees. Current Opinion in Plant Biology. 2009, 12 (2): 148-156. 10.1016/j.pbi.2008.12.008.View ArticlePubMed
- DeCook R, Lall S, Nettleton D, Howell SH: Genetic regulation of gene expression during shoot development in Arabidopsis. Genetics. 2006, 172 (2): 1155-PubMed CentralView ArticlePubMed
- Druka A, Potokina E, Luo Z, Bonar N, Druka I, Zhang L, Marshall DF, Steffenson BJ, Close TJ, Wise RP: Exploiting regulatory variation to identify genes underlying quantitative resistance to the wheat stem rust pathogen Puccinia graminis f. sp. tritici in barley. Theoretical and Applied Genetics. 2008, 117 (2): 261-272. 10.1007/s00122-008-0771-x.View ArticlePubMed
- Kloosterman B, Oortwijn M, UitdeWilligen J, America T, de Vos R, Visser RGF, Bachem CWB: From QTL to candidate gene: Genetical genomics of simple and complex traits in potato using a pooling strategy. BMC Genomics. 2010, 11: 158-10.1186/1471-2164-11-158.PubMed CentralView ArticlePubMed
- Li Y, Breitling R, Jansen RC: Generalizing genetical genomics: getting added value from environmental perturbation. Trends in Genetics. 2008, 24 (10): 518-524. 10.1016/j.tig.2008.08.001.View ArticlePubMed
- Jansen RC, Nap JP: Genetical genomics: the added value from segregation. Trends in Genetics. 2001, 17 (7): 388-391. 10.1016/S0168-9525(01)02310-1.View ArticlePubMed
- Keurentjes JJB, Fu J, Terpstra IR, Garcia JM, Van Den Ackerveken G, Snoek LB, Peeters AJM, Vreugdenhil D, Koornneef M, Jansen RC: Regulatory network construction in Arabidopsis by using genome-wide gene expression quantitative trait loci. Proceedings of the National Academy of Sciences USA. 2007, 104 (5): 1708-10.1073/pnas.0610429104.View Article
- West MAL, Kim K, Kliebenstein DJ, Van Leeuwen H, Michelmore RW, Doerge R, St Clair DA: Global eQTL mapping reveals the complex genetic architecture of transcript-level variation in Arabidopsis. Genetics. 2007, 175 (3): 1441-PubMed CentralView ArticlePubMed
- Chan EKF, Rowe HC, Corwin JA, Joseph B, Kliebenstein DJ: Combining genome-wide association mapping and transcriptional networks to identify novel genes controlling glucosinolates in Arabidopsis thaliana. PLoS Biology. 2011, 9 (8): e1001125-10.1371/journal.pbio.1001125.PubMed CentralView ArticlePubMed
- Shi C, Uzarowska A, Ouzunova M, Landbeck M, Wenzel G, Lübberstedt T: Identification of candidate genes associated with cell wall digestibility and eQTL (expression quantitative trait loci) analysis in a Flint × Flint maize recombinant inbred line population. BMC Genomics. 2007, 8 (1): 22-10.1186/1471-2164-8-22.PubMed CentralView ArticlePubMed
- Jordan MC, Somers DJ, Banks TW: Identifying regions of the wheat genome controlling seed development by mapping expression quantitative trait loci†. Plant Biotechnology Journal. 2007, 5 (3): 442-453. 10.1111/j.1467-7652.2007.00253.x.View ArticlePubMed
- Street NR, Skogström O, Sjödin A, Tucker J, Rodríguez Acosta M, Nilsson P, Jansson S, Taylor G: The genetics and genomics of the drought response in Populus. The Plant Journal. 2006, 48 (3): 321-341. 10.1111/j.1365-313X.2006.02864.x.View ArticlePubMed
- Drost DR, Benedict CI, Berg A, Novaes E, Novaes CRDB, Yu Q, Dervinis C, Maia JM, Yap J, Miles B, Kirst M: Diversification in the genetic architecture of gene expression and transcriptional networks in organ differentiation of Populus. Proceedings of the National Academy of Sciences USA. 2010, 107 (18): 8492-10.1073/pnas.0914709107.View Article
- Wingfield MJ, Crous PW, Coutinho TA: A serious canker disease of Eucalyptus in South Africa caused by a new species of Coniothyrium. Mycopathologia. 1996, 136 (3): 139-145. 10.1007/BF00438919.View ArticlePubMed
- Vigneron P, Bouvet J: Eucalypt hybrid breeding in Congo. Hybrid Breeding and Genetics of Forest Trees. Proceedings of QFRI/CRC-SPF Symposium, 9–14th April 2000 Noosa, Queensland. Edited by: Dungey HS, Dieters MJ, Nikles DG. 2000, Australia. Department of Primary Industries, Brisbane, 14-26.
- Campinhos E, Ikemori YK: Selection and management of the basic population Eucalyptus grandis and E. urophylla established at Aracruz for the long term breeding programme. Proceedings of the IUFRO Conference: Breeding Tropical Trees: Population Structure and Genetic Improvement Strategies in Clonal and Seedling Forestry. Pattaya. Edited by: Gibson GI, Griffin AR, Matheson AC. 1989, Thailand. Oxford Forestry Institute, Oxford, 169-175.
- Kullan ARK, van Dyk MM, Jones N, Kanzler A, Bayley A, Myburg AA: High-density genetic linkage maps with over 2,400 sequence-anchored DArT markers for genetic dissection in an F2 pseudo-backcross of Eucalyptus grandis × E. urophylla. Tree Genetics & Genomes. 2011, 8 (1): 163-175.View Article
- Brondani RPV, Williams ER, 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 Biology. 2006, 6 (1): 20-10.1186/1471-2229-6-20.PubMed CentralView ArticlePubMed
- Schadt EE, Lamb J, Yang X, Zhu J, Edwards S, GuhaThakurta D, Sieberts SK, Monks S, Reitman M, Zhang C: An integrative genomics approach to infer causal associations between gene expression and disease. Nature Genetics. 2005, 37 (7): 710-10.1038/ng1589.PubMed CentralView ArticlePubMed
- 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.View Article
- 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.View Article
- Retief ECL, Stanger TK: Genetic parameters of pure and hybrid populations of Eucalyptus grandis and E. urophylla and implications for hybrid breeding strategy. Southern Forests: a Journal of Forest Science. 2009, 71 (2): 133-140.View Article
- Bouvet JM, Vigneron P: Variance structure in Eucalyptus hybrid populations. Silvae Genetica. 1996, 45: 171-176.
- De Assis TF, Rezende G, Aguiar AM: Current status of breeding and deployment for clonal forestry with tropical eucalypt hybrids in Brazil. International Forestry Review. 2005, 7: 61-
- True JR, Liu J, Stam LF, Zeng ZB, Laurie CC: Quantitative genetic analysis of divergence in male secondary sexual traits between Drosophila simulans and Drosophila mauritiana. Evolution. 1997, 51 (3): 816-832. 10.2307/2411157.View Article
- McDonald A, Borralho N, Potts B: Genetic variation for growth and wood density in Eucalyptus globulus ssp. globulus in Tasmania (Australia). Silvae Genetica. 1997, 46 (4): 236-241.
- Muneri A, Raymond CA: Genetic parameters and genotype-by-environment interactions for basic density, pilodyn penetration and stem diameter in Eucalyptus globulus. Forest Genetics. 2000, 7 (4): 317-328.
- López GA, Potts BM, Dutkowski GW, Apiolaza LA, Gelid P: Genetic variation and inter-trait correlations in Eucalyptus globulus base population trials in Argentina. Forest Genetics. 2002, 9 (3): 217-231.
- Harrand L, Hernandez JJV, Upton JL, Valverde GR: Genetic parameters of growth traits and wood density in Eucalyptus grandis progenies planted in Argentina. Silvae Genetica. 2009, 58 (1–2): 11-19.
- Costa e Silva J, Borralho NMG, Araújo JA, Vaillancourt RE, Potts BM: Genetic parameters for growth, wood density and pulp yield in Eucalyptus globulus. Tree Genetics & Genomes. 2009, 5 (2): 291-305. 10.1007/s11295-008-0174-9.View Article
- Marques C, Brondani R, Grattapaglia D, Sederoff R: Conservation and synteny of SSR loci and QTLs for vegetative propagation in four Eucalyptus species. Theoretical and Applied Genetics. 2002, 105 (2): 474-478. 10.1007/s00122-002-0899-z.PubMed
- Hudson CJ, Kullan ARK, Freeman JS, Faria DA, Grattapaglia D, Kilian A, Myburg AA, Potts BM, Vaillancourt RE: High synteny and colinearity among Eucalyptus genomes revealed by high-density comparative genetic mapping. Tree Genetics & Genomes. 2011, 8 (2): 339-352.View Article
- Xu Y, Crouch JH: Marker-assisted selection in plant breeding: from publications to practice. Crop Science. 2008, 48 (2): 391-10.2135/cropsci2007.04.0191.View Article
- Kirst M, Basten CJ, Myburg AA, Zeng ZB, Sederoff RR: Genetic architecture of transcript-level variation in differentiating xylem of a Eucalyptus hybrid. Genetics. 2005, 169 (4): 2295-2303. 10.1534/genetics.104.039198.PubMed CentralView ArticlePubMed
- Clough SJ, Fengler KA, Yu I: The Arabidopsis dnd1 “defense, no death” gene encodes a mutated cyclic nucleotide-gated ion channel. Proceedings of the National Academy of Sciences USA. 2000, 97 (16): 9323-10.1073/pnas.150005697.View Article
- Nishimura C, Ohashi Y, Sato S, Kato T, Tabata S, Ueguchi C: Histidine kinase homologs that act as cytokinin receptors possess overlapping functions in the regulation of shoot and root growth in Arabidopsis. The Plant Cell Online. 2004, 16 (6): 1365-1377. 10.1105/tpc.021477.View Article
- Mok DWS, Mok MC: Cytokinin metabolism and action. Annual Review of Plant Biology. 2001, 52 (1): 89-118. 10.1146/annurev.arplant.52.1.89.View Article
- Ueguchi C, Sato S, Kato T, Tabata S: The AHK4 gene involved in the cytokinin-signaling pathway as a direct receptor molecule in Arabidopsis thaliana. Plant and Cell Physiology. 2001, 42 (7): 751-755. 10.1093/pcp/pce094.View ArticlePubMed
- Weng JK, Chapple C: The origin and evolution of lignin biosynthesis. New Phytologist. 2010, 187 (2): 273-285. 10.1111/j.1469-8137.2010.03327.x.View ArticlePubMed
- Friml J, Yang X, Michniewicz M, Weijers D, Quint A, Tietz O, Benjamins R, Ouwerkerk PBF, Ljung K, Sandberg G: A PINOID-dependent binary switch in apical-basal PIN polar targeting directs auxin efflux. Science. 2004, 306 (5697): 862-865. 10.1126/science.1100618.View ArticlePubMed
- Taylor LP, Grotewold E: Flavonoids as developmental regulators. Current Opinion in Plant Biology. 2005, 8 (3): 317-323. 10.1016/j.pbi.2005.03.005.View ArticlePubMed
- Blilou I, Xu J, Wildwater M, Willemsen VA, Paponov I, Friml J, Heidstra R, Aida M, Palme K, Scheres BJG: The PIN auxin efflux facilitator network controls growth and patterning in Arabidopsis roots. Nature. 2005, 433: 39-44. 10.1038/nature03184.View ArticlePubMed
- Besseau S, Hoffmann L, Geoffroy P, Lapierre C, Pollet B, Legrand M: Flavonoid accumulation in Arabidopsis repressed in lignin synthesis affects auxin transport and plant growth. The Plant Cell Online. 2007, 19 (1): 148-162. 10.1105/tpc.106.044495.View Article
- Roberts AW, Haigler CH: Tracheary-element differentiation in suspension-cultured cells of Zinnia requires uptake of extracellular Ca2+. Planta. 1990, 180 (4): 502-509. 10.1007/BF02411447.View ArticlePubMed
- Kobayashi H, Fukuda H: Involvement of calmodulin and calmodulin-binding proteins in the differentiation of tracheary elements in Zinnia cells. Planta. 1994, 194 (3): 388-394.View Article
- TAPPI: Basic density and moisture content of pulp wood. 2006, TAPPI No. T258 om-06
- Grattapaglia D, Sederoff R: Genetic linkage maps of Eucalyptus grandis and Eucalyptus urophylla using a pseudotestcross mapping strategy and RAPD markers. Genetics. 1994, 137: 1121-1137.PubMed CentralPubMed
- Basten CJ, Weir BS, Zeng ZB: QTL Cartographer. 2001, Department of Statistics, North Carolina State University, Raleigh, NC
- Darvasi A, Soller M: Optimum spacing of genetic markers for determining linkage between marker loci and quantitative trait loci. Theoretical and Applied Genetics. 1994, 89 (2): 351-357.PubMed
- Van Ooijen JW: JoinMap 4. Software for the calculation of genetic linkage maps in experimental populations. 2006, Plant Research International Wageningen, The Netherlands
- Zeng ZB: Precision mapping of quantitative trait loci. Genetics. 1994, 136 (4): 1457-1468.PubMed CentralPubMed
- Wang S, Basten CJ, Zeng ZB: Windows QTL Cartographer 2.5. Department of Statistics, North Carolina State University. 2005, Raleigh, NC
- Paterson AH, Damon S, Hewitt JD, Zamir D, Rabinowitch HD, Lincoln SE, Lander ES, Tanksley SD: Mendelian factors underlying quantitative traits in tomato: comparison across species, generations, and environments. Genetics. 1991, 127 (1): 181-197.PubMed CentralPubMed
- Yang J, Hu CC, Ye XZ, Zhu J: QTLNetwork 2.0. 2005, Institute of Bioinformatics, Zhejiang University, Hangzhou, China, Available at http://ibi.zju.edu.cn/software/qtlnetwork.
- Voorrips RE: MapChart: software for the graphical presentation of linkage maps and QTLs. Journal of Heredity. 2002, 93 (1): 77-78. 10.1093/jhered/93.1.77.View ArticlePubMed
- Chang S, Puryear J, Cairney J: A simple and efficient method for isolating RNA from pine trees. Plant Molecular Biology Reporter. 1993, 11 (2): 113-116. 10.1007/BF02670468.View Article
- Trapnell C, Pachter L, Salzberg SL: TopHat: discovering splice junctions with RNA-Seq. Bioinformatics. 2009, 25 (9): 1105-10.1093/bioinformatics/btp120.PubMed CentralView ArticlePubMed
- Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, Van Baren MJ, Salzberg SL, Wold BJ, Pachter L: Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nature Biotechnology. 2010, 28 (5): 511-515. 10.1038/nbt.1621.PubMed CentralView ArticlePubMed
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.