Anopheles coluzzii stearoyl -CoA desaturase is essential for adult female survival and reproduction upon blood feeding

Abstract

Inside the mosquito midgut lumen, red blood cell membranes are broken down, and hemoglobin and other proteins of the blood and serum are digested into amino acids, chiefly through the action of bloodmeal-induced trypsin-like and other proteases. Human red blood cells provided by the National Blood Service of the United Kingdom National Health Service were obtained from healthy donors upon taking a written informed consent. Briefly, SCD1-specific oligonucleotide primers flanked at their 5′ end by T7 RNA polymerase promoter sequence (SCD1-RNAi-F: taatacgactcactataggg-AAATGTGATTGCCTTCGGT; SCD1-RNAi-R: taatacgactcactataggg-GCGAGAAGAAGAAGCCAC) were used in PCR reactions on total RNA extracted from 10 sugar-fed adult female mosquitoes using the TRIzol reagent (Invitrogen). Gene expression profiling Control (dsLacZ-injected) and SCD1 KD female mosquitoes were fed on human blood though a standard membrane feeding assay 3 days after recovery from dsRNA injection, and 10 fully-engorged female mosquitoes from both groups were sampled at each of the following time points: 0h (just prior to BF), 6h, 12h, 18h and 24h PBM.

Introduction

In preparation for reproduction, anautogenous female mosquitoes must ingest vertebrate blood to improve their energy status and to stimulate vitellogenesis. This induces massive mobilization of lipids from the fat body to the ovaries where they serve as the principal energy source for the maturing oocytes and the developing embryos [1–3]. As bloodmeals are required in every gonotrophic cycle, this has made mosquitoes ideal vectors of blood-borne pathogens such as the malaria parasites.

Inside the mosquito midgut lumen, red blood cell membranes are broken down, and hemoglobin and other proteins of the blood and serum are digested into amino acids, chiefly through the action of bloodmeal-induced trypsin-like and other proteases. Amino acids are then used for protein and fatty acid biosynthesis from reduced carbon atoms [3,4]. Specifically, carbon skeletons of amino acids are converted to pyruvate, acetyl CoA, or one of TCA cycle intermediates, which in turn are used for glucose synthesis through gluconeogenesis or used as intermediates for the synthesis of dietary triacylglycerols (TAGs). At the same time, blood-derived TAGs are converted into monoacylglycerols (MAGs) and free fatty acids (FFAs) in the midgut lumen through the action of lipases before they are transported into the midgut cells and converted into diacylglycerides (DAGs), TAGs and phospholipids [1]. TAGs are stored intracellularly within lipid droplets and are transported principally to the fat body from where they provide energy upon demand. TAGs are also synthesized from dietary carbohydrates through glycolysis and gluconeogenesis. The loading and unloading of lipids into and from the fat body is performed by a shuttling mechanism that involves lipophorin (Lp) and other lipid transfer molecules. Lp also caries DAGs from the midgut directly to the ovarian tissues where they convert into TAGs [5]. In addition to their importance as energy resource, lipids are crucial for signal transduction and maintenance of membrane fluidity of the midgut cells [6,7].

Saturated (SFAs) and unsaturated (UFAs) fatty acids, the building blocks of TAGs and other fats, are carboxylic acids consisting of an aliphatic tail and a terminal carboxyl group. The introduction of the first double bond at the Δ9 position of the aliphatic chain (between carbons 9 and 10) of SFAs, the most critical commitment step in the biosynthesis of mono UFAs (MUFAs), is catalyzed by SCD, an iron-containing, microsomal enzyme conserved across eukaryotes. The reaction requires NADH, cytochrome b5 and molecular oxygen [8]. Major MUFAs such as oleic acid and palmitoleic acid are key precursors of membrane phospholipids, cholesterol esters and triglycerides, and are synthesized from SFAs such as stearic acid and palmitic acid, respectively [8]. The ratio of saturated to unsaturated phospholipids is critical for maintenance of cell membrane fluidity and homeostasis.

SCDs have been investigated for their potential use as chemotherapeutic targets in human metabolic disorders, cancer and infections. For example, the Mycobacterium tuberculosis SCD is a known target of the thiourea drug isoxyl [9], and Sterculic oil (SO) is the best recognized natural inhibitor of SCD [10]. The importance of lipid biosynthesis during mosquito blood feeding offers an opportunity for exploitation toward development of novel vector control interventions. Here, we investigate the function of SCD1 in the malaria vector mosquito Anopheles coluzzii upon blood feeding. A. coluzzii was until recently known as A. gambiae M form [11] and is the species used in the majority of experimental studies, most of which still referred to as A. gambiae. Silencing SCD1 through RNA interference and inhibition of SCD1 enzymatic activity by supplementing bloodmeals with SA allowed us to investigate the transcriptional, metabolic and phenotypic impacts of the enzyme on mosquito physiology. The data has led us to conclude that SCD1 is essential for mosquito lipid metabolism, survival and reproduction following a bloodmeal and identify SCD1 as a potential target of novel vector control interventions.

Methods

Ethics statement

All animal procedures were reviewed and approved by the Imperial College Animal Welfare and Ethical Review Body and carried out in accordance with the Animal Scientifics Procedures Act 1986 under the United Kingdom Home Office licenses PLL70/7185 and PPL70/8788. Human red blood cells provided by the National Blood Service of the United Kingdom National Health Service were obtained from healthy donors upon taking a written informed consent.

Mosquito colony and maintenance

A. coluzzii mosquitoes used in all experiments were of the N’gousso colony established from A. gambiae M molecular form and the Forest chromosomal form mosquitoes collected in Cameroon in 2006. Adult mosquitoes were held in netted cages and maintained on 5% fructose solution at 28°C (±1°C), 70–80% relative humidity and a 12-hour dark/ light cycle. Larvae were fed on tetraamine fish food.

Bioinformatics analysis

