Article

Up-regulated pathways in response to Deformed Wing Virus infection in Apis mellifera (Hymenoptera: Apidae)

Vías de señalización activadas en respuesta a la infección por el virus de las alas deformes en Apis mellifera (Hymenoptera: Apidae)

 

QUINTANA, Silvina^1,2,3, BRASESCO, Constanza^2,3*, NEGRI, Pedro^2,3, MARIN, Maia³, PAGNUCO, Inti³, SZAWARSKI, Nicolás^2,3, REYNALDI, Francisco^3,4, LARSEN, Alejandra⁵, EGUARAS, Martín^2,3 & MAGGI, Matías^2,3

¹ Laboratorio de Biología Molecular, Instituto de Análisis Fares Taie. Mar del Plata, Argentina.
² Centro de Investigación en Abejas Sociales, Laboratorio de Artrópodos, Facultad de Ciencias Exactas y Naturales, Universidad Nacional de Mar del Plata. Mar del Plata, Argentina. * E-mail: cobrasesco@gmail.com
³ Consejo Nacional de Investigaciones Científicas y Técnicas. Ciudad Autónoma de Buenos Aires, Argentina.
⁴ Laboratorio de Virología (LAVIR), Facultad de Ciencias Veterinarias, Universidad Nacional de La Plata. La Plata, Argentina.
⁵ Inmunología Aplicada, Facultad de Ciencias Veterinarias, Universidad Nacional de La Plata. La Plata, Argentina.

Received 08 - VI - 2018 | Accepted 21 - II - 2019 | Published 28 - III - 2019 https://doi.org/10.25085/rsea.780101

 

RESUMEN. La transmisión de virus por Varroa destructorAnderson & Trueman a las abejas ha sido identificada como una de las principales amenazas para las colonias de Apis melliferaL. También se ha sugerido que la interacción sinérgica entre parásitos y patógenos mediada por la inmunidad de A. mellifera puede conducir al colapso de las colonias de abejas melíferas. Analizamos la relación entre la carga viral del virus de las alas deformes (DWV) y los niveles de expresión génica de Toll Wheeler, cactus, domeless, TEPA, óxido nítrico sintasa, himenoptecina, vitellogenina, defensina, abaecina en muestras de abeja invernal de A. mellifera con y sin alas deformadas que se obtuvieron de colmenas en el campo con alta infestación de V. destructor. En las abejas con alas deformadas encontramos aumento en las cargas virales de DWV y regulación positiva de los genes de la vía de señalización Jak-STAT.

PALABRAS CLAVE. Abeja melífera. Virus. Vías del sistema inmune de insectos. Expresión de inmunogenes. Cargas virales.

ABSTRACT. Transmission of viruses by Varroa destructorAnderson & Trueman to bees has been identifed as one of the main threats to Apis mellifera L. colonies. It has also been suggested that synergistic parasite-pathogen interaction mediated by A. mellifera immunity can drive the collapse of honey bee colonies. We analyzed the relationship between deformed wing virus (DWV) viral load and gene expression levels of Toll Wheeler, cactus, domeless, TEPA, nitric oxide synthase, hymenoptaecin, vitellogenin, defensin, abaecin in winter bee samples of A. melliferawith and without deformed wings. Samples were collected from hives in the feld with high V. destructor infestation. In bees with deformed wings we found an increase in DWV titre and upregulation of genes in the Jak-STAT signalling pathway.

KEYWORDS. Honey bees. Virus. Insect immune system pathways. Immunogene expression. Viral load.

 

Alarming increases in depopulation and the loss of bee colonies have been reported over the last few years globally (Le Conte et al., 2010; Potts et al., 2010). These losses are of great ecological, agricultural and economic relevance as honey bees pollinate natural and agricultural ecosystems worldwide. This critical situation has major impacts on the fragile balance between bees and plants, forcing many scientists to redirect their research areas to tackle the possible factors which cause this problem. Most researchers agree that there is no single explanation for the large losses of colonies observed, but they are the result of a synergistic interaction between different stressors (Hendriksma et al., 2011). Argentina is no exception, with beekeepers across the country reporting cases of depopulation of hives and colony losses. Unfortunately, there is no nationwide data for this available (Maggi et al., 2013).

Transmission of viruses by the mite Varroa destructor Anderson & Trueman is considered one of the main threats to Apis mellifera colonies (Martin et al., 2012). Viruses replicate within mites (without causing the death of the parasite) and consequently infect the bees with high viral loads, as is the case for the deformed wing virus (DWV) (Gisder et al., 2009; Martin et al., 2012). So far, 24 different viruses have been isolated from the honey bee and it has been demonstrated that some can be vectored by V. destructor (Genersch & Aubert, 2010; Runckel et al., 2011;de Miranda et al., 2013; Dalmon et al., 2017). These include Kashmir bee virus, Sacbrood virus (SBV), acute bee paralysis virus (ABPV), Israeli acute paralysis virus (IAPV) and DWV (Genersch & Aubert, 2010). ABPV and DWV have been implicated in colony losses in Germany and in the colony collapse disorder (CCD) observed in the USA (Cox-Foster et al., 2007; vanEngelsdorp et al., 2007, 2008, 2009; Genersch, 2010).

In Argentina, the information about bee viruses is scarce and fragmented. However, Reynaldi et al. (2010) reported the molecular identification of chronic bee paralysis virus (CBPV), ABPV and SBV in Buenos Aires, the main honey-producing province of Argentina. In the study, a lower rate of infection and some cases of co-infection with more than one virus were found when compared to other countries of South America (Antúnez et al., 2006). In 2011, the same authors reported the presence of IAPV by RT-PCR in samples taken from several provinces. These data indicated a high frequency of IAPV in asymptomatic hives in Argentina (Reynaldi et al., 2011).

