Strong balancing selection at HLA loci: evidence from segregation in South Amerindian families

Admixture has obscured signals of historical hard sweeps in humans

The role of natural selection in shaping biological diversity is an area of intense interest in modern biology. To date, studies of positive selection have primarily relied on genomic datasets from contemporary populations, which are susceptible to confounding factors associated with complex and often unknown aspects of population history. In particular, admixture between diverged populations can distort or hide prior selection events in modern genomes, though this process is not explicitly accounted for in most selection studies despite its apparent ubiquity in humans and other species. Through analyses of ancient and modern human genomes, we show that previously reported Holocene-era admixture has masked more than 50 historic hard sweeps in modern European genomes. Our results imply that this canonical mode of selection has probably been underappreciated in the evolutionary history of humans and suggest that our current understanding of the tempo and mode of selection in natural populations may be inaccurate.

Challenging ancient DNA results about putative HLA protection or susceptibility to Yersina pestis

In a recent article, Immel et al. (Immel A, Key FM, Szolek A, Barquera R, Robinson MK, Harrison GF, Palmer WH, Spyrou MA, Susat J, Krause-Kyora B, et al. 2021. Analysis of genomic DNA from medieval plague victims suggests long-term effect of Yersinia pestis on human immunity genes. Mol Biol Evol. 38:4059–4076) extracted DNA from 36 individuals dead from plague in Ellwangen, Southern Germany, during the 16th century. By comparing their human leukocyte antigen (HLA) genotypes with those of 50 present-day Ellwangen inhabitants, the authors reported a significant decrease of HLA-B*51:01 and HLA-C*06:02 and a significant increase of HLA-DRB1*13:01/13:02 frequencies from ancient to modern populations. After comparing these frequencies with a larger sample of 8,862 modern Germans and performing simulations of natural selection, they concluded that these changes had been driven by natural selection. In an attempt to provide more evidence on such stimulating results, we explored the HLA frequency patterns over all of Europe, we predicted binding affinities of HLA-B/C/DRB1 alleles to 106,515 Yersinia pestis-derived peptides, and we performed forward simulations of HLA genetic profiles under neutrality. Our analyses do not sustain the conclusions of HLA protection or susceptibility to plague based on ancient DNA.

Sequencing of 53,831 diverse genomes from the NHLBI TOPMed Program

The Trans-Omics for Precision Medicine (TOPMed) programme seeks to elucidate the genetic architecture and biology of heart, lung, blood and sleep disorders, with the ultimate goal of improving diagnosis, treatment and prevention of these diseases. The initial phases of the programme focused on whole-genome sequencing of individuals with rich phenotypic data and diverse backgrounds. Here we describe the TOPMed goals and design as well as the available resources and early insights obtained from the sequence data. The resources include a variant browser, a genotype imputation server, and genomic and phenotypic data that are available through dbGaP (Database of Genotypes and Phenotypes)1. In the first 53,831 TOPMed samples, we detected more than 400 million single-nucleotide and insertion or deletion variants after alignment with the reference genome. Additional previously undescribed variants were detected through assembly of unmapped reads and customized analysis in highly variable loci. Among the more than 400 million detected variants, 97% have frequencies of less than 1% and 46% are singletons that are present in only one individual (53% among unrelated individuals). These rare variants provide insights into mutational processes and recent human evolutionary history. The extensive catalogue of genetic variation in TOPMed studies provides unique opportunities for exploring the contributions of rare and noncoding sequence variants to phenotypic variation. Furthermore, combining TOPMed haplotypes with modern imputation methods improves the power and reach of genome-wide association studies to include variants down to a frequency of approximately 0.01%.

Fetal HLA-G mediated immune tolerance and interferon response in preeclampsia

BACKGROUND

Fetal immune tolerance is crucial for pregnancy success. We studied the link between preeclampsia, a severe pregnancy disorder with uncertain pathogenesis, and fetal human leukocyte antigen G (HLA-G) and other genes regulating maternal immune responses.

METHODS

We assessed sex ratios and regulatory HLA-G haplotypes in population cohorts and series of preeclampsia and stillbirth. We studied placental mRNA expression of 136 genes by sequencing and HLA-G and interferon alpha (IFNα) protein expression by immunohistochemistry.

FINDINGS

We found underrepresentation of males in preeclamptic births, especially those delivered preterm or small for gestational age. Balancing selection at HLA-G associated with the sex ratio, stillbirth, and preeclampsia. We observed downregulation of HLA-G, its receptors, and many other tolerogenic genes, and marked upregulation of IFNA1 in preeclamptic placentas.

INTERPRETATION

These findings indicate that an evolutionary trade-off between immune tolerance and protection against infections at the maternal-fetal interface promotes genetic diversity in fetal HLA-G, thereby affecting survival, preeclampsia, and sex ratio. We highlight IFNA1 as a potential mediator of preeclampsia and a target for therapeutic trials.

FUNDING

Finnish Medical Foundation, Päivikki and Sakari Sohlberg Foundation, Karolinska Institutet Research Foundation, Scandinavia-Japan Sasakawa Foundation, Japan Eye Bank Association, Astellas Foundation for Research on Metabolic Disorders, Japan Society for the Promotion of Science, Knut and Alice Wallenberg Foundation, Swedish Research Council, Medical Society Liv och Hälsa, Sigrid Jusélius Foundation, Helsinki University Hospital and University of Helsinki, Jane and Aatos Erkko Foundation, Academy of Finland, Finska Läkaresällskapet, Novo Nordisk Foundation, Finnish Foundation for Pediatric Research, and Emil Aaltonen Foundation.

Divergent Allele Advantage Provides a Quantitative Model for Maintaining Alleles with a Wide Range of Intrinsic Merits

