Haplotype variation, recombination, and gene conversion within the turkey MHC-B locus

The major histocompatibility complex (MHC) is a gene dense region with profound effects on the disease phenotype. In many species, characterizations of MHC polymorphisms have focused on identifying allelic haplotypes of the highly polymorphic class I and class II loci through direct immunological approaches such as monoclonal antibodies specific for the major antigens or indirectly through DNA sequence-based approaches. Invariably, these studies fail to assess the broader range of variation at the other loci within the MHC. This study examines variation in the turkey MHC by resequencing 15 interspersed amplicons (∼14 kb) spaced across the MHC-B locus in a representative sampling of 52 commercial birds. Over 200 single nucleotide polymorphisms (SNPs) were identified with high levels of polymorphism (1 SNP/70 bp) and heterozygosity (average minor allele frequency of 0.15). SNP genotypes were used to identify the major haplotypes segregating in the commercial lines. Sequencing of the peptide binding region (PBR, exon 2) of the class IIB loci of select individuals identified 10 PBR alleles/isotypes among the major MHC haplotypes. Examination of pedigreed families provides direct evidence of gene conversion and recombination within the B locus. Results of this study demonstrate the MHC diversity available in commercial flocks and provide genomic resources for studying the effect of this diversity (alleles and/or haplotypes) on disease susceptibility and resistance.     

Genetic Variation at the MHC in a Population of Introduced Wild Turkeys

Genetic variation in the major histocompatibility complex (MHC) is known to affect disease resistance in many species. Investigations of MHC diversity in populations of wild species have focused on the antigen presenting class IIβ molecules due to the known polymorphic nature of these genes and the role these molecules play in pathogen recognition. Studies of MHC haplotype variation in the turkey (Meleagris gallopavo) are limited. This study was designed to examine MHC diversity in a group of Eastern wild turkeys (Meleagris gallopavo silvestris) collected during population expansion following reintroduction of the species in southern Wisconsin, USA. Southern blotting with BG and class IIβ probes and single nucleotide polymorphism (SNP) genotyping was used to measure MHC variation. SNP analysis focused on single copy MHC genes flanking the highly polymorphic class IIβ genes. Southern blotting identified 27 class IIβ phenotypes, whereas SNP analysis identified 13 SNP haplotypes occurring in 28 combined genotypes. Results show that genetic diversity estimates based on RFLP (Southern blot) analysis underestimate the level of variation detected by SNP analysis. Sequence analysis of the mitochondrial D-loop identified 7 mitochondrial haplotypes (mitotypes) in the sampled birds. Results show that wild turkeys located in southern Wisconsin have a genetically diverse MHC and originate from several maternal lineages.     

Extended sequence of the turkey MHC <Emphasis Type="Italic">B</Emphasis>-locus and sequence variation in the highly polymorphic B-G loci

Genetic variation in the major histocompatibility complex (MHC) is directly correlated to differences in disease resistance. Immunity is greatly dependent on highly polymorphic genes in the MHC, such as class I, class II, and class III complement genes. Preliminary studies of wild turkey populations show extreme polymorphisms in a family of genes exclusive to the avian MHC, the class IV or B-G genes. Significance of this variation is unclear as there are few and conflicting studies of the expression of these genes. Confounding understanding of B-G variation is the lack of a complete delineation of the number of loci in the turkey genome. Direct 454 sequencing of a clone from the CHORI-260 BAC library was used to extend the turkey MHC B-locus sequence, identifying five additional complete B-locus genes including two B-G loci. Sequences of the new B-G genes were compared with those of other turkey gene (BG1–3) and sequences available for other galliformes. Phylogenetic analysis shows species-specific gene evolution supporting a birth–death model of evolution for the B-G gene family. Analysis of variation within the signal peptide sequence (exon 1) found two clusters of polymorphism among the turkey B-G genes. Resequencing of exon 1 in a diverse sample including wild, heritage, and commercial turkeys confirmed multiple alleles at each B-G gene. Future studies aim to correlate B-G variation with group and individual immunological differences.     

