Serine protease-related proteins in the malaria mosquito, Anopheles gambiae

A single-cell atlas of the Culex tarsalis midgut during West Nile virus infection

The mosquito midgut functions as a key interface between pathogen and vector. However, studies of midgut physiology and associated virus infection dynamics are scarce, and in Culex tarsalis - an extremely efficient vector of West Nile virus (WNV) - nonexistent. We performed single-cell RNA sequencing on Cx. tarsalis midguts, defined multiple cell types, and determined whether specific cell types are more permissive to WNV infection. We identified 20 cell states comprised of 8 distinct cell types, consistent with existing descriptions of Drosophila and Aedes aegypti midgut physiology. Most midgut cell populations were permissive to WNV infection. However, there were higher levels of WNV RNA (vRNA) in enteroendocrine cells and cells enriched for mitochondrial genes, suggesting enhanced replication in these populations. In contrast, proliferating intestinal stem cell (ISC) populations had the lowest levels of vRNA, a finding consistent with studies suggesting ISC proliferation in the midgut is involved in viral control. Notably, we did not detect significant WNV-infection induced upregulation of canonical mosquito antiviral immune genes (e.g., AGO2, R2D2, etc.) at the whole-midgut level. Rather, we observed a significant positive correlation between immune gene expression levels and vRNA in individual cells, suggesting that within midgut cells, high levels of vRNA may trigger antiviral responses. Our findings establish a Cx. tarsalis midgut cell atlas, and provide insight into midgut infection dynamics of WNV by characterizing cell-type specific enhancement/restriction of, and immune response to, infection at the single-cell level.

The expanded immunoregulatory protease network in mosquitoes is governed by gene co-expression

Serine protease cascades regulate key innate immune responses. In mosquitoes, these cascades involve clip-domain serine proteases and their non-catalytic homologs (CLIPs), forming a complex network whose make-up and structural organization is not fully understood. This study assessed the impact of 85 CLIPs on humoral immunity in Anopheles gambiae. By coupling RNAi with assays measuring antimicrobial activity and melanization, we identified 27 CLIPs as immunoregulators that together form two distinct subnetworks. CLIPs regulating antimicrobial activity were found to control infection resistance, as knockdowns reduced bacterial load and improved survival. Furthermore, our analysis of CLIP gene expression unveiled a novel immunoregulatory mechanism reliant on protease baseline co-expression rather than infection-induced upregulation. These findings underscore that despite its complexity mosquito immune regulation may be targeted for malaria interventions.

The Set of Serine Peptidases of the Tenebrio molitor Beetle: Transcriptomic Analysis on Different Developmental Stages

Serine peptidases (SPs) of the chymotrypsin S1A subfamily are an extensive group of enzymes found in all animal organisms, including insects. Here, we provide analysis of SPs in the yellow mealworm Tenebrio molitor transcriptomes and genomes datasets and profile their expression patterns at various stages of ontogeny. A total of 269 SPs were identified, including 137 with conserved catalytic triad residues, while 125 others lacking conservation were proposed as non-active serine peptidase homologs (SPHs). Seven deduced sequences exhibit a complex domain organization with two or three peptidase units (domains), predicted both as active or non-active. The largest group of 84 SPs and 102 SPHs had no regulatory domains in the propeptide, and the majority of them were expressed only in the feeding life stages, larvae and adults, presumably playing an important role in digestion. The remaining 53 SPs and 23 SPHs had different regulatory domains, showed constitutive or upregulated expression at eggs or/and pupae stages, participating in regulation of various physiological processes. The majority of polypeptidases were mainly expressed at the pupal and adult stages. The data obtained expand our knowledge on SPs/SPHs and provide the basis for further studies of the functions of proteins from the S1A subfamily in T. molitor.

Characteristics of Trypsin genes and their roles in insecticide resistance based on omics and functional analyses in the malaria vector Anopheles sinensis

Trypsin is one of the most diverse and widely studied protease hydrolases. However, the diversity and characteristics of the Trypsin superfamily of genes have not been well understood, and their role in insecticide resistance is yet to be investigated. In this study, a total of 342 Trypsin genes were identified and classified into seven families based on homology, characteristic domains and phylogenetics in Anopheles sinensis, and the LY-Domain and CLECT-Domain families are specific to the species. Four Trypsin genes, (Astry2b, Astry43a, Astry90, Astry113c) were identified to be associated with pyrethroid resistance based on transcriptome analyses of three field resistant populations and qRT-PCR validation, and the knock-down of these genes significantly decrease the pyrethroid resistance of Anopheles sinensis based on RNAi. The activity of Astry43a can be reduced by five selected insecticides (indoxacarb, DDT, temephos, imidacloprid and deltamethrin); and however, the Astry43a could not directly metabolize these five insecticides, like the trypsin NYD-Tr did in earlier reports. This study provides the overall information frame of Trypsin genes, and proposes the role of Trypsin genes to insecticide resistance. Further researches are necessary to investigate the metabolism function of these trypsins to insecticides.

Comprehensive Proteomics Analysis of the Hemolymph Composition of Sugar-Fed Aedes aegypti Female and Male Mosquitoes

In arthropods, hemolymph carries immune cells and solubilizes and transports nutrients, hormones, and other molecules that are involved in diverse physiological processes including immunity, metabolism, and reproduction. However, despite such physiological importance, little is known about its composition. We applied mass spectrometry-based label-free quantification approaches to study the proteome of hemolymph perfused from sugar-fed female and male Aedes aegypti mosquitoes. A total of 1403 proteins were identified, out of which 447 of them were predicted to be extracellular. In both sexes, almost half of these extracellular proteins were predicted to be involved in defense/immune response, and their relative abundances (based on their intensity-based absolute quantification, iBAQ) were 37.9 and 33.2%, respectively. Interestingly, among them, 102 serine proteases/serine protease-homologues were identified, with almost half of them containing CLIP regulatory domains. Moreover, proteins belonging to families classically described as chemoreceptors, such as odorant-binding proteins (OBPs) and chemosensory proteins (CSPs), were also highly abundant in the hemolymph of both sexes. Our data provide a comprehensive catalogue of A. aegypti hemolymph basal protein content, revealing numerous unexplored targets for future research on mosquito physiology and disease transmission. It also provides a reference for future studies on the effect of blood meal and infection on hemolymph composition.

Wolbachia infection-responsive immune genes suppress Plasmodium falciparum infection in Anopheles stephensi

Wolbachia, a maternally transmitted symbiotic bacterium of insects, can suppress a variety of human pathogens in mosquitoes, including malaria-causing Plasmodium in the Anopheles vector. However, the mechanistic basis of Wolbachia-mediated Plasmodium suppression in mosquitoes is not well understood. In this study, we compared the midgut and carcass transcriptomes of stably infected Anopheles stephensi with Wolbachia wAlbB to uninfected mosquitoes in order to discover Wolbachia infection-responsive immune genes that may play a role in Wolbachia-mediated anti-Plasmodium activity. We show that wAlbB infection upregulates 10 putative immune genes and downregulates 14 in midguts, while it upregulates 31 putative immune genes and downregulates 15 in carcasses at 24 h after blood-fed feeding, the time at which the Plasmodium ookinetes are traversing the midgut tissue. Only a few of these regulated immune genes were also significantly differentially expressed between Wolbachia-infected and non-infected midguts and carcasses of sugar-fed mosquitoes. Silencing of the Wolbachia infection-responsive immune genes TEP 4, TEP 15, lysozyme C2, CLIPB2, CLIPB4, PGRP-LD and two novel genes (a peritrophin-44-like gene and a macro domain-encoding gene) resulted in a significantly greater permissiveness to P. falciparum infection. These results indicate that Wolbachia infection modulates mosquito immunity and other processes that are likely to decrease Anopheles permissiveness to Plasmodium infection.

Insight into the structural hierarchy of the protease cascade that regulates the mosquito melanization response

