Rolling down that mountain: microgeographical adaptive divergence during a fast population expansion along a steep environmental gradient in European beech.

Forest tree populations harbour high genetic diversity thanks to large effective population sizes and strong gene flow, allowing them to diversify through adaptation to local environmental pressures within dispersal distance. Many tree populations also experienced historical demographic fluctuations, including spatial population contraction or expansions at various temporal scales, which may constrain their ability to adapt to environmental variations. Our aim is to investigate how recent contraction and expansion events interfere with local adaptation, by studying patterns of adaptive divergence between closely related stands undergoing environmentally contrasted conditions, and having or not recently expanded. To investigate genome-wide signatures of local adaptation while accounting for demography, we analysed divergence in a European beech population by testing pairwise differentiation among four tree stands at ~35k Single Nucleotide Polymorphisms from ~9k genomic regions. We applied three divergence outlier search methods resting on different assumptions and targeting either single SNPs or contiguous genomic regions, while accounting for the effect of population size variations on genetic divergence. We found 27 signals of selective signatures in 19 target regions. Putatively adaptive divergence involved all stand pairs. We retrieved signals both when comparing old-growth stands and recently colonised areas and when comparing stands within the old-growth area. Therefore, adaptive divergence processes have taken place both over short time spans, under strong environmental contrasts, and over short ecological gradients, in populations that have been stable in the long term. This suggests that standing genetic variation supports local, microgeographic divergence processes, which can maintain genetic diversity at the landscape level.

Common microgeographical selection patterns revealed in four European conifers.

Microgeographical adaptation occurs when the effects of directional selection persist despite gene flow. Traits and genetic loci under selection can then show adaptive divergence, against the backdrop of little differentiation at other traits or loci. How common such events are and how strong the selection is that underlies them remain open questions. Here, we discovered and analysed microgeographical patterns of genomic divergence in four European and Mediterranean conifers with widely differing life-history traits and ecological requirements (Abies alba MIll., Cedrus atlantica [Endl.] Manetti, Pinus halepensis Mill. and Pinus pinaster Aiton) by screening pairs from geographically close forest stands sampled along steep ecological gradients. We inferred patterns of genomic divergence by applying a combination of divergence outlier detection methods, demographic modelling, Approximate Bayesian Computation inferences and genomic annotation to genomic data. Surprisingly for such small geographical scales, we showed that selection is strong in all species but generally affects different loci in each. A clear signature of selection was systematically detected on a fraction of the genome, of the order of 0.1%-1% of the loci depending on the species. The novel modelling method we designed for estimating selection coefficients showed that the microgeographical selection coefficient scaled by population size (Ns) was 2-30. Our results convincingly suggest that selection maintains within-population diversity at microgeographical scales in spatially heterogeneous environments. Such genetic diversity is likely to be a major reservoir of adaptive potential, helping populations to adapt under fluctuating environmental conditions.

Interactions between microenvironment, selection and genetic architecture drive multiscale adaptation in a simulation experiment.

When environmental conditions differ both within and among populations, multiscale adaptation results from processes at both scales and interference across scales. We hypothesize that within-population environmental heterogeneity influences the chance of success of migration events, both within and among populations, and maintains within-population adaptive differentiation. We used a simulation approach to analyse the joint effects of environmental heterogeneity patterns, selection intensity and number of QTL controlling a selected trait on local adaptation in a hierarchical metapopulation design. We show the general effects of within-population environmental heterogeneity: (i) it increases occupancy rate at the margins of distribution ranges, under extreme environments and high levels of selection; (ii) it increases the adaptation lag in all environments; (iii) it impacts the genetic variance in each environment, depending on the ratio of within- to between-populations environmental heterogeneity; (iv) it reduces the selection-induced erosion of adaptive gene diversity. Most often, the smaller the number of QTL involved, the stronger are these effects. We also show that both within- and between-populations phenotypic differentiation (Q ST ) mainly results from covariance of QTL effects rather than QTL differentiation (F STq ), that within-population QTL differentiation is negligible, and that stronger divergent selection is required to produce adaptive differentiation within populations than among populations. With a high number of QTL, when the difference between environments within populations exceeds the smallest difference between environments across populations, high levels of within-population differentiation can be reached, reducing differentiation among populations. Our study stresses the need to account for within-population environmental heterogeneity when investigating local adaptation.

