Structural, evolutionary and functional analysis of APN genes in the Lepidoptera Bombyx mori

A chromosome-level genome assembly of the soybean pod borer: insights into larval transcriptional response to transgenic soybean expressing the pesticidal Cry1Ac protein

BACKGROUND

Genetically modified (GM) crop plants with transgenic expression of Bacillus thuringiensis (Bt) pesticidal proteins are used to manage feeding damage by pest insects. The durability of this technology is threatened by the selection for resistance in pest populations. The molecular mechanism(s) involved in insect physiological response or evolution of resistance to Bt is not fully understood.

RESULTS

To investigate the response of a susceptible target insect to Bt, the soybean pod borer, Leguminivora glycinivorella (Lepidoptera: Tortricidae), was exposed to soybean, Glycine max, expressing Cry1Ac pesticidal protein or the non-transgenic parental cultivar. Assessment of larval changes in gene expression was facilitated by a third-generation sequenced and scaffolded chromosome-level assembly of the L. glycinivorella genome (657.4 Mb; 27 autosomes + Z chromosome), and subsequent structural annotation of 18,197 RefSeq gene models encoding 23,735 putative mRNA transcripts. Exposure of L. glycinivorella larvae to transgenic Cry1Ac G. max resulted in prediction of significant differential gene expression for 204 gene models (64 up- and 140 down-regulated) and differential splicing among isoforms for 10 genes compared to unexposed cohorts. Differentially expressed genes (DEGs) included putative peritrophic membrane constituents, orthologs of Bt receptor-encoding genes previously linked or associated with Bt resistance, and those involved in stress responses. Putative functional Gene Ontology (GO) annotations assigned to DEGs were significantly enriched for 36 categories at GO level 2, respectively. Most significantly enriched cellular component (CC), biological process (BP), and molecular function (MF) categories corresponded to vacuolar and microbody, transport and metabolic processes, and binding and reductase activities. The DEGs in enriched GO categories were biased for those that were down-regulated (≥ 0.783), with only MF categories GTPase and iron binding activities were bias for up-regulation genes.

CONCLUSIONS

This study provides insights into pathways and processes involved larval response to Bt intoxication, which may inform future unbiased investigations into mechanisms of resistance that show no evidence of alteration in midgut receptors.

Culex quinquefasciatus alpha-glucosidase serves as a putative receptor of the Cry48Aa toxin from Lysinibacillus sphaericus

The Cry48Aa/Cry49Aa toxin of Lysinibacillus sphaericus shows specific toxicity towards larvae of Culex spp. Individual Cry48Aa and Cry49Aa subunits interact with distinct target sites in the larval midgut and overcome the resistance of Culex to the Bin toxin. However, the toxin-binding proteins have not yet been identified. The present study aimed to identify Cry48Aa-binding proteins in Culex quinquefasciatus. Pulldown assays using C. quinquefasciatus midgut brush-border membrane fractions (BBMFs) identified a class of proteins, including aminopeptidases (APNs), protease m1 zinc metalloproteases, alkaline phosphatases (ALPs), and maltases, that could be potentially involved in the mode of action of this toxin. RNA interference analysis showed that silenced larvae treated with dsRNA of the alpha-glucosidase (named Glu71) gene were more tolerant of the Cry48Aa/Cry49Aa toxin, which induced less than 20% mortality. The amino acid sequence of Glu71 exhibited 42% identity with Cqm1/Cpm1, which acted as a Bin toxin receptor. Toxin binding assays showed that Cry48Aa had a high specific binding capacity for the Glu71 protein, whereas Cry49Aa exhibited no specific binding. Overall, our results showed that Glu71 is a Cry48-binding protein involved in Cry48Aa/Cry49Aa toxicity.

Genome-wide identification and comparative analysis of Cry toxin receptor families in 7 insect species with a focus on Spodoptera litura

