Recombinant DRB sequences produced by mismatch repair of heteroduplexes during cloning in Escherichia coli

New data from basal Australian songbird lineages show that complex structure of MHC class II β genes has early evolutionary origins within passerines

Background

The major histocompatibility complex (MHC) plays a crucial role in the adaptive immune system and has been extensively studied across vertebrate taxa. Although the function of MHC genes appears to be conserved across taxa, there is great variation in the number and organisation of these genes. Among avian species, for instance, there are notable differences in MHC structure between passerine and non-passerine lineages: passerines typically have a high number of highly polymorphic MHC paralogs whereas non-passerines have fewer loci and lower levels of polymorphism. Although the occurrence of highly polymorphic MHC paralogs in passerines is well documented, their evolutionary origins are relatively unexplored. The majority of studies have focussed on the more derived passerine lineages and there is very little empirical information on the diversity of the MHC in basal passerine lineages. We undertook a study of MHC diversity and evolutionary relationships across seven species from four families (Climacteridae, Maluridae, Pardalotidae, Meliphagidae) that comprise a prominent component of the basal passerine lineages. We aimed to determine if highly polymorphic MHC paralogs have an early evolutionary origin within passerines or are a more derived feature of the infraorder Passerida.

Results

We identified 177 alleles of the MHC class II β exon 2 in seven basal passerine species, with variation in numbers of alleles across individuals and species. Overall, we found evidence of multiple gene loci, pseudoalleles, trans-species polymorphism and high allelic diversity in these basal lineages. Phylogenetic reconstruction of avian lineages based on MHC class II β exon 2 sequences strongly supported the monophyletic grouping of basal and derived passerine species.

Conclusions

Our study provides evidence of a large number of highly polymorphic MHC paralogs in seven basal passerine species, with strong similarities to the MHC described in more derived passerine lineages rather than the simpler MHC in non-passerine lineages. These findings indicate an early evolutionary origin of highly polymorphic MHC paralogs in passerines and shed light on the evolutionary forces shaping the avian MHC.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0681-5) contains supplementary material, which is available to authorized users.

Characterization of class II β chain major histocompatibility complex genes in a family of Hawaiian honeycreepers: 'amakihi (Hemignathus virens)

Hawaiian honeycreepers (Drepanidinae) have evolved in the absence of mosquitoes for over five million years. Through human activity, mosquitoes were introduced to the Hawaiian archipelago less than 200 years ago. Mosquito-vectored diseases such as avian malaria caused by Plasmodium relictum and Avipoxviruses have greatly impacted these vulnerable species. Susceptibility to these diseases is variable among and within species. Due to their function in adaptive immunity, the role of major histocompatibility complex genes (Mhc) in disease susceptibility is under investigation. In this study, we evaluate gene organization and levels of diversity of Mhc class II β chain genes (exon 2) in a captive-reared family of Hawaii 'amakihi (Hemignathus virens). A total of 233 sequences (173 bp) were obtained by PCR+1 amplification and cloning, and 5720 sequences were generated by Roche 454 pyrosequencing. We report a total of 17 alleles originating from a minimum of 14 distinct loci. We detected three linkage groups that appear to represent three distinct haplotypes. Phylogenetic analysis revealed one variable cluster resembling classical Mhc sequences (DAB) and one highly conserved, low variability cluster resembling non-classical Mhc sequences (DBB). High net evolutionary divergence values between DAB and DBB resemble that seen between chicken BLB system and YLB system genes. High amino acid identity among non-classical alleles from 12 species of passerines (DBB) and four species of Galliformes (YLB) was found, suggesting that these non-classical passerine sequences may be related to the Galliforme YLB sequences.

Selective pressures on MHC class II genes in the guppy (Poecilia reticulata) as inferred by hierarchical analysis of population structure

