Applications of codon and rate variation models in molecular phylogeny

The current article illustrates the practical advantages of some new models and statistical algorithms for codon substitution and spatial rate variation in molecular phylogeny. Our companion paper in this issue discusses at length the mathematical properties of these models for nucleotide and codon substitution, for site-to-site and branch-to-branch heterogeneity in rates of evolution, and for spatial correlation in the assignment of rates. In this study we summarize the theoretical background and apply the models and algorithms to data on beta-globin, the complete HIV genome, and the mitochondrial genome. Our complex but realistic models enhance biological interpretation of sequence data and show substantial improvements in model fit over existing models. All the new statistical algorithms applied are incorporated in our phylogeny software LINNAEUS, which is tuned for performance and modeling flexibility.     

Long-branch attraction phenomenon and the impact of among-site rate variation on rodent phylogeny

The phylogenetic relationships among major lineages of rodents is one of the issues most debated by both paleontologists and molecular biologists. In the present study, we have analyzed all complete mammalian mitochondrial genomes available in the databases, including five rodent species (rat, mouse, dormouse, squirrel and guinea-pig). Phylogenetic analyses were performed on H-strand amino acid sequences by means of maximum-likelihood and on H-strand protein-coding and ribosomal genes by means of distance methods. Also, log-likelihood ratio tests were applied to different tree topologies under the assumption of rodent monophyly, paraphyly or polyphyly. The analyses significantly rejected rodent monophyly and showed the existence of two differentiated clades, one containing non-murids (dormouse, squirrel and guinea-pig) and the other containing murids (rat and mouse). Long-branch attraction between murids and the outgroups could not be responsible for the existence of two different rodent clades, as no significant differences in evolutionary rate have been observed, except in the case of the squirrel, which shows a lower rate. The impact of among-site rate variation models on the phylogeny of rodents has been evaluated using the gamma distribution model. Results have shown that relationships among rodents remained unchanged, and the general topology of the tree was not affected, even though some branches were not properly resolved, most likely due to a lack of fit between estimated and real rate heterogeneity parameters.     

Phylogenetic estimates of diversification rate are affected by molecular rate variation

Molecular phylogenies are increasingly being used to investigate the patterns and mechanisms of macroevolution. In particular, node heights in a phylogeny can be used to detect changes in rates of diversification over time. Such analyses rest on the assumption that node heights in a phylogeny represent the timing of diversification events, which in turn rests on the assumption that evolutionary time can be accurately predicted from DNA sequence divergence. But there are many influences on the rate of molecular evolution, which might also influence node heights in molecular phylogenies, and thus affect estimates of diversification rate. In particular, a growing number of studies have revealed an association between the net diversification rate estimated from phylogenies and the rate of molecular evolution. Such an association might, by influencing the relative position of node heights, systematically bias estimates of diversification time. We simulated the evolution of DNA sequences under several scenarios where rates of diversification and molecular evolution vary through time, including models where diversification and molecular evolutionary rates are linked. We show that commonly used methods, including metric‐based, likelihood and Bayesian approaches, can have a low power to identify changes in diversification rate when molecular substitution rates vary. Furthermore, the association between the rates of speciation and molecular evolution rate can cause the signature of a slowdown or speedup in speciation rates to be lost or misidentified. These results suggest that the multiple sources of variation in molecular evolutionary rates need to be considered when inferring macroevolutionary processes from phylogenies.     

Codon usage determines translation rate in Escherichia coli

We wish to determine whether differences in translation rate are correlated with differences in codon usage or with differences in mRNA secondary structure. We therefore inserted a small DNA fragment in the lacZ gene either directly or flanked by a few frame-shifting bases, leaving the reading frame of the lacZ gene unchanged. The fragment was chosen to have "infrequent" codons in one reading frame and "common" codons in the other. The insert in these constructs does not seem to give mRNAs that are able to form extensive secondary structures. The translation time for these modified lacZ mRNAs was measured with a reproducibility better than plus or minus one second. We found that the mRNA with infrequent codons inserted has an approximately three-seconds longer translation time than the one with common codons. In another set of experiments we constructed two almost identical lacZ genes in which the lacZ mRNAs have the potential to generate stem structures with stabilities of about -75 kcal/mol. In this way we could investigate the influence of mRNA structure on translation rate. This type of modified gene was generated in two reading frames with either common or infrequent codons similar to our first experiments. We find that the yield of protein from these mRNAs is reduced, probably due to the action in vivo of an RNase. Nevertheless, the data do not indicate that there is any effect of mRNA secondary structure on translation rate. In contrast, our data persuade us that there is a difference in translation rate between infrequent codons and common codons that is of the order of sixfold.     

