Population biology of antigen presentation by MHC class I molecules

Disentangling the effects of recombination, selection, and demography on the genetic variation at a major histocompatibility complex class II gene in the alpine chamois

The major histocompatibility complex (MHC) harbours some of the most polymorphic loci in vertebrate genomes. MHC genes are thought to be subject to some form of balancing selection, most likely pathogen-mediated selection. Hence, MHC genes are excellent candidates for exploring adaptive processes. In this study, we investigated the genetic variation at exon 2 of the DRB class II MHC locus in 191 alpine chamois (Rupicapra rupicapra) from 10 populations in the eastern Alps of Italy. In particular, we were interested in distinguishing and estimating the relative impact of selective and demographic factors, while taking into account the confounding effect of recombination. The extremely high d(n)/d(s) ratio and the presence of trans-species polymorphisms suggest that a strong long-term balancing selection effect has been operating at this locus throughout the evolutionary history of this species. We analysed patterns of genetic variation within and between populations, and the mitochondrial D-loop polymorphism patterns were analysed to provide a baseline indicator of the effects of demographic processes. These analyses showed that (i) the chamois experienced a demographic decline in the last 5000-30 000 years, most likely related to the postglacial elevation in temperature; (ii) this demographic process can explain the results of neutrality tests applied to MHC variation within populations, but cannot justify the much weaker divergence between populations implied by MHC as opposed to mitochondrial DNA; (iii) similar sets of divergent alleles are probably maintained with similar frequencies by balancing selection in different populations, and this mechanism is also operating in small isolated populations, which are strongly affected by drift.

Population genetic models of duplicated genes

Various population genetic models of duplicated genes are introduced. The problems covered in this review include the fixation process of a duplicated copy, copy number polymorphism, the fates of duplicated genes and single nucleotide polymorphism in duplicated genes. Because of increasing evidence for concerted evolution by gene conversion, this review introduces recently developed gene conversion models. In the first half, models assuming independent evolution of duplicated genes are introduced, and then the effect of gene conversion is considered in the second half.

HLA-A, -B, -C, and -DRB1 allele and haplotype frequencies distinguish Eastern European Americans from the general European American population

