Transcriptomic analysis of developmental features of Bombyx mori wing disc during metamorphosis

Non-gonadal somatic piRNA pathways ensure sexual differentiation, larval growth, and wing development in silkworms

PIWI-interacting RNAs (piRNAs) guide PIWI proteins to target transposons in germline cells, thereby suppressing transposon activity to preserve genome integrity in metazoans' gonadal tissues. Piwi, one of three Drosophila PIWI proteins, is expressed in the nucleus and suppresses transposon activity by forming heterochromatin in an RNA cleavage-independent manner. Recently, Piwi was reported to control cell metabolism in Drosophila fat body, providing an example of piRNAs acting in non-gonadal somatic tissues. However, mutant flies of the other two PIWI proteins, Aubergine (Aub) and Argonaute3 (Ago3), show no apparent phenotype except for infertility, blurring the importance of the piRNA pathway in non-gonadal somatic tissues. The silkworm, Bombyx mori, possesses two PIWI proteins, Siwi (Aub homolog) and BmAgo3 (Ago3 homolog), whereas B. mori does not have a Piwi homolog. Siwi and BmAgo3 are mainly expressed in gonadal tissues and play a role in repressing transposon activity by cleaving transposon RNA in the cytoplasm. Here, we generated Siwi and BmAgo3 loss-of-function mutants of B. mori and found that they both showed delayed larval growth and failed to become adult moths. They also exhibited defects in wing development and sexual differentiation. Transcriptome analysis revealed that loss of somatic piRNA biogenesis pathways results in abnormal expression of not only transposons but also host genes, presumably causing severe growth defects. Our results highlight the roles of non-gonadal somatic piRNAs in B. mori development.

Cloak Scavenges the Reactive Oxygen Species around the Larvae of Drino inconspicuoides (Diptera: Tachinidae)

Drino inconspicuoides (Diptera: Tachinidae) is an endoparasitoid that develops inside the lepidopteran host. When the larva of D. inconspicuoides penetrates into the host, Mythimna separata (Lepidoptera: Noctuidae), the larva creates a cap-like structure, called the funnel, by using host hemocytes, forming a respiratory attachment to permit efficient respiration. A newly described cloudy and cottony structure, called the "cloak", is formed outside the funnel within 24 h of parasitism. The cloak contains the host fat body and hemocytes. In this study, we aimed to clarify the function of the cloak, which has to date remained unknown. We hypothesized that the funnel generates reactive oxygen species (ROS) through melanization, and that the cloak detoxifies them. We confirmed that the black pigments of the funnel were caused by melanization, which inevitably generates ROS that are potentially harmful to the D. inconspicuoides larva inside the funnel. The cloak showed high activities of antioxidant enzymes, including superoxide dismutase, glutathione peroxidase, and catalase. These results suggest that the cloak scavenged the ROS from the melanized funnel through the diversion of antioxidant enzymes in the fat body, thereby protecting the D. inconspicuoides larva from oxidative damage.

Comparative transcriptomics of the irradiated melon fly (Zeugodacus cucurbitae) reveal key developmental genes

Irradiation can be used as an insect pest management technique to reduce post-harvest yield losses. It causes major physiological changes, impairing insect development and leading to mortality. This technique is used to control the melon fly Zeugodacus cucurbitae, a major pest of Cucurbitaceae in Asia. Here, we applied irradiation to melon fly eggs, and the larvae emerged from irradiated eggs were used to conduct comparative transcriptomics and thereby identify key genes involved in the development and survival. We found 561 upregulated and 532 downregulated genes in irradiated flies compared to non-irradiated flies. We also observed abnormal small-body phenotypes in irradiated flies. By screening the 532 downregulated genes, we selected eight candidate genes putatively involved in development based in described functions in public databases and in the literature. We first established the expression profile of each candidate gene. Using RNA interference (RNAi), we individually knocked down each gene in third instar larvae and measured the effects on development. The knockdown of ImpE2 ecdysone-inducible gene controlling life stage transitions-led to major body size reductions in both pupae and adults. The knockdown of the tyrosine-protein kinase-like tok (Tpk-tok) caused severe body damage to larvae, characterized by swollen and black body parts. Adults subject to knockdown of the eclosion hormone (Eh_1) failed to shed their old cuticle which remained attached to their bodies. However, no obvious developmental defects were observed following the knockdown of the heat shock protein 67B1-like (Hsp67), the insulin receptor (Insr), the serine/threonine-protein kinase Nek4 (Nek4), the tyrosine-protein kinase transmembrane receptor Ror (Ror_1) and the probable insulin-like peptide 1 (Insp_1). We argue that irradiation can be successfully used not only as a pest management technique but also for the screening of essential developmental genes in insects via comparative transcriptomics. Our results demonstrate that ImpE2 and Eh_1 are essential for the development of melon fly and could therefore be promising candidates for the development of RNAi-based pest control strategies.

RNA Sequencing Reveals the Potential Adaptation Mechanism to Different Hosts of Grapholita molesta

Grapholita molesta is an important fruit tree worldwide pest which feeds on hosts extensively and does serious harm. In this paper, the growth and development parameters and protease activities of G. molesta fed on different hosts were compared. Using Illumina RNA sequencing technology, 18 midgut samples from five different hosts (apple, pear, plum, peach and peach shoots) and artificial diet were sequenced and compared with the reference genome, resulting in 15269 genes and 2785 predicted new genes. From 15 comparative combinations, DEGs were found from 286 to 4187 in each group, with up-regulated genes from 107 to 2395 and down-regulated genes from 83 to 2665. KEGG pathway analysis showed that DEGs were associated with amino acid metabolism, starch and sucrose metabolism, carbohydrate metabolism, and hydrolase activity. A total of 31 co-expression gene modules of different hosts were identified by WGCNA. qRT-PCR showed that the expression pattern of the trypsin gene was consistent with RNA sequencing. In this study, growth and development parameters, protease activity, DEGs, enrichment analysis and qRT-PCR were combined to reveal the adaptation process to different hosts of G. molesta in many aspects. The results of this study provide a basis for further exploration of the molecular mechanism of host adaptation of G. molesta.

DNA G-quadruplex structure participates in regulation of lipid metabolism through acyl-CoA binding protein

G-quadruplex structure (G4) is a type of DNA secondary structure that widely exists in the genomes of many organisms. G4s are believed to participate in multiple biological processes. Acyl-CoA binding protein (ACBP), a ubiquitously expressed and highly conserved protein in eukaryotic cells, plays important roles in lipid metabolism by transporting and protecting acyl-CoA esters. Here, we report the functional identification of a G4 in the promoter of the ACBP gene in silkworm and human cancer cells. We found that G4 exists as a conserved element in the promoters of ACBP genes in invertebrates and vertebrates. The BmACBP G4 bound with G4-binding protein LARK regulated BmACBP transcription, which was blocked by the G4 stabilizer pyridostatin (PDS) and G4 antisense oligonucleotides. PDS treatment with fifth instar silkworm larvae decreased the BmACBP expression and triacylglycerides (TAG) level, resulting in reductions in fat body mass, body size and weight and growth and metamorphic rates. PDS treatment and knocking out of the HsACBP G4 in human hepatic adenocarcinoma HepG2 cells inhibited the expression of HsACBP and decreased the TAG level and cell proliferation. Altogether, our findings suggest that G4 of the ACBP genes is involved in regulation of lipid metabolism processes in invertebrates and vertebrates.

U1 small nuclear ribonucleoprotein is essential for early larval development in silkworm, Bombyx mori

U1 small nuclear ribonucleoproteins (U1 snRNP) associates with 5' splice sites in the form of ribonucleoprotein particles and is highly conserved among species. The physiological functions of U1 snRNP in a lepidopteran model insect Bombyx mori is mostly unknown. Here, we showed that U1 snRNP plays an important role in the development of silkworm. Knockout of U1 snRNP in silkworm showed either delayed or stationary 1st instar larva development compared with the wild-type group. U1 snRNP deletion mutants exhibited abnormal cellular phenotypes with enlarged cell nucleus, scant cytoplasm and enlarged nuclei. RNA-seq analysis revealed that genes involved in metabolic pathway, biosynthesis of secondary metabolites and steroid hormone biosynthesis were significantly affected by U1 snRNP depletion. Taken together, our study suggests that U1 snRNP homeostasis plays an important role in silkworm development.

Genome-Wide Identification of Aldehyde Oxidase Genes in Moths and Butterflies Suggests New Insights Into Their Function as Odorant-Degrading Enzymes

