mosGILT controls innate immunity and germ cell development in Anopheles gambiae

Gene‑edited mosquitoes lacking a gamma‑interferon‑inducible lysosomal thiol reductase‑like protein, namely ( mosGILT null ) have lower Plasmodium infection, which is linked to impaired ovarian development and immune activa‑ tion. The transcriptome of mosGILT null Anopheles gambiae was therefore compared to wild type (WT) mosquitoes by RNA‑sequencing to delineate mosGILT‑dependent pathways. Compared to WT mosquitoes, mosGILT null A. gambiae demonstrated altered expression of genes related to oogenesis, 20‑hydroxyecdysone synthesis, as well as immune‑ related genes. Serendipitously, the zero population growth gene, zpg , an essential regulator of germ cell development was found to be one of the most downregulated genes in mosGILT null mosquitoes. These results provide a crucial missing link between two previous studies on the role of zpg and mosGILT in ovarian development. This study further demonstrates that mosGILT has the potential to serve as a target for the biological control of mosquito vectors and to influence the Plasmodium life cycle within the vector.


Introduction
Malaria remains one of the deadliest vector-borne diseases, with over 247 million cases and approximately 619 thousand deaths reported globally in 2021 [1].It is caused by apicomplexan protozoa of the genus Plasmodium and transmitted by female Anopheles mosquitoes.During its journey through the vector, Plasmodium is profoundly influenced by mosquito factors that either facilitate or suppress parasite development [2][3][4][5].In the mosquito midgut, Plasmodium development begins with the differentiation of gametocytes into gametes, followed by the formation of zygotes, ookinetes, and then oocysts.Parasite numbers suffer great losses, especially during traversal of the midgut epithelium [6,7].This midgut bottleneck of the Plasmodium life cycle is, therefore, an optimal target for blocking malaria transmission [8][9][10][11][12][13].Considering the increase of insecticide resistance in Anopheles [5,[14][15][16], a better knowledge of how mosquito-derived factors influence midgut-stage parasites is important for developing novel strategies for malaria intervention.
Blood-feeding is important for the female mosquito to obtain the nutrients required for egg development.Studies suggest that Plasmodium maturation in the midgut is influenced by concurrent mosquito reproductive processes, yet the underlying mechanisms remain poorly understood [17,18].Through CRISPR/Cas9 gene editing, our group demonstrated that disruption of a mosquito gamma-interferon-inducible lysosomal thiol reductaselike protein-mosGILT-caused underdeveloped ovaries with a reduction in egg development and the production of ecdysteroid hormone 20-hydroxyecdysone (20E), and vitellogenin (Vg) [19].Interestingly, mosGILT-deficient (mosGILT null ) Anopheles gambiae also showed significant suppression of both Plasmodium berghei and Plasmodium falciparum infection, as measured by reduced oocyst intensity and prevalence [19].In parallel, another study showed that Anopheles mosquitoes with a disrupted zero-population growth-zpg-gene also exhibited underdeveloped ovaries that produced low amounts of 20E and were less permissive to P. falciparum infection [20].Though both mosGILT null and zpg null mosquitoes exhibited somewhat similar phenotypes, it is unclear whether their mechanisms of action are similar or distinct.In this study, we elucidate the mosGILT-dependent global transcriptional responses in A. gambiae.The information here may provide novel insights into the A. gambiae reproductive system network, as well as innate responses associated with Plasmodium infection.

Transcriptome analysis of mosGILT null A. gambiae
To examine whether the molecular responses to a blood meal were altered in mosGILT null A. gambiae, RNA sequencing (RNA-seq) was performed.mosGILT null and Wild type (WT, control) mosquitoes were fed on naïve mice and then separated into cardboard containers.The whole body, ovaries, and fat bodies were collected from blood-fed mosquitoes 24 h after completion of the meal.RNA-seq was carried out with total RNA extracted from the whole body, ovaries, or fat bodies from A. gambiae using Illumina sequencing technology.A total of 16 samples (from control and mosGILT null A. gambiae) were sequenced (Supplementary Fig. 1).After sequencing, each library generated around 10-28 million reads and approximately 78.0-99.7% of the reads in each library mapped to the A. gambiae genome (Supplementary Figs. 2, 3, and 4).
Principal component analysis of mosquito libraries examined the clustering of data/samples after a blood meal.All biological replicates of control and mosGILT null samples were distributed into two distinct groups (Fig. 1A).Volcano plots indicated that there were differentially expressed genes (DEGs) in all three groups (Fig. 1B), with whole-body samples having the maximum number of DEGs (n = 5749), followed by the ovaries and fat bodies.In all the three samples analyzed, 108 DEGs were common to all three groups (Fig. 2A).Of particular interest were 4991 DEGs in the ovaries of mosGILT null mosquitoes which are potentially helpful to elucidate the role of mosGILT in A. gambiae development and Plasmodium infection.

DEG analysis reveals that mosGILT has a role in nucleocytoplasmic transport
Since female mosGILT null mosquitoes display defects in ovarian development and show refractoriness to Plasmodium, we first focused on the nature of genes that showed differential expression patterns in the ovarian tissue.Among the 4991 DEGs, 117 genes had a more than a 100-fold increase in mRNA abundance compared to the WT controls, and 97 genes showed at least a 100fold lower mRNA abundance.Several groups of genes showed altered expression patterns in oxidative phosphorylation and nucleocytoplasmic transport based on Gene Ontology (GO) annotation (Figs.2B, and 3A, B, C, D).Nucleocytoplasmic transport (75 genes, p-value 2.60E-13) is a process crucial for cell function (Supplementary Table 1 and Supplementary Fig. 5).Depletion  3 A, B and C Bubble chart showing enrichment of Gene Ontology (GO) terms of the differentially expressed genes in the whole body, fat body, and ovary of mosGILT null mosquitoes.The biological processes, cellular components, and molecular functions are shown in blue, magenta and orange respectively.The -Log10 (p-value) for each process is shown on the X-axis.D The figure presents the Manhattan plot showcasing enrichment analysis of differentially expressed genes in the ovary of mosGILT null mosquitoes performed using g:Profiler.The x-axis represents functional GO terms grouped and color-coded based on data sources.The y-axis displays adjusted enrichment p-values on a negative log10 scale of mammalian GILT (IFI30) in immune cells can inhibit cytokine production and the release of high mobility group box 1 (HMGB1), a nonhistone nuclear protein.
The cytosolic translocation of the nuclear HMGB1 protein is then activated and associated with increased autophagy in GILT-depleted fibroblasts [21,22].These findings suggest that mammalian GILT is involved in the homeostatic regulation of oxidative stress.In the absence of mosGILT, nucleocytoplasmic transport pathways are upregulated which can influence diverse cellular responses, such as autophagy (Fig. 2B).Indeed, autophagy negatively regulates Vg production in gonotrophic cycles in Aedes mosquitoes, which is consistent with the changes in Vg production as shown by our group previously [23].

Whole body analysis of A. gambiae confirms the role of mosGILT in development and metabolism
Our previous studies demonstrated that mosGILT is implicated in embryonic and ovarian development and shows refractoriness to Plasmodium infection [19].To fully understand how female mosGILT null mosquitoes show resistance to Plasmodium, we also examined the genes that showed differential expression patterns in the whole body.Among the 5749 DEGs, 76 genes had a more than a 100-fold decrease in expression compared to the WT control, and 7 genes showed at least 100fold upregulation.Similar to the ovaries, nucleocytoplasmic transport pathways were affected in the whole body of mosGILT null mosquitoes.The most downregulated proteins included AGAP007823 (meiosis arrest female protein 1, fold change: -1000000, p-value < 0.0001) and AGAP008738 (eukaryotic translation initiation factor 4E transporter, fold change -14625, p-value 0.0001) which is involved in protein translation and oogenesis (Fig. 4A).
The most overexpressed genes included AGAP010658 (fold change 169, p-value < 0.00001) which encodes a homolog of the hexameric 2 beta protein or larval serum protein 1 which serves as a major storage protein.AGAP010658 antagonizes Plasmodium berghei oocyst development and silencing of this gene by RNAi results in a significant increase in oocyst infection intensities and melanization of ookinetes [24].Other highly expressed genes were AGAP000046 (a sugar transporter, fold change 160, p-value < 00001) and AGAP002518 (delta l-pyrroline-5-carboxylate synthetase, fold change 102, p-value < 0.00001).While transcripts of AGAP004203 were decreased (fold change -343), another Vg-like protein was overexpressed (AGAP008369 fold change 102, p-value < 0.00001).Among the major downregulated genes were AGAP002134 (vitelline membrane Fig. 4 Heat Map analysis of key differentially expressed genes in the ovary (A), and fat body (B) of mosGILT null mosquitoes.C Heat map analysis of key differentially expressed genes that are related to germ cell development in mosGILT null mosquitoes.D STRING functional gene analysis of zpg-related genes: Each node represents a protein product of an investigated gene.Lines denote protein-protein interactions, with each edge representing an interaction including either physical or functional associations component which is required for the proper assembly of the eggshell) AGAP003545 (maternal effect protein oskar), and AGAP007767 (mosquito-specific bZIP1 transcription factor) (Fig. 4A).Indeed, silencing of bZIP1 in the embryos of Aedes aegypti resulted in embryonic death [25].These findings indicate that mosGILT is involved in both development and metabolism.

Effect of mosGILT on metabolic responses in the A. gambiae fat body
Our previous study showed that mosGILT is expressed in the ovaries and fat body [19].In insects, the fat body plays a central role in energy metabolism.After a blood meal, the fat body, which is the mosquito fat storage organ, initiates Vg synthesis.The fat body releases Vg into the hemolymph, and it is then taken up by developing oocytes via receptor-mediated endocytosis.Since mosGILT null mosquitoes showed lower Vg mRNA levels, we examined the transcriptional response in the fat body.Our analysis shows that in the fat body, gene transcripts of metabolic pathways (gene count: 93, p-value 6.30E-14) and pathways related to carbohydrate metabolism (pentose and glucuronate interconversions, fructose and mannose metabolism, galactose metabolism, biosynthesis of nucleotide sugars, glycolysis/gluconeogenesis and starch, and sucrose metabolism) demonstrated major changes in the mosGILT null mosquitoes.Among the most altered DEGs, AGAP029067 (lipid-binding protein) transcripts were enriched by a 194-fold while transcripts of AGAP008547, a gene with a predicted function in carbohydrate metabolism showed nearly 160-fold enrichment in the fat body of mosGILT null mosquitoes (Fig. 4A).Blood-meal ingestion, therefore, not only induces the expression of genes related to egg development but also increased the expression of genes related to oxidative stress.
In mammals, GILT is expressed in lysosomes and is required for optimal activity of cysteine proteases, such as cathepsin in macrophages [26,27].In mosquitoes, the fat body plays a central role in metabolism and nutrient storage, mirroring the functions of the liver in vertebrates.In response to blood-feeding, the Aedes fat body secretes a cathepsin-b-like thiol protease [28].Interestingly, in our analysis of mosGILT null A. gambiae, the most downregulated genes include proteases such as AGAP004534 (cathepsin B -like serine protease, fold change -18, p-value 0.0001), and AGAP004700 (trypsin-like serine protease, fold change -20, p-value 0.0003) (Fig. 4B).

mosGILT regulates germ cell maintenance by controlling the expression of zero population growth (zpg)
A. gambiae females express zero population growth (zpg, AGAP006241), an essential gene for germ cell development.Genetic ablation of zpg in mosquitoes leads to severely underdeveloped ovaries and impaired egg production after a blood meal [20,29,30] similar to the mosGILT null mosquitoes.zpg and mosGILT are important for Plasmodium infection since mosquitoes lacking either gene show fewer oocysts after feeding on P. falciparum-infected blood compared to sibling controls.In the ovaries of mosGILT null A. gambiae, zpg was strongly downregulated (fold change -8130, p-value 0.0001) which suggests mosGILT acts upstream of zpg or influences the molecular pathway that regulates zpg expression (Fig. 4C).We further performed protein network analysis (Fig. 4D, Supplementary Fig. 6) to determine whether mosGILT affected the expression of other genes that are associated with zpg.STRING analysis showed the genes that are in the zpg network may be involved in germ cell development.Lack of mosGILT affected the expression of AGAP000427 (fold change -26814, p-value > 0.00001), which is in the zpg network and is a putative Vg receptor.In the ovaries of mosGILT null mosquitoes, another zpg-related protein, oskar (AGAP003545) was strongly downregulated (fold change -8162, p-value > 0.00001) (Fig. 4C).Oskar is required for the formation of the germplasm in Drosophila [31].Germplasm promotes the formation and specification of the primordial germ cells, the first cell lineage to form in the fertilized embryo [32].Our genetic analyses of mosGILT null mosquitoes reveal underdeveloped ovaries and a decrease in 20E production.Here, we have analyzed genes involved in insect hormone biosynthesis and found that five genes (AGAP001038, AGAP001039, AGAP004665, AGAP005992, AGAP000284) involved in 20E biosynthesis were downregulated in mosGILT null mosquitoes (Fig. 5).Signaling of 20E is initiated by binding of 20E to the heterodimer ecdysone receptor (EcR)/ultraspiracle (USP) complex, which triggers a transcriptional response that promotes egg development.Our results show that mosGILT null also decreases the expression of ultraspiracle (USP, AGAP002095, fold change -4, p-value < 0.00001) (Fig. 4C).The complex formed after 20E binding to EcR/USP upregulates the expression of an array of 20E-regulated genes which in turn controls the Vg transcriptional program [33].These results show that mosGILT regulates 20E levels, as well as the downstream transcriptional response.

Immune pathway alteration in mosGILT null mosquitoes
GILT was first identified as an IFNγ-inducible, lysosomal thiol reductase in mammals [37].GILT is also key to antigen processing in mammals but in other organisms it is part of the innate immune response [38].Thioester-containing protein 1 (TEP1) is an important component of the A. gambiae innate immune system and plays a major role against Plasmodium parasites [39,40].TEP is an antimicrobial protein family that acts like the human complement protein C3 and can damage the cell membranes of pathogens.In our study, we observed an increase in TEP1 activity in mosGILT null mosquitoes [19].We also showed that strong refractoriness to Plasmodium infection in mosGILT null mosquitoes was related to increased TEP1 activity.Therefore, we examined immune genes that are differentially expressed in mosGILT null mosquitoes to better understand how TEP1 activity is influenced.
In our analysis, we found that AGAP028725 (SPCLIP1) was upregulated in the ovary (fold change 73, p-value 0.0001) and whole body (fold change 12, p-value 0.0001) of mosGILT null mosquitoes (Fig. 6).Clip domain serine proteases (CLIPs) play an important role in innate immunity in mosquitoes.These proteases are present in hemolymph and are involved in Plasmodium immunity, antimicrobial peptide synthesis and melanization.Previously it was shown that complement C3-like protein TEP1 mediated killing of Plasmodium is dependent on the CLIP-domain serine protease homolog SPCLIP1 [41,42].SPCLIP1 is required for the accumulation of TEP1 on microbial surfaces, a reaction that leads to lysis of malaria parasites [41,42].
In addition, we also observed increase in defensin transcripts, an anti-microbial peptide (AGAP011294, DEF1).Def1 expression was upregulated in the ovary (fold change 255, p-value 0.0001) and whole body (fold change 3, p-value 0.0029) of mosGILT null mosquitoes (Fig. 6).Transcripts of another anti-microbial peptide gene, cercopin 1 (AGAP000693) was also increased in the ovary (fold change 300, p-value 0.0004) and whole body of the mosGILT null mosquitoes (fold change 2, p-value 0.0235), and cercopin 3 (AGAP000694) expression was upregulated in the whole body (fold change 3, p-value 0.0016).Toll-like receptors (TLRs) are evolutionarily conserved pattern recognition receptor proteins with diverse biological functions.In Drosophila melanogaster, Toll plays an important role in formation during embryogenesis and in immunity.In this study we observed that TLR9 expression the mosGILT null mosquito was upregulated in the whole body (AGAP006974, fold change 2, p-value 0.0095) Though there is not much information on this gene in A. gambiae, in the silk moth, Bombyx mori, TLR9 is involved in innate immunity and functions as a pattern recognition receptor that binds lipopolysaccharide.

Discussion
To further understand Plasmodium infection of the vector, and develop new mosquito and malaria control strategies, it is important to decipher the molecular mechanisms related to mosquito reproduction and immunity.Recently our lab uncovered a critical role for mosGILT in A. gambiae ovary development and Plasmodium infection [19].The mosGILT null mosquitoes displayed underdeveloped ovaries and reproductive defects including impaired production of 20E and Vg while showing an increased TEP1-mediated immune response against Plasmodium ookinetes [19].However, there is a lack of knowledge regarding how mosGILT functions in A. gambiae, its interacting partners, and the implicated signaling pathways.Here we characterize the transcriptional response in mosGILT null mosquitoes, including the ovaries, whole body, and fat body, and delineate how mosGILT is regulating reproduction and immune responses.
GILT is a thioredoxin-related oxidoreductase which is the only known lysosomal thiol reductase [27].Mammalian GILT is related to the thioredoxin family of oxidoreductases, characterized by a CXXC motif in the active site and an enzymatic mechanism in which the pair of active site cysteine residues cooperate to reduce substrate disulfide bonds [45][46][47][48].In Anopheles mosquitoes and other insect genera, the active motif is CXXS, where the second cysteine is replaced with serine (Fig. 7A  and B).This change in the active site of mosGILT hinders the thioredoxin reaction of the mixed intermolecular disulfide.The CXXS containing thioredoxins can interact with their substrate proteins in a more stable manner since there is a covalent complex formed between thioredoxins and their target.However, some studies report that at least some CXXS-containing forms of thioredoxin-fold proteins retain catalytic activity albeit at a Fig. 6 Heat map analysis of key differentially expressed genes related to immune pathways in the whole body (WB), fat body (FB) and ovary (OV) of mosGILT null mosquitoes lower level [49,50].It is therefore possible that mosGILT has partial thioredoxin activity that is critical to the transcriptional response that leads to oogenesis and ovarian development.
Vector reproduction is closely associtaed with immune defense pathways, as these two processes are simultaneously initiated after a blood meal [51][52][53].Identification of mosquito proteins that regulate cross-talk between reproduction and could serve as targets for vector-control strategies.Our previous study demonstrated that mosGILT plays an essential role in the reproduction of A. gambiae [19].In this study we show that mosGILT influences complement activity and upregulates defense-related genes, which are associated with anti-Plasmodium responses.In addition, mosGILT regulates TEP1 activity by altering SPCLIP1 expression, which can influence Plasmodium survival in mosquitoes.mosGILT null mosquito refractoriness to microbes was also suggested by the increase in genes encoding antimicrobial peptides and innate immune molecules.mosGILT may therefore be a potential target for new malaria control strategies, such as gene drive -propelled transgenes targeting mosGILT, that would induce mosquito sterility and also reduce the overall malaria parasite burden.
Development and immunity are highly demanding processes and, thus, a balance must be maintained between the two [54].Previously it was suggested that there is a trade-off between reproductive fitness and immune defense in mosquitoes including A. gambiae and activation of immune pathways can cause a reduction in fecundity by inducing the resorption of ovarian follicles during oocyte development [55][56][57][58].In another example, activation of immune signaling leads to decreased in reproductive capacity of D. melanogaster [59].These observations suggest that in the absence of mosGILT and a reduction in ovary development, there may be an increase in the immune response.This is also supported by a recent study where it was shown that 20E binding to the ecdysone receptor (EcR) directly induces a Pirklike inhibitor, affecting the immune deficiency (IMD) pathway in Aedes mosquitoes [57].Pirk downregulates the IMD pathway by inhibitng the nuclear translocation of the NF-κB transcriptional factor Rel2 and this downregulation of the immune response is vital for maintaining mosquito fertility.We hypothesize that the effects of downregulating ecdysone on mosquito immunity can be complex.We acknowledge that in some studies 20E is shown to have immunomodulatory effects, and it can enhance the mosquitoes immune response [52,[60][61][62].Our results with mosGILT null mosquitoes show that the influence of 20E on immunity may be part of a broader regulatory network.Overall our results suggest that mos-GILT may be important for resource allocation in respect to both development and immunity.In the absence of mosGILT, resources are moved to immune pahways and there is reduction in reproductive capacity.Overall, our results suggest tha mosquito immunity is multilayered and hierarchical, where downregulation of the mosGILT in Anopheles mosquitoes is sufficient to activate the immune system even in the absence of 20E.
mosGILT mutants have underdeveloped ovaries that produced low amounts of 20E, and were significantly less permissive to P. falciparum infection [19].This is somewhat similar to Plasmodium infection in zpg mutant mosquitoes, a gene that is essential for germ cell development [20,30].Though both mutant mosquitoes provide an association between 20E signaling and intensity of Plasmodium infection, it is not known if they act together or independently.Our studies now suggest that mos-GILT may affect several genes that are associated with zpg.Since zpg and other germ cell development genes are considered a prime target for gene drive [63][64][65], our results suggest mosGILT may be more attractive target for gene drive of A. gambiae.
mosGILT plays an important role in reproduction and host defenses in A. gambiae.The elucidation of mos-GILT as a regulator of zpg and involvement in cross-talk between reproduction and immunity could provide new targets for the development of disease-control strategies (Fig. 8).These findings can potentially lead to new preventive strategies for malaria based on Anopheles mos-GILT and also potentially suggest a role for arthropod GILT-like proteins for other vector-borne diseases.

Generation of transgenic mosGILT mutant mosaic
Anopheles mosquitoes CRISPR/Cas9 mediated A. gambiae mosGILT somatic KO mutants were generated using the A. gambiae docking line X1 as described in Yang et al., 2020 [19].In brief, the mosGILT-gRNA female mosquitoes were crossed with the homozygous male Vasa Cas9 transgenic strain (vasa2-Cas9; with green/yellow fluorescence in the eyes;) [67,68].The progeny of this cross was screened at the larvae stage with both RFP (red fluorescence marker) and YFP (appeared as green fluorescence through GFPchannel under the fluorescence microscope) signals was gRNA/Cas9 trans-heterozygous.Molecular screening of these trans-heterozygotes by PCR and sequencing was done as described to confirm the mosGILT null [19].

RNA isolation
Total RNA was extracted from the whole body, ovaries, or fat bodies from mosquitoes.Trizol was added to the tissues and RNA was isolated using a using a combination method which utilizes TRIzol (Thermo Fisher Scientific) and RNeasy Mini Kit (QIAGEN) as described before [69].

RNA-seq
RNA was submitted for library preparation using TruSeq (Illumina, San Diego, CA, USA) and sequenced using Illumina HiSeq 2500 by paired-end sequencing at the Yale Center for Genome Analysis (YCGA).

Differential gene expression analysis
All the RNA-seq analyses including alignment, quantitation, normalization, and differential gene expression analyses were performed using Partek Genomics Flow software (St.Louis, MO, USA).Specifically, RNA-seq data were trimmed and aligned to the Anopheles PEST genome, version P4.14. with associated annotation file using STAR (v2.7.3a) [70].The aligned reads were quantified using the Partek E/M algorithm [71] and the subsequent steps were performed on gene-level annotation followed by total count normalization.The gene-level data were normalized by dividing the gene counts by the total number of reads followed by the addition of a small offset (0.0001).

Principal components analysis
Principal components analysis (PCA) was performed using default parameters for the determination of the component number, with all components contributing equally in Partek Flow.Volcano plot Hierarchal clustering was performed on the genes that were differentially expressed across the conditions (P < 0.05, fold change ≥ 2 for each comparison).

PPI network and Cluster analysis
PPI network was generated by STRING online database and visualized by Cytoscape software [75,76].K means clustering with 10 clusters was selected to generate cluster with edges.
• fast, convenient online submission • thorough peer review by experienced researchers in your field • rapid publication on acceptance • support for research data, including large and complex data types • gold Open Access which fosters wider collaboration and increased citations maximum visibility for your research: over 100M website views per year

•
At BMC, research is always in progress.

Learn more biomedcentral.com/submissions
Ready to submit your research Ready to submit your research ?Choose BMC and benefit from: ? Choose BMC and benefit from:

Fig. 1 A
Fig. 1 A Principal component analysis of the RNA-Seq data.The whole body (WB), fat body (FB) and ovary (OV) samples were collected at 24 h from blood-fed (BF) wild-type (WT) and mosGILT null Anopheles gambiae.B Volcano plot analysis of differentially expressed genes (DEGs) of whole body (WB), fat body (FB) and ovary (OV) between wild type (control) and mosGILT null Anopheles gambiae.DEGs with False Discovery Rate (FDR) p-values of < 0.05 and log 2 fold changes > 1 are presented

Fig. 2 A
Fig. 2 A Venn diagrams showing the number of differentially expressed genes (DEGs) between whole body (WB), fat body (FB) and ovary (OV) tissues of blood-fed Anopheles gambiae.The numbers in the overlapping areas indicate genes that were common to both or to all three different time points.The number of DEGs was highest between the whole body and ovary samples.B Bubble chart and circle graph showing enriched pathways in the whole body, fat body, and ovary of mosGILT null mosquitoes by KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis.Gene ontology annotation and pathway enrichment analysis of the top differentially expressed genes (DEGs) in the whole body, fat body, and ovary are shown.Enriched functional ontology terms are shown on the left.The whole body, fat body, and ovary are represented by blue, red, and orange colors respectively.The Log 10 (p-value) is shown on the X-axis

Fig. 7
Fig. 7 Structural and sequence analysis of mosGILT (A) The mosGILT structural model is created by Alpha Fold.Structural alignment of mosGILT homologs from Anopheles, Aedes, Drosophila and human GILT.The active site motif is shown in the yellow in the inset.B Sequence alignment of Anopheles, Aedes, Drosophila and Ixodes GILT.Protein sequence of mosGILT homologs from Anopheles gambiae (AGAP004551), Aedes aegypti (XP_021707151.1 GILT-like protein 1), Drosophila melanogaster (NP_650287.3Gamma-interferon-inducible lysosomal thiol reductase 1), Ixodes scapularis (XP_029851151.3GILT-like protein 1).Active site residues are shown in the inset box

Fig. 8
Fig. 8 Diagrammatic representation of mosGILT's role in immunity and development