Forest-tree population genomics and adaptive evolution

Local adaptation of Pinus leiophylla under climate and land use change models in the Avocado Belt of Michoacán

Climate change and land use change are two main drivers of global biodiversity decline, decreasing the genetic diversity that populations harbour and altering patterns of local adaptation. Landscape genomics allows measuring the effect of these anthropogenic disturbances on the adaptation of populations. However, both factors have rarely been considered simultaneously. Based on a set of 3660 SNPs from which 130 were identified as outliers by a genome-environment association analysis (LFMM), we modelled the spatial turnover of allele frequencies in 19 localities of Pinus leiophylla across the Avocado Belt in Michoacán state, Mexico. Then, we evaluated the effect of climate change and land use change scenarios, in addition to evaluating assisted gene flow strategies and connectivity metrics across the landscape to identify priority conservation areas for the species. We found that localities in the centre-east of the Avocado Belt would be more vulnerable to climate change, while localities in the western area are more threatened by land conversion to avocado orchards. Assisted gene flow actions could aid in mitigating both threats. Connectivity patterns among forest patches will also be modified by future habitat loss, with central and eastern parts of the Avocado Belt maintaining the highest connectivity. These results suggest that areas with the highest priority for conservation are in the eastern part of the Avocado Belt, including the Monarch Butterfly Biosphere Reserve. This work is useful as a framework that incorporates distinct layers of information to provide a more robust representation of the response of tree populations to anthropogenic disturbances.

Signatures of local adaptation to current and future climate in phenology-related genes in natural populations of Quercus robur

BACKGROUND

Local adaptation is a key evolutionary process that enhances the growth of plants in their native habitat compared to non-native habitats, resulting in patterns of adaptive genetic variation across the entire geographic range of the species. The study of population adaptation to local environments and predicting their response to future climate change is important because of climate change.

RESULTS

Here, we explored the genetic diversity of candidate genes associated with bud burst in pedunculate oak individuals sampled from 6 populations in Poland. Single nucleotide polymorphism (SNP) diversity was assessed in 720 candidate genes using the sequence capture technique, yielding 18,799 SNPs. Using landscape genomic approaches, we identified 8 FST outliers and 781 unique SNPs in 389 genes associated with geography, climate, and phenotypic variables (individual/family spring and autumn phenology, family diameter at breast height (DBH), height, and survival) that are potentially involved in local adaptation. Then, using a nonlinear multivariate model, Gradient Forests, we identified vulnerable areas of the pedunculate oak distribution in Poland that are at risk from climate change.

CONCLUSIONS

The model revealed that pedunculate oak populations in the eastern part of the analyzed geographical region are the most sensitive to climate change. Our results might offer an initial evaluation of a potential management strategy for preserving the genetic diversity of pedunculate oak.

Strategies to mitigate shifts in red oak (Quercus sect. Lobatae) distribution under a changing climate

Red oaks (Quercus sect. Lobatae) are a taxonomic group of hardwood trees, which occur in swamp forests, subtropical chaparral and savannahs from Columbia to Canada. They cover a wide range of ecological niches, and many species are thought to be able to cope with current trends in climate change. Genus Quercus encompasses ca. 500 species, of which ca. 80 make up sect. Lobatae. Species diversity is greatest within the southeastern USA and within the northern and eastern regions of Mexico. This review discusses the weak reproductive barriers between species of red oaks and the effects this has on speciation and niche range. Distribution and diversity have been shaped by drought adaptations common to the species of sect. Lobatae, which enable them to fill various xeric niches across the continent. Drought adaptive traits of this taxonomic group include deciduousness, deep tap roots, ring-porous xylem, regenerative stump sprouting, greater leaf thickness and smaller stomata. The complex interplay between these anatomical and morphological traits has given red oaks features of drought tolerance and avoidance. Here, we discuss physiological and genetic components of these adaptations to address how many species of sect. Lobatae reside within xeric sites and/or sustain normal metabolic function during drought. Although extensive drought adaptation appears to give sect. Lobatae a resilience to climate change, aging tree stands, oak life history traits and the current genetic structures place many red oak species at risk. Furthermore, oak decline, a complex interaction between abiotic and biotic agents, has severe effects on red oaks and is likely to accelerate species decline and fragmentation. We suggest that assisted migration can be used to avoid species fragmentation and increase climate change resilience of sect. Lobatae.

Evaluation of environmental factors affecting the genetic diversity, genetic structure, and the potential distribution of Rhododendron aureum Georgi under changing climate

Understanding genetic variation and structure, adaptive genetic variation, and its relationship with environmental factors is of great significance to understand how plants adapt to climate change and design effective conservation and management strategies. The objective of this study was to (I) investigate the genetic diversity and structure by AFLP markers in 36 populations of R. aureum from northeast China, (Ⅱ) reveal the relative contribution of geographical and environmental impacts on the distribution and genetic differentiation of R. aureum, (Ⅲ) identify outlier loci under selection and evaluate the association between outlier loci and environmental factors, and (Ⅳ) exactly calculate the development trend of population of R. aureum, as it is confronted with severe climate change and to provide information for designing effective conservation and management strategies. We found high genetic variation (I = 0.584) and differentiation among populations (ΦST  = 0.703) and moderate levels of genetic diversity within populations of R. aureum. A significant relationship between genetic distance and environmental distance was identified, which suggested that the differentiation of different populations was caused by environmental factors. Using BayeScan and Dfdist, 42 outlier loci are identified and most of the outlier loci are associated with climate or relief factors, suggesting that these loci are linked to genes that are involved in the adaptability of R. aureum to the environment. Species distribution models (SDMs) showed that climate warming will cause a significant reduction in suitable areas for R. aureum, especially under the RCP 85 scenario. Our results help to understand the potential response of R. aureum to climatic changes and provide new perspectives for R. aureum resource management and conservation strategies.

Genetic Architecture and Genome-Wide Adaptive Signatures Underlying Stem Lenticel Traits in Populus tomentosa

The stem lenticel is a highly specialized tissue of woody plants that has evolved to balance stem water retention and gas exchange as an adaptation to local environments. In this study, we applied genome-wide association studies and selective sweeping analysis to characterize the genetic architecture and genome-wide adaptive signatures underlying stem lenticel traits among 303 unrelated accessions of P. tomentosa, which has significant phenotypic and genetic variations according to climate region across its natural distribution. In total, we detected 108 significant single-nucleotide polymorphisms, annotated to 88 candidate genes for lenticel, of which 9 causative genes showed significantly different selection signatures among climate regions. Furthermore, PtoNAC083 and PtoMYB46 showed significant association signals and abiotic stress response, so we overexpressed these two genes in Arabidopsis thaliana and found that the number of stem cells in all three overexpression lines was significantly reduced by PtoNAC083 overexpression but slightly increased by PtoMYB46 overexpression, suggesting that both genes are involved in cell division and expansion during lenticel formation. The findings of this study demonstrate the successful application of an integrated strategy for dissecting the genetic basis and landscape genetics of complex adaptive traits, which will facilitate the molecular design of tree ideotypes that may adapt to future climate and environmental changes.