Chloroplast DNA variation and phylogeny of theRanunculaceae

Restriction site mapping of chloroplast DNA from 31 species representing 26 genera of theRanunculaceae was performed using eleven restriction endonucleases. The chloroplast genome varies in length from approximately 152 to 160 kb. Length variants are frequent in theRanunculaceae and range from usually less than 300 bp to rarely 1.5 kb. The inverted repeat is extended into the large single copy (LSC) region by 4–4.5 kb inAnemone, Clematis, Clematopsis, Hepatica, Knowltonia, andPulsatilla. Several inversions are present in the LSC-region of the cpDNA in all these genera and inAdonis. The frequency of restriction site mutations varies within the chloroplast genome in theRanunculaceae between 4 and 32 mutations per kilobase, and is lowest in the inverted repeat and the regions containing the ATPase-genes and the genespsaA, psaB, psbA, rpoB, andrbcL. A total of 547 phylogenetically informative restriction sites was utilized in cladistic analyses of the family using Wagner, Dollo, and weighted parsimony. These three parsimony analyses result in different tree topologies. Four, six, and one equally most parsimonious trees were obtained with Wagner, Dollo, and weighted parsimony, respectively. The amount of support for the monophyletic groups was evaluated using bootstrapping and decay analysis. All three parsimony methods suggest thatHydrastis is the sister group to the remainder of theRanunculaceae, and that theAnemone-Clematis group, which shares several derived cpDNA rearrangements, is monophyletic. Only a few of the traditional groups in theRanunculaceae are supported by cpDNA restriction side data. Only Dollo parsimony provides support for the hypothesis thatThalictroideae andRanunculoideae are monophyletic.     

Intraspecific molecular phylogeny, genetic variation and phylogeography of Reticulitermes speratus (Isoptera: Rhinotermitidae)

Population structure was investigated in Reticulitermes speratus populations in the Korean Peninsula and the Japanese Archipelago. All trees derived from analyses of the combined sequence dataset of two mitochondrial genes, COII and COIII, showed that R. speratus populations cluster into two major clades comprising the Korean/southern Japanese populations and the north-ern Japanese populations. Analysis of population ge-netic structure showed strong genetic partitioning between populations of the two clades. To understand historical migration routes and current distributions, the phylogeographic history of R. speratus was inferred from intra-/interspecific phylogeny and diver-gence times estimated between the clades of the phylogenetic tree. The estimated migration route and divergence time of ancestral R. speratus are congruent with recent paleogeographic hypotheses involving land-bridge connections between the Asian continent and the Japanese Archipelago. We suggest that ancestral R. speratus separated into northern and southern Japanese populations after its migration into the Japanese main islands from East China during the early Pleistocene via the East China Sea basin, which may have been exposed during that period. The Korean populations seem to have diverged recently from southern Japanese populations; this may explain the current distribution of R. speratus in the Japanese Arachipelago, and account for why it is restricted to northern areas of the Tokara Strait.     

Mitochondrial DNA as effective molecular markers for the genetic variation and phylogeny of the family Osteoglossidae

The present study examined the genetic variation of the family Osteoglossidae from different geographical locations based on the mitochondrial NADH dehydrogenase subunit 2 (ND2) and ATPase subunit 6 (ATPase6) genes; we then re-constructed the phylogenetic relationships using the two sequences in combination. The results showed that the partial sequences of mitochondrial ND2 and ATPase6 of the family Osteoglossidae were 813 bp and 669 bp, respectively. A total of 42 species-specific nucleotide positions of the family Osteoglossidae were found to be useful for molecular identification. The sequence variation showed greater differences (8.3% ~ 28.1% for the combined sequences, 8.3% ~ 26.7% for the ND2 gene, and 9.3% ~ 28.7% for the ATPase6 gene) among the different species of Osteoglossidae, and there was a significant association between the genetic difference and geographical location. Phylogenetic analyses using neighbor-joining, Bayesian inference, and maximum parsimony (MP) methods based on the combined sequences of the two genes were able to distinguish the different species and were in agreement with the existing taxonomy based on morphological characters and in association with the geographical distribution among seven species of the family Osteoglossidae.     

