球孢白僵菌线粒体序列分析及系统进化树构建

线粒体是真核细胞的重要的细胞器,在细胞能量代谢和细胞衰老方面扮演了重要角色,同时也是研究物种进化关系的重要资源。本文应用球孢白僵菌的线粒体序列分析了白僵菌与粪壳菌纲其他11属的13种真菌间的进化关系,为阐明球孢白僵菌的进化地位提供了新的论据。线粒体编码的rnl,rns,25个tRNAs和14个蛋白质基因均由同一条链编码,大部分的tRNA分布为三个tRNA簇,这与粪壳菌纲的其他物种相似。与NCBI上公布的序列比对发现球孢白僵菌线粒体序列与蝇蚧霉线粒体序列最相似,相似率达73%,用14个蛋白编码基因通过邻接法(Neighbor-joint)和最大简约法(Maximum parsimony)构建的系统进化树也得到相同结论,其支持率均为100%。     

The complete mitochondrial genome of the entomopathogenic fungus Metarhizium anisopliae var. anisopliae: gene order and trn gene clusters reveal a common evolutionary course for all Sordariomycetes, while intergenic regions show variation

The mitochondrial genome (mtDNA) of the entomopathogenic fungus Metarhizium anisopliae var. anisopliae, with a total size of 24,673 bp, was one of the smallest known mtDNAs of Pezizomycotina. It contained the 14 typical genes coding for proteins related to oxidative phosphorylation, the two rRNA genes, a single intron that harbored an intronic ORF coding for a putative ribosomal protein (rps) within the large rRNA gene (rnl), and a set of 24 tRNA genes which recognized codons for all amino acids, except proline and valine. Gene order comparison with all known mtDNAs of Sordariomycetes illustrated a highly conserved genome organization for all the protein- and rRNA-coding genes, as well as three clusters of tRNA genes. By considering all mitochondrial essential protein-coding genes as one unit a phylogenetic study of these small genomes strongly supported the common evolutionary course of Sordariomycetes (100% bootstrap support) and highlighted the advantages of analyzing small genomes (mtDNA) over single genes. In addition, comparative analysis of three intergenic regions demonstrated sequence variability that can be exploited for intra- and inter-specific identification of Metarhizium. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

The analysis of the complete mitochondrial genome of Lecanicillium muscarium (synonym Verticillium lecanii) suggests a minimum common gene organization in mtDNAs of Sordariomycetes: phylogenetic implications

The mitochondrial genome (mtDNA) of the entomopathogenic fungus Lecanicillium muscarium (synonym Verticillium lecanii) with a total size of 24,499-bp has been analyzed. So far, it is the smallest known mitochondrial genome among Pezizomycotina, with an extremely compact gene organization and only one group-I intron in its large ribosomal RNA (rnl) gene. It contains the 14 typical genes coding for proteins related to oxidative phosphorylation, the two rRNA genes, one intronic ORF coding for a possible ribosomal protein (rps), and a set of 25 tRNA genes which recognize codons for all amino acids, except alanine and cysteine. All genes are transcribed from the same DNA strand. Gene order comparison with all available complete fungal mtDNAs-representatives of all four Phyla are included-revealed some characteristic common features like uninterrupted gene pairs, overlapping genes, and extremely variable intergenic regions, that can all be exploited for the study of fungal mitochondrial genomes. Moreover, a minimum common mtDNA gene order could be detected, in two units, for all known Sordariomycetes namely nad1-nad4-atp8-atp6 and rns-cox3-rnl, which can be extended in Hypocreales, to nad4L-nad5-cob-cox1-nad1-nad4-atp8-atp6 and rns-cox3-rnl nad2-nad3, respectively. Phylogenetic analysis of all fungal mtDNA essential protein-coding genes as one unit, clearly demonstrated the superiority of small genome (mtDNA) over single gene comparisons.     

Complete mitochondrial genome of the longtooth grouper Epinephelus bruneus (Perciformes, Serranidae)

