Heteroplasmy in the mitochondrial genomes of human lice and ticks revealed by high throughput sequencing

The mitochondrial genome of Bottapotamon fukienense (Brachiura: Potamidae) is fragmented into two chromosomes

BACKGROUND

China is the hotspot of global freshwater crab diversity, but their wild populations are facing severe pressures associated with anthropogenic factors, necessitating the need to map their taxonomic and genetic diversity and design conservation policies.

RESULTS

Herein, we sequenced the mitochondrial genome of a Chinese freshwater crab species Bottapotamon fukienense, and found that it is fragmented into two chromosomes. We confirmed that fragmentation was not limited to a single specimen or population. Chromosome 1 comprised 15,111 base pairs (bp) and there were 26 genes and one pseudogene (pseudo-nad1) encoded on it. Chromosome 2 comprised 8,173 bp and there were 12 genes and two pseudogenes (pseudo-trnL2 and pseudo-rrnL) encoded on it. Combined, they comprise the largest mitogenome (23,284 bp) among the Potamidae. Bottapotamon was the only genus in the Potamidae dataset exhibiting rearrangements of protein-coding genes. Bottapotamon fukienense exhibited average rates of sequence evolution in the dataset and did not differ in selection pressures from the remaining Potamidae.

CONCLUSIONS

This is the first experimentally confirmed fragmentation of a mitogenome in crustaceans. While the mitogenome of B. fukienense exhibited multiple signs of elevated mitogenomic architecture evolution rates, including the exceptionally large size, duplicated genes, pseudogenisation, rearrangements of protein-coding genes, and fragmentation, there is no evidence that this is matched by elevated sequence evolutionary rates or changes in selection pressures.

Unparalleled mitochondrial heteroplasmy and Wolbachia co-infection in the non-model bee, Amphylaeus morosus

Mitochondrial heteroplasmy is the occurrence of more than one type of mitochondrial DNA within a single individual. Although generally reported to occur in a small subset of individuals within a species, there are some instances of widespread heteroplasmy across entire populations. Amphylaeus morosus is an Australian native bee species in the diverse and cosmopolitan bee family Colletidae. This species has an extensive geographical range along the eastern Australian coast, from southern Queensland to western Victoria, covering approximately 2,000 km. Seventy individuals were collected from five localities across this geographical range and sequenced using Sanger sequencing for the mitochondrial cytochrome c oxidase subunit I (COI) gene. These data indicate that every individual had the same consistent heteroplasmic sites but no other nucleotide variation, suggesting two conserved and widespread heteroplasmic mitogenomes. Ion Torrent shotgun sequencing revealed that heteroplasmy occurred across multiple mitochondrial protein-coding genes and is unlikely explained by transposition of mitochondrial genes into the nuclear genome (NUMTs). DNA sequence data also demonstrated a consistent co-infection of Wolbachia across the A. morosus distribution with every individual infected with both bacterial strains. Our data are consistent with the presence of two mitogenomes within all individuals examined in this species and suggest a major divergence from standard patterns of mitochondrial inheritance. Because the host's mitogenome and the Wolbachia genome are genetically linked through maternal inheritance, we propose three possible hypotheses that could explain maintenance of the widespread and conserved co-occurring bacterial and mitochondrial genomes in this species.

Comparative hologenomics of two Ixodes scapularis tick populations in New Jersey

