Advertisement

  • Loading metrics

Open Access

Peer-reviewed

Research Article

Dengue virus reduces AGPAT1 expression to alter phospholipids and enhance infection in Aedes aegypti

Abstract

More than half of the world population is at risk of dengue virus (DENV) infection because of the global distribution of its mosquito vectors. DENV is an envelope virus that relies on host lipid membranes for its life-cycle. Here, we characterized how DENV hijacks the mosquito lipidome to identify targets for novel transmission-blocking interventions. To describe metabolic changes throughout the mosquito DENV cycle, we deployed a Liquid chromatography–high resolution mass spectrometry (LC-HRMS) workflow including spectral similarity annotation in cells, midguts and whole mosquitoes at different times post infection. We revealed a major aminophospholipid reconfiguration with an overall early increase, followed by a reduction later in the cycle. We phylogenetically characterized acylglycerolphosphate acyltransferase (AGPAT) enzyme isoforms to identify those that catalyze a rate-limiting step in phospholipid biogenesis, the acylation of lysophosphatidate to phosphatidate. We showed that DENV infection decreased AGPAT1, but did not alter AGPAT2 expression in cells, midguts and mosquitoes. Depletion of either AGPAT1 or AGPAT2 increased aminophospholipids and partially recapitulated DENV-induced reconfiguration before infection in vitro. However, only AGPAT1 depletion promoted infection by maintaining high aminophospholipid concentrations. In mosquitoes, AGPAT1 depletion also partially recapitulated DENV-induced aminophospholipid increase before infection and enhanced infection by maintaining high aminophospholipid concentrations. These results indicate that DENV inhibition of AGPAT1 expression promotes infection by increasing aminophospholipids, as observed in the mosquito’s early DENV cycle. Furthermore, in AGPAT1-depleted mosquitoes, we showed that enhanced infection was associated with increased consumption/redirection of aminophospholipids. Our study suggests that DENV regulates aminophospholipids, especially phosphatidylcholine and phosphatidylethanolamine, by inhibiting AGPAT1 expression to increase aminophospholipid availability for virus multiplication.

Author summary

Dengue is endemic in tropical and subtropical regions, and has now encroached onto temperate regions because of the geographic expansion of its vector, Aedes aegypti. In the absence of effective vaccine and curative drug, the sole intervention relies on containment strategies using insecticide. However, occurrence of insecticide resistance diminishes vector control efficacy. Here, we explore the nascent field of mosquito metabolomics as part of our discovery effort for new transmission-blocking targets. Dengue virus (DENV) relies on host metabolome, specifically the lipid membrane to complete its life-cycle. However, little is known about how DENV subverts the mosquito physiology. Using high-resolution mass spectrometry, we described metabolic changes incurred by DENV throughout the mosquito cycle, from cellular replication onset to systemic infection. Membrane phospholipids were highly reconfigured and were associated with reduced expression of AGPAT1, an enzyme involved in their biogenesis. AGPAT1 depletion partially recapitulated DENV-induced metabolic reconfiguration and enhanced infection by maintaining high phospholipid concentrations. These phospholipids were then consumed/redirected later in the mosquito DENV cycle. Our work comprehensively describes metabolic changes associated with DENV infection. In addition, we reveal how DENV subdues the lipidome for its benefit by demonstrating the role of phospholipids in mosquito infection.

Citation: Vial T, Tan W-L, Wong Wei Xiang B, Missé D, Deharo E, Marti G, et al. (2019) Dengue virus reduces AGPAT1 expression to alter phospholipids and enhance infection in Aedes aegypti. PLoS Pathog 15(12): e1008199. https://doi.org/10.1371/journal.ppat.1008199

Editor: Glenn Randall, The University of Chicago, UNITED STATES

Received: May 30, 2019; Accepted: November 7, 2019; Published: December 9, 2019

Copyright: © 2019 Vial et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: Raw data were deposited in the MassIVE data repository under number MSV000083868; ftp://massive.ucsd.edu/MSV000083868/.

Funding: This work was supported by a grant from National Medical Research Council, Singapore, (NMRC/ZRRF/0007/2017) awarded to JP, a grant from the Ministry of Education, Singapore, (MOE2015-T3-1-003) partially awarded to JP, and by the Duke-NUS Signature Research Programme funded by the Agency for Science, Technology and Research (A*STAR), Singapore, and the Ministry of Health, Singapore. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

Increased global distribution of dengue virus (DENV) is driven by the expansion of its mosquito vectors, mainly Aedes aegypti [1]. An estimated 400 million infections occur yearly in over 100 countries [2] and cause a range of symptoms from flu-like illness to potentially lethal complication called severe dengue. Without approved antiviral drug, treatment is limited to supportive care. Although dengue vaccine is now licensed in several countries, it has variable efficacy against all dengue serotypes [3], and is only suitable for dengue-seropositive patients [4]. To curb dengue epidemics, containment strategy mainly relies on vector control that includes the use of insecticides [5]. However, insecticide resistance is rapidly developing, compromising the efficacy of the only available intervention [6]. Characterization of viral metabolic requirements in mosquitoes will identify targets for novel chemical-based control strategies [7].

As obligate and intracellular parasites, viruses rely on the host to fulfill their metabolite requirements. Glycolysis, amino acid and lipid pathways provide the energy and structural compounds necessary for multiplication. In human cells, DENV alters energy [8,9], glycolysis [10], nucleic acid [11], mitochondrial [12] and lipid metabolisms [13–15]. As an envelope virus, DENV is particularly dependent on host-derived lipid membranes, with which it interacts for entry, replication, translation, assembly and egress [7,16]. In mosquito cells and midgut, DENV reconfigures the lipid profile as indicated by lipidomics [17,18] and transcriptomics [19], particularly altering the membrane lipids such as phospholipids (PL) and sphingolipids. Chemical inhibition of lipid synthesis confirmed the DENV requirements for lipids in human and mosquito cells [20]. Several mechanisms related to elevated autophagy [13] and recruitment of lipogenesis enzymes to the replication complex [21] have been elucidated in the mammalian host. However, how DENV reconfigures lipids in mosquitoes remains unknown.

PL de novo biogenesis is initiated by two types of acyl-transferases that sequentially add two acyls to one glycerol-3-phosphate (G3P) [22,23]. The second addition is catalyzed by 1-acyl-sn-glycerol-3-phosphate O-acyltransferases (AGPAT) that transform lysophospatidate (lysoPA) in phosphatidate (PA) [24]. PAs are then used to produce all PLs, positioning AGPATs as rate-limiting enzymes of PL biogenesis. PLs are produced within the endoplasmic reticulum and subject to swift reconfiguration to meet the cell needs [25,26]. DENV translation, replication and assembly harness the endoplasmic reticulum membranes, suggesting potential alteration of PL biogenesis.

Following an infectious blood meal, DENV infects the mosquito midgut, multiplies, and propagates to the whole mosquito body, including salivary glands, from where it is expectorated during subsequent blood feeding [27]. Throughout various infected tissues, the virus modifies its metabolic environment. In this study, we aim to understand how DENV modifies mosquitoes’ metabolome. Using high-resolution mass spectrometry, we explored the metabolic changes throughout the DENV cycle, in Ae. aegypti Aag2 cell line, midguts and whole mosquitoes at different times post infection. Infection-induced metabolic alteration mostly affected PLs and was associated with AGPAT expression regulations. Combining metabolomic profiling with RNAi-mediated AGPAT depletions, we partially recapitulated the infection-induced PL reconfiguration and demonstrated its pro-viral impact. Eventually, PL profiling upon infection both in vitro and in vivo in PL-altered environments indicated that increased PL redirection/consumption favored virus multiplication. Our study reveals how DENV reconfigures the metabolome and identifies PLs as important components in the virus life-cycle.

Results

Mosquito phospholipidome is reconfigured throughout DENV infection

To describe the metabolic changes in Ae. aegypti, we designed an untargeted multidimensional approach that covered the different stages of DENV cycle (Fig 1A). Extracts from Aag2 cells at 6, 12, 24 and 48 h post infection (hpi) represented changes caused by entry, initiation of replication, replication and virion production (S1A and S1B Fig) [28]. Extracts from midguts at 1 and 7 days post oral infection (dpi) represented infection onset and replication peak, respectively (S1C and S1CD Fig) [27]. Extracts from whole mosquitoes at 1, 7 and 14 dpi represented the different dissemination stages in the mosquito body (S1C and S1CD Fig) [27]. The extracts were analyzed using a Liquid Chromatography-High Resolution Mass Spectrometry (LC-HRMS) metabolomic workflow that detects polar and nonpolar metabolites (S2 Fig). We detected 667, 486 and 1121 compounds in the cells, midgut and whole mosquitoes, respectively (Table 1). Comparison of MS spectra with available non-mosquito databases enabled the annotation of 77% of the peaks. To annotate the remaining mosquito-specific metabolites, we deployed MS spectral similarity network (S3 Fig; S1 Table) [29]. Based on the spectral similarity with database-identified features, we identified the class of 20 MS features without homologues in the databases, increasing the identification coverage to 77, 70 and 82% in the cells, midguts and whole mosquitoes, respectively (Table 1). Among all annotated metabolites, we detected 29% of PLs, 21% of non-PL lipids, 8% of amino acids and peptides, 5% of organic acids, 4% of carbohydrates and 33% of other minor classes (fewer than 3 occurrences). The distribution of metabolite classes indicates that the metabolomic workflow had a broad coverage, although slightly biased towards lipids.

thumbnail
Download:
Fig 1. Aminophospholipid composition is altered throughout DENV infection in mosquitoes.

