Two circular chromosomes of unequal copy number make up the mitochondrial genome of the rotifer Brachionus plicatilis

Multigenerational resilience of the marine rotifer Brachionus plicatilis to high temperature after additive exposure to high salinity and nanoplastics

To study multigenerational resilience to high temperature (HT) conditions, we exposed Brachionus plicatilis marine rotifers to HT, high salinity (HS), and nanoplastics (NPs), and measured reproductive and life-cycle endpoints. After exposure to HT, rotifer lifespans were reduced, but daily production of offspring increased. However, both combined HT/HS and HT/HS/NP exposure led to additional decreases in longevity and reproductive ability; the antioxidant defense mechanisms of the rotifers were also notably upregulated as measured by reactive oxygen species levels. Fatty-acid profiles were reduced in all conditions. In multigenerational experiments, the negative effects of HT dissipated rapidly; however, the effects of HT/HS and HT/HS/NPs required four generations to disappear completely. The findings indicated that B. plicatilis were able to recover from these environmental stressors. This study demonstrated the resilience of aquatic organisms in response to changing environmental conditions and provides insights into the complex interactions of different abiotic stressors.

Single and combined effects of increased temperature and methylmercury on different stages of the marine rotifer Brachionus plicatilis

Rapid, anthropogenic activity-induced global warming is a severe problem that not only raises water temperatures but also shifts aquatic environments by increasing the bioavailability of heavy metals (HMs), with potentially complicated effects on aquatic organisms, including small aquatic invertebrates. For this paper, we investigated the combined effects of temperature (23 and 28 °C) and methylmercury (MeHg) by measuring physiological changes, bioaccumulation, oxidative stress, antioxidants, and the mitogen-activated protein kinase signaling pathway in the marine rotifer Brachionus plicatilis. High temperature and MeHg adversely affected the survival rate, lifespan, and population of rotifers, and bioaccumulation, oxidative stress, and biochemical reactions depended on the developmental stage, with neonates showing higher susceptibility than adults. These findings demonstrate that increased temperature enhances potentially toxic effects from MeHg, and susceptibility differs with the developmental stage. This study provides a comprehensive understanding of the combined effects of elevated temperature and MeHg on rotifers. ENVIRONMENTAL IMPLICATION: Methylmercury (MeHg) is a widespread and harmful heavy metal that can induce lethal effects on aquatic organisms in even trace amounts. The toxicity of metals can vary depending on various environmental conditions. In particular, rising temperatures are considered a major factor affecting bioavailability and toxicity by changing the sensitivity of organisms. However, there are few studies on the combinational effects of high temperatures and MeHg on aquatic animals, especially invertebrates. Our research would contribute to understanding the actual responses of aquatic organisms to complex aquatic environments.

Two-Fold ND5 Genes, Three-Fold Control Regions, lncRNA, and the "Missing" ATP8 Found in the Mitogenomes of Polypedates megacephalus (Rhacophridae: Polypedates)

In prior research on the mitochondrial genome (mitogenome) of Polypedates megacephalus, the one copy of ND5 gene was translocated to the control region (CR) and the ATP8 gene was not found. Gene loss is uncommon among vertebrates. However, in this study, we resequenced the mitogenomes of P. megacephalus from different regions using a "primer bridging" approach with Sanger sequencing technologies, which revealed the "missing" ATP8 gene in P. megacephalus as well as three other previously published Polypedates. The mitogenome of this species was found to contain two copies of the ND5 genes and three copies of the control regions. Furthermore, multiple tandem repeats were identified in the control regions. Notably, we observed that there was no correlation between genetic divergence and geographic distance. However, using the mitogenome, gene expression analysis was performed via RT-qPCR of liver samples and it was thus determined that COIII, ND2, ND4, and ND6 were reduced to 0.64 ± 0.24, 0.55 ± 0.34, 0.44 ± 0.21 and 0.65 ± 0.17, respectively, under low-temperature stress (8 °C) as compared with controls (p < 0.05). Remarkably, the transcript of long non-coding RNA (lncRNA) between positions 8029 and 8612 decreased significantly with exposure to low-temperature stress (8 °C). Antisense ND6 gene expression showed a downward trend, but this was not significant. These results reveal that modulations of protein-coding mitochondrial genes and lncRNAs of P. megacephalus play a crucial role in the molecular response to cold stress.

Physiological effects and molecular response in the marine rotifer Brachionus plicatilis after combined exposure to nanoplastics and copper

Because nanoplastics (NPs) can transport pollutants, the absorption of surrounding pollutants into NPs and their effects are important environmental issues. This study shows a combined effect of high concentrations of NPs and copper (Cu) in the marine rotifer Brachionus plicatilis. Co-exposure decreased the growth rate, reproduction, and lifespan. The highest level of NP ingestion was detected in the co-treated group, but the Cu concentration was higher in the Cu single-exposure group. ERK activation played a key role in the downstream cell signaling pathway activated by the interaction of NPs and Cu. The increased sensitivity of B. plicatilis to Cu could be due to the impairment of MXR function caused by a high concentration of NPs, which supports our in vivo experiment results. Our results show that exposure to NPs could induce the dysfunction of several critical molecular responses, weakening resistance to Cu and thereby increasing its physiological toxicity in B. plicatilis.

Origin of minicircular mitochondrial genomes in red algae

Eukaryotic organelle genomes are generally of conserved size and gene content within phylogenetic groups. However, significant variation in genome structure may occur. Here, we report that the Stylonematophyceae red algae contain multipartite circular mitochondrial genomes (i.e., minicircles) which encode one or two genes bounded by a specific cassette and a conserved constant region. These minicircles are visualized using fluorescence microscope and scanning electron microscope, proving the circularity. Mitochondrial gene sets are reduced in these highly divergent mitogenomes. Newly generated chromosome-level nuclear genome assembly of Rhodosorus marinus reveals that most mitochondrial ribosomal subunit genes are transferred to the nuclear genome. Hetero-concatemers that resulted from recombination between minicircles and unique gene inventory that is responsible for mitochondrial genome stability may explain how the transition from typical mitochondrial genome to minicircles occurs. Our results offer inspiration on minicircular organelle genome formation and highlight an extreme case of mitochondrial gene inventory reduction.

Variation in heat shock protein 40 kDa relates to divergence in thermotolerance among cryptic rotifer species

Genetic divergence and the frequency of hybridization are central for defining species delimitations, especially among cryptic species where morphological differences are merely absent. Rotifers are known for their high cryptic diversity and therefore are ideal model organisms to investigate such patterns. Here, we used the recently resolved Brachionus calyciflorus species complex to investigate whether previously observed between species differences in thermotolerance and gene expression are also reflected in their genomic footprint. We identified a Heat Shock Protein gene (HSP 40 kDa) which exhibits cross species pronounced sequence variation. This gene exhibits species-specific fixed sites, alleles, and sites putatively under positive selection. These sites are located in protein binding regions involved in chaperoning and may therefore reflect adaptive diversification. By comparing three genetic markers (ITS, COI, HSP 40 kDa), we revealed hybridization events between the cryptic species. The low frequency of introgressive haplotypes/alleles suggest a tight, but not fully impermeable boundary between the cryptic species.

