Molecular genetics of growth and development in populus. IV. Mapping QTLs with large effects on growth, form, and phenology traits in a forest tree

Climate-driven convergent evolution in riparian ecosystems on sky islands

Climate-induced evolution will determine population persistence in a changing world. However, finding natural systems in which to study these responses has been a barrier to estimating the impact of global change on a broad scale. We propose that isolated sky islands (SI) and adjacent mountain chains (MC) are natural laboratories for studying long-term and contemporary climatic pressures on natural populations. We used greenhouse common garden trees to test whether populations on SI exposed to hot and dry climates since the end of the Pleistocene have phenotypically diverged from populations on MC, and if SI populations have converged in these traits. We show: (1) populations of Populus angustifolia from SI have diverged from MC, and converged across SI, in reproductive and productivity traits, (2) these traits (cloning and aboveground biomass, respectively) are significantly correlated, suggesting a genetic linkage between them, and (3) the trait variation is driven by both natural selection and genetic drift. These shifts represent potentially beneficial phenotypes for population persistence in a changing world. These results suggest that the SI-MC comparison is a natural laboratory, as well as a predictive framework, for studying long-term responses to climate change across the globe.

QTL mapping of drought-related traits in the hybrids of Populus deltoides 'Danhong'×Populus simonii 'Tongliao1'

BACKGROUND

Poplar trees provide a large amount of wood material, but many parts of the world are arid or semi-arid areas because of insufficient annual precipitation, which seriously affects the growth of poplar trees. Populus simonii 'Tongliao1' shows strong tolerance to stress environments, and Populus deltoides 'Danhong' shows a stronger growth rate in a suitable environment. To identify drought tolerance-related QTLs and genes, an F₁ population derived from the cross between the 'Danhong' and 'Tongliao 1' Populus was assessed under drought stress.

RESULTS

We measured drought-related traits such as the relative height growth, relative diameter growth, leaf senescence number, specific leaf area, and leaf relative water content in the population under control and drought environments. The results showed that drought stress reduced the plant height relative growth, ground diameter relative growth, specific leaf area and leaf relative water content and increased the number of leaf drops. A total of 208 QTLs were identified by QTL mapping analysis, and they consisted of 92, 63 and 53 QTLs under control, drought stress treatment and drought index conditions, respectively. A molecular identification marker for drought tolerance, np2841, which was associated with a QTL (qDLRWC-LG10-1) for relative leaf water content, was initially developed. We mined 187 candidate genes for QTL regions of five traits under a drought environment. The reference genome annotation for Populus trichocarpa and a homologous gene analysis of Arabidopsis thaliana identified two candidate genes, Potri.003G171300 and Potri.012G123900, with significant functions in response to drought stress. We identified five key regulatory genes (Potri.006G273500, Potri.007G111500, Potri.007G111600, Potri.007G111700, and Potri.007G111800) related to drought tolerance through the poplar coexpression network.

CONCLUSION

In this study, our results indicate that the QTLs can effectively enhance the drought tolerance of poplar. It is a step closer towards unravelling the genetic basis of poplar drought tolerance-related traits, and to providing validated candidate genes and molecular markers for future genetic improvement.

Time-sequential detection of quantitative trait loci and candidate genes underlying the dynamic growth of Salix suchowensis

Elucidating the genetic factors underlying long-term biological processes remains challenging since the relevant genes and their effects may vary across different developmental stages. In this study, we carried out a large-scale field trial of the progeny of an F1 full-sib pedigree of Salix suchowensis and measured plant height and ground diameter periodically over a time course of 240 days. With the obtained data, we characterized plant growth rhythms and performed time-sequential analyses of quantitative trait loci underlying the dynamic growth of the plants. The dynamic mapping of quantitative trait loci revealed that stem height and ground diameter were under the control of four quantitative trait loci, and the effects of these quantitative trait loci varied greatly throughout the growth process, in which two quantitative trait loci were found to exert a pleiotropic effect determining the correlation between stem height and ground diameter. The analysis of candidate genes in the target genetic intervals showed that the pleiotropic effect of the two quantitative trait loci arises from the colocalization of genes with independent effects on stem height and ground diameter. Further examination of the expression patterns of the candidate genes indicated that height and circumference growth involve different activities of leaf and cambium tissues. This study provides unprecedented information to help us understand the dynamic growth of plants and presents an applicable strategy for elucidating the genetic mechanism underlying a long-term biological process by using plant growth as an example.

Simulating selection and evolution at the community level using common garden data

A key issue in evolutionary biology is whether selection acting at levels higher than the individual can cause evolutionary change. If it can, then conceptual and empirical studies must consider how selection operates at multiple levels of biological organization. Here, we test the hypothesis that estimates of broad‐sense community heritability, HC2, can be used to predict the evolutionary response by community‐level phenotypes when community‐level selection is imposed. Using an approach informed by classic quantitative genetics, we made three predictions. First, when we imposed community‐level selection, we expected a significant change in the average phenotype of arthropod communities associated with individual tree genotypes [we imposed selection by favoring high and low NMDS (nonmetric multidimensional scaling) scores that reflected differences in arthropod species richness, abundance and composition]. Second, we expected HC2 to predict the magnitude of the community‐level response. Third, we expected no significant change in average NMDS scores with community‐level selection imposed at random. We tested these hypotheses using three years of common garden data for 102 species comprising the arthropod communities, associated with nine clonally replicated Populus angustifolia genotypes. Each of our predictions were met. We conclude that estimates of HC2 account for the resemblance among communities sharing common ancestry, the persistence of community composition over time, and the outcome of selection when it occurs at the community level. Our results provide a means for exploring how this process leads to large‐scale community evolutionary change, and they identify the circumstances in which selection may routinely act at the community level.

Spatial and temporal shifts in photoperiod with climate change

Climate change causes both temporal (e.g. advancing spring phenology) and geographic (e.g. range expansion poleward) species shifts, which affect the photoperiod experienced at critical developmental stages ('experienced photoperiod'). As photoperiod is a common trigger of seasonal biological responses - affecting woody plant spring phenology in 87% of reviewed studies that manipulated photoperiod - shifts in experienced photoperiod may have important implications for future plant distributions and fitness. However, photoperiod has not been a focus of climate change forecasting to date, especially for early-season ('spring') events, often assumed to be driven by temperature. Synthesizing published studies, we find that impacts on experienced photoperiod from temporal shifts could be orders of magnitude larger than from spatial shifts (1.6 h of change for expected temporal vs 1 min for latitudinal shifts). Incorporating these effects into forecasts is possible by leveraging existing experimental data; we show that results from growth chamber experiments on woody plants often have data relevant for climate change impacts, and suggest that shifts in experienced photoperiod may increasingly constrain responses to additional warming. Further, combining modeling approaches and empirical work on when, where and how much photoperiod affects phenology could rapidly advance our understanding and predictions of future spatio-temporal shifts from climate change.

Multivariate linear mixed model enhanced the power of identifying genome-wide association to poplar tree heights in a randomized complete block design

