Population genetic consequences of geographic disjunction: a prairie plant isolated on Great Lakes alvars

Species may often exhibit geographic variation in population genetic structure due to contemporary and historical variation in population size and gene flow. Here, we test the predictions that populations on the margins of a species' distribution contain less genetic variation and are more differentiated than populations towards the core of the range by comparing patterns of genetic variation at five microsatellite loci between disjunct and core populations of the perennial, allohexaploid herb Geum triflorum. We sampled nine populations isolated on alvar habitat within the eastern Great Lakes region in North America, habitats that include disjunct populations of several plant species, and compared these to 16 populations sampled from prairie habitat throughout the core of the species' distribution in midwestern Canada and the USA. Alvar populations exhibited much lower within-population diversity and contained only a subset of alleles found in prairie populations. We detected isolation by distance across the species' range and within alvar and prairie regions separately. As predicted, genetic differentiation was higher among alvar populations than among prairie populations, even after controlling for the geographic distance between sampled populations. Low diversity and high differentiation can be accounted for by the greater contemporary spatial isolation of alvar populations. However, the genetic structure of alvar populations may also have been influenced by postglacial range expansion and contraction. Our results are consistent with alvar populations being founded during an expansion of prairie habitat during the warmer, hypsithermal period approximately 5000 bp and subsequently becoming stranded on isolated alvar habitat as the climate grew cooler and wetter.     

Disjunct habitats as islands: genetic variability in the Caucasian rock lizard Lacerta portschinskii

Genetic diversity at 37 allozyme loci was surveyed in Lacerta portschinskii from contiguous populations and from a disjunct population. Indices of genetic diversity (heterozygosity, number of alleles per locus, and percentage of loci polymorphic) were greater in contiguous populations than in the smaller disjunct population. In this regard, disjunct populations appear to be similar to island populations. Indices of genetic diversity in Caucasian Lacerta are less than those reported from vagile lizard taxa and more similar to those of sit-and-wait predators. This revised version was published online in July 2006 with corrections to the Cover Date.     

The role of propagule pressure in the invasion success of bluegill sunfish, Lepomis macrochirus, in Japan

The bluegill sunfish, Lepomis macrochirus, is a widespread exotic species in Japan that is considered to have originated from 15 fish introduced from Guttenberg, Iowa, in 1960. Here, the genetic and phenotypic traits of Japanese populations were examined, together with 11 native populations of the USA using 10 microsatellite markers and six meristic traits. Phylogenetic analysis reconfirmed a single origin of Japanese populations, among which populations established in the 1960s were genetically close to Guttenberg population, keeping high genetic diversity comparable to the ancestral population. In contrast, genetic diversity of later-established populations significantly declined with genetic divergence from the ancestral population. Among the 1960s established populations, that from Lake Biwa showed a significant isolation-by-distance pattern with surrounding populations in which genetic bottlenecks increased with geographical distance from Lake Biwa. Although phenotypic divergence among populations was recognized in both neutral and adaptive traits, P(ST)-F(ST) comparisons showed that it is independent of neutral genetic divergence. Divergent selection was suggested in some populations from reservoirs with unstable habitats, while stabilizing selection was dominant. Accordingly, many Japanese populations of L. macrochirus appear to have derived from Lake Biwa population, expanding their distribution with population bottlenecks. Despite low propagule pressure, the invasion success of L. macrochirus is probably because of its drastic population growth in Lake Biwa shortly after its introduction, together with artificial transplantations. It not only enabled the avoidance of a loss in genetic diversity but also formed a major gene pool that supported local adaptation with high phenotypic plasticity.     

Genetic variation of the major histocompatibility complex (MHC class II β gene) in the threatened Gila trout, <Emphasis Type="Italic">Oncorhynchus gilae gilae</Emphasis>

Gila trout (Oncorhynchus gilae gilae) was federally protected in 1973 because of severe declines in abundance and geographic range size. At present, four relict genetic lineages of the species remain in mountain streams of New Mexico and Arizona, USA. Management actions aimed at species recovery, including hatchery production and restocking of formerly occupied streams, have been guided by information from non-functional genetic markers. In this study, we investigated genetic variation at exon 2 of the major histocompatibility complex (MHC) class II β gene that is involved in pathogen resistance and thus presumably under natural selection. Phylogenetic analysis revealed trans-species polymorphism and a significantly high ratio of non-synonymous to synonymous amino acid changes consistent with the action of historical balancing selection that maintained diversity at this locus in the past. However, Gila trout exhibited low allelic diversity (five alleles from 142 individuals assayed) compared to some other salmonid fishes, and populations that originated exclusively from hatcheries possessed three or fewer MHC alleles. Comparative analysis of genetic variation at MHC and six presumably neutrally evolving microsatellite loci revealed that genetic drift cannot be rejected as a primary force governing evolution of MHC in contemporary populations of Gila trout. Maintenance of diversity at MHC will require careful implementation of hatchery breeding protocols and continued protection of wild populations to prevent loss of allelic diversity due to drift.     

