Tick-box for 3'-end formation of mitochondrial transcripts in Ixodida, basal chelicerates and Drosophila

Quantification of rearrangements and evolution of mitochondrial gene order of Acari (Chelicerata: Arachnida)

Invertebrate mitogenomes are generally fixed with formal 37 genes: 13 PCGs encoded subunits of OXPHOS, 2 ribosomal RNA (rRNA) functional in the translation of these PCGs and 22 transfer RNA (tRNA) genes. The order of these genes varies greatly among organisms and named rearrangement. Rearrangement patterns of mitochondrial genomes may shed light on mutation processes and evolutionary relationships of organisms. Mitochondrial gene organization is highly variable among Acari, so rearrangement is a very common mitogenomic pattern in this group. In this study, 258 unique Acari (Acariformes + Parasitiformes) mitogenomes were downloaded from NCBI and studied about rearrangement patterns. Sixty-seven mitotypes were determined among Acari and the most rearranged genes were trnL1 and nad2. Following that, trnI, trnS1, trnN, trnE, trnT, and trnP genes are remarkably mobile (RF > 95%). Conversely, atp6, cox3, trnG, and cytb genes also appears to be quite stable (RF < 20%). Within Acari, mean distance calculations are varied from 1.210 in atp8 to 0.155 in rrnS. Contrary to expectations, among Acari mobile tRNA genes appear to be conserved in sequences, whereas PCGs have higher distance values and seem to be mutated. Consistently, tRNA genes seem saturated, but some PCGs (atp6, cox genes, cytb, nad1, and nad6) are not saturated. These values do not correlate with each other (p > 0.005). This discrepancy may indicate that the genes were rearranged after mutation load; consistent with this, DAMBE saturation values are also not correlated with RF values. Parasitiformes mitogenomes are more mobile than Acariformes mitogenomes and may be under the effect of selective sweeping.

Characterization of the complete mitochondrial genome and phylogenetic analyses of Haemaphysalis tibetensis Hoogstraal, 1965 (Acari: Ixodidae)

Ticks are specialized ectoparasites that feed on blood, causing physical harm to the host and facilitating pathogen transmission. The genus Haemaphysalis contains vectors for numerous infectious agents. These agents cause various diseases in humans and animals. Mitochondrial genome sequences serve as reliable molecular markers, forming a crucial basis for evolutionary analyses, studying species origins, and exploring molecular phylogeny. We extracted mitochondrial genome from the enriched mitochondria of Haemaphysalis tibetensis and obtained a 14,714-bp sequence. The mitochondrial genome consists of 13 protein-coding genes (PCGs), two ribosomal RNA, 22 transfer RNAs (tRNAs), and two control regions. The nucleotide composition of H. tibetensis mitochondrial genome was 38.38 % for A, 9.61 % for G, 39.32 % for T, and 12.69 % for C. The A + T content of H. tibetensis mitochondrial genome was 77.7 %, significantly higher than the G + C content. The repeat units of H. tibetensis exhibited two identical repeat units of 33 bp in length, positioned downstream of nad1 and rrnL genes. Furthermore, phylogenetic analyses based on the 13 PCGs indicated that Haemaphysalis tibetensis (subgenus Allophysalis) formed a monophyletic clade with Haemaphysalis nepalensis (subgenus Herpetobia) and Haemaphysalis danieli (subgenus Allophysalis). Although the species Haemaphysalis inermis, Haemaphysalis kitaokai, Haemaphysalis kolonini, and Haemaphysalis colasbelcouri belong to the subgenus Alloceraea, which were morphologically primitive hemaphysalines just like H. tibetensis, these four tick species cannot form a single clade with H. tibetensis. In this study, the whole mitochondrial genome sequence of H. tibetensis from Tibet was obtained, which enriched the mitochondrial genome data of ticks and provided genetic markers to study the population heredity and molecular evolution of the genus Haemaphysalis.

Seventy-eight entire mitochondrial genomes and nuclear rRNA genes provide insight into the phylogeny of the hard ticks, particularly the Haemaphysalis species, Africaniella transversale and Robertsicus elaphensis

