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

Characterization of the fragmented mitochondrial genome of domestic pig louse Haematopinus suis (Insecta: Haematopinidae) from China

The domestic pig louse Haematopinus suis (Linnaeus, 1758) (Phthiraptera: Anoplura) is a common ectoparasite of domestic pigs, which can act as a vector of various infectious disease agents. Despite its significance, the molecular genetics, biology and systematics of H. suis from China have not been studied in detail. In the present study, the entire mitochondrial (mt) genome of H. suis isolate from China was sequenced and compared with that of H. suis isolate from Australia. We identified 37 mt genes located on nine circular mt minichromosomes, 2.9 kb-4.2 kb in size, each containing 2-8 genes and one large non-coding region (NCR) (1,957 bp-2,226 bp). The number of minichromosomes, gene content, and gene order in H. suis isolates from China and Australia are identical. Total sequence identity across coding regions was 96.3% between H. suis isolates from China and Australia. For the 13 protein-coding genes, sequence differences ranged from 2.8%-6.5% consistent nucleotides with amino acids. Our result is H. suis isolates from China and Australia being the same H. suis species. The present study determined the entire mt genome of H. suis from China, providing additional genetic markers for studying the molecular genetics, biology and systematics of domestic pig louse.

The fragmented mitochondrial genomes of two Linognathus lice reveal active minichromosomal recombination and recombination hotspots

Evidence for recombination between mitochondrial (mt) minichromosomes has been reported in sucking lice, but it is still not clear how frequent mt minichromosomal recombination occurs. We sequenced the mt genomes of the cattle louse Linognathus vituli and the goat louse L. africanus. Both Linognathus species have 10 mt minichromosomes, and seven of them have the same gene content and gene arrangement. Comparison of mt karyotypes revealed numerous inter-minichromosomal recombination events in the evolution of Linognathus species. Minichromosome merger, gene duplication and gene translocation occurred in the lineage leading to Linognathus lice. After the divergence of L. vituli and L. africanus, duplication, degeneration, deletion and translocation of genes also occurred independently in each species. Most of the recombination events in the Linognathus species occurred upstream of either cox3 or nad2, indicating these two locations were hotspots for inter-minichromosomal recombination. Our results provide an important perspective on mt genome evolution in metazoans.

Mitochondrial Genome Fragmentation Occurred Multiple Times Independently in Bird Lice of the Families Menoponidae and Laemobothriidae

Mitochondrial (mt) genome fragmentation has been discovered in all five parvorders of parasitic lice (Phthiraptera). To explore whether minichromosomal characters derived from mt genome fragmentation are informative for phylogenetic studies, we sequenced the mt genomes of 17 species of bird lice in Menoponidae and Laemobothriidae (Amblycera). Four species of Menoponidae (Actornithophilus sp. 1 ex [pied oystercatcher], Act. sp. 2 ex [masked lapwing], Austromenopon sp. 2 ex [sooty tern and crested tern], Myr. sp. 1 ex [satin bowerbird]) have fragmented mt genomes, whereas the other 13 species retain the single-chromosome mt genomes. The two Actornithophilus species have five and six mt minichromosomes, respectively. Aus. sp. 2 ex [sooty tern and crested tern] has two mt minichromosomes, in contrast to Aus. sp. 1 ex [sooty shearwater], which has a single mt chromosome. Myr. sp. 1 ex [satin bowerbird] has four mt minichromosomes. When mapped on the phylogeny of Menoponidae and Laemobothriidae, it is evident that mt genome fragmentation has occurred multiple times independently among Menoponidae and Laemobothriidae species. We found derived mt minichromosomal characters shared between Myrsidea species, between Actornithophilus species, and between and among different ischnoceran genera, respectively. We conclude that while mt genome fragmentation as a general feature does not unite all the parasitic lice that have this feature, each independent mt genome fragmentation event does produce minichromosomal characters that can be informative for phylogenetic studies of parasitic lice at different taxonomic levels.

A Novel Mitochondrial Genome Fragmentation Pattern in the Buffalo Louse Haematopinus tuberculatus (Psocodea: Haematopinidae)