(A) Schematic of the multidimensional strategy deployed to profile the mosquito metabolome after DENV infection. Aag2 mosquito cells were infected with DENV at an MOI of 5 and collected at 6, 12, 24 and 48 h post infection (hpi). Mosquitoes were orally infected with 107 pfu/ml of DENV. Midguts and whole mosquitoes were collected at 1, 7 and 14 days post infection (dpi). Metabolic extracts were analyzed by LC-HRMS as in S2 Fig. (B) Venn diagrams show regulated metabolites at the different time points and bars indicate class distribution within tissue levels. Triangles indicate direction of regulation and are color-coded with regards to the time point. Uninfected condition was used as control. (C) Fold changes of annotated and significantly regulated metabolites (|log2 fold change| > 1 and p-value < 0.05) by DENV in cells, midguts and mosquitoes as compared to mock infection. Only metabolites from the general metabolism (i.e., lipid, carbohydrate, amino acid and peptide, nucleotide and nucleoside, sialic acid) are shown. *, indicates metabolite annotations determined by spectral similarity. Carbo., carbohydrate; SA, sialic acid; PL, phospholipid; PE, phosphatidylethanolamine; PC, phosphatidylcholine; PS, phosphatidylserine; PA, phosphatidic acid; PI, phosphatidylinositol; PG, phosphatidylglycerol; LysoPS, lysophosphatidylserine; LysoPC, lysophosphatidylcholine; LysoPE, lysophosphatidylethanolamine; LysoPG, lysophosphatidylglycerol; LysoPI, lysophosphatidylinositol; SM, Sphingomyelin; DAG, Diacylglycerol; MAG, Monoacylglycerol; FAHFA, Fatty Acid ester of Hydroxyl Fatty Acid; NAE, N-acylethanolamine; HEA, Heneicosanoic acid; pep., peptides.

https://doi.org/10.1371/journal.ppat.1008199.g001

thumbnail
Download:
Table 1. Summary of metabolites detected across mosquito tissues.

MS Finder (HMDB, ChEBI, LipidMAPS, LipidBlast) score annotation ≥ 5; * p-value < 0.05 as indicated by unpaired t-test and |log2 Fold Change| ≥ 1.

https://doi.org/10.1371/journal.ppat.1008199.t001

Infection regulated 67, 24 and 181 unique metabolites in cells, midguts and whole mosquitoes, respectively (Fig 1B; S2 Table). The PL class had the highest number of regulated metabolites in the midgut and mosquito, and the third most regulated in cells. Fatty acyls were also differentially regulated in the three tissue levels, while amino acids and carbohydrates were regulated in cells and mosquitoes.

Several classes of PLs were regulated in all three tissues (Fig 1C, S4 Fig). Aminophospholipids (aminoPL) are the major constituents of membranes and are synthesized in the endoplasmic reticulum [30]. AminoPLs include phosphatidylethanolamine (PE), phosphatidylcholine (PC) and phosphatidylserine (PS) [31]. Although infection altered different species of aminoPLs in the different tissue-time combinations, we observed a general increase at the beginning of DENV cycle, followed by a reduction at the end. In cells, the majority of regulated aminoPLs (3 out of 5) were upregulated at 48 hpi. In midguts and mosquitoes a total of 9 aminoPLs were increased at 1 dpi, whereas 15 and 14 were reduced at 7 and 14 dpi, respectively.

Phosphatidylglycerols (PG) are another group of PL. Although PGs have a lower abundance than aminoPLs, they are also constituents of membranes and synthesized in the mitochondria [32]. Although only detected in mosquito extracts, PGs responded differently to infection when compared with aminoPLs (Fig 1C). Two shorter PGs decreased at 1 dpi and increased at 14 dpi, and three longer PGs were upregulated at 1 dpi.

Lysophopholipids (lysoPL) are produced from fatty acid remodeling of the different PL classes [33]. In cells, one lysoPG was downregulated at 48 hpi (Fig 1C). In mosquitoes, the lysoPCs followed the same trend as aminoPLs, with 6 increasing at 1 dpi and decreasing at 14 dpi.

Different types of PL precursors were regulated. PA, the direct precursor of all PLs [24,31], was upregulated at 1 dpi in midgut, similar to aminoPLs. Two diacylglycerols (DAG), intermediates for PE and PC productions, were up and downregulated in mosquitoes at 7 and 14 dpi, respectively. Fatty acyls can be incorporated into a glycerol head to produce PL and were regulated. Fatty acyls with an ethanolamine group can be related to PE either via degradation or as a precursor. Similar to aminoPLs, two acyl-ethanolamine increased at 1 dpi and decreased at 14 dpi in mosquitoes. In accordance with previous studies [34], one sphingomyelin was regulated in mosquitoes and one cholesterol compound was depleted in midgut.

Our comprehensive metabolomic profiling revealed that DENV infection profoundly reconfigures the phospholidome. The aminoPLs including PE, PC and PS increased at the beginning of DENV cycle and decreased later on. A large phospholipidome reconfiguration was previously observed in DENV-infected midguts at 3, 7 and 11 dpi [18]. However, in this previous kinetic study, PLs were mostly up-regulated throughout the infection cycle. They reported a decrease in lysoPL abundance, which we observed at 14 dpi in mosquitoes. The authors also noted an increase in sphingolipids, which we saw in mosquitoes but only for one species. Those variations between the studies can stem from methodological differences in extraction and analyses, or from biological differences in mosquito colony and virus strain. Taken together, our results and those of others indicate that DENV infection reconfigures lipid membrane composition.

DENV infection modulates expression of AGPAT1 that is involved in PL biogenesis

To determine how PLs are reorganized upon DENV infection, we first characterized the rate-limiting AGPAT enzymes in Ae. aegypti (Fig 2A). In humans, there are five AGPAT isoforms with different activities depending on four motif sequences [35]. Motifs I and IV bind to acyl-CoA and catalyze lysoPA to PA acylation, while motifs II and III bind to lysoPA. Human AGPAT (hAGPAT) 1 and 2 are localized in the endoplasmic reticulum and have the highest acyltransferase activity and lysoPA affinity [36]. hAGPAT3-5 have lower transferase activity and target different substrates that lysoPA such as lysoPL, participating in PL remodeling [37]. Aedes aegypti also has five AGPAT isoforms: AGPAT1 (AAEL011898), AGPAT2 (AAEL001000), AGPAT3 (AAEL011902), AGPAT4 (AAEL014026) and AGPAT5 (AAEL011901). Based on amino acid similarity, AGPAT1, 2, 3 and 5 cluster with hAGPAT1 and 2, whereas AGPAT4 clusters with hAGPAT3-5 (Fig 2B; S5 Table). Further classification based on functional motifs that are define the biochemical activity [38] establishes AGPAT1, 2, 3 as homologues of hAGPAT1 and 2 (Table 2), suggesting that they all share the acyltransferase activity and lysoPA affinity.

thumbnail
Download:
Fig 2. DENV infection decreases AGPAT1 but not AGPAT2 expression.

(A) An overview of phospholipid biogenesis. Sequential additions of two acyl-coA to one G3P produce a LysoPA and a PA. Several AGPATs mediate the second addition. PA forms either DAG or CDP-DAG, each generating a different set of phospholipids. DAG produces aminophospholipids by addition of a Cho or Eth group. Aminophospholipids are also produced by acylation of lysoPL and base group modifications. CDP-DAG produces PI and PG, the latter being transformed in CL by combination with a second PG. Phospholipases cleave off an acyl chain from PL to produce lysoPL. G3P, glycerol-3-phosphate; LysoPA, lysophosphatidate; PA, phosphatidate; AGPAT, acyl-sn-glycerol-3-phosphate acyltransferases; DAG, diacylglycerol; CDP-DAG, cytidine diphosphate diacylglycerol; PC, phosphatidylcholine; LPC, lysophosphatidylcholine; Cho, choline; PE, phosphatidylethanolamine; LPE, lysophosphatidylethanolamine; Eth, ethanolamine; PS, phosphatidylserine; LPS, lysophosphatidylserine; PI, phosphatidylinositol; LPI, lysophosphatidylinositol; PG, phosphatidylglycerol; LPG, lysophosphatidylglycerol; CL, cardiolipin; LCL, lysocardiolipin. (B) Maximum likelihood tree between AGPATs from Ae. aegypti and humans. (C) AGPAT1 and AGPAT2 expressions in DENV-infected cells at 24 and 48 hpi with an MOI of 5. Expressions of AGPAT1 and AGPAT2 (D) in midguts at 1, 7 dpi and (E) in whole mosquitoes at 1, 7 and 14 dpi with DENV at 107 pfu/ml. (C-E) Actin expression was used for normalization. Bars show means ± s.e.m from 4 independent wells or 3 pools of 5 midguts or 5 mosquitoes. *, p-value < 0.05; **, p-value<0.01; as indicated by unpaired t-test.

https://doi.org/10.1371/journal.ppat.1008199.g002

thumbnail
Download:
Table 2. Acyltransferase motif comparison between human and Ae. aegypti AGPAT homologues.

Purple indicates amino acid residues highly conserved across all hAGPATs and yellow in hAGPAT1 and 2 only.

https://doi.org/10.1371/journal.ppat.1008199.t002

To test whether infection-induced phospholipidome reconfiguration is associated with AGAPT regulation, we quantified AGPAT1 and 2 expressions. Interestingly, we observed that AGPAT1 was downregulated at 24 and 48 hpi in cells, at 1 and 7 dpi in midguts and at 1, 7 and 14 dpi in mosquitoes (Fig 2C–2E). In contrast, AGPAT2 was not significantly regulated in cells, midguts and mosquitoes (Fig 2C–2E). This partially corroborates a transcriptomic study [39] that shows AGPAT1 downregulation and AGPAT2 upregulation at 1, 2 and 7 days post DENV-inoculation in mosquitoes (S5 Fig). To test whether virus replication was required for AGPAT1 down-regulation, we incubated cells with UV-inactivated virus (S6A Fig). AGAPT1 expression did not vary between mock and UV-inactivated DENV at 24 and 48 hpi (S6B Fig), indicating that active infection is required.