Multi-Approach Analysis Reveals Local Adaptation in a Widespread Forest Tree of Reunion Island

Detecting processes of local adaptation in forest trees and identifying environmental selective drivers are of primary importance for forest management and conservation. Transplant experiments, functional genomics and population genomics are complementary tools to efficiently characterize heritable phenotypic traits and to decipher the genetic bases of adaptive traits. Using an integrative approach combining phenotypic assessment in common garden, transcriptomics and landscape genomics, we investigated leaf adaptive traits in Coffea mauritiana, a forest tree endemic to Reunion Island. Eight populations of C. mauritiana originating from sites with contrasted environmental conditions were sampled in common garden to assess several leaf morphological traits, to analyze the leaf transcriptome and leaf cuticular wax composition. The relative alkane content of cuticular waxes was significantly correlated with major climatic gradients, paving the way for further transcriptome-based analyses. The expression pattern of cuticle biosynthetic genes was consistent with a modulation of alkane accumulation across the population studied, supporting the hypothesis that the composition of cuticular wax is involved in the local adaptation of C. mauritiana. Association tests in landscape genomics performed using RNA-seq-derived single-nucleotide polymorphisms revealed that genes associated with cell wall remodeling also likely play an adaptive role. By combining these different approaches, this study efficiently identified local adaptation processes in a non-model species. Our results provide the first evidence for local adaptation in trees endemic to Reunion Island and highlight the importance of cuticle composition for the adaptation of trees to the high evaporative demand in warm climates.

The GenTree Platform: growth traits and tree-level environmental data in 12 European forest tree species

BACKGROUND

Progress in the field of evolutionary forest ecology has been hampered by the huge challenge of phenotyping trees across their ranges in their natural environments, and the limitation in high-resolution environmental information.

FINDINGS

The GenTree Platform contains phenotypic and environmental data from 4,959 trees from 12 ecologically and economically important European forest tree species: Abies alba Mill. (silver fir), Betula pendula Roth. (silver birch), Fagus sylvatica L. (European beech), Picea abies (L.) H. Karst (Norway spruce), Pinus cembra L. (Swiss stone pine), Pinus halepensis Mill. (Aleppo pine), Pinus nigra Arnold (European black pine), Pinus pinaster Aiton (maritime pine), Pinus sylvestris L. (Scots pine), Populus nigra L. (European black poplar), Taxus baccata L. (English yew), and Quercus petraea (Matt.) Liebl. (sessile oak). Phenotypic (height, diameter at breast height, crown size, bark thickness, biomass, straightness, forking, branch angle, fructification), regeneration, environmental in situ measurements (soil depth, vegetation cover, competition indices), and environmental modeling data extracted by using bilinear interpolation accounting for surrounding conditions of each tree (precipitation, temperature, insolation, drought indices) were obtained from trees in 194 sites covering the species' geographic ranges and reflecting local environmental gradients.

CONCLUSION

The GenTree Platform is a new resource for investigating ecological and evolutionary processes in forest trees. The coherent phenotyping and environmental characterization across 12 species in their European ranges allow for a wide range of analyses from forest ecologists, conservationists, and macro-ecologists. Also, the data here presented can be linked to the GenTree Dendroecological collection, the GenTree Leaf Trait collection, and the GenTree Genomic collection presented elsewhere, which together build the largest evolutionary forest ecology data collection available.

Resin ducts as resistance traits in conifers: linking dendrochronology and resin-based defences

Conifers have evolved different chemical and anatomical defences against a wide range of antagonists. Resin ducts produce, store and translocate oleoresin, a complex terpenoid mixture that acts as both a physical and a chemical defence. Although resin duct characteristics (e.g., number, density, area) have been positively related to biotic resistance in several conifer species, the literature reporting this association remains inconclusive. Axial resin ducts recorded in annual growth rings are an archive of annual defensive investment in trees. This whole-life record of defence investment can be analysed using standard dendrochronological procedures, which allows us to assess interannual variability and the effect of understudied drivers of phenotypic variation on resin-based defences. Understanding the sources of phenotypic variation in defences, such as genetic differentiation and environmental plasticity, is essential for assessing the adaptive potential of forest tree populations to resist pests under climate change. Here, we reviewed the evidence supporting the importance of resin ducts in conifer resistance, and summarized current knowledge about the sources of variation in resin duct production. We propose a standardized methodology to measure resin duct production by means of dendrochronological procedures. This approach will illuminate the roles of resin ducts in tree defence across species, while helping to fill pivotal knowledge gaps in plant defence theory, and leading to a robust understanding of the patterns of variation in resin-based defences throughout the tree's lifespan.

High Photosynthetic Rates in a Solanum pennellii Chromosome 2 QTL Is Explained by Biochemical and Photochemical Changes

Enhanced photosynthesis is strictly associated with to productivity and it can be accomplished by genetic approaches through identification of genetic variation. By using a Solanum pennellii introgression lines (ILs) population, it was previously verified that, under normal (CO2), IL 2-5 and 2-6 display increased photosynthetic rates by up to 20% in comparison with their parental background (M82). However, the physiological mechanisms involved in the enhanced CO2 assimilation exhibited by these lines remained unknown, precluding their use for further biotechnological applications. Thereby, here we attempted to uncover the physiological factors involved in the upregulation of photosynthesis in ILs 2-5 and 2-6 under normal (CO2) as well as under elevated (CO2). The results provide evidence for increased biochemical capacity (higher maximum carboxylation velocity and maximum electron transport rate) in plants from IL 2-5 and 2-6, whereas the diffusive components (stomatal and mesophyll conductances) were unaltered in these ILs in comparison to M82. Our analyses revealed that the higher photosynthetic rate observed in these ILs was associated with higher levels of starch as well as total protein levels, specially increased RuBisCO content. Further analyses performed in plants under high (CO2) confirmed that biochemical properties are involved in genetic variation on chromosome 2 related to enhanced photosynthesis.

Testing the Effect of Mountain Ranges as a Physical Barrier to Current Gene Flow and Environmentally Dependent Adaptive Divergence in Cunninghamia konishii (Cupressaceae)

Populations can be genetically isolated by differences in their ecology or environment that hampered efficient migration, or they may be isolated solely by geographic distance. Moreover, mountain ranges across a species’ distribution area might have acted as barriers to gene flow. Genetic variation was quantified using amplified fragment length polymorphism (AFLP) and 13 selective amplification primer combinations used generated a total of 482 fragments. Here, we tested the barrier effects of mountains on gene flow and environmentally dependent local adaptation of Cunninghamia konishii occur in Taiwan. A pattern of genetic isolation by distance was not found and variation partitioning revealed that environment explained a relatively larger proportion of genetic variation than geography. The effect of mountains as barriers to genetic exchange, despite low population differentiation indicating a high rate of gene flow, was found within the distribution range of C. konishii. Twelve AFLP loci were identified as potential selective outliers using genome-scan methods (BAYESCAN and DFDIST) and strongly associated with environmental variables using regression approaches (LFMM, Samβada, and rstanarm) demonstrating adaptive divergence underlying local adaptation. Annual mean temperature, annual precipitation, and slope could be the most important environmental factors causally associated with adaptive genetic variation in C. konishii. The study revealed the existence of physical barriers to current gene flow and environmentally dependent adaptive divergence, and a significant proportion of the rate of gene flow may represent a reflection of demographic history.