With the advances in high-throughput sequencing technologies, it is not difficult to extract tens of thousands of single-nucleotide polymorphisms (SNPs) across many individuals in a fast and cheap way, making it possible to perform genome-wide association studies (GWAS) of quantitative traits in outbred forest trees. It is very valuable to apply traditional breeding experiments in GWAS for identifying genome variants associated with ecologically and economically important traits in Populus. Here, we reported a GWAS of tree height measured at multiple time points from a randomized complete block design (RCBD), which was established with clones from an F₁ hybrid population of Populus deltoides and Populus simonii. A total of 22,670 SNPs across 172 clones in the RCBD were obtained with restriction site-associated DNA sequencing (RADseq) technology. The multivariate mixed linear model was applied by incorporating the pedigree relationship matrix of individuals to test the association of each SNP to the tree heights over 8 time points. Consequently, 41 SNPs were identified significantly associated with the tree height under the P-value threshold determined by Bonferroni correction at the significant level of 0.01. These SNPs were distributed on all but two chromosomes (Chr02 and Chr18) and explained the phenotypic variance ranged from 0.26% to 2.64%, amounting to 63.68% in total. Comparison with previous mapping studies for poplar height as well as the candidate genes of these detected SNPs were also investigated. We therefore showed that the application of multivariate linear mixed model to the longitudinal phenotypic data from the traditional breeding experimental design facilitated to identify far more genome-wide variants for tree height in poplar. The significant SNPs identified in this study would enhance understanding of molecular mechanism for growth traits and would accelerate marker-assisted breeding programs in Populus.

High-Quality SNP Linkage Maps Improved QTL Mapping and Genome Assembly in Populus

With the advances in high-throughput sequencing technologies and the development of new software for extracting single nucleotide polymorphisms (SNPs) across a mapping population, it is possible to construct high-quality genetic maps with thousands of SNPs in outbred forest trees. Two parent-specific linkage maps were constructed with restriction site-associated DNA sequencing data from an F₁ hybrid population derived from Populus deltoides and Populus simonii, and applied in QTL mapping and genome assembly. The female P. deltoides map contained 4018 SNPs, which were divided into 19 linkage groups under a wide range of LOD thresholds from 7 to 55. The male P. simonii map showed similar characteristics, consisting of 2097 SNPs, which also belonged to 19 linkage groups under LOD thresholds of 7 to 29. The SNP order of each linkage group was optimal among different ordering results from several available software. Moreover, the linkage maps allowed the detection of 39 QTLs underlying tree height and 47 for diameter at breast height. In addition, the linkage maps improved the anchoring of 689 contigs of P. simonii to chromosomes. The 2 parental genetic maps of Populus are of high quality, especially in terms of SNP data quality, the SNP order within linkage groups, and the perfect match between the number of linkage groups and the karyotype of Populus, as well as the excellent performances in QTL mapping and genome assembly. Both approaches for extracting and ordering SNPs could be applied to other species for constructing high-quality genetic maps.

Evolutionary network genomics of wood formation in a phylogenetic survey of angiosperm forest trees

Wood formation was present in early angiosperms, but has been highly modified through evolution to generate the anatomical diversity seen in extant angiosperm lineages. In this project, we modeled changes in gene coexpression relationships associated with the evolution of wood formation in a phylogenetic survey of 13 angiosperm tree species. Gravitropic stimulation was used as an experimental treatment to alter wood formation and also perturb gene expression. Gene transcript abundances were determined using RNA sequencing of developing wood tissues from upright trees, and from the top (tension wood) and bottom (opposite wood) tissues of gravistimulated trees. A network-based approach was employed to align gene coexpression networks across species based on orthologous relationships. A large-scale, multilayer network was modeled that identified both lineage-specific gene coexpression modules and modules conserved across multiple species. Functional annotation and analysis of modules identified specific regulatory processes associated with conserved modules, including regulation of hormones, protein phosphorylation, meristem development and epigenetic processes. Our results provide novel insights into the evolution and development of wood formation, and demonstrate the ability to identify biological processes and genes important for the evolution of a foundational trait in nonmodel, undomesticated forest trees.

Endophyte-Promoted Phosphorus Solubilization in Populus

Phosphorus is one of the essential nutrients for plant growth, but it may be relatively unavailable to plants because of its chemistry. In soil, the majority of phosphorus is present in the form of a phosphate, usually as metal complexes making it bound to minerals or organic matter. Therefore, inorganic phosphate solubilization is an important process of plant growth promotion by plant associated bacteria and fungi. Non-nodulating plant species have been shown to thrive in low-nutrient environments, in some instances by relying on plant associated microorganisms called endophytes. These microorganisms live within the plant and help supply nutrients for the plant. Despite their potential enormous environmental importance, there are a limited number of studies looking at the direct molecular impact of phosphate solubilizing endophytic bacteria on the host plant. In this work, we studied the impact of two endophyte strains of wild poplar (Populus trichocarpa) that solubilize phosphate. Using a combination of x-ray imaging, spectroscopy methods, and proteomics, we report direct evidence of endophyte-promoted phosphorus uptake in poplar. We found that the solubilized phosphate may react and become insoluble once inside plant tissue, suggesting that endophytes may aid in the re-release of phosphate. Using synchrotron x-ray fluorescence spectromicroscopy, we visualized the nutrient phosphorus inside poplar roots inoculated by the selected endophytes and found the phosphorus in both forms of organic and inorganic phosphates inside the root. Tomography-based root imaging revealed a markedly different root biomass and root architecture for poplar samples inoculated with the phosphate solubilizing bacteria strains. Proteomics characterization on poplar roots coupled with protein network analysis revealed novel proteins and metabolic pathways with possible involvement in endophyte enriched phosphorus uptake. These findings suggest an important role of endophytes for phosphorus acquisition and provide a deeper understanding of the critical symbiotic associations between poplar and the endophytic bacteria.

Quantitative genetic architecture of adaptive phenology traits in the deciduous tree, Populus trichocarpa (Torr. and Gray)

In a warming climate, the ability to accurately predict and track shifting environmental conditions will be fundamental for plant survival. Environmental cues define the transitions between growth and dormancy as plants synchronise development with favourable environmental conditions, however these cues are predicted to change under future climate projections which may have profound impacts on tree survival and growth. Here, we use a quantitative genetic approach to estimate the genetic basis of spring and autumn phenology in Populus trichocarpa to determine this species capacity for climate adaptation. We measured bud burst, leaf coloration, and leaf senescence traits across two years (2017–2018) and combine these observations with measures of lifetime growth to determine how genetic correlations between phenology and growth may facilitate or constrain adaptation. Timing of transitions differed between years, although we found strong cross year genetic correlations in all traits, suggesting that genotypes respond in consistent ways to seasonal cues. Spring and autumn phenology were correlated with lifetime growth, where genotypes that burst leaves early and shed them late had the highest lifetime growth. We also identified substantial heritable variation in the timing of all phenological transitions (h² = 0.5–0.8) and in lifetime growth (h² = 0.8). The combination of additive variation and favourable genetic correlations in phenology traits suggests that populations of cultivated varieties of P. Trichocarpa may have the capability to adapt their phenology to climatic changes without negative impacts on growth.

Marginal Lands to Grow Novel Bio-Based Crops: A Plant Breeding Perspective

The biomass demand to fuel a growing global bio-based economy is expected to tremendously increase over the next decades, and projections indicate that dedicated biomass crops will satisfy a large portion of it. The establishment of dedicated biomass crops raises huge concerns, as they can subtract land that is required for food production, undermining food security. In this context, perennial biomass crops suitable for cultivation on marginal lands (MALs) raise attraction, as they could supply biomass without competing for land with food supply. While these crops withstand marginal conditions well, their biomass yield and quality do not ensure acceptable economic returns to farmers and cost-effective biomass conversion into bio-based products, claiming genetic improvement. However, this is constrained by the lack of genetic resources for most of these crops. Here we first review the advantages of cultivating novel perennial biomass crops on MALs, highlighting management practices to enhance the environmental and economic sustainability of these agro-systems. Subsequently, we discuss the preeminent breeding targets to improve the yield and quality of the biomass obtainable from these crops, as well as the stability of biomass production under MALs conditions. These targets include crop architecture and phenology, efficiency in the use of resources, lignocellulose composition in relation to bio-based applications, and tolerance to abiotic stresses. For each target trait, we outline optimal ideotypes, discuss the available breeding resources in the context of (orphan) biomass crops, and provide meaningful examples of genetic improvement. Finally, we discuss the available tools to breed novel perennial biomass crops. These comprise conventional breeding methods (recurrent selection and hybridization), molecular techniques to dissect the genetics of complex traits, speed up selection, and perform transgenic modification (genetic mapping, QTL and GWAS analysis, marker-assisted selection, genomic selection, transformation protocols), and novel high-throughput phenotyping platforms. Furthermore, novel tools to transfer genetic knowledge from model to orphan crops (i.e., universal markers) are also conceptualized, with the belief that their development will enhance the efficiency of plant breeding in orphan biomass crops, enabling a sustainable use of MALs for biomass provision.

