Maintenance of multiallelic polymorphism at the MHC region

Diversity of major histocompatibility complex of II B gene and mate choice in a monogamous and long-lived seabird, the Little Auk (Alle alle)

The major histocompatibility complex (MHC) plays a key role in the adaptive immune system of vertebrates, and is known to influence mate choice in many species. In birds, the MHC has been extensively examined but mainly in galliforms and passerines while other taxa that represent specific ecological and evolutionary life-histories, like seabirds, are underexamined. Here, we characterized diversity of MHC Class II B exon 2 in a colonial pelagic seabird, the Little Auk (or Dovekie Alle alle). We further examined whether MHC variation could be maintained through balancing selection and disassortative mating. We found high polymorphism at the genotyped MHC fragment, characterizing 99 distinct alleles across 140 individuals from three populations. The alleles frequencies exhibited a similar skewed distribution in both sexes, with the four most commonly occurring alleles representing approximately 35% of allelic variation. The results of a Bayesian site-by-site selection analysis suggest evidence of balancing selection and no direct evidence for MHC-dependent disassortative mating preferences in the Little Auk. The latter result might be attributed to the high overall polymorphism of the examined fragment, which itself may be maintained by the large population size of the species.

High Resolution HLA ∼A, ∼B, ∼C, ∼DRB1, ∼DQA1, and ∼DQB1 Diversity in South African Populations

Background: Lack of HLA data in southern African populations hampers disease association studies and our understanding of genetic diversity in these populations. We aimed to determine HLA diversity in South African populations using high resolution HLA ∼A, ∼B, ∼C, ∼DRB1, ∼DQA1 and ∼DQB1 data, from 3005 previously typed individuals.

Methods: We determined allele and haplotype frequencies, deviations from Hardy-Weinberg equilibrium (HWE), linkage disequilibrium (LD) and neutrality test. South African HLA class I data was additionally compared to other global populations using non-metrical multidimensional scaling (NMDS), genetic distances and principal component analysis (PCA).

Results: All loci strongly (p < 0.0001) deviated from HWE, coupled with excessive heterozygosity in most loci. Two of the three most frequent alleles, HLA ∼DQA1*05:02 (0.2584) and HLA ∼C*17:01 (0.1488) were previously reported in South African populations at lower frequencies. NMDS showed genetic distinctness of South African populations. Phylogenetic analysis and PCA clustered our current dataset with previous South African studies. Additionally, South Africans seem to be related to other sub-Saharan populations using HLA class I allele frequencies.

Discussion and Conclusion: Despite the retrospective nature of the study, data missingness, the imbalance of sample sizes for each locus and haplotype pairs, and induced methodological difficulties, this study provides a unique and large HLA dataset of South Africans, which might be a useful resource to support anthropological studies, disease association studies, population based vaccine development and donor recruitment programs. We additionally provide simulated high resolution HLA class I data to augment the mixed resolution typing results generated from this study.

Flexible Mixture Model Approaches That Accommodate Footprint Size Variability for Robust Detection of Balancing Selection

Long-term balancing selection typically leaves narrow footprints of increased genetic diversity, and therefore most detection approaches only achieve optimal performances when sufficiently small genomic regions (i.e., windows) are examined. Such methods are sensitive to window sizes and suffer substantial losses in power when windows are large. Here, we employ mixture models to construct a set of five composite likelihood ratio test statistics, which we collectively term B statistics. These statistics are agnostic to window sizes and can operate on diverse forms of input data. Through simulations, we show that they exhibit comparable power to the best-performing current methods, and retain substantially high power regardless of window sizes. They also display considerable robustness to high mutation rates and uneven recombination landscapes, as well as an array of other common confounding scenarios. Moreover, we applied a specific version of the B statistics, termed B2, to a human population-genomic data set and recovered many top candidates from prior studies, including the then-uncharacterized STPG2 and CCDC169-SOHLH2, both of which are related to gamete functions. We further applied B2 on a bonobo population-genomic data set. In addition to the MHC-DQ genes, we uncovered several novel candidate genes, such as KLRD1, involved in viral defense, and SCN9A, associated with pain perception. Finally, we show that our methods can be extended to account for multiallelic balancing selection and integrated the set of statistics into open-source software named BalLeRMix for future applications by the scientific community.

Footprints of natural selection at the mannose-6-phosphate isomerase locus in barnacles

The rocky intertidal is a natural laboratory to study how natural selection acts on the genes and proteins responsible for organismal survival and reproduction. Alternative forms of enzymes that differ across the intertidal have been known for decades and have provided examples of selection, but the genetic basis of such enzyme variation is known in only a few cases. In this paper, we present molecular evidence of natural selection at the Mpi gene, a key enzyme in energy metabolism that alters survival of barnacles living across the stress gradient imposed by the intertidal. Our study demonstrates how natural selection can facilitate survival in highly heterogeneous environments through the maintenance of multiple molecular solutions to ecological stresses.

