Xiaolong Cao

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.

Analysis of chitin-binding proteins from Manduca sexta provides new insights into evolution of peritrophin A-type chitin-binding domains in insects

In insects, chitin is a major structural component of the cuticle and the peritrophic membrane (PM). In nature, chitin is always associated with proteins among which chitin-binding proteins (CBPs) are the most important for forming, maintaining and regulating the functions of these extracellular structures. In this study, a genome-wide search for genes encoding proteins with ChtBD2-type (peritrophin A-type) chitin-binding domains (CBDs) was conducted. A total of 53 genes encoding 56 CBPs were identified, including 15 CPAP1s (cuticular proteins analogous to peritrophins with 1 CBD), 11 CPAP3s (CPAPs with 3 CBDs) and 17 PMPs (PM proteins) with a variable number of CBDs, which are structural components of cuticle or of the PM. CBDs were also identified in enzymes of chitin metabolism including 6 chitinases and 7 chitin deacetylases encoded by 6 and 5 genes, respectively. RNA-seq analysis confirmed that PMP and CPAP genes have differential spatial expression patterns. The expression of PMP genes is midgut-specific, while CPAP genes are widely expressed in different cuticle forming tissues. Phylogenetic analysis of CBDs of proteins in insects belonging to different orders revealed that CPAP1s from different species constitute a separate family with 16 different groups, including 6 new groups identified in this study. The CPAP3s are clustered into a separate family of 7 groups present in all insect orders. Altogether, they reveal that duplication events of CBDs in CPAP1s and CPAP3s occurred prior to the evolutionary radiation of insect species. In contrast to the CPAPs, all CBDs from individual PMPs are generally clustered and distinct from other PMPs in the same species in phylogenetic analyses, indicating that the duplication of CBDs in each of these PMPs occurred after divergence of insect species. Phylogenetic analysis of these three CBP families showed that the CBDs in CPAP1s form a clearly separate family, while those found in PMPs and CPAP3s were clustered together in the phylogenetic tree. For chitinases and chitin deacetylases, most of phylogenetic analysis performed with the CBD sequences resulted in similar clustering to the one obtained by using catalytic domain sequences alone, suggesting that CBDs were incorporated into these enzymes and evolved in tandem with the catalytic domains before the diversification of different insect orders. Based on these results, the evolution of CBDs in insect CBPs is discussed to provide a new insight into the CBD sequence structure and diversity, and their evolution and expression in insects.

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.

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(2+)-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(2+)-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 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.

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 sexta is 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 calculated Mrs, suggestive of posttranslational modifications. In addition, some low Mr proteins were detected in the range of 80 to >300 kDa on a reducing SDS-polyacrylamide gel, indicating the existence of high Mr covalent 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 in M. 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.

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

Abstract The Asian citrus psyllid (Diaphorina citri Kuwayama) is the insect vector of the bacterium Candidatus Liberibacter asiaticus (CLas), the pathogen associated with citrus Huanglongbing (HLB, citrus greening). HLB threatens citrus production worldwide. Suppression or reduction of the insect vector using chemical insecticides has been the primary method to inhibit the spread of citrus greening disease. Accurate structural and functional annotation of the Asian citrus psyllid genome, as well as a clear understanding of the interactions between the insect and CLas, are required for development of new molecular-based HLB control methods. A draft assembly of the D. citri genome has been generated and annotated with automated pipelines. However, knowledge transfer from well-curated reference genomes such as that of Drosophila melanogaster to newly sequenced ones is challenging due to the complexity and diversity of insect genomes. To identify and improve gene models as potential targets for pest control, we manually curated several gene families with a focus on genes that have key functional roles in D. citri biology and CLas interactions. This community effort produced 530 manually curated gene models across developmental, physiological, RNAi regulatory and immunity-related pathways. As previously shown in the pea aphid, RNAi machinery genes putatively involved in the microRNA pathway have been specifically duplicated. A comprehensive transcriptome enabled us to identify a number of gene families that are either missing or misassembled in the draft genome. In order to develop biocuration as a training experience, we included undergraduate and graduate students from multiple institutions, as well as experienced annotators from the insect genomics research community. The resulting gene set (OGS v1.0) combines both automatically predicted and manually curated gene models.

Manduca sexta serpin-12 controls the prophenoloxidase activation system in larval hemolymph

Insect prophenoloxidase activation is coordinated by a serine protease network, which is regulated by serine protease inhibitors of the serpin superfamily. The enzyme system also leads to proteolytic processing of a Spätzle precursor. Binding of Spätzle to a Toll receptor turns on a signaling pathway to induce the synthesis of defense proteins. Previous studies of the tobacco hornworm Manduca sexta have revealed key members of the protease cascade, which generates phenoloxidase for melanogenesis and Spätzle to induce immunity-related genes. Here we provide evidence that M. sexta serpin-12 regulates hemolymph protease-14 (HP14), an initiating protease of the cascade. This inhibitor, unlike the other serpins characterized in M. sexta, has an amino-terminal extension rich in hydrophilic residues and an unusual P1 residue (Leu429) right before the scissile bond cleaved by a target protease. Serpins with similarities to serpin-12, including Drosophila Necrotic, were identified in a wide range of insects including flies, moths, wasps, beetles, and two hemimetabolous species. The serpin-12 mRNA is present at low, constitutive levels in larval fat body and hemocytes and becomes more abundant after an immune challenge. We produced the serpin-12 core domain (serpin-12ΔN) in insect cells and in Escherichia coli and demonstrated its inhibition of human cathepsin G, bovine α-chymotrypsin, and porcine pancreatic elastase. MALDI-TOF analysis of the reaction mixtures confirmed the predicted P1 residue of Leu429. Supplementation of larval plasma samples with the serpin-12ΔN decreased prophenoloxidase activation elicited by microbial cells and reduced the proteolytic activation of the protease precursors of HP6, HP8, PAPs, and other serine protease-related proteins. After incubation of plasma stimulated with peptidoglycan, a 72 kDa protein appeared, which was recognized by polyclonal antibodies against both serpin-12 and HP14, suggesting that a covalent serpin-protease complex formed when serpin-12 inhibited HP14. Together, these data suggest that M. sexta serpin-12 inhibits HP14 to regulate melanization and antimicrobial peptide induction.