Relationships of Nematodirus species and Nematodirus battus isolates (Nematoda: Trichostrongyloidea) based on nuclear ribosomal DNA sequences

Nuclear ribosomal sequence data from the internal transcribed spacers (ITS-1 and ITS-2), 5.8S subunit, and regions of the 18S and 28S genes were used to investigate sequence diversity among geographic samples of Nematodirus battus, and to infer phylogenetic relationships among Nematodirus species. Phylogenetic analysis of these data yielded strong support for relationships among species, depicting Nematodirus helvetianus and Nematodirus spathiger as sister-taxa and a clade of these 2 species and Nematodirus filicollis. This tree is consistent with caprine bovids as ancestral hosts, with a subsequent host shift to Bovinae in N. helvetianus. Eleven of 14 N. battus sequences were unique, with 19 variable sites among sequences representing 5 geographic samples. The lowest number of variable nucleotide sites was observed in samples representing apparently recent introductions to the United States and Canada, which is consistent with a population bottleneck concomitant with translocation. Comparison of directly sequenced polymerase chain reaction products and clones revealed evidence for intraindividual variation at some of the sequence sites, and this pattern of variation and that within geographic samples indicates incomplete rDNA repeat homogenization within species. This pattern of variation is not conducive for inferring phylogenetic relationships among sequences representing N. battus or addressing the putative history of introduction.     

Genetic Diversity of Moose (Alces alces L.) in Eurasia

Polymorphism of nucleotide sequence of D-loop fragment of the mitochondrial DNA was studied in 20 moose from several local populations on the territory of Eurasia. Three main haplotype variants of D-loop were detected by molecular phylogenetic method, which formed three clusters named European, Asian, and American. Intraspecies variation in the length of HVSI of D-loop of the mitochondrial DNA of moose was revealed. In the Far Eastern and Yakutian moose, haplotypes with a 75-bp deletion were found, which were most similar with haplotypes (also with the deletion), earlier observed in North American moose [1]. The highest diversity of the haplotypes of mitochondrial DNA is characteristic of Yakutia and the Far East (where three haplotype variants were found), which demonstrates the probable role of the region as the center of the species origin or as the region of ancient population mixture. The geographic region might be considered as a probable source of ancient moose migrations from Asia to America, basing on the data of distribution of mitochondrial haplotypes of D-loop and alleles of MhcAlal-DRB1. Divergence of nucleotide sequences of haplotypes with the 75-bp deletion (forming the American cluster on the phylogenetic tree) was the lowest (0.4%), which evidences respectively recent origin of the group of haplotypes. In Europe, only haplotypes of mitochondrial DNA referred to European variant were observed. Basing on analysis of variation of nucleotide sequences of D-loop, exon 4 of -Casein and exon 2 of MhcAlal-DRB1, we demonstrated that Eurasian moose studied belong to the unique species, which has probably passed through a bottle neck. The time of the origin of modern diversity of D-loop haplotypes of the species was estimated as 0.075–0.15 Myr ago.     

African-derived mitochondria in South American native cattle breeds (Bos taurus): evidence of a new taurine mitochondrial lineage

This article reports the nucleotide diversity within the control region of 42 mitochondrial chromosomes belonging to five South American native cattle breeds (Bos taurus). Analysis of these data in conjunction with B. taurus and B. indicus sequences from Africa, Europe, the Near East, India, and Japan allowed the recognition of eight new mitochondrial haplotypes and their relative positions in a phylogenetic network. The structure of genetic variation among different hypothetical groupings was tested through the molecular variance decomposition, which was best explained by haplotype group components. Haplotypes surveyed were classified as European-related and African-related. Unexpectedly, two haplotypes within the African cluster were more divergent from the African consensus than the latter from the European consensus. A neighbor-joining tree shows the position of two haplotypes compared to European/African mitochondrial lineage splitting. This different and putatively ancestral mitochondrial lineage (AA) is supported by the calibration of sequence divergence based on the Bos-Bison separation. The European/African mitochondria divergence might be subsequent (67,100 years before present) to that between AA and Africans (84,700 years before present), also preceding domestication times. These genetic data could reflect the haplotype distribution of Iberian cattle five centuries ago.     

