Ostrinia furnacalis serpin-3 regulates melanization cascade by inhibiting a prophenoloxidase-activating protease

Serine protease cascade-mediated prophenolxidase activation is a prominent innate immune response in insect defense against the invading pathogens. Serpins regulate this reaction to avoid excessive activation. However, the function of serpins in most insect species, especially in some non-model agriculture insect pests, is largely unknown. We here cloned a full-length cDNA for a serpin, named as serpin-3, from Asian corn borer, Ostrinia furnacalis (Guenée). The open reading frame of serpin-3 encodes 462-amino acid residue protein with a 19-residue signal peptide. It contains a reactive center loop strikingly similar to the proteolytic activation site in prophenoloxidase. Sequence comparison indicates that O. furnacalis serpin-3 is an apparent ortholog of Manduca sexta serpin-3, a defined negative regulator of melanization reaction. Serpin-3 mRNA and protein levels significantly increase after a bacterial or fungal injection. Recombinant serpin-3 efficiently blocks prophenoloxidase activation in larval plasma in a concentration-dependent manner. It forms SDS-stable complexes with serine protease 13 (SP13), and prevents SP13 from cleaving prophenoloxidase. Injection of recombinant serpin-3 into larvae results in decreased fungi-induced melanin synthesis and reduced the expression of attacin, cecropin, gloverin, and peptidoglycan recognition protein-1 genes in the fat body. Altogether, serpin-3 plays important roles in the regulation of prophenoloxidase activation and antimicrobial peptide production in O. furnacalis larvae.     

Mechanisms of nodule-specific melanization in the hemocoel of the silkworm,Bombyx mori

In the insect immune system, nodules are known to be a product of the cellular response against microorganisms and may be a preferential target for melanization. However, the mechanism of nodule-preferential melanization remains to be explored. In this study, we identified several mechanisms of nodule-preferential melanization by analyzing congregation and the activation of several factors involved in the prophenoloxidase (proPO)-activating system in the silkworm, Bombyx mori. Microorganism-binding assays revealed that B. mori larval plasma have an effective invading microorganism-surveillance network consisting of at least six pattern-recognition receptors (PRRs). We also found that a hemolymph serine proteinase, BmHP14, can bind to Saccharomyces cerevisiae. Pull-down assays showed that PRR C-type lectins form protein complexes with serine proteinase homologs, BmSPH1 and BmSPH2, which leads to the activated forms of BmSPH1 and BmSPH2 being gathered on microorganisms and trapped in nodules. Immunostaining analysis revealed that most factors in the proPO-activating system and some factors in the triggering system for antimicrobial peptide production exist in the granules of hemocytes which can gather in nodules. Western blot analysis showed that factors in the proPO-activating system are congregated in formed nodules by their concentration in plasma and aggregating hemocytes.     

Overview of chitin metabolism enzymes in Manduca sexta: Identification,domain organization,phylogenetic analysis and gene expression

Chitin is one of the most abundant biomaterials in nature. The biosynthesis and degradation of chitin in insects are complex and dynamically regulated to cope with insect growth and development. Chitin metabolism in insects is known to involve numerous enzymes, including chitin synthases (synthesis of chitin), chitin deacetylases (modification of chitin by deacetylation) and chitinases (degradation of chitin by hydrolysis). In this study, we conducted a genome-wide search and analysis of genes encoding these chitin metabolism enzymes in Manduca sexta. Our analysis confirmed that only two chitin synthases are present in M. sexta as in most other arthropods. Eleven chitin deacetylases (encoded by nine genes) were identified, with at least one representative in each of the five phylogenetic groups that have been described for chitin deacetylases to date. Eleven genes encoding for family 18 chitinases (GH18) were found in the M. sexta genome. Based on the presence of conserved sequence motifs in the catalytic sequences and phylogenetic relationships, two of the M. sexta chitinases did not cluster with any of the current eight phylogenetic groups of chitinases: two new groups were created (groups IX and X) and their characteristics are described. The result of the analysis of the Lepidoptera-specific chitinase-h (group h) is consistent with its proposed bacterial origin. By analyzing chitinases from fourteen species that belong to seven different phylogenetic groups, we reveal that the chitinase genes appear to have evolved sequentially in the arthropod lineage to achieve the current high level of diversity observed in M. sexta. Based on the sequence conservation of the catalytic domains and on their developmental stage- and tissue-specific expression, we propose putative functions for each group in each category of enzymes.     