Aldehyde oxidases (AOXs) are common detoxifying enzymes in several organisms. In insects, AOXs act in xenobiotic metabolism and as odorant-degrading enzymes (ODEs). These last appear as crucial enzymes in the life cycle of insects, helping to reset their olfactory system, particularly in lepidopterans, which fulfill important ecological roles (e.g., pollination or destructive life cycles). A comprehensive understanding of their olfactory system has provided opportunities to study key chemosensory proteins. However, no significant advance has been made around lepidopteran AOXs research, and even less around butterflies, a recently evolved lineage. In this study we have identified novel AOX gene families in moths and butterflies in order to understand their role as ODEs. Eighteen genomes from both moths and butterflies were used for phylogenetics, molecular evolution and sequence analyses. We identified 164 AOXs, from which 91 are new. Their phylogeny showed two main clades that are potentially related to odorant-degrading function, where both moths and butterflies have AOXs. A first ODE-related clade seems to have a non-ditrysian origin, likely related to plant volatiles. A second ODE-related clade could be more pheromone-biased. Molecular evolution analysis suggests a slight purifying selection process, though a number of sites appeared under positive selection. ODE-related AOXs have changed a phenylalanine residue by proline in the active site. Finally, this study could serve as a reference for further evolutionary and functional studies around Lepidopteran AOXs.

The Hox gene Antennapedia is essential for wing development in insects

A long-standing view in the field of evo-devo is that insect forewings develop without any Hox gene input. The Hox gene Antennapedia (Antp), despite being expressed in the thoracic segments of insects, has no effect on wing development. This view has been obtained from studies in two main model species: Drosophila and Tribolium. Here, we show that partial loss of function of Antp resulted in reduced and malformed adult wings in Bombyx, Drosophila and Tribolium. Antp mediates wing growth in Bombyx by directly regulating the ecdysteriod biosynthesis enzyme gene (shade) in the wing tissue, which leads to local production of the growth hormone 20-hydroxyecdysone. Additional targets of Antp are wing cuticular protein genes CPG24, CPH28 and CPG9, which are essential for wing development. We propose, therefore, that insect wing development occurs in an Antp-dependent manner. This article has an associated 'The people behind the papers' interview.

Transcriptome dynamics during metamorphosis of imaginal discs into wings and thoracic dorsum in Apis mellifera castes

BACKGROUND

Much of the complex anatomy of a holometabolous insect is built from disc-shaped epithelial structures found inside the larva, i.e., the imaginal discs, which undergo a rapid differentiation during metamorphosis. Imaginal discs-derived structures, like wings, are built through the action of genes under precise regulation.

RESULTS

We analyzed 30 honeybee transcriptomes in the search for the gene expression needed for wings and thoracic dorsum construction from the larval wing discs primordia. Analyses were carried out before, during, and after the metamorphic molt and using worker and queen castes. Our RNA-seq libraries revealed 13,202 genes, representing 86.2% of the honeybee annotated genes. Gene Ontology analysis revealed functional terms that were caste-specific or shared by workers and queens. Genes expressed in wing discs and descendant structures showed differential expression profiles dynamics in premetamorphic, metamorphic and postmetamorphic developmental phases, and also between castes. At the metamorphic molt, when ecdysteroids peak, the wing buds of workers showed maximal gene upregulation comparatively to queens, thus underscoring differences in gene expression between castes at the height of the larval-pupal transition. Analysis of small RNA libraries of wing buds allowed us to build miRNA-mRNA interaction networks to predict the regulation of genes expressed during wing discs development.

CONCLUSION

Together, these data reveal gene expression dynamics leading to wings and thoracic dorsum formation from the wing discs, besides highlighting caste-specific differences during wing discs metamorphosis.

Genome and transcriptome analysis of the beet armyworm Spodoptera exigua reveals targets for pest control

The genus Spodoptera (Lepidoptera: Noctuidae) includes some of the most infamous insect pests of cultivated plants including Spodoptera frugiperda, Spodoptera litura, and Spodoptera exigua. To effectively develop targeted pest control strategies for diverse Spodoptera species, genomic resources are highly desired. To this aim, we provide the genome assembly and developmental transcriptome comprising all major life stages of S. exigua, the beet armyworm. Spodoptera exigua is a polyphagous herbivore that can feed on > 130 host plants, including several economically important crops. The 419 Mb beet armyworm genome was sequenced from a female S. exigua pupa. Using a hybrid genome sequencing approach (Nanopore long-read data and Illumina short read), a high-quality genome assembly was achieved (N50 = 1.1 Mb). An official gene set (18,477 transcripts) was generated by automatic annotation and by using transcriptomic RNA-seq datasets of 18 S. exigua samples as supporting evidence. In-depth analyses of developmental stage-specific expression combined with gene tree analyses of identified homologous genes across Lepidoptera genomes revealed four potential genes of interest (three of them Spodoptera-specific) upregulated during first- and third-instar larval stages for targeted pest-outbreak management. The beet armyworm genome sequence and developmental transcriptome covering all major developmental stages provide critical insights into the biology of this devastating polyphagous insect pest species worldwide. In addition, comparative genomic analyses across Lepidoptera significantly advance our knowledge to further control other invasive Spodoptera species and reveals potential lineage-specific target genes for pest control strategies.

RNA-seq analysis of gene expression changes in cuticles during the larval-pupal metamorphosis of Plutella xylostella

The diamondback moth, Plutella xylostella (Lepidoptera: Plutellidae) is a holometabolous insect that its cuticles must undergo the significant changes during the larval-pupal metamorphosis development. To elucidate these changes at molecular levels, RNA-seq analysis of cuticles from LLS (later fourth instar larval stage), PPS (prepupal stage) and PS (pupal stage) were performed in P. xylostella. In this paper, a total of 17,710 transcripts were obtained in the larval-pupal transition of P. xylostella, and out of which 2293 (881 up-regulated and 1412 down-regulated) and 2989 transcripts (2062 up-regulated and 927 down-regulated) were identified to be differentially expressed between LLS and PPS, as well as PPS and PS, respectively. The further GO and KEGG analysis of differentially expressed genes (DEGs) revealed that the 'structural constituent of cuticle', 'chitin metabolic process', 'chitin binding', 'tyrosine metabolism' and 'insect hormone biosynthesis' pathways were significantly enriched, indicating these pathways might be involved in the process of larval pupation in P. xylostella. Then, we found some genes that encoded cuticular proteins, chitinolytic enzymes, chitin synthesis enzymes, and cuticle tanning proteins changed their expression levels remarkably, indicating these genes might play important roles in the restruction (degradation and biosynthesis) of insect cuticles during the larval metamorphosis. Additionally, the significant changes in the mRNA levels of 20-hydroxyecdysone (20E) and juvenile hormone (JH) related genes suggested their crucial roles in regulating cuticle remodeling during the larval metamorphosis of P. xylostella. In conclusion, the present study provide us the comprehensive gene expression profiles to explore the molecular mechanisms of cuticle metamorphosis in P. xylostella, which laid a molecular basis to study roles of specific pathways and genes in insect development.

Molecular characterization and function of chitin deacetylase-like from the Chinese mitten crab, Eriocheir sinensis

Chitin deacetylases are essential enzymes in the chitin-modifying process and play vital roles in arthropod molting. In this study, we identified and characterized a chitin deacetylase-like (EsCDA-l) gene in the Chinese mitten crab, Eriocheir sinensis. The open reading frame of EsCDA-l was 2555 bp and encoded 554 amino acid residues that contained typical domain structure of carbohydrate esterase family 4. Phylogenetic analysis reveal that EsCDA-l belongs to the group I chitin deacetylase family. Quantitative real-time PCR analyses showed that EsCDA-l was highly expressed in exoskeletal tissues and megalopa stages. During the molting cycle, EsCDA-l was up-regulated periodically in the post-molt stage. Knockdown of EsCDA-l resulted in the abnormal ultrastructure of cuticle, prevented molting to high mortality suggesting EsCDA-l is indispensable for molting. The characterization and function analysis of the EsCDA-l should provide useful reference for further research on the utility of key genes involved in the chitin metabolic pathway in the molting process of the Chinese mitten crab as well as other crustaceans.

Multilayer fibroin/chitosan oligosaccharide lactate and pullulan immunomodulatory patch for treatment of hernia and prevention of intraperitoneal adhesion