Hoogstraal and Kim (1985) proposed from morphology, three groups of Haemaphysalis subgenera: (i) the "structurally advanced"; (ii) the "structurally intermediate"; and (iii) the "structurally primitive" subgenera. Nuclear gene phylogenies, however, did not indicate monophyly of these morphological groups but alas, only two mitochondrial (mt) genomes from the "structurally intermediate" subgenera had been sequenced. The phylogeny of Haemaphysalis has not yet been resolved. We aimed to resolve the phylogeny of the genus Haemaphysalis, with respect to the subgenus Alloceraea. We presented 15 newly sequenced and annotated mt genomes from 15 species of ticks, five species of which have not been sequenced before, and four new 18S rRNA and 28S rRNA nuclear gene sequences. Our datasets were constructed from 10 mt protein-coding genes, cox1, and the 18S and 28S nuclear rRNA genes. We found a 132-bp insertion between tRNA-Glu (E) gene and the nad1 gene in the mt genome of Haemaphysalis (Alloceraea) inermis that resembles insertions in H. (Alloceraea) kitaokai and Rhipicephalus (Boophilus) geigyi. Our mt phylogenies had the three species of Amblyomma (Aponomma) we sequenced embedded in the main clade of Amblyomma: Am. (Aponomma) fimbriatum, Am. (Aponomma) gervaisi and Am. (Aponomma) latum. This is further support for the hypothesis that the evolution of eyes appears to have occurred in the most-recent-common-ancestor of Amblyocephalus (i.e. Amblyomminae plus Rhipicephalinae) and that eyes were subsequently lost in the most-recent-common-ancestor of the subgenus Am. (Aponomma). Either Africaniella transversale or Robertsicus elaphensis, or perhaps Af. transversale plus Ro. elaphensis, appear to be the sister-group to the rest of the metastriate Ixodida. Our cox1 phylogenies did not indicate monophyly of the "structurally primitive", "structurally intermediate" nor the "structurally advanced" groups of Haemaphysalis subgenera. Indeed, the subgenus Alloceraea may be the only monophyletic subgenus of the genus Haemaphysalis sequenced thus far. All of our mt genome and cox1 phylogenies had the subgenus Alloceraea in a clade that was separate from the rest of the Haemaphysalis ticks. If Alloceraea is indeed the sister to the rest of the Haemaphysalis subgenera this would resonate with the argument of Hoogstraal and Kim (1985), that Alloceraea was a subgenus of "primitive" Haemaphysalis. Alectorobius capensis from Japan had a higher genetic-identity to A. sawaii, which was also from Japan, than to the A. capensis from South Africa. This indicates that A. capensis from Japan may be a cryptic species with respect to the A. capensis from South Africa.

Complete Mitogenomes of Ticks Ixodes acutitarsus and Ixodes ovatus Parasitizing Giant Panda: Deep Insights into the Comparative Mitogenomic and Phylogenetic Relationship of Ixodidae Species

Ticks rank second in the world as vectors of disease. Tick infestation is one of the factors threatening the health and survival of giant pandas. Here, we describe the mitogenomes of Ixodes acutitarsus and Ixodes ovatus parasitizing giant pandas, and perform comparative and phylogenetic genomic analyses on the newly sequenced and other available mitogenomes of hard ticks. All six newly determined mitogenomes contain a typical gene component and share an ancient Arthropoda gene arrangement pattern. Our study suggests that I. ovatus is a species complex with high genetic divergence, indicating that different clades of I. ovatus represent distinct species. Comparative mitogenomic analyses show that the average A + T content of Ixodidae mitogenomes is 78.08%, their GC-skews are strongly negative, while AT-skews fluctuate around 0. A large number of microsatellites are detected in Ixodidae mitogenomes, and the main microsatellite motifs are mononucleotide A and trinucleotide AAT. We summarize five gene arrangement types, and identify the trnY-COX1-trnS1-COX2-trnK-ATP8-ATP6-COX3-trnG fragment is the most conserved region, whereas the region near the control region is the rearrangement hotspot in Ixodidae mitogenomes. The phylogenetic trees based on 15 genes provide a very convincing relationship (Ixodes + (Robertsicus + ((Bothriocroton + Haemaphysalis) + (Amblyomma + (Dermacentor + (Rhipicentor + (Hyalomma + Rhipicephalus))))))) with very strong supports. Remarkably, Archaeocroton sphenodonti is embedded in the Haemaphysalis clade with strong supports, resulting in paraphyly of the Haemaphysalis genus, so in-depth morphological and molecular studies are essential to determine the taxonomic status of A. sphenodonti and its closely related species. Our results provide new insights into the molecular phylogeny and evolution of hard ticks, as well as basic data for population genetics assessment and efficient surveillance and control for the giant panda-infesting ticks.

