Variation in Coding (NADH Dehydrogenase Subunits 2, 3, and 6) and Noncoding Intergenic Spacer Regions of the Mitochondrial Genome in Octocorallia (Cnidaria: Anthozoa)

Low rates of evolution in cnidarian mitochondrial genes such as COI and 16S rDNA have hindered molecular systematic studies in this important invertebrate group. We sequenced fragments of 3 mitochondrial protein-coding genes (NADH dehydrogenase subunits ND2, ND3 and ND6) as well as the COI-COII intergenic spacer, the longest noncoding region found in the octocoral mitochondrial genome, to determine if any of these regions contain levels of variation sufficient for reconstruction of phylogenetic relationships among genera of the anthozoan subclass Octocorallia. Within and between the soft coral families Alcyoniidae and Xeniidae, sequence divergence in the genes ND2 (539 bp), ND3 (102 bp), and ND6 (444 bp) ranged from 0.5% to 12%, with the greatest pairwise distances between the 2 families. The COI-COII intergenic spacer varied in length from 106 to 122 bp, and pairwise sequence divergence values ranged from 0% to 20.4%. Phylogenetic trees constructed using each region separately were poorly resolved. Better phylogenetic resolution was obtained in a combined analysis using all 3 protein-coding regions (1085 bp total). Although relationships among some pairs of species and genera were well supported in the combined analysis, the base of the alcyoniid family tree remained an unresolved polytomy. We conclude that variation in the NADH subunit coding regions is adequate to resolve phylogenetic relationships among families and some genera of Octocorallia, but insufficient for most species - or population-level studies. Although the COI-COII intergenic spacer exhibits greater variability than the protein-coding regions and may contain useful species-specific markers, its short length limits its phylogenetic utility.     

Complete mitochondrial genome of Tubulipora flabellaris (Bryozoa: Stenolaemata): the first representative from the class Stenolaemata with unique gene order

Mitochondrial genomes play a significant role in the reconstruction of phylogenetic relationships within metazoans. There are still many controversies concerning the phylogenetic position of the phylum Bryozoa. In this research, we have finished the complete mitochondrial genome of one bryozoan (Tubulipora flabellaris), which is the first representative from the class Stenolaemata. The complete mitochondrial genome of T. flabellaris is 13,763 bp in length and contains 36 genes, which lacks the atp8 gene in contrast to the typical metazoan mitochondrial genomes. Gene arrangement comparisons indicate that the mitochondrial genome of T. flabellaris has unique gene order when compared with other metazoans. The four known bryozoans complete mitochondrial genomes also have very different gene arrangements, indicates that bryozoan mitochondrial genomes have experienced drastic rearrangements. To investigate the phylogenetic relationship of Bryozoa, phylogenetic analyses based on amino acid sequences of 11 protein coding genes (excluding atp6 and atp8) from 26 metazoan complete mitochondrial genomes were made utilizing Maximum Likelihood (ML) and Bayesian methods, respectively. The results indicate the monopoly of Lophotrochozoa and a close relationship between Chaetognatha and Bryozoa. However, more evidences are needed to clarify the relationship between two groups. Lophophorate appeared to be polyphyletic according to our analyses. Meanwhile, neither analysis supports close relationship between Branchiopod and Phoronida. Four bryozoans form a clade and the relationship among them is T. flabellaris + (F. hispida + (B. neritina + W. subtorquata)), which is in coincidence with traditional classification system.     

The complete mitochondrial genome of the Sichuan taimen (Hucho bleekeri): repetitive sequences in the control region and phylogenetic implications for Salmonidae

The complete mitochondrial DNA genome of the Sichuan taimen (Hucho bleekeri) was determined by the long and accurate polymerase chain reaction (LA-PCR) and primer walking sequence method. The entire mitochondrial genome of this species is 16,997 bp in length, making it the longest among the completely sequenced Salmonidae mitochondrial genomes. It consists of two ribosomal RNA (rRNA) genes, 13 protein-coding genes, 22 transfer RNA (tRNA) genes, and one control region (CR). The gene arrangement, nucleotide composition, and codon usage pattern of the mitochondrial genome are similar to those of other teleosts. A T-type mononucleotide microsatellite and an 82 bp tandem repeat were identified in the control region, which were almost identical among the three H. bleekeri individuals examined. Both phylogenetic analyses based on 12 concatenated protein-coding genes of the heavy strand and on just the control region show that H. bleekeri is a basal species in Salmoninae. In addition, Salmo, Salvelinus and Oncorhynchus all represent monophyletic groups, respectively. All freshwater species occupied basal phylogenetic positions, and also possessed various tandem repeats in their mitochondrial control regions. These results support established phylogenetic relationships among genera in Salmonidae based on morphological and molecular analyses, and are consistent with the hypothesis that Salmonidae evolved from freshwater species.     