Using molecular DNA markers in phylogeny and systematics

The review considers data on the use of the main evolutionary markers (ribosomal, mitochondrial, and RAPD markers; dispersed and tandem repeats). Some circumstances impending analysis of these data are discussed.     

A molecular phylogeny of annelids

We present parsimony analyses of annelids based on the largest taxon sample and most extensive molecular data set yet assembled, with two nuclear ribosomal genes (18S rDNA and the D1 region of 28S rDNA), one nuclear protein coding‐gene (Histone H3) and one mitochondrial ribosomal gene (16S rDNA) from 217 terminal taxa. Of these, 267 sequences are newly sequenced, and the remaining were obtained from GenBank. The included taxa are based on the criteria that the taxon must have 18S rDNA or at least two other loci. Our analyses show that 68% of annelid family ranked taxa represented by more than one taxon in our study are supported by a jackknife value > 50%. In spite of the size of our data set, the phylogenetic signal in the deepest part of the tree remains weak and the majority of the currently recognized major polychaete clades (except Amphinomida and Aphroditiformia) could not be recovered. Terbelliformia is monophyletic (with the exclusion of Pectinariidae, for which only 18S data were available), whereas members of taxa such as Phyllodocida, Cirratuliformia, Sabellida and Scolecida are scattered over the trees. Clitellata is monophyletic, although Dinophilidae should possibly be included, and Clitellata has a sister group within the polychaetes. One major problem is the current lack of knowledge on the closest relatives to annelids and the position of the annelid root. We suggest that the poor resolution in the basal parts of the trees presented here may be due to lack of signal connected to incomplete data sets both in terms of terminal and gene sampling, rapid radiation events and/or uneven evolutionary rates and long‐branch attraction. © The Willi Hennig Society 2006.     

A new method for characterizing replacement rate variation in molecular sequences. Application of the Fourier and wavelet models to Drosophila and mammalian proteins

We propose models for describing replacement rate variation in genes and proteins, in which the profile of relative replacement rates along the length of a given sequence is defined as a function of the site number. We consider here two types of functions, one derived from the cosine Fourier series, and the other from discrete wavelet transforms. The number of parameters used for characterizing the substitution rates along the sequences can be flexibly changed and in their most parameter-rich versions, both Fourier and wavelet models become equivalent to the unrestricted-rates model, in which each site of a sequence alignment evolves at a unique rate. When applied to a few real data sets, the new models appeared to fit data better than the discrete gamma model when compared with the Akaike information criterion and the likelihood-ratio test, although the parametric bootstrap version of the Cox test performed for one of the data sets indicated that the difference in likelihoods between the two models is not significant. The new models are applicable to testing biological hypotheses such as the statistical identity of rate variation profiles among homologous protein families. These models are also useful for determining regions in genes and proteins that evolve significantly faster or slower than the sequence average. We illustrate the application of the new method by analyzing human immunoglobulin and Drosophilid alcohol dehydrogenase sequences.     

A molecular phylogeny of bryozoans

