Mitochondrial heteroplasmy and the evolution of insecticide resistance: non-Mendelian inheritance in action

Novel protein genes in animal mtDNA: a new sex determination system in freshwater mussels (Bivalvia: Unionoida)?

Mitochondrial (mt) function depends critically on optimal interactions between components encoded by mt and nuclear DNAs. mitochondrial DNA (mtDNA) inheritance (SMI) is thought to have evolved in animal species to maintain mito-nuclear complementarity by preventing the spread of selfish mt elements thus typically rendering mtDNA heteroplasmy evolutionarily ephemeral. Here, we show that mtDNA intraorganismal heteroplasmy can have deterministic underpinnings and persist for hundreds of millions of years. We demonstrate that the only exception to SMI in the animal kingdom, that is, the doubly uniparental mtDNA inheritance system in bivalves, with its three-way interactions among egg mt-, sperm mt- and nucleus-encoded gene products, is tightly associated with the maintenance of separate male and female sexes (dioecy) in freshwater mussels. Specifically, this mother-through-daughter and father-through-son mtDNA inheritance system, containing highly differentiated mt genomes, is found in all dioecious freshwater mussel species. Conversely, all hermaphroditic species lack the paternally transmitted mtDNA (=possess SMI) and have heterogeneous macromutations in the recently discovered, novel protein-coding gene (F-orf) in their maternally transmitted mt genomes. Using immunoelectron microscopy, we have localized the F-open reading frame (ORF) protein, likely involved in specifying separate sexes, in mitochondria and in the nucleus. Our results support the hypothesis that proteins coded by the highly divergent maternally and paternally transmitted mt genomes could be directly involved in sex determination in freshwater mussels. Concomitantly, our study demonstrates novel features for animal mt genomes: the existence of additional, lineage-specific, mtDNA-encoded proteins with functional significance and the involvement of mtDNA-encoded proteins in extra-mt functions. Our results open new avenues for the identification, characterization, and functional analyses of ORFs in the intergenic regions, previously defined as "noncoding," found in a large proportion of animal mt genomes.

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.

Acaricide resistance mechanisms in the two-spotted spider mite Tetranychus urticae and other important Acari: a review

The two-spotted spider mite Tetranychus urticae Koch is one of the economically most important pests in a wide range of outdoor and protected crops worldwide. Its control has been and still is largely based on the use of insecticides and acaricides. However, due to its short life cycle, abundant progeny and arrhenotokous reproduction, it is able to develop resistance to these compounds very rapidly. As a consequence, it has the dubious reputation to be the"most resistant species" in terms of the total number of pesticides to which populations have become resistant, and its control has become problematic in many areas worldwide. Insecticide and acaricide resistance has also been reported in the ectoparasite Sarcoptes scabiei, the causative organism of scabies, and other economically important Acari, such as the Southern cattle tick Rhipicephalus microplus, one of the biggest arthropod threats to livestock, and the parasitic mite Varroa destructor, a major economic burden for beekeepers worldwide. Although resistance research in Acari has not kept pace with that in insects, a number of studies on the molecular mechanisms responsible for the resistant phenotype has been conducted recently. In this review, state-of-the-art information on T. urticae resistance, supplemented with data on other important Acari has been brought together. Considerable attention is given to the underlying resistance mechanisms that have been elucidated at the molecular level. The incidence of bifenazate resistance in T. urticae is expanded as an insecticide resistance evolutionary paradigm in arthropods.

Local adaptation of aboveground herbivores towards plant phenotypes induced by soil biota

BACKGROUND

Soil biota may trigger strong physiological responses in plants and consequently induce distinct phenotypes. Plant phenotype, in turn, has a strong impact on herbivore performance. Here, we tested the hypothesis that aboveground herbivores are able to adapt to plant phenotypes induced by soil biota.

METHODOLOGY AND PRINCIPAL FINDINGS

We bred spider mites for 15 generations on snap beans with three different belowground biotic interactions: (i) no biota (to serve as control), (ii) arbuscular mycorrhizal fungi and (ii) root-feeding nematodes. Subsequently, we conducted a reciprocal selection experiment using these spider mites, which had been kept on the differently treated plants. Belowground treatments induced changes in plant biomass, nutrient composition and water content. No direct chemical defence through cyanogenesis was detected in any of the plant groups. Growth rates of spider mites were higher on the ecotypes on which they were bred for 15 generations, although the statistical significance disappeared for mites from the nematode treatment when corrected for all multiple comparisons.