Tick-borne diseases, such as those transmitted by the blacklegged tick Ixodes scapularis, are a significant and growing public health problem in the US. There is mounting evidence that co-occurring non-pathogenic microbes can also impact tick-borne disease transmission. Shotgun metagenome sequencing enables sampling of the complete tick hologenome—the collective genomes of the tick and all of the microbial species contained therein, whether pathogenic, commensal or symbiotic. This approach simultaneously uncovers taxonomic composition and allows the detection of intraspecific genetic variation, making it a useful tool to compare spatial differences across tick populations. We evaluated this approach by comparing hologenome data from two tick samples (N = 6 ticks per location) collected at a relatively fine spatial scale, approximately 23 km apart, within a single US county. Several intriguing variants in the data between the two sites were detected, including polymorphisms in both in the tick’s own mitochondrial DNA and that of a rickettsial endosymbiont. The two samples were broadly similar in terms of the microbial species present, including multiple known tick-borne pathogens (Borrelia burgdorferi, Babesia microti, and Anaplasma phagocytophilum), filarial nematodes, and Wolbachia and Babesia species. We assembled the complete genome of the rickettsial endosymbiont (most likely Rickettsia buchneri) from both populations. Our results provide further evidence for the use of shotgun metagenome sequencing as a tool to compare tick hologenomes and differentiate tick populations across localized spatial scales.

Genetic variation and heteroplasmy of Varroa destructor inferred from ND4 mtDNA sequences

Varroa destructor, a parasitic mite of the western honey bee, Apis mellifera L., is a serious threat to colonies and beekeeping worldwide. Population genetics studies of the mite have provided information on two mitochondrial haplotypes infecting honey bee colonies, named K and J (after Korea and Japan, respectively, where they were originally identified). On the American continent, the K haplotype is much more prevalent, with the J haplotype only detected in some areas of Brazil. The aims of the present study were to assess the genetic diversity of V. destructor populations in the major beekeeping region of Argentina and to evaluate the presence of heteroplasmy at the nucleotide level. Phoretic mites were collected from managed A. mellifera colonies in ten localities, and four mitochondrial DNA (mtDNA) regions (COXI, ND4, ND4L, and ND5) were analyzed. Based on cytochrome oxidase subunit I (COXI) sequencing, exclusively the K haplotype of V. destructor was detected. Furthermore, two sub-haplotypes (KArg-N1 and KArg-N2) were identified from a variation in ND4 sequences and the frequency of these sub-haplotypes was found to significantly correlate with geographical latitude. The occurrence of site heteroplasmy was also evident for this gene. Therefore, ND4 appears to be a sensitive marker for detecting genetic variability in mite populations. Site heteroplasmy emerges as a phenomenon that could be relatively frequent in V. destructor.

Investigating the genetic diversity and presence of forensically informative nucleotide sequences in Indian antelope (Antilope cervicapra) using multiple genes of the mitochondrial genome

Indian antelope or Blackbuck (Antilope cervicapra) is one of the widely distributed endemic species in India among wild bovids and a majority of preferred habitats are in human-dominated landscapes. Poaching threats and habitat degradation are major factors for the decline in Blackbuck population from its distribution range. Till date, there is no detailed study using molecular techniques in India on Blackbuck, except a few studies entailing phylogenetic scenario based on inadequate sampling and DNA sequences restricted over limited geographic areas. In view of this, the present study is aimed to screen the Blackbuck samples from a large part of its distribution range and to investigate the genetic diversity as well as to identify the forensically informative nucleotide sequences (FINS) for species identification. We relied on multi-genes approach using three genes of mtDNA genome viz. Cytochrome Oxidase I, Cytochrome b and 16S rRNA and identified the FINS in the Blackbuck population along with conspecific sequences divergence and genetic diversity indices. In all three genes, we observed 8 to 17 haplotypes with the intra-species sequence divergence of 0.004-0.016. Inter-species sequence divergence with the other closely related species of the Blackbuck was 0.0225-0.033. We report the presence of FINS across three genes from 12 to 18 and found more informative nucleotide sites using Cytochrome Oxidase I genes compared to Cytochrome b and 16S rRNA gene. We did not observe the presence of geographic-specific FINS amongst Blackbuck population that can be used to assign individuals to geographic origin. Besides, in the phylogenetic tree, samples from different locations did not cluster into geographic-specific clade and exhibited mixed homology for these sequences. We suggest exploring the feasibility of using nuclear markers for population assignment.