Analysis of MHC class I genes across horse MHC haplotypes

The genomic sequences of 15 horse major histocompatibility complex (MHC) class I genes and a collection of MHC class I homozygous horses of five different haplotypes were used to investigate the genomic structure and polymorphism of the equine MHC. A combination of conserved and locus-specific primers was used to amplify horse MHC class I genes with classical and nonclassical characteristics. Multiple clones from each haplotype identified three to five classical sequences per homozygous animal and two to three nonclassical sequences. Phylogenetic analysis was applied to these sequences, and groups were identified which appear to be allelic series, but some sequences were left ungrouped. Sequences determined from MHC class I heterozygous horses and previously described MHC class I sequences were then added, representing a total of ten horse MHC haplotypes. These results were consistent with those obtained from the MHC homozygous horses alone, and 30 classical sequences were assigned to four previously confirmed loci and three new provisional loci. The nonclassical genes had few alleles and the classical genes had higher levels of allelic polymorphism. Alleles for two classical loci with the expected pattern of polymorphism were found in the majority of haplotypes tested, but alleles at two other commonly detected loci had more variation outside of the hypervariable region than within. Our data indicate that the equine major histocompatibility complex is characterized by variation in the complement of class I genes expressed in different haplotypes in addition to the expected allelic polymorphism within loci.     

Polymorphism,haplotype composition,and selection in the <Emphasis Type="Italic">Mhc-DRB</Emphasis> of wild baboons

General patterns of organization in the major histocompatibility complex (MHC) have been successfully explained by the model of birth-and-death evolution, but understanding why certain MHC genes are maintained together into specific haplotypes remains challenging. The haplotype configurations of the functionally important class II DR region have been described in few primates and display important interspecific variability with respect to the extent of allelic variation, the number of loci and/or combinations of loci present. Understanding the evolutionary mechanisms driving such variation is conditional upon characterizing haplotypes in new species and identifying the selective pressures acting on haplotypes. This study explores the variability of haplotype configurations in the Mhc-DRB region (exon 2) for the first time in wild non-human primates, chacma baboons (Papio ursinus). Paur-DRB haplotypes were characterized through segregation studies and linkage disequilibrium. 23 Paur-DRB sequences and 15 haplotype configurations were identified in 199 animals. The Paur-DRB exon 2 is shown to be subjected to intense positive selection and frequent recombination. An approach recently developed for human vaccine studies was used to classify Paur-DRB sequences into supertypes, based on the physico-chemical properties of amino acids that are positively selected, thus most probably involved in antigen recognition. Sequences grouped into the same supertype (thus presumably sharing antigen-binding affinities) are non-randomly distributed within haplotypes, leading to an increased individual diversity of supertypes. Our results suggest that selection favoring haplotypes with complementary sets of DRB supertypes shapes functionally tuned haplotypes in this natural baboon population.     

Characterization of the turkey MHC chromosome through genetic and physical mapping

Previous studies in the chicken have identified a single microchromosome (GGA16) containing the ribosomal DNA (rDNA) and two genetically unlinked MHC regions, MHC-B and MHC-Y. Chicken DNA sequence from these loci was used to develop PCR primers for amplification of homologous fragments from the turkey (Meleagris gallopavo). PCR products were sequenced and overgo probes were designed to screen the CHORI 260 turkey BAC library. BAC clones corresponding to the turkey rDNA, MHC-B and MHC-Y were identified. BAC end and subclone sequencing confirmed identity and homology of the turkey BAC clones to the respective chicken loci. Based on subclone sequences, single-nucleotide polymorphisms (SNPs) segregating within the UMN/NTBF mapping population were identified and genotyped. Analysis of SNP genotypes found the B and Y to be genetically unlinked in the turkey. Silver staining of metaphase chromosomes identified a single pair of microchromosomes with nucleolar organizer regions (NORs). Physical locations of the rDNA and MHC loci were determined by fluorescence in situ hybridization (FISH) of the BAC clones to metaphase chromosomes. FISH clearly positioned the rDNA distal to the Y locus on the q-arm of the MHC chromosome and the MHC-B on the p-arm. An internal telomere array on the MHC chromosome separates the B and Y loci.     

