Genetic divergence and evolutionary relationships in six species of genera Hoplobatrachus and Euphlyctis (Amphibia: Anura) from Bangladesh and other Asian countries revealed by mitochondrial gene sequences

To elucidate the species composition, genetic divergence, evolutionary relationships, and divergence time of Hoplobatrachus and Euphlyctis frogs (subfamily Dicroglossinae, family Ranidae) in Bangladesh and other Asian countries, we analyzed the mitochondrial Cyt b, 12S, and 16S rRNA genes of 252 specimens. Our phylogenetic analyses showed 13 major clades corresponding to several cryptic species as well as to nominal species in the two genera. The results suggested monophyly of Asian Hoplobatrachus species, but the position of African Hoplobatrachus occipitalis was not clarified. Nucleotide divergence and phylogenetic data suggested the presence of allopatric cryptic species allied to Euphlyctis hexadactylus in Sundarban, Bangladesh and several parapatric cryptic species in the Western Ghats, India. The presence of at least two allopatric cryptic species among diverged Euphlyctis cyanophlyctis in Bangladesh, India, and Sri Lanka was also suggested. In some cases, our estimated divergence times matched the paleogeological events of South and Southeast Asian regions that may have led to the divergence of Hoplobatrachus and Euphlyctis taxa. Especially, land formation at Bangladesh (15–10 Ma) may have allowed the spread of these frog taxa to Southeast Asian areas, and the aridification of central India (5.1–1.6 Ma) might have affected the gene flow of widely distributed species. The present study revealed prior underestimation of the richness of the amphibian fauna in this region, indicating the possible occurrence of many cryptic species among these groups.     

虎纹蛙及其它几种蛙科动物线粒体16S rRNA的空间结构分析

测定了虎纹蛙16S rRNA基因全序列,并从GeneBank上下载了5种蛙科动物同源序列,利用CLUSTAL与MEGA3.1软件比较了种间的序列差异,用RNAdraw软件推测分析二级结构的空间构象,并根据空间结构对这6种蛙科动物的分类地位进行讨论分析。结果表明：虎纹蛙与黑斑蛙的二级结构最为相似,与泽蛙结构大体一致,而与黑斑蛙同属的金线蛙则表现出较大的差异;凹耳蛙与上述几种差别较大;福建大头娃与这几种蛙科动物差距最为明显。空间结构分析与序列比较结果存在一定的差异。     

The complete mitochondrial genomes of two schizothoracine fishes (Teleostei,Cypriniformes): A novel minisatellite in fish mitochondrial genomes

Minisatellites, a class of variable number tandem repeats (VNTRs), are abundant throughout the control region in animal mitochondrial DNA (mtDNA) but rare in other regions of animal mtDNA. Here, we reported a novel minisatellite in fish mitochondrial genomes. We first determined the complete mitochondrial genomes of two schizothoracine fishes (Herzensteinia microcephalus and Schizopygopsis pylzovi) and found a type of minisatellites in a novel region between the tRNA‐Thr and tRNA‐Pro genes in their mtDNA. To explore the origin and evolution of the minisatellites in different schizothoracine and closely related fishes, we analyzed the available 80 fish mitogenomes which represent five closely related tribes of cyprinine fishes. The results from the phylogenetic analyses show that the schizothoracine fishes sensu stricto is not a monophyletic group and is divided into two clades (Schizothoracini and Schizopygopsini); and the minisatellite is only present in Schizopygopsini distributed in the region between the two tRNA genes (tRNA‐Thr and tRNA‐Pro) of the mtDNA. This is the first record of a minisatellite in a non‐control region of fish mitogenome.     

Multiple rearrangements in mitochondrial genomes of Isopoda and phylogenetic implications

In this study, we analyse the evolutionary dynamics and phylogenetic implications of gene order rearrangements in five newly sequenced mitochondrial (mt) genomes and four published mt genomes of isopod crustaceans. The sequence coverage is nearly complete for four of the five newly sequenced species, with only the control region and some tRNA genes missing, while in Janira maculosa only two thirds of the genome could be determined. Mitochondrial gene order in isopods seems to be more plastic than that in other crustacean lineages, making all nine known mt gene orders different. Especially the asellote Janira is characterized by many autapomorphies. The following inferred ancestral isopod mt gene order exists slightly modified in modern isopods: nad1, tnrL1, rrnS, control region, trnS1, cob, trnT, nad5, trnF. We consider the inferred gene translocation events leading to gene rearrangements as valuable characters in phylogenetic analyses. In this first study covering major isopod lineages, potential apomorphies were identified, e.g., a shared relative position of trnR in Valvifera. We also report one of the first findings of homoplasy in mitochondrial gene order, namely a shared relative position of trnV in unrelated isopod lineages. In addition to increased taxon sampling secondary structure, modification in tRNAs and GC-skew inversion may be potentially fruitful subjects for future mt genome studies in a phylogenetic context.     

