Genome-wide annotation and comparative analysis of cuticular protein genes in the noctuid pest Spodoptera litura.

Insect cuticle is considered an adaptable and versatile building material with roles in the construction and function of exoskeleton. Its physical properties are varied, as the biological requirements differ among diverse structures and change during the life cycle of the insect. Although the bulk of cuticle consists basically of cuticular proteins (CPs) associated with chitin, the degree of cuticular sclerotization is an important factor in determining its physical properties. Spodoptera litura, the tobacco cutworm, is an important agricultural pest in Asia. Compared to the domestic silkworm, Bombyx mori, another lepidopteran whose CP genes have been well annotated, S. litura has a shorter life cycle, hides in soil during daytime beginning in the 5th instar and is exposed to soil in the pupal stage without the protection of a cocoon. In order to understand how the CP genes may have been adapted to support the characteristic life style of S. litura, we searched its genome and found 287 putative cuticular proteins that can be classified into 9 CP families (CPR with three groups (RR-1, RR-2, RR-3), CPAP1, CPAP3, CPF, CPFL, CPT, CPG, CPCFC and CPLCA), and a collection of unclassified CPs named CPH. There were also 112 cuticular proteins enriched in Histidine residues with content varying from 6% to 30%, comprising many more His-rich cuticular proteins than B. mori. A phylogenetic analysis between S. litura, M. sexta and B. mori uncovered large expansions of RR-1 and RR-2 CPs, forming large gene clusters in different regions of S. litura chromosome 9. We used RNA-seq analysis to document the expression profiles of CPs in different developmental stages and tissues of S. litura. The comparative genomic analysis of CPs between S. litura and B. mori integrated with the unique behavior and life cycle of the two species offers new insights into their contrasting ecological adaptations.

Genome-wide patterns of copy number variations in Spodoptera litura.

Spodoptera litura is a polyphagous pest and can feed on more than 100 species of plants, causing great damage to agricultural production. The SNP results showed that there were gene exchanges between different regions. To explore the variations of larger segments in S. litura genome, we used genome resequencing samples from 14 regions of China, India, and Japan to study the copy number variations (CNVs). We identified 3976 CNV events and 1581 unique copy number variation regions (CNVRs) occupying the 108.5 Mb genome of S. litura. A total of 5527 genes that overlapped with CNVRs were detected. Selection signal analysis identified 19 shared CNVRs and 105 group-specific CNVRs, whose related genes were involved in various biological processes in S. litura. We constructed the first CNVs map in S. litura genome, and our findings will be valuable for understanding the genomic variations and population differences of S. litura.

Enzymatic characterization of two epsilon-class glutathione S-transferases of Spodoptera litura.

Two cDNAs encoding glutathione S-transferase (GST) of the tobacco cutworm, Spodoptera litura, were cloned by reverse transcriptase-polymerase chain reaction. The deduced amino acid sequences of the resulting clones revealed 32-51% identities to the epsilon-class GSTs from other organisms. The recombinant proteins were functionally overexpressed in Escherichia coli cells in soluble form and were purified to homogeneity. The enzymes were capable of catalyzing the bioconjugation of glutathione with 1-chloro-2,4-dinitrobenzene, 1,2-epoxy-3-(4-nitrophenoxy)-propane, and ethacrynic acid. A competition assay revealed that the GST activity was inhibited by insecticides, suggesting that it could be conducive to insecticide tolerance in the tobacco cutworm.

Genomic adaptation to polyphagy and insecticides in a major East Asian noctuid pest.

The tobacco cutworm, Spodoptera litura, is among the most widespread and destructive agricultural pests, feeding on over 100 crops throughout tropical and subtropical Asia. By genome sequencing, physical mapping and transcriptome analysis, we found that the gene families encoding receptors for bitter or toxic substances and detoxification enzymes, such as cytochrome P450, carboxylesterase and glutathione-S-transferase, were massively expanded in this polyphagous species, enabling its extraordinary ability to detect and detoxify many plant secondary compounds. Larval exposure to insecticidal toxins induced expression of detoxification genes, and knockdown of representative genes using short interfering RNA (siRNA) reduced larval survival, consistent with their contribution to the insect's natural pesticide tolerance. A population genetics study indicated that this species expanded throughout southeast Asia by migrating along a South India-South China-Japan axis, adapting to wide-ranging ecological conditions with diverse host plants and insecticides, surviving and adapting with the aid of its expanded detoxification systems. The findings of this study will enable the development of new pest management strategies for the control of major agricultural pests such as S. litura.