Multitrait genome-wide association analysis of Populus trichocarpa identifies key polymorphisms controlling morphological and physiological traits

Genome-wide association studies (GWAS) have great promise for identifying the loci that contribute to adaptive variation, but the complex genetic architecture of many quantitative traits presents a substantial challenge. We measured 14 morphological and physiological traits and identified single nucleotide polymorphism (SNP)-phenotype associations in a Populus trichocarpa population distributed from California, USA to British Columbia, Canada. We used whole-genome resequencing data of 882 trees with more than 6.78 million SNPs, coupled with multitrait association to detect polymorphisms with potentially pleiotropic effects. Candidate genes were validated with functional data. Broad-sense heritability (H² ) ranged from 0.30 to 0.56 for morphological traits and 0.08 to 0.36 for physiological traits. In total, 4 and 20 gene models were detected using the single-trait and multitrait association methods, respectively. Several of these associations were corroborated by additional lines of evidence, including co-expression networks, metabolite analyses, and direct confirmation of gene function through RNAi. Multitrait association identified many more significant associations than single-trait association, potentially revealing pleiotropic effects of individual genes. This approach can be particularly useful for challenging physiological traits such as water-use efficiency or complex traits such as leaf morphology, for which we were able to identify credible candidate genes by combining multitrait association with gene co-expression and co-methylation data.

Neutral and Climate-Driven Adaptive Processes Contribute to Explain Population Variation in Resin Duct Traits in a Mediterranean Pine Species

Resin ducts are important anatomical defensive traits related to biotic resistance in conifers. Previous studies have reported intraspecific genetic variation in resin duct characteristics. However, little is currently known about the micro-evolutionary patterns and adaptive value of these defensive structures. Here, we quantified inter-population genetic variation in resin duct features and their inducibility in Pinus pinaster and assessed whether such variation was associated with climate gradients. To that end, we characterized the resin duct system of 2-year-old saplings from 10 populations across the species' distribution range. We measured axial resin duct features (density, mean size, and percentage conductive area of resin ducts) and their inducibility in response to methyl jasmonate. Genotyping of single nucleotide polymorphisms allowed to account for the population genetic structure in our models in order to avoid spurious correlations between resin duct characteristics and climate. We found large inter-population variation in resin duct density and conductive area, but not in their inducibility. Our results suggest that population variation in the percentage conductive area of resin ducts likely arise from adaptation to local climate conditions. This study highlights the adaptive relevance of resin ducts and helps to shed light on the micro-evolutionary patterns of resin-based defenses in conifers.

Unexpectedly high genetic diversity and divergence among populations of the apomictic Neotropical tree Miconia albicans

Since tropical trees often have long generation times and relatively small reproductive populations, breeding systems and genetic variation are important for population viability and have consequences for conservation. Miconia albicans is an obligate, diplosporous, apomictic species widespread in the Brazilian Cerrado, the savanna areas in central Brazil and elsewhere in the Neotropics. The genetic variability would be, theoretically, low within these male-sterile and possibly clonal populations, although some variation would be expected due to recombination during restitutional meiosis. We used ISSR markers to assess genetic diversity of M. albicans and to compare with other tropical trees, including invasive species of Melastomataceae. A total of 120 individuals from six populations were analysed using ten ISSR primers, which produced 153 fully reproducible fragments. The populations of M. albicans presented mean Shannon's information index (I) of 0.244 and expected heterozygosity (H_e ) of 0.168. Only two pairs of apparently clonal trees were identified, and genetic diversity was relatively high. A hierarchical amova for all ISSR datasets showed that 74% of the variance was found among populations, while only 26% of the variance was found within populations of this species. Multivariate and Bayesian analyses indicated marked separation between the studied populations. The genetic diversity generated by restitutional meiosis, polyploidy and possibly other genome changes may explain the morpho-physiological plasticity and the ability of these plants to differentiate and occupy such a wide territory and different environmental conditions. Producing enormous amounts of bird-dispersed fruits, M. albicans possess weedy potential that may rival other Melastomataceae alien invaders.

The genetics of drought tolerance in conifers

Contents 1034 I. 1034 II. 1035 III. 1037 IV. 1038 V. 1042 VI. 1043 VII. 1045 References 1045 SUMMARY: As temperatures warm and precipitation patterns shift as a result of climate change, interest in the identification of tree genotypes that will thrive under more arid conditions has grown. In this review, we discuss the multiple definitions of 'drought tolerance' and the biological processes involved in drought responses. We describe the three major approaches taken in the study of genetic variation in drought responses, the advantages and shortcomings of each, and what each of these approaches has revealed about the genetic basis of adaptation to drought in conifers. Finally, we discuss how a greater knowledge of the genetics of drought tolerance may aid forest management, and provide recommendations for how future studies may overcome the limitations of past approaches. In particular, we urge a more direct focus on survival, growth and the traits that directly predict them (rather than on proxies, such as water use efficiency), combining research approaches with complementary strengths and weaknesses, and the inclusion of a wider range of taxa and life stages.

Influence of altitude on local adaptation in upland tree species from central Argentina

Steep climatic gradients boost morphological and physiological adjustments in plants, with consequences on performance. The three principal woody species of the Sierras Grandes Mountains of central Argentina have marked differences in sapling performance along their altitudinal distribution. We hypothesize that the steep gradient of climatic conditions across the species' altitudinal distribution promotes trait differences between populations of different altitudes that are inherited by the following generation. Seeds from different altitudes were exposed to three temperature regimes to assess differential germination responses. Saplings were then transplanted to a greenhouse to assess possible variations in attributes and performance after 18 months. The three species showed differences in germination responses to temperature among altitudes and/or in sapling attributes and performance. In Maytenus boaria and Escallonia cordobensis, germination success was higher under high temperatures for the highest-altitude, whereas lower temperatures boosted germination of the lowest altitudes. Polylepis australis showed no differences in germination among temperature treatments. In the greenhouse, saplings of the three species from intermediate altitudes showed high performance, whereas the upper and lower populations seemed to be adjusted to tolerating more stressful conditions (i.e., lower temperatures at the upper end and water stress at the lower end), showing lower performance toward both altitudinal limits. These patterns agree with those described for saplings growing under field conditions, suggesting adjustments in response to environmental changes undergone by populations along the altitudinal range. The marked adjustments of populations to the local environment suggest a potentially high impact of climatic change on species distribution.

Forest genetic monitoring: an overview of concepts and definitions

