Mitochondrial evidence on the phylogenetic position of caecilians (Amphibia: Gymnophiona)

Unexpected Discovery and Expression of Amphibian Class II Endogenous Retroviruses

Endogenous retroviruses (ERVs) are the remnants of past retroviral infections. Fossil records of class II retroviruses have been discovered in a range of vertebrates, with the exception of amphibians, which are known only to possess class I and class III-like ERVs. Through genomic mining of all available amphibian genomes, we identified, for the first time, class II ERVs in amphibians. The class II ERVs were found only in Gymnophiona (caecilians) and not in the genomes of the other amphibian orders, Anura (frogs and toads) and Caudata (salamanders and newts), which are phylogenetically closely related. Therefore, the ERV endogenization occurred after the split of Gymnophiona, Anura, and Caudata (323 million years ago). Investigation of phylogenetic relationship and genomic structure revealed that the ERVs may originate from alpha- or betaretroviruses. We offer evidence that class II ERVs infiltrated amphibian genomes recently and may still have infectious members. Remarkably, certain amphibian class II ERVs can be expressed in diverse tissues. This discovery closes the major gap in the retroviral fossil record of class II ERVs and provides important insights into the evolution of class II ERVs in vertebrates.IMPORTANCE Class II retroviruses, largely distributed among mammals and birds, are of particular importance for medicine and economics. Class II ERVs have been discovered in a range of vertebrates, with the exception of amphibians, which are known only to possess class I and class III-like ERVs. Here, for the first time, we discovered class II ERVs in amphibians. We also revealed that the ERVs may originate from alpha- or betaretroviruses. We revealed that class II ERVs were integrated into amphibian genomes recently and certain amphibian class II ERVs can be expressed in diverse tissues. Our discovery closes the major gap in the retroviral fossil record of class II ERVs, and also indicates that amphibians may be still infected by class II retroviruses.

The first complete mitochondrial genome sequence of Nanorana parkeri and Nanorana ventripunctata (Amphibia: Anura: Dicroglossidae), with related phylogenetic analyses

Members of the Nanorana genus (family Dicroglossidae) are often referred to as excellent model species with which to study amphibian adaptations to extreme environments and also as excellent keystone taxa for providing insights into the evolution of the Dicroglossidae. However, a complete mitochondrial genome is currently only available for Nanorana pleskei. Thus, we analyzed the complete mitochondrial genomes of Nanorana parkeri and Nanorana ventripunctata to investigate their evolutionary relationships within Nanorana and their phylogenetic position in the family Dicroglossidae. Our results showed that the genomes of N. parkeri (17,837 bp) and N. ventripunctata (18,373 bp) encode 13 protein‐coding genes (PCGs), two ribosomal RNA genes, 23 transfer RNA (tRNA) genes, and a noncoding control region. Overall sequences and genome structure of the two species showed high degree of similarity with N. pleskei, although the motif structures and repeat sequences of the putative control region showed clear differences among these three Nanorana species. In addition, a tandem repeat of the tRNA‐Met gene was found located between the tRNA‐Gln and ND2 genes. On both the 5′ and 3′‐sides, the control region possessed distinct repeat regions; however, the CSB‐2 motif was not found in N. pleskei. Based on the nucleotide sequences of 13 PCGs, our phylogenetic analyses, using Bayesian inference and maximum‐likelihood methods, illustrate the taxonomic status of Nanorana with robust support showing that N. ventripunctata and N. pleskei are more closely related than they are to N. parkeri. In conclusion, our analyses provide a more robust and reliable perspective on the evolutionary history of Dicroglossidae than earlier analyses, which used only a single species (N. pleskei).

Next-Generation Mitogenomics: A Comparison of Approaches Applied to Caecilian Amphibian Phylogeny