Sucking lice are obligate ectoparasites of mammalian hosts, causing serious public health problems and economic losses worldwide. It is well known that sucking lice have fragmented mitochondrial (mt) genomes, but many remain undetermined. To better understand patterns of mt genome fragmentation in the sucking lice, we sequenced the mt genome of the buffalo louse Haematopinus tuberculatus using next-generation sequencing (NGS). The mt genome of H. tuberculatus has ten circular minichromosomes containing a total of 37 genes. Each minichromosome is 2.9-5.0 kb long and carries one to eight genes plus one large non-coding region. The number of mt minichromosomes of H. tuberculatus (ten) is different from those of congeneric species (horse louse H. asini, domestic pig louse H. suis and wild pig louse H. apri) and other sucking lice. Two events (gene translocation and merger of mt minichromosome) are observed in Haematopinus. Compared to other studies, our phylogeny generated from mt genome datasets showed a different topology, suggesting that inclusion of data other than mt genomes would be required to resolve phylogeny of sucking lice. To our knowledge, this is the first report of a ten mt minichromosomes genome in sucking lice, which opens a new outlook into unexplored mt genome fragmentation patterns in sucking lice.

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.

Fragmented mitochondrial genomes of seal lice (family Echinophthiriidae) and gorilla louse (family Pthiridae): frequent minichromosomal splits and a host switch of lice between seals

Background

The mitochondrial (mt) genomes of 15 species of sucking lice from seven families have been studied to date. These louse species have highly dynamic, fragmented mt genomes that differ in the number of minichromosomes, the gene content, and gene order in a minichromosome between families and even between species of the same genus.

Results

In the present study, we analyzed the publicly available data to understand mt genome fragmentation in seal lice (family Echinophthiriidae) and gorilla louse, Pthirus gorillae (family Pthiridae), in particular the role of minichromosome split and minichromosome merger in the evolution of fragmented mt genomes. We show that 1) at least three ancestral mt minichromosomes of sucking lice have split in the lineage leading to seal lice, 2) one minichromosome ancestral to primate lice has split in the lineage to the gorilla louse, and 3) two ancestral minichromosomes of seal lice have merged in the lineage to the northern fur seal louse. Minichromosome split occurred 15-16 times in total in the lineages leading to species in six families of sucking lice investigated. In contrast, minichromosome merger occurred only four times in the lineages leading to species in three families of sucking lice. Further, three ancestral mt minichromosomes of sucking lice have split multiple times independently in different lineages of sucking lice. Our analyses of mt karyotypes and gene sequences also indicate the possibility of a host switch of crabeater seal louse to Weddell seals.

Conclusions

We conclude that: 1) minichromosome split contributes more than minichromosome merger in mt genome fragmentation of sucking lice, and 2) mt karyotype comparison helps understand the phylogenetic relationships between sucking louse species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08530-8.

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.

Frequent tRNA gene translocation towards the boundaries with control regions contributes to the highly dynamic mitochondrial genome organization of the parasitic lice of mammals

Background

The typical single-chromosome mitochondrial (mt) genome of animals has fragmented into multiple minichromosomes in the lineage Mitodivisia, which contains most of the parasitic lice of eutherian mammals. These parasitic lice differ from each other even among congeneric species in mt karyotype, i.e. the number of minichromosomes, and the gene content and gene order in each minichromosome, which is in stark contrast to the extremely conserved single-chromosome mt genomes across most animal lineages. How fragmented mt genomes evolved is still poorly understood. We use Polyplax sucking lice as a model to investigate how tRNA gene translocation shapes the dynamic mt karyotypes.

Results

We sequenced the full mt genome of the Asian grey shrew louse, Polyplax reclinata. We then inferred the ancestral mt karyotype for Polyplax lice and compared it with the mt karyotypes of the three Polyplax species sequenced to date. We found that tRNA genes were entirely responsible for mt karyotype variation among these three species of Polyplax lice. Furthermore, tRNA gene translocation observed in Polyplax lice was only between different types of minichromosomes and towards the boundaries with the control region. A similar pattern of tRNA gene translocation can also been seen in other sucking lice with fragmented mt genomes.

Conclusions

