Response to cold acclimation in diapause pupae of Hyles euphorbiae (Lepidoptera: Sphingidae): candidate biomarker identification using proteomics

Genetic constraints in genes exhibiting splicing plasticity in facultative diapause

Phenotypic plasticity is produced and maintained by processes regulating the transcriptome. While differential gene expression is among the most important of these processes, relatively little is known about other sources of transcriptional variation. Previous work suggests that alternative splicing plays an extensive and functionally unique role in transcriptional plasticity, though plastically spliced genes may be more constrained than the remainder of expressed genes. In this study, we explore the relationship between expression and splicing plasticity, along with the genetic diversity in those genes, in an ecologically consequential polyphenism: facultative diapause. Using 96 samples spread over two tissues and 10 timepoints, we compare the extent of differential splicing and expression between diapausing and direct developing pupae of the butterfly Pieris napi. Splicing differs strongly between diapausing and direct developing trajectories but alters a smaller and functionally unique set of genes compared to differential expression. We further test the hypothesis that among these expressed loci, plastically spliced genes are likely to experience the strongest purifying selection to maintain seasonally plastic phenotypes. Genes with unique transcriptional changes through diapause consistently had the lowest nucleotide diversity, and this effect was consistently stronger among genes that were differentially spliced compared to those with just differential expression through diapause. Further, the strength of negative selection was higher in the population expressing diapause every generation. Our results suggest that maintenance of the molecular mechanisms involved in diapause progression, including post-transcriptional modifications, are highly conserved and likely to experience genetic constraints, especially in northern populations of P. napi.

Transcriptional Regulation of Reproductive Diapause in the Convergent Lady Beetle, Hippodamia convergens

Simple Summary

Diapause is a dormant period typically controlled by daylength that ensures an insect’s survival through harsh environmental conditions. The convergent lady beetle, Hippodamia convergens, undergoes a reproductive diapause in winter, where female ovaries remain immature and no eggs are laid. This species is an important biological control agent, but during diapause, beetles are less likely to eat pest insects. Thus, knowledge of diapause mechanisms may facilitate manipulation thereof to improve biological control. Further, molecular studies of adult diapause and diapause in Coleoptera are relatively lacking. Here, we assembled and annotated a transcriptome for this species and quantified transcript expression changes during diapause. Female beetles were sampled at three times in diapause (early, mid, and late diapause), which allowed us to characterize the molecular processes occurring at distinct transitions throughout diapause. We found that transcripts involved in flight were consistently upregulated during diapause, which is consistent with dispersal flights at this stage, while transcripts involved in ovarian development were downregulated, which is consistent with the shutdown of reproduction in diapausing females. These findings identify key regulators of diapause in H. convergens and contribute to a growing body of literature on the molecular mechanisms of diapause across the insect phylogeny.

Abstract

Diapause is an alternate development program that synchronizes an insect’s life cycle with seasonally abundant resources and ensures survival in unfavorable conditions. The physiological basis of diapause has been well characterized, but the molecular mechanisms regulating it are still being elucidated. Here, we present a de novo transcriptome and quantify transcript expression during diapause in the convergent lady beetle Hippodamia convergens. H. convergens is used as an augmentative biocontrol agent, and adult females undergo reproductive diapause that is regulated by photoperiod. We sampled females at three stages (early, mid, and late diapause) and compared transcript expression to non-diapausing individuals. Based on principle component analysis, the transcriptomes of diapausing beetles were distinct from non-diapausing beetles, and the three diapausing points tended to cluster together. However, there were still classes of transcripts that differed in expression across distinct phases of diapause. In general, transcripts involved in muscle function and flight were upregulated during diapause, likely to support dispersal flights that occur during diapause, while transcripts involved in ovarian development were downregulated. This information could be used to improve biological control by manipulating diapause. Additionally, our data contribute to a growing understanding of the genetic regulation of diapause across diverse insects.

The Importance of ATM and ATR in Physcomitrella patens DNA Damage Repair, Development, and Gene Targeting

