Effects of spatially varying selection on nucleotide diversity and linkage disequilibrium: insights from deer mouse globin genes

Balancing selection on an MYB transcription factor maintains the twig trichome color variation in Melastoma normale

BACKGROUND

The factors that maintain phenotypic and genetic variation within a population have received long-term attention in evolutionary biology. Here the genetic basis and evolution of the geographically widespread variation in twig trichome color (from red to white) in a shrub Melastoma normale was investigated using Pool-seq and evolutionary analyses.

RESULTS

The results show that the twig trichome coloration is under selection in different light environments and that a 6-kb region containing an R2R3 MYB transcription factor gene is the major region of divergence between the extreme red and white morphs. This gene has two highly divergent groups of alleles, one of which likely originated from introgression from another species in this genus and has risen to high frequency (> 0.6) within each of the three populations under investigation. In contrast, polymorphisms in other regions of the genome show no sign of differentiation between the two morphs, suggesting that genomic patterns of diversity have been shaped by homogenizing gene flow. Population genetics analysis reveals signals of balancing selection acting on this gene, and it is suggested that spatially varying selection is the most likely mechanism of balancing selection in this case.

CONCLUSIONS

This study demonstrate that polymorphisms on a single transcription factor gene largely confer the twig trichome color variation in M. normale, while also explaining how adaptive divergence can occur and be maintained in the face of gene flow.

Genetic variation in haemoglobin is associated with evolved changes in breathing in high-altitude deer mice

Physiological systems often have emergent properties but the effects of genetic variation on physiology are often unknown, which presents a major challenge to understanding the mechanisms of phenotypic evolution. We investigated whether genetic variants in haemoglobin (Hb) that contribute to high-altitude adaptation in deer mice (Peromyscus maniculatus) are associated with evolved changes in the control of breathing. We created F2 inter-population hybrids of highland and lowland deer mice to test for phenotypic associations of α- and β-globin variants on a mixed genetic background. Hb genotype had expected effects on Hb-O2 affinity that were associated with differences in arterial O2 saturation in hypoxia. However, high-altitude genotypes were also associated with breathing phenotypes that should contribute to enhancing O2 uptake in hypoxia. Mice with highland α-globin exhibited a more effective breathing pattern, with highland homozygotes breathing deeper but less frequently across a range of inspired O2, and this difference was comparable to the evolved changes in breathing pattern in deer mouse populations native to high altitude. The ventilatory response to hypoxia was augmented in mice that were homozygous for highland β-globin. The association of globin variants with variation in breathing phenotypes could not be recapitulated by acute manipulation of Hb-O2 affinity, because treatment with efaproxiral (a synthetic drug that acutely reduces Hb-O2 affinity) had no effect on breathing in normoxia or hypoxia. Therefore, adaptive variation in Hb may have unexpected effects on physiology in addition to the canonical function of this protein in circulatory O2 transport.

Intermolecular interactions drive protein adaptive and co-adaptive evolution at both species and population levels

Proteins are the building blocks for almost all the functions in cells. Understanding the molecular evolution of proteins and the forces that shape protein evolution is essential in understanding the basis of function and evolution. Previous studies have shown that adaptation frequently occurs at the protein surface, such as in genes involved in host–pathogen interactions. However, it remains unclear whether adaptive sites are distributed randomly or at regions associated with particular structural or functional characteristics across the genome, since many proteins lack structural or functional annotations. Here, we seek to tackle this question by combining large-scale bioinformatic prediction, structural analysis, phylogenetic inference, and population genomic analysis of Drosophila protein-coding genes. We found that protein sequence adaptation is more relevant to function-related rather than structure-related properties. Interestingly, intermolecular interactions contribute significantly to protein adaptation. We further showed that intermolecular interactions, such as physical interactions, may play a role in the coadaptation of fast-adaptive proteins. We found that strongly differentiated amino acids across geographic regions in protein-coding genes are mostly adaptive, which may contribute to the long-term adaptive evolution. This strongly indicates that a number of adaptive sites tend to be repeatedly mutated and selected throughout evolution in the past, present, and maybe future. Our results highlight the important roles of intermolecular interactions and coadaptation in the adaptive evolution of proteins both at the species and population levels.

Broad Concordance in the Spatial Distribution of Adaptive and Neutral Genetic Variation across an Elevational Gradient in Deer Mice

When species are continuously distributed across environmental gradients, the relative strength of selection and gene flow shape spatial patterns of genetic variation, potentially leading to variable levels of differentiation across loci. Determining whether adaptive genetic variation tends to be structured differently than neutral variation along environmental gradients is an open and important question in evolutionary genetics. We performed exome-wide population genomic analysis on deer mice sampled along an elevational gradient of nearly 4,000 m of vertical relief. Using a combination of selection scans, genotype-environment associations, and geographic cline analyses, we found that a large proportion of the exome has experienced a history of altitude-related selection. Elevational clines for nearly 30% of these putatively adaptive loci were shifted significantly up- or downslope of clines for loci that did not bear similar signatures of selection. Many of these selection targets can be plausibly linked to known phenotypic differences between highland and lowland deer mice, although the vast majority of these candidates have not been reported in other studies of highland taxa. Together, these results suggest new hypotheses about the genetic basis of physiological adaptation to high altitude, and the spatial distribution of adaptive genetic variation along environmental gradients.

Limited Evidence for Parallel Evolution Among Desert-Adapted Peromyscus Deer Mice

Warming climate and increasing desertification urge the identification of genes involved in heat and dehydration tolerance to better inform and target biodiversity conservation efforts. Comparisons among extant desert-adapted species can highlight parallel or convergent patterns of genome evolution through the identification of shared signatures of selection. We generate a chromosome-level genome assembly for the canyon mouse (Peromyscus crinitus) and test for a signature of parallel evolution by comparing signatures of selective sweeps across population-level genomic resequencing data from another congeneric desert specialist (Peromyscus eremicus) and a widely distributed habitat generalist (Peromyscus maniculatus), that may be locally adapted to arid conditions. We identify few shared candidate loci involved in desert adaptation and do not find support for a shared pattern of parallel evolution. Instead, we hypothesize divergent molecular mechanisms of desert adaptation among deer mice, potentially tied to species-specific historical demography, which may limit or enhance adaptation. We identify a number of candidate loci experiencing selective sweeps in the P. crinitus genome that are implicated in osmoregulation (Trypsin, Prostasin) and metabolic tuning (Kallikrein, eIF2-alpha kinase GCN2, APPL1/2), which may be important for accommodating hot and dry environmental conditions.

High-Altitude Adaptation: Mechanistic Insights from Integrated Genomics and Physiology

