Distribution of genomic regions differentiating oak species assessed by QTL detection

Genotype-by-environment and QTL-by-environment interactions in sweet cherry (Prunus avium L.) for flowering date

In sweet cherry (Prunus avium L.), flowering date is strongly dependent on the environment conditions and, therefore, is a trait of major interest for adaptation to climate change. Such trait can be influenced by genotype-by-environment interaction (G×E), that refers to differences in the response of genotypes to different environments. If not taken into account, G×E can reduce selection accuracy and overall genetic gain. However, little is known about G×E in fruit tree species. Flowering date is a highly heritable and polygenic trait for which many quantitative trait loci (QTLs) have been identified. As for the overall genetic performance, differential expression of QTLs in response to environment (QTL-by-environment interaction, QTL×E) can occur. The present study is based on the analysis of a multi-environment trial (MET) suitable for the study of G×E and QTL×E in sweet cherry. It consists of a sweet cherry F₁ full-sib family (n = 121) derived from the cross between cultivars 'Regina' and 'Lapins' and planted in two copies in five locations across four European countries (France, Italy, Slovenia and Spain) covering a large range of climatic conditions. The aim of this work was to study the effect of the environment on flowering date and estimate G×E, to carry QTL detection in different environments in order to study the QTL stability across environments and to estimate QTL×E. A strong effect of the environment on flowering date and its genetic control was highlighted. Two large-effect and environment-specific QTLs with significant QTL×E were identified on linkage groups (LGs) 1 and 4. This work gives new insights into the effect of the environment on a trait of main importance in one of the most economically important fruit crops in temperate regions. Moreover, molecular markers were developed for flowering date and a strategy consisting in using specific markers for warm or cold regions was proposed to optimize marker-assisted selection (MAS) in sweet cherry breeding programs.

Quercus Conservation Genetics and Genomics: Past, Present, and Future

Quercus species (oaks) have been an integral part of the landscape in the northern hemisphere for millions of years. Their ability to adapt and spread across different environments and their contributions to many ecosystem services is well documented. Human activity has placed many oak species in peril by eliminating or adversely modifying habitats through exploitative land usage and by practices that have exacerbated climate change. The goal of this review is to compile a list of oak species of conservation concern, evaluate the genetic data that is available for these species, and to highlight the gaps that exist. We compiled a list of 124 Oaks of Concern based on the Red List of Oaks 2020 and the Conservation Gap Analysis for Native U.S. Oaks and their evaluations of each species. Of these, 57% have been the subject of some genetic analysis, but for most threatened species (72%), the only genetic analysis was done as part of a phylogenetic study. While nearly half (49%) of published genetic studies involved population genetic analysis, only 16 species of concern (13%) have been the subject of these studies. This is a critical gap considering that analysis of intraspecific genetic variability and genetic structure are essential for designing conservation management strategies. We review the published population genetic studies to highlight their application to conservation. Finally, we discuss future directions in Quercus conservation genetics and genomics.

Quercus species divergence is driven by natural selection on evolutionarily less integrated traits

Functional traits are organismal attributes that can respond to environmental cues, thereby providing important ecological functions. In addition, an organism's potential for adaptation is defined by the patterns of covariation among groups of functionally related traits. Whether an organism is evolutionarily constrained or has the potential for adaptation is based on the phenotypic integration or modularity of these traits. Here, we revisited leaf morphology in two European sympatric white oaks (Quercus petraea (Matt.) Liebl. and Quercus robur L.), sampling 2098 individuals, across much of their geographical distribution ranges. At the phenotypic level, leaf morphology traditionally encompasses discriminant attributes among different oak species. Here, we estimated in situ heritability, genetic correlation, and integration across such attributes. Also, we performed Selection Response Decomposition to test these traits for potential differences in oak species' evolutionary responses. Based on the uncovered functional units of traits (modules) in our study, the morphological module "leaf size gradient" was highlighted among functionally integrated traits. Equally, this module was defined in both oaks as being under "global regulation" in vegetative bud establishment and development. Lamina basal shape and intercalary veins' number were not, or, less integrated within the initially defined leaf functional unit, suggesting more than one module within the leaf traits' ensemble. Since these traits generally show the greatest species discriminatory power, they potentially underwent effective differential response to selection among oaks. Indeed, the selection of these traits could have driven the ecological preferences between the two sympatric oaks growing under different microclimates.

Association and linkage mapping to unravel genetic architecture of phenological traits and lateral bearing in Persian walnut (Juglans regia L.)

BACKGROUND

Unravelling the genetic architecture of agronomic traits in walnut such as budbreak date and bearing habit, is crucial for climate change adaptation and yield improvement. A Genome-Wide Association Study (GWAS) using multi-locus models was conducted in a panel of 170 walnut accessions genotyped using the Axiom™ J. regia 700 K SNP array, with phenological data from 2018, 2019 and legacy data. These accessions come from the INRAE walnut germplasm collection which is the result of important prospecting work performed in many countries around the world. In parallel, an F1 progeny of 78 individuals segregating for phenology-related traits, was genotyped with the same array and phenotyped for the same traits, to construct linkage maps and perform Quantitative Trait Loci (QTLs) detection.

RESULTS

Using GWAS, we found strong associations of SNPs located at the beginning of chromosome 1 with both budbreak and female flowering dates. These findings were supported by QTLs detected in the same genomic region. Highly significant associated SNPs were also detected using GWAS for heterodichogamy and lateral bearing habit, both on chromosome 11. We developed a Kompetitive Allele Specific PCR (KASP) marker for budbreak date in walnut, and validated it using plant material from the Walnut Improvement Program of the University of California, Davis, demonstrating its effectiveness for marker-assisted selection in Persian walnut. We found several candidate genes involved in flowering events in walnut, including a gene related to heterodichogamy encoding a sugar catabolism enzyme and a cell division related gene linked to female flowering date.

CONCLUSIONS

This study enhances knowledge of the genetic architecture of important agronomic traits related to male and female flowering processes and lateral bearing in walnut. The new marker available for budbreak date, one of the most important traits for good fruiting, will facilitate the selection and development of new walnut cultivars suitable for specific climates.

