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

Genome-wide identification, classification, and expression profiling of serine esterases and other esterase-related proteins in the tobacco hornworm, Manduca sexta

Serine esterases (SEs) are hydrolases that catalyze the conversion of carboxylic esters into acids and alcohols. Lipases and carboxylesterases constitute two major groups of SEs. Although over a hundred of insect genomes are known, systematic identification and classification of SEs are rarely performed, likely due to large size and complex composition of the gene family in each species. Considering their key roles in lipid metabolism and other physiological processes, we have categorized 144 M. sexta SEs and SE homologs (SEHs), 114 of which contain a motif of GXSXG. Multiple sequence alignment and phylogenetic tree analysis have revealed 39 neutral lipases (NLs), 3 neutral lipase homologs (NLHs), 11 acidic lipases (ALs), 3 acidic lipase homologs (ALHs), a lipase-3, a triglyceride lipase, a monoglyceride lipase, a hormone-sensitive lipase, and a GDSL lipase. Eighty-three carboxylesterase genes encode 29 α-esterases (AEs), 12 AEHs (e.g., SEH4-1-3), 20 feruloyl esterases (FEs), 2 FEHs, 2 β-esterases (BEs), 2 integument esterases (IEs), 1 IEH, 4 juvenile hormone esterases, 2 acetylcholinesterases, gliotactin, 6 neuroligins, neurotactin, and an uncharacteristic esterase homolog. In addition to these GXSXG proteins, we have identified 26 phospholipases and 13 thioesterases. Expression profiling of these genes in specific tissues and stages has provided insights into their functions including digestion, detoxification, hormone processing, neurotransmission, reproduction, and developmental regulation. In summary, we have established a framework of information on SEs and related proteins in M. sexta to stimulate their research in the model species and comparative investigations in agricultural pests or disease vectors.

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.

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.

Digestion-related proteins in the tobacco hornworm, Manduca sexta

Food digestion is vital for the survival and prosperity of insects. Research on insect digestive enzymes yields knowledge of their structure and function, and potential targets of antifeedants to control agricultural pests. While such enzymes from pest species are more relevant for inhibitor screening, a systematic analysis of their counterparts in a model insect has broader impacts. In this context, we identified a set of 122 digestive enzyme genes from the genome of Manduca sexta, a lepidopteran model related to some major agricultural pests. These genes encode hydrolases of proteins (85), lipids (20), carbohydrates (16), and nucleic acids (1). Gut serine proteases (62) and their noncatalytic homologs (11) in the S1A subfamily are encoded by abundant transcripts whose levels correlate well with larval feeding stages. Aminopeptidases (10), carboxypeptidases (10), and other proteases (3) also participate in dietary protein digestion. A large group of 11 lipases as well as 9 esterases are probably responsible for digesting lipids in diets. The repertoire of carbohydrate hydrolases (16) is relatively small, including two amylases, three maltases, two sucrases, two α-glucosidases, and others. Lysozymes, peptidoglycan amidases, and β-1,3-glucanase may hydrolyze peptidoglycans and glucans to harvest energy and defend the host from microbes on plant leaves. One alkaline nuclease is associated with larval feeding, which is likely responsible for hydrolyzing denatured DNA and RNA undergoing autolysis at a high pH of midgut. Proteomic analysis of the ectoperitrophic fluid from feeding larvae validated at least 131 or 89% of the digestive enzymes and their homologs. In summary, this study provides for the first time a holistic view of the digestion-related proteins in a lepidopteran model insect and clues for comparative research in lepidopteran pests and beyond.

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.

The Manduca sexta serpinome: Analysis of serpin genes and proteins in the tobacco hornworm