Patterns of selection and allele diversity of class I and class II major histocompatibility loci across the species range of sockeye salmon (Oncorhynchus nerka)

The major histocompatibility complex (MHC), an important component of the vertebrate immune system, provides an important suite of genes to examine the role of genetic diversity at non‐neutral loci for population persistence. We contrasted patterns of diversity at the two classical MHC loci in sockeye salmon (Oncorhynchus nerka), MHC class I (UBA) and MHC class II (DAB), and neutral microsatellite loci across 70 populations spanning the species range from Washington State to Japan. There was no correlation in allelic richness or heterozygosity between MHC loci or between MHC loci and microsatellites. The two unlinked MHC loci may be responding to different selective pressures; the distribution of F_ST values for the two loci was uncorrelated, and evidence for both balancing and directional selection on alleles and lineages of DAB and UBA was observed in populations throughout the species range but rarely on both loci within a population. These results suggest that fluctuating selection has resulted in the divergence of MHC loci in contemporary populations.     

MHC Variability of a Small Lemur in the Littoral Forest Fragments of Southeastern Madagascar

Habitat fragmentation inhibits gene flow between populations often resulting in a loss of genetic diversity with possible negative effects on fitness parameters. In vertebrates, growing evidence suggests that such genetic diversity is particularly important at the level of the major histocompatibility complex (MHC) because its gene products play an important role in immune functions. Diversity in the MHC is assumed to improve population viability. Here, we investigated the impact of forest fragmentation on the genetic variability of one of the functionally important parts of the MHC, DRB exon 2, of the endemic mouse lemur Microcebus murinus by comparing populations inhabiting two littoral forest fragments of different size in southeastern Madagascar. Twelve different alleles of DRB exon 2 were found in 145 individuals of M. murinus with high levels of sequence divergence between alleles. In both subpopulations, levels of genetic diversity were high, and the genetic analyses revealed only limited effects of fragmentation. Significantly more non-synonymous than synonymous substitutions were found in the functionally important antigen recognition and binding sites indicating selection processes maintaining MHC polymorphism. This is the first study on MHC variation in a free-ranging Malagasy lemur population.     

Range-wide phylogeography of the European temperate-montane herbaceous plant Meum athamanticum Jacq.: evidence for periglacial persistence

Aim  The aim of this study is to analyse the genetic population structure of Meum athamanticum Jacq. in order to explore the alternative hypotheses (1) that the central and northern highland populations are the result of post-glacial recolonization from southern refugia, and the disjunct distribution of M. athamanticum can be explained by modern ecological conditions, or (2) that extant populations north of the Alps and Pyrenees persisted in situ during glacial periods.
Location  Europe.
Methods  Variation of amplified fragment length polymorphisms (AFLPs) was analysed for 23 populations from the entire range of the species. We used band-based approaches and methods based on allele frequencies to measure genetic diversity and to identify population structure.
Results  Our analyses reveal a north–south differentiation within M. athamanticum . High levels of genetic diversity, as well as private fragments, are found in populations both north and south of the Alps. Differentiation among populations is lower in the northern than in the southern population group, and significant isolation-by-distance is found only in the latter group.
Main conclusions  Our results indicate that M. athamanticum survived the last ice age in multiple refugia throughout its contemporary range and did not expand into areas north of the Alps from southern refugia. We found evidence that regional-scale migration in northern, formerly periglacial, parts of the species range has resulted in the intermingling of populations. In contrast, southern populations are characterized by long-term isolation. The south-west Alps represent an area where immigration and mixing of populations from northern and southern refugia appears to have taken place.     

MHC diversity in bottlenecked populations: a simulation model

