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

Differential transcriptomic responses to heat stress in surface and subterranean diving beetles

Subterranean habitats are generally very stable environments, and as such evolutionary transitions of organisms from surface to subterranean lifestyles may cause considerable shifts in physiology, particularly with respect to thermal tolerance. In this study we compared responses to heat shock at the molecular level in a geographically widespread, surface-dwelling water beetle to a congeneric subterranean species restricted to a single aquifer (Dytiscidae: Hydroporinae). The obligate subterranean beetle Paroster macrosturtensis is known to have a lower thermal tolerance compared to surface lineages (CT_max 38 °C cf. 42–46 °C), but the genetic basis of this physiological difference has not been characterized. We experimentally manipulated the thermal environment of 24 individuals to demonstrate that both species can mount a heat shock response at high temperatures (35 °C), as determined by comparative transcriptomics. However, genes involved in these responses differ between species and a far greater number were differentially expressed in the surface taxon, suggesting it can mount a more robust heat shock response; these data may underpin its higher thermal tolerance compared to subterranean relatives. In contrast, the subterranean species examined not only differentially expressed fewer genes in response to increasing temperatures, but also in the presence of the experimental setup employed here alone. Our results suggest P. macrosturtensis may be comparatively poorly equipped to respond to both thermally induced stress and environmental disturbances more broadly. The molecular findings presented here have conservation implications for P. macrosturtensis and contribute to a growing narrative concerning weakened thermal tolerances in obligate subterranean organisms at the molecular level.

Transcriptomic analysis of pre-diapause larvae of Chilo suppressalis (Walker) (Lepidoptera: Pyralidae) in natural populations

Chilo suppressalis Walker is a devastating pest of rice in Asia and exhibits facultative diapause in the larval stage. Most prior experiments on diapausing and non-diapausing C. suppressalis were conducted in the laboratory. In this study, transcriptome analyses were performed on pre-diapausing larvae collected from field populations of C. suppressalis and compared to laboratory populations. Among 2674 differentially expressed genes (DEGs), 32 DEGs related to pre-diapause and 239 universally expressed genes were screened; these were primarily enriched in "neuroactive ligand-receptor interaction", "lysosome" and "glycerolipid metabolism" in KEGG pathway analysis. With respect to clusters of orthologous genes (COG), DEGs were assigned to "posttranslational modification, protein turnover, chaperones", "carbohydrate transport and metabolism", and "secondary metabolite biosynthesis, transport and catabolism" categories. Further analysis also revealed that a key "circadian clock-controlled protein" gene is sensitive to photoperiod and significantly decreased during the pre-diapause phase. Genes encoding two small heat shock proteins, hsp21.4 and hsp27.2, were significantly expressed on August 15 as compared to three other sampling times in August 2018. Eight DEGs were randomly chosen and evaluated by real-time quantitative PCR (RT-qPCR) to validate the accuracy of the transcriptome data. The expression of six DEGs (gene-evm_000752, gene-evm_006486, gene-evm_008626, gene-evm_002485, gene-evm_011981 and Chilo_suppressalis_newGene_18103) showed significant same patterns of differential expression in both the RT-qPCR and RNA-Seq analyses. This study increases our understanding of the complex physiological and molecular mechanisms involved in C. suppressalis at the pre-diapause phase.

The role of N6-methyladenosine modification on diapause in silkworm (Bombyx mori) strains that exhibit different voltinism

N6-methyladenosine (m⁶ A) plays a key role in regulating gene expression in myriad organisms. Diapause is an important plastic phenotype that allows insects to survive under specific environmental conditions. However, the diapause molecular mechanism remains unknown. In this study, we analyzed the phylogenetics of genes related to the m⁶ A modification complex in the silkworm (Bombyx mori) based on identified sequences from other organisms. We detected the expression of these genes during different developmental phases from four strains with different voltinism. We also determined total m⁶ A content in cells treated with different diapause hormone concentrations or eggs exposed to hydrochloric acid. Our data revealed that m⁶ A-modification-related gene expression and m⁶ A content were greater in diapause-destinated compared to nondiapause-destined strains. Our findings suggest that m⁶ A modification may provide significant epigenetic regulation of diapause-related genes in the silkworm.

Expression Analysis Reveals the Association of Several Genes with Pupal Diapause in Bactrocera minax (Diptera: Tephritidae)

The Chinese citrus fly, Bactrocera minax, is a devastating pest of citrus, which enters the obligatory diapause in overwintering pupae to resist harsh environmental conditions. However, little is known about the molecular mechanisms underlying pupal diapause. The previous transcriptomic analysis revealed that a large number of genes were regulated throughout the pupal stage. Of these genes, 12 and six ones that are remarkably up- and downregulated, respectively, specifically in intense diapause were manually screened out in present study. To validate the expression of these genes throughout the pupal stage, the quantitative real-time PCR (qRT-PCR) was conducted, and the genes displaying different expression patterns with those of previous study were excluded. Then, the expressions of remaining genes were compared between diapause-destined and non-diapause-destined pupae to reveal their association with diapause using qRT-PCR and semiquantitative PCR. Finally, five genes, TTLL3B, Cyp6a9, MSTA, Fru, and UC2, and two genes, KSPI and LYZ1, were demonstrated to be positively and negatively associated with diapause, respectively. These findings provide a solid foundation for the further investigation of molecular mechanisms underlying B. minax pupal diapause.

