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

Genomic and metagenomic analyses of the domestic mite Tyrophagus putrescentiae identify it as a widespread environmental contaminant and a host of a basal, mite-specific Wolbachia lineage (supergroup Q)

Tyrophagus putrescentiae (mould mite) is a global, microscopic trophic generalist that commonly occurs in various human-created habitats, causing allergies and damaging stored food. Its ubiquity and extraordinary ability to penetrate research samples or cultures through air currents or by active walking through tights spaces (such as treads of screw caps) may lead to sample contamination and introduction of its DNA to research materials in the laboratory. This prompts a thorough investigation into potential sequence contamination in public genomic databases. The trophic success of T. putrescentiae is primarily attributed to the symbiotic bacteria housed in specialized internal mite structures, facilitating adaptation to varied nutritional niches. However, recent work suggests that horizontal transfer of bacterial/fungal genes related to nutritional functionality may also contribute to the mite's trophic versatility. This aspect requires independent confirmation. Additionally, T. putrescentiae harbors an uncharacterized and genetically divergent bacterium, Wolbachia, displaying blocking and microbiome-modifying effects. The phylogenomic position and supergroup assignment of this bacterium are unknown. Here, we sequenced and assembled the T. putrescentiae genome, analyzed its microbiome, and performed detailed phylogenomic analyses of the mite-specific Wolbachia. We show that T. putrescentiae DNA is a substantial source of contamination of research samples. Its DNA may inadvertently be co-extracted with the DNA of the target organism, eventually leading to sequence contamination in public databases. We identified a diversity of bacterial species associated with T. putrescentiae, including those capable of rapidly developing antibiotic resistance, such as Escherichia coli. Despite the presence of diverse bacterial communities in T. putrescentiae, we did not detect any recent horizontal gene transfers in this mite species and/or in astigmatid (domestic) mites in general. Our phylogenomic analysis of Wolbachia recovered a basal, mite-specific lineage (supergroup Q) represented by two Wolbachia spp. from the mould mite and a gall-inducing plant mite. Fluorescence in situ hybridization confirmed the presence of Wolbachia inside the mould mite. The discovery of an early derivative Wolbachia lineage (supergroup Q) in two phylogenetically unrelated and ecologically dissimilar mites suggests that this endosymbiotic bacterial lineage formed a long-term association with mites. This finding provides a unique insight into the early evolution and host associations of Wolbachia. Further discoveries of Wolbachia diversity in acariform mites are anticipated.

Mitochondrial metagenomics reveal the independent colonization of the world's coasts by intertidal oribatid mites (Acari, Oribatida, Ameronothroidea)

Oribatid mites are an ancient group that already roamed terrestrial ecosystems in the early and middle Devonian. The superfamily of Ameronothroidea, a supposedly monophyletic lineage, represents the only group of oribatid mites that has successfully invaded the marine coastal environment. By using mitogenome data and nucleic ribosomal RNA genes (18S, 5.8S, 28S), we show that Ameronothroidea are a paraphyletic assemblage and that the land-to-sea transition happened three times independently. Common ancestors of the tropical Fortuyniidae and Selenoribatidae were the first to colonize the coasts and molecular calibration of our phylogeny dates this event to a period in the Triassic and Jurassic era (225-146 mya), whereas present-day distribution indicates that this event might have happened early in this period during the Triassic, when the supercontinent Pangaea still existed. The cold temperate northern hemispheric Ameronothridae colonized the marine littoral later in the late Jurassic-Early Cretaceous and had an ancient distribution on Laurasian coasts. The third and final land-to-sea transition happened in the same geological period, but approx. 30 my later when ancestors of Podacaridae invaded coastal marine environments of the Gondwanan landmasses.

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.

Diversity and Distribution of Mites (ACARI) Revealed by Contamination Survey in Public Genomic Databases

Acari (mites and ticks) are a biodiverse group of microarthropods within the Arachnida. Because of their diminutive size, mites are often overlooked. We hypothesized that mites, like other closely related microorganisms, could also contaminate public genomic database. Here, using a strategy based on DNA barcodes previously reported, we scanned contaminations related to mites (Acari, exclusive of Ixodida) in Genbank WGS/TSA database. In 22,114 assemblies (17,845 animal and 4269 plant projects), 1717 contigs in 681 assemblies (3.1%) were detected as mite contaminations. Additional taxonomic analysis showed the following: (1) most of the contaminants (1445/1717) were from the specimens of Magnoliopsida, Insecta and Pinopsida; (2) the contamination rates were higher in plant or TSA projects; (3) mite distribution among different classes of hosts varied considerably. Additional phylogenetic analysis of these contaminated contigs further revealed complicated mite-host associations. Overall, we conducted a first systemic survey and analysis of mite contaminations in public genomic database, and these DNA barcode related mite contigs will provide a valuable resource of information for understanding the diversity and phylogeny of mites.

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.

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.

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.

Adaptation of the Romanomermis culicivorax CCA-Adding Enzyme to Miniaturized Armless tRNA Substrates