Genes of the major histocompatibility complex, which are the most polymorphic of all vertebrate genes, are a pre-eminent system for the study of selective pressures that arise from host-pathogen interactions. Balancing selection capable of maintaining high polymorphism should lead to the homogenization of MHC allele frequencies among populations, but there is some evidence to suggest that diversifying selection also operates on the MHC. However, the pattern of population structure observed at MHC loci is likely to depend on the spatial and/or temporal scale examined. Here, we investigated selection acting on MHC genes at different geographic scales using Venezuelan guppy populations inhabiting four regions. We found a significant correlation between MHC and microsatellite allelic richness across populations, which suggests the role of genetic drift in shaping MHC diversity. However, compared to microsatellites, more MHC variation was explained by differences between populations within larger geographic regions and less by the differences between the regions. Furthermore, among proximate populations, variation in MHC allele frequencies was significantly higher compared to microsatellites, indicating that selection acting on MHC may increase population structure at small spatial scales. However, in populations that have significantly diverged at neutral markers, the population-genetic signature of diversifying selection may be eradicated in the long term by that of balancing selection, which acts to preserve rare alleles and thus maintain a common pool of MHC alleles.

Stepwise threshold clustering: a new method for genotyping MHC loci using next-generation sequencing technology

Genes of the vertebrate major histocompatibility complex (MHC) are of great interest to biologists because of their important role in immunity and disease, and their extremely high levels of genetic diversity. Next generation sequencing (NGS) technologies are quickly becoming the method of choice for high-throughput genotyping of multi-locus templates like MHC in non-model organisms. Previous approaches to genotyping MHC genes using NGS technologies suffer from two problems:1) a "gray zone" where low frequency alleles and high frequency artifacts can be difficult to disentangle and 2) a similar sequence problem, where very similar alleles can be difficult to distinguish as two distinct alleles. Here were present a new method for genotyping MHC loci--Stepwise Threshold Clustering (STC)--that addresses these problems by taking full advantage of the increase in sequence data provided by NGS technologies. Unlike previous approaches for genotyping MHC with NGS data that attempt to classify individual sequences as alleles or artifacts, STC uses a quasi-Dirichlet clustering algorithm to cluster similar sequences at increasing levels of sequence similarity. By applying frequency and similarity based criteria to clusters rather than individual sequences, STC is able to successfully identify clusters of sequences that correspond to individual or similar alleles present in the genomes of individual samples. Furthermore, STC does not require duplicate runs of all samples, increasing the number of samples that can be genotyped in a given project. We show how the STC method works using a single sample library. We then apply STC to 295 threespine stickleback (Gasterosteus aculeatus) samples from four populations and show that neighboring populations differ significantly in MHC allele pools. We show that STC is a reliable, accurate, efficient, and flexible method for genotyping MHC that will be of use to biologists interested in a variety of downstream applications.

Extensive variation at MHC DRB in the New Zealand sea lion (Phocarctos hookeri) provides evidence for balancing selection

Marine mammals are often reported to possess reduced variation of major histocompatibility complex (MHC) genes compared with their terrestrial counterparts. We evaluated diversity at two MHC class II B genes, DQB and DRB, in the New Zealand sea lion (Phocarctos hookeri, NZSL) a species that has suffered high mortality owing to bacterial epizootics, using Sanger sequencing and haplotype reconstruction, together with next-generation sequencing. Despite this species' prolonged history of small population size and highly restricted distribution, we demonstrate extensive diversity at MHC DRB with 26 alleles, whereas MHC DQB is dimorphic. We identify four DRB codons, predicted to be involved in antigen binding, that are evolving under adaptive evolution. Our data suggest diversity at DRB may be maintained by balancing selection, consistent with the role of this locus as an antigen-binding region and the species' recent history of mass mortality during a series of bacterial epizootics. Phylogenetic analyses of DQB and DRB sequences from pinnipeds and other carnivores revealed significant allelic diversity, but little phylogenetic depth or structure among pinniped alleles; thus, we could neither confirm nor refute the possibility of trans-species polymorphism in this group. The phylogenetic pattern observed however, suggests some significant evolutionary constraint on these loci in the recent past, with the pattern consistent with that expected following an epizootic event. These data may help further elucidate some of the genetic factors underlying the unusually high susceptibility to bacterial infection of the threatened NZSL, and help us to better understand the extent and pattern of MHC diversity in pinnipeds.