CONCLUSION/SIGNIFICANCE

These results demonstrate that belowground biota may indeed impose selection on the aboveground insect herbivores mediated by the host plant. The observed adaptation was driven by variable quantitative changes of the different separately studied life history traits (i.e. fecundity, longevity, sex-ratio, time to maturity).

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

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

Milbemectin resistance in Tetranychus urticae (Acari: Tetranychidae): selection, stability and cross-resistance to abamectin

Studies on artificial laboratory selections with milbemectin, stability of milbemectin resistance and possible cross-resistance with abamectin were carried out with Tetranychus urticae Koch to provide basic information for a milbemectin resistance management program. Selections for resistance and susceptibility to milbemectin were performed in a population of T. urticae, collected from a commercial chrysanthemum field in the State of São Paulo, Brazil. After six selections for resistance and five selections for susceptibility, susceptible (S) and resistant (R) strains of T. urticae to milbemectin were obtained. The resistance ratio (R/S) at the LC(50) reached 409-fold value. The stability of milbemectin resistance was also studied under laboratory conditions, using a population with initial frequency of 75% of resistant mites. The frequencies of milbemectin resistance were evaluated monthly for a period of 7 months. In order to observe possible correlation between milbemectin and abamectin resistance, the frequencies of abamectin resistance were also evaluated for that population, during the same period. The frequency of milbemectin resistance decreased from 75 to 14.5%, while the percentage of abamectin resistant mites decreased from 57 to 9.1%, in 7 months. The frequencies of milbemectin and abamectin resistance were also evaluated in 25 field populations of T. urticae, collected from several crops in the State of São Paulo. The frequencies of milbemectin resistance varied from 4.1 to 89.5%, and of abamectin, from 7.0 to 90.5%. A positive and significant correlation was observed between the frequencies of milbemectin and abamectin resistance, indicating positive cross-resistance between these acaricides. The results indicate that abamectin should be avoided for managing milbemectin resistance in T. urticae. This is the first report on milbemectin resistance in T. urticae in Brazil.

Development of acaricide resistance in Pacific spider mite (Tetranychus pacificus) from California vineyards

In recent years, grape growers in California reported failures of acaricides against Tetranychus pacificus McGregor. We collected T. pacificus populations from four vineyards and tested them for resistance to bifenazate, propargite and pyridaben. In addition, we sequenced part of the cytochrome b gene of bifenazate-resistant and -susceptible T. pacificus to test for the presence of mutations reported to confer resistance to the congeneric T. urticae. None of the mutations conferring resistance to bifenazate in T. urticae were present in resistant T. pacificus. Resistance levels ranged from full susceptibility to statistically significant 11-fold resistance to pyridaben, sevenfold resistance to bifenazate and fourfold resistance to propargite compared to a susceptible population. Despite the relatively low levels of resistance detected, we estimated that under the conditions of our study the highest field rates of bifenazate and pyridaben application would cause less than 58 and 66% mortality of adult females in the most resistant populations, respectively. In contrast, field rates of propargite application would cause close to 100% mortality in the least susceptible population. These results highlight a potential link between resistance development and reduced field effectiveness for bifenazate and pyridaben. Finally, T. pacificus may be more tolerant to bifenazate and propargite than T. urticae, since the LC(50) values for the susceptible population of T. pacificus were several times higher than LC(50)'s reported for susceptible T. urticae.

Acetylcholinesterase point mutations in European strains of Tetranychus urticae (Acari: Tetranychidae) resistant to organophosphates

BACKGROUND

In Tetranychus urticae Koch, acetylcholinesterase insensitivity is often involved in organophosphate (OP) and carbamate (CARB) resistance. By combining toxicological, biochemical and molecular data from three reference laboratory and three OP selected strains (OP strains), the AChE1 mutations associated with resistance in T. urticae were characterised.

RESULTS

The resistance ratios of the OP strains varied from 9 to 43 for pirimiphos-methyl, from 78 to 586 for chlorpyrifos, from 8 to 333 for methomyl and from 137 to 4164 for dimethoate. The insecticide concentration needed to inhibit 50% of the AChE1 activity was, in the OP strains, at least 2.7, 55, 58 and 31 times higher for the OP pirimiphos-methyl, chlorpyrifos oxon, paraoxon and omethoate respectively, and 87 times higher for the CARB carbaryl. By comparing the AChE1 sequence, four amino acid substitutions were detected in the OP strains: (1) F331W (Torpedo numbering) in all the three OP strains; (2) T280A found in the three OP strains but not in all clones; (3) G328A, found in two OP strains; (4) A201S found in only one OP strain.

