Trans-species polymorphism and evidence of selection on class II MHC loci in tuco-tucos (Rodentia: Ctenomyidae)

Evolutionary variation in gene conversion at the avian MHC is explained by fluctuating selection, gene copy numbers and life history

The major histocompatibility complex (MHC) multigene family encodes key pathogen-recognition molecules of the vertebrate adaptive immune system. Hyper-polymorphism of MHC genes is de novo generated by point mutations, but new haplotypes may also arise by re-shuffling of existing variation through intra- and inter-locus gene conversion. Although the occurrence of gene conversion at the MHC has been known for decades, we still have limited understanding of its functional importance. Here, I took advantage of extensive genetic resources (~9000 sequences) to investigate broad scale macroevolutionary patterns in gene conversion processes at the MHC across nearly 200 avian species. Gene conversion was found to constitute a universal mechanism in birds, as 83% of species showed footprints of gene conversion at either MHC class and 25% of all allelic variants were attributed to gene conversion. Gene conversion processes were stronger at MHC-II than MHC-I, but inter-specific variation at both MHC classes was explained by similar evolutionary scenarios, reflecting fluctuating selection towards different optima and drift. Gene conversion showed uneven phylogenetic distribution across birds and was driven by gene copy number variation, supporting significant role of inter-locus gene conversion processes in the evolution of the avian MHC. Finally, MHC gene conversion was stronger in species with fast life histories (high fecundity) and in long-distance migrants, likely reflecting variation in population sizes and host-pathogen coevolutionary dynamics. The results provide a robust comparative framework for understanding macroevolutionary variation in gene conversion at the avian MHC and reinforce important contribution of this mechanism to functional MHC diversity.

Population-wide genetic analysis of Ovar-DQA1 and DQA2 loci across sheep breeds in India revealed their evolutionary importance and fitness of sheep in a tropical climate

Genetic variability at the major histocompatibility complex (MHC) is important in any species due to significant role played by MHC for antigen presentation. DQA locus has not been studied for its genetic variability across sheep population in India. In the present study, MHC of sheep at DQA1 and DQA2 loci were evaluated across 17 Indian sheep breeds. Results revealed high degree of heterozygosity (10.34% to 100% for DQA1 and 37.39 to 100% for DQA2). 18 DQA1 alleles and 22 DQA2 alleles were isolated in different breeds. Nucleotide content for DQA region revealed richness of AT content (54.85% for DQA1 and 53.89% for DQA2). DQA1 and DQA2 sequences clustered independently. We could see evidence of divergence of DQA as DQA1 and DQA2 across sheep breeds. Wu-Kabat variability index revealed vast genetic variation across DQA1 and DQA2, specifically at peptide binding sites (PBS) that consisted 21 residues for DQA1 and 17 residues for DQA2. Evolutionary analysis revealed the presence of positive and balancing selection for DQA1 locus, however DQA2 was under purifying selection across sheep breeds. Higher heterozygosity and large diversity at both loci especially at PBS indicated the fitness of the sheep population for evading pathogens and adapt to the harsh tropical climate.

Cross-population genetic analysis revealed genetic variation and selection in the Ovar-DRB1 gene of Indian sheep breeds

