Neutral and adaptive drivers of microgeographic genetic divergence within continuous populations: the case of the neotropical tree Eperua falcata (Aubl.)

Environmental selection, rather than neutral processes, best explain regional patterns of diversity in a tropical rainforest fish

To conserve the high functional and genetic variation in hotspots such as tropical rainforests, it is essential to understand the forces driving and maintaining biodiversity. We asked to what extent environmental gradients and terrain structure affect morphological and genomic variation across the wet tropical distribution of an Australian rainbowfish, Melanotaenia splendida splendida. We used an integrative riverscape genomics and morphometrics framework to assess the influence of these factors on both putative adaptive and non-adaptive spatial divergence. We found that neutral genetic population structure was largely explainable by restricted gene flow among drainages. However, environmental associations revealed that ecological variables had a similar power to explain overall genetic variation, and greater power to explain body shape variation, than the included neutral covariables. Hydrological and thermal variables were the strongest environmental predictors and were correlated with traits previously linked to heritable habitat-associated dimorphism in rainbowfishes. In addition, climate-associated genetic variation was significantly associated with morphology, supporting heritability of shape variation. These results support the inference of evolved functional differences among localities, and the importance of hydroclimate in early stages of diversification. We expect that substantial evolutionary responses will be required in tropical rainforest endemics to mitigate local fitness losses due to changing climates.

Topography drives microgeographic adaptations of closely related species in two tropical tree species complexes

Closely related tree species that grow in sympatry are abundant in rainforests. However, little is known of the eco-evolutionary processes that govern their niches and local coexistence. We assessed genetic species delimitation in closely related sympatric species belonging to two Neotropical tree species complexes and investigated their genomic adaptation to a fine-scale topographic gradient with associated edaphic and hydrologic features. Combining LiDAR-derived topography, tree inventories, and single nucleotide polymorphisms (SNPs) from gene capture experiments, we explored genome-wide population genetic structure, covariation of environmental variables, and genotype-environment association to assess microgeographic adaptations to topography within the species complexes Symphonia (Clusiaceae), and Eschweilera (Lecythidaceae) with three species per complex and 385 and 257 individuals genotyped, respectively. Within species complexes, closely related tree species had different realized optima for topographic niches defined through the topographic wetness index or the relative elevation, and species displayed genetic signatures of adaptations to these niches. Symphonia species were genetically differentiated along water and nutrient distribution particularly in genes responding to water deprivation, whereas Eschweilera species were genetically differentiated according to soil chemistry. Our results suggest that varied topography represents a powerful driver of processes modulating tropical forest biodiversity with differential adaptations that stabilize local coexistence of closely related tree species.

Mechanisms driving plant functional trait variation in a tropical forest

Plant functional trait variation in tropical forests results from taxonomic differences in phylogeny and associated genetic differences, as well as, phenotypic plastic responses to the environment. Accounting for the underlying mechanisms driving plant functional trait variation is important for understanding the potential rate of change of ecosystems since trait acclimation via phenotypic plasticity is very fast compared to shifts in community composition and genetic adaptation. We here applied a statistical technique to decompose the relative roles of phenotypic plasticity, genetic adaptation, and phylogenetic constraints. We examined typically obtained plant functional traits, such as wood density, plant height, specific leaf area, leaf area, leaf thickness, leaf dry mass content, leaf nitrogen content, and leaf phosphorus content. We assumed that genetic differences in plant functional traits between species and genotypes increase with environmental heterogeneity and geographic distance, whereas trait variation due to plastic acclimation to the local environment is independent of spatial distance between sampling sites. Results suggest that most of the observed trait variation could not be explained by the measured environmental variables, thus indicating a limited potential to predict individual plant traits from commonly assessed parameters. However, we found a difference in the response of plant functional traits, such that leaf traits varied in response to canopy-light regime and nutrient availability, whereas wood traits were related to topoedaphic factors and water availability. Our analysis furthermore revealed differences in the functional response of coexisting neotropical tree species, which suggests that endemic species with conservative ecological strategies might be especially prone to competitive exclusion under projected climate change.

Ecological Factors Generally Not Altitude Related Played Main Roles in Driving Potential Adaptive Evolution at Elevational Range Margin Populations of Taiwan Incense Cedar (Calocedrus formosana)