Toward understanding MHC disease associations: partial resequencing of 46 distinct HLA haplotypes

We carried out a resequencing project that examined 552 kb of sequence from each of 46 individual HLA haplotypes representing a diversity of HLA allele types, generating nearly 27 Mb of fully phased genomic sequence. Haplotype blocks were defined extending from telomeric of HLA-F to centromeric of HLA-DP including in total 5186 MHC SNPs. To investigate basic questions about the evolutionary origin of common HLA haplotypes, and to obtain an estimate of rare variation in the MHC, we similarly examined two additional sets of samples. In 19 independent HLA-A1, B8, DR3 chromosomes, the most common HLA haplotype in Northern European Caucasians, variation was found at 11 SNP positions in the 3600-kb region from HLA-A to DR. Partial resequencing of 282 individuals in the gene-dense class III region identified significant variability beyond what could have been detected by linkage to common SNPs.     

Multiple divergent haplotypes express completely distinct sets of class I MHC genes in zebrafish

The zebrafish is an important animal model for stem cell biology, cancer, and immunology research. Histocompatibility represents a key intersection of these disciplines; however, histocompatibility in zebrafish remains poorly understood. We examined a set of diverse zebrafish class I major histocompatibility complex (MHC) genes that segregate with specific haplotypes at chromosome 19, and for which donor-recipient matching has been shown to improve engraftment after hematopoietic transplantation. Using flanking gene polymorphisms, we identified six distinct chromosome 19 haplotypes. We describe several novel class I U lineage genes and characterize their sequence properties, expression, and haplotype distribution. Altogether, ten full-length zebrafish class I genes were analyzed, mhc1uba through mhc1uka. Expression data and sequence properties indicate that most are candidate classical genes. Several substitutions in putative peptide anchor residues, often shared with deduced MHC molecules from additional teleost species, suggest flexibility in antigen binding. All ten zebrafish class I genes were uniquely assigned among the six haplotypes, with dominant or codominant expression of one to three genes per haplotype. Interestingly, while the divergent MHC haplotypes display variable gene copy number and content, the different genes appear to have ancient origin, with extremely high levels of sequence diversity. Furthermore, haplotype variability extends beyond the MHC genes to include divergent forms of psmb8. The many disparate haplotypes at this locus therefore represent a remarkable form of genomic region configuration polymorphism. Defining the functional MHC genes within these divergent class I haplotypes in zebrafish will provide an important foundation for future studies in immunology and transplantation.     

A fast and accurate enumeration-based algorithm for haplotyping a triploid individual

Background

Haplotype assembly, reconstructing haplotypes from sequence data, is one of the major computational problems in bioinformatics. Most of the current methodologies for haplotype assembly are designed for diploid individuals. In recent years, genomes having more than two sets of homologous chromosomes have attracted many research groups that are interested in the genomics of disease, phylogenetics, botany and evolution. However, there is still a lack of methods for reconstructing polyploid haplotypes.

Results

In this work, the minimum error correction with genotype information (MEC/GI) model, an important combinatorial model for haplotyping a single individual, is used to study the triploid individual haplotype reconstruction problem. A fast and accurate enumeration-based algorithm enumeration haplotyping triploid with least difference (EHTLD) is proposed for solving the MEC/GI model. The EHTLD algorithm tries to reconstruct the three haplotypes according to the order of single nucleotide polymorphism (SNP) loci along them. When reconstructing a given SNP site, the EHTLD algorithm enumerates three kinds of SNP values in terms of the corresponding site’s genotype value, and chooses the one, which leads to the minimum difference between the reconstructed haplotypes and the sequenced fragments covering that SNP site, to fill the SNP loci being reconstructed.