Safeguarding sustainability of forest ecosystems with their habitat variability and all their functions is of highest priority. Therefore, the long-term adaptability of forest ecosystems to a changing environment must be secured, e.g., through sustainable forest management. High adaptability is based on biological variation starting at the genetic level. Thus, the ultimate goal of the Convention on Biological Diversity (CBD) to halt the ongoing erosion of biological variation is of utmost importance for forest ecosystem functioning and sustainability. Monitoring of biological diversity over time is needed to detect changes that threaten these biological resources. Genetic variation, as an integral part of biological diversity, needs special attention, and its monitoring can ensure its effective conservation. We compare forest genetic monitoring to other biodiversity monitoring concepts. Forest genetic monitoring (FGM) enables early detection of potentially harmful changes of forest adaptability before these appear at higher biodiversity levels (e.g., species or ecosystem diversity) and can improve the sustainability of applied forest management practices and direct further research. Theoretical genetic monitoring concepts developed up to now need to be evaluated before being implemented on a national and international scale. This article provides an overview of FGM concepts and definitions, discusses their advantages and disadvantages, and provides a flow chart of the steps needed for the optimization and implementation of FGM. FGM is an important module of biodiversity monitoring, and we define an effective FGM scheme as consisting of an assessment of a forest population's capacity to survive, reproduce, and persist under rapid environmental changes on a long-term scale.

AFLP Genome Scanning Reveals Divergent Selection in Natural Populations of Liriodendron chinense (Magnoliaceae) along a Latitudinal Transect

Understanding adaptive genetic variation and its relation to environmental factors are important for understanding how plants adapt to climate change and for managing genetic resources. Genome scans for the loci exhibiting either notably high or low levels of population differentiation (outlier loci) provide one means of identifying genomic regions possibly associated with convergent or divergent selection. In this study, we combined Amplified Fragment Length Polymorphism (AFLP) genome scan and environmental association analysis to test for signals of natural selection in natural populations of Liriodendron chinense (Chinese Tulip Tree; Magnoliaceae) along a latitudinal transect. We genotyped 276 individuals from 11 populations of L. chinense using 987 AFLP markers. Both frequency-based (Dfdist and BayeScan) and correlation-based (MLM) methods were applied to detect outlier loci. Our analyses recovered both neutral and potentially adaptive genetic differentiation among populations of L. chinense. We found moderate genetic diversity within populations and high genetic differentiation among populations with reduced genetic diversity toward the periphery of the species ranges. Nine AFLP marker loci showed evidence of being outliers for population differentiation for both detection methods. Of these, six were strongly associated with at least one climate factor. Temperature, precipitation, and radiation were found to be three important factors influencing local adaptation of L. chinense. The outlier AFLP loci are likely not the target of natural selection, but the neighboring genes of these loci might be involved in local adaptation. Hence, these candidates should be validated by further studies.

Local adaptation to temperature and precipitation in naturally fragmented populations of Cephalotaxus oliveri, an endangered conifer endemic to China

Cephalotaxus oliveri is an endangered tertiary relict conifer endemic to China. The species survives in a wide range from west to east with heterogeneous climatic conditions. Precipitation and temperature are main restrictive factors for distribution of C. oliveri. In order to comprehend the mechanism of adaptive evolution to climate variation, we employed ISSR markers to detect adaptive evolution loci, to identify the association between variation in temperature and precipitation and adaptive loci, and to investigate the genetic structure for 22 C. oliveri natural populations. In total, 14 outlier loci were identified, of which five were associated with temperature and precipitation. Among outlier loci, linkage disequilibrium (LD) was high (42.86%), which also provided strong evidence for selection. In addition, C. oliveri possessed high genetic variation (93.31%) and population differentiation, which may provide raw material to evolution and accelerate local adaptation, respectively. Ecological niche modeling showed that global warming will cause a shift for populations of C. oliveri from south to north with a shrinkage of southern areas. Our results contribute to understand the potential response of conifers to climatic changes, and provide new insights for conifer resource management and conservation strategies.

Landscape genomic analysis of candidate genes for climate adaptation in a California endemic oak, Quercus lobata

PREMISE OF THE STUDY

The ability of California tree populations to survive anthropogenic climate change will be shaped by the geographic structure of adaptive genetic variation. Our goal is to test whether climate-associated candidate genes show evidence of spatially divergent selection in natural populations of valley oak, Quercus lobata, as preliminary indication of local adaptation.

METHODS

Using DNA from 45 individuals from 13 localities across the species' range, we sequenced portions of 40 candidate genes related to budburst/flowering, growth, osmotic stress, and temperature stress. Using 195 single nucleotide polymorphisms (SNPs), we estimated genetic differentiation across populations and correlated allele frequencies with climate gradients using single-locus and multivariate models.

RESULTS

The top 5% of FST estimates ranged from 0.25 to 0.68, yielding loci potentially under spatially divergent selection. Environmental analyses of SNP frequencies with climate gradients revealed three significantly correlated SNPs within budburst/flowering genes and two SNPs within temperature stress genes with mean annual precipitation, after controlling for multiple testing. A redundancy model showed a significant association between SNPs and climate variables and revealed a similar set of SNPs with high loadings on the first axis. In the RDA, climate accounted for 67% of the explained variation, when holding climate constant, in contrast to a putatively neutral SSR data set where climate accounted for only 33%.

CONCLUSIONS

Population differentiation and geographic gradients of allele frequencies in climate-associated functional genes in Q. lobata provide initial evidence of adaptive genetic variation and background for predicting population response to climate change.

Natural Selection and Recombination Rate Variation Shape Nucleotide Polymorphism Across the Genomes of Three Related Populus Species

A central aim of evolutionary genomics is to identify the relative roles that various evolutionary forces have played in generating and shaping genetic variation within and among species. Here we use whole-genome resequencing data to characterize and compare genome-wide patterns of nucleotide polymorphism, site frequency spectrum, and population-scaled recombination rates in three species of Populus: Populus tremula, P. tremuloides, and P. trichocarpa. We find that P. tremuloides has the highest level of genome-wide variation, skewed allele frequencies, and population-scaled recombination rates, whereas P. trichocarpa harbors the lowest. Our findings highlight multiple lines of evidence suggesting that natural selection, due to both purifying and positive selection, has widely shaped patterns of nucleotide polymorphism at linked neutral sites in all three species. Differences in effective population sizes and rates of recombination largely explain the disparate magnitudes and signatures of linked selection that we observe among species. The present work provides the first phylogenetic comparative study on a genome-wide scale in forest trees. This information will also improve our ability to understand how various evolutionary forces have interacted to influence genome evolution among related species.

Transcriptional regulation of the paper mulberry under cold stress as revealed by a comprehensive analysis of transcription factors

BACKGROUND

Several studies have focused on cold tolerance in multiple regulated levels. However, a genome-scale molecular analysis of the regulated network under the control of transcription factors (TFs) is still lacking, especially for trees. To comprehensively identify the TFs that regulate cold stress response in the paper mulberry and understand their regulatory interactions, transcriptomic data was used to assess changes in gene expression induced by exposure to cold.

RESULTS