Genomic Porosity between Invasive Chondrostoma nasus and Endangered Endemic Parachondrostoma toxostoma (Cyprinidae): The Evolution of MHC IIB Genes

Two cyprinid species, Parachondrostoma toxostoma, an endemic threatened species, and Chondrostoma nasus, an invasive species, live in sympatry in southern France and form two sympatric zones where the presence of intergeneric hybrids is reported. To estimate the potential threat to endemic species linked to the introduction of invasive species, we focused on the DAB genes (functional MHC IIB genes) because of their adaptive significance and role in parasite resistance. More specifically, we investigated (1) the variability of MHC IIB genes, (2) the selection pattern shaping MHC polymorphism, and (3) the extent to which trans-species evolution and intergeneric hybridization affect MHC polymorphism. In sympatric areas, the native species has more diversified MHC IIB genes when compared to the invasive species, probably resulting from the different origins and dispersal of both species. A similar level of MHC polymorphism was found at population level in both species, suggesting similar mechanisms generating MHC diversity. In contrast, a higher number of DAB-like alleles per specimen were found in invasive species. Invasive species tended to express the alleles of two DAB lineages, whilst native species tended to express the alleles of only the DAB3 lineage. Hybrids have a pattern of MHC expression intermediate between both species. Whilst positive selection acting on peptide binding sites (PBS) was demonstrated in both species, a slightly higher number of positively selected sites were identified in C. nasus, which could result from parasite-mediated selection. Bayesian clustering analysis revealed a similar pattern of structuring for the genetic variation when using microsatellites or the MHC approach. We confirmed the importance of trans-species evolution for MHC polymorphism. In addition, we demonstrated bidirectional gene flow for MHC IIB genes in sympatric areas. The positive significant correlation between MHC and microsatellites suggests that demographic factors may contribute to MHC variation on a short time scale.

High diversity at PRDM9 in chimpanzees and bonobos

Background

The PRDM9 locus in mammals has increasingly attracted research attention due to its role in mediating chromosomal recombination and possible involvement in hybrid sterility and hence speciation processes. The aim of this study was to characterize sequence variation at the PRDM9 locus in a sample of our closest living relatives, the chimpanzees and bonobos.

Methodology/Principal Findings

PRDM9 contains a highly variable and repetitive zinc finger array. We amplified this domain using long-range PCR and determined the DNA sequences using conventional Sanger sequencing. From 17 chimpanzees representing three subspecies and five bonobos we obtained a total of 12 alleles differing at the nucleotide level. Based on a data set consisting of our data and recently published Pan PRDM9 sequences, we found that at the subspecies level, diversity levels did not differ among chimpanzee subspecies or between chimpanzee subspecies and bonobos. In contrast, the sample of chimpanzees harbors significantly more diversity at PRDM9 than samples of humans. Pan PRDM9 shows signs of rapid evolution including no alleles or ZnFs in common with humans as well as signals of positive selection in the residues responsible for DNA binding.

Conclusions and Significance

The high number of alleles specific to the genus Pan, signs of positive selection in the DNA binding residues, and reported lack of conservation of recombination hotspots between chimpanzees and humans suggest that PRDM9 could be active in hotspot recruitment in the genus Pan. Chimpanzees and bonobos are considered separate species and do not have overlapping ranges in the wild, making the presence of shared alleles at the amino acid level between the chimpanzee and bonobo species interesting in view of the hypothesis that PRDM9 plays a universal role in interspecific hybrid sterility.

Identification of MhcMafa-DRB alleles in a cohort of cynomolgus macaques of Vietnamese origin