The Major Histocompatibility Complex (MHC) is the most genetically diverse region of the genome in most vertebrates. Some form of balancing selection is necessary to account for the extreme diversity, but the precise mechanism of balancing selection is unknown. Due to the way MHC molecules determine immune recognition, overdominance (also referred to as heterozygote advantage) has been suggested as the main driving force behind this unrivalled diversity. However, both theoretical results and simulation models have shown that overdominance in its classical form cannot maintain large numbers of alleles unless all alleles confer unrealistically similar levels of fitness. There is increasing evidence that heterozygotes containing genetically divergent alleles allow for broader antigen presentation to immune cells, providing a selective mechanism for MHC polymorphism. By framing competing models of overdominance within a general framework, we show that a model based on Divergent Allele Advantage (DAA) provides a superior mechanism for maintaining alleles with a wide range of intrinsic merits, as intrinsically less-fit MHC alleles that are more divergent can survive under DAA. Specifically, our results demonstrate that a quantitative mechanism built from the DAA hypothesis is able to maintain polymorphism in the MHC. Applying such a model to both livestock breeding and conservation could provide a better way of identifying superior heterozygotes, and quantifying the advantages of genetic diversity at the MHC.

HLA-A promoter, coding, and 3'UTR sequences in a Brazilian cohort, and their evolutionary aspects

HLA-A is the second most polymorphic locus of the human leucocyte antigen (HLA) complex encoding a key molecule for antigen presentation and NK cell modulation. Many studies have evaluated HLA-A variability in worldwide populations, focusing mainly on exons, but the regulatory segments have been poorly characterized. HLA-A variability is particularly high in the segment encoding the peptide-binding groove (exons 2 and 3), which is related to the antigen presentation function and the balancing selection in these segments. Here we evaluate the genetic diversity of the HLA-A gene considering a continuous segment encompassing the extended promoter (1.5 kb upstream of the first translated ATG), all exons and introns, and the entire 3' untranslated region, by using massively parallel sequencing. To achieve this goal, we used a freely available bioinformatics workflow that optimizes read mapping for HLA genes and defines complete sequences using either the phase among variable sites directly observed in sequencing data and probabilistic models. The HLA-A variability detected in a highly admixed population sample from Brazil shows that the HLA-A regulatory segments present few, but divergent sequences. The regulatory segments are in close association with the coding alleles. Both exons and introns are highly variable. Moreover, patterns of molecular diversity suggest that the promoter, in addition to the coding region, might be under the same selective pressure, but a different scenario arises when it comes to exon 4 and the 3'UTR segment.

Analysis of type 2 diabetes and obesity genetic variants in Mexican Pima Indians: Marked allelic differentiation among Amerindians at HLA

Prevalence of diabetes and obesity in Mexican Pima Indians is low, while prevalence in US Pima Indians is high. Although lifestyle likely accounts for much of the difference, the role of genetic factors is not well explored. To examine this, we genotyped 359 single nucleotide polymorphisms, including established type 2 diabetes and obesity variants from genome-wide association studies (GWAS) and 96 random markers, in 342 Mexican Pimas. A multimarker risk score of obesity variants was associated with body mass index (BMI; β = 0.81 kg/m² per SD, P = 0.0066). The mean value of the score was lower in Mexican Pimas than in US Pimas (P = 4.3 × 10^-11 ), and differences in allele frequencies at established loci could account for approximately 7% of the population difference in BMI; however, the difference in risk scores was consistent with evolutionary neutrality given genetic distance. To identify loci potentially under recent natural selection, allele frequencies at 283 variants were compared between US and Mexican Pimas, accounting for genetic distance. The largest differences were seen at HLA markers (e.g., rs9271720, difference = 0.75, P = 8.7 × 10^-9 ); genetic distances at HLA were greater than at random markers (P = 1.6 × 10^-46 ). Analyses of GWAS data in 937 US Pimas also showed sharing of alleles identical by descent at HLA that exceeds its genomic expectation (P = 7.0 × 10^-10 ). These results suggest that, in addition to the widely recognized balancing selection at HLA, recent directional selection may also occur, resulting in marked allelic differentiation between closely related populations.

Population Genetics and Natural Selection in Rheumatic Disease

Human genetic diversity is the result of population genetic forces. This genetic variation influences disease risk and contributes to health disparities. Natural selection is an important influence on human genetic variation. Because immune and inflammatory function genes are enriched for signals of positive selection, the prevalence of rheumatic disease-risk alleles seen in different populations is partially the result of differing selective pressures (eg, due to pathogens). This review summarizes the genetic regions associated with susceptibility to different rheumatic diseases and concomitant evidence for natural selection, including known agents of selection exerting selective pressure in these regions.

Inference from the stationary distribution of allele frequencies in a family of Wright-Fisher models with two levels of genetic variability

The distribution of allele frequencies obtained from diffusion approximations to Wright-Fisher models is useful in developing intuition about the population level effects of evolutionary processes. The statistical properties of the stationary distributions of K-allele models have been extensively studied under neutrality or under selection. Here, we introduce a new family of Wright-Fisher models in which there are two hierarchical levels of genetic variability. The genotypes composed of alleles differing from each other at the selected level have fitness differences with respect to each other and evolve under selection. The genotypes composed of alleles differing from each other only at the neutral level have the same fitness and evolve under neutrality. We show that with an appropriate scaling of the mutation parameter with respect to the number of alleles at each level, the frequencies of alleles at the selected and the neutral level are conditionally independent of each other, conditional on knowing the number of alleles at all levels. This conditional independence allows us to simulate from the joint stationary distribution of the allele frequencies. We use these simulated frequencies to perform inference on parameters of the model with two levels of genetic variability using Approximate Bayesian Computation.

