Natural selection at major histocompatibility complex loci of vertebrates

Evidence of balancing selection at the HLA-G promoter region

HLA-G is a class Ib HLA gene with unique tissue expression pattern and immunomodulatory properties. Polymorphisms in the HLA-G promoter region have been associated with miscarriage and asthma, whereas expression levels have been associated with a wide range of pathologic conditions as well as survival of embryos after in vitro fertilization and of organs after transplantation. Here, we characterize the sequence variation and haplotype structure of the HLA-G promoter and flanking sequences in 44 African Americans, 47 European Americans and 43 Han Chinese by haplotype-specific PCR and sequencing. In all three populations, we observed high levels of nucleotide variation, an excess of intermediate-frequency alleles, and a genealogy with two common haplotypes separated by deep branches, features that are suggestive of balancing selection acting in this region. Comparisons to HLA-A and a pseudogene, HLA-J, suggested that the observed pattern of sequence variation in the HLA-G promoter region is not likely due to other selected HLA genes. We suggest that the mechanism for this selection is related to the highly regulated expression pattern of HLA-G and that high- and low-expressing promoters may be favored under temporally and/or spatially varying selective pressures.

Associations between malaria and MHC genes in a migratory songbird

Malaria parasites are a widespread and species-rich group infecting many wild populations of mammals, birds and reptiles. Studies on humans have demonstrated that genetic factors play a key role in the susceptibility and outcome of malaria infections. Until the present study, it has not been examined whether genetic variation in hosts is important for the outcome of malaria infections in natural avian populations. We investigated associations between major histocompatibility complex (MHC) genes and prevalence of three different avian malaria parasites (Haemoproteus payevskyi (GRW1), Plasmodium sp. (GRW2) and Plasmodium sp. (GRW4)) in a long-term study of great reed warblers Acrocephalus arundinaceus. We hypothesized that the MHC genes could either give full protection against a malaria infection, or confer protection against lethal malaria and direct the infection towards being milder. We found a positive association between numbers of MHC class I alleles (a measure of level of heterozygosity) and prevalence of the GRW2 parasite, suggesting the latter scenario. There was also a positive association between a specific MHC allele (B4b), previously shown to be under frequency-dependent selection in the study population, and prevalence of GRW2. These associations suggest that individuals carrying either a large number of MHC alleles or a specific MHC allele are protected against lethal malaria infections.

Pathogen-driven selection and worldwide HLA class I diversity

The human leukocyte antigen (HLA; known as MHC in other vertebrates) plays a central role in the recognition and presentation of antigens to the immune system and represents the most polymorphic gene cluster in the human genome [1]. Pathogen-driven balancing selection (PDBS) has been previously hypothesized to explain the remarkable polymorphism in the HLA complex, but there is, as yet, no direct support for this hypothesis [2 and 3]. A straightforward prediction coming out of the PDBS hypothesis is that populations from areas with high pathogen diversity should have increased HLA diversity in relation to their average genomic diversity. We tested this prediction by using HLA class I genetic diversity from 61 human populations. Our results show that human colonization history explains a substantial proportion of HLA genetic diversity worldwide. However, between-population variation at the HLA class I genes is also positively correlated with local pathogen richness (notably for the HLA B gene), thus providing support for the PDBS hypothesis. The proportion of variations explained by pathogen richness is higher for the HLA B gene than for the HLA A and HLA C genes. This is in good agreement with both previous immunological and genetic data suggesting that HLA B could be under a higher selective pressure from pathogens.

Evolution by recombination and transspecies polymorphism in the MHC class I gene of Xenopus laevis

The patterns of major histocompatibility complex (MHC) evolution involve duplications, deletions, and independent divergence of loci during episodes punctuated by natural selection. Major differences in MHC evolution among taxa have previously been attributed to variation in linkage patterns of class I and class II MHC genes. Here we characterize patterns of evolution in the MHC class Ia gene of Xenopus laevis in terms of polymorphism, recombination, and extent of transspecies polymorphism. We also compare these patterns to see if a correlation exists with linkage or separation of the MHC class I and class II regions as seen in amphibians and teleost fishes. In X. laevis, we find high levels of polymorphism. Also, genetic exchange is relatively frequent and occurs in intron II, reshuffling allelic forms of exons 2 and 3. Evolutionary relationships among class I alleles show an intermingling of alleles from divergent Xenopus species rather than a species-specific clustering. Results indicate that the patterns of evolution are similar to those found in salmonid fishes and are different from the mode of evolution seen in primates. Similar patterns of class Ia evolution in salmonid fishes and X. laevis suggest that nonlinkage of class I and class II regions alone is insufficient to explain some patterns of MHC evolution in salmonids.

The evolutionary ecology of the major histocompatibility complex

The major histocompatibility complex (MHC) has become a paradigm for how selection can act to maintain adaptively important genetic diversity in natural populations. Here, we review the contribution of studies on the MHC in non-model species to our understanding of how selection affects MHC diversity, emphasising how ecological and ethological processes influence the tempo and mode of evolution at the MHC, and conversely, how variability at the MHC affects individual fitness, population dynamics and viability. We focus on three main areas: the types of information that have been used to detect the action of selection on MHC genes; the relative contributions of parasite-mediated and sexual selection on the maintenance of MHC diversity; and possible future lines of research that may help resolve some of the unanswered issues associated with MHC evolution.

Serological and molecular diversity in the cattle MHC class I region

