Genome-wide identification and comparative analysis of lipocalin families in Lepidoptera with an emphasis on Bombyx mori

Lipocalins exhibit functional diversity, including roles in retinol transport, invertebrate cryptic coloration, and stress response. However, genome-wide identification and characterization of lipocalin in the insect lineage have not been thoroughly explored. Here, we found that a lineage-specific expansion of the lipocalin genes in Lepidoptera occurred in large part due to tandem duplication events and several lipocalin genes involving insect coloration were expanded more via tandem duplication in butterflies. A comparative analysis of conserved motifs showed both conservation and divergence of lepidopteran lipocalin family protein structures during evolution. We observe dynamic changes in tissue expression preference of paralogs in Bombyx mori, suggesting differential contribution of paralogs to specific organ functions during evolution. Subcellular localization experiments revealed that lipocalins localize to the cytoplasm, nuclear membrane, or nucleus in BmN cells. Moreover, several lipocalin genes exhibited divergent responses to abiotic and biotic stresses, and 1 lipocalin gene was upregulated by 300 fold in B. mori. These results suggest that lipocalins act as signaling components in defense responses by mediating crosstalk between abiotic and biotic stress responses. This study deepens our understanding of the comprehensive characteristics of lipocalins in insects.     

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.     

Transcriptome analysis of the response of silkworm to drastic changes in ambient temperature

Bombyx mori is a poikilothermic insect and is economically important for silk production. Drastic changes in the ambient temperature have a negative impact on sericulture. However, the reason as to why high temperature is associated with the occurrence of diseases in silkworm and the response of silkworm to low temperature remain unclear and were the focus of the present study. Dazao silkworm exposed to 13 °C (DZ-13), 25 °C (DZ-25), and 37 °C (DZ-37) were used for RNA-seq analysis. There were 478 and 194 upregulated differentially expressed genes (DEGs) in DZ-13 and DZ-37 while 49 and 273 downregulated DEGs in DZ-13 and DZ-37, respectively. Eight DEGs were co-upregulated, in which seven genes were for heat shock proteins (Hsps), implying that Hsps play important roles in the tolerance of silkworm to high and low temperature. Gene ontology analysis revealed that the developmental process was downregulated in DZ-13. All the DEGs in the oxidative phosphorylation and insulin signaling pathways were upregulated in DZ-13. Several cuticular proteins and ATP synthesis-related genes were upregulated in DZ-13, suggesting that thickening of the cuticle and increase in the ATPase expression would help silkworms to protect themselves from low temperature-induced stress. Several immune-related genes, such as BmRel and BmSerpin-2, were downregulated in DZ-37, revealing that the resistance of silkworm is decreased under high temperature shock resulting in susceptibility to pathogens. Thus, the increase in the thermo-tolerance of silkworm should be related to the enhancement in the pathogen resistance.

     

Conserved roles of Osiris genes in insect development,polymorphism and protection

Much of the variation among insects is derived from the different ways that chitin has been moulded to form rigid structures, both internal and external. In this study, we identify a highly conserved expression pattern in an insect‐only gene family, the Osiris genes, that is essential for development, but also plays a significant role in phenotypic plasticity and in immunity/toxicity responses. The majority of Osiris genes exist in a highly syntenic cluster, and the cluster itself appears to have arisen very early in the evolution of insects. We used developmental gene expression in the fruit fly, Drosophila melanogaster, the bumble bee, Bombus terrestris, the harvester ant, Pogonomyrmex barbatus, and the wood ant, Formica exsecta, to compare patterns of Osiris gene expression both during development and between alternate caste phenotypes in the polymorphic social insects. Developmental gene expression of Osiris genes is highly conserved across species and correlated with gene location and evolutionary history. The social insect castes are highly divergent in pupal Osiris gene expression. Sets of co‐expressed genes that include Osiris genes are enriched in gene ontology terms related to chitin/cuticle and peptidase activity. Osiris genes are essential for cuticle formation in both embryos and pupae, and genes co‐expressed with Osiris genes affect wing development. Additionally, Osiris genes and those co‐expressed seem to play a conserved role in insect toxicology defences and digestion. Given their role in development, plasticity, and protection, we propose that the Osiris genes play a central role in insect adaptive evolution.     

Genome-wide comparison of genes involved in the biosynthesis,metabolism, and signaling of juvenile hormone between silkworm and other insects

