Sequence and Characterization of Six Mitochondrial Subgenomes from Globodera rostochiensis: Multipartite Structure Is Conserved Among Close Nematode Relatives

High-Throughput Sequencing Enables Rapid Analyses of Nematode Mitochondrial Genomes from an Environmental Sample

Mitochondrial genomes serve as essential tools in evolutionary biology, phylogenetics, and population genetics due to their maternal inheritance, lack of recombination, and conserved structure. Traditional morphological methods for identifying nematodes are often insufficient for distinguishing cryptic species complexes. This study highlights recent advancements in nematode mitochondrial genome research, particularly the impact of long-read sequencing technologies such as Oxford Nanopore. These technologies have facilitated the assembly of mitochondrial genomes from mixed soil samples, overcoming challenges associated with designing specific primers for long PCR amplification across different groups of parasitic nematodes. In this study, we successfully recovered and assembled eleven nematode mitochondrial genomes using long-read sequencing, including those of two plant-parasitic nematode species. Notably, we detected Heterodera cruciferae in Victoria, expanding its known geographic range within Australia. Additionally, short-read sequencing data from a previous draft genome study revealed the presence of the mitochondrial genome of Heterodera filipjevi. Comparative analyses of Heterodera mitogenomes revealed conserved protein-coding genes essential for oxidative phosphorylation, as well as gene rearrangements and variations in transfer RNA placement, which may reflect adaptations to parasitic lifestyles. The consistently high A+T content and strand asymmetry observed across species align with trends reported in related genera. This study demonstrates the utility of long-read sequencing for identifying coexisting nematode species in agricultural fields, providing a rapid, accurate, and comprehensive alternative to traditional diagnostic methods. By incorporating non-target endemic species into public databases, this approach enhances biodiversity records and informs biosecurity strategies. These findings reinforce the potential of mitochondrial genomics to strengthen Australia's as well as the global biosecurity framework against plant-parasitic nematode threats.

One Ring does not rule them all: Linear mtDNA in Metazoa

Recent advances in genome sequencing technologies have facilitated the exploration of the architecture of genomes, including mitochondrial genomes (mtDNA). In particular, whole genome sequencing has provided easier access to mitochondrial genomes with unusual organizations, which were difficult to obtain using traditional PCR-based approaches. As a consequence, there has been a steep increase in complete mtDNA sequences, particularly for Metazoa. The popular view of metazoan mtDNA is that of a small gene-dense circular chromosome. This view clashes with discoveries of a number of linear mtDNAs, particularly in non-bilaterian animals. Here, we review the distribution of linear mtDNA in Metazoa, namely in isopods, cnidarians, and sponges. We discuss the multiple origins of linear mitogenomes in these clades, where linearity has been linked to the likely insertion of a linear plasmid in cnidarians and the demosponge Acanthella acuta, while fixation of a heteroplasmy in the anticodon site of a tRNA might be responsible for the monolinear form of the mtDNA in some isopods. We also summarize our current knowledge of mechanisms that maintain the integrity of linear mitochromosomes, where a recurrent theme is the presence of terminal repeats that likely play the role of telomeres. We caution in defining a linear chromosome as complete, particularly when coding sequences and key features of linear DNA are missing. Finally, we encourage authors interested in mitogenome science to utilize all available data for linear mtDNA, including those tagged as "incomplete" or "unverified" in public databases, as they can still provide useful information such as phylogenetic characters and gene order.

Insect Mitochondrial Genomics: A Decade of Progress

The past decade has seen the availability of insect genomic data explode, with mitochondrial (mt) genome data seeing the greatest growth. The widespread adoption of next-generation sequencing has solved many earlier methodological limitations, allowing the routine sequencing of whole mt genomes, including from degraded or museum specimens and in parallel to nuclear genomic projects. The diversity of available taxa now allows finer-scale comparisons between mt and nuclear phylogenomic analyses; high levels of congruence have been found for most orders, with some significant exceptions (e.g., Odonata, Mantodea, Diptera). The evolution of mt gene rearrangements and their association with haplodiploidy have been tested with expanded taxonomic sampling, and earlier proposed trends have been largely supported. Multiple model systems have been developed based on findings unique to insects, including mt genome fragmentation (lice and relatives) and control region duplication (thrips), allowing testing of hypothesized evolutionary drivers of these aberrant genomic phenomena. Finally, emerging research topics consider the contributions of mt genomes to insect speciation and habitat adaption, with very broad potential impacts. Integration between insect mt genomic research and other fields within entomology continues to be our field's greatest opportunity and challenge.

Mitochondrial genome fragmentation is correlated with increased rates of molecular evolution

While mitochondrial genome content and organization is quite diverse across all Eukaryotes, most bilaterian animal mitochondrial genomes (mitogenomes) exhibit highly conserved gene content and organisation, with genes typically encoded on a single circular chromosome. However, many species of parasitic lice (Insecta: Phthiraptera) are among the notable exceptions, having mitogenomes fragmented into multiple circular chromosomes. To better understand the process of mitogenome fragmentation, we conducted a large-scale genomic study of a major group of lice, Amblycera, with extensive taxon sampling. Analyses of the evolution of mitogenome structure across a phylogenomic tree of 90 samples from 53 genera revealed evidence for multiple independent origins of mitogenome fragmentation, some inferred to have occurred less than five million years ago. We leveraged these many independent origins of fragmentation to compare the rates of DNA substitution and gene rearrangement, specifically contrasting branches with fragmented and non-fragmented mitogenomes. We found that lineages with fragmented mitochondrial genomes had significantly higher rates of mitochondrial sequence evolution. In addition, lineages with fragmented mitochondrial genomes were more likely to have mitogenome gene rearrangements than those with single-chromosome mitochondrial genomes. By combining phylogenomics and mitochondrial genomics we provide a detailed portrait of mitogenome evolution across this group of insects with a remarkably unstable mitogenome structure, identifying processes of molecular evolution that are correlated with mitogenome fragmentation.

