The complete nucleotide sequence of the mitochondrial DNA of the agnathan Lampetra fluviatilis: bearings on the phylogeny of cyclostomes

There are two competing theories about the interrelationships of craniates: the cyclostome theory assumes that lampreys and hagfishes are a clade, the cyclostomes, whose sister group is the jawed vertebrates (gnathostomes); the vertebrate theory assumes that lampreys and gnathostomes are a clade, the vertebrates, whose sister group is hagfishes. The vertebrate theory is best supported by a number of unique anatomical and physiological characters. Molecular sequence data from 18S and 28S rRNA genes rather support the cyclostome theory, but mtDNA sequence of Myxine glutinosa rather supports the vertebrate theory. Additional molecular data are thus needed to elucidate this three-taxon problem. We determined the complete nucleotide sequence of the mtDNA of the lamprey Lampetra fluviatilis. The mtDNA of L. fluviatilis possesses the same genomic organization as Petromyzon marinus, which validates this gene order as a synapomorphy of lampreys. The mtDNA sequence of L. fluviatilis was used in combination with relevant mtDNA sequences for an approach to the hagfish/lamprey relationships using the maximum-parsimony, neighbor-joining, and maximum-likelihood methods. Although trees compatible with our present knowledge of the phylogeny of craniates can be reconstructed by using the three methods, the data collected do not support the vertebrate or the cyclostome hypothesis. The present data set does not allow the resolution of this three-taxon problem, and new kinds of data, such as nuclear DNA sequences, need to be collected.     

Monophyly of Lampreys and Hagfishes Supported by Nuclear DNA–Coded Genes

The phylogenetic position of hagfishes in vertebrate evolution is currently controversial. The 18S and 28S rRNA trees support the monophyly of hagfishes and lampreys. In contrast, the mitochondrial DNAs suggest the close association of lampreys and gnathostomes. To clarify this controversial issue, we have conducted cloning and sequencing of the four nuclear DNA–coded single-copy genes encoding the triose phosphate isomerase, calreticulin, and the largest subunit of RNA polymerase II and III. Based on these proteins, together with the Mn superoxide dismutase for which hagfish and lamprey sequences are available in database, phylogenetic trees have been inferred by the maximum likelihood (ML) method of protein phylogeny. It was shown that all the five proteins prefer the monophyletic tree of cyclostomes, and the total log-likelihood of the five proteins significantly supports the cyclostome monophyly at the level of ±1 SE. The ML trees of aldolase family comprising three nonallelic isoforms and the complement component group comprising C3, C4, and C5, both of which diverged during vertebrate evolution by gene duplications, also suggest the cyclostome monophyly. Received: 28 April 1999 / Accepted: 30 June 1999     

The mitochondrial DNA molecule of the hagfish (myxine glutinosa) and vertebrate phylogeny

The vertebrates are traditionally classified into two distinct groups, Agnatha (jawless vertebrates) and Gnathostomata (jawed vertebrates). Extant agnathans are represented by hagfishes (Myxiniformes) and lampreys (Petromyzontiformes), frequently grouped together within the Cyclostomata. Whereas the recognition of the Gnathostomata as a clade is commonly acknowledged, a consensus has not been reached regarding whether or not Cyclostomata represents a clade. In the present study we have used newly established sequences of the protein-coding genes of the mitochondrial DNA molecule of the hagfish to explore agnathan and gnathostome relationships. The phylogenetic analysis of Pisces, using echinoderms as outgroup, placed the hagfish as a sister group of Vertebrata sensu stricto, i.e., the lamprey and the gnathostomes. The phylogenetic analysis of the Gnathostomata identified a basal divergence between gnathostome fishes and a branch leading to birds and mammals, i.e., between ``Anamnia' and Amniota. The lungfish has a basal position among gnathostome fishes with the teleosts as the most recently evolving lineage. The findings portray a hitherto unrecognized polarity in the evolution of bony fishes. The presently established relationships are incompatible with previous molecular studies. Received: 15 August 1997 / Accepted: 1 October 1997     

Hagfish (cyclostomata, vertebrata): searching for the ancestral developmental plan of vertebrates