Complete Mitogenomes of Polypedates Tree Frogs Unveil Gene Rearrangement and Concerted Evolution within Rhacophoridae

Simple Summary

Duplicated control regions have been reported several times in the tree frog family Rhacophoridae, and previous studies have mostly relied on sequence analysis to reconstruct their evolution. This is the first study to employ a phylogenetic method to demonstrate the existence of concerted and parallel evolution succinctly and intuitively in the duplicated control regions of the family Rhacophoridae. Phylogenetic relationships were also used to illustrate the parallel evolution of ATP8 loss of function in the genus Polypedates. In general, this study elucidated the evolutionary patterns and pathways of mitochondrial gene rearrangement of the family Rhacophoridae from a phylogenetic perspective, which aids in understanding the evolutionary history of this fascinating tree frog taxon from a molecular evolution standpoint.

Abstract

New developments in sequencing technology and nucleotide analysis have allowed us to make great advances in reconstructing anuran phylogeny. As a clade of representative amphibians that have radiated from aquatic to arboreal habitats, our understanding of the systematic status and molecular biology of rhacophorid tree frogs is still limited. We determined two new mitogenomes for the genus Polypedates (Rhacophoridae): P. impresus and P. mutus. We conducted comparative and phylogenetic analyses using our data and seven other rhacophorid mitogenomes. The mitogenomes of the genera Polypedates, Buergeria, and Zhangixalus were almost identical, except that the ATP8 gene in Polypedates had become a non-coding region; Buergeria maintained the legacy “LTPF” tRNA gene cluster compared to the novel “TLPF” order in the other two genera; and B. buergeri and Z. dennysi had no control region (CR) duplication. The resulting phylogenetic relationship supporting the above gene rearrangement pathway suggested parallel evolution of ATP8 gene loss of function (LoF) in Polypedates and CR duplication with concerted evolution of paralogous CRs in rhacophorids. Finally, conflicting topologies in the phylograms of 185 species reflected the advantages of phylogenetic analyses using multiple loci.

Oxidative stress-mediated synergistic deleterious effects of nano- and microplastics in the hypoxia-conditioned marine rotifer Brachionus plicatilis

While pollution due to nano- and micro-sized plastics (NMPs) and hypoxic conditions both occur in coastal areas, the deleterious potential of co-exposure to hypoxia and NMPs (hypoxia and micro-sized plastics, HMPs; hypoxia and nano-sized plastics, HNPs) is largely unclear. Here, we provide evidence for multigenerational effects of HMP and HNP in the marine rotifer Brachionus plicatilis by investigating changes in its life traits, antioxidant system, and hypoxia-inducible factor (HIF) pathway using an orthogonal experimental design, with nanoscale and microscale particles measuring 0.05 μm and 6.0 μm in diameter, respectively, and hypoxic conditions of 0.5 mg/L for six generations. Combined exposure to NMPs and hypoxia caused a significant decrease in fecundity and overproduction of reactive oxygen species (ROS). The HIF pathway and circadian clock genes were also significantly upregulated in response to HMP and HNP exposure. In particular, synergistic deleterious effects of HNP were evident, suggesting that size-dependent toxicity can be a major driver of the effects of hypoxia and NMP co-exposure. After several generations of exposure, ROS levels returned to basal levels and transcriptomic resilience was observed, although rotifer reproduction remained suppressed. These findings help eluciating the underlying molecular mechanisms involved in responses to plastic pollution in hypoxic conditions.

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome’s size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

Comparative analysis reveals within-population genome size variation in a rotifer is driven by large genomic elements with highly abundant satellite DNA repeat elements

BACKGROUND

Eukaryotic genomes are known to display an enormous variation in size, but the evolutionary causes of this phenomenon are still poorly understood. To obtain mechanistic insights into such variation, previous studies have often employed comparative genomics approaches involving closely related species or geographically isolated populations within a species. Genome comparisons among individuals of the same population remained so far understudied-despite their great potential in providing a microevolutionary perspective to genome size evolution. The rotifer Brachionus asplanchnoidis represents one of the most extreme cases of within-population genome size variation among eukaryotes, displaying almost twofold variation within a geographic population.

RESULTS

Here, we used a whole-genome sequencing approach to identify the underlying DNA sequence differences by assembling a high-quality reference genome draft for one individual of the population and aligning short reads of 15 individuals from the same geographic population including the reference individual. We identified several large, contiguous copy number variable regions (CNVs), up to megabases in size, which exhibited striking coverage differences among individuals, and whose coverage overall scaled with genome size. CNVs were of remarkably low complexity, being mainly composed of tandemly repeated satellite DNA with only a few interspersed genes or other sequences, and were characterized by a significantly elevated GC-content. CNV patterns in offspring of two parents with divergent genome size and CNV patterns in several individuals from an inbred line differing in genome size demonstrated inheritance and accumulation of CNVs across generations.

CONCLUSIONS

By identifying the exact genomic elements that cause within-population genome size variation, our study paves the way for studying genome size evolution in contemporary populations rather than inferring patterns and processes a posteriori from species comparisons.

Observations of linear aggregation behavior in rotifers (Brachionus calyciflorus)

Linear aggregation is present in some animals, such as the coordinated movement of ants and the migration of caterpillars and spinylobsters, but none has been reported on rotifers. The rotifers were collected and clone cultured in the laboratory at 25 ± 1°C, under natural light (light intensity ~130 lx, L:D = 14:10). The culture medium(pH = 7.3) was formulated as described by Suga et al., and rotifers were fed on the micro algae Scenedesmus obliquus grown in HB-4 medium to the exponential growth stage. When density was high (150 individuals ml-1), the behavior of rotifers was observed using a stereo microscope (Motic ES-18TZLED). In this paper, linear aggregation in Brachionus calyciflorus was found for the first time, and experiments were carried out to verify the correlation between linear aggregation and culture density of B. calyciflorus. With the increase of density, the number of aggregations increase, the number of individuals in the aggregation increased, and the maintenance time of the aggregation was also increased. Therefore, we speculate that the formation of aggregates is related to density and may be a behavioral signal of density increase, which may transmit information between density increase and formation of dormant eggs.

Mitochondrial Genomic Landscape: A Portrait of the Mitochondrial Genome 40 Years after the First Complete Sequence

