The nucleotide sequence of a microvitellogenin encoding gene from the tobacco hornworm, Manduca sexta

LP30K protein manifested in hemocytes of Bombyx mori larva on Nosema bombycis infection and showed functional evolution based on glucose- binding domain

Infection by microsporidian Nosema bombycis induced appearance of exclusive protein conjugate of 190 kDa in hemocytes of silkworm Bombyx mori L (Lepidoptera: Bombycidae). Mass spectrometry of the band showed peptides of low molecular weight 30 kDa lipoprotein (LP30K). Six accessions of LP30K identified from the hemocytes comprised 30 K lipoprotein 1, 30 K protein 1, 2, 6, 7 and 11. Two uncharacterised proteins (UCP) identified from the hemocytes showed 100% similarity with LP30K sequence, altogether showed abundance after the infection. The LP30K accessions H9J4F6 (Q00802), E5EVW2 and the UCP accessions D4QGC0 and D4QGB9 showed presence of glucose binding protein I domain "ADSDVPNDILEEQLYNSIVVADYDSAVEK" that binds with fungal glucans to inhibit infection. However glucose binding protein II domain "TLAPRTDDVLAEQLYMSVVIGEYETAIAK" is absent in LP30K accessions from hemocytes showed loss of DNA sequences encoding the domain. The accessions H9J4F5, H9B440, A7LIK7 and H9B444 showed 92% identity with B. mori LP30K protein (NP_001095198.2) however the glucose binding domain I is absent in these accessions suggesting isoform- specific restricted fungal defense activity. Phylogeny tree of the LP30K homologues showed four groups including microvitellogenin and 30 kDa proteins showing functional diversity endorsed with evolutionary diversity. LP30K accessions with glucose binding domain diverged from that without glucose binding domain exemplify co-evolution for domain- dependent functional roles like storage and immune reactions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03685-x.

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.

Bom m 9 from Bombyx mori is a novel protein related to asthma

The silkworm (Bombyx mori) can cause severe IgE-mediated allergic disease, however, the mechanism remains unclear. The aim of this study was to investigate the immunologic mechanism by which silkworms induce allergy. Whole silkworm pupa proteins were separated by SDS-PAGE and 2-D PAGE. Then, IgE-binding proteins were detected by immunoblotting with sera of patients having an allergy to Bombyx mori. After tryptic digestion, the peptides of IgE-binding proteins were analyzed by matrix-assisted laser desorption ionization tandem time-of-flight mass spectrometry or tandem mass spectrometry. Database searches were used to identify allergens in silkworm pupa, after which Bom m 9 was to construct an asthma model. Thus, in the current study, a mouse asthma model was constructed with Bom m 9.

Identification and Characterization of a Novel Microvitellogenin from the Chinese Oak Silkworm Antheraea pernyi

Microvitellogenin (mVg) is a relatively small vitellogenic protein only characterized in the eggs of the lepidopteran insects Manduca sexta and Bombyx mori. In the present study, we report a novel mVg (ApmVg) isolated from the Chinese oak silkworm Antheraea pernyi. The obtained ApmVg cDNA sequence contains an open reading frame of 783 bp encoding a protein of 260 amino acids with a predicted molecular weight of 29.96 kDa. This gene does not contain introns. Structural analysis revealed that this protein shares putative conserved domains with the lepidopteran low-molecular weight lipoprotein, which belongs to the lipoprotein_11 superfamily. The protein sequence of ApmVg exhibits 48% sequence identity with mVg from M. sexta and 40-47% sequence identity with the 30K lipoproteins from B. mori. Phylogenetic analysis suggests that ApmVg is a novel member of the lepidopteran low-molecular weight lipoproteins. Transcriptional analysis indicated that ApmVg mRNA is mainly expressed in the fat body (both female and male) during post-diapause development of the pupal stage, and it was also detected in ovaries and spermaries in smaller amounts. RT-PCR and Western blot analyses revealed that ApmVg is synthesized by the fat body and secreted into hemolymph and ultimately accumulates in eggs. The ApmVg transcript can be detected in the fat bodies of female pupae four days after treatment with 20-hydroxyecdysone and shows an expression pattern distinct from that of vitellogenin (Vg), which is detectable throughout diapausing and in post-diapause development. ApmVg decreased dramatically during embryonic development. These results represent the first study of mVg outside M. sexta and B. mori and provide insight into the physiological role and evolution of mVgs.

