Novel Mitochondrial Gene Content and Gene Arrangement Indicate Illegitimate Inter-mtDNA Recombination in the Chigger Mite, Leptotrombidium pallidum

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.

Let me know your name: a study of chigger mites (Acariformes: Trombiculidae) associated with the edible dormouse (Glis glis) in the Carpathian-Balkan distribution gradient

Trombiculid mites were collected from the edible dormouse (Glis glis) within the Carpathian-Balkan distribution gradient of host species. Representatives of five genera (Leptotrombidium, Neotrombicula, Brunehaldia, Hirsutiella, Schoutedenichia) and 10 species of chiggers were discovered in the material, based on morphological and/or molecular data. Brunehaldia, new to the fauna of Greece, was recorded for the first time from the edible dormouse. Neotrombicula talmiensis was new to the fauna of Greece and Neotrombicula vulgaris was new to the fauna of North Macedonia. Successful amplification and sequencing of COI was carried out in relation to three genera and six species. The intraspecific variation of taxa hitherto distinguished based on morphological criteria was juxtaposed with molecular data, using the distance method and the phylogenetic approach. The molecular methods indicated wider than hitherto recognized, intraspecific morphological variation for Leptotrombidium europaeum and N. talmiensis. On the other hand, an inference limited to morphology proved to be insufficient for species delineation, which was confirmed by the relatively low identity (%) of examined COI sequences as well as the size of inter-/intraspecific K2P distance threshold. Our study provides support for integrative taxonomy that combines different sources of evidence and contributes to recognition of the scope of intraspecific variation. The high degree of hidden diversity revealed with the application of molecular tools, votes for a careful approach to the identification of chiggers. The confirmed cases of co-invasion, including the representatives of various genera (Leptotrombidium and Neotrombicula, Brunehaldia and Neotrombicula, Neotrombicula and Schoutedenichia, Hirsutiella and Schoutedenichia) additionally support the need to include all larvae found on a given host specimen in the identification process.

The complete mitochondrial genome and novel gene arrangement in Nesodiprion zhejiangensis Zhou & Xiao (Hymenoptera: Diprionidae)

The complete mitochondrial genome (mitogenome) of the sawfly, Nesodiprion zhejiangensis Zhou & Xiao, was sequenced, assembled, and deposited in GenBank (Accession Number: OM501121). The 15,660 bp N. zhejiangensis mitogenome encodes for 2 ribosomal RNAs (rrnL and rrnS), 22 transfer RNAs (tRNAs), 13 protein-coding genes (PCGs), and an AT-rich region of 450 bp in length. The nucleotide composition is biased toward adenine and thymine (A + T = 81.8%). Each PCG is initiated by an ATN codon, except for cox2, which starts with a TTG. Of 13 PCGs, 9 have a TAA termination codon, while the remainder terminate with a TAG or a single T. All tRNAs have the classic cloverleaf structure, except for the dihydrouridine (DHU) arm of tRNA^val, which forms a simple loop. There are 49 helices belonging to 6 domains in rrnL and 30 helices belonging to 4 domains in rrnS. In comparison to the ancestral architecture, N. zhejiangensis has the most rearranged mitogenome in Symphyta, in which rearrangement events of local inversion and transposition are identified in three gene clusters. Specifically, the main hotspot of gene rearrangement occurred between rrnS and trnY, and rearranged from rrnS-(AT-rich region)-I-Q-M-nd2-W-C-Y to rrnS-Q-W-C-nd2-I-M-(AT-rich region)-Y, involving a local inversion event of a large gene cluster and transposition events of some tRNAs. Transposition of trnA and trnR (rearranged from A-R to R-A) was observed at the nd3-nd5 gene junction while shuffling of trnP and trnT (rearranged from T-P to P-T) occurred at the nd4l-nd6 gene junction. While illegitimate inter-mtDNA recombination might explain the opposite orientations of transcription between rrnS and trnY, transposition events of tRNA in some gene blocks can be accounted for by the tandem duplication/random loss (TDRL) model. Our phylogenetic analysis suggests that N. zhejiangensis is closely related to congeneric species N. biremis and N. japonicus, which together form a sister lineage with the European pine sawfly, Neodiprion sertifer.

Complete mitochondrial genome of Hygrobates turcicus Pešić, Esen & Dabert, 2017 (Acari, Hydrachnidia, Hygrobatoidea)

The aim of the study was sequencing of the mitogenome of Hygrobates turcicus Pešić, Esen & Dabert, 2017 to expand knowledge of the polymorphism and cryptic or pseudocryptic diversity within Hydrachnidia. The samples originated from Bulgaria, Vidima River near Debnewo, 42°56'41.4''N, 24°48'44.6''E, depth 0.4 m, stones on the bottom, water flow 0.71 m/s, temperature 10 °C, pH 8.53, oxygen 110%, conductivity 279 µS/cm, hardness 121 CaO mg/l; 11 males, 27 females, 2 deutonymphs 12.x.2019 leg. Zawal, Michoński & Bańkowska; one male and one female dissected and slides mounted. The study was carried out using the following methods: DNA extraction, sequencing, assembly and annotation, comparison with other populations of H. turcicus, and multigene phylogeny. As a result of the study, it was determined that the mitogenome is 15,006 bp long and encodes for 13 proteins, 2 rRNAs, and 22 tRNAs. The genome is colinear with those of H. longiporus and H. taniguchii, the difference in size originating from a non-coding region located between protein-coding genes ND4L and ND3. Five genes have alternative start-codon, and four display premature termination. The multigene phylogeny obtained using all mitochondrial protein-coding genes unambiguously associates H. turcicus with the cluster formed by H. longiporus and H. taniguchii.

Biology, Systematics, Microbiome, Pathogen Transmission and Control of Chiggers (Acari: Trombiculidae, Leeuwenhoekiidae) with Emphasis on the United States

Chiggers are the larval stage of Trombiculidae and Leeuwenhoekiidae mites of medical and veterinary importance. Some species in the genus Leptotrombidium and Herpetacarus vector Orientia species, the bacteria that causes scrub typhus disease in humans. Scrub typhus is a life-threatening, febrile disease. Chigger bites can also cause dermatitis. There were 248 chigger species reported from the US from almost every state. However, there are large gaps in our knowledge of the life history of other stages of development. North American wide morphological keys are needed for better species identification, and molecular sequence data for identification are minimal and not clearly matched with morphological data. The role of chiggers in disease transmission in the US is especially understudied, and the role of endosymbionts in Orientia infection are suggested in the scientific literature but not confirmed. The most common chiggers in the eastern United States were identified as Eutrombicula alfreddugesi but were likely misidentified and should be replaced with Eutrombicula cinnabaris. Scrub typhus was originally believed to be limited to the Tsutsugamushi Triangle and the chigger genus, Leptotrombidium, but there is increasing evidence this is not the case. The potential of Orientia species establishing in the US is high. In addition, several other recognized pathogens to infect humans, namely Hantavirus, Bartonella, Borrelia, and Rickettsia, were also detected in chiggers. The role that chiggers play in these disease transmissions in the US needs further investigation. It is possible some of the tick-borne diseases and red meat allergies are caused by chiggers.

Genome Organization and Comparative Evolutionary Mitochondriomics of Brown Planthopper, Nilaparvata lugens Biotype 4 Using Next Generation Sequencing (NGS)

