Novel isoforms of heat shock transcription factor 1 are induced by hypoxia in the Pacific oyster Crassostrea gigas

Interactive abiotic and biotic stressor impacts on a stream-dwelling amphibian

Organisms within freshwater and marine environments are subject to a diverse range of often co-occurring abiotic and biotic stressors. Despite growing awareness of the complex multistress systems at play in aquatic ecosystems, many questions remain regarding how simultaneous stressors interact with one another and jointly impact aquatic species. We looked at multistress interactions in a protected stream ecosystem in Mendocino County, California. Specifically, we examined how diurnal temperature variation, turbidity, and predator cues altered the movement speed of larval Pacific giant salamanders (Dicamptodon tenebrosus). In a second experiment, we looked at how simulated low-flow summer conditions impact the expression of heat-shock proteins (HSPs) in the same species. Larvae moved almost one and a half times faster in the presence of chemical cues from trout and suspended sediment, and almost two times faster when both sediment and trout cues were present but were only marginally affected by temperature and visual cues from conspecifics. Interestingly, the order of stressor exposure also appeared to influence larval speed, where exposure to sediment and trout in earlier trials tended to lead to faster speeds in later trials. Additionally, larvae exposed to low-flow conditions had more variable, but not statistically significantly higher, expression of HSPs. Our findings highlight the potential interactive effects of an abiotic stressor, sedimentation, and a biotic stressor, and predator chemical cues on an ecologically important trait: movement speed. Our findings also demonstrate the likely role of HSPs in larval salamander survival in challenging summer conditions. Taken together, these findings show that larval D. tenebrosus responds behaviorally to biotic and abiotic stressors and suggests a possible pathway for physiological tolerance of environmental stress. Consideration of multistress systems and their effects is important for understanding the full effects of co-occurring stressors on aquatic organisms to guide appropriate conservation and management efforts based on ecologically relevant responses of organisms within an environment.

Contribution of HIF-1α to Heat Shock Response by Transcriptional Regulation of HSF1/HSP70 Signaling Pathway in Pacific Oyster, Crassostrea gigas

Ocean temperature rising drastically threatens the adaptation and survival of marine organisms, causing serious ecological impacts and economic losses. It is crucial to understand the adaptive mechanisms of marine organisms in response to high temperature. In this study, a novel regulatory mechanism that is mediated by hypoxia-inducible factor-1α (HIF-1α) was revealed in Pacific oyster (Crassostrea gigas) in response to heat stress. We identified a total of six HIF-1α genes in the C. gigas genome, of which HIF-1α and HIF-1α-like5 were highly induced under heat stress. We found that the HIF-1α and HIF-1α-like5 genes played critical roles in the heat shock response (HSR) through upregulating the expression of heat shock protein (HSP). Knocking down of HIF-1α via RNA interference (RNAi) inhibited the expression of heat shock factor 1 (HSF1) and HSP70 genes in C. gigas under heat stress. Both HIF-1α and HIF-1α-like5 promoted the transcriptional activity of HSF1 by binding to hypoxia response elements (HREs) within the promoter region. Furthermore, the survival of C. gigas under heat stress was significantly decreased after knocking down of HIF-1α. This work for the first time revealed the involvement of HIF-1α/HSF1/HSP70 pathway in response to heat stress in the oyster and provided an insight into adaptive mechanism of bivalves in the face of ocean warming.

Alternative Splicing of Heat Shock Transcription Factor 2 Regulates the Expression of Laccase Gene Family in Response to Copper in Trametes trogii

