Mhc class II diversity and balancing selection in greater prairie-chickens

Selection and demography drive range-wide patterns of MHC-DRB variation in mule deer

BACKGROUND

Standing genetic variation is important especially in immune response-related genes because of threats to wild populations like the emergence of novel pathogens. Genetic variation at the major histocompatibility complex (MHC), which is crucial in activating the adaptive immune response, is influenced by both natural selection and historical population demography, and their relative roles can be difficult to disentangle. To provide insight into the influences of natural selection and demography on MHC evolution in large populations, we analyzed geographic patterns of variation at the MHC class II DRB exon 2 locus in mule deer (Odocoileus hemionus) using sequence data collected across their entire broad range.

RESULTS

We identified 31 new MHC-DRB alleles which were phylogenetically similar to other cervid MHC alleles, and one allele that was shared with white-tailed deer (Odocoileus virginianus). We found evidence for selection on the MHC including high dN/dS ratios, positive neutrality tests, deviations from Hardy-Weinberg Equilibrium (HWE) and a stronger pattern of isolation-by-distance (IBD) than expected under neutrality. Historical demography also shaped variation at the MHC, as indicated by similar spatial patterns of variation between MHC and microsatellite loci and a lack of association between genetic variation at either locus type and environmental variables.

CONCLUSIONS

Our results show that both natural selection and historical demography are important drivers in the evolution of the MHC in mule deer and work together to shape functional variation and the evolution of the adaptive immune response in large, well-connected populations.

Long-Read Genome Assemblies Reveal Extraordinary Variation in the Number and Structure of MHC Loci in Birds

Our knowledge of the Major Histocompatibility Complex (MHC) in birds is limited because it often consists of numerous duplicated genes within individuals that are difficult to assemble with short read sequencing technology. Long-read sequencing provides an opportunity to overcome this limitation because it allows the assembly of long regions with repetitive elements. In this study, we used genomes based on long-read sequencing to predict the number and location of MHC loci in a broad range of bird taxa. From the long-read-based genomes of 34 species, we found that there was extremely large variation in the number of MHC loci between species. Overall, there were greater numbers of both class I and II loci in passerines than nonpasserines. The highest numbers of loci (up to 193 class II loci) were found in manakins (Pipridae), which had previously not been studied at the MHC. Our results provide the first direct evidence from passerine genomes of this high level of duplication. We also found different duplication patterns between species. In some species, both MHC class I and II genes were duplicated together, whereas in most species they were duplicated independently. Our study shows that the analysis of long-read-based genomes can dramatically improve our knowledge of MHC structure, although further improvements in chromosome level assembly are needed to understand the evolutionary mechanisms producing the extraordinary interspecific variation in the architecture of the MHC region.

Historical gene flow and profound spatial genetic structure among golden pheasant populations suggested by multi-locus analysis

Major histocompatibility complex (MHC) is a good marker system for geographical genetics since they are functional genes in the immune system that are likely to affect the fitness of the individual, and the survival and evolutionary potential of a population in a changing environment. Golden pheasant (Chrysolophus pictus) is a wild Phasianidae distributed in central and north China. In this study, we used a locus-specific genotyping technique for MHC IIB genes of golden pheasant. Combining with microsatellites (simple sequence repeat, SSR) and mitochondrial DNA (mtDNA) D-loop region, we investigated the demographic history and illuminate genetic structure of this bird in detail. SYR (south of Yangtze river) - NYR (north of Yangtze river) lineages, separated by Yangtze River, were defined in genetic structure of MHC IIB. NYR was supposed as refuge during glacial period, suggested by diversity parameters and more ancient alleles in this region. Based on this hypothesis, there was gene flow from NYR to SYR, which was proved by three pieces of evidence: (1) distinct demographic histories of SYR (kept stable) and NYR (experienced expansion); (2) specific affiliation of LC in genetic structure of SSR and MHC genes; (3) significant gene flow from NYR to SYR. Moreover, we also found balancing selection by combination of three Grouping A2's regions (SC, QL and North) into one in Grouping B4 (NYR) and no pattern of isolation by distance (IBD) found in MHC IIB, whereas for SSR we found a relatively strong and significant IBD. Several mechanisms in the evolution of MHC IIB genes, including recombination, historically positive selection, trans-species evolution and concerted evolution, were shown by molecular and phylogenetic analysis. Overall these results suggest the Yangtze River was inferred to be a geological barrier for this avian and NYR might experience population expansion, which invaded into a neighboring region. This study contributes to the understanding of the effects of geographic features on contemporary patterns of genetic variation in the golden pheasant in China, and helps us to define the adaptive unite (AU) for this avian.

