Comparative Analysis of Mitochondrial Genomes in <Emphasis Type="Italic">Bombina</Emphasis> (Anura; Bombinatoridae)

The complete mitochondrial genomes of two basal anurans, Bombina bombina and B. variegata (Anura; Bombinatoridae), were sequenced. The gene order of their mitochondrial DNA (mtDNA) is identical to that of canonical vertebrate mtDNA. In contrast, we show that there are structural differences in regulatory regions and protein coding genes between the mtDNA of these two closely related species. Corrected sequence divergence between the mtDNA of B. bombina and B. variegata amounts to 8.7% (2.3% divergence in amino acids). Comparisons with two East Asian congeners show that the control region contains two repeat regions, LV1 and LV2, present in all species except for B. bombina, in which LV2 has been secondarily lost. The rRNAs and tRNAs are characterized by low nucleotide divergence. The protein coding genes are considerably more disparate, although functional constraint is high but variable among genes, as evidenced by dN/dS ratios. A mtDNA phylogeny established the distribution of autapomorphic nonsynonomous substitutions in the mitogenomes of B. bombina and B. variegata. Nine of 98 nonsynonomous substitutions led to radical amino acid replacements that may alter mitochondrial protein function. Most radical substitutions were found in ND2, ND4, or ND5, encoding mitochondrial subunits of complex I of the electron transport system. The extensive divergence between the mitogenomes of B. bombina and B. variegata is discussed in terms of its possible role in impeding gene flow in natural hybrid zones between these two species.     

A low rate of replacement substitutions in two major Ovis aries mitochondrial genomes

Two mitochondrial DNA molecules which represent major Ovis aries mtDNA haplogroups were cloned and comparatively sequenced to assess the degree of intraspecific variation. A total of 9623 bp that correspond to 58% of both mitochondrial genomes were determined. The control region, the Cyt b , ND2, ND3, ND4L, COIII and 12 tRNA genes, including the origin of L-strand replication, were completely characterized. Partial sequence information was obtained from the 12S and 16S rRNA and an additional six protein coding and six tRNA genes. The control regions of the two mtDNAs showed a nucleotide divergence of 4·34% while coding regions differed by 0·44%. The number of sheep coding region substitutions was similar to values observed in intraspecific comparisons of mitochondrial DNAs that represent remote points in genealogical trees of mice and humans. However, replacement substitutions were only observed at ∼30% of the rate in mice and ∼20% of the rate in humans. Nucleotide substitutions with a potential for phenotypic effects were found in the 12S and 16S rRNA and in the ND1 and COIII genes.     

中华卵索线虫线粒体基因组多态性分析

为完善昆虫病原索科线虫线粒体基因组全序列数据库, 更系统地研究其基因组特征和系统演化规律, 进而为发挥该线虫生防潜力打下基础, 我们开展了中华卵索线虫Ovomermis sinensis线粒体全基因组的研究。该研究通过线粒体基因组滚环复制及酶切图谱, 揭示了中华卵索线虫线粒体基因组具有种内遗传多态性, 即群体中单体线虫具有独特的酶切条带, 且条带累加之和变化范围较大, 为16.5～24.5 kb。为进一步了解线粒体基因组多态性特征及产生的分子机制, 采用两步长PCR方法对2条代表性成虫线粒体基因组进行了测序及拼接, 得其全长分别为18 864和16 777 bp。对这2条序列的比对表明, 线粒体基因组中位于ND2和ND4之间的可变区域, 不仅基因排列顺序不同, 且存在ND3基因重复现象, 这是导致中华卵索线虫线粒体基因组呈现多态性的主要原因。通过对以上研究结果的分析及与GenBank中已有的6种索科线虫线粒体基因组序列进行比对, 概括出其线粒体基因组基本特点: ①线粒体基因排列顺序各不相同;②部分线虫线粒体基因存在重复现象, 且重复次数不同;③线粒体基因组大小存在很大差异。     

The mitochondrial genome sequence of Enterobius vermicularis (Nematoda: Oxyurida)--an idiosyncratic gene order and phylogenetic information for chromadorean nematodes

