A deep-sea slant on the molecular phylogeny of the Scleractinia

Lophelia pertusa and Madrepora oculata are azooxanthellate corals with nearly cosmopolitan distributions. They form cold-water reefs in the upper bathyal zone on continental margins and offshore banks [A.D. Rogers, Int. Rev. Hydrobiol. 84 (1999) 315]. Lophelia is classified in the family Caryophylliidae and Madrepora in the family Oculinidae, both on the basis of skeletal morphology. Recent molecular studies of the scleractinians have given a new insight into the evolutionary history of this group. This study was aimed at clarifying the phylogenetic relationships of Lophelia and Madrepora, through the analysis of partial sequences of the mitochondrial 16S rDNA. Sequences were obtained for samples of L. pertusa collected in the northeast Atlantic and off Brazil, M. oculata, four other deep-sea and eight tropical coral species from the Réunion island in the Indian Ocean. The sequences were aligned with 69 homologous sequences of Scleractinia. Maximum parsimony and Bayesian analyses support previously published molecular topologies. The two specimens of L. pertusa grouped with two caryophylliids, confirming the existing classification of the species, but the large genetic distance between the two Lophelia samples suggests that these populations are genetically isolated from one another. M. oculata did not cluster with oculinids, but formed a monotypic clade lying between the families Pocilloporidae and Caryophyliidae. Phylogenetic analysis also suggested cryptic speciation within the tropical taxa Pocillopora meandriana and possibly Acropora humilis.     

Highly conserved gene arrangement of the mitochondrial genomes of 23 Plasmodium species

Mitochondrial (mt) genomes from diverse phylogenetic groups vary considerably in size, structure and organization. The genus Plasmodium, the causative agent of malaria, has the smallest mt genome in the form of a tandemly repeated, linear element of 6 kb. The Plasmodium mt genome encodes only three protein genes (cox1, cox3 and cob) and large- and small-subunit ribosomal RNA (rRNA) genes, which are highly fragmented with 19 identified rRNA pieces. The complete mt genome sequences of 21 Plasmodium species have been published but a thorough investigation of the arrangement of rRNA gene fragments has been undertaken for only Plasmodium falciparum, the human malaria parasite. In this study, we determined the arrangement of mt rRNA gene fragments in 23 Plasmodium species, including two newly determined mt genome sequences from P. gallinaceum and P. vinckei vinckei, as well as Leucocytozoon caulleryi, an outgroup of Plasmodium. Comparative analysis reveals complete conservation of the arrangement of rRNA gene fragments in the mt genomes of all the 23 Plasmodium species and L. caulleryi. Surveys for a new rRNA gene fragment using hidden Markov models enriched with recent mt genome sequences led us to suggest the mtR-26 sequence as a novel candidate LSU rRNA fragment in the mt genomes of the 24 species. Additionally, we found 22-25 bp-inverted repeat sequences, which may be involved in the generation of lineage-specific mt genome arrangements after divergence from a common ancestor of the genera Eimeria and Plasmodium/Leucocytozoon.     

The mitogenome of Paphia euglypta (Bivalvia: Veneridae) and comparative mitogenomic analyses of three venerids

Extraordinary variation has been found in mitochondrial (mt) genome inheritance, gene content and arrangement among bivalves. However, only few bivalve mt genomes have been comparatively analyzed to infer their evolutionary scenarios. In this study, the complete mt genome of the venerid Paphia euglypta (Bivalvia: Veneridae) was firstly studied and, secondly, it was comparatively analyzed with other venerids (e.g., Venerupis philippinarum and Meretrix petechialis) to better understand the mt genome evolution within a family. Though several common features such as the AT content, codon usage of protein-coding genes, and AT/GC skew are shared by the three venerids, a high level of variability is observed in genome size, gene content, gene order, arrangements and primary sequence of nucleotides or amino acids. Most of the gene rearrangement can be explained by the "tandem duplication and random loss" model. From the observed rearrangement patterns, we speculate that block interchange between adjacent genes may be common in the evolution of mt genomes in venerids. Furthermore, this study presents several new findings in mt genome annotation of V. philippinarum and M. petechialis, and hence we have reannotated the genome of these two species as: (1) the ORF of the formerly annotated cox2 gene in V. philippinarum is deduced by using a truncated "T" codon and a second cox2 gene is identified; (2) the trnS-AGN gene is identified and marked in the mt genome of both venerids. Thus, this study demonstrated a high variability of mt genomes in the Veneridae, and showed the importance of comparative mt genome analysis to interpret the evolution of the bivalve mt genome.     