Exome resequencing and GWAS for growth, ecophysiology, and chemical and metabolomic composition of wood of Populus trichocarpa

Background

Populus trichocarpa is an important forest tree species for the generation of lignocellulosic ethanol. Understanding the genomic basis of biomass production and chemical composition of wood is fundamental in supporting genetic improvement programs. Considerable variation has been observed in this species for complex traits related to growth, phenology, ecophysiology and wood chemistry. Those traits are influenced by both polygenic control and environmental effects, and their genome architecture and regulation are only partially understood. Genome wide association studies (GWAS) represent an approach to advance that aim using thousands of single nucleotide polymorphisms (SNPs). Genotyping using exome capture methodologies represent an efficient approach to identify specific functional regions of genomes underlying phenotypic variation.

Results

We identified 813 K SNPs, which were utilized for genotyping 461 P. trichocarpa clones, representing 101 provenances collected from Oregon and Washington, and established in California. A GWAS performed on 20 traits, considering single SNP-marker tests identified a variable number of significant SNPs (p-value < 6.1479E-8) in association with diameter, height, leaf carbon and nitrogen contents, and δ¹⁵N. The number of significant SNPs ranged from 2 to 220 per trait. Additionally, multiple-marker analyses by sliding-windows tests detected between 6 and 192 significant windows for the analyzed traits. The significant SNPs resided within genes that encode proteins belonging to different functional classes as such protein synthesis, energy/metabolism and DNA/RNA metabolism, among others.

Conclusions

SNP-markers within genes associated with traits of importance for biomass production were detected. They contribute to characterize the genomic architecture of P. trichocarpa biomass required to support the development and application of marker breeding technologies.

Examining Changes in Stem Taper and Volume Growth with Two-Date 3D Point Clouds

Exact knowledge over tree growth is valuable information for decision makers when considering the purposes of sustainable forest management and planning or optimizing the use of timber, for example. Terrestrial laser scanning (TLS) can be used for measuring tree and forest attributes in very high detail. The study aims at characterizing changes in individual tree attributes (e.g., stem volume growth and taper) during a nine year-long study period in boreal forest conditions. TLS-based three-dimensional (3D) point cloud data were used for identifying and quantifying these changes. The results showed that observing changes in stem volume was possible from TLS point cloud data collected at two different time points. The average volume growth of sample trees was 0.226 m3 during the study period, and the mean relative change in stem volume was 65.0%. In addition, the results of a pairwise Student’s t-test gave strong support (p-value 0.0001) that the used method was able to detect tree growth within the nine-year period between 2008–2017. The findings of this study allow the further development of enhanced methods for TLS-based single tree and forest growth modeling and estimation, which can thus improve the accuracy of forest inventories and offer better tools for future decision-making processes.

Ecological genomics of variation in bud-break phenology and mechanisms of response to climate warming in Populus trichocarpa

Spring bud-break phenology is a critical adaptive feature common to temperate perennial woody plants. Understanding the molecular underpinnings of variation in bud-break is important for elucidating adaptive evolution and predicting outcomes relating to climate change. Field and controlled growth chamber tests were used to assess population-wide patterns in bud-break from wild-sourced black cottonwood (Populus trichocarpa) genotypes. We conducted a genome-wide association study (GWAS) with single nucleotide polymorphisms (SNPs) derived from whole genome sequencing to test for loci underlying variation in bud-break. Bud-break had a quadratic relationship with latitude, where southern- and northern-most provenances generally broke bud earlier than those from central parts of the species' range. Reduced winter chilling increased population-wide variation in bud-break, whereas greater chilling decreased variation. GWAS uncovered 16 loci associated with bud-break. Phenotypic changes connected with allelic variation were replicated in an independent set of P. trichocarpa trees. Despite phenotypic similarities, genetic profiles between southern- and northern-most genotypes were dissimilar based on our GWAS-identified SNPs. We propose that the GWAS-identified loci underpin the geographical pattern in P. trichocarpa and that variation in bud-break reflects different selection for winter chilling and heat sum accumulation, both of which can be affected by climate warming.

High density SNP mapping and QTL analysis for time of leaf budburst in Corylus avellana L

The growing area of European hazelnut (Corylus avellana L.) is increasing, as well as the number of producing countries, and there is a pressing need for new improved cultivars. Hazelnut conventional breeding process is slow, due to the length of juvenile phase and the high heterozygosity level. The development of genetic linkage maps and the identification of molecular markers tightly linked to QTL (quantitative trait loci) of agronomic interest are essential tools for speeding up the selection of seedlings carrying desired traits through marker-assisted selection. The objectives of this study were to enrich a previous linkage map and confirm QTL related to time of leaf budburst, using an F₁ population obtained by crossing Tonda Gentile delle Langhe with Merveille de Bollwiller. Genotyping-by-Sequencing was used to identify a total of 9,999 single nucleotide polymorphism markers. Well saturated linkage maps were constructed for each parent using the double pseudo-testcross mapping strategy. A reciprocal translocation was detected in Tonda Gentile delle Langhe between two non-homologous chromosomes. Applying a bioinformatic approach, we were able to disentangle ‘pseudo-linkage’ between markers, removing markers around the translocation breakpoints and obtain a linear order of the markers for the two chromosomes arms, for each linkage group involved in the translocation. Twenty-nine QTL for time of leaf budburst were identified, including a stably expressed region on LG_02 of the Tonda Gentile delle Langhe map. The stability of these QTL and their coding sequence content indicates promise for the identification of specific chromosomal regions carrying key genes involved in leaf budburst.

MVQTLCIM: composite interval mapping of multivariate traits in a hybrid F1 population of outbred species

Background

With the plummeting cost of the next-generation sequencing technologies, high-density genetic linkage maps could be constructed in a forest hybrid F₁ population. However, based on such genetic maps, quantitative trait loci (QTL) mapping cannot be directly conducted with traditional statistical methods or tools because the linkage phase and segregation pattern of molecular markers are not always fixed as in inbred lines.

Results

We implemented the traditional composite interval mapping (CIM) method to multivariate trait data in forest trees and developed the corresponding software, mvqtlcim. Our method not only incorporated the various segregations and linkage phases of molecular markers, but also applied Takeuchi’s information criterion (TIC) to discriminate the QTL segregation type among several possible alternatives. QTL mapping was performed in a hybrid F₁ population of Populus deltoides and P. simonii, and 12 QTLs were detected for tree height over 6 time points. The software package allowed many options for parameters as well as parallel computing for permutation tests. The features of the software were demonstrated with the real data analysis and a large number of Monte Carlo simulations.

Conclusions

We provided a powerful tool for QTL mapping of multiple or longitudinal traits in an outbred F₁ population, in which the traditional software for QTL mapping cannot be used. This tool will facilitate studying of QTL mapping and thus will accelerate molecular breeding programs especially in forest trees. The tool package is freely available from https://github.com/tongchf /mvqtlcim.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1908-1) contains supplementary material, which is available to authorized users.