Notwithstanding the initial claims of general conservation, mitochondrial genomes are a largely heterogeneous set of organellar chromosomes which displays a bewildering diversity in terms of structure, architecture, gene content, and functionality. The mitochondrial genome is typically described as a single chromosome, yet many examples of multipartite genomes have been found (for example, among sponges and diplonemeans); the mitochondrial genome is typically depicted as circular, yet many linear genomes are known (for example, among jellyfish, alveolates, and apicomplexans); the chromosome is normally said to be “small”, yet there is a huge variation between the smallest and the largest known genomes (found, for example, in ctenophores and vascular plants, respectively); even the gene content is highly unconserved, ranging from the 13 oxidative phosphorylation-related enzymatic subunits encoded by animal mitochondria to the wider set of mitochondrial genes found in jakobids. In the present paper, we compile and describe a large database of 27,873 mitochondrial genomes currently available in GenBank, encompassing the whole eukaryotic domain. We discuss the major features of mitochondrial molecular diversity, with special reference to nucleotide composition and compositional biases; moreover, the database is made publicly available for future analyses on the MoZoo Lab GitHub page.

Complete mitochondrial genomes of two marine monogonont rotifer Brachionus manjavacas strains

The complete mitochondrial genomes of Brachionus manjavacas German strain were 10,721 bp (mitochondrial DNA I) and 12,274 bp (mitochondrial DNA II) in size, while the complete mitochondrial genomes of B. manjavacas Australian strain were 10,889 bp (mitochondrial DNA I) and 12,443 bp (mitochondrial DNA II) in size. Of 12 protein-coding genes (PCGs), 99.6% of amino acid sequences were identical between the two strains. Of 12 PCGs of both B. manjavacas strains, three genes (ND1, ATP6, and ND5) had incomplete stop codon T. Furthermore, ATA was the start codon for ND4, ND5, and CO3 genes, whereas that for other PCGs was ATG in both strains. The base compositions of 12 PCGs in the B. manjavacas strains were similar, indicating that the mitochondrial genome of the two strains was structurally conserved over evolution. The gene structure and its orientation of 12 PGCs of B. manjavacas strains were identical, as shown in other marine Brachionus rotifers and the freshwater Brachionus rotifers, while the freshwater rotifer B. calyciflorus had an additional cytochrome b gene in the mitochondrial DNA I.

Complete mitochondrial genome of Brachionus rubens from wuhu, China (Rotifera, Brachionidae)

The complete mitochondrial genome of Brachionus rubens was sequenced using primers design, clone culture, DNA extraction, LONG-PCR amplification, purification and clone sequencing. We found that it is composed of two circular chromosomes, designated mtDNA I (11,398 bp) and mtDNA II (12,820 bp). The gene content of the B. rubens mitochondrial genome was similar to that of the previously reported mitochondrial genome of B. plicatilis. It contained 22 tRNA genes, 2 rRNA genes and 12 protein-coding genes (PCGs). Four of the 12 PCGs had an incomplete stop codons, TA（cob, atp6, nd3）or T(cox3). The A + T content of B. rubens mitochondrial genome was apparently higher (mtDNA-I 70.2% and mtDNA II 70.4%) than that of the mitochondrial genome of B. plicatilis (mtDNA-I 63.9% and mtDNA-II 62.9%).

Complete mitochondrial genome of the freshwater monogonont rotifer Brachionus angularis (Rotifera, Brachionidae)

The two complete mitochondrial genomes were sequenced from the freshwater monogonont rotifer Brachionus angularis. The mitochondrial genome sequences were 10,764 bp (mitochondrial DNA I) and 12,238 bp (mitochondrial DNA II) in size, respectively. The gene structure and its orientation of 12 protein-coding genes (PCGs) of complete mitochondrial genomes of B. angularis was identical to those shown in other marine rotifers and the freshwater rotifer Brachionus rubens, but was different from the freshwater rotifer Brachionus calyciflorus. Of 12 PCGs, one gene (ND5) had incomplete stop codon. Furthermore, the start codon for CO1, ND4L, ND5, and CO2 was GTG, while the start codon for ND3 and other PCGs was ATA and ATG, respectively. The base composition of 12 PCGs in B. angularis mitogenome showed 20.4% for A, 47.3% for T, 17.5% for C, and 14.8% for G, respectively. The mitochondrial genome A + T base composition (67.7%) of 12 PCGs was higher than G + C (32.3%), while the complete mitochondrial genome A + T base composition (66.3%) was higher than G + C (33.7%).

Revisiting Ceriantharian (Anthozoa) Mitochondrial Genomes: Casting Doubts about Their Structure and Size

Recently, Stampar et al. (2019. Linear mitochondrial genome in Anthozoa (Cnidaria): a case study in. Sci Rep. 9(1):6094.) uncovered highly atypical mitochondrial genome structures in the cnidarian species Pachycerianthus magnus and Isarachnanthus nocturnus (Anthozoa, Ceriantharia). These two mitochondrial DNAs assembled as linear fragmented genomes, comprising eight and five chromosomes, respectively—architectures unlike any other anthozoan mitogenome described to date. What’s more, they have cumulative lengths of 77.8 (P. magnus) and 80.9 kb (I. nocturnus), making them the largest animal mitochondrial DNAs on record, a finding which garnered significant attention by various news media. Here, I take a closer look at the work of Stampar et al. and question their key results. I provide evidence that the currently available mitogenome sequences for I. nocturnus and P. magnus, including their structures, sizes, and chromosome numbers, should be treated with caution. More work must be done on these genomes before one can say with any certainty that they are linear, fragmented, or the largest animal mitogenomes observed to date.

The genome, transcriptome, and proteome of the fish parasite Pomphorhynchus laevis (Acanthocephala)

Thorny-headed worms (Acanthocephala) are endoparasites exploiting Mandibulata (Arthropoda) and Gnathostomata (Vertebrata). Despite their world-wide occurrence and economic relevance as a pest, genome and transcriptome assemblies have not been published before. However, such data might hold clues for a sustainable control of acanthocephalans in animal production. For this reason, we present the first draft of an acanthocephalan nuclear genome, besides the mitochondrial one, using the fish parasite Pomphorhynchus laevis (Palaeacanthocephala) as a model. Additionally, we have assembled and annotated the transcriptome of this species and the proteins encoded. A hybrid assembly of long and short reads resulted in a near-complete P. laevis draft genome of ca. 260 Mb, comprising a large repetitive portion of ca. 63%. Numbers of transcripts and translated proteins (35,683) were within the range of other members of the Rotifera-Acanthocephala clade. Our data additionally demonstrate a significant reorganization of the acanthocephalan gene repertoire. Thus, more than 20% of the usually conserved metazoan genes were lacking in P. laevis. Ontology analysis of the retained genes revealed many connections to the incorporation of carotinoids. These are probably taken up via the surface together with lipids, thus accounting for the orange coloration of P. laevis. Furthermore, we found transcripts and protein sequences to be more derived in P. laevis than in rotifers from Monogononta and Bdelloidea. This was especially the case in genes involved in energy metabolism, which might reflect the acanthocephalan ability to use the scarce oxygen in the host intestine for respiration and simultaneously carry out fermentation. Increased plasticity of the gene repertoire through the integration of foreign DNA into the nuclear genome seems to be another underpinning factor of the evolutionary success of acanthocephalans. In any case, energy-related genes and their proteins may be considered as candidate targets for the acanthocephalan control.