Comparative mitogenomics of Brachiopods reveals conservatism in articulate species and unusualness in inarticulate species

BACKGROUND

Brachiopods are a phylum of marine invertebrates with over 10,000 fossil species. Today, there are fewer than 500 extant species assigned to the class Articulata or Inarticulata and for which knowledge of evolutionary genetics and genomics is still poor. Until now, complete mitogenome sequences of two inarticulate species and four articulate species were available.

METHODS AND RESULTS

The complete mitogenome of the inarticulate brachiopod species Lingula reevii (20,778 bp) was obtained by using next generation sequencing. It contains 12 protein-coding genes (the annotation of atp8 is unsure), two ribosomal RNA genes, 26 transfer RNA genes, and one supernumerary ORF that is also conserved in the inarticulate species Lingula anatina. It is hypothesized that this ORF could represent a Lingula-specific mtORFan gene (without obvious homology to other genes). Comparative mitogenomics indicate the mitochondrial gene order of L. reevii is unique among brachiopods, and that compared to articulate species, inarticulate species exhibit massive mitogenome rearrangements, deviant ATP8 protein sequences and supernumerary ORFs, possibly representing species- or lineage-specific mtORFan genes.

CONCLUSION

The results of this study enrich genetics knowledge of extant brachiopods, which may eventually help to test hypotheses about their decline.

Independent evolution of highly variable, fragmented mitogenomes of parasitic lice

The mitochondrial genomes (mitogenomes) of bilaterian animals are highly conserved structures that usually consist of a single circular chromosome. However, several species of parasitic lice (Insecta: Phthiraptera) possess fragmented mitogenomes, where the mitochondrial genes are present on separate, circular chromosomes. Nevertheless, the extent, causes, and consequences of this structural variation remain poorly understood. Here, we combined new and existing data to better understand the evolution of mitogenome fragmentation in major groups of parasitic lice. We found strong evidence that fragmented mitogenomes evolved many times within parasitic lice and that the level of fragmentation is highly variable, including examples of heteroplasmic arrangements. We also found a significant association between mitochondrial fragmentation and signatures of relaxed selection. Mitochondrial fragmentation was also associated with changes to a lower AT%, possibly due to differences in mutation biases. Together, our results provide a significant advance in understanding the process of mitogenome fragmentation and provide an important perspective on mitochondrial evolution in eukaryotes.

An Integrated Approach for Synonymization of Rotylenchus rhomboides with R. goodeyi (Nematoda: Hoplolaimidae) Reveals High Intraspecific Mitogenomic Variation

Rotylenchus is a widely distributed, economically important plant-parasitic nematode group whose species-level identification relies largely on limited morphological characters, including character-based tabular keys and molecular data of ribosomal and mitochondrial genes. In this study, a combined morphological and molecular analysis of three populations of Rotylenchus goodeyi from Belgium, Poland, and the Netherlands revealed important character variations of this species, leading to synonymization of R. rhomboides with R. goodeyi and a high nucleotide variation within cox1 gene sequences in these populations. Additional Illumina sequencing of DNA from individuals of the Dutch population revealed two variants of mitogenomes, each approximately 23 Kb in size, differing by approximately 9% and containing 11 protein-coding genes, 2 ribosomal RNA genes, and as many as 29 transfer RNA genes. In addition to the first representative whole-genome shotgun sequence datasets of the genus Rotylenchus, this study also provides the full-length mitogenome and the ribosomal DNA sequences of R. goodeyi.

Fragmentation in mitochondrial genomes in relation to elevated sequence divergence and extreme rearrangements

Background

A single circular mitochondrial (mt) genome is a common feature across most metazoans. The mt-genome includes protein-coding genes involved in oxidative phosphorylation, as well as RNAs necessary for translation of mt-RNAs, whose order and number are highly conserved across animal clades, with few known exceptions of alternative mt-gene order or mt-genome architectures. One such exception consists of the fragmented mitochondrial genome, a type of genome architecture where mt-genes are split across two or more mt-chromosomes. However, the origins of mt-genome fragmentation and its effects on mt-genome evolution are unknown. Here, we investigate these origin and potential mechanisms underlying mt-genome fragmentation, focusing on a genus of booklice, Liposcelis, which exhibits elevated sequence divergence, frequent rearrangement of mt-gene order, and fragmentation of the mt genome, and compare them to other Metazoan clades.

Results

We found this genus Liposcelis exhibits very low conservation of mt-gene order across species, relative to other metazoans. Levels of gene order rearrangement were, however, unrelated to whether or not mt-genomes were fragmented or intact, suggesting mitochondrial genome fragmentation is not affecting mt-gene order directly. We further investigated possible mechanisms underpinning these patterns and revealed very high conservation of non-coding sequences at the edges of multiple recombination regions across populations of one particular Liposcelis species, supportive of a hypothesis that mt-fragmentation arises from recombination errors between mt-genome copies. We propose these errors may arise as a consequence of a heightened mutation rate in clades exhibiting mt-fragmentation. Consistent with this, we observed a striking pattern across three Metazoan phyla (Arthropoda, Nematoda, Cnidaria) characterised by members exhibiting high levels of mt-gene order rearrangement and cases of mt-fragmentation, whereby the mt-genomes of species more closely related to species with fragmented mt-genomes diverge more rapidly despite experiencing strong purifying selection.

Conclusions

We showed that contrary to expectations, mt-genome fragmentation is not correlated with the increase in mt-genome rearrangements. Furthermore, we present evidence that fragmentation of the mt-genome may be part of a general relaxation of a natural selection on the mt-genome, thus providing new insights into the origins of mt-genome fragmentation and evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-021-01218-7.

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.