Members of the serpin superfamily of proteins occur in animals, plants, bacteria, archaea and some viruses. They adopt a variety of physiological functions, including regulation of immune system, modulation of apoptosis, hormone transport and acting as storage proteins. Most members of the serpin family are inhibitors of serine proteinases. In this study, we searched the genome of Manduca sexta and identified 32 serpin genes. We analyzed the structure of these genes and the sequences of their encoded proteins. Three M. sexta genes (serpin-1, serpin-15, and serpin-28) have mutually exclusive alternatively spliced exons encoding the carboxyl-terminal reactive center loop of the protein, which is the site of interaction with target proteases. We discovered that MsSerpin-1 has 14 splicing isoforms, including two undiscovered in previous studies. Twenty-eight of the 32 M. sexta serpins include a putative secretion signal peptide and are predicted to be extracellular proteins. Phylogenetic analysis of serpins in M. sexta and Bombyx mori indicates that 17 are orthologous pairs, perhaps carrying out essential physiological functions. Analysis of the reactive center loop and hinge regions of the protein sequences indicates that 16 of the serpin genes encode proteins that may lack proteinase inhibitor activity. Our annotation and analysis of these serpin genes and their transcript profiles should lead to future advances in experimental study of their functions in insect biochemistry.

WaspBase: a genomic resource for the interactions among parasitic wasps, insect hosts and plants

Insect pests reduce yield and cause economic losses, which are major problems in agriculture. Parasitic wasps are the natural enemies of many agricultural pests and thus have been widely used as biological control agents. Plants, phytophagous insects and parasitic wasps form a tritrophic food chain. Understanding the interactions in this tritrophic system should be helpful for developing parasitic wasps for pest control and deciphering the mechanisms of parasitism. However, the genomic resources for this tritrophic system are not well organized. Here, we describe the WaspBase, a new database that contains 573 transcriptomes of 35 parasitic wasps and the genomes of 12 parasitic wasps, 5 insect hosts and 8 plants. In addition, we identified long non-coding RNA, untranslated regions and 25 widely studied gene families from the genome and transcriptome data of these species. WaspBase provides conventional web services such as Basic Local Alignment Search Tool, search and download, together with several widely used tools such as profile hidden Markov model, Multiple Alignment using Fast Fourier Transform, automated alignment trimming and JBrowse. We also present a collection of active researchers in the field of parasitic wasps, which should be useful for constructing scientific networks in this field.

An analysis of 67 RNA-seq datasets from various tissues at different stages of a model insect, Manduca sexta

Background

Manduca sexta is a large lepidopteran insect widely used as a model to study biochemistry of insect physiological processes. As a part of its genome project, over 50 cDNA libraries have been analyzed to profile gene expression in different tissues and life stages. While the RNA-seq data were used to study genes related to cuticle structure, chitin metabolism and immunity, a vast amount of the information has not yet been mined for understanding the basic molecular biology of this model insect. In fact, the basic features of these data, such as composition of the RNA-seq reads and lists of library-correlated genes, are unclear. From an extended view of all insects, clear-cut tempospatial expression data are rarely seen in the largest group of animals including Drosophila and mosquitoes, mainly due to their small sizes.

Results

We obtained the transcriptome data, analyzed the raw reads in relation to the assembled genome, and generated heatmaps for clustered genes. Library characteristics (tissues, stages), number of mapped bases, and sequencing methods affected the observed percentages of genome transcription. While up to 40% of the reads were not mapped to the genome in the initial Cufflinks gene modeling, we identified the causes for the mapping failure and reduced the number of non-mappable reads to <8%. Similarities between libraries, measured based on library-correlated genes, clearly identified differences among tissues or life stages. We calculated gene expression levels, analyzed the most abundantly expressed genes in the libraries. Furthermore, we analyzed tissue-specific gene expression and identified 18 groups of genes with distinct expression patterns.

Conclusion

We performed a thorough analysis of the 67 RNA-seq datasets to characterize new genomic features of M. sexta. Integrated knowledge of gene functions and expression features will facilitate future functional studies in this biochemical model insect.

Electronic supplementary material

The online version of this article doi: (10.1186/s12864-017-4147-y) contains supplementary material, which is available to authorized users.

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

