Genetics of gene expression responses to temperature stress in a sea urchin gene network

High temperature influences DNA methylation and transcriptional profiles in sea urchins (Strongylocentrotus intermedius)

BACKGROUND

DNA methylation plays an important role in life processes by affecting gene expression, but it is still unclear how DNA methylation is controlled and how it regulates gene transcription under high temperature stress conditions in Strongylocentrotus intermedius. The potential link between DNA methylation variation and gene expression changes in response to heat stress in S. intermedius was investigated by MethylRAD-seq and RNA-seq analysis. We screened DNA methylation driver genes in order to comprehensively elucidate the regulatory mechanism of its high temperature adaptation at the DNA/RNA level.

RESULTS

The results revealed that high temperature stress significantly affected not only the DNA methylation and transcriptome levels of S. intermedius (P < 0.05), but also growth. MethylRAD-seq analysis revealed 12,129 CG differential methylation sites and 966 CWG differential methylation sites, and identified a total of 189 differentially CG methylated genes and 148 differentially CWG methylated genes. Based on KEGG enrichment analysis, differentially expressed genes (DEGs) are mostly enriched in energy and cell division, immune, and neurological damage pathways. Further RNA-seq analysis identified a total of 1968 DEGs, of which 813 genes were upregulated and 1155 genes were downregulated. Based on the joint MethylRAD-seq and RNA-seq analysis, metabolic processes such as glycosaminoglycan degradation, oxidative phosphorylation, apoptosis, glutathione metabolism, thermogenesis, and lysosomes are regulated by DNA methylation.

CONCLUSIONS

High temperature affected the DNA methylation and expression levels of genes such as MOAP-1, GGT1 and RDH8, which in turn affects the metabolism of HPSE, Cox, glutathione, and retinol, thereby suppressing the immune, energy metabolism, and antioxidant functions of the organism and finally manifesting as stunted growth. In summary, the observations in the present study improve our understanding of the molecular mechanism of the response to high temperature stress in sea urchin.

Marine heatwave temperatures enhance larval performance but are meditated by paternal thermal history and inter-individual differences in the purple sea urchin, Strongylocentrotus purpuratus

Marine heatwave (MHW) events, characterized by periods of anomalous temperatures, are an increasingly prevalent threat to coastal marine ecosystems. Given the seasonal phenology of MHWs, the full extent of their biological consequences may depend on how these thermal stress events align with an organism's reproductive cycle. In organisms with more complex life cycles (e.g., many marine invertebrate species) the alignment of adult and larval environments may be an important factor determining offspring success, setting the stage for MHW events to influence reproduction and development in situ. Here, the influence of MHW-like temperatures on the early development of the California purple sea urchin, Strongylocentrotus purpuratus, were explored within the context of paternal thermal history. Based on temperature data collected during MHW events seen in Southern California from 2014-2020, adult urchins were acclimated to either MHW or non-MHW temperatures for 28 days before their sperm was used to produce embryos that were subsequently raised under varying thermal conditions. Once offspring reached an early larval stage, the impact of paternal and offspring environments were assessed on two aspects of offspring performance: larval size and thermal tolerance. Exposure to elevated temperatures during early development resulted in larger, more thermally tolerant larvae, with further influences of paternal identity and thermal history, respectively. The alignment of paternal and offspring exposure to MHW temperatures had additional positive benefits on larval thermal tolerance, but this tolerance significantly decreased when their thermal experience mismatched. As the highest recorded temperatures within past MHW events have occurred during the gametogenesis of many kelp forest benthic marine invertebrate species, such as the purple sea urchin, such parental mediated impacts may represent important drivers of future recruitment and population composition for these species.

Without a pinch of salt: effect of low salinity on eggs and nauplii of the salmon louse (Lepeophtheirus salmonis)

The salmon louse is an economically important parasite on Atlantic salmon and poses a major threat to aquaculture. Several treatment methods have lost their effect due to resistance development in the lice. A rather new method for combatting sea lice is freshwater treatment where the various life stages of lice are differently affected by this treatment. In this study, we analyzed the effect of freshwater on the egg strings. A 3-h treatment with freshwater had a detrimental effect on the egg strings. First, the water penetrated the string, widening it, then entering the eggs and enlarging them. Finally, the ordered structure of the egg strings collapsed, and no alive animals hatched. Shorter treatments had a lower effectivity, and treatments with brackish water also showed milder effects. The egg strings were found to have a protective effect against low salinities, as hatched nauplii died rapidly under conditions that embryos survived. We also found that embryos react to low salinity on a molecular level by changing gene expression of several genes, when incubated in brackish water. Additionally, the hatching of embryos treated with brackish water was delayed in comparison to seawater controls.

How thermal challenges change gene regulation in the songbird brain and gonad: implications for sexual selection in our changing world

In a rapidly warming world, exposure to high temperatures may impact fitness, but the gene regulatory mechanisms that link sublethal heat to sexually selected traits are not well understood, particularly in endothermic animals. Our experiment used zebra finches (Taeniopygia guttata), songbirds that experience extreme temperature fluctuations in their native Australia. We exposed captive males to an acute thermal challenge (43°C) compared with thermoneutral (35°C) and lower (27°C) temperatures. We found significantly more heat dissipation behaviors at 43°C, a temperature previously shown to reduce song production and fertility, and more heat retention behaviors at 27°C. Next, we characterized transcriptomic responses in tissues important for mating effort - the posterior telencephalon, for its role in song production, and the testis, for its role in fertility and hormone production. Differential expression of hundreds of genes in the testes, but few in the brain, suggest the brain is less responsive to extreme temperatures. Nevertheless, gene network analyses revealed that expression related to dopaminergic signaling in the brain co-varied with heat dissipation behaviors, providing a mechanism by which temporary thermal challenges may alter motivational circuits for song production. In both brain and testis, we observed correlations between thermally sensitive gene networks and individual differences in thermoregulatory behavior. Although we cannot directly relate these gene regulatory changes to mating success, our results suggest that individual variation in response to thermal challenges could impact sexually selected traits in a warming world.

Gene expression patterns of sea urchins (Strongylocentrotus intermedius) exposed to different combinations of temperature and hypoxia

Strongylocentrotus intermedius is one of the most economically valuable sea urchin species in China, and its growth and survival are severely constrained by ocean warming and the hypoxia that often accompanies high water temperatures. To elucidate the molecular mechanisms of S. intermedius that regulate gene expression in response to multi-causal environmental stresses. We performed a de novo transcriptome analysis of coelomocyte from S. intermedius to heat (25 °C), hypoxia (2 mg/L), and the combined stress. We identified 35,635, 29,107, and 29,440 differentially expressed genes (DEGs) in S. intermedius cultured under high temperature, low oxygen, and combined stress, respectively. Further Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyses revealed that these DEGs mainly enriched the functional categories of "Protein processing in endoplasmic reticulum," and "Glutathione metabolism" by heat stress, such as HSP70, GSTO1, PDIA4. After hypoxic stress, "Notch signaling pathway" and metabolism-related pathways such as "Glycerolipid metabolism", "Pyruvate metabolism" were significantly enriched. Exposure to combined stress resulted in a two-factor additive effect at the transcriptome level and have a more extensive impact on the immune correlated pathways in S. intermedius than single stress, the expression of related immune genes (C3, C5, and AIFM2) were up-regulated. Quantitative real-time PCR (qRT-PCR) analysis of the expression of 18 DEGs confirmed the RNA-Seq results. Observations in the present study will improve the understanding of the molecular mechanism of S. intermedius in response to multi-causal environmental stress.

