Stéphane Joost

A spatial analysis method (SAM) to detect candidate loci for selection: towards a landscape genomics approach to adaptation

The detection of adaptive loci in the genome is essential as it gives the possibility of understanding what proportion of a genome or which genes are being shaped by natural selection. Several statistical methods have been developed which make use of molecular data to reveal genomic regions under selection. In this paper, we propose an approach to address this issue from the environmental angle, in order to complement results obtained by population genetics. We introduce a new method to detect signatures of natural selection based on the application of spatial analysis, with the contribution of geographical information systems (GIS), environmental variables and molecular data. Multiple univariate logistic regressions were carried out to test for association between allelic frequencies at marker loci and environmental variables. This spatial analysis method (SAM) is similar to current population genomics approaches since it is designed to scan hundreds of markers to assess a putative association with hundreds of environmental variables. Here, by application to studies of pine weevils and breeds of sheep we demonstrate a strong correspondence between SAM results and those obtained using population genetics approaches. Statistical signals were found that associate loci with environmental parameters, and these loci behave atypically in comparison with the theoretical distribution for neutral loci. The contribution of this new tool is not only to permit the identification of loci under selection but also to establish hypotheses about ecological factors that could exert the selection pressure responsible. In the future, such an approach may accelerate the process of hunting for functional genes at the population level.

Perspectives on the use of landscape genetics to detect genetic adaptive variation in the field

Understanding the genetic basis of species adaptation in the context of global change poses one of the greatest challenges of this century. Although we have begun to understand the molecular basis of adaptation in those species for which whole genome sequences are available, the molecular basis of adaptation is still poorly understood for most non‐model species. In this paper, we outline major challenges and future research directions for correlating environmental factors with molecular markers to identify adaptive genetic variation, and point to research gaps in the application of landscape genetics to real‐world problems arising from global change, such as the ability of organisms to adapt over rapid time scales. High throughput sequencing generates vast quantities of molecular data to address the challenge of studying adaptive genetic variation in non‐model species. Here, we suggest that improvements in the sampling design should consider spatial dependence among sampled individuals. Then, we describe available statistical approaches for integrating spatial dependence into landscape analyses of adaptive genetic variation.

Objectives, criteria and methods for using molecular genetic data in priority setting for conservation of animal genetic resources

The genetic diversity of the worlds livestock populations is decreasing, both within and across breeds. A wide variety of factors has contributed to the loss, replacement or genetic dilution of many local breeds. Genetic variability within the more common commercial breeds has been greatly decreased by selectively intense breeding programmes. Conservation of livestock genetic variability is thus important, especially when considering possible future changes in production environments. The world has more than 7500 livestock breeds and conservation of all of them is not feasible. Therefore, prioritization is needed. The objective of this article is to review the state of the art in approaches for prioritization of breeds for conservation, particularly those approaches that consider molecular genetic information, and to identify any shortcomings that may restrict their application. The Weitzman method was among the first and most well-known approaches for utilization of molecular genetic information in conservation prioritization. This approach balances diversity and extinction probability to yield an objective measure of conservation potential. However, this approach was designed for decision making across species and measures diversity as distinctiveness. For livestock, prioritization will most commonly be performed among breeds within species, so alternatives that measure diversity as co-ancestry (i.e. also within-breed variability) have been proposed. Although these methods are technically sound, their application has generally been limited to research studies; most existing conservation programmes have effectively primarily based decisions on extinction risk. The development of user-friendly software incorporating these approaches may increase their rate of utilization.

Adaptations to Climate-Mediated Selective Pressures in Sheep