New features of Asian Crassostrea oyster mitochondrial genomes: A novel alloacceptor tRNA gene recruitment and two novel ORFs

A feasible way to perform evolutionary analyses is to compare characters divergent enough to observe significant differences, but sufficiently similar to exclude saturation of the differences that occurred. Thus, comparisons of invertebrate mitochondrial (mt) genomes at low taxonomic levels can be extremely helpful in investigating patterns of variation and evolutionary dynamics of genomes, as intermediate stages of the process may be identified. Fortunately, in this study, we newly sequenced the mt genome of the eighth member of Asian Crassostrea oysters which can provide necessary intermediate characters for us to believe that the variation of Crassostrea mt genomes is considerably greater than previously acknowledged. Several new features of Asian Crassostrea oyster mitochondrial genomes were revealed, and our results are particularly significant as they 1) suggest a novel model of alloacceptor tRNA gene recruitment, namely "vertical" tRNA gene recruitment, which can be successfully used to explain the origination of the unusually additional trnK and trnQ genes (annotated as trnK(2) and trnQ(2) respectively) in the mt genomes of the five Asian oysters, and we speculate that this recruitment progress may be a common phenomenon in the evolution of the tRNA multigene family; 2) reveal the existence of two additional, lineage-specific, mtDNA-encoded genes that may originate from duplication of nad2 followed by rapid evolutionary change. Each of these two genes encodes a unique amino terminal signal peptide, thus each might possess an unknown function; and 3) identify for the first time the atp8 gene in oysters. The present study thus gives further credence to the comparison of congeneric bivalves as a meaningful strategy to investigate mt genomic evolutionary trends in genome organization, tRNA multigene family, and gene loss and/or duplication that are difficult to undertake at higher taxonomic levels. In particular, our study provides new evidence for the identification and characterization of ORFs in the "non-coding region" of animal mt genomes.     

A New Species of Euphlyctis (Anura: Dicroglossidae) from Barisal,Bangladesh

A new species of the genus Euphlyctis is described from the Barisal district of Bangladesh and compared with its morphologically similar and geographically proximate congeners. The new species is highly divergent in comparison to other congeneric species on basis of sequence divergence in mitochondrial DNA gene sequences (ranging from 5.5% to 17.8% divergence). Euphlyctis kalasgramensis sp. nov. can be readily diagnosed by having the following combination of characters: snout-vent length (SVL) 30.44 – 37.88 mm, absence of mid-dorsal line, nostril–snout length 3% of SVL, nostril much closer to snout tip than eye, nostril–snout length 48% of distance from front of eyes to nostril, relative length of fingers (shortest to longest: 1 = 2 < 4 < 3), tibia length 59% of SVL, foot length 55% of SVL.     

Complete sequences of mitochondrial genomes from the Baltic mussel Mytilus trossulus

Marine mussels Mytilus possess two mitochondrial (mt) genomes, which undergo doubly uniparental inheritance (DUI). Female (F) and male (M) genomes are usually highly diverged at the sequence level. Both genomes contain the same set of metazoan genes (for 12 proteins, 2 rRNAs and 23 tRNAs), both lack the atp8 gene and have two tRNAs for methionine. However, recently recombination between those variants has been reported. Both original F and M mt genomes of M. trossulus were replaced by M. edulis mtDNA in the Baltic populations. Highly diverged M genome occurs rarely in the Baltic mussels. Full sequences of the M genome identified in males (sperm) and F genome in females (eggs) were obtained. Both genomes were diverged by 24% in nucleotide sequence, but had similar nucleotide composition and codon usage bias. Constant domain (CD) of the control region (CR), the tRNA and rRNA genes were the most conserved. The most diverged was the variable domain 1 (VD1) of the control region. The F genome was longer than M by 147 bp. and the main difference was localised in the VD1 region. No recombination was observed in whole mtDNA of both studied variants. Nuclear mitochondrial pseudogenes (numts) have not been found by hybridisation with probes complementary to several fragments of the Baltic M. trossulus mtDNA.     

