Draft genome of the honey bee ectoparasitic mite, Tropilaelaps mercedesae, is shaped by the parasitic life history

Establishment and validation of a simple and accurate qPCR detection method for Haemophilus parasuis

As an infectious disease that poses a significant threat to the rapidly growing pig breeding industry, the detection of Haemophilus parasuis (HPS) is often compromised by various interfering substances present in the test sample during quantitative real-time PCR (qPCR). The rapid detection of HPS is important for the isolation of infectious pigs and their treatment. We designed and optimized a rapid qPCR test to detect the INFB gene of HPS in clinical and environmental samples on pig farms. The method was evaluated for its specificity, sensitivity, repeatability, anti-interference capability, and its ability to detect HPS in clinical samples. The results indicated that the method was specific for the detection of HPS when evaluated against pathogens and intestinal probiotics found in pig farms. By using a seven-fold dilution series of the recombinant plasmid DNA in triplicate, it was determined that the lowest limit of detection (LOD) for this method was less than 10 copies/µL. The results of inter-batch and intra-batch repeatability tests showed that the coefficient of variation (CV) was consistently below 1%. Furthermore, the impact of 14 endogenous and exogenous interfering substances on the Ct values detected by the HPS qPCR was found to be less than 5% when compared to the Ct values obtained in the absence of interfering substances. A total of 248 clinical samples were analyzed using the HPS qPCR, commercial kits, and corresponding national standards, yielding positive rates of 9.27%, 6.05%, and 9.27%, respectively. Notably, the positive and negative percent agreement between the detection method developed in this research and the national standard was 100%. These findings demonstrate that the established detection method is suitable for epidemiological research on HPS and for diagnosing clinical samples containing interfering substances, thereby providing essential technical support for the prevention and control of HPS.

High-quality reference genome of predatory mite Neoseiulus californicus McGregor (Acari: Phytoseiidae) provides insights into its biological traits and potential RNAi off-target effects

BACKGROUND

Neoseiulus californicus is a predatory mite that can control various spider mites and other small arthropods. Despite its acknowledged effectiveness in the natural enemy market, a crucial knowledge gap exists in understanding the genomic features related to its predatory traits and adaptation. With the increasing emphasis on modern pest management strategies and dynamic environmental changes in plant production trends, constructing a reliable genomic resource for N. californicus becomes imperative.

RESULTS

In this study, we provided a high-quality genome assembly and annotation of N. californicus, with a size of 188.43 Mb and 12 946 predicted genes. We identified genomic features and traits related to its detoxification, stress response, sensory system, mobility, secretory toxins, digestive enzymes and horizontally transferred genes (HGTs) by comparative genomics. We also predicted the potential off-target effects of double-stranded RNA (dsRNA) that could be used to control spider mites, based on public data and the N. californicus genome.

CONCLUSION

Our study provided a valuable genomic resource for a commercialized predatory mite, offering useful insights for the design of integrated pest management (IPM) strategies in the new era. Further studies are needed to explore the functional roles of key gene families in predatory mites, as well as their interactions with their prey and IPM strategies. © 2024 Society of Chemical Industry.

Using minor variant genomes and machine learning to study the genome biology of SARS-CoV-2 over time

In infected individuals, viruses are present as a population consisting of dominant and minor variant genomes. Most databases contain information on the dominant genome sequence. Since the emergence of SARS-CoV-2 in late 2019, variants have been selected that are more transmissible and capable of partial immune escape. Currently, models for projecting the evolution of SARS-CoV-2 are based on using dominant genome sequences to forecast whether a known mutation will be prevalent in the future. However, novel variants of SARS-CoV-2 (and other viruses) are driven by evolutionary pressure acting on minor variant genomes, which then become dominant and form a potential next wave of infection. In this study, sequencing data from 96 209 patients, sampled over a 3-year period, were used to analyse patterns of minor variant genomes. These data were used to develop unsupervised machine learning clusters to identify amino acids that had a greater potential for mutation than others in the Spike protein. Being able to identify amino acids that may be present in future variants would better inform the design of longer-lived medical countermeasures and allow a risk-based evaluation of viral properties, including assessment of transmissibility and immune escape, thus providing candidates with early warning signals for when a new variant of SARS-CoV-2 emerges.

Life-history stage determines the diet of ectoparasitic mites on their honey bee hosts

Ectoparasitic mites of the genera Varroa and Tropilaelaps have evolved to exclusively exploit honey bees as food sources during alternating dispersal and reproductive life history stages. Here we show that the primary food source utilized by Varroa destructor depends on the host life history stage. While feeding on adult bees, dispersing V. destructor feed on the abdominal membranes to access to the fat body as reported previously. However, when V. destructor feed on honey bee pupae during their reproductive stage, they primarily consume hemolymph, indicated by wound analysis, preferential transfer of biostains, and a proteomic comparison between parasite and host tissues. Biostaining and proteomic results were paralleled by corresponding findings in Tropilaelaps mercedesae, a mite that only feeds on brood and has a strongly reduced dispersal stage. Metabolomic profiling of V. destructor corroborates differences between the diet of the dispersing adults and reproductive foundresses. The proteome and metabolome differences between reproductive and dispersing V. destructor suggest that the hemolymph diet coincides with amino acid metabolism and protein synthesis in the foundresses while the metabolism of non-reproductive adults is tuned to lipid metabolism. Thus, we demonstrate within-host dietary specialization of ectoparasitic mites that coincides with life history of hosts and parasites.

