The complete mitochondrial genome analysis of Haemaphysalis hystricis Supino, 1897 (Ixodida: Ixodidae) and its phylogenetic implications

In order to study the sequence characteristics, gene order, and codon usage of the mitochondrial genome of Haemaphysalis hystricis, and to explore its phylogenetic relationship, a total of 36 H. hystricis isolated from dogs were used as sample in this study. The mitochondrial genome of a H. hystricis was amplified with several pairs of specific primers by PCR, and was sequenced by first generation sequencing. The mitochondrial genome of H. hystricis was 14,719 bp in size, and it contained 37 genes including 13 protein coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and AT-rich region. Each PCG sequence had different lengths, the sequence longest and shortest gene were nad5 (1,652 bp) and atp8 (155 bp), respectively, among the 13 PCGs. All PCGs used ATN as their initiation codon, 10 of 13 PCGs used TAN as their termination codon, and 3 of which had incomplete termination codon (TA/T). Most of the 22 tRNAs with different sizes could form the classical cloverleaf structures expect for tRNA-Ala, tRNA-Ser1, tRNA-Ser2, and tRNA-Glu, and there were base mismatch (U-U and U-G) in all the 22 tRNAs sequences. Two rRNAs, namely rrnL and rrnS, had different lengths, rrnL located between tRNA-Leu1 and tRNA-Val, and rrnS located between tRNA-Val and tRNA-Ile, respectively. Two AT (D-loop) control areas with different lengths were in the mitochondrial genome, the NCRL was located between tRNA-Leu2 and tRNA-Cys, and the NCRS was located between rrnS and tRNA-Ile. The complete mitochondrial genome sequence of H. hystricis was AT preferences, and the gene order is the same as that of other Haemaphysalis family ticks. However, phylogenetic analysis showed that H. hystricis was most closely related to Haemaphysalis longicornis among the selected ticks. The mitochondrial genome not only enriches the genome database, provides more novel genetic markers for identifying tick species, and studying its molecular epidemiology, population genetics, systematics, but also have implications for the diagnosis, prevention, and control of ticks and tick-borne diseases in animals and humans.

Rapid evolution of mitochondrion-related genes in haplodiploid arthropods

BACKGROUND

Mitochondrial genes and nuclear genes cooperate closely to maintain the functions of mitochondria, especially in the oxidative phosphorylation (OXPHOS) pathway. However, mitochondrial genes among arthropod lineages have dramatic evolutionary rate differences. Haplodiploid arthropods often show fast-evolving mitochondrial genes. One hypothesis predicts that the small effective population size of haplodiploid species could enhance the effect of genetic drift leading to higher substitution rates in mitochondrial and nuclear genes. Alternatively, positive selection or compensatory changes in nuclear OXPHOS genes could lead to the fast-evolving mitochondrial genes. However, due to the limited number of arthropod genomes, the rates of evolution for nuclear genes in haplodiploid species, besides hymenopterans, are largely unknown. To test these hypotheses, we used data from 76 arthropod genomes, including 5 independently evolved haplodiploid lineages, to estimate the evolutionary rates and patterns of gene family turnover of mitochondrial and nuclear genes.

RESULTS

We show that five haplodiploid lineages tested here have fast-evolving mitochondrial genes and fast-evolving nuclear genes related to mitochondrial functions, while nuclear genes not related to mitochondrion showed no significant evolutionary rate differences. Among hymenopterans, bees and ants show faster rates of molecular evolution in mitochondrial genes and mitochondrion-related nuclear genes than sawflies and wasps. With genome data, we also find gene family expansions and contractions in mitochondrion-related genes of bees and ants.

CONCLUSIONS

Our results reject the small population size hypothesis in haplodiploid species. A combination of positive selection and compensatory changes could lead to the observed patterns in haplodiploid species. The elevated evolutionary rates in OXPHOS complex 2 genes of bees and ants suggest a unique evolutionary history of social hymenopterans.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