Following domestication, sheep (Ovis aries) have become essential farmed animals across the world through adaptation to a diverse range of environments and varied production systems. Climate-mediated selective pressure has shaped phenotypic variation and has left genetic “footprints” in the genome of breeds raised in different agroecological zones. Unlike numerous studies that have searched for evidence of selection using only population genetics data, here, we conducted an integrated coanalysis of environmental data with single nucleotide polymorphism (SNP) variation. By examining 49,034 SNPs from 32 old, autochthonous sheep breeds that are adapted to a spectrum of different regional climates, we identified 230 SNPs with evidence for selection that is likely due to climate-mediated pressure. Among them, 189 (82%) showed significant correlation (P ≤ 0.05) between allele frequency and climatic variables in a larger set of native populations from a worldwide range of geographic areas and climates. Gene ontology analysis of genes colocated with significant SNPs identified 17 candidates related to GTPase regulator and peptide receptor activities in the biological processes of energy metabolism and endocrine and autoimmune regulation. We also observed high linkage disequilibrium and significant extended haplotype homozygosity for the core haplotype TBC1D12-CH1 of TBC1D12. The global frequency distribution of the core haplotype and allele OAR22_18929579-A showed an apparent geographic pattern and significant (P ≤ 0.05) correlations with climatic variation. Our results imply that adaptations to local climates have shaped the spatial distribution of some variants that are candidates to underpin adaptive variation in sheep.

Landscape genomics and biased FST approaches reveal single nucleotide polymorphisms under selection in goat breeds of North-East Mediterranean

BackgroundIn this study we compare outlier loci detected using a FST based method with those identified by a recently described method based on spatial analysis (SAM). We tested a panel of single nucleotide polymorphisms (SNPs) previously genotyped in individuals of goat breeds of southern areas of the Mediterranean basin (Italy, Greece and Albania). We evaluate how the SAM method performs with SNPs, which are increasingly employed due to their high number, low cost and easy of scoring.ResultsThe combined use of the two outlier detection approaches, never tested before using SNP polymorphisms, resulted in the identification of the same three loci involved in milk and meat quality data by using the two methods, while the FST based method identified 3 more loci as under selection sweep in the breeds examined.ConclusionData appear congruent by using the two methods for FST values exceeding the 99% confidence limits. The methods of FST and SAM can independently detect signatures of selection and therefore can reduce the probability of finding false positives if employed together. The outlier loci identified in this study could indicate adaptive variation in the analysed species, characterized by a large range of climatic conditions in the rearing areas and by a history of intense trade, that implies plasticity in adapting to new environments.

Integrating landscape genomics and spatially explicit approaches to detect loci under selection in clinal populations

Uncovering the genetic basis of adaptation hinges on the ability to detect loci under selection. However, population genomics outlier approaches to detect selected loci may be inappropriate for clinal populations or those with unclear population structure because they require that individuals be clustered into populations. An alternate approach, landscape genomics, uses individual‐based approaches to detect loci under selection and reveal potential environmental drivers of selection. We tested four landscape genomics methods on a simulated clinal population to determine their effectiveness at identifying a locus under varying selection strengths along an environmental gradient. We found all methods produced very low type I error rates across all selection strengths, but elevated type II error rates under “weak” selection. We then applied these methods to an AFLP genome scan of an alpine plant, Campanula barbata, and identified five highly supported candidate loci associated with precipitation variables. These loci also showed spatial autocorrelation and cline patterns indicative of selection along a precipitation gradient. Our results suggest that landscape genomics in combination with other spatial analyses provides a powerful approach for identifying loci potentially under selection and explaining spatially complex interactions between species and their environment.

Detecting spatial genetic signatures of local adaptation in heterogeneous landscapes

The spatial structure of the environment (e.g. the configuration of habitat patches) may play an important role in determining the strength of local adaptation. However, previous studies of habitat heterogeneity and local adaptation have largely been limited to simple landscapes, which poorly represent the multiscale habitat structure common in nature. Here, we use simulations to pursue two goals: (i) we explore how landscape heterogeneity, dispersal ability and selection affect the strength of local adaptation, and (ii) we evaluate the performance of several genotype–environment association (GEA) methods for detecting loci involved in local adaptation. We found that the strength of local adaptation increased in spatially aggregated selection regimes, but remained strong in patchy landscapes when selection was moderate to strong. Weak selection resulted in weak local adaptation that was relatively unaffected by landscape heterogeneity. In general, the power of detection methods closely reflected levels of local adaptation. False‐positive rates (FPRs), however, showed distinct differences across GEA methods based on levels of population structure. The univariate GEA approach had high FPRs (up to 55%) under limited dispersal scenarios, due to strong isolation by distance. By contrast, multivariate, ordination‐based methods had uniformly low FPRs (0–2%), suggesting these approaches can effectively control for population structure. Specifically, constrained ordinations had the best balance of high detection and low FPRs and will be a useful addition to the GEA toolkit. Our results provide both theoretical and practical insights into the conditions that shape local adaptation and how these conditions impact our ability to detect selection.