Population genomic analyses of high-altitude humans and other vertebrates have identified numerous candidate genes for hypoxia adaptation, and the physiological pathways implicated by such analyses suggest testable hypotheses about underlying mechanisms. Studies of highland natives that integrate genomic data with experimental measures of physiological performance capacities and subordinate traits are revealing associations between genotypes (e.g., hypoxia-inducible factor gene variants) and hypoxia-responsive phenotypes. The subsequent search for causal mechanisms is complicated by the fact that observed genotypic associations with hypoxia-induced phenotypes may reflect second-order consequences of selection-mediated changes in other (unmeasured) traits that are coupled with the focal trait via feedback regulation. Manipulative experiments to decipher circuits of feedback control and patterns of phenotypic integration can help identify causal relationships that underlie observed genotype–phenotype associations. Such experiments are critical for correct inferences about phenotypic targets of selection and mechanisms of adaptation.

Evolution of Hemoglobin Genes in a Subterranean Rodent Species (Lasiopodomys mandarinus)

The Mandarin vole (Lasiopodomys mandarinus), a typical subterranean rodent, has undergone hematological adaptations to tolerate the hypoxic/hypercapnic underground environment. Hemoglobin (Hb) genes encode respiratory proteins functioning principally in oxygen binding and transport to various tissues and organs. To investigate the evolution of α- and β-hemoglobin (Hb) in subterranean rodent species, we sequenced Hb genes of the Mandarin vole and the related aboveground Brandt's vole (L. brandtii). Sequencing showed that in both voles, α-globin was encoded by a cluster of five functional genes in the following linkage order: HBZ, HBA-T1, HBQ-T1, HBA-T2, and HBQ-T2; among these, HBQ-T2 is a pseudogene in both voles. The β-globin gene cluster in both voles also included five functional genes in the following linkage order: HBE, HBE/HBG, HBG, HBB-T1, and HBB-T2. Phylogenetic analysis revealed that the Mandarin vole underwent convergent evolution with its related aboveground species (Brandt's vole) but not with other subterranean rodent species. Selection pressure analyses revealed that α- and β-globin genes are under strong purifying selection (ω < 1), and branch-site analyses identified positive selection sites on HBAT-T1 and HBB-T1 in different subterranean rodent species. This suggests that the adaptive evolution of these genes enhanced the ability of Hb to store and transport oxygen in subterranean rodent species. Our findings highlight the critical roles of Hb genes in the evolution of hypoxia tolerance in subterranean rodent species.

Evolutionary Diversity in the Intracellular Microsporidian Parasite Nosema sp. Infecting Wild Silkworm Revealed by IGS Nucleotide Sequence Diversity

Intracellular microsporidian Nosema mylitta infects Indian wild silkworm Antheraea mylitta causing pebrine disease. Genetic structure and phylogeny of N. mylitta are analysed using nucleotide variability in 5S ribosomal DNA and intergenic spacer (IGS) sequence from 20 isolates collected from Southern, Northern and Central regions of Jharkhand State. Nucleotide diversity (π) and genetic differentiation Gst were highest in the Central isolates whereas lowest in the North. Among the isolates, absence of nucleotides, transitions and transversions were observed. Haplotyping showed nucleotide variability at 83 positions in IGS and 13 positions in 5S rDNA. Haplotype-based genetic differentiation was 0.96 to 0.97 whereas nucleotide sequence-based genetic differentiation was higher (Ks = 22.29) between Southern and Central isolates. Bottleneck analysis showed negative value for Tajima's D and other summary statistics revealing induction of loss of rare alleles and population explosion. From IGS, 17 ancestral sequences were inferred by Network algorithm. Core of nine closely related nodes having ancient nucleotides and peripheral nodes with highly divergent nucleotides were derived. Most diverged peripheral haplotype was Bero (H11) from the Central region whereas Deoghar (H3) of the Northern region diverged early. Phylogeny of N. mylitta grouped Southern and Northern isolates together revealed weak phylogenetic signal for these locations. Phylogeny of N. mylitta with Nosema sp. infecting other lepidopterans clustered N. mylitta isolates with N. antheraea and N. philosamiae of China indicating genetic similarity whereas other species were dissimilar showing diversity irrespective of country of origin.

Ontogenesis of evolved changes in respiratory physiology in deer mice native to high altitude

High-altitude environments are cold and hypoxic, and many high-altitude natives have evolved changes in respiratory physiology that improve O2 uptake in hypoxia as adults. Altricial mammals undergo a dramatic metabolic transition from ectothermy to endothermy in early post-natal life, which may influence the ontogenetic development of respiratory traits at high altitude. We examined the developmental changes in respiratory and haematological traits in deer mice (Peromyscus maniculatus) native to high altitude, comparing the respiratory responses to progressive hypoxia between highland and lowland deer mice. Among adults, highlanders exhibited higher total ventilation and a more effective breathing pattern (relatively deeper tidal volumes), for mice that were caught and tested at their native altitudes and those lab-raised in normoxia. Lab-raised progeny of each population were also tested at post-natal day (P)7, 14, 21 and 30. Highlanders developed an enhanced hypoxic ventilatory response by P21, concurrent with the full maturation of the carotid bodies, and their more effective breathing pattern arose by P14; these ages correspond to critical benchmarks in the full development of homeothermy in highlanders. However, highlanders exhibited developmental delays in ventilatory sensitivity to hypoxia, hyperplasia of type I cells in the carotid body and increases in blood haemoglobin content compared with lowland mice. Nevertheless, highlanders maintained consistently higher arterial O2 saturation in hypoxia across development, in association with increases in blood-O2 affinity that were apparent from birth. We conclude that evolved changes in respiratory physiology in high-altitude deer mice become expressed in association with the post-natal development of endothermy.

The fitness challenge of studying molecular adaptation

Advances in bioinformatics and high-throughput genetic analysis increasingly allow us to predict the genetic basis of adaptive traits. These predictions can be tested and confirmed, but the molecular-level changes - i.e. the molecular adaptation - that link genetic differences to organism fitness remain generally unknown. In recent years, a series of studies have started to unpick the mechanisms of adaptation at the molecular level. In particular, this work has examined how changes in protein function, activity, and regulation cause improved organismal fitness. Key to addressing molecular adaptations is identifying systems and designing experiments that integrate changes in the genome, protein chemistry (molecular phenotype), and fitness. Knowledge of the molecular changes underpinning adaptations allow new insight into the constraints on, and repeatability of adaptations, and of the basis of non-additive interactions between adaptive mutations. Here we critically discuss a series of studies that examine the molecular-level adaptations that connect genetic changes and fitness.

Evolution of physiological performance capacities and environmental adaptation: insights from high-elevation deer mice (Peromyscus maniculatus)