Identification and Characterization of 30 K Protein Genes Found in Bombyx mori (Lepidoptera: Bombycidae) Transcriptome

The 30 K proteins, the major group of hemolymph proteins in the silkworm, Bombyx mori (Lepidoptera: Bombycidae), are structurally related with molecular masses of ∼30 kDa and are involved in various physiological processes, e.g., energy storage, embryonic development, and immune responses. For this report, known 30 K protein gene sequences were used as Blastn queries against sequences in the B. mori transcriptome (SilkTransDB). Twenty-nine cDNAs (Bm30K-1-29) were retrieved, including four being previously unidentified in the Lipoprotein_11 family. The genomic structures of the 29 genes were analyzed and they were mapped to their corresponding chromosomes. Furthermore, phylogenetic analysis revealed that the 29 genes encode three types of 30 K proteins. The members increased in each type is mainly a result of gene duplication with the appearance of each type preceding the differentiation of each species included in the tree. Real-Time Quantitative Polymerase Chain Reaction (Q-PCR) confirmed that the genes could be expressed, and that the three types have different temporal expression patterns. Proteins from the hemolymph was separated by SDS-PAGE, and those with molecular mass of ∼30 kDa were isolated and identified by mass spectrometry sequencing in combination with searches of various databases containing B. mori 30K protein sequences. Of the 34 proteins identified, 13 are members of the 30 K protein family, with one that had not been found in the SilkTransDB, although it had been found in the B. mori genome. Taken together, our results indicate that the 30 K protein family contains many members with various functions. Other methods will be required to find more members of the family.

Two crystal structures of Bombyx mori lipoprotein 3 - structural characterization of a new 30-kDa lipoprotein family member

The 30-kDa family of lipoproteins from insect hemolymph has been the focus of a number of studies over the last few years. Recently, four crystal structures of Bombyx mori lipoprotein 7 have been determined. Here we report two crystal structures of another member of the 30-kDa lipoprotein family, Bombyx mori lipoprotein 3 (Bmlp3). The protein was isolated from its natural source, mulberry silkworm hemolymph. It crystallized in two different crystal forms, Bmlp3-p21 (space group P21) and Bmlp3-c2 (space group C2). The crystal structures were solved by molecular replacement using the coordinates of Bmlp7 as a starting model. The crystals of Bmlp3-p21 diffracted X-rays to 2.4 Å resolution and of Bmlp3-c2 to 2.1 Å resolution. Bmlp3 has an overall fold characteristic of 30-kDa lipoproteins, with a VHS-type N-terminal domain and β-trefoil C-terminal domain. Structural comparison of Bmlp3 and Bmlp7 shows that the loops present in the C-terminal domain are flexible and participate in dimer formation. Additionally, new putative binding sites of Bmlp3 have been analyzed in detail and the electrostatic potential of the protein surface at physiological pH 7.4 conditions has been calculated. The results of these calculations are the starting point for an explanation of the recently reported cell-penetrating properties of the 30-kDa lipoproteins.

A proteomic view on the developmental transfer of homologous 30 kDa lipoproteins from peripheral fat body to perivisceral fat body via hemolymph in silkworm, Bombyx mori

BACKGROUND

A group of abundant proteins of ~30 kDa is synthesized in silkworm larval peripheral fat body (PPFB) tissues and transported into the open circulatory system (hemolymph) in a time-depended fashion to be eventually stored as granules in the pupal perivisceral fat body (PVFB) tissues for adult development during the non-feeding stage. These proteins have been shown to act anti-apoptotic besides being assigned roles in embryogenesis and defense. However, detailed protein structural information for individual PPFB and PVFB tissues during larval and pupal developmental stages is still missing. Gel electrophoresis and chromatography were used to separate the 30 kDa proteins from both PPFB and PVFB as well as hemolymph total proteomes. Mass spectrometry (MS) was employed to elucidate individual protein sequences. Furthermore, 30 kDa proteins were purified and biochemically characterized.

RESULTS

