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

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.

Aberrant Mitochondrial tRNA Genes Appear Frequently in Animal Evolution

Mitochondrial tRNAs have acquired a diverse portfolio of aberrant structures throughout metazoan evolution. With the availability of more than 12,500 mitogenome sequences, it is essential to compile a comprehensive overview of the pattern changes with regard to mitochondrial tRNA repertoire and structural variations. This, of course, requires reanalysis of the sequence data of more than 250,000 mitochondrial tRNAs with a uniform workflow. Here, we report our results on the complete reannotation of all mitogenomes available in the RefSeq database by September 2022 using mitos2. Based on the individual cases of mitochondrial tRNA variants reported throughout the literature, our data pinpoint the respective hotspots of change, i.e. Acanthocephala (Lophotrochozoa), Nematoda, Acariformes, and Araneae (Arthropoda). Less dramatic deviations of mitochondrial tRNAs from the norm are observed throughout many other clades. Loss of arms in animal mitochondrial tRNA clearly is a phenomenon that occurred independently many times, not limited to a small number of specific clades. The summary data here provide a starting point for systematic investigations into the detailed evolutionary processes of structural reduction and loss of mitochondrial tRNAs as well as a resource for further improvements of annotation workflows for mitochondrial tRNA annotation.

Microbiome and mitogenomics of the chigger mite Pentidionis agamae: potential role as an Orientia vector and associations with divergent clades of Wolbachia and Borrelia

BACKGROUND

Trombiculid mites are globally distributed, highly diverse arachnids that largely lack molecular resources such as whole mitogenomes for the elucidation of taxonomic relationships. Trombiculid larvae (chiggers) parasitise vertebrates and can transmit bacteria (Orientia spp.) responsible for scrub typhus, a zoonotic febrile illness. Orientia tsutsugamushi causes most cases of scrub typhus and is endemic to the Asia-Pacific Region, where it is transmitted by Leptotrombidium spp. chiggers. However, in Dubai, Candidatus Orientia chuto was isolated from a case of scrub typhus and is also known to circulate among rodents in Saudi Arabia and Kenya, although its vectors remain poorly defined. In addition to Orientia, chiggers are often infected with other potential pathogens or arthropod-specific endosymbionts, but their significance for trombiculid biology and public health is unclear.

RESULTS

Ten chigger species were collected from rodents in southwestern Saudi Arabia. Chiggers were pooled according to species and screened for Orientia DNA by PCR. Two species (Microtrombicula muhaylensis and Pentidionis agamae) produced positive results for the htrA gene, although Ca. Orientia chuto DNA was confirmed by Sanger sequencing only in P. agamae. Metagenomic sequencing of three pools of P. agamae provided evidence for two other bacterial associates: a spirochaete and a Wolbachia symbiont. Phylogenetic analysis of 16S rRNA and multi-locus sequence typing genes placed the spirochaete in a clade of micromammal-associated Borrelia spp. that are widely-distributed globally with no known vector. For the Wolbachia symbiont, a genome assembly was obtained that allowed phylogenetic localisation in a novel, divergent clade. Cytochrome c oxidase I (COI) barcodes for Saudi Arabian chiggers enabled comparisons with global chigger diversity, revealing several cases of discordance with classical taxonomy. Complete mitogenome assemblies were obtained for the three P. agamae pools and almost 50 SNPs were identified, despite a common geographic origin.

CONCLUSIONS

P. agamae was identified as a potential vector of Ca. Orientia chuto on the Arabian Peninsula. The detection of an unusual Borrelia sp. and a divergent Wolbachia symbiont in P. agamae indicated links with chigger microbiomes in other parts of the world, while COI barcoding and mitogenomic analyses greatly extended our understanding of inter- and intraspecific relationships in trombiculid mites.

The complete mitochondrial genome of Echinolaelaps fukienensis provide insights into phylogeny and rearrangement in the superfamily Dermanyssoidea

BACKGROUND

Echinolaelaps fukienensis is the dominant mite species parasitic on the body surface of the genus Niviventer. The mitochondrial genome (mitogenome) has its own independent genetic material and genetic system, and is now widely used in population genetics, genealogical biogeography, phylogeny and molecular evolution studies. Species diversity of the superfamily Dermanyssoidea is very rich, but its mitogenomes AT content is high, and it is difficult to amplify the complete mitogenome by routine PCR. To date, we have only obtained the mitogenomes of 6 species, scarcity on sequence data has greatly impeded the studies in the superfamily Dermanyssoidea.

METHODS

Echinolaelaps fukienensis were collected in 2019 from the body surface of Niviventer confucianus (Rodentia, Muridae) in Yunnan Province. The E. fukienensis mitogenome was determined and analyzed for the first time using the Illumina Novoseq 6000 platform. Phylogenetic analyses of the superfamily Dermanyssoidea were conducted based on the entire mitogenome sequences.

RESULTS

The E. fukienensis mitogenome was 14,402 bp, which is known the smallest genome of the superfamily Dermanyssoidea, encoding a total of 37 genes, including 13 PCGs, 22 tRNAs, 2 rRNAs and 1 control region. Most protein-coding genes use ATN as the start codon and TAN as the stop codon. AT and GC skew of atp8 genes in E. fukienensis were both 0. The average length of 22 tRNA genes of E. fukienensis was 64 bp, and secondary structures of tRNAs showed base mismatches and missing D-arms in many places. Compared with gene arrangement pattern of the hypothetical ancestor of arthropods, the E. fukienensis mitogenome shows a novel arrangement pattern. Phylogenetic tree supported the monophyly of the superfamily Dermanyssoidea. Echinolaelaps fukienensis being the least genetic distant (0.2762) and most closely related to Varroa destructor.