Coordinated by ataxia-telangiectasia-mutated (ATM) and ATM and Rad3-related (ATR), two highly conserved kinases, DNA damage repair ensures genome integrity and survival in all organisms. The Arabidopsis thaliana (A. thaliana) orthologues are well characterized and exhibit typical mammalian characteristics. We mutated the Physcomitrellapatens (P. patens) PpATM and PpATR genes by deleting functionally important domains using gene targeting. Both mutants showed growth abnormalities, indicating that these genes, particularly PpATR, are important for normal vegetative development. ATR was also required for repair of both direct and replication-coupled double-strand breaks (DSBs) and dominated the transcriptional response to direct DSBs, whereas ATM was far less important, as shown by assays assessing resistance to DSB induction and SuperSAGE-based transcriptomics focused on DNA damage repair genes. These characteristics differed significantly from the A. thaliana genes but resembled those in yeast (Saccharomyces cerevisiae). PpATR was not important for gene targeting, pointing to differences in the regulation of gene targeting and direct DSB repair. Our analysis suggests that ATM and ATR functions can be substantially diverged between plants. The differences in ATM and ATR reflect the differences in DSB repair pathway choices between A. thaliana and P. patens, suggesting that they represent adaptations to different demands for the maintenance of genome stability.

Functional characterization of the Hyles euphorbiae hawkmoth transcriptome reveals strong expression of phorbol ester detoxification and seasonal cold hardiness genes

Background

The European spurge hawkmoth, Hyles euphorbiae (Lepidoptera, Sphingidae), has been intensively studied as a model organism for insect chemical ecology, cold hardiness and evolution of species delineation. To understand species isolation mechanisms at a molecular level, this study aims at determining genetic factors underlying two adaptive ecological trait candidates, phorbol ester (TPA) detoxification and seasonal cold acclimation.

Method

A draft transcriptome of H. euphorbiae was generated using Illumina sequencing, providing the first genomic resource for the hawkmoth subfamily Macroglossinae. RNA expression levels in tissues of experimental TPA feeding larvae and cooled pupae was compared to levels in control larvae and pupae using 26 bp RNA sequence tag libraries (DeepSuperSAGE). Differential gene expression was assessed by homology searches of the tags in the transcriptome.

Results

In total, 389 and 605 differentially expressed transcripts for detoxification and cold hardiness, respectively, could be identified and annotated with proteins. The majority (22 of 28) of differentially expressed detox transcripts of the four 'drug metabolism' enzyme groups (cytochrome P450 (CYP), carboxylesterases (CES), glutathione S-transferases (GST) and lipases) are up-regulated. Triacylglycerol lipase was significantly over proportionally annotated among up-regulated detox transcripts. We record several up-regulated lipases, GSTe2, two CESs, CYP9A21, CYP6BD6 and CYP9A17 as candidate genes for further H. euphorbiae TPA detoxification analyses. Differential gene expression of the cold acclimation treatment is marked by metabolic depression with enriched Gene Ontology terms among down-regulated transcripts almost exclusively comprising metabolism, aerobic respiration and dissimilative functions. Down-regulated transcripts include energy expensive respiratory proteins like NADH dehydrogenase, cytochrome oxidase and ATP synthase. Gene expression patterns show shifts in carbohydrate metabolism towards cryoprotectant production. The Glycolysis enzymes, G1Pase, A1e, Gpi and an Akr isoform are up-regulated. Glycerol, an osmolyte which lowers the body liquid supercooling point, appears to be the predominant polyol cryoprotectant in H. euphorbiae diapause pupae. Several protein candidates involved in glucose, glycerol, myo-inositol and potentially sorbitol and trehalose synthesis were identified.

Conclusions

A majority of differently expressed transcripts unique for either detoxification or cold hardiness indicates highly specialized functional adaptation which may have evolved from general cell metabolism and stress response.The transcriptome and extracted candidate biomarkers provide a basis for further gene expression studies of physiological processes and adaptive traits in H. euphorbiae.

Differences in protein expression among five species of stream stonefly (Plecoptera) along a latitudinal gradient in Japan

