A comparison of two methods for constructing evolutionary distances from a weighted contribution of transition and transversion differences

Since the initial work of Jukes and Cantor (1969), a number of procedures have been developed to estimate the expected number of nucleotide substitutions corresponding to a given observed level of nucleotide differentiation assuming particular evolutionary models. Unlike the proportion of different sites, the expected number of substitutions that would have occurred grows linearly with time and therefore has had great appeal as an evolutionary distance. Recently, however, a number of authors have tried to develop improved statistical approaches for generating and evaluating evolutionary distances (Schoniger and von Haeseler 1993; Goldstein and Polock 1994; Tajima and Takezaki 1994). These studies clearly show that the estimated number of nucleotide substitutions is generally not the best estimator for use in reconstruction of phylogenetic relationships. The reason for this is that there is often a large error associated with the estimation of this number. Therefore, even though its expectation is correct (i.e., on average the expected number of substitutions is proportional to time- -but see Tajima 1993), it is not expected to be as useful as estimators designed to have a lower variance.      

Molecular phylogeny of the Chinese ranids inferred from nuclear and mitochondrial DNA sequences

The phylogeny of representative species of Chinese ranids was reconstructed using two nuclear (tyrosinase and rhodopsin) and two mitochondrial (12S rRNA, 16S rRNA) DNA fragments. Maximum parsimony, Bayesian, and maximum likelihood analyses were employed. In comparison with the results from nuclear and mitochondrial data, we used nuclear gene data as our preferred phylogenetic hypothesis. We proposed two families (Ranidae, Dicroglossidae) for Chinese ranids, with the exception of genus Ingerana. Within Dicroglossidae, four tribes were supported including Dicroglossini, Paini, Limnonectini, and Occidozygini. A broader sampling strategy and evidence from additional molecular markers are required to decisively evaluate the evolutionary history of Chinese ranids.     

Molecular phylogeny of Moenkhausia (Characidae) inferred from mitochondrial and nuclear DNA evidence

Moenkhausia is one of the most speciose genera in Characidae, currently composed of 75 nominal species of small fishes distributed across South American hydrographic basins, primarily the Amazon and Guyanas. Despite the large number of described species, studies involving a substantial number of its species designed to better understand their relationships and putative monophyly are still lacking. In this study, we analysed a large number of species of Moenkhausia to test the monophyly of the genus based on the phylogenetic analysis of DNA sequences of two mitochondrial and three nuclear genes. The in‐group included 29 species of Moenkhausia, and the out‐group was composed of representatives of Characidae and other members of Characiformes. All species of Moenkhausia belong to the same clade (Clade C); however, they appear distributed in five monophyletic groups along with other different genera, which means that Moenkhausia is polyphyletic and indicates the necessity of an extensive revision of the group.     

Molecular phylogeny of Clupeiformes (Actinopterygii) inferred from nuclear and mitochondrial DNA sequences

The taxonomy of clupeiforms has been extensively studied, yet phylogenetic relationships among component taxa remain controversial or unresolved. Here we test current and new hypotheses of relationships among clupeiforms using mitochondrial rRNA genes (12S and 16S) and nuclear RAG1 and RAG2 sequences (total of 4749bp) for 37 clupeiform taxa representing all five extant families and all subfamilies of Clupeiformes, except Pristigasterinae, plus seven outgroups. Our results, based on maximum parsimony, maximum likelihood, and Bayesian analyses of these data, show that some traditional hypotheses are supported. These include the monophyly of the families Engraulidae, consisting of two monophyletic subfamilies, Engraulinae (Engraulis and Anchoa) and Coilinae (Coilia and Setipinna), and Pristigasteridae (here represented only by Ilisha and Pellona). The basal position of Denticeps among clupeiforms is consistent with the molecular data when base compositional biases are accounted for. However, the monophyly of Clupeidae was not supported. Some clupeids were more closely related to taxa assigned to Pristigasteridae and Chirocentridae (Chirocentrus). These results suggest that a major revision in the classification of clupeiform fishes may be necessary, but should await a more complete taxonomic sampling and additional data.     

Snake phylogeny: evidence from nuclear and mitochondrial genes