The mitochondrial genome of the nematode Romanomermis culicivorax encodes for miniaturized hairpin-like tRNA molecules that lack D- as well as T-arms, strongly deviating from the consensus cloverleaf. The single tRNA nucleotidyltransferase of this organism is fully active on armless tRNAs, while the human counterpart is not able to add a complete CCA-end. Transplanting single regions of the Romanomermis enzyme into the human counterpart, we identified a beta-turn element of the catalytic core that-when inserted into the human enzyme-confers full CCA-adding activity on armless tRNAs. This region, originally identified to position the 3'-end of the tRNA primer in the catalytic core, dramatically increases the enzyme's substrate affinity. While conventional tRNA substrates bind to the enzyme by interactions with the T-arm, this is not possible in the case of armless tRNAs, and the strong contribution of the beta-turn compensates for an otherwise too weak interaction required for the addition of a complete CCA-terminus. This compensation demonstrates the remarkable evolutionary plasticity of the catalytic core elements of this enzyme to adapt to unconventional tRNA substrates.

Naturally Occurring tRNAs With Non-canonical Structures

Transfer RNA (tRNA) is the central molecule in genetically encoded protein synthesis. Most tRNA species were found to be very similar in structure: the well-known cloverleaf secondary structure and L-shaped tertiary structure. Furthermore, the length of the acceptor arm, T-arm, and anticodon arm were found to be closely conserved. Later research discovered naturally occurring, active tRNAs that did not fit the established 'canonical' tRNA structure. This review discusses the non-canonical structures of some well-characterized natural tRNA species and describes how these structures relate to their role in translation. Additionally, we highlight some newly discovered tRNAs in which the structure-function relationship is not yet fully understood.

Comparative analysis of mitochondrial genomes among the subfamily Sarcophaginae (Diptera: Sarcophagidae) and phylogenetic implications

The subfamily Sarcophaginae is extremely diverse in morphology, habit and geographical distribution, and usually considered to be of significant ecological, medical, and forensic significance. In the present study, 18 mitochondrial genomes (mitogenomes) of sarcophagid flies were first obtained. The rearrangement and orientation of genes were identical with that of ancestral insects. The degrees of compositional heterogeneity in the datasets were extremely low. Furthermore, 13 protein-coding genes were evolving under purifying selection. The phylogenic relationship of the genus-group taxa Boettcheria + (Sarcophaga + (Peckia + (Ravinia + Oxysarcodexia))) was strongly supported. Four subgenera were recovered as monophyletic (Liopygia, Liosarcophaga, Pierretia, Heteronychia) in addition to Parasarcophaga as polyphyletic. The sister-relationships between S. dux and S. aegyptiaca, S. pingi and S. kawayuensis were recovered, respectively. Moreover, the molecular phylogenetic relationships among the subgenera Helicophagella, Kozlovea, Kramerea, Pandelleisca, Phallocheira, Pseudothyrsocnema, Sinonipponia and Seniorwhitea were rarely put forward prior to this study. This study provides insight into the population genetics, molecular biology, and phylogeny for the subfamily Sarcophaginae, especially for the subgeneric classification of Sarcophaga. However, compared with the enormous species diversity of flesh flies, the available mitogenomes are still limited for recovering the phylogeny of Sarcophaginae.

De novo sequence of the mitochondrial genome of Tyrophagus putrescentiae (Acari: Sarcoptiformes) including 22 tRNA sequences and the largest non-coding region

In this study, we de novo sequenced and analyzed the circular mitochondrial genome (mitogenome) of Tyrophagus putrescentiae. It was 14,156 bp long and contained a complete set of 37 genes, contrary to the initial published sequences; it included 22 tRNA sequences and the largest non-coding region. The mtDNA gene order of T. putrescentiae was found to be identical to that of Aleuroglyphus ovatus, Caloglyphus berlesei, and Rhizoglyphus robini (all Acaroidea). Most tRNAs of T. putrescentiae lack at least a D-arm or T-arm. Tyrophagus putrescentiae tRNAs also shared considerable structural and sequence similarity with the tRNAs of other reported Acaroidea species that have the full set of tRNAs. The largest non-coding region was located between trnF and trnS1, and it contained a microsatellite-like (AT)_n sequence, short palindromic sequences, and several hairpin loops, as observed in other reported Acaroidea species (excepting Tyrophagus longior).

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.

Small but large enough: structural properties of armless mitochondrial tRNAs from the nematode Romanomermis culicivorax

As adapter molecules to convert the nucleic acid information into the amino acid sequence, tRNAs play a central role in protein synthesis. To fulfill this function in a reliable way, tRNAs exhibit highly conserved structural features common in all organisms and in all cellular compartments active in translation. However, in mitochondria of metazoans, certain dramatic deviations from the consensus tRNA structure are described, where some tRNAs lack the D- or T-arm without losing their function. In Enoplea, this miniaturization comes to an extreme, and functional mitochondrial tRNAs can lack both arms, leading to a considerable size reduction. Here, we investigate the secondary and tertiary structure of two such armless tRNAs from Romanomermis culicivorax. Despite their high AU content, the transcripts fold into a single and surprisingly stable hairpin structure, deviating from standard tRNAs. The three-dimensional form is boomerang-like and diverges from the standard L-shape. These results indicate that such unconventional miniaturized tRNAs can still fold into a tRNA-like shape, although their length and secondary structure are very unusual. They highlight the remarkable flexibility of the protein synthesis apparatus and suggest that the translational machinery of Enoplea mitochondria may show compensatory adaptations to accommodate these armless tRNAs for efficient translation.

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.

