Combining niche modelling and landscape genetics to study local adaptation: A novel approach illustrated using alpine plants

Understanding the factors that shape adaptive genetic variation across species niches has become of paramount importance in evolutionary ecology, especially to understand how adaptation to changing climate affects the geographic range of species. The distribution of adaptive alleles in the ecological niche is determined by the emergence of novel mutations, their fitness consequences and gene flow that connects populations across species niches. Striking demographical differences and source sink dynamics of populations between the centre and the margin of the niche can play a major role in the emergence and spread of adaptive alleles. Although some theoretical predictions have long been proposed, the origin and distribution of adaptive alleles within species niches remain untested. In this paper, we propose and discuss a novel empirical approach that combines landscape genetics with species niche modelling, to test whether alleles that confer local adaptation are more likely to occur in either marginal or central populations of species niches. We illustrate this new approach by using a published data set of 21 alpine plant species genotyped with a total of 2483 amplified fragment length polymorphisms (AFLP), distributed over more than 1733 sampling sites across the Alps. Based on the assumption that alleles that were statistically associated with environmental variables were adaptive, we found that adaptive alleles in the margin of a species niche were also present in the niche centre, which suggests that adaptation originates in the niche centre. These findings corroborate models of species range evolution, in which the centre of the niche contributes to the emergence of novel adaptive alleles, which diffuse towards niche margins and facilitate niche and range expansion through subsequent local adaptation. Although these results need to be confirmed via fitness measurements in natural populations and functionally characterised genetic sequences, this study provides a first step towards understanding how adaptive genetic variation emerges and shapes species niches and geographic ranges along environmental gradients.

Broad-scale adaptive genetic variation in alpine plants is driven by temperature and precipitation.

Identifying adaptive genetic variation is a challenging task, in particular in non-model species for which genomic information is still limited or absent. Here, we studied distribution patterns of amplified fragment length polymorphisms (AFLPs) in response to environmental variation, in 13 alpine plant species consistently sampled across the entire European Alps. Multiple linear regressions were performed between AFLP allele frequencies per site as dependent variables and two categories of independent variables, namely Moran's eigenvector map MEM variables (to account for spatial and unaccounted environmental variation, and historical demographic processes) and environmental variables. These associations allowed the identification of 153 loci of ecological relevance. Univariate regressions between allele frequency and each environmental factor further showed that loci of ecological relevance were mainly correlated with MEM variables. We found that precipitation and temperature were the best environmental predictors, whereas topographic factors were rarely involved in environmental associations. Climatic factors, subject to rapid variation as a result of the current global warming, are known to strongly influence the fate of alpine plants. Our study shows, for the first time for a large number of species, that the same environmental variables are drivers of plant adaptation at the scale of a whole biome, here the European Alps.

Automated scoring of AFLPs using RawGeno v 2.0, a free R CRAN library.

Amplified Fragment Length Polymorphisms (AFLPs) are a cheap and efficient protocol for generating large sets of genetic markers. This technique has become increasingly used during the last decade in various fields of biology, including population genomics, phylogeography, and genome mapping. Here, we present RawGeno, an R library dedicated to the automated scoring of AFLPs (i.e., the coding of electropherogram signals into ready-to-use datasets). Our program includes a complete suite of tools for binning, editing, visualizing, and exporting results obtained from AFLP experiments. RawGeno can either be used with command lines and program analysis routines or through a user-friendly graphical user interface. We describe the whole RawGeno pipeline along with recommendations for (a) setting the analysis of electropherograms in combination with PeakScanner, a program freely distributed by Applied Biosystems; (b) performing quality checks; (c) defining bins and proceeding to scoring; (d) filtering nonoptimal bins; and (e) exporting results in different formats.

Integrating species distribution models (SDMs) and phylogeography for two species of Alpine Primula.