The complete mitochondrial genome of Eulaelaps huzhuensis (Mesostigmata: Haemogamasidae)

Some mites of the family Haemogamasidae can transmit a variety of zoonotic diseases and have important public health and safety implications. Currently, however, little attention has been paid to molecular data of Haemogamasidae species, limiting our understanding of their evolutionary and phylogenetic relationships. In this study, the complete mitochondrial genome of Eulaelaps huzhuensis was determined for the first time, and its genomic information was analyzed in detail. The mitochondrial genome of E. huzhuensis is 14,872 bp in length with 37 genes and two control regions. The base composition showed a distinct AT preference. Twelve protein-coding genes have a typical ATN as the start codon, and three protein-coding genes have incomplete stop codons. During the folding of tRNA genes, a total of 30 mismatches occurred, and three tRNA genes had an atypical cloverleaf secondary structure. The order of the E. huzhuensis mitochondrial genome arrangement is a new type of rearrangement in Mesostigmata. The phylogenetic analysis confirmed that the family Haemogamasidae is a monophyletic branch and does not belong to a subfamily of the Laelapidae. Our results lay the foundation for subsequent studies on the phylogeny and evolutionary history of the family Haemogamasidae.

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

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

The Mitochondrial Genome of a Freshwater Pelagic Amphipod Macrohectopus branickii Is among the Longest in Metazoa

There are more than 350 species of amphipods (Crustacea) in Lake Baikal, which have emerged predominantly through the course of endemic radiation. This group represents a remarkable model for studying various aspects of evolution, one of which is the evolution of mitochondrial (mt) genome architectures. We sequenced and assembled the mt genome of a pelagic Baikalian amphipod species Macrohectopus branickii. The mt genome is revealed to have an extraordinary length (42,256 bp), deviating significantly from the genomes of other amphipod species and the majority of animals. The mt genome of M. branickii has a unique gene order within amphipods, duplications of the four tRNA genes and Cox2, and a long non-coding region, that makes up about two thirds of the genome’s size. The extension of the mt genome was most likely caused by multiple duplications and inversions of regions harboring ribosomal RNA genes. In this study, we analyzed the patterns of mt genome length changes in amphipods and other animal phyla. Through a statistical analysis, we demonstrated that the variability in the mt genome length may be a characteristic of certain phyla and is primarily conferred by expansions of non-coding regions.

Massive gene rearrangement in mitogenomes of phytoseiid mites

Mitochondrial (mt) gene sequences have been widely used to infer phylogeny in animals. The relative order of mt genes in the mitogenome can also be a useful marker for evolution, but the propensity of mt gene rearrangements vary tremendously among taxa. Ticks and mites in Acari exemplify this trend as some families retain the ancestral arthropod gene order, while others show highly divergent gene orders. Mites in Phytoseiidae, many of which are effective biological control agents, show some of the most divergent gene orders. However, the diversity of mitogenome order within this family is little known. We thus sequenced three mt genomes of phytoseiid mites from two of the most speciose genera: Amblyseius swirskii (Athias-Henriot), Amblyseius tsugawai (Ehara) and Neoseiulus womersleyi (Schicha). We find differences in mt GC skew and nucleotide composition, especially between N. womersleyi and the two Amblyseius species. Each species within Phytoseiidae (including three previously available sequences) present a unique gene order. Phytoseiid mitogenomes show some of the highest numbers of breakpoints when compared to the ancestral arthropod order (up to 33), as well as high numbers of breakpoints within the family (14-30). This suggests a history of massive, ongoing mitogenome rearrangements in the family. Phylogenetic analyses of mt sequences confirm that the degree of gene rearrangements follows phylogenetic relatedness. We discuss possible causes for the high degree of mt gene rearrangement within phytoseiid mites as well as selection in the mt and nuclear genome tied to the independent evolution of many diverse feeding strategies in the family. Finally, we suggest N. womersleyi should be used instead of the synonym Amblyseius pseudolongispinosus.

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.

