Molecular characterization of MHC class II in a nonmodel anuran species, the fire-bellied toad Bombina bombina

Southern introgression increases adaptive immune gene variability in northern range margin populations of Fire-bellied toad

Northern range margin populations of the European fire-bellied toad (Bombina bombina) have rapidly declined during recent decades. Extensive agricultural land use has fragmented the landscape, leading to habitat disruption and loss, as well as eutrophication of ponds. In Northern Germany (Schleswig-Holstein) and Southern Sweden (Skåne), this population decline resulted in decreased gene flow from surrounding populations, low genetic diversity, and a putative reduction in adaptive potential, leaving populations vulnerable to future environmental and climatic changes. Previous studies using mitochondrial control region and nuclear transcriptome-wide SNP data detected introgressive hybridization in multiple northern B. bombina populations after unreported release of toads from Austria. Here, we determine the impact of this introgression by comparing the body conditions (proxy for fitness) of introgressed and nonintrogressed populations and the genetic consequences in two candidate genes for putative local adaptation (the MHC II gene as part of the adaptive immune system and the stress response gene HSP70 kDa). We detected regional differences in body condition and observed significantly elevated levels of within individual MHC allele counts in introgressed Swedish populations, associated with a tendency toward higher body weight, relative to regional nonintrogressed populations. These differences were not observed among introgressed and nonintrogressed German populations. Genetic diversity in both MHC and HSP was generally lower in northern than Austrian populations. Our study sheds light on the potential benefits of translocations of more distantly related conspecifics as a means to increase adaptive genetic variability and fitness of genetically depauperate range margin populations without distortion of local adaptation.

Genetic variation and selection of MHC class I loci differ in two congeneric frogs

Major histocompatibility complex (MHC) genes encode proteins in the acquired immune response pathway that often show distinctive selection-driven patterns in wild vertebrate populations. We examined genetic variation and signatures of selection in the MHC class I alpha 1 (A1)- and alpha 2 (A2)-domain encoding exons of two frog congeners [Agalychnis callidryas (n = 20) and A. lemur (n = 20)] from a single locality in Panama. We also investigated how historical demographic processes may have impacted MHC genetic diversity by analyzing a neutral mitochondrial marker. We found that both MHC domains were highly variable in both species, with both species likely expressing three loci. Our analyses revealed different signatures of selection between the two species, most notably that the A. callidryas A2 domain had experienced positive selection while the A2 domain of A. lemur had not. Diversifying selection acted on the same number of A1 and A2 allelic lineages, but on a higher percentage of A1 sites compared to A2 sites. Neutrality tests of mitochondrial haplotypes predominately indicated that the two species were at genetic equilibrium when the samples were collected. In addition, two historical tests of demography indicated both species have had relatively stable population sizes over the past 100,000 years; thus large population size changes are unlikely to have greatly influenced MHC diversity in either species during this time period. In conclusion, our results suggest that the impact of selection on MHC diversity varied between these two closely related species, likely due to a combination of distinct ecological conditions and past pathogenic pressures.

Limited polymorphism of the functional MHC class II B gene in the black-spotted frog (Pelophylax nigromaculatus) identified by locus-specific genotyping

Amphibians can be more vulnerable to environmental changes and diseases than other species because of their complex life cycle and physiological requirements. Therefore, understanding the adaptation of amphibians to environmental changes is crucial for their conservation. Major histocompatibility complex (MHC) presents an excellent tool for the investigation of adaptive variations and the assessment of adaptive potential because it can be under strong diversifying selection. Here, we isolated the MHC class II B (MHCIIB) gene from cDNA sequences of the black-spotted frog (Pelophylax nigromaculatus), a widespread amphibian species in China, and designed locus-specific primers to characterize adaptive variability of this amphibian. Ten alleles were identified from 67 individual frogs of three populations and no more than two alleles were present in each individual animal. Furthermore, none of the sequences had indels or/and stop codons, which is in good agreement with locus-specific amplification of a functional gene. However, we found low polymorphism at both nucleotide and amino acid levels, even in the antigen-binding region. Purifying selection acting at this locus was supported by the findings that the d_N/d_S ratio across all alleles was below 1 and that negatively selected sites were detected by different methods. Allele frequency distributions were significantly different among geographic populations, indicating that physiographic factors may have strong effect on the genetic structure of the black-spotted frog. This study revealed limited polymorphism of three adjacent black-spotted frog populations at the functional MHCIIB locus, which may be attributed to region-specific differences. The locus-specific genotyping technique developed in this study would provide a foundation for future studies on adaptive divergence among different frog populations.