The major intention of the present study was to investigate whether an approach combining the use of niche-based palaeodistribution modeling and phylo-geography would support or modify hypotheses about the Quaternary distributional history derived from phylogeographic methods alone. Our study system comprised two closely related species of Alpine Primula. We used species distribution models based on the extant distribution of the species and last glacial maximum (LGM) climate models to predict the distribution of the two species during the LGM. Phylogeographic data were generated using amplified fragment length polymorphisms (AFLPs). In Primula hirsuta, models of past distribution and phylogeographic data are partly congruent and support the hypothesis of widespread nunatak survival in the Central Alps. Species distribution models (SDMs) allowed us to differentiate between alpine regions that harbor potential nunatak areas and regions that have been colonized from other areas. SDMs revealed that diversity is a good indicator for nunataks, while rarity is a good indicator for peripheral relict populations that were not source for the recolonization of the inner Alps. In P. daonensis, palaeo-distribution models and phylogeographic data are incongruent. Besides the uncertainty inherent to this type of modeling approach (e.g., relatively coarse 1-km grain size), disagreement of models and data may partly be caused by shifts of ecological niche in both species. Nevertheless, we demonstrate that the combination of palaeo-distribution modeling with phylogeographical approaches provides a more differentiated picture of the distributional history of species and partly supports (P. hirsuta) and partly modifies (P. daonensis and P. hirsuta) hypotheses of Quaternary distributional history. Some of the refugial area indicated by palaeodistribution models could not have been identified with phylogeographic data.

Cadmium hyperaccumulation and genetic differentiation of Thlaspi caerulescens populations

Although the knowledge on heavy metal hyperaccumulation mechanisms is increasing, the genetic basis of cadmium (Cd) hyperaccurnulation remains to be elucidated. Thlaspi caerulescens is an attractive model since Cd accumulation polymorphism observed in this species suggests genetic differences between populations with low versus high Cd hyperaccumulation capacities. In our study, a methodology is proposed to analyse at a regional scale the genetic differentiation of T. caerulescens natural populations in relation to Cd hyperaccumulation capacity while controlling for different environmental, soil, plant parameters and geographic origins of populations. Twenty-two populations were characterised with AFLP markers and cpDNA polymorphism. Over all loci, a partial Mantel test showed no significant genetic structure with regard to the Cd hyperaccumulation capacity. Nevertheless, when comparing the marker variation to a neutral model, seven AFLP fragments (9% of markers) were identified as presenting particularly high genetic differentiation between populations with low and high Cd hyperaccurnulation capacity. Using simulations, the number of outlier loci was showed to be significantly higher than expected at random. These loci presented a genetic structure linked to Cd hyperaccumulation capacity independently of the geography, environment, soil parameters and Zn, Pb, Fe and Cu concentrations in plants. Using a canonical correspondence analysis, we identified three of them as particularly related to the Cd hyperaccumutation capacity. This study demonstrates that populations with low and high hyperaccurnulation capacities can be significantly distinguished based on molecular data. Further investigations with candidate genes and mapped markers may allow identification and characterization of genomic regions linked to factors involved in Cd hyperaccumulation.

Morphological, ecological and genetic aspects associated with endemism in the Fly Orchid group.

The European genus Ophrys (Orchidaceae) is famous for its insect-like floral morphology, an adaptation for a pseudocopulatory pollination strategy involving Hymenoptera males. A large number of endemic Ophrys species have recently been described, especially within the Mediterranean Basin, which is one of the major species diversity hotspots. Subtle morphological variation and specific pollinator dependence are the two main perceptible criteria for describing numerous endemic taxa. However, the degree to which endemics differ genetically remains a challenging question. Additionally, knowledge regarding the factors underlying the emergence of such endemic entities is limited. To achieve new insights regarding speciation processes in Ophrys, we have investigated species boundaries in the Fly Orchid group (Ophrys insectifera sensu lato) by examining morphological, ecological and genetic evidence. Classically, authors have recognized one widespread taxon (O. insectifera) and two endemics (O. aymoninii from France and O. subinsectifera from Spain). Our research has identified clear morphological and ecological factors segregating among these taxa; however, genetic differences were more ambiguous. Insights from cpDNA sequencing and amplified fragment length polymorphisms genotyping indicated a recent diversification in the three extant Fly Orchid species, which may have been further obscured by active migration and admixture across the European continent. Our genetic results still indicate weak but noticeable phylogeographic clustering that partially correlates with the described species. Particularly, we report several isolated haplotypes and genetic clusters in central and southeastern Europe. With regard to the morphological, ecological and genetic aspects, we discuss the endemism status within the Fly Orchid group from evolutionary, taxonomical and conservation perspectives.

