Multilocus Phylogeny of Lycodon and the Taxonomic Revision of Oligodon multizonatum

Classification of the Asian snake genera Lycodon and Oligodon has proven challenging. We conducted a molecular phylogenetic analysis to estimate the phylogenetic relationships in the genus of Lycodon and clarify the taxonomic status of Oligodon multizonatum using mitochondrial（cyt b, ND4） and nuclear（c-mos） genes. Phylogenetic trees estimated using Maximum Likelihood and Bayesian Inference indicated that O. multizonatum is actually a species of Lycodon. Comparing morphological data from O. multizonatum and its closest relatives also supported this conclusion. Our results imply that a thorough review of the evolutionary relationships in the genus of Lycodon is strong suggested.     

Revaluation of deuterostome phylogeny and evolutionary relationships among chordate subphyla using mitogenome data

The traditional knowledge in textbooks indicated that cephalochordates were the closest relatives to vertebrates among all extant organisms. However, this opinion was challenged by several recent phylogenetic studies using hundreds of nuclear genes. The researchers suggested that urochordates, but not cephalochordates, should be the closest living relatives to vertebrates. In the present study, by using data generated from hundreds of mtDNA sequences, we revalue the deuterostome phylogeny in terms of whole mitochondrial genomes （mitogenomes）. Our results firmly demonstrate that each of extant deuterostome phyla and chordate subphyla is monophyletic. But the results present several alternative phylogenetic trees depending on different sequence datasets used in the analysis. Although no clear phylogenetic relationships are obtained, those trees indicate that the ancient common ancestor diversified rapidly soon after their appearance in the early Cambrian and generated all major deuterostome lineages during a short historical period, which is consistent with ＂Cambrian explosion＂ revealed by paleontologists. It was the 520-million-year＇s evolution that obscured the phylogenetic relationships of extant deuterostomes. Thus, we conclude that an integrative analysis approach rather than simply using more DNA sequences should be employed to address the distant evolutionary relationship.     

基于28S rRNA基因D2序列的优茧蜂亚科分子系统发育（膜翅目：茧蜂科）

The D2 variable region of 28S ribosomal RNA was sequenced from ethanol specimens or obtained from the literature to provide the first phylogenetie reconstruction of the subfamily Euphorinae (Hymenoptera;Braconidae). Phylogenetic relationships were established by comparing the results using two different methods (distance-based neighbor-joining, NJ; and maximum parsimony, MP) and three different outgroups. The monophyly of the Euphorinae is well supported by all trees generated from molecular data. All phylogenetic reconstructions yielded trees with very similar topologies that only partially resolved the morphologically defined tribes and the relationships within the subfamily. We found no evidence for the monophyletic natures of the tribes Euphorinl, Dinocampini,Perilitini, Syntretini, Comsophorini and Centisitini, but we did find some evidence for the tribes Meteorini and Microctonini. The monophyletic nature of the tribe Meteodnl was well-supported in all trees. We also found the clade containing the LecythodeUa,Microctonus, Orionis and Streblocera to be a monophyletic group, which corresponded to the tribe Microtonini, with Orionis transferred from the tribe Eupholini into Microtonini.Among the genera of Euphorini our results showed strong support for a paraphyletic nature of this group, which can be roughly divided into two clades, one consisting of Aridelus Wesmaelia, the other of Leiophron Peristenus, suggesting both of which may be given tribal rank. The placement of the genus Chrysopophorus is largely uncertain. Two clades,Dinocampus Perilitus and Cosmophorus Rhopalophorus, were constantly resolved in our analyses, with 42-96 and 97-100 bootstrap value support, respectively, suggesting that both of them form monophyletic groups. For members of the Centistini, Pygostolus may be removed and included in Microctonini or other relative tribe.     

Bacterial microbiome of Coptotermes curvignathus （Isoptera： Rhinotermitidae） reflects the coevolution of species and dietary pattern