Serine protease cascades regulate important insect immune responses, including melanization and Toll pathway activation. In the context of melanization, central components of these cascades are clip domain serine proteases (CLIPs) including the catalytic, clip domain serine proteases (cSPs) and their non-catalytic homologs (cSPHs). Here, we define partially the structural hierarchy of An. gambiae cSPs of the CLIPB family, central players in melanization, and characterize their relative contributions to bacterial melanization and to mosquito susceptibility to bacterial infections. Using in vivo genetic analysis we show that the protease cascade branches downstream of the cSPs CLIPB4 and CLIPB17 into two branches one converging on CLIPB10 and the second on CLIPB8. We also show that the contribution of key cSPHs to melanization in vivo in response to diverse microbial challenges is more significant than any of the individual cSPs, possibly due to partial functional redundancy among the latter. Interestingly, we show that the key cSPH CLIPA8 which is essential for the efficient activation cleavage of CLIPBs in vivo is efficiently cleaved itself by several CLIPBs in vitro, suggesting that cSPs and cSPHs regulate signal amplification and propagation in melanization cascades by providing positive reinforcement upstream and downstream of each other.

Serine protease homolog pairs CLIPA4-A6, A4-A7Δ, and A4-A12 act as cofactors for proteolytic activation of prophenoloxidase-2 and -7 in Anopheles gambiae

Phenoloxidase (PO) catalyzed melanization and other insect immune responses are mediated by serine proteases (SPs) and their noncatalytic homologs (SPHs). Many of these SP-like proteins have a regulatory clip domain and are called CLIPs. In most insects studied so far, PO precursors are activated by a PAP (i.e., PPO activating protease) and its cofactor of clip-domain SPHs. Although melanotic encapsulation is a well-known refractory mechanism of mosquitoes against malaria parasites, it is unclear if a cofactor is required for PPO activation. In Anopheles gambiae, CLIPA4 is 1:1 orthologous to Manduca sexta SPH2; CLIPs A5-7, A12-14, A26, A31, A32, E6, and E7 are 11:4 orthologous to M. sexta SPH1a, 1b, 4, and 101, SPH2 partners in the cofactors. Here we produced proCLIPs A4, A6, A7Δ, A12, and activated them with CLIPB9 or M. sexta PAP3. A. gambiae PPO2 and PPO7 were expressed in Escherichia coli for use as PAP substrates. CLIPB9 was mutated to CLIPB9Xa by including a Factor Xa cleavage site. CLIPA7Δ was a deletion mutant with a low complexity region removed. After PAP3 or CLIPB9Xa processing, CLIPA4 formed a high Mr complex with CLIPA6, A7Δ or A12, which assisted PPO2 and PPO7 activation. High levels of specific PO activity (55-85 U/μg for PO2 and 1131-1630 U/μg for PO7) were detected in vitro, indicating that cofactor-assisted PPO activation also occurs in this species. The cleavage sites and mechanisms for complex formation and cofactor function are like those reported in M. sexta and Drosophila melanogaster. In conclusion, these data suggest that the three (and perhaps more) SPHI-II pairs may form cofactors for CLIPB9-mediated activation of PPOs for melanotic encapsulation in A. gambiae.

Insight into the structural hierarchy of the protease cascade that regulates the mosquito melanization response

Serine protease cascades regulate important insect immune responses, including melanization and Toll pathway activation. In the context of melanization, central components of these cascades are clip domain serine proteases (CLIPs) including the catalytic, clip domain serine proteases (cSPs) and their non-catalytic homologs (cSPHs). Here, we define partially the structural hierarchy of An. gambiae cSPs of the CLIPB family, central players in melanization, and characterize their relative contributions to bacterial melanization and to mosquito susceptibility to bacterial infections. Using in vivo genetic analysis we show that the protease cascade branches downstream of the cSPs CLIPB4 and CLIPB17 into two branches one converging on CLIPB10 and the second on CLIPB8. We also show that the contribution of key cSPHs to melanization in vivo in response to diverse microbial challenges is more significant than any of the individual cSPs, possibly due to partial functional redundancy among the latter. Interestingly, we show that the key cSPH CLIPA8 which is essential for the efficient activation cleavage of CLIPBs in vivo is efficiently cleaved itself by several CLIPBs in vitro, suggesting that cSPs and cSPHs regulate signal amplification and propagation in melanization cascades by providing positive reinforcement upstream and downstream of each other.

Two clip-domain serine protease homologs, cSPH35 and cSPH242, act as a cofactor for prophenoloxidase-1 activation in Drosophila melanogaster

Insect phenoloxidases (POs) catalyze phenol oxygenation and o-diphenol oxidation to form reactive intermediates that kill invading pathogens and form melanin polymers. To reduce their toxicity to host cells, POs are produced as prophenoloxidases (PPOs) and activated by a serine protease cascade as required. In most insects studied so far, PPO activating proteases (PAPs) generate active POs in the presence of a high Mr cofactor, comprising two serine protease homologs (SPHs) each with a Gly residue replacing the catalytic Ser of an S1A serine protease (SP). These SPHs have a regulatory clip domain at the N-terminus, like most of the SP cascade members including PAPs. In Drosophila, PPO activation and PO-catalyzed melanization have been examined in genetic analyses but it is unclear if a cofactor is required for PPO activation. In this study, we produced the recombinant cSPH35 and cSPH242 precursors, activated them with Manduca sexta PAP3, and confirmed their predicted role as a cofactor for Drosophila PPO1 activation by MP2 (i.e., Sp7). The cleavage sites and mechanisms for complex formation and cofactor function are highly similar to those reported in M. sexta. In the presence of high Mr complexes of the cSPHs, PO at a high specific activity of 260 U/μg was generated in vitro. To complement the in vitro analysis, we measured hemolymph PO activity levels in wild-type flies, cSPH35, and cSPH242 RNAi lines. Compared with the wild-type flies, only 4.4% and 18% of the control PO level (26 U/μl) was detected in the cSPH35 and cSPH242 knockdowns, respectively. Consistently, percentages of adults with a melanin spot at the site of septic pricking were 82% in wild-type, 30% in cSPH35 RNAi, and 53% in cSPH242 RNAi lines; the survival rate of the control (45%) was significantly higher than those (30% and 15%) of the two RNAi lines. These data suggest that Drosophila cSPH35 and cSPH242 are components of a cofactor for MP2-mediated PPO1 activation, which are indispensable for early melanization in adults.

CLIPB4 Is a Central Node in the Protease Network that Regulates Humoral Immunity in Anopheles gambiae Mosquitoes

Insect humoral immune responses are regulated in part by protease cascades, whose components circulate as zymogens in the hemolymph. In mosquitoes, these cascades consist of clip-domain serine proteases (cSPs) and/or their non-catalytic homologs, which form a complex network, whose molecular make-up is not fully understood. Using a systems biology approach, based on a co-expression network of gene family members that function in melanization and co-immunoprecipitation using the serine protease inhibitor (SRPN)2, a key negative regulator of the melanization response in mosquitoes, we identify the cSP CLIPB4 from the African malaria mosquito Anopheles gambiae as a central node in this protease network. CLIPB4 is tightly co-expressed with SRPN2 and forms protein complexes with SRPN2 in the hemolymph of immune-challenged female mosquitoes. Genetic and biochemical approaches validate our network analysis and show that CLIPB4 is required for melanization and antibacterial immunity, acting as a prophenoloxidase (proPO)-activating protease, which is inhibited by SRPN2. In addition, we provide novel insight into the structural organization of the cSP network in An. gambiae, by demonstrating that CLIPB4 is able to activate proCLIPB8, a cSP upstream of the proPO-activating protease CLIPB9. These data provide the first evidence that, in mosquitoes, cSPs provide branching points in immune protease networks and deliver positive reinforcement in proPO activation cascades.

CLIPB4 is a central node in the protease network that regulates humoral immunity in Anopheles gambiae mosquitoes

Insect humoral immune responses are regulated in part by protease cascades, whose components circulate as zymogens in the hemolymph. In mosquitoes, these cascades consist of clip domain serine proteases (cSPs) and/or their non-catalytic homologs (cSPHs), which form a complex network, whose molecular make-up is not fully understood. Using a systems biology approach, based on a co-expression network of gene family members that function in melanization and co-immunoprecipitation using the serine protease inhibitor (SRPN)2, a key negative regulator of the melanization response in mosquitoes, we identify the cSP CLIPB4 from the African malaria mosquito Anopheles gambiae as a central node in this protease network. CLIPB4 is tightly co-expressed with SRPN2 and forms protein complexes with SRPN2 in the hemolymph of immune-challenged female mosquitoes. Genetic and biochemical approaches validate our network analysis and show that CLIPB4 is required for melanization and antibacterial immunity, acting as a prophenoloxidase (proPO)-activating protease, which is inhibited by SRPN2. In addition, we provide novel insight into the structural organization of the cSP network in An. gambiae, by demonstrating that CLIPB4 is able to activate proCLIPB8, a cSP upstream of the proPO-activating protease CLIPB9. These data provide the first evidence that, in mosquitoes, cSPs provide branching points in immune protease networks and deliver positive reinforcement in proPO activation cascades.