Population diversification can be shaped by a combination of environmental factors as well as geographic isolation interacting with gene flow. We surveyed genetic variation of 243 samples from 12 populations of Calocedrus formosana using amplified fragment length polymorphism (AFLP) and scored a total of 437 AFLP fragments using 11 selective amplification primer pairs. The AFLP variation was used to assess the role of gene flow on the pattern of genetic diversity and to test environments in driving population adaptive evolution. This study found the relatively lower level of genetic diversity and the higher level of population differentiation in C. formosana compared with those estimated in previous studies of conifers including Cunninghamia konishii, Keteleeria davidiana var. formosana, and Taiwania cryptomerioides occurring in Taiwan. BAYESCAN detected 26 F_ST outlier loci that were found to be associated strongly with various environmental variables using multiple univariate logistic regression, latent factor mixed model, and Bayesian logistic regression. We found several environmentally dependent adaptive loci with high frequencies in low- or high-elevation populations, suggesting their involvement in local adaptation. Ecological factors, including relative humidity and sunshine hours, that are generally not altitude related could have been the most important selective drivers for population divergent evolution in C. formosana. The present study provides fundamental information in relation to adaptive evolution and can be useful for assisted migration program of C. formosana in the future conservation of this species.

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.

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.

Looking for Local Adaptation: Convergent Microevolution in Aleppo Pine (Pinus halepensis)

Finding outlier loci underlying local adaptation is challenging and is best approached by suitable sampling design and rigorous method selection. In this study, we aimed to detect outlier loci (single nucleotide polymorphisms, SNPs) at the local scale by using Aleppo pine (Pinus halepensis), a drought resistant conifer that has colonized many habitats in the Mediterranean Basin, as the model species. We used a nested sampling approach that considered replicated altitudinal gradients for three contrasting sites. We genotyped samples at 294 SNPs located in genomic regions selected to maximize outlier detection. We then applied three different statistical methodologies-Two Bayesian outlier methods and one latent factor principal component method-To identify outlier loci. No SNP was an outlier for all three methods, while eight SNPs were detected by at least two methods and 17 were detected only by one method. From the intersection of outlier SNPs, only one presented an allelic frequency pattern associated with the elevational gradient across the three sites. In a context of multiple populations under similar selective pressures, our results underline the need for careful examination of outliers detected in genomic scans before considering them as candidates for convergent adaptation.

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.

Adaptive Genetic Divergence Despite Significant Isolation-by-Distance in Populations of Taiwan Cow-Tail Fir (Keteleeria davidiana var. formosana)

Double digest restriction site-associated DNA sequencing (ddRADseq) is a tool for delivering genome-wide single nucleotide polymorphism (SNP) markers for non-model organisms useful in resolving fine-scale population structure and detecting signatures of selection. This study performs population genetic analysis, based on ddRADseq data, of a coniferous species, Keteleeria davidiana var. formosana, disjunctly distributed in northern and southern Taiwan, for investigation of population adaptive divergence in response to environmental heterogeneity. A total of 13,914 SNPs were detected and used to assess genetic diversity, FST outlier detection, population genetic structure, and individual assignments of five populations (62 individuals) of K. davidiana var. formosana. Principal component analysis (PCA), individual assignments, and the neighbor-joining tree were successful in differentiating individuals between northern and southern populations of K. davidiana var. formosana, but apparent gene flow between the southern DW30 population and northern populations was also revealed. Fifteen of 23 highly differentiated SNPs identified were found to be strongly associated with environmental variables, suggesting isolation-by-environment (IBE). However, multiple matrix regression with randomization analysis revealed strong IBE as well as significant isolation-by-distance. Environmental impacts on divergence were found between populations of the North and South regions and also between the two southern neighboring populations. BLASTN annotation of the sequences flanking outlier SNPs gave significant hits for three of 23 markers that might have biological relevance to mitochondrial homeostasis involved in the survival of locally adapted lineages. Species delimitation between K. davidiana var. formosana and its ancestor, K. davidiana, was also examined (72 individuals). This study has produced highly informative population genomic data for the understanding of population attributes, such as diversity, connectivity, and adaptive divergence associated with large- and small-scale environmental heterogeneity in K. davidiana var. formosana.

Altitudinal gradients, biogeographic history and microhabitat adaptation affect fine-scale spatial genetic structure in African and Neotropical populations of an ancient tropical tree species