Balancing selection and recombination as evolutionary forces caused population genetic variations in golden pheasant MHC class I genes

BACKGROUND

The major histocompatibility complex (MHC) genes are vital partners in the acquired immune processes of vertebrates. MHC diversity may be directly associated with population resistance to infectious pathogens. Here, we screened for polymorphisms in exons 2 and 3 of the IA1 and IA2 genes in 12 golden pheasant populations across the Chinese mainland to characterize their genetic variation levels, to understand the effects of historical positive selection and recombination in shaping class I diversity, and to investigate the genetic structure of wild golden pheasant populations.

RESULTS

Among 339 individual pheasants, we identified 14 IA1 alleles in exon 2 (IA1-E2), 11 IA1-E3 alleles, 27 IA2-E2 alleles, and 28 IA2-E3 alleles. The non-synonymous substitution rate was significantly greater than the synonymous substitution rate at sequences in the IA2 gene encoding putative peptide-binding sites but not in the IA1 gene; we also found more positively selected sites in IA2 than in IA1. Frequent recombination events resulted in at least 9 recombinant IA2 alleles, in accordance with the intermingling pattern of the phylogenetic tree. Although some IA alleles are widely shared among studied populations, large variation occurs in the number of IA alleles across these populations. Allele frequency analysis across 2 IA loci showed low levels of genetic differentiation among populations on small geographic scales; however, significant genetic differentiation was observed between pheasants from the northern and southern regions of the Yangtze River. Both STRUCTURE analysis and F-statistic (F ST ) value comparison classified those populations into 2 major groups: the northern region of the Yangtze River (NYR) and the southern region of the Yangtze River (SYR).

CONCLUSIONS

More extensive polymorphisms in IA2 than IA1 indicate that IA2 has undergone much stronger positive-selection pressure during evolution. Moreover, the recombination events detected between the genes and the intermingled phylogenetic pattern indicate that interlocus recombination accounts for much of the allelic variation in IA2. Analysis of the population differentiation implied that homogenous balancing selection plays an important part in maintaining an even distribution of MHC variations. The natural barrier of the Yangtze River and heterogeneous balancing selection might help shape the NYR-SYR genetic structure in golden pheasants.

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.

MyD88 Polymorphisms and Association with Susceptibility to Salmonella Pullorum

Myeloid differentiation primary response gene 88 (MYD88), a universal adapter protein, plays an important role in activating the nuclear factor-κB (NF-κB) and regulating the expression of proinflammatory genes like tumor necrosis factor (TNF) and interleukin-1 (IL-1), which were highly involved in Salmonella Pullorum infection. To detect the relationship between polymorphisms of the MyD88 gene and Salmonella Pullorum disease, we screened the coding region (CDS) of the MYD88 gene by DNA pool construction and sequencing based on case-control study. Eight single nucleotide polymorphisms (SNPs) in the sequenced fragment (5 exons), 7 known loci and one novel mutation named G4810372T (SNP8), were found in the fifth exon. In addition, we found 7 nonsynonymous substitutions. The allele frequency of only one SNP, g.4810191C > T (SNP1), was significantly different (P < 0.05) between case and control groups. The genotype frequencies of SNP1 (g.4810191C > T) and SNP3 (g.4810257G > T) were of significant difference between the case and the control groups (P < 0.05). Collectively, SNPs of the MyD88 gene were significantly associated with susceptibility to Salmonella Pullorum infection, which can be used as a disease-resistant marker in chicken. These results provided a theoretical basis for future research on chicken breeding by marker-assisted selection.

Contrasting patterns of selection and drift between two categories of immune genes in prairie-chickens