Expression and evolution of hexamerins from the tobacco hornworm,Manduca sexta,and other Lepidoptera

Hexamerins are large hemolymph-proteins that accumulate during the late larval stages of insects. Hexamerins have emerged from hemocyanin, but have lost the ability to bind oxygen. Hexamerins are mainly considered as storage proteins for non-feeding stages, but may also have other functions, e.g. in cuticle formation, transport and immune response. The genome of the hornworm Manduca sexta harbors six hexamerin genes. Two of them code for arylphorins (Msex2.01690, Msex2.15504) and two genes correspond to a methionine-rich hexamerin (Msex2.10735) and a moderately methionine-rich hexamerin (Msex2.01694), respectively. Two other genes do not correspond to any known hexamerin and distantly resemble the arylphorins (Msex2.01691, Msex2.01693). Five of the six hexamerin genes are clustered within ∼45 kb on scaffold 00023, which shows conserved synteny in various lepidopteran genomes. The methionine-rich hexamerin gene is located at a distinct site. M. sexta and other Lepidoptera have lost the riboflavin-binding hexamerin. With the exception of Msex2.01691, which displays low mRNA levels throughout the life cycle, all hexamerins are most highly expressed during pre-wandering phase of the 5th larval instar of M. sexta, supporting their role as storage proteins. Notably, Msex2.01691 is most highly expressed in the brain, suggesting a divergent function. Phylogenetic analyses showed that hexamerin evolution basically follows insect systematics. Lepidoptera display an unparalleled diversity of hexamerins, which exceeds that of other hexapod orders. In contrast to previous analyses, the lepidopteran hexamerins were found monophyletic. Five distinct types of hexamerins have been identified in this order, which differ in terms of amino acid composition and evolutionary history: i. the arylphorins, which are rich in aromatic amino acids (∼20% phenylalanine and tyrosine), ii. the distantly related arylphorin-like hexamerins, iii. the methionine-rich hexamerins, iv. the moderately methionine rich hexamerins, and v. the riboflavin-binding hexamerins.     

Integrated modeling of protein-coding genes in the Manduca sexta genome using RNA-Seq data from the biochemical model insect

The genome sequence of Manduca sexta was recently determined using 454 technology. Cufflinks and MAKER2 were used to establish gene models in the genome assembly based on the RNA-Seq data and other species' sequences. Aided by the extensive RNA-Seq data from 50 tissue samples at various life stages, annotators over the world (including the present authors) have manually confirmed and improved a small percentage of the models after spending months of effort. While such collaborative efforts are highly commendable, many of the predicted genes still have problems which may hamper future research on this insect species. As a biochemical model representing lepidopteran pests, M. sexta has been used extensively to study insect physiological processes for over five decades. In this work, we assembled Manduca datasets Cufflinks 3.0, Trinity 4.0, and Oases 4.0 to assist the manual annotation efforts and development of Official Gene Set (OGS) 2.0. To further improve annotation quality, we developed methods to evaluate gene models in the MAKER2, Cufflinks, Oases and Trinity assemblies and selected the best ones to constitute MCOT 1.0 after thorough crosschecking. MCOT 1.0 has 18,089 genes encoding 31,666 proteins: 32.8% match OGS 2.0 models perfectly or near perfectly, 11,747 differ considerably, and 29.5% are absent in OGS 2.0. Future automation of this process is anticipated to greatly reduce human efforts in generating comprehensive, reliable models of structural genes in other genome projects where extensive RNA-Seq data are available.     

An independent occurrence of the chimeric P450 enzyme CYP337B3 of Helicoverpa armigera confers cypermethrin resistance in Pakistan