The complete mitochondrial genome sequence was determined for the human pinworm Enterobius vermicularis (Oxyurida: Nematoda) and used to infer its phylogenetic relationship to other major groups of chromadorean nematodes. The E. vermicularis genome is a 14,010-bp circular DNA molecule that encodes 36 genes (12 proteins, 22 tRNAs, and 2 rRNAs). This mtDNA genome lacks atp8, as reported for almost all other nematode species investigated. Phylogenetic analyses (maximum parsimony, maximum likelihood, neighbor joining, and Bayesian inference) of nucleotide sequences for the 12 protein-coding genes of 25 nematode species placed E. vermicularis, a representative of the order Oxyurida, as sister to the main Ascaridida+Rhabditida group. Tree topology comparisons using statistical tests rejected an alternative hypothesis favoring a closer relationship among Ascaridida, Spirurida, and Oxyurida, which has been supported from most studies based on nuclear ribosomal DNA sequences. Unlike the relatively conserved gene arrangement found for most chromadorean taxa, E. vermicularis mtDNA gene order is very unique, not sharing similarity to any other nematode species reported to date. This lack of gene order similarity may represent idiosyncratic gene rearrangements unique to this specific lineage of the oxyurids. To more fully understand the extent of gene rearrangement and its evolutionary significance within the nematode phylogenetic framework, additional mitochondrial genomes representing a greater evolutionary diversity of species must be characterized.     

Mitochondrial DNA sequences of primates: Tempo and mode of evolution

Summary We cloned and sequenced a segment of mitochondrial DNA from human, chimpanzee, gorilla, orangutan, and gibbon. This segment is 896 bp in length, contains the genes for three transfer RNAs and parts of two proteins, and is homologous in all 5 primates. The 5 sequences differ from one another by base substitutions at 283 positions and by a deletion of one base pair. The sequence differences range from 9 to 19% among species, in agreement with estimates from cleavage map comparisons, thus confirming that the rate of mtDNA evolution in primates is 5 to 10 times higher than in nuclear DNA. The most striking new finding to emerge from these comparisons is that transitions greatly outnumber transversions. Ninety-two percent of the differences among the most closely related species (human, chimpanzee, and gorilla) are transitions. For pairs of species with longer divergence times, the observed percentage of transitions falls until, in the case of comparisons between primates and non-primates, it reaches a value of 45. The time dependence is probably due to obliteration of the record of transitions by multiple substitutions at the same nucleotide site. This finding illustrates the importance of choosing closely related species for analysis of the evolutionary process. The remarkable bias toward transitions in mtDNA evolution necessitates the revision of equations that correct for multiple substitutions at the same site. With revised equations, we calculated the incidence of silent and replacement substitutions in the two protein-coding genes. The silent substitution rate is 4 to 6 times higher than the replacement rate, indicating strong functional constraints at replacement sites. Moreover, the silent rate for these two genes is about 10% per million years, a value 10 times higher than the silent rate for the nuclear genes studied so far. In addition, the mean substitution rate in the three mitochondrial tRNA genes is at least 100 times higher than in nuclear tRNA genes. Finally, genealogical analysis of the sequence differences supports the view that the human lineage branched off only slightly before the gorilla and chimpanzee lineages diverged and strengthens the hypothesis that humans are more related to gorillas and chimpanzees than is the orangutan.Abbreviations mtDNA mitochondrial DNA - bp base pair - URF unidentified reading frame     

Mode and Tempo of Molecular Evolution in the Nematode Caenorhabditis: Cytochrome Oxidase II and Calmodulin Sequences

Through direct sequencing methods, the mitochondrial gene for cytochrome oxidase subunit two (CO II) and the single-copy nuclear gene for calmodulin were compared among strains of Caenorhabidits elegans and two other Caenorhabditis species (C. remanei and C. briggsae). In addition the CO II sequence was determined from a distantly related nematode, Steinernema intermedii. Among the 11 strains of C. elegans tested, there are four types of CO II gene, arising from two major lineages. Levels of intraspecific difference in the CO II gene are low (less than 2.0%) compared to the extraordinary divergence between congeneric species, which is about 50% when corrected for multiple hits. Concordant with the increase in divergence between taxa is a change in the pattern of substitution from a strong transition bias (24 transitions compared to two transversions) within species to a substitution pattern that appears to reflect the base composition of the mitochondrial genome when more divergent nematodes are compared. The base composition of the Caenorhabditis CO II gene is strongly biased toward A + T at all three positions of codons and appears to constrain the amino acid composition of the protein. Both the CO II and calmodulin genes show extreme conservation of amino acid sequences. When the accumulation of changes at silent sites in the two genes is compared among strains, it becomes evident that the mitochondrial gene is changing faster than the nuclear gene.     