Parallel evolution in radiation ofOhomopterus ground beetles inferred from mitochondrial ND5 gene sequences

Molecular phylogenetic analyses using mitochondrial NADH dehydrogenase subunit 5 (ND5) gene sequences representing all 15 species and the majority of subspecies or races of theOhomopterus ground beetles from all over the Japanese archipelago have uncovered a remarkable evolutionary history. Clustering of the species in the molecular phylogenetic tree is linked to their geographic distribution and does not correlate with morphological characters. Taxonomically the same species or the members belonging to the same species-group fall out in more than two different places on the ND5 tree. Evidence has been presented against a possible participation of ancestral polymorphism and random lineage sorting or of hybrid individuals for the observed distribution of mitochondrial DNA haplotypes. The most plausible explanation of our results is that parallel evolution took place in different lineages. Most notably,O. dehaanii, O. yaconinus, andO. japonicus in a lineage reveal almost identical morphology with those of the same species (or subspecies) but belonging to the phylogenetically remote lineages.The nucleotide sequence data reported in this paper will appear in the DDBJ, EMBL, and GenBank nucleotide sequence databases with accession numbers D50711-DD-50733 and D87131-D87186.     

Genetic study of the Paleolithic and Neolithic Southeast Asians.

DNA samples were extracted from six prehistoric human remains, found on the Malay Peninsula, dating to the Paleolithic and the Neolithic periods. Nucleotide sequences of mitochondrial DNA were determined by the polymerase chain reaction-direct sequencing method. A phylogenetic tree between prehistoric and present humans was constructed based on the nucleotide sequence data. Mitochondrial DNA phylogenetic relationships and ethnoarchaeological evidence suggest that there is a continuity beetween the pre-Neolithic humans and the present Semang and that the Neolithic humans in this area might be an ancestral group of the Senoi.     

Ancestral sequence reconstruction in primate mitochondrial DNA: compositional bias and effect on functional inference

Reconstruction of ancestral DNA and amino acid sequences is an important means of inferring information about past evolutionary events. Such reconstructions suggest changes in molecular function and evolutionary processes over the course of evolution and are used to infer adaptation and convergence. Maximum likelihood (ML) is generally thought to provide relatively accurate reconstructed sequences compared to parsimony, but both methods lead to the inference of multiple directional changes in nucleotide frequencies in primate mitochondrial DNA (mtDNA). To better understand this surprising result, as well as to better understand how parsimony and ML differ, we constructed a series of computationally simple "conditional pathway" methods that differed in the number of substitutions allowed per site along each branch, and we also evaluated the entire Bayesian posterior frequency distribution of reconstructed ancestral states. We analyzed primate mitochondrial cytochrome b (Cyt-b) and cytochrome oxidase subunit I (COI) genes and found that ML reconstructs ancestral frequencies that are often more different from tip sequences than are parsimony reconstructions. In contrast, frequency reconstructions based on the posterior ensemble more closely resemble extant nucleotide frequencies. Simulations indicate that these differences in ancestral sequence inference are probably due to deterministic bias caused by high uncertainty in the optimization-based ancestral reconstruction methods (parsimony, ML, Bayesian maximum a posteriori). In contrast, ancestral nucleotide frequencies based on an average of the Bayesian set of credible ancestral sequences are much less biased. The methods involving simpler conditional pathway calculations have slightly reduced likelihood values compared to full likelihood calculations, but they can provide fairly unbiased nucleotide reconstructions and may be useful in more complex phylogenetic analyses than considered here due to their speed and flexibility. To determine whether biased reconstructions using optimization methods might affect inferences of functional properties, ancestral primate mitochondrial tRNA sequences were inferred and helix-forming propensities for conserved pairs were evaluated in silico. For ambiguously reconstructed nucleotides at sites with high base composition variability, ancestral tRNA sequences from Bayesian analyses were more compatible with canonical base pairing than were those inferred by other methods. Thus, nucleotide bias in reconstructed sequences apparently can lead to serious bias and inaccuracies in functional predictions.     