Homoplasy or plesiomorphy? Reconstruction of the evolutionary history of mitochondrial gene order rearrangements in the subphylum Neodermata

Recent mitogenomic studies have exposed a gene order (GO) shared by two classes, four orders and 31 species ('common GO') within the flatworm subphylum Neodermata. There are two possible hypotheses for this phenomenon: convergent evolution (homoplasy) or shared ancestry (plesiomorphy). To test those, we conducted a meta-analysis on all available mitogenomes to infer the evolutionary history of GO in Neodermata. To improve the resolution, we added a newly sequenced mitogenome that exhibited the common GO, Euryhaliotrema johni (Ancyrocephalinae), to the dataset. Phylogenetic analyses conducted on two datasets (nucleotides of all 36 genes and amino acid sequences of 12 protein coding genes) and four algorithms (MrBayes, RAxML, IQ-TREE and PhyloBayes) produced topology instability towards the tips, so ancestral GO reconstructions were conducted using TreeREx and MLGO programs using all eight obtained topologies, plus three unique topologies from previous studies. The results consistently supported the second hypothesis, resolving the common GO as a plesiomorphic ancestral GO for Neodermata, Cestoda, Monopisthocotylea, Cestoda + Trematoda and Cestoda + Trematoda + Monopisthocotylea. This allowed us to trace the evolutionary GO scenarios from each common ancestor to its descendants amongst the Monogenea and Cestoda classes, and propose that the common GO was most likely retained throughout all of the common ancestors, leading to the extant species possessing the common GO. Neodermatan phylogeny inferred from GOs was largely incongruent with all 11 topologies described above, but it did support the mitogenomic dataset in resolving Polyopisthocotylea as the earliest neodermatan branch. Although highly derived GOs might be of some use in resolving isolated taxonomic and phylogenetic uncertainties, we conclude that, due to the discontinuous nature of their evolution, they tend to produce artefactual phylogenetic relationships, which makes them unsuitable for phylogenetic reconstruction in Neodermata. Wider and denser sampling of neodermatan mitogenomic sequences will be needed to infer the evolutionary pathways leading to the observed diversity of GOs with confidence.

The Complete Mitochondrial Genome and Expression Profile of Mitochondrial Protein-Coding Genes in the Bisexual and Parthenogenetic Haemaphysalis longicornis

The tick Haemaphysalis longicornis is widely distributed in eastern Asia, New Zealand and Australia, and is well-known as a vector of multiple zoonotic pathogens. This species exhibits two reproductive strategies, bisexual and obligate parthenogenetic reproduction. Hence, in the current study, the complete mitochondrial genomes of the bisexual and parthenogenetic populations were assembled and analyzed, and the expression of the mitochondrial protein-coding genes was evaluated and compared between the two reproductive populations. The results indicated that the length of the mitochondrial genomes of the two reproductive populations is 14,694 and 14,693 bp in the bisexual and parthenogenetic populations, respectively. The AT content in the mitochondrial genome of the bisexual and obligate parthenogenetic population reached 77.22 and 77.34%, respectively. The phylogenetic tree was constructed combining 13 protein-coding genes, which showed that the genetic distance between the bisexual and parthenogenetic populations was less than that between the subspecies. The expression of the mitochondrial protein-coding genes was quantitatively analyzed at different feeding status for the bisexual and parthenogenetic populations, and the results showed significant differences in the expression patterns of these genes, suggesting that they might trigger specific energy utilization mechanisms due to their different reproductive strategies and environmental pressures.

Paternal leakage and mtDNA heteroplasmy in Rhipicephalus spp. ticks

