A dynamic history of admixture from Mediterranean and Carpathian glacial refugia drives genomic diversity in the bank vole

Understanding the historical contributions of differing glacial refugia is key to evaluating the roles of microevolutionary forces, such as isolation, introgression, and selection in shaping genomic diversity in present‐day populations. In Europe, where both Mediterranean and extra‐Mediterranean (e.g., Carpathian) refugia of the bank vole (Clethrionomys glareolus) have been identified, mtDNA indicates that extra‐Mediterranean refugia were the main source of colonization across the species range, while Mediterranean peninsulas harbor isolated, endemic lineages. Here, we critically evaluate this hypothesis using previously generated genomic data (>6,000 SNPs) for over 800 voles, focusing on genomic contributions to bank voles in central Europe, a key geographic area in considering range‐wide colonization. The results provide clear evidence that both extra‐Mediterranean (Carpathian) and Mediterranean (Spanish, Calabrian, and Balkan) refugia contributed to the ancestry and genomic diversity of bank vole populations across Europe. Few strong barriers to dispersal and frequent admixture events in central Europe have led to a prominent mid‐latitude peak in genomic diversity. Although the genomic contribution of the centrally located Carpathian refugium predominates, populations in different parts of Europe have admixed origins from Mediterranean (28%–47%) and the Carpathian (53%–72%) sources. We suggest that the admixture from Mediterranean refugia may have provisioned adaptive southern alleles to more northern populations, facilitating the end‐glacial spread of the admixed populations and contributing to increased bank vole diversity in central Europe. This study adds critical details to the complex end‐glacial colonization history of this well‐studied organism and underscores the importance of genomic data in phylogeographic interpretation.     

Unraveling the Effects of Selection and Demography on Immune Gene Variation in Free-Ranging Plains Zebra (Equus quagga) Populations

Demography, migration and natural selection are predominant processes affecting the distribution of genetic variation among natural populations. Many studies use neutral genetic markers to make inferences about population history. However, the investigation of functional coding loci, which directly reflect fitness, is critical to our understanding of species'' ecology and evolution. Immune genes, such as those of the Major Histocompatibility Complex (MHC), play an important role in pathogen recognition and provide a potent model system for studying selection. We contrasted diversity patterns of neutral data with MHC loci, ELA-DRA and -DQA, in two southern African plains zebra (Equus quagga) populations: Etosha National Park, Namibia, and Kruger National Park, South Africa. Results from neutrality tests, along with observations of elevated diversity and low differentiation across populations, supported previous genus-level evidence for balancing selection at these loci. Despite being low, MHC divergence across populations was significant and may be attributed to drift effects typical of geographically separated populations experiencing little to no gene flow, or alternatively to shifting allele frequency distributions driven by spatially variable and fluctuating pathogen communities. At the DRA, zebra exhibited geographic differentiation concordant with microsatellites and reduced levels of diversity in Etosha due to highly skewed allele frequencies that could not be explained by demography, suggestive of spatially heterogeneous selection and local adaptation. This study highlights the complexity in which selection affects immune gene diversity and warrants the need for further research on the ecological mechanisms shaping patterns of adaptive variation among natural populations.     

The population genomic signature of environmental selection in the widespread insect-pollinated tree species Frangula alnus at different geographical scales

The evaluation of the molecular signatures of selection in species lacking an available closely related reference genome remains challenging, yet it may provide valuable fundamental insights into the capacity of populations to respond to environmental cues. We screened 25 native populations of the tree species Frangula alnus subsp. alnus (Rhamnaceae), covering three different geographical scales, for 183 annotated single-nucleotide polymorphisms (SNPs). Standard population genomic outlier screens were combined with individual-based and multivariate landscape genomic approaches to examine the strength of selection relative to neutral processes in shaping genomic variation, and to identify the main environmental agents driving selection. Our results demonstrate a more distinct signature of selection with increasing geographical distance, as indicated by the proportion of SNPs (i) showing exceptional patterns of genetic diversity and differentiation (outliers) and (ii) associated with climate. Both temperature and precipitation have an important role as selective agents in shaping adaptive genomic differentiation in F. alnus subsp. alnus, although their relative importance differed among spatial scales. At the ‘intermediate'' and ‘regional'' scales, where limited genetic clustering and high population diversity were observed, some indications of natural selection may suggest a major role for gene flow in safeguarding adaptability. High genetic diversity at loci under selection in particular, indicated considerable adaptive potential, which may nevertheless be compromised by the combined effects of climate change and habitat fragmentation.     