Large-scale genome-wide study reveals climate adaptive variability in a cosmopolitan pest

Understanding the genetic basis of climatic adaptation is essential for predicting species' responses to climate change. However, intraspecific variation of these responses arising from local adaptation remains ambiguous for most species. Here, we analyze genomic data from diamondback moth (Plutella xylostella) collected from 75 sites spanning six continents to reveal that climate-associated adaptive variation exhibits a roughly latitudinal pattern. By developing an eco-genetic index that combines genetic variation and physiological responses, we predict that most P. xylostella populations have high tolerance to projected future climates. Using genome editing, a key gene, PxCad, emerged from our analysis as functionally temperature responsive. Our results demonstrate that P. xylostella is largely capable of tolerating future climates in most of the world and will remain a global pest beyond 2050. This work improves our understanding of adaptive variation along environmental gradients, and advances pest forecasting by highlighting the genetic basis for local climate adaptation.

PAHs and PCBs Affect Functionally Intercorrelated Genes in the Sea Urchin Paracentrotus lividus Embryos

Polycyclic aromatic hydrocarbons (PAHs) and polychlorinated biphenyls (PCBs) represent the most common pollutants in the marine sediments. Previous investigations demonstrated short-term sublethal effects of sediments polluted with both contaminants on the sea urchin Paracentrotus lividus after 2 months of exposure in mesocosms. In particular, morphological malformations observed in P. lividus embryos deriving from adults exposed to PAHs and PCBs were explained at molecular levels by de novo transcriptome assembly and real-time qPCR, leading to the identification of several differentially expressed genes involved in key physiological processes. Here, we extensively explored the genes involved in the response of the sea urchin P. lividus to PAHs and PCBs. Firstly, 25 new genes were identified and interactomic analysis revealed that they were functionally connected among them and to several genes previously defined as molecular targets of response to the two pollutants under analysis. The expression levels of these 25 genes were followed by Real Time qPCR, showing that almost all genes analyzed were affected by PAHs and PCBs. These findings represent an important further step in defining the impacts of slight concentrations of such contaminants on sea urchins and, more in general, on marine biota, increasing our knowledge of molecular targets involved in responses to environmental stressors.

Resurvey of sea urchins and mussels at protected and harvested shores a decade after: A beyond-BACI approach

Comparing temporal patterns of distribution and abundance of target organisms between protected and harvested shores is essential to assess the extant effectiveness of marine protected areas (MPAs) and whether it is maintained through time. By means of an adapted Beyond-BACI approach, we compared the short- and long-term patterns of variation in the abundance of the sea urchin Paracentrotus lividus and the mussel Mytilus galloprovincialis at a protected shore (within the Parque Litoral Norte MPA, Portugal) and at three adjacent shores subject to intense harvesting over a decadal interval. Despite the existence of the MPA for more than 30 years, we did not obtain clear evidence of its persistent or recent effectiveness on intertidal species of commercial interest. We suggest the need for refining management options along the northern Portuguese coast, possibly by better enforcing current regulations and reconsidering the present design of protection schemes. Moreover, the adopted analytical approach may represent a methodological reference for similar investigations in systems where the perturbation of interest (protection or disturbance) would not occur at a given time during the course of the study, but has been operating since before the first sampling occasion and maintained until subsequent surveys.

Comparative CpG methylation kinetic patterns of cis-regulatory regions of heat stress-related genes in Sahiwal and Frieswal cattle upon persistent heat stress

The kinetic patterns of CpG methylation of the cis-regulatory region of heat stress-related genes on exposed to heat stress (at 42 °C) between the Sahiwal and Frieswal cattle was compared in the present study. Using an in vitro whole blood culture model, cells were continuously exposed to heat stress (at 42 °C) for 6 h. Methylation levels of five genes, viz., GPX1, HSP70, HSP90, c-FOS, and JUN were estimated by SyberGreen-based quantitative methylation-specific PCR (qMSP) assay. CpG methylation kinetics at different time points of heat stress (0.5, 1, 2, 4, 6 h) were analyzed using mixed ANOVA. The initial methylation level, estimated at 37 °C, of HSP70 was significantly high in the Sahiwal breed. A significant (p<0.001) time-dependent hypomethylation of an antioxidant gene (GPX1) CpG islands was detected at the acute phase of the stress. Heat shock protein gene (HSP70) showed a similar CpG methylation kinetics where the hypomethylation was prominent from 1 h and persisted up to 4 h. The heat stress responses of both Sahiwal and Frieswal cattle were identical as there was no distinctiveness in the methylation kinetics of CpG islands of studied genes. The acclimatization of Frieswal cattle-a breed developed in India over the years to the tropical climatic conditions, maybe one of the reasons for this similarity. Thus, the present study results could pave a path to understand the molecular mechanism of heat stress and adaptation of indigenous and crossbred cattle populations to the changing scenario in tropical climate conditions.

Gene expression patterns of red sea urchins (Mesocentrotus franciscanus) exposed to different combinations of temperature and pCO2 during early development

BACKGROUND

The red sea urchin Mesocentrotus franciscanus is an ecologically important kelp forest herbivore and an economically valuable wild fishery species. To examine how M. franciscanus responds to its environment on a molecular level, differences in gene expression patterns were observed in embryos raised under combinations of two temperatures (13 °C or 17 °C) and two pCO2 levels (475 μatm or 1050 μatm). These combinations mimic various present-day conditions measured during and between upwelling events in the highly dynamic California Current System with the exception of the 17 °C and 1050 μatm combination, which does not currently occur. However, as ocean warming and acidification continues, warmer temperatures and higher pCO2 conditions are expected to increase in frequency and to occur simultaneously. The transcriptomic responses of the embryos were assessed at two developmental stages (gastrula and prism) in light of previously described plasticity in body size and thermotolerance under these temperature and pCO2 treatments.

RESULTS

Although transcriptomic patterns primarily varied by developmental stage, there were pronounced differences in gene expression as a result of the treatment conditions. Temperature and pCO2 treatments led to the differential expression of genes related to the cellular stress response, transmembrane transport, metabolic processes, and the regulation of gene expression. At each developmental stage, temperature contributed significantly to the observed variance in gene expression, which was also correlated to the phenotypic attributes of the embryos. On the other hand, the transcriptomic response to pCO2 was relatively muted, particularly at the prism stage.

CONCLUSIONS

M. franciscanus exhibited transcriptomic plasticity under different temperatures, indicating their capacity for a molecular-level response that may facilitate red sea urchins facing ocean warming as climate change continues. In contrast, the lack of a robust transcriptomic response, in combination with observations of decreased body size, under elevated pCO2 levels suggest that this species may be negatively affected by ocean acidification. High present-day pCO2 conditions that occur due to coastal upwelling may already be influencing populations of M. franciscanus.

Different amplitudes of temperature fluctuation induce distinct transcriptomic and metabolomic responses in the dung beetle Phanaeus vindex