Genetics of phenotypic plasticity and biomass traits in hybrid willows across contrasting environments and years

BACKGROUND AND AIMS

Phenotypic plasticity can affect the geographical distribution of taxa and greatly impact the productivity of crops across contrasting and variable environments. The main objectives of this study were to identify genotype-phenotype associations in key biomass and phenology traits and the strength of phenotypic plasticity of these traits in a short-rotation coppice willow population across multiple years and contrasting environments to facilitate marker-assisted selection for these traits.

METHODS

A hybrid Salix viminalis  × ( S. viminalis × Salix schwerinii ) population with 463 individuals was clonally propagated and planted in three common garden experiments comprising one climatic contrast between Sweden and Italy and one water availability contrast in Italy. Several key phenotypic traits were measured and phenotypic plasticity was estimated as the trait value difference between experiments. Quantitative trait locus (QTL) mapping analyses were conducted using a dense linkage map and phenotypic effects of S. schwerinii haplotypes derived from detected QTL were assessed.

KEY RESULTS

Across the climatic contrast, clone predictor correlations for biomass traits were low and few common biomass QTL were detected. This indicates that the genetic regulation of biomass traits was sensitive to environmental variation. Biomass QTL were, however, frequently shared across years and across the water availability contrast. Phenology QTL were generally shared between all experiments. Substantial phenotypic plasticity was found among the hybrid offspring, that to a large extent had a genetic origin. Individuals carrying influential S. schwerinii haplotypes generally performed well in Sweden but less well in Italy in terms of biomass production.

CONCLUSIONS

The results indicate that specific genetic elements of S. schwerinii are more suited to Swedish conditions than to those of Italy. Therefore, selection should preferably be conducted separately for such environments in order to maximize biomass production in admixed S. viminalis × S. schwerinii populations.

Quantitative trait loci influencing forking defects in an outbred pedigree of loblolly pine

Background

The use of wood as an industrial raw material has led to development of plantation forestry, in which trees are planted, managed, and harvested as crops. The productivity of such plantations often exceeds that of less-intensively-managed forests, and land managers have the option of choosing specific planting stock to produce specific types of wood for industrial use. Stem forking, or division of the stem into two or more stems of roughly equal size, is a character trait important in determining the quality of the stem for production of solid wood products. This trait typically has very low individual-tree heritability, but can be more accurately assessed in clonally-replicated plantings where each genotype is represented by several individual trees. We report results from a quantitative trait mapping experiment in a clonally-replicated full-sibling family of loblolly pine (Pinus taeda L.).

Results

Quantitative trait loci influencing forking defects were identified in an outbred full-sibling family of loblolly pine, using single-nucleotide polymorphism markers. Genetic markers in this family segregated either in 1:2:1 (F2 intercross-like segregation) or 1:1 ratio (backcross-like segregation). An integrated linkage map combining markers with different segregation ratios was assembled for this full-sib family, and a total of 409 SNP markers were mapped on 12 linkage groups, covering 1622 cM. Two and three trait loci were identified for forking and ramicorn branch traits, respectively, using the interval mapping method. Three trait loci were detected for both traits using multiple-trait analysis.

Conclusions

The detection of three loci for forking and ramicorn branching in a multiple-trait analysis could mean that there are genes with pleiotropic effects on both traits, or that separate genes affecting different traits are clustered together. The detection of genetic loci associated with variation in stem quality traits in this study supports the hypothesis that marker-assisted selection can be used to decrease the rate of stem defects in breeding populations of loblolly pine.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-016-0446-6) contains supplementary material, which is available to authorized users.

Integrating mixed-effect models into an architectural plant model to simulate inter- and intra-progeny variability: a case study on oil palm (Elaeis guineensis Jacq.)

Three-dimensional (3D) reconstruction of plants is time-consuming and involves considerable levels of data acquisition. This is possibly one reason why the integration of genetic variability into 3D architectural models has so far been largely overlooked. In this study, an allometry-based approach was developed to account for architectural variability in 3D architectural models of oil palm (Elaeis guineensis Jacq.) as a case study. Allometric relationships were used to model architectural traits from individual leaflets to the entire crown while accounting for ontogenetic and morphogenetic gradients. Inter- and intra-progeny variabilities were evaluated for each trait and mixed-effect models were used to estimate the mean and variance parameters required for complete 3D virtual plants. Significant differences in leaf geometry (petiole length, density of leaflets, and rachis curvature) and leaflet morphology (gradients of leaflet length and width) were detected between and within progenies and were modelled in order to generate populations of plants that were consistent with the observed populations. The application of mixed-effect models on allometric relationships highlighted an interesting trade-off between model accuracy and ease of defining parameters for the 3D reconstruction of plants while at the same time integrating their observed variability. Future research will be dedicated to sensitivity analyses coupling the structural model presented here with a radiative balance model in order to identify the key architectural traits involved in light interception efficiency.

Ecological Consequences of Clonal Integration in Plants

Clonal plants are widespread throughout the plant kingdom and dominate in diverse habitats. Spatiotemporal heterogeneity of environment is pervasive at multiple scales, even at scales relevant to individual plants. Clonal integration refers to resource translocation and information communication among the ramets of clonal plants. Due to clonal integration, clonal plant species possess a series of peculiar attributes: plasticity in response to local and non-local conditions, labor division with organ specialization for acquiring locally abundant resources, foraging behavior by selective placement of ramets in resource-rich microhabitats, and avoidance of intraclonal competition. Clonal integration has very profound ecological consequences for clonal plants. It allows them to efficiently cope with environmental heterogeneity, by alleviating local resource shortages, buffering environmental stresses and disturbances, influencing competitive ability, increasing invasiveness, and altering species composition and invasibility at the community level. In this paper, we present a comprehensive review of research on the ecological consequences of plant clonal integration based on a large body of literature. We also attempt to propose perspectives for future research.

Ecological Consequences of Clonal Integration in Plants

Clonal plants are widespread throughout the plant kingdom and dominate in diverse habitats. Spatiotemporal heterogeneity of environment is pervasive at multiple scales, even at scales relevant to individual plants. Clonal integration refers to resource translocation and information communication among the ramets of clonal plants. Due to clonal integration, clonal plant species possess a series of peculiar attributes: plasticity in response to local and non-local conditions, labor division with organ specialization for acquiring locally abundant resources, foraging behavior by selective placement of ramets in resource-rich microhabitats, and avoidance of intraclonal competition. Clonal integration has very profound ecological consequences for clonal plants. It allows them to efficiently cope with environmental heterogeneity, by alleviating local resource shortages, buffering environmental stresses and disturbances, influencing competitive ability, increasing invasiveness, and altering species composition and invasibility at the community level. In this paper, we present a comprehensive review of research on the ecological consequences of plant clonal integration based on a large body of literature. We also attempt to propose perspectives for future research.

QTLs for Woolly Poplar Aphid (Phloeomyzus passerinii L.) Resistance Detected in an Inter-Specific Populus deltoides x P. nigra Mapping Population