BLAST searches were conducted using the National Centre for Biotechnology Information (NCBI) and Vector Base sequences and tools. Sequence alignments were performed using locally installed ClustalW and protein prediction was done using the SMART Package (http://smart.embl-heidelberg.de/). Pathway analysis of transcriptomic and metabolomic data was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG).

DsRNA preparation and gene silencing

Gene-specific double-stranded RNA (dsRNA) was synthesized and injected as described [12]. Briefly, SCD1-specific oligonucleotide primers flanked at their 5′ end by T7 RNA polymerase promoter sequence (SCD1-RNAi-F: taatacgactcactataggg-AAATGTGATTGCCTTCGGT; SCD1-RNAi-R: taatacgactcactataggg-GCGAGAAGAAGAAGCCAC) were used in PCR reactions on total RNA extracted from 10 sugar-fed adult female mosquitoes using the TRIzol reagent (Invitrogen). DsRNA was synthesized using the T7 MEGAscript kit (Invitrogen) and the PCR product was purified using the RNeasy kit (QIAGEN). DsRNA for the control LacZ gene was synthesized as reported [13]. Gene silencing was accomplished by injection of 69 nL (3 ug/uL) of dsRNA into the mesothoracic spiracle of each carbon dioxide-anesthetized 0-2-day old female mosquito. Mosquitoes were allowed to recover for 3 days before any further treatment.

RNA extraction and qRT PCR

Total RNA was extracted from 10 dsSCD1 and dsLacZ injected mosquitoes, respectively, using the using the TRIzol reagent (Invitrogen) and the QIAGEN RNeasy kit. cDNA synthesis was performed from total 1 μg RNA following the manufacturer’s instructions for the QuantiTect Reverse Transcription kit (QIAGEN). Reactions of 20ul consisted of 1ul of Quantiscript Reverse Transcriptase, 4 ul of Quantiscript RT Buffer and 1ul of RT primer mix. Transcript abundance / gene silencing efficiency were determined using the SYBR Premix Ex Taq kit (Takara) according to the manufacturer’s instructions and the ABI Prism 7500 Fast Real-Time thermocycler (Applied Biosystems). The A. gambiae ribosomal gene S7 was used as a reference for constitutive expression. The primers used were SCD1-qRT-F (GGTGTCGAAGGAGATCGTGG), SCD1-qRT-R (CTGGTTGAGAATGGTGGCG), S7A (GGCGATCATCATCTACGT), S7B (GTAGCTGCTGCAAACTTCGG), CEC1F (GGAAGCGGGACGCCTGAA) and CEC1R (CCTTGACACCTGCCACCACC), as reported previously [13–15]. In brief, each 25 μl reaction consisted of 2.5 ng template RNA, 0.2 μM forward and reverse primers and 12.5 μl of the SYBR premix. Amplification protocol is as follows: 95°C for 5–10 seconds, at (94°C for 3 seconds, 60°C for 20 seconds annealing and 72°C for 30 seconds extension for 30 cycles).

Transmission electron microscopy

For Transmission Electron Microscopy (TEM), midguts from dsLacZ or dsSCD1 injected A. gambiae females were dissected just before blood feeding (0h) and 24h post bloodmeal (PBM) through a standard membrane feeding apparatus (Hemotech). Dissected midguts were fixed for 2h at room temperature (2.5% E.M. Grade Glutaraldehyde, 0.1 M sodium cacodylate buffer containing 2 mM calcium chloride), transferred to 1% osmium tetroxide for 1 h, washed in deionized water, and stained with 2% aqueous uranyl acetate for 30 min. Specimens were dehydrated with ethanol dilutions, infiltrated with Spurr’s resin, and flat-embedded at 60°C. Transverse sections (70 nm) were prepared using a Leica ultracut UCT ultramicrotome and placed onto 200 thin mesh bar copper grids followed by counterstaining with lead citrate. Midgut sections were viewed with a JEOL-JEM-1230 electron microscope operated at 80 kV and TIFF images (8 bit) were captured using an AMT 4M pixel camera.

Gene expression profiling

Control (dsLacZ-injected) and SCD1 KD female mosquitoes were fed on human blood though a standard membrane feeding assay 3 days after recovery from dsRNA injection, and 10 fully-engorged female mosquitoes from both groups were sampled at each of the following time points: 0h (just prior to BF), 6h, 12h, 18h and 24h PBM. Total RNA was extracted using the TRIzol reagent (Invitrogen) and purified using the RNeasy kit (Qiagen). RNA quantification was performed using the Nanodrop 1000 spectrophotometer (Thermo Scientific) and RNA quality and integrity was assessed though gel electrophoresis (1% agarose gel) and the Agilent Bioanalyzer 2100 (Agilent Technologies). Labelling and hybridization were performed, following the instructions in the Low Input Quick Amp Labeling kit, for two-color microarray-based expression analysis (Agilent). The Pfalcip_Agamb2009 microarray platform (designed by Dr. Vlachou, personal communication) was used for all of the hybridizations presented. This platform encompasses a total of 15,424 A. gambiae annotated transcripts of the AgamP3.8 according to the 2014 VectorBase release along with P. falciparum probes. A. coluzzii genome assembly suggests 14,712 gene transcripts and shows minor differences from A. gambiae Slides were scanned using the Genepix 4000B scanner equipped with Genepix Pro 6.1 software (Axon instruments). Data were normalized using the Genespring 11.0 GX software (Agilent). The Lowess Normalization Method with the threshold of raw signals set to 5 which was enough to eliminate background regulation of P. falciparum probes. Significant differences in gene expression between experiment and control groups were evaluated by combining an expression ratio cut-off of 1.0 on a log2 scale and one-way ANOVA statistical tests (with P values of ≤0.025) across different time points. Expression profiles were clustered using the Cluster software, version 3.0, according to the Pearson correlation score. The data were visualized with Java TreeView version 1.1.6r4. All subsequent data analysis was done in Microsoft Excel (Microsoft Corp., Redmond, WA).

Preparation of SA solution and mosquito treatment

SA was freshly synthesized and kindly provided by Professor Mark Baird in (Bangor University, UK), and 1M stock solution was prepared by dissolving SA in 1M Tris buffer and kept at -200 C. For experiments, SA stock solution was diluted further into Tris buffer and human blood to a final concentration of 1mM. Four-day old female mosquitoes were fed on human blood supplemented with a 1mM SA (final concentration in blood) in 1mM Tris buffer or with 1mM Tris buffer alone. Mosquitoes were maintained on 5% sucrose solution throughout the experiments. For survival analysis, dead mosquitoes from each group were counted every 6h post bloodmeal (PMB). Before performing experiments with 1mM SA concentration, different concentration of SA including 0.4, 0.6, 0.8, 1, 1.4, 1.6 and 1.8 mM were tested to determine the IC50.

Metabolomic profiling

SCD1 KD and dsLacZ-injected (control), and SA-treated and 1mM Tris-treated (control) mosquitoes were sampled at 0h (just prior to blood feeding), 18h and 36h PBM. Sample preparation for metabolomic profiling and data analysis were performed as previously described [16]. In brief, 3 mosquitoes, from each time point, were extracted in a 1 ml ice-cold solution of methanol: water (8:2 v/v) and centrifuged (14000 rpm, 40 C, 15 min) to separate the cell debris from the supernatant. The supernatant was then transferred to a silanized 1.5 ml glass vial (Agilent Technologies UK Ltd) and dried in a SpeedVac concentrator (Eppendorf). Methoxymation, in combination with trimethylsilylation (methoxy-TMS), was accomplished for derivatization. Sample analysis was performed on an Agilent 7890 GC coupled to a 5975c mass spectrometer, using the Fiehnlib settings, running a pooled biological quality control (pbQC) sample every six samples [17]. Myristic-d27 acid was used as the endogenous reference. Deconvolution and integration of the extracted metabolites was performed using an in-house work-flow as described [18]. All experiments were carried out in 5 independent biological replicates with 3 technical replicates per each biological replicate. Statistical analysis was performed using a student’s t-test. To produce heat maps, metabolite data was normalized to the internal standard (Myristic-d27 acid) and pbQC levels, then the averages of each metabolite from five replicates for control or treated group at different time points were calculated. Finally, each metabolite was calibrated to the highest value for that metabolite across all time points and groups. Significant differences in metabolite signals between experimental and control groups were evaluated by combining metabolite ratio cut-off 0.1 on a log2 scale. MATLAB was used to plot the data.

Antibiotic treatment

To reduce the natural microbiota load in the midgut, mosquitoes were treated with antibiotic solutions as described previously [13]. In brief, mosquitoes were collected just after emergence and kept on a cocktail of 25 μg/ml gentamycin, 10 μg/ml penicillin and 10 units/ml streptomycin, diluted in 5% sucrose. The treatment continued for 5 days, with the antibiotic solution refreshed every 24 hours. Total RNA was extracted from 10 dsSCD1 and dsLacZ-treated mosquitoes per replicate, followed by cDNA synthesis from total RNA was performed using the QuantiTect Reverse Transcription kit (Qiagen). 16S rDNA primers used were as reported previously [13].

Results and discussion

Putative stearoyl-CoA desaturases

Bioinformatics searches of the A. gambiae and A. coluzzii genomes identified AGAP001713 as a putative stearoyl-CoA Δ9-desaturase, designated as SCD1. The gene encodes a 355 amino acid-long protein with a single transmembrane domain and a putative acyl-CoA desaturase domain (IPR015876; 68–309 aa) showing 45.1%, 92.2% and 71.9% identity to equivalent domains of SCD homologues in human, Drosophila melanogaster and Plasmodium falciparum, respectively (S1 Fig). The desaturase domain contains three Histidine (His) boxes designated as region Ia, Ib and II, which together contain 8 conserved His residues that are essential for SCD1 activity.

To determine the relationships of SCD1 with other putative stearoyl-CoA desaturases encoded in the A. gambiae/A. coluzzii genomes we first searched VectorBase for proteins matching the acyl-CoA desaturase domain. This identified six additional genes on chromosome 2R and one on the UNKN chromosome (S1 Table). The maximum percent identity of amino acid matches between SCD1 and each of these proteins was 55% (AGAP003418), indicating that none of these are likely to be recently duplicated paralogues of SCD1. Secondly, to check for any potential homologues not currently annotated in the genome we used the SCD1 protein sequence to perform a tBLASTn search (default parameters) of the A. gambiae/A. coluzzii genomes. No additional loci were identified, the hit with the highest identity (60%) encoded the AGAP003418 homologue of SCD1 (S1 Table) We then queried VectorBase for orthologues of SCD1, and identified 1-to-1 orthologues in 19 other mosquitoes that were orthologous to Desat1 and Desat2 from Drosophila melanogaster (S2 Table). Target and query percent identities with mosquito 1-to-1 orthologues ranged from 100% to 70% (mean = 96%) and with the fly orthologues from 75%-63% (mean = 69%). Therefore, the orthologues of SCD1 across mosquitoes, and furthermore in flies, all showed higher sequence identities than SCD1 has with any of the putative other putative stearoyl-CoA desaturases encoded in the A. gambiae/A. coluzzii genomes.

SCD1 is essential for adult female survival upon blood feeding

Newly emerged, age-matched (0-2-days old) female mosquitoes were injected with control LacZ and SCD1 dsRNAs, respectively, kept separately in netted cups and maintained on 5% w/v glucose solution. SCD1 silencing efficiency 3 days post recovery was 78% as assessed by qRT-PCR (Fig 1A). At day 3 post recovery, the mosquitoes were fed on human blood through a membrane feeder. Fully engorged mosquitoes were separated, counted and maintained on 5% w/v glucose solution throughout the rest of the experiment. Dead mosquitoes from each group were counted every 6 hours post bloodmeal (PBM). The results revealed a dramatic increase of mortality in SCD1 KD blood-fed mosquitoes starting at about 24h PBM and reaching 100% by 48h PBM (Fig 1B). The median survival time was 42h. Only 17% mortality was observed in control mosquitoes at 48h PBM, which was significantly higher than that of SCD1 KD mosquitoes (p<0.0001).

[Figure omitted. See PDF.]

Fig 1. Blood and sugar meal induced mortality in SCD1 KD and SA treated female A. coluzzii.

(A) Efficiency of RNAi-mediated silencing of SCD1. Silencing was measured with qRT-PCR 3 days post recovery following dsRNA injection of freshly hatched adult females. Control (CTRL) mosquitoes were injected with dsRNA against the LacZ gene. (B-D) Survival analysis of adult female mosquitoes that were: (B) injected with SCD1 or control (CTRL) LacZ dsRNAs following feeding on human blood (post blood feeding; PBF) at day 3 post injection, (C) treated with 1mM SA or control Tris buffer, added to human blood provided through a membrane feeder, or (D) injected with SCD1 KD or CTRL LacZ dsRNA and maintained on 5% sucrose for 25 days post injection recovery that started at day 3 post injection (PIR). Survival was recorded every 6 h. Data were analyzed using Kaplan-Meier survival analysis and log-rank statistics (Mantel-Cox test), showing highly statistical differences between treatments and controls in all three sets (P<0.0001). The mean ± SEM of data obtained from 3 independent biological experiments and the median survival time (dashed lines) are shown.

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

Survival assays were repeated, this time following supplementation of the human blood with 1 mM SA in 1 mM Tris buffer instead of SCD1 silencing. Human blood supplemented with 1 mM Tris buffer alone was used as a control. Similar to the SCD1 KD assays, mosquitoes fed on blood supplemented with SA also exhibited increased mortality rates compared to controls, albeit not as rapid, which reached 52% at 42h PBM (p<0.0001; Fig 1C). Interestingly, the mortality rate was reduced after 42 h, presumably due to quick SA break down.

We examined whether the observed survival phenotype was specific to blood-fed mosquitoes by monitoring the survival of SCD1 KD and dsLacZ control mosquitoes that were not blood-fed but maintained on 5% sucrose solution throughout the experiment. Whilst mortality rates of SCD1 KD sugar-fed mosquitoes were evidently slower than the SCD1 KD blood-fed mosquitoes, these rates were significantly higher than their sugar-fed controls (p<0.0001; Fig 1D). The median survival time for the SCD1 KD group was 6 days and for the control group was 11.5 days.

SCD1 silencing affects midgut epithelial integrity and egg development

Light microscopy 24h PBM revealed a blood-filled thoracic cavity in SCD1 KD mosquitoes compared to controls, indicating blood perfusion of the hemocoel compartment (Fig 2A, top panel). In addition, SCD1 KD midguts were extremely fragile and appeared swollen, with bloodmeal not contained within a well-defined blood bolus (Fig 2A, middle panel). These phenotypes suggested that SCD1 KD may affect the midgut physiology and cell membrane fluidity, directly or indirectly compromising the ability of mosquitoes to cope with the mechanical and metabolic stress of blood feeding.

[Figure omitted. See PDF.]

Fig 2. Phenotypic characterization of SCD1 KD mosquitoes.

(A) Representative pictures of thoraces, midguts and ovaries of blood-fed control (left) and SCD1 KD (right) female mosquitoes at 24h PBM. (B) Representative TEM pictures of midgut epithelial cells dissected from control (left) and SCD1 KD (right) female mosquitoes at 0h before blood meal and 24h PBM, respectively. LCM denotes lateral borders between adjacent midgut cells. Thickness of LCM is the mean value of 10 different points on the same LCM. LD denotes lipid droplets and MV denotes microvilli. Scale bar is set at 1 μm.

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

We investigated the subcellular basis of the SCD1 KD phenotype using transmission electron microscopy (TEM) on midgut tissues at 0h (just before bloodmeal) and 24h PBM. Imaging revealed that midgut cell membranes were thicker in SCD1 KD (mean width 60 nm) compared to control (31nm) mosquitoes prior to bloodmeal (Fig 2B, top panel), suggesting loss of cell membrane rigidity in agreement with previously reported induction of negative cell membrane curvature [19,20]. At 24h PBM, cell membranes were substantially thicker in SCD1 KD (55 nm) than control (14 nm) midguts. Noteworthily, cell membrane thickness in SCD1 KD midguts was only marginally reduced after bloodmeal (60 nm compared to 55 nm respectively). Reduction of cell membrane thickness is a known adjustment to the expansion of the epithelial cell surface in order to provide structural support to the engorged midgut epithelium [21]. We conclude that SCD1 loss-of-function alters the ability of midgut epithelial cells to undergo adjustments needed to offset the increased pressure within the lumen, causing rupturing of the epithelial cell layer upon blood feeding and blood perfusion of the thoracic cavity.

TEM analysis also revealed a great deficiency in lipid droplets in SCD1 KD compared to control midgut cells at 24h PBM (Fig 2B, bottom panel). Therefore, blockage of SCD1-catalyzed MUFA synthesis appears to inhibit TAG and phospholipid biosynthesis from bloodmeal-derived carbon in the midgut cells as reported previously [4]. These data prompted us to examine the effect of SCD1 KD on egg development. Ovaries are a major target tissue of Lp-carried lipids synthesized in the midgut from bloodmeal-derived carbon. Indeed, SCD1 KD mosquitoes showed no signs of egg maturation following a bloodmeal (Fig 2A, bottom panel), similar to the phenotype previously observed upon silencing Apolipophorin II/I (ApoII/I), the main component of the insect Lp [22]. We conclude that compromised lipid synthesis in the midgut cells due to inhibition of SCD1 is responsible for reduced lipid shuttling in the hemolymph leading to underdeveloped ovaries.

SCD1 silencing greatly impacts on metabolic homeostasis

We investigated the molecular basis of the SCD1 KD phenotype in female mosquitoes in response to blood feeding using DNA microarray expression profiling of SCD1 KD versus dsLacZ-injected control mosquitoes. Five time points from three independent biological replicates were examined: 0h (just prior to blood feeding), 6h, 12h, 18h and 24h PBM. In total, 1432 genes registered expression values in at least one of the time points examined. Genes with at least 2-fold expression changes (p<0.05 in ANOVA t-test) in SCD1 KD versus control mosquitoes, in any of the time points, were analyzed further. To derive a general picture of expression changes, genes were functionally classified into 9 major functional categories as done previously [22]. Overall, upregulation dominated downregulation in most of the time points; however, there was a marked difference between the different categories (Fig 3 and S3 Table). SCD1 KD considerably disturbed metabolic homeostasis that is essential for mosquito survival and reproduction, as well as other gene functional categories, all of which are presented in this section as they provide an integrated picture of the disruption of mosquito physiology. The impact of SCD1 KD on the immune response is analyzed later.

[Figure omitted. See PDF.]

Fig 3. Gene transcription profiles upon blood feeding of SCD1 KD mosquitoes.

(A) Heat map showing the transcription profiles of a total of 137 genes exhibiting at least 2-fold differential expression in SCD1 KD vs. dsLacZ-injected control across all time points. The time points examined are 0h (just prior to blood meal), 6h, 12h, 18h and 24h post bloodmeal (PBM). Blue color indicates downregulation and yellow color indicates upregulation, as shown in the color bar below the heat map. The data shown is the mean from 3 independent experiments. The number of genes grouped into each of 8 functional categories or are of unknown function is shown. (B) The graphs show the percentage of genes that are upregulated or downregulated in each time point and belonging to each of the 8 functional categories or having unknown function (grey colored bars). The color code of bars representing each of the functional categories is the same as the color code shown in panel A.

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

Carbohydrate and lipid metabolism.

A total of 79 genes involved in “carbohydrate and lipid metabolism” were differentially regulated in at least one time point; of these, 73% were downregulated and 19 were differentially regulated across the time course. Approximately 39% and 29% of all the genes differentially regulated at 6h and 24h PBM belonged to this category.

Differential expression of lipid metabolism genes was noticed (21/28), with several genes involved in lipid biosynthesis (6/10), lipolysis (5/6) and lipid transportation (4/9) being downregulated. Genes encoding acyl-CoA carboxylases (ACC) that catalyze the irreversible conversion of acetyl-CoA to malonyl-CoA were also downregulated. As silencing SCD1 increases SFA accumulation by inhibiting SFA to MUFA conversion, it is hypothesized that SFA buildup over time hinders ACC activity through a well-known feedback mechanism resulting in a drop of malonyl-CoA and increased transport of acetyl-CoA into the mitochondria for β-oxidation [23]. 3-hydroxy-3-methylglutaryl-CoA lyase (HMG-CoA lyase) that releases acetyl-CoA by breaking down leucine and fats [24] was downregulated across all time points. We speculate that HMG-CoA lyase downregulation is due to the increased flux of acetyl-CoA caused by over-activation of β-oxidation. UDP glucuronosyltransferase (UGT) that catalyzes glucuronidation removing toxic compounds such as fatty acid derivatives from body was also downregulated. Insect UGTs play important roles in cuticle formation, pigmentation and olfaction, besides their major role in inactivation and excretion of endogenous and exogenous compounds [25,26]. This may reflect a feedback mechanism to prevent the removal of fatty acids.

Importantly, Lp was markedly downregulated (>3 fold, at 6h and 12h PBM) in SCD1 KD mosquitoes, which may highlight a feedback regulatory mechanism caused by the depletion of lipid droplets in the midgut cells due to the blockade of SFA to MUFA conversion. Silencing ApoII/I is previously shown to inhibit egg development, which is consistent with the SCD1 KD phenotype reported here [27,28]. Surprisingly, phospholipid scramblase (PLS1), a plasma membrane-anchored enzyme that shuttles phospholipids between the inner and outer leaflets of the cell membrane to maintain membrane integrity, was upregulated (>2 folds, at 6h PBM). This response may represent an effort to sustain the SFA/ MUFA ratio in SCD1 KD cells.

Genes encoding TCA cycle enzymes including isocitrate dehydrogenase (IDH) that catalyzes the reversible conversion of isocitrate to alpha-ketoglutarate, and succinate-CoA synthetase (SUCS-GDP forming) and succinate-CoA synthetase (SUCS-ADP forming), which catalyze the conversion of succinate and GTP or ATP to succinyl-CoA and ADP or GDP, respectively, were also downregulated at 6h PBM. This suggests higher succinyl-CoA involvement in the TCA cycle [29] in SCD1 KD mosquitoes. Limitation of free CoA in SCD1 KD mosquitoes may also be an important factor for the slowdown of the forward reaction of TCA cycle and accumulation of TCA cycle intermediates.

A large number of carbohydrate metabolism genes (41) including many involved in glycolysis and gluconeogenesis (11) as well as pyruvate/acetyl-CoA metabolism (8) were differentially regulated in SCD1 KD mosquitoes. Interestingly, α-glucosidase (AGL) and α-amylase (AAM), which convert glycogen/starch and their intermediates to glucose were downregulated (>2 fold at 0h, 6h and 18h PBM), which could imply lower levels of glucose and other downstream metabolites of glucose. However, downregulation of two key gluconeogenesis enzymes, namely pyruvate carboxylase (PC) and phosphoenol pyruvate carboxykinase (PEPCK), was also observed. At the same time, galactokinase (GALK), a key enzyme of the Leloir pathway, was upregulated >2 fold at 0h PBM). Since the Leloir pathway converts galactose to UDP-glucose in order to maintain homeostasis during glucose scarcity, it is possible that activation of this pathway together with PC and PEPCK upregulation is part of an effort to compensate for glucose deficiency.