Complete mitochondrial genome of the euryhaline monogonont rotifer Brachionus paranguensis (Rotifera, Brachionidae)

The two complete mitochondrial genomes were sequenced from the euryhaline monogonont rotifer Brachionus paranguensis. The mitochondrial genome sequences were 11,603 bp and 12,901 bp in size, and the gene order of 12 protein-coding genes (PCGs) were identical to those of the marine rotifer Brachionus plicatilis, but the positions of some tRNAs (e.g. tRNA-Ile, tRNA-Leu[TTA], tRNA-Phe, and tRNA-Leu[CTC]) of mitochondrial DNA I were different between B. paranguensis and B. plicatilis mitochondrial genomes. Of 12 PCGs, four genes (ND1, ATPase 6, ND5, and ND3) had incomplete stop codons. Furthermore, the start codon of ND4, ND5, and CO3 genes was ATT, while the start codon of other PCGs was ATG. The base composition of 12 PCGs in B. paranguensis mitochondrial genomes was 26.6% for A, 43.0% for T, 17.65% for C, and 12.75% for G, respectively.

Complete mitochondrial genome of the freshwater monogonont rotifer Brachionus rubens (Rotifera, Brachionidae)

The two complete mitochondrial genomes were sequenced from the freshwater monogonont rotifer Brachionus rubens. The genome sequences were 12,041 bp and 13,793 bp in size, and the gene order and contents were identical to those of the freshwater rotifer B. rubens China, but were different in three tRNA-Arg, tRNA-Ile, and tRNA-Leu between both B. rubens mitochondrial genomes, while B. calyciflorus had peculiar gene order in mitochondrial DNA I. Of 12 protein-coding genes (PCGs), one gene (ND5) had incomplete stop codons. Furthermore, the start codon of ND4 and CO2 gene was ATT, while the start codon of other PCGs was ATG. The base composition of 12 PCGs in B. rubens mitogenome showed 22.5% for A, 46.5% for T, 16.3% for C, and 14.7% for G, respectively.

Complete mitochondrial genome of the freshwater monogonont rotifer Brachionus calyciflorus (Rotifera, Brachionidae)

The two complete mitochondrial genomes were sequenced from the Netherlands strain of the freshwater monogonont rotifer Brachionus calyciflorus. The mitochondrial genome sequences were 27,698 bp and 9,906 bp in size, respectively. The gene order and contents of the two B. calyciflorus strains were mostly identical to one another, except for the additional identification and translocation of several tRNAs in mitochondrial DNA I and II. Of 13 protein-coding genes (PCGs), three genes (ND1, ND5, and ND3) had incomplete stop codons. Furthermore, the start codon of ND2, CO2, and CO3 and ND4 genes was ATT, GTG, and ATA, respectively, while the start codon of other PCGs was ATG. The base composition of 13 PCGs of B. calyciflorus (the Netherlands strain) mitogenome showed 31.1% for A, 37.6% for T, 16.5% for C, and 14.8% for G, respectively.

The revised complete mitogenome sequence of the tree frog Polypedatesmegacephalus (Anura, Rhacophoridae) by next-generation sequencing and phylogenetic analysis

The mitochondrial genome (mitogenome) sequence of the tree frog Polypedates megacephalus (16,473 bp) was previously reported as having the unusual characteristic of lacking the ND5 gene. In this study, a new mitogenome of P. megacephalus (19,952 bp) was resequenced using the next-generation sequencing (NGS) and standard Sanger sequencing technologies. It was discovered that the ND5 gene was not lost but translocated to the control region (CR) from its canonical location between the ND4 and ND6 genes. In addition, a duplicated control region was found in the new mitogenome of this species. Conservative region identification of the ND5 gene and phylogenetic analysis confirmed that the ND5 gene was located between two control regions. The phylogenetic relationship among 20 related species of anura revealed a rearrangement of the ND5 gene during the evolutionary process. These results also highlighted the advantages of next-generation sequencing. It will not only decrease the time and cost of sequencing, but also will eliminate the errors in published mitogenome databases.

The bipartite mitochondrial genome of Ruizia karukerae (Rhigonematomorpha, Nematoda)

Mitochondrial genes and whole mitochondrial genome sequences are widely used as molecular markers in studying population genetics and resolving both deep and shallow nodes in phylogenetics. In animals the mitochondrial genome is generally composed of a single chromosome, but mystifying exceptions sometimes occur. We determined the complete mitochondrial genome of the millipede-parasitic nematode Ruizia karukerae and found its mitochondrial genome consists of two circular chromosomes, which is highly unusual in bilateral animals. Chromosome I is 7,659 bp and includes six protein-coding genes, two rRNA genes and nine tRNA genes. Chromosome II comprises 7,647 bp, with seven protein-coding genes and 16 tRNA genes. Interestingly, both chromosomes share a 1,010 bp sequence containing duplicate copies of cox2 and three tRNA genes (trnD, trnG and trnH), and the nucleotide sequences between the duplicated homologous gene copies are nearly identical, suggesting a possible recent genesis for this bipartite mitochondrial genome. Given that little is known about the formation, maintenance or evolution of abnormal mitochondrial genome structures, R. karukerae mtDNA may provide an important early glimpse into this process.

The complete mitochondrial genome of parasitic nematode Camallanus cotti: extreme discontinuity in the rate of mitogenomic architecture evolution within the Chromadorea class

BACKGROUND

Complete mitochondrial genomes are much better suited for the taxonomic identification and phylogenetic studies of nematodes than morphology or traditionally-used molecular markers, but they remain unavailable for the entire Camallanidae family (Chromadorea). As the only published mitogenome in the Camallanina suborder (Dracunculoidea superfamily) exhibited a unique gene order, the other objective of this research was to study the evolution of mitochondrial architecture in the Spirurida order. Thus, we sequenced the complete mitogenome of the Camallanus cotti fish parasite and conducted structural and phylogenomic comparative analyses with all available Spirurida mitogenomes.

RESULTS

The mitogenome is exceptionally large (17,901 bp) among the Chromadorea and, with 46 (pseudo-) genes, exhibits a unique architecture among nematodes. Six protein-coding genes (PCGs) and six tRNAs are duplicated. An additional (seventh) tRNA (Trp) was probably duplicated by the remolding of tRNA-Ser2 (missing). Two pairs of these duplicated PCGs might be functional; three were incomplete and one contained stop codons. Apart from Ala and Asp, all other duplicated tRNAs are conserved and probably functional. Only 19 unique tRNAs were found. Phylogenomic analysis included Gnathostomatidae (Spirurina) in the Camallanina suborder.