Conclusion

Extensive experimental comparisons were performed between the EHTLD algorithm and the well known HapCompass and HapTree. Compared with algorithms HapCompass and HapTree, the EHTLD algorithm can reconstruct more accurate haplotypes, which were proven by a number of experiments.

     

Genetic analysis of completely sequenced disease-associated MHC haplotypes identifies shuffling of segments in recent human history

The major histocompatibility complex (MHC) is recognised as one of the most important genetic regions in relation to common human disease. Advancement in identification of MHC genes that confer susceptibility to disease requires greater knowledge of sequence variation across the complex. Highly duplicated and polymorphic regions of the human genome such as the MHC are, however, somewhat refractory to some whole-genome analysis methods. To address this issue, we are employing a bacterial artificial chromosome (BAC) cloning strategy to sequence entire MHC haplotypes from consanguineous cell lines as part of the MHC Haplotype Project. Here we present 4.25 Mb of the human haplotype QBL (HLA-A26-B18-Cw5-DR3-DQ2) and compare it with the MHC reference haplotype and with a second haplotype, COX (HLA-A1-B8-Cw7-DR3-DQ2), that shares the same HLA-DRB1, -DQA1, and -DQB1 alleles. We have defined the complete gene, splice variant, and sequence variation contents of all three haplotypes, comprising over 259 annotated loci and over 20,000 single nucleotide polymorphisms (SNPs). Certain coding sequences vary significantly between different haplotypes, making them candidates for functional and disease-association studies. Analysis of the two DR3 haplotypes allowed delineation of the shared sequence between two HLA class II-related haplotypes differing in disease associations and the identification of at least one of the sites that mediated the original recombination event. The levels of variation across the MHC were similar to those seen for other HLA-disparate haplotypes, except for a 158-kb segment that contained the HLA-DRB1, -DQA1, and -DQB1 genes and showed very limited polymorphism compatible with identity-by-descent and relatively recent common ancestry (<3,400 generations). These results indicate that the differential disease associations of these two DR3 haplotypes are due to sequence variation outside this central 158-kb segment, and that shuffling of ancestral blocks via recombination is a potential mechanism whereby certain DR-DQ allelic combinations, which presumably have favoured immunological functions, can spread across haplotypes and populations.     

Diversity and locus specificity of chicken MHC B class I sequences

The major histocompatibility complex B (MHC B) region in a standard haplotype of Leghorn chickens contains two closely linked class I loci, B-FI and B-FIV. Few sequences of B-FI alleles are available, and therefore alleles of the two loci have not been compared with regard to sequence diversity or locus specificity. Here, we report eight new B-F alpha 1/alpha 2-coding sequences from broiler chicken MHC B haplotypes, and a unique recombinant between the two B-F loci. The new sequences were combined with existing B-F sequences from Leghorn and broiler haplotypes for analysis. On the basis of phylogenetic analysis and conserved sequence motifs, B-F sequences separated into two groups (Groups A and B), corresponding to B-FIV and B-FI locus, respectively. Every broiler haplotype had one B-F sequence in Group A and the second B-F sequence, if it existed, clustered in Group B. Group B (presumptive B-FI locus) sequences identified in broiler haplotypes resembled the human MHC class I HLA-C locus in their distinctive pattern of allelic polymorphism. Compared with B-FIV, B-FI alleles were less polymorphic and possessed a conserved locus-specific motif in the alpha1 helix, but nevertheless demonstrated evidence of diversifying selection. One B-FI alpha 1/alpha 2-coding nucleotide sequence was completely conserved in four different broiler haplotypes, but each allele differed in the exon encoding the alpha 3 domain.     

DNA sequence and haplotype variation in two candidate genes for dilated cardiomyopathy in the turkey Meleagris gallopavo