We constructed phylogenies of snakes from the c-mos and cytochrome b genes using conventional phylogenetic methods as well as the relatively new method of Bayesian inference. For all methods, there was excellent congruence between the c-mos and cytochrome b genes, implying a high level of support for the shared clades. Our results agree with previous studies in two important respects: first, that the scolecophidians and alethinophidians are monophyletic sister clades; and second, that the Colubroidea is a monophyletic group with the Acrochordidae as its sister clade. However, our results differ from previous studies in the finding that Loxocemus and Xenopeltis cluster with pythons. An additional noteworthy result from our data is that the genera Exiliboa and Ungaliophis, often placed with Tropidophis (and Trachyboa, not included in the present study) in the Tropidophiidae, are in reality boids.     

Molecular phylogeny of elasmobranchs inferred from mitochondrial and nuclear markers

The elasmobranchs (sharks, rays and skates) being the extant survivors of one of the earliest offshoots of the vertebrate evolutionary tree are good model organisms to study the primitive vertebrate conditions. They play a significant role in maintaining the ecological balance and have high economic value. Due to over-exploitation and illegal fishing worldwide, the elasmobranch stocks are being decimated at an alarming rate. Appropriate management measures are necessary for restoring depleted elasmobranch stocks. One approach for restoring stocks is implementation of conservation measures and these measures can be formulated effectively by knowing the evolutionary relationship among the elasmobranchs. In this study, a total of 30 species were chosen for molecular phylogeny studies using mitochondrial cytochrome c oxidase subunit I, 12S ribosomal RNA gene and nuclear Internal Transcribed Spacer 2. Among different genes, the combined dataset of COI and 12S rRNA resulted in a well resolved tree topology with significant bootstrap/posterior probabilities values. The results supported the reciprocal monophyly of sharks and batoids. Within Galeomorphii, Heterodontiformes (bullhead sharks) formed as a sister group to Lamniformes (mackerel sharks): Orectolobiformes (carpet sharks) and to Carcharhiniformes (ground sharks). Within batoids, the Myliobatiformes formed a monophyly group while Pristiformes (sawfishes) and Rhinobatiformes (guitar fishes) formed a sister group to all other batoids.     

Transition and transversion rate in the evolution of animal mitochondrial DNA

We present a further application of the stochastic model previously described (Lanave et al., 1984, 1985) for measuring the nucleotide substitution rate in the mammalian evolution of the mitochondrial DNA (mtDNA). The applicability of this method depends on the validity of "stationarity conditions" (equal nucleotide frequencies at first, second and third silent codon positions in homologous protein coding genes). In the comparison of homologous sequences satisfying the stationarity condition at the silent sites, only the four codon families (quartets) for which both transitions and transversions are silent at the third position are considered here. This has allowed us to estimate the transition and transversion rates for any pair of species. We have analyzed the third silent codon position of the triplet rat-mouse-cow, of a series of slightly divergent primates and of two Drosophila species. In terms of two external dating input we have then determined the phylogenetic trees for rat, mouse, and cow as well as for a number of primates including man. The phylogenetic tree that we have derived for the triplet rat, mouse and cow agrees with that we had previously determined by analyzing the first, second and third silent codon positions (in both duets and quartets) of mt genes (Lanave et al., 1985). For primates our method leads to the following branching order from the oldest to the most recent: Gibbon, Orangutan, Gorilla, Chimpanzee and Man. In absolute time, fixing the distance Chimpanzee-Man as 5 million years (Myr) we estimate the dating of the divergence nodes as: Gorilla 7 Myr; Orangutan 16 Myr; Gibbon 20 Myr. In all cases analyzed, the transition rate has been found to be substantially higher than the transversion rate. Moreover we have found that the transition/transversion ratio is different in the various lineages. We suggest that this fact is probably related to the nucleotide frequencies at the third silent codon position.     

Molecular phylogeny of a red-snow-crab species complex using mitochondrial and nuclear DNA markers

Nucleotide sequence variation of mitochondrial DNA COI and nuclear rRNA gene regions was used to reconstruct phylogenetic relationships for the red-snow-crab species complex, including the red snow crab, Chionoecetes japonicus, its nominal subspecies, C. japonicus pacificus, and the triangle tanner crab, C. angulatus. The topologies of the Bayesian and neighbor-joining (NJ) trees of the COI and of NJ trees of rRNA sequences placed C. japonicus and C. angulatus in a single clade. The net sequence divergence between these taxa was d(net) = 0.000 in COI, and strongly suggests that these taxa represent a single species. In contrast, haplotypes in C. j. pacificus clustered separately from the C. japonicus - C. angulatus clade. Net sequence divergence from C. japonicus - C. angulatus to C. j. pacificus was d(net) = 0.026 in COI, indicating that C. j. pacificus should be elevated to a separate species, C. pacificus. A 165 bp insert appeared in the rRNA gene of C. j. pacificus, but was absent in the remaining species of Chionoecetes. This autapomorphic condition in C. j. pacificus adds support for an independent evolution of this taxon. Evolutionary divergences between these taxa may reflect contrasting evolutionary process influenced by ocean bathymetry.     

