Studying sources of incongruence in arthropod molecular phylogenies: Sea spiders (Pycnogonida) as a case study

In this report, we analyze the phylogeny of Pycnogonida using the three nuclear and three mitochondrial markers currently sequenced for studying inter- and intrafamilial relationships within Arthropoda: 18S and 28S rRNA genes, Histone H3, cytochrome c oxidase subunit 1 (CO1), 12S and 16S rRNA genes. We identify several problems in previous studies, due to the use of inappropriate sequences (taxonomic misidentification, DNA contamination, sequencing errors, missing data) or taxa (outgroup choice). Our analyses show that most markers are not powerful to study the phylogeny of sea spiders. The results suggest however a recent diversification of the group (Mesozoic rather than Paleozoic) and the early divergence of Austrodecidae, followed by Colossendeidae, Pycnogonidae and Rhynchothoracidae. Except Ammotheidae and Callipallenidae, all other families were recovered as monophyletic. Analyses of synonymous sites in CO1 sequences reveal an extreme heterogeneity of nucleotide composition within sea spiders, as six unrelated species show a reverse strand-specific bias. We therefore suggest that several independent reversals of asymmetric mutational constraints occurred during the evolution of Pycnogonida, as a consequence of genomic inversions involving either the control region or a fragment containing the CO1 gene. These hypotheses are supported by the comparison of two complete mitochondrial genomes of sea spiders (Achelia bituberculata and Nymphon gracile) with that of Limulus.     

The mitochondrial genome sequence of the scorpion Centruroides limpidus (Karsch 1879) (Chelicerata; Arachnida)

The mitochondrial genome of the scorpion Centruroides limpidus (Chelicerata; Arachnida) has been completely sequenced and is 14519 bp long. The genome contains 13 protein-encoding genes, two ribosomal RNA genes, 21 transfer RNA genes and a large non-coding region related to the control region. The overall A + T composition is the lowest among the complete mitochondrial sequences published within the Chelicerata subphylum. Gene order and gene content differ slightly from that of Limulus polyphemus (Chelicerata: Xiphosura): i.e., the lack of the trnD gene, and the translocation–inversion of the trnI gene. Preliminary phylogenetic analysis of some Chelicerata shows that scorpions (C. limpidus and Mesobuthus gibbosus) make a tight cluster with the spiders (Arachnida; Araneae). Our analysis does not support that Scorpiones order is the sister group to all Arachnida Class, since it is closer to Araneae than to Acari orders.     

The evolution of mitochondrial genomes in modern frogs (Neobatrachia): nonadaptive evolution of mitochondrial genome reorganization

Background

Although mitochondrial (mt) gene order is highly conserved among vertebrates, widespread gene rearrangements occur in anurans, especially in neobatrachians. Protein coding genes in the mitogenome experience adaptive or purifying selection, yet the role that selection plays on genomic reorganization remains unclear. We sequence the mitogenomes of three species of Glandirana and hot spots of gene rearrangements of 20 frog species to investigate the diversity of mitogenomic reorganization in the Neobatrachia. By combing these data with other mitogenomes in GenBank, we evaluate if selective pressures or functional constraints act on mitogenomic reorganization in the Neobatrachia. We also look for correlations between tRNA positions and codon usage.

Results

Gene organization in Glandirana was typical of neobatrachian mitogenomes except for the presence of pseudogene trnS (AGY). Surveyed ranids largely exhibited gene arrangements typical of neobatrachian mtDNA although some gene rearrangements occurred. The correlation between codon usage and tRNA positions in neobatrachians was weak, and did not increase after identifying recurrent rearrangements as revealed by basal neobatrachians. Codon usage and tRNA positions were not significantly correlated when considering tRNA gene duplications or losses. Change in number of tRNA gene copies, which was driven by genomic reorganization, did not influence codon usage bias. Nucleotide substitution rates and d_N/d_S ratios were higher in neobatrachian mitogenomes than in archaeobatrachians, but the rates of mitogenomic reorganization and mt nucleotide diversity were not significantly correlated.

Conclusions

