key: cord-0799734-p95vzcoy authors: Schläppi, Daniel; Lattrell, Patrick; Yañez, Orlando; Chejanovsky, Nor; Neumann, Peter title: Foodborne Transmission of Deformed Wing Virus to Ants (Myrmica rubra) date: 2019-11-07 journal: Insects DOI: 10.3390/insects10110394 sha: ba29366173f97f54a22e5c410b3d05e9a9649d28 doc_id: 799734 cord_uid: p95vzcoy Virus host shifts occur frequently, but the whole range of host species and the actual transmission pathways are often poorly understood. Deformed wing virus (DWV), an RNA virus described from honeybees (Apis mellifera), has been shown to have a broad host range. Since ants are often scavenging on dead honeybees, foodborne transmission of these viruses may occur. However, the role of the ant Myrmica rubra as an alternative host is not known and foodborne transmission to ants has not been experimentally addressed yet. Here, we show with a 16-week feeding experiment that foodborne transmission enables DWV type-A and -B to infect M. rubra and that these ants may serve as a virus reservoir. However, the titers of both plus- and minus-sense viral RNA strands decreased over time. Since the ants were fed with highly virus-saturated honeybee pupae, this probably resulted in initial viral peaks, then approaching lower equilibrium titers in infected individuals later. Since DWV infections were also found in untreated field-collected M. rubra colonies, our results support the wide host range of DWV and further suggest foodborne transmission as a so far underestimated spread mechanism. Emerging infectious diseases (EIDs) can cause a significant impact on human-and animal health [1] [2] [3] and are often the product of a host shift, where a pathogen jumps from its original host to new species [4] . Prominent examples include severe acute respiratory syndrome (SARS) coronavirus or the still ongoing human immunodeficiency virus (HIV-1) pandemic [3] [4] [5] . Pathogens and their novel hosts lack a co-evolutionary history, which can lead to drastic effects on host populations [6] . In particular, RNA viruses have a high potential to cross species barriers because of high mutation rates that enable fast adaptive changes [7, 8] , which make them prominent among EIDs [9] . Such RNA viruses have also been suggested to contribute to the recent decline of wild pollinators and losses of managed ones, thereby potentially endangering valuable pollination services [10] [11] [12] [13] . Thus, it is not surprising that virus transmission between managed honeybees, Apis mellifera, and wild bees has recently been speculated [14, 15] . The transmission potential of honeybee viruses between bees and other ecosystem service providing insects has received considerably less attention [16] . It appears as if cross-species transmission Fifteen M. rubra colonies, collected in their founding stage around Berlin, Germany, were purchased from https://www.antstore.net/shop/de [50] . These colonies were kept in nesting tubes (155 mm length, 14 inner diameters) and fed with Drosophila hydei and sugar-water ad libitum (40% mass fraction of sugar). For the experiment, the colonies were transferred into boxes (60 × 70 × 15 mm) attached to a foraging arena (200 × 150 × 120 mm). Humidity in the nest was maintained with a moist cotton ball. Red foil covered the nest boxes to reduce stress [51] . The foraging arenas had an opening covered with removable insect nets and the edges were covered with paraffin oil to minimize escape. The colonies were maintained at RT (19-23 • C) with a natural day/night cycle. White-eyed honeybee pupae were obtained from sealed worker brood frames of two local A. mellifera colonies and then microinjected laterally between the second and third segment of the abdomen with 2 µL of a solution containing DWV [52] . The pupae were then incubated at 34.5 • C, ≥50% relative humidity and darkness for seven days [53] and afterward frozen at −80 • C. The DWV solution was prepared by homogenizing a DWV symptomatic honeybee in PBS and chloroform, followed by centrifugation (15800 rcf for 10 min). The supernatant was collected, diluted 1:1000, and frozen at −80 • C before injection into the pupae [52] . For the injected honeybee pupae used as food and positive controls for the virus detection, the logarithmic mean number of viral copies was 11.49 ± 0.29 and DWV type-specific primers showed a ratio of 5.88/1 for DWV-B to DWV-A, but no DWV-C. After four weeks, the colonies were randomly assigned to three treatments, differing in DWV infected protein supply over eight-weeks: (1) The controls (C, N = 3) were fed with D. hydei ad libitum; (2) treatment 1 (T1) colonies (N = 6) were fed for four weeks with DWV-infected honeybee pupae and then for eight weeks with D. hydei; (3) treatment 2 (T2) colonies (N = 6) were fed for eight weeks with DWV infected honeybee pupae and for 4 weeks with D. hydei. D. hydei were chosen as food for the controls, because honeybee pupae taken from local Swiss colonies inevitably contain background levels of DWV [54] . Sugar water (40% mass fraction of sugar) supply ad libitum was maintained. Then, ten adult workers were randomly sampled from each colony over four consecutive weeks until the end of the experiment (week [13] [14] [15] [16] and stored at −80 • C until further processing. The different feeding regimes were chosen to see if there is an impact of feeding duration on virus uptake and the time gap since the last feeding of the virus until the first sampling should reduce the risk of detecting DWV in the ant's gut content [55, 56] . Some colonies were not large enough to provide adults, thereby decreasing sample size at the last sampling (C -N = 1; T1 -N = 5; T2 -N = 6). For the quantification of viral titers a pooled sample of five workers per sampling week was crushed in 100 µL TN buffer (100 mM Tris, 100 mM NaCl, pH 7.6) with a 5-mm metal bead in 2-mL Eppendorf ® (Basel, Switzerland) tubes using a Retsch ® (Haan, Germany) MM 300 mixer mill for 1 min at the frequency 25 1/s [57] . Injected honeybee pupae extracted as well as positive controls using 200 µL TN buffer. Fifty microliters of the homogenate was transferred to 1.5-mL Eppendorf ® tubes and the RNA extraction was done using a NucleoSpin ® RNA II kit (Macherey-Nagel, Oensingen, Switzerland) following the manufacturer's recommendations. In the final step, RNA was eluted in 30 µL of elution buffer and stored at −80 • C [57] . To monitor the efficiency of the RNA purification and cDNA synthesis, tobacco mosaic virus (TMV, 1 µL TMV solution (1:1000)) was used as an exogenous internal reference added to each sample at the first step of the RNA extraction [58] . For the reverse transcription, we used the M-MLV RT (Promega, Dübendorf, Switzerland) Kit and followed the manufacturer's recommendations: Template RNA (50-500 µg, Table S1 ) was incubated with 0.75 µL of a random hexamer oligonucleotide (100 µM) and H 2 O for 5 min at 70 • C in a thermocycler (Biometra; Analytik Jena, Jena, Germany). For the cDNA synthesis, 5 µL 5× buffer, 1.125 µL nucleoside triphosphate (dNTP) 10 mM and 1 µL reverse transcriptase (M-MLV) was added to a final reaction volume of 25 µL, which was incubated at 37 • C for 60 min. The resulting cDNA was diluted 1/5 and stored at −25 • C until further processing. The RT-qPCR was performed in duplicate for each sample using the KAPA SYBR ® FAST Universal qPCR kit (KAPA Biosystems, Wilmington, DE, USA) with 12 µL volumes containing 6 µL KAPA SYBR ® (Wilmington, DE, USA) green reaction mix, 0.24 µL each of the forward and reverse primers for DWV or TMV (Table 1) , 2.52 µL H 2 O and 3 µL diluted cDNA [52] . An ECO™ real-time PCR machine (Illumina, San Diego, CA, USA) processed the reaction at the following qPCR cycling profile: 3 min incubation at 95 • C and 40 cycles of 3 s at 95 • C for denaturation, 30 s at 57 • C for annealing and extension, and data collection. By reading the fluorescence at 0.5 • C intervals between 55 • C and 95 • C, a melting curve analysis was performed after the amplification to verify the specificity of the PCR products. On each plate, a ten-folds serial dilution of purified PCR products that served as standard curves for DWV and TMV and two no-template negatives were included [59] . From the q-PCR output data, the standard curves and the experimental dilution factors viral titers respectively estimated viral copies per sample were derived [60] . To account for the exponential distribution of the data, virus titers were log-transformed and throughout the manuscript, the logarithmic values of the viral titers are reported. Because it is not possible to log transform zero values of negative samples we assigned to them a hypothetical Cq value of 36, which was transformed to virus titers as above. From the averaged value of these "negative titers," we obtained the titer detection threshold [61] . Table 1 . Primers used during qPCR for the detection and quantification of deformed wing virus (DWV) and minus strand-specific PCR as a proxy for viral replication as well as the DWV type-specific (A, B, and C) primers. The detection of the negative-sense strand RNA was used as a token of DWV replication [48, 49] . First, the negative-sense strand RNA was tagged and then amplified: The cDNA of all T1 and T2 colonies sampled at Week 13 and 16 were obtained by using Superscript ® III reverse transcriptase (Invitrogen, Carlsbad, CA, USA) following the manufacturer's recommendations (1 µL of DWV 3F tagged primer, Table 1 , 1 µL dNTP 10 mM, 4 µL of 5× first strand buffer, 2 µL 0.1 M DTT, 1 µL M-MLV reverse transcriptase, H 2 O and sample in a final reaction volume of 20 µL) [52] . The reaction was processed in a thermocycler (Biometra; Analytik Jena, Jena, Germany) according to the following protocol: 5 min at 65 • C, 10 min at 25 • C, 60 min at 50 • C followed by 15 min at 70 • C. The obtained cDNA then was purified using the NucleoSpin ® Gel & PCR Clean-up kit (Macherey-Nagel) and eluted in 30 µL elution buffer to avoid false-positive results. A five-fold dilution of the purified cDNA was used for the subsequent amplification with conventional PCR (MyTay™ kit by Bioline, London, UK). DWV-R1 was used as the reverse primer, while a Tag oligonucleotide (Table 1 ) was used as forward primer. For each sample, a second reaction was carried out without forward primer to check for efficient removal of DWV 3F tagged primer after the reverse transcription. The thermal cycling was done with a profile that consisted of 2 min incubation at 95 • C and 35 cycles of 20 s at 95 • C for denaturation, 20 s at 42 • C for annealing, and 30 s at 72 • C for extension. By electrophoresis on a 1.5% agarose gel, the PCR products were analyzed and visually checked under UV light. If a clear band was present at 221 bp and the associated negative control was indeed negative, a sample was considered positive. All statistical analyses were performed using R version 3.5.1 [64] and R Studio version 1.1.456 [65] . Data were checked for normality with the Shapiro-Wilk test and homogeneity of variances with the Levene's test and subsequent statistical tests chosen accordingly. We used the lmer function from the package lme4 [66] to perform a linear mixed effect analysis to test for differences of viral titers over time between the treatments. The logarithm of the viral titers entered the model as a response variable. As fixed effects, we used sampling week and treatment including their interaction term. As a random effect, we added the colony identity to control for multiple testing of the same colonies (week [13] [14] [15] [16] . Visual inspection of residual plots did not reveal deviations from homoscedasticity or normality. p-values were obtained by comparing the model with the effect in question against the model without it, using an ANOVA as a likelihood ratio test. The post-hoc pairwise comparison was then performed using the Bonferroni correction. All the samples from T1 (colonies fed for four weeks with DWV-infected honeybee pupae) and T2 (colonies fed for eight weeks with DWV-infected honeybee pupae) were strongly positive, with two exceptions (one pooled sample each during week 13 and 15). With 4.17/1 (T1) and 5.46/1 (T2), the ratio of DWV type-B to type-A was similar to the one in the feeding pupae (5.88/1 for DWV-B to DWV-A). Sequencing of the PCR products confirmed the identity of both DWV strains found in M. rubra (GenBank accession: DWV-A, JF346565.1, 99.73% identity, 100% query cover; DWV-B, AY251269.2, 98.84% identity, 99% query cover). The viral titers found in ants were in the range of 6.46-9.66 with a log mean of 7.95 ± 0.72 for T1, 6.17-10.18, 8.26 ± 1.02 for T2, and 3.17-5.88, 4.22 ± 0.85 in the control group seen over all weeks ( Figure 1 , Table 2 ). Viral titers in ants were significantly affected by the factors treatment (X 2 (2) = 42.39, p < 0.001) and sampling week (X 2 (1) = 38.43, p < 0.001), including their interaction term (X 2 (2) = 8.12, p < 0.017). Viral titers of both treatment groups were significantly higher than in the control group (p < 0.0001, Bonferroni post hoc test), but there was no significant difference between the two treatment groups (p = 0.53, Bonferroni post hoc test). The significant interaction between treatment and sampling week indicates that the slope of the viral titer decrease was dependent on treatment. 9.47 ± 0.51 [6] 8.47 ± 0.48 [6] 8.01 ± 0.65 [6] 7.12 ± 0.65 [6] 3. The presence of the negative-sense strand of DWV, as indicated by a band at 221 bp of the agarose gel electrophoresis [49, 52] was confirmed for all colonies of T1 at week 13, but only for two colonies at week 16 ( Figure 2) . Similarly, the band at 221 bp was clearly visible in all but one colony of T2 at the first sampling point. Four weeks later, the band was still well visible in one colony, less pronounced in two others, and not visible in three colonies. Our results suggest that the ant M. rubra can be an alternative host for both DWV A and B as a result of foodborne transmission because of the presence of both high viral titers and the negative-sense RNA strand for up to 13 weeks after the last virus-infected food has been consumed. These results further support that DWV has a wide host range including other predatory insects [14, 15, [17] [18] [19] [20] [21] [29] [30] [31] [32] 37] . Interestingly, there were only two ant species from the family Formicidae among the more than 60 species in which DWV has been detected so far [21] . It is obvious that recently consumed highly virus infected food like honeybee pupae, containing the virus negative-sense strand, might lead to false positive results simply because of the presence of viral particles in ant's gut. Because of the delayed sampling, we can however expect that the gut content of the ants has been cleared out at least once before the virus analyses [55] . The experimental data show that foodborne transmission, e.g., via consumption of infected honeybees, can be an underlying mechanism, which enables DWV to infect M. rubra and potentially other arthropod predators and scavengers [29] [30] [31] [32] 37] . Indeed, both treatment groups that were fed with virus-infected honeybee pupae had consistently higher viral titers than the controls. Since symptomatic honeybees with an overt infection have 10 10+ genomic equivalents of DWV and at lower infection levels 10 3 -10 9 genomic copies [67] , the titers observed in M. rubra seem to reflect covert infections. However, the average honeybee worker (~120 mg) weights about 60 times more than a M. rubra worker [68, 69] . When controlled for these bodyweight differences, the DWV A and B titers found in the treatment groups therefore almost reach infection levels per mg ant-tissue analogous to an overt infection in bees. Further, the viral titers were very low in the control colonies, thereby indicating that DWV A and B were only replicating in the treatments. Interestingly, the ratios of DWV type-B to type-A were very similar in the ants, regardless of the treatment (4.17/1 (T1), 5.46/1 (T2), as well as in the injected honeybee pupae (5.88/1 for DWV-B to DWV-A), thereby suggesting that both strains had similar reproductive capacities in both hosts. Over time, a decrease of both the viral titers and the presence of the negative-sense strand was observed in the treatment colonies. It may well be that the treatment colonies would eventually reach a low viral level such as the controls, display low level replication and at that both A and B strains are present. However, we have no data to address those questions. What could explain the observed decrease of viral titers over time? Neither altruistic self-removal of infected nestmates [70] nor increased mortality in the treatment colonies were observed. Therefore, the higher initial viral titers and the more common minus-strand presence could reflect virus particles consumed with the infected honeybee food. In light of ant physiology [71] , this is nevertheless unlikely to explain the comparatively high viral titers and especially the presence of the minus strand for up to 13 weeks after the last virus-infected food has been consumed. The most parsimonious explanation therefore is that the ants were fed with highly virus saturated honeybee pupae resulting in initial viral peaks, then approaching lower equilibrium titers in infected ant individuals later. In light of these results, it appears as if controlled laboratory experiments are required before drawing general conclusions on the role of predators/scavengers as alternative hosts of any honeybee viruses, i.e., the timing of the last virus-infected food consumption should be known. Similarly, the detection of the negative-sense strand in V. destructor suggested DWV replication in the mites until recently, when a non-propagative manner DWV transmission has been proposed [72] . Given the prevalence and high titers of DWV and its wide host range [21] , it is not surprising that the field-collected control colonies were also tested positive for DWV at low levels. Since our experimental data clearly show that DWV transmission to ants is feasible via food, we assume that either (1) DWV-contaminated food must have been given to the ant colonies in captivity prior to our experiments or alternatively, but not mutually exclusive, (2) the D. hydei flies used as control food were DWV positive, which we regard as unlikely, or (3) the control colonies had consumed DWV-infected food prior to field-sampling by the company. Given the latter case, this would support the view that DWV transmission to M. rubra also occurs in the field, which would be in line with other studies showing the presence of viruses, originally described from honeybees, in ants [29, 37] . A closer look at the controls also reveals that the viruses are not equally distributed within these naturally infected colonies. It would be interesting to see if the uneven distribution is due to age-related division of labor and non-stochastic interaction patterns [73] , i.e., foragers might have a higher exposure risk compared to in-house workers. Considering the low viral titers in the controls, it is possible that the viruses persist in M. rubra as a covert infection that might turn into an overt infection under suitable conditions, e.g., parasite infections as in honeybees [26] . Although the data suggest that DWV prevails and replicates in M. rubra and that it also occurs in field-collected samples, it provides no insights on pathogenicity. Neither abnormalities in colony development nor individual clinical symptoms in workers were observed. Additional studies are therefore required to understand if and how honeybee viruses might affect M. rubra or other predators/scavengers. The data suggest for the first time that M. rubra is another biological host for DWV type-A and -B known from honeybees and that foodborne transmission is an underlying transmission mechanism. Since ants are ubiquitous and important in terrestrial ecosystems, it appears as if there is a considerable potential for virus transmission between wild insects and managed bees with possible consequences for ecosystem functioning. It has been suggested that viruses might offer novel options for the control of invasive ant species [37] . However, we consider such biocontrol suboptimal due to the evident wide host range of viruses and the here demonstrated considerable chances of host-switches, with possibly dramatic consequences [5] . Supplementary Materials: The following are available online at http://www.mdpi.com/2075-4450/10/11/394/s1, Table S1 : DWV titers and positive/negative assignment of colonies will be available online. Emerging infectious diseases of wildlife-Threats to biodiversity and human health The challenge of emerging and re-emerging infectious diseases Global trends in emerging infectious diseases The evolution and genetics of virus host shifts AIDS as a zoonosis: Scientific and public health implications Biological and biomedical implications of the co-evolution of pathogens and their hosts Origins of major human infectious diseases González-Candelas, F. The population genetics and evolutionary epidemiology of RNA viruses Emerging pathogens: The epidemiology and evolution of species jumps Global pollinator declines: Trends, impacts and drivers Honey bee colony losses Emerging viral disease risk to pollinating insects: Ecological, evolutionary and anthropogenic factors Worldwide decline of the entomofauna: A review of its drivers Disease associations between honeybees and bumblebees as a threat to wild pollinators A sting in the spit: Widespread cross-infection of multiple RNA viruses across wild and managed bees The economic value of ecological services provided by insects Detection of deformed wing virus, a honey bee viral pathogen, in bumble bees (Bombus terrestris and Bombus pascuorum) with wing deformities New evidence that deformed wing virus and black queen cell virus are multi-host pathogens Cross-species transmission of honey bee viruses in associated arthropods Widespread occurrence of honey bee pathogens in solitary bees Deformed wing virus in Honeybees and other Insects Deformed wing virus Predictive markers of honey bee colony collapse Clinical signs of deformed wing virus infection are predictive markers for honey bee colony losses Global honey bee viral landscape altered by a parasitic mite Varroa invasion and virus adaptation Deformed wing virus is a recent global epidemic in honeybees driven by Varroa mites Knock-on community impacts of a novel vector: Spillover of emerging DWV-B from Varroa-infested honeybees to wild bumblebees Detection of Chronic bee paralysis virus (CBPV) genome and its replicative RNA form in various hosts and possible ways of spread A scientific note on Israeli acute paralysis virus infection of Eastern honeybee Apis cerana and vespine predator Vespa velutina Detection of deformed wing virus in Vespa crabro First detection of replicative deformed wing virus (DWV) in Vespa velutina nigrithorax Socialized medicine: Individual and communal disease barriers in honey bees Field assays for hygienic behavior in honey bees (Hymenoptera: Apidae) Hygienic behaviour of honey bees and its application for control of brood diseases and Varroa: Part I. Hygienic behaviour and resistance to American foulbrood Evidence for antiseptic behaviour towards sick adult bees in honey bee colonies Invasive ants carry novel viruses in their new range and form reservoirs for a honeybee pathogen A metatranscriptomic survey of the invasive yellow crazy ant, Anoplolepis gracilipes, identifies several potential viral and bacterial pathogens and mutualists Single-stranded RNA viruses infecting the invasive Argentine ant Soil invertebrates and ecosystem services Ecosystem engineering and predation: The multi-trophic impact of two ant species The little things that run the world revisited: A review of ant-mediated ecosystem services and disservices (Hymenoptera: Formicidae) Ant biodiversity and its relationship to ecosystem functioning: A review Now you see them, now you don't!-population crashes of established introduced species The evolution of epidemic influenza Die Ameisen Mittel-und Nordeuropas Distribution of an invasive ant, Myrmica rubra (Hymenoptera: Formicidae), in Maine Development of a strand-specific RT-PCR based assay to detect the replicative form of hepatitis C virus RNA RT-PCR analysis of Deformed wing virus in honeybees (Apis mellifera) and mites (Varroa destructor) The influence of red light on the aggregation of two castes of the ant, Lasius niger Standard methods for virus research in Apis mellifera Standard methods for maintaining adult Apis mellifera in cages under in vitro laboratory conditions Protein nutrition governs within-host race of honey bee pathogens The multiple recruitment systems of the African weaver ant Oecophylla longinoda (Latreille) (Hymenoptera: Formicidae) Nest etiquette-where ants go when nature calls Standard methods for molecular research in Apis mellifera Comparative analysis of deformed wing virus (DWV) RNA in Apis mellifera and Varroa destructor The MIQE guidelines: Minimum information for publication of quantitative real-time PCR experiments Deformed wing virus and drone mating flights in the honey bee (Apis mellifera): Implications for sexual transmission of a major honey bee virus On the front line: Quantitative virus dynamics in honeybee (Apis mellifera L.) colonies along a new expansion front of the parasite Varroa destructor Method development and application to quantify the role of three DWV master variants in overwinter colony losses of European honey bees Viral load estimation in asymptomatic honey bee colonies using the quantitative RT-PCR technique R: A Language and Environment for Statistical Computing, R version 3 R Foundation for Statistical Computing RStudio: Integrated Development Environment for R Fitting linear mixed-effects models using lme4 Replication of Varroa destructor virus 1 (VDV-1) and a Varroa destructor virus 1-deformed wing virus recombinant (VDV-1-DWV) in the head of the honey bee Hydrolases of developing worker brood and newly emerged worker of Apis mellifera carnica Feeding and growth in the ant Myrmica A review on self-destructive defense behaviors in social insects The Ants Deformed wing virus type A, a major honey bee pathogen, is vectored by the mite Varroa destructor in a non-propagative manner We wish to thank Kaspar Roth for technical support and Laura Bosco for fruitful discussions. The authors declare no conflict of interest.