Computer simulation analysis suggests weak balancing selection operative at the MICA locus

A high degree of polymorphism has been reported at the major histocompatibility class I chain-related gene A (MICA) locus, which is located 46 kb away from HLA-Bin the human major histocompatibility complex (MHC) class I region. Although it is known that the polymorphisms at the conventional MHC class I loci have been maintained by balancing selection, it is unclear whether positive natural selection is also operative in maintaining the polymorphism at the MICA locus. In order to explain the degree of polymorphism at the MICA locus, a computer simulation study was carried out. The high degree of polymorphism at the MICA locus (heterozygosity and number of polymorphic residues) could not be explained solely by balancing selection at the HLA-B locus even if no recombination was assumed between MICA and HLA-B. Although there is no definite evidence indicating that balancing selection is operative at the MICA locus, our results suggest that the MICA gene is subject to weak balancing selection.     

Genetics of cell-mediated lympholysis in man.

Cell-mediated lympholysis (CML) was studied in a family containing two siblings in who genetic recombinaiton had occurred in the HLA comples. In one sibling, recombination occurred between the HLA-A locus and the HLA-B locus. In the second sibling recombination occurred between the HLA-B locus and the HLA-D locus. Strong CML activity was generated in mixed lymphocyte cultures (MLC) when stimulator and responder cells differed in HLA-A, B, and D antigens. MLC involving HLA-D differences alone did not generate CML. Weak, but definite CML activity was generated during MLC with cells differing at HLA-A and HLA-B but sharing HLA-D. HLA-B antigens were good targets for lysis in all combinations studied. HLA-A antigens were poor targets in some but not in all combinations. However, combinations where HLA-A antigens seemed to be good targets could have involved HLA-B differences due to polymorphism of HLA-B7 antigens each inherited from a different parent. HLA-D antigens did not serve as targets for lysis. In three cell experiments, excellent CML activity was generated when responder cells were stimulated by HLA-D antigens and by HLA-A and B antigens present on separate stimulator cells.     

Sojourn times and substitution rate at overdominant and linked neutral loci

The sojourn times until fixation of an overdominant allele were investigated based on the diffusion equation. Furthermore, the rate of accumulation of mutations, or the substitution rate, was predicted from the mean extinction time of a common overdominant allele. The substitution rate calculated theoretically agreed well with that determined by computer simulation. Overdominant selection enhances the polymorphism at linked loci, while its effect on the sojourn times and the substitution rate at linked loci has not been studied yet. To solve these problems, a model that assumed two linked loci, each with infinite alleles, was examined by computer simulation. A decrease in the recombination rate between two loci markedly changed the distribution of sojourn times of a neutral allele. Although overdominant selection obviously increased the sojourn times and the polymorphism at a linked locus, the rate of nucleotide substitution at the neutral locus was not influenced significantly even if complete linkage was assumed. These results suggest that, in regions containing overdominant genes, linked neutral loci will exhibit elevated levels of polymorphism, but their rate of molecular evolution remains that predicted by neutral theory.     

Factors affecting polymorphism at microsatellite loci in bread wheat [<Emphasis Type="Italic">Triticum aestivum</Emphasis> (L.) Thell]: effects of mutation processes and physical distance from the centromere

The effects of factors known to influence the level of polymorphism at microsatellite loci were studied using 99 markers and seven lines of bread wheat. Mutational factors as well as indirect selective events shape diversity at these loci. Theory predicts that the selection of favorable alleles should reduce polymorphism at neutral neighboring loci in genomic areas with low recombination rates. In wheat, local recombination rate is positively correlated with physical distance from the centromere. Seventy four loci among the 99 used could be physically located on the chromosome. We studied how the following affected the diversity among a set of inbred lines: the length of the alleles, the motif (CA versus CT), the structure of the loci (perfect versus imperfect) and the chromosomal position of the loci. For each locus, we determined whether the polymorphism observed at a locus was compatible with the Stepwise Mutation Model (SMM) or the Two-Phase Model (TPM). Both the mutation rate and the compatibility with the SMM or the TPM were shown to be variable between loci. Wheat microsatellite loci were found to be more variable when segregating alleles were perfect and had long motifs (composed of many repetitions). Diversity observed at 19 loci was not compatible with the SMM. Loci located in distal regions, with presumably high recombination rates, had longer allele sizes and were more polymorphic than loci located in proximal regions. We conclude that both mutation factors and indirect selective events vary according to the local recombination rate and therefore jointly influence the level of polymorphism at microsatellite loci in wheat.Communicated by J. Dvorak     