This study aims to develop a novel intraperitoneal two- or three-layered patch with immunomodulatory property for treatment of hernia, regeneration of abdominal wall and prevention of intraperitoneal adhesions. Polypropylene (PP) mesh, middle layer, was intended to provide mechanical support whereas pullulan (PUL) hydrogel coating layer was designed to prevent intraperitoneal adhesions. Fibroin/chitosan oligosaccharide lactate (F/COS) layer electrospun on one side of pullulan was chosen for immunomodulation and abdominal wall regeneration. Physical and mechanical properties and regenerative capacity of intraperitoneal patches were determined. Immunomodulatory property of electrospun layer and whole patch was studied by determining nitric oxide amount produced by RAW 264.7 macrophages. 25 % (w/v) PUL hydrogel and F/COS with 90:10 (w/w) ratio yielded optimal results. Here, we report that fabricated intraperitoneal patches successfully prevented cell adhesion on one side and increased cell viability and proliferation on other side, along with immunomodulation, in vitro.

Glycosyltransferases: the multifaceted enzymatic regulator in insects

Glycosyltransferases (GTs) catalyse the reaction of glyco-conjugation of various biomolecules by transferring the saccharide moieties from an activated nucleotide sugar to nucleophilic glycosyl acceptor. In insects, GTs show diverse temporal and site-specific expression patterns and thus play significant roles in forming the complex biomolecular structures that are necessary for insect survival, growth and development. Several insects exhibit GT-mediated detoxification as a key defence strategy against plant allelochemicals and xenobiotic compounds, as well as a mechanism for pesticide cross-resistance. Also, these enzymes act as crucial effectors and modulators in various developmental processes of insects such as eye development, UV shielding, cuticle formation, epithelial development and other specialized functions. Furthermore, many of the known insect GTs have been shown to play a fundamental role in other physiological processes like body pigmentation, cuticular tanning, chemosensation and stress response. This review provides a detailed overview of the multifaceted functionality of insect GTs and summarizes numerous case studies associated with it.

DIA-based proteome reveals the involvement of cuticular proteins and lipids in the wing structure construction in the silkworm

Wing discs of Bombyx mori (B. mori) are transformed into wings during metamorphosis via dramatic morphological and structural changes. Mutations in genes related to the wings cause the adults to have altered wing shapes or abnormal wing colour. At present, there are more than 20 wing mutants recorded in the silkworm. However, the key factors that influence B. mori wing development are still unclear. Here, we used the strains +^Wes/+^Wes and Wes/+^Wes that are typical for the normal wing and shriveled wing phenotypes, respectively, to identify differentially expressed proteins by label-free data-independent acquisition (DIA). Ten enriched GO terms and 9 KEGG pathways were identified based on the 3993 proteins in the wings. Among the identified and quantified proteins, 370 differentially expressed proteins (DEPs) were detected (P-value <0.01, |log2FC| > 0.58). Mapping of the DEPs to the reference canonical pathways in KEGG showed that the top 20% of the pathways were related to fatty acid, cutin, suberin and wax biosynthesis, protein processing in endoplasmic reticulum, protein export, etc. Of the 370 DEPs, 238 were down-regulated, and 132 were up-regulated of Wes/+^Wes compared with +^Wes/+^Wes. Numerous cuticular proteins were down-regulated, and fatty metabolism enzymes were up-regulated, in Wes/+^Wes compared with +^Wes/+^Wes. SIGNIFICANCE: The comparative analysis of proteomes suggested that cuticular proteins and fatty metabolism enzymes are the main abnormally expressed proteins in the pupal wings of Wes/+^Wes, leading to curly and shrunken wings after moth transformation. Our results also identify the substances affecting the development of silkworm wings from the perspective of proteins. The information from this study is important for further research on the molecular mechanisms of wing development in lepidopteran insects, and these differentially expressed genes may be targets for Lepidoptera pest control.

Mechanisms underlying the control of dynamic regulatory element activity and chromatin accessibility during metamorphosis

Cis-regulatory modules of metazoan genomes determine the when and where of gene expression during development. Here we discuss insights into the genetic and molecular mechanisms behind cis-regulatory module usage that have come from recent application of genomics assays to insect metamorphosis. Assays including FAIRE-seq, ATAC-seq, and CUT&RUN indicate that sequential changes in chromatin accessibility play a key role in mediating stage-specific cis-regulatory module activity and gene expression. We review the current understanding of what controls precisely coordinated changes in chromatin accessibility during metamorphosis and describe evidence that points to systemic hormone signaling as a primary signal to trigger genome-wide shifts in accessibility patterns and cis-regulatory module usage.

The Developmental Transcriptome of Bagworm, Metisa plana (Lepidoptera: Psychidae) and Insights into Chitin Biosynthesis Genes

Bagworm, Metisa plana (Lepidoptera: Psychidae) is a ubiquitous insect pest in the oil palm plantations. M. plana infestation could reduce the oil palm productivity by 40% if it remains untreated over two consecutive years. Despite the urgency to tackle this issue, the genome and transcriptome of M. plana have not yet been fully elucidated. Here, we report a comprehensive transcriptome dataset from four different developmental stages of M. plana, comprising of egg, third instar larva, pupa and female adult. The de novo transcriptome assembly of the raw data had produced a total of 193,686 transcripts, which were then annotated against UniProt, NCBI non-redundant (NR) database, Gene Ontology, Cluster of Orthologous Group, and Kyoto Encyclopedia of Genes and Genomes databases. From this, 46,534 transcripts were annotated and mapped to 146 known metabolic or signalling KEGG pathways. The paper further identified 41 differentially expressed transcripts encoding seven genes in the chitin biosynthesis pathways, and their expressions across each developmental stage were further analysed. The genetic diversity of M. plana was profiled whereby there were 21,516 microsatellite sequences and 379,895 SNPs loci found in the transcriptome of M. plana. These datasets add valuable transcriptomic resources for further study of developmental gene expression, transcriptional regulations and functional gene activities involved in the development of M. plana. Identification of regulatory genes in the chitin biosynthesis pathway may also help in developing an RNAi-mediated pest control management by targeting certain pathways, and functional studies of the genes in M. plana.

Molecular identification, expression and function analysis of peroxidasin in Chilo suppressalis

Peroxidasin plays a unique role in the formation and stability of extracellular matrix (ECM) in the animal kingdom; however, it was only characterized in Diptera, not in other insect orders. In this study peroxidasin (CsPxd) was first identified and characterized from Chilo suppressalis, a lepidopteran pest. CsPxd complementary DNA with a 4080 bp open reading frame encodes a peptide of 1359 amino acids; the derived amino acid sequence of CsPxd harbors the typical structural characteristics of peroxidasin family in heme-peroxidase superfamily, including the signal peptide at N-terminal, leucine-rich repeat domain, Ig-loop motifs and peroxidase domain, signifying the extracellular location of protein and the involvement in ECM formation. Eukaryotic expression reveals CsPxd protein displays peroxidase activity on H2 O2 , justifying the membership of peroxidase. Phyletic analysis shows the monophyletic evolution pattern of peroxidasin in insect phyle, and moreover only one peroxidasin is present in each species of insects, suggesting its evolutionary conservation on function. Peroxidasin messenger RNA is mainly expressed in egg and the final instar larvae stage. Injection of peroxidasin double-stranded RNA into the final instar larvae impacts the cuticle sclerotization during the metamorphosis from larvae to pupa, and eventually lead to lethality of larvae and pupa. These results suggest the presence of collagen crosslink in chorion and cuticle of insects, and indicate peroxidasin plays a role in the development of chorion and cuticle; furthermore peroxidasin might be the one of potential target genes for pest control using RNA interference.

Transcriptome and gene expression analysis of Rhynchophorus ferrugineus (Coleoptera: Curculionidae) during developmental stages

Background

Red palm weevil, Rhynchophorus ferrugineus Olivier, is one of the most destructive pests harming palm trees. However, genomic resources for R. ferrugineus are still lacking, limiting the ability to discover molecular and genetic means of pest control.

Methods

In this study, PacBio Iso-Seq and Illumina RNA-seq were used to generate transcriptome from three developmental stages of R. ferrugineus (pupa, 7th-instar larva, adult) to increase the understanding of the life cycle and molecular characteristics of the pest.

Results

Sequencing generated 625,983,256 clean reads, from which 63,801 full-length transcripts were assembled with N50 of 3,547 bp. Expression analyses revealed 8,583 differentially expressed genes (DEGs). Moreover, gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that these DEGs were mainly related to the peroxisome pathway which associated with metabolic pathways, material transportation and organ tissue formation. In summary, this work provides a valuable basis for further research on the growth and development, gene expression and gene prediction, and pest control of R. ferrugineus.

Pigment Dispersing Factors and Their Cognate Receptors in a Crustacean Model, With New Insights Into Distinct Neurons and Their Functions