Expression analysis and functional identification of several genes related to diapause in Bombyx mori

Diapause is an important characteristic of insects used to adapt to extreme changes in environmental conditions. Embryonic diapause of the bivoltine silkworm (Bombyx mori) is determined by environmental conditions experienced by the mother while in the embryo stage. If they are incubated at 25°C with natural light, their progenies will be diapause-destined. If they are incubated at 17°C in darkness, their progenies will be non-diapause-destined. The molecular mechanism of diapause remains unknown. In the present study, we analyzed two downregulated genes (BGIBMGA003835, BGIBMGA012335) and two upregulated genes (BGIBMG012996, BGIBMG002426) related to carbohydrate metabolism, verified differentially expressed in ovaries and heads of 1-day-old fifth-instar larvae to 6-day pupae by quantitative real-time reverse transcription polymerase chain reaction (qRT-PCR). In line with published data, the expression level of these genes in larvae were generally lower than in pupae. We further analyzed the expression levels of the four genes in BmN cells that had been treated with various concentrations of diapause hormone (DH). It demonstrated that the expression of these genes was affected by DH. Knockdown of the selected genes in non-diapause-destined female pupae changed the fate of the progeny from non-diapause- to daipause-destined, as seen by the appearance of diapause eggs. Our study provides insight into the molecular mechanism of diapause in B. mori.

Molecular cloning of heat shock protein 10 (Hsp10) and 60 (Hsp60) cDNAs from Galeruca daurica (Coleoptera: Chrysomelidae) and their expression analysis

Galeruca daurica (Joannis) is a new outbreak pest in the Inner Mongolia grasslands in northern China. Heat shock protein 10 and 60 (Hsp10 and Hsp60) genes of G. daurica, designated as GdHsp10 and GdHsp60, were cloned by rapid amplification of cDNA ends techniques. Sequence analysis showed that GdHsp10 and GdHsp60 encoded polypeptides of 104 and 573 amino acids, respectively. Sequence alignment and phylogenetic analysis clearly revealed that the amino acids of GdHsp10 and GdHsp60 had high homology and were clustered with other Hsp10 and Hsp60 genes in insects which are highly relative with G. daurica based on morphologic taxonomy. The mRNA expression analysis by real-time PCR revealed that GdHsp10 and GdHsp60 were expressed at all development stages and in all tissues examined, but expressed highest in eggs and in adults' abdomen; both heat and cold stresses could induce mRNA expression of GdHsp10 and GdHsp60 in the 2nd instar larvae; the two Hsp genes were expressed from high to low with the extension of treatment time in G. daurica eggs exposed to freezing point. Overall, our study provides useful information to understand temperature stress responses of Hsp60 and Hsp10 in G. daurica, and provides a basis to further study functions of Hsp60/Hsp10 relative to thermotolerance and cold hardiness mechanism.

Insect heat shock proteins during stress and diapause

Insect heat shock proteins include ATP-independent small heat shock proteins and the larger ATP-dependent proteins, Hsp70, Hsp90, and Hsp60. In concert with cochaperones and accessory proteins, heat shock proteins mediate essential activities such as protein folding, localization, and degradation. Heat shock proteins are synthesized constitutively in insects and induced by stressors such as heat, cold, crowding, and anoxia. Synthesis depends on the physiological state of the insect, but the common function of heat shock proteins, often working in networks, is to maintain cell homeostasis through interaction with substrate proteins. Stress-induced expression of heat shock protein genes occurs in a background of protein synthesis inhibition, but in the course of diapause, a state of dormancy and increased stress tolerance, these genes undergo differential regulation without the general disruption of protein production. During diapause, when ATP concentrations are low, heat shock proteins may sequester rather than fold proteins.

Cloning of cDNA encoding a Bombyx mori homolog of human oxidation resistance 1 (OXR1) protein from diapause eggs, and analyses of its expression and function

To better understand the molecular mechanisms of diapause initiation, we used the sensitive cDNA subtraction (selective amplification via biotin- and restriction-mediated enrichment) method and isolated a novel gene expressed abundantly in diapause eggs of the silkworm, Bombyx mori, which encodes a homolog of the human oxidation resistance 1 (OXR1) protein. Quantitative real-time polymerase chain reaction (qRT-PCR) and Western blotting analyses confirmed that BmOXR1 mRNA and its 140-kDa protein were differentially expressed in diapause eggs compared to non-diapause eggs. OXR1 double-stranded RNA (dsRNA) was injected into diapause-destined eggs before the cellular blastoderm stage, and 4days later, when untreated eggs reached the diapause stage, the OXR1 protein disappeared; however, these eggs remained in diapause, suggesting that BmOXR1 is not essential for diapause initiation and/or maintenance. To further investigate the in vivo function of BmOXR1 apart from its role in diapause, we overexpressed BmOXR1 in Drosophila melanogaster. The fruit fly male adult life-span was significantly extended in the 50%-survival time when adults were reared on diets both with and without H2O2 solution under 25°C incubation. These results suggest that BmOXR1 functions in D. melanogaster via a possible antioxidant effect. As BmOXR1 was expressed mainly in the nuclei of D. melanogaster cells, the mechanism underlying its antioxidation effect appears to be different from that in humans where it is expressed mainly in the mitochondria. Taken together, these results suggest that BmOXR1 might serve as an antioxidant regulator during the early diapause stage.