Identification of anhydrobiosis-related genes from an expressed sequence tag database in the cryptobiotic midge Polypedilum vanderplanki (Diptera; Chironomidae).

Some organisms are able to survive the loss of almost all their body water content, entering a latent state known as anhydrobiosis. The sleeping chironomid (Polypedilum vanderplanki) lives in the semi-arid regions of Africa, and its larvae can survive desiccation in an anhydrobiotic form during the dry season. To unveil the molecular mechanisms of this resistance to desiccation, an anhydrobiosis-related Expressed Sequence Tag (EST) database was obtained from the sequences of three cDNA libraries constructed from P. vanderplanki larvae after 0, 12, and 36 h of desiccation. The database contained 15,056 ESTs distributed into 4,807 UniGene clusters. ESTs were classified according to gene ontology categories, and putative expression patterns were deduced for all clusters on the basis of the number of clones in each library; expression patterns were confirmed by real-time PCR for selected genes. Among up-regulated genes, antioxidants, late embryogenesis abundant (LEA) proteins, and heat shock proteins (Hsps) were identified as important groups for anhydrobiosis. Genes related to trehalose metabolism and various transporters were also strongly induced by desiccation. Those results suggest that the oxidative stress response plays a central role in successful anhydrobiosis. Similarly, protein denaturation and aggregation may be prevented by marked up-regulation of Hsps and the anhydrobiosis-specific LEA proteins. A third major feature is the predicted increase in trehalose synthesis and in the expression of various transporter proteins allowing the distribution of trehalose and other solutes to all tissues.

Molecular evolution of lepidopteran silk proteins: insights from the ghost moth, Hepialus californicus.

Silk production has independently evolved in numerous arthropod lineages, such as Lepidoptera, the moths and butterflies. Lepidopteran larvae (caterpillars) synthesize silk proteins in modified salivary glands and spin silk fibers into protective tunnels, escape lines, and pupation cocoons. Molecular sequence data for these proteins are necessary to determine critical features of their function and evolution. To this end, we constructed an expression library from the silk glands of the ghost moth, Hepialus californicus, and characterized light chain fibroin and heavy chain fibroin gene transcripts. The predicted H. californicus silk fibroins share many elements with other lepidopteran and trichopteran fibroins, such as conserved placements of cysteine, aromatic, and polar amino acid residues. Further comparative analyses were performed to determine site-specific signatures of selection and to assess whether fibroin genes are informative as phylogenetic markers. We found that purifying selection has constrained mutation within the fibroins and that light chain fibroin is a promising molecular marker. Thus, by characterizing the H. californicus fibroins, we identified key functional amino acids and gained insight into the evolutionary processes that have shaped these adaptive molecules.

Conservation of silk genes in Trichoptera and Lepidoptera.

Larvae of the sister orders Trichoptera and Lepidoptera are characterized by silk secretion from a pair of labial glands. In both orders the silk filament consists of heavy (H)- and light (L)-chain fibroins and in Lepidoptera it also includes a P25 glycoprotein. The L-fibroin and H-fibroin genes of Rhyacophila obliterata and Hydropsyche angustipennis caddisflies have exon/intron structuring (seven exons in L-fibroin and two in H-fibroin) similar to that in their counterparts in Lepidoptera. Fibroin cDNAs are also known in Limnephilus decipiens, representing the third caddisfly suborder. Amino acid sequences of deduced L-fibroin proteins and of the terminal H-fibroin regions are about 50% identical among the three caddisfly species but their similarity to lepidopteran fibroins is <25%. Positions of some residues are conserved, including cysteines that were shown to link the L-fibroin and H-fibroin by a disulfide bridge in Lepidoptera. The long internal part of H-fibroins is composed of short motifs arranged in species-specific repeats. They are extremely uniform in R. obliterata. Motifs (SX)(n), GGX, and GPGXX occur in both Trichoptera and Lepidoptera. The trichopteran H-fibroins further contain charged amphiphilic motifs but lack the strings of alanines or alanine-glycine dipeptides that are typical lepidopteran motifs. On the other hand, sequences composed of a motif similar to ERIVAPTVITR surrounded by the (SX)(4-6) strings and modifications of the GRRGWGRRG motif occur in Trichoptera and not in Lepidoptera.