Sequence-based typing was used to identify human leukocyte antigen (HLA)-A, -B, -C, and -DRB1 alleles from 558 consecutively recruited US volunteers with Eastern European ancestry for an unrelated hematopoietic stem cell registry. Four of 31 HLA-A alleles, 29 HLA-C alleles, 59 HLA-B alleles, and 42 HLA-DRB1 alleles identified (A*0325, B*440204, Cw*0332, and *0732N) are novel. The HLA-A*02010101g allele was observed at a frequency of 0.28. Two-, three-, and four-locus haplotypes were estimated using the expectation-maximization algorithm. The highest frequency extended haplotypes (A*010101g-Cw*070101g-B*0801g-DRB1*0301 and A*03010101g-Cw*0702-B*0702-DRB1*1501) were observed at frequencies of 0.04 and 0.03, respectively. Linkage disequilibrium values (Dij') of the constituent two-locus haplotypes were highly significant for both extended haplotypes (P values were less than 8 x 10(-10)) but were consistently higher for the more frequent haplotype. Balancing selection was inferred to be acting on all the four loci, with the strongest evidence of balancing selection observed for the HLA-C locus. Comparisons of the A-C-B haplotypes and DRB1 frequencies in this population with those for African, European, and western Asian populations showed high degrees of identity with Czech, Polish, and Slovenian populations and significant differences from the general European American population.

Is mate choice in humans MHC-dependent?

In several species, including rodents and fish, it has been shown that the Major Histocompatibility Complex (MHC) influences mating preferences and, in some cases, that this may be mediated by preferences based on body odour. In humans, the picture has been less clear. Several studies have reported a tendency for humans to prefer MHC-dissimilar mates, a sexual selection that would favour the production of MHC-heterozygous offspring, who would be more resistant to pathogens, but these results are unsupported by other studies. Here, we report analyses of genome-wide genotype data (from the HapMap II dataset) and HLA types in African and European American couples to test whether humans tend to choose MHC-dissimilar mates. In order to distinguish MHC-specific effects from genome-wide effects, the pattern of similarity in the MHC region is compared to the pattern in the rest of the genome. African spouses show no significant pattern of similarity/dissimilarity across the MHC region (relatedness coefficient, R = 0.015, p = 0.23), whereas across the genome, they are more similar than random pairs of individuals (genome-wide R = 0.00185, p<10⁻³). We discuss several explanations for these observations, including demographic effects. On the other hand, the sampled European American couples are significantly more MHC-dissimilar than random pairs of individuals (R = −0.043, p = 0.015), and this pattern of dissimilarity is extreme when compared to the rest of the genome, both globally (genome-wide R = −0.00016, p = 0.739) and when broken into windows having the same length and recombination rate as the MHC (only nine genomic regions exhibit a higher level of genetic dissimilarity between spouses than does the MHC). This study thus supports the hypothesis that the MHC influences mate choice in some human populations.

There has been a longstanding hypothesis that selection may have led to mating patterns that encourage heterozygosity at Major Histocompatibility Complex (MHC) loci because of improved immune response to pathogens in the offspring of such matings, and, indeed, this has been observed in several model systems. However, in humans, previous studies regarding the role of the MHC in mate choice or preference, both directly in couples and also indirectly in “sweaty T-shirts” experiments, have reported conflicting results. Here, by using genome-wide genotype data and HLA types in African and European American couples, we test whether humans tend to choose MHC-dissimilar mates. This approach allows us to distinguish MHC-specific effects from genome-wide effects. In the African sample, the patterns at MHC loci is confounded by genome-wide effects, possibly resulting from demographic processes relating to the social organization of this population, and no tendency to choose MHC-dissimilar mates is detected. On the other hand, the sampled European Americans appear to have favoured MHC-dissimilar mates, supporting the hypothesis that MHC influences mate choice in some human populations. Thus, this study suggests that, in some cases, humans may rely on biological factors, in addition to social factors, when choosing a mate.

Preservation of a pseudogene by gene conversion and diversifying selection

Interlocus gene conversion is considered a crucial mechanism for generating novel combinations of polymorphisms in duplicated genes. The importance of gene conversion between duplicated genes has been recognized in the major histocompatibility complex and self-incompatibility genes, which are likely subject to diversifying selection. To theoretically understand the potential role of gene conversion in such situations, forward simulations are performed in various two-locus models. The results show that gene conversion could significantly increase the number of haplotypes when diversifying selection works on both loci. We find that the tract length of gene conversion is an important factor to determine the efficacy of gene conversion: shorter tract lengths can more effectively generate novel haplotypes given the gene conversion rate per site is the same. Similar results are also obtained when one of the duplicated genes is assumed to be a pseudogene. It is suggested that a duplicated gene, even after being silenced, will contribute to increasing the variability in the other locus through gene conversion. Consequently, the fixation probability and longevity of duplicated genes increase under the presence of gene conversion. On the basis of these findings, we propose a new scenario for the preservation of a duplicated gene: when the original donor gene is under diversifying selection, a duplicated copy can be preserved by gene conversion even after it is pseudogenized.

Maternal HLA homozygosity and mother-child HLA concordance increase the risk of vertical transmission of HIV-1

BACKGROUND

Mother-child human leukocyte antigen (HLA) concordance and maternal HLA homozygosity may increase the risk of vertical transmission of human immunodeficiency virus type 1 (HIV-1) risk by reducing infant immune responses.

METHODS

We analyzed mother-child HLA concordance and maternal HLA homozygosity in a Kenyan perinatal cohort receiving antenatal zidovudine. HLA concordance was scored as the number of shared class I alleles, and relative risk estimates were adjusted for maternal HIV-1 load.

RESULTS

Among 277 mother-infant pairs, HIV-1 transmission occurred in 58 infants (21%), with in utero transmission in 21 (36%), peripartum transmission in 26 (45%), and transmission via breast-feeding in 11 (19%). With increased concordance, we observed a significant increase in the risk of transmission overall (adjusted hazard ratio [aHR], 1.3 [95% confidence interval {CI}, 1.0-1.7]; P = .04 in utero (adjusted odds ratio, 1.72 [95% CI, 1.0-1.7]; P = .04), and via breast-feeding (aHR, 1.6 [95% CI, 1.0-2.5]; P = .04). Women with homozygosity had higher plasma HIV-1 RNA levels at 32 weeks of gestation (5.1 vs. 4.8 log(10) copies/mL; P = .03) and an increased risk of transmission overall (aHR, 1.7 [95% CI, 1.1-2.7]; P = .03) and via breast-feeding (aHR, 5.8 [95% CI, 1.9-17.7]; P = .002).

CONCLUSION

The risks of overall, in utero, and breast milk HIV-1 transmission increased with HLA concordance and homozygosity. The increased risk may be due to reduced alloimmunity or less diverse protective immune responses.

Selection, diversity and evolutionary patterns of the MHC class II DAB in free-ranging Neotropical marsupials

BACKGROUND

Research on the genetic architecture and diversity of the MHC has focused mainly on eutherian mammals, birds and fish. So far, studies on model marsupials used in laboratory investigations indicated very little or even no variation in MHC class II genes. However, natural levels of diversity and selection are unknown in marsupials as studies on wild populations are virtually absent. We used two endemic South American mouse opossums, Gracilinanus microtarsus and Marmosops incanus, to investigate characteristic features of MHC selection. This study is the first investigation of MHC selection in free-ranging Neotropical marsupials. In addition, the evolutionary history of MHC lineages within the group of marsupials was examined.

RESULTS

G. microtarsus showed extensive levels of MHC diversity within and among individuals as 47 MHC-DAB alleles and high levels of sequence divergence were detected at a minimum of four loci. Positively selected codon sites were identified, of which most were congruent with human antigen binding sites. The diversity in M. incanus was rather low with only eight observed alleles at presumably two loci. However, these alleles also revealed high sequence divergence. Again, positive selection was identified on specific codon sites, all congruent with human ABS and with positively selected sites observed in G. microtarsus. In a phylogenetic comparison alleles of M. incanus interspersed widely within alleles of G. microtarsus with four alleles being present in both species.

CONCLUSION

Our investigations revealed extensive MHC class II polymorphism in a natural marsupial population, contrary to previous assumptions. Furthermore, our study confirms for the first time in marsupials the presence of three characteristic features common at MHC loci of eutherian mammals, birds and fish: large allelic sequence divergence, positive selection on specific sites and trans-specific polymorphism.

Pinpointing a selective sweep to the chimpanzee MHC class I region by comparative genomics

Chimpanzees experienced a reduction of the allelic repertoire at the major histocompatibility complex (MHC) class I A and B loci, which may have been caused by a retrovirus belonging to the simian immunodeficiency virus (SIV) family. Extended MHC haplotypes were defined in a pedigreed chimpanzee colony. Comparison of genetic variation at microsatellite markers mapping inside and outside the Mhc region was carried out in humans and chimpanzees to investigate the genomic extent of the repertoire reduction. Multilocus demographic analyses underscored that chimpanzees indeed experienced a selective sweep that mainly targeted the chromosomal segment carrying the Mhc class I region. Probably due to genetic linkage, the sweep also affected other polymorphic loci, mapping in the close vicinity of the Mhc class I region genes. Nevertheless, although the allelic repertoire at particular Mhc class I and II loci appears to be limited, naturally occurring recombination events allowed the establishment of haplotype diversity after the sweep. However, recombination did not have sufficient time to erase the signal of the selective sweep.

Evaluating linkage disequilibrium and recombination provides a haplotype-tagging SNP panel of the major histocompatibility complex in African Americans

The major histocompatibility complex (MHC) (Chromosome 6p21.3) is a dynamic, immune gene-rich region that is associated with multiple diseases. Haplotype-tagging single-nucleotide polymorphism (htSNP) panels for the MHC can aid association studies but have only been reported for African, Asian and Caucasian populations to date. We genotyped 2154 SNPs spanning a 3.8-Mb region of the classical MHC in 94 healthy African Americans using Illumina BeadArray technology. We describe the haplotype structure of the MHC in African Americans, calculate the recombination rate (0.35 cM Mb(-1)) across the region, identify recombination hot spots and develop a panel of htSNPs for future genetic association studies in this population. We conclude that while patterns of LD and recombination are similar within the MHC to that reported in other populations, differences in minor allele frequency at specific markers necessitates an htSNP panel unique to African Americans, which we provide here for use in future genetic association studies.

Selective downregulation of rhesus macaque and sooty mangabey major histocompatibility complex class I molecules by Nef alleles of simian immunodeficiency virus and human immunodeficiency virus type 2

Human immunodeficiency virus type 1 (HIV-1) Nef downregulates HLA-A and -B molecules, but not HLA-C or -E molecules, based on amino acid differences in their cytoplasmic domains to simultaneously evade cytotoxic T lymphocyte (CTL) and natural killer cell surveillance. Rhesus macaques and sooty mangabeys express orthologues of HLA-A, -B, and -E, but not HLA-C, and many of these molecules have unique amino acid differences in their cytoplasmic tails. We found that these differences also resulted in differential downregulation by primary simian immunodeficiency virus (SIV) SIV(smm/mac) and HIV-2 Nef alleles. Thus, selective major histocompatibility complex class I downregulation is a conserved mechanism of immune evasion for pathogenic SIV infection of rhesus macaques and nonpathogenic SIV infection of sooty mangabeys.

Reanalysis of sequence-based HLA-A, -B and -Cw typings: how ambiguous is today?s SBT typing tomorrow

The permanently increasing number of human leukocyte antigen (HLA)-alleles and the growing list of ambiguities require continuous updating of high-resolution HLA typing results. Two different kinds of ambiguities exist: the first, when two or more allele combinations have identical heterozygous sequences, and the second, when differences are located outside the analyzed region. The number of HLA-A, B and C alleles recognized in 1999 was almost tripled in 2006. Two hundred individuals, sequence-based typing (SBT) typed in the period from 1999 to 2002, were reanalyzed using the 2006 database. A final allele typing result of at least four digits was obtained for HLA-A, -B and -C by heterozygous sequencing of exons 2 and 3 and, if necessary, additional exons and/or allele-specific sequencing. Storage of the individual sequences in a specially developed database enabled reanalysis with all present and future HLA releases. In the 5-year period HLA-A, -B and -C typing results became ambiguous in 37%, 46% and 41% of the cases. Most were because of differences outside the analyzed region; ambiguities because of different allele combinations with identical heterozygous sequences were present in 7%, 8% and 13% of the HLA-A, -B and -C typings. These results indicate that within 5 years, approximately half of the HLA SBT typings become ambiguous.

HLA genes of Aleutian Islanders living between Alaska (USA) and Kamchatka (Russia) suggest a possible southern Siberia origin

Aleuts HLA profile has been compared with that of neighboring and worldwide populations. Thirteen thousand one hundred and sixty-four chromosomes have been used for this study. Computer programs have obtained HLA allele frequencies, genetic distances between populations, NJ relatedness dendrograms, correspondence analysis and most frequent HLA extended haplotypes. Aleuts have inhabited Aleutian Islands since about 9000 years BP according to fossil and genetic (mtDNA) records. They are genetically different to Eskimo, Amerindian and Na-Dene speakers according to their HLA profile; this correlates with cultural and anthropological Aleut distinctiveness. No typical Amerindian HLA alleles have been found in a significant frequency. Their HLA relatedness to Saami (or Lapps, northern Scandinavians), Finns and Pomors (North-West Russia) indicates an ancient possible origin from the Baikal Lake Area (southern Siberia) around the present day Buryat peopling area; other origins are not discarded. Aleuts characteristic HLA profile may influence future transplantation programs in the region and be useful to study diseases linked to HLA epidemiology.

Phase-dependent immune evasion of herpesviruses

Viruses employ various modes to evade immune detection. Two possible evasion modes are a reduction of the number of epitopes presented and the mimicry of host epitopes. The immune evasion efforts are not uniform among viral proteins. The number of epitopes in a given viral protein and the similarity of the epitopes to host peptides can be used as a measure of the viral attempts to hide this protein. Using bioinformatics tools, we here present a genomic analysis of the attempts of four human herpesviruses (herpes simplex virus type 1-human herpesvirus 1, Epstein-Barr virus-human herpesvirus 4, human cytomegalovirus-human herpesvirus 5, and Kaposi's sarcoma-associated herpesvirus-human herpesvirus 8) and one murine herpesvirus (murine herpesvirus 68) to escape from immune detection. We determined the full repertoire of CD8 T-lymphocyte epitopes presented by each viral protein and show that herpesvirus proteins present many fewer epitopes than expected. Furthermore, the epitopes that are presented are more similar to host epitopes than are random viral epitopes, minimizing the immune response. We defined a score for the size of the immune repertoire (the SIR score) based on the number of epitopes in a protein. The numbers of epitopes in proteins expressed in the latent and early phases of infection were significantly smaller than those in proteins expressed in the lytic phase in all tested viruses. The latent and immediate-early epitopes were also more similar to host epitopes than were lytic epitopes. A clear trend emerged from the analysis. In general, herpesviruses demonstrated an effort to evade immune detection. However, within a given herpesvirus, proteins expressed in phases critical to the fate of infection (e.g., early lytic and latent) evaded immune detection more than all others. The application of the SIR score to specific proteins allows us to quantify the importance of immune evasion and to detect optimal targets for immunotherapy and vaccine development.

MHC adaptive divergence between closely related and sympatric African cichlids

Background

The haplochromine cichlid species assemblages of Lake Malawi and Victoria represent some of the most important study systems in evolutionary biology. Identifying adaptive divergence between closely-related species can provide important insights into the processes that may have contributed to these spectacular radiations. Here, we studied a pair of sympatric Lake Malawi species, Pseudotropheus fainzilberi and P. emmiltos, whose reproductive isolation depends on olfactory communication. We tested the hypothesis that these species have undergone divergent selection at MHC class II genes, which are known to contribute to olfactory-based mate choice in other taxa.

Methodology/Principal Findings

Divergent selection on functional alleles was inferred from the higher genetic divergence at putative antigen binding sites (ABS) amino acid sequences than at putatively neutrally evolving sites at intron 1, exon 2 synonymous sequences and exon 2 amino acid residues outside the putative ABS. In addition, sympatric populations of these fish species differed significantly in communities of eukaryotic parasites.

Conclusions/Significance

We propose that local host-parasite coevolutionary dynamics may have driven adaptive divergence in MHC alleles, influencing odor-mediated mate choice and leading to reproductive isolation. These results provide the first evidence for a novel mechanism of adaptive speciation and the first evidence of adaptive divergence at the MHC in closely related African cichlid fishes.

Sensory genes and mate choice: Evidence that duplications, mutations, and adaptive evolution alter variation in mating cue genes and their receptors

Fascinating new data, revealed through gene sequencing, comparative genomics, and genetic engineering, precisely establish which genes are involved in mate choice and mating activity--behaviors that are surprisingly understudied from a genetic perspective. Discussed here are some of the recently identified visual and chemosensory genes that are involved in mate choice and mating behavior. These genes' products are involved in the production, transmission, and receipt of crucial sensory mate-choice cues that affect fitness. This review exposes newfound evidence that alternative splicing, gene-expression pattern changes, and molecular genetic variation in sensory genes are crucial for both intra- and interspecific mate choice and mating success. Many sensory genes have arisen through gene duplications, and data amassed from studies conducted at scales ranging from individual genes to genomic comparisons show that strong, positive Darwinian selection acts on several mating-related genes and that these genes evolve rapidly.

Major histocompatibility complex heterozygosity reduces fitness in experimentally infected mice

It is often suggested that heterozygosity at major histocompatibility complex (MHC) loci confers enhanced resistance to infectious diseases (heterozygote advantage, HA, hypothesis), and overdominant selection should contribute to the evolution of these highly polymorphic genes. The evidence for the HA hypothesis is mixed and mainly from laboratory studies on inbred congenic mice, leaving the importance of MHC heterozygosity for natural populations unclear. We tested the HA hypothesis by infecting mice, produced by crossbreeding congenic C57BL/10 with wild ones, with different strains of Salmonella, both in laboratory and in large population enclosures. In the laboratory, we found that MHC influenced resistance, despite interacting wild-derived background loci. Surprisingly, resistance was mostly recessive rather than dominant, unlike in most inbred mouse strains, and it was never overdominant. In the enclosures, heterozygotes did not show better resistance, survival, or reproductive success compared to homozygotes. On the contrary, infected heterozygous females produced significantly fewer pups than homozygotes. Our results show that MHC effects are not masked on an outbred genetic background, and that MHC heterozygosity provides no immunological benefits when resistance is recessive, and can actually reduce fitness. These findings challenge the HA hypothesis and emphasize the need for studies on wild, genetically diverse species.

Identification of a novel HLA-B*56 allele, B*5618 and an extension of B*2736 by sequence-based typing

We report the identification of a novel human leukocyte antigen (HLA)-B*56 allele, B*5618 and an extension of B*2736 that were found during routine high-resolution sequence-based typing in Chinese Han individual. The B*5618 allele has 4nt changes from B*5610 in exon 3, The B*2736 allele has 10nt changes from B*270401 in exons 3-4.

The impact of HLA-B micropolymorphism outside primary peptide anchor pockets on the CTL response to CMV

The factors controlling epitope selection in the T cell response to persistent viruses are not fully understood, and we have examined this issue in the context of four HLA-B*35-binding peptides from the pp65 antigen of human cytomegalovirus, two of which are previously undescribed. Striking differences in the hierarchy of immunodominance between these four epitopes were observed in healthy virus carriers expressing HLA-B*3501 versus B*3508, two HLA-B allotypes that differ by a single amino acid at position 156 (HLA-B*3501, (156)Leucine; HLA-B*3508, (156)Arginine) that projects from the alpha2 helix into the centre of the peptide-binding groove. While HLA-B*3501(+) individuals responded most strongly to the (123)IPSINVHHY(131) and (366)HPTFTSQY(373) epitopes, HLA-B*3508(+) individuals responded preferentially to (103)CPSQEPMSIYVY(114) and (188)FPTKDVAL(195). By comparing peptide-MHC association and disassociation rates with peptide immunogenicity, it was clear that dissociation rates correlate more closely with the hierarchy of immunodominance among the four pp65 peptides. These findings demonstrate that MHC micropolymorphism at positions outside the primary anchor residue binding pockets can have a major impact on determinant selection in antiviral T cell responses. Such influences may provide the evolutionary pressure that maintains closely related MHC molecules in diverse human populations.

A novel HLA-A2 variant, A*9203, identified by sequence-based typing

HLA-A*9203, found in Taiwan using sequence-based typing method, was identical to HLA-A*0207 in exon 3 but differed in exon 2 by five nucleotide substitutions at positions 240-282 corresponding to three amino acid changes at codons 62, 66 and 70. Since this substitution motif is also seen in A*11 and A*03, it is likely that a gene conversion event from A*11 or A*03 to a A*0207 backbone may have been the process used in generating HLA-A*9203.

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.

Specific combinations of HLA-DR2 and DR3 class II haplotypes contribute graded risk for disease susceptibility and autoantibodies in human SLE

The human leukocyte antigen (HLA) Class II antigen presentation alleles DR and DQ are associated with susceptibility to systemic lupus erythematosus (SLE) and the production of lupus-related autoantibodies. Here, we explore the effect of different combinations of Class II risk haplotypes in a large, multi-center collection of 780 SLE families. Haplotypes bearing the DRB1*1501/DQB1*0602 (DR2) and DRB1*0301/DQB1*0201 (DR3) alleles were present in nearly two-thirds of SLE cases and were significantly associated with disease susceptibility in both family-based and case-control study designs. DR3-containing haplotypes conferred higher risk for disease than DR2, and individual homozygous for DR3 or compound heterozygous for DR3 and DR2 showed the highest risk profile. DR2 haplotypes were also found to be associated with antibodies to the nuclear antigen Sm, and, as previously observed, DR3 genotypes were associated with Ro and La autoantibodies. Interestingly, SLE cases and unaffected family members who were DR2/DR3 compound heterozygotes showed particularly strong risk of developing antibodies to Ro, and were enriched for La and Sm. These data provide convincing evidence that particular combinations of HLA Class II DR2 and DR3 haplotypes are key determinants of autoantibody production and disease susceptibility in human SLE.

HLA-A, -B, -C, -DRB1 allele and haplotype frequencies in an African American population

Sequence-based typing was used to identify human leukocyte antigen (HLA)-A, -B, -C, and -DRB1 alleles from 564 consecutively recruited African American volunteers for an unrelated hematopoietic stem cell registry. The number of known alleles identified at each locus was 42 for HLA-A, HLA-B 67, HLA-C 33, and HLA-DRB1 44. Six novel alleles (A*260104, A*7411, Cw*0813, Cw*1608, Cw*1704, and DRB1*130502) not observed in the initial sequence-specific oligonucleotide probe testing were characterized. The action of balancing selection, shaping more 'even' than expected allele frequency distributions, was inferred for all four loci and significantly so for the HLA-A and DRB1 loci. Two-, three-, and four-locus haplotypes were estimated using the expectation maximization algorithm. Comparisons with other populations from Africa and Europe suggest that the degree of European admixture in the African American population described here is lower than that in other African American populations previously reported, although HLA-A:B haplotype frequencies similar to those in previous studies of African American individuals were also noted.

DC-SIGN and DC-SIGNR genetic diversity among different ethnic populations: potential implications for pathogen recognition and disease susceptibility

Dendritic cell–specific intracellular adhesion molecule-3–grabbing nonintegrin (DC-SIGN) and DC-SIGNR are C-type lectins that serve both as cell adhesion and pathogen recognition receptors. Because of the essential role of the these molecules in the immune response, the implication of their alleles in human disease states, and the possible genetic variation at these loci among ethnically diverse populations, we undertook a study to analyze the full extent of DC-SIGN and DC-SIGNR polymorphisms in Caucasian Canadian and indigenous African populations. We report several novel nucleotide variants within regulatory 5′- and 3′-untranslated regions of the genes that could affect their transcription and translation. There were significant differences in the distribution of DC-SIGN and DC-SIGNR alleles among African and non-African populations. Finally, our study clearly demonstrates that Africans show greater genetic diversity at these two closely-related immune loci than observed in other major population groups. The differences may reflect evolutionary pressures generated by environmental factors, such as prevalent pathogens in these geographically distinct regions. Further studies will be needed to determine the net impact of DC-SIGN and DC-SIGNR genetic variants on the expression, translation, and function of the proteins and to understand how these functional polymorphisms may affect immune responses or immune escape.

Gene duplication, allelic diversity, selection processes and adaptive value of MHC class II DRB genes of the bank vole, Clethrionomys glareolus

The generation and maintenance of allelic polymorphism in genes of the major histocompatibility complex (MHC) is a central issue in evolutionary genetics. Recently, the focus has changed from ex situ to in situ populations to understand the mechanisms that determine adaptive MHC polymorphism under natural selection. Birth-and-death evolution and gene conversion events are considered to generate sequence diversity in MHC genes, which subsequently is maintained by balancing selection through parasites. The ongoing arms race between the host and parasites leads to an adaptive selection pressure upon the MHC, evident in high rates of non-synonymous vs synonymous substitution rates. We characterised the MHC class II DRB exon 2 of free living bank voles, Clethrionomys glareolus by single-strand conformation polymorphism and direct sequencing. Unlike other arvicolid species, the DRB locus of the bank vole is at least quadruplicated. No evidence for gene conversion events in the Clgl-DRB sequences was observed. We found not only high allelic polymorphism with 26 alleles in 36 individuals but also high rates of silent polymorphism. Exceptional for MHC class II genes is a purifying selection pressure upon the majority of MHC-DRB sequences. Further, we analysed the association between certain DRB alleles and the parasite burden with gastrointestinal trichostrongyle nematodes Heligmosomum mixtum and Heligmosomoides glareoli and found significant quality differences between specific alleles with respect to infection intensity. Our findings suggest a snapshot in an evolutionary process of ongoing birth-and-death evolution. One allele cluster has lost its function and is already silenced, another is loosing its adaptive value in terms of gastrointestinal nematode resistance, while a third group of alleles indicates all signs of classical functional MHC alleles.

Direct class I HLA antigen discovery to distinguish virus-infected and cancerous cells

Class I human leukocyte antigen molecules are nature's proteome-scanning chips, presenting thousands of endogenously loaded peptides on the surface of virtually every cell in the body. Cytotoxic T cells survey the class I human leukocyte antigen peptide cargo presented, recognize peptides unique to unhealthy cells and destroy diseased cells. A precise understanding of how class I molecules distinguish diseased cells is positioned to drive immune-based diagnostics, therapies and vaccines. When identifying epitopes unique to unhealthy cells, the most experimentally direct approach is to examine the class I-presented peptides of infected/cancerous cells. Here we discuss the strategies adapted for protein production, protein/peptide purification, peptide separation and for maintaining experimental reproducibility during the direct characterization of class I human leukocyte antigen peptides.

Genetic basis for MHC-dependent mate choice

Genes in the major histocompatibility complex (MHC), best known for their role in immune recognition and transplantation success, are also involved in modulating mate choice in mice. Early studies with inbred, congenic mouse lines showed that mate choice tended to favor nonself MHC types. A similar phenomenon was demonstrated with semi-wild mice as well. Subsequent studies showed that, rather than nonself choices, it was more accurate to say that mice chose nonparental MHC types for mates since preferences for nonself could be reversed if mice were fostered from birth on parents with nonself MHC types. Other studies have demonstrated that parent-offspring recognition is also regulated by MHC-determined signals suggesting that this system is one of general importance for mouse behavior. Many studies have now demonstrated that volatile mouse body odors are regulated by MHC genes and it is presumably these odor differences that underlie mate choice and familial recognition. Recent studies have shown that many odorants are controlled by the MHC but the mechanism by which MHC genes exert their influence has not been identified. Surprisingly, not only are volatile body odors influenced by MHC genes but so too are nonvolatile signals. Peptides bound to the MHC protein may also function in individual recognition. The extent to which this system is involved in mate choice of other species is unclear although there are some suggestive studies. Indeed, there is tentative evidence that MHC differences, presumably acting via odor changes, may influence human partner selection. Further studies should clarify both the mechanism underlying MHC influence on body odors as well as the generality of their importance in mate selection.

Beta2-microglobulin promotes the growth of human renal cell carcinoma through the activation of the protein kinase A, cyclic AMP-responsive element-binding protein, and vascular endothelial growth factor axis

PURPOSE

Beta(2)-microglobulin (beta2M), a soluble protein secreted by cancer and host inflammatory cells, has various biological functions, including antigen presentation. Because aberrant expression of beta2M has been reported in human renal cell carcinoma, we investigated the effects of beta2M overexpression on cancer cell growth and analyzed its molecular signaling pathway.

EXPERIMENTAL DESIGN

We established clonal cell lines that overexpressed beta2M in human renal cell carcinoma (SN12C) cells and then examined cell growth in vitro and in vivo and studied the beta2M-mediated downstream cell signaling pathway.

RESULTS

Our results showed that beta2M expression positively correlates with (a) in vitro growth on plastic dishes and as Matrigel colonies, (b) cell invasion and migration in Boyden chambers, and (c) vascular endothelial growth factor (VEGF) expression and secretion by cells. We found, in addition, that beta2M mediates its action through increased phosphorylation of cyclic AMP-responsive element-binding protein (CREB) via the protein kinase A-CREB axis, resulting in increased VEGF expression and secretion. In convergence with this signal axis, beta2M overexpression also activated both phosphatidylinositol 3-kinase/Akt and mitogen-activated protein kinase pathways. Beta2M overexpression induced accelerated growth of SN12C in mouse subcutis and bone. Interrupting the beta2M signaling pathway using small interfering RNA led to apoptosis with increased activation of caspase-3 and caspase-9 and cleaved poly(ADP-ribose) polymerase.

CONCLUSIONS

Our results showed for the first time that the beta2M-protein kinase A-CREB-VEGF signaling axis plays a crucial role in support of renal cell carcinoma growth and progression and reveals a novel therapeutic target.

Balancing selection and the evolution of functional polymorphism in Old World monkey TRIM5alpha

Retroviral restriction factor TRIM5alpha exhibits a high degree of sequence variation among primate species. It has been proposed that this diversity is the cumulative result of ancient, lineage-specific episodes of positive selection. Here, we describe the contribution of within-species variation to the evolution of TRIM5alpha. Sampling within two geographically distinct Old World monkey species revealed extensive polymorphism, including individual polymorphisms that predate speciation (shared polymorphism). In some instances, alleles were more closely related to orthologues of other species than to one another. Both silent and nonsynonymous changes clustered in two domains. Functional assays revealed consequences of polymorphism, including differential restriction of a small panel of retroviruses by very similar alleles. Together, these features indicate that the primate TRIM5alpha locus has evolved under balancing selection. Except for the MHC there are few, if any, examples of long-term balancing selection in primates. Our results suggest a complex evolutionary scenario, in which fixation of lineage-specific adaptations is superimposed on a subset of critical polymorphisms that predate speciation events and have been maintained by balancing selection for millions of years.

Receptor-ligand analyses define minimal killer cell Ig-like receptor (KIR) in humans

Interactions between inhibitory killer cell immunoglobulin-like receptors (iKIR) and human leukocyte antigen (HLA) class I molecules regulate natural killer (NK) cell responses to eliminate infected and transformed cells while maintaining tolerance to healthy cells. Unlinked polymorphic gene families encode KIR receptors and HLA class I ligands and their independent segregation results in a variable number and type of iKIR + HLA pairs inherited in individuals. The diversity in the co-inheritance of iKIR + HLA pairs and activating KIR (aKIR) genes in 759 unrelated individuals from four ethnic populations was analyzed. Every individual studied inherited a minimum of one iKIR + HLA pair; suggesting that major histocompatibility complex class I-dependent inhibitory KIR signaling is essential for human NK cell function. In contrast, 13.4% of the study group lacked all aKIR genes. Twenty percent of the study group carried only one of the four iKIR + HLA pairs. Interestingly, 3% of the study group carrying only KIR2DL3 + HLA-C1 as an iKIR + HLA pair lacked aKIR genes. These data suggest that a single iKIR can constitute the minimal KIR repertoire for human NK cells. Genotypes carrying an equal number of iKIR + HLA pairs and aKIR genes represented 20% of the study group. The remaining individuals had either a dominant inhibitory KIR genotype (iKIR + HLA > aKIR) or a dominant activating KIR genotype (iKIR + HLA < aKIR). Genotypes encoding these imbalanced inhibitory and activating interactions may contribute to susceptibility or resistance to human diseases.

Major histocompatibility complex class IIB allele polymorphism and its association with resistance/susceptibility to Vibrio anguillarum in Japanese flounder (Paralichthys olivaceus)

The full length of major histocompatibility complex (MHC) class IIB cDNA was cloned from a Chinese population of Paralichthys olivaceus by homology cloning and rapid amplification of cDNA ends-polymerase chain reaction (RACE-PCR). The MHC IIB genomic sequence is 1,864 bp long and consists of 34-bp 5'UTR, 741-bp open reading frame, 407-bp 3'UTR, 96-bp intron1, 392-bp intron2, 85-bp intron3, and 109-bp intron4. Phylogenetic analysis showed that the putative MHC class IIB amino acid of the Chinese P. olivaceus shared 28.3% to 85.4% identity with that of the reported MHC class IIB in other species. A significant association between MHC IIB polymorphism and disease resistance/susceptibility was found in Chinese P. olivaceus. Thirteen different MHC IIB alleles were identified among 411 clones from 84 individuals. Among the 280 (268) nucleotides, 32 (11.4%) nucleotide positions were variable. Most alleles such as alleles a, b, c, d, e, f, j, k, i, m were commonly found in both resistant and susceptible stock. Via chi2 test, allele d was significantly more prevalent in individuals from susceptible stock than from resistant stock, and their percentages were 23.80% and 7.14%, respectively. In addition, allele g occurred in 9 and allele h in 4 of 42 resistant individuals that were not present in the susceptible stock; their percentages were 21.4% and 9.52%, respectively. Although allele l was found only in 8 individuals from the susceptible stock, its percentage is 19.05%.

KIR3DL3 allelic diversity: six new alleles exhibit both conservative and non-conservative substitutions

KIR3DL3 alleles were characterized in two families and one unrelated individual. Based on exon 2-9 nucleotide sequences, six novel alleles, 3DL3*00402, *005, *006, *007, *00801, *00802, were identified bringing the total number of known alleles to 11. Compared with 3DL3*001, the six new alleles differ by from three to nine nucleotides and from three to four amino acids. The new alleles double the number of known polymorphic positions to 18 with variation in exons encoding the extracellular domains, transmembrane region, and a portion of the cytoplasmic tail. Many of the nucleotide substitutions are shared among alleles of 3DL3 or other KIR loci, but five were found only in single 3DL3 alleles. Comparison of intron sequences among individuals carrying the same allele showed a modest number of substitutions with the exception of 3DL3*001 which differed substantially in its intron sequences. Two alleles sharing coding region sequences, 3DL3*00201 and 3DL3*00202, were also substantially different in intron sequences.

Generation and maintenance of diversity in the cattle MHC class I region

Major histocompatibility complex (MHC) class I genes play a crucial role in the immune defence against intracellular pathogens. An important evolutionary strategy is to generate and maintain a high level of diversity in these genes. Humans express three highly polymorphic classical MHC class I genes (HLA-A, HLA-B and HLA-C). In contrast, some species, for example rat and rhesus macaque, maintain diversity by generation of haplotypes that vary considerably with regard to the number and combination of transcribed genes. Cattle appear to use both strategies. We show that various combinations of six apparently classical genes, three of which are highly polymorphic, are transcribed on different haplotypes. Although additional sequences were identified in both cDNA and gDNA, it was not possible to assign them to any of these defined genes. Most were highly divergent or were non-classical class I genes. Thus, we found little evidence for frequent duplication and deletion of classical class I genes as reported in some other species. However, the maintenance of class I diversity in cattle may involve limited gene shuffling and deletion, possibly as a result of unequal crossing-over within the class I region.

The human brain is a detector of chemosensorily transmitted HLA-class I-similarity in same- and opposite-sex relations

Studies on subjective body odour ratings suggest that humans exhibit preferences for human leucocyte antigen (HLA)-dissimilar persons. However, with regard to the extreme polymorphism of the HLA gene loci, the behavioural impact of the proposed HLA-related attracting signals seems to be minimal. Furthermore, the role of HLA-related chemosignals in same- and opposite-sex relations in humans has not been specified so far. Here, we investigate subjective preferences and brain evoked responses to body odours in males and females as a function of HLA similarity between odour donor and smeller. We show that pre-attentive processing of body odours of HLA-similar donors is faster and that late evaluative processing of these chemosignals activates more neuronal resources than the processing of body odours of HLA-dissimilar donors. In same-sex smelling conditions, HLA-associated brain responses show a different local distribution in male (frontal) and female subjects (parietal). The electrophysiological results are supported by significant correlations between the odour ratings and the amplitudes of the brain potentials. We conclude that odours of HLA-similar persons function as important social warning signals in inter- and intrasexual human relations. Such HLA-related chemosignals may contribute to female and male mate choice as well as to male competitive behaviour.

Polymorphism, natural selection, and structural modeling of class Ia MHC in the African clawed frog (Xenopus laevis)

In the African clawed frog (Xenopus laevis), two deeply divergent allelic lineages of multiple genes of the class I MHC region have been discovered. For the MHC class I UAA locus, functional differences and the molecular basis for lineages maintenance are unknown. Alleles of linked class I region genes also exhibit strong disequilibrium with specific MHC alleles, but the underlying cause is not clear. We use MHC class Ia sequence data to estimate substitution rates and investigate structural differences between allelic lineages from protein models. Results indicate the operation of natural selection, and differences in the steric properties in the F pocket of the peptide-binding region among lineages. Variability in this pocket likely enables allelic lineages to bind very different sets of peptides and to interact differently with MHC chaperones in the endoplasmic reticulum. These results constitute evidence of the molecular evolutionary basis for 1) the maintenance of allelic lineages, 2) functional differences among lineages, and 3) strong linkage disequilibrium of allelic variants of class I region genes in X. laevis.

