Cloning of cDNAs encoding sorbitol dehydrogenase-2a and b, enzymatic characterization, and up-regulated expression of the genes in Bombyx mori diapause eggs exposed to 5 °C

Identification of genes associated with the silk gland size using multi-omics in silkworm (Bombyx mori)

Silk gland size in silkworms (Bombyx mori) affects silk output. However, the molecular mechanisms by which genes regulate silk gland size remain unclear. In this study, silk glands from three pure silkworm strains (A798, A306 and XH) with different silk gland weight phenotypes were compared using transcriptomics and proteomics to identify differentially expressed genes (DEGs) and proteins (DEPs). When comparing A798 to A306 and A798 to XH, 830 and 469 DEGs were up-regulated, respectively. These genes were related to the gene ontology terms, metabolic process, transport activity and biosynthesis process. In addition, 372 and 302 up-regulated differentially expressed proteins were detected in A798 to A306 and A798 to XH, respectively, related to the gene ontology terms, ribosome and protein export, ribosome and polypeptide biosynthesis processes. Moreover, combined transcriptomics, proteomics and weighted correlation network analyses showed that five genes (BGIBMGA002524, BGIBMGA002629, BGIBMGA005659, BGIBMGA005711 and BGIBMGA010889) were significantly associated with the silk gland weight. Reverse Transcription-quantitative real-time Polymerase Chain Reaction (RT-qPCR) and Enzyme linked immunosorbent assay (ELISA) were used to verify the mRNA and protein expression of five genes in the silk glands and tissues of 18 silkworm strains. The results showed that four genes have higher expression levels in heavier silk glands. These genes are associated with glycogen metabolism, fatty acid synthesis and branched chain amino acid metabolism, thus potentially promoting growth and silk protein synthesis. These findings provide valuable insights into the molecular mechanisms underlying the relationship between silk gland weight and silk yield in silkworms.

Noncoding RNA Regulation of Dormant States in Evolutionarily Diverse Animals

Dormancy is evolutionarily widespread and can take many forms, including diapause, dauer formation, estivation, and hibernation. Each type of dormancy is characterized by distinct features; but accumulating evidence suggests that each is regulated by some common processes, often referred to as a common "toolkit" of regulatory mechanisms, that likely include noncoding RNAs that regulate gene expression. Noncoding RNAs, especially microRNAs, are well-known regulators of biological processes associated with numerous dormancy-related processes, including cell cycle progression, cell growth and proliferation, developmental timing, metabolism, and environmental stress tolerance. This review provides a summary of our current understanding of noncoding RNAs and their involvement in regulating dormancy.

Comparative transcriptomic analysis of Tibetan Gynaephora to explore the genetic basis of insect adaptation to divergent altitude environments

Adaptation of insects to different altitudes remain largely unknown, especially those endemic to the Tibetan Plateau (TP). Here, we generated the transcriptomes of Gynaephora menyuanensis and G. alpherakii, inhabiting different high altitudes on the TP, and used these and the previously available transcriptomic and genomic sequences from low-altitude insects to explore potential genetic basis for divergent high-altitude adaptation in Gynaephora. An analysis of 5,869 orthologous genes among Gynaephora and other three low-altitude insects uncovered that fast-evolving genes and positively selected genes (PSGs) in the two Gynaephora species were enriched in energy metabolism and hypoxia response categories (e.g. mitochondrion, oxidation-reduction process, and response to oxidative stress). Particularly, mTOR signaling pathway involving hypoxia was enriched by PSGs, indicating this well-known pathway in mammal hypoxia adaptation may be an important signaling system in Gynaephora. Furthermore, some PSGs were associated with response to hypoxia (e.g. cytochrome proteins), cold (e.g. dehydrogenase) and DNA repair (e.g. DNA repair proteins). Interestingly, several insect-specific genes that were associated with exoskeleton and cuticle development (e.g. chitinase and ecdysteroids) had experienced positive selection, suggesting the specific adaptive mechanisms in insects. This study is favourable for understanding the adaptive evolution of Gynaephora and even TP insects to divergent altitudes.

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.