Kourami: graph-guided assembly for novel human leukocyte antigen allele discovery

Accurate typing of human leukocyte antigen (HLA) is important because HLA genes play important roles in immune responses and disease genesis. Previously available computational methods are database-matching approaches and their outputs are inherently limited by the completeness of already known types, making them unsuitable for discovery of novel alleles. We have developed a graph-guided assembly technique for classical HLA genes, which can construct allele sequences given high-coverage whole-genome sequencing data. Our method delivers highly accurate HLA typing, comparable to the current state-of-the-art methods. Using various data, we also demonstrate that our method can type novel alleles.

Balancing selection on the number of repeats in the ribosomal intergenic spacer present in naturally occurring yellow perch (Perca flavescens) populations

The ribosomal intergenic spacer (IGS), responsible for the rate of transcription of rRNA genes, is associated with the growth and fecundity of individuals. A previous study of IGS length variants in a yellow perch (Perca flavescens) population revealed the presence of two predominant alleles differing by 1 kb due to variation in the number of repeat units. This study aims to assess whether length variation of IGS is the result of selection in natural populations. Length variation of IGS and 11 neutral microsatellite loci were assessed in geographically distant yellow perch populations. Most populations displayed the very same IGS alleles; they did not differ in frequencies among populations and the F_ST was not significantly different from zero. In contrast, diversity at microsatellite loci was high and differed among populations (F_ST = 0.18). Selection test based on F_ST identified IGS as a significant outlier from neutral expectations for population differentiation. Heterozygote excess was also detected in one specific cohort, suggesting temporal variation in the selection regime. While the exact mechanism remains to be specified, together the results of this study support the contention that balancing selection is acting to maintain two distinct IGS alleles in natural fish populations.

No evidence of prenatal diversifying selection at locus or supertype levels in the dog MHC class II loci

BACKGROUND

Despite decades of studying, the mechanisms maintaining high diversity in the genes of the Major Histocompatibility Complex (MHC) are still puzzling scientists. In addition to pathogen recognition and other functions, MHC molecules may act prenatally in mate choice and in maternal-foetal interactions. These interactions are potential selective mechanisms that increase genetic diversity in the MHC. During pregnancy, immune response has a dual role: the foetus represents foreign tissue compared to mother, but histo-incompatibility is required for successful pregnancy. We have studied the prenatal selection in MHC class II loci (DLA-DQA1, DLA-DQB1 and DLA-DRB1) in domestic dogs by comparing the observed and expected offspring genotype proportions in 110 dog families. Several potential selection targets were addressed, including the peptide-binding site, the MHC locus, three-locus haplotype and supertype levels. For the supertype analysis, the first canine supertype classification was created based on in silico analysis of peptide-binding amino-acid polymorphism.

RESULTS

In most loci and levels, no deviation from the expected genotype frequencies was observed. However, one peptide-binding site in DLA-DRB1 had an excess of heterozygotes among the offspring. In addition, if the father shared a DLA-DRB1 allele with the mother, that allele was inherited by the offspring more frequently than expected, suggesting the selective advantage of a histo-compatible foetus, in contrast to our expectations.

CONCLUSIONS

We conclude that there is some evidence of post-copulatory selection at nucleotide site level in the MHC loci of pet dogs. But due to no indication of selection at locus, three-locus, or supertype levels, we estimated that the prenatal selection coefficient is less than 0.3 in domestic dogs and very likely other factors are more important in maintaining the genetic diversity in MHC loci.

An Ancient Fecundability-Associated Polymorphism Switches a Repressor into an Enhancer of Endometrial TAP2 Expression

Variation in female reproductive traits, such as fertility, fecundity, and fecundability, is heritable in humans, but identifying and functionally characterizing genetic variants associated with these traits has been challenging. Here, we explore the functional significance and evolutionary history of a T/C polymorphism of SNP rs2071473, which we have previously shown is an eQTL for TAP2 and significantly associated with fecundability (time to pregnancy). We replicated the association between the rs2071473 genotype and TAP2 expression by using GTEx data and demonstrated that TAP2 is expressed by decidual stromal cells at the maternal-fetal interface. Next, we showed that rs2071473 is located within a progesterone-responsive cis-regulatory element that functions as a repressor with the T allele and an enhancer with the C allele. Remarkably, we found that this polymorphism arose before the divergence of modern and archaic humans, segregates at intermediate to high frequencies across human populations, and has genetic signatures of long-term balancing selection. This variant has also previously been identified in genome-wide association studies of immune-related disease, suggesting that both alleles are maintained as a result of antagonistic pleiotropy.

Genetics of autoimmune diseases: insights from population genetics

Human genetic diversity is the result of population genetic forces. This genetic variation influences disease risk and contributes to health disparities. Autoimmune diseases (ADs) are a family of complex heterogeneous disorders with similar underlying mechanisms characterized by immune responses against self. Collectively, ADs are common, exhibit gender and ethnic disparities, and increasing incidence. As natural selection is an important influence on human genetic variation, and immune function genes are enriched for signals of positive selection, it is thought that the prevalence of AD risk alleles seen in different population is partially the result of differing selective pressures (for example, due to pathogens). With the advent of high-throughput technologies, new analytical methodologies and large-scale projects, evidence for the role of natural selection in contributing to the heritable component of ADs keeps growing. This review summarizes the genetic regions associated with susceptibility to different ADs and concomitant evidence for selection, including known agents of selection exerting selective pressure in these regions. Examples of specific adaptive variants with phenotypic effects are included as an evidence of natural selection increasing AD susceptibility. Many of the complexities of gene effects in different ADs can be explained by population genetics phenomena. Integrating AD susceptibility studies with population genetics to investigate how natural selection has contributed to genetic variation that influences disease risk will help to identify functional variants and elucidate biological mechanisms. As such, the study of population genetics in human population holds untapped potential for elucidating the genetic causes of human disease and more rapidly focusing to personalized medicine.