Phylogeographic dynamics of the arthropod vector, the blacklegged tick (Ixodes scapularis)

BACKGROUND

The emergence of vector-borne pathogens in novel geographic areas is regulated by the migration of their arthropod vectors. Blacklegged ticks (Ixodes scapularis) and the pathogens they vector, including the causative agents of Lyme disease, babesiosis and anaplasmosis, continue to grow in their population sizes and to expand in geographic range. Migration of this vector over the previous decades has been implicated as the cause of the re-emergence of the most prevalent infectious diseases in North America.

METHODS

We systematically collected ticks from across New York State (hereafter referred to as New York) from 2004 to 2017 as part of routine tick-borne pathogen surveillance in the state. This time frame corresponds with an increase in range and incidence of tick-borne diseases within New York. We randomly sampled ticks from this collection to explore the evolutionary history and population dynamics of I. scapularis. We sequenced the mitochondrial genomes of each tick to characterize their current and historical spatial genetic structure and population growth using phylogeographic methods.

RESULTS

We sequenced whole mitochondrial genomes from 277 ticks collected across New York between 2004 and 2017. We found evidence of population genetic structure at a broad geographic scale due to differences in the relative abundance, but not the composition, of haplotypes among sampled ticks. Ticks were often most closely related to ticks from the same and nearby collection sites. The data indicate that both short- and long-range migration events shape the population dynamics of blacklegged ticks in New York.

CONCLUSIONS

We detailed the population dynamics of the blacklegged tick (Ixodes scapularis) in New York during a time frame in which tick-borne diseases were increasing in range and incidence. Migration of ticks occurred at both coarse and fine scales in the recent past despite evidence of limits to gene flow. Past and current tick population dynamics have implications for further range expansion as habitat suitability for ticks changes due to global climate change. Analyses of mitochondrial genome sequencing data will expound upon previously identified drivers of tick presence and abundance as well as identify additional drivers. These data provide a foundation on which to generate testable hypotheses on the drivers of tick population dynamics occurring at finer scales.

The complete mitochondrial genome of Ixodes vespertilionis (Acari: Ixodidae)

Ixodes vespertilionis is a tick parasitizing on the bodies of bats. In our study, the complete mitogenome of I. vespertilionis was determined by using Illumina sequencing technology. The mitogenome was 14,559 bp in size and was predicted to encode 37 genes including 13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes, and one control region. The gene order of the mitogenome is identical to Argasidae and non-Australasian Prostriata. The phylogenetic analysis by the Maximum-likelihood method reveals that I. vespertilionis is phylogenetically closest to Ixodes simplex. These data provide novel reference for further studies on the population genetics and phylogenetics of ticks.

Complete Mitochondrial Genome Sequence of the Gulf Coast Tick (Amblyomma maculatum)

The complete circularized mitochondrial genome sequence of Amblyomma maculatum is 14,803 bp long. It encodes 13 protein coding genes, 2 rRNA genes, 22 tRNA genes, 2 tick box motifs, and 2 control regions. The gene arrangement and content are consistent with those of previously reported Metastriata tick mitochondrial genomes.

Complete mitochondrial genome of Asian longhorned tick, Haemaphysalis longicornis, Neumann, 1901 (Acari: Ixodida: Ixodidae) identified in the United States

Haemaphysalis longicornis (Ixodida: Ixodidae), the Asian longhorned tick, which is native to temperate East Asia, has been recently detected in the northeastern region of the United States, drawing concerns about its potential impact on the US animal and public health sectors. Knowledge about the genetic features of H. longicornis found in the US is limited. Therefore, we sequenced the complete mitochondrial genome (mt-genome) from two H. longicornis ticks recently collected in the State of New York, USA, in 2020. These ticks were morphologically identified and tested for tick-borne pathogens at the Connecticut Veterinary Medical Diagnostic Laboratory (Storrs, CT). The mt-genome was 14,694 bp in length and encoded 37 genes, including 13 protein-coding genes, 22 transfer RNAs, and two ribosomal RNAs. Phylogenetic analysis showed that the mt-genome clustered with those of other H. longicornis identified in China. The mt-genome sequence was 99.7% identical to a H. longicornis mt-genome (GenBank: MK439888) collected in China. The cox1 gene haplotype in these ticks belonged to the H1 type, which is the dominant haplotype present in central NJ and Staten Island, NY. The complete mt-genome data are needed to provide insights into genetic changes and phylogenetic studies of H. longicornis ticks.