Protein metabolism.

As many as 46 genes involved in the “protein degradation” category was differentially regulated in SCD1 KD compared to control mosquitoes, of which 67% were downregulated, mostly at 24h PBM. Seven (7) of these genes were differentially regulated across the experimental time course. Mosquitoes require a variety of proteolytic enzymes to digest the bloodmeal. Among downregulated genes were several exopeptidases, including alanine aminopeptidase, serine protease S28, astacin, carboxypeptidase A1 and various trypsins. Nevertheless, some proteolytic enzymes were upregulated including peptidase M19 (6h PBM), peptidase S33 (12h PBM), trypsin 4 (18h PBM) and peptidase S60 (24h PBM). It is reported that enzymes expressed early after a bloodmeal such as trypsin 4 may indirectly activate the transcription of the major endoproteolytic enzyme trypsin 1 [30]. Taken together, these data suggest that a complex network of proteolytic processes involved in bloodmeal metabolism are differentially affected in SCD1 KD mosquitoes contributing to the dramatic phenotype observed.

Arginosuccinate synthase (ASS) and arginosuccinate lyase (ASL), two key enzymatic components of the arginine biosynthetic pathway, whereby citrulline is converted to arginine releasing fumarate, were downregulated in SCD1 KD mosquitoes (18h and 24h PBM). It is plausible that accumulation of fumarate leads to inhibition of ASS and ASL expression [31]. The upregulation of ornithine decarboxylase (ODC), which converts ornithine to putrescine may suggest further metabolism of putrescine initially to spermidine and then through the glutathione detoxification pathways.

