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.     

从细胞色素b基因序列探讨笛鲷属的分子系统发生关系

测定了9种中国南海的笛鲷属鱼类的细胞色素b基因的部分序列,结合来自GenBank中1种分布于菲律宾和9种分布于美国大西洋的笛鲷属鱼类的相应同源序列,用邻接法和最大简约法构建分子系统树。结果显示:红鳍笛鲷(Lutjanuserythropterus)与红笛鲷(L.sanguineus)之间的同源序列碱基差异百分率只有0.32%,支持二者是同种异名的观点;中国南海的笛鲷属鱼类间的平均碱基差异要高于美国大西洋笛鲷属鱼类。在MP和NJ树中,美国大西洋笛鲷表现为亲缘关系较近,来源于中国南海的笛鲷鱼类相对集中在树的基部,分歧较大。这与所研究的笛鲷地理分布和地理隔离基本相一致,同时也说明中国南海笛鲷分化较早并且分歧较大。     

Phylogenetic relationships among Pneumocystis from Asian macaques inferred from mitochondrial rRNA sequences

The presence of Pneumocystis organisms was detected by nested-PCR at mitochondrial large subunit (mtLSU) rRNA gene in 23 respiratory samples from Asian macaques representing two species: Macaca mulatta and M. fascicularis. A very high level of sequence heterogeneity was detected with 18 original sequence types. Two genetic groups of Pneumocystis could be distinguished from the samples. Within each group, the extent of genetic divergence was low (2.5+/-1.4% in group 1 and 2.3+/-1.7% in group 2). Genetic divergences were systematically higher when macaque-derived sequence types were compared with Pneumocystis mtLSU sequences from other primate species (from 5.3+/-2.7% to 19.3+/-3.0%). The two macaque-derived groups may be considered as distinct Pneumocystis species. Surprisingly, these Pneumocystis species were recovered from both M. mulatta and M. fascicularis suggesting that host-species restriction may not systematically occur in the genus Pneumocystis. Alternatively, these observations question about the species concept in macaques.     

Higher-level relationships of snakes inferred from four nuclear and mitochondrial genes

Higher-level snake relationships are inferred from sequence analyses of one nuclear gene (C-mos) and three mitochondrial genes (12S rRNA, 16S rRNA and cytochrome b). Extant snakes belong to two lineages: the fossorial Scolecophidia, which feed on small prey on a frequent basis, and the ecologically diverse Alethinophidia ('typical' snakes), which feed on large prey on an infrequent basis. The vast majority of Alethinophidia, if not all of them, belong to two clades, corresponding to two distinct prey neutralization modes: unimodal constriction for the Henophidia (locomotor and feeding systems coupled) and injection of toxic saliva, in addition (or not) to diverse alternate modes of constriction, for the Caenophidia (locomotor and feeding systems uncoupled). Within Alethinophidia, non-macrostomatan (small gape) Aniliidae (genus Anilius) and macrostomatan (large gape) Tropidophiidae (genera Trachyboa and Tropidophis), both from the Neotropics, are closest relatives. Although our data are insufficient to robustly infer the ancestral mode of life of snakes, we find evidence of plasticity in the basic ecological and trophic modes of snakes. Consequently, the macrostomatan condition should not be treated a priori as a derived character state devoid of homoplasy.     

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

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

Molecular phylogeny of gibbons inferred from mitochondrial DNA sequences: Preliminary report

We analyzed the 896 base-pair (bp) mitochondrial DNA (mtDNA) sequences for seven gibbons, representative of three out of four subgenera. The result from our molecular analysis is consistent with previous studies as to the monophyly of subgenus Hylobates species, yet the relationship among subgenera remains slightly ambiguous. A striking result of the analysis is the phylogenetic location of Kloss's gibbon (H. klossii). Kloss's gibbon has been considered to be an initial off-shoot of the subgenus Hylobates because of its morphological primitiveness. However, our molecular data strongly suggest that Kloss's gibbon speciated most recently within the subgenus Hylobates. Correspondence to: S. Horai     

Higher-level relationships of caenophidian snakes inferred from four nuclear and mitochondrial genes

Higher-level caenophidian snake relationships are inferred from sequence analyses of one nuclear gene (C-mos) and three mitochondrial genes (12S rRNA, 16S rRNA and ND4). Caenophidians, which are haenophidian closest relatives, have an Asiatic origin. An African clade comprising atractaspidids, psammophiines, 'lamprophiines' and 'pseudoxyrhophiines' is identified. We discern no evolutionary trend such as an improvement of the venom apparatus with a linear progression from the absence of a venom system to the presence of a front-fanged one. The venom apparatus is contemporary with the origin of colubroids and its absence in a few lineages results from secondary losses. The front-fanged venom system appeared three times independently. The active diurnal foraging mode (associated with a high metabolic rate) appears in a derived position among colubroids.     

