N-glycosylation influences the catalytic activity of mosquito α-glucosidases associated with susceptibility or refractoriness to Lysinibacillus sphaericus

Culex quinquefasciatus Resistant to the Binary Toxin from Lysinibacillus sphaericus Displays a Consistent Downregulation of Pantetheinase Transcripts

Culex quinquefasciatus resistance to the binary (Bin) toxin, the major larvicidal component from Lysinibacillus sphaericus, is associated with mutations in the cqm1 gene, encoding the Bin-toxin receptor. Downregulation of the cqm1 transcript was found in the transcriptome of larvae resistant to the L. sphaericus IAB59 strain, which produces both the Bin toxin and a second binary toxin, Cry48Aa/Cry49Aa. Here, we investigated the transcription profiles of two other mosquito colonies having Bin resistance only. These confirmed the cqm1 downregulation and identified transcripts encoding the enzyme pantetheinase as the most downregulated mRNAs in both resistant colonies. Further quantification of these transcripts reinforced their strong downregulation in Bin-resistant larvae. Multiple genes were found encoding this enzyme in Cx. quinquefasciatus and a recombinant pantetheinase was then expressed in Escherichia coli and Sf9 cells, with its presence assessed in the midgut brush border membrane of susceptible larvae. The pantetheinase was expressed as a ~70 kDa protein, potentially membrane-bound, which does not seem to be significantly targeted by glycosylation. This is the first pantetheinase characterization in mosquitoes, and its remarkable downregulation might reflect features impacted by co-selection with the Bin-resistant phenotype or potential roles in the Bin-toxin mode of action that deserve to be investigated.

A Comparative Analysis of Bombyx mori (Lepidoptera: Bombycidae) β-fructofuranosidase Homologs Reveals Different Post-Translational Regulations in Glyphodes pyloalis Walker (Lepidoptera: Pyralidae)

Simple Summary

The β-fructofuranosidase (β-FFase) encoding gene BmSuc1 regulates the glycometabolism of silkworm larvae, and it participates in the resistance of mulberry alkaloids. However, there is no molecular or biochemical information available about the mulberry pest Glyphodespyloalis Walker β-FFase homologs. In this paper, we have obtained five β-FFase homologous genes in G. pyloalis and characterized the expression and the localization of GpSUC1a in the midgut. The β-FFase activity in the midgut of G. pyloalis larvae and GpSUC1a were both confirmed, while recombinant GpSUC1a displayed little activity as compared with the higher activity of BmSUC1. Some putative N-glycosylation sites were found in GpSUC1a but none in BmSUC1, while there was more methylation in BmSUC1 than in GpSUC1a. The results indicate that such post-translational modifications (PTMs) are differentially supporting that β-FFase are active in these two mulberry feeding caterpillars, and the activation of GpSUC1a may be controlled by a more complex post-translational regulatory system in G. pyloalis larvae. This is the first report on the characterization of β-FFase genes from G. pyloalis and the first comparison of expression regulation between two mulberry feeding insects B. mori and G. pyloalis. Moreover, this research may provide new ideas for the management of mulberry borers.

Abstract

The silk-spinning and Lepidopteran model insect Bombyx mori (Bombycidae) is a mulberry specialist. The BmSuc1 gene is the first β-fructofuranosidase (β-FFase) encoding gene identified in animals, and β-FFase acts as an essential sucrase for glycometabolism modulation in the silkworm larvae, involved in resistance to mulberry alkaloids. Glyphodes pyloalis Walker (Lepidoptera: Pyralidae) is an important mulberry pest leading to heavy economic loss of sericulture. However, no molecular or biochemical information is available about G. pyloalis β-FFase homologs. In this study, five β-FFase homologous genes in G. pyloalis were obtained. The genes GpSuc1a and GpSuc2c were expressed in the midgut; GpSuc2c encodes a truncated polypeptide. The expression and the localization of GpSUC1a in the midgut was characterized. Whereas recombinant GpSUC1a expressed in both Escherichia coli and BmN cells displayed little activity as compared with higher activity of BmSUC1, β-FFase activity in the larval midgut of G. pyloalis and GpSUC1a purified from the midgut were both confirmed. The data suggested that the activation of GpSUC1a is probably controlled by a more complicated post-translational regulation system in G. pyloalis larvae than that of BmSUC1 in B. mori. To study post-translational modifications (PTMs), GpSUC1a and BmSUC1 were purified from larval midguts using immunoprecipitation and subjected to LC-MS to perform PTMs analysis. Some putative N-glycosylated sites were found in GpSUC1a but none in BmSUC1, while there was more methylation in BmSUC1 than in GpSUC1a, indicating that such PTMs were supporting the differential β-FFases activities in these two mulberry feeding caterpillars.

