Table of Contents

0878-B1

Quantitative Trait Loci (QTL) Mapping of Growth and Wood Property Traits in Populus

Zhang Deqiang[1], Zhang Zhiyi, Yang Kai and Li Bailian


Abstract

The genetic control of wood volume growth and wood chemical contents traits was studied in an interspecific backcross progeny derived from the cross between clone “LM50” (P. tomentosa) and its hybrid clone “TB01” (P. tomentosa, P. bolleana). Genetic maps were constructed for each parent with AFLP markers according to Pseudo-test-cross mapping strategy, based on 120 offspring. A total of 247 and 146 markers distributed among 19 linkage groups covered 3265.1 cM and 1992 cM for the clone “LM50” and clone “TB01” parents, respectively. Using genetic linkage maps constructed in this research, 32 putative QTLs were detected that are associated with the five main economic traits by internal mapping approach. They were dispersed on 16 linkage groups for clone “LM50” maps and 10 for clone “TB01” maps. The phenotypic variance explained by each QTL ranged from 7.0% to 14.6%. QTLs controlling correlated traits between sylleptic branch numbers and stem volume were co-localized in two genomic regions on linkage groups TLG6 and TLG8, respectively. The favorable alleles were derived from the parent P. tomentosa, which was phenotypically superior over parent P. bolleana for sylleptic branch angle, stem volume and wood chemical content traits, conversely for the sylleptic branch number trait. These preliminary results should be verified in larger population, different environmental backgrounds and years, and eventually the stable QTLs were detected for future development of marker-assisted selection for growth and wood properties in P. tomentosa.


Introduction

Populus tomentosa is a native species of main economic importance among the forest trees cultivated in China. It has played a key role in timber production and ecological environment protection along the Yellow River. Over past several years, genetic improvements in P. tomentosa by hybridization, chromosome doubling and clone breeding have been developed, and great progress has been achieved. Owing to poor research of genetic control system for this indigenous species, however, the existed problems for propagated rooting, yields of wood and quality have not been resolved in essence to date. These commercially economic traits are of a quantitative nature, and controlled by its corresponding quantitative trait loci (QTLs). The use of molecular markers enables one to identify and map QTLs that are involved in the variation of such traits in the segregated population. Knowledge of the number, genomic regions and individual effects of genes controlling these traits will provide a better understanding of the genetic architecture of quantitative traits and expedite breeding progress through opportunities for developing marker-assisted breeding strategies in forest trees.

For the last decade, QTLs mapping with molecular markers has been received extensive attention in forest trees. QTLs were detected in poplar for many important traits such as growth, form, phenology and disease resistance, based on F2 generation population of P. trichocarpa?P. deltoides (Bradshaw and Stettler 1995; Frewen et al. 2000). Tree growth, wood quality and vegetative propagation in Eucalyptus (Grattapaglia et al. 1996; Marques et al. 1999), growth, adaptation, wood specific gravity, wood physical and chemical properties in Pinus (Groover et al. 1994; Sewell et al. 2000, 2002) and many more quantitative traits have already been studied in other forest tree species. However, the QTLs controlling growth and wood chemical content traits remain unknown in P. tomentosa. To detect the QTLs affecting stem growth and wood chemical contents traits, the parent-specific moderate resolution maps with AFLP markers were constructed for two cultivars clone “LM50” (P. tomentosa) and clone “TB01” (P. tomentosa´P. bolleana) based on 120 progeny. In the present study, the separated AFLP genetic maps were used with objectives: (1) identify QTLs associated with five economic importance traits for sylleptic branch numbers, sylleptic branch angle, stem volume, wood cellulose content and wood lignin content; (2) estimate the number, genomic location and individual effect of gene affecting these traits, using the interval mapping approach. To our knowledge, this is the first report on QTL analysis based on an interspecific BC1 population in P. tomentosa.

Materials and methods

Mapping population

A subset of 120 progeny, chosen randomly from inbred population of 696 lines derived from a cross between two elite poplar species, “TB01” clone and “LM50”clone, was used in this research.

Linkage map construction

Separated genetic linkage maps have been constructed with AFLP markers according to Pseudo-test-cross mapping strategy (Grattapaglia and sedroff 1994) using the Kosambi’s mapping function.

QTL analysis

QTL detection was performed with software package MAPMAKER-QTL (Lander and Botstein 1989) underling the backcross model, based on the individual linkage maps for clone “TB01” and clone “LM50”.

Results

QTL analysis