The Complete Mitochondrial Genome of Aleurocanthus camelliae: Insights into Gene Arrangement and Genome Organization within the Family Aleyrodidae

There are numerous gene rearrangements and transfer RNA gene absences existing in mitochondrial (mt) genomes of Aleyrodidae species. To understand how mt genomes evolved in the family Aleyrodidae, we have sequenced the complete mt genome of Aleurocanthus camelliae and comparatively analyzed all reported whitefly mt genomes. The mt genome of A. camelliae is 15,188 bp long, and consists of 13 protein-coding genes, two rRNA genes, 21 tRNA genes and a putative control region (GenBank: KU761949). The tRNA gene, trnI, has not been observed in this genome. The mt genome has a unique gene order and shares most gene boundaries with Tetraleurodes acaciae. Nineteen of 21 tRNA genes have the conventional cloverleaf shaped secondary structure and two (trnS₁ and trnS₂) lack the dihydrouridine (DHU) arm. Using ARWEN and homologous sequence alignment, we have identified five tRNA genes and revised the annotation for three whitefly mt genomes. This result suggests that most absent genes exist in the genomes and have not been identified, due to be lack of technology and inference sequence. The phylogenetic relationships among 11 whiteflies and Drosophila melanogaster were inferred by maximum likelihood and Bayesian inference methods. Aleurocanthus camelliae and T. acaciae form a sister group, and all three Bemisia tabaci and two Bemisia afer strains gather together. These results are identical to the relationships inferred from gene order. We inferred that gene rearrangement plays an important role in the mt genome evolved from whiteflies.

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.

Complete Mitochondrial Genome of Euseius nicholsi (Ehara et Lee) (Acari:Phytoseiidae)

In this study, the complete mitochondrial genome (mitogenome) sequence of Euseius nicholsi (Acarinae: Phytoseiidae) was determined by long PCR and primer walking methods. The complete mitochondrial genome is 15,561 bp in length and contains 13 protein-coding genes, 2 ribosomal RNA genes, 21 transfer RNA genes. The overall base composition of the genome is A (34.69%), T (43.14%), C (8.06%) and G (14.10%). Stop codon was missed for the COXIII gene in E.nicholsi. The mitogenome would contribute to resolving phylogenetic position and interrelationships of Euseius.

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.

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.

A resourceful genome: updating the functional repertoire and evolutionary role of animal mitochondrial DNAs

Recent data from mitochondrial genomics and proteomics research demonstrate the existence of several atypical mitochondrial protein-coding genes (other than the standard set of 13) and the involvement of mtDNA-encoded proteins in functions other than energy production in several animal species including humans. These results are of considerable importance for evolutionary and cellular biology because they indicate that animal mtDNAs have a larger functional repertoire than previously believed. This review summarizes recent studies on animal species with a non-standard mitochondrial functional repertoire and discusses how these genetic novelties represent promising candidates for studying the role of the mitochondrial genome in speciation.

Parasitic mites of medical and veterinary importance--is there a common research agenda?

There are an estimated 0.5-1 million mite species on earth. Among the many mites that are known to affect humans and animals, only a subset are parasitic but these can cause significant disease. We aim here to provide an overview of the most recent work in this field in order to identify common biological features of these parasites and to inform common strategies for future research. There is a critical need for diagnostic tools to allow for better surveillance and for drugs tailored specifically to the respective parasites. Multi-'omics' approaches represent a logical and timely strategy to identify the appropriate mite molecules. Recent advances in sequencing technology enable us to generate de novo genome sequence data, even from limited DNA resources. Consequently, the field of mite genomics has recently emerged and will now rapidly expand, which is a particular advantage for parasitic mites that cannot be cultured in vitro. Investigations of the microbiota associated with mites will elucidate the link between parasites and pathogens, and define the role of the mite in transmission and pathogenesis. The databases generated will provide the crucial knowledge essential to design novel diagnostic tools, control measures, prophylaxes, drugs and immunotherapies against the mites and associated secondary infections.