Information on major histocompatibility complex (MHC) diversity in cattle is important to aid our understanding of immune responses and may contribute to maintenance of healthy cattle populations. Equally, understanding the mechanisms involved in generating this diversity may shed light on the complex nature of mammalian MHC evolution. The aim of this study was to assess molecular and serological variation within cattle MHC class I molecules and to study the mechanisms generating diversity. To address this aim, sequence variation was examined in 12 serologically assigned alleles from three putative loci and correlated with monoclonal antibody (mAb) binding data. The results demonstrate that both alloantisera and mAbs often fail to distinguish gene products that differ by a significant number of amino acids. Conversely, some mAbs could distinguish alleles differing by only one or two amino acids. Examination of the sequences demonstrates sharing of motifs between alleles, some encoded at distinct loci, supporting the occurrence of interlocus recombination within the cattle MHC class I region. The implications of this for MHC sequence diversity, and functional capability, are discussed.

MHC diversity and the association to nematode parasitism in the yellow-necked mouse (Apodemus flavicollis)

In vertebrates, the genes of the major histocompatibility complex (MHC) are among the most debated candidates accounting for co-evolutionary processes of host-parasite interaction at the molecular level. The exceptionally high allelic polymorphism found in MHC loci is believed to be maintained by pathogen-driven selection, mediated either through heterozygous advantage or rare allele advantage (= frequency dependent selection). While investigations under natural conditions are still very rare, studies on humans or mice under laboratory conditions revealed support for both hypotheses. We investigated nematode burden and allelic diversity of a functional important MHC class II gene (DRB exon2) in free-ranging yellow-necked mice (Apodemus flavicollis). Twenty-seven distinct Apfl-DRB alleles were detected in 146 individuals with high levels of amino acid sequence divergence, especially at the antigen binding sites (ABS), indicating selection processes acting on this locus. Heterozygosity had no influence on the infection status (being infected or not), the number of different nematode infections (NNI) or the intensity of infection, measured as the individual faecal egg count (FEC). However, significant associations of specific Apfl-DRB alleles to both nematode susceptibility and resistance were found, for all nematodes as well as in separate analyses of the two most common nematodes. Apodemus flavicollis individuals carrying the alleles Apfl-DRB*5 or Apfl-DRB*15 revealed significantly higher FEC than individuals with other alleles. In contrast, the allele Apfl-DRB*23 showed a significant association to low FEC of the most common nematode. Thus, our results provide evidence for pathogen-driven selection acting through rare allele advantage under natural conditions.

A genomics approach to the detection of positive selection in cattle: adaptive evolution of the T-cell and natural killer cell-surface protein CD2

The detection of adaptive evolution at the molecular level is of interest not only as an insight into the process of evolution but also because of its functional implications for genes of interest. Here, we present the first genomics approach to detecting positive selection operating on the Bos taurus lineage, an important domestic species. This analysis led to the identification of the T-cell and natural killer (NK) cell receptor cluster of differentiation 2 (CD2) as having a strong signal of selection. Further detailed investigation of CD2 revealed that this gene was subject to positive selection during the evolution of a number of mammalian lineages. Moreover, we show that selection has operated primarily on the extracellular domain of CD2 and discuss the implications of this for an important regulator of the adaptive immune response.

MHC-based patterns of social and extra-pair mate choice in the Seychelles warbler

The existence and nature of indirect genetic benefits to mate choice remain contentious. Major histocompatibility complex (MHC) genes, which play a vital role in determining pathogen resistance in vertebrates, may be the link between mate choice and the genetic inheritance of vigour in offspring. Studies have shown that MHC-dependent mate choice can occur in mammal and fish species, but little work has focused on the role of the MHC in birds. We tested for MHC-dependent mating patterns in the Seychelles warbler (Acrocephalus sechellensis). There was no influence of MHC class I exon 3 variation on the choice of social mate. However, females were more likely to obtain extra-pair paternity (EPP) when their social mate had low MHC diversity, and the MHC diversity of the extra-pair male was significantly higher than that of the cuckolded male. There was no evidence that females were mating disassortatively, or that they preferred males with an intermediate number of MHC bands. Overall, the results are consistent with the 'good genes' rather than the 'genetic compatibility' hypothesis. As female choice will result in offspring of higher MHC diversity, MHC-dependent EPP may provide indirect benefits in the Seychelles warbler if survival is positively linked to MHC diversity.

The complexity of selection at the major primate beta-defensin locus

BACKGROUND

We have examined the evolution of the genes at the major human beta-defensin locus and the orthologous loci in a range of other primates and mouse. For the first time these data allow us to examine selective episodes in the more recent evolutionary history of this locus as well as the ancient past. We have used a combination of maximum likelihood based tests and a maximum parsimony based sliding window approach to give a detailed view of the varying modes of selection operating at this locus.

RESULTS

We provide evidence for strong positive selection soon after the duplication of these genes within an ancestral mammalian genome. Consequently variable selective pressures have acted on beta-defensin genes in different evolutionary lineages, with episodes both of negative, and more rarely positive selection, during the divergence of primates. Positive selection appears to have been more common in the rodent lineage, accompanying the birth of novel, rodent-specific beta-defensin genes. These observations allow a fuller understanding of the evolution of mammalian innate immunity. In both the rodent and primate lineages, sites in the second exon have been subject to positive selection and by implication are important in functional diversity. A small number of sites in the mature human peptides were found to have undergone repeated episodes of selection in different primate lineages. Particular sites were consistently implicated by multiple methods at positions throughout the mature peptides. These sites are clustered at positions predicted to be important for the specificity of the antimicrobial or chemoattractant properties of beta-defensins. Surprisingly, sites within the prepropeptide region were also implicated as being subject to significant positive selection, suggesting previously unappreciated functional significance for this region.