Cicada Endosymbionts Have tRNAs That Are Correctly Processed Despite Having Genomes That Do Not Encode All of the tRNA Processing Machinery

The smallest bacterial genomes, in the range of about 0.1 to 0.5 million base pairs, are commonly found in the nutritional endosymbionts of insects. These tiny genomes are missing genes that encode proteins and RNAs required for the translation of mRNAs, one of the most highly conserved and important cellular processes. In this study, we found that the bacterial endosymbionts of cicadas have genomes which encode incomplete tRNA sets and lack genes required for tRNA processing. Nevertheless, we found that endosymbiont tRNAs are correctly processed at their 5′ and 3′ ends and, surprisingly, that mostly exist as tRNA halves. We hypothesize that the cicada host must supply its symbionts with these missing tRNA processing activities.

Gene loss and genome reduction are defining characteristics of endosymbiotic bacteria. The most highly reduced endosymbiont genomes have lost numerous essential genes related to core cellular processes such as replication, transcription, and translation. Computational gene predictions performed for the genomes of the two bacterial symbionts of the cicada Diceroprocta semicincta, “Candidatus Hodgkinia cicadicola” (Alphaproteobacteria) and “Ca. Sulcia muelleri” (Bacteroidetes), have found only 26 and 16 tRNA genes and 15 and 10 aminoacyl tRNA synthetase genes, respectively. Furthermore, the original “Ca. Hodgkinia cicadicola” genome annotation was missing several essential genes involved in tRNA processing, such as those encoding RNase P and CCA tRNA nucleotidyltransferase as well as several RNA editing enzymes required for tRNA maturation. How these cicada endosymbionts perform basic translation-related processes remains unknown. Here, by sequencing eukaryotic mRNAs and total small RNAs, we show that the limited tRNA set predicted by computational annotation of “Ca. Sulcia muelleri” and “Ca. Hodgkinia cicadicola” is likely correct. Furthermore, we show that despite the absence of genes encoding tRNA processing activities in the symbiont genomes, symbiont tRNAs have correctly processed 5′ and 3′ ends and seem to undergo nucleotide modification. Surprisingly, we found that most “Ca. Hodgkinia cicadicola” and “Ca. Sulcia muelleri” tRNAs exist as tRNA halves. We hypothesize that “Ca. Sulcia muelleri” and “Ca. Hodgkinia cicadicola” tRNAs function in bacterial translation but require host-encoded enzymes to do so.

Population and Culture Age Influence the Microbiome Profiles of House Dust Mites

Interactions with microorganisms might enable house dust mites (HDMs) to derive nutrients from difficult-to-digest structural proteins and to flourish in human houses. We tested this hypothesis by investigating the effects of changes in the mite culture growth and population of two HDM species on HDM microbiome composition and fitness. Growing cultures of laboratory and industrial allergen-producing populations of Dermatophagoides farinae (DFL and DFT, respectively) and Dermatophagoides pteronyssinus (DPL and DPT, respectively) were sampled at four time points. The symbiotic microorganisms of the mites were characterized by DNA barcode sequencing and quantified by qPCR using universal/specific primers. The population growth of mites and nutrient contents of mite bodies were measured and correlated with the changes in bacteria in the HDM microbiome. The results showed that both the population and culture age significantly influenced the microbiome profiles. Cardinium formed 93% and 32% of the total sequences of the DFL and DFT bacterial microbiomes, respectively, but this bacterial species was less abundant in the DPL and DPT microbiomes. Staphylococcus abundance was positively correlated with increased glycogen contents in the bodies of mites, and increased abundances of Aspergillus, Candida, and Kocuria were correlated with increased lipid contents in the bodies of mites. The xerophilic fungus Wallemia accounted for 39% of the fungal sequences in the DPL microbiome, but its abundance was low in the DPT, DFL, and DFT microbiomes. With respect to the mite culture age, we made three important observations: the mite population growth from young cultures was 5-8-fold higher than that from old cultures; specimens from old cultures had greater abundances of fungi and bacteria in their bodies; and yeasts predominated in the gut contents of specimens from young cultures, whereas filamentous mycelium prevailed in specimens from old cultures. Our results are consistent with the hypothesis that mites derive nutrients through associations with microorganisms.

Comparative analysis of mitochondrial genomes among four species of muscid flies (Diptera: Muscidae) and its phylogenetic implications

