Genetic variation and the mitogen-activated protein kinase (MAPK) signaling pathway

Non-synonymous single nucleotide polymorphisms (nsSNPs) are known to alter protein function, contributing to disease susceptibility. This report explores the nature of nsSNPs in the gene products of the highly conserved mitogen-activated protein kinase (MAPK) signaling pathways already implicated in cancer development. MAPK signaling pathways regulate cellular processes such as proliferation, differentiation, apoptosis, and survival mediated through interconnected signaling cascades. Using the dbSNP database, we have identified 25 nsSNPs in 17 out of 98 MAPK genes studied. Computational algorithms were used to predict whether the amino acid substitutions were evolutionarily tolerated, or affected putative functional units such as phosphorylation sites, protein motifs and domains. This study predicts that 36% of nsSNPs are likely to have functional consequences, based on evolutionary conservation analysis, and 36% based on phosphorylation prediction analysis. All such nsSNPs represent potentially functional and disease-causing/modifying alleles. More interestingly, the epistatic relationships discussed in this report represent potential synergistic/ antagonistic/additive effects of nsSNP combinations found within the same protein, or within members of the same protein complex and cascades. This strategy can effectively determine which nsSNPs potentially alter protein function, and can be utilized to study the genetic architecture and disease association of other biological protein complexes and networks.     

Association of the I264T Variant in the Sulfide Quinone Reductase-Like (SQRDL) Gene with Osteoporosis in Korean Postmenopausal Women

To identify novel susceptibility variants for osteoporosis in Korean postmenopausal women, we performed a genome-wide association analysis of 1180 nonsynonymous single nucleotide polymorphisms (nsSNPs) in 405 individuals with osteoporosis and 722 normal controls of the Korean Association Resource cohort. A logistic regression analysis revealed 72 nsSNPs that showed a significant association with osteoporosis (p<0.05). The top 10 nsSNPs showing the lowest p-values (p = 5.2×10^-4–8.5×10^-3) were further studied to investigate their effects at the protein level. Based on the results of an in silico prediction of the protein’s functional effect based on amino acid alterations and a sequence conservation evaluation of the amino acid residues at the positions of the nsSNPs among orthologues, we selected one nsSNP in the SQRDL gene (rs1044032, SQRDL I264T) as a meaningful genetic variant associated with postmenopausal osteoporosis. To assess whether the SQRDL I264T variant played a functional role in the pathogenesis of osteoporosis, we examined the in vitro effect of the nsSNP on bone remodeling. Overexpression of the SQRDL I264T variant in the preosteoblast MC3T3-E1 cells significantly increased alkaline phosphatase activity, mineralization, and the mRNA expression of osteoblastogenesis markers, Runx2, Sp7, and Bglap genes, whereas the SQRDL wild type had no effect or a negative effect on osteoblast differentiation. Overexpression of the SQRDL I264T variant did not affect osteoclast differentiation of the primary-cultured monocytes. The known effects of hydrogen sulfide (H₂S) on bone remodeling may explain the findings of the current study, which demonstrated the functional role of the H₂S-catalyzing enzyme SQRDL I264T variant in osteoblast differentiation. In conclusion, the results of the statistical and experimental analyses indicate that the SQRDL I264T nsSNP may be a significant susceptibility variant for osteoporosis in Korean postmenopausal women that is involved in osteoblast differentiation.     

Predicting deleterious nsSNPs: an analysis of sequence and structural attributes

Background  

There has been an explosion in the number of single nucleotide polymorphisms (SNPs) within public databases. In this study we focused on non-synonymous protein coding single nucleotide polymorphisms (nsSNPs), some associated with disease and others which are thought to be neutral. We describe the distribution of both types of nsSNPs using structural and sequence based features and assess the relative value of these attributes as predictors of function using machine learning methods. We also address the common problem of balance within machine learning methods and show the effect of imbalance on nsSNP function prediction. We show that nsSNP function prediction can be significantly improved by 100% undersampling of the majority class. The learnt rules were then applied to make predictions of function on all nsSNPs within Ensembl.     

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.     