Selective pressures on MHC class II genes in the guppy (Poecilia reticulata) as inferred by hierarchical analysis of population structure

Genes of the major histocompatibility complex, which are the most polymorphic of all vertebrate genes, are a pre‐eminent system for the study of selective pressures that arise from host–pathogen interactions. Balancing selection capable of maintaining high polymorphism should lead to the homogenization of MHC allele frequencies among populations, but there is some evidence to suggest that diversifying selection also operates on the MHC. However, the pattern of population structure observed at MHC loci is likely to depend on the spatial and/or temporal scale examined. Here, we investigated selection acting on MHC genes at different geographic scales using Venezuelan guppy populations inhabiting four regions. We found a significant correlation between MHC and microsatellite allelic richness across populations, which suggests the role of genetic drift in shaping MHC diversity. However, compared to microsatellites, more MHC variation was explained by differences between populations within larger geographic regions and less by the differences between the regions. Furthermore, among proximate populations, variation in MHC allele frequencies was significantly higher compared to microsatellites, indicating that selection acting on MHC may increase population structure at small spatial scales. However, in populations that have significantly diverged at neutral markers, the population‐genetic signature of diversifying selection may be eradicated in the long term by that of balancing selection, which acts to preserve rare alleles and thus maintain a common pool of MHC alleles.     

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.     