A review of the molecular mechanisms of acaricide resistance in mites and ticks

The Arachnida subclass of Acari comprises many harmful pests that threaten agriculture as well as animal health, including herbivorous spider mites, the bee parasite Varroa, the poultry mite Dermanyssus and several species of ticks. Especially in agriculture, acaricides are often used intensively to minimize the damage they inflict, promoting the development of resistance. Beneficial predatory mites used in biological control are also subjected to acaricide selection in the field. The development and use of new genetic and genomic tools such as genome and transcriptome sequencing, bulked segregant analysis (QTL mapping), and reverse genetics via RNAi or CRISPR/Cas9, have greatly increased our understanding of the molecular genetic mechanisms of resistance in Acari, especially in the spider mite Tetranychus urticae which emerged as a model species. These new techniques allowed to uncover and validate new resistance mutations in a larger range of species. In addition, they provided an impetus to start elucidating more challenging questions on mechanisms of gene regulation of detoxification associated with resistance.

Do mites eat and run? A systematic review of feeding and dispersal strategies

Dispersal is an important process affecting the survival of organisms and the structure and dynamics of communities and ecosystems in space and time. It is a multiphase phenomenon influenced by many internal and external factors. Dispersal syndromes can be complicated, but they are vital to our knowledge of the biology of any organism. We analysed dispersal ability in mites (Acariformes and Parasitiformes), a highly diverse group of wingless arthropods, taking into consideration various modes of dispersal, feeding strategies, body size and the number of articles published for each species. Based on 174 articles summarized for this study, it appears that mites are opportunistic when it comes to dispersal, regardless of their feeding habits, and are often able to adopt several different strategies as needs arise. Moreover, we find a significant positive relationship between the amount of research effort that was put into studying a given species and the number of modes of dispersal that were described. The most salient conclusion to be drawn from this positive correlation is that additional studies are needed, especially on a broader set of mite taxa, until the aforementioned correlation is no longer demonstrably significant.

Identification of potential sex determination genes and functional analyses in Neoseiulus californicus under prey stress

BACKGROUND

Phytoseiid mites are important natural enemies of spider mites. Sex-determination mechanism are important basic scientific issues in the reproduction and evolution of predatory mites. Clarifying sex-determination mechanism may provide reference for exploring genetic approach to have the phytoseiid mites produce more female offspring, which could improve their effectiveness as a biological control agent.

RESULTS

We used transcriptome sequencing to identify and characterize 20 putative sex-determination genes in the phytoseiid mite Neoseiulus californicus, a species with uncommon pseudo-arrhenotoky, including doublesex-like (dsx1-like), transformer-2 (tra-2), intersex (ix), and fruitless-like (BTB2). A significant negative correlation was found between prey stress and offspring sex ratio. But the most genes identified showed no difference in expression between the groups with lowest and highest female offspring ratios. The hatching rate and sex ratio of female offspring were reduced when the ix gene was silenced, and the oviposition days and fecundity were reduced when the BTB2 gene was silenced. The fecundity was reduced when the tra2 gene was silenced and the snf gene is essential for oviposition in female. There was no effect on reproduction and female sex determination when silencing the dsx1-like and dsx2-like gene.

CONCLUSION

The genes BTB2, tra2 and snf are involved in oviposition, and ix may be involved in female sex determination and egg formation in Neoseiulus californicus. The results are conductive to further understanding molecular regulatory mechanism of sex determination in predatory mites and may provide a reference for better use of this predatory by producing more females. © 2022 Society of Chemical Industry.

Effects of Niemann-Pick type C2 (NPC2) gene transcripts silencing on behavior of Varroa destructor and molecular changes in the putative olfactory gene networks

To understand the role of two Niemann-Pick type C2 (NPC2) transcripts, Vd40090 (NP1) and Vd74517 (NP5), in the chemosensing pathway of Varroa destructor, we evaluated the impact of NP5 silencing on mites behavior and compared the effect of silencing of either transcripts on the interaction between chemosensory transcripts. In contrast to silencing NP1, which reduced feeding and reproduction by the mite (Nganso et al., 2021), silencing of NP5 reduced significantly the host reaching ability, but it did not affect the feeding on nurse bee. However, silencing of either transcript changed dramatically the co-expression patterns among the putative chemosensory genes, binding proteins and receptors. The results suggest the role of gustatory receptors in the detection of long-range chemical cues in the chemosensory cascade of the Varroa mite.