We present the most comprehensive molecular phylogeny of bryozoans to date. Our concatenated alignment of two nuclear ribosomal and five mitochondrial genes includes 95 taxa and 13,292 nucleotide sites, of which 8297 were included. The number of new sequences generated during this project are for each gene:ssrDNA (32), lsrDNA (22), rrnL (38), rrnS (35), cox1 (37), cox3 (34), and cytb (44). Our multi-gene analysis provides a largely stable topology across the phylum. The major groups were unambiguously resolved as (Phylactolaemata (Cyclostomata (Ctenostomata, Cheilostomata))), with Ctenostomata paraphyletic. Within Phylactolaemata, (Stephanellidae, Lophopodidae) form the earliest divergent clade. Fredericellidae is not resolved as a monophyletic family and forms a clade together with Plumatellidae, Cristatellidae and Pectinatellidae, with the latter two as sister taxa. Hyalinella and Gelatinella nest within the genus Plumatella. Cyclostome taxa fall into three major clades: i. (Favosipora (Plagioecia, Rectangulata)); ii. (Entalophoroecia ((Diplosolen, Cardioecia) (Frondipora, Cancellata))); and iii. (Articulata ((Annectocyma, Heteroporidae) (Tubulipora (Tennysonia, Idmidronea)))), with suborders Tubuliporina and Cerioporina, and family Plagioeciidae each being polyphyletic. Ctenostomata is composed of three paraphyletic clades to the inclusion of Cheilostomata: ((Alcyonidium, Flustrellidra) (Paludicella (Anguinella, Triticella)) (Hislopia (Bowerbankia, Amathia)) Cheilostomata); Flustrellidra nests within the genus Alcyonidium, and Amathia nests within the genus Bowerbankia. Suborders Carnosa and Stolonifera are not monophyletic. Within the cheilostomes, Malacostega is paraphyletic to the inclusion of all other cheilostomes. Conopeum is the most early divergent cheilostome, forming the sister group to ((Malacostega, Scrupariina, Inovicellina) ((Hippothoomorpha, Flustrina) (Lepraliomorpha, Umbonulomorpha))); Flustrina is paraphyletic to the inclusion of the hippothoomorphs; neither Lepraliomorpha nor Umbonulomorpha is monophyletic. Ascophorans are polyphyletic, with hippothoomorphs grouping separately from lepraliomorphs and umbonulomorphs; no cribrimorphs were included in the analysis. Results are discussed in the light of molecular and morphological evidence. Ancestral state reconstruction of larval strategy in Gymnolaemata revealed planktotrophy and lecithotrophy as equally parsimonious solutions for the ancestral condition. More comprehensive taxon sampling is expected to clarify this result. We discuss the extent of non-bryozoan contaminant sequences deposited in GenBank and their impact on the reconstruction of metazoan phylogenies and those of bryozoan interrelationships.     

An examination of phylogenetic models of substitution rate variation among lineages

Molecular evolutionary rates can show significant variation among lineages, complicating the task of estimating substitution rates and divergence times using phylogenetic methods. Accordingly, relaxed molecular clock models have been developed to accommodate such rate heterogeneity, but these often make the assumption of rate autocorrelation among lineages. In this paper, I examine the validity of this assumption.     

Nucleotide substitution models and estimation of phylogeny

The nucleotide substitution matrix inferred from avian data sets using cytochrome b differs considerably from the models commonly used in phylogenetic analyses. To analyze the possible effects of this particular pattern of change in phylogeny estimation we performed a computer simulation in which we started with a real sequence and used the inferred model of change to produce a tree of 10 species. Maximum parsimony (MP), maximum likelihood (ML), and various distance methods were then used to recover the topology and the branch lengths. We used two kinds of data with varying levels of variation. In addition, we tested with the removal of third positions and different weighting schemes. At low levels of variation, MP was outstanding in recovering the topology (90% correct), while unweighted pair-group method, arithmetic average (UPGMA), regardless of distances used, was poor (40%). At the higher level, most methods had a chance of around 40%-58% of finding the true tree. However, in most cases, the trees found were only slightly wrong, with only one or a few branches misplaced. On the other hand, the use of a "wrong" model had serious effects on the estimation of branch lengths (distances). Although precision was high, accuracy was poor with most methods, giving branch lengths that were biased downward. When seeded with the true distance matrix, Fitch and NJ always found the true tree, while UPGMA frequently failed to do so. The effect of removing third positions was dramatic at low levels of variation, because only one MP program was able to find a true tree at all, albeit rarely, while none of the others ever did so. At higher levels, the situation was better, but still much worse than with the whole data set.     

A molecular phylogeny and classification of Bignoniaceae