Diapause, signal and molecular characteristics of overwintering Chilo suppressalis (Insecta: Lepidoptera: Pyralidae)

Diapause is a complex and dynamic process. Chilo suppressalis, an important rice pest in Asia enters facultative diapause as larvae. Our results demonstrated in Yangzhou, China, diapause was initiated between September 4 and 12, 2010. After diapause termination, C. suppressalis remained in quiescence in the field for as long as three months. The average time between collection of field larvae of C. suppressalis and their pupation decreased as the season progressed from fall to next spring. Unexpectedly, the pupated ratio of female to male in the initiation of diapause was 0.22. The abundance of hsp90, hsp70, hsp60 and CsAQP1 all peaked on January 8 or 15, 2011. Nitric oxide (NO) is a secondary messenger that is positively correlated with the diapause of C. suppressalis. Among several geographically separated populations of C. suppressalis, there are no significant differences in the mRNA levels of hsp70, hsp60 or CsAQP1.

Global transcriptome analysis of orange wheat blossom midge, Sitodiplosis mosellana (Gehin) (Diptera: Cecidomyiidae) to identify candidate transcripts regulating diapause

BACKGROUND

Many insects enter a developmental arrest (diapause) that allows them to survive harsh seasonal conditions. Despite the well-established ecological significance of diapause, the molecular basis of this crucial adaptation remains largely unresolved. Sitodiplosis mosellana (Gehin), the orange wheat blossom midge (OWBM), causes serious damage to wheat throughout the northern hemisphere, and sporadic outbreaks occur in the world. Traits related to diapause appear to be important factors contributing to their rapid spread and outbreak. To better understand the diapause mechanisms of OWBM, we sequenced the transcriptome and determined the gene expression profile of this species.

METHODOLOGY/PRINCIPAL FINDINGS

In this study, we performed de novo transcriptome analysis using short-read sequencing technology (Illumina) and gene expression analysis with a tag-based digital gene expression (DGE) system. The sequencing results generated 89,117 contigs, and 45,713 unigenes. These unigenes were annotated by Blastx alignment against the NCBI non-redundant (nr), Clusters of orthologous groups (COG), gene orthology (GO), and the Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. 20,802 unigenes (45.5% of the total) matched with protein in the NCBI nr database. Two digital gene expression (DGE) libraries were constructed to determine differences in gene expression profiles during diapause and non-diapause developmental stages. Genes related to diapause were analyzed in detail and in addition, nine diapause-related genes were analyzed by real time PCR.

CONCLUSIONS/SIGNIFICANCE

The OWBM transcriptome greatly improves our genetic understanding and provides a platform for functional genomics research of this species. The DGE profiling data provides comprehensive information at the transcriptional level that facilitates our understanding of the molecular mechanisms of various physiological aspects including development and diapause stages in OWBM. From this study it is evident that various genes coding metabolic enzymes are crucial for diapause and metamorphosis.

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.

Transcript profiling reveals mechanisms for lipid conservation during diapause in the mosquito, Aedes albopictus

The Asian tiger mosquito, Aedes albopictus, is a medically important invasive species whose geographic distribution has expanded dramatically during the past 20 years, and one of the key elements of its success is its capacity to survive long distance transport as a diapausing pharate first instar larva, encased within the chorion of the egg. We report that pharate larvae entering diapause are larger and contain 30% more lipid than their nondiapausing counterparts. To improve our understanding of the molecular regulation of lipid metabolism during diapause, we assessed the relative mRNA abundance of 21 genes using qRT-PCR. Elevated expression of lipid storage droplet protein 2 during embryonic development likely contributes to the higher amounts of lipid we noted in diapausing individuals. The conservation of lipids during diapause is reflected in downregulation of genes involved in lipid catabolism, including lipase 2, lipase 3, lipase 4, acyl-CoA dehydrogenase 4, and isovaleryl-CoA dehydrogenase. Two genes involved in fatty acid synthesis and modification, Δ(9)-desaturase, and fatty acyl-CoA elongase, were both upregulated in diapausing pharate larvae, suggesting roles for their gene products in generating unsaturated fatty acids to enhance membrane fluidity at low temperatures and generating precursors to the surface hydrocarbons needed to resist desiccation, respectively. Together, the results point to substantial distinctions in lipid metabolism within the embryo as a consequence of the diapause program, and these differences occur both before the actual onset of diapause as well as during the diapause state.