Class II genes of miniature swine. IV. Characterization and expression of two allelic class II DQB cDNA clones

Two cDNA clones coding for allelic miniature swine MHC class II Ag DQB chains have been isolated, characterized, and shown to be expressed after transfection into mouse fibroblasts. The two alleles differ at the nucleotide level by an overwhelming proportion of replacement substitutions, suggesting the influence of selection for polymorphism. Most of the resulting predicted amino acid replacements are in regions commonly polymorphic in mouse Ab and human DQB sequences, corresponding to the predicted Ag recognition site. Nucleotide and amino acid sequence comparisons to homologous mouse and human sequences show more similarity between swine and man than between either swine and mouse or man and mouse. This tendency is most pronounced when comparing the 3' untranslated regions. However, an examination of unique cross-species sharing of amino acid residues suggests a closer relationship between both man and miniature swine and man and mouse than between miniature swine and mouse. The simplest explanation we can envision for these findings is that the mouse DQB gene homologue (Ab) has been subject to a higher substitution rate than either swine or human DQB genes. An additional cytoplasmic exon expressed in mouse Ab gene products and in putative human DQB2 gene products is lacking in both swine and human DQB cDNA clones. Its absence suggests either that the expression of this exon in mouse Ab genes was activated after mammalian speciation or that the expression of this exon was independently inactivated in swine DQB and human DQB1 genes. Alternatively, the mouse Ab gene may be derived from the same primordial gene as human DQB2, whereas the pig DQB gene may be derived from the same primordial gene as the human DQB1 gene.     

Polymorphism and gene organization of water buffalo MHC-DQB genes show homology to the BoLA DQB region

In cattle (Bos taurus), there is evidence of more than 50 alleles of BoLA-DQB (bovine lymphocyte antigen DQB) that are distributed across at least five DQB loci, making this region one of the most complex in the BoLA gene family. In this study, DQB alleles were analysed for the water buffalo (Bubalus bubalis), another economically important bovine species. Twelve alleles for Bubu-DQB (Bubalis bubalis DQB) were determined by nucleotide sequence analysis. A phylogenetic analysis revealed numerous trans-species polymorphisms, with alleles from water buffalo assigned to at least three different loci (BoLA-DQB1, BoLA-DQB3 and BoLA-DQB4) that are also found in cattle. These presumptive loci were analysed for patterns of synonymous (d(S)) and non-synonymous (d(N)) substitution. Like BoLA-DQB1, Bubu-DQB1 was observed to be under strong positive selection for polymorphism. We conclude that water buffalo and cattle share the current arrangement of their DQB region because of their common ancestry.     

High similarity at three MHC loci between the baiji and finless porpoise: trans-species or convergent evolution?

Two DRA alleles and six MHC-I alleles were identified from a group of 15 baiji (Lipotes vexillifer), the most threatened cetacean in the world. Little sequence variation was detected at the DRA locus but extensive variation at the MHC-I locus. In combination with data at the DQB locus previously reported, three MHC loci exon 2 of the baiji all revealed striking similarity with those of the finless porpoise. Especially, some identical alleles shared by both species at the MHC-I and DQB loci suggested the convergent evolution as a consequence of common adaptive solutions to similar environmental pressures in the Yangtze River. As for DRA locus, the identity alleles were shared not only by baiji and finless porpoise but by some other cetacean species of the families Phocoenidae and Delphinidae, suggesting trans-species evolution on this gene.     

A second locus and new alleles in the major histocompatibility complex class II (ELA-DQB) region in the horse

More than two nucleotide sequences of the second exon of the ELA-DQB region retrieved from a single animal and two different sequences isolated from horses homozygous in the major histocompatibility complex (MHC) region by descent indicated the existence of at least two ELA-DQB loci at the genomic level. New alleles detected by polymerase chain reaction single strand conformation polymorphism (SSCP) and defined by nucleotide sequencing of the second exon of the DQB gene(s) were described. Based on the level of nucleotide sharing, at least two groups of alleles were shown to exist. The newly defined alleles belonged preferentially to one of the groups. However, their specific locus assignment was not possible from the data collected. At least one of these alleles was shown to be transcribed. No frame-shift mutations were identified among the new alleles, although one pseudoallele containing a stop codon was identified at the genomic DNA level.     

Mouse models and the interpretation of human GWAS in type 2 diabetes and obesity