Bacterial Toxins Active against Mosquitoes: Mode of Action and Resistance

Larvicides based on the bacteria Bacillus thuringiensis svar. israelensis (Bti) and Lysinibacillus sphaericus are effective and environmentally safe compounds for the control of dipteran insects of medical importance. They produce crystals that display specific and potent insecticidal activity against larvae. Bti crystals are composed of multiple protoxins: three from the three-domain Cry type family, which bind to different cell receptors in the midgut, and one cytolytic (Cyt1Aa) protoxin that can insert itself into the cell membrane and act as surrogate receptor of the Cry toxins. Together, those toxins display a complex mode of action that shows a low risk of resistance selection. L. sphaericus crystals contain one major binary toxin that display an outstanding persistence in field conditions, which is superior to Bti. However, the action of the Bin toxin based on its interaction with a single receptor is vulnerable for resistance selection in insects. In this review we present the most recent data on the mode of action and synergism of these toxins, resistance issues, and examples of their use worldwide. Data reported in recent years improved our understanding of the mechanism of action of these toxins, showed that their combined use can enhance their activity and counteract resistance, and reinforced their relevance for mosquito control programs in the future years.

Molecular and biological features of Culex quinquefasciatus homozygous larvae for two cqm1 alleles that confer resistance to Lysinibacillus sphaericus larvicides

BACKGROUND

Culex quinquefasciatus resistance to the binary toxin from Lysinibacillus sphaericus larvicides can occur because of mutations in the cqm1 gene that prevents the expression of the toxin receptor, Cqm1 α-glucosidase. In a resistant laboratory-selected colony maintained for more than 250 generations, cqm1_REC and cqm1_REC-2 resistance alleles were identified. The major allele initially found, cqm1_REC , became minor and was replaced by cqm1_REC-2 . This study aimed to investigate the features associated with homozygous larvae for each allele to understand the reasons for the allele replacement and to generate knowledge on resistance to microbial larvicides.

RESULTS

Homozygous larvae for each allele were compared. Both larvae displayed the same level of resistance to the binary toxin (3500-fold); therefore, a change in phenotype was not the reason for the replacement observed. The lack of Cqm1 expression did not reduce the total specific α-glucosidase activity for homozygous cqm1_REC and cqm1_REC-2 larvae, which were statistically similar to the susceptible strain, using artificial or natural substrates. The expression of eight Cqm1 paralog α-glucosidases was demonstrated in resistant and susceptible larvae. Bioassays in which cqm1_REC or cqm1_REC-2 homozygous larvae were reared under stressful conditions showed that most adults produced were cqm1_REC-2 homozygous (69%). Comparatively, in the offspring of a heterozygous sub-colony reared under optimal conditions for 20 generations, the cqm1_REC allele assumed a higher frequency (0.72).

CONCLUSION

Homozygous larvae for each allele exhibited a similar resistant phenotype. However, they presented specific advantages that might favor their selection and can be used in designing resistance management practices. © 2021 Society of Chemical Industry.

Characterization of Glycoside Hydrolase Families 13 and 31 Reveals Expansion and Diversification of α-Amylase Genes in the Phlebotomine Lutzomyia longipalpis and Modulation of Sandfly Glycosidase Activities by Leishmania Infection