HERMES: An improved method to test mitochondrial genome molecular synapomorphies among clades

Mitochondrial chromosomes have diversified among eukaryotes and many different architectures and features are now acknowledged for this genome. Here we present the improved HERMES index, which can measure and quantify the amount of molecular change experienced by mitochondrial genomes. We test the improved approach with ten molecular phylogenetic studies based on complete mitochondrial genomes, representing six bilaterian Phyla. In most cases, HERMES analysis spotted out clades or single species with peculiar molecular synapomorphies, allowing to identify phylogenetic and ecological patterns. The software presented herein handles linear, circular, and multi-chromosome genomes, thus widening the HERMES scope to the complete eukaryotic domain.

Mitochondrial genomes of Columbicola feather lice are highly fragmented, indicating repeated evolution of minicircle-type genomes in parasitic lice

Most animals have a conserved mitochondrial genome structure composed of a single chromosome. However, some organisms have their mitochondrial genes separated on several smaller circular or linear chromosomes. Highly fragmented circular chromosomes (“minicircles”) are especially prevalent in parasitic lice (Insecta: Phthiraptera), with 16 species known to have between nine and 20 mitochondrial minicircles per genome. All of these species belong to the same clade (mammalian lice), suggesting a single origin of drastic fragmentation. Nevertheless, other work indicates a lesser degree of fragmentation (2–3 chromosomes/genome) is present in some avian feather lice (Ischnocera: Philopteridae). In this study, we tested for minicircles in four species of the feather louse genus Columbicola (Philopteridae). Using whole genome shotgun sequence data, we applied three different bioinformatic approaches for assembling the Columbicola mitochondrial genome. We further confirmed these approaches by assembling the mitochondrial genome of Pediculus humanus from shotgun sequencing reads, a species known to have minicircles. Columbicola spp. genomes are highly fragmented into 15–17 minicircles between ∼1,100 and ∼3,100 bp in length, with 1–4 genes per minicircle. Subsequent annotation of the minicircles indicated that tRNA arrangements of minicircles varied substantially between species. These mitochondrial minicircles for species of Columbicola represent the first feather lice (Philopteridae) for which minicircles have been found in a full mitochondrial genome assembly. Combined with recent phylogenetic studies of parasitic lice, our results provide strong evidence that highly fragmented mitochondrial genomes, which are otherwise rare across the Tree of Life, evolved multiple times within parasitic lice.

Dicyemida and Orthonectida: Two Stories of Body Plan Simplification

Two enigmatic groups of morphologically simple parasites of invertebrates, the Dicyemida (syn. Rhombozoa) and the Orthonectida, since the 19th century have been usually considered as two classes of the phylum Mesozoa. Early molecular evidence suggested their relationship within the Spiralia (=Lophotrochozoa), however, high rates of dicyemid and orthonectid sequence evolution led to contradicting phylogeny reconstructions. Genomic data for orthonectids revealed that they are highly simplified spiralians and possess a reduced set of genes involved in metazoan development and body patterning. Acquiring genomic data for dicyemids, however, remains a challenge due to complex genome rearrangements including chromatin diminution and generation of extrachromosomal circular DNAs, which are reported to occur during the development of somatic cells. We performed genomic sequencing of one species of Dicyema, and obtained transcriptomic data for two Dicyema spp. Homeodomain (homeobox) transcription factors, G-protein-coupled receptors, and many other protein families have undergone a massive reduction in dicyemids compared to other animals. There is also apparent reduction of the bilaterian gene complements encoding components of the neuromuscular systems. We constructed and analyzed a large dataset of predicted orthologous proteins from three species of Dicyema and a set of spiralian animals including the newly sequenced genome of the orthonectid Intoshia linei. Bayesian analyses recovered the orthonectid lineage within the Annelida. In contrast, dicyemids form a separate clade with weak affinity to the Rouphozoa (Platyhelminthes plus Gastrotricha) or (Entoprocta plus Cycliophora) suggesting that the historically proposed Mesozoa is a polyphyletic taxon. Thus, dramatic simplification of body plans in dicyemids and orthonectids, as well as their intricate life cycles that combine metagenesis and heterogony, evolved independently in these two lineages.

Mitochondrial Genome Fragmentation Unites the Parasitic Lice of Eutherian Mammals

Organelle genome fragmentation has been found in a wide range of eukaryotic lineages; however, its use in phylogenetic reconstruction has not been demonstrated. We explored the use of mitochondrial (mt) genome fragmentation in resolving the controversial suborder-level phylogeny of parasitic lice (order Phthiraptera). There are approximately 5000 species of parasitic lice in four suborders (Amblycera, Ischnocera, Rhynchophthirina, and Anoplura), which infest mammals and birds. The phylogenetic relationships among these suborders are unresolved despite decades of studies. We sequenced the mt genomes of eight species of parasitic lice and compared them with 17 other species of parasitic lice sequenced previously. We found that the typical single-chromosome mt genome is retained in the lice of birds but fragmented into many minichromosomes in the lice of eutherian mammals. The shared derived feature of mt genome fragmentation unites the eutherian mammal lice of Ischnocera (family Trichodectidae) with Anoplura and Rhynchophthirina to the exclusion of the bird lice of Ischnocera (family Philopteridae). The novel clade, namely Mitodivisia, is also supported by phylogenetic analysis of mt genome and cox1 gene sequences. Our results demonstrate, for the first time, that organelle genome fragmentation is informative for resolving controversial high-level phylogenies.

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 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.

Mitochondrial genome diversity in dagger and needle nematodes (Nematoda: Longidoridae)