Nilaparvata lugens is the main rice pest in India. Until now, the Indian N. lugens mitochondrial genome has not been sequenced, which is a very important basis for population genetics and phylogenetic evolution studies. An attempt was made to sequence two examples of the whole mitochondrial genome of N. lugens biotype 4 from the Indian population for the first time. The mitogenomes of N. lugens are 16,072 and 16,081 bp long with 77.50% and 77.45% A + T contents, respectively, for both of the samples. The mitochondrial genome of N. lugens contains 37 genes, including 13 protein-coding genes (PCGs) (cox1-3, atp6, atp8, nad1-6, nad4l, and cob), 22 transfer RNA genes, and two ribosomal RNA (rrnS and rrnL) subunits genes, which are typical of metazoan mitogenomes. However, both samples of N. lugens mitogenome in the present study retained one extra copy of the trnC gene. Additionally, we also found 93 bp lengths for the atp8 gene in both of the samples, which were 60–70 bp less than that of the other sequenced mitogenomes of hemipteran insects. The phylogenetic analysis of the 19 delphacids mitogenome dataset yielded two identical topologies when rooted with Ugyops sp. in one clade, and the remaining species formed another clade with P. maidis and M. muiri being sisters to the remaining species. Further, the genus Nilaparvata formed a separate subclade with the other genera (Sogatella, Laodelphax, Changeondelphax, and Unkanodes) of Delphacidae. Additionally, the relationship among the biotypes of N. lugens was recovered as the present study samples (biotype-4) were separated from the three biotypes reported earlier. The present study provides the reference mitogenome for N. lugens biotype 4 that may be utilized for biotype differentiation and molecular-aspect-based future studies of N. lugens.

Functional Characterization of Non-Ankyrin Repeat Domains of Orientia tsutsugamushi Ank Effectors Reveals Their Importance for Molecular Pathogenesis

Orientia tsutsugamushi is a genetically intractable obligate intracellular bacterium, causes scrub typhus, and has one of the largest known armamentariums of ankyrin repeat-containing effectors (Anks). Most have a C-terminal F-box presumed to interact with the SCF ubiquitin ligase complex primarily based on their ability to bind overexpressed Skp1. Whether all F-box-containing Anks bind endogenous SCF components and the F-box residues essential for such interactions has gone unexplored. Many O. tsutsugamushi Ank F-boxes occur as part of a PRANC (pox protein repeats of ankyrin-C-terminal) domain. Roles of the non-F-box portion of the PRANC and intervening sequence region (ISR) that links the ankyrin repeat and F-box/PRANC domains are unknown. The functional relevance of these effectors' non-ankyrin repeat domains was investigated. The F-box was necessary for Flag-tagged versions of most F-box-containing Anks to precipitate endogenous Skp1, Cul1, and/or Rbx1, while the ISR and PRANC were dispensable. Ank toxicity in yeast was predominantly F-box dependent. Interrogations of Ank1, Ank5, and Ank6 established that L1, P2, E4, I9, and D17 of the F-box consensus are key for binding native SCF components and for Ank1 and Ank6 to inhibit NF-κB. The ISR is also essential for Ank1 and Ank6 to impair NF-κB. Ectopically expressed Ank1 and Ank6 lacking the ISR or having a mutagenized F-box incapable of binding SCF components performed as dominant-negative inhibitors to block O. tsutsugamushi NF-κB modulation. This study advances knowledge of O. tsutsugamushi Ank functional domains and offers an approach for validating their roles in infection.

Parasitism, seasonality, and diversity of trombiculid mites (Trombidiformes: Parasitengona, Trombiculidae) infesting bats (Chiroptera) in Poland

The study aims to ascertain the diversity of trombiculid species associated with Chiroptera in Poland, and for the first time in the case of research on Central European Trombiculidae, we use both DNA and morphology in an integrative taxonomic approach to determine species identities of trombiculids. The research was carried out from 2015 to 2019. In total, 2725 larvae were collected from 300 specimens of bats belonging to 11 species. Deutonymphs were obtained through laboratory rearing of larvae; few larvae and deutonymphs were collected also from bats' daily roosts. The presence of trombiculid larvae on hosts was observed between July and April of the following year, with the highest numbers recorded in autumn, during bat swarming. Male bats were infested more often than females (16.4 vs. 6.6%). The highest infestation rate was recorded for Barbastella barbastellus, Myotis nattereri and Plecotus auritus, and the highest prevalence of chiggers (> 30%) for Myotis bechsteinii and P. auritus. The larvae found on bats occupied the areas with free access to the host's skin: auricles, tragus, and snout. Morphological identification of specimens to the species level was hindered by the mosaic distribution of diagnostic traits. Morphological analyses indicated the presence of Leptotrombidium russicum and Leptotrombidium spp. in the examined material, whereas molecular analyses additionally suggested three other potential species assigned to the same genus based on the assessed scope of intrageneric variation (ASAP method). We argue that the identification of the parasitic larvae (chiggers) using morphological characters does not address the question of actual species boundaries, which, in turn, affects the inferences about host specificity and host range.

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.

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.

The challenge of Coccidae (Hemiptera: Coccoidea) mitochondrial genomes: The case of Saissetia coffeae with novel truncated tRNAs and gene rearrangements

There have been few reports of complete mitochondrial genomes (mitogenomes) of scale insects, and it has been indicated that complex and novel structures in their mitogenomes may lead to difficulties in sequencing, assembly and annotation. Transfer RNAs (tRNAs) usually possess typical cloverleaf secondary structures, and truncated tRNAs are rarely found in insect mitogenomes. Here, we report a complete Saissetia coffeae mitogenome (15,389 bp) with high A+T content (84.7%) sequenced by next-generation sequencing (NGS) methods. Genes in the mitogenome were annotated, and nine tRNAs were not found using MITOS. Most of the detected tRNAs were significantly truncated without the dihydrouridine (DHU) arm or the TΨC (T) arm. In addition, the 9 "lost" tRNAs containing mismatched base pairs were retrieved based on the tRNA annotation workflow for Coccidae described in our study. The gene arrangement in the Saissetia coffeae mitogenome was significantly different from that in other hemipteran insects. Additionally, Bayesian and maximum likelihood trees based on the mitochondrial genes showed a long branch of the Saissetia lineage, indicating significant nonsynonymous substitutions or high evolutionary rates in the Saissetia lineage. We provide a reference mitogenome for the assembly and annotation of the Coccidae mitogenome and offer insights into the evolution of scale insects.

Double control regions of some flatfish mitogenomes evolve in a concerted manner

Mitochondrial genomes (mitogenomes) typically contain 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and a single control region (CR). Flatfish mitochondrial genomes (mitogenomes) from three different genera in Bothidae (bothids) contain double CRs that evolved in a concerted manner. How these double CRs maintained identical sequences throughout the evolutionary process is an interesting issue. In the present study, over four hundred arrays of the double CRs of mitogenomes from three bothids (Arnoglossus tenuis, Lophonectes gallus and Psettina iijimae) were performed. Interesting variations between double CRs were observed in P. iijimae mitogenomes, and the networks of CR sequences from P. iijimae indicated a high possibility of genetic information exchange between CRs. No recombination product was detected in our results, indicating that the mechanism of the concerted evolution between the double CRs of P. iijimae was not recombination. We speculate that mismatch repair, a mitochondrial DNA repair mechanism, is a potential explanation for the concerted evolution between these double CRs.

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.

