Yanyan Pan?Shuchun Li?Chenglu Wang?Wenjun Ma?Guiyou Xu?Lili Shao?Kailong Li?Xiyang Zhao?Tingbo Jiang

Introduction

Larix kaempferi(Larix Mill)has become one of the most important timber and afforestation species in China since it was introduced from Japan 100 years ago(Toda and Mikami 1976).Because of its strong adaptability,wide distribution,and rapid growth rate(Kurinobu 2005),L.kaempferi is used for many purposes,such as building,furniture making,paper making,and biofuel(Burczyk et al.1997;Fukatsu et al.2013).Previous studies on L.kaempferi have focused on the selection of provenances(Teruyoshi et al.2014;Nagasaka et al.2011),families(Hou et al.2012),and clones(Liu et al.2009),as well as its wood properties(Pa?ques et al.2010;Apiolaza 2009),photosynthetic characters(Fukatsu et al.2008;Crane et al.1983),genetic diversity(Achere 2004;Arcade et al.2000),and molecular markers(San et al.2009;Moriguchi et al.2008).

Among the several studies that focused on genetic variation in growth traits of L.kaempferi,most were based on trait measurements in a single year(Xu et al.2014;Zhang et al.2013;Takata et al.1992).However,for selecting superior materials,genetic variations in the chosen variables need to be monitored over several growth years(Zhao et al.2013;Xia et al.2016).Early selection can be conducted if growth parameters are correlated among different growth years(Weng et al.2007).In previous studies,Jiao et al.(2005)determined that the early selection time for L.kaempferi clones should be the 6th growth year,in contrast to the 10th year reported by Ma(2006).Different materials and growth environments may lead to different phenotypes,resulting in a range of estimates for the best early selection time.In the present study,for determining the most appropriate time for early selection of L.kaempferi seedlings,we measured tree height and diameter at breast height of 57 L.kaempferi clones over many years.Different clones were selected as excellent materials in different years under various selection ratios,and a suitable early selected time was proposed.This research provides a theoretical basis for breeding improvement of L.kaempferi.

Materials and methods

Materials and experimental design

The experimental field was located at Wudaogou Forestry Farm,Liuhe County(N 41°54′,E 125°17′),Jilin Province,China.This area has a continental monsoon climate.The annual average temperature,precipitation,frost-free period,and duration of sunshine of this area are 5.5°C,736 mm,145 days,and 2479 h,respectively.

The 57 clones of L.kaempferi(L1–L26,L60–L90)used for this study were propagated by grafting,and experimental plantations were established with 1-year-old plants in the spring of 1988 in a completely randomized block design(Marron and Ceulemans 2006)with 4 replication blocks and 10-tree row plot and 3.0×4.0 m spacing.

Data collection

The tree height(H)of all living and unbroken trees was measured from 1988(2nd growth year)to 1994(8th growth year),and in 1996(10th growth year),2003(17th growth year)and 2016(30th growth year),respectively.The diameter at breast height(DBH)of each tree was measured in 2003(17th growth year)and in 2016(30th growth year),respectively.

Statistical analyses

All data were analyzed using SPSS software,version 19.0(IBM,Armonk,NY,USA)(Shi 2012).The significance offixed effects was tested by analysis of variance(One-way ANOVA)F tests.Variations among families in different years were determined by ANOVA according to Hansen and Roulund(1997),as follows:

where yijis the performance of an individual of clone i in year j,μ is the overall mean,αiis the clone effect,βjis the year effect,αβijis the random effect of clone i in year j,and εijis the random error.

Variations among clones in the same year were determined by ANOVA according to Xu(2006)using the following equation:

where yijis the performance of individual clone i,μ is the overall mean,αiis the clone effect(i=1,…,50),and εijis the random error.

The coefficient of phenotypic variation(PCV)was calculated using the formula(Hai et al.2008):

where SD is the standard deviation andˉX is the phenotypic mean of the trait,respectively.

Repeatability(R)was calculated according to Hansen and Roulund(1997)as:

whereis the additive genetic variance component among clones,is the block variance,andis the error variance component.