Results indicated that 794 TFs, belonging to 47 families and comprising more than 59% of the total TFs of this plant, were involved in the cold stress response. They were clustered into three groups, namely early, intermediate and late responsive groups which contained 95, 550 and 149 TFs, respectively. Among of these differentially expressed TFs, one bHLH, two ERFs and three CAMTAs were considered to be the key TFs functioning in the primary signal transduction. After that, at the intermediate stage of cold stress, there were mainly two biological processes that were regulated by TFs, namely cold stress resistance (including 5 bHLH, 14 ERFs, one HSF, 4 MYBs, 3 NACs, 11 WRKYs and so on) and growth and development of lateral organ or apical meristem (including ARR-B, B3, 5 bHLHs, 2 C2H2, 4 CO-like, 2 ERF, 3 HD-ZIP, 3 YABBYs, G2-like, GATA, GRAS and TCP). In late responsive group, 3 ARR-B, C3H, 6 CO-like, 2 G2-like, 2 HSFs, 2 NACs and TCP. Most of them presented the up-regulated expression at 12 or 24 hours after cold stress implied their important roles for the new growth homeostasis under cold stress.

CONCLUSIONS

Our study identified the key TFs that function in the regulatory cascades mediating the activation of downstream genes during cold tress tolerance in the paper mulberry. Based on the analysis, we found that the AP2/ERF, bHLH, MYB, NAC and WRKY families might play the central and significant roles during cold stress response in the paper mulberry just as in other species. Meanwhile, many other TF families previously reported as involving in regulation of growth and development, including ARF, DBB, G2-like, GRF, GRAS, LBD, WOX and YAABY exhibited their important potential function in growth regulation under cold stress.

Genomic sequencing and microsatellite marker development for Boswellia papyrifera, an economically important but threatened tree native to dry tropical forests

Microsatellite (or simple sequence repeat, SSR) markers are highly informative DNA markers often used in conservation genetic research. Next-generation sequencing enables efficient development of large numbers of SSR markers at lower costs. Boswellia papyrifera is an economically important tree species used for frankincense production, an aromatic resinous gum exudate from bark. It grows in dry tropical forests in Africa and is threatened by a lack of rejuvenation. To help guide conservation efforts for this endangered species, we conducted an analysis of its genomic DNA sequences using Illumina paired-end sequencing. The genome size was estimated at 705 Mb per haploid genome. The reads contained one microsatellite repeat per 5.7 kb. Based on a subset of these repeats, we developed 46 polymorphic SSR markers that amplified 2-12 alleles in 10 genotypes. This set included 30 trinucleotide repeat markers, four tetranucleotide repeat markers, six pentanucleotide markers and six hexanucleotide repeat markers. Several markers were cross-transferable to Boswellia pirrotae and B. popoviana. In addition, retrotransposons were identified, the reads were assembled and several contigs were identified with similarity to genes of the terpene and terpenoid backbone synthesis pathways, which form the major constituents of the bark resin.

Genetic differentiation and evolutionary adaptation in Cryptomeria japonica

Local adaptation of plant species is a central issue for survival during global climate change, especially for long-lived forest trees, with their lengthy regeneration time and spatially limited gene flow. Identification of loci and/or genomic regions associated with local adaptation is necessary for knowledge of both evolution and molecular breeding for climate change. Cryptomeria japonica is an important species for forestry in Japan; it has a broad natural distribution and can survive in a range of different environments. The genetic structure of 14 natural populations of this species was investigated using 3930 SNP markers. Populations on the Pacific Ocean side of Japan are clearly different from those on the Japan Sea side, as discussed in previous studies. Structure analysis and population network trees show that peripheral populations, including the most northerly and southerly ones, have unique features. We found that the genetic differentiation coefficient is low, F_ST = 0.05, although it must account for the presence of important genes associated with adaptation to specific environments. In total, 208 outlier loci were detected, of which 43 were associated with environmental variables. Four clumped regions of outlier loci were detected in the genome by linkage analysis. Linkage disequilibrium (LD) was quite high in these clumps of outlier loci, which were found in linkage groups (LGs) 2, 7, 10, and 11, especially between populations of two varieties, and when interchromosomal LD was also detected. The LG7 region is characteristic of the Yakushima population, which is a large, isolated, peripheral population occupying a specific environment resulting from isolation combined with volcanic activity in the region. The detected LD may provide strong evidence for selection between varieties.

Genetic control of functional traits related to photosynthesis and water use efficiency in Pinus pinaster Ait. drought response: integration of genome annotation, allele association and QTL detection for candidate gene identification

BACKGROUND

Understanding molecular mechanisms that control photosynthesis and water use efficiency in response to drought is crucial for plant species from dry areas. This study aimed to identify QTL for these traits in a Mediterranean conifer and tested their stability under drought.

RESULTS

High density linkage maps for Pinus pinaster were used in the detection of QTL for photosynthesis and water use efficiency at three water irrigation regimes. A total of 28 significant and 27 suggestive QTL were found. QTL detected for photochemical traits accounted for the higher percentage of phenotypic variance. Functional annotation of genes within the QTL suggested 58 candidate genes for the analyzed traits. Allele association analysis in selected candidate genes showed three SNPs located in a MYB transcription factor that were significantly associated with efficiency of energy capture by open PSII reaction centers and specific leaf area.

CONCLUSIONS

The integration of QTL mapping of functional traits, genome annotation and allele association yielded several candidate genes involved with molecular control of photosynthesis and water use efficiency in response to drought in a conifer species. The results obtained highlight the importance of maintaining the integrity of the photochemical machinery in P. pinaster drought response.

Intraspecific divergence and evolution of a life-history trade-off along a successional gradient in Hawaii's Metrosideros polymorpha

The importance of environmental gradients in the diversification of long-lived tree species is poorly understood. Two morphologically distinct varieties of the endemic Hawaiian tree, 'ōhi'a lehua (Metrosideros polymorpha), are the canopy dominants at alternate extremes of a successional gradient formed by the recurring disturbance of lava flows on east Hawai'i Island. The maintenance of these varieties despite hybridization may be due to disruptive selection at either end of the successional gradient. To test this hypothesis, seeds from three, replicate monotypic stands of each variety on east Hawai'i Island were germinated and the resulting seedlings grown under four combinations of light and nitrogen levels in a greenhouse, and at early- and late-successional field sites. Growth and survivorship measures revealed differential fitness of these varieties in high- and low-light environments in the greenhouse with corresponding differential fitness in early- and late-successional field sites. Unique light-by-nitrogen interaction effects on growth were observed in each variety, and only the late-successional variety appeared to be nitrogen limited. These two varieties exhibit the classic plant life-history trade-off between fast growth in high light and high survivorship in shade, but notably within a single tree species. These findings strongly implicate a role for Hawaii's striking environmental heterogeneity in the emergence of at least two endemic forms of this woody genus.

In situ genetic association for serotiny, a fire-related trait, in Mediterranean maritime pine (Pinus pinaster)