SAGE analysis of early oogenesis in the silkworm, Bombyx mori.

To identify genes involved in the differentiation of Bombyx cystoblast, we constructed two 3' long serial analysis of gene expression (Long SAGE) libraries from stage 1-3 or stage 2-3 egg chambers and compared their gene expression profiles. In both libraries, the most frequent tags were derived from the same novel transcript. The transcript does not have any open reading frame capable of encoding a protein with over 100 amino acids in length. RNA blot analysis revealed that this transcript is specifically and abundantly expressed in the Bombyx ovary, mainly the germ line cells in the ovarioles. These results suggest that Bombyx oogenesis may be regulated by a previously unidentified non-coding RNA. Comparison of the gene expression profiles between the stage 1-3 and stage 2-3 egg chamber libraries revealed that 272 tags were significantly more abundant in stage 1-3 egg chambers (p<0.05 and at least two-fold change) than in library 2. Among the differentially expressed transcripts were the sequences that correspond to ATP synthase subunit d (3.1-fold enriched) and ATP synthase coupling factor 6 (9.1-fold enriched), suggesting that they are involved in regulation of cell cycle of cystocytes.

Isolation and Characterization of a Novel Human Radiosusceptibility Gene, NP95.

The murine nuclear protein Np95 has been shown to underlie resistance to ionizing radiation and other DNA insults or replication arrests in embryonic stem (ES) cells. Using the databases for expressed sequenced tags and a two-step PCR procedure, we isolated human NP95, the full-length human homologue of the murine Np95 cDNA, which consists of 4,327 bp with a single open reading frame (ORF) encoding a polypeptide of 793 amino acids and 73.3% homology to Np95. The ORF of human NP95 cDNA is identical to the UHRF1 (ubiquitin-like protein containing PHD and RING domain 1). The NP95 gene, assigned to 19p13.3, consists of 18 exons, spanning 60 kb. Several stable transformants from HEK293 and WI-38 cells that had been transfected with the antisense NP95 cDNA were, like the murine Np95-knockout ES cells, more sensitive to X rays, UV light and hydroxyurea than the corresponding parental cells. In HEK293 cells, the lack of NP95 did not affect the activities of topoisomerase IIalpha, whose expression had been demonstrated to be regulated by the inverted CCAAT box binding protein of 90 kDa (ICBP90) that closely resembles NP95 in amino acid sequence and in cDNA but differs greatly in genomic organization. These findings collectively indicate that the human NP95 gene is the functional orthologue of the murine Np95 gene.

Lepidopteran ortholog of Drosophila breathless is a receptor for the baculovirus fibroblast growth factor.

The Bombyx mori nucleopolyhedrovirus (BmNPV) encodes a gene homologous to the mammalian fibroblast growth factor (FGF) family. We report the cloning of B. mori and Spodoptera frugiperda orthologous genes (Bmbtl and Sfbtl, respectively) of Drosophila melanogaster breathless (btl) encoding a receptor for Branchless/FGF and show that these genes encode the receptor for a baculovirus-encoded FGF (vFGF). Sequence analysis showed that BmBtl is composed of 856 amino acid residues, which potentially encodes a 97.3-kDa polypeptide and shares structural features and sequence similarities with the FGF receptor family. Reverse transcription-PCR experiments showed that Bmbtl was abundantly expressed in the trachea and midgut in B. mori larvae, with moderate expression observed in the hemocytes and the B. mori cultured cell line BmN. We generated Sf-9 cells that stably expressed His-tagged BmBtl. Western blot analysis revealed that BmBtl was an approximately 110-kDa protein. Immunoprecipitation experiments showed that BmNPV vFGF markedly phosphorylated BmBtl in Sf-9 cells. In addition, we found that BmBtl overexpression enhanced the migration activity for BmNPV vFGF. Furthermore, we generated Sf-9 cells in which Sfbtl was knocked down by transfection with double-strand RNA-expressing plasmids. In these cells, cell motility triggered by vFGF was markedly reduced. These results strongly suggest that the Btl orthologs, BmBtl and SfBtl, are the receptors for vFGF, which mediate vFGF-induced host cell chemotaxis.