The phenotypic correlation,rA(xy),of traits x and y was calculated according to Bi et al.(2000),as:

where covp（x，y）is the covariance between trait x and y,σp（x）is the variance component for trait x,and σp（y）is the variance component for trait y.

Comprehensive evaluations of different clones were calculated as described by Zhao et al.(2016),as:

where ai=Xij/Xjmax,Qiis a comprehensive evaluation value for each clone,Xijis the average value of one trait,Xjmaxis the maximum value of the trait;and n is the number of traits.

Genetic gain was estimated using the formula of Francisco et al.(2008):where W is the selection difference,andˉX and h2are the mean value and heritability of a given trait,respectively.

Results

Analyses of variance for H and DBH data

The H and DBH of 57 L.kaempferi clones were compared among different growth years by ANOVA(Table 1).All of the variance sources were significant(P＜0.01)except for the age×clone interaction(P=0.741)for DBH.

The ANOVA results for H and DBH among clones in the same growth year are shown in Table 2.For both H and DBH,there were highly significant differences among different clones in the same growth year(P＜0.01).

Genetic and phenotypic variations in different traits

The genetic and phenotypic variations in H and DBH among different clones are shown in Table 3.The average H of the 57 clones ranged from 0.53 to 19.09 m(from the 2nd to the 30th growth year).The average DBH of all clones was 17.56 cm in the 17th growth year and 29.03 cm in the 30th growth year.The PCVs of H and DBH in different growth years ranged from 14.89 to 35.65%and from 19.17 to 23.86%,respectively.The PCVs of H decreased as the trees grew.The PCVs of DBH were higher than those of H in the 17th and the 30th growth year.The R of H increased as the trees grew,and the R values were all higher than 0.6000.In the 17th and the 30th growth year,the R of DBH was 0.8497 and 0.8443,respectively,indicating high repeatability.

Correlation analysis

Correlations among H and DBH in different growth years are shown in Table 4.There were significant positive correlations between all pairs of growth traits except for H(in the 2nd growth year)×H(in the 30th growth year);and H(in the 2nd growth year)×DBH(in the 30th growth year).The correlation coefficients between H at different growth years ranged from 0.018(between 2nd and 30th growth year)to 0.974(between 7th and 8th growth year)and the values decreased as the age interval increased.The correlation coefficients between H in consecutive years increased from the 2nd growth year(0.749)to the 8th growth year(0.974).The correlation coefficient between H and DBH in the same growth year was extremely signif icant(0.483 in the 17th growth year;0.644 in the 30th growth year).

Comprehensive evaluation

In this study,excellent clones were selected based on H(in the 2nd,4th,7th,and 10th growth year)and Qivalues(in the 17th and the 30th growth year)with different selection ratios in different growth years.The results are shown in Tables 5 and 6.Thirty-five excellent clones were selected in the 2nd growth year under a 60%selection ratio.The average H of the selected clones was 0.56 m,which was 6.59%higher than the overall average value,and the genetic gain was 4.08%.Twenty-nine excellent clones were selected in the 4th growth year under a 50%selection ratio.The average H of selected clones was 2.06 m,which was 7.20%higher than the overall average value,and the genetic gain was 5.45%.In the 7th growth year,23 excellent clones were selected under a 40%selection ratio.The average H of selected clones was 4.29 m,which was 7.56%higher than overall average value,and the genetic gain was 5.54%.In the 10th growth year,18 excellent clones were selected under a 30%selection ratio.The average H of selected clones was 5.92 m,which was 7.74%higher than the overall average value,and the genetic gain was 5.69%.

Table 1 ANOVA of H and DBH of 57 L.kaempferi clones

Table 2 ANOVA of different traits in the same year of 57 L.kaempferi clones

Table 3 Variation parameters of different clones in different years

Table 4 Correlation analysis of different traits in different growth years in Larix clones