Paternal leakage of mitochondrial DNA (mtDNA) and heteroplasmy have been recently described in several animal species. In arthropods, by searching in the Scopus database, we found only 23 documented cases of paternal leakage. Therefore, although arthropods represent a large fraction of animal biodiversity, this phenomenon has been investigated only in a paucity of species in this phylum, thus preventing a reliable estimate of its frequency. Here, we investigated the occurrence of paternal leakage and mtDNA heteroplasmy in ticks belonging to one of the most significant tick species complexes, the so-called Rhipicephalus sanguineussensu lato. By developing a multiplex allele-specific PCR assay targeting a fragment of the 12S rRNA ribosomal region of the mtDNA, we showed the occurrence of paternal leakage and mtDNA heteroplasmy in R. sanguineuss.l. ticks originated from experimental crosses, as well as in individuals collected from the field. Our results add a new evidence of paternal leakage in arthropods and document for the first time this phenomenon in ticks. Furthermore, they suggest the importance of using allele-specific assays when searching for paternal leakage and/or heteroplasmy, as standard sequencing methods may fail to detect the rare mtDNA molecules.

Mitochondrial genomes of two diplectanids (Platyhelminthes: Monogenea) expose paraphyly of the order Dactylogyridea and extensive tRNA gene rearrangements

BACKGROUND

Recent mitochondrial phylogenomics studies have reported a sister-group relationship of the orders Capsalidea and Dactylogyridea, which is inconsistent with previous morphology- and molecular-based phylogenies. As Dactylogyridea mitochondrial genomes (mitogenomes) are currently represented by only one family, to improve the phylogenetic resolution, we sequenced and characterized two dactylogyridean parasites, Lamellodiscus spari and Lepidotrema longipenis, belonging to a non-represented family Diplectanidae.

RESULTS

The L. longipenis mitogenome (15,433 bp) contains the standard 36 flatworm mitochondrial genes (atp8 is absent), whereas we failed to detect trnS1, trnC and trnG in L. spari (14,614 bp). Both mitogenomes exhibit unique gene orders (among the Monogenea), with a number of tRNA rearrangements. Both long non-coding regions contain a number of different (partially overlapping) repeat sequences. Intriguingly, these include putative tRNA pseudogenes in a tandem array (17 trnV pseudogenes in L. longipenis, 13 trnY pseudogenes in L. spari). Combined nucleotide diversity, non-synonymous/synonymous substitutions ratio and average sequence identity analyses consistently showed that nad2, nad5 and nad4 were the most variable PCGs, whereas cox1, cox2 and cytb were the most conserved. Phylogenomic analysis showed that the newly sequenced species of the family Diplectanidae formed a sister-group with the Dactylogyridae + Capsalidae clade. Thus Dactylogyridea (represented by the Diplectanidae and Dactylogyridae) was rendered paraphyletic (with high statistical support) by the nested Capsalidea (represented by the Capsalidae) clade.

CONCLUSIONS

Our results show that nad2, nad5 and nad4 (fast-evolving) would be better candidates than cox1 (slow-evolving) for species identification and population genetics studies in the Diplectanidae. The unique gene order pattern further suggests discontinuous evolution of mitogenomic gene order arrangement in the Class Monogenea. This first report of paraphyly of the Dactylogyridea highlights the need to generate more molecular data for monogenean parasites, in order to be able to clarify their relationships using large datasets, as single-gene markers appear to provide a phylogenetic resolution which is too low for the task.

Species Concepts: What about Ticks?

Since ancient times, philosophers and taxonomists have tried to classify forms of life. This is what taxonomy is about: the science of identifying, naming, classifying, and describing organisms. In this article I address the issue of the species concept in tick taxonomy. While the typological species concept is still the most widely used, the biological and phylogenetic species concepts are growing in popularity among tick taxonomists. The integrative approach is increasingly being used, but the question is how to define a tick species when using this approach, particularly if data are incongruent. The adoption of an integrative species concept is discussed, in light of recent advances in our understanding of the genetics, morphology, and biology of ticks.

Haplotype identification and detection of mitochondrial DNA heteroplasmy in Varroa destructor mites using ARMS and PCR-RFLP methods

