A satellite repeat-derived piRNA controls embryonic development of Aedes.

Tandem repeat elements such as the diverse class of satellite repeats occupy large parts of eukaryotic chromosomes, mostly at centromeric, pericentromeric, telomeric and subtelomeric regions1. However, some elements are located in euchromatic regions throughout the genome and have been hypothesized to regulate gene expression in cis by modulating local chromatin structure, or in trans via transcripts derived from the repeats2-4. Here we show that a satellite repeat in the mosquito Aedes aegypti promotes sequence-specific gene silencing via the expression of two PIWI-interacting RNAs (piRNAs). Whereas satellite repeats and piRNA sequences generally evolve extremely quickly5-7, this locus was conserved for approximately 200 million years, suggesting that it has a central function in mosquito biology. piRNA production commenced shortly after egg laying, and inactivation of the more abundant piRNA resulted in failure to degrade maternally deposited transcripts in the zygote and developmental arrest. Our results reveal a mechanism by which satellite repeats regulate global gene expression in trans via piRNA-mediated gene silencing that is essential for embryonic development.

Pan-flavivirus analysis reveals sfRNA-independent, 3' UTR-biased siRNA production from an insect-specific flavivirus.

RNA interference (RNAi) plays an essential role in mosquito antiviral immunity, but it is not known whether viral small interfering RNA (siRNA) profiles differ between mosquito-borne and mosquito-specific viruses. A pan-Orthoflavivirus analysis in Aedes albopictus cells revealed that viral siRNAs were evenly distributed across the viral genome of most representatives of the Flavivirus genus. In contrast, siRNA production was biased toward the 3' untranslated region (UTR) of the genomes of classical insect-specific flaviviruses (cISF), which was most pronounced for Kamiti River virus (KRV), a virus with a unique, 1.2 kb long 3' UTR. KRV-derived siRNAs were produced in high quantities and almost exclusively mapped to the 3' UTR. We mapped the 5' end of KRV subgenomic flavivirus RNAs (sfRNAs), products of the 5'-3' exoribonuclease XRN1/Pacman stalling on secondary RNA structures in the 3' UTR of the viral genome. We found that KRV produces high copy numbers of a long, 1,017 nt sfRNA1 and a short, 421 nt sfRNA2, corresponding to two predicted XRN1-resistant elements. Expression of both sfRNA1 and sfRNA2 was reduced in Pacman-deficient Aedes albopictus cells; however, this did not correlate with a shift in viral siRNA profiles. We suggest that cISFs, particularly KRV, developed a unique mechanism to produce high amounts of siRNAs as a decoy for the antiviral RNAi response in an sfRNA-independent manner.IMPORTANCEThe Flavivirus genus contains diverse mosquito viruses ranging from insect-specific viruses circulating exclusively in mosquito populations to mosquito-borne viruses that cause disease in humans and animals. Studying the mechanisms of virus replication and antiviral immunity in mosquitoes is important to understand arbovirus transmission and may inform the development of disease control strategies. In insects, RNA interference (RNAi) provides broad antiviral activity and constitutes a major immune response against viruses. Comparing diverse members of the Flavivirus genus, we found that all flaviviruses are targeted by RNAi. However, the insect-specific Kamiti River virus was unique in that small interfering RNAs are highly skewed toward its uniquely long 3' untranslated region. These results suggest that mosquito-specific viruses have evolved unique mechanisms for genome replication and immune evasion.

Arbovirus-vector protein interactomics identifies Loquacious as a co-factor for dengue virus replication in Aedes mosquitoes.