Evolutionary dynamics of an expressed MHC class IIβ locus in the Ranidae (Anura) uncovered by genome walking and high-throughput amplicon sequencing

The Major Histocompatibility Complex (MHC) is a genomic region encoding immune loci that are important and frequently used markers in studies of adaptive genetic variation and disease resistance. Given the primary role of infectious diseases in contributing to global amphibian declines, we characterized the hypervariable exon 2 and flanking introns of the MHC Class IIβ chain for 17 species of frogs in the Ranidae, a speciose and cosmopolitan family facing widespread pathogen infections and declines. We find high levels of genetic variation concentrated in the Peptide Binding Region (PBR) of the exon. Ten codons are under positive selection, nine of which are located in the mammal-defined PBR. We hypothesize that the tenth codon (residue 21) is an amphibian-specific PBR site that may be important in disease resistance. Trans-species and trans-generic polymorphisms are evident from exon-based genealogies, and co-phylogenetic analyses between intron, exon and mitochondrial based reconstructions reveal incongruent topologies, likely due to different locus histories. We developed two sets of barcoded adapters that reliably amplify a single and likely functional locus in all screened species using both 454 and Illumina based sequencing methods. These primers provide a resource for multiplexing and directly sequencing hundreds of samples in a single sequencing run, avoiding the labour and chimeric sequences associated with cloning, and enabling MHC population genetic analyses. Although the primers are currently limited to the 17 species we tested, these sequences and protocols provide a useful genetic resource and can serve as a starting point for future disease, adaptation and conservation studies across a range of anuran taxa.

Multiple major histocompatibility complex class I genes in Asian anurans: Ontogeny and phylogeny

Amphibians, as the first terrestrial vertebrates, offer a window into early major histocompatibility complex (MHC) evolution. We characterized the MHC class I of two Korean amphibians, the Asiatic toad (Bufo gargarizans) and the Japanese tree frog (Hyla japonica). We found at least four transcribed MHC class I (MHC I) loci, the highest number confirmed in any anuran to date. Furthermore, we identified MHC I transcripts in terrestrial adults, and possibly in aquatic larvae, of both species. We conducted a phylogenetic analysis based on MHC I sequence data and found that B. gargarizans and H. japonica cluster together in the superfamily Nobleobatrachia. We further identified three supertypes shared by the two species. Our results reveal substantial variation in the number of MHC I loci in anurans and suggest that certain supertypes have particular physiochemical properties that may confer pathogen resistance.

Adaptive tolerance to a pathogenic fungus drives major histocompatibility complex evolution in natural amphibian populations

Amphibians have been affected globally by the disease chytridiomycosis, caused by the fungus Batrachochytrium dendrobatidis (Bd), and we are just now beginning to understand how immunogenetic variability contributes to disease susceptibility. Lineages of an expressed major histocompatibility complex (MHC) class II locus involved in acquired immunity are associated with chytridiomycosis susceptibility in controlled laboratory challenge assays. Here, we extend these findings to natural populations that vary both in exposure and response to Bd. We find that MHC alleles and supertypes associated with Bd survival in the field show a molecular signal of positive selection, while those associated with susceptibility do not, supporting the hypothesis that heritable Bd tolerance is rapidly evolving. We compare MHC supertypes to neutral loci to demonstrate where selection versus demography is shaping MHC variability. One population with Bd tolerance in nature shows a significant signal of directional selection for the same allele (allele Q) that was significantly associated with survival in an earlier laboratory study. Our findings indicate that selective pressure for Bd survival drives rapid immunogenetic adaptation in some natural populations, despite differences in environment and demography. Our field-based analysis of immunogenetic variation confirms that natural amphibian populations have the evolutionary potential to adapt to chytridiomycosis.

Susceptibility of amphibians to chytridiomycosis is associated with MHC class II conformation

The pathogenic chytrid fungus Batrachochytrium dendrobatidis (Bd) can cause precipitous population declines in its amphibian hosts. Responses of individuals to infection vary greatly with the capacity of their immune system to respond to the pathogen. We used a combination of comparative and experimental approaches to identify major histocompatibility complex class II (MHC-II) alleles encoding molecules that foster the survival of Bd-infected amphibians. We found that Bd-resistant amphibians across four continents share common amino acids in three binding pockets of the MHC-II antigen-binding groove. Moreover, strong signals of selection acting on these specific sites were evident among all species co-existing with the pathogen. In the laboratory, we experimentally inoculated Australian tree frogs with Bd to test how each binding pocket conformation influences disease resistance. Only the conformation of MHC-II pocket 9 of surviving subjects matched those of Bd-resistant species. This MHC-II conformation thus may determine amphibian resistance to Bd, although other MHC-II binding pockets also may contribute to resistance. Rescuing amphibian biodiversity will depend on our understanding of amphibian immune defence mechanisms against Bd. The identification of adaptive genetic markers for Bd resistance represents an important step forward towards that goal.