White-rot fungi, especially Trametes strains, are the primary source of industrial laccases in bioenergy and bioremediation. Trametes strains express members of the laccase gene family with different physicochemical properties and expression patterns. However, the literature on the expression pattern of the laccase gene family in T. trogii S0301 and the response mechanism to Cu²⁺, a key laccase inducer, in white-rot fungal strains is scarce. In the present study, we found that Cu²⁺ could induce the mRNAs and proteins of the two alternative splicing variants of heat shock transcription factor 2 (TtHSF2). Furthermore, the overexpression of alternative splicing variants TtHSF2α and TtHSF2β-I in the homokaryotic T. trogii S0301 strain showed opposite effects on the extracellular total laccase activity, with the maximum laccase activity of approximately 0.6 U mL^-1 and 3.0 U mL^-1, respectively, on the eighth day, which is 0.4 and 2.3 times that of the wild type strain. Similarly, TtHSF2α and TtHSF2β-I play opposite roles in the oxidation tolerance to H₂O₂ In addition, the direct binding of TtHSF2α to the promoter regions of the representative laccase isoenzymes (TtLac1 and TtLac13) and protein-protein interactions between TtHSF2α and TtHSF2β-I were detected. Our results demonstrate the crucial roles of TtHSF2 and its alternative splicing variants in response to Cu²⁺ We believe that these findings will deepen our understanding of alternative splicing of HSFs and their regulatory mechanism of the laccase gene family in white-rot fungi.Importance The members of laccase gene family in Trametes strains are the primary source of industrial laccase and have gained widespread attention. Increasing the yield and enzymatic properties of laccase through various methods has always been a topic worthy of attention, and there is no report on the regulation of laccase expression through HSF transcription factor engineering. Here, we found that two alternative splicing variants of TtHSF2 functioned oppositely in regulating the expression of laccase genes, and copper can induce the expression of almost all members of the laccase gene family. Most importantly, our study suggested that TtHSF2 and its alternative splicing variants are vital for copper-induced production of laccases in T. trogii S0301.

Regulation Between HSF1 Isoforms and HSPs Contributes to the Variation in Thermal Tolerance Between Two Oyster Congeners

Heat shock transcription factor 1 (HSF1) plays an important role in regulating heat shock, which can activate heat shock proteins (HSPs). HSPs can protect organisms from thermal stress. Oysters in the intertidal zone can tolerate thermal stress. The Pacific oyster (Crassostrea gigas gigas) and Fujian oyster (C. gigas angulata)—allopatric subspecies with distinct thermal tolerances—make good study specimens for analyzing and comparing thermal stress regulation. We cloned and compared HSF1 isoforms, which is highly expressed under heat shock conditions in the two subspecies. The results revealed that two isoforms (HSF1a and HSF1d) respond to heat shock in both Pacific and Fujian oysters, and different heat shock conditions led to various combinations of isoforms. Subcellular localization showed that isoforms gathered in the nucleus when exposed to heat shock. The co-immunoprecipitation revealed that HSF1d can be a dimer. In addition, we selected HSPs that are expressed under the heat shock response, according to the RNA-seq and proteomic analyses. For the HSPs, we analyzed the coding part and the promoter sequences. The result showed that the domains of HSPs are conserved in two subspecies, but the promoters are significantly different. The Dual-Luciferase assay showed that the induced expression isoform HSF1d had the highest activity in C. gigas gigas, while the constitutively-expressed HSF1a was most active in C. gigas angulata. In addition, variation in the level of HSP promoters appeared to be correlated with gene expression. We argue that this gene is regulated based on the different expression levels between the two subspecies’ responses to heat shock. In summary, various stress conditions can yield different HSF1 isoforms and respond to heat shock in both oyster subspecies. Differences in how the isoforms and promoter are activated may contribute to their differential expressions. Overall, the results comparing C. gigas gigas and C. gigas angulata suggest that these isoforms have a regulatory relationship under heat shock, providing valuable information on the thermal tolerance mechanism in these commercially important oyster species.

Identification of HSF1 Target Genes Involved in Thermal Stress in the Pacific Oyster Crassostrea gigas by ChIP-seq