Coptotermes curvignathus Holmgren is capable of feeding on living trees. This ability is attributed to their effective digestive system that is furnished by the termite＇s own cellulolytic enzymes and cooperative enzymes produced by their gut microbes. In this study, the identity of an array of diverse microbes residing in the gut of C. curvignathus was revealed by sequencing the near-full-length 16S rRNA genes. A total of 154 bacterial phylotypes were found. The Bacteroidetes was the most abundant phylum and accounted for about 65% of the gut microbial profile. This is followed by Firmicutes, Actinobacteria, Spirochetes, Proteobacteria, TM7, Deferribacteres, Planctomycetes, Ferrucomicrobia, and Termite Group 1. Based on the phylogenetic study, this symbiosis can be a result of long coevolution of gut enterotypes with the phylogenic distribution, strong selection pressure in the gut, and other speculative pressures that determine bacterial biome to follow. The phylogenetic distribution of cloned rRNA genes in the bacterial domain that was considerably different from other termite reflects the strong selection pressures in the gut where a proportional composition of gut microbiome of C. curvignathus has established. The selection pressures could be linked to the unique diet preference of C. curvignathus that profoundly feeds on living trees. The delicate gut microbiome composition may provide available nutrients to the host as well as potential protection against opportunistic pathogen.     

Whole genome duplication of intra- and inter-chromosomes in the tomato genome

Whole genome duplication(WGD)events have been proven to occur in the evolutionary history of most angiosperms.Tomato is considered a model species of the Solanaceae family.In this study,we describe the details of the evolutionary process of the tomato genome by detecting collinearity blocks and dating the WGD events on the tree of life by combining two different methods:synonymous substitution rates(Ks)and phylogenetic trees.In total,593 collinearity blocks were discovered out of 12 pseudo-chromosomes constructed. It was evident that chromosome 2 had experienced an intra-chromosomal duplication event.Major inter-chromosomal duplication occurred among all the pseudo-chromosome.We calculated the Ks value of these collinearity blocks.Two peaks of Ks distribution were found,corresponding to two WGD events occurring approximately 36-82 million years ago(MYA)and 148-205 MYA.Additionally, the results of phylogenetic trees suggested that the more recent WGD event may have occurred after the divergence of the rosidasterid clade,but before the major diversification in Solanaceae.The older WGD event was shown to have occurred before the divergence of the rosid-asterid clade and after the divergence of rice-Arabidopsis(monocot-dicot).     

Phylogenomic analyses reveal intractable evolutionary history of a temperate bamboo genus (Poaceae: Bambusoideae)

Shibataea is a genus of temperate bamboos(Poaceae:Bambusoideae)endemic to China,but little is known about its phylogenetic position and interspecific relationships.To elucidate the phylogenetic relationship of the bamboo genus Shibataea,we performed genome-scale phylogenetic analysis of all seven species and one variety of the genus using double digest restriction-site associated DNA sequencing(dd RAD-seq)and whole plastid genomes generated using genome skimming.Our phylogenomic analyses based on dd RAD-seq and plastome data congruently recovered Shibataea as monophyletic.The nuclear data resolved S.hispida as the earliest diverged species,followed by S.chinensis,while the rest of Shibataea can be further divided into two clades.However,the plastid and nuclear topologies conflict significantly.By comparing the results of network analysis and topologies reconstructed from different datasets,we identify S.kumasasa as the most admixed species,which may be caused by incomplete lineage sorting(ILS)or interspecific gene flow with four sympatric species.This study highlights the power of dd RAD and plastome data in resolving complex relationships in the intractable bamboo genus.     

Analysis of plastid and nuclear DNA data in plant phylogenetics——evaluation and improvement

Correct combination of plastid(cp)and nuclear(nr)DNA data for plant phylogenetic reconstructions is not a new issue,but with an increasing number of nrDNA loci being used,it is of ever greater practical concern.For accurately reconstructing the phylogeny and evolutionary history of plant groups,correct treatment of phylogenetic incongruence is a vital step in the proper analysis of cpDNA and nrDNA data.We first evaluated the current status of analyzing cpDNA and nrDNA data by searching all articles published in the journal Systematic Botany between 2005 and 2011.Many studies combining cpDNA and nrDNA data did not rigorously assess the combinability of the data sets,or did not address in detail possible reasons for incongruence between the two data sets.By reviewing various methods,we outline a procedure to more accurately analyze and/or combine cpDNA and nrDNA data,which includes four steps:identifying significant incongruence,determining conflicting taxa,providing possible interpretations for incongruence,and reconstructing the phylogeny after treating incongruence.Particular attention is given to explanation of the cause of incongruence.We hope that our procedure will help raise awareness of the importance of rigorous analysis and help identify the cause of incongruence before combining cpDNA and nrDNA data.     