Cadherin, aminopeptidase N (APN) and alkaline phosphatase (ALP) have been characterized as Cry receptors. In this study, comparative genomic analysis of the 3 receptor families was performed in 7 insects. ALPs and APNs are divided into three and eight clades in phylogenetic trees, respectively. ALPs in clade 3 and APNs in clade 1 contain multiple paralogs within each species and most paralogs are located closely in chromosomes. Drosophila melanogaster has expanded APNs in clade 5 and were lowly expressed in midgut. Cadherins are divided into 16 clades; they may diverge before holometabolous insect speciation except for BtR and Cad89D-like clades. Eight insects from different orders containing BtR orthologs are sensitive to Cry1A or Cry3A, while five species without BtR are insensitive to both toxins. Most APNs in clade 1, several ALPs in clade 3, BtR and Cad89D-like genes were highly or moderately expressed in larval midgut of Spodoptera litura and the other six species, and several members in these clades have been identified as Cry receptors. Expressions of putative S. litura Cry receptors in the midgut after exposing to Bt toxins were also analyzed.

Cytotoxicity of Vibrio parahaemolyticus AHPND toxin on shrimp hemocytes, a newly identified target tissue, involves binding of toxin to aminopeptidase N1 receptor

Acute hepatopancreatic necrosis disease (AHPND) caused by PirAB^VP-producing strain of Vibrio parahaemolyticus, VP_AHPND, has seriously impacted the shrimp production. Although the VP_AHPND toxin is known as the VP_AHPND virulence factor, a receptor that mediates its action has not been identified. An in-house transcriptome of Litopenaeus vannamei hemocytes allows us to identify two proteins from the aminopeptidase N family, LvAPN1 and LvAPN2, the proteins of which in insect are known to be receptors for Cry toxin. The membrane-bound APN, LvAPN1, was characterized to determine if it was a VP_AHPND toxin receptor. The increased expression of LvAPN1 was found in hemocytes, stomach, and hepatopancreas after the shrimp were challenged with either VP_AHPND or the partially purified VP_AHPND toxin. LvAPN1 knockdown reduced the mortality, histopathological signs of AHPND in the hepatopancreas, and the number of virulent VP_AHPND bacteria in the stomach after VP_AHPND toxin challenge. In addition, LvAPN1 silencing prevented the toxin from causing severe damage to the hemocytes and sustained both the total hemocyte count (THC) and the percentage of living hemocytes. We found that the rLvAPN1 directly bound to both rPirA^VP and rPirB^VP toxins, supporting the notion that silencing of LvAPN1 prevented the VP_AHPND toxin from passing through the cell membrane of hemocytes. We concluded that the LvAPN1 was involved in AHPND pathogenesis and acted as a VP_AHPND toxin receptor mediating the toxin penetration into hemocytes. Besides, this was the first report on the toxic effect of VP_AHPND toxin on hemocytes other than the known target tissues, hepatopancreas and stomach.

A specific strain of Vibrio parahaemolyticus causing acute hepatopancreatic necrosis disease (AHPND) in shrimp or VP_AHPND produces a binary toxin (PirAB^vp toxin) that is previously known to induce cell death of stomach and hepatopancreas but the molecular mechanism has not been defined. Similar to Cry toxin receptor in insects, a novel aminopeptidase N1 protein from L. vannamei (LvAPN1) was identified as a putative receptor of VP_AHPND toxin. Suppression of LvAPN1 reduced the number of AHPND virulence plasmids in stomach and occurrence of AHPND clinical sign, sustained the number of total hemocyte count, and elevated the number of viable hemocyte. We demonstrated that VP_AHPND toxin challenge induces hemocyte cell damage and it interacts with LvAPN1 in vitro. Collectively, our finding suggested that not only stomach and hepatopancreas but also hemocyte are the VP_AHPND target tissues where LvAPN1 serves as a VP_AHPND toxin receptor. This study provides novel insight into the contributions of LvAPN1 receptor towards the AHPND pathogenesis in shrimp and may extend to the development of AHPND preventive measure in shrimp.

iTRAQ-Based Quantitative Proteomic Analysis of Digestive Juice across the First 48 Hours of the Fifth Instar in Silkworm Larvae