CONCLUSIONS

Identification of these putatively functional sites has important implications for our understanding of beta-defensin function and for novel antibiotic design.

Complex interactions: The immunogenetics of human leukocyte antigen and killer cell immunoglobulin-like receptors

The killer cell immunoglobulin-like receptors (KIR) for human leukocyte antigen (HLA) modulate innate and adaptive immunity by controlling effector cells. HLA and KIR are encoded in genomic regions that have complex organization and exhibit exceptional diversity within and among human population groups. This diversity is likely to have arisen to combat a constantly evolving pathogen challenge. Numerous variations influence the expression level or function of KIR molecules and can affect their interaction with HLA, with important implications for the immune response. The functional variety of natural immune responses that are controlled by HLA and KIR interactions is genetically determined and maintained by natural selection.

Recombination and the origin of sequence diversity in the DRB MHC class II locus in chamois (Rupicapra spp.)

We examined the evolutionary processes contributing to genetic diversity at the major histocompatibility complex (MHC) class II DRB locus in chamois (Rupicapra spp., subfamily Caprinae). We characterised the pattern of intragenic recombination (or homologous gene conversion) and quantified the amount of recombination in the genealogical history of the two chamois species, Pyrenean chamois (Rupicapra pyrenaica) and Alpine chamois (Rupicapra rupicapra). We found evidence for intragenic recombination, and the estimated amount of population recombination suggests that recombination has been a significant process in generating DRB allelic diversity in the genealogical history of the genus Rupicapra. Moreover, positive selection appears to act on the same peptide-binding residues in both analysed chamois species, but not in identical intensity. Recombination coupled with positive selection drives the rapid evolution at the peptide-binding sites in the MHC class II DRB gene. Many chamois MHC class II DRB alleles are thus much younger than previously assumed.

Comparative evolutionary histories of chitinase genes in the Genus zea and Family poaceae

Patterns of DNA sequence diversity vary widely among genes encoding proteins that protect plants against pathogens and herbivores. Comparative studies may help determine whether these differences are due to the strength of selection acting on different types of defense, in different evolutionary lineages, or both. I analyzed sequence diversity at three chitinases, a well-studied component of defense, in two species of Zea and several Poaceae taxa. Although the Zea species are closely related and these genes code for proteins with similar biochemical function, patterns of diversity vary widely within and among species. Intraspecific diversity at chiB, chiI, and Z. mays ssp. parviglumis chiA are consistent with a neutral-equilibrium model whereas chiA had no segregating sites within Z. diploperennis--consistent with a recent and strong selective sweep. Codons identified as having diverged among Poaceae taxa in response to positive selection were significantly overrepresented among targets of selection in Arabis, suggesting common responses to selection in distantly related plant taxa. Divergence of the recent duplicates chiA and chiB is consistent with positive selection but relaxed constraint cannot be rejected. Weak evidence for adaptive divergence of these duplicated downstream components of defense contrasts with strong evidence for adaptive divergence of genes involved in pathogen recognition.

A temporal analysis shows major histocompatibility complex loci in the Scandinavian wolf population are consistent with neutral evolution

The major histocompatibility complex (MHC) has an integral role in the immune system, and hence diversity at its genes may be of particular importance for the health of populations. In large populations, balancing selection maintains diversity in MHC genes, but theoretical expectations indicate that this form of selection is absent or inefficient in small populations. We examine the level of diversity at three MHC class II loci in the wolf population of Scandinavia, a population naturally recolonized with a genetic contribution from as few as three founders, and in four neighbouring wolf populations. In the Scandinavian wolf population, two alleles were found for each locus and the distribution of alleles is compatible with their linkage into two haplotypes. Changes in the level of heterozygosity over time since recolonization demonstrate the effects of the proposed arrival of an immigrant wolf. The maintenance of diversity is shown to be compatible with a neutral, random allocation of alleles, in conjunction with crossing between packs. A total of 15 DRB1, seven DQA and 10 DQB1 alleles are found in four neighbouring wolf populations, with substantial sharing across populations. Even in these larger populations, bottlenecks and fragmentation with consequent genetic drift are likely to have resulted in few indicators for balancing selection and significant differentiation of populations.

Variability and evolution of bovine β-defensin genes

Defensins comprise an important family of antimicrobial peptides. Among vertebrates numerous defensin genes have been detected, but their evolutionary background is still discussed. We investigated the molecular evolution of bovine defensins via screening of different bovine species including the extinct ancestor of domestic cattle (Bos primigenius) for beta-defensin encoding genes. We detected a large variability of new defensin encoding sequences similar to previously published bovine neutrophil beta-defensin (bnbd), neutrophil beta-defensin 12 (nbd12), enteric beta-defensin (ebd), lingual antimicrobial peptide (lap), and tracheal antimicrobial peptide (tap). Our data suggest that variants of the same so-called subfamily (tap, lap, ebd, and nbd12) each share a common ancestry independent of their species origin, implicating several duplication events of tap, lap, ebd, and nbd12 before the different bovine lineages diverged. Variants of bovine neutrophil beta-defensins bnbd5 and bnbd9 were detected exclusively in domestic cattle and aurochs. Values of synonymous and nonsynonymous substitutions demonstrated lap, bnbd5, bnbd9 and nbd12 evolving under positive selection, whereas amino-acid altering substitutions among variants of ebd and tap are purified. Comparison of the amino-acid sequences with available structures of human and murine defensins suggested conservation of the typical secondary elements of defensins in the absence of high sequence similarity.