Altogether, our results and those of others suggest that AGPAT1 down-regulation correlates with DENV phospholipidome reconfiguration.

Depletion of AGPAT1 but not AGPAT2 promotes DENV infection by increasing aminoPL concentrations in cells

Based on the association between the infection-induced phospholipidome reconfiguration and AGPAT1 down-regulation, we hypothesized that AGPAT1 mediates the phospholipidome reconfiguration that promotes DENV infection. First, we described how both AGPAT1 and 2, the latter used as control, regulate PL biogenesis in non-infected mosquito cells. LC-HRMS polar mode detection was used to target phospholipid metabolites (S2 Fig). Depletion of either AGPAT1 or 2 increased the concentrations of aminoPLs (Fig 3A–3C; S6 Table), as previously observed in human cells [40]. Strikingly, PC (34:1), PC (38:5) and G3P were upregulated by both AGPAT1 and 2 depletion. AGPAT2 depletion also increased two other PCs and decreased two forms of sphinganine, a precursor of sphingolipid biosynthesis [41] (Fig 3D; S7A Fig), and six amino acids or nucleoside (arginine, cysteine, methionine, glutamic acid, oxidized glutathione and adenosine) (S6 Table). Of note, depletion of one of the AGPATs did not alter expression of the other (S8 Fig), indicating the enzyme-specificity of the aminoPL alterations. These results confirm the roles of AGPAT1 and 2 in PL biogenesis, and reveal differences between the two enzymes.

thumbnail
Download:
Fig 3. AGPAT1 but not AGPAT2 depletion increases DENV multiplication and aminoPL in cells.

Aag2 cells were transfected with either dsRNA against AGPAT1 or 2 (dsAGPAT1 or 2) or with dsRNA control (dsControl). At 24 h post transfection, cells were infected with DENV at MOI of 1 or mock infected. Supernatant was collected at 48 hpi. Expressions of (A) AGPAT1 and (B) AGPAT2 in mock-infected cells at 72 h post transfection. Actin expression was used for normalization. Bars show mean ± s.e.m. (C) Number of phospholipid-related metabolites significantly regulated. (D) Impact of AGPAT1 or 2 depletion on the lipidome of mock-infected cells at 72 h post transfection. *, p-value <0.05 and |log2 fold change| > 1. (E) Impact of AGPAT1 or 2 depletion on DENV production at 48 hpi as determined by plaque forming unit (pfu) assay. Bars show mean ± s.e.m. and each point represents independent wells. (F) Impact of AGPAT1 or 2 depletion on the lipidome of infected cells at 48 hpi. *, p-value < 0.05 and |log2 fold change| > 1. (A, B, D and F) result from three biological replicates. (A, B) ***; p-value < 0.001; ****, p-value < 0.0001, as indicated by unpaired t-test. (E) **, p-value < 0.01, as indicated by Dunnett’s test. PE, phosphatidylethanolamine; PC, phosphatidylcholine; PS, phosphatidylserine; LPE, lysophosphatidylethanolamine.

https://doi.org/10.1371/journal.ppat.1008199.g003

To test whether AGPAT-mediated reconfiguration of aminoPLs promotes DENV infection, we quantified infection in cells depleted of either AGPAT1 or 2. While DENV titer was not altered by AGPAT2 depletion, it increased 2.09 ± 0.36 fold (p-value = 0.0012) following AGPAT1 depletion (Fig 3E). Next, we described how DENV infection altered the lipidome in cells depleted of either AGPAT1 or 2. Infection in AGPAT1-depleted cells, but not in AGPAT2-depleted cells, mostly increased aminoPLs as compared to infected wild-type cells (Fig 3C, S7 Table). Specifically, two PEs, two PCs, two PSs, and one lysoPE were upregulated (Fig 3F; S7B Fig). These results indicate that: (i) depletion of AGPAT1 prior infection amplifies the aminoPL increase observed in early mosquito DENV infection (Fig 1C), and (ii) AGPAT1-mediated aminoPL reconfiguration is associated with increased DENV production. Therefore, AGPAT1 downregulation by infection generates a pro-viral environment.

To test whether AGPAT1 effect on DENV was related to aminoPLs, we modified the AGPAT1-induced reconfiguration of aminoPLs by media supplementation. In PL biogenesis, PA produces DAG, which is transformed into PE by addition of ethanolamine (Fig 2A). Extracellular source of ethanolamine influences phospholipid metabolism [42], especially PEs that were altered by AGPAT1 depletion (Fig 3F). Therefore, we measured the impact on DENV gRNA of ethanolamine supplementation upon AGPAT1 depletion. Controls were non-depleted cells without supplementation (standard media that does not contain ethanolamine), non-depleted cells with ethanolamine supplementation and AGPAT1-depleted cells without supplementation. While ethanolamine supplementation did not alter gRNA in non-depleted cells, the increase observed upon AGPAT1-depletion was reverted to non-depleted non-supplemented control when ethanolamine was supplemented (S9 Fig). These results confirm the role of AminoPLs in AGPAT1 increase of DENV multiplication and incriminate metabolites downstream of DAG as important for DENV.

AGPAT1 depletion promotes DENV infection by amplifying aminoPL reconfiguration in mosquitoes

To further characterize how AGPAT1-mediated PL alteration increases DENV infection, we depleted AGPAT1 in mosquitoes (Fig 4A). We first described how AGAPT1 regulates PLs in mosquitoes that fed on a non-infectious blood meal (Fig 4B; S10 Fig; S8 Table). As in cells, AGPAT1 mostly regulated aminoPLs and their derivatives (S8 Table). At 2 days post oral feeding, one DAG and one PA, two key precursors of PLs, and one PS were increased. At 7 days post oral feeding, eight lysoPLs (one lysoPI, two lysoPGs, one lysoPS, two lysoPEs and two lysoPCs) were downregulated, similarly to those in DENV-infected wild-type mosquitoes at 14 dpi (Fig 1C). Strikingly, some metabolites were similarly regulated by either AGPAT1-depletion or DENV infection. PS (34:2) was upregulated by either AGPAT1-depletion at 2 days post oral feeding or DENV infection at 1 dpi. LysoPI (20:4), lysoPE (18:1), lysoPC (16:1) and lysoPC (20:4) were reduced by either AGPAT1-depletion at 7 days post oral feeding or DENV infection at 14 dpi. At 7 days post oral feeding in AGPAT1-depleted mosquitoes, one fatty acyl, two acylcarnitines, oleoylcarnitine (18:1) and acetyl-L-carnitine, were decreased. Of note, one of the regulated acylcarnitines (i.e., oleoylcarnitine) was inversely regulated by DENV infection at 14 dpi (Fig 1C). These results show that AGPAT1 depletion partially reproduces the aminoPL reconfiguration caused by DENV infection.

thumbnail
Download:
Fig 4. AGPAT1 depletion increases DENV multiplication and consumption of aminoPLs in mosquitoes.

Mosquitoes were injected with either dsRNA against AGPAT1 (dsAGPAT1) or dsRNA control (dsControl). Two days later, mosquitoes were orally fed with non-infectious blood or infectious blood containing DENV at 107 pfu/ml. (A) Validation of AGPAT1 silencing in non-infected mosquitoes at two days post dsRNA injection. Actin expression was used for normalization. Bars show mean ± s.e.m. from 3 pools of 5 mosquitoes each. ****, p-value < 0.0001 as indicated by unpaired t-test. (B) Impact of AGPAT1 depletion on the lipidome of mosquitoes at 2 and 7 days post non-infectious blood feeding (dpb). Impact of AGPAT1 depletion on (C) DENV gRNA copies and (D) viral load measured as pfu/ml at 7 days post oral infection (dpi). Bars indicate mean (C) or geometric means ± 95% CI (D) with each dot representing one mosquito. ***, p-value < 0.001 ****, p-value < 0.0001 as indicated by unpaired t-test (C) or Mann-Whitney test (D). (E) Impact of AGPAT1 depletion on the lipidome of infected mosquitoes at 7 dpi. (F) Impact of infection on the lipidome of AGPAT1-depleted mosquitoes at 7 dpi. (B, E, F) Only regulated metabolites are shown (p-value < 0.05 and |log2 fold change| > 1). PL, phospholipid; FA; fatty acyl; AC, acylcarnitine; DAG, diacylglycerol; PA, phosphatidic acid; PS, phosphatidylserine; LPE, lysophosphatidylethanolamine; LPC, lysophosphatidylcholine; LPS, lysophosphatidylserine; LPG, lysophosphatidylglycerol; LPI, lysophosphatidylinositol; PE, phosphatidylethanolamine; PC, phosphatidylcholine; PI, phosphatidylinositol; GPEtn, glycerophosphoethanolamine.

https://doi.org/10.1371/journal.ppat.1008199.g004