Insect serine proteases (SPs) and serine protease homologs (SPHs) participate in digestion, defense, development, and other physiological processes. In mosquitoes, some clip-domain SPs and SPHs (i.e. CLIPs) have been investigated for possible roles in antiparasitic responses. In a recent test aimed at improving quality of gene models in the Anopheles gambiae genome using RNA-seq data, we observed various discrepancies between gene models in AgamP4.5 and corresponding sequences selected from those modeled by Cufflinks, Trinity and Bridger. Here we report a comparative analysis of the 337 SP-related proteins in A. gambiae by examining their domain structures, sequence diversity, chromosomal locations, and expression patterns. One hundred and ten CLIPs contain 1 to 5 clip domains in addition to their protease domains (PDs) or non-catalytic, protease-like domains (PLDs). They are divided into five subgroups: CLIPAs (22) are clip_1-5-PLD; CLIPBs (29), CLIPCs (12) and CLIPDs (14) are mainly clip-PD; most CLIPEs (33) have a domain structure of PD/PLD-PLD-clip-PLD_0-1. While expression of the CLIP genes in group-1 is generally low and detected in various tissue- and stage-specific RNA-seq libraries, some putative GPs/GPHs (i.e. single domain gut SPs/SPHs) in group-2 are highly expressed in midgut, whole larva or whole adult libraries. In comparison, 46 SPs, 26 SPHs, and 37 multi-domain SPs/SPHs (i.e. PD/PLD-PLD_≥1) in group-3 do not seem to be specifically expressed in digestive tract. There are 16 SPs and 2 SPH containing other types of putative regulatory domains (e.g. LDLa, CUB, Gd). Of the 337 SP and SPH genes, 159 were sorted into 46 groups (2-8 members/group) based on similar phylogenetic tree position, chromosomal location, and expression profile. This information and analysis, including improved gene models and protein sequences, constitute a solid foundation for functional analysis of the SP-related proteins in A. gambiae.

Transcriptome Analysis of Male Drosophila melanogaster Exposed to Ethylparaben Using Digital Gene Expression Profiling

Ethylparaben (EP) has been shown to have estrogenic effects and can affect the normal development, longevity, and reproductive system of some animals. In this study, we investigated the effects of EP in male Drosophila melanogaster using transcriptome analysis or digital gene expression profiling. We then screened differentially expressed genes (DEGs) between the two groups (EP-treated and control group) of Drosophila, and performed clustering analysis, gene ontology (GO) function annotation, kyoto encyclopedia of gene and genomes metabolic pathway analysis. We found that EP enriched GO in three processes: cellular component, molecular function, and biological process. Consequently, we detected 13,959 genes and among them, 18 genes were identified to be significantly expressed between the EP-treated and control samples. Of these, seven genes were down-regulated, and eleven genes were up-regulated in EP-treated samples. Furthermore, four DEGs including two down-regulated genes (CG9465, CG9468) and two up-regulated genes (TotA, Sqz) were verified by real-time quantitative PCR. This study revealed the impact of EP on gene expression in fruit fly and provided new insight into the mechanisms of this response, which is helpful for understanding EP toxicity to humans.

Improved annotation of the insect vector of citrus greening disease: biocuration by a diverse genomics community

Database URL

https://citrusgreening.org/.

Multifaceted biological insights from a draft genome sequence of the tobacco hornworm moth, Manduca sexta

Manduca sexta, known as the tobacco hornworm or Carolina sphinx moth, is a lepidopteran insect that is used extensively as a model system for research in insect biochemistry, physiology, neurobiology, development, and immunity. One important benefit of this species as an experimental model is its extremely large size, reaching more than 10 g in the larval stage. M. sexta larvae feed on solanaceous plants and thus must tolerate a substantial challenge from plant allelochemicals, including nicotine. We report the sequence and annotation of the M. sexta genome, and a survey of gene expression in various tissues and developmental stages. The Msex_1.0 genome assembly resulted in a total genome size of 419.4 Mbp. Repetitive sequences accounted for 25.8% of the assembled genome. The official gene set is comprised of 15,451 protein-coding genes, of which 2498 were manually curated. Extensive RNA-seq data from many tissues and developmental stages were used to improve gene models and for insights into gene expression patterns. Genome wide synteny analysis indicated a high level of macrosynteny in the Lepidoptera. Annotation and analyses were carried out for gene families involved in a wide spectrum of biological processes, including apoptosis, vacuole sorting, growth and development, structures of exoskeleton, egg shells, and muscle, vision, chemosensation, ion channels, signal transduction, neuropeptide signaling, neurotransmitter synthesis and transport, nicotine tolerance, lipid metabolism, and immunity. This genome sequence, annotation, and analysis provide an important new resource from a well-studied model insect species and will facilitate further biochemical and mechanistic experimental studies of many biological systems in insects.