Catarrhine primate divergence dates estimated from complete mitochondrial genomes: concordance with fossil and nuclear DNA evidence

Accurate divergence date estimates improve scenarios of primate evolutionary history and aid in interpretation of the natural history of disease-causing agents. While molecule-based estimates of divergence dates of taxa within the superfamily Hominoidea (apes and humans) are common in the literature, few such estimates are available for the Cercopithecoidea (Old World monkeys), the sister taxon of the hominoids in the primate infraorder Catarrhini. To help fill this gap, we have sequenced the entire mitochondrial DNA (mtDNA) genomes from a representative of three cercopithecoid tribes, Cercopithecini (Chlorocebus aethiops), Colobini (Colobus guereza), and Presbytini (Trachypithecus obscurus), and analyzed these new data together with other catarrhine mtDNA genomes available in public databases. Molecular divergence date estimates are dependent on calibration points gleaned from the paleontological record. We defined criteria for the selection of good calibration points and identified three points meeting these criteria: Homo-Pan, 6.0 Ma; Pongo-hominines, 14.0 Ma; hominoid/cercopithecoid, 23.0 Ma. Because a uniform molecular clock does not fit the catarrhine mtDNA data, we estimated divergence dates using a penalized likelihood and a Bayesian method, both of which take into account the effects of rate differences on lineages, phylogenetic tree structure, and multiple calibration points. The penalized likelihood method applied to the coding regions of the mtDNA genome yielded the following divergence date estimates, with approximate 95% confidence intervals: cercopithecine-colobine, 16.2 (14.4-17.9) Ma; colobin-presbytin, 10.9 (9.6-12.3) Ma; cercopithecin-papionin, 11.6 (10.3-12.9) Ma; and Macaca-Papio, 9.8 (8.6-10.9) Ma. Within the hominoids, the following dates were inferred: hylobatid-hominid, 16.8 (15.0-18.5) Ma; Gorilla-Homo+Pan, 8.1 (7.1-9.0) Ma; Pongo pygmaeus pygmaeus-P. p. abelii, 4.1 (3.5-4.7) Ma; and Pan troglodytes-P. paniscus, 2.4 (2.0-2.7) Ma. These dates were similar to those found using penalized likelihood on other subsets of the data, but slightly younger than several of the Bayesian estimates.     

Taenia asiatica and Taenia saginata: genetic divergence estimated from their mitochondrial genomes

We conducted a differential identification of Taenia asiatica and Taenia saginata, through the mapping of mitochondrial genomes and the sequencing of the cox1 and cob genes. The entire mitochondrial genomes of T. asiatica and T. saginata were amplified by long-extension PCR and cloned; each was approximately 14 kb in size. Restriction maps of T. asiatica and T. saginata mitochondrial genomes were then constructed using 13 restriction enzymes. The resulting restriction patterns enable us to estimate their genetic divergence at 4.8%. The actual sequence divergence was computed 4.5% from the cox1 gene, and 4.1% from the cob gene. These results support the designation of T. asiatica as a separate species from T. saginata.     

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.     

Phylogenetic relationships of Nembrothinae (Mollusca: Doridacea: Polyceridae) inferred from morphology and mitochondrial DNA

Within the Polyceridae, Nembrothinae includes some of the most striking and conspicuous sea slugs known, although several features of their biology and phylogenetic relationships remain unknown. This paper reports a phylogenetic analysis based on partial sequences of two mitochondrial genes (cytochrome c oxidase subunit I and 16S rRNA) and morphology for most species included in Nembrothinae. Our phylogenetic reconstructions using both molecular and combined morphological and molecular data support the taxonomic splitting of Nembrothinae into several taxa. Excluding one species (Tambja tentaculata), the monophyly of Roboastra was supported by all the phylogenetic analyses of the combined molecular data. Nembrotha was monophyletic both in the morphological and molecular analyses, always with high support. However, Tambja was recovered as para- or polyphyletic, depending on the analysis performed. Our study also rejects the monophyly of "phanerobranch" dorids based on molecular data.     

乌龟线粒体全基因组序列和结构分析