Pigment dispersing factors (PDFs, or PDHs in crustaceans) form a structurally related group of neuropeptides found throughout the Ecdysozoa and were first discovered as pigmentary effector hormones in crustaceans. In insects PDFs fulfill crucial neuromodulatory roles, most notably as output regulators of the circadian system, underscoring their central position in physiological and behavioral organization of arthropods. Intriguingly, decapod crustaceans express multiple isoforms of PDH originating from separate genes, yet their differential functions are still to be determined. Here, we functionally define two PDH receptors in the crab Carcinus maenas and show them to be selectively activated by four PDH isoforms: PDHR 43673 was activated by PDH-1 and PDH-2 at low nanomolar doses whilst PDHR 41189 was activated by PDH-3 and an extended 20 residue e-PDH. Detailed examination of the anatomical distribution of all four peptides and their cognate receptors indicate that they likely perform different functions as secreted hormones and/or neuromodulators, with PDH-1 and its receptor 43,673 implicated in an authentic hormonal axis. PDH-2, PDH-3, and e-PDH were limited to non-neurohemal interneuronal sites in the CNS; PDHR 41189 was largely restricted to the nervous system suggesting a neuromodulatory function. Notably PDH-3 and e-PDH were without chromatophore dispersing activity. This is the first report which functionally defines a PDHR in an endocrine system in a crustacean and to indicate this and other putative roles of this physiologically pivotal peptide group in these organisms. Thus, our findings present opportunities to further examine the endocrine and circadian machinery in this important arthropod phylum.

Alcohol dehydrogenase 5 of Helicoverpa armigera interacts with the CYP6B6 promoter in response to 2-tridecanone

Alcohol dehydrogenase 5 (ADH5) is a member of medium-chain dehydrogenase/reductase family and takes part in cellular formaldehyde and S-nitrosoglutathione metabolic network. 2-tridecanone (2-TD) is a toxic compound in many Solanaceae crops to defend against a variety of herbivory insects. In the broader context of insect development and pest control strategies, this study investigates how a new ADH5 from Helicoverpa armigera (HaADH5) regulates the expression of CYP6B6, a gene involved in molting and metamorphosis, in response to 2-TD treatment. Cloning of the HaADH5 complementary DNA sequence revealed that its 1002 bp open reading frame encodes 334 amino acids with a predicted molecular weight of 36.5 kD. HaADH5 protein was purified in the Escherichia coli Transetta (pET32a-HaADH5) strain using a prokaryotic expression system. The ability of HaADH5 protein to interact with the 2-TD responsive region within the promoter of CYP6B6 was confirmed by an in vitro electrophoretic mobility shift assay and transcription activity validation in yeast. Finally, the expression levels of both HaADH5 and CYP6B6 were found to be significantly decreased in the midgut of 6th instar larvae after 48 h of treatment with 10 mg/g 2-TD artificial diet. These results indicate that upon 2-TD treatment of cotton bollworm, HaADH5 regulates the expression of CYP6B6 by interacting with its promoter. As HaADH5 regulation of CYP6B6 expression may contribute to the larval xenobiotic detoxification, molting and metamorphosis, HaADH5 is a candidate target for controlling the growth and development of cotton bollworm.

RNA interference of eclosion hormone gene reveals its roles in the control of ecdysis behavior in Heortia vitessoides Moore

Eclosion hormone (EH) is an important neuropeptide that regulates growth and development. This study predicted the EH gene (HvEH) of Heortia vitessoides Moore (Lepidoptera: Crambidae) from the transcriptome database and its expression patterns were determined using quantitative real-time polymerase chain reaction. HvEH was expressed in all developmental stages and especially in the head area. RNA interference-mediated silencing of HvEH (2 μg/individual) with double-stranded HvEH RNA (dsHvEH) was achieved within 48 hr. Abnormal phenotypes appeared in the pupa and adult stages. dsHvEH injection suppressed pupation and eclosion rates. HvEH expression increased upon treatment with 20-hydroxyecdysone but decreased at extreme temperatures. These results suggest that HvEH plays an essential role in ecdysis and wing formation in H. vitessoides.

Potential targets for controlling Bactrocera dorsalis using cuticle- and hormone-related genes revealed by a developmental transcriptome analysis

BACKGROUND

The oriental fruit fly, Bactrocera dorsalis (Hendel), is an important agricultural pest and has developed resistance to many insecticides. To investigate vital genes participating in metamorphosis for development of additional control strategies, a comprehensive transcriptome analysis covering ten developmental stages of B. dorsalis was performed.

RESULTS

There were 2132, 952, 1062, 2301 and 1333 differentially expressed genes identified during hatching, 1st-instar larval molting, 2nd-instar larval molting, pupariation and emergence, respectively. Further expression analyses indicated that genes in hormone- (20-hydroxyecdysone and juvenile hormone) and cuticle- (chitin and cuticle protein) related pathways were essential for metamorphosis in B. dorsalis. Among chitinase (Cht) genes, BdCht-5, -8 and -10 were differentially expressed during larval-larval, larval-pupal and pupal-adult moltings. However, BdCht7 was differentially expressed during egg-larval and larval-larval moltings. Knockdown of BdCht7 at the 1st-instar larval stage disrupted normal development of larvae and was lethal to B. dorsalis. Among cuticle protein (CP) genes, 15 genes (BdCPLCG-1, BdCPLCP-2, BdCPAP1-B2, BdRR1-21, BdRR1-31, BdRR2-15, BdRR2-26, BdRR2-30, BdRR2-32, BdTweedle-9, BdTweedle-24, BdRR2-10, BdCPAP3-C1, BdRR1-34 and BdRR1-41) were differentially expressed during four of five types of moltings. Among hormone-relative genes, BdJHBP-4, -9 and -13 were differentially expressed during all five types of moltings, whereas BdJHBP-5, -12 and BdHR4 were differentially expressed during four of five types of moltings.

CONCLUSION

This study reveals critical genes involved in development and metamorphosis of B. dorsaslis, and BdCht7 is dispensable for larval survival. It also provides comprehensive transcriptome information for finding more molecular targets to control this pest. © 2020 Society of Chemical Industry.

miR-34 regulates larval growth and wing morphogenesis by directly modulating ecdysone signalling and cuticle protein in Bombyx mori

microRNAs (miRNA) are small non-coding RNAs that modulate the myriad biological activities by targeting genes, and many studies showed that miRNAs played a pivotal role in insect development. Here, we find that Bm-miRNA (miR-34) controls larval growth and wing morphology by targeting BmE74 and BmCPG4. Overexpression of miR-34 in the whole body caused a smaller body size, partially displays deformed wings and venation defects in adults. Ablation of miR-34 by transgenic CRISPR/Cas9 technology resulted in a severe developmental delay during the larval stage. Moreover, we confirmed that miR-34 directly targeted BmE74 and BmCPG4 by using a dual luciferase reporter assay in HEK293T cells. Remarkably, loss-of-function of BmCPG4 caused wing defects, which was similar to the phenotype of miR-34 overexpression in animals. In addition, our analysis revealed that ecdysone strongly inhibited miR-34 expression in vivo. Taken together, our study identifies miR-34 as a modulator that regulates larval growth and wing morphogenesis by directly modulating ecdysone signalling and cuticle protein in Bombyx mori.

Effects of Osiris9a on Silk Properties in Bombyx mori Determined by Transgenic Overexpression

Osiris is an insect-specific gene family with multiple biological roles in development, phenotypic polymorphism, and protection. In the silkworm, we have previously identified twenty-five Osiris genes with high evolutionary conservation and remarkable synteny among several insects. Bombxy moriOsiris9a (BmOsi9a) is expressed only in the silk gland, particularly in the middle silk gland (MSG). However, the biological function of BmOsi9a is still unknown. In this study, we overexpressed BmOsi9a in the silk gland by germline transgene expression. BmOsi9a was overexpressed not only in the MSG but also in the posterior silk gland (PSG). Interestingly, BmOsi9a could be secreted into the lumen in the MSG but not in the PSG. In the silk fiber, overexpressed BmOsi9a interacted with Sericin1 in the MSG, as confirmed by a co-immunoprecipitation assay. The overexpression of BmOsi9a altered the secondary structure and crystallinity of the silk fiber, thereby changing the mechanical properties. These results provide insight into the mechanisms underlying silk proteins secretion and silk fiber formation.

Transcriptomic analysis of the testicular fusion in Spodoptera litura

Background

Lepidoptera is one group of the largest plant-feeding insects and Spodoptera litura (Lepidoptera: Noctuidae) is one of the most serious agricultural pests in Asia countries. An interesting and unique phenomenon for gonad development of Lepidoptera is the testicular fusion. Two separated testes fused into a single one during the larva-to-pupa metamorphosis, which is believed to contribute to sperm production and the prevalence in field. To study the molecular mechanism of the testicular fusion, RNA sequencing (RNA-seq) experiments of the testes from 4-day-old sixth instar larvae (L6D4) (before fusion), 6-day-old sixth instar larvae (L6D6, prepupae) (on fusing) and 4-day-old pupae (P4D) (after fusion) of S. litura were performed.