We determined the complete nucleotide sequence of the mitochondrial (mt) genome for the longtooth grouper, Epinephelus bruneus (Perciformes, Serranidae). This mt genome, consisting of 16,686 base pairs (bp), encoded genes for 13 protein-coding genes, 2 ribosomal RNAs, 22 transfer RNAs, and a noncoding control region as those found in other vertebrates, with the gene order identical to that of typical vertebrates. A major noncoding region between the trnP and trnF genes (991 bp) was considered to be the control region (D-loop). Within this sequence, 22 copies of a 17-bp tandem repeat element, 5'-TGATATTACATATATGC-3', were identified in the control region unlike previous reported Epinephelus species.     

Concatenated mitochondrial DNA of the coccidian parasite Eimeria tenella

Apicomplexan parasites of the genus Plasmodium, pathogens causing malaria, and the genera Babesia and Theileria, aetiological agents of piroplasmosis, are closely related. However, their mitochondrial (mt) genome structures are highly divergent: Plasmodium has a concatemer of 6-kb unit and Babesia/Theileria a monomer of 6.6- to 8.2-kb with terminal inverted repeats. Fragmentation of ribosomal RNA (rRNA) genes and gene arrangements are remarkably distinctive. To elucidate the evolutionary origin of this structural divergence, we determined the mt genome of Eimeria tenella, pathogens of coccidiosis in domestic fowls. Analysis revealed that E. tenella mt genome was concatemeric with similar protein-coding genes and rRNA gene fragments to Plasmodium. Copy number was 50-fold of the nuclear genome. Evolution of structural divergence in the apicomplexan mt genomes is discussed.     

基于太平洋甲胁虱裂化线粒体基因组推测甲胁虱属祖先线粒体核型

[目的]动物典型的单一染色体线粒体基因组在甲胁虱属Hoplopleura已裂化成多个线粒体微环染色体.本研究旨在通过测定太平洋甲胁虱Hoplopleura pacifica的线粒体基因组来推测甲胁虱属祖先线粒体核型.[方法]利用Illumina HiSeq X Ten高通量测序技术对太平洋甲胁虱裂化线粒体基因组进行测定...     

Oesophagostomum dentatum and Oesophagostomum quadrispinulatum: characterization of the complete mitochondrial genome sequences of the two pig nodule worms

In the present study, the complete mitochondrial DNA (mtDNA) sequences of the pig nodule worm Oesophagostomum quadrispinulatum were determined for the first time, and the mt genome of Oesophagostomum dentatum from China was also sequenced for comparative analysis of their gene contents and genome organizations. The mtDNA sequences of O. dentatum China isolate and O. quadrispinulatum were 13,752 and 13,681 bp in size, respectively. Each of the two mt genomes comprises 36 genes, including 12 protein-coding genes, two ribosomal RNA and 22 transfer RNA genes, but lacks the ATP synthetase subunit 8 gene. All genes are transcribed in the same direction and have a nucleotide composition high in A and T. The contents of A+T are 75.79% and 77.52% for the mt genomes of O. dentatum and O. quadrispinulatum, respectively. Phylogenetic analyses using concatenated amino acid sequences of the 12 protein-coding genes, with three different computational algorithms (maximum likelihood, maximum parsimony and Bayesian inference), all revealed that O. dentatum and O. quadrispinulatum represent distinct but closely-related species. These data provide novel and useful markers for studying the systematics, population genetics and molecular diagnosis of the two pig nodule worms.     

Numerous gene rearrangements in the mitochondrial genome of the wallaby louse, Heterodoxus macropus (Phthiraptera)