Cytoskeleton, cell adhesion and structural components.

Several genes (27) belonging to the “cytoskeleton, cell adhesion and structural components” category showed differential expression throughout the time course; the majority of them were upregulated. Importantly, upregulation of an actin-like gene was detected at 6h PBM. It was previously reported that increased levels of actin mRNA at 3-4h PBM correlates with epithelial remodeling to accommodate distention of the midgut upon blood ingestion [21]. As SCD1 KD causes rigidity of the epithelial cell membranes, this is suggestive of a stress response to overcome the mechanical pressure caused by the rigid cell membranes. E-Cadherin, a Calcium-dependent cell adhesion protein involved in cell-cell interactions was previously shown to be overexpressed in SCD1 silenced breast cancer cells [32]. Consistent with this finding, an upregulation of a gene encoding a cadherin-like protein was observed 24h PBM. This protein is expressed in the midgut of A. gambiae larvae, where it plays a key role in detoxification of the nonchemical larvicide-Bti [33,34]. Thus, its overexpression could serve in detoxifying excess SFAs.

A gene encoding a chitin-binding domain protein of the peritrophic matrix (PM) was consistently downregulated throughout the time course but 18h PBM, with a peak >5-fold decrease at 24h PBM. A. gambiae PM is detectable as early as 12h PBM and fully formed by 48h PBM [35]. We hypothesize that PM synthesis is compromised in SCD1 KD mosquitoes owing to the downregulation of this gene. Indeed, the juvenile hormone (JH) binding protein (JHBP) was also upregulated at 24h PBM. JH plays important roles in larval development, reproduction and stress management and is known to suppress PM synthesis and function in adult Calliphora erythocephala [36]. As JH levels drop with increasing JHBP levels [37,38], JHBP upregulation in SCD1 KD mosquitoes may play part in suppression of ovary development and PM synthesis.

Signaling/ATPases/GTPases.

JHBP belongs to the “signaling/ATPase/GTPase” gene category. Of this category, 52 genes in total were differentially regulated throughout the time course, mostly at 6h and 24h PBM representing a quarter and a third of all upregulated genes at these time points, respectively; the majority of them upregulated and 15 of them differentially regulated across the time course.

Insulin-like peptides (ILPs) that can have a multitude of functions including lipid and glycogen processing, immunity, reproduction and longevity [39,40] were significantly downregulated across 5 time points. In Ae. aegypti, ILP3 is known to stimulate ovaries to synthesize ecdysteroid hormones including 2-hydroxyecdysone that induces production of vitellogenin in the fat body [41,42]. A possible reduction of ecdysteroids could also lead to lower levels of JH [37,43]. These data could suggest interruption of insulin mediated signaling in SCD1 KD mosquitoes and perhaps defective ecdysteroid hormone synthesis.

The insulin/mTOR (mammalian target of rapamycin) signaling plays a central role in regulating cell growth and proliferation by sensing and integrating a variety of inputs arising from amino acids, cellular stresses, energy status and growth factors [41]. The phosphatidylinositol 5-phosphate 4-kinase (PIP4K), one of the essential mediators of many membrane signaling events and a modulator of mTOR Complex 1 (mTORC1) activity in insects was upregulated. Overexpression of PIP4K in Drosophila upregulates the activity of mTORC1 and increases the levels of mTOR downstream [44,45]. We hypothesize that the increased levels of amino acids and substrates of TCA cycle in response to blood feeding in SCD1 KD mosquitoes, induces the expression of PIP4K to enhance the mTORC1 function [41].

The gene encoding the alpha subunit of the Na/K-transporting ATPase, which catalyzes the cation transfer across the membrane by converting ATP to ADP, was highly upregulated at 6h, 12h and 24h PBM further suggesting disruption of the epithelial cell membrane integrity and homeostasis.

Genes encoding odorant binding proteins OBP13 and OBP26 were markedly upregulated across all 5 time points in SCD1 KD mosquitoes. OBPs are shown to have strong affinity for SFAs like palmitate, stearate, octanoic acid and hexanoic acid [46]. Therefore, the accumulation of SFAs in SCD1 KD mosquitoes may trigger overexpression of these OBPs.

Replication, transcription and translation.

Significant differential regulation was observed at 24h PBM for 64% of genes (29/ 48) in this category, with 35% of them (18/48) exhibiting differential expression across all 5 time points. Nuclear export factor RNA Binding Motif protein 15 (RBM15) that is required for efficient mRNA export from the nucleus [47] was consistently downregulated in SCD1 KD mosquitoes across the time course. Cyclin A, an important regulator of cell cycle progression, was also downregulated at 24h PBM. The link between SCD1 and cell cycle regulators is largely unknown. However, recent data indicate that SCD1 blockage reduces AKT activity by dephosphorylation, with subsequent activation of glycogen synthase-kinase (GSK3-β), which in turn degrades cyclin D1 thereby inhibiting cell proliferation. AKT is a key component of the survival signaling pathway, and GSK3β is a downstream element of AKT signaling that is regulated by insulin [48]. Together these data suggest that SCD1 KD disrupts cell cycle progression in A. coluzzii.

Redox, apoptosis and detoxification.

A total of 37 genes in this category were differentially regulated in the SCD1 KD mosquitoes. The majority of them (83%) were downregulated, mostly at 24h PBM, and a third of them (12/37) were differentially regulated across all time points.

The detoxification gene GstD7, which plays an important role in protecting against cellular damage by controlling oxidative stress [49,50], was consistently upregulated in SCD1 KD mosquitoes across all time points. GstD7 is highly expressed in Drosophila suzukii upon exposure to the pesticide malathion [51]. Cytochrome c oxidase polypeptide 7A (COX7A), the terminal component of the mitochondrial respiratory chain that catalyzes the reduction of oxygen to water was downregulated across all time points. COX7A is shown to be downregulated during acute inflammation [52], supporting our hypothesis of a potent inflammatory response in SCD1 KD mosquitoes.

Genes encoding cytochrome P450s (CY314A1, CY302A1 CYP12F1, CYP6P1 & CYP9K1) were significantly downregulated at different time points. Previous studies demonstrate that cytochrome P450s play significant role in lipid homeostasis by assisting lipid degradation [53]. Hence, it is possible that scarcity of TAGs leads to the downregulation of cytochrome P450s as a feedback response. Glutathione-S-transferase D7 (GSTD7), a major detoxifying enzyme in mosquitoes, was upregulated, which may be due to the SFA:MUFA ratio distortion and represent a homeostatic response to hinder lipid peroxidation.

Transport/vesicle mediated transport.

As many as 26 genes involved in transport/ vesicle mediated transport system were differentially regulated in SCD1 KD compared to control mosquitoes, of which 15 were upregulated, mostly at 24h PBM. Trehalose transporter 1 (TRET1) that transports trehalose disaccharide from the fat body to hemolymph [54] was downregulated across the time course. Silencing trehalose transporter is shown to reduce trehalose content in hemolymph and shorten the lifespan of A. gambiae under stress conditions like desiccation and high temperature [55,56].

Metabolomic profile of SCD1 KD mosquitoes

We used GC-MS-based metabolomic analysis of deproteinized extracts to examine the metabolic changes induced in female mosquitoes prior to (0h PBM) and upon feeding on human blood (18h and 36h PBM) following SCD1 silencing or SA (in 1 mM Tris) treatment. Mosquitoes injected with dsLacZ or treated with 1 mM Tris were used as controls, respectively. Differences in metabolite signals between experimental and control groups were evaluated.

Desaturase indices, the FFA product to substrate, serves as a biomarker for the desaturase activity. The results showed that SCD1 KD causes a significant reduction in desaturase indices 16:1/16:0 (palmitoleate to palmitate) at 18h and 36h PBM (Fig 4A) and 18:1/18:0 (oleate to stearate) at all time points (Fig 4B). A statistically significant reduction in both 16:1/16:0 and 18:1/18:0 indices was also observed upon SA treatment albeit only at 36h PBM. These data confirm that silencing SCD1 or inhibition of the desaturase activity blocks the conversion of SFAs to UFAs. They are also in agreement with other studies correlating decreased SCD1 activity with increased SFA content, caused by administration of the cyclopropenoic fatty acids SA and malvalic acid [57].

[Figure omitted. See PDF.]

Fig 4. Metabolic responses in blood fed SCD1 KD and SA treated mosquitoes.

Desaturase indices (A) 16:1/16:0 (Palmitoleate to Palmitate) and (B) 18:1/18:0 (Oleate to Stearate) in SCD1 KD (left) and SA-treated (right) adult female mosquitoes prior to (0h) and 18h and 26h post bloodmeal (PBM). DsLacZ injected and 1mM Tris treated mosquitoes were used as controls for SCD1 KD and SA-treatment, respectively. Reported values are mean ± SE of 5 independent biological replicates, where n = 3 mosquitoes per group per replicate. *p<0.05, *p<0.001, ***p< 0.0001 using the student’s t-test. (C) Metabolite profiles in SCD1 KD and SA-treated mosquitoes and their respective dsLacZ and 1mM Tris treated controls. Only metabolites which showed statistically significant (*p<0.01) changes in the percentage of metabolite signals in the experimental group compared to the control group in at least one of the conditions are shown in the heat map. Red bars represent the mean percentage of metabolite signals in SCD1 KD or SA treated mosquitoes from 5 independent biological repeats. Black bars show mean percentage of metabolite signals in respective controls.

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