Mitochondrial genomes of two Polydora (Spionidae) species provide further evidence that mitochondrial architecture in the Sedentaria (Annelida) is not conserved

Contrary to the early evidence, which indicated that the mitochondrial architecture in one of the two major annelida clades, Sedentaria, is relatively conserved, a handful of relatively recent studies found evidence that some species exhibit elevated rates of mitochondrial architecture evolution. We sequenced complete mitogenomes belonging to two congeneric shell-boring Spionidae species that cause considerable economic losses in the commercial marine mollusk aquaculture: Polydora brevipalpa and Polydora websteri. The two mitogenomes exhibited very similar architecture. In comparison to other sedentarians, they exhibited some standard features, including all genes encoded on the same strand, uncommon but not unique duplicated trnM gene, as well as a number of unique features. Their comparatively large size (17,673 bp) can be attributed to four non-coding regions larger than 500 bp. We identified an unusually large (putative) overlap of 14 bases between nad2 and cox1 genes in both species. Importantly, the two species exhibited completely rearranged gene orders in comparison to all other available mitogenomes. Along with Serpulidae and Sabellidae, Polydora is the third identified sedentarian lineage that exhibits disproportionally elevated rates of mitogenomic architecture rearrangements. Selection analyses indicate that these three lineages also exhibited relaxed purifying selection pressures.

Complete mitochondrial genome characterization and phylogenetic analyses of the main vector of Crimean-Congo haemorrhagic fever virus: Hyalomma marginatum Koch, 1844

The Mediterranean tick, Hyalomma marginatum, is the most important vector of Crimean-Congo haemorrhagic fever virus and several pathogens that cause animal and human diseases and economic losses to livestock production. Given the medical and veterinary importance of this tick species, we sequenced and characterized its mitochondrial genome (mitogenome) for the first time. We designed two new primer sets and combined long-range PCR with next generation sequencing to generate complete mitogenomes with deep coverage from 10 H. marginatum adults. The mitogenomes contained 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal subunits, two control regions, and three tick-box motifs. The nucleotide composition of the H. marginatum mitogenomes were A+T biased (79.76%) and exhibited negative AT- and GC- skews across most PCGs. All PCGs were initiated by ATK codons and two truncated termination codons were seen in the COX2 and COX3 genes. All tRNAs exhibited typical cloverleaf structures, except for tRNA^Cys and tRNA^Ser1. A total of 62 polymorphic sites defined ten unique haplotypes. Phylogenetic analyses based on the 13 PCGs of 56 tick species revealed that four Hyalomma species (H. marginatum, H. asiaticum, H. rufipes, and H. truncatum) formed a monophyletic clade with strong support. The results of this study provide a comprehensive resource for further studies on the systematics, population genetics, molecular epidemiology, and evolution of ticks.

The complete mitochondrial genome of Dermacentor (Indocentor) auratus (Acari, Ixodidae)

Dermacentor (Indocentor) auratus Supino, 1897 is a prominent ixodid vector of numerous pathogens of public health and veterinary importance. Using long-range PCR of two overlapping regions sequenced on an Illumina MiSeq machine, the complete mitochondrial genome of D. auratus is reported here. The resulting contigs were able to be assembled into a complete and circularised genome which had the general organisation of the mitochondrial genomes of the Metastriates. It had a total length of 14,766 bp and contained 37 genes, including 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes, as well as 2 non-coding control regions and 3 tick-boxes. The phylogenetic analysis on the whole mitogenome confirmed the position of D. auratus within the Dermacentor clade.

Phylogenies from mitochondrial genomes of 120 species of ticks: Insights into the evolution of the families of ticks and of the genus Amblyomma