Elevated rates of nonsynonymous substitution in island birds

Slightly deleterious mutations are expected to fix at relatively higher rates in small populations than in large populations. Support for this prediction of the nearly-neutral theory of molecular evolution comes from many cases in which lineages inferred to differ in long-term average population size have different rates of nonsynonymous substitution. However, in most of these cases, the lineages differ in many other ways as well, leaving open the possibility that some factor other than population size might have caused the difference in substitution rates. We compared synonymous and nonsynonymous substitutions in the mitochondrial cyt b and ND2 genes of nine closely related island and mainland lineages of ducks and doves. We assumed that island taxa had smaller average population sizes than those of their mainland sister taxa for most of the time since they were established. In all nine cases, more nonsynonymous substitutions occurred on the island branch, but synonymous substitutions showed no significant bias. As in previous comparisons of this kind, the lineages with smaller populations might differ in other respects that tend to increase rates of nonsynonymous substitution, but here such differences are expected to be slight owing to the relatively recent origins of the island taxa. An examination of changes to apparently "preferred" and "unpreferred" synonymous codons revealed no consistent difference between island and mainland lineages.     

Synonymous substitution rates in Drosophila: Mitochondrial versus nuclear genes

Synonymous substitution rates in mitochondrial and nuclear genes of Drosophila were compared. To make accurate comparisons, we considered the following: (1) relative synonymous rates, which do not require divergence time estimates, should be used; (2) methods estimating divergence should take into account base composition; (3) only very closely related species should be used to avoid effects of saturation; (4) the heterogeneity of rates should be examined. We modified the methods estimating synonymous substitution numbers to account for base composition bias. By using these methods, we found that mitochondrial genes have 1.7–3.4 times higher synonymous substitution rates than the fastest nuclear genes or 4.5–9.0 times higher rates than the average nuclear genes. The average rate of synonymous transversions was 2.7 (estimated from the melanogaster species subgroup) or 2.9 (estimated from the obscura group) times higher in mitochondrial genes than in nuclear genes. Synonymous transversions in mitochondrial genes occurred at an approximately equivalent rate to those in the fastest nuclear genes. This last result is not consistent with the hypothesis that the difference in turnover rates between mitochondrial and nuclear genomes is the major factor determining higher synonymous substitution rates in mtDNA. We conclude that the difference in synonymous substitution rates is due to a combination of two factors: a higher transitional mutation rate in mtDNA and constraints on nuclear genes due to selection for codon usage. Received: 27 November 1996 / Accepted: 8 May 1997     

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

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

Sequences from 14 mitochondrial genes provide a well-supported phylogeny of the Charadriiform birds congruent with the nuclear RAG-1 tree

Because of the difficulties of constructing a robust phylogeny for Charadriiform birds using morphological characters, recent studies have turned to DNA sequences to resolve the systematic uncertainties of family-level relationships in this group. However, trees constructed using nuclear genes or the mitochondrial Cytochrome b gene suggest deep-level relationships of shorebirds that differ from previous studies based on morphology or DNA-DNA hybridization distances. To test phylogenetic hypotheses based on nuclear genes (RAG-1, myoglobin intron-2) and single mitochondrial genes (Cytochrome b), approximately 13,000 bp of mitochondrial sequence was collected for one exemplar species of 17 families of Charadriiformes plus potential outgroups. Maximum likelihood and Bayesian analyses show that trees constructed from long mitochondrial sequences are congruent with the nuclear gene topologies [Chardrii (Lari, Scolopaci)]. Unlike short mitochondrial sequences (such as Cytochrome b alone), longer sequences yield a well-supported phylogeny for shorebirds across various taxonomic levels. Examination of substitution patterns among mitochondrial genes reveals specific genes (especially ND5, ND4, ND2, and COI) that are better suited for phylogenetic analyses among shorebird families because of their relatively homogeneous nucleotide composition among lineages, slower accumulation of substitutions at third codon positions, and phylogenetic utility in both closely and distantly related lineages. For systematic studies of birds in which family and generic levels are examined simultaneously, we recommend the use of both nuclear and mitochondrial sequences as the best strategy to recover relationships that most likely reflect the phylogenetic history of these lineages.     

Neutrality tests on mtDNA: unusual results from nematodes