Genomic and phenotypic divergence unveil microgeographic adaptation in the Amazonian hyperdominant tree Eperua falcata Aubl. (Fabaceae).

Plant populations can undergo very localized adaptation, allowing widely distributed populations to adapt to divergent habitats in spite of recurrent gene flow. Neotropical trees-whose large and undisturbed populations often span a variety of environmental conditions and local habitats-are particularly good models to study this process. Here, we explore patterns of adaptive divergence from large (i.e., regional) to small (i.e., microgeographic) spatial scales in the hyperdominant Amazonian tree Eperua falcata Aubl. (Fabaceae) under a replicated design involving two microhabitats (~300 m apart) in two study sites (~300 km apart). A three-year reciprocal transplant illustrates that, beyond strong maternal effects and phenotypic plasticity, genetically driven divergence in seedling growth and leaf traits was detected both between seedlings originating from different regions, and between seedlings from different microhabitats. In parallel, a complementary genome scan for selection was carried out through whole-genome sequencing of tree population pools. A set of 290 divergence outlier SNPs was detected at the regional scale (between study sites), while 185 SNPs located in the vicinity of 106 protein-coding genes were detected as replicated outliers between microhabitats within regions. Outlier-surrounding genomic regions are involved in a variety of physiological processes, including plant responses to stress (e.g., oxidative stress, hypoxia and metal toxicity) and biotic interactions. Together with evidence of microgeographic divergence in functional traits, the discovery of genomic candidates for microgeographic adaptive divergence represents a promising advance in our understanding of local adaptation, which probably operates across multiple spatial scales and underpins divergence and diversification in Neotropical trees.

Genetic signatures of divergent selection in European beech (Fagus sylvatica L.) are associated with the variation in temperature and precipitation across its distribution range.

High genetic variation and extensive gene flow may help forest trees with adapting to ongoing climate change, yet the genetic bases underlying their adaptive potential remain largely unknown. We investigated range-wide patterns of potentially adaptive genetic variation in 64 populations of European beech (Fagus sylvatica L.) using 270 SNPs from 139 candidate genes involved either in phenology or in stress responses. We inferred neutral genetic structure and processes (drift and gene flow) and performed differentiation outlier analyses and gene-environment association (GEA) analyses to detect signatures of divergent selection. Beech range-wide genetic structure was consistent with the species' previously identified postglacial expansion scenario and recolonization routes. Populations showed high diversity and low differentiation along the major expansion routes. A total of 52 loci were found to be putatively under selection and 15 of them turned up in multiple GEA analyses. Temperature and precipitation related variables were equally represented in significant genotype-climate associations. Signatures of divergent selection were detected in the same proportion for stress response and phenology-related genes. The range-wide adaptive genetic structure of beech appears highly integrated, suggesting a balanced contribution of phenology and stress-related genes to local adaptation, and of temperature and precipitation regimes to genetic clines. Our results imply a best-case scenario for the maintenance of high genetic diversity during range shifts in beech (and putatively other forest trees) with a combination of gene flow maintaining within-population neutral diversity and selection maintaining between-population adaptive differentiation.

Tropical rainforests that persisted: inferences from the Quaternary demographic history of eight tree species in the Guiana shield.