In the present study, amplification refractory mutation system (ARMS) and polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) methods were used for identification of recently described Serbia 1 (S1) and Peshter 1 (P1) mitochondrial haplotypes of Varroa destructor. Based on single nucleotide polymorphisms (SNPs) within cytochrome oxidase 1 (cox1) and cytochrome b (cytb) gene sequences, a total of 64 adult V. destructor females were analyzed from locations where the S1 and P1 haplotypes had been detected previously. Results of haplotype identification obtained by ARMS and PCR-RFLP methods were completely consistent with the sequencing data. Furthermore, in some analyzed samples the occurrence of site heteroplasmy at haplotype-defining sites was detected, as it was confirmed by double peaks in the sequence chromatograms. Neither mites with simultaneous nucleotide variability, nor those with combined SNP and heteroplasmy in cox1 and cytb were found. Given that this is the first occurrence of site heteroplasmy in V. destructor, the origin of this phenomenon and possible specific traits of heteroplasmic mites have yet to be determined.

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.

Comparative population genetics of two invading ticks: Evidence of the ecological mechanisms underlying tick range expansions

Two species of ixodid tick, Ixodes affinis Neumann and Amblyomma maculatum Koch, are simultaneously expanding their ranges throughout the mid-Atlantic region of the US. Although we have some understanding of the ecology and life history of these species, the ecological mechanisms governing where and how new populations establish and persist are unclear. To assess population connectivity and ancestry, we sequenced a fragment of the 16S mitochondrial rRNA gene from a representative sample of individuals of both species from populations throughout the eastern US. We found that despite overlapping host preferences throughout ontogeny, each species exhibited very different genetic and geographic patterns of population establishment and connectivity. I. affinis was of two distinct mitochondrial clades, with a clear geographic break separating northern and southern populations. Both I. affinis populations showed evidence of recent expansion, although the southern population was more genetically diverse, indicating a longer history of establishment. A. maculatum exhibited diverse haplotypes that showed no significant relationship with geographic patterns and little apparent connectivity between sites. Heteroplasmy was also observed in the 16S mitochondrial rRNA gene in 3.5% of A. maculatum individuals. Genetic evidence suggests that these species rely on different key life stages to successfully disperse into novel environments, and that host vagility, habitat stability and habitat connectivity all play critical roles in the establishment of new tick populations.

Population and demographic structure of Ixodes scapularis Say in the eastern United States

Introduction

The most significant vector of tick-borne pathogens in the United States is Ixodes scapularis Say (the blacklegged tick). Previous studies have identified significant genetic, behavioral and morphological differences between northern vs. southern populations of this tick. Because tick-borne pathogens are dependent on their vectors for transmission, a baseline understanding of the vector population structure is crucial to determining the risks and epidemiology of pathogen transmission.

Methods

We investigated population genetic variation of I. scapularis populations in the eastern United States using a multilocus approach. We sequenced and analyzed the mitochondrial COI and 16S genes and three nuclear genes (serpin2, ixoderin B and lysozyme) from wild specimens.

Results

We identified a deep divergence (3–7%) in I. scapularis COI gene sequences from some southern specimens, suggesting we had sampled a different Ixodes species. Analysis of mitochondrial 16S rRNA sequences did not support this hypothesis and indicated that all specimens were I. scapularis. Phylogenetic analysis and analysis of molecular variance (AMOVA) supported significant differences between northern vs. southern populations. Demographic analysis suggested that northern populations had experienced a bottleneck/expansion event sometime in the past, possibly associated with Pleistocene glaciation events.

Conclusions

Similar to other studies, our data support the division of northern vs. southern I. scapularis genetic lineages, likely due to differences in the demographic histories between these geographic regions. The deep divergence identified in some COI gene sequences highlights a potential hazard of relying solely on COI for species identification (“barcoding”) and population genetics in this important vector arthropod.