Spatial analysis of gene regulation reveals new insights into the molecular basis of upper thermal limits

Expression of Truncated Products at the 5'-Terminal Region of RIPK2 and Evolutive Aspects that Support Their Biological Importance

Alternative splicing is the process of generating different mRNAs from the same primary transcript, which contributes to increase the transcriptome and proteome diversity. Abnormal splicing has been associated with the development of several diseases including cancer. Given that mutations and abnormal levels of the RIPK2 transcript and RIP-2 protein are frequent in tumors, and that RIP-2 modulates immune and inflammatory responses, we investigated alternative splicing events that result in partial deletions of the kinase domain at the N-terminus of RIP-2. We also investigated the structure and expression of the RIPK2 truncated variants and isoforms in different environments. In addition, we searched data throughout Supraprimates evolution that could support the biological importance of RIPK2 alternatively spliced products. We observed that human variants and isoforms were differentially regulated following temperature stress, and that the truncated transcript was more expressed than the long transcript in tumor samples. The inverse was found for the longer protein isoform. The truncated variant was also detected in chimpanzee, gorilla, hare, pika, mouse, rat, and tree shrew. The fact that the same variant has been preserved in mammals with divergence times up to 70 million years raises the hypothesis that it may have a functional significance.

Transcriptomic Analysis Pipeline (TAP) for quality control and functional assessment of transcriptomes

Background

RNA-sequencing (RNA-seq) has revolutionized the exploration of biological mechanisms, shedding light on the roles of non-coding RNAs, including long non-coding RNAs (lncRNAs), across various biological processes, including stress responses. Despite these advancements, there remains a gap in our understanding of the implications of different RNA-seq library protocols on comprehensive lncRNA expression analysis, particularly in non-mammalian organisms.

Results

In this study, we sought to bridge this knowledge gap by investigating lncRNA expression patterns in Drosophila melanogaster under thermal stress conditions. To achieve this, we conducted a comparative analysis of two RNA-seq library protocols: polyA + RNA capture and rRNA-depletion. Our approach involved the development and application of a Transcriptome Analysis Pipeline (TAP) designed to systematically assess both the technical and functional dimensions of RNA-seq, facilitating a robust comparison of these library protocols. Our findings underscore the efficacy of the polyA + protocol in capturing the majority of expressed lncRNAs within the Drosophila melanogaster transcriptome. In contrast, rRNA-depletion exhibited limited advantages in the context of D. melanogaster studies. Notably, the polyA + protocol demonstrated superior performance in terms of usable read yield and the accurate detection of splice junctions.

Conclusions

Our study introduces a versatile transcriptomic analysis pipeline, TAP, designed to uniformly process RNA-seq data from any organism with a reference genome. It also highlights the significance of selecting an appropriate RNA-seq library protocol tailored to the specific research context.

Background

Advances in next generation sequencing (NGS) technologies enable the comprehensive analysis of genetic sequences of organisms in a relatively cost-effective manner [1, 2]. Among these technologies, RNA-sequencing (RNA-seq) has emerged as a preeminent method to study fundamental biological mechanisms at the level of cells, tissues, and whole organisms. RNA-seq enables the detection and quantification of various RNA populations, including messenger RNA (mRNA) and various species of non-coding RNA, such as long non-coding RNA (lncRNA), as well as an assessment of features including splice junctions in RNA.

Assessing the role of family level variation and heat shock gene expression in the thermal stress response of the mosquito Aedes aegypti

The geographical range of the mosquito vector for many human disease-causing viruses, Aedes aegypti, is expanding, in part owing to changing climate. The capacity of this species to adapt to thermal stress will affect its future distributions. It is unclear how much heritable genetic variation may affect the upper thermal limits of mosquito populations over the long term. Nor are the genetic pathways that confer thermal tolerance fully understood. In the short term, cells induce a plastic, protective response known as 'heat shock'. Using a physiological 'knockdown' assay, we investigated mosquito thermal tolerance to characterize the genetic architecture of the trait. While families representing the extreme ends of the distribution for knockdown time differed from one another, the trait exhibited low but non-zero broad-sense heritability. We then explored whether families representing thermal performance extremes differed in their heat shock response by measuring gene expression of heat shock protein-encoding genes Hsp26, Hsp83 and Hsp70. Contrary to prediction, the families with higher thermal tolerance demonstrated less Hsp expression. This pattern may indicate that other mechanisms of heat tolerance, rather than heat shock, may underpin the stress response, and the costly production of HSPs may instead signal poor adaptation. This article is part of the theme issue 'Infectious disease ecology and evolution in a changing world'.