Muscidae, commonly known as house flies and their close relatives, is one of the dipteran insects of recognized medical, veterinary, and ecological importance. Mitochondrial genomes (Mitogenomes) have been widely used for exploring phylogenetic analysis and taxonomic diagnosis due to the difficulty in distinguishing them morphologically. In this study, our complete mitogenomes of muscid flies were sequenced and aligned, which ranged from 15,117 bp (Synthesiomyia nudiseta) to 16,089 bp (Musca sorbens) in length, and contained a typical circular molecule comprising 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs) and a non-coding control region. The order and orientation of genes were identical with that from the ancestral insects. The phylogenetic analysis based on the mitochondrial genes indicated that the subfamily relationships within Muscidae were reconstructed as (Mydaeinae (Muscinae (Reinwardtiinae + Azeliinae))). Similar tree topologies were recovered from both Maximum Likelihood (ML) and Bayesian Inference (BI) analysis. Furthermore, we compared the phylogenetic analyses that were constructed using internal transcribed spacer 2 (ITS2), elongation factor-1α (EF-1α), 13 PCGs and 13 PCGs + ITS2 + EF-1α, respectively. Combined analysis of nuclear gene partitions improved support and resolution for resulting topologies but the positions of branches were obviously inconsistent due to limited species. More mitogenomes should be sequenced representing various taxonomic levels, especially close related species, which will enhance our understanding of phylogenetic relationships among muscids.

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 complete mitochondrial genome of the feather mite Trouessartia rubecula Jablonska, 1968 (Astigmata: Analgoidea: Trouessartiidae)

We assembled and annotated the complete mitochondrial genome of Trouessartia rubecula, the first feather mite complete mitochondrial genome from the largest feather mite superfamily Analgoidea (ca. 1150 spp). The mitogenome was composed of 13 protein, 17 tRNA, and 2 rRNA-coding genes and was 14,125 bp in length.

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.

Populations of Stored Product Mite Tyrophagus putrescentiae Differ in Their Bacterial Communities

Background:Tyrophagus putrescentiae colonizes different human-related habitats and feeds on various post-harvest foods. The microbiota acquired by these mites can influence the nutritional plasticity in different populations. We compared the bacterial communities of five populations of T. putrescentiae and one mixed population of T. putrescentiae and T. fanetzhangorum collected from different habitats.

Material: The bacterial communities of the six mite populations from different habitats and diets were compared by Sanger sequencing of cloned 16S rRNA obtained from amplification with universal eubacterial primers and using bacterial taxon-specific primers on the samples of adults/juveniles or eggs. Microscopic techniques were used to localize bacteria in food boli and mite bodies. The morphological determination of the mite populations was confirmed by analyses of CO1 and ITS fragment genes.

Results: The following symbiotic bacteria were found in compared mite populations: Wolbachia (two populations), Cardinium (five populations), Bartonella-like (five populations), Blattabacterium-like symbiont (three populations), and Solitalea-like (six populations). From 35 identified OTUs₉₇, only Solitalea was identified in all populations. The next most frequent and abundant sequences were Bacillus, Moraxella, Staphylococcus, Kocuria, and Microbacterium. We suggest that some bacterial species may occasionally be ingested with food. The bacteriocytes were observed in some individuals in all mite populations. Bacteria were not visualized in food boli by staining, but bacteria were found by histological means in ovaria of Wolbachia-infested populations.

Conclusion: The presence of Blattabacterium-like, Cardinium, Wolbachia, and Solitalea-like in the eggs of T. putrescentiae indicates mother to offspring (vertical) transmission. Results of this study indicate that diet and habitats influence not only the ingested bacteria but also the symbiotic bacteria of T. putrescentiae.

Mitochondrial Genome Sequence of the Scabies Mite Provides Insight into the Genetic Diversity of Individual Scabies Infections

The scabies mite, Sarcoptes scabiei, is an obligate parasite of the skin that infects humans and other animal species, causing scabies, a contagious disease characterized by extreme itching. Scabies infections are a major health problem, particularly in remote Indigenous communities in Australia, where co-infection of epidermal scabies lesions by Group A Streptococci or Staphylococcus aureus is thought to be responsible for the high rate of rheumatic heart disease and chronic kidney disease. We collected and separately sequenced mite DNA from several pools of thousands of whole mites from a porcine model of scabies (S. scabiei var. suis) and two human patients (S. scabiei var. hominis) living in different regions of northern Australia. Our sequencing samples the mite and its metagenome, including the mite gut flora and the wound micro-environment. Here, we describe the mitochondrial genome of the scabies mite. We developed a new de novo assembly pipeline based on a bait-and-reassemble strategy, which produced a 14 kilobase mitochondrial genome sequence assembly. We also annotated 35 genes and have compared these to other Acari mites. We identified single nucleotide polymorphisms (SNPs) and used these to infer the presence of six haplogroups in our samples, Remarkably, these fall into two closely-related clades with one clade including both human and pig varieties. This supports earlier findings that only limited genetic differences may separate some human and animal varieties, and raises the possibility of cross-host infections. Finally, we used these mitochondrial haplotypes to show that the genetic diversity of individual infections is typically small with 1-3 distinct haplotypes per infestation.

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.

First complete mitochondrial genome of the South American annual fish Austrolebias charrua (Cyprinodontiformes: Rivulidae): peculiar features among cyprinodontiforms mitogenomes

Background