The Pacific oyster Crassostrea gigas, a commercially important species inhabiting the intertidal zone, facing enormous temperature fluctuations. Therefore, it is important to identify candidate genes and key regulatory relationships associated with thermal tolerance, which can aid the molecular breeding of oysters. Heat shock transcription factor 1 (HSF1) plays an important role in the thermal stress resistance. However, the regulatory relationship between the expansion of heat shock protein (HSP) HSP 70 and HSF1 is not yet clear in C. gigas. In this study, we analyzed genes regulated by HSF1 in response to heat shock by chromatin immunoprecipitation followed by sequencing (ChIP-seq), determined the expression patterns of target genes by qRT-PCR, and validated the regulatory relationship between one HSP70 and HSF1. We found 916 peaks corresponding to HSF1 binding sites, and these peaks were annotated to the nearest genes. In Gene Ontology analysis, HSF1 target genes were related to signal transduction, energy production, and response to biotic stimulus. Four HSP70 genes, two HSP40 genes, and one small HSP gene exhibited binding to HSF1. One HSP70 with a binding site in the promoter region was validated to be regulated by HSF1 under heat shock. These results provide a basis for future studies aimed at determining the mechanisms underlying thermal tolerance and provide insights into gene regulation in the Pacific oyster.

RNAi based transcriptome suggests genes potentially regulated by HSF1 in the Pacific oyster Crassostrea gigas under thermal stress

Background

The Pacific oyster Crassostrea gigas is an important fishery resource that is sensitive to temperature fluctuations. Thus, it has evolved a protection mechanism against heat stress by increasing the expression of the gene coding for heat shock protein (HSP) 70 under elevated temperatures. In other animals, heat shock response is a transcriptional response driven by the heat shock transcription factor 1 (HSF1) and thermal stress can trigger HSP70 expression to protect the organism via HSF1. However, the regulatory relationship between HSF1 and HSP remains unclear in Pacific oyster. Therefore, in the present study, we examined the transcriptomic response of several to thermal stress following HSF1 interference.

Results

We identified 150 genes responsive to heat shock including seven HSP genes, six of which belonging to the group of 17 HSP genes enriched in response to heat shock, according to weighted gene co-expression network analysis (WGCNA). The other gene was enriched in the module correlated with HSF1 interference. In addition, we found 48 and 47 genes that were upregulated and downregulated by HSF1 in response to heat shock, respectively. In the upregulated genes, we identified one HSP70 potentially regulated by HSF1 in response to heat shock. Furthermore, based on differentially expressed genes and WGCNA analyses, we found that the hypoxia signaling pathway was enriched under heat shock conditions. Five genes were then selected to detect dynamic changes through time. The results suggested that gene expression was correlated with HSF1 expression. The regulation of HSP70 by HSF1 was preliminarily confirmed by binding site predictions and by a dual luciferase assay.

Conclusions

Our results revealed that the expression of HSP70 and HSP20 was initially triggered after 2 h of heat shock, and one of the HSP70 genes was potentially regulated by HSF1. From these results, it is evident that not all heat-inducible genes were triggered simultaneously in response to heat shock stress. Overall, the results revealed a possible HSF1–HSP regulatory relationship in Pacific oyster, providing valuable information on the mechanisms of thermal tolerance in this commercially important oyster.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-6003-8) contains supplementary material, which is available to authorized users.

Analysis of transcriptional responses of normalizing genes on Crassostrea brasiliana under different experimental conditions

Bivalves show remarkable plasticity to environmental changes and have been proposed as sentinel organisms in biomonitoring. Studies related to transcriptional analysis using quantitative real-time polymerase chain reaction (qRT-PCR) in these organisms have notably increased, imposing a need to identify and validate adequate reference genes for an accurate and reliable analysis. In the present study, 9 reference genes were selected from transcriptome data of Crassostrea brasiliana to identify their suitability as qRT-PCR normalizer genes. The transcriptional patterns were analyzed in gills of oysters under 3 different conditions: different temperatures (18, 24, or 32 °C) and phenanthrene (100 µg L^-1 ) combined exposure; different salinities (10, 25, or 35‰) and phenanthrene combined exposure; and 10% of diesel fuel water-accommodated fraction (diesel-WAF) exposure. Reference gene stability was calculated using 5 algorithms (geNorm, NormFinder, BestKeeper, ΔCt, RefFinder). Transcripts of ankyrin-like (ANK), glyceraldehyde 3-phosphate dehydrogenase-like (GAPDH), and α-tubulin-like (TUBA) genes showed minor changes in different temperature/phenanthrene treatment. Transcripts of ANK, β-actin-like, and β-tubulin-like genes showed better stability at salinity/phenanthrene treatment, and ANK, TUBA, and 28S ribosomal protein-like genes showed the most stable transcription pattern in oysters exposed to diesel-WAF exposure. The present study constitutes the first systematic analysis of reference gene selection for qRT-PCR normalization in C. brasiliana. These genes could be employed in studies using qRT-PCR analysis under similar experimental conditions. Environ Toxicol Chem 2017;36:2190-2198. © 2017 SETAC.

