Conserving marine biodiversity: insights from life-history trait candidate genes in Atlantic cod (Gadus morhua)

Recent technological developments have facilitated an increased focus on identifying genomic regions underlying adaptive trait variation in natural populations, and it has been advocated that this information should be important for designating population units for conservation. In marine fishes, phenotypic studies have suggested adaptation through divergence of life-history traits among natural populations, but the distribution of adaptive genetic variation in these species is still relatively poorly known. In this study, we extract information about the geographical distribution of genetic variation for 33 single nucleotide polymorphisms (SNPs) associated with life-history trait candidate genes, and compare this to variation in 70 putatively neutral SNPs in Atlantic cod (Gadus morhua). We analyse samples covering the major population complexes in the eastern Atlantic and find strong evidence for non-neutral levels and patterns of population structuring for several of the candidate gene-associated markers, including two SNPs in the growth hormone 1 gene. Thus, this study aligns with findings from phenotypic studies, providing molecular data strongly suggesting that these or closely linked genes are under selection in natural populations of Atlantic cod. Furthermore, we find that patterns of variation in outlier markers do not align with those observed at selectively neutral markers, and that outlier markers identify conservation units on finer geographical scales than those revealed when analysing only neutral markers. Accordingly, results also suggest that information about adaptive genetic variation will be useful for targeted conservation and management in this and other marine species.     

Pedigrees, MHC and microsatellites: an integrated approach for genetic management of captive African wild dogs (Lycaon pictus)

Captive breeding programmes aim to provide an insurance against extinction in the wild and a source for re-introductions making it essential to minimise genetic threats, and maximise representation of wild adaptive genetic diversity. As such, genetic assessments of captive breeding programmes are increasingly common. However, these rarely include comprehensive comparisons with wild populations and typically neutral, rather than adaptive, genetic diversity is assayed. Moreover, genetic data are rarely integrated with studbook information, which enables the most robust assessments. Here we use the European captive African wild dog (Lycaon pictus) population to demonstrate the utility of this combined approach. Specifically, we combined studbook pedigree information with genetic assessments of captive and wild samples at both neutral markers and a locus thought to be important for adaptation (a gene at the Major Histocompatibility Complex, MHC). With these data we were able to evaluate founder origin and representation, as well as the distribution and origin of genetic variation within the captive population. We found discrepancies between diversity metrics derived from neutral and adaptive markers and pedigree versus genetic derived inbreeding estimates. Overall, however, we found a large proportion of genetic diversity from wild populations to be conserved in the captive population, much of which can be attributed to recent imports from outside of the European breeding programme. Nonetheless, we also found a high incidence of inbreeding and very skewed founder contributions. Based on these results, we proposed and implemented a genetic management plan to prevent further losses of diversity and reduce inbreeding.     

On the relative roles of selection and genetic drift in shaping MHC variation

Genes of the major histocompatibility complex (MHC) have provided some of the clearest examples of how natural selection generates discordances between adaptive and neutral variation in natural populations. The type and intensity of selection as well as the strength of genetic drift are believed to be important in shaping the resulting pattern of MHC diversity. However, evaluating the relative contribution of multiple microevolutionary forces is challenging, and empirical studies have reported contrasting results. For instance, balancing selection has been invoked to explain high levels of MHC diversity and low population differentiation in comparison with other nuclear markers. Other studies have shown that genetic drift can sometimes overcome selection and then patterns of genetic variation at adaptive loci cannot be discerned from those occurring at neutral markers. Both empirical and simulated data also indicate that loss of genetic diversity at adaptive loci can occur faster than at neutral loci when selection and population bottlenecks act simultaneously. Diversifying selection, on the other hand, explains accelerated MHC divergence as the result of spatial variation in pathogen‐mediated selective regimes. Because of all these possible scenarios and outcomes, collecting information from as many study systems as possible, is crucial to enhance our understanding about the evolutionary forces driving MHC polymorphism. In this issue, Miller and co‐workers present an illuminating contribution by combining neutral markers (microsatellites) and adaptive MHC class I loci during the investigation of genetic differentiation across island populations of tuatara Sphenodon punctatus. Their study of geographical variation reveals a major role of genetic drift in shaping MHC variation, yet they also discuss some support for diversifying selection.     