Insights into gene manipulation techniques for Acari functional genomics

Functional genomics is an essential tool for elucidating the structure and function of genes in any living organism. Here, we review the use of different gene manipulation techniques in functional genomics of Acari (mites and ticks). Some of these Acari species inflict severe economic losses to managed crops and health problems to humans, wild and domestic animals, but many also provide important ecosystem services worldwide. Currently, RNA interference (RNAi) is the leading gene expression manipulation tool followed by gene editing via the bacterial type II Clustered Regularly Interspaced Short Palindromic Repeats and associated protein 9 system (CRISPR-Cas9). Whilst RNAi, via siRNA, does not always lead to expected outcomes, the exploitations of the CRISPR systems in Acari are still in their infancy and are limited only to CRISP/Cas9 to date. In this review, we discuss the advantages and disadvantages of RNAi and CRISPR-Cas9 and the technical challenges associated with their exploitations. We also compare the biochemical machinery of RNAi and CRISPR-Cas9 technologies. We highlight some potential solutions for experimental optimization of each mechanism in gene function studies. The potential benefits of adopting various CRISPR-Cas9 systems for expanding on functional genomics experiments in Acari are also discussed.

De novo transcriptome sequencing of the northern fowl mite, Ornithonyssus sylviarum, shed light on parasitiform poultry mites evolution and its chemoreceptor repertoires

The northern fowl mite (NFM), Ornithonyssus sylviarum, and the poultry red mite (PRM), Dermanyssus gallinae, are the most serious pests of poultry, both of which have an expanding global prevalence. Research on NFM has been constrained by a lack of genomic and transcriptomic data. Here, we report and analyze the first global transcriptome data across all mite live stages and sexes. A total of 28,999 unigenes were assembled, of which 19,750 (68.10%) were annotated using seven functional databases. The biological function of these unigenes was classified using the GO, KOG, and KEGG databases. To gain insight into the chemosensory receptor-based system of parasitiform mites, we furthermore assessed the gene repertoire of gustatory receptors (GRs) and ionotropic receptors (IRs), both of which encode putative ligand-gated ion channel proteins. While these receptors are well characterized in insect model species, our understanding of chemosensory detection in mites and ticks is in its infancy. To address this paucity of data, we identified 9 IR/iGluRs and 2 GRs genes by analyzing transcriptome data in the NFM, while 9 GRs and 41 IR/iGluRs genes were annotated in the PRM genome. Taken together, the transcriptomic and genomic annotation of these two species provide a valuable reference for studies of parasitiform mites and also help to understand how chemosensory gene family expansion/contraction events may have been reshaped by an obligate parasitic lifestyle compared with their free-living closest relatives. Future studies should include additional species to validate this observation and functional characterization of the identified proteins as a step forward in identifying tools for controlling these poultry pests.

Towards disrupting Varroa -honey bee chemosensing: A focus on a Niemann-Pick type C2 transcript

We focused our study on the 12 recently identified putative odorant carrier proteins in the ectoparasitic mite, Varroa destructor. Here we show, via an exclusion of the chemosensory appendages (forelegs and gnathosoma) that transcripts of five of the 12 genes were significantly lower, suggesting that they are likely involved in carrying host volatiles. Specifically, three transcripts were found to be foreleg-specific while the other two transcripts were expressed in both the forelegs and gnathosoma. We focused on one of the highly expressed and foreleg-specific transcript Vd40090, which encodes a Niemann-Pick disease protein type C2 (NPC2) protein. Effects of dsRNA-mediated silencing of Vd40090 were first measured by quantifying the transcript levels of genes that encode other putative odorant carrier proteins as well as reproduction related proteins. In addition, the impact of silencing on mites behaviour and survival was tested. Silencing of Vd40090 effectively disrupted Varroa host selection, acceptance and feeding and significantly impaired the expression of genes that regulate its reproduction in brood cells, resulting in reduced reproduction and survival.

Comparative analysis of mite genomes reveals positive selection for diet adaptation

Diet is a powerful evolutionary force for species adaptation and diversification. Acari is one of the most abundant clades of Arachnida, exhibiting diverse dietary types, while the underlying genetic adaptive mechanisms are not fully understood. Based on comparative analyses of 15 Acari genomes, we found genetic bases for three specialized diets. Herbivores experienced stronger selection pressure than other groups; the olfactory genes and gene families involving metabolizing toxins showed strong adaptive signals. Genes and gene families related to anticoagulation, detoxification, and haemoglobin digestion were found to be under strong selection pressure or significantly expanded in the blood-feeding species. Lipid metabolism genes have a faster evolutionary rate and been subjected to greater selection pressures in fat-feeding species; one positively selected site in the fatty-acid amide hydrolases 2 gene was identified. Our research provides a new perspective for the evolution of Acari and offers potential target loci for novel pesticide development.

A genome-wide screening for RNAi pathway proteins in Acari