CONCLUSIONS

This study analyzed comprehensive the structure and evolution of the E. fukienensis mitogenome for the first time, enriches molecular data of the genus Echinolaelaps, which will contribute to further understand phylogeny and rearrangement patterns of the superfamily Dermanyssoidea.

The evolutionary characterization of Gamasida based on mitochondrial genes codon usage pattern

Mites belonging to the suborder Gamasida are species-rich and habitat-diverse, with a worldwide distribution. To adapt to the environment and obtain better living conditions, all species of the suborder Gamasida have been undergoing constant evolution. The complete mitochondrial genome (mitogenome) is an invaluable molecular marker for studying the origin of species, genetic differentiation between closely related species, and between intraspecific groups. In some species of the suborder Gamasida, mitochondrial tRNA genes are truncated and carried unstable genetic information. This study presents a comparative analysis of codon usage pattern and preference of 13 protein-coding genes of 24 species in 17 genera and 10 families of the suborder Gamasida. Results showed that have an obvious AT preference (0.664-0.829) for codon usage in the suborder Gamasida. Most of the optimal and high-frequency codons also end in A/T. The degree of natural selection varies between the same protein-coding genes of different gamasid mites or among different protein-coding genes within the same gamasid mites. Base and codon usage pattern and preference are very similar between the same species and genus, namely the closer species, the more similar their bases and codons usage patterns and preference are. T bases and C bases were the preference bases for codon usage of 24 species in the suborder Gamasida. Evolution of the suborder Gamasida was dominated by natural selection (64.1%). This study provides the first comprehensive analysis of codon usage in the suborder Gamasida, which will greatly improve our understanding of codon usage patterns and preference, genetics, and evolution of the suborder Gamasida. It will help to evaluate the degree of molecular adaptation in the suborder Gamasida and to further explore evolutionary features of the suborder Gamasida.

Atypically Shaped Setae in Gall Mites (Acariformes, Eriophyoidea) and Mitogenomics of the Genus Leipothrix Keifer (Eriophyidae)

The setae in Eriophyoidea are filiform, slightly bent and thickened near the base. Confocal microscopy indicates that their proximal and distal parts differ in light reflection and autofluorescence. Approximately 50 genera have atypically shaped setae: bifurcated, angled or swollen. These modifications are known in the basal part of prosomal setae u', ft', ft″, d, v, bv, ve, sc and caudal setae h2. We assessed the distribution of atypically shaped setae in Eriophyoidea and showed that they are scattered in different phylogenetic lineages. We hypothesized that the ancestral setae of eriophyoid mites were bifurcated before later simplifying into filiform setae. We also proposed that hypo-furcating setae are a synapomorphy that unites Eriophyoidea with Nematalycidae. We analyzed four new mitochondrial genomes of Leipothrix, the largest genus with bifurcated d, and showed that it is monophyletic and has a unique mitochondrial gene order with translocated trnK. We exclude Cereusacarus juniperensisn. comb. Xue and Yin, 2020 from Leipothrix and transfer five Epitrimerus spp. to Leipothrix: L. aegopodii (Liro 1941) n. comb., L. femoralis (Liro 1941) n. comb., L. geranii (Liro 1941) n. comb., L. ranunculi (Liro 1941) n. comb., and L. triquetra (Meyer 1990) n. comb.

Where Eriophyoidea (Acariformes) Belong in the Tree of Life

Over the past century and a half, the taxonomic placement of Eriophyoidea has been in flux. For much of this period, this group has been treated as a subtaxon within Trombidiformes. However, the vast majority of recent phylogenetic analyses, including almost all phylogenomic analyses, place this group outside Trombidiformes. The few studies that still place Eriophyoidea within Trombidiformes are likely to be biased by incomplete taxon/gene sampling, long branch attraction, the omission of RNA secondary structure in sequence alignment, and the inclusion of hypervariable expansion-contraction rRNA regions. Based on the agreement among a number of independent analyses that use a range of different datasets (morphology; multiple genes; mitochondrial/whole genomes), Eriophyoidea are almost certain to be closely related to Nematalycidae, a family of vermiform mites within Endeostigmata, a basal acariform grade. Much of the morphological evidence in support of this relationship was apparent after the discovery of Nematalycidae in the middle of the 20th century. However, this evidence has largely been disregarded until very recently, perhaps because of overconfidence in the placement of Eriophyoidea within Trombidiformes. Here, we briefly review and identify a number of biases, both molecular- and morphology-based, that can lead to erroneous reconstructions of the position of Eriophyoidea in the tree of life.

Integrative Taxonomy of the Gall Mite Nothopoda todeica n. sp. (Eriophyidae) from the Disjunct Afro-Australasian Fern Todea barbara: Morphology, Phylogeny, and Mitogenomics

