A conceptual framework for the spatial analysis of landscape genetic data

Targeting environmental adaptation in the monocot model Brachypodium distachyon: a multi-faceted approach

BACKGROUND

The local environment plays a major role in the spatial distribution of plant populations. Natural plant populations have an extremely poor displacing capacity, so their continued survival in a given environment depends on how well they adapt to local pedoclimatic conditions. Genomic tools can be used to identify adaptive traits at a DNA level and to further our understanding of evolutionary processes. Here we report the use of genotyping-by-sequencing on local groups of the sequenced monocot model species Brachypodium distachyon. Exploiting population genetics, landscape genomics and genome wide association studies, we evaluate B. distachyon role as a natural probe for identifying genomic loci involved in environmental adaptation.

RESULTS

Brachypodium distachyon individuals were sampled in nine locations with different ecologies and characterized with 16,697 SNPs. Variations in sequencing depth showed consistent patterns at 8,072 genomic bins, which were significantly enriched in transposable elements. We investigated the structuration and diversity of this collection, and exploited climatic data to identify loci with adaptive significance through i) two different approaches for genome wide association analyses considering climatic variation, ii) an outlier loci approach, and iii) a canonical correlation analysis on differentially sequenced bins. A linkage disequilibrium-corrected Bonferroni method was applied to filter associations. The two association methods jointly identified a set of 15 genes significantly related to environmental adaptation. The outlier loci approach revealed that 5.7% of the loci analysed were under selection. The canonical correlation analysis showed that the distribution of some differentially sequenced regions was associated to environmental variation.

CONCLUSIONS

We show that the multi-faceted approach used here targeted different components of B. distachyon adaptive variation, and may lead to the discovery of genes related to environmental adaptation in natural populations. Its application to a model species with a fully sequenced genome is a modular strategy that enables the stratification of biological material and thus improves our knowledge of the functional loci determining adaptation in near-crop species. When coupled with population genetics and measures of genomic structuration, methods coming from genome wide association studies may lead to the exploitation of model species as natural probes to identify loci related to environmental adaptation.

Local topography shapes fine-scale spatial genetic structure in the Arkansas Valley evening primrose, Oenothera harringtonii (Onagraceae)

Identifying factors that shape the spatial distribution of genetic variation is crucial to understanding many population- and landscape-level processes. In this study, we explore fine-scale spatial genetic structure in Oenothera harringtonii (Onagraceae), an insect-pollinated, gravity-dispersed herb endemic to the grasslands of south-central and southeastern Colorado, USA. We genotyped 315 individuals with 11 microsatellite markers and utilized a combination of spatial autocorrelation analyses and landscape genetic models to relate life history traits and landscape features to dispersal processes. Spatial genetic structure was consistent with theoretical expectations of isolation by distance, but this pattern was weak (Sp = 0.00374). Anisotropic analyses indicated that spatial genetic structure was markedly directional, in this case consistent with increased dispersal along prominent slopes. Landscape genetic models subsequently confirmed that spatial genetic variation was significantly influenced by local topographic heterogeneity, specifically that geographic distance, elevation and aspect were important predictors of spatial genetic structure. Among these variables, geographic distance was ~68% more important than elevation in describing spatial genetic variation, and elevation was ~42% more important than aspect after removing the effect of geographic distance. From these results, we infer a mechanism of hydrochorous seed dispersal along major drainages aided by seasonal monsoon rains. Our findings suggest that landscape features may shape microevolutionary processes at much finer spatial scales than typically considered, and stress the importance of considering how particular dispersal vectors are influenced by their environmental context.

Modelling the dispersal of the two main hosts of the raccoon rabies variant in heterogeneous environments with landscape genetics