Characterization of microsatellite loci and reliable genotyping in a polyploid plant, Mercurialis perennis (Euphorbiaceae).

For many applications in population genetics, codominant simple sequence repeats (SSRs) may have substantial advantages over dominant anonymous markers such as amplified fragment length polymorphisms (AFLPs). In high polyploids, however, allele dosage of SSRs cannot easily be determined and alleles are not easily attributable to potentially diploidized loci. Here, we argue that SSRs may nonetheless be better than AFLPs for polyploid taxa if they are analyzed as effectively dominant markers because they are more reliable and more precise. We describe the transfer of SSRs developed for diploid Mercurialis huetii to the clonal dioecious M. perennis. Primers were tested on a set of 54 male and female plants from natural decaploid populations. Eight of 65 tested loci produced polymorphic fragments. Binary profiles from 4 different scoring routines were used to define multilocus lineages (MLLs). Allowing for fragment differences within 1 MLL, all analyses revealed the same 14 MLLs without conflicting with merigenet, sex, or plot assignment. For semiautomatic scoring, a combination of as few as 2 of the 4 most polymorphic loci resulted in unambiguous discrimination of clones. Our study demonstrates that microsatellite fingerprinting of polyploid plants is a cost efficient and reliable alternative to AFLPs, not least because fewer loci are required than for diploids.

Fine-scale genetic structure and marginal processes in an expanding population of Biscutella laevigata L. (Brassicaceae).

Evolutionary processes acting at the expanding margins of a species' range are still poorly understood. Genetic drift is considered prevalent in marginal populations, and the maintenance of genetic diversity during recolonization might seem puzzling. To investigate such processes, a fine-scale investigation of 219 individuals was performed within a population of Biscutella laevigata (Brassicaceae), located at the leading edge of its range. The survey used amplified fragment length polymorphisms (AFLPs). As commonly reported across the whole species distribution range, individual density and genetic diversity decreased along the local axis of recolonization of this expanding population, highlighting the enduring effect of the historical colonization on present-day diversity. The self-incompatibility system of the plant may have prevented local inbreeding in newly found patches and sustained genetic diversity by ensuring gene flow from established populations. Within the more continuously populated region, spatial analysis of genetic structure revealed restricted gene flow among individuals. The distribution of genotypes formed a mosaic of relatively homogenous patches within the continuous population. This pattern could be explained by a history of expansion by long-distance dispersal followed by fine-scale diffusion (that is, a stratified dispersal combination). The secondary contact among expanding patches apparently led to admixture among differentiated genotypes where they met (that is, a reshuffling effect). This type of dynamics could explain the maintenance of genetic diversity during recolonization.

Combining palaeodistribution modelling and phylogeographical approaches for identifying glacial refugia in Alpine Primula

Aim We investigated the late Quaternary history of two closely related and partly sympatric species of Primula from the south-western European Alps, P. latifolia Lapeyr. and P. marginata Curtis, by combining phylogeographical and palaeodistribution modelling approaches. In particular, we were interested in whether the two approaches were congruent and identified the same glacial refugia. Location South-western European Alps. Methods For the phylogeographical analysis we included 353 individuals from 28 populations of P. marginata and 172 individuals from 15 populations of P. latifolia and used amplified fragment length polymorphisms (AFLPs). For palaeodistribution modelling, species distribution models (SDMs) were based on extant species occurrences and then projected to climate models (CCSM, MIROC) of the Last Glacial Maximum (LGM), approximately 21 ka. Results The locations of the modelled LGM refugia were confirmed by various indices of genetic variation. The refugia of the two species were largely geographically isolated, overlapping only 6% to 11% of the species' total LGM distribution. This overlap decreased when the position of the glacial ice sheet and the differential elevational and edaphic distributions of the two species were considered. Main conclusions The combination of phylogeography and palaeodistribution modelling proved useful in locating putative glacial refugia of two alpine species of Primula. The phylogeographical data allowed us to identify those parts of the modelled LGM refugial area that were likely source areas for recolonization. The use of SDMs predicted LGM refugial areas substantially larger and geographically more divergent than could have been predicted by phylogeographical data alone

Combining molecular evolution and environmental genomics to unravel adaptive processes of MHC class IIB diversity in European minnows (Phoxinus phoxinus)