The complete mitochondrial genome of Clostera anastomosis (Lepidoptera: Notodontidae) and implication for the phylogenetic relationships of Noctuoidea species

In the present study, the complete mitochondrial genome (mitogenome) of Clostera anastomosis (C. anastomosis) has been determined for the first time. The mitogenome is 15,390 base pairs (bp) in length, comprised of 13 protein-coding genes (PCGs), 2 ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs) and one non-coding control region (CR). The gene order shows a typical trnM rearrangement (trnM-trnI-trnQ) compared to ancestral insects (trnI-trnQ-trnM). Almost all the PCGs have the same start codon (ATN) except for cox1 (CGA), and almost all tRNAs have a typical cloverleaf secondary structure except for trnS1. At the beginning of the CR, we found a conserved motif "ATAGA + poly-T" as found in other lepidopteran insects. There are 20 intergenic regions and 11 overlapping regions, ranging from 1 to 53 bp and 1 to 9 bp, respectively. The A + T content is relatively high across the whole mitogenome. The optimal tree topologies of Noctuoidea were given by the dataset consisting of all 13 PCGs from five families (exclude Oenosandridae). Our trees suggested a topology of (Notodontidae + (Erebidae + (Nolidae + (Euteliidae + Noctuidae)))) and identified that C. anastomosis belongs to Notodontidae.

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.

Next-Generation Sequencing of Two Mitochondrial Genomes from Family Pompilidae (Hymenoptera: Vespoidea) Reveal Novel Patterns of Gene Arrangement

Animal mitochondrial genomes have provided large and diverse datasets for evolutionary studies. Here, the first two representative mitochondrial genomes from the family Pompilidae (Hymenoptera: Vespoidea) were determined using next-generation sequencing. The sequenced region of these two mitochondrial genomes from the species Auplopus sp. and Agenioideus sp. was 16,746 bp long with an A + T content of 83.12% and 16,596 bp long with an A + T content of 78.64%, respectively. In both species, all of the 37 typical mitochondrial genes were determined. The secondary structure of tRNA genes and rRNA genes were predicted and compared with those of other insects. Atypical trnS1 using abnormal anticodons TCT and lacking D-stem pairings was identified. There were 49 helices belonging to six domains in rrnL and 30 helices belonging to three domains in rrns present. Compared with the ancestral organization, four and two tRNA genes were rearranged in mitochondrial genomes of Auplopus and Agenioideus, respectively. In both species, trnM was shuffled upstream of the trnI-trnQ-trnM cluster, and trnA was translocated from the cluster trnA-trnR-trnN-trnS1-trnE-trnF to the region between nad1 and trnL1, which is novel to the Vespoidea. In Auplopus, the tRNA cluster trnW-trnC-trnY was shuffled to trnW-trnY-trnC. Phylogenetic analysis within Vespoidea revealed that Pompilidae and Mutillidae formed a sister lineage, and then sistered Formicidae. The genomes presented in this study have enriched the knowledge base of molecular markers, which is valuable in respect to studies about the gene rearrangement mechanism, genomic evolutionary processes and phylogeny of Hymenoptera.

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.

Characterization of the complete mitochondrial genome of the storage mite pest Tyrophagus longior (Gervais) (Acari: Acaridae) and comparative mitogenomic analysis of four acarid mites

Mites of the genus Tyrophagus are economically important polyphagous pest commonly living on stored products and also responsible for allergic reactions to humans. Complete mitochondrial genomes (mitogenomes) and the gene features therein are widely used as molecular markers in the study of population genetics, phylogenetics as well as molecular evolution. However, scarcity on the sequence data has greatly impeded the studies in these areas pertaining to the Acari (mites and ticks). Information on the Tyrophagus mitogenomes is quite critical for phylogenetic evaluation and molecular evolution of the mitogenomes within Acariformes. Herein, we reported the complete mitogenome of the allergenic acarid storage mite Tyrophagus longior (Astigmata: Acaridae), an important member of stored food pests, and compared with those of other three acarid mites. The complete mitogenome of T. longior was a circular molecule of 13,271 bp. Unexpectedly, only 19 transfer RNA genes (tRNAs) were present, lacking trnF, trnS1 and trnQ. Furthermore, it also contained 13 protein-coding genes (PCGs) and 2 genes for rRNA (rrnS and rrnL) commonly detected in metazoans. The four mitogenomes displayed similar characteristics with respect to the gene content, nucleotide comparison, and codon usages. Yet, the gene order of T. longior was different from that in other Acari. The J-strands of the four mitogenomes possessed high A+T content (67.4-70.0%), and exhibited positive GC-skews and negative AT-skews. Most inferred tRNAs of T. longior were extremely truncated, lacking either a D- or T-arm, as found in other acarid mites. In T. longior mitogenome the A+T-rich region was just 50 bp in length and can be folded as a stable stem-loop structure, whereas in the region some structures of microsatellite-like (AT)n and palindromic sequences was not present. Besides, reconstructing of the phylogenetic relationship based on concatenated amino acid sequences of 13 PCGs supported that monophyly of the family Acaridae and the order Astigmata, to which the former belongs. Our results were consistent with the traditional classifications.

The complete mitogenomes of Calameuta filiformis (Eversmann, 1847) and Calameuta idolon (Rossi, 1794) (Hymenoptera: Cephidae): The remarkable features of the elongated A+T rich region in Cephini

Two complete mitogenomes of the stem borers, Calameuta filiformis and Calameuta idolon, and the complete A+T-rich region of Trachelus iudaicus (Hymenoptera: Cephidae), are reported. The mitogenomes of these species are the longest reported from hymenopterans to date. A remarkable increase in length of the A+T-rich region, the longest for Hymenoptera, was found and compared across the tribe Cephini. The presence of the tRNA- and rRNA-like sequences were reported in the A+T-rich region of sawflies and they were suggested to play a role in replication and/or transcription. The long and short tandem repeats were orderly located in both sides of the A+T-rich region producing stable secondary structures. We suggest that the short tandem repeats are likely to function as a replication fork barrier.

Extreme variation in patterns of tandem repeats in mitochondrial control region of yellow-browed tits (Sylviparus modestus, Paridae)

To investigate the evolutionary pattern and origins of tandem repeats in the mitochondrial control region of the yellow-browed tit (Sylviparus modestus), the control region and another four mitochondrial loci from fifteen individuals were analyzed. A 117-bp tandem repeat unit that repeated once, twice or three times in different individuals was found, and a rarely reported arrangement for this tandem repeats region that a 5' imperfect copy at its downstream and a 3' imperfect copy at its upstream was observed. The haplotype network, phylogenetic trees, and ancestral state reconstruction of the combined dataset of five loci suggested multiple origins of the same repeat number. The turnover model via slipped-strand mispairing was introduced to interpret the results, because mispairing occurred so frequently that multiple origins of certain repeat number were observed. Insertion via recombination should be a better explanation for the origin of this tandem repeat unit, considering characteristics of the combined sequence of the 3' and 5' imperfect copy, including identification of its homolog in other passerines and its predicted secondary structure.