Copy number variation and genetic diversity of MHC Class IIb alleles in an alien population of Xenopus laevis

Xenopus laevis (the African clawed frog), which originated through hybridisation and whole genome duplication, has been used as a model for genetics and development for many years, but surprisingly little is known about immune gene variation in natural populations. The purpose of this study was to use an isolated population of X. laevis that was introduced to Wales, UK in the past 50 years to investigate how variation at the MHC compares to that at other loci, following a severe population bottleneck. Among 18 individuals, we found nine alleles based on exon 2 sequences of the Class IIb region (which includes the peptide binding region). Individuals carried from one to three of the loci identified from previous laboratory studies. Genetic variation was an order of magnitude higher at the MHC compared with three single-copy nuclear genes, but all loci showed high levels of heterozygosity and nucleotide diversity and there was not an excess of homozygosity or decrease in diversity over time that would suggest extensive inbreeding in the introduced population. Tajima’s D was positive for all loci, which is consistent with a bottleneck. Moreover, comparison with published sequences identified the source of the introduced population as the Western Cape region of South Africa, where most commercial suppliers have obtained their stocks. These factors suggest that despite founding by potentially already inbred individuals, the alien population in Wales has maintained substantial genetic variation at both adaptively important and neutral genes.

Genetic polymorphism of major histocompatibility complex class IIB alleles and pathogen resistance in the giant spiny frog Quasipaa spinosa

Major histocompatibility complex (MHC) genes are candidates for determining disease susceptibility due to their pivotal role in both innate and adaptive immune responses. Accordingly, the association between the genetic variation of MHC genes and the pathogen resistance has been investigated in numerous vertebrates. To date, however, little is reported in amphibians. In this study, we investigate the genetic variation at the MHC class IIB gene in the giant spiny frog Quasipaa spinosa, which has high commercial value in China. The full length of MHC class IIB cDNA was cloned from Q. spinosa by homology cloning and rapid amplification of cDNA end-polymerase chain reaction (RACE-PCR). Two MHC class IIB loci were identified in Q. spinosa. We also developed PCR primers for a portion of the second exon of the MHC class IIB gene. A total of 26 MHC class IIB alleles were identified. The dN rate was significantly higher than the dS rate in the putative peptide-binding region, thereby proving the positive selection hypothesis. In addition, individuals intraperitoneally injected with Aeromonas hydrophila were used to study the association between MHC class IIB alleles and pathogen resistance/susceptibility, to explore the specific alleles in balancing selection. Eighty frogs were used after exposure to A. hydrophila infection. Nine alleles were used to study the association between the alleles and disease resistance. Two alleles, namely, Pasa-DAB(∗)1301 and Pasa-DAB(∗)0901, were significantly associated with resistance against A. hydrophila. This study provides valuable information on the structure of the MHC class IIB gene and confirms the association between MHC class IIB gene alleles and disease resistance to bacterial infection in Q. spinosa. Moreover, pathogen resistance-related MHC markers can be used for the selective breeding of the giant spiny frog.

Drift rather than selection dominates MHC class II allelic diversity patterns at the biogeographical range scale in natterjack toads Bufo calamita

Study of major histocompatibility complex (MHC) loci has gained great popularity in recent years, partly due to their function in protecting vertebrates from infections. This is of particular interest in amphibians on account of major threats many species face from emergent diseases such as chytridiomycosis. In this study we compare levels of diversity in an expressed MHC class II locus with neutral genetic diversity at microsatellite loci in natterjack toad (Bufo (Epidalea) calamita) populations across the whole of the species' biogeographical range. Variation at both classes of loci was high in the glacial refugium areas (REF) and much lower in postglacial expansion areas (PGE), especially in range edge populations. Although there was clear evidence that the MHC locus was influenced by positive selection in the past, congruence with the neutral markers suggested that historical demographic events were the main force shaping MHC variation in the PGE area. Both neutral and adaptive genetic variation declined with distance from glacial refugia. Nevertheless, there were also some indications from differential isolation by distance and allele abundance patterns that weak effects of selection have been superimposed on the main drift effect in the PGE zone.

Conservation and divergence in the frog immunome: pyrosequencing and de novo assembly of immune tissue transcriptomes

BACKGROUND