Determining variation in genes is fundamental to understanding their function in the disease state. Cardiac troponin T (cTnT) and phospholamban (PLN) genes have been implicated in dilated cardiomyopathy (DCM) in human and model species. To investigate the role of these 2 candidate genes in DCM in the turkey Meleagris gallopavo, understanding sequence variants and map position distribution is necessary. To this end, a total of 1854 and 1771 bp of cTnT and PLN gene sequences, respectively, were scanned for single nucleotide polymorphisms (SNPs) in a randomly bred population. A total of 15 SNPs was identified in the cTnT and PLN genomic sequences. Nine haplotypes, 5 in cTnT and 4 in PLN, were identified. Observed heterozygosities (0.02-0.39) in the turkey population were low for both genes. Within each gene, 1 SNP corresponding to a restriction enzyme site was identified and used to develop a PCR-restriction fragment length polymorphism (RFLP) genotyping assay. The PLN gene was genetically mapped to turkey chromosome 2, equivalent to Gallus gallus chromosome 3, and cTnT mapped to a turkey microchromosome. Although limited because of the relatively small sample size of 55 birds, the data from this SNP analysis of PLN and cTnT provide a foundation from which to evaluate the function of cTnT and PLN in the turkey. Information about the distribution of the SNPs and haplotypes will facilitate future association and linkage studies.     

Haplotype diversity in 11 candidate genes across four populations

Analysis of haplotypes based on multiple single-nucleotide polymorphisms (SNP) is becoming common for both candidate gene and fine-mapping studies. Before embarking on studies of haplotypes from genetically distinct populations, however, it is important to consider variation both in linkage disequilibrium (LD) and in haplotype frequencies within and across populations, as both vary. Such diversity will influence the choice of "tagging" SNPs for candidate gene or whole-genome association studies because some markers will not be polymorphic in all samples and some haplotypes will be poorly represented or completely absent. Here we analyze 11 genes, originally chosen as candidate genes for oral clefts, where multiple markers were genotyped on individuals from four populations. Estimated haplotype frequencies, measures of pairwise LD, and genetic diversity were computed for 135 European-Americans, 57 Chinese-Singaporeans, 45 Malay-Singaporeans, and 46 Indian-Singaporeans. Patterns of pairwise LD were compared across these four populations and haplotype frequencies were used to assess genetic variation. Although these populations are fairly similar in allele frequencies and overall patterns of LD, both haplotype frequencies and genetic diversity varied significantly across populations. Such haplotype diversity has implications for designing studies of association involving samples from genetically distinct populations.     

Molecular genotype identification of the Gallus gallus major histocompatibility complex

The chicken major histocompatibility complex (MHC) is commonly defined by serologic reactions of erythrocytes with antibodies specific to the highly polymorphic MHC class I (BF) and MHC class IV (BG) antigens. The microsatellite marker LEI0258 is known to be physically located within the MHC, between the BG and BF regions. DNA from various serologically defined MHC haplotypes was amplified by polymerase chain reaction with primers surrounding this marker. Twenty-six distinctive allele sizes were identified. Some serologically well-defined MHC haplotypes shared a common LEI0258 allele size but could be distinguished either by the addition of information from another nearby marker (MCW0371) or by small indels or single nucleotide polymorphism (SNP) differences between the alleles. The association between LEI0258 allele and serologically defined MHC haplotype was very consistent for the same haplotype from multiple sources. Sequence information for the region defined by LEI0258 was obtained for 51 different haplotypes. Two internal repeats whose lengths were 13 and 12 bp, respectively, are the primary basis for allelic variability. Allele size variation ranges from 182 to 552 bp. Four indels and five SNPs in the surrounding sequence provide additional means for distinguishing alleles. Typing with LEI0258 and MCW0371 will be useful in identifying MHC haplotypes in outbred populations of chickens particularly for the initial development of serological reagents.     

Balancing selection, random genetic drift, and genetic variation at the major histocompatibility complex in two wild populations of guppies (Poecilia reticulata)