A reanalysis of the ancient mitochondrial DNA sequences recovered from Neandertal bones

Recent reports analyzing mitochondrial DNA sequences from Neandertal bones have claimed that Neandertals and modern humans are different species. The phylogenetic analyses carried out in these articles did not take into account the high substitution rate variation among sites observed in the human mitochondrial D-loop region and also lack an estimation of the parameters of the nucleotide substitution model. The separate phylogenetic position of Neandertals is not supported when these factors are considered. Our analysis shows that Neandertal-Human and Human-Human pairwise distance distributions overlap more than what previous studies suggested. We also show that the most ancient Neandertal HVI region is the most divergent when compared with modern human sequences. However, the opposite would be expected if the sequence had not been modified since the death of the specimen. Such incongruence is discussed in the light of diagenetic modifications in ancient Neandertal DNA sequences.     

Estimating the pattern of nucleotide substitution

Knowledge of the pattern of nucleotide substitution is important both to our understanding of molecular sequence evolution and to reliable estimation of phylogenetic relationships. The method of parsimony analysis, which has been used to estimate substitution patterns in real sequences, has serious drawbacks and leads to results difficult to interpret. In this paper a model-based maximum likelihood approach is proposed for estimating substitution patterns in real sequences. Nucleotide substitution is assumed to follow a homogeneous Markov process, and the general reversible process model (REV) and the unrestricted model without the reversibility assumption are used. These models are also applied to examine the adequacy of the model of Hasegawa et al. (J. Mol. Evol. 1985;22:160–174) (HKY85). Two data sets are analyzed. For the -globin pseudogenes of six primate species, the REV model fits the data much better than HKY85, while, for a segment of mtDNA sequences from nine primates, REV cannot provide a significantly better fit than HKY85 when rate variation over sites is taken into account in the models. It is concluded that the use of the REV model in phylogenetic analysis can be recommended, especially for large data sets or for sequences with extreme substitution patterns, while HKY85 may be expected to provide a good approximation. The use of the unrestricted model does not appear to be worthwhile.     

PHYMYCO-DB: A Curated Database for Analyses of Fungal Diversity and Evolution

Background

In environmental sequencing studies, fungi can be identified based on nucleic acid sequences, using either highly variable sequences as species barcodes or conserved sequences containing a high-quality phylogenetic signal. For the latter, identification relies on phylogenetic analyses and the adoption of the phylogenetic species concept.Such analysis requires that the reference sequences are well identified and deposited in public-access databases. However, many entries in the public sequence databases are problematic in terms of quality and reliability and these data require screening to ensure correct phylogenetic interpretation.

Methods and Principal Findings

To facilitate phylogenetic inferences and phylogenetic assignment, we introduce a fungal sequence database. The database PHYMYCO-DB comprises fungal sequences from GenBank that have been filtered to satisfy stringent sequence quality criteria. For the first release, two widely used molecular taxonomic markers were chosen: the nuclear SSU rRNA and EF1-α gene sequences. Following the automatic extraction and filtration, a manual curation is performed to remove problematic sequences while preserving relevant sequences useful for phylogenetic studies. As a result of curation, ∼20% of the automatically filtered sequences have been removed from the database. To demonstrate how PHYMYCO-DB can be employed, we test a set of environmental Chytridiomycota sequences obtained from deep sea samples.

Conclusion

PHYMYCO-DB offers the tools necessary to: (i) extract high quality fungal sequences for each of the 5 fungal phyla, at all taxonomic levels, (ii) extract already performed alignments, to act as ‘reference alignments’, (iii) launch alignments of personal sequences along with stored data. A total of 9120 SSU rRNA and 672 EF1-α high-quality fungal sequences are now available.The PHYMYCO-DB is accessible through the URL http://phymycodb.genouest.org/.     

