- Open Access
RNA interference-based antiviral immune response against the salivary gland hypertrophy virus in Glossina pallidipes
BMC Microbiology volume 18, Article number: 170 (2018)
Glossina pallidipes salivary gland hypertrophy virus (GpSGHV; Hytrosaviridae) is a non-occluded dsDNA virus that specifically infects the adult stages of the hematophagous tsetse flies (Glossina species, Diptera: Glossinidae). GpSGHV infections are usually asymptomatic, but unknown factors can result to a switch to acute symptomatic infection, which is characterized by the salivary gland hypertrophy (SGH) syndrome associated with decreased fecundity that can ultimately lead to a colony collapse. It is uncertain how GpSGHV is maintained amongst Glossina spp. populations but RNA interference (RNAi) machinery, a conserved antiviral defense in insects, is hypothesized to be amongst the host’s mechanisms to maintain the GpSGHV in asymptomatic (persistent or latent) infection state. Here, we investigated the involvement of RNAi during GpSGHV infections by comparing the expression of three key RNAi machinery genes, Dicer (DCR), Argonaute (AGO) and Drosha, in artificially virus injected, asymptomatic and symptomatic infected G. pallidipes flies compared to PBS injected (controls) individuals. We further assessed the impact of AGO2 knockdown on virus infection by RT-qPCR quantification of four selected GpSGHV genes, i.e. odv-e66, dnapol, maltodextrin glycosyltransferase (a tegument gene) and SGHV091 (a capsid gene).
We show that in response to hemocoelic injections of GpSGHV into G. pallidipes flies, increased virus replication was accompanied by significant upregulation of the expression of three RNAi key genes; AGO1, AGO2 and DCR2, and a moderate increase in the expression of Drosha post injection compared to the PBS-injected controls. Furthermore, compared to asymptomatically infected individuals, symptomatic flies showed significant downregulation of AGO1, AGO2 and Drosha, but a moderate increase in the expression of DCR2. Compared to the controls, knockdown of AGO2 did not have a significant impact on virus infection in the flies as evidenced by unaltered transcript levels of the selected GpSGHV genes.
The upregulation of the expression of the RNAi genes implicate involvement of this machinery in controlling GpSGHV infections and the establishment of symptomatic GpSGHV infections in Glossina. These findings provide a strategic foundation to understand GpSGHV infections and to control latent (asymptomatic) infections in Glossina spp. and thereby control SGHVs in insect production facilities.
Tsetse flies (Diptera: Glossinidae) are naturally infected by the Glossina pallidipes salivary gland hypertrophy virus (GpSGHV, family Hytrosaviridae), a large double-stranded DNA (dsDNA) virus pathogenic specifically to Glossina spp., [1, 2]. In G. pallidipes, although the majority of cultured and wild tsetse fly species are asymptomatically infected by GpSGHV (low virus titers), some unknown factors can trigger symptomatic infections (high virus titers) . This, in turn, is associated with the occurrence of overt salivary gland hypertrophy (SGH) symptoms [4, 5]. In the mass rearing of G. pallidipes, SGH epizootics reduce fly survival and productivity, and have caused the collapse of three colonies; two in the Insect Pest Control Laboratory (IPCL) in Seibersdorf, Austria (in 1987 and 2001), and one in the mass rearing facility in Kality, Ethiopia (in 2012) . These GpSGHV-induced effects have significantly compromised the implementation of the sterile insect technique (SIT), a component of area-wide integrated pest management (AW-IPM) strategies designed for the eradication of G. pallidipes from the Southern Rift Valley of Ethiopia .
Although it is uncertain how the virus is maintained within tsetse populations in nature and in laboratory colonies, three hypothetical scenarios may account for the maintenance of asymptomatic GpSGHV infection state. The first is a persistent infection whereby the virus remains in specific host cells with low-level production of virions, but without causing substantial cell damage . The second is a latent infection state, during which viral genomes and maybe some viral proteins are present in the infected host cells of certain organs, but without detectable production of infectious viral particles . In the third, the virus can exist in both persistent and latent infection states at the same time, but in different tissues . Persistent infection in the salivary gland (SG) cells is accompanied by a low number of virions (102 viral genome copies/fly) released by asymptomatic flies via saliva during feeding . In addition, detection of viral DNA in other tissues such as the tracheal cells without detectable viral gene transcripts  may reflect a latent infection state. In any case, the persistent or latent GpSGHV infection in G. pallidipes potentially represents a homeostatic equilibrium between the host’s immune system and the viral escape strategies. Consequently, the viral infection is kept under control (asymptomatic state), but is not completely eliminated from the fly.
Amongst the possible host mechanisms that keep GpSGHV under control is the insect’s RNA interference (RNAi) machinery, which regulates both host and viral gene expression by use of small RNAs that bind to their complementary messenger RNA (mRNA) targets [12, 13]. This hypothesis is based on evidence from various studies indicating that the RNAi machinery is a conserved antiviral defense mechanism for several groups of large dsDNA viruses infecting insects, including the related baculoviruses and nudiviruses, as well as ascoviruses and iridoviruses [14, 15]. RNAi is mediated through three pathways: short interfering RNA (siRNA), microRNA (miRNA) and Piwi-interacting RNA (piRNA) pathways . The siRNAs are processed in the cytoplasm by the ribonuclease III enzyme Dicer-2 (DCR2) from exogenous double stranded (ds) RNAs (e.g. dsRNAs that arise as viral replication intermediates or from overlapping transcripts). These siRNAs destroy (viral or cellular) single-stranded RNAs (ssRNAs) in a sequence-specific manner . The miRNAs, on the other hand, are processed by DCR1 from cellular or viral pre-miRNAs, originating from DNA components of nuclear replicating viruses that are processed by RNase III enzyme Drosha in the nucleus . The miRNAs are then exported to the cytoplasm where they post-transcriptionally regulate cellular or viral protein expression, thereby modulating developmental and physiological processes of the host as well as virus infection [19, 20]. The piRNAs are processed via diverse pathways independent of DCR proteins , and are involved in the regulation of cellular genes and activities of transposons . They may also have a role in antiviral strategies, as has been suggested for arboviruses . In the RNAi biogenesis pathways, the siRNAs, miRNAs and piRNAs are loaded into Argonaute proteins 2, 1 and 3 (AGO2, 1 and 3), respectively, that mediate the RNAi by either cleavage or degradation of target RNAs (AGO2), translation repression (AGO1), or epigenetic modifications (AGO3) . However, some of the above-mentioned enzymes may participate in two or more of these pathways. For instance, in Drosophila melanogaster, DCR1 is involved in both siRNA and miRNA pathways, while DCR2 is only involved in the siRNA pathway. Additionally, Drosophila AGO1 and AGO2 enzymes may participate in both the siRNA and miRNA pathways .
The siRNA-mediated RNAi pathway is a potent antiviral immune pathway in insects [21, 26] and is implicated in controlling the replication of RNA and DNA viruses [15, 27,28,29,30]. In addition, it has been shown for several viruses that the knockdown of RNAi pathway components leads to increased viral replication. For instance, loss-of-function mutations in DCR2 enhanced the susceptibility of Helicoverpa armigera to infection by H. armigera single nucleopolyhedrovirus (HearNPV) . Similar observations were made in Drosophila during infection by Flock House virus (FHV), Drosophila C virus (DCV), and Sindbis virus (SINV) . In addition to controlling viral replication, the siRNA pathways have also been implicated in establishing persistent virus infection . To establish persistent infection, viral fragments generated during viral DNA genome replication are transcribed or reverse transcribed for RNA viruses and integrated into the host genome. The generated transcripts are processed by DCR2 into virus-derived siRNAs (vsiRNAs); the vsiRNAs are then loaded into AGO2 to mediate specific cleavage of viral mRNAs, leading to persistent infection [32,33,34]. In the case of GpSGHV infections, the outcome of RNAi-based immune responses would hypothetically be restriction of viral replication and prevention of the development of overt SGH. If this is indeed the case, then virus and host would progress into a stable equilibrium of a persistent or latent infection state, which may account for the widespread chronic asymptomatic GpSGHV infections in many tsetse species, particularly in colonized flies .
In the current study, we evaluated whether GpSGHV infection induces an RNAi response in G. pallidipes and whether this would downregulate the development of SGH and instead induce a covert infection state (persistent or latent). To accomplish this, we performed comparative analyses of the expression of AGO, DCR and Drosha between artificially (intra-hemocoelic) virus injected and uninfected (PBS injected) individuals, and between asymptomatic and symptomatic infected flies (with overt SGH symptoms). It should be noted that artificial injection of the virus does not result into overt SGH in the same (parental) generations but in the progeny flies . We complemented these bioassays by testing the impact of downregulation of a key component in the siRNA pathway (AGO2), on GpSGHV infection in G. pallidipes. The data obtained in this study offer a rationale for similar studies on other Hytrosaviridae family members and may open novel strategies to manage SGHVs in insect production facilities.
The Argonaute family in Glossina species
The analyses of the genomes of G. pallidipes, G. m. morsitans, G. f. fuscipes, G. p. palpalis, G. austeni and G. brevipalpis resulted in the identification of AGO 1, 2 and 3 in all these species (See Table 1), key components of the RNA induced silencing complex (RISC); AGOs activate and cleave target mRNA within the RISC complex . The identification of AGOs in the six tsetse species underscores the conservation of the RNAi machinery in Glossina species. Phylogenetically, the three AGO proteins segregated into distinct clusters with their orthologs in D. melanogaster, which corresponded to the siRNA, microRNA and piRNA pathways of the RNAi machinery (Fig. 1a). The phylogenetic clustering was supported by robust bootstrap values. Additionally, similar to the D. melanogaster AGO family proteins, their orthologs in Glossina species contained the critical functional domains, i.e. the PAZ domain (for dsRNA binding) and the PIWI domain (executioner of the RNase activity) (Fig. 1b). These results strongly suggest that the three RNAi machinery pathways are functional in Glossina spp.
The dicer family in Glossina species
The bioinformatics analysis of the DCR protein family did not result in the identification of the homolog to the Drosophila DCR1 protein in G. pallidipes and G. brevipalpis, but the other Glossina species included in this study contained a DCR1 protein homolog (Fig. 2a). However, homologs to the Drosophila DCR2 protein were present in both G. pallidipes and G. brevipalpis, suggesting that DCR2 might be involved in both siRNA and miRNA pathways, at least in these species. However, DCR2 was lacking in G. p. palpalis and G. austeni (Fig. 2a), implying that in these two species DCR1 might be involved in both siRNA and miRNA pathways. Homologs to both DCR1 and DCR2 were found only in G. m. morsitans and G. f. fuscipes suggesting that in these species they may be involved in two separate pathways (e.g. miRNA and siRNA) as reported in Drosophila . Drosha was present in all six tsetse species investigated (Fig. 2a). Our bioinformatics analysis revealed the presence of all the functional motifs in the identified DCR (N-terminal helicases, DCR- dsRBF, PAZ, two C-terminal RNA III, and dsRBD) and Drosha (C-terminal RNA III and the dsRBD) protein homologs, which were organized as in their orthologs in D. melanogaster (Fig. 2b). Notably, no major differences (phylogenetic and domain architecture) were observed between the DCR1 or DCR2 protein sequences of the tsetse species containing one or both DCR proteins (Fig. 2a and b). The presence of the functional domains in DCR and Drosha homologs could imply the conservation and functionality of both siRNA and miRNA pathways in tsetse. As indicated above, the piRNA pathway is independent of DCR implying that function of this pathway may not be affected by the presence or function of this protein.
Expression levels of AGO, DCR and Drosha in virus-injected G. pallidipes
Having identified the core genes involved in the RNAi machinery pathways, we determined whether GpSGHV infection induces an RNAi response. Using G. pallidipes as the model species, and due to the absence of DCR1 in this Glossina species (See Fig. 2), we quantified the expression levels of the AGO family (AGO1, AGO2 and AGO3), DCR2, and Drosha in flies injected with the virus suspensions compared to the PBS-injected control flies. Additionally, we correlated the expression levels of these RNAi-related genes to the level of virus replication by quantifying the expression levels of GpSGHV odv-e66, a conserved and late viral gene. Compared to the control (PBS-injected) fly group, the expression of odv-e66 increased significantly with time (t = 8.657; d.f. = 44; P < 0.001) (Fig. 3a), implying active replication and late gene expression of the virus in the injected flies. This increased virus replication was accompanied by a significant increase in the expression level of AGO1 (t = 2.306; d.f. = 44; P = 0.026) and AGO2 (t = 3.334; d.f. = 44; P = 0.00174) but not AGO3 (t = 1.651; d.f. = 44; P = 0.106), of which the AGO2 (involved in siRNA pathway) was the most upregulated (compare panels B, C and D in Fig. 3). Similarly, to the AGO genes DCR2, which may be involved in both siRNA and miRNA pathways in G. pallidipes, was also found to be significantly upregulated (t = 3.968; d.f. = 44; P < 0.001) in response to the virus injection (Fig. 3e). However, unlike AGO and DCR, the expression levels of Drosha showed no significant increase (t = 0.601; d.f. = 44; P = 0.551) in the virus-injected flies compared to the levels observed in the PBS-injected flies (Fig. 3f). This suggests that Drosha, part of the miRNA pathway, is not involved in the immune response against a lytic infection by GpSGHV. Full statistics for each regression are given in Additional file 1.
Expression levels of AGO, DCR and Drosha in symptomatic G. pallidipes
We also compared the expression levels of the DCR and AGO family genes in the virus-injected G. pallidipes flies described above with the expression levels in symptomatic (flies with overt SGH symptoms and high virus titers) and asymptomatically infected individuals (flies with low virus titers) (t = 16.72; d.f. = 10; P < 0.001) (Fig. 4a). We found that AGO1 (t = − 5.454; d.f. = 10; P < 0.001), AGO2 (t = − 3.899; d.f. = 10; P = 0.00363) and Drosha (t = − 3.549; d.f. = 10; P = 0.00623) were significantly downregulated in symptomatic G. pallidipes flies as compared to asymptomatically infected flies (Fig. 4b). There was no difference in DCR2 (t = 1.318; d.f. = 10; P = 0.2202), or AGO3 (t = − 0.858; d.f. = 10; P = 0.413) expression between the asymptomatic and symptomatic infected flies (Fig. 4b).
Effect of AGO2 knockdown on virus infection in G. pallidipes
We also assessed the impact of AGO2 knockdown on GpSGHV infection. The AGO2 gene, which is involved in the siRNA pathway, was chosen for knockdown largely because its expression levels were significantly modulated in both the virus-injected flies (upregulated; see Fig. 3c) and symptomatic infected flies (downregulated; See Fig. 4b). Notably, compared to the flies injected with nuclease-free water, injection of flies with the dsRNAs did not cause any difference in mortality rate. The injection of dsRNAs specific for AGO2 and tsetse EP in addition to PBS or virus injection resulted in a significant decrease in the expression levels of both AGO2 (AGO2dsRNA/PBS cf. water/PBS: t = − 4.265; d.f. = 42; P < 0.001, AGO2dsRNA/virus cf. water/virus: t = − 3.543 d.f. = 42; P < 0.001) (Fig. 5a) and tsetse EP (TsetseEPdsRNA/PBS cf. water/PBS: t = − 5.392; d.f. = 40; P < 0.001, TsetseEPdsRNA/virus cf. water/virus: t = − 6.798; d.f. = 40; P < 0.001) (Fig. 5b), respectively, compared to the water-injected control flies. We then assessed the effect of AGO2 knockdown on virus infection by quantifying expression levels of the selected viral genes, odv-e66, DNApol, SGHV038 and SGHV091. AGO2 knockdown did not have a significant impact on the transcript levels of any of the selected viral genes; odv-e66 (t = − 1.861; d.f. = 119; P = 0.391), DNApol (t = − 0.422; d.f. = 119; P = 0.674), SGHV038 (t = − 0.179; d.f. = 119; P = 0.858) and SGHV091 (t = − 0.877; d.f. = 119; P = 0.382) compared to the controls (Fig. 6a, b, c and d). As expected, knockdown of tsetse EP, which is not associated with the RNAi machinery, did not affect the expression levels of these selected viral genes. The transformations and parameters for all regression analyses are given in Additional file 1.
In the current study, we have investigated the potential involvement of the RNAi machinery during GpSGHV infections in G. pallidipes by quantifying the expression of both the host (AGO, DCR and Drosha) and four viral genes. The finding that AGO and DCR proteins in Glossina species contain the functional domains or motifs known to mediate the RNAi response strongly indicate that the RNAi machinery is functional in tsetse, presumably in a similar version as reported in Drosophila . It is known that the presence or absence of these functional domains in AGO or DCR proteins affects the efficiency of the RNAi response. For instance, Gu et al.  discovered that although AGO2 lacking the PAZ domain interacts with duplex siRNAs, the truncated protein was unable to unwind the siRNAs or eject the passenger RNA strands. The passenger strand should be degraded or ejected from the siRNA duplex to allow the guide strand to be incorporated into the RISC complex and target the mRNA . Moreover, the results from our phylogenetic analysis confirmed that the predicted AGO proteins from the various tsetse species analyzed clustered with those of Drosophila.
The presence of only DCR2 in G. pallidipes and G. brevipalpis, and only DCR1 in G. p. palpalis and G. austeni suggest that in these species only one of the respective proteins is involved in the siRNA and miRNA pathways. In contrast, genes for both DCR1 and DCR2 proteins were found in G. m. morsitans and G. f. fuscipes. Given that in Drosophila DCR1 and DCR2 are involved in the miRNA and siRNA pathways, respectively , it is possible that these two proteins are involved in these pathways in these two Glossina species as well. Alternatively, they might be involved in both pathways, as suggested for the other Glossina species as there were no differences in DCR1 or DCR2 protein sequences of species with one or both DCR proteins, which could be linked to the RNAi pathways. The DCR proteins can be involved in both siRNA and miRNA pathways  or be involved in the separate pathways . Taken together, the identification of genes for AGO and DCR proteins may, in principle, be an indication of a robust RNAi silencing response in Glossina species [40, 41].
Currently, the mechanism(s) permitting the covert (asymptomatic) GpSGHV infection in G. pallidipes, and the reactivation from viral persistence/latency to overt symptomatic infection that is associated with overt SGH symptoms , are poorly understood. Our finding of significant upregulation of the expression of AGO2 in GpSGHV-injected flies suggests that the virus infection induces the host’s siRNA-mediated response, presumably to inhibit the virus infection. Notably DCR2, which might be involved in both siRNA and miRNA pathways in G. pallidipes, was similarly upregulated during virus infection. The upregulated expression of both AGO2 and DCR2 post virus injection in G. pallidipes, which are key components in the dsRNA-mediated gene silencing in several insects, suggests a functional RNAi-mediated innate immunity response in Glossina species. However, more work is required to elucidate the precise details of this pathway in Glossina as well as the involvement of the RNAi machinery in other members of the Hytrosaviridae family.
In contrast to the above-mentioned increased levels of AGO2 and DCR2 in virus-injected flies, the comparative analysis of the expression of the two genes between asymptomatic and symptomatic infections showed a significant downregulation of the expression of AGO2, but insignificant upregulation of DCR2 in the flies exhibiting diagnostic SGH symptoms. The high expression of siRNA pathway genes in the virus-injected flies suggests a tight control of the virus via the RNAi response during asymptomatic GpSGHV infections. However, during symptomatic infections as evidenced by increased virus titers, the siRNA pathway may be compromised (as supported by the low expression of AGO2) enabling the virus to escape the RNAi-mediated innate immunity, thereby increasing virus titers and in turn causing the detectable SGH symptoms. A similar outcome has been documented in the case of the African malaria mosquito, Anopheles gambiae, where dsRNA-mediated silencing of AGO2, which functions in conjunction with DCR2 in this mosquito, resulted in increased O’nyong-nyong virus (ONNV) viral loads . It should be noted that, due to its involvement in both the siRNA and miRNA pathways, DCR2 was not considered a suitable candidate for the knockdown assays in this study. Therefore, we assessed the potential involvement of dsRNA-mediated gene silencing (siRNA pathway) in G. pallidipes by knockdown of AGO2.
Several examples of the function of RNAi in insects have been clearly demonstrated, including for species from the orders Diptera, Dictyoptera, Isoptera, Hymenoptera and Orthoptera [43, 44]. In the case of tsetse flies, Walshe et al.,  showed that micro-injection of dsRNA into 6–8 day-old G. m. morsitans flies could persistently silence expression of tsetse EP, a gene that is demonstrated to protect the fly against establishment of trypanosome infections in the midgut . In the current study, the knockdown of AGO2 in G. pallidipes did not alter the transcript levels of the selected GpSGHV genes implying that AGO2 knockdown had no effect on the GpSGHV infection. These findings are in contrast to previous results; for instance in Drosophila melanogaster in which flies that were deficient in the DCR2 protein showed increased susceptibility to infection by members of three different RNA virus families; i.e. FHV (Nodaviridae), DCV (Dicistroviridae), and SINV (Togaviridae) . In general, many studies have led to the conclusion that flies that contain mutations in genes that encode components of the siRNA pathway (including DCR2 and AGO2) or the Janus kinase/signal transducers and activators of transcription (JAK-STAT) pathway, are not only more sensitive to infection by several viruses, but also harbor higher viral titers than their wild-type counterparts. The JAK-STAT pathway is also a conserved insect innate immune antiviral response [13, 26, 31, 38, 47]. In another report, D. melanogaster defective for the AGO2 were found to be hypersensitive to infections by DCV, an infection which also supported a 1000-fold increased production of progeny virus .
The results obtained from the current study showed that reduction of RNAi efficiency in G. pallidipes did not cause a significant impact on the fly immunocompetence. Notably, in addition to the RNAi pathway investigated in this study, the JAK-STAT, immune deficiency (imd) and Toll immune pathways are also involved in elimination of viruses [48,49,50]. The presence of these innate immune pathways may indicate that following the interruption of the RNAi pathway in G. pallidipes flies in our study alternative antiviral pathways in these flies were able to control the virus infections.
Given the high diversity of pathogens, their elimination by host organisms is challenging and therefore many organisms, including insects, employ multiple mechanisms to deal with them. The RNAi machinery, particularly the siRNA pathway, plays a central role in insects by specifically recognizing and eliminating invading pathogens and other invasive elements such as transposons. This study has elucidated important milestones in the infection of GpSGHV in Glossina spp. We detected RNAi key genes in all the analyzed Glossina species, which may indicate a functional antiviral role of RNAi machinery in tsetse flies. In G. pallidipes, the model Glossina species in this study, the siRNA pathway genes AGO2 and DCR2 were upregulated during virus infection, which confirmed the involvement of the RNAi response in the flies’ defense against GpSGHV. We also noted that in G. pallidipes, the siRNA pathway was compromised during symptomatic infection as evidenced by the low expression of AGO2. However, although knockdown of AGO2 in G. pallidipes did not have an impact on virus infection, it would be worthwhile to further investigate the long-term effects of the gene knockdowns on GpSGHV transmission and the expression of SGH symptoms in F1 progeny produced by the parents with the knockdown. The F1 generation is of interest since induction of SGH symptoms does not occur in the parental generation of flies that are intra-hemocoelically injected with the virus but are observed in the subsequent F1 generation produced by injected mothers . It may therefore be important to determine whether the dsRNA-mediated gene silencing is a heritable trait, as has been demonstrated in other studies . The next question would be to determine which factors (exogenous or endogenous) trigger the transition from the asymptomatic to symptomatic state in colonized and wild flies. Together, the approach and results of the current study offer a starting point for further investigations into Hytrosaviridae family members and may open novel strategies to control SGHVs in insect production facilities.
Tsetse fly and virus injections
The G. pallidipes flies used in this study were obtained from the colony maintained at the Joint FAO/IAEA IPCL, Seibersdorf, Austria. G pallidipes was selected for this study because, unlike other species, it shows both asymptomatic and symptomatic virus infections . The experimental G. pallidipes flies were maintained in controlled insectaria with 70–80% relative humidity, 24 ± 1 °C temperature and 12 h photo-phase. The flies were fed for 10–15 min, 3 times per week on defibrinated bovine blood using an in vitro membrane feeding system . The virus inoculum was prepared from a one pair of hypertrophied salivary glands (with overt SGH symptoms) dissected from G. pallidipes male flies; viral titers in the gland homogenates (in PBS) were estimated by qPCR as described previously [3, 35]. The experimental flies were injected with 2 μl of the virus inoculum estimated to contain ~ 106 virus genome copies per μl.
Identification of core RNAi gene orthologs in Glossina
To determine whether the G. pallidipes genome contains the key RNAi pathway genes, AGO, DCR and Drosha sequences were retrieved from the VectorBase database  using as query sequences the annotated homologous gene sequences of G. morsitans morsitans  (BLASTp; e-value ≤10− 2). To determine the conservation of RNAi in Glossina species, sequences of these core genes for G. fuscipes fuscipes, G. palpalis palpalis, G. austeni and G. brevipalpis were similarly retrieved from VectorBase. The functional domain architecture of the retrieved AGO, DCR and Drosha sequences was analyzed using the ScanProsite tool . Homologous search of the retrieved sequences was performed to determine the relatedness of these sequences with those of D. melanogaster for which RNAi mechanisms and pathways have been demonstrated. Multiple alignments of protein-coding loci of the identified gene sequences were performed in BioEdit . Phylogenetic analysis was performed with MEGA6 using default settings for Maximum-Likelihood (ML) based on the General Time Reversible model with gamma distributed substitution rates with 1000 bootstrap replications .
Analysis of the expression of core RNAi genes in virus-injected G. pallidipes by RT-qPCR
To investigate the impact of GpSGHV on the core RNAi genes in G. pallidipes, two groups of teneral flies (newly-eclosed and non-fed; 50 females and 50 males per group) were injected with either the virus inoculum as described above, or phosphate buffered saline (PBS) as control. Following the injections, four females and four males were sampled within 1 h post injection and at 7, 14 and 21 days post injection (dpi). Total RNA was extracted from individual whole bodies of the sampled flies using Trizol reagent (Invitrogen, Paisley UK) according to the manufacturer’s instructions. Contaminating DNA was removed from the extracted RNA by treating the samples with DNase 1 (Invitrogen, Paisley, UK), after which the concentration of the RNA was quantified using a Nanodrop ND-1000 spectrophotometer (Thermo Fisher Scientific, Wilmington, DE). Complementary DNAs (cDNAs) were synthesised using the SuperScript® III Reverse Transcriptase kit (Invitrogen, Paisley UK) following the manufacturer’s instructions. The iQ SYBR green supermix (Bio-Rad laboratories, Hercules, CA) was used for RT-qPCR analysis. The viral infection was assessed by quantifying the expression of GpSGHV odv-e66, a conserved, late viral gene (highly expressed upon viral genome replication), followed by expression analysis of AGO, DCR and Drosha transcripts, using the PCR cycling conditions: 95 °C for 3 min, followed by forty cycles of 95 °C for 10 s, 60 °C for 1 min, then 95 °C for 1 min and 55 °C for 1 min, using the primers shown in Table 2. The tsetse housekeeping gene β-tubulin was used to normalize gene expression.
Differential expression of RNAi genes in asymptomatic and symptomatic flies
To determine the differential expression of AGO, DCR and Drosha between asymptomatic and symptomatic infected flies, 10-day old F1 progeny flies produced by virus injected mothers were screened under a stereo microscope for the occurence of diagnostic SGH symptoms. Total RNA was extracted from whole bodies of 8 asymptomatic and 8 symptomatic infected flies (4 females and 4 males). The viral infection was estimated by quantification of the GpSGHV odv-e66 transcripts, followed by expression analysis of the three genes as described above. These expression analyses were replicated three times (biological replicates).
Design of dsRNA constructs and prediction of off-targets
The optimal regions on the AGO2 mRNA for the synthesis of dsRNA constructs were determined by siRNA design software (default setting) , which uses three predictive steps; (i) selection of functional siRNA sequence, (ii) selection of siRNA sequence with reduced off-target effects and (iii) elimination of near-perfect matched off-target genes. The identification of off-targets was performed by BLAST (BLASTn; e-value ≤10− 2) search of VectorBase. Following identification of the siRNAs, primers for dsRNAs synthesis were designed to flank the most effective siRNAs based on the above-described steps and a T7 promoter sequence added on each primer (See Additional files 2 and 3).
Synthesis of dsRNAs
To generate dsRNA to knockdown AGO2, total DNA was isolated from G. pallidipes using the Qiagen DNeasy Blood and Tissue kit (QIAGEN Inc., Valencia, CA). The extracted DNA was subsequently used to produce T7 tailed PCR amplicons of AGO2 using primers designed to contain 5′-T7 promotor sequences (See Table 2). These primers allowed dsRNAs transcription using the Hiscribe T7 Quick high yield RNA synthesis kit (New England Biolabs, UK) according to the manufacturer’s instructions. Template DNA was removed from the transcription reaction by DNase treatment, as described in the transcription kit. The synthesized dsRNAs were purified using MEGAclear columns (Ambion, ThermoFisher Scientific, USA) and eluted in 50 μl nuclease free water. The tsetse EP gene, an immune response gene with extensive glutamic acid-proline dipeptide repeats, that has been successfully knocked down in tsetse, [45, 46] was used to assess the efficiency of the knockdown treatment (i.e. by measuring the expression of the tsetse EP gene).
Injections of flies with GpSGHV and dsRNAs
To investigate the impact of AGO2 knockdown on GpSGHV infection, teneral G. pallidipes flies were divided into three groups each consisting of 40 females and 40 males and offered one blood meal. After 48 h, two of three groups were injected with 4 μl of either AGO-2 or tsetse EP dsRNAs (2.5 μg/μl dsRNA in RNase free water) (See the procedure in Additional file 4). The selection of this dsRNA dose was optimized for effective knockdown based on previous bioassays on dsRNA-mediated gene knockdown in tsetse flies [45, 46]. The third group of flies (an additional negative control) was injected with RNase-free water. For the injections, flies were anaesthetized by chilling (~ 5 min) on ice, and subsequently injected in the dorsolateral surface of the thorax. Five days after the dsRNAs/RNase-free water injections, half of the injected flies (20 females and 20 males) from each group were injected with 2 μl of the virus suspension as described above, while the other half were injected with PBS. This time point (i.e. 5 days post dsRNA injection) was selected because successful gene knockdown in tsetse has been shown to occur after ~ 3 dpi , implying that the 5 dpi in our case ensured that the virus was injected after successful knockdown. To monitor the impact of AGO2 knockdown on GpSGHV infection, 3 females and 3 males were sampled from each of the above described treatment groups at 1 h post injection, and at 7, 14 and 21 dpi. The samples were stored at − 20 °C until further analysis as described below.
Analysis of the impacts of AGO2 knockdown on GpSGHV replication
To assess the effect of AGO2 knockdown of GpSGHV replication, total RNA was extracted and cDNA synthesized as described above from the frozen fly samples collected from different time points post dsRNA and virus injection. The efficiency of gene knockdowns was assessed by quantifying (by RT-qPCR) AGO2 and tsetse EP gene transcripts using the qPCR primer sets listed in Table 2. The impact of AGO2 knockdown on GpSGHV infection was assessed by RT-qPCR quantification of mRNA transcripts of the selected conserved GpSGHV genes; the per os infectivity odv-e66 (SGHV005) gene, DNApol (SGHV079) gene involved in DNA replication, a tegument gene (SGHV038) and capsid gene (SGHV091) [59, 60]. Note that a clear correlation between the GpSGHVodv-e66 gene transcripts and the total virus copy numbers has been previously reported , which may demonstrate the impact of AGO2 knockdown on virus replication.
All quantitative RT-qPCR results were representative of at least three independent biological experiments, each with three technical replicates. Statistical differences in the expression of the above described host and viral genes between the different treatments and the controls were performed with RStudio v1.0.143  (R v3.4.0 ) using the packages lattice v0.20–35  and MASS v7.3.47  The obtained data were visualized using the ggplot2 v2.2.1 package  available within the RStudio platform. Data was checked for normality and transformed where necessary using the Box-Cox routine. The data was log transformed where the confidence interval of lambda includes 0 and transformed with (xλ-1)/λ in other cases. T-tests were used for the comparison of RT-qPCR data.
area-wide integrated pest management
Drosophila C virus
double stranded RNA
Flock house virus
Glossina pallidipes salivary gland hypertrophy virus
Helicoverpa armigera single nucleopolyhedrovirus
Insect Pest Control Laboratory
Janus kinase/signal transducers and activators of transcription
RNA induced silencing complex
salivary gland hypertrophy
short interfering RNA
Abd-Alla AM, Vlak JM, Bergoin M, Maruniak JE, Parker A, Burand JP, Jehle JA, Boucias DG. Hytrosaviridae: a proposal for classification and nomenclature of a new insect virus family. Arch Virol. 2009;154(6):909–18.
Abd-Alla AMM, Boucias DG, Bergoin M. Hytrosaviruses: structure and genomic properties. In: Virology I, editor. Asgari SJ, K.N.: caister academic press: Norfolk; 2010. p. 103–21.
Abd-Alla AMM, Cousserans F, Parker AG, Jridi C, Bergoin M, Robinson AS. Quantitative PCR analysis of the salivary gland hypertrophy virus (GpSGHV) in a laboratory colony of Glossina pallidipes. Virus Res. 2009;139(1):48–53.
Kariithi HM, van Oers MM, Vlak JM, Vreysen MJ, Parker AG, Abd-Alla AM. Virology, epidemiology and pathology of Glossina hytrosavirus, and its control prospects in laboratory colonies of the tsetse fly, Glossina pallidipes (Diptera; Glossinidae). Insects. 2013;4(3):287–319.
Kariithi HM, Ahmadi M, Parker AG, Franz G, Ros VID, Haq I, Elashry AM, Vlak JM, Bergoin M, Vreysen MJB, et al. Prevalence and genetic variation of salivary gland hypertrophy virus in wild populations of the tsetse fly Glossina pallidipes from southern and eastern Africa. J Invertebr Pathol. 2013;112(Supplement 1):S123–32.
Abd-Alla AMM, Kariithi HM, Mohamed AH, Lapiz E, Parker AG, Vreysen MJB. Managing hytrosavirus infections in Glossina pallidipes colonies: feeding regime affects the prevalence of salivary gland hypertrophy syndrome. PLoS One. 2013;8(5):e61875.
Alemu T, Kapitano B, Mekonnen S, Aboset G, Kiflom M, Bancha B, Woldeyes G, Bekele K, Feldmann U. Area-wide control of tsetse and trypanosomosis: Ethiopian experience in the southern Rift Valley. In: Vreysen MJB, Robinson AS, Hendrichs J, editors. Area-wide control of insect pests: from research to field implementation. Dordrecht: Springer Netherlands; 2007. p. 325–35.
Boldogh I, Albrecht T, Porter DD. Persistent viral infections. In: Baron S, editor. Med Microbiol. 4th ed: Galveston, University of Texas Medical Branch at Galveston; 1996. (Chapter 46).
Goic B, Saleh M-C. Living with the enemy: viral persistent infections from a friendly viewpoint. Curr Opin Microbiol. 2012;15(4):531–7.
Abd-Alla AMM, Kariithi HM, Parker AG, Robinson AS, Kiflom M, Bergoin M, Vreysen MJB. Dynamics of the salivary gland hypertrophy virus in laboratory colonies of Glossina pallidipes (Diptera: Glossinidae). Virus Res. 2010;150(1–2):103–10.
Kariithi HM, Meki IK, Boucias DG, Abd-Alla AM. Hytrosaviruses: current status and perspective. Curr. Opin. Insect Sci. 2017;22:71–8.
van Rij RP. Virus meets RNAi. Symposium on antiviral applications of RNA interference. EMBO Rep. 2008;9(8):725–9.
Kingsolver MB, Huang Z, Hardy RW. Insect antiviral innate immunity: pathways, effectors, and connections. J Mol Biol. 2013;425(24):4921–36.
Burand JP, Hunter WB. RNAi: future in insect management. J Invertebr Pathol. 2013;112:S68–74.
Bronkhorst AW, van Cleef KWR, Vodovar N, İnce İA, Blanc H, Vlak JM, Saleh M-C, van Rij RP. The DNA virus invertebrate iridescent virus 6 is a target of the Drosophila RNAi machinery. Proc Natl Acad Sci U S A. 2012;109(51):3604–13.
Ding SW. RNA-based antiviral immunity. Nat Rev Immunol. 2010;10(9):632–44.
Lan H, Chen H, Liu Y, Jiang C, Mao Q, Jia D, Chen Q, Wei T. Small interfering RNA pathway modulates initial viral infection in midgut epithelium of insect after ingestion of virus. J Virol. 2016;90(2):917–29.
Skalsky RL, Cullen BR. Viruses, microRNAs, and host interactions. Annu Rev Microbiol. 2010;64:123–41.
Asgari S. MicroRNA functions in insects. Insect Biochem Mol Biol. 2013;43(4):388–97.
Kincaid RP, Sullivan CS. Virus-encoded microRNAs: an overview and a look to the future. PLoS Pathog. 2012;8(12):e1003018.
Nandety RS, Kuo Y-W, Nouri S, Falk BW. Emerging strategies for RNA interference (RNAi) applications in insects. Bioengineered. 2015;6(1):8–19.
Weick E-M, Miska EA. piRNAs: from biogenesis to function. Development. 2014;141(18):3458.
Miesen P, Joosten J, van Rij RP. PIWIs go viral: arbovirus-derived piRNAs in vector mosquitoes. PLoS Pathog. 2016;12(12):e1006017.
Carmell MA, Hannon GJ. RNase III enzymes and the initiation of gene silencing. Nat Struct Mol Biol. 2004;11(3):214–8.
Tomoyasu Y, Miller SC, Tomita S, Schoppmeier M, Grossmann D, Bucher G. Exploring systemic RNA interference in insects: a genome-wide survey for RNAi genes in Tribolium. Genome Biol. 2008;9(1):R10.
Zambon RA, Vakharia VN, Wu LP. RNAi is an antiviral immune response against a dsRNA virus in Drosophila melanogaster. Cell Microbiol. 2006;8(5):880–9.
Kemp C, Mueller S, Goto A, Barbier V, Paro S, Bonnay F, Dostert C, Troxler L, Hetru C, Meignin C, et al. Broad RNA interference-mediated antiviral immunity and virus-specific inducible responses in Drosophila. J Immunol. 2013;190(2):650–8.
Jayachandran B, Hussain M, Asgari S. RNA interference as a cellular defense mechanism against the DNA virus baculovirus. J Virol. 2012;86(24):13729–34.
Bronkhorst AW, van Cleef KWR, Venselaar H, van Rij RP. A dsRNA-binding protein of a complex invertebrate DNA virus suppresses the Drosophila RNAi response. Nucleic Acids Res. 2014;42(19):12237–48.
Sabin LR, Zheng Q, Thekkat P, Yang J, Hannon GJ, Gregory BD, Tudor M, Cherry S. Dicer-2 processes diverse viral RNA species. PLoS One. 2013;8(2):e55458.
Galiana-Arnoux D, Dostert C, Schneemann A, Hoffmann JA, Imler J-L. Essential function in vivo for Dicer-2 in host defense against RNA viruses in Drosophila. Nat Immunol. 2006;7(6):590–7.
Goic B, Vodovar N, Mondotte JA, Monot C, Frangeul L, Blanc H, Gausson V, Vera-Otarola J, Cristofari G, Saleh M-C. RNA-mediated interference and reverse transcription control the persistence of RNA viruses in the insect model Drosophila. Nat Immunol. 2013;14(4):396–403.
Gammon DB, Mello CC. RNA interference-mediated antiviral defense in insects. Curr Opin Insect Sci. 2015;8:111–20.
Bronkhorst AW, van Rij RP: The long and short of antiviral defense: small RNA-based immunity in insects. Curr. Opin. Virol. 2014, 7:19–28.
Boucias DG, Kariithi HM, Bourtzis K, Schneider DI, Kelley K, Miller WJ, Parker AG, Abd-Alla AMM. Transgenerational transmission of the Glossina pallidipes hytrosavirus depends on the presence of a functional symbiome. PLoS One. 2013;8(4):e61150.
Carmell MA, Xuan Z, Zhang MQ, Hannon GJ. The Argonaute family: tentacles that reach into RNAi, developmental control, stem cell maintenance, and tumorigenesis. Genes Dev. 2002;16(21):2733–42.
Lee YS, Nakahara K, Pham JW, Kim K, He Z, Sontheimer EJ, Carthew RW. Distinct roles for Drosophila Dicer-1 and Dicer-2 in the siRNA/miRNA silencing pathways. Cell. 2004;117(1):69–81.
Wang X-H, Aliyari R, Li W-X, Li H-W, Kim K, Carthew R, Atkinson P, Ding S-W. RNA interference directs innate immunity against viruses in adult Drosophila. Science. 2006;312(5772):452–4.
Gu S, Jin L, Huang Y, Zhang F, Kay MA. Slicing-independent RISC activation requires the argonaute PAZ domain. Curr Biol. 2012;22(16):1536–42.
Meister G, Tuschl T. Mechanisms of gene silencing by double-stranded RNA. Nature. 2004;431(7006):343–9.
Parker JS, Barford D. Argonaute: a scaffold for the function of short regulatory RNAs. Trends Biochem Sci. 2006;31(11):622–30.
Keene KM, Foy BD, Sanchez-Vargas I, Beaty BJ, Blair CD, Olson KE. RNA interference acts as a natural antiviral response to O'nyong-nyong virus (Alphavirus; Togaviridae) infection of Anopheles gambiae. Proc Natl Acad Sci U S A. 2004;101(49):17240–5.
Huvenne H, Smagghe G. Mechanisms of dsRNA uptake in insects and potential of RNAi for pest control: a review. J Insect Physiol. 2010;56(3):227–35.
Swevers L, Smagghe G. Use of RNAi for control of insect crop pests. In: Arthropod-plant interactions: novel insights and approaches for IPM. Edited by Smagghe G, Diaz I, vol. 14. Dordrecht: Springer Netherlands; 2012. p. 177–97.
Walshe DP, Lehane SM, Lehane MJ, Haines LR. Prolonged gene knockdown in the tsetse fly Glossina by feeding double stranded RNA. Insect Mol Biol. 2009;18(1):11–9.
Haines LR, Lehane SM, Pearson TW, Lehane MJ. Tsetse EP protein protects the fly midgut from trypanosome establishment. PLoS Pathog. 2010;6(3):e1000793.
van Rij RP, Saleh MC, Berry B, Foo C, Houk A, Antoniewski C, Andino R. The RNA silencing endonuclease Argonaute 2 mediates specific antiviral immunity in Drosophila melanogaster. Genes Dev. 2006;20(21):2985–95.
Costa A, Jan E, Sarnow P, Schneider D. The IMD pathway is involved in antiviral immune responses in Drosophila. PLoS One. 2009;4(10):e7436.
Sabin LR, Hanna SL, Cherry S. Innate antiviral immunity in Drosophila. Curr Opin Immunol. 2010;22(1):4–9.
Fullaondo A, Lee SY. Regulation of Drosophila-virus interaction. Dev Comp Immunol. 2012;36(2):262–6.
Rechavi O. “Guest list” or “black list”? Heritable small RNAs as immunogenic memories. Trends Cell Biol. 2014;24(4):212–20.
Feldmann U. Guidelines for the rearing of tsetse flies using the membrane feeding technique. In: Ochieng'-Odero JPR, editor. Techniques of insect rearing for the development of integrated pest and vector management strategies. Nairobi: ICIPE Science Press; 1994. p. 449–71.
Giraldo-Calderón GI, Emrich SJ, MacCallum RM, Maslen G, Dialynas E, Topalis P, Ho N, Gesing S, the VectorBase C, Madey G, et al. VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases. Nucleic Acids Res. 2015;43(Database issue):D707–13.
International Glossina Genome Initiative. Genome sequence of the tsetse fly (Glossina morsitans): vector of African trypanosomiasis. Science. 2014;344(6182):380–6.
De Castro E, Sigrist CJA, Gattiker A, Bulliard V, Langendijk-Genevaux PS, Gasteiger E, Bairoch A, Hulo N. ScanProsite: detection of PROSITE signature matches and ProRule-associated functional and structural residues in proteins. Nucleic Acids Res. 2006;34(Web Server):W362–5.
Hall T. BioEdit: an important software for molecular biology. GERF Bulletin of Biosciences. 2011;2:30–61.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.
Naito Y, Ui-Tei K. siRNA design software for a target gene-specific RNA interference. Front Genet. 2012;3:102.
Abd-Alla AMM, Cousserans F, Parker AG, Jehle JA, Parker NJ, Vlak JM, Robinson AS, Bergoin M. Genome analysis of a Glossina pallidipes salivary gland hypertrophy virus reveals a novel, large, double-stranded circular DNA virus. J Virol. 2008;82(9):4595–611.
Abd-Alla AMM, Kariithi HM, Cousserans F, Parker NJ, İnce İA, Scully ED, Boeren S, Geib SM, Mekonnen S, Vlak JM, et al. Comprehensive annotation of Glossina pallidipes salivary gland hypertrophy virus from Ethiopian tsetse flies: a proteogenomics approach. J Gen Virol. 2016;97(4):1010–31.
RStudio T. RStudio: integrated development environment for R. RStudio, Inc. Boston. In: MA; 2015.
R-Core T. R: a language and environment for statistical computing. Vienna, Austria: R Foundation for statistical computing; 2015.
Sarkar D. Lattice: multivariate data visualization with R. New York: Springer; 2008.
Venables WN, Ripley BD. Modern applied statistics with S. 4th ed. New York: Springer; 2002.
Wickham H. Ggplot2: elegant graphics for data analysis. New York: Springer; 2009.
Caljon G, Broos K, De Goeyse I, De Ridder K, Sternberg JM, Coosemans M, De Baetselier P, Guisez Y, Van Den Abbeele J. Identification of a functional Antigen5-related allergen in the saliva of a blood feeding insect, the tsetse fly. Insect Biochem Mol Biol. 2009;39(5–6):332–41.
The authors acknowledge Prof. Serap Aksoy of Yale University for insightful discussions during the conception of the study and Hasim Mohammed, Carmen Marin and Henry Adun of the Joint FAO/IAEA Insect Pest Control Laboratory (IPCL), Seibersdorf, Austria for technical assistance.
This work was financially supported by the Joint FAO/IAEA Division of Nuclear Techniques in Food and Agriculture, Vienna, Austria (CRP #D42015) and Wageningen University (WU). Both FAO/IAEA Division and WU participated in the design of the study, provided the resources for running of the experiments and writing of the manuscript.
Availability of data and materials
Materials described in the manuscript, including all relevant raw data, are deposited in this link https://dataverse.harvard.edu/dataset.xhtml?persistentId=doi:10.7910/DVN/ZN20Q6.
About this supplement
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Regression parameters applied in the statistical analysis. For each regression line (column 2) found in the figures (column 1), the confidence interval of lambda values (column 3) obtained based on the Box-Cox routine to determine the method of transformation (column 4) of the data, the intercept values (column 5), the t value (column 6) the P values (column 7) and the degrees of freedom (D.F) in column 8. (DOCX 23 kb)
Result of AGO2 dsRNA design and off-target prediction. A). AGO2 mRNA sequence; the primers flanking the knockdown sequence are highlighted in cyan, predicted siRNAs in grey and qPCR primers in yellow. (B) Shows a graphical view of effective siRNAs candidates. (TIF 9175 kb)
Result of tsetse EP dsRNA design and off-target prediction. A). Tsetse EP mRNA sequence; the primers flanking knockdown sequence are highlighted in cyan, predicted siRNAs in grey and qPCR primers in yellow. (B) Shows a graphical view of effective siRNAs candidates. (TIF 8491 kb)
Step by step procedure to determine the effect of Knock-down of AGO2 on virus infection in G. pallidipes. Teneral G. pallidipes flies were collected and offered one blood meal. 48 h later the flies were divided onto three groups and injected with 4 μl AGO2 dsRNAs, tsetse EP dsRNAs or RNase-free water. Five days later, each of the 3 groups was divided into two and each injected with either 2 μl of the virus suspension or with PBS. Three females and 3 males were sampled from each of the described treatments at 1 h post injection, and at 7, 14 and 21 days post injection to determine the effect of AGO2 knockdown on virus infection. (TIF 1030 kb)
About this article
Cite this article
Meki, I.K., Kariithi, H.M., Parker, A.G. et al. RNA interference-based antiviral immune response against the salivary gland hypertrophy virus in Glossina pallidipes. BMC Microbiol 18, 170 (2018). https://doi.org/10.1186/s12866-018-1298-1
- Sterile insect technique
- Symptomatic and asymptomatic infection
- Covert infections