ddRAD Sequencing-Based Identification of Genomic Boundaries and Permeability in Quercus ilex and Q. suber Hybrids

Hybridization and its relevance is a hot topic in ecology and evolutionary biology. Interspecific gene flow may play a key role in species adaptation to environmental change, as well as in the survival of endangered populations. Despite the fact that hybridization is quite common in plants, many hybridizing species, such as Quercus spp., maintain their integrity, while precise determination of genomic boundaries between species remains elusive. Novel high throughput sequencing techniques have opened up new perspectives in the comparative analysis of genomes and in the study of historical and current interspecific gene flow. In this work, we applied ddRADseq technique and developed an ad hoc bioinformatics pipeline for the study of ongoing hybridization between two relevant Mediterranean oaks, Q. ilex and Q. suber. We adopted a local scale approach, analyzing adult hybrids (sensu lato) identified in a mixed stand and their open-pollinated progenies. We have identified up to 9,435 markers across the genome and have estimated individual introgression levels in adults and seedlings. Estimated contribution of Q. suber to the genome is higher, on average, in hybrid progenies than in hybrid adults, suggesting preferential backcrossing with this parental species, maybe followed by selection during juvenile stages against individuals with higher Q. suber genomic contribution. Most discriminating markers seem to be scattered throughout the genome, suggesting that a large number of small genomic regions underlie boundaries between these species. A noticeable proportion of the markers (26%) showed allelic frequencies in adult hybrids very similar to one of the parental species, and very different from the other; a finding that seems relevant for understanding the hybridization process and the occurrence of adaptive introgression. Candidate marker databases developed in this study constitute a valuable resource to design large scale re-sequencing experiments in Mediterranean sclerophyllous oak species and could provide insight in species boundaries and on adaptive introgression between Q. suber and Q. ilex.

Heritability and genetic architecture of reproduction-related traits in a temperate oak species

Reproduction, one of the main components of plant fitness, is highly variable in response to environmental cues, but little is known about the genetic determinism underlying reproduction-related traits in forest tree species. There is therefore an urgent need to characterize the genetic architecture of those traits if we are to predict the evolutionary trajectories of forest populations facing rapidly changing environment and mitigate their impacts. Using a full-sib family of pedunculate oak (Quercus robur), we investigated the within population variability of seed production and mean seed mass during four consecutive years. Reproductive traits were highly variable between trees and between years. The high narrow sense heritability and evolvability estimated underline the important genetic effect on the variability in seed production and mean seed mass. Despite a large variability over years, reproductive traits show significant genetic correlation between years. Furthermore, for the first time in forest tree species, quantitative trait loci (QTLs) associated with seed production and mean mass of a seed have been identified. While it is commonly assumed and observed that fitness-traits have low narrow sense heritabilities, our findings show that reproduction-related traits may undergo evolutionary changes under selective pressure and may be determinant for tree adaptation.

Oak genome reveals facets of long lifespan

Oaks are an important part of our natural and cultural heritage. Not only are they ubiquitous in our most common landscapes¹ but they have also supplied human societies with invaluable services, including food and shelter, since prehistoric times². With 450 species spread throughout Asia, Europe and America³, oaks constitute a critical global renewable resource. The longevity of oaks (several hundred years) probably underlies their emblematic cultural and historical importance. Such long-lived sessile organisms must persist in the face of a wide range of abiotic and biotic threats over their lifespans. We investigated the genomic features associated with such a long lifespan by sequencing, assembling and annotating the oak genome. We then used the growing number of whole-genome sequences for plants (including tree and herbaceous species) to investigate the parallel evolution of genomic characteristics potentially underpinning tree longevity. A further consequence of the long lifespan of trees is their accumulation of somatic mutations during mitotic divisions of stem cells present in the shoot apical meristems. Empirical⁴ and modelling⁵ approaches have shown that intra-organismal genetic heterogeneity can be selected for⁶ and provides direct fitness benefits in the arms race with short-lived pests and pathogens through a patchwork of intra-organismal phenotypes⁷. However, there is no clear proof that large-statured trees consist of a genetic mosaic of clonally distinct cell lineages within and between branches. Through this case study of oak, we demonstrate the accumulation and transmission of somatic mutations and the expansion of disease-resistance gene families in trees.

Exploring the potential of gametic reconstruction of parental genotypes by F1 hybrids as a bridge for rapid introgression

Interspecific hybridization and genetic introgression are commonly observed in natural populations of many species, especially trees. Among oaks, gene flow between closely related species has been well documented. And yet, hybridization does not lead to a "melting pot", i.e., the homogenization of phenotypic traits. Here, we explore how the combination of several common reproductive and genomic traits could create an avenue for interspecific gene flow that partially explains this apparent paradox. During meiosis, F₁ hybrids will produce approximately (½)ⁿ "reconstructed" parental gametes, where n equals the number of chromosomes. Crossing over would introduce a small amount of introgressive material. The resulting parental-type gametophytes would probably possess a similar fertilization advantage as conspecific pollen. The resulting "backcross" would actually be the genetic equivalent of a conspecific out-cross, with a small amount of heterospecific DNA captured through crossing over. Even with detailed genomic analysis, the resulting offspring would not appear to be a backcross. This avenue for rapid introgression between species through the F₁ hybrid will be viable for organisms that meet certain conditions: low base chromosome number, conserved genomic structure and size, production of billions of gametes/gametophytes during each reproductive event, and conspecific fertilization advantage.

Decoding the oak genome: public release of sequence data, assembly, annotation and publication strategies