Immune-receptor genes of the adaptive immune system, such as the major histocompatibility complex (MHC), are involved in recognizing specific pathogens and are known to have high rates of adaptive evolution, presumably as a consequence of rapid co-evolution between hosts and pathogens. In contrast, many 'mediating' genes of the immune system do not interact directly with specific pathogens and are involved in signalling (e.g. cytokines) or controlling immune cell growth. As a consequence, we might expect stronger selection at immune-receptor than mediating genes, but these two types of genes have not been compared directly in wild populations. Here, we tested the hypothesis that selection differs between MHC (class I and II) and mediating genes by comparing levels of population differentiation across the range of greater prairie-chickens (Tympanuchus cupido). As predicted, there was stronger population differentiation and isolation by distance at immune receptor (MHC) than at either mediating genes or neutral microsatellites, suggesting a stronger role of local adaptation at the MHC. In contrast, mediating genes displayed weaker differentiation between populations than neutral microsatellites, consistent with selection favouring similar alleles across populations for mediating genes. In addition to selection, drift also had a stronger effect on immune receptor (MHC) than mediating genes as indicated by the stronger decline of MHC variation in relation to population size. This is the first study in the wild to show that the effects of selection and drift on immune genes vary across populations depending on their functional role.

High major histocompatibility complex class I polymorphism despite bottlenecks in wild and domesticated populations of the zebra finch (Taeniopygia guttata)

BACKGROUND

Two subspecies of zebra finch, Taeniopygia guttata castanotis and T. g. guttata are native to Australia and the Lesser Sunda Islands, respectively. The Australian subspecies has been domesticated and is now an important model system for research. Both the Lesser Sundan subspecies and domesticated Australian zebra finches have undergone population bottlenecks in their history, and previous analyses using neutral markers have reported reduced neutral genetic diversity in these populations. Here we characterize patterns of variation in the third exon of the highly variable major histocompatibility complex (MHC) class I α chain. As a benchmark for neutral divergence, we also report the first mitochondrial NADH dehydrogenase 2 (ND2) sequences in this important model system.

RESULTS

Despite natural and human-mediated population bottlenecks, we find that high MHC class I polymorphism persists across all populations. As expected, we find higher levels of nucleotide diversity in the MHC locus relative to neutral loci, and strong evidence of positive selection acting on important residues forming the peptide-binding region (PBR). Clear population differentiation of MHC allele frequencies is also evident, and this may be due to adaptation to new habitats and associated pathogens and/or genetic drift. Whereas the MHC Class I locus shows broad haplotype sharing across populations, ND2 is the first locus surveyed to date to show reciprocal monophyly of the two subspecies.

CONCLUSIONS

Despite genetic bottlenecks and genetic drift, all surveyed zebra finch populations have maintained high MHC Class I diversity. The diversity at the MHC Class I locus in the Lesser Sundan subspecies contrasts sharply with the lack of diversity in previously examined neutral loci, and may thus be a result of selection acting to maintain polymorphism. Given uncertainty in historical population demography, however, it is difficult to rule out neutral processes in maintaining the observed diversity. The surveyed populations also differ in MHC Class I allele frequencies, and future studies are needed to assess whether these changes result in functional immune differences.

Genetic variation of the major histocompatibility complex (MHC class II B gene) in the threatened Hume's pheasant, Syrmaticus humiae

Major histocompatibility complex (MHC) genes are the most polymorphic genes in vertebrates and encode molecules that play a crucial role in pathogen resistance. As a result of their diversity, they have received much attention in the fields of evolutionary and conservation biology. Here, we described the genetic variation of MHC class II B (MHCIIB) exon 2 in a wild population of Hume's pheasant (Syrmaticus humiae), which has suffered a dramatic decline in population over the last three decades across its ranges in the face of heavy exploitation and habitat loss. Twenty-four distinct alleles were found in 73 S. humiae specimens. We found seven shared alleles among four geographical groups as well as six rare MHCIIB alleles. Most individuals displayed between one to five alleles, suggesting that there are at least three MHCIIB loci of the Hume's pheasant. The dN ⁄ dS ratio at putative antigen-binding sites (ABS) was significantly greater than one, indicating balancing selection is acting on MHCIIB exon 2. Additionally, recombination and gene conversion contributed to generating MHCIIB diversity in the Hume's pheasant. One to three recombination events and seventy-five significant gene conversion events were observed within the Hume's pheasant MHCIIB loci. The phylogenetic tree and network analysis revealed that the Hume's pheasant alleles do not cluster together, but are scattered through the tree or network indicating a trans-species evolutionary mode. These findings revealed the evolution of the Hume's pheasant MHC after suffering extreme habitat fragmentation.

