Annotation and expression analysis of cuticular proteins from the tobacco hornworm,Manduca sexta

The insect cuticle is a unique material that covers the exterior of the animal as well as lining the foregut, hindgut, and tracheae. It offers protection from predators and desiccation, defines body shape, and serves as an attachment site for internal organs and muscle. It has demonstrated remarkable variations in hardness, flexibility and elasticity, all the while being light weight, which allows for ease of movement and flight. It is composed primarily of chitin, proteins, catecholamines, and lipids. Proteomic analyses of cuticle from different life stages and species of insects has allowed for a more detailed examination of the protein content and how it relates to cuticle mechanical properties. It is now recognized that several groups of cuticular proteins exist and that they can be classified according to conserved amino acid sequence motifs. We have annotated the genome of the tobacco hornworm, Manduca sexta, for genes that encode putative cuticular proteins that belong to seven different groups: proteins with a Rebers and Riddiford motif (CPR), proteins analogous to peritrophins (CPAP), proteins with a tweedle motif (CPT), proteins with a 44 amino acid motif (CPF), proteins that are CPF-like (CPFL), proteins with an 18 amino acid motif (18 aa), and proteins with two to three copies of a C-X₅-C motif (CPCFC). In total we annotated 248 genes, of which 207 belong to the CPR family, the most for any insect genome annotated to date. Additionally, we discovered new members of the CPAP family and determined that orthologous genes are present in other insects. We established orthology between the M. sexta and Bombyx mori genes and identified duplication events that occurred after separation of the two species. Finally, we utilized 52 RNAseq libraries to ascertain gene expression profiles that revealed commonalities and differences between different tissues and developmental stages.     

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

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

The immune properties of Manduca sexta transferrin

Transferrins are secreted proteins that bind iron. The well-studied transferrins are mammalian serum transferrin, which is involved in iron transport, and mammalian lactoferrin, which functions as an immune protein. Lactoferrin and lactoferrin-derived peptides have bactericidal activity, and the iron-free form of lactoferrin has bacteriostatic activity due to its ability to sequester iron. Insect transferrin is similar in sequence to both serum transferrin and lactoferrin, and its functions are not well-characterized; however, many studies of insect transferrin indicate that it has some type of immune function. The goal of this study was to determine the specific immune functions of transferrin from Manduca sexta (tobacco hornworm). We verified that transferrin expression is upregulated in response to infection in M. sexta larvae and determined that the concentration of transferrin in hemolymph increases from 2 μM to 10 μM following an immune challenge. It is also present in molting fluid and prepupal midgut fluid, two extracellular fluids with immune capabilities. No immune-induced proteolytic cleavage of transferrin in hemolymph was observed; therefore, M. sexta transferrin does not appear to be a source of antimicrobial peptides. Unlike iron-saturated lactoferrin, iron-saturated transferrin had no detectable antibacterial activity. In contrast, 1 μM iron-free transferrin inhibited bacterial growth, and this inhibition was blocked by supplementing the culture medium with 1 μM iron. Our results suggest that M. sexta transferrin does not have bactericidal activity, but that it does have a bacteriostatic function that depends on its iron sequestering ability. This study supports the hypothesis that insect transferrin participates in an iron withholding strategy to protect insects from infectious bacteria.     

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

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

Extensive Conserved Synteny of Genes between the Karyotypes of Manduca sexta and Bombyx mori Revealed by BAC-FISH Mapping

Background

Genome sequencing projects have been completed for several species representing four highly diverged holometabolous insect orders, Diptera, Hymenoptera, Coleoptera, and Lepidoptera. The striking evolutionary diversity of insects argues a need for efficient methods to apply genome information from such models to genetically uncharacterized species. Constructing conserved synteny maps plays a crucial role in this task. Here, we demonstrate the use of fluorescence in situ hybridization with bacterial artificial chromosome probes as a powerful tool for physical mapping of genes and comparative genome analysis in Lepidoptera, which have numerous and morphologically uniform holokinetic chromosomes.

Methodology/Principal Findings

We isolated 214 clones containing 159 orthologs of well conserved single-copy genes of a sequenced lepidopteran model, the silkworm, Bombyx mori, from a BAC library of a sphingid with an unexplored genome, the tobacco hornworm, Manduca sexta. We then constructed a BAC-FISH karyotype identifying all 28 chromosomes of M. sexta by mapping 124 loci using the corresponding BAC clones. BAC probes from three M. sexta chromosomes also generated clear signals on the corresponding chromosomes of the convolvulus hawk moth, Agrius convolvuli, which belongs to the same subfamily, Sphinginae, as M. sexta.

Conclusions/Significance

Comparison of the M. sexta BAC physical map with the linkage map and genome sequence of B. mori pointed to extensive conserved synteny including conserved gene order in most chromosomes. Only a few rearrangements, including three inversions, three translocations, and two fission/fusion events were estimated to have occurred after the divergence of Bombycidae and Sphingidae. These results add to accumulating evidence for the stability of lepidopteran genomes. Generating signals on A. convolvuli chromosomes using heterologous M. sexta probes demonstrated that BAC-FISH with orthologous sequences can be used for karyotyping a wide range of related and genetically uncharacterized species, significantly extending the ability to develop synteny maps for comparative and functional genomics.     