Temporal specific coevolution of Hsp70 and co-chaperone stv expression in Drosophila melanogaster under selection for heat tolerance

Heat shock proteins (Hsps) have long been candidates for ecological adaptation given their unequivocal role in mitigating cell damage from heat stress, but linking Hsps to heat tolerance has proven difficult given the complexity of thermal adaptation. Experimental evolution has been utilized to examine direct and correlated responses to selection for increased heat tolerance in Drosophila, often focusing on the major Hsp family Hsp70 and/or the master regulator HSF as a selection response, but rarely on other aspects of the heat shock complex. We examined Hsp70 and co-chaperone stv isoform transcript expression in Australian D. melanogaster lines selected for static heat tolerance, and observed a temporal and stv isoform specific, coordinated transcriptional selection response with Hsp70, suggesting that increased chaperone output accompanied increased heat tolerance. We hypothesize that the coordinated evolutionary response of Hsp70 and stv may have arisen as a correlated response resulting from a shared regulatory hierarchy. Our work highlights the complexity and specificity of the heat shock response in D. melanogaster. The selected lines examined also showed correlated responses for other measures of heat tolerance, and the coevolution of Hsp70 and stv provide new avenues to examine the common mechanisms underpinning direct and correlated phenotypic responses to selection for heat tolerance.

Interactions between developmental and adult acclimation have distinct consequences for heat tolerance and heat stress recovery

Developmental and adult thermal acclimation can have distinct, even opposite, effects on adult heat resistance in ectotherms. Yet, their relative contribution to heat-hardiness of ectotherms remains unclear despite the broad ecological implications thereof. Furthermore, the deterministic relationship between heat knockdown and recovery from heat stress is poorly understood but significant for establishing causal links between climate variability and population dynamics. Here, using Drosophila melanogaster in a full-factorial experimental design, we assessed the heat tolerance of flies in static stress assays, and document how developmental and adult acclimation interact with a distinct pattern to promote survival to heat stress in adults. We show that warmer adult acclimation is the initial factor enhancing survival to constant stressful high temperatures in flies, but also that the interaction between adult and developmental acclimation becomes gradually more important to ensure survival as the stress persists. This provides an important framework revealing the dynamic interplay between these two forms of acclimation that ultimately enhance thermal tolerance as a function of stress duration. Furthermore, by investigating recovery rates post-stress, we also show that the process of heat-hardening and recovery post-heat knockdown are likely to be based on set of (at least partially) divergent mechanisms. This could bear ecological significance as a trade-off may exist between increasing thermal tolerance and maximizing recovery rates post-stress, constraining population responses when exposed to variable and stressful climatic conditions.

Evolution of Plasticity in Response to Ethanol between Sister Species with Different Ecological Histories (Drosophila melanogaster and D. simulans)

AbstractWhen populations evolve adaptive reaction norms in response to novel environments, it can occur through a process termed genetic accommodation. Under this model, the initial response to the environment is widely variable between genotypes as a result of cryptic genetic variation, which is then refined by selection to a single adaptive response. Here, I empirically test these predictions from genetic accommodation by measuring reaction norms in individual genotypes and across several time points. I compare two species of Drosophila that differ in their adaptation to ethanol (D. melanogaster and D. simulans). Both species are human commensals with a recent cosmopolitan expansion, but only D. melanogaster is adapted to ethanol exposure. Using gene expression as a phenotype and an approach that combines information about expression and alternative splicing, I find that D. simulans exhibits cryptic genetic variation in the response to ethanol, while D. melanogaster has almost no genotype-specific variation in reaction norm. This is evidence for adaptation to ethanol through genetic accommodation, suggesting that the evolution of phenotypic plasticity could be an important contributor to the ability to exploit novel resources.

Transcriptomic signatures of cold adaptation and heat stress in the winter ant (Prenolepis imparis)