Sequence variations of the MHC class I gene exon 2 and exon 3 between infected and uninfected chickens challenged with Marek's disease virus

The major histocompatibility complex (MHC) among chickens has been well established as being associated with disease resistance and pathogens infection, but the genetic differences in MHC between chickens susceptible to certain infections and those chickens that remain uninfected have not been sufficiently determined. In this study, we sought the genetic basis that may underlie differences in susceptibility to infection among chickens by challenging four groups of broilers with Marek's disease virus (MDV). Over the course of the experiment, lesions began to appear between 21 and 35 days post challenge (dpc), and commercial broilers were not necessarily better than indigenous chickens in terms of disease resistance. The four groups showed neutral resistance to MDV infection validated by challenge results and evolutionary analysis of exons 2 and 3 of the MHC class I region. Several variable sites in exon 2 and exon 3 were exclusively appeared in infected chickens. Exon 3 was likely more crucial than exon 2 in disease resistance. Our observations offered a support for a potential association between promiscuous pathogens and conspicuous genetic diversity in the MHC class I region.

MHC class IIB exon 2 polymorphism in the Grey partridge (Perdix perdix) is shaped by selection, recombination and gene conversion

Among bird species, the most studied major histocompatibility complex (MHC) is the chicken MHC. Although the number of studies on MHC in free-ranging species is increasing, the knowledge on MHC variation in species closely related to chicken is required to understand the peculiarities of bird MHC evolution. Here we describe the variation of MHC class IIB (MHCIIB) exon 2 in a population of the Grey partridge (Perdix perdix), a species of high conservation concern throughout Europe and an emerging galliform model in studies of sexual selection. We found 12 alleles in 108 individuals, but in comparison to other birds surprisingly many sites show signatures of historical positive selection. Individuals displayed between two to four alleles both on genomic and complementary DNA, suggesting the presence of two functional MHCIIB loci. Recombination and gene conversion appear to be involved in generating MHCIIB diversity in the Grey partridge; two recombination breakpoints and several gene conversion events were detected. In phylogenetic analysis of galliform MHCIIB, the Grey partridge alleles do not cluster together, but are scattered through the tree instead. Thus, our results indicate that the Grey partridge MHCIIB is comparable to most other galliforms in terms of copy number and population polymorphism.

MHC class I variation in a natural blue tit population (Cyanistes caeruleus)

The major histocompatibility complex (MHC) is central to the vertebrate immune system and its highly polymorphic genes are considered to influence several life-history traits of individuals. To characterize the MHC in a natural population of blue tits (Cyanistes caeruleus) we investigated the class I exon 3 diversity of more than 900 individuals. We designed two pairs of motif-specific primers that reliably amplify independent subsets of MHC alleles. Applying denaturing gradient gel electrophoresis (DGGE) we obtained 48 independently inherited units of unique band patterns (DGGE-haplogroups), which were validated in a segregation analysis within 105 families. In a second approach, we extensively sequenced 6 unrelated individuals to confirm that DGGE-haplogroup composition reflects individual allelic variation. The highest number of different DGGE-haplogroups in a single individual corresponded in 19 MHC exon 3 sequences, suggesting a minimum of 10 amplified MHC class I loci in the blue tit. In total, we identified 50 unique functional and 3 non-functional sequences. Functional sequences showed high levels of recombination and strong positive selection in the antigen binding region, whereas nucleotide diversity was comparatively low in the range of all passerine species. Finally, in a phylogenetic comparison of passerine MHC class I exon 3 sequences we discuss conflicting evolutionary signals possibly due to recent gene duplication, recombination events and concerted evolution. Our results indicate that the described method is suitable to effectively explore the MHC diversity and its ecological impacts in blue tits in future studies.

Signatures of balancing selection are maintained at disease resistance loci following mating system evolution and a population bottleneck in the genus Capsella

Background

Population bottlenecks can lead to a loss of variation at disease resistance loci, which could have important consequences for the ability of populations to adapt to pathogen pressure. Alternatively, current or past balancing selection could maintain high diversity, creating a strong heterogeneity in the retention of polymorphism across the genome of bottlenecked populations. We sequenced part of the LRR region of 9 NBS-LRR disease resistance genes in the outcrossing Capsella grandiflora and the recently derived, bottlenecked selfing species Capsella rubella, and compared levels and patterns of nucleotide diversity and divergence with genome-wide reference loci.