McDonald-Kreitman tests of neutrality on mitochondrial DNA (mtDNA) of butterflies, Drosophila, and a variety of vertebrates usually show excess (over the neutral expectation) intraspecific polymorphism at nonsilent sites. These results are of great interest because they are the opposite of what is usually found for nuclear genes, in which the neutral pattern or evidence of adaptive divergence between species is usually observed. However, only vertebrates and insects have been tested so far, so it is not clear whether this intriguing pattern is typical for mtDNA in all taxa. Here I tested three pairs of nematode species and found that they all show a deficit of replacement polymorphism. Taken at face value, this result suggests that adaptive evolution proceeds more efficiently in nematode mtDNA than in mtDNA of vertebrates or insects. An alternate explanation is that the nematode pattern is an artifact of silent-site saturation that results from the rapid and composition-biased way in which nematode mtDNA evolves. Further studies are needed to distinguish between these two hypotheses.     

Molecular systematics of cytochrome oxidase I and 16S from Neochlamisus leaf beetles and the importance of sampling.

If a gene tree is to be judiciously used for inferring the histories of closely related taxa, (1) its topology must be sufficiently resolved and robust that noteworthy phylogenetic patterns can be confidently documented, and (2) sampling of species, populations, and pertinent biological variation must be sufficiently broad that otherwise misleading sources of genetic variation can be detected. These principles are illustrated by the complex gene tree of Neochlamisus leaf beetles that I reconstructed using 90,000 bp of cytochrome oxidase I (COI) and 16S mitochondrial DNA (mtDNA) sequences from over 100 specimens. Cytochrome oxidase I haplotypes varied up to 25.1% within Neochlamisus and up to 11.1% within the gibbosus species group, while exhibiting very low A + T bias for insect mtDNA (63%), low transition saturation, and conservative patterns of amino acid variation. 16S exhibited lower sequence divergences and greater A + T bias and transition saturation than COI, and substitutions were more constrained in stems than in loops. Comparisons with an earlier study of Ophraella leaf beetles highlighted conservative and labile elements of molecular evolution across genes and taxa. Cytochrome oxidase I parsimony and neighbor-joining analyses strongly supported a robust mtDNA genealogy that revealed the monophyly of Neochlamisus and of the gibbosus species group. Phylogeographic relationships suggested that the eastern U.S. gibbosus group derives from southwestern velutinus group ancestors. Haplotypes from individual velutinus group species clustered monophyletically, as expected. However, haplotypes from each of several gibbosus group taxa were polyphyletically distributed, appearing in divergent parts of the tree. 16S provided a less-resolved gibbosus group topology that was congruent with the COI tree and corroborated patterns of mitochondrial polyphyly. By subsampling haplotypes corresponding to particular species, populations, and ecological variants of gibbosus group taxa, I demonstrate that recovered topologies and genetic distances vary egregiously according to sampling regime. This study thus documents the potentially dire consequences of inadequate sampling when inferring the evolutionary history of closely related and mitochondrially polyphyletic taxa.     

Excess amino acid polymorphism in mitochondrial DNA: contrasts among genes from Drosophila, mice, and humans