Can balancing selection on MHC loci counteract genetic drift in small fragmented populations of black grouse?

The ability of natural populations to adapt to new environmental conditions is crucial for their survival and partly determined by the standing genetic variation in each population. Populations with higher genetic diversity are more likely to contain individuals that are better adapted to new circumstances than populations with lower genetic diversity. Here, we use both neutral and major histocompatibility complex (MHC) markers to test whether small and highly fragmented populations hold lower genetic diversity than large ones. We use black grouse as it is distributed across Europe and found in populations with varying degrees of isolation and size. We sampled 11 different populations; five continuous, three isolated, and three small and isolated. We tested patterns of genetic variation in these populations using three different types of genetic markers: nine microsatellites and 21 single nucleotide polymorphisms (SNPs) which both were found to be neutral, and two functional MHC genes that are presumably under selection. The small isolated populations displayed significantly lower neutral genetic diversity compared to continuous populations. A similar trend, but not as pronounced, was found for genotypes at MHC class II loci. Populations were less divergent at MHC genes compared to neutral markers. Measures of genetic diversity and population genetic structure were positively correlated among microsatellites and SNPs, but none of them were correlated to MHC when comparing all populations. Our results suggest that balancing selection at MHC loci does not counteract the power of genetic drift when populations get small and fragmented.     

Computational refinement of functional single nucleotide polymorphisms associated with ATM gene

Background

Understanding and predicting molecular basis of disease is one of the major challenges in modern biology and medicine. SNPs associated with complex disorders can create, destroy, or modify protein coding sites. Single amino acid substitutions in the ATM gene are the most common forms of genetic variations that account for various forms of cancer. However, the extent to which SNPs interferes with the gene regulation and affects cancer susceptibility remains largely unknown.

Principal findings

We analyzed the deleterious nsSNPs associated with ATM gene based on different computational methods. An integrative scoring system and sequence conservation of amino acid residues was adapted for a priori nsSNP analysis of variants associated with cancer. We further extended our approach on SNPs that could potentially influence protein Post Translational Modifications in ATM gene.

Significance

In the lack of adequate prior reports on the possible deleterious effects of nsSNPs, we have systematically analyzed and characterized the functional variants in both coding and non coding region that can alter the expression and function of ATM gene. In silico characterization of nsSNPs affecting ATM gene function can aid in better understanding of genetic differences in disease susceptibility.     

Inbreeding depression and low between-population heterosis in recently diverged experimental populations of a selfing species

In fragmented populations, genetic drift and selection reduce genetic diversity, which in turn results in a loss of fitness or in a loss of evolvability. Genetic rescue, that is, controlled input of diversity from distant populations, may restore evolutionary potential, whereas outbreeding depression might counteract the positive effect of this strategy. We carried out self-pollination and crosses within and between populations in an experimental subdivided population of a selfing species, Triticum aestivum L., to estimate the magnitude of these two phenomena. Surprisingly, for a self-fertilizing species, we found significant inbreeding depression within each population for four of the six traits studied, indicating that mildly deleterious mutations were still segregating in these populations. The progeny of within- and between-population crosses was very similar, indicating low between-population heterosis and little outbreeding depression. We conclude that relatively large population effective sizes prevented fixation of a high genetic load and that local adaptation was limited in these recently diverged populations. The kinship coefficient estimated between the parents using 20 neutral markers was a poor predictor of the progeny phenotypic values, indicating that there was a weak link between neutral diversity and genes controlling fitness-related traits. These results show that when assessing the viability of natural populations and the need for genetic rescue, the use of neutral markers should be complemented with information about the presence of local adaptation in the subdivided population.     

Knowledge-based computational mutagenesis for predicting the disease potential of human non-synonymous single nucleotide polymorphisms