The silkworm is an oligophagous insect for which mulberry leaves are the sole diet. The nutrients needed for vital activities of the egg, pupal, and adult stages, and the proteins formed in the cocoon, are all derived from the larval stages. The silkworm feeds and grows quickly during the larval stages. In particular, the amount of leaf ingested and digested quickly increases from the ecdysis to the gluttonous stage in the fifth instar period. In this study, we used the iTRAQ proteomic technique to identify and analyze silkworm larval digestive juice proteins during this period. A total of 227 proteins were successfully identified. These were primarily serine protease activity, esterase activity, binding, and serine protease inhibitors, which were mainly involved in the digestion and overcoming the detrimental effects of mulberry leaves. Moreover, 30 genes of the identified proteins were expressed specifically in the midgut. Temporal proteomic analysis of digestive juice revealed developmental dynamic features related to molecular mechanisms of the principal functions of digesting, resisting pathogens, and overruling the inhibitory effects of mulberry leaves protease inhibitors (PIs) with a dynamic strategy, although overruling the inhibitory effects has not yet been confirmed by previous study. These findings will help address the potential functions of digestive juice in silkworm larvae.

Newly identified APN splice isoforms suggest novel splicing mechanisms may underlie circRNA circularization in moth

Circular RNA (circRNA) have long been considered by-products of splicing errors, but the coordination of RNA transcription and exon circularization events remains poorly understood. Here, we investigated this question using genes encoding aminopeptidases N (APNs), which are receptors of Bacillus thuringiensis toxins, in the cotton bollworm, Helicoverpa armigera. We cloned and sequenced the cDNA of ten APN genes (HaAPN1-10) located in the same APN gene cluster, and detected 20 and 14 novel splicing isoforms with exon skipping in HaAPN1 and HaAPN3, respectively, whereas no or very few variants were found in the remaining genes. Further study identified 14 and 6 circular RNA (circRNA) in HaAPN1 and HaAPN3, respectively. Neither novel splicing isoforms nor circRNA were detected in HaAPN2 and HaAPN5. Distinct from the conventional GT/AG splicing signal, short co-directional repeats were involved in the splicing of the linear and circular isoforms of HaAPN1 and HaAPN3. Identification of the splice sites revealed that the linear isoforms may be related in some way to the circularization. Moreover, phylogenetic analysis and detection of circRNA of the APN gene of the diamondback moth, Plutella xylostella (PxAPN3), showed that circRNA formation is relatively conserved during the lepidopteran evolutionary process. These results contribute to an improved understanding of lepidopteran APNs and this novel class of insect circRNA.

Gene Expression Pattern and Regulatory Network of α-Toxin Treatment in Bombyx mori

Bacillus bombyseptieus is a pathogen of Bombyx mori; it can cause bacterial septicemia in silkworm. One of the components of the parasporal crystal toxin of B. bombyseptieus, α-toxin, plays an important role in the process of infection in silkworm. In this study, we investigated the immune response of silkworm induced by α-toxin by using RNA-seq. We compared the changes in gene expression in the midgut, fatbody, and hemocytes of silkworm and in the B. mori embryonic cell line (BmE) after treatment with α-toxin and identified 952 differentially expressed genes and 353 differentially expressed long noncoding RNAs (lncRNAs). These regulated genes in different tissues were found to be enriched in different pathways. The upregulated genes in the midgut were mainly involved in peptidoglycan catabolic process and tyrosine kinase signaling pathway, whereas the downregulated genes were mainly involved in chitin metabolic pathways. The upregulated genes in fatbody were also involved in peptidoglycan catabolic process, but they were for a different peptidoglycan subtype. Further, genes encoding cecropins were enriched in the fatbody. The downregulated genes were mainly involved in the metabolic pathways of fundamental substances such as cellular protein metabolic process and nucleobase-containing compound metabolic process. These results suggest that α-toxin can induce various immune responses in silkworm, and further studies are warranted to understand the mechanism of α-toxin action in silkworm. Further, lncRNAs and differentially expressed genes were correlated using coexpression network analysis. Our findings revealed potential candidate genes and lncRNAs that might play important physiological functions in the immune response to α-toxins in silkworm.