Dagger and needle nematodes included in the family Longidoridae (viz. Longidorus, Paralongidorus, and Xiphinema) are highly polyphagous plant-parasitic nematodes in wild and cultivated plants and some of them are plant-virus vectors (nepovirus). The mitochondrial (mt) genomes of the dagger and needle nematodes, Xiphinema rivesi, Xiphinema pachtaicum, Longidorus vineacola and Paralongidorus litoralis were sequenced in this study. The four circular mt genomes have an estimated size of 12.6, 12.5, 13.5 and 12.7 kb, respectively. Up to date, the mt genome of X. pachtaicum is the smallest genome found in Nematoda. The four mt genomes contain 12 protein-coding genes (viz. cox1-3, nad1-6, nad4L, atp6 and cob) and two ribosomal RNA genes (rrnL and rrnS), but the atp8 gene was not detected. These mt genomes showed a gene arrangement very different within the Longidoridae species sequenced, with the exception of very closely related species (X. americanum and X. rivesi). The sizes of non-coding regions in the Longidoridae nematodes were very small and were present in a few places in the mt genome. Phylogenetic analysis of all coding genes showed a closer relationship between Longidorus and Paralongidorus and different phylogenetic possibilities for the three Xiphinema species.

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.

Selfish Mitochondrial DNA Proliferates and Diversifies in Small, but not Large, Experimental Populations of Caenorhabditis briggsae

Evolutionary interactions across levels of biological organization contribute to a variety of fundamental processes including genome evolution, reproductive mode transitions, species diversification, and extinction. Evolutionary theory predicts that so-called “selfish” genetic elements will proliferate when the host effective population size (N_e) is small, but direct tests of this prediction remain few. We analyzed the evolutionary dynamics of deletion-containing mitochondrial DNA (ΔmtDNA) molecules, previously characterized as selfish elements, in six different natural strains of the nematode Caenorhabditis briggsae allowed to undergo experimental evolution in a range of population sizes (N = 1, 10, 100, and 1,000) for a maximum of 50 generations. Mitochondrial DNA (mtDNA) was analyzed for replicate lineages at each five-generation time point. Ten different ΔmtDNA molecule types were observed and characterized across generations in the experimental populations. Consistent with predictions from evolutionary theory, lab lines evolved in small-population sizes (e.g., nematode N = 1) were more susceptible to accumulation of high levels of preexisting ΔmtDNA compared with those evolved in larger populations. New ΔmtDNA elements were observed to increase in frequency and persist across time points, but almost exclusively at small population sizes. In some cases, ΔmtDNA levels decreased across generations when population size was large (nematode N = 1,000). Different natural strains of C. briggsae varied in their susceptibilities to ΔmtDNA accumulation, owing in part to preexisting compensatory mtDNA alleles in some strains that prevent deletion formation. This analysis directly demonstrates that the evolutionary trajectories of ΔmtDNA elements depend upon the population-genetic environments and molecular-genetic features of their hosts.

First comparative insight into the architecture of COI mitochondrial minicircle molecules of dicyemids reveals marked inter-species variation

Dicyemids, poorly known parasites of benthic cephalopods, are one of the few phyla in which mitochondrial (mt) genome architecture departs from the typical ~16 kb circular metazoan genome. In addition to a putative circular genome, a series of mt minicircles that each comprises the mt encoded units (I-III) of the cytochrome c oxidase complex have been reported. Whether the structure of the mt minicircles is a consistent feature among dicyemid species is unknown. Here we analyse the complete cytochrome c oxidase subunit I (COI) minicircle molecule, containing the COI gene and an associated non-coding region (NCR), for ten dicyemid species, allowing for first time comparisons between species of minicircle architecture, NCR function and inferences of minicircle replication. Divergence in COI nucleotide sequences between dicyemid species was high (average net divergence = 31.6%) while within species diversity was lower (average net divergence = 0.2%). The NCR and putative 5' section of the COI gene were highly divergent between dicyemid species (average net nucleotide divergence of putative 5' COI section = 61.1%). No tRNA genes were found in the NCR, although palindrome sequences with the potential to form stem-loop structures were identified in some species, which may play a role in transcription or other biological processes.

Mitochondrial genome plasticity among species of the nematode genus Meloidogyne (Nematoda: Tylenchina)

The mitochondrial (mt) genomes of the plant-parasitic root-knot nematodes Meloidogyne arenaria, Meloidogyne enterolobii and Meloidogyne javanica were sequenced and compared with those of three other root-knot nematode species in order to explore the mt genome plasticity within Meloidogyne. The mt genomes of M. arenaria, M. enterolobii and M. javanica are circular, with an estimated size of 18.8, 18.9 and 19.6 kb, respectively. Compared to other nematodes these mt genomes are larger, due to the presence of large non-coding regions. The mt genome architecture within the genus Meloidogyne varied in the position of trn genes and in the position, length and nucleotide composition of non-coding regions. These variations were observed independent of the species' natural environments or reproductive modes. M. enterolobii showed three main non-coding regions whereas Meloidogyne chitwoodi, Meloidogyne incognita, M. javanica and M. arenaria had two non-coding regions, and Meloidogyne graminicola had a unique large non-coding region interrupted by two trn genes. trn genes were positioned in different regions of the mt genomes in M. chitwoodi, M. enterolobii and M. graminicola, whereas the trn gene order was identical between M. arenaria, M. incognita and M. javanica. Importantly, M. graminicola had extra copies of trnV and trnS2. High divergence levels between the two copies of each trn might indicate duplication events followed by random loss and mutations in the anticodon. Tree-based methods based on amino acid sequences of 12 mt protein-coding genes support the monophyly for the tropical and mitotic parthenogenetic species, M. arenaria, M. enterolobii, M. incognita and M. javanica and for a clade that includes the meiotic parthenogenetic species, M. chitwoodi and M. graminicola. A comparison of the mt genome architecture in plant-parasitic nematodes and phylogenetic analyses support that Pratylenchus is the most recent ancestor of root-knot nematodes.

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 Meloidogyne graminicola (Tylenchina): a unique gene arrangement and its phylogenetic implications