The 1.5 Gbp/2C genome of pedunculate oak (Quercus robur) has been sequenced. A strategy was established for dealing with the challenges imposed by the sequencing of such a large, complex and highly heterozygous genome by a whole-genome shotgun (WGS) approach, without the use of costly and time-consuming methods, such as fosmid or BAC clone-based hierarchical sequencing methods. The sequencing strategy combined short and long reads. Over 49 million reads provided by Roche 454 GS-FLX technology were assembled into contigs and combined with shorter Illumina sequence reads from paired-end and mate-pair libraries of different insert sizes, to build scaffolds. Errors were corrected and gaps filled with Illumina paired-end reads and contaminants detected, resulting in a total of 17,910 scaffolds (>2 kb) corresponding to 1.34 Gb. Fifty per cent of the assembly was accounted for by 1468 scaffolds (N50 of 260 kb). Initial comparison with the phylogenetically related Prunus persica gene model indicated that genes for 84.6% of the proteins present in peach (mean protein coverage of 90.5%) were present in our assembly. The second and third steps in this project are genome annotation and the assignment of scaffolds to the oak genetic linkage map. In accordance with the Bermuda and Fort Lauderdale agreements and the more recent Toronto Statement, the oak genome data have been released into public sequence repositories in advance of publication. In this presubmission paper, the oak genome consortium describes its principal lines of work and future directions for analyses of the nature, function and evolution of the oak genome.

Genetic mapping of natural variation in schooling tendency in the threespine stickleback

Although there is a heritable basis for many animal behaviors, the genetic architecture of behavioral variation in natural populations remains mostly unknown, particularly in vertebrates. We sought to identify the genetic basis for social affiliation in two populations of threespine sticklebacks (Gasterosteus aculeatus) that differ in their propensity to school. Marine sticklebacks from Japan school strongly whereas benthic sticklebacks from a lake in Canada are more solitary. Here, we expanded on our previous efforts to identify quantitative trait loci (QTL) for differences in schooling tendency. We tested fish multiple times in two assays that test different aspects of schooling tendency: 1) the model school assay, which presents fish with a school of eight model sticklebacks; and 2) the choice assay, in which fish are given a choice between the model school and a stationary artificial plant. We found low-to-moderate levels of repeatability, ranging from 0.1 to 0.5, in schooling phenotypes. To identify the genomic regions that contribute to differences in schooling tendency, we used QTL mapping in two types of crosses: benthic × marine backcrosses and an F2 intercross. We found two QTL for time spent with the school in the model school assay, and one QTL for number of approaches to the school in the choice assay. These QTL were on three different linkage groups, not previously linked to behavioral differences in sticklebacks. Our results highlight the importance of using multiple crosses and robust behavioral assays to uncover the genetic basis of behavioral variation in natural populations.

A linkage disequilibrium perspective on the genetic mosaic of speciation in two hybridizing Mediterranean white oaks

We analyzed the genetic mosaic of speciation in two hybridizing Mediterranean white oaks from the Iberian Peninsula (Quercus faginea Lamb. and Quercus pyrenaica Willd.). The two species show ecological divergence in flowering phenology, leaf morphology and composition, and in their basic or acidic soil preferences. Ninety expressed sequence tag-simple sequence repeats (EST-SSRs) and eight nuclear SSRs were genotyped in 96 trees from each species. Genotyping was designed in two steps. First, we used 69 markers evenly distributed over the 12 linkage groups (LGs) of the oak linkage map to confirm the species genetic identity of the sampled genotypes, and searched for differentiation outliers. Then, we genotyped 29 additional markers from the chromosome bins containing the outliers and repeated the multilocus scans. We found one or two additional outliers within four saturated bins, thus confirming that outliers are organized into clusters. Linkage disequilibrium (LD) was extensive; even for loosely linked and for independent markers. Consequently, score tests for association between two-marker haplotypes and the 'species trait' showed a broad genomic divergence, although substantial variation across the genome and within LGs was also observed. We discuss the influence of several confounding effects on neutrality tests and review the evolutionary processes leading to extensive LD. Finally, we examine how LD analyses within regions that contain outlier clusters and quantitative trait loci can help to identify regions of divergence and/or genomic hitchhiking in the light of predictions from ecological speciation theory.

Challenges and prospects in genome-wide quantitative trait loci mapping of standing genetic variation in natural populations

A considerable challenge in evolutionary genetics is to understand the genetic mechanisms that facilitate or impede evolutionary adaptation in natural populations. For this, we must understand the genetic loci contributing to trait variation and the selective forces acting on them. The decreased costs and increased feasibility of obtaining genotypic data on a large number of individuals have greatly facilitated gene mapping in natural populations, particularly because organisms whose genetics have been historically difficult to study are now within reach. Here we review the methods available to evolutionary ecologists interested in dissecting the genetic basis of traits in natural populations. Our focus lies on standing genetic variation in outbred populations. We present an overview of the current state of research in the field, covering studies on both plants and animals. We also draw attention to particular challenges associated with the discovery of quantitative trait loci and discuss parallels to studies on crops, livestock, and humans. Finally, we point to some likely future developments in genetic mapping studies.

Genetic determinism of phenological traits highly affected by climate change in Prunus avium: flowering date dissected into chilling and heat requirements

The present study investigated the genetic determinism of flowering date (FD), dissected into chilling (CR) and heat (HR) requirements. Elucidation of the genetic determinism of flowering traits is crucial to anticipate the increasing of ecological misalignment of adaptative traits with novel climate conditions in most temperate-fruit species. CR and HR were evaluated over 3 yr and FD over 5 yr in an intraspecific sweet cherry (Prunus avium) F1 progeny, and FD over 6 yr in a different F1 progeny. One quantitative trait locus (QTL) with major effect and high stability between years of evaluation was detected for CR and FD in the same region of linkage group (LG) 4. For HR, no stable QTL was detected. Candidate genes underlying the major QTL on LG4 were investigated and key genes were identified for CR and FD. Phenotypic dissection of FD and year repetitions allowed us to identify CR as the high heritable component of FD and a high genotype × environment interaction for HR. QTLs for CR reported in this study are the first described in this species. Our results provide a foundation for the identification of genes involved in CR and FD in sweet cherry which could be used to develop ideotypes adapted to future climatic conditions.

QTLs influencing carbohydrate and fat choice in a LOU/CxFischer 344 F2 rat population

OBJECTIVE

Individual differences in macronutrient selection, particularly fat and carbohydrate, and associated body weight gain are partly inherited as polygenic traits, but the potential genetic pathways are unknown. To give an overview of the Quantitative Trait Loci (QTLs) and candidate gene pathways influencing these differences in rat was aimed in this study.