One- and two-dimensional gel electrophoresis (1/2D-PAGE) was used to visualize the relative changes of abundance of the 30 kDa proteins in PPFB and PVFB as well as hemolymph from day 1 of V instar larval stage to day 6 of pupal stage. Their concentrations were markedly increased in hemolymph and PVFB up to the first two days of pupal development and these proteins were consumed during development of the adult insect. Typically, three protein bands were observed (~29, 30, 31 kDa) in 1D-PAGE, which were subjected to MS-based protein identification along with spots excised from 2D-gels run for those proteomes. Gas phase fragmentation was used to generate peptide sequence information, which was matched to the available nucleotide data pool of more than ten highly homologous insect 30 kDa lipoproteins. Phylogenetic and similarity analyses of those sequences were performed to assist in the assignment of experimentally identified peptides to known sequences. Lipoproteins LP1 to LP5 and L301/302 could be matched to peptides extracted from all bands suggesting the presence of full length and truncated or modified protein forms in all of them. The individual variants could not be easily separated by classical means of purification such as 2D-PAGE because of their high similarity. They even seemed to aggregate as was indicated by native gel electrophoresis. Multistep chromatographic procedures eventually allowed purification of an LP3-like protein. The protein responded to lipoprotein-specific staining.

CONCLUSIONS

In B. mori larvae and pupae, 30 kDa lipoproteins LP1 to LP5 and L301/302 were detected in PPFB and PVFB tissue as well as in hemolymph. The concentration of these proteins changed progressively during development from their synthesis in PPFB, transport in hemolymph to storage in PVFB. While the 30 kDa proteins could be reproducibly separated in three bands electrophoretically, the exact nature of the individual protein forms present in those bands remained partially ambiguous. The amino acid sequences of all known 30 kDa proteins showed very high homology. High-resolution separation techniques will be necessary before MS and other structural analysis can shed more light on the complexity of the 30 kDa subproteome in B. mori. A first attempt to that end allowed isolation of a B. mori LP3-like protein, the complete structure, properties and function of which will now be elucidated in detail.

The homeodomain protein PBX participates in JH-related suppressive regulation on the expression of major plasma protein genes in the silkworm, Bombyx mori

In the silkworm, Bombyx mori, major plasma proteins referred to as 30K proteins are the most abundant proteins in the hemolymph of final (fifth) instar larvae. Surgical extirpation of corpora allata, the source of a juvenile hormone (JH), causes rapid accumulation of 30K proteins in the hemolymph of fourth instar larvae. The 30K protein 6G1 (30K6G1) gene was repressed in primary cultured fat body cells treated with a JH analog (JHA), methoprene. To identify the JH response element present in the promoter region of the 30K6G1 gene, we performed transfection analyses of the 5'-deletion mutants of the 30K6G1 gene using primary cultured fat body cells, gel retardation assays and in vivo footprinting analysis. The results from those analyses revealed that a JH response element exists in the sequence between positions -147 and -140. When the promoter construct mutated at positions -143, -142, and -141 was transfected to fat body primary cultured cells, the suppression effect on the reporter gene expression caused by JHA was reduced. Gel retardation assay using specific antibody revealed that a PBX protein binds to the JH response element. Northern blot analysis revealed that the gene expression of Bombyx PBX is enhanced in the fat body cells by JHA treatment. These results indicate that PBX proteins are involved in the JH signaling pathway and play an important role in suppressing 30K protein gene expression in the fat body of B. mori.

Organization of serpin gene-1 from Manduca sexta. Evolution of a family of alternate exons encoding the reactive site loop

Manduca sexta serpin gene-1 encodes a family of serpins whose amino acid sequences are identical in their amino-terminal 336 residues but variable in their carboxyl-terminal 39-46 residues, which includes the reactive site loop (Jiang, H., Wang, Y., and Kanost, M. R. (1994) J. Biol. Chem. 269, 55-58). Here, we report the gene's complete nucleotide sequence and exon-intron structure. A unique characteristic of this gene is its exon 9, which is present in 12 alternate forms between exons 8 and 10. Isolation and characterization of cDNA clones containing exons 9C, 9H, and 9I, which were not found previously, indicate that all 12 alternate forms of exon 9 can be utilized to generate 12 different serpins. The splicing pathway apparently allows inclusion of only one exon 9 per molecule of mature serpin-1 mRNA. Analysis of exon-intron border sequences reveals unique features that may be involved in regulation of RNA splicing. The exon 9 region has apparently evolved through rounds of exon duplication and sequence divergence. The exons near the center of the region may have evolved recently, whereas the outermost exons are the most ancient. Exons 9G and 9H were duplicated as a pair from exons 9E and 9F, an event that may have occurred more than once in the history of this gene.