BACKGROUND

RNA interference (RNAi) is a highly conserved, sequence-specific gene silencing mechanism present in Eukaryotes. Three RNAi pathways are known, namely micro-RNA (miRNA), piwi-interacting RNA (piRNA) and short interfering RNA (siRNA). However, little knowledge exists about the proteins involved in these pathways in Acari. Moreover, variable successes has been obtained in gene knockdown via siRNA pathway in their functional genomics and management. We hypothesized that the clue may be in the variability of the composition and the efficacy of siRNA machinery among Acari.

RESULTS

Both comparative genomic analyses and domain annotation suggest that all the analyzed species have homologs of putative core proteins that mediate cleaving of targeted genes via the three RNAi pathways. We identified putative homologs of Caenorhabditis elegans RNA-dependent RNA polymerase (RdRP) protein in all species though no secondary Argonaute homologs that operate with this protein in siRNA amplification mechanism were found, suggesting that the siRNA amplification mechanism present in Acari may be distinct from that described in C. elegans. Moreover, the genomes of these species do not encode homologs of C. elegans systemic RNAi defective-1 (Sid-1) protein that mediate silencing of the mRNA target throughout the treated organisms suggesting that the phenomena of systemic RNAi that has been reported in some Acari species probably occur through a different mechanism. However, homologs of putative RNAi spreading defective-3 (Rsd-3) protein and scavenger receptors namely Eater and SR-CI that mediate endocytosis cellular update of dsRNA in C. elegans and Drosophila melanogaster were found in Acari genomes. This result suggests that cellular dsRNA uptake in Acari is endocytosis-dependent. Detailed phylogenetic analyses of core RNAi pathway proteins in the studied species revealed that their evolution is compatible with the proposed monophyletic evolution of this group.

CONCLUSIONS

Our analyses have revealed the potential activity of all three pathways in Acari. Still, much experimental work remains to be done to confirm the mechanisms behind these pathways in particular those that govern systemic/parental RNAi and siRNA amplification in Acari. Disclosure of these mechanisms will facilitate the development of new and specific management tools for the harmful species and enrichment of the beneficial species.

Honey Bee Suppresses the Parasitic Mite Vitellogenin by Antimicrobial Peptide

The negative effects of honey bee parasitic mites and deformed wing virus (DWV) on honey bee and colony health have been well characterized. However, the relationship between DWV and mites, particularly viral replication inside the mites, remains unclear. Furthermore, the physiological outcomes of honey bee immune responses stimulated by DWV and the mite to the host (honey bee) and perhaps the pathogen/parasite (DWV/mite) are not yet understood. To answer these questions, we studied the tripartite interactions between the honey bee, Tropilaelaps mercedesae, and DWV as the model. T. mercedesae functioned as a vector for DWV without supporting active viral replication. Thus, DWV negligibly affected mite fitness. Mite infestation induced mRNA expression of antimicrobial peptides (AMPs), Defensin-1 and Hymenoptaecin, which correlated with DWV copy number in honey bee pupae and mite feeding, respectively. Feeding T. mercedesae with fruit fly S2 cells heterologously expressing honey bee Hymenoptaecin significantly downregulated mite Vitellogenin expression, indicating that the honey bee AMP manipulates mite reproduction upon feeding on bee. Our results provide insights into the mechanism of DWV transmission by the honey bee parasitic mite to the host, and the novel role of AMP in defending against mite infestation.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Characterisation of three novel β-1,3 glucanases from the medically important house dust mite Dermatophagoides pteronyssinus (airmid)

The European house dust mite, Dermatophagoides pteronyssinus is a major source of airborne allergens worldwide and is found in half of European homes. Interactions between microbes and house dust mites (HDM) are considered important factors that allow them to persist in the home. Laboratory studies indicate the European HDM, D. pteronyssinus is a mycophagous mite, capable of utilising a variety of fungi for nutrients, however specific mycolytic digestive enzymes are unknown. Our previous work identified a number of putative glycosyl hydrolases present in the predicted proteome of D. pteronyssinus airmid and validated the expression of 42 of these. Of note, three GH16 proteins with predicted β-1,3 glucanase activity were found to be consistently present in the mite body and excretome. Here, we performed an extensive bioinformatic, proteomic and biochemical study to characterize three-novel β-1,3 glucanases from this medically important house dust mite. The genes encoding novel β-1,3 glucanases designated Glu1, Glu2 and Glu3 were identified in D. pteronyssinus airmid, each exhibited more than 59% amino acid identity to one another. These enzymes are encoded by Glu genes present in a tri-gene cluster and protein homologs are found in other acari. The patchy phyletic distribution of Glu proteins means their evolutionary history remains elusive, however horizontal gene transfer cannot be completely excluded. Recombinant Glu1 and Glu2 exhibit hydrolytic activity toward laminarin, pachyman and barley glucan. Excreted β-1,3 glucanase activity was increased in response to D. pteronyssinus airmid feeding on baker's yeast. Active β-1,3 glucanases are expressed and excreted in the faeces of D. pteronyssinus airmid indicating they are digestive enzymes capable of breaking down β-1,3 glucans of fungi present in house dust.