Wildfire is a major ecological driver of plant evolution. Understanding the genetic basis of plant adaptation to wildfire is crucial, because impending climate change will involve fire regime changes worldwide. We studied the molecular genetic basis of serotiny, a fire-related trait, in Mediterranean maritime pine using association genetics. A single nucleotide polymorphism (SNP) set was used to identify genotype : phenotype associations in situ in an unstructured natural population of maritime pine (eastern Iberian Peninsula) under a mixed-effects model framework. RR-BLUP was used to build predictive models for serotiny in this region. Model prediction power outside the focal region was tested using independent range-wide serotiny data. Seventeen SNPs were potentially associated with serotiny, explaining approximately 29% of the trait phenotypic variation in the eastern Iberian Peninsula. Similar prediction power was found for nearby geographical regions from the same maternal lineage, but not for other genetic lineages. Association genetics for ecologically relevant traits evaluated in situ is an attractive approach for forest trees provided that traits are under strong genetic control and populations are unstructured, with large phenotypic variability. This will help to extend the research focus to ecological keystone non-model species in their natural environments, where polymorphisms acquired their adaptive value.

The evolutionary genetics of the genes underlying phenotypic associations for loblolly pine (Pinus taeda, Pinaceae)

A primary goal of evolutionary genetics is to discover and explain the genetic basis of fitness-related traits and how this genetic basis evolves within natural populations. Unprecedented technological advances have fueled the discovery of genetic variants associated with ecologically relevant phenotypes in many different life forms, as well as the ability to scan genomes for deviations from selectively neutral models of evolution. Theoretically, the degree of overlap between lists of genomic regions identified using each approach is related to the genetic architecture of fitness-related traits and the strength and type of natural selection molding variation at these traits within natural populations. Here we address for the first time in a plant the degree of overlap between these lists, using patterns of nucleotide diversity and divergence for >7000 unique amplicons described from the extensive expressed sequence tag libraries generated for loblolly pine (Pinus taeda L.) in combination with the >1000 published genetic associations. We show that loci associated with phenotypic traits are distinct with regard to neutral expectations. Phenotypes measured at the whole plant level (e.g., disease resistance) exhibit an approximately twofold increase in the proportion of adaptive nonsynonymous substitutions over the genome-wide average. As expected for polygenic traits, these signals were apparent only when loci were considered at the level of functional sets. The ramifications of this result are discussed in light of the continued efforts to dissect the genetic basis of quantitative traits.

Highly local environmental variability promotes intrapopulation divergence of quantitative traits: an example from tropical rain forest trees

BACKGROUND AND AIMS

In habitat mosaics, plant populations face environmental heterogeneity over short geographical distances. Such steep environmental gradients can induce ecological divergence. Lowland rainforests of the Guiana Shield are characterized by sharp, short-distance environmental variations related to topography and soil characteristics (from waterlogged bottomlands on hydromorphic soils to well-drained terra firme on ferralitic soils). Continuous plant populations distributed along such gradients are an interesting system to study intrapopulation divergence at highly local scales. This study tested (1) whether conspecific populations growing in different habitats diverge at functional traits, and (2) whether they diverge in the same way as congeneric species having different habitat preferences.

METHODS

Phenotypic differentiation was studied within continuous populations occupying different habitats for two congeneric, sympatric, and ecologically divergent tree species (Eperua falcata and E. grandiflora, Fabaceae). Over 3000 seeds collected from three habitats were germinated and grown in a common garden experiment, and 23 morphological, biomass, resource allocation and physiological traits were measured.

KEY RESULTS

In both species, seedling populations native of different habitats displayed phenotypic divergence for several traits (including seedling growth, biomass allocation, leaf chemistry, photosynthesis and carbon isotope composition). This may occur through heritable genetic variation or other maternally inherited effects. For a sub-set of traits, the intraspecific divergence associated with environmental variation coincided with interspecific divergence.

CONCLUSIONS

The results indicate that mother trees from different habitats transmit divergent trait values to their progeny, and suggest that local environmental variation selects for different trait optima even at a very local spatial scale. Traits for which differentiation within species follows the same pattern as differentiation between species indicate that the same ecological processes underlie intra- and interspecific variation.

The genomic architecture and association genetics of adaptive characters using a candidate SNP approach in boreal black spruce

BACKGROUND

The genomic architecture of adaptive traits remains poorly understood in non-model plants. Various approaches can be used to bridge this gap, including the mapping of quantitative trait loci (QTL) in pedigrees, and genetic association studies in non-structured populations. Here we present results on the genomic architecture of adaptive traits in black spruce, which is a widely distributed conifer of the North American boreal forest. As an alternative to the usual candidate gene approach, a candidate SNP approach was developed for association testing.

RESULTS

A genetic map containing 231 gene loci was used to identify QTL that were related to budset timing and to tree height assessed over multiple years and sites. Twenty-two unique genomic regions were identified, including 20 that were related to budset timing and 6 that were related to tree height. From results of outlier detection and bulk segregant analysis for adaptive traits using DNA pool sequencing of 434 genes, 52 candidate SNPs were identified and subsequently tested in genetic association studies for budset timing and tree height assessed over multiple years and sites. A total of 34 (65%) SNPs were significantly associated with budset timing, or tree height, or both. Although the percentages of explained variance (PVE) by individual SNPs were small, several significant SNPs were shared between sites and among years.

CONCLUSIONS

The sharing of genomic regions and significant SNPs between budset timing and tree height indicates pleiotropic effects. Significant QTLs and SNPs differed quite greatly among years, suggesting that different sets of genes for the same characters are involved at different stages in the tree's life history. The functional diversity of genes carrying significant SNPs and low observed PVE further indicated that a large number of polymorphisms are involved in adaptive genetic variation. Accordingly, for undomesticated species such as black spruce with natural populations of large effective size and low linkage disequilibrium, efficient marker systems that are predictive of adaptation should require the survey of large numbers of SNPs. Candidate SNP approaches like the one developed in the present study could contribute to reducing these numbers.

The UDP-glucuronate decarboxylase gene family in Populus: structure, expression, and association genetics

In woody crop plants, the oligosaccharide components of the cell wall are essential for important traits such as bioenergy content, growth, and structural wood properties. UDP-glucuronate decarboxylase (UXS) is a key enzyme in the synthesis of UDP-xylose for the formation of xylans during cell wall biosynthesis. Here, we isolated a multigene family of seven members (PtUXS1-7) encoding UXS from Populus tomentosa, the first investigation of UXSs in a tree species. Analysis of gene structure and phylogeny showed that the PtUXS family could be divided into three groups (PtUXS1/4, PtUXS2/5, and PtUXS3/6/7), consistent with the tissue-specific expression patterns of each PtUXS. We further evaluated the functional consequences of nucleotide polymorphisms in PtUXS1. In total, 243 single-nucleotide polymorphisms (SNPs) were identified, with a high frequency of SNPs (1/18 bp) and nucleotide diversity (π_T = 0.01033, θ_w =0.01280). Linkage disequilibrium (LD) analysis showed that LD did not extend over the entire gene (r²<0.1, P<0.001, within 700 bp). SNP- and haplotype-based association analysis showed that nine SNPs (Q <0.10) and 12 haplotypes (P<0.05) were significantly associated with growth and wood property traits in the association population (426 individuals), with 2.70% to 12.37% of the phenotypic variation explained. Four significant single-marker associations (Q <0.10) were validated in a linkage mapping population of 1200 individuals. Also, RNA transcript accumulation varies among genotypic classes of SNP10 was further confirmed in the association population. This is the first comprehensive study of the UXS gene family in woody plants, and lays the foundation for genetic improvements of wood properties and growth in trees using genetic engineering or marker-assisted breeding.