Meloidogyne graminicola is one of the most economically important plant parasitic-nematodes (PPNs). In the present study, we determined the complete mitochondrial (mt) DNA genome sequence of this plant pathogen. Compared with other PPNs genera, this genome (19,589 bp) is only slightly smaller than that of Pratylenchus vulnus (21,656 bp). The nucleotide composition of the whole mtDNA sequence of M. graminicola is significantly biased toward A and T, with T being the most favored nucleotide and C being the least favored. The A+T content of the entire genome is 83.51%. The mt genome of M. graminicola contains 36 genes (lacking atp8) that are transcribed in the same direction. The gene arrangement of the mt genome of M. graminicola is unique. A total of 21 out of 22 tRNAs possess a DHU loop only, while tRNASer(AGN) lacks a DHU loop. The two large noncoding regions (2,031 bp and 5,063 bp) are disrupted by tRNASer(UCN). Phylogenetic analysis based on concatenated amino acid sequences of 12 protein-coding genes support the monophylies of the three orders Rhabditida, Mermithida and Trichinellida, the suborder Rhabditina and the three infraorders Spiruromorpha, Oxyuridomorpha and Ascaridomorpha, but do not support the monophylies of the two suborders Spirurina and Tylenchina, and the three infraorders Rhabditomorpha, Panagrolaimomorpha and Tylenchomorpha. The four Tylenchomorpha species including M. graminicola, P. vulnus, H. glycines and R. similis from the superfamily Tylenchoidea are placed within a well-supported monophyletic clade, but far from the other two Tylenchomorpha species B. xylophilus and B. mucronatus of Aphelenchoidea. In the clade of Tylenchoidea, M. graminicola is sister to P. vulnus, and H. glycines is sister to R. similis, which suggests root-knot nematodes has a closer relationship to Pratylenchidae nematodes than to cyst nematodes.

Mitochondrial genomes of Meloidogyne chitwoodi and M. incognita (Nematoda: Tylenchina): comparative analysis, gene order and phylogenetic relationships with other nematodes

Root-knot nematodes (Meloidogyne spp.) are among the most important plant pathogens. In this study, the mitochondrial (mt) genomes of the root-knot nematodes, M. chitwoodi and M. incognita were sequenced. PCR analyses suggest that both mt genomes are circular, with an estimated size of 19.7 and 18.6-19.1kb, respectively. The mt genomes each contain a large non-coding region with tandem repeats and the control region. The mt gene arrangement of M. chitwoodi and M. incognita is unlike that of other nematodes. Sequence alignments of the two Meloidogyne mt genomes showed three translocations; two in transfer RNAs and one in cox2. Compared with other nematode mt genomes, the gene arrangement of M. chitwoodi and M. incognita was most similar to Pratylenchus vulnus. Phylogenetic analyses (Maximum Likelihood and Bayesian inference) were conducted using 78 complete mt genomes of diverse nematode species. Analyses based on nucleotides and amino acids of the 12 protein-coding mt genes showed strong support for the monophyly of class Chromadorea, but only amino acid-based analyses supported the monophyly of class Enoplea. The suborder Spirurina was not monophyletic in any of the phylogenetic analyses, contradicting the Clade III model, which groups Ascaridomorpha, Spiruromorpha and Oxyuridomorpha based on the small subunit ribosomal RNA gene. Importantly, comparisons of mt gene arrangement and tree-based methods placed Meloidogyne as sister taxa of Pratylenchus, a migratory plant endoparasitic nematode, and not with the sedentary endoparasitic Heterodera. Thus, comparative analyses of mt genomes suggest that sedentary endoparasitism in Meloidogyne and Heterodera is based on convergent evolution.

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.

Recombination and evolution of duplicate control regions in the mitochondrial genome of the Asian big-headed turtle, Platysternon megacephalum

Complete mitochondrial (mt) genome sequences with duplicate control regions (CRs) have been detected in various animal species. In Testudines, duplicate mtCRs have been reported in the mtDNA of the Asian big-headed turtle, Platysternon megacephalum, which has three living subspecies. However, the evolutionary pattern of these CRs remains unclear. In this study, we report the completed sequences of duplicate CRs from 20 individuals belonging to three subspecies of this turtle and discuss the micro-evolutionary analysis of the evolution of duplicate CRs. Genetic distances calculated with MEGA 4.1 using the complete duplicate CR sequences revealed that within turtle subspecies, genetic distances between orthologous copies from different individuals were 0.63% for CR1 and 1.2% for CR2app:addword:respectively, and the average distance between paralogous copies of CR1 and CR2 was 4.8%. Phylogenetic relationships were reconstructed from the CR sequences, excluding the variable number of tandem repeats (VNTRs) at the 3' end using three methods: neighbor-joining, maximum likelihood algorithm, and Bayesian inference. These data show that any two CRs within individuals were more genetically distant from orthologous genes in different individuals within the same subspecies. This suggests independent evolution of the two mtCRs within each P. megacephalum subspecies. Reconstruction of separate phylogenetic trees using different CR components (TAS, CD, CSB, and VNTRs) suggested the role of recombination in the evolution of duplicate CRs. Consequently, recombination events were detected using RDP software with break points at ≈290 bp and ≈1,080 bp. Based on these results, we hypothesize that duplicate CRs in P. megacephalum originated from heterological ancestral recombination of mtDNA. Subsequent recombination could have resulted in homogenization during independent evolutionary events, thus maintaining the functions of duplicate CRs in the mtDNA of P. megacephalum.

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.

The multipartite mitochondrial genome of Liposcelis bostrychophila: insights into the evolution of mitochondrial genomes in bilateral animals