Genes associated with SLE are targets of recent positive selection

The reasons for the ethnic disparities in the prevalence of systemic lupus erythematosus (SLE) and the relative high frequency of SLE risk alleles in the population are not fully understood. Population genetic factors such as natural selection alter allele frequencies over generations and may help explain the persistence of such common risk variants in the population and the differential risk of SLE. In order to better understand the genetic basis of SLE that might be due to natural selection, a total of 74 genomic regions with compelling evidence for association with SLE were tested for evidence of recent positive selection in the HapMap and HGDP populations, using population differentiation, allele frequency, and haplotype-based tests. Consistent signs of positive selection across different studies and statistical methods were observed at several SLE-associated loci, including PTPN22, TNFSF4, TET3-DGUOK, TNIP1, UHRF1BP1, BLK, and ITGAM genes. This study is the first to evaluate and report that several SLE-associated regions show signs of positive natural selection. These results provide corroborating evidence in support of recent positive selection as one mechanism underlying the elevated population frequency of SLE risk loci and supports future research that integrates signals of natural selection to help identify functional SLE risk alleles.

A signature of balancing selection in the region upstream to the human UGT2B4 gene and implications for breast cancer risk

UDP-glucuronosyltransferase 2 family, polypeptide B4 (UGT2B4) is an important metabolizing enzyme involved in the clearance of many xenobiotics and endogenous substrates, especially steroid hormones and bile acids. The HapMap data show that numerous SNPs upstream of UGT2B4 are in near-perfect linkage disequilibrium with each other and occur at intermediate frequency, indicating that this region might contain a target of natural selection. To investigate this possibility, we chose three regions (4.8 kb in total) for resequencing and observed a striking excess of intermediate-frequency alleles that define two major haplotypes separated by many mutation events and with little differentiation across populations, thus suggesting that the variation pattern upstream UGT2B4 is highly unusual and may be the result of balancing selection. We propose that this pattern is due to the maintenance of a regulatory polymorphism involved in the fine tuning of UGT2B4 expression so that heterozygous genotypes result in optimal enzyme levels. Considering the important role of steroid hormones in breast cancer susceptibility, we hypothesized that variation in this region could predispose to breast cancer. To test this hypothesis, we genotyped tag SNP rs13129471 in 1,261 patients and 825 normal women of African ancestry from three populations. The frequency comparison indicated that rs13129471 was significantly associated with breast cancer after adjusting for ethnicity [P = 0.003; heterozygous odds ratio (OR) 1.02, 95% confidence interval (CI) 0.81-1.28; homozygous OR 1.50, 95% CI 1.15-1.95]. Our results provide new insights into UGT2B4 sequence variation and indicate that a signal of natural selection may lead to the identification of disease susceptibility variants.

In the nose of the beholder: are olfactory influences on human mate choice driven by variation in immune system genes or sex hormone levels?

The human leukocyte antigen (HLA) is the most polymorphic region of the genome, coding for proteins that mediate human immune response. This polymorphism may be maintained by balancing selection and certain populations show deviations from expected gene frequencies. Supporting this hypothesis, studies into olfactory preferences have suggested that females prefer the scent of males with dissimilar HLA to their own. However, it has also been proposed that androstenones play a role in female mate choice, and as these molecules inhibit the immune system, this has implications for the theory of HLA-driven mate preference. This review will critically analyze the findings of studies investigating olfactory preference in humans, and their implications for these two contrasting theories of mate choice.

Estimation of selection intensity under overdominance by Bayesian methods

A balanced pattern in the allele frequencies of polymorphic loci is a potential sign of selection, particularly of overdominance. Although this type of selection is of some interest in population genetics, there exists no likelihood based approaches specifically tailored to make inference on selection intensity. To fill this gap, we present Bayesian methods to estimate selection intensity under k-allele models with overdominance. Our model allows for an arbitrary number of loci and alleles within a locus. The neutral and selected variability within each locus are modeled with corresponding k-allele models. To estimate the posterior distribution of the mean selection intensity in a multilocus region, a hierarchical setup between loci is used. The methods are demonstrated with data at the Human Leukocyte Antigen loci from world-wide populations.

Evolutionary history of an MHC gene in two leporid species: characterisation of Mhc-DQA in the European brown hare and comparison with the European rabbit

We surveyed the genetic diversity of the expressed major histocompatibility complex class II DQA locus in natural populations of European brown hares, Lepus europaeus, from Austria and Belgium (267 individuals in total). Based on cDNA sequences, we designed hare-specific primers to amplify the highly variable second exon of the DQA gene. Using cloning-sequencing methodology and capillary electrophoresis single-strand conformation polymorphism, we found ten alleles of the DQA exon 2 locus across these two European regions, of which eight are described for the first time. To search for signals of selection and recombination in the evolution of the DQA gene within the leporids, we augmented our sample with orthologous DQA alleles from the European rabbit, Oryctolagus cuniculus, in order to carry out a species level, species pairwise comparison. We found evidence of recombination in the history of the DQA sequences in leporids with some recombinant alleles bridging the species divide. In both species, selection on peptide binding site codons can be detected, though stronger for the rabbit. This result suggests that there may be a differential selection pressure in the deeper evolutionary history of these two species due to differences in several demographic and ecological traits likely subjecting them to differential selection by parasites. Finally, evolutionary relationships show a widespread and statistically significant intermingling of alleles from the two species. The many macroparasites shared between hares and rabbits may explain this pattern of trans-species polymorphism.