Analysis of Homologs of Cry-toxin Receptor-Related Proteins in the Midgut of a Non-Bt Target, Nilaparvata lugens (Stål) (Hemiptera: Delphacidae)

The brown planthopper (BPH) Nilaparvata lugens is one of the most destructive insect pests in the rice fields of Asia. Like other hemipteran insects, BPH is not susceptible to Cry toxins of Bacillus thuringiensis (Bt) or transgenic rice carrying Bt cry genes. Lack of Cry receptors in the midgut is one of the main reasons that BPH is not susceptible to the Cry toxins. The main Cry-binding proteins (CBPs) of the susceptible insects are cadherin, aminopeptidase N (APN), and alkaline phosphatase (ALP). In this study, we analyzed and validated de novo assembled transcripts from transcriptome sequencing data of BPH to identify and characterize homologs of cadherin, APN, and ALP. We then compared the cadherin-, APN-, and ALP-like proteins of BPH to previously reported CBPs to identify their homologs in BPH. The sequence analysis revealed that at least one cadherin, one APN, and two ALPs of BPH contained homologous functional domains identified from the Cry-binding cadherin, APN, and ALP, respectively. Quantitative real-time polymerase chain reaction used to verify the expression level of each putative Cry receptor homolog in the BPH midgut indicated that the CBPs homologous APN and ALP were expressed at high or medium-high levels while the cadherin was expressed at a low level. These results suggest that homologs of CBPs exist in the midgut of BPH. However, differences in key motifs of CBPs, which are functional in interacting with Cry toxins, may be responsible for insusceptibility of BPH to Cry toxins.

Identification and Characterization of Hyphantria cunea Aminopeptidase N as a Binding Protein of Bacillus thuringiensis Cry1Ab35 Toxin

The fall webworm, Hyphantria cunea (Drury) is a major invasive pest in China. Aminopeptidase N (APN) isoforms in lepidopteran larvae midguts are known for their involvement in the mode of action of insecticidal crystal (Cry) proteins from Bacillus thuringiensis. In the present work, we identified a putative Cry1Ab toxin-binding protein, an APN isoform designated HcAPN3, in the midgut of H. cunea by ligand blot and mass spectrometry. HcAPN3 was highly expressed throughout all larval developmental stages and was abundant in the midgut and hindgut tissues. HcAPN3 was down-regulated at 6 h, then was up-regulated significantly at 12 h and 24 h after Cry1Ab toxin treatment. We expressed HcAPN3 in insect cells and detected its interaction with Cry1Ab toxin by ligand blot assays. Furthermore, RNA interference (RNAi) against HcAPN3 using oral delivery and injection of double-stranded RNA (dsRNA) resulted in a 61–66% decrease in transcript level. Down-regulating of the expression of HcAPN3 was closely associated with reduced susceptibility of H. cunea to Cry1Ab. In addition, the HcAPN3E fragment peptide expressed in Escherichia coli enhanced Cry1Ab toxicity against H. cunea larvae. This work represents the first evidence to suggest that an APN in H. cunea is a putative binding protein involved in Cry1Ab susceptibility.

Down-regulation of aminopeptidase N and ABC transporter subfamily G transcripts in Cry1Ab and Cry1Ac resistant Asian corn borer, Ostrinia furnacalis (Lepidoptera: Crambidae)