Gene coexpression network during ontogeny in the yellow fever mosquito, Aedes aegypti

BACKGROUND

The behaviors and ontogeny of Aedes aegypti are closely related to the spread of diseases caused by dengue (DENV), chikungunya (CHIKV), Zika (ZIKV), and yellow fever (YFV) viruses. During the life cycle, Ae. aegypti undergoes drastic morphological, metabolic, and functional changes triggered by gene regulation and other molecular mechanisms. Some essential regulatory factors that regulate insect ontogeny have been revealed in other species, but their roles are still poorly investigated in the mosquito.

RESULTS

Our study identified 6 gene modules and their intramodular hub genes that were highly associated with the ontogeny of Ae. aegypti in the constructed network. Those modules were found to be enriched in functional roles related to cuticle development, ATP generation, digestion, immunity, pupation control, lectins, and spermatogenesis. Additionally, digestion-related pathways were activated in the larvae and adult females but suppressed in the pupae. The integrated protein‒protein network also identified cilium-related genes. In addition, we verified that the 6 intramodular hub genes encoding proteins such as EcKinase regulating larval molt were only expressed in the larval stage. Quantitative RT‒PCR of the intramodular hub genes gave similar results as the RNA-Seq expression profile, and most hub genes were ontogeny-specifically expressed.

CONCLUSIONS

The constructed gene coexpression network provides a useful resource for network-based data mining to identify candidate genes for functional studies. Ultimately, these findings will be key in identifying potential molecular targets for disease control.

Characterization and transcriptomic analysis of a native fungal pathogen against the rice pest Nilaparvata lugens

The brown planthopper (BPH), Nilaparvata lugens, is one of the most destructive pests of rice. Given the threats posed by insecticide resistance to its control, eco-friendly strategies based on microbial pathogens emerged as a promising biocontrol alternative. In the present study, we isolated a native fungal pathogen against BPH from infected BPH cadavers and preliminarily identified as a strain of Aspergillus fumigatus based on morphological and molecular methods. Laboratory bioassay revealed that this fungal strain was highly virulent to BPH both at nymphal and adult stages, with the median lethal times (LT₅₀) of 7.5 and 5.8 days under high conidial concentration of 1 × 10⁹ conidia mL^-1. A genome-wide view of gene expressions in BPH against fungal attack was analyzed by transcriptomic sequencing and consequently a large number of differentially expressed genes that mainly involved in host immune defense and cell detoxification were found. RNAi-mediated knockdown of an upregulated gene encoding a serine protease (NlSPN) could cause a significant decrease in BPH survival. Combination of dsRNA injection and fungal infection showed an additive effect on BPH mortality, which provided clues to develop new pest management strategies against BPH.

Identification and gene expression analysis of serine proteases and their homologs in the Asian corn borer Ostrinia furnacalis

Serine proteases (SPs) and their homologs (SPHs) are among the best-characterized gene families. They are involved in several physiological processes, including digestion, embryonic development and immunity. In the current study, a total of 177 SPs-related genes were characterized in the genome of Ostrinia furnacalis. The activation site of SPs/SPHs and enzyme specificity of SPs were identified, and the findings showed that most of the SPs analyzed possessed trypsin substrate specificity. Several SPs/SPHs with similar simple gene structures had tandem repeat-like distributions on the scaffold, indicated that gene expansion has occurred in this large family. Furthermore, we constructed 30 RNA sequencing libraries including four with developmental stage and four middle larval stage tissues to study the transcript levels of these genes. Differentially upregulated and downregulated genes were obtained via data analysis. More than one-quarter of the genes were specifically identified as highly expressed in the midgut in compared to the other three tissues evaluated. In the current study, the domain structure, gene location and phylogenetic relationship of genes in O. furnacalis were explored. Orthologous comparisons of SPs/SPHs between model insects and O. furnacalis indicated their possible functions. This information provides a basis for understanding the functional roles of this large family.

Immune gene expression in the mosquito vector Culex quinquefasciatus during an avian malaria infection

Plasmodium relictum is the most widespread avian malaria parasite in the world. It is listed as one of the 100 most dangerous invasive species, having been responsible for the extinction of several endemic bird species, and the near-demise of several others. Here we present the first transcriptomic study focused on the effect of P. relictum on the immune system of its vector (the mosquito Culex quinquefasciatus) at different times post-infection. We show that over 50% of immune genes identified as being part of the Toll pathway and 30%-40% of the immune genes identified within the Imd pathway are overexpressed during the critical period spanning the parasite's oocyst and sporozoite formation (8-12 days), revealing the crucial role played by both these pathways in this natural mosquito-Plasmodium combination. Comparison of infected mosquitoes with their uninfected counterparts also revealed some unexpected immune RNA expression patterns earlier and later in the infection: significant differences in expression of several immune effectors were observed as early as 30 min after ingestion of the infected blood meal. In addition, in the later stages of the infection (towards the end of the mosquito lifespan), we observed an unexpected increase in immune investment in uninfected, but not in infected, mosquitoes. In conclusion, our work extends the comparative transcriptomic analyses of malaria-infected mosquitoes beyond human and rodent parasites and provides insights into the degree of conservation of immune pathways and into the selective pressures exerted by Plasmodium parasites on their vectors.

Molting-related proteases in the brown planthopper, Nilaparvata lugens

Digestion and absorption of old cuticles during insect molting are necessary for new cuticle formation, during which complicated enzyme catalysis is essential. To date, a few carboxypeptidases, aminopeptidases and serine proteases (mostly trypsins) connected with cuticle digestion, zymogen activation and histological differentiation during the ecdysis of lepidopteran, dipteran and hymenopteran insects have been identified. However, little is known about these proteins in hemimetabolous insects. In this study, we identified 33 candidate trypsin and trypsin-like homologs, 14 metallocarboxypeptidase and 32 aminopeptidase genes in the brown planthopper Nilaparvata lugens, a hemipteran rice pest. Among the proteins encoded by these genes, 9 trypsin-like proteases, 3 metallocarboxypeptidases and 1 aminopeptidase were selected as potential procuticle hydrolases by bioinformatics analysis and in vivo validation. RNA interference targeting these genes demonstrated that 3 trypsin-like proteases (NlTrypsin-8, NlTrypsin-29 and NlTrypsin-32) genes and 1 metallocarboxypeptidase (NlCpB) gene were found to be essential for ecdysis in N. lugens; specifically, gene silencing led to incomplete cuticle degradation and arrested ecdysis, causing lethal morphological phenotype acquisition. Spatiotemporal expression profiling by quantitative PCR and western blotting revealed their specific expression in the integument and their periodic expression during each stadium, with a peak before ecdysis and eclosion. Transmission electron microscopy demonstrated corresponding ultrastructural defects after RNAi targeting, with NlCpB-silenced specimens having the most undigested old procuticles. Immunohistochemical staining revealed that NlTrypsin-8, NlTrypsin-29 and NlCpB were predominantly located in the exuvial space. This research further adds to our understanding of proteases and its potential role in insect ecdysis.

Population genomics reveal distinct and diverging populations of An. minimus in Cambodia

Anopheles minimus is an important malaria vector throughout its wide geographic range across Southeast Asia. Genome sequencing could provide important insights into the unique malaria transmission dynamics in this region, where many vector species feed and rest outdoors. We describe results from a study using Illumina deep whole-genome sequencing of 302 wild-caught An. minimus collected from three Cambodian provinces over several years (2010, 2014, 2016) and seasons to examine the level of population structure and genetic diversity within this species. These specimens cluster into four distinct populations of An. minimus s.s., with two populations overlapping geographically. We describe the underlying genetic diversity and divergence of these populations and investigated the genetic variation in genes known to be involved in insecticide resistance. We found strong signals of selection within these An. minimus populations, most of which were present in the two Northeastern Cambodian populations and differ from those previously described in African malaria vectors. Cambodia is the focus of the emergence and spread of drug-resistant malaria parasites, so understanding the underlying genetic diversity and resilience of the vectors of these parasites is key to implementing effective malaria control and elimination strategies. These data are publicly available as part of the MalariaGEN Vector Observatory, an open access resource of genome sequence data.