Efficient virus replication in Aedes vector mosquitoes is essential for the transmission of arboviral diseases such as dengue virus (DENV) in human populations. Like in vertebrates, virus-host protein-protein interactions are essential for viral replication and immune evasion in the mosquito vector. Here, 79 mosquito host proteins interacting with DENV non-structural proteins NS1 and NS5 were identified by label-free mass spectrometry, followed by a functional screening. We confirmed interactions with host factors previously observed in mammals, such as the oligosaccharyltransferase complex, and we identified protein-protein interactions that seem to be specific for mosquitoes. Among the interactors, the double-stranded RNA (dsRNA) binding protein Loquacious (Loqs), an RNA interference (RNAi) cofactor, was found to be essential for efficient replication of DENV and Zika virus (ZIKV) in mosquito cells. Loqs did not affect viral RNA stability or translation of a DENV replicon and its proviral activity was independent of its RNAi regulatory activity. Interestingly, Loqs colocalized with DENV dsRNA replication intermediates in infected cells and directly interacted with high affinity with DENV RNA in the 3' untranslated region in vitro (KD = 48-62 nM). Our study provides an interactome for DENV NS1 and NS5 and identifies Loqs as a key proviral host factor in mosquitoes. We propose that DENV hijacks a factor of the RNAi mechanism for replication of its own RNA.

The impact of BNT162b2 mRNA vaccine on adaptive and innate immune responses.

The mRNA-based BNT162b2 protects against severe disease and mortality caused by SARS-CoV-2 via induction of specific antibody and T-cell responses. Much less is known about its broad effects on immune responses against other pathogens. Here, we investigated the adaptive immune responses induced by BNT162b2 vaccination against various SARS-CoV-2 variants and its effects on the responsiveness of immune cells upon stimulation with heterologous stimuli. BNT162b2 vaccination induced effective humoral and cellular immunity against SARS-CoV-2 that started to wane after six months. We also observed long-term transcriptional changes in immune cells after vaccination. Additionally, vaccination with BNT162b2 modulated innate immune responses as measured by inflammatory cytokine production after stimulation - higher IL-1/IL-6 release and decreased IFN-α production. Altogether, these data expand our knowledge regarding the overall immunological effects of this new class of vaccines and underline the need for additional studies to elucidate their effects on both innate and adaptive immune responses.

A piRNA-lncRNA regulatory network initiates responder and trailer piRNA formation during mosquito embryonic development.

PIWI-interacting (pi)RNAs are small silencing RNAs that are crucial for the defense against transposable elements in germline tissues of animals. In Aedes aegypti mosquitoes, the piRNA pathway also contributes to gene regulation in somatic tissues, illustrating additional roles for piRNAs and PIWI proteins besides transposon repression. Here, we identify a highly abundant endogenous piRNA (propiR1) that associates with both Piwi4 and Piwi5. PropiR1-mediated target silencing requires base-pairing in the seed region with supplemental base-pairing at the piRNA 3' end. Yet, propiR1 represses a limited set of targets, among which is the lncRNA AAEL027353 (lnc027353). Slicing of lnc027353 initiates production of responder and trailer piRNAs from the cleavage fragment. Expression of propiR1 commences early during embryonic development and mediates degradation of maternally provided lnc027353 Both propiR1 and its lncRNA target are conserved in the closely related Aedes albopictus mosquito, underscoring the importance of this regulatory network for mosquito development.

Endogenous piRNA-guided slicing triggers responder and trailer piRNA production from viral RNA in Aedes aegypti mosquitoes.

In the germline of animals, PIWI interacting (pi)RNAs protect the genome against the detrimental effects of transposon mobilization. In Drosophila, piRNA-mediated cleavage of transposon RNA triggers the production of responder piRNAs via ping-pong amplification. Responder piRNA 3' end formation by the nuclease Zucchini is coupled to the production of downstream trailer piRNAs, expanding the repertoire of transposon piRNA sequences. In Aedes aegypti mosquitoes, piRNAs are generated from viral RNA, yet, it is unknown how viral piRNA 3' ends are formed and whether viral RNA cleavage gives rise to trailer piRNA production. Here we report that in Ae. aegypti, virus- and transposon-derived piRNAs have sharp 3' ends, and are biased for downstream uridine residues, features reminiscent of Zucchini cleavage of precursor piRNAs in Drosophila. We designed a reporter system to study viral piRNA 3' end formation and found that targeting viral RNA by abundant endogenous piRNAs triggers the production of responder and trailer piRNAs. Using this reporter, we identified the Ae. aegypti orthologs of Zucchini and Nibbler, two nucleases involved in piRNA 3' end formation. Our results furthermore suggest that autonomous piRNA production from viral RNA can be triggered and expanded by an initial cleavage event guided by genome-encoded piRNAs.