Next, we found that AGPAT1-depletion increased DENV gRNA (p-value = 0.0004) and titer (p-value < 0.0001) at 7 dpi in whole mosquito (Fig 4C and 4D), confirming that AGPAT1 depletion induces a pro-viral environment in vivo. However, the pro-viral impact of AGPAT1 depletion was not observed on infection rate (S11A Fig) and at 2 dpi on gRNA (S11B Fig). We repeated the experiment using a lower inoculum for oral infection and similarly observed no effect at 2 dpi but a moderate increase of gRNA (p = 0.0421) at 7 dpi (S11C–S11E Fig). It is intriguing that AGPAT1 depletion only increased infection at 7 dpi. This may indicate that DENV infection is necessary to amplify AGPAT1-mediated alteration in establishing a pro-viral environment. To characterize the AGPAT1-induced proviral environment, we compared the lipidome of infected mosquitoes that were AGPAT1 depleted or not. AGPAT1-depletion increased aminoPL concentrations at 7 dpi (Fig 4E; S9 Table). Specifically, one PC, four PEs and one glycerophosphoethanolamine were upregulated. By artificially modulating the metabolome, we revealed that AGPAT1 depletion favors DENV by increasing aminoPL concentrations.

When profiling metabolic changes throughout mosquito DENV cycle, we observed that aminoPLs were reduced at the end of DENV cycle (Fig 1C). To test whether this reduction occurred when aminoPL concentrations were increased by AGPAT1 depletion (Fig 4E), we examined how DENV modifies the lipidome in AGPAT1-depleted mosquitoes. Control mosquitoes were depleted of AGPAT1 and fed a non-infectious blood meal. At 7 dpi, aminoPLs were drastically decreased by DENV despite AGPAT1-depletion enhancement (Fig 4F; S10 Table). Specifically, three PEs and five PCs were reduced more than 100 folds. These results suggest that, by reproducing infection-induced aminoPL reconfiguration, AGPAT1 depletion amplifies the DENV reconfiguration, thereby promoting infection. Altogether, our results indicate that DENV inhibits AGPAT1 expression to increase the amplitude of aminoPLs consumption/redirection for virus multiplication.

Discussion

By combining observational and manipulative approaches in vitro and in vivo, we deciphered how DENV hijacks the mosquitoes’ metabolome. We deployed metabolomic profiling throughout the mosquito DENV cycle, and revealed an overall increase in aminoPLs, followed by a reduction at the end of the cycle. We next showed that aminoPL reconfiguration is partially mediated by infection-induced AGPAT1 down-regulation, which increases aminoPL concentrations. Because AGPAT1 depletion promotes virus multiplication, our study discovers a mechanism whereby DENV reconfigures the metabolome to its benefit. Furthermore, we show that in an environment richer in aminoPLs (AGPAT1 depletion), infection reduces aminoPL concentrations at a higher and earlier rate than in a wild-type organism, while increasing virus production. This suggests that DENV consumes/redirects aminoPLs, although this is based on correlations between indirect lipidome alteration and virus multiplication. Altogether, we propose a model whereby DENV regulates aminoPL enzymes to increase aminoPL concentrations early during the cycle and consumes/redirects them for its multiplication.

We revealed that DENV-induced phospholipidome reconfiguration is partially mediated by lowering AGPAT1 expression. Indeed, AGPAT1 was downregulated upon infection and its depletion partially recapitulated DENV-induced aminoPL reconfiguration both in vitro and in vivo. That we only partially recapitulated the DENV-induced phospholipidome reconfiguration suggests that other enzymes, such as those directly producing and those hydrolyzing aminophospholipids [43], play a role. AGPAT enzymes can influence PL composition at several levels. By catalyzing the second acylation that produces PA (Fig 2A), AGPAT enzymes are responsible for a rate-limiting step in PL biogenesis [35]. However, we did not observe PA reduction upon AGPAT1 depletion and instead reported an alteration of aminoPLs. This pattern was previously observed in hAGPAT-depleted human cells [40] and was attributed to compensation by other AGPAT isoforms. Lack of AGPAT1 may unbalance substrate competition with other isoforms with different lysoPA and acyl affinity [44]. Consequently, this can influence PA structure and subsequent PL composition. In addition, certain AGPAT isoforms can transfer acyls to lysoPL, thus, participating in PL remodeling independently of de novo synthesis [45]. In mosquito cells, AGPAT2 depletion also altered PL profile. However, AGPAT2 expression is not reduced by DENV infection and its depletion does not impact virus multiplication. The correlation between expression of AGPAT isoforms and their impacts on infection indicates a fine regulation of metabolism by DENV, and identifies a new mechanism by which the virus hijacks the host phospholipidome.

The expression of AGPATs is regulated either directly through their transcription factors or indirectly by altering the endoplasmic reticulum topology. Sterol regulatory element binding proteins (SREBP) coordinate fatty acid, sterol and phospholipid metabolisms by transcriptionally regulating enzymes, such as acyl transferases [46]. DENV protein interaction with SREBP, as described with other transcription factors [47], could alter AGPAT expressions. In support of this, chemical inhibition of SREBP blocks DENV replication [48]. Alternatively, AGPAT activity is influenced by substrate accessibility, which depends on membrane topology [36]. Several DENV proteins are embedded in the endoplasmic reticulum membrane and modify its topology [49,50]. This can alter endoplasmic reticulum-located AGPAT activity, PL profile and activate acyl-transferase expressions to restore homeostasis [51].

DENV intricately interacts with aminoPL-containing membranes and a change as we report can alter several stages of its life-cycle, such as entry, replication, translation, assembly and egress [7,16]. DENV interacts with the plasma membrane for entry and fusion, and with the endoplasmic reticulum membrane for translation, replication and assembly [28]. These membranes are mostly composed of PE and PC, and PS in lower proportions [52]. Their fluidity and topology are determined by the aminoPL structure and concentration. The cylindrical shape of PCs stabilizes lipid bilayers, whereas the conical shape of PEs induces curvature. Several flavivirus-aminoPL interactions have been described in human cells. PLs present in the DENV envelope are recognized by cellular ligands and mediate entry [53,54]. Replication of flaviviruses induce the invagination of the endoplasmic reticulum by altering PE [55,56] and lysoPC compositions [43]. Alternatively, AGPAT expressions could alter anti-viral immune response [57]. In mammalian cells, AGPAT overexpressions amplify cellular signaling of cytokine [58] that can reduce DENV infection in mosquito cells [59]. Although the precise function of AGPAT1-regulated aminoPLs is unknown and may be multifactorial, previous studies and ours indicate that membrane aminoPL reconfiguration influences DENV multiplication.

In conclusion, our study determines the importance of aminoPL reconfiguration for DENV infection. We also reveal the underlying mechanism for viral phospholipidome reconfiguration. The intricate metabolic interactions between DENV and mosquitoes represents a target to control transmission.

Materials and methods

Aedes aegypti mosquitoes and cell line

The Aedes aegypti colony was established in 2010 from Singapore and was reared at 28°C and 60% relative humidity with 12h:12h light:dark cycle. Eggs hatched in milliQ water were fed with a mix of fish food (TetraMin fish flakes), yeast and liver powder (MP Biomedicals). Adults were held in rearing cages (Bioquip) supplemented with water and 10% sucrose solution. Aedes aegypti Aag2 cells [60] were grown in RPMI-1640 medium (Gibco) with 10% filtered fetal bovine serum (FBS) (Hyclone) and 1% Penicillin-Streptomycin (Gibco). For media supplementation, 2 mM of ethanolamine (Sigma) was added to Aag2 growth medium. Cells were maintained in vented culture flasks in a humidified incubator with 5% CO2 at 28°C. BHK-21 (baby hamster kidney) (ATCC CCL-10) cells were grown in the same media and maintained at 37°C with 5% CO2.

Dengue virus

Dengue virus serotype-2 strain ST (DENV) was collected from the Singapore General Hospital in 1997 [61]. DENV was propagated alternatively in Vero (ATCC CCL-81) and C6/36 (ATCC CRL-1660) cells. Virus titer was determined by plaque assay using BHK-21 cells. DENV supernatant was exposed to UV light from biosafety cabinet (Sterilgard III advance, the Baker Company) for 1h at room temperature for inactivation.

Oral infection of mosquitoes

Two- to four-day-old adult female mosquitoes were starved for 24h before oral feeding on a blood meal containing 40% volume of washed erythrocytes from SPF pig’s blood (PWG Genetics, Singapore), 5% of 100 mM ATP (Thermo Fisher Scientific), 5% of human serum (Sigma-Aldrich) and 50% of DENV-2 in RPMI media (Gibco). The virus titer in the blood meal was 2 x 107 pfu/ml and validated by plaque assay. Blood was maintained at 37°C using hemotek membrane feeder system (Discovery Workshops) with sausage casing for 1.5 h. A control group was allowed to feed on the same mix of SPF pig blood meal without virus. Engorged mosquitoes were visually selected and maintained at 28°C with water and 10% sucrose solution.

Cell inoculation

5 x 106 cells were inoculated with DENV at an MOI of 5 in serum-free RPMI media for 1h. The inoculum was then replaced with 2% FBS RPMI media. Mock infection was used as negative control.

Metabolite extraction from mosquitoes, midguts and cells

At 1, 7 and 14 days post-oral feeding, 10 mosquitoes in 500 μl of ice-cold methanol and water ratio of 80:20 (LCMS grade, Thermo Fisher) were homogenized with bead Mill homogenizer (FastPrep-24, MP Biomedicals) and sonicated for 15 min in an ultrasonic bath (J.R. Selecta) at 4°C. Homogenates were centrifuged at 10,000 rpm for 1 min at 4°C to collect 400 μl supernatant. Pellets were further extracted twice by addition of 500 μl of methanol:water (80:20), followed by centrifugation. Supernatants were combined and vacuum-dried (Speed-Vac, Thermo-Scientific) before storage at -20°C. At 1 and 7 days post-oral feeding, 10 midguts were homogenized in 200 μl of the methanol:water (80:20) solution and 120 μl of supernatant was collected three times with the same protocol. At 6, 12, 24 and 48 h post-inoculation, cells were washed with room temperature 0.9% NaCl and collected in 2 ml of ice-cold methanol:water (80:20) by scraping. Cells were homogenized by ultrasound and extracted thrice as detailed above by adding 500 μl of ice-cold methanol:water (80:20). Three biological replicates were conducted per condition.