The evolution and phylogenetic relationships of the ticks at both the family and genus levels are contested. The genus Amblyomma and its subgenera are in a state of flux; moreover, the relationships among the three tick families are controversial due to conflicting phylogenetic support for different arrangements of the three families of living ticks. With 18 newly sequenced mitochondrial (mt) genomes of ticks included, we executed the largest mt genome phylogenetic study of ticks so far. Phylogenetic trees were inferred from one sea spider mt genome, one horseshoe crab, five mite mt genomes and 146 tick mt genomes from 120 species: 153 mt genomes in total. Sixteen phylogenetic trees were inferred from 10 datasets using both maximum likelihood and Bayesian inference methods. We describe the first novel mt gene-arrangement for the metastriate Ixodidae in Amblyomma (Africaniella) transversale. Also, three unusual partial 16S rRNA gene inserts were found in the mt genome of Haemaphysalis (Alloceraea) kitaokai: we consider the possible role of past genome translocation events in the formation of these inserts. Our phylogenies revealed evidence that: (i) the genus Amblyomma is polyphyletic with respect to Amblyomma (Africaniella) transversale; (ii) the subgenus Aponomma is apparently embedded in the genus Amblyomma; (iii) Haemaphysalis (Segalia) parva and Haemaphysalis (Alloceraea) kitaokai form a clade to the exclusion of other Haemaphysalis species; and (iv) the phylogenetic position of the family Nuttalliellidae is unstable among phylogenies from different datasets.

Molecular evidence confirms occurrence of Rhipicephalus microplus Clade A in Kenya and sub-Saharan Africa

Background

The tick vector Rhipicephalus microplus which transmits Babesia spp. and rickettsial pathogens has not been reported in Kenya since 1998. More recently, the pathogenic Babesia bovis has been detected in cattle blood DNA. The status of R. microplus in Kenya remains unknown. This study employed morphological and molecular tools to characterize R. microplus originating from Kenya and assess the genetic relationships between Kenyan and other African R. microplus genotypes.

Methods

Ticks were collected in south-eastern Kenya (Kwale County) from cattle and characterized to investigate the existence of R. microplus. Genetic and phylogenetic relationships between the Kenyan and other annotated R. microplus reference sequences was investigated by analysis of the cytochrome c oxidase subunit 1 (cox1) gene. To further characterize Kenyan ticks, we generated low coverage whole genome sequences of two R. microplus, one R. decoloratus and R. appendiculatus. A B. bovis specific TaqMan probe qPCR assay was used to detect B. bovis in gDNA from R. microplus ticks.

Results

Occurrence of R. microplus was confirmed in Kwale County, Kenya. The Kenyan R. microplus cox1 sequences showed very high pairwise identities (> 99%) and clustered very closely with reference African R. microplus sequences. We found a low genetic variation and lack of geographical sub-structuring among the African cox1 sequences of R. microplus. Four complete mitochondrial (mt) genomes for two R. microplus, one R. decoloratus and one R. appendiculatus were assembled from next generation sequence data. The mitochondrial genome sequences of the two Kenyan R. microplus ticks clustered closely with reference genome sequences from Brazil, USA, Cambodia and India forming R. microplus Clade A. No B. bovis was detected in the Kwale R. microplus DNA.

Conclusions

These findings confirm the presence of R. microplus in Kenya and suggest that R. microplus Clade A is prevalent in cattle in sub-Saharan Africa. These and other recent findings of widespread occurrence of R. microplus in Africa provide a strong justification for urgent surveillance to determine and monitor the spread of R. microplus and vector competence of Boophilus ticks for B. bovis in Africa, with the ultimate goal of strategic control.

Precise annotation of tick mitochondrial genomes reveals multiple copy number variation of short tandem repeats and one transposon-like element

Background

In the present study, we used long-PCR amplification coupled with Next-Generation Sequencing (NGS) to obtain complete mitochondrial (mt) genomes of individual ticks and unprecedently performed precise annotation of these mt genomes. We aimed to: (1) develop a simple, cost-effective and accurate method for the study of extremely high AT-content mt genomes within an individual animal (e.g. Dermacentor silvarum) containing miniscule DNA; (2) provide a high-quality reference genome for D. silvarum with precise annotation and also for future studies of other tick mt genomes; and (3) detect and analyze mt DNA variation within an individual tick.

Results

These annotations were confirmed by the PacBio full-length transcriptome data to cover both entire strands of the mitochondrial genomes without any gaps or overlaps. Moreover, two new and important findings were reported for the first time, contributing fundamental knowledge to mt biology. The first was the discovery of a transposon-like element that may eventually reveal much about mechanisms of gene rearrangements in mt genomes. Another finding was that Copy Number Variation (CNV) of Short Tandem Repeats (STRs) account for mitochondrial sequence diversity (heterogeneity) within an individual tick, insect, mouse or human, whereas SNPs were not detected. The CNV of STRs in the protein-coding genes resulted in frameshift mutations in the proteins, which can cause deleterious effects. Mitochondria containing these deleterious STR mutations accumulate in cells and can produce deleterious proteins.