CONCLUSIONS

Within the Nematoda, comparable PCG duplications were observed only in the enoplean Mermithidae family, but those result from mitochondrial recombination, whereas characteristics of the studied mitogenome suggest that likely rearrangement mechanisms are either a series of duplications, transpositions and random loss events, or duplication, fragmentation and subsequent reassembly of the mitogenome. We put forward a hypothesis that the evolution of mitogenomic architecture is extremely discontinuous, and that once a long period of stasis in gene order and content has been punctuated by a rearrangement event, such a destabilised mitogenome is much more likely to undergo subsequent rearrangement events, resulting in an exponentially accelerated evolutionary rate of mitogenomic rearrangements. Implications of this model are particularly important for the application of gene order similarity as an additive source of phylogenetic information. Chromadorean nematodes, and particularly Camallanina clade (with C. cotti as an example of extremely accelerated rate of rearrangements), might be a good model to further study this discontinuity in the dynamics of mitogenomic evolution.

The Multipartite Mitochondrial Genome of Enteromyxum leei (Myxozoa): Eight Fast-Evolving Megacircles

Myxozoans are a large group of poorly characterized cnidarian parasites. To gain further insight into their evolution, we sequenced the mitochondrial (mt) genome of Enteromyxum leei and reevaluate the mt genome structure of Kudoa iwatai. Although the typical animal mt genome is a compact, 13-25 kb, circular chromosome, the mt genome of E. leei was found to be fragmented into eight circular chromosomes of ∼23 kb, making it the largest described animal mt genome. Each chromosome was found to harbor a large noncoding region (∼15 kb), nearly identical between chromosomes. The protein coding genes show an unusually high rate of sequence evolution and possess little similarity to their cnidarian homologs. Only five protein coding genes could be identified and no tRNA genes. Surprisingly, the mt genome of K. iwatai was also found to be composed of two chromosomes. These observations confirm the remarkable plasticity of myxozoan mt genomes.

Untangling Heteroplasmy, Structure, and Evolution of an Atypical Mitochondrial Genome by PacBio Sequencing

The highly compact mitochondrial (mt) genome of terrestrial isopods (Oniscidae) presents two unusual features. First, several loci can individually encode two tRNAs, thanks to single nucleotide polymorphisms at anticodon sites. Within-individual variation (heteroplasmy) at these loci is thought to have been maintained for millions of years because individuals that do not carry all tRNA genes die, resulting in strong balancing selection. Second, the oniscid mtDNA genome comes in two conformations: a ∼14 kb linear monomer and a ∼28 kb circular dimer comprising two monomer units fused in palindrome. We hypothesized that heteroplasmy actually results from two genome units of the same dimeric molecule carrying different tRNA genes at mirrored loci. This hypothesis, however, contradicts the earlier proposition that dimeric molecules result from the replication of linear monomers-a process that should yield totally identical genome units within a dimer. To solve this contradiction, we used the SMRT (PacBio) technology to sequence mirrored tRNA loci in single dimeric molecules. We show that dimers do present different tRNA genes at mirrored loci; thus covalent linkage, rather than balancing selection, maintains vital variation at anticodons. We also leveraged unique features of the SMRT technology to detect linear monomers closed by hairpins and carrying noncomplementary bases at anticodons. These molecules contain the necessary information to encode two tRNAs at the same locus, and suggest new mechanisms of transition between linear and circular mtDNA. Overall, our analyses clarify the evolution of an atypical mt genome where dimerization counterintuitively enabled further mtDNA compaction.

Complete mitochondrial genome of the monogonont rotifer Brachionus rotundiformis (Rotifera, Brachionidae)

The two complete mitochondrial genomes were sequenced from the monogonont rotifer Brachionus rotundiformis. The genome sequences were 10,268 bp and 11,703 bp in size, and the gene order and contents were identical with those of B. koreanus but were different in tRNA-Cys with B. plicatilis mitochondrial genomes. Of the 12 protein-coding genes (PCGs), five genes (ND1, ATP6, ND5, CO3, ND3) had incomplete stop codons. Furthermore, the start codon of ND4 and CO3 gene was ATA, while the start codon of other PCGs was ATG. The base composition of B. rotundiformis mitogenome shows an anti-G bias (12.05% and 10.24%) on the second and third position of the PCGs, respectively.

Monogonont Rotifer, Brachionus calyciflorus, Possesses Exceptionally Large, Fragmented Mitogenome

In contrast to the highly conserved mitogenomic structure and organisation in most animals (including rotifers), the two previously sequenced monogonont rotifer mitogenomes were fragmented into two chromosomes similar in size, each of which possessed one major non-coding region (mNCR) of about 4–5 Kbp. To further explore this phenomenon, we have sequenced and analysed the mitogenome of one of the most studied monogonont rotifers, Brachionus calyciflorus. It is also composed of two circular chromosomes, but the chromosome-I is extremely large (27 535 bp; 3 mNCRs), whereas the chromosome-II is relatively small (9 833 bp; 1 mNCR). With the total size of 37 368 bp, it is one of the largest metazoan mitogenomes ever reported. In comparison to other monogononts, gene distribution between the two chromosomes and gene order are different and the number of mNCRs is doubled. Atp8 was not found (common in rotifers), and Cytb was present in two copies (the first report in rotifers). A high number (99) of SNPs indicates fast evolution of the Cytb-1 copy. The four mNCRs (5.3–5.5 Kb) were relatively similar. Publication of this sequence shall contribute to the understanding of the evolutionary history of the unique mitogenomic organisation in this group of rotifers.

The Complete Mitochondrial Genome of Aleurocanthus camelliae: Insights into Gene Arrangement and Genome Organization within the Family Aleyrodidae

There are numerous gene rearrangements and transfer RNA gene absences existing in mitochondrial (mt) genomes of Aleyrodidae species. To understand how mt genomes evolved in the family Aleyrodidae, we have sequenced the complete mt genome of Aleurocanthus camelliae and comparatively analyzed all reported whitefly mt genomes. The mt genome of A. camelliae is 15,188 bp long, and consists of 13 protein-coding genes, two rRNA genes, 21 tRNA genes and a putative control region (GenBank: KU761949). The tRNA gene, trnI, has not been observed in this genome. The mt genome has a unique gene order and shares most gene boundaries with Tetraleurodes acaciae. Nineteen of 21 tRNA genes have the conventional cloverleaf shaped secondary structure and two (trnS₁ and trnS₂) lack the dihydrouridine (DHU) arm. Using ARWEN and homologous sequence alignment, we have identified five tRNA genes and revised the annotation for three whitefly mt genomes. This result suggests that most absent genes exist in the genomes and have not been identified, due to be lack of technology and inference sequence. The phylogenetic relationships among 11 whiteflies and Drosophila melanogaster were inferred by maximum likelihood and Bayesian inference methods. Aleurocanthus camelliae and T. acaciae form a sister group, and all three Bemisia tabaci and two Bemisia afer strains gather together. These results are identical to the relationships inferred from gene order. We inferred that gene rearrangement plays an important role in the mt genome evolved from whiteflies.