How Quaternary climatic and geological disturbances influenced the composition of Neotropical forests is hotly debated. Rainfall and temperature changes during and/or immediately after the last glacial maximum (LGM) are thought to have strongly affected the geographical distribution and local abundance of tree species. The paucity of the fossil records in Neotropical forests prevents a direct reconstruction of such processes. To describe community-level historical trends in forest composition, we turned therefore to inferential methods based on the reconstruction of past demographic changes. In particular, we modelled the history of rainforests in the eastern Guiana Shield over a timescale of several thousand generations, through the application of approximate Bayesian computation and maximum-likelihood methods to diversity data at nuclear and chloroplast loci in eight species or subspecies of rainforest trees. Depending on the species and on the method applied, we detected population contraction, expansion or stability, with a general trend in favour of stability or expansion, with changes presumably having occurred during or after the LGM. These findings suggest that Guiana Shield rainforests have globally persisted, while expanding, through the Quaternary, but that different species have experienced different demographic events, with a trend towards the increase in frequency of light-demanding, disturbance-associated species.

High-throughput transcriptome sequencing and preliminary functional analysis in four Neotropical tree species.

BACKGROUND: The Amazonian rainforest is predicted to suffer from ongoing environmental changes. Despite the need to evaluate the impact of such changes on tree genetic diversity, we almost entirely lack genomic resources. RESULTS: In this study, we analysed the transcriptome of four tropical tree species (Carapa guianensis, Eperua falcata, Symphonia globulifera and Virola michelii) with contrasting ecological features, belonging to four widespread botanical families (respectively Meliaceae, Fabaceae, Clusiaceae and Myristicaceae). We sequenced cDNA libraries from three organs (leaves, stems, and roots) using 454 pyrosequencing. We have developed an R and bioperl-based bioinformatic procedure for de novo assembly, gene functional annotation and marker discovery. Mismatch identification takes into account single-base quality values as well as the likelihood of false variants as a function of contig depth and number of sequenced chromosomes. Between 17103 (for Symphonia globulifera) and 23390 (for Eperua falcata) contigs were assembled. Organs varied in the numbers of unigenes they apparently express, with higher number in roots. Patterns of gene expression were similar across species, with metabolism of aromatic compounds standing out as an overrepresented gene function. Transcripts corresponding to several gene functions were found to be over- or underrepresented in each organ. We identified between 4434 (for Symphonia globulifera) and 9076 (for Virola surinamensis) well-supported mismatches. The resulting overall mismatch density was comprised between 0.89 (S. globulifera) and 1.05 (V. surinamensis) mismatches/100 bp in variation-containing contigs. CONCLUSION: The relative representation of gene functions in the four transcriptomes suggests that secondary metabolism may be particularly important in tropical trees. The differential representation of transcripts among tissues suggests differential gene expression, which opens the way to functional studies in these non-model, ecologically important species. We found substantial amounts of mismatches in the four species. These newly identified putative variants are a first step towards acquiring much needed genomic resources for tropical tree species.

Estimation of contemporary effective population size in plant populations: Limitations of genomic datasets.

Effective population size (N e) is a pivotal evolutionary parameter with crucial implications in conservation practice and policy. Genetic methods to estimate N e have been preferred over demographic methods because they rely on genetic data rather than time-consuming ecological monitoring. Methods based on linkage disequilibrium (LD), in particular, have become popular in conservation as they require a single sampling and provide estimates that refer to recent generations. A software program based on the LD method, GONE, looks particularly promising to estimate contemporary and recent-historical N e (up to 200 generations in the past). Genomic datasets from non-model species, especially plants, may present some constraints to the use of GONE, as linkage maps and reference genomes are seldom available, and SNP genotyping is usually based on reduced-representation methods. In this study, we use empirical datasets from four plant species to explore the limitations of plant genomic datasets when estimating N e using the algorithm implemented in GONE, in addition to exploring some typical biological limitations that may affect N e estimation using the LD method, such as the occurrence of population structure. We show how accuracy and precision of N e estimates potentially change with the following factors: occurrence of missing data, limited number of SNPs/individuals sampled, and lack of information about the location of SNPs on chromosomes, with the latter producing a significant bias, previously unexplored with empirical data. We finally compare the N e estimates obtained with GONE for the last generations with the contemporary N e estimates obtained with the programs currentNe and NeEstimator.

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.

Fire frequency, as well as stress response and developmental gene control serotiny level variation in a widespread pioneer Mediterranean conifer, Pinus halepensis.