Cynomolgus macaques have been used widely to build a research model of infectious and chronic diseases, as well as in transplantation studies, where disease susceptibility and/or resistance are associated with the major histocompatibility complex (MHC). To better elucidate polymorphisms and genetic differences in the Mafa-DRB locus, and facilitate the experimental use of cynomolgus macaques, we used pool screening combined with cloning and direct sequencing of polymerase chain reaction products to characterize MhcMafa-DRB gene alleles in 153 Vietnamese cynomolgus macaques. We identified 30 Mafa-DRB alleles belonging to 17 allelic lineages, including four novel sequences that had not been documented in earlier reports. The highest frequency allele was Mafa-DRB*W27:04, which was present in 7 of 35 (20%) monkeys. The next most frequent alleles were Mafa-DRB*3:07 and Mafa-DRB*W7:01, which were detected in 5 of 35 (14.3%) and 4 of 35 (11.4%) of the monkeys, respectively. The high-frequency alleles in this Vietnamese population may be high priority targets for additional characterization of immune functions. Only the DRB1*03 and DRB1*10 lineages were also present in humans, whereas the remaining alleles were monkey-specific lineages. We found 25 variable sites by aligning the deduced amino acid sequences of 29 identified alleles. Evolutionary and population analyses based on these sequences showed that human, rhesus, and cynomolgus macaques share several Mhc-DRB lineages and the shared polymorphisms in the DRB region may be attributable to the existence of interbreeding between rhesus and cynomolgus macaques. This information will promote the understanding of MHC diversity and polymorphism in cynomolgus macaques and increase the value of this species as a model for biomedical research.

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.

Evaluation of two approaches to genotyping major histocompatibility complex class I in a passerine-CE-SSCP and 454 pyrosequencing

Genes of the highly dynamic major histocompatibility complex (MHC) are directly linked to individual fitness and are of high interest in evolutionary ecology and conservation genetics. Gene duplication and positive selection usually lead to high levels of polymorphism in the MHC region, making genotyping of MHC a challenging task. Here, we compare the performance of two methods for MHC class I genotyping in a passerine with highly duplicated MHC class I genes: capillary electrophoresis-single-strand conformation polymorphism (CE-SSCP) analysis and 454 GS FLX Titanium pyrosequencing. According to our findings, the number of MHC variants (called alleles for simplicity) detected by CE-SSCP is significantly lower than detected by 454. To resolve discrepancies between the two methods, we cloned and Sanger sequenced a MHC class I amplicon for an individual with high number of alleles. We found a perfect congruence between cloning/Sanger sequencing results and 454. Thus, in case of multi-locus amplification, CE-SSCP considerably underestimates individual MHC diversity. However, numbers of alleles detected by both methods are significantly correlated, although the correlation is weak (r = 0.32). Thus, in systems with highly duplicated MHC, 454 provides more reliable information on individual diversity than CE-SSCP.

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.

Characterization and locus-specific typing of MHC class I genes in the red-billed gull (Larus scopulinus) provides evidence for major, minor, and nonclassical loci

A major challenge facing studies of major histocompatibility complex (MHC) evolution in birds is the difficulty in genotyping alleles at individual loci, and the consequent inability to investigate sequence variation and selection pressures for each gene. In this study, four MHC class I loci were isolated from the red-billed gull (Larus scopulinus), representing both the first characterized MHCI genes within Charadriiformes (shorebirds, gulls, and allies) and the first full-length MHCI sequences described outside Galloanserae (gamebirds + waterfowl). Complete multilocus genotypes were obtained for 470 individuals using a combination of reference-strand conformation analysis and direct sequencing of gene-specific amplification products, and variation of peptide-binding region (PBR) exons was surveyed for all loci. Each gene is transcribed and has conserved sequence features characteristic of antigen-presenting MHCI molecules. However, higher allelic variation, a more even allele frequency distribution, and evidence of positive selection acting on a larger number of PBR residues suggest that only one locus (Lasc-UAA) functions as a major classical MHCI gene. Lasc-UBA, with more limited variation and PBR motifs that encompass a subset of Lasc-UAA diversity, was assigned a putative minor classical function, whereas the divergent and largely invariant binding-groove motifs of Lasc-UCA and -UDA are suggestive of nonclassical loci with specialized ligand-binding roles.

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.

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