The relationship between microsatellite polymorphism and recombination hot spots in the human genome

Although previous studies have failed to detect an association between microsatellite polymorphism and broadscale recombination rates in the human genome, there are several possible reasons why such a relationship could exist. For instance, there might be a direct link if recombination is mutagenic to microsatellite sequences or if polymorphic microsatellites act as recombination signals. Alternatively, recombination could exert an indirect effect by uncoupling of natural selection at linked loci, promoting polymorphism. As recombination is concentrated in narrow hotspot regions in the human genome, we investigated the relationship between microsatellite polymorphism and recombination hot spots. By using data from a common allele frequency database, we found several polymorphism estimates to be similar for hot spots and the genomic average. However, this is likely explained by an ascertainment bias because markers with high polymorphism information content are usually selected for genotyping in human populations and pedigrees. In contrast, by using an unbiased set of shotgun sequence data, we found an excess of microsatellite polymorphism in recombination hot spots of 14%. However, when other genomic variables are taken into account in a generalized model and using wavelet analysis, the effect is no longer detectable and the only firm predictor of microsatellite polymorphism is the incidence of SNPs and indels. One possible neutral explanation to these observations is that there is a common denominator affecting the local rate of mutation in unique as well as in repetitive DNA, for example, base composition.     

Pervasive Hitchhiking at Coding and Regulatory Sites in Humans

Much effort and interest have focused on assessing the importance of natural selection, particularly positive natural selection, in shaping the human genome. Although scans for positive selection have identified candidate loci that may be associated with positive selection in humans, such scans do not indicate whether adaptation is frequent in general in humans. Studies based on the reasoning of the MacDonald–Kreitman test, which, in principle, can be used to evaluate the extent of positive selection, suggested that adaptation is detectable in the human genome but that it is less common than in Drosophila or Escherichia coli. Both positive and purifying natural selection at functional sites should affect levels and patterns of polymorphism at linked nonfunctional sites. Here, we search for these effects by analyzing patterns of neutral polymorphism in humans in relation to the rates of recombination, functional density, and functional divergence with chimpanzees. We find that the levels of neutral polymorphism are lower in the regions of lower recombination and in the regions of higher functional density or divergence. These correlations persist after controlling for the variation in GC content, density of simple repeats, selective constraint, mutation rate, and depth of sequencing coverage. We argue that these results are most plausibly explained by the effects of natural selection at functional sites—either recurrent selective sweeps or background selection—on the levels of linked neutral polymorphism. Natural selection at both coding and regulatory sites appears to affect linked neutral polymorphism, reducing neutral polymorphism by 6% genome-wide and by 11% in the gene-rich half of the human genome. These findings suggest that the effects of natural selection at linked sites cannot be ignored in the study of neutral human polymorphism.     

Molecular evolution near a two-locus balanced polymorphism

Balancing selection at one locus can increase the amount of selectively neutral variation within neighboring genomic regions. Discrete phenotypic polymorphisms studied in natural populations are frequently determined by sets of interacting genes instead of alternative alleles at single loci. We extend coalescent theory to investigate balancing selection on combinations of linked genes. We find that variation at neutral sites is increased across a much larger genomic region relative to the single-locus models: the entire region lying between the two loci in balanced combination is affected to some degree. Epistatic selection maintains these high levels of neutral variation because it directly opposes the homogenizing effect of recombination. The results of the theory are discussed in relation to published gene sequence data, primarily from Drosophila.     

A Test of Neutral Molecular Evolution Based on Nucleotide Data

The neutral theory of molecular evolution predicts that regions of the genome that evolve at high rates, as revealed by interspecific DNA sequence comparisons, will also exhibit high levels of polymorphism within species. We present here a conservative statistical test of this prediction based on a constant-rate neutral model. The test requires data from an interspecific comparison of at least two regions of the genome and data on levels of intraspecific polymorphism in the same regions from at least one species. The model is rejected for data from the region encompassing the Adh locus and the 5' flanking sequence of Drosophila melanogaster and Drosophila sechellia. The data depart from the model in a direction that is consistent with the presence of balanced polymorphism in the coding region.     

