Sequencing of complete mitochondrial genomes confirms synonymization of Hyalomma asiaticum asiaticum and kozlovi, and advances phylogenetic hypotheses for the Ixodidae

Unraveling the phylogenetics of genetically closely related species, Haemaphysalis japonica and Haemaphysalis megaspinosa, using entire tick mitogenomes and microbiomes

Ticks have a profound impact on public health. Haemaphysalis is one of the most widespread genera in Asia, including Japan. The taxonomy and genetic differentiation of Haemaphysalis spp. is challenging. For instance, previous studies struggled to distinguish Haemaphysalis japonica and Haemaphysalis megaspinosa due to the dearth of nucleotide sequence polymorphisms in widely used barcoding genes. The classification of H. japonica japonica and its related sub-species Haemaphysalis japonica douglasi or Haemaphysalis jezoensis is also confused due to their high morphological similarity and a lack of molecular data that support the current classification. We used mitogenomes and microbiomes of H. japonica and H. megaspinosa to gain deeper insights into the phylogenetic relationships and genetic divergence between two species. Phylogenetic analyses of concatenated nucleotide sequences of protein-coding genes and ribosomal DNA genes distinguished H. japonica and H. megaspinosa as monophyletic clades, with further subdivision within the H. japonica clade. The 16S rRNA and NAD5 genes were valuable markers for distinguishing H. japonica and H. megaspinosa. Population genetic structure analyses indicated that genetic variation within populations accounted for a large proportion of the total variation compared to variation between populations. Microbiome analyses revealed differences in alpha and beta diversity between H. japonica and H. megaspinosa: H. japonica had the higher diversity. Coxiella sp., a likely endosymbiont, was found in both Haemaphysalis species. The abundance profiles of likely endosymbionts, pathogens, and commensals differed between H. japonica and H. megaspinosa: H. megaspinosa was more diverse.

Comparative analysis of the mitochondrial genome of Dermacentor steini from different regions in China

Ticks are a group of blood-sucking ectoparasites that play an important role in human health and livestock production development as vectors of zoonotic diseases. The phylogenetic tree of single genes cannot accurately reflect the true kinship between species. Based on the complete mitochondrial genome analysis one can help to elucidate the phylogenetic relationships among species. In this study, the complete mitochondrial genome of Dermacentor steini (isolate Longyan) was sequenced and compared with the mitochondrial genes of 3 other Chinese isolates (Nanchang, Jinhua and Yingtan). In Dermacentor steini 4 isolates had identical or similar mitochondrial genome lengths and an overall variation of 0.76% between sequences. All nucleotide compositions showed a distinct AT preference. The most common initiation and stop codons were ATG and TAA, respectively. Fewer base mismatches were found in the tRNA gene of D. steini (isolate Longyan), and the vicinity of the control region and tRNA gene was a hot rearrangement region of the genus Dermacentor. Maximum likelihood trees and Bayesian trees indicate that D. steini is most closely related to Dermacentor auratus. The results enrich the mitochondrial genomic data of species in the genus Dermacentor and provide novel insights for further studies on the phylogeographic classification and molecular evolution of ticks.

Comparative mitogenomics elucidates the population genetic structure of Amblyomma testudinarium in Japan and a closely related Amblyomma species in Myanmar

Ticks are the second most important vector capable of transmitting diseases affecting the health of both humans and animals. Amblyomma testudinarium Koch 1844 (Acari: Ixodidae), is a hard tick species having a wide geographic distribution in Asia. In this study, we analyzed the composition of A. testudinarium whole mitogenomes from various geographical regions in Japan and investigated the population structure, demographic patterns, and phylogeographic relationship with other ixodid species. In addition, we characterized a potentially novel tick species closely related to A. testudinarium from Myanmar. Phylogeographic inference and evolutionary dynamics based on the 15 mitochondrial coding genes supported that A. testudinarium population in Japan is resolved into a star-like haplogroup and suggested a distinct population structure of A. testudinarium from Amami island in Kyushu region. Correlation analysis using Mantel test statistics showed that no significant correlation was observed between the genetic and geographic distances calculated between the A. testudinarium population from different localities in Japan. Finally, demographic analyses, including mismatch analysis and Tajima's D test, suggested a possibility of recent population expansion occurred within Japanese haplogroup after a bottleneck event. Although A. testudinarium has been considered widespread and common in East and Southeast Asia, the current study suggested that potentially several cryptic Amblyomma spp. closely related to A. testudinarium are present in Asia.

Prediction and validation of cross-protective candidate antigen of Hyalomma asiaticum cathepsin L between H. asiaticum and H. anatolicum

Hyalomma asiaticum and H. anatolicum are tick species in Eurasia and Africa with major medical and veterinary significance. Beside their direct pathogenic effects, H. asiaticum and H. anatolicum are vectors of important diseases of livestock and in some instances of zoonoses. In search of ways to address the increasing incidence of global acaricide resistance, tick control through vaccination is regarded as a sustainable alternative approach. Cathepsin L-like cysteine protease (CPL) is a potent hemoglobinase, and plays important roles in the digestion of blood acquired from a host. CPL from H. anatolicum (HanCPL) with high similarity (> 90%) for H. asiaticum CPL (HasCPL) were aligned by in silico analysis. After further in vitro validation, the anti-HasCPL sera have cross-reactivity between the different total native protein of life stages and tissues for H. asiaticum and H. anatolicum. Furthermore, we further confirmed that recombinant HasCPL (rHasCPL) immunized rabbits were partially cross-protected (54.8%) by H. anatolicum infestation.