The histone methyltransferase G9a regulates tolerance to oxidative stress-induced energy consumption.

Stress responses are crucial processes that require activation of genetic programs that protect from the stressor. Stress responses are also energy consuming and can thus be deleterious to the organism. The mechanisms coordinating energy consumption during stress response in multicellular organisms are not well understood. Here, we show that loss of the epigenetic regulator G9a in Drosophila causes a shift in the transcriptional and metabolic responses to oxidative stress (OS) that leads to decreased survival time upon feeding the xenobiotic paraquat. During OS exposure, G9a mutants show overactivation of stress response genes, rapid depletion of glycogen, and inability to access lipid energy stores. The OS survival deficiency of G9a mutants can be rescued by a high-sugar diet. Control flies also show improved OS survival when fed a high-sugar diet, suggesting that energy availability is generally a limiting factor for OS tolerance. Directly limiting access to glycogen stores by knocking down glycogen phosphorylase recapitulates the OS-induced survival defects of G9a mutants. We propose that G9a mutants are sensitive to stress because they experience a net reduction in available energy due to (1) rapid glycogen use, (2) an inability to access lipid energy stores, and (3) an overinduced transcriptional response to stress that further exacerbates energy demands. This suggests that G9a acts as a critical regulatory hub between the transcriptional and metabolic responses to OS. Our findings, together with recent studies that established a role for G9a in hypoxia resistance in cancer cell lines, suggest that G9a is of wide importance in controlling the cellular and organismal response to multiple types of stress.

A DNA virus-encoded immune antagonist fully masks the potent antiviral activity of RNAi in Drosophila.

Coevolution of viruses and their hosts may lead to viral strategies to avoid, evade, or suppress antiviral immunity. An example is antiviral RNA interference (RNAi) in insects: the host RNAi machinery processes viral double-stranded RNA into small interfering RNAs (siRNAs) to suppress viral replication, whereas insect viruses encode suppressors of RNAi, many of which inhibit viral small interfering RNA (vsiRNA) production. Yet, many studies have analyzed viral RNAi suppressors in heterologous systems, due to the lack of experimental systems to manipulate the viral genome of interest, raising questions about in vivo functions of RNAi suppressors. To address this caveat, we generated an RNAi suppressor-defective mutant of invertebrate iridescent virus 6 (IIV6), a large DNA virus in which we previously identified the 340R protein as a suppressor of RNAi. Loss of 340R did not affect vsiRNA production, indicating that 340R binds siRNA duplexes to prevent RNA-induced silencing complex assembly. Indeed, vsiRNAs were not efficiently loaded into Argonaute 2 during wild-type IIV6 infection. Moreover, IIV6 induced a limited set of mature microRNAs in a 340R-dependent manner, most notably miR-305-3p, which we attribute to stabilization of the miR-305-5p:3p duplex by 340R. The IIV6 340R deletion mutant did not have a replication defect in cells, but was strongly attenuated in adult Drosophila This in vivo replication defect was completely rescued in RNAi mutant flies, indicating that 340R is a bona fide RNAi suppressor, the absence of which uncovers a potent antiviral immune response that suppresses virus accumulation ∼100-fold. Together, our work indicates that viral RNAi suppressors may completely mask antiviral immunity.

Cationic Geminoid Peptide Amphiphiles Inhibit DENV2 Protease, Furin, and Viral Replication.

