Comprehensive Quantitative Proteome Analysis of Aedes aegypti Identifies Proteins and Pathways Involved in Wolbachia pipientis and Zika Virus Interference Phenomenon

Abstract

Zika virus (ZIKV) is a global public health emergency due to its association with microcephaly, Guillain-Barré syndrome, neuropathy, and myelitis in children and adults. A total of 87 countries have had evidence of autochthonous mosquito-borne transmission of ZIKV, distributed across four continents, and no antivirus therapy or vaccines are available. Therefore, several strategies have been developed to target the main mosquito vector, Aedes aegypti, to reduce the burden of different arboviruses. Among such strategies, the use of the maternally-inherited endosymbiont Wolbachia pipientis has been applied successfully to reduce virus susceptibility and decrease transmission. However, the mechanisms by which Wolbachia orchestrate resistance to ZIKV infection remain to be elucidated. In this study, we apply isobaric labeling quantitative mass spectrometry (MS)-based proteomics to quantify proteins and identify pathways altered during ZIKV infection; Wolbachia infection; co-infection with Wolbachia/ZIKV in the A. aegypti heads and salivary glands. We show that Wolbachia regulates proteins involved in reactive oxygen species production, regulates humoral immune response, and antioxidant production. The reduction of ZIKV polyprotein in the presence of Wolbachia in mosquitoes was determined by MS and corroborates the idea that Wolbachia helps to block ZIKV infections in A. aegypti. The present study offers a rich resource of data that may help to elucidate mechanisms by which Wolbachia orchestrate resistance to ZIKV infection in A. aegypti, and represents a step further on the development of new targeted methods to detect and quantify ZIKV and Wolbachia directly in complex tissues.

Highlights

• -

  The abundance of ZIKV polyprotein is reduced in the presence of Wolbachia.

• -

  Shotgun proteomics quantifies ZIKV and Wolbachia proteins directly in tissues.

• -

  Wolbachia regulates proteins involved in ROS production.

• -

  Wolbachia regulates humoral immune response and antioxidant production.

• -

  Metabolism and detoxification processes were associated with mono infections.

Introduction

Zika virus (ZIKV) is a single-stranded RNA virus that belongs to the Flavivirus genus and Flaviviridae family. Originally discovered in a primate from the Zika forest in Uganda in 1947, remained unimportant for public health until the beginning of the 20th century (Dick et al., 1952). In 2014, ZIKV emerged in the Pacific islands and a few years later invaded the Americas, leading the World Health Organization to claim a global public health emergency due to its association with microcephaly in newborns and Guillain-Barré syndrome, neuropathy, and myelitis in children and adults (Campos et al., 2015; Barreto et al., 2016). As of May 2017, around 75% of inhabitants of the Yap Islands were infected with ZIKV, but no microcephaly was recorded. Whilst in Brazil more than 220,000 cases were recorded and more than 2,600 newborns had birth defects like microcephaly and other neuropathies (Baud et al., 2017).

The mosquito Aedes aegypti is the primary vector of ZIKV worldwide (Chouin-Carneiro et al., 2016; Ferreira-de-Brito et al., 2016). To date, no antivirus therapy or vaccines are available to mitigate ZIKV transmission, which increases the urgency for effective methods targeting the mosquito vector to reduce the burden of arboviruses as ZIKV. Among such strategies, the use of the maternally inherited endosymbiont Wolbachia pipientis has gained attention recently. The bacterium Wolbachia is naturally found in around 60% of arthropod species. Although Wolbachia does not naturally infect A. aegypti, some strains were transinfected by microinjection into this mosquito species with the further discovery of their pathogen-blocking properties (Moreira et al., 2009; Aliota et al., 2016; Dutra et al., 2016). The rationale underlying a mass release of Wolbachia infected mosquitoes is to replace the local A. aegypti population. Usually, it has high vector competence to arboviruses, like ZIKV, while the Wolbachia transinfected mosquitoes have reduced susceptibility and decreased transmission (Hoffmann et al., 2011). This approach is currently ongoing in at least 12 countries and the first communications showing its impact on disease transmission decrease are being revealed (Nazni et al., 2019; Indriani et al., 2020; Ryan et al., 2020).

The mechanisms which are orchestrated by Wolbachia, conferring resistance to arboviruses infection remain to be elucidated, however, it has been assigned to the activation of the mosquito immune system, down-regulation of genes that encode proteins or receptors necessary to the virus binding and depletion of host’s resources essentials to the virus life cycle (Ford et al., 2020). When infected with Wolbachia, the mosquito Toll pathway is activated through the increase of reactive oxygen species (ROS), resulting in the production of antimicrobial peptides like cecropin (Pan et al., 2012). The gene expression modulation of tubulin, insulin receptors, and other genes seems to hinder the virus infection (Haqshenas et al., 2019; Paingankar et al., 2020). Also, it has been proposed that the lipids are relevant to Wolbachia infection. While one study implies it competes with the host’s lipid resources (mainly cholesterol), others suggest it modulates the lipid production which can have a role blocking virus replication (Geoghegan et al., 2017; Koh et al., 2020; Manokaran et al., 2020).

Protein quantification on a genome-wide scale allows a systemic bird’s-eye view of protein variation and pathway regulation under different conditions and samples (Choudhary and Mann, 2010). Therefore, it helps to highlight some features and emergent properties of complex systems, which could be obscured when the analysis is performed only from a reductionist point of view. In this way, mass spectrometry (MS)-based proteomics, together with robust bioinformatic tools, allows the formulation of more ambitious hypotheses due to the possibility of generating large and detailed datasets. Reliable large-scale protein quantification can be achieved using isobaric labeling and shotgun proteomics (Pappireddi et al., 2019). Several quantitative methods have already been successfully applied in insects to access differentially regulated proteins (De Mandal et al., 2020; García-Robles et al., 2020; Serteyn et al., 2020). The A. aegypti infection with Mayaro virus or microsporidian parasites have already been studied (Duncan et al., 2012; Vasconcellos et al., 2020) and a previous descriptive work was based on comparing the proteome of A. aegypti female heads during blood and sugar meals conditions (Nunes et al., 2016).

Here, we use an isobaric labeling-based quantitative proteomic strategy to interrogate the interaction between the Wolbachia wMel strain and ZIKV infection in A. aegypti heads and salivary glands. Wolbachia is present in several tissues and organs in A. aegypti mosquitoes, with higher density in ovaries and tissues on the head and anterior part of the thorax, like the ommatidia cells, brain, salivary glands, and cardia (Moreira et al., 2009). To the best of our knowledge, this is the most complete proteome of A. aegypti reported so far, with a total of 3,790 proteins identified, corresponding to 25.75% of total protein-coding genes. This work includes the identification and quantification of several peptides from the ZIKV and 323 Wolbachia proteins in a complex sample tissue, paving the way to the development of sensitive approaches to detect and quantify the ZIKV and Wolbachia directly in the A. aegypti tissue via MS. Moreover, we describe and discuss proteins and pathways altered in A. aegypti during ZIKV infections, Wolbachia infections, co-infection Wolbachia/ZIKV, and compared with no infection.

Materials and Methods

Mosquitoes

To access the effects of ZIKV and Wolbachia infection on A. aegypti proteome, we used mosquitoes from two different sites in Rio de Janeiro (Rio de Janeiro State, Brazil): Tubiacanga (22°47′06″S; 43°13′32″W) and Porto (22°53′43″S, 43°11′03″W), distant 13 km apart from each other. Tubiacanga was selected to represent an area in which mosquitoes have Wolbachia (wMel strain) since it is the first site in Latin America with an established invasion (>90% frequency) (Garcia et al., 2019). Porto represents a Wolbachia free area. Eggs were collected through 80 ovitraps (Codeço et al., 2015) installed every 25–50 m each other in both sites. Ovitraps were placed over an extensive geographic area to ensure we captured the local A. aegypti genetic variability, collecting at least 500 eggs per site. The eggs were hatched and the mosquitoes were maintained at the insectary under a relative humidity of 80 ± 5% and a temperature of 25 ± 3°C, with ad libitum access to a 10% sucrose solution. The experimental infection was performed with mosquitoes from the F1 generation.

Viral Strain

Aedes aegypti females were orally infected with the ZIKV strain Asian genotype isolated from a patient in Rio de Janeiro (GenBank accession number KU926309). Previous reports have shown high infectivity of this ZIKV isolate to Rio de Janeiro A. aegypti populations, including Porto, which had 100% infected females after receive an infective blood meal (Fernandes et al., 2016; da Silveira et al., 2018; Petersen et al., 2018). Viral titers in supernatants were previously determined by serial dilutions in Vero cells, expressed in plaque-forming units per milliliters (PFU/ml). All the assays were performed with samples containing 3.55 × 10⁶ PFU/ml. Viral stocks were maintained at −80°C until its use.