Within the last 3 years, genome-wide association studies (GWAS) have had unprecedented success in identifying loci that are involved in common diseases. For example, more than 35 susceptibility loci have been identified for type 2 diabetes and 32 for obesity thus far. However, the causal gene and variant at a specific linkage disequilibrium block is often unclear. Using a combination of different mouse alleles, we can greatly facilitate the understanding of which candidate gene at a particular disease locus is associated with the disease in humans, and also provide functional analysis of variants through an allelic series, including analysis of hypomorph and hypermorph point mutations, and knockout and overexpression alleles. The phenotyping of these alleles for specific traits of interest, in combination with the functional analysis of the genetic variants, may reveal the molecular and cellular mechanism of action of these disease variants, and ultimately lead to the identification of novel therapeutic strategies for common human diseases. In this Commentary, we discuss the progress of GWAS in identifying common disease loci for metabolic disease, and the use of the mouse as a model to confirm candidate genes and provide mechanistic insights.     

Genetic variability of MHC class II DQB exon 2 alleles in yak (Bos grunniens)

The major histocompatibility class (MHC) DQ molecules are dimeric glycoproteins revealing antigen presentation to CD⁴⁺ T cells. In the present study, the exon 2 of the MHC class II DQB gene from 32 yaks (Bos grunniens) was cloned, sequenced and compared with previously reported patterns for other bovidae. It was revealed by sequence analyses that there are 25 DQB exon 2 alleles among 32 yaks, all alleles are found to belong to DQB1 loci. These alleles exhibited a high degree of nucleotide and amino acid polymorphisms with most amino acid variations occurring at positions forming the peptide-binding sites. The DQB loci were analyzed for patterns of synonymous (d _S) and non-synonymous (d _N) substitution. The yak was observed to be under strong positive selection in the DQB exon 2 peptide-binding sites (d _N = 0.15, P < 0.001). It appears that this variability among yaks confers the ability to mount immune responses to a wide variety of peptides or pathogens.     

Pyrosequence-based typing of alleles of the HLA-DQB1 gene

DNA typing of alleles of the highly polymorphic HLA-DQBI gene was performed by Pyrosequencing using purified DNA from the 11th International Histocompatibility Workshop human cell lines and samples from the Children's Hospital of Pittsburgh registry of diabetics and their first-degree relatives. Pyrosequencing was optimized for genotyping exon 2 of the HLA-DQB1 gene, but the procedure should be applicable to other HLA loci. The 47 HLA-DQB1 alleles were readily identifiable, as were the 1,128 potential allelic heterozygous combinations. The method required PCR conditions that specifically amplified DQB1 but not the pseudogene, DQB2. The new method of pyrosequence-based typing can be performed in 96- or 384-well format. The 61 polymorphic residues of DQB1 exon 2 were identified within four pyrosequencing reactions, obtained by a 70-nucleotide read length in each reaction, in about an hour's time. Allelic combinations of HLA-DQB1 most frequentlyfound in the population of diabetics and their immediate family members were analyzed and successfully compared to typing of the DQB1 alleles by sequence-specific oligonucleotide probe protocols. Pyrosequence-based typing is compatible with genotyping of allelic combinations expected from heterozygous individuals, resulting in nucleotide resolution of the highly polymorphic HLA system. Using pyrosequencing, more than 750 sample wells can be processed in a working day, resulting in the identification of more than 50,000 bases.     

Development and characterization of eleven polymorphic microsatellite loci from South China field mouse (Apodemus draco)

One population distributed in Yunnan of China was regarded as a new species based on mitochondrial cytochrome b sequences. However, the usefulness of mitochondrial sequence data in determining species boundaries is not universally agreed upon, the frequency data from multiple nuclear gene loci is necessary in determining species boundaries. So, we describe in this paper the isolation and characterization of eleven microsatellite loci in the South China field mouse from genomic DNA-enriched libraries. The eleven loci were tested in 24 individuals from two populations in Southwest China. These loci were highly polymorphic with numbers of alleles per locus ranging from 9 to 24 and expected heterozygosities from 0.898 to 0.967. Eight loci followed Hardy–Weinberg expectations after Bonferroni correction for multiple comparisons. No significant linkage association was found among all these loci. The eleven polymorphic microsatellite loci will be useful in determining species boundaries of the South China field mouse.     

Gene duplications and sequence polymorphism of bovine class II DQB genes