Adaptive and neutral genetic differentiation among Scottish and endangered Irish red grouse (<Emphasis Type="Italic">Lagopus lagopus scotica</Emphasis>)

Studying patterns of intra-specific genetic variation among populations allows for a better understanding of population structure and local adaptation. However, those patterns may differ according to the genetic markers applied, as neutral genetic markers reflect demographic processes and random genetic drift, whereas adaptive markers also carry the footprint of selection. In combination, neutral and adaptive genetic markers permit to assess the relative roles of drift and selection in shaping population structure. Among the best understood adaptive genetic loci are the genes of the major histocompatibility complex (MHC). We here study variation and differentiation at neutral SNP markers and MHC class II genes in red grouse (Lagopus lagopus scotica) from Ireland and Scotland. Irish red grouse populations are fragmented and drastically declining, but red grouse are abundant in Scotland. We find evidence for positive selection acting on the MHC genes and variation in MHC gene copy numbers among Irish individuals. Furthermore, there was significant population differentiation among red grouse from Ireland and Scotland at the neutral SNP markers (F_ST = 0.084) and the MHC-BLB genes (F_ST: BLB1 = 0.116, BLB2 = 0.090, BLB3 = 0.104). Differentiation at the MHC-BLB1 was significantly higher than at the neutral SNP markers, suggesting that selection plays an important role in shaping MHC variation, in addition to genetic drift. We speculate that the observed differentiation pattern might be due to local adaptation to different parasite regimes. These findings have strong conservation implications and we advise against the introduction of Scottish red grouse to supplement Irish populations.     

Adaptive and neutral markers both show continent‐wide population structure of mountain pine beetle (Dendroctonus ponderosae)

Assessments of population genetic structure and demographic history have traditionally been based on neutral markers while explicitly excluding adaptive markers. In this study, we compared the utility of putatively adaptive and neutral single‐nucleotide polymorphisms (SNPs) for inferring mountain pine beetle population structure across its geographic range. Both adaptive and neutral SNPs, and their combination, allowed range‐wide structure to be distinguished and delimited a population that has recently undergone range expansion across northern British Columbia and Alberta. Using an equal number of both adaptive and neutral SNPs revealed that adaptive SNPs resulted in a stronger correlation between sampled populations and inferred clustering. Our results suggest that adaptive SNPs should not be excluded prior to analysis from neutral SNPs as a combination of both marker sets resulted in better resolution of genetic differentiation between populations than either marker set alone. These results demonstrate the utility of adaptive loci for resolving population genetic structure in a nonmodel organism.     

Spatially explicit and multi-sourced genetic information is critical for conservation of an endangered plant species,San Diego thornmint (<Emphasis Type="Italic">Acanthomintha ilicifolia</Emphasis>)

San Diego thornmint Acanthomintha ilicifolia (Gray) Gray (Lamiaceae) is a winter herb restricted to San Diego county in the United States and Baja California Norte in Mexico. Historic records document 80 occurrences of this species, with 55 extant occurrences in San Diego County currently known. We compared three measures of genetic variation to inform ongoing conservation efforts: putatively neutral genetic structure revealed from isozyme markers, apparent cytogenetic variation confirmed using flow cytometry, and potentially adaptive morphological variation quantified in a common-garden study. Together, these data indicated that this rare endemic is genetically complex, revealing significant differentiation of neutral and potentially adaptive genetic variation among populations, and possessing at least two cytotypes, sometimes even within the same population. While additional study is required to resolve the extent of potential local adaptation in this species, conservation plans should limit the movement of germplasm among occurrences and monitor populations in order to limit potential long-term impacts to population viability. Given that these findings challenge the canonical model of genetic structure in rare plants (low genetic variation and limited genetic structure), we recommend guidelines to apply genetic information to conservation strategies.     

The gap between the concept and definitions in the Evolutionarily Significant Unit: the need to integrate neutral genetic variation and adaptive variation