ZIKV Infection

Sugar supply was removed 36-h before mosquitoes were challenged with the infective blood meal to increase female’s avidity. Six–seven days old inseminated A. aegypti females from each of the two populations (Tubiacanga and Porto) were separated in cylindrical plastic cages for blood-feeding. Oral infection procedures were performed through a membrane feeding system (Hemotek, Great Harwood, United Kingdom), adapted with a pig-gut covering, which gives access to human blood. The infective blood meals consisted of 1 ml of the supernatant of infected Vero cell culture, 2 ml of washed rabbit erythrocytes, and 0.5 mM of adenosine triphosphate (ATP) as phagostimulant. The same procedure and membrane feeding apparatus were used to feed control mosquitoes, i.e., control insects received a noninfectious blood meal, with 1 ml of Vero cell culture medium replacing the viral supernatant (therefore, with no ZIKV), 2 ml of washed rabbit erythrocytes, and 0.5 mM of ATP as phagostimulant. After the experimental infection, we had a total of 152 A. aegypti females.

Ethical Approval

The maintenance of A. aegypti colonies with Wolbachia in the lab were achieved by offering blood obtained from anonymous donors from the blood bank of the Rio de Janeiro State University, whose blood bags would be discarded due to small volume. We have no information on donors, including sex, age, and clinical condition, but the blood bank discards those bags positive for Hepatitis B, Hepatitis C, Chagas disease, syphilis, human immunodeficiency virus, and human T-cell lymphotropic virus. Before offering the blood for mosquitoes, it was screened for Dengue virus (DENV) using the Dengue NS1 Ag STRIP (Bio-Rad). The use of human blood was approved by the Fiocruz Ethical Committee (CAAE 53419815.9.0000.5248).

Protein Extraction for Proteomics

Proteins from a total of 152 A. aegypti females were processed, of which 39 were coinfected with Wolbachia and ZIKV (WZ); 36 non-infected (A); 35 infected with Wolbachia (W) only and 42 infected with ZIKV (Z). On 14 days post-infection, each mosquito head plus the salivary gland was separated from the body using needles and a scalpel blade which were sterilized after every single use (Schmid et al., 2017). Proteins were extracted by lysis with a buffer containing 7 M urea, 2 M thiourea, 50 mM HEPES pH 8, 75 mM NaCl, 1 mM EDTA, 1 mM PMSF, and protease/phosphatase inhibitor cocktails (Roche). Then, samples were vortexed, sonicated in a cold bath for 10 min and sonicated by the probe for 10 min. Lysates were centrifuged at 10,000 g for 10 min at 4°C and the supernatants were carefully transferred to new tubes for further steps.

Protein Digestion and iTRAQ Labeling

The protein concentration was determined by the Qubit Protein Assay Kit^® fluorometric (Invitrogen), following the manufacturer’s instructions. A total of 100 μg of proteins from each condition was processed. Protein samples were incubated with dithiothreitol (DTT; GE Healthcare) at 10 mM and 30°C for 1 h. Subsequently, iodoacetamide solution (GE Healthcare) was added, with a final concentration of 40 mM. The reaction was performed at room temperature, in the dark, for 30 min. The samples were diluted 10× with 50 mM HEPES buffer to reduce the concentration of urea/thiourea. The solution was incubated with trypsin (Promega) in a 1:50 (w/w, enzyme/protein) ratio at 37°C for 18 h. The resulting peptides were desalted with a C-18 macro spin column (Harvard Apparatus) and then vacuum dried. The dried peptides were labeled with isobaric tags for relative and absolute quantitation (iTRAQ) 4-plex (ABSciex). The labeling process was performed according to Murillo et al. (2017). Briefly, dried peptides were dissolved in 50 mM TEAB and iTRAQ tags were dissolved in dry ethanol, labeling was performed by mixing peptide and tag solutions for 1 h at room temperature. The four tag reactions were combined in a single tube (ratio 1:1:1:1) and dried in a speed vac prior to fractionation. Each condition was labeled as follows: (i) Tag 114 corresponds to sample W (Wolbachia infected); (ii) Tag 115 to sample A (none infection); (iii) Tag 116 to sample WZ (Wolbachia and ZIKV co-infection); and (iv) Tag 117 to sample Z (ZIKV infection).

Fractionation by Hydrophilic Interaction Chromatography and LC-MS/MS Analysis

The iTRAQ labeled peptide mixture was fractionated off-line by HILIC (hydrophilic interaction liquid chromatography) before LC-MS/MS analysis. The dried samples were resuspended in acetonitrile (ACN) 90% / trifluoroacetic acid (TFA) 0.1% and injected into Shimadzu UFLC chromatography using a TSKGel Amide-80 column (15 cm × 2 mm i.d. × 3 μm – Supelco), flow rate 0.2 ml/min; mobile phases A (ACN 85%/TFA 0.1%) and B (TFA 0.1%); gradient: 0% of phase B in 5 min; 0–10% of phase B from 5 to 10 min; 10–20% of phase B from 10 to 55 min; 20–100% of phase B from 55 to 60 min; and 100–0% of phase B from 60 to 65 min. For each pooled sample, eight fractions were collected and combined according to the separation and intensity of the peaks. The pools of fractions were dried in a speed vac and resuspend in 0.1% formic acid (FA). The samples were analyzed in Easy-nLC 1000 coupled to a Q-Exactive Plus mass spectrometer (Thermo Fisher Scientific). For each fraction (or pool of fractions), a linear gradient of 5–40% of B em 105 min, 40–95% in 8 min and 95% B in 5 min was performed in 120 min; phase A (0.1% FA), B (ACN 95%, 0.1% FA). Ionization was performed in an electrospray source with the acquisition of spectra in positive polarity by data-dependent acquisition (DDA) mode, spray voltage of 2.5 kV, and temperature of 200°C in the heated capillary. The acquisition was set as follows: full scan or MS1 in the range of 375–1,800 m/z, resolution of 70,000 (m/z 200), fragmentation of the 10 most intense ions in the HCD collision cell, with standardized collision energy (NCE) of 30, resolution of 17,000 in the acquisition of MS/MS spectra, the first mass of 110 m/z, isolation window of 2.0 m/z and dynamic exclusion of 45 s.

Data Analysis

The raw data were processed using Proteome Discoverer 2.4 Software (Thermo Fisher Scientific). Peptide identification was performed with the Sequest HT search engine against A. aegypti (genome version/assembly ID: INSDC: GCA_002204515.1), ZIKV, and Wolbachia databases provided by VectorBase (Giraldo-Calderón et al., 2015), ViPR (Pickett et al., 2012), and UniProt (The Uniprot Consortium, 2019), respectively. Searches were performed with peptide mass tolerance of 10 ppm, MS/MS tolerance of 0.1 Da, tryptic cleavage specificity, two maximum missed cleavage sites allowed, fixed modification of carbamidomethyl (Cys), variable modification of iTRAQ 4-plex (Tyr, Lys, and peptide N-terminus), phosphate (Ser, Thr, and Tyr), and oxidation (Met); peptides with high confidence were considered. False discovery rates (FDR) were obtained using Target Decoy PSM selecting identifications with a q-value equal or to less than 0.01. Data of all technical replicates were log2-transformed and normalized by subtracting the column median. Relative iTRAQ quantification analysis was carried out in the Perseus software, version 1.6.12.0 (Tyanova et al., 2016), based on the intensity of the fragmented reporter peaks in MS/MS.

Statistics and Gene Enrichment

For ZIKV peptides abundance comparison, statistical analysis was performed in experimental groups using T-test analysis (95% confidence values, p < 0.05), meanwhile for differential proteins was used one-way analysis of variance (ANOVA) and Tukey’s multiple comparisons post-test (95% confidence values, p < 0.05) in Graphpad Prism 8.0.0 for Windows, GraphPad Software, San Diego, CA, United States, www.graphpad.com. Biological processes enrichment for differentially expressed proteins were performed using Fisher’s Exact Test (p < 0.05) in VectorBase website (bit.ly/38OmEX0) (Ashburner et al., 2000; Giraldo-Calderón et al., 2015). The resulting Gene Ontology (GO) lists were summarized and represented in a network interaction format in REVIGO (Supek et al., 2011).

Results and Discussion

Proteome Analysis