In the 17th growth year,12 excellent clones were selected under a 20% selection ratio with the comprehensive evaluation method.The average H and DBH of selected clones was 10.19 m and 20.67 cm,which were 12.78 and 17.71%higher than their respective overall average values,and the genetic gain was 10.91 and15.04%,respectively.In the 30th growth year,six excellent clones were selected under a 10%selection ratio with the comprehensive evaluation method.The H and DBH of selected clones was 20.96 m and 33.84 cm,respectively,9.81 and 16.57%higher than their overall average values.The genetic gain in H and DBH in the 30th growth year was 6.46 and 13.99%,respectively.

Table 5 Excellent clones selected under different selected rate in different growth years

Discussion

Analysis of variance

Among the various statistical analyses,ANOVA is the most important method to estimate the extent and signif icance of variability(Zhao et al.2014).In this study,the ANOVA results showed that all sources of variance were highly significant(P＜0.01),except for the age×clone interaction for DBH(P=0.741).These results were consistent with those reported for Pinus taeda(Isik et al.2003)and L.kaempferi(Nakada et al.2005),indicating that the selection of excellent clones is practical and feasible.

Table 6 Average and variation parameters of selected clones in different year

Genetic variation analysis

Genetic and phenotypic variations are the main factors in tree genetic and breeding(Mwase et al.2008),and both are very important for the selection and evaluation of excellent clones(Safavi et al.2010).In the process of tree growth,genetic and phenotypic variations in growth traits undergo dynamic changes and show certain patterns(Karacic and Weih 2006).Understanding the dynamic changes in genetic variables is important for early selection(Pliura et al.2007)and the genetic improvement of forest trees(Sixto et al.2011;Neale 2007).In this research,the PCVs of H ranged from 35.65 to 14.89%,and the PCVs decreased as the trees grew.A similar phenomenon was reported in other studies(Shigematsu 1991;Oshima 1998);that is,the differences in H among clones decreased as the trees grew.In the 30th growth year,the PCV value was 14.89%,suggesting that tree height was close to the critical value.

The size of the R value indicates the reliability with which the genotype will be recognized by its phenotypic expression(Zhao et al.2013).The R values of H and DBH ranged from 0.6181 to 0.8531 and from 0.8443 to 0.8497,respectively,in different growth years.These high values indicated that the variation in H and DBH was strongly controlled by genetic factors,consistent with the results of a previous study(Du et al.2015).For H,the R values increased as the trees grew.Similar findings were also reported for Picea abies(Isik et al.2010)and Pinus sylvestris(Haapanen 2001),and indicated that the effect of the environment decreased as the trees grew.Overall,high R values indicate that the selection for excellent traits would be effective(Maniee et al.2009).The high variation and R values of growth traits in this study were useful for the selection and evaluation of excellent clones.

Correlation analysis

In previous studies,age–age correlations of growth traits were important parameters for early selection(Goncalves et al.2005).Age–age correlations for growth have been reported for Pinus radiata(Matheson et al.1994;King and Burdon 1991),but the relationship between H and DBH was rather complex(Sumida et al.2013).In this study,there were significant positive correlations between all pairs of growth traits except for H(in the 2nd growth year)×H(in the 30th growth year);and H(in the 2nd growth year)×DBH(in the 30th growth year).These results indicated that there was a close relationship between early and late growth traits.

The results indicated that early selection of trees was feasible (Pa?ques 2004).The correlation coefficients between DBH(in the 17th and 30th growth year)and H(in the 2nd to 10th growth year)were higher than those between H(in the 17th and 30th growth year)and H(in the 2nd to 10th growth year).Therefore,although H would be a good selection variable for clones with fast growth at an early stage,DBH would be a better to select clones with strong growth at a later stage,consistent with the findings of another study(Zas et al.2004).The high correlation coefficients also suggested that excellent clones could be selected early based on their measured traits.

Selection of excellent clones