No evidence suggests that adaptive selection drove the reorganization of neobatrachian mitogenomes. In contrast, protein-coding genes that function in metabolism showed evidence for purifying selection, and some functional constraints appear to act on the organization of rRNA and tRNA genes. As important nonadaptive forces, genetic drift and mutation pressure may drive the fixation and evolution of mitogenomic reorganizations.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-691) contains supplementary material, which is available to authorized users.     

Complete mitochondrial DNA sequence of Aulopus japonicus (Teleostei: Aulopiformes), a basal Eurypterygii: longer DNA sequences and higher-level relationships

For the first step toward resolution of the higher-level relationships of the order Aulopiformes (Teleostei: Eurypterygii) using longer DNA sequences, we determined the complete mitochondrial DNA sequence for Aulopus japonicus (Aulopodidae). The entire genome was purified by gene amplification using a long PCR technique, and the products were subsequently used as templates for PCR with 63 fish-versatile and 3 species-specific primers that amplify contiguous, overlapping segments of the entire genome. Direct sequencing of the PCR products demonstrated that the genome (16 653 base pairs [bp]) contained the same 37 mitochondrial genes (2 ribosomal RNA, 22 transfer RNA, and 13 protein-coding genes) as found in other vertebrates, with the gene order identical to that in typical vertebrates. Maximum-parsimony analysis using nucleotide sequences from the concatenated 12 protein-coding genes (no third codon positions and excluding the ND6 gene) plus 22 tRNA genes (stem regions only) from eight teleosts placed A. japonicus in a reasonable phylogenetic position; those from individual protein-coding genes and the concatenated 22 tRNA genes alone, however, did not reproduce the expected phylogeny with few exceptions, probably owing to insufficient phylogenetic information in these smaller data sets. This result suggests that further taxonomic sampling and sequencing efforts may clarify limits and intra- and interrelationships of this morphologically and ecologically diverse group of fishes using mitochondrial genomic (mitogenomic) data. Received: August 31, 2000 / Revised: December 20, 2000 / Accepted: January 23, 2001     

The complete mitochondrial genome of the onychophoran Epiperipatus biolleyi reveals a unique transfer RNA set and provides further support for the ecdysozoa hypothesis

Onychophora (velvet worms) play a crucial role in current discussions on position of arthropods. The ongoing Articulata/Ecdysozoa debate is in need of additional ground pattern characters for Panarthropoda (Arthropoda, Tardigrada, and Onychophora). Hence, Onychophora is an important outgroup taxon in resolving the relationships among arthropods, irrespective of whether morphological or molecular data are used. To date, there has been a noticeable lack of mitochondrial genome data from onychophorans. Here, we present the first complete mitochondrial genome sequence of an onychophoran, Epiperipatus biolleyi (Peripatidae), which shows several characteristic features. Specifically, the gene order is considerably different from that in other arthropods and other bilaterians. In addition, there is a lack of 9 tRNA genes usually present in bilaterian mitochondrial genomes. All these missing tRNAs have anticodon sequences corresponding to 4-fold degenerate codons, whereas the persisting 13 tRNAs all have anticodons pairing with 2-fold degenerate codons. Sequence-based phylogenetic analysis of the mitochondrial protein-coding genes provides a robust support for a clade consisting of Onychophora, Priapulida, and Arthropoda, which confirms the Ecdysozoa hypothesis. However, resolution of the internal ecdysozoan relationships suffers from a cluster of long-branching taxa (including Nematoda and Platyhelminthes) and a lack of data from Tardigrada and further nemathelminth taxa in addition to nematodes and priapulids.     

Hox genes in sea spiders (Pycnogonida) and the homology of arthropod head segments