Juvenile hormone (JH) contributes to the regulation of larval molting and metamorphosis in insects. Herein, we comprehensively identified 55 genes involved in JH biosynthesis, metabolism and signaling in the silkworm (Bombyx mori) as well as 35 in Drosophila melanogaster, 35 in Anopheles gambiae, 36 in Apis mellifera, 47 in Tribolium castaneum, and 44 in Danaus plexippus. Comparative analysis showed that each gene involved in the early steps of the mevalonate (MVA) pathway, in the neuropeptide regulation of JH biosynthesis, or in JH signaling is a single copy in B. mori and other surveyed insects, indicating that these JH-related pathways or steps are likely conserved in all surveyed insects. However, each gene participating in the isoprenoid branch of JH biosynthesis and JH metabolism, together with the FPPS genes for catalyzing the final step of the MVA pathway of JH biosynthesis, exhibited an obvious duplication in Lepidoptera, including B. mori and D. plexippus. Microarray and real-time RT-PCR analysis revealed that different copies of several JH-related genes presented expression changes that correlated with the dynamics of JH titer during larval growth and metamorphosis. Taken together, the findings suggest that duplication-derived copy variation of JH-related genes might be evolutionarily associated with the variation of JH types between Lepidoptera and other insect orders. In conclusion, our results provide useful clues for further functional analysis of JH-related genes in B. mori and other insects.     

The G protein-coupled receptors in the silkworm,Bombyx mori

G protein-coupled receptors (GPCRs) are the largest and most versatile family of transmembrane receptors in the cell, occupying the highest hierarchical positions in the regulation of many physiological processes. Although they have been extensively studied in a number of model insects, there have been few investigations of GPCRs in large Lepidopterans, such as Bombyx mori, an organism that provides a means to perform detailed tissue expression analyses, which may help to characterize GPCRs and their ligands. In addition, B. mori, also known as the silkworm, is an insect of substantial economic importance, due to its use in silk production and traditional medicines. In this work, we computationally identified 90 putative GPCRs in B. mori, 33 of which represent novel proteins. These GPCRs were annotated and compared with their homologs in Drosophila melanogaster and Anopheles gambiae. Phylogenetics analyses of the GPCRs from these three insects showed that GPCRs may easily duplicate or disappear during insect evolution, especially in the neuropeptide and protein hormone receptor subfamily. Interestingly, we observed a decrease in the quantity and diversity of the stress-tolerance gene, Methuselah, in B. mori, which may be related to its long history of domestication. Moreover, the presence of many Bombyx-specific GPCRs suggests that neither Drosophila nor Anopheles is good representatives for the GPCRs in the Class Insecta.     

Characterization of multiple <Emphasis Type="Italic">CYP9A</Emphasis> genes in the silkworm, <Emphasis Type="Italic">Bombyx mori</Emphasis>

Based on the advances in the silkworm genome project, a new genome-wide analysis of cytochrome P450 genes was performed, focusing mainly on gene duplication. All four CYP9A subfamily members from the silkworm, Bombyx mori, were cloned by RT-PCR and designated CYP9A19–CYP9A22 by the P450 Nomenclature Committee. They each contain an open reading frame of 1,593 bp in length and encode a putative polypeptide of 531 amino acids. Both nucleic acid and amino acid sequences share very high identities with one another. The typical motifs of insect cytochrome P450, including the heme-binding region, helix-C, helix-I, helix-K, and PERF, show high sequence conservation among the multiple proteins. Alignment with their cDNA sequences revealed that these paralogues share identical gene structures, each comprising ten exons and nine introns of variable sizes. The locations of their introns (all nine introns follow the GT–AG rule) are absolutely conserved. CYP9A19, CYP9A20, and CYP9A21 form a tandem cluster on chromosome 17, whereas CYP9A22 is separated from the cluster by four tandem alcohol-dehydrogenase-like genes. Their phylogenetic relationships and structural comparisons indicated that these paralogues arose as the results of gene duplication events. RT-PCR detected their mRNAs in different “first line of defense” tissues, as well as in several other organs, suggesting diverse functions. Tissue-selective expression also indicates their functional divergence. The identified CYP9A genes have not yet been found outside the Lepidoptera, and are probably unique to the Lepidoptera. They show high sequence and structural similarities to each other, indicating that the Lepidoptera-specific P450s may be of functional importance. This analysis constitutes the first report of the clustering, spatial organization, and functional divergence of P450 in the silkworm.     

Identification of MicroRNAs in Helicoverpa armigera and Spodoptera litura Based on Deep Sequencing and Homology Analysis

