Prospects for Inferring Very Large Phylogenies by Using the Neighbor-joining Method

Overview

Journal Proc Natl Acad Sci U S A

Specialty Science

Date 2004 Jul 20

PMID 15258291

Citations 1613

Authors

Koichiro Tamura

Masatoshi Nei

Sudhir Kumar

Affiliations

Soon will be listed here.

Abstract

Current efforts to reconstruct the tree of life and histories of multigene families demand the inference of phylogenies consisting of thousands of gene sequences. However, for such large data sets even a moderate exploration of the tree space needed to identify the optimal tree is virtually impossible. For these cases the neighbor-joining (NJ) method is frequently used because of its demonstrated accuracy for smaller data sets and its computational speed. As data sets grow, however, the fraction of the tree space examined by the NJ algorithm becomes minuscule. Here, we report the results of our computer simulation for examining the accuracy of NJ trees for inferring very large phylogenies. First we present a likelihood method for the simultaneous estimation of all pairwise distances by using biologically realistic models of nucleotide substitution. Use of this method corrects up to 60% of NJ tree errors. Our simulation results show that the accuracy of NJ trees decline only by approximately 5% when the number of sequences used increases from 32 to 4,096 (128 times) even in the presence of extensive variation in the evolutionary rate among lineages or significant biases in the nucleotide composition and transition/transversion ratio. Our results encourage the use of complex models of nucleotide substitution for estimating evolutionary distances and hint at bright prospects for the application of the NJ and related methods in inferring large phylogenies.

Citing Articles

Staphylococcus borealis causing subclinical mastitis in dromedary camels.

Ranjan R, Kuniyal A, Ranjan A, Singh M, Jyotsana B, Sawal R Trop Anim Health Prod. 2025; 57(2):117.

PMID: 40085299 DOI: 10.1007/s11250-025-04372-0.

First report of the emerging pathogenic yeast in commercially relevant fish from southeastern Brazil.

Oliveira M, Lima V, da Costa G, Baptista B, Augusto J, Hauser-Davis R Front Fungal Biol. 2025; 6:1545572.

PMID: 40066366 PMC: 11891357. DOI: 10.3389/ffunb.2025.1545572.

Investigation of anti-microbial and cytotoxic potential of Streptomyces werraensis GRS9 derived from the sediments of river Ganga.

Malik D, Panwar H, Rani S, Gupta K Braz J Microbiol. 2025; .

PMID: 40064810 DOI: 10.1007/s42770-025-01642-9.

Assessing the potential for fall armyworm exchanges between the two American continents across the Mexico-Central America land bridge.

Nagoshi R, Malo E, Cruz-Esteban S, Rosas-Garcia N, Herrera-Mayorga V, Meagher R PLoS One. 2025; 20(3):e0308501.

PMID: 40029859 PMC: 11875382. DOI: 10.1371/journal.pone.0308501.

Oral Microbiota and Inflammatory Bowel Diseases: Detection of Emerging Fungal Pathogens and Herpesvirus.

Oliveira M, Campos L, Brito F, de Carvalho F, Silva-Junior G, da Costa G Biomedicines. 2025; 13(2).

PMID: 40002893 PMC: 11852465. DOI: 10.3390/biomedicines13020480.

References

Desper R, Gascuel O . Fast and accurate phylogeny reconstruction algorithms based on the minimum-evolution principle. J Comput Biol. 2002; 9(5):687-705. DOI: 10.1089/106652702761034136. View

Joost P, Methner A . Phylogenetic analysis of 277 human G-protein-coupled receptors as a tool for the prediction of orphan receptor ligands. Genome Biol. 2002; 3(11):RESEARCH0063. PMC: 133447. DOI: 10.1186/gb-2002-3-11-research0063. View

Irving J, Askew D, Whisstock J . Computational analysis of evolution and conservation in a protein superfamily. Methods. 2003; 32(2):73-92. DOI: 10.1016/s1046-2023(03)00200-7. View

Ruiz-Pesini E, Mishmar D, Brandon M, Procaccio V, Wallace D . Effects of purifying and adaptive selection on regional variation in human mtDNA. Science. 2004; 303(5655):223-6. DOI: 10.1126/science.1088434. View

Desper R, Gascuel O . Theoretical foundation of the balanced minimum evolution method of phylogenetic inference and its relationship to weighted least-squares tree fitting. Mol Biol Evol. 2003; 21(3):587-98. DOI: 10.1093/molbev/msh049. View

Kumar S, Tamura K, Nei M . MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004; 5(2):150-63. DOI: 10.1093/bib/5.2.150. View

Kimura M . A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. J Mol Evol. 1980; 16(2):111-20. DOI: 10.1007/BF01731581. View

Hasegawa M, Kishino H, Yano T . Dating of the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985; 22(2):160-74. DOI: 10.1007/BF02101694. View

Saitou N, Nei M . The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987; 4(4):406-25. DOI: 10.1093/oxfordjournals.molbev.a040454. View

10.

Jin L, Nei M . Limitations of the evolutionary parsimony method of phylogenetic analysis. Mol Biol Evol. 1990; 7(1):82-102. DOI: 10.1093/oxfordjournals.molbev.a040588. View

11.

Tamura K . Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G+C-content biases. Mol Biol Evol. 1992; 9(4):678-87. DOI: 10.1093/oxfordjournals.molbev.a040752. View

12.

Tamura K, Nei M . Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol. 1993; 10(3):512-26. DOI: 10.1093/oxfordjournals.molbev.a040023. View

13.

Tajima F . Unbiased estimation of evolutionary distance between nucleotide sequences. Mol Biol Evol. 1993; 10(3):677-88. DOI: 10.1093/oxfordjournals.molbev.a040031. View

14.

Gu X, Fu Y, Li W . Maximum likelihood estimation of the heterogeneity of substitution rate among nucleotide sites. Mol Biol Evol. 1995; 12(4):546-57. DOI: 10.1093/oxfordjournals.molbev.a040235. View

15.

Kumar S . A stepwise algorithm for finding minimum evolution trees. Mol Biol Evol. 1996; 13(4):584-93. DOI: 10.1093/oxfordjournals.molbev.a025618. View

16.

Rzhetsky A, Sitnikova T . When is it safe to use an oversimplified substitution model in tree-making?. Mol Biol Evol. 1996; 13(9):1255-65. DOI: 10.1093/oxfordjournals.molbev.a025691. View

17.

Nei M . Phylogenetic analysis in molecular evolutionary genetics. Annu Rev Genet. 1996; 30:371-403. DOI: 10.1146/annurev.genet.30.1.371. View

18.

Nei M, Kumar S, Takahashi K . The optimization principle in phylogenetic analysis tends to give incorrect topologies when the number of nucleotides or amino acids used is small. Proc Natl Acad Sci U S A. 1998; 95(21):12390-7. PMC: 22842. DOI: 10.1073/pnas.95.21.12390. View

19.

Takahashi K, Nei M . Efficiencies of fast algorithms of phylogenetic inference under the criteria of maximum parsimony, minimum evolution, and maximum likelihood when a large number of sequences are used. Mol Biol Evol. 2000; 17(8):1251-8. DOI: 10.1093/oxfordjournals.molbev.a026408. View

20.

Kumar S, Gadagkar S . Efficiency of the neighbor-joining method in reconstructing deep and shallow evolutionary relationships in large phylogenies. J Mol Evol. 2000; 51(6):544-53. DOI: 10.1007/s002390010118. View