六种灵猫科动物的分子系统关系研究

对六种灵猫科物种线粒体12 S rRNA基因及其中四种的Cytb基因部分序列进行了测定,并从Gen-Bank获得斑林狸(Prionodon pardicolor)、熊狸(Arctictis binturong)的Cytb基因同源序列。两基因整合序列比对后长755 bp,12 S rRNA基因序列中有70个变异位点,31个简约信息位点,在Cytb基因序列中,共有120个位点呈现变异,60个简约信息位点,Cytb基因的碱基变异百分比高于12 S rRNA基因的碱基变异百分比。使用邻接法(NJ)、最大似然法(ML)重建的分子系统树显示:斑林狸从灵猫亚科中分离出来,支持灵猫亚科的多系起源,而且斑林狸可能是中国起源最早且最特化的灵猫科动物。另外,同属于灵猫亚科的大灵猫(Viverra zibe-tha)、小灵猫(Viverricula indica)聚为一支,同属于棕榈狸亚科的果子狸(Viverricula indica)、熊狸聚为姐妹群,这些与传统形态学分类观点一致。     

Phylogenetic Utility of the External Transcribed Spacer (ETS) of 18S–26S rDNA: Congruence of ETS and ITS Trees ofCalycadenia(Compositae)

The 3′ region of the external transcribed spacer (ETS) of 18S–26S nuclear ribosomal DNA was sequenced in 19 representatives ofCalycadenia/Osmadeniaand two outgroup species (Compositae) to assess its utility for phylogeny reconstruction compared to rDNA internal transcribed spacer (ITS) data. Universal primers based on plant, fungal, and animal sequences were designed to amplify the intergenic spacer (IGS) and an angiosperm primer was constructed to sequence the 3′ end of the ETS in members of tribe Heliantheae. Based on these sequences, an internal ETS primer useful across Heliantheaesensu latowas designed to amplify and sequence directly the 3′ ETS region in the study taxa, which were the subjects of an earlier phylogenetic investigation based on ITS sequences. Size variation in the amplified ETS region varied across taxa of Heliantheaesensu latofrom approximately 350 to 700 bp, in part attributable to an approximately 200-bp tandem duplication in a common ancestor ofCalycadenia/Osmadenia.Phylogenetic analysis of the 200-bp subrepeats and examination of apomorphic changes in the duplicated region demonstrate that the subrepeats inCalycadenia/Osmadeniahave evolved divergently. Phylogenetic analyses of the entire amplified ETS region yielded a highly resolved strict consensus tree that is nearly identical in topology to the ITS tree, with strong bootstrap and decay support on most branches. Parsimony analyses of combined ETS and ITS data yielded a strict consensus tree that is better resolved and generally better supported than trees based on either data set analyzed separately. We calculated an approximately 1.3- to 2.4-fold higher rate of sequence evolution by nucleotide substitution in the ETS region studied than in ITS-1 + ITS-2. A similar disparity in the proportion of variable (1.3 ETS:1 ITS) and potentially informative (1.5 ETS:1 ITS) sites was observed for the ingroup. Levels of homoplasy are similar in the ETS and ITS data. We conclude that the ETS holds great promise for augmenting ITS data for phylogenetic studies of young lineages.     

Complete mitochondrial genome of the <Emphasis Type="Italic">Pleuronichthys lighti</Emphasis> (Pleuronectiformes,Pleuronectidae) with phylogenetic consideration