The current identification of microRNAs (miRNAs) in insects is largely dependent on genome sequences. However, the lack of available genome sequences inhibits the identification of miRNAs in various insect species. In this study, we used a miRNA database of the silkworm Bombyx mori as a reference to identify miRNAs in Helicoverpa armigera and Spodoptera litura using deep sequencing and homology analysis. Because all three species belong to the Lepidoptera, the experiment produced reliable results. Our study identified 97 and 91 conserved miRNAs in H. armigera and S. litura, respectively. Using the genome of B. mori and BAC sequences of H. armigera as references, 1 novel miRNA and 8 novel miRNA candidates were identified in H. armigera, and 4 novel miRNA candidates were identified in S. litura. An evolutionary analysis revealed that most of the identified miRNAs were insect-specific, and more than 20 miRNAs were Lepidoptera-specific. The investigation of the expression patterns of miR-2a, miR-34, miR-2796-3p and miR-11 revealed their potential roles in insect development. miRNA target prediction revealed that conserved miRNA target sites exist in various genes in the 3 species. Conserved miRNA target sites for the Hsp90 gene among the 3 species were validated in the mammalian 293T cell line using a dual-luciferase reporter assay. Our study provides a new approach with which to identify miRNAs in insects lacking genome information and contributes to the functional analysis of insect miRNAs.     

Extreme convergence in stick insect evolution: phylogenetic placement of the Lord Howe Island tree lobster

The ‘tree lobsters’ are an enigmatic group of robust, ground-dwelling stick insects (order Phasmatodea) from the subfamily Eurycanthinae, distributed in New Guinea, New Caledonia and associated islands. Its most famous member is the Lord Howe Island stick insect Dryococelus australis (Montrouzier), which was believed to have become extinct but was rediscovered in 2001 and is considered to be one of the rarest insects in the world. To resolve the evolutionary position of Dryococelus, we constructed a phylogeny from approximately 2.4 kb of mitochondrial and nuclear sequence data from representatives of all major phasmatodean lineages. Our data placed Dryococelus and the New Caledonian tree lobsters outside the New Guinean Eurycanthinae as members of an unrelated Australasian stick insect clade, the Lanceocercata. These results suggest a convergent origin of the ‘tree lobster’ body form. Our reanalysis of tree lobster characters provides additional support for our hypothesis of convergent evolution. We conclude that the phenotypic traits leading to the traditional classification are convergent adaptations to ground-living behaviour. Our molecular dating analyses indicate an ancient divergence (more than 22 Myr ago) between Dryococelus and its Australian relatives. Hence, Dryococelus represents a long-standing separate evolutionary lineage within the stick insects and must be regarded as a key taxon to protect with respect to phasmatodean diversity.     

从形态学不同到基因水平差异: 家蚕与果蝇早期胚胎发育比较

昆虫的分化多样性一直是研究的一个热点, 对家蚕胚胎发育调控机制的研究将有助于加深人们对其他昆虫早期胚胎发育机制的了解。文章综述了鳞翅目昆虫家蚕和双翅目昆虫果蝇的早期胚胎发育在形态学及基因水平的研究结果。通过比较, 发现二者的早期胚胎发育存在较大差异。这一结果不但可以为家蚕胚胎发育调控机制的分子水平研究提供重要参考, 也可为昆虫分化多样性的分子机制研究提供线索。     

Effects of transgenic overexpression of diapause hormone and diapause hormone receptor genes on non-diapause silkworm

Diapause is a state of developmental arrest that is most often observed in arthropods, especially insects. The domesticated silkworm, Bombyx mori, is a typical insect that enters diapause at an early embryonic stage. Previous studies have revealed that the diapause hormone (DH) signaling molecules, especially the core members DH and DH receptor 1 (DHR1), are crucial for the determination of embryonic diapause in diapause silkworm strains. However, whether they function in non-diapause silkworm strains remains largely unknown. Here, we generated two transgenic lines overexpressing DH or DHR1 genes in a non-diapause silkworm strain, Nistari. Our results showed that developmental expression patterns of DH and DHR1 are quite similar in transgenic silkworms: both genes are highly expressed in the mid to late stages of pupae and are most highly expressed in day-6 pupae but are expressed at very low levels in other developmental stages. Moreover, the overexpression of DH or DHR1 can affect the expression of diapause-related genes but is not sufficient to induce embryonic diapause in their offspring. This study provides new insights into the function of DH and DHR1 in a non-diapause silkworm strain.     

Constitutive Activation of the Midgut Response to Bacillus thuringiensis in Bt-Resistant Spodoptera exigua