Natural Selection and Neutral Evolution Jointly Drive Population Divergence between Alpine and Lowland Ecotypes of the Allopolyploid Plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074–0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (F_ST  = 0.041–0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.     

Characterisation of Adaptive Genetic Diversity in Environmentally Contrasted Populations of Eucalyptus camaldulensis Dehnh. (River Red Gum)

As an increasing number of ecosystems face departures from long standing environmental conditions under climate change, our understanding of the capacity of species to adapt will become important for directing conservation and management of biodiversity. Insights into the potential for genetic adaptation might be gained by assessing genomic signatures of adaptation to historic or prevailing environmental conditions. The river red gum (Eucalyptus camaldulensis Dehnh.) is a widespread Australian eucalypt inhabiting riverine and floodplain habitats which spans strong environmental gradients. We investigated the effects of adaptation to environment on population level genetic diversity of E. camaldulensis, examining SNP variation in candidate gene loci sampled across 20 climatically diverse populations approximating the species natural distribution. Genetic differentiation among populations was high (F_ST = 17%), exceeding previous estimates based on neutral markers. Complementary statistical approaches identified 6 SNP loci in four genes (COMT, Dehydrin, ERECTA and PIP2) which, after accounting for demographic effects, exhibited higher than expected levels of genetic differentiation among populations and whose allelic variation was associated with local environment. While this study employs but a small proportion of available diversity in the eucalyptus genome, it draws our attention to the potential for application of wide spread eucalypt species to test adaptive hypotheses.     

Balancing selection versus allele and supertype turnover in MHC class II genes in guppies

Selection pressure from parasites is thought to be a major force shaping the extreme polymorphism of the major histocompatibility complex (MHC) genes, but the modes and consequences of selection remain unclear. Here, we analyse MHC class II and microsatellite diversity in 16 guppy (Poecilia reticulata) populations from two islands (Trinidad and Tobago) that have been separated for at least 10 ky. Within-population MHC diversity was high, but allele sharing was limited within islands and even lower between islands, suggesting relatively fast turnover of alleles. Allelic lineages strongly supported in phylogenetic analyses tended to be island-specific, suggesting rapid lineage sorting, and an expansion of an allelic lineage private to Tobago was observed. New alleles appear to be generated locally at a detectably high frequency. We did not detect a consistent signature of local adaptation, but F_ST outlier analysis suggested that balancing selection may be the more general process behind spatial variation in MHC allele frequencies in this system, particularly within Trinidad. We found no evidence for divergent allele advantage within populations, or for decreased genetic structuring of MHC supertypes compared to MHC alleles. The dynamic and complex nature of MHC evolution we observed in guppies, coupled with some evidence for balancing selection shaping MHC allele frequencies, are consistent with Red Queen processes of host-parasite coevolution.Subject terms: Population genetics, Evolutionary genetics     

High MHC DQB Variation and Asymmetric Allelic Distribution in the Endangered Yangtze Finless Porpoise, Neophocaena phocaenoides asiaeorientalis

The endangered Yangtze finless porpoise is found in the middle and lower reaches of the Yangtze River and its adjoining big lakes. To explore the major histocompatibility complex (MHC) genetic diversity and allelic distribution patterns across its range, we investigated variation at DQB exon 2. From 76 porpoises, we identified 18 DQB sequences. The freshwater Yangtze populations had much higher allelic diversity than marine populations. Among these freshwater populations, the middle-reach population had higher allelic diversity than the lower-reach population. The high DQB diversity level, relative to that of a neutral mtDNA locus, suggests that balancing selection is acting at the DQB gene and that rapid evolution and local positive selection play critical roles in generating and retaining high MHC diversity in the freshwater population. As the balancing selection might be driven by environmental pathogens, we suggest that maintaining MHC variation should be a high priority in the conservation and management of this endangered population, especially as an ex situ conservation strategy.     

Environmental and climatic determinants of molecular diversity and genetic population structure in a coenagrionid damselfly

Identifying environmental factors that structure intraspecific genetic diversity is of interest for both habitat preservation and biodiversity conservation. Recent advances in statistical and geographical genetics make it possible to investigate how environmental factors affect geographic organisation and population structure of molecular genetic diversity within species. Here we present a study on a common and wide ranging insect, the blue tailed damselfly Ischnuraelegans, which has been the target of many ecological and evolutionary studies. We addressed the following questions: (i) Is the population structure affected by longitudinal or latitudinal gradients?; (ii) Do geographic boundaries limit gene flow?; (iii) Does geographic distance affect connectivity and is there a signature of past bottlenecks?; (iv) Is there evidence of a recent range expansion and (vi) what is the effect of geography and climatic factors on population structure? We found low to moderate genetic sub-structuring between populations (mean F_ST = 0.06, D_est = 0.12), and an effect of longitude, but not latitude, on genetic diversity. No significant effects of geographic boundaries (e.g. water bodies) were found. F_ST-and D_est-values increased with geographic distance; however, there was no evidence for recent bottlenecks. Finally, we did not detect any molecular signatures of range expansions or an effect of geographic suitability, although local precipitation had a strong effect on genetic differentiation. The population structure of this small insect has probably been shaped by ecological factors that are correlated with longitudinal gradients, geographic distances, and local precipitation. The relatively weak global population structure and high degree of genetic variation within populations suggest that I. elegans has high dispersal ability, which is consistent with this species being an effective and early coloniser of new habitats.     

Neutral and adaptive genomic signatures of rapid poleward range expansion

Many species are expanding their range polewards, and this has been associated with rapid phenotypic change. Yet, it is unclear to what extent this reflects rapid genetic adaptation or neutral processes associated with range expansion, or selection linked to the new thermal conditions encountered. To disentangle these alternatives, we studied the genomic signature of range expansion in the damselfly Coenagrion scitulum using 4950 newly developed genomic SNPs and linked this to the rapidly evolved phenotypic differences between core and (newly established) edge populations. Most edge populations were genetically clearly differentiated from the core populations and all were differentiated from each other indicating independent range expansion events. In addition, evidence for genetic drift in the edge populations, and strong evidence for adaptive genetic variation in association with the range expansion was detected. We identified one SNP under consistent selection in four of the five edge populations and showed that the allele increasing in frequency is associated with increased flight performance. This indicates collateral, non‐neutral evolutionary changes in independent edge populations driven by the range expansion process. We also detected a genomic signature of adaptation to the newly encountered thermal regimes, reflecting a pattern of countergradient variation. The latter signature was identified at a single SNP as well as in a set of covarying SNPs using a polygenic multilocus approach to detect selection. Overall, this study highlights how a strategic geographic sampling design and the integration of genomic, phenotypic and environmental data can identify and disentangle the neutral and adaptive processes that are simultaneously operating during range expansions.     

The Population Genetics of Cultivation: Domestication of a Traditional Chinese Medicine,Scrophularia ningpoensis Hemsl. (Scrophulariaceae)

Background

Domestic cultivation of medicinal plants is an important strategy for protecting these species from over harvesting. Some species of medicinal plants have been brought into cultivation for more than hundreds years. Concerns about severe loss of genetic diversity and sustainable cultivation can potentially limit future use of these valuable plants. Genetic studies with comprehensive sampling of multiple medicinal species by molecular markers will allow for assessment and management of these species. Here we examine the population genetic consequences of cultivation and domestication in Scrophularia ningpoensis Hemsl. We used chloroplast DNA and genomic AFLP markers to clarify not only the effects of domestication on genetic diversity, but also determine the geographic origins of cultivars and their genetic divergence from native populations. These results will allow both better management of cultivated populations, but also provide insights for crop improvement.

Results

Twenty-one cpDNA haplotypes of S. ningpoensis were identified. Wild populations contain all haplotypes, whereas only three haplotypes were found in cultivated populations with wild populations having twice the haplotype diversity of cultivated populations. Genetic differentiation between cultivated populations and wild populations was significant. Genomic AFLP markers revealed similar genetic diversity patterns. Furthermore, Structure analysis grouped all wild populations into two gene pools; two of which shared the same gene pool with cultivated S. ningpoensis. The result of Neighbor-Joining analysis was consistent with the structure analysis. In principal coordinate analysis, three cultivated populations from Zhejiang Province grouped together and were separated from other cultivated populations.

Conclusions

These results suggest that cultivated S. ningpoensis has experienced dramatic loss of genetic diversity under anthropogenic influence. We postulate that strong artificial selection for medicinal quality has resulted in genetic differentiation between cultivated and wild populations. Furthermore, it appears that wild populations in Jiangxi-Hunan area were involved in the origin of cultivated S. ningpoensis.     

Genetic divergence and isolation by thermal environment in geothermal populations of an aquatic invertebrate

Temperature is one of the most influential forces of natural selection impacting all biological levels. In the face of increasing global temperatures, studies over small geographic scales allowing investigations on the effects of gene flow are of great value for understanding thermal adaptation. Here, we investigated genetic population structure in the freshwater gastropod Radix balthica originating from contrasting thermal habitats in three areas of geothermal activity in Iceland. Snails from 32 sites were genotyped at 208 AFLP loci. Five AFLPs were identified as putatively under divergent selection in Lake Mývatn, a geothermal lake with an almost 20 °C difference in mean temperature across a distance of a few kilometres. In four of these loci, variation across all study populations was correlated with temperature. We found significant population structure in neutral markers both within and between the areas. Cluster analysis using neutral markers classified the sites mainly by geography, whereas analyses using markers under selection differentiated the sites based on temperature. Isolation by distance was stronger in the neutral than in the outlier loci. Pairwise differences based on outlier F_ST were significantly correlated with temperature at different spatial scales, even after correcting for geographic distance or neutral pairwise F_ST differences. In general, genetic variation decreased with increasing environmental temperature, possibly suggesting that natural selection had reduced the genetic diversity in the warm origin sites. Our results emphasize the influence of environmental temperature on the genetic structure of populations and suggest local thermal adaptation in these geothermal habitats.     

Fluctuating Environments Maintain Genetic Diversity through Neutral Fitness Effects and Balancing Selection

Genetic variation is the raw material upon which selection acts. The majority of environmental conditions change over time and therefore may result in variable selective effects. How temporally fluctuating environments impact the distribution of fitness effects and in turn population diversity is an unresolved question in evolutionary biology. Here, we employed continuous culturing using chemostats to establish environments that switch periodically between different nutrient limitations and compared the dynamics of selection to static conditions. We used the pooled Saccharomyces cerevisiae haploid gene deletion collection as a synthetic model for populations comprising thousands of unique genotypes. Using barcode sequencing, we find that static environments are uniquely characterized by a small number of high-fitness genotypes that rapidly dominate the population leading to dramatic decreases in genetic diversity. By contrast, fluctuating environments are enriched in genotypes with neutral fitness effects and an absence of extreme fitness genotypes contributing to the maintenance of genetic diversity. We also identified a unique class of genotypes whose frequencies oscillate sinusoidally with a period matching the environmental fluctuation. Oscillatory behavior corresponds to large differences in short-term fitness that are not observed across long timescales pointing to the importance of balancing selection in maintaining genetic diversity in fluctuating environments. Our results are consistent with a high degree of environmental specificity in the distribution of fitness effects and the combined effects of reduced and balancing selection in maintaining genetic diversity in the presence of variable selection.     

Differing Patterns of Selection and Geospatial Genetic Diversity within Two Leading Plasmodium vivax Candidate Vaccine Antigens

Although Plasmodium vivax is a leading cause of malaria around the world, only a handful of vivax antigens are being studied for vaccine development. Here, we investigated genetic signatures of selection and geospatial genetic diversity of two leading vivax vaccine antigens – Plasmodium vivax merozoite surface protein 1 (pvmsp-1) and Plasmodium vivax circumsporozoite protein (pvcsp). Using scalable next-generation sequencing, we deep-sequenced amplicons of the 42 kDa region of pvmsp-1 (n = 44) and the complete gene of pvcsp (n = 47) from Cambodian isolates. These sequences were then compared with global parasite populations obtained from GenBank. Using a combination of statistical and phylogenetic methods to assess for selection and population structure, we found strong evidence of balancing selection in the 42 kDa region of pvmsp-1, which varied significantly over the length of the gene, consistent with immune-mediated selection. In pvcsp, the highly variable central repeat region also showed patterns consistent with immune selection, which were lacking outside the repeat. The patterns of selection seen in both genes differed from their P. falciparum orthologs. In addition, we found that, similar to merozoite antigens from P. falciparum malaria, genetic diversity of pvmsp-1 sequences showed no geographic clustering, while the non-merozoite antigen, pvcsp, showed strong geographic clustering. These findings suggest that while immune selection may act on both vivax vaccine candidate antigens, the geographic distribution of genetic variability differs greatly between these two genes. The selective forces driving this diversification could lead to antigen escape and vaccine failure. Better understanding the geographic distribution of genetic variability in vaccine candidate antigens will be key to designing and implementing efficacious vaccines.     

Balancing Selection on a Regulatory Region Exhibiting Ancient Variation That Predates Human–Neandertal Divergence

Ancient population structure shaping contemporary genetic variation has been recently appreciated and has important implications regarding our understanding of the structure of modern human genomes. We identified a ∼36-kb DNA segment in the human genome that displays an ancient substructure. The variation at this locus exists primarily as two highly divergent haplogroups. One of these haplogroups (the NE1 haplogroup) aligns with the Neandertal haplotype and contains a 4.6-kb deletion polymorphism in perfect linkage disequilibrium with 12 single nucleotide polymorphisms (SNPs) across diverse populations. The other haplogroup, which does not contain the 4.6-kb deletion, aligns with the chimpanzee haplotype and is likely ancestral. Africans have higher overall pairwise differences with the Neandertal haplotype than Eurasians do for this NE1 locus (p<10⁻¹⁵). Moreover, the nucleotide diversity at this locus is higher in Eurasians than in Africans. These results mimic signatures of recent Neandertal admixture contributing to this locus. However, an in-depth assessment of the variation in this region across multiple populations reveals that African NE1 haplotypes, albeit rare, harbor more sequence variation than NE1 haplotypes found in Europeans, indicating an ancient African origin of this haplogroup and refuting recent Neandertal admixture. Population genetic analyses of the SNPs within each of these haplogroups, along with genome-wide comparisons revealed significant F_ST (p = 0.00003) and positive Tajima''s D (p = 0.00285) statistics, pointing to non-neutral evolution of this locus. The NE1 locus harbors no protein-coding genes, but contains transcribed sequences as well as sequences with putative regulatory function based on bioinformatic predictions and in vitro experiments. We postulate that the variation observed at this locus predates Human–Neandertal divergence and is evolving under balancing selection, especially among European populations.     

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.     

Genomic scanning using AFLP to detect loci under selection in the moss Funaria hygrometrica along a climate gradient in the Sierra Nevada Mountains,Spain

The common cord moss Funaria hygrometrica has a worldwide distribution and thrives in a wide variety of environments. Here, we studied the genetic diversity in F. hygrometrica along an abiotic gradient in the Mediterranean high mountain of Sierra Nevada (Spain) using a genome scan method. Eighty‐four samples from 17 locations from 24 to 2700 m were fingerprinted based on their amplified fragment length polymorphism (AFLP) banding pattern. Using PCA and Bayesian inference we found that the genetic diversity was structured in three or four clusters, respectively. Using a genome scan method we identified 13 outlier loci, which showed a signature of positive selection. Partial Mantel tests were performed between the Euclidean distance matrices of geographic and climatic variables, versus the pair‐wise genetic distance of the AFLP dataset and AFLP‐positive outliers dataset. AFLP‐positive outlier data were significantly correlated with the gradient of the climatic variables, suggesting adaptive variation among populations of F. hygrometrica along the Sierra Nevada Mountains. We highlight the additional analyses necessary to identify the nature of these loci, and their biological role in the adaptation process.     

Contrasting patterns of diversity and population differentiation at the innate immunity gene toll-like receptor 2 (TLR2) in two sympatric rodent species

Comparing patterns of diversity and divergence between populations at immune genes and neutral markers can give insights into the nature and geographic scale of parasite-mediated selection. To date, studies investigating such patterns of selection in vertebrates have primarily focused on the acquired branch of the immune system, whereas it remains largely unknown how parasite-mediated selection shapes innate immune genes both within and across vertebrate populations. Here, we present a study on the diversity and population differentiation at the innate immune gene Toll-like receptor 2 (TLR2) across nine populations of yellow-necked mice (Apodemus flavicollis) and bank voles (Myodes glareolus) in southern Sweden. In yellow-necked mice, TLR2 diversity was very low, as was TLR2 population differentiation compared to neutral loci. In contrast, several TLR2 haplotypes co-occurred at intermediate frequencies within and across bank vole populations, and pronounced isolation by distance between populations was observed. The diversity and differentiation at neutral loci was similar in the two species. These results indicate that parasite-mediated selection has been acting in dramatically different ways on a given immune gene in ecologically similar and sympatric species. Furthermore, the finding of TLR2 population differentiation at a small geographical scale in bank voles highlights that vertebrate innate immune defense may be evolutionarily more dynamic than has previously been appreciated.     

Moving beyond heritability in the search for coral adaptive potential

Global environmental change is happening at unprecedented rates. Coral reefs are among the ecosystems most threatened by global change. For wild populations to persist, they must adapt. Knowledge shortfalls about corals' complex ecological and evolutionary dynamics, however, stymie predictions about potential adaptation to future conditions. Here, we review adaptation through the lens of quantitative genetics. We argue that coral adaptation studies can benefit greatly from “wild” quantitative genetic methods, where traits are studied in wild populations undergoing natural selection, genomic relationship matrices can replace breeding experiments, and analyses can be extended to examine genetic constraints among traits. In addition, individuals with advantageous genotypes for anticipated future conditions can be identified. Finally, genomic genotyping supports simultaneous consideration of how genetic diversity is arrayed across geographic and environmental distances, providing greater context for predictions of phenotypic evolution at a metapopulation scale.