Association between major histocompatibility complex class II DRB alleles and parasite load in the hairy-footed gerbil, Gerbillurus paeba, in the southern Kalahari

We investigated the importance of the MHC-constitution (major histocompatibility complex-constitution) on the endoparasite load in free-range hairy-footed gerbils (Gerbillurus paeba) in the southern Kalahari Desert. While the number of alleles of the duplicated DRB exon 2 gene had no significant effects on the individual status of being 'not infected' or 'infected' and on the number of helminth morphotype infections per individual, it significantly affected the faecal egg count values. One allele (Gepa-DRB*15) was only found in uninfected mice. Our results support the hypotheses that MHC polymorphism in G. paeba is maintained by pathogen-driven selection. The present study is the first investigation on associations between duplicated DRB gene loci and the parasite load in mammals.

The Evolution of Alloimmunity and the Genesis of Adaptive Immunity

Infectious agents select for host immune responses that destroy infectious nonself yet maintain tolerance to self. Here we propose that retroviruses and other host-antigen associated pathogens (HAAPs) select for the genetic, biochemical, and cell biological properties of alloimmunity, also known as the histocompatibility or tissue rejection response. This hypothesis predicts the major observations regarding histocompatibility responses, including: (i) their existence in animals as diverse as sponges and humans; (ii) extreme polymorphism and balanced allele frequencies at histocompatibility loci, including the human MHC and blood group loci; (iii) the frequency dependent selection of histocompatibility alleles; (iv) the ancient age of many alloantigenic polymorphisms; (v) the high ratio of nonsynonymous mutations to synonymous mutations at histocompatibility loci; (vi) disassortative mating based on MHC alleles; (vii) the inability to explain the existence and continuing selection of histocompatibility alleles by other more conventional biochemical and genetic paradigms; and (viii) the susceptibility of HAAPs, particularly retroviruses such as HIV (human immunodeficiency virus), to histocompatibility reactions. In addition, the hypothesis that HAAPs select the forms and molecules of alloimmunity offers simple explanations for the evolution of histocompatibility systems over time, the initial selection of hypervariable immune mechanisms, and the genesis of adaptive immunity.

Identification and sequence analysis of chicken Toll-like receptors

Toll-like receptors (TLRs) play an important role in the recognition of microbial components. Only chicken TLR2 and -4 have been reported in the literature. The objectives of this study were to identify new chicken TLRs and to evaluate evolutionary significance of these receptors. Searching chicken genomic databases and DNA sequencing revealed five new TLRs, TLR1 (type 1 and 2), -3, -5, and -7. No chicken orthologues of mammalian TLR8, -9, or -10 were found. As in mammals, all chicken TLRs (chTLRs) share identical protein secondary structure that consists of several leucine-rich domains, a transmembrane domain, and Toll/Interleukin-1 receptor domain(s). Phylogenetic analyses indicate that the identified chTLR genes are the orthologues of TLRs in mammals. Analyses of the number of synonymous substitutions per synonymous site and nonsynonymous substitutions per nonsynonymous site indicate that the nucleotide sequences coding for the leucine-rich repeats of chicken TLR1 type 1 and type 2 were significantly under positive Darwinian selection. In contrast, the sequences of other TLRs were under purifying selection. These results support the hypothesis that one of the major evolutionary strategies of the innate immune system is to recognize a few highly conserved microbial components with several conserved TLRs. The results also indicate that the sequence changes in the ligand-binding domains of TLR1 in chickens provide adaptive advantages during evolution.

Evolution of MHC-DRB class II polymorphism in the genus Apodemus and a comparison of DRB sequences within the family Muridae (Mammalia: Rodentia)

Allelic diversity at major histocompatibility complex (MHC) genes is thought to be maintained by balancing selection over long periods of time, even across multiple speciation events. Trans-species sharing of MHC alleles among genera has been supported by many studies on mammals and fish, but in rodents, the results are ambiguous. We investigated natural levels of MHC- DRB variability and evolutionary processes in the wood mouse (Apodemus sylvaticus) and the yellow-necked mouse (Apodemus flavicollis), which are common, sympatric murid rodents in European forests. Using single-strand conformation polymorphism analysis and DNA sequencing, 38 DRB exon 2 alleles were detected among 162 A. sylvaticus from nine different locations in Germany and Switzerland, and 15 DRB exon 2 alleles were detected among 60 A. flavicollis from three different locations in northern Germany. There was evidence for balancing selection in both species. Phylogenetic analysis, including additional murid taxa, showed that the DRB exon 2 sequences did not separate according to species, consistent with trans-species evolution of the MHC in these taxa.

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

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

Mouse biodiversity in the genomic era

Comparative genomics has developed by comparison of distantly related genomes, for which the link between the reported evolutionary changes and species development/physiology/ecology is not obvious. It is argued that the mouse (genus Mus) is an optimal model for microevolutionary genomics in vertebrates. This is because the mouse genome sequence, physical and genetic map have been completed, because mouse genetics, morpho-anatomy, pathology, behavior and ecology are well-studied, and because the Mus genus is a diverse, well- documented taxon, allowing comparative studies at the level of individual, population, subspecies, and species. The potential of the interaction between mouse genome and mouse biodiversity is illustrated by recent studies of speciation in the house mouse Mus musculus, and studies about the evolution of isochores, the peculiar pattern of GC-content variation across mammalian genomes.

Low diversity in the major histocompatibility complex class II DRB1 gene of the Spanish ibex, Capra pyrenaica