The depletion of variation at MHC loci, which play a crucial role in pathogen recognition, has been postulated to be one of important extinction risk factors for endangered populations. Thus, it is important to understand how selection affects the level of polymorphism in these genes when populations undergo a reduction in size. We followed MHC diversity in computer simulations of population bottlenecks. The fates of MHC alleles in the simulations were determined either by drift, or by balancing selection resulting from host–parasite coevolution. We found that the impact of selection on MHC polymorphism in bottlenecked populations was dependent upon the timescales involved. Initially, selection maintained lower number of alleles than drift, but after ~40 generations of hosts selection maintained higher MHC diversity, as compared to drift. The adverse effects of decreased MHC polymorphism on population viability may be, to some extent, compensated for if selection helps to retain MHC alleles which show high functional diversity, which should allow protection against a broader range of pathogens. Our simulation shows, however, that the mean divergence of alleles retained under selection in bottlenecked populations is not, on average, significantly higher than the divergence due to drift.     

Population subdivision in westslope cutthroat trout (Oncorhynchus clarki lewisi) at the northern periphery of its range: evolutionary inferences and conservation implications

Westslope cutthroat trout (Oncorhynchus clarki lewisi, Salmonidae) are native to the upper Columbia, Missouri, and South Saskatchewan river drainages of western North America and are at the northern periphery of their range in southeastern British Columbia, Canada. We examined geographical variation in allele frequencies at eight microsatellite loci in 36 samples of westslope cutthroat trout from British Columbia to assess levels of population subdivision and to test the hypothesis that different habitat types (principally mainstem vs. above migration barrier habitats) would influence levels of genetic diversity, genetic divergence among populations, and attainment of equilibrium between gene flow and genetic drift. Across all samples, the mean number of alleles per locus was 3.9 and mean expected heterozygosity was 0.56. Population subdivision was extensive with an overall Fst (theta) of 0.32. Populations sampled above migration barriers had significantly fewer alleles, lower expected heterozygosity, but greater average pairwise Fst than populations sampled from mainstem localities. We found evidence for isolation-by-distance from a significant correlation between genetic distance and geographical distance (r = 0.31), but the pattern was much stronger (r = 0.51) when above barrier populations and a population that may have been involved in headwater exchanges were removed. By contrast, isolation-by-distance was not observed when only above barrier populations were tested among themselves. Our data support the maintenance of separate demographic management strategies for westslope cutthroat trout inhabiting different river systems and illustrate how differing habitat structure (e.g. presence of migration barriers) may influence patterns of biodiversity and gene flow-drift equilibrium.     

High genetic diversity in Sarracenia leucophylla(Sarraceniaceae), a carnivorous wetland herb

Eighteen allozyme loci were used to examine genetic diversity in 10 natural populations of Sarracenia leucophylla Raf., a pitcher plant restricted to the southeastern United States. One ex situ population propagated for restoration in Georgia was also analyzed. S. leucophylla is an insect-pollinated, outcrossing perennial wetland herb that is threatened over much of its geographic range. Fifteen loci (83.3%) were polymorphic, with a mean number of alleles of 3.33. Compared to species having similar life-history traits and to previously analyzed Sarracenia species, S. leucophylla displayed unexpectedly high genetic diversity. For example, genetic diversity within the species (Hes) was 0.224 and mean population genetic diversity (Hep) was 0.183. Although small S. leucophylla populations maintained less genetic diversity than larger ones, these differences were not statistically significant. Nonetheless, this suggests that small populations may have lost rare alleles. Statistically significant genetic differentiation among populations was found (theta = 0.192, P < .01), although it was not atypical considering the species' life-history characteristics. A significant correlation (P < .01) between genetic and geographic distance was found, indicating an isolation-by-distance effect. However, the correlation coefficient for this relationship was low (r = 0.46), suggesting that factors other than gene flow play a prominent role in the geographic distribution of genetic diversity within the species. The ex situ population captured most of the allozyme variation found in its source population.     

At the interface of phylogenetics and population genetics, the phylogeography of Dirca occidentalis (Thymelaeaceae)

Dirca occidentalis is a rare shrub indigenous to only six counties near the San Francisco Bay in California, United States. We used intersimple sequence repeat (ISSR) markers and automated genotyping to probe 29 colonies of D. occidentalis from four geographically disjunct populations (East Bay, North Bay, Salmon Creek, and Peninsula) and used methods of phylogenetics and population genetics to model variation across the species. Results show that the four disjunct populations are genetically isolated and have undergone divergence. Phylogenetic analyses indicate that the East Bay population was the first to diverge, followed by the North Bay, then the Salmon Creek and Peninsula populations. This order of divergence suggests an intriguing natural history for D. occidentalis that is explained by the dynamic geological and climatic history of the Bay Area. Spatial genetic structure detected for the species suggests an interaction of four factors: limited seed dispersal, clonal regeneration, distances traveled by pollinators, and genetic isolation of the four populations. Genetic diversity within the North Bay and Salmon Creek populations is low, indicating poor ecological fitness and risk of decline. ISSRs resolved phylogeographic structure within D. occidentalis, results unattainable with ITS methods, and the integration of tools of phylogenetics and population biology led to an enhanced understanding of this endemic species.     

