Abstract
MHC genes are key components in disease resistance and an excellent system for studying selection acting on genetic variation in natural populations. Current patterns of variation in MHC genes are likely to be influenced by past and ongoing selection as well as demographic fluctuations in population size such as those imposed by post-glacial recolonization processes. Here, we investigated signatures of historical selection and demography on an MHC class II gene in 12 moor frog populations along a 1700-km latitudinal gradient. Sequences were obtained from 207 individuals and consecutively assigned into two different clusters (northern and southern clusters, respectively) in concordance with a previously described dual post-glacial colonization route. Selection analyses comparing the relative rates of non-synonymous to synonymous substitutions (dN/dS) suggested evidence of different selection patterns in the northern and the southern clusters, with divergent selection prevailing in the south but uniform positive selection predominating in the north. Also, models of codon evolution revealed considerable differences in the strength of selection: The southern cluster appeared to be under strong selection while the northern cluster showed moderate signs of selection. Our results indicate that the MHC alleles in the north diverged from southern MHC alleles as a result of differential selection patterns.
The major histocompatibility complex class II plays a pivotal role in binding and presenting peptides derived from extracellular pathogens for recognition on CD4+ lymphocytes and T cell receptors (Janeway et al. 2017). The antigen presentation will elicit an appropriate immune response in the host based on the type of MHC class II molecule coded by different class II alleles (Wood 2006). The MHC class II exon 2 is by far the most polymorphic and diverse coding genetic region in most vertebrates (Hughes 1999, Klein 1986), and variation is especially pronounced in the peptide binding regions (PBRs; Bernatchez and Landry 2003; Piertney and Oliver 2006). This high MHC variability is beneficial for the host in the face of newly emerging diseases (Schemske et al. 2009). As a survival strategy, pathogens evolve to escape being recognized by the most frequent host MHC molecules, which might lead to infections in large parts of the population and ultimately increase the risk of extinction (Hedrick 2002; Wood 2006). As a consequence of this arms race, host populations may express hundreds of different MHC variants, making it more difficult for pathogens to fully avoid detection (Kohn et al. 2006; Wood 2006).
The molecular structure of the MHC has been extensively studied in humans and many other mammalian species (Janeway 2001). However, the MHC structure in other vertebrates, such as amphibians, is less well known. To study the MHC structure in other taxa, researchers commonly compare the amino acid positions of PBRs of the human leukocyte antigen (HLA; the human MHC; Bondinas et al. 2007; Brown et al. 1993; Tong et al. 2004) to infer sites involved in peptide binding. The variability at the PBR is believed to be subject to strong selective pressure and contains signatures of both positive and negative selection (Meyer and Thomson 2001).
In Scandinavian moor frogs (Rana arvalis), genetic variation at the MHC class II exon 2 is shaped by a complex pattern of past and ongoing selection, drift, and/or historical demographic events including post-glacial recolonization history where selection is a relatively more important process in the south whereas drift is more important in the north (Cortázar-Chinarro et al. 2017). These processes can increase or decrease genetic variation and explain the observed patterns of current genetic variation in nature (Hedrick 1998). Several studies have found that MHC II exon 2 can be under divergent selection among different populations, suggesting adaptation to differences in the local parasite fauna (e.g., Ekblom et al. 2007; Meyer-Lucht et al. 2016). Others have found patterns indicating long-term balancing selection maintaining trans-species polymorphism (e.g., Vlček et al. 2016). Furthermore, other authors showed positive selection and trans-species polymorphism working at the same time (Surridge et al. 2008). However, none of the studies have explicitly examined the relationships between the patterns of historical selection at MHC at latitudinal scales, and how MHC variation can be shaped by drift, selection, and colonization processes.
Amphibian populations are suffering worldwide declines and are the most threatened vertebrate taxon on the planet (Stuart et al. 2004). Habitat fragmentation in combination with infectious diseases, such as chytridiomycosis caused by the fungus Batrachochytrium dendrobatidis (Bd), has been identified as the most important factor in the amphibian decline (Berger et al. 1998; Garner et al. 2006; Lips et al. 2006; Morgan et al. 2007). A number of studies in amphibians suggest that specific MHC class II molecules are involved in immunity against Bd and other disease agents (Babik et al. 2008, 2009; Lillie et al. 2015; Zeisset and Beebee 2009, 2013), suggesting a link between MHC II variability and the susceptibility or resistance to fungi and other disease agents (Bataille et al. 2015; Savage and Zamudio 2011).
In this study, we explore the mode of historical selection in 12 R. arvalis populations along a 1700-km latitudinal gradient from northern Germany to northern Sweden. Previous phylogeographic studies using microsatellites, mtDNA, and allelic variation at the MHC II exon 2 gene in R. arvalis indicate two post-glacial recolonization routes into Scandinavia (Babik et al. 2004; Cortázar-Chinarro et al. 2017; Knopp and Merilä 2009). Under this scenario, populations were assigned belonging to either of two different clusters (northern and southern clusters, respectively) suggesting that individuals from the two lineages evolved in allopatry until they came into secondary contact ca. 5000 years ago. Our ultimate goal is to investigate the mechanisms underlying the historical formation and maintenance of MHC II exon 2 variation along the two post-glacial colonization routes.
We used the MHC class II exon 2 data set described in Cortázar-Chinarro et al. (2017) for analyzing signatures of historical selection in natural populations of the moor frog (R. arvalis) along the latitudinal gradient (see supplementary material; Table S1). DNA extraction methods, primer design, MHC genotyping, and data filtration are detailed in Cortázar-Chinarro et al. (2017). In this previous study, we assessed contemporary selection primarily based on general allelic frequencies and a classical FST outlier approach. We identified the MHC II exon 2 (corresponding to the β-2 domain) to be subject to diversifying selection, while five microsatellite loci showed signals of stabilizing selection among populations. In general, contemporary signs of selection acting on MHC were weaker in the north, possibly as a consequence of the effects of demography. In the present study, we conduct calculations for inferring signatures of historical selection on the MHC sequences, dN/dS, ratios and Tajima’s D, as well as dN/dS per codon in a sliding window approach on each position, looking at selection from a molecular point of view.
We used MEGA v.7.0 to visualize the data and align haplotype sequences of the MHC class II gene exon 2β chain (Kumar et al. 2016). To determine whether the sequences constitute classical functional class II alleles, we looked for the presence of indels causing a shift in the reading frame and/or stop codons. The peptide binding region was defined according to the human HLA (Bondinas et al. 2007). Some of the codons are generally highly conserved among species in classical MHC class II molecules (Bondinas et al. 2007; Brown et al. 1993). We used the MHC class II exon 2 (272 bp) data set from Cortázar-Chinarro et al. (2017) for analyzing natural populations of the moor frog (R. arvalis; total 207 individuals, 12 populations) along a latitudinal gradient (77 individuals from the north, Umeå and Luleå, and 130 from the south, Germany, Skåne, and Uppsala, Table S1). We amplified and sequenced the complete exon 2 (272 bp) at a single MHC II locus. We detected 57 unique haplotypes over the entire gradient; 4 haplotypes were shared between the northern and southern cluster (Table 1; Fig. 1; see Cortázar-Chinarro et al. 2017). Each of the 57 haplotypes translated into 90 amino acids, none of the allele sequences contained an indel or a stop codon, and we thus assume that they all are functional.
Nucleotide diversity (П), number of segregating sites (S), and average number of pairwise nucleotide differences (theta K) were calculated in DNAsp v5.0 (Librado and Rozas 2009). Observed heterozygosity (HO), expected heterozygosity (HE), and allelic richness (AR) were calculated in Fstat (Goudet 2001). Deviations from Hardy-Weinberg equilibrium were tested for in Genepop (Rousset 2008). All the genetic diversity measures were calculated for the entire gradient, and for the northern and the southern cluster separately. The number of segregating sites (S) was 31 (11% of the 272 nucleotide sites) in the entire gradient and in the southern cluster, and 19 in the northern cluster (7% of the 272 nucleotide sites). MHC II exon 2 revealed a higher nucleotide diversity (П) and higher average number of pairwise differences (theta K) in the entire gradient and in the south compared to the north (Table 1).
We constructed a phylogenetic tree to illustrate the phylogenetic relationship among R. arvalis MHC sequences using the Neighbor joining method with bootstrapping (1000 replicates, Fig. S1) implemented in MEGA v7.0 (Kumar et al. 2016). We also reconstructed an unrooted phylogenetic network in the software SplitTree4, using the neighbor net method (Huson and Bryant 2005, Fig. 1). In both analyses, we included a natterjack toad (Epidalea calamita; GenBank: HQ388291.1) sequence as an outgroup. We further constructed a haplotype network in the Hapstar software (Teacher and Griffiths 2011), which shows the number of base-pair changes over the sequences (Fig. S2). The input file was generated by selecting the minimum spanning tree option in ARLEQUIN software (Excoffier and Lischer 2010). Statistical network analyses were run with the “igraph” package in R (Csardi and Nepusz 2006). All three approaches showed that the MHC II exon 2 sequences from the northern groups are clearly non-randomly distributed within the network, indicating their sequence similarity among each other (Figs. 1, S1, and S2). Only four out of the 57 haplotypes were shared between the northern and southern clusters (Raar_DAB*15, Raar_DAB*18, Raar_DAB*20, Raar_DAB*22). Among these four haplotypes, Raar_DAB*15 was found in nine out of 12 populations along the gradient and it is central in the haplotype network (Fig. S2).
We used a number of approaches to investigate whether regions of the R. arvalis MHC class II sequences carry signatures of selection (i.e., over evolutionary time), either as positive or negative selection (Piertney and Oliver 2006). Positive selection is the process by which new advantageous alleles are spread to fixation and therefore increase the fitness of individuals in a population. In contrast, negative selection is the selective removal of deleterious alleles which ultimately reduces variation in a population. In MHC genes, positive selection may suggest adaptive responses to pathogen diversity while negative selection or purifying selection may also suggest that the MHC fragment plays a functional role independent of diseases or microbial resistance by removal of deleterious alleles in a population and reducing the frequency of mutations in a pool of linked genes (see e.g., Charlesworth 2006).
We used DNAsp (Librado and Rozas 2009) to calculate Tajima’s D, a statistical method for testing for deviation from neutral evolution of DNA sequences, on the MHC II exon 2 locus. While we detected negative Tajima’s D values in the northern cluster, and positive values of Tajima’s D in the entire gradient and southern clusters (Table 1), the lack of significance (p > 0.10) indicates that deviation from neutrality could not be detected by this approach.
To detect historical selection (i.e., over evolutionary time), we used a common method based on the ratio of non-synonymous to synonymous nucleotide substitutions (dN/dS or ω). Ratios of dN/dS > 1 indicate positive selection and dN/dS < 1 indicate negative purifying selection (Garrigan and Hedrick 2003; Hughes 1999, 2007). We used two different methods to calculate ratios of synonymous to non-synonymous substitutions. In the first approach, we estimated the average synonymous (dS) and non-synonymous (dN) substitutions per synonymous and non-synonymous site as implemented in MEGA v.7.0 (Kumar et al. 2016). We used the Nei-Gojobori/Jukes-Cantor method with 5000 bootstrap replicates to calculate the overall average of dN/dS (ω) for all sites, the peptide binding region (PBR), and non-PBR according to Bondinas et al. (2007). Deviations from neutrality (dN ≠ dS) were assessed by a Z test. Differences in dN and dS between PBR and non-PBR were tested for significance using Mann-Whitney U tests. In addition, we identified signatures of positive selection using OmegaMap as implemented by Wilson and McVean (2006).
We estimated dN and dS at the PBR and the non-PBR (Bondinas et al. 2007) for the entire gradient, the northern and the southern clusters (Table 2; Fig. S3). For the PBR, the values for the ratio (ω) were larger than 1 in all the groups (entire gradient, northern and southern subgroups). The Z test for selection indicated that ω was significantly different from neutral expectations in all groups (entire gradient Z = 3.894; northern cluster; Z = 2.672 and southern cluster = 3.621; p < 0.001) (Table 2; Fig. 2). This is considered as evidence of positive selection acting on the PBR. Codons outside the putative PBR regions showed a ω value < 1 in the entire gradient and in the south, suggesting negative selection; however, neither of these were significant (Table 2). In the north, ω was > 1, but not significant (Table 2). Values of ω > 1 may indicate positive selection and/or effects of demography.
A signature of natural selection was also detected using the dN/dS (ω) ratio and a signature of recombination (Rho; ρ) from the patterns of linkage disequilibrium by Bayesian interference in OmegaMap (Wilson and McVean 2006). As recommended, the model was run twice and then two runs were combined with “burn-in” of 50,000 iterations using the Summarize module provided by OmegaMap. In each simulation, we run 500,000 iterations and thinned every 1000 iterations to obtain the posterior distribution. We allowed for ω and ρ to vary across the codons as suggested by the manual, set the ω prior and ρ prior = inverse (0.010, 100), and ω and ρ model = independent. The remaining priors were set to improver inverse. In addition, we performed a sliding window analysis to visualize the position of the codons under selection (ω > 1) in MHC II exon 2 in concordance to OmegaMap results. All selection analyses were run separately for the northern and southern clusters. Analyses provided significant evidence of positive selection (ω > 1) in several positions (Table 2; Fig. 2 and Fig. S3). We found 11 codons in the southern cluster (positions 11, 13, 27, 28, 34, 47 60, 63, 70, 74, 78) and six codons in the northern cluster (positions 13, 37, 47, 60, 63, 74) showing signatures of positive selection with significant posterior probabilities (prob. of selection > 95%). Codons 11, 13, 27, 28, 37, 47, 70, 74, and 78 are part of the putative PBR of the MHC class II (Bondinas et al. 2007; marked in Table 2). While some of the codons exhibiting evidence of positive selection were not identified as part of the putative PBR (Bondinas et al. 2007), these codons were located close to the putative PBR sites (Table 2). Different signatures of natural selection between the northern and the southern clusters were also detected by the SELECTON server (Stern et al. 2007) and CodeML method (See supplementary Table S2). Both approaches are implemented within the Phylogenetic Analysis by Maximum Likelihood software package (PAML, version 4.9d) (Yang 2007). We followed a conservative criterion to select the codons identified at least by two methods to consider them as positive selective codons (Supplementary Table S2; see, e.g., Wlasiuk and Nachman 2010). All the codons selected for positive selection in OmegaMap were also highlighted in the two other methods. According to OmegaMap, there was no evidence of recombination in the MHC class II sequences (lower 95% CI ρ > 1). Moreover, as the power to detect recombination varies with different detection methods, we used an additional array of six methods from the RDP package (Martin et al. 2010) to analyze recombination (Supplementary Table S3). These revealed very weak indications of recombination, as four methods did not find recombination signals at all and two methods revealed a single recombination event in the south and the entire gradient.
Based on earlier studies on microsatellites and mitochondrial markers indicating dual recolonization routes to Scandinavia from southeastern Europe glacial refugia (Babik et al. 2004; Cortázar-Chinarro et al. 2017; Knopp and Merilä 2009), we assigned R. arvalis populations to two different clusters: a northern cluster consisting of populations from Umeå and Luleå and a southern cluster with populations from Germany, Skåne, and Uppsala. Consistent with this hypothesis, both the phylogenetic tree and the unrooted phylogenetic network in this study illustrated contrasting evolutionary histories in the northern and southern clusters, and we found striking differences in the signature of historical selection in the MHC II exon 2 between the northern and southern parts of the latitudinal gradient.
We found that almost all MHC II exon 2 haplotype sequences from the north clustered in a single clade, with the exception of the Raar_DAB*32, suggesting that they are more similar to each other than to southern alleles. On the contrary, alleles from the southern populations are scattered throughout the tree showing high allelic divergence (see Fig. 1). Alleles from northern and southern populations might have a different demographic origin, as predicted by the dual colonization hypothesis.
In terms of genetic diversity, we found that the MHC II exon 2 locus was more diverse in the southern cluster than in the northern in R. arvalis. As a result of post-glacial recolonization, genetically less diverse populations are often found at northern latitudes, differing from the southern populations due to their demographic history (Hewitt 1996, 1999; Taberlet et al. 1998). We found a total of 47 MHC II haplotypes in the south, 14 in the north and only four haplotypes which were shared between the north and south (Raar_DAB*18, Raar_DAB*22, Raar_DAB*9, Raar_DAB*15). The haplotype Raar_DAB*15 is present in Umeå but not further north in Luleå and in nine out of 12 populations along the gradient. It is a central haplotype in the network and connects the northern and southern clusters. This haplotype, very common in the southern cluster, could have arrived to the northern populations via introgression from the south across a hybrid zone. The average number of pairwise nucleotide differences (theta K) between the alleles is very high, indicating a high allelic divergence in the southern cluster (theta K = 11.7). It is considerably lower (theta K = 4.4) in the northern cluster, indicating that a majority of MHC II exon 2 alleles are very similar to each other in the north. This might suggest that MHC exon 2 is under uniform positive selection in the north, decreasing genetic variation.
Biotic interactions such as host-parasite interactions are believed to be of great importance for the origin and maintenance of genetic variation in latitudinal gradients. The number of parasite species increases towards the tropics, and there is evidence that parasites are more prevalent and virulent at low latitudes (Schemske et al. 2009). Smaller and more isolated R. arvalis populations at northern latitudes displaying a lower MHC II diversity might be more vulnerable to new emerging pathogens and, consequently, at higher risk of extinction as suggested in previous studies on mammals confronted with new pathogens due to global warming (e.g., Siddle et al. 2007). There is a need for further studies linking patterns of genetic variation and divergence with occurrence of infectious disease across latitudinal and other environmental gradients to better infer evolutionary processes shaping adaptive genetic variation in wildlife populations.
The high levels of polymorphism observed at most genes of the MHC are generated and maintained by means of balancing selection (Hughes and Nei 1988). In this study, we found deviating signals of positive selection at the MHC II exon 2 in the entire gradient and in the southern cluster compared to the northern cluster. While the high levels of genetic variation and allelic divergence in the southern cluster represent a classical pattern for MHC variation, the pattern of uniform positive selection in the north is deviating from this classical case and rather uncommon. Elevated and highly significant dN/dS ratios were found at the PBR at MHC II exon 2, though signals of historical selection in northern populations were weaker compared to those of southern populations. Consistent with previous studies (Bos and DeWoody 2005; Huang et al. 2016; Kiemnec-Tyburczy et al. 2012), we did not find any evidence that MHC II exon 2 locus has been undergoing recombination. Also, several positions along the MHC class II exon 2 showed evidence of positive selection, as revealed by OmegaMap. However, these positively selected positions deviated among the southern and the northern clusters. For the northern cluster, we found six codons under positive selection, whereas 11 codons were under positive selection for the southern cluster. This could mean that in spite of a strong signal of historical selection in the north, the effects of positive selection/drift could be acting on MHC and continue to act if the populations remain small and isolated, reflecting a differential selection patterns among northern and southern clusters.
Our earlier study (Cortázar-Chinarro et al. 2017) disentangling contemporary evolutionary forces acting on MHC II exon 2 along the same latitudinal gradient suggested that the MHC II exon 2 shows signs of contemporary diversifying selection in the southern populations but drift and uniform selection are dominant forces in the northern populations. In a similar example, a study on two mouse lemur species showed contrasting selection patterns acting on MHC II exon 2 locus, providing strong evidence of historical selection in Microcebus murinus, while selection was more relaxed for M. berthae (Pechouskova et al. 2015). The different selection patterns found in our study could be mediated by differential pathogen-driven selection processes together with the effects of the effective population size and the demographic history along the gradient. Although our results support differential selection patterns found in the previous studies on MHC II exon 2 variation (Cortázar-Chinarro et al. 2017; Pechouskova et al. 2015), we cannot rule out the possibility that our results may be partly due to lower sample sizes and fewer population contrasts in the north. Consequently, there is a need for additional studies with stronger sampling effort to further our understanding how genetic variation in R. arvalis is geographically partitioned and distributed from the Balkans to the northern part of Finland.
The results of the current study suggest that selection patterns differ from northern and southern clusters in R. arvalis. The MHC class II exon 2 has fewer sites subject to positive selection in the northern populations compared to the southern populations. We believe that the differential selection patterns observed in our study could be explained by the historical post-glacial recolonization processes and by the different selective pressures from parasites in the different populations along the present gradient and/or lower parasite diversity at northern latitudes, which should be investigated in more detail in future studies.
References
Babik W, Branicki W, Sandera M, Litvinchuk S, Borkin L, Irwin J, Rafiński J (2004) Mitochondrial phylogeography of the moor frog, Rana arvalis. Mol Ecol 13(6):1469–1480. https://doi.org/10.1111/j.1365-294X.2004.02157.x
Babik W, Pabijan M, Radwan J (2008) Contrasting patterns of variation in MHC loci in the Alpine newt. Mol Ecol 17(10):2339–2355. https://doi.org/10.1111/j.1365-294X.2008.03757.x
Babik W, Pabijan M, Arntzen J, Cogalniceanu D, Durka W, Radwan J (2009) Long-term survival of a urodele amphibian despite depleted major histocompatibility complex variation. Mol Ecol 18(5):769–781. https://doi.org/10.1111/j.1365-294X.2008.04057.x
Bataille A, Cashins SD, Grogan L, Skerratt LF, Hunter D, McFadden M, Scheele B, Brannelly LA, Macris A, Harlow PS (2015) Susceptibility of amphibians to chytridiomycosis is associated with MHC class II conformation. Proc R Soc B 282:1805
Berger L, Speare R, Daszak P, Green DE, Cunningham AA, Goggin CL, Slocombe R, Ragan MA, Hyatt AD, McDonald KR (1998) Chytridiomycosis causes amphibian mortality associated with population declines in the rain forests of Australia and Central America. Proc Natl Acad Sci U S A 95(15):9031–9036. https://doi.org/10.1073/pnas.95.15.9031
Bernatchez L, Landry C (2003) MHC studies in nonmodel vertebrates: what have we learned about natural selection in 15 years? J Evol Biol 16(3):363–377. https://doi.org/10.1046/j.1420-9101.2003.00531.x
Bondinas GP, Moustakas AK, Papadopoulos GK (2007) The spectrum of HLA-DQ and HLA-DR alleles, 2006: a listing correlating sequence and structure with function. Immunogenetics 59(7):539–553. https://doi.org/10.1007/s00251-007-0224-8
Bos DH, DeWoody JA (2005) Molecular characterization of major histocompatibility complex class II alleles in wild tiger salamanders (Ambystoma tigrinum). Immunogenetics 57(10):775–781. https://doi.org/10.1007/s00251-005-0038-5
Brown MG, Driscoll J, Monaco J (1993) MHC-linked low-molecular mass polypeptide subunits define distinct subsets of proteasomes. Implications for divergent function among distinct proteasome subsets. J Immunol 51:1193–1204
Charlesworth D (2006) Balancing selection and its effects on sequences in nearby genome regions. PLoS Genet 2(4):e64. https://doi.org/10.1371/journal.pgen.0020064
Cortázar-Chinarro M, Lattenkamp EZ, Meyer-Lucht Y, Luquet E, Laurila A, Höglund J (2017) Drift, selection, or migration? Processes affecting genetic differentiation and variation along a latitudinal gradient in an amphibian. BMC Evol Biol 17(1):189. https://doi.org/10.1186/s12862-017-1022-z
Csardi G, Nepusz T (2006) The igraph software package for complex network research. Complex Syst Inter 1695:1–9
Ekblom R, Saether SA, Jacobsson P, Fiske P, Sahlman T, Grahn M, Kålås JA, Höglund J (2007) Spatial pattern of MHC class II variation in the great snipe (Gallinago media). Mol Ecol 16(7):1439–1451. https://doi.org/10.1111/j.1365-294X.2007.03281.x
Excoffier L, Lischer HE (2010) Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour 10(3):564–567. https://doi.org/10.1111/j.1755-0998.2010.02847.x
Garner TW, Perkins MW, Govindarajulu P, Seglie D, Walker S, Cunningham AA, Fisher MC (2006) The emerging amphibian pathogen Batrachochytrium dendrobatidis globally infects introduced populations of the North American bullfrog, Rana catesbeiana. Biol Lett 2(3):455–459. https://doi.org/10.1098/rsbl.2006.0494
Garrigan D, Hedrick PW (2003) Perspective: detecting adaptive molecular polymorphism: lessons from the MHC. Evolution 57(8):1707–1722. https://doi.org/10.1111/j.0014-3820.2003.tb00580.x
Goudet J (2001) FSTAT; a program to estimate and test gene diversities and fixation indices version 2.9. 3. http://www.unil.ch/izea/softwares/fstat
Hedrick PW (1998) Balancing selection and MHC. Genetica 104(3):207–214. https://doi.org/10.1023/A:1026494212540
Hedrick PW (2002) Pathogen resistance and genetic variation at MHC loci. Evolution 56(10):1902–1908. https://doi.org/10.1111/j.0014-3820.2002.tb00116.x
Hewitt GM (1996) Some genetic consequences of ice ages, and their role in divergence and speciation. Biol J Linn Soc 58(3):247–276. https://doi.org/10.1111/j.1095-8312.1996.tb01434.x
Hewitt GM (1999) Post-glacial re-colonization of European biota. Biol J Linn Soc 68(1-2):87–112. https://doi.org/10.1111/j.1095-8312.1999.tb01160.x
Huang L, Zhao M, Luo Z, Wu H (2016) Allelic polymorphism, gene duplication and balancing selection of MHC class IIB genes in the Omei treefrog (Rhacophorus omeimontis). Asian Herpetol Res 7:1–11
Hughes AL (1999) Adaptive evolution of genes and genomes. Oxford University Press, USA
Hughes AL (2007) Looking for Darwin in all the wrong places: the misguided quest for positive selection at the nucleotide sequence level. Heredity 99(4):364–373. https://doi.org/10.1038/sj.hdy.6801031
Hughes AL, Nei M (1988) Pattern of nucleotide substitution at major histocompatibility complex class I loci reveals overdominant selection. Nature 335(6186):167–170. https://doi.org/10.1038/335167a0
Huson DH, Bryant D (2005) Application of phylogenetic networks in evolutionary studies. Mol Biol Evol 23:254–267
Janeway CA (2001) How the immune system works to protect the host from infection: A personal view. Proc Natl Acad Sci 98(13):7461–7468
Janeway CA, Travers P, Walport M, Shlomchik M (2017) Immunobiology: the immune system in health and disease. 2005. Garland Science 6, New York
Kiemnec-Tyburczy K, Richmond J, Savage A, Lips K, Zamudio K (2012) Genetic diversity of MHC class I loci in six non-model frogs is shaped by positive selection and gene duplication. Heredity 109(3):146–155. https://doi.org/10.1038/hdy.2012.22
Klein J (1986) Natural history of the major histocompatibility complex. Wiley, Hoboken
Knopp T, Merilä J (2009) The postglacial recolonization of northern Europe by Rana arvalis as revealed by microsatellite and mitochondrial DNA analyses. Heredity 102(2):174–181. https://doi.org/10.1038/hdy.2008.91
Kohn MH, Murphy WJ, Ostrander EA, Wayne RK (2006) Genomics and conservation genetics. Trends Ecol Evol 21(11):629–637. https://doi.org/10.1016/j.tree.2006.08.001
Kumar S, Stecher G, Tamura K (2016) MEGA7: Molecular Evolutionary Genetics Analysis version 7.0 for bigger datasets. Mol Biol Evol 33(7):1870–1874. https://doi.org/10.1093/molbev/msw054
Librado P, Rozas J (2009) DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics 25(11):1451–1452. https://doi.org/10.1093/bioinformatics/btp187
Lillie M, Grueber CE, Sutton JT, Howitt R, Bishop PJ, Gleeson D, Belov K (2015) Selection on MHC class II supertypes in the New Zealand endemic Hochstetter’s frog. BMC Evol Biol 15(1):63. https://doi.org/10.1186/s12862-015-0342-0
Lips KR, Brem F, Brenes R, Reeve JD, Alford RA, Voyles J, Carey C, Livo L, Pessier AP, Collins JP (2006) Emerging infectious disease and the loss of biodiversity in a Neotropical amphibian community. Proc Natl Acad Sci U S A 103(9):3165–3170. https://doi.org/10.1073/pnas.0506889103
Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P (2010) RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics 26(19):2462–2463. https://doi.org/10.1093/bioinformatics/btq467
Meyer D, Thomson G (2001) How selection shapes variation of the human major histocompatibility complex: a review. Ann Hum Genet 65(1):1–26. https://doi.org/10.1046/j.1469-1809.2001.6510001.x
Meyer-Lucht Y, Mulder KP, James MC, McMahon BJ, Buckley K, Piertney SB, Höglund J (2016) Adaptive and neutral genetic differentiation among Scottish and endangered Irish red grouse (Lagopus lagopus scotica). Conserv Genet 17(3):615–630. https://doi.org/10.1007/s10592-016-0810-0
Morgan JA, Vredenburg VT, Rachowicz LJ, Knapp RA, Stice MJ, Tunstall T, Bingham RE, Parker JM, Longcore JE, Moritz C (2007) Population genetics of the frog-killing fungus Batrachochytrium dendrobatidis. Proc Natl Acad Sci U S A 104(34):13845–13850. https://doi.org/10.1073/pnas.0701838104
Pechouskova E, Dammhahn M, Brameier M, Fichtel C, Kappeler PM, Huchard E (2015) MHC class II variation in a rare and ecological specialist mouse lemur reveals lower allelic richness and contrasting selection patterns compared to a generalist and widespread sympatric congener. Immunogenetics 67(4):229–245. https://doi.org/10.1007/s00251-015-0827-4
Piertney S, Oliver M (2006) The evolutionary ecology of the major histocompatibility complex. Heredity 96(1):7–21. https://doi.org/10.1038/sj.hdy.6800724
Rousset F (2008) genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour 8(1):103–106. https://doi.org/10.1111/j.1471-8286.2007.01931.x
Savage AE, Zamudio KR (2011) MHC genotypes associate with resistance to a frog-killing fungus. Proc Natl Acad Sci U S A 108(40):16705–16710. https://doi.org/10.1073/pnas.1106893108
Schemske DW, Mittelbach GG, Cornell HV, Sobel JM, Roy K (2009) Is there a latitudinal gradient in the importance of biotic interactions? Annu Rev Ecol Evol Syst 40(1):245–269. https://doi.org/10.1146/annurev.ecolsys.39.110707.173430
Siddle HV, Kreiss A, Eldridge MD, Noonan E, Clarke CJ, Pyecroft S, Belov K (2007) Transmission of a fatal clonal tumor by biting occurs due to depleted MHC diversity in a threatened carnivorous marsupial. Proc Natl Acad Sci 104(41):16221–16226. https://doi.org/10.1073/pnas.0704580104
Stern A, Doron-Faigenboim A, Erez E, Martz E, Bacharach E, Pupko T (2007) Selecton: advanced models for detecting positive and purifying selection using a Bayesian inference approach. Nucleic Acids Res 35:506–511
Stuart SN, Chanson JS, Cox NA, Young BE, Rodrigues AS, Fischman DL, Waller RW (2004) Status and trends of amphibian declines and extinctions worldwide. Science 306(5702):1783–1786. https://doi.org/10.1126/science.1103538
Surridge AK, Van der Loo W, Abrantes J, Carneiro M, Hewitt GM, Esteves PJ (2008) Diversity and evolutionary history of the MHC DQA gene in leporids. Immunogenetics 60(9):515–525. https://doi.org/10.1007/s00251-008-0309-z
Taberlet P, Fumagalli L, Wust-Saucy AG, Cosson J-F (1998) Comparative phylogeography and postglacial colonization routes in Europe. Mol Ecol 7(4):453–464. https://doi.org/10.1046/j.1365-294x.1998.00289.x
Teacher A, Griffiths D (2011) HapStar: automated haplotype network layout and visualization. Mol Ecol Resour 11(1):151–153. https://doi.org/10.1111/j.1755-0998.2010.02890.x
Tong JC, Tan TW, Ranganathan S (2004) Modeling the structure of bound peptide ligands to major histocompatibility complex. Protein Sci 13(9):2523–2532. https://doi.org/10.1110/ps.04631204
Vlček J, Hoeck PE, Keller LF, Wayhart JP, Dolinová I, Štefka J (2016) Balancing selection and genetic drift create unusual patterns of MHCIIβ variation in Galápagos mockingbirds. Mol Ecol 25(19):4757–4772. https://doi.org/10.1111/mec.13807
Wilson DJ, McVean G (2006) Estimating diversifying selection and functional constraint in the presence of recombination. Genetics 172(3):1411–1425. https://doi.org/10.1534/genetics.105.044917
Wlasiuk G, Nachman MW (2010) Adaptation and constraint at toll-like receptors in primates. Mol Biol Evol 27(9):2172–2186. https://doi.org/10.1093/molbev/msq104
Wood PJ (2006) Understanding immunology. Pearson Education, London
Yang Z (2007) PAML 4: a program package for phylogenetic analysis by maximum likelihood. Mol Biol Evol 24(8):1586–1591. https://doi.org/10.1093/molbev/msm088
Zeisset I, Beebee T (2009) Molecular characterization of major histocompatibility complex class II alleles in the common frog, Rana temporaria. Mol Ecol Resour 9(3):738–745. https://doi.org/10.1111/j.1755-0998.2009.02535.x
Zeisset I, Beebee TJ (2013) Bufo MHC class II loci with conserved introns flanking exon 2: cross-species amplification with common primers. Conserv Genet Resour 5(1):211–213. https://doi.org/10.1007/s12686-012-9770-y
Acknowledgments
We thank Alex Richter-Boix for invaluable help with the haplotype network reconstruction, the statistical analyses, and for the helpful comments and discussions.
Funding
Funding was provided by Swedish Research Council Formas (grant 215-2014-594 to JH), Stiftelsen Oscar och Lili Lamms Minne (DO2013-0013 to JH), the Swedish Research Council (621-2013-4503 to AL), and Stiftelsen för Zoologisk Forskning.
Author information
Authors and Affiliations
Corresponding author
Electronic supplementary material
Figure S1
Neighbor-joining tree for Rana arvalis MHC II exon 2 representing the 57 different alleles. Alleles present in the north are marked with an orange triangle and shared alleles between the northern and southern cluster are represented with a green rhombus. A natterjack toad sequence [Genbank HQ388291.1] from MHC II exon 2 was used as an outgroup. The scale bar represents substitution per site. (PDF 949 KB)
Figure S2
Haplotype network for the 57 MHC II exon 2 sequences. Every colored circle represents a different haplotype, the dots represent single nucleotide changes. Circles in blue represents haplotypes from the southern cluster, circles in orange represents haplotypes from the northern cluster and green circles represent shared haplotypes between the northern and southern cluster. (PDF 938 KB)
Table S1
(PDF 104 kb)
Table S2
(PDF 29.1 kb)
Table S3
(PDF 89.1 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Cortázar-Chinarro, M., Meyer-Lucht, Y., Laurila, A. et al. Signatures of historical selection on MHC reveal different selection patterns in the moor frog (Rana arvalis). Immunogenetics 70, 477–484 (2018). https://doi.org/10.1007/s00251-017-1051-1
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00251-017-1051-1