The major histocompatibility complex (Mhc) of domestic chickens has been characterized as small and relatively simple compared with that of mammals. However, there is growing evidence that the Mhc of many bird lineages may be more complex, even within the Order Galliformes. In this study, we measured genetic variation and balancing selection at Mhc loci in another galliform, the greater prairie-chicken. We cloned and sequenced a 239 bp fragment of Mhc Class II beta-chain (BLB) exon 2 in 14 individuals. There was a total of 10 unique sequences and a minimum of four BLB loci. The d(N)/d(S) ratio at peptide-binding codons was significantly greater than one, suggesting balancing selection is acting on the BLB. We also recovered two YLB sequences, which clustered tightly with YLB sequences from three other species: domestic chicken, black grouse and common quail. The relatively large number of loci revealed in our study suggests that even closely related galliforms differ in the level of Mhc variation and structure.

Targeted single nucleotide polymorphism (SNP) discovery in a highly polyploid plant species using 454 sequencing

Discovering single nucleotide polymorphisms (SNPs) in specific genes in a heterozygous polyploid plant species, such as sugarcane, is challenging because of the presence of a large number of homologues. To discover SNPs for mapping genes of interest, 454 sequencing of 307 polymerase chain reaction (PCR) amplicons (> 59 kb of sequence) was undertaken. One region of a four-gasket sequencing run, on a 454 Genome Sequencer FLX, was used for pooled PCR products amplified from each parent of a quantitative trait locus (QTL) mapping population (IJ76-514 x Q165). The sequencing yielded 96,755 (IJ76-514) and 86,241 (Q165) sequences with perfect matches to a PCR primer used in amplification, with an average sequence depth of approximately 300 and an average read length of 220 bases. Further analysis was carried out on amplicons whose sequences clustered into a single contig using an identity of 80% with the program cap3. In the more polymorphic sugarcane parent (Q165), 94% of amplicons (227/242) had evidence of a reliable SNP--an average of one every 35 bases. Significantly fewer SNPs were found in the pure Saccharum officinarum parent--with one SNP every 58 bases and SNPs in 86% (213/247) of amplicons. Using automatic SNP detection, 1632 SNPs were detected in Q165 sequences and 1013 in IJ76-514. From 225 candidate SNP sites tested, 209 (93%) were validated as polymorphic using the Sequenom MassARRAY system. Amplicon re-sequencing using the 454 system enables cost-effective SNP discovery that can be targeted to genes of interest and is able to perform in the highly challenging area of polyploid genomes.

Structure and evolution of a new avian MHC class II B gene in a sub-Antarctic seabird, the thin-billed prion (Procellariiformes: Pachyptila belcheri)

The major histocompatibility complex encodes molecules that present foreign peptides to T cells of the immune system. The peptide binding region (PBR) of these molecules is among the most polymorphic regions found in vertebrate taxa. Genomic cloning approaches are improving our understanding of the evolution of this multigene family in nonmodel avian groups. By building a cosmid library, a new MHC class II B gene, Pabe-DAB1, was isolated and characterized at the genomic level in a sub-Antarctic seabird, the thin-billed prion (Pachyptila belcheri). Pabe-DAB1 exhibits the hallmark structural features of functional MHC class II loci. Direct sequencing of the PBR encoding exon in a panel of prions revealed significantly higher levels of genetic diversity compared to two noncoding neutral loci, with most alleles differing by at least one replacement substitution in the peptide binding codons. We estimated evolutionary dynamics for Pabe-DAB1 using a variety of Bayesian and other approaches. Evidence for balancing selection comes from a spatially variable ratio of nonsynonymous-to-synonymous substitutions (mean d (N)/d (S) = 2.87) in the PBR, with sites predicted to be functionally relevant exhibiting the highest omega values. We estimate the population recombination rate to be approximately 0.3 per site per generation, indicating an important role for recombination in generating polymorphism at this locus. Pabe-DAB1 is among the few avian class II loci characterized at the genomic level and with a known intron-exon structure, a feature that greatly facilitated the amplification and sequencing of a single MHC locus in this species.