Uncovering the genetic basis of adaptive change: on the intersection of landscape genomics and theoretical population genetics

A workshop recently held at the École Polytechnique Fédérale de Lausanne (EPFL, Switzerland) was dedicated to understanding the genetic basis of adaptive change, taking stock of the different approaches developed in theoretical population genetics and landscape genomics and bringing together knowledge accumulated in both research fields. Indeed, an important challenge in theoretical population genetics is to incorporate effects of demographic history and population structure. But important design problems (e.g. focus on populations as units, focus on hard selective sweeps, no hypothesis‐based framework in the design of the statistical tests) reduce their capability of detecting adaptive genetic variation. In parallel, landscape genomics offers a solution to several of these problems and provides a number of advantages (e.g. fast computation, landscape heterogeneity integration). But the approach makes several implicit assumptions that should be carefully considered (e.g. selection has had enough time to create a functional relationship between the allele distribution and the environmental variable, or this functional relationship is assumed to be constant). To address the respective strengths and weaknesses mentioned above, the workshop brought together a panel of experts from both disciplines to present their work and discuss the relevance of combining these approaches, possibly resulting in a joint software solution in the future.

Gene–obesogenic environment interactions in the UK Biobank study

Abstract Background: Previous studies have suggested that modern obesogenic environments accentuate the genetic risk of obesity. However, these studies have proven controversial as to which, if any, measures of the environment accentuate genetic susceptibility to high body mass index (BMI). Methods: We used up to 120 000 adults from the UK Biobank study to test the hypothesis that high-risk obesogenic environments and behaviours accentuate genetic susceptibility to obesity. We used BMI as the outcome and a 69-variant genetic risk score (GRS) for obesity and 12 measures of the obesogenic environment as exposures. These measures included Townsend deprivation index (TDI) as a measure of socio-economic position, TV watching, a ‘Westernized’ diet and physical activity. We performed several negative control tests, including randomly selecting groups of different average BMIs, using a simulated environment and including sun-protection use as an environment. Results: We found gene–environment interactions with TDI (Pinteraction = 3 × 10–10), self-reported TV watching (Pinteraction = 7 × 10–5) and self-reported physical activity (Pinteraction = 5 × 10–6). Within the group of 50% living in the most relatively deprived situations, carrying 10 additional BMI-raising alleles was associated with approximately 3.8 kg extra weight in someone 1.73 m tall. In contrast, within the group of 50% living in the least deprivation, carrying 10 additional BMI-raising alleles was associated with approximately 2.9 kg extra weight. The interactions were weaker, but present, with the negative controls, including sun-protection use, indicating that residual confounding is likely. Conclusions: Our findings suggest that the obesogenic environment accentuates the risk of obesity in genetically susceptible adults. Of the factors we tested, relative social deprivation best captures the aspects of the obesogenic environment responsible.

High performance computation of landscape genomic models including local indicators of spatial association.

With the increasing availability of both molecular and topo‐climatic data, the main challenges facing landscape genomics – that is the combination of landscape ecology with population genomics – include processing large numbers of models and distinguishing between selection and demographic processes (e.g. population structure). Several methods address the latter, either by estimating a null model of population history or by simultaneously inferring environmental and demographic effects. Here we present samβada, an approach designed to study signatures of local adaptation, with special emphasis on high performance computing of large‐scale genetic and environmental data sets. samβada identifies candidate loci using genotype–environment associations while also incorporating multivariate analyses to assess the effect of many environmental predictor variables. This enables the inclusion of explanatory variables representing population structure into the models to lower the occurrences of spurious genotype–environment associations. In addition, samβada calculates local indicators of spatial association for candidate loci to provide information on whether similar genotypes tend to cluster in space, which constitutes a useful indication of the possible kinship between individuals. To test the usefulness of this approach, we carried out a simulation study and analysed a data set from Ugandan cattle to detect signatures of local adaptation with samβada, bayenv, lfmm and an FST outlier method (FDIST approach in arlequin) and compare their results. samβada – an open source software for Windows, Linux and Mac OS X available at http://lasig.epfl.ch/sambada – outperforms other approaches and better suits whole‐genome sequence data processing.