Eriophyoidea is a group of phytoparasitic mites with poorly resolved phylogeny. Previous studies inferred Eriophyidae s.l. as the largest molecular clade of Eriophyoidea, and Nothopodinae as the basal divergence of Eriophyidae s.l. We investigate the morphology and molecular phylogeny of Nothopoda todeican. sp. (Nothopodinae, Nothopodini), associated with a disjunct Afro-Australasian fern Todea barbara (Osmundaceae) from South Africa. Our analyses (1) determine new erroneous sequences (KF782375, KF782475, KF782586) wrongly assigned to Nothopodinae instead of Phyllocoptinae, (2) confirm the basal position of Nothopodinae in Eriophyoidea s.l., (3) question the monophyly of the Colopodacini and Nothopodini tribes, and (4) show the nested position of African fern-associated Nothopoda within a clade dominated by Asian nothopodines from angiosperms, which implies (a) a secondary association of nothopodines with ferns and (b) no relation between geography (continents) and the phylogenetic relationships of Nothopodinae species. Finally, we obtained a first complete mitochondrial genome for Nothopodinae and revealed a new gene order in the mitogenome of N. todeican. sp., notably deviating from those in other investigated eriophyoids. Our results contribute to resolving the phylogeny of Eriophyoidea and provide an example of an integrative study of a new taxon belonging to an economically important group of acariform mites.

One-way ticket to the blue: A large-scale, dated phylogeny revealed asymmetric land-to-water transitions in acariform mites (Acari: Acariformes)

Acariform mites are an ancient and megadiverse lineage that may have experienced a complex pattern of invasions into terrestrial and aquatic habitats. These among-realm transitions may relate to periods of turmoil in Earth's history or be simply results of uneven biodiversity patterns across habitats. Here, we inferred a dated, representative acariform phylogeny (five genes, 9,200 bp aligned, 367 terminals belonging to 150 ingroup plus 15 outgroup families, 23 fossil calibration points) which was used to infer transitions between marine/freshwater/terrestrial habitats. We detected four unambiguous transitions from terrestrial to freshwater habitats (Hydrozetes, Naiadacarus, Fusohericia, Afronothrus, Homocaligus); one from freshwater to marine (Pontarachnidae), and four from marine to brackish or freshwater transitions (all among Halacaridae: Acarothrix; Halacarellus petiti; Copidognathus sp.; clade Limnohalacarus + Soldanellonyx + Porohalacarus + Porolohmannella). One transition to the sea was inferred ambiguously with respect to the ancestor being either terrestrial or freshwater (Hyadesiidae), and another must be most carefully examined by adding potential related taxa (Selenoribatidae + Fortuyniidae). Finally, we inferred a single, remarkable transition from aquatic to terrestrial habitats involving early evolution of the large and ecologically diverse lineage: the ancestor of the Halacaridae + Parasitengona clade was probably freshwater given our dataset, thus making terrestrial Parasitengona secondarily terrestrial. Overall, our results suggested a strong asymmetry in environmental transitions: the majority occurred from terrestrial to aquatic habitats. This asymmetry is probably linked to mites' biological properties and uneven biodiversity patterns across habitats rather than Earth's geological history. Since the land holds more acariform diversity than water habitats, a shift from the former is more likely than from the latter. We inferred the following relationships: alicid endeostigmatid + eriophyoid (Alycidae, (Nanorchestidae, (Nematalycidae, Eriophyoidea))) being sister group to the remaining Acariformes: (proteonematalycid Endeostigmata, alicorhagiid Endeostigmata, Trombidiformes, Oribatida (including Astigmata)). Trombidiform relationships had several novel rearrangements: (i) traditional Eupodina lacked support for the inclusion of Bdelloidea; (ii) Teneriffidae, traditionally placed among Anystina, was consistently recovered in a clade including Heterostigmata in Eleutherengona; (iii) several lineages, such as Adamystidae, Paratydeidae, Caeculidae and Erythracaridae, were recovered in a large clade along other Anystina and Eleutherengona, suggesting single origins of several fundamental character states, such as the reduction of the cheliceral fixed digit and development of the palpal thumb-claw complex.

Complete mitochondrial genomes of Thyreophagus entomophagus and Acarus siro (Sarcoptiformes: Astigmatina) provide insight into mitogenome features, evolution, and phylogeny among Acaroidea mites

Mites from the Acaroidea (Sarcoptiformes: Astigmatina) are important pests of various stored products, posing potential threats to preserved foods. In addition, mites can cause allergic diseases. Complete mitochondrial genomes (mitogenomes) are valuable resources for different research fields, including comparative genomics, molecular evolutionary analysis, and phylogenetic inference. We sequenced and annotated the complete mitogenomes of Thyreophagus entomophagus and Acarus siro. A comparative analysis was made between mitogenomic sequences from 10 species representing nine genera within Acaroidea. The mitogenomes of T. entomophagus and A. siro contained 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region. In Acaroidea species, mitogenomes have highly conserved gene size and order, and codon usage. Among Acaroidea mites, most PCGs were found to be under purifying selection, implying that most PCGs might have evolved slowly. Our findings showed that nad4 evolved most rapidly, whereas cox1 and cox3 evolved most slowly. The evolutionary rates of Acaroidea vary considerably across families. In addition, selection analyses were also performed in 23 astigmatid mite species, and the evolutionary rate of the same genes in different superfamilies exhibited large differences. Phylogenetic results are mostly consistent with those identified by previous phylogenetic studies on astigmatid mites. The monophyly of Acaroidea was rejected, and the Suidasiidae and Lardoglyphidae appeared to deviate from the Acaroidea branch. Our research proposed a review of the current Acaroidea classification system.