DESIGN AND METHODS

To that end, F2 rats obtained from the crossbreeding between LOU/C and Fischer 344 rat strains to diet self-selection during 3 weeks were submitted. A genome scan was conducted with microsatellite markers covering evenly the whole genome. Genotypes and phenotypes were analyzed separately in male and female F2 rats by multiple interval mapping. Then, lists of candidate genes were treated by the Ingenuity Pathway software to propose gene pathways involved in our phenotypes.

RESULTS

Among numerous others, a QTL on chromosome 12 that influences body weight gain, and fat and carbohydrate choices in the LOU/C x Fischer 344 F2 rat population was found. This locus contains notably the acyl-co-A dehydrogenase gene.

CONCLUSION

A strong genetic determinism and complex pathways involving numerous candidate genes and processes, notably in accordance with the metabolic theory of feeding behavior control were found.

High rates of gene flow by pollen and seed in oak populations across Europe

Gene flow is a key factor in the evolution of species, influencing effective population size, hybridisation and local adaptation. We analysed local gene flow in eight stands of white oak (mostly Quercus petraea and Q. robur, but also Q. pubescens and Q. faginea) distributed across Europe.

Adult trees within a given area in each stand were exhaustively sampled (range [239, 754], mean 423), mapped, and acorns were collected ([17,147], 51) from several mother trees ([3], [47], 23). Seedlings ([65,387], 178) were harvested and geo-referenced in six of the eight stands. Genetic information was obtained from screening distinct molecular markers spread across the genome, genotyping each tree, acorn or seedling. All samples were thus genotyped at 5–8 nuclear microsatellite loci. Fathers/parents were assigned to acorns and seedlings using likelihood methods. Mating success of male and female parents, pollen and seed dispersal curves, and also hybridisation rates were estimated in each stand and compared on a continental scale.

On average, the percentage of the wind-borne pollen from outside the stand was 60%, with large variation among stands (21–88%). Mean seed immigration into the stand was 40%, a high value for oaks that are generally considered to have limited seed dispersal. However, this estimate varied greatly among stands (20–66%). Gene flow was mostly intraspecific, with large variation, as some trees and stands showed particularly high rates of hybridisation.

Our results show that mating success was unevenly distributed among trees. The high levels of gene flow suggest that geographically remote oak stands are unlikely to be genetically isolated, questioning the static definition of gene reserves and seed stands.

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.

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.

Comparative mapping in the Fagaceae and beyond with EST-SSRs

BACKGROUND

Genetic markers and linkage mapping are basic prerequisites for comparative genetic analyses, QTL detection and map-based cloning. A large number of mapping populations have been developed for oak, but few gene-based markers are available for constructing integrated genetic linkage maps and comparing gene order and QTL location across related species.

RESULTS

We developed a set of 573 expressed sequence tag-derived simple sequence repeats (EST-SSRs) and located 397 markers (EST-SSRs and genomic SSRs) on the 12 oak chromosomes (2n = 2x = 24) on the basis of Mendelian segregation patterns in 5 full-sib mapping pedigrees of two species: Quercus robur (pedunculate oak) and Quercus petraea (sessile oak). Consensus maps for the two species were constructed and aligned. They showed a high degree of macrosynteny between these two sympatric European oaks. We assessed the transferability of EST-SSRs to other Fagaceae genera and a subset of these markers was mapped in Castanea sativa, the European chestnut. Reasonably high levels of macrosynteny were observed between oak and chestnut. We also obtained diversity statistics for a subset of EST-SSRs, to support further population genetic analyses with gene-based markers. Finally, based on the orthologous relationships between the oak, Arabidopsis, grape, poplar, Medicago, and soybean genomes and the paralogous relationships between the 12 oak chromosomes, we propose an evolutionary scenario of the 12 oak chromosomes from the eudicot ancestral karyotype.

CONCLUSIONS

This study provides map locations for a large set of EST-SSRs in two oak species of recognized biological importance in natural ecosystems. This first step toward the construction of a gene-based linkage map will facilitate the assignment of future genome scaffolds to pseudo-chromosomes. This study also provides an indication of the potential utility of new gene-based markers for population genetics and comparative mapping within and beyond the Fagaceae.

Comparison of the genetic determinism of two key phenological traits, flowering and maturity dates, in three Prunus species: peach, apricot and sweet cherry

The present study investigates the genetic determinism of flowering and maturity dates, two traits highly affected by global climate change. Flowering and maturity dates were evaluated on five progenies from three Prunus species, peach, apricot and sweet cherry, during 3-8 years. Quantitative trait locus (QTL) detection was performed separately for each year and also by integrating data from all years together. High heritability estimates were obtained for flowering and maturity dates. Several QTLs for flowering and maturity dates were highly stable, detected each year of evaluation, suggesting that they were not affected by climatic variations. For flowering date, major QTLs were detected on linkage groups (LG) 4 for apricot and sweet cherry and on LG6 for peach. QTLs were identified on LG2, LG3, LG4 and LG7 for the three species. For maturity date, a major QTL was detected on LG4 in the three species. Using the peach genome sequence data, candidate genes underlying the major QTLs on LG4 and LG6 were investigated and key genes were identified. Our results provide a basis for the identification of genes involved in flowering and maturity dates that could be used to develop cultivar ideotypes adapted to future climatic conditions.

Assortative mating and gene flow generate clinal phenological variation in trees

BACKGROUND

On-going climate change is shifting the timing of bud burst (TBB) of broad leaf and conifer trees in temperate areas, raising concerns about the abilities of natural populations to respond to these shifts. The level of expected evolutionary change depends on the level and distribution of genetic variation of TBB. While numerous experimental studies have highlighted the role of divergent selection in promoting clinal TBB differentiation, we explored whether the observed patterns of variation could be generated by the joint effects of assortative mating for TBB and gene flow among natural populations. We tested this hypothesis using an in silico approach based on quantitative genetic models.

RESULTS