Our understanding of the evolution of genes of the major histocompatibility complex (MHC) is rapidly increasing, but there are still enigmatic questions remaining, particularly regarding the maintenance of high levels of MHC polymorphisms in small, isolated populations. Here, we analyze the genetic variation at eight microsatellite loci and sequence variation at exon 2 of the MHC class IIB (DAB) genes in two wild populations of the Trinidadian guppy, Poecilia reticulata. We compare the genetic variation of a small (Ne, 100) and relatively isolated upland population to that of its much larger (Ne approximately 2400) downstream counterpart. As predicted, microsatellite diversity in the upland population is significantly lower and highly differentiated from the population further downstream. Surprisingly, however, these guppy populations are not differentiated by MHC genetic variation and show very similar levels of allelic richness. Computer simulations indicate that the observed level of genetic variation can be maintained with overdominant selection acting at three DAB loci. The selection coefficients differ dramatically between the upland (s > or = 0.2) and lowland (s < or = 0.01) populations. Parasitological analysis on wild-caught fish shows that parasite load is significantly higher on upland than on lowland fish, which suggests that large differences in selection intensity may indeed exist between populations. Based on the infection intensity, a substantial proportion of the upland fish would have suffered direct or indirect fitness consequences as a result of their high parasite loads. Selection by parasites plays a particularly important role in the evolution of guppies in the upland habitat, which has resulted in high levels of MHC diversity being maintained in this population despite considerable genetic drift.     

Considerable haplotypic diversity in the RT1-CE class I gene region of the rat major histocompatibility complex

The major histocompatibility complex (MHC) class I region extending between the Bat1 and Pou5f1 genes shows considerable genomic plasticity in mouse and rhesus macaque but not in human haplotypes. In the rat, this region is known as the RT1-CE region. The recently published rat MHC sequence gave rise to a complete set of class I gene sequences in a single MHC haplotype, namely the RT1ⁿ haplotype of the widely used BN inbred strain. To study the degree of genetic diversity, we compared the RT1-CE region-derived class I genes of the RT1ⁿ haplotype with class I sequences of other rat haplotypes. By using phylogenetic tree analyses, we obtained evidence for extensive presence and absence polymorphisms of single loci and even small subfamilies of class I genes in the rat. Alleles of RT1-CE region class I genes could also be identified, but the rate of allelic nucleotide substitutions appeared rather low, indicating that the diversity in the RT1-CE region is mainly based on genomic plasticity.     

Estimation of DNA sequence diversity in bovine cytokine genes

DNA sequence variation provides the fundamental material for improving livestock through selection. In cattle, single nucleotide polymorphisms and small insertions/deletions (collectively referred to here as SNPs) have been identified in cytokine genes and scored in a reference population to determine linkage map positions. The aim of the present study was twofold: first, to estimate the SNP frequency in a reference population of beef cattle, and second, to determine cytokine haplotypes in a group of sires from commercial populations. Forty-five SNP markers in DNA segments from nine cytokine gene loci were analyzed in 26 reference parents. Comparison of all 52 haploid genomes at each PCR amplicon locus revealed an average of one SNP per 143 bp of sequence, whereas comparison of any two chromosomes identified heterozygous sites, on average, every 443 bp. The combination of these 45 SNP genotypes was sufficient to uniquely identify each of the 26 animals. The average number of haplotype alleles (4.4) per PCR amplicon (688 bp) and the percentage heterozygosity among founding parents (50%) were similar to those for microsatellite markers in the same population. For 49 sires from seven common breeds of beef cattle, SNP genotypes (1225 total) were obtained by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) at three amplicon loci. All three of the amplicon haplotypes were correctly deduced for each sire without the use of parent or progeny genotypes. The latter allows a wide range of genetic studies in commercial populations of cattle where genotypic information from relatives may not be available. Received: 16 June 2000 / Accepted: 23 August 2000