The Evolutionarily Significant Unit (ESU) was conceptualized in 1986 as a conservation unit below the species level, theoretically applicable to a wide range of taxa. The concept has gained support, and various definitions or criteria, some of which are inconsistent with each other, have since been proposed. Recent critiques of the ESU have pointed out the dominance of definitions biased to the identification of long-term isolation or neutral genetic variation, which has largely ignored the adaptive components. We present here the validity of such claims and show how the ESU definitions have actually been applied in research. We surveyed scientific journals for original papers supporting ESU designations and determined who among the proponents of ESU definitions have gained wider support. Our results indicate that indeed there are inconsistencies with the original concept and with the existing definitions. Although the original concept recommended both ecological and genetic data as the basis for identification of ESUs, which reflect true evolutionary variation, recent definitions have become biased to either neutral genetic variation or adaptive variation. The definition which uses genetic data to assess neutral genetic variation (long-term isolation) has gained major support, and therefore validates the earlier claims. To bridge the gap between the original concept and the practical application, we propose the use of partial ESU and full ESU designations. The application of full ESU should be limited solely to when both information about neutral genetic variation and adaptive variation are available. In other cases, in which only a part of the variation is examined, we should use the term partial ESU (e.g., molecular-based ESU) and continue to investigate focal populations from other aspects of variations to designate full ESU.     

Genetic basis of adaptation in Arabidopsis thaliana: local adaptation at the seed dormancy QTL DOG1

Local adaptation provides an opportunity to study the genetic basis of adaptation and investigate the allelic architecture of adaptive genes. We study delay of germination 1 (DOG1), a gene controlling natural variation in seed dormancy in Arabidopsis thaliana and investigate evolution of dormancy in 41 populations distributed in four regions separated by natural barriers. Using F(ST) and Q(ST) comparisons, we compare variation at DOG1 with neutral markers and quantitative variation in seed dormancy. Patterns of genetic differentiation among populations suggest that the gene DOG1 contributes to local adaptation. Although Q(ST) for seed dormancy is not different from F(ST) for neutral markers, a correlation with variation in summer precipitation supports that seed dormancy is adaptive. We characterize dormancy variation in several F(2) -populations and show that a series of functionally distinct alleles segregate at the DOG1 locus. Theoretical models have shown that the number and effect of alleles segregatin at quantitative trait loci (QTL) have important consequences for adaptation. Our results provide support to models postulating a large number of alleles at quantitative trait loci involved in adaptation.     

Ecological and life history characteristics predict population genetic divergence of two salmonids in the same landscape