Among teleosts, the South American genus Austrolebias (Cyprinodontiformes: Rivulidae) includes 42 taxa of annual fishes divided into five different species groups. It is a monophyletic genus, but morphological and molecular data do not resolve the relationship among intrageneric clades and high rates of substitution have been previously described in some mitochondrial genes. In this work, the complete mitogenome of a species of the genus was determined for the first time. We determined its structure, gene order and evolutionary peculiar features, which will allow us to evaluate the performance of mitochondrial genes in the phylogenetic resolution at different taxonomic levels.

Results

Regarding gene content and order, the circular mitogenome of A. charrua (17,271 pb) presents the typical pattern of vertebrate mitogenomes. It contains the full complement of 13 proteins-coding genes, 22 tRNA, 2 rRNA and one non-coding control region. Notably, the tRNA-Cys was only 57 bp in length and lacks the D-loop arm. In three full sibling individuals, heteroplasmatic condition was detected due to a total of 12 variable sites in seven protein-coding genes. Among cyprinodontiforms, the mitogenome of A. charrua exhibits the lowest G+C content (37 %) and GCskew, as well as the highest strand asymmetry with a net difference of T over A at 1st and 3rd codon positions. Considering the 12 coding-genes of the H strand, correspondence analyses of nucleotide composition and codon usage show that A and T at 1st and 3rd codon positions have the highest weight in the first axis, and segregate annual species from the other cyprinodontiforms analyzed. Given the annual life-style, their mitogenomes could be under different selective pressures. All 13 protein-coding genes are under strong purifying selection and we did not find any significant evidence of nucleotide sites showing episodic selection (dN >dS) at annual lineages. When fast evolving third codon positions were removed from alignments, the “supergene” tree recovers our reference species phylogeny as well as the Cytb, ND4L and ND6 genes. Therefore, third codon positions seem to be saturated in the aforementioned coding regions at intergeneric Cyprinodontiformes comparisons.

Conclusions

The complete mitogenome obtained in present work, offers relevant data for further comparative studies on molecular phylogeny and systematics of this taxonomic controversial endemic genus of annual fishes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2090-3) contains supplementary material, which is available to authorized users.

Precise mapping and dynamics of tRNA-derived fragments (tRFs) in the development of Triops cancriformis (tadpole shrimp)

Background

In a deep sequencing analysis of small RNAs prepared from a living fossil, the tadpole shrimp Triops cancriformis, a 32-nt small RNA was specifically detected in the adult stage. A nucleotide sequence comparison between the 32-nt small RNA and predicted tRNA sequences in the draft nuclear genomic DNA showed that the small RNA was derived from tRNA^Gly(GCC). To determine the overall features of the tRNA-derived fragments (tRFs) of T. cancriformis, the small RNA sequences in each of the six developmental stages (egg, 1st − 4th instar larvae, and adult) were compared with the mitochondrial and nuclear tRNA sequences.

Results

We found that the tRFs were derived from mitochondrial and nuclear tRNAs corresponding to 16 and 39 anticodons, respectively. The total read number of nuclear tRFs was approximately 400 times larger than the number of mitochondrial tRFs. Interestingly, the main regions in each parental tRNA from which these tRFs were derived differed, depending on the parental anticodon. Mitochondrial tRF^Ser(GCU)s were abundantly produced from the 5’ half regions of the parental tRNA, whereas mitochondrial tRF^Val(UAC)s were mainly produced from the 3’ end regions. Highly abundant nuclear tRFs, tRF^Gly(GCC)s, tRF^Gly(CCC)s, tRF^Glu(CUC)s, and tRF^Lys(CUU)s were derived from the 5’ half regions of the parental tRNAs. Further analysis of the tRF read counts in the individual developmental stages suggested that the expression of mitochondrial and nuclear tRFs differed during the six stages. Based on these data, we precisely summarized the positions of the tRFs in their parental tRNAs and their expression changes during development.

Conclusions

Our results reveal the entire dynamics of the tRFs from both the nuclear and mitochondrial genomes of T. cancriformis and indicate that the majority of tRFs in the cell are derived from nuclear tRNAs. This study provides the first examples of developmentally expressed mitochondrial tRFs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0245-5) contains supplementary material, which is available to authorized users.

Comparative mitogenomic analysis of the superfamily Pentatomoidea (Insecta: Hemiptera: Heteroptera) and phylogenetic implications

BACKGROUND

Insect mitochondrial genomes (mitogenomes) are the most extensively used genetic marker for evolutionary and population genetics studies of insects. The Pentatomoidea superfamily is economically important and the largest superfamily within Pentatomomorpha with over 7,000 species. To better understand the diversity and evolution of pentatomoid species, we sequenced and annotated the mitogenomes of Eurydema gebleri and Rubiconia intermedia, and present the first comparative analysis of the 11 pentatomoid mitogenomes that have been sequenced to date.

RESULTS