Exploring the utility of an indel-rich, mitochondrial intergenic region as a molecular barcode for bamboo corals (Octocorallia: Isididae)

The DNA barcoding initiative has advocated the use of the 5'-end (～658bp) of mitochondrial (mt) cytochrome c oxidase subunit 1 (cox1) to genetically distinguish species. However, this has proven difficult within the subclass Octocorallia due to extraordinarily low substitution rates within mt protein-coding genes. Intergenic regions (IGRs), which have been little examined among octocorals, may be subject to high mutation rates and have proven useful target regions at both the interspecific and population levels of metazoans. Herein we examine a mt IGR (igr4) between the cytochrome b (cob) and NADH dehydrogenase subunit 6 (nad6) genes among species of the bamboo coral subfamily Keratoisidinae to evaluate its utility for barcoding and phylogenetic studies. Among 77 keratoisidin specimens, we found igr4 to vary in length between either 42bp (Acanella Gray, 1870 and Orstomisis Bayer, 1990) or 302-605bp (Isidella Gray, 1857, Lepidisis Verrill, 1883, Keratoisis Wright, 1869, and two undescribed genera). We interpreted the short igr4 sequence of Acanella eburnea (Pourtalès, 1868) as potentially indicative of additional mt genome-related novelties and thus sequenced its entire mt genome; gene content and gene order were the same as in a previously-sequenced bamboo coral mt genome. Alignment of the longer igr4 sequences included 108 parsimony-informative characters, as well as numerous indels ranging from 2-262bp in length. Uncorrected pairwise 'p' distances indicated sequence variation of 0-27.2%, as compared to 0-4.8% among the same specimens for the MutS homolog (msh1), currently the most widely sequenced octocorallian mt gene, and <0.4% for cox1 for a subset of the taxa. Despite the greater levels of variation, fewer unique haplotypes were observed at igr4 compared to msh1; however, in combination, the two gene regions revealed increased mt haplotype diversity relative to either gene region on their own.     

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.     

The Mitochondrial Genome of a Deep-Sea Bamboo Coral (Cnidaria,Anthozoa, Octocorallia,Isididae): Genome Structure and Putative Origins of Replication Are Not Conserved Among Octocorals

Octocoral mitochondrial (mt) DNA is subject to an exceptionally low rate of substitution, and it has been suggested that mt genome content and structure are conserved across the subclass, an observation that has been supported for most octocorallian families by phylogenetic analyses using PCR products spanning gene boundaries. However, failure to recover amplification products spanning the nad4L-msh1 gene junction in species from the family Isididae (bamboo corals) prompted us to sequence the complete mt genome of a deep-sea bamboo coral (undescribed species). Compared to the "typical" octocoral mt genome, which has 12 genes transcribed on one strand and 5 genes on the opposite (cox2, atp8, atp6, cox3, trnM), in the bamboo coral genome a contiguous string of 5 genes (msh1, rnl, nad2, nad5, nad4) has undergone an inversion, likely in a single event. Analyses of strand-specific compositional asymmetry suggest that (i) the light-strand origin of replication was also inverted and is adjacent to nad4, and (ii) the orientation of the heavy-strand origin of replication (OriH) has reversed relative to that of previously known octocoral mt genomes. Comparative analyses suggest that intramitochondrial recombination and errors in replication at OriH may be responsible for changes in gene order in octocorals and hexacorals, respectively. Using primers flanking the regions at either end of the inverted set of five genes, we examined closely related taxa and determined that the novel gene order is restricted to the deep-sea subfamily Keratoisidinae; however, we found no evidence for strand-specific mutational biases that may influence phylogenetic analyses that include this subfamily of bamboo corals.     