We conclude that inter-minichromosomal tRNA gene translocation orientated towards the boundaries with the control region is a major contributing factor to the highly dynamic mitochondrial genome organization in the parasitic lice of mammals.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07859-w.

Mitochondrial genome of Geomydoecus aurei, a pocket-gopher louse

Parasitic lice demonstrate an unusual array of mitochondrial genome architectures and gene arrangements. We characterized the mitochondrial genome of Geomydoecus aurei, a chewing louse (Phthiraptera: Trichodectidae) found on pocket gophers (Rodentia: Geomyidae) using reads from both Illumina and Oxford Nanopore sequencing coupled with PCR, cloning, and Sanger sequencing to verify structure and arrangement for each chromosome. The genome consisted of 12 circular mitochondrial chromosomes ranging in size from 1,318 to 2,088 nucleotides (nt). Total genome size was 19,015 nt. All 37 genes typical of metazoans (2 rRNA genes, 22 tRNA genes, and 13 protein-coding genes) were present. An average of 26% of each chromosome was composed of non-gene sequences. Within the non-gene region of each chromosome, there was a 79-nt nucleotide sequence that was identical among chromosomes and a conserved sequence with secondary structure that was always followed by a poly-T region. We hypothesize that these regions may be important in the initiation of transcription and DNA replication, respectively. The G. aurei genome shares 8 derived gene clusters with other chewing lice of mammals, but in G. aurei, genes on several chromosomes are not contiguous.

Highly rearranged mitochondrial genome in Falcolipeurus lice (Phthiraptera: Philopteridae) from endangered eagles

Background

Fragmented mitochondrial (mt) genomes and extensive mt gene rearrangements have been frequently reported from parasitic lice (Insecta: Phthiraptera). However, relatively little is known about the mt genomes from the family Philopteridae, the most species-rich family within the suborder Ischnocera.

Methods

Herein, we use next-generation sequencing to decode the mt genome of Falcolipeurus suturalis and compare it with the mt genome of F. quadripustulatus. Phylogenetic relationships within the family Philopteridae were inferred from the concatenated 13 protein-coding genes of the two Falcolipeurus lice and members of the family Philopteridae using Bayesian inference (BI) and maximum likelihood (ML) methods.

Results

The complete mt genome of F. suturalis is a circular, double-stranded DNA molecule 16,659bp in size that contains 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and three non-coding regions. The gene order of the F. suturalis mt genome is rearranged relative to that of F. quadripustulatus, and is radically different from both other louse species and the putative ancestral insect. Phylogenetic analyses revealed clear genetic distinctiveness between F. suturalis and F. quadripustulatus (Bayesian posterior probabilities=1.0 and bootstrapping frequencies=100), and that the genus Falcolipeurus is sister to the genus Ibidoecus (Bayesian posterior probabilities=1.0 and bootstrapping frequencies=100).

Conclusions

These datasets help to better understand gene rearrangements in lice and the phylogenetic position of Falcolipeurus and provide useful genetic markers for systematic studies of bird lice.

Graphic abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04776-5.

A novel fragmented mitochondrial genome in the protist pathogen Toxoplasma gondii and related tissue coccidia

Mitochondrial genome content and structure vary widely across the eukaryotic tree of life, with protists displaying extreme examples. Apicomplexan and dinoflagellate protists have evolved highly reduced mitochondrial genome sequences, mtDNA, consisting of only three cytochrome genes and fragmented rRNA genes. Here, we report the independent evolution of fragmented cytochrome genes in Toxoplasma and related tissue coccidia and evolution of a novel genome architecture consisting minimally of 21 sequence blocks (SBs) totaling 5.9 kb that exist as nonrandom concatemers. Single-molecule Nanopore reads consisting entirely of SBs ranging from 0.1 to 23.6 kb reveal both whole and fragmented cytochrome genes. Full-length cytochrome transcripts including a divergent coxIII are detected. The topology of the mitochondrial genome remains an enigma. Analysis of a cob point mutation reveals that homoplasmy of SBs is maintained. Tissue coccidia are important pathogens of man and animals, and the mitochondrion represents an important therapeutic target. The mtDNA sequence has been elucidated, but a definitive genome architecture remains elusive.