Booklice (order Psocoptera) in the genus Liposcelis are major pests to stored grains worldwide and are closely related to parasitic lice (order Phthiraptera). We sequenced the mitochondrial (mt) genome of Liposcelis bostrychophila and found that the typical single mt chromosome of bilateral animals has fragmented into and been replaced by two medium-sized chromosomes in this booklouse; each of these chromosomes has about half of the genes of the typical mt chromosome of bilateral animals. These mt chromosomes are 8,530 bp (mt chromosome I) and 7,933 bp (mt chromosome II) in size. Intriguingly, mt chromosome I is twice as abundant as chromosome II. It appears that the selection pressure for compact mt genomes in bilateral animals favors small mt chromosomes when small mt chromosomes co-exist with the typical large mt chromosomes. Thus, small mt chromosomes may have selective advantages over large mt chromosomes in bilateral animals. Phylogenetic analyses of mt genome sequences of Psocodea (i.e. Psocoptera plus Phthiraptera) indicate that: 1) the order Psocoptera (booklice and barklice) is paraphyletic; and 2) the order Phthiraptera (the parasitic lice) is monophyletic. Within parasitic lice, however, the suborder Ischnocera is paraphyletic; this differs from the traditional view that each suborder of parasitic lice is monophyletic.

Widespread atypical mitochondrial DNA structure in isopods (Crustacea, Peracarida) related to a constitutive heteroplasmy in terrestrial species

Metazoan mitochondrial DNA (mtDNA) is generally composed of circular monomeric molecules. However, a few exceptions do exist and among them two terrestrial isopods Armadillidium vulgare and Porcellionides pruinosus have an atypical mtDNA composed of linear monomers associated with circular "head-to-head" dimers: a very unusual structure for animal mtDNA genome. To assess the distribution of this atypical mtDNA among isopods, we performed RFLP and Southern blot analyses on mtDNA of 16 terrestrial (Oniscidea family) and two aquatic isopod species: the marine Sphaeroma serratum (suborder Flabellifera, sister group of Oniscidea) and the freshwater Asellus aquaticus (Asellota, early derived taxon of isopod). The atypical mtDNA structure was observed in 15 terrestrial isopod species and A. aquaticus, suggesting a wide distribution of atypical mtDNA among isopods. However, a typical metazoan mtDNA structure was detected in the marine isopod S. serratum and the Oniscidea Ligia oceanica . Our results suggest two possible scenarios: an early origin of the atypical mtDNA in isopods followed by reversion to the typical ancestral mtDNA structure for several species, or a convergent appearance of the atypical mtDNA structure in two isopod suborders. We compare this distribution of the atypical mtDNA structure with the presence of a heteroplasmy also observed in the mtDNA of several terrestrial isopod species. We discuss if this transmitted heteroplasmy is vectored by the atypical mtDNA and its impact on the maintenance of the atypical mtDNA in isopods.

Mitochondrial genome deletions and minicircles are common in lice (Insecta: Phthiraptera)

Background

The gene composition, gene order and structure of the mitochondrial genome are remarkably stable across bilaterian animals. Lice (Insecta: Phthiraptera) are a major exception to this genomic stability in that the canonical single chromosome with 37 genes found in almost all other bilaterians has been lost in multiple lineages in favour of multiple, minicircular chromosomes with less than 37 genes on each chromosome.

Results

Minicircular mt genomes are found in six of the ten louse species examined to date and three types of minicircles were identified: heteroplasmic minicircles which coexist with full sized mt genomes (type 1); multigene chromosomes with short, simple control regions, we infer that the genome consists of several such chromosomes (type 2); and multiple, single to three gene chromosomes with large, complex control regions (type 3). Mapping minicircle types onto a phylogenetic tree of lice fails to show a pattern of their occurrence consistent with an evolutionary series of minicircle types. Analysis of the nuclear-encoded, mitochondrially-targetted genes inferred from the body louse, Pediculus, suggests that the loss of mitochondrial single-stranded binding protein (mtSSB) may be responsible for the presence of minicircles in at least species with the most derived type 3 minicircles (Pediculus, Damalinia).

Conclusions

Minicircular mt genomes are common in lice and appear to have arisen multiple times within the group. Life history adaptive explanations which attribute minicircular mt genomes in lice to the adoption of blood-feeding in the Anoplura are not supported by this expanded data set as minicircles are found in multiple non-blood feeding louse groups but are not found in the blood-feeding genus Heterodoxus. In contrast, a mechanist explanation based on the loss of mtSSB suggests that minicircles may be selectively favoured due to the incapacity of the mt replisome to synthesize long replicative products without mtSSB and thus the loss of this gene lead to the formation of minicircles in lice.

The mitochondrial genome of the soybean cyst nematode, Heterodera glycines

We sequenced the entire coding region of the mitochondrial genome of Heterodera glycines. The sequence obtained comprised 14.9 kb, with PCR evidence indicating that the entire genome comprised a single, circular molecule of approximately 21-22 kb. The genome is the most T-rich nematode mitochondrial genome reported to date, with T representing over half of all nucleotides on the coding strand. The genome also contains the highest number of poly(T) tracts so far reported (to our knowledge), with 60 poly(T) tracts ≥ 12 Ts. All genes are transcribed from the same mitochondrial strand. The organization of the mitochondrial genome of H. glycines shows a number of similarities compared with Radopholus similis, but fewer similarities when compared with Meloidogyne javanica. Very few gene boundaries are shared with Globodera pallida or Globodera rostochiensis. Partial mitochondrial genome sequences were also obtained for Heterodera cardiolata (5.3 kb) and Punctodera chalcoensis (6.8 kb), and these had identical organizations compared with H. glycines. We found PCR evidence of a minicircular mitochondrial genome in P. chalcoensis, but at low levels and lacking a noncoding region. Such circularised genome fragments may be present at low levels in a range of nematodes, with multipartite mitochondrial genomes representing a shift to a condition in which these subgenomic circles predominate.