龟鳖类同其它类群脊椎动物的系统进化关系一直存在争论。为进一步从分子水平上探讨这一问题,本文参照近源物种的线粒体基因组,设计了16对特异引物,采用PCR产物直接测序法测得了乌龟线粒体基因组全序列。结果表明:乌龟线粒体基因组序列全长16576bp,包括2个rRNA基因、22个tRNA基因、13个蛋白质编码基因和1个非编码控制区。乌龟线粒体基因组结构和基因排列顺序与其它龟鳖类相同,在“WANCY区”包含一个“stemloop”结构,ND3基因174位点存在一个额外插入的腺苷酸(A)。本文通过比较分析结构基因在主要脊椎动物类群中的排列顺序,探讨了龟鳖类与其它主要脊椎动物类群的系统进化关系     

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.     

Complete mitochondrial genome of Saturnia jonasii (Lepidoptera: Saturniidae): Genomic comparisons and phylogenetic inference among Bombycoidea

The complete mitochondrial genome (mitogenome) of Saturnia jonasii (Lepidoptera: Saturniidae) was sequenced and compared to those of 19 other bombycoid species. Furthermore, the mitogenome sequences were used to infer phylogenetic relationships among bombycoid species. The 15,261-bp Saturnia jonasii mitogenome contained the typical sets of genes and gene arrangements found in majority of Lepidoptera. All Bombycoidea species, including Saturnia jonasii, have a 15–33-bp spacer sequence at the trnS₂-ND1 junction. The phylogenetic reconstruction of bombycoid species consistently and strongly supported monophylies of the families, Saturniidae, Bombycidae, and Sphingidae, based on Bayesian inference (BI) and maximum-likelihood (ML) methods. Among these families, the Bombycidae and Sphingidae species consistently showed a sister relationship, regardless of data partitions; the BI method strongly supported this relationship, whereas it was moderately supported using the ML method.     

The evolutionary history of cockatoos (Aves: Psittaciformes: Cacatuidae)

Cockatoos are the distinctive family Cacatuidae, a major lineage of the order of parrots (Psittaciformes) and distributed throughout the Australasian region of the world. However, the evolutionary history of cockatoos is not well understood. We investigated the phylogeny of cockatoos based on three mitochondrial and three nuclear DNA genes obtained from 16 of 21 species of Cacatuidae. In addition, five novel mitochondrial genomes were used to estimate time of divergence and our estimates indicate Cacatuidae diverged from Psittacidae approximately 40.7 million years ago (95% CI 51.6-30.3 Ma) during the Eocene. Our data shows Cacatuidae began to diversify approximately 27.9 Ma (95% CI 38.1-18.3 Ma) during the Oligocene. The early to middle Miocene (20-10 Ma) was a significant period in the evolution of modern Australian environments and vegetation, in which a transformation from mainly mesic to xeric habitats (e.g., fire-adapted sclerophyll vegetation and grasslands) occurred. We hypothesize that this environmental transformation was a driving force behind the diversification of cockatoos. A detailed multi-locus molecular phylogeny enabled us to resolve the phylogenetic placements of the Palm Cockatoo (Probosciger aterrimus), Galah (Eolophus roseicapillus), Gang-gang Cockatoo (Callocephalon fimbriatum) and Cockatiel (Nymphicus hollandicus), which have historically been difficult to place within Cacatuidae. When the molecular evidence is analysed in concert with morphology, it is clear that many of the cockatoo species' diagnostic phenotypic traits such as plumage colour, body size, wing shape and bill morphology have evolved in parallel or convergently across lineages.     

Molecular phylogeny of the Chinese ranids inferred from nuclear and mitochondrial DNA sequences

The phylogeny of representative species of Chinese ranids was reconstructed using two nuclear (tyrosinase and rhodopsin) and two mitochondrial (12S rRNA, 16S rRNA) DNA fragments. Maximum parsimony, Bayesian, and maximum likelihood analyses were employed. In comparison with the results from nuclear and mitochondrial data, we used nuclear gene data as our preferred phylogenetic hypothesis. We proposed two families (Ranidae, Dicroglossidae) for Chinese ranids, with the exception of genus Ingerana. Within Dicroglossidae, four tribes were supported including Dicroglossini, Paini, Limnonectini, and Occidozygini. A broader sampling strategy and evidence from additional molecular markers are required to decisively evaluate the evolutionary history of Chinese ranids.     

Phylogenetic estimation of Mytilidae in the East China Sea inferred from mitochondrial genes and nuclear DNA sequence variation