Differentially expressed microRNAs in diapausing versus HCl-treated Bombyx embryos

Differentially expressed microRNAs were detected to explore the molecular mechanisms of diapause termination. The total small RNA of diapause-destined silkworm eggs and HCl-treated eggs was extracted and then sequenced using HiSeq high-throughput method. 44 novel miRNAs were discovered. Compared to those in the diapause-destined eggs, 61 miRNAs showed significant changes in the acid-treated eggs, with 23 being up-regulated and 38 being down-regulated. The potential target genes of differentially expressed miRNAs were predicted by miRanda. Gene Ontology and KEGG pathway enrichment analysis of these potential target genes revealed that they were mainly located within cells and organelles, involved in cellular and metabolic processes, and participated in protein production, processing and transportation. Two differentially expressed genes, Bombyx mori SDH and Bmo-miR-2761-3p, were further analyzed with qRT-PCR. BmSDH was significantly up-regulated in the HCl-treated eggs, while Bmo-miR-2761-3p was down-regulated. These results suggested that these two genes were well coordinated in silkworm eggs. Dual luciferase reporter assay demonstrated that Bmo-miR-2761-3p inhibited the expression of BmSDH.

Energy and lipid metabolism during direct and diapause development in a pierid butterfly

Diapause is a fundamental component of the life-cycle in the majority of insects living in environments characterized by strong seasonality. The present study addresses poorly understood associations and trade-offs between endogenous diapause duration, thermal sensitivity of development, energetic cost of development and cold tolerance. Diapause intensity, metabolic rate trajectories and lipid profiles of directly developing and diapausing animals were studied using pupae and adults of Pieris napi butterflies from a population for which endogenous diapause is well studied. Endogenous diapause was terminated after 3 months and termination required chilling. Metabolic and postdiapause development rates increased with diapause duration, while the metabolic cost of postdiapause development decreased, indicating that once diapause is terminated development proceeds at a low rate even at low temperature. Diapausing pupae had larger lipid stores than the directly developing pupae and lipids constituted the primary energy source during diapause. However, during diapause lipid stores did not decrease. Thus, despite lipid catabolism meeting the low energy costs of the diapausing pupae, primary lipid store utilization did not occur until the onset of growth and metamorphosis in spring. In line with this finding, diapausing pupae contained low amounts of mitochondria-derived cardiolipins, which suggests a low capacity for fatty acid β-oxidation. While ontogenic development had a large effect on lipid and fatty acid profiles, only small changes in these were seen during diapause. The data therefore indicate that the diapause lipidomic phenotype is built early, when pupae are still at high temperature, and retained until diapause post-diapause development.

A de novo transcriptome and valid reference genes for quantitative real-time PCR in Colaphellus bowringi

Background

The cabbage beetle Colaphellus bowringi Baly is a serious insect pest of crucifers and undergoes reproductive diapause in soil. An understanding of the molecular mechanisms of diapause regulation, insecticide resistance, and other physiological processes is helpful for developing new management strategies for this beetle. However, the lack of genomic information and valid reference genes limits knowledge on the molecular bases of these physiological processes in this species.

Results

Using Illumina sequencing, we obtained more than 57 million sequence reads derived from C. bowringi, which were assembled into 39,390 unique sequences. A Clusters of Orthologous Groups classification was obtained for 9,048 of these sequences, covering 25 categories, and 16,951 were assigned to 255 Kyoto Encyclopedia of Genes and Genomes pathways. Eleven candidate reference gene sequences from the transcriptome were then identified through reverse transcriptase polymerase chain reaction. Among these candidate genes, EF1α, ACT1, and RPL19 proved to be the most stable reference genes for different reverse transcriptase quantitative polymerase chain reaction experiments in C. bowringi. Conversely, aTUB and GAPDH were the least stable reference genes.

Conclusion

The abundant putative C. bowringi transcript sequences reported enrich the genomic resources of this beetle. Importantly, the larger number of gene sequences and valid reference genes provide a valuable platform for future gene expression studies, especially with regard to exploring the molecular mechanisms of different physiological processes in this species.