The Mitochondrial Genome of Xiphinema americanum sensu stricto (Nematoda: Enoplea): Considerable Economization in the Length and Structural Features of Encoded Genes

The complete sequence of the mitochondrial genome of the plant parasitic nematode Xiphinema americanum sensu stricto has been determined. At 12626bp it is the smallest metazoan mitochondrial genome reported to date. Genes are transcribed from both strands. Genes coding for 12 proteins, 2 rRNAs and 17 putative tRNAs (with the tRNA-C, I, N, S1, S2 missing) are predicted from the sequence. The arrangement of genes within the X. americanum mitochondrial genome is unique and includes gene overlaps. Comparisons with the mtDNA of other nematodes show that the small size of the X. americanum mtDNA is due to a combination of factors. The two mitochondrial rRNA genes are considerably smaller than those of other nematodes, with most of the protein encoding and tRNA genes also slightly smaller. In addition, five tRNAs genes are absent, lengthy noncoding regions are not present in the mtDNA, and several gene overlaps are present. [Reviewing Editor: Dr. Yues van de Peer] F. Lamberti: Deceased, 2004     

Highly conserved D-loop-like nuclear mitochondrial sequences (Numts) in tiger (Panthera tigris)

Using oligonucleotide primers designed to match hypervariable segments I (HVS-1) ofPanthera tigris mitochondrial DNA (mtDNA), we amplified two different PCR products (500 bp and 287 bp) in the tiger (Panthera tigris), but got only one PCR product (287 bp) in the leopard (Panthera pardus). Sequence analyses indicated that the sequence of 287 bp was a D-loop-like nuclear mitochondrial sequence (Numts), indicating a nuclear transfer that occurred approximately 4.8–17 million years ago in the tiger and 4.6–16 million years ago in the leopard. Although the mtDNA D-loop sequence has a rapid rate of evolution, the 287-bp Numts are highly conserved; they are nearly identical in tiger subspecies and only 1.742% different between tiger and leopard. Thus, such sequences represent molecular ‘fossils’ that can shed light on evolution of the mitochondrial genome and may be the most appropriate outgroup for phylogenetic analysis. This is also proved by comparing the phylogenetic trees reconstructed using the D-loop sequence of snow leopard and the 287-bp Numts as outgroup.     

Novel Repetitive Structures, Deviant Protein-Encoding Sequences andUnidentified ORFs in the Mitochondrial Genome of the BrachiopodLingula anatina

Complete sequence determination of the brachiopod Lingula anatina mtDNA (28,818 bp) revealed an organization that is remarkably atypical for an animal mt-genome. In addition to the usual set of 37 animal mitochondrial genes, which make up only 57% (16,555 bp) of the entire sequence, the genome contains lengthy unassigned sequences. All the genes are encoded in the same DNA strand, generally in a compact way, whereas the overall gene order is highly divergent in comparison with known animal mtDNA. Individual genes are generally longer and deviate considerably in sequence from their homologues in other animals. The genome contains two major repeat regions, in which 11 units of unassigned sequences and six genes (atp8, trnM, trnQ, trnV, and part of cox2 and nad2) are found in repetition, in the form of nested direct repeats of unparalleled complexity. One of the repeat regions contains unassigned repeat units dispersed among several unique sequences, novel repetitive structure for animal mtDNAs. Each of those unique sequences contains an open reading frame for a polypeptide between 80 and 357 amino acids long, potentially encoding a functional molecule, but none of them has been identified with known proteins. In both repeat regions, tRNA genes or tRNA gene-like sequences flank major repeated units, supporting the view that those structures play a role in the mitochondrial gene rearrangements. Although the intricate repeated organization of this genome can be explained by recurrent tandem duplications and subsequent deletions mediated by replication errors, other mechanisms, such as nonhomologous recombinations, appear to explain certain structures more easily.     