The mannose-6-phosphate isomerase (Mpi) locus in Semibalanus balanoides has been studied as a candidate gene for balancing selection for more than two decades. Previous work has shown that Mpi allozyme genotypes (fast and slow) have different frequencies across Atlantic intertidal zones due to selection on postsettlement survival (i.e., allele zonation). We present the complete gene sequence of the Mpi locus and quantify nucleotide polymorphism in S. balanoides, as well as divergence to its sister taxon Semibalanus cariosus. We show that the slow allozyme contains a derived charge-altering amino acid polymorphism, and both allozyme classes correspond to two haplogroups with multiple internal haplotypes. The locus shows several footprints of balancing selection around the fast/slow site: an enrichment of positive Tajima’s D for nonsynonymous mutations, an excess of polymorphism, and a spike in the levels of silent polymorphism relative to silent divergence, as well as a site frequency spectrum enriched for midfrequency mutations. We observe other departures from neutrality across the locus in both coding and noncoding regions. These include a nonsynonymous trans-species polymorphism and a recent mutation under selection within the fast haplogroup. The latter suggests ongoing allelic replacement of functionally relevant amino acid variants. Moreover, predicted models of Mpi protein structure provide insight into the functional significance of the putatively selected amino acid polymorphisms. While footprints of selection are widespread across the range of S. balanoides, our data show that intertidal zonation patterns are variable across both spatial and temporal scales. These data provide further evidence for heterogeneous selection on Mpi.

LOWER MITE INFESTATIONS IN AN ASEXUAL GECKO COMPARED WITH ITS SEXUAL ANCESTORS

What advantage do sexually reproducing organisms gain from their mode of reproduction that compensates for their twofold loss in reproductive rate relative to their asexual counterparts? One version of the Red Queen hypothesis suggests that selective pressure from parasites is strongest on the most common genotype in a population, and thus genetically identical clonal lineages are more vulnerable to parasitism over time than genetically diverse sexual lineages. Our surveys of the ectoparasites of an asexual gecko and its two sexual ancestral species show that the sexuals have a higher prevalence, abundance, and mean intensity of mites than asexuals sharing the same habitat. Our experimental data indicate that in one sexual/asexual pair this pattern is at least partly attributable to higher attachment rates of mites to sexuals. Such a difference may occur as a result of exceptionally high susceptibility of the sexuals to mites because of their low genetic diversity (relative to other more-outbred sexual species) and their potentially high stress levels, or as a result of exceptionally low susceptibility of the asexuals to mites because of their high levels of heterozygosity.

MHC allele frequency distributions under parasite-driven selection: A simulation model

BACKGROUND

The extreme polymorphism that is observed in major histocompatibility complex (MHC) genes, which code for proteins involved in recognition of non-self oligopeptides, is thought to result from a pressure exerted by parasites because parasite antigens are more likely to be recognized by MHC heterozygotes (heterozygote advantage) and/or by rare MHC alleles (negative frequency-dependent selection). The Ewens-Watterson test (EW) is often used to detect selection acting on MHC genes over the recent history of a population. EW is based on the expectation that allele frequencies under balancing selection should be more even than under neutrality. We used computer simulations to investigate whether this expectation holds for selection exerted by parasites on host MHC genes under conditions of heterozygote advantage and negative frequency-dependent selection acting either simultaneously or separately.

RESULTS

In agreement with simple models of symmetrical overdominance, we found that heterozygote advantage acting alone in populations does, indeed, result in more even allele frequency distributions than expected under neutrality, and this is easily detectable by EW. However, under negative frequency-dependent selection, or under the joint action of negative frequency-dependent selection and heterozygote advantage, distributions of allele frequencies were less predictable: the majority of distributions were indistinguishable from neutral expectations, while the remaining runs resulted in either more even or more skewed distributions than under neutrality.

CONCLUSIONS

Our results indicate that, as long as negative frequency-dependent selection is an important force maintaining MHC variation, the EW test has limited utility in detecting selection acting on these genes.

Reintroducing group selection to the human behavioral sciences

In both biology and the human sciences, social groups are sometimes treated as adaptive units whose organization cannot be reduced to individual interactions. This group-level view is opposed by a more individualistic one that treats social organization as a byproduct of self-interest. According to biologists, group-level adaptations can evolve only by a process of natural selection at the group level. Most biologists rejected group selection as an important evolutionary force during the 1960s and 1970s but a positive literature began to grow during the 1970s and is rapidly expanding today. We review this recent literature and its implications for human evolutionary biology. We show that the rejection of group selection was based on a misplaced emphasis on genes as “replicators” which is in fact irrelevant to the question of whether groups can be like individuals in their functional organization. The fundamental question is whether social groups and other higher-level entities can be “vehicles” of selection. When this elementary fact is recognized, group selection emerges as an important force in nature and what seem to be competing theories, such as kin selection and reciprocity, reappear as special cases of group selection. The result is a unified theory of natural selection that operates on a nested hierarchy of units.