Changes in the Plasma Proteome of Manduca sexta Larvae in Relation to the Transcriptome Variations after an Immune Challenge: Evidence for High Molecular Weight Immune Complex Formation

Manduca sextais a lepidopteran model widely used to study insect physiological processes, including innate immunity. In this study, we explored the proteomes of cell-free hemolymph from larvae injected with a sterile buffer (C for control) or a mixture of bacteria (I for induced). Of the 654 proteins identified, 70 showed 1.67 to >200-fold abundance increases after the immune challenge; 51 decreased to 0-60% of the control levels. While there was no strong parallel between plasma protein levels and their transcript levels in hemocytes or fat body, the mRNA level changes (i.e.I/C ratios of normalized read numbers) in the tissues concurred with their protein level changes (i.e.I/C ratios of normalized spectral counts) with correlation coefficients of 0.44 and 0.57, respectively. Better correlations support that fat body contributes a more significant portion of the plasma proteins involved in various aspects of innate immunity. Consistently, ratios of mRNA and protein levels were better correlated for immunity-related proteins than unrelated ones. There is a set of proteins whose apparent molecular masses differ considerably from the calculatedMr's, suggestive of posttranslational modifications. In addition, some lowMrproteins were detected in the range of 80 to >300 kDa on a reducing SDS-polyacrylamide gel, indicating the existence of highMrcovalent complexes. We identified 30 serine proteases and their homologs, 11 of which are known members of an extracellular immune signaling network. Along with our quantitative transcriptome data, the protein identification, inducibility, and association provide leads toward a focused exploration of humoral immunity inM. sexta.

Phylogenetic analysis and expression profiling of the pattern recognition receptors: Insights into molecular recognition of invading pathogens in Manduca sexta

Pattern recognition receptors (PRRs) detect microbial pathogens and trigger innate immune responses. Previous biochemical studies have elucidated the physiological functions of eleven PRRs in Manduca sexta but our understanding of the recognition process is still limited, lacking genomic perspectives. While 34 C-type lectin-domain proteins and 16 Toll-like receptors are reported in the companion papers, we present here 120 other putative PRRs identified through the genome annotation. These include 76 leucine-rich repeat (LRR) proteins, 14 peptidoglycan recognition proteins, 6 EGF/Nim-domain proteins, 5 β-1,3-glucanase-related proteins, 4 galectins, 4 fibrinogen-related proteins, 3 thioester proteins, 5 immunoglobulin-domain proteins, 2 hemocytins, and 1 Reeler. Sequence alignment and phylogenetic analysis reveal the evolution history of a diverse repertoire of proteins for pathogen recognition. While functions of insect LRR proteins are mostly unknown, their structure diversification is phenomenal: In addition to the Toll homologs, 22 LRR proteins with a signal peptide are expected to be secreted; 18 LRR proteins lacking signal peptides may be cytoplasmic; 36 LRRs with a signal peptide and a transmembrane segment may be non-Toll receptors on the surface of cells. Expression profiles of the 120 genes in 52 tissue samples reflect complex regulation in various developmental stages and physiological states, including some likely by Rel family transcription factors via κB motifs in the promoter regions. This collection of information is expected to facilitate future biochemical studies detailing their respective roles in this model insect.

The immune signaling pathways of Manduca sexta