Making the most of mitochondrial genomes--markers for phylogeny, molecular ecology and barcodes in Schistosoma (Platyhelminthes: Digenea)

An increasing number of complete sequences of mitochondrial (mt) genomes provides the opportunity to optimise the choice of molecular markers for phylogenetic and ecological studies. This is particularly the case where mt genomes from closely related taxa have been sequenced; e.g., within Schistosoma. These blood flukes include species that are the causative agents of schistosomiasis, where there has been a need to optimise markers for species and strain recognition. For many phylogenetic and population genetic studies, the choice of nucleotide sequences depends primarily on suitable PCR primers. Complete mt genomes allow individual gene or other mt markers to be assessed relative to one another for potential information content, prior to broad-scale sampling. We assess the phylogenetic utility of individual genes and identify regions that contain the greatest interspecific variation for molecular ecological and diagnostic markers. We show that variable characters are not randomly distributed along the genome and there is a positive correlation between polymorphism and divergence. The mt genomes of African and Asian schistosomes were compared with the available intraspecific dataset of Schistosoma mansoni through sliding window analyses, in order to assess whether the observed polymorphism was at a level predicted from interspecific comparisons. We found a positive correlation except for the two genes (cox1 and nad1) adjoining the putative control region in S. mansoni. The genes nad1, nad4, nad5, cox1 and cox3 resolved phylogenies that were consistent with a benchmark phylogeny and in general, longer genes performed better in phylogenetic reconstruction. Considering the information content of entire mt genome sequences, partial cox1 would not be the ideal marker for either species identification (barcoding) or population studies with Schistosoma species. Instead, we suggest the use of cox3 and nad5 for both phylogenetic and population studies. Five primer pairs designed against Schistosoma mekongi and Schistosoma malayensis were tested successfully against Schistosoma japonicum. In combination, these fragments encompass 20-27% of the variation amongst the genomes (average total length approximately 14,000bp), thus providing an efficient means of encapsulating the greatest amount of variation within the shortest sequence. Comparative mitogenomics provides the basis of a rational approach to molecular marker selection and optimisation.     

Separate origins of group I introns in two mitochondrial genes of the katablepharid Leucocryptos marina

Mitochondria are descendants of the endosymbiotic α-proteobacterium most likely engulfed by the ancestral eukaryotic cells, and the proto-mitochondrial genome should have been severely streamlined in terms of both genome size and gene repertoire. In addition, mitochondrial (mt) sequence data indicated that frequent intron gain/loss events contributed to shaping the modern mt genome organizations, resulting in the homologous introns being shared between two distantly related mt genomes. Unfortunately, the bulk of mt sequence data currently available are of phylogenetically restricted lineages, i.e., metazoans, fungi, and land plants, and are insufficient to elucidate the entire picture of intron evolution in mt genomes. In this work, we sequenced a 12 kbp-fragment of the mt genome of the katablepharid Leucocryptos marina. Among nine protein-coding genes included in the mt genome fragment, the genes encoding cytochrome b and cytochrome c oxidase subunit I (cob and cox1) were interrupted by group I introns. We further identified that the cob and cox1 introns host open reading frames for homing endonucleases (HEs) belonging to distantly related superfamilies. Phylogenetic analyses recovered an affinity between the HE in the Leucocryptos cob intron and two green algal HEs, and that between the HE in the Leucocryptos cox1 intron and a fungal HE, suggesting that the Leucocryptos cob and cox1 introns possess distinct evolutionary origins. Although the current intron (and intronic HE) data are insufficient to infer how the homologous introns were distributed to distantly related mt genomes, the results presented here successfully expanded the evolutionary dynamism of group I introns in mt genomes.     

Mitochondrial genome reorganization characterizes various lineages of mesostigmatid mites (Acari: Parasitiformes)