Phylogenetic utility, and variability in structure and content, of complete mitochondrial genomes among genetic lineages of the Hawaiian anchialine shrimp Halocaridina rubra Holthuis 1963 (Atyidae:Decapoda)

The Atyidae are caridean shrimp possessing hair-like setae on their claws and are important contributors to ecological services in tropical and temperate fresh and brackish water ecosystems. Complete mitochondrial genomes have only been reported from five of the 449 species in the family, thus limiting understanding of mitochondrial genome evolution and the phylogenetic utility of complete mitochondrial sequences in the Atyidae. Here, comparative analyses of complete mitochondrial genomes from eight genetic lineages of Halocaridina rubra, an atyid endemic to the anchialine ecosystem of the Hawaiian Archipelago, are presented. Although gene number, order, and orientation were syntenic among genomes, three regions were identified and further quantified where conservation was substantially lower: (1) high length and sequence variability in the tRNA-Lys and tRNA-Asp intergenic region; (2) a 317-bp insertion between the NAD6 and CytB genes confined to a single lineage and representing a partial duplication of CytB; and (3) the putative control region. Phylogenetic analyses utilizing complete mitochondrial sequences provided new insights into relationships among the H. rubra genetic lineages, with the topology of one clade correlating to the geologic sequence of the islands. However, deeper nodes in the phylogeny lacked bootstrap support. Overall, our results from H. rubra suggest intra-specific mitochondrial genomic diversity could be underestimated across the Metazoa since the vast majority of complete genomes are from just a single individual of a species.

The complete mitochondrial genomes of six species of Tetranychus provide insights into the phylogeny and evolution of spider mites

Many spider mites belonging to the genus Tetranychus are of agronomical importance. With limited morphological characters, Tetranychus mites are usually identified by a combination of morphological characteristics and molecular diagnostics. To clarify their molecular evolution and phylogeny, the mitochondrial genomes of the green and red forms of Tetranychus urticae as well as T. kanzawai, T. ludeni, T. malaysiensis, T. phaselus, T. pueraricola were sequenced and compared. The seven mitochondrial genomes are typical circular molecules of about 13,000 bp encoding and they are composed of the complete set of 37 genes that are usually found in metazoans. The order of the mitochondrial (mt) genes is the same as that in the mt genomes of Panonychus citri and P. ulmi, but very different from that in other Acari. The J-strands of the mitochondrial genomes have high (∼ 84%) A+T contents, negative GC-skews and positive AT-skews. The nucleotide sequence of the cox1 gene, which is commonly used as a taxon barcode and molecular marker, is more highly conserved than the nucleotide sequences of other mitochondrial genes in these seven species. Most tRNA genes in the seven genomes lose the D-arm and/or the T-arm. The functions of these tRNAs need to be evaluated. The mitochondrial genome of T. malaysiensis differs from the other six genomes in having a slightly smaller genome size, a slight difference in codon usage, and a variable loop in place of the T-arm of some tRNAs by a variable loop. A phylogenic analysis shows that T. malaysiensis first split from other Tetranychus species and that the clade of the family Tetranychoidea occupies a basal position in the Trombidiformes. The mt genomes of the green and red forms of T. urticae have limited divergence and short evolutionary distance.

Evolution of multipartite mitochondrial genomes in the booklice of the genus Liposcelis (Psocoptera)

BACKGROUND

The genus Liposcelis (Psocoptera: Troctomorpha) has more than 120 species with a worldwide distribution and they pose a risk for global food security. The organization of mitochondrial (mt) genomes varies between the two species of booklice investigated in the genus Liposcelis. Liposcelis decolor has its mt genes on a single chromosome, like most other insects; L. bostrychophila, however, has a multipartite mt genome with genes on two chromosomes.

RESULTS

To understand how multipartite mt genome organization evolved in the genus Liposcelis, we sequenced the mt genomes of L. entomophila and L. paeta in this study. We found that these two species of booklice also have multipartite mt genomes, like L. bostrychophila, with the mt genes we identified on two chromosomes. Numerous pseudo mt genes and non-coding regions were found in the mt genomes of these two booklice, and account for 30% and 10% respectively of the entire length we sequenced. In L. bostrychophila, the mt genes are distributed approximately equally between the two chromosomes. In L. entomophila and L. paeta, however, one mt chromosome has most of the genes we identified whereas the other chromosome has largely pseudogenes and non-coding regions. L. entomophila and L. paeta differ substantially from each other and from L. bostrychophila in gene content and gene arrangement in their mt chromosomes.

CONCLUSIONS

Our results indicate unusually fast evolution in mt genome organization in the booklice of the genus Liposcelis, and reveal different patterns of mt genome fragmentation among L. bostrychophila, L. entomophila and L. paeta.

The complete mitochondrial genome of the brown leg mite, Aleuroglyphus ovatus (Acari: Sarcoptiformes): evaluation of largest non-coding region and unique tRNAs

The circular mitochondrial genome (mitogenome) of Aleuroglyphus ovatus was sequenced. It was 14,328 bp long, and consisted of 37 coding genes including 13 protein-coding genes, two ribosomal RNA genes, and 22 transfer RNA genes. This is the first description of the complete mitogenome of a species in the Acaridae (Acari: Sarcoptiformes). The mtDNA gene order for A. ovatus is identical to those of Dermatophagoides farinae and D. pteronyssinus, but distinctly different from the mtDNA of other Acari. Most inferred tRNA genes of A. ovatus are extremely truncated (48-62 bp), lack stem-loops on either the T- or D-arm (except the trnK), and are unable to fold into the canonical tRNA cloverleaf structure. The largest non-coding region (378 bp) contained several conserved sequences involved in the regulation of mitogenome replication, including one core sequence (ACAT) associated with termination of the J-strand replication and several hypothetical stem-loop structures. The microsatellite-like (AT)n sequence in the largest non-coding region was observed in two other Astigmata species, but it has not been found in other Acari.

The complete mitochondrial genomes of two rice planthoppers, Nilaparvata lugens and Laodelphax striatellus: conserved genome rearrangement in Delphacidae and discovery of new characteristics of atp8 and tRNA genes

Background

Nilaparvata lugens (the brown planthopper, BPH) and Laodelphax striatellus (the small brown planthopper, SBPH) are two of the most important pests of rice. Up to now, there was only one mitochondrial genome of rice planthopper has been sequenced and very few dependable information of mitochondria could be used for research on population genetics, phylogeographics and phylogenetic evolution of these pests. To get more valuable information from the mitochondria, we sequenced the complete mitochondrial genomes of BPH and SBPH. These two planthoppers were infected with two different functional Wolbachia (intracellular endosymbiont) strains (wLug and wStri). Since both mitochondria and Wolbachia are transmitted by cytoplasmic inheritance and it was difficult to separate them when purified the Wolbachia particles, concomitantly sequencing the genome of Wolbachia using next generation sequencing method, we also got nearly complete mitochondrial genome sequences of these two rice planthoppers. After gap closing, we present high quality and reliable complete mitochondrial genomes of these two planthoppers.

Results