During 2013 our research group collected and

updated the distribution of SBV, DWV, IAPV and ABPV in colonies of the South-eastern region of Buenos Aires province. Viruses were detected in 33.3% of all samples tested. SBV was identifed in pupae, IAPV and DWV was identifed in adult bees and in V. destructormites. ABPV was not detected in any of the samples analyzed (Brasesco et al., 2013). High loads of DWV in bees with deformed wings and lower viral loads in bees with normal wing appearance were observed, as previously reported (Yang & Cox-Foster, 2005; Martin et al., 2012). In this manner, the information was expanded to provide new molecular records about the presence of such bee viruses in the country (Brasesco et al., 2013).

Martin and colleagues studied viral load, prevalence and genetic diversity of DWV in four Hawaiian Islands, finding differences in viral load and also variations in DWV strain diversity (Martin et al., 2012). Ryabov et al. (2014) also described the existence of variants of DWV associated with differing virulence.

Some reports have described changes in the expression levels of genes of the immune system of A. mellifera in response to pathogen infection, like ABPV and Paenibacillus larvae(Siede et al., 2012; Boncristiani et al., 2013). It has also been suggested that synergistic parasite-pathogen interaction mediated by A. mellifera immunity can drive the collapse of honey bee colonies (Nazzi et al., 2012). There is considerable interest in understanding the mechanisms mediating non-self-recognition, activation of signal transduction pathways that controls the expression of immune genes, and the identity and mode of action of the biocidal molecules used by insects (Nappi et al., 2000).

Vitellogenin (Vg) is a yolk precursor as well as a pathogen pattern recognition receptor (Zhang et al., 2011). It is a nutritious lipoprotein synthesized by the invertebrate fat body or vertebrate liver, secreted to the hemolymph/blood and taken up by nurse cells and eggs through receptor-mediated endocytosis (Finn, 2007). Vg concentration varies between the members of the honey bee colony, from almost undetectable to 40% of the total hemolymph protein fraction in the functionally sterile helper females -called workers- and it constitutes up to 70% of the hemolymph protein fraction in the egg-laying queens (Barchuk et al., 2002). Vg has an important role in bee's immune response being the carrier of immune elicitors; it binds to bacteria, both the gram-positive P. larvae, causing American foulbrood disease, and the gram-negative Escherichia coli(Salmela et al., 2015). It also binds to pathogen-associated molecular patterns; lipopolysaccharide, peptidoglycan and zymosan, as demonstrated by Salmela et al. (2015). Furthermore, there have also been previously described changes in Vg expression in response to microsporidia, for example a decrease after Nosema ceranae infection (Antúnez et al., 2009). Therefore, it is interesting to study the implications of this molecule during virus infections in bees.

The major mechanism of antiviral defence in insects is the small interfering RNA pathway that responds through the detection of virus-derived double-stranded RNA to suppress virus replication. However, other innate antimicrobial pathways such as Imd, Toll, and Jak-STAT and the autophagy pathway have also been shown to play important roles in antiviral immunity (Evans et al., 2006; Kingsolver et al., 2013). In insects, Toll receptors are necessary for sensing the invasion by pathogens (Takeda & Akira, 2003). Recent reports have implicated Drosophila Toll and Jak-STAT signalling pathways in controlling RNA viral infection (Merkling & van Rij, 2013). Initially characterized for its role in development and hemocyte proliferation, it has been shown that the Jak-STAT pathway also responds to bacterial and viral infections by regulating the production of downstream effector molecules including antimicrobial peptides (AMPs). This pathway is composed of one signal transducer and activator of transcription (STAT92E; hereafter STAT), a Janus kinase (Hopscotch), a receptor (Domeless), and three unpaired-related ligands (Kingsolver et al., 2013).

Also, data from various research groups, using both RNA and DNA viruses in experimental systems, have shown the role for nitric oxide (NO) in insect immune response. The data suggests that NO may inhibit an early stage of viral replication and thus prevent viral spread, promoting viral clearance and recovery of the host (Reiss & Komatsu, 1998). Further, it has been found that in insects, nitric oxide synthase (NOS) expression is also regulated by the STAT pathway (Chen et al., 2002).

A deeper knowledge of the response of selected immune genes to viral infection may help to better understand the complex interactions between host and pathogen, illuminating the phenomenon of colony bee losses. In order to examine this interaction between immune genes and virus, bees with deformed wings that normally carry high viral load (Martin et al., 2012) and bees with normal wings that are expected to show lower viral levels were used to study viral load of the DWV and to analyse the changes in the expression levels of Toll Wheeler, cactus, domeless, TEPA, NOS, hymenoptaecin, Vg, defensin, and abaecin in response to the virus infection. This study was conducted in bee samples with and without deformed wings which were collected during winter from hives in the feld with high V. destructorinfestation.

Collection of honey bees