Phylogenetic positions of several amitochondriate protozoa

Several groups of parasitic protozoa, as represented by Giardia, Trichomonas, Entamoeba and Microsporida, were once widely considered to be the most primitive extant eukaryotic group―Archezoa. The main evidence for this is their ‘lacking mitochondria’ and possessing some other primitive features between prokaryotes and eukaryotes, and being basal to all eukaryotes with mitochondria in phylogenies inferred from many molecules. Some authors even proposed that these organisms diverged before the endosymbiotic origin of mitochondria within eukaryotes. This view was once considered to be very significant to the study of origin and evolution of eukaryotic cells (eukaryotes). However, in recent years this has been challenged by accumulating evidence from new studies. Here the sequences of DNA topoisomerase II in G. lamblia, T. vaginalis and E. histolytica were identified first by PCR and sequencing, then combining with the sequence data of the microsporidia Encephalitozoon cunicul and other eukaryotic groups of different evolutionary positions from GenBank, phylogenetic trees were constructed by various methods to investigate the evolutionary positions of these amitochondriate protozoa. Our results showed that since the characteristics of DNA topoisomerase II make it avoid the defect of ‘long-branch attraction’ appearing in the previous phylogenetic analyses, our trees can not only reflect effectively the relationship of different major eukaryotic groups, which is widely accepted, but also reveal phylogenetic positions for these amitochondriate protozoa, which is different from the previous phylogenetic trees. They are not the earliest-branching eukaryotes, but diverged after some mitochondriate organisms such as kinetoplastids and mycetozoan; they are not a united group but occupy different phylogenetic positions. Combining with the recent cytological findings of mitochondria-like organelles in them, we think that though some of them (e.g. diplo-monads, as represented by Giardia) may occupy a very low evolutionary position, generally these organisms are not as extremely primitive as was thought before; they should be poly-phyletic groups diverging after the endosymbiotic origin of mitochondrion to adapt themselves to anaerobic parasitic life.     

Estimating ancestral distributions of lineages with uncertain sister groups： a statistical approach to Dispersal-Vicariance Analysis and a case using Aesculus L. （Sapindaceae） including fossils

We propose a simple statistical approach for using Dispersal-Vicariance Analysis （DIVA） software to infer biogeographic histories without fully bifurcating trees. In this approach, ancestral ranges are first optimized for a sample of Bayesian trees. The probability P of an ancestral range r at a node is then calculated as P（rY） = ∑t^n=1 F（rY）t Pt where Y is a node, and F（rY） is the frequency of range r among all the optimal solutions resulting from DIVA optimization at node Y, t is one of n topologies optimized, and Pt is the probability of topology t. Node Y is a hypothesized ancestor shared by a specific crown lineage and the sister of that lineage ＂x＂, where x may vary due to phylogenetic uncertainty （polytomies and nodes with posterior probability 〈 100%）. Using this method, the ancestral distribution at Y can be estimated to provide inference of the geographic origins of the specific crown group of interest. This approach takes into account phylogenetic uncertainty as well as uncertainty from DIVA optimization. It is an extension of the previously described method called Bayes-DIVA, which pairs Bayesian phylogenetic analysis with biogeographic analysis using DIVA. Further, we show that the probability P of an ancestral range at Y calculated using this method does not equate to pp＊F（rY） on the Bayesian consensus tree when both variables are 〈 100%, where pp is the posterior probability and F（rY） is the frequency of range r for the node containing the specific crown group. We tested our DIVA-Bayes approach using Aesculus L., which has major lineages unresolved as a polytomy. We inferred the most probable geographic origins of the five traditional sections of Aesculus and ofAesculus californica Nutt. and examined range subdivisions at parental nodes of these lineages. Additionally, we used the DIVA-Bayes data from Aesculus to quantify the effects on biogeographic inference of including two wildcard fossil taxa in phylogenetic analysis. Our analysis resolved the geographic     

Molecular phylogeny of Soims-laubachia （Brassicaceae） s.l., based on multiple nuclear and plastid DNA sequences,and its biogeographic implications

