PubMed Health⌕ Search

Biomedical subjects

M Nei

Publications and source records attributed to M Nei.

At least 19 recordsLinked to original sources
PubMed1992

Statistical properties of the ordinary least-squares, generalized least-squares, and minimum-evolution methods of phylogenetic inference.

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.

Biological Evolution↗
PubMed1991

Relative efficiencies of the maximum-parsimony and distance-matrix methods of phylogeny construction for restriction data.

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.

Computer Simulation↗
PubMed1990

Allelic genealogy under overdominant and frequency-dependent selection and polymorphism of major histocompatibility complex loci.

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.

Alleles↗
PubMed1990

A simple method for estimating average number of nucleotide substitutions within and between populations from restriction data.

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.

Analysis of Variance↗
PubMed1990

Limitations of the evolutionary parsimony method of phylogenetic analysis.

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.

Base Sequence↗
PubMed1990

Evolutionary relationships of class II major-histocompatibility-complex genes in mammals.

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.

Amino Acids↗
PubMed1990

Positive Darwinian selection promotes charge profile diversity in the antigen-binding cleft of class I major-histocompatibility-complex molecules.

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.

Alleles↗
PubMed1989

Nucleotide substitution at major histocompatibility complex class II loci: evidence for overdominant selection.

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.

Alleles↗
PubMed1989

Evolution of the major histocompatibility complex: independent origin of nonclassical class I genes in different groups of mammals.

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.

Animals↗
PubMed1989

[Study on molecular evolution for both variable segments of immunoglobulin heavy chain and T cell receptor].

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.

Animals↗
PubMed1979

Mathematical model for studying genetic variation in terms of restriction endonucleases.

A mathematical model for the evolutionary change of restriction sites in mitochondrial DNA is developed. Formulas based on this model are presented for estimating the number of nucleotide substitutions between two populations or species. To express the degree of polymorphism in a population at the nucleotide level, a measure called "nucleotide diversity" is proposed.

Base Sequence↗
PubMed1978

Nonrandom amino acid substitution and estimation of the number of nucleotide substitutions in evolution.

A method of estimating the number of nucleotide substitutions from amino acid sequence data is developed by using Dayhoff's mutation probability matrix. This method takes into account the effect of nonrandom amino acid substitutions and gives an estimate which is similar to the value obtained by Fitch's counting method, but larger than the estimate obtained under the assumption of random substitutions (Jukes and Cantor's formula). Computer simulations based on Dayhoff's mutation probability matrix have suggested that Jukes and Holmquist's method of estimating the number of nucleotide substitutions gives an overestimate when amino acid substitution is not random and the variance of the estimate is generally very large. It is also shown that when the number of nucleotide substitutions is small, this method tends to give an overestimate even when amino acid substitution is purely at random.

Amino Acid Sequence↗
PubMed1978

Goodman et al.'s method for augmenting the number of nucleotide substitutions.

Statistical properties of Goodman et al.'s (1974) method of compensating for undetected nucleotide substitutions in evolution are investigated by using computer simulation. It is found that the method tends to overcompensate when the stochastic error of the number of nucleotide substitutions is large. Furthermore, the estimate of the number of nucleotide substitutions obtained by this method has a large variance. However, in order to see whether this method gives overcompensation when applied together with the maximum parsimony method, a much larger scale of simulation seems to be necessary.

Biological Evolution↗