Animal Mitochondrial DNA as We Do Not Know It: mt-Genome Organization and Evolution in Nonbilaterian Lineages

Animal mitochondrial DNA (mtDNA) is commonly described as a small, circular molecule that is conserved in size, gene content, and organization. Data collected in the last decade have challenged this view by revealing considerable diversity in animal mitochondrial genome organization. Much of this diversity has been found in nonbilaterian animals (phyla Cnidaria, Ctenophora, Placozoa, and Porifera), which, from a phylogenetic perspective, form the main branches of the animal tree along with Bilateria. Within these groups, mt-genomes are characterized by varying numbers of both linear and circular chromosomes, extra genes (e.g. atp9, polB, tatC), large variation in the number of encoded mitochondrial transfer RNAs (tRNAs) (0-25), at least seven different genetic codes, presence/absence of introns, tRNA and mRNA editing, fragmented ribosomal RNA genes, translational frameshifting, highly variable substitution rates, and a large range of genome sizes. This newly discovered diversity allows a better understanding of the evolutionary plasticity and conservation of animal mtDNA and provides insights into the molecular and evolutionary mechanisms shaping mitochondrial genomes.

The mitochondrial genome of Globodera ellingtonae is composed of two circles with segregated gene content and differential copy numbers

Background

The evolution of animal mitochondrial (mt) genomes has resulted in a highly conserved structure: a single compact circular chromosome approximately 14 to 20 kb long. Within the last two decades exceptions to this conserved structure, such as the division of the genome into multiple chromosomes, have been reported in a diverse set of metazoans. We report on the two circle multipartite mt genome of a newly described cyst nematode, Globodera ellingtonae.

Results

The G. ellingtonae mt genome was found to be comprised of two circles, each larger than any other multipartite circular mt chromosome yet reported, and both were larger than the single mt circle of the model nematode Caenorhabditis elegans. The genetic content of the genome was disproportionately divided between the two circles, although they shared a ~6.5 kb non-coding region. The 17.8 kb circle (mtDNA-I) contained ten protein-coding genes and two tRNA genes, whereas the 14.4 kb circle (mtDNA-II) contained two protein-coding genes, 20 tRNA genes and both rRNA genes. Perhaps correlated with this division of genetic content, the copy number of mtDNA-II was more than four-fold that of mtDNA-I in individual nematodes. The difference in copy number increased between second-stage and fourth-stage juveniles.

Conclusions

The segregation of gene types to different mt circles in G. ellingtonae could provide benefit by localizing gene functional types to independent transcriptional units. This is the first report of both two-circle and several-circle mt genomes within a single genus. The differential copy number associated with this multipartite mt organization could provide a model system for deconstructing mechanisms regulating mtDNA copy number both in somatic cells and during germline development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3047-x) contains supplementary material, which is available to authorized users.

The massive mitochondrial genome of the angiosperm Silene noctiflora is evolving by gain or loss of entire chromosomes

Across eukaryotes, mitochondria exhibit staggering diversity in genomic architecture, including the repeated evolution of multichromosomal structures. Unlike in the nucleus, where mitosis and meiosis ensure faithful transmission of chromosomes, the mechanisms of inheritance in fragmented mitochondrial genomes remain mysterious. Multichromosomal mitochondrial genomes have recently been found in multiple species of flowering plants, including Silene noctiflora, which harbors an unusually large and complex mitochondrial genome with more than 50 circular-mapping chromosomes totaling ∼7 Mb in size. To determine the extent to which such genomes are stably maintained, we analyzed intraspecific variation in the mitochondrial genome of S. noctiflora. Complete genomes from two populations revealed a high degree of similarity in the sequence, structure, and relative abundance of mitochondrial chromosomes. For example, there are no inversions between the genomes, and there are only nine SNPs in 25 kb of protein-coding sequence. Remarkably, however, these genomes differ in the presence or absence of 19 entire chromosomes, all of which lack any identifiable genes or contain only duplicate gene copies. Thus, these mitochondrial genomes retain a full gene complement but carry a highly variable set of chromosomes that are filled with presumably dispensable sequence. In S. noctiflora, conventional mechanisms of mitochondrial sequence divergence are being outstripped by an apparently nonadaptive process of whole-chromosome gain/loss, highlighting the inherent challenge in maintaining a fragmented genome. We discuss the implications of these findings in relation to the question of why mitochondria, more so than plastids and bacterial endosymbionts, are prone to the repeated evolution of multichromosomal genomes.

Multiple Conserved Heteroplasmic Sites in tRNA Genes in the Mitochondrial Genomes of Terrestrial Isopods (Oniscidea)

Mitochondrial genome structure and organization are relatively conserved among metazoans. However, in many isopods, especially the terrestrial isopods (Oniscidea), the mitochondrial genome consists of both ∼14-kb linear monomers and ∼28-kb circular dimers. This unusual organization is associated with an ancient and conserved constitutive heteroplasmic site. This heteroplasmy affects the anticodon of a tRNA gene, allowing this single locus to function as a “dual” tRNA gene for two different amino acids. Here, we further explore the evolution of these unusual mitochondrial genomes by assembling complete mitochondrial sequences for two additional Oniscidean species, Trachelipus rathkei and Cylisticus convexus. Strikingly, we find evidence of two additional heteroplasmic sites that also alter tRNA anticodons, creating additional dual tRNA genes, and that are conserved across both species. These results suggest that the unique linear/circular organization of isopods’ mitochondrial genomes may facilitate the evolution of stable mitochondrial heteroplasmies, and, conversely, once such heteroplasmies have evolved, they constrain the multimeric structure of the mitochondrial genome in these species. Finally, we outline some possible future research directions to identify the factors influencing mitochondrial genome evolution in this group.

The mitochondrial genome of the glomeromycete Rhizophagus sp. DAOM 213198 reveals an unusual organization consisting of two circular chromosomes