The Hengduan Mountains region of south-west China is a noted biodiversity,hotspot, but the geographic origins and historical assembly of its rich endemic flora, including the sky-island species ofSolms-laubachia Muschl. （Brassicaceae）, have been little studied. Previous molecular studies on the phylogeny of Solms-laubachia showed it to be paraphyletic, leading to considerable expansion not only of its taxonomic limits, but also its geographic range, with the inclusion of taxa from outside the Hengduan region. However, these studies provided little resolution of interspecific relationships, preventing inferences about historical biogeography within the clade. In the present study, new sequence data from two nuclear genes （LEAFY and G3pdh） and two chloroplast intergenic spacers （petN-psbM and psbM-trnD） were combined with existing markers to increase phylogenetic signals. Phaeonychium villosum （Maxim.） Al-Shehbaz was found to be nested within Solms-laubachia s.l. In general, phylogenetic relationships appear to be a good predictor of geography, with the Hengduan Mountain endemics embedded in a paraphyletic grade of species from the western Himalayas and central Asia, but they also imply morphological homoplasy, lncongruence was detected between the nuclear and chloroplast gene trees, perhaps resulting from incomplete lineage sorting of ancestral polymorphisms. The crown age ofSolms-laubachia s.l. was estimated to be approximately 1.42-3.68 mya, using Bayesian relaxed molecular clock analysis. Historical biogeographic analysis using a parametric dispersalextinction-cladogenesis model inferred central Asia and the western Himalayas as most probable ancestral range of Solms-laubachia s.l., and estimated higher rates of eastward expansion than westward during the diversification of descendant lineages. In summary, our results suggest that Solms-laubachia s.l. originated during the Pliocene in central Asia, and subsequently migrated eastward into the Hengduan Mountains, colonizing sky-island, alpine scre     

BIMLR: A method for constructing rooted phylogenetic networks from rooted phylogenetic trees

Rooted phylogenetic trees constructed from different datasets (e.g. from different genes) are often conflicting with one another, i.e. they cannot be integrated into a single phylogenetic tree. Phylogenetic networks have become an important tool in molecular evolution, and rooted phylogenetic networks are able to represent conflicting rooted phylogenetic trees. Hence, the development of appropriate methods to compute rooted phylogenetic networks from rooted phylogenetic trees has attracted considerable research interest of late. The C_ASS algorithm proposed by van Iersel et al. is able to construct much simpler networks than other available methods, but it is extremely slow, and the networks it constructs are dependent on the order of the input data. Here, we introduce an improved C_ASS algorithm, BIMLR. We show that BIMLR is faster than C_ASS and less dependent on the input data order. Moreover, BIMLR is able to construct much simpler networks than almost all other methods. BIMLR is available at http://nclab.hit.edu.cn/wangjuan/BIMLR/.     

SuperFine: fast and accurate supertree estimation

Many research groups are estimating trees containing anywhere from a few thousands to hundreds of thousands of species, toward the eventual goal of the estimation of a Tree of Life, containing perhaps as many as several million leaves. These phylogenetic estimations present enormous computational challenges, and current computational methods are likely to fail to run even on data sets in the low end of this range. One approach to estimate a large species tree is to use phylogenetic estimation methods (such as maximum likelihood) on a supermatrix produced by concatenating multiple sequence alignments for a collection of markers; however, the most accurate of these phylogenetic estimation methods are extremely computationally intensive for data sets with more than a few thousand sequences. Supertree methods, which assemble phylogenetic trees from a collection of trees on subsets of the taxa, are important tools for phylogeny estimation where phylogenetic analyses based upon maximum likelihood (ML) are infeasible. In this paper, we introduce SuperFine, a meta-method that utilizes a novel two-step procedure in order to improve the accuracy and scalability of supertree methods. Our study, using both simulated and empirical data, shows that SuperFine-boosted supertree methods produce more accurate trees than standard supertree methods, and run quickly on very large data sets with thousands of sequences. Furthermore, SuperFine-boosted matrix representation with parsimony (MRP, the most well-known supertree method) approaches the accuracy of ML methods on supermatrix data sets under realistic conditions.     

DendroBLAST: Approximate Phylogenetic Trees in the Absence of Multiple Sequence Alignments