The metabolomic analysis of 43 aqueous metabolites showed significant differences (P<0.01 or less) of several metabolite concentrations between the SCD1 KD compared to control mosquitoes at one or more of the examined time points (Fig 4C). In SCD1 KD mosquitoes, the levels of substrates of the TCA cycle including citrate (36h PBM), malate (18h and 36h PBM) and fumarate (36h PBM) increased in response to blood feeding, while the intermediates of glycolysis, glycerol-3P (also 36h PBM), and the pentose phosphate pathway including ribose 5-P (36h PBM) decreased in response to blood feeding. The levels of D-glucose and D-glucose-6P remained unchanged.

Female mosquitoes largely convert protein-rich liquid meals into amino acids after blood ingestion. Consistent with this, only traces of most amino acids were observed at 0h PBM which significantly increased after bloodmeal in SCD1 KD mosquitoes at 36h PBM compared to that of the control. During oogenesis, phenylalanine is preferentially channeled towards late protein synthesis in order to meet the need for egg cuticle hardening in A. gambiae [16]. Interestingly, decreased levels of phenylalanine and–even more so–of its derivative tyrosine were observed in SCD1 KD mosquitoes at 18h and 36h PBM. It has been previously shown that interference with A. gambiae phenylalanine hydroxylase resulting in reduced tyrosine formation leads to fewer eggs and retarded vitellogenesis [16]. It can be therefore inferred that SCD1 KD leads to a perturbed conversion of phenylalanine to tyrosine and a subsequent suppression of egg development. Indeed, decreased levels of histidine may be linked to the role of histidine as a major intracellular buffering agent in animal cells [58]. Buffering agents maintain cellular pH and little changes in pH can cause either acidosis or alkalosis in the body with deleterious effects on metabolism. Hence, decreased levels of histidine may indicate imbalance in the pH system of SCD1 KD mosquitoes. Increased concentration of essential amino acids, such as isoleucine, valine, glutamic acid/glutamine, alanine, lysine and methionine may be the result of a feedback mechanism to improve amino acid uptake.

The levels of ornithine, a central component of the urea cycle, were higher in SCD1 KD at 18h and 36h PBM compared to control whereas the levels of putrescine were increased at 36h PBM.

Integrative metabolic pathway analysis

To better understand metabolic perturbations caused by SCD1 KD, 17 genes and 21 metabolites from the transcriptomic and metabolomic analyses, respectively, were mapped onto cellular metabolic processes (Fig 5). Pathway analysis provided insights into SCD1 KD-induced cellular toxicity dynamics in mosquitoes and revealed that silencing SCD1 significantly disrupts the TCA cycle and its interwind pathways with increasing accumulation of substrates and feedback regulation of expression of metabolic enzymes.

[Figure omitted. See PDF.]

Fig 5. Integrative analysis of metabolic pathways.

Metabolite and transcript abundance differences between dsSCD1 (KD) and dsLaZ (control) injected mosquitoes are mapped on reconstructed KEGG pathways encompassing glycolysis, gluconeogenesis, TCA cycle and Urea cycle. Transcript differences in any of five time points are indicated as upregulated (yellow) or downregulated (blue) when 2-fold or more expression changes are documented in SCD1 KD versus control mosquitoes (p<0.05 in ANOVA t-test). The metabolite data were normalized for each metabolite basis and plotted as circle fragments of which the radius is directly proportional to the highest value for the metabolite in any of the three time points and two conditions. AGL, α-glucosidase; AAM, α-amylase; PC, pyruvate carboxylase; PEPCK, phosphoenol pyruvate carboxykinase; LDH, lactate dehydrogenase; GALK, galactokinase; LDH, lactate dehydrogenase; FBPA, fructose-bisphosphate aldolase; SUCL-GDP, succinate-CoA ligase; SUCS-ADP, succinate-CoA synthetase; ACC, acetyl-CoA carboxylases; ASS, arginosuccinate synthase; ASL, arginosuccinate lyase; ODC, ornithine decarboxylase; CPT, carnitine-palmitoyl transferase; Glu-6P, glucose 6-phospate; Glu-1P, glucose 1-phospate; Gal-6P, galactose 6-phospate; Gal-1P, galactose 1-phospate; Rib-5P, ribose 5-phospate; PEP; phosphoenol pyruvate; Fru-6P, Fructose 6-phospate; Fru-1,6BP, fructose 1,6-biphospate; GAP, glyceraldehyde-3-phosphate; 2PG, 2-phosphoglycerate; 3PG, 2-phosphoglycerate; PEP, phosphoenol pyruvate.

https://doi.org/10.1371/journal.ppat.1009486.g005

Silencing SCD1 results an increased accumulation of SFA over time by modulating SFA/MUFA ratios that hinder acetyl-CoA carboxylases (ACC) activity through a well-known feedback mechanism resulting in a drop of malonyl-CoA that down-regulates carnitine-palmitoyl transferase (CPT) and in turn increases transport of acyl-CoA into mitochondria for β-oxidation [59]. This upregulated β-oxidation in turn increases the reaction rate of many of the steps in the TCA cycle, thereby increasing flux throughout the pathway and resulting in higher accretion of TCA cycle intermediates over time, i.e. citrate, succinate, fumarate and malate. The significantly higher levels of amino acids that can fuel the TCA cycle at various steps, such as methionine, isoleucine and valine, which is thought to be part of a feedback mechanism to increase amino acid uptake could indicate a major cause of accretion of these metabolites in SCD1 KD mosquitoes. Therefore, accumulation of TCA cycle intermediates is the likely cause of the decreased expression of rate limiting enzymes such as SUCS-GDP forming and SUCS-ADP forming enzymes, as well as those that feed into the TCA cycle including arginosuccinate synthase (ASS) and arginosuccinate lyase (ASL) observed in SCD1 KD mosquitoes in order to achieve homeostasis. The much higher levels of asparagine and glutamate may also be part of the feedback mechanism to suppress the TCA cycle through depletion of intermediate metabolites.

The low intracellular levels of the polyamine putrescine, produced through the urea cycle, in SCD1 KD compared to control mosquitoes could be a direct result of the low ASS and ASL levels. However, it may also be due to the parallel upregulation of ornithine decarboxylase (ODC), which converts ornithine to putrescine, and suggest increased metabolism of putrescine initially to spermidine, which is also low in SCD1 KD mosquitoes, and then though the glutathione detoxification pathways. Hence the levels of ornithine appear to be unaffected between SCD1 KD and control mosquitoes. These changes may also indicate a reduction in the speed of the urea cycle and possibly reduced rates of amino acid degradation.

Lower levels of sugar and downstream substrates of the glucose metabolism, such as ribose-5P (Rib-5P, which controls the rate-limiting purine synthesis), trehalose, 3-phosphoglycerate (3PG) and–to a lesser extent–glucose-6P (Glu-6P), in SCD1 KD mosquitoes compared to controls may indicate suppression of gluconeogenesis as a feedback mechanism, also suggested by the strong PEPCK downregulation. Interestingly, SCD1 KD mosquitoes appear to activate compensatory pathways, such as the Leloir gluconeogenic pathway, by upregulating galactokinase (GALK) to replenish the glucose flux. Downregulation of fructose-bisphosphate aldolase (FBPA) that converts fructose 1, 6-BP to glyceraldehyde-3-phosphate (GAP), both precursors of dihydroxyacetone phosphate (DHAP), cause low levels of glycerol-3P which would lead to reduced glycerol-based phospholipids, the main component of biological membranes.

In conclusion, the metabolic pathway analysis suggests that inhibition of SCD1 in female mosquitoes reprograms the interwinding pathways for cellular metabolism, leading to accumulation of various metabolites causing cytotoxicity and death.

SCD1 KD triggers a robust immune response

Immunity-related genes became increasingly and proportionally more numerous in the SCD1 KD upregulated gene set as time progressed (Fig 3 and S4 Table). A total of 58 genes were significantly differentially regulated in at least one of the time points, 83% of them upregulated and 14 genes differentially regulated across the time course. Immune induction was particularly prominent at 18h and 24h PBM.

One fifth of differentially regulated genes (19%) encode clip-domain serine proteases and serine protease homologs, including CLIPA5, CLIPA8, CLIPB11-12, CLIPB14-20 and CLIPB20; all but CLIPB20 upregulated. Many of these proteins were previously shown to be involved in microbial lysis and melanization in the mosquito hemolymph [60–63], suggesting activation of a systemic immune response in SCD1 KD mosquitoes. The melanization catalyzing enzyme PPO3 was also upregulated at 24h PBM.

TEP3, a member of the thioester-containing protein family, was amongst upregulated genes. TEP3 is paralogous to the main complement pathway effector, TEP1, but does not have a thioester motif and therefore is presumed to have a regulatory rather than effector function. Like TEP1, TEP3 binds to the putative receptor/adaptor complex LRIM1/APL1C and is involved in responses against bacteria and malaria parasites [64,65]. Another component of the complement pathway upregulated in SCD1 KD mosquitoes was CTLMA2, a member of the C-type lectin protein family. CTLMA2 is transcriptionally induced in response to infection and is thought to promote pathogen lysis and/or phagocytosis and negatively regulate melanization [66].