Positive Correlation of the Gene Rearrangements and Evolutionary Rates in the Mitochondrial Genomes of Thrips (Insecta: Thysanoptera)

Simple Summary

Aeolothrips, commonly known as banded thrips, is the largest genus of the family Aeolothripidae (predatory thrips). In the current study, we sequenced the mitochondrial genome (mitogenome) of the banded thrip species Aeolothrips xinjiangensis. We found a novel gene arrangement in this mitogenome that has not been reported in Thysanoptera. By comparing the gene order and rearrangement patterns, we found seven identical gene blocks and three identical rearrangement events in two mitogenomes of banded thrips. There was marked variation in the mitochondrial gene order across thrip species, with only two conserved gene blocks shared by all 14 thrips. In addition, we found a positive correlation between the degree of gene rearrangement and evolutionary rate. Our results suggested that the mitogenomes of thrips have tended to be stable since their massive rearrangement.

Abstract

Extensive gene rearrangement is characteristic in the mitogenomes of thrips (Thysanoptera), but the historical process giving rise to the contemporary gene rearrangement pattern remains unclear. To better understand the evolutionary processes of gene rearrangement in the mitogenomes of thrips, we sequenced the mitogenome of the banded thrip species Aeolothrips xinjiangensis. First, we found a novel mitochondrial gene order in this species. This mitogenome is 16,947 bp in length and encodes the typical 37 coding genes (13 protein-coding genes, 22 tRNA genes, and two rRNA genes) of insects. The gene arrangement was dramatically different from the putative ancestral mitogenome, with 26 genes being translocated, eight of which were inverted. Moreover, we found a novel, conserved gene block, trnC-trnY, which has not been previously reported in the mitogenomes of thrips. With this newly assembled mitogenome, we compared mitogenome sequences across Thysanoptera to assess the evolutionary processes giving rise to the current gene rearrangement pattern in thrips. Seven identical gene blocks were shared by two sequenced banded thrip mitogenomes, while the reversal of ND2 combined with TDRL events resulted in the different gene orders of these two species. In phylogenetic analysis, the monophyly of the suborders and families of Thysanoptera was well supported. Across the gene orders of 14 thrips, only two conserved gene blocks, ATP8-ATP6 and ND4-ND4L, could be found. Correlation analysis showed that the degree of gene rearrangement was positively correlated with the non-synonymous substitution rate in thrips. Our study suggests that the mitogenomes of thrips remain stable over long evolutionary timescales after massive rearrangement during early diversification.

Symbiotic bacteria of the gall-inducing mite Fragariocoptes setiger (Eriophyoidea) and phylogenomic resolution of the eriophyoid position among Acari

Eriophyoid mites represent a hyperdiverse, phytophagous lineage with an unclear phylogenetic position. These mites have succeeded in colonizing nearly every seed plant species, and this evolutionary success was in part due to the mites' ability to induce galls in plants. A gall is a unique niche that provides the inducer of this modification with vital resources. The exact mechanism of gall formation is still not understood, even as to whether it is endogenic (mites directly cause galls) or exogenic (symbiotic microorganisms are involved). Here we (i) investigate the phylogenetic affinities of eriophyoids and (ii) use comparative metagenomics to test the hypothesis that the endosymbionts of eriophyoid mites are involved in gall formation. Our phylogenomic analysis robustly inferred eriophyoids as closely related to Nematalycidae, a group of deep-soil mites belonging to Endeostigmata. Our comparative metagenomics, fluorescence in situ hybridization, and electron microscopy experiments identified two candidate endosymbiotic bacteria shared across samples, however, it is unlikely that they are gall inducers (morphotype1: novel Wolbachia, morphotype2: possibly Agrobacterium tumefaciens). We also detected an array of plant pathogens associated with galls that may be vectored by the mites, and we determined a mite pathogenic virus (Betabaculovirus) that could be tested for using in biocontrol of agricultural pest mites.

The mitogenome of Halotydeus destructor (Tucker) and its relationships with other trombidiform mites as inferred from nucleotide sequences and gene arrangements

The redlegged earth mite, Halotydeus destructor (Tucker, 1925: Trombidiformes, Eupodoidea, Penthaleidae), is an invasive mite species. In Australia, this mite has become a pest of winter pastures and grain crops. We report the complete mitogenome for H. destructor, the first to represent the family Penthaleidae, superfamily Eupodoidea. The mitogenome of H. destructor is 14,691 bp in size, and has a GC content of 27.87%, 13 protein-coding genes, two rRNA genes, and 22 tRNA genes. We explored evolutionary relationships of H. destructor with other members of the Trombidiformes using phylogenetic analyses of nucleotide sequences and the order of protein-coding and rRNA genes. We found strong, consistent support for the superfamily Tydeoidea being the sister taxon to the superfamily Eupodoidea based on nucleotide sequences and gene arrangements. Moreover, the gene arrangements of Eupodoidea and Tydeoidea are not only identical to each other but also identical to that of the hypothesized arthropod ancestor, showing a high level of conservatism in the mitogenomic structure of these mite superfamilies. Our study illustrates the utility of gene arrangements for providing complementary information to nucleotide sequences with respect to inferring the evolutionary relationships of species within the order Trombidiformes. The mitogenome of H. destructor provides a valuable resource for further population genetic studies of this important agricultural pest. Given the co-occurrence of closely related, morphologically similar Penthaleidae mites with H. destructor in the field, a complete mitogenome provides new opportunities to develop metabarcoding tools to study mite diversity in agro-ecosystems. Moreover, the H. destructor mitogenome fills an important taxonomic gap that will facilitate further study of trombidiform mite evolution.