Certain genetic variations in the human population are associated with heritable diseases, and single nucleotide polymorphisms (SNPs) represent the most common form of such differences in DNA sequence. In particular, substantial interest exists in determining whether a non-synonymous SNP (nsSNP), leading to a single residue replacement in the translated protein product, is neutral or disease-related. The nature of protein structure-function relationships suggests that nsSNP effects, either benign or leading to aberrant protein function possibly associated with disease, are dependent on relative structural changes introduced upon mutation. In this study, we characterize a representative sampling of 1790 documented neutral and disease-related human nsSNPs mapped to 243 diverse human protein structures, by quantifying environmental perturbations in the associated proteins with the use of a computational mutagenesis methodology that relies on a four-body, knowledge-based, statistical contact potential. These structural change data are used as attributes to generate a vector representation for each nsSNP, in combination with additional features reflecting sequence and structure of the corresponding protein. A trained model based on the random forest supervised classification algorithm achieves 76% cross-validation accuracy. Our classifier performs at least as well as other methods that use significantly larger datasets of nsSNPs for model training, and the novelty of our attributes differentiates the model as an orthogonal approach that can be utilized in conjunction with other techniques. A dedicated server for obtaining predictions, as well as supporting datasets and documentation, is available at http://proteins.gmu.edu/automute.     

Geographical distribution of selected and putatively neutral SNPs in Southeast Asian malaria parasites

Loci targeted by directional selection are expected to show elevated geographical population structure relative to neutral loci, and a flurry of recent papers have used this rationale to search for genome regions involved in adaptation. Studies of functional mutations that are known to be under selection are particularly useful for assessing the utility of this approach. Antimalarial drug treatment regimes vary considerably between countries in Southeast Asia selecting for local adaptation at parasite loci underlying resistance. We compared the population structure revealed by 10 nonsynonymous mutations (nonsynonymous single-nucleotide polymorphisms [nsSNPs]) in four loci that are known to be involved in antimalarial drug resistance, with patterns revealed by 10 synonymous mutations (synonymous single-nucleotide polymorphisms [sSNPs]) in housekeeping genes or genes of unknown function in 755 Plasmodium falciparum infections collected from 13 populations in six Southeast Asian countries. Allele frequencies at known nsSNPs underlying resistance varied markedly between locations (F(ST) = 0.18-0.66), with the highest frequencies on the Thailand-Burma border and the lowest frequencies in neighboring Lao PDR. In contrast, we found weak but significant geographic structure (F(ST) = 0-0.14) for 8 of 10 sSNPs. Importantly, all 10 nsSNPs showed significantly higher F(ST) (P < 8 x 10(-5)) than simulated neutral expectations based on observed F(ST) values in the putatively neutral sSNPs. This result was unaffected by the methods used to estimate allele frequencies or the number of populations used in the simulations. Given that dense single-nucleotide polymorphism (SNP) maps and rapid SNP assay methods are now available for P. falciparum, comparing genetic differentiation across the genome may provide a valuable aid to identifying parasite loci underlying local adaptation to drug treatment regimes or other selective forces. However, the high proportion of polymorphic sites that appear to be under balancing selection (or linked to selected sites) in the P. falciparum genome violates the central assumption that selected sites are rare, which complicates identification of outlier loci, and suggests that caution is needed when using this approach.     

Habitat fragmentation differentially shapes neutral and immune gene variation in a tropical bird species

