Maintenance of MHC polymorphism.
Explore the source record for details and available documents.
Biomedical subjects
Publications and source records attributed to M Nei.
Explore the source record for details and available documents.
Statistical properties of the ordinary least-squares (OLS), generalized least-squares (GLS), and minimum-evolution (ME) methods of phylogenetic inference were studied by considering the case of four DNA sequences. Analytical study has shown that all three methods are statistically consistent in the sense that as the number of nucleotides examined (m) increases they tend to choose the true tree as long as the evolutionary distances used are unbiased. When evolutionary distances (dij's) are large and sequences under study are not very long, however, the OLS criterion is often biased and may choose an incorrect tree more often than expected under random choice. It is also shown that the variance-covariance matrix of dij's becomes singular as dij's approach zero and thus the GLS may not be applicable when dij's are small. The ME method suffers from neither of these problems, and the ME criterion is statistically unbiased. Computer simulation has shown that the ME method is more efficient in obtaining the true tree than the OLS and GLS methods and that the OLS is more efficient than the GLS when dij's are small, but otherwise the GLS is more efficient.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The relative efficiencies of the maximum-parsimony (MP), UPGMA, and neighbor-joining (NJ) methods in obtaining the correct tree (topology) for restriction-site and restriction-fragment data were studied by computer simulation. In this simulation, six DNA sequences of 16,000 nucleotides were assumed to evolve following a given model tree. The recognition sequences of 20 different six-base restriction enzymes were used to identify the restriction sites of the DNA sequences generated. The restriction-site data and restriction-fragment data thus obtained were used to reconstruct a phylogenetic tree, and the tree obtained was compared with the model tree. This process was repeated 300 times. The results obtained indicate that when the rate of nucleotide substitution is constant the probability of obtaining the correct tree (Pc) is generally higher in the NJ method than in the MP method. However, if we use the average topological deviation from the model tree (dT) as the criterion of comparison, the NJ and MP methods are nearly equally efficient. When the rate of nucleotide substitution varies with evolutionary lineage, the NJ method is better than the MP method, whether Pc or dT is used as the criterion of comparison. With 500 nucleotides and when the number of nucleotide substitutions per site was very small, restriction-site data were, contrary to our expectation, more useful than sequence data. Restriction-fragment data were less useful than restriction-site data, except when the sequence divergence was very small. UPGMA seems to be useful only when the rate of nucleotide substitution is constant and sequence divergence is high.
To explain the long-term persistence of polymorphic alleles (trans-specific polymorphism) at the major histocompatibility complex (MHC) loci in rodents and primates, a computer simulation study was conducted about the coalescence time of different alleles sampled under various forms of selection. At the same time, average heterozygosity, the number of alleles in a sample, and the rate of codon substitution were examined to explain the mechanism of maintenance of polymorphism at the MHC loci. The results obtained are as follows. (1) The coalescence time for neutral alleles is too short to explain the trans-specific polymorphism at the MHC loci. (2) Under overdominant selection, the coalescence time can be tens of millions of years, depending on the parameter values used. The average heterozygosity and the number of alleles observed are also high enough to explain MHC polymorphism. (3) The pathogen adaptation model proposed by Snell is incapable of explaining MHC polymorphism, since the coalescence time for this model is too short and the expected heterozygosity and the expected number of alleles are too small. (4) From the mathematical point of view, the minority advantage model of frequency-dependent selection is capable of explaining a high degree of polymorphism and trans-specific polymorphism. (5) The molecular mimicry hypothesis also gives a sufficiently long coalescence time when the mutation rate is low in the host but very high in the parasite. However, the expected heterozygosity and the expected number of alleles tend to be too small. (6) Consideration of the molecular mechanism of the function of MHC molecules and other biological observations suggest that the most important factor for the maintenance of MHC polymorphism is overdominant selection. However, some experiments are necessary to distinguish between the overdominance and frequency-dependent selection hypotheses.
A simple method is proposed for estimating the average number of nucleotide substitutions per site within and between populations for the case where a large number of individuals are examined for many restriction enzymes. This method gives essentially the same results as those obtained by Nei and Li's method but saves a large amount of computer time. The variances of the quantities estimated can be obtained by the jackknife method, and these variances are very similar to those obtained by Nei and Jin's more sophisticated method. A similar method can also be applied to DNA sequence data.
Lake's evolutionary parsimony (EP) method of constructing a phylogenetic tree is primarily applied to four DNA sequences. In this method, three quantities--X, Y, and Z--that correspond to three possible unrooted trees are computed, and an invariance property of these quantities is used for choosing the best tree. However, Lake's method depends on a number of unrealistic assumptions. We therefore examined the theoretical basis of his method and reached the following conclusions: (1) When the rates of two transversional changes from a nucleotide are unequal, his invariance property breaks down. (2) Even if the rates of two transversional changes are equal, the invariance property requires some additional conditions. (3) When Kimura's two-parameter model of nucleotide substitution applies and the rate of nucleotide substitution varies greatly with branch, the EP method is generally better than the standard maximum-parsimony (MP) method in recovering the correct tree but is inferior to the neighbor-joining (NJ) and a few other distance matrix methods. (4) When the rate of nucleotide substitution is the same or nearly the same for all branches, the EP method is inferior to the MP method even if the proportion of transitional changes is high. (5) When Lake's assumptions fail, his chi2 test may identify an erroneous tree as the correct tree. This happens because the test is not for comparing different trees. (6) As long as a proper distance measure is used, the NJ method is better than the EP and MP methods whether there is a transition/transversion bias or whether there is variation in substitution rate among different nucleotide sites.
The major histocompatibility complex (MHC) class II molecule consists of noncovalently associated alpha and beta chains. In mammals studied so far, the class II MHC can be divided into a number of regions, each containing one or more alpha-chain genes (A genes) and beta-chain genes (B genes), and it has been known for some time that orthologous relationships exist between genes in corresponding regions from different mammalian species. A phylogenetic analysis of DNA sequences of class II A and B genes confirmed these relationships; but no such orthologous relationship was observed between the B genes of mammals and those of birds. Thus, the class II regions have diverged since the separation of birds and mammals (approximately 300 Mya) but before the radiation of the placental mammalian orders (60-80 Mya). Comparison of the phylogenetic trees for A and B genes revealed an unexpected characteristic of DP-region genes: DPB genes are most closely related to DQB genes, whereas DPA chain genes are most closely related to DRA-chain genes. Thus, the DP region seems to have originated through a recombinational event which brought together a DQB gene and a DRA gene (perhaps approximately 120 Mya). The 5' untranslated region of all class II genes includes sequences which are believed to be important in regulating class II gene expression but which are not conserved in known pseudogenes. These sequences are conserved to an extraordinary degree in the human DQB1 gene and its mouse homologue A beta 1, suggesting that regulation of expression of this locus may play a key role in expression of the entire class II MHC.
Certain major-histocompatibility-complex (MHC) loci are highly polymorphic, and the mechanism of maintenance of this polymorphism remains controversial. Recent studies of the pattern of nucleotide substitution at MHC loci have produced strong evidence that this polymorphism is maintained mainly by positive Darwinian selection that operates on the antigen recognition site (ARS) of the MHC molecule. The ARS of the class I MHC consists of three subregions: (1) the binding cleft, (2) T-cell-receptor-directed residues, and (3) outward-directed residues. Here we report that the rate of nonsynonymous nucleotide substitution is much higher in the binding cleft than in the other ARS subregions. Furthermore, nonsynonymous nucleotide substitutions that result in a change of residue side-chain charge occur significantly more frequently than expected by chance. We conclude that the main target of positive selection on the class I MHC molecules is the binding cleft of the ARS and that this selection acts primarily to promote diversity among alleles with respect to the pattern of residue side-chain charges (charge profile) in the binding cleft. These results provide additional support for the hypothesis that MHC polymorphism is maintained by overdominant selection relating to antigen-binding capacity and thus to disease resistance.
To study the mechanism of maintenance of polymorphism at major histocompatibility complex (MHC) loci, synonymous and nonsynonymous (amino acid-altering) nucleotide substitutions in the putative antigen-recognition site (included in the first domain of the MHC molecule) and other regions of human and mouse class II genes were examined. In the putative antigen-recognition site, the rate of nonsynonymous substitution was found to exceed that of synonymous substitution, whereas in the second domain the former was significantly lower than the latter. In light of a previous theoretical study and parallel findings in class I MHC loci, we conclude that the unusually high degree of polymorphism at class II MHC loci is caused mainly by overdominant selection (heterozygote advantage) operating in the antigen-recognition site.
The major histocompatibility complex (MHC) in humans and chimpanzees includes three classical class I loci, A, B and C, which encode glycoproteins expressed on the surface of all nucleated cells. There are also several nonclassical class I loci including E, which have more limited expression. By analyzing published sequences, we have shown that in exons 4 and 5, A locus alleles from both humans and chimpanzees are much more similar to E than to B or C alleles, whereas in exons 2 and 3 alleles from all three classical class I loci are much more similar to each other than any one is to E. We propose that some 20 million years ago, interlocus recombination led to the formation of a hybrid gene in which exons 2 and 3 were derived from the original A locus and exons 4 and 5 were derived from the E locus. The fact that such an ancient event can still be detected suggests that interlocus recombination is rare in the MHC and does not significantly contribute to MHC polymorphism, which is known to be extremely high. The present finding, however, supports Gilbert's idea that exons in a gene may occasionally be replaced by those from another gene in the evolutionary process.
Statistical methods for computing the variances of nucleotide diversity within populations and of nucleotide divergence between populations are developed. Both variances are computed by finding the phylogenetic relationships of the DNA sequences studied through the unweighted pair-group method or some other tree-making method. The methods developed are applicable to both DNA sequence and restriction-site map data.
To examine Gojobori and Nei's hypothesis that the immunoglobulin heavy-chain variable-region (VH) genes in mammals are subject to diversity-enhancing selection, we studied the rates of synonymous and nonsynonymous nucleotide substitution in the complementarity-determining regions (CDRs) and in the framework regions (FRs) of mouse and human VH genes. The results obtained indicate that the non-synonymous rate is higher than the synonymous rate in CDRs, whereas the reverse is true in FRs. This observation supports Gojobori and Nei's hypothesis and suggests that diversity-enhancing selection (similar to overdominant selection) operates mainly in CDRs and that this is one of the evolutionary factors that increase antibody diversity.
The class I major histocompatibility complex genes are composed of classical and nonclassical genes, the latter being largely nonfunctional. To understand the evolutionary relationships of the two groups of class I genes, a phylogenetic analysis of DNA sequences was conducted using 45 genes from six mammalian and one avian species. The results indicate that nonclassical genes in one species are more closely related to classical genes from the same species than to nonclassical genes from a species belonging to a different order or family. This indicates that the differentiation of classical and nonclassical genes occurs rather rapidly in the genome. Classical genes are apparently duplicated with a high frequency in the evolutionary process, and many of the duplicated genes seem to degenerate into nonclassical genes as a result of deleterious mutation. The nonclassical Qa genes in the mouse have sequences homologous to regulatory sequences involved in the universal expression of classical class I genes, but they have accumulated numerous nucleotide substitutions in these sequences. The pattern of nucleotide substitution in nonclassical genes is different from that in classical genes. In nonclassical genes, the rate of nonsynonymous substitution is higher in the antigen recognition site than in other gene regions, as is true of classical genes. However, unlike the case of classical genes, the nonsynonymous rate does not always exceed the synonymous rate in the antigen recognition site. Nonclassical proteins further differ from classical proteins in having amino acid replacements in conserved antigen recognition site positions. These observations are consistent with the hypothesis that nonclassical genes have originated from classical genes but have lost classical class I function because of deleterious mutation.
To study the evolutionary relationships of the three major groups of humans, Europeans, Asians and Africans, the genetic distances between them were computed by using 4 different sets of genetic loci (84 protein loci, 33 blood group loci, 8 HLA and immunoglobulin loci, and 61 DNA markers). The results obtained indicate that the overall genetic distance between Europeans and Asians is significantly lower than that between Europeans and Africans of that between Asians and Africans and support the hypothesis of an African origin of modern humans. This seems to be the first study to establish the evolutionary relationships of the three major groups of humans at a statistically significant level.
In order to explain difference and similarity in producing antibody diversity between immunoglobulin (Ig) and T cell receptor (TCR), authors compared both codon substitution and concerted evolution rate between the variable segment of Ig heavy (Ig VH) and that of TCR (TCR V). The protein sequences of TCR V alpha (including 8 gene segments from mouse and 3 from human), TCR V beta (including 11 from mouse and one from human) and T cell V gamma (including 2 from mouse and 4 from human) were compiled, as well as the protein sequences of Ig VH (3 from human, 11 from mouse, 3 from caiman and one from shark) were collected. It is shown that: (1) the nucleotide substitution of TCR V segment is 2.4 times as large as that of Ig VH in coding region; (2) as for concerted evolution, gene duplicate rates in TCR V and Ig VH are 1.7 X 10(-8) and 1.6 X 10(-8)/gene/year, respectively. The number of TCR V(V alpha equals to 100 and V beta equals to 30) is less than the one of Ig VH (VH equals to 300), for TCR V is subject to negative selection of major histocompatibility complex according to the neutral theory. We discussed that is somatic mutation or DNA rearrangement the main force in producing antibody diversity and are there pseudogenes in TCR V or not.
The extent and pattern of protein and DNA polymorphisms are discussed with emphasis on the mechanism of maintenance of the polymorphisms. Statistical studies suggest that a large proportion of genetic variability at the molecular level is maintained by a mutation-drift balance. At some loci, such as those for histocompatibility in mammals, however, a form of overdominant selection seems to be involved. In the presence of overdominant selection, polymorphic alleles may be maintained for tens of millions of years, so that the number of nucleotide differences between alleles is often very large, as in the case of self-incompatibility alleles in plants. There are also an increasing number of examples in which an adaptive change of a morphological or physiological character is caused by a single nucleotide substitution. Nevertheless, these mutations seem to be a small proportion of the total nucleotide changes that contribute to genetic variability and evolution. Although there are many examples of frequency-dependent selection, this form of selection is apparently unimportant for the maintenance of genetic variability except in some special cases. Observations on the evolutionary change of DNA suggest that the driving force of evolution is mutation rather than selection.