The vehicle-based theory makes it clear that group selection is an important force to consider in human evolution. Humans can facultatively span the full range from self-interested individuals to “organs” of group-level “organisms.” Human behavior not only reflects the balance between levels of selection but it can also alter the balance through the construction of social structures that have the effect of reducing fitness differences within groups, concentrating natural selection (and functional organization) at the group level. These social structures and the cognitive abilities that produce them allow group selection to be important even among large groups of unrelated individuals.

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.

MHC polymorphism under host-pathogen coevolution

The genes encoding major histocompatibility (MHC) molecules are among the most polymorphic genes known for vertebrates. Since MHC molecules play an important role in the induction of immune responses, the evolution of MHC polymorphism is often explained in terms of increased protection of hosts against pathogens. Two selective pressures that are thought to be involved are (1) selection favoring MHC heterozygous hosts, and (2) selection for rare MHC alleles by host-pathogen coevolution. We have developed a computer simulation of coevolving hosts and pathogens to study the relative impact of these two mechanisms on the evolution of MHC polymorphism. We found that heterozygote advantage per se is insufficient to explain the high degree of polymorphism at the MHC, even in very large host populations. Host-pathogen coevolution, on the other hand, can easily account for realistic polymorphisms of more than 50 alleles per MHC locus. Since evolving pathogens mainly evade presentation by the most common MHC alleles in the host population, they provide a selective pressure for a large variety of rare MHC alleles. Provided that the host population is sufficiently large, a large set of MHC alleles can persist over many host generations under host-pathogen coevolution, despite the fact that allele frequencies continuously change.

Genomic structure and function in the MHC

The major histocompatibility complex (MHC) is a fascinating region of the human genome. More is known about this 4 Mb of DNA (0.1% of the genome) on the short arm of chromosome 6 than about any other region of similar size. Among the 80 or so MHC genes found so far are several clusters with related functions in antigen processing and presentation. In addition to its importance in immunology, the MHC is a useful model for investigating gene organization, polymorphism, linkage disequilibrium and recombination. A large number of diseases, many of the autoimmune type, are associated with the region.

Diallelic polymorphism may explain variations of the blood concentration of mannan-binding protein in Eskimos, but not in black Africans

Mannan-binding protein (MBP) is a lectin which, upon binding to certain carbohydrates, activates the classical pathway of complement without the involvement of antibody or C1q. Deficiency of the MBP is associated with an opsonic defect and recurrent infections during early life. An amino acid substitution in the exon 1 at codon 54 in the MBP gene (GGC [glycine] to GAC [aspartic acid]) has been shown to be closely associated with low MBP concentration in Caucasoids. The gene frequency of the mutant allele in this population has been estimated at 0.13. In the study described here, we investigated the association between the mutant allele and MBP protein concentration in Eskimos from East-Greenland and black Africans from the Baringo District in Kenya. The frequency of the GAC allele was identical in Eskimos and Caucasoids (0.13). No overlap with regard to MBP concentration between the genotypes was found in the Eskimos. In contrast, the Africans revealed a low frequency of the GAC allele (0.009). However, the median MBP protein concentration was approximately 5 times lower among the Africans than the Eskimos. In 12.6% of the Africans and in 2.5% of the Eskimos, MBP was undetectable. Thus, MBP deficiency is the most frequent immunodeficiency so far described. The high prevalence of MBP deficiency among healthy individuals indicates that MBP deficiency also confers some selective advantages. We advance the hypothesis that MBP deficiency is maintained in populations because MBP deficiency decreases the infectivity of some intracellular micro-organisms which are dependent on opsonization.

Unusual HLA-B alleles in two tribes of Brazilian Indians

The Kaingang and Guarani are culturally and linguistically distinct tribes of southern Brazil. Like all Amerindian groups they show limited HLA polymorphism, which probably reflects the small founder populations that colonized America by overland migration from Asia 11,000-40,000 years ago. We find the nucleotide sequences of HLA-B alleles from the Kaingang and Guarani to be distinct from those characterized in caucasian, oriental and other populations. By comparison, the HLA-A and C alleles are familiar. These results and those reported in the accompanying paper on the Waorani of Ecuador reveal that a marked evolution of HLA-B has occurred since humans first entered South America. New alleles have been formed through recombination between pre-existing alleles, not by point mutation, giving rise to distinctive diversification of HLA-B in different South American Indian tribes.