Bignoniaceae are woody, trees, shrubs, and lianas found in all tropical floras of the world with lesser representation in temperate regions. Phylogenetic analyses of chloroplast sequences (rbcL, ndhF, trnL-F) were undertaken to infer evolutionary relationships in Bignoniaceae and to revise its classification. Eight clades are recognized as tribes (Bignonieae, Catalpeae, Coleeae, Crescentieae, Jacarandeae, Oroxyleae, Tecomeae, Tourrettieae); additional inclusive clades are named informally. Jacarandeae and Catalpeae are resurrected; the former is sister to the rest of the family, and the latter occupies an unresolved position within the "core" Bignoniaceae. Tribe Eccremocarpeae is included in Tourrettieae. Past classifications recognized a large Tecomeae, but this tribe is paraphyletic with respect to all other tribes. Here Tecomeae are reduced to a clade of approximately 12 genera with a worldwide distribution in both temperate and tropical ecosystems. Two large clades, Bignonieae and Crescentiina, account for over 80% of the species in the family. Coleeae and Crescentieae are each included in larger clades, the Paleotropical alliance and Tabebuia alliance, respectively; each alliance includes a grade of taxa assigned to the traditional Tecomeae. Parsimony inference suggests that the family originated in the neotropics, with at least five dispersal events leading to the Old World representatives.     

A molecular phylogeny and classification of Verbenaceae

? Premise of the study: Verbenaceae consist of trees, shrubs, lianas, and herbs distributed primarily in Latin America, where they occur in a wide array of ecosystems. A second center of diversity exists in Africa. Competing morphology-based classifications that rely on different traits conflict in significant ways. A broad phylogenetic study was undertaken to assess those classifications and to examine the historical geography of the family. ? Methods: Analysis of seven chloroplast DNA regions for 109 species, representing all genera except one monotypic genus, provide inference into evolutionary relationships in Verbenaceae. ? Key results: The phylogeny shows that none of the traditional classifications reflect phylogenetic relationships very well. Eight clades are recognized as tribes (Casselieae, Citharexyleae, Duranteae, Lantaneae, Neospartoneae trib. nov., Petreeae, Priveae, and Verbeneae). Two genera, Dipyrena and Rhaphithamnus, remain unplaced in these larger clades. Petreeae, which consist of Neotropical lianas, are sister to the rest of the family. Lantaneae and Verbeneae together form a derived clade that comprises approximately two-thirds of the species in Verbenaceae. ? Conclusions: We present a new tribal classification, including one new tribe, Neospartoneae trib. nov., to accommodate three small genera of Argentine species (Diostea, Neosparton, and Lampaya). Phylogenetic inference suggests a South American origin for Verbenaceae, with approximately six colonization events having given rise to the Old World species.     

Invariable sites models and their use in phylogeny reconstruction

Phylogenetic inference is well known to be problematic if both long and short branches occur together in the underlying tree. With biological data, correcting for this problem may require simultaneous consideration for both substitution biases and rate heterogeneity between lineages and across sequence positions. A particular form of the latter is the presence of invariable sites, which are well known to mislead estimation of genetic divergences. Here we describe a capture-recapture method to estimate the proportion of invariable sites in an alignment of amino acids or nucleotides. We use it to investigate phylogenetic signals in 18S ribosomal DNA sequences from Holometabolus insects. Our results suggest that, as taxa diverged, their 18S rDNA sequences have altered in both their distribution of sites that can vary as well as in their base compositions.     

Stochastic models of molecular evolution and the estimation of phylogeny and rates of nucleotide substitution in the hominoid primates

We have estimated phylogenetic patterns and rates of nucleotide substitution in the hominoid primates using two different probabilistic models of molecular evolution as applied to three different data sets of nucleic acid sequences. The orang-utan was found to be the out-group of the other hominoids examined. Within the African apes and human clade the sister-group relationship of chimpanzee and human was found to be statistically the best, although the magnitude of the error estimates (a reflection of random statistical fluctuations) makes this conclusion tentative. The ψν-globin data sets were found to be statistically the most consistent and gave estimates of the times of divergence of chimpanzee and human from gorilla and of chimpanzee from human as 7·7 ± 1·5 Ma (Millions of years ago) and 7·4 ± 1·5 Ma respectively, although the speculative nature of these estimates is emphasized. In all cases the calibration point was the assumed divergence of the orang-utan from the remaining hominoids at 14·5 Ma. There was no statistically significant evidence of a slowdown in nucleotide substitution rate for the human lineage, or among the hominoids as a whole with respect to the Old and New World monkeys. We advocate the continued use and development of stochastic models of molecular evolution as a basis for phylogenetic estimation. On this basis one can choose between competing hypotheses of relationship in a statistical manner and can provide estimates of the errors involved in such estimations. The assumptions of all stochastic models are open to test and future refinement.     

Gene copy number variation and its significance in cyanobacterial phylogeny

