The Mitochondrial Subgenomes of the Nematode Globodera pallida Are Mosaics: Evidence of Recombination in an Animal Mitochondrial Genome

Invited Review Beyond parasitic convergence: unravelling the evolution of the organellar genomes in holoparasites

BACKGROUND

The molecular evolution of organellar genomes in angiosperms has been studied extensively, with some lineages, such as parasitic ones, displaying unique characteristics. Parasitism has emerged 12 times independently in angiosperm evolution. Holoparasitism is the most severe form of parasitism, and is found in ~10 % of parasitic angiosperms. Although a few holoparasitic species have been examined at the molecular level, most reports involve plastomes instead of mitogenomes. Parasitic plants establish vascular connections with their hosts through haustoria to obtain water and nutrients, which facilitates the exchange of genetic information, making them more susceptible to horizontal gene transfer (HGT). HGT is more prevalent in the mitochondria than in the chloroplast or nuclear compartments.

SCOPE

This review summarizes current knowledge on the plastid and mitochondrial genomes of holoparasitic angiosperms, compares the genomic features across the different lineages, and discusses their convergent evolutionary trajectories and distinctive features. We focused on Balanophoraceae (Santalales), which exhibits extraordinary traits in both their organelles.

CONCLUSIONS

Apart from morphological similarities, plastid genomes of holoparasitic plants also display other convergent features, such as rampant gene loss, biased nucleotide composition and accelerated evolutionary rates. In addition, the plastomes of Balanophoraceae have extremely low GC and gene content, and two unexpected changes in the genetic code. Limited data on the mitochondrial genomes of holoparasitic plants preclude thorough comparisons. Nonetheless, no obvious genomic features distinguish them from the mitochondria of free-living angiosperms, except for a higher incidence of HGT. HGT appears to be predominant in holoparasitic angiosperms with a long-lasting endophytic stage. Among the Balanophoraceae, mitochondrial genomes exhibit disparate evolutionary paths with notable levels of heteroplasmy in Rhopalocnemis and unprecedented levels of HGT in Lophophytum. Despite their differences, these Balanophoraceae share a multichromosomal mitogenome, a feature also found in a few free-living angiosperms.

The minicircular and extremely heteroplasmic mitogenome of the holoparasitic plant Rhopalocnemis phalloides

The plastid and nuclear genomes of parasitic plants exhibit deeply altered architectures,^1-13 whereas the few examined mitogenomes range from deeply altered to conventional.^14-20 To provide further insight on mitogenome evolution in parasitic plants, we report the highly modified mitogenome of Rhopalocnemis phalloides, a holoparasite in Balanophoraceae. Its mitogenome is uniquely arranged in 21 minicircular chromosomes that vary in size from 4,949 to 7,861 bp, with a total length of only 130,713 bp. All chromosomes share an identical 896 bp conserved region, with a large stem-loop that acts as the origin of replication, flanked on each side by hypervariable and semi-conserved regions. Similar minicircular structures with shared and unique regions have been observed in parasitic animals and free-living protists,^21-24 suggesting convergent structural evolution. Southern blots confirm both the minicircular structure and the replication origin of the mitochondrial chromosomes. PacBio reads provide evidence for chromosome recombination and rolling-circle replication for the R. phalloides mitogenome. Despite its small size, the mitogenome harbors a typical set of genes and introns within the unique regions of each chromosome, yet introns are the smallest among seed plants and ferns. The mitogenome also exhibits extreme heteroplasmy, predominantly involving short indels and more complex variants, many of which cause potential loss-of-function mutations for some gene copies. All heteroplasmic variants are transcribed, and functional and nonfunctional protein-coding variants are spliced and RNA edited. Our findings offer a unique perspective into how mitogenomes of parasitic plants can be deeply altered and shed light on plant mitogenome replication.

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.

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.

Clonorchis sinensis and Clonorchiasis: The Relevance of Exploring Genetic Variation

Parasitic trematodes (flukes) cause substantial mortality and morbidity in humans. The Chinese liver fluke, Clonorchis sinensis, is one of the most destructive parasitic worms in humans in China, Vietnam, Korea and the Russian Far East. Although C. sinensis infection can be controlled relatively well using anthelmintics, the worm is carcinogenic, inducing cholangiocarcinoma and causing major suffering in ~15 million people in Asia. This chapter provides an account of C. sinensis and clonorchiasis research-covering aspects of biology, epidemiology, pathogenesis and immunity, diagnosis, treatment and control, genetics and genomics. It also describes progress in the area of molecular biology (genetics, genomics, transcriptomics and proteomics) and highlights challenges associated with comparative genomics and population genetics. It then reviews recent advances in the sequencing and characterisation of the mitochondrial and nuclear genomes for a Korean isolate of C. sinensis and summarises salient comparative genomic work and the implications thereof. The chapter concludes by considering how advances in genomic and informatics will enable research on the genetics of C. sinensis and related parasites, as well as the discovery of new fluke-specific intervention targets.