Ecological and life history characteristics such as population size, dispersal pattern, and mating system mediate the influence of genetic drift and gene flow on population subdivision. Bull trout (Salvelinus confluentus) and mountain whitefish (Prosopium williamsoni) differ markedly in spawning location, population size and mating system. Based on these differences, we predicted that bull trout would have reduced genetic variation within and greater differentiation among populations compared with mountain whitefish. To test this hypothesis, we used microsatellite markers to determine patterns of genetic divergence for each species in the Clark Fork River, Montana, USA. As predicted, bull trout had a much greater proportion of genetic variation partitioned among populations than mountain whitefish. Among all sites, FST was seven times greater for bull trout (FST = 0.304 for bull trout, 0.042 for mountain whitefish. After removing genetically differentiated high mountain lake sites for each species FST, was 10 times greater for bull trout (FST = 0.176 for bull trout; FST = 0.018 for mountain whitefish). The same characteristics that affect dispersal patterns in these species also lead to predictions about the amount and scale of adaptive divergence among populations. We provide a theoretical framework that incorporates variation in ecological and life history factors, neutral divergence, and adaptive divergence to interpret how neutral and adaptive divergence might be correlates of ecological and life history factors.     

Patterns of differentiation in a colour polymorphism and in neutral markers reveal rapid genetic changes in natural damselfly populations

The existence and mode of selection operating on heritable adaptive traits can be inferred by comparing population differentiation in neutral genetic variation between populations (often using F(ST) values) with the corresponding estimates for adaptive traits. Such comparisons indicate if selection acts in a diversifying way between populations, in which case differentiation in selected traits is expected to exceed differentiation in neutral markers [F(ST )(selected) > F(ST )(neutral)], or if negative frequency-dependent selection maintains genetic polymorphisms and pulls populations towards a common stable equilibrium [F(ST) (selected) < F(ST) (neutral)]. Here, we compared F(ST) values for putatively neutral data (obtained using amplified fragment length polymorphism) with estimates of differentiation in morph frequencies in the colour-polymorphic damselfly Ischnura elegans. We found that in the first year (2000), population differentiation in morph frequencies was significantly greater than differentiation in neutral loci, while in 2002 (only 2 years and 2 generations later), population differentiation in morph frequencies had decreased to a level significantly lower than differentiation in neutral loci. Genetic drift as an explanation for population differentiation in morph frequencies could thus be rejected in both years. These results indicate that the type and/or strength of selection on morph frequencies in this system can change substantially between years. We suggest that an approach to a common equilibrium morph frequency across all populations, driven by negative frequency-dependent selection, is the cause of these temporal changes. We conclude that inferences about selection obtained by comparing F(ST) values from neutral and adaptive genetic variation are most useful when spatial and temporal data are available from several populations and time points and when such information is combined with other ecological sources of data.     

Spatial and temporal scales of adaptive divergence in marine fishes and the implications for conservation

Knowledge of geographic and temporal scales of adaptive genetic variation is crucial to species conservation, yet understanding of these phenomena, particularly in marine systems, is scant. Until recently, the belief has been that because most marine species have highly dispersive or mobile life stages, local adaptation could occur only on broad geographic scales. This view is supported by comparatively low levels of genetic variation among populations as detected by neutral markers. Similarly, the time scale of adaptive divergence has also been assumed to be very long, requiring thousands of generations. Recent studies of a variety of species have challenged these beliefs. First, there is strong evidence of geographically structured local adaptation in physiological and morphological traits. Second, the proportion of quantitative trait variation at the among-population level ( Q _ST) is much higher than it is for neutral markers ( F _ST) and these two metrics of genetic variation are poorly correlated. Third, evidence that selection is a potent evolutionary force capable of sustaining adaptive divergence on contemporary time scales is summarized. The differing spatial and temporal scales of adaptive v. neutral genetic divergence call for a new paradigm in thinking about the relationship between phenogeography (the geography of phenotypic variation) and phylogeography (the geography of lineages) in marine species. The idea that contemporary selective processes can cause fine-scale spatial and temporal divergence underscores the need for a new emphasis on Darwinian fishery science.     

LANDSCAPE GENOMICS IN ATLANTIC SALMON (SALMO SALAR): SEARCHING FOR GENE–ENVIRONMENT INTERACTIONS DRIVING LOCAL ADAPTATION

A growing number of studies are examining the factors driving historical and contemporary evolution in wild populations. By combining surveys of genomic variation with a comprehensive assessment of environmental parameters, such studies can increase our understanding of the genomic and geographical extent of local adaptation in wild populations. We used a large‐scale landscape genomics approach to examine adaptive and neutral differentiation across 54 North American populations of Atlantic salmon representing seven previously defined genetically distinct regional groups. Over 5500 genome‐wide single nucleotide polymorphisms were genotyped in 641 individuals and 28 bulk assays of 25 pooled individuals each. Genome scans, linkage map, and 49 environmental variables were combined to conduct an innovative landscape genomic analysis. Our results provide valuable insight into the links between environmental variation and both neutral and potentially adaptive genetic divergence. In particular, we identified markers potentially under divergent selection, as well as associated selective environmental factors and biological functions with the observed adaptive divergence. Multivariate landscape genetic analysis revealed strong associations of both genetic and environmental structures. We found an enrichment of growth‐related functions among outlier markers. Climate (temperature–precipitation) and geological characteristics were significantly associated with both potentially adaptive and neutral genetic divergence and should be considered as candidate loci involved in adaptation at the regional scale in Atlantic salmon. Hence, this study significantly contributes to the improvement of tools used in modern conservation and management schemes of Atlantic salmon wild populations.     

Nuclear microsatellite loci for Arabidopsis halleri (Brassicaceae), a model species to study plant adaptation to heavy metals

? Premise of the study: Arabidopsis halleri is a model species to study the adaptation of plants to soils contaminated by zinc, cadmium, and lead. To provide a neutral genetic background with which adaptive genetic markers could be compared, we developed highly polymorphic neutral microsatellite markers. ? Methods and Results: Using a microsatellite-enriched library method, we identified 120 microsatellite loci for quantitative trait locus (QTL) mapping analysis, of which eight primer pairs were developed in a single multiplex for population genetic studies. Analyses were performed on 508 individuals from 26 populations. All loci were polymorphic with six to 23 alleles per locus. Genetic diversity varied between 0.56 and 0.76. ? Conclusions: Our results demonstrated the value of these eight microsatellite markers to investigate neutral population genetic structure in A. halleri. To increase the resolution of population genetic analyses, we suggest adding them to the 11 markers previously developed independently.     

Conserving adaptive genetic diversity in dynamic landscapes

Genetic variation supplies the raw material for adaptation, evolution and survival of populations and has therefore been a key focus of conservation biology ever since its foundation (Soulé 1985). In previous decades, the neutral component of genetic diversity (generated by mutation and shaped by drift) has been the subject of intense scientific research, fuelled by the increasing availability of molecular markers. On the other hand, the adaptive component of genetic diversity, which is shaped by the action of natural selection, has long remained elusive and difficult to assess, especially at small spatial or temporal scales (Ouborg et al. 2010). Fortunately, new technological and methodological developments now make it possible to identify loci in the genome that are influenced by selection, and thus to get a more complete view of genetic diversity. One article featured in this issue of Molecular Ecology is a good example of this recent breakthrough. Richter-Boix et al. (2011) examined a network of moor frog populations breeding in contrasting habitats in order to understand how landscape features influence patterns of genetic variation. They combined information from both neutral markers and loci putatively under selection to quantify the relative roles of selection and isolation in the evolution of fine-scale local adaptations in these populations. This study nicely illustrates how data on polymorphisms of neutral and adaptive loci can now be judiciously synthesized to help identify the best strategies for preserving adaptive variation, and more generally to enlighten conservation and population-management plans.     

Genetic isolation of fragmented populations is exacerbated by drift and selection

Reduced genetic variation at marker loci in small populations has been well documented, whereas the relationship between quantitative genetic variation and population size has attracted little empirical investigation. Here we demonstrate that both neutral and quantitative genetic variation are reduced in small populations of a fragmented plant metapopulation, and that both drift and selective change are enhanced in small populations. Measures of neutral genetic differentiation (F(ST)) and quantitative genetic differentiation (Q(ST)) in two traits were higher among small demes, and Q(ST) between small populations exceeded that expected from drift alone. This suggests that fragmented populations experience both enhanced genetic drift and divergent selection on phenotypic traits, and that drift affects variation in both neutral markers and quantitative traits. These results highlight the need to integrate natural selection into conservation genetic theory, and suggests that small populations may represent reservoirs of genetic variation adaptive within a wide range of environments.     

Genetic diversity and differentiation at MHC genes in island populations of tuatara (Sphenodon spp.)

Neutral genetic markers are commonly used to understand the effects of fragmentation and population bottlenecks on genetic variation in threatened species. Although neutral markers are useful for inferring population history, the analysis of functional genes is required to determine the significance of any observed geographical differences in variation. The genes of the major histocompatibility complex (MHC) are well‐known examples of genes of adaptive significance and are particularly relevant to conservation because of their role in pathogen resistance. In this study, we survey diversity at MHC class I loci across a range of tuatara populations. We compare the levels of MHC variation with that observed at neutral microsatellite markers to determine the relative roles of balancing selection, diversifying selection and genetic drift in shaping patterns of MHC variation in isolated populations. In general, levels of MHC variation within tuatara populations are concordant with microsatellite variation. Tuatara populations are highly differentiated at MHC genes, particularly between the northern and Cook Strait regions, and a trend towards diversifying selection across populations was observed. However, overall our results indicate that population bottlenecks and isolation have a larger influence on patterns of MHC variation in tuatara populations than selection.     

Proportioning whole-genome single-nucleotide-polymorphism diversity for the identification of geographic population structure and genetic ancestry

The identification of geographic population structure and genetic ancestry on the basis of a minimal set of genetic markers is desirable for a wide range of applications in medical and forensic sciences. However, the absence of sharp discontinuities in the neutral genetic diversity among human populations implies that, in practice, a large number of neutral markers will be required to identify the genetic ancestry of one individual. We showed that it is possible to reduce the amount of markers required for detecting continental population structure to only 10 single-nucleotide polymorphisms (SNPs), by applying a newly developed ascertainment algorithm to Affymetrix GeneChip Mapping 10K SNP array data that we obtained from samples of globally dispersed human individuals (the Y Chromosome Consortium panel). Furthermore, this set of SNPs was able to recover the genetic ancestry of individuals from all four continents represented in the original data set when applied to an independent, much larger, worldwide population data set (Centre d'Etude du Polymorphisme Humain-Human Genome Diversity Project Cell Line Panel). Finally, we provide evidence that the unusual patterns of genetic variation we observed at the respective genomic regions surrounding the five most informative SNPs is in agreement with local positive selection being the explanation for the striking SNP allele-frequency differences we found between continental groups of human populations.     

The mode of evolution of molecular markers in populations of house mice under artificial selection for locomotor behavior

A complete understanding of the mode of evolution of molecular markers is important for making inferences about different population genetic parameters, especially because a number of studies have reported patterns of allelic variation at molecular markers that are not in agreement with neutral evolutionary expectations. In the present study, house mice (Mus domesticus) from the fourteenth generation of a selection experiment for increased voluntary wheel-running activity were used to test how selection on a complex behavior affects the distribution of allelic variation by examining patterns of variation at six microsatellite and four allozyme loci. This population had a hierarchical structure that allowed for simultaneous testing of the effects of selection and genetic drift on the distribution of allelic variation by comparing observed patterns of allele frequencies and estimates of genetic divergence at multiple hierarchical levels to expectations under models of neutral evolution. The levels of genetic divergence among replicate lines and between selection groups, estimated from microsatellite data or pooled microsatellite and allozyme data, were not significantly different from expectations under neutral evolution. Furthermore, the pattern of change of allele frequencies between the base population and generation 14 was largely in agreement with expectations under neutral evolution (although the PGM locus exhibited a pattern of change within populations that was difficult to explain under neutral evolution). Overall the results generally provide support for the neutral evolution of molecular markers.     

Genomics and conservation genetics

In large part, the relevance of genetics to conservation rests on the premise that neutral marker variation in populations reflects levels of detrimental and adaptive genetic variation. Despite its prominence, this tenet has been difficult to evaluate, until now. As we discuss here, genome sequence information and new technological and bioinformatics platforms now enable comprehensive surveys of neutral variation and more direct inferences of detrimental and adaptive variation in species with sequenced genomes and in 'genome-enabled' endangered taxa. Moreover, conservation schemes could begin to consider specific pathological genetic variants. A new conservation genetic agenda would utilize data from enhanced surveys of genomic variation in endangered species to better manage functional genetic variation.     

Machine‐learning‐based detection of adaptive divergence of the stream mayfly Ephemera strigata populations

Adaptive divergence is a key mechanism shaping the genetic variation of natural populations. A central question linking ecology with evolutionary biology is how spatial environmental heterogeneity can lead to adaptive divergence among local populations within a species. In this study, using a genome scan approach to detect candidate loci under selection, we examined adaptive divergence of the stream mayfly Ephemera strigata in the Natori River Basin in northeastern Japan. We applied a new machine‐learning method (i.e., random forest) besides traditional distance‐based redundancy analysis (dbRDA) to examine relationships between environmental factors and adaptive divergence at non‐neutral loci. Spatial autocorrelation analysis based on neutral loci was employed to examine the dispersal ability of this species. We conclude the following: (a) E. strigata show altitudinal adaptive divergence among the populations in the Natori River Basin; (b) random forest showed higher resolution for detecting adaptive divergence than traditional statistical analysis; and (c) separating all markers into neutral and non‐neutral loci could provide full insight into parameters such as genetic diversity, local adaptation, and dispersal ability.