Long‐term ‘islands’ in the landscape: low gene flow,effective population size and genetic divergence in the shrub Hakea oldfieldii (Proteaceae)

The genetic structure of disjunct populations is determined by founding genetic properties, demographic processes, gene flow, drift and local selection. We aim to identify the genetic consequences of natural population disjunction at regional and local scales in Hakea oldfieldii using nuclear and plastid markers to investigate long‐term effective population sizes and gene flow, and patterns of diversity and divergence, among populations. Regional divergence was significant as shown by a consistent pattern in principal coordinates, neighbor‐joining and Bayesian analyses, but divergence at the local level was also significant with localized distribution of plastid haplotypes and populations clustering separately in Bayesian analyses. Historical, recent and first‐generation gene flow was low, suggesting that recent habitat fragmentation has not reduced gene migration significantly. Genetic bottlenecks were detected in three populations. Long‐term effective population size was significantly correlated with the number of alleles/locus and observed heterozygosity, but not with census population size, suggesting that the loss of diversity is associated with long‐term changes rather than recent fragmentation. Inbreeding coefficients were significant in only three populations, suggesting that the loss of diversity is linked to drift and bottlenecks associated with demographic processes (local extinction by fires) rather than inbreeding. Historical disjunction as a result of specific ecological requirements, contraction of habitats following drying during the Pleistocene, low gene flow and changes in population size are likely to have been important forces driving divergence through isolation by distance and drift. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 179 , 319–334.     

Population genomic evidence for plant glacial survival in Scandinavia

Quaternary glaciations have played a major role in shaping the genetic diversity and distribution of plant species. Strong palaeoecological and genetic evidence supports a postglacial recolonization of most plant species to northern Europe from southern, eastern and even western glacial refugia. Although highly controversial, the existence of small in situ glacial refugia in northern Europe has recently gained molecular support. We used genomic analyses to examine the phylogeography of a species that is critical in this debate. Carex scirpoidea Michx subsp. scirpoidea is a dioecious, amphi‐Atlantic arctic–alpine sedge that is widely distributed in North America, but absent from most of Eurasia, apart from three extremely disjunct populations in Norway, all well within the limits of the Weichselian ice sheet. Range‐wide population sampling and variation at 5,307 single nucleotide polymorphisms show that the three Norwegian populations comprise unique evolutionary lineages divergent from Greenland with high between‐population divergence. The Norwegian populations have low within‐population genetic diversity consistent with having experienced genetic bottlenecks in glacial refugia, and host private alleles that probably accumulated in long‐term isolated populations. Demographic analyses support a single, pre‐Weichselian colonization into Norway from East Greenland, and subsequent divergence of the three populations in separate refugia. Other refugial areas are identified in North‐east Greenland, Minnesota/Michigan, Colorado and Alaska. Admixed populations in British Columbia and West Greenland indicate postglacial contact. Taken together, evidence from this study strongly indicates in situ glacial survival in Scandinavia.     

Genetic wealth,population health: Major histocompatibility complex variation in captive and wild ring‐tailed lemurs (Lemur catta)

Across species, diversity at the major histocompatibility complex (MHC) is critical to individual disease resistance and, hence, to population health; however, MHC diversity can be reduced in small, fragmented, or isolated populations. Given the need for comparative studies of functional genetic diversity, we investigated whether MHC diversity differs between populations which are open, that is experiencing gene flow, versus populations which are closed, that is isolated from other populations. Using the endangered ring‐tailed lemur (Lemur catta) as a model, we compared two populations under long‐term study: a relatively “open,” wild population (n = 180) derived from Bezà Mahafaly Special Reserve, Madagascar (2003–2013) and a “closed,” captive population (n = 121) derived from the Duke Lemur Center (DLC, 1980–2013) and from the Indianapolis and Cincinnati Zoos (2012). For all animals, we assessed MHC‐DRB diversity and, across populations, we compared the number of unique MHC‐DRB alleles and their distributions. Wild individuals possessed more MHC‐DRB alleles than did captive individuals, and overall, the wild population had more unique MHC‐DRB alleles that were more evenly distributed than did the captive population. Despite management efforts to maintain or increase genetic diversity in the DLC population, MHC diversity remained static from 1980 to 2010. Since 2010, however, captive‐breeding efforts resulted in the MHC diversity of offspring increasing to a level commensurate with that found in wild individuals. Therefore, loss of genetic diversity in lemurs, owing to small founder populations or reduced gene flow, can be mitigated by managed breeding efforts. Quantifying MHC diversity within individuals and between populations is the necessary first step to identifying potential improvements to captive management and conservation plans.     