Results

RNA-seq data of the testes showed that totally 12,339 transcripts were expressed at L6D4, L6D6 and P4D stages. A large number of differentially expressed genes (DEGs) were up-regulated from L6D4 to L6D6, and then more genes were down-regulated from L6D6 to P4D. The DEGs mainly belongs to the genes related to the 20E signal transduction pathway, transcription factors, chitin metabolism related enzymes, the families of cytoskeleton proteins, extracellular matrix (ECM) components, ECM-related protein, its receptor integrins and ECM-remodeling enzymes. The expression levels of these genes that were up-regulated significantly during the testicular fusion were verified by qRT-PCR. The matrix metalloproteinases (MMPs) were found to be the main enzymes related to the ECM degradation and contribute to the testicular fusion. The testis was not able to fuse if MMPs inhibitor GM6001 was injected into the 5th abdomen region at L6D6 early stage.

Conclusions

The transcriptome and DEGs analysis of the testes at L6D4, L6D6 and P4D stages provided genes expression information related to the testicular fusion in S. litura. These results indicated that cytoskeleton proteins, ECM-integrin interaction genes and ECM-related proteins were involved in cell migration, adhesion and fusion during the testicular fusion. The ECM degradation enzymes MMPs probably play a critical role in the fusion of testis.

Induced Thermotolerance and Expression of Three Key Hsp Genes (Hsp70, Hsp21, and sHsp21) and Their Roles in the High Temperature Tolerance of Agasicles hygrophila

Thermal adaptation plays a fundamental role in the expansion and distribution of insects, and heat shock proteins (Hsps) play important roles in the temperature adaptation of various organisms. To determine the roles of Hsp genes (Hsp70, Hsp21, and sHsp21) on the high temperature tolerance of Agasicles hygrophila, we obtained complete cDNA (complementary DNA) sequences for Hsp70, Hsp21, and sHsp21 by rapid amplification of cDNA ends (RACE), analyzed their expression profiles under different high temperature treatments by quantitative reverse transcription polymerase chain reaction (RT-qPCR), and performed functional verification by RNA interference (RNAi). The open reading frames of Hsp70, Hsp21, and sHsp21 were 1940, 543, and 567 bp, encoding 650, 180, and 188 amino acids, respectively. Their molecular weights (MWs) were 71.757, 20.879, and 21.510 kDa, and the isoelectric points were 5.63, 6.45, and 6.24, respectively. Phylogenetic tree analysis showed that the Hsp70, Hsp21, and sHsp21 genes of A. hygrophila were relatively conserved in evolution. The Hsp70 and Hsp21 genes in A. hygrophila were homologous to those in Leptinotarsa decemlineata (87 and 79% similarity, respectively), and the sHsp21 gene in A. hygrophila was homologous to that in Lissorhoptrus oryzophilus (74% similarity). The amino acid polypeptide chain had highly conserved sequences of DLGGGTFD, VLVGGSTR, and GPTIEEVD. The sequence of EEVD was the characteristic motif of cytoplasmic Hsp70, and the highly conserved sequences of MALFR and MSLLP were characteristic sequences of Hsp2 and sHsp21, respectively. Relative quantitative real time PCR showed that the three Hsps could be induced by 4-h treatment at high temperatures. Significant upregulation of these Hsps was observed when the temperature was further increased. The RNAi results showed that the injection of the three Hsps' dsRNA could suppress the expression at the gene level significantly. Compared with the control group, high temperature heat shock reduced the fecundity of A. hygrophila significantly, and the fecundity decreased with the increase in temperature. Our results suggest that Hsp70, Hsp21, and sHsp21 might play key roles in high temperature adaptation of A. hygrophila and help improve our understanding of their mechanism of thermotolerance.

BR-C Z4 and FoxJ interact to regulate expression of a chitin synthase gene CHSA-2b in the pupal wing discs of the silkworm, Bombyx mori

Elaborate regulation of tissue- and stage-specific expression of genes is prerequisite for insect development. The hormone 20-hydroxyecdysone (20E) initiates metamorphosis by regulating the expression of a series of genes. However, how 20E orderly regulates the pupa-specific expression of genes remains unclear. In this study, we report a regulatory mechanism for the pupa-specific expression of chitin synthase A 2b (CHSA-2b) in Bombyx mori. We found that Broad-Complex Z4 (BR-C Z4) was up-regulated by 20E just before pupation, while transcription factor FoxJ and CHSA-2b were up-regulated during the pupal stage. There is a Fox cis-regulatory element in the CHSA-2b promoter region, and FoxJ protein bound to this element, enhancing the CHSA-2b transcription during the pupal stage. In addition to CHSA-2b, FoxJ also up-regulated the expression of 16 out of 19 pupa-specific genes tested. However, at the prepupal stage, 20E-induced BR-C Z4 inhibited the FoxJ transcription, indirectly inhibiting the CHSA-2b transcription. These data suggest that at the pre-pupation stage, 20E-induced BR-C Z4 inhibited the expression of pupa-stage genes like CHSA-2b by inhibiting the expression of FoxJ; by the pupal stage, the expression of BR-C Z4 decreased, releasing its inhibition on FoxJ, which then up-regulated the expression of the pupa-specific genes. This study explains the elaborate regulation of the pupa-specific gene expression during metamorphosis in B. mori.

BmBlimp-1 gene encoding a C2H2 zinc finger protein is required for wing development in the silkworm Bombyx mori

Cys2-His2 zinc finger (C2H2-ZF) proteins represent the most common class of transcription factors. These factors have great potential for the management of developmental progression by regulating the specific spatiotemporal expression of genes. In this study, we cloned one C2H2-ZF protein gene of Bombyx mori, BGIBMGA000319, that is orthologous to B-lymphocyte-induced maturation protein-1 (Blimp-1); we thus named it as Bombyx mori Blimp-1 (BmBlimp-1). In the silkworm, the BmBlimp-1 gene is specifically upregulated during day 2 of the pupal to adult stage and is highly expressed in wing discs on day 3 of the pupa. Knockdown of its expression level in the pupal stage results in a crumpled-winged silkworm moth. Using the predicted DNA-binding sequences of BmBlimp-1 to search the silkworm genome to screen target genes of BmBlimp-1, 7049 genes were identified to have at least one binding site of BmBlimp-1 on their 1 kb upstream and downstream genome regions. Comparisons of those genes with a reported pupal wing disc transcriptome data resulted in 4065 overlapping genes being retrieved. GO enrichment analysis of the overlapping genes showed that most of the genes were enriched in the binding term. Combining functional annotation and real-time quantitative PCR, 15 genes were identified as the candidate target genes of BmBlimp-1, including several wing cuticular protein genes, chitin synthase A, and wing disc development genes, such as Wnt1, cubitus interruptus (ci) and engrailed (en). Moreover, the amino acid sequence of the zinc finger motif of Blimp-1 gene was highly conserved among the 15 insect species. We propose that BmBlimp-1 is an important regulatory factor in silkworm wing development.

CRISPR Disruption of BmOvo Resulted in the Failure of Emergence and Affected the Wing and Gonad Development in the Silkworm Bombyx mori

The domesticated silkworm is an economically important insect that is widely used as a lepidopteran insect model. Although somatic sex determination in the silkworm is well characterized, germline sex determination is not. Here, we used the transgenic-based CRISPR/Cas9 genome editing system to study the function of the Ovo gene in Bombyx mori. BmOvo is the homolog of a factor important in germline sex determination in Drosophila melanogaster. BmOvo mutants had abnormally shaped eggs that were disordered in the ovarioles, and gonad development was abnormal. Interestingly, wing discs and wings did not develop properly, and most of the mutants failed to eclose. Gene expression analyses by qRT-PCR showed that BmOvo gene was highly expressed in the wing disc and epidermis. Genes involved in the WNT signaling pathway and wing development genes BmWCP10 and BmE74 were downregulated in the BmOvo mutants when compared with wild-type animals. These results demonstrate that the BmOvo gene product plays an important role in wing metamorphosis. Thus, this study provides new insights into the multiple functions of BmOvo beyond germline sex determination.

High-quality genome assembly of the silkworm, Bombyx mori