Comparative genetics of MHC polymorphisms in different primate species: duplications and deletions

Gene products of the major histocompatibility complex (MHC) play a crucial role in the activation of adaptive (antigen-dependent) immune responses. In this paper similarities and dissimilarities among the MHCs of different primate species and their functional implications are reviewed. The human HLA system represents the most thoroughly investigated MHC of any contemporary living primate species, and so it will serve as a reference.

Diversity and duplication of DQB and DRB-like genes of the MHC in baleen whales (suborder: Mysticeti)

The molecular diversity and phylogenetic relationships of two class II genes of the baleen whale major histocompatibility complex were investigated and compared to toothed whales and out-groups. Amplification of the DQB exon 2 provided sequences showing high within-species and between-species nucleotide diversity and uninterrupted reading frames consistent with functional class II loci found in related mammals (e.g., ruminants). Cloning of amplified products indicated gene duplication in the humpback whale and triplication in the southern right whale, with average nucleotide diversity of 5.9 and 6.3%, respectively, for alleles of each species. Significantly higher nonsynonymous divergence at sites coding for peptide binding (32% for humpback and 40% for southern right) suggested that these loci were subject to positive (overdominant) selection. A population survey of humpback whales detected 23 alleles, differing by up to 21% of their inferred amino acid sequences. Amplification of the DRB exon 2 resulted in two groups of sequences. One was most similar to the DRB3 of the cow and present in all whales screened to date, including toothed whales. The second was most similar to the DRB2 of the cow and was found only in the bowhead and right whales. Both loci showed low diversity among species and apparent loss of function or altered function including interruption of reading frames. Finally, comparison of inferred protein sequence of the DRB3-like locus suggested convergence with the DQB, perhaps resulting from intergenic conversion or recombination.