A comparison of biallelic markers and microsatellites for the estimation of population and conservation genetic parameters in Atlantic salmon (Salmo salar)

Biallelic markers such as single nucleotide polymorphisms (SNPs) and insertion/deletion polymorphisms have become increasingly popular markers for various population genetics applications. However, the effort required to develop biallelic markers in nonmodel organisms is still substantial. In this study, we compared the estimation of various population genetic parameters (genetic divergence and structuring, isolation-by-distance, genetic diversity) using a limited number of biallelic markers (in total 7 loci) to those estimated with 14 microsatellite loci in 21 Atlantic salmon (Salmo salar) populations from northern Europe. Pairwise FST values were significantly correlated between biallelic loci and microsatellite datasets, as was overall heterozygosity when both anadromous and nonanadromous populations were analyzed together. However, when the anadromous and nonanadromous samples were analyzed separately, only genetic divergence correlations remained significant. Biallelic markers alone were not sufficient for reliable neighbor-joining clustering of populations but gave highly similar isolation-by-distance signals when compared with microsatellites. Finally, although several population prioritization measures for conservation exhibited significant correlation between different marker types, the specific populations highlighted as being most valuable for conservation purposes varied depending on the marker type and conservation criteria applied. This study demonstrates that a relatively small set of biallelic markers can be sufficient for obtaining concordant results in most of the analyses compared with microsatellites, although estimates of genetic distance are generally more concordant than estimates of genetic diversity. This suggests that a relatively small number of biallelic markers can provide useful information for various population genetic applications. However, we emphasize that the use of much higher number of loci is preferable, especially when the genetic differences between populations are subtle or individual multilocus genotype-based analyses are to be performed.     

Spatial and temporal patterns of genetic diversity in an endangered Hawaiian honeycreeper,the Hawaii Akepa (<Emphasis Type="Italic">Loxops coccineus coccineus</Emphasis>)

As a result of disease, habitat destruction, and other anthropogenic factors, the Hawaii Akepa (Loxops coccineus coccineus) currently occupies <10% of its original range and exists in five disjunct populations, raising concerns about what effect such reduction and fragmentation has had on the connectivity and diversity of Akepa populations. In this study, we used both historical and contemporary samples to assess genetic diversity and structure in this endangered Hawaiian honeycreeper. We generated sequence data from two mtDNA regions (ND2, control region) and two nuclear introns for contemporary samples representing three of the five current populations. We also generated control region sequence data for museum specimens collected over 100 years ago from throughout the historical range of the bird. Results indicate that despite recent declines and fragmentation, genetic diversity has not been lost. We detected a modest level of genetic differentiation, which followed a combined pattern of isolation-by-barriers and isolation-by-distance, across the historical range of Akepa. The similarly low level of differentiation observed between the contemporary populations indicates that not much divergence, if any, has occurred post-fragmentation. Rather, the present structure seen likely reflects the historical pattern of distribution. Ironically, this declining species exhibits the genetic signal of an expanding population, demonstrating that earlier demographic events are outweighing the effects of recent changes in population size, and genetic estimates of N_e, though crude, suggest Hawaii Akepa were at least an order of magnitude more abundant prior to the decline.     

The STR polymorphism (AAAAT)n within the intron 1 of the tumor protein 53 (TP53) locus in 17 populations of different ethnic groups of Africa, America, Asia and Europe

The STR (AAAAT)n within intron 1 of the TP53 locus was screened in 17 populations from 3 main ethnic groups: Europeans, Asiatics, and Africans, and from the hybrid population of Costa Rica (1968 samples). Three alleles, 126/7 (bp/copies of the repeat), 131/8 and 136/9 were the most prevalent in all populations. Other alleles rarely reached frequencies of 10% or higher. Observed heterozygosities ranged between 0.351 and 0.829. Patterns of diversity fit well with both the geographic origin of the samples and the history of the populations screened. A statistical test suggests that single-step mutational events have been the main mechanism producing new alleles at this locus. Fixation indexes (R(ST)) for this marker showed an effect of population subdivision on divergence only within the Asiatic group; they were insensitive at the level of major ethnic groups as well as within Africans and within Europeans.     