The genetic diversity of bovine class II DQB genes was investigated by polymerase chain reaction amplification and DNA sequencing. The first domain exon was amplified from genomic DNA samples representing 14 class II haplotypes, defined by restriction fragment length polymorphism (RFLP) analysis. The presence of a polymorphism in the copy number of DQB genes was confirmed since two DQB sequences were isolated from certain haplotypes. Four subtypes of bovine DQB genes were found. DQB1 is the major type and was found in almost all haplotypes. DQB2 is very similar to DQB1 but was found only in the duplicated haplotypes DQ9 to 12. DQB3 and DQB4 are two quite divergent genes only present in certain duplicated haplotypes. The bovine DQB complexity thus resembles that in the human DRB region. Bovine DQB genes were found to be highly polymorphic as ten DQB1 alleles and four DQB2 alleles were identified. The observed sequence polymorphism correlated well with previously defined DQB RFLPs. Bovine and human DQB alleles show striking similarities at the amino acid level. In contrast, the frequency of silent substitutions is much higher in comparisons of DQB alleles between species than within species ruling out the possibility that any of the contemporary DQB alleles have been maintained since the divergence of humans and cattle. The frequency of silent substitutions between DQB alleles was markedly lower in cattle than in humans, in agreement with a previous comparison of human and bovine DRB alleles.     

Molecular characterization of expressed DQA and DQB genes in the California sea lion ( Zalophus californianus)

To date, there are no published MHC sequences from the California sea lion (Zalophus californianus), a thriving species that, by feeding high on the marine food web, could be a sentinel for disturbances in marine and coastal ecosystems. In this study, degenerate primers and RACE technology were used to amplify near-full-length (MhcZaca- DQB) and full-length (MhcZaca- DQA) expressed class II MHC gene products from the peripheral blood mononuclear cells of two California sea lions in rehabilitation. Five unique Zaca- DQA sequences and eight unique Zaca- DQB sequences, all encoding functional proteins, were identified in the two animals, indicating the presence of multiple DQ- loci in this species. An additional three Zaca- DQB sequences containing features compatible with pseudogenes or null alleles were also identified. Despite the identification of multiple DQA and DQB sequences, the degree of heterogeneity between them was extremely low. To confirm the limited degree of Zaca-DQ nucleotide variation between individuals, we used denaturing gradient gel electrophoresis to examine putative peptide binding region sequences from the peripheral blood leukocyte-derived RNAs of 19 wild-caught California sea lions from physically distinct populations. The pattern of Zaca-DQ sequence migration was identical between individuals and independent of geographical region. This apparent Zaca-DQ sequence identity between sea lions was confirmed by direct sequencing of individual bands. In combination, these findings raise important questions regarding immunogenetic diversity within this thriving species, and should prompt further research into the existence of a highly polymorphic sea lion class II MHC molecule with sequence features that support traditional peptide binding functions.     

Natural killer gene complex (Nkc) allelic variability in inbred mice: evidence for Nkc haplotypes

Allelic variability for mouse Chromosome 6 Nkc loci was assessed in 22 common laboratory strains of mice using selected natural killer gene complex (Nkc)-linked sequence tagged site markers. Most Nkc markers distinguished three or more alleles for a particular locus in the assessed mouse strains. Nkc locus alleles were highly conserved among genealogically related inbred strains, whereas far less similarity was observed among unrelated strains. Concurrent strain-to-strain comparisons for all Nkc-linked loci revealed common and uncommon Nkc haplotypes, including some that were likely recombinant. Nkc allele and haplotype assignments in inbred mouse strains and correlation with phenotypic traits should facilitate positional gene cloning strategies for unknown Nkc-linked trait modification loci.     

Specific amino acid residues in the second hypervariable region of HLA-DQA1 and DQB1 chain genes promote the Ro (SS-A)/La (SS-B) autoantibody responses

In order to define the HLA-DR and DQ alleles, as well as the specific DQA1 and DQB1 chain genes involved in the anti-Ro/La autoantibody responses, RFLP analysis and sequence-specific oligonucleotide typing was carried out on 58 Caucasians and 48 American blacks with SLE or Sj?gren's syndrome and anti-Ro antibodies. Among both Caucasian and black patients, the highest relative risk for the anti-Ro response (both with and without accompanying anti-La) was conferred by heterozygosity for the DQw2.1 (in linkage disequilibrium with HLA-DR3) and DQw6 (a subtype of DQw1) alleles compared with either 269 normal race-matched controls or 80 anti-Ro negative SLE/Sj?gren's syndrome patients. Analysis of individual DQA1 and DQB1 chain alleles revealed that DQA1*0501 and DQB1*0201 were most frequent, followed by DQA1 and DQB1 alleles comprising DQw6. In patients not possessing DQw2.1 and/or DQw6 alleles, HLA-DQB1*0302 and HLA-DQA1*0401 (especially in blacks) were significantly increased. Nucleotide sequence analysis of these associated alleles showed that 100% of patients with anti-Ro had a glutamine residue at position 34 of the outermost domain of the DQA1 chain and/or a leucine at position 26 of the outermost domain of the DQB1 chain. Patients with anti-Ro plus La were more likely to have all four of their DQA1/DQB1 chains containing these amino acid residues than either anti-Ro-negative SLE patients or controls. These data implicate specific amino acid residues on both DQA1 and DQB1 chains located in the floor of the Ag binding cleft of the HLA-DQA1:B1 heterodimer and further suggest a role for "gene dosage" in the anti-Ro (+/- La) autoantibody response.     