Proteome variation among natural populations along an environmental gradient may provide insights into how the biological functions of species are related to their local adaptation. We investigated protein expression in five stream stonefly species from four geographic regions along a latitudinal gradient in Japan with varying climatic conditions. The extracted proteins were separated by two-dimensional gel electrophoresis and identified by matrix-assisted laser desorption/ionization of time-of-flight (MALDI TOF/TOF), yielding 446 proteins. Low interspecies variation in the proteome profiles was observed among five species within geographical regions, presumably due to the co-occurring species sharing the environments. However, large spatial variations in protein expression were found among four geographic regions, suggesting strong regulation of protein expression in heterogeneous environments, where the spatial variations were positively correlated with water temperature. We identified 21 unique proteins expressed specifically in a geographical region and six common proteins expressed throughout all regions. In warmer regions, metabolic proteins were upregulated, whereas proteins related to cold stress, the photoperiod, and mating were downregulated. Oxygen-related and energy-production proteins were upregulated in colder regions with higher altitudes. Thus, our proteomic approach is useful for identifying and understanding important biological functions related to local adaptations by populations of stoneflies.

Arabidopsis glucosinolates trigger a contrasting transcriptomic response in a generalist and a specialist herbivore

Phytophagous insects have to deal with toxic defense compounds from their host plants. Although it is known that insects have evolved genes and mechanisms to detoxify plant allochemicals, how specialist and generalist precisely respond to specific secondary metabolites at the molecular level is less understood. Here we studied the larval performance and transcriptome of the generalist moth Heliothis virescens and the specialist butterfly Pieris brassicae feeding on Arabidopsis thaliana genotypes with different glucosinolate (GS) levels. H. virescens larvae gained significantly more weight on the GS-deficient mutant quadGS compared to wild-type (Col-0) plants. On the contrary, P. brassicae was unaffected by the presence of GS and performed equally well on both genotypes. Strikingly, there was a considerable differential gene expression in H. virescens larvae feeding on Col-0 compared to quadGS. In contrast, compared to H. virescens, P. brassicae displayed a much-reduced transcriptional activation when fed on both plant genotypes. Transcripts coding for putative detoxification enzymes were significantly upregulated in H. virescens, along with digestive enzymes and transposable elements. These data provide an unprecedented view on transcriptional changes that are specifically activated by GS and illustrate differential molecular responses that are linked to adaptation to diet in lepidopteran herbivores.

Transcriptome and proteome analyses reveal complex mechanisms of reproductive diapause in the two-spotted spider mite, Tetranychus urticae

Although a variety of factors underlying diapause have been identified in arthropods and other organisms, the molecular mechanisms regulating diapause are still largely unknown. Here, to better understand this process, we examined diapause-associated genes in the two-spotted spider mite, Tetranychus urticae, by comparing the transcriptomes and proteomes of early diapausing and reproductive adult females. Amongst genes underlying diapause revealed by the transcriptomic and proteomic data sets, we described the noticeable change in Ca²⁺ -associated genes, including 65 Ca²⁺ -binding protein genes and 23 Ca²⁺ transporter genes, indicating that Ca²⁺ signalling has a substantial role in diapause regulation. Other interesting changes in diapause included up-regulation of (1) glutamate receptors that may be involved in synaptic plasticity changes, (2) genes involved in cytoskeletal reorganization including genes encoding each of the components of thick and thin filaments, tubulin and members of integrin signalling and (3) genes involved in anaerobic energy metabolism, which reflects a shift to anaerobic energy metabolism in early diapausing mites.

The Genetic Basis of Natural Variation in Drosophila (Diptera: Drosophilidae) Virgin Egg Retention

Drosophila melanogaster is able to thrive in harsh northern climates through adaptations in life-history traits and physiological mechanisms that allow for survival through the winter. We examined the genetic basis of natural variation in one such trait, female virgin egg retention, which was previously shown to vary clinally and seasonally. To further our understanding of the genetic basis and evolution of virgin egg retention, we performed a genome-wide association study (GWAS) using the previously sequenced Drosophila Genetic Reference Panel (DGRP) mapping population. We found 29 single nucleotide polymorphisms (SNPs) associated with virgin egg retention and assayed 6 available mutant lines, each harboring a mutation in a candidate gene, for effects on egg retention time. We found that four out of the six mutant lines had defects in egg retention time as compared with the respective controls: mun, T48, Mes-4, and Klp67A Surprisingly, none of these genes has a recognized role in ovulation control, but three of the four genes have known effects on fertility or have high expression in the ovaries. We also found that the SNP set associated with egg retention time was enriched for clinal SNPs. The majority of clinal SNPs had alleles associated with longer egg retention present at higher frequencies in higher latitudes. Our results support previous studies that show higher frequency of long retention times at higher latitude, providing evidence for the adaptive value of virgin egg-retention.