The Asian corn borer (ACB), Ostrinia furnacalis (Lepidoptera: Crambidae), is a highly destructive pest of cultivated maize throughout East Asia. Bacillus thuringiensis (Bt) crystalline protein (Cry) toxins cause mortality by a mechanism involving pore formation or signal transduction following toxin binding to receptors along the midgut lumen of susceptible insects, but this mechanism and mutations therein that lead to resistance are not fully understood. In the current study, quantitative comparisons were made among midgut expressed transcripts from O. furnacalis susceptible (ACB-BtS) and laboratory selected strains resistant to Cry1Ab (ACB-AbR) and Cry1Ac toxins (ACB-AcR) when feeding on non-Bt diet. From a combined de novo transcriptome assembly of 83,370 transcripts, ORFs of ≥ 100 amino acids were predicted and annotated for 28,940 unique isoforms derived from 12,288 transcripts. Transcriptome-wide expression estimated from RNA-seq read depths predicted significant down-regulation of transcripts for previously known Bt resistance genes, aminopeptidase N1 (apn1) and apn3, as well as a putative ATP binding cassette transporter group G (abcg) gene in both ACB-AbR and -AcR (log2[fold-change] ≥ 1.36; P < 0.0001). The transcripts that were most highly differentially regulated in both ACB-AbR and -AcR compared to ACB-BtS (log2[fold-change] ≥ 2.0; P < 0.0001) included up- and down-regulation of serine proteases, storage proteins and cytochrome P450 monooxygenases, as well as up-regulation of genes with predicted transport function. This study predicted the significant down-regulation of transcripts for previously known Bt resistance genes, aminopeptidase N1 (apn1) and apn3, as well as abccg gene in both ACB-AbR and -AcR. These data are important for the understanding of systemic differences between Bt resistant and susceptible genotypes.

Transcriptome analysis of interactions between silkworm and cytoplasmic polyhedrosis virus

Bombyx mori cytoplasmic polyhedrosis virus (BmCPV) specifically infects silkworm midgut (MG) and multiplication occurs mainly in posterior midgut (PM). In this study, MG and fat body (FB) were extracted at 0, 3, 24, and 72 h after BmCPV infection. The total sequence reads of each sample were more than 1510000, and the mapping ratio exceeded 95.3%. Upregulated transcripts increased in MG during the infection process. Gene ontology (GO) categories showed that antioxidants were all upregulated in FB but not in MG. BGI001299, BGI014434, BGI012068, and BGI009201 were MG-specific genes with transmembrane transport function, the expression of which were induced by BmCPV. BGI001299, BGI014434, and BGI012068 expressed in entire MG and may be involved in BmCPV invasion. BGI009201 expressed only in PM and may be necessary for BmCPV proliferation. BmPGRP-S2 and BGI012452 (a putative serine protease) were induced by BmCPV and may be involved in immune defense against BmCPV. The expression level of BmCPV S1, S2, S3, S6, and S7 was high and there was no expression of S9 in MG 72 h, implying that the expression time of structural protein coding genes is earlier. These results provide insights into the mechanism of BmCPV infection and host defense.

Bacillus bombysepticus α-Toxin Binding to G Protein-Coupled Receptor Kinase 2 Regulates cAMP/PKA Signaling Pathway to Induce Host Death

Bacterial pathogens and their toxins target host receptors, leading to aberrant behavior or host death by changing signaling events through subversion of host intracellular cAMP level. This is an efficient and widespread mechanism of microbial pathogenesis. Previous studies describe toxins that increase cAMP in host cells, resulting in death through G protein-coupled receptor (GPCR) signaling pathways by influencing adenylyl cyclase or G protein activity. G protein-coupled receptor kinase 2 (GRK2) has a central role in regulation of GPCR desensitization. However, little information is available about the pathogenic mechanisms of toxins associated with GRK2. Here, we reported a new bacterial toxin-Bacillus bombysepticus (Bb) α-toxin that was lethal to host. We showed that Bb α-toxin interacted with BmGRK2. The data demonstrated that Bb α-toxin directly bound to BmGRK2 to promote death by affecting GPCR signaling pathways. This mechanism involved stimulation of G_αs, increase level of cAMP and activation of protein kinase A (PKA). Activated cAMP/PKA signal transduction altered downstream effectors that affected homeostasis and fundamental biological processes, disturbing the structural and functional integrity of cells, resulting in death. Preventing cAMP/PKA signaling transduction by inhibitions (NF449 or H-89) substantially reduced the pathogenicity of Bb α-toxin. The discovery of a toxin-induced host death specifically linked to GRK2 mediated signaling pathway suggested a new model for bacterial toxin action. Characterization of host genes whose expression and function are regulated by Bb α-toxin and GRK2 will offer a deeper understanding of the pathogenesis of infectious diseases caused by pathogens that elevate cAMP.