The increasing resistance level of insect pest species is a major concern to agriculture worldwide. The cotton bollworm, Helicoverpa armigera, is one of the most important pest species due to being highly polyphagous, geographically widespread, and resistant towards many chemical classes of insecticides. We previously described the mechanism of fenvalerate resistance in Australian populations conferred by the chimeric cytochrome P450 monooxygenase CYP337B3, which arose by unequal crossing-over between CYP337B1 and CYP337B2. Here, we show that this mechanism is also present in the cypermethrin-resistant FSD strain from Pakistan. The Pakistani and the Australian CYP337B3 alleles differ by 18 synonymous and three nonsynonymous SNPs and additionally in the length and sequence of the intron. Nevertheless, the activity of both CYP337B3 proteins is comparable. We demonstrate that CYP337B3 is capable of metabolizing cypermethrin (trans- and especially cis-isomers) to the main metabolite 4'-hydroxycypermethrin, which exhibits no intrinsic toxicity towards susceptible larvae. In a bioassay, CYP337B3 confers a 7-fold resistance towards cypermethrin in FSD larvae compared to susceptible larvae from the Australian TWB strain lacking CYP337B3. Linkage analysis shows that presence of CYP337B3 accounts for most of the cypermethrin resistance in the FSD strain; up-regulation of other P450s in FSD plays no detectable role in resistance. The presence or absence of CYP337B3 can be easily detected by a simple PCR screen, providing a powerful tool to rapidly distinguish resistant from susceptible individuals in the field and to determine the geographical distribution of this resistance gene. Our results suggest that CYP337B3 evolved twice independently by unequal crossing-over between CYP337B2 and two different CYP337B1 alleles.     

Structural features,evolutionary relationships,and transcriptional regulation of C-type lectin-domain proteins in Manduca sexta

C-type lectins (CTLs) are a large family of Ca²⁺-dependent carbohydrate-binding proteins recognizing various glycoconjugates and functioning primarily in immunity and cell adhesion. We have identified 34 CTLDP (for CTL-domain protein) genes in the Manduca sexta genome, which encode proteins with one to three CTL domains. CTL-S1 through S9 (S for simple) have one or three CTL domains; immulectin-1 through 19 have two CTL domains; CTL-X1 through X6 (X for complex) have one or two CTL domains along with other structural modules. Nine simple CTLs and seventeen immulectins have a signal peptide and are likely extracellular. Five complex CTLs have both an N-terminal signal peptide and a C-terminal transmembrane region, indicating that they are membrane anchored. Immulectins exist broadly in Lepidoptera and lineage-specific gene duplications have generated three clusters of fourteen genes in the M. sexta genome, thirteen of which have similar expression patterns. In contrast to the family expansion, CTL-S1∼S6, S8, and X1∼X6 have 1:1 orthologs in at least four lepidopteran/dipteran/coleopteran species, suggestive of conserved functions in a wide range of holometabolous insects. Structural modeling suggests the key residues for Ca²⁺-dependent or independent binding of certain carbohydrates by CTL domains. Promoter analysis identified putative κB motifs in eighteen of the CTL genes, which did not have a strong correlation with immune inducibility in the mRNA or protein levels. Together, the gene identification, sequence comparisons, structure modeling, phylogenetic analysis, and expression profiling establish a solid foundation for future studies of M. sexta CTL-domain proteins.     

The silkworm glutathione S-transferase gene noppera-bo is required for ecdysteroid biosynthesis and larval development

Insect molting and metamorphosis are tightly controlled by ecdysteroids, which are important steroid hormones that are synthesized from dietary sterols in the prothoracic gland. One of the ecdysteroidogenic genes in the fruit fly Drosophila melanogaster is noppera-bo (nobo), also known as GSTe14, which encodes a member of the epsilon class of glutathione S-transferases. In D. melanogaster, nobo plays a crucial role in utilizing cholesterol via regulating its transport and/or metabolism in the prothoracic gland. However, it is still not known whether the orthologs of nobo from other insects are also involved in ecdysteroid biosynthesis via cholesterol transport and/or metabolism in the prothoracic gland. Here we report genetic evidence showing that the silkworm Bombyx mori ortholog of nobo (nobo-Bm; GSTe7) is essential for silkworm development. nobo-Bm is predominantly expressed in the prothoracic gland. To assess the functional importance of nobo-Bm, we generated a B. mori genetic mutant of nobo-Bm using TALEN-mediated genome editing. We show that loss of nobo-Bm function causes larval arrest and a glossy cuticle phenotype, which are rescued by the application of 20-hydroxyecdysone. Moreover, the prothoracic gland cells isolated from the nobo-Bm mutant exhibit an abnormal accumulation of 7-dehydrocholesterol, a cholesterol metabolite. These results suggest that the nobo family of glutathione S-transferases is essential for development and for the regulation of sterol utilization in the prothoracic gland in not only the Diptera but also the Lepidoptera. On the other hand, loss of nobo function mutants of D. melanogaster and B. mori abnormally accumulates different sterols, implying that the sterol utilization in the PG is somewhat different between these two insect species.     

