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.     

Limitations of mitochondrial gene barcoding in Octocorallia

The widespread assumption that COI and other mitochondrial genes will be ineffective DNA barcodes for anthozoan cnidarians has not been well tested for most anthozoans other than scleractinian corals. Here we examine the limitations of mitochondrial gene barcoding in the sub-class Octocorallia, a large, diverse, and ecologically important group of anthozoans. Pairwise genetic distance values (uncorrected p) were compared for three candidate barcoding regions: the Folmer region of COI; a fragment of the octocoral-specific mitochondrial protein-coding gene, msh1; and an extended barcode of msh1 plus COI with a short, adjacent intergenic region (igr1). Intraspecific variation was <0.5%, with most species exhibiting no variation in any of the three gene regions. Interspecific divergence was also low: 18.5% of congeneric morphospecies shared identical COI barcodes, and there was no discernible barcoding gap between intra- and interspecific p values. In a case study to assess regional octocoral biodiversity, COI and msh1 barcodes each identified 70% of morphospecies. In a second case study, a nucleotide character-based analysis correctly identified 70% of species in the temperate genus Alcyonium. Although interspecific genetic distances were 2× greater for msh1 than COI, each marker identified similar numbers of species in the two case studies, and the extended COI + igr1 + msh1 barcode more effectively discriminated sister taxa in Alcyonium. Although far from perfect for species identification, a COI + igr1 + msh1 barcode nonetheless represents a valuable addition to the depauperate set of characters available for octocoral taxonomy.     

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.     

Complete mitochondrial genomes of two Japanese precious corals, Paracorallium japonicum and Corallium konojoi (Cnidaria, Octocorallia, Coralliidae): notable differences in gene arrangement

Precious coral are taxonomically a group of corals that belong to the family Coralliidae within the order Alcyonacea, subclass Octocorallia, and class Anthozoa, whose skeletal axes are used for jewelry. They are distributed in the Mediterranean Sea and in waters adjacent to Japan, Taiwan, Midway Island and the Hawaiian Islands. The genus Corallium of the family Coralliidae was recently divided into two genera, Corallium and Paracorallium, based on morphological observations, but insufficient molecular evidence to support this classification has been presented to date. We determined for the first time the complete mitochondrial genome sequence of two precious corals P. japonicum and C. konojoi, in order to clarify their systematic positions. The circular mitochondrial genomes of P. japonicum and C. konojoi are 18,913bp and 18,969bp in length, respectively, and encode 13 typical energy pathway protein coding genes (nad1-6, nad4L, cox1-3, cob, atp6 and atp8), two ribosomal RNA genes (rns and rnl), a transfer RNA (trnM) and a mismatch repair gene homologue msh1. The two genomes have an overall nucleotide sequence identity of 97.5%, which is comparable to that between Acanella eburnea and Keratoisidinae sp. belonging to Octocorallia. Surprisingly, however, their gene arrangements were not identical. Phylogenetic analyses using seven complete mitochondrial genome sequences belonging to species in the subclass Octocorallia indicated that within the subclass, at least three gene order rearrangement events occurred during evolution. Our results support the validity of the morphological classification that separated the family Coralliidae into two genera, Corallium and Paracorallium.     

Assessing the use of the mitochondrial cox1 marker for use in DNA barcoding of red algae (Rhodophyta)

The red algae, a remarkably diverse group of organisms, are difficult to identify using morphology alone. Following the proposal to use the mitochondrial cytochrome c oxidase subunit I (cox1) for DNA barcoding animals, we assessed the use of this gene in the identification of red algae using 48 samples plus 31 sequences obtained from GenBank. The data set spanned six orders of red algae: the Bangiales, Ceramiales, Corallinales, Gigartinales, Gracilariales and Rhodymeniales. The results indicated that species could be discriminated. Intraspecific variation was between 0 and 4 bp over 539 bp analyzed except in Mastocarpus stellatus (0-14 bp) and Gracilaria gracilis (0-11 bp). Cryptic diversity was found in Bangia fuscopurpurea, Corallina officinalis, G. gracilis, M. stellatus, Porphyra leucosticta and P. umbilicalis. Interspecific variation across all taxa was between 28 and 148 bp, except for G. gracilis and M. stellatus. A comparison of cox1 with the plastid Rubisco spacer for Porphyra species revealed that it was a more sensitive marker in revealing incipient speciation and cryptic diversity. The cox1 gene has the potential to be used for DNA barcoding of red algae, although a good taxonomic foundation coupled with extensive sampling of taxa is essential for the development of an effective identification system.     

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.     

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.     

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.     

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.     

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.     

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.     

Aborjinia corallicola sp. n., a new nematode species (Nematoda: Marimermithidae) associated with the bamboo coral Acanella arbuscula (Johnson)

Systematic Parasitology - A new species of Aborjinia Özdikmen, 2010 is described from the tissues of the cold-water bamboo coral Acanella arbuscula (Johnson) from the northwest Atlantic....     

Novel organization of the mitochondrial genome in the deep-sea coral, Madrepora oculata (Hexacorallia, Scleractinia, Oculinidae) and its taxonomic implications