The mitogenomes of N. lugens (BPH) and L. striatellus (SBPH) are 17, 619 bp and 16, 431 bp long with A + T contents of 76.95% and 77.17%, respectively. Both species have typical circular mitochondrial genomes that encode the complete set of 37 genes which are usually found in metazoans. However, the BPH mitogenome also possesses two additional copies of the trnC gene. In both mitochondrial genomes, the lengths of the atp8 gene were conspicuously shorter than that of all other known insect mitochondrial genomes (99 bp for BPH, 102 bp for SBPH). That two rearrangement regions (trnC-trnW and nad6-trnP-trnT) of mitochondrial genomes differing from other known insect were found in these two distantly related planthoppers revealed that the gene order of mitochondria might be conservative in Delphacidae. The large non-coding fragment (the A+T-rich region) putatively corresponding responsible for the control of replication and transcription of mitochondria contained a variable number of tandem repeats (VNTRs) block in different natural individuals of these two planthoppers. Comparison with a previously sequenced individual of SBPH revealed that the mitochondrial genetic variation within a species exists not only in the sequence and secondary structure of genes, but also in the gene order (the different location of trnH gene).

Conclusion

The mitochondrial genome arrangement pattern found in planthoppers was involved in rearrangements of both tRNA genes and protein-coding genes (PCGs). Different species from different genera of Delphacidae possessing the same mitochondrial gene rearrangement suggests that gene rearrangements of mitochondrial genome probably occurred before the differentiation of this family. After comparatively analyzing the gene order of different species of Hemiptera, we propose that except for some specific taxonomical group (e.g. the whiteflies) the gene order might have diversified in family level of this order. The VNTRs detected in the control region might provide additional genetic markers for studying population genetics, individual difference and phylogeographics of planthoppers.

Evolution of extensively fragmented mitochondrial genomes in the lice of humans

Bilateral animals are featured by an extremely compact mitochondrial (mt) genome with 37 genes on a single circular chromosome. The human body louse, Pediculus humanus, however, has its mt genes on 20 minichromosomes. We sequenced the mt genomes of two other human lice: the head louse, P. capitis, and the pubic louse, Pthirus pubis. Comparison among the three human lice revealed the presence of fragmented mt genomes in their most recent common ancestor, which lived ∼7 Ma. The head louse has exactly the same set of mt minichromosomes as the body louse, indicating that the number of minichromosomes, and the gene content and gene arrangement in each minichromosome have remained unchanged since the body louse evolved from the head louse ∼107,000 years ago. The pubic louse has the same pattern of one protein-coding or rRNA gene per minichromosome (except one minichromosome with two protein-coding genes, atp6 and atp8) as the head louse and the body louse. This pattern is apparently ancestral to all human lice and has been stable for at least 7 Myr. Most tRNA genes of the pubic louse, however, are on different minichromosomes when compared with their counterparts in the head louse and the body louse. It is evident that rearrangement of four tRNA genes (for leucine, arginine and glycine) was due to gene-identity switch by point mutation at the third anticodon position or by homologous recombination, whereas rearrangement of other tRNA genes was by gene translocation between minichromosomes, likely caused by minichromosome split via gene degeneration and deletion.

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

According to the tRNA punctuation model, the mitochondrial genome (mtDNA) of mammals and arthropods is transcribed as large polycistronic precursors that are maturated by endonucleolytic cleavage at tRNA borders and RNA polyadenylation. Starting from the newly sequenced mtDNA of Ixodes ricinus and using a combination of mitogenomics and transcriptional analyses, we found that in all currently-sequenced tick lineages (Prostriata, Metastriata and Argasidae) the 3'-end of the polyadenylated nad1 and rrnL transcripts does not follow the tRNA punctuation model and is located upstream of a degenerate 17-bp DNA motif. A slightly different motif is also present downstream the 3'-end of nad1 transcripts in the primitive chelicerate Limulus polyphemus and in Drosophila species, indicating the ancient origin and the evolutionary conservation of this motif in arthropods. The transcriptional analyses suggest that this motif directs the 3'-end formation of the nad1/rrnL mature RNAs, likely working as a transcription termination signal or a processing signal of precursor transcripts. Moreover, as most regulatory elements, this motif is characterized by a taxon-specific evolution. Although this signal is not exclusive of ticks, making a play on words it has been named "Tick-Box", since it is a check mark that has to be verified for the 3'-end formation of some mt transcripts, and its consensus sequence has been here carefully characterized in ticks. Indeed, in the whole mtDNA of all ticks, the Tick-Box is always present downstream of nad1 and rrnL, mainly in non-coding regions (NCRs) and occasionally within trnL(CUN). However, some metastriates present a third Tick-Box at an intriguing site--inside the small NCR located at one end of a 3.4 kb translocated region, the other end of which exhibits the nad1 Tick-Box--hinting that this motif could have been involved in metastriate gene order rearrangements.

Pseudoscorpion mitochondria show rearranged genes and genome-wide reductions of RNA gene sizes and inferred structures, yet typical nucleotide composition bias

BACKGROUND

Pseudoscorpions are chelicerates and have historically been viewed as being most closely related to solifuges, harvestmen, and scorpions. No mitochondrial genomes of pseudoscorpions have been published, but the mitochondrial genomes of some lineages of Chelicerata possess unusual features, including short rRNA genes and tRNA genes that lack sequence to encode arms of the canonical cloverleaf-shaped tRNA. Additionally, some chelicerates possess an atypical guanine-thymine nucleotide bias on the major coding strand of their mitochondrial genomes.

RESULTS

We sequenced the mitochondrial genomes of two divergent taxa from the chelicerate order Pseudoscorpiones. We find that these genomes possess unusually short tRNA genes that do not encode cloverleaf-shaped tRNA structures. Indeed, in one genome, all 22 tRNA genes lack sequence to encode canonical cloverleaf structures. We also find that the large ribosomal RNA genes are substantially shorter than those of most arthropods. We inferred secondary structures of the LSU rRNAs from both pseudoscorpions, and find that they have lost multiple helices. Based on comparisons with the crystal structure of the bacterial ribosome, two of these helices were likely contact points with tRNA T-arms or D-arms as they pass through the ribosome during protein synthesis.The mitochondrial gene arrangements of both pseudoscorpions differ from the ancestral chelicerate gene arrangement. One genome is rearranged with respect to the location of protein-coding genes, the small rRNA gene, and at least 8 tRNA genes. The other genome contains 6 tRNA genes in novel locations. Most chelicerates with rearranged mitochondrial genes show a genome-wide reversal of the CA nucleotide bias typical for arthropods on their major coding strand, and instead possess a GT bias. Yet despite their extensive rearrangement, these pseudoscorpion mitochondrial genomes possess a CA bias on the major coding strand. Phylogenetic analyses of all 13 mitochondrial protein-coding gene sequences consistently yield trees that place pseudoscorpions as sister to acariform mites.

CONCLUSION

The well-supported phylogenetic placement of pseudoscorpions as sister to Acariformes differs from some previous analyses based on morphology. However, these two lineages share multiple molecular evolutionary traits, including substantial mitochondrial genome rearrangements, extensive nucleotide substitution, and loss of helices in their inferred tRNA and rRNA structures.

Mitochondrial genome sequence of Unionicola parkeri (Acari: Trombidiformes: Unionicolidae): molecular synapomorphies between closely-related Unionicola gill mites