Early selection is an important method to shorten the breeding cycle and accelerate the selection process(Kroon et al.2011;Kumar and Singh 2001).The optimum age for selection depends on the correlation coefficients between various variables and on genetic stability(Gwaze and Bridgwater 2002).The optimum age for early selection of L.kaempferi differed depending on the selected breeding targets(Fujimoto et al.2006;Diao et al.2016).In terms of growth traits,some studies found that H was the best factor for early selection(Sun et al.2004;Lai et al.2014).The optimum age for early selection was reported to be the 6th to the 8th growth year(Sun 2003;Balocchi et al.1993).Other studies suggested that DBH was a more appropriate than H for early selection(Vasquez and Dvorak 1996)because it is harder to measure H accurately as trees grow and the DBH measurements were more accurate(Zhou et al.2008).In this research,H was selected as an early evaluation index,and different clones were selected in different growth years under different selection ratios.Under a selection ratio of 10%,the same excellent clones were selected in the 7th and the 30th growth year,indicating that the early selection of these materials could be conducted from the 6th to the 8th growth year.

Some studies reported that selection based on dual growth traits(H and DBH)was more effective than that based on a single trait(Collet and Chenost 2006;Chen et al.2003),which has important implications for breeding programs.In this study,the combination of H and DBH was used to select excellent clones in the 30th growth year,and L65,L1,L62,L9,L15,and L78 were selected as the excellent clones.Among these clones,L65 and L62 showed fast growth in the later growth stage,but clones L1,L9,L15,and L78 grew fast all the time.Thus,the clones had different growth patterns,consistent with the findings of another study(Jansson et al.2003).In the process of selecting clones,the use of appropriate selection ratios and the choice of appropriate selection years are vital to select the best clones(Leksono et al.2006).

Genetic gain

Genetic gain was one of the most important variables for selecting excellent materials.The level of genetic gain depends on the repeatability,amplitude of genetic variation,and selection rate(Zhou et al.2014).Most breeding programs aim to increase genetic gain via selection(Montes et al.2008).In our research,the genetic gains of H increased from the 2nd to the 17th growth year(4.08 to 10.91%),but decreased from the 17th to the 30th growth year(from 10.91 to 6.46%),possibly because of smaller annual increases in height that led to smaller differences among clones,as reflected by the low PCV values.

Conclusions

Early selection is very important to shorten the breeding cycle and improve breeding efficiency.In this research,there were highly significant differences in all the traits among 57 L.kaempferi clones,and the high PCV and R values indicated that it was feasible to select excellent clones.There were significant positive age–age correlations in growth traits among different growth years.Under the same selection ratio,the clones selected in the 7th growth year were the same as those selected in the 30th growth year,suggesting that the optimum early selection age is the 6th–8th growth year.The selected clones could serve as materials for the genetic improvement of L.kaempferi.This study provides a theoretical foundation for the early selection of L.kaempferi.

Author’s contributionsYP,SL and CW have contributed equally to this work.

Achere V(2004)Chloroplast and mitochondrial molecular tests identify European×Japanese larch hybrids.Theor Appl Genet 108:1643–1649

Apiolaza LA(2009)Very early selection for solid wood quality:screening for early winners.Ann For Sci.doi:10.1051/forest/2009047

Arcade A,Anselin F,Faivre RP(2000)Application of AFLP,RAPD and ISSR markers to genetic mapping of European and Japanese larch.Theor Appl Genet 100:299–307

Balocchi C,Bridgwater F,Zobel B,Jahromi S(1993)Age trends in genetic parameters for tree height in a nonselected population of loblolly pine.For Sci 39(2):231–251

Bi CX,Guo JZ,Wang HW,Shu QY(2000)The correlation and path analysis on the quantity characters of Chinese pine.J Northwest For Univ 15(2):7–12

Burczyk J,Nikkanen T,Lewandowski A(1997)Evidence of an unbalanced mating pattern in a seed orchard composed of two larch species.Silvae Genet 46:176–181

Chen XY,Hawkins B,Xie CY(2003)Age trends in genetic parameters and early selection of lodgepole pine provenances with particular reference to the lambeth model.For Genet 10(3):249–258

Collet C,Chenost C(2006)Using competition and light estimates to predict diameter and height growth of naturally regenerated beech seedlings growing under changing canopy conditions.Forestry 79(5):489–502