Mesostigmata is an extremely diverse group of mites with more than 11,000 described species in 109 families. The complete mitochondrial (mt) genomes of five species of mesostigmatid mites from three families (Varroidae, Ologamasidae, Phytoseiidae) have been reported previously; all of them are rearranged or highly rearranged in gene order. However, it is unclear when mt genome reorganization occurred and how common it is in mesostigmatid mites. We sequenced the mt genomes of ten species of mesostigmatid mites from five more families (Blattisociidae, Diplogyniidae, Laelapidae, Macrochelidae, Parasitidae). We found that species in the families Diplogyniidae and Parasitidae have retained the ancestral mt genome organization of arthropods, which is in stark contrast to the highly rearranged mt genomes in the Phytoseiidae species. As in the Varroidae and Ologamasidae species, the mt genomes of the Blattisociidae, Macrochelidae and Laelapidae species are also rearranged but are less rearranged than in the Phytoseiidae species. Each of the six mesostigmatid families that have rearranged mt genomes is characterized by unique gene order not seen in other mesostigmatid families. Furthermore, the mt genome organization also differs among three genera of the Phytoseiidae, between two genera of the Laelapidae, and among three Macrocheles species of the Macrochelidae. Our results indicate that: (a) the most recent common ancestor of mesostigmatid mites likely retained the ancestral mt genome organization of arthropods; and (b) mt genome organization characterizes various lineages of mesostigmatid mites and provides a valuable source of information for understanding their phylogeny and evolution.     

Molecular variability in the Celleporella hyalina (Bryozoa; Cheilostomata) species complex: evidence for cryptic speciation from complete mitochondrial genomes

The bryozoan Celleporella has been shown to be composed of multiple, often cryptic, lineages. We sequenced two complete mitochondrial (mt) genomes of the Celleporella hyalina species complex from Wales, UK and Norway (i) to determine genetic divergence at the complete mt genome level, and (ii) to design new molecular markers for examining the interrelationships amongst the major lineages. In addressing (i), we estimated genetic divergence at three levels: (a) nucleotide diversity (π), (b) genome size, and (c) gene order. Genes nad4L, nad6, and atp8 showed the highest levels of divergence, and rrnL, rrnS, and cox1 showed the lowest levels. Inter-genome nucleotide divergence of protein-coding and ribosomal RNA genes, measured as π, was 0.21. The two genomes differed substantially in size, with the Norwegian genome being 2,573 base pairs (bp) longer than the Welsh genome, 17,265 and 14,692 bp, respectively. This difference in size is attributable to long non-coding regions present in the Norwegian genome. Both genomes exhibit similar gene orders, except for the translocation of one transfer RNA (trnA). Considering the high nucleotide diversity, genome size difference and change in gene order, these mt genomes are considered sufficiently divergent to have originated from two distinct species. In addressing (ii) we designed PCR primers that flank the most conserved regions of the genome: 1,300 bp of cox1 and a contiguous 2,000 bp fragment of rrnL + rrnS. The primers have yielded products for tissue from Wales, Norway, New Zealand, Alaska and Chile and should provide useful tools in establishing species- and population-level diversity within the Celleporella complex.     

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

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

Pyganodon (Bivalvia: Unionoida: Unionidae) phylogenetics: a male- and female-transmitted mitochondrial DNA perspective

Species boundaries, evolutionary relationships and geographic distributions of many unionoid bivalve species, like those in the genus Pyganodon, remain unresolved in Eastern North America. Because unionoid bivalves are one of the most imperiled groups of animals in the world, understanding the genetic variation within and among populations as well as among species is crucial for effective conservation planning. Conservation of unionoid species is indispensable from a freshwater habitat perspective but also because they possess a unique mitochondrial inheritance system where distinct gender-associated mitochondrial DNA lineages coexist: a female-transmitted (F) mt genome and a male-transmitted (M) mt genome that are involved in the maintenance of separate sexes (=dioecy). In this study, 42 populations of Pyganodon sp. were sampled across a large geographical range and fragments of two mitochondrial genes (cox1 and cox2) were sequenced from both the M- and F-transmitted mtDNA genomes. Our results support the recency of the divergence between P. cataracta and P. fragilis. We also found two relatively divergent F and M lineages within P. grandis. Surprisingly, the relationships among the P. grandis specimens in the F and M sequence trees are not congruent. We found that a single haplotype in P. lacustris has recently swept throughout the M genotype space leading to an unexpectedly low diversity in the M lineage in that species. Our survey put forward some challenging results that force us to rethink hybridization and species boundaries in the genus Pyganodon. As the M and F genomes do not always display the same phylogeographic story in each species, we also discuss the importance of being careful in the interpretation of molecular data based solely on maternal transmitted mtDNA genomes. The involvement of F and M genomes in unionoid bivalve sex determination likely played a role in the genesis of the unorthodox phylogeographic patterns reported herein.     