Mitochondrial genome has been used to shed light on many fields of both basic and applied research, including the study of molecular evolution. The complete mitochondrial genome sequence of 17368 bp nucleotides from the Pleuronichthys lighti was determined. It was a circular double-stranded DNA molecule with identical set of 22 transfer RNA genes, 2 ribosomal RNA genes, 13 protein-coding genes as well as a non-coding control region. Stand asymmetry in the nucleotide composition was reflected in the codon usage of genes oriented in opposite directions. In the control region, we identified the extended termination associated sequence domain, the central conserved sequence block domain and the conserved sequence block domain, and two complete repeat region. They were “TTACAATA” and “TGTTGTAA”, respectively. All known 12 mitochondrial genomes of Pleuronectinae fishes were downloaded and analyzed; there were 5570 variable sites in the consensus sequences of 15241 base pairs, calculation of total sites were 35.5%. The highest sequence divergence was 50% (ATP8) and the Kimura-2-parameter genetic distance was 0.235 (ND6), whereas the COIII had the lowest sequence divergence (28.8%) and genetic distance (0.128); the protein coding genes were mainly acted by purifying selection which was detected by selection tests. Analysis of confidence and the information content for per nucleotide revealed ND5, ATP6, COI and ND4 genes were suitable molecular markers for phylogenetic study of Pleuronectinae fishes. Phylogenetic analysis using Bayesian computational algorithms based on COI genes provided support for the taxonomic status of P. lighti, which was consistent with the traditional taxonomy.     

Monophyly of clade III nematodes is not supported by phylogenetic analysis of complete mitochondrial genome sequences

Background

The orders Ascaridida, Oxyurida, and Spirurida represent major components of zooparasitic nematode diversity, including many species of veterinary and medical importance. Phylum-wide nematode phylogenetic hypotheses have mainly been based on nuclear rDNA sequences, but more recently complete mitochondrial (mtDNA) gene sequences have provided another source of molecular information to evaluate relationships. Although there is much agreement between nuclear rDNA and mtDNA phylogenies, relationships among certain major clades are different. In this study we report that mtDNA sequences do not support the monophyly of Ascaridida, Oxyurida and Spirurida (clade III) in contrast to results for nuclear rDNA. Results from mtDNA genomes show promise as an additional independently evolving genome for developing phylogenetic hypotheses for nematodes, although substantially increased taxon sampling is needed for enhanced comparative value with nuclear rDNA. Ultimately, topological incongruence (and congruence) between nuclear rDNA and mtDNA phylogenetic hypotheses will need to be tested relative to additional independent loci that provide appropriate levels of resolution.

Results

For this comparative phylogenetic study, we determined the complete mitochondrial genome sequences of three nematode species, Cucullanus robustus (13,972 bp) representing Ascaridida, Wellcomia siamensis (14,128 bp) representing Oxyurida, and Heliconema longissimum (13,610 bp) representing Spirurida. These new sequences were used along with 33 published nematode mitochondrial genomes to investigate phylogenetic relationships among chromadorean orders. Phylogenetic analyses of both nucleotide and amino acid sequence datasets support the hypothesis that Ascaridida is nested within Rhabditida. The position of Oxyurida within Chromadorea varies among analyses; in most analyses this order is sister to the Ascaridida plus Rhabditida clade, with representative Spirurida forming a distinct clade, however, in one case Oxyurida is sister to Spirurida. Ascaridida, Oxyurida, and Spirurida (the sampled clade III taxa) do not form a monophyletic group based on complete mitochondrial DNA sequences. Tree topology tests revealed that constraining clade III taxa to be monophyletic, given the mtDNA datasets analyzed, was a significantly worse result.

Conclusion

The phylogenetic hypotheses from comparative analysis of the complete mitochondrial genome data (analysis of nucleotide and amino acid datasets, and nucleotide data excluding 3^rd positions) indicates that nematodes representing Ascaridida, Oxyurida and Spirurida do not share an exclusive most recent common ancestor, in contrast to published results based on nuclear ribosomal DNA. Overall, mtDNA genome data provides reliable support for nematode relationships that often corroborates findings based on nuclear rDNA. It is anticipated that additional taxonomic sampling will provide a wealth of information on mitochondrial genome evolution and sequence data for developing phylogenetic hypotheses for the phylum Nematoda.

     

Corrections to the human mitochondrial ribosomal RNA sequences