During the last two centuries, the Spanish ibex (Capra pyrenaica) has shown a significant demographic decline as a result of the progressive destruction of its natural habitat, disease epidemics, and uncontrolled hunting. Partial sequencing of the class II MHC DRB1 gene revealed that the Spanish ibex has remarkably low levels of genetic variation at this locus, with only six different DRB1 alleles and an observed heterozygosity of 0.429-0.579. The rates of nonsynonymous vs synonymous substitutions were significantly different in the peptide-binding region (dN/dS=5.347, P=0.002), a feature that indicates that the DRB1 gene is under positive selection. A phylogenetic analysis of the Spanish ibex and a set of domestic goat DRB1 alleles revealed that the reported sequences represent four major allelic lineages. The limited allelic repertoire of the DRB1 gene in the Spanish ibex is likely the direct result of the recent history of population bottlenecks and marked demographic decline of this species. A genetic survey of 13 microsatellite loci was consistent with this idea. The Spanish ibex subspecies C. p. hispanica and C. p. victoriae consistently showed considerably lower levels of microsatellite heterozygosity (Ho=0.184-0.231) and allelic diversity (mean number of alleles per locus=2-2.4) than those reported in other wild ruminants. This study demonstrates the significance of both natural selection and the demographic history of populations in determining patterns of genetic variation at MHC loci. In addition, our results emphasize the importance of locally adapted populations for the preservation of genetic diversity.

Diversity of the ovine DQA2 gene1

Variation in the ovine DQA2 gene was investigated in approximately 2,000 sheep from six breeds. Fragments of DNA containing the ovine DQA2 exon 2 were amplified using PCR. Single-strand conformational polymorphism analysis and DNA sequence analysis were employed to detect genetic variation. Twenty-three nucleic acid sequences, encoding 22 DQA2 amino acid sequences, were identified. This increases the number of alleles identified from 10 to 23. In some cases, three or four unique sequences were isolated from individual sheep, suggesting that these DQA2 sequences may represent two loci. Phylogenetic tree analysis revealed that 5 of these 23 sequences were more closely related to cattle DQA3 or DQA4 sequences than to other sheep DQA2 sequences. These sequences clustered together and were called DQA2-like to differentiate them from other DQA2 sequences. There was no evidence of DQA5-like sequences in sheep. Information theory-based analysis indicated that some of the DQA2-like sequences had low information content at splice sites, suggesting that these alleles may have low functional activity. Allelic lineages were observed not only at the DQA2 locus, but also at the DQA2-like locus, supporting the trans-species mode of evolution of MHC genes. Comparison of the allelic sequences suggests that polymorphism seems to have arisen largely by point mutation and gene conversion, and a recent gene conversion event seems to have occurred between the DQA2 and DQA2-like loci. The high level of sequence polymorphism detected and varied number of loci demonstrate the extensive diversity of the ovine DQA2 gene.

Human genes that limit AIDS

Discernable genetic variation among people and populations has an important role in infectious disease epidemics, including that of acquired immune deficiency syndrome (AIDS). Genetic association analysis of several large AIDS cohorts implicate 14 AIDS restriction genes, polymorphic variants in loci that regulate HIV-1 cell entry, acquired and innate immunity, and cytokine defenses to HIV-1. The influence and translational impact of these genes on individual and population sensitivity to AIDS is considerable.

The Evolution of SMC Proteins: Phylogenetic Analysis and Structural Implications

The SMC proteins are found in nearly all living organisms examined, where they play crucial roles in mitotic chromosome dynamics, regulation of gene expression, and DNA repair. We have explored the phylogenetic relationships of SMC proteins from prokaryotes and eukaryotes, as well as their relationship to similar ABC ATPases, using maximum-likelihood analyses. We have also investigated the coevolution of different domains of eukaryotic SMC proteins and attempted to account for the evolutionary patterns we have observed in terms of available structural data. Based on our analyses, we propose that each of the six eukaryotic SMC subfamilies originated through a series of ancient gene duplication events, with the condensins evolving more rapidly than the cohesins. In addition, we show that the SMC5 and SMC6 subfamily members have evolved comparatively rapidly and suggest that these proteins may perform redundant functions in higher eukaryotes. Finally, we propose a possible structure for the SMC5/SMC6 heterodimer based on patterns of coevolution.

Hypervariability, suppressed recombination and the genetics of individuality

We define 'genetic individuality' as intraspecies variation that has substantial heritability and involves traits that are sufficiently common that they can be observed in any modest-sized sampling of individuals. We propose that genetic individuality is largely shaped by the combinatory shuffling of a modest number of genes, each of which exists as a family of functionally and structurally diverged alleles. Unequivocal examples of such allele families are found at the O-antigen-biosynthetic locus in Pseudomonas aeruginosa and the human leucocyte antigen locus in humans. We examine characteristic features of these allele families and explore the possibility that genetic loci with similar characteristics can be recognized in a whole-genome scan of human genetic variation.

Heterozygote advantage fails to explain the high degree of polymorphism of the MHC