Point mutations with positive selection were a major force during the evolution of a receptor-kinase resistance gene family of rice

The rice (Oryza sativa) Xa26 gene, which confers resistance to bacterial blight disease and encodes a leucine-rich repeat (LRR) receptor kinase, resides at a locus clustered with tandem homologous genes. To investigate the evolution of this family, four haplotypes from the two subspecies of rice, indica and japonica, were analyzed. Comparative sequence analysis of 34 genes of 10 types of paralogs of the family revealed haplotype polymorphisms and pronounced paralog diversity. The orthologs in different haplotypes were more similar than the paralogs in the same haplotype. At least five types of paralogs were formed before the separation of indica and japonica subspecies. Only 7% of amino acid sites were detected to be under positive selection, which occurred in the extracytoplasmic domain. Approximately 74% of the positively selected sites were solvent-exposed amino acid residues of the LRR domain that have been proposed to be involved in pathogen recognition, and 73% of the hypervariable sites detected in the LRR domain were subject to positive selection. The family is formed by tandem duplication followed by diversification through recombination, deletion, and point mutation. Most variation among genes in the family is caused by point mutations and positive selection.

Cytotoxic T cells

The immune system is a complex arrangement of cells and molecules that preserve the integrity of the organism by elimination of all elements judged dangerous. Within the immune system, a humoral and a cellular as well as an innate and an adaptive arm can be differentiated. The key players of adaptive cellular immune responses are T lymphocytes in general and, for the effector function, cytotoxic T lymphocytes (CTLs) in particular. T lymphocytes arise in the bone marrow and migrate to the thymus for maturation. During this process, T cells somatically rearrange gene segments, eventually leading to the expression of a unique antigen-binding molecule, the T-cell receptor (TCR). This receptor allows them to monitor all cells of the body, ready to destroy any cell posing a threat to the organism. Cytotoxicity is exerted directly through the Fas or perforin pathway and/or indirectly by the release of cytokines. Obviously, the activity of such a potent cell is tightly regulated. Indeed, a predominance of stimulatory over inhibitory signals is required for effective immune responses to pathogens, and a predominance of inhibitory over stimulatory signals is required for maintenance of self-tolerance. Still, several situations occur in which an inappropriate CTL response leads to either autoimmune disease or persistence of pathogens.