Human endogenous retrovirus (HERVK9) structural polymorphism with haplotypic HLA-A allelic associations

The frequency and HLA-A allelic associations of a HERVK9 DNA structural polymorphism located in close proximity to the highly polymorphic HLA-A gene within the major histocompatibility complex (MHC) genomic region were determined in Japanese, African Americans, and Australian Caucasians to better understand its human population evolutionary history. The HERVK9 insertion or deletion was detected as a 3' LTR or a solo LTR, respectively, by separate PCR assays. The average insertion frequency of the HERVK9.HG was significantly different (P < 1.083e(-6)) between the Japanese (0.59) and the African Americans (0.34) or Australian Caucasians (0.37). LD analysis predicted a highly significant (P < 1.0e(-5)) linkage between the HLA-A and HERVK9 alleles, probably as a result of hitchhiking (linkage). Evolutionary time estimates of the solo, 5' and 3' LTR nucleotide sequence divergences suggest that the HERVK9 was inserted 17.3 MYA with the first structural deletion occurring 15.1 MYA. The LTR/HLA-A haplotypes appear to have been formed mostly during the past 3.9 MY. The HERVK9 insertion and deletion, detected by a simple and economical PCR method, is an informative genetic and evolutionary marker for the study of HLA-A haplotype variations, human migration, the origins of contemporary populations, and the possibility of disease associations.

A new test for detecting ongoing selection

Much effort has been made to search for signatures of past natural selection in DNA sequences. However, currently acting selection is rarely detected in natural populations because of its rarity, low detection power of available methods, or both. Here, we develop a new test to detect viability selection over a single generation. In this test, one specific type of chromosomes is chosen as a reference, while all other chromosomes are designated as "focal". The test compares measures of variation between two groups of "focal" chromosomes: those found in reference/focal heterozygous individuals and those found in focal/focal homozygous individuals. In the absence of selection, we do not expect differences between these two groups as long as mating is random. On the other hand, currently acting selection can cause differences in some measures of variation. We applied this test to typing data for In(2L)t inversion polymorphism in a Drosophila melanogaster population, using "standard" (non-inverted) chromosomes as the focal class. Although the frequencies of In(2L)t and standard chromosomes did not deviate from the Hardy-Weinberg equilibrium, we found differences in allele frequency and the number of haplotypes between the two groups of standard chromosomes. This new test, in conjunction with the Hardy-Weinberg test, may shed light on how often strong selection is operating in extant populations.

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.

Amazonian Amerindians exhibit high variability of KIR profiles

Natural killer cell immunoglobulin-like receptors (KIRs) mediate cell lysis through the recognition of human leukocyte antigen class I complexes in target cells, playing an important role in innate immune response. In this context, disease-based selective pressures could be relevant, leaving signatures detected by population studies. However, most population studies on KIR variability have focused on Europe and Asia, while Americas, Oceania, and Africa remain poorly studied. The aim of this study was to analyze the variability of KIR genes in Amerindian tribes from the Amazon region to infer about their evolutionary history. KIR profiles were estimated in 40 individuals from six Amazonian Amerindian tribes using single specific primer polymerase chain reaction. Twenty-five different profiles were identified, and surprisingly, the haplogroup A frequency was the lowest observed in human populations (16%). Results showed also that KIR variability was higher in this group in contrast to Venezuelan Amerindians. Principal components analysis evidenced that Amerindians formed a separated group from other worldwide populations and showed a higher intraethnic differentiation in comparison to other ethnic groups. Such pattern may reflect small effective size and intense genetic drift. However, because of the role of KIR in immune response, selective pressures cannot be entirely ruled out.

Selection on MHC-linked microsatellite loci in sheep populations

Microsatellites within the major histocompatibility complex (MHC) region have received increasing attention as proxy measures of the level of polymorphism at the Mhc genes themselves. We assessed the diversity of microsatellite loci within or in close proximity of the Mhc genes in several breeds of domestic sheep (Ovis aries) and the wild Mouflon (Ovis orientalis musimon). This was compared to variation at other microsatellite loci scattered throughout the sheep genome. Significantly higher number of alleles were observed at the MHC microsatellites. The sheep breeds studied fell into high- and low-diversity group. This grouping is not related to the agricultural use of the breeds, whether for milk, meat or wool. It is, however, correlated with the geographic origins of the breeds. Southern breeds are genetically more diverse than northern breeds. The observed heterozygosity was in most cases lower than Hardy-Weinberg expectations. The potential impact of selective breeding by man on this is discussed. Neutrality tests indicated that for most of the breeds, the distribution of alleles at the MHC-linked microsatellites are more even than would be expected if the genes were neutral and sampled from populations under drift-mutation equilibrium. Hitchhiking due to tight linkage with alleles at the MHC loci that are under balancing selection is proposed as a possible explanation for this pattern.

Genomic variation in rice: genesis of highly polymorphic linkage blocks during domestication

Genomic regions that are unusually divergent between closely related species or racial groups can be particularly informative about the process of speciation or the operation of natural selection. The two sequenced genomes of cultivated Asian rice, Oryza sativa, reveal that at least 6% of the genomes are unusually divergent. Sequencing of ten unlinked loci from the highly divergent regions consistently identified two highly divergent haplotypes with each locus in nearly complete linkage disequilibrium among 25 O. sativa cultivars and 35 lines from six wild species. The existence of two highly divergent haplotypes in high divergence regions in species from all geographical areas (Africa, Asia, and Oceania) was in contrast to the low polymorphism and low linkage disequilibrium that were observed in other parts of the genome, represented by ten reference loci. While several natural processes are likely to contribute to this pattern of genomic variation, domestication may have greatly exaggerated the trend. In this hypothesis, divergent haplotypes that were adapted to different geographical and ecological environments migrated along with humans during the development of domesticated varieties. If true, these high divergence regions of the genome would be enriched for loci that contribute to the enormous range of phenotypic variation observed among domesticated breeds.