Host-pathogen interactions are a major evolutionary force promoting local adaptation. Genes of the major histocompatibility complex (MHC) represent unique candidates to investigate evolutionary processes driving local adaptation to parasite communities. The present study aimed at identifying the relative roles of neutral and adaptive processes driving the evolution of MHC class IIB (MHCIIB) genes in natural populations of European minnows (Phoxinus phoxinus). To this end, we isolated and genotyped exon 2 of two MHCIIB gene duplicates (DAB1 and DAB3) and 1665 amplified fragment length polymorphism (AFLP) markers in nine populations, and characterized local bacterial communities by 16S rDNA barcoding using 454 amplicon sequencing. Both MHCIIB loci exhibited signs of historical balancing selection. Whereas genetic differentiation exceeded that of neutral markers at both loci, the populations' genetic diversities were positively correlated with local pathogen diversities only at DAB3. Overall, our results suggest pathogen-mediated local adaptation in European minnows at both MHCIIB loci. While at DAB1 selection appears to favor different alleles among populations, this is only partially the case in DAB3, which appears to be locally adapted to pathogen communities in terms of genetic diversity. These results provide new insights into the importance of host-pathogen interactions in driving local adaptation in the European minnow, and highlight that the importance of adaptive processes driving MHCIIB gene evolution may differ among duplicates within species, presumably as a consequence of alternative selective regimes or different genomic context.

Ecological temporal stability of Staphylococcus aureus nasal carriage.

We described the colonization dynamics of Staphylococcus aureus in a group of 266 healthy carriers over a period of approximately 1 year. We used precise genotyping methods, i.e., amplified fragment length polymorphism (AFLP), spa typing, and double-locus sequence typing (DLST), to detect changes in strain identity. Strain change took place rather rarely: out of 89 carriers who had initially been colonized, only 7 acquired a strain different from the original one. Approximately one-third of the carriers eliminated the colonization, and a similar number became newly colonized. Some of these events probably represent detection failure rather than genuine colonization loss or acquisition. Lower bacterial counts were associated with increased probability of eliminating the colonization. We have confirmed a high mutation rate in the spa locus: 6 out of 53 strains underwent mutation in the spa locus. There was no overall change in S. aureus genotype composition.

Link between genotype and antimicrobial resistance in bovine mastitis-related Staphylococcus aureus strains, determined by comparing Swiss and French isolates from the Rhône Valley.

Staphylococcus aureus is a major bovine mastitis pathogen. Although the reported antimicrobial resistance was generally low, the emergence of new genetic clusters in bovine mastitis requires examination of the link between antimicrobial resistance and genotypes. Here, amplified fragment length polymorphism (AFLP) profiles and standard antimicrobial resistance profiles were determined in order to characterize a total of 343 S. aureus cow mastitis isolates from two geographically close regions of Switzerland and France. AFLP profiles revealed similar population compositions in the two regions, with 4 major clusters (C8, C20, C97, and C151), but the proportions of isolates in each cluster significantly diverged between the two countries (P = 9.2 × 10⁻⁹). Antimicrobial resistance was overall low (< 5% resistance to all therapeutically relevant molecules), with the exception of penicillin resistance, which was detected in 26% of the isolates. Penicillin resistance proportions differed between clusters, with only 1 to 2% of resistance associated with C20 and C151 and up to 70% associated with bovine C97. The prevalence of C20 and C8 was unexpectedly high and requires further investigation into the mechanism of adaptation to the bovine host. The strong association of penicillin resistance with few clusters highlights the fact that the knowledge of local epidemiology is essential for rational choices of antimicrobial treatment in the absence of susceptibility testing. Taken together, these observations argue in favor of more routine scrutiny of antimicrobial resistance and antibiotic-resistant clones in cattle and the farm environment.

Antimicrobial resistance of Staphylococcus aureus strains acquired by pig farmers from pigs