The phylogenetic position of the hagfish remains enigmatic. In contrast to molecular data that suggest monophyly of the cyclostomes, several morphological features imply a more ancestral state of this animal compared with the lampreys. To resolve this question requires an understanding of the embryology of the hagfish, especially of the neural crest. The early development of the hagfish has long remained a mystery. We collected a shallow-water-dwelling hagfish, Eptatretus burgeri, set up an aquarium tank designed to resemble its habitat, and successfully obtained several embryos. By observing the histology and expression of genes known to play fundamental roles in the neural crest, we found that the hagfish crest develops as delaminating migratory cells, as in other vertebrates. We conclude that the delaminating neural crest is a vertebrate synapomorphy that seems to have appeared from the beginning of their evolutionary history, before the splitting away of the hagfish lineage.     

Molecular phylogeny of early vertebrates: monophyly of the agnathans as revealed by sequences of 35 genes

Extant vertebrates are divided into three major groups: hagfishes (Hyperotreti, myxinoids), lampreys (Hyperoartia, petromyzontids), and jawed vertebrates (Gnathostomata). The phylogenetic relationships among the groups and within the jawed vertebrates are controversial, for both morphological and molecular studies have rendered themselves to conflicting interpretations. Here, we use the sequences of 35 nuclear protein-encoding genes to provide definitive evidence for the monophyly of the Agnatha (jawless vertebrates, a group encompassing the hagfishes and lampreys). Our analyses also give a strong support for the separation of Chondrichthyes (cartilaginous fishes) before the divergence of Osteichthyes (bony fishes) from the other gnathostomes.     

Time scale for cyclostome evolution inferred with a phylogenetic diagnosis of hagfish and lamprey cDNA sequences

The Cyclostomata consists of the two orders Myxiniformes (hagfishes) and Petromyzoniformes (lampreys), and its monophyly has been unequivocally supported by recent molecular phylogenetic studies. Under this updated vertebrate phylogeny, we performed in silico evolutionary analyses using currently available cDNA sequences of cyclostomes. We first calculated the GC-content at four-fold degenerate sites (GC(4)), which revealed that an extremely high GC-content is shared by all the lamprey species we surveyed, whereas no striking pattern in GC-content was observed in any of the hagfish species surveyed. We then estimated the timing of diversification in cyclostome evolution using nucleotide and amino acid sequences. We obtained divergence times of 470-390 million years ago (Mya) in the Ordovician-Silurian-Devonian Periods for the interordinal split between Myxiniformes and Petromyzoniformes; 90-60 Mya in the Cretaceous-Tertiary Periods for the split between the two hagfish subfamilies, Myxininae and Eptatretinae; 280-220 Mya in the Permian-Triassic Periods for the split between the two lamprey subfamilies, Geotriinae and Petromyzoninae; and 30-10 Mya in the Tertiary Period for the split between the two lamprey genera, Petromyzon and Lethenteron. This evolutionary configuration indicates that Myxiniformes and Petromyzoniformes diverged shortly after the common ancestor of cyclostomes split from the future gnathostome lineage. Our results also suggest that intra-subfamilial diversification in hagfish and lamprey lineages (especially those distributed in the northern hemisphere) occurred in the Cretaceous or Tertiary Periods.     

Bayesian phylogenetic analysis supports monophyly of ambulacraria and of cyclostomes

Vertebrates are part of the phylum Chordata, itself part of a three-phylum group known as the deuterostomes. Despite extensive phylogenetic analysis of the deuterostome animals, several unresolved relationships remain. These include the relationship between the three deuterostome phyla (chordates, echinoderms and hemichordates), and the monophyletic or paraphyletic origin of the cyclostomes (hagfish and lampreys). Using robust Bayesian statistical analysis of 18S ribosomal DNA, mitochondrial genes and nuclear protein-coding DNA, we find strong support for a hemichordate-echinoderm clade, and for monophyly of the cyclostomes.     

Phylogenetic analysis of 48 gene families revealing relationships between Hagfishes, Lampreys, and Gnathostomata