LC-HRMS metabolism profiling

Dry extracts were normalized at 2 mg/mL in methanol:water 80:20 solution and metabolites were detected using two methods. Compounds from medium range polarity to lipophilic substances were detected using a UPLC-UV-QTOF-MSE instrument (Xevo G2 QTof, Waters) mounted on an electrospray ionization (ESI) source, with a UPLC BEH C18 Acquity column (100 × 2.1 mm i.d., 1.7 μm, Waters) equipped with a guard column. Mobile phase A was of 0.1% formic acid in water, B was 0.1% formic acid in acetonitrile. The flow rate was 400 μl/min and the gradient ran with 98% A for 0.5 min to 20% B over 3.5 min, 98% B for 8 min, held at 98% B for 3 min, and returned in 0.5 min to initial conditions (98% A), finally held for 3.5 min to assure equilibration before the subsequent analysis. Detection was performed at 254 nm by TOF-MS in both electrospray (ESI) negative mode with voltage at 2.5 kV and positive mode with voltage at 3.0 kV. The m/z range was 100–1200 Da with a scan time of 0.1 s. All detected ions were fragmented using MSE scan with an energy collision ramp from 20 to 50 eV. All analyses were acquired using leucine enkephalin as the lock mass at a concentration of 400 pg/μl and flow rate 7 μl/min. The injection volume was 2 μl and samples were kept at 10°C during the whole analysis.

Polar metabolites were profiled using a UPLC-LTQ Orbitrap XL instrument (Ultimate 3000, Thermo Fisher Scientific, Hemel Hempstead, UK) set at 15,000 resolution, with a Zic-pHilic column (150 × 2.1 mm i.d., 5 μm, SeQuant, Merck). Mobile phase A was 20 mM ammonium acetate buffered at pH 9 and B was acetonitrile. The flow rate was 250 μl/min and the gradient ran from 90% B for 0.5 min to 40% B over 18 min, held at 40% B for a further 3 min, and then returned in 0.5 min to initial conditions (90% B) finally held for 5 min before subsequent analysis. The m/z range was 100–1500 and ISpray voltage at 4.2 kV (positive mode) and 3.0 kV (negative mode). Each full MS scan was followed by data dependent MS/MS on the two most intense ions using stepped CID fragmentation mode at 35% normalized collision energy, isolation width of 2 u and activation Q set at 0.250.

Data analysis and visualization

Peak detection and alignment were performed using MS-DIAL (ver. 3.12) [62]. Peak annotation was done using MS-finder (ver. 3.04) [63] with HMDB, ChEBI, LipidMAPS and LipidBlast databases, allowing a level 2.2 of metabolite identification [64,65]. Data were normalized by total ion chromatogram (TIC) and features lower than 2-fold average blank were removed. Each LC-HRMS condition was analyzed separately before concatenation. Data were normalized by auto-scaling before selecting regulated metabolites with more than 2-fold intensity change and a p-value < 0.05 with FDR adjustment as indicate by unpaired t-test using MetaboAnalyst (ver. 4.0) [66]. To account for physiological variations between the different time-tissue combinations, the t-tests were done by comparing the same tissue in infected and uninfected conditions within each time. Uninfected condition was used as control. PCA for quality control was performed with MetaboAnalyst (ver. 4.0). MS/MS similarity metabolic networks with cut-off > 60% were generated with MS-Finder for each LC-HRMS mode-tissue combination.

dsRNA-mediated RNAi

Templates for dsRNA against AAEL011898 and AAEL001000 were PCR amplified with primers flanked with a T7 promoter (S3 Table) from mosquito cDNA. dsRNA was synthetized with megaScript T7 transcription kit (Thermo Fisher Scientific), extracted in DEPC-treated water and annealed by slowly cooling down from 95°C. Control dsRNA targeting LacZ was produced [67]. Two- to five-day-old cold-anesthetized female mosquitoes were intrathoracically injected with 69 nl of 3 mg/ml of dsRNA by using Nanoject II injector (Drummond Scientific). Mosquitoes were then maintained at 28°C with water and 10% sucrose solution before oral infection as detailed above. Cells were seeded at 2 x 105 per 24-well plate and transfected after 24h with 1 μg of dsRNA by using TransIT-mRNA Transfection kit (Mirusbio). Infection with DENV was done one day post transfection.

Quantification of DENV genomic RNA

Single mosquitoes or tissues were homogenized in 350μl of TRK lysis buffer (Omega Bio-tek) using a bead Mill homogenizer (FastPrep-24, MP Biomedicals). Total RNA was extracted using E.Z.N.A. Total RNA kit I (Omega Bio-tek) and eluted in 30μl of DEPC-treated water. Genomic RNA (gRNA) was quantified with one-step RT-qPCR using iTaq Universal probe kit (Bio-Rad) and primers and probes targeting the DENV Envelope [68]. The 12.5 μl reaction mix contained 1 μM of forward and reverse primers, 0.125 μM of probe and 4 μl of RNA extract. Quantification was conducted on a CFX96 Touch Real-Time PCR Detection System (Bio-Rad). Thermal profile was 50˚C for 10 min, 95˚C for 1 min and 40 cycles of 95˚C for 10 sec and 60˚C for 15 sec.

An absolute standard curve was generated by amplifying fragments containing the qPCR target using a forward primer tagged with T7 promoter; forward: 5’-CAGGATAAGAGGTTCGTCTG-3’ and reverse: 5’-TTGACTCTTGTTTATCCGCT-3’, resulting in a 453bp fragment. The fragment was reverse transcribed using MegaScript T7 transcription kit (Ambion) and purified using E.Z.N. A. Total RNA kit I. The total amount of RNA was quantified using a Nanodrop (Thermo Fisher Scientific) to estimate copy number. Ten times serial dilutions were made and used to generate absolute standard equation for gRNA. In each subsequent RT-qPCR plate, five standards were added to adjust for threshold variation between plates. The infection rate was calculated by dividing the number of samples with detectable gRNA over total number of samples.

Titration

Titration was conducted by plaque assay with BHK-21 cells as described previously [69]. Briefly, 80–90% confluent cells were inoculated with serial 10-fold dilutions of samples for 1h. Cells were then incubated with 1% carboxyl-methyl cellulose (CMC) (Merck) for 5 days, fixed with 4% formaldehyde (Merck) -PBS and stained with 1% crystal violet (Sigma-Aldrich) solution to count plaque forming units (pfu).

Quantification of gene expression

Total RNA from five mosquitoes was extracted using E.Z.N.A. Total RNA kit I, treated with RapidOut DNA Removal kit (Thermo Fisher Scientific) and reverse transcribed with iScript cDNA Synthesis Kit (Bio-Rad). Gene expression was quantified using iTaq Universal SYBR Green Supermix (Bio-Rad) and primers detailed in S4 Table. Actin expression was used for normalization. Quantification was conducted in a CFX96 Touch Real-Time PCR Detection System (Bio-Rad). Thermal profile was 95°C for 1 min and 40 cycles of 95°C for 10 sec and 60°C for 15 sec. Three biological replicates were conducted.

Statistical analysis

Differences in gRNA copies per infected mosquito were tested on values using unpaired t-test or Mann-Whitney test depending on the normal distribution estimated with D’agostino and Pearson normality test. Differences in percentages were tested using χ2 test. Tests were performed with GraphPad PRISM software (ver. 6.01).

AGPAT sequence alignment

Amino acid sequence homology was determined by MEGA X software (ver. 10.0.5), using Maximum Likelihood and bootstrapping. FASTA sequences (S5 Table) were retrieved from ncbi.nlm.nih.gov.

Supporting information

S1 Fig. Quantification of DENV infection in Aag2 cells, A. aegypti mosquito and midgut.

(A) DENV gRNA copies in Aag2 cells at 6, 12, 24 and 48h post-infection with DENV at MOI = 5. Points from 4 repeat and standard errors show geometric mean ± 95% CI. (B) Plaque titer (plaque forming unit—pfu) in supernatant from Aag2 cells at 6, 12, 24 and 48h post-infection. Aag2 cells were infected with DENV at MOI = 5 and virus titer was calculated using plaque assay. Each point represents one well. (C) DENV gRNA copies per infected mosquitoes and dissected midguts at 1, 7 and 14 days post-oral infection with 107 pfu/ml. Each point represents one mosquito or midgut. Bars show geometric mean ± 95% CI. (D) Infection rate in whole mosquito and midgut at 1, 7 and 14 days post-oral infection. Bars represent percentages ± s.e.

https://doi.org/10.1371/journal.ppat.1008199.s001

(TIF)

S2 Fig. LC-HMRS analytical pipeline.

Liquid Chromatography-High Resolution Mass Spectrometry (LC-HRMS) pipeline used to detect polar and nonpolar metabolites.

https://doi.org/10.1371/journal.ppat.1008199.s002

(TIF)

S3 Fig. Spectral similarity network from mosquito MS features.

Example of a molecular spectral network for Ae. aegypti mosquito at 14 days post-infection using MS features detected with the non-polar LC condition and MS negative mode. Line length represents the MS/MS score similarity. Ontology for unknown features was determined based on the proximity with database-identified features.

https://doi.org/10.1371/journal.ppat.1008199.s003

(TIF)

S4 Fig. Ion intensity of regulated metabolites in cells, midguts and mosquitoes infected with DENV and mock.