Polymorphism at the ovine major histocompatibility complex class II loci

Southern hybridization analysis of the ovine major histocompatibility complex (MHC) ( MhcOvar ) class II region, using sheep-specific probes for the DQA1, DQA2, DQB and DRA loci, has revealed extensive polymorphism. DQA1 and DQAP had eight and 16 alleles respectively, DQB had six and DRA had three alleles. Little information was derived from the DRB locus owing to extensive cross-hybridization between the DRB probe and the DQB locus. Differences in allele frequency between breeds were revealed. At the DQA1 locus a null allele (DQA1-N) was observed with a frequency of between 27% and 45%, making this the most common DQA1 allele in all breeds examined. The frequency of DQA1-N homozygotes was between 11% and 18%, raising questions as to the functional significance of the DQA1 gene. Linkage analysis between the DQA1, DQA2, DQB and DRA loci did not reveal any recombination.     

HLA-AR, an inactivated antigen-presenting locus related to HLA-A. Implications for the evolution of the MHC

The MHC contains many class I genes other than those known to present peptides to T lymphocytes. These additional class I genes vary between species and their functions are unknown. Genes involved in Ag presentation, HLA-A,B,C in humans, are highly diverse whereas other class I genes are of much more limited diversity. We have studied alleles of a gene, HLA-AR, that is closely linked and structurally related to HLA-A; properties consistent with these two loci having been formed by a gene duplication. Compared to HLA-A the diversity in HLA-AR is much less, and does not focus on residues of a putative Ag recognition site. However, the structure of HLA-AR alleles closely resembles those encoding Ag-presenting molecules, although the presence of one or two deleterious mutations prevents these alleles being active in Ag presentation. These results suggest HLA-AR derives from an Ag-presenting locus that became inactivated, possibly as a result of positive natural selection due to changing demands on T cell immunity. Thus absence of diversity may sometimes correlate with loss rather than preservation of function in class I MHC genes.     

The ABCs of MHC

The major histocompatibility complex (MHC) contains the most diverse genes known in vertebrates. These genes encode cell‐surface molecules that play a central role in controlling immunological activity and, as a consequence, in tissue rejection, autoimmunity, and immune responses to infectious diseases. In vertebrates, there are many different MHC genes, most with many alleles. This is true for all primates studied thus far. Multiple loci and alleles allow for an increased peptide‐binding repertoire; their variety has a profound impact on an organism's ability to battle constantly evolving pathogens. The argument that infectious disease is a driving force for MHC variability is supported by observations that most of the allelic variation centers on the amino acid residues that directly interact with foreign peptides. However, while MHC diversity could be maintained through heterozygote advantage, frequency‐dependent selection, or both, the direct evidence that natural selection enhances diversity is limited. Indeed, it is not wholly clear whether selection operates only with respect to disease resistance or if behavioral and biological mechanisms also contribute to the extreme variation that has been observed for many species. Furthermore, reproductive behavior and biology may also help to maintain genetic variability at MHC loci.     

HLA class II linkage disequilibrium and haplotype evolution in the Cayapa Indians of Ecuador.

DNA-based typing of the HLA class II loci in a sample of the Cayapa Indians of Ecuador reveals several lines of evidence that selection has operated to maintain and to diversify the existing level of polymorphism in the class II region. As has been noticed for other Native American groups, the overall level of polymorphism at the DRB1, DQA1, DQB1, and DPB1 loci is reduced relative to that found in other human populations. Nonetheless, the relative evenness in the distribution of allele frequencies at each of the four loci points to the role of balancing selection in the maintenance of the polymorphism. The DQA1 and DQB1 loci, in particular, have near-maximum departures from the neutrality model, which suggests that balancing selection has been especially strong in these cases. Several novel DQA1-DQB1 haplotypes and the discovery of a new DRB1 allele demonstrate an evolutionary tendency favoring the diversification of class II alleles and haplotypes. The recombination interval between the centromeric DPB1 locus and the other class II loci will, in the absence of other forces such as selection, reduce disequilibrium across this region. However, nearly all common alleles were found to be part of DR-DP haplotypes in strong disequilibrium, consistent with the recent action of selection acting on these haplotypes in the Cayapa.     