The genus Populus represents one of the most economically important groups of forest trees. It is composed by approximately 30 species used for wood and non-wood products, phytoremediation and biomass. Poplar is subjected to several biological and environmental threats although, compared to annual crops, we know far less about the genetic bases of biotic stress resistance. Woolly poplar aphid (Phloeomyzus passerinii) is considered a main pest of cultivated poplars in European and American countries. In this work we present two high density linkage maps in poplar obtained by a genotyping by sequencing (GBS) approach and the identification of QTLs involved in Ph. passerinii resistance. A total of 5,667 polymorphic markers (5,606 SNPs and 61 SSRs) identified on expressed sequences have been used to genotype 131 plants of an F1 population P ×canadensis obtained by an interspecific mate between Populus deltoides (resistant to woolly poplar aphid) and Populus nigra (susceptible to woolly poplar aphid). The two linkage maps, obtained following the two-way pseudo-testcross mapping strategy, have been used to investigate the genetic bases of woolly poplar aphid resistance. One major QTL and two QTLs with minor effects (mapped on LGV, LGXVI and LG XIX) explaining the 65.8% of the genetic variance observed in the progeny in response to Ph. passerinii attack were found. The high density coverage of functional markers allowed the identification of three genes belonging to disease resistance pathway as putative candidates for P. deltoides resistance to woolly poplar aphid. This work is the first report on genetic of woolly poplar aphid genetic resistance and the resistant loci associated markers identified represent a valuable tool in resistance poplar breeding programs.

Biomass traits and candidate genes for bioenergy revealed through association genetics in coppiced European Populus nigra (L.)

BACKGROUND

Second generation (2G) bioenergy from lignocellulosic feedstocks has the potential to develop as a sustainable source of renewable energy; however, significant hurdles still remain for large-scale commercialisation. Populus is considered as a promising 2G feedstock and understanding the genetic basis of biomass yield and feedstock quality are a research priority in this model tree species.

RESULTS

We report the first coppiced biomass study for 714 members of a wide population of European black poplar (Populus nigra L.), a native European tree, selected from 20 river populations ranging in latitude and longitude between 40.5 and 52.1°N and 1.0 and 16.4°E, respectively. When grown at a single site in southern UK, significant Site of Origin (SO) effects were seen for 14 of the 15 directly measured or derived traits including biomass yield, leaf area and stomatal index. There was significant correlation (p < 0.001) between biomass yield traits over 3 years of harvest which identified leaf size and cell production as strong predictors of biomass yield. A 12 K Illumina genotyping array (constructed from 10,331 SNPs in 14 QTL regions and 4648 genes) highlighted significant population genetic structure with pairwise FST showing strong differentiation (p < 0.001) between the Spanish and Italian subpopulations. Robust associations reaching genome-wide significance are reported for main stem height and cell number per leaf; two traits tightly linked to biomass yield. These genotyping and phenotypic data were also used to show the presence of significant isolation by distance (IBD) and isolation by adaption (IBA) within this population.

CONCLUSIONS

The three associations identified reaching genome-wide significance at p < 0.05 include a transcription factor; a putative stress response gene and a gene of unknown function. None of them have been previously linked to bioenergy yield; were shown to be differentially expressed in a panel of three selected genotypes from the collection and represent exciting, novel candidates for further study in a bioenergy tree native to Europe and Euro-Asia. A further 26 markers (22 genes) were found to reach putative significance and are also of interest for biomass yield, leaf area, epidermal cell expansion and stomatal patterning. This research on European P. nigra provides an important foundation for the development of commercial native trees for bioenergy and for advanced, molecular breeding in these species.

Vulnerability to drought-induced cavitation in poplars: synthesis and future opportunities

Vulnerability to drought-induced cavitation is a key trait of plant water relations. Here, we summarize the available literature on vulnerability to drought-induced cavitation in poplars (Populus spp.), a genus of agronomic, ecological and scientific importance. Vulnerability curves and vulnerability parameters (including the water potential inducing 50% loss in hydraulic conductivity, P50) were collected from 37 studies published between 1991 and 2014, covering a range of 10 species and 12 interspecific hybrid crosses. Results of our meta-analysis confirm that poplars are among the most vulnerable woody species to drought-induced cavitation (mean P50  = -1.44 and -1.55 MPa across pure species and hybrids, respectively). Yet, significant variation occurs among species (P50 range: 1.43 MPa) and among hybrid crosses (P50 range: 1.12 MPa), within species and hybrid crosses (max. P50 range reported: 0.8 MPa) as well as in response to environmental factors including nitrogen fertilization, irradiance, temperature and drought (max. P50 range reported: 0.75 MPa). Potential implications and gaps in knowledge are discussed in the context of poplar cultivation, species adaptation and climate modifications. We suggest that poplars represent a valuable model for studies on drought-induced cavitation, especially to elucidate the genetic and molecular basis of cavitation resistance in Angiosperms.

Geographical and environmental gradients shape phenotypic trait variation and genetic structure in Populus trichocarpa

• Populus trichocarpa is widespread across western North America spanning extensive variation in photoperiod, growing season and climate. We investigated trait variation in P. trichocarpa using over 2000 trees from a common garden at Vancouver, Canada, representing replicate plantings of 461 genotypes originating from 136 provenance localities. • We measured 40 traits encompassing phenological events, biomass accumulation, growth rates, and leaf, isotope and gas exchange-based ecophysiology traits. With replicated plantings and 29,354 single nucleotide polymorphisms (SNPs) from 3518 genes, we estimated both broad-sense trait heritability (H(2)) and overall population genetic structure from principal component analysis. • Populus trichocarpa had high phenotypic variation and moderate/high H(2) for many traits. H(2) ranged from 0.3 to 0.9 in phenology, 0.3 to 0.8 in biomass and 0.1 to 0.8 in ecophysiology traits. Most traits correlated strongly with latitude, maximum daylength and temperature of tree origin, but not necessarily with elevation, precipitation or heat : moisture indices. Trait H(2) values reflected trait correlation strength with geoclimate variables. The population genetic structure had one significant principal component (PC1) which correlated with daylength and showed enrichment for genes relating to circadian rhythm and photoperiod. • Robust relationships between traits, population structure and geoclimate in P. trichocarpa reflect patterns which suggest that range-wide geographical and environment gradients have shaped its genotypic and phenotypic variability.

A molecular genetic linkage map of Eucommia ulmoides and quantitative trait loci (QTL) analysis for growth traits

Eucommia ulmoides is an economically important tree species for both herbal medicine and organic chemical industry. Effort to breed varieties with improved yield and quality is limited by the lack of knowledge on the genetic basis of the traits. A genetic linkage map of E. ulmoides was constructed from a full-sib family using sequence-related amplified polymorphism, amplified fragment length polymorphism, inter-simple sequence repeat and simple sequence repeat markers. In total, 706 markers were mapped in 25 linkage groups covering 2133 cM. The genetic linkage map covered approximately 89% of the estimated E. ulmoides genome with an average of 3.1 cM between adjacent markers. The present genetic linkage map was used to identify quantitative trait loci (QTL) affecting growth-related traits. Eighteen QTLs were found to explain 12.4%–33.3% of the phenotypic variance. This genetic linkage map provides a tool for marker-assisted selection and for studies of genome in E. ulmoides.

Genetic architecture of spring and autumn phenology in Salix

BACKGROUND

In woody plants from temperate regions, adaptation to the local climate results in annual cycles of growth and dormancy, and optimal regulation of these cycles are critical for growth, long-term survival, and competitive success. In this study we have investigated the genetic background to growth phenology in a Salix pedigree by assessing genetic and phenotypic variation in growth cessation, leaf senescence and bud burst in different years and environments. A previously constructed linkage map using the same pedigree and anchored to the annotated genome of P. trichocarpa was improved in target regions and used for QTL analysis of the traits. The major aims in this study were to map QTLs for phenology traits in Salix, and to identify candidate genes in QTL hot spots through comparative mapping with the closely related Populus trichocarpa.

RESULTS

All traits varied significantly among genotypes and the broad-sense heritabilities ranged between 0.5 and 0.9, with the highest for leaf senescence. In total across experiment and years, 80 QTLs were detected. For individual traits, the QTLs explained together from 21.5 to 56.5% of the variation. Generally each individual QTL explained a low amount of the variation but three QTLs explained above 15% of the variation with one QTL for leaf senescence explaining 34% of the variation. The majority of the QTLs were recurrently identified across traits, years and environments. Two hotspots were identified on linkage group (LG) II and X where narrow QTLs for all traits co-localized.