Mitogenomic analysis of Montipora cactus and Anacropora matthai (cnidaria; scleractinia; acroporidae) indicates an unequal rate of mitochondrial evolution among Acroporidae corals

The complete nucleotide sequence of the mitochondrial (mt) genome was determined for specimens of the coral species Montipora cactus (Bernard 1897) and Anacropora matthai (Pillai 1973), representing two morphologically distinct genera of the family Acroporidae. These sequences were compared with the published mt genome sequence for the confamilial species, Acropora tenuis (Dana 1846). The size of the mt genome was 17,887 bp and 17,888 bp for M. cactus and A. matthai. Gene content and organization was found to be very similar among the three Acroporidae mt genomes with a group I intron occurring in the NADH dehyrogenase 5 (nad5) gene. The intergenic regions were also similar in length among the three corals. The control region located between the small ribosomal RNA (ms) and the cytochrome oxidase 3 (cox3) gene was significantly smaller in M. cactus and A. matthai (both 627 bp) than in A. tenuis (1086 bp). Only one set of repeated sequences was identified at the 3′-end of the control regions in M. cactus and A. matthai. A lack of the abundant repetitive elements which have been reported for A. tenuis, accounts for the relatively short control regions in M. cactus and A. matthai. Pairwise distances and relative rate analyses of 13 protein coding genes, the group I intron and the largest intergenic region, igr3, revealed significant differences in the rate of molecular evolution of the mt genome among the three species, with an extremely slow rate being seen between Montipora and Anacropora. It is concluded that rapid mt genome evolution is taking place in genus Acropora relative to the confamilial genera Montipora and Anacropora although all are within the relatively slow range thought to be typical of Anthozoa.     

The Complete Maternally and Paternally Inherited Mitochondrial Genomes of the Endangered Freshwater Mussel Solenaia carinatus (Bivalvia: Unionidae) and Implications for Unionidae Taxonomy

Doubly uniparental inheritance (DUI) is an exception to the typical maternal inheritance of mitochondrial (mt) DNA in Metazoa, and found only in some bivalves. In species with DUI, there are two highly divergent gender-associated mt genomes: maternal (F) and paternal (M), which transmit independently and show different tissue localization. Solenaia carinatus is an endangered freshwater mussel species exclusive to Poyang Lake basin, China. Anthropogenic events in the watershed greatly threaten the survival of this species. Nevertheless, the taxonomy of S. carinatus based on shell morphology is confusing, and the subfamilial placement of the genus Solenaia remains unclear. In order to clarify the taxonomic status and discuss the phylogenetic implications of family Unionidae, the entire F and M mt genomes of S. carinatus were sequenced and compared with the mt genomes of diverse freshwater mussel species. The complete F and M mt genomes of S. carinatus are 16716 bp and 17102 bp in size, respectively. The F and M mt genomes of S. carinatus diverge by about 40% in nucleotide sequence and 48% in amino acid sequence. Compared to F counterparts, the M genome shows a more compact structure. Different gene arrangements are found in these two gender-associated mt genomes. Among these, the F genome cox2-rrnS gene order is considered to be a genome-level synapomorphy for female lineage of the subfamily Gonideinae. From maternal and paternal mtDNA perspectives, the phylogenetic analyses of Unionoida indicate that S. carinatus belongs to Gonideinae. The F and M clades in freshwater mussels are reciprocal monophyly. The phylogenetic trees advocate the classification of sampled Unionidae species into four subfamilies: Gonideinae, Ambleminae, Anodontinae, and Unioninae, which is supported by the morphological characteristics of glochidia.     