Evolution of HLA class II molecules: Allelic and amino acid site variability across populations

Analysis of the highly polymorphic beta1 domains of the HLA class II molecules encoded by the DRB1, DQB1, and DPB1 loci reveals contrasting levels of diversity at the allele and amino acid site levels. Statistics of allele frequency distributions, based on Watterson's homozygosity statistic F, reveal distinct evolutionary patterns for these loci in ethnically diverse samples (26 populations for DQB1 and DRB1 and 14 for DPB1). When examined over all populations, the DQB1 locus allelic variation exhibits striking balanced polymorphism (P < 10(-4)), DRB1 shows some evidence of balancing selection (P < 0.06), and while there is overall very little evidence for selection of DPB1 allele frequencies, there is a trend in the direction of balancing selection (P < 0.08). In contrast, at the amino acid level all three loci show strong evidence of balancing selection at some sites. Averaged over polymorphic amino acid sites, DQB1 and DPB1 show similar deviation from neutrality expectations, and both exhibit more balanced polymorphic amino acid sites than DRB1. Across ethnic groups, polymorphisms at many codons show evidence for balancing selection, yet data consistent with directional selection were observed at other codons. Both antigen-binding pocket- and non-pocket-forming amino acid sites show overall deviation from neutrality for all three loci. Only in the case of DRB1 was there a significant difference between pocket- and non-pocket-forming amino acid sites. Our findings indicate that balancing selection at the MHC occurs at the level of polymorphic amino acid residues, and that in many cases this selection is consistent across populations.     

Patterns of Adaptive and Neutral Diversity Identify the Xiaoxiangling Mountains as a Refuge for the Giant Panda

Genetic variation plays a significant role in maintaining the evolutionary potential of a species. Comparing the patterns of adaptive and neutral diversity in extant populations is useful for understanding the local adaptations of a species. In this study, we determined the fine-scale genetic structure of 6 extant populations of the giant panda (Ailuropoda melanoleuca) using mtDNA and DNA fingerprints, and then overlaid adaptive variations in 6 functional Aime-MHC class II genes (DRA, DRB3, DQA1, DQA2, DQB1, and DQB2) on this framework. We found that: (1) analysis of the mtDNA and DNA fingerprint-based networks of the 6 populations identified the independent evolutionary histories of the 2 panda subspecies; (2) the basal (ancestral) branches of the fingerprint-based Sichuan-derived network all originated from the smallest Xiaoxiangling (XXL) population, suggesting the status of a glacial refuge in XXL; (3) the MHC variations among the tested populations showed that the XXL population exhibited extraordinary high levels of MHC diversity in allelic richness, which is consistent with the diversity characteristics of a glacial refuge; (4) the phylogenetic tree showed that the basal clades of giant panda DQB sequences were all occupied by XXL-specific sequences, providing evidence for the ancestor-resembling traits of XXL. Finally, we found that the giant panda had many more DQ alleles than DR alleles (33∶13), contrary to other mammals, and that the XXL refuge showed special characteristics in the DQB loci, with 7 DQB members of 9 XXL-unique alleles. Thus, this study identified XXL as a glacial refuge, specifically harboring the most number of primitive DQB alleles.     

Duplicated DQ haplotypes increase the complexity of restriction element usage in cattle

The MHC of cattle encodes two distinct isotypes of class II molecules, DR and DQ. Unlike humans, cattle lack the DP locus and about half the common haplotypes express duplicated DQ genes. The number and frequency of DQA and DQB alleles means that most cattle are heterozygous. If inter- and/or intrahaplotype pairing of DQA and DQB molecules occurs, cattle carrying DQ-duplicated haplotypes may express more restriction elements than would be predicted by the number of expressed alleles. We are investigating whether duplicated haplotypes cause differences in immune response, particularly in terms of generating protective immunity. We have analyzed the Ag-presenting function of DQ molecules in two heterozygous animals, one of which carries a duplicated haplotype. We compared the class II isotype specificity of T cell clones recognizing a putative vaccinal peptide from foot-and-mouth disease virus (FMDV15). We show for the first time that bovine T cells can recognize Ag in the context of DQ molecules. We also present evidence that interhaplotype pairings of DQA and DQB molecules form functional restriction elements. Both animals showed distinct biases to usage of particular restriction elements. Mainly DQ-restricted clones were derived from the animal with duplicated DQ genes, whereas the majority of clones from the animal with a single DQ gene pair were DR restricted. Furthermore, haplotype bias was observed with both animals. These experiments show that understanding of class II chain pairing in addition to knowledge of the genotype may be important in vaccine design where effective epitope selection is essential.