Many plants undergo adaptation to fire. Yet, as global change is increasing fire frequency worldwide, our understanding of the genetics of adaptation to fire is still limited. We studied the genetic basis of serotiny (the ability to disseminate seeds exclusively after fire) in the widespread, pioneer Mediterranean conifer Pinus halepensis Mill., by linking individual variation in serotiny presence and level to fire frequency and to genetic polymorphism in natural populations. After filtering steps, 885 single nucleotide polymorphisms (SNPs) out of 8000 SNPs used for genotyping were implemented to perform an in situ association study between genotypes and serotiny presence and level. To identify serotiny-associated loci, we performed random forest analyses of the effect of SNPs on serotiny levels, while controlling for tree size, frequency of wildfires, and background environmental parameters. Serotiny showed a bimodal distribution, with serotinous trees more frequent in populations exposed to fire in their recent history. Twenty-two SNPs found in genes involved in stress tolerance were associated with the presence-absence of serotiny while 37 found in genes controlling for flowering were associated with continuous serotiny variation. This study shows the high potential of P. halepensis to adapt to changing fire regimes, benefiting from a large and flexible genetic basis of trait variation.

Signature of mid-Pleistocene lineages in the European silver fir (Abies alba Mill.) at its geographic distribution margin.

In a conservation and sustainable management perspective, we identify the ecological, climatic, and demographic factors responsible for the genetic diversity patterns of the European silver fir (Abies alba Mill.) at its southwestern range margin (Pyrenees Mountains, France, Europe). We sampled 45 populations throughout the French Pyrenees and eight neighboring reference populations in the Massif Central, Alps, and Corsica. We genotyped 1,620 individuals at three chloroplast and ten nuclear microsatellite loci. We analyzed within- and among-population genetic diversity using phylogeographic reconstructions, tests of isolation-by-distance, Bayesian population structure inference, modeling of demographic scenarios, and regression analyses of genetic variables with current and past environmental variables. Genetic diversity decreased from east to west suggesting isolation-by-distance from the Alps to the Pyrenees and from the Eastern to the Western Pyrenees. We identified two Pyrenean lineages that diverged from a third Alpine-Corsica-Massif Central lineage 0.8 to 1.1 M years ago and subsequently formed a secondary contact zone in the Central Pyrenees. Population sizes underwent contrasted changes, with a contraction in the west and an expansion in the east. Glacial climate affected the genetic composition of the populations, with the western genetic cluster only observed in locations corresponding to the coldest past climate and highest elevations. The eastern cluster was observed over a larger range of temperatures and elevations. All demographic events shaping the current spatial structure of genetic diversity took place during the Mid-Pleistocene Transition, long before the onset of the Holocene. The Western Pyrenees lineage may require additional conservation efforts, whereas the eastern lineage is well protected in in situ gene conservation units. Due to past climate oscillations and the likely emergence of independent refugia, east-west oriented mountain ranges may be important reservoir of genetic diversity in a context of past and ongoing climate change in Europe.

A Chromosome-Level Genome Assembly of the European Beech (Fagus sylvatica) Reveals Anomalies for Organelle DNA Integration, Repeat Content and Distribution of SNPs.

The European Beech is the dominant climax tree in most regions of Central Europe and valued for its ecological versatility and hardwood timber. Even though a draft genome has been published recently, higher resolution is required for studying aspects of genome architecture and recombination. Here, we present a chromosome-level assembly of the more than 300 year-old reference individual, Bhaga, from the Kellerwald-Edersee National Park (Germany). Its nuclear genome of 541 Mb was resolved into 12 chromosomes varying in length between 28 and 73 Mb. Multiple nuclear insertions of parts of the chloroplast genome were observed, with one region on chromosome 11 spanning more than 2 Mb which fragments up to 54,784 bp long and covering the whole chloroplast genome were inserted randomly. Unlike in Arabidopsis thaliana, ribosomal cistrons are present in Fagus sylvatica only in four major regions, in line with FISH studies. On most assembled chromosomes, telomeric repeats were found at both ends, while centromeric repeats were found to be scattered throughout the genome apart from their main occurrence per chromosome. The genome-wide distribution of SNPs was evaluated using a second individual from Jamy Nature Reserve (Poland). SNPs, repeat elements and duplicated genes were unevenly distributed in the genomes, with one major anomaly on chromosome 4. The genome presented here adds to the available highly resolved plant genomes and we hope it will serve as a valuable basis for future research on genome architecture and for understanding the past and future of European Beech populations in a changing climate.