Bird mitochondrial gene order: insight from 3 warbler mitochondrial genomes

Two main gene orders exist in birds: the ancestral gene order and the remnant control region (CR) 2 gene order. These gene orders differ by the presence of 1 or 2 copies of the CR, respectively. Among songbirds, Oscines were thought to follow the ancestral gene order, with the exception of the lyrebird and Phylloscopus warblers. Here, we determined the complete mitochondrial genome sequence of 3 non-Phylloscopus warblers species and found that the blackcap (Sylvia atricapilla) and the reed warbler (Acrocephalus scirpaceus) have 2 almost identical copies of the CR, whereas the eastern orphean warbler (Sylvia crassirostris) follows the remnant CR 2 gene order. Our results contradict previous studies suggesting that Acrocephalus and most sylvioid warblers exhibit the ancestral gene order. We were able to trace this contradiction to a misidentification of gene order from polymerase chain reaction length determination. We thus suggest that passerine gene order evolution needs to be revised.     

Massive gene rearrangements of mitochondrial genomes and implications for the phylogeny of Trichoptera (Insecta)

Mitochondrial genomes have been widely used for phylogenetic reconstruction and evolutionary analysis in various groups of Insecta. Gene rearrangements in the mitogenome can be informative characters for phylogenetic reconstruction and adaptive evolution. Trichoptera is one of the most important groups of aquatic insects. Prior to this study, complete mitogenomes from Trichoptera were restricted to eight families, resulting in a biased view of their mitogenome structure and evolution. Here, we assemble new mitogenomes for 66 species by high-throughput sequencing. The mitogenomes of 19 families and 47 genera are documented for the first time. Combined with 16 previously published mitogenomes of Trichoptera, we find 14 kinds of gene rearrangement patterns novel for Trichoptera, including rearrangement of protein-coding genes, tRNAs and control regions. Simultaneously, we provide evidence for the occurrence of tandem duplication and non-random loss events in the mitogenomes of three families. Phylogenetic analyses show that Hydroptilidae was recovered as a sister group to Annulipalpia. The increased nucleotide substitution rate and adaptive evolution may have affected the mitochondrial gene rearrangements in Trichoptera. Our study offers new insights into the mechanisms and patterns of mitogenome rearrangements in Insecta at large and into the usefulness of mitogenomic gene order as a phylogenetic marker within Trichoptera.     

Gene rearrangements in gekkonid mitochondrial genomes with shuffling,loss, and reassignment of tRNA genes

Background

Vertebrate mitochondrial genomes (mitogenomes) are 16–18 kbp double-stranded circular DNAs that encode a set of 37 genes. The arrangement of these genes and the major noncoding region is relatively conserved through evolution although gene rearrangements have been described for diverse lineages. The tandem duplication-random loss model has been invoked to explain the mechanisms of most mitochondrial gene rearrangements. Previously reported mitogenomic sequences for geckos rarely included gene rearrangements, which we explore in the present study.

Results

We determined seven new mitogenomic sequences from Gekkonidae using a high-throughput sequencing method. The Tropiocolotes tripolitanus mitogenome involves a tandem duplication of the gene block: tRNA^Arg, NADH dehydrogenase subunit 4L, and NADH dehydrogenase subunit 4. One of the duplicate copies for each protein-coding gene may be pseudogenized. A duplicate copy of the tRNA^Arg gene appears to have been converted to a tRNA^Gln gene by a C to T base substitution at the second anticodon position, although this gene may not be fully functional in protein synthesis. The Stenodactylus petrii mitogenome includes several tandem duplications of tRNA^Leu genes, as well as a translocation of the tRNA^Ala gene and a putative origin of light-strand replication within a tRNA gene cluster. Finally, the Uroplatus fimbriatus and U. ebenaui mitogenomes feature the apparent loss of the tRNA^Glu gene from its original position. Uroplatus fimbriatus appears to retain a translocated tRNA^Glu gene adjacent to the 5’ end of the major noncoding region.

Conclusions