Mitochondrial analysis of oribatid mites provides insights into their atypical tRNA annotation, genome rearrangement and evolution

Background

The mitochondrial (mt) genomes of Sarcoptiformes mites typically contain 37 genes. Although the loss of genes is rare in Sarcoptiformes mite mitogenomes, two of the six previously reported oribatid mites (Acariforms: Sarcoptiformes) are reported to have lost parts of their tRNA genes. To confirm whether the tRNA genes were indeed lost and whether the loss is universal, we re-annotated the available oribatid mite sequences and sequenced the mitogenome of Oribatula sakamorii.

Methods

The mitogenome of O. sakamorii was sequenced using an Illumina HiSeq sequencer. The mt tRNA gene was annotated using multi-software combined with a manual annotation approach. Phylogenetic analyses were performed using the maximum likelihood and Bayesian inference methods with concatenated nucleotide and amino acid sequences.

Results

The mitogenomes of O. sakamorii contained 37 genes, including 22 tRNA genes. We identified all mt tRNA genes that were reported as “lost” in Steganacarus magnus and Paraleius leontonychus and revealed certain atypical tRNA annotation errors in oribatid mite sequences. Oribatid mite mitogenomes are characterized by low rates of genetic rearrangement, with six or seven gene blocks conserved between the mitogenome of all species and that of ancestral arthropods. Considering the relative order of the major genes (protein-coding genes and rRNAs), only one or two genes were rearranged with respect to their positions in the ancestral genome. We explored the phylogenetic relationships among the available oribatid mites, and the results confirmed the systematic position of Hermannia in the Crotonioidea superfamily. This was also supported by the synapomorphic gene-derived boundaries.

Conclusions

The tRNA “lost” phenomenon is not universal in oribatid mites. Rather, highly atypical secondary structure of the inferred mt tRNA genes made them unidentifiable using a single type of tRNA search program. The use of multi-software combined with a manual annotation approach can improve the accuracy of tRNA gene annotation. In addition, we identified the precise systematic position of Hermannia and validated that Astigmata is nested in Oribatida.

Graphic Abstract

Supplementary Information

The online version contains supplementary material available at 10.1186/s13071-021-04719-0.

Hopeful monsters: unintended sequencing of famously malformed mite mitochondrial tRNAs reveals widespread expression and processing of sense-antisense pairs

Although tRNA structure is one of the most conserved and recognizable shapes in molecular biology, aberrant tRNAs are frequently found in the mitochondrial genomes of metazoans. The extremely degenerate structures of several mitochondrial tRNAs (mt-tRNAs) have led to doubts about their expression and function. Mites from the arachnid superorder Acariformes are predicted to have some of the shortest mt-tRNAs, with a complete loss of cloverleaf-like shape. While performing mitochondrial isolations and recently developed tRNA-seq methods in plant tissue, we inadvertently sequenced the mt-tRNAs from a common plant pest, the acariform mite Tetranychus urticae, to a high enough coverage to detect all previously annotated T. urticae tRNA regions. The results not only confirm expression, CCA-tailing and post-transcriptional base modification of these highly divergent tRNAs, but also revealed paired sense and antisense expression of multiple T. urticae mt-tRNAs. Mirrored expression of mt-tRNA genes has been hypothesized but not previously demonstrated to be common in any system. We discuss the functional roles that these divergent tRNAs could have as both decoding molecules in translation and processing signals in transcript maturation pathways, as well as how sense–antisense pairs add another dimension to the bizarre tRNA biology of mitochondrial genomes.

Mitochondrial Metagenomics Reveals the Ancient Origin and Phylodiversity of Soil Mites and Provides a Phylogeny of the Acari

High-throughput DNA methods hold great promise for phylogenetic analysis of lineages that are difficult to study with conventional molecular and morphological approaches. The mites (Acari), and in particular the highly diverse soil-dwelling lineages, are among the least known branches of the metazoan Tree-of-Life. We extracted numerous minute mites from soils in an area of mixed forest and grassland in southern Iberia. Selected specimens representing the full morphological diversity were shotgun sequenced in bulk, followed by genome assembly of short reads from the mixture, which produced >100 mitochondrial genomes representing diverse acarine lineages. Phylogenetic analyses in combination with taxonomically limited mitogenomes available publicly resulted in plausible trees defining basal relationships of the Acari. Several critical nodes were supported by ancestral-state reconstructions of mitochondrial gene rearrangements. Molecular calibration placed the minimum age for the common ancestor of the superorder Acariformes, which includes most soil-dwelling mites, to the Cambrian-Ordovician (likely within 455-552 Ma), whereas the origin of the superorder Parasitiformes was placed later in the Carboniferous-Permian. Most family-level taxa within the Acariformes were dated to the Jurassic and Triassic. The ancient origin of Acariformes and the early diversification of major extant lineages linked to the soil are consistent with a pioneering role for mites in building the earliest terrestrial ecosystems.