Analysis of variation in whole-animal performance can shed light on causal connections between specific traits, integrated physiological capacities, and Darwinian fitness. Here, we review and synthesize information on naturally occurring variation in physiological performance capacities and how it relates to environmental adaptation in deer mice (Peromyscus maniculatus). We discuss how evolved changes in aerobic exercise capacity and thermogenic capacity have contributed to adaptation to high elevations. Comparative work on deer mice at high and low elevations has revealed evolved differences in aerobic performance capacities in hypoxia. Highland deer mice have consistently higher aerobic performance capacities under hypoxia relative to lowland natives, consistent with the idea that it is beneficial to have a higher maximal metabolic rate (as measured by the maximal rate of O₂ consumption, VO_2max) in an environment characterized by lower air temperatures and lower O₂ availability. Observed differences in aerobic performance capacities between highland and lowland deer mice stem from changes in numerous subordinate traits that alter the flux capacity of the O₂-transport system, the oxidative capacity of tissue mitochondria, and the relationship between O₂ consumption and ATP synthesis. Many such changes in physiological phenotype are associated with hypoxia-induced changes in gene expression. Research on natural variation in whole-animal performance forms a nexus between physiological ecology and evolutionary biology that requires insight into the natural history of the study species.

Water availability drives signatures of local adaptation in whitebark pine (Pinus albicaulis Engelm.) across fine spatial scales of the Lake Tahoe Basin, USA

Patterns of local adaptation at fine spatial scales are central to understanding how evolution proceeds, and are essential to the effective management of economically and ecologically important forest tree species. Here, we employ single and multilocus analyses of genetic data (n = 116 231 SNPs) to describe signatures of fine-scale adaptation within eight whitebark pine (Pinus albicaulis Engelm.) populations across the local extent of the environmentally heterogeneous Lake Tahoe Basin, USA. We show that despite highly shared genetic variation (F_ST  = 0.0069), there is strong evidence for adaptation to the rain shadow experienced across the eastern Sierra Nevada. Specifically, we build upon evidence from a common garden study and find that allele frequencies of loci associated with four phenotypes (mean = 236 SNPs), 18 environmental variables (mean = 99 SNPs), and those detected through genetic differentiation (n = 110 SNPs) exhibit significantly higher signals of selection (covariance of allele frequencies) than could be expected to arise, given the data. We also provide evidence that this covariance tracks environmental measures related to soil water availability through subtle allele frequency shifts across populations. Our results replicate empirical support for theoretical expectations of local adaptation for populations exhibiting strong gene flow and high selective pressures and suggest that ongoing adaptation of many P. albicaulis populations within the Lake Tahoe Basin will not be constrained by the lack of genetic variation. Even so, some populations exhibit low levels of heritability for the traits presumed to be related to fitness. These instances could be used to prioritize management to maintain adaptive potential. Overall, we suggest that established practices regarding whitebark pine conservation be maintained, with the additional context of fine-scale adaptation.

Population structure and local selection yield high genomic variation in Mimulus guttatus

Across western North America, Mimulus guttatus exists as many local populations adapted to site-specific environmental challenges. Gene flow between locally adapted populations will affect genetic diversity both within demes and across the larger metapopulation. Here, we analyse 34 whole-genome sequences from the intensively studied Iron Mountain population (IM) in conjunction with sequences from 22 Mimulus individuals sampled from across western North America. Three striking features of these data address hypotheses about migration and selection in a locally adapted population. First, we find very high levels of intrapopulation polymorphism (synonymous π = 0.033). Variation outside of genes is likely even higher but difficult to estimate because excessive divergence reduces the efficiency of read mapping. Second, IM exhibits a significantly positive genomewide average for Tajima's D. This indicates allele frequencies are typically more intermediate than expected from neutrality, opposite the pattern observed in many other species. Third, IM exhibits a distinctive haplotype structure with a genomewide excess of positive associations between rarer alleles at linked loci. This suggests an important effect of gene flow from other Mimulus populations, although a residual effect of population founding might also contribute. The combination of multiple analyses, including a novel tree-based analytic method, illustrates how the balance of local selection, limited dispersal and metapopulation dynamics manifests across the genome. The overall genomic pattern of sequence diversity suggests successful gene flow of divergent immigrant genotypes into IM. However, many loci show patterns indicative of local adaptation, particularly at SNPs associated with chromosomal inversions.

Single-Locus versus Multilocus Patterns of Local Adaptation to Climate in Eastern White Pine (Pinus strobus, Pinaceae)

Natural plant populations are often adapted to their local climate and environmental conditions, and populations of forest trees offer some of the best examples of this pattern. However, little empirical work has focused on the relative contribution of single-locus versus multilocus effects to the genetic architecture of local adaptation in plants/forest trees. Here, we employ eastern white pine (Pinus strobus) to test the hypothesis that it is the inter-genic effects that primarily drive climate-induced local adaptation. The genetic structure of 29 range-wide natural populations of eastern white pine was determined in relation to local climatic factors using both a reference set of SSR markers, and SNPs located in candidate genes putatively involved in adaptive response to climate. Comparisons were made between marker sets using standard single-locus outlier analysis, single-locus and multilocus environment association analyses and a novel implementation of Population Graphs. Magnitudes of population structure were similar between the two marker sets. Outlier loci consistent with diversifying selection were rare for both SNPs and SSRs. However, genetic distances based on the multilocus among population covariances (cGD) were significantly more correlated to climate, even after correcting for spatial effects, for SNPs as compared to SSRs. Coalescent simulations confirmed that the differences in mutation rates between SSRs and SNPs did not affect the topologies of the Population Graphs, and hence values of cGD and their correlations with associated climate variables. We conclude that the multilocus covariances among populations primarily reflect adaptation to local climate and environment in eastern white pine. This result highlights the complexity of the genetic architecture of adaptive traits, as well as the need to consider multilocus effects in studies of local adaptation.

A haplotype method detects diverse scenarios of local adaptation from genomic sequence variation

Identifying genomic targets of population-specific positive selection is a major goal in several areas of basic and applied biology. However, it is unclear how often such selection should act on new mutations versus standing genetic variation or recurrent mutation, and furthermore, favoured alleles may either become fixed or remain variable in the population. Very few population genetic statistics are sensitive to all of these modes of selection. Here, we introduce and evaluate the Comparative Haplotype Identity statistic (χMD ), which assesses whether pairwise haplotype sharing at a locus in one population is unusually large compared with another population, relative to genomewide trends. Using simulations that emulate human and Drosophila genetic variation, we find that χMD is sensitive to a wide range of selection scenarios, and for some very challenging cases (e.g. partial soft sweeps), it outperforms other two-population statistics. We also find that, as with FST , our haplotype approach has the ability to detect surprisingly ancient selective sweeps. Particularly for the scenarios resembling human variation, we find that χMD outperforms other frequency- and haplotype-based statistics for soft and/or partial selective sweeps. Applying χMD and other between-population statistics to published population genomic data from D. melanogaster, we find both shared and unique genes and functional categories identified by each statistic. The broad utility and computational simplicity of χMD will make it an especially valuable tool in the search for genes targeted by local adaptation.