CONCLUSIONS

Four AChE1 mutations were found in resistant strains of T. urticae, and three of them, F331W, G328A and A201S, are possibly involved in resistance to OP and CARB insecticides. Among them, F331W is probably the most important and the most common in T. urticae. It can be easily detected by the diagnostic PCR-RLFP assay developed in this study.

Susceptibility of an organophosphate resistant strain of the two-spotted spider mite (Tetranychus urticae) to mixtures of bifenazate with organophosphate and carbamate insecticides

Bifenazate, a new and frequently used carbazate, is a pro-acaricide which needs to be activated by carboxylesterases. We evaluated the possible antagonism of organophosphate and carbamate insecticides on bifenazate toxicity in Tetranychus urticae applied in mixtures. Two organophosphate resistant strains were used (WI and MR-VL) and several organophosphate (chlorpyrifos, azinphosmethyl and phosmet) and carbamate (carbaryl and methomyl) insecticides were evaluated. Mixing chlorpyrifos with bifenazate decreased bifenazate toxicity in both tested strains. However, in the strain with a higher esterase activity, antagonism decreased after 2 days. Of all other tested chemicals, only methomyl displayed an antagonistic effect 1 day after treatment. These findings indicate that mixing organophosphate and carbamate insecticides with bifenazate may inhibit bifenazate efficacy under field conditions, especially when resistant strains are present.

Identification of pyrethroid resistance associated mutations in the para sodium channel of the two-spotted spider mite Tetranychus urticae (Acari: Tetranychidae)

We investigated pyrethroid resistance mechanisms in Tetranychus urticae strains from Greece. Combined bioassay, biochemical and synergistic data indicated that although P450 mono-oxygenase activities were associated with the trait, target site insensitivity was the major resistance component. A 3.3 kb cDNA fragment of the T. urticae para sodium channel gene encompassing segment 4 of domain II to segment 6 of domain IV was obtained by a degenerate PCR strategy. The T. urticae sequence showed highest identity (56%) to the scabies mite, Sarcoptes scabiei, and was phylogenetically classified within the divergent group of Arachnida. Comparison of resistant and susceptible strains identified the point mutation F1538I in segment 6 of domain III, which is known to confer strong resistance to pyrethroids, along with a second mutation (A1215D) in the intracellular linker connecting domains II and III with an unknown role. Three transcripts were identified corresponding to the k and l alternative exons. The mode of inheritance of resistance was confirmed as incompletely recessive, which is consistent with a target site mechanism for pyrethroids.

Wolbachia induces strong cytoplasmic incompatibility in the predatory bug Macrolophus pygmaeus

Macrolophus pygmaeus is a heteropteran predator that is widely used in European glasshouses for the biological control of whiteflies, aphids, thrips and spider mites. We have demonstrated that the insect is infected with the endosymbiotic bacterium Wolbachia pipientis. Several gene fragments of the endosymbiont were sequenced and subsequently used for phylogenetic analysis, revealing that it belongs to the Wolbachia supergroup B. The endosymbiont was visualized within the ovarioles using immunolocalization. Tetracycline treatments were used to cure M. pygmaeus from its infection. Although a completely cured line could not be obtained by this approach, the application of a constant antibiotic pressure over 13 generations resulted in a line with a significantly reduced Wolbachia concentration. Crosses performed with this tetracycline-treated line revealed that the endosymbiont causes severe cytoplasmic incompatibility. This is the first report of a reproductive effect induced by Wolbachia in an economically important heteropteran predator that may have vital implications for its commercial production and use in biological control.

Genetic and biochemical analysis of a laboratory-selected spirodiclofen-resistant strain of Tetranychus urticae Koch (Acari: Tetranychidae)

BACKGROUND

Spirodiclofen is a selective, non-systemic acaricide from the new chemical class of tetronic acid derivatives. In order to develop strategies to minimise resistance in the field, a laboratory-selected spirodiclofen-resistant strain of the two-spotted spider mite, Tetranychus urticae Koch, was used to determine genetic, toxicological, biochemical and cross-resistance data.

RESULTS