Conclusions

We proposed that the accumulation of CNV of STRs in mitochondria may cause aging or diseases. Future tests of the CNV of STRs hypothesis help to ultimately reveal the genetic basis of mitochondrial DNA variation and its consequences (e.g., aging and diseases) in animals. Our study will lead to the reconsideration of the importance of STRs and a unified study of CNV of STRs with longer and shorter repeat units (particularly polynucleotides) in both nuclear and mt genomes.

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.

Identification of closely related Ixodes species by protein profiling with MALDI-TOF mass spectrometry

Ticks are vectors of infectious diseases of major importance in human and veterinary medicine. For epidemiological studies, accurate identification of ticks is crucial to define their potential role as vectors and to develop control and prevention strategies. Although morphological and molecular methods are widely used to identify ticks, an innovative approach using MALDI-TOF MS technology recently emerged as an alternative tool. Previous works showed that MALDI-TOF MS was highly effective in identifying ticks, but these works mainly tested tick specimens of different genera. To confirm the accuracy of this new tool for tick identification, nine closely related tick species belonging to the Ixodes genus were analysed, specimens of the Dermacentor reticulatus species were also included in the analysis as an outer group. Three of the species used for the present study belonged to the I. ricinus species complex, which are known to transmit Borrelia burgdorferi sensu lato, the causative agent of Lyme borreliosis. A total of 246 tick specimens were submitted to MALDI-TOF MS analysis, and two body parts (half-idiosoma and four legs) were individually investigated. For each body part, intraspecies reproducibility and interspecies specificity of the MS profiles were determined. The profile analysis revealed that the main determinant for spectra clustering was the tick species for both legs and half-idiosoma. For each body part, a reference database of spectra was set up including 2 to 5 specimens per species randomly selected, and genotyped using 16s rDNA and COI genes to confirm their morphological identification. Both created spectral databases were individually blind tested with their respective body part using the remaining specimens, which were correctly identified in 98.5% of the cases. MALDI-TOF MS is a reliable tool for tick identification, including specimens belonging to closely related species and hardly distinguishable using morphology. The 4-legs as well as the half-idiosoma of ticks can now be applied for specimen identification using two different databases. The combined use of these two body parts improves the rate of tick identification and their confidence level.

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.

A transcriptome-based phylogenetic study of hard ticks (Ixodidae)

Hard ticks are widely distributed across temperate regions, show strong variation in host associations, and are potential vectors of a diversity of medically important zoonoses, such as Lyme disease. To address unresolved issues with respect to the evolutionary relationships among certain species or genera, we produced novel RNA-Seq data sets for nine different Ixodes species. We combined this new data with 18 data sets obtained from public databases, both for Ixodes and non-Ixodes hard tick species, using soft ticks as an outgroup. We assembled transcriptomes (for 27 species in total), predicted coding sequences and identified single copy orthologues (SCO). Using Maximum-likelihood and Bayesian frameworks, we reconstructed a hard tick phylogeny for the nuclear genome. We also obtained a mitochondrial DNA-based phylogeny using published genome sequences and mitochondrial sequences derived from the new transcriptomes. Our results confirm previous studies showing that the Ixodes genus is monophyletic and clarify the relationships among Ixodes sub-genera. This work provides a baseline for studying the evolutionary history of ticks: we indeed found an unexpected acceleration of substitutions for mitochondrial sequences of Prostriata, and for nuclear and mitochondrial genes of two species of Rhipicephalus, which we relate with patterns of genome architecture and changes of life-cycle, respectively.

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.

Argasid and ixodid systematics: Implications for soft tick evolution and systematics, with a new argasid species list

