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.     

Spatially and temporally fluctuating selection at non-MHC immune genes: evidence from TAP polymorphism in populations of brown trout (Salmo trutta, L.)

Temporal samples of Danish brown trout (Salmo trutta) from populations representing varying geographical scales were analysed using eight putatively neutral microsatellite loci and two microsatellite loci embedded in TAP genes (Transporter associated with Antigen Processing). These genes encode molecules that are central to the major histocompatibility complex (MHC) class I restricted antigen presentation and thus integral components in the adaptive immune system. As such, they could be influenced by selection, driven by pathogens and parasites in a manner similar to MHC genes. Analysis of allele frequencies at presumably neutral microsatellite loci revealed a temporally unstable population structure within regions, while the population structure was stable over time among regions. Analyses of the two TAP markers indicated an effect of selection at both a regional and micro-geographical spatial scale. Moreover, signals of divergent selection among temporal samples within localities suggest that selection also might fluctuate at a temporal scale. These results suggest that immune genes other than the classical MHC class I and II might be subject to selection and warrant further studies of functional polymorphism of such genes in natural populations.     

Is MHC diversity a better marker for conservation than neutral genetic diversity? A case study of two contrasting dolphin populations

Genetic diversity is essential for populations to adapt to changing environments. Measures of genetic diversity are often based on selectively neutral markers, such as microsatellites. Genetic diversity to guide conservation management, however, is better reflected by adaptive markers, including genes of the major histocompatibility complex (MHC). Our aim was to assess MHC and neutral genetic diversity in two contrasting bottlenose dolphin (Tursiops aduncus) populations in Western Australia—one apparently viable population with high reproductive output (Shark Bay) and one with lower reproductive output that was forecast to decline (Bunbury). We assessed genetic variation in the two populations by sequencing the MHC class II DQB, which encompasses the functionally important peptide binding regions (PBR). Neutral genetic diversity was assessed by genotyping twenty‐three microsatellite loci. We confirmed that MHC is an adaptive marker in both populations. Overall, the Shark Bay population exhibited greater MHC diversity than the Bunbury population—for example, it displayed greater MHC nucleotide diversity. In contrast, the difference in microsatellite diversity between the two populations was comparatively low. Our findings are consistent with the hypothesis that viable populations typically display greater genetic diversity than less viable populations. The results also suggest that MHC variation is more closely associated with population viability than neutral genetic variation. Although the inferences from our findings are limited, because we only compared two populations, our results add to a growing number of studies that highlight the usefulness of MHC as a potentially suitable genetic marker for animal conservation. The Shark Bay population, which carries greater adaptive genetic diversity than the Bunbury population, is thus likely more robust to natural or human‐induced changes to the coastal ecosystem it inhabits.     

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.     

Diversity at the major histocompatibility complex Class II in the platypus, Ornithorhynchus anatinus

The platypus (Ornithorhynchus anatinus) is the sole survivor of a previously widely distributed and diverse lineage of ornithorhynchid monotremes. Its dependence on healthy water systems imposes an inherent sensitivity to habitat degradation and climate change. Here, we compare genetic diversity at the major histocompatibility complex (MHC) Class II-DZB gene and 3 MHC-associated microsatellite markers with diversity at 6 neutral microsatellite markers in 70 platypuses from across their range, including the mainland of Australia and the isolated populations of Tasmania, King Island, and Kangaroo Island. Overall, high DZB diversity was observed in the platypus, with 57 DZB β1 alleles characterized. Significant positive selection was detected within the DZB peptide-binding region, promoting variation in this domain. Low levels of genetic diversity were detected at all markers in the 2 island populations, King Island (endemic) and Kangaroo Island (introduced), with the King Island platypuses monomorphic at the DZB locus. Loss of MHC diversity on King Island is of concern, as the population may have compromised immunological fitness and reduced ability to resist changing environmental conditions.     

Variation in MHC genotypes in two populations of house sparrow (Passer domesticus) with different population histories