In sheep, MHC variability is studied widely to explore disease association. The aim of the current study was to explore the genetic diversity of Ovar-DRB diversity across sheep breeds of India. Here, Ovar-DRB1 locus was studied across 20 sheep breeds. DRB1 was amplified (301 bp) and sequenced using a PCR-sequence-based typing approach. Results revealed a high degree of heterozygosity across breeds (mean: 73.99%). Overall mean distance for DRB1 was highest in Sangamneri (0.18) and lowest in Madgyal sheep (0.10). There was a higher rate of transition, across breeds. Further, 39 alleles were isolated in different breeds, out of which 10 were new. To allow easy access and use of the immune-polymorphic database, an online database management system was launched (http://www.mhcdbms.in/). Nucleotide content across breeds for the DRB1 region revealed the richness of GC content (59.26%). Wu-Kabat index revealed vast genetic variation across peptide binding sites (PBS) of DRB1. Residues 6, 66, 69, 52, and 81, were polymorphic showing utility for antigen presentation. All breeds were under positive selection for DRB1 locus (dN > dS). Study revealed the importance of DRB locus diversity for beta chain specifically at PBS across sheep breeds of the Indian subcontinent and presented evidence of positive selection for DRB owing to its evolutionary significance.

Appearances are deceptive: a cryptic lineage within the assumed distributional boundaries of Ctenomys talarum (Rodentia, Ctenomyidae)

The coastal dunes of the southeast of Buenos Aires province, Argentina, present two different described species of tuco-tucos: Ctenomys talarum and C. australis. Formerly, C. talarum was subdivided into three subspecies (C. t. talarum, C. t. recessus, and C. t. occidentalis), mainly based on its geographic distribution and phenotypic variation in characters of external morphology (e.g., body size and pelage color). This study assesses the phylogenetic relationships of C. talarum, focusing on the populations at the western end of its coastal distribution (localities of Pehuen-Có and Sauce Grande), which have been previously identified as highly genetically divergent. In this regard, populations distributed throughout the range of the species were sampled. Complete DNA sequences of the mitochondrial cytochrome b gene (1,140 bp), partial sequences of the mitochondrial D-loop region (426 bp), and partial sequences of the nuclear DNA intron 8 of the β-fibrinogen (about 870 bp) gene were used for the analyses. Phylogenetic inferences based on mitochondrial and nuclear markers were performed separately or combined to obtain a species tree. Populations distributed at the western end of the coastal dunes (between Pehuen-Có and Sauce Grande), previously assumed as C. talarum, were found to belong to an independent lineage relative to the other populations from the Pampas region. The average genetic distance between these two lineages is within the order of the genetic distances observed between different species of the genus. Also, our results show that this lineage of Ctenomys presents a high affinity with the magellanicus group, which is distributed further south, in Patagonia. In conclusion, tuco-tuco populations occurring in the coastal expanse between these two localities should be considered a possible distinct cryptic species, highly differentiated from C. talarum.

Adaptive evolution of major histocompatibility complex class I immune genes and disease associations in coastal juvenile sea turtles

Characterizing polymorphism at the major histocompatibility complex (MHC) genes is key to understanding the vertebrate immune response to disease. Despite being globally afflicted by the infectious tumour disease fibropapillomatosis (FP), immunogenetic variation in sea turtles is minimally explored. We sequenced the α 1 peptide-binding region of MHC class I genes (162 bp) from 268 juvenile green (Chelonia mydas) and 88 loggerhead (Caretta caretta) sea turtles in Florida, USA. We recovered extensive variation (116 alleles) and trans-species polymorphism. Supertyping analysis uncovered three functional MHC supertypes corresponding to the three well-supported clades in the phylogeny. We found significant evidence of positive selection at seven amino acid sites in the class I exon. Random forest modelling and risk ratio analysis of Ch. mydas alleles uncovered one allele weakly associated with smooth FP tumour texture, which may be associated with disease outcome. Our study represents the first characterization of MHC class I diversity in Ch. mydas and the largest sample of sea turtles used to date in any study of adaptive genetic variation, revealing tremendous genetic variation and high adaptive potential to viral pathogen threats. The novel associations we identified between MHC diversity and FP outcomes in sea turtles further highlight the importance of evaluating genetic predictors of disease, including MHC and other functional markers.

Trans‐specific polymorphism and the convergent evolution of supertypes in major histocompatibility complex class II genes in darters ( Etheostoma )

Major Histocompatibility Complex (MHC) genes are one of the most polymorphic gene groups known in vertebrates. MHC genes also exhibit allelic variants that are shared among taxa, referred to as trans‐specific polymorphism (TSP). The role that selection plays in maintaining such high diversity within species, as well as TSP, is an ongoing discussion in biology. In this study, we used deep‐sequencing techniques to characterize MHC class IIb gene diversity in three sympatric species of darters. We found at least 5 copies of the MHC gene in darters, with 126 genetic variants encoding 122 unique amino acid sequences. We identified four supertypes based on the binding properties of proteins encoded by the sequences. Although each species had a unique pool of variants, many variants were shared between species pairs and across all three species. Phylogenetic analysis showed that the variants did not group together monophyletically based on species identity or on supertype. An expanded phylogenetic analysis showed that some darter alleles grouped together with alleles from other percid fishes. Our findings show that TSP occurs in darters, which suggests that balancing selection is acting at the genotype level. Supertypes, however, are most likely evolving convergently, as evidenced by the fact that alleles do not form monophyletic groups based on supertype. Our research demonstrates that selection may be acting differently on MHC genes at the genotype and supertype levels, selecting for the maintenance of high genotypic diversity while driving the convergent evolution of similar MHC phenotypes across different species.

Major histocompatibility complex class II DR and DQ evolution and variation in wild capuchin monkey species (Cebinae)

The major histocompatibility complex (MHC) is an important gene complex contributing to adaptive immunity. Studies of platyrrhine MHC have focused on identifying experimental models of immune system function in the equivalent Human Leukocyte Antigen (HLA). These genes have thus been explored primarily in captive platyrrhine individuals from research colonies. However, investigations of standing MHC variation and evolution in wild populations are essential to understanding its role in immunity, sociality and ecology. Capuchins are a promising model group exhibiting the greatest habitat diversity, widest diet breadth and arguably the most social complexity among platyrrhines, together likely resulting in varied immunological challenges. We use high-throughput sequencing to characterize polymorphism in four Class II DR and DQ exons for the first time in seven capuchin species. We find evidence for at least three copies for DQ genes and at least five for DRB, with possible additional unrecovered diversity. Our data also reveal common genotypes that are inherited across our most widely sampled population, Cebus imitator in Sector Santa Rosa, Costa Rica. Notably, phylogenetic analyses reveal that platyrrhine DQA sequences form a monophyletic group to the exclusion of all Catarrhini sequences examined. This result is inconsistent with the trans-species hypothesis for MHC evolution across infraorders in Primates and provides further evidence for the independent origin of current MHC genetic diversity in Platyrrhini. Identical allele sharing across cebid species, and more rarely genera, however, does underscore the complexity of MHC gene evolution and the need for more comprehensive assessments of allelic diversity and genome structure.

Distinct evolutionary trajectories of MHC class I and class II genes in Old World finches and buntings

Major histocompatibility complex (MHC) genes code for key proteins of the adaptive immune system, which present antigens from intra-cellular (MHC class I) and extra-cellular (MHC class II) pathogens. Because of their unprecedented diversity, MHC genes have long been an object of scientific interest, but due to methodological difficulties in genotyping of duplicated loci, our knowledge on the evolution of the MHC across different vertebrate lineages is still limited. Here, we compared the evolution of MHC class I and class II genes in three sister clades of common passerine birds, finches (Fringillinae and Carduelinae) and buntings (Emberizidae) using a uniform methodological (genotyping and data processing) approach and uniform sample sizes. Our analyses revealed contrasting evolutionary trajectories of the two MHC classes. We found a stronger signature of pervasive positive selection and higher allele diversity (allele numbers) at the MHC class I than class II. In contrast, MHC class II genes showed greater allele divergence (in terms of nucleotide diversity) and a much stronger recombination (gene conversion) signal. Gene copy numbers at both MHC class I and class II evolved via fluctuating selection and drift (Brownian Motion evolution), but the evolutionary rate was higher at class I. Our study constitutes one of few existing examples, where evolution of MHC class I and class II genes was directly compared using a multi-species approach. We recommend that re-focusing MHC research from single-species and single-class approaches towards multi-species analyses of both MHC classes can substantially increase our understanding MHC evolution in a broad phylogenetic context.

Spatiotemporal adaptive evolution of an MHC immune gene in a frog-fungus disease system

Genetic diversity of major histocompatibility complex (MHC) genes is linked to reduced pathogen susceptibility in amphibians, but few studies also examine broad spatial and temporal patterns of MHC and neutral genetic diversity. Here, we characterized range-wide MHC diversity in the Northern leopard frog, Rana pipiens, a species found throughout North America that is experiencing disease-related declines. We used previously sequenced neutral markers (mitochondrial DNA and microsatellites), sequenced an expressed MHC class IIß gene fragment, and measured infection prevalence and intensity of the global fungal pathogen Batrachochytrium dendrobatidis (Bd) across 14 populations. Four populations were sampled across two decades, enabling temporal comparisons of selection and demography. We recovered 37 unique MHC alleles, including 17 that were shared across populations. Phylogenetic and population genetic patterns between MHC and neutral markers were incongruent, and five MHC codon positions associated with peptide binding were under positive selection. MHC heterozygosity, but not neutral marker heterozygosity, was a significant factor explaining spatial patterns of Bd prevalence, whereas only environmental variables predicted Bd intensity. MHC allelic richness (AR) decreased significantly over time but microsatellite-based AR did not, highlighting a loss of functional immunogenetic diversity that may be associated with Bd selective pressures. MHC supertype 4 was significantly associated with an elevated risk of Bd infection, whereas one supertype 2 allele was associated with a nearly significant reduced risk of Bd. Taken together, these results provide evidence that positive selection contributes to MHC class IIß evolution in R. pipiens and suggest that functional MHC differences across populations may contribute to disease adaptation.

Genetic variation and selection in the major histocompatibility complex Class II gene in the Guizhou pony

The Guizhou pony (GZP) is an indigenous species of equid found in the mountains of the Guizhou province in southwest China. We selected four regions of the equine leukocyte antigen (ELA), including DQA, DRA, DQB, and DRB, and used them to assess the diversity of the major histocompatibility complex (MHC) class II gene using direct sequencing technology. DRA had the lowest d_N/d_S ratio (0.560) compared with the other three loci, indicating that DRA was conserved and could be conserved after undergoing selective processes. Nine DQA, five DQB, nine DRA, and seven DRB codons were under significant positive selection at the antigen binding sites (ABS), suggesting that the selected residues in ABS may play a significant role in the innate immune system of the GZP. Two GZP alleles were shared with Przewalski’s horse, and six older GZP haplotypes had a better relationship with other horse species by one or two mutational steps, indicating that the GZP may be a natural ancient variety of equid. The specific diversity of ABS and the numbers of unique haplotypes in the evolutionary process affords this species a better genetic fitness and ability to adapt to the native environment.

A global analysis of selection at the avian MHC

Recent advancements in sequencing technology have resulted in rapid progress in the study of the major histocompatibility complex (MHC) in non-model avian species. Here, we analyze a global dataset of avian MHC class I and class II sequences (ca. 11,000 sequences from over 250 species) to gain insight into the processes that govern macroevolution of MHC genes in birds. Analysis of substitution rates revealed striking differences in the patterns of diversifying selection between passerine and non-passerine birds. Non-passerines showed stronger selection at MHC class II, which is primarily involved in recognition of extracellular pathogens, while passerines showed stronger selection at MHC class I, which is involved in recognition of intracellular pathogens. Positions of positively selected amino-acid residues showed marked discrepancies with peptide-binding residues (PBRs) of human MHC molecules, suggesting that using a human classification of PBRs to assess selection patterns at the avian MHC may be unjustified. Finally, our analysis provided evidence that indel mutations can make a substantial contribution to adaptive variation at the avian MHC.

Selection on MHC in a Context of Historical Demographic Change in 2 Closely Distributed Species of Tuco-tucos (Ctenomys australis and C. talarum)

Selection necessarily acts within the same current and historical demographic framework as neutral evolutionary processes, and the outcome of the interplay between these forces may vary according to their relative strength. In this study, we compare the variation at a major histocompatibility complex (MHC) locus (DRB exon 2), typically subject to strong diversifying selection, and mitochondrial diversity (D-loop) across populations encompassing the entire distribution of 2 species of South American subterranean rodents: Ctenomys australis and C. talarum (tuco-tucos). Although these species are parapatric along most of their distribution, historically they have followed distinct demographic trajectories associated with sea level changes during the Quaternary. We surveyed 8 populations of C. australis and 15 of C. talarum, from which we analyzed 70 and 212 D-loop haplotypes and 91 and 346 DRB genotypes, respectively. Both species have gone through a recent demographic expansion; however, the signal of this process only encompasses the entire distribution of one of the species: C. australis. While balancing selection on MHC in C. talarum-enhanced DRB diversity at the local level compared to D-loop, although not promoting divergence among populations, in C. australis local diversifying selection may have driven higher population differentiation at DRB than at D-loop. Our findings reinforce the idea that the relative strength of selection acting on MHC genes varies spatially and temporally within and among species, even between species using the same macrohabitat and exposed to similar immune challenges.

Extensive shared polymorphism at non-MHC immune genes in recently diverged North American prairie grouse

Gene polymorphisms shared between recently diverged species are thought to be widespread and most commonly reflect introgression from hybridization or retention of ancestral polymorphism through incomplete lineage sorting. Shared genetic diversity resulting from incomplete lineage sorting is usually maintained for a relatively short period of time, but under strong balancing selection it may persist for millions of years beyond species divergence (balanced trans-species polymorphism), as in the case of the major histocompatibility complex (MHC) genes. However, balancing selection is much less likely to act on non-MHC immune genes. The aim of this study was to investigate the patterns of shared polymorphism and selection at non-MHC immune genes in five grouse species from Centrocercus and Tympanuchus genera. For this purpose, we genotyped five non-MHC immune genes that do not interact directly with pathogens, but are involved in signaling and regulate immune cell growth. In contrast to previous studies with MHC, we found no evidence for balancing selection or balanced trans-species polymorphism among the non-MHC immune genes. No haplotypes were shared between genera and in most cases more similar allelic variants sorted by genus. Between species within genera, however, we found extensive shared polymorphism, which was most likely attributable to introgression or incomplete lineage sorting following recent divergence and large ancestral effective population size (i.e., weak genetic drift). Our study suggests that North American prairie grouse may have attained relatively low degree of reciprocal monophyly at nuclear loci and reinforces the rarity of balancing selection in non-MHC immune genes.

A genome-wide scan for genes under balancing selection in Drosophila melanogaster

BACKGROUND

In the history of population genetics balancing selection has been considered as an important evolutionary force, yet until today little is known about its abundance and its effect on patterns of genetic diversity. Several well-known examples of balancing selection have been reported from humans, mice, plants, and parasites. However, only very few systematic studies have been carried out to detect genes under balancing selection. We performed a genome scan in Drosophila melanogaster to find signatures of balancing selection in a derived (European) and an ancestral (African) population. We screened a total of 34 genomes searching for regions of high genetic diversity and an excess of SNPs with intermediate frequency.

RESULTS

In total, we found 183 candidate genes: 141 in the European population and 45 in the African one, with only three genes shared between both populations. Most differences between both populations were observed on the X chromosome, though this might be partly due to false positives. Functionally, we find an overrepresentation of genes involved in neuronal development and circadian rhythm. Furthermore, some of the top genes we identified are involved in innate immunity.

CONCLUSION

Our results revealed evidence of genes under balancing selection in European and African populations. More candidate genes have been found in the European population. They are involved in several different functions.

Genomic structure and expression pattern of MHC IIα and IIβ genes reveal an unusual immune trait in lined seahorse Hippocampus erectus

The major histocompatibility complex (MHC) genes are crucial in the adaptive immune system, and the gene duplication of MHC in animals can generally result in immune flexibility. In this study, we found that the lined seahorse (Hippocampus erectus) has only one gene copy number (GCN) of MHC IIα and IIβ, which is different from that in other teleosts. Together with the lack of spleen and gut-associated lymphatic tissue (GALT), the seahorse may be referred to as having a partial but natural "immunodeficiency". Highly variable amino acid residues were found in the IIα and IIβ domains, especially in the α1 and β1 domains with 9.62% and 8.43% allelic variation, respectively. Site models revealed seven and ten positively selected positions in the α1 and β1 domains, respectively. Real-time PCR experiments showed high expression levels of the MHC II genes in intestine (In), gill (Gi) and trunk kidney (TK) and medium in muscle (Mu) and brood pouch (BP), and the expression levels were significantly up-regulated after bacterial infection. Specially, relative higher expression level of both MHC IIα and IIβ was found in Mu and BP when compared with other fish species, in which MHC II is expressed negligibly in Mu. These results indicate that apart from TK, Gi and In, MU and BP play an important role in the immune response against pathogens in the seahorse. In conclusion, high allelic variation and strong positive selection in PBR and relative higher expression in MU and BP are speculated to partly compensate for the immunodeficiency.

Contrasting patterns of selection between MHC I and II across populations of Humboldt and Magellanic penguins

The evolutionary and adaptive potential of populations or species facing an emerging infectious disease depends on their genetic diversity in genes, such as the major histocompatibility complex (MHC). In birds, MHC class I deals predominantly with intracellular infections (e.g., viruses) and MHC class II with extracellular infections (e.g., bacteria). Therefore, patterns of MHC I and II diversity may differ between species and across populations of species depending on the relative effect of local and global environmental selective pressures, genetic drift, and gene flow. We hypothesize that high gene flow among populations of Humboldt and Magellanic penguins limits local adaptation in MHC I and MHC II, and signatures of selection differ between markers, locations, and species. We evaluated the MHC I and II diversity using 454 next-generation sequencing of 100 Humboldt and 75 Magellanic penguins from seven different breeding colonies. Higher genetic diversity was observed in MHC I than MHC II for both species, explained by more than one MHC I loci identified. Large population sizes, high gene flow, and/or similar selection pressures maintain diversity but limit local adaptation in MHC I. A pattern of isolation by distance was observed for MHC II for Humboldt penguin suggesting local adaptation, mainly on the northernmost studied locality. Furthermore, trans-species alleles were found due to a recent speciation for the genus or convergent evolution. High MHC I and MHC II gene diversity described is extremely advantageous for the long-term survival of the species.

Contrasting evolutionary histories of MHC class I and class II loci in grouse--effects of selection and gene conversion

Genes of the major histocompatibility complex (MHC) encode receptor molecules that are responsible for recognition of intracellular and extracellular pathogens (class I and class II genes, respectively) in vertebrates. Given the different roles of class I and II MHC genes, one might expect the strength of selection to differ between these two classes. Different selective pressures may also promote different rates of gene conversion at each class. Despite these predictions, surprisingly few studies have looked at differences between class I and II genes in terms of both selection and gene conversion. Here, we investigated the molecular evolution of MHC class I and II genes in five closely related species of prairie grouse (Centrocercus and Tympanuchus) that possess one class I and two class II loci. We found striking differences in the strength of balancing selection acting on MHC class I versus class II genes. More than half of the putative antigen-binding sites (ABS) of class II were under positive or episodic diversifying selection, compared with only 10% at class I. We also found that gene conversion had a stronger role in shaping the evolution of MHC class II than class I. Overall, the combination of strong positive (balancing) selection and frequent gene conversion has maintained higher diversity of MHC class II than class I in prairie grouse. This is one of the first studies clearly demonstrating that macroevolutionary mechanisms can act differently on genes involved in the immune response against intracellular and extracellular pathogens.

Genetic Variation at Exon 2 of the MHC Class II DQB Locus in Blue Whale (Balaenoptera musculus) from the Gulf of California

The genes of the Major Histocompatibility Complex (MHC) play an important role in the vertebrate immune response and are among the most polymorphic genes known in vertebrates. In some marine mammals, MHC genes have been shown to be characterized by low levels of polymorphism compared to terrestrial taxa; this reduction in variation is often explained as a result of lower pathogen pressures in marine habitats. To determine if this same reduction in variation applies to the migratory population of blue whales (Balaenoptera musculus) that occurs in the Gulf of California, we genotyped a 172 bp fragment of exon 2 of the MHC Class II DQB locus for 80 members of this population. Twenty-two putatively functional DQB allotypes were identified, all of which were homologous with DQB sequences from other cetacean species. Up to 5 putative alleles per individual were identified, suggesting that gene duplication has occurred at this locus. Rates of non-synonymous to synonymous substitutions (ω) and maximum likelihood analyses of models of nucleotide variation provided potential evidence of ongoing positive selection at this exon. Phylogenetic analyses of DQB alleles from B. musculus and 16 other species of cetaceans revealed trans-specific conservation of MHC variants, suggesting that selection has acted on this locus over prolonged periods of time. Collectively our findings reveal that immunogenic variation in blue whales is comparable to that in terrestrial mammals, thereby providing no evidence that marine taxa are subject to reduced pathogen-induced selective pressures.

Genetic diversity comparison of the DQA gene in European rabbit (Oryctolagus cuniculus) populations

The European rabbit (Oryctolagus cuniculus) natural populations within the species native region, the Iberian Peninsula, are considered a reservoir of genetic diversity. Indeed, the Iberia was a Pleistocene refuge to the species and currently two subspecies are found in the peninsula (Oryctolagus cuniculus cuniculus and Oryctolagus cuniculus algirus). The genes of the major histocompatibility complex (MHC) have been substantially studied in wild populations due to their exceptional variability, believed to be pathogen driven. They play an important function as part of the adaptive immune system affecting the individual fitness and population viability. In this study, the MHC variability was assessed by analysing the exon 2 of the DQA gene in several European rabbit populations from Portugal, Spain and France and in domestic breeds. Twenty-eight DQA alleles were detected, among which 18 are described for the first time. The Iberian rabbit populations are well differentiated from the French population and domestic breeds. The Iberian populations retained the higher allelic diversity with the domestic breeds harbouring the lowest; in contrast, the DQA nucleotide diversity was higher in the French population. Signatures of positive selection were detected in four codons which are putative peptide-binding sites and have been previously detected in other mammals. The evolutionary relationships showed instances of trans-species polymorphism. Overall, our results suggest that the DQA in European rabbits is evolving under selection and genetic drift.

Disparate independent genetic events disrupt the secondary metabolism gene perA in certain symbiotic Epichloë species

Peramine is an insect-feeding deterrent produced by Epichloë species in symbiotic association with C3 grasses. The perA gene responsible for peramine synthesis encodes a two-module nonribosomal peptide synthetase. Alleles of perA are found in most Epichloë species; however, peramine is not produced by many perA-containing Epichloë isolates. The genetic basis of these peramine-negative chemotypes is often unknown. Using PCR and DNA sequencing, we analyzed the perA genes from 72 Epichloë isolates and identified causative mutations of perA null alleles. We found nonfunctional perA-ΔR* alleles, which contain a transposon-associated deletion of the perA region encoding the C-terminal reductase domain, are widespread within the Epichloë genus and represent a prevalent mutation found in nonhybrid species. Disparate phylogenies of adjacent A2 and T2 domains indicated that the deletion of the reductase domain (R*) likely occurred once and early in the evolution of the genus, and subsequently there have been several recombinations between those domains. A number of novel point, deletion, and insertion mutations responsible for abolishing peramine production in full-length perA alleles were also identified. The regions encoding the first and second adenylation domains (A1 and A2, respectively) were common sites for such mutations. Using this information, a method was developed to predict peramine chemotypes by combining PCR product size polymorphism analysis with sequencing of the perA adenylation domains.

Characterization and genotyping of the DRB1 gene of the major histocompatibility complex (MHC) in the Marmota monax, animal model of hepatitis B

The major histocompatibility complex (MHC)-containing genes are among the most polymorphic in vertebrates. MHC genes code for proteins that are critical in the immune system response. In this study, the polymorphism of the second exon of the MHC class II DRB gene was characterized in the Eastern woodchuck (Marmota monax). Woodchucks chronically infected with the woodchuck hepatitis virus (WHV) represent the best available animal model for the study of chronic hepatitis B infection in humans. In the genotyped animals we found fifteen alleles, which were expressed in two independent loci and that were named DRB1A and DRB1B in this work. The 15 alleles investigated showed an elevated divergence. A significant excess of non-synonymous substitutions was detected, which could indicate that a historical positive selection is acting in the woodchuck DRB1 genes. This hypothesis was confirmed in our study by the high variability in or near the antigen binding sites (ABS) and by the results obtained in sequence variability analyses. This analysis identified the presence of a microsatellite sequence that is located at the start of the second intron, which could further allow the development of a fast and cheap semiautomatic sequencing method.

Advantageous diversity maintained by balancing selection in humans

Most human polymorphisms are neutral or slightly deleterious, but some genetic variation is advantageous and maintained in populations by balancing selection. Considered a rarity and overlooked for years, balanced polymorphisms have recently received renewed attention with several lines of evidence showing their relevance in human evolution. From theoretical work on its role in adaptation to empirical studies that identify its targets, recent developments have showed that balancing selection is more prevalent than previously thought. Here we review these developments and discuss their implications in our understanding of the influence of balancing selection in human evolution. We also review existing evidence on the biological functions that benefit most from advantageous diversity, and the functional consequences of these variants. Overall, we argue that balancing selection must be considered an important selective force in human evolution.

Characterization of class I- and class II-like major histocompatibility complex loci in pedigrees of North Atlantic right whales

North Atlantic right whales have one of the lowest levels of genetic variation at minisatellite loci, microsatellite loci, and mitochondrial control region haplotypes among mammals. Here, adaptive variation at the peptide binding region of class I and class II DRB-like genes of the major histocompatibility complex was assessed. Amplification of a duplicated region in 222 individuals revealed at least 11 class II alleles. Six alleles were assigned to the locus Eugl-DRB1 and 5 alleles were assigned to the locus Eugl-DRB2 by assessing segregation patterns of alleles from 81 parent/offspring pedigrees. Pedigree analysis indicated that these alleles segregated into 12 distinct haplotypes. Genotyping a smaller subset of unrelated individuals (n = 5 and 10, respectively) using different primer sets revealed at least 2 class II pseudogenes (with ≥ 4 alleles) and at least 3 class I loci (with ≥ 6 alleles). Class II sequences were significantly different from neutrality at peptide binding sites suggesting loci may be under the influence of balancing selection. Trans-species sharing of alleles was apparent for class I and class II sequences. Characterization of class II loci represents the first step in determining the relationship between major histocompatibility complex variability and factors affecting health and reproduction in this species.

Nucleotide variability at its limit? Insights into the number and evolutionary dynamics of the sex-determining specificities of the honey bee Apis mellifera

Deciphering the evolutionary processes driving nucleotide variation in multiallelic genes is limited by the number of genetic systems in which such genes occur. The complementary sex determiner (csd) gene in the honey bee Apis mellifera is an informative example for studying allelic diversity and the underlying evolutionary forces in a well-described model of balancing selection. Acting as the primary signal of sex determination, diploid individuals heterozygous for csd develop into females, whereas csd homozygotes are diploid males that have zero fitness. Examining 77 of the functional heterozygous csd allele pairs, we established a combinatorical criteria that provide insights into the minimum number of amino acid differences among those pairs. Given a data set of 244 csd sequences, we show that the total number of csd alleles found in A. mellifera ranges from 53 (locally) to 87 (worldwide), which is much higher than was previously reported (20). Using a coupon-collector model, we extrapolate the presence of in total 116–145 csd alleles worldwide. The hypervariable region (HVR) is of particular importance in determining csd allele specificity, and we provide for this region evidence of high evolutionary rate for length differences exceeding those of microsatellites. The proportion of amino acids driven by positive selection and the rate of nonsynonymous substitutions in the HVR-flanking regions reach values close to 1 but differ with respect to the HVR length. Using a model of csd coalescence, we identified the high originating rate of csd specificities as a major evolutionary force, leading to an origin of a novel csd allele every 400,000 years. The csd polymorphism frequencies in natural populations indicate an excess of new mutations, whereas signs of ancestral transspecies polymorphism can still be detected. This study provides a comprehensive view of the enormous diversity and the evolutionary forces shaping a multiallelic gene.

Diversity and selection of MHC class IIb gene exon3 in Chinese alligator

Our study used MHC class IIb gene exon3 complete sequence as markers to investigate genetic variability, selection and population differentiation in Chinese alligator. In this study, 282 bp MHC IIb exon3 complete sequence was got, none of the sequences contained insertions/deletions or stop codons, suggesting that all sequences might come from functional molecules in the genome. The neighbor-joining (NJ) tree revealed that Xuangzhou and Changxing populations were genetically close related, while Wild population showed the most diverse from the other. Gene flow (Nm) was very higher than one, suggesting that inter-group gene flow may have occurred. Furthermore, the results showed that MHC IIb gene might be a good molecular marker, we think that this technology could be used for Chinese alligator breeding and releasing in future.

Evolution of the MHC-DQB exon 2 in marine and terrestrial mammals

On the basis of a general low polymorphism, several studies suggest that balancing selection in the class II major histocompatibility complex (MHC) is weaker in marine mammals as compared with terrestrial mammals. We investigated such differential selection among Cetacea, Artiodactyla, and Primates at exon 2 of MHC-DQB gene by contrasting indicators of molecular evolution such as occurrence of transpecific polymorphisms, patterns of phylogenetic branch lengths by codon position, rates of nonsynonymous and synonymous substitutions as well as accumulation of variable sites on the sampling of alleles. These indicators were compared between the DQB and the mitochondrial cytochrome b gene (cytb) as a reference of neutral expectations and differences between molecular clocks resulting from life history and historical demography. All indicators showed that the influence of balancing selection on the DQB is more variable and overall weaker for cetaceans. In our sampling, ziphiids, the sperm whale, monodontids and the finless porpoise formed a group with lower DQB polymorphism, while mysticetes exhibited a higher DQB variation similar to that of terrestrial mammals as well as higher occurrence of transpecific polymorphisms. Different dolphins appeared in the two groups. Larger variation of selection on the cetacean DQB could be related to greater stochasticity in their historical demography and thus, to a greater complexity of the general ecology and disease processes of these animals.

Adaptive molecular evolution of the Major Histocompatibility Complex genes, DRA and DQA, in the genus Equus

BACKGROUND

Major Histocompatibility Complex (MHC) genes are central to vertebrate immune response and are believed to be under balancing selection by pathogens. This hypothesis has been supported by observations of extremely high polymorphism, elevated nonsynonymous to synonymous base pair substitution rates and trans-species polymorphisms at these loci. In equids, the organization and variability of this gene family has been described, however the full extent of diversity and selection is unknown. As selection is not expected to act uniformly on a functional gene, maximum likelihood codon-based models of selection that allow heterogeneity in selection across codon positions can be valuable for examining MHC gene evolution and the molecular basis for species adaptations.

RESULTS

We investigated the evolution of two class II MHC genes of the Equine Lymphocyte Antigen (ELA), DRA and DQA, in the genus Equus with the addition of novel alleles identified in plains zebra (E. quagga, formerly E. burchelli). We found that both genes exhibited a high degree of polymorphism and inter-specific sharing of allele lineages. To our knowledge, DRA allelic diversity was discovered to be higher than has ever been observed in vertebrates. Evidence was also found to support a duplication of the DQA locus. Selection analyses, evaluated in terms of relative rates of nonsynonymous to synonymous mutations (dN/dS) averaged over the gene region, indicated that the majority of codon sites were conserved and under purifying selection (dN <dS). However, the most likely evolutionary codon models allowed for variable rates of selection across codon sites at both loci and, at the DQA, supported the hypothesis of positive selection acting on specific sites.

CONCLUSIONS

Observations of elevated genetic diversity and trans-species polymorphisms supported the conclusion that balancing selection may be acting on these loci. Furthermore, at the DQA, positive selection was occurring at antigen binding sites, suggesting that a few selected residues may play a significant role in equid immune function. Future studies in natural equid populations will be valuable for understanding the functional significance of the uniquely diverse DRA locus and for elucidating the mechanism maintaining diversity at these MHC loci.

MHC variation, multiple simultaneous infections and physiological condition in the subterranean rodent Ctenomys talarum

Parasites and pathogens can play a significant role in shaping the genetic diversity of host populations, particularly at genes associated with host immune response. To explore this relationship in a natural population of vertebrates, we characterized Major Histocompatibility Complex (MHC) variation in the subterranean rodent Ctenomys talarum (the talas tuco-tuco) as a function of parasite load and ability to mount an adaptive immune response against a novel antigen. Specifically, we quantified genotypic diversity at the MHC class II DRB locus in relation to (1) natural variation in infection by multiple genera of parasites (potential agents of selection on MHC genes) and (2) antibody production in response to injection with sheep red blood cells (a measure of immunocompetence). Data were analyzed using co-inertia multivariate statistics, with epidemiological proxies for individual condition (hematocrit, leukocyte profile, body weight) and risk of parasite exposure (season of capture, sex). A significant excess of DRB heterozygotes was evident in the study population. Co-inertia analyses revealed significant associations between specific DRB alleles and both parasite load and intensity of humoral immune response against sheep red blood cells. The presence of specific DRB aminoacid sequences appeared to be more strongly associated with parasite load and response to a novel antigen than was heterozygosity at the DRB locus. These data suggest a role for parasite-driven balancing selection in maintaining MHC variation in natural populations of C. talarum. At the same time, these findings underscore the importance of using diverse parameters to study interactions among physiological conditions, immunocompetence, and MHC diversity in free-living animals that are confronted with multiple simultaneous immune challenges.

Allelic diversity of the MHC class II DRB genes in brown bears (Ursus arctos) and a comparison of DRB sequences within the family Ursidae

The allelic diversity of the DRB locus in major histocompatibility complex (MHC) genes was analyzed in the brown bear (Ursus arctos) from the Hokkaido Island of Japan, Siberia, and Kodiak of Alaska. Nineteen alleles of the DRB exon 2 were identified from a total of 38 individuals of U. arctos and were highly polymorphic. Comparisons of non-synonymous and synonymous substitutions in the antigen-binding sites of deduced amino acid sequences indicated evidence for balancing selection on the bear DRB locus. The phylogenetic analysis of the DRB alleles among three genera (Ursus, Tremarctos, and Ailuropoda) in the family Ursidae revealed that DRB allelic lineages were not separated according to species. This strongly shows trans-species persistence of DRB alleles within the Ursidae.

Genetic diversity and differentiation at MHC genes in island populations of tuatara (Sphenodon spp.)

Neutral genetic markers are commonly used to understand the effects of fragmentation and population bottlenecks on genetic variation in threatened species. Although neutral markers are useful for inferring population history, the analysis of functional genes is required to determine the significance of any observed geographical differences in variation. The genes of the major histocompatibility complex (MHC) are well-known examples of genes of adaptive significance and are particularly relevant to conservation because of their role in pathogen resistance. In this study, we survey diversity at MHC class I loci across a range of tuatara populations. We compare the levels of MHC variation with that observed at neutral microsatellite markers to determine the relative roles of balancing selection, diversifying selection and genetic drift in shaping patterns of MHC variation in isolated populations. In general, levels of MHC variation within tuatara populations are concordant with microsatellite variation. Tuatara populations are highly differentiated at MHC genes, particularly between the northern and Cook Strait regions, and a trend towards diversifying selection across populations was observed. However, overall our results indicate that population bottlenecks and isolation have a larger influence on patterns of MHC variation in tuatara populations than selection.

Sequence polymorphism and geographical variation at a positively selected MHC-DRB gene in the finless porpoise (Neophocaena phocaenoides): implication for recent differentiation of the Yangtze finless porpoise?

Sequence polymorphism at the MHC class II DRB locus was investigated in three finless porpoise (Neophocaena phocaenoides) populations in Chinese waters. Intragenic recombination and strong positive selection were the main forces in generating sequence diversity in the DRB gene. MHC sequence diversity changed significantly along the study period. Significant decrease in heterozygosity and lost alleles have been detected in the Yangtze River population and South China Sea population since 1990. Furthermore, there is a trend of increasing population differentiation over time. Especially, the genetic differentiation between the Yangtze River population and the Yellow Sea population was very low prior to 1990 (F (ST) = 0.036, P = 0.009), but became very significant after 1990 (F (ST) = 0.134, P < 0.001), suggesting a recent augmentation of genetic differentiation between both populations probably in a relatively short-term period. Porpoises from the Yangtze River displayed divergent frequencies of shared and private alleles from those displayed by two marine populations, which suggest that the former riverine population has been under a different selection regime (characteristic of a fresh water environment) than that of its marine counterparts.

Sequence polymorphism and evolution of three cetacean MHC genes

Sequence variability at three major histocompatibility complex (MHC) genes (DQB, DRA, and MHC-I) of cetaceans was investigated in order to get an overall understanding of cetacean MHC evolution. Little sequence variation was detected at the DRA locus, while extensive and considerable variability were found at the MHC-I and DQB loci. Phylogenetic reconstruction and sequence comparison revealed extensive sharing of identical MHC alleles among different species at the three MHC loci examined. Comparisons of phylogenetic trees for these MHC loci with the trees reconstructed only based on non-PBR sites revealed that allelic similarity/identity possibly reflected common ancestry and were not due to adaptive convergence. At the same time, trans-species evolution was also evidenced that the allelic diversity of the three MHC loci clearly pre-dated species divergence events according to the relaxed molecular clock. It may be the forces of balancing selection acting to maintain the high sequence variability and identical alleles in trans-specific manner at the MHC-I and DQB loci.

Polymorphism and selection in the major histocompatibility complex DRA and DQA genes in the family Equidae

The major histocompatibility complex genes coding for antigen binding and presenting molecules are the most polymorphic genes in the vertebrate genome. We studied the DRA and DQA gene polymorphism of the family Equidae. In addition to 11 previously reported DRA and 24 DQA alleles, six new DRA sequences and 13 new DQA alleles were identified in the genus Equus. Phylogenetic analysis of both DRA and DQA sequences provided evidence for trans-species polymorphism in the family Equidae. The phylogenetic trees differed from species relationships defined by standard taxonomy of Equidae and from trees based on mitochondrial or neutral gene sequence data. Analysis of selection showed differences between the less variable DRA and more variable DQA genes. DRA alleles were more often shared by more species. The DQA sequences analysed showed strong amongst-species positive selection; the selected amino acid positions mostly corresponded to selected positions in rodent and human DQA genes.

Trans-species polymorphism of the Mhc class II DRB-like gene in banded penguins (genus Spheniscus)

The Major Histocompatibility Complex (Mhc) class II DRB locus of vertebrates is highly polymorphic and some alleles may be shared between closely related species as a result of balancing selection in association with resistance to parasites. In this study, we developed a new set of PCR primers to amplify, clone, and sequence overlapping portions of the Mhc class II DRB-like gene from the 5'UTR end to intron 3, including exons 1, 2, and 3 and introns 1 and 2 in four species (20 Humboldt, six African, five Magellanic, and three Galapagos penguins) of penguin from the genus Spheniscus (Sphe). Analysis of gene sequence variation by the neighbor-joining method of 21 Sphe sequences and 20 previously published sequences from four other penguin species revealed overlapping clades within the Sphe species, but species-specific clades for the other penguin species. The overlap of the DRB-like gene sequence variants between the four Sphe species suggests that, despite their allopatric distribution, the Sphe species are closely related and that some shared DRB1 alleles may have undergone a trans-species inheritance because of balancing selection and/or recent rapid speciation. The new primers and PCR assays that we have developed for the identification of the DRB1 DNA and protein sequence variations appear to be useful for the characterization of the molecular evolution of the gene in closely related Penguin species and might be helpful for the assessment of the genetic health and the management of the conservation and captivity of these endangered species.

Evolutionary genetics: changed sex determination in honeybees

In several insects and fish, and probably some mammals, the gene controlling the male-female switch has changed during evolution. It now seems that this has also happened in honeybees, where the sex-determining gene has now been shown to be a duplicate of another Hymenopteran sex-determining gene.