Sugar-rich food sources are essential for sandflies to meet their energy demands, achieving more prolonged survival. The digestion of carbohydrates from food is mainly realized by glycoside hydrolases (GH). To identify genes coding for α-glycosidases and α-amylases belonging to Glycoside Hydrolase Family 13 (GH13) and Glycoside Hydrolase Family 31 (GH31) in Lutzomyia longipalpis, we performed an HMMER search against its genome using known sequences from other dipteran species. The sequences retrieved were classified based on BLASTP best hit, analysis of conserved regions by alignment with sequences of proteins with known structure, and phylogenetic analysis comparing with orthologous proteins from other dipteran species. Using RT-PCR analysis, we evaluated the expression of GH13 and GH31 genes, in the gut and rest of the body of females, in four different conditions: non-fed, sugar-fed, blood-fed, and Leishmania mexicana infected females. L. longipalpis has GH13/31 genes that code for enzymes involved in various aspects of sugar metabolism, as carbohydrate digestion, storage, and mobilization of glycogen reserves, proteins involved in transport, control of N-glycosylation quality, as well as others with a putative function in the regulation of myogenesis. These proteins are representatives of GH13 and GH31 families, and their roles seem to be conserved. Most of the enzymes seem to be active with conserved consense sequences, including the expected catalytic residues. α-amylases also demonstrated the presence of calcium and chloride binding sites. L. longipalpis genome shows an expansion in the α-amylase gene family, with two clusters. In contrast, a retraction in the number of α-glucosidases occurred. The expansion of α-amylases is probably related to the specialization of these proteins for different substrates or inhibitors, which might correlate with the higher diversity of plant foods available in the natural habitat of L. longipalpis. The expression of α-glucosidase genes is higher in blood-fed females, suggesting their role in blood digestion. Besides that, in blood-fed females infected with the parasite Leishmania mexicana, these genes were also modulated. Glycoside Hydrolases from families 13 and 31 are essential for the metabolism of L. longipalpis, and GH13 enzymes seem to be involved in the interaction between sandflies and Leishmania.

Structure-function analysis of silkworm sucrose hydrolase uncovers the mechanism of substrate specificity in GH13 subfamily 17 exo-α-glucosidases

The domestic silkworm Bombyx mori expresses two sucrose-hydrolyzing enzymes, BmSUH and BmSUC1, belonging to glycoside hydrolase family 13 subfamily 17 (GH13_17) and GH32, respectively. BmSUH has little activity on maltooligosaccharides, whereas other insect GH13_17 α-glucosidases are active on sucrose and maltooligosaccharides. Little is currently known about the structural mechanisms and substrate specificity of GH13_17 enzymes. In this study, we examined the crystal structures of BmSUH without ligands; in complexes with substrates, products, and inhibitors; and complexed with its covalent intermediate at 1.60-1.85 Å resolutions. These structures revealed that the conformations of amino acid residues around subsite -1 are notably different at each step of the hydrolytic reaction. Such changes have not been previously reported among GH13 enzymes, including exo- and endo-acting hydrolases, such as α-glucosidases and α-amylases. Amino acid residues at subsite +1 are not conserved in BmSUH and other GH13_17 α-glucosidases, but subsite -1 residues are absolutely conserved. Substitutions in three subsite +1 residues, Gln¹⁹¹, Tyr²⁵¹, and Glu⁴⁴⁰, decreased sucrose hydrolysis and increased maltase activity of BmSUH, indicating that these residues are key for determining its substrate specificity. These results provide detailed insights into structure-function relationships in GH13 enzymes and into the molecular evolution of insect GH13_17 α-glucosidases.

A differential transcriptional profile by Culex quinquefasciatus larvae resistant to Lysinibacillus sphaericus IAB59 highlights genes and pathways associated with the resistance phenotype

Background