Small populations are likely to have a low genetic ability for disease resistance due to loss of genetic variation through inbreeding and genetic drift. In vertebrates, the highest genetic diversity of the immune system is located at genes within the major histocompatibility complex (MHC). Interestingly, parasite‐mediated selection is thought to potentially maintain variation at MHC loci even in populations that are monomorphic at other loci. Therefore, general loss of genetic variation in the genome may not necessarily be associated with low variation at MHC loci. We evaluated inter‐ and intrapopulation variation in MHC genotypes between an inbred (Aldra) and a relatively outbred population (Hestmannøy) of house sparrows (Passer domesticus) in a metapopulation at Helgeland, Norway. Genomic (gDNA) and transcribed (cDNA) alleles of functional MHC class I and IIB loci, along with neutral noncoding microsatellite markers, were analyzed to obtain relevant estimates of genetic variation. We found lower allelic richness in microsatellites in the inbred population, but high genetic variation in MHC class I and IIB loci in both populations. This suggests that also the inbred population could be under balancing selection to maintain genetic variation for pathogen resistance.     

Adaptive and Neutral Genetic Variation and Colonization History of Atlantic Salmon, Salmo salar

Synopsis I combined neutral microsatellite markers with the major histocompatibility complex (MHC) class IIB to study genetic differentiation and colonization history in Atlantic salmon, Salmo salar, in the Baltic Sea and in the north-eastern Atlantic. Baltic salmon populations have lower levels of microsatellite genetic variation, in terms of heterozygosity and allelic richness than Atlantic populations, confirming earlier findings with other genetic markers, suggesting that the Baltic Sea populations have been exposed to genetic bottlenecks, most likely at a founding event. On the other hand, the level of MHC variation was similar in the Baltic and in the north-eastern Atlantic, indicating that positive balancing selection has increased the level of MHC-variation. Both microsatellite and MHC class IIB genetic variation give strong support to the hypothesis that the Baltic salmon are of a biphyletic origin, the southern population in this study is strongly differentiated from both the northern Baltic salmon populations and from the north-eastern Atlantic populations. Salmon may have colonized the northern Baltic Sea either from the south, via the so called “N?rke strait” or from the north, via a proposed historical connection between the White Sea and the northern Baltic. At microsatellites, no significant isolation-by distance was found at either colonization route. At the MHC, populations were significantly isolated by distance when assuming that colonization occurred via the “N?rke strait”.     

Genetic Variation of Major Histocompatibility Complex and Microsatellite Loci: A Comparison in Bighorn Sheep

Examining and comparing genetic variation for major histocompatibility complex (MHC) and microsatellite (MS) loci in the same individuals provides an opportunity to understand the forces influencing genetic variation. We examined five MHC and three MS loci in 235 bighorn sheep (Ovis canadensis) from 14 populations and found that both types of loci were highly variable and were in Hardy-Weinberg proportions. Mean F(ST) values for both markers were very similar and MHC and MS genetic variability was predominantly distributed within rather than among populations. However, analyses of genetic distances and tree topologies revealed different spatial patterns of variation for the two types of loci. Collectively, these results indicated that neutral forces substantially influenced MS and MHC variation, and they provided limited evidence for selection acting on the MHC.     

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.     

Detecting the signature of selection on immune genes in highly structured populations of wild sheep (Ovis dalli)

The confounding effects of population structure complicate efforts to identify regions of the genome under the influence of selection in natural populations. Here we test for evidence of selection in three genes involved in vertebrate immune function - the major histocompatibility complex (MHC), interferon gamma (IFNG) and natural resistance associated macrophage polymorphism (NRAMP) - in highly structured populations of wild thinhorn sheep (Ovis dalli). We examined patterns of variation at microsatellite loci linked to these gene regions and at the DNA sequence level. Simple Watterson's tests indicated balancing selection at all three gene regions. However, evidence for selection was confounded by population structure, as the Watterson's test statistics from linked markers were not outside of the range of values from unlinked and presumably neutral microsatellites. The translated coding sequences of thinhorn IFNG and NRAMP are fixed and identical to those of domestic sheep (Ovis aries). In contrast, the thinhorn MHC DRB locus shows significant evidence of overdominance through both an excess of nonsynonymous substitution and trans-species polymorphism. The failure to detect balancing selection at microsatellite loci linked to the MHC is likely the result of recombination between the markers and expressed gene regions.     

Less can be more: loss of MHC functional diversity can reflect adaptation to novel conditions during fish invasions