The systematics of the genera and subgenera within the soft tick family Argasidae is not adequately resolved. Different classification schemes, reflecting diverse schools of scientific thought that elevated or downgraded groups to genera or subgenera, have been proposed. In the most recent classification scheme, Argas and Ornithodoros are paraphyletic and the placement of various subgenera remains uncertain because molecular data are lacking. Thus, reclassification of the Argasidae is required. This will enable an understanding of soft tick systematics within an evolutionary context. This study addressed that knowledge gap using mitochondrial genome and nuclear (18S and 28S ribosomal RNA) sequence data for representatives of the subgenera Alectorobius, Argas, Chiropterargas, Ogadenus, Ornamentum, Ornithodoros, Navis (subgen. nov.), Pavlovskyella, Persicargas, Proknekalia, Reticulinasus and Secretargas, from the Afrotropical, Nearctic and Palearctic regions. Hard tick species (Ixodidae) and a new representative of Nuttalliella namaqua (Nuttalliellidae), were also sequenced with a total of 83 whole mitochondrial genomes, 18S rRNA and 28S rRNA genes generated. The study confirmed the utility of next-generation sequencing to retrieve systematic markers. Paraphyly of Argas and Ornithodoros was resolved by systematic analysis and a new species list is proposed. This corresponds broadly with the morphological cladistic analysis of Klompen and Oliver (1993). Estimation of divergence times using molecular dating allowed dissection of phylogeographic patterns for argasid evolution. The discovery of cryptic species in the subgenera Chiropterargas, Ogadenus and Ornithodoros, suggests that cryptic speciation is common within the Argasidae. Cryptic speciation has implications for past biological studies of soft ticks. These are discussed in particular for the Ornithodoros (Ornithodoros) moubata and Ornithodoros (Ornithodoros) savignyi groups.

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

Phylogeny of hard ticks (Ixodidae) remains unresolved. Mitochondrial genomes (mitogenomes) are increasingly used to resolve phylogenetic controversies, but remain unavailable for the entire large Hyalomma genus. Hyalomma asiaticum is a parasitic tick distributed throughout the Asia. As a result of great morphological variability, two subspecies have been recognised historically; until a morphological data-based synonymization was proposed. However, this hypothesis was never tested using molecular data. Therefore, objectives of this study were to: 1. sequence the first Hyalomma mitogenome; 2. scrutinise the proposed synonymization using molecular data, i.e. complete mitogenomes of both subspecies: H. a. asiaticum and kozlovi; 3. conduct phylogenomic and comparative analyses of all available Ixodidae mitogenomes. Results corroborate the proposed synonymization: the two mitogenomes are almost identical (99.6%). Genomic features of both mitogenomes are standard for Metastriata; which includes the presence of two control regions and all three "Tick-Box" motifs. Gene order and strand distribution are perfectly conserved for the entire Metastriata group. Suspecting compositional biases, we conducted phylogenetic analyses (29 almost complete mitogenomes) using homogeneous and heterogeneous (CAT) models of substitution. The results were congruent, apart from the deep-level topology of prostriate ticks (Ixodes): the homogeneous model produced a monophyletic Ixodes, but the CAT model produced a paraphyletic Ixodes (and thereby Prostriata), divided into Australasian and non-Australasian clades. This topology implies that all metastriate ticks have evolved from the ancestor of the non-Australian branch of prostriate ticks. Metastriata was divided into three clades: 1. Amblyomminae and Rhipicephalinae (Rhipicephalus, Hyalomma, Dermacentor); 2. Haemaphysalinae and Bothriocrotoninae, plus Amblyomma sphenodonti; 3. Amblyomma elaphense, basal to all Metastriata. We conclude that mitogenomes have the potential to resolve the long-standing debate about the evolutionary history of ticks, but heterogeneous evolutionary models should be used to alleviate the effects of compositional heterogeneity on deep-level relationships.

Mitochondrial genome evolution and tRNA truncation in Acariformes mites: new evidence from eriophyoid mites

The subclass Acari (mites and ticks) comprises two super-orders: Acariformes and Parasitiformes. Most species of the Parasitiformes known retained the ancestral pattern of mitochondrial (mt) gene arrangement of arthropods, and their mt tRNAs have the typical cloverleaf structure. All of the species of the Acariformes known, however, have rearranged mt genomes and truncated mt tRNAs. We sequenced the mt genomes of two species of Eriophyoidea: Phyllocoptes taishanensis and Epitrimerus sabinae. The mt genomes of P. taishanensis and E. sabinae are 13,475 bp and 13,531 bp, respectively, are circular and contain the 37 genes typical of animals; most mt tRNAs are highly truncated in both mites. On the other hand, these two eriophyoid mites have the least rearranged mt genomes seen in the Acariformes. Comparison between eriophyoid mites and other Aacariformes mites showed that: 1) the most recent common ancestor of Acariformes mites retained the ancestral pattern of mt gene arrangement of arthropods with slight modifications; 2) truncation of tRNAs for cysteine, phenylalanine and histidine occurred once in the most recent common ancestor of Acariformes mites whereas truncation of other tRNAs occurred multiple times; and 3) the placement of eriophyoid mites in the order Trombidiformes needs to be reviewed.