Can long-range PCR be used to amplify genetically divergent mitochondrial genomes for comparative phylogenetics? A case study within spiders (Arthropoda: Araneae)

The development of second generation sequencing technology has resulted in the rapid production of large volumes of sequence data for relatively little cost, thereby substantially increasing the quantity of data available for phylogenetic studies. Despite these technological advances, assembling longer sequences, such as that of entire mitochondrial genomes, has not been straightforward. Existing studies have been limited to using only incomplete or nominally intra-specific datasets resulting in a bottleneck between mitogenome amplification and downstream high-throughput sequencing. Here we assess the effectiveness of a wide range of targeted long-range PCR strategies, encapsulating single and dual fragment primer design approaches to provide full mitogenomic coverage within the Araneae (Spiders). Despite extensive rounds of optimisation, full mitochondrial genome PCR amplifications were stochastic in most taxa, although 454 Roche sequencing confirmed the successful amplification of 10 mitochondrial genomes out of the 33 trialled species. The low success rates of amplification using long-Range PCR highlights the difficulties in consistently obtaining genomic amplifications using currently available DNA polymerases optimised for large genomic amplifications and suggests that there may be opportunities for the use of alternative amplification methods.

Transcriptome sequencing and annotation of the predatory mite Metaseiulus occidentalis (Acari: Phytoseiidae): a cautionary tale about possible contamination by prey sequences

Next-generation sequencing was applied to the transcriptome of the phytoseiid Metaseiulus occidentalis to characterize gene expression in all life stages reared under different conditions to optimize the recovery of as many genes as possible. One production and one titration run produced a total of 862,069 reads (average size: 314.87 bp), which generated 255.6 Mbp of sequences on the GS-FLX Titanium sequencing platform. After removal of putative prey sequences 850,543 reads were used in NewBler and PTA assemblies to produce 74,172 non-redundant sequences, including 30,691 contigs and 43,481 singlets with 11,994 contigs consisting of more than 500 bp and 37,278 sequences >300 bp, constituting 48.7 % of all sequences. There were 25,888 hits with the NCBI non-redundant database and 15,376 unique transcripts. There were 26,225 hits with the Ixodes scapularis genome and 6,634 unique transcripts. There were 22,225 hits with the RefSeq of Homo sapiens with 6,465 unique transcripts, and 23,656 hits with the RefSeq of Drosophila melanogaster with 9,216 unique transcripts. Selected ESTs corresponding to genes of interest were analyzed including those related to transposable elements, GPCRs, Sox transcription factors, diapause and foraging behavior, and pesticide resistances. Novel and important genes appear to have been discovered that provide insight into the evolution, biology, and physiology of this important predator of pest mites in agriculture and will be useful in analyzing complete genome sequences of this natural enemy.

Obtaining molecular data for all life stages of Typhlodromus (Typhlodromus) exhilaratus (Mesostigmata: Phytoseiidae): consequences for species identification

Several species of the family Phytoseiidae are known to control mite pests in many crops worldwide. However, biological control success greatly depends on the accurate identification of these predatory mites. Species diagnostics is essentially based on the morphological characters of females. Thus, when only immature stages and/or males are collected, their identification is poorly supported. Molecular tools could be of great help to overcome these difficulties, as molecular sequences are assumed to be identical for the life stage considered. However, one of the essential points is to extract a sufficient DNA amount from a single specimen of immature stages (eggs, protonymphs, deutonymphs) and males (less than 300 μm in length) to amplify and sequence DNA. The markers used were two mitochondrial DNA fragments (12S rRNA and Cytb mtDNA) and the species studied were Typhlodromus (Typhlodromus) exhilaratus and T. (T.) phialatus, two cryptic species, reported to control mite pests in crops of southern Europe and commonly found on the same plants. Despite a low quantity of DNA extracted, particularly for the egg, larva and protonymph stages, DNA was amplified and sequences were obtained from all the life stages considered with the two mtDNA fragments. Sequences from all the developmental stages of T. (T.) exhilaratus were identical and well differentiated from those of its sister-species. However, contaminations were observed especially for eggs and DNA amplified with the Cytb mt marker. Utility of the present results are discussed and protocol improvements are proposed.