Habitat fragmentation is a major cause of biodiversity loss, responsible for an alteration of intraspecific patterns of neutral genetic diversity and structure. Although neutral genetic variation can be informative for demographic inferences, it may be a poor predictor of adaptive genetic diversity and thus of the consequences of habitat fragmentation on selective evolutionary processes. In this context, we contrasted patterns of genetic diversity and structure of neutral loci (microsatellites) and immune genes (i.e., toll-like receptors) in an understorey bird species, the wedge-billed woodcreeper Glyphorynchus spirurus. The objectives were (1) to investigate forest fragmentation effects on population genetic diversity, (2) to disentangle the relative role of demography (genetic drift and migration) and selection, and (3) to assess whether immunogenetic patterns could be associated with variation of ectoparasite (i.e., ticks) pressures. Our results revealed an erosion of neutral genetic diversity and a substantial genetic differentiation among fragmented populations, resulting from a decrease in landscape connectivity and leading to the divergence of distinct genetic pools at a small spatial scale. Patterns of genetic diversity observed for TLR4 and TLR5 were concordant with neutral genetic patterns, whereas those observed for TLR3 and TLR21 were discordant. This result underlines that the dominant evolutionary force shaping immunogenetic diversity (genetic drift vs. selection) may be different depending on loci considered. Finally, tick prevalence was higher in fragmented environments. We discussed the hypothesis that pathogen selective pressures may contribute to maintain adaptive genetic diversity despite the negative demographic effect of habitat fragmentation on neutral genetic diversity.Subject terms: Tropical ecology, Genetic variation     

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake,the Eastern Massasauga (Sistrurus catenatus)

Genetic diversity is fundamental to maintaining the long‐term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.     

Computational SNP Analysis: Current Approaches and Future Prospects

The computational approaches in determining disease-associated Non-synonymous single nucleotide polymorphisms (nsSNPs) have evolved very rapidly. Large number of deleterious and disease-associated nsSNP detection tools have been developed in last decade showing high prediction reliability. Despite of all these highly efficient tools, we still lack the accuracy level in determining the genotype–phenotype association of predicted nsSNPs. Furthermore, there are enormous questions that are yet to be computationally compiled before we might talk about the prediction accuracy. Earlier we have incorporated molecular dynamics simulation approaches to foster the accuracy level of computational nsSNP analysis roadmap, which further helped us to determine the changes in the protein phenotype associated with the computationally predicted disease-associated mutation. Here we have discussed on the present scenario of computational nsSNP characterization technique and some of the questions that are crucial for the proper understanding of pathogenicity level for any disease associated mutations.     

Demographic and Genetic Patterns of Variation among Populations of Arabidopsis thaliana from Contrasting Native Environments

Background

Understanding the relationship between environment and genetics requires the integration of knowledge on the demographic behavior of natural populations. However, the demographic performance and genetic composition of Arabidopsis thaliana populations in the species'' native environments remain largely uncharacterized. This information, in combination with the advances on the study of gene function, will improve our understanding on the genetic mechanisms underlying adaptive evolution in A. thaliana.

Methodology/Principal Findings

We report the extent of environmental, demographic, and genetic variation among 10 A. thaliana populations from Mediterranean (coastal) and Pyrenean (montane) native environments in northeast Spain. Geographic, climatic, landscape, and soil data were compared. Demographic traits, including the dynamics of the soil seed bank and the attributes of aboveground individuals followed over a complete season, were also analyzed. Genetic data based on genome-wide SNP markers were used to describe genetic diversity, differentiation, and structure. Coastal and montane populations significantly differed in terms of environmental, demographic, and genetic characteristics. Montane populations, at higher altitude and farther from the sea, are exposed to colder winters and prolonged spring moisture compared to coastal populations. Montane populations showed stronger secondary seed dormancy, higher seedling/juvenile mortality in winter, and initiated flowering later than coastal populations. Montane and coastal regions were genetically differentiated, montane populations bearing lower genetic diversity than coastal ones. No significant isolation-by-distance pattern and no shared multilocus genotypes among populations were detected.

Conclusions/Significance

Between-region variation in climatic patterns can account for differences in demographic traits, such as secondary seed dormancy, plant mortality, and recruitment, between coastal and montane A. thaliana populations. In addition, differences in plant mortality can partly account for differences in the genetic composition of coastal and montane populations. This study shows how the interplay between variation in environmental, demographic, and genetic parameters may operate in natural A. thaliana populations.     

Standing genetic diversity and selection at functional gene loci are associated with differential invasion success in two non‐native fish species