New Records of Bartonella spp. and Rickettsia spp. in Lice Collected from Small Rodents

Lice are blood-sucking insects that are of medical and veterinary significance as parasites and vectors for various infectious agents. More than half of described blood-sucking lice species are found on rodents. Rodents are important hosts of several Bartonella and Rickettsia species, and some of these bacteria are characterized as human pathogens in Europe. Rodent ectoparasites, such as fleas and ticks, are important vectors of Bartonella spp. and Rickettsia spp., but knowledge about the presence of these bacteria in lice is limited. The aim of this study was to determine the prevalence of Bartonella and Rickettsia bacteria in lice collected from rodents in Slovakia. The ectoparasites were collected from small rodents captured from 2010 to 2015 at four different sites in eastern Slovakia. The presence of Bartonella and Rickettsia species in lice samples was screened by real-time PCR, targeting ssrA and gltA genes, respectively. The molecular characterization of the Bartonella strains was based on sequence analysis of partial rpoB and intergenic spacer (ITS) genes, and of the Rickettsia species on sequence analysis of the gltA gene. A total of 1074 lice of seven species were collected from six rodent species. Bartonella DNA was detected in Hoplopleura affinis (collected from Apodemus agrarius, Apodemus flavicollis, and Myodes glareolus), Polyplax serrata (from A. agrarius), and Hoplopleura sp. (from A. flavicollis). Sequence analysis revealed that the Bartonella strains belonged to the Bartonella coopersplainsensis, Bartonella tribocorum, and Bartonella taylorii genogroups. Rickettsia DNA was detected in H. affinis and P. serrata collected from A. agrarius. Sequence analysis revealed two Rickettsia species: Rickettsia helvetica and Rickettsia sp. The results of the study confirm the presence of Bartonella spp. and Rickettsia spp. in lice collected from rodents.

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.

Structure, gene order, and nucleotide composition of mitochondrial genomes in parasitic lice from Amblycera

Parasitic lice have unique mitochondrial (mt) genomes characterized by rearranged gene orders, variable genome structures, and less AT content compared to most other insects. However, relatively little is known about the mt genomes of Amblycera, the suborder sister to all other parasitic lice. Comparing among nine different genera (including representative of all seven families), we show that Amblycera have variable and highly rearranged mt genomes. Some genera have fragmented genomes that vary considerably in length, whereas others have a single mt chromosome. Notably, these genomes are more AT-biased than most other lice. We also recover genus-level phylogenetic relationships among Amblycera that are consistent with those reported from large nuclear datasets, indicating that mt sequences are reliable for reconstructing evolutionary relationships in Amblycera. However, gene order data cannot reliably recover these same relationships. Overall, our results suggest that the mt genomes of lice, already know to be distinctive, are even more variable than previously thought.

Variation of mitochondrial minichromosome composition in Hoplopleura lice (Phthiraptera: Hoplopleuridae) from rats

Background

The family Hoplopleuridae contains at least 183 species of blood-sucking lice, which widely parasitize both mice and rats. Fragmented mitochondrial (mt) genomes have been reported in two rat lice (Hoplopleura kitti and H. akanezumi) from this family, but some minichromosomes were unidentified in their mt genomes.

Methods

We sequenced the mt genome of the rat louse Hoplopleura sp. with an Illumina platform and compared its mt genome organization with H. kitti and H. akanezumi.

Results

Fragmented mt genome of the rat louse Hoplopleura sp. contains 37 genes which are on 12 circular mt minichromosomes. Each mt minichromosome is 1.8–2.7 kb long and contains 1–5 genes and one large non-coding region. The gene content and arrangement of mt minichromosomes of Hoplopleura sp. (n = 3) and H. kitti (n = 3) are different from those in H. akanezumi (n = 3). Phylogenetic analyses based on the deduced amino acid sequences of the eight protein-coding genes showed that the Hoplopleura sp. was more closely related to H. akanezumi than to H. kitti, and then they formed a monophyletic group.

Conclusions