The pycnogonids (or sea spiders) are an enigmatic group of arthropods, classified in recent phylogenies as a sister-group of either euchelicerates (horseshoe crabs and arachnids), or all other extant arthropods. Because of their bizarre morpho-anatomy, homologies with other arthropod taxa have been difficult to assess. We review the main morphology-based hypotheses of correspondence between anterior segments of pycnogonids, arachnids and mandibulates. In an attempt to provide new relevant data to these controversial issues, we performed a PCR survey of Hox genes in two pycnogonid species, Endeis spinosa and Nymphon gracile, from which we could recover nine and six Hox genes, respectively. Phylogenetic analyses allowed to identify their orthology relationships. The Deformed gene from E. spinosa and the abdominal-A gene from N. gracile exhibit unusual sequence divergence in their homeodomains, which, in the latter case, may be correlated with the extreme reduction of the posterior region in pycnogonids. Expression patterns of two Hox genes (labial and Deformed) in the E. spinosa protonymphon larva are discussed. The anterior boundaries of their expression domains favour homology between sea spider chelifores, euchelicerates chelicerae and mandibulate (first) antennae, in contradistinction with previously proposed alternative schemes such as the protocerebral identity of sea spider chelifores or the absence of a deutocerebrum in chelicerates. In addition, while anatomical and embryological evidences suggest the possibility that the ovigers of sea spiders could be a duplicated pair of pedipalps, the Hox data support them as modified anterior walking legs, consistent with the classical views.Supplementary material is available for this article at and is accessible for authorized users.Guest editors Jean Deutsch and Gerhard Scholtz     

Pleistocene diversification of living squirrel monkeys (Saimiri spp.) inferred from complete mitochondrial genome sequences

In order to enhance our understanding of the evolutionary history of squirrel monkeys (Saimiri spp.), we newly sequenced and analyzed data from seven complete mitochondrial genomes representing six squirrel monkey taxa. While previous studies have lent insights into the taxonomy and phylogeny of the genus, phylogenetic relationships and divergence date estimates among major squirrel monkey clades remain unclear. Using maximum likelihood and Bayesian procedures, we inferred a highly resolved phylogenetic tree with strong support for a sister relationship between Saimiri boliviensis and all other Saimiri, for monophyly of Saimiri oerstedii and Saimiri sciureus sciureus, and for Saimiri sciureus macrodon as the sister lineage to the S. oerstedii/S. s. sciureus clade. We inferred that crown lineages for extant squirrel monkeys diverged around 1.5 million years ago (MYA) in the Pleistocene Epoch, with other major clades diverging between 0.9 and 1.1 MYA. Our results suggest a relatively recent timeline of squirrel monkey evolution and challenge previous conceptions about the diversification of the genus and its expansion into Central America.     

The compact mitochondrial genome of Zorotypus medoensis provides insights into phylogenetic position of Zoraptera

Background

Zoraptera, generally regarded as a member of Polyneoptera, represents one of the most enigmatic insect orders. Although phylogenetic analyses based on a wide array of morphological and/or nuclear data have been performed, the position of Zoraptera is still under debate. Mitochondrial genome (mitogenome) information is commonly considered to be preferable to reconstruct phylogenetic relationships, but no efforts have been made to incorporate it in Zorapteran phylogeny. To characterize Zoraptera mitogenome features and provide insights into its phylogenetic placement, here we sequenced, for the first time, one complete mitogenome of Zoraptera and reconstructed the phylogeny of Polyneoptera.

Results

The mitogenome of Zorotypus medoensis with an A + T content of 72.50% is composed of 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and a noncoding A + T-rich region. The gene content and arrangement are identical to those considered ancestral for insects. This mitogenome shows a number of very unusual features. First, it is very compact, comprising 14,572 bp, and is the smallest among all known polyneopteran mitogenomes. Second, both noncoding sequences and coding genes exhibit a significant decrease in size compared with those of other polyneopterans. Third, Z. medoensis mitogenome has experienced an accelerated substitution rate. Fourth, truncated secondary structures of tRNA genes occur with loss of dihydrouridine (DHU) arm in trnC, trnR, and trnS(AGN) and loss of TΨC arm in trnH and trnT. The phylogenetic analyses based on the mitogenome sequence information indicate that Zoraptera, represented by Z. medoensis, is recovered as sister to Embioptera. However, both Zoraptera and Embioptera exhibit very long branches in phylogenetic trees.

Conclusions