Thirty-two putative QTLs associated with the 5 main economic traits (sylleptic angle, sylleptic number, stem volume, cellulose content and lignin content) were identified. These putative QTLs were dispersed on 16 linkage groups for clone “LM50” maps and 10 for clone “TB01” maps, respectively. The number of QTLs detected for each trait ranged from four to eight. The proportion of phenotypic variation explained by each of these QTLs varied from 7.0% (TBLG17) to 14.6% (TBLG11). More than 50.0% of total phenotypic variation was explained by the QTLs for each trait except for stem volume trait. Table 1 summarized the genomic location and associated additive effect of magnitude for all traits, together with their corresponding confidence interval and the percentage of the phenotypic variation explained by each QTL.

Sylleptic angle

Six QTLs detected were fell into 6 genomic regions corresponding, respectively, to linkage groups TLG2, TLG5, TLG7, TLG10, TLG16 and TBLG14 (Table 1). Among which 3 QTLs increased phenotypic values and another 3 have negative effects on the Sylleptic angles. The percentage of phenotypic variance explained by each QTLs ranged from 7.4% (TLG16) to 12.8% (TLG10). Total phenotypic variation explained by the six putative QTL was 57.7%.

Sylleptic number

A total of eight genomic regions were found to have effects on Sylleptic number (Table 1) and seven QTLs originated from both parents increased phenotypic value. The proportion of the total variance explained by a single QTLs varied between 7.0% (TBLG17) and 9.2% (TLG9). The total phenotypic variation explained by these seven putative QTLs was 61.7%, but individual QTL explained a small proportion of the phenotypic variance (less than 10.0%).

Stem volume

Four putative QTLs were identified that influenced stem volume (Table 1). The percentage of trait variation explained by each individual QTL ranged from 8.0% (TLG1) to 11.0% (TLG6) with most of them accounting for less than 10.0% of the total phenotypic variation. At these four loci, stem volume was increased by the alleles from or a closely linked QTL had additive effects between 0.936 and 1.198 m3 for the phenotypic value. QTLs controlling correlated traits between sylleptic branch numbers and stem volume were co-localized in two genomic regions on linkage groups TLG6 and TLG8, respectively.

Wood cellulose content

Eight unique quantitative trait loci affecting cellulose content were detected on linkage group TLG1, TLG5, TLG12, TLG13, TLG18, TBLG9, TBLG11 and TBLG14, respectively (Table 1). The respective putative QTLs accounted for between 8.7% (TLG13) and 14.6% (TBLG11) of phenotypic variation, with LOD scores ranging from 2.02 to 3.2. These eight QTLs could account for 83.7% of the total variation. One QTL of these was located in an identical region to the QTL for sylleptic number on TLG5 at interval of E35M52350r-E44M6097r.

Klason lignin content

For this trait, six putative QTLs were identified (Table 1). They were positioned on linkage group TLG2 and TLG15 in linkage map of “LM50” and on linkage group TBLG3, TBLG6, TBLG8 and TBLG5 in linkage map of “TB01”, respectively. Together these explained 54.9% of the additive genetic variation and all but one each explained less than 10.0% of the total phenotypic variation for the trait analyzed.

Discussion

QTL identification in inbred population

In order to understand the genetic control of growth and wood chemical components, backcrossing was used as a method to generate segregating population for QTL mapping. The use of genetic maps based on molecular markers enables one to identify and map QTLs that are involved in the variation of such traits in the segregated population. These could expedite efficiency of early selection of main economic traits through marker-assisted breeding strategies contrast to classical quantitative genetics. In this paper we identified 32 putative QTLs on 26 genomic regions controlling wood growth and chemical content traits. The number of QTLs identified ranged from 4 to 8 per trait and each QTL accounted only for a small proportion of the phenotypic variation (7.0%-14.6%). A previous QTL detection in the F2 progeny of P. trichocarpa ´P. deltoids was achieved with RFLP and RAPD markers by Bradshaw et al.(1995). Their QTL mapping results suggest that 24.4-33.4% phenotypic variation for growth and development trait was explained by single QTL. There are many reasons why no more than 14.6% of the phenotypic variation explained by single QTL for any trait was found in our study. The reason for this may be due to types and the sample size of mapping population, types and density of marker used in genetic maps, and QTL mapping software. First, the efficiency of QTL detection obtained from backcross is less than that obtained from F2 population because one, rather than two, recombinant gametes are sampled in plant. The small size of mapping population decreases the sensitivity of detecting QTLs, and only major QTLs could be detected, while smaller ones might remain unrevealed. Further, the dominant markers only used by us lead to lower LOD scores and also to an under-estimation of the magnitude of explained variation if multiple alleles are segregating at a locus in offspring of out-breeding species. Furthermore, our medium-density parentic maps constructed have several gaps and have not full Populus genome coverage. As a result, some QTLs of large effect have gone undetected. Therefore, There is much additive genetic variation unaccounted for which should require further field trials in order to resolve, or attribute to the QTLs in larger population.