The separate and combined effects of MHC genotype, parasite clone, and host gender on the course of malaria in mice

Background

The link between host MHC (major histocompatibility complex) genotype and malaria is largely based on correlative data with little or no experimental control of potential confounding factors. We used an experimental mouse model to test for main effects of MHC-haplotypes, MHC heterozygosity, and MHC × parasite clone interactions. We experimentally infected MHC-congenic mice (F2 segregants, homo- and heterozygotes, males and females) with one of two clones of Plasmodium chabaudi and recorded disease progression.

Results

We found that MHC haplotype and parasite clone each have a significant influence on the course of the disease, but there was no significant host genotype by parasite genotype interaction. We found no evidence for overdominance nor any other sort of heterozygote advantage or disadvantage.

Conclusion

When tested under experimental conditions, variation in the MHC can significantly influence the course of malaria. However, MHC heterozygote advantage through overdominance or dominance of resistance cannot be assumed in the case of single-strain infections. Future studies might focus on the interaction between MHC heterozygosity and multiple-clone infections.

Animal breeding and disease

Single-locus disorders in domesticated animals were among the first Mendelian traits to be documented after the rediscovery of Mendelism, and to be included in early linkage maps. The use of linkage maps and (increasingly) comparative genomics has been central to the identification of the causative gene for single-locus disorders of considerable practical importance. The 'score-card' in domestic animals is now more than 100 disorders for which the molecular lesion has been identified and hence for which a DNA test is available. Because of the limited lifespan of any such test, a cost-effective and hence popular means of protecting the intellectual property inherent in a DNA test is not to publish the discovery. While understandable, this practice creates a disconcerting precedent. For multifactorial disorders that are scored on an all-or-none basis or into many classes, the effectiveness of control schemes could be greatly enhanced by selection on estimated breeding values for liability. Genetic variation for resistance to pathogens and parasites is ubiquitous. Selection for resistance can therefore be successful. Because of the technical and welfare challenges inherent in the requirement to expose animals to pathogens or parasites in order to be able to select for resistance, there is a very active search for DNA markers for resistance. The first practical fruits of this research were seen in 2002, with the launch of a national scrapie control programme in the UK.

Genetic variability in a genomic region with long-range linkage disequilibrium reveals traces of a bottleneck in the history of the European population

The inference of the demographic history of populations from genetic variability data is not only of academic interest. It also provides background information for the identification of genes which may have played a role in human evolution or in the aetiology of human disease. To obtain a clear picture of this background, it is necessary to compare data obtained from a number of genomic loci. Due to its very low recombination rate, the NF1 gene region can be regarded as a further suitable locus. A combined resequencing and SNP typing project in a European population disclosed the presence of only two well separated subgroups of NF1 sequences. Statistical analysis revealed a bimodal distribution of the pairwise differences, a positive value of Tajima's D and a TMRCA of 700,000 years for the whole sample, and pairwise differences indicative for a growing population and TMRCAs of 130,000 to 150,000 years for the subgroups. Together, the data lead to a model that the recent European population went through a bottleneck during the last 150,000 years of its history. Regarding the given timeframe, this bottleneck could either reflect a speciation event which led to the anatomically modern human (AMH), or a severe reduction of the population size during the emigration of AMHs out of Africa or the immigration into Europe.

Detecting site-specific physicochemical selective pressures: applications to the Class I HLA of the human major histocompatibility complex and the SRK of the plant sporophytic self-incompatibility system

Models of codon substitution are developed that incorporate physicochemical properties of amino acids. When amino acid sites are inferred to be under positive selection, these models suggest the nature and extent of the physicochemical properties under selection. This is accomplished by first partitioning the codons on the basis of some property of the encoded amino acids. This partition is used to parametrize the rates of property-conserving and property-altering base substitutions at the codon level by means of finite mixtures of Markov models that also account for codon and transition:transversion biases. Here, we apply this method to two positively selected receptors involved in ligand-recognition: the class I alleles of the human major histocompatibility complex (MHC) of known structure and the S-locus receptor kinase (SRK) of the sporophytic self-incompatibility system (SSI) in cruciferous plants (Brassicaceae), whose structure is unknown. Through likelihood ratio tests we demonstrate that at some sites, the positively selected MHC and SRK proteins are under physicochemical selective pressures to alter polarity, volume, polarity and/or volume, and charge to various extents. An empirical Bayes approach is used to identify sites that may be important for ligand recognition in these proteins.

The course of malaria in mice: major histocompatibility complex (MHC) effects, but no general MHC heterozygote advantage in single-strain infections

A general MHC-heterozygote advantage in parasite-infected organisms is often assumed, although there is little experimental evidence for this. We tested the response of MHC-congenic mice (F2 segregants) to malaria and found the course of infection to be significantly influenced by MHC haplotype, parasite strain, and host gender. However, the MHC heterozygotes did worse than expected from the average response of the homozygotes.

Major histocompatibility complex and microsatellite variation in two populations of wild gorillas