Morphological, molecular and cross-breeding analysis of geographic populations of coconut-mite associated predatory mites identified as Neoseiulus baraki: evidence for cryptic species?

Surveys were conducted in Brazil, Benin and Tanzania to collect predatory mites as candidates for control of the coconut mite Aceria guerreronis Keifer, a serious pest of coconut fruits. At all locations surveyed, one of the most dominant predators on infested coconut fruits was identified as Neoseiulus baraki Athias-Henriot, based on morphological similarity with regard to taxonomically relevant characters. However, scrutiny of our own and published descriptions suggests that consistent morphological differences may exist between the Benin population and those from the other geographic origins. In this study, we combined three methods to assess whether these populations belong to one species or a few distinct, yet closely related species. First, multivariate analysis of 32 morphological characters showed that the Benin population differed from the other three populations. Second, DNA sequence analysis based on the mitochondrial cytochrome oxidase subunit I (COI) showed the same difference between these populations. Third, cross-breeding between populations was unsuccessful in all combinations. These data provide evidence for the existence of cryptic species. Subsequent morphological research showed that the Benin population can be distinguished from the others by a new character (not included in the multivariate analysis), viz. the number of teeth on the fixed digit of the female chelicera.

The complete mitochondrial genome of Pseudochauhanea macrorchis (Monogenea: Chauhaneidae) revealed a highly repetitive region and a gene rearrangement hot spot in Polyopisthocotylea

The complete mitochondrial genome of Pseudochauhanea macrorchis was determined and compared with other monogenean mitochondrial genomes from GenBank. The circular genome was 15,031 bp in length and encoded 36 genes (12 protein-coding genes, two ribosomal RNAs, and 22 transfer RNAs) typically found in flatworms. Structures of the mitochondrial genome were mostly concordant with that known for Microcotyle sebastis and Polylabris halichoeres, but also contained two noted features-a gene rearrangement hot spot and the highly repetitive region (HRR) in major non-coding region (NCR). The gene rearrangement hot spot located between the cox3 and nad5 genes, including a cluster of tRNA genes, nad6 gene and one major NCR. The HRR seemed to be a unique feature of the polyopisthocotylean mitochondrial genomes. In conclusion, the present study provided new molecular data for future studies of the comparative mitochondrial genomics and also served as a resource of markers for the studies of species populations and monogenean phylogenetics.

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.

Nucleotide composition of CO1 sequences in Chelicerata (Arthropoda): detecting new mitogenomic rearrangements

Here we study the evolution of nucleotide composition in third codon-positions of CO1 sequences of Chelicerata, using a phylogenetic framework, based on 180 taxa and three markers (CO1, 18S, and 28S rRNA; 5,218 nt). The analyses of nucleotide composition were also extended to all CO1 sequences of Chelicerata found in GenBank (1,701 taxa). The results show that most species of Chelicerata have a positive strand bias in CO1, i.e., in favor of C nucleotides, including all Amblypygi, Palpigradi, Ricinulei, Solifugae, Uropygi, and Xiphosura. However, several taxa show a negative strand bias, i.e., in favor of G nucleotides: all Scorpiones, Opisthothelae spiders and several taxa within Acari, Opiliones, Pseudoscorpiones, and Pycnogonida. Several reversals of strand-specific bias can be attributed to either a rearrangement of the control region or an inversion of a fragment containing the CO1 gene. Key taxa for which sequencing of complete mitochondrial genomes will be necessary to determine the origin and nature of mtDNA rearrangements involved in the reversals are identified. Acari, Opiliones, Pseudoscorpiones, and Pycnogonida were found to show a strong variability in nucleotide composition. In addition, both mitochondrial and nuclear genomes have been affected by higher substitution rates in Acari and Pseudoscorpiones. The results therefore indicate that these two orders are more liable to fix mutations of all types, including base substitutions, indels, and genomic rearrangements.

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.