The rapidly growing availability of genome information has created considerable demand for both fast and accurate phylogenetic inference algorithms. We present a novel method called DendroBLAST for reconstructing phylogenetic dendrograms/trees from protein sequences using BLAST. This method differs from other methods by incorporating a simple model of sequence evolution to test the effect of introducing sequence changes on the reliability of the bipartitions in the inferred tree. Using realistic simulated sequence data we demonstrate that this method produces phylogenetic trees that are more accurate than other commonly-used distance based methods though not as accurate as maximum likelihood methods from good quality multiple sequence alignments. In addition to tests on simulated data, we use DendroBLAST to generate input trees for a supertree reconstruction of the phylogeny of the Archaea. This independent analysis produces an approximate phylogeny of the Archaea that has both high precision and recall when compared to previously published analysis of the same dataset using conventional methods. Taken together these results demonstrate that approximate phylogenetic trees can be produced in the absence of multiple sequence alignments, and we propose that these trees will provide a platform for improving and informing downstream bioinformatic analysis. A web implementation of the DendroBLAST method is freely available for use at http://www.dendroblast.com/.     

Recovering speciation and extinction dynamics based on phylogenies

Phylogenetic trees of only extant species contain information about the underlying speciation and extinction pattern. In this review, I provide an overview over the different methodologies that recover the speciation and extinction dynamics from phylogenetic trees. Broadly, the methods can be divided into two classes: (i) methods using the phylogenetic tree shapes (i.e. trees without branch length information) allowing us to test for speciation rate variation and (ii) methods using the phylogenetic trees with branch length information allowing us to quantify speciation and extinction rates. I end the article with an overview on limitations, open questions and challenges of the reviewed methodology.     

Encoding phylogenetic trees in terms of weighted quartets

One of the main problems in phylogenetics is to develop systematic methods for constructing evolutionary or phylogenetic trees. For a set of species X, an edge-weighted phylogenetic X-tree or phylogenetic tree is a (graph theoretical) tree with leaf set X and no degree 2 vertices, together with a map assigning a non-negative length to each edge of the tree. Within phylogenetics, several methods have been proposed for constructing such trees that work by trying to piece together quartet trees on X, i.e. phylogenetic trees each having four leaves in X. Hence, it is of interest to characterise when a collection of quartet trees corresponds to a (unique) phylogenetic tree. Recently, Dress and Erdös provided such a characterisation for binary phylogenetic trees, that is, phylogenetic trees all of whose internal vertices have degree 3. Here we provide a new characterisation for arbitrary phylogenetic trees.     

Student construction of phylogenetic trees in an introductory biology course

Background

Phylogenetic trees have become increasingly essential across biology disciplines. Consequently, learning about phylogenetic trees has become an important component of biology education and an area of interest for biology education research. Construction tasks, in which students generate phylogenetic trees from some type of data, are often used for instruction. However, the impact of these exercises on student learning is uncertain, in part due to our fragmented knowledge of what students construct during the tasks. The goal of this project was to develop a more robust method for describing student-generated phylogenetic trees, which will support future investigations that attempt to link construction tasks with student learning.

Results

Through iterative examination of data from an introductory biology course, we developed a method for describing student-generated phylogenetic trees in terms of style, conventionality, and accuracy. Students used the diagonal style more often than the bracket style for construction tasks. The majority of phylogenetic trees were constructed conventionally, and variable orientation of branches was the most common unconventional feature. In addition, the majority of phylogenetic trees were generated correctly (no errors) or adequately (minor errors only) in terms of accuracy. Suggesting extant taxa are descended from other extant taxa was the most common major error, while empty branches and extra nodes were very common minor errors.

Conclusions

The method we developed to describe student-constructed phylogenetic trees uncovered several trends that warrant further investigation. For example, while diagonal and bracket phylogenetic trees contain equivalent information, student preference for using the diagonal style could impact comprehension. In addition, despite a lack of explicit instruction, students generated phylogenetic trees that were largely conventional and accurate. Surprisingly, accuracy and conventionality were also dependent on each other. Our method for describing phylogenetic trees constructed by students is based on data from one introductory biology course at one institution, and the results are likely limited. We encourage researchers to use our method as a baseline for developing a more generalizable tool, which will support future investigations that attempt to link construction tasks with student learning.

     

Multilocus analyses of an Antarctic fish species flock (Teleostei, Notothenioidei, Trematominae): phylogenetic approach and test of the early-radiation event