CLIPA7 Exhibits Pleiotropic Roles in the Anopheles gambiae Immune Response

Clip domain serine proteases and clip domain serine protease homologs (cSPHs) are key components of serine protease cascades that drive the melanization response. Despite lacking catalytic activity, cSPHs play essential roles in regulating melanization, but the spectrum of functions they catalyze within and outside these cascades is not fully understood. Aside from their classical role as cofactors for PPO activation, we have previously revealed an unprecedented complexity in the function and molecular organization of these cSPHs in the immune response of the malaria vector Anopheles gambiae. Here, we add yet another dimension to the complex roles underpinning the contributions of cSPHs to mosquito immunity by showing that CLIPA7, a member of the expanded cSPH family, defines a novel branch within the cSPH network that is essential for the melanization of Escherichia coli but not Plasmodium ookinetes or Gram-positive bacteria. Despite its dispensability for the melanization of Gram-positive bacteria, we show that CLIPA7 is required for the clearance of systemic infections with Staphylococcus aureus. CLIPA7 is produced by hemocytes and associates with the surfaces of live E. coli and S. aureus cells in vivo as well as with those of melanized cells. Based on its RNAi phenotypes and its unique domain architecture among A. gambiae cSPHs including the presence of an RGD motif, we propose that CLIPA7 exhibits pleiotropic roles in mosquito immunity that extend beyond the regulation of melanization to microbial clearance.

Identification and expression profiling of serine protease-related genes in Tenebrio molitor

In insects, serine proteases and serine protease homologs (SPs/SPHs) are involved in a variety of physiological processes including digestion, development, and immunity. Here, we identified 112 SP and 88 SPH genes in the genome of the yellow mealworm, Tenebrio molitor. Based on the features of domain structure, they were divided into "S" group containing single Tryp-SPc or Tryp-SPHc domain, "C" group containing 1-4 CLIP domain (CLIPA-D) and "M" group containing the CBD, CUB, EGF, Fz, Gd, LDLa, PAN, SEA, SR, Sushi, and TSP domains, and have 115, 48, and 37 gene members, respectively. According to the active sites in the catalytic triad, the putative trypsin, chymotrypsin, or elastase-like enzyme specificity of the identified SPs/SPHs were predicted. Phylogenetic and genomic location analyses revealed that gene duplication exists in the large amount of SPs/SPHs. Gene expression profiling using RNA-seq data along with real time reverse transcription-polymerase chain reaction analysis showed that most SP/SPH genes display life stage specific expression patterns, indicating their important roles in development. Many SP/SPH genes are specifically or highly expressed in the gut, salivary gland, fat body, hemocyte, ovary, and testis, suggesting that they participate in digestion, immunity, and reproduction. The findings lay the foundation for further functional characterization of SPs/SPHs in T. molitor.

Transcriptome Analysis Reveals Early Hemocyte Responses upon In Vivo Stimulation with LPS in the Stick Insect Bacillus rossius (Rossi, 1788)

Simple Summary

Non-model insect species such as B. rossius suffer from a profound gap of knowledge regarding the temporal progression of physiological responses following the challenge with bacterial pathogens or cell wall components thereof. The reason for this mostly lies in the lack of genomic/transcriptomic resources, which would provide an unparalleled in-depth capacity in the analysis of molecular, biochemical, and metabolic mechanisms. We present a high-quality transcriptome obtained from high-coverage sequencing of hemocytes harvested from adult stick insect specimens both pre- and post-LPS stimulation. Such a resource served as the basis for a stringent differential gene expression and functional enrichment analyses, the results of which were characterized and discussed in depth. Selected transcripts encoding for C-type lectins and ML-domain containing proteins were further investigated from a phylogenetic perspective. Overall, these findings shed light on the physiological responses driven by a short-term LPS stimulation in the European stick insect.

Abstract

Despite a growing number of non-model insect species is being investigated in recent years, a greater understanding of their physiology is prevented by the lack of genomic resources. This is the case of the common European stick insect Bacillus rossius (Rossi, 1788): in this species, some knowledge is available on hemocyte-related defenses, but little is known about the physiological changes occurring in response to natural or experimental challenges. Here, the transcriptional signatures of adult B. rossius hemocytes were investigated after a short-term (2 h) LPS stimulation in vivo: a total of 2191 differentially expressed genes, mostly involved in proteolysis and carbohydrate and lipid metabolic processes, were identified in the de novo assembled transcriptome and in-depth discussed. Overall, the significant modulation of immune signals—such as C-type lectins, ML domain-containing proteins, serpins, as well as Toll signaling-related molecules—provide novel information on the early progression of LPS-induced responses in B. rossius.

Transcriptomic and small RNA response to Mayaro virus infection in Anopheles stephensi mosquitoes

Mayaro virus (MAYV) is an arboviral pathogen in the genus Alphavirus that is circulating in South America with potential to spread to naïve regions. MAYV is also one of the few viruses with the ability to be transmitted by mosquitoes in the genus Anopheles, as well as the typical arboviral transmitting mosquitoes in the genus Aedes. Few studies have investigated the infection response of Anopheles mosquitoes. In this study we detail the transcriptomic and small RNA responses of An. stephensi to infection with MAYV via infectious bloodmeal at 2, 7, and 14 days post infection (dpi). 487 unique transcripts were significantly regulated, 78 putative novel miRNAs were identified, and an siRNA response is observed targeting the MAYV genome. Gene ontology analysis of transcripts regulated at each timepoint shows a number of proteases regulated at 2 and 7 dpi, potentially representative of Toll or melanization pathway activation, and repression of pathways related to autophagy and apoptosis at 14 dpi. These findings provide a basic understanding of the infection response of An. stephensi to MAYV and help to identify host factors which might be useful to target to inhibit viral replication in Anopheles mosquitoes.

Mayaro virus (MAYV) is a mosquito-borne Alphavirus responsible for outbreaks in South America and the Caribbean. In this study we infected Anopheles stephensi with MAYV and sequenced mRNA and small RNA to understand how MAYV infection impacts gene transcription and the expression of small RNAs in the mosquito vector. Genes involved with innate immunity and signaling pathways related to cell death are regulated in response to MAYV infection of An. stephensi, we also discovered novel miRNAs and describe the expression patterns of miRNAs, siRNAs, and piRNAs following bloodmeal ingestion. These results suggest that MAYV does induce a molecular response to infection in its mosquito vector species.

Host-Parasite Interactions: Regulation of Leishmania Infection in Sand Fly

PURPOSE

Sand flies are the only proven vectors of leishmaniases, a tropical neglected disease endemic in at least 92 countries. Vector-parasite interactions play a significant role in vector-borne disease transmission. There are various bottlenecks to Leishmania colonization of the sand fly midgut. Such bottlenecks include the production of innate immune-related molecules, digestive proteases, parasite impermeable peritrophic membrane, and resident gut microbiota. These barriers determine the parasite load transmitted and, consequently, the disease outcome in mammalian host. Therefore, it is important to understand the molecular responses of both sand fly and Leishmania during infection.

METHOD

Here, we reviewed the published literature on sand fly-Leishmania interactions bringing together earlier and current findings to highlight new developments and research gaps in the field.

CONCLUSION

Recent research studies on sand fly-Leishmania interaction have revealed contrasting observations to past studies. However, how Leishmania parasites evade the sand fly immune response still needs further research. Sand fly response to Leishmania infection can be best understood by analyzing its tissue transcriptome. Better characterization of the role of midgut components could be a game changer in development of transmission-blocking strategies for leishmaniasis.

A global [Formula: see text] gene co-expression network constructed from hundreds of experimental conditions with missing values

BACKGROUND

Gene co-expression networks (GCNs) can be used to determine gene regulation and attribute gene function to biological processes. Different high throughput technologies, including one and two-channel microarrays and RNA-sequencing, allow evaluating thousands of gene expression data simultaneously, but these methodologies provide results that cannot be directly compared. Thus, it is complex to analyze co-expression relations between genes, especially when there are missing values arising for experimental reasons. Networks are a helpful tool for studying gene co-expression, where nodes represent genes and edges represent co-expression of pairs of genes.

RESULTS