The structure of allozyme and leaf shape variation in isolated, range-margin populations of the shrub Hippocrepis emerus (Leguminosae)

We investigated allozyme variation and phenotypic variation in leaf shape in 15 populations of the shrub Hippocrepis emerus (Leguminosae) from the three isolated, regional populations in Scandinavia and analysed patterns of differentiation and the hierarchical structuring of diversity on different geographic scales There are pronounced geographic differences between the Norwegian and Swedish isolates of the species and most of the polymorphic allozyme loci show reciprocal fixation in the two isolates The Scandinavian populations of H emerus are not only disjunct on a macrogeographic scale but also show considerable disjunction within regions Within the Oland regional population, a central group of populations shows low levels of Inter-population differentiation the Gotland group of populations is related to this core group of Oland populations The geographically marginal populations on Oland are spatially isolated and show a higher degree of divergence between populations than does the central group of Oland populations We interpret genetic divergence between the marginal populations in terms of genetic drift - as a result of historical fluctuations in habitat availability and population size     

Lost but not forgotten: MHC genotypes predict overwinter survival despite depauperate MHC diversity in a declining frog

The amphibian disease chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), has contributed to the decline of Chiricahua leopard frogs (Rana chiricahuensis), a federally threatened species native to the Southwestern United States. We characterized immunogenetic variability in R. chiricahuensis by sequencing an expressed Major Histocompatibility Complex (MHC) class IIβ gene across 13 natural populations in Arizona, USA, as well as 283 individuals that were captive reared from two egg masses. We recovered a total of five class IIβ MHC alleles compared to 84 alleles previously characterized in eight natural populations of the Arizona congener R. yavapaiensis, demonstrating reduced MHC diversity in R. chiricahuensis. One allele was fixed in five populations but none of the R. chiricahuensis alleles were closely related to R. yavapaiensis allele Q, which is significantly associated with chytridiomycosis resistance in laboratory trials. Nine of 13 R. chiricahuensis population localities were Bd positive, and bearing allele RachDRB*04 was the best genetic predictor of an individual being infected with Bd. A total of three class IIβ alleles were recovered from captive reared individuals, which were released to two natural population localities followed by recapture surveys to assess MHC-based survival over winter, the time when chytridiomycosis outbreaks are most severe. At one site, all released animals were fixed for a single allele and MHC-based survival could not be assessed. At the second site, fewer than half of the released but all of the recaptured individuals were homozygous for RachDRB*05, indicating that MHC genotype is important in determining Bd survival under natural field conditions. We conclude that the limited MHC variation in R. chiricahuensis is likely the consequence rather than the cause of natural selection favoring alleles that promote survival in the face of Bd. Our study highlights that preserving even low levels of functional genetic variation may be essential for population persistence, and that local disease adaptation may present as a reduction in genetic diversity. These finding also suggest that for populations that have declined due to a specific infectious pathogen, MHC-based genetically-informed reintroduction approaches may enhance species recovery efforts.     

Plasmodium relictum infection and MHC diversity in the house sparrow (Passer domesticus)

Antagonistic coevolution between hosts and parasites has been proposed as a mechanism maintaining genetic diversity in both host and parasite populations. In particular, the high level of genetic diversity usually observed at the major histocompatibility complex (MHC) is generally thought to be maintained by parasite-driven selection. Among the possible ways through which parasites can maintain MHC diversity, diversifying selection has received relatively less attention. This hypothesis is based on the idea that parasites exert spatially variable selection pressures because of heterogeneity in parasite genetic structure, abundance or virulence. Variable selection pressures should select for different host allelic lineages resulting in population-specific associations between MHC alleles and risk of infection. In this study, we took advantage of a large survey of avian malaria in 13 populations of the house sparrow (Passer domesticus) to test this hypothesis. We found that (i) several MHC alleles were either associated with increased or decreased risk to be infected with Plasmodium relictum, (ii) the effects were population specific, and (iii) some alleles had antagonistic effects across populations. Overall, these results support the hypothesis that diversifying selection in space can maintain MHC variation and suggest a pattern of local adaptation where MHC alleles are selected at the local host population level.