CONCLUSIONS

This study provides the most detailed analysis of QTL detection for phenology in Salix conducted so far. Several hotspot regions were found where QTLs for different traits and QTLs for the same trait but identified during different years co-localised. Many QTLs co-localised with QTLs found in poplar for similar traits that could indicate common pathways for these traits in Salicaceae. This study is an important first step in identifying QTLs and candidate genes for phenology traits in Salix.

A molecular framework for seasonal growth-dormancy regulation in perennial plants

The timing of the onset and release of dormancy impacts the survival, productivity and spatial distribution of temperate horticultural and forestry perennials and is mediated by at least three main regulatory programs involving signal perception and processing by phytochromes (PHYs) and PHY-interacting transcription factors (PIFs). PIF4 functions as a key regulator of plant growth in response to both external and internal signals. In poplar, the expression of PIF4 and PIF3-LIKE1 is upregulated in response to short days, while PHYA and PHYB are not regulated at the transcriptional level. Integration of light and environmental signals is achieved by gating the expression and transcriptional activity of PIF4. During this annual cycle, auxin promotes the degradation of Aux/IAA transcriptional repressors through the SKP-Cullin-F-boxTIR1 complex, relieving the repression of auxin-responsive genes by allowing auxin response factors (ARFs) to activate the transcription of auxin-responsive genes involved in growth responses. Analyses of transcriptome changes during dormancy transitions have identified MADS-box transcription factors associated with endodormancy induction. Previous studies show that poplar dormancy-associated MADS-box (DAM) genes PtMADS7 and PtMADS21 are differentially regulated during the growth-dormancy cycle. Endodormancy may be regulated by internal factors, which are specifically localized in buds. PtMADS7/PtMADS21 may function as an internal regulator in poplar. The control of flowering time shares certain regulatory hierarchies with control of the dormancy/growth cycle. However, the particularities of different stages of the dormancy/growth cycle warrant comprehensive approaches to identify the causative genes for the entire cycle. A growing body of knowledge also indicates epigenetic regulation plays a role in these processes in perennial horticultural and forestry plants. The increased knowledge contributes to better understanding of the dormancy process and consequently to precise manipulation of dormancy-related horticultural traits, such as flowering time.

Insight into the genetic components of community genetics: QTL mapping of insect association in a fast-growing forest tree

Identifying genetic sequences underlying insect associations on forest trees will improve the understanding of community genetics on a broad scale. We tested for genomic regions associated with insects in hybrid poplar using quantitative trait loci (QTL) analyses conducted on data from a common garden experiment. The F2 offspring of a hybrid poplar (Populus trichocarpa x P. deltoides) cross were assessed for seven categories of insect leaf damage at two time points, June and August. Positive and negative correlations were detected among damage categories and between sampling times. For example, sap suckers on leaves in June were positively correlated with sap suckers on leaves (P<0.001) but negatively correlated with skeletonizer damage (P<0.01) in August. The seven forms of leaf damage were used as a proxy for seven functional groups of insect species. Significant variation in insect association occurred among the hybrid offspring, including transgressive segregation of susceptibility to damage. NMDS analyses revealed significant variation and modest broad-sense heritability in insect community structure among genets. QTL analyses identified 14 genomic regions across 9 linkage groups that correlated with insect association. We used three genomics tools to test for putative mechanisms underlying the QTL. First, shikimate-phenylpropanoid pathway genes co-located to 9 of the 13 QTL tested, consistent with the role of phenolic glycosides as defensive compounds. Second, two insect association QTL corresponded to genomic hotspots for leaf trait QTL as identified in previous studies, indicating that, in addition to biochemical attributes, leaf morphology may influence insect preference. Third, network analyses identified categories of gene models over-represented in QTL for certain damage types, providing direction for future functional studies. These results provide insight into the genetic components involved in insect community structure in a fast-growing forest tree.

Deciphering structural and temporal interplays during the architectural development of mango trees

Plant architecture is commonly defined by the adjacency of organs within the structure and their properties. Few studies consider the effect of endogenous temporal factors, namely phenological factors, on the establishment of plant architecture. This study hypothesized that, in addition to the effect of environmental factors, the observed plant architecture results from both endogenous structural and temporal components, and their interplays. Mango tree, which is characterized by strong phenological asynchronisms within and between trees and by repeated vegetative and reproductive flushes during a growing cycle, was chosen as a plant model. During two consecutive growing cycles, this study described vegetative and reproductive development of 20 trees submitted to the same environmental conditions. Four mango cultivars were considered to assess possible cultivar-specific patterns. Integrative vegetative and reproductive development models incorporating generalized linear models as components were built. These models described the occurrence, intensity, and timing of vegetative and reproductive development at the growth unit scale. This study showed significant interplays between structural and temporal components of plant architectural development at two temporal scales. Within a growing cycle, earliness of bud burst was highly and positively related to earliness of vegetative development and flowering. Between growing cycles, flowering growth units delayed vegetative development compared to growth units that did not flower. These interplays explained how vegetative and reproductive phenological asynchronisms within and between trees were generated and maintained. It is suggested that causation networks involving structural and temporal components may give rise to contrasted tree architectures.

Transcriptional profiling of bud dormancy induction and release in oak by next-generation sequencing

Background

In temperate regions, the time lag between vegetative bud burst and bud set determines the duration of the growing season of trees (i.e. the duration of wood biomass production). Dormancy, the period during which the plant is not growing, allows trees to avoid cold injury resulting from exposure to low temperatures. An understanding of the molecular machinery controlling the shift between these two phenological states is of key importance in the context of climatic change. The objective of this study was to identify genes upregulated during endo- and ecodormancy, the two main stages of bud dormancy. Sessile oak is a widely distributed European white oak species. A forcing test on young trees was first carried out to identify the period most likely to correspond to these two stages. Total RNA was then extracted from apical buds displaying endo- and ecodormancy. This RNA was used for the generation of cDNA libraries, and in-depth transcriptome characterization was performed with 454 FLX pyrosequencing technology.

Results

Pyrosequencing produced a total of 495,915 reads. The data were cleaned, duplicated reads removed, and sequences were mapped onto the oak UniGene data. Digital gene expression analysis was performed, with both R statistics and the R-Bioconductor packages (edgeR and DESeq), on 6,471 contigs with read numbers ≥ 5 within any contigs. The number of sequences displaying significant differences in expression level (read abundance) between endo- and ecodormancy conditions ranged from 75 to 161, depending on the algorithm used. 13 genes displaying significant differences between conditions were selected for further analysis, and 11 of these genes, including those for glutathione-S-transferase (GST) and dehydrin xero2 (XERO2) were validated by quantitative PCR.

Conclusions

The identification and functional annotation of differentially expressed genes involved in the “response to abscisic acid”, “response to cold stress” and “response to oxidative stress” categories constitutes a major step towards characterization of the molecular network underlying vegetative bud dormancy, an important life history trait of long-lived organisms.

Genome anchored QTLs for biomass productivity in hybrid Populus grown under contrasting environments