It has become clear that the extant vertebrates are divided into three major groups, that is, hagfishes, lampreys, and jawed vertebrates.Morphological and molecular studies, however, have resulted in conflicting views with regard m their interrelationships. To clarify the phylogenetic relationships between them, 48 orthologous protein-coding gene families were analyzed. Even as the analysis of 34 nuclear gene families supported the monophyly of cyclostomes, the analysis of 14 mitochondrial gene families suggested a closer relationship between lampreys and gnathostomes compared to hagfishes. Lampreys were sister group of gnathostomes. The results of this study sup-ported the eyclostomes. Choice of outgroup, tree-making methods, and software may affect the phylogenetic prediction, which may have caused much debate over the subject. Development of new methods for tackling such problems is still necessary.     

The origin of developmental mechanisms underlying vertebral elements: implications from hagfish evo-devo

The origins of the vertebral elements and the underlying developmental mechanisms have so far remained unclear, largely due to the unusual axial skeletal morphology of hagfish, one of two extant jawless vertebrate clades. Hagfish axial supporting tissue is generally believed to consist of the notochord and cartilaginous fin rays only. However, careful investigations of whether vertebral elements are truly absent in hagfish are scarce, and it is also unclear whether the axial skeletal morphology of the hagfish is an ancestral or a derived condition. To address these questions, we re-examined the axial skeletal morphology of the Japanese inshore hagfish (Eptatretus burgeri). Based on a report published a century ago which implied the existence of vertebral elements in hagfish, we conducted anatomical and histological analyses of the hagfish axial skeletal systems and their development. Through this analysis, we demonstrate that hagfish possesses sclerotome-derived cartilaginous vertebral elements at the ventral aspect of the notochord. Based on (i) molecular phylogenetic evidence in support of the monophyly of cyclostomes (hagfish and lampreys) and jawed vertebrates (gnathostomes), and (ii) the morphology of the vertebral elements in extant gnathostomes and cyclostomes, we propose that the embryos of the common ancestor of all vertebrates would have possessed sclerotomal cells that formed the segmentally arranged vertebral elements attached to the notochord. We also conclude that the underlying developmental mechanisms are likely to have been conserved among extinct jawless vertebrates and modern gnathostomes.     

Phylogenetic patterns in populations of Chilean species of the genus Orestias (Teleostei: Cyprinodontidae): results of mitochondrial DNA analysis

Patterns of molecular genetic differentiation among taxa of the "agassii species complex" (Parenti, 1984) were analysed based on partial mtDNA control region sequences. Special attention has been paid to Chilean populations of Orestias agassii and species from isolated lakes of northern Chile, e.g., O. agassii, Orestias chungarensis, Orestias parinacotensis, Orestias laucaensis, and Orestias ascotanensis. Orestias tschudii, Orestias luteus, and Orestias ispi were analysed comparatively. Our findings support the utility of mtDNA control region sequences for phylogenetic studies within the "agassii species complex" and confirmed the monophyly of this particular lineage, excluding O. luteus. However, the monophyly of further morphologically defined lineages within the "agassii complex" appears doubtful. No support was found for the utility of these data sets for inferring phylogenetic relationships between more distantly related taxa originating from Lake Titicaca.     

Complete mitochondrial DNA replacement in a Lake Tanganyika cichlid fish

We used nuclear and mitochondrial DNA (mtDNA) sequences from specimens collected throughout Lake Tanganyika to clarify the evolutionary relationship between Lamprologus callipterus and Neolamprologus fasciatus . The nuclear data support the reciprocal monophyly of these two shell-breeding lamprologine cichlids. However, mtDNA sequences show that (i) L. callipterus includes two divergent and geographically disjunct (North–South) mtDNA lineages; and that (ii) N. fasciatus individuals cluster in a lineage sister group to the northern lineage of L. callipterus . The two mtDNA lineages of L. callipterus diverged c . 684 kya to 1.2 Ma, coinciding with a major water level low stand in Lake Tanganyika, which divided the lake into isolated sub-lakes. This suggests that the two mtDNA lineages originated as the result of the separation of L. callipterus populations in different sub-basins. The incongruent phylogenetic position of N. fasciatus can best be explained by an ancient unidirectional introgression from L. callipterus into N. fasciatus. Remarkably, our data indicate that this event resulted in the complete mtDNA replacement in N. fasciatus . Our data suggest that hybridization occurred soon after the divergence of the two L. callipterus mtDNA lineages, probably still during the water level low stand, and that subsequently the invading mtDNA lineage spread throughout the lake.     

Resolving deep phylogenetic relationships in salamanders: analyses of mitochondrial and nuclear genomic data