Bohr effect and temperature sensitivity of hemoglobins from highland and lowland deer mice

An important means of physiological adaptation to environmental hypoxia is an increased oxygen (O2) affinity of the hemoglobin (Hb) that can help secure high O2 saturation of arterial blood. However, the trade-off associated with a high Hb-O2 affinity is that it can compromise O2 unloading in the systemic capillaries. High-altitude deer mice (Peromyscus maniculatus) have evolved an increased Hb-O2 affinity relative to lowland conspecifics, but it is not known whether they have also evolved compensatory mechanisms to facilitate O2 unloading to respiring tissues. Here we investigate the effects of pH (Bohr effect) and temperature on the O2-affinity of high- and low-altitude deer mouse Hb variants, as these properties can potentially facilitate O2 unloading to metabolizing tissues. Our experiments revealed that Bohr factors for the high- and low-altitude Hb variants are very similar in spite of the differences in O2-affinity. The Bohr factors of deer mouse Hbs are also comparable to those of other mammalian Hbs. In contrast, the high- and low-altitude variants of deer mouse Hb exhibited similarly low temperature sensitivities that were independent of red blood cell anionic cofactors, suggesting an appreciable endothermic allosteric transition upon oxygenation. In conclusion, high-altitude deer mice have evolved an adaptive increase in Hb-O2 affinity, but this is not associated with compensatory changes in sensitivity to changes in pH or temperature. Instead, it appears that the elevated Hb-O2 affinity in high-altitude deer mice is compensated by an associated increase in the tissue diffusion capacity of O2 (via increased muscle capillarization), which promotes O2 unloading.

Functional Genomic Insights into Regulatory Mechanisms of High-Altitude Adaptation

Recent studies of indigenous human populations at high altitude have provided proof-of-principle that genome scans of DNA polymorphism can be used to identify candidate loci for hypoxia adaptation. When integrated with experimental analyses of physiological phenotypes, genome-wide surveys of DNA polymorphism and tissue-specific transcriptional profiles can provide insights into actual mechanisms of adaptation. It has been suggested that adaptive phenotypic evolution is largely mediated by cis-regulatory changes in genes that are located at integrative control points in regulatory networks. This hypothesis can be tested by conducting transcriptomic analyses of hypoxic signaling pathways in conjunction with experimental measures of vascular oxygen supply and metabolic pathway flux. Such studies may reveal whether the architecture of gene regulatory networks can be used to predict which loci (and which types of loci) are likely to be "hot spots" for adaptive physiological evolution. Functional genomic studies of deer mice (Peromyscus maniculatus) demonstrate how the integrated analysis of variation in tissue-specific transcriptomes, whole-animal physiological performance, and various subordinate traits can yield insights into the mechanistic underpinnings of high-altitude adaptation.

Using the variability of linkage disequilibrium between subpopulations to infer sweeps and epistatic selection in a diverse panel of chickens

A whole-genome scan for identifying selection acting on pairs of linked loci is proposed and implemented. The scan is based on , one of Ohta's 1982 measures of between-population linkage disequilibrium (LD). An approximate empirical null distribution for the statistic is suggested. Although the partitioning of LD into between-population components was originally used to investigate epistatic selection, we demonstrate that values of may also be influenced by single-locus selective sweeps with linkage but no epistasis. The proposed scan is implemented in a diverse panel of chickens including 72 distinct breeds genotyped at 538 298 single-nucleotide polymorphisms. In all, 1723 locus pairs are identified as putatively corresponding to a selective sweep or epistatic selection. These pairs of loci generally cluster to form overlapping or neighboring signals of selection. Known variants that were expected to have been under selection in the panel are identified, as well as an assortment of novel regions that have putatively been under selection in chickens. Notably, a promising pair of genes located 8 MB apart on chromosome 9 are identified based on as demonstrating strong evidence of dispersive epistatic selection between populations.

Genetic approaches in comparative and evolutionary physiology

Whole animal physiological performance is highly polygenic and highly plastic, and the same is generally true for the many subordinate traits that underlie performance capacities. Quantitative genetics, therefore, provides an appropriate framework for the analysis of physiological phenotypes and can be used to infer the microevolutionary processes that have shaped patterns of trait variation within and among species. In cases where specific genes are known to contribute to variation in physiological traits, analyses of intraspecific polymorphism and interspecific divergence can reveal molecular mechanisms of functional evolution and can provide insights into the possible adaptive significance of observed sequence changes. In this review, we explain how the tools and theory of quantitative genetics, population genetics, and molecular evolution can inform our understanding of mechanism and process in physiological evolution. For example, lab-based studies of polygenic inheritance can be integrated with field-based studies of trait variation and survivorship to measure selection in the wild, thereby providing direct insights into the adaptive significance of physiological variation. Analyses of quantitative genetic variation in selection experiments can be used to probe interrelationships among traits and the genetic basis of physiological trade-offs and constraints. We review approaches for characterizing the genetic architecture of physiological traits, including linkage mapping and association mapping, and systems approaches for dissecting intermediary steps in the chain of causation between genotype and phenotype. We also discuss the promise and limitations of population genomic approaches for inferring adaptation at specific loci. We end by highlighting the role of organismal physiology in the functional synthesis of evolutionary biology.

Intraspecific polymorphism, interspecific divergence, and the origins of function-altering mutations in deer mouse hemoglobin

Major challenges for illuminating the genetic basis of phenotypic evolution are to identify causative mutations, to quantify their functional effects, to trace their origins as new or preexisting variants, and to assess the manner in which segregating variation is transduced into species differences. Here, we report an experimental analysis of genetic variation in hemoglobin (Hb) function within and among species of Peromyscus mice that are native to different elevations. A multilocus survey of sequence variation in the duplicated HBA and HBB genes in Peromyscus maniculatus revealed that function-altering amino acid variants are widely shared among geographically disparate populations from different elevations, and numerous amino acid polymorphisms are also shared with closely related species. Variation in Hb-O₂ affinity within and among populations of P. maniculatus is attributable to numerous amino acid mutations that have individually small effects. One especially surprising feature of the Hb polymorphism in P. maniculatus is that an appreciable fraction of functional standing variation in the two transcriptionally active HBA paralogs is attributable to recurrent gene conversion from a tandemly linked HBA pseudogene. Moreover, transpecific polymorphism in the duplicated HBA genes is not solely attributable to incomplete lineage sorting or introgressive hybridization; instead, it is mainly attributable to recurrent interparalog gene conversion that has occurred independently in different species. Partly as a result of concerted evolution between tandemly duplicated globin genes, the same amino acid changes that contribute to variation in Hb function within P. maniculatus also contribute to divergence in Hb function among different species of Peromyscus. In the case of function-altering Hb mutations in Peromyscus, there is no qualitative or quantitative distinction between segregating variants within species and fixed differences between species.