The complete arrangement of genes in the mitochondrial (mt) genome is known for 12 species of insects, and part of the gene arrangement in the mt genome is known for over 300 other species of insects. The arrangement of genes in the mt genome is very conserved in insects studied, since all of the protein-coding and rRNA genes and most of the tRNA genes are arranged in the same way. We sequenced the entire mt genome of the wallaby louse, Heterodoxus macropus, which is 14,670 bp long and has the 37 genes typical of animals and some noncoding regions. The largest noncoding region is 73 bp long (93% A+T), and the second largest is 47 bp long (92% A+T). Both of these noncoding regions seem to be able to form stem-loop structures. The arrangement of genes in the mt genome of this louse is unlike that of any other animal studied. All tRNA genes have moved and/or inverted relative to the ancestral gene arrangement of insects, which is present in the fruit fly Drosophila yakuba. At least nine protein-coding genes (atp6, atp8, cox2, cob, nad1-nad3, nad5, and nad6) have moved; moreover, four of these genes (atp6, atp8, nad1, and nad3) have inverted. The large number of gene rearrangements in the mt genome of H. macropus is unprecedented for an arthropod.     

Phylogenetic relationships of discoglossid frogs (Amphibia:Anura:Discoglossidae) based on complete mitochondrial genomes and nuclear genes

The complete nucleotide sequence of the mitochondrial (mt) genome was determined for three species of discoglossid frogs (Amphibia:Anura:Discoglossidae), representing three of the four recognized genera: Alytes obstetricans, Bombina orientalis, and Discoglossus galganoi. The organization and size of these newly determined mt genomes are similar to those previously reported for other vertebrates. Phylogenetic analyses (maximum likelihood, Bayesian inference, minimum evolution, and maximum parsimony) of mt protein-coding genes at the amino acid level were performed in combination with already published mt genome sequence data of three species of Neobatrachia, one of Pipoidea, and four of Caudata. Phylogenetic analyses based on the deduced amino acid sequences of all mt protein-coding genes arrived at the same topology. The monophyly of Discoglossidae is strongly supported. Within the Discoglossidae, Alytes is consistently recovered as sister group of Discoglossus, to the exclusion of Bombina. The three species representing Neobatrachia exhibited extremely long branches irrespective of the phylogenetic inference method used, and hence their relative position with respect to Discoglossidae and Xenopus may be artefactual due to a severe long branch attraction effect. To further investigate the phylogenetic intrarelationships of discoglossids, nucleotide sequences of four nuclear protein-coding genes (CXCR4, RAG1, RAG2, and Rhodopsin) with sequences available for the three discoglossid genera and Xenopus were retrieved from GenBank, and together with a concatenated nucleotide sequence data set containing all mt protein-coding genes except ND6 were subjected to separate and combined phylogenetic analyses. In all cases, a sister group relationship between Alytes and Discoglossus was recovered with high statistical support.     

Sequencing and analysis of the complete mitochondrial genome of the giant dobsonfly Acanthacorydalis orientalis(McLachlan)(Insecta:Megaloptera:Corydalidae)

The complete mitochondrial genome of Acanthacorydalis orientalis(McLachlan)was determined and analyzed(GenBank accession number:KF840564).This paper represents the first mitochondrial(mt) genome of the dobsonfly genus Acanthacorydalis.The mt genome is a typical circular DNA of 15 753 bp composed of 37 genes with an A+T content of 76.7%.It has an ancestral gene arrangement of the insect mt genomes.Eleven of the 13 PCGs start with codon ATT and ATG,while several exceptions such as ATA and TTG respectively for atpS and nad\ are also present.Five protein-coding genes end with a single T,while others have a termination codon of TAA or TAG.Most tRNAs are folded into the typical clover-leaf structure except for the trnS 1 whose dihydrouridine arm was a simple loop.The secondary structure of rrnl consists of five structural domains and 50 helices,while the rrns includes three domains and 34 helices.The control region has a stretches of Ts with a length of 22 bp but lacks obvious tandem repeat region.Both Bayesian inference and maximum likelihood(ML) analyses,based on all 13 protein-coding genes and two rRNA genes of the mt genomes,confirm the monophyly of Corydalinae and suggest that Acanthacorydalis,together with Corydalus,which is an endemic dobsonfly genus from the New World,belong to a monophyletic lineage.     