Normalized ion intensity was calculated after total ion chromatography normalization and auto scaling from three biological replicates. Conditions with significantly regulated metabolites (p-value <0.05 and |log2 fold change| >1) were indicated with an asterisk. Only metabolites from the general metabolism (i.e., lipid, carbohydrate, amino acid and peptide, nucleotide and nucleoside, sialic acid) are shown. †, indicates metabolite annotated by spectral similarity. Carbo., carbohydrate; SA, sialic acid; PL, phospholipid; PE, phosphatidylethanolamine; PC, phosphatidylcholine; PS, phosphatidylserine; PA, phosphatidic acid; PI, phosphatidylinositol; PG, phosphatidylglycerol; LysoPS, lysophosphatidylserine; LysoPC, lysophosphatidylcholine; LysoPE, lysophosphatidylethanolamine; LysoPG, lysophosphatidylglycerol; LysoPI, lysophosphatidylinositol; SM, Sphingomyelin; DAG, Diacylglycerol; MAG, Monoacylglycerol; FAHFA, Fatty Acid ester of Hydroxyl Fatty Acid; NAE, N-acylethanolamine; HEA, Heneicosanoic acid; pep., peptides.

https://doi.org/10.1371/journal.ppat.1008199.s004

(TIF)

S5 Fig. AGPAT genes are regulated by DENV infection, from Colpitts transcriptomic data.

Ae. aegypti female mosquitoes were inoculated with DENV serotype 2 and collected at 1, 2 and 7 days post inoculation for transcriptomic analysis using microarray. Data from 3 separate infections. AGAPAT 1–5 were found significantly regulated by DENV infection. Data retrieved from (Colpitts et al., 2011, PMID: 21909258).

https://doi.org/10.1371/journal.ppat.1008199.s005

(TIF)

S6 Fig. UV-inactivated DENV does not regulate AGPAT1 expression.

Aag2 cells were infected with an MOI of 5 of DENV (DENV-WT), UV-inactivated DENV (DENV-UV) or mock. Cells were analyzed at 24 and 48 hours post-infection (hpi). (A) DENV gRNA copies. Bars show geometric means ± 95% C.I. (B) AGPAT1 expression relative to Actin level. Bars show arithmetic means ± s.e.m. (A-B) Each point represents an independent well. ***, p-value < 0.001 as determined by unpaired t-test.

https://doi.org/10.1371/journal.ppat.1008199.s006

(TIF)

S7 Fig. Metabolomic impact of AGPAT1 and AGPAT2 depletion and infection in cells as measured by ion intensity.

Aag2 cells were transfected with dsRNA against AGPAT1 (dsAGPAT1) or AGPAT2 (dsAGPAT2) or a dsRNA control (dsControl). 24h later, cells were infected with DENV at MOI of 1. (A) Ion intensity of regulated metabolites in mock cells at 72h post transfection. (B) Ion intensity of regulated metabolites in infected cells at 24 hpi. Normalized ion intensity was calculated after total ion chromatography normalization and auto scaling from three replicates. Conditions with significantly regulated metabolites (p-value <0.05 and |log2 fold change| >1) were indicated with an asterisk. PE, phosphatidylethanolamine; PC, phosphatidylcholine; PS, phosphatidylserine; LPE, lysophosphatidylethanolamine.

https://doi.org/10.1371/journal.ppat.1008199.s007

(TIF)

S8 Fig. AGPAT1 and 2 expression in cells after the other AGPAT depletion.

Aag2 cells were transfected with dsRNA against AGPAT1 or 2 (dsAGPAT1 or 2). Control cells were transfected with dsRNA control (dsControl). Cells were collected 72h post dsRNA. (A) AGPAT1 expression in AGPAT2-depleted cells. (B) AGPAT2 expression in AGPAT1-depleted cells. Bars show mean ± s.e.m from 3 biological replicates. ns, non-significant, as indicated by unpaired t-test.

https://doi.org/10.1371/journal.ppat.1008199.s008

(TIF)

S9 Fig. Ethanolamine supplementation partially rescued infection increase upon AGPAT1 depletion.

24h before infection, Aag2 cells were transfected with dsRNA against AGPAT1 (dsAGPAT1) or with dsRNA control (dsControl) and reared in standard growth media or the same media supplemented with 2mM ethanolamine. Cells were infected with DENV at MOI of 1 and gRNA copy was quantified 48h later. Bars show geometric means ± 95% C.I. Each point represents an independent well. *, p-value < 0.05; **, p-value < 0.01 as determined by unpaired t-test.

https://doi.org/10.1371/journal.ppat.1008199.s009

(TIF)

S10 Fig. Metabolomic impact of AGPAT1 and AGPAT2 depletion in uninfected mosquitoes as measured by ion intensity.

Two days post dsRNA injection against AGPAT1 (dsAGPAT1) or control (dsControl), mosquitoes were orally infected with DENV at 107 pfu/ml. Metabolomic analyses were performed at 2 and 7 dpi. Normalized ion intensity was calculated after total ion chromatography normalization and auto scaling from three replicates. Conditions with significantly regulated metabolites (p-value <0.05 and |log2 fold change| >1) were indicated with an asterisk.

https://doi.org/10.1371/journal.ppat.1008199.s010

(TIF)

S11 Fig. Impact of AGPAT1-depletion in mosquitoes on DENV infection rate and gRNA copies.

Mosquitoes were injected with either dsRNA against AGPAT1 (dsAGPAT1) or dsRNA control (dsControl). Two days post injection, mosquitoes were orally fed with either non-infectious blood or DENV infectious blood. Impact of AGPAT1 depletion on (A) infection rate and (B) DENV gRNA copies at 2 days post oral infection (dpi) with 107 pfu/ml. (C-E) Impact of AGPAT1 depletion on infection rate (C) and DENV gRNA copies at 2 (D) and 7 (E) dpi with 106 pfu/ml. Bars indicate percentage ± s.e. (A, C) or geometric means ± 95% C.I. (B, D, E) with each dot representing one mosquito. *, p-value < 0.05 as indicated by Mann-Whitney test.

https://doi.org/10.1371/journal.ppat.1008199.s011

(TIF)

S1 Table. Identification of mosquito specific metabolites by spectral similarity.

https://doi.org/10.1371/journal.ppat.1008199.s012

(DOCX)

S2 Table. Metabolites detected from cell, midgut and whole mosquito with differential regulation upon DENV infection.

The three tabs contain compound detected on cell, midgut and mosquito with the following information: ionization mode (positive and negative), phase detection (polar and non-polar), the mass average m/z between replicate, the retention time average Rt in minutes between replicates, the MS/MS spectrum fragmentation and intensity, adducts [M+H]+ and [M-H]-, metabolites of importance in Fig 1C, regulated compound with abundance between DENV-infected and uninfected samples (p-value < 0.05 and |log2 fold change| ≥ 1) and annotation classes with the first 3 ranks of identification by MS-Finder.

https://doi.org/10.1371/journal.ppat.1008199.s013

(XLSX)

S3 Table. Primers for dsRNA.

https://doi.org/10.1371/journal.ppat.1008199.s014

(DOCX)

S4 Table. Primers for Real-Time qPCR.

https://doi.org/10.1371/journal.ppat.1008199.s015

(DOCX)

S5 Table. AGPAT FASTA protein sequences.

https://doi.org/10.1371/journal.ppat.1008199.s016

(DOCX)

S6 Table. Metabolites detected from uninfected cell, after AGPAT1 or AGPAT2 depletion.

The table contains compound detected on cell with the following information: ionization mode (positive and negative), the mass average m/z between replicate, the retention time average Rt in minutes between replicates, the MS/MS spectrum fragmentation and intensity, adducts [M+H]+ and [M-H]-, metabolites of importance in Fig 3C, regulated compound with abundance between dsAGPAT1 or 2 and dsControl samples (p-value < 0.05, |log2 fold change| ≥ 1) and annotation classes with the first 3 ranks of identification by MS-Finder.

https://doi.org/10.1371/journal.ppat.1008199.s017

(XLSX)

S7 Table. Metabolites detected from DENV-infected cell at 48 hpi, after AGPAT1 or AGPAT2 depletion.

The table contains compound detected on cell with the following information: ionization mode (positive and negative), the mass average m/z between replicate, the retention time average Rt in minutes between replicates, the MS/MS spectrum fragmentation and intensity, adducts [M+H]+ and [M-H]-, metabolites of importance in Fig 3E, regulated compound with abundance between dsAGPAT1 or 2 and dsControl samples (p-value < 0.05, |log2 fold change| ≥ 1) and annotation classes with the first 3 ranks of identification by MS-Finder.

https://doi.org/10.1371/journal.ppat.1008199.s018

(XLSX)

S8 Table. Metabolites detected from uninfected mosquito after AGPAT1 depletion.

The two tabs contain compound detected on cell after 2 or 7 days post blood feeding with the following information: ionization mode (positive and negative), the mass average m/z between replicate, the retention time average Rt in minutes between replicates, the MS/MS spectrum fragmentation and intensity, adducts [M+H]+ and [M-H]-, metabolites of importance in Fig 4B, regulated compound with abundance between dsAGPAT1 and dsControl samples (p-value < 0.05, |log2 fold change| ≥ 1) and annotation classes with the first 3 ranks of identification by MS-Finder.

https://doi.org/10.1371/journal.ppat.1008199.s019

(XLSX)

S9 Table. Metabolites detected from DENV-infected mosquito after dsAGPAT1 depletion or dsControl.

The table contains compound detected on cell with the following information: ionization mode (positive and negative), the mass average m/z between replicate, the retention time average Rt in minutes between replicates, the MS/MS spectrum fragmentation and intensity, adducts [M+H]+ and [M-H]-, metabolites of importance in Fig 4E, regulated compound with abundance between dsAGPAT1 and dsControl samples (p-value < 0.05, |log2 fold change| ≥ 1) and annotation classes with the first 3 ranks of identification by MS-Finder.