Balancing and directional selection at exon-2 of the MHC DQB1 locus among populations of odontocete cetaceans

The diversity of exon-2 (peptide-binding region) of the DQB1 locus (Class II, major histocompatibility complex, MHC) was investigated on an extended sample of populations of three focal cetacean species (two sibling delphinid species and another in the same family). We tested the hypothesis that dolphin populations with a worldwide distribution across different habitats and geographic regions will be under differential selective pressure by comparing DQB1 variation with variation at neutral markers and by investigating putative functional residues within the exon-2 sequence at the population level. Variation at the DQB1 locus was not correlated to neutral differentiation (assessed by comparison with microsatellite DNA markers), and overall F(ST) values were significantly lower for the MHC locus, consistent with expectations for balancing selection. Measures of heterozygosity and d(n)/d(s) ratios were also consistent with balancing selection. However, outliers in the F(ST) comparisons and the analysis of putative functional residues suggested incidences of directional selection in local populations.

Reducing cloning artifacts for recovery of allelic sequences by T7 endonuclease I cleavage and single re-extension of PCR products--a benchmark

Occurrence of chimeric sequences and related artifacts in PCR cloning procedures gives us risks of over-estimation of haplotypes or alleles. Recombination among haplotypes occurs through template switching during PCR cycles or through random repair of mismatch sites on heteroduplex DNA by the host cell. To eliminate the chimeric cloning artifacts, we tested two alternative protocols using T7 endonuclease I cleavage of mismatch sites and re-extension of nascent strands. Though T7 endonuclease I effectively eliminated chimeric clones in some cases, it produced many short fragments. Protocol with single re-extension of PCR products successfully recovered non-recombinant clones with fewer short fragments. In spite of the improvement of allelic recovery through these two protocols, there were still a few recombinants that remained in both reaction mixtures, and thus interpretation of the results for haplotype diversity in a PCR-amplified DNA population should be cautionary. Because re-extension in a diluted reaction mixture is quick, inexpensive and effective, it is advisable to use this procedure for recovery of chromosomal alleles with PCR cloning.

Low MHC variation in the endangered Galápagos penguin (Spheniscus mendiculus)

The major histocompatibility complex (MHC) is one of the most polymorphic regions of the genome, likely due to balancing selection acting to maintain alleles over time. Lack of MHC variability has been attributed to factors such as genetic drift in small populations and relaxed selection pressure. The Galápagos penguin (Spheniscus mendiculus), endemic to the Galápagos Islands, is the only penguin that occurs on the equator. It relies upon cold, nutrient-rich upwellings and experiences severe population declines when ocean temperatures rise during El Niño events. These bottlenecks, occurring in an already small population, have likely resulted in reduced genetic diversity in this species. In this study, we used MHC class II exon 2 sequence data from a DRB1-like gene to characterize the amount of genetic variation at the MHC in 30 Galápagos penguins, as well as one Magellanic penguin (S. magellanicus) and two king penguins (Aptenodytes patagonicus), and compared it to that in five other penguin species for which published data exist. We found that the Galápagos penguin had the lowest MHC diversity (as measured by number of polymorphic sites and average divergence among alleles) of the eight penguin species studied. A phylogenetic analysis showed that Galápagos penguin MHC sequences are most closely related to Humboldt penguin (Spheniscus humboldti) sequences, its putative sister species based on other loci. An excess of non-synonymous mutations and a pattern of trans-specific evolution in the neighbor-joining tree suggest that selection is acting on the penguin MHC.

Study of Cynomolgus monkey (Macaca fascicularis) MhcDRB (Mafa-DRB) polymorphism in two populations