Carriage of animal-associated methicillin-resistant Staphylococcus aureus (MRSA) clonal complex 398 (CC398) is common among pig farmers. This study was conducted (i) to investigate whether pig farmers are colonized with pig-specific S. aureus genotypes other than CC398 and (ii) to survey antimicrobial resistance of S. aureus isolates from pigs and pig farmers. Forty-eight S. aureus isolates from pig farmers and veterinarians and 130 isolates from pigs collected in Western Switzerland were genotyped by spa typing and amplified fragment length polymorphism (AFLP). Antimicrobial resistance profiles were determined for representative sample of the isolates. The data obtained earlier on healthy S. aureus carriers without exposure to agriculture were used for comparison. The genotype composition of S. aureus isolates from pig farmers and veterinarians was similar to isolates from pigs with predominant AFLP clusters CC398, CC9, and CC49. The resistance to tetracycline and macrolides (clarithromycin) was common among the isolates from farmers and veterinarians (52 and 21%, respectively) and similar to resistance levels in isolates from pigs (39 and 23%, respectively). This was in contrast to isolates from persons without contact with agriculture, where no (0/128) isolates were resistant to tetracycline and 3% of the isolates were resistant to clarithromycin. MRSA CC398 was isolated from pigs (n = 11) and pig farmers (n = 5). These data imply that zoonotic transmission of multidrug-resistant S. aureus from pigs to farmers is frequent, and well-known MRSA transmission merely represents the tip of the iceberg for this phenomenon. We speculate that the relatively low frequency of MRSA isolation is related to lower antimicrobial use in Switzerland compared to, for example, the Netherlands.

Origin and expansion of the allotetraploid Aegilops geniculata, a wild relative of wheat.

*This study reconstructs the phylogeography of Aegilops geniculata, an allotetraploid relative of wheat, to discuss the impact of past climate changes and recent human activities (e.g. the early expansion of agriculture) on the genetic diversity of ruderal plant species. *We combined chloroplast DNA (cpDNA) sequencing, analysed using statistical parsimony network, with nonhierarchical K-means clustering of amplified fragment length polymorphism (AFLP) genotyping, to unravel patterns of genetic structure across the native range of Ae. geniculata. The AFLP dataset was further explored by measurement of the regional genetic diversity and the detection of isolation by distance patterns. *Both cpDNA and AFLP suggest an eastern Mediterranean origin of Ae. geniculata. Two lineages have spread independently over northern and southern Mediterranean areas. Northern populations show low genetic diversity but strong phylogeographical structure among the main peninsulas, indicating a major influence of glacial cycles. By contrast, low genetic structuring and a high genetic diversity are detected in southern Mediterranean populations. Finally, we highlight human-mediated dispersal resulting in substantial introgression between resident and migrant populations. *We have shown that the evolutionary trajectories of ruderal plants can be similar to those of wild species, but are interfered by human activities, promoting range expansions through increased long-distance dispersal and the creation of suitable habitats.

Staphylococcus aureus host range and human-bovine host shift.

Staphylococcus aureus is a major agent of bovine mastitis. The concomitant emergence of pig-associated methicillin-resistant S. aureus (MRSA) in human carriage and infection requires a reexamination of the host range and specificity of human- and cow-associated S. aureus strains, something which has not been systematically studied previously. The genetic relatedness of 500 S. aureus isolates from bovine mastitis cases, 57 isolates from nasal carriage of farmers, and 133 isolates from nonfarmers was determined by amplified fragment length polymorphism (AFLP) analysis and spa typing. Multilocus sequence typing (MLST) was conducted on a subset of isolates to match AFLP clusters with MLST clonal complexes (CCs). This data set allowed us to study host range and host specificity and to estimate the extent of bovine-to-human transmission. The genotype compositions of S. aureus isolates from farmers and nonfarmers were very similar, while the mastitis isolates were quite distinct. Overall, transmission was low, but specific genotypes did show increased cow-to-human transmission. Unexpectedly, more than one-third of mastitis isolates belonged to CC8, a lineage which has not been considered to be bovine mastitis associated, but it is well known from human carriage and infection (i.e., USA300). Despite the fact that we did detect some transmission of other genotypes from cows to farmers, no transmission of CC8 isolates to farmers was detected, except for one tentative case. This was despite the close genetic relatedness of mastitis CC8 strains to nonfarmer carriage strains. These results suggest that the emergence of the new bovine-adapted genotype was due to a recent host shift from humans to cows concurrent with a loss of the ability to colonize humans. More broadly, our results indicate that host specificity is a lineage-specific trait that can rapidly evolve.