Frogs are a diverse group of vertebrates for which limited genomic resources are available. Natural frog populations face a multitude of threats, including habitat degradation, infectious disease, and environmental change. Characterizing the functional genomics of anuran tissues in general - and the immune system in particular - will enhance our knowledge of genetic and epigenetic responses to environmental threats and inform conservation and recovery efforts.

RESULTS

To increase the number of species with genomic datasets and characterize gene expression in immune-related tissues, we sequenced the transcriptomes of three tissues from two frogs (Espadarana prosoblepon and Lithobates yavapaiensis) on the Roche 454 GS FLX platform. Our sequencing produced 8881 E. prosoblepon and 5428 L. yavapaiensis annotated gene products after de novo assembly and Gene Ontology classification. Transcripts of the innate and acquired immune system were expressed in all three tissues. Inflammatory response and acquired immunity transcripts were significantly more diverged between E. prosoblepon and L. yavapaiensis compared to innate immunity and immune system development transcripts. Immune-related transcripts did not show an overall elevated rate of functional evolution, with the exception of glycosyl proteases, which include lysozymes, central bacterial and fungal-killing enzymes of the innate immune system.

CONCLUSIONS

The three frog transcriptomes provide more than 600 Mbp of new genomic data, and will serve as a valuable framework for future comparative studies of non-model anurans. Additionally, we show that immune gene divergence varies by functional group and that transcriptome studies can be useful in comparing rates of evolutionary change across gene families.

Genetic diversity of MHC class I loci in six non-model frogs is shaped by positive selection and gene duplication

Comparative studies of major histocompatibility complex (MHC) genes across vertebrate species can reveal the evolutionary processes that shape the structure and function of immune regulatory proteins. In this study, we characterized MHC class I sequences from six frog species representing three anuran families (Hylidae, Centrolenidae and Ranidae). Using cDNA from our focal species, we amplified a total of 79 unique sequences spanning exons 2-4 that encode the extracellular domains of the functional alpha chain protein. We compared intra- and interspecific nucleotide and amino-acid divergence, tested for recombination, and identified codon sites under selection by estimating the rate of non-synonymous to synonymous substitutions with multiple codon-based maximum likelihood methods. We determined that positive (diversifying) selection was acting on specific amino-acid sites located within the domains that bind pathogen-derived peptides. We also found significant signals of recombination across the physical distance of the genes. Finally, we determined that all the six species expressed two or three putative classical class I loci, in contrast to the single locus condition of Xenopus laevis. Our results suggest that MHC evolution in anurans is a dynamic process and that variation in numbers of loci and genetic diversity can exist among taxa. Thus, the accumulation of genetic data for more species will be useful in further characterizing the relative importance of processes such as selection, recombination and gene duplication in shaping MHC loci among amphibian lineages.

Interspecific hybridization increases MHC class II diversity in two sister species of newts

Our understanding of the evolutionary mechanisms generating variation within the highly polymorphic major histocompatibility complex (MHC) genes remains incomplete. Assessing MHC variation across multiple populations, of recent and ancient divergence, may facilitate understanding of geographical and temporal aspects of variation. Here, we applied 454 sequencing to perform a large-scale, comprehensive analysis of MHC class II in the closely related, hybridizing newts, Lissotriton vulgaris (Lv) and Lissotriton montandoni (Lm). Our study revealed an extensive (299 alleles) geographically structured polymorphism. Populations at the southern margin of the Lv distribution, inhabited by old and distinct lineages (southern Lv), exhibited moderate MHC variation and strong population structure, indicating little gene flow or extensive local adaptation. Lissotriton vulgaris in central Europe and the northern Balkans (northern Lv) and almost all Lm populations had a high MHC variation. A much higher proportion of MHC alleles was shared between Lm and northern Lv than between Lm and southern Lv. Strikingly, the average pairwise F(ST) between northern Lv and Lm was significantly lower than between northern and southern Lv for MHC, but not for microsatellites. Thus, high MHC variation in Lm and northern Lv may result from gene flow between species. We hypothesize that the interspecific exchange of MHC genes may be facilitated by frequency-dependent selection. A marginally significant correlation between the MHC and microsatellite allelic richness indicates that demographic factors may have contributed to the present-day pattern of MHC variation, but unequivocal signatures of adaptive evolution in MHC class II sequences emphasize the role of selection on a longer timescale.

MHC genotypes associate with resistance to a frog-killing fungus