Aquaporins in the wild: natural genetic diversity and selective pressure in the PIP gene family in five Neotropical tree species.

BACKGROUND: Tropical trees undergo severe stress through seasonal drought and flooding, and the ability of these species to respond may be a major factor in their survival in tropical ecosystems, particularly in relation to global climate change. Aquaporins are involved in the regulation of water flow and have been shown to be involved in drought response; they may therefore play a major adaptive role in these species. We describe genetic diversity in the PIP sub-family of the widespread gene family of Aquaporins in five Neotropical tree species covering four botanical families. RESULTS: PIP Aquaporin subfamily genes were isolated, and their DNA sequence polymorphisms characterised in natural populations. Sequence data were analysed with statistical tests of standard neutral equilibrium and demographic scenarios simulated to compare with the observed results. Chloroplast SSRs were also used to test demographic transitions. Most gene fragments are highly polymorphic and display signatures of balancing selection or bottlenecks; chloroplast SSR markers have significant statistics that do not conform to expectations for population bottlenecks. Although not incompatible with a purely demographic scenario, the combination of all tests tends to favour a selective interpretation of extant gene diversity. CONCLUSIONS: Tropical tree PIP genes may generally undergo balancing selection, which may maintain high levels of genetic diversity at these loci. Genetic variation at PIP genes may represent a response to variable environmental conditions.

Genomic patterns of adaptive divergence between chromosomally differentiated sunflower species.

Understanding the genetic mechanisms of speciation and basis of species differences is among the most important challenges in evolutionary biology. Two questions of particular interest are what roles divergent selection and chromosomal differentiation play in these processes. A number of recently proposed theories argue that chromosomal rearrangements can facilitate the development and maintenance of reproductive isolation and species differences by suppressing recombination within rearranged regions. Reduced recombination permits the accumulation of alleles contributing to isolation and adaptive differentiation and protects existing differences from the homogenizing effects of introgression between incipient species. Here, we examine patterns of genetic diversity and divergence in rearranged versus collinear regions in two widespread, extensively hybridizing sunflower species, Helianthus annuus and Helianthus petiolaris, using sequence data from 77 loci distributed throughout the genomes of the two species. We find weak evidence for increased genetic divergence near chromosomal break points but not within rearranged regions overall. We find no evidence for increased rates of adaptive divergence on rearranged chromosomes; in fact, collinear chromosomes show a far greater excess of fixed amino acid differences between the two species. A comparison with a third sunflower species indicates that much of the nonsynonymous divergence between H. annuus and H. petiolaris probably occurred during or soon after their formation. Our results suggest a limited role for chromosomal rearrangements in genetic divergence, but they do document substantial adaptive divergence and provide further evidence of how species integrity and genetic identity can be maintained at many loci in the face of extensive hybridization and gene flow.

Mechanosensing of stem bending and its interspecific variability in five neotropical rainforest species.