Allodeterminants and evolution of a novel HLA-B5 CREG antigen, HLA-B SNA.

A novel HLA-B5 CREG gene, HLA-B SNA was cloned and the primary structure was determined. The sequence data showed that HLA-B SNA was identical to HLA-B51 except the alpha 1 domain in which one amino acid substitution at residue 74 and 5 amino acid substitutions associated with the Bw4/Bw6 epitopes were observed between these Ag. The comparison with other HLA-B locus genes suggested that HLA-B SNA evolved from HLA-B51 by gene exchange or recombination at the exon 2 between HLA-B51 and B8. A total of 10 of 14 HLA-B51-specific CTL clones showed significantly weak or no recognition of HLA-B SNA Ag. They also gave the same degree of a lysis of Hmy2CIR cells expressing the HLA-B35/51 chimeric Ag composed of the alpha 1 domain of HLA-B35 and other domains of HLA-B51 as that of Hmy2CIR cells expressing the HLA-B SNA Ag. These results demonstrated that amino acid substitutions within positions 77-83 associated with the HLA-Bw4/Bw6 epitopes have an influence on recognition of the HLA-B SNA antigen by HLA-B51-specific CTL.     

Species and recombination effects on DNA variability in the tomato genus

Population genetics theory predicts that strong selection for rare, beneficial mutations or against frequent, deleterious mutations reduces polymorphism at linked neutral (or weakly selected) sites. The reduction of genetic variation is expected to be more severe when recombination rates are lower. In outbreeding species, low recombination rates are usually confined to certain chromosomal regions, such as centromeres and telomeres. In contrast, in predominantly selfing species, the rarity of double heterozygotes leads to a reduced effective recombination rate in the whole genome. We investigated the effects of restricted recombination on DNA polymorphism in these two cases, analyzing five Lycopersicon species with contrasting mating systems: L. chilense, L. hirsutum, L. peruvianum, L. chmielewskii, and L. pimpinellifolium, of which only the first three species have self-incompatibility alleles. In each species, we determined DNA sequence variation of five single-copy genes located in chromosomal regions with either high or low recombination rate. We found that the mating system has a highly significant effect on the level of polymorphism, whereas recombination has only a weak influence. The effect of recombination on levels of polymorphism in Lycopersicon is much weaker than in other well-studied species, including Drosophila. To explain these observations, we discuss a number of hypotheses, invoking selection, recombination, and demographic factors associated with the mating system. We also provide evidence that L. peruvianum, showing a level of polymorphism (almost 3%) that is comparable to the level of divergence in the whole genus, is the ancestral species from which the other species of the genus Lycopersicon have originated relatively recently.     

Human MHC class III genes, Bf and C4. Polymorphism, complotypes and association with MHC class I genes in the Finnish population

Electrophoretically detected genetic polymorphism of human MHC class III genes, factor B (Bf) and complement C4A and C4B, was studied in the Finnish population. Bf alleles were determined in a panel of sera from 70 unrelated individuals. The common Bf alleles, Bf*S and Bf*F, had frequencies of 73% and 26%, respectively. Only in 1 individual was another allele, Bf*F1, detected. The frequencies of the C4A and C4B alleles were based on studies of 254 unrelated individuals. In this panel, five different alleles were detected at the C4A locus and four at the C4B locus. At both loci an allele without a gene product, i.e. a 'null' allele, was observed with high frequency, 11% for C4A 'null' and 17% for C4B 'null'. The association of complotypes to HLA haplotypes was analyzed in 70 chromosomes. The most common combination, defined by class I and class III alleles, was HLA-B7-S31 (13%), followed by HLA-B35-F20 (8.4%) and HLA-B8-S03 (7.1%). Some HLA-B specificities, for example B15, B27 and B40, were associated with a variety of complotypes. The importance of complotyping in HLA genetics is discussed.     

Detecting the signature of selection on immune genes in highly structured populations of wild sheep (Ovis dalli)