Characterization of the complete mitochondrial genomes of five Eimeria species from domestic chickens

Chicken coccidiosis caused by members of the genus Eimeria causes significant economic losses worldwide. In the present study we sequenced the complete mitochondrial DNA (mtDNA) sequences of six Eimeria species and analyzed features of their gene contents and genome organizations. The complete mt genomes of E. acervulina, E. brunetti, E. maxima, E. necatrix, E. tenella and E. praecox were 6179bp, 6148bp, 6169bp, 6214bp, 6213bp and 6174bp in size, respectively. All of the mt genomes consist of 3 genes for proteins (cox1, cox3, and cytb), 12 gene fragments for the large subunit (LSU) rRNA, and 7 gene fragments for the small subunit (SSU) rRNA, but no transfer RNA genes. The organization of the mt genomes is similar to that of Plasmodium, but distinct from Babesia and Theileria. The putative direction of translation for 3 genes (cox1, cox3, and cytb) was the same in all six Eimeria species. The contents of A+T of the mt genomes were 65.35% for E. acervulina, 65.43% for E. brunetti, 64.53% for E. maxima, 65.04% for E. necatrix, 64.98% for E. tenella and 65.59% for E. praecox. The AT bias has a significant effect on both the codon usage pattern and amino acid composition of proteins. Phylogenetic analyses using concatenated nucleotide sequences of the 2 protein-coding genes (cytb and cox1), with three different computational algorithms (Bayesian analysis, maximum parsimony and maximum likelihood), all revealed distinct groups with high statistical support, indicating that the six Eimeria spp. represent six distinct but closely-related species. These data provide novel mtDNA markers for studying the molecular epidemiology and population genetics of the six Eimeria spp., and should have implications for the molecular diagnosis, prevention and control of coccidiosis in domestic chickens.     

The mitochondrial genome of Moniliophthora roreri, the frosty pod rot pathogen of cacao

In this study, we report the sequence of the mitochondrial (mt) genome of the Basidiomycete fungus Moniliophthora roreri, which is the etiologic agent of frosty pod rot of cacao (Theobroma cacao L.). We also compare it to the mtDNA from the closely-related species Moniliophthora perniciosa, which causes witches' broom disease of cacao. The 94 Kb mtDNA genome of M. roreri has a circular topology and codes for the typical 14 mt genes involved in oxidative phosphorylation. It also codes for both rRNA genes, a ribosomal protein subunit, 13 intronic open reading frames (ORFs), and a full complement of 27 tRNA genes. The conserved genes of M. roreri mtDNA are completely syntenic with homologous genes of the 109 Kb mtDNA of M. perniciosa. As in M. perniciosa, M. roreri mtDNA contains a high number of hypothetical ORFs (28), a remarkable feature that make Moniliophthoras the largest reservoir of hypothetical ORFs among sequenced fungal mtDNA. Additionally, the mt genome of M. roreri has three free invertron-like linear mt plasmids, one of which is very similar to that previously described as integrated into the main M. perniciosa mtDNA molecule. Moniliophthora roreri mtDNA also has a region of suspected plasmid origin containing 15 hypothetical ORFs distributed in both strands. One of these ORFs is similar to an ORF in the mtDNA gene encoding DNA polymerase in Pleurotus ostreatus. The comparison to M. perniciosa showed that the 15 Kb difference in mtDNA sizes is mainly attributed to a lower abundance of repetitive regions in M. roreri (5.8 Kb vs 20.7 Kb). The most notable differences between M. roreri and M. perniciosa mtDNA are attributed to repeats and regions of plasmid origin. These elements might have contributed to the rapid evolution of mtDNA. Since M. roreri is the second species of the genus Moniliophthora whose mtDNA genome has been sequenced, the data presented here contribute valuable information for understanding the evolution of fungal mt genomes among closely-related species.     

啮总目昆虫的线粒体基因组多样性及系统发育研究进展