Fibrinogen related proteins, FREPs, play key roles in mosquito humoral and systemic responses and are thought to have complementary and synergistic functions, some seemingly acting as pattern recognition receptors (PRRs) [67]. Strong upregulation in SCD1 KD mosquitoes was observed for FREP19, FREP21, FREP24, FREP32, FREP59 and ficolin A. Another putative PRR, the scavenger receptor SCRBQ3, was also upregulated in SCD1 KD mosquitoes. SCRBQ3 is involved in melanization and induced in response to injury and associated opportunistic infections as well as upon injection of Sephadex beads [68,69].

The lipid-binding MD2-like proteins (MLs), ML4 and ML7, were strongly upregulated at 18h PBM in SCD1 KD. MLs are thought to bind lipids, and the founding member of the family, MD2, synergizes with TLR4 to trigger expression of pro-inflammatory cytokines in response to bacterial lipopolysaccharide (LPS). We hypothesize that mosquitoes upregulate ML4 and ML7 in response to abnormal and excessive saturated fatty acids, activating systemic reactions. Indeed, SFAs are shown to be recognized by the TLR4-MD2 complex, triggering trigger inflammatory pathways, similar to LPS [70].

There are two possible hypotheses for the observed upregulation of immune-related genes, most of which are involved in systemic immune reactions of the hemolymph, in SCD1 KD mosquitoes. The first is that over-activation of the immune system is caused by the gut microbiota (or their immune inducers) that gain increased access to immune receptors on the midgut epithelium or the hemolymph due to the loss of the peritrophic membrane and epithelium integrity. The second hypothesis is that the increased concentrations of SFAs and/or purinergic molecules such as ATP and ADP trigger an auto-inflammatory signaling. In invertebrates, humoral substances, often called danger associated molecular patterns (DAMPs), released by damaged cells or tissues, activate signaling pathways leading to inflammation.

We investigated the microbiota-induced response hypothesis by examining the expression of the antimicrobial peptide (AMP) gene Cecropin 1 (CEC1) in SCD1 KD and control (dsLacZ-injected) mosquitoes. Cecropins are some of the best characterized AMPs in insects, shown to provide first line defense against viruses, bacteria, and parasites. In the mosquito gut, CEC1 shows strong microbiota-dependent expression compared to other AMPs [15], and in our microarray experiments it was upregulated across all time points in SCD1 KD mosquitoes. Both, SCD1 KD and dsLacZ-injected mosquitoes were treated with antibiotics provided though the sugar meal since adult emergence to rid their gut from resident microbiota. As a control, mosquitoes were fed only on sugar. After 5 days, mosquitoes were either sampled (0 h) or fed on human blood though a membrane feeder and sampled 24 hours later. The load of bacteria in the mosquito gut was quantified by qPCR of the 16S rRNA at 0h PBM. Three independent biological replicates were performed.

The results showed that CEC1 was significantly upregulated both at 0h PBM and, more so, at 24h PBM in SCD1 KD compared to control mosquitoes (Fig 6). Surprisingly, the microbiota load prior to bloodmeal was significantly lower in SCD1 KD mosquitoes compared to controls. These results indicate that induction of CEC1 expression in SCD1 KD mosquitoes is independent of the microbiota load, which is decreased in SCD1 KD mosquitoes either due to immune activation or due to the physiological and metabolic conditions in the gut. Indeed, a much stronger upregulation in SCD1 KD mosquitoes was detected in mosquitoes treated with antibiotics. As expected, the bacterial load was substantially lower in antibiotic-treated compared to non-treated mosquitoes at 0h PBM, and even lower in SCD1 KD antibiotic-treated mosquitoes. These data reject the hypothesis that over-activation of the immune system in SCD1 KD mosquitoes is caused by the gut microbiota and suggest that these mosquitoes may suffer from an auto-inflammatory condition caused by high concentrations of SFAs, purinergic molecules and/or TCA cycle intermediates. They additionally suggest that the microbiota may be playing an inhibitory role to immune signaling in SCD1 KD mosquitoes, possibly triggered by the distorted SFA:MUFA ratio and mediated though microbiota tolerance pathways.

[Figure omitted. See PDF.]

Fig 6. Impact of gut microbiota on immune activation in SCD1 KD mosquitoes.

(A) Relative expression of the AMP CEC1 expression in dsLacZ-injected (C) and SCD1 KD (KD) adult mosquitoes fed on sugar (-) or sugar plus antibiotic (Abx; +) solution since emergence. After five, mosquitoes were sampled either without further treatment (0h PBM) or after they have been provided a blood meal (24h PBM). (B) Bacterial load in the midguts of mosquitoes treated as above, measured by qRT-PCR of 16S rRNA. Reported values are the mean ± SE of three independent biological experiments of 10 mosquitoes per group. ***p< 0.0001, *p<0.01.

https://doi.org/10.1371/journal.ppat.1009486.g006

Conclusion

Stearoyl-CoA desaturase is a rate limiting key enzyme regulating the ratio of SFAs to UFAs that is essential in the maintenance of cell membrane fluidity, proliferation, lipid accumulation and signaling. Its loss of function is known to alter lipid homeostasis due to the overloading with SFAs, inducing a cellular toxic response known as lipotoxicity, which eventually modulates the global metabolic profile and leads to various physiological abnormalities in animal cells [4]. Here we demonstrate that loss-of-function of SCD1 in adult female A. coluzzii mosquitoes, the main vector of human malaria in West and Central sub-Saharan Africa, leads, upon blood feeding, to a highly skewed SFA:UFA ratio, loss of cell membrane integrity, transcriptional and metabolic reprogramming and various other detrimental physiological effects, including an acute auto-inflammatory condition and compromised reproduction. We therefore conclude that SCD1 may be a good target of interventions to control vector abundance and block malaria transmission, including small molecule inhibitors (insecticides and drugs) and transgenesis methods aiming at mosquito population suppression or replacement. Indeed, targeting SCD1 specifically in the adult mosquito midgut after or at the time of the first blood meal could lead to a reduction of the mosquito lifespan without compromising reproduction in the first gonotrophic cycle, this providing a powerful means to impact on malaria transmission. SCD1 is highly conserved among insects with 1-to-1 orthologues found in all mosquito genomes sequenced to date, suggesting that any such intervention could be applicable to other mosquito vectors of human diseases.

Supporting information

S1 Fig. Sequence analysis of the A. coluzzii Stearoyl- CoA desaturase1.

Alignment of the central domain of SCD1 (amino acid residues 88 to 303) and the equivalent domain of SCD1 orthologues in human, Drosophila melanogaster and Plasmodium falciparum, performed using CLUTALW. The three His boxes designated as region Ia, Ib and II, respectively, and the eight conserved His residues are shown.

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

(DOCX)

S1 Table. Putative stearoyl-CoA desaturases in the A. gambiae genome.

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

(DOCX)

S2 Table. Putative orthologues of SCD1 in other mosquitoes and flies.

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

(DOCX)

S3 Table. List of genes showing altered expression across timepoints in SCD1 KD mosquitoes.

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

(DOCX)

S4 Table. List of immune genes showing differential expression in SCD1 KD mosquitoes.

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

(DOCX)

Acknowledgments

We thank Professor Mark S. Baird for kindly providing us with Sterculic Acid and Katarzyna Sala for her assistance with mosquito rearing.

Citation: Ferdous Z, Fuchs S, Behrends V, Trasanidis N, Waterhouse RM, Vlachou D, et al. (2021) Anopheles coluzzii stearoyl-CoA desaturase is essential for adult female survival and reproduction upon blood feeding. PLoS Pathog 17(5): e1009486. https://doi.org/10.1371/journal.ppat.1009486

References

1. Arrese EL, Soulages JL. Insect fat body: energy, metabolism, and regulation. Annu Rev Entomol. 2010;55:207–25. Epub 2009/09/04. pmid:19725772; PubMed Central PMCID: PMC3075550.

2. Clifton ME, Noriega FG. The fate of follicles after a blood meal is dependent on previtellogenic nutrition and juvenile hormone in Aedes aegypti. J Insect Physiol. 2012;58(7):1007–19. Epub 2012/05/26. pmid:22626792; PubMed Central PMCID: PMC3389259.

3. Lampe L, Jentzsch M, Kierszniowska S, Levashina EA. Metabolic balancing by miR-276 shapes the mosquito reproductive cycle and Plasmodium falciparum development. Nat Commun. 2019;10(1):5634. Epub 2019/12/12. pmid:31822677; PubMed Central PMCID: PMC6904670.

4. Alabaster A, Isoe J, Zhou G, Lee A, Murphy A, Day WA, et al. Deficiencies in acetyl-CoA carboxylase and fatty acid synthase 1 differentially affect eggshell formation and blood meal digestion in Aedes aegypti. Insect Biochem Mol Biol. 2011;41(12):946–55. Epub 2011/10/06. pmid:21971482; PubMed Central PMCID: PMC3210400.

5. Wang X, Hou Y, Saha TT, Pei G, Raikhel AS, Zou Z. Hormone and receptor interplay in the regulation of mosquito lipid metabolism. Proc Natl Acad Sci U S A. 2017;114(13):E2709–E18. Epub 2017/03/16. pmid:28292900; PubMed Central PMCID: PMC5380040.

6. Athenstaedt K, Daum G. The life cycle of neutral lipids: synthesis, storage and degradation. Cell Mol Life Sci. 2006;63(12):1355–69. Epub 2006/05/02. pmid:16649142.

7. 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(2):e1006853. Epub 2018/02/16. pmid:29447265; PubMed Central PMCID: PMC5814098.

8. Ntambi JM. Regulation of stearoyl-CoA desaturase by polyunsaturated fatty acids and cholesterol. J Lipid Res. 1999;40(9):1549–58. Epub 1999/09/14. pmid:10484602.

9. Phetsuksiri B, Jackson M, Scherman H, McNeil M, Besra GS, Baulard AR, et al. Unique mechanism of action of the thiourea drug isoxyl on Mycobacterium tuberculosis. J Biol Chem. 2003;278(52):53123–30. Epub 2003/10/16. pmid:14559907; PubMed Central PMCID: PMC4747054.