Traits related to biomass production were analyzed for the presence of quantitative trait loci (QTLs) in a Populus trichocarpa × P. deltoides F(2) population. A genetic linkage map composed of 841 SSR, AFLP, and RAPD markers and phenotypic data from 310 progeny were used to identify genomic regions harboring biomass QTLs. Twelve intervals were identified, of which BM-1, BM-2, and BM-7 were identified in all three years for both height and diameter. One putative QTL, BM-7, and one suggestive QTL exhibited significant evidence of over-dominance in all three years for both traits. Conversely, QTLs BM-4 and BM-6 exhibited evidence of under-dominance in both environments for height and diameter. Seven of the nine QTLs were successfully anchored, and QTL peak positions were estimated for each one on the P. trichocarpa genome assembly using flanking SSR markers with known physical positions. Of the 3,031 genes located in genome-anchored QTL intervals, 1,892 had PFAM annotations. Of these, 1,313, representing 255 unique annotations, had at least one duplicate copy in a QTL interval identified on a separate scaffold. This observation suggests that some QTLs identified in this study may have shared the same ancestral sequence prior to the salicoid genome duplication in Populus.

Integrating genome annotation and QTL position to identify candidate genes for productivity, architecture and water-use efficiency in Populus spp

BACKGROUND

Hybrid poplars species are candidates for biomass production but breeding efforts are needed to combine productivity and water use efficiency in improved cultivars. The understanding of the genetic architecture of growth in poplar by a Quantitative Trait Loci (QTL) approach can help us to elucidate the molecular basis of such integrative traits but identifying candidate genes underlying these QTLs remains difficult. Nevertheless, the increase of genomic information together with the accessibility to a reference genome sequence (Populus trichocarpa Nisqually-1) allow to bridge QTL information on genetic maps and physical location of candidate genes on the genome. The objective of the study is to identify QTLs controlling productivity, architecture and leaf traits in a P. deltoides x P. trichocarpa F1 progeny and to identify candidate genes underlying QTLs based on the anchoring of genetic maps on the genome and the gene ontology information linked to genome annotation. The strategy to explore genome annotation was to use Gene Ontology enrichment tools to test if some functional categories are statistically over-represented in QTL regions.

RESULTS

Four leaf traits and 7 growth traits were measured on 330 F1 P. deltoides x P. trichocarpa progeny. A total of 77 QTLs controlling 11 traits were identified explaining from 1.8 to 17.2% of the variation of traits. For 58 QTLs, confidence intervals could be projected on the genome. An extended functional annotation was built based on data retrieved from the plant genome database Phytozome and from an inference of function using homology between Populus and the model plant Arabidopsis. Genes located within QTL confidence intervals were retrieved and enrichments in gene ontology (GO) terms were determined using different methods. Significant enrichments were found for all traits. Particularly relevant biological processes GO terms were identified for QTLs controlling number of sylleptic branches: intervals were enriched in GO terms of biological process like 'ripening' and 'adventitious roots development'.

CONCLUSION

Beyond the simple identification of QTLs, this study is the first to use a global approach of GO terms enrichment analysis to fully explore gene function under QTLs confidence intervals in plants. This global approach may lead to identification of new candidate genes for traits of interest.

Integrating genome annotation and QTL position to identify candidate genes for productivity, architecture and water-use efficiency in Populus spp

BACKGROUND

Hybrid poplars species are candidates for biomass production but breeding efforts are needed to combine productivity and water use efficiency in improved cultivars. The understanding of the genetic architecture of growth in poplar by a Quantitative Trait Loci (QTL) approach can help us to elucidate the molecular basis of such integrative traits but identifying candidate genes underlying these QTLs remains difficult. Nevertheless, the increase of genomic information together with the accessibility to a reference genome sequence (Populus trichocarpa Nisqually-1) allow to bridge QTL information on genetic maps and physical location of candidate genes on the genome. The objective of the study is to identify QTLs controlling productivity, architecture and leaf traits in a P. deltoides x P. trichocarpa F1 progeny and to identify candidate genes underlying QTLs based on the anchoring of genetic maps on the genome and the gene ontology information linked to genome annotation. The strategy to explore genome annotation was to use Gene Ontology enrichment tools to test if some functional categories are statistically over-represented in QTL regions.

RESULTS

Four leaf traits and 7 growth traits were measured on 330 F1 P. deltoides x P. trichocarpa progeny. A total of 77 QTLs controlling 11 traits were identified explaining from 1.8 to 17.2% of the variation of traits. For 58 QTLs, confidence intervals could be projected on the genome. An extended functional annotation was built based on data retrieved from the plant genome database Phytozome and from an inference of function using homology between Populus and the model plant Arabidopsis. Genes located within QTL confidence intervals were retrieved and enrichments in gene ontology (GO) terms were determined using different methods. Significant enrichments were found for all traits. Particularly relevant biological processes GO terms were identified for QTLs controlling number of sylleptic branches: intervals were enriched in GO terms of biological process like 'ripening' and 'adventitious roots development'.

CONCLUSION

Beyond the simple identification of QTLs, this study is the first to use a global approach of GO terms enrichment analysis to fully explore gene function under QTLs confidence intervals in plants. This global approach may lead to identification of new candidate genes for traits of interest.

Multiallelic epistatic model for an out-bred cross and mapping algorithm of interactive quantitative trait loci

BACKGROUND

Genetic mapping has proven to be powerful for studying the genetic architecture of complex traits by characterizing a network of the underlying interacting quantitative trait loci (QTLs). Current statistical models for genetic mapping were mostly founded on the biallelic epistasis of QTLs, incapable of analyzing multiallelic QTLs and their interactions that are widespread in an outcrossing population.

RESULTS

Here we have formulated a general framework to model and define the epistasis between multiallelic QTLs. Based on this framework, we have derived a statistical algorithm for the estimation and test of multiallelic epistasis between different QTLs in a full-sib family of outcrossing species. We used this algorithm to genomewide scan for the distribution of multiallelic epistasis for a rooting ability trait in an outbred cross derived from two heterozygous poplar trees. The results from simulation studies indicate that the positions and effects of multiallelic QTLs can well be estimated with a modest sample and heritability.

CONCLUSIONS

The model and algorithm developed provide a useful tool for better characterizing the genetic control of complex traits in a heterozygous family derived from outcrossing species, such as forest trees, and thus fill a gap that occurs in genetic mapping of this group of important but underrepresented species.

Association genetics of wood physical traits in the conifer white spruce and relationships with gene expression

Marker-assisted selection holds promise for highly influencing tree breeding, especially for wood traits, by considerably reducing breeding cycles and increasing selection accuracy. In this study, we used a candidate gene approach to test for associations between 944 single-nucleotide polymorphism markers from 549 candidate genes and 25 wood quality traits in white spruce. A mixed-linear model approach, including a weak but nonsignificant population structure, was implemented for each marker-trait combination. Relatedness among individuals was controlled using a kinship matrix estimated either from the known half-sib structure or from the markers. Both additive and dominance effect models were tested. Between 8 and 21 single-nucleotide polymorphisms (SNPs) were found to be significantly associated (P ≤ 0.01) with each of earlywood, latewood, or total wood traits. After controlling for multiple testing (Q ≤ 0.10), 13 SNPs were still significant across as many genes belonging to different families, each accounting for between 3 and 5% of the phenotypic variance in 10 wood characters. Transcript accumulation was determined for genes containing SNPs associated with these traits. Significantly different transcript levels (P ≤ 0.05) were found among the SNP genotypes of a 1-aminocyclopropane-1-carboxylate oxidase, a β-tonoplast intrinsic protein, and a long-chain acyl-CoA synthetase 9. These results should contribute toward the development of efficient marker-assisted selection in an economically important tree species.

QTL mapping in white spruce: gene maps and genomic regions underlying adaptive traits across pedigrees, years and environments

BACKGROUND

The genomic architecture of bud phenology and height growth remains poorly known in most forest trees. In non model species, QTL studies have shown limited application because most often QTL data could not be validated from one experiment to another. The aim of our study was to overcome this limitation by basing QTL detection on the construction of genetic maps highly-enriched in gene markers, and by assessing QTLs across pedigrees, years, and environments.