The novel transcriptional factor HP1BP3 negatively regulates Hsp70 transcription in Crassostrea hongkongensis

ChHP1BP3, a chromatin complex-related protein known with dynamic features, was identified as a ChHsp70 promoter-associated factor in Crassostrea hongkongensis by DNA-affinity purification and mass spectrometry analysis. Direct interaction between purified ChHP1BP3 and the ChHsp70 promoter region was demonstrated using EMSA. ChHp1bp3 depletion led to clear enhancements in ChHsp70 mRNA expression in C. hongkongensis hemocytes. However, ChHp1bp3 overexpression in heterologous HEK293T cells correlated with fluctuations in ChHsp70 transcription. Quantitative RT-PCR analysis showed that both ChHsp70 and ChHp1bp3 transcription were responsive to external physical/chemical stresses by heat, CdCl₂ and NP. This indicated a plausible correlation between ChHsp70 and ChHp1bp3 in the stress-induced genetic regulatory pathway. While, the distinctive ChHp1bp3 expression patterns upon physical and chemical stresses suggest that the mechanisms that mediate ChHp1bp3 induction might be stress-specific. This study discovered a novel role for HP1BP3 as a negative regulator in controlling Hsp70 transcription in C. hongkongensis, and contributed to better understanding the complex regulatory mechanisms governing Hsp70 transcription.

Transcriptomic evaluation of the American oyster, Crassostrea virginica, deployed during the Deepwater Horizon oil spill: Evidence of an active hydrocarbon response pathway

Estuarine organisms were impacted by the Deepwater Horizon oil spill which released ∼5 million barrels of crude oil into the Gulf of Mexico in the spring and summer of 2010. Crassostrea virginica, the American oyster, is a keystone species in these coastal estuaries and is routinely used for environmental monitoring purposes. However, very little is known about their cellular and molecular responses to hydrocarbon exposure. In response to the spill, a monitoring program was initiated by deploying hatchery-reared oysters at three sites along the Alabama and Mississippi coast (Grand Bay, MS, Fort Morgan, AL, and Orange Beach, AL). Oysters were deployed for 2-month periods at five different time points from May 2010 to May 2011. Gill and digestive gland tissues were harvested for gene expression analysis and determination of aliphatic and polycyclic aromatic hydrocarbon (PAH) concentrations. To facilitate identification of stress response genes that may be involved in the hydrocarbon response, a nearly complete transcriptome was assembled using Roche 454 and Illumina high-throughput sequencing from RNA samples obtained from the gill and digestive gland tissues of deployed oysters. This effort resulted in the assembly and annotation of 27,227 transcripts comprised of a large assortment of stress response genes, including members of the aryl hydrocarbon receptor (AHR) pathway, Phase I and II biotransformation enzymes, antioxidant enzymes and xenobiotic transporters. From this assembly several potential biomarkers of hydrocarbon exposure were chosen for expression profiling, including the AHR, two cytochrome P450 1A genes (CYP1A-like 1 and CYP1A-like 2), Cu/Zn superoxide dismutase (CuZnSOD), glutathione S-transferase theta (GST theta) and multidrug resistance protein 3 (MRP3). Higher expression levels of GST theta and MRP3 were observed in gill tissues from all three sites during the summer to early fall 2010 deployments. Linear regression analysis indicated a statistically significant relationship between total PAH levels in digestive gland tissue samples with CYP1A-like 2, CuZnSOD, GST theta and MRP3 induction. These observations provide evidence of a potentially conserved AHR pathway in invertebrates and yield new insight into the development of novel biomarkers for use in environmental monitoring activities.

High mobility group protein DSP1 negatively regulates HSP70 transcription in Crassostrea hongkongensis