The confounding effects of population structure complicate efforts to identify regions of the genome under the influence of selection in natural populations. Here we test for evidence of selection in three genes involved in vertebrate immune function - the major histocompatibility complex (MHC), interferon gamma (IFNG) and natural resistance associated macrophage polymorphism (NRAMP) - in highly structured populations of wild thinhorn sheep (Ovis dalli). We examined patterns of variation at microsatellite loci linked to these gene regions and at the DNA sequence level. Simple Watterson's tests indicated balancing selection at all three gene regions. However, evidence for selection was confounded by population structure, as the Watterson's test statistics from linked markers were not outside of the range of values from unlinked and presumably neutral microsatellites. The translated coding sequences of thinhorn IFNG and NRAMP are fixed and identical to those of domestic sheep (Ovis aries). In contrast, the thinhorn MHC DRB locus shows significant evidence of overdominance through both an excess of nonsynonymous substitution and trans-species polymorphism. The failure to detect balancing selection at microsatellite loci linked to the MHC is likely the result of recombination between the markers and expressed gene regions.     

Molecular evolution of recombination hotspots and highly recombining pseudoautosomal regions in hominoids

We examined the effects of recombination on the molecular evolution of noncoding regions in pseudoautosomal regions (PARs) and recombination hotspots in hominoids. The PAR-linked regions analyzed had on average longer branch lengths than those of the recombination hotspots. Moreover, contrary to previous observations, we found no correlation between recombination rate and silent site divergence in our data set and little change in the GC content during recent hominoid evolution. This suggests that the current rate of recombination is not a good indicator of the past rates of recombination for these highly recombining regions. Furthermore, human recombination hotspots show increased AT to GC substitutions in the human lineage, while no such pattern is detected for PAR-linked regions. Together, these observations suggest that recombination hotspots in hominoids are transient in the evolutionary time-scale. Interestingly, the 16p13.3 recombination hotspot locus violates a local molecular clock, though the locus appears to be noncoding and should evolve neutrally. We hypothesize that sudden changes in recombination rate have caused the changes in substitution rate at this locus.     

Balancing selection and low recombination affect diversity near the self-incompatibility loci of the plant Arabidopsis lyrata

The self-incompatibility (S-) locus region of plants in the Brassica family is a small genome region. In Arabidopsis lyrata, the S-genes, SRK and SCR, encode the functional female and pollen recognition proteins, which must be coadapted to maintain correct associations between the two component genes, and thus self-incompatibility (SI). Recombinants would be self-compatible and thus probably disadvantageous in self-incompatible species. Therefore, tight linkage between the two genes in incompatibility systems is predicted to evolve to avoid producing such recombinant haplotypes. The evolution of low recombination in S-locus regions has not been rigorously tested. To test whether these regions' per-nucleotide recombination rates differ from those elsewhere in the genome, and to investigate whether the A. lyrata S-loci have the predicted effect on diversity in their immediate genome region, we studied diversity in genes that are linked to the S-loci but are not involved in incompatibility and are not under balancing selection. Compared with other A. lyrata loci, genes linked to the S-loci have extraordinarily high polymorphism. Our estimated recombination in this region, from fitting a model of the effects of S-allele polymorphism on linked neutral sites, supports the hypothesis of locally suppressed recombination around the S-locus.     

Global Analysis of ATM Polymorphism Reveals Significant Functional Constraint

ATM, the gene that is mutated in ataxia-telangiectasia, is associated with cerebellar degeneration, abnormal proliferation of small blood vessels, and cancer. These clinically important manifestations have stimulated interest in defining the sequence variation in the ATM gene. Therefore, we undertook a comprehensive survey of sequence variation in ATM in diverse human populations. The protein-encoding exons of the gene (9,168 bp) and the adjacent intron and untranslated sequences (14,661 bp) were analyzed in 93 individuals from seven major human populations. In addition, the coding sequence was analyzed in one chimpanzee, one gorilla, one orangutan, and one Old World monkey. In human ATM, 88 variant sites were discovered by denaturing high-performance liquid chromatography, which is 96%-100% sensitive for detection of DNA sequence variation. ATM was compared to 14 other autosomal genes for nucleotide diversity. The noncoding regions of ATM had diversity values comparable to other genes, but the coding regions had very low diversity, especially in the last 29% of the protein sequence. A test of the neutral evolution hypothesis, through use of the Hudson/Kreitman/Aguadé statistic, revealed that this region of the human ATM gene was significantly constrained relative to that of the orangutan, the Old World monkey, and the mouse, but not relative to that of the chimpanzee or the gorilla. ATM displayed extensive linkage disequilibrium, consistent with suppression of meiotic recombination at this locus. Seven haplotypes were defined. Two haplotypes accounted for 82% of all chromosomes analyzed in all major populations; two others carrying the same D126E missense polymorphism accounted for 33% of chromosomes in Africa but were never observed outside of Africa. The high frequency of this polymorphism may be due either to a population expansion within Africa or to selective pressure.     