We obtained the complete mitogenome of Eurydema gebleri (16,005 bp) and a nearly complete mitogenome of Rubiconia intermedia (14,967 bp). Our results show that gene content, gene arrangement, base composition, codon usage, and mitochondrial transcription termination factor sequences are highly conserved in pentatomoid species, especially for species in the same family. Evolutionary rate analyses of protein-coding genes reveal that the highest and lowest rates are found in atp8 and cox1 and distinctive evolutionary patterns are significantly correlated with the G + C content of genes. We inferred the secondary structures for two rRNA genes for eleven pentatomoid species, and identify some conserved motifs of RNA structures in Pentatomidea. All tRNA genes in pentatomoid mitogenomes have a canonical cloverleaf secondary structure, except for two tRNAs (trnS1 and trnV) which appear to lack the dihydrouridine arm. Regions that are A + T-rich have several distinct characteristics (e.g. size variation and abundant tandem repeats), and have potential as species or population level molecular markers. Phylogenetic analyses based on mitogenomic data strongly support the monophyly of Pentatomoidea, and the estimated phylogenetic relationships are: (Urostylididae + (Plataspidae + (Pentatomidae + (Cydnidae + (Dinidoridae + Tessaratomidae))))).

CONCLUSIONS

This comparative mitogenomic analysis sheds light on the architecture and evolution of mitogenomes in the superfamily Pentatomoidea. Mitogenomes can be effectively used to resolve phylogenetic relationships of pentatomomorphan insects at various taxonomic levels. Sequencing more mitogenomes at various taxonomic levels, particularly from closely related species, will improve the annotation accuracy of mitochondrial genes, as well as greatly enhance our understanding of mitogenomic evolution and phylogenetic relationships in pentatomoids.

The identity of the discriminator base has an impact on CCA addition

CCA-adding enzymes synthesize and maintain the C-C-A sequence at the tRNA 3'-end, generating the attachment site for amino acids. While tRNAs are the most prominent substrates for this polymerase, CCA additions on non-tRNA transcripts are described as well. To identify general features for substrate requirement, a pool of randomized transcripts was incubated with the human CCA-adding enzyme. Most of the RNAs accepted for CCA addition carry an acceptor stem-like terminal structure, consistent with tRNA as the main substrate group for this enzyme. While these RNAs show no sequence conservation, the position upstream of the CCA end was in most cases represented by an adenosine residue. In tRNA, this position is described as discriminator base, an important identity element for correct aminoacylation. Mutational analysis of the impact of the discriminator identity on CCA addition revealed that purine bases (with a preference for adenosine) are strongly favoured over pyrimidines. Furthermore, depending on the tRNA context, a cytosine discriminator can cause a dramatic number of misincorporations during CCA addition. The data correlate with a high frequency of adenosine residues at the discriminator position observed in vivo. Originally identified as a prominent identity element for aminoacylation, this position represents a likewise important element for efficient and accurate CCA addition.

Complete mitochondrial genomes of the human follicle mites Demodex brevis and D. folliculorum: novel gene arrangement, truncated tRNA genes, and ancient divergence between species

Background

Follicle mites of the genus Demodex are found on a wide diversity of mammals, including humans; surprisingly little is known, however, about the evolution of this association. Additional sequence information promises to facilitate studies of Demodex variation within and between host species. Here we report the complete mitochondrial genome sequences of two species of Demodex known to live on humans—Demodex brevis and D. folliculorum—which are the first such genomes available for any member of the genus. We analyzed these sequences to gain insight into the evolution of mitochondrial genomes within the Acariformes. We also used relaxed molecular clock analyses, based on alignments of mitochondrial proteins, to estimate the time of divergence between these two species.

Results

Both Demodex genomes shared a novel gene order that differs substantially from the ancestral chelicerate pattern, with transfer RNA (tRNA) genes apparently having moved much more often than other genes. Mitochondrial tRNA genes of both species were unusually short, with most of them unable to encode tRNAs that could fold into the canonical cloverleaf structure; indeed, several examples lacked both D- and T-arms. Finally, the high level of sequence divergence observed between these species suggests that these two lineages last shared a common ancestor no more recently than about 87 mya.

Conclusions

Among Acariformes, rearrangements involving tRNA genes tend to occur much more often than those involving other genes. The truncated tRNA genes observed in both Demodex species would seem to require the evolution of extensive tRNA editing capabilities and/or coevolved interacting factors. The molecular machinery necessary for these unusual tRNAs to function might provide an avenue for developing treatments of skin disorders caused by Demodex. The deep divergence time estimated between these two species sets a lower bound on the time that Demodex have been coevolving with their mammalian hosts, and supports the hypothesis that there was an early split within the genus Demodex into species that dwell in different skin microhabitats.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1124) contains supplementary material, which is available to authorized users.

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 genome of the scab mite Psoroptes cuniculi (Arthropoda: Arachnida) provides insights into Acari phylogeny

Background

Limited available sequence information has greatly impeded population genetics, phylogenetics and systematics studies in the subclass Acari (mites and ticks). Mitochondrial (mt) DNA is well known to provide genetic markers for investigations in these areas, but complete mt genomic data have been lacking for many Acari species. Herein, we present the complete mt genome of the scab mite Psoroptes cuniculi.

Methods