The emerging amphibian disease chytridiomycosis is caused by the fungal pathogen Batrachochytrium dendrobatidis (Bd). Amphibian populations and species differ in susceptibility to Bd, yet we know surprisingly little about the genetic basis of this natural variation. MHC loci encode peptides that initiate acquired immunity in vertebrates, making them likely candidates for determining disease susceptibility. However, MHC genes have never been characterized in the context of chytridiomycosis. Here, we performed experimental Bd infections in laboratory-reared frogs collected from five populations that show natural variation in Bd susceptibility. We found that alleles of an expressed MHC class IIB locus associate with survival following Bd infection. Across populations, MHC heterozygosity was a significant predictor of survival. Within populations, MHC heterozygotes and individuals bearing MHC allele Q had a significantly reduced risk of death, and we detected a significant signal of positive selection along the evolutionary lineage leading to allele Q. Our findings demonstrate that immunogenetic variation affects chytridiomycosis survival under controlled experimental conditions, confirming that host genetic polymorphisms contribute to chytridiomycosis resistance.

Towards the simplification of MHC typing protocols: targeting classical MHC class II genes in a passerine, the pied flycatcher Ficedula hypoleuca

Background

Major Histocompatibility Complex (MHC) has drawn the attention of evolutionary biologists due to its importance in crucial biological processes, such as sexual selection and immune response in jawed vertebrates. However, the characterization of classical MHC genes subjected to the effects of natural selection still remains elusive in many vertebrate groups. Here, we have tested the suitability of flanking intron sequences to guide the selective exploration of classical MHC genes driving the co-evolutionary dynamics between pathogens and their passerine (Aves, Order Passeriformes) hosts.

Findings

Intronic sequences flanking the usually polymorphic exon 2 were isolated from different species using primers sitting on conserved coding regions of MHC class II genes (β chain). Taking the pied flycatcher Ficedula hypoleuca as an example, we demonstrate that careful primer design can evade non-classical MHC gene and pseudogene amplification. At least four polymorphic and expressed loci were co-replicated using a single pair of primers in five non-related individuals (N = 28 alleles). The cross-amplification and preliminary inspection of similar MHC fragments in eight unrelated songbird taxa suggests that similar approaches can also be applied to other species.

Conclusions

Intron sequences flanking the usually polymorphic exon 2 may assist the specific investigation of classical MHC class II B genes in species characterized by extensive gene duplication and pseudogenization. Importantly, the evasion of non-classical MHC genes with a more specific function and non-functional pseudogenes may accelerate data collection and diminish lab costs. Comprehensive knowledge of gene structure, polymorphism and expression profiles may be useful not only for the selective examination of evolutionarily relevant genes but also to restrict chimera formation by minimizing the number of co-amplifying loci.

Characterisation of MHC class II DRB genes in the northern tree shrew (Tupaia belangeri)

Genes of the major histocompatibility complex (MHC) mainly code for proteins of the immune system of jawed vertebrates. In particular, MHC class I and II cell surface proteins are crucial for the self/non-self discrimination of the adaptive immune system and are the most polymorphic genes in vertebrates. Positive selection, gene duplications and pseudogenes shape the face of the MHC and reflect a highly dynamic evolution. Here, we present for the first time data of the highly polymorphic MHC class II DRB exon 2 of a representative of the mammalian order scandentia, the northern tree shrew Tupaia belangeri. We found up to eight different alleles per individual and determined haplotype constitution by intensively studying their inheritance. The alleles were assigned to four putative loci, all of which were polymorphic. Only the most polymorphic locus was subject to positive selection within the antigen binding sites and only alleles of this locus were transcribed.

Species-specific evolution of class I MHC genes in iguanas (order: Squamata; subfamily: Iguaninae)

Over the last few decades, the major histocompatibility complex (MHC) has emerged as a model for understanding the influence of natural selection on genetic diversity in populations as well as for investigating the genetic basis of host resistance to pathogens. However, many vertebrate taxa remain underrepresented in the field of MHC research, preventing its application to studies of disease, evolution, and conservation genetics in these groups. This is particularly true for squamates, which are by far the most diversified order of non-avian reptiles but have not been the subject of any recent MHC studies. In this paper, we present MHC class I complementary DNA data from three squamate species in the subfamily Iguaninae (iguanas): the Galápagos marine iguana (Amblyrhynchus cristatus), the Galápagos land iguana (Conolophus subcristatus), and the green iguana (Iguana iguana). All sequences obtained are related to the few published class I genes from other squamates. There is evidence for multiple loci in each species, and the conserved alpha-3 domain appears to be evolving in a species-specific manner. Conversely, there is some indication of shared polymorphism between species in the peptide-binding alpha-1 and alpha-2 domains, suggesting that these two regions have different phylogenetic histories. The great similarity between alpha-3 sequences in marine iguanas in particular suggests that concerted evolution is acting to homogenize class I loci within species. However, while less likely, the data are also compatible with a birth and death model of evolution.