The HKA test revisited: a maximum-likelihood-ratio test of the standard neutral model

We present a maximum-likelihood-ratio test of the standard neutral model, using multilocus data on polymorphism within species and divergence between species. The model is based on the Hudson-Kreitman-Aguade (HKA) test, but allows for an explicit test of selection at individual loci in a multilocus framework. We use coalescent simulations to show that the likelihood-ratio test statistic is conservative, particularly when the assumption of no recombination is violated. Application of the method to polymorphism data from 18 loci from a population of Arabidopsis lyrata provides significant evidence for a balanced polymorphism at a candidate locus thought to be linked to the centromere. The method is also applied to polymorphism data in maize, providing support for the hypothesis of directional selection on genes in the starch pathway.     

Selection at linked sites in the partial selfer Caenorhabditis elegans

Natural selection can produce a correlation between local recombination rates and levels of neutral DNA polymorphism as a consequence of genetic hitchhiking and background selection. Theory suggests that selection at linked sites should affect patterns of neutral variation in partially selfing populations more dramatically than in outcrossing populations. However, empirical investigations of selection at linked sites have focused primarily on outcrossing species. To assess the potential role of selection as a determinant of neutral polymorphism in the context of partial self-fertilization, we conducted a multivariate analysis of single-nucleotide polymorphism (SNP) density throughout the genome of the nematode Caenorhabditis elegans. We based the analysis on a published SNP data set and partitioned the genome into windows to calculate SNP densities, recombination rates, and gene densities across all six chromosomes. Our analyses identify a strong, positive correlation between recombination rate and neutral polymorphism (as estimated by noncoding SNP density) across the genome of C. elegans. Furthermore, we find that levels of neutral polymorphism are lower in gene-dense regions than in gene-poor regions in some analyses. Analyses incorporating local estimates of divergence between C. elegans and C. briggsae indicate that a mutational explanation alone is unlikely to explain the observed patterns. Consequently, we interpret these findings as evidence that natural selection shapes genome-wide patterns of neutral polymorphism in C. elegans. Our study provides the first demonstration of such an effect in a partially selfing animal. Explicit models of genetic hitchhiking and background selection can each adequately describe the relationship between recombination rate and SNP density, but only when they incorporate selfing rate. Clarification of the relative roles of genetic hitchhiking and background selection in C. elegans awaits the development of specific theoretical predictions that account for partial self-fertilization and biased sex ratios.     

Diversity and diversification of HLA-A,B,C alleles

The nucleotide sequences encoding 14 HLA-A,B,C and 5 ChLA-A,B,C molecules have been determined. Combining these sequences with published data has enabled the polymorphism in 40 HLA-A,B,C and 9 ChLA-A,B,C alleles to be analyzed. Diversity is generated through assortment of point mutations by recombinational mechanisms including gene and allelic conversions. The distribution and frequency of silent and replacement substitutions indicate that there has been positive selection for allelic diversity in the 5' part of the gene (exons 1 to 3) and for allelic homogenization and locus specificity in the 3' part of the gene (exons 4 to 8). These differences may correlate with the lengths of converted sequences in the two parts of the gene and frequency of the CpG dinucleotide. Locus-specific divergence of HLA-A,B, and C demonstrates that recombinational events involving alleles of a locus have been more important than conversion between loci. This contrasts with the predominance of gene conversion events in the evolution of mutants of the H-2Kb gene. However, a striking example of gene conversion involving HLA-B and C alleles of an oriental haplotype has been found. Comparison of human and chimpanzee alleles reveals extensive sharing of polymorphisms, confirming that diversification is a slow process, and that much of contemporary polymorphism originated in ancestral primate species before the emergence of Homo sapiens. There is less polymorphism at the HLA-A locus compared to HLA-B, with greater similarity also being seen between HLA-A and ChLA-A alleles than between HLA-B and ChLA-B alleles. Although greater diversity is seen in the 5' "variable" exons of HLA-B compared to HLA-A, there is increased heterogeneity in the 3' "conserved" exons of HLA-A compared to HLA-B.