长爪栘[木衣]叶绿体基因组特征系统发育及密码子偏好性分析

长爪栘[木衣](Docynia longiunguis Q.Luo & J.L.Liu)是我国特有的栘[木衣]属植物,具有较高的食药用价值.对其叶绿体基因组进行分析,有助于阐明栘[木衣]属内的系统发育关系,为长爪栘[木衣]资源的开发利用及进一步研究奠定基础.结合其近缘种云南移[木衣]叶绿体基因组数据,在进行全序列比对后...     

Mitochondrial Genome of <Emphasis Type="Italic">Ciona savignyi</Emphasis> (Urochordata,Ascidiacea, Enterogona): Comparison of Gene Arrangement and tRNA Genes with <Emphasis Type="Italic">Halocynthia roretzi</Emphasis> Mitochondrial Genome

The complete nucleotide sequence of the urochordate Ciona savignyi (Ascidiacea, Enterogona) mitochondrial (mt) genome (14,737 bp) was determined. The Ciona mt genome does not encode a gene for ATP synthetase subunit 8 but encodes an additional tRNA^Gly gene (anticodon UCU), as is the case in another urochordate, Halocynthia roretzi (Ascidiacea, Pleurogona), mt genome. In addition, the Ciona mt genome encodes two tRNA^Met genes; anticodon CAT and anticodon TAT. The tRNA^Cys gene is thought to lack base pairs at the D-stem. Thus, the Ciona mt genome encodes 12 protein, 2 rRNA, and 24 tRNA genes. The gene arrangement of the Ciona mt genome differs greatly from those of any other metazoan mt genomes reported to date. Only three gene boundaries are shared between the Halocynthia and the Ciona mt genomes. Molecular phylogenetic analyses based on amino acid sequences of mt protein genes failed to demonstrate the monophyly of the chordates.     

Incongruence between primary sequence data and the distribution of a mitochondrial atp1 group II intron among ferns and horsetails

Using DNA sequence data from multiple genes (often from more than one genome compartment) to reconstruct phylogenetic relationships has become routine. Augmenting this approach with genomic structural characters (e.g., intron gain and loss, changes in gene order) as these data become available from comparative studies already has provided critical insight into some long-standing questions about the evolution of land plants. Here we report on the presence of a group II intron located in the mitochondrial atp1 gene of leptosporangiate and marattioid ferns. Primary sequence data for the atp1 gene are newly reported for 27 taxa, and results are presented from maximum likelihood-based phylogenetic analyses using Bayesian inference for 34 land plants in three data sets: (1) single-gene mitochondrial atp1 (exon+intron sequences); (2) five combined genes (mitochondrial atp1 [exon only]; plastid rbcL, atpB, rps4; nuclear SSU rDNA); and (3) same five combined genes plus morphology. All our phylogenetic analyses corroborate results from previous fern studies that used plastid and nuclear sequence data: the monophyly of euphyllophytes, as well as of monilophytes; whisk ferns (Psilotidae) sister to ophioglossoid ferns (Ophioglossidae); horsetails (Equisetopsida) sister to marattioid ferns (Marattiidae), which together are sister to the monophyletic leptosporangiate ferns. In contrast to the results from the primary sequence data, the genomic structural data (atp1 intron distribution pattern) would seem to suggest that leptosporangiate and marattioid ferns are monophyletic, and together they are the sister group to horsetails--a topology that is rarely reconstructed using primary sequence data.     

Complete DNA sequence of the mitochondrial genome ofCepaea nemoralis (Gastropoda: Pulmonata)