The ability of invasive species to adapt to novel conditions depends on population size and environmental mismatch, but also on genetic variation. Away from their native range, invasive species confronted with novel selective pressures may display different levels of neutral versus functional genetic variation. However, the majority of invasion studies have only examined genetic variation at neutral markers, which may reveal little about how invaders adapt to novel environments. Salmonids are good model systems to examine adaptation to novel pressures because they have been translocated all over the world and represent major threats to freshwater biodiversity in the Southern Hemisphere, where they have become invasive. We examined patterns of genetic differentiation at seven putatively neutral (microsatellites) loci and one immune‐related major histocompatibility complex (MHC class II‐β) locus among introduced rainbow trout living in captivity (farmed) or under natural conditions (naturalized) in Chilean Patagonia. A significant positive association was found between differentiation at neutral and functional markers, highlighting the role of neutral evolutionary forces in shaping genetic variation at immune‐related genes in salmonids. However, functional (MHC) genetic diversity (but not microsatellite diversity) decreased with time spent in the wild since introduction, suggesting that there was selection against alleles associated with captive rearing of donor populations that do not provide an advantage in the wild. Thus, although high genetic diversity may initially enhance fitness in translocated populations, it does not necessarily reflect invasion success, as adaptation to novel conditions may result in rapid loss of functional MHC diversity.     

TEMPORAL VARIATION AT THE MHC CLASS IIB IN WILD POPULATIONS OF THE GUPPY (POECILIA RETICULATA)

Understanding genetic diversity in natural populations is a fundamental objective of evolutionary biology. The immune genes of the major histocompatibility complex (MHC) are excellent candidates to study such diversity because they are highly polymorphic in populations. Although balancing selection may be responsible for maintaining diversity at these functionally important loci, temporal variation in selection pressure has rarely been examined. We examine temporal variation in MHC class IIB diversity in nine guppy (Poecilia reticulata) populations over two years. We found that five of the populations changed significantly more at the MHC than at neutral (microsatellite) loci as measured by F_ST, which suggests that the change at the MHC was due to selection and not neutral processes. Additionally, pairwise population differentiation measures at the MHC were higher in 2007 than in 2006, with the signature of selection changing from homogenizing to diversifying selection or neutral evolution. Interestingly, within the populations the magnitude of the change at the MHC between years was related to the change in the proportion of individuals infected by a common parasite, indicating a link between genetic structure and the parasite. Our data thereby implicate temporal variation in selective pressure as an important mechanism maintaining diversity at the MHC in wild populations.     

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

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

Is genetic variability so important? Non-native salmonids in South America

Three salmonid species introduced in Patagonian national parks in Argentine have experienced different degrees of expansion. Atlantic salmon Salmo salar is restricted to a few river-lake systems and its populations have been declining over recent years. Both rainbow Oncorhynchus mykiss and brown trout Salmo trutta populations have expanded from their introduction sites and now occupy a wide range of freshwater ecosystems. Genetic variation at the same neutral markers (microsatellite loci) was examined for different populations of the three species acclimatized to the same areas, and compared with that of native populations. Founder effects denoted as reduced variability and great differentiation with respect to the native populations were detected. Significant reduction in variability has not been an obstacle for successful adaptation of rainbow and brown trout, indicating that genetic variability per se cannot be claimed as the reason for their different outcomes in the new habitats.     

Extensive polymorphism and geographical variation at a positively selected MHC class II B gene of the lesser kestrel (Falco naumanni)

Understanding the selective forces that shape genetic variation in natural populations remains a high priority in evolutionary biology. Genes at the major histocompatibility complex (MHC) have become excellent models for the investigation of adaptive variation and natural selection because of their crucial role in fighting off pathogens. Here we present one of the first data sets examining patterns of MHC variation in wild populations of a bird of prey, the lesser kestrel, Falco naumanni . We report extensive polymorphism at the second exon of a putatively functional MHC class II gene, Fana- DAB*1. Overall, 103 alleles were isolated from 121 individuals sampled from Spain to Kazakhstan. Bayesian inference of diversifying selection suggests that several amino acid sites may have experienced strong positive selection (ω = 4.02 per codon). The analysis also suggests a prominent role of recombination in generating and maintaining MHC diversity (ρ = 4 Nc  = 0.389 per codon, θ = 0.017 per codon). Both the Fana -DAB*1 locus and a set of eight polymorphic microsatellite markers revealed an isolation-by-distance pattern across the Western Palaearctic ( r  = 0.67; P  = 0.01 and r  = 0.50; P  = 0.04, respectively). Nonetheless, geographical variation at the MHC contrasts with relatively uniform distributions in the frequencies of microsatellite alleles. In addition, we found lower fixation rates in the MHC than those predicted by genetic drift after controlling for neutral mitochondrial sequences. Our results therefore underscore the role of balancing selection as well as spatial variations in parasite-mediated selection regimes in shaping MHC diversity when gene flow is limited.     

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.     