Our simulations showed that genetic clines can develop even without divergent selection. Assortative mating in association with environmental gradients substantially shifted the mean genetic values of populations. Owing to assortative mating, immigrant alleles were screened for proximal or distant populations depending on the strength of the environmental cline. Furthermore, we confirmed that assortative mating increases the additive genetic variance within populations. However, we observed also a rapid decline of the additive genetic variance caused by restricted gene flow between neighboring populations resulting from preferential matings between phenologically-matching phenotypes.

CONCLUSIONS

We provided evidence that the patterns of genetic variation of phenological traits observed in forest trees can be generated solely by the effects of assortative mating and gene flow. We anticipate that predicted temperature increases due to climate change will further enhance genetic differentiation across the landscape. These trends are likely to be reinforced or counteracted by natural selection if phenological traits are correlated to fitness.

What can patterns of differentiation across plant genomes tell us about adaptation and speciation?

Genome scans have become a common approach to identify genomic signatures of natural selection and reproductive isolation, as well as the genomic bases of ecologically relevant phenotypes, based on patterns of polymorphism and differentiation among populations or species. Here, we review the results of studies taking genome scan approaches in plants, consider the patterns of genomic differentiation documented and their possible causes, discuss the results in light of recent models of genomic differentiation during divergent adaptation and speciation, and consider assumptions and caveats in their interpretation. We find that genomic regions of high divergence generally appear quite small in comparisons of both closely and more distantly related populations, and for the most part, these differentiated regions are spread throughout the genome rather than strongly clustered. Thus, the genome scan approach appears well-suited for identifying genomic regions or even candidate genes that underlie adaptive divergence and/or reproductive barriers. We consider other methodologies that may be used in conjunction with genome scan approaches, and suggest further developments that would be valuable. These include broader use of sequence-based markers of known genomic location, greater attention to sampling strategies to make use of parallel environmental or phenotypic transitions, more integration with approaches such as quantitative trait loci mapping and measures of gene flow across the genome, and additional theoretical and simulation work on processes related to divergent adaptation and speciation.

Strength, diversity and plasticity of postmating reproductive barriers between two hybridizing oak species (Quercus robur L. and Quercus petraea (Matt) Liebl.)

Very little is known about the nature and strength of reproductive isolation (RI) in Quercus species, despite extensive research on the estimation and evolutionary significance of hybridization rates. We characterized postmating pre- and postzygotic RI between two hybridizing oak species, Quercus robur and Quercus petraea, using a large set of controlled crosses between different genotypes. Various traits potentially associated with reproductive barriers were quantified at several life history stages, from pollen-pistil interactions to seed set and progeny fitness-related traits. Results indicate strong intrinsic postmating prezygotic barriers, with significant barriers also at the postzygotic level, but relatively weaker extrinsic barriers on early hybrid fitness measures assessed in controlled conditions. Using general linear modelling of common garden data with clonal replicates, we showed that most traits exhibited important genotypic differences, as well as different levels of sensitivity to micro-environmental heterogeneity. These new findings suggest a large potential genetic diversity and plasticity of reproductive barriers and are confronted with hybridization evidence in these oak species.

A fast and cost-effective approach to develop and map EST-SSR markers: oak as a case study

Background

Expressed Sequence Tags (ESTs) are a source of simple sequence repeats (SSRs) that can be used to develop molecular markers for genetic studies. The availability of ESTs for Quercus robur and Quercus petraea provided a unique opportunity to develop microsatellite markers to accelerate research aimed at studying adaptation of these long-lived species to their environment. As a first step toward the construction of a SSR-based linkage map of oak for quantitative trait locus (QTL) mapping, we describe the mining and survey of EST-SSRs as well as a fast and cost-effective approach (bin mapping) to assign these markers to an approximate map position. We also compared the level of polymorphism between genomic and EST-derived SSRs and address the transferability of EST-SSRs in Castanea sativa (chestnut).

Results

A catalogue of 103,000 Sanger ESTs was assembled into 28,024 unigenes from which 18.6% presented one or more SSR motifs. More than 42% of these SSRs corresponded to trinucleotides. Primer pairs were designed for 748 putative unigenes. Overall 37.7% (283) were found to amplify a single polymorphic locus in a reference full-sib pedigree of Quercus robur. The usefulness of these loci for establishing a genetic map was assessed using a bin mapping approach. Bin maps were constructed for the male and female parental tree for which framework linkage maps based on AFLP markers were available. The bin set consisting of 14 highly informative offspring selected based on the number and position of crossover sites. The female and male maps comprised 44 and 37 bins, with an average bin length of 16.5 cM and 20.99 cM, respectively. A total of 256 EST-SSRs were assigned to bins and their map position was further validated by linkage mapping. EST-SSRs were found to be less polymorphic than genomic SSRs, but their transferability rate to chestnut, a phylogenetically related species to oak, was higher.

Conclusion

We have generated a bin map for oak comprising 256 EST-SSRs. This resource constitutes a first step toward the establishment of a gene-based map for this genus that will facilitate the dissection of QTLs affecting complex traits of ecological importance.

Genetic and genomic approaches to assess adaptive genetic variation in plants: forest trees as a model

With the increasing availability of sequence information at putatively important genes or regulatory regions, the characterization of adaptive genetic diversity and their association with phenotypic trait variation becomes feasible for many non-model organisms such as forest trees. Especially in predominantly outcrossing forest tree populations with large effective size, a high genetic variation in relevant genes is maintained, that is the raw material for the adaptation to changing and variable environments, and likewise for plant breeding. Oaks (Quercus spp.) are excellent model species to study the adaptation of forest trees to changing environments. They show a wide geographic distribution in Europe as dominant tree species in many forests and grow under a wide range of climatic and edaphic conditions. With the availability of a growing amount of functional and expressional candidate genes, we are now able to test the functional importance of single nucleotide polymorphisms (SNPs) by associating nucleotide variation in these genes with phenotypic variation in adaptive traits in segregating or natural populations. Here, we report on quantitative trait locus (QTL), candidate gene and association mapping approaches that are applicable to characterize gene markers and SNPs associated with variation in adaptive traits, such as bud burst, drought resistance and other traits showing selective responses to environmental change and stress. Because genome-wide association mapping studies are not feasible because of the enormous amount of SNP markers required in outcrossing trees with high recombination rates, the success of such an approach depends largely on the reasonable selection of candidate genes.