Table 1 Putative QTLs detected for growth and wood chemical content traits by interval mapping via MapMaker/QTL.

Traits/QTL

Linkage group

Interval

Length a (cM)

Position b (cM)

LOD

Variance c %

Additive effect

SBAd

SBATL2

TLG2

E60M34268-E63M39436

16.0

2.0

2.74

11.5

-6.05

SBATL5

TLG5

E35M52350r-E44M6097r

16.7

16.2

2.0

7.6

-4.06

SBATL7

TLG7

E63M32396-E44M60413

11.2

3.9

2.02

7.8

+5.08

SBATL10

TLG10

E44M40100r-E34M47391r

12.1

12.0

3.1

12.8

+6.31

SBATL16

TLG16

E61M41145-E61M41144

4.7

1.3

2.0

7.4

+4.37

SBATBL14

TBLG14

E33M48277-E65M3181r

12.5

6.3

2.12

10.6

-5.82






Total

57.7


SBNe

SBNTL3

TLG3

E61M4152r-E63M33357

11.3

0.1

2.02

7.3

+8.76

SBNTL6

TLG6

E33M62439-E61M3192

30.1

0.1

2.05

7.6

+9.43

SBNTL8

TLG8

E60M33270r-E33M44290r

12.2

0.1

2.0

7.2

+8.06

SBNTL9

TLG9

E63M36218-E63M39248

16.9

15.2

2.15

9.2

+9.87

SBNTL14

TLG14

E35M52154-E61M3170

19.1

0.2

2.0

7.3

+8.76

SBNTL19

TLG19

E60M33227r-E33M40307

20.3

6.0

2.03

7.5


SBNTBL1

TBLG1

E33M44298-E34M44490r

11.1

5.4

2.1

8.6

+8.89

SBNTBL17

TBLG17

E63M52396-E60M3267

22.0

2.9

2.0

7.0

-9.55






Total

61.7

+8.16

SVf

SVTL1

TLG1

E33M32273-E33M32258r

6.2

3.5

2.0

8.0


SVTL6

TLG6

E33M62439-E61M3192

30.1

0.2

2.16

11.0

+0.936

SVTL8

TLG8

E60M33270r-E33M44290r

12.2

1.4

2.04

8.2

+0.940

SVTBLG2

TBLG2

E44M46517-E3342134r

33.0

7.2

2.08

8.6

+1.198






Total

35.8

+0.950

WCCg

WCCTL1

TLG1

E61M31154-E61M31157

8.0

0.7

2.05

9.0


WCCTL5

TLG5

E35M52350r-E44M6097r

16.7

6.8

2.15

9.8

+1.48

WCCTL12

TLG12

E34M48492r-E44M44664

25.3

0.4

2.63

10.1

+1.60

WCCTL13

TLG13

E63M36140r-E3338126r

36.7

12.6

2.02

8.7

-1.63

WCCTL18

TLG18

E63M3396-E65M3474

17.1

9.1

2.92

12.9

-1.44

WCCTBL9

TBLG9

E33M79222-E63M36283

28.0

4.0

2.10

9.4

+1.83

WCCTBL11

TBLG11

E33M79660-E44M40293

40.6

9.2

3.20

14.6

-1.56

WCCTBL14

TBLG14

E35M50220-E60M34338

26.7

23.1

2.08

9.2

+2.41






Total

83.7

+1.52

WLCh

WLCTL2

TLG2

E63M39436-E35M52185

16.0

0.4

2.06

9.2


WLCTL15

TLG15

E33M32134-E63M6198

24.9

1.7

2.20

10.0

+1.52

WLCTBL3

TBLG3

E63M3672-E61M31125

13.9

0.3

2.04

8.8

-1.64

WLCTBL5

TBLG5

E33M82207-E44M60644

12.4

12.3

2.00

8.5

+1.46

WLCTBL6

TBLG6

E33M40134-E65M34504

14.3

1.2

2.18

9.8

-1.60

WLCTBL8

TBLG8

E33M4266-E33M41151

35.3

17.7

2.02

8.6

+1.44






Total

54.9

+1.42

aInterval between the two flanking markers (cM) where a QTL is located
b QTL position from the first marker (cM)
c Phenotypic variation explained by each QTL
dSylleptic brance angle
eSylleptic brance number
fStem volume
gWood cellulose content
hwood lignin content

Trait correlation and pleiotropic effects of QTLs