Comparison among the three rat lice revealed variation in the composition of mt minichromosomes within the genus Hoplopleura. Hoplopleura sp. is the first species from the family Hoplopleuridae for which a complete fragmented mt genome has been sequenced. The new data provide useful genetic markers for studying the population genetics, molecular systematics and phylogenetics of blood-sucking lice.

"Parasite turnover zone" at secondary contact: A new pattern in host-parasite population genetics

We describe here a new pattern of population genetic structure in a host-parasite system that can arise after secondary contact of previously isolated populations. Due to different generation times, and therefore different tempos of molecular evolution, the host and parasite populations reach different degrees of genetic differentiation during their separation (e.g., in refugia). Consequently, upon secondary contact, the host populations are able to re-establish a single panmictic population across the area of contact, while the parasite populations stop their dispersal at the secondary contact zone and create a narrow hybrid zone. From the host's perspective, the parasite's hybrid zone functions on a microevolutionary scale as a "parasite turnover zone": while the hosts are passing from area A to area B, their parasites turn genetically from the area A genotypes to the area B genotypes. We demonstrate this novel pattern with a model composed of Apodemus mice and Polyplax lice by comparing maternally inherited markers (complete mitochondrial genomes, and complete genomes of the vertically transmitted symbiont Legionella polyplacis) with single nucleotide polymorphisms derived from louse genomic data. We discuss the circumstances that may lead to this pattern and possible reasons why it has been overlooked in studies of host-parasite population genetics.

Fragmented mitochondrial genomes evolved in opposite directions between closely related macaque louse Pedicinus obtusus and colobus louse Pedicinus badii

We report for the first time the fragmented mitochondrial (mt) genomes of two Pedicinus species: Pedicinus obtusus and Pedicinus badii, and compared them with the lice of humans and chimpanzees. Despite being congeneric, the two monkey lice are distinct from each other in mt karyotype. The variation in mt karyotype between the two Pedicinus lice is the most pronounced among the congeneric species of sucking lice observed to date and is attributable to the opposite directions between them in mt karyotype evolution. Two of the inferred ancestral mt minichromosomes of the higher primate lice merged as one in the macaque louse whereas one of the ancestral minichromosomes split into two in the colobus louse after these two species diverged from their most recent common ancestor. Our results showed that mt genome fragmentation was a two-way process in the higher primate lice, and minichromosome merger was more common than previously thought.

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.

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 Mitochondrial Genome of the Guanaco Louse, Microthoracius praelongiceps: Insights into the Ancestral Mitochondrial Karyotype of Sucking Lice (Anoplura, Insecta)

Fragmented mitochondrial (mt) genomes have been reported in 11 species of sucking lice (suborder Anoplura) that infest humans, chimpanzees, pigs, horses, and rodents. There is substantial variation among these lice in mt karyotype: the number of minichromosomes of a species ranges from 9 to 20; the number of genes in a minichromosome ranges from 1 to 8; gene arrangement in a minichromosome differs between species, even in the same genus. We sequenced the mt genome of the guanaco louse, Microthoracius praelongiceps, to help establish the ancestral mt karyotype for sucking lice and understand how fragmented mt genomes evolved. The guanaco louse has 12 mt minichromosomes; each minichromosome has 2-5 genes and a non-coding region. The guanaco louse shares many features with rodent lice in mt karyotype, more than with other sucking lice. The guanaco louse, however, is more closely related phylogenetically to human lice, chimpanzee lice, pig lice, and horse lice than to rodent lice. By parsimony analysis of shared features in mt karyotype, we infer that the most recent common ancestor of sucking lice, which lived ∼75 Ma, had 11 minichromosomes; each minichromosome had 1-6 genes and a non-coding region. As sucking lice diverged, split of mt minichromosomes occurred many times in the lineages leading to the lice of humans, chimpanzees, and rodents whereas merger of minichromosomes occurred in the lineage leading to the lice of pigs and horses. Together, splits and mergers of minichromosomes created a very complex and dynamic mt genome organization in the sucking lice.

The Highly Divergent Mitochondrial Genomes Indicate That the Booklouse, Liposcelis bostrychophila (Psocoptera: Liposcelididae) Is a Cryptic Species