The mitochondrial genome of Unionicola parkeri is a 14,734 bp circular DNA molecule. The sequence and annotation revealed a unique gene order, related to but distinct from the gene order in the closely related species U. foili. Mitochondrial tRNA sequences annotated in this genome predict non-canonical secondary structures for these molecules. The continuing patterns of unique gene orders and unusual tRNA structures in the Trombidiformes in general and Unionicola in particular support the use of phylogenetic approaches that use these types of molecular features as shared, derived character states. Further progress in using these molecular character states to reconstruct phylogeny will depend on careful annotation, especially of tRNA genes.

The complete mitochondrial genome of the citrus red mite Panonychus citri (Acari: Tetranychidae): high genome rearrangement and extremely truncated tRNAs

Background

The family Tetranychidae (Chelicerata: Acari) includes ~1200 species, many of which are of agronomic importance. To date, mitochondrial genomes of only two Tetranychidae species have been sequenced, and it has been found that these two mitochondrial genomes are characterized by many unusual features in genome organization and structure such as gene order and nucleotide frequency. The scarcity of available sequence data has greatly impeded evolutionary studies in Acari (mites and ticks). Information on Tetranychidae mitochondrial genomes is quite important for phylogenetic evaluation and population genetics, as well as the molecular evolution of functional genes such as acaricide-resistance genes. In this study, we sequenced the complete mitochondrial genome of Panonychus citri (Family Tetranychidae), a worldwide citrus pest, and provide a comparison to other Acari.

Results

The mitochondrial genome of P. citri is a typical circular molecule of 13,077 bp, and contains the complete set of 37 genes that are usually found in metazoans. This is the smallest mitochondrial genome within all sequenced Acari and other Chelicerata, primarily due to the significant size reduction of protein coding genes (PCGs), a large rRNA gene, and the A + T-rich region. The mitochondrial gene order for P. citri is the same as those for P. ulmi and Tetranychus urticae, but distinctly different from other Acari by a series of gene translocations and/or inversions. The majority of the P. citri mitochondrial genome has a high A + T content (85.28%), which is also reflected by AT-rich codons being used more frequently, but exhibits a positive GC-skew (0.03). The Acari mitochondrial nad1 exhibits a faster amino acid substitution rate than other genes, and the variation of nucleotide substitution patterns of PCGs is significantly correlated with the G + C content. Most tRNA genes of P. citri are extremely truncated and atypical (44-65, 54.1 ± 4.1 bp), lacking either the T- or D-arm, as found in P. ulmi, T. urticae, and other Acariform mites.

Conclusions

The P. citri mitochondrial gene order is markedly different from those of other chelicerates, but is conserved within the family Tetranychidae indicating that high rearrangements have occurred after Tetranychidae diverged from other Acari. Comparative analyses suggest that the genome size, gene order, gene content, codon usage, and base composition are strongly variable among Acari mitochondrial genomes. While extremely small and unusual tRNA genes seem to be common for Acariform mites, further experimental evidence is needed.

Mosaic gene conversion after a tandem duplication of mtDNA sequence in Diomedeidae (albatrosses)

Although the tandem duplication of mitochondrial (mt) sequences, especially those of the control region (CR), has been detected in metazoan species, few studies have focused on the features of the duplicated sequence itself, such as the gene conversion rate, distribution patterns of the variation, and relative rates of evolution between the copies. To investigate the features of duplicated mt sequences, we partially sequenced the mt genome of 16 Phoebastria albatrosses belonging to three species (P. albatrus, P. nigripes, and P. immutabilis). More than 2,300 base pairs of tandemly-duplicated sequence were shared by all three species. The observed gene arrangement was shared in the three Phoebastria albatrosses and suggests that the duplication event occurred in the common ancestor of the three species. Most of the copies in each individual were identical or nearly identical, and were maintained through frequent gene conversions. By contrast, portions of CR domains I and III had different phylogenetic signals, suggesting that gene conversion had not occurred in those sections after the speciation of the three species. Several lines of data, including the heterogeneity of the rate of molecular evolution, nucleotide differences, and putative secondary structures, suggests that the two sequences in CR domain I are maintained through selection; however, additional studies into the mechanisms of gene conversion and mtDNA synthesis are required to confirm this hypothesis.

Mitochondrial genome analysis of the predatory mite Phytoseiulus persimilis and a revisit of the Metaseiulus occidentalis mitochondrial genome

In this study we sequenced and analysed the complete mitochondrial (mt) genome of the Chilean predatory mite Phytoseiulus persimilis Athias-Henriot (Chelicerata: Acari: Mesostigmata: Phytoseiidae: Amblyseiinae). The 16 199 bp genome (79.8% AT) contains the standard set of 13 protein-coding and 24 RNA genes. Compared with the ancestral arthropod mtDNA pattern, the gene order is extremely reshuffled (35 genes changed position) and represents a novel arrangement within the arthropods. This is probably related to the presence of several large noncoding regions in the genome. In contrast with the mt genome of the closely related species Metaseiulus occidentalis (Phytoseiidae: Typhlodrominae) — which was reported to be unusually large (24 961 bp), to lack nad6 and nad3 protein-coding genes, and to contain 22 tRNAs without T-arms — the genome of P. persimilis has all the features of a standard metazoan mt genome. Consequently, we performed additional experiments on the M. occidentalis mt genome. Our preliminary restriction digests and Southern hybridization data revealed that this genome is smaller than previously reported. In addition, we cloned nad3 in M. occidentalis and positioned this gene between nad4L and 12S-rRNA on the mt genome. Finally, we report that at least 15 of the 22 tRNAs in the M. occidentalis mt genome can be folded into canonical cloverleaf structures similar to their counterparts in P. persimilis.

Improved tRNA prediction in the American house dust mite reveals widespread occurrence of extremely short minimal tRNAs in acariform mites

BACKGROUND

Atypical tRNAs are functional minimal tRNAs, lacking either the D- or T-arm. They are significantly shorter than typical cloverleaf tRNAs. Widespread occurrence of atypical tRNAs was first demonstrated for secernentean nematodes and later in various arachnids. Evidence started to accumulate that tRNAs of certain acariform mites are even shorter than the minimal tRNAs of nematodes, raising the possibility that tRNAs lacking both D- and T-arms might exist in these organisms. The presence of cloverleaf tRNAs in acariform mites, particularly in the house dust mite genus Dermatophagoides, is still disputed.

RESULTS

Mitochondrial tRNAs of Dermatophagoides farinae are minimal, atypical tRNAs lacking either the T- or D-arm. The size (49-62, 54.4 +/- 2.86 nt) is significantly (p = 0.019) smaller than in Caenorhabditis elegans (53-63, 56.3 +/- 2.30 nt), a model minimal tRNA taxon. The shortest tRNA (49 nt) in Dermatophagoides is approaching the length of the shortest known tRNAs (45-49 nt) described in other acariform mites. The D-arm is absent in these tRNAs, and the inferred T-stem is small (2-3 bp) and thermodynamically unstable, suggesting that it may not exist in reality. The discriminator nucleotide is probably not encoded and is added postranscriptionally in many Dermatophagoides tRNAs.

CONCLUSIONS

Mitochondrial tRNAs of acariform mites are largely atypical, non-cloverleaf tRNAs. Among them, the shortest known tRNAs with no D-arm and a short and unstable T-arm can be inferred. While our study confirmed seven tRNAs in Dermatophagoides by limited EST data, further experimental evidence is needed to demonstrate extremely small and unusual tRNAs in acariform mites.