Most studies exploring molecular and physiological responses to temperature have focused on constant temperature treatments. To gain a better understanding of the impact of fluctuating temperatures, we investigated the effects of increased temperature variation on Phanaeus vindex dung beetles across levels of biological organization. Specifically, we hypothesized that increased temperature variation is energetically demanding. We predicted that thermal sensitivity of metabolic rate and energetic reserves would be reduced with increasing fluctuation. To test this, we examined the responses of dung beetles to constant (20°C), low fluctuation (20±5°C), or high fluctuation (20±12°C) temperature treatments using respirometry, assessment of energetic reserves and HPLC-MS-based metabolomics. We found no significant differences in metabolic rate or energetic reserves, suggesting increased fluctuations were not energetically demanding. To understand why there was no effect of increased amplitude of temperature fluctuation on energetics, we assembled and annotated a de novo transcriptome, finding non-overlapping transcriptomic and metabolomic responses of beetles exposed to different fluctuations. We found that 58 metabolites increased in abundance in both fluctuation treatments, but 15 only did so in response to high-amplitude fluctuations. We found that 120 transcripts were significantly upregulated following acclimation to any fluctuation, but 174 were upregulated only in beetles from the high-amplitude fluctuation treatment. Several differentially expressed transcripts were associated with post-translational modifications to histones that support a more open chromatin structure. Our results demonstrate that acclimation to different temperature fluctuations is distinct and may be supported by increasing transcriptional plasticity. Our results indicate for the first time that histone modifications may underlie rapid acclimation to temperature variation.

Post-ejaculation thermal stress causes changes to the RNA profile of sperm in an external fertilizer

Sperm cells experience considerable post-ejaculation environmental variation. However, little is known about whether this affects their molecular composition, probably owing to the assumption that sperm are transcriptionally quiescent. Nevertheless, recent evidence shows sperm have distinct RNA profiles that affect fertilization and embryo viability. Moreover, RNAs are expected to be highly sensitive to extracellular changes. One such group of RNAs are heat shock protein (hsp) transcripts, which function in stress responses and are enriched in sperm. Here, we exploit the experimental tractability of the mussel Mytilus galloprovincialis by exposing paired samples of ejaculated sperm to ambient (19°C) and increased (25°C) temperatures, then measure (i) sperm motility phenotypes, and (ii) messenger RNA (mRNA) levels of two target genes (hsp70 and hsp90) and several putative reference genes. We find no phenotypic changes in motility, but reduced mRNA levels for hsp90 and the putative reference gene gapdh at 25°C. This could reflect either decay of specific RNAs, or changes in translation and degradation rates of transcripts to maintain sperm function under stress. These findings represent, to our knowledge, the first evidence for changes in sperm RNA profiles owing to post-ejaculation environments, and suggest that sperm may be more vulnerable to stress from rising temperatures than currently thought.

Transcriptomes shed light on transgenerational and developmental effects of ocean warming on embryos of the sea urchin Strongylocentrotus intermedius

Ocean warming increasingly endangers the fitness of marine invertebrates. Transgenerational effects (TE) potentially mitigate the impacts of environmental stress on the embryos of marine invertebrates. The molecular mechanisms, however, remain largely unknown. Using high-throughput RNA sequencing technology, we investigated the gene expression patterns of embryos (the gastrula stage) of the sea urchin Strongylocentrotus intermedius at different developmental temperatures, whose parents were exposed to long-term (15 months) elevated temperature (A) or not (B). The temperatures at which adults were held for ~4 weeks prior to the start of the experiment (21 °C for A and 18 °C for B) were also used for the development of offspring (high: 21 °C and ambient (laboratory): 18 °C) resulting in four experimental groups (HA and HB at 21 °C, and LA and LB at 18 °C). The embryos were sampled ~24 h after fertilization. All samples were in the gastrula stage. Twelve mRNA libraries (groups HA, HB, LA, LB, 3 replicates for each group) were established for the following sequencing. Embryos whose parents were exposed to elevated temperatures or not showed 1891 significantly different DEGs (differentially expressed genes) at the ambient developmental temperature (LB vs LA, LB as control) and 2203 significantly different DEGs at the high developmental temperature (HB vs HA, HB as control), respectively. This result indicates complex molecular mechanisms of transgenerational effects of ocean warming, in which a large number of genes are involved. With the TE, we found 904 shared DEGs in both LB vs LA (LB as control) and HB vs HA (HB as control) changed in the same direction of expression (i.e., up- or down-regulated), indicating that parental exposed temperatures affect the expression of these genes in the same manner regardless of the development temperature. With developmental exposure, we found 198 shared DEGs in both HB vs LB (HB as control) and HA vs LA (HA as control) changed in the same direction of expression. Of the 198 DEGs, more genes were up-regulated at high developmental temperature. Interestingly, embryos whose parents were exposed to high temperature showed fewer differently expressed DEGs between high and low developmental temperatures than the individuals whose parents were exposed to ambient temperature. The results indicate that gene expressions are probably depressed by the transgenerational effect of ocean warming. The roles of hsp70 and hnf6 in thermal acclimation are highlighted for future studies. The present study provides new insights into the molecular mechanisms of the transgenerational and developmental effects of ocean warming on the embryos of sea urchins.

Cooler temperatures slow the repair of DNA damage in tadpoles exposed to ultraviolet radiation: Implications for amphibian declines at high altitude

Ultraviolet B radiation (UVBR) damages the DNA of exposed cells, causing dimers to form between adjacent pyrimidine nucleotides. These dimers block DNA replication, causing mutations and apoptosis. Most organisms utilize biochemical or biophysical DNA repair strategies to restore DNA structure; however, as with most biological reactions, these processes are likely to be thermally sensitive. Tadpoles exposed to elevated UVBR at low environmental temperatures have significantly higher rates of mortality and developmental deformities compared with tadpoles exposed to the same levels of UVBR at higher environmental temperatures. We hypothesized that low environmental temperatures impair the primary enzymatic (photolyase) DNA repair pathway in amphibians, leading to the accumulation of DNA damage. To test this hypothesis, we compared DNA repair rates and photolyase gene expression patterns in Limnodynastes peronii. Tadpoles were acutely exposed to UVBR for 1 hr at either 20 or 30°C, and we measured DNA damage and photolyase expression levels at intervals following this exposure. Temperature had a significant effect on the rate of DNA repair, with repair at 30°C occurring twice as fast as repair at 20°C. Photolyase gene expression (6-4 PP and CPD) was significantly upregulated by UVBR exposure, with expression levels increasing within 6 hr of UVBR exposure. CPD expression levels were not significantly affected by temperature, but 6-4 PP expression was significantly higher in tadpoles in the 30°C treatment within 12 hr of UVBR exposure. These data support the hypothesis that DNA repair rates are thermally sensitive in tadpoles and may explain why enigmatic amphibian declines are higher in montane regions where UVBR levels are naturally elevated and environmental temperatures are lower.

Morphological and molecular responses of the sea urchin Paracentrotus lividus to highly contaminated marine sediments: The case study of Bagnoli-Coroglio brownfield (Mediterranean Sea)