The booklouse, Liposcelis bostrychophila is an important storage pest worldwide. The mitochondrial (mt) genome of an asexual strain (Beibei, China) of the L. bostrychophila comprises two chromosomes; each chromosome contains approximate half of the 37 genes typically found in bilateral animals. The mt genomes of two sexual strains of L. bostrychophila, however, comprise five and seven chromosomes, respectively; each chromosome contains one to six genes. To understand mt genome evolution in L. bostrychophila, and whether L. bostrychophila is a cryptic species, we sequenced the mt genomes of six strains of asexual L. bostrychophila collected from different locations in China, Croatia, and the United States. The mt genomes of all six asexual strains of L. bostrychophila have two chromosomes. Phylogenetic analysis of mt genome sequences divided nine strains of L. bostrychophila into four groups. Each group has a distinct mt genome organization and substantial sequence divergence (48.7–87.4%) from other groups. Furthermore, the seven asexual strains of L. bostrychophila, including the published Beibei strain, are more closely related to two other species of booklice, L. paeta and L. sculptilimacula, than to the sexual strains of L. bostrychophila. Our results revealed highly divergent mt genomes in the booklouse, L. bostrychophila, and indicate that L. bostrychophila is a cryptic species.

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.

The first complete genome sequences of the acI lineage, the most abundant freshwater Actinobacteria, obtained by whole-genome-amplification of dilution-to-extinction cultures

The acI lineage of the phylum Actinobacteria is the most abundant bacterial group in most freshwater lakes. However, due to difficulties in laboratory cultivation, only two mixed cultures and some incomplete single-amplified or metagenome-derived genomes have been reported for the lineage. Here, we report the initial cultivation and complete genome sequences of four novel strains of the acI lineage from the tribes acI-A1, -A4, -A7, and -C1. The acI strains, initially isolated by dilution-to-extinction culturing, eventually failed to be maintained as axenic cultures. However, the first complete genomes of the acI lineage were successfully obtained from these initial cultures through whole genome amplification applied to more than hundreds of cultured acI cells. The genome sequences exhibited features of genome streamlining and showed that the strains are aerobic chemoheterotrophs sharing central metabolic pathways, with some differences among tribes that may underlie niche diversification within the acI lineage. Actinorhodopsin was found in all strains, but retinal biosynthesis was complete in only A1 and A4 tribes.

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

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

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

The Mitochondrial Cytochrome Oxidase Subunit I Gene Occurs on a Minichromosome with Extensive Heteroplasmy in Two Species of Chewing Lice, Geomydoecus aurei and Thomomydoecus minor

In animals, mitochondrial DNA (mtDNA) typically occurs as a single circular chromosome with 13 protein-coding genes and 22 tRNA genes. The various species of lice examined previously, however, have shown mitochondrial genome rearrangements with a range of chromosome sizes and numbers. Our research demonstrates that the mitochondrial genomes of two species of chewing lice found on pocket gophers, Geomydoecus aurei and Thomomydoecus minor, are fragmented with the 1,536 base-pair (bp) cytochrome-oxidase subunit I (cox1) gene occurring as the only protein-coding gene on a 1,916-1,964 bp minicircular chromosome in the two species, respectively. The cox1 gene of T. minor begins with an atypical start codon, while that of G. aurei does not. Components of the non-protein coding sequence of G. aurei and T. minor include a tRNA (isoleucine) gene, inverted repeat sequences consistent with origins of replication, and an additional non-coding region that is smaller than the non-coding sequence of other lice with such fragmented mitochondrial genomes. Sequences of cox1 minichromosome clones for each species reveal extensive length and sequence heteroplasmy in both coding and noncoding regions. The highly variable non-gene regions of G. aurei and T. minor have little sequence similarity with one another except for a 19-bp region of phylogenetically conserved sequence with unknown function.

Fragmented mitochondrial genomes in two suborders of parasitic lice of eutherian mammals (Anoplura and Rhynchophthirina, Insecta)