Characterization of Z. medoensis mitogenome contributes to our understanding of the enigmatic Zoraptera. Mitogenome data demonstrate an overall strong resolution of deep-level phylogenies of Polyneoptera but not Insecta. It is preferable to expand taxon sampling of Zoraptera and other poorly represented orders in future to break up long branches.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1156) contains supplementary material, which is available to authorized users.     

The complete mitochondrial genome of the yeast Pichia sorbitophila

Here, we report the complete nucleotide sequence of the 39 107-bp mitochondrial genome of the yeast Pichia sorbitophila . This genome is closely related to those of Candida parapsilosis and Debaryomyces hansenii , as judged from sequence similarities and synteny conservation. It encodes three subunits of cytochrome oxidase ( COX1, COX2 and COX3 ), three subunits of ATP synthase ( ATP6, ATP8 and ATP9 ), the seven subunits of NADH dehydrogenase ( NAD1-6 and NAD4L ), the apocytochrome b ( COB ), the large and small rRNAs and a complete set of tRNAs. Although the mitochondrial genome of P. sorbitophila contains the same core of mitochondrial genes observed in the ascomycetous yeasts, those coding for the RNAse P and the ribosomal protein VAR1p are missing. Moreover, the mtDNA of P. sorbitophila contains several introns in its genes and has the particularity of possessing an intron, which is not linked to any upstream exon.     

Making the most of mitochondrial genomes--markers for phylogeny, molecular ecology and barcodes in Schistosoma (Platyhelminthes: Digenea)

An increasing number of complete sequences of mitochondrial (mt) genomes provides the opportunity to optimise the choice of molecular markers for phylogenetic and ecological studies. This is particularly the case where mt genomes from closely related taxa have been sequenced; e.g., within Schistosoma. These blood flukes include species that are the causative agents of schistosomiasis, where there has been a need to optimise markers for species and strain recognition. For many phylogenetic and population genetic studies, the choice of nucleotide sequences depends primarily on suitable PCR primers. Complete mt genomes allow individual gene or other mt markers to be assessed relative to one another for potential information content, prior to broad-scale sampling. We assess the phylogenetic utility of individual genes and identify regions that contain the greatest interspecific variation for molecular ecological and diagnostic markers. We show that variable characters are not randomly distributed along the genome and there is a positive correlation between polymorphism and divergence. The mt genomes of African and Asian schistosomes were compared with the available intraspecific dataset of Schistosoma mansoni through sliding window analyses, in order to assess whether the observed polymorphism was at a level predicted from interspecific comparisons. We found a positive correlation except for the two genes (cox1 and nad1) adjoining the putative control region in S. mansoni. The genes nad1, nad4, nad5, cox1 and cox3 resolved phylogenies that were consistent with a benchmark phylogeny and in general, longer genes performed better in phylogenetic reconstruction. Considering the information content of entire mt genome sequences, partial cox1 would not be the ideal marker for either species identification (barcoding) or population studies with Schistosoma species. Instead, we suggest the use of cox3 and nad5 for both phylogenetic and population studies. Five primer pairs designed against Schistosoma mekongi and Schistosoma malayensis were tested successfully against Schistosoma japonicum. In combination, these fragments encompass 20-27% of the variation amongst the genomes (average total length approximately 14,000bp), thus providing an efficient means of encapsulating the greatest amount of variation within the shortest sequence. Comparative mitogenomics provides the basis of a rational approach to molecular marker selection and optimisation.     

Cloning and characterization of the platypus mitochondrial genome

The vertebrate mitochondrial genome is highly conserved in size and gene content. Among the chordates there appears to be one basic gene arrangement, but rearrangements in the mitochondrial gene order of the avian lineages have indicated that the mitochondrial genome may be more variable than once thought. Different gene orders in marsupials and eutherian mammals leave the ancestral mammalian order in some doubt. We have investigated the mitochondrial gene order in the platypus (Ornithorhynchus anatinus), a representative of the third major group of mammals, to determine which mitochondrial gene arrangement is ancestral in mammals. We have found that the platypus mtDNA conforms to the basic chordate gene arrangement, common to fish, amphibians, and eutherian mammals, indicating that this arrangement was the original mammalian arrangement, and that the unusual rearrangements observed in the avians and marsupials are probably lineage-specific. Correspondence to: N.J. Gemmell     