Major histocompatibility (MHC) molecules are encoded by extremely polymorphic genes and play a crucial role in vertebrate immunity. Natural selection favors MHC heterozygous hosts because individuals heterozygous at the MHC can present a larger diversity of peptides from infectious pathogens than homozygous individuals. Whether or not heterozygote advantage is sufficient to account for a high degree of polymorphism is controversial, however. Using mathematical models we studied the degree of MHC polymorphism arising when heterozygote advantage is the only selection pressure. We argue that existing models are misleading in that the fitness of heterozygotes is not related to the MHC alleles they harbor. To correct for this, we have developed novel models in which the genotypic fitness of a host directly reflects the fitness contributions of its MHC alleles. The mathematical analysis suggests that a high degree of polymorphism can only be accounted for if the different MHC alleles confer unrealistically similar fitnesses. This conclusion was confirmed by stochastic simulations, including mutation, genetic drift, and a finite population size. Heterozygote advantage on its own is insufficient to explain the high population diversity of the MHC.

Co-duplication of olfactory receptor and MHC class I genes in the mouse major histocompatibility complex

We report the 897 kb sequence of a cluster of olfactory receptor (OR) genes located at the distal end of the major histocompatibility complex (MHC) class I region on mouse chromosome 17 of strain 129/SvJ (H2bc). With additional information from the mouse genome draft sequence, we identified 59 OR loci (approximately 20% pseudogenes) in contrast to only 25 OR loci (approximately 50% pseudogenes) in the corresponding centromeric OR cluster that is part of the 'extended MHC class I region' on human chromosome 6. Comparative analysis leads to three major observations: (i) most of the OR subfamilies have evolved independently in the two species, expanding more in the mouse, and resulting in co-orthologs--subfamilies of highly similar paralogs that keep orthologous relationships with their human counterparts; (ii) three of the mouse OR subfamilies have no orthologs in humans; and (iii) MHC class I loci are interspersed in the OR cluster in mouse but not in human, and were subjected to co-duplication with OR genes. Screening of our sequence against the available sequences of other strains/haplotypes revealed that most of the OR loci are polymorphic and that the number of OR loci may vary among strains/haplotypes. Our findings that MHC-linked OR loci share duplication with MHC class I loci, have duplicated extensively and are polymorphic revives questions about potential reciprocal influences acting on the dynamics and evolution of the H2 region and the H2-linked OR loci.

Parasitic exploitation as an engine of diversity

Parasitic exploitation occurs within and between a wide variety of taxa in a plethora of diverse contexts. Theoretical and empirical analyses indicate that parasitic exploitation can generate substantial genetic and phenotypic polymorphism within species. Under some circumstances, parasitic exploitation may also be an important factor causing reproductive isolation and promoting speciation. Here we review research relevant to the relationship between parasitic exploitation, within species-polymorphism, and speciation in some of the major arenas in which such exploitation has been studied. This includes research on the vertebrate major histocompatibility loci, plant-pathogen interactions, the evolution of sexual reproduction, intragenomic conflict, sexual conflict, kin mimicry and social parasitism, tropical forest diversity and the evolution of language. We conclude by discussing some of the issues raised by comparing the effect of parasitic exploitation on polymorphism and speciation in different contexts.

Detecting excess radical replacements in phylogenetic trees

There are a few instances in which positive Darwinian selection has been convincingly demonstrated at the molecular level. In this study, we present a novel test for detecting excess of radical amino-acid replacements. Such excess is usually indicative of positive Darwinian selection, but may also be due to relaxed functional constraints or model misspecification. In our test, each amino-acid replacement is characterized in terms of a physicochemical distance, i.e., the degree of dissimilarity between the exchanged amino-acid residues. By using phylogenetic trees based on protein sequences, our test identifies statistically significant deviations of the mean physicochemical distance from the random expectation, either along a taxonomic lineage or across a subtree. The mean inferred distance is calculated as the average physicochemical distance over all possible ancestral sequence reconstructions weighted by their likelihood. Our method substantially improves over previous approaches by taking into account the stochastic process, tree phylogeny, among-site rate variation, and alternative ancestral reconstructions. We provide a fast linear time algorithm for applying this test to all branches and all subtrees of a given phylogenetic tree. We validate this approach by applying it to two well-studied datasets: the MHC class I glycoproteins serving as a positive control, and the house-keeping gene carbonic anhydrase I serving as a negative control.

MHC diversity in two Acrocephalus species: the outbred Great reed warbler and the inbred Seychelles warbler

The Great reed warbler (GRW) and the Seychelles warbler (SW) are congeners with markedly different demographic histories. The GRW is a normal outbred bird species while the SW population remains isolated and inbred after undergoing a severe population bottleneck. We examined variation at Major Histocompatibility Complex (MHC) class I exon 3 using restriction fragment length polymorphism, denaturing gradient gel electrophoresis and DNA sequencing. Although genetic variation was higher in the GRW, considerable variation has been maintained in the SW. The ten exon 3 sequences found in the SW were as diverged from each other as were a random sub-sample of the 67 sequences from the GRW. There was evidence for balancing selection in both species, and the phylogenetic analysis showing that the exon 3 sequences did not separate according to species, was consistent with transspecies evolution of the MHC.

Geographic patterns of functional categories of HLA-DRB1 alleles: a new approach to analyse associations between HLA-DRB1 and disease