Two new errors and one consensus change were identified in the human mitochondrial Cambridge consensus sequence. The errors are an A to G substitution at nucleotide 750 in the 12S rRNA gene and a single nucleotide deletion at nt 3107 in the 16S rRNA gene. The consensus change is nt 2706 AG in the 16S rRNA gene.     

Identification of astigmatid mites using the second internal transcribed spacer (ITS2) region and its application for phylogenetic study

The second internal transcribed spacer (ITS2) of nuclear ribosomal DNA from 73 specimens of Astigmata was analyzed by PCR amplification and DNA sequencing. The length of the ITS2 region varied from 282 to 592 bp. The interspecific variation based on consensus sequences was more than 4.1%, while the intraspecific or intra-individual variation was from 0 to 5.7%. The variation between geographically separated populations (0–3.2%) was almost the same as the variation within strains. The sequences of the ITS2 region of Astigmata were concluded to be species-specific. The phylogenetic tree inferred from the ITS2 region supported Zachvatkins morphological classification in the subfamily Rhizoglyphinae. The species-specific ITS2 sequence is useful for the species identification of astigmatid mites and for studying low-level phylogenetic relationships.Chemical Ecology of Astigmatid Mites LXXVThis revised version was published online in May 2005 with a corrected cover date.     

FASconCAT-G: extensive functions for multiple sequence alignment preparations concerning phylogenetic studies

Background

Phylogenetic and population genetic studies often deal with multiple sequence alignments that require manipulation or processing steps such as sequence concatenation, sequence renaming, sequence translation or consensus sequence generation. In recent years phylogenetic data sets have expanded from single genes to genome wide markers comprising hundreds to thousands of loci. Processing of these large phylogenomic data sets is impracticable without using automated process pipelines. Currently no stand-alone or pipeline compatible program exists that offers a broad range of manipulation and processing steps for multiple sequence alignments in a single process run.

Results

Here we present FASconCAT-G, a system independent editor, which offers various processing options for multiple sequence alignments. The software provides a wide range of possibilities to edit and concatenate multiple nucleotide, amino acid, and structure sequence alignment files for phylogenetic and population genetic purposes. The main options include sequence renaming, file format conversion, sequence translation between nucleotide and amino acid states, consensus generation of specific sequence blocks, sequence concatenation, model selection of amino acid replacement with ProtTest, two types of RY coding as well as site exclusions and extraction of parsimony informative sites. Convieniently, most options can be invoked in combination and performed during a single process run. Additionally, FASconCAT-G prints useful information regarding alignment characteristics and editing processes such as base compositions of single in- and outfiles, sequence areas in a concatenated supermatrix, as well as paired stem and loop regions in secondary structure sequence strings.

Conclusions

FASconCAT-G is a command-line driven Perl program that delivers computationally fast and user-friendly processing of multiple sequence alignments for phylogenetic and population genetic applications and is well suited for incorporation into analysis pipelines.

     

Base sequence studies of 300 nucleotide renatured repeated human DNA clones

A band of 300 nucleotide long duplex DNA is released by treating renatured repeated human DNA with the single strand-specific endonuclease S₁. Since many of the interspersed repeated sequences in human DNA are 300 nucleotides long, this band should be enriched in such repeats. We have determined the nucleotide sequences of 15 clones constructed from these 300 nucleotide S₁-resistant repeats. Ten of these cloned sequences are members of the Alu family of interspersed repeats. These ten sequences share a recognizable consensus sequence from which individual clones have an average divergence of 12.8%. The 300 nucleotide Alu family consensus sequence has a dimeric structure and was evidently formed from a head to tail duplication of an ancestral monomeric sequence. Three of the remaining clones are variations on a simple pentanucleotide sequence previously reported for human satellite III DNA. Two of the 15 clones have distinct and complex sequences and may represent other families of interspersed repeated sequences.     

A Complete Mitochondrial Genome Sequence from a Mesolithic Wild Aurochs (Bos primigenius)