10. Pelaez R, Pariente A, Perez-Sala A, Larrayoz IM. Sterculic Acid: The Mechanisms of Action beyond Stearoyl-CoA Desaturase Inhibition and Therapeutic Opportunities in Human Diseases. Cells. 2020;9(1). Epub 2020/01/16. pmid:31936134; PubMed Central PMCID: PMC7016617.

11. Coetzee M, Hunt RH, Wilkerson R, Della Torre A, Coulibaly MB, Besansky NJ. Anopheles coluzzii and Anopheles amharicus, new members of the Anopheles gambiae complex. Zootaxa. 2013;3619:246–74. Epub 2013/01/01. pmid:26131476.

12. Blandin S, Moita LF, Kocher T, Wilm M, Kafatos FC, Levashina EA. Reverse genetics in the mosquito Anopheles gambiae: targeted disruption of the Defensin gene. EMBO Rep. 2002;3(9):852–6. Epub 2002/08/22. pmid:12189180; PubMed Central PMCID: PMC1084233.

13. Stathopoulos S, Neafsey DE, Lawniczak MK, Muskavitch MA, Christophides GK. Genetic dissection of Anopheles gambiae gut epithelial responses to Serratia marcescens. PLoS Pathog. 2014;10(3):e1003897. Epub 2014/03/08. pmid:24603764; PubMed Central PMCID: PMC3946313.

14. Habtewold T, Povelones M, Blagborough AM, Christophides GK. Transmission blocking immunity in the malaria non-vector mosquito Anopheles quadriannulatus species A. PLoS Pathog. 2008;4(5):e1000070. Epub 2008/05/24. pmid:18497855; PubMed Central PMCID: PMC2374904.

15. Rodgers FH, Gendrin M, Wyer CAS, Christophides GK. Microbiota-induced peritrophic matrix regulates midgut homeostasis and prevents systemic infection of malaria vector mosquitoes. PLoS Pathog. 2017;13(5):e1006391. Epub 2017/05/26. pmid:28545061; PubMed Central PMCID: PMC5448818.

16. Fuchs S, Behrends V, Bundy JG, Crisanti A, Nolan T. Phenylalanine metabolism regulates reproduction and parasite melanization in the malaria mosquito. PLoS One. 2014;9(1):e84865. Epub 2014/01/11. pmid:24409310; PubMed Central PMCID: PMC3883676.

17. Kind T, Wohlgemuth G, Lee DY, Lu Y, Palazoglu M, Shahbaz S, et al. FiehnLib: mass spectral and retention index libraries for metabolomics based on quadrupole and time-of-flight gas chromatography/mass spectrometry. Anal Chem. 2009;81(24):10038–48. Epub 2009/11/26. pmid:19928838; PubMed Central PMCID: PMC2805091.

18. Behrends V, Tredwell GD, Bundy JG. A software complement to AMDIS for processing GC-MS metabolomic data. Anal Biochem. 2011;415(2):206–8. Epub 2011/05/18. pmid:21575589.

19. Gullberg RC, Steel JJ, Pujari V, Rovnak J, Crick DC, Perera R. Stearoly-CoA desaturase 1 differentiates early and advanced dengue virus infections and determines virus particle infectivity. PLoS Pathog. 2018;14(8):e1007261. Epub 2018/08/18. pmid:30118512; PubMed Central PMCID: PMC6114894.

20. Jarsch IK, Daste F, Gallop JL. Membrane curvature in cell biology: An integration of molecular mechanisms. J Cell Biol. 2016;214(4):375–87. Epub 2016/08/17. pmid:27528656; PubMed Central PMCID: PMC4987295.

21. Sodja A, Fujioka H, Lemos FJ, Donnelly-Doman M, Jacobs-Lorena M. Induction of actin gene expression in the mosquito midgut by blood ingestion correlates with striking changes of cell shape. J Insect Physiol. 2007;53(8):833–9. Epub 2007/06/01. pmid:17537455; PubMed Central PMCID: PMC2756156.

22. Mendes AM, Awono-Ambene PH, Nsango SE, Cohuet A, Fontenille D, Kafatos FC, et al. Infection intensity-dependent responses of Anopheles gambiae to the African malaria parasite Plasmodium falciparum. Infect Immun. 2011;79(11):4708–15. Epub 2011/08/17. pmid:21844236; PubMed Central PMCID: PMC3257935.

23. Carta G, Murru E, Banni S, Manca C. Palmitic Acid: Physiological Role, Metabolism and Nutritional Implications. Front Physiol. 2017;8:902. Epub 2017/11/24. pmid:29167646; PubMed Central PMCID: PMC5682332.

24. Arnedo M, Latorre-Pellicer A, Lucia-Campos C, Gil-Salvador M, Antonanzas-Perez R, Gomez-Puertas P, et al. More Than One HMG-CoA Lyase: The Classical Mitochondrial Enzyme Plus the Peroxisomal and the Cytosolic Ones. Int J Mol Sci. 2019;20(24). Epub 2019/12/11. pmid:31817290; PubMed Central PMCID: PMC6941031.

25. Bock KW. The UDP-glycosyltransferase (UGT) superfamily expressed in humans, insects and plants: Animal-plant arms-race and co-evolution. Biochem Pharmacol. 2016;99:11–7. Epub 2015/10/11. pmid:26453144.

26. Huang FF, Chai CL, Zhang Z, Liu ZH, Dai FY, Lu C, et al. The UDP-glucosyltransferase multigene family in Bombyx mori. BMC Genomics. 2008;9:563. Epub 2008/11/29. pmid:19038024; PubMed Central PMCID: PMC2633020.

27. Mendes AM, Schlegelmilch T, Cohuet A, Awono-Ambene P, De Iorio M, Fontenille D, et al. Conserved mosquito/parasite interactions affect development of Plasmodium falciparum in Africa. PLoS Pathog. 2008;4(5):e1000069. Epub 2008/05/17. pmid:18483558; PubMed Central PMCID: PMC2373770.

28. Vlachou D, Schlegelmilch T, Christophides GK, Kafatos FC. Functional genomic analysis of midgut epithelial responses in Anopheles during Plasmodium invasion. Curr Biol. 2005;15(13):1185–95. Epub 2005/07/12. pmid:16005290.

29. Laurenti G, Tennant DA. Isocitrate dehydrogenase (IDH), succinate dehydrogenase (SDH), fumarate hydratase (FH): three players for one phenotype in cancer? Biochem Soc Trans. 2016;44(4):1111–6. Epub 2016/08/17. pmid:27528759.

30. Muller HM, Catteruccia F, Vizioli J, della Torre A, Crisanti A. Constitutive and blood meal-induced trypsin genes in Anopheles gambiae. Exp Parasitol. 1995;81(3):371–85. Epub 1995/11/01. pmid:7498434.

31. Schmidt C, Sciacovelli M, Frezza C. Fumarate hydratase in cancer: A multifaceted tumour suppressor. Semin Cell Dev Biol. 2020;98:15–25. Epub 2019/05/16. pmid:31085323; PubMed Central PMCID: PMC6974395.

32. Angelucci C, D’Alessio A, Iacopino F, Proietti G, Di Leone A, Masetti R, et al. Pivotal role of human stearoyl-CoA desaturases (SCD1 and 5) in breast cancer progression: oleic acid-based effect of SCD1 on cell migration and a novel pro-cell survival role for SCD5. Oncotarget. 2018;9(36):24364–80. Epub 2018/06/01. pmid:29849946; PubMed Central PMCID: PMC5966257.

33. Chen J, Aimanova KG, Gill SS. Aedes cadherin receptor that mediates Bacillus thuringiensis Cry11A toxicity is essential for mosquito development. PLoS Negl Trop Dis. 2020;14(2):e0007948. Epub 2020/02/06. pmid:32012156; PubMed Central PMCID: PMC7018227.

34. Hua G, Zhang R, Abdullah MA, Adang MJ. Anopheles gambiae cadherin AgCad1 binds the Cry4Ba toxin of Bacillus thuringiensis israelensis and a fragment of AgCad1 synergizes toxicity. Biochemistry. 2008;47(18):5101–10. Epub 2008/04/15. pmid:18407662; PubMed Central PMCID: PMC2662694.

35. Dana AN, Hong YS, Kern MK, Hillenmeyer ME, Harker BW, Lobo NF, et al. Gene expression patterns associated with blood-feeding in the malaria mosquito Anopheles gambiae. BMC Genomics. 2005;6:5. Epub 2005/01/18. pmid:15651988; PubMed Central PMCID: PMC546002.

36. Hegedus D, Erlandson M, Gillott C, Toprak U. New insights into peritrophic matrix synthesis, architecture, and function. Annu Rev Entomol. 2009;54:285–302. Epub 2008/12/11. pmid:19067633.

37. Noriega FG. Juvenile Hormone Biosynthesis in Insects: What Is New, What Do We Know, and What Questions Remain? Int Sch Res Notices. 2014;2014:967361. Epub 2014/01/01. pmid:27382622; PubMed Central PMCID: PMC4897325.

38. Tauchman SJ, Lorch JM, Orth AP, Goodman WG. Effects of stress on the hemolymph juvenile hormone binding protein titers of Manduca sexta. Insect Biochem Mol Biol. 2007;37(8):847–54. Epub 2007/07/14. pmid:17628283.

39. Brown MR, Clark KD, Gulia M, Zhao Z, Garczynski SF, Crim JW, et al. An insulin-like peptide regulates egg maturation and metabolism in the mosquito Aedes aegypti. Proc Natl Acad Sci U S A. 2008;105(15):5716–21. Epub 2008/04/09. pmid:18391205; PubMed Central PMCID: PMC2311378.

40. Ling L, Raikhel AS. Serotonin signaling regulates insulin-like peptides for growth, reproduction, and metabolism in the disease vector Aedes aegypti. Proc Natl Acad Sci U S A. 2018;115(42):E9822–E31. Epub 2018/10/03. pmid:30275337; PubMed Central PMCID: PMC6196551.