Proteomics has improved in recent years, especially considering its coverage and sensitivity, which offer more opportunities to study changes in the global proteome of complex tissues (Aebersold and Mann, 2016). Therefore, the use of isobaric-labeled quantitative proteomics measurements provides efficiency and increased depth coverage of proteomics analysis (Rauniyar and Yates, 2014). Peptides from the four experimental groups were labeled with iTRAQ 4-plex, enabling multiplexing of samples prior to off-line fractionation and LC-MS/MS analysis (Figure 1). We were able to identify 4,117 protein groups (Table 1), 26,898 peptides and 909,046 MS/MS spectra. The protein search was performed using A. aegypti, ZIKV, and Wolbachia databases simultaneously. Of all identified proteins, 3,790 belong to A. aegypti, 323 to Wolbachia and three unique peptides with a total of 10 PSMs to the ZIKV polyprotein (Figure 2 and Supplementary Figures 1, 2).

FIGURE 1

FIGURE 2

Quantitative Proteomics

The quantitative proteomics analysis was performed by comparing the four experimental groups using the normalized iTRAQ reporter ion intensity of identified peptides. Replicates were highly correlated (Pearson correlation mean of 0.97) and normalization methods were applied (Supplementary Figure 3). The 196 proteins with significant differences were statistically determined by the ANOVA test. Using the Tukey post-test ANOVA, we defined pairs of proteins with significant differences between the groups (Supplementary Table 2). The comparisons between the samples were: single-infected mosquitoes versus non-infected mosquitoes (W versus A and Z versus A); coinfected versus non-infected mosquitoes (WZ versus A); coinfected versus single-infected mosquitoes (WZ versus W and WZ versus Z).

Abundance of ZIKV Polyprotein Is Reduced in the Presence of Wolbachia

The quantitative analysis using the iTRAQ reporter ions from ZIKV peptides revealed that ZIKV proteins were more abundant in the absence of Wolbachia in contrast to mosquitoes coinfected with ZIKV and Wolbachia (Figure 2). This latter condition showed reporter ions intensity compared to the background signal of uninfected mosquitoes (Figure 2A). This result indicates Wolbachia presence likely reduces virus replication in the female head and salivary gland (Figure 2B). Noteworthy, the presence of the wMel strain in organs like salivary glands and ovaries of A. aegypti females induce the antiviral response that supports Wolbachia deployment to mitigate disease transmission.

The detection of ZIKV in the mosquito head + salivary gland was possible even considering that virus proteins were a minute fraction of all proteins identified. Peptides LITANPVITESTENS, SHTLWTDGIEESDLIIP, and LNQMSALEFYSYKK showed a robust increase in ZIKV infected tissue and this data represents a resource to further develop targeted proteomics methods to quantify ZIKV in infected mosquitoes (MS/MS spectra shown in Supplementary Figure 1).

Identification of Salivary Gland Proteins, Neuropeptides, and Hormones in Head Samples

During a blood meal, the female mosquitoes inject saliva into the host skin. The salivary glands produce a huge amount of molecules to avoid host responses against blood loss and to evade the host immune system (Ribeiro et al., 2016). For transmission of arboviroses, the virus must be secreted with the saliva during a blood meal in a vertebrate host (Linthicum et al., 1996). In this work, we took advantage of the sensitivity of protein detection by MS to detect and quantify proteins from the salivary gland, which are located on the anterior part of the thorax of A. aegypti. We compared the head proteome with two proteomes and one transcriptome (Wasinpiyamongkol et al., 2010; Ribeiro et al., 2016; Dhawan et al., 2017) from salivary glands of female mosquitoes to identify the salivary gland related proteins (Figure 3). We observed that among the 3,790 proteins identified in A. aegypti heads, 830 proteins were previously known to be expressed in the salivary glands. Moreover, we identified 77% of the proteins described by Wasinpiyamongkol et al. (2010) and 65% by Dhawan et al. (2017) in previous salivary gland proteome analysis. The 16 salivary proteins that were detected in all datasets probably are the most expressed and include inhibitors of platelet aggregation (salivary apyrase, AAEL006333-PA; apyrase precursor, AAEL006347PA), clotting inhibitors (SRPN23, AAEL002704-PB; SRPN26, AAEL003182-PA), aegyptins (AAEL010228-PA and AAEL010235-PA), angiopoietins (AAEL000749-PA and AAEL000726-PA), D7 family proteins (AAEL006417-PA and AAEL006424-PA), among others. We compared the 830 proteins previously known to be expressed in the salivary glands with the 196 differentially expressed proteins, to understand if these proteins could play a role during ZIKV and Wolbachia infections and/or co-infections. The modulation of saliva proteins is discussed below where the systemic changes in W, Z, and WZ will be addressed.

FIGURE 3

As a class of neuronal signal molecules, neuropeptides are produced by the neurosecretory cells, that are mainly located in the brain, suboesophageal ganglion, among others (Li et al., 2020). In insects, neuropeptides and their receptors play important physiological processes such as development, reproduction, behavior, and feeding (Holt et al., 2002). A previous study that investigated the presence of neuropeptidomics in A. aegypti proteome carried out in the central nervous system (CNS) and midgut samples identified 43 neuropeptides and hormones (Predel et al., 2010). Of these, we were able to identify 13 in our head proteomics analysis (Figure 4): pyrokinin (PK) (AAEL012060), short neuropeptide F (sNPF) (AAEL019691), neuropeptide-like precursor-1(NPLP) (AAEL012640), calcitonin-like diuretic hormone (DH-31) (AAEL008070), SIFamide (AAEL009858), crustacean cardioactive peptide (CCAP) (AAEL000630), allatotropin (AAEL009541), CAPA (AAEL005444), sulfakinin (SK) (AAEL006451), ovary ecdysteroidogenic hormone (OEH; neuroparsin A homolog) (AAEL004155), insulin-like peptide 1 and 7 (AAEL000937; AAEL024251), and adipokinetic hormone 1 (AKH) (AAEL011996). A few neuropeptides and hormones were modulated by co-infection and will be discussed afterward.

FIGURE 4

ZIKV Infection Modulates Aerobic Metabolism, Lipids Pathways, and Immune System Response

By comparing mosquitoes infected and non-infected with ZIKV (Figure 5 and Supplementary Figure 4A), we observe four main enriched processes: nicotinamide nucleotide metabolism, lipid transport, humoral innate response, and carbohydrate catabolism. Nicotinamide mononucleotide (NMN) is known for being an intermediator in NAD+ biosynthesis (Poddar et al., 2019; Hong et al., 2020), which participates actively in the aerobic metabolism process. This enrichment is mediated by the phosphoglycerate kinase (AAEL004988) and transaldolase (AAEL009389) upregulation, leading to carbohydrate catabolism elevation. The pentose phosphate pathway (PPP) was also up-regulated in our study as a response to transaldolase (AAEL009389) protein. The PPP is known to generate NADPH and pentoses (5-carbon sugars) as well as ribose 5-phosphate, a precursor for the synthesis of nucleotides, this led us to relate the activation of PPP with the probation of nucleotide biosynthesis for virus replication (Guo et al., 2019). Both increased fluxes have been observed following Human Cytomegalovirus, Mayaro, and Kaposi-sarcoma viruses in human cells (El-Bacha and Da Poian, 2013).

FIGURE 5

Viruses are capable of modulating mitochondrial bioenergetics aiming to synthesize macromolecular precursors for self-replication, assembly and egress (El-Bacha and Da Poian, 2013). Thus, the enhancement of the lipid transport process should drive a redirection of these energetic resources to the viral capsids formation (Kuzmina et al., 2020), an event that was also described in A. aegypti coinfected with DENV and Wolbachia (Koh et al., 2020). Koh et al. (2020) compared lipid abundances between non-infected and DENV-3 infected mosquitoes and found that most modulated lipids are glycerophospholipids, as described earlier (Perera et al., 2012), with a higher level of unsaturated fatty acids and triacylglycerols. Glycerophospholipids modulation can be related to a virus attempt to build new virus capsids, as these lipids are a compound of this structure (Martín-Acebes et al., 2016). In our results, three genes related to lipid localization and transport were up-regulated: apolipophorin-III (AAEL008789) and two vitellogenin-A1 precursors (AAEL006126 and AAEL006138). Apolipophorins are protein components of lipoprotein particles that are essential for lipid transportation in the insect body, but apolipophorin-III can work as a pattern recognition receptor for insect immune system response (Wang et al., 2019), including in A. aegypti (Phillips and Clark, 2017). Lipid carrier protein lipophorin (Lp) and its lipophorin receptor (LpRfb) were significantly increased in A. aegypti after infections with Gram (+) bacteria and fungi, suggest that lipid metabolism is involved in the mosquito systemic immune responses to pathogens and parasites, as seen in studies before using fat body tissue (Cheon et al., 2006). The apolipophorin-III upregulation can be rolled to lipid transport and also to the immune response as ZIKV is expected to regulate them and the protein participates in both processes.