啮总目包括啮虫目（皮虱和书虱）和虱目（羽虱和吸虱）,是农业和医学等领域具有重要经济意义和研究价值的类群,目前已鉴定和描述的物种超过10 000个。啮总目昆虫线粒体基因组的变异性在昆虫各类群中最为剧烈,这些变异包括基因组的结构、基因排序、基因含量和链上分布等诸多方面。本文全面分析和总结了啮总目昆虫裂化线粒体基因组的进化属性,并结合两侧对称动物线粒体基因组的裂化特征重构了线粒体基因组环裂化的过程。引入“线粒体基因组核型”的概念来描述动物线粒体基因组丰富的变异程度。动物线粒体的染色体有减小的趋势,而线粒体基因组的裂化正是体现这种趋势的一种重要策略。同时,总结和探讨了目前具有争议的啮总目主要类群间的系统发育关系。本综述为啮总目昆虫线粒体基因组学、啮总目系统发生关系以及两侧对称动物线粒体基因组进化模式的研究提供一个新的视角。     

Heteroplasmy in the Mitochondrial Genomes of Human Lice and Ticks Revealed by High Throughput Sequencing

The typical mitochondrial (mt) genomes of bilateral animals consist of 37 genes on a single circular chromosome. The mt genomes of the human body louse, Pediculus humanus, and the human head louse, Pediculus capitis, however, are extensively fragmented and contain 20 minichromosomes, with one to three genes on each minichromosome. Heteroplasmy, i.e. nucleotide polymorphisms in the mt genome within individuals, has been shown to be significantly higher in the mt cox1 gene of human lice than in humans and other animals that have the typical mt genomes. To understand whether the extent of heteroplasmy in human lice is associated with mt genome fragmentation, we sequenced the entire coding regions of all of the mt minichromosomes of six human body lice and six human head lice from Ethiopia, China and France with an Illumina HiSeq platform. For comparison, we also sequenced the entire coding regions of the mt genomes of seven species of ticks, which have the typical mitochondrial genome organization of bilateral animals. We found that the level of heteroplasmy varies significantly both among the human lice and among the ticks. The human lice from Ethiopia have significantly higher level of heteroplasmy than those from China and France (P_t<0.05). The tick, Amblyomma cajennense, has significantly higher level of heteroplasmy than other ticks (P_t<0.05). Our results indicate that heteroplasmy level can be substantially variable within a species and among closely related species, and does not appear to be determined by single factors such as genome fragmentation.     

Complete nucleotide sequence of the mitochondrial genome of a Malagasy poison frog Mantella madagascariensis: evolutionary implications on mitochondrial genomes of higher anuran groups

We determined the complete nucleotide sequence of the mitochondrial (mt) genome of a Malagasy poison frog, Mantella madagascariensis (family Mantellidae), and partial sequences of two Mantella (M. baroni and M. bernhardi) and two additional mantellid species (Boophis madagascariensis and Mantidactylus cf. ulcerosus). The M. madagascariensis genome was shown to be the largest (23kbp) of all vertebrate mtDNAs investigated so far. Furthermore, the following unique features were revealed: (1) the positions of some genes and gene regions were rearranged compared to mitochondrial genomes typical for vertebrates and other anuran groups, (2) two distinct genes and a pseudogene corresponding to transfer RNA gene for methionine (tRNA-Met) were encoded, and (3) two control regions with very high sequence homology were present. These features were shared by the two other Mantella species but not the other mantellid species, indicating dynamic genome reorganization in a common ancestor linage before divergence of the Mantella genus. The reorganization pathway could be explained by a model of gene duplication and deletion. Duplication and deletion events also seem to have been responsible for concerted sequence evolution of the control regions in Mantella mt genomes. It is also suggested that the pseudo tRNA-Met gene sustained for a long time in Mantella mt genomes possibly functions as a punctuation marker for NADH dehydrogenase subunit (ND) 2 mRNA processing. Phylogenetic analyses employing a large sequence data set of mt genes supported the monophyly of Mantellidae and Rhacophoridae and other recent phylogenetic views for ranoid frogs. The resultant phylogenetic relationship also suggested parallel occurrence of two tRNA-Met genes, duplicated control regions, and ND5 gene translocation in independent ranoid lineages.