Marine sediments store complex mixtures of compounds, including heavy metals, organotins and a large array of other contaminants. Sediment quality monitoring, characterization and management are priorities, due to potential impacts of the above compounds on coastal waters and their biota, especially in cases of pollutants released during dredging activities. Harbours and marinas, as well as estuaries and bays, where limited exchanges of water occurr, the accumulation of toxic compounds poses major concerns for human and environmental health. Here we report the effects of highly contaminated sediments from the site of national interest Bagnoli-Coroglio (Tyrrhenian Sea, Western Mediterranean) on the sea urchin Paracentrotus lividus, considered a good model for ecotoxicological studies. Adult sea urchins were reared one month in aquaria in the presence of contaminated sediment that was experimentally subject to different patterns of re-suspension events (mimicking the effect of natural storms occurring in the field), crossed with O₂ enrichment versus natural gas exchanges in the water. The development of embryos deriving from adult urchins exposed to such experimental conditions was followed until the pluteus stage, checking the power of contaminated sediment to induce morphological malformations and its eventual buffering by high oxygenation. Real-Time qPCR analysis revealed that the expression of several genes (among the fifty analyzed, involved in different functional processes) was targeted by contaminated sediments more than those exposed in oxygen-enriched condition. Our findings have biological and ecological relevance in terms of assessing the actual impact on local organisms of chronic environmental contamination by heavy metals and polycyclic aromatic hydrocarbons affecting the Bagnoli-Coroglio area, and of exploring enhanced sediment and water oxygenation as a promising tool to mitigate the effects of contamination in future environmental restoration actions.

Constitutive and Plastic Gene Expression Variation Associated with Desiccation Resistance Differences in the Drosophila americana Species Group

Stress response mechanisms are ubiquitous and important for adaptation to heterogenous environments and could be based on constitutive or plastic responses to environmental stressors. Here we quantify constitutive and plastic gene expression differences under ambient and desiccation stress treatments, in males and females of three species of Drosophila known to differ in desiccation resistance. Drosophila novamexicana survives desiccation trials significantly longer than the two subspecies of Drosophila americana, consistent with its natural species range in the desert southwest USA. We found that desiccation stress reduces global expression differences between species—likely because many general stress response mechanisms are shared among species—but that all species showed plastic expression changes at hundreds of loci during desiccation. Nonetheless, D. novamexicana had the fewest genes with significant plastic expression changes, despite having the highest desiccation resistance. Of the genes that were significantly differentially expressed between species—either within each treatment (>200 loci), constitutively regardless of treatment (36 loci), or with different species-specific plasticity (26 loci)—GO analysis did not find significant enrichment of any major gene pathways or broader functions associated with desiccation stress. Taken together, these data indicate that if gene expression changes contribute to differential desiccation resistance between species, these differences are likely shaped by a relatively small set of influential genes rather than broad genome-wide differentiation in stress response mechanisms. Finally, among the set of genes with the greatest between-species plasticity, we identified an interesting set of immune-response genes with consistent but opposing reaction norms between sexes, whose potential functional role in sex-specific mechanisms of desiccation resistance remains to be determined.

Genomics of Developmental Plasticity in Animals

Developmental plasticity refers to the property by which the same genotype produces distinct phenotypes depending on the environmental conditions under which development takes place. By allowing organisms to produce phenotypes adjusted to the conditions that adults will experience, developmental plasticity can provide the means to cope with environmental heterogeneity. Developmental plasticity can be adaptive and its evolution can be shaped by natural selection. It has also been suggested that developmental plasticity can facilitate adaptation and promote diversification. Here, we summarize current knowledge on the evolution of plasticity and on the impact of plasticity on adaptive evolution, and we identify recent advances and important open questions about the genomics of developmental plasticity in animals. We give special attention to studies using transcriptomics to identify genes whose expression changes across developmental environments and studies using genetic mapping to identify loci that contribute to variation in plasticity and can fuel its evolution.

The Impact of Chronic Heat Stress on the Growth, Survival, Feeding, and Differential Gene Expression in the Sea Urchin Strongylocentrotus intermedius

To explore the impact of chronic heat stress on commercial echinoderms, the present study assessed the effects of chronic high temperature on the growth, survival, feeding, and differential gene expression in the sea urchin Strongylocentrotus intermedius cultured in northern Yellow Sea in China. One suitable seawater condition (20°C) and one laboratory-controlled high temperature condition (25°C) were set up. After 28 days incubation, our results showed that: (1) The specific growth, survival, and ingestion rates of S. intermedius reared under high temperature (25°C) decreased compared to those reared under optimal temperature (20°C) conditions; (2) comparative transcriptome analysis identified 2,125 differentially expressed genes (DEGs) in S. intermedius reared under high temperature (25°C) compared to those subjected to optimal temperature condition (20°C), which included 1,015 upregulated and 1,100 downregulated genes. The accuracy of the transcriptome profiles was verified by quantitative real-time PCR (qRT-PCR). Further Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways analyses revealed that these DEGs mainly enriched the functional categories of ribosome, protein processing in endoplasmic reticulum, and prion diseases. A total of 732 temperature-induced expressed genes, such as ATP5, heat shock protein 70, and heat shock protein 90, were identified as candidates that were closely correlated with heat resistance in S. intermedius. Differentially expressed transcription factors (TFs), such as AP-1, Fos, CREB, and ZNF, were also identified as potential regulators that regulate the molecular network that was associated with responses to heat stress in sea urchins. Observations in the present study provide additional information that improves our understanding of the molecular mechanism of temperate echinoid species in response to heat stress, as well as theoretical basis for the molecular-assisted breeding of heat-resistant sea urchins.

Aberrant gene expression profiles in Mediterranean sea urchin reproductive tissues after metal exposures

Marine organisms are simultaneously exposed to numerous pollutants, among which metals probably represent the most abundant in marine environments. In order to evaluate the effects of metal exposure at molecular level in reproductive tissues, we profiled the sea urchin transcriptional response after non-lethal exposures using pathway-focused mRNA expression analyses. Herein, we show that exposures to relatively high concentrations of both essential and toxic metals hugely affected the gonadic expression of several genes involved in stress-response, detoxification, transcriptional and post-transcriptional regulation, without significant changes in gonadosomatic indices. Even though treatments did not result in reproductive tissues visible alterations, metal exposures negatively affected the main mechanisms of stress-response, detoxification and survival of adult P. lividus. Additionally, transcriptional changes observed in P. lividus gonads may cause altered gametogenesis and maintenance of heritable aberrant epigenetic effects. This study leads to the conclusion that exposures to metals, as usually occurs in polluted coastal areas, may affect sea urchin gametogenesis, thus supporting the hypothesis that parental exposure to environmental stressors affects the phenotype of the offspring.

Toxicity of diatom-derived polyunsaturated aldehyde mixtures on sea urchin Paracentrotus lividus development

Diatom-derived polyunsaturated aldehydes (PUAs), decadienal, heptadienal and octadienal, derive from the oxidation of fatty acids and have cytotoxic and anticancer effects. PUAs, tested separately, induce malformations in sea urchin Paracentrotus lividus embryos. Decadienal induces the worst malformations and lowest survival rates. Interestingly, decadienal, heptadienal and octadienal place in motion several genes to counteract their negative effects. To date, no studies are available reporting on the effects of PUA mixtures on marine invertebrates. Here we test binary and ternary mixtures on embryonic development of P. lividus. Our findings demonstrate that mixtures of PUAs act (i) at morphological level in synergistic way, being much more severe compared to individual PUAs; (ii) at molecular level also reveal an additive effect, affecting almost all fifty genes, previously tested using individual PUAs. This study is relevant from an ecological point of view since diatoms are a major food source for both pelagic and benthic organisms. This work opens new perspectives for understanding the molecular mechanisms that marine organisms use in reacting to environmental natural toxin mixtures such as diatom PUAs.