Another group of up-regulated biological processes includes an innate immune response in host cells due to virus infection (antimicrobial humoral response, humoral immune response, and response to biotic stimulus). Xi et al. (2008) found a significant role for the Toll pathway in DENV resistance regulation, and it was possible to identify differentially expressed antimicrobial peptide cecropin (AAEL029038) related to this pathway in our samples, confirming that Toll-pathway plays a major part in virus response (Xi et al., 2008; Pan et al., 2012; Kingsolver et al., 2013).

We observed down-regulation of tomosyn (AAEL006948), adenylate cyclase (AAEL009314), putative syntaxin-7 (AAEL009398), putative sodium/dicarboxylate cotransporter (AAEL009863), cytochrome c oxidase subunit II (AAEL018664), major facilitator superfamily domain-containing protein (AAEL006002), nuclear pore complex protein Nup188 (AAEL008557). and talin-1 (AAEL024235) related to cation transport processes. Also, three proteins were down-regulated: calcium-transporting ATPase sarcoplasmic/endoplasmic reticulum type (AAEL006582), Na+/K+ ATPase alpha subunit (AAEL012062), and calmodulin (AAEL012326-PA) related to cation transport processes as well. Similarly cytoskeleton, vesicle-mediated transport and ion transport-associated proteins down-regulation were described in CHIKV infection (Cui et al., 2020). Also, a recent study showed that Na+/K+ ATPase inhibition using two distinct drugs in mice helped to block ZIKV replication (Guo et al., 2020), suggesting that this down-regulation could be a mosquito cellular response to fight infection.

Adenosine triphosphate-citrate synthase (AAEL004297), sphingosine phosphate lyase (AAEL003188) and probable cysteine-tRNA synthetase (AAEL022886), proteins associated with amino acids insertion in the citric acid cycle, were also down-regulated and may lead to a decrease in citric acid cycle and oxidative phosphorylation processes. This phenomenon was also observed in A. aegypti salivary glands infected with DENV-2 (Chisenhall et al., 2014; Shrinet et al., 2018), suggesting the occurrence of cellular biomolecule and energy production control by virus infection to avoid mitochondrial exhaustion.

Wolbachia Endosymbiosis Interfere in Glycoconjugates Production Pathways, Aerobic Metabolism, Immune Response Induction, and Blood-Feeding Process

It is possible to notice in our data (Figure 6 and Supplementary Figures 3B, 4B) an enrichment of GDP-mannose metabolism comparing mosquitoes non-infected with Wolbachia-mono infected. The enrichment was based on the upregulation of gdp mannose-4,6-dehydratase (AAEL014961), which converts GDP-mannose to GDP-4-dehydro-6-deoxy-D-mannose, the first of three steps for the conversion of GDP-mannose to GDP-fucose. This result suggests an induction caused by the presence of Wolbachia in insect cells to generate glycoconjugates, as GDP-mannose is a precursor for these molecules (Stick and Williams, 2009; Rexer et al., 2020), that may be used as elements to form new bacterial structures (García-del Portillo, 2020). An important observation was the enrichment of processes downregulated related to the citric acid cycle in mosquito cells infected with Wolbachia. As shown by da Rocha Fernandes et al. (2014), Wolbachia modulated glycogen metabolism in Aedes fluviatilis embryos in favor of its aerobic metabolism (da Rocha Fernandes et al., 2014).

FIGURE 6

Up-regulation of 4-hydroxyphenylpyruvate dioxygenase (AAEL014600) can be related to plastoquinol and vitamin E (tocopherol) synthesis. Both molecules are described in the literature as participants of antioxidant response to different stimuli in organisms (MaitiDutta et al., 2020; Nowicka et al., 2020), including oxidative stress (Kumar et al., 2020; Minter et al., 2020). In this way, we consider that genes correlated to these pathways were regulated during Wolbachia infection because this microorganism can induce the expression of antioxidant proteins (Brennan et al., 2008). Moreover, 4-hydroxyphenylpyruvate dioxygenase is related to L-tyrosine degradation I and L-phenylalanine degradation IV, which are essential traits for blood-feeding arthropods, as they consume an amount of blood proteins and need to deal with a high quantity of free amino acids (Sterkel et al., 2016).

The Wolbachia infection up-regulated 10 salivary gland proteins, such as endoplasmin (AAEL012827-PA), which is required for proper folding of Toll-like receptors in mammals (Fitzgerald and Kagan, 2020), and FKBP12 (AAEL010491-PB) that has a prolyl isomerase activity and bidins to immunosuppressive molecules (Caminati and Procacci, 2020). Twenty salivary proteins were down-regulated by Wolbachia infection. Interestingly, the same ion related proteins down-regulated by ZIKV are also down-regulated by Wolbachia, moreover, proteins related with contraction such as paramyosin (AAEL010975-PA), tropomyosin (AAEL002761-PA), and Actin-1 (AAEL001928-PA) were also down-regulated. Also, a serine protease inhibitor 25 (SRPN25) (AAEL007420-PB) was also identified in our analysis. Its orthologue is Alboserpin, the major salivary gland anticoagulant in Aedes albopictus, which prevents coagulation in an atypical reversible interaction with factor Xa (coagulation activation factor), important for blood intake by female mosquitoes (Stark and James, 1998; Calvo et al., 2011). SRPN25 was up-regulated in W, which may infer that blood takes longer to coagulate, increasing mosquito enzymes access to digest blood. The consequence of these salivary proteins modulation suggests an impairment in salivary gland injection into the host and a possible reduction in the amount of blood ingestion during Wolbachia presence in the mosquito (Turley et al., 2009).

Co-infection Modulated Pathways and Processes Show an Induction in Reactive Oxygen Species Production That Can Lead to an Immune System Response to Virus Infection

This analysis was based on a comparison of the ZW with Z (Figure 7 and Supplementary Figure 4C) and W (Figure 8 and Supplementary Figure 4D) mosquito heads samples.

FIGURE 7

FIGURE 8

Zika virus and Wolbachia use mosquitoes cell machinery respectively to synthesize lipids and glycoconjugates molecules for their structural use, as mentioned before. According to WZ results, the enriched lipid synthesis biological processes remain up-regulated despite the Wolbachia presence, similarly to ZIKV infection. Koh et al. (2020) confirmed such behavior during co-infection-induced intracellular events with DENV and Wolbachia, suggesting that they do not compete by lipids in host cell and virus-driven modulation dominates over that of Wolbachia (Koh et al., 2020), which may have occurred as well as in our experiment. In contrast, there is a down-regulation in GDP-mannose proteins genes synthesis, which may suggest that host cells no longer produce glycoconjugates for bacteria use.

Within the up-regulated genes in WZ compared to Z, ATP synthesis (ATPase subunit O; AAEL010823-PA), and cytoskeleton dynamics like microtubule binding and associated-proteins (AAEL009375-PK, AAEL004176-PB, and AAEL009847-PB) suggest activation of the ubiquitin-proteasome system. Indeed, it was previously reported in Drosophila and mosquito cell lines transfected, that Wolbachia relies on host proteolysis through ubiquitination to acquire essential amino acid for the infection (Fallon and Witthuhn, 2009; White et al., 2017). Flavivirus also need a functioning ubiquitin-proteasome system to complete their life cycle (Choy et al.,2015a,b), however, Ub3881 (AAEL003881), a ubiquitin-protein, is highly down-regulated in mosquitoes infected with DENV and its overexpression was able to control the virus infection (Troupin et al., 2016). In our analysis, we identified a paralog of the Ub3881 ubiquitin protein (AAEL003888), which was up-regulated in WZ versus Z and down-regulated in Z, suggesting that Wolbachia have a role in AAEL003888 expression enhancement which might target the envelope proteins of ZIKV for degradation, therefore hindering the virus assembly.