Climate change is a serious threat to biodiversity; it is therefore important to understand how animals will react to this stress. Ectotherms, such as ants, are especially sensitive to the climate as the environmental temperature influences myriad aspects of their biology, from optimal foraging time to developmental rate. In this study, we conducted an RNA-seq analysis to identify stress-induced genes in the winter ant (Prenolepis imparis). We quantified gene expression during heat and cold stress relative to a control temperature. From each of our conditions, we sequenced the transcriptome of three individuals. Our de novo assembly included 13,324 contigs that were annotated against the nr and SwissProt databases. We performed gene ontology and enrichment analyses to gain insight into the physiological processes involved in the stress response. We identified a total of 643 differentially expressed genes across both treatments. Of these, only seven genes were differentially expressed in the cold-stressed ants, which could indicate that the temperature we chose for trials did not induce a strong stress response, perhaps due to the cold adaptations of this species. Conversely, we found a strong response to heat: 426 upregulated genes and 210 downregulated genes. Of these, ten were expressed at a greater than ten-fold change relative to the control. The transcripts we could identify included those encoding for protein folding genes, heat shock proteins, histones, and Ca2+ ion transport. One of these transcripts, hsc70-4L was found to be under positive selection. We also characterized the functional categories of differentially expressed genes. These candidate genes may be functionally conserved and relevant for related species that will deal with rapid climate change.

Validating measurements of acclimation for climate change adaptation

Acclimation and other forms of plasticity that can increase stress resistance feature strongly in discussions surrounding climate change impacts or vulnerability projections of insects and other ectotherms. There is interest in compiling databases for assessing the adequacy of acclimation for dealing with climate change. Here, we argue that the nature of acclimation is context dependent and therefore that estimates summarised across studies, especially those that have assayed stress using diverse methods, are limited in their utility when applied as a standardized metric or to a single general context such as average climate warming. Moreover, the dynamic nature of tolerances and acclimation drives important variation that is quickly obscured through many summary statistics or even in effect size analyses; retaining a strong focus on the temporal-level, population-level and treatment-level variance in forecasting climate change impacts on insects is essential. We summarise recent developments within the context of climate change and propose how future studies might validate the role of acclimation by integration across field studies and mechanistic modelling. Despite arguments to the contrary, to date no studies have convincingly demonstrated an important role for acclimation in recent climate change adaptation of insects. Paramount to these discussions is i) developing a strong conceptual framework for acclimation in the focal trait(s), ii) obtaining novel empirical data dissecting the fitness benefits and consequences of acclimation across diverse contexts and timescales, with iii) better coverage of under-represented geographic regions and taxa.

Highly dynamic transcriptional reprogramming and shorter isoform shifts under acute stresses during biological invasions

Phenotypic plasticity has been increasingly recognized for its importance in adaptation to novel environments, and initial rapid plastic response to acute stresses usually serves as the stepping stone for future adaptation. Differential gene expression and alternative splicing have been proposed as two underlying mechanisms for rapid plastic response to environmental stresses. Here, we used an invasive model species, Ciona savignyi, to investigate the temporary plastic changes under temperature stresses on gene expression and alternative splicing. Our results revealed rapid and highly dynamic gene expression reprogramming and alternative splicing switch under acute stresses. Distinct transcriptional response profiles were triggered by two types of temperature stresses, showing resilience recovery and increasing divergence under heat and cold challenges, respectively. Interestingly, alternative exons were more inclined to be skipped under both heat and cold stresses, leading to shorter isoforms but with maintained Open Reading Frames (ORFs). Although similar response patterns were observed between differential gene expression and alternative splicing, low overlap between Differentially Expressed Genes (DEGs) and Differentially Alternative Spliced Genes (DASGs) suggests that distinct gene sets and associated functions should be involved in temperature challenges. Thus, alternative splicing should offer an additional layer of plastic response to environmental challenges. Finally, we identified key plastic genes involved in both gene expression regulation and alternative splicing. The results obtained here shed light on adaptation and accommodation mechanisms during biological invasions, particularly for acute environmental changes at early stages of biological invasions such as transport and introduction.

Adaptation to the abiotic environment in insects: the influence of variability on ecophysiology and evolutionary genomics