Cynomolgus monkey is one of the macaque species currently used as an animal model for experimental surgery and medicine, in particular, to experiment new drugs or therapy protocols designed for the prevention of allograft rejection. In this field, it is of utmost importance to select histoincompatible recipient-donor pairs. One way to ensure incompatibility between donor and recipient is to check their major histocompatibility complex (MHC) genotypes at the loci playing a determinant role in histocompatibility. We report in this paper on the cynomolgus monkey DRB polymorphism evidenced by sequencing of amplified exon 2 separated either by denaturing gradient gel electrophoresis (DGGE), or by cloning. By the study of 253 unrelated animals from two populations (Mauritius and The Philippines), we characterized 50 exon 2 sequences among which 28 were identical to sequences already reported in Macaca fascicularis or other macaque species (Macaca mulatta, Macaca nemestrina). By cloning and sequencing DRB cDNA, we revealed two additional DRB alleles. Out of the 20 haplotypes that we defined here, only two were found in both populations. The functional impact of DR incompatibility was studied in vitro by mixed lymphocyte culture.

Major histocompatibility complex and microsatellite variation in two populations of wild gorillas

In comparison to their close relatives the chimpanzees and humans, very little is known concerning the amount and structure of genetic variation in gorillas. Two species of gorillas are recognized and while the western gorillas number in the tens of thousands, only several hundred representatives of the mountain gorilla subspecies of eastern gorillas survive. To analyse the possible effects of these different population sizes, this study compares the variation observed at microsatellite and major histocompatibility complex (MHC) loci in samples of wild western and mountain gorillas, collected using a sampling scheme that targeted multiple social groups within defined geographical areas. Noninvasive samples proved a viable source of DNA for sequence analysis of the second exon of the DRB loci of the MHC. Observed levels of variation at the MHC locus were similar between the two gorilla species and were comparable to those in other primates. Comparison of results from analysis of variation at multiple microsatellite loci found only a slight reduction in heterozygosity for the mountain gorillas despite the relatively smaller population size.

Denaturing gradient gel electrophoresis and its use in the detection of major histocompatibility complex polymorphism

The major histocompatibility complex (MHC) has been studied extensively in humans and in mice and many methods are available for MHC typing of these well-characterized species. Studies of MHC variation in other species are ever increasing and researchers can choose one of a number of approaches for MHC typing of their species of interest. DNA sequencing is regarded as the 'gold standard' and it is frequently used for MHC typing. However, DNA sequencing is impractical when many individuals must be typed. Denaturing gradient gel electrophoresis (DGGE) offers a flexible and sensitive method for identifying and characterizing MHC alleles in any vertebrate species. This article reviews the theory and the practice of DGGE and examines the use of DGGE for MHC identification in various species. DGGE is compared to other similar techniques for MHC typing, such as single-stranded conformational polymorphism and reference strand-mediated conformational analysis. The advantages, problems, pitfalls and limitations of DGGE are considered and future perspectives on the use of DGGE for MHC typing are discussed.

Natural selection of the major histocompatibility complex (Mhc) in Hawaiian honeycreepers (Drepanidinae)

The native Hawaiian honeycreepers represent a classic example of adaptive radiation and speciation, but currently face one the highest extinction rates in the world. Although multiple factors have likely influenced the fate of Hawaiian birds, the relatively recent introduction of avian malaria is thought to be a major factor limiting honeycreeper distribution and abundance. We have initiated genetic analyses of class II beta chain Mhc genes in four species of honeycreepers using methods that eliminate the possibility of sequencing mosaic variants formed by cloning heteroduplexed polymerase chain reaction products. Phylogenetic analyses group the honeycreeper Mhc sequences into two distinct clusters. Variation within one cluster is high, with dN > dS and levels of diversity similar to other studies of Mhc (B system) genes in birds. The second cluster is nearly invariant and includes sequences from honeycreepers (Fringillidae), a sparrow (Emberizidae) and a blackbird (Emberizidae). This highly conserved cluster appears reminiscent of the independently segregating Rfp-Y system of genes defined in chickens. The notion that balancing selection operates at the Mhc in the honeycreepers is supported by transpecies polymorphism and strikingly high dN/dS ratios at codons putatively involved in peptide interaction. Mitochondrial DNA control region sequences were invariant in the i'iwi, but were highly variable in the 'amakihi. By contrast, levels of variability of class II beta chain Mhc sequence codons that are hypothesized to be directly involved in peptide interactions appear comparable between i'iwi and 'amakihi. In the i'iwi, natural selection may have maintained variation within the Mhc, even in the face of what appears to a genetic bottleneck.