The nucleotide sequence of a mitochondrial genome of the pulmonate gastropod molluscCepaea nemoralis has been determined. Contained within the 14,100 basepairs (bp) are the two ribosomal RNA genes and 13 protein coding genes typical of metazoan mitochondrial genomes. TheCepaea mtDNA does contain a gene for ATPase subunit 8, like the clausiliid pulmonate,Albinaria, and the chiton,Katharina, but unlike the bivalve mollusc,Mytilus. The mitochondrial genetic code ofCepaea is proposed to be the same as that ofMytilus, Katharina, andDrosophila. Only 14 putative tRNA genes are presented, although there is sufficient unassigned sequence to encode the remainder of the expected total of 22 tRNA genes. These 14 tRNA genes are a mixture of standard cloverleaf structures and nonstandard structures containing TV replacement loops as seen in nematode and mosquito mitochondrial genomes. If the eight unidentified tRNA genes are indeed present, very little unassigned sequence would remain to serve as a control region. Genes are transcribed from both strands of the molecule. Base composition is the least biased for any reported animal mitochondrial genome and is also very little skewed between strands using measures independent of base composition. TheCepaea mitochondrial gene order is quite unlike that of any other reported metazoan mtDNA, with the exception of the recently reported partial sequences ofAlbinaria. No gene bound-aries are shared among all the reported molluscan taxa, demonstrating a complete lack of conservation of mitochondrial gene order across the phylum Mollusca.     

极小种群植物川柿叶绿体基因组特征与系统发育分析（附录）

川柿(Diospyros sutchuensis)为极小种群和国家重点保护野生植物,分布范围狭窄,种群数量极少。目前,川柿基因组信息缺乏,在柿属(Diospyros)中的系统亲缘关系不明确。该研究通过Illumina平台对川柿叶绿体基因组进行测序,应用Getorganellev1.7.3.4和PGA软件对基因组进行组装和注释,使用DnaSP6.12.03软件进行多序列对比分析,并使用REPuter、Tandem Reapeats Finder和MISA软件进行重复序列分析,使用CodonW1.4和EasyCodemL软件分别进行密码子偏好性和选择压力分析。同时,基于4个不同的叶绿体基因组序列数据集,使用IQtree软件分析川柿与11个柿属物种的系统发育关系。结果表明:(1)川柿叶绿体基因组全长157 917 bp,包含1对26 111 bp的反向重复区、大单拷贝区(87 303 bp)和小单拷贝区(18 392 bp),GC碱基含量为37.4%。(2)川柿叶绿体基因组共注释到113个基因,包括79个蛋白编码基因、30个tRNA基因和4个rRNA基因; 共检测到49个长重复序列、27个串联重复序列和34个简单重复序列; 蛋白编码基因中高频密码子31个,多数密码子末位碱基为A或U,编码亮氨酸的密码子使用最多; 基因组编码区比非编码区更为保守,10个高变热点区域可作为潜在的分子标记; 蛋白编码基因中有8个基因(ndhB、ndhG、ndhI、rbcL、rpoB、petB、petD、rps12)受到正选择压力。(3)系统发育分析显示,川柿与老鸦柿(D. rhombifolia)和乌柿(D. cathayensis)亲缘关系最为密切,它们与海南柿(D. hainanensis)共同形成一个单系分支。该研究结果既为川柿及柿属种质资源鉴定、遗传多样性保护以及种群恢复等提供了叶绿体基因组资源,也为阐明川柿的系统进化提供了重要的分子信息。     

Complete Mitochondrial Genomes of Carcinoscorpius rotundicauda and Tachypleus tridentatus (Xiphosura,Arthropoda) and Implications for Chelicerate Phylogenetic Studies

Horseshoe crabs (order Xiphosura) are often referred to as an ancient order of marine chelicerates and have been considered as keystone taxa for the understanding of chelicerate evolution. However, the mitochondrial genome of this order is only available from a single species, Limulus polyphemus. In the present study, we analyzed the complete mitochondrial genomes from two Asian horseshoe crabs, Carcinoscorpius rotundicauda and Tachypleus tridentatus to offer novel data for the evolutionary relationship within Xiphosura and their position in the chelicerate phylogeny. The mitochondrial genomes of C. rotundicauda (15,033 bp) and T. tridentatus (15,006 bp) encode 13 protein-coding genes, two ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes. Overall sequences and genome structure of two Asian species were highly similar to that of Limulus polyphemus, though clear differences among three were found in the stem-loop structure of the putative control region. In the phylogenetic analysis with complete mitochondrial genomes of 43 chelicerate species, C. rotundicauda and T. tridentatus were recovered as a monophyly, while L. polyphemus solely formed an independent clade. Xiphosuran species were placed at the basal root of the tree, and major other chelicerate taxa were clustered in a single monophyly, clearly confirming that horseshoe crabs composed an ancestral taxon among chelicerates. By contrast, the phylogenetic tree without the information of Asian horseshoe crabs did not support monophyletic clustering of other chelicerates. In conclusion, our analyses may provide more robust and reliable perspective on the study of evolutionary history for chelicerates than earlier analyses with a single Atlantic species.     