Complete mitochondrial genome of Periplaneta brunnea (Blattodea: Blattidae) and phylogenetic analyses within Blattodea

The complete mitochondrial genome (mitogenome) of Periplaneta brunnea was sequenced in this study and used to reconstruct the phylogenetic relationship of Blattodea. The circular mitogenome was 15,604?bp long and exhibited typical gene organization and order, consistent with other sequenced Periplaneta mitogenomes. The initiation codon of the P. brunnea COX1 gene was unusual in that no typical ATN or TTG start codon was found. The two longest intergenic spacer sequences found in the P. brunnea mitogenome were 21 and 17?bp long. Twenty-one base spacer had a 4?bp motif (TATT) between tRNA-Glu and tRNA-Met that conservatively displayed in 9 sequenced blattarian mitogenomes. The second spacer was between tRNA-Ser (UCN) and NAD1 containing a 7?bp motif (WACTTAA) that was highly conserved in 14 blattarian mitogenomes. The control region showed a relatively fixed motif present in 6 Blattidae mitogenomes, with a big stem-loop structure. Phylogenetic analyses were conducted using site-homogeneous models based on 13 protein-coding genes (PCGs) and two RNA genes. The trees derived from Bayesian inference and maximum likelihood analyses and recovered a relatively stable relationship among major lineages except for the position of Polyphagidae and inter-family relationships of Blaberidae. Analyses supported the monophyly of Blattidae, Blaberidae, Blattellidae, Polyphagidae, Dictyoptera, and the paraphyly of Blattaria. We also found Mantodea was the sister clade to (Blattaria?+?Isoptera), being the basal position of Dictyoptera in all topologies. Meanwhile, our results also consistently supported that Isoptera should be clustered with Blattaria of Blattodea.     

The complete mitochondrial genome of the Antarctic sea spider Ammothea carolinensis (Chelicerata; Pycnogonida)

Mitochondria are responsible for the oxidative phosphorylation process. Accordingly, putatively adaptive changes in their genomic features have been variously associated with major eco-physiological shifts in animal evolution, including increased metabolic rates and heat adaptation. Antarctic pycnogonids offer an interesting system to test whether the selective pressure for heat production and increased aerobic metabolism may be driving genomic changes like: (a) unusual compositional biases at the nucleotide and amino acid level, possibly related to cold adaptation; (b) an accelerated rate of mutations/genomic rearrangements, possibly related to the mutagenic effects of oxygen intermediates. The complete mitochondrial genome (mtDNA) of the Antarctic sea spider Ammothea carolinensis Leach, 1814 (Arthropoda: Pycnogonida), the type species for the genus Ammothea, has been determined and is here compared to known genomes from Antarctic and temperate species. We describe a marked heterogeneity in base composition skewness parameters as well as a strong signature of purifying selection toward an increase in thymines at second codon positions, possibly associated with an increased stability of hydrophobic inter-membrane domains. We further observe a fairly high rate of genomic changes, including a possible hot spot of recombination at the level of tRNA-Q. Nevertheless, these features do not seem to be restricted to the two Antarctic pycnogonids analyzed, as to suggest a causal relationship between cold adaptation and genomic changes, and are better interpreted as basal features shared by the entire group. The relevance of the newly determined sequence for the phylogeny of pycnogonids, including its base composition and genomic rearrangements, is further discussed.     

The mitochondrial genome of the spinyhead croaker Collichthys lucida: genome organization and phylogenetic consideration