Dengue is an important arboviral infectious disease for which there is currently no specific cure. We report gemini-like (geminoid) alkylated amphiphilic peptides containing lysines in combination with glycines or alanines (C15H31C(O)-Lys-(Gly or Ala)nLys-NHC16H33, shorthand notation C16-KXnK-C16 with X = A or G, and n = 0-2). The representatives with 1 or 2 Ala inhibit dengue protease and human furin, two serine proteases involved in dengue virus infection that have peptides with cationic amino acids as their preferred substrates, with IC50 values in the lower µM range. The geminoid C16-KAK-C16 combined inhibition of DENV2 protease (IC50 2.3 µM) with efficacy against replication of wildtype DENV2 in LLC-MK2 cells (EC50 4.1 µM) and an absence of toxicity. We conclude that the lysine-based geminoids have activity against dengue virus infection, which is based on their inhibition of the proteases involved in viral replication and are therefore promising leads to further developing antiviral therapeutics, not limited to dengue.

Desialylation of platelets induced by Von Willebrand Factor is a novel mechanism of platelet clearance in dengue.

Thrombocytopenia and platelet dysfunction are commonly observed in patients with dengue virus (DENV) infection and may contribute to complications such as bleeding and plasma leakage. The etiology of dengue-associated thrombocytopenia is multifactorial and includes increased platelet clearance. The binding of the coagulation protein von Willebrand factor (VWF) to the platelet membrane and removal of sialic acid (desialylation) are two well-known mechanisms of platelet clearance, but whether these conditions also contribute to thrombocytopenia in dengue infection is unknown. In two observational cohort studies in Bandung and Jepara, Indonesia, we show that adult patients with dengue not only had higher plasma concentrations of plasma VWF antigen and active VWF, but that circulating platelets had also bound more VWF to their membrane. The amount of platelet-VWF binding correlated well with platelet count. Furthermore, sialic acid levels in dengue patients were significantly reduced as assessed by the binding of Sambucus nigra lectin (SNA) and Maackia amurensis lectin II (MAL-II) to platelets. Sialic acid on the platelet membrane is neuraminidase-labile, but dengue virus has no known neuraminidase activity. Indeed, no detectable activity of neuraminidase was present in plasma of dengue patients and no desialylation was found of plasma transferrin. Platelet sialylation was also not altered by in vitro exposure of platelets to DENV nonstructural protein 1 or cultured DENV. In contrast, induction of binding of VWF to glycoprotein 1b on platelets using the VWF-activating protein ristocetin resulted in the removal of platelet sialic acid by translocation of platelet neuraminidase to the platelet surface. The neuraminidase inhibitor oseltamivir reduced VWF-induced platelet desialylation. Our data demonstrate that excessive binding of VWF to platelets in dengue results in neuraminidase-mediated platelet desialylation and platelet clearance. Oseltamivir might be a novel treatment option for severe thrombocytopenia in dengue infection.

Population genomics in the arboviral vector Aedes aegypti reveals the genomic architecture and evolution of endogenous viral elements.

Horizontal gene transfer from viruses to eukaryotic cells is a pervasive phenomenon. Somatic viral integrations are linked to persistent viral infection whereas integrations into germline cells are maintained in host genomes by vertical transmission and may be co-opted for host functions. In the arboviral vector Aedes aegypti, an endogenous viral element from a nonretroviral RNA virus (nrEVE) was shown to produce PIWI-interacting RNAs (piRNAs) to limit infection with a cognate virus. Thus, nrEVEs may constitute a heritable, sequence-specific mechanism for antiviral immunity, analogous to piRNA-mediated silencing of transposable elements. Here, we combine population genomics and evolutionary approaches to analyse the genomic architecture of nrEVEs in A. aegypti. We conducted a genome-wide screen for adaptive nrEVEs and searched for novel population-specific nrEVEs in the genomes of 80 individual wild-caught mosquitoes from five geographical populations. We show a dynamic landscape of nrEVEs in mosquito genomes and identified five novel nrEVEs derived from two currently circulating viruses, providing evidence of the environmental-dependent modification of a piRNA cluster. Overall, our results show that virus endogenization events are complex with only a few nrEVEs contributing to adaptive evolution in A. aegypti.

The Tudor protein Veneno assembles the ping-pong amplification complex that produces viral piRNAs in Aedes mosquitoes.