Results

In strong contrast with reference loci, average diversity at resistance loci was comparable between C. rubella and C. grandiflora, primarily due to two loci with highly elevated diversity indicative of past or present balancing selection. Average between-species differentiation was also reduced at the set of R-genes compared with reference loci, which is consistent with the maintenance of ancestral polymorphism.

Conclusions

Historical or ongoing balancing selection on plant disease resistance genes is a likely contributor to the retention of ancestral polymorphism in some regions of the bottlenecked Capella rubella genome.

Polymorphism and Balancing Selection of MHC Class II DAB Gene in 7 Selective Flounder (Paralichthys olivaceus) Families

In order to determine the genetic variation of the MHC class IIB exon2 allele in the offspring, 700 fry from seven families of Japanese flounder challenged with V. anguillarum were studied, and different mortality rates were found in those families. Five to ten surviving and dead fry from each of the seven families were selected to study the MHC class II B exon2 gene with PCR and a direct sequencing method. One hundred and sixteen different exon2 sequences were found and 116 different alleles were identified, while a minimum of four loci were revealed in the MHC class II B exon2 gene. The ratio (d_N/d_S) of nonsynonymous substitution (d_N) to synonymous substitutions (d_S) in the peptide-binding region (PBR) of the MHC class IIB gene was 6.234, which indicated that balancing selection is acting on the MHC class IIB genes. The MHC IIB alleles were thus being passed on to their progeny. Some alleles were significantly more frequent in surviving than dead individuals. All together our data suggested that the alleles Paol-DAB*4301, Paol-DAB*4601, Paol-DAB*4302, Paol-DAB*3803, and Paol-DAB*4101 were associated with resistance to V. anguillarum in flounder.

Allelic polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of Cynoglossus semilaevis (Cynoglossidae)

Major histocompatibility complex (MHC) genes play an important role in the immune response of vertebrates. Allelic polymorphism and evolutionary mechanism of MHC genes have been investigated in many mammals, but much less is known in teleosts. We examined the polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of the half-smooth tongue sole (Cynoglossus semilaevis); 23 alleles were found in this species. Gene duplication manifested as three to six distinct sequences at each domain in the same individuals. Non-synonymous substitutions occurred at a significantly higher frequency than synonymous substitutions in the PBR domain, suggesting balancing selection for maintaining polymorphisms at the MHC II DAB locus. Many positive selection sites were found to act very intensely on antigen-binding sites of MHC class II DAB gene.

454 sequencing reveals extreme complexity of the class II Major Histocompatibility Complex in the collared flycatcher

Background

Because of their functional significance, the Major Histocompatibility Complex (MHC) class I and II genes have been the subject of continuous interest in the fields of ecology, evolution and conservation. In some vertebrate groups MHC consists of multiple loci with similar alleles; therefore, the multiple loci must be genotyped simultaneously. In such complex systems, understanding of the evolutionary patterns and their causes has been limited due to challenges posed by genotyping.

Results

Here we used 454 amplicon sequencing to characterize MHC class IIB exon 2 variation in the collared flycatcher, an important organism in evolutionary and immuno-ecological studies. On the basis of over 152,000 sequencing reads we identified 194 putative alleles in 237 individuals. We found an extreme complexity of the MHC class IIB in the collared flycatchers, with our estimates pointing to the presence of at least nine expressed loci and a large, though difficult to estimate precisely, number of pseudogene loci. Many similar alleles occurred in the pseudogenes indicating either a series of recent duplications or extensive concerted evolution. The expressed alleles showed unambiguous signals of historical selection and the occurrence of apparent interlocus exchange of alleles. Placing the collared flycatcher's MHC sequences in the context of passerine diversity revealed transspecific MHC class II evolution within the Muscicapidae family.

Conclusions

454 amplicon sequencing is an effective tool for advancing our understanding of the MHC class II structure and evolutionary patterns in Passeriformes. We found a highly dynamic pattern of evolution of MHC class IIB genes with strong signals of selection and pronounced sequence divergence in expressed genes, in contrast to the apparent sequence homogenization in pseudogenes. We show that next generation sequencing offers a universal, affordable method for the characterization and, in perspective, genotyping of MHC systems of virtually any complexity.