Ten adult winter bees with deformed wings (DW) and 10 adult winter bees with normal wings (NW) were collected from feld hives highly infested (5.45% of phoretic Varroaprevalence in the apiary), naturally, with V. destructor in Santa Paula apiary (37° 56' 00.69" S; 57° 40' 40.53" W). Overwintering losses associated with V. destructor have been detected in this region in the past (Eguaras & Ruffinengo, 2006; Invernizzi et al., 2011; Maggi et al., 2013); thus, it is expected to find hives highly infested in winter.

RNA extraction, DNAse treatment and reverse transcription

Total RNA from individual adult A. mellifera was isolated using TRI-reagent (Biobasic Inc., USA) according to the manufacturer's protocol. DNase I Amplification Grade (Invitrogen, Carlsbad, CA, USA) was used to degrade genomic DNA carried over in the RNA extraction by digestion at 37 °C for 30 min. The quantity of the resulting RNA was determined in a Rotor Gene 6000 cycler (Qiagen, USA) using a Quant-iTT RiboGreen® RNA Assay Kit (Invitrogen, USA). Complementary DNA (cDNA) was synthesized using a reaction mixture containing 1 µg of total RNA, random hexamers (12 ng/μl) and Moloney murine leukemia virus reverse transcriptase (Invitrogen, USA) according to the procedures suggested by the manufacturer. Negative controls, omitting either the RNA or the reverse transcriptase were included during reverse transcription.

Primers

A list of the primers for qPCR (synthesized by Eurofins MWG Operon, USA) used in PCR protocols, specific annealing temperature and melting temperature are summarised in Table I. Primers for the study of Toll wheeler mRNA levels were designed with Primer Premier Software (PREMIER Biosoft International, Palo Alto, USA).

RT-qPCR

To determine whether DWV infection has an effect o n A. mellifera immune related genes expression, the levels of Toll Wheeler, cactus, domeless, TEPA, NOS, hymenoptaecin, Vg, defensin and abaecin mRNA were quantifed. mRNA expression levels of all analyzed genes from DW bees (10) were compared to their expression in NW bees (10).

RT-qPCR reactions were run with KAPA HRM FAST qPCR kit, with EvaGreen as intercalating dye (KAPA Biosystems, USA) and 1 µl of cDNA sample in a final volume of 20 µl. The amplification and detection of the specific products were carried out using a Rotor Gene 6000 cycler (Qiagen, Germany), with the following amplification conditions: 2 min at 95 °C, 40 cycles of 10 s at 95 °C, 15 s at specific annealing temperature and 15 s at 72 °C. After amplification, a melting curve analysis was performed, which resulted in single product-specific melting curve. In all cases, the experiments were performed in duplicate. Negative controls for cDNA synthesis and qPCR procedures were included in all cases. The amplification efficiency was determined for each gene using 10-fold dilutions of cDNA.

The results are reported as the mean fold change of DWV and immune system genes transcription levels in DW bees over levels detected in NW bees, which served as the control group (Pfaff et al., 2002). Apis mellifera β-actin gene was used as a reference / housekeeping gene.

DWV analysis

For the determination of DWV levels and strain diversity, primers that amplify a ~250 bp fragment of the RNA-dependent RNA polymerase (RdRp) gene of DWV were used (Kukielka et al., 2008). PCR amplifications were performed in duplicate and a High-Resolution Melting (HRM) analysis was carried out in all 20 samples immediately after the DWV RT-qPCR at temperature increment steps of 0.1 °C from 75 °C to 95 °C. The results were analyzed with Rotor Gene Q software, version 1.7.94. The genotype identification confdence cut off value of 90% was used as the threshold in the genotyping module of the software and PCR reactions of each sample were performed in triplicate.

Statistical analysis

The relative expression analysis of DWV viral levels and immune related genes was performed using the Relative Expression Software Tool (REST, Qiagen Inc., Valencia, CA, USA), which compares the expression of the target gene in a sample relative to the control and tests the group differences for significance with a pair wise fxed reallocation randomization test (Pfaff et al., 2002). RT-qPCR efficiency for each gene was determined by a linear regression model, according to the equation: E = 10[-1/slope].

RESULTS

Viral load

Relative levels of DWV RNA were analyzed and DW bees showed a 178-fold higher viral load compared to NW bees (p<0.01).

Sequencing

After HRM analysis of DWV amplifications, the PCR amplicons were resolved on agarose gels and purifed using Accuprep Gel Purification Kit (Bioneer, South Korea). Purifed fragments were sequenced using BigDye Terminator v1.1 Cycle Sequencing Kit (Applied Biosystems) with the ABI Prism 310 Genetic Analyzer (Applied Biosystems). To verify the differences identifed by HRM analysis, all sequences obtained were aligned and compared by use of MAFFT software on-line (Katoh et al., 2017). Using ExPASy Translate tool online (Artimo et al., 2012) we obtained and analyzed the amino acid coding sequence of the fragments of the non-structural polyprotein of the DWV, of each of the nucleotide sequences obtained.

Table I. Primer sequences (forward-fw and reverse-rv), annealing temperature (AT), and melting temperature (MT) used in this study (see details on main text).

Gene expression studies

No statistically significant differences (p >0.05) were found for the expression levels of Toll Wheeler, Vg, defensin and abaecin mRNA in DW bees when compared to NW (Fig. 1). However, significant up regulation of cactus, domeless, TEPA, NOS and hymenoptecin mRNA levels (p<0.05) was found in DW bees when compared to NW (Fig. 1).

DWV analysis

HRM analysis of DWV RT-qPCR products showed that amplifications of DW bees were primarily dominated by a single genotype cluster, which showed the same HRM profle in all analyzed samples.

In the case of NW samples, half of them (5) exhibited different melting profles, indicating the presence of a range of DWV variant sequences. The differences of fve DW and fve NW bee samples are shown in a Normalized temperature-shifted difference graph (Fig. 2), in this set of samples all DW bees have the same melting profle, while two NW samples show a different melting profle. T o verify the differences identifed by HRM analysis, all sequences were aligned and compared by use of MAFFT software on-line (Katoh et al., 2017) (Fig. 2). In this way the strain diversity present in NW samples was corroborated at sequence level observing that sample A9 had two SNPs in the middle of the sequence, and sample A10 had one SNP. In the case of A9 a change in the amino acid coding sequence of the RdRp protein could be observed (Fig. 3).

Blast search (Altschul et al., 1990) of these sequences share over 98% of identity to those reported in countries such as United Kingdom, Brazil, and Syria (GenBank accession: DQ434932, DQ434942, KT733631, KX530470). Using MAFFT software on-line (Katoh et al., 2017) a phylogenetic tree was created to visualize the clusters of the sequences obtained with the top hits from the GenBank (Fig. 4).

DISCUSSION

DWV is considered the most widespread bee virus, often associated with V. destructor infestation, and the role of this parasite in the transmission of the virus has already been demonstrated (Allen & Ball, 1996). In the present work we studied DWV viral load, strain diversity and A. mellifera immunogene expression in samples of bees, with and without deformed wings, collected from feld hives highly infected by V. destructor.

DW bees showed a 178-fold higher viral load in comparison to NW bees (p<0.01). These results agree with those reported by Yang and Cox-Foster who also described that DW bees had higher virus titers than NW bees (Yang & Cox-Foster, 2005). It is worth mentioning that virus detection does not implicate an established infection (Lanzi et al., 2006) as this virus, as well as several others RNA bee viruses, usually persist as covert infections in the colony without the vectoring capacity of Varroa (Locke et al., 2012).

DWV detection was performed amplifying a ~250 bp fragment of the RNA-dependent RNA polymerase (RdRp) gene. RdRp is an essential protein encoded in the genomes of all RNA viruses with no DNA intermediate stage. RdRp catalyzes synthesis of the complementary RNA strand, is highly conserved (Lanzi et al., 2006) and vital for viral replication. Substitutions in its nucleotide sequence could contribute to variations in the amino acid coding sequence of this protein, hence in viral epidemiology, since any variation in this gene could suggest an alteration in the genetic shift of the rest of the viral genome (Ryabov et al., 2014). Here, we found that the portion of RdRp sequence analyzed from the local viruses share over 98% of identity to those reported in countries such as United Kingdom, Brazil, and Syria (Berényi et al., 2007; Wang et al., 2013).

High-resolution melting (HRM) analysis has been previously used to monitor the relationship between Varroainfestation and virus diversity (Martin et al., 2012). We have also used this technology as a virus identification method (Marin et al., 2016). In the present work we analyzed a region of the RdRp gene to exploit the possibility of differentiating virus strains by their slight differences in specific nucleotides, as observed in the alignment of the amplifed fragments of DWV RdRp gene fragment. The alignment revealed differences in the 250 bp amplifed product, with distinct melting peaks (Fig. 2A); all DWV amplifications of DW bees showed the same melting profle, revealing a unique viral strain in these samples, while fve NW samples showed different melting profles, showing strain diversity in bees without clinical symptoms. Sample A9 (NW bee) had two SNPs in the middle of the sequence, and sample A10 (NW bee) had one SNP. In the case of A9, an alteration of the amino acid coding sequence of the fragment of the RdRp protein analyzed was found (Fig. 3).

The DW honey bees analyzed here were sampled from feld colonies that were parasitized by V. destructor for the necessary period of time (more than three years of infestation) required for the selection of virus variants adapted to mite transmission (Martin et al., 2012).

Fig. 1. Relative mRNA expression levels of Toll Wheeler, cactus, abaecin, defensin, hymenoptaecin, domeless, NOS, TEPA and Vg in NW and DW samples

Asterisks indicate statistically significant differences (P < 0.05) in relation to NW samples (control). mRNA relative levels and statistical significance of the differences in mRNA expression were obtained with the software REST (Qiagen Inc, USA).

Fig. 2. A. HRM analysis of DWV-RT-qPCR products. Normalized temperature-shifted difference graph, DWV amplicons of 5 DW samples (blue) and 5 NW samples (red) analyzed in duplicate are shown. The melting profile of a DWV RT-qPCR product of a DW bee sample has been chosen as a horizontal baseline. B. Partial alignment of DWV- RT-qPCR products between 2 NW samples (A9 and A10) which showed different melting profiles and 2 DW samples (AD6 and AD7).

 

A large number of polymorphisms were described in a number of regions in the deep sequencing data of DWV by Wang et al. (2013) suggesting that sequence conservation might not necessarily be maintained at the population level. Further studies including complete viral genome sequencing and experimental infections of honey bees with different DWV strains are required for testing virulence.

Some studies have found no correlation between the presence of Varroa and changes in host immune responses (Gregory et al., 2005; Navajas et al., 2008). Yet, Marringa et al. (2014) demonstrated changes in the type of hemocytes freed in the hemolymph, suggesting a cellular response to Varroaattack and could be the link to induction of the cellular responses pathways. Also, Yang and Cox-Foster (2005) suggest that the increase in DWV titre and reduction in variant diversity cannot be solely explained by Varroa-induced immunosuppression of honey bees. Considering this background, here we studied the expression levels of genes of various effector pathways of the immune system of A. mellifera i n samples of DW and NW bees to make a contribution to the studies focussed in this host-pathogen interaction.

Antimicrobial peptides expression levels have been previously studied in A. mellifera in relation to viral infections. It has been reported that the expression of these peptides is induced by activation of Toll and Imd pathways (Brutscher et al., 2015). In this study, three antimicrobial peptides were analyzed, abaecin, defensin-1 and hymenoptaecin. Only hymenoptaecin mRNA expression levels were affected in DW bees, showing a 38-fold increase in comparison with NW bees. Hymenoptaecin expression was not correlated to either abaecin or defensin, this result could indicate that hymenoptaecin uses a different regulatory pathway than abaecin and defensin as previously suggested (Yang & Cox-Foster, 2005). It was proposed that hymenoptaecin is regulated by the IMD pathway (Schlüns & Crozier, 2007) and the Toll pathway (Evans et al., 2006), whereas - Blast identification of the complete sequence of the non-structural polyprotein gene of the DWV used for the analysis of the amino acid coding sequence of the RdRp protein fragment.

 

Fig. 3. Amino acid coding sequence analysis of the fragments of the RdRp protein, of the non-structural polyprotein of the DWV, obtained from the nucleotide sequences of AD6, AD7, A9 & A10.

defensin has been reported to be regulated by the Toll pathway (Lourenço et al., 2018). In contrast to our results, Yang and Cox-Foster (2005) described a decrease in the expression of abaecin, defensin-1 and hymenoptaecin in DW bees, relating this immunosuppression with viral amplification. Interestingly, Azzami et al. (2012) reported that bees infected with ABPV produced neither elevated levels of specific antimicrobial peptides, such as hymenoptaecin and defensin, nor any general antimicrobial activity.

In our study no differences were found in NW and DW bees in Vg mRNA levels, which could indicate that this glycolipoprotein is not involved in the immune response induced by DWV.

We also studied the expression of two genes of the Toll pathway by analyzing mRNA levels of Toll wheeler, previously studied in antimicrobial immune defence in A. mellifera(Aronstein & Saldivar, 2005) and cactus, the NF

kappa B inhibitor (IκB). According to Evans et al. (2006), following conformational changes of the activated Toll-like receptor, several intracellular death-domain (DD) containing proteins are recruited to form a receptor complex. Activation of this complex leads to the degradation of the NF kappa B inhibitor (IκB) cactus and subsequent nuclear translocation of the NF-κB transcription factor Dorsal. In this study no differences in Toll wheeler mRNA levels were found, but, in concordance to Nazzi et al. (2012) who showed a down-regulation of a honey bee NF-kB transcription factor, a 3.3-fold induction of cactus mRNA levels was detected in DW bees in comparison with NW bees. These results suggest that this pathway is inactivated against high levels of Varroa infestation and high viral loads of the DWV. This could indicate that the overexpression of hymenoptaecin discussed above could be related with other signalling pathway other than Toll (see Evans et al., 2006) or that other Toll receptors (besides Toll wheeler) should be evaluated.

 

Fig. 4. Phylogenetic tree created using MAFFT software on-line (Katoh et al., 2017) to visualize the clusters of the sequences obtained (AD6, AD7, A9 & A10) with the top hits from the GenBank.

 

The Jak-STAT pathway has been shown to respond to bacterial and viral infections in insects (Kingsolver et al., 2013;Chen et al., 2014). In this work, two members of the Jak-STAT pathway were studied: the receptor domeless and the effector protein TEPA showing a significant up-regulation in DW bees. These results demonstrate that the Jak-STAT pathway is activated in response to Varroa infestation and DWV infection in A. mellifera. The Jak-STAT pathway represents an extremely rapid membrane-to-nucleus signalling system. Interferon-g (Ifin-g) is well known to play crucial roles in several aspects of the immune responses and one action responsible for its infammatory function results from activation of the nitric oxide synthase (NOS) promoter and NOS gene expression. To achieve pleiotropic biologic activities, including NOS induction, signalling cascades initiated from receptors associated with Jak-STAT pathway are critical (Chen et al., 2002). Interestingly, in this study the expression levels of the NOS gene were detected in high level at the DW group.

Nitric oxide (NO, the signalling molecule that results from NOS activity) has been described recently as an effector that mediates cellular and humoral immune response in insects (Sadekuzzaman & Kim, 2018). The involvement of NO in the immune response of A. mellifera has been previously reported (Negri et al., 2013, 2014) and a role for NO in the immune response to viruses has also been described (Reiss & Komatsu, 1998). Increased NOS has also been shown to be related to bacterial infection (Hillyer & Estévez-Lao, 2010) and induced after blood feeding (Lim et al., 2005) in mosquitoes that transmit intracellular pathogens like malarial protozoans. In other insects as Bombyx mori (Lepidoptera) (Ishii et al., 2013) NOS mediates the expression of AMP (Sadekuzzaman & Kim, 2018), and this particularity could be related to the hymenoptaecin up-regulation observed from our results.

In DW bees we found higher titers of DWV compared with NW bees, and an up-regulation of A. melliferaJak-STAT signalling pathway genes and cactus expression level, while the abaecin and defensin expression levels were not significantly different. In insects, the connection between the RNAi and Jak-STAT pathway has been demonstrated with the NF-κB pathways (Toll and Imd). The explanation of AMPs expression decrease could be explained according to Jak-STAT and the Imd/Toll pathways as a negative interaction, as seen in other insects (Kingsolver et al., 2013). This was related to the inhibitory effect of a complex integrated by STAT that competes with Relish (the NF-κB transcription factor from Imd pathway) inducing the negative regulation of certain AMPs (Kim et al., 2007).

This study shows that DW bees are infected with a virulent DWV strain that not only induces changes in immunogene expression indicating activation of the Jak-STAT pathway, but also shows an increased viral load, in comparison with NW bees.Varroa-adapted DWV variants present in our country should continue evolving, and investigations of virus strain differences may explain the different pathologies currently seen in honey bee colonies. Such variants may interact with other pests, pathogens, environmental factors, and regional beekeeping practices, resulting in recent largescale losses of honey bee colonies (Cox- Foster et al., 2007; Maggi et al., 2013). In concordance with this and all the results found in the current study more analyses are required to target the possible cross talk between Jak-STAT and Toll/Imd pathways, and also the possibility of a relation between NOS expression and AMP production. Analysis of the differences in the modulation of immunogens in response to DWV different variants and viral load between summer and winter bees, with different infestation levels o f Varroa, is also needed.

ACKNOWLEDGEMENTS

This research was supported by grant and PICT-2011 No. 0162 to M.M. The authors would like to thank the Institute of Analysis Fares Taie, CONICET and the UNMdP for financial support.

 

LITERATURE CITED

Allen, M.F., & Ball, B.V. (1996) The incidence and world distribution of honeybee viruses. Bee World, 77, 141-162.         [ Links ]

Altschul, S.F., Gish, W., Miller, W., Myers, E.W., & Lipman, D.J. (1990) Basic local alignment search tool. Journal of Molecular Biology, 215, 403-410.         [ Links ]

Antúnez, K., D'Alessandro, B., Corbella, E., Ramallo, G., & Zunino, P. (2006) Honeybee viruses in Uruguay. Journal of Invertebrate Pathology, 93(1), 67-70.         [ Links ]

Antúnez, K., Martín-Hernández, R., Prieto, L., Meana, A., Zunino, P., & Higes, M. (2009) Immune suppression in the honey bee (Apis mellifera) following infection by Nosema ceranae(Microsporidia). Environmental Microbiology, 11(9), 2284-2290.         [ Links ]

Aronstein, K., & Saldivar, E. (2005) Characterization of a honey bee Toll related receptor gene Am18w and its potential involvement in antimicrobial immune defense. Apidologie, 36, 3-14.         [ Links ]

Artimo, P., Jonnalagedda, M., Arnold, K., Baratin, D., Csardi, G., de Castro, E., Duvaud, S., Flegel, V., Fortier, A., et al. (2012) ExPASy: SIB bioinformatics resource portal. Nucleic Acids Research, 40(W1), W597-W603        [ Links ]

Azzami, K., Ritter, W., Tautz, J., & Beier, H. (2012) Infection of honey bees with acute bee paralysis virus does not trigger humoral or cellular immune responses. Archives of Virology, 157, 689-702.         [ Links ]

Barchuk, A.R., Bitondi, M.M.G., & Simões, Z.L.P. (2002) Effects of juvenile hormone and ecdysone on the timing of vitellogenin appearance in hemolymph of queen and worker pupae of Apis mellifera. Journal of Insect Science, 2(1).         [ Links ]

Berényi, O., Bakonyi, T., Derakhshifar, I., Köglberger, H., Topolska, G., Ritter, W., Pechhacker, H., & Nowotny, N. (2007) Phylogenetic Analysis of Deformed Wing Virus Genotypes from Diverse Geographic Origins Indicates Recent Global Distribution of the Virus. Applied and Environmental Microbiology, 73(11), 3605-3611.         [ Links ]

Boncristiani, H.F., Evans, J.D., Chen, Y., Pettis, J., Murphy, C., Lopez, D.L., Simone-Finstrom, M., Strand, M., Tarpy, D.R., & Rueppell, O. (2013) In Vitro Infection of Pupae with Israeli Acute Paralysis Virus Suggests Disturbance of Transcriptional Homeostasis in Honey Bees (Apis mellifera). PLoS ONE, 8(9), e73429.         [ Links ]

Brasesco, C., Quintana, S., Negri, P., Medici, S., Ruffinengo, S., Eguaras, M., & Maggi, M. (2013) Detección mediante PCR en Tiempo Real de virus patógenos de Apis mellifera en el sudeste de la provincia de Buenos Aires. In: Acta Resúmenes del VIII Encuentro Anual Biólogos en Red, 2013, p. 70.         [ Links ]

Brutscher, L.M., Daughenbaugh, K.F., & Flenniken, M.L. (2015) Antiviral Defense Mechanisms in Honey Bees. Current Opinion in Insect Science, 10, 71-82.         [ Links ]

Chen, C., Chao, Y., Chang, Y., Hsu, M., & Lin, W. (2002) Inhibition of cytokine-induced Jak-STAT signalling pathways by an endonuclease inhibitor aurintricarboxylic acid. British Journal of Pharmacology, 137, 1011-1020.         [ Links ]

Chen, Y.P., Pettis, J.S., Corona, M., Chen, W.P., Li, C.J., Spivak, M., Visscher, P.K., Degrandi-Hoffman, G., Boncristiani, et al. (2014) Israeli acute paralysis virus: Epidemiology. Pathogenesis and implications for honey bee health. PLoS Pathogens, 10, e1004261.         [ Links ]

Cox-Foster, D.L., Conlan, S., Holmes, E.C., Palacios, G., Evans, J.D., Moran, N.A., Quan, P.L., Briese, T., Hornig, M., et al. (2007) A metagenomic survey of microbes in honey bee colony collapse disorder. Science, 318, 283-286.         [ Links ]

Dalmon, A., Desbiez, C., Coulon, M., Thomasson, M., Le Conte, Y., Alaux, C., Vallon, J., & Moury, B. (2017) Evidence for positive selection and recombination hotspots in Deformed wing virus (DWV). Scientific reports, 7, 41045.         [ Links ]

de Miranda, J.R., Bailey, L., Ball, B.V., Blanchard, P., Budge, G.E., Chejanovsky, N., Chen, Y.P., Gauthier, L., Genersch, E., et al. (2013) Standard methods for virus research in Apis mellifera. Journal of Apicultural Research, 52(4), 1-56.         [ Links ]

Eguaras, M., & Ruffinengo, S. (2006) Estrategias para el control de Varroa. Editorial Martin. Mar del Plata, Argentina.

Evans, J.D. (2006) Beepath: an ordered quantitative-PCR array for exploring honey bee immunity and disease. Journal of Invertebrate Pathology, 93(2), 135-139.

Evans, J.D., Aronstein, K., Chen, Y.P., Hetru, C., Imler, J-L., Jiang, H., Kanost, M., Thompson, G.J., Zou, Z., & Hultmark, D. (2006) Immune pathways and defence mechanisms in honey bees Apis mellifera. Insect Molecular Biology, 15(5), 645-656.

Finn, R.N. (2007) Vertebrate yolk complexes and the functional implications of phosvitins and other subdomains in vitellogenins. Biology of Reproduction, 76, 926-935.

Genersch, E. (2010) Honey bee pathology: current threats to honey bees and beekeeping. Applied Microbiology and Biotechnology, 87, 87-97.

Genersch, E., & Aubert, M. (2010) Emerging and re-emerging viruses of the honey bee (Apis mellifera L.). Veterinary Research, 41(6), 54.

Gisder, S., Aumeier, P., & Genersch, E. (2009) Deformed wing virus (DWV): viral load and replication in mites (Varroa destructor). Journal of General Virology, 90, 463-467.

Gregorc, A., Evans, J.D., Scharf, M., & Ellis, J.D. (2012) Gene expression in honey bee (Apis mellifera) larvae exposed to pesticides and Varroa mites (Varroa destructor). Journal of Insect Physiology, 58, 1042-1049.

Gregory, P.G., Evans, J.D., Rinderer, T., & de Guzman, L. (2005) Conditional immune-gene suppression of honeybees parasitized by Varroamites. Journal of Insect Science, 5(1).

Hendriksma, H.P., Härtel, S., & Steffan‐Dewenter, I. (2011) Honey bee risk assessment: new approaches for in vitro larvae rearing and data analyses. Methods in Ecology and Evolution, 2(5), 509-517.

Hillyer, J.F., & Estévez-Lao, T.Y. (2010) Nitric oxide is an essential component of the hemocyte-mediated mosquito immune response against bacteria. Developmental & Comparative Immunology, 34(2), 141-149.

Invernizzi, C., Antúnez, K., Campa, J.P., Harriet, J., Mendoza, Y., Santos, E., & Zunino, P. (2011) Situación sanitaria de las abejas melíferas en Uruguay. Veterinaria, 47, 15-27

Ishii, K., Adachi, T., Hamamoto, H., Oonishi, T., Kamimura, M., Imamura, K., & Sekimizu K. (2013) Insect cytokine paralytic peptide activates innate immunity via nitric oxide production in the silkworm Bombyxmori. Developmental&Comparative Immunology, 39, 147-153.

Katoh, K., Rozewicki, J., & Yamada, K. D. (2017) MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Briefings in bioinformatics, bbx108.

Kim, L.K., Choi, U.Y., Cho, H.S., Lee, J.S., Lee, W.B., Kim, J., Jeong, K., Shim, J., Kim-Ha, J., & Kim, Y.J. (2007) Downregulation of NF-kappaB target genes by the AP-1 and STAT complex during the innate immune response in Drosophila. PLoS biology, 5, e238.

Kingsolver, M.B., Huang, Z., & Hardy, R.W. (2013) Insect Antiviral Innate Immunity: Pathways, Effectors, and Connections. Journal of Molecular Biology, 425(24), 4921-4936.

Kukielka, D., Esperón, F., Higes, M., & Sánchez-Vizcaíno, J.M. (2008) A sensitive one-step real-time RT-PCR method for detection of deformed wing virus and black queen cell virus in honeybee Apis mellifera. Journal of Virological Methods, 147(2), 275-281.

Lanzi, G., de Miranda, J., Boniotti, M.B., Cameron, C.E., Lavazza, A., Capucci, L., Camazine, S.M., & Rossi, C. (2006) Molecular and Biological Characterization of Deformed Wing Virus of Honeybees (Apis mellifera L.). Journal of Virology, 80(10), 4998-5009.

Le Conte, Y., Ellis, M., & Ritter, W. (2010) Varroa mites and honey bee health: can Varroa explain part of the colony losses? Apidologie, 41, 353-363.

Lim, J., Gowda, D.C., Krishnegowda, G., & Luckhart, S. (2005) Induction of nitric oxide synthase in Anopheles stephensi by Plasmodium falciparum: mechanism of signaling and the role of parasite glycosylphosphatidylinositols. Infection and Immunity, 73, 2778-2789.

Locke, B., Forsgren, E., Fries, I., & de Miranda, J.R. (2012) Acaricide treatment affects viral dynamics in Varroa destructor-infested honey bee colonies via both host physiology and mite control. Applied and Environmental Microbiology, 78(1), 227-235.

Lourenço, A.P., Florecki, M.M., Simões, Z.L.P., & Evans, J.D. (2018) Silencing of Apis melliferadorsal genes reveals their role in expression of the antimicrobial peptide defensin‐1. Insect molecular biology, 27(5), 577-589.

Maggi, M., Ruffinengo, S., Negri, P., Brasesco, C., Medici, S., Quintana, S., Szawarski, N., Gimenez Martinez, P., De Piano, F., et al. (2013) The status of bee health and colony losses in Argentina. Honeybees: Foraging Behavior, Reproductive Biology and Diseases (ed. Cameron Molley), pp. 212-234. Nova Publishing Group.

Marin, M.S., Quintana, S., Leunda, M.R., Recavarren, M., Pagnuco, I., Späth, E., Pérez, S., & Odeón, A. (2016) A new method for simultaneous detection and discrimination of Bovine herpesvirus types 1 (BoHV-1) and 5 (BoHV-5) using real time PCR with high resolution melting (HRM) analysis. Journal of Virological Methods, 227, 14-22.

Marringa, W.J., Krueger, M.J., Burritt, N.L., & Burritt, J.B. (2014) Honey bee hemocyte profling by fow cytometry. PloS ONE, 9(10), e108486.

Martin, S., Highfeld, A., Brettell, L., Villalobos, E., Budge, G., Powell, M., Nikaido, S., & Schroeder, D. (2012) Global Honey Bee Viral Landscape Altered by a Parasitic Mite. Science, 336(6086), 1304-1306.

Merkling, S.H., & van Rij, R.P. (2013) Beyond RNAi: antiviral defense strategies in Drosophila and mosquito. Journal of Insect Physiology, 59(2), 159-170.

Nappi, A.J., Vass, E., Frey, F., & Carton, Y. (2000) Nitric oxide involvement in Drosophila immunity. Nitric Oxide, 4(4), 423-430.

Navajas, M., Migeon, A., Alaux, C., Martin-Magniette, M.L., Robinson, G.E., Evans, J.D., Cros-Arteil, S., Crauser, D., & Le Conte, Y. (2008) Differential gene expression of the honey bee Apis melliferaassociated with Varroa destructor infection. BMC Genomics, 9(1), 301-311.

Nazzi, F., Brown, S.P., Annoscia, D., Del Piccolo, F., Di Prisco, G., Varricchio, P., Della Vedova, G., Cattonaro, F., Caprio, E., & Pennacchio, F. (2012) Synergistic parasite pathogen interactions mediated by host immunity can drive the collapse of honeybee colonies. PLoS Pathogens, 8(6), e1002735.

Negri, P., Maggi, M., Correa-Aragunde, N., Brasesco, C., Eguaras, M., & Lamattina, L. (2013) Nitric oxide participates at the frst steps of Apis melliferacellular immune activation in response to non-self recognition. Apidologie, 44(5), 575-585.

Negri, P., Quintana, S., Maggi, M., Szawarski, N., Lamattina, L., & Eguaras, M. (2014) Apis melliferahemocytes generate increased amounts of nitric oxide in response to wounding/ encapsulation. Apidologie, 45(5), 610-617.

Pfaff, M., Horgan, G., & Demplfe, L. (2002) Relative expression software tool (REST) for group-wise comparison and statistical analysis of relative expression results in real time PCR. Nucleic Acids Research, 30(9), e36-e36.

Potts, S., Biesmeijer, J., Kremen, C., Neumann, P., Schweiger, O., & Kunin, W. (2010) Global pollinator declines: trends, impacts and drivers. Trends in Ecology&Evolution, 25(6), 345-353.

Reiss, C.S., & Komatsu, T. (1998) Does Nitric Oxide Play a Critical Role in Viral Infections? Journal of Virology, 72(6), 4547-4551.

Reynaldi, F.J., Sguazza, G.H., Pecoraro, M.R., Tizzano, M.A., & Galosi, C.M. (2010) First report of viral infections that affect Argentine honey bees. Environmental Microbiology Reports, 2, 749-751.

Reynaldi, F.J., Sguazza, G.H., Tizzano, M.A., Fuentealba, N.A., Galosi, C.M., & Pecoraro, M.R. (2011) First report of Israeli acute paralysis virus in asymptomatic hives of Argentina. Revista Argentina de Microbiología, 43, 84-86.

Runckel, C., Flenniken, M.L., Engel, J.C., Ruby, J.G., Ganem, D., Andino, R., & DeRisi, J.L. (2011) Temporal analysis of the honey bee microbiome reveals four novel viruses and seasonal prevalence of known viruses, Nosema, and Crithidia. PloS ONE, 6(6), e20656.

Ryabov, E.V., Wood, G.R., Fannon, J.M., Moore, J.D., Bull, J.C., Chandler, D., Mead, A., Burroughs, N., & Evans, D.J. (2014) A virulent strain of deformed wing virus (DWV) of honeybees (Apis mellifera) prevails after Varroa destructor-mediated, or in vitro, transmission. PLoS Pathogens, 10(6), e1004230.

Sadekuzzaman, M., & Kim, Y. (2018) Nitric oxide mediates antimicrobial peptide gene expression by activating eicosanoid signaling. PLoS ONE, 13(2), e0193282.

Salmela, H., Amdam, G.V., & Freitak, D. (2015) Transfer of Immunity from Mother to Offspring Is Mediated via Egg-Yolk Protein Vitellogenin. PLoS Pathogens, 11(7), e1005015.

Schlüns, H., & Crozier, R.H. (2007) Relish regulates expression of antimicrobial peptide genes in the honeybee, Apis mellifera, shown by RNA interference. Insect molecular biology, 16(6), 753-759.

Siede, R., Meixner, M., & Büchler, R. (2012) Comparison of transcriptional changes of immune genes to experimental challenge in the honey bee (Apis mellifera). Journal of Apicultural Research, 51(4), 320-328.

Simone, M., Evans, J.D., & Spivak, M. (2009) Resin collection and social immunity in honey bees. Evolution, 11, 3016-3022.

Takeda, K., & Akira, S. (2003) Toll receptors and pathogen resistance. Cellular Microbiology, 5(3), 143-153.

vanEngelsdorp, D., Underwood, R., Caron, D., & Hayes Jr., J. (2007) An estimate of managed colony losses in the winter of 2006-2007: a report commissioned by the Apiary Inspectors of America. The American Bee Journal, 147, 599-603.

vanEngelsdorp, D., Hayes Jr., J., Underwood, R.M., & Pettis, J. (2008) A survey of honey bee colony losses in the U.S., fall 2007 to spring 2008. PLoS ONE, 3, e4071.

vanEngelsdorp, D., Evans, J.D., Saegerman, C., Mullin, C., Haubruge, E., Nguyen, B.K., Frazier, M., Frazier, J., Cox-Foster, D., et al. (2009) Colony Collapse Disorder: A Descriptive Study. PLoS ONE, 4(8), e6481.

Wang, H., Xie, J., Shreeve, T.G., Ma, J., Pallett, D.W., King, L.A., & Possee, R.D. (2013) Sequence recombination and conservation of Varroa destructor virus-1 and deformed wing virus in feld collected honey bees (Apis mellifera). PloS ONE, 8(9), e74508.

Yang, X., & Cox-Foster, D.L. (2005) Impact of an ectoparasite on the immunity and pathology of an invertebrate: Evidence for host immunosuppression and viral amplification. Proceedings of the National Academy of Sciences of the United States of America, 102(21), 7470-7475.

Zhang, S., Wang, S., Li, H., & Li, L. (2011) Vitellogenin, a multivalent sensor and an antimicrobial effector. International Journal of Biochemistry&Cell Biology, 43, 303-305.