Characterization of the complete mitochondrial genomes of two Ixodes ticks, I. nipponensis and Ixodes (Pholeoixodes) sp

In this study, the authors sequenced and characterized the complete mitochondrial (mt) genomes of two hard ticks of the genus Ixodes, I. nipponensis and Ixodes (Pholeoixodes) sp., which were 14 505 and 14 543 bp in length, respectively. Their mt genomes encoded 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNA genes and two ribosomal RNA genes, and have only one non-coding region. The gene order in their mt genomes was the same as that of other Ixodes spp. mt genomes. The average sequence identity, combined nucleotide diversity, non-synonymous/synonymous substitutions ratio analyses consistently demonstrated that cox1, rrnS, cox2, cox3 and cytb were the most conserved and atp8, nad6 and nad2 were the most variable genes across Ixodes mitogenomes. Phylogeny of the present Ixodes spp., and other selected hard tick species, based on concatenated amino acid sequences of PCGs, confirmed their position within the genus Ixodes and sub-family Ixodinae. The novel mt markers described herein will be useful for further studies of the population genetics, molecular epidemiology and systematics of hard ticks.

Characterization of the complete mitochondrial genome of Amblyomma ovale, comparative analyses and phylogenetic considerations

In this study, we sequenced two complete mitochondrial genomes of Amblyomma ovale, a tick of public health importance. Sequencing two distinct individuals, the resulting mitochondrial genomes were 14,756 and 14,760 bp in length and maintained the same gene order previously reported in Amblyomma. These were combined with RNA-seq derived mitochondrial sequences from three additional species, Amblyomma aureolatum, Amblyomma maculatum, and Amblyomma moreliae, to carry out mitogenome comparative and evolutionary analyses against all previously published tick mitochondrial genomes. We described a derivative genome rearrangement that isolates Ixodes from the remaining Ixodidae and consists of both a reverse translocation as well as an event of Tandem Duplication Random Loss. Genetic distance analyses indicated that cox2, nd1, nd5, and 16S are good candidates for future population studies in A. ovale. The phylogenetic analyses corroborated the utility of complete mitochondrial genomes as phylogenetic markers within the group. This study further supplements the genome information available for Amblyomma and facilitates future evolutionary and population genetic studies within the genus.

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.

Disentangling the interplay of positive and negative selection forces that shaped mitochondrial genomes of Gammarus pisinnus and Gammarus lacustris

We hypothesized that the mitogenome of Gammarus lacustris (GL), native to the Qinghai-Tibet Plateau, might exhibit genetic adaptations to the extreme environmental conditions associated with high altitudes (greater than 3000 m). To test this, we also sequenced the mitogenome of Gammarus pisinnus (GP), whose native range is close to the Tibetan plateau, but at a much lower altitude (200-1500 m). The two mitogenomes exhibited conserved mitochondrial architecture, but low identity between genes (55% atp8 to 76.1% cox1). Standard (homogeneous) phylogenetic models resolved Gammaridae as paraphyletic, but 'heterogeneous' CAT-GTR model as monophyletic. In indirect support of our working hypothesis, GL, GP and Gammarus fossarum exhibit evidence of episodic diversifying selection within the studied Gammaroidea dataset. The mitogenome of GL generally evolves under a strong purifying selection, whereas GP evolves under directional (especially pronounced in atp8) and/or relaxed selection. This is surprising, as GP does not inhabit a unique ecological niche compared to other gammarids. We propose that this rapid evolution of the GP mitogenome may be a reflection of its relatively recent speciation and heightened non-adaptive (putatively metabolic rate-driven) mutational pressures. To test these hypotheses, we urge sequencing mitogenomes of remaining Gammarus species populating the same geographical range as GP.

Tick mitochondrial genomes: structural characteristics and phylogenetic implications

Ticks are obligate blood-sucking arachnid ectoparasites from the order Acarina, and many are notorious as vectors of a wide variety of zoonotic pathogens. However, the systematics of ticks in several genera is still controversial. The mitochondrial genome (mt-genome) has been widely used in arthropod phylogeny, molecular evolution and population genetics. With the development of sequencing technologies, an increasing number of tick mt-genomes have been sequenced and annotated. To date, 63 complete tick mt-genomes are available in the NCBI database, and these genomes have become an increasingly important genetic resource and source of molecular markers in phylogenetic studies of ticks in recent years. The present review summarizes all available complete mt-genomes of ticks in the NCBI database and analyses their characteristics, including structure, base composition and gene arrangement. Furthermore, a phylogenetic tree was constructed using mitochondrial protein-coding genes (PCGs) and ribosomal RNA (rRNA) genes from ticks. The results will provide important clues for deciphering new tick mt-genomes and establish a foundation for subsequent taxonomic research.

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.