Stressful conditions reveal decrease in size, modification of shape but relatively stable asymmetry in bumblebee wings

Human activities can generate a wide variety of direct and indirect effects on animals, which can manifest as environmental and genetic stressors. Several phenotypic markers have been proposed as indicators of these stressful conditions but have displayed contrasting results, depending, among others, on the phenotypic trait measured. Knowing the worldwide decline of multiple bumblebee species, it is important to understand these stressors and link them with the drivers of decline. We assessed the impact of several stressors (i.e. natural toxin-, parasite-, thermic- and inbreeding- stress) on both wing shape and size and their variability as well as their directional and fluctuating asymmetries. The total data set includes 650 individuals of Bombus terrestris (Hymenoptera: Apidae). Overall wing size and shape were affected by all the tested stressors. Except for the sinigrin (e.g. glucosinolate) stress, each stress implies a decrease of wing size. Size variance was affected by several stressors, contrary to shape variance that was affected by none of them. Although wing size directional and fluctuating asymmetries were significantly affected by sinigrin, parasites and high temperatures, neither directional nor fluctuating shape asymmetry was significantly affected by any tested stressor. Parasites and high temperatures led to the strongest phenotype modifications. Overall size and shape were the most sensitive morphological traits, which contrasts with the common view that fluctuating asymmetry is the major phenotypic marker of stress.

Comparative Studies of Gene Expression Kinetics: Methodologies and Insights on Development and Evolution

Across the animal kingdom, embryos of closely related species show high morphological similarity despite genetic and environmental distances. Deciphering the molecular mechanisms that underlie morphological conservation and those that support embryonic adaptation are keys to understand developmental robustness and evolution. Comparative studies of developmental gene regulatory networks can track the genetic changes that lead to evolutionary novelties. However, these studies are limited to a relatively small set of genes and demand extensive experimental efforts. An alternative approach enabled by next-generation sequencing, is to compare the expression kinetic of large sets of genes between different species. The advantages of these comparisons are that they can be done relatively easily, for any species and they provide information of all expressed genes. The challenge in these experiments is to compare the kinetic profiles of thousands of genes between species that develop in different rates. Here we review recent comparative studies that tackled the challenges of accurate staging and large-scale analyses using different computational approaches. These studies reveal how correct temporal scaling exposes the striking conservation of developmental gene expression between morphologically similar species. Different clustering approaches are used to address various comparative questions and identify the conservation and divergence of large gene sets. We discuss the unexpected contribution of housekeeping genes to the interspecies correlations and how this contribution distorts the hourglass pattern generated by developmental genes. Overall, we demonstrate how comparative studies of gene expression kinetics can provide novel insights into the developmental constraints and plasticity that shape animal body plans.

Interactions between genetic variation and cellular environment in skeletal muscle gene expression

From whole organisms to individual cells, responses to environmental conditions are influenced by genetic makeup, where the effect of genetic variation on a trait depends on the environmental context. RNA-sequencing quantifies gene expression as a molecular trait, and is capable of capturing both genetic and environmental effects. In this study, we explore opportunities of using allele-specific expression (ASE) to discover cis-acting genotype-environment interactions (GxE)—genetic effects on gene expression that depend on an environmental condition. Treating 17 common, clinical traits as approximations of the cellular environment of 267 skeletal muscle biopsies, we identify 10 candidate environmental response expression quantitative trait loci (reQTLs) across 6 traits (12 unique gene-environment trait pairs; 10% FDR per trait) including sex, systolic blood pressure, and low-density lipoprotein cholesterol. Although using ASE is in principle a promising approach to detect GxE effects, replication of such signals can be challenging as validation requires harmonization of environmental traits across cohorts and a sufficient sampling of heterozygotes for a transcribed SNP. Comprehensive discovery and replication will require large human transcriptome datasets, or the integration of multiple transcribed SNPs, coupled with standardized clinical phenotyping.

Parallel embryonic transcriptional programs evolve under distinct constraints and may enable morphological conservation amidst adaptation

Embryonic development evolves by balancing stringent morphological constraints with genetic and environmental variation. The design principle that allows developmental transcriptional programs to conserve embryonic morphology while adapting to environmental changes is still not fully understood. To address this fundamental challenge, we compare developmental transcriptomes of two sea urchin species, Paracentrotus lividus and Strongylocentrotus purpuratus, that shared a common ancestor about 40 million years ago and are geographically distant yet show similar morphology. We find that both developmental and housekeeping genes show highly dynamic and strongly conserved temporal expression patterns. The expression of other gene sets, including homeostasis and response genes, show divergent expression which could result from either evolutionary drift or adaptation to local environmental conditions. The interspecies correlations of developmental gene expressions are highest between morphologically similar developmental time points whereas the interspecies correlations of housekeeping gene expression are high between all the late zygotic time points. Relatedly, the position of the phylotypic stage varies between these two groups of genes: developmental gene expression shows highest conservation at mid-developmental stage, in agreement with the hourglass model while the conservation of housekeeping genes keeps increasing with developmental time. When all genes are combined, the relationship between conservation of gene expression and morphological similarity is partially masked by housekeeping genes and genes with diverged expression. Our study illustrates various transcriptional programs that coexist in the developing embryo and evolve under different constraints. Apparently, morphological constraints underlie the conservation of developmental gene expression while embryonic fitness requires the conservation of housekeeping gene expression and the species-specific adjustments of homeostasis gene expression. The distinct evolutionary forces acting on these transcriptional programs enable the conservation of similar body plans while allowing adaption.

New inter-correlated genes targeted by diatom-derived polyunsaturated aldehydes in the sea urchin Paracentrotus lividus

The marine environment is continually subjected to the action of stressors (including natural toxins), which represent a constant danger for benthic communities. In the present work using network analysis we identified ten genes on the basis of associated functions (FOXA, FoxG, GFI-1, nodal, JNK, OneCut/Hnf6, TAK1, tcf4, TCF7, VEGF) in the sea urchin Paracentrotus lividus, having key roles in different processes, such as embryonic development and asymmetry, cell fate specification, cell differentiation and morphogenesis, and skeletogenesis. These genes are correlated with three HUB genes, Foxo, Jun and HIF1A. Real Time qPCR revealed that during sea urchin embryonic development the expression levels of these genes were modulated by three diatom-derived polyunsaturated aldehydes (PUAs), decadienal, heptadienal and octadienal. Our findings show how changes in gene expression levels may be used as an early indicator of stressful conditions in the marine environment. The identification of key genes and the molecular pathways in which they are involved represents a fundamental tool in understanding how marine organisms try to afford protection against toxicants, to avoid deleterious consequences and irreversible damages. The genes identified in this work as targets for PUAs can be considered as possible biomarkers to detect exposure to different environmental pollutants.

Genomic Characterization of the Evolutionary Potential of the Sea Urchin Strongylocentrotus droebachiensis Facing Ocean Acidification