The mitochondrial genome of Priapulus caudatus Lamarck (Priapulida: Priapulidae)

We sequenced and annotated the complete mitochondrial (mt) genome of the priapulid Priapulus caudatus in order to provide a source of phylogenetic characters including an assessment of gene order arrangement. The genome was 14,919 bp in its entirety with few, short non-coding regions. A number of protein-coding and tRNA genes overlapped, making the genome relatively compact. The gene order was: cox1, cox2, trnK, trnD, atp8, atp6, cox3, trnG, nad3, trnA, trnR, trnN, rrnS, trnV, rrnL, trnL(yaa), trnL(nag), nad1, -trnS(nga), -cob, -nad6, trnP, -trnT, nad4L, nad4, trnH, nad5, trnF, -trnE, -trnS(nct), trnI, -trnQ, trnM, nad2, trnW, -trnC, -trnY; where '-' indicates genes transcribed on the opposite strand. The gene order, although unique amongst Metazoa, shared the greatest number of gene boundaries and the longest contiguous fragments with the chelicerate Limulus polyphemus. The mt genomes of these taxa differed only by a single inversion of 18 contiguous genes bounded by rrnS and trnS(nct). Other arthropods and nematodes shared fewer gene boundaries but considerably more than the most similar non-ecdysozoan.     

The complete mitochondrial genome of Arctic Calanus hyperboreus (Copepoda,Calanoida) reveals characteristic patterns in calanoid mitochondrial genome

Copepoda is the most diverse and abundant group of crustaceans, but its phylogenetic relationships are ambiguous. Mitochondrial (mt) genomes are useful for studying evolutionary history, but only six complete Copepoda mt genomes have been made available and these have extremely rearranged genome structures. This study determined the mt genome of Calanus hyperboreus, making it the first reported Arctic copepod mt genome and the first complete mt genome of a calanoid copepod. The mt genome of C. hyperboreus is 17,910 bp in length and it contains the entire set of 37 mt genes, including 13 protein-coding genes, 2 rRNAs, and 22 tRNAs. It has a very unusual gene structure, including the longest control region reported for a crustacean, a large tRNA gene cluster, and reversed GC skews in 11 out of 13 protein-coding genes (84.6%). Despite the unusual features, comparing this genome to published copepod genomes revealed retained pan-crustacean features, as well as a conserved calanoid-specific pattern. Our data provide a foundation for exploring the calanoid pattern and the mechanisms of mt gene rearrangement in the evolutionary history of the copepod mt genome.     

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.     

The complete mitochondrial genome of the Asiatic cavity-nesting honeybee Apis cerana (Hymenoptera: Apidae)

In the present study, we determined the complete mitochondrial DNA (mtDNA) sequence of Apis cerana, the Asiatic cavity-nesting honeybee. We present here an analysis of features of its gene content and genome organization in comparison with Apis mellifera to assess the variation within the genus Apis and among main groups of Hymenoptera. The size of the entire mt genome of A. cerana is 15,895 bp, containing 2 ribosomal RNA genes, 13 protein-coding genes, 22 transfer RNA (tRNA) genes and one control region. These genes are transcribed from both strands and have a nucleotide composition high in A and T. The contents of A+T of the complete genomes are 83.96% for A. cerana. The AT bias had a significant effect on both the codon usage pattern and amino acid composition of proteins. There are a total of 3672 codons in all 13 protein-coding genes, excluding termination codons. The most frequently used amino acid is Leu (15.52%), followed by Ile (12.85%), Phe (10.10%), Ser (9.15%) and Met (8.96%). Intergenic regions in the mt genome of A. cerana are 705 bp in total. The order and orientation of the gene arrangement pattern is identical to that of A. mellifera, except for the position of the tRNA-Ser(AGN) gene. Phylogenetic analyses using concatenated amino acid sequences of 13 protein-coding genes, with three different computational algorithms (NJ, MP and ML), all revealed two distinct groups with high statistical support, indicating that A. cerana and A. mellifera are two separate species, consistent with results of previous morphological and molecular studies. The complete mtDNA sequence of A. cerana provides additional genetic markers for studying population genetics, systematics and phylogeographics of honeybees.     

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 mitochondrial genome of Strongyloides stercoralis (Nematoda) - idiosyncratic gene order and evolutionary implications