https://doi.org/10.1371/journal.ppat.1008199.s020

(XLSX)

S10 Table. Metabolites detected from DENV-infected or uninfected mosquito after AGPAT1 depletion.

The table contains compound detected on cell with the following information: ionization mode (positive and negative), the mass average m/z between replicate, the retention time average Rt in minutes between replicates, the MS/MS spectrum fragmentation and intensity, adducts [M+H]+ and [M-H]-, metabolites of importance in Fig 4F, regulated compound in dsAGPAT1 condition only, regulated compound with abundance between DENV-infected and Mock samples in dsControl or dsAGPAT1 condition (p-value < 0.05, |log2 fold change| ≥ 1) and annotation classes with the first 3 ranks of identification by MS-Finder.

https://doi.org/10.1371/journal.ppat.1008199.s021

(XLSX)

Acknowledgments

We are grateful to Dr. Mariano Garcia-Blanco for constant support and the Pompon/Garcia-Blanco’s team for comments on a previous version. We thank Dr. Eng Eong Ooi for providing the ST virus. We also thank the Mass spectrometry platform, Université de Toulouse, ICT, UPS, Toulouse, France.

References

  1. 1. Gubler DJ. The Global Emergence/Resurgence of Arboviral Diseases As Public Health Problems. Arch Med Res. 2002;33: 330–342. pmid:12234522
  2. 2. Bhatt S, Gething PW, Brady OJ, Messina JP, Farlow AW, Moyes CL, et al. The global distribution and burden of dengue. Nature. 2013;496: 504–507. pmid:23563266
  3. 3. Capeding MR, Tran NH, Hadinegoro SRS, Ismail HIHM, Chotpitayasunondh T, Chua MN, et al. Clinical efficacy and safety of a novel tetravalent dengue vaccine in healthy children in Asia: a phase 3, randomised, observer-masked, placebo-controlled trial. The Lancet. 2014;384: 1358–1365.
  4. 4. WHO. Dengue vaccine: WHO position paper–September 2018. 2018. Available: http://apps.who.int/iris/bitstream/handle/10665/274315/WER9336.pdf?ua=1
    • 5. WHO SEARO, Regional Office for South-East Asia. Comprehensive guidelines for prevention and control of dengue and dengue haemorrhagic fever. New Delhi, India: World Health Organization Regional Office for South-East Asia; 2011. Available: http://www.searo.who.int/entity/vector_borne_tropical_diseases/documents/SEAROTPS60/en/
      • 6. Liu N. Insecticide Resistance in Mosquitoes: Impact, Mechanisms, and Research Directions. Annu Rev Entomol. 2015;60: 537–559. pmid:25564745
      • 7. Villareal VA, Rodgers MA, Costello DA, Yang PL. Targeting host lipid synthesis and metabolism to inhibit dengue and hepatitis C viruses. Antiviral Res. 2015;124: 110–121. pmid:26526588
      • 8. Pando-Robles V, Oses-Prieto JA, Rodríguez-Gandarilla M, Meneses-Romero E, Burlingame AL, Batista CVF. Quantitative proteomic analysis of Huh-7 cells infected with Dengue virus by label-free LC–MS. J Proteomics. 2014;111: 16–29. pmid:25009145
      • 9. Allonso D, Andrade IS, Conde JN, Coelho DR, Rocha DCP, da Silva ML, et al. Dengue Virus NS1 Protein Modulates Cellular Energy Metabolism by Increasing Glyceraldehyde-3-Phosphate Dehydrogenase Activity. J Virol. 2015;89: 11871–11883. pmid:26378175
      • 10. Fontaine KA, Sanchez EL, Camarda R, Lagunoff M. Dengue Virus Induces and Requires Glycolysis for Optimal Replication. Sandri-Goldin RM, editor. J Virol. 2015;89: 2358–2366. pmid:25505078
      • 11. Birungi G, Chen SM, Loy BP, Ng ML, Li SFY. Metabolomics Approach for Investigation of Effects of Dengue Virus Infection Using the EA.hy926 Cell Line. J Proteome Res. 2010;9: 6523–6534. pmid:20954703
      • 12. El-Bacha T, Midlej V, Pereira da Silva AP, Silva da Costa L, Benchimol M, Galina A, et al. Mitochondrial and bioenergetic dysfunction in human hepatic cells infected with dengue 2 virus. Biochim Biophys Acta BBA—Mol Basis Dis. 2007;1772: 1158–1166. pmid:17964123
      • 13. Heaton NS, Randall G. Dengue Virus-Induced Autophagy Regulates Lipid Metabolism. Cell Host Microbe. 2010;8: 422–432. pmid:21075353
      • 14. Diop F, Vial T, Ferraris P, Wichit S, Bengue M, Hamel R, et al. Zika virus infection modulates the metabolomic profile of microglial cells. PLOS ONE. 2018;13: e0206093. pmid:30359409
      • 15. Cui L, Lee YH, Kumar Y, Xu F, Lu K, Ooi EE, et al. Serum Metabolome and Lipidome Changes in Adult Patients with Primary Dengue Infection. Michael SF, editor. PLoS Negl Trop Dis. 2013;7: e2373. pmid:23967362
      • 16. Martín-Acebes MA, Vázquez-Calvo Á, Saiz J-C. Lipids and flaviviruses, present and future perspectives for the control of dengue, Zika, and West Nile viruses. Prog Lipid Res. 2016;64: 123–137. pmid:27702593
      • 17. Perera R, Riley C, Isaac G, Hopf-Jannasch AS, Moore RJ, Weitz KW, et al. Dengue Virus Infection Perturbs Lipid Homeostasis in Infected Mosquito Cells. PLoS Pathog. 2012;8. pmid:22457619
      • 18. Chotiwan N, Andre BG, Sanchez-Vargas I, Islam MN, Grabowski JM, Hopf-Jannasch A, et al. Dynamic remodeling of lipids coincides with dengue virus replication in the midgut of Aedes aegypti mosquitoes. PLOS Pathog. 2018;14: e1006853. pmid:29447265
      • 19. Xi Z, Ramirez JL, Dimopoulos G. The Aedes aegypti Toll Pathway Controls Dengue Virus Infection. PLoS Pathog. 2008;4. pmid:18604274
      • 20. Martín-Acebes MA, Blázquez A-B, Jiménez de Oya N, Escribano-Romero E, Saiz J-C. West Nile Virus Replication Requires Fatty Acid Synthesis but Is Independent on Phosphatidylinositol-4-Phosphate Lipids. Wang T, editor. PLoS ONE. 2011;6: e24970. pmid:21949814
      • 21. Heaton NS, Perera R, Berger KL, Khadka S, LaCount DJ, Kuhn RJ, et al. Dengue virus nonstructural protein 3 redistributes fatty acid synthase to sites of viral replication and increases cellular fatty acid synthesis. Proc Natl Acad Sci U S A. 2010;107: 17345–17350. pmid:20855599
      • 22. Fagone P, Jackowski S. Membrane phospholipid synthesis and endoplasmic reticulum function. J Lipid Res. 2009;50: S311–S316. pmid:18952570
      • 23. Takeuchi K, Reue K. Biochemistry, physiology, and genetics of GPAT, AGPAT, and lipin enzymes in triglyceride synthesis., Biochemistry, physiology, and genetics of GPAT, AGPAT, and lipin enzymes in triglyceride synthesis. Am J Physiol Endocrinol Metab Am J Physiol—Endocrinol Metab. 2009;296, 296: E1195, E1195–209. pmid:19336658
      • 24. Athenstaedt K, Daum G. Phosphatidic acid, a key intermediate in lipid metabolism. Eur J Biochem. 1999;266: 1–16. pmid:10542045
      • 25. Gillespie LK, Hoenen A, Morgan G, Mackenzie JM. The Endoplasmic Reticulum Provides the Membrane Platform for Biogenesis of the Flavivirus Replication Complex. J Virol. 2010;84: 10438–10447. pmid:20686019
      • 26. Aktepe TE, Mackenzie JM. Shaping the flavivirus replication complex: It is curvaceous! Cell Microbiol. 2018;20: e12884. pmid:29933527
      • 27. Salazar MI, Richardson JH, Sánchez-Vargas I, Olson KE, Beaty BJ. Dengue virus type 2: replication and tropisms in orally infected Aedes aegypti mosquitoes. BMC Microbiol. 2007;7: 9. pmid:17263893
      • 28. Barrows Nicholas J., Mariano A. Garcia-Blanco. Biochemistry and Molecular Biology of Flaviviruses. Chem Rev. 2018;118: 4448–4482. pmid:29652486
      • 29. Tsugawa H. Advances in computational metabolomics and databases deepen the understanding of metabolisms. Curr Opin Biotechnol. 2018;54: 10–17. pmid:29413746
      • 30. van Meer G, Voelker DR, Feigenson GW. Membrane lipids: where they are and how they behave. Nat Rev Mol Cell Biol. 2008;9: 112–124. pmid:18216768
      • 31. Vance JE. Phospholipid Synthesis and Transport in Mammalian Cells. Traffic. 2015;16: 1–18. pmid:25243850
      • 32. Lagace TA, Ridgway ND. The role of phospholipids in the biological activity and structure of the endoplasmic reticulum. Biochim Biophys Acta BBA—Mol Cell Res. 2013;1833: 2499–2510. pmid:23711956
      • 33. Yamashita A, Hayashi Y, Nemoto-Sasaki Y, Ito M, Oka S, Tanikawa T, et al. Acyltransferases and transacylases that determine the fatty acid composition of glycerolipids and the metabolism of bioactive lipid mediators in mammalian cells and model organisms. Prog Lipid Res. 2014;53: 18–81. pmid:24125941
      • 34. Mackenzie JM, Khromykh AA, Parton RG. Cholesterol Manipulation by West Nile Virus Perturbs the Cellular Immune Response. Cell Host Microbe. 2007;2: 229–239. pmid:18005741
      • 35. Yamashita A, Hayashi Y, Matsumoto N, Nemoto-Sasaki Y, Oka S, Tanikawa T, et al. Glycerophosphate/Acylglycerophosphate Acyltransferases. Biology. 2014;3: 801–830. pmid:25415055
      • 36. Dircks LK, Ke J, Sul HS. A Conserved Seven Amino Acid Stretch Important for Murine Mitochondrial Glycerol-3-phosphate Acyltransferase Activity. J Biol Chem. 1999;274: 34728–34734. pmid:10574940
      • 37. Lu B, Jiang YJ, Zhou Y, Xu FY, Hatch GM, Choy PC. Cloning and characterization of murine 1-acyl-sn-glycerol 3-phosphate acyltransferases and their regulation by PPARalpha in murine heart. Biochem J. 2005;385: 469–477. pmid:15367102
      • 38. Yamashita A, Nakanishi H, Suzuki H, Kamata R, Tanaka K, Waku K, et al. Topology of acyltransferase motifs and substrate specificity and accessibility in 1-acyl-sn-glycero-3-phosphate acyltransferase 1. Biochim Biophys Acta BBA—Mol Cell Biol Lipids. 2007;1771: 1202–1215. pmid:17707131
      • 39. Colpitts TM, Cox J, Vanlandingham DL, Feitosa FM, Cheng G, Kurscheid S, et al. Alterations in the Aedes aegypti Transcriptome during Infection with West Nile, Dengue and Yellow Fever Viruses. Rice CM, editor. PLoS Pathog. 2011;7: e1002189. pmid:21909258
      • 40. Gale SE, Frolov A, Han X, Bickel PE, Cao L, Bowcock A, et al. A Regulatory Role for 1-Acylglycerol-3-phosphate-O-acyltransferase 2 in Adipocyte Differentiation. J Biol Chem. 2006;281: 11082–11089. pmid:16495223
      • 41. Merrill AH. Sphingolipid and Glycosphingolipid Metabolic Pathways in the Era of Sphingolipidomics. Chem Rev. 2011;111: 6387–6422. pmid:21942574
      • 42. Lipton BA, Davidson EP, Ginsberg BH, Yorek MA. Ethanolamine metabolism in cultured bovine aortic endothelial cells. J Biol Chem. 1990;265: 7195–7201. pmid:2110161
      • 43. Liebscher S, Ambrose RL, Aktepe TE, Mikulasova A, Prier JE, Gillespie LK, et al. Phospholipase A2 activity during the replication cycle of the flavivirus West Nile virus. PLoS Pathog. 2018;14. pmid:29709018
      • 44. Hollenback D, Bonham L, Law L, Rossnagle E, Romero L, Carew H, et al. Substrate specificity of lysophosphatidic acid acyltransferase β—evidence from membrane and whole cell assays. J Lipid Res. 2006;47: 593–604. pmid:16369050
      • 45. Prasad SS, Garg A, Agarwal AK. Enzymatic activities of the human AGPAT isoform 3 and isoform 5: localization of AGPAT5 to mitochondria. J Lipid Res. 2011;52: 451–462. pmid:21173190
      • 46. Coleman RA, Lee DP. Enzymes of triacylglycerol synthesis and their regulation. Prog Lipid Res. 2004;43: 134–176. pmid:14654091
      • 47. Ashour J, Laurent-Rolle M, Shi P-Y, García-Sastre A. NS5 of Dengue Virus Mediates STAT2 Binding and Degradation. J Virol. 2009;83: 5408–5418. pmid:19279106
      • 48. Merino-Ramos T, Jiménez de Oya N, Saiz J-C, Martín-Acebes MA. Antiviral Activity of Nordihydroguaiaretic Acid and Its Derivative Tetra-O-Methyl Nordihydroguaiaretic Acid against West Nile Virus and Zika Virus. Antimicrob Agents Chemother. 2017;61. pmid:28507114
      • 49. Miller S, Kastner S, Krijnse-Locker J, Bühler S, Bartenschlager R. The Non-structural Protein 4A of Dengue Virus Is an Integral Membrane Protein Inducing Membrane Alterations in a 2K-regulated Manner. J Biol Chem. 2007;282: 8873–8882. pmid:17276984
      • 50. Fajardo-Sánchez E, Galiano V, Villalaín J. Spontaneous membrane insertion of a dengue virus NS2A peptide. Arch Biochem Biophys. 2017;627: 56–66. pmid:28666739
      • 51. Caviglia JM, Dumm INT de G, Coleman RA, Igal RA. Phosphatidylcholine deficiency upregulates enzymes of triacylglycerol metabolism in CHO cells. J Lipid Res. 2004;45: 1500–1509. pmid:15175356
      • 52. Meer G van, Kroon AIPM de. Lipid map of the mammalian cell. J Cell Sci. 2011;124: 5–8. pmid:21172818
      • 53. Richard AS, Zhang A, Park S-J, Farzan M, Zong M, Choe H. Virion-associated phosphatidylethanolamine promotes TIM1-mediated infection by Ebola, dengue, and West Nile viruses. Proc Natl Acad Sci U S A. 2015;112: 14682–14687. pmid:26575624
      • 54. Carnec X, Meertens L, Dejarnac O, Perera-Lecoin M, Hafirassou ML, Kitaura J, et al. The Phosphatidylserine and Phosphatidylethanolamine Receptor CD300a Binds Dengue Virus and Enhances Infection. J Virol. 2016;90: 92–102. pmid:26468529
      • 55. Xu K, Nagy PD. RNA virus replication depends on enrichment of phosphatidylethanolamine at replication sites in subcellular membranes. Proc Natl Acad Sci U S A. 2015;112: E1782–E1791. pmid:25810252
      • 56. Belov G. Less Grease, Please. Phosphatidylethanolamine Is the Only Lipid Required for Replication of a (+)RNA Virus. Viruses. 2015;7: 3500–3505. pmid:26131959
      • 57. Sim S, Jupatanakul N, Dimopoulos G. Mosquito Immunity against Arboviruses. Viruses. 2014;6: 4479–4504. pmid:25415198
      • 58. West J, Tompkins CK, Balantac N, Nudelman E, Meengs B, White T, et al. Cloning and expression of two human lysophosphatidic acid acyltransferase cDNAs that enhance cytokine-induced signaling responses in cells. DNA Cell Biol. 1997;16: 691–701. pmid:9212163
      • 59. Kanthong N, Laosutthipong C, Flegel TW. Response to Dengue virus infections altered by cytokine-like substances from mosquito cell cultures. BMC Microbiol. 2010;10: 290. pmid:21078201
      • 60. Barletta ABF, Silva MCLN, Sorgine MHF. Validation of Aedes aegypti Aag-2 cells as a model for insect immune studies. Parasit Vectors. 2012;5: 148. pmid:22827926
      • 61. Schreiber MJ, Holmes EC, Ong SH, Soh HSH, Liu W, Tanner L, et al. Genomic Epidemiology of a Dengue Virus Epidemic in Urban Singapore. J Virol. 2009;83: 4163–4173. pmid:19211734
      • 62. Tsugawa H, Cajka T, Kind T, Ma Y, Higgins B, Ikeda K, et al. MS-DIAL: data-independent MS/MS deconvolution for comprehensive metabolome analysis. Nat Methods. 2015;12: 523–526. pmid:25938372
      • 63. Tsugawa H, Kind T, Nakabayashi R, Yukihira D, Tanaka W, Cajka T, et al. Hydrogen Rearrangement Rules: Computational MS/MS Fragmentation and Structure Elucidation Using MS-FINDER Software. Anal Chem. 2016;88: 7946–7958. pmid:27419259
      • 64. Sumner LW, Amberg A, Barrett D, Beale MH, Beger R, Daykin CA, et al. Proposed minimum reporting standards for chemical analysis Chemical Analysis Working Group (CAWG) Metabolomics Standards Initiative (MSI). Metabolomics Off J Metabolomic Soc. 2007;3: 211–221. pmid:24039616
      • 65. Tsugawa H, Nakabayashi R, Mori T, Yamada Y, Takahashi M, Rai A, et al. A cheminformatics approach to characterize metabolomes in stable-isotope-labeled organisms. Nat Methods. 2019;16: 295. pmid:30923379
      • 66. Chong J, Soufan O, Li C, Caraus I, Li S, Bourque G, et al. MetaboAnalyst 4.0: towards more transparent and integrative metabolomics analysis. Nucleic Acids Res. 2018;46: W486–W494. pmid:29762782
      • 67. Fraiture M, Baxter RHG, Steinert S, Chelliah Y, Frolet C, Quispe-Tintaya W, et al. Two Mosquito LRR Proteins Function as Complement Control Factors in the TEP1-Mediated Killing of Plasmodium. Cell Host Microbe. 2009;5: 273–284. pmid:19286136
      • 68. Johnson BW, Russell BJ, Lanciotti RS. Serotype-Specific Detection of Dengue Viruses in a Fourplex Real-Time Reverse Transcriptase PCR Assay. J Clin Microbiol. 2005;43: 4977–4983. pmid:16207951
      • 69. Manokaran G, Finol E, Wang C, Gunaratne J, Bahl J, Ong EZ, et al. Dengue subgenomic RNA binds TRIM25 to inhibit interferon expression for epidemiological fitness. Science. 2015;350: 217. pmid:26138103