In comparison to their close relatives the chimpanzees and humans, very little is known concerning the amount and structure of genetic variation in gorillas. Two species of gorillas are recognized and while the western gorillas number in the tens of thousands, only several hundred representatives of the mountain gorilla subspecies of eastern gorillas survive. To analyse the possible effects of these different population sizes, this study compares the variation observed at microsatellite and major histocompatibility complex (MHC) loci in samples of wild western and mountain gorillas, collected using a sampling scheme that targeted multiple social groups within defined geographical areas. Noninvasive samples proved a viable source of DNA for sequence analysis of the second exon of the DRB loci of the MHC. Observed levels of variation at the MHC locus were similar between the two gorilla species and were comparable to those in other primates. Comparison of results from analysis of variation at multiple microsatellite loci found only a slight reduction in heterozygosity for the mountain gorillas despite the relatively smaller population size.

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.

High MHC diversity maintained by balancing selection in an otherwise genetically monomorphic mammal

The San Nicolas Island fox (Urocyon littoralis dickeyi) is genetically the most monomorphic sexually reproducing animal population yet reported and has no variation in hypervariable genetic markers. Such low levels of variation imply lower resistance to pathogens, reduced fitness, and problems in distinguishing kin from non-kin. In vertebrates, the MHC contains genes that influence disease resistance and kin recognition and may be under intense balancing selection in some populations. Hence, genetic variation at the MHC might persist despite the extreme monomorphism shown by neutral markers. We examine variation of five loci within the MHC of San Nicolas Island foxes and find remarkably high levels of variation. Further, we show by simulation that genetic monomorphism at neutral loci and high MHC variation could arise only through an extreme population bottleneck of <10 individuals, approximately 10-20 generations ago, accompanied by unprecedented selection coefficients of >0.5 on MHC loci. These results support the importance of balancing selection as a mechanism to maintain variation in natural populations and expose the difficulty of using neutral markers as surrogates for variation in fitness-related loci.

Perspective: detecting adaptive molecular polymorphism: lessons from the MHC

In the 1960s, when population geneticists first began to collect data on the amount of genetic variation in natural populations, balancing selection was invoked as a possible explanation for how such high levels of molecular variation are maintained. However, the predictions of the neutral theory of molecular evolution have since become the standard by which cases of balancing selection may be inferred. Here we review the evidence for balancing selection acting on the major histocompatibility complex (MHC) of vertebrates, a genetic system that defies many of the predictions of neutrality. We apply many widely used tests of neutrality to MHC data as a benchmark for assessing the power of these tests. These tests can be categorized as detecting selection in the current generation, over the history of populations, or over the histories of species. We find that selection is not detectable in MHC datasets in every generation, population, or every evolutionary lineage. This suggests either that selection on the MHC is heterogeneous or that many of the current neutrality tests lack sufficient power to detect the selection consistently. Additionally, we identify a potential inference problem associated with several tests of neutrality. We demonstrate that the signals of selection may be generated in a relatively short period of microevolutionary time, yet these signals may take exceptionally long periods of time to be erased in the absence of selection. This is especially true for the neutrality test based on the ratio of nonsynonymous to synonymous substitutions. Inference of the nature of the selection events that create such signals should be approached with caution. However, a combination of tests on different time scales may overcome such problems.

Role of viruses in human evolution

The study of viral molecular genetics has produced a considerable body of research into the sequences and phylogenetic relationships of human and animal viruses. A review of this literature suggests that humans have been afflicted by viruses throughout their evolutionary history, although the number and types have changed. Some viruses show evidence of long-standing intimate relationship and cospeciation with hominids, while others are more recently acquired from other species, including African monkeys and apes while our line was evolving in that continent, and domesticated animals and rodents since the Neolithic. Viral selection for specific resistance polymorphisms is unlikely, but in conjunction with other parasites, viruses have probably contributed to selection pressure maintaining major histocompatibility complex (MHC) diversity and a strong immune response. They may also have played a role in the loss in our lineage of N-glycolylneuraminic acid (Neu5Gc), a cell-surface receptor for many infectious agents. Shared viruses could have affected hominid species diversity both by promoting divergence and by weeding out less resistant host populations, while viruses carried by humans and other animals migrating out of Africa may have contributed to declines in other populations. Endogenous retroviral insertions since the divergence between humans and chimpanzees were capable of directly affecting hominid evolution through changes in gene expression and development.

Balanced polymorphism selected by genetic versus infectious human disease

The polymorphisms within the human genome include several functional variants that cause debilitating inherited diseases. An elevated frequency of some of these deleterious mutations can be explained by a beneficial effect that confers a selective advantage owing to disease resistance in carriers of such mutations during an infectious disease outbreak. We here review plausible examples of balanced functional polymorphisms and their roles in the defense against pathogens. The genome organization of the chemokine receptor and HLA gene clusters and their influence on the HIV/AIDS epidemic provides compelling evidence for the interaction of infectious and genetic diseases in recent human history.

Pathogen resistance and genetic variation at MHC loci

Balancing selection in the form of heterozygote advantage, frequency-dependent selection, or selection that varies in time and/or space, has been proposed to explain the high variation at major histocompatibility complex (MHC) genes. Here the effect of variation of the presence and absence of pathogens over time on genetic variation at multiallelic loci is examined. In the basic model, resistance to each pathogen is conferred by a given allele, and this allele is assumed to be dominant. Given that s is the selective disadvantage for homozygotes (and heterozygotes) without the resistance allele and the proportion of generations, which a pathogen is present, is e, fitnesses for homozygotes become (1 - s)((n - 1)e) and the fitnesses for heterozygotes become (1 - s)((n - 2)e), where n is the number of alleles. In this situation, the conditions for a stable, multiallelic polymorphism are met even though there is no intrinsic heterozygote advantage. The distribution of allele frequencies and consequently heterozygosity are a function of the autocorrelation of the presence of the pathogen in subsequent generations. When there is a positive autocorrelation over generations, the observed heterozygosity is reduced. In addition, the effects of lower levels of selection and dominance and the influence of genetic drift were examined. These effects were compared to the observed heterozygosity for two MHC genes in several South American Indian samples. Overall, resistance conferred by specific alleles to temporally variable pathogens may contribute to the observed polymorphism at MHC genes and other similar host defense loci.