The Discriminatory Transfer Routes of tRNA Genes Among Organellar and Nuclear Genomes in Flowering Plants: A Genome-Wide Investigation of <Emphasis Type="Italic">indica</Emphasis> Rice

The transfer and integration of tRNA genes from organellar genomes to the nuclear genome and between organellar genomes occur extensively in flowering plants. The routes of the genetic materials flowing from one genome to another are biased, limited largely by compatibility of DNA replication and repair systems differing among the organelles and nucleus. After thoroughly surveying the tRNA gene transfer among organellar genomes and the nuclear genome of a domesticated rice (Oryza sativa L. ssp. indica), we found that (i) 15 mitochondrial tRNA genes originate from the plastid; (ii) 43 and 80 nuclear tRNA genes are mitochondrion-like and plastid-like, respectively; and (iii) 32 nuclear tRNA genes have both mitochondrial and plastid counterparts. Besides the native (or genuine) tRNA gene sets, the nuclear genome contains organelle-like tRNA genes that make up a complete set of tRNA species capable of transferring all amino acids. More than 97% of these organelle-like nuclear tRNA genes flank organelle-like sequences over 20 bp. Nearly 40% of them colocalize with two or more other organelle-like tRNA genes. Twelve of the 15 plastid-like mitochondrial tRNA genes possess 5′- and 3′-flanking sequences over 20 bp, and they are highly similar to their plastid counterparts. Phylogenetic analyses of the migrated tRNA genes and their original copies suggest that intergenomic tRNA gene transfer is an ongoing process with noticeable discriminatory routes among genomes in flowering plants. Electronic Supplementary Material Electronic Supplementary material is available for this article at and accessible for authorised users. Reviewing Editor: Dr. David Guttman     

A comparative analysis of the complete mitochondrial genome of the Eurasian otter <Emphasis Type="Italic">Lutra lutra</Emphasis> (Carnivora; Mustelidae)

Otter populations are declining throughout the world and most otter species are considered endangered. Molecular methods are suitable tools for population genetic research on endangered species. In the present study, we analyzed the complete mitochondrial genome (mitogenome) sequence of the Eurasian otter Lutra lutra. The mitochondrial DNA sequence of the Eurasian otter is 16,505 bp in length and consists of 13 protein-coding genes, 22 tRNAs, 2 rRNAs, and a control region (CR). The CR sequence of otters from Europe and Asia showed nearly identical numbers and nucleotide sequences of minisatellites. Phylogenetic analysis of Mustelidae mitogenomes, including individual genes, revealed that Lutrinae and Mustelinae form a clade, and that L. lutra and Enhydra lutris are sister taxa within the Lutrinae. Phylogenetic analyses revealed that of the 13 mitochondrial protein-coding genes, ND5 is the most reliable marker for analysis of phylogenetic relationships within the Mustelidae.     

Molecular phylogeny and systematics of the Lilium carniolicum group (Liliaceae) based on nuclear ITS sequences

The Lilium carniolicum group consists of several taxonomically dubious taxa endemic to the European flora. Internal transcribed spacer (ITS) sequences of nuclear ribosomal DNA (nrDNA) were used to clarify both the delineation of, and relationships among, taxa in the group as well as to provide insight on the phylogenetic position of the group within the genus. Maximum parsimony, maximum likelihood and Bayesian analyses were in general agreement, with all taxa in the group being very closely related, and the entire group being monophyletic. L. pyrenaicum and L. pomponium are placed at the basal position in the group, while L. chalcedonicum is shown to be more closely related to L. carniolicum than previously thought. Our analyses suggested that L. albanicum and L. jankae are distinct from L. carniolicum, while no evidence was found to support the same separation for L. bosniacum.     

食线虫真菌洛斯里被毛孢线粒体基因组的再分析