Crane J,Dickmann D,Flore J(1983)Photosynthesis and transpiration by young Larix kaempferi trees:C3 responses to light and temperature.Physiol Plant 59(4):635–640

Diao S,Hou YM,Xie YH,Sun XM(2016)Age trends of genetic parameters,early selection and family by site interactions for growth traits in Larix kaempferi open-pollinated families.BMC Genet 17:104–115

Du CQ,Xu YZ,Sun XH,Liu YH,Xie PS(2015)Variation of growth traits and early selection of Larix Kaempferi clones in sub-alpine area of western Hubei Province.J Huazhong Agric Univ 34(3):19–23

Francisco F,Messias G,Helaine C,Pedro C,Telma N,Ana P,Alexandre P,Geraldo A(2008)Selection and estimation of the genetic gain in segregating generations of papaya(Carica papaya L.).Crop Breed Appl Biotechnol 8:1–8

Fujimoto T,Kita K,Uchiyama K(2006)Age trends in the genetic parameters of wood density and the relationship with growth rates in hybrid larch(Larix gmelinii var.japonica×L.kaemferi)F1.Jpn For Soc Springer Verlag Tokyo 11:157–163

Fukatsu E,Fukuda Y,Takahashi M,Nakada R(2008)Clonal variation of carbon content in wood of Larix kaempferi(Japanese larch).J Wood Sci 54:247–251

Fukatsu E,Tsubomura M,Fujisawa Y,Nakada R(2013)Genetic improvement of wood density and radial growth in Larix kaempferi:results from a diallel mating test.Ann For Sci 70:451–459

Goncalves P,Bortoletto N,Cardinal A,Gouvea L,Costa R,Moraes M(2005)Age–age correlation for early selection of rubber tree genotypes in Sao Paulo State,Brazil.Genet Mol Biol 28:758–764

Gwaze DP,Bridgwater FE(2002)Determining optimum selection age for diameter and height in loblolly pine(Pinus taeda).For Genet 9:159–165

Haapanen M(2001)Time trends in genetic parameter estimates and selection efficiency for Scots pine in relation to field testing method.For Genet 8:129–144

Hai PH,Jansson G,Harwood C,Hannrup B,Thinh HH(2008)Genetic variation in growth,stem straightness and branch thickness in clonal trials of Acacia auriculiformis at three contrasting sites in Vietnam.For Ecol Manag 255:156–167

Hansen J,Roulund H(1997)Genetic parameters for spiral grain,stem form,pilodyn and growth in 13 years old clones of Sitka Spruce(Picea sitchensis(Bong.)Carr.).Silvae Genet 46:107–113

Hou YM,Quan YS,Cao J,Guan LH(2012)Study on superior family selection of Larix kaempferi.J Cent South Univ For Technol 3(9):6–10

Isik F,Li BL,Frampton J(2003)Estimates of additive,dominance and epistatic genetic variances from a clonally replicated test of loblolly pine.For Sci 49(1):77–88

Isik K,Kleinschmit J,Steiner W(2010)Age–age correlations and early selection for height in a clonal genetic test of Norway spruce.For Sci 56(2):212–221

Jansson G,Li BL,Hannrup B(2003)Time trends in genetic parameters for height and optimum age for parental selection in Scots pine.For Sci 49:696–705

Jiao YD,Zhao K,Qin YF,Hu ZY,Liu HF(2005)Study on superior clones and early selection of Larix kaempferi.J Henan For Sci Technol 25(2):9–11

Karacic A,Weih M(2006)Variation in growth and resource utilization among eitht poplar clones grown under different irrigation and fertilization regimes in Sweden.Biomass Bioenerg 30:115–124

King JN,Burdon RD(1991)Time trends in inheritance and projected efficiencies of early selection in a large 17-year-old progeny test of Pinus radiata.Can J For Res 21(8):1200–1207

Kroon J,Ericsson T,Jansson G,Andersson B(2011)Patterns of genetic parameters for height in field genetic tests of Picea abies and Pinus sylvestris in Sweden.Tree Genet Genomes 7(6):1099–1111