The present study describes several new mitochondrial gene rearrangements from Gekkonidae. The loss and reassignment of tRNA genes is not very common in vertebrate mitogenomes and our findings raise new questions as to how missing tRNAs are supplied and if the reassigned tRNA gene is fully functional. These new examples of mitochondrial gene rearrangements in geckos should broaden our understanding of the evolution of mitochondrial gene arrangements.

Electronic supplementary material

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

Complete mitochondrial genomes of two Japanese precious corals, Paracorallium japonicum and Corallium konojoi (Cnidaria, Octocorallia, Coralliidae): notable differences in gene arrangement

Precious coral are taxonomically a group of corals that belong to the family Coralliidae within the order Alcyonacea, subclass Octocorallia, and class Anthozoa, whose skeletal axes are used for jewelry. They are distributed in the Mediterranean Sea and in waters adjacent to Japan, Taiwan, Midway Island and the Hawaiian Islands. The genus Corallium of the family Coralliidae was recently divided into two genera, Corallium and Paracorallium, based on morphological observations, but insufficient molecular evidence to support this classification has been presented to date. We determined for the first time the complete mitochondrial genome sequence of two precious corals P. japonicum and C. konojoi, in order to clarify their systematic positions. The circular mitochondrial genomes of P. japonicum and C. konojoi are 18,913bp and 18,969bp in length, respectively, and encode 13 typical energy pathway protein coding genes (nad1-6, nad4L, cox1-3, cob, atp6 and atp8), two ribosomal RNA genes (rns and rnl), a transfer RNA (trnM) and a mismatch repair gene homologue msh1. The two genomes have an overall nucleotide sequence identity of 97.5%, which is comparable to that between Acanella eburnea and Keratoisidinae sp. belonging to Octocorallia. Surprisingly, however, their gene arrangements were not identical. Phylogenetic analyses using seven complete mitochondrial genome sequences belonging to species in the subclass Octocorallia indicated that within the subclass, at least three gene order rearrangement events occurred during evolution. Our results support the validity of the morphological classification that separated the family Coralliidae into two genera, Corallium and Paracorallium.     

An insight into molecular taxonomy of bufonids,microhylids, and dicroglossid frogs: First genetic records from Pakistan

The current study was focused on documentation of amphibian assemblage in North Punjab and Islamabad Capital Territory, Pakistan, by using mitochondrial gene sequences of 16S rRNA. Our study entailed 37% of the known amphibian species of the country. We provided a phylogenetic analysis based on 74 newly generated mitochondrial 16S rRNAs from nine species of genus Microlyla, Duttaphrynus, Allopaa, Nanorana, Sphaerotheca, Minervarya, Hoplobatrachus, and Euphlyctis. We employed the maximum‐likelihood inference and Bayesian analysis to assess the taxonomic status of the samples obtained from Pakistan, with respect to other congeneric species from surrounding regions. Our findings confirmed the taxonomic status of South Asian anuran species Duttaphrynus stomaticus, Duttaphrynus melanostictus, Microhyla nilphamariensis, Allopaa hazarensis, Nanorana vicina, Sphaerotheca maskeyi (synonym: S. pashchima), Minervarya pierrei, Hoplobatrachus tigerinus, and Euphlyctis kalasgramensis in Pakistan. We have reported new country records of genus Minervarya ( M. pierrei). Minervarya pierrei was previously misidentified as Fejervarya limnocharis, due to dearth of genetic information. We provided the first genetic records of our endemic species N. vicina. The results revealed the taxonomic placement of N. vicina with respect to its congeners and validated the taxonomic status of N. vicina from its type locality (Murree) for the first time. The findings of the present study also indicated the paraphyletic relationship of A.‐ hazarensis with Nanorana species. So, based on our phylogenetic inferences, morphological characters, and habitat preferences, validity of generic status of A. hazarensis is undecided. As our data were not enough to resolve this issue, we suggest sequencing of additional mitochondrial and nuclear genes in the future studies to get a better resolution. We recommend carrying out extensive surveys throughout the country for proper scientific documentation of amphibians of Pakistan. Many new species, some of them might be endemic to Pakistan, are expected to be discovered, and taxonomic status of other species would be resolved.     

The complete mitochondrial genomes of sixteen ardeid birds revealing the evolutionary process of the gene rearrangements

Background