Mitochondrial genome sequence of Unionicola foili (Acari: Unionicolidae): a unique gene order with implications for phylogenetic inference

The mitochondrial genome of Unionicola foili is circular, 14,738 bp in length, and contains several notable features. The sequence and annotation revealed a unique gene order, continuing a pattern of highly rearranged mitochondrial genomes in the Trombidiformes. U. foili mitochondrial tRNA sequences predict non-canonical secondary structures for these molecules, and our annotation suggests an in-frame fusion between the nad4L and nad5 genes in this genome. The unique gene order and unusual tRNA structures could serve as idiosyncratic characters and have the potential to be phylogenetically informative.

Mitochondrial rRNA secondary structures and genome arrangements distinguish chelicerates: comparisons with a harvestman (Arachnida: Opiliones: Phalangium opilio)

Arachnids are a highly diverse group of arthropods, and many of the mitochondrial genomes that have been sequenced from arachnids possess unusual features in their inferred gene structures and genome organization. The first complete sequence of a mitochondrial genome from the arachnid order Opiliones (harvestmen) is presented here. Secondary structures of the two mitochondrial ribosomal subunits of Phalangium opilio are inferred and compared to mitochondrial rRNA structures of a hexapod and a chelicerate. The large subunit rRNA of P. opilio is found to have more helices conserved than in other arthropods, while the small subunit rRNA shows a complexity similar to that of other arthropods. These comparisons suggest that a reduction in rRNA complexity occurred in Pancrustacea after the divergence of Pancrustacea and Chelicerata from a common ancestor. The gene arrangement of the mitochondrial genome of P. opilio is compared with the gene order of taxa from all seven other orders of arachnids for which representative mitochondrial genomes have been sequenced. Taxa from five of these seven orders possess gene arrangements identical to that of Limulus polyphemus, and P. opilio is found to have a similar arrangement. However, in P. opilio, some genes near the putative control region are rearranged, with the suite of genes encoding tRNA(Gln), the control region, and tRNA(Ile) located downstream of the two ribosomal RNA genes, and upstream of where they are typically located in chelicerates. The genome encodes only 21 of the typical 22 mitochondrial tRNA genes and lacks the gene for tRNA(Leu(CUN)). The protein-coding genes in the mitochondrial genome of P. opilio show a significantly decreased use of codons recognized by tRNA(Leu(CUN)), likely due to selection to utilize the more specific tRNA(Leu(UUR)) anticodon. The gene arrangement and lack of a tRNA(Leu(CUN)) gene in P. opilio is most parsimoniously explained by the occurrence of at least two translocation events, one of which probably destroyed the function of the tRNA(Leu(CUN)) gene. Phylogenetic relationships among the major orders of arachnids are inferred, using all 13 mt protein-coding genes, and gene rearrangements are mapped onto the phylogeny. The phylogenetic analyses are unable to resolve the placement of P. opilio but are generally consistent with an early divergence of members of the Dromopoda (harvestmen, scorpions, and solifuges) from the Micruran arachnids (spiders, whip spiders, vinegaroons, ricinuleids, and mites). However, unlike some morphologically based phylogenetic analyses, the existence of a clade of Dromopoda is not supported. While data on genome arrangement and gene loss do not provide further information to help resolve relationships among the arachnid orders, they distinguish some groups of arachnids, distinguish chelicerates from other arthropods, and further clarify the ancestral gene order of this diverse group of arthropods.

The complete mitochondrial genome of the house dust mite Dermatophagoides pteronyssinus (Trouessart): a novel gene arrangement among arthropods

BACKGROUND

The apparent scarcity of available sequence data has greatly impeded evolutionary studies in Acari (mites and ticks). This subclass encompasses over 48,000 species and forms the largest group within the Arachnida. Although mitochondrial genomes are widely utilised for phylogenetic and population genetic studies, only 20 mitochondrial genomes of Acari have been determined, of which only one belongs to the diverse order of the Sarcoptiformes. In this study, we describe the mitochondrial genome of the European house dust mite Dermatophagoides pteronyssinus, the most important member of this largely neglected group.

RESULTS

The mitochondrial genome of D. pteronyssinus is a circular DNA molecule of 14,203 bp. It contains the complete set of 37 genes (13 protein coding genes, 2 rRNA genes and 22 tRNA genes), usually present in metazoan mitochondrial genomes. The mitochondrial gene order differs considerably from that of other Acari mitochondrial genomes. Compared to the mitochondrial genome of Limulus polyphemus, considered as the ancestral arthropod pattern, only 11 of the 38 gene boundaries are conserved. The majority strand has a 72.6% AT-content but a GC-skew of 0.194. This skew is the reverse of that normally observed for typical animal mitochondrial genomes. A microsatellite was detected in a large non-coding region (286 bp), which probably functions as the control region. Almost all tRNA genes lack a T-arm, provoking the formation of canonical cloverleaf tRNA-structures, and both rRNA genes are considerably reduced in size. Finally, the genomic sequence was used to perform a phylogenetic study. Both maximum likelihood and Bayesian inference analysis clustered D. pteronyssinus with Steganacarus magnus, forming a sistergroup of the Trombidiformes.

CONCLUSION

Although the mitochondrial genome of D. pteronyssinus shares different features with previously characterised Acari mitochondrial genomes, it is unique in many ways. Gene order is extremely rearranged and represents a new pattern within the Acari. Both tRNAs and rRNAs are truncated, corroborating the theory of the functional co-evolution of these molecules. Furthermore, the strong and reversed GC- and AT-skews suggest the inversion of the control region as an evolutionary event. Finally, phylogenetic analysis using concatenated mt gene sequences succeeded in recovering Acari relationships concordant with traditional views of phylogeny of Acari.

First divergence time estimate of spiders, scorpions, mites and ticks (subphylum: Chelicerata) inferred from mitochondrial phylogeny

Spiders, scorpions, mites and ticks (chelicerates) form one of the most diverse groups of arthropods on land, but their origin and times of diversification are not yet established. We estimated, for the first time, the molecular divergence times for these chelicerates using complete mitochondrial sequences from 25 taxa. All mitochondrial genes were evaluated individually or after concatenation. Sequences belonging to three missing genes (ND3, 6, and tRNA-Asp) from three taxa, as well as the faster-evolving ribosomal RNAs (12S and 16S), tRNAs, and the third base of each codon from 11 protein-coding genes (PCGs) (COI-III, CYTB, ATP8, 6, ND1-2, 4L, and 4-5), were identified and removed. The remaining concatenated sequences from 11 PCGs produced a completely resolved phylogenetic tree and confirmed that all chelicerates are monophyletic. Removing the third base from each codon was essential to resolve the phylogeny, which allowed deep divergence times to be calculated using three nodes calibrated with upper and lower priors. Our estimates indicate that the orders and classes of spiders, scorpions, mites, and ticks diversified in the late Paleozoic, much earlier than previously reported from fossil date estimates. The divergence time estimated for ticks suggests that their first land hosts could have been amphibians rather than reptiles. Using molecular data, we separated the spider-scorpion clades and estimated their divergence times at 397 +/- 23 million years ago. Algae, fungi, plants, and animals, including insects, were well established on land when these chelicerates diversified. Future analyses, involving mitochondrial sequences from additional chelicerate taxa and the inclusion of nuclear genes (or entire genomes) will provide a more complete picture of the evolution of the Chelicerata, the second most abundant group of animals on earth.