BACKGROUND AND AIMS: In rain forests, sapling survival is highly dependent on the regulation of trunk slenderness (height/diameter ratio): shade-intolerant species have to grow in height as fast as possible to reach the canopy but also have to withstand mechanical loadings (wind and their own weight) to avoid buckling. Recent studies suggest that mechanosensing is essential to control tree dimensions and stability-related morphogenesis. Differences in species slenderness have been observed among rainforest trees; the present study thus investigates whether species with different slenderness and growth habits exhibit differences in mechanosensitivity. METHODS: Recent studies have led to a model of mechanosensing (sum-of-strains model) that predicts a quantitative relationship between the applied sum of longitudinal strains and the plant's responses in the case of a single bending. Saplings of five different neotropical species (Eperua falcata, E. grandiflora, Tachigali melinonii, Symphonia globulifera and Bauhinia guianensis) were subjected to a regimen of controlled mechanical loading phases (bending) alternating with still phases over a period of 2 months. Mechanical loading was controlled in terms of strains and the five species were subjected to the same range of sum of strains. The application of the sum-of-strain model led to a dose-response curve for each species. Dose-response curves were then compared between tested species. KEY RESULTS: The model of mechanosensing (sum-of-strain model) applied in the case of multiple bending as long as the bending frequency was low. A comparison of dose-response curves for each species demonstrated differences in the stimulus threshold, suggesting two groups of responses among the species. Interestingly, the liana species B. guianensis exhibited a higher threshold than other Leguminosae species tested. CONCLUSIONS: This study provides a conceptual framework to study variability in plant mechanosensing and demonstrated interspecific variability in mechanosensing.

Discrepancies between genetic and ecological divergence patterns suggest a complex biogeographic history in a Neotropical genus.

Phylogenetic patterns and the underlying speciation processes can be deduced from morphological, functional, and ecological patterns of species similarity and divergence. In some cases, though, species retain multiple similarities and remain almost indistinguishable; in other cases, evolutionary convergence can make such patterns misleading; very often in such cases, the "true" picture only emerges from carefully built molecular phylogenies, which may come with major surprises. In addition, closely related species may experience gene flow after divergence, thus potentially blurring species delimitation. By means of advanced inferential methods, we studied molecular divergence between species of the Virola genus (Myristicaceae): widespread Virola michelii and recently described, endemic V. kwatae, using widespread V. surinamensis as a more distantly related outgroup with different ecology and morphology-although with overlapping range. Contrary to expectations, we found that the latter, and not V. michelii, was sister to V. kwatae. Therefore, V. kwatae probably diverged from V. surinamensis through a recent morphological and ecological shift, which brought it close to distantly related V. michelii. Through the modeling of the divergence process, we inferred that gene flow between V. surinamensis and V. kwatae stopped soon after their divergence and resumed later, in a classical secondary contact event which did not erase their ecological and morphological differences. While we cannot exclude that initial divergence occurred in allopatry, current species distribution and the absence of geographical barriers make complete isolation during speciation unlikely. We tentatively conclude that (a) it is possible that divergence occurred in allopatry/parapatry and (b) secondary contact did not suppress divergence.

Microgeographic local adaptation and ecotype distributions: The role of selective processes on early life-history traits in sympatric, ecologically divergent Symphonia populations.

Trees are characterized by the large number of seeds they produce. Although most of those seeds will never germinate, plenty will. Of those which germinate, many die young, and eventually, only a minute fraction will grow to adult stage and reproduce. Is this just a random process? Do variations in germination and survival at very young stages rely on variations in adaptations to microgeographic heterogeneity? and do these processes matter at all in determining tree species distribution and abundance? We have studied these questions with the Neotropical Symphonia tree species. In the Guiana shield, Symphonia are represented by at least two sympatric taxa or ecotypes, Symphonia globulifera found almost exclusively in bottomlands, and a yet undescribed more generalist taxon/ecotype, Symphonia sp1. A reciprocal transplantation experiment (510 seeds, 16 conditions) was set up and followed over the course of 6 years to evaluate the survival and performance of individuals from different ecotypes and provenances. Germination, survival, growth, and herbivory showed signs of local adaptation, with some combinations of ecotypes and provenances growing faster and surviving better in their own habitat or provenance region. S. globulifera was strongly penalized when planted outside its home habitat but showed the fastest growth rates when planted in its home habitat, suggesting it is a specialist of a high-risk high-gain strategy. Conversely, S. sp1 behaved as a generalist, performing well in a variety of environments. The differential performance of seeds and seedlings in the different habitats matches the known distribution of both ecotypes, indicating that environmental filtering at the very early stages can be a key determinant of tree species distributions, even at the microgeographic level and among very closely related taxa. Furthermore, such differential performance also contributes to explain, in part, the maintenance of the different Symphonia ecotypes living in intimate sympatry despite occasional gene flow.