The animal mitochondrial genome is generally considered to be under selection for both compactness and gene order conservation. As more mitochondrial genomes are sequenced, mitochondrial duplications and gene rearrangements have been frequently identified among diverse animal groups. Although several mechanisms of gene rearrangement have been proposed thus far, more observational evidence from major taxa is needed to validate specific mechanisms. In the current study, the complete mitochondrial DNA of sixteen bird species from the family Ardeidae was sequenced and the evolution of mitochondrial gene rearrangements was investigated. The mitochondrial genomes were then used to review the phylogenies of these ardeid birds.

Results

The complete mitochondrial genome sequences of the sixteen ardeid birds exhibited four distinct mitochondrial gene orders in which two of them, named as “duplicate tRNA^Glu–CR” and “duplicate tRNA^Thr–tRNA^Pro and CR”, were newly discovered. These gene rearrangements arose from an evolutionary process consistent with the tandem duplication - random loss model (TDRL). Additionally, duplications in these gene orders were near identical in nucleotide sequences within each individual, suggesting that they evolved in concert. Phylogenetic analyses of the sixteen ardeid species supported the idea that Ardea ibis, Ardea modesta and Ardea intermedia should be classified as genus Ardea, and Ixobrychus flavicollis as genus Ixobrychus, and indicated that within the subfamily Ardeinae, Nycticorax nycticorax is closely related to genus Egretta and that Ardeola bacchus and Butorides striatus are closely related to the genus Ardea.

Conclusions

The duplicate tRNAThr–CR gene order is found in most ardeid lineages, suggesting this gene order is the ancestral pattern within these birds and persisted in most lineages via concerted evolution. In two independent lineages, when the concerted evolution stopped in some subsections due to the accumulation of numerous substitutions and deletions, the duplicate tRNAThr–CR gene order was transformed into three other gene orders. The phylogenetic trees produced from concatenated rRNA and protein coding genes have high support values in most nodes, indicating that the mitochondrial genome sequences are promising markers for resolving the phylogenetic issues of ardeid birds when more taxa are added.

Electronic supplementary material

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

Numerous small rearrangements of gene content,order and orientation differentiate grass genomes

Comparative genetic mapping has indicated that the grass family (Poaceae) exhibits extensive chromosomal collinearity. In order to investigate microcollinearity in these genomes, several laboratories have begun to undertake comparative DNA sequence analyses of orthologous chromosome segments from various grass species. Five different regions have now been investigated in detail, with four regions sequenced for maize, rice and sorghum, plus two for wheat and one for barley. In all five of these segments, gene rearrangements were observed in at least one of the comparisons. Most of the detected rearrangements are small, involving the inversion, duplication, translocation or deletion of DNA segments that contain only 1-3 genes. Even closely related species, like barley and wheat or maize and sorghum, exhibit approximately 20% alterations in gene content or orientation. These results indicate that thousands of small genetic rearrangements have occurred in several grass lineages since their divergence from common ancestors. These rearrangements have largely been missed by genetic mapping and will both complicate and enrich the use of comparative genetics in the grasses.     

Complete mitochondrial genomes of the human follicle mites Demodex brevis and D. folliculorum: novel gene arrangement,truncated tRNA genes,and ancient divergence between species

Background

Follicle mites of the genus Demodex are found on a wide diversity of mammals, including humans; surprisingly little is known, however, about the evolution of this association. Additional sequence information promises to facilitate studies of Demodex variation within and between host species. Here we report the complete mitochondrial genome sequences of two species of Demodex known to live on humans—Demodex brevis and D. folliculorum—which are the first such genomes available for any member of the genus. We analyzed these sequences to gain insight into the evolution of mitochondrial genomes within the Acariformes. We also used relaxed molecular clock analyses, based on alignments of mitochondrial proteins, to estimate the time of divergence between these two species.

Results

Both Demodex genomes shared a novel gene order that differs substantially from the ancestral chelicerate pattern, with transfer RNA (tRNA) genes apparently having moved much more often than other genes. Mitochondrial tRNA genes of both species were unusually short, with most of them unable to encode tRNAs that could fold into the canonical cloverleaf structure; indeed, several examples lacked both D- and T-arms. Finally, the high level of sequence divergence observed between these species suggests that these two lineages last shared a common ancestor no more recently than about 87 mya.

Conclusions