Because specific amino acids found within the peptide-binding cleft of human leukocyte antigen (HLA) molecules have been implicated in HLA/disease associations, an approach which consists in grouping the alleles according to their functional properties at the protein level may enable us to better understand HLA associations than the conventional allelic classification. In this study, we applied this methodology to investigate the associations between HLA-DRB1 and rheumatoid arthritis. The alleles were first classified into seven functional categories [restrictive supertype patterns (RSPs)], among which three were known to be significantly associated with susceptibility (one category) or resistance (two categories) to rheumatoid arthritis. The frequencies of these categories were then estimated in 104 population samples previously tested for HLA-DRB1, and their variability was analysed spatially on a worldwide scale by applying an original methodology for detecting discontinuities in geographically patterned data. RSP frequencies were also compared to known values of rheumatoid arthritis prevalence in some populations. The results indicated that the three RSP frequency distributions were geographically structured, and that these patterns could generally be explained by the history of human migrations. However, the peculiar pattern observed for RSP 'A' (conferring susceptibility to rheumatoid arthritis) indicated a possible association with some latitude-dependent disease. Furthermore, the very high correlation coefficient found between RSP 'A' frequencies and rheumatoid arthritis prevalence confirmed the significant disease association of this functional category. In contrast, the putative protective effect of the other RSPs ('De' and 'Q') was not detectable at the worldwide level, but may be significant in specific geographic areas. This study shows that population genetic diversity analyses based on a functional grouping of HLA alleles provide an efficient way to explore the mutual influence of HLA genetic variation and disease.

Effects of habitat fragmentation and changes of dispersal behaviour after a recent population decline on the genetic variability of noncoding and coding DNA of a monogamous Malagasy rodent