Clades that have undergone episodes of rapid cladogenesis are challenging from a phylogenetic point of view. They are generally characterised by short or missing internal branches in phylogenetic trees and by conflicting topologies among individual gene trees. This may be the case of the subfamily Trematominae, a group of marine teleosts of coastal Antarctic waters, which is considered to have passed through a period of rapid diversification. Despite much phylogenetic attention, the relationships among Trematominae species remain unclear. In contrast to previous studies that were mostly based on concatenated datasets of mitochondrial and/or single nuclear loci, we applied various single-locus and multilocus phylogenetic approaches to sequences from 11 loci (eight nuclear) and we also used several methods to assess the hypothesis of a radiation event in Trematominae evolution. Diversification rate analyses support the hypothesis of a period of rapid diversification during Trematominae history and only a few nodes in the hypothetical species tree were consistently resolved with various phylogenetic methods. We detected significant discrepancies among trees from individual genes of these species, most probably resulting from incomplete lineage sorting, suggesting that concatenation of loci is not the most appropriate way to investigate Trematominae species interrelationships. These data also provide information about the possible effects of historic climate changes on the diversification rate of this group of fish.     

Understanding phylogenetic incongruence: lessons from phyllostomid bats

All characters and trait systems in an organism share a common evolutionary history that can be estimated using phylogenetic methods. However, differential rates of change and the evolutionary mechanisms driving those rates result in pervasive phylogenetic conflict. These drivers need to be uncovered because mismatches between evolutionary processes and phylogenetic models can lead to high confidence in incorrect hypotheses. Incongruence between phylogenies derived from morphological versus molecular analyses, and between trees based on different subsets of molecular sequences has become pervasive as datasets have expanded rapidly in both characters and species. For more than a decade, evolutionary relationships among members of the New World bat family Phyllostomidae inferred from morphological and molecular data have been in conflict. Here, we develop and apply methods to minimize systematic biases, uncover the biological mechanisms underlying phylogenetic conflict, and outline data requirements for future phylogenomic and morphological data collection. We introduce new morphological data for phyllostomids and outgroups and expand previous molecular analyses to eliminate methodological sources of phylogenetic conflict such as taxonomic sampling, sparse character sampling, or use of different algorithms to estimate the phylogeny. We also evaluate the impact of biological sources of conflict: saturation in morphological changes and molecular substitutions, and other processes that result in incongruent trees, including convergent morphological and molecular evolution. Methodological sources of incongruence play some role in generating phylogenetic conflict, and are relatively easy to eliminate by matching taxa, collecting more characters, and applying the same algorithms to optimize phylogeny. The evolutionary patterns uncovered are consistent with multiple biological sources of conflict, including saturation in morphological and molecular changes, adaptive morphological convergence among nectar‐feeding lineages, and incongruent gene trees. Applying methods to account for nucleotide sequence saturation reduces, but does not completely eliminate, phylogenetic conflict. We ruled out paralogy, lateral gene transfer, and poor taxon sampling and outgroup choices among the processes leading to incongruent gene trees in phyllostomid bats. Uncovering and countering the possible effects of introgression and lineage sorting of ancestral polymorphism on gene trees will require great leaps in genomic and allelic sequencing in this species‐rich mammalian family. We also found evidence for adaptive molecular evolution leading to convergence in mitochondrial proteins among nectar‐feeding lineages. In conclusion, the biological processes that generate phylogenetic conflict are ubiquitous, and overcoming incongruence requires better models and more data than have been collected even in well‐studied organisms such as phyllostomid bats.     

Evolutionary change of restriction cleavage sites and phylogenetic inference for man and apes

A mathematical theory for the evolutionary change of restriction endonuclease cleavage sites is developed, and the probabilities of various types of restriction-site changes are evaluated. A computer simulation is also conducted to study properties of the evolutionary change of restriction sites. These studies indicate that parsimony methods of constructing phylogenetic trees often make erroneous inferences about evolutionary changes of restriction sites unless the number of nucleotide substitutions per site is less than 0.01 for all branches of the tree. This introduces a systematic error in estimating the number of mutational changes for each branch and, consequently, in constructing phylogenetic trees. Therefore, parsimony methods should be used only in cases where nucleotide sequences are closely related. Reexamination of Ferris et al.'s data on restriction-site differences of mitochondrial DNAs does not support Templeton's conclusions regarding the phylogenetic tree for man and apes and the molecular clock hypothesis. Templeton's claim that Nei and Li's method of estimating the number of nucleotide substitutions per site is seriously affected by parallel losses and loss-gains of restriction sites is also unsupported.