Evolution and inheritance of animal mitochondrial DNA: rules and exceptions

Mitochondrial DNA (mtDNA) has been studied intensely for “its own” merit. Its role for the function of the cell and the organism remains a fertile field, its origin and evolution is an indispensable part of the evolution of life and its interaction with the nuclear DNA is among the most important cases of genome synergism and co-evolution. Also, mtDNA was proven one of the most useful tools in population genetics and molecular phylogenetics. In this article we focus on animal mtDNA and discuss briefly how our views about its structure, function and transmission have changed, how these changes affect the information we have accumulated through its use in the fields of phylogeny and population structure and what are the most important questions that remain open for future research.

A metagenetic approach to determine the diversity and distribution of cyst nematodes at the level of the country, the field and the individual

Distinct populations of the potato cyst nematode (PCN) Globodera pallida exist in the UK that differ in their ability to overcome various sources of resistance. An efficient method for distinguishing between populations would enable pathogen‐informed cultivar choice in the field. Science and Advice for Scottish Agriculture (SASA) annually undertake national DNA diagnostic tests to determine the presence of PCN in potato seed and ware land by extracting DNA from soil floats. These DNA samples provide a unique resource for monitoring the distribution of PCN and further interrogation of the diversity within species.

We identify a region of mitochondrial DNA descriptive of three main groups of G. pallida present in the UK and adopt a metagenetic approach to the sequencing and analysis of all SASA samples simultaneously. Using this approach, we describe the distribution of G. pallida mitotypes across Scotland with field‐scale resolution. Most fields contain a single mitotype, one‐fifth contain a mix of mitotypes, and less than 3% contain all three mitotypes. Within mixed fields, we were able to quantify the relative abundance of each mitotype across an order of magnitude. Local areas within mixed fields are dominated by certain mitotypes and indicate towards a complex underlying ‘pathoscape’. Finally, we assess mitotype distribution at the level of the individual cyst and provide evidence of ‘hybrids’. This study provides a method for accurate, quantitative and high‐throughput typing of up to one thousand fields simultaneously, while revealing novel insights into the national genetic variability of an economically important plant parasite.

The genomic basis of parasitism in the Strongyloides clade of nematodes

Soil-transmitted nematodes, including the Strongyloides genus, cause one of the most prevalent neglected tropical diseases. Here we compare the genomes of four Strongyloides species, including the human pathogen Strongyloides stercoralis, and their close relatives that are facultatively parasitic (Parastrongyloides trichosuri) and free-living (Rhabditophanes sp. KR3021). A significant paralogous expansion of key gene families--families encoding astacin-like and SCP/TAPS proteins--is associated with the evolution of parasitism in this clade. Exploiting the unique Strongyloides life cycle, we compare the transcriptomes of the parasitic and free-living stages and find that these same gene families are upregulated in the parasitic stages, underscoring their role in nematode parasitism.

Mitogenomics reveals high synteny and long evolutionary histories of sympatric cryptic nematode species

Species with seemingly identical morphology but with distinct genetic differences are abundant in the marine environment and frequently co-occur in the same habitat. Such cryptic species are typically delineated using a limited number of mitochondrial and/or nuclear marker genes, which do not yield information on gene order and gene content of the genomes under consideration. We used next-generation sequencing to study the composition of the mitochondrial genomes of four sympatrically distributed cryptic species of the Litoditis marina species complex (PmI, PmII, PmIII, and PmIV). The ecology, biology, and natural occurrence of these four species are well known, but the evolutionary processes behind this cryptic speciation remain largely unknown. The gene order of the mitochondrial genomes of the four species was conserved, but differences in genome length, gene length, and codon usage were observed. The atp8 gene was lacking in all four species. Phylogenetic analyses confirm that PmI and PmIV are sister species and that PmIII diverged earliest. The most recent common ancestor of the four cryptic species was estimated to have diverged 16 MYA. Synonymous mutations outnumbered nonsynonymous changes in all protein-encoding genes, with the Complex IV genes (coxI-III) experiencing the strongest purifying selection. Our mitogenomic results show that morphologically similar species can have long evolutionary histories and that PmIII has several differences in genetic makeup compared to the three other species, which may explain why it is better adapted to higher temperatures than the other species.

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.

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

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.