Utility of multilocus genotypes for taxon assignment in stands of closely related European white oaks from Switzerland

BACKGROUND AND AIMS

European white oaks (Quercus petraea, Q. pubescens, Q. robur) have long puzzled plant biologists owing to disputed species differentiation. Extensive hybridization or shared ancestry have been proposed as alternative hypotheses to explain why genetic differentiation between these oak species is low. Species delimitation is usually weak and often shows gradual transitions in leaf morphology. Hence, individual identification may be difficult, but remains a critical step for both scientific work and practical management.

METHODS

Multilocus genotype data (five nuclear microsatellites) were used from ten Swiss oak stands for taxon identification without a priori grouping of individuals or populations, using model-based Bayesian assignment tests.

KEY RESULTS

Three groups best structured the data, indicating that the taxonomical signal was stronger than the spatial signal. Most individuals showed high posterior probabilities for either of three genetic groups that were best circumscribed as taxonomical units. The assignment of a subset of trees, whose taxonomic status had been previously characterized in detail, supported this classification scheme.

CONCLUSIONS

Molecular-genetic assignment tests are useful in the identification of species status in critical taxon complexes such as the European white oaks. Such an approach is of practical importance for forest management, e.g. for stand certification or in seed trade to trace the origin of forest products.

QTL analysis of leaf morphological characters in a Quercus robur full-sib family (Q. robur x Q. robur ssp. slavonica)

The distinction between white oak species (section Quercus sensu stricto) is largely based on leaf morphological characters. There is, however, considerable within-species variation and no single species-diagnostic character, possibly due to phenotypic plasticity and/or underlying genetic variation. The aim of the present study was to identify quantitative trait loci (QTL) underlying the high within-species variation for leaf morphological characters in an F(1) full-sib family derived from a cross between Q. robur and Q. robur ssp. slavonica. In accordance with an earlier QTL mapping study in an intraspecific Q. robur full-sib family, polygenic inheritance was detected for leaf morphological characters that are used to discriminate between the species Quercus robur and Q. petraea. QTLs were distributed over ten linkage groups, showed a moderate effect in terms of phenotypic variance explained (PVE) in the mapping pedigree (3.6-9.6%), but accounted for a considerable amount of the parental differences. Co-localisation of QTLs on the same linkage group in different genetic backgrounds was found for the number and percentage of intercalary veins (NV, PV) on linkage group 3 and for NV on linkage group 5, revealing a high congruence in the relative QTL positions. The generally low correspondence of the other QTLs in the different mapping pedigrees may be an effect of the genetic background and of the environment. In conclusion, leaf morphological characters were found to be under polygenic control, and a comparison to earlier published results led to the identification of two QTLs that were stable across different genetic backgrounds.

A population genomic analysis of species boundaries: neutral processes, adaptive divergence and introgression between two hybridizing plant species

Little is known about the nature of species boundaries between closely related plant species and about the extent of introgression as a consequence of hybridization upon secondary contact. To address these topics we analyzed genome-wide differentiation between two closely related Silene species, Silene latifolia and S. dioica, and assessed the strength of introgression in sympatry. More than 300 AFLP markers were genotyped in three allopatric and three sympatric populations of each species. Outlier analyses were performed separately for sympatric and allopatric populations. Both positive and negative outlier loci were found, indicating that divergent and balancing selection, respectively, have shaped patterns of divergence between the two species. Sympatric populations of the two species were found to be less differentiated genetically than allopatric populations, indicating that hybridization has led to gene introgression. We conclude that differentiation between S. latifolia and S. dioica has been shaped by a combination of introgression and selection. These results challenge the view that species differentiation is a genome-wide phenomenon, and instead support the idea that genomes can be porous and that species differentiation has a genic basis.

Evidence for hybridization and introgression within a species-rich oak (Quercus spp.) community

Background

Analysis of interspecific gene flow is crucial for the understanding of speciation processes and maintenance of species integrity. Oaks (genus Quercus, Fagaceae) are among the model species for the study of hybridization. Natural co-occurrence of four closely related oak species is a very rare case in the temperate forests of Europe. We used both morphological characters and genetic markers to characterize hybridization in a natural community situated in west-central Romania and which consists of Quercus robur, Q. petraea, Q. pubescens, and Q. frainetto, respectively.

Results

On the basis of pubescence and leaf morphological characters ~94% of the sampled individuals were assigned to pure species. Only 16 (~6%) individual trees exhibited intermediate morphologies or a combination of characters of different species. Four chloroplast DNA haplotypes were identified in the study area. The distribution of haplotypes within the white oak complex showed substantial differences among species. However, the most common haplotypes were present in all four species. Furthermore, based on a set of 7 isozyme and 6 microsatellite markers and using a Bayesian admixture analysis without any a priori information on morphology we found that four genetic clusters best fit the data. There was a very good correspondence of each species with one of the inferred genetic clusters. The estimated introgression level varied markedly between pairs of species ranging from 1.7% between Q. robur and Q. frainetto to 16.2% between Q. pubescens and Q. frainetto. Only nine individuals (3.4%) appeared to be first-generation hybrids.

Conclusion

Our data indicate that natural hybridization has occurred at relatively low rates. The different levels of gene flow among species might be explained by differences in flowering time and spatial position within the stand. In addition, a partial congruence between phenotypically and genetically intermediate individuals was found, suggesting that intermediate appearance does not necessarily mean hybridization. However, it appears that natural hybridization did not seriously affect the species identity in this area of sympatry.

Influences of gene flow on adaptive speciation in theDubautia arborea-D. ciliolatacomplex