Detecting short spatial scale local adaptation and epistatic selection in climate-related candidate genes in European beech (Fagus sylvatica) populations

Detecting signatures of selection in tree populations threatened by climate change is currently a major research priority. Here, we investigated the signature of local adaptation over a short spatial scale using 96 European beech (Fagus sylvatica L.) individuals originating from two pairs of populations on the northern and southern slopes of Mont Ventoux (south-eastern France). We performed both single and multilocus analysis of selection based on 53 climate-related candidate genes containing 546 SNPs. FST outlier methods at the SNP level revealed a weak signal of selection, with three marginally significant outliers in the northern populations. At the gene level, considering haplotypes as alleles, two additional marginally significant outliers were detected, one on each slope. To account for the uncertainty of haplotype inference, we averaged the Bayes factors over many possible phase reconstructions. Epistatic selection offers a realistic multilocus model of selection in natural populations. Here, we used a test suggested by Ohta based on the decomposition of the variance of linkage disequilibrium. Overall populations, 0.23% of the SNP pairs (haplotypes) showed evidence of epistatic selection, with nearly 80% of them being within genes. One of the between gene epistatic selection signals arose between an FST outlier and a nonsynonymous mutation in a drought response gene. Additionally, we identified haplotypes containing selectively advantageous allele combinations which were unique to high or low elevations and northern or southern populations. Several haplotypes contained nonsynonymous mutations situated in genes with known functional importance for adaptation to climatic factors.

AFLPs and mitochondrial haplotypes reveal local adaptation to extreme thermal environments in a freshwater gastropod

The way environmental variation shapes neutral and adaptive genetic variation in natural populations is a key issue in evolutionary biology. Genome scans allow the identification of the genetic basis of local adaptation without previous knowledge of genetic variation or traits under selection. Candidate loci for divergent adaptation are expected to show higher F_ST than neutral loci influenced solely by random genetic drift, migration and mutation. The comparison of spatial patterns of neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection among populations living in contrasting environments. Using the gastropod Radix balthica as a system, we analyzed 376 AFLP markers and 25 mtDNA COI haplotypes for candidate loci and associations with local adaptation among contrasting thermal environments in Lake Mývatn, a volcanic lake in northern Iceland. We found that 2% of the analysed AFLP markers were under directional selection and 12% of the mitochondrial haplotypes correlated with differing thermal habitats. The genetic networks were concordant for AFLP markers and mitochondrial haplotypes, depicting distinct topologies at neutral and candidate loci. Neutral topologies were characterized by intense gene flow revealed by dense nets with edges connecting contrasting thermal habitats, whereas the connections at candidate loci were mostly restricted to populations within each thermal habitat and the number of edges decreased with temperature. Our results suggest microgeographic adaptation within Lake Mývatn and highlight the utility of genome scans in detecting adaptive divergence.

Ecological genomics of local adaptation

It is increasingly important to improve our understanding of the genetic basis of local adaptation because of its relevance to climate change, crop and animal production, and conservation of genetic resources. Phenotypic patterns that are generated by spatially varying selection have long been observed, and both genetic mapping and field experiments provided initial insights into the genetic architecture of adaptive traits. Genomic tools are now allowing genome-wide studies, and recent theoretical advances can help to design research strategies that combine genomics and field experiments to examine the genetics of local adaptation. These advances are also allowing research in non-model species, the adaptation patterns of which may differ from those of traditional model species.

Epistasis among adaptive mutations in deer mouse hemoglobin

Epistatic interactions between mutant sites in the same protein can exert a strong influence on pathways of molecular evolution. We performed protein engineering experiments that revealed pervasive epistasis among segregating amino acid variants that contribute to adaptive functional variation in deer mouse hemoglobin (Hb). Amino acid mutations increased or decreased Hb-O2 affinity depending on the allelic state of other sites. Structural analysis revealed that epistasis for Hb-O2 affinity and allosteric regulatory control is attributable to indirect interactions between structurally remote sites. The prevalence of sign epistasis for fitness-related biochemical phenotypes has important implications for the evolutionary dynamics of protein polymorphism in natural populations.

Molecular divergence in tropical tree populations occupying environmental mosaics

Unveiling the genetic basis of local adaptation to environmental variation is a major goal in molecular ecology. In rugged landscapes characterized by environmental mosaics, living populations and communities can experience steep ecological gradients over very short geographical distances. In lowland tropical forests, interspecific divergence in edaphic specialization (for seasonally flooded bottomlands and seasonally dry terra firme soils) has been proven by ecological studies on adaptive traits. Some species are nevertheless capable of covering the entire span of the gradient; intraspecific variation for adaptation to contrasting conditions may explain the distribution of such ecological generalists. We investigated whether local divergence happens at small spatial scales in two stands of Eperua falcata (Fabaceae), a widespread tree species of the Guiana Shield. We investigated Single Nucleotide Polymorphisms (SNP) and sequence divergence as well as spatial genetic structure (SGS) at four genes putatively involved in stress response and three genes with unknown function. Significant genetic differentiation was observed among sub-populations within stands, and eight SNP loci showed patterns compatible with disruptive selection. SGS analysis showed genetic turnover along the gradients at three loci, and at least one haplotype was found to be in repulsion with one habitat. Taken together, these results suggest genetic differentiation at small spatial scale in spite of gene flow. We hypothesize that heterogeneous environments may cause molecular divergence, possibly associated to local adaptation in E. falcata.

Phenotypic plasticity in blood-oxygen transport in highland and lowland deer mice