Molecular evolution and expression of the CRAL_TRIO protein family in insects

CRAL_TRIO domain proteins are known to bind small lipophilic molecules such as retinal, inositol and Vitamin E and include such gene family members as PINTA, α-tocopherol transfer (ATT) proteins, retinoid binding proteins, and clavesins. In insects, very little is known about either the molecular evolution of this family of proteins or their ligand specificity. Here we characterize insect CRAL_TRIO domain proteins and present the first insect CRAL_TRIO protein phylogeny constructed by performing reciprocal BLAST searches of the reference genomes of Drosophila melanogaster, Anopheles gambiae, Apis mellifera, Tribolium castaneum, Bombyx mori, Manduca sexta and Danaus plexippus. We find several highly conserved amino acid residues in the CRAL_TRIO domain-containing genes across insects and a gene expansion resulting in more than twice as many gene family members in lepidopterans than in other surveyed insect species, but no lepidopteran homolog of the PINTA gene in Drosophila. In addition, we examined the expression pattern of CRAL_TRIO domain genes in Manduca sexta heads using RNA-Seq data. Of the 42 gene family members found in the M. sexta reference genome, we found 30 expressed in the head tissue with similar expression profiles between males and females. Our results suggest this gene family underwent a large expansion in Lepidoptera, making the lepidopteran CRAL_TRIO domain family distinct from other holometabolous insect lineages.     

Molecular cloning,expression pattern and phylogenetic analysis of myosin light chain 2 gene from Antheraea pernyi: A potential marker for phylogenetic inference

Myosin light chain 2 (MLC-2) gene was isolated and characterized from Antheraea pernyi, a well-known wild silkmoth. The isolated cDNA sequence is 905 bp in length with an open reading frame of 612 bp encoding a polypeptide of 203 amino acids. Semi-quantitative RT-PCR analysis showed that the MLC-2 gene was transcribed during four developmental stages (egg, larva, pupa, and moth), and present in all tissues tested. Alignment analysis revealed that the deduced protein sequence has over 95% identity to myosin light chain 2 of lepidopteran species, and 57–88% identity to other insect species, suggesting that insect MLC-2 proteins are highly conserved throughout evolution. The protein sequence was used to construct phylogenetic trees with other known vertebrate and invertebrate MLC-2 sequences, and the obtained trees demonstrated similar topology with the classical systematics, indicating the potential value of MLC-2 gene in phylogenetic study.     

昆虫基因组及数据库研究进展

基因组序列为昆虫分子生物学研究提供丰富的数据资源,推动系统生物学在古老的昆虫学中蓬勃发展。昆虫基因组学研究已经成为当前的研究热点,目前在NCBI登录注册的昆虫基因组测序计划有494项,其中已提交原始测序数据的昆虫有225种,完成基因组拼接的有215种,具有基因注释的有65种,公开发表的昆虫基因组有43篇。本文综述了测序技术发展的历史及其对昆虫基因组研究的推动作用、昆虫基因组的组装和注释及其存在的问题、昆虫基因组测序进展、昆虫基因组数据库的发展及基因数据挖掘利用的基本思路和对策,以及昆虫基因大数据在害虫防治和资源昆虫利用中的应用前景。     

FISH identification of Helicoverpa armigera and Mamestra brassicae chromosomes by BAC and fosmid probes

Since the Bombyx mori genome sequence was published, conserved synteny between B. mori and some other lepidopteran species has been revealed by either FISH (fluorescence in situ hybridization) with BAC (bacterial artificial chromosome) probes or linkage analysis. However, no species belonging to the Noctuidae, the largest lepidopteran family which includes serious polyphagous pests, has been analyzed so far with respect to genome-wide conserved synteny and gene order. For that purpose, we selected the noctuid species Helicoverpa armigera and Mamestra brassicae, both with n = 31 chromosomes. Gene-defined fosmid clones from M. brassicae and BAC clones from a closely related species of H. armigera, Heliothis virescens, were used for a FISH analysis on pachytene chromosomes. We recognized all H. armigera chromosomes from specific cross-hybridization signals of 146 BAC probes. With 100 fosmid clones we identified and characterized all 31 bivalents of M. brassicae. Synteny and gene order were well conserved between the two noctuid species. The comparison with the model species B. mori (n = 28) showed the same phenomenon for 25 of the 28 chromosomes. Three chromosomes (#11, #23 and #24) had two counterparts each in H. armigera and M. brassicae. Since n = 31 is the modal chromosome number in Lepidoptera, the noctuid chromosomes probably represent an ancestral genome organization of Lepidoptera. This is the first identification of a full karyotype in Lepidoptera by means of BAC cross-hybridization between species. The technique shows the potential to expand the range of analyzed species efficiently.