Signal transduction pathways and their coordination are critically important for proper functioning of animal immune systems. Our knowledge of the constituents of the intracellular signaling network in insects mainly comes from genetic analyses in Drosophila melanogaster. To facilitate future studies of similar systems in the tobacco hornworm and other lepidopteran insects, we have identified and examined the homologous genes in the genome of Manduca sexta. Based on 1:1 orthologous relationships in most cases, we hypothesize that the Toll, Imd, MAPK-JNK-p38 and JAK-STAT pathways are intact and operative in this species, as are most of the regulatory mechanisms. Similarly, cellular processes such as autophagy, apoptosis and RNA interference probably function in similar ways, because their mediators and modulators are mostly conserved in this lepidopteran species. We have annotated a total of 186 genes encoding 199 proteins, studied their domain structures and evolution, and examined their mRNA levels in tissues at different life stages. Such information provides a genomic perspective of the intricate signaling system in a non-drosophiline insect.

A genome-wide analysis of antimicrobial effector genes and their transcription patterns in Manduca sexta

Antimicrobial proteins/peptides (AMPs) are effectors of innate immune systems against pathogen infection in multicellular organisms. Over half of the AMPs reported so far come from insects, and these effectors act in concert to suppress or kill bacteria, fungi, viruses, and parasites. In this work, we have identified 86 AMP genes in the Manduca sexta genome, most of which seem likely to be functional. They encode 15 cecropins, 6 moricins, 6 defensins, 3 gallerimycins, 4 X-tox splicing variants, 14 diapausins, 15 whey acidic protein homologs, 11 attacins, 1 gloverin, 4 lebocins, 6 lysozyme-related proteins, and 4 transferrins. Some of these genes (e.g. attacins, cecropins) constitute large clusters, likely arising after rounds of gene duplication. We compared the amino acid sequences of M. sexta AMPs with their homologs in other insects to reveal conserved structural features and phylogenetic relationships. Expression data showed that many of them are synthesized in fat body and midgut during the larval-pupal molt. Certain genes contain one or more predicted κB binding sites and other regulatory elements in their promoter regions, which may account for the dramatic mRNA level increases in fat body and hemocytes after an immune challenge. Consistent with these strong mRNA increases, many AMPs become highly abundant in the larval plasma at 24 h after the challenge, as demonstrated in our previous peptidomic study. Taken together, these data suggest the existence of a large repertoire of AMPs in M. sexta, whose expression is up-regulated via immune signaling pathways to fight off invading pathogens in a coordinated manner.

Sequence conservation, phylogenetic relationships, and expression profiles of nondigestive serine proteases and serine protease homologs in Manduca sexta

Serine protease (SP) and serine protease homolog (SPH) genes in insects encode a large family of proteins involved in digestion, development, immunity, and other processes. While 68 digestive SPs and their close homologs are reported in a companion paper (Kuwar et al., in preparation), we have identified 125 other SPs/SPHs in Manduca sexta and studied their structure, evolution, and expression. Fifty-two of them contain cystine-stabilized structures for molecular recognition, including clip, LDLa, Sushi, Wonton, TSP, CUB, Frizzle, and SR domains. There are nineteen groups of genes evolved from relatively recent gene duplication and sequence divergence. Thirty-five SPs and seven SPHs contain 1, 2 or 5 clip domains. Multiple sequence alignment and molecular modeling of the 54 clip domains have revealed structural diversity of these regulatory modules. Sequence comparison with their homologs in Drosophila melanogaster, Anopheles gambiae and Tribolium castaneum allows us to classify them into five subfamilies: A are SPHs with 1 or 5 group-3 clip domains, B are SPs with 1 or 2 group-2 clip domains, C, D1 and D2 are SPs with a single clip domain in group-1a, 1b and 1c, respectively. We have classified into six categories the 125 expression profiles of SP-related proteins in fat body, brain, midgut, Malpighian tubule, testis, and ovary at different stages, suggesting that they participate in various physiological processes. Through RNA-Seq-based gene annotation and expression profiling, as well as intragenomic sequence comparisons, we have established a framework of information for future biochemical research of nondigestive SPs and SPHs in this model species.