Kumar D,Singh NB(2001)Age–age correlation for early selection of clones of Populus in India.Silvae Genet 50:103–108

Kurinobu S(2005)Forest tree breeding for Japanese larch.Eurasian J For Res 8:127–134

Lai M,Sun XM,Chen DS,Xie YH,Zhang SG(2014)Age-related trends in genetic parameters for Larix kaempferi and their implications for early selection.BMC Genet 15(Suppl 1):S10

Leksono B,Kurinobu S,Ide Y(2006)Optimum age for selection base on a time trend of genetic parameters related to diameter growth in seedling seed orchards of Eucalyptus pellita in Indonesia.J For Res 11:359–363

Liu LX,Bai YR,Wu ZY,Zhao A,Ma ZG(2009)Experiment on choosing excellent clones of Larix kaempferi Carr.J Inner Mong For Technol 35(2):10–13

Ma SX(2006)Studies on genetic variation and early selection of Japanese larch clone.Master Thesis of Henan Agricultural University

Maniee MD,Kahrizi MD,Mohammadi R(2009)Genetic variability of some morphophysiological traits in durum wheat(Triticum turgidum var.durum).J Appl Sci 9:1383–1387

Marron N,Ceulemans R(2006)Genetic variation of leaf traits related to productivity in a Populus deltoides×Populus nigra family.Can J For Res 36:390–400

Matheson A,Spencer DJ,Magnussen D(1994)Optimum age for selection in pinus radiata using basal area under bark for age–age correlations.Silvae Genet 43:352–357

Montes CS,Hernandez RE,Beaulieu J,Weber JC(2008)Genetic variation in wood color and its correlations with tree growth and wood density of Calycophyllum spruceanum at an early age in the Peruvian Amazon.New For 35:57–73

Moriguchi Y,Kita K,Uchiyama K,Kuromaru M,Tsumura Y(2008)Enhanced hybridization rates in a Larix gmelinii var.japonica×L.kaempferi interspecific seed orchard with a single maternal clone revealed by cytoplasmic DNA markers.Tree Genet Genomes 4:637–645

Mwase WF,Savill PS,Hemery G(2008)Genetic parameter estimates for growth and form traits in common ash(Fraxinus excelsior,L.)in a breeding seedling orchard at Little Wittenham in England.New For 36:225–238

Nagasaka K,Yoshimura K,Akashi T,Arai K,Yamamoto C(2011)Assessment of characteristics and grouping of Larix kaempferi provenances in provenance test sites in Nagano prefecture in Japan.J Jap For Soc 93:179–186