Characterization and transcriptional profiles of three Spodoptera frugiperda genes encoding cysteine-rich peptides. A new class of defensin-like genes from lepidopteran insects?

The present work describes sequence and transcription of three Spodoptera frugiperda genes encoding 6-cysteine-rich peptides. Sequence alignments indicate that the predicted peptides belong to the insect defensin family, although phylogenetic analyses suggest they form a cluster distinct from that of other neopteran insect defensins. The three genes were identified in a non-immune-challenged Sf9 cells cDNA (DNA complementary to RNA) library (Landais et al., Bioinformatics, in press) and were named spodoptericin, Sf-gallerimycin and Sf-cobatoxin. Spodoptericin is a novel defensin-like gene that appears to be weakly up-regulated following injection of bacteria and fungi. Interestingly, no sequence motif clearly homologous to cis regulatory element involved in the regulation of antimicrobial genes was found. An homologue of the spodoptericin gene was identified in the SilkBase Bombyx mori cDNA library. Sf-gallerimycin is related to the Galleria mellonella gallerimycin gene and is induced after immune challenge by injection of bacteria in the larval fat body as well as in hemocytes. In silico analysis of the sequence upstream from the cDNA reveals the presence of at least one motif homologous to a nuclear factor kappaB (NF-kappaB) binding site. Finally, Sf-cobatoxin is related to the G. mellonella cobatoxin-like gene. Despite high levels of constitutive expression compared to the two previous genes, transcription of Sf-cobatoxin is increased after immune, in particular, bacterial challenge. We therefore confirm that these three genes encode potential candidate molecules involved in S. frugiperda innate humoral response.

Spectrum of Znfn1a1 (Ikaros) inactivation and its association with loss of heterozygosity in radiogenic T-cell lymphomas in susceptible B6C3F1 mice.

Ikaros (now known as Znfn1a1), a Krüppel-type zinc-finger transcription factor that plays a critical role in both lineage commitment and differentiation of lymphoid cells, has recently been shown to function as a tumor suppressor gene. We have previously reported a high frequency of LOH (approximately 50%) at the Znfn1a1 locus in radiation-induced T-cell lymphoma in susceptible B6C3F1 mice. The aim of the present study was to delineate the types of Znfn1a1 inactivation, with special reference to the LOH status, and to determine the relative contribution of each type of Znfn1a1 inactivation in radiation-induced T-cell lymphomas in B6C3F1 mice. We demonstrated that Znfn1a1 was frequently altered (in approximately 50% of T-cell lymphomas), and that its inactivation was caused by a variety of mechanisms, which came under one of the following four categories: (1) null expression (14%); (2) expression of unusual dominant-negative isoforms (11%); (3) amino acid substitutions in the N-terminal zinc-finger domain for DNA binding caused by point mutations (22%); (4) lack of the Znfn1a1 isoform 1 due to the creation of a stop codon by insertion of a dinucleotide in exon 3 (3%). The null expression, amino acid substitutions, and dinucleotide insertion inactivation types were well correlated with LOH at the Znfn1a1 allele (86%) and were consistent with Knudson's two-hit theory. On the other hand, T-cell lymphomas expressing dominant-negative Znfn1a1 isoforms retained both alleles. These results indicate that Znfn1a1 inactivation takes place by a variety of mechanisms in radiation-induced murine T-cell lymphomas and is frequently associated with LOH, this association depending on the type of inactivation.