Aquatic Hemiptera in Southwest Cameroon: Biodiversity of Potential Reservoirs of Mycobacterium ulcerans and Multiple Wolbachia Sequence Types Revealed by Metagenomics

Buruli ulcer (BU), caused by Mycobacterium ulcerans, is a neglected tropical disease associated with freshwater habitats. A variety of limnic organisms harbor this pathogen, including aquatic bugs (Hemiptera: Heteroptera), which have been hypothesized to be epidemiologically important reservoirs. Aquatic Hemiptera exhibit high levels of diversity in the tropics, but species identification remains challenging. In this study, we collected aquatic bugs from emerging foci of BU in the Southwest Region of Cameroon, which were identified using morphological and molecular methods. The bugs were screened for mycobacterial DNA and a selection of 20 mycobacteria-positive specimens from the families Gerridae and Veliidae were subjected to next-generation sequencing. Only one individual revealed putative M. ulcerans DNA, but all specimens contained sequences from the widespread alpha-proteobacterial symbiont, Wolbachia. Phylogenetic analysis placed the Wolbachia sequences into supergroups A, B, and F. Circularized mitogenomes were obtained for seven gerrids and two veliids, the first from these families for the African continent. This study suggests that aquatic Hemiptera may have a minor role (if any) in the spread of BU in Southwest Cameroon. Our metagenomic analysis provides new insights into the incursion of Wolbachia into aquatic environments and generated valuable resources to aid molecular taxonomic studies of aquatic Hemiptera.

High Resolution Analysis of Respiratory Syncytial Virus Infection In Vivo

Human respiratory syncytial virus (HRSV) is a major cause of pediatric infection and also causes disease in the elderly and those with underlying respiratory problems. There is no vaccine for HRSV and anti-viral therapeutics are not broadly applicable. To investigate the effect of HRSV biology in children, nasopharyngeal aspirates were taken from children with different viral loads and a combined high throughput RNAseq and label free quantitative proteomics approach was used to characterize the nucleic acid and proteins in these samples. HRSV proteins were identified in the nasopharyngeal aspirates from infected children, and their abundance correlated with viral load (Ct value), confirming HRSV infection. Analysis of the HRSV genome indicated that the children were infected with sub-group A virus and that minor variants in nucleotide frequency occurred in discrete clusters along the HRSV genome, and within a patient clustered distinctly within the glycoprotein gene. Data from the samples were binned into four groups; no-HRSV infection (control), high viral load (Ct < 20), medium viral load (Ct = 20-25), and low viral load (Ct > 25). Cellular proteins associated with the anti-viral response (e.g., ISG15) were identified in the nasopharyngeal aspirates and their abundance was correlated with viral load. These combined approaches have not been used before to study HRSV biology in vivo and can be readily applied to the study the variation of virus host interactions.

Divergent evolutionary trajectories following speciation in two ectoparasitic honey bee mites

Multispecies host-parasite evolution is common, but how parasites evolve after speciating remains poorly understood. Shared evolutionary history and physiology may propel species along similar evolutionary trajectories whereas pursuing different strategies can reduce competition. We test these scenarios in the economically important association between honey bees and ectoparasitic mites by sequencing the genomes of the sister mite species Varroa destructor and Varroa jacobsoni. These genomes were closely related, with 99.7% sequence identity. Among the 9,628 orthologous genes, 4.8% showed signs of positive selection in at least one species. Divergent selective trajectories were discovered in conserved chemosensory gene families (IGR, SNMP), and Halloween genes (CYP) involved in moulting and reproduction. However, there was little overlap in these gene sets and associated GO terms, indicating different selective regimes operating on each of the parasites. Based on our findings, we suggest that species-specific strategies may be needed to combat evolving parasite communities.

Feeding by Tropilaelaps mercedesae on pre- and post-capped brood increases damage to Apis mellifera colonies