Invasive species are expected to experience a unique combination of high genetic drift due to demographic factors while also experiencing strong selective pressures. The paradigm that reduced genetic diversity should limit the evolutionary potential of invasive species, and thus, their potential for range expansion has received little empirical support, possibly due to the choice of genetic markers. Our goal was to test for effects of genetic drift and selection at functional genetic markers as they relate to the invasion success of two paired invasive goby species, one widespread (successful) and one with limited range expansion (less successful). We genotyped fish using two marker types: single nucleotide polymorphisms (SNPs) in known‐function, protein‐coding genes and microsatellites to contrast the effects of neutral genetic processes. We identified reduced allelic variation in the invaded range for the less successful tubenose goby. SNPs putatively under selection were responsible for the observed differences in population structure between marker types for round goby (successful) but not tubenose goby (less successful). A higher proportion of functional loci experienced divergent selection for round goby, suggesting increased evolutionary potential in invaded ranges may be associated with round goby's greater invasion success. Genes involved in thermal tolerance were divergent for round goby populations but not tubenose goby, consistent with the hypothesis that invasion success for fish in temperate regions is influenced by capacity for thermal tolerance. Our results highlight the need to incorporate functional genetic markers in studies to better assess evolutionary potential for the improved conservation and management of species.     

Exploring the structural and functional effect of pRB by significant nsSNP in the coding region of <Emphasis Type="Italic">RB1</Emphasis> gene causing retinoblastoma

In this study, we identified the most deleterious nsSNP in RB1 gene through structural and functional properties of its protein (pRB) and investigated its binding affinity with E2F-2. Out of 956 SNPs, we investigated 12 nsSNPs in coding region in which three of them (SNPids rs3092895, rs3092903 and rs3092905) are commonly found to be damaged by I-Mutant 2.0, SIFT and PolyPhen programs. With this effort, we modeled the mutant pRB proteins based on these deleterious nsSNPs. From a comparison of total energy, stabilizing residues and RMSD of these three mutant proteins with native pRB protein, we identified that the major mutation is from Glutamic acid to Glycine at the residue position of 746 of pRB. Further, we compared the binding efficiency of both native and mutant pRB (E746G) with E2F-2. We found that mutant pRB has less binding affinity with E2F-2 as compared to native type. This is due to sixteen hydrogen bonding and two salt bridges that exist between native type and E2F-2, whereas mutant type makes only thirteen hydrogen bonds and one salt bridge with E2F-2. Based on our investigation, we propose that the SNP with an id rs3092905 could be the most deleterious nsSNP in RB1 gene causing retinoblastoma.     

Urban colonization through multiple genetic lenses: The city‐fox phenomenon revisited

Urbanization is driving environmental change on a global scale, creating novel environments for wildlife to colonize. Through a combination of stochastic and selective processes, urbanization is also driving evolutionary change. For instance, difficulty in traversing human‐modified landscapes may isolate newly established populations from rural sources, while novel selective pressures, such as altered disease risk, toxicant exposure, and light pollution, may further diverge populations through local adaptation. Assessing the evolutionary consequences of urban colonization and the processes underlying them is a principle aim of urban evolutionary ecology. In the present study, we revisited the genetic effects of urbanization on red foxes (Vulpes vulpes) that colonized Zurich, Switzerland. Through use of genome‐wide single nucleotide polymorphisms and microsatellite markers linked to the major histocompatibility complex (MHC), we expanded upon a previous neutral microsatellite study to assess population structure, characterize patterns of genetic diversity, and detect outliers associated with urbanization. Our results indicated the presence of one large evolutionary cluster, with substructure evident between geographic sampling areas. In urban foxes, we observed patterns of neutral and functional diversity consistent with founder events and reported increased differentiation between populations separated by natural and anthropogenic barriers. We additionally reported evidence of selection acting on MHC‐linked markers and identified outlier loci with putative gene functions related to energy metabolism, behavior, and immunity. We concluded that demographic processes primarily drove patterns of diversity, with outlier tests providing preliminary evidence of possible urban adaptation. This study contributes to our overall understanding of urban colonization ecology and emphasizes the value of combining datasets when examining evolutionary change in an increasingly urban world.     