Ocean acidification (OA) is increasing due to anthropogenic CO₂ emissions and poses a threat to marine species and communities worldwide. To better project the effects of acidification on organisms’ health and persistence, an understanding is needed of the 1) mechanisms underlying developmental and physiological tolerance and 2) potential populations have for rapid evolutionary adaptation. This is especially challenging in nonmodel species where targeted assays of metabolism and stress physiology may not be available or economical for large-scale assessments of genetic constraints. We used mRNA sequencing and a quantitative genetics breeding design to study mechanisms underlying genetic variability and tolerance to decreased seawater pH (-0.4 pH units) in larvae of the sea urchin Strongylocentrotus droebachiensis. We used a gene ontology-based approach to integrate expression profiles into indirect measures of cellular and biochemical traits underlying variation in larval performance (i.e., growth rates). Molecular responses to OA were complex, involving changes to several functions such as growth rates, cell division, metabolism, and immune activities. Surprisingly, the magnitude of pH effects on molecular traits tended to be small relative to variation attributable to segregating functional genetic variation in this species. We discuss how the application of transcriptomics and quantitative genetics approaches across diverse species can enrich our understanding of the biological impacts of climate change.

Loggerhead sea turtle embryos (Caretta caretta) regulate expression of stress response and developmental genes when exposed to a biologically realistic heat stress

Oviparous reptile embryos are expected to breach their critical thermal maxima if temperatures reach those predicted under current climate change models due to the lack of the maternal buffering processes and parental care. Heat-shock proteins (HSPs) are integral in the molecular response to thermal stress, and their expression is heritable, but the roles of other candidate families such as the heat-shock factors (HSFs) have not been determined in reptiles. Here, we subject embryonic sea turtles (Caretta caretta) to a biologically realistic thermal stress and employ de novo transcriptomic profiling of brain tissue to investigate the underlying molecular response. From a reference transcriptome of 302 293 transcripts, 179 were identified as differentially expressed between treatments. As anticipated, genes enriched in the heat-shock treatment were primarily associated with the Hsp families, or were genes whose products play similar protein editing and chaperone functions (e.g. bag3, MYOC and serpinh1). Unexpectedly, genes encoding the HSFs were not significantly upregulated under thermal stress, indicating their presence in unstressed cells in an inactive state. Genes that were downregulated under thermal stress were less well functionally defined but were associated with stress response, development and cellular organization, suggesting that developmental processes may be compromised at realistically high temperatures. These results confirm that genes from the Hsp families play vital roles in the thermal tolerance of developing reptile embryos and, in addition with a number of other genes, should be targets for evaluating the capacity of oviparous reptiles to respond adaptively to the effects of climate change.

Gene expression profiling during the embryo-to-larva transition in the giant red sea urchin Mesocentrotus franciscanus

In echinoderms, major morphological transitions during early development are attributed to different genetic interactions and changes in global expression patterns that shape the regulatory program for the specification of embryonic territories. In order more thoroughly to understand these biological and molecular processes, we examined the transcriptome structure and expression profiles during the embryo‐to‐larva transition of a keystone species, the giant red sea urchin Mesocentrotus franciscanus. Using a de novo assembly approach, we obtained 176,885 transcripts from which 60,439 (34%) had significant alignments to known proteins. From these transcripts, ~80% were functionally annotated allowing the identification of ~2,600 functional, structural, and regulatory genes involved in developmental process. Analysis of expression profiles between gastrula and pluteus stages of M. franciscanus revealed 791 differentially expressed genes with 251 GO overrepresented terms. For gastrula, up‐regulated GO terms were mainly linked to cell differentiation and signal transduction involved in cell cycle checkpoints. In the pluteus stage, major GO terms were associated with phosphoprotein phosphatase activity, muscle contraction, and olfactory behavior, among others. Our evolutionary comparative analysis revealed that several of these genes and functional pathways are highly conserved among echinoids, holothuroids, and ophiuroids.

Modulation of yellow expression contributes to thermal plasticity of female abdominal pigmentation in Drosophila melanogaster

Phenotypic plasticity describes the ability of a given genotype to produce distinct phenotypes in different environments. We use the temperature sensitivity of abdominal pigmentation in Drosophila melanogaster females as a model to analyse the effect of the environment on development. We reported previously that thermal plasticity of abdominal pigmentation in females involves the pigmentation gene tan (t). However, the expression of the pigmentation gene yellow (y) was also modulated by temperature in the abdominal epidermis of pharate females. We investigate here the contribution of y to female abdominal pigmentation plasticity. First, we show that y is required for the production of black Dopamine-melanin. Then, using in situ hybridization, we show that the expression of y is strongly modulated by temperature in the abdominal epidermis of pharate females but not in bristles. Interestingly, these two expression patterns are known to be controlled by distinct enhancers. However, the activity of the y-wing-body epidermal enhancer only partially mediates the effect of temperature suggesting that additional regulatory sequences are involved. In addition, we show that y and t co-expression is needed to induce strong black pigmentation indicating that y contributes to female abdominal pigmentation plasticity.

Ocean acidification has little effect on developmental thermal windows of echinoderms from Antarctica to the tropics

As the ocean warms, thermal tolerance of developmental stages may be a key driver of changes in the geographical distributions and abundance of marine invertebrates. Additional stressors such as ocean acidification may influence developmental thermal windows and are therefore important considerations for predicting distributions of species under climate change scenarios. The effects of reduced seawater pH on the thermal windows of fertilization, embryology and larval morphology were examined using five echinoderm species: two polar (Sterechinus neumayeri and Odontaster validus), two temperate (Fellaster zelandiae and Patiriella regularis) and one tropical (Arachnoides placenta). Responses were examined across 12-13 temperatures ranging from -1.1 °C to 5.7 °C (S. neumayeri), -0.5 °C to 10.7 °C (O. validus), 5.8 °C to 27 °C (F. zelandiae), 6.0 °C to 27.1 °C (P. regularis) and 13.9 °C to 34.8 °C (A. placenta) under present-day and near-future (2100+) ocean acidification conditions (-0.3 pH units) and for three important early developmental stages 1) fertilization, 2) embryo (prehatching) and 3) larval development. Thermal windows for fertilization were broad and were not influenced by a pH decrease. Embryological development was less thermotolerant. For O. validus, P. regularis and A. placenta, low pH reduced normal development, albeit with no effect on thermal windows. Larval development in all five species was affected by both temperature and pH; however, thermal tolerance was not reduced by pH. Results of this study suggest that in terms of fertilization and development, temperature will remain as the most important factor influencing species' latitudinal distributions as the ocean continues to warm and decrease in pH, and that there is little evidence of a synergistic effect of temperature and ocean acidification on the thermal control of species ranges.

Population-dependent effects of ocean acidification

Elevated carbon dioxide levels and the resultant ocean acidification (OA) are changing the abiotic conditions of the oceans at a greater rate than ever before and placing pressure on marine species. Understanding the response of marine fauna to this change is critical for understanding the effects of OA. Population-level variation in OA tolerance is highly relevant and important in the determination of ecosystem resilience and persistence, but has received little focus to date. In this study, whether OA has the same biological consequences in high-salinity-acclimated population versus a low-salinity-acclimated population of the same species was investigated in the marine isopod Idotea balthica.The populations were found to have physiologically different responses to OA. While survival rate was similar between the two study populations at a future CO2 level of 1000 ppm, and both populations showed increased oxidative stress, the metabolic rate and osmoregulatory activity differed significantly between the two populations. The results of this study demonstrate that the physiological response to OA of populations from different salinities can vary. Population-level variation and the environment provenance of individuals used in OA experiments should be taken into account for the evaluation and prediction of climate change effects.