DNA barcoding of genus Metapenaeopsis (Decapoda: Penaeidae) and molecular phylogeny inferred from mitochondrial and nuclear DNA sequences

Metapenaeopsis Bouvier, 1905 is the most diverse genus within Penaeidae. Metapenaeopsis shrimps exhibited subtle morphological differences, which make identification a difficult task based on taxonomic keys alone. In this study, we carried out DNA barcoding and phylogenetic analyses to examine taxonomy and phylogeny of genus Metapenaeopsis based on mitochondrial (COI) and nuclear (PEPCK and NaK) genes. Despite limited performance of DNA barcoding in delineating Metapenaeopsis shrimps, it questioned the taxonomic status of the two subspecies, Metapenaeopsis mogiensis intermedia and Metapenaeopsis mogiensis mogiensis, as well as three separate species: Metapenaeopsis provocatoria longirostris, Metapenaeopsis quinquedentata and Metapenaeopsis velutina. The major pattern of relationships between all studied taxa of Metapenaeopsis was similar across all analytical methods in which species with one-valved petasma were genetically distinct from those with two-valved petasma. As expected from morphology, the remaining species with stridulating organ constituted a strongly supported clade. In contrast, a paraphyletic clade was resolved for species without stridulating organ which contradicts Crosnier's morphological classification scheme for Metapenaeopsis. Overall, the present molecular data indicated that the shape of petasma and stridulating organ were both phylogenetically significant morphological characters for this genus, adding further evidence for the Crosnier's proposal.     

Molecular phylogeny of New World Myotis (Chiroptera, Vespertilionidae) inferred from mitochondrial and nuclear DNA genes

Recent studies have shown that species in the genus Myotis have evolved a number of convergent morphological traits, many of which are more related to their mode of food procurement than to their phylogeny. Surprisingly, the biogeographic origins of these species are a much better predictor of phylogenetic relationships, than their morphology. In particular, a monophyletic clade that includes all New World species was apparent, but only a third of the 38 species have been analysed. In order to better understand the evolution of this clade, we present phylogenetic reconstructions of 17 Nearctic and 13 Neotropical species of Myotis compared to a number of Old World congeners. These reconstructions are based on mitochondrial cytochrome b (1140 bp), and nuclear Rag 2 genes (1148 bp). Monophyly of the New World clade is strongly supported in all analyses. Two Palaearctic sister species, one from the west (M. brandtii) and one from the east (M. gracilis), are embedded within the New World clade, suggesting that they either moved across the Bering Strait, or that they descended from the same ancestor that reached the New World. An emerging feature of these phylogenetic reconstructions is that limited faunal exchanges have occurred, including between the North and South American continents, further emphasizing the importance of biogeography in the radiation of Myotis. A fossil-calibrated, relaxed molecular-clock model was used to estimate the divergence time of New World lineages to 12.2+/-2.0 MYA. Early diversification of New World Myotis coincides with the sharp global cooling of the Middle Miocene. Radiation of the temperate-adapted Myotis may have been triggered by these climatic changes. The relative paucity of species currently found in South America might result from a combination of factors including the early presence of competitors better adapted to tropical habitats.     

Towards the phylogeny of the Curculionoidea (Coleoptera): Reconstructions from mitochondrial and nuclear ribosomal DNA sequences