Exploring neuropeptides and hormones, neuropeptide sulfakinin (SK) (AAEL006451-PB) and insulin-like peptide (ILP) 1 (AAEL000937-PA) were down-regulated comparing WZ with Z. SK are a family of neuropeptides, homologous to mammalian gastrin/cholecystokinin (CCK), responsible for food uptake control in insects (Zels et al., 2015) while ILPs are known for signaling carbohydrates intake in cells (Nässel and Broeck, 2016). A recent study correlated SK influence in ILPs concentration, where a rise in SK concentration promotes ILPs production (Słocińska et al., 2020). SK down-regulation, observed in our data, may decrease in mosquito feeding satisfaction. This fact allied to ILP down-regulation, leading to a decrease in glycolysis cell assimilation, can be understood as a virus-infection influence to increase mosquito blood uptake, elevating ZIKV transmission success.

There are several processes up-regulated after enrichment comparison between WZ and Z. Cellular respiration, aerobic respiration, small molecule metabolic process, tricarboxylic acid cycle, oxidative phosphorylation, ATP synthesis coupled electron transport, respiratory electron transport chain, generation of precursor metabolites and energy, tricarboxylic acid metabolic process, nicotinamide nucleotide metabolic process, and cellular process were enriched up-regulated biological process by proteins genes ATP-citrate synthase (AAEL004297), glutamate semialdehyde dehydrogenase (AAEL006834), ATP synthase delta chain (AAEL010823), phosphatidylinositol-4,5-bisphosphate 4-phosphatase (AAEL011216), NADH-ubiquinone oxidoreductase (AAEL012552), mitochondrial aconitase (AAEL012897), 4-hydroxyphenylpyruvate dioxygenase (AAEL014600), titin (AAEL002565), and myosin (AAEL000596) up-regulation. The increase in aerobic respiration metabolism in Wolbachia infection proved to stimulate ROS production in insect cells as part of the host immune response, but it counterbalances with antioxidant pathways activation (Zug and Hammerstein, 2015). However, it was observed in WZ that the aerobic metabolism increases together with transition metal ions biological process and pathways, mediated by putative secreted ferritin G subunit precursor (AAEL007383) and transferrin (AAEL015458), which can also generate ROS (Liu et al., 2016). ROS likely activates the Toll-immune pathway to confront ZIKV, similar to DENV infection (Selivanov et al., 2008). Also, two eukaryotic elongation factor 1-alpha (eEF1A) proteins (AAEL017096; AAEL017301) were up-regulated comparing WZ and Z. eEF1A increase was described as a response to endoplasmic reticulum stress caused by ROS in Chinese hamster ovary (Borradaile et al., 2005), which may have occurred in the mosquito. Another observed data was SRPN25 down-regulated in WZ, which was up-regulated in W. SRPN25 down-regulation in WZ highlights its function related to the insect immune system. Serpins can block the Toll signaling cascade and β-1,3-glucan-mediated melanin biosynthesis (Jiang et al., 2009), so its downregulation can infer that the insect immune system is activated.

This led us to understand insects defense mechanisms, induced by the presence of endosymbiotic bacteria, that may have been acting against virus infection: melanization activation mediated by SRPN25 down-regulation and Toll-pathway with the influence of SPRN25 down-regulation and the exacerbated production of ROS by aerobic metabolism increase. ROS occurs inside the mitochondria, after the oxidative phosphorylation stage, which has oxygen as the final electron acceptor, forming metabolic water. When oxygen is prematurely and incompletely reduced, it gives rise to superoxide anion, classified as ROS, together with hydrogen peroxide, hypochlorous acid, hydroxyl radical, singlet oxygen, and ozone (Lambert and Brand, 2009). Pan et al. described the relation between Wolbachia presence in A. aegypti and Toll-pathway activation by ROS to control DENV (Pan et al., 2012). This study not only confirms oxidative stress caused by Wolbachia, but also antioxidant genes expression as described earlier in this article. The role of symbiotic microorganisms in arbovirus infection of arthropods vectors is widely discussed, highlighting the importance of a ROS-mediated stimulation of the Toll-pathway and AMPs expression (Yin et al., 2020). As SPRN25 is down-regulated, serpins activity is up-regulated, which lead to serpins key response in the defense mechanism of insects, activating especially the Toll pathway, as well as ROS, and prophenoloxidase (PPO) cascade (Meekins et al., 2017). As for PPO cascade, it mediates melanization immune response in insects, including mosquito vectors (Christensen et al., 2005), It was already investigated that melanization suppression by SRPN activation can lead to a increase in viral infection (Toufeeq et al., 2019), calling attention to SRPN importance in mosquito immune response.

Conclusion

Proteomics analysis of female’s head of A. aegypti during mono-infection with ZIKV or Wolbachia highlighted how those microorganisms use cell host for their benefits. Our data support that ZIKV may induces lipid synthesis and transport allied with glycolysis pathways while Wolbachia increases glycoconjugates production. During co-infection, Wolbachia probably helps A. aegypti to prevent virus infection by stimulating ROS production leading to Toll-pathway humoral immune response, together with antioxidant production to control cell homeostasis. This mechanism seems to be efficient since we have shown that peptides coming from the ZIKV polyprotein are reduced in female’s heads of A. aegypti in the presence of Wolbachia. Our study has provided important insights into the differential manner in which A. aegypti proteome is systemically regulated during mono- and co-infections of ZIKV and Wolbachia. Ultimately, this work represents a rich resource for the insect research community to help elucidate the mechanisms by which Wolbachia orchestrate resistance to ZIKV infection in A. aegypti and it allows the generation of specific hypothesis-driven experiments in future studies.

References

• 1

  AebersoldR.MannM. (2016). Mass-spectrometric exploration of proteome structure and function.Nature537347–355. 10.1038/nature19949

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 2

  AliotaM. T.WalkerE. C.YepesA. U.VelezI. D.ChristensenB. M.OsorioJ. E. (2016). The wMel strain of Wolbachia reduces transmission of Chikungunya virus in Aedes aegypti.PLoS Negl. Trop. Dis.10:e0004677. 10.1371/journal.pntd.0004677

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 3

  AshburnerM.BallC. A.BlakeJ. A.BotsteinD.ButlerH.CherryJ. M.et al (2000). Gene ontology: tool for the unification of biology.Nat. Genet.2525–29. 10.1038/75556

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 4

  BarretoM. L.Barral-NettoM.StabeliR.Almeida-FilhoN.VasconcelosP. F. C.TeixeiraM.et al (2016). Zika virus and microcephaly in Brazil: a scientific agenda.Lancet387919–921. 10.1016/S0140-6736(16)00545-6

  • CrossRef
  • Google Scholar

• 5

  BaudD.GublerD. J.SchaubB.LanteriM. C.MussoD. (2017). An update on Zika virus infection.Lancet3902099–2109. 10.1016/S0140-6736(17)31450-2

  • CrossRef
  • Google Scholar