Integrative taxonomy demonstrates the unexpected synonymy between two predatory mite species: Cydnodromus idaeus and C. picanus (Acari:Phytoseiidae)

Species of the mite family Phytoseiidae are well known as predators of mite pests all over the world. Their identification is thus of great interest for biological control. The specimens examined in this study belong to the species Cydnodromus idaeus (described from Brazil) and C. picanus (described from Chile). They were collected together on the same plant in Argentina. These species are mainly differentiated by the presence/absence of the dorsal solenostomes (gland openings) gd2. Some morphometric differences were observed between the Argentinian specimens and the type material of C. idaeus and C. picanus; however, they were not sufficient to support a species separation. Morphological and molecular analyses with two mtDNA (12S rRNA, CytB) and a nuclear DNA (ITS) marker showed no difference between Argentinian specimens with and without the solenostome gd2. Also, molecular differentiation between these latter specimens, C. picanus and C. idaeus, was very low, suggesting that all these individuals belong to the same species. Biological experiments confirm these results, as 75% of the progeny resulting from females with gd2 had gd2 present whereas all the descendants resulting from the females without gd2 had gd2 present too. The presence/absence of this solenostome seems thus to be variable within a species and cannot be used diagnostically. We conclude that all specimens from Argentina belong to C. idaeus and that C. picanus is a junior synonym of C. idaeus. The present results were unexpected as the presence/absence of solenostomes is used in Phytoseiidae identification keys and as one of a suite of diagnostic characters for discriminating between species. The consequences of these findings for the taxonomy of Phytoseiidae are discussed.

The puzzle of DNA sequences of Phytoseiidae (Acari:Mesostigmata) in the public GenBank database

The public database GenBank is an increasingly important source of sequence data for diagnostic and phylogenetic research; however, not all deposited sequences are necessarily correctly ascribed to a source species. We considered the example of the mite family Phytoseiidae to determine how the corresponding sequences could be accurately exploited. Phytoseiidae mites are well known worldwide for their ability to control certain mite and insect pests. The number of molecular approaches, especially for diagnostic purposes, has increased over the past decade, leading to an increase in the number of sequences registered in the GenBank database. The aim of the present paper was to evaluate the validity of the DNA sequences presently assigned to Phytoseiidae species in this database. Three hundred and fifty-one sequences, corresponding to the four most frequently registered DNA fragments (ITS, COI, Cytb and 12S rRNA), were considered. DNA extraction, amplification and sequencing were performed for the fragments 12S rRNA and ITS for Amblyseius andersoni, A. swirskii, Iphiseius degenerans, Euseius ovalis, E. stipulatus, Neoseiulus cucumeris and Typhlodromus pyri, as some identifications were questionable. Numerous problems were evident based on genetic distance analyses of these sequences. First, nomenclatural problems were encountered, preventing the correct identification of the taxa sequenced in one case. Suspected misidentifications were frequent, stressing the importance of voucher specimen availability. For the 12S rRNA fragment, sequences assigned to three Phytoseiidae species were those of their prey (Astigmata), underlining the care that must be taken when manipulating the DNA of such predators (sterile conditions and specific PCR primers). Finally, sequences of two regions of the COI mtDNA were encountered, leading to alignment problems between sequences of a same gene and same species. These results are discussed in relation to responsibilities of authors in terms of taxon identification and the great utility of open access DNA sequence databases, such as GenBank, for improving taxonomic identifications and advancing scientific research.

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.