Poly(T) variation in heteroderid nematode mitochondrial genomes is predominantly an artefact of amplification

We assessed the rate of in vitro polymerase errors at polythymidine [poly(T)] tracts in the mitochondrial DNA (mtDNA) of a heteroderid nematode (Heterodera cajani). The mtDNA of these nematodes contain unusually high numbers of poly(T) tracts, and have previously been suggested to contain biological poly(T) length variation. However, using a cloned molecule, we observed that poly(T) variation was generated in vitro at regions containing more than six consecutive Ts. This artefactual error rate was estimated at 7.3 × 10(-5) indels/poly(T) tract >6 Ts/cycle. This rate was then compared to the rate of poly(T) variation detected after the amplification of a biological sample, in order to estimate the 'biological + artefactual' rate of poly(T) variation. There was no significant difference between the artefactual and the artefactual + biological rates, suggesting that the majority of poly(T) variation in the biological sample was artefactual. We then examined the generation of poly(T) variation in a range of templates with tracts up to 16 Ts long, utilizing a range of Heteroderidae species. We observed that T deletions occurred five times more frequently than insertions, and a trend towards increasing error rates with increasing poly(T) tract length. These findings have significant implications for studies involving genomes with many homopolymer tracts.

MITOCHONDRIAL DNA IN THE OOGAMOCHLAMYS CLADE (CHLOROPHYCEAE): HIGH GC CONTENT AND UNIQUE GENOME ARCHITECTURE FOR GREEN ALGAE(1)

Most mitochondrial genomes in the green algal phylum Chlorophyta are AT-rich, circular-mapping DNA molecules. However, mitochondrial genomes from the Reinhardtii clade of the Chlorophyceae lineage are linear and sometimes fragmented into subgenomic forms. Moreover, Polytomella capuana, from the Reinhardtii clade, has an elevated GC content (57.2%). In the present study, we examined mitochondrial genome conformation and GC bias in the Oogamochlamys clade of the Chlorophyceae, which phylogenetic data suggest is closely related to the Reinhardtii clade. Total DNA from selected Oogamochlamys taxa, including four Lobochlamys culleus (H. Ettl) Pröschold, B. Marin, U. G. Schlöss. et Melkonian strains, Lobochlamys segnis (H. Ettl) Pröschold, B. Marin, U. G. Schlöss. et Melkonian, and Oogamochlamys gigantea (O. Dill) Pröschold, B. Marin, U. G. Schlöss. et Melkonian, was subjected to Southern blot analyses with cob and cox1 probes, and the results suggest that the mitochondrial genome of these taxa is represented by multiple-sized linear DNA fragments with overlapping homologies. On the basis of these data, we propose that linear mitochondrial DNA with a propensity to become fragmented arose in an ancestor common to the Reinhardtii and Oogamochlamys clades or even earlier in the evolutionary history of the Chlorophyceae. Analyses of partial cob and cox1 sequences from these Oogamochlamys taxa revealed an unusually high GC content (49.9%-65.1%) and provided evidence for the accumulation of cob and cox1 pseudogenes and truncated sequences in the mitochondrial genome of all L. culleus strains examined.

Complete mitochondrial genome sequences of the three pelagic chaetognaths Sagitta nagae, Sagitta decipiens and Sagitta enflata

The complete nucleotide sequences of the mitochondrial genomes were determined for the three pelagic chaetognaths, Sagitta nagae, Sagitta decipiens, and Sagitta enflata. The mitochondrial genomes of these species which were 11,459, 11,121, and 12,631bp in length, respectively, contained 14 genes (11 protein-coding genes, one transfer RNA gene, and two ribosomal RNA genes), and were found to have lost 23 genes that are present in the typical metazoan mitochondrial genome. The same mitochondrial genome contents have been reported from the benthic chaetognaths belonging to the family Spadellidae, Paraspadella gotoi and Spadella cephaloptera. Within the phylum Chaetognatha, Sagitta and Spadellidae are distantly related, suggesting that the gene loss occurred in the ancestral species of the phylum. The gene orders of the three Sagitta species are markedly different from those of the other non-Chaetognatha metazoans. In contrast to the region with frequent gene rearrangements, no gene rearrangements were observed in the gene cluster encoding COII-III, ND1-3, srRNA, and tRNA(met). Within this conserved gene cluster, gene rearrangements were not observed in the three Sagitta species or between the Sagitta and Spadellidae species. The gene order of this cluster was also assumed to be the ancestral state of the phylum.

A unique genetic code change in the mitochondrial genome of the parasitic nematode Radopholus similis

BACKGROUND

Mitochondria (mt) contain their own autonomously replicating DNA, constituted as a small circular genome encoding essential subunits of the respiratory chain. Mt DNA is characterized by a genetic code which differs from the standard one. Interestingly, the mt genome of nematodes share some peculiar features, such as small transfer RNAs, truncated ribosomal RNAs and - in the class of Chromadorean nematodes - unidirectional transcription.

FINDINGS

We present the complete mt genomic sequence (16,791 bp) of the plant-parasitic nematode Radopholus similis (class Chromadorea). Although it has a gene content similar to most other nematodes, many idiosyncrasies characterize the extremely AT-rich mt genome of R. similis (85.4% AT). The secondary structure of the large (16S) rRNA is further reduced, the gene order is unique, the large non-coding region contains two large repeats, and most interestingly, the UAA codon is reassigned from translation termination to tyrosine. In addition, 7 out of 12 protein-coding genes lack a canonical stop codon and analysis of transcriptional data showed the absence of polyadenylation. Northern blot analysis confirmed that only one strand is transcribed and processed. Furthermore, using nucleotide content bias methods, regions for the origin of replication are suggested.