The complete mitochondrial genome of the spiny head croaker Collichthys lucida was determined in the present study. The mitochondrial DNA was 16,442 base pairs in length, and contained 13 protein coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and one major non-coding control region, with the content and order of genes being similar to those in typical teleosts. Most of the genes of C. lucida were encoded on the H-strand, while the ND6 and eight tRNA (Gln, Ala, Asn, Cys, Tyr, Ser (UCN), Glu and Pro) genes were encoded on the L-strand. The reading frames of two pairs of genes overlapped: ATPase 8 and 6 and ND4L and ND4 by ten and seven nucleotides, respectively. The control region was unusually short at only 768bp, and absence of typical conserved blocks (CSB-D, CSB-E, and CSB-F). Phylogenetic analyses indicated that C. lucida was located in the cluster of fish species from the family Sciaenidae, supporting the traditional taxonomic classification of fish, and in the cluster of Serranidae, the divergence time in Plectropomus leopardus is longer than that among its coordinal species. On the other hand, phylogenetic analyses do not support the monophyletic of family Centracanthidae and genera Larimichthys and Collichthys, which is against the morphological results.     

Variation in Nymphon brevirostre (Hodge), and the status of N. rubrum (Hodge) (Arthropoda: Pycnogonida)

A series of morphometric analyses conducted on specimens of Nymphon brevirostre Hodge and Nymphon rubrum Hodge from British waters has shown that the characteristics variously used in the past to distinguish these two species are invalid. In this material these characteristics range from the extreme of the typical N. brevirostre form to that of the typical N. rubrum form, with most of the specimens exhibiting an intermediate morphology. Multivariate analyses showed no significant splitting of the material into two groups which might correlate to Hodge's two species. Protonymphon larvae from males of both typical extreme forms were morphologically identical. It is concluded that N. rubrum is a junior synonym of N. brevirostre, a variable species whose morphology tends towards the rubrum form as the animal grows.     

Variation in Nymphon brevirostre (Hodge), and the status of N. rubrum (Hodge) (Arthropoda: Pycnogonida)

A series of morphometric analyses conducted on specimens of Nymphon brevirostre Hodge and Nymphon rubrum Hodge from British waters has shown that the characteristics variously used in the past to distinguish these two species are invalid. In this material these characteristics range from the extreme of the typical N. brevirostre form to that of the typical N. rubrum form, with most of the specimens exhibiting an intermediate morphology. Multivariate analyses showed no significant splitting of the material into two groups which might correlate to Hodge's two species. Protonymphon larvae from males of both typical extreme forms were morphologically identical. It is concluded that N. rubrum is a junior synonym of N. brevirostre , a variable species whose morphology tends towards the rubrum form as the animal grows.     

尼罗鳄线粒体基因组全序列分析及鳄类系统发生关系的探讨

大多数脊椎动物的线粒体基因组（约16—18kb）的组成是相对较稳定的,但在不同类群中,线粒体基因组在基因结构和基因排列方式等方面均显示了极大的多样性,这种多样性可能反映了真核细胞不同的进化路线（Saccone et al．,1999）。就目前的研究而言,线粒体基因组是惟一一个能够从基因组水平上来分析动物系统发生的分子标记,可以从线粒体基因组序列信息、基因组成及基因排列方式等进行多方位的分子进化研究,因而线粒体基因组全序列将成为动物分子系统发生最有力的证据（Saccone et al．,1999）。     

The partial mitochondrial genome of the Cephalothrix rufifrons (Nemertea, Palaeonemertea): characterization and implications for the phylogenetic position of Nemertea

A continuous 10.1kb fragment of the Cephalothrix rufifrons (Nemertea, Palaeonemertea) mitochondrial genome was sequenced and characterized to further assess organization of protostome mitochondrial genomes and evaluate the phylogenetic potential of gene arrangement and amino acid characters. The genome is A-T rich (72%), and this biased base composition is partly reflected in codon usage. Inferred tRNA secondary structures are typical of those reported for other metazoan mitochondrial DNAs. The arrangement of the 26 genes contained in the fragment exhibits marked similarity to those of many protostome taxa, most notably molluscs with highly conserved arrangements and a phoronid. Separate and simultaneous phylogenetic analyses of inferred amino acid sequences and gene adjacencies place the nemertean within the protostomes among coelomate lophotrochozoan taxa, but do not find a well-supported sister taxon link.