Structure of the Bombyx sorbitol dehydrogenase gene: a possible alternative use of the promoter

In an initial effort to understand the molecular mechanism of how low temperature induces sorbitol dehydrogenase gene expression in diapause eggs of the silkworm, the sorbitol dehydrogenase gene was isolated from a Bombyx genomic library using a cDNA encoding the Bombyx homologue of mammalian sorbitol dehydrogenase as a probe. The gene extended for about 10 kb, consisting of eight exons and seven introns. Four TATA motifs were found in the 5' upstream region of the gene, without CCAAT. AATTAA, instead of AATAAA, was localized in the upstream region of the polyadenylation site. Although a single copy of this gene was present per haploid genome, 1.2 kb and 1.1 kb transcripts were found from yolk cells in diapause eggs and from larval fat-body cells, respectively. The two major transcription initiation sites corresponding to both transcripts were localized at 355 and 226 base pairs upstream from the transition start site, indicating an alternative use of promoter. The 5'-upstream region of the gene contained a consensus sequence, TGA(A/T)AA(A/G/T), that has been found in insect genes expressed mainly in larval and pupal fat bodies. It also contained three kinds of sequences similar to cis-elements recognized by members of the steroid receptor superfamily, such as chicken ovalbumin upstream promoter transcription factor (COUP-TF)/Drosophila Seven up (SVP), Drosophila hormone receptor 39 (DHR39) and Bombyx fushi tarazu transcriptional factor 1 (BmFTZ-F1).

Isolation, characterization and complete nucleotide sequence of a Galleria mellonella cDNA coding for the follicle cell-specific yolk protein, YP4

An ovariole-specific yolk polypeptide cDNA from Galleria mellonella was isolated from a cDNA library constructed using vitellogenic ovariole poly (A)+ RNA, by differential screening with the ovariole and day-0 male pupal poly (A)+ RNA. The Ov11 cDNA hybridized to a 1.1 kb message from the ovariole but not from male pharate adults, ovariectomized pharate adults or day-5 last instar larvae. Developmental Northern analysis showed that this message is expressed only during vitellogenic stages. In situ hybridization of digoxigenin-dUTP labeled Ov11 cDNA to the day-8 pharate adult ovariole showed that only follicle cells express the Ov11 gene. The 960 bp Ov11 cDNA is nearly full-length, containing a single open reading frame coding for a 286 amino acid polypeptide with a hydrophobic signal sequence and two potential N-glycosylation and tyrosine sulfation sites. The fact that Ov11 cDNA codes for the 37 kDa ovariole-specific yolk protein, YP4, was confirmed by the N-terminal amino acid sequence analysis of purified YP4. Genomic Southern analysis suggests that the yp4 gene is a single copy gene and PCR analysis using the 5' and 3' end primers of the cDNA indicates that the yp4 gene is intronless within the coding region.

Juvenile hormone regulation and developmental expression of a Tenebrio desiccation stress protein gene

Levels of a 28 kDa hemolymph protein, desiccation stress protein (dsp28), in the mealworm beetle, Tenebrio molitor, increase in response to desiccation and cold stress and are also developmentally regulated under nonstress conditions. Dsp28 mRNA is produced in the fat body, and its abundance changes dramatically throughout development. Transcript abundance increases throughout larval development, drops at pupation, and increases again in adults, the highest levels being found in females. The juvenile hormone (JH) analogue, methoprene, increased dsp28 transcript accumulation in pupae, suggesting that changing JH titres have a role in developmental expression and also perhaps contribute to dsp28 regulation during conditions of environmental stress. Genomic DNA, containing the entire dsp28 coding region plus 1.3 kb of upstream sequence, was isolated and potential regulatory sequences, including putative JH response elements, were identified.

Sex-specific and hormone-controlled expression of a vitellogenin-encoding gene in the gypsy moth