Advances in tools to gather environmental, phenotypic, and molecular data have accelerated our ability to detect abiotic drivers of variation across the genome-to-phenome spectrum in model and non-model insects. However, differences in the spatial and temporal resolution of these data sets may create gaps in our understanding of linkages between environment, genotype, and phenotype that yield missed or misleading results about adaptive variation. In this review we highlight sources of variability that might impact studies of phenotypic and 'omic environmental adaptation, challenges to collecting data at relevant scales, and possible solutions that link intensive fine-scale reductionist studies of mechanisms to large-scale biogeographic patterns.

How important is thermal history? Evidence for lasting effects of developmental temperature on upper thermal limits in Drosophila melanogaster

A common practice in thermal biology is to take individuals directly from the field and estimate a range of thermal traits. These estimates are then used in studies aiming to understand broad scale distributional patterns, understanding and predicting the evolution of phenotypic plasticity, and generating predictions for climate change risk. However, the use of field-caught individuals in such studies ignores the fact that many traits are phenotypically plastic and will be influenced by the thermal history of the focal individuals. The current study aims to determine the extent to which estimates of upper thermal limits (CTmax), a frequently used measure for climate change risk, are sensitive to developmental and adult acclimation temperatures and whether these two forms of plasticity are reversible. Examining a temperate and tropical population of Drosophila melanogaster we show that developmental acclimation has a larger and more lasting effect on CTmax than adult acclimation. We also find evidence for an interaction between developmental and adult acclimation, particularly when flies are acclimated for a longer period, and that these effects can be population specific. These results suggest that thermal history can have lasting effects on estimates of CTmax. In addition, we provide evidence that developmental and/or adult acclimation are unlikely to contribute to substantial shifts in CTmax and that acclimation capacity may be constrained at higher temperatures.

Transcriptional profiles of plasticity for desiccation stress in Drosophila

We examined the transcriptional responses of desiccation resistance candidate genes in populations of Drosophila melanogaster divergent for desiccation resistance and in capacity to improve resistance via phenotypic plasticity. Adult females from temperate and tropical eastern Australian populations were exposed to a rapid desiccation hardening (RDH) treatment, and groups without RDH to acute desiccation stress, and the transcript expression of 12 candidate desiccation genes were temporally profiled during, and in recovery from stress. We found that desiccation exposure resulted in largely transitory, stress-specific transcriptional changes in all but one gene. However linking the expression profiles to the population-level phenotypic divergence was difficult given subtle, and time-point specific population expression variation. Nonetheless, rapid desiccation hardening had the largest effect on gene expression, resulting in distinct molecular profiles. We report a hitherto uncharacterised desiccation molecular hardening response where prior exposure essentially 'primes' genes to respond to subsequent stress without discernible transcript changes prior to stress. This, taken together with some population gene expression variation of several bona fide desiccation candidates associated with different water balance strategies speaks of the complexity of natural desiccation resistance and plasticity and provides new avenues for understanding the molecular basis of a trait of ecological significance.

RNA sequencing reveals differential thermal regulation mechanisms between sexes of Glanville fritillary butterfly in the Tianshan Mountains, China

The Glanville fritillary butterfly (Melitaea cinxia; Nymphalidae) has been extensively studied as a model species in metapopulation ecology. We investigated in the earlier studies that female butterflies exhibit higher thermal tolerance than males in the Tianshan Mountains of China. We aim to understand the molecular mechanism of differences of thermal responses between sexes. We used RNA-seq approach and performed de novo assembly of transcriptome to compare the gene expression patterns between two sexes after heat stress. All the reads were assembled into 84,376 transcripts and 72,701 unigenes. The number of differential expressed genes (DEGs) between control and heat shock samples was 175 and 268 for males and females, respectively. Heat shock proteins genes (hsps) were up-regulated in response to heat stress in both males and females. Most of the up-regulated hsps showed higher fold changes in males than in females. Females expressed more ribosomal subunit protein genes, transcriptional elongation factor genes, and methionine-rich storage protein genes, participating in protein synthesis. It indicated that protein synthesis is needed for females to replace the damaged proteins due to heat shock. In addition, aspartate decarboxylase might contribute to thermal tolerance in females. These differences in gene expression may at least partly explain the response to high temperature stress, and the fact that females exhibit higher thermal tolerance.

Thermal plasticity in Drosophila melanogaster populations from eastern Australia: quantitative traits to transcripts