Mechanisms of reproductive isolation during plant speciation are often unclear because distinct species often experience high levels of gene flow and hybridization. Adaptive radiations such as the Hawaiian silversword alliance (HSA) provide unique opportunities to study the interactions of selection, gene flow and isolating mechanisms during the speciation process. We examined patterns of phenotypic and genetic differentiation in Dubautia arborea and Dubautia ciliolata, two parapatric HSA taxa that show marked morphological divergence but evidence of weak molecular differentiation, in order to estimate genome-wide differentiation and gene flow patterns. We scored 166 amplified fragment length polymorphism markers in a set of 89 plants from two populations each of D. arborea and D. ciliolata and phenotypically D. arborea-like and D. ciliolata-like plants from a natural hybrid zone. Analyses of population subdivision showed low levels of differentiation between the two species (F(ST) = 0.089) and evidence that the phenotypically parental hybrid zone plants were largely of parental species rather than of hybrid origin. A Bayesian analysis of population ancestry identified a number of plants with admixed D. arborea and D. ciliolata ancestry, even in nonhybrid-zone populations. These results suggest that genome-wide low levels of differentiation between D. arborea and D. ciliolata are in part due to gene flow, and favour models of genic speciation and collective evolution in which gene flow has different effects on selected loci vs. nonselected genomic regions. We discuss ecological and climatic factors that may have shaped patterns of differentiation in this species complex.

Linkage maps of the dwarf and Normal lake whitefish (Coregonus clupeaformis) species complex and their hybrids reveal the genetic architecture of population divergence

Elucidating the genetic architecture of population divergence may reveal the evolution of reproductive barriers and the genomic regions implicated in the process. We assembled genetic linkage maps for the dwarf and Normal lake whitefish species complex and their hybrids. A total of 877 AFLP loci and 30 microsatellites were positioned. The homology of mapped loci between families supported the existence of 34 linkage groups (of 40n expected) exhibiting 83% colinearity among linked loci between these two families. Classes of AFLP markers were not randomly distributed among linkage groups. Both AFLP and microsatellites exhibited deviations from Mendelian expectations, with 30.4% exhibiting significant segregation distortion across 28 linkage groups of the four linkage maps in both families (P < 0.00001). Eight loci distributed over seven homologous linkage groups were significantly distorted in both families and the level of distortion, when comparing homologous loci of the same phase between families, was correlated (Spearman R = 0.378, P = 0.0021). These results suggest that substantial divergence incurred during allopatric glacial separation and subsequent sympatric ecological specialization has resulted in several genomic regions that are no longer complementary between dwarf and Normal populations issued from different evolutionary glacial lineages.

QTL analysis of intraspecific differences between two Silene vulgaris ecotypes

BACKGROUND AND AIMS

Serpentine soils provide a highly selective substrate for plant colonization and growth and represent an ideal system for studying the evolution of plant-ecotypes. In the present study the aim was to identify the genetic architecture of morphological traits distinguishing serpentine and non-serpentine ecotypes of Silene vulgaris.

METHODS

Using an F(2) mapping population derived from an intraspecific cross between a serpentine and a non-serpentine ecotype of S. vulgaris, the genetic architecture of 12 morphological traits was explored using a quantitative trait locus (QTL) analysis.

KEY RESULTS

The QTL analysis identified a total of 49 QTLs, of which 24 were classified as major QTLs. The mean number of QTLs per trait category was found to correspond well with numbers reported in the literature for similar crosses. Clustering of QTLs for different traits was found on several linkage groups.

CONCLUSIONS

Morphological traits that differentiate the two ecotypes are strongly correlated, presumably as a consequence of the joint effects of extensive linkage of QTLs for different traits and directional selection. The signature of consistent directional selection was found for leaf and shoot trait divergence. Intraspecific ecotype differences in S. vulgaris were found to be distributed across the entire genome. The study shows that QTL analyses on non-model organisms can provide novel insights into the genetic basis of plant diversification.

Barrier to gene flow between two ecologically divergent Populus species, P. alba (white poplar) and P. tremula (European aspen): the role of ecology and life history in gene introgression

The renewed interest in the use of hybrid zones for studying speciation calls for the identification and study of hybrid zones across a wide range of organisms, especially in long-lived taxa for which it is often difficult to generate interpopulation variation through controlled crosses. Here, we report on the extent and direction of introgression between two members of the "model tree" genus Populus: Populus alba (white poplar) and Populus tremula (European aspen), across a large zone of sympatry located in the Danube valley. We genotyped 93 hybrid morphotypes and samples from four parental reference populations from within and outside the zone of sympatry for a genome-wide set of 20 nuclear microsatellites and eight plastid DNA restriction site polymorphisms. Our results indicate that introgression occurs preferentially from P. tremula to P. alba via P. tremula pollen. This unidirectional pattern is facilitated by high levels of pollen vs. seed dispersal in P. tremula (pollen/seed flow = 23.9) and by great ecological opportunity in the lowland floodplain forest in proximity to P. alba seed parents, which maintains gene flow in the direction of P. alba despite smaller effective population sizes (N(e)) in this species (P. alba N(e)c. 500-550; P. tremula N(e)c. 550-700). Our results indicate that hybrid zones will be valuable tools for studying the genetic architecture of the barrier to gene flow between these two ecologically divergent Populus species.

Results on quantitative trait loci for flushing date in oaks can be transferred to different segregating progenies

Flushing date (bud burst) is one of the most important traits for the adaptation to different environments and climates in the temperate zone. Because of their wide geographic distribution, Quercus robur L. and Q. petraea (Matt.) Liebl. are suitable as model plants to study the genetic basis of bud burst. QTLs (Quantitative Trait Loci) with comparatively large effects have been mapped in a former study in a Q. robur x Q. robur full-sib family (French cross). In the present study, we performed a Bulked Segregant Analysis (BSA) in the F (1) progeny comprising 144 seedlings derived from a cross between a single Q. robur tree as common seed parent and five different pollen donors both from Q. robur and Q. petraea (Q. robur x Q. spp., Diekholzen crosses). In addition, markers linked to two bud burst QTLs with comparably strong effect in the above-mentioned full-sib family (French cross) were tested for their association with bud burst in the Q. robur x Q. spp. (Diekholzen) progeny. Using three microsatellite markers as anchor points, we could map QTLs on linkage group 7 and on linkage group 2, together explaining 16.2 % of the total phenotypic variance (PVE) in 1999 and 38.1 % in 2003. Out of 10 markers that segregated in both mapping progenies, four markers including the two microsatellite markers, showed a significant effect on bud burst in both materials. At microsatellite loci ssrQpZAG1/5 (linkage group 7) and ssrQpZAG119 (linkage group 2) alleles associated with early (allele 166 bp in ssrQpZAG1/5) and late bud burst (allele 57 bp in ssrQpZAG119) in the Q. robur x Q. robur full-sib family (French cross) showed a highly significant association with the same polarity of the effect in the Q. robur x Q. spp. (Diekholzen) progeny. The usefulness of these markers for marker-assisted selection in full-sib and half-sib families is discussed.