P. cuniculi was collected from a naturally infected New Zealand white rabbit from China and identified by morphological criteria. The complete mt genome of P. cuniculi was amplified by PCR and then sequenced. The relationships of this scab mite with selected members of the Acari were assessed by phylogenetic analysis of concatenated amino acid sequence datasets by Bayesian inference (BI), maximum likelihood (ML) and maximum parsimony (MP).

Results

This mt genome (14,247 bp) is circular and consists of 37 genes, including 13 genes for proteins, 22 genes for tRNA, 2 genes for rRNA. The gene arrangement in mt genome of P. cuniculi is the same as those of Dermatophagoides farinae (Pyroglyphidae) and Aleuroglyphus ovatus (Acaridae), but distinct from those of Steganacarus magnus (Steganacaridae) and Panonychus citri (Tetranychidae). Phylogenetic analyses using concatenated amino acid sequences of 12 protein-coding genes, with three different computational algorithms (BI, ML and MP), showed the division of subclass Acari into two superorders, supported the monophylies of the both superorders Parasitiformes and Acariformes; and the three orders Ixodida and Mesostigmata and Astigmata, but rejected the monophyly of the order Prostigmata.

Conclusions

The mt genome of P. cuniculi represents the first mt genome of any member of the family Psoroptidae. Analysis of mt genome sequences in the present study has provided new insights into the phylogenetic relationships among several major lineages of Acari species.

Fast assembly of the mitochondrial genome of a plant parasitic nematode (Meloidogyne graminicola) using next generation sequencing

Little is known about the variations of nematode mitogenomes (mtDNA). Sequencing a complete mtDNA using a PCR approach remains a challenge due to frequent genome reorganizations and low sequence similarities between divergent nematode lineages. Here, a genome skimming approach based on HiSeq sequencing (shotgun) was used to assemble de novo the first complete mtDNA sequence of a root-knot nematode (Meloidogyne graminicola). An AT-rich genome (84.3%) of 20,030 bp was obtained with a mean sequencing depth superior to 300. Thirty-six genes were identified with a semi-automated approach. A comparison with a gene map of the M. javanica mitochondrial genome indicates that the gene order is conserved within this nematode lineage. However, deep genome rearrangements were observed when comparing with other species of the superfamily Hoplolaimoidea. Repeat elements of 111 bp and 94 bp were found in a long non-coding region of 7.5 kb, as similarly reported in M. javanica and M. hapla. This study points out the power of next generation sequencing to produce complete mitochondrial genomes, even without a reference sequence, and possibly opening new avenues for species/race identification, phylogenetics and population genetics of nematodes.

Biological evidence for the world's smallest tRNAs

Due to their function as adapters in translation, tRNA molecules share a common structural organization in all kingdoms and organelles with ribosomal protein biosynthesis. A typical tRNA has a cloverleaf-like secondary structure, consisting of acceptor stem, D-arm, anticodon arm, a variable region, and T-arm, with an average length of 73 nucleotides. In several mitochondrial genomes, however, tRNA genes encode transcripts that show a considerable deviation of this standard, having reduced D- or T-arms or even completely lack one of these elements, resulting in tRNAs as small as 66 nts. An extreme case of such truncations is found in the mitochondria of Enoplea. Here, several tRNA genes are annotated that lack both the D- and the T-arm, suggesting even shorter transcripts with a length of only 42 nts. However, direct evidence for these exceptional tRNAs, which were predicted by purely computational means, has been lacking so far. Here, we demonstrate that several of these miniaturized armless tRNAs consisting only of acceptor- and anticodon-arms are indeed transcribed and correctly processed by non-encoded CCA addition in the mermithid Romanomermis culicivorax. This is the first direct evidence for the existence and functionality of the smallest tRNAs ever identified so far. It opens new possibilities towards exploration/assessment of minimal structural motifs defining a functional tRNA and their evolution.

Arthropod phylogenetics in light of three novel millipede (myriapoda: diplopoda) mitochondrial genomes with comments on the appropriateness of mitochondrial genome sequence data for inferring deep level relationships

BACKGROUND

Arthropods are the most diverse group of eukaryotic organisms, but their phylogenetic relationships are poorly understood. Herein, we describe three mitochondrial genomes representing orders of millipedes for which complete genomes had not been characterized. Newly sequenced genomes are combined with existing data to characterize the protein coding regions of myriapods and to attempt to reconstruct the evolutionary relationships within the Myriapoda and Arthropoda.

RESULTS

The newly sequenced genomes are similar to previously characterized millipede sequences in terms of synteny and length. Unique translocations occurred within the newly sequenced taxa, including one half of the Appalachioria falcifera genome, which is inverted with respect to other millipede genomes. Across myriapods, amino acid conservation levels are highly dependent on the gene region. Additionally, individual loci varied in the level of amino acid conservation. Overall, most gene regions showed low levels of conservation at many sites. Attempts to reconstruct the evolutionary relationships suffered from questionable relationships and low support values. Analyses of phylogenetic informativeness show the lack of signal deep in the trees (i.e., genes evolve too quickly). As a result, the myriapod tree resembles previously published results but lacks convincing support, and, within the arthropod tree, well established groups were recovered as polyphyletic.

CONCLUSIONS