The complete mitochondrial genome of the sexual oribatid mite Steganacarus magnus: genome rearrangements and loss of tRNAs

Background

Complete mitochondrial (mt) genomes and the gene rearrangements therein are increasingly used as molecular markers for investigating phylogenetic relationships, especially for elucidating deep splits. Contributing to the complete mt genomes of arthropods, especially Arachnida, available so far, we provide the first complete mt genome of a sarcoptiform mite species, the sexually reproducing oribatid mite Steganacarus magnus (Acari, Oribatida) which was determined by sequencing of long PCR products.

Results

The mt genome of S. magnus lacks 16 tRNAs, only those for leucine, histidine, proline, tryptophan, glutamine and serine are present. Within those tRNAs only tRNA-His and tRNA-Pro have kept their original position, the others are translocated. Furthermore, the mt genome of S. magnus consists of 13,818 bp and it is composed of 13 protein-coding genes and two genes for the ribosomal RNA subunits that are typically found in metazoan mt genomes. The gene order in S. magnus differs from the hypothetical ancestral chelicerate arrangement as conserved in Limulus polyphemus: instead of nad1-rrnL-rrnS-LNR-nad2 (tRNAs excluded) S. magnus is nad2-rrnL-nad1-rrnS-LNR. Phylogenetic analyses of a concatenated amino acid dataset of all mt protein-coding genes of 28 arthropod species suggest a sister-group relationship of sarcoptiform and prostigmatid mites (S. magnus and Leptotrombidium).

Conclusion

The mt gene arrangement of S. magnus differs from the hypothetical ground plan of arthropods and from that of other mites further contributing to the variety of mt gene arrangements found in Arachnida. The unexpected lack of tRNAs is enigmatic, probably showing that the loss of mt genes is an ongoing evolutionary process. For solving phylogenetic relationships of oribatid mite lineages and their position within Acari further complete mt genomes are needed.

Evolution of the mitochondrial genome of Metazoa as exemplified by comparison of congeneric species

The mitochondrial genome (mtDNA) of Metazoa is a good model system for evolutionary genomic studies and the availability of more than 1000 sequences provides an almost unique opportunity to decode the mechanisms of genome evolution over a large phylogenetic range. In this paper, we review several structural features of the metazoan mtDNA, such as gene content, genome size, genome architecture and the new parameter of gene strand asymmetry in a phylogenetic framework. The data reviewed here show that: (1) the plasticity of Metazoa mtDNA is higher than previously thought and mainly due to variation in number and location of tRNA genes; (2) an exceptional trend towards stabilization of genomic features occurred in deuterostomes and was exacerbated in vertebrates, where gene content, genome architecture and gene strand asymmetry are almost invariant. Only tunicates exhibit a very high degree of genome variability comparable to that found outside deuterostomes. In order to analyse the genomic evolutionary process at short evolutionary distances, we have also compared mtDNAs of species belonging to the same genus: the variability observed in congeneric species significantly recapitulates the evolutionary dynamics observed at higher taxonomic ranks, especially for taxa showing high levels of genome plasticity and/or fast nucleotide substitution rates. Thus, the analysis of congeneric species promises to be a valuable approach for the assessment of the mtDNA evolutionary trend in poorly or not yet sampled metazoan groups.

The complete mitochondrial genome of Pseudocellus pearsei (Chelicerata: Ricinulei) and a comparison of mitochondrial gene rearrangements in Arachnida

Background

Mitochondrial genomes are widely utilized for phylogenetic and population genetic analyses among animals. In addition to sequence data the mitochondrial gene order and RNA secondary structure data are used in phylogenetic analyses. Arachnid phylogeny is still highly debated and there is a lack of sufficient sequence data for many taxa. Ricinulei (hooded tickspiders) are a morphologically distinct clade of arachnids with uncertain phylogenetic affinities.

Results

The first complete mitochondrial DNA genome of a member of the Ricinulei, Pseudocellus pearsei (Arachnida: Ricinulei) was sequenced using a PCR-based approach. The mitochondrial genome is a typical circular duplex DNA molecule with a size of 15,099 bp, showing the complete set of genes usually present in bilaterian mitochondrial genomes. Five tRNA genes (trnW, trnY, trnN, trnL(CUN), trnV) show different relative positions compared to other Chelicerata (e.g. Limulus polyphemus, Ixodes spp.). We propose that two events led to this derived gene order: (1) a tandem duplication followed by random deletion and (2) an independent translocation of trnN. Most of the inferred tRNA secondary structures show the common cloverleaf pattern except tRNA-Glu where the TψC-arm is missing. In phylogenetic analyses (maximum likelihood, maximum parsimony, Bayesian inference) using concatenated amino acid and nucleotide sequences of protein-coding genes the basal relationships of arachnid orders remain unresolved.

Conclusion

Phylogenetic analyses (ML, MP, BI) of arachnid mitochondrial genomes fail to resolve interordinal relationships of Arachnida and remain in a preliminary stage because there is still a lack of mitogenomic data from important taxa such as Opiliones and Pseudoscorpiones. Gene order varies considerably within Arachnida – only eight out of 23 species have retained the putative arthropod ground pattern. Some gene order changes are valuable characters in phylogenetic analysis of intraordinal relationships, e.g. in Acari.

The mitochondrial genome of Phallusia mammillata and Phallusia fumigata (Tunicata, Ascidiacea): high genome plasticity at intra-genus level

Background

Within Chordata, the subphyla Vertebrata and Cephalochordata (lancelets) are characterized by a remarkable stability of the mitochondrial (mt) genome, with constancy of gene content and almost invariant gene order, whereas the limited mitochondrial data on the subphylum Tunicata suggest frequent and extensive gene rearrangements, observed also within ascidians of the same genus.

Results

To confirm this evolutionary trend and to better understand the evolutionary dynamics of the mitochondrial genome in Tunicata Ascidiacea, we have sequenced and characterized the complete mt genome of two congeneric ascidian species, Phallusia mammillata and Phallusia fumigata (Phlebobranchiata, Ascidiidae). The two mtDNAs are surprisingly rearranged, both with respect to one another and relative to those of other tunicates and chordates, with gene rearrangements affecting both protein-coding and tRNA genes. The new data highlight the extraordinary variability of ascidian mt genome in base composition, tRNA secondary structure, tRNA gene content, and non-coding regions (number, size, sequence and location). Indeed, both Phallusia genomes lack the trnD gene, show loss/acquisition of DHU-arm in two tRNAs, and have a G+C content two-fold higher than other ascidians. Moreover, the mt genome of P. fumigata presents two identical copies of trnI, an extra tRNA gene with uncertain amino acid specificity, and four almost identical sequence regions. In addition, a truncated cytochrome b, lacking a C-terminal tail that commonly protrudes into the mt matrix, has been identified as a new mt feature probably shared by all tunicates.

Conclusion

The frequent occurrence of major gene order rearrangements in ascidians both at high taxonomic level and within the same genus makes this taxon an excellent model to study the mechanisms of gene rearrangement, and renders the mt genome an invaluable phylogenetic marker to investigate molecular biodiversity and speciation events in this largely unexplored group of basal chordates.