PIWI-interacting RNAs (piRNAs) comprise a class of small RNAs best known for suppressing transposable elements in germline tissues. The vector mosquito Aedes aegypti encodes seven PIWI genes, four of which are somatically expressed. This somatic piRNA pathway generates piRNAs from viral RNA during infection with cytoplasmic RNA viruses through ping-pong amplification by the PIWI proteins Ago3 and Piwi5. Yet, additional insights into the molecular mechanisms mediating non-canonical piRNA production are lacking. TUDOR-domain containing (Tudor) proteins facilitate piRNA biogenesis in Drosophila melanogaster and other model organisms. We thus hypothesized that Tudor proteins are required for viral piRNA production and performed a knockdown screen targeting all A. aegypti Tudor genes. Knockdown of the Tudor genes AAEL012437, Vreteno, Yb, SMN and AAEL008101-RB resulted in significantly reduced viral piRNA levels, with AAEL012437-depletion having the strongest effect. This protein, which we named Veneno, associates directly with Ago3 in an sDMA-dependent manner and localizes in cytoplasmic foci reminiscent of piRNA processing granules of Drosophila. Veneno-interactome analyses reveal a network of co-factors including the orthologs of the Drosophila piRNA pathway components Vasa and Yb, which in turn interacts with Piwi5. We propose that Veneno assembles a multi-protein complex for ping-pong dependent piRNA production from viral RNA.

Agua Salud alphavirus defines a novel lineage of insect-specific alphaviruses discovered in the New World.

The genus Alphavirus harbours mostly insect-transmitted viruses that cause severe disease in humans, livestock and wildlife. Thus far, only three alphaviruses with a host range restricted to insects have been found in mosquitoes from the Old World, namely Eilat virus (EILV), Taï Forest alphavirus (TALV) and Mwinilunga alphavirus (MWAV). In this study, we found a novel alphavirus in one Culex declarator mosquito sampled in Panama. The virus was isolated in C6/36 mosquito cells, and full genome sequencing revealed an 11 468 nt long genome with maximum pairwise nucleotide identity of 62.7 % to Sindbis virus. Phylogenetic analyses placed the virus as a solitary deep rooting lineage in a basal relationship to the Western equine encephalitis antigenic complex and to the clade comprising EILV, TALV and MWAV, indicating the detection of a novel alphavirus, tentatively named Agua Salud alphavirus (ASALV). No growth of ASALV was detected in vertebrate cell lines, including cell lines derived from ectothermic animals, and replication of ASALV was strongly impaired above 31 °C, suggesting that ASALV represents the first insect-restricted alphavirus of the New World.

Induction and Suppression of NF-κB Signalling by a DNA Virus of Drosophila.

Interactions between the insect immune system and RNA viruses have been extensively studied in Drosophila, in which RNA interference, NF-κB, and JAK-STAT pathways underlie antiviral immunity. In response to RNA interference, insect viruses have convergently evolved suppressors of this pathway that act by diverse mechanisms to permit viral replication. However, interactions between the insect immune system and DNA viruses have received less attention, primarily because few Drosophila-infecting DNA virus isolates are available. In this study, we used a recently isolated DNA virus of Drosophila melanogaster, Kallithea virus (KV; family Nudiviridae), to probe known antiviral immune responses and virus evasion tactics in the context of DNA virus infection. We found that fly mutants for RNA interference and immune deficiency (Imd), but not Toll, pathways are more susceptible to Kallithea virus infection. We identified the Kallithea virus-encoded protein gp83 as a potent inhibitor of Toll signalling, suggesting that Toll mediates antiviral defense against Kallithea virus infection but that it is suppressed by the virus. We found that Kallithea virus gp83 inhibits Toll signalling through the regulation of NF-κB transcription factors. Furthermore, we found that gp83 of the closely related Drosophila innubila nudivirus (DiNV) suppresses D. melanogaster Toll signalling, suggesting an evolutionarily conserved function of Toll in defense against DNA viruses. Together, these results provide a broad description of known antiviral pathways in the context of DNA virus infection and identify the first Toll pathway inhibitor in a Drosophila virus, extending the known diversity of insect virus-encoded immune inhibitors.IMPORTANCE Coevolution of multicellular organisms and their natural viruses may lead to an intricate relationship in which host survival requires effective immunity and virus survival depends on evasion of such responses. Insect antiviral immunity and reciprocal virus immunosuppression tactics have been well studied in Drosophila melanogaster, primarily during RNA, but not DNA, virus infection. Therefore, we describe interactions between a recently isolated Drosophila DNA virus (Kallithea virus [KV]) and immune processes known to control RNA viruses, such as RNA interference (RNAi) and Imd pathways. We found that KV suppresses the Toll pathway and identified gp83 as a KV-encoded protein that underlies this suppression. This immunosuppressive ability is conserved in another nudivirus, suggesting that the Toll pathway has conserved antiviral activity against DNA nudiviruses, which have evolved suppressors in response. Together, these results indicate that DNA viruses induce and suppress NF-κB responses, and they advance the application of KV as a model to study insect immunity.