Phylogenetic relationships among salamander families illustrate analytical challenges inherent to inferring phylogenies in which terminal branches are temporally very long relative to internal branches. We present new mitochondrial DNA sequences, approximately 2,100 base pairs from the genes encoding ND1, ND2, COI, and the intervening tRNA genes for 34 species representing all 10 salamander families, to examine these relationships. Parsimony analysis of these mtDNA sequences supports monophyly of all families except Proteidae, but yields a tree largely unresolved with respect to interfamilial relationships and the phylogenetic positions of the proteid genera Necturus and Proteus. In contrast, Bayesian and maximum-likelihood analyses of the mtDNA data produce a topology concordant with phylogenetic results from nuclear-encoded rRNA sequences, and they statistically reject monophyly of the internally fertilizing salamanders, suborder Salamandroidea. Phylogenetic simulations based on our mitochondrial DNA sequences reveal that Bayesian analyses outperform parsimony in reconstructing short branches located deep in the phylogenetic history of a taxon. However, phylogenetic conflicts between our results and a recent analysis of nuclear RAG-1 gene sequences suggest that statistical rejection of a monophyletic Salamandroidea by Bayesian analyses of our mitochondrial genomic data is probably erroneous. Bayesian and likelihood-based analyses may overestimate phylogenetic precision when estimating short branches located deep in a phylogeny from data showing substitutional saturation; an analysis of nucleotide substitutions indicates that these methods may be overly sensitive to a relatively small number of sites that show substitutions judged uncommon by the favored evolutionary model.     

Stochasticity overrules the "three-times rule": genetic drift,genetic draft,and coalescence times for nuclear loci versus mitochondrial DNA

Abstract.— Palumbi et al. (2001) proposed a "three-times rule" that uses mitochondrial DNA (mtDNA) sequences to predict probabilities of monophyly for nuclear loci (i.e., whether the alleles within a taxon coalesce with one another before they coalesce with alleles from a sister taxon). They use neutral coalescent theory to infer these probabilities from the ratio of interspecific divergence to intraspecific variation of mtDNA. We show that the estimated probabilities have very wide confidence intervals because of the inherent stochasticity of the mtDNA coalescent process. Under neutrality, the true probability of monophyly can be much higher, or much lower, than predicted by the three-times rule. We also review recent empirical and theoretical studies that refute neutrality-based predictions concerning mtDNA variation and divergence. We conclude that the three-times rule is neither a useful test for neutral molecular evolution nor a reliable guide to genealogical species.     

Phylogeny of Anophelinae (Diptera: Culicidae) based on nuclear ribosomal and mitochondrial DNA sequences

Abstract Phylogenetic relationships among thirty-two species of mosquitoes in subfamily Anophelinae are inferred from portions of the mitochondrial genes COI and COII, the nuclear 18S small subunit rRNA gene and the expansion D2 region of the nuclear large subunit 28S rRNA gene. Sequences were obtained from the genera Anopheles , Bironella and Chagasia . Representatives of all six subgenera of Anopheles were included: Anopheles , Cellia , Kerteszia , Lophopodomyia , Nyssorhynchus and Stethomyia. Using parsimony and maximum likelihood methods, various combinations of these DNA sequence data were analysed separately: 18S, 28S, combined 18S and 28S, combined COI and COII, and combined 18S, 28S, COI and COII ('total evidence'). The combined rDNA data contain strong phylogenetic signal, moderately to strongly supporting most clades in MP and ML analyses; however, the mtDNA data (analysed as either nucleotide or amino acid sequences) contain little phylogenetic signal, except for relationships of very recently derived groups of species and, at the deepest level, for the monophyly of Anophelinae. The paraphyly of Anopheles relative to Bironella is confirmed by most analyses and statistical tests. Support for the monophyly of subgenera Anopheles , Cellia , Kerteszia and Nyssorhynchus is indicated by most analyses. Subgenus Lophopodomyia is reconstructed as the sister to Bironella , nested within a clade also containing Nyssorhynchus and Kerteszia . The most basal relationships within genus Anopheles are not well resolved by any of the data partitions, although the results of statistical analyses of the rDNA data (S-H-tests, likelihood ratio tests for monophyly and Bayesian MCMC analyses) suggest that the clade consisting of Bironella , Lophopodomyia , Nyssorhynchus and Kerteszia is the sister to the clade containing Cellia and Anopheles .     