Nakada R,Fujisawa Y,Taniguchi T(2005)Variations of wood properties between plus-tree clones in Larix kaempferi(Lamb.)Carrie`re.Bull For Tree Breed Cent 21:85–105

Neale DB(2007)Genomics to tree breeding and forest health.Curr Opin Genet Dev 17:539–544

Oshima T(1998)The age trends in the combining ability and heritability in the wood density in Japanese larch.Trans Meet Hokkaido Branch Jpn For Soc 46:157–159

Pa?ques LE(2004)Roles of European and Japanese larch in the genetic control of growth,architecture and wood quality traits in interspecific hybrids(Larix×eurolepis Henry).Ann For Sci 61:25–33

Pa?ques LE,Millier F,Rozenberg P(2010)Selection perspectives for genetic improvement of wood stiffness in hybrid larch(Larix × eurolepis Henry).Tree Genet Genomes 6:83–92

Pliura A,Zhang SY,Mackay J,Bousquet J(2007)Genotypic variation in wood density and growth traits of poplar hybrids at four clonal trials.For Ecol Manag 238:92–106

Safavi SA,Pourdad SA,Mohmmad T,Mahmoud K(2010)Assessment of genetic variation among safflower(Carthamus tinctorius L.)accessions using agro-morphological traits and molecular markers.J Food Agric Environ 8(1):616–625

San J,Uchida K,Tomaru N(2009)Mitochondrial DNA variation in natural populations of Japanese larch(Larix kaempferi).Silvae Genet 58:5–6

Shi LW(2012)Statistical analysis from entry to the master of SPSS 19.0.Tsinghua University Press,Beijing

Shigematsu Y(1991)The wood quality of planting trees for Japanese larch relating with growth.II.The grade of wood for planting trees of Japanese larch.Mokuzai Kogyo 46:9–16

Sixto H,Salvia J,Barrio M,Ciria MP,Canellas I(2011)Genetic variation and genotype-environment interactions in short rotation Populus plantations in southern Europe.New For 42:163–177

Sumida A,Miyaura T,Torii H(2013)Relationships of tree height and diameter at breast height revisited:analyses of stem growth using 20-year data of an even-aged Chamaecyparis obtusa stand.Tree Physiol 33:106–118

Sun XM(2003)Selection of superior Larix kaempferi families for pulpwood purpose and family based growth simulation.Ph.D.Thesis of Chinese Academy of Forestry

Sun XM,Zhang SG,Hou YM,Li SY(2004)Age trends of genetic parameters for growth traits in short rotation Larix kaempferi families.Sci Silvae Sin 40(6):68–74

Takata K,Koizumi A,Ueda K(1992)Variations of radial growth and wood quality among provenances in Japanese larch.Mokuzai Gakkaishi 38:1082–1088

Teruyoshi N,Kazutoshi N,Hiroshi Y,Yoshihiko T(2014)Provenance tests for survival and growth of 50-year-old Japanese larch(Larix kaempfer i)trees related to climatic conditions in central Japan.Tree Genet Genomes 10:87–99

Toda R,Mikami S(1976)The provenance trials of Japanese larch established in Japan and the tentative achievements.Silvae Genet 25:209–216

Vasquez J,Dvorak WS(1996)Trend in variance and heritabilities with stand development of tropical pines.Can J For Res 26:1473–1480

Weng YH,Tosh KJ,Park YS,Fullarton MS(2007)Age-related trends in genetic parameters for Jack Pine and their implications for early selection.Silvae Genet 56(5):242–251

Xia H,Zhao GH,Zhang LS,Sun XY,Yin SP,Liang DY,Li Y,Zheng M,Zhao XY(2016)Genetic and variation analyses of growth traits of half-sib Larix olgensis families in northeastern China.Euphytica.doi:10.1007/s10681-016-1765-4

Xu JR (2006)Quantitative genetics of forest.China Forestry Publishing,Beijing

Xu SD,Sun XM,Chen JL,Li DX,Wang GF,Huang XF,Qi WY(2014)Primary research on fine family and clone of Larix kaempferi.Hubei For Sci Technol 43(3):12–16

Zas R,Merlo E,Ferna′ndez-Lo′pez J(2004)Genetic parameter estimates for Maritime pine in the Atlantic coast of North-west Spain.For Genet 11(1):45–53

Zhang ZG,Ma JW,Jin XC,Wang HC,Ding PJ,Hu MH,Pan CL,Lei SG(2013)Analysis and selection of progeny test forest of the open pollinated family of Japanese larch.J Northwest For Univ 28(4):74–79

Zhao XY,Hou W,Zheng HQ,Zhang ZY(2013)Analyses of genotypic variation in white poplar at four sites in China.Silvae Genet 62(4–5):187–195

Zhao XY,Bian XY,Liu MR,Li ZX,Li Y,Zheng M,Teng WH,Jiang J,Liu GF(2014)Analysis of genetic effects on complete diallel cross test of Betula platyphylla.Euphytica 200:221–229

Zhao XY,Xia H,Wang XW,Wang C,Liang DY,Li KL,Liu GF(2016)Variance and stability analyses of growth characters in halfsib Betula platyphylla families at three different sites in China.Euphtyica 208:173–186

Zhou GX,Hu X,Zhu Q,Ning SH(2008)Research on the relationship between stand productivity and site condition of Larix kaempferi.Ecol Econ 4:216–226

Zhou H,Fries A,Wu HX(2014)High negative genetic correlations between growth traits and wood properties suggest incorporating multiple traits selection including economic weights for the future Scots pine breeding programs.Ann For Sci 71:463–472