The use of reference strand-mediated conformational analysis for the study of cheetah (Acinonyx jubatus) feline leucocyte antigen class II DRB polymorphisms

There is now considerable evidence to suggest the cheetah (Acinonyx jubatus) has limited genetic diversity. However, the extent of this and its significance to the fitness of the cheetah population, both in the wild and captivity, is the subject of some debate. This reflects the difficulty associated with establishing a direct link between low variability at biologically significant loci and deleterious aspects of phenotype in this, and other, species. Attempts to study one such region, the feline leucocyte antigen (FLA), are hampered by a general reliance on cloning and sequencing which is expensive, labour-intensive, subject to PCR artefact and always likely to underestimate true variability. In this study we have applied reference strand-mediated conformational analysis (RSCA) to determine the FLA-DRB phenotypes of 25 cheetahs. This technique was rapid, repeatable and less prone to polymerase chain reaction (PCR)-induced sequence artefacts associated with cloning. Individual cheetahs were shown to have up to three FLA-DRB genes. A total of five alleles were identified (DRB*ha14-17 and DRB*gd01) distributed among four genotypes. Fifteen cheetahs were DRB*ha14/ha15/ha16/ha17, three were DRB*ha15/ha16/ha17, six were DRB*ha14/ha16/ha17 and one was DRB*ha14/ha15/ha16/ha17/gd01. Sequence analysis of DRB*gd01 suggested it was a recombinant of DRB*ha16 and DRB*ha17. Generation of new alleles is difficult to document, and the clear demonstration of such an event is unusual. This study confirms further the limited genetic variability of the cheetah at a biologically significant region. RSCA will facilitate large-scale studies that will be needed to correlate genetic diversity at such loci with population fitness in the cheetah and other species.

Resolution of complex feline leukocyte antigen DRB loci by reference strand-mediated conformational analysis (RSCA)

The DRB genes of the domestic cat are highly polymorphic. Studies based on clonal sequence analysis have suggested the existence of two distinct loci within individual animals and good evidence for 24 distinct FLA-DRB alleles. This variability, the complexity of clonal sequence analysis and its susceptibility to PCR-induced artefacts has represented a bottleneck to further progress. In this study we have applied reference strand-mediated conformational analysis (RSCA) to FLA-DRB. This protocol has been shown to be highly reproducible. Using five reference strands including two derived from non-domestic felines, we could distinguish 23 FLA-DRB alleles. We used RSCA to explore genetic polymorphism of FLA-DRB in 71 cats including 31 for which clonal sequence analysis was also available. On average, RSCA identified 0.9 more alleles within cats than clonal sequence analysis. Reference strand-mediated conformational analysis was also able to identify animals containing new alleles that could be targeted for sequence analysis. Analysis of allele patterns showed clear evidence for different allele distributions between breeds of cats, and suggested the Burmese breed may have highly restricted FLA-DRB polymorphism. Results from two families provided clear evidence for variation in the number of DRB genes on different haplotypes, with some haplotypes carrying two genes and some containing three. This study highlights the utility of RSCA for the resolution of complex amplicons containing up to six distinct alleles. A simple, rapid method for characterizing FLA-DRB makes possible studies on vaccine response and susceptibility/resistance to viral infections, which are a significant clinical problem in cats.