Tropilaelaps mercedesae parasitism can cause Apis mellifera colony mortality in Asia. Here, we report for the first time that tropilaelaps mites feed on both pre- and post-capped stages of honey bees. Feeding on pre-capped brood may extend their survival outside capped brood cells, especially in areas where brood production is year-round. In this study, we examined the types of injury inflicted by tropilaelaps mites on different stages of honey bees, the survival of adult honey bees, and level of honey bee viruses in 4th instar larvae and prepupae. The injuries inflicted on different developing honey bee stages were visualised by staining with trypan blue. Among pre-capped stages, 4th instar larvae sustained the highest number of wounds (4.6 ± 0.5/larva) while 2nd-3rd larval instars had at least two wounds. Consequently, wounds were evident on uninfested capped brood (5th-6th instar larvae = 3.91 ± 0.64 wounds; prepupae = 5.25 ± 0.73 wounds). Tropilaelaps mite infestations resulted in 3.4- and 6-fold increases in the number of wounds in 5th-6th instar larvae and prepupae as compared to uninfested capped brood, respectively. When wound-inflicted prepupae metamorphosed to white-eyed pupae, all wound scars disappeared with the exuviae. This healing of wounds contributed to the reduction of the number of wounds (≤10) observed on the different pupal stages. Transmission of mite-borne virus such as Deformed Wing Virus (DWV) was also enhanced by mites feeding on early larval stages. DWV and Black Queen Cell Virus (BQCV) were detected in all 4th instar larvae and prepupae analysed. However, viral levels were more pronounced in scarred 4th instar larvae and infested prepupae. The remarkably high numbers of wounds and viral load on scarred or infested developing honey bees may have caused significant weight loss and extensive injuries observed on the abdomen, wings, legs, proboscis and antennae of adult honey bees. Together, the survival of infested honey bees was significantly compromised. This study demonstrates the ability of tropilaelaps mites to inflict profound damage on A. mellifera hosts. Effective management approaches need to be developed to mitigate tropilaelaps mite problems.

Varroa chemosensory proteins: some are conserved across Arthropoda but others are arachnid specific

The tight synchronization between the life cycle of the obligatory parasitic mite Varroa destructor (Varroa) and its host, the honeybee, is mediated by honeybee chemical stimuli. These stimuli are mainly perceived by a pit organ located on the distal part of the mite's foreleg. In the present study, we searched for Varroa chemosensory molecular components by comparing transcriptomic and proteomic profiles between forelegs from different physiological stages, and rear legs. In general, a comparative transcriptomic analysis showed a clear separation of the expression profiles between the rear legs and the three groups of forelegs (phoretic, reproductive and tray-collected mites). Most of the differentially expressed transcripts and proteins in the mite's foreleg were previously uncharacterized. Using a conserved domain approach, we identified 45 transcripts with known chemosensory domains belonging to seven chemosensory protein families, of which 14 were significantly upregulated in the mite's forelegs when compared to rear legs. These are soluble and membrane bound proteins, including the somewhat ignored receptors of degenerin/epithelial Na+ channels and transient receptor potentials. Phylogenetic clustering and expression profiles of the putative chemosensory proteins suggest their role in chemosensation and shed light on the evolution of these proteins in Chelicerata.

Honey Bee Parasitic Mite Contains the Sensilla-Rich Sensory Organ on the Foreleg Tarsus Expressing Ionotropic Receptors With Conserved Functions

Honey bee parasitic mites (Tropilaelaps mercedesae and Varroa destructor) detect temperature, humidity, and odor but the underlying sensory mechanisms are poorly understood. To uncover how T. mercedesae responds to environmental stimuli inside a hive, we first identified the sensilla-rich sensory organ on the foreleg tarsus. The organ appeared to correspond to Haller's organ in ticks and contained four types of sensilla, which may respond to different stimuli based on their morphology. We searched for differentially expressed genes between the forelegs and hindlegs to identify mRNAs potentially associated with the sensory organ. The forelegs were enriched with mRNAs encoding sensory proteins such as ionotropic receptors (IRs) and gustatory receptors, as well as proteins involved in ciliary transport. We also found that T. mercedesae IR25a and IR93a were capable of rescuing temperature and humidity preference defects in Drosophila melanogaster IR25a and IR93a mutants. These results demonstrate that the structures and physiological functions of ancient IRs have been conserved during arthropod evolution. Our study provides insight into the sensory mechanisms of honey bee parasitic mites, as well as potential targets for methods to control the most serious honey bee pest.

Genomes of trombidid mites reveal novel predicted allergens and laterally transferred genes associated with secondary metabolism

Background

Trombidid mites have a unique life cycle in which only the larval stage is ectoparasitic. In the superfamily Trombiculoidea ("chiggers"), the larvae feed preferentially on vertebrates, including humans. Species in the genus Leptotrombidium are vectors of a potentially fatal bacterial infection, scrub typhus, that affects 1 million people annually. Moreover, chiggers can cause pruritic dermatitis (trombiculiasis) in humans and domesticated animals. In the Trombidioidea (velvet mites), the larvae feed on other arthropods and are potential biological control agents for agricultural pests. Here, we present the first trombidid mites genomes, obtained both for a chigger, Leptotrombidium deliense, and for a velvet mite, Dinothrombium tinctorium.

Results

Sequencing was performed using Illumina technology. A 180 Mb draft assembly for D. tinctorium was generated from two paired-end and one mate-pair library using a single adult specimen. For L. deliense, a lower-coverage draft assembly (117 Mb) was obtained using pooled, engorged larvae with a single paired-end library. Remarkably, both genomes exhibited evidence of ancient lateral gene transfer from soil-derived bacteria or fungi. The transferred genes confer functions that are rare in animals, including terpene and carotenoid synthesis. Thirty-seven allergenic protein families were predicted in the L. deliense genome, of which nine were unique. Preliminary proteomic analyses identified several of these putative allergens in larvae.