Predicting the geographic spread of wildlife epidemics requires knowledge about the movement patterns of disease hosts or vectors. The field of landscape genetics provides valuable approaches to study dispersal indirectly, which in turn may be used to understand patterns of disease spread. Here, we applied landscape genetic analyses and spatially explicit models to identify the potential path of raccoon rabies spread in a mesocarnivore community. We used relatedness estimates derived from microsatellite genotypes of raccoons and striped skunks to investigate their dispersal patterns in a heterogeneous landscape composed predominantly of agricultural, forested and residential areas. Samples were collected in an area covering 22 000 km² in southern Québec, where the raccoon rabies variant (RRV) was first detected in 2006. Multiple regressions on distance matrices revealed that genetic distance among male raccoons was strictly a function of geographic distance, while dispersal in female raccoons was significantly reduced by the presence of agricultural fields. In skunks, our results suggested that dispersal is increased in edge habitats between fields and forest fragments in both males and females. Resistance modelling allowed us to identify likely dispersal corridors used by these two rabies hosts, which may prove especially helpful for surveillance and control (e.g. oral vaccination) activities.

Climate-Driven Reshuffling of Species and Genes: Potential Conservation Roles for Species Translocations and Recombinant Hybrid Genotypes

Comprising 50%-75% of the world's fauna, insects are a prominent part of biodiversity in communities and ecosystems globally. Biodiversity across all levels of biological classifications is fundamentally based on genetic diversity. However, the integration of genomics and phylogenetics into conservation management may not be as rapid as climate change. The genetics of hybrid introgression as a source of novel variation for ecological divergence and evolutionary speciation (and resilience) may generate adaptive potential and diversity fast enough to respond to locally-altered environmental conditions. Major plant and herbivore hybrid zones with associated communities deserve conservation consideration. This review addresses functional genetics across multi-trophic-level interactions including "invasive species" in various ecosystems as they may become disrupted in different ways by rapid climate change. "Invasive genes" (into new species and populations) need to be recognized for their positive creative potential and addressed in conservation programs. "Genetic rescue" via hybrid translocations may provide needed adaptive flexibility for rapid adaptation to environmental change. While concerns persist for some conservationists, this review emphasizes the positive aspects of hybrids and hybridization. Specific implications of natural genetic introgression are addressed with a few examples from butterflies, including transgressive phenotypes and climate-driven homoploid recombinant hybrid speciation. Some specific examples illustrate these points using the swallowtail butterflies (Papilionidae) with their long-term historical data base (phylogeographical diversity changes) and recent (3-decade) climate-driven temporal and genetic divergence in recombinant homoploid hybrids and relatively recent hybrid speciation of Papilio appalachiensis in North America. Climate-induced "reshuffling" (recombinations) of species composition, genotypes, and genomes may become increasingly ecologically and evolutionarily predictable, but future conservation management programs are more likely to remain constrained by human behavior than by lack of academic knowledge.

Mantel test in population genetics

The comparison of genetic divergence or genetic distances, estimated by pairwise FST and related statistics, with geographical distances by Mantel test is one of the most popular approaches to evaluate spatial processes driving population structure. There have been, however, recent criticisms and discussions on the statistical performance of the Mantel test. Simultaneously, alternative frameworks for data analyses are being proposed. Here, we review the Mantel test and its variations, including Mantel correlograms and partial correlations and regressions. For illustrative purposes, we studied spatial genetic divergence among 25 populations of Dipteryx alata ("Baru"), a tree species endemic to the Cerrado, the Brazilian savannas, based on 8 microsatellite loci. We also applied alternative methods to analyze spatial patterns in this dataset, especially a multivariate generalization of Spatial Eigenfunction Analysis based on redundancy analysis. The different approaches resulted in similar estimates of the magnitude of spatial structure in the genetic data. Furthermore, the results were expected based on previous knowledge of the ecological and evolutionary processes underlying genetic variation in this species. Our review shows that a careful application and interpretation of Mantel tests, especially Mantel correlograms, can overcome some potential statistical problems and provide a simple and useful tool for multivariate analysis of spatial patterns of genetic divergence.