In this paper, we establish a method for constructing a gene co-expression network for the Anopheles gambiae transcriptome from 257 unique studies obtained with different methodologies and experimental designs. We introduce the sliding threshold approach to select node pairs with high Pearson correlation coefficients. The resulting network, which we name AgGCN1.0, is robust to random removal of conditions and has similar characteristics to small-world and scale-free networks. Analysis of network sub-graphs revealed that the core is largely comprised of genes that encode components of the mitochondrial respiratory chain and the ribosome, while different communities are enriched for genes involved in distinct biological processes.

CONCLUSION

Analysis of the network reveals that both the architecture of the core sub-network and the network communities are based on gene function, supporting the power of the proposed method for GCN construction. Application of network science methodology reveals that the overall network structure is driven to maximize the integration of essential cellular functions, possibly allowing the flexibility to add novel functions.

Cleavage activation and functional comparison of Manduca sexta serine protease homologs SPH1a, SPH1b, SPH4, and SPH101 in conjunction with SPH2

Phenoloxidase (PO) is a crucial component of the insect immune response against microbial infection. In the tobacco hornworm Manduca sexta, PO is generated from its precursor proPO by prophenoloxidase activating proteases (PAPs) in the presence of two noncatalytic serine protease homologs (SPHs). cDNA cloning and genome analysis indicate that SPH1a (formerly known as SPH1), SPH1b, SPH4, SPH101, and SPH2 contain a clip domain, a linker, and a protease-like domain (PLD). The first 22 residues of the SPH1b, SPH4, and SPH101 PLDs are identical, and differ from SPH1a only at position 4, Thr¹⁵⁴ substituted with Asn¹⁵⁴ in SPH1a. While the sequence from Edman degradation was used to establish PAP cofactor as a high M_r complex of SPH1a and SPH2, this assignment needed further validation, especially because SPH1b mRNA levels are much higher than SPH1a's and better correlate with SPH2 transcription. Thus, here we determined expression profiles of these SPH genes in different tissues from various developmental stages using highly specific primers. High levels of SPH1b and SPH2 proteins, low SPH4, and no SPH1a or SPH101 were detected in hemolymph from larvae in the feeding, wandering and bar stages, pupae, and adults by targeted LC-MS/MS analysis, based on unique peptides from the trypsin-treated SPHs. We expressed the five proSPHs in baculovirus-infected Sf9 cells for use as standards to identify and quantify their counterparts in plasma samples. Moreover, we tested their cleavage by PAP3 and efficacy of the SPH1a, 1b, 4, and 101 as SPH2 partners in PAP3-mediated proPO activation. PAP3 processed proSPH1b and 101 more readily than proSPH1a and 4; PAP3 activated proPO more efficiently in the presence of SPH2 with SPH101 or SPH1b than with SPH1a or SPH4. These results generally agree with their order of appearance or sequence similarity: SPH101 > SPH1b (98%) > SPH1a (90%) > SPH4 (83%). In other words, likely due to positive selection, products of the newly duplicated genes (SPH1b and SPH101) are more favorable substrates of PAP3 and better SPH2 partners in forming a high M_r cofactor than SPH1a or SPH4 is. Electrophoresis on native gel and immunoblot analysis further indicated that SPH101 or 1b form high M_r complexes more readily than SPH1a or 4 does. In comparison, SPH2 showed a small mobility decrease and then increase on native gel after PAP3 cleavage at the first site. Since the natural cofactor in bar-stage hemolymph is complexes of SPH1 and 2 with an average M_r of 790 kDa, PAP3-activated SPH2 may associate with the higher M_r SPH1b scaffolds to form super-complexes. Their structures and formation in relation to cleavage of SPH1b at different sites await further exploration.

Homeodomain proteins POU-M2, antennapedia and abdominal-B are involved in regulation of the segment-specific expression of the clip-domain serine protease gene CLIP13 in the silkworm, Bombyx mori

Clip-domain serine proteases (CLIPs) play important roles in insect innate immunity and development. Our previous studies indicated that CLIP13, an epidermis-specific gene, was involved in cuticle remodeling during molting and metamorphosis in the silkworm, Bombyx mori. However, the transcriptional regulatory mechanism and regulatory pathways of CLIP13 remained unclear. In the present study, we investigated CLIP13 expression and the regulation pathway controlled by 20-hydroxyecdysone (20E) in the silkworm. At the transcriptional level, expression of CLIP13 exhibited pronounced spatial and temporal specificity in different regions of the epidermis; homeodomain transcription factors POU-M2, antennapedia (Antp), and abdominal-B (Abd-B) showed similar expression change trends as CLIP13 in the head capsule, thorax, and abdomen, respectively. Furthermore, results of cell transfection assays, electrophoretic mobility shift assays, and chromatin immunoprecipitation demonstrated that POU-M2, Antp, and Abd-B were involved in the transcriptional regulation of CLIP13 by directly binding to their cis-response elements in CLIP13 promoter. RNA interference-mediated silencing of POU-M2, Antp, and Abd-B led to a decrease of CLIP13 expression in the head capsule, the epidermis of the 1st to 3rd thoracic segments and the 7th to 10th abdominal segments, respectively. Consistent with CLIP13, 20E treatment significantly upregulated expression of POU-M2, Antp, and Abd-B in the silkworm epidermis. Taken together, these data suggest that 20E positively regulates transcription of CLIP13 via homeodomain proteins POU-M2, Antp, and Abd-B in different regions of the silkworm epidermis during metamorphosis, thus affecting the molting process. Our findings provide new insight into the functions of homeodomain transcription factors in insect molting.

Serine pseudoproteases in physiology and disease

Serine proteases (SPs) constitute a very important family of enzymes, both physiologically and pathologically. The effects produced by these proteins have been explained by their proteolytic activity. However, the discovery of pharmacologically active SP molecules that show no enzymatic activity, as the so-called pseudo SPs or SP homologs (SPHs), has exposed a profoundly neglected possibility of nonenzymatic functions of these SP molecules. In this review, the most thoroughly described SPHs are presented. The main physiological domains in which SPHs operate appear to be in reproduction, embryonic development, immune response, host defense, and hemostasis. Hitherto unexplained actions of SPs should therefore be considered also as the result of the ligand-like attributes of SPs. The gain of a novel function by an SPH is a consequence of specific amino acid replacements that have resulted in a novel interaction interface or a 'catalytic trap'. Unraveling the SP/SPH interactome will provide a description of previously unknown physiological functions of SPs/SPHs, aiding the creation of innovative medical approaches.

Transposable element variants and their potential adaptive impact in urban populations of the malaria vector Anopheles coluzzii

Anopheles coluzzii is one of the primary vectors of human malaria in sub-Saharan Africa. Recently, it has spread into the main cities of Central Africa threatening vector control programs. The adaptation of An. coluzzii to urban environments partly results from an increased tolerance to organic pollution and insecticides. Some of the molecular mechanisms for ecological adaptation are known, but the role of transposable elements (TEs) in the adaptive processes of this species has not been studied yet. As a first step toward assessing the role of TEs in rapid urban adaptation, we sequenced using long reads six An. coluzzii genomes from natural breeding sites in two major Central Africa cities. We de novo annotated TEs in these genomes and in an additional high-quality An. coluzzii genome, and we identified 64 new TE families. TEs were nonrandomly distributed throughout the genome with significant differences in the number of insertions of several superfamilies across the studied genomes. We identified seven putatively active families with insertions near genes with functions related to vectorial capacity, and several TEs that may provide promoter and transcription factor binding sites to insecticide resistance and immune-related genes. Overall, the analysis of multiple high-quality genomes allowed us to generate the most comprehensive TE annotation in this species to date and identify several TE insertions that could potentially impact both genome architecture and the regulation of functionally relevant genes. These results provide a basis for future studies of the impact of TEs on the biology of An. coluzzii.

Trained Immunity in Anopheles gambiae: Antibacterial Immunity Is Enhanced by Priming via Sugar Meal Supplemented With a Single Gut Symbiotic Bacterial Strain