• 6

  BorradaileN. M.BuhmanK. K.ListenbergerL. L.MageeC. J.MorimotoE. T. A.OryD. S.et al (2005). A critical role for eukaryotic elongation factor 1A-1 in lipotoxic cell death.MBoC17770–778. 10.1091/mbc.e05-08-0742

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 7

  BrennanL. J.KeddieB. A.BraigH. R.HarrisH. L. (2008). The endosymbiont Wolbachia pipientis induces the expression of host antioxidant proteins in an Aedes albopictus cell line.PLoS One3:e2083. 10.1371/journal.pone.0002083

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 8

  CalvoE.MizuriniD. M.Sá-NunesA.RibeiroJ. M. C.AndersenJ. F.MansB. J.et al (2011). Alboserpin, a factor Xa inhibitor from the mosquito vector of yellow fever, binds heparin and membrane phospholipids and exhibits antithrombotic activity.J. Biol. Chem.28627998–28010. 10.1074/jbc.M111.247924

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 9

  CaminatiG.ProcacciP. (2020). Mounting evidence of FKBP12 implication in neurodegeneration.Neural Regen. Res.152195–2202. 10.4103/1673-5374.284980

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 10

  CamposG.BandeiraA.SardiS. (2015). Zika virus outbreak, Bahia, Brazil.Emerg. Infect. Dis. J.21:1885. 10.3201/eid2110.150847

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 11

  CheonH.-M.ShinS. W.BianG.ParkJ.-H.RaikhelA. S. (2006). Regulation of lipid metabolism genes, lipid carrier protein lipophorin, and its receptor during immune challenge in the mosquito Aedes aegypti.J. Biol. Chem.2818426–8435. 10.1074/jbc.M510957200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 12

  ChisenhallD. M.LondonoB. L.ChristoffersonR. C.McCrackenM. K.MoresC. N. (2014). Effect of dengue-2 virus infection on protein expression in the salivary glands of Aedes aegypti mosquitoes.Am. J. Trop. Med. Hyg.90431–437. 10.4269/ajtmh.13-0412

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 13

  ChoudharyC.MannM. (2010). Decoding signalling networks by mass spectrometry-based proteomics.Nat. Rev. Mol. Cell Biol.11427–439. 10.1038/nrm2900

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 14

  Chouin-CarneiroT.Vega-RuaA.VazeilleM.YebakimaA.GirodR.GoindinD.et al (2016). Differential susceptibilities of Aedes aegypti and Aedes albopictus from the Americas to Zika Virus.PLoS Negl. Trop. Dis.10:e0004543. 10.1371/journal.pntd.0004543

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 15

  ChoyM. M.SessionsO. M.GublerD. J.OoiE. E. (2015a). Production of infectious dengue virus in Aedes aegypti is dependent on the ubiquitin proteasome pathway.PLoS Negl. Trop. Dis.9:e0004227. 10.1371/journal.pntd.0004227

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 16

  ChoyM. M.ZhangS. L.CostaV. V.TanH. C.HorrevortsS.OoiE. E. (2015b). Proteasome inhibition suppresses dengue virus egress in antibody dependent infection.PLoS Negl. Trop. Dis.9:e0004058. 10.1371/journal.pntd.0004058

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 17

  ChristensenB. M.LiJ.ChenC.-C.NappiA. J. (2005). Melanization immune responses in mosquito vectors.Trends Parasitol.21192–199. 10.1016/j.pt.2005.02.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 18

  CodeçoC. T.LimaA. W. S.AraújoS. C.LimaJ. B. P.Maciel-de-FreitasR.HonórioN. A.et al (2015). Surveillance of Aedes aegypti: comparison of house index with four alternative traps.PLoS Negl. Trop. Dis.9:e0003475. 10.1371/journal.pntd.0003475

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 19

  CuiY.LiuP.MooneyB. P.FranzA. W. E. (2020). Quantitative proteomic analysis of Chikungunya virus-infected Aedes aegypti reveals proteome modulations indicative of persistent infection.J. Proteome Res.192443–2456. 10.1021/acs.jproteome.0c00173

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 20

  da Rocha FernandesM.MartinsR.Pessoa CostaE.PacidônioE. C.Araujoet al (2014). The modulation of the symbiont/host interaction between Wolbachia pipientis and Aedes fluviatilis embryos by glycogen metabolism.PLoS One9:e98966. 10.1371/journal.pone.0098966

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 21

  da SilveiraI. D.PetersenM. T.SylvestreG.GarciaG. A.DavidM. R.PavanM. G.et al (2018). Zika virus infection produces a reduction on Aedes aegypti lifespan but no effects on mosquito fecundity and oviposition success.Front. Microbiol.9:3011. 10.3389/fmicb.2018.03011

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 22

  De MandalS.LinB.ShiM.LiY.XuX.JinF. (2020). iTRAQ-based comparative proteomic analysis of larval midgut from the beet armyworm, Spodoptera exigua (Hübner) (Lepidoptera: Noctuidae) challenged with the entomopathogenic bacteria Serratia marcescens.Front. Physiol.11:442. 10.3389/fphys.2020.00442

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 23

  DhawanR.MohantyA. K.KumarM.DeyG.AdvaniJ.PrasadT. S. K.et al (2017). Data from salivary gland proteome analysis of female Aedes aegypti Linn.Data Brief.13274–277. 10.1016/j.dib.2017.05.034

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 24

  DickG. W. A.KitchenS. F.HaddowA. J. (1952). Zika Virus (I). Isolations and serological specificity.Trans. R. Soc. Trop. Med. Hyg.46509–520. 10.1016/0035-9203(52)90042-4

  • CrossRef
  • Google Scholar

• 25

  DuncanA. B.AgnewP.NoelV.DemettreE.SevenoM.BrizardJ.-P.et al (2012). Proteome of Aedes aegypti in response to infection and coinfection with microsporidian parasites.Ecol. Evol.2681–694. 10.1002/ece3.199

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 26

  DutraH. L. C.RochaM. N.DiasF. B. S.MansurS. B.CaragataE. P.MoreiraL. A. (2016). Wolbachia blocks currently circulating Zika virus isolates in Brazilian Aedes aegypti mosquitoes.Cell Host Microbe19771–774. 10.1016/j.chom.2016.04.021

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 27

  El-BachaT.Da PoianA. T. (2013). Virus-induced changes in mitochondrial bioenergetics as potential targets for therapy.Int. J. Biochem. Cell Biol.4541–46. 10.1016/j.biocel.2012.09.021

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 28

  El-GebaliS.MistryJ.BatemanA.EddyS. R.LucianiA.PotterS. C.et al (2018). The Pfam protein families database in 2019.Nucleic Acids Res.47427–432. 10.1093/nar/gky995

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 29

  FallonA. M.WitthuhnB. A. (2009). Proteasome activity in a naïve mosquito cell line infected with Wolbachia pipientis wAlbB.Vitro Cell. Dev. Biol. Anim.45460–466. 10.1007/s11626-009-9193-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 30

  FernandesR. S.CamposS. S.Ferreira-de-BritoA.MirandaR. M.Barbosa, da SilvaK. A.et al (2016). Culex quinquefasciatus from Rio de janeiro is not competent to transmit the local Zika Virus.PLoS Negl. Trop. Dis.10:e0004993. 10.1371/journal.pntd.0004993

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 31

  Ferreira-de-BritoA.RibeiroI. P.MirandaR. M.FernandesR. S.CamposS. S.SilvaK. A.et al (2016). First detection of natural infection of Aedes aegypti with Zika virus in Brazil and throughout South America.Memórias Instit. Oswaldo Cruz111655–658. 10.1590/0074-02760160332

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 32

  FitzgeraldK. A.KaganJ. C. (2020). Toll-like receptors and the control of immunity.Cell1801044–1066. 10.1016/j.cell.2020.02.041

  • CrossRef
  • Google Scholar

• 33

  FordS. A.AlbertI.AllenS. L.ChenowethS. F.JonesM.KohC.et al (2020). Artificial selection finds new hypotheses for the mechanism of Wolbachia-Mediated dengue blocking in mosquitoes.Front. Microbiol.11:1456. 10.3389/fmicb.2020.01456

  • CrossRef
  • Google Scholar

• 34

  GarciaG. A.SylvestreG.AguiarR.da CostaG. B.MartinsA. J.LimaJ. B. P.et al (2019). Matching the genetics of released and local Aedes aegypti populations is critical to assure Wolbachia invasion.PLoS Negl. Trop. Dis.13:e0007023. 10.1371/journal.pntd.0007023

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 35

  García-RoblesI.De LomaJ.CapillaM.RogerI.Boix-MontesinosP.CarriónP.et al (2020). Proteomic insights into the immune response of the Colorado potato beetle larvae challenged with Bacillus thuringiensis.Dev. Comp. Immunol.104:103525. 10.1016/j.dci.2019.103525

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 36

  García-del PortilloF. (2020). Building peptidoglycan inside eukaryotic cells: a view from symbiotic and pathogenic bacteria.Mol. Microbiol.113613–626. 10.1111/mmi.14452

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 37

  GeogheganV.StaintonK.RaineyS. M.AntT. H.DowleA. A.LarsonT.et al (2017). Perturbed cholesterol and vesicular trafficking associated with dengue blocking in Wolbachia -infected Aedes aegypti cells.Nat. Commun.8:526. 10.1038/s41467-017-00610-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 38

  Giraldo-CalderónG. I.EmrichS. J.MacCallumR. M.MaslenG.DialynasE.TopalisP.et al (2015). VectorBase: an updated bioinformatics resource for invertebrate vectors and other organisms related with human diseases.Nucleic Acids Res.43707–713. 10.1093/nar/gku1117

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 39

  GuoJ.JiaX.LiuY.WangS.CaoJ.ZhangB.et al (2020). Inhibition of Na + /K + ATPase blocks Zika virus infection in mice.Commun. Biol.3380. 10.1038/s42003-020-1109-8

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 40

  GuoX.WuS.LiN.LinQ.LiuL.LiangH.et al (2019). Accelerated metabolite levels of aerobic glycolysis and the pentose phosphate pathway are required for efficient replication of infectious spleen and kidney necrosis virus in chinese perch brain cells.Biomolecules9:440. 10.3390/biom9090440

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 41

  HaqshenasG.TerradasG.ParadkarP. N.DucheminJ.-B.McGrawE. A.DoerigC. (2019). A role for the insulin receptor in the suppression of dengue virus and Zika virus in wolbachia-infected mosquito cells.Cell Rep.26529.e3–535.e3. 10.1016/j.celrep.2018.12.068

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 42

  HoffmannA. A.MontgomeryB. L.PopoviciJ.Iturbe-OrmaetxeI.JohnsonP. H.MuzziF.et al (2011). Successful establishment of Wolbachia in Aedes populations to suppress dengue transmission.Nature476454–457. 10.1038/nature10356

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 43

  HoltR. A.SubramanianG. M.HalpernA.SuttonG. G.CharlabR.NusskernD. R.et al (2002). The genome sequence of the malaria mosquito Anopheles gambiae.Science298129–149. 10.1126/science.1076181

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 44

  HongW.MoF.ZhangZ.HuangM.WeiX. (2020). Nicotinamide mononucleotide: a promising molecule for therapy of diverse diseases by targeting NAD+ metabolism.Front. Cell Dev. Biol.8:246. 10.3389/fcell.2020.00246

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 45

  IndrianiC.TantowijoyoW.RancèsE.AndariB.PrabowoE.YusdiD.et al (2020). Reduced dengue incidence following deployments of Wolbachia-infected Aedes aegypti in Yogyakarta, Indonesia: a quasi-experimental trial using controlled interrupted time series analysis.Gates Open Res.4:50. 10.12688/gatesopenres.13122.1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 46

  JiangR.KimE.-H.GongJ.-H.KwonH.-M.KimC.-H.RyuK.-H.et al (2009). Three pairs of protease-serpin complexes cooperatively regulate the insect innate immune responses.J. Biol. Chem.28435652–35658. 10.1074/jbc.M109.071001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 47

  KingsolverM. B.HuangZ.HardyR. W. (2013). Insect antiviral innate immunity: pathways, effectors, and connections.J. Mol. Biol.4254921–4936. 10.1016/j.jmb.2013.10.006

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 48

  KohC.IslamM. N.YeY. H.ChotiwanN.GrahamB.BelisleJ. T.et al (2020). Dengue virus dominates lipid metabolism modulations in Wolbachia -coinfected Aedes aegypti.Commun. Biol.3:e0007443. 10.1038/s42003-020-01254-z

  • CrossRef
  • Google Scholar