In 2008, the genome assembly and gene models for the domestic silkworm, Bombyx mori, were published by a Japanese and Chinese collaboration group. However, the genome assembly contains a non-negligible number of misassembled and gap regions due to the presence of many repetitive sequences within the silkworm genome. The erroneous genome assembly occasionally causes incorrect gene prediction. Here we performed hybrid assembly based on 140 × deep sequencing of long (PacBio) and short (Illumina) reads. The remaining gaps in the initial genome assembly were closed using BAC and Fosmid sequences, giving a new total length of 460.3 Mb, with 30 gap regions and an N50 comprising 16.8 Mb in scaffolds and 12.2 Mb in contigs. More RNA-seq and piRNA-seq reads were mapped on the new genome assembly compared with the previous version, indicating that the new genome assembly covers more transcribed regions, including repetitive elements. We performed gene prediction based on the new genome assembly using available mRNA and protein sequence data. The number of gene models was 16,880 with an N50 of 2154 bp. The new gene models reflected more accurate coding sequences and gene sets than old ones. The proportion of repetitive elements was also reestimated using the new genome assembly, and was calculated to be 46.8% in the silkworm genome. The new genome assembly and gene models are provided in SilkBase (http://silkbase.ab.a.u-tokyo.ac.jp).

Candidate olfactory genes identified in Heortia vitessoides (Lepidoptera: Crambidae) by antennal transcriptome analysis

Heortia vitessoides Moore is the most severe defoliating pest of Aquilaria sinensis (Lour.) Gilg (Thymelaeaceae) forests. Olfaction in insects is essential for host identification, mating, and oviposition, in which olfactory proteins, including odorant-binding proteins (OBPs), chemosensory proteins (CSPs), olfactory receptors (ORs), ionotropic receptors (IRs), and sensory neuron membrane proteins (SNMPs), are responsible for chemical signaling. Here, we determined the transcriptomes of male and female adult antennae of H. vitessoides. We assembled 52,383 unigenes and annotated their putative gene functions based on the gene ontology (GO), eukaryotic ortholog groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. Overall, 61 olfactory-related transcripts, including nine OBPs, 10 CSPs, 28 ORs, 12 IRs, and two SNMPs, were identified. Expression patterns of OBPs and CSPs in the female antennae, male antennae, and legs were performed using reverse transcription quantitative PCR (RT-qPCR). The results revealed that HvitOBP1, HvitOBP6, and HvitGOBP1 were enriched in the female antennae, while HvitOBP2, HvitOBP3, HvitOBP5, HvitGOBP2, and HvitPBP1 were enriched in the male antennae. HvitOBP4 was expressed at nearly the same level in the antennae of both males and females. Four CSPs (HvitCSP3, HvitCSP5, HvitCSP7, and HvitCSP10) and two CSPs (HvitCSP1 and HvitCSP4) were expressed at higher levels in the female and male antennae, respectively. HvitCSP6 was expressed at higher levels both in the female antennae and legs. Three CSP genes (HvitCSP2, HvitCSP8, and HvitCSP9) were expressed at higher levels in the legs. These results provide a basis for further studies on the molecular olfactory mechanisms of H. vitessoides.

RNA-seq analysis of gene expression changes during pupariation in Bactrocera dorsalis (Hendel) (Diptera: Tephritidae)

Background

The oriental fruit fly, Bactrocera dorsalis (Hendel) has been considered to be one of the most important agricultural pest around the world. As a holometabolous insect, larvae must go through a metamorphosis process with dramatic morphological and structural changes to complete their development. To better understand the molecular mechanisms of these changes, RNA-seq of B. dorsalis from wandering stage (WS), late wandering stage (LWS) and white puparium stage (WPS) were performed.

Results

In total, 11,721 transcripts were obtained, out of which 1914 genes (578 up-regulated and 1336 down-regulated) and 2047 genes (655 up-regulated and 1392 down-regulated) were found to be differentially expressed between WS and LWS, as well as between WS and WPS, respectively. Of these DEGs, 1862 and 1996 genes were successfully annotated in various databases. The analysis of RNA-seq data together with qRT-PCR validation indicated that during this transition, the genes in the oxidative phosphorylation pathway, and genes encoding P450s, serine protease inhibitor, and cuticular proteins were down-regulated, while the serine protease genes were up-regulated. Moreover, we found some 20-hydroxyecdysone (20E) biosynthesis and signaling pathway genes had a higher expression in the WS, while the genes responsible for juvenile hormone (JH) synthesis, degradation, signaling and transporter pathways were down-regulated, suggesting these genes might be involved in the process of larval pupariation in B. dorsalis. For the chitinolytic enzymes, the genes encoding chitinases (chitinase 2, chitinase 5, chitinase 8, and chitinase 10) and chitin deacetylase might play the crucial role in the degradation of insect chitin with their expressions significantly increased during the transition. Here, we also found that chitin synthase 1A might be involved in the chitin synthesis of cuticles during the metamorphosis in B. dorsalis.

Conclusions

Significant changes at transcriptional level were identified during the larval pupariation of B. dorsalis. Importantly, we also obtained a vast quantity of RNA-seq data and identified metamorphosis associated genes, which would all help us to better understand the molecular mechanism of metamorphosis process in B. dorsalis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5077-z) contains supplementary material, which is available to authorized users.

Dual role of the serine protease homolog BmSPH-1 in the development and immunity of the silkworm Bombyx mori

Serine proteases and serine protease homologs are involved in the prophenoloxidase (proPO)-activating system leading to melanization. The Bombyx mori serine protease homolog BmSPH-1 regulates nodule melanization. Here, we show the dual role of BmSPH-1 in the development and immunity of B. mori. BmSPH-1 was expressed in hemocytes after molting and during the larval-pupal transformation in normal development. In contrast, following infection, BmSPH-1 was expressed in hemocytes and cleaved in the hemolymph, which resulted in the induction of PO activity. Moreover, BmSPH-1 was cleaved in the cuticle during the larval-pupal transformation and early pupal stages. In BmSPH-1 RNAi-treated silkworms, the reduced BmSPH-1 mRNA levels during the spinning stage or the prepupal stage resulted in the arrest of pupation or pupal cuticular melanization, respectively. The binding assays revealed that BmSPH-1 interacts with B. mori immulectin, proPO, and proPO-activating enzyme. Our findings demonstrate that BmSPH-1 paticipates larval-pupal transformation, pupal cuticular melanization and innate immunity of silkworms, illustrating the dual role of BmSPH-1 in development and immunity.

Functional Identification and Characterization of the Diuretic Hormone 31 (DH31) Signaling System in the Green Shore Crab, Carcinus maenas

The functional characterization of crustacean neuropeptides and their cognate receptors has not kept pace with the recent advances in sequence determination and, therefore, our understanding of the physiological roles of neuropeptides in this important arthropod sub-phylum is rather limited. We identified a candidate receptor-ligand pairing for diuretic hormone 31 (DH31) in a neural transcriptome of the crab, Carcinus maenas. In insects, DH31 plays species -specific but central roles in many facets of physiology, including fluid secretion, myoactivity, and gut peristalsis but little is known concerning its functions in crustaceans. The C. maenas DH31 transcript codes for a 147 amino acid prepropeptide, and a single receptor transcript translates to a secretin-like (Class B1) G protein-coupled receptor (GPCR). We used an in vitro aequorin luminescence Ca2+ mobilization assay to demonstrate that this candidate DH31R is activated byCarcinus and insect DH31s in a dose-dependent manner (EC50 15-30 nM). Whole mount immunohistochemical and in situ hybridization localization revealed extensive DH31 expressing neurons throughout the central nervous system, most notably in the abdominal ganglion where large, unpaired cells give rise to medial nerves, which terminate in extensive DH31 immunopositive dendritic fields intimately associated with oesophageal musculature. This system constitutes a large and hitherto undescribed neurohemal area adjacent to key muscle groups associated with the gastric system. DH31 expressing neurons were also seen in the cardiac, commissural, oesophageal, and stomatogastric ganglia and intense labeling was seen in dendrites innervating fore- and hindgut musculature but not with limb muscles. These labeling patterns, together with measurement of DH31R mRNA in the heart and hindgut, prompted us test the effects of DH31 on semi-isolated heart preparations. Cardiac superfusion with peptide evoked increased heart rates (10-100 nM). The neuroanatomical distribution of DH31 and its receptor transcripts, particularly that associated with gastric and cardiac musculature, coupled with the cardio- acceleratory effects of the peptide implicate this peptide in key myoactive roles, likely related to rhythmic coordination.

DNA methylation mediates BmDeaf1-regulated tissue- and stage-specific expression of BmCHSA-2b in the silkworm, Bombyx mori

BACKGROUND

Accurate regulation of tissue- and stage-specific expression of genes is prerequisite for normal development in organisms. DNA methylation plays an important role in modulating gene expression in mammals and plants. However, there is no direct evidence showing how DNA methylation regulates gene transcription in insects.