In vertebrates living at high altitude, arterial hypoxemia may be ameliorated by reversible changes in the oxygen-carrying capacity of the blood (regulated by erythropoiesis) and/or changes in blood-oxygen affinity (regulated by allosteric effectors of hemoglobin function). These hematological traits often differ between taxa that are native to different elevational zones, but it is often unknown whether the observed physiological differences reflect fixed, genetically based differences or environmentally induced acclimatization responses (phenotypic plasticity). Here, we report measurements of hematological traits related to blood-O2 transport in populations of deer mice (Peromyscus maniculatus) that are native to high- and low-altitude environments. We conducted a common-garden breeding experiment to assess whether altitude-related physiological differences were attributable to developmental plasticity and/or physiological plasticity during adulthood. Under conditions prevailing in their native habitats, high-altitude deer mice from the Rocky Mountains exhibited a number of pronounced hematological differences relative to low-altitude conspecifics from the Great Plains: higher hemoglobin concentrations, higher hematocrits, higher erythrocytic concentrations of 2,3-diphosphoglycerate (an allosteric regulator of hemoglobin-oxygen affinity), lower mean corpuscular hemoglobin concentrations and smaller red blood cells. However, these differences disappeared after 6 weeks of acclimation to normoxia at low altitude. The measured traits were also indistinguishable between the F1 progeny of highland and lowland mice, indicating that there were no persistent differences in phenotype that could be attributed to developmental plasticity. These results indicate that the naturally occurring hematological differences between highland and lowland mice are environmentally induced and are largely attributable to physiological plasticity during adulthood.

Extensive linkage disequilibrium and parallel adaptive divergence across threespine stickleback genomes

Population genomic studies are beginning to provide a more comprehensive view of dynamic genome-scale processes in evolution. Patterns of genomic architecture, such as genomic islands of increased divergence, may be important for adaptive population differentiation and speciation. We used next-generation sequencing data to examine the patterns of local and long-distance linkage disequilibrium (LD) across oceanic and freshwater populations of threespine stickleback, a useful model for studies of evolution and speciation. We looked for associations between LD and signatures of divergent selection, and assessed the role of recombination rate variation in generating LD patterns. As predicted under the traditional biogeographic model of unidirectional gene flow from ancestral oceanic to derived freshwater stickleback populations, we found extensive local and long-distance LD in fresh water. Surprisingly, oceanic populations showed similar patterns of elevated LD, notably between large genomic regions previously implicated in adaptation to fresh water. These results support an alternative biogeographic model for the stickleback radiation, one of a metapopulation with appreciable bi-directional gene flow combined with strong divergent selection between oceanic and freshwater populations. As predicted by theory, these processes can maintain LD within and among genomic islands of divergence. These findings suggest that the genomic architecture in oceanic stickleback populations may provide a mechanism for the rapid re-assembly and evolution of multi-locus genotypes in newly colonized freshwater habitats, and may help explain genetic mapping of parallel phenotypic variation to similar loci across independent freshwater populations.

Genomic insights into adaptation to high-altitude environments

Elucidating the molecular genetic basis of adaptive traits is a central goal of evolutionary genetics. The cold, hypoxic conditions of high-altitude habitats impose severe metabolic demands on endothermic vertebrates, and understanding how high-altitude endotherms cope with the combined effects of hypoxia and cold can provide important insights into the process of adaptive evolution. The physiological responses to high-altitude stress have been the subject of over a century of research, and recent advances in genomic technologies have opened up exciting opportunities to explore the molecular genetic basis of adaptive physiological traits. Here, we review recent literature on the use of genomic approaches to study adaptation to high-altitude hypoxia in terrestrial vertebrates, and explore opportunities provided by newly developed technologies to address unanswered questions in high-altitude adaptation at a genomic scale.

Estimating population divergence time and phylogeny from single-nucleotide polymorphisms data with outgroup ascertainment bias

The inference of population divergence times and branching patterns is of fundamental importance in many population genetic analyses. Many methods have been developed for estimating population divergence times, and recently, there has been particular attention towards genome-wide single-nucleotide polymorphisms (SNP) data. However, most SNP data have been affected by an ascertainment bias caused by the SNP selection and discovery protocols. Here, we present a modification of an existing maximum likelihood method that will allow approximately unbiased inferences when ascertainment is based on a set of outgroup populations. We also present a method for estimating trees from the asymmetric dissimilarity measures arising from pairwise divergence time estimation in population genetics. We evaluate the methods by simulations and by applying them to a large SNP data set of seven East Asian populations.

BAC-based sequencing of behaviorally-relevant genes in the prairie vole

The prairie vole (Microtus ochrogaster) is an important model organism for the study of social behavior, yet our ability to correlate genes and behavior in this species has been limited due to a lack of genetic and genomic resources. Here we report the BAC-based targeted sequencing of behaviorally-relevant genes and flanking regions in the prairie vole. A total of 6.4 Mb of non-redundant or haplotype-specific sequence assemblies were generated that span the partial or complete sequence of 21 behaviorally-relevant genes as well as an additional 55 flanking genes. Estimates of nucleotide diversity from 13 loci based on alignments of 1.7 Mb of haplotype-specific assemblies revealed an average pair-wise heterozygosity (8.4×10⁻³). Comparative analyses of the prairie vole proteins encoded by the behaviorally-relevant genes identified >100 substitutions specific to the prairie vole lineage. Finally, our sequencing data indicate that a duplication of the prairie vole AVPR1A locus likely originated from a recent segmental duplication spanning a minimum of 105 kb. In summary, the results of our study provide the genomic resources necessary for the molecular and genetic characterization of a high-priority set of candidate genes for regulating social behavior in the prairie vole.

Population genomic analysis of model and nonmodel organisms using sequenced RAD tags

The evolutionary processes of mutation, migration, genetic drift, and natural selection shape patterns of genetic variation among individuals, populations, and species, and they can do so differentially across genomes. The field of population genomics provides a comprehensive genome-scale view of these processes, even beyond traditional model organisms. Until recently, genome-wide studies of genetic variation have been prohibitively expensive. However, next-generation sequencing (NGS) technologies are revolutionizing the field of population genomics, allowing for genetic analysis at scales not previously possible even in organisms for which few genomic resources presently exist. To speed this revolution in evolutionary genetics, we and colleagues developed Restriction site Associated DNA (RAD) sequencing, a method that uses Illumina NGS to simultaneously type and score tens to hundreds of thousands of single nucleotide polymorphism (SNP) markers in hundreds of individuals for minimal investment of resources. The core molecular protocol is described elsewhere in this volume, which can be modified to suit a diversity of evolutionary genetic questions. In this chapter, we outline the conceptual framework of population genomics, relate genomic patterns of variation to evolutionary processes, and discuss how RAD sequencing can be used to study population genomics. In addition, we discuss bioinformatic considerations that arise from unique aspects of NGS data as compared to traditional marker based approaches, and we outline some general analytical approaches for RAD-seq and similar data, including a computational pipeline that we developed called Stacks. This software can be used for the analysis of RAD-seq data in organisms with and without a reference genome. Nonetheless, the development of analytical tools remains in its infancy, and further work is needed to fully quantify sampling variance and biases in these data types. As data-gathering technology continues to advance, our ability to understand genomic evolution in natural populations will be limited more by conceptual and analytical weaknesses than by the amount of molecular data.

Experimental approaches to evaluate the contributions of candidate protein-coding mutations to phenotypic evolution