【目的】鉴定洛斯里被毛孢OWVT-1菌株的线粒体基因组,验证公布的USA-87-5菌株线粒体基因组中的错误,对洛斯里被毛孢正确的线粒体基因组序列进行注释并开展不同被毛孢物种间的比较线粒体基因组学分析。【方法】借助DNA高通量测序数据并通过必要的Sanger测序组装OWVT-1的线粒体基因组。通过PCR验证OWVT-1与公布的USA-87-5线粒体基因组序列差异的真实性。利用多种生物信息方法分析和注释洛斯里被毛孢的线粒体基因组。【结果】公布的洛斯里被毛孢USA-87-5菌株的线粒体基因组存在几处序列错误,包括3处长片段的插入缺失和多处短片段的插入缺失。实际上,洛斯里被毛孢USA-87-5与OWVT-1菌株的线粒体基因组序列完全相同。该菌的线粒体基因组全长62949 bp,在7个基因中共插入13个内含子,部分内含子和基因间区显现出序列退化的特征。洛斯里被毛孢、明尼苏达被毛孢、线虫被毛孢的线粒体基因组具有较强的共线性关系。除一些独立的ORF外,核心蛋白编码基因、rRNA基因和tRNA基因的排列顺序非常保守。基因间区的长短是影响3种被毛孢线粒体基因组大小最主要的因素。【结论】公布的洛斯里被毛孢USA-87-5菌株线粒体基因组中存在序列错误。本文新报道了OWVT-1菌株的线粒体基因组,并进行注释和比较线粒体基因组学分析。     

多花海棠叶绿体基因组分析

多花海棠（Malus floribunda Siebold.）是世界范围内广泛栽培的苹果属物种,具有较高的观赏价值和育种意义。对其进行叶绿体基因组比较分析,有利于完善苹果属系统进化以及种质利用的研究内容。基于全基因组测序数据,组装获得一个完整的具有四分体结构的多花海棠叶绿体基因组。该基因组包括大单拷贝区（88 142 bp）、反向重复区B （26 353 bp）、小单拷贝区（19 189 bp）与反向重复区A （26 353 bp）,共计160 037 bp。多花海棠叶绿体全基因组共注释到111个基因,包括78个蛋白编码基因、29个tRNA基因和4个rRNA基因。此外,在其基因组中识别到大量的重复序列,与三叶海棠和变叶海棠略有差异。通过计算相对同义密码子使用度,发现其高频密码子共30种,并且密码子具有偏向A/T结尾的使用模式。种间序列比对、边界分析的结果表明,大单拷贝区序列变异较大,8种苹果属植物SC区与IR区扩张收缩情况整体上较为相似。基于叶绿体基因组序列的系统进化分析,将多花海棠、湖北海棠和变叶海棠聚为一类。多花海棠叶绿体基因组的研究可为今后遗传标记开发与种质资源利用等提供数据支持。     

Evolution of fragmented mitochondrial ribosomal RNA genes inChlamydomonas

The fragmented mitochondrial ribosomal RNAs (rRNAs) of the green algaeChlamydomonas eugametos andChlamydomonas reinhardtii are discontinuously encoded in subgenic modules that are scrambled in order and interspersed with protein coding and tRNA genes. The mitochondrial rRNA genes of these two algae differ, however, in both the distribution and organization of rRNA coding information within their respective genomes. The objectives of this study were (1) to examine the phylogenetic relationships between the mitochondrial rRNA gene sequences ofC. eugametos andC. reinhardtii and those of the conventional mitochondrial rRNA genes of the green alga,Prototheca wickerhamii, and land plants and (2) to attempt to deduce the evolutionary pathways that gave rise to the unusual mitochondrial rRNA gene structures in the genusChlamydomonas. Although phylogenetic analysis revealed an affiliation between the mitochondrial rRNA gene sequences of the twoChlamydomonas taxa to the exclusion of all other mitochondrial rRNA gene sequences tested, no specific affiliation was noted between theChlamydomonas sequences andP. wickerhamii or land plants. Calculations of the minimal number of transpositions required to convert hypothetical ancestral rRNA gene organizations to the arrangements observed forC. eugametos andC. reinhardtii mitochondrial rRNA genes, as well as a limited survey of the size of mitochondrial rRNAs in other members of the genus, lead us to propose that the last common ancestor ofChlamydomonas algae contained fragmented mitochondrial rRNA genes that were nearly co-linear with conventional rRNA genes.