Conclusions

Trombidid mite genomes appear to be more dynamic than those of other acariform mites. A priority for future research is to determine the biological function of terpene synthesis in this taxon and its potential for exploitation in disease control.

Multiple combinations of RDL subunits diversify the repertoire of GABA receptors in the honey bee parasite Varroa destructor

In insects, γ-aminobutyric acid (GABA) is the major inhibitory neurotransmitter, and GABA-gated ion channels are the target of different classes of insecticides, including fipronil. We report here the cloning of six subunits (four RDL, one LCCH3, and one GRD) that constitute the repertoire of the GABA-gated ion channel family of the Varroa mite (Varroa destructor), a honey bee ectoparasite. We also isolated a truncated GRD subunit with a premature stop codon. We found that when expressed in Xenopus laevis oocytes, three of the four RDL subunits (VdesRDL1, VdesRDL2, and VdesRDL3) formed functional, homomultimeric anionic receptors, whereas GRD and LCCH3 produced heteromultimeric cationic receptors. These receptors displayed specific sensitivities toward GABA and fipronil, and VdesRDL1 was the most resistant to the insecticide. We identified specific residues in the VdesRDL1 pore-lining region that explain its high resistance to fipronil. VdesRDL4 did not form a functional receptor when expressed alone, but it assembled with VdesRDL1 to form a heteromultimeric receptor with properties distinct from those of the VdesRDL1 homomultimeric receptor. Moreover, VdesRDL1 physically interacted with VdesRDL3, generating a heteromultimeric receptor combining properties of both subunits. On the other hand, we did not detect any functional interaction between VdesLCCH3 and the VdesRDL subunits, an observation that differed from what was previously reported for Drosophila melanogaster In conclusion, this study provides insights relevant to improve our understanding of the precise role of GABAergic signaling in insects and new tools for the development of Varroa mite-specific insecticidal agents that do not harm honey bees.

Transcriptome profiling of the honeybee parasite Varroa destructor provides new biological insights into the mite adult life cycle

Background

The parasite Varroa destructor represents a significant threat to honeybee colonies. Indeed, development of Varroa infestation within colonies, if left untreated, often leads to the death of the colony. Although its impact on bees has been extensively studied, less is known about its biology and the functional processes governing its adult life cycle and adaptation to its host. We therefore developed a full life cycle transcriptomic catalogue in adult Varroa females and included pairwise comparisons with males, artificially-reared and non-reproducing females (10 life cycle stages and conditions in total).

Results

Extensive remodeling of the Varroa transcriptome was observed, with an upregulation of energetic and chitin metabolic processes during the initial and final phases of the life cycle (e.g. phoretic and post-oviposition stages), whereas during reproductive stages in brood cells genes showing functions related to transcriptional regulation were overexpressed. Several neurotransmitter and neuropeptide receptors involved in behavioural regulation, as well as active compounds of salivary glands, were also expressed at a higher level outside the reproductive stages. No difference was detected between artificially-reared phoretic females and their counterparts in colonies, or between females who failed to reproduce and females who successfully reproduced, indicating that phoretic individuals can be reared outside host colonies without impacting their physiology and that mechanisms underlying reproductive failure occur before oogenesis.

Conclusions

We discuss how these new findings reveal the remarkable adaptation of Varroa to its host biology and notably to the switch from living on adults to reproducing in sealed brood cells. By spanning the entire adult life cycle, our work captures the dynamic changes in the parasite gene expression and serves as a unique resource for deciphering Varroa biology and identifying new targets for mite control.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4668-z) contains supplementary material, which is available to authorized users.

Successful reproduction of unmated Tropilaelaps mercedesae and its implication on mite population growth in Apis mellifera colonies

Successful reproduction by unmated Tropilaelaps mercedesae is reported here for the first time. Of the eight mature daughters that did not have male mates within their natal cells, four produced both mature sons and daughters, and four produced mature daughters only. Overall, 78% of the new daughters that had no egg-laying experience, and 84% of the foundresses that had or had not laid previously reproduced. Both inoculum daughter and foundress mites were collected from tan-bodied pupae and inoculated immediately. Therefore, our results suggest that phoresy is not required for reproduction in tropilaelaps mites. The ability of virgin females to lay both males and females (deuterotoky), and to reproduce without spending a phoretic period on adult bees may play major roles in tropilaelaps mites' competitive advantage over varroa mites in Apis mellifera colonies.

Chemosensing of honeybee parasite, Varroa destructor: Transcriptomic analysis