We performed a phylogenetic estimation of the family Mytilidae in the East China Sea based on nuclear internal transcribed spacer (ITS) genes and two mitochondrial genes (COI and 16S RNA). Analysis of five mytilid species based on each of the three genes resulted in mostly congruent trees, although there were some discrepancies in the classification of these species. We combine the results obtained from the three separate analyses to provide a phylogenetic estimation of Mytilidae. We found that the Mytilidae was divided into two major lineages: in one clade, Mytilus galloprovincialis was grouped with Mytilus coruscus; in the second clade, Septifer bilocularis was placed at the basal position in an individual clade, and Perna viridis and Musculista senhousia were recovered as a monophyletic group. Although these finding provide important insights into the taxonomic relationships among the Mytilidae, many aspects of Mytilidae phylogeny remain unresolved. Further analysis based on more molecular information and extensive taxon sampling is necessary to elucidate the phylogenetic relationships among the major lineages within the Mytilidae.     

Evolutionary and taxonomic relationships among Far-Eastern salmonid fishes inferred from mitochondrial DNA divergence

Mitochondrial DNA (mtDNA) restriction analysis was used to examine the evolutionary and taxonomic relationships among 11 taxa of the subfamily Salmoninae. The genera Brachymystax and Hucho were closely related, diverging by sequence divergence estimates of 3.1%. Because the mtDNA sequence divergence between blunt- and sharp-snouted forms of Brachymystax (2.24%) was similar to divergence level of Brachymystax and Hucho , then taking into account the distinct morphological, ecological and allozyme differences between them, it is possible to recognize these forms as two separate species. The subgenus Parahucho formed a very distinct group differing by 6.35–7.08% (sequence divergence estimate) from both Brachymystax and Hucho and must be considered as a valid genus. The UPGMA and neighbour-joined phenograms showed that the five genera studied are divided into two main groupings: (1) Hucho, Brachymystax and Salvelinus ; and (2) Oncorhynchus and Parahucho species. The mtDNA sequence divergence estimates between these groupings were about 8.1%. However, the subsequent bootstrap analysis of mtDNA RFLP data did not support the monophyly of the latter grouping. The concordance of morphological and mtDNA phylogenetic patterns is discussed.     

One-step PCR amplification of complete arthropod mitochondrial genomes

A new PCR primer set which enables one-step amplification of complete arthropod mitochondrial genomes was designed from two conserved 16S rDNA regions for the long PCR technique. For this purpose, partial 16S rDNAs amplified with universal primers 16SA and 16SB were newly sequenced from six representative arthropods: Armadillidium vulgare and Macrobrachium nipponense (Crustacea), Anopheles sinensis (Insecta), Lithobius forficatus and Megaphyllum sp. (Myriapoda), and Limulus polyphemus (Chelicerata). The genomic locations of two new primers, HPK16Saa and HPK16Sbb, correspond to positions 13314-13345 and 12951-12984, respectively, in the Drosophila yakuba mitochondrial genome. The usefulness of the primer set was experimentally examined and confirmed with five of the representative arthropods, except for A. vulgare, which has a linearized mitochondrial genome. With this set, therefore, we could easily and rapidly amplify complete mitochondrial genomes with small amounts of arthropod DNA. Although the primers suggested here were examined only with arthropod groups, a possibility of successful application to other invertebrates is very high, since the high degree of sequence conservation is shown on the primer sites in other invertebrates. Thus, this primer set can serve various research fields, such as molecular evolution, population genetics, and molecular phylogenetics based on DNA sequences, RFLP, and gene rearrangement of mitochondrial genomes in arthropods and other invertebrates.     

Evolutionary landscape of amphibians emerging from ancient freshwater fish inferred from complete mitochondrial genomes

It is very interesting that the only extant marine amphibian is the marine frog, Fejervarya cancrivora. This study investigated the reasons for this apparent rarity by conducting a phylogenetic tree analysis of the complete mitochondrial genomes from 14 amphibians, 67 freshwater fishes, four migratory fishes, 35 saltwater fishes, and one hemichordate. The results showed that amphibians, living fossil fishes, and the common ancestors of modern fishes are phylogenetically separated. In general, amphibians, living fossil fishes, saltwater fishes, and freshwater fishes are clustered in different clades. This suggests that the ancestor of living amphibians arose from a type of primordial freshwater fish, rather than the coelacanth, lungfish, or modern saltwater fish. Modern freshwater fish and modern saltwater fish were probably separated from a common ancestor by a single event, caused by crustal movement.