• 49

  KumarA.PrasadA.SedlářováM.KsasB.HavauxM.PospíšilP. (2020). Interplay between antioxidants in response to photooxidative stress in Arabidopsis.Free Radic. Biol. Med.160894–907. 10.1016/j.freeradbiomed.2020.08.027

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 50

  KuzminaN. V.LoshkarevaA. S.ShilovaL. A.ShtykovaE. V.KnyazevD. G.ZimmerbergJ.et al (2020). Protein-lipid interactions in formation of viral envelopes.Biophys. J.118:523a. 10.1016/j.bpj.2019.11.2874

  • CrossRef
  • Google Scholar

• 51

  LambertA. J.BrandM. D. (2009). “Reactive oxygen species production by mitochondria,” in Mitochondrial DNA: Methods and Protocols Methods in Molecular Biology^TM, ed.StuartJ. A. (Totowa, NJ: Humana Press), 165–181. 10.1007/978-1-59745-521-3_11

  • CrossRef
  • Google Scholar

• 52

  LiF.ZhaoX.ZhuS.WangT.LiT.WoolfleyT.et al (2020). Identification and expression profiling of neuropeptides and neuropeptide receptor genes in Atrijuglans hetaohei.Gene743:144605. 10.1016/j.gene.2020.144605

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 53

  LinthicumK. J.PlattK.MyintK. S.LerdthusneeK.InniB. L.GhnD. W. V. (1996). Dengue 3 virus distribution in the mosquito Aedes aegypti: an immunocytochemical study.Medic. Vet. Entomol.1087–92. 10.1111/j.1365-2915.1996.tb00086.x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 54

  LiuK.GaoY.LiuJ.WenY.ZhaoY.ZhangK.et al (2016). Photoreactivity of metal–organic frameworks in aqueous solutions: metal dependence of reactive oxygen species production.Environ. Sci. Technol.503634–3640. 10.1021/acs.est.5b06019

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 55

  MaitiDuttaS.ChenG.MaitiS. (2020). Tocopherol moderately induces the expressions of some human sulfotransferases, which are activated by oxidative stress.Cell Biochem. Biophys.78439–446. 10.1007/s12013-020-00938-x

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 56

  ManokaranG.FloresH. A.DicksonC. T.NarayanaV. K.KanojiaK.DayalanS.et al (2020). Modulation of acyl-carnitines, the broad mechanism behind Wolbachia-mediated inhibition of medically important flaviviruses in Aedes aegypti.PNAS11724475–24483. 10.1073/pnas.1914814117

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 57

  Martín-AcebesM. A.Vázquez-CalvoÁSaizJ.-C. (2016). Lipids and flaviviruses, present and future perspectives for the control of dengue, Zika, and West Nile viruses.Prog. Lipid Res.64123–137. 10.1016/j.plipres.2016.09.005

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 58

  MeekinsD. A.KanostM. R.MichelK. (2017). Serpins in arthropod biology.Semin. Cell Dev. Biol.62105–119. 10.1016/j.semcdb.2016.09.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 59

  MinterB. E.LowesD. A.WebsterN. R.GalleyH. F. (2020). Differential effects of MitoVitE, α-Tocopherol and Trolox on oxidative stress, mitochondrial function and inflammatory signalling pathways in endothelial cells cultured under conditions mimicking sepsis.Antioxidants9:195. 10.3390/antiox9030195

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 60

  MoreiraL. A.Iturbe-OrmaetxeI.JefferyJ. A.LuG.PykeA. T.HedgesL. M.et al (2009). A Wolbachia symbiont in Aedes aegypti limits infection with dengue, chikungunya, and plasmodium.Cell1391268–1278. 10.1016/j.cell.2009.11.042

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 61

  MurilloJ. R.Goto-SilvaL.SánchezA.NogueiraF. C. S.DomontG. B.JunqueiraM. (2017). Quantitative proteomic analysis identifies proteins and pathways related to neuronal development in differentiated SH-SY5Y neuroblastoma cells.EuPA Open Proteom.161–11. 10.1016/j.euprot.2017.06.001

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 62

  NässelD. R.BroeckJ. V. (2016). Insulin/IGF signaling in Drosophila and other insects: factors that regulate production, release and post-release action of the insulin-like peptides.Cell. Mol. Life Sci.73271–290. 10.1007/s00018-015-2063-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 63

  NazniW. A.HoffmannA. A.NoorAfizahA.CheongY. L.ManciniM. V.GoldingN.et al (2019). Establishment of Wolbachia strain wAlbB in Malaysian populations of Aedes aegypti for dengue control.Curr. Biol.294241.e5–4248.e5. 10.1016/j.cub.2019.11.007

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 64

  NowickaB.FesenkoT.WalczakJ.KrukJ. (2020). The inhibitor-evoked shortage of tocopherol and plastoquinol is compensated by other antioxidant mechanisms in Chlamydomonas reinhardtii exposed to toxic concentrations of cadmium and chromium ions.Ecotoxicol. Environ. Saf.191:110241. 10.1016/j.ecoenv.2020.110241

  • CrossRef
  • Google Scholar

• 65

  NunesA. T.BritoN. F.OliveiraD. S.AraujoG. D. T.NogueiraF. C. S.DomontG. B.et al (2016). Comparative proteome analysis reveals that blood and sugar meals induce differential protein expression in Aedes aegypti female heads.Proteomics162582–2586. 10.1002/pmic.201600126

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 66

  PaingankarM. S.GokhaleM. D.DeobagkarD. D.DeobagkarD. N. (2020). Drosophila melanogaster as a model host to study arbovirus–vector interaction.bioRxiv10.1101/2020.09.03.282350

  • CrossRef
  • Google Scholar