Identifying mechanisms of molecular adaptation can provide important insights into the process of phenotypic evolution, but it can be exceedingly difficult to quantify the phenotypic effects of specific mutational changes. To verify the adaptive significance of genetically based changes in protein function, it is necessary to document functional differences between the products of derived and wild-type alleles and to demonstrate that such differences impinge on higher-level physiological processes (and ultimately, fitness). In the case of metabolic enzymes, this requires documenting in vivo differences in reaction rate that give rise to differences in flux through the pathway in which the enzymes function. These measured differences in pathway flux should then give rise to differences in cellular or systemic physiology that affect fitness-related variation in whole-organism performance. Efforts to establish these causal connections between genotype, phenotype, and fitness require experiments that carefully control for environmental variation and background genetic variation. Here, we discuss experimental approaches to evaluate the contributions of amino-acid mutations to adaptive phenotypic change. We discuss conceptual and methodological issues associated with in vitro and in vivo studies of protein function, and the evolutionary insights that can be gleaned from such studies. We also discuss the importance of isolating the effects of individual mutations to distinguish between positively selected substitutions that directly contribute to improvements in protein function versus positively selected, compensatory substitutions that mitigate negative pleiotropic effects of antecedent changes.

Multilocus coalescent analysis of haemoglobin differentiation between low- and high-altitude populations of crested ducks (Lophonetta specularioides)

Hypoxia is a key factor determining survival, and haemoglobins are targets of selection in species native to high-altitude regions. We studied population genetic structure and evaluated evidence for local adaptation in the crested duck (Lophonetta specularioides). Differentiation, gene flow and time since divergence between highland and lowland populations were assessed for three haemoglobin genes (α(A) , α(D) , β(A) ) and compared to seven reference loci (six autosomal introns and mtDNA). Four derived amino acid replacements were found in the globin genes that had elevated Φ(ST) values between the Andean highlands and Patagonian lowlands. A single β(A) -globin polymorphism at a site known to influence O(2) affinity was fixed for different alleles in the two populations, whereas three α(A) - and α(D) -globin polymorphisms exhibited high heterozygosity in the highlands but not in the lowlands. Coalescent analyses supported restricted gene flow for haemoglobin alleles and mitochondrial DNA but nonzero gene flow for the introns. Simulating genetic data under a drift-migration model of selective neutrality, the β(A) -globin fell outside the 95% confidence limit of simulated data, suggesting that directional selection is maintaining different variants in the contrasting elevational environments, thereby restricting migration of β(A) -globin alleles. The α(A) - and α(D) -globins, by contrast, did not differ from the simulated values, suggesting that variants in these genes are either selectively neutral, or that the effects of selection could not be differentiated from background levels of population structure and linkage disequilibrium. This study illustrates the combined effects of selection and population history on inferring levels of population divergence for a species distributed across an altitudinal gradient in which selection for hypoxia resistance has likely played an important role.

Altitudinal variation at duplicated β-globin genes in deer mice: effects of selection, recombination, and gene conversion

Spatially varying selection on a given polymorphism is expected to produce a localized peak in the between-population component of nucleotide diversity, and theory suggests that the chromosomal extent of elevated differentiation may be enhanced in cases where tandemly linked genes contribute to fitness variation. An intriguing example is provided by the tandemly duplicated β-globin genes of deer mice (Peromyscus maniculatus), which contribute to adaptive differentiation in blood-oxygen affinity between high- and low-altitude populations. Remarkably, the two β-globin genes segregate the same pair of functionally distinct alleles due to a history of interparalog gene conversion and alleles of the same functional type are in perfect coupling-phase linkage disequilibrium (LD). Here we report a multilocus analysis of nucleotide polymorphism and LD in highland and lowland mice with different genetic backgrounds at the β-globin genes. The analysis of haplotype structure revealed a paradoxical pattern whereby perfect LD between the two β-globin paralogs (which are separated by 16.2 kb) is maintained in spite of the fact that LD within both paralogs decays to background levels over physical distances of less than 1 kb. The survey of nucleotide polymorphism revealed that elevated levels of altitudinal differentiation at each of the β-globin genes drop away quite rapidly in the external flanking regions (upstream of the 5' paralog and downstream of the 3' paralog), but the level of differentiation remains unexpectedly high across the intergenic region. Observed patterns of diversity and haplotype structure are difficult to reconcile with expectations of a two-locus selection model with multiplicative fitness.

Molecular spandrels: tests of adaptation at the genetic level

Although much progress has been made in identifying the genes (and, in rare cases, mutations) that contribute to phenotypic variation, less is known about the effects that these genes have on fitness. Nonetheless, genes are commonly labelled as 'adaptive' if an allele has been shown to affect a phenotype with known or suspected functional importance or if patterns of nucleotide variation at the locus are consistent with positive selection. In these cases, the 'adaptive' designation may be premature and may lead to incorrect conclusions about the relationships between gene function and fitness. Experiments to test targets and agents of natural selection within a genomic context are necessary for identifying the adaptive consequences of individual alleles.

Molecular signatures of selection on reproductive character displacement of flower color in Phlox drummondii

Character displacement, which arises when species diverge in sympatry to decrease competition for resources or reproductive interference, has been observed in a wide variety of plants and animals. A classic example of reproductive character displacement, presumed to be caused by reinforcing selection, is flower-color variation in the native Texas wildflower Phlox drummondii. Here, we use population genetic analyses to investigate molecular signatures of selection on flower-color variation in this species. First, we quantify patterns of neutral genetic variation across the range of P. drummondii to demonstrate that restricted gene flow and genetic drift cannot explain the pattern of flower-color divergence in this species. There is evidence of extensive gene flow across populations with different flower colors, suggesting selection caused flower-color divergence. Second, analysis of sequence variation in the genes underlying this divergence reveals a signature of a selective sweep in one of the two genes, further indicating selection is responsible for divergence in sympatry. The lack of a signature of selection at the second locus does not necessarily indicate a lack of selection on this locus but instead brings attention to the uncertainty in depending on molecular signatures to identify selection.

Local selection across a latitudinal gradient shapes nucleotide diversity in balsam poplar, Populus balsamifera L

Molecular studies of adaptive evolution often focus on detecting selective sweeps driven by positive selection on a species-wide scale; however, much adaptation is local, particularly of ecologically important traits. Here, we look for evidence of range-wide and local adaptation at candidate genes for adaptive phenology in balsam poplar, Populus balsamifera, a widespread forest tree whose range extends across environmental gradients of photoperiod and growing season length. We examined nucleotide diversity of 27 poplar homologs of the flowering-time network-a group of genes that control plant developmental phenology through interactions with environmental cues such as photoperiod and temperature. Only one gene, ZTL2, showed evidence of reduced diversity and an excess of fixed replacement sites, consistent with a species-wide selective sweep. Two other genes, LFY and FRI, harbored high levels of nucleotide diversity and exhibited elevated differentiation between northern and southern accessions, suggesting local adaptation along a latitudinal gradient. Interestingly, FRI has also been identified as a target of local selection between northern and southern accessions of Arabidopsis thaliana, indicating that this gene may be commonly involved in ecological adaptation in distantly related species. Our findings suggest an important role for local selection shaping molecular diversity and reveal limitations of inferring molecular adaptation from analyses designed only to detect species-wide selective sweeps.