Parasitic lice (order Phthiraptera) infest birds and mammals. The typical animal mitochondrial (mt) genome organization, which consists of a single chromosome with 37 genes, was found in chewing lice in the suborders Amblycera and Ischnocera. The sucking lice (suborder Anoplura) known, however, have fragmented mt genomes with 9–20 minichromosomes. We sequenced the mt genome of the elephant louse, Haematomyzus elephantis – the first species of chewing lice investigated from the suborder Rhynchophthirina. We identified 33 mt genes in the elephant louse, which were on 10 minichromosomes. Each minichromosome is 3.5–4.2 kb in size and has 2–6 genes. Phylogenetic analyses of mt genome sequences confirm that the elephant louse is more closely related to sucking lice than to the chewing lice in the Amblycera and Ischnocera. Our results indicate that mt genome fragmentation is shared by the suborders Anoplura and Rhynchophthirina. Nine of the 10 mt minichromosomes of the elephant louse differ from those of the sucking lice (Anoplura) known in gene content and gene arrangement, indicating that distinct mt karyotypes have evolved in Anoplura and Rhynchophthirina since they diverged ~92 million years ago.

Review: can diet influence the selective advantage of mitochondrial DNA haplotypes?

This review explores the potential for changes in dietary macronutrients to differentially influence mitochondrial bioenergetics and thereby the frequency of mtDNA haplotypes in natural populations. Such dietary modification may be seasonal or result from biogeographic or demographic shifts. Mechanistically, mtDNA haplotypes may influence the activity of the electron transport system (ETS), retrograde signalling to the nuclear genome and affect epigenetic modifications. Thus, differential provisioning by macronutrients may lead to selection through changes in the levels of ATP production, modulation of metabolites (including AMP, reactive oxygen species (ROS) and the NAD⁺/NADH ratio) and potentially complex epigenetic effects. The exquisite complexity of dietary influence on haplotype frequency is further illustrated by the fact that macronutrients may differentially influence the selective advantage of specific mutations in different life-history stages. In Drosophila, complex I mutations may affect larval growth because dietary nutrients are fed through this complex in immaturity. In contrast, the majority of electrons are provided to complex III in adult flies. We conclude the review with a case study that considers specific interactions between diet and complex I of the ETS. Complex I is the first enzyme of the mitochondrial ETS and co-ordinates in the oxidation of NADH and transfer of electrons to ubiquinone. Although the supposition that mtDNA variants may be selected upon by dietary macronutrients could be intuitively consistent to some and counter intuitive to others, it must face a multitude of scientific hurdles before it can be recognized.

The mitochondrial genome of the chimpanzee louse, Pediculus schaeffi: insights into the process of mitochondrial genome fragmentation in the blood-sucking lice of great apes

Background

Blood-sucking lice in the genera Pediculus and Pthirus are obligate ectoparasites of great apes. Unlike most bilateral animals, which have 37 mitochondrial (mt) genes on a single circular chromosome, the sucking lice of humans have extensively fragmented mt genomes. The head louse, Pediculus capitis, and the body louse, Pe. humanus, have their 37 mt genes on 20 minichromosomes. The pubic louse, Pthirus pubis, has its 34 mt genes known on 14 minichromosomes. To understand the process of mt genome fragmentation in the sucking lice of great apes, we sequenced the mt genome of the chimpanzee louse, Pe. schaeffi, and compared it with the three human lice.

Results

We identified all of the 37 mt genes typical of bilateral animals in the chimpanzee louse; these genes are on 18 types of minichromosomes. Seventeen of the 18 minichromosomes of the chimpanzee louse have the same gene content and gene arrangement as their counterparts in the human head louse and the human body louse. However, five genes, cob, trnS₁, trnN, trnE and trnM, which are on three minichromosomes in the human head louse and the human body louse, are together on one minichromosome in the chimpanzee louse.

Conclusions

Using the human pubic louse, Pt. pubis, as an outgroup for comparison, we infer that a single minichromosome has fragmented into three in the lineage leading to the human head louse and the human body louse since this lineage diverged from the chimpanzee louse ~6 million years ago. Our results provide insights into the process of mt genome fragmentation in the sucking lice in a relatively fine evolutionary scale.

Electronic supplementary material

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

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.