The complete mitochondrial genome sequence of the parasitic nematode Strongyloides stercoralis was determined, and its organisation and structure compared with other nematodes for which complete mitochondrial sequence data were available. The mitochondrial genome of S. stercoralis is 13,758 bp in size and contains 36 genes (all transcribed in the clockwise direction) but lacks the atp8 gene. This genome has a high T content (55.9%) and a low C content (8.3%). Corresponding to this T content, there are 16 (poly-T) tracts of >/=12 Ts distributed across the genome. In protein-coding genes, the T bias is greatest (76.4%) at the third codon position compared with the first and second codon positions. Also, the C content is higher at the first (9.3%) and second (13.4%) codon positions than at the third (2%) position. These nucleotide biases have a significant effect on predicted codon usage patterns and, hence, on amino acid compositions of the mitochondrial proteins. Interestingly, six of the 12 protein-coding genes are predicted to employ a unique initiation codon (TTT), which has not yet been reported for any other animal mitochondrial genome. The secondary structures predicted for the 22 transfer RNA (trn) genes and the two ribosomal RNA (rrn) genes are similar to those of other nematodes. In contrast, the gene arrangement in the mitochondrial genome of S. stercoralis is different from all other nematodes studied to date, revealing only a limited number of shared gene boundaries (atp6-nad2 and cox2-rrnL). Evolutionary analyses of mitochondrial nucleotide and amino acid sequence data sets for S. stercoralis and seven other nematodes demonstrate that the mitochondrial genome provides a rich source of phylogenetically informative characters. In conclusion, the S. stercoralis mitochondrial genome, with its unique gene order and characteristics, should provide a resource for comparative mitochondrial genomics and systematics studies of parasitic nematodes.     

Complete sequence of the mitochondrial genome of Taenia saginata: comparison with T. solium and T. asiatica

The complete sequence of the Taenia saginata mitochondrial genome was determined, and its organization and structure were compared to other human-tropic Taenia tapeworms for which complete mitochondrial sequence data were available. The mitochondrial genome was 13,670 bp long, contained 12 protein-coding genes, two ribosomal RNAs (rRNAs, a small and a large subunit), and 22 transfer RNAs (tRNAs). It did not encode the atp8 gene. Overlapping regions were found between nad4L and nad4, nad1 and trnN, and cox1 and trnT. The ATG initiation codon was used for 10 protein-coding genes, and the GTG initiation codon was used for the remaining 2 genes (nad4 and atp6). The size of the protein-coding genes of the three human Taenia tapeworms did not vary, except for Taenia solium nad1 (891 aa) and nad4 (1212 aa) and Taenia asiatica cox2 (576 aa). The tRNA genes were 57-75 bp long, and the predicted secondary structures of 18 of these genes had typical clover-leaf shapes with paired dihydrouridine (DHU) arms. The genes in all human Taenia tapeworms for the two mitochondrial rRNA subunits rrnL and rrnS are separated by trnC. The putative T. saginata rrnL and rrnS are 972 and 732 bp long, respectively. The non-coding regions of the mt genome of T. saginata consisted of 2 regions: a short non-coding region (SNR, 66 nucleotides) and a long non-coding region (LNR, 159 nucleotides). The overall sequence difference in the full mitochondrial genome between T. saginata and T. asiatica was 4.6%, while T. solium differed by 11%. In conclusion, the complete sequence of the T. saginata mitochondrial genome will serve as a resource for comparative mitochondrial genomics and systematic studies of the parasitic cestodes.     

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hemipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one another. We found that the correlation was positive and statistically significant (R2 = 0.73, P = 0.01; Rs = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.