Characterization of the complete mitochondrial genome of the predatory mite Stratiolaelaps scimitus (Acari: Laelapidae)

In this study, we recovered the complete mitochondrial genome of Stratiolaelaps scimitus through Illumina sequencing data. The circularized mitogenome is 16,009 in length, which consists of 37 genes (13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes). The overall base composition is 35.4% for A, 40.1% for T, 16.8% for G, 7.7% for C, demonstrating an extreme bias of high AT content (75.5%). The whole mitogenome of S. scimitus and other Acari mitogenomes (10 species, in total) were used for phylogenetic analysis, and the result showed that the relationship of S. scimitus was close to Varroa destructor in the same Superfamily Dermanyssoidea.

Genomic insights into mite phylogeny, fitness, development, and reproduction

BACKGROUND

Predatory mites (Acari: Phytoseiidae) are the most important beneficial arthropods used in augmentative biological pest control of protected crops around the world. However, the genomes of mites are far less well understood than those of insects and the evolutionary relationships among mite and other chelicerate orders are contested, with the enigmatic origin of mites at one of the centres in discussion of the evolution of Arachnida.

RESULTS

We here report the 173 Mb nuclear genome (from 51.75 Gb pairs of Illumina reads) of the predatory mite, Neoseiulus cucumeris, a biocontrol agent against pests such as mites and thrips worldwide. We identified nearly 20.6 Mb (~ 11.93% of this genome) of repetitive sequences and annotated 18,735 protein-coding genes (a typical gene 2888 bp in size); the total length of protein-coding genes was about 50.55 Mb (29.2% of this assembly). About 37% (6981) of the genes are unique to N. cucumeris based on comparison with other arachnid genomes. Our phylogenomic analysis supported the monophyly of Acari, therefore rejecting the biphyletic origin of mites advocated by other studies based on limited gene fragments or few taxa in recent years. Our transcriptomic analyses of different life stages of N. cucumeris provide new insights into genes involved in its development. Putative genes involved in vitellogenesis, regulation of oviposition, sex determination, development of legs, signal perception, detoxification and stress-resistance, and innate immune systems are identified.

CONCLUSIONS

Our genomics and developmental transcriptomics analyses of N. cucumeris provide invaluable resources for further research on the development, reproduction, and fitness of this economically important mite in particular and Arachnida in general.

Unravelling the phylogeny, cryptic diversity and morphological evolution of Diptilomiopus mites (Acari: Eriophyoidea)

The Eriophyoidea, notable for specific morphological characters (four-legged mites) and gall-formation in host plants (gall mites), is one of the most species-rich superfamilies of Acari. Monophyly of the superfamily Eriophyoidea is accepted by all acarologists; however, monophyly of most genera has not been evaluated in a molecular phylogenetic network. Furthermore, most eriophyoid mites, especially species in the genus Diptilomiopus, are morphologically similar, challenging their identification. Here we test the phylogeny and cryptic diversity of Diptilomiopus species using fragments of two mitochondrial (COI and 12S) and two nuclear (18S and 28S) genes. Our results revealed the monophyly of Diptilomiopus. Sequence distance, barcode gap, and species delimitation analyses of the COI gene allowed us to resolve cryptic diversity of Diptilomiopus species. Additionally, we supposed that characteristics of genu fused with femur on both legs and seta ft' absent on leg II evolved only once within Diptilomiopus, which are potential morphological synapomorphies. In contrast, characteristics of both setae ft' and ft″ divided into a short branch and a long branch were supposed evolving multiple times independently. Our findings contribute to the understanding of phylogeny and morphological evolution of Diptilomiopus species and provide a DNA-based approach for species delimitation of Diptilomiopus mites.

Arm-less mitochondrial tRNAs conserved for over 30 millions of years in spiders

BACKGROUND

In recent years, Next Generation Sequencing (NGS) has accelerated the generation of full mitogenomes, providing abundant material for studying different aspects of molecular evolution. Some mitogenomes have been observed to harbor atypical sequences with bizarre secondary structures, which origins and significance could only be fully understood in an evolutionary framework.

RESULTS

Here we report and analyze the mitochondrial sequences and gene arrangements of six closely related spiders in the sister genera Parachtes and Harpactocrates, which belong to the nocturnal, ground dwelling family Dysderidae. Species of both genera have compacted mitogenomes with many overlapping genes and strikingly reduced tRNAs that are among the shortest described within metazoans. Thanks to the conservation of the gene order and the nucleotide identity across close relatives, we were able to predict the secondary structures even on arm-less tRNAs, which would be otherwise unattainable for a single species. They exhibit aberrant secondary structures with the lack of either DHU or TΨC arms and many miss-pairings in the acceptor arm but this degeneracy trend goes even further since at least four tRNAs are arm-less in the six spider species studied.

CONCLUSIONS

The conservation of at least four arm-less tRNA genes in two sister spider genera for about 30 myr suggest that these genes are still encoding fully functional tRNAs though they may be post-transcriptionally edited to be fully functional as previously described in other species. We suggest that the presence of overlapping and truncated tRNA genes may be related and explains why spider mitogenomes are smaller than those of other invertebrates.

Mitochondrial genome reorganization provides insights into the relationship between oribatid mites and astigmatid mites (Acari: Sarcoptiformes: Oribatida)