The study of the mechanisms by which larvae of the Culex quinquefasciatus mosquito survive exposure to the entomopathogen Lysinibacillus sphaericus has benefited substantially from the generation of laboratory-selected colonies resistant to this bacterium. One such colony, RIAB59, was selected after regular long-term exposure of larvae to the L. sphaericus IAB59 strain. This strain is characterized by its ability to produce the well known Binary (Bin) toxin, and the recently characterized Cry48Aa/Cry49Aa toxin, able to kill Bin-resistant larvae. Resistance to Bin is associated with the depletion of its receptor, Cqm1 α-glucosidase, from the larvae midgut. This study aimed to identify novel molecules and pathways associated with survival of the RIAB59 larvae and the resistance phenotype.

Methods

A transcriptomic approach and bioinformatic tools were used to compare the profiles derived from the midguts of larvae resistant and susceptible to L. sphaericus IAB59.

Results

The RNA-seq profiles identified 1355 differentially expressed genes (DEGs), with 673 down- and 682 upregulated transcripts. One of the most downregulated DEGs was cqm1, which validates the approach. Other strongly downregulated mRNAs encode the enzyme pantetheinase, apolipoprotein D, lipases, heat-shock proteins and a number of lesser known and hypothetical polypeptides. Among the upregulated DEGs, the top most encodes a peroxisomal enzyme involved in lipid metabolism, while others encode enzymes associated with juvenile hormone synthesis, ion channels, DNA binding proteins and defense polypeptides. Further analyses confirmed a strong downregulation of several enzymes involved in lipid catabolism while the assignment of DEGs into metabolic pathways highlighted the upregulation of those related to DNA synthesis and maintenance, confirmed by their clustering into related protein networks. Several other pathways were also identified with mixed profiles of down- and upregulated transcripts. Quantitative RT-PCR confirmed the changes in levels seen for selected mRNAs.

Conclusions

Our transcriptome-wide dataset revealed that the RIAB59 colony, found to be substantially more resistant to Bin than to the Cry48Aa/Cry49Aa toxin, developed a differential expression profile as well as metabolic features co-selected during the long-term adaptation to IAB59 and that are most likely linked to Bin resistance.

Electronic supplementary material

The online version of this article (10.1186/s13071-019-3661-y) contains supplementary material, which is available to authorized users.

Receptor protein of Lysinibacillus sphaericus mosquito-larvicidal toxin displays amylomaltase activity

The activated binary toxin (BinAB) from Lysinibacillus sphaericus binds to surface receptor protein (Cqm1) on the midgut cell membrane and kills Culex quinquefasciatus larvae on internalization. Cqm1 is attached to cells via a glycosyl-phosphatidylinositol (GPI) anchor. It has been classified as a member of glycoside hydrolase family 13 of the CAZy database. Here, we report characterization of the ordered domain (residues 23-560) of Cqm1. Gene expressing Cqm1 of BinAB susceptible mosquito was chemically synthesized and the protein was purified using E. coli expression system. Values for the Michaelis-Menten kinetics parameters towards 4-nitrophenyl α-D-glucopyranoside (α-pNPG) substrate were estimated to be 0.44 mM (Km) and 1.9 s^-1 (kcat). Thin layer chromatography experiments established Cqm1 as α-glucosidase competent to cleave α-1,4-glycosidic bonds of maltose and maltotriose with high glycosyltransferase activity to form glucose-oligomers. The observed hydrolysis and synthesis of glucose-oligomers is consistent with open and accessible active-site in the structural model. The protein also hydrolyses glycogen and sucrose. These activities suggest that Cqm1 may be involved in carbohydrate metabolism in mosquitoes. Further, toxic BinA component does not inhibit α-glucosidase activity of Cqm1, while BinB reduced the activity by nearly 50%. The surface plasmon resonance study reveals strong binding of BinB with Cqm1 (Kd, 9.8 nM). BinA interaction with Cqm1 however, is 1000-fold weaker. Notably the estimated Kd values match well with dissociation constants reported earlier with larvae brush border membrane fractions. The Cqm1 protein forms a stable dimer that is consistent with its apical localization in lipid rafts. Its melting temperature (T_m) as observed by thermofluor-shift assay is 51.5 °C and Ca²⁺ provides structural stability to the protein.