RESULTS

During the development of Bombyx mori wing, the expression level of DNA methyltransferase 1 (BmDnmt1) gradually declined and became stationary at pupal stage, resulting in a lower methylation rate of the intragenic promoter of the mid-pupal wing-specific gene BmCHSA-2b, an epidermal chitin synthase controlling mid-pupal wing development in B. mori. The higher methylation rate of the promoter in the pupal epidermis was decreased and BmCHSA-2b transcription was significantly increased by the treatment with the DNA methylation inhibitor, 5-azacytidine-2'-deoxycytidine, suggesting that DNA methylation regulates the tissue-specific expression of BmCHSA-2b. Pupa-specific transcription factor BmDEAF1 bound to the unmethylated intragenic promoter and activated the BmCHSA-2b transcription in the mid-pupal wing. BmDnmt1 and BmDeaf1 influenced the BmCHSA-2b transcription by binding competitively to the CpG island in the promoter.

CONCLUSIONS

All the data together demonstrate that the cooperation between the down-regulation of BmDnmt1 and increased stage-specific expression of BmDeaf1 enhances BmCHSA-2b tissue- and stage-specific transcription to ensure mid-wing development in B. mori. This study highlights an elaborate regulation mechanism how tissue- and stage-specific gene expression is regulated through promoter methylation in insect development.

Innate immune responses in the Chinese oak silkworm, Antheraea pernyi

Innate immunity, the evolutionarily conserved defense system, has been extensively analyzed in insect models over recent decades. The significant progress in this area has formed our dominant conceptual framework of the innate immune system, but critical advances in other insects have had a profound impact on our insights into the mystery of innate immunity. In recent years, we focused on the immune responses in Antheraea pernyi, an important commercial silkworm species reared in China. Here, we review the immune responses of A. pernyi based on immune-related gene-encoded proteins that are divided into five categories, namely pattern recognition receptors, hemolymph proteinases and their inhibitors, prophenoloxidase, Toll pathway factors and antimicrobial peptides, and others. Although the summarized information is limited since the research on A. pernyi immunity is in its infancy, we hope to provide evidence for further exploration of innate immune mechanisms.

Dynamic regulation of mRNA and miRNA associated with the developmental stages of skin pigmentation in Japanese ornamental carp

The Japanese ornamental carp (Cyprinus carpio var. Koi) is famous for multifarious colors and patterns, making it commonly culture and trade across the world. Although functional genes and inheritance of color traits have been commonly studied, seldom attentions were focused on the genetic regulation during the developmental process of pigmentation. To better understand the mechanism of skin color development, we observed the morphogenesis of pigment cells during the post-embryonic stages and analysed the temporal expression pattern of mRNAs/miRNAs profiles in four distinct developmental stages. 59 and 103 differentially expressed genes/miRNAs (DEGs/DEMs) associated with pigmentation and skin were identified, including pax7, mitf, tyr, tyrp1, etc., and the highest DEGs were detected at 11 days post hatching (dph). In addition, the functional characteristics of mRNAs/miRNAs associated with pteridine and carotenoid pathway were also examined. Furthermore, 65 miRNA-mRNA interaction pairs related to pigmentation, pteridines and carotenoids metabolism were detected between different stages. Interestingly, the largest pairs appeared in the transition from 11 dph to 48 dph, which had the similar trend with DEGs further manifesting the importance of 11 dph. This study produced a comprehensive programme of DEGs/DEMs during color development, which will provide resources to understand the regulation mechanism in color formation. The understanding of genetic basis in color formation might promote the production and breeding of the Koi carp.

Transcriptome profiling with focus on potential key genes for wing development and evolution in Megaloprepus caerulatus, the damselfly species with the world's largest wings

The evolution, development and coloration of insect wings remains a puzzling subject in evolutionary research. In basal flying insects such as Odonata, genomic research regarding bauplan evolution is still rare. Here we focus on the world's largest odonate species-the "forest giant" Megaloprepus caerulatus, to explore its potential for looking deeper into the development and evolution of wings. A recently discovered cryptic species complex in this genus previously considered monotypic is characterized by morphological differences in wing shape and color patterns. As a first step toward understanding wing pattern divergence and pathways involved in adaptation and speciation at the genomic level, we present a transcriptome profiling of M. caerulatus using RNA-Seq and compare these data with two other odonate species. The de novo transcriptome assembly consists of 61,560 high quality transcripts and is approximately 93% complete. For almost 75% of the identified transcripts a possible function could be assigned: 48,104 transcripts had a hit to an InterPro protein family or domain, and 28,653 were mapped to a Gene Ontology term. In particular, we focused on genes related to wing development and coloration. The comparison with two other species revealed larva-specific genes and a conserved 'core' set of over 8,000 genes forming orthologous clusters with Ischnura elegans and Ladona fulva. This transcriptome may provide a first point of reference for future research in odonates addressing questions surrounding the evolution of wing development, wing coloration and their role in speciation.

Osiris9a is a major component of silk fiber in lepidopteran insects

In a previous high-throughput proteomics study, it was found that the silkworm cocoon contains hundreds of complex proteins, many of which have unknown functions, in addition to fibroins, sericins, and some protease inhibitors. Osiris was one of the proteins with no known function. In this study, we identified the Osiris gene family members and constructed a phylogenetic tree based on the sequences from different species. Our results indicate that the Osiris9 gene subfamily contains six members; it is specifically expressed in lepidopteran insects and has evolved by gene duplication. An Osiris gene family member from Bombyx mori was designated as BmOsiris9a (BmOsi9a) on the basis of its homology to Drosophila melanogaster Osiris9. The expression pattern of BmOsi9a showed that it was highly expressed only in the middle silk gland of silkworm larvae, similar to Sericin1 (Ser1). BmOsi9a was visualized as two bands in western blot analysis, implying that it probably undergoes post-translational modifications. Immunohistochemistry analysis revealed that BmOsi9a was synthesized and secreted into the lumen of the middle silk gland, and was localized in the sericin layer in the silk fiber. BmOsi9a was found in the silk fibers of not only three Bombycidae species, viz. B. mori, B. mandarina, and B. huttoni, but also in the fibers collected from Saturniidae species, including Antheraea assama, Antheraea mylitta, and Samia cynthia. Although the exact biological function of Osi9a in the silk fibers is unknown, our results are important because they demonstrate that Osi9a is a common structural component of silk fiber and is expressed widely among the silk-producing Bombycidae and Saturniidae insects. Our results should help in understanding the role of Osi9a in silk fibers.

RNA-seq of Rice Yellow Stem Borer Scirpophaga incertulas Reveals Molecular Insights During Four Larval Developmental Stages

The yellow stem borer (YSB), Scirpophaga incertulas, is a prominent pest in rice cultivation causing serious yield losses. The larval stage is an important stage in YSB, responsible for maximum infestation. However, limited knowledge exists on the biology and mechanisms underlying the growth and differentiation of YSB. To understand and identify the genes involved in YSB development and infestation, so as to design pest control strategies, we performed de novo transcriptome analysis at the first, third, fifth, and seventh larval developmental stages employing Illumina Hi-seq. High-quality reads (HQR) of ∼229 Mb were assembled into 24,775 transcripts with an average size of 1485 bp. Genes associated with various metabolic processes, i.e., detoxification mechanism [CYP450, GSTs, and carboxylesterases (CarEs)], RNA interference (RNAi) machinery (Dcr-1, Dcr-2, Ago-1, Ago-2, Sid-1, Sid-2, Sid-3, and Sid-1-related gene), chemoreception (CSPs, GRs, OBPs, and ORs), and regulators [transcription factors (TFs) and hormones] were differentially regulated during the developmental stages. Identification of stage-specific transcripts made it possible to determine the essential processes of larval development. Comparative transcriptome analysis revealed that YSB has not evolved much with respect to the detoxification mechanism, but showed the presence of distinct RNAi machinery. The presence of strong specific visual recognition coupled with chemosensory mechanisms supports the monophagous nature of YSB. Designed expressed sequenced tags-simple-sequence repeats (EST-SSRs) will facilitate accurate estimation of the genetic diversity of YSB. This is the first report on characterization of the YSB transcriptome and the identification of genes involved in key processes, which will help researchers and industry to devise novel pest control strategies. This study also opens up a new avenue to develop next-generation resistant rice using RNAi or genome editing approaches.

BmCHSA-2b, a Lepidoptera specific alternative splicing variant of epidermal chitin synthase, is required for pupal wing development in Bombyx mori