RESULTS

Four saturated individual linkage maps representing two unrelated mapping populations of 260 and 500 clonally replicated progeny were assembled from 471 to 570 markers, including from 283 to 451 gene SNPs obtained using a multiplexed genotyping assay. Thence, a composite linkage map was assembled with 836 gene markers.For individual linkage maps, a total of 33 distinct quantitative trait loci (QTLs) were observed for bud flush, 52 for bud set, and 52 for height growth. For the composite map, the corresponding numbers of QTL clusters were 11, 13, and 10. About 20% of QTLs were replicated between the two mapping populations and nearly 50% revealed spatial and/or temporal stability. Three to four occurrences of overlapping QTLs between characters were noted, indicating regions with potential pleiotropic effects. Moreover, some of the genes involved in the QTLs were also underlined by recent genome scans or expression profile studies.Overall, the proportion of phenotypic variance explained by each QTL ranged from 3.0 to 16.4% for bud flush, from 2.7 to 22.2% for bud set, and from 2.5 to 10.5% for height growth. Up to 70% of the total character variance could be accounted for by QTLs for bud flush or bud set, and up to 59% for height growth.

CONCLUSIONS

This study provides a basic understanding of the genomic architecture related to bud flush, bud set, and height growth in a conifer species, and a useful indicator to compare with Angiosperms. It will serve as a basic reference to functional and association genetic studies of adaptation and growth in Picea taxa. The putative QTNs identified will be tested for associations in natural populations, with potential applications in molecular breeding and gene conservation programs. QTLs mapping consistently across years and environments could also be the most important targets for breeding, because they represent genomic regions that may be least affected by G × E interactions.

Forest tree genomics: growing resources and applications

Over the past two decades, research in forest tree genomics has lagged behind that of model and agricultural systems. However, genomic research in forest trees is poised to enter into an important and productive phase owing to the advent of next-generation sequencing technologies, the enormous genetic diversity in forest trees and the need to mitigate the effects of climate change. Research on long-lived woody perennials is extending our molecular knowledge of complex life histories and adaptations to the environment - enriching a field that has traditionally drawn biological inference from a few short-lived herbaceous species.

Climate-driven local adaptation of ecophysiology and phenology in balsam poplar, Populus balsamifera L. (Salicaceae)

PREMISE OF THE STUDY

During past episodes of climate change, many plant species experienced large-scale range expansions. Expanding populations likely encountered strong selection as they colonized new environments. In this study we examine the extent to which populations of the widespread forest tree Populus balsamifera L. have become locally adapted as the species expanded into its current range since the last glaciation.

METHODS

We tested for adaptive variation in 13 ecophysiology and phenology traits on clonally propagated genotypes originating from a range-wide sample of 20 subpopulations. The hypothesis of local adaption was tested by comparing among-population variation at ecologically important traits (Q(ST)) to expected variation based on demographic history (F(ST)) estimated from a large set of nuclear single nucleotide polymorphism loci.

KEY RESULTS

Evidence for divergence in excess of neutral expectations was present for eight of 13 traits. Bud phenology, petiole length, and leaf nitrogen showed the greatest divergence (all Q(ST) > 0.6), whereas traits related to leaf water usage showed the least (all Q(ST) ≤ 0.30) and were not different from neutrality. Strong correlations were present between traits, geography, and climate, and they revealed a general pattern of northern subpopulations adapted to shorter, drier growing seasons compared with populations in the center or eastern regions of the range.

CONCLUSIONS

Our study demonstrates pronounced adaptive variation in ecophysiology and phenology among balsam poplar populations. These results suggest that as this widespread forest tree species expanded its range since the end of the last glacial maximum, it evolved rapidly in response to geographically variable selection.

Functional mapping of quantitative trait loci associated with rice tillering

Several biologically significant parameters that are related to rice tillering are closely associated with rice grain yield. Although identification of the genes that control rice tillering and therefore influence crop yield would be valuable for rice production management and genetic improvement, these genes remain largely unidentified. In this study, we carried out functional mapping of quantitative trait loci (QTLs) for rice tillering in 129 doubled haploid lines, which were derived from a cross between IR64 and Azucena. We measured the average number of tillers in each plot at seven developmental stages and fit the growth trajectory of rice tillering with the Wang-Lan-Ding mathematical model. Four biologically meaningful parameters in this model--the potential maximum for tiller number (K), the optimum tiller time (t(0)), and the increased rate (r), or the reduced rate (c) at the time of deviation from t(0)--were our defined variables for multi-marker joint analysis under the framework of penalized maximum likelihood, as well as composite interval mapping. We detected a total of 27 QTLs that accounted for 2.49-8.54% of the total phenotypic variance. Nine common QTLs across multi-marker joint analysis and composite interval mapping showed high stability, while one QTL was environment-specific and three were epistatic. We also identified several genomic segments that are associated with multiple traits. Our results describe the genetic basis of rice tiller development, enable further marker-assisted selection in rice cultivar development, and provide useful information for rice production management.

Unconditional and conditional QTL mapping for the developmental behavior of tiller number in rice (Oryza sativa L.)

A single segment substitution population of 26 lines and their recipient parent Hua-jing-xian 74 (HJX74) were selected as experimental materials for analyzing the developmental behavior of tiller number in rice. By the unconditional QTL (quantitative trait locus) mapping method, a total number of 14 SSSLs were detected with QTLs controlling rice tiller number. The number of QTLs significantly affecting tiller number and their effect values estimated differed across measuring stages. More QTLs could be detected based on time-dependent measures of different stages. By the conditional QTL mapping method, it is possible to reveal net expression of gene in a time interval. 14 QTLs on tiller number expressed their effects in dynamic patterns of themselves during whole ontogeny. They exhibited mainly negative effects within 7 days after transplanting. During 7-21 days, QTLs were in active status and expressed larger positive effects. In the mid-period of 21-35 days, they had opposite genetic effects to wither tillers. Since then these QTLs expressed positive effects again to cause the appearance of noneffective tillers. The dynamics of QTL effects was in agreement with the actual change of tillers. Mapping QTL combining unconditional with conditional analysis for time-dependent measures is helpful to understand roundly the genetic bases for the development of quantitative traits.

Productivity, water-use efficiency and tolerance to moderate water deficit correlate in 33 poplar genotypes from a Populus deltoides x Populus trichocarpa F1 progeny

Genotypic variability for productivity, water-use efficiency and leaf traits in 33 genotypes selected from an F1 progeny of Populus deltoides Bartr. ex Marsh x Populus trichocarpa L. was explored under optimal and moderate water-deficit conditions. Saplings of the 33 genotypes were grown in a two-plot open field at INRA Orléans (France) and coppiced every year. A moderate water deficit was induced during two successive years on one plot by withholding irrigation, while the second one remained irrigated (control). Stem biomass and leaf structure (e.g., specific leaf area and leaf area) were measured in 2004 and 2005 and functional leaf traits (e.g., carbon isotope discrimination, Delta) were measured only in 2004. Tolerance to water deficit was estimated at genotype level as the ability to limit losses in biomass production in water deficit versus control trees. Stem biomass, leaf structure and Delta displayed a significant genotypic variability whatever the irrigation regime. For all traits, genotype ranks remained stable across years for similar irrigation conditions. Carbon isotope discrimination scaled negatively with productivity and leaf nitrogen content in controls. The most productive genotypes were the least tolerant to moderate water deficit. No relationship was evidenced between Delta and the level of tolerance to water deficit. The relationships between traits evidenced in this collection of P. deltoides x P. trichocarpa F1 genotypes contrast with the ones that were previously detected in a collection of P. deltoides x Populus nigra L. cultivars tested in the same field trial.