Molecular divergence in tropical tree populations occupying environmental mosaics

Unveiling the genetic basis of local adaptation to environmental variation is a major goal in molecular ecology. In rugged landscapes characterized by environmental mosaics, living populations and communities can experience steep ecological gradients over very short geographical distances. In lowland tropical forests, interspecific divergence in edaphic specialization (for seasonally flooded bottomlands and seasonally dry terra firme soils) has been proven by ecological studies on adaptive traits. Some species are nevertheless capable of covering the entire span of the gradient; intraspecific variation for adaptation to contrasting conditions may explain the distribution of such ecological generalists. We investigated whether local divergence happens at small spatial scales in two stands of Eperua falcata (Fabaceae), a widespread tree species of the Guiana Shield. We investigated Single Nucleotide Polymorphisms (SNP) and sequence divergence as well as spatial genetic structure (SGS) at four genes putatively involved in stress response and three genes with unknown function. Significant genetic differentiation was observed among sub-populations within stands, and eight SNP loci showed patterns compatible with disruptive selection. SGS analysis showed genetic turnover along the gradients at three loci, and at least one haplotype was found to be in repulsion with one habitat. Taken together, these results suggest genetic differentiation at small spatial scale in spite of gene flow. We hypothesize that heterogeneous environments may cause molecular divergence, possibly associated to local adaptation in E. falcata.

Targeted enrichment of the black cottonwood (Populus trichocarpa) gene space using sequence capture

Background

High-throughput re-sequencing is rapidly becoming the method of choice for studies of neutral and adaptive processes in natural populations across taxa. As re-sequencing the genome of large numbers of samples is still cost-prohibitive in many cases, methods for genome complexity reduction have been developed in attempts to capture most ecologically-relevant genetic variation. One of these approaches is sequence capture, in which oligonucleotide baits specific to genomic regions of interest are synthesized and used to retrieve and sequence those regions.

Results

We used sequence capture to re-sequence most predicted exons, their upstream regulatory regions, as well as numerous random genomic intervals in a panel of 48 genotypes of the angiosperm tree Populus trichocarpa (black cottonwood, or ‘poplar’). A total of 20.76Mb (5%) of the poplar genome was targeted, corresponding to 173,040 baits. With 12 indexed samples run in each of four lanes on an Illumina HiSeq instrument (2x100 paired-end), 86.8% of the bait regions were on average sequenced at a depth ≥10X. Few off-target regions (>250bp away from any bait) were present in the data, but on average ~80bp on either side of the baits were captured and sequenced to an acceptable depth (≥10X) to call heterozygous SNPs. Nucleotide diversity estimates within and adjacent to protein-coding genes were similar to those previously reported in Populus spp., while intergenic regions had higher values consistent with a relaxation of selection.

Conclusions

Our results illustrate the efficiency and utility of sequence capture for re-sequencing highly heterozygous tree genomes, and suggest design considerations to optimize the use of baits in future studies.

Genome resequencing reveals multiscale geographic structure and extensive linkage disequilibrium in the forest tree Populus trichocarpa

• Plant population genomics informs evolutionary biology, breeding, conservation and bioenergy feedstock development. For example, the detection of reliable phenotype-genotype associations and molecular signatures of selection requires a detailed knowledge about genome-wide patterns of allele frequency variation, linkage disequilibrium and recombination. • We resequenced 16 genomes of the model tree Populus trichocarpa and genotyped 120 trees from 10 subpopulations using 29,213 single-nucleotide polymorphisms. • Significant geographic differentiation was present at multiple spatial scales, and range-wide latitudinal allele frequency gradients were strikingly common across the genome. The decay of linkage disequilibrium with physical distance was slower than expected from previous studies in Populus, with r(2) dropping below 0.2 within 3-6 kb. Consistent with this, estimates of recent effective population size from linkage disequilibrium (N(e) ≈ 4000-6000) were remarkably low relative to the large census sizes of P. trichocarpa stands. Fine-scale rates of recombination varied widely across the genome, but were largely predictable on the basis of DNA sequence and methylation features. • Our results suggest that genetic drift has played a significant role in the recent evolutionary history of P. trichocarpa. Most importantly, the extensive linkage disequilibrium detected suggests that genome-wide association studies and genomic selection in undomesticated populations may be more feasible in Populus than previously assumed.

Annotated genetic linkage maps of Pinus pinaster Ait. from a Central Spain population using microsatellite and gene based markers

BACKGROUND

Pinus pinaster Ait. is a major resin producing species in Spain. Genetic linkage mapping can facilitate marker-assisted selection (MAS) through the identification of Quantitative Trait Loci and selection of allelic variants of interest in breeding populations. In this study, we report annotated genetic linkage maps for two individuals (C14 and C15) belonging to a breeding program aiming to increase resin production. We use different types of DNA markers, including last-generation molecular markers.

RESULTS

We obtained 13 and 14 linkage groups for C14 and C15 maps, respectively. A total of 211 and 215 markers were positioned on each map and estimated genome length was between 1,870 and 2,166 cM respectively, which represents near 65% of genome coverage. Comparative mapping with previously developed genetic linkage maps for P. pinaster based on about 60 common markers enabled aligning linkage groups to this reference map. The comparison of our annotated linkage maps and linkage maps reporting QTL information revealed 11 annotated SNPs in candidate genes that co-localized with previously reported QTLs for wood properties and water use efficiency.

CONCLUSIONS

This study provides genetic linkage maps from a Spanish population that shows high levels of genetic divergence with French populations from which segregating progenies have been previously mapped. These genetic maps will be of interest to construct a reliable consensus linkage map for the species. The importance of developing functional genetic linkage maps is highlighted, especially when working with breeding populations for its future application in MAS for traits of interest.

Predicting adaptive phenotypes from multilocus genotypes in Sitka spruce (Picea sitchensis) using random forest

Climate is the primary driver of the distribution of tree species worldwide, and the potential for adaptive evolution will be an important factor determining the response of forests to anthropogenic climate change. Although association mapping has the potential to improve our understanding of the genomic underpinnings of climatically relevant traits, the utility of adaptive polymorphisms uncovered by such studies would be greatly enhanced by the development of integrated models that account for the phenotypic effects of multiple single-nucleotide polymorphisms (SNPs) and their interactions simultaneously. We previously reported the results of association mapping in the widespread conifer Sitka spruce (Picea sitchensis). In the current study we used the recursive partitioning algorithm ‘Random Forest’ to identify optimized combinations of SNPs to predict adaptive phenotypes. After adjusting for population structure, we were able to explain 37% and 30% of the phenotypic variation, respectively, in two locally adaptive traits—autumn budset timing and cold hardiness. For each trait, the leading five SNPs captured much of the phenotypic variation. To determine the role of epistasis in shaping these phenotypes, we also used a novel approach to quantify the strength and direction of pairwise interactions between SNPs and found such interactions to be common. Our results demonstrate the power of Random Forest to identify subsets of markers that are most important to climatic adaptation, and suggest that interactions among these loci may be widespread.