Recent studies of mitochondrial DNA (mtDNA) variation in mammals and Drosophila have shown an excess of amino acid variation within species (replacement polymorphism) relative to the number of silent and replacement differences fixed between species. To examine further this pattern of nonneutral mtDNA evolution, we present sequence data for the ND3 and ND5 genes from 59 lines of Drosophila melanogaster and 29 lines of D. simulans. Of interest are the frequency spectra of silent and replacement polymorphisms, and potential variation among genes and taxa in the departures from neutral expectations. The Drosophila ND3 and ND5 data show no significant excess of replacement polymorphism using the McDonald-Kreitman test. These data are in contrast to significant departures from neutrality for the ND3 gene in mammals and other genes in Drosophila mtDNA (cytochrome b and ATPase 6). Pooled across genes, however, both Drosophila and human mtDNA show very significant excesses of amino acid polymorphism. Silent polymorphisms at ND5 show a significantly higher variance in frequency than replacement polymorphisms, and the latter show a significant skew toward low frequencies (Tajima's D = -1.954). These patterns are interpreted in light of the nearly neutral theory where mildly deleterious amino acid haplotypes are observed as ephemeral variants within species but do not contribute to divergence. The patterns of polymorphism and divergence at charge-altering amino acid sites are presented for the Drosophila ND5 gene to examine the evolution of functionally distinct mutations. Excess charge-altering polymorphism is observed at the carboxyl terminal and excess charge-altering divergence is detected at the amino terminal. While the mildly deleterious model fits as a net effect in the evolution of nonrecombining mitochondrial genomes, these data suggest that opposing evolutionary pressures may act on different regions of mitochondrial genes and genomes.      

Phylogeny, rate variation, and genome size evolution of Pelargonium (Geraniaceae)

The phylogeny of 58 Pelargonium species was estimated using five plastid markers (rbcL, matK, ndhF, rpoC1, trnL-F) and one mitochondrial gene (nad5). The results confirmed the monophyly of three major clades and four subclades within Pelargonium but also indicate the need to revise some sectional classifications. This phylogeny was used to examine karyotype evolution in the genus: plotting chromosome sizes, numbers and 2C-values indicates that genome size is significantly correlated with chromosome size but not number. Accelerated rates of nucleotide substitution have been previously detected in both plastid and mitochondrial genes in Pelargonium, but sparse taxon sampling did not enable identification of the phylogenetic distribution of these elevated rates. Using the multigene phylogeny as a constraint, we investigated lineage- and locus-specific heterogeneity of substitution rates in Pelargonium for an expanded number of taxa and demonstrated that both plastid and mitochondrial genes have had accelerated substitution rates but with markedly disparate patterns. In the plastid, the exons of rpoC1 have significantly accelerated substitution rates compared to its intron and the acceleration was mainly due to nonsynonymous substitutions. In contrast, the mitochondrial gene, nad5, experienced substantial acceleration of synonymous substitution rates in three internal branches of Pelargonium, but this acceleration ceased in all terminal branches. Several lineages also have dN/dS ratios significantly greater than one for rpoC1, indicating that positive selection is acting on this gene, whereas the accelerated synonymous substitutions in the mitochondrial gene are the result of elevated mutation rates.     

c-mos variation in songbirds: molecular evolution, phylogenetic implications, and comparisons with mitochondrial differentiation

Nucleotide sequences from the c-mos proto-oncogene have previously been used to reconstruct the phylogenetic relationships between distantly related vertebrate taxa. To explore c-mos variation at shallower levels of avian divergence, we compared c-mos sequences from representative passerine taxa that span a range of evolutionary differentiation, from basal passerine lineages to closely allied genera. Phylogenetic reconstructions based on these c-mos sequences recovered topologies congruent with previous DNA-DNA hybridization-based reconstructions, with many nodes receiving high support, as indicated by bootstrap and reliability values. One exception was the relationship of Acanthisitta to the remaining passerines, where the c-mos-based searches indicated a three-way polytomy involving the Acanthisitta lineage and the suboscine and oscine passerine clades. We also compared levels of c-mos and mitochondrial differentiation across eight oscine passerine taxa and found that c-mos nucleotide substitutions accumulate at a rate similar to that of transversion substitutions in mitochondrial protein-coding genes. These comparisons suggest that nuclear-encoded loci such as c-mos provide a temporal window of phylogenetic resolution that overlaps the temporal range where mitochondrial protein-coding sequences have their greatest utility and that c-mos substitutions and mtDNA transversions can serve as complementary, informative, and independent phylogenetic markers for the study of avian relationships.     

Analysis of nucleotide substitutions of mitochondrial DNAs in Drosophila melanogaster and its sibling species

To study the rate and pattern of nucleotide substitution in mitochondrial DNA (mtDNA), we cloned and sequenced a 975-bp segment of mtDNA from Drosophila melanogaster, D. simulans, and D. mauritiana containing the genes for three transfer RNAs and parts of two protein- coding genes, ND2 and COI. Statistical analysis of synonymous substitutions revealed a predominance of transitions over transversions among the three species, a finding differing from previous results obtained from a comparison of D. melanogaster and D. yakuba. The number of transitions observed was nearly the same for each species comparison, including D. yakuba, despite the differences in divergence times. However, transversions seemed to increase steadily with increasing divergence time. By contrast, nonsynonymous substitutions in the ND2 gene showed a predominance of transversions over transitions. Most transversions were between A and T and seemed to be due to some kind of mutational bias to which the A + T-rich mtDNA of Drosophila species may be subject. The overall rate of nucleotide substitution in Drosophila mtDNA appears to be slightly faster (approximately 1.4 times) than that of the Adh gene. This contrasts with the result obtained for mammals, in which the mtDNA evolves approximately 10 times faster than single-copy nuclear DNA. We have also shown that the start codon of the COI gene is GTGA in D. simulans and GTAA in D. mauritiana. These codons are different from that of D. melanogaster (ATAA).      

Human Mitochondrial DNA Variation and Evolution: Analysis of Nucleotide Sequences from Seven Individuals

We have analyzed nucleotide sequence variation in an approximately 900-base pair region of the human mitochondrial DNA molecule encompassing the heavy strand origin of replication and the D-loop. Our analysis has focused on nucleotide sequences available from seven humans. Average nucleotide diversity among the sequences is 1.7%, several-fold higher than estimates from restriction endonuclease site variation in mtDNA from these individuals and previously reported for other humans. This disparity is consistent with the rapidly evolving nature of this noncoding region. However, several instances of convergent or parallel gain and loss of restriction sites due to multiple substitutions were observed. In addition, other results suggest that restriction site (as well as pairwise sequence) comparisons may underestimate the total number of substitutions that have occurred since the divergence of two mtDNA sequences from a common ancestral sequence, even at low levels of divergence. This emphasizes the importance of recognizing the large standard errors associated with estimates of sequence variability, particularly when constructing phylogenies among closely related sequences. Analysis of the observed number and direction of substitutions revealed several significant biases, most notably a strand dependence of substitution type and a 32-fold bias favoring transitions over transversions. The results also revealed a significantly nonrandom distribution of nucleotide substitutions and sequence length variation. Significantly more multiple substitutions were observed than expected for these closely related sequences under the assumption of uniform rates of substitution. The bias for transitions has resulted in predominantly convergent or parallel changes among the observed multiple substitutions. There is no convincing evidence that recombination has contributed to the mtDNA sequence diversity we have observed.     

A molecular phylogeny of the New World orioles (Icterus): the importance of dense taxon sampling.

We sequenced 2005 bp of the mitochondrial ND2 and cytochrome b genes from the 25 recognized species of New World orioles (Icterus). Our data confirmed the monophyly of Icterus and produced a well-resolved phylogeny showing three main clades of orioles. We also sequenced multiple subspecies for most polytypic taxa. Our findings demonstrated the importance of dense taxon sampling below the species level in two ways. First, we found evidence that two species are polyphyletic, I. galbula (Northern oriole) and I. dominicensis (Black-cowled oriole). Choosing different subspecies from either of these taxa would lead to different species-level phylogenies. Second, adding subspecies even to monophyletic groups produced a bootstrap tree with significantly more support. Of the two genes that we used, ND2 provided more resolution than did cytochrome b. ND2 evolved up to 40% faster than cytochrome b, yet shows a higher saturation threshold. Our findings suggest that ND2 may be superior to cytochrome b for resolving species-level phylogenies in passerine birds.     

Nematode mitochondrial metagenomics: A new tool for biodiversity analysis

DNA barcoding approaches have greatly increased our understanding of biodiversity on the planet, and metabarcoding is widely used for classifying members of the phylum Nematoda. However, loci typically utilized in metabarcoding studies are often unable to resolve closely related species or are unable to recover all taxa present in a sample due to inadequate PCR primer binding. Mitochondrial metagenomics (mtMG) is an alternative approach utilizing shotgun sequencing of total DNA to recover the mitochondrial genomes of all species present in samples. However, this approach requires a comprehensive reference database for identification and currently available mitochondrial sequences for nematodes are highly dominated by sequences from the order Rhabditida, and excludes many clades entirely. Here, we analysed the efficacy of mtMG for the recovery of nematode taxa and the generation of mitochondrial genomes. We first developed a curated reference database of nematode mitochondrial sequences and expanded it with 40 newly sequenced taxa. We then tested the mito-metagenomics approach using a series of nematode mock communities consisting of morphologically identified nematode species representing various feeding traits, life stages, and phylogenetic relationships. We were able to identify all but two species through the de novo assembly of COX1 genes. We were also able to recover additional mitochondrial protein coding genes (PCGs) for 23 of the 24 detected species including a full array of 12 PCGs from five of the species. We conclude that mtMG offers a potential for the effective recovery of nematode biodiversity but remains limited by the breadth of the reference database.     

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.