Mitochondrial genome (mitogenome) sequences are being generated with increasing speed due to the advances of next-generation sequencing (NGS) technology and associated analytical tools. However, detailed comparisons to explore the utility of alternative NGS approaches applied to the same taxa have not been undertaken. We compared a 'traditional' Sanger sequencing method with two NGS approaches (shotgun sequencing and non-indexed, multiplex amplicon sequencing) on four different sequencing platforms (Illumina's HiSeq and MiSeq, Roche's 454 GS FLX, and Life Technologies' Ion Torrent) to produce seven (near-) complete mitogenomes from six species that form a small radiation of caecilian amphibians from the Seychelles. The fastest, most accurate method of obtaining mitogenome sequences that we tested was direct sequencing of genomic DNA (shotgun sequencing) using the MiSeq platform. Bayesian inference and maximum likelihood analyses using seven different partitioning strategies were unable to resolve compellingly all phylogenetic relationships among the Seychelles caecilian species, indicating the need for additional data in this case.

Life-history evolution and mitogenomic phylogeny of caecilian amphibians

We analyze mitochondrial genomes to reconstruct a robust phylogenetic framework for caecilian amphibians and use this to investigate life-history evolution within the group. Our study comprises 45 caecilian mitochondrial genomes (19 of them newly reported), representing all families and 27 of 32 currently recognized genera, including some for which molecular data had never been reported. Support for all relationships in the inferred phylogenetic tree is high to maximal, and topology tests reject all investigated alternatives, indicating an exceptionally robust molecular phylogenetic framework of caecilian evolution consistent with current morphology-based supraspecific classification. We used the mitogenomic phylogenetic framework to infer ancestral character states and to assess correlation among three life-history traits (free-living larvae, viviparity, specialized pre-adult or vernal teeth), each of which occurs only in some caecilian species. Our results provide evidence that an ancestor of the Seychelles caecilians abandoned direct development and re-evolved a free-living larval stage. This study yields insights into the concurrent evolution of direct development and of vernal teeth in an ancestor of Teresomata that likely gave rise to skin-feeding (maternal dermatophagy) behavior and subsequently enabled evolution of viviparity, with skin feeding possibly a homologous precursor of oviduct feeding in viviparous caecilians.

Molecular phylogenetic relationships of pangasiid and schilbid catfishes in Thailand

In this study, the phylogenetic relationships among 13 pangasiids and six schilbids of Thailand were reconstructed based on the almost complete mitochondrial cytochrome b (cyt b), 12S rRNA, tRNA-Val and 16S rRNA, as well as the partial nuclear recombination-activating gene 1 (rag1) sequences by using the maximum likelihood and the Bayesian inference methods of phylogenetic reconstruction. The reconstructed phylogeny based on the concatenated sequence data set recovered Pangasiidae and Schilbidae as reciprocally monophyletic groups. Within Pangasiidae, four major clades were recovered, which according to the cyt b genetic distances can be categorized into four genera: Pangasius, Pseudolais, Helicophagus and Pangasianodon. The genus Pangasianodon was strongly supported as the most basal taxon within pangasiids, whereas Pseudolais + Helicophagus were recovered as a sister group of Pangasius. Within the latter, the giant pangasius Pangasius sanitwongsei was recovered as a sister group of the spot pangasius Pangasius larnaudii, Pangasius krempfi as a sister group of Pangasius nasutus + Pangasius conchophilus and Pangasius polyuranodon as a sister species of Pangasius macronema. Other internal phylogenetic relationships, however, were unresolved. Within Schilbidae, Pseudeutropius was supported as the most basal lineage. Eutropiichthys was recovered as a sister group of Clupisoma. The enigmatic Clupisoma sinense was recognized as more closely related to Laides longibarbis than to Clupisoma prateri. Thus, based on the cyt b genetic distances, a recategorization of C. sinense to the genus Laides is suggested. On the basis of a relaxed clock fossil calibration, the divergence of pangasiids and schilbids was estimated to have occurred 14·93 million years before present (b.p.) during the Miocene epoch. The separation between Pangasiidae and Schilbidae took place c. 13·12 Mb.p. during the early middle Miocene. The estimated divergence time of pangasiids is similar to the age of the calibrated fossil, Cetopangasius chaetobranchus, which was discovered in north-central Thailand. This suggests that the oldest pangasiid ancestor diverged into diverse genera in the area.

Complete nucleotide sequence and gene organization of the mitochondrial genome of Paa spinosa (Anura: Ranoidae)

The mt genome of Paa spinosa (Anura: Ranoidae) is a circular molecule of 18,012 bp in length, containing 38 genes (including an extra copy of tRNA-Met gene). This mt genome is characterized by three distinctive features: a cluster of rearranged tRNA genes (LTPF tRNA gene cluster), a tandem duplication of tRNA-Met gene (Met1 and Met2), and distinct repeat regions at both 5' and 3'-sides in the control region. Comparing the locations and the sequences of all tRNA-Met genes among Ranoidae, and constructing NJ tree of the nucleotide of those tRNA-Met genes, we suggested a tandem duplication of tRNA-Met gene can be regarded as a synapomorphy of Dicroglossinae. To further investigate the phylogenetic relationships of anurans, phylogenetic analyses (BI, ML and MP) based on the nucleotide dataset and the corresponding amino acid dataset of 11 protein-coding genes (except ND5 and ATP8) arrived at the similar topology.

The complete mitochondrial genome of the large-headed frog, Limnonectes bannaensis (Amphibia: Anura), and a novel gene organization in the vertebrate mtDNA

We determined the complete nucleotide sequence of the mitochondrial (mt) genome of the large-headed frog, Limnonectes bannaensis (Amphibia, Anura) by using polymerase chain reaction (PCR). The entire mtDNA sequence is 16,867 bp in length with a novel case of tRNAs in vertebrates. This mt genome is characterized by three distinctive features: (1) a tandem duplication of tRNA(Met) gene is observed, (2) the tRNA(Ala), tRNA(Asn), tRNA(Cys) and tRNA(Glu) genes coded on the L-strand are absent from the L. bannaensis mtDNA, the tRNA(Cys) and tRNA(Glu) genes change into tRNA pseudogenes by reason of degenerative anticodon, and a noncoding sequence of 206 nt long (NC1) has replaced the original position of other two tRNAs, (3) besides NC1, another three noncoding spacers (NC2-4) longer than 50 bp are found in the broken WANCY region and the region NC3-ND5-NC4-ND6-PsiE-Cytb-CR of the new sequence. These features could be explained by a model of gene duplication and deletion. The new sequence data was used to assess the phylogenetic relationships among 25 species of Anura using neighbor-joining, Bayesian, and maximum likelihood methods, and the phylogenetic tree shows the rice frog Fejervarya limnocharis is closest to L. bannaensis in the study.

Gene expression and tissue distribution of cytoglobin and myoglobin in the Amphibia and Reptilia: possible compensation of myoglobin with cytoglobin in skeletal muscle cells of anurans that lack the myoglobin gene

Cytoglobin (Cygb), a recently discovered vertebrate cytoplasmic heme-binding globin, is considered to be in a clade with vertebrate myoglobin (Mb), which is exclusively distributed in the cytoplasm of cardiac and skeletal muscles as an oxygen storage protein. GenBank databases (NCBI and JGI) and gene synteny analyses showed the absence of the Mb gene (mb) in two anuran amphibians, Xenopus laevis and X. tropicalis. Here we conducted comparative studies on the gene expression and tissue distribution of Cygb and Mb in anuran and reptilian tissues. Cygb and Mb genes were cloned from a reptile, iguana (Iguana iguana). Two types of cygb (cygb-1 and -2) were cloned, with lengths of 1066 and 1034 bp, and 196 and 193 amino acid residues, respectively. Their nucleotide and amino acid sequence identities were 90 and 87%, respectively. The Mb gene covered 1416 bp with an open reading frame of 465 bp, giving rise to a 154 amino acid protein. The distal ligand-binding histidine at E7, the proximal heme-binding histidine at F8, and the phenylalanine residue at CD1 were conserved in Mb and Cygb. The nucleotide and amino acid sequence identity of I. iguana cygb-1 and cygb-2 against X. laevis cygb were approximately 67% and 65%, respectively. RT-PCR demonstrated that X. laevis cygb was uniquely expressed in the heart and skeletal muscles, and faintly in the liver and spleen, which was quite contrasted with Iguana and the other vertebrates, where mb is exclusively expressed in the heart and skeletal muscles. Immunohistochemical analyses showed the distribution of Cygb in the cytoplasm of skeletal muscle cells. Interestingly, Cygb in the heart was localized in the nuclei. Considering the absence of mb in the Anura, we hypothesize that Cygb in muscle cells of anurans compensates for the lack of Mb for the storage and intracellular transportation of oxygen.

The complete mitochondrial genome of a basal teleost, the Asian arowana (Scleropages formosus, Osteoglossidae)

BACKGROUND

Mitochondrial DNA-derived sequences have become popular markers for evolutionary studies, as their comparison may yield significant insights into the evolution of both the organisms and their genomes. From the more than 24,000 teleost species, only 254 complete mtDNA sequences are available (GenBank status on 06 Sep 2006). In this paper, we report the complete mitochondrial genome sequence of Asian arowana, a basal bonytongue fish species, which belongs to the order of Osteoglossiformes.

RESULTS

The complete mitochondrial genomic sequence (mtDNA) of Asian arowana (Scleropages formosus) was determined by using shotgun sequencing method. The length of Asian arowana mtDNA is ca. 16,650 bp (its variation is due to polymorphic repeats in the control region), containing 13 protein-coding genes, 22 tRNA and 2 rRNA genes. Twelve of the thirteen protein coding genes were found to be encoded by the heavy strand in the order typically observed for vertebrate mitochondrial genomes, whereas only nad6 was located on the light strand. An interesting feature of Asian arowana mitogenome is that two different repeat arrays were identified in the control region: a 37 bp tandem repeat at the 5' end and an AT-type dinucleotide microsatellite at the 3' end. Both repeats show polymorphism among the six individuals tested; moreover the former one is present in the mitochondrial genomes of several other teleost groups. The TACAT motif described earlier only from mammals and lungfish was found in the tandem repeat of several osteoglossid and eel species. Phylogenetic analysis of fish species representing Actinopterygii and Sarcopterygii taxa has shown that the Asian arowana is located near the baseline of the teleost tree, confirming its status among the ancestral teleost lineages.

CONCLUSION

The mitogenome of Asian arowana is very similar to the typical vertebrate mitochondrial genome in terms of gene arrangements, codon usage and base composition. However its control region contains two different types of repeat units at both ends, an interesting feature that to our knowledge has never been reported before for other vertebrate mitochondrial control regions. Phylogenetic analysis using the complete mtDNA sequence of Asian arowana confirmed that it belongs to an ancestral teleost lineage.

Structure and evolution of the mitochondrial control region and flanking sequences in the European cave salamander Proteus anguinus

The European cave salamander Proteus anguinus Laurenti 1768 is one of the best-known subterranean animals, yet its evolutionary history and systematic relationships remain enigmatic. This is the first comprehensive study on molecular evolution within the taxon, using an mtDNA segment containing the control region (CR) and adjacent sequences. Two to seven tandem repeats of 24-32 bp were found in the intergenic spacer region (VNTR1), and three, four or six repeats, 59-77 bp each, in the 3' end of the CR (VNTR2). Different molecular mechanisms account for VNTR2 formation in different lineages of Proteus. The overall CR variation was lower than that of the spacer region, the 3' end of the cytb gene, or the tRNA genes. Individual genes and the concatenated non-repetitive sequences produced similar, well resolved maximum likelihood, Bayesian inference and parsimony trees. The numbers of repeat elements as well as the genealogy of the VNTR2 repeat units were mostly inconsistent with the groupings of the non-repetitive sequences. Different degrees of repeat array homogenization were detected in all major groups. Orthology was established for the first and the second VNTR2 elements of some populations. These two copies may therefore be used for analyses at the population level. The pattern of CR sequence variation points to strong genetic isolation of hydrographically separated populations. Genetic separation of the major groups of populations is incongruent with the current division into subspecies.

Complete nucleotide sequence of the mitochondrial genome of Schlegel's tree frog Rhacophorus schlegelii (family Rhacophoridae): duplicated control regions and gene rearrangements

The complete nucleotide sequence (21,359 bp) of the mitochondrial DNA of the rhacophorid frog Rhacophorus schlegelii was determined. The gene content, nucleotide composition, and codon usage of this genome corresponded to those typical of vertebrates. However, the Rh. schlegelii genome was unusually large due to the inclusion of two control regions and the accumulation of lengthy repetitive sequences in these regions. The two control regions had 97% sequence similarity over 1,510 bp, suggesting the occurrence of concerted sequence evolution. Comparison of the gene organizations among anuran species revealed that the mitochondrial gene arrangement of Rh. schlegelii diverged from that of typical vertebrates but was similar to that of Buergeria buergeri. The positions of the tRNA-Leu(CUN) and tRNA-Thr genes were exchanged between Rh. schlegelii and B. buergeri. Based on parsimonious consideration and the basal phylogenetic position of B. buergeri, these genes seemed to have been rearranged in an ancestral lineage leading to Rh. schlegelii.

Cranial musculature in the larva of the caecilian,Ichthyophis kohtaoensis (Lissamphibia: Gymnophiona)

Within the Gymnophiona (caecilians) oviparous species with biphasic life-cycles possess a free living semi-aquatic larval stage that feeds in aquatic habitats. The larvae pass through a metamorphosis to a purely terrestrial adult stage. It is likely that the cranial morphology of caecilian larvae has specializations for aquatic feeding. However, little is known about the cranial morphology, and the cranial musculature is especially neglected in the literature. This study provides a detailed description of the jaw and hyobranchial musculature in larval stages of a caecilian. We studied late embryonic and early larval specimens of Ichthyophis kohtaoensis. Furthermore, we compared and homologized the cranial muscles found in larval I. kohtaoensis with the muscles described for adult caecilians. Most cranial muscles of larval I. kohtaoensis are also present in the adult, except for the m. levator mandibulae externus and the m. subarcualis obliquus II. Our results were compared with the data available for larval frogs and salamanders in order to hypothesize the cranial musculature in the larva of the most recent common ancestor of the Lissamphibia. Larval caecilians, frog tadpoles, and salamander larvae share many characters in their cranial musculature, which, consequently, can be assigned to the lissamphibian ground pattern. However, the m. pterygoideus and the m. levator quadrati are unique to the Gymnophiona.

A Hotspot of Gene Order Rearrangement by Tandem Duplication and Random Loss in the Vertebrate Mitochondrial Genome

Most reported examples of change in vertebrate mitochondrial (mt) gene order could be explained by a tandem duplication followed by random loss of redundant genes (tandem duplication-random loss [TDRL] model). Under this model of evolution, independent loss of genes arising from a single duplication in an ancestral species are predicted, and remnant pseudogenes expected, intermediate states that may remain in rearranged genomes. However, evidence for this is rare and largely scattered across vertebrate lineages. Here, we report new derived mt gene orders in the vertebrate "WANCY" region of four closely related caecilian amphibians. The novel arrangements found in this genomic region (one of them is convergent with the derived arrangement of marsupials), presence of pseudogenes, and positions of intergenic spacers fully satisfy predictions from the TDRL model. Our results, together with comparative data for the available vertebrate complete mt genomes, provide further evidence that the WANCY genomic region is a hotspot for gene order rearrangements and support the view that TDRL is the dominant mechanism of gene order rearrangement in vertebrate mt genomes. Convergent gene rearrangements are not unlikely in hotspots of gene order rearrangement by TDRL.

Mitogenomic perspectives on the origin and phylogeny of living amphibians

Establishing the relationships among modern amphibians (lissamphibians) and their ancient relatives is necessary for our understanding of early tetrapod evolution. However, the phylogeny is still intractable because of the highly specialized anatomy and poor fossil record of lissamphibians. Paleobiologists are still not sure whether lissamphibians are monophyletic or polyphyletic, and which ancient group (temnospondyls or lepospondyls) is most closely related to them. In an attempt to address these problems, eight mitochondrial genomes of living amphibians were determined and compared with previously published amphibian sequences. A comprehensive molecular phylogenetic analysis of nucleotide sequences yields a highly resolved tree congruent with the traditional hypotheses (Batrachia). By using a molecular clock-independent approach for inferring dating information from molecular phylogenies, we present here the first molecular timescale for lissamphibian evolution, which suggests that lissamphibians first emerged about 330 million years ago. By observing the fit between molecular and fossil times, we suggest that the temnospondyl-origin hypothesis for lissamphibians is more credible than other hypotheses. Moreover, under this timescale, the potential geographic origins of the main living amphibian groups are discussed: (i) advanced frogs (neobatrachians) may possess an Africa-India origin; (ii) salamanders may have originated in east Asia; (iii) the tropic forest of the Triassic Pangaea may be the place of origin for the ancient caecilians. An accurate phylogeny with divergence times can be also helpful to direct the search for "missing" fossils, and can benefit comparative studies of amphibian evolution.

Complete Mitochondrial DNA Sequences of Six Snakes: Phylogenetic Relationships and Molecular Evolution of Genomic Features

Complete mitochondrial DNA (mtDNA) sequences were determined for representative species from six snake families: the acrochordid little file snake, the bold boa constrictor, the cylindrophiid red pipe snake, the viperid himehabu, the pythonid ball python, and the xenopeltid sunbeam snake. Thirteen protein-coding genes, 22 tRNA genes, 2 rRNA genes, and 2 control regions were identified in these mtDNAs. Duplication of the control region and translocation of the tRNALeu gene were two notable features of the snake mtDNAs. The duplicate control regions had nearly identical nucleotide sequences within species but they were divergent among species, suggesting concerted sequence evolution of the two control regions. In addition, the duplicate control regions appear to have facilitated an interchange of some flanking tRNA genes in the viperid lineage. Phylogenetic analyses were conducted using a large number of sites (9570 sites in total) derived from the complete mtDNA sequences. Our data strongly suggested a new phylogenetic relationship among the major families of snakes: ((((Viperidae, Colubridae), Acrochordidae), (((Pythonidae, Xenopeltidae), Cylindrophiidae), Boidae)), Leptotyphlopidae). This conclusion was distinct from a widely accepted view based on morphological characters in denying the sister-group relationship of boids and pythonids, as well as the basal divergence of nonmacrostomatan cylindrophiids. These results imply the significance to reconstruct the snake phylogeny with ample molecular data, such as those from complete mtDNA sequences.

Phylogeny of caecilian amphibians (Gymnophiona) based on complete mitochondrial genomes and nuclear RAG1

We determined the complete nucleotide sequence of the mitochondrial (mt) genome of five individual caecilians (Amphibia: Gymnophiona) representing five of the six recognized families: Rhinatrema bivittatum (Rhinatrematidae), Ichthyophis glutinosus (Ichthyophiidae), Uraeotyphlus cf. oxyurus (Uraeotyphlidae), Scolecomorphus vittatus (Scolecomorphidae), and Gegeneophis ramaswamii (Caeciliidae). The organization and size of these newly determined mitogenomes are similar to those previously reported for the caecilian Typhlonectes natans (Typhlonectidae), and for other vertebrates. Nucleotide sequences of the nuclear RAG1 gene were also determined for these six species of caecilians, and the salamander Mertensiella luschani atifi. RAG1 (both at the amino acid and nucleotide level) shows slower rates of evolution than almost all mt protein-coding genes (at the amino acid level). The new mt and nuclear sequences were compared with data for other amphibians and subjected to separate and combined phylogenetic analyses (Maximum Parsimony, Minimum Evolution, Maximum Likelihood, and Bayesian Inference). All analyses strongly support the monophyly of the three amphibian Orders. The Batrachia hypothesis (Gymnophiona, (Anura, Caudata) receives moderate or good support depending on the method of analysis. Within Gymnophiona, the optimal tree (Rhinatrema, (Ichthyophis, Uraeotyphlus), (Scolecomorphus, (Gegeneophis Typhlonectes) agrees with the most recent morphological and molecular studies. The sister group relationship between Rhinatrematidae and all other caecilians, that between Ichthyophiidae and Uraeotyphlidae, and the monophyly of the higher caecilians Scolecomorphidae+Caeciliidae+Typhlonectidae, are strongly supported, whereas the relationships among the higher caecilians are less unambiguously resolved. Analysis of RAG1 is affected by a spurious local rooting problem and associated low support that is ameliorated when outgroups are excluded. Comparisons of trees using the non-parametric Templeton, Kishino-Hasegawa, Approximately Unbiased, and Shimodaira-Hasegawa tests suggest that the latter may be too conservative.

The complete mitochondrial genome of a tree frog, Polypedates megacephalus (Amphibia: Anura: Rhacophoridae), and a novel gene organization in living amphibians

In this study, we have determined the complete sequence of the mitochondrial genome of an Old World tree frog Polypedates megacephalus (Anura: Rhacophoridae) by using a long polymerase chain reaction (PCR) technique and shotgun strategy of sequencing. The entire mtDNA sequence is 16,473 nt long with a novel mitogenomic gene organization in amphibians. Unlike other neobatrachian frogs, the transfer ribonucleic acid (tRNA)-Leu(CUN) and tRNA-Thr genes exchange their positions in P. megacephalus and form a Thr-Leu(CUN)-Pro-Phe tRNA gene tetrad. Moreover, we found that the ATP8 gene was replaced by a noncoding sequence of 853 nt long and that the ND5 gene was absent in the new mitogenome. These peculiar features of P. megacephalus mtDNA were further studied among related anuran species by PCR amplification. The new sequence data was used to assess the phylogenetic relationships of the three living amphibian orders using neighbor-joining, maximum likelihood, and Bayesian methods. In agreement with most morphological studies, phylogenetic analyses of a whole mitochondrial genome data set suggest a close relationship between salamanders and frogs. Moreover, using a molecular clock-independent Bayesian approach for inferring dating information from molecular phylogenies, we have provided a rough timescale for living amphibian evolution. This timescale provides a working framework for future paleontological researches on amphibian evolution and improves our understanding of the evolutionary history of modern amphibians.

The mitochondrial genome organization of the rice frog, Fejervarya limnocharis (Amphibia: Anura): a new gene order in the vertebrate mtDNA

The mitochondrial DNA of the rice frog, Fejervarya limnocharis (Amphibia, Anura), was obtained using long-and-accurate polymerase chain reaction (LA-PCR) combining with subcloning method. The complete nucleotide sequence (17,717 bp) of mitochondrial genome was determined subsequently. This mitochondrial genome is characterized by four distinctive features: the translocation of ND5 gene, a cluster of rearranged tRNA genes (tRNAThr, tRNAPro, tRNA(Leu CUN)), a tandem duplication of tRNAMet gene, and eight large 89-bp tandem repeats in the control region, as well as three short noncoding regions containing two repeated motifs existing in the gene cluster of ND5/tRNAThr/tRNAPro/tRNALeu/tRNAPhe. The tandem duplication of gene regions followed by deletions of supernumerary genes can be invoked to explain the shuffling of tRNAMet and a cluster of tRNA and ND5 genes, as observed in this study. Both ND5 gene translocation and tandem duplication of tRNAMet were first observed in the vertebrate mitochondrial genomes.

Integrating reptilian herpesviruses into the family herpesviridae

The phylogeny of reptilian herpesviruses (HVs) relative to mammalian and avian HVs was investigated by using available gene sequences and by alignment of encoded amino acid sequences and derivation of trees by maximum-likelihood and Bayesian methods. Phylogenetic loci were obtained for green turtle HV (GTHV) primarily on the basis of DNA polymerase (POL) and DNA binding protein sequences, and for lung-eye-trachea disease-associated HV (LETV) primarily from its glycoprotein B sequence; both have nodes on the branch leading to recognized species in the Alphaherpesvirinae subfamily and should be regarded as new members of that subfamily. A similar but less well defined locus was obtained for an iguanid HV based on a partial POL sequence. On the basis of short POL sequences (around 60 amino acid residues), it appeared likely that GTHV and LETV belong to a private clade and that three HVs of gerrhosaurs (plated lizards) are associated with the iguanid HV. Based on phylogenetic branching patterns for mammalian HV lineages that mirror those of host lineages, we estimated a date for the HV tree's root of around 400 million years ago. Estimated dates for branching events in the development of reptilian, avian, and mammalian Alphaherpesvirinae lineages could plausibly be accounted for in part but not completely by ancient coevolution of these virus lines with reptilian lineages and with the development of birds and mammals from reptilian progenitors.

Complete Mitochondrial Genome Sequences of the South American and the Australian Lungfish: Testing of the Phylogenetic Performance of Mitochondrial Data Sets for Phylogenetic Problems in Tetrapod Relationships

We determined the complete nucleotide sequences (16403 and 16572 base pairs, respectively) of the mitochondrial genomes of the South American lungfish, Lepidosiren paradoxa, and the Australian lungfish, Neoceratodus forsteri (Sarcopterygii, Dipnoi). The mitochondrial DNA sequences were established in an effort to resolve the debated evolutionary positions of the lungfish and the coelacanth relative to land vertebrates. Previous molecular phylogenetic studies based on complete mtDNA sequences, including only the African lungfish, Protopterus dolloi, sequence were able to strongly reject the traditional textbook hypothesis that coelacanths are the closest relatives of land vertebrates. However, these studies were unable to statistically significantly distinguish between the two remaining scenarios: lungfish as the closest relatives to land vertebrates and lungfish and coelacanths jointly as their sister group (Cao et al. 1998; Zardoya et al. 1998; Zardoya and Meyer 1997a). Lungfish, coelacanths, and the fish ancestors of the tetrapod lineage all originated within a short time window of about 20 million years, back in the early Devonian (about 380 to 400 million years ago). This short divergence time makes the determination of the phylogenetic relationships among these three lineages difficult. In this study, we attempted to break the long evolutionary branch of lungfish, in an effort to better resolve the phylogenetic relationships among the three extant sarcopterygian lineages. The gene order of the mitochondrial genomes of the South American and Australian lungfish conforms to the consensus gene order among gnathostome vertebrates. The phylogenetic analyses of the complete set of mitochondrial proteins (without ND6) suggest that the lungfish are the closest relatives of the tetrapods, although the support in favor of this scenario is not statistically significant. The two other smaller data sets (tRNA and rRNA genes) give inconsistent results depending on the different reconstruction methods applied and cannot significantly rule out any of the three alternative hypotheses. Nuclear protein-coding genes, which might be better phylogenetic markers for this question, support the lungfish-tetrapod sister-group relationship (Brinkmann et al. 2004).

Molecular systematics of primary reptilian lineages and the tuatara mitochondrial genome

We provide phylogenetic analyses for primary Reptilia lineages including, for the first time, Sphenodon punctatus (tuatara) using data from whole mitochondrial genomes. Our analyses firmly support a sister relationship between Sphenodon and Squamata, which includes lizards and snakes. Using Sphenodon as an outgroup for select squamates, we found evidence indicating a sister relationship, among our study taxa, between Serpentes (represented by Dinodon) and Varanidae. Our analyses support monophyly of Archosauria, and a sister relationship between turtles and archosaurs. This latter relationship is congruent with a growing set of morphological and molecular analyses placing turtles within crown Diapsida and recognizing them as secondarily anapsid (lacking a skull fenestration). Inclusion of Sphenodon, as the only surviving member of Sphenodontia (with fossils from the mid-Triassic), helps to fill a sampling gap within previous analyses of reptilian phylogeny. We also report a unique configuration for the mitochondrial genome of Sphenodon, including two tRNA(Lys) copies and an absence of ND5, tRNA(His), and tRNA(Thr) genes.

Palaeoecology of triassic stem turtles sheds new light on turtle origins

Competing hypotheses of early turtle evolution contrast sharply in implying very different ecological settings-aquatic versus terrestrial-for the origin of turtles. We investigate the palaeoecology of extinct turtles by first demonstrating that the forelimbs of extant turtles faithfully reflect habitat preferences, with short-handed turtles being terrestrial and long-handed turtles being aquatic. We apply this metric to the two successive outgroups to all living turtles with forelimbs preserved, Proganochelys quenstedti and Palaeochersis talampayensis, to discover that these earliest turtle outgroups were decidedly terrestrial. We then plot the observed distribution of aquatic versus terrestrial habits among living turtles onto their hypothesized phylogenies. Both lines of evidence indicate that although the common ancestor of all living turtles was aquatic, the earliest turtles clearly lived in a terrestrial environment. Additional anatomical and sedimentological evidence favours these conclusions. The freshwater aquatic habitat preference so characteristic of living turtles cannot, consequently, be taken as positive evidence for an aquatic origin of turtles, but must rather be considered a convergence relative to other aquatic amniotes, including the marine sauropterygians to which turtles have sometimes been allied.

Complete nucleotide sequence and gene rearrangement of the mitochondrial genome of the bell-ring frog, Buergeria buergeri (family Rhacophoridae)

In this study we determined the complete nucleotide sequence (19,959 bp) of the mitochondrial DNA of the rhacophorid frog Buergeria buergeri. The gene content, nucleotide composition, and codon usage of B. buergeri conformed to those of typical vertebrate patterns. However, due to an accumulation of lengthy repetitive sequences in the D-loop region, this species possesses the largest mitochondrial genome among all the vertebrates examined so far. Comparison of the gene organizations among amphibian species (Rana, Xenopus, salamanders and caecilians) revealed that the positioning of four tRNA genes and the ND5 gene in the mtDNA of B. buergeri diverged from the common vertebrate gene arrangement shared by Xenopus, salamanders and caecilians. The unique positions of the tRNA genes in B. buergeri are shared by ranid frogs, indicating that the rearrangements of the tRNA genes occurred in a common ancestral lineage of ranids and rhacophorids. On the other hand, the novel position of the ND5 gene seems to have arisen in a lineage leading to rhacophorids (and other closely related taxa) after ranid divergence. Phylogenetic analysis based on nucleotide sequence data of all mitochondrial genes also supported the gene rearrangement pathway.

Complete nucleotide sequence of the mitochondrial genome of a salamander, Mertensiella luschani

The complete nucleotide sequence (16,650 bp) of the mitochondrial genome of the salamander Mertensiella luschani (Caudata, Amphibia) was determined. This molecule conforms to the consensus vertebrate mitochondrial gene order. However, it is characterized by a long non-coding intervening sequence with two 124-bp repeats between the tRNA(Thr) and tRNA(Pro) genes. The new sequence data were used to reconstruct a phylogeny of jawed vertebrates. Phylogenetic analyses of all mitochondrial protein-coding genes at the amino acid level recovered a robust vertebrate tree in which lungfishes are the closest living relatives of tetrapods, salamanders and frogs are grouped together to the exclusion of caecilians (the Batrachia hypothesis) in a monophyletic amphibian clade, turtles show diapsid affinities and are placed as sister group of crocodiles+birds, and the marsupials are grouped together with monotremes and basal to placental mammals. The deduced phylogeny was used to characterize the molecular evolution of vertebrate mitochondrial proteins. Amino acid frequencies were analyzed across the main lineages of jawed vertebrates, and leucine and cysteine were found to be the most and least abundant amino acids in mitochondrial proteins, respectively. Patterns of amino acid replacements were conserved among vertebrates. Overall, cartilaginous fishes showed the least variation in amino acid frequencies and replacements. Constancy of rates of evolution among the main lineages of jawed vertebrates was rejected.

The complete mitochondrial genome of the Chinese giant salamander, Andrias davidianus (Amphibia: Caudata)

The mitochondrial genome of the Chinese giant salamander Andrias davidianus was isolated using the long-and-accurate polymerase chain reaction (LA PCR) method. The sequencing work adopted the shotgun strategy accompanying with seven internal primers to cover the gaps where overlapping clones were not available. The entire mtDNA sequence is 16,503 bp long, with a gene content of 13 protein-coding, two ribosomal RNA and 22 transfer RNA genes, and order identical to that observed in most other vertebrates except for an additional 318 bp non-coding sequence between tRNA-Thr and tRNA-Pro genes. In order to carry out molecular phylogenetic analyses, all 13 protein sequences deduced from whole mitochondrial genomes for eight vertebrate species (six amphibians, two lobe-finned fishes) were combined to a single data set. This data set was refined by the program Gblocks using a stringent parameter setting and then subjected to MP, ML and NJ analyses. Thus phylogenetic relationships among living amphibians were discussed.

Complete mitochondrial DNA of the hagfish, Eptatretus burgeri: the comparative analysis of mitochondrial DNA sequences strongly supports the cyclostome monophyly

The phylogenetic position of cyclostomes, i.e., the relationships between hagfishes, lampreys, and jawed vertebrates is an unresolved problem. Anatomical data support the paraphyly of cyclostomes, whereas nuclear genes data support monophyly of cyclostomes. Previous results obtained using mitochondrial DNA are ambiguous, presumably due to a lack of informative sequences. By adding the complete mtDNA of a hagfish, Eptatretus burgeri, we have generated a novel data set for sequences of hagfishes and of lampreys. The addition of this mtDNA sequence to the 12 taxa we have already used becomes sufficient to obtain unambiguous results. This data set, which includes sequences of mtDNA of animals closely related to the lamprey/hagfish node, was used in a phylogenetic analysis with two independent statistical approaches and unequivocally supported the monophyly of cyclostomes. Thus molecular data, i.e., our results and those obtained using nuclear genes, conclude that hagfishes and lampreys form a clade.

Molecular phylogeny of French Guiana Hylinae: implications for the systematic and biodiversity of the Neotropical frogs

In this study we used nucleotide sequences from a segment of mitochondrial 16S ribosomal DNA gene to investigate the evolutionary relationships of some French Guiana Hylinae. New sequences, representing the members of different French Guiana frogs-five specimens of the Scinax genus, two Hyla, one Osteocephalus, one Hyalinobatrachium and two Rana as out-group-were examined. In addition, 26 sequences available from GenBank database representing the other subfamilies of the Hylidae were added to our study. This work allowed us to clarify relationships within the four hylids subfamilies (Pelodryadinae, Phyllomedusinae, Hemiphractinae and Hylinae) and the phylogenetic placement of the enigmatic Scinax genus within the Hylidae. We found that: (1) the Scinax genus displays a high level of differentiation in comparison to two other genera (Litoria and Hyla) belonging to 'Hylidae' family; (2) the Hylinae are paraphyletic given the position of the Litoria, which was the sister-group of the Hyla and the Osteocephalus genera; (3) the anterior works and our results (based on two different data sets) showed the paraphyly of the Hylidae questioning the validity of this family; (4) the reassessment of these different taxonomic groups will induce a huge implication on the estimation (past, present and future) of the biodiversity (in Neotropical frogs).

Evolutionary dynamics in a novel L2 clade of non-LTR retrotransposons in Deuterostomia

The evolution of the novel L2 clade of non-long terminal repeat (LTR) retrotransposons and their evolutionary dynamics in Deuterostomia has been examined. The short-term evolution of long interspersed nuclear element 2s (LINE2s) has been studied in 18 reptilian species by analysis of a PCR amplified 0.7-kb fragment encoding the palm/fingers subdomain of reverse transcriptase (RT). Most of the reptilian LINE2s examined are inactive since they contain multiple stop codons, indels, or frameshift mutations that disrupt the RT. Analysis of reptilian LINE2s has shown a high degree of sequence divergence and an unexpectedly large number of deletions. The evolutionary dynamics of LINE2s in reptiles has been found to be complex. LINE2s are shown to form a novel clade of non-LTR retrotransposons that is well separated from the CR1 clade. This novel L2 clade is more widely distributed than previously thought, and new representatives have been discovered in echinoderms, insects, teleost fishes, Xenopus, Squamata, and marsupials. There is an apparent absence of LINE2s from different vertebrate classes, such as cartilaginous fishes, Archosauria (birds and crocodiles), and turtles. Whereas the LINE2s are present in echinoderms and teleost fishes in a conserved form, in most tetrapods only highly degenerated pseudogenes can be found. The predominance of inactive LINE2s in Tetrapoda indicates that, in the host genomes, only inactive copies are still present. The present data indicate that the vertical inactivation of LINE2s might have begun at the time of Tetrapoda origin, 400 MYA. The evolutionary dynamics of the L2 clade in Deuterostomia can be described as a gradual vertical inactivation in Tetrapoda, stochastic loss in Archosauria and turtles, and strict vertical transmission in echinoderms and teleost fishes.

Complete nucleotide sequence and gene rearrangement of the mitochondrial genome of the Japanese pond frog Rana nigromaculata

In this study, we determined the complete nucleotide sequence of the mitochondrial genome of the Japanese pond frog Rana nigromaculata. The length of the sequence of the frog was 17,804 bp, though this was not absolute due to length variation caused by differing numbers of repetitive units in the control regions of individual frogs. The gene content, base composition, and codon usage of the Japanese pond frog conformed to those of typical vertebrate patterns. However, the comparison of gene organization between three amphibian species (Rana, Xenopus and caecilian) provided evidence that the gene arrangement of Rana differs by four tRNA gene positions from that of Xenopus or caecilian, a common gene arrangement in vertebrates. These gene rearrangements are presumed to have occurred by the tandem duplication of a gene region followed by multiple deletions of redundant genes. It is probable that the rearrangements start and end at tRNA genes involved in the initial production of a tandemly duplicated gene region. Putative secondary structures for the 22 tRNAs and the origin of the L-strand replication (OL) are described. Evolutionary relationships were estimated from the concatenated sequences of the 12 proteins encoded in the H-strand of mtDNA among 37 vertebrate species. A quartet-puzzling tree showed that three amphibian species form a monophyletic clade and that the caecilian is a sister group of the monophyletic Anura.

Molecular synapomorphies resolve evolutionary relationships of extant jawed vertebrates

The evolutionary relationships of gnathostomes (jawed vertebrates), which comprise chondrichthyans (cartilaginous fishes), lobe-finned fishes (coelacanths and lungfishes), tetrapods, and actinopterygians (ray-finned fishes), have been debated for almost a century. Phylogenetic analyses based on fossils, morphology, and molecular sequences have generated different models of relationships that remain unresolved. We identified 13 derived shared molecular markers (synapomorphies) that define clades in the vertebrate lineage and used them to resolve the phylogenetic relationships of extant jawed vertebrates. Our markers include the presence or absence of insertions and deletions in coding sequences, nuclear introns, and alternatively spliced transcripts. The synapomorphies identified by us are congruent with each other and give rise to a single phylogenetic tree. This tree confirms that chondrichthyans are basal to all living gnathostomes, that lungfishes (Dipnoi) are the closest living relatives of tetrapods, and that bichirs (Cladistia) are the living members of the most ancient family of ray-finned fishes. Our study also provides molecular evidence to support the monophyly of living tetrapods and teleosts.

On the origin of and phylogenetic relationships among living amphibians

The phylogenetic relationships among the three orders of modern amphibians (Caudata, Gymnophiona, and Anura) have been estimated based on both morphological and molecular evidence. Most morphological and paleontological studies of living and fossil amphibians support the hypothesis that salamanders and frogs are sister lineages (the Batrachia hypothesis) and that caecilians are more distantly related. Previous interpretations of molecular data based on nuclear and mitochondrial rRNA sequences suggested that salamanders and caecilians are sister groups to the exclusion of frogs. In an attempt to resolve this apparent conflict, the complete mitochondrial genomes of a salamander (Mertensiella luschani) and a caecilian (Typhlonectes natans) were determined (16,656 and 17,005 bp, respectively) and compared with previously published sequences from a frog (Xenopus laevis) and several other groups of vertebrates. Phylogenetic analyses of the mitochondrial data supported with high bootstrap values the monophyly of living amphibians with respect to other living groups of tetrapods, and a sister group relationship of salamanders and frogs. The lack of phylogenetically informative sites in the previous rRNA data sets (because of its shorter size and higher among-site rate variation) likely explains the discrepancy between our results and those based on previous molecular data. Strong support of the Batrachia hypothesis from both molecule- and morphology-based studies provides a robust phylogenetic framework that will be helpful to comparative studies among the three living orders of amphibians and will permit better understanding of the considerably divergent vertebral, brain, and digit developmental patterns found in frogs and salamanders.