Coexistence of minicircular and a highly rearranged mtDNA molecule suggests that recombination shapes mitochondrial genome organization

Recombination has been proposed as a possible mechanism to explain mitochondrial (mt) gene rearrangements, although the issue of whether mtDNA recombination occurs in animals has been controversial. In this study, we sequenced the entire mt genome of the megaspilid wasp Conostigmus sp., which possessed a highly rearranged mt genome. The sequence of the A+T-rich region contained a number of different types of repeats, similar to those reported previously in the nematode Meloidogyne javanica, in which recombination was discovered. In Conostigmus, we detected the end products of recombination: a range of minicircles. However, using isolated (cloned) fragments of the A+T-rich region, we established that some of these minicircles were found to be polymerase chain reaction (PCR) artifacts. It appears that regions with repeats are prone to PCR template switching or PCR jumping. Nevertheless, there is strong evidence that one minicircle is real, as amplification primers that straddle the putative breakpoint junction produce a single strong amplicon from genomic DNA but not from the cloned A+T-rich region. The results provide support for the direct link between recombination and mt gene rearrangement. Furthermore, we developed a model of recombination which is important for our understanding of mtDNA evolution.

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.

Comparison of transcript profiles in different life stages of the nematode Globodera pallida under different host potato genotypes

The potato cyst nematodes (PCNs) Globodera pallida and Globodera rostochiensis are important parasites of potato. PCNs undergo complex biotrophic interactions with their hosts that involve gene expression changes in both the nematode and the host plant. The aim of this study was to determine key genes that are differentially expressed in Globodera pallida life cycle stages and during the initiation of the feeding site in susceptible and partially resistant potato genotypes. For this purpose, two microarray experiments were designed: (i) a comparison of eggs, infective second-stage juveniles (J2s) and sedentary parasitic-stage J2s (SJ2); (ii) a comparison of SJ2s at 8 days after inoculation (DAI) in the susceptible cultivar (Desirée) and two partially resistant lines. The results showed differential expression of G. pallida genes during the stages studied, including previously characterized effectors. In addition, a large number of genes changed their expression between SJ2s in the susceptible cultivar and those infecting partially resistant lines; the number of genes with modified expression was lower when the two partially resistant lines were compared. Moreover, a histopathological study was performed at several time points (7, 14 and 30 DAI) and showed the similarities between both partially resistant lines with a delay and degeneration in the formation of the syncytia in comparison with the susceptible cultivar. Females at 30 DAI in partially resistant lines showed a delay in their development in comparison with those in the susceptible cultivar.

Paternal leakage of mitochondrial DNA in experimental crosses of populations of the potato cyst nematode Globodera pallida

Animal mtDNA is typically assumed to be maternally inherited. Paternal mtDNA has been shown to be excluded from entering the egg or eliminated post-fertilization in several animals. However, in the contact zones of hybridizing species and populations, the reproductive barriers between hybridizing organisms may not be as efficient at preventing paternal mtDNA inheritance, resulting in paternal leakage. We assessed paternal mtDNA leakage in experimental crosses of populations of a cyst-forming nematode, Globodera pallida. A UK population, Lindley, was crossed with two South American populations, P5A and P4A. Hybridization of these populations was supported by evidence of nuclear DNA from both the maternal and paternal populations in the progeny. To assess paternal mtDNA leakage, a ~3.4 kb non-coding mtDNA region was analyzed in the parental populations and in the progeny. Paternal mtDNA was evident in the progeny of both crosses involving populations P5A and P4A. Further, paternal mtDNA replaced the maternal mtDNA in 22 and 40 % of the hybrid cysts from these crosses, respectively. These results indicate that under appropriate conditions, paternal leakage occurs in the mtDNA of parasitic nematodes, and supports the hypothesis that hybrid zones facilitate paternal leakage. Thus, assumptions of strictly maternal mtDNA inheritance may be frequently violated, particularly when divergent populations interbreed.

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.

Distinctive mitochondrial genome of Calanoid copepod Calanus sinicus with multiple large non-coding regions and reshuffled gene order: useful molecular markers for phylogenetic and population studies

Background

Copepods are highly diverse and abundant, resulting in extensive ecological radiation in marine ecosystems. Calanus sinicus dominates continental shelf waters in the northwest Pacific Ocean and plays an important role in the local ecosystem by linking primary production to higher trophic levels. A lack of effective molecular markers has hindered phylogenetic and population genetic studies concerning copepods. As they are genome-level informative, mitochondrial DNA sequences can be used as markers for population genetic studies and phylogenetic studies.

Results