Genome scanning for interspecific differentiation between two closely related oak species [Quercus robur L. and Q. petraea (Matt.) Liebl.]

Interspecific differentiation values (G(ST)) between two closely related oak species (Quercus petraea and Q. robur) were compiled across different studies with the aim to explore the distribution of differentiation at the genome level. The study was based on a total set of 389 markers (isozymes, AFLPs, SCARs, microsatellites, and SNPs) for which allelic frequencies were estimated in pairs of populations sampled throughout the sympatric distribution of the two species. The overall distribution of G(ST) values followed an L-shaped curve with most markers exhibiting low species differentiation (G(ST) < 0.01) and only a few loci reaching >10% levels. Twelve percent of the loci exhibited significant G(ST) deviations to neutral expectations, suggesting that selection contributed to species divergence. Coding regions expressed higher differentiation than noncoding regions. Among the 389 markers, 158 could be mapped on the 12 linkage groups of the existing Q. robur genetic map. Outlier loci with large G(ST) values were distributed over 9 linkage groups. One cluster of three outlier loci was found within 0.51 cM; but significant autocorrelation of G(ST) was observed at distances <2 cM. The size and distribution of genomic regions involved in species divergence are discussed in reference to hitchhiking effects and disruptive selection.

Evidence for shared ancestral polymorphism rather than recurrent gene flow at microsatellite loci differentiating two hybridizing oaks (Quercus spp.)

Quercus petraea and Quercus robur are two closely related oak species, considered to hybridize. Genetic markers, however, indicate that despite sharing most alleles, the two species remain separate genetic units. Analysis of 20 microsatellite loci in multiple populations from both species suggested a genome-wide differentiation. Thus, the allele sharing between both species could be explained either by low rates of gene flow or shared ancestral variation. We performed further analyses of population differentiation in a biogeographical setting and an admixture analysis in mixed oak stands to distinguish between both hypotheses. Based on our results we propose that the low genetic differentiation among these species results from shared ancestry rather than high rates of gene flow.

Detection of quantitative trait loci controlling bud burst and height growth in Quercus robur L

Genetic variation of bud burst and early growth components was estimated in a full-sib family of Quercus robur L. comprising 278 offspring. The full sibs were vegetatively propagated, and phenotypic assessments were made in three field tests. This two-generation pedigree was also used to construct a genetic linkage map (12 linkage groups, 128 markers) and locate quantitative trait loci (QTLs) controlling bud burst and growth components. In each field test, the date of bud burst extended over a period of 20 days from the earliest to the latest clone. Bud burst exhibited higher heritability (0.15-0.51) than growth components (0.04-0.23) and also higher correlations across field tests. Over the three tests there were 32 independent detected QTLs ( P<or=5% at the chromosome level) controlling bud burst, which likely represent at least 12 unique genes or chromosomal regions controlling this trait. QTLs explained from 3% to 11% of the variance of the clonal means. The number of QTLs controlling height growth components was lower and varied between two and four. However the contribution of each QTL to the variance of the clonal mean was higher (from 4% to 19%). These results indicate that the genetic architecture of two important fitness-related traits are quite different. On the one hand, bud burst is controlled by several QTLs with rather low to moderate effects, but contributing to a high genetic (additive) variance. On the other hand, height growth depends on fewer QTLs with moderate to strong effects, resulting in lower heritabilities of the trait.

Hybrid zones as a tool for identifying adaptive genetic variation in outbreeding forest trees: lessons from wild annual sunflowers (Helianthus spp.)

The identification and study of adaptively important genes in forest trees represents a formidable challenge because of their long generation spans. In annual or perennial herbs, formal genetic studies can be employed to identify the quantitative trait loci (QTLs) and/or candidate genes that underlie important traits, and the segregating populations can be transplanted into natural populations to measure the strength and direction of selection. However, the application of these methods to forest trees is difficult, because the creation of appropriate genetic material is extremely time-consuming in long-lived, woody plants, and lifetime fitness estimates are difficult or impossible to obtain. Although QTL mapping should in principle be feasible in wild intraspecific populations (as an alternative to artificial crosses), this approach is less likely to be successful in trees because LD (linkage disequilibrium) will decay quickly in large outbreeding plant populations. Within the present paper, we discuss a modified approach based on natural hybrid zones. We describe the use of wild annual sunflowers (Helianthus spp.) as a model for exploring the hybrid zone approach. Transplanted experimental hybrids allowed us to assess the adaptive value of individual chromosomal blocks in nature, and data on natural Helianthus hybrids suggest that similar approaches are possible in natural hybrid zones. Our results allowed us to test the role of hybridization in the origin of ecological divergence in wild sunflowers. In addition, they have practical implications for identifying adaptively important genes or QTLs in trees. This is exemplified by three temperate forest taxa, Populus (poplars, aspens, cottonwoods), Fraxinus (ash), and Quercus (oak). All three are diploid and important genomic tools are under development. Moreover, all three offer extensive hybrid zones whose likely age can be inferred from fossil data. Age data enables estimates of the size and frequency of chromosomal blocks in hybrids, thereby providing guidance in designing marker-based experiments. We predict that natural hybrid zones will be valuable tools for identifying the QTLs and/or candidate genes responsible for adaptive traits in forest trees.