Oribatida s.l. represents one of the most species-rich mite lineages, including two recognized groups: oribatid mites (Oribatida s.s., non-astigmatan oribatids) and astigmatid mites (Astigmata). However, the relationship between these two groups has been debated. Here, we sequenced the complete mitochondrial (mt) genome of one oribatid mite and one astigmatid mite, retrieved complete mt genomes of three oribatid mites, and compared them with two other oribatid mites and 12 astigmatid mites sequenced previously. We find that gene orders in the mt genomes of both oribatid mites and astigmatid mites are rearranged relative to the hypothetical ancestral arrangement of the arthropods. Based on the shared derived gene clusters in each mt genome group, rearranged mt genomes are roughly divided into two groups corresponding to each mite group (oribatid mites or astigmatid mites). Phylogenetic results show that Astigmata nested in Oribatida. The monophyly of Astigmata is recovered, while paraphyly of Oribatida s.s. is observed. Our results show that rearranged gene orders in the mt genomes characterize various lineages of oribatid mites and astigmatid mites, and have potential phylogenetic information for resolving the high-level (cohort or supercohort) phylogeny of Oribatida.

A new, sensitive and efficient method for taxonomic placement in the Eriophyoidea and virus detection in individual eriophyoids

Eriophyoids affect crops around the globe directly or indirectly as virus vectors. Eriophyoid systematics initiated over a century ago, yet more than 90% of their fauna remain undescribed. Morphological identification is challenging because of a limited number of traits, cryptic speciation and complex life cycle reported for many species in the group. Nucleic acids extraction for mite identification is challenging due to their microscopic size with researchers using pooled samples leading to polymorphisms and inconclusive results. Identification of mite virus vectors is a tiresome task that could be simplified with a protocol that allows for the detection of viruses in the individual specimen. This communication describes an innovative, highly efficient extraction and detection pipeline. Direct Reverse Transcriptase - Polymerase Chain Reaction (Drt-PCR) assays were implemented in the molecular identification of eriophyoids and detection of viruses present in their bodies. The reverse transcription step allows for amplification from a single mite or egg, as in addition to the genomic DNA, it incorporates the abundant transcripts of targeted genes, whereas it also allows for the amplification of viruses. This communication provides an efficient, sensitive and cost-effective alternative that can be implemented in pest identification and detection as well as biological and ecological studies.

Going Deeper into High and Low Phylogenetic Relationships of Protura

Proturans are small, wingless, soil-dwelling arthropods, generally associated with the early diversification of Hexapoda. Their bizarre morphology, together with conflicting results of molecular studies, has nevertheless made their classification ambiguous. Furthermore, their limited dispersal capability (due to the primarily absence of wings) and their euedaphic lifestyle have greatly complicated species-level identification. Mitochondrial and nuclear markers have been applied herein to investigate and summarize proturan systematics at different hierarchical levels. Two new mitochondrial genomes are described and included in a phylum-level phylogenetic analysis, but the position of Protura could not be resolved with confidence due to an accelerated rate of substitution and extensive gene rearrangements. Mitochondrial and nuclear loci were also applied in order to revise the intra-class systematics, recovering three proturan orders and most of the families/subfamilies included as monophyletic, with the exception of the subfamily Acerentominae. At the species level, most morphologically described species were confirmed using molecular markers, with some exceptions, and the advantages of including nuclear, as well as mitochondrial, markers and morphology are discussed. At all levels, an enlarged taxon sampling and the integration of data from different sources may be of significant help in solving open questions that still persist on the evolutionary history of Protura.

The mitochondrial genomes of sarcoptiform mites: are any transfer RNA genes really lost?

BACKGROUND

Mitochondrial (mt) genomes of animals typically contain 37 genes for 13 proteins, two ribosomal RNA (rRNA) genes and 22 transfer RNA (tRNA) genes. In sarcoptiform mites, the entire set of mt tRNA genes is present in Aleuroglyphus ovatus, Caloglyphus berlesei, Dermatophagoides farinae, D. pteronyssinus, Histiostoma blomquisti and Psoroptes cuniculi. Loss of 16 mt tRNA genes, however, was reported in Steganacarus magnus; loss of 2-3 tRNA genes was reported in Tyrophagus longior, T. putrescentiae and Sarcoptes scabiei. Nevertheless, convincing evidence for mt gene loss is lacking in these mites.

RESULTS

We sequenced the mitochondrial genomes of two sarcoptiform mites, Histiostoma feroniarum (13,896 bp) and Rhizoglyphus robini (14,244 bp). Using tRNAScan and ARWEN programs, we identified 16 and 17 tRNA genes in the mt genomes of H. feroniarum and R. robini, respectively. The other six mt tRNA genes in H. feroniarum and five mt tRNA genes in R. robini can only be identified manually by sequence comparison when alternative anticodons are considered. We applied this manual approach to other mites that were reported previously to have lost mt tRNA genes. We were able to identify all of the 16 mt tRNA genes that were reported as lost in St. magnus, two of the three mt tRNA genes that were reported as lost in T. longior and T. putrescentiae, and the two mt tRNA genes that were reported as lost in Sa. scabiei. All of the tRNA genes inferred from these manually identified genes have truncation in the arms and mismatches in the stems.

CONCLUSIONS

Our results reveal very unconventional tRNA structures in sarcoptiform mites and do not support the loss of mt tRNA genes in these mites. The functional implication of the drastic structural changes in these tRNA genes remains to be investigated.