Selecting for spirodiclofen resistance in the laboratory yielded a strain (SR-VP) with a resistance ratio of 274, determined on the larval stage. The egg stage remained far more susceptible. No cross-resistance was found against other established acaricides, except for spiromesifen. Based on synergist experiments and enzyme assays, it appeared that especially P450 monooxygenases, but also esterases and glutathione-S-transferases, could be involved in the metabolic detoxification of spirodiclofen. Genetic analysis showed that the resistance is inherited as an intermediate trait under control of more than one gene.

CONCLUSIONS

Resistance to spirodiclofen exceeded by far the recommended field rate. A good acaricide resistance management programme is necessary to prevent fast resistance build-up in the field. Spirodiclofen can be used in alternation with most established acaricides, except for other tetronic acid derivatives. Without selection pressure, resistance tends to be unstable and can decrease in the presence of susceptible individuals owing to the intermediate, polygenic inheritance mode.

Mutations in the mitochondrial cytochrome b of Tetranychus urticae Koch (Acari: Tetranychidae) confer cross-resistance between bifenazate and acequinocyl

BACKGROUND

Resistance of Tetranychus urticae Koch to bifenazate was recently linked with mutations in the mitochondrial cytochrome b Q(o) pocket, suggesting that bifenazate acts as a Q(o) inhibitor (Q(o)I). Since these mutations might cause cross-resistance to the known acaricidal Q(o)I acequinocyl and fluacrypyrim, resistance levels and inheritance patterns were investigated in several bifenazate-susceptible and bifenazate-resistant strains with different mutations in the cd1 and ef helices aligning the Q(o) pocket.

RESULTS

Cross-resistance to acequinocyl in two bifenazate-resistant strains was shown to be maternally inherited and caused by the combination of two specific mutations in the cytochrome b Q(o) pocket. Although most investigated strains were resistant to fluacrypyrim, resistance was not inherited maternally, but as a monogenic autosomal highly dominant trait. As a consequence, there was no correlation between cytochrome b genotype and fluacrypyrim resistance.

CONCLUSIONS

Although there is no absolute cross-resistance between bifenazate, acequinocyl and fluacrypyrim, some bifenazate resistance mutations confer cross-resistance to acequinocyl. In the light of resistance development and management, high prudence is called for when alternating bifenazate and acequinocyl in the same crop. Maternally inherited cross-resistance between bifenazate and acequinocyl reinforces the likelihood of bifenazate acting as a mitochondrial complex III inhibitor at the Q(o) site.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

The symbiotic mites of some Appalachian Xystodesmidae (Diplopoda:Polydesmida) and the complete mitochondrial genome sequence of the mite Stylochyrus rarior (Berlese) (Acari:Mesostigmata:Ologamasidae)

Millipedes of the family Xystodesmidae (Polydesmida) are often host to several symbiotic mite species, but very little work has been done to identify these acarines or to understand their relationship to the millipedes. In an attempt to better understand these associations, mites found on xystodesmid millipedes, a group for which a species phylogeny has been proposed, were collected in the Appalachian Mountains of Kentucky, Virginia, Tennessee and North Carolina. Mites in the genera Stylochyrus Canestrini & Canestrini, 1882 (Mesostigmata: Ologamasidae) and Schwiebea Oudemans, 1916 (Sarcoptiformes: Acaridae) were prevalent among millipedes in the genera Apheloria Chamberlin, 1921, Appalachioria Marek & Bond, 2006, Boraria Chamberlin, 1943, Brachoria Chamberlin, 1939, Dixioria Chamberlin, 1947, Nannaria Chamberlin, 1918, Pleuroloma Rafinesque, 1820, Prionogonus Shelley, 1982, Rudiloria Causey, 1955 and Sigmoria Chamberlin, 1939. Of the mite taxa collected, the species Stylochyrus rarior (Berlese, 1916) was found on the greatest number of sampled millipede taxa. To enhance future coevolutionary studies of xystodesmid millipedes and their mite symbionts, the complete mitochondrial genome of S. rarior associated with the millipede genus Apheloria (Polydesmida: Xystodesmidae) was sequenced. The genome is 14 899 nucleotides in length, has all the typical genes of an arthropod mitochondrion, differs in gene arrangement from that of the ancestral arthropod, and has a gene order that is unique among mites and ticks. The major difference in S. rarior is the placement of the protein-coding gene nad1, which is positioned between the rRNA gene 12S and the protein-coding gene nad2 (tRNA genes and non-coding regions excluded). There are also two non-coding control regions within this mitochondrial genome.