Mosquitoes have evolved an effective innate immune system. The mosquito gut accommodates various microbes, which play a crucial role in shaping the mosquito immune system during evolution. The resident bacteria in the gut microbiota play an essential role in priming basal immunity. In this study, we show that antibacterial immunity in Anopheles gambiae can be enhanced by priming via a sugar meal supplemented with bacteria. Serratia fonticola S1 and Enterobacter sp. Ag1 are gut bacteria in mosquitoes. The intrathoracic injection of the two bacteria can result in an acute hemocoelic infection in the naïve mosquitoes with mortality of ∼40% at 24 h post-infection. However, the Enterobacter orSerratia primed mosquitoes showed a better 24 h survival upon the bacterial challenge. The priming confers the protection with a certain degree of specificity, the Enterobacter primed mosquitoes had a better survival upon the Enterobacter but not Serratia challenge, and the Serratia primed mosquitoes had a better survival upon the Serratia but not Enterobacter challenge. To understand the priming-mediated immune enhancement, the transcriptomes were characterized in the mosquitoes of priming as well as priming plus challenges. The RNA-seq was conducted to profile 10 transcriptomes including three samples of priming conditions (native microbiota, Serratia priming, and Enterobacter priming), six samples of priming plus challenges with the two bacteria, and one sample of injury control. The three priming regimes resulted in distinctive transcriptomic profiles with about 60% of genes affected by both bacteria. Upon challenges, different primed mosquitoes displayed different transcriptomic patterns in response to different bacteria. When a primed cohort was challenged with a heterogenous bacterium, more responsive genes were observed than when challenged with a homogenous bacterium. As expected, many canonical immune genes were responsive to the priming and challenge, but much more non-immune genes with various functions were also responsive in the contexts, which implies that the prior priming triggers a delicately coordinated systemic regulation that results in an enhanced immunity against the subsequent challenge. Besides the participation of typical immune pathways, the transcriptome data suggest the involvement of lysosome and metabolism in the context. Overall, this study demonstrated a trained immunity via priming with bacteria in diet.

Regulators and signalling in insect antimicrobial innate immunity: Functional molecules and cellular pathways

Insects possess an immune system that protects them from attacks by various pathogenic microorganisms that would otherwise threaten their survival. Immune mechanisms may deal directly with the pathogens by eliminating them from the host organism or disarm them by suppressing the synthesis of toxins and virulence factors that promote the invasion and destructive action of the intruder within the host. Insects have been established as outstanding models for studying immune system regulation because innate immunity can be explored as an integrated system at the level of the whole organism. Innate immunity in insects consists of basal immunity that controls the constitutive synthesis of effector molecules such as antimicrobial peptides, and inducible immunity that is activated after detection of a microbe or its product(s). Activation and coordination of innate immune defenses in insects involve evolutionary conserved immune factors. Previous research in insects has led to the identification and characterization of distinct immune signalling pathways that modulate the response to microbial infections. This work has not only advanced the field of insect immunology, but it has also rekindled interest in the innate immune system of mammals. Here we review the current knowledge on key molecular components of insect immunity and discuss the opportunities they present for confronting infectious diseases in humans.

CLIPB10 is a Terminal Protease in the Regulatory Network That Controls Melanization in the African Malaria Mosquito Anopheles gambiae

Humoral immune responses in animals are often tightly controlled by regulated proteolysis. This proteolysis is exerted by extracellular protease cascades, whose activation culminates in the proteolytic cleavage of key immune proteins and enzymes. A model for such immune system regulation is the melanization reaction in insects, where the activation of prophenoxidase (proPO) leads to the rapid formation of eumelanin on the surface of foreign entities such as parasites, bacteria and fungi. ProPO activation is tightly regulated by a network of so-called clip domain serine proteases, their proteolytically inactive homologs, and their serpin inhibitors. In Anopheles gambiae, the major malaria vector in sub-Saharan Africa, manipulation of this protease network affects resistance to a wide range of microorganisms, as well as host survival. However, thus far, our understanding of the molecular make-up and regulation of the protease network in mosquitoes is limited. Here, we report the function of the clip domain serine protease CLIPB10 in this network, using a combination of genetic and biochemical assays. CLIPB10 knockdown partially reversed melanotic tumor formation induced by Serpin 2 silencing in the absence of infection. CLIPB10 was also partially required for the melanization of ookinete stages of the rodent malaria parasite Plasmodium berghei in a refractory mosquito genetic background. Recombinant serpin 2 protein, a key inhibitor of the proPO activation cascade in An. gambiae, formed a SDS-stable protein complex with activated recombinant CLIPB10, and efficiently inhibited CLIPB10 activity in vitro at a stoichiometry of 1.89:1. Recombinant activated CLIPB10 increased PO activity in Manduca sexta hemolymph ex vivo, and directly activated purified M. sexta proPO in vitro. Taken together, these data identify CLIPB10 as the second protease with prophenoloxidase-activating function in An. gambiae, in addition to the previously described CLIPB9, suggesting functional redundancy in the protease network that controls melanization. In addition, our data suggest that tissue melanization and humoral melanization of parasites are at least partially mediated by the same proteases.

Changes in composition and levels of hemolymph proteins during metamorphosis of Manduca sexta

The tobacco hornworm, Manduca sexta, is a lepidopteran model species widely used to study insect biochemical processes. Some of its larval hemolymph proteins are well studied, and a detailed proteomic analysis of larval plasma proteins became available in 2016, revealing features such as correlation with transcriptome data, formation of immune complexes, and constitution of an immune signaling system in hemolymph. It is unclear how the composition of these proteins may change in other developmental stages. In this paper, we report the proteomes of cell-free hemolymph from prepupae, pupae on day 4 and day 13, and young adults. Of the 1824 proteins identified, 907 have a signal peptide and 410 are related to immunity. Drastic changes in abundance of the storage proteins, lipophorins and vitellogenin, for instance, reflect physiological differences among prepupae, pupae, and adults. Considerably more proteins lacking signal peptide are present in the late pupae, suggesting that plasma contains relatively low concentrations of intracellular components released from remodeling tissues during metamorphosis. The defense proteins detected include 43 serine proteases and 11 serine protease homologs. Some of these proteins are members of the extracellular immune signaling network found in feeding larvae, and others may play additional roles and hence confer new features in the later life stages. In summary, the proteins and their levels revealed in this study, together with their transcriptome data, are expected to stimulate focused explorations of humoral immunity and other physiological systems in wandering larvae, pupae, and adults of M. sexta and shed light upon functional and comparative genomic research in other holometabolous insects.

Identification of Immune Response to Sacbrood Virus Infection in Apis cerana Under Natural Condition

Chinese sacbrood virus (CSBV) is a serious threat to eastern honeybees (Apis cerana), especially larvae. However, the pathological mechanism of this deadly disease remains unclear. Here, we employed mRNA and small RNA (sRNA) transcriptome approach to investigate the microRNAs (miRNAs) and small interfering RNAs (siRNAs) expression changes of A. cerana larvae infected with CSBV under natural condition. We found that serine proteases involved in immune response were down-regulated, while the expression of siRNAs targeted to serine proteases were up-regulated. In addition, CSBV infection also affected the expression of larvae cuticle proteins such as larval cuticle proteins A1A and A3A, resulting in increased susceptibility to CSBV infection. Together, our results provide insights into sRNAs that they are likely to be involved in regulating honeybee immune response.

The CLIP-domain serine protease CLIPC9 regulates melanization downstream of SPCLIP1, CLIPA8, and CLIPA28 in the malaria vector Anopheles gambiae

The arthropod melanization immune response is activated by extracellular protease cascades predominantly comprised of CLIP-domain serine proteases (CLIP-SPs) and serine protease homologs (CLIP-SPHs). In the malaria vector, Anopheles gambiae, the CLIP-SPHs SPCLIP1, CLIPA8, and CLIPA28 form the core of a hierarchical cascade downstream of mosquito complement that is required for microbial melanization. However, our understanding of the regulatory relationship of the CLIP-SPH cascade with the catalytic CLIP-SPs driving melanization is incomplete. Here, we report on the development of a novel screen to identify melanization pathway components based on the quantitation of melanotic mosquito excreta, eliminating the need for microdissections or hemolymph enzymatic assays. Using this screen, we identified CLIPC9 and subsequent functional analyses established that this protease is essential for the melanization of both Escherichia coli and the rodent malaria parasite Plasmodium berghei. Mechanistically, septic infection with E. coli promotes CLIPC9 cleavage and both full-length and cleaved CLIPC9 localize to this bacterium in a CLIPA8-dependent manner. The steady state level of CLIPC9 in the hemolymph is regulated by thioester-containing protein 1 (TEP1), suggesting it functions downstream of mosquito complement. In support, CLIPC9 cleavage is inhibited following SPCLIP1, CLIPA8, and CLIPA28 knockdown positioning it downstream of the CLIP-SPH cascade. Moreover, like CLIPA8 and CLIPA28, CLIPC9 processing is negatively regulated by serine protease inhibitor 2 (SRPN2). This report demonstrates how our novel excretion-based approach can be utilized to dissect the complex protease networks regulating mosquito melanization. Collectively, our findings establish that CLIPC9 is required for microbial melanization in An. gambiae and shed light on how the CLIP-SPH cascade regulates this potent immune response.