Conventional quantitative genetics suggest that trait correlations may be attributable to either pleiotropic effects of single genes or to tight linkage of several genes that individually influence specific traits. Thus, it is expected that correlated traits have QTLs mapping to the same genomic regions (Veldboom et al. 1994; Xiao et al. 1996). The same trend was observed in our study. A significant positive correlation was found between sylleptic branch number and stem volume, and two QTLs (on linkage group TLG6 and TLG8, respectively) affecting both traits were mapped to two same genomic regions. The presence of QTLs for both sylleptic number and stem volume are detected in the same genomic regions may be due to the presence of two genes tightly linked or to one gene with a pleiotrophic effect on these two traits. In a previous study by Bradshaw et al.(1995), pleiotropy was suggested for QTLs on linkage group E that were simultaneously associated with spring flush, stem height, basal area and sylleptic branch number. In fact, it is observed that the close trait correlations cannot be explained by QTLs mapped in the same genomic region, or non-significant correlations were found to locate in same position were been observed in some cases (Damerval et al. 1994; Wu et al. 1997). For example, in this study, the QTLs located on same position were found to associate with sylleptic angle and cellulose content although non-correlation both traits have been observed. The real reasons for this phenomenon have not well understood but are generally believed to include that some genomic regions controlling these phenotypic traits have not been detected due to not saturated linkage maps and small experimental population or genetic factors such as repulsion and epistasis, and statistical methods or sampling error.

References

Bradshaw HD, Stettler RF (1995) Molecular genetics of growth and development in Populus. IV. Mapping QTLs with large effects on growth, form, and phenology traits in a forest tree. Genetics 139: 963-973

Damerval C, Maurice A, Josse JM, De-Vienne D (1994) Quantitative trait loci underlying gene product variation: a novel perspective for analyzing regulation of genome expression. Genetics 137: 289-301

Frewen Barbara E, Chen Tony HH, Howe Glenn T, Davis Joel, Rohde Antje, Boerjan Wout, Bradshaw HD (2000) Quantitative trait loci and candidate gene mapping of bud set and bud flush in Populus. Genetics 154: 837-845

Grattapaglia, D. and Sederoff, R., 1994. Genetic linkage maps of Eucalyptus grandis and Eucalyptus urophylla using a pseudo-test-cross: mapping strategy and RAPD markers. Genetics 137: 1121-1137

Grattapaglia D, Bertolucci FLG, Penchel R, Sederoff RR (1996) 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 144: 1205-1214

Groover, A, Devey M, Lee J, Megraw R, Mitchell-olds T, Sherman B, Vujcic S, Williams C, Neale D (1994) Identification of quantitative trait loci influencing wood specific gravity in an outbred pedigree of loblolly pine. Genetics 138: 1293-1300

Lander ES and Botstein D (1989) Mapping Mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics 121: 185-199

Marques CM, Vasquez-Kool J, Carocha, VJ, Ferreira, JG, O’Malley, DM, Liu, BH, Sederoff, R (1999) Genetic dissection of vegetative propagation traits in Eucalyptus tereticornis and E. globules. Theor Appl Genet 99: 936-946

Sewell MM, Bassoni DL, Megraw RA, Wheeler NC, Neale DB (2000) Identification of QTLs influencing wood property traits in loblolly pine (pinus taeda L.)£®I£®Physical wood properties. Theor Appl Genet 101: 1273-1281

Sewell MM, Davis MF, Tuskan GA, Wheeler NC, Elam CC, Bassoni DL, Neale DB (2002) Identification of QTLs influencing wood property traits in loblolly pine (pinus taeda L.). II. Chemical wood properties. Theor Appl Genet 104: 214-222

Veldboom LR, Lee M, Woodman WL (1994) Molecular marker-facilitated studies of morphological traits in an elite maize population. II. Determination of QTLs for grain yield and yield components. Theor Appl Genet 89: 451-458

Vos P, Hogers R, Bleeker M, Reijans M, van de Lee T, Hornes M, Frijters A, Pot J, Peleman J, Kuiper M, Zabeau M (1995) AFLP: a new technique for DNA fingerprinting. Nucleic Acids Res 23: 4407-4414

Wu R, Bradshaw HD, Stettler RF (1997) Molecular genetics of growth and development in Populus (Salicaceae). Mapping quantitative trait loci affecting leaf variation. American Journal of Botany 84: 143-153

Xiao J, Li J, Tanksley SD (1996) Identification of QTLs affecting traits of agronomic importance in a recombinant inbred population derived from a subspecific rice cross. Theor Appl Genet 92: 230-244


[1] Corresponding author: Institute of Populus tomentosa, Beijing Forestry University, Beijing 100083. Email: [email protected]