• 67

  PanX.ZhouG.WuJ.BianG.LuP.RaikhelA. S.et al (2012). Wolbachia induces reactive oxygen species (ROS)-dependent activation of the Toll pathway to control dengue virus in the mosquito Aedes aegypti.PNAS109E23–E31. 10.1073/pnas.1116932108

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 68

  PappireddiN.MartinL.WührM. (2019). A review on quantitative multiplexed proteomics.ChemBioChem201210–1224. 10.1002/cbic.201800650

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 69

  PereraR.RileyC.IsaacG.Hopf-JannaschA. S.MooreR. J.WeitzK. W.et al (2012). Dengue virus infection perturbs lipid homeostasis in infected mosquito cells.PLoS Pathog.8:e1002584. 10.1371/journal.ppat.1002584

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 70

  PetersenM. T.da SilveiraI. D.Tátila-FerreiraA.DavidM. R.Chouin-CarneiroT.Van den WouwerL.et al (2018). The impact of the age of first blood meal and Zika virus infection on Aedes aegypti egg production and longevity.PLoS One13:e0200766. 10.1371/journal.pone.0200766

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 71

  PhillipsD. R.ClarkK. D. (2017). Bombyx mori and Aedes aegypti form multi-functional immune complexes that integrate pattern recognition, melanization, coagulants, and hemocyte recruitment.PLoS One12:e0171447. 10.1371/journal.pone.0171447

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 72

  PickettB. E.SadatE. L.ZhangY.NoronhaJ. M.SquiresR. B.HuntV.et al (2012). ViPR: an open bioinformatics database and analysis resource for virology research.Nucleic Acids Res.40D593–D598. 10.1093/nar/gkr859

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 73

  PoddarS. K.SifatA. E.HaqueS.NahidN. A.ChowdhuryS.MehediI. (2019). Nicotinamide mononucleotide: exploration of diverse therapeutic applications of a potential molecule.Biomolecules9:34. 10.3390/biom9010034

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 74

  PredelR.NeupertS.GarczynskiS. F.CrimJ. W.BrownM. R.RussellW. K.et al (2010). Neuropeptidomics of the mosquito Aedes aegypti.J. Proteome Res.92006–2015. 10.1021/pr901187p

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 75

  RauniyarN.YatesJ. R. (2014). Isobaric labeling-based relative quantification in shotgun proteomics.Proteome Res.135293–5309. 10.1021/pr500880b

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 76

  RexerT.LaafD.GottschalkJ.FrohnmeyerH.RappE.EllingL. (2020). “Enzymatic synthesis of glycans and glycoconjugates,” in Advances in Biochemical Engineering/Biotechnology, (Berlin: Springer), 1–50. 10.1007/10_2020_148

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 77

  RibeiroJ. M. C.Martin-MartinI.ArcàB.CalvoE. (2016). A deep insight into the sialome of male and female Aedes aegypti mosquitoes.PLoS One11:e0151400. 10.1371/journal.pone.0151400

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 78

  RyanP. A.TurleyA. P.WilsonG.HurstT. P.RetzkiK.Brown-KenyonJ.et al (2020). Establishment of wMel Wolbachia in Aedes aegypti mosquitoes and reduction of local dengue transmission in cairns and surrounding locations in northern Queensland, Australia.Gates Open Res.3:1547. 10.12688/gatesopenres.13061.2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 79

  SchmidM. A.KauffmanE.PayneA.HarrisE.KramerL. D. (2017). Preparation of mosquito salivary gland extract and intradermal inoculation of mice.Bio Protocol7:e2407.

  • Google Scholar

• 80

  SelivanovV. A.ZeakJ. A.RocaJ.CascanteM.TruccoM.VotyakovaT. V. (2008). The role of external and matrix pH in mitochondrial reactive oxygen species generation.J. Biol. Chem.28329292–29300. 10.1074/jbc.M801019200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 81

  SerteynL.PonnetL.SaiveM.FauconnierM.-L.FrancisF. (2020). Changes of feeding behavior and salivary proteome of Brown Marmorated Stink Bug when exposed to insect-induced plant defenses.Arthropod-Plant Interact.14101–112. 10.1007/s11829-019-09718-8

  • CrossRef
  • Google Scholar

• 82

  ShrinetJ.BhaveshN. S.SunilS. (2018). Understanding oxidative stress in aedes during chikungunya and dengue virus infections using integromics analysis.Viruses10:314. 10.3390/v10060314

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 83

  SłocińskaM.ChowańskiS.MarciniakP. (2020). Identification of sulfakinin receptors (SKR) in Tenebrio molitor beetle and the influence of sulfakinins on carbohydrates metabolism.J. Comp. Physiol. B190669–679. 10.1007/s00360-020-01300-6

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 84

  StarkK. R.JamesA. A. (1998). Isolation and characterization of the gene encoding a novel factor xa-directed anticoagulant from the yellow fever mosquito, Aedes aegypti.J. Biol. Chem.27320802–20809. 10.1074/jbc.273.33.20802

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 85

  SterkelM.PerdomoH. D.GuizzoM. G.BarlettaA. B. F.NunesR. D.DiasF. A.et al (2016). Tyrosine detoxification is an essential trait in the life history of blood-feeding arthropods.Curr. Biol.262188–2193. 10.1016/j.cub.2016.06.025

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 86

  StickR. V.WilliamsS. J. (2009). “Chapter 11 - Glycoproteins and proteoglycans,” in Carbohydrates: The Essential Molecules of Life (Second Edition), edsStickR. V.WilliamsS. J. (Oxford: Elsevier), 369–412.

  • Google Scholar

• 87

  SupekF.BošnjakM.ŠkuncaN.ŠmucT. (2011). REVIGO summarizes and visualizes long lists of gene ontology terms.PLoS One6:e21800. 10.1371/journal.pone.0021800

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 88

  The Uniprot Consortium (2019). UniProt: a worldwide hub of protein knowledge.Nucleic Acids Res.47D506–D515. 10.1093/nar/gky1049

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 89

  ToufeeqS.WangJ.ZhangS.-Z.LiB.HuP.ZhuL.-B.et al (2019). Bmserpin2 is involved in BmNPV infection by suppressing melanization in Bombyx mori.Insects10:399. 10.3390/insects10110399

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 90

  TroupinA.Londono-RenteriaB.ConwayM. J.ClohertyE.JamesonS.HiggsS.et al (2016). A novel mosquito ubiquitin targets viral envelope protein for degradation and reduces virion production during dengue virus infection.Biochim. Biophys. Acta Gen. Subj.18601898–1909. 10.1016/j.bbagen.2016.05.033

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 91

  TurleyA. P.MoreiraL. A.O’NeillS. L.McGrawE. A. (2009). Wolbachia infection reduces blood-feeding success in the dengue fever mosquito, Aedes aegypti.PLoS Negl. Trop. Dis.3:e516. 10.1371/journal.pntd.0000516

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 92

  TyanovaS.TemuT.SinitcynP.CarlsonA.HeinM. Y.GeigerT.et al (2016). The Perseus computational platform for comprehensive analysis of (prote)omics data.Nat. Methods13731–740. 10.1038/nmeth.3901

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 93

  VasconcellosA. F.MandacaruS. C.de OliveiraA. S.FontesW.MeloR. M.de SousaM. V.et al (2020). Dynamic proteomic analysis of Aedes aegypti Aag-2 cells infected with Mayaro virus.Parasit. Vect.13:297. 10.1186/s13071-020-04167-2

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 94

  WangX.ZhangY.ZhangR.ZhangJ. (2019). The diversity of pattern recognition receptors (PRRs) involved with insect defense against pathogens.Curr. Opin. Insect Sci.33105–110. 10.1016/j.cois.2019.05.004

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 95

  WasinpiyamongkolL.PatramoolS.LuplertlopN.SurasombatpattanaP.DoucoureS.MouchetF.et al (2010). Blood-feeding and immunogenic Aedes aegypti saliva proteins.Proteomics101906–1916. 10.1002/pmic.200900626

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 96

  WhiteP. M.SerbusL. R.DebecA.CodinaA.BrayW.GuichetA.et al (2017). Reliance of Wolbachia on high rates of host proteolysis revealed by a genome-wide RNAi screen of Drosophila cells.Genetics2051473–1488. 10.1534/genetics.116.198903

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 97

  XiZ.RamirezJ. L.DimopoulosG. (2008). The Aedes aegypti toll pathway controls dengue virus infection.PLoS Pathog.4:e1000098. 10.1371/journal.ppat.1000098

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 98

  YinC.SunP.YuX.WangP.ChengG. (2020). Roles of symbiotic microorganisms in arboviral infection of arthropod vectors.Trends Parasitol.36607–615. 10.1016/j.pt.2020.04.009

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 99

  ZelsS.DillenS.CrabbéK.SpitJ.NachmanR. J.Vanden BroeckJ. (2015). Sulfakinin is an important regulator of digestive processes in the migratory locust, Locusta migratoria.Insect Biochem. Mol. Biol.618–16. 10.1016/j.ibmb.2015.03.008

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 100

  ZugR.HammersteinP. (2015). Wolbachia and the insect immune system: what reactive oxygen species can tell us about the mechanisms of Wolbachia–host interactions.Front. Microbiol.6:1201. 10.3389/fmicb.2015.01201

  • Pubmed Abstract
  • CrossRef
  • Google Scholar