Genetic Structure of a Worldwide Germplasm Collection of Prunus armeniaca L. Reveals Three Major Diffusion Routes for Varieties Coming From the Species' Center of Origin.

The characterization of the largest worldwide representative data set of apricot (Prunus armeniaca L.) germplasm was performed using molecular markers. Genetic diversity and structure of the cultivated apricot genetic resources were analyzed to decipher the history of diffusion of this species around the world. A common set of 25 microsatellite markers was used for genotyping a total of 890 apricot accessions in different collections from the center of origin to the more recent regions of apricot culture. Using a Bayesian model-based clustering approach, the apricot genotypes can be structured into five different genetic clusters (FST = 0.174), correlated with the geographical regions of origin of the accessions. Accessions from China and Central Asia were clustered together and exhibited the highest levels of diversity, confirming an origin in this region. A loss of genetic diversity was observed from the center of origin to both western and eastern zones of recent apricot culture. Altogether, our results revealed that apricot spread from China and Central Asia, defined as the center of origin, following three major diffusion routes with a decreasing gradient of genetic variation in each geographical group. The identification of specific alleles outside the center of origin confirmed the existence of different secondary apricot diversification centers. The present work provides more understanding of the worldwide history of apricot species diffusion as well as the field of conservation of the available genetic resources. Data have been used to define an apricot core collection based on molecular marker diversity which will be useful for further identification of genomic regions associated with commercially important horticultural traits through genome-wide association studies to sustain apricot breeding programs.

Genetic Variation and Differentiation of Hylesia metabus (Lepidoptera: Saturniidae): Moths of Public Health Importance in French Guiana and in Venezuela.

Hylesia moths impact human health in South America, inducing epidemic outbreaks of lepidopterism, a puriginous dermatitis caused by the urticating properties of females' abdominal setae. The classification of the Hylesia genus is complex, owing to its high diversity in Amazonia, high intraspecific morphological variance, and lack of interspecific diagnostic traits which may hide cryptic species. Outbreaks of Hylesia metabus have been considered responsible for the intense outbreaks of lepidopterism in Venezuela and French Guiana since the C20, however, little is known about genetic variability throughout the species range, which is instrumental for establishing control strategies on H. metabus. Seven microsatellites and mitochondrial gene markers were analyzed from Hylesia moths collected from two major lepidopterism outbreak South American regions. The mitochondrial gene sequences contained significant genetic variation, revealing a single, widespread, polymorphic species with distinct clusters, possibly corresponding to populations evolving in isolation. The microsatellite markers validated the mitochondrial results, and suggest the presence of three populations: one in Venezuela, and two in French Guiana. All moths sampled during outbreak events in French Guiana were assigned to a single coastal population. The causes and implications of this finding require further research.

Chloroplast DNA variation in a hyperdiverse tropical tree community.

We investigate chloroplast DNA variation in a hyperdiverse community of tropical rainforest trees in French Guiana, focusing on patterns of intraspecific and interspecific variation. We test whether a species genetic diversity is higher when it has congeners in the community with which it can exchange genes and if shared haplotypes are more frequent in genetically diverse species, as expected in the presence of introgression.We sampled a total of 1,681 individual trees from 472 species corresponding to 198 genera and sequenced them at a noncoding chloroplast DNA fragment.Polymorphism was more frequent in species that have congeneric species in the study site than in those without congeners (30% vs. 12%). Moreover, more chloroplast haplotypes were shared with congeners in polymorphic species than in monomorphic ones (44% vs. 28%).Despite large heterogeneities caused by genus-specific behaviors in patterns of hybridization, these results suggest that the higher polymorphism in the presence of congeners is caused by local introgression rather than by incomplete lineage sorting. Our findings suggest that introgression has the potential to drive intraspecific genetic diversity in species-rich tropical forests.