While interactions among demography, behaviour and genetic structure are well-documented for neutral genetic markers, the role of these parameters and the effects of genetic drift and selection are considerably less well understood in functional genes, such as the major histocompatibility complex (MHC). In this study, the consequences of habitat fragmentation and the effects of a current population decline on noncoding (mitochondrial DNA) and two coding MHC loci (DQA, DRB) with different functional importance were investigated in the small remnant subdivided population of the endangered Malagasy giant jumping rat (Hypogeomys antimena). Both neutral and selective markers revealed a significant genetic differentiation between the two remnant subpopulations. The FST values were much lower in the MHC DQA and DRB genes than in the mitochondrial data. The MHC DRB loci display the effects of both balancing selection (high sequence diversity, four times higher nonsynonymous than synonymous substitutions in the functionally important antigen-binding site positions, twice the average heterozygosity of individual amino acids at the positions identified as part of the antigen-binding site (ABS) than those outside the ABS and nonselective forces including genetic drift. Simultaneously with a current population decline offspring reduced their dispersal distances. No substantial effects were detected within the first 6 years of reduced gene flow in either mitochondrial or MHC markers.

The influence of HLA genotype on AIDS

Genetic resistance to infectious diseases is likely to involve a complex array of immune-response and other genes with variants that impose subtle but significant consequences on gene expression or protein function. We have gained considerable insight into the genetic determinants of HIV-1 disease, and the HLA class I genes appear to be highly influential in this regard. Numerous reports have identified a role for HLA genotype in AIDS outcomes, implicating many HLA alleles in various aspects of HIV disease. Here we review the HLA associations with progression to AIDS that have been consistently affirmed and discuss the underlying mechanisms behind some of these associations based on functional studies of immune cell recognition.

Comparative genome organization of human, murine, and feline MHC class II region

To study comparative molecular dynamics in the genesis of the major histocompatibility complex (MHC), we determined a complete nucleotide sequence spanning 758,291 bp of the domestic cat (Felis catus) extended and classical class II region. The feline class II MHC includes 44 genes (31 predicted to be expressed) which display DNA sequence homology and ordered gene synteny with human HLA and mouse H2, in extended class II and centromere proximal regions (DM to DO) of the classical class II region. However, remarkable genomic alterations including gene gain and loss plus size differentials of 250 kb are evident in comparisons of the cat class II with those of human and mouse. The cat MHC lacks the entire DQ region and retains only relict pseudogene homologs of DP genes, compensated by expansion and reorganization of seven modern DR genes. Repetitive gene families within the feline MHC comprise 35% of the feline MHC with very different density and abundance of GC levels, SINES, LINES, STRs, and retro-elements from the same repeats in human and mouse MHC. Comparison of the feline MHC with the murine and human MHC offers a detailed view of the consequences of genome organization in three mammalian lineages.

Effects of MHC class I on HIV/SIV disease in primates

Data indicate that resistance to HIV-1 disease involves an array of contrasting HLA genotypic effects that are subtle, but significant, particularly when these genetic effects are considered as a whole. Numerous reports attributing a role for HLA genotype in AIDS outcomes have been reported, and a few of these have been affirmed in multiple studies. Functional studies of immune cell recognition have provided clues to the underlying mechanisms behind some of the strongest HLA associations, suggesting the means by which relative resistance or susceptibility to the virus may occur. SIV infection in non-human primates has served as an invaluable model for understanding AIDS pathogenesis (in rhesus monkeys) and viral resistance (in chimpanzee). The effect of rhesus MHC class I molecules on the evolution of SIV has been convincingly described [19], and a recent study in humans has suggested that selection pressure conferred by HLA molecules is responsible for specific genetic variation in HIV-1 [114]. HIV-1 may eventually have conspicuous evolutionary effects on HLA and other AIDS restriction genes, a prolonged process that could have occurred in chimpanzee [92]. To prevent such an outcome, it will be necessary to approach the disease from many perspectives, andapply comprehensively the knowledge gained to the successful control of the virus.

Evolution of MHC class II E beta diversity within the genus Peromyscus

Progress in understanding the evolution of variation at the MHC has been slowed by an inability to assess the relative roles of mutation vs. intragenic recombination in contributing to observed polymorphism. Recent theoretical advances now permit a quantitative treatment of the problem, with the result that the amount of recombination is at least an order of magnitude greater than that of mutation in the history of class II genes. We suggest that this insight allows progress in evaluating the importance of other factors affecting the evolution of the MHC. We investigated the evolution of MHC class II E beta sequence diversity in the genus Peromyscus. We find evidence for extensive recombination in the history of these sequences. Nevertheless, it appears that intragenic recombination alone is insufficient to account for evolution of MHC diversity in Peromyscus. Significant differences in silent variation among subgenera arose over a relatively short period of time, with little subsequent change. We argue that these observations are consistent with the effects of historical population bottleneck(s). Population restrictions may explain general features of MHC evolution, including the large amount of recombination in the history of MHC genes, because intragenic recombination may efficiently regenerate allelic polymorphism following a population constriction.

Duplication and selection in the evolution of primate beta-defensin genes

BACKGROUND

Innate immunity is the first line of defense against microorganisms in vertebrates and acts by providing an initial barrier to microorganisms and triggering adaptive immune responses. Peptides such as beta-defensins are an important component of this defense, providing a broad spectrum of antimicrobial activity against bacteria, fungi, mycobacteria and several enveloped viruses. Beta-defensins are small cationic peptides that vary in their expression patterns and spectrum of pathogen specificity. Disruptions in beta-defensin function have been implicated in human diseases, including cystic fibrosis, and a fuller understanding of the variety, function and evolution of human beta-defensins might form the basis for novel therapies. Here we use a combination of laboratory and computational techniques to characterize the main human beta-defensin locus on chromosome 8p22-p23.

RESULTS

In addition to known genes in the region we report the genomic structures and expression patterns of four novel human beta-defensin genes and a related pseudogene. These genes show an unusual pattern of evolution, with rapid divergence between second exon sequences that encode the mature beta-defensin peptides matched by relative stasis in first exons that encode signal peptides.

CONCLUSIONS

We conclude that the 8p22-p23 locus has evolved by successive rounds of duplication followed by substantial divergence involving positive selection, to produce a diverse cluster of paralogous genes established before the human-baboon divergence more than 23 million years ago. Positive selection, disproportionately favoring alterations in the charge of amino-acid residues, is implicated as driving second exon divergence in these genes.

Genetic variation and population structure of moose (Alces alces) at neutral and functional DNA loci

Genetic variation was examined for moose (Alces alces) from Riding Mountain, Isle Royale, and Pukaskwa national parks; northwestern, Nipigon, northeastern, and central Ontario; New Brunswick; and Newfoundland. The national parks were identified as maintaining potentially different local selection pressures due to the absence of hunting and the presence or absence of the parasite Parelaphostrongylus tenuis. Genetic variation was estimated using neutral DNA markers, assessed by multilocus DNA fingerprinting and five microsatellite loci, and the functional antigen binding region (ARS) (exon 2) of the major histocompatibility complex (MHC) gene DRB. There was discordance in the allelic diversity observed at the neutral loci compared with the MHC DRB locus in a number of populations. Ontario populations demonstrated higher levels of variability at the neutral loci and relatively low levels at the DRB locus. Conversely, the Isle Royale population has the lowest genetic variability, consistent with a historic small founding event, at the neutral DNA markers and relatively high variability at the MHC gene. Relatively high levels of genetic variation at the DRB locus were observed in protected park populations concomitant with the absence of white-tailed deer (Odocoileus virginianus) or the parasite P. tenuis and an absence of hunting. Gene flow was observed among the neighboring geographic regions within Ontario, including Pukaskwa National Park, with evidence of isolation-by-distance among more distant regions within Ontario. The discordant patterns between DNA markers suggest that neutral DNA markers may not accurately reflect adaptive variation present at functional loci.

Efficient discovery of single-nucleotide polymorphisms in coding regions of human genes

Single nucleotide polymorphisms in protein coding regions (cSNPs) are of great interest for their effects on phenotype and potential for mapping disease genes. We have identified 5,400 novel exonic SNPs from alignments of public EST data to the draft human genome sequence, and approximately 12,000 more novel exonic SNPs from EST cluster alignments. We found 82% of the genomic-aligned SNPs and 63% of the EST-only SNPs to be detectably polymorphic in 20 Finnish DNA samples. 37% of the SNPs mapped to known protein coding regions, yielding 6,500 distinct, novel cSNPs from the two datasets. These data reveal selection against mutations that alter protein structure, and distinct classes of genes under strongly positive vs. negative pressure from natural selection for amino acid replacement (detected by K(A)/K(S)ratio). We have searched these cSNPs for compatibility with the amino acid profile at each site and structural impact on protein core stability.

Multiple parasites are driving major histocompatibility complex polymorphism in the wild

Parasite mediated selection may result in arms races between host defence and parasite virulence. In particular, simultaneous infections from multiple parasite species should cause diversification (i.e. balancing selection) in resistance genes both at the population and the individual level. Here, we tested these ideas in highly polymorphic major histocompatibility complex (MHC) genes from three-spined sticklebacks (Gasterosteus aculeatus L.). In eight natural populations, parasite diversity (15 different species), and MHC class IIB diversity varied strongly between habitat types (lakes vs. rivers vs. estuaries) with lowest values in rivers. Partial correlation analysis revealed an influence of parasite diversity on MHC class IIB variation whereas general genetic diversity assessed at seven microsatellite loci was not significantly correlated with parasite diversity. Within individual fish, intermediate, rather than maximal allele numbers were associated with minimal parasite load, supporting theoretical models of self-reactive T-cell elimination. The optimal individual diversity matched those values female fish try to achieve in their offspring by mate choice. We thus present correlative evidence supporting the 'allele counting' strategy for optimizing the immunocompetence in stickleback offspring.