Background

The derivation of domestic cattle from the extinct wild aurochs (Bos primigenius) has been well-documented by archaeological and genetic studies. Genetic studies point towards the Neolithic Near East as the centre of origin for Bos taurus, with some lines of evidence suggesting possible, albeit rare, genetic contributions from locally domesticated wild aurochsen across Eurasia. Inferences from these investigations have been based largely on the analysis of partial mitochondrial DNA sequences generated from modern animals, with limited sequence data from ancient aurochsen samples. Recent developments in DNA sequencing technologies, however, are affording new opportunities for the examination of genetic material retrieved from extinct species, providing new insight into their evolutionary history. Here we present DNA sequence analysis of the first complete mitochondrial genome (16,338 base pairs) from an archaeologically-verified and exceptionally-well preserved aurochs bone sample.

Methodology

DNA extracts were generated from an aurochs humerus bone sample recovered from a cave site located in Derbyshire, England and radiocarbon-dated to 6,738±68 calibrated years before present. These extracts were prepared for both Sanger and next generation DNA sequencing technologies (Illumina Genome Analyzer). In total, 289.9 megabases (22.48%) of the post-filtered DNA sequences generated using the Illumina Genome Analyzer from this sample mapped with confidence to the bovine genome. A consensus B. primigenius mitochondrial genome sequence was constructed and was analysed alongside all available complete bovine mitochondrial genome sequences.

Conclusions

For all nucleotide positions where both Sanger and Illumina Genome Analyzer sequencing methods gave high-confidence calls, no discrepancies were observed. Sequence analysis reveals evidence of heteroplasmy in this sample and places this mitochondrial genome sequence securely within a previously identified aurochsen haplogroup (haplogroup P), thus providing novel insights into pre-domestic patterns of variation. The high proportion of authentic, endogenous aurochs DNA preserved in this sample bodes well for future efforts to determine the complete genome sequence of a wild ancestor of domestic cattle.     

Estimating evolution of temporal sequence changes: a practical approach to inferring ancestral developmental sequences and sequence heterochrony

Developmental biology often yields data in a temporal context. Temporal data in phylogenetic systematics has important uses in the field of evolutionary developmental biology and, in general, comparative biology. The evolution of temporal sequences, specifically developmental sequences, has proven difficult to examine due to the highly variable temporal progression of development. Issues concerning the analysis of temporal sequences and problems with current methods of analysis are discussed. We present here an algorithm to infer ancestral temporal sequences, quantify sequence heterochronies, and estimate pseudoreplicate consensus support for sequence changes using Parsimov-based genetic inference [PGi]. Real temporal developmental sequence data sets are used to compare PGi with currently used approaches, and PGi is shown to be the most efficient, accurate, and practical method to examine biological data and infer ancestral states on a phylogeny. The method is also expandable to address further issues in developmental evolution, namely modularity.     

Nonneutral Mitochondrial DNA Variation in Humans and Chimpanzees

We sequenced the NADH dehydrogenase subunit 3 (ND3) gene from a sample of 61 humans, five common chimpanzees, and one gorilla to test whether patterns of mitochondrial DNA (mtDNA) variation are consistent with a neutral model of molecular evolution. Within humans and within chimpanzees, the ratio of replacement to silent nucleotide substitutions was higher than observed in comparisons between species, contrary to neutral expectations. To test the generality of this result, we reanalyzed published human RFLP data from the entire mitochondrial genome. Gains of restriction sites relative to a known human mtDNA sequence were used to infer unambiguous nucleotide substitutions. We also compared the complete mtDNA sequences of three humans. Both the RFLP data and the sequence data reveal a higher ratio of replacement to silent nucleotide substitutions within humans than is seen between species. This pattern is observed at most or all human mitochondrial genes and is inconsistent with a strictly neutral model. These data suggest that many mitochondrial protein polymorphisms are slightly deleterious, consistent with studies of human mitochondrial diseases.