Microvitellogenin and vitellogenin cDNA from Manduca sexta (tobacco hornworm) were tested for use as molecular probes to investigate the expression of genes coding for vitellogenins in Spodoptera frugiperda (fall armyworm) and Lymantria dispar (gypsy moth). Cross-hybridization was not observed between the M. sexta cDNAs and S. frugiperda DNA and mRNA. Vitellogenin cDNA from M. sexta did not hybridize to L. dispar DNA or mRNA. However, the 834 bp microvitellogenin cDNA from M. sexta hybridized to an approximately 850 bp transcript in L. dispar mRNA. A 2.5 kb cDNA clone, pz64, was isolated from late last instar larvae of female L. dispar by differential screening. This clone has 38% amino acid sequence (deduced) and 55% nucleic acid sequence similarities with the 3'-end of high molecular weight vitellogenin in Bombyx mori (silkworm). When used as a probe in northern analysis of L. dispar mRNA, this cDNA hybridized to a 5.3 kb transcript in female last instar larvae, pupae, and adults, but not to male last instar larvae and adults. This cDNA did not hybridize to mRNA from M. sexta or S. frugiperda. Expression of the 5.3 kb vitellogenin transcript hybridizing to the 2.5 kb cDNA clone was suppressed in 5-day-old last instar larvae of female L. dispar treated on day 2 with doses of the juvenile hormone analog, methoprene, greater than 10 nmol. Apparently, the high in vivo titer of juvenile hormone during the first 2 days of the last instar represses the transcription of vitellogenin mRNA.

Codon usage patterns among genes for lepidopteran hemolymph proteins

Patterns in codon usage were examined for the coding regions of the 23 known lepidopteran hemolymph proteins. Coding triplets are GC rich at the third position and a significant linear relationship between GC content of silent and nonsilent (replacement) sites was demonstrated. Intron GC content was significantly lower than in coding regions and no relationship between intron GC content and the same at silent and nonsilent sites was found. Though hemolymph proteins are all produced by the same tissue--fat body--significantly less bias was observed when all moth sequences were pooled than when sequences of the two major species were analyzed separately, as predicted by the genome hypothesis. In cases where no statistically significant bias was observed, polar or acidic/basic amino acids were almost exclusively involved. Calculation of codon adaptation indices (CAI) was of limited value in quantifying the degree of codon bias and probably reflects the complexity of multicellular-organism life cycles and the changing patterns of gene expression over different developmental stages.

An extraovarian protein accumulated in mosquito oocytes is a carboxypeptidase activated in embryos

We report a phenomenon previously unknown for oviparous animals; in Aedes aegypti mosquitoes a serine carboxypeptidase is synthesized extraovarially and then internalized by oocytes. The cDNA encoding mosquito vitellogenic carboxypeptidase (VCP) was cloned and sequenced. The VCP cDNA hybridizes to a 1.5-kilobase mRNA present only in the fat body of vitellogenic females. The deduced amino acid sequence of VCP shares significant homology with members of the serine carboxypeptidase family. Binding assays using a serine protease inhibitor, [3H]diisopropyl fluorophosphate, showed that VCP is activated in eggs at the onset of embryonic development. Activation of VCP is associated with the reduction in its size from 53 kDa (inactive proenzyme) to 48 kDa (active enzyme). The active, 48-kDa, form of VCP is maximally present at the middle of embryonic development and disappears by the end.

Structures and organization of major plasma protein genes of the silkworm Bombyx mori

In the silkworm, Bombyx mori, a group of structurally related proteins, termed 30K proteins, accumulate in the hemolymph of the last instar larvae. We have isolated and characterized three genes, each of which encodes a distinct 30K protein component. Each 30K protein gene is composed of a short first exon and a protein-coding second exon interspersed by a single intron. The transcription initiation site of the 30K protein mRNA was identified at the nucleotide level. A typical TATA box exists some 30 base-pairs upstream from the transcription initiation site. The 5'-flanking region of each gene also contains octamer-like sequences. Restriction mapping analyses revealed that the cloned 46 x 10(3) base-pair region of the chromosomal DNA bears three 30K protein genes. Several copies of highly reiterated retrotransposon-like sequences are present around the 30K protein genes. S1 nuclease protection analysis provided evidence that the biosynthesis of 30K protein is regulated in a stage-specific manner at the transcriptional level in the fat body.