Viral suppressors of RNAi employ a rapid screening mode to discriminate viral RNA from cellular small RNA.

RNA interference (RNAi) is an indispensable mechanism for antiviral defense in insects, including mosquitoes that transmit human diseases. To escape this antiviral defense system, viruses encode suppressors of RNAi that prevent elimination of viral RNAs, and thus ensure efficient virus accumulation. Although the first animal Viral Suppressor of RNAi (VSR) was identified more than a decade ago, the molecular basis of RNAi suppression by these viral proteins remains unclear. Here, we developed a single-molecule fluorescence assay to investigate how VSRs inhibit the recognition of viral RNAs by Dcr-2, a key endoribonuclease enzyme in the RNAi pathway. Using VSRs from three insect RNA viruses (Culex Y virus, Drosophila X virus and Drosophila C virus), we reveal bimodal physical interactions between RNA molecules and VSRs. During initial interactions, these VSRs rapidly discriminate short RNA substrates from long dsRNA. VSRs engage nearly irreversible binding with long dsRNAs, thereby shielding it from recognition by Dcr-2. We propose that the length-dependent switch from rapid screening to irreversible binding reflects the main mechanism by which VSRs distinguish viral dsRNA from cellular RNA species such as microRNAs.

Histone-derived piRNA biogenesis depends on the ping-pong partners Piwi5 and Ago3 in Aedes aegypti.

The piRNA pathway is of key importance in controlling transposable elements in most animal species. In the vector mosquito Aedes aegypti, the presence of eight PIWI proteins and the accumulation of viral piRNAs upon arbovirus infection suggest additional functions of the piRNA pathway beyond genome defense. To better understand the regulatory potential of this pathway, we analyzed in detail host-derived piRNAs in A. aegypti Aag2 cells. We show that a large repertoire of protein-coding genes and non-retroviral integrated RNA virus elements are processed into genic piRNAs by different combinations of PIWI proteins. Among these, we identify a class of genes that produces piRNAs from coding sequences in an Ago3- and Piwi5-dependent fashion. We demonstrate that the replication-dependent histone gene family is a genic source of ping-pong dependent piRNAs and that histone-derived piRNAs are dynamically expressed throughout the cell cycle, suggesting a role for the piRNA pathway in the regulation of histone gene expression. Moreover, our results establish the Aag2 cell line as an accessible experimental model to study gene-derived piRNAs.

PIWIs Go Viral: Arbovirus-Derived piRNAs in Vector Mosquitoes.

Vector mosquitoes are responsible for transmission of the majority of arthropod-borne (arbo-) viruses. Virus replication in these vectors needs to be sufficiently high to permit efficient virus transfer to vertebrate hosts. The mosquito immune response therefore is a key determinant for arbovirus transmission. Mosquito antiviral immunity is primarily mediated by the small interfering RNA pathway. Besides this well-established antiviral machinery, the PIWI-interacting RNA (piRNA) pathway processes viral RNA into piRNAs. In recent years, significant progress has been made in characterizing the biogenesis and function of these viral piRNAs. In this review, we discuss these developments, identify knowledge gaps, and suggest directions for future research.