Multiple nuclear-gene phylogenies: application to pinnipeds and comparison with a mitochondrial DNA gene phylogeny

Phylogenetic analyses of closely related species should use information from multiple, independent genes with relatively high rates of sequence evolution. To investigate species for which there are few prior sequence data for single-copy nuclear (scnDNA) genes, primers for gene amplification can be designed to highly conserved regions of exons in order to amplify both coding (exons) and noncoding (introns) sequences. We have explored this approach in a phylogenetic analysis of six species of pinnipeds that, together with terrestrial carnivore outgroups, encompass divergence times < or = 40-50 Mya. We sequenced one intron from each of the aldolase A (ALD-A), aldolase C (ALD-C), and histone H2AF genes; one exon from the major-histocompatibility-complex DQA gene; a H2AF processed pseudogene (psi H2AF); and, for comparison with the nuclear genes, the 5' portion of the mitochondrial DNA (mtDNA) control region. The pinniped psi H2AF genes were found to be of limited use because they were paralogous with the gene in the outgroup. The rate of silent substitution in scnDNA (primarily introns) was 5-10-fold lower than that for mtDNA control region I, and scnDNA sequence divergence increased linearly with time < or = 40-50 Mya. Alleles at three polymorphic scnDNA loci (ALD-A, H2AF, and DQA) in the southern elephant seal were paraphyletic with respect to the allele from the closely related northern elephant seal, while the more numerous mtDNA alleles were monophyletic. This we attribute to the consequences of a higher mutation rate rather than to a lower effective population size of mtDNA compared with scnDNA. Within the short (i.e., < 500-bp) sequences of individual scnDNA sequences, phylogenetically informative variation was insufficient to obtain robust phylogenies. However, the combined scnDNA sequences produced a well-supported phylogeny congruent with that derived from mtDNA. This analysis illustrates the high resolution of mtDNA sequences compared with a similar length of scnDNA sequence, but it also demonstrates the utility of combining information from multiple short scnDNA sequences obtained using broadly applicable primers.      

Phylogenetic relationships, character evolution, and biogeography of the subfamily Lygosominae (Reptilia: Scincidae) inferred from mitochondrial DNA sequences

Phylogenetic relationships among the lygosomine skinks were inferred from 1249 base positions of mitochondrial DNA sequences of 12S and 16S rRNA genes. The monophyly of this subfamily was confirmed and the presence of five distinct infrasubfamilial lineages detected. Of these, the Sphenomorphus group appears to have diverged first, followed by the Lygosoma and Egernia groups in order, leaving the Eugongylus and Mabuya groups as sister groups. Our results did not support monophyly of the Mabuya group sensu lato (i.e., an assemblage of the Lygosoma, Egernia, and Mabuya groups), for which a number of morphological and karyological studies demonstrated a considerable similarity. Our results also contradict the previous hypothesis, formulated on the basis of morphological and immunological data, which argued for the sister relationship between the Egernia and the Eugongylus groups. Morphological and karyological characters used to define the Mabuya group (sensu lato) may actually represent plesiomorphic states. The phylogenetic diversity of lygosomine skinks in the Australian region appears to have increased through multiple colonizations from Southeast Asia.     

Evolutionary implications of lactate dehydrogenases (LDHs) of hagfishes compared to lampreys: LDH cDNA sequences from Eptatretus burgeri, Paramyxine atami and Eptatretus okinoseanus

Heart muscles of hagfishes Paramyxine atami and Eptatretus okinoseanus express the B4 isozyme of lactate dehydrogenase [L-LDH: NAD oxidoreductase, EC 1.1.1.27] (LDH-B4) whereas their skeletal muscles express LDH-A4. To examine the relationship of hagfish LDHs to lamprey and other vertebrate LDHs, we determined the cDNA sequences of LDH-A from three hagfishes and compared them with previously published sequences. A phylogenic tree shows that hagfishes diverged just after lampreys. The deduced amino acid sequences showed ten regions common to all vertebrate LDHs examined, i.e., the active site, the pocket recognizing the substrate-coenzyme complex, part of a loop at the surface, and the substrate binding site. The cyclostomate-specific regions (S1, S2) were located in the neighborhood of the active site loop. Three regions, IGS1, IGS2 and IGS3, seem to have altered their structures during the differentiation of LDH isozymes, and the regions remain in LDH-B of vertebrates hitherto examined. IGS2 and IGS3, which are in the neighborhood of the active site, may regulate catalytic activity. There were differences in six amino acid residues (6, 10, 20, 156, 269, and 341) in LDHs of hagfishes. These differences might reflect the tolerance to high pressure and low temperature of LDHs from hagfishes at different habitat depths.     