Bacillus thuringiensis is the most effective microbial control agent for controlling numerous species from different insect orders. The main threat for the long term use of B. thuringiensis in pest control is the ability of insects to develop resistance. Thus, the identification of insect genes involved in conferring resistance is of paramount importance. A colony of Spodoptera exigua (Lepidoptera: Noctuidae) was selected for 15 years in the laboratory for resistance to Xentari™, a B. thuringiensis-based insecticide, reaching a final resistance level of greater than 1,000-fold. Around 600 midgut ESTs were analyzed by DNA-macroarray in order to find differences in midgut gene expression between susceptible and resistant insects. Among the differentially expressed genes, repat and arylphorin were identified and their increased expression was correlated with B. thuringiensis resistance. We also found overlap among genes that were constitutively over-expressed in resistant insects with genes that were up-regulated in susceptible insects after exposure to Xentari™, suggesting a permanent activation of the response to Xentari™ in resistant insects. Increased aminopeptidase activity in the lumen of resistant insects in the absence of exposure to Xentari™ corroborated the hypothesis of permanent activation of response genes. Increase in midgut proliferation has been proposed as a mechanism of response to pathogens in the adult from several insect species. Analysis of S. exigua larvae revealed that midgut proliferation was neither increased in resistant insects nor induced by exposure of susceptible larvae to Xentari™, suggesting that mechanisms other than midgut proliferation are involved in the response to B. thuringiensis by S. exigua larvae.     

Ecdysteroid promotes cell cycle progression in the Bombyx wing disc through activation of c-Myc

Developmental switching from growth to metamorphosis in imaginal primordia is an essential process of adult body planning in holometabolous insects. Although it is disciplined by a sequential action of the ecdysteroid, molecular mechanisms linking to cell proliferation are poorly understood. In the present study, we investigated the expression control of cell cycle–related genes by the ecdysteroid using the wing disc of the final-instar larvae of the silkworm, Bombyx mori. We found that the expression level of c-myc was remarkably elevated in the post-feeding cell proliferation phase, which coincided with a small increase in ecdysteroid titer. An in vitro wing disc culture showed that supplementation of the moderate level of the ecdysteroid upregulated c-myc expression within an hour and subsequently increased the expression of cell cycle core regulators, including A-, B-, D-, and E-type cyclin genes, Cdc25 and E2F1. We demonstrated that c-myc upregulation by the ecdysteroid was not inhibited in the presence of a protein synthesis inhibitor, suggesting a possibility that the ecdysteroid directly stimulates c-myc expression. Finally, results from the administration of a c-Myc inhibitor demonstrated that c-Myc plays an essential role in 20E-inducible cell proliferation. These findings suggested a novel pathway for ecdysteroid-inducible cell proliferation in insects, and it is likely to be conserved between insects and mammals in terms of steroid hormone regulation.     

Evolutionary genomics analysis reveals gene expansion and functional diversity of arylalkylamine N-acetyltransferases in the Culicinae subfamily of mosquitoes

Arylalkylamine N-acetyltransferase (aaNAT), considered a potential new insecticide target, catalyzes the acetylation of arylalkylamine substrates such as serotonin and dopamine and, hence, mediates diverse functions in insects. However, the origin of insect aaNATs (iaaNATs) and the evolutionary process that generates multiple aaNATs in mosquitoes remain largely unknown. Here, we have analyzed the genomes of 33 species to explore and expand our understanding of the molecular evolution of this gene family in detail. We show that aaNAT orthologs are present in Bacteria, Cephalochordata, Chondrichthyes, Cnidaria, Crustacea, Mammalia, Placozoa, and Teleoste, as well as those from a number of insects, but are absent in some species of Annelida, Echinozoa, and Mollusca as well as Arachnida. Particularly, more than 10 aaNATs were detected in the Culicinae subfamily of mosquitoes. Molecular evolutionary analysis of aaNAT/aaNAT-like genes in mosquitoes reveals that tandem duplication events led to gene expansion in the Culicinae subfamily of mosquitoes more than 190 million years ago. Further selection analysis demonstrates that mosquito aaNATs evolved under strongly positive pressures that generated functional diversity following gene duplication events. Overall, this study may provide novel insights into the molecular evolution of the aaNAT family in mosquitoes.     

A chromosome-level genome assembly of the beet armyworm Spodoptera exigua

BackgroundThe beet armyworm Spodoptera exigua is a polyphagous caterpillar that causes serious damage to many species of crops and vegetables. To gain insight into how this polyphagous insect differs from less harmful oligophagous species, we generated a chromosome-level assembly and compared it to closely related species with the same or different feeding habits.ResultsBased on Illumina and Pacific Biosciences data and Hi-C technology, 425.6 Mb of genome sequences were anchored and oriented into 31 linkage groups, with an N50 length of 14.8 Mb. A total of 24,649 gene models were predicted, of which 97.4% were identified in the genome assembly. Chemosensory genes are vital for locating food: of the four main families, odorant-binding proteins, chemosensory proteins and olfactory receptors showed little difference, whereas gustatory receptors are greatly expanded in S. exigua. Examination of other polyphagous insects confirmed this difference from oligophagous congeners and further identified the bitter receptor subfamily as being particularly affected.ConclusionOur high-quality genome sequence for beet armyworm identified a key expansion of the bitter gustatory receptor subfamily in this and other pests that differs crucially from more benign relatives and offers insight into the biology and possible future means of control for these economically important insects.