Among Acariformes, rearrangements involving tRNA genes tend to occur much more often than those involving other genes. The truncated tRNA genes observed in both Demodex species would seem to require the evolution of extensive tRNA editing capabilities and/or coevolved interacting factors. The molecular machinery necessary for these unusual tRNAs to function might provide an avenue for developing treatments of skin disorders caused by Demodex. The deep divergence time estimated between these two species sets a lower bound on the time that Demodex have been coevolving with their mammalian hosts, and supports the hypothesis that there was an early split within the genus Demodex into species that dwell in different skin microhabitats.

Electronic supplementary material

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

Complete mitochondrial genomes of prasinophyte algae Pyramimonas parkeae and Cymbomonas tetramitiformis

Mitochondria are archetypal eukaryotic organelles that were acquired by endosymbiosis of an ancient species of alpha‐proteobacteria by the last eukaryotic common ancestor. The genetic information contained within the mitochondrial genome has been an important source of information for resolving relationships among eukaryotic taxa. In this study, we utilized mitochondrial and chloroplast genomes to explore relationships among prasinophytes. Prasinophytes are represented by diverse early‐diverging green algae whose physical structures and genomes have the potential to elucidate the traits of the last common ancestor of the Viridiplantae (or Chloroplastida). We constructed de novo mitochondrial genomes for two prasinophyte algal species, Pyramimonas parkeae and Cymbomonas tetramitiformis, representing the prasinophyte clade. Comparisons of genome structure and gene order between these species and to those of other prasinophytes revealed that the mitochondrial genomes of P. parkeae and C. tetramitiformis are more similar to each other than to other prasinophytes, consistent with other molecular inferences of the close relationship between these two species. Phylogenetic analyses using the inferred amino acid sequences of mitochondrial and chloroplast protein‐coding genes resolved a clade consisting of P. parkeae and C. tetramitiformis; and this group (representing the prasinophyte clade I) branched with the clade II, consistent with previous studies based on the use of nuclear gene markers.     

Complete mitochondrial genomes of Ceratobaeus sp. and Idris sp. (Hymenoptera: Scelionidae): shared gene rearrangements as potential phylogenetic markers at the tribal level

We sequenced the complete mitochondrial genomes of two sceliond taxa (Ceratobaeus sp. and Idris sp.). An atypical tRNA-Arg which lacks a D-stem was identified in both taxa, and represents a potentially derived character of sceliond wasps. A number of tRNA genes have rearranged in the two mitochondrial genomes compared with the ancestral organization. Some of these derived genome organizations are shared, and thus have much potential as phylogenetic markers at the tribal level in the subfamily Scelioninae. We test the influence of third codon inclusion/exclusion, alignment methods and partition schemes on the reconstruction of phylogenetic relationships. The results show that inclusion of third codon positions does not appear to be problematic when investigating the phylogeny of closely related taxa. Muscle and PartitionFinder schemes significantly improve the likelihood scores.     

Complete nucleotide sequences of mitochondrial genomes of two solitary entoprocts, Loxocorone allax and Loxosomella aloxiata: implications for lophotrochozoan phylogeny

The complete nucleotide sequences of the mitochondrial (mt) genomes of the entoprocts Loxocorone allax and Loxosomella aloxiata were determined. Both species carry the typical gene set of metazoan mt genomes and have similar organizations of their mt genes. However, they show differences in the positions of two tRNA(Leu) genes. Additionally, the tRNA(Val) gene, and half of the long non-coding region, is duplicated and inverted in the Loxos. aloxiata mt genome. The initiation codon of the Loxos. aloxiata cytochrome oxidase subunit I gene is expected to be ACG rather than AUG. The mt gene organizations in these two entoproct species most closely resemble those of mollusks such as Katharina tunicata and Octopus vulgaris, which have the most evolutionarily conserved mt gene organization reported to date in mollusks. Analyses of the mt gene organization in the lophotrochozoan phyla (Annelida, Brachiopoda, Echiura, Entoprocta, Mollusca, Nemertea, and Phoronida) suggested a close phylogenetic relationship between Brachiopoda, Annelida, and Echiura. However, Phoronida was excluded from this grouping. Molecular phylogenetic analyses based on the sequences of mt protein-coding genes suggested a possible close relationship between Entoprocta and Phoronida, and a close relationship among Brachiopoda, Annelida, and Echiura.