A combined immuno-informatics and structure-based modeling approach for prediction of T cell epitopes of secretory proteins of Mycobacterium tuberculosis

The role of secretory proteins of Mycobacterium tuberculosis in pathogenesis and stimulation of specific host responses is well documented. They are also shown to activate different cell types, which subsequently present mycobacterial antigens to T cells. Therefore identification of T cell epitopes from this set of proteins may serve to define candidate antigens with vaccine potential. Fifty-two secretory proteins of M. tuberculosis H37Rv were analyzed computationally for the presence of HLA class I binding nonameric peptides. All possible overlapping nonameric peptide sequences from 52 secretory proteins were generated in silico and analyzed for their ability to bind to 33 alleles belonging to A, B and C loci of HLA class I. Fifteen percent of generated peptides are predicted to bind to HLA with halftime of dissociation T(1/2) >or=100 min and 73% of the peptides predicted to bind are mono-allelic in their binding. The structural basis for recognition of no-namers by different HLA molecules was studied employing structural modeling of HLA class I-peptide complexes and there exists a good correlation between structural analysis and binding prediction. Pathogen peptides that could behave as self- or partially self-peptides in the host were eliminated using a comparative study with the human proteome, thus reducing the number of peptides for analysis. The implications of the finding for vaccine development are discussed vis-à-vis the limitations of the use of subunit vaccine and DNA vaccine.