Expression and purification of recombinant hemoglobin in Escherichia coli

Background

Recombinant DNA technologies have played a pivotal role in the elucidation of structure-function relationships in hemoglobin (Hb) and other globin proteins. Here we describe the development of a plasmid expression system to synthesize recombinant Hbs in Escherichia coli, and we describe a protocol for expressing Hbs with low intrinsic solubilities. Since the α- and β-chain Hbs of different species span a broad range of solubilities, experimental protocols that have been optimized for expressing recombinant human HbA may often prove unsuitable for the recombinant expression of wildtype and mutant Hbs of other species.

Methodology/Principal Findings

As a test case for our expression system, we produced recombinant Hbs of the deer mouse (Peromyscus maniculatus), a species that has been the subject of research on mechanisms of Hb adaptation to hypoxia. By experimentally assessing the combined effects of induction temperature, induction time and E. coli expression strain on the solubility of recombinant deer mouse Hbs, we identified combinations of expression conditions that greatly enhanced the yield of recombinant protein and which also increased the efficiency of post-translational modifications.

Conclusion/Significance

Our protocol should prove useful for the experimental study of recombinant Hbs in many non-human animals. One of the chief advantages of our protocol is that we can express soluble recombinant Hb without co-expressing molecular chaperones, and without the need for additional reconstitution or heme-incorporation steps. Moreover, our plasmid construct contains a combination of unique restriction sites that allows us to produce recombinant Hbs with different α- and β-chain subunit combinations by means of cassette mutagenesis.

Genetic evaluation of eri silkworm Samia cynthia ricini: ISSR loci specific to high and low altitude regimes and quantitative attributes

Genetic structure of populations is under constant pressure from varying geographical conditions that induce phenotypic plasticity in insects. Spatial distribution of 15 populations of Indian eri silkworm, Samia cynthia ricini originated at various altitudes of sub Himalayas based on Euclidean distance realized from yield attributes showed two population clusters irrespective of their place of origin and altitude. However, DNA amplification profile by inter SSR (ISSR) markers showed genetic variations among the populations depend on low and high altitudes. One ISSR locus each specific to high and low altitude population was identified. The locus from high altitude showed deviation from Hardy-Weinberg equilibrium but that from low altitude was in neutrality suggests that the high altitude loci could be under pressure from the altitudinal variations. In association with different yield traits, 18 loci were identified. Of which, three markers showed association with more than one trait indicative of pleiotropic influence. Stepwise addition of markers enhanced the correlation between markers and the associated trait pointed to polygenic influence. Association of markers with altitude and yield traits suggests an imperative relation of rare genetic loci with gene-environment interaction and phenotypic variability in S. c. ricini.

Integrating evolutionary and functional approaches to infer adaptation at specific loci

Inferences about adaptation at specific loci are often exclusively based on the static analysis of DNA sequence variation. Ideally,population-genetic evidence for positive selection serves as a stepping-off point for experimental studies to elucidate the functional significance of the putatively adaptive variation. We argue that inferences about adaptation at specific loci are best achieved by integrating the indirect, retrospective insights provided by population-genetic analyses with the more direct, mechanistic insights provided by functional experiments. Integrative studies of adaptive genetic variation may sometimes be motivated by experimental insights into molecular function, which then provide the impetus to perform population genetic tests to evaluate whether the functional variation is of adaptive significance. In other cases, studies may be initiated by genome scans of DNA variation to identify candidate loci for recent adaptation. Results of such analyses can then motivate experimental efforts to test whether the identified candidate loci do in fact contribute to functional variation in some fitness-related phenotype. Functional studies can provide corroborative evidence for positive selection at particular loci, and can potentially reveal specific molecular mechanisms of adaptation.

Enhanced AFLP genome scans detect local adaptation in high-altitude populations of a small rodent (Microtus arvalis)

Adaptation to adverse environmental conditions such as high altitude requires physiological and/or morphological changes. Genome scans provide a means to identify the genetic basis of such adaptations without previous knowledge about the particular genetic variants or traits under selection. In this study, we scanned 3027 amplified fragment length polymorphisms (AFLP) in four populations of the common vole Microtus arvalis for loci associated with local adaptation and high altitude. We investigated voles from two populations at high elevation (~2000 m a.s.l.) representing the upper limit of the altitudinal distribution of the species and two geographically close low-altitude populations (<600 m a.s.l.). Statistical analysis incorporated a new Bayesian F(ST) outlier approach specifically developed for AFLP markers, which considers the intensity of AFLP bands instead of mere presence/absence and allows to derive population-based estimates of allele frequencies and F(IS) values. Computer simulations showed that this approach increases the statistical power of the detection of AFLP markers under selection almost to the power of single nucleotide polymorphism (SNP) data without compromising specificity. Our enhanced genome scan resulted in 20 prime candidate markers for positive selection, which show mostly extremely high allele frequency differences between the low- and high-altitude populations. The comparison of global- and pairwise-enhanced genome scans demonstrated further that very strong selective signatures may also be associated with single populations suggesting the importance of local adaptation in alpine populations of common voles.

A tale of two genomes: contrasting patterns of phylogeographic structure in a widely distributed bat

One of the most widely distributed bats in the New World, the big brown bat (Eptesicus fuscus) exhibits well-documented geographic variation in morphology and life history traits, suggesting the potential for significant phylogeographic structure as well as adaptive differentiation among populations. In a pattern broadly consistent with morphologically defined subspecies, we found deeply divergent mitochondrial lineages restricted to different geographic regions. In contrast, sequence data from two nuclear loci suggest a general lack of regional genetic structure except for peripheral populations in the Caribbean and Mexico/South America. Coalescent analyses suggest that the striking difference in population structure between genomes cannot be attributed solely to different rates of lineage sorting, but is likely due to male-mediated gene flow homogenizing nuclear genetic diversity across most of the continental range. Despite this ongoing gene flow, selection has apparently been effective in producing and maintaining adaptive differentiation among populations, while strong female site fidelity, maintained over the course of millions of years, has produced remarkably deep divergence among geographically isolated matrilines. Our results highlight the importance of evaluating multiple genetic markers for a more complete understanding of population structure and history.