The flexibility afforded to genotypes in different environments by phenotypic plasticity is of interest to biologists studying thermal adaptation because of the thermal lability of many traits. Differences in thermal performance and reaction norms can provide insight into the evolution of thermal adaptation to explore broader questions such as species distributions and persistence under climate change. One approach is to study the effects of temperature on fitness, morphological and more recently gene expression traits in populations from different climatic origins. The diverse climatic conditions experienced by Drosophila melanogaster along the eastern Australian temperate-tropical gradient are ideal given the high degree of continuous trait differentiation, but reaction norm variation has not been well studied in this system. Here, we reared a tropical and temperate population from the ends of the gradient over six developmental temperatures and examined reaction norm variation for five quantitative traits including thermal performance for fecundity, and reaction norms for thermotolerance, body size, viability and 23 transcript-level traits. Despite genetic variation for some quantitative traits, we found no differentiation between the populations for fecundity thermal optima and breadth, and the reaction norms for the other traits were largely parallel, supporting previous work suggesting that thermal evolution occurs by changes in trait means rather than by reaction norm shifts. We examined reaction norm variation in our expanded thermal regime for a gene set shown to previously exhibit GxE for expression plasticity in east Australian flies, as well as key heat-shock genes. Although there were differences in curvature between the populations suggesting a higher degree of thermal plasticity in expression patterns than for the quantitative traits, we found little evidence to support a role for genetic variation in maintaining expression plasticity.

The Effect of In Vitro Cultivation on the Transcriptome of Adult Brugia malayi

BACKGROUND

Filarial nematodes cause serious and debilitating infections in human populations of tropical countries, contributing to an entrenched cycle of poverty. Only one human filarial parasite, Brugia malayi, can be maintained in rodents in the laboratory setting. It has been a widely used model organism in experiments that employ culture systems, the impact of which on the worms is unknown.

METHODOLOGY/PRINCIPAL FINDINGS

Using Illumina RNA sequencing, we characterized changes in gene expression upon in vitro maintenance of adult B. malayi female worms at four time points: immediately upon removal from the host, immediately after receipt following shipment, and after 48 h and 5 days in liquid culture media. The dramatic environmental change and the 24 h time lapse between removal from the host and establishment in culture caused a globally dysregulated gene expression profile. We found a maximum of 562 differentially expressed genes based on pairwise comparison between time points. After an initial shock upon removal from the host and shipping, a few stress fingerprints remained after 48 h in culture and until the experiment was stopped. This was best illustrated by a strong and persistent up-regulation of several genes encoding cuticle collagens, as well as serpins.

CONCLUSIONS/SIGNIFICANCE

These findings suggest that B. malayi can be maintained in culture as a valid system for pharmacological and biological studies, at least for several days after removal from the host and adaptation to the new environment. However, genes encoding several stress indicators remained dysregulated until the experiment was stopped.

Constraints, independence, and evolution of thermal plasticity: probing genetic architecture of long- and short-term thermal acclimation

Seasonal and daily thermal variation can limit species distributions because of physiological tolerances. Low temperatures are particularly challenging for ectotherms, which use both basal thermotolerance and acclimation, an adaptive plastic response, to mitigate thermal stress. Both basal thermotolerance and acclimation are thought to be important for local adaptation and persistence in the face of climate change. However, the evolutionary independence of basal and plastic tolerances remains unclear. Acclimation can occur over longer (seasonal) or shorter (hours to days) time scales, and the degree of mechanistic overlap is unresolved. Using a midlatitude population of Drosophila melanogaster, we show substantial heritable variation in both short- and long-term acclimation. Rapid cold hardening (short-term plasticity) and developmental acclimation (long-term plasticity) are positively correlated, suggesting shared mechanisms. However, there are independent components of these traits, because developmentally acclimated flies respond positively to short-term acclimation. A strong negative correlation between basal cold tolerance and developmental acclimation suggests that basal cold tolerance may constrain developmental acclimation, whereas a weaker negative correlation between basal cold tolerance and short-term acclimation suggests less constraint. Using genome-wide association mapping, we show the genetic architecture of rapid cold hardening and developmental acclimation responses are nonoverlapping at the SNP and corresponding gene level. However, genes associated with each trait share functional similarities, including genes involved in apoptosis and autophagy, cytoskeletal and membrane structural components, and ion binding and transport. These results indicate substantial opportunity for short-term and long-term acclimation responses to evolve separately from each other and for short-term acclimation to evolve separately from basal thermotolerance.