Analysis of blood-induced Anopheles gambiae midgut proteins and sexual stage Plasmodium falciparum interaction reveals mosquito genes important for malaria transmission

Plasmodium invasion of mosquito midguts is a mandatory step for malaria transmission. The roles of mosquito midgut proteins and parasite interaction during malaria transmission are not clear. This study aims to identify mosquito midgut proteins that interact with and affect P. falciparum invasion. Based on gene expression profiles and protein sequences, 76 mosquito secretory proteins that are highly expressed in midguts and up-regulated by blood meals were chosen for analysis. About 61 candidate genes were successfully cloned from Anopheles gambiae and expressed in insect cells. ELISA analysis showed that 25 of the insect cell-expressed recombinant mosquito proteins interacted with the P. falciparum-infected cell lysates. Indirect immunofluorescence assays confirmed 17 of them interacted with sexual stage parasites significantly stronger than asexual stage parasites. Knockdown assays found that seven candidate genes significantly changed mosquitoes' susceptibility to P. falciparum. Four of them (AGAP006268, AGAP002848, AGAP006972, and AGAP002851) played a protective function against parasite invasion, and the other three (AGAP008138, FREP1, and HPX15) facilitated P. falciparum transmission to mosquitoes. Notably, AGAP008138 is a unique gene that only exists in Anopheline mosquitoes. These gene products are ideal targets to block malaria transmission.

In vitro shared transcriptomic responses of Aedes aegypti to arboviral infections: example of dengue and Rift Valley fever viruses

Background

Arthropod borne virus infections are the cause of severe emerging diseases. Among the diseases due to arboviruses, dengue (DEN) and Rift Valley fever (RVF) are in the top ten in the list of diseases responsible of severe human cases worldwide. Understanding the effects of viral infection on gene expression in competent vectors is a challenge for the development of early diagnostic tools and may enable researchers and policy makers to better anticipate outbreaks in the next future.

Methods

In this study, alterations in gene expression across the entire Aedes aegypti genome during infection with DENV and RVFV were investigated in vitro at two time points of infection, the early phase (24 h) and the late phase (6 days) of infection using the RNA sequencing approach

Results

A total of 10 upregulated genes that share a similar expression profile during infection with both viruses at early and late phases of infection were identified. Family B and D clip-domain serine proteases (CLIP) were clearly overrepresented as well as C-type lectins and transferrin.

Conclusions

Our data highlight the presence of 10 viral genes upregulated in Ae. aegypti during infection. They may also be targeted in the case of the development of broad-spectrum anti-viral diagnostic tools focusing the mosquito vectors rather than the mammalian hosts as they may predict the emergence of outbreaks.

Quantitative Proteomic Analysis of Chikungunya Virus-Infected Aedes aegypti Reveals Proteome Modulations Indicative of Persistent Infection