Interference with regulation of host signaling by pathogens can alter gene expression, leading to functional disarray in the host cells that causes abnormal division or death. Here, we propose a previously undescribed model for how bacterial toxins subvert host processes via interaction with GRK2 that influences cAMP/PKA signaling. Our findings provide new fundamental information about how bacterial pathogens regulate host signal transduction to cause death, which offers additional perspectives in host-pathogen systems. These findings will help to advance our understanding of bacteria pathogenic mechanism. Furthermore, these might extend to other microbial pathogenesis and assist in designing new or safer strategies against pathogens.

The 5'-UTR intron of the midgut-specific BmAPN4 gene affects the level and location of expression in transgenic silkworms

Introns are important for regulating gene expression. BmAPN4, which has a 5'-UTR upstream intron (5 UI), is specifically expressed in the entire silkworm midgut. In our previous study, the promoter region upstream of the 5 UI of BmAPN4 was cloned and identified as the P3 promoter (P3P) with activity only in the anterior midgut. In this study, the sequence consisting of the P3P and the 5 UI was cloned and named as P3P+5 UI. A transgenic vector was constructed in which EGFP was controlled by P3P+5 UI. Transgenic P3+5 UI silkworms were generated by embryo microinjection. RT-PCR showed P3P+5 UI activity throughout the larval stage. Intense green fluorescence was seen only in the entire midgut of P3+5 UI silkworms and expression was confirmed by RT-PCR. qPCR revealed that expression of EGFP in the anterior midgut of P3+5 UI silkworms was 64% higher than in P3 silkworms, indicating the 5 UI sustained intron-mediated enhancement of gene expression. These results suggested that the BmAPN4 5 UI affected the level and site of expression. The 5 UI was cloned and added behind P2P, another specific promoter with activity only in the anterior midgut of silkworm, to construct the P2P+5 UI and transgenic P2+5 UI silkworms. Expression patterns were the same for P2P+5 UI and P2P, suggesting that the 5UI of BmAPN4 did not affect P2P. This study found that the BmAPN4 5 UI affected the amount and location of gene expression. Its influence appeared to be dependent on a specific promoter.

PC, a Novel Oral Insecticidal Toxin from Bacillus bombysepticus Involved in Host Lethality via APN and BtR-175

Insect pests have developed resistance to chemical insecticides, insecticidal toxins as bioinsecticides or genetic protection built into crops. Consequently, novel, orally active insecticidal toxins would be valuable biological alternatives for pest control. Here, we identified a novel insecticidal toxin, parasporal crystal toxin (PC), from Bacillus bombysepticus (Bb). PC shows oral pathogenic activity and lethality towards silkworms and Cry1Ac-resistant Helicoverpa armigera strains. In vitro assays, PC after activated by trypsin binds to BmAPN4 and BtR-175 by interacting with CR7 and CR12 fragments. Additionally, trypsin-activated PC demonstrates cytotoxicity against Sf9 cells expressing BmAPN4, revealing that BmAPN4 serves as a functional receptor that participates in Bb and PC pathogenicity. In vivo assay, knocking out BtR-175 increased the resistance of silkworms to PC. These data suggest that PC is the first protein with insecticidal activity identified in Bb that is capable of causing silkworm death via receptor interactions, representing an important advance in our understanding of the toxicity of Bb and the contributions of interactions between microbial pathogens and insects to disease pathology. Furthermore, the potency of PC as an insecticidal protein makes it a good candidate for inclusion in integrated agricultural pest management systems.

Complete Genome Sequence of Bacillus bombysepticus, a Pathogen Leading to Bombyx mori Black Chest Septicemia

Bacillus bombysepticus is a Gram-positive spore-forming bacterium. Here, we announce the first complete genome sequence of this organism isolated from the cadavers of silkworm larvae that had been sick. The genome contains a single circular chromosome and a circular plasmid. Analyses of the B. bombysepticus genome will provide insights into its pathomechanisms and biology.