CONCLUSION

The extraordinary mt genome of R. similis with its unique genetic code appears to contain exceptional features correlated to DNA decoding. Therefore the genome may provide an incentive to further elucidate these barely understood processes in nematodes. This comprehension may eventually lead to parasitic nematode-specific control targets as healthy mitochondria are imperative for organism survival. In addition, the presented genome is an interesting exceptional event in genetic code evolution.

A fragmented metazoan organellar genome: the two mitochondrial chromosomes of Hydra magnipapillata

Background

Animal mitochondrial (mt) genomes are characteristically circular molecules of ~16–20 kb. Medusozoa (Cnidaria excluding Anthozoa) are exceptional in that their mt genomes are linear and sometimes subdivided into two to presumably four different molecules. In the genus Hydra, the mt genome comprises one or two mt chromosomes. Here, we present the whole mt genome sequence from the hydrozoan Hydra magnipapillata, comprising the first sequence of a fragmented metazoan mt genome encoded on two linear mt chromosomes (mt1 and mt2).

Results

The H. magnipapillata mt chromosomes contain the typical metazoan set of 13 genes for respiratory proteins, the two rRNA genes and two tRNA genes. All genes are unidirectionally oriented on mt1 and mt2, and several genes overlap. The gene arrangement suggests that the two mt chromosomes originated from one linear molecule that separated between nd5 and rns. Strong correlations between the AT content of rRNA genes (rns and rnl) and the AT content of protein-coding genes among 24 cnidarian genomes imply that base composition is mainly determined by mt genome-wide constraints. We show that identical inverted terminal repeats (ITR) occur on both chromosomes; these ITR contain a partial copy or part of the 3' end of cox1 (54 bp). Additionally, both mt chromosomes possess identical oriented sequences (IOS) at the 5' and 3' ends (5' and 3' IOS) adjacent to the ITR. The 5' IOS contains trnM and non-coding sequences (119 bp), whereas the 3' IOS comprises a larger part (mt2) with a larger partial copy of cox1 (243 bp).

Conclusion

ITR are also documented in the two other available medusozoan mt genomes (Aurelia aurita and Hydra oligactis). In H. magnipapillata, the arrangement of ITR and 5' IOS and 3' IOS suggest that these regions are crucial for mt DNA replication and/or transcription initiation. An analogous organization occurs in a highly fragmented ichthyosporean mt genome. With our data, we can reject a model of mt replication that has previously been proposed for Hydra. This raises new questions regarding replication mechanisms probably employed by all medusozoans, and also has general implications for the expected organization of fragmented linear mt chromosomes of other taxa.

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

The monogonont rotifer Brachionus plicatilis is an emerging model system for a diverse array of questions in limnological ecosystem dynamics, the evolution of sexual recombination, cryptic speciation, and the phylogeny of basal metazoans. We sequenced the complete mitochondrial genome of B. plicatilis sensu strictu NH1L and found that it is composed of 2 circular chromosomes, designated mtDNA-I (11,153 bp) and mtDNA-II (12,672 bp). Hybridization to DNA isolated from mitochondria demonstrated that mtDNA-I is present at 4 times the copy number of mtDNA-II. The only nucleotide similarity between the 2 chromosomes is a 4.9-kbp region of 99.5% identity including a transfer RNA (tRNA) gene and an extensive noncoding region that contains putative D-loop and control sequence. The mtDNA-I chromosome encodes 4 proteins (ATP6, COB, NAD1, and NAD2), 13 tRNAs, and the large and small subunit ribosomal RNAs; mtDNA-II encodes 8 proteins (COX1-3, NAD3-6, and NAD4L) and 9 tRNAs. Gene order is not conserved between B. plicatilis and its closest relative with a sequenced mitochondrial genome, the acanthocephalan Leptorhynchoides thecatus, or other sequenced mitochondrial genomes. Polymerase chain reaction assays and Southern hybridization to DNA from 18 strains of Brachionus suggest that the 2-chromosome structure has been stable for millions of years. The novel organization of the B. plicatilis mitochondrial genome into 2 nearly equal chromosomes of 4-fold different copy number may provide insight into the evolution of metazoan mitochondria and the phylogenetics of rotifers and other basal animal phyla.

Poly(T) variation within mitochondrial protein-coding genes in Globodera (Nematoda: Heteroderidae)

We sequenced a mitochondrial subgenome from the nematode Globodera rostochiensis, in two overlapping pieces. The subgenome was 9210 bp and contained four protein-coding genes (ND4, COIII, ND3, Cytb) and two tRNA genes (tRNA(Thr), tRNA(Gln)). Genome organization was similar to that of Globodera pallida, which is multipartite. Together with the small number of genes on this subgenome, this suggests that the mitochondrial genome of G. rostochiensis is also multipartite. In the initial clones sequenced, COIII and ND3 were full-length, while ND4 and Cytb were interrupted by premature stop codons and contained point indels that disrupted the reading frame. However, sequencing of multiple clones, from DNA extracted both from multiple individuals and from single cysts, revealed a predominant source of variation-in the length of polythymidine tracts. Comparison of our genomic sequences with ESTs similarly revealed variation in the length of polythymidine tracts. We subsequently sequenced both genomic DNA and mRNA from populations of G. pallida. In each case, variation in the length of polythymidine tracts was observed. The levels of expression of mitochondrial genes in G. pallida were representative of the subgenomes present: little evidence of differential expression was observed. These observations are consistent with the operation of posttranscriptional editing in Globodera mitochondria, although this is difficult to show conclusively in the presence of intraindividual gene sequence variation. Further, alternative explanations cannot be discounted; these include the operation of slippage during translation or that genomic copies of most genes are pseudogenes with a small proportion of full-length sequences able to maintain mitochondrial function.