Increased heterozygosity for MHC class II lineages in newborn males

In plants, fungi and marine invertebrates, there are genetic compatibility systems to ensure diversity in the offspring. The importance of genetic compatibility in gametic union and selective abortion in vertebrate animals has also been appreciated recently. There have been suggestions that the major histocompatibility complex (HLA in humans) may be a compatibility system in vertebrates. HLA class II haplotypes often contain a second expressed DRB locus which can be either DRB3, DRB4 or DRB5. These encode the supertypical specificities and mark the ancestral lineages. The members of each lineage have related DNA sequences at the main class II locus HLA-DRB1. We analysed 415 newborns at all expressed DRB loci by PCR analysis to seek evidence for sex-specific prenatal selection events. While there was no significant change in heterozygosity rates between males and females at DRB1, the proportion of males carrying two DRB1 specificities from different ancestral lineages was significantly increased (53.7% in males vs 39.3% in females, P = 0.003). The genotypes consisting of phylogenetically most distinct ones, namely the DRB3 and DRB4 haplotypes, showed the most striking difference between sexes (P = 0.007). These results suggested a more favourable outcome for male concepti heterozygous for supertypical haplotypes. Heterozygosity for most divergent haplotypical families ensures the highest degree of functional heterozygosity at the main HLA class II locus DRB1 while increasing the likelihood of heterozygosity also at other MHC loci. Our observations agree with the previously reported heterozygote excess in male newborn rats and mice. Correlations between MHC class II heterozygosity and advertised male quality in deer and pheasant as well as increased reproductive success in MHC class II heterozygous male macaques are examples of postnatal benefits of heterozygosity in males that may be behind the development of prenatal selection mechanisms. The MHC-mediated prenatal selection of males may also be one of the selective events suggested by the very high primary (male-to-female) sex ratio at fertilization reaching close to unity at birth in humans. These results provide an appealing working hypothesis for further studies in humans and other vertebrates.

Immunoglobulin GM and KM allotypes and vaccine immunity

Immunoglobulin (Ig) GM and KM allotypes are genetic markers of gamma and kappa-type light chains, respectively. The striking qualitative and quantitative differences in the distribution of these determinants among different races raise questions concerning the nature of the evolutionary selective mechanism that maintains this variation. Associations between Ig allotypes and specific antibody responses could be a selective force for the maintenance of various haplotypes and their frequencies. Data from several studies reporting significant associations between certain GM and KM allotypes and immune responsiveness to polysaccharide vaccines and to particular infectious pathogens support this hypothesis. Possible ways in which constant (C)-region allotypes could contribute to the antibody specificity include the following: (i) certain alleles coding for allotypes may be in linkage disequilibrium with particular variable (V)-region determinants associated with immune responsiveness; (ii) they could directly contribute to the formation of specific idiotypes, as shown for the T15 system in mice; and (iii) allotype-associated structural variability in the C-region could modulate the kinetic competence of the antigen binding sites.

Do infectious diseases drive MHC diversity?

The primary function of the major histocompatibility complex (MHC) is to allow the immune system to identify infectious pathogens and eliminate them. Infectious diseases are now thought to be the main selection force that drives and maintains the extraordinary diversity of the MHC.

Congenital adrenal hyperplasia due to 21-hydroxylase deficiency

More than 90% of cases of congenital adrenal hyperplasia (CAH, the inherited inability to synthesize cortisol) are caused by 21-hydroxylase deficiency. Females with severe, classic 21-hydroxylase deficiency are exposed to excess androgens prenatally and are born with virilized external genitalia. Most patients cannot synthesize sufficient aldosterone to maintain sodium balance and may develop potentially fatal "salt wasting" crises if not treated. The disease is caused by mutations in the CYP21 gene encoding the steroid 21-hydroxylase enzyme. More than 90% of these mutations result from intergenic recombinations between CYP21 and the closely linked CYP21P pseudogene. Approximately 20% are gene deletions due to unequal crossing over during meiosis, whereas the remainder are gene conversions--transfers to CYP21 of deleterious mutations normally present in CYP21P. The degree to which each mutation compromises enzymatic activity is strongly correlated with the clinical severity of the disease in patients carrying it. Prenatal diagnosis by direct mutation detection permits prenatal treatment of affected females to minimize genital virilization. Neonatal screening by hormonal methods identifies affected children before salt wasting crises develop, reducing mortality from this condition. Glucocorticoid and mineralocorticoid replacement are the mainstays of treatment, but more rational dosing and additional therapies are being developed.

HLA and HIV-1: heterozygote advantage and B*35-Cw*04 disadvantage

A selective advantage against infectious disease associated with increased heterozygosity at the human major histocompatibility complex [human leukocyte antigen (HLA) class I and class II] is believed to play a major role in maintaining the extraordinary allelic diversity of these genes. Maximum HLA heterozygosity of class I loci (A, B, and C) delayed acquired immunodeficiency syndrome (AIDS) onset among patients infected with human immunodeficiency virus-type 1 (HIV-1), whereas individuals who were homozygous for one or more loci progressed rapidly to AIDS and death. The HLA class I alleles B*35 and Cw*04 were consistently associated with rapid development of AIDS-defining conditions in Caucasians. The extended survival of 28 to 40 percent of HIV-1-infected Caucasian patients who avoided AIDS for ten or more years can be attributed to their being fully heterozygous at HLA class I loci, to their lacking the AIDS-associated alleles B*35 and Cw*04, or to both.