The novel genome sequences described herein provide useful genomic information concerning millipede groups that had not been investigated. Taken together with existing sequences, the variety of compositions and evolution of myriapod mitochondrial genomes are shown to be more complex than previously thought. Unfortunately, the use of mitochondrial protein-coding regions in deep arthropod phylogenetics appears problematic, a result consistent with previously published studies. Lack of phylogenetic signal renders the resulting tree topologies as suspect. As such, these data are likely inappropriate for investigating such ancient relationships.

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.

Improved systematic tRNA gene annotation allows new insights into the evolution of mitochondrial tRNA structures and into the mechanisms of mitochondrial genome rearrangements

Transfer RNAs (tRNAs) are present in all types of cells as well as in organelles. tRNAs of animal mitochondria show a low level of primary sequence conservation and exhibit 'bizarre' secondary structures, lacking complete domains of the common cloverleaf. Such sequences are hard to detect and hence frequently missed in computational analyses and mitochondrial genome annotation. Here, we introduce an automatic annotation procedure for mitochondrial tRNA genes in Metazoa based on sequence and structural information in manually curated covariance models. The method, applied to re-annotate 1876 available metazoan mitochondrial RefSeq genomes, allows to distinguish between remaining functional genes and degrading 'pseudogenes', even at early stages of divergence. The subsequent analysis of a comprehensive set of mitochondrial tRNA genes gives new insights into the evolution of structures of mitochondrial tRNA sequences as well as into the mechanisms of genome rearrangements. We find frequent losses of tRNA genes concentrated in basal Metazoa, frequent independent losses of individual parts of tRNA genes, particularly in Arthropoda, and wide-spread conserved overlaps of tRNAs in opposite reading direction. Direct evidence for several recent Tandem Duplication-Random Loss events is gained, demonstrating that this mechanism has an impact on the appearance of new mitochondrial gene orders.

The mitochondrial genome of Sinentomon erythranum (Arthropoda: Hexapoda: Protura): an example of highly divergent evolution

BACKGROUND

The phylogenetic position of the Protura, traditionally considered the most basal hexapod group, is disputed because it has many unique morphological characters compared with other hexapods. Although mitochondrial genome information has been used extensively in phylogenetic studies, such information is not available for the Protura. This has impeded phylogenetic studies on this taxon, as well as the evolution of the arthropod mitochondrial genome.

RESULTS

In this study, the mitochondrial genome of Sinentomon erythranum was sequenced, as the first proturan species to be reported. The genome contains a number of special features that differ from those of other hexapods and arthropods. As a very small arthropod mitochondrial genome, its 14,491 nucleotides encode 37 typical mitochondrial genes. Compared with other metazoan mtDNA, it has the most biased nucleotide composition with T = 52.4%, an extreme and reversed AT-skew of -0.351 and a GC-skew of 0.350. Two tandemly repeated regions occur in the A+T-rich region, and both could form stable stem-loop structures. Eighteen of the 22 tRNAs are greatly reduced in size with truncated secondary structures. The gene order is novel among available arthropod mitochondrial genomes. Rearrangements have involved in not only small tRNA genes, but also PCGs (protein-coding genes) and ribosome RNA genes. A large block of genes has experienced inversion and another nearby block has been reshuffled, which can be explained by the tandem duplication and random loss model. The most remarkable finding is that trnL2(UUR) is not located between cox1 and cox2 as observed in most hexapod and crustacean groups, but is between rrnL and nad1 as in the ancestral arthropod ground pattern. The "cox1-cox2" pattern was further confirmed in three more representative proturan species. The phylogenetic analyses based on the amino acid sequences of 13 mitochondrial PCGs suggest S. erythranum failed to group with other hexapod groups.

CONCLUSIONS

The mitochondrial genome of S. erythranum shows many different features from other hexapod and arthropod mitochondrial genomes. It underwent highly divergent evolution. The "cox1-cox2" pattern probably represents the ancestral state for all proturan mitogenomes, and suggests a long evolutionary history for the Protura.

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.

Repeated parallel evolution of minimal rRNAs revealed from detailed comparative analysis

The concept of a minimal ribosomal RNA-containing ribosome, a structure with a minimal set of elements capable of providing protein biosynthesis, is essential for understanding this fundamental cellular process. Nematodes and trypanosomes have minimal mitochondrial rRNAs and detailed reconstructions of their secondary structures indicate that certain conserved helices have been lost in these taxa. In contrast, several recent studies on acariform mites have argued that minimal rRNAs may evolve via shortening of secondary structure elements but not the loss of these elements as shown for trypanosomes and nematodes. Based on extensive structural analysis of chelicerate arthropods, we demonstrate that extremely short rRNAs of acariform mites share certain structural modifications with nematodes and trypanosomes: loss of helices of the GTPase region and divergence in the evolutionarily conserved connecting loop between helices H1648 and H1764 of the large subunit rRNA. These highly concerted parallel modifications indicate that minimal rRNAs were generated under the strong selection that favored or tolerated reductions of helices in particular locations while maintaining the functionality of the rRNA molecules throughout evolution. We also discuss potential evolution of minimal rRNAs and atypical transfer RNAs.