Spatial and temporal patterns of neutral and adaptive genetic variation in the endangered African wild dog (Lycaon pictus)

Deciphering patterns of genetic variation within a species is essential for understanding population structure, local adaptation and differences in diversity between populations. Whilst neutrally evolving genetic markers can be used to elucidate demographic processes and genetic structure, they are not subject to selection and therefore are not informative about patterns of adaptive variation. As such, assessments of pertinent adaptive loci, such as the immunity genes of the major histocompatibility complex (MHC), are increasingly being incorporated into genetic studies. In this study, we combined neutral (microsatellite, mtDNA) and adaptive (MHC class II DLA‐DRB1 locus) markers to elucidate the factors influencing patterns of genetic variation in the African wild dog (Lycaon pictus); an endangered canid that has suffered extensive declines in distribution and abundance. Our genetic analyses found all extant wild dog populations to be relatively small (N_e < 30). Furthermore, through coalescent modelling, we detected a genetic signature of a recent and substantial demographic decline, which correlates with human expansion, but contrasts with findings in some other African mammals. We found strong structuring of wild dog populations, indicating the negative influence of extensive habitat fragmentation and loss of gene flow between habitat patches. Across populations, we found that the spatial and temporal structure of microsatellite diversity and MHC diversity were correlated and strongly influenced by demographic stability and population size, indicating the effects of genetic drift in these small populations. Despite this correlation, we detected signatures of selection at the MHC, implying that selection has not been completely overwhelmed by genetic drift.     

The Pattern of Neutral Molecular Variation under the Background Selection Model

Stochastic simulations of the infinite sites model were used to study the behavior of genetic diversity at a neutral locus in a genomic region without recombination, but subject to selection against deleterious alleles maintained by recurrent mutation (background selection). In large populations, the effect of background selection on the number of segregating sites approaches the effct on nucleotide site diversity, i.e., the reduction in genetic variability caused by background selection resembles that caused by a simple reduction in effective population size. We examined, by coalescence-based methods, the power of several tests for the departure from neutral expectation of the frequency spectra of alleles in samples from randomly mating populations (TAJIMA's, FU and LI's, and WATTERSON's tests). All of the tests have low power unless the selection against mutant alleles is extremely weak. In Drosophila, significant TAJIMA's tests are usually not obtained with empirical data sets from loci in genomic regions with restricted recombination frequencies and that exhibit low genetic diversity. This is consistent with the operation of background selection as opposed to selective sweeps. It remains to be decided whether background selection is sufficient to explain the observed extent of reduction in diversity in regions of restricted recombination.     

Selection maintains MHC diversity through a natural population bottleneck

A perceived consequence of a population bottleneck is the erosion of genetic diversity and concomitant reduction in individual fitness and evolutionary potential. Although reduced genetic variation associated with demographic perturbation has been amply demonstrated for neutral molecular markers, the effective management of genetic resources in natural populations is hindered by a lack of understanding of how adaptive genetic variation will respond to population fluctuations, given these are affected by selection as well as drift. Here, we demonstrate that selection counters drift to maintain polymorphism at a major histocompatibility complex (MHC) locus through a population bottleneck in an inbred island population of water voles. Before and after the bottleneck, MHC allele frequencies were close to balancing selection equilibrium but became skewed by drift when the population size was critically low. MHC heterozygosity generally conformed to Hardy-Weinberg expectations except in one generation during the population recovery where there was a significant excess of heterozygous genotypes, which simulations ascribed to strong differential MHC-dependent survival. Low allelic diversity and highly skewed frequency distributions at microsatellite loci indicated potent genetic drift due to a strong founder affect and/or previous population bottlenecks. This study is a real-time examination of the predictions of fundamental evolutionary theory in low genetic diversity situations. The findings highlight that conservation efforts to maintain the genetic health and evolutionary potential of natural populations should consider the genetic basis for fitness-related traits, and how such adaptive genetic diversity will vary in response to both the demographic fluctuations and the effects of selection.