Phylogenetic analysis of four nuclear protein-encoding genes largely corroborates the traditional classification of Bivalvia (Mollusca)

Revived interest in molluscan phylogeny has resulted in a torrent of molecular sequence data from phylogenetic, mitogenomic, and phylogenomic studies. Despite recent progress, basal relationships of the class Bivalvia remain contentious, owing to conflicting morphological and molecular hypotheses. Marked incongruity of phylogenetic signal in datasets heavily represented by nuclear ribosomal genes versus mitochondrial genes has also impeded consensus on the type of molecular data best suited for investigating bivalve relationships. To arbitrate conflicting phylogenetic hypotheses, we evaluated the utility of four nuclear protein-encoding genes-ATP synthase β, elongation factor-1α, myosin heavy chain type II, and RNA polymerase II-for resolving the basal relationships of Bivalvia. We sampled all five major lineages of bivalves (Archiheterodonta, Euheterodonta [including Anomalodesmata], Palaeoheterodonta, Protobranchia, and Pteriomorphia) and inferred relationships using maximum likelihood and Bayesian approaches. To investigate the robustness of the phylogenetic signal embedded in the data, we implemented additional datasets wherein length variability and/or third codon positions were eliminated. Results obtained include (a) the clade (Nuculanida+Opponobranchia), i.e., the traditionally defined Protobranchia; (b) the monophyly of Pteriomorphia; (c) the clade (Archiheterodonta+Palaeoheterodonta); (d) the monophyly of the traditionally defined Euheterodonta (including Anomalodesmata); and (e) the monophyly of Heteroconchia, i.e., (Palaeoheterodonta+Archiheterodonta+Euheterodonta). The stability of the basal tree topology to dataset manipulation is indicative of signal robustness in these four genes. The inferred tree topology corresponds closely to those obtained by datasets dominated by nuclear ribosomal genes (18S rRNA and 28S rRNA), controverting recent taxonomic actions based solely upon mitochondrial gene phylogenies.     

Molecular phylogeny of a circum-global, diverse gastropod superfamily (Cerithioidea: Mollusca: Caenogastropoda): pushing the deepest phylogenetic limits of mitochondrial LSU rDNA sequences

The Cerithioidea is a very diverse group of gastropods with ca. 14 extant families and more than 200 genera occupying, and often dominating, marine, estuarine, and freshwater habitats. While the composition of Cerithioidea is now better understood due to recent anatomical and ultrastructural studies, the phylogenetic relationships among families remain chaotic. Morphology-based studies have provided conflicting views of relationships among families. We generated a phylogeny of cerithioideans based on mitochondrial large subunit rRNA and flanking tRNA gene sequences (total aligned data set 1873 bp). Nucleotide evidence and the presence of a unique pair of tRNA genes (i.e., threonine + glycine) between valine-mtLSU and the mtSSU rRNA gene support conclusions based on ultrastructural data that Vermetidae and Campanilidae are not Cerithioidea, certain anatomical similarities being due to convergent evolution. The molecular phylogeny shows support for the monophyly of the marine families Cerithiidae [corrected], Turritellidae, Batillariidae, Potamididae, and Scaliolidae as currently recognized. The phylogenetic data reveal that freshwater taxa evolved on three separate occasions; however, all three recognized freshwater families (Pleuroceridae, Melanopsidae, and Thiaridae) are polyphyletic. Mitochondrial rDNA sequences provide valuable data for testing the monophyly of cerithioidean [corrected] families and relationships within families, but fail to provide strong evidence for resolving relationships among families. It appears that the deepest phylogenetic limits for resolving caenogastropod relationships is less than about 245--241 mya, based on estimates of divergence derived from the fossil record.