Madrepora is one of the most ecologically important genera of reef-building scleractinians in the deep sea, occurring from tropical to high-latitude regions. Despite this, the taxonomic affinities and relationships within the genus Madrepora remain unclear. To clarify these issues, we sequenced the mitochondrial (mt) genome of the most widespread Madrepora species, M. oculata, and compared this with data for other scleractinians. The architecture of the M. oculata mt genome was very similar to that of other scleractinians, except for a novel gene rearrangement affecting only cox2 and cox3. This pattern of gene organization was common to four geographically distinct M. oculata individuals as well as the congeneric species M. minutiseptum, but was not shared by other genera that are closely related on the basis of cox1 sequence analysis nor other oculinids, suggesting that it might be unique to Madrepora.     

The complete mitochondrial genome of the verongid sponge Aplysina cauliformis: implications for DNA barcoding in demosponges

DNA “barcoding,” the determination of taxon-specific genetic variation typically within a fragment of the mitochondrial cytochrome oxidase 1 (cox1) gene, has emerged as a useful complement to morphological studies, and is routinely used by expert taxonomists to identify cryptic species and by non-experts to better identify samples collected during field surveys. The rate of molecular evolution in the mitochondrial genomes (mtDNA) of nonbilaterian animals (sponges, cnidarians, and placozoans) is much slower than in bilaterian animals for which DNA barcoding strategies were developed. If sequence divergence among nonbilaterian mtDNA and specifically cox1 is too slow to generate diagnostic variation, alternative genes for DNA barcoding and species-level phylogenies should be considered. Previous study across the Aplysinidae (Demospongiae, Verongida) family of sponges demonstrated no nucleotide substitutions in the traditional cox1 barcoding fragment among the Caribbean species of Aplysina. As the mitochondrial genome of Aplysina fulva has previously been sequenced, we are now able to make the first comparisons between complete mtDNA of congeneric demosponges to assess whether potentially informative variation exists in genes other than cox1. In this article, we present the complete mitochondrial genome of Aplysina cauliformis, a circular molecule 19620 bp in size. The mitochondrial genome of A. cauliformis is the same length as is A. fulva and shows six confirmed nucleotide differences and an additional 11 potential SNPs. Of the six confirmed SNPs, NADH dehydrogenase subunit 5 (nad5) and nad2 each contain two, and in nad2 both yield amino acid substitutions, suggesting balancing selection may act on this gene. Thus, while the low nucleotide diversity in Caribbean aplysinid cox1 extends to the entire mitochondrial genome, some genes do display variation. If these represent interspecific differences, then they may be useful alternative markers for studies in recently diverged sponge clades.     

Deciphering amphibian diversity through DNA barcoding: chances and challenges

Amphibians globally are in decline, yet there is still a tremendous amount of unrecognized diversity, calling for an acceleration of taxonomic exploration. This process will be greatly facilitated by a DNA barcoding system; however, the mitochondrial population structure of many amphibian species presents numerous challenges to such a standardized, single locus, approach. Here we analyse intra- and interspecific patterns of mitochondrial variation in two distantly related groups of amphibians, mantellid frogs and salamanders, to determine the promise of DNA barcoding with cytochrome oxidase subunit I (cox1) sequences in this taxon. High intraspecific cox1 divergences of 7-14% were observed (18% in one case) within the whole set of amphibian sequences analysed. These high values are not caused by particularly high substitution rates of this gene but by generally deep mitochondrial divergences within and among amphibian species. Despite these high divergences, cox1 sequences were able to correctly identify species including disparate geographic variants. The main problems with cox1 barcoding of amphibians are (i) the high variability of priming sites that hinder the application of universal primers to all species and (ii) the observed distinct overlap of intraspecific and interspecific divergence values, which implies difficulties in the definition of threshold values to identify candidate species. Common discordances between geographical signatures of mitochondrial and nuclear markers in amphibians indicate that a single-locus approach can be problematic when high accuracy of DNA barcoding is required. We suggest that a number of mitochondrial and nuclear genes may be used as DNA barcoding markers to complement cox1.     

Characterization of the complete mitochondrial genome of Bombyx mori strain H9 (Lepidoptera: Bombycidae)

The complete mitochondrial genome (mitogenome) of Bombyx mori strain H9 (Lepidoptera: Bombycidae) is 15,670 base pairs (bp) in length, encoding 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and a control region. The nucleotide composition of the genome is highly A + T biased, accounting for 81.31%, with a slightly positive AT skewness (0.059). The arrangement of 13 PCGs is similar to that of other sequenced lepidopterans. All the PCGs are initiated by ATN codons, except for the cytochrome c oxidase subunit 1 (cox1) gene, which is proposed by the TTAG sequence as observed in other lepidopterans. Unlike the other PCGs, the cox1 and cytochrome c oxidase subunit 2 (cox2) genes have incomplete stop codons consisting of just a T. All tRNAs have typical structures of insect mitochondrial tRNAs, which is different from other sequenced lepidopterans. The structure of A + T-rich region is similar to that of other sequenced lepidopterans, including non-repetitive sequences, the ATAGA binding domain, a 18 bp poly-T stretch and a poly-A element upstream of transfer RNA M (trnM) gene. Phylogenetic analysis shows that the domesticated silkmoth B. mori originated from the Chinese Bombyx mandarina.