41. Sharma V, Pandey AK, Kumar A, Misra S, Gupta HPK, Gupta S, et al. Functional male accessory glands and fertility in Drosophila require novel ecdysone receptor. PLoS Genet. 2017;13(5):e1006788. Epub 2017/05/12. pmid:28493870; PubMed Central PMCID: PMC5444863.

42. Vogel KJ, Brown MR, Strand MR. Ovary ecdysteroidogenic hormone requires a receptor tyrosine kinase to activate egg formation in the mosquito Aedes aegypti. Proc Natl Acad Sci U S A. 2015;112(16):5057–62. Epub 2015/04/08. pmid:25848040; PubMed Central PMCID: PMC4413300.

43. Liu S, Li K, Gao Y, Liu X, Chen W, Ge W, et al. Antagonistic actions of juvenile hormone and 20-hydroxyecdysone within the ring gland determine developmental transitions in Drosophila. Proc Natl Acad Sci U S A. 2018;115(1):139–44. Epub 2017/12/20. pmid:29255055; PubMed Central PMCID: PMC5776822.

44. Gupta A, Toscano S, Trivedi D, Jones DR, Mathre S, Clarke JH, et al. Phosphatidylinositol 5-phosphate 4-kinase (PIP4K) regulates TOR signaling and cell growth during Drosophila development. Proc Natl Acad Sci U S A. 2013;110(15):5963–8. Epub 2013/03/27. pmid:23530222; PubMed Central PMCID: PMC3625332.

45. Sharma S, Mathre S, Ramya V, Shinde D, Raghu P. Phosphatidylinositol 5 Phosphate 4-Kinase Regulates Plasma-Membrane PIP3 Turnover and Insulin Signaling. Cell Rep. 2019;27(7):1979–90 e7. Epub 2019/05/16. pmid:31091438; PubMed Central PMCID: PMC6591132.

46. Wang J, Murphy EJ, Nix JC, Jones DNM. Aedes aegypti Odorant Binding Protein 22 selectively binds fatty acids through a conformational change in its C-terminal tail. Sci Rep. 2020;10(1):3300. Epub 2020/02/26. pmid:32094450; PubMed Central PMCID: PMC7039890.

47. Zolotukhin AS, Uranishi H, Lindtner S, Bear J, Pavlakis GN, Felber BK. Nuclear export factor RBM15 facilitates the access of DBP5 to mRNA. Nucleic Acids Res. 2009;37(21):7151–62. Epub 2009/09/30. pmid:19786495; PubMed Central PMCID: PMC2790900.

48. Sharma A, Nuss AB, Gulia-Nuss M. Insulin-Like Peptide Signaling in Mosquitoes: The Road Behind and the Road Ahead. Front Endocrinol (Lausanne). 2019;10:166. Epub 2019/04/16. pmid:30984106; PubMed Central PMCID: PMC6448002.

49. Misra JR, Horner MA, Lam G, Thummel CS. Transcriptional regulation of xenobiotic detoxification in Drosophila. Genes Dev. 2011;25(17):1796–806. Epub 2011/09/08. pmid:21896655; PubMed Central PMCID: PMC3175716.

50. Zanotto-Filho A, Dashnamoorthy R, Loranc E, de Souza LH, Moreira JC, Suresh U, et al. Combined Gene Expression and RNAi Screening to Identify Alkylation Damage Survival Pathways from Fly to Human. PLoS One. 2016;11(4):e0153970. Epub 2016/04/23. pmid:27100653; PubMed Central PMCID: PMC4839732.

51. Hamby KA, Kwok RS, Zalom FG, Chiu JC. Integrating circadian activity and gene expression profiles to predict chronotoxicity of Drosophila suzukii response to insecticides. PLoS One. 2013;8(7):e68472. Epub 2013/07/19. pmid:23861907; PubMed Central PMCID: PMC3702611.

52. Huttemann M, Helling S, Sanderson TH, Sinkler C, Samavati L, Mahapatra G, et al. Regulation of mitochondrial respiration and apoptosis through cell signaling: cytochrome c oxidase and cytochrome c in ischemia/reperfusion injury and inflammation. Biochim Biophys Acta. 2012;1817(4):598–609. Epub 2011/07/21. pmid:21771582; PubMed Central PMCID: PMC3229836.

53. Veith A, Moorthy B. Role of Cytochrome P450s in the Generation and Metabolism of Reactive Oxygen Species. Curr Opin Toxicol. 2018;7:44–51. Epub 2018/03/13. pmid:29527583; PubMed Central PMCID: PMC5841237.

54. Kikawada T, Saito A, Kanamori Y, Nakahara Y, Iwata K, Tanaka D, et al. Trehalose transporter 1, a facilitated and high-capacity trehalose transporter, allows exogenous trehalose uptake into cells. Proc Natl Acad Sci U S A. 2007;104(28):11585–90. Epub 2007/07/04. pmid:17606922; PubMed Central PMCID: PMC1905927.

55. Hagan RW, Didion EM, Rosselot AE, Holmes CJ, Siler SC, Rosendale AJ, et al. Dehydration prompts increased activity and blood feeding by mosquitoes. Sci Rep. 2018;8(1):6804. Epub 2018/05/03. pmid:29717151; PubMed Central PMCID: PMC5931509.

56. Liu K, Dong Y, Huang Y, Rasgon JL, Agre P. Impact of trehalose transporter knockdown on Anopheles gambiae stress adaptation and susceptibility to Plasmodium falciparum infection. Proc Natl Acad Sci U S A. 2013;110(43):17504–9. Epub 2013/10/09. pmid:24101462; PubMed Central PMCID: PMC3808592.

57. Ortinau LC, Pickering RT, Nickelson KJ, Stromsdorfer KL, Naik CY, Haynes RA, et al. Sterculic Oil, a Natural SCD1 Inhibitor, Improves Glucose Tolerance in Obese ob/ob Mice. ISRN Endocrinol. 2012;2012:947323. Epub 2012/12/05. pmid:23209931; PubMed Central PMCID: PMC3504409.

58. Sasaki S, Kobayashi M, Futagi Y, Ogura J, Yamaguchi H, Iseki K. Involvement of Histidine Residue His382 in pH Regulation of MCT4 Activity. PLoS One. 2014;10(4):e0122738. Epub 2015/04/29. pmid:25919709; PubMed Central PMCID: PMC4412719.

59. Imai N, Cohen DE. Trimming the Fat: Acetyl-CoA Carboxylase Inhibition for the Management of NAFLD. Hepatology. 2018;68(6):2062–5. Epub 2018/08/05. pmid:30076622.

60. Nakhleh J, Christophides GK, Osta MA. The serine protease homolog CLIPA14 modulates the intensity of the immune response in the mosquito Anopheles gambiae. J Biol Chem. 2017;292(44):18217–26. Epub 2017/09/21. pmid:28928218; PubMed Central PMCID: PMC5672044.

61. Schnitger AK, Kafatos FC, Osta MA. The melanization reaction is not required for survival of Anopheles gambiae mosquitoes after bacterial infections. J Biol Chem. 2007;282(30):21884–8. Epub 2007/06/01. pmid:17537726.

62. Volz J, Osta MA, Kafatos FC, Muller HM. The roles of two clip domain serine proteases in innate immune responses of the malaria vector Anopheles gambiae. J Biol Chem. 2005;280(48):40161–8. Epub 2005/09/29. pmid:16188883.

63. Yassine H, Kamareddine L, Osta MA. The mosquito melanization response is implicated in defense against the entomopathogenic fungus Beauveria bassiana. PLoS Pathog. 2012;8(11):e1003029. Epub 2012/11/21. pmid:23166497; PubMed Central PMCID: PMC3499577.

64. Povelones M, Upton LM, Sala KA, Christophides GK. Structure-function analysis of the Anopheles gambiae LRIM1/APL1C complex and its interaction with complement C3-like protein TEP1. PLoS Pathog. 2011;7(4):e1002023. Epub 2011/05/03. pmid:21533217; PubMed Central PMCID: PMC3077365.

65. Yan Y, Hillyer JF. Complement-like proteins TEP1, TEP3 and TEP4 are positive regulators of periostial hemocyte aggregation in the mosquito Anopheles gambiae. Insect Biochem Mol Biol. 2019;107:1–9. Epub 2019/01/29. pmid:30690067.

66. Schnitger AK, Yassine H, Kafatos FC, Osta MA. Two C-type lectins cooperate to defend Anopheles gambiae against Gram-negative bacteria. J Biol Chem. 2009;284(26):17616–24. Epub 2009/04/22. pmid:19380589; PubMed Central PMCID: PMC2719400.

67. Dong Y, Dimopoulos G. Anopheles fibrinogen-related proteins provide expanded pattern recognition capacity against bacteria and malaria parasites. J Biol Chem. 2009;284(15):9835–44. Epub 2009/02/06. pmid:19193639; PubMed Central PMCID: PMC2665105.

68. Gonzalez-Lazaro M, Dinglasan RR, Hernandez-Hernandez Fde L, Rodriguez MH, Laclaustra M, Jacobs-Lorena M, et al. Anopheles gambiae Croquemort SCRBQ2, expression profile in the mosquito and its potential interaction with the malaria parasite Plasmodium berghei. Insect Biochem Mol Biol. 2009;39(5–6):395–402. Epub 2009/04/16. pmid:19366631; PubMed Central PMCID: PMC4138513.

69. Warr E, Lambrechts L, Koella JC, Bourgouin C, Dimopoulos G. Anopheles gambiae immune responses to Sephadex beads: involvement of anti-Plasmodium factors in regulating melanization. Insect Biochem Mol Biol. 2006;36(10):769–78. Epub 2006/10/10. pmid:17027843.

70. Rocha DM, Caldas AP, Oliveira LL, Bressan J, Hermsdorff HH. Saturated fatty acids trigger TLR4-mediated inflammatory response. Atherosclerosis. 2016;244:211–5. Epub 2015/12/22. pmid:26687466.