Mitochondrial (mt) genomes are intensively studied in Ascomycota and Basidiomycota, but they are poorly documented in basal fungal lineages. In this study, we sequenced the complete mtDNA of Rhizophagus sp. DAOM 213198, a close relative to Rhizophagus irregularis, a widespread, ecologically and economical relevant species belonging to Glomeromycota. Unlike all other known taxonomically close relatives harboring a full-length circular chromosome, mtDNA of Rhizophagus sp. reveals an unusual organization with two circular chromosomes of 61,964 and 29,078 bp. The large chromosome contained nine protein-coding genes (atp9, nad5, cob, nad4, nad1, nad4L, cox1, cox2, and atp8), small subunit rRNA gene (rns), and harbored 20 tRNA-coding genes and 10 orfs, while the small chromosome contained five protein-coding genes (atp6, nad2, nad3, nad6, and cox3), large subunit rRNA gene (rnl) in addition to 5 tRNA-coding genes, and 8 plasmid-related DNA polymerases (dpo). Although structural variation of plant mt genomes is well documented, this study is the first report of the presence of two circular mt genomes in arbuscular mycorrhizal fungi. Interestingly, the presence of dpo at the breakage point in intergenes cox1-cox2 and rnl-atp6 for large and small mtDNAs, respectively, could be responsible for the conversion of Rhizophagus sp. mtDNA into two chromosomes. Using quantitative real-time polymerase chain reaction, we found that both mtDNAs have an equal abundance. This study reports a novel mtDNA organization in Glomeromycota and highlights the importance of studying early divergent fungal lineages to describe novel evolutionary pathways in the fungal kingdom.

Evolution of multipartite mitochondrial genomes in the booklice of the genus Liposcelis (Psocoptera)

BACKGROUND

The genus Liposcelis (Psocoptera: Troctomorpha) has more than 120 species with a worldwide distribution and they pose a risk for global food security. The organization of mitochondrial (mt) genomes varies between the two species of booklice investigated in the genus Liposcelis. Liposcelis decolor has its mt genes on a single chromosome, like most other insects; L. bostrychophila, however, has a multipartite mt genome with genes on two chromosomes.

RESULTS

To understand how multipartite mt genome organization evolved in the genus Liposcelis, we sequenced the mt genomes of L. entomophila and L. paeta in this study. We found that these two species of booklice also have multipartite mt genomes, like L. bostrychophila, with the mt genes we identified on two chromosomes. Numerous pseudo mt genes and non-coding regions were found in the mt genomes of these two booklice, and account for 30% and 10% respectively of the entire length we sequenced. In L. bostrychophila, the mt genes are distributed approximately equally between the two chromosomes. In L. entomophila and L. paeta, however, one mt chromosome has most of the genes we identified whereas the other chromosome has largely pseudogenes and non-coding regions. L. entomophila and L. paeta differ substantially from each other and from L. bostrychophila in gene content and gene arrangement in their mt chromosomes.

CONCLUSIONS

Our results indicate unusually fast evolution in mt genome organization in the booklice of the genus Liposcelis, and reveal different patterns of mt genome fragmentation among L. bostrychophila, L. entomophila and L. paeta.

The mitochondrial genome of Frankliniella intonsa: insights into the evolution of mitochondrial genomes at lower taxonomic levels in Thysanoptera

Thrips is an ideal group for studying the evolution of mitochondrial (mt) genomes in the genus and family due to independent rearrangements within this order. The complete sequence of the mitochondrial DNA (mtDNA) of the flower thrips Frankliniella intonsa has been completed and annotated in this study. The circular genome is 15,215bp in length with an A+T content of 75.9% and contains the typical 37 genes and it has triplicate putative control regions. Nucleotide composition is A+T biased, and the majority of the protein-coding genes present opposite CG skew which is reflected by the nucleotide composition, codon and amino acid usage. Although the known thrips have massive gene rearrangements, it showed no reversal of strand asymmetry. Gene rearrangements have been found in the lower taxonomic levels of thrips. Three tRNA genes were translocated in the genus Frankliniella and eight tRNA genes in the family Thripidae. Although the gene arrangements of mt genomes of all three thrips species differ massively from the ancestral insect, they are all very similar to each other, indicating that there was a large rearrangement somewhere before the most recent common ancestor of these three species and very little genomic evolution or rearrangements after then. The extremely similar sequences among the CRs suggest that they are ongoing concerted evolution. Analyses of the up and downstream sequence of CRs reveal that the CR2 is actually the ancestral CR. The three CRs are in the same spot in each of the three thrips mt genomes which have the identical inverted genes. These characteristics might be obtained from the most recent common ancestor of this three thrips. Above observations suggest that the mt genomes of the three thrips keep a single massive rearrangement from the common ancestor and have low evolutionary rates among them.

The complete mitochondrial genome of the booklouse, Liposcelis decolor: insights into gene arrangement and genome organization within the genus Liposcelis

Booklice in the genus Liposcelis are pests of stored grain products. They pose a considerable economic threat to global food security and safety. To date, the complete mitochondrial genome has only been determined for a single booklouse species Liposcelis bostrychophila. Unlike most bilateral animals, which have their 37 mt genes on one circular chromosome, ≈15 kb in size, the mt genome of L. bostrychophila has two circular chromosomes, 8 and 8.5 kb in size. Here, we report the mt genome of another booklouse, Liposcelis decolor. The mt genome of L. decolor has the typical mt chromosome of bilateral animals, 14,405 bp long with 37 genes (13 PCGs, 22 tRNAs and 2 rRNAs). However, the arrangement of these genes in L. decolor differs substantially from that observed in L. bostrychophila and other insects. With the exception of atp8-atp6, L. decolor differs from L. bostrychophila in the arrangement of all of the other 35 genes. The variation in the mt genome organization and mt gene arrangement between the two Liposcelis species is unprecedented for closely related animals in the same genus. Furthermore, our results indicate that the two-chromosome mt genome organization observed in L. bostrychophila likely evolved recently after L. bostrychophila and L. decolor split from their most recent common ancestor.

Fragmented mitochondrial genomes of the rat lice, Polyplax asiatica and Polyplax spinulosa: intra-genus variation in fragmentation pattern and a possible link between the extent of fragmentation and the length of life cycle

BACKGROUND

Blood-sucking lice (suborder Anoplura) parasitize eutherian mammals with 67% of the 540 described species found on rodents. The five species of blood-sucking lice that infest humans and pigs have fragmented mitochondrial genomes and differ substantially in the extent of fragmentation. To understand whether, or not, any life-history factors are linked to such variation, we sequenced the mt genomes of Polyplax asiatica and Polyplax spinulosa, collected from the greater bandicoot rat, Bandicota indica, and the Asian house rat, Rattus tanezumi, respectively.

RESULTS

We identified all of the 37 mitochondrial genes common to animals in Polyplax asiatica and Polyplax spinulosa. The mitochondrial genes of these two rat lice are on 11 circular minichromosomes; each minichromosome is 2-4 kb long and has 2-7 genes. The two rat lice share the same pattern for the distribution of the protein-coding genes and ribosomal RNA genes over the minichromosomes, but differ in the pattern for the distribution of 8 of the 22 transfer RNA genes. The mitochondrial genomes of the Polyplax rat lice have 3.4 genes, on average, on each minichromosome and, thus, are less fragmented than those of the human lice (2.1 and 2.4 genes per minichromosome), but are more fragmented than those of the pig lice (4.1 genes per minichromosome).

CONCLUSIONS