HSP70 acts mostly as a molecular chaperone and plays important roles in facilitating the folding of nascent peptides as well as the refolding or degradation of the denatured proteins. Under stressed conditions, the expression level of HSP70 is upregulated significantly and rapidly, as is known to be achieved by various regulatory factors controlling the transcriptional level. In this study, a high mobility group protein DSP1 was identified by DNA-affinity purification from the nuclear extracts of Crassostrea hongkongensis using the ChHSP70 promoter as a bait. The specific interaction between the prokaryotically expressed ChDSP1 and the FITC-labeled ChHSP70 promoter was confirmed by EMSA analysis. ChDSP1 was shown to negatively regulate ChHSP70 promoter expression by Luciferase Reporter Assay in the heterologous HEK293T cells. Both ChHSP70 and ChDSP1 transcriptions were induced by either thermal or CdCl2 stress, while the accumulated expression peaks of ChDSP1 were always slightly delayed when compared with that of ChHSP70. This indicates that ChDSP1 is involved, very likely to exert its suppressive role, in the recovery of the ChHSP70 expression from the induced level to its original state. This study is the first to report negative regulator of HSP70 gene transcription, and provides novel insights into the mechanisms controlling heat shock protein expression.

Molecular Basis for Adaptation of Oysters to Stressful Marine Intertidal Environments

Oysters that occupy estuarine and intertidal habitats have well-developed stress tolerance mechanisms to tolerate harsh and dynamically changing environments. In this review, we summarize common pathways and genomic features in oyster that are responsive to environmental stressors such as temperature, salinity, hypoxia, air exposure, pathogens, and anthropogenic pollutions. We first introduce the key genes involved in several pathways, which constitute the molecular basis for adaptation to stress. We use genome analysis to highlight the strong cellular homeostasis system, a unique adaptive characteristic of oysters. Next, we provide a global view of features of the oyster genome that contribute to stress adaptation, including oyster-specific gene expansion, highly inducible expression, and functional divergence. Finally, we review the consequences of interactions between oysters and the environment from ecological and evolutionary perspectives by discussing mass mortality and adaptive divergence among populations and related species of the genus Crassostrea. We conclude with prospects for future study.

Seasonal changes in gene expression and polymorphism of hsp70 in cultivated oysters (Crassostrea gigas) at extreme temperatures

The HSP70 proteins are an important element of the response against thermal stress and infectious diseases, and they are highly conserved and ubiquitous. In some species, variations on the hsp70 encoding sequence resulted in intraspecific differential expression, which leads to variations on thermo-tolerance among individuals. This phenomenon has not been described in the Pacific oyster Crassostrea gigas, which is cultivated in Mexico under temperature conditions highly above the optimal for this species. The present study was aimed to identify associations between hsp70 genotypes and their expression levels in C. gigas. By analyzing a 603 bp fragment from the 3' end of the hsp70 gene, 21 different genotypes with 60 nucleotide polymorphic sites were detected, of which 34 sites were found in heterozygous condition. Although no correlation was found between genotype-expression-season, a minimum expression threshold that should be taken into account as an important feature for a future breeding program is proposed.

Novel isoforms of heat shock transcription factor 1, HSF1γα and HSF1γβ, regulate chaperone protein gene transcription

The heat shock response, resulting in the production of heat shock proteins or molecular chaperones, is triggered by elevated temperature and a variety of other stressors. Its master regulator is heat shock transcription factor 1 (HSF1). Heat shock factors generally exist in multiple isoforms. The two known isoforms of HSF1 differ in the inclusion (HSF1α) or exclusion (HSF1β) of exon 11. Although there are some data concerning the differential expression patterns and transcriptional activities of HSF2 isoforms during development, little is known about the distinct properties of the HSF1 isoforms. Here we present evidence for two novel HSF1 isoforms termed HSF1γα and HSF1γβ, and we show that the HSF1 isoform ratio differentially regulates heat shock protein gene transcription. Hsf1γ isoforms are expressed in various mouse tissues and are translated into protein. Furthermore, after heat shock, HSF1γ isoforms are exported from the nucleus more rapidly or degraded more quickly than HSF1α or HSF1β. We also show that each individual HSF1 isoform is sufficient to induce the heat shock response and that expression of combinations of HSF1 isoforms, in particular HSF1α and HSF1β, results in a synergistic enhancement of the transcriptional response. In addition, HSF1γ isoforms potentially suppress the synergistic effect of HSF1α and HSF1β co-expression. Collectively, our observations suggest that the expression of HSF1 isoforms in a specific ratio provides an additional layer in the regulation of heat shock protein gene transcription.