The immunogenicity of a viral cytotoxic T cell epitope is controlled by its MHC-bound conformation

Thousands of potentially antigenic peptides are encoded by an infecting pathogen; however, only a small proportion induce measurable CD8(+) T cell responses. To investigate the factors that control peptide immunogenicity, we have examined the cytotoxic T lymphocyte (CTL) response to a previously undefined epitope ((77)APQPAPENAY(86)) from the BZLF1 protein of Epstein-Barr virus (EBV). This peptide binds well to two human histocompatibility leukocyte antigen (HLA) allotypes, HLA-B*3501 and HLA-B*3508, which differ by a single amino acid at position 156 ((156)Leucine vs. (156)Arginine, respectively). Surprisingly, only individuals expressing HLA-B*3508 show evidence of a CTL response to the (77)APQPAPENAY(86) epitope even though EBV-infected cells expressing HLA-B*3501 process and present similar amounts of peptide for CTL recognition, suggesting that factors other than peptide presentation levels are influencing immunogenicity. Functional and structural analysis revealed marked conformational differences in the peptide, when bound to each HLA-B35 allotype, that are dictated by the polymorphic HLA residue 156 and that directly affected T cell receptor recognition. These data indicate that the immunogenicity of an antigenic peptide is influenced not only by how well the peptide binds to major histocompatibility complex (MHC) molecules but also by its bound conformation. It also illustrates a novel mechanism through which MHC polymorphism can further diversify the immune response to infecting pathogens.