A new gene regulatory network model based on BP algorithm for interrogating differentially expressed genes of Sea Urchin

Background

Computer science and mathematical theories are combined to analyze the complex interactions among genes, which are simplified to a network to establish a theoretical model for the analysis of the structure, module and dynamic properties. In contrast, traditional model of gene regulatory networks often lack an effective method for solving gene expression data because of high durational and spatial complexity. In this paper, we propose a new model for constructing gene regulatory networks using back propagation (BP) neural network based on predictive function and network topology.

Results

Combined with complex nonlinear mapping and self-learning, the BP neural network was mapped into a complex network. Network characteristics were obtained from the parameters of the average path length, average clustering coefficient, average degree, modularity, and map’s density to simulate the real gene network by an artificial network. Through the statistical analysis and comparison of network parameters of Sea Urchin mRNA microarray data under different temperatures, the value of network parameters was observed. Differentially expressed Sea Urchin genes associated with temperature were determined by calculating the difference in the degree of each gene from different networks.

Conclusion

The new model we developed is suitable to simulate gene regulatory network and has capability of determining differentially expressed genes.

Subtle reproductive impairment through nitric oxide-mediated mechanisms in sea urchins from an area affected by harmful algal blooms

The health of the sea urchin Paracentrotus lividus, a key species in the Mediterranean Sea, is menaced by several pressures in coastal environments. Here, we aimed at assessing the reproductive ability of apparently healthy P. lividus population in a marine protected area affected by toxic blooms of Ostreospsis cf. ovata. Wide-ranging analyses were performed in animals collected prior to and during the bloom, as well as at several times thereafter, during the reproductive season. Adults showed a low fertilization rate, along with high nitric oxide (NO) levels in the gonads and the nitration of the major yolk protein toposome, which is an important player in sea urchin development. Serious developmental anomalies were observed in the progeny, which persist several months after the bloom. NO levels were high in the different developmental stages, which also showed variations in the transcription of several genes that were found to be directly or indirectly modulated by NO. These results highlight subtle but important reproductive flaws transmitted from the female gonads to the offspring with the NO involvement. Despite a recovery along time after the bloom, insidious damages can be envisaged in the local sea urchin population, with possible reverberation on the whole benthic system.

Robustness and Accuracy in Sea Urchin Developmental Gene Regulatory Networks

Developmental gene regulatory networks robustly control the timely activation of regulatory and differentiation genes. The structure of these networks underlies their capacity to buffer intrinsic and extrinsic noise and maintain embryonic morphology. Here I illustrate how the use of specific architectures by the sea urchin developmental regulatory networks enables the robust control of cell fate decisions. The Wnt-βcatenin signaling pathway patterns the primary embryonic axis while the BMP signaling pathway patterns the secondary embryonic axis in the sea urchin embryo and across bilateria. Interestingly, in the sea urchin in both cases, the signaling pathway that defines the axis controls directly the expression of a set of downstream regulatory genes. I propose that this direct activation of a set of regulatory genes enables a uniform regulatory response and a clear cut cell fate decision in the endoderm and in the dorsal ectoderm. The specification of the mesodermal pigment cell lineage is activated by Delta signaling that initiates a triple positive feedback loop that locks down the pigment specification state. I propose that the use of compound positive feedback circuitry provides the endodermal cells enough time to turn off mesodermal genes and ensures correct mesoderm vs. endoderm fate decision. Thus, I argue that understanding the control properties of repeatedly used regulatory architectures illuminates their role in embryogenesis and provides possible explanations to their resistance to evolutionary change.

Temperature-induced variation in gene expression burst size in metazoan cells

Background

Gene expression is an inherently stochastic process, owing to its dynamic molecular nature. Protein amount distributions, which can be acquired by cytometry using a reporter gene, can inform about the mechanisms of the underlying microscopic molecular system.

Results

By using different clones of chicken erythroid progenitor cells harboring different integration sites of a CMV-driven mCherry protein, we investigated the dynamical behavior of such distributions. We show that, on short term, clone distributions can be quickly regenerated from small population samples with a high accuracy. On longer term, on the contrary, we show variations manifested by correlated fluctuation in the Mean Fluorescence Intensity. In search for a possible cause of this correlation, we demonstrate that in response to small temperature variations cells are able to adjust their gene expression rate: a modest (2 °C) increase in external temperature induces a significant down regulation of mean expression values, with a reverse effect observed when the temperature is decreased. Using a two-state model of gene expression we further demonstrate that temperature acts by modifying the size of transcription bursts, while the burst frequency of the investigated promoter is less systematically affected.

Conclusions

For the first time, we report that transcription burst size is a key parameter for gene expression that metazoan cells from homeotherm animals can modify in response to an external thermal stimulus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12867-015-0048-2) contains supplementary material, which is available to authorized users.

Parental Effect of Long Acclimatization on Thermal Tolerance of Juvenile Sea Cucumber Apostichopus japonicus

To evaluate the thermal resistance of marine invertebrates to elevated temperatures under scenarios of future climate change, it is crucial to understand parental effect of long acclimatization on thermal tolerance of offspring. To test whether there is parental effect of long acclimatization, adult sea cucumbers (Apostichopus japonicus) from the same broodstock were transplanted southward and acclimatized at high temperature in field mesocosms. Four groups of juvenile sea cucumbers whose parents experienced different durations of high temperature acclimatization were established. Upper thermal limits, oxygen consumption and levels of heat shock protein mRNA of juveniles was determined to compare thermal tolerance of individuals from different groups. Juvenile sea cucumbers whose parents experienced high temperature could acquire high thermal resistance. With the increase of parental exposure duration to high temperature, offspring became less sensitive to high temperature, as indicated by higher upper thermal limits (LT₅₀), less seasonal variations of oxygen consumption, and stable oxygen consumption rates between chronic and acute thermal stress. The relatively high levels of constitutive expression of heat-shock proteins should contribute to the high thermal tolerance. Together, these results indicated that the existence of a parental effect of long acclimatization would increase thermal tolerance of juveniles and change the thermal sensitivity of sea cucumber to future climate change.

Temperature-Related Reaction Norms of Gene Expression: Regulatory Architecture and Functional Implications

The environment has profound effects on the expression of many traits and reaction norms describe the expression dynamics of a trait across a broad range of environmental conditions. Here, we analyze gene expression in Drosophila melanogaster across four different developmental temperatures (13-29 °C). Gene expression is highly plastic with 83.3% of the genes being differentially expressed. We distinguished three components of plasticity: 1) Dynamics of gene expression intensity (sum of change), 2) direction of change, and 3) curvature of the reaction norm (linear vs. quadratic). Studying their regulatory architecture we found that all three plasticity components were most strongly affected by the number of different transcription factors (TFs) binding to the target gene. More TFs were found in genes with less expression changes across temperatures. Although the effect of microRNAs was weaker, we consistently noted a trend in the opposite direction. The most plastic genes were regulated by fewer TFs and more microRNAs than less plastic genes. Different patterns of plasticity were also reflected by their functional characterization based on gene ontology. Our results suggest that reaction norms provide an important key to understand the functional requirements of natural populations exposed to variable environmental conditions.