Fragmented mitochondrial genomes are present in both major clades of the blood-sucking lice (suborder Anoplura): evidence from two Hoplopleura rodent lice (family Hoplopleuridae)

Background

The suborder Anoplura contains 540 species of blood-sucking lice that parasitize over 840 species of eutherian mammals. Fragmented mitochondrial (mt) genomes have been found in the lice of humans, pigs, horses and rats from four families: Pediculidae, Pthiridae, Haematopinidae and Polyplacidae. These lice, eight species in total, are from the same major clade of the Anoplura. The mt genomes of these lice consist of 9–20 minichromosomes; each minichromosome is 1.5–4 kb in size and has 1–8 genes. To understand mt genome fragmentation in the other major clade of the Anoplura, we sequenced the mt genomes of two species of rodent lice in the genus Hoplopleura (family Hoplopleuridae).

Results

We identified 28 mt genes on 10 minichromosomes in the mouse louse, Ho. akanezumi; each minichromosome is 1.7–2.7 kb long and has 1–6 genes. We identified 34 mt genes on 11 minichromosomes in the rat louse, Ho. kitti; each minichromosome is 1.8–2.8 kb long and has 1–5 genes. Ho. akanezumi also has a chimeric minichromosome with parts of two rRNA genes and a full-length tRNA gene for tyrosine. These two rodent lice share the same pattern for the distribution of all of the protein-coding and rRNA genes but differ in tRNA gene content and gene arrangement in four minichromosomes. Like the four genera of blood-sucking lice that have been investigated in previous studies, the Hoplopleura species have four minichromosomes that are only found in this genus.

Conclusions

Our results indicate that fragmented mt genomes were present in the most recent common ancestor of the two major clades of the blood-sucking lice, which lived ~75 million years ago. Intra-genus variation in the pattern of mt genome fragmentation is common in the blood-sucking lice (suborder Anoplura) and genus-specific minichromosomes are potential synapomorphies. Future studies should expand into more species, genera and families of blood-sucking lice to explore further the phylogenetic utility of the novel features associated with fragmented mt genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-751) contains supplementary material, which is available to authorized users.

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.

Variation in mitochondrial minichromosome composition between blood-sucking lice of the genus Haematopinus that infest horses and pigs

BACKGROUND

The genus Haematopinus contains 21 species of blood-sucking lice, parasitizing both even-toed ungulates (pigs, cattle, buffalo, antelopes, camels and deer) and odd-toed ungulates (horses, donkeys and zebras). The mitochondrial genomes of the domestic pig louse, Haematopinus suis, and the wild pig louse, Haematopinus apri, have been sequenced recently; both lice have fragmented mitochondrial genomes with 37 genes on nine minichromosomes. To understand whether the composition of mitochondrial minichromosomes and the gene content and gene arrangement of each minichromosome are stable within the genus, we sequenced the mitochondrial genome of the horse louse, Haematopinus asini.

METHODS

We used a PCR-based strategy to amplify four mitochondrial minichromosomes in near full-length, and then amplify the entire coding regions of all of the nine mitochondrial minichromosomes of the horse louse. These amplicons were sequenced with an Illumina Hiseq platform.

RESULTS

We identified all of the 37 mitochondrial genes typical of bilateral animals in the horse louse, Haematopinus asini; these genes are on nine circular minichromosomes. Each minichromosome is 3.5-5.0 kb in size and consists of a coding region and a non-coding region except R-nad4L-rrnS-C minichromosome, which contains two coding regions and two non-coding regions. Six of the nine minichromosomes of the horse louse have their counterparts in the pig lice with the same gene content and gene arrangement. However, the gene content and arrangement of the other three minichromosomes of the horse louse, including R-nad4L-rrnS-C, are different from that of the other three minichromosomes of the pig lice.

CONCLUSIONS

Comparison between the horse louse and the pig lice revealed variation in the composition of mitochondrial minichromosomes within the genus Haematopinus, which can be accounted for by gene translocation events between minichromosomes. The current study indicates that inter-minichromosome recombination plays a major role in generating the variation in the composition of mitochondrial minichromosomes and provides novel insights into the evolution of fragmented mitochondrial genomes in the blood-sucking lice.