The mitochondrial genome of C. sinicus is distinct from other arthropods owing to the concurrence of multiple non-coding regions and a reshuffled gene arrangement. Further particularities in the mitogenome of C. sinicus include low A + T-content, symmetrical nucleotide composition between strands, abbreviated stop codons for several PCGs and extended lengths of the genes atp6 and atp8 relative to other copepods. The monophyletic Copepoda should be placed within the Vericrustacea. The close affinity between Cyclopoida and Poecilostomatoida suggests reassigning the latter as subordinate to the former. Monophyly of Maxillopoda is rejected. Within the alignment of 11 C. sinicus mitogenomes, there are 397 variable sites harbouring three 'hotspot' variable sites and three microsatellite loci.

Conclusion

The occurrence of the circular subgenomic fragment during laboratory assays suggests that special caution should be taken when sequencing mitogenomes using long PCR. Such a phenomenon may provide additional evidence of mitochondrial DNA recombination, which appears to have been a prerequisite for shaping the present mitochondrial profile of C. sinicus during its evolution. The lack of synapomorphic gene arrangements among copepods has cast doubt on the utility of gene order as a useful molecular marker for deep phylogenetic analysis. However, mitochondrial genomic sequences have been valuable markers for resolving phylogenetic issues concerning copepods. The variable site maps of C. sinicus mitogenomes provide a solid foundation for population genetic studies.

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.

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

Parasitism genes and host range disparities in biotrophic nematodes: the conundrum of polyphagy versus specialisation

This essay considers biotrophic cyst and root-knot nematodes in relation to their biology, host-parasite interactions and molecular genetics. These nematodes have to face the biological consequences of the physical constraints imposed by the soil environment in which they live while their hosts inhabit both above and below ground environments. The two groups of nematodes appear to have adopted radically different solutions to these problems with the result that one group is a host specialist and reproduces sexually while the other has an enormous host range and reproduces by mitotic parthenogenesis. We consider what is known about the modes of parasitism used by these nematodes and how it relates to their host range, including the surprising finding that parasitism genes in both nematode groups have been recruited from bacteria. The nuclear and mitochondrial genomes of these two nematode groups are very different and we consider how these findings relate to the biology of the organisms.

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.

Thermodynamic characterization of specific interactions between the human Lon protease and G-quartet DNA

Lon is an ATP-powered protease that binds DNA. However, the function of DNA binding by Lon remains elusive. Studies suggest that human Lon (hLon) binds preferentially to a G-rich single-stranded DNA (ssDNA) sequence overlapping the light strand promoter of mitochondrial DNA. This sequence is contained within a 24-base oligonucleotide referred to as LSPas. Here, we use biochemical and biophysical approaches to elucidate the structural properties of ssDNAs bound by hLon, as well as the thermodynamics of DNA binding by hLon. Electrophoretic mobility shift assay and circular dichroism show that ssDNAs with a propensity for forming parallel G-quartets are specifically bound by hLon. Isothermal titration calorimetry demonstrates that hLon binding to LSPas is primarily driven by enthalpy change associated with a significant reduction in heat capacity. Differential scanning calorimetry pinpoints an excess heat capacity upon hLon binding to LSPas. By contrast, hLon binding to an 8-base G-rich core sequence is entropically driven with a relatively negligible change in heat capacity. A considerable enhancement of thermal stability accompanies hLon binding to LSPas as compared to the G-rich core. Taken together, these data support the notion that hLon binds G-quartets through rigid-body binding and that binding to LSPas is coupled with structural adaptation.

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

Recently, a multipartite mitochondrial genome was characterized in the potato cyst nematode, Globodera pallida. Six subgenomic circles were detectable by PCR, while full-length genomes were not. We investigate here whether this subgenomic organization occurs in a close relative of G. pallida. We amplified and sequenced one entire mitochondrial subgenome from the cyst-forming nematode, Globodera rostochiensis. Comparison of the noncoding region of this subgenome with those reported previously for G. pallida facilitated the design of amplification primers for a range of subgenomes from G. rostochiensis. We then randomly sequenced five subgenomic fragments, each representative of a unique subgenome. This study indicates that the multipartite structure reported for G. pallida is conserved in G. rostochiensis. A comparison of subgenomic organization between these two Globodera species indicates a considerable degree of overlap between them. Indeed, we identify two subgenomes with an organization identical with that reported for G. pallida. However, other subgenomes are unique to G. rostochiensis, although some of these have blocks of genes comparable to those in G. pallida. Dot-plot comparisons of pairs of subgenomes from G. rostochiensis indicate that the different subgenomes share fragments with high sequence identity. We interpret this as evidence that recombination is operating in the mitochondria of G. rostochiensis.