Comparative genomics shows that viral integrations are abundant and express piRNAs in the arboviral vectors Aedes aegypti and Aedes albopictus.

BACKGROUND: Arthropod-borne viruses (arboviruses) transmitted by mosquito vectors cause many important emerging or resurging infectious diseases in humans including dengue, chikungunya and Zika. Understanding the co-evolutionary processes among viruses and vectors is essential for the development of novel transmission-blocking strategies. Episomal viral DNA fragments are produced from arboviral RNA upon infection of mosquito cells and adults. Additionally, sequences from insect-specific viruses and arboviruses have been found integrated into mosquito genomes. RESULTS: We used a bioinformatic approach to analyse the presence, abundance, distribution, and transcriptional activity of integrations from 425 non-retroviral viruses, including 133 arboviruses, across the presently available 22 mosquito genome sequences. Large differences in abundance and types of viral integrations were observed in mosquito species from the same region. Viral integrations are unexpectedly abundant in the arboviral vector species Aedes aegypti and Ae. albopictus, in which they are approximately ~10-fold more abundant than in other mosquito species analysed. Additionally, viral integrations are enriched in piRNA clusters of both the Ae. aegypti and Ae. albopictus genomes and, accordingly, they express piRNAs, but not siRNAs. CONCLUSIONS: Differences in the number of viral integrations in the genomes of mosquito species from the same geographic area support the conclusion that integrations of viral sequences is not dependent on viral exposure, but that lineage-specific interactions exist. Viral integrations are abundant in Ae. aegypti and Ae. albopictus, and represent a thus far underappreciated component of their genomes. Additionally, the genome locations of viral integrations and their production of piRNAs indicate a functional link between viral integrations and the piRNA pathway. These results greatly expand the breadth and complexity of small RNA-mediated regulation and suggest a role for viral integrations in antiviral defense in these two mosquito species.

Unity in defence: honeybee workers exhibit conserved molecular responses to diverse pathogens.

BACKGROUND: Organisms typically face infection by diverse pathogens, and hosts are thought to have developed specific responses to each type of pathogen they encounter. The advent of transcriptomics now makes it possible to test this hypothesis and compare host gene expression responses to multiple pathogens at a genome-wide scale. Here, we performed a meta-analysis of multiple published and new transcriptomes using a newly developed bioinformatics approach that filters genes based on their expression profile across datasets. Thereby, we identified common and unique molecular responses of a model host species, the honey bee (Apis mellifera), to its major pathogens and parasites: the Microsporidia Nosema apis and Nosema ceranae, RNA viruses, and the ectoparasitic mite Varroa destructor, which transmits viruses. RESULTS: We identified a common suite of genes and conserved molecular pathways that respond to all investigated pathogens, a result that suggests a commonality in response mechanisms to diverse pathogens. We found that genes differentially expressed after infection exhibit a higher evolutionary rate than non-differentially expressed genes. Using our new bioinformatics approach, we unveiled additional pathogen-specific responses of honey bees; we found that apoptosis appeared to be an important response following microsporidian infection, while genes from the immune signalling pathways, Toll and Imd, were differentially expressed after Varroa/virus infection. Finally, we applied our bioinformatics approach and generated a gene co-expression network to identify highly connected (hub) genes that may represent important mediators and regulators of anti-pathogen responses. CONCLUSIONS: Our meta-analysis generated a comprehensive overview of the host metabolic and other biological processes that mediate interactions between insects and their pathogens. We identified key host genes and pathways that respond to phylogenetically diverse pathogens, representing an important source for future functional studies as well as offering new routes to identify or generate pathogen resilient honey bee stocks. The statistical and bioinformatics approaches that were developed for this study are broadly applicable to synthesize information across transcriptomic datasets. These approaches will likely have utility in addressing a variety of biological questions.