Abstract
Lipid metabolism is an essential component in reproductive physiology. While lipid mobilization has been implicated in the growth of Plasmodium falciparum malaria parasites in their Anopheles vectors, the role of this process in the reproductive biology of these mosquitoes remains elusive. Here, we show that impairing lipolysis in Anopheles gambiae, the major malaria vectors, leads to embryonic lethality. Embryos derived from females in which we silenced the triglyceride lipase AgTL2 or the lipid storage droplet AgLSD1 develop normally during early embryogenesis but fail to hatch due to severely impaired metabolism. Embryonic lethality is efficiently recapitulated by exposing adult females to broad-spectrum lipase inhibitors prior to blood feeding, unveiling lipolysis as a potential target for inducing mosquito sterility. Our findings provide mechanistic insights into the importance of maternal lipid mobilization in embryonic health that may inform studies on human reproduction.
Citation: Itoe MA, Shaw WR, Stryapunina I, Vidoudez C, Peng D, Du EW, et al. (2024) Maternal lipid mobilization is essential for embryonic development in the malaria vector Anopheles gambiae. PLoS Biol 22(12):
e3002960.
https://doi.org/10.1371/journal.pbio.3002960
Academic Editor: Louis Lambrechts, Institut Pasteur, FRANCE
Received: April 29, 2024; Accepted: November 29, 2024; Published: December 17, 2024
Copyright: © 2024 Itoe 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: The authors confirm that all data underlying the findings are fully available without restriction. The numerical data are accessible from the Harvard Dataverse online repository using the link https://doi.org/10.7910/DVN/ULTW1K. Metabolomics and Lipidomics data are accessible from the EMBL-EBI MetaboLights online database using the link https://www.ebi.ac.uk/metabolights/MTBLS9881. The RNAseq data are accessible from the Gene Expression Omnibus (GEO) online database using the link https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE263712.
Funding: F.C. is funded by Howard Hughes Medical Institute (HHMI, www.hhmi.org) as an HHMI Investigator and by the National Institute of Health (NIH) grants R01AI148646, R01AI153404. M.A.I. was funded by the Charles A. King Trust postdoctoral research fellowship in Basic Sciences from Health Resources in Action (HRiA, www.hria.org). I.S. was funded by Natural Sciences and Engineering Research Council of Canada (NSERC, www.nserc-crsng.gc.ca) as a scholarship recipient (PGSD3 – 545866 – 2020). The funders were not involved in study design, data collection, analysis, or interpretation.
Competing interests: I have read the journal’s policy and the authors of this manuscript have the following competing interest: A patent application covering the concept of the application of lipase inhibitors to induce sterility in mosquitoes is being processed on behalf of F.C. and M.A.I. by the President and Fellows of Harvard University. The authors state that they have no other competing interests.
Abbreviations:
ACTH,
adrenocorticotropic hormone; AKH,
adipokinetic hormone; ATGL,
adipose triglyceride lipase; SA,
bovine serum albumin; cAMP,
cyclic adenosine monophosphate; CE,
cholesterol ester; DAG,
diacylglyceride; FDR,
false discovery rate; HSL,
hormone-sensitive lipase; LD,
lipid droplet; LPE,
lysophosphatidylethanolamine; MB,
maleate buffer; MPA,
mobile phase A; MPB,
mobile phase B; PBM,
post blood meal; PBS,
phosphate-buffered saline; PFA,
paraformaldehyde; PKA,
protein kinase A; TAG,
triglyceride; TCA,
tricarboxylic acid; TEM,
transmission electron microscopy
Introduction
Reproduction is a highly energy-demanding process that requires tight metabolic control coordinated by hormones, metabolites, and immune factors [1,2]. In mammals, lipid stores in the form of triglycerides (TAGs) are deposited in the cells of adipose tissue, where they serve as a dynamic energy reservoir. As TAGs are unable to cross cell membranes, for transport to other tissues, they need to be either incorporated into vesicles or remobilized. Regulation of these processes is critically important during gestation as they coordinate the necessary energy supply to the growing fetus, and their perturbations can have negative effects on pregnancy outcomes [3,4]. In reproductive organs, lipid mobilization allows essential reproductive functions such as the development of ovarian follicles and the synthesis of steroid hormones, which prepare the endometrium for implantation [3].
In simplified terms (reviewed in [5]), TAG remobilization follows a hydrolytic process requiring their breakdown into free fatty acids (FFAs) and glycerol, regulated by hormonal stimuli that activate lipolytic enzymes like hormone-sensitive lipase (HSL) and adipose triglyceride lipase (ATGL) [6–8]. The canonical pathway for lipolysis activation in adipocytes is driven by catecholamines, adrenocorticotropic hormone (ACTH), and secretin, which upon binding to β-adrenergic G-protein-coupled receptors stimulate cyclic adenosine monophosphate (cAMP) synthesis, activating protein kinase A (PKA). PKA activation causes the phosphorylation of HSL, which induces its translocation onto the lipid droplet (LD) surface aided by PKA-mediated phosphorylation of perilipin-1 [9,10]. Perilipin-1 phosphorylation also results in the dissociation of the co-activator comparative gene identification-58 (CGI-58) [11,12], which can then bind to ATGL and activate lipolysis [13]. Upon completion of hydrolysis, FFAs will bind to carrier proteins that transport them to target tissues in the body for their use in energy production through ß-oxidation, hormone synthesis, or other cellular metabolic needs [14].
Although the last common ancestor of insects and mammals is estimated to have lived over 500 million years ago, insects possess a similar metabolic cascade to mobilize lipids. Lipolysis of TAGs is hormonally stimulated by a neuropeptide hormone, adipokinetic hormone (AKH), that has an analogous role to glucagon in mammals and is released from the corpora cardiaca [15,16]. Once released, AKH binds to its G protein-coupled receptor (AKHR) triggering an increase in intracellular calcium levels and cAMP synthesis, followed by PKA activation [17]. Activated PKA in turn phosphorylates lipid storage droplet protein 1 (LSD1), a homologue of perilipin 1-like protein, on the surface of LDs, ultimately inducing the activation of TAG lipases to catalyze the rate-limiting step in the hydrolysis of fatty acids from TAGs to generate FFA and diacylglycerides (DAGs) [18–22]. In hematophagous insects, after blood feeding, complex lipids such as TAGs and phospholipids are broken down into FFAs in the midgut lumen for transport across the luminal membranes [23], where they are reutilized to synthesize DAGs, TAGs, esters, and phospholipids [24]. These lipids are then loaded—in concert with lipid transfer proteins—onto circulating lipoproteins, principally lipophorin (Lp), on the basal side of the midgut epithelium for delivery to key tissues such as flight muscles and ovaries, or directed to the fat body for storage or further mobilization [19,20,25,26]. Dynamic changes occur in the lipid profile of the fat body during a reproductive cycle [27], with substantial proportions of stored lipids being mobilized for egg production [28,29]. Lp-transported lipids, as well as the yolk protein vitellogenin (Vg), are incorporated into developing oocytes [30,31] whereupon lipid is organized into LDs and Vg crystallizes into vitellin yolk granules [32] to be degraded after oviposition to nourish the developing embryo [33]. Coordinated lipolysis regulation is mediated by multiple conserved signaling pathways, including 20-hydroxyecdysone [34,35], Hedgehog [36], Decapentaplegic [37], NF-κB [38], as well as the neuropeptides FMRFamide [39] and AKH [40]. The components of the lipolytic machinery are shared by important insect vectors of human diseases including malaria, Chagas disease, and African Trypanosomiasis [23,41–44], and could potentially be targeted to reduce the size of vector populations and aid in disease control efforts. Orthologues of TAG lipase and AKHR are significantly increased after blood feeding in the fat body of Aedes aegypti mosquitoes [34,40], and TAG stores in the fat body, which are sharply reduced after blood feeding due to their mobilization and incorporation into developing eggs, are not mobilized when the AKH receptor AKHR is silenced by RNAi [40]. Similarly, disruption of lipid mobilization through targeting AKH signaling, ATGL or Perilipin/LSD1 in several insects decreases fecundity [43,45] and causes embryonic lethality [21,46]. However, information regarding lipid mobilization in Anopheles is relatively scarce.
In one of the most important vectors of human malaria, the Anopheles gambiae mosquito, lipid transport by Lp has been shown to impact the development of Plasmodium falciparum and Plasmodium berghei, respectively, a human and rodent malaria parasite. Indeed, silencing Lp before an infectious blood meal reduces oocyst numbers in both parasite species [47–49], while it impairs oocyst growth and sporozoite infectivity in P. berghei but not P. falciparum [48–50]. These studies also revealed that Lp, and to a lesser extent Vg, which transports approximately 5% of lipids within oocytes [51], are essential for egg development as their silencing drastically reduces the number of developing ovarian oocytes [31,47,48,50]. A role of upstream lipolytic pathways in the reproduction of this or other human malaria vectors has yet to be determined, and their function has only previously been studied in the context of pathogen infections. Silencing of the TAG lipase AgTL2 (a pancreatic lipase that is the functional homologue of ATGL and predominantly induced in the fat body after blood feeding [52]) by RNAi leads to accelerated P. falciparum oocyst growth, possibly due to the accumulation of TAGs, DAGs, and phosphatidylcholine observed in midguts and fat bodies [48], while depleting AKH1 and AKHR leads to reduced number of sporozoites in the salivary glands [53]. Whether impairing lipolysis impacts reproduction in this anopheline is unknown, and addressing this question could uncover key reproductive targets to reduce mosquito populations and curb malaria incidence in disease endemic countries.
Here, we demonstrate a key role of the An. gambiae maternal lipolytic machinery in shaping the development and survival of progeny. Targeted depletion of AgTL2 in mothers induces considerable perturbations to energy metabolism in the progeny, preventing embryos from hatching and causing embryonic lethality. Treating adult females with orlistat, a broad-spectrum hydrolase inhibitor, mimics the effects of AgTL2 silencing, causing dose-dependent embryonic mortality without, however, accelerating P. falciparum development. These data provide a proof of principle that targeting lipolysis can aid mosquito control and increase our understanding of the importance of lipid metabolism during pregnancy and embryonic development.
Results
Maternal triglyceride (TAG) lipolysis is critical for the survival of An. gambiae embryos
To determine the role of TAG mobilization in An. gambiae reproduction, we silenced the expression of the TAG lipase AgTL2 and found that its depletion reduced the number of eggs laid by females compared to controls (Fig 1A and 1B). Silencing of AgTL2, which is induced after blood feeding predominantly in the fat body (S1A and S1B Fig) [52], led to an increase in glyceride levels in midguts and fat bodies after blood feeding, mostly evident at 48 hours (h) post blood meal (PBM), paralleled by a large decrease of glycerides in ovaries, where accumulation was severely impaired at both 24 h and 48 h PBM (Fig 1C). This effect was confirmed by reduced staining of neutral lipids in ovaries of AgTL2-depleted females, while midguts and fat bodies were laden with large LDs not observed in controls (S1D Fig).
Fig 1. Impairing lipolysis in An. gambiae females induces severe embryonic lethality.
(A) Schematic representation of the experimental setup, showing that females are first injected with dsAgTL2 or control dsGFP, then mated and blood fed, and then analyzed for number of eggs and survival of the progeny. d = day. (B, C) Females injected with dsAgTL2 (B) laid fewer eggs after blood feeding (Mann–Whitney test) and (C) have increased glyceride levels in midguts and fat body and decreased levels in ovaries, compared to dsGFP-injected controls. hPBM = hours post blood meal; NBF = non-blood fed (Least square means model, followed by FDR-corrected t test). (D) Glyceride levels are significantly reduced in embryos from AgTL2-depleted females throughout development, compared to control dsGFP-injected group (Ordinary 2-way ANOVA, Šídák’s multiple comparisons correction). (E) Eggs oviposited by AgTL2-silenced females have significantly impaired embryo survival and produce very few larvae (Mann–Whitney test). (F, G) Injection of dsAgLSD1 (F) significantly reduces the number of eggs laid and (G) severely decreases survival rates (Mann–Whitney test). n = number of mosquitoes (except for glyceride assay where n is the number of pools of tissues, including 3 insects each); 3–4 biological replicates for A–E, 2 biological replicates for F and G. Not all significant comparisons are shown, for clarity. Numerical data supporting this figure is available on the Harvard Dataverse online repository at https://doi.org/10.7910/DVN/ULTW1K.
Embryos from AgTL2-depleted mothers had a strong reduction in glyceride levels (Fig 1D), as expected given the low levels observed in oocytes. Strikingly, only a handful of larvae emerged from several egg batches oviposited by these females (Fig 1E). We confirmed that these effects were indeed due to impaired TAG lipolysis by silencing the positive regulator AgLSD1 (S1C Fig), which similarly reduced the number of eggs laid by females and induced almost complete infertility (Fig 1F and 1G). Embryonic development appeared comparable through gastrulation, differentiation, and into late embryogenesis in both groups (Fig 2A), but embryos derived from AgTL2-depleted females (henceforth also referred to as dsAgTL2 embryos) failed to hatch. Further characterization by transmission electron microscopy (TEM) confirmed normal morphology at early stages, while at later time points the gut of dsAgTL2 embryos was filled with large vesicles not observed in controls (Fig 2B and 2C). Altogether, these data demonstrate that the maternal lipolytic machinery is essential for embryonic development in these malaria vectors.
Fig 2. Progeny of AgTL2-deficient females appear to develop normally during embryogenesis.
(A) Nile Red (neutral lipids) and DAPI (DNA) staining show apparent normal development in embryos from AgTL2-knockdown and control females at 4–7, 21–24, 35–38 h after oviposition (scale bar = 50 μm). (B) TEM reveals normal granular content at early time points in guts of embryos from AgTL2-knockdown females (scale bar = 2 μm). (C) Brightfield (lower panel) and TEM (upper panel) reveals large granules in the midgut of embryos from AgTL2-depleted females at 38 h after oviposition (Brightfield, scale bar = 50 μm; TEM, scale bar = 2 μm). Representative images from 3 experiments are shown. TEM, transmission electron microscopy.
Key metabolic transitions are altered in embryos from AgTL2-deficient females
To reconstruct a global view of the changes induced in embryos by TAG lipase silencing in mothers, we performed multiomics analyses during embryogenesis (at 6 h, 24 h, and 38 h after oviposition), collecting lipidomics, metabolomics and RNAseq data sets. dsAgTL2 embryos inherited significantly reduced levels of maternal lipids than controls, driven mainly by reduced TAGs (Fig 3A and S1 Data), consistent with our glyceride quantifications, while cholesterol ester (CE) and lysophosphatidylethanolamine (LPE) levels were increased (Figs 3A, S2A, and S2B). The changes in TAG and CE levels were not specific for a particular fatty acid composition (S2C and S2D Fig).
Fig 3. AgTL2 embryos have altered lipids, metabolites, and transcriptional profiles during development.
(A) LC-MS analysis shows reduced levels of total lipids and TAGs and increased levels of CEs in embryos from AgTL2-depleted females compared to controls throughout development (Ordinary 2-way ANOVA, Šídák’s multiple comparisons correction). (B) Principal component analysis shows distinct clustering of metabolites in AgTL2 embryos relative to controls. Small circles represent individual replicates. Numbers indicate time after oviposition. (C) Heatmap of metabolites analyzed by LC-MS showing levels of the 60 most highly dysregulated metabolites in AgTL2 embryos compared to controls throughout development (values of t test statistic, range: blue to red = significant decrease to increase) (3 biological replicates represented by triplicate columns at each time point). (D) Volcano plots of RNAseq data of embryos from control and AgTL2-depleted females at 4–6 h (6 h), 22–24 h (24 h), and 36–38 h (38 h) after oviposition (4 biological replicates). Few genes, including FAS, ABC transporter, trehalose transporter, and tRNA-Ser, are altered at early time points, while 785 genes are up-regulated and 220 genes are down-regulated at the final time point. At 38 h, top 5 up-regulated genes (log2 FC ≥ 2.0) are: 1: AGAP007650 (Growth Arrest and DNA-damage-inducible protein); 2: AGAP001610 (unknown); 3: AGAP029366 (Retrotrans-gag domain-containing protein); 4: AGAP009460 (c-Jun N-terminal Kinase JNK); 5: AGAP028615 (unknown). Top 5 down-regulated genes (log2 FC ≤ −2.0) are: 6: AGAP009591 (formyltetrahydrofolate dehydrogenase); 7: AGAP003620 (unknown); 8: AGAP007903 (Excitatory amino acid transporter 2 isoform 2); 9: AGAP009896 (Proton-coupled amino acid transporter); and 10: AGAP001742 (unknown). See Tables 1 and S1 and S2 Data and S3 Data. Numerical data supporting this figure is available in the Harvard Dataverse online repository at https://doi.org/10.7910/DVN/ULTW1K. CE, cholesterol ester; FAS, fatty acid synthase; TAG, triglyceride.
Additionally, embryos from AgTL2-depleted mothers showed highly different metabolic profiles (Fig 3B) indicative of substantial breakdown in key metabolic processes. Indeed, while in controls we observed a steady increase in major metabolites over time—including sugars (trehalose, lactose, and glucose), intermediates of the tricarboxylic acid (TCA) cycle and ß-oxidation, and nucleotides and their derivatives—dsAgTL2 embryos had significantly reduced levels of most of these metabolites at one or more time points, while hallmarks of purine and pyrimidine catabolism (ureidopropionic acid, hypoxanthine, xanthine, and uric acid) were notably increased at 38 h post oviposition (Figs 3C and S3 and S2 Data).
While changes in lipid and metabolite composition started early during embryogenesis, we observed limited differences in the transcriptional profile of the 2 groups at the 2 earlier time points (Fig 3D and S3 Data). One exception was the up-regulation of genes involved in fatty acid synthesis in dsAgTL2 embryos at 24 h post oviposition, likely to compensate for the severe energetic deficits derived from the reduction in fatty acids and ß-oxidation (Table 1 and S4A and S4B Fig). However, by 38 h post oviposition, there were large differences in the dsAgTL2 group characterized by hundreds of up- and down-regulated genes. Among the down-regulated genes, we observed a large decrease in neuropeptides and G-protein coupled receptors involved in sensory perceptions (Table 1 and S4A, S4B, and S5 Figs). Indeed, corazonin and eclosion hormone, 2 neuropeptides produced by neurosecretory cells in the brain that regulate ecdysis [54,55], and peptide hormones such as cchamide-1-related, secreted by endocrine cells localized in the adult gut and possibly regulating gut-brain communication [56,57], were all substantially down-regulated (Table 1 and S4B Fig). Several cuticular protein genes were also suppressed, suggesting inhibition of the second cuticle that is usually formed before hatching in Diptera (58). Up-regulated genes showed a strong signature of stress response and protein folding, including heat shock proteins, growth arrest and DNA-damage inducible protein, and a c-Jun N-terminal Kinase (JNK), possibly indicative of cellular death (Table 1 and S4A and S4B Fig).
These data reveal that maternal TAGs are essential for driving key metabolic processes during embryogenesis, such that when these lipids are limiting, embryos cannot complete development. A summary of the major metabolic changes and how they may be connected is provided in S5 Fig.
A hydrolase inhibitor induces severe embryonic lethality in An. gambiae
The remarkable phenotype observed when silencing AgTL2 and AgLSD1 makes TAG hydrolysis an appealing target for reducing the size of mosquito populations in field settings. Indeed, strategies that rely on chemically interfering with the reproductive biology of Anopheles mosquitoes are already in use in malaria endemic regions to mitigate the negative effects of widespread insecticide resistance, as in the case of insecticide-treated bed nets also containing the growth regulator pyriproxyfen [59–61], a juvenile hormone analogue that reduces the mosquito reproductive output [62]. We hypothesized that using lipase inhibitors may also induce embryonic lethality when administered to adult females before blood feeding. Consistent with our hypothesis, injections of orlistat—a broad-spectrum hydrolase inhibitor that covalently binds to Ser152 in the active site of mammalian lipases [63]—into the hemolymph of females prior to blood feeding significantly impaired egg development and induced strong embryonic lethality in a concentration-dependent manner (Fig 4A–4C).
Fig 4. Treatment with orlistat induces severe embryonic mortality but does not affect P. falciparum growth.
(A) Experimental scheme showing how 4-day-old adult An. gambiae females are injected with control (DMSO) or orlistat and then blood fed with either infected or uninfected blood. (B, C) Orlistat treatment decreases both (B) the number of eggs laid and (C) the survival of the embryos in a dose-dependent manner (4 biological replicates). (D) Experimental scheme showing how An. gambiae females are exposed for 1 h to orlistat-treated glass or control plates before blood feeding to measure oviposition and fertility. (E, F) Tarsal exposure to orlistat (E) decreases the number of eggs laid and (F) reduces embryo survival (2 biological replicates). (G) Orlistat treatment has no effect on P. falciparum oocyst size (two-tailed Unpaired t test) (2 biological replicates). d = day; n = number of individual mosquitoes; pi = post injection. (Ordinary 1-way ANOVA (B and E) and Kruskal–Wallis test (C and F), Dunn’s multiple comparisons correction). Numerical data supporting this figure is available in the Harvard Dataverse online repository at https://doi.org/10.7910/DVN/ULTW1K.
In field settings, mosquitoes can absorb insecticides and sterilizing compounds via their legs when they land on bed nets incorporating these ingredients, as in the case of pyriproxyfen nets [59–62]. When we tested this method of delivery by allowing females to rest on a surface coated with orlistat before blood feeding, their reproductive output was also strongly impaired in a dose-dependent manner (Fig 4D–4F). These data reveal that the development of Anopheles-specific lipase inhibitors may represent a promising tool for field interventions aimed at reducing the size of field mosquito populations.
Orlistat does not affect P. falciparum development in the mosquito
Our previous work had shown that the growth of P. falciparum oocysts is accelerated in females depleted for AgTL2, leading to early sporozoite invasion of the salivary glands [48]. This effect would have obvious negative consequence when using lipase inhibitors to induce sterility in field populations, as it would lead to mosquitoes that can better support parasite transmission. We first confirmed that TAG lipolysis affects parasite growth in An. gambiae by silencing AgLSD1 prior to infection with P. falciparum. Paralleling results previously observed in AgTL2-silenced females, AgLSD1 depletion accelerated oocyst growth without affecting prevalence or intensity of infection (S6A and S6B Fig), leading to increased sporozoite numbers at early time points (12 to 13 days) after infection relative to controls (S6C Fig). Orlistat treatment, however, did not induce an increase in parasite growth rates, and oocyst size was indistinguishable from controls (Fig 4G), providing reassurance that this strategy for inducing sterility would not have unwanted consequences on parasite transmission. Interestingly, orlistat-treated mosquitoes had elevated glyceride levels in the midgut but not in the fat body (S6D Fig), as instead observed in AgTL2-silenced females, possibly explaining the lack of effects on oocyst growth and showing this inhibitor may preferentially affect midgut lipases rather than AgTL2. When combined, these data further support the development of Anopheles-specific lipase inhibitors for controlling mosquito populations that transmit malaria.
Discussion
Our data reveal an unappreciated role for maternal TAG remobilization in shaping the development of An. gambiae embryos. We show that impairing lipolysis via AgTL2 depletion in the fat body, the fat storage tissue, leads to a substantial decrease in TAGs in the ovaries of adult females, the destination of lipids shuttled via the function of lipoproteins. In turn, we demonstrate that inheriting vastly reduced TAG levels triggers a cascade of nutritional dysregulation that prevents embryos from hatching and leads to their death. The almost complete penetrance of this phenotype led us to test the possibility of targeting lipolysis in mosquitoes to reduce the size of field populations. We were pleasantly surprised by the efficacy in proof of principle experiments of the lipase inhibitor orlistat, which induced significant embryonic lethality when applied to females. Orlistat is a drug used in humans to counter obesity and therefore it would not be suitable for a mosquito-targeting strategy. However, in the future, it may be feasible to generate Anopheles-specific inhibitors that target the function of AgTL2, AgLSD1, or other factors critical to mosquito lipolysis, for incorporation in bed nets (or other vector control strategies like indoor residual sprays [64]). Developing these inhibitors for mosquito control may prove challenging given the close homology between mosquito and human lipases [21], although it is possible that allosteric inhibitors may show higher specificity against mosquito enzymes [65]. Potential impacts on nontarget organisms would need to be evaluated through toxicological testing in mammals and insect surveys during efficacy trials of the intervention. If necessary, effects could be mitigated through usage methods, dosages, or formulations that reduce exposure of other insects to these inhibitors, such as physical incorporation into bed nets, which generally only target night-biting species contacting the net. Importantly, as mentioned above, compounds that interfere with mosquito reproduction are already incorporated into bed nets utilized in malaria control strategies. Indeed, dual-ingredient nets including the insect growth regulator pyriproxyfen show promise in malaria endemic countries [59–61], as they act as secondary active agents in the event of mosquito resistance to the primary insecticide. This is particularly relevant at a time when insecticide-treated nets, after contributing to averting millions of malaria cases worldwide [66], are showing worrying signs of inefficacy [67], representing one of the principal causes behind the increase in malaria cases after years of decline [68,69]. Once developed, Anopheles-specific inhibitors of the lipolytic machinery could therefore contribute to feed the pipeline for the next generation of dual-ingredient bed net products.
Besides Lp-transported lipids, the yolk protein Vg conveys approximately 5% of total oocyte lipids [51], and so yolk deposition could potentially be affected in AgTL2 embryos and contribute to their abnormal development. Interestingly, embryonic lethality was reported when maternally deposited enzymes required for yolk usage in the embryo were silenced in Rhodnius prolixus embryos [70], suggesting that such phenotypes are due to an inability to mobilize yolk effectively. Our studies in An. gambiae have however shown that embryonic lethality occurs much earlier when Vg is silenced. Indeed, Vg depletion in mothers, which led to reduced yolk protein deposition into the oocyte, arrested embryonic development within the first few hours of oviposition [31], suggesting the late mortality observed in dsAgTL2 embryos is not due to defects in yolk deposition. The effects of AgTL2 silencing on Anopheles reproduction are in line with phenotypes observed in insects after manipulation of other components of the lipolytic machinery. In Drosophila, knockout of brummer (bmm), the orthologue of the mammalian adipose TAG lipase ATGL, in both parents leads to excessive fat accumulation and impaired lipid mobilization from fat bodies, and severely affects hatching of mutant embryos by unknown mechanisms [21], and similar results were obtained in recent studies in Bombyx mori [46] and R. prolixus [45]. Although these studies across insects did not determine the mechanisms triggered by inhibition of lipolysis, they highlight a degree of conservation in the importance of TAG mobilization in insect reproduction that could well be exploited to target populations of insect vectors of human and animal pathogens, as well as agricultural pests.
Our multi-omics analyses revealed that embryos from AgTL2-depleted mothers have considerable deficits in key metabolic processes. Carbohydrate metabolism and ß-oxidation are shut down in early developmental phases, followed by disruption of the TCA cycle. While these metabolic defects are not apparent in the histological samples and are not reflected by changes in transcriptomic profiles at the early RNAseq time points, it is likely that their accumulation leads eventually to the mis-regulation of several genes observed at the last time point analyzed (38 h post oviposition), a time when control embryos are close to hatching. Up-regulated genes are enriched in stress response and protein folding genes, revealing embryos that are under severe stress. This is paralleled with a strong down-regulation in key sensory signaling receptors and neuropeptides, eclosion hormone [54,55], and cuticular proteins that form the second cuticle needed before hatching [58]. LDs are known to be multifunctional storage compartments in the developing embryo [71], providing a source of not only TAGs and phospholipids but other lipids classes as well, such as the waxy saturated hydrocarbons required for cuticle synthesis [72]. Along with a significant decrease in expression of 26 annotated cuticle-related genes (S3 Data), a lack of LD breakdown in dsAgTL2 embryos may sequester these structural lipids away from cuticle synthesis enzymes, preventing the production of adequate cuticle for hatching [58]. Lipids stored within LDs may also serve as the source of several key signaling molecules, such as eicosanoids, hormones derived from polyunsaturated FFAs with roles in egg and eggshell formation [73]. Limited eicosanoid signaling may also explain embryonic lethality as eggshell defects and low hatching rates were recently reported in R. prolixus after bmm silencing [45], and eggshell defects have been observed after inhibition of eicosanoid signaling in Aedes mosquitoes [74].
We also observed a highly significant GO term enrichment for neural functions in down-regulated genes at 38 h post oviposition (top 4 terms), including neuropeptide and GPCR signaling pathways. CCHamide1-related proteins, induced in control but not dsAgTL2 embryos, are arthropod neuropeptides secreted by enteroendocrine cells in the gut upon sensing of different dietary conditions in the intestinal lumen. They regulate activity of dopaminergic neurons in the brain to control gut motility and other physiological behaviors [75–77], although interestingly they have not previously been implicated in embryogenesis and hatching. Lipids and their derivatives are required for neuronal signaling through their roles in synapse function and axonal wrapping [78,79]. Indeed, we observed increased levels of LPE and decreased ceramide phosphoethanolamine (CerPE) in dsAgTL2 embryos at 38 h post oviposition, changes that resemble Drosophila glial mutants unable to ensheathe axons [80], again suggestive of neuronal defects. In mammals, deficiencies in lipid metabolism and lipid signaling have been associated with developmental and cognitive problems as well as with neurodegenerative diseases. Imbalances in fatty acids during gestation causes altered brain development in the fetus, while severe fatty acid deprivation can impair cognitive functions (reviewed in [78]). Although we could not pinpoint the precise mechanisms behind embryonic lethality, our results show how reduced energy stored caused by impaired lipid remobilization and low TAG levels leads to the collapse of normal metabolic functions, eventually causing severe stress and cellular death.
Finally, our data confirm that P. falciparum parasites are affected by the mosquito lipolytic machinery. Depletion of AgTL2 (previous studies [48]) and AgLSD1 (this study) leads to accelerated oocyst growth rates and faster sporozoite invasion of salivary glands. Interestingly, however, such acceleration is not recapitulated by orlistat treatment, an observation that implicates fat body (rather than midgut) lipids in parasite growth. Indeed, while both AgTL2 depletion and orlistat treatment induce glyceride accumulation in the midgut (this study and [48]), differences in glyceride levels in the fat body are only detected after lipase silencing. This difference suggests that orlistat, while closely phenocopying the effects of AgTL2 depletion in terms of oogenesis and embryonic development, may inhibit different or additional aspects of the female lipolytic machinery. Regardless of the mechanism, when combined with the observed effects on fertility, these data support the idea that if Anopheles-specific lipase inhibitors were to be generated, they may represent promising tools for field interventions aimed at reducing the size of field populations of malaria vectors, aiding efforts to control this devastating disease.
Materials and methods
Rearing of An. gambiae mosquitoes
Anopheles gambiae (the standard laboratory wild-type strain G3) were reared in cages at 27°C, 80% humidity, and on a 12 h light-dark cycle. Adults were fed 10% glucose ad libitum and weekly on purchased human blood (Research Blood Components) for maintenance.
Culturing of P. falciparum parasites
Plasmodium falciparum (NF54 strain) parasites were maintained as asexual stages in human erythrocytes and gametocytes were induced as described in [48].
Embryo collection
Glass bowls containing circular Whatman filter papers submerged in deionized water were provided for 2 h to mated control and AgTL2-knockdown females at 4 days (d) PBM. Embryos were aged outside cages and collected for an assay of interest at 4 to 6 h, 22 to 24 h, and 36 to 38 h after oviposition.
Measurement of triglyceride (TAG) levels
Non-blood fed and blood fed control and AgTL2-silenced adult female An. gambiae mosquitoes were cooled on ice and decapitated in 1× phosphate-buffered saline (PBS). Ovaries, midgut, and fat body were dissected and transferred into separate tubes containing two 2 mm beads, diluted NP-40 substitute assay buffer (Triglyceride Assay Colorimetric kit, Cayman Chemicals), and 1× cocktail of protease inhibitor (cOmplete Protease Inhibitor Cocktail, EDTA-Free, Roche). For orlistat experiments, tissues were dissected from control or orlistat-treated mated wild-type blood fed mosquitoes at 48 h PBM. Each tube contained tissues pooled from 3 animals. For TAG determination in embryos, 200 embryos laid by control and AgTL2-silenced females were collected at 6, 24, and 30 h after oviposition. TAG quantification in tissues and embryos was performed according to the manufacturer’s instructions in the Triglyceride Assay Colorimetric kit (Cayman Chemicals). The assay quantifies glycerol content derived after release of fatty acids by lipoprotein lipase hydrolytic activity. Protein content in each pool was determined by Bradford assay and used for normalization of TAG levels.
Embryo fixation, immunofluorescence assays, and transmission electron microscopy (TEM)
Embryos were fixed as described by [81]. Briefly, embryos were rinsed with distilled water onto a filter paper at the indicated time point and bleached for 75 s with 25% bleach (2% sodium hypochlorite). After rinsing 6 to 8 times with distilled water, they were transferred into scintillation vials with 1:1 solution of 9% paraformaldehyde (PFA) and heptane and rotated on a roller shaker for 25 min at 70 rpm at room temperature. The water phase was replaced once and then washed with distilled water for 30 min on a shaker. Following incubation with boiled and iced-cold distilled water for 30 s and 15 s, respectively, the water phase was removed, and heptane was replaced. To remove the endochorion, an equal volume of methanol was added, and the vials were strongly swirled once and allowed to stand for 15 min at room temperature. The heptane-methanol mixture was removed, and embryos were washed 3 to 4 times with methanol and stored in fresh methanol at −20°C until the day of removal of eggshell. To remove eggshells, embryos were transferred onto a double-sided tape and peeled according to [82]. For immunofluorescence assays, embryos were permeabilized and blocked with a 2% bovine serum albumin (BSA) in 0.1% Tween-20 in 1× PBS solution for 30 min at room temperature.
For TEM, peeled embryos were collected in 200 μl of fixative (2.5% PFA, 5% glutaraldehyde, 0.06% picric acid in 0.2 M cacodylate buffer) and submitted to the Harvard Medical School Electron Microscopy Core. There, they were washed once with 0.1 M cacodylate buffer, twice with water, and then postfixed with 1% osmium tetroxide and 1.5% potassium ferrocyanide in water for 1 h at room temperature. Samples were then washed twice in water followed by once in 50 mM maleate buffer (pH 5.15) (MB). Next, the samples were incubated for 1 h in 1% uranyl acetate in MB, followed by 1 wash in MB, and 2 washes in water. Samples were then dehydrated via an ethanol gradient: increasing concentrations of ethanol were used with samples placed successively in 50%, 70%, 90%, 100%, 100% ethanol for 10 min each. After dehydration, samples were placed in propylene oxide for 1 h and then infiltrated overnight in a 1:1 mixture of propylene oxide and TAAB 812 Resin (www.taab.co.uk, #T022). The following day, samples were embedded in TAAB 812 Resin and polymerized at 60°C for 48 h. Ultrathin sections (roughly 80 nm) were cut on a Reichert Ultracut-S microtome, sections were picked up onto copper grids, stained with lead citrate, and imaged in a JEOL 1200EX transmission electron microscope equipped with an AMT 2K CCD camera.
Staining of lipids with LD540 or Nile Red
Ovaries, midguts, and fat body tissues were dissected from control dsGFP, TL2-silenced, control DMSO-, or orlistat-injected females at 48 h PBM and fixed in 4% PFA for 45 min at room temperature. Tissues were washed twice with 1× PBS for 15 min and incubated in permeabilization and blocking solution (1% BSA, 0.1% Tween 20, 1× PBS) for 30 min with shaking at room temperature. Samples were stained with LD540 in 1× PBS (0.5 μg/ml) according to [83]. Nuclei were stained with 4′,6-diamidino-2-phenylindole (DAPI) in 1× PBS (1 μg/ml) for 5 min and mounted on microscope slides with Vectorshield (Vector laboratories). Peeled embryos were fixed, permeabilized, and blocked as above but stained with Nile Red (10 μg/ml) for 30 min with shaking at room temperature. Images were captured on an Inverted Zeiss microscope Observer Z1.
Lipidomic and metabolomic analyses
All solvents were HPLC-MS grade from Sigma-Aldrich.
For lipidomic and metabolomic analyses of An. gambiae embryos, 200 embryos oviposited by control and AgTL2-knockdown blood fed mosquitoes in 3 biological replicates were collected at 4 to 6 h, 22 to 24 h, and 36 to 38 h after laying with 1 ml of methanol into a vial containing five 2 mm steel beads. Embryos were homogenized in a bead-beating homogenizer (Mini-Beadbeater 96, BioSpec Products) for 5 min at 2,400 rpm in a cold block and transferred to 8 ml glass vials (Duran Wheaton Kimble Life Sciences). After rinsing the bead-beater tubes with an additional 1 ml of methanol, 4 ml of cold chloroform was added to each glass vial and vortexed for 1 min, and 2 ml of ultrapure water was added to each vial and vortexed for another 1 min. Glass vials were centrifuged at 3,000 g for 10 min at 4°C, and 3 ml of the upper aqueous phase and 3 ml of the lower chloroform phase were transferred into separate glass vials for metabolomic and lipidomic analyses, respectively.
For metabolomics, samples were dried under nitrogen flow and resuspended in 50 μl of acetonitrile 25% in water. A pool was created by combining 10 μl of each sample. Samples were analyzed by LC-MS on a Vanquish LC coupled to an ID-X MS (Thermo Fisher Scientific); 5 μl of sample was injected on a ZIC-pHILIC PEEK-coated column (150 mm × 2.1 mm, 5 μm particles, maintained at 40°C, Sigma-Aldrich). Buffer A was 20 mM ammonium carbonate, 0.1% ammonium hydroxide in water and Buffer B was acetonitrile 97% in water. The LC program was as follow: starting at 93% B, to 40% B in 19 min, then to 0% B in 9 min, maintained at 0% B for 5 min, then back to 93% B in 3 min and re-equilibrated at 93% B for 9 min. The flow rate was maintained at 0.15 ml min-1, except for the first 30 s where the flow rate was uniformly ramped from 0.05 to 0.15 ml min-1. Data was acquired on the ID-X in switching polarities at 120,000 resolution, with an automatic gain control target of 1e5, and an m/z range of 65 to 1,000. MS1 data was acquired in switching polarities for all samples. MS2 and MS3 data was acquired on the pool sample using the AquirX DeepScan function, with 5 reinjections, separately in positive and negative ion mode. Data was analyzed in Compound Discoverer 3.2 (Thermo Fisher Scientific). Identification was based on MS2/MS3 matching with a local mzvault library and corresponding retention time built with pure standards (level 1) or on mzcloud match (level 2). Each match was manually inspected.
For lipidomics, samples were dried under nitrogen flow and resuspended in 60 μl of chloroform. Each sample was split into 2 equal aliquots, one for each polarity analysis. LC-MS analyses were modified from [84] and were performed on an Orbitrap Exactive plus MS (Thermo Fisher Scientific) in line with an Ultimate 3000 LC (Thermo Fisher Scientific). Each sample was analyzed in positive and negative modes, in top 5 automatic data-dependent MS/MS mode. Column hardware consisted of a Biobond C4 column (4.6 × 50 mm, 5 μm, Dikma Technologies). Flow rate was set to 100 μl min−1 for 5 min with 0% mobile phase B (MPB), then switched to 400 μl min−1 for 50 min, with a linear gradient of MPB from 20% to 100%. The column was then washed at 500 μl min−1 for 8 min at 100% MPB before being re-equilibrated for 7 min at 0% MPB and 500 μl min−1. For positive mode runs, buffers consisted for mobile phase A (MPA) of 5 mM ammonium formate, 0.1% formic acid, and 5% methanol in water, and for MPB of 5 mM ammonium formate, 0.1% formic acid, 5% water, 35% methanol in isopropanol. For negative runs, buffers consisted for MPA of 0.03% ammonium hydroxide, 5% methanol in water, and for MPB of 0.03% ammonium hydroxide, 5% water, 35% methanol in isopropanol. Lipids were identified and their signal integrated using the Lipidsearch software (version 4.2.27, Mitsui Knowledge Industry, University of Tokyo). Integrations and peak quality were curated manually before exporting and analyzing the data in Microsoft Excel. Heatmaps and principal component analysis plot were generated using MetaboAnalyst 5.0 [85].
Egg development and fertility assays
To determine the number of eggs laid and fertility after dsRNA injections or orlistat treatment, treated females were put into individual oviposition cups containing distilled water and Whatman filter paper at 48 to 72 h PBM. One week after oviposition, the number of eggs laid and the fertility was determined using a Nikon SMZ1000 dissecting microscope. Specifically, fertility was assessed by the degree of melanization of laid eggs (black for fertile and brown for infertile), followed by needle pricking of unhatched eggs to assess the presence (or absence) of fully developed larvae.
Knockdown of gene expression using dsRNA
Creation of templates for dsRNA production
cDNA from female whole bodies without heads was used for RT-qPCR amplification of fragments of interest. AgTL2 (AGAP000211, 406 bp) and AgLSD1 (AGAP002890, 485 bp) fragments were amplified using the following primers with the T7 promoter sequence added to their 5′ ends: AgTL2-Fwd: 5′-taatacgactcactatagggCTTTAGTCGGGTGGAGGACA-3′, AgTL2-Rev: 5′-taatacgactcactatagggGATGATCGGTGTTCGGCTT-3′; AgLSD1-Fwd: 5′-taatacgactcactatagggCGCAGATGGAGTCGCTGT-3′, AgLSD1-Rev: 5′-taatacgactcactatagggGCACGAGTCTGGCTGATCTT-3′. GFP fragments were amplified from plasmid pCR2.1-eGFP plasmid as in [48]. The PCR products were confirmed by gel electrophoresis and dsRNA was transcribed using the Megascript T7 transcription kit (Thermo Fisher Scientific) previously described [86,87].
dsRNA injections.
For single knockdown experiments, 690 ng of dsRNA (dsGFP, dsAgTL2, dsAgLSD1) was injected (Nanoject III, Drummond) into 1-day-old females at a concentration of 10 ng/nl. Gene knockdown levels were determined in at least 3 biological replicates by quantitative real-time PCR (RT-qPCR).
RNA extraction, cDNA synthesis, and RT-qPCR
For RNA extraction, 7 to 10 female mosquitoes were decapitated in 1× PBS and transferred to Tri-reagent (Thermo Fisher Scientific). Samples were stored at −80°C until the day of extraction. RNA was extracted according to the manufacturer’s instructions, with minor modifications. RNA pellet was washed 3 times with 75% ethanol in RNase-free water and air-dried for 7 to 8 min. Pellets were resuspended in heated water for 10 min and DNase-treated using Turbo DNase according to manufacturer’s instructions (Thermo Fisher Scientific). Purified RNA was quantified using a NanoDrop Spectrophotometer 2000c (Thermo Fisher Scientific) and cDNA synthesis was performed as in [87] using 2 μg of RNA in 100 μl cDNA reactions.
For quantification of target genes, RT-qPCR was performed in duplicate either in 7.5 or 15 μl reactions containing 1× PowerUp SYBR Green Master mix for RT-qPCR (Thermo Fisher Scientific), primers, and 2.5 or 5 μl of 10-fold diluted cDNA. Reactions were run on a QuantStudio 6 Pro real time-PCR thermocycler system (Thermo Fisher Scientific). Primers used for RT-qPCR had been previously used or designed spanning exon-exon junctions using Primer-BLAST [88]. Relative quantities of each gene in 3 or more biological replicates were determined using a standard curve with normalization against the ribosomal gene Rpl19 (AGAP004422). With the exception of RT-qPCR primers for AgTL2 that have previously been published [48], the following additional primer sets were used: AgLSD1: Fwd: 5′-CATTCAGCATACCTACGAC-3′, Rev: 5′-CACTGGTGGTGATCTTTT-3′; AgFAS: Fwd: 5′-GAGCTTGCGTAGAGATAGAT-3′, Rev: 5′-GCTAAGATTCTGTTGGAATG-3′; AGAP012959: Fwd: 5′-CTACTCTGTTCTTCCTGGTC-3′, Rev: 5′-TACCCTGTACTTCGGACTT-3′; AGAP013072: Fwd: 5′-TGAAGGCTACTCTGTTCTTC-3′, Rev: 5′-CCTGTACTTCGGACTTTC-3′; Hsp70 (AGAP004581): Fwd: 5′-CGTACTGTTTCAGTCTCAAG-3′, Rev: 5′-ATCGTACTCTTCCTTCTCTG-3′.
RNA sequencing and differential gene expression analysis
For RNAseq, RNA was extracted from embryos collected at different time points after oviposition from control and AgTL2-depleted blood fed An. gambiae females in 4 biological replicates as above, and then purified using PureLink RNA purification mini kit according to manufacturer’s instructions (Invitrogen). Libraries were prepared using a SciClone G3 NGSx workstation (Perkin Elmer) using the Kapa mRNA HyperPrep kit (Roche Applied Science). Polyadenylated mRNAs were captured using oligo-dT-conjugated magnetic beads (Kapa mRNA HyperPrep kit, Roche Sequencing) from 300 ng of total RNA on a Perkin Elmer SciClone G3 NGSx automated workstation. Poly-adenylated mRNA samples were immediately fragmented to 300 to 400 bp using heat and magnesium. First-strand synthesis was completed using random priming followed by second-strand synthesis and A-tailing. A dUTP was incorporated into the second strand to allow strand-specific sequencing of the library. Libraries were enriched and indexed using 11 cycles of amplification (Kapa mRNA HyperPrep kit, Roche Sequencing) with PCR primers, which included dual 8 bp index sequences to allow for multiplexing (IDT for Illumina unique dual 8 bp indexes). Excess PCR reagents were removed through magnetic bead-based cleanup using Kapa Pure magnetic beads on a SciClone G3 NGSx workstation (Perkin Elmer). The resulting libraries were assessed using a 4200 TapeStation (Agilent Technologies) and quantified by qPCR (Roche Sequencing). Libraries were pooled and sequenced on one lane of an Illumina NovaSeq S4 flow cell using paired-end, 100 bp reads.
Sequencing reads were aligned to the An. gambiae genome (PEST strain, version 4.14, VectorBase release 62) using HISAT2 (version 2.2.1) with the default parameters. Reads with mapping quality scores
Orlistat treatment
Approximately 5.0 to 6.7 mg orlistat (Sigma-Aldrich) was dissolved in 50 μl pure dimethyl sulfoxide (DMSO) to give a 200 to 270 mM stock solution, which was then diluted as required with 1× PBS to a 10 mM working solution, resulting in 3.7% to 5.0% DMSO in 1× PBS. The 5 mM and 2.5 mM orlistat solutions were prepared by 2-fold serial dilution with 1× PBS. Mated wild-type females were cooled on ice and pre-injected with control or diluted amounts of orlistat (1.0, 0.5, 0.25 μg/mosquito). The amount of orlistat injected was achieved by injecting 207 nl of either the 10, 5, or 2.5 mM solutions, respectively. The control group was injected with the same DMSO concentration as the 10 mM orlistat group (i.e., 3.7% to 5.0% DMSO in 1× PBS). After 3 to 4 h of recovery at 27°C, some mosquitoes were blood fed and ovaries, midgut, and fat body tissues were dissected from non-blood fed and blood fed females at 24 and 48 h after feeding for TAG analysis. For neutral lipid staining, tissues dissected at 48 h PBM were fixed and stained with the neutral lipid dye LD540. At 3 d PBM, the remaining blood fed females were put into individual oviposition cups to assess the number of eggs developed and fertility. For tarsal delivery, mated wild-type females were allowed to rest on their tarsi for 1 h on glass plates previously coated with orlistat at 1 or 2 mmol/m2 or control plates treated with the acetone vehicle alone. After 30 min of recovery, mosquitoes were blood fed and egg development, and fertility were determined after 7 d of oviposition.
P. falciparum infections of An. gambiae mosquitoes
Treated mosquitoes aged 4 to 6 d were housed in small cages and provided about 400 μl of P. falciparum gametocyte cultures for 1 h through heated membrane feeders. After blood feeding, poorly fed mosquitoes were removed and cages were maintained for the duration of the experiment in a glove box: 7 to 8 d post infectious feed for oocyst quantification or 13 d for sporozoite quantification. While in the glove box, mosquitoes were provided 10% glucose solution ad libitum.
Statistical analysis
Data were analyzed using GraphPad Prism 10 and JMP 16 Pro statistical software. GraphPad Prism 10 was used to calculate Fisher’s exact test, Student’s t tests, and 1- or 2-way analysis of variance (ANOVA) with Dunnett’s or Šídák’s multiple comparisons correction. Prior to ANOVA, data were first assessed for normal distribution and were transformed if found not be normally distributed. JMP Pro 16 was used to generate Least Square Means models incorporating treatment, time, and their interactions as fixed effects and replicate as a random effect. A post hoc t test between groups was done to assess significance and p-values corrected using a false discovery rate (FDR) of 0.05.
Supporting information
S1 Fig. Gene silencing efficacy and lipid staining in AgTL2-depleted and control females after blood feeding.
(A) AgTL2 is significantly induced in the fat body but not the ovaries or midgut of females at 2 time points after blood feeding (hPBM = hours post blood meal) compared to non-blood fed (NBF) controls. (Ordinary 2-way ANOVA, Šídák’s multiple comparisons correction) (B) AgTL2 levels after blood feeding are significantly reduced upon RNAi injections in the fat body but not the ovaries or midgut (Ordinary 2-way ANOVA, Šídák’s multiple comparisons correction). (C) AgLSD1 levels are reduced upon RNAi injections at 3 d post injection in whole body minus head samples (Rest of body) (unpaired t test). Each dot in A–C represents RT-qPCR analysis from a pool of 10 mosquitoes per replicate (3–5 biological replicates). (D) Neutral lipid staining with LD540 (red) shows intense staining in midgut and fat body tissues and weak staining in ovaries from AgTL2-depleted females at 48 hPBM (scale bar = 10 μm for midguts, 50 μm for ovaries and fat body images). Blue = DAPI (DNA). Numerical data supporting this figure is available in the Harvard Dataverse online repository at https://doi.org/10.7910/DVN/ULTW1K.
https://doi.org/10.1371/journal.pbio.3002960.s001
(TIF)
S2 Fig. dsAgTL2 embryos have altered lipid profiles.
(A) Heatmap of major lipids analyzed by LC-MS reveals dysregulated levels in dsAgTL2 embryos compared to controls throughout development. (B) Lysophosphatidylethanolamine (LPE) levels are significantly increased in dsAgTL2 embryos (Ordinary 2-way ANOVA, Šídák’s multiple comparisons correction, * = p t test statistic, range: blue to red = significant decrease to increase. Three biological replicates, represented by triplicate columns at each timepoint. See S1 Table and S1 Data. Numerical data supporting this figure is available in the Harvard Dataverse online repository at https://doi.org/10.7910/DVN/ULTW1K.
https://doi.org/10.1371/journal.pbio.3002960.s002
(TIF)
S3 Fig. Key metabolites are significantly altered in embryos from AgTL2-depleted mothers.
LC-MS analysis reveals dysregulated levels of key metabolites in embryos from AgTL2-deficient females at one or multiple time points after oviposition. Statistical significance was assessed by Least Square Means models testing effect of treatment, time points, and replicate. Adjusted p-values were calculated using an FDR = 0.05. Not all significant comparisons are shown for clarity. See S2 Data. (ns = not significant, * = p p p p https://doi.org/10.7910/DVN/ULTW1K.
https://doi.org/10.1371/journal.pbio.3002960.s003
(TIF)
S4 Fig. Gene ontology of dysregulated genes and validation by RT-qPCR.
(A) At 38 h post oviposition, genes involved in neuropeptide signaling pathway, sensory perception, and nervous system processes are most significantly represented among down-regulated genes in dsAgTL2 embryos compared to controls, while genes involved in protein translation, protein folding, and metabolic processes are most significantly represented in up-regulated genes (numbers in brackets indicate the number of genes in each data set associated with a biological process). (B) RT-qPCR analyses confirm significantly down-regulation of 2 Cchamide-1-related genes (neuropeptide signaling pathway) at 38 h post oviposition in dsAgTL2 embryos compared to controls, while Heat shock protein 70 (Hsp70) (protein folding) and fatty acid synthase (FAS) (metabolic process) are significantly up-regulated. Kruskal–Wallis test with Dunn’s multiple comparisons correction. See S3 Data. n = number of pools of at least 50 embryos per female from 2 biological replicates. Numerical data supporting this figure is available in the Harvard Dataverse online repository at https://doi.org/10.7910/DVN/ULTW1K.
https://doi.org/10.1371/journal.pbio.3002960.s004
(TIF)
S5 Fig. Overall representation of metabolic pathways affected in embryos from AgTL2-silenced mothers, based on lipidomics, metabolomics, and transcriptomics analyses.
dsAgTL2 embryos have significantly reduced levels of main insect sugars and the glycolytic intermediate fructose-bisphosphate, reduced levels of early TCA and ß-oxidation intermediates citrate, aconitate and Acetyl L-carnitine, altered levels of amino acids, and increased levels of nucleotide degradation products (xanthine, hypoxanthine, uric acid, and 3-ureidopropionic acid). Only data from the 38 h time point are shown (but Least Square Means models testing effect of treatment, time point, and replicate were built on the entire time course). Metabolites highlighted in blue are significantly reduced and those highlighted in red are significantly increased. Adjusted p-values were calculated using an FDR = 0.05. Not all significant comparisons are shown, for clarity. See S3 Fig and S2 Data. (* = p p p p https://doi.org/10.7910/DVN/ULTW1K.
https://doi.org/10.1371/journal.pbio.3002960.s005
(TIF)
S6 Fig. AgLSD1 accelerates parasite growth and orlistat treatment does not affect lipid levels in fat body.
(A–C) AgLSD1 knockdown has (A) no effect on P. falciparum infection intensity (Mann–Whitney) and prevalence (P) (Pie charts, Fisher’s exact, ns; p > 0.05) compared to dsGFP controls but leads to (B) increased oocyst size on 7 d PBM and (C) higher sporozoite numbers in salivary glands at 12–13 d PBM (Unpaired t test). (D) Orlistat treatment leads to increased glyceride levels in midguts and decreased levels in ovaries compared to control-injected mosquitoes, but no changes are observed in the fat body. All data at 48 h PBM (Ordinary 1-way ANOVA, Dunnett’s multiple comparisons correction). Three–four biological replicates are represented in A–C, and 4 in D. n = number of individual mosquitoes analyzed (except for glyceride assay, where n represents the number of pools of 3 insects). Numerical data supporting this figure is available in the Harvard Dataverse online repository at https://doi.org/10.7910/DVN/ULTW1K.
https://doi.org/10.1371/journal.pbio.3002960.s006
(TIF)
S1 Data. Major lipid levels in embryos from control and AgTL2-depleted females at different time points after oviposition.
Lipid levels (as peak area) determined after LC-MS of control and AgTL2-depleted embryos at 4–6 h, 22–24 h, and 36–38 h after oviposition. Values are given for each class and subclass of lipid by fatty acid composition. Total values are summed under each major lipid class in bold. The raw lipidomics data are available from the EMBL-EBI Metabolights online repository using the link https://www.ebi.ac.uk/metabolights/MTBLS9881.
https://doi.org/10.1371/journal.pbio.3002960.s008
(XLSX)
S2 Data. Intensities of metabolites determined by LC-MS in embryos from control and AgTL2-depleted mothers at different time points after oviposition.
Metabolite levels (as peak area) determined after LC-MS of control and AgTL2-depleted embryos at 4–6 h, 22–24 h, and 36–38 h after oviposition. Values are given for each metabolite and confidence of metabolite identification (ID). The raw metabolomics data are available from the EMBL-EBI Metabolights online repository using the link https://www.ebi.ac.uk/metabolights/MTBLS9881.
https://doi.org/10.1371/journal.pbio.3002960.s009
(XLSX)
S3 Data. RNAseq data sets for embryos derived from control and AgTL2-depleted An. gambiae mosquitoes at different time points after oviposition.
Transcriptomics of control and AgTL2-depleted embryos at 4–6 h, 22–24 h, and 36–38 h after oviposition. Each time point has 3 associated spreadsheets: all genes, up-regulated genes and down-regulated genes. Values are given for mean ± SE log2 (fold change on TL2 depletion), test statistic, and unadjusted and adjusted p-values, with gene description and name if available. The raw RNAseq data are available from the GEO database using the link https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE263712.
https://doi.org/10.1371/journal.pbio.3002960.s010
(XLSX)
Acknowledgments
We are particularly grateful to Aaron Stanton, Emily Selland, and Elizabeth Nelson for mosquito rearing and all the members of the Catteruccia laboratory for helpful suggestions. We are also grateful to The Bauer Core Facility at Harvard University for RNA sequencing and the Harvard Medical School Electron Microscopy Facility for TEM.
The findings and conclusions within this publication are those of the authors and do not reflect the policies or positions of the HHMI or the NIH.
References
- 1.
Liang L, Rasmussen MH, Piening B, Shen X, Chen S, Rost H, et al. Metabolic Dynamics and Prediction of Gestational Age and Time to Delivery in Pregnant Women. Cell. 2020;181(7):1680–92 e15. pmid:32589958 - 2.
Wang Q, Wurtz P, Auro K, Makinen VP, Kangas AJ, Soininen P, et al. Metabolic profiling of pregnancy: cross-sectional and longitudinal evidence. BMC Med. 2016;14(1):205. pmid:27955712 - 3.
Dunning KR, Russell DL, Robker RL. Lipids and oocyte developmental competence: the role of fatty acids and beta-oxidation. Reproduction. 2014;148(1):R15–R27. - 4.
Sturmey RG, Reis A, Leese HJ, McEvoy TG. Role of fatty acids in energy provision during oocyte maturation and early embryo development. Reprod Domest Anim. 2009;44(Suppl 3):50–58. pmid:19660080 - 5.
Grabner GF, Xie H, Schweiger M, Zechner R. Lipolysis: cellular mechanisms for lipid mobilization from fat stores. Nat Metab. 2021;3(11):1445–1465. pmid:34799702 - 6.
Walther TC, Farese RV Jr. The life of lipid droplets. Biochim Biophys Acta. 2009;1791(6):459–466. pmid:19041421 - 7.
Egan JJ, Greenberg AS, Chang MK, Wek SA, Moos MC Jr, Londos C. Mechanism of hormone-stimulated lipolysis in adipocytes: translocation of hormone-sensitive lipase to the lipid storage droplet. Proc Natl Acad Sci U S A. 1992;89(18):8537–8541. pmid:1528859 - 8.
Zimmermann R, Strauss JG, Haemmerle G, Schoiswohl G, Birner-Gruenberger R, Riederer M, et al. Fat mobilization in adipose tissue is promoted by adipose triglyceride lipase. Science. 2004;306(5700):1383–1386. pmid:15550674 - 9.
Sztalryd C, Xu G, Dorward H, Tansey JT, Contreras JA, Kimmel AR, et al. Perilipin A is essential for the translocation of hormone-sensitive lipase during lipolytic activation. J Cell Biol. 2003;161 (6):1093–1103. pmid:12810697 - 10.
Zhang HH, Souza SC, Muliro KV, Kraemer FB, Obin MS, Greenberg AS. Lipase-selective functional domains of perilipin A differentially regulate constitutive and protein kinase A-stimulated lipolysis. J Biol Chem. 2003;278(51):51535–51542. pmid:14527948 - 11.
Yamaguchi T, Omatsu N, Matsushita S, Osumi T. CGI-58 interacts with perilipin and is localized to lipid droplets. Possible involvement of CGI-58 mislocalization in Chanarin-Dorfman syndrome. J Biol Chem. 2004;279(29):30490–30497. pmid:15136565 - 12.
Yamaguchi T, Omatsu N, Morimoto E, Nakashima H, Ueno K, Tanaka T, et al. CGI-58 facilitates lipolysis on lipid droplets but is not involved in the vesiculation of lipid droplets caused by hormonal stimulation. J Lipid Res. 2007;48(5):1078–1089. pmid:17308334 - 13.
Lass A, Zimmermann R, Haemmerle G, Riederer M, Schoiswohl G, Schweiger M, et al. Adipose triglyceride lipase-mediated lipolysis of cellular fat stores is activated by CGI-58 and defective in Chanarin-Dorfman Syndrome. Cell Metab. 2006;3(5):309–319. - 14.
Cho CH, Patel S, Rajbhandari P. Adipose tissue lipid metabolism: lipolysis. Curr Opin Genet Dev. 2023;83:102114. pmid:37738733 - 15.
Kim SK, Rulifson EJ. Conserved mechanisms of glucose sensing and regulation by Drosophila corpora cardiaca cells. Nature. 2004;431(7006):316–320. pmid:15372035 - 16.
Mayer RJ, Candy DJ. Changes in energy reserves during flight of the desert locust, Schistocerca gregaria. Comp Biochem Physiol. 1969;31(3):409–418. pmid:5351972 - 17.
Arrese EL, Flowers MT, Gazard JL, Wells MA. Calcium and cAMP are second messengers in the adipokinetic hormone-induced lipolysis of triacylglycerols in Manduca sexta fat body. J Lipid Res. 1999;40(3):556–564. pmid:10064744 - 18.
Arrese EL, Mirza S, Rivera L, Howard AD, Chetty PS, Soulages JL. Expression of lipid storage droplet protein-1 may define the role of AKH as a lipid mobilizing hormone in Manduca sexta. Insect Biochem Mol Biol. 2008;38(11):993–1000. pmid:18793726 - 19.
Ziegler R, Eckart K, Law JH. Adipokinetic hormone controls lipid metabolism in adults and carbohydrate metabolism in larvae of Manduca sexta. Peptides. 1990;11(5):1037–1040. pmid:2284195 - 20.
Arrese EL, Soulages JL. Insect fat body: energy, metabolism, and regulation. Annu Rev Entomol. 2010;55:207–225. pmid:19725772 - 21.
Gronke S, Mildner A, Fellert S, Tennagels N, Petry S, Muller G, et al. Brummer lipase is an evolutionary conserved fat storage regulator in Drosophila. Cell Metab. 2005;1(5):323–330. - 22.
Dos Santos LV, Silva E, Caiado MS, Rezende S, de Carvalho MG, Pontes EG. Differential expression of brummer and levels of TAG in different developmental stages Aedes aegypti (Diptera: Culicidae), including fasted adults. Arch Insect Biochem Physiol. 2024;115(1):e22084. pmid:38288494 - 23.
Geering K, Freyvogel TA. Lipase activity and stimulation mechanism of esterases in the midgut of female Aedes aegypti. J Insect Physiol. 1975;21(6):1251–1256. pmid:1133492 - 24.
Weintraub H, Tietz A. Triglyceride digestion and absorption in the locust, Locusta migratoria. Biochim Biophys Acta. 1973;306(1):31–41. pmid:4703573 - 25.
Kaufmann C, Briegel H. Flight performance of the malaria vectors Anopheles gambiae and Anopheles atroparvus. J Vector Ecol. 2004;29(1):140–153. pmid:15266751 - 26.
Ziegler R, Isoe J, Moore W, Riehle MA, Wells MA. The putative AKH receptor of the tobacco hornworm, Manduca sexta, and its expression. J Insect Sci. 2011;11:40. pmid:21529255 - 27.
Pinch M, Mitra S, Rodriguez SD, Li Y, Kandel Y, Dungan B, et al. Fat and Happy: Profiling Mosquito Fat Body Lipid Storage and Composition Post-blood Meal. Front Insect Sci. 2021;1:693168. pmid:38468893 - 28.
Briegel H, Waltert A, Kuhn AR. Reproductive physiology of Aedes (Aedimorphus) vexans (Diptera: Culicidae) in relation to flight potential. J Med Entomol. 2001;38(4):557–565. pmid:11476336 - 29.
Briegel H, Hefti M, DiMarco E. Lipid metabolism during sequential gonotrophic cycles in large and small female Aedes aegypti. J Insect Physiol. 2002;48(5):547–554. pmid:12770082 - 30.
Raikhel AS, Dhadialla TS. Accumulation of yolk proteins in insect oocytes. Annu Rev Entomol. 1992;37:217–251. pmid:1311540 - 31.
Stryapunina I, Itoe MA, Trinh Q, Vidoudez C, Du E, Mendoza L, et al. Precise coordination between nutrient transporters ensures fertility in the malaria mosquito Anopheles gambiae. PLoS Genet. 2024;20(1):e1011145. pmid:38285728 - 32.
Raikhel AS, Lea AO. Internalized proteins directed into accumulative compartments of mosquito oocytes by the specific ligand, vitellogenin. Tissue Cell. 1986;18(4):559–574. pmid:3750316 - 33.
Cho WL, Tsao SM, Hays AR, Walter R, Chen JS, Snigirevskaya ES, et al. Mosquito cathepsin B-like protease involved in embryonic degradation of vitellin is produced as a latent extraovarian precursor. J Biol Chem. 1999;274(19):13311–13321. pmid:10224092 - 34.
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–E2718. pmid:28292900 - 35.
Hossain MS, Liu Y, Zhou S, Li K, Tian L, Li S. 20-Hydroxyecdysone-induced transcriptional activity of FoxO upregulates brummer and acid lipase-1 and promotes lipolysis in Bombyx fat body. Insect Biochem Mol Biol. 2013;43(9):829–838. - 36.
Zhang J, Liu Y, Jiang K, Jia J. Hedgehog signaling promotes lipolysis in adipose tissue through directly regulating Bmm/ATGL lipase. Dev Biol. 2020;457(1):128–139. pmid:31550483 - 37.
Qian W, Guo M, Peng J, Zhao T, Li Z, Yang Y, et al. Decapentaplegic retards lipolysis during metamorphosis in Bombyx mori and Drosophila melanogaster. Insect Biochem Mol Biol. 2023;155:103928. pmid:36870515 - 38.
Molaei M, Vandehoef C, Karpac J. NF-kappaB Shapes Metabolic Adaptation by Attenuating Foxo-Mediated Lipolysis in Drosophila. Dev Cell. 2019;49(5):802–10 e6. - 39.
Song T, Qin W, Lai Z, Li H, Li D, Wang B, et al. Dietary cysteine drives body fat loss via FMRFamide signaling in Drosophila and mouse. Cell Res. 2023;33(6):434–447. pmid:37055592 - 40.
Dou X, Chen K, Brown MR, Strand MR. Multiple endocrine factors regulate nutrient mobilization and storage in Aedes aegypti during a gonadotrophic cycle. Insect Sci. 2023;30(2):425–442. pmid:36056560 - 41.
Grillo LA, Majerowicz D, Gondim KC. Lipid metabolism in Rhodnius prolixus (Hemiptera: Reduviidae): role of a midgut triacylglycerol-lipase. Insect Biochem Mol Biol. 2007;37(6):579–588. pmid:17517335 - 42.
Alves-Bezerra M, De Paula IF, Medina JM, Silva-Oliveira G, Medeiros JS, Gade G, et al. Adipokinetic hormone receptor gene identification and its role in triacylglycerol metabolism in the blood-sucking insect Rhodnius prolixus. Insect Biochem Mol Biol. 2016;69:51–60. pmid:26163435 - 43.
Attardo GM, Benoit JB, Michalkova V, Yang G, Roller L, Bohova J, et al. Analysis of lipolysis underlying lactation in the tsetse fly, Glossina morsitans. Insect Biochem Mol Biol. 2012;42(5):360–370. pmid:22509523 - 44.
Kaufmann C, Brown MR. Adipokinetic hormones in the African malaria mosquito, Anopheles gambiae: identification and expression of genes for two peptides and a putative receptor. Insect Biochem Mol Biol. 2006;36(6):466–481. pmid:16731343 - 45.
Aredes DS, Rios T, Carvalho-Kelly LF, Braz V, Araripe LO, Bruno RV, et al. Deficiency of Brummer lipase disturbs lipid mobilization and locomotion, and impairs reproduction due to defects in the eggshell ultrastructure in the insect vector Rhodnius prolixus. Biochim Biophys Acta Mol Cell Biol Lipids. 2024;1869(2):159442. pmid:38042331 - 46.
Yang F, Xu X, Hu B, Zhang Z, Chen K, Yu Y, et al. Lipid homeostasis is essential for oogenesis and embryogenesis in the silkworm, Bombyx mori. Cell Mol Life Sci. 2024;81(1):127. pmid:38472536 - 47.
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–1195. pmid:16005290 - 48.
Werling K, Shaw WR, Itoe MA, Westervelt KA, Marcenac P, Paton DG, et al. Steroid Hormone Function Controls Non-competitive Plasmodium Development in Anopheles. Cell. 2019;177(2):315–25 e14. pmid:30929905 - 49.
Shaw WR, Holmdahl IE, Itoe MA, Werling K, Marquette M, Paton DG, et al. Multiple blood feeding in mosquitoes shortens the Plasmodium falciparum incubation period and increases malaria transmission potential. PLoS Pathog. 2020;16(12):e1009131. pmid:33382824 - 50.
Costa G, Gildenhard M, Eldering M, Lindquist RL, Hauser AE, Sauerwein R, et al. Non-competitive resource exploitation within mosquito shapes within-host malaria infectivity and virulence. Nat Commun. 2018;9(1):3474. pmid:30150763 - 51.
Kawooya JK, Law JH. Role of lipophorin in lipid transport to the insect egg. J Biol Chem. 1988;263(18):8748–8753. pmid:3379043 - 52.
Marinotti O, Calvo E, Nguyen QK, Dissanayake S, Ribeiro JM, James AA. Genome-wide analysis of gene expression in adult Anopheles gambiae. Insect Mol Biol. 2006;15(1):1–12. pmid:16469063 - 53.
Nyasembe VO, Hamerly T, Lopez-Gutierrez B, Leyte-Vidal AM, Coatsworth H, Dinglasan RR. Adipokinetic hormone signaling in the malaria vector Anopheles gambiae facilitates Plasmodium falciparum sporogony. Commun Biol. 2023;6(1):171. pmid:36782045 - 54.
Kubrak OI, Lushchak OV, Zandawala M, Nassel DR. Systemic corazonin signalling modulates stress responses and metabolism in Drosophila. Open Biol. 2016;6(11). pmid:27810969 - 55.
Ewer J. Behavioral actions of neuropeptides in invertebrates: insights from Drosophila. Horm Behav. 2005;48(4):418–429. pmid:15996666 - 56.
Li S, Torre-Muruzabal T, Sogaard KC, Ren GR, Hauser F, Engelsen SM, et al. Expression patterns of the Drosophila neuropeptide CCHamide-2 and its receptor may suggest hormonal signaling from the gut to the brain. PLoS ONE. 2013;8(10):e76131. pmid:24098432 - 57.
Ida T, Takahashi T, Tominaga H, Sato T, Sano H, Kume K, et al. Isolation of the bioactive peptides CCHamide-1 and CCHamide-2 from Drosophila and their putative role in appetite regulation as ligands for G protein-coupled receptors. Front Endocrinol (Lausanne). 2012;3:177. pmid:23293632 - 58.
Moussian B. Recent advances in understanding mechanisms of insect cuticle differentiation. Insect Biochem Mol Biol. 2010;40(5):363–375. pmid:20347980 - 59.
Tiono AB, Ouedraogo A, Ouattara D, Bougouma EC, Coulibaly S, Diarra A, et al. Efficacy of Olyset Duo, a bednet containing pyriproxyfen and permethrin, versus a permethrin-only net against clinical malaria in an area with highly pyrethroid-resistant vectors in rural Burkina Faso: a cluster-randomised controlled trial. Lancet. 2018;392(10147):569–580. pmid:30104047 - 60.
Accrombessi M, Cook J, Dangbenon E, Yovogan B, Akpovi H, Sovi A, et al. Efficacy of pyriproxyfen-pyrethroid long-lasting insecticidal nets (LLINs) and chlorfenapyr-pyrethroid LLINs compared with pyrethroid-only LLINs for malaria control in Benin: a cluster-randomised, superiority trial. Lancet. 2023;401(10375):435–446. pmid:36706778 - 61.
Mosha JF, Kulkarni MA, Lukole E, Matowo NS, Pitt C, Messenger LA, et al. Effectiveness and cost-effectiveness against malaria of three types of dual-active-ingredient long-lasting insecticidal nets (LLINs) compared with pyrethroid-only LLINs in Tanzania: a four-arm, cluster-randomised trial. Lancet. 2022;399(10331):1227–1241. pmid:35339225 - 62.
Aiku AO, Yates A, Rowland M. Laboratory evaluation of pyriproxifen treated bednets on mosquito fertility and fecundity. A preliminary study. West Afr J Med. 2006;25(1):22–26. pmid:16722354 - 63.
Hadvary P, Sidler W, Meister W, Vetter W, Wolfer H. The lipase inhibitor tetrahydrolipstatin binds covalently to the putative active site serine of pancreatic lipase. J Biol Chem. 1991;266(4):2021–2027. pmid:1899234 - 64.
Sherrard-Smith E, Griffin JT, Winskill P, Corbel V, Pennetier C, Djenontin A, et al. Systematic review of indoor residual spray efficacy and effectiveness against Plasmodium falciparum in Africa. Nat Commun. 2018;9(1):4982. pmid:30478327 - 65.
Wagner JR, Lee CT, Durrant JD, Malmstrom RD, Feher VA, Amaro RE. Emerging Computational Methods for the Rational Discovery of Allosteric Drugs. Chem Rev. 2016;116(11):6370–6390. pmid:27074285 - 66.
Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The effect of malaria control on Plasmodium falciparum in Africa between 2000 and 2015. Nature. 2015;526(7572):207–211. pmid:26375008 - 67.
Churcher TS, Lissenden N, Griffin JT, Worrall E, Ranson H. The impact of pyrethroid resistance on the efficacy and effectiveness of bednets for malaria control in Africa. Elife. 2016;5:e16090. pmid:27547988 - 68.
Ranson H, Lissenden N. Insecticide Resistance in African Anopheles Mosquitoes: A Worsening Situation that Needs Urgent Action to Maintain Malaria Control. Trends Parasitol. 2016;32(3):187–196. pmid:26826784 - 69.
WHO. World Malaria Report. Geneva: World Health Organization; 2023. Report No.: 978-92-4-008617-3. - 70.
de Almeida E, Dittz U, Pereira J, Walter-Nuno AB, Paiva-Silva GO, Lacerda-Abreu MA, et al. Functional characterization of maternally accumulated hydrolases in the mature oocytes of the vector Rhodnius prolixus reveals a new protein phosphatase essential for the activation of the yolk mobilization and embryo development. Front Physiol. 2023;14:1142433. pmid:36923285 - 71.
Li T, Jin Y, Wu J, Ren Z. Beyond energy provider: multifunction of lipid droplets in embryonic development. Biol Res. 2023;56(1):38. pmid:37438836 - 72.
Fruttero LL, Leyria J, Canavoso LE. Lipids in Insect Oocytes: From the Storage Pathways to Their Multiple Functions. Results Probl Cell Differ. 2017;63:403–434. pmid:28779328 - 73.
Vrablik TL, Watts JL. Polyunsaturated fatty acid derived signaling in reproduction and development: insights from Caenorhabditis elegans and Drosophila melanogaster. Mol Reprod Dev. 2013;80(4):244–259. pmid:23440886 - 74.
Choi DY, Kim Y. Prostaglandin E(2) mediates chorion formation of the Asian tiger mosquito, Aedes albopictus, at late oogenesis. Insect Mol Biol. 2023;32(5):484–509. pmid:37158315 - 75.
Titos I, Juginovic A, Vaccaro A, Nambara K, Gorelik P, Mazor O, et al. A gut-secreted peptide suppresses arousability from sleep. Cell. 2023;186(7):1382–97 e21. pmid:36958331 - 76.
Psichas A, Reimann F, Gribble FM. Gut chemosensing mechanisms. J Clin Invest. 2015;125(3):908–917. pmid:25664852 - 77.
Nassel DR, Zandawala M. Recent advances in neuropeptide signaling in Drosophila, from genes to physiology and behavior. Prog Neurobiol. 2019;179:101607. pmid:30905728 - 78.
Falomir-Lockhart LJ, Cavazzutti GF, Gimenez E, Toscani AM. Fatty Acid Signaling Mechanisms in Neural Cells: Fatty Acid Receptors. Front Cell Neurosci. 2019;13:162. pmid:31105530 - 79.
Murate M, Tomishige N, Kobayashi T. Wrapping axons in mammals and Drosophila: Different lipids, same principle. Biochimie. 2020;178:39–48. pmid:32800899 - 80.
Kunduri G, Turner-Evans D, Konya Y, Izumi Y, Nagashima K, Lockett S, et al. Defective cortex glia plasma membrane structure underlies light-induced epilepsy in cpes mutants. Proc Natl Acad Sci U S A. 2018;115(38):E8919–E8928. pmid:30185559 - 81.
Goltsev Y, Hsiong W, Lanzaro G, Levine M. Different combinations of gap repressors for common stripes in Anopheles and Drosophila embryos. Dev Biol. 2004;275(2):435–446. pmid:15501229 - 82.
Juhn J, James AA. Hybridization in situ of salivary glands, ovaries, and embryos of vector mosquitoes. J Vis Exp. 2012;64. pmid:22781778 - 83.
Spandl J, White DJ, Peychl J, Thiele C. Live cell multicolor imaging of lipid droplets with a new dye, LD540. Traffic. 2009;10(11):1579–1584. pmid:19765264 - 84.
Miraldi ER, Sharfi H, Friedline RH, Johnson H, Zhang T, Lau KS, et al. Molecular network analysis of phosphotyrosine and lipid metabolism in hepatic PTP1b deletion mice. Integr Biol (Camb). 2013;5(7):940–963. pmid:23685806 - 85.
Pang Z, Zhou G, Ewald J, Chang L, Hacariz O, Basu N, et al. Using MetaboAnalyst 5.0 for LC-HRMS spectra processing, multi-omics integration and covariate adjustment of global metabolomics data. Nat Protoc. 2022;17(8):1735–1761. pmid:35715522 - 86.
Blandin S, Shiao SH, Moita LF, Janse CJ, Waters AP, Kafatos FC, et al. Complement-like protein TEP1 is a determinant of vectorial capacity in the malaria vector Anopheles gambiae. Cell. 2004;116(5):661–670. pmid:15006349 - 87.
Rogers DW, Baldini F, Battaglia F, Panico M, Dell A, Morris HR, et al. Transglutaminase-mediated semen coagulation controls sperm storage in the malaria mosquito. PLoS Biol. 2009;7(12):e1000272. pmid:20027206 - 88.
Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, Madden TL. Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics. 2012;13:134. pmid:22708584
ADVERTISEMENT:
Halo, para pengemar slots Pernah denger istilah “slot gaco”? Kalau belum, siap-siap jatuh cinta sama konsep ini. slot gacor adalah mesin slots yang sering memberi win. Ya, slot-slot ini bisa disebut sebagai andalannya buat membawa come back hasil. but, gimana sih
tekniknya nemuin slot gaco yang benar? Santai Bro and Sis, kita bahas tenang saja di tempat ini
Games terpercaya saat sekarang satu-satunya berada Indonesia yaitu yang memberikan ROI terbaik
Daftar hanya di :
Informasi mengenai KING SLOT, Segera Daftar Bersama king selot terbaik dan terpercaya no satu di Indonesia. Boleh mendaftar melalui sini king slot serta memberikan hasil kembali yang paling tinggi saat sekarang ini hanyalah KING SLOT atau Raja slot paling gacor, gilak dan gaco saat sekarang di Indonesia melalui program return tinggi di kingselot serta pg king slot
slot demo gacor
slot demo gacor permainan paling top dan garansi imbal balik hasil besar bersama kdwapp.com
akun demo slot gacor
akun demo slot gacor permainan paling top dan garansi imbal balik hasil besar bersama kdwapp.com
akun slot demo gacor
akun slot demo gacor permainan paling top dan garansi imbal balik hasil besar bersama kdwapp.com
akun demo slot pragmatic
akun demo slot pragmatic permainan paling top dan garansi imbal balik hasil besar bersama kdwapp.com
akun slot demo pragmatic
akun slot demo pragmatic permainan paling top dan garansi imbal balik hasil besar bersama kdwapp.com
akun slot demo
akun slot demo permainan paling top dan garansi imbal balik hasil besar bersama kdwapp.com
akun demo slot
akun demo slot permainan paling top dan garansi imbal balik hasil besar bersama kdwapp.com
slot demo gacor
slot demo gacor permainan paling top dan garansi imbal balik hasil besar bersama jebswagstore.com
akun demo slot gacor
akun demo slot gacor permainan paling top dan garansi imbal balik hasil besar bersama jebswagstore.com
akun slot demo gacor
akun slot demo gacor permainan paling top dan garansi imbal balik hasil besar bersama jebswagstore.com
akun demo slot pragmatic
akun demo slot pragmatic permainan paling top dan garansi imbal balik hasil besar bersama jebswagstore.com
akun slot demo pragmatic
akun slot demo pragmatic permainan paling top dan garansi imbal balik hasil besar bersama jebswagstore.com
akun slot demo
akun slot demo permainan paling top dan garansi imbal balik hasil besar bersama jebswagstore.com
akun demo slot
akun demo slot permainan paling top dan garansi imbal balik hasil besar bersama jebswagstore.com
slot demo gacor
slot demo gacor permainan paling top dan garansi imbal balik hasil besar bersama demoslotgacor.pro
akun demo slot gacor
akun demo slot gacor permainan paling top dan garansi imbal balik hasil besar bersama demoslotgacor.pro
akun slot demo gacor
akun slot demo gacor permainan paling top dan garansi imbal balik hasil besar bersama demoslotgacor.pro
akun demo slot pragmatic
akun demo slot pragmatic permainan paling top dan garansi imbal balik hasil besar bersama demoslotgacor.pro
akun slot demo pragmatic
akun slot demo pragmatic permainan paling top dan garansi imbal balik hasil besar bersama demoslotgacor.pro
akun slot demo
akun slot demo permainan paling top dan garansi imbal balik hasil besar bersama demoslotgacor.pro
akun demo slot
akun demo slot permainan paling top dan garansi imbal balik hasil besar bersama demoslotgacor.pro
slot demo gacor
slot demo gacor permainan paling top dan garansi imbal balik hasil besar bersama situsslotterbaru.net
akun demo slot gacor
akun demo slot gacor permainan paling top dan garansi imbal balik hasil besar bersama situsslotterbaru.net
akun slot demo gacor
akun slot demo gacor permainan paling top dan garansi imbal balik hasil besar bersama situsslotterbaru.net
akun demo slot pragmatic
akun demo slot pragmatic permainan paling top dan garansi imbal balik hasil besar bersama situsslotterbaru.net
akun slot demo pragmatic
akun slot demo pragmatic permainan paling top dan garansi imbal balik hasil besar bersama situsslotterbaru.net
akun slot demo
akun slot demo permainan paling top dan garansi imbal balik hasil besar bersama situsslotterbaru.net
akun demo slot
akun demo slot permainan paling top dan garansi imbal balik hasil besar bersama situsslotterbaru.net
situs slot terbaru
situs slot terbaru permainan paling top dan garansi imbal balik hasil besar bersama situsslotterbaru.net
slot terbaru
slot terbaru permainan paling top dan garansi imbal balik hasil besar bersama situsslotterbaru.net
suara88 permainan paling top dan garansi imbal balik hasil besar bersama suara88.biz
sumo7777 permainan paling top dan garansi imbal balik hasil besar bersama sumo7777.com
supermoney888 permainan paling top dan garansi imbal balik hasil besar bersama supermoney888.biz
teratai88 permainan paling top dan garansi imbal balik hasil besar bersama teratai88.biz
thor88 permainan paling top dan garansi imbal balik hasil besar bersama thor88.biz
togelhk88 permainan paling top dan garansi imbal balik hasil besar bersama togelhk88.net
topjitu88 permainan paling top dan garansi imbal balik hasil besar bersama topjitu88.net
totosloto88 permainan paling top dan garansi imbal balik hasil besar bersama totosloto88.com
trisula888 permainan paling top dan garansi imbal balik hasil besar bersama trisula888.biz
udangbet88 permainan paling top dan garansi imbal balik hasil besar bersama udangbet88.net
via88 permainan paling top dan garansi imbal balik hasil besar bersama via88.biz
virusjp88 permainan paling top dan garansi imbal balik hasil besar bersama virusjp88.net
warga888 permainan paling top dan garansi imbal balik hasil besar bersama warga888.biz
waw88 permainan paling top dan garansi imbal balik hasil besar bersama waw88.biz
winjitu88 permainan paling top dan garansi imbal balik hasil besar bersama winjitu88.net
wisdom88 permainan paling top dan garansi imbal balik hasil besar bersama wisdom88.biz
wnitogel88 permainan paling top dan garansi imbal balik hasil besar bersama wnitogel88.com
yoyo888 permainan paling top dan garansi imbal balik hasil besar bersama yoyo888.biz
validtoto88 permainan paling top dan garansi imbal balik hasil besar bersama validtoto88.com
sule999 permainan paling top dan garansi imbal balik hasil besar bersama sule999.com
sule88 permainan paling top dan garansi imbal balik hasil besar bersama sule88.org
ss888bet permainan paling top dan garansi imbal balik hasil besar bersama ss888bet.com
sia77 permainan paling top dan garansi imbal balik hasil besar bersama sia77.info
seluang88 permainan paling top dan garansi imbal balik hasil besar bersama seluang88.com
satu88 permainan paling top dan garansi imbal balik hasil besar bersama satu88.biz
satu777 permainan paling top dan garansi imbal balik hasil besar bersama satu777.asia
rp88 permainan paling top dan garansi imbal balik hasil besar bersama rp88.biz
rp88 permainan paling top dan garansi imbal balik hasil besar bersama rp88.asia
rp88 permainan paling top dan garansi imbal balik hasil besar bersama rp77.live
qiuqiu88 permainan paling top dan garansi imbal balik hasil besar bersama qiuqiu88.biz
pt88 permainan paling top dan garansi imbal balik hasil besar bersama pt88.org
pt77 permainan paling top dan garansi imbal balik hasil besar bersama pt77.info
produk88 permainan paling top dan garansi imbal balik hasil besar bersama produk88.asia
mt88 permainan paling top dan garansi imbal balik hasil besar bersama mt88.org
mt77 permainan paling top dan garansi imbal balik hasil besar bersama mt77.biz
menang66 permainan paling top dan garansi imbal balik hasil besar bersama menang66.biz
latobet888 permainan paling top dan garansi imbal balik hasil besar bersama latobet888.org
kedai96 permainan paling top dan garansi imbal balik hasil besar bersama kedai96.org
kedai188 permainan paling top dan garansi imbal balik hasil besar bersama kedai188.biz
ids88 permainan paling top dan garansi imbal balik hasil besar bersama ids88.biz
hp88 permainan paling top dan garansi imbal balik hasil besar bersama hp88.org
hp77 permainan paling top dan garansi imbal balik hasil besar bersama hp77.org
gm88 permainan paling top dan garansi imbal balik hasil besar bersama gm88.asia
gm77 permainan paling top dan garansi imbal balik hasil besar bersama gm77.net
final888 permainan paling top dan garansi imbal balik hasil besar bersama final888.org
duit88 permainan paling top dan garansi imbal balik hasil besar bersama duit88.asia
duit168 permainan paling top dan garansi imbal balik hasil besar bersama duit168.biz
divisi88 permainan paling top dan garansi imbal balik hasil besar bersama divisi88.org
dewi500 permainan paling top dan garansi imbal balik hasil besar bersama dewi500.biz
devil88 permainan paling top dan garansi imbal balik hasil besar bersama devil88.info
cuputoto88 permainan paling top dan garansi imbal balik hasil besar bersama cuputoto88.com
cukongbet88 permainan paling top dan garansi imbal balik hasil besar bersama cukongbet88.asia
bom888 permainan paling top dan garansi imbal balik hasil besar bersama bom888.biz
bintaro888 permainan paling top dan garansi imbal balik hasil besar bersama bintaro888.info
askasino88 permainan paling top dan garansi imbal balik hasil besar bersama askasino88.org
999aset permainan paling top dan garansi imbal balik hasil besar bersama 999aset.com
afb77 permainan paling top dan garansi imbal balik hasil besar bersama afb77.biz
aset99 permainan paling top dan garansi imbal balik hasil besar bersama aset99.biz
bendera77 permainan paling top dan garansi imbal balik hasil besar bersama bendera77.biz
bendera888 permainan paling top dan garansi imbal balik hasil besar bersama bendera888.com
coco88 permainan paling top dan garansi imbal balik hasil besar bersama coco88.org
cuma77 permainan paling top dan garansi imbal balik hasil besar bersama cuma77.biz
cuma88 permainan paling top dan garansi imbal balik hasil besar bersama cuma88.org
dwv88 permainan paling top dan garansi imbal balik hasil besar bersama dwv88.org
fafajp88 permainan paling top dan garansi imbal balik hasil besar bersama fafajp88.com
gemar88 permainan paling top dan garansi imbal balik hasil besar bersama gemar88.biz
gocap88 permainan paling top dan garansi imbal balik hasil besar bersama gocap88.info
gocaptoto permainan paling top dan garansi imbal balik hasil besar bersama gocaptoto.asia
hakabet88 permainan paling top dan garansi imbal balik hasil besar bersama hakabet88.com
hwtoto88 permainan paling top dan garansi imbal balik hasil besar bersama hwtoto88.org
ina77 permainan paling top dan garansi imbal balik hasil besar bersama ina77.biz
ina88 permainan paling top dan garansi imbal balik hasil besar bersama ina88.info
jingga8888 permainan paling top dan garansi imbal balik hasil besar bersama jingga8888.com
juragan777 permainan paling top dan garansi imbal balik hasil besar bersama juragan777.asia
kastil77 permainan paling top dan garansi imbal balik hasil besar bersama kastil77.info
kebo888 permainan paling top dan garansi imbal balik hasil besar bersama kebo888.biz
kkwin77 permainan paling top dan garansi imbal balik hasil besar bersama kkwin77.com
kokoslot88 permainan paling top dan garansi imbal balik hasil besar bersama kokoslot88.asia
luckydf88 permainan paling top dan garansi imbal balik hasil besar bersama luckydf88.org
microstar888 permainan paling top dan garansi imbal balik hasil besar bersama microstar888.biz
monperatoto88 permainan paling top dan garansi imbal balik hasil besar bersama monperatoto88.com
mpo1122 permainan paling top dan garansi imbal balik hasil besar bersama mpo1122.biz
mpo122 permainan paling top dan garansi imbal balik hasil besar bersama mpo122.biz
mpopelangi88 permainan paling top dan garansi imbal balik hasil besar bersama mpopelangi88.com
pamanslot88 permainan paling top dan garansi imbal balik hasil besar bersama pamanslot88.biz
panel88 permainan paling top dan garansi imbal balik hasil besar bersama panel88.org
paragon77 permainan paling top dan garansi imbal balik hasil besar bersama paragon77.biz
paragon888 permainan paling top dan garansi imbal balik hasil besar bersama paragon888.info
pion77 permainan paling top dan garansi imbal balik hasil besar bersama pion77.biz
prada88 permainan paling top dan garansi imbal balik hasil besar bersama prada88.asia
prada888 permainan paling top dan garansi imbal balik hasil besar bersama prada888.com
qqslot88slot permainan paling top dan garansi imbal balik hasil besar bersama qqslot88slot.com
rejekibet88 permainan paling top dan garansi imbal balik hasil besar bersama rejekibet88.com
rezekibet88 permainan paling top dan garansi imbal balik hasil besar bersama rezekibet88.org
sensa77 permainan paling top dan garansi imbal balik hasil besar bersama sensa77.biz
sensa888 permainan paling top dan garansi imbal balik hasil besar bersama sensa888.biz
singajp88 permainan paling top dan garansi imbal balik hasil besar bersama singajp88.com
sr77 permainan paling top dan garansi imbal balik hasil besar bersama sr77.org
sr88 permainan paling top dan garansi imbal balik hasil besar bersama sr88.org
surya77 permainan paling top dan garansi imbal balik hasil besar bersama surya77.biz
surya88 permainan paling top dan garansi imbal balik hasil besar bersama surya88.asia
tajir77 permainan paling top dan garansi imbal balik hasil besar bersama tajir77.info
tajir88 permainan paling top dan garansi imbal balik hasil besar bersama tajir88.biz
toto122 permainan paling top dan garansi imbal balik hasil besar bersama toto122.com
toto123 permainan paling top dan garansi imbal balik hasil besar bersama toto123.biz
uangvip88 permainan paling top dan garansi imbal balik hasil besar bersama uangvip88.com
wajik77 permainan paling top dan garansi imbal balik hasil besar bersama wajik77.asia
777neko permainan paling top dan garansi imbal balik hasil besar bersama 777neko.org
88judi permainan paling top dan garansi imbal balik hasil besar bersama 88judi.net
99judi permainan paling top dan garansi imbal balik hasil besar bersama 99judi.org
abcslot88 permainan paling top dan garansi imbal balik hasil besar bersama abcslot88.asia