With about 60,000 described species, Curculionoidea represent the most species-rich superfamily in the animal kingdom. The immense diversity apparently creates difficulties in the reconstruction of the phylogenetic relationships. Independent morphological studies have led to very different classifications. This study is based on molecular data from two independent molecular sources, the 16S and 18S rDNA. Sensitivity analyses were conducted for the sequence alignment (gap costs were varied) as well as the phylogenetic reconstruction algorithms and some of their parameters. The higher-level relationships reconstructed within Curculionoidea are sensitive to alignment and reconstruction method. Nemonychidae or Oxycorynidae+Belidae were found to be sister to all remaining Curculionoidea in many analyses. The 16S rDNA sequence data (obtained from 157 species) corroborate many tribes and genera as monophyletic. It is observed that the phylogenetic reconstruction of genera with specific genetic features such as polyploidy and parthenogenetic reproduction is difficult in weevils. The curculionid subfamily Lixinae appears monophyletic. A new monophylum consisting of Entiminae, Hyperinae, Cyclominae, Myllorhinus plus possibly the Cossoninae is distinguished and we call it Entiminae s.l. For most other subfamilies and families homoplasy concealed the phylogenetic signal (due to saturation of the 16S sequences), or the species sampling was insufficient, although our sampling scheme was rather broad. We observed that although data from one source can easily be misleading (16S) or hardly informative (18S), the combination of the two independent data sets can result in useful information for such a speciose group of organisms. Our study represents the most thorough analysis of molecular sequence data of the Curculionoidea to date and although the phylogenetic results appear less stable than expected, they reflect the information content of these sequence data realistically and thus contribute to the total knowledge about the phylogeny of the Curculionoidea.     

Introgressive hybridization and the evolutionary history of the herring gull complex revealed by mitochondrial and nuclear DNA

Background  

Based on extensive mitochondrial DNA (mtDNA) sequence data, we previously showed that the model of speciation among species of herring gull (Larus argentatus) complex was not that of a ring species, but most likely due more complex speciation scenario's. We also found that two species, herring gull and glaucous gull (L. hyperboreus) displayed an unexpected biphyletic distribution of their mtDNA haplotypes. It was evident that mtDNA sequence data alone were far from sufficient to obtain a more accurate and detailed insight into the demographic processes that underlie speciation of this complex, and that extensive autosomal genetic analysis was warranted.     

Re-examination of the phylogeny of Rhacophoridae (Anura) based on mitochondrial and nuclear DNA

The phylogenetic relationships among rhacophorid frogs are under dispute. We use partial sequences of three mitochondrial (12S rRNA, 16S rRNA, and cytochrome b) and three nuclear protein-coding (Rag-1, rhodopsin exon 1, and tyrosinase exon 1) genes from 57 ingroup taxa and eight outgroup taxa to propose a hypothesis for phylogenetic relationships within Rhacophoridae. Our results support recognition of the genus Feihyla, and Chiromantis is the sister taxon to the clade formed by Feihyla, Polypedates and Rhacophorus. We place Aquixalus odontotarsus within Kurixalus, and the remaining species of Aquixalus and Philautus jinxiuensis into the genus Gracixalus. We give Philautus (Kirtixalus) the rank of genus and place Philautus menglaensis within it. The division of species groups among Chinese Rhacophorus needs revision, and a cryptic species is revealed within Rhacophorus nigropunctatus. Rhacophorus pingbianensis is considered a synonym of Rhacophorus omeimontis. The validity of Rhacophorus hui is confirmed by present molecular evidence.     

Global phylogeny of Hadrosauridae (Dinosauria: Ornithopoda) using parsimony and Bayesian methods

The late Cretaceous hadrosaurids were the most specialized and diverse clade of ornithopod dinosaurs. Parsimony and Bayesian methods were implemented to elucidate the phylogenetic relationships of all hadrosaurid species. Traditional and geometric morphometrics were applied to discover patterns of variation containing phylogenetic information. In total, 286 phylogenetically informative characters (196 cranial and 90 postcranial) were defined and documented: the most extensive character data set ever constructed for hadrosaurid dinosaurs. Of these, 136 characters were used for the first time in phylogenetic analysis of these ornithopods, and 93 were modified from those of other authors. Parsimony and the Bayesian analysis (using the Mk model without the gamma parameter) confirmed the split of hadrosaurids into Saurolophinae and Lambeosaurinae. Saurolophines included a major clade composed of the Prosaurolophus–Saurolophus and the Kritosaurus–Gryposaurus–Secernosaurus subclades. Edmontosaurus and Shantungosaurus were recovered outside the major clade of saurolophines. The Brachylophosaurus clade was recovered as the most basal clade of saurolophines in the parsimony analysis, whereas following the Bayesian analysis it was recovered as the sister clade to the Kritosaurus–Gryposaurus–Secernosaurus clade. These two analyses resulted in a Lambeosaurinae composed of a succession of Eurasian sister taxa to two major clades: the Parasaurolophus clade and the Hypacrosaurs altispinus–Corythosaurus clade. In contrast, the Bayesian analysis using the Mk model with the gamma parameter included, resulted in an unbalanced hadrosauroid tree, with a paraphyletic Saurolophinae, and with the Prosaurolophus clade, Edmontosaurus, and Shantungosaurus as successively closer sister taxa to Lambeosaurinae. Based on the strict reduced consensus tree derived from the parsimony analysis, Hadrosauridae was redefined as the clade stemming from the most recent common ancestor of Hadrosaurus foulkii and Parasaurolophus walkeri. © 2010 The Linnean Society of London, Zoological Journal of the Linnean Society, 2010, 159 , 435–502.     