Chemosensing is a primary sense in nature, however little is known about its mechanism in Chelicerata. As a model organism we used the mite Varroa destructor, a key parasite of honeybees. Here we describe a transcriptomic analysis of two physiological stages for the Varroa foreleg, the site of primary olfactory organ. The transcriptomic analysis revealed transcripts of chemosensory related genes belonging to several groups. These include Niemann-Pick disease protein, type C2 (NPC2), gustatory receptors (GRs), ionotropic receptors (IRs), sensory neuron membrane proteins (SNMPs) and odorant binding proteins (OBP). However, no insect odorant receptors (ORs) and odorant co-receptors (ORcos) were found. In addition, we identified a homolog of the most ancient IR co-receptor, IR25a, in Varroa as well as in other members of Acari. High expression of this transcript in the mite's forelegs, while not detectable in the other pairs of legs, suggests a function for this IR25a-like in Varroa chemosensing.

Gene Duplication and Protein Evolution in Tick-Host Interactions

Ticks modulate their hosts' defense responses by secreting a biopharmacopiea of hundreds to thousands of proteins and bioactive chemicals into the feeding site (tick-host interface). These molecules and their functions evolved over millions of years as ticks adapted to blood-feeding, tick lineages diverged, and host-shifts occurred. The evolution of new proteins with new functions is mainly dependent on gene duplication events. Central questions around this are the rates of gene duplication, when they occurred and how new functions evolve after gene duplication. The current review investigates these questions in the light of tick biology and considers the possibilities of ancient genome duplication, lineage specific expansion events, and the role that positive selection played in the evolution of tick protein function. It contrasts current views in tick biology regarding adaptive evolution with the more general view that neutral evolution may account for the majority of biological innovations observed in ticks.

Characterization of the Copy Number and Variants of Deformed Wing Virus (DWV) in the Pairs of Honey Bee Pupa and Infesting Varroa destructor or Tropilaelaps mercedesae

Recent honey bee colony losses, particularly during the winter, have been shown to be associated with the presence of both ectoparasitic mites and Deformed Wing Virus (DWV). Whilst the role of Varroa destructor mites as a viral vector is well established, the role of Tropilaelaps mercedesae mites in viral transmission has not been fully investigated. In this study, we tested the effects that V. destructor and T. mercedesae infestation have on fluctuation of the DWV copy number and alteration of the virus variants in honey bees by characterizing individual pupae and their infesting mites. We observed that both mite species were associated with increased viral copy number in honey bee pupae. We found a positive correlation between DWV copy number in pupae and copy number in infesting mites, and the same DWV type A variant was present in either low or high copy number in both honey bee pupae and infesting V. destructor. These data also suggest that variant diversity is similar between honey bee pupae and the mites that infest them. These results support a previously proposed hypothesis that DWV suppresses the honey bee immune system when virus copy number reaches a specific threshold, promoting greater replication.

BlobTools: Interrogation of genome assemblies

The goal of many genome sequencing projects is to provide a complete representation of a target genome (or genomes) as underpinning data for further analyses. However, it can be problematic to identify which sequences in an assembly truly derive from the target genome(s) and which are derived from associated microbiome or contaminant organisms.  We present BlobTools, a modular command-line solution for visualisation, quality control and taxonomic partitioning of genome datasets. Using guanine+cytosine content of sequences, read coverage in sequencing libraries and taxonomy of sequence similarity matches, BlobTools can assist in primary partitioning of data, leading to improved assemblies, and screening of final assemblies for potential contaminants.  Through simulated paired-end read dataset,s containing a mixture of metazoan and bacterial taxa, we illustrate the main BlobTools workflow and suggest useful parameters for taxonomic partitioning of low-complexity metagenome assemblies.

Draft genome of the honey bee ectoparasitic mite, Tropilaelaps mercedesae, is shaped by the parasitic life history

The number of managed honey bee colonies has considerably decreased in many developed countries in recent years and ectoparasitic mites are considered as major threats to honey bee colonies and health. However, their general biology remains poorly understood. We sequenced the genome of Tropilaelaps mercedesae, the prevalent ectoparasitic mite infesting honey bees in Asia, and predicted 15 190 protein-coding genes that were well supported by the mite transcriptomes and proteomic data. Although amino acid substitutions have been accelerated within the conserved core genes of two mites, T. mercedesae and Metaseiulus occidentalis, T. mercedesae has undergone the least gene family expansion and contraction between the seven arthropods we tested. The number of sensory system genes has been dramatically reduced, but T. mercedesae contains all gene sets required to detoxify xenobiotics. T. mercedesae is closely associated with a symbiotic bacterium (Rickettsiella grylli-like) and Deformed Wing Virus, the most prevalent honey bee virus. T. mercedesae has a very specialized life history and habitat as the ectoparasitic mite strictly depends on the honey bee inside a stable colony. Thus, comparison of the genome and transcriptome sequences with those of a tick and free-living mites has revealed the specific features of the genome shaped by interaction with the honey bee and colony environment. Genome and transcriptome sequences of T. mercedesae, as well as Varroa destructor (another globally prevalent ectoparasitic mite of honey bee), not only provide insights into the mite biology, but may also help to develop measures to control the most serious pests of the honey bee.