Validation of housekeeping genes as internal controls for studying gene expression during Pacific oyster (Crassostrea gigas) development by quantitative real-time PCR

Hatchery-reared larvae of the Pacific oyster (Crassostrea gigas) often suffer from massive mortality induced by Ostreid herpesvirus 1 (OsHV-1) infection, indicating the importance of better understanding of oyster immune defense systems. The accuracy of measurements of gene expression levels based on quantitative real-time PCR assays relies on the use of housekeeping genes as internal controls; however, few studies have focused on the selection of such internal controls. In this study, we conducted a comprehensive investigation of internal control genes during oyster development in virus-infected and uninfected samples. Transcriptome data for 38 developmental stages were downloaded and the gene expression patterns were classified into 30 clusters. A total of 317 orthologs of classical housekeeping genes in the oyster genome were annotated. After combining the expression profiles and oyster housekeeping gene dataset, 14 candidate internal controls were selected for further investigation: Elongation factor-1α (EF-1α), 18S rRNA (18S), 28S rRNA (28S), Glyceraldehyde 3-phosphate dehydrogenase (GAPDH), β-actin (ACT), Ribosomal protein L7 (RL7), Ribosomal protein L27 (RL27), Ribosomal protein L36 (RL36), Ribosomal protein S18 (RS18), Heterogeneous nuclear ribonucleoprotein A2/B1 (RO21), Eukaryotic translation elongation factor 2 (EF2), Ubiquitin-conjugating enzyme E2D2 (UBCD1), S-phase kinase-associated protein 1 (SKP1) and Heterogeneous nuclear ribonucleoprotein Q (HNRPQ). RNA was extracted from oyster larvae infected with OsHV-1 (group A; GA), and OsHV-1 free larvae (group B; GB). The expression levels of the 14 candidate internal controls were studied in GA and GB larvae by real-time PCR. Their expression stabilities were further analyzed using the GeNorm program. RL7 and RS18 were the most stable genes in both OsHV-1 infected (GA) and uninfected (GB) larvae. These results suggest that RL7 and RS18 could be used as internal controls for studying gene expression in normal growing oyster larvae and in OsHV-1 infected larvae. These high quality internal controls will be a valuable resource in future studies of oyster larval mortality.

Regulation of a truncated isoform of AMP-activated protein kinase α (AMPKα) in response to hypoxia in the muscle of Pacific oyster Crassostrea gigas

AMP-activated protein kinase α (AMPKα) is a key regulator of energy balance in many model species during hypoxia. In a marine bivalve, the Pacific oyster Crassostrea gigas, we analyzed the protein content of adductor muscle in response to hypoxia during 6 h. In both smooth and striated muscles, the amount of full-length AMP-activated protein kinase α (AMPKα) remained unchanged during hypoxia. However, hypoxia induced a rapid and muscle-specific response concerning truncated isoforms of AMPKα. In the smooth muscle, a truncated isoform of AMPKα was increased from 1 to 6 h of hypoxia, and was linked with accumulation of AKT kinase, a key enzyme of the insulin signaling pathway which controls intracellular glucose metabolism. In this muscle, aerobic metabolism was maintained over the 6 h of hypoxia, as mitochondrial citrate synthase activity remained constant. In contrast, in striated muscle, hypoxia did not induce any significant modification of neither truncated AMPKα nor AKT protein content, and citrate synthase activity was altered after 6 h of hypoxia. Together, our results demonstrate that hypoxia response is specific to muscle type in Pacific oyster, and that truncated AMPKα and AKT proteins might be involved in maintaining aerobic metabolism in smooth muscle. Such regulation might occur in vivo during tidal intervals that cause up to 6 h of hypoxia.