The analysis of fine-scale spatial genetic structure (FSGS) within populations can provide insights into eco-evolutionary processes. Restricted dispersal and locally occurring genetic drift are the primary causes for FSGS at equilibrium, as described in the isolation by distance (IBD) model. Beyond IBD expectations, spatial, environmental or historical factors can affect FSGS. We examined FSGS in seven African and Neotropical populations of the late-successional rain forest tree Symphonia globulifera L. f. (Clusiaceae) to discriminate the influence of drift-dispersal vs. landscape/ecological features and historical processes on FSGS. We used spatial principal component analysis and Bayesian clustering to assess spatial genetic heterogeneity at SSRs and examined its association with plastid DNA and habitat features. African populations (from Cameroon and São Tomé) displayed a stronger FSGS than Neotropical populations at both marker types (mean Sp = 0.025 vs. Sp = 0.008 at SSRs) and had a stronger spatial genetic heterogeneity. All three African populations occurred in pronounced altitudinal gradients, possibly restricting animal-mediated seed dispersal. Cyto-nuclear disequilibria in Cameroonian populations also suggested a legacy of biogeographic history to explain these genetic patterns. Conversely, Neotropical populations exhibited a weaker FSGS, which may reflect more efficient wide-ranging seed dispersal by Neotropical bats and other dispersers. The population from French Guiana displayed an association of plastid haplotypes with two morphotypes characterized by differential habitat preferences. Our results highlight the importance of the microenvironment for eco-evolutionary processes within persistent tropical tree populations.

Postglacial range expansion and the role of ecological factors in driving adaptive evolution of Musa basjoo var. formosana

Genetic variation evolves during postglacial range expansion of a species and is important for adapting to varied environmental conditions. It is crucial for the future survival of a species. We investigate the nuclear DNA sequence variation to provide evidence of postglacial range expansion of Musa basjoo var. formosana, a wild banana species, and test for adaptive evolution of amplified fragment length polymorphic (AFLP) loci underlying local adaptation in association with environmental variables. Postglacial range expansion was suggested by phylogeographical analyses based on sequence variation of the second intron of copper zinc superoxide dismutase 2 gene. Two glacial refugia were inferred by the average F_ST parameter (mean F_ST of a population against the remaining populations). Using variation partitioning by redundancy analysis, we found a significant amount of explained AFLP variation attributed to environmental and spatially-structured environmental effects. By combining genome scan methods and multiple univariate logistic regression, four AFLP loci were found to be strongly associated with environmental variables, including temperature, precipitation, soil moisture, wet days, and surface coverage activity representing vegetation greenness. These environmental variables may have played various roles as ecological drivers for adaptive evolution of M. basjoo var. formosana during range expansion after the last glacial maximum.

Pollen-mediated gene flow and fine-scale spatial genetic structure in Olea europaea subsp. europaea var. sylvestris

BACKGROUND AND AIMS

Wild olive ( Olea europaea subsp. europaea var. sylvestris ) is important from an economic and ecological point of view. The effects of anthropogenic activities may lead to the genetic erosion of its genetic patrimony, which has high value for breeding programmes. In particular, the consequences of the introgression from cultivated stands are strongly dependent on the extent of gene flow and therefore this work aims at quantitatively describing contemporary gene flow patterns in wild olive natural populations.

METHODS

The studied wild population is located in an undisturbed forest, in southern Spain, considered one of the few extant hotspots of true oleaster diversity. A total of 225 potential father trees and seeds issued from five mother trees were genotyped by eight microsatellite markers. Levels of contemporary pollen flow, in terms of both pollen immigration rates and within-population dynamics, were measured through paternity analyses. Moreover, the extent of fine-scale spatial genetic structure (SGS) was studied to assess the relative importance of seed and pollen dispersal in shaping the spatial distribution of genetic variation.

KEY RESULTS

The results showed that the population under study is characterized by a high genetic diversity, a relatively high pollen immigration rate (0·57), an average within-population pollen dispersal of about 107 m and weak but significant SGS up to 40 m. The population is a mosaic of several intermingled genetic clusters that is likely to be generated by spatially restricted seed dispersal. Moreover, wild oleasters were found to be self-incompatible and preferential mating between some genotypes was revealed.

CONCLUSIONS

Knowledge of the within-population genetic structure and gene flow dynamics will lead to identifying possible strategies aimed at limiting the effect of anthropogenic activities and improving breeding programmes for the conservation of olive tree forest genetic resources.