The mosquito-borne chikungunya virus (CHIKV) poses a threat to human health in tropical countries throughout the world. The molecular interactions of CHIKV with its mosquito vector Aedes aegypti are not fully understood. Following oral acquisition of CHIKV via salinemeals, we analyzed changes in the proteome of Ae. aegypti in 12 h intervals by label-free quantification using a timsTOF Pro mass spectrometer. For each of the seven time points, between 2647 and 3167 proteins were identified among CHIKV-infected and noninfected mosquito samples, and fewer than 6% of those identified proteins were affected by the virus. Functional enrichment analysis revealed that the three pathways, Endocytosis, Oxidative phosphorylation, and Ribosome biogenesis, were enriched during CHIKV infection. On the other hand, three pathways of the cellular RNA machinery and five metabolism related pathways were significantly attenuated in the CHIKV-infected samples. Furthermore, proteins associated with cytoskeleton and vesicular transport, as well as various serine-type endopeptidases and metallo-proteinases, were modulated in the presence of CHIKV. Our study reveals biological pathways and novel proteins interacting with CHIKV in the mosquito. Overall, CHIKV infection caused minor changes to the mosquito proteome demonstrating a high level of adaption between the vector and the virus, essentially coexisting in a nonpathogenic relationship. The mass spectrometry data have been deposited to the MassIVE repository (https://massive.ucsd.edu/ProteoSAFe/dataset.jsp?task=abfd14f7015243c69854731998d55df1) with the data set identifier MSV000085115.

Venom serine proteinase homolog of the ectoparasitoid Scleroderma guani impairs host phenoloxidase cascade

The ant-like bethylid ectoparasitoid Scleroderma guani (Hymenoptera: Bethylidae) envenomates host to suppress immune response. Yet, the roles of its venom in inhibiting melanization of the host hemolymph have not been fully characterized. Here, we demonstrated that S. guani envenomation induced strong inhibition of melanization of the hemolymph from Tenebrio molitor (Coleoptera: Tenebrionidae), permitting the successful development of parasitoid offspring. To reveal venom component associated with such function, a serine proteinase homolog (SguaSPH) rich in the venom of S. guani was characterized. It was found that one of the catalytic triad residues for serine proteinase is absent in the amino acid sequence of SguaSPH. This venom component was abundantly expressed in venom apparatus and adult stages. By enzymatic assays, SguaSPH displayed low trypsin and no chymotrypsin activity, and was able to inhibit phenoloxidase activity in the hemolymph of Ostrinia furnacalis (Lepidoptera: Crambidae). The findings suggest that SguaSPH is essential for interfering with hemolymph melanization of S. guani envenomated host via phenoloxidase cascade disruption.

Identification, characterization, and expression analysis of clip-domain serine protease genes in the silkworm, Bombyx mori

Clip-domain serine proteases (CLIPs), characterized by regulatory module clip domains, constitute an important serine protease family identified in insects and other arthropods. They participate in host immune response and embryonic development in a cascade-activated manner. Here, we present a genome-wide identification and expression analysis of CLIP genes in the silkworm, Bombyx mori. A total of 26 CLIP genes were identified in the silkworm genome. Bioinformatics analysis indicated that these CLIPs clustered into four subfamilies (CLIPA-D), and exhibit a close evolutionary relationship with CLIPs of Manduca sexta. Tissue expression profiling revealed that silkworm CLIP genes are mainly expressed in the integument, head, fat body, and hemocytes. Temporal expression profiles showed that 15 CLIP genes were predominantly expressed during the fifth-instar larval stage, early and later period of the pupal stage, and adult stage, whereas 10 CLIP genes were mainly expressed in the wandering stage and middle to later period of the pupal stage in the integument. Pathogens and 20-hydroxyecdysone (20E) induction analysis indicated that 14 CLIP genes were positively regulated by 20E, 9 were negatively regulated by 20E but positively regulated by pathogens, and 5 were positively regulated by both factors in the integument. Together, these results suggested that silkworm CLIP genes may play multiple functions in integument development, including melanization of new cuticle, molting and immune defense. Our data provide a comprehensive understanding of CLIP genes in the silkworm integument and lays a foundation for further functional studies of CLIP genes in the silkworm.

Inhibition of immune pathway-initiating hemolymph protease-14 by Manduca sexta serpin-12, a conserved mechanism for the regulation of melanization and Toll activation in insects

A network of serine proteases (SPs) and their non-catalytic homologs (SPHs) activates prophenoloxidase (proPO), Toll pathway, and other insect immune responses. However, integration and conservation of the network and its control mechanisms have not yet been fully understood. Here we present evidence that these responses are initiated through a conserved serine protease and negatively regulated by serpins in two species, Manduca sexta and Anopheles gambiae. We have shown that M. sexta serpin-12 reduces the proteolytic activation of HP6, HP8, proPO activating proteases (PAPs), SPHs, and POs in larval hemolymph, and we hypothesized that these effects are due to the inhibition of the immune pathway-initiating protease HP14. To test whether these changes are due to HP14 inhibition, we isolated a covalent complex of HP14 with serpin-12 from plasma using polyclonal antibodies against the HP14 protease domain or against serpin-12, and confirmed formation of the complex by 2D-electrophoresis, immunoblotting, and mass spectrometry. Upon recognition of bacterial peptidoglycans or fungal β-1,3-glucan, the zymogen proHP14 became active HP14, which formed an SDS-stable complex with serpin-12 in vitro. Activation of proHP21 by HP14 was suppressed by serpin-12, consistent with the decrease in steps downstream of HP21, proteolytic activation of proPAP3, proSPH1/2 and proPO in hemolymph. Guided by the results of phylogenetic analysis, we cloned and expressed A. gambiae proSP217 (an ortholog of HP14) and core domains of A. gambiae serpin-11 and -17. The recombinant SP217 zymogen became active during expression, with cleavage between Tyr394 and Ile395. Both MsHP14 and AgSP217 cleaved MsSerpin-12 and AgSRPN11 at Leu*Ser (P1*P1') and formed complexes in vitro. ProPO activation in M. sexta plasma increased after recombinant AgSP217 had been added, indicating that it may function in a similar manner as the endogenous initiating protease HP14. Based on these data, we propose that inhibition of an initiating modular protease by a serpin may be a common mechanism in holometabolous insects to regulate proPO activation and other protease-induced immune responses.

The mosquito melanization response requires hierarchical activation of non-catalytic clip domain serine protease homologs

Serine protease cascades regulate important insect immune responses namely melanization and Toll pathway activation. An important component of these cascades are clip-domain serine protease homologs (cSPHs), which are non-catalytic, but essential for activating the enzyme prophenoloxidase (PPO) in the melanization response during septic infections. The activation of cSPHs requires their proteolytic cleavage, yet factors that control their activation and the complexity of their interactions within these cascades remain unclear. Here, we report the identification of CLIPA28 as a novel immune-related cSPH in the malaria vector Anopheles gambiae. Functional genetic analysis using RNA interference (RNAi) revealed that CLIPA28 is essential for the melanization of Plasmodium berghei parasites in refractory mosquitoes, and for mosquito resistance to fungal infections. We further show, using combined biochemical and genetic approaches, that CLIPA28 is member of a network of at least four cSPHs, whereby members are activated in a hierarchical manner following septic infections. Depletion of the complement-like protein TEP1 abolished the activation of this network after septic infections, whereas, depletion of the serine protease inhibitor 2 (SRPN2) triggered enhanced network activation, even in naïve mosquitoes, culminating in a dramatic reduction in cSPHs hemolymph levels, which paralleled that of PPO. Our data suggest that cSPHs are engaged in complex and multilayered interactions within serine protease cascades that regulate melanization, and identify TEP1 and SRPN2 as two master regulators of the cSPH network.

Melanization is a spectacular immune response of insects that culminates in the deposition of melanin on microbial surfaces leading to their death. Despite more than half a century of investigation of this immune reaction in several insect models, important knowledge gaps, specifically concerning the regulation of this response, remain, especially that its output is tightly regulated by complex genetic and biochemical interactions between members of structurally diverse immune gene families. Here, we identify a new hemolymph protein that is involved in the melanization of malaria parasites and show that it forms, together with other proteins of the same gene family, a network, whereby members are activated in an orderly manner during mosquito septic infections. Furthemore, we show that this network is controlled by two structurally distinct immune proteins which act as master regulators of its activation. Dissecting the complexity of the regulatory mechanisms mediating insect melanization responses, may facilitate fine manipulation of mosquito immunity to reduce vectorial capacity.

Building a platform for predicting functions of serine protease-related proteins in Drosophila melanogaster and other insects

Serine proteases (SPs) and serine protease homologs (SPHs) play essential roles in insect physiological processes including digestion, defense and development. Studies of insect genomes, transcriptomes and proteomes have generated a vast amount of information on these proteins, dwarfing the biological data acquired from a few model species. The large number and high diversity of homologous sequences makes it a challenge to use the limited functional information for making predictions across a broad taxonomic group of insects. In this work, we have extensively updated the framework of knowledge on the SP-related proteins in Drosophila melanogaster by identifying 52 new SPs/SPHs, classifying the 257 proteins into four groups (CLIP, gut, single- and multi-domain SPs/SPHs), and detecting inherent connections among phylogenetic relationships, genomic locations and expression profiles for 99 of the genes. Information on the existence of specific proteins in eggs, larvae, pupae and adults is presented to facilitate future research. More importantly, we have developed an approach to reveal close homologous or orthologous relationships among SPs/SPHs from D. melanogaster, Anopheles gambiae, Apis mellifera, Manduca sexta, and Tribolium castaneum thus inspiring functional studies in these and other holometabolous insects. This approach is useful for tackling similar problems on large and diverse protein families in other groups of organisms.

Manduca sexta hemolymph protease-2 (HP2) activated by HP14 generates prophenoloxidase-activating protease-2 (PAP2) in wandering larvae and pupae

Melanization is a universal defense mechanism of insects against microbial infection. During this response, phenoloxidase (PO) is activated from its precursor by prophenoloxidase activating protease (PAP), the terminal enzyme of a serine protease (SP) cascade. In the tobacco hornworm Manduca sexta, hemolymph protease-14 (HP14) is autoactivated from proHP14 to initiate the protease cascade after host proteins recognize invading pathogens. HP14, HP21, proHP1*, HP6, HP8, PAP1-3, and non-catalytic serine protease homologs (SPH1 and SPH2) constitute a portion of the extracellular SP-SPH system to mediate melanization and other immune responses. Here we report the expression, purification, and functional characterization of M. sexta HP2. The HP2 precursor is synthesized in hemocytes, fat body, integument, nerve and trachea. Its mRNA level is low in fat body of 5th instar larvae before wandering stage; abundance of the protein in hemolymph displays a similar pattern. HP2 exists as an active enzyme in plasma of the wandering larvae and pupae in the absence of an infection. HP14 cleaves proHP2 to yield active HP2. After incubating active HP2 with larval hemolymph, we detected higher levels of PO activity, i.e. an enhancement of proPO activation. HP2 cleaved proPAP2 (but not proPAP3 or proPAP1) to yield active PAP2, responsible for a major increase in IEARpNA hydrolysis. PAP2 activates proPOs in the presence of a cofactor of SPH1 and SPH2. In summary, we have identified a new member of the proPO activation system and reconstituted a pathway of HP14-HP2-PAP2-PO. Since high levels of HP2 mRNA were present in integument and active HP2 in plasma of wandering larvae, HP2 likely plays a role in cuticle melanization during pupation and protects host from microbial infection in a soil environment.

Unraveling the Plasmodium vivax sporozoite transcriptional journey from mosquito vector to human host

Malaria parasites transmitted by mosquito bite are remarkably efficient in establishing human infections. The infection process requires roughly 30 minutes and is highly complex as quiescent sporozoites injected with mosquito saliva must be rapidly activated in the skin, migrate through the body, and infect the liver. This process is poorly understood for Plasmodium vivax due to low infectivity in the in vitro models. To study this skin-to-liver-stage of malaria, we used quantitative bioassays coupled with transcriptomics to evaluate parasite changes linked with mammalian microenvironmental factors. Our in vitro phenotyping and RNA-seq analyses revealed key microenvironmental relationships with distinct biological functions. Most notable, preservation of sporozoite quiescence by exposure to insect-like factors coupled with strategic activation limits untimely activation of invasion-associated genes to dramatically increase hepatocyte invasion rates. We also report the first transcriptomic analysis of the P. vivax sporozoite interaction in salivary glands identifying 118 infection-related differentially-regulated Anopheles dirus genes. These results provide important new insights in malaria parasite biology and identify priority targets for antimalarial therapeutic interventions to block P. vivax infection.

The immune strategies of mosquito Aedes aegypti against microbial infection

Yellow fever mosquito Aedes aegypti transmits many devastating arthropod-borne viruses (arboviruses), such as dengue virus, yellow fever virus, Chikungunya virus, and Zika virus, which cause great concern to human health. Mosquito control is an effective method to block the spread of infectious diseases. Ae. aegypti uses its innate immune system to fight against arboviruses, parasites, and fungi. In this review, we briefly summarize the recent findings in the immune response of Ae. aegypti against arboviral and entomopathogenic infections. This review enriches our understanding of the mosquito immune system and provides evidence to support the development of novel mosquito control strategies.