Sequence and secondary structure of the mitochondrial 16S ribosomal RNA gene of Ixodes scapularis

The complete DNA sequences and secondary structure of the mitochondrial (mt) 16S ribosomal (r) RNA gene were determined for six Ixodes scapularis adults. There were 44 variable nucleotide positions in the 1252 bp sequence alignment. Most (95%) nucleotide alterations did not affect the integrity of the secondary structure of the gene because they either occurred at unpaired positions or represented compensatory changes that maintained the base pairing in helices. A large proportion (75%) of the intraspecific variation in DNA sequence occurred within Domains I, II and VI of the 16S gene. Therefore, several regions within this gene may be highly informative for studies of the population genetics and phylogeography of I. scapularis, a major vector of pathogens of humans and domestic animals in North America.

Complete mitochondrial genomes of Carcinoscorpius rotundicauda and Tachypleus tridentatus (Xiphosura, Arthropoda) and implications for chelicerate phylogenetic studies

Horseshoe crabs (order Xiphosura) are often referred to as an ancient order of marine chelicerates and have been considered as keystone taxa for the understanding of chelicerate evolution. However, the mitochondrial genome of this order is only available from a single species, Limulus polyphemus. In the present study, we analyzed the complete mitochondrial genomes from two Asian horseshoe crabs, Carcinoscorpius rotundicauda and Tachypleus tridentatus to offer novel data for the evolutionary relationship within Xiphosura and their position in the chelicerate phylogeny. The mitochondrial genomes of C. rotundicauda (15,033 bp) and T. tridentatus (15,006 bp) encode 13 protein-coding genes, two ribosomal RNA (rRNA) genes, and 22 transfer RNA (tRNA) genes. Overall sequences and genome structure of two Asian species were highly similar to that of Limulus polyphemus, though clear differences among three were found in the stem-loop structure of the putative control region. In the phylogenetic analysis with complete mitochondrial genomes of 43 chelicerate species, C. rotundicauda and T. tridentatus were recovered as a monophyly, while L. polyphemus solely formed an independent clade. Xiphosuran species were placed at the basal root of the tree, and major other chelicerate taxa were clustered in a single monophyly, clearly confirming that horseshoe crabs composed an ancestral taxon among chelicerates. By contrast, the phylogenetic tree without the information of Asian horseshoe crabs did not support monophyletic clustering of other chelicerates. In conclusion, our analyses may provide more robust and reliable perspective on the study of evolutionary history for chelicerates than earlier analyses with a single Atlantic species.

"Candidatus Midichloriaceae" fam. nov. (Rickettsiales), an ecologically widespread clade of intracellular alphaproteobacteria

"Candidatus Midichloria mitochondrii" is an intramitochondrial bacterium of the order Rickettsiales associated with the sheep tick Ixodes ricinus. Bacteria phylogenetically related to "Ca. Midichloria mitochondrii" (midichloria and like organisms [MALOs]) have been shown to be associated with a wide range of hosts, from amoebae to a variety of animals, including humans. Despite numerous studies focused on specific members of the MALO group, no comprehensive phylogenetic and statistical analyses have so far been performed on the group as a whole. Here, we present a multidisciplinary investigation based on 16S rRNA gene sequences using both phylogenetic and statistical methods, thereby analyzing MALOs in the overall framework of the Rickettsiales. This study revealed that (i) MALOs form a monophyletic group; (ii) the MALO group is structured into distinct subgroups, verifying current genera as significant evolutionary units and identifying several subclades that could represent novel genera; (iii) the MALO group ranks at the level of described Rickettsiales families, leading to the proposal of the novel family "Candidatus Midichloriaceae." In addition, based on the phylogenetic trees generated, we present an evolutionary scenario to interpret the distribution and life history transitions of these microorganisms associated with highly divergent eukaryotic hosts: we suggest that aquatic/environmental protista have acted as evolutionary reservoirs for members of this novel family, from which one or more lineages with the capacity of infecting metazoa have evolved.