Seasonal phenotype-specific transcriptional reprogramming during metamorphosis in the European map butterfly Araschnia levana

The European map butterfly (Araschnia levana) is a classic example of seasonal polyphenism because the spring and summer imagoes display two distinct morphological phenotypes. The light regime and temperature during larval and prepupal development determine whether or not the pupae commit to diapause and overwintering and thus whether spring or summer imagoes emerge. We used suppression subtractive hybridization to experimentally screen for genes that are differentially expressed in prepupae committed either to accelerated metamorphosis and egg production or diapause and overwintering. The range and ontology of the differentially expressed genes in prepupae developing from larvae exposed either to long‐day (LD) or short‐day (SD) conditions revealed fundamental differences. The SD prepupae preferentially expressed genes related to cuticle formation and immunity, reflecting the formation of a robust pupal exoskeleton and the upregulation of antimicrobial peptides as preparations for overwintering. One protein preferentially expressed in SD prepupae has a counterpart in Bombyx mori that functions as a diapause duration clock. The differentially expressed genes in LD prepupae included several members of the dusky and osiris families. We also observed the strong induction of different yellow‐like genes under SD and LD conditions which suggest a role in the developmental choice between seasonal phenotypes. Our transcriptomic data will facilitate the more detailed analysis of molecular mechanisms underlying seasonal polyphenism.

PROTEOMICS ANALYSIS OF OVEREXPRESSED PLASMA PROTEINS IN RESPONSE TO COLD ACCLIMATION IN Ostrinia furnacalis

Many insects in temperate regions overwinter in diapause. In these insects, one of the metabolic adaptations to cold stress is the synthesis of responsive proteins. Using proteomic analysis, an investigation aimed to a better understanding of the molecular adaptation mechanisms to cold stress was carried out in Ostrinia furnacalis larva. Proteins were extracted from the larval hemolymph collected from both control and overwintering larva. By polyethylene glycol precipitation, approximately 560 protein spots were separated and visualized on two-dimensional (2D) gels after silver staining. Eighteen protein spots were found to be upregulated in overwinter larval plasma in different patterns. As an initial work, 13 of these proteins were identified using MALDI TOF/TOF MS. The differentially overexpressed proteins include heat shock 70 kDa cognate protein, small heat shock protein (sHSP), putative aliphatic nitrilase, arginine kinase, phosphoglyceromutase, triosephosphateisomerase, and glutathione transferase. Alterations in the levels of these proteins were further confirmed by qPCR. This study is the first analysis of differentially expressed plasma proteins in O. furnacalis diapause larvae under extremely low temperature conditions and gives new insights into the acclimation mechanisms responsive to cold stress. Our results also support the idea that energy metabolism, alanine and proline metabolism, and antioxidative reaction act in the cold acclimation of O. furnacalis diapause larvae.

Comparative proteomic analysis of Bactrocera dorsalis (Hendel) in response to thermal stress

Temperature is one of the most important environmental variables affecting growth, reproduction and distribution of insects. The rise of comparative proteomics provides a powerful tool to explore the response in proteins to thermal stress. As an important worldwide pest, the oriental fruit fly Bactrocera dorsalis causes severe economic losses to crops. To understand the response of B. dorsalis to thermal stress, we performed a comparative proteome analysis of this insect after exposure to extreme low and high temperatures using two-dimensional electrophoresis. Among the separated proteins, 51 diverse protein spots were present differently in response to extreme temperatures. Using tandem mass spectrometry sequencing analysis 39 proteins were successfully identified, which included 13 oxidoreductases, 10 binding proteins, 5 transferases, and 2 each of lyases, isomerases, ligases, and developmental proteins. Subsequently, the expression of these protein transcripts was studied by RT-qPCR to validate the proteomic results. In conclusion, this study provides a first look into the thermal stress response of B. dorsalis at the protein level, and thus it paves the way for further functional studies in the physiological mechanism related to thermal stress.