Cyanide detoxification in an insect herbivore: Molecular identification of β-cyanoalanine synthases from Pieris rapae

Cyanogenic compounds occur widely in the plant kingdom. Therefore, many herbivores are adapted to the presence of these compounds in their diet by either avoiding cyanide release or by efficient cyanide detoxification mechanisms. The mechanisms of adaptation are not fully understood. Larvae of Pieris rapae (Lepidoptera: Pieridae) are specialist herbivores on glucosinolate-containing plants. They are exposed to cyanide during metabolism of phenylacetonitrile, a product of benzylglucosinolate breakdown catalyzed by plant myrosinases and larval nitrile-specifier protein (NSP) in the gut. Cyanide is metabolized to β-cyanoalanine and thiocyanate in the larvae. Here, we demonstrate that larvae of P. rapae possess β-cyanoalanine activity in their gut. We have identified three gut-expressed cDNAs designated PrBSAS1-PrBSAS3 which encode proteins with similarity to β-substituted alanine synthases (BSAS). Characterization of recombinant PrBSAS1-PrBSAS3 shows that they possess β-cyanoalanine activity. In phylogenetic trees, PrBSAS1-PrBSAS3, the first characterized insect BSAS, group together with a characterized mite β-cyanoalanine synthase and bacterial enzymes indicating a similar evolutionary history.     

Potential detoxification of gossypol by UDP-glycosyltransferases in the two Heliothine moth species Helicoverpa armigera and Heliothis virescens

The cotton bollworm Helicoverpa armigera and the tobacco budworm Heliothis virescens are closely related generalist insect herbivores and serious pest species on a number of economically important crop plants including cotton. Even though cotton is well defended by its major defensive compound gossypol, a toxic sesquiterpene dimer, larvae of both species are capable of developing on cotton plants. In spite of severe damage larvae cause on cotton plants, little is known about gossypol detoxification mechanisms in cotton-feeding insects. Here, we detected three monoglycosylated and up to five diglycosylated gossypol isomers in the feces of H. armigera and H. virescens larvae fed on gossypol-supplemented diet. Candidate UDP-glycosyltransferase (UGT) genes of H. armigera were selected by microarray studies and in silico analyses and were functionally expressed in insect cells. In enzymatic assays, we show that UGT41B3 and UGT40D1 are capable of glycosylating gossypol mainly to the diglycosylated gossypol isomer 5 that is characteristic for H. armigera and is absent in H. virescens feces. In conclusion, our results demonstrate that gossypol is partially metabolized by UGTs via glycosylation, which might be a crucial step in gossypol detoxification in generalist herbivores utilizing cotton as host plant.     

Midgut aminopeptidase N isoforms from Ostrinia nubilalis: Activity characterization and differential binding to Cry1Ab and Cry1Fa proteins from Bacillus thuringiensis

Aminopeptidase N (APN) isoforms from Lepidoptera are known for their involvement in the mode of action of insecticidal Cry proteins from Bacillus thuringiensis. These enzymes belong to a protein family with at least eight different members that are expressed simultaneously in the midgut of lepidopteran larvae. Here, we focus on the characterization of the APNs from Ostrinia nubilalis (OnAPNs) to identify potential Cry receptors. We expressed OnAPNs in insect cells using a baculovirus system and analyzed their enzymatic activity by probing substrate specificity and inhibitor susceptibility. The interaction with Cry1Ab and Cry1Fa proteins (both found in transgenic insect-resistant maize) was evaluated by ligand blot assays and immunocytochemistry. Ligand blots of brush border membrane proteins showed that both Cry proteins bound mainly to a 150 kDa-band, in which OnAPNs were greatly represented. Binding analysis of Cry proteins to the cell-expressed OnAPNs showed that OnAPN1 interacted with both Cry1Ab and Cry1Fa, whereas OnAPN3a and OnAPN8 only bound to Cry1Fa. Two isoforms, OnAPN2 and OnAPN3b, did not interact with any of these two proteins. This work provides the first evidence of a differential role of OnAPN isoforms in the mode of action of Cry proteins in O. nubilalis.