Molecular phylogeny of Oryzeae (Poaceae) based on DNA sequences from chloroplast, mitochondrial, and nuclear genomes

The phylogeny and evolutionary history of the rice tribe (Oryzeae) were explored using sequences of five DNA fragments (matK, trnL, nad1, Adh2, and GPA1) from chloroplast, mitochondrial, and nuclear genomes. Results indicate that (1) Oryzeae is monophyletic and falls into two main clades corresponding to the traditionally recognized subtribes; (2) previous recognition of three monotypic genera (Hydrochloa, Porteresia, and Prosphytochloa) is not justified; and (3) close affinities of the monoecious genera are not supported, suggesting the possibility of multiple origins of unisexual florets. Based on the magnitude of matK and GPA1 sequence divergence, we suggest that Oryza and Leersia branched off from the remaining genera of Oryzeae ～20 million years ago (mya), and separated from each other ～14 mya. A divergence time of ～9 mya is obtained for the most basal split within Oryza. These estimates suggest that Oryzeae diverged during the Miocene, and thus imply that long-distance dispersal appears to be one of the important factors in the diversification of the tribe.     

Molecular phylogeny and evolution of primate mitochondrial DNA

We determined nucleotide sequences of homologous 0.9-kb fragments of mitochondrial DNAs (mtDNAs) derived from four species of old-world monkeys, one species of new-world monkeys, and two species of prosimians. With these nucleotide sequences and homologous sequences for five species of hominoids, we constructed a phylogenetic tree for the four groups of primates. The phylogeny obtained is generally consistent with evolutionary trees constructed in previous studies. Our results also suggest that the rate of nucleotide substitution for mtDNAs in hominines (human, chimpanzee, and gorilla) may have slowed down compared with that for old-world monkeys. This evolutionary feature of mitochondrial genes is similar to one found in nuclear genes.      

A comparison of methods for estimating the transition:transversion ratio from DNA sequences

Estimation of the ratio of the rates of transitions to transversions (TI:TV ratio) for a collection of aligned nucleotide sequences is important because it provides insight into the process of molecular evolution and because such estimates may be used to further model the evolutionary process for the sequences under consideration. In this paper, we compare several methods for estimating the TI:TV ratio, including the pairwise method [TREE 11 (1996) 158], a modification of the pairwise method due to Ina [J. Mol. Evol. 46 (1998) 521], a method based on parsimony (TREE 11 (1996) 158), a method due to Purvis and Bromham [J. Mol. Evol. 44 (1997) 112] that uses phylogenetically independent pairs of sequences, the maximum likelihood method, and a Bayesian method [Bioinformatics 17 (2001) 754]. We examine the performance of each estimator under several conditions using both simulated and real data.     

Correlation between strand asymmetry and phylogeny in mitochondrial DNA

An evolutionary distance is introduced in order to propose an efficient and feasible procedure for phylogeny studies. Our analysis are based on the strand asymmetry property of mitochondrial DNA, but can be applied to other genomes. Comparison of our results with those reported in conventional phylogenetic trees, gives confidence about our approximation. Our findings support the hypotheses about the origin of the skew and its dependence upon evolutionary pressures, and improves previous efforts on using the strand asymmetry property of genomes for phylogeny inference. For the evolutionary distance introduced here, we observe that the more adequate technique for tree reconstructions correspond to an average link method which employs a sequential clustering algorithm.     

Analysis of the mitochondrial DNA diversity in Yukaghirs in the evolutionary context

Based on the mtDNA first hypervariable segment sequence variation data, statistical analysis of the diversity in Yukaghirs in comparison with the other indigenous populations of Siberia, was carried out. The level of the Yukaghir mtDNA gene diversity (GD) constituted 0.920, which was only slightly different from the corresponding estimate for the other Siberian populations. Integral estimates of the genetic structure of Siberian populations (k, S, θ _S , and π) are presented. Phylogenetic analysis, performed using the neighbor-joining method, showed that the Siberian populations clustered irrespectively to their language affiliation. Negative F _s values found in Yukaghirs pointed to the possible influence of adaptive selection.