Between‐year variation of MHC allele frequencies in great reed warblers: selection or drift?

The major histocompatibility complex (MHC) genes are extremely polymorphic and this variation is assumed to be maintained by balancing selection. Cyclic interactions between pathogens and their hosts could generate such selection, and specific MHC alleles or heterozygosity at certain MHC loci have been shown to confer resistance against particular pathogens. Here we compare the temporal variation in allele frequencies of 23 MHC class I alleles with that of 23 neutral microsatellite markers in adult great reed warblers (a passerine bird) in nine successive cohorts. Overall, the MHC alleles showed a significantly higher variation in allele frequencies between cohorts than the microsatellite alleles, using a multi-variate genetic analysis (amova). The frequency of two specific MHC alleles, A3e (P = 0.046) and B4b (P = 0.0018), varied more between cohorts than expected from random, whereas none of the microsatellite alleles showed fluctuations exceeding the expectation from stochastic variation. These results imply that the variation in MHC allele frequencies between cohorts is not a result of demographic events, but rather an effect of selection favouring different MHC alleles in different years.     

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.     

Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes

Microsatellites and mitochondrial DNA (mtDNA) have traditionally been used in population genetics because of their variability and presumed neutrality, whereas genes of the major histocompatibility complex (MHC) are increasingly of interest because strong selective pressures shape their standing variation. Despite the potential for MHC genes, microsatellites, and mtDNA sequences to complement one another in deciphering population history and demography, the three are rarely used in tandem. Here we report on MHC, microsatellite, and mtDNA variability in a single large population of the eastern tiger salamander (Ambystoma tigrinum tigrinum). We use the mtDNA mismatch distribution and, on microsatellite data, the imbalance index and bottleneck tests to infer aspects of population history and demography. Haplotype and allelic variation was high at all loci surveyed, and heterozygosity was high at the nuclear loci. We find concordance among neutral molecular markers that suggests our study population originated from post-Pleistocene expansions of multiple, fragmented sources that shared few migrants. Differences in N(e) estimates derived from haploid and diploid genetic markers are potentially attributable to secondary contact among source populations that experienced rapid mtDNA divergence and comparatively low levels of nuclear DNA divergence. We find strong evidence of natural selection acting on MHC genes and estimate long-term effective population sizes (N(e)) that are very large, making small selection intensities significant evolutionary forces in this population.     

Genetic diversity at neutral and adaptive loci determines individual fitness in a long-lived territorial bird

There is compelling evidence about the manifest effects of inbreeding depression on individual fitness and populations' risk of extinction. The majority of studies addressing inbreeding depression on wild populations are generally based on indirect measures of inbreeding using neutral markers. However, the study of functional loci, such as genes of the major histocompatibility complex (MHC), is highly recommended. MHC genes constitute an essential component of the immune system of individuals, which is directly related to individual fitness and survival. In this study, we analyse heterozygosity fitness correlations of neutral and adaptive genetic variation (22 microsatellite loci and two loci of the MHC class II, respectively) with the age of recruitment and breeding success of a decimated and geographically isolated population of a long-lived territorial vulture. Our results indicate a negative correlation between neutral genetic diversity and age of recruitment, suggesting that inbreeding may be delaying reproduction. We also found a positive correlation between functional (MHC) genetic diversity and breeding success, together with a specific positive effect of the most frequent pair of cosegregating MHC alleles in the population. Globally, our findings demonstrate that genetic depauperation in small populations has a negative impact on the individual fitness, thus increasing the populations' extinction risk.