Genome wide microarray based expression profiles during early embryogenesis in diapause induced and non-diapause eggs of polyvoltine silkworm Bombyx mori

Diapause was induced in polyvoltine silkworm B. mori eggs and the molecular mechanism involved in diapause was investigated using a genome wide microarray. In diapause eggs, 638 and 675 genes were upregulated, while, in non-diapause eggs 1136 and 595 genes were upregulated at 18 h and 30 h, respectively after oviposition. Real-time qPCR analysis confirmed the expression of 20 genes, and the relative expression levels of the Aquaporin gene was highest among the 20 genes, followed by Sorbitol dehydrogenase-2 and Cytochrome b5 in diapause eggs, while, Kruppel homolog, Period and Relish were higher in non-diapause eggs. The upregulation of SDH-2 and cytochrome b5 indicates increased metabolic rate in diapause-destined embryos prior to the onset of diapause within 36 h as a preparatory phase. This study provides an insight into the early molecular events for the induction and maintenance of diapause in B. mori.

Diapause specific gene expression in the eggs of multivoltine silkworm Bombyx mori, identified by suppressive subtractive hybridization

Molecular mechanism controlling egg diapause remains obscure in silkworm, Bombyx mori. An attempt is made to decipher various molecular events occurring during embryonic diapause in multivoltine silkworm, B. mori. Using suppressive subtractive hybridization (SSH), 186 cDNA clones isolated from both diapause and nondiapause eggs were sequenced. Of the sequenced clones, 29 matched with silkbase entries and these identified putative genes were classified into six functional groups such as regulatory, food utilization, stress response, metabolic, ribosomal and transposable elements. Among these genes, twelve belonged to regulatory group while, one taste receptor type 2 member 117 gene was related to food utilization. One heat shock cognate 70 kDa protein and 3 of the ubiquitin family were identified under stress response category. Similarly, four genes were identified as metabolic genes, 3 belonging to chitin family and one propanediol utilization protein. Of the seven genes identified in ribosomal groups, most of them were 60s ribosomal protein subunits. However, one negative regulation of transcription gene identified was a transposable element. The qPCR analysis confirmed the expression of 21 of the above genes, wherein, 6 were upregulated during diapause, 12 during non-diapause, while, 3 remained unchanged.

Heat shock factor binds to heat shock elements upstream of heat shock protein 70a and Samui genes to confer transcriptional activity in Bombyx mori diapause eggs exposed to 5°C

To understand the molecular mechanisms of how 5°C-incubation activates mRNA expression of Hsp70a and Samui genes in Bombyx mori diapause eggs, we first searched the 5'-upstream regions of the Hsp70a and Samui genes for heat shock elements (HSEs) and found two regions [Hsp70aHSE-1 (-95 to -58) and -2 (-145 to -121), and SamuiHSE-1 (-84 to -55) and -2 (-304 to -290)] corresponding to HSEs (repeats of nGAAn and/or nTTCn). We cloned four cDNAs encoding heat shock factor (HSF)-a2 (627 amino acids), -b (685 aa), -c (682 aa) and -d (705 aa), which were produced by alternative splicing. When we exposed diapause eggs to 5°C beginning at 2 day post-oviposition to break diapause, HSFd mRNA only increased after chilling for 6-8 days, a pattern very similar to those of Hsp70a and Samui mRNAs. To examine further whether HSFd binds to the respective HSEs, we carried out a gel shift assay using HSFd protein expressed in a cell-free system and the isolated HSEs; migration of the respective digoxigenin(DIG)-labeled HSE-1 and -2 of Hsp70a and Samui was retarded by addition of HSFd; the retarded bands disappeared after addition of the corresponding unlabeled HSE-1 and -2 as competitors, but were not affected by addition of the respective mutated unlabeled HSE-1 and -2. These results indicated that HSFd protein binds to the respective HSEs and may activate mRNA expression of Hsp70a and Samui genes upon exposure of diapause eggs to 5°C.