Genome scanning for detecting adaptive genes along environmental gradients in the Japanese conifer, Cryptomeria japonica

Local adaptation is important in evolutionary processes and speciation. We used multiple tests to identify several candidate genes that may be involved in local adaptation from 1026 loci in 14 natural populations of Cryptomeria japonica, the most economically important forestry tree in Japan. We also studied the relationships between genotypes and environmental variables to obtain information on the selective pressures acting on individual populations. Outlier loci were mapped onto a linkage map, and the positions of loci associated with specific environmental variables are considered. The outlier loci were not randomly distributed on the linkage map; linkage group 11 was identified as a genomic island of divergence. Three loci in this region were also associated with environmental variables such as mean annual temperature, daily maximum temperature, maximum snow depth, and so on. Outlier loci identified with high significance levels will be essential for conservation purposes and for future work on molecular breeding.

Spatially variable natural selection and the divergence between parapatric subspecies of lodgepole pine (Pinus contorta, Pinaceae)

PREMISE OF THE STUDY

Plant populations arrayed across sharp environmental gradients are ideal systems for identifying the genetic basis of ecologically relevant phenotypes. A series of five uplifted marine terraces along the northern coast of California represents one such system where morphologically distinct populations of lodgepole pine (Pinus contorta) are distributed across sharp soil gradients ranging from fertile soils near the coast to podzolic soils ca. 5 km inland.

METHODS

A total of 92 trees was sampled across four coastal marine terraces (N = 10-46 trees/terrace) located in Mendocino County, California and sequenced for a set of 24 candidate genes for growth and responses to various soil chemistry variables. Statistical analyses relying on patterns of nucleotide diversity were employed to identify genes whose diversity patterns were inconsistent with three null models.

KEY RESULTS

Most genes displayed patterns of nucleotide diversity that were consistent with null models (N = 19) or with the presence of paralogs (N = 3). Two genes, however, were exceptional: an aluminum responsive ABC-transporter with F(ST) = 0.664 and an inorganic phosphate transporter characterized by divergent haplotypes segregating at intermediate frequencies in most populations.

CONCLUSIONS

Spatially variable natural selection along gradients of aluminum and phosphate ion concentrations likely accounted for both outliers. These results shed light on some of the genetic components comprising the extended phenotype of this ecosystem, as well as highlight ecotones as fruitful study systems for the detection of adaptive genetic variants.

Allele discovery of ten candidate drought-response genes in Austrian oak using a systematically informatics approach based on 454 amplicon sequencing

Background

Rise of temperatures and shortening of available water as result of predicted climate change will impose significant pressure on long-lived forest tree species. Discovering allelic variation present in drought related genes of two Austrian oak species can be the key to understand mechanisms of natural selection and provide forestry with key tools to cope with future challenges.

Results

In the present study we have used Roche 454 sequencing and developed a bioinformatic pipeline to process multiplexed tagged amplicons in order to identify single nucleotide polymorphisms and allelic sequences of ten candidate genes related to drought/osmotic stress from sessile oak (Quercus robur) and sessile oak (Q. petraea) individuals. Out of these, eight genes of 336 oak individuals growing in Austria have been detected with a total number of 158 polymorphic sites. Allele numbers ranged from ten to 52 with observed heterozygosity ranging from 0.115 to 0.640. All loci deviated from Hardy-Weinberg equilibrium and linkage disequilibrium was found among six combinations of loci.

Conclusions

We have characterized 183 alleles of drought related genes from oak species and detected first evidences of natural selection. Beside the potential for marker development, we have created an expandable bioinformatic pipeline for the analysis of next generation sequencing data.

Functional genomics of a living fossil tree, Ginkgo, based on next-generation sequencing technology

Ginkgo biloba is monotypic species native to China and has old, dioecious, medicinally important characteristics. The functional genes related to these characteristics have not been effectively explored due to a limited number of expressed sequence tags (ESTs) from Ginkgo. To discover novel functional genes efficiently and to understand the development of a living fossil tree, Ginkgo, we used massive parallel pyrosequencing on the Roche 454 GS FLX Titanium platform to generate 64 057 ESTs. The ESTs combined with the 21 590 Ginkgo ESTs in genbank were assembled into 22 304 unique putative transcripts, in which 13 922 novel unique putative transcripts were identified by 454 sequencing. After being assigned to putative functions with Gene Ontology terms, a detailed view of the Ginkgo biological systems was displayed, including characterization of unique putative transcripts with homology to known key enzymes and transcription factors involved in ginkgolide/bilobalide and flavonoid biosynthetic pathways, as well as unique putative transcripts related to development, response to disease and defence. The fact that three full-length Ginkgo genes encoding key enzymes were found and cloned, suggests that high-throughput sequencing technology is superior to traditional gene-by-gene approach in discovery of genes. Additionally, a total of 204 simple sequence repeat motifs were detected. Our study not only lays the foundations for transcriptome-led studies in biosynthetic mechanisms, but also contributes significantly to the understanding of functional genomics and development in non-model plants.

The heterogeneous levels of linkage disequilibrium in white spruce genes and comparative analysis with other conifers

In plants, knowledge about linkage disequilibrium (LD) is relevant for the design of efficient single-nucleotide polymorphism arrays in relation to their use in population and association genomics studies. Previous studies of conifer genes have shown LD to decay rapidly within gene limits, but exceptions have been reported. To evaluate the extent of heterogeneity of LD among conifer genes and its potential causes, we examined LD in 105 genes of white spruce (Picea glauca) by sequencing a panel of 48 haploid megagametophytes from natural populations and further compared it with LD in other conifer species. The average pairwise r(2) value was 0.19 (s.d.=0.19), and LD dropped quickly with a half-decay being reached at a distance of 65 nucleotides between sites. However, LD was significantly heterogeneous among genes. A first group of 29 genes had stronger LD (mean r(2)=0.28), and a second group of 38 genes had weaker LD (mean r(2)=0.12). While a strong relationship was found with the recombination rate, there was no obvious relationship between LD and functional classification. The level of nucleotide diversity, which was highly heterogeneous across genes, was also not significantly correlated with LD. A search for selection signatures highlighted significant deviations from the standard neutral model, which could be mostly attributed to recent demographic changes. Little evidence was seen for hitchhiking and clear relationships with LD. When compared among conifer species, on average, levels of LD were similar in genes from white spruce, Norway spruce and Scots pine, whereas loblolly pine and Douglas fir genes exhibited a significantly higher LD.