Our results revealed distinct patterns of mitochondrial genome fragmentation within the genus Polyplax and, furthermore, indicated a possible inverse link between the extent of mitochondrial genome fragmentation and the length of life cycle of the blood-sucking lice.

Insect mitochondrial genomics: implications for evolution and phylogeny

The mitochondrial (mt) genome is, to date, the most extensively studied genomic system in insects, outnumbering nuclear genomes tenfold and representing all orders versus very few. Phylogenomic analysis methods have been tested extensively, identifying compositional bias and rate variation, both within and between lineages, as the principal issues confronting accurate analyses. Major studies at both inter- and intraordinal levels have contributed to our understanding of phylogenetic relationships within many groups. Genome rearrangements are an additional data type for defining relationships, with rearrangement synapomorphies identified across multiple orders and at many different taxonomic levels. Hymenoptera and Psocodea have greatly elevated rates of rearrangement offering both opportunities and pitfalls for identifying rearrangement synapomorphies in each group. Finally, insects are model systems for studying aberrant mt genomes, including truncated tRNAs and multichromosomal genomes. Greater integration of nuclear and mt genomic studies is necessary to further our understanding of insect genomic evolution.

A new intertidal Brachionus and intrageneric phylogenetic relationships among Brachionus as revealed by allometry and CO1-ITS1 gene analysis

Background

The rotifer genus Brachionus is distributed worldwide along a gradient from freshwater to seawater habitats. This genus is recognized as a suitable organism for testing environmental and evolutionary hypotheses. Here we present the description of a new species and the phylogenetic relationship within the genus Brachionus as revealed by morphometric allometry of two representatives and CO1-ITS1 gene analysis of five representatives belonging to Brachionus distributed in the fresh and seawaters of Korea, Japan, and China.

Results

Similarities of populations were studied using nuclear rDNA ITS1 sequences from Brachionus spp. collected from different geographical areas of Far East Asia. The phylogeographic analysis of nuclear DNA ITS1 and mitochondrial CO1 sequences showed that Brachionus from South Korea formed five distinct clades according to their geographic origin. Interspecific differences suggest that Brachionus species established on conventional morphological characters also forms five separate clades. When ITS1-CO1 of the Asian specimens was compared with representatives of Brachionus worldwide, a high genotypic similarity was found. However, they slightly differed between localities.

Conclusions

This study sets a first step for an integrative morphological and molecular characterization of the diversity contained within the ecologically and economically important rotifer genus Brachionus.

A commentary on the XIII(th) International Rotifer Symposium (Shillong, 2012)

Rotifers have attracted the attention of biologists for well over 200 years. Interest in these exquisite animals rests in their diverse morphology, short generation time resulting in high growth rates, ability to withstand desiccation, and wide distribution, coupled with evidence of cryptic speciation. Moreover, three modes of reproduction are present in the phylum: obligatory sexuality, cyclical parthenogenesis, and obligatory ameiotic parthenogenesis. Thus, this phylum offers a rich field of study. Recognizing the need to share advances in knowledge, a triennial meeting, the International Rotifer Symposium (IRS), was begun in 1976. The most recent symposium (13(th) IRS) was held at Shillong (India) from 18-24, November 2012. In this commentary we considered the development of rotifer research as viewed through the lens of more than 35 years of IRS. Initially papers presented at the IRS focused on ecology, morphology, and pure taxonomic problems, with little applied work being reported. However, after more than three decades, the emphasis has swung to a balance of both basic (e.g., aging, ecology, genetics, and taxonomy) and applied (aquaculture and ecotoxicology) research.

The complete mitochondrial genome of Solemya velum (Mollusca: Bivalvia) and its relationships with conchifera

Background

Bivalve mitochondrial genomes exhibit a wide array of uncommon features, like extensive gene rearrangements, large sizes, and unusual ways of inheritance. Species pertaining to the order Solemyida (subclass Opponobranchia) show many peculiar evolutionary adaptations, f.i. extensive symbiosis with chemoautotrophic bacteria. Despite Opponobranchia are central in bivalve phylogeny, being considered the sister group of all Autobranchia, a complete mitochondrial genome has not been sequenced yet.

Results

In this paper, we characterized the complete mitochondrial genome of the Atlantic awning clam Solemya velum: A-T content, gene arrangement and other features are more similar to putative ancestral mollusks than to other bivalves. Two supranumerary open reading frames are present in a large, otherwise unassigned, region, while the origin of replication could be located in a region upstream to the cox3 gene.

Conclusions

We show that S. velum mitogenome retains most of the ancestral conchiferan features, which is unusual among bivalve mollusks, and we discuss main peculiarities of this first example of an organellar genome coming from the subclass Opponobranchia. Mitochondrial genomes of Solemya (for bivalves) and Haliotis (for gastropods) seem to retain the original condition of mollusks, as most probably exemplified by Katharina.

Evolution of extensively fragmented mitochondrial genomes in the lice of humans

Bilateral animals are featured by an extremely compact mitochondrial (mt) genome with 37 genes on a single circular chromosome. The human body louse, Pediculus humanus, however, has its mt genes on 20 minichromosomes. We sequenced the mt genomes of two other human lice: the head louse, P. capitis, and the pubic louse, Pthirus pubis. Comparison among the three human lice revealed the presence of fragmented mt genomes in their most recent common ancestor, which lived ∼7 Ma. The head louse has exactly the same set of mt minichromosomes as the body louse, indicating that the number of minichromosomes, and the gene content and gene arrangement in each minichromosome have remained unchanged since the body louse evolved from the head louse ∼107,000 years ago. The pubic louse has the same pattern of one protein-coding or rRNA gene per minichromosome (except one minichromosome with two protein-coding genes, atp6 and atp8) as the head louse and the body louse. This pattern is apparently ancestral to all human lice and has been stable for at least 7 Myr. Most tRNA genes of the pubic louse, however, are on different minichromosomes when compared with their counterparts in the head louse and the body louse. It is evident that rearrangement of four tRNA genes (for leucine, arginine and glycine) was due to gene-identity switch by point mutation at the third anticodon position or by homologous recombination, whereas rearrangement of other tRNA genes was by gene translocation between minichromosomes, likely caused by minichromosome split via gene degeneration and deletion.

Genetic aspects of mitochondrial genome evolution

Many years of extensive studies of metazoan mitochondrial genomes have established differences in gene arrangements and genetic codes as valuable phylogenetic markers. Understanding the underlying mechanisms of replication, transcription and the role of the control regions which cause e.g. different gene orders is important to assess the phylogenetic signal of such events. This review summarises and discusses, for the Metazoa, the general aspects of mitochondrial transcription and replication with respect to control regions as well as several proposed models of gene rearrangements. As whole genome sequencing projects accumulate, more and more observations about mitochondrial gene transfer to the nucleus are reported. Thus occurrence and phylogenetic aspects concerning nuclear mitochondrial-like sequences (NUMTS) is another aspect of this review.