Insect chitin synthase A (CHSA) is an epidermis-specific enzyme that plays an essential role in insect development. In this study, the function and regulation of CHSA-2b, an alternative splicing variant of Bombxy mori CHSA that is discovered only in Lepidopteran insects, were investigated. Analysis of mRNA level showed that BmCHSA-2b was responsive to 20-hydroxyecdysone (20E) in pupal wing unlike BmCHSA-2a, which shares almost the identical sequence as BmCHSA-2b except the first 31 amino acids, suggesting that the expression of these two alternative splicing variants is driven by different promoters of CHSA gene. Reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR) analysis showed that BmCHSA-2b was up-regulated in the wing of mid-pupa unlike BmCHSA-2a, which was up-regulated in epidermis and wing disc at the beginning and end of pupal stage. Further analysis reveals that the up-regulations of BmCHSA-2a and BmCHSA-2b in pupal wing were consistent with the increase of chitin content and wing area at the same stages, respectively. Furthermore, the higher transcription level of BmCHSA-2b in the mid-pupal wing of male than that in female was consistent with the chitin content of pupal wing between genders. Injection of double-stranded RNAs of BmCHSA-2b resulted in the decrease in the area and chitin content of the wing, and irregular and crimpled vein. All these results together suggest that B. mori evolves an extra promoter in CHSA gene to activate BmCHSA-2b expression in the wing of mid-pupal stage in response to 20E, and BmCHSA-2b is required for the wing development in the mid-pupa of B. mori.

De novo transcriptome of the mayfly Cloeon viridulum and transcriptional signatures of Prometabola

Mayflies (Ephemeroptera) display many primitive characters and a unique type of metamorphosis (Prometabola). However, information on the genomes and transcriptomes of this insect group is limited. The RNA sequencing study presented here generated the first de novo transcriptome assembly of Cloeon viridulum (Ephemeroptera: Baetidae), and compared gene expression signatures among the young larva (YL), mature larva (ML), subimago (SI), and imago (IM) stages of this mayfly. The transcriptome, based on 88 Gb of sequence data, comprised a set of 81,185 high quality transcripts. The number of differentially expressed genes (DEGs) in YL vs. ML, ML vs. SI, and SI vs. IM, was 4,825, 1,584, and 1,278, respectively, according to the reads per kilobase of transcript per million mapped reads analysis, assuming a false discovery rate <0.05 and a fold change >2. Gene enrichment analysis revealed that these DEGs were enriched in the "chitin metabolic process", "germ cell development", "steroid hormone biosynthesis", and "cutin, suberine, and wax biosynthesis" pathways. Finally, the expression pattern of a selected group of candidate signature genes for Prometabola, including vestigial, methoprene-tolerant, wingless, and broad-complex were confirmed by quantitative real time-PCR analysis. The Q-PCR analysis of larval, subimaginal, and imaginal stages of C. viridulum suggests that the development of mayflies more closely resembles hemimetamorphosis than holometamorphosis.

Comparative transcriptomes analysis of the wing disc between two silkworm strains with different size of wings

Wings of Bombyx mori (B. mori) develop from the primordium, and different B. mori strains have different wing types. In order to identify the key factors influencing B. mori wing development, we chose strains P50 and U11, which are typical for normal wing and minute wing phenotypes, respectively. We dissected the wing disc on the 1^st-day of wandering stage (P50D1 and U11D1), 2^nd-day of wandering stage (P50D2 and U11D2), and 3^rd-day of wandering stage (P50D3 and U11D3). Subsequently, RNA-sequencing (RNA-Seq) was performed on both strains in order to construct their gene expression profiles. P50 exhibited 628 genes differentially expressed to U11, 324 up-regulated genes, and 304 down-regulated genes. Five enriched gene ontology (GO) terms were identified by GO enrichment analysis based on these differentially expressed genes (DEGs). KEGG enrichment analysis results showed that the DEGs were enriched in five pathways; of these, we identified three pathways related to the development of wings. The three pathways include amino sugar and nucleotide sugar metabolism pathway, proteasome signaling pathway, and the Hippo signaling pathway. The representative genes in the enrichment pathways were further verified by quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). The RNA-Seq and qRT-PCR results were largely consistent with each other. Our results also revealed that the significantly different genes obtained in our study might be involved in the development of the size of B. mori wings. In addition, several KEGG enriched pathways might be involved in the regulation of the pathways of wing formation. These results provide a basis for further research of wing development in B. mori.

Large-scale identification of differentially expressed genes during pupa development reveals solute carrier gene is essential for pupal pigmentation in Chilo suppressalis

Insects undergo metamorphosis, involving an abrupt change in body structure through cell growth and differentiation. Rice stem stripped borer (SSB), Chilo suppressalis, is one of the most destructive rice pests. However, little is known about the regulation mechanism of metamorphosis development in this notorious insect pest. Here, we studied the expression of 22,197 SSB genes at seven time points during pupa development with a customized microarray, identifying 622 differentially expressed genes (DEG) during pupa development. Gene ontology (GO) analysis of these DEGs indicated that the genes related to substance metabolism were highly expressed in the early pupa, which participate in the physiological processes of larval tissue disintegration at these stages. In comparison, highly expressed genes in the late pupal stages were mainly associated with substance biosynthesis, consistent with adult organ formation at these stages. There were 27 solute carrier (SLC) genes that were highly expressed during pupa development. We knocked down SLC22A3 at the prepupal stage, demonstrating that silencing SLC22A3 induced a deficiency in pupa stiffness and pigmentation. The RNAi-treated individuals had white and soft pupa, suggesting that this gene has an essential role in pupal development.

Knockdown of two trehalose-6-phosphate synthases severely affects chitin metabolism gene expression in the brown planthopper Nilaparvata lugens

BACKGROUND

RNA interference combined with digital gene expression (DGE) analysis can be used to study gene function. Trehalose-6-phosphate synthase (TPS) plays a key role in the synthesis of trehalose and insect development.

RESULTS

DGE analysis revealed that the expression of nine or four chitinase genes was reduced significantly 48 h after NlTPS1 and NlTPS2 knockdown by RNAi, respectively. Additionally, abnormal phenotypes were noted, and approximately 30% of insects died. HK and G6PI2 expression decreased significantly whereas GFAT, GNPNA and UAP expression increased significantly 72 h after NlTPS1 and NlTPS2 knockdown. PGM1 expression decreased significantly after TPS2 knockdown, whereas PGM2 expression increased significantly and the expression of three CHS genes decreased 48 h after TPS1 knockdown. The mRNA expression of all 12 chitin degradation genes decreased 48 h after NlTPS1 and NlTPS2 treatment, and Cht2, Cht3, Cht6, Cht7, Cht10 and ENGase levels remained significantly decreased up to 72 h after NlTPS1 and NlTPS2 knockdown.

CONCLUSIONS

These results demonstrate that silencing of TPS genes can lead to increased moulting deformities and mortality rates owing to the misregulation of genes involved in chitin metabolism, and TPS genes are potential pest control targets in the future. © 2016 Society of Chemical Industry.

Analysis of expression and chitin-binding activity of the wing disc cuticle protein BmWCP4 in the silkworm, Bombyx mori

The insect exoskeleton is mainly composed of chitin filaments linked by cuticle proteins. When insects molt, the cuticle of the exoskeleton is renewed by degrading the old chitin and cuticle proteins and synthesizing new ones. In this study, chitin-binding activity of the wing disc cuticle protein BmWCP4 in Bombyx mori was studied. Sequence analysis showed that the protein had a conservative hydrophilic "R&R" chitin-binding domain (CBD). Western blotting showed that BmWCP4 was predominately expressed in the wing disc-containing epidermis during the late wandering and early pupal stages. The immunohistochemistry result showed that the BmWCP4 was mainly present in the wing disc tissues containing wing bud and trachea blast during day 2 of wandering stage. Recombinant full-length BmWCP4 protein, "R&R" CBD peptide (CBD), non-CBD peptide (BmWCP4-CBD^- ), four single site-directed mutated peptides (M₁ , M₂ , M₃ and M₄ ) and four-sites-mutated peptide (M_F ) were generated and purified, respectively, for in vitro chitin-binding assay. The results indicated that both the full-length protein and the "R&R" CBD peptide could bind with chitin, whereas the BmWCP4-CBD^- could not bind with chitin. The single residue mutants M₁ , M₂ , M₃ and M₄ reduced but did not completely abolish the chitin-binding activity, while four-sites-mutated protein M_F completely lost the chitin-binding activity. These data indicate that BmWCP4 protein plays a critical role by binding to the chitin filaments in the wing during larva-to-pupa transformation. The conserved aromatic amino acids are critical in the interaction between chitin and the cuticle protein.