ABSTRACT: BACKGROUND: In eukaryotes, variation in gene copy numbers is often associated with deleterious effects, but may also have positive effects. For prokaryotes, studies on gene copy number variation are rare. Previous studies have suggested that high numbers of rRNA gene copies can be advantageous in environments with changing resource availability, but further association of gene copies and phenotypic traits are not documented. We used one of the morphologically most diverse prokaryotic phyla to test whether numbers of gene copies are associated with levels of cell differentiation. RESULTS: We implemented a search algorithm that identified 44 genes with highly conserved copies across 22 fully sequenced cyanobacterial taxa. For two very basal cyanobacterial species, Gloeobacter violaceus and a thermophilic Synechococcus species, distinct phylogenetic positions previously found were supported by identical protein coding gene copy numbers. Furthermore, we found that increased ribosomal gene copy numbers showed a strong correlation to cyanobacteria capable of terminal cell differentiation. Additionally, we detected extremely low variation of 16S rRNA sequence copies within the cyanobacteria. We compared our results for 16S rRNA to three other eubacterial phyla (Chroroflexi, Spirochaetes and Bacteroidetes). Based on Bayesian phylogenetic inference and the comparisons of genetic istances, we could confirm that cyanobacterial 16S rRNA paralogs and orthologs show significantly stronger conservation than found in other eubacterial phyla. Conclusions: A higher number of ribosomal operons could potentially provide an advantage to terminally differentiated cyanobacteria. Furthermore, we suggest that 16S rRNA gene copies in cyanobacteria are homogenized by both concerted evolution and purifying selection. In addition, the small ribosomal subunit in cyanobacteria appears to evolve at extraordinary slow evolutionary rates, an observation that has been made previously for morphological characteristics of cyanobacteria.     

Codon substitution models based on residue similarity and their applications

Codon models are now widely used to draw evolutionary inferences from alignments of homologous sequence data. Incorporating physicochemical properties of amino acids into codon models, two novel codon substitution models describing the evolution of protein-coding DNA sequences are presented based on the similarity scores of amino acids. To describe substitutions between codons a continue-time Markov process is used. Transition/transversion rate bias and nonsynonymous codon usage bias are allowed in the models. In our implementation, the parameters are estimated by maximum-likelihood (ML) method as in previous studies. Furthermore, instantaneous mutations involving more than one nucleotide position of a codon are considered in the second model. Then the two suggested models are applied to five real data sets. The analytic results indicate that the new codon models considering physicochemical properties of amino acids can provide a better fit to the data comparing with existing codon models, and then produce more reliable estimates of certain biologically important measures than existing methods.     

Morphological and molecular phylogeny of Pistacia species in Turkey

This study addresses the taxonomic relationships and genetic variation of wild Pistacia germplasm in Turkey using morphological data and RAPD analysis. P. atlantica, P. terebinthus and P. eurycarpa are common wild species in the flora of Turkey, and their phenotypic appearance and productivity are highly variable. Understanding such variation would facilitate their use in rootstock breeding programs as rootstock for edible pistachio. We have sampled and characterized a total of 40 wild Pistacia genotypes from different parts of Turkey for this study. These included 10 P. eurycarpa (locally identified as P. khinjuk) genotypes from Siirt and Gaziantep provinces and 20 P. atlantica and 10 P. terebinthus genotypes from Adana, Aydin and Manisa provinces. In addition, two local P. vera varieties, cvs. Kirmizi and Siirt, were added for comparison. Cluster analysis based on morphological data revealed that the closest species to P. vera is P. eurycarpa, followed by P. atlantica and P. terebinthus. Ten polymorphic RAPD primers, yielding a total of 138 scorable bands, were selected and used for DNA fingerprinting of these genotypes. In the resulting molecular phylogeny, the four Pistacia species are clearly separated from each other. P. terebinthus appears to be the most diverged species, and the closest pair of species was found to be P. atlantica and P. eurycarpa. This supported the classification of trees that had been identified by local growers as P. khinjuk, as P. eurycarpa. Comparison of these samples with a reference P. khinjuk, obtained from a germplasm collection in the USA, also supported such identification. Several wild genotypes were found to be inter-specific hybrids, and the RAPD patterns revealed their probable origin. Species-specific markers were identified for each of the four species, and these may aid in future classification of new germplasm materials.