The mitochondrial genome of the oribatid mite Paraleius leontonychus: new insights into tRNA evolution and phylogenetic relationships in acariform mites

Bilaterian mitochondrial (mt) genomes are circular molecules that typically contain 37 genes. To date, only a single complete mitogenome sequence is available for the species-rich sarcoptiform mite order Oribatida. We sequenced the mitogenome of Paraleius leontonychus, another species of this suborder. It is 14,186 bp long and contains 35 genes, including only 20 tRNAs, lacking tRNA^Gly and tRNA^Tyr. Re-annotation of the mitogenome of Steganacarus magnus increased the number of mt tRNAs for this species to 12. As typical for acariform mites, many tRNAs are highly truncated in both oribatid species. The total number of tRNAs and the number of tRNAs with a complete cloverleaf-like structure in P. leontonychus, however, clearly exceeds the numbers previously reported for Sarcoptiformes. This indicates, contrary to what has been previously assumed, that reduction of tRNAs is not a general characteristic for sarcoptiform mites. Compared to other Sarcoptiformes, the two oribatid species have the least rearranged mt genome with respect to the pattern observed in Limulus polyphemus, a basal arachnid species. Phylogenetic analysis of the newly sequenced mt genome and previously published data on other acariform mites confirms paraphyly of the Oribatida and an origin of the Astigmata within the Oribatida.

Complete mitochondrial genome of the feather mite Ardeacarus ardeae (Acari, Sarcoptiformes, Pterolichidae)

In this study, we determined the mitogenome sequence of Ardeacarus ardeae (Canestrini, 1878) in the family Pterolichidae (Acari, Sarcoptiformes), which is the first complete mitogenome sequence in feather mite. The mitogenome of A. ardeae is 14,069 bp in length and contains 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), and a control region (CR). The phylogenetic tree shows that A. ardeae belong to the supercohort Desmonomatides within the order Sarcoptiformes.

The phylogenetic position of eriophyoid mites (superfamily Eriophyoidea) in Acariformes inferred from the sequences of mitochondrial genomes and nuclear small subunit (18S) rRNA gene

Eriophyoid mites (superfamily Eriophyoidea) comprise >4400 species worldwide. Despite over a century of study, the phylogenetic position of these mites within Acariformes is still poorly resolved. Currently, Eriophyoidea is placed in the order Trombidiformes. We inferred the high-level phylogeny of Acari with the mitochondrial (mt) genome sequences of 110 species including four eriophyoid species, and the nuclear small subunit (18S) rRNA gene sequences of 226 species including 25 eriophyoid species. Maximum likelihood (ML), Bayesian inference (BI) and Maximum parsimony (MP) methods were used to analyze the sequence data. Divergence times were estimated for major lineages of Acari using Bayesian approaches. Our analyses consistently recovered the monophyly of Eriophyoidea but rejected the monophyly of Trombidiformes. The eriophyoid mites were grouped with the sarcoptiform mites, or were the sister group of sarcoptiform mites+non-eriophyoid trombidiform mites, depending on data partition strategies. Eriophyoid mites diverged from other mites in the Devonian (384Mya, 95% HPD, 352-410Mya). The origin of eriophyoid mites was dated to the Permian (262Mya, 95% HPD 230-307Mya), mostly prior to the radiation of gymnosperms (Triassic-Jurassic) and angiosperms (early Cretaceous). We propose that the placement of Eriophyoidea in the order Trombidiformes under the current classification system should be reviewed.

Suboral fork: a newly discerned gnathosomal structure from the proboscis of eriophyoid mites (Acari, Eriophyoidea)

The infracapitulum of eriophyoid mites comprises a cone-like basal infracapitulum, containing a pharynx, and a distal infracapitulum, forming a proboscis ensheathing a bunch of stylets. A well-developed basal labral section was observed in all studied specimens. A newly discerned structure, the suboral fork, situated in the ventral part of the proboscis was discovered. It is larger in diptilomiopids and Nalepella and notably smaller in eriophyids and phytoptids. This structure presumably determines the site of piercing and functions in a similar way to the pressure foot of a sewing machine which controls the movements of a needle. In diptilomiopids the suboral fork might have an additional function: it is a stopper which prevents the proboscis from further penetrating into plant tissues. It is possible that the suboral fork is homologous with the labium of early derivative acariform mites. The proboscis might be a fusion product of the infracapitular lateral lips, malapophyses and the labium. The proboscis serves as a feeding structure in eriophyoids; two ways of sucking plant cell sap, depending on shapes of proboscis and labrum, are hypothesized. Further work is needed to draw conclusions on homologies and the function of all gnathosomal structures in eriophyoids.

The complete mitochondrial genome of Histiostoma blomquisti (Acari: Histiostomatidae)

The mite Histiostoma blomquisti is a microorganism feeder that uses the red imported fire ant (Solenopsis invicta) as a phoretic carrier for dispersal. We sequenced the H. blomquisti mitogenome using next-generation sequencing methods. The circular mitogenome of H. blomquisti is 15,892 bp and is composed of 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNAs, and 6 non-coding regions >100 bp. Most tRNAs are highly reduced, like those found in other Acariformes. Phylogenetic analysis based on the concatenated nucleotide sequence of the 13 protein-coding genes supports Histiostomatid mites forming the basal-most lineage in Astigmata.