Comparative Study of Regulatory Circuits in Two Sea Urchin Species Reveals Tight Control of Timing and High Conservation of Expression Dynamics

Accurate temporal control of gene expression is essential for normal development and must be robust to natural genetic and environmental variation. Studying gene expression variation within and between related species can delineate the level of expression variability that development can tolerate. Here we exploit the comprehensive model of sea urchin gene regulatory networks and generate high-density expression profiles of key regulatory genes of the Mediterranean sea urchin, Paracentrotus lividus (Pl). The high resolution of our studies reveals highly reproducible gene initiation times that have lower variation than those of maximal mRNA levels between different individuals of the same species. This observation supports a threshold behavior of gene activation that is less sensitive to input concentrations. We then compare Mediterranean sea urchin gene expression profiles to those of its Pacific Ocean relative, Strongylocentrotus purpuratus (Sp). These species shared a common ancestor about 40 million years ago and show highly similar embryonic morphologies. Our comparative analyses of five regulatory circuits operating in different embryonic territories reveal a high conservation of the temporal order of gene activation but also some cases of divergence. A linear ratio of 1.3-fold between gene initiation times in Pl and Sp is partially explained by scaling of the developmental rates with temperature. Scaling the developmental rates according to the estimated Sp-Pl ratio and normalizing the expression levels reveals a striking conservation of relative dynamics of gene expression between the species. Overall, our findings demonstrate the ability of biological developmental systems to tightly control the timing of gene activation and relative dynamics and overcome expression noise induced by genetic variation and growth conditions.

Maternal Exposure to Cadmium and Manganese Impairs Reproduction and Progeny Fitness in the Sea Urchin Paracentrotus lividus

Metal contamination represents one of the major sources of pollution in marine environments. In this study we investigated the short-term effects of ecologically relevant cadmium and manganese concentrations (10(-6) and 3.6 x 10(-5) M, respectively) on females of the sea urchin Paracentrotus lividus and their progeny, reared in the absence or presence of the metal. Cadmium is a well-known heavy metal, whereas manganese represents a potential emerging contaminant, resulting from an increased production of manganese-containing compounds. The effects of these agents were examined on both P. lividus adults and their offspring following reproductive state, morphology of embryos, nitric oxide (NO) production and differential gene expression. Here, we demonstrated that both metals differentially impaired the fertilization processes of the treated female sea urchins, causing modifications in the reproductive state and also affecting NO production in the ovaries. A detailed analysis of the progeny showed a high percentage of abnormal embryos, associated to an increase in the endogenous NO levels and variations in the transcriptional expression of several genes involved in stress response, skeletogenesis, detoxification, multi drug efflux processes and NO production. Moreover, we found significant differences in the progeny from females exposed to metals and reared in metal-containing sea water compared to embryos reared in non-contaminated sea water. Overall, these results greatly expanded previous studies on the toxic effects of metals on P. lividus and provided new insights into the molecular events induced in the progeny of sea urchins exposed to metals.

Ocean acidification impacts on sperm mitochondrial membrane potential bring sperm swimming behaviour near its tipping point

Broadcast spawning marine invertebrates are susceptible to environmental stressors such as climate change, as their reproduction depends on the successful meeting and fertilization of gametes in the water column. Under near-future scenarios of ocean acidification, the swimming behaviour of marine invertebrate sperm is altered. We tested whether this was due to changes in sperm mitochondrial activity by investigating the effects of ocean acidification on sperm metabolism and swimming behaviour in the sea urchin Centrostephanus rodgersii. We used a fluorescent molecular probe (JC-1) and flow cytometry to visualize mitochondrial activity (measured as change in mitochondrial membrane potential, MMP). Sperm MMP was significantly reduced in ΔpH −0.3 (35% reduction) and ΔpH −0.5 (48% reduction) treatments, whereas sperm swimming behaviour was less sensitive with only slight changes (up to 11% decrease) observed overall. There was significant inter-individual variability in responses of sperm swimming behaviour and MMP to acidified seawater. We suggest it is likely that sperm exposed to these changes in pH are close to their tipping point in terms of physiological tolerance to acidity. Importantly, substantial inter-individual variation in responses of sperm swimming to ocean acidification may increase the scope for selection of resilient phenotypes, which, if heritable, could provide a basis for adaptation to future ocean acidification.

The environmental genomics of metazoan thermal adaptation

Continued and accelerating change in the thermal environment places an ever-greater priority on understanding how organisms are going to respond. The paradigm of ‘move, adapt or die', regarding ways in which organisms can respond to environmental stressors, stimulates intense efforts to predict the future of biodiversity. Assuming that extinction is an unpalatable outcome, researchers have focussed attention on how organisms can shift in their distribution to stay in the same thermal conditions or can stay in the same place by adapting to a changing thermal environment. How likely these respective outcomes might be depends on the answer to a fundamental evolutionary question, namely what genetic changes underpin adaptation to the thermal environment. The increasing access to and decreasing costs of next-generation sequencing (NGS) technologies, which can be applied to both model and non-model systems, provide a much-needed tool for understanding thermal adaptation. Here we consider broadly what is already known from non-NGS studies about thermal adaptation, then discuss the benefits and challenges of different NGS methodologies to add to this knowledge base. We then review published NGS genomics and transcriptomics studies of thermal adaptation to heat stress in metazoans and compare these results with previous non-NGS patterns. We conclude by summarising emerging patterns of genetic response and discussing future directions using these increasingly common techniques.

Temperature stress mediates decanalization and dominance of gene expression in Drosophila melanogaster

The regulatory architecture of gene expression remains an area of active research. Here, we studied how the interplay of genetic and environmental variation affects gene expression by exposing Drosophila melanogaster strains to four different developmental temperatures. At 18°C we observed almost complete canalization with only very few allelic effects on gene expression. In contrast, at the two temperature extremes, 13°C and 29°C a large number of allelic differences in gene expression were detected due to both cis- and trans-regulatory effects. Allelic differences in gene expression were mainly dominant, but for up to 62% of the genes the dominance swapped between 13 and 29°C. Our results are consistent with stabilizing selection causing buffering of allelic expression variation in non-stressful environments. We propose that decanalization of gene expression in stressful environments is not only central to adaptation, but may also contribute to genetic disorders in human populations.

Transgenerational acclimation of fishes to climate change and ocean acidification

There is growing concern about the impacts of climate change and ocean acidification on marine organisms and ecosystems, yet the potential for acclimation and adaptation to these threats is poorly understood. Whereas many short-term experiments report negative biological effects of ocean warming and acidification, new studies show that some marine species have the capacity to acclimate to warmer and more acidic environments across generations. Consequently, transgenerational plasticity may be a powerful mechanism by which populations of some species will be able to adjust to projected climate change. Here, I review recent advances in understanding transgenerational acclimation in fishes. Research over the past 2 to 3 years shows that transgenerational acclimation can partially or fully ameliorate negative effects of warming, acidification, and hypoxia in a range of different species. The molecular and cellular pathways underpinning transgenerational acclimation are currently unknown, but modern genetic methods provide the tools to explore these mechanisms. Despite the potential benefits of transgenerational acclimation, there could be limitations to the phenotypic traits that respond transgenerationally, and trade-offs between life stages, that need to be investigated. Future studies should also test the potential interactions between transgenerational plasticity and genetic evolution to determine how these two processes will shape adaptive responses to environmental change over coming decades.