Transcriptome Analysis Reveals Candidate Genes for Cold Tolerance in Drosophila ananassae

From whole bodies to single cells: A guide to transcriptomic approaches for ecology and evolutionary biology

RNA sequencing (RNAseq) methodology has experienced a burst of technological developments in the last decade, which has opened up opportunities for studying the mechanisms of adaptation to environmental factors at both the organismal and cellular level. Selecting the most suitable experimental approach for specific research questions and model systems can, however, be a challenge and researchers in ecology and evolution are commonly faced with the choice of whether to study gene expression variation in whole bodies, specific tissues, and/or single cells. A wide range of sometimes polarised opinions exists over which approach is best. Here, we highlight the advantages and disadvantages of each of these approaches to provide a guide to help researchers make informed decisions and maximise the power of their study. Using illustrative examples of various ecological and evolutionary research questions, we guide the readers through the different RNAseq approaches and help them identify the most suitable design for their own projects.

Temporal regulation of temperature tolerances and gene expression in an arctic insect

Terrestrial arthropods in the Arctic are exposed to highly variable temperatures that frequently reach cold and warm extremes. Yet, ecophysiological studies on arctic insects typically focus on the ability of species to tolerate low temperatures, whereas studies investigating physiological adaptations of species to periodically warm and variable temperatures are few. In this study, we investigated temporal changes in thermal tolerances and the transcriptome in the Greenlandic seed bug Nysius groenlandicus, collected in the field across different times and temperatures in Southern Greenland. We found that plastic changes in heat and cold tolerances occurred rapidly (within hours) and at a daily scale in the field, and that these changes are correlated with diurnal temperature variation. Using RNA sequencing, we provide molecular underpinnings of the rapid adjustments in thermal tolerance across ambient field temperatures and in the laboratory. We show that transcriptional responses are sensitive to daily temperature changes, and days characterized by high temperature variation induced markedly different expression patterns than thermally stable days. Further, genes associated with laboratory-induced heat responses, including expression of heat shock proteins and vitellogenins, were shared across laboratory and field experiments, but induced at time points associated with lower temperatures in the field. Cold stress responses were not manifested at the transcriptomic level.

Tropical super flies: Integrating Cas9 into Drosophila ananassae and its phenotypic effects

Ectotherms such as insects are animals whose body temperature largely depends on ambient temperature and temperature variations provide a selection pressure affecting the geographical distribution of these species. However, over the course of evolution, some insect species managed to colonize environments characterized by various temperature ranges. Therefore, insects provide an excellent study system to investigate the basis of adaptation to temperature changes and extremes. We are generally using the vinegar fly Drosophila ananassae as a model system to investigate the genetic basis of cold tolerance. This species has expanded from its tropical ancestral range to more temperate regions resulting in a cosmopolitan, domestic distribution. Previously, we identified candidate genes significantly associated with cold tolerance in this species. We now established molecular genetic tools to assess the function of these genes. Using CRISPR/Cas9 methodology for genome editing and the PiggyBac system, the Cas9 enzyme was successfully integrated into the genome of three fly strains with different levels of cold tolerance. We further report on preliminary findings that the Cas9 integration itself did not have a consistent effect on tolerance to cold. In conclusion, we offer with our study the molecular tools that allow studying stress-related candidate genes in D. ananassae in the future. In addition, we point out and provide guidance on the challenges that come with genome editing in a non-model species.

Transcriptomic and functional genetic evidence for distinct ecophysiological responses across complex life cycle stages

Organisms with complex life cycles demonstrate a remarkable ability to change their phenotypes across development, presumably as an evolutionary adaptation to developmentally variable environments. Developmental variation in environmentally sensitive performance, and thermal sensitivity in particular, has been well documented in holometabolous insects. For example, thermal performance in adults and juvenile stages exhibit little genetic correlation (genetic decoupling) and can evolve independently, resulting in divergent thermal responses. Yet, we understand very little about how this genetic decoupling occurs. We tested the hypothesis that genetic decoupling of thermal physiology is driven by fundamental differences in physiology between life stages, despite a potentially conserved Cellular Stress Response. We used RNAseq to compare transcript expression in response to a cold stressor in Drosophila melanogaster larvae and adults and used RNAi (RNA interference) to test whether knocking down nine target genes differentially affected larval and adult cold tolerance. Transcriptomic responses of whole larvae and adults during and following exposure to -5°C were largely unique both in identity of responding transcripts and in temporal dynamics. Further, we analyzed the tissue-specificity of differentially-expressed transcripts from FlyAtlas 2 data, and concluded that stage-specific differences in transcription were not simply driven by differences in tissue composition. In addition, RNAi of target genes resulted in largely stage-specific and sometimes sex-specific effects on cold tolerance. The combined evidence suggests that thermal physiology is largely stage-specific at the level of gene expression, and thus natural selection may be acting on different loci during the independent thermal adaptation of different life stages.

Identification and Characterization of UDP-Glycosyltransferase Genes in a Cerambycid Beetle, Pharsalia antennata Gahan, 1894 (Coleoptera: Cerambycidae)

The cerambycid beetle, Pharsalia antennata Gahan, 1894 (Coleoptera: Cerambycidae), is a wood-boring pest that spends most of its life cycle in the trunks or under the bark of trees. These distinctive biological characteristics make it likely that this beetle will encounter a number of plant defensive compounds, coupled with a broad range of host plants, possibly resulting in the overexpression or expansion of uridine diphosphate (UDP)-glycosyltransferase (UGT) genes. Here, we identified and characterized the UGT gene family in P. antennata through transcriptome data, sequence and phylogenetic analyses, and PCR and homology modeling approaches. In total, 59 transcripts encoding UGTs were identified, 34 of which harbored full-length sequences and shared high conservation with the UGTs of Anoplophora glabripennis. Of the 34 PantUGTs, only 31.78% amino acid identity was observed on average, but catalytic and sugar binding residues were highly conserved. Phylogenetic analyses revealed four Cerambycidae-specific clades, including 30 members from P. antennata. Combining the transcriptome and PCR data showed that PantUGTs had a wide tissue expression, and the majority of the genes were presented mainly in antennae or abdomens, suggesting their putative roles in olfaction and detoxification. This study provides, for the first time, information on the molecular and genetic basis of P. antennata, greatly enhancing our knowledge of the detoxification-related UGT gene family.

Phylogenomic analyses of the genus Drosophila reveals genomic signals of climate adaptation

Many Drosophila species differ widely in their distributions and climate niches, making them excellent subjects for evolutionary genomic studies. Here, we have developed a database of high‐quality assemblies for 46 Drosophila species and one closely related Zaprionus. Fifteen of the genomes were newly sequenced, and 20 were improved with additional sequencing. New or improved annotations were generated for all 47 species, assisted by new transcriptomes for 19. Phylogenomic analyses of these data resolved several previously ambiguous relationships, especially in the melanogaster species group. However, it also revealed significant phylogenetic incongruence among genes, mainly in the form of incomplete lineage sorting in the subgenus Sophophora but also including asymmetric introgression in the subgenus Drosophila. Using the phylogeny as a framework and taking into account these incongruences, we then screened the data for genome‐wide signals of adaptation to different climatic niches. First, phylostratigraphy revealed relatively high rates of recent novel gene gain in three temperate pseudoobscura and five desert‐adapted cactophilic mulleri subgroup species. Second, we found differing ratios of nonsynonymous to synonymous substitutions in several hundred orthologues between climate generalists and specialists, with trends for significantly higher ratios for those in tropical and lower ratios for those in temperate‐continental specialists respectively than those in the climate generalists. Finally, resequencing natural populations of 13 species revealed tropics‐restricted species generally had smaller population sizes, lower genome diversity and more deleterious mutations than the more widespread species. We conclude that adaptation to different climates in the genus Drosophila has been associated with large‐scale and multifaceted genomic changes.

A cis-regulatory element promoting increased transcription at low temperature in cultured ectothermic Drosophila cells

Background

Temperature change affects the myriad of concurrent cellular processes in a non-uniform, disruptive manner. While endothermic organisms minimize the challenge of ambient temperature variation by keeping the core body temperature constant, cells of many ectothermic species maintain homeostatic function within a considerable temperature range. The cellular mechanisms enabling temperature acclimation in ectotherms are still poorly understood. At the transcriptional level, the heat shock response has been analyzed extensively. The opposite, the response to sub-optimal temperature, has received lesser attention in particular in animal species. The tissue specificity of transcriptional responses to cool temperature has not been addressed and it is not clear whether a prominent general response occurs. Cis-regulatory elements (CREs), which mediate increased transcription at cool temperature, and responsible transcription factors are largely unknown.

Results

The ectotherm Drosophila melanogaster with a presumed temperature optimum around 25 °C was used for transcriptomic analyses of effects of temperatures at the lower end of the readily tolerated range (14–29 °C). Comparative analyses with adult flies and cell culture lines indicated a striking degree of cell-type specificity in the transcriptional response to cool. To identify potential cis-regulatory elements (CREs) for transcriptional upregulation at cool temperature, we analyzed temperature effects on DNA accessibility in chromatin of S2R+ cells. Candidate cis-regulatory elements (CREs) were evaluated with a novel reporter assay for accurate assessment of their temperature-dependency. Robust transcriptional upregulation at low temperature could be demonstrated for a fragment from the pastrel gene, which expresses more transcript and protein at reduced temperatures. This CRE is controlled by the JAK/STAT signaling pathway and antagonizing activities of the transcription factors Pointed and Ets97D.

Conclusion

Beyond a rich data resource for future analyses of transcriptional control within the readily tolerated range of an ectothermic animal, a novel reporter assay permitting quantitative characterization of CRE temperature dependence was developed. Our identification and functional dissection of the pst_E1 enhancer demonstrate the utility of resources and assay. The functional characterization of this CoolUp enhancer provides initial mechanistic insights into transcriptional upregulation induced by a shift to temperatures at the lower end of the readily tolerated range.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08057-4.

Natural variation in the transcriptional response of Drosophila melanogaster to oxidative stress

Broadly distributed species must cope with diverse and changing environmental conditions, including various forms of stress. Cosmopolitan populations of Drosophila melanogaster are more tolerant to oxidative stress than those from the species’ ancestral range in sub-Saharan Africa, and the degree of tolerance is associated with an insertion/deletion polymorphism in the 3′ untranslated region of the Metallothionein A (MtnA) gene that varies clinally in frequency. We examined oxidative stress tolerance and the transcriptional response to oxidative stress in cosmopolitan and sub-Saharan African populations of D. melanogaster, including paired samples with allelic differences at the MtnA locus. We found that the effect of the MtnA polymorphism on oxidative stress tolerance was dependent on the genomic background, with the deletion allele increasing tolerance only in a northern, temperate population. Genes that were differentially expressed under oxidative stress included MtnA and other metallothioneins, as well as those involved in glutathione metabolism and other genes known to be part of the oxidative stress response or the general stress response. A gene coexpression analysis revealed further genes and pathways that respond to oxidative stress including those involved in additional metabolic processes, autophagy, and apoptosis. There was a significant overlap among the genes induced by oxidative and cold stress, which suggests a shared response pathway to these two stresses. Interestingly, the MtnA deletion was associated with consistent changes in the expression of many genes across all genomic backgrounds, regardless of the expression level of the MtnA gene itself. We hypothesize that this is an indirect effect driven by the loss of microRNA binding sites within the MtnA 3′ untranslated region.

Thermal tolerance and preference are both consistent with the clinal distribution of house fly proto-Y chromosomes

Selection pressures can vary within localized areas and across massive geographical scales. Temperature is one of the best studied ecologically variable abiotic factors that can affect selection pressures across multiple spatial scales. Organisms rely on physiological (thermal tolerance) and behavioral (thermal preference) mechanisms to thermoregulate in response to environmental temperature. In addition, spatial heterogeneity in temperatures can select for local adaptation in thermal tolerance, thermal preference, or both. However, the concordance between thermal tolerance and preference across genotypes and sexes within species and across populations is greatly understudied. The house fly, Musca domestica, is a well-suited system to examine how genotype and environment interact to affect thermal tolerance and preference. Across multiple continents, house fly males from higher latitudes tend to carry the male-determining gene on the Y chromosome, whereas those from lower latitudes usually have the male determiner on the third chromosome. We tested whether these two male-determining chromosomes differentially affect thermal tolerance and preference as predicted by their geographical distributions. We identify effects of genotype and developmental temperature on male thermal tolerance and preference that are concordant with the natural distributions of the chromosomes, suggesting that temperature variation across the species range contributes to the maintenance of the polymorphism. In contrast, female thermal preference is bimodal and largely independent of congener male genotypes. These sexually dimorphic thermal preferences suggest that temperature-dependent mating dynamics within populations could further affect the distribution of the two chromosomes. Together, the differences in thermal tolerance and preference across sexes and male genotypes suggest that different selection pressures may affect the frequencies of the male-determining chromosomes across different spatial scales.

Different ways to play it cool: Transcriptomic analysis sheds light on different activity patterns of three amphipod species under long-term cold exposure

Species of littoral freshwater environments in regions with continental climate experience pronounced seasonal temperature changes. Coping with long cold winters and hot summers requires specific physiological and behavioural adaptations. Endemic amphipods of Lake Baikal, Eulimnogammarus verrucosus and Eulimnogammarus cyaneus, show high metabolic activity throughout the year; E. verrucosus even reproduces in winter. In contrast, the widespread Holarctic amphipod Gammarus lacustris overwinters in torpor. This study investigated the transcriptomic hallmarks of E. verrucosus, E. cyaneus and G. lacustris exposed to low water temperatures. Amphipods were exposed to 1.5°C and 12°C (corresponding to the mean winter and summer water temperatures, respectively, in the Baikal littoral) for one month. At 1.5°C, G. lacustris showed upregulation of ribosome biogenesis and mRNA processing genes, as well as downregulation of genes related to growth, reproduction and locomotor activity, indicating enhanced energy allocation to somatic maintenance. Our results suggest that the mitogen-activated protein kinase (MAPK) signalling pathway is involved in the preparation for hibernation; downregulation of the actin cytoskeleton pathway genes could relate to the observed low locomotor activity of G. lacustris at 1.5°C. The differences between the transcriptomes of E. verrucosus and E. cyaneus from the 1.5°C and 12°C exposures were considerably smaller than for G. lacustris. In E. verrucosus, cold-exposure triggered reproductive activity was indicated by upregulation of respective genes, whereas in E. cyaneus, genes related to mitochondria functioning were upregulated, indicating cold compensation in this species. Our data elucidate the molecular characteristics behind the different adaptations of amphipod species from the Lake Baikal area to winter conditions.

The UDP-Glycosyltransferase Family in Drosophila melanogaster: Nomenclature Update, Gene Expression and Phylogenetic Analysis

UDP-glycosyltransferases (UGTs) are important conjugation enzymes found in all kingdoms of life, catalyzing a sugar conjugation with small lipophilic compounds and playing a crucial role in detoxification and homeostasis. The UGT gene family is defined by a signature motif in the C-terminal domain where the uridine diphosphate (UDP)-sugar donor binds. UGTs have been identified in a number of insect genomes over the last decade and much progress has been achieved in characterizing their expression patterns and molecular functions. Here, we present an update of the complete repertoire of UGT genes in Drosophila melanogaster and provide a brief overview of the latest research in this model insect. A total of 35 UGT genes are found in the D. melanogaster genome, localized to chromosomes 2 and 3 with a high degree of gene duplications on the chromosome arm 3R. All D. melanogaster UGT genes have now been named in FlyBase according to the unified UGT nomenclature guidelines. A phylogenetic analysis of UGT genes shows lineage-specific gene duplications. Analysis of anatomical and induced gene expression patterns demonstrate that some UGT genes are differentially expressed in various tissues or after environmental treatments. Extended searches of UGT orthologs from 18 additional Drosophila species reveal a diversity of UGT gene numbers and composition. The roles of Drosophila UGTs identified to date are briefly reviewed, and include xenobiotic metabolism, nicotine resistance, olfaction, cold tolerance, sclerotization, pigmentation, and immunity. Together, the updated genomic information and research overview provided herein will aid further research in this developing field.

Small heat shock protein Hsp67Bc plays a significant role in Drosophila melanogaster cold stress tolerance

Hsp67Bc in Drosophila melanogaster is a member of the small heat shock protein family, the main function of which is to prevent the aggregation of misfolded or damaged proteins. Hsp67Bc interacts with Starvin and Hsp23, which are known to be a part of the cold stress response in the fly during the recovery phase. In this study, we investigated the role of the Hsp67Bc gene in the cold stress response. We showed that in adult Drosophila, Hsp67Bc expression increases after cold stress and decreases after 1.5 h of recovery, indicating the involvement of Hsp67Bc in short-term stress recovery. We also implemented a deletion in the D. melanogaster Hsp67Bc gene using imprecise excision of a P-element, and analysed the cold tolerance of Hsp67Bc-null mutants at different developmental stages. We found that Hsp67Bc-null homozygous flies are viable and fertile but display varying cold stress tolerance throughout the stages of ontogenesis: the survival after cold stress is slightly impaired in late third instar larvae, unaffected in pupae, and notably affected in adult females. Moreover, the recovery from chill coma is delayed in Hsp67Bc-null adults of both sexes. In addition, the deletion in the Hsp67Bc gene caused more prominent up-regulation of Hsp70 following cold stress, suggesting the involvement of Hsp70 in compensation of the lack of the Hsp67Bc protein. Taken together, our results suggest that Hsp67Bc is involved in the recovery of flies from a comatose state and contributes to the protection of the fruit fly from cold stress.

Integrating GWAS and Transcriptomics to Identify the Molecular Underpinnings of Thermal Stress Responses in Drosophila melanogaster

Thermal tolerance of an organism depends on both the ability to dynamically adjust to a thermal stress and preparatory developmental processes that enhance thermal resistance. However, the extent to which standing genetic variation in thermal tolerance alleles influence dynamic stress responses vs. preparatory processes is unknown. Here, using the model species Drosophila melanogaster, we used a combination of Genome Wide Association mapping (GWAS) and transcriptomic profiling to characterize whether genes associated with thermal tolerance are primarily involved in dynamic stress responses or preparatory processes that influence physiological condition at the time of thermal stress. To test our hypotheses, we measured the critical thermal minimum (CTmin) and critical thermal maximum (CTmax) of 100 lines of the Drosophila Genetic Reference Panel (DGRP) and used GWAS to identify loci that explain variation in thermal limits. We observed greater variation in lower thermal limits, with CTmin ranging from 1.81 to 8.60°C, while CTmax ranged from 38.74 to 40.64°C. We identified 151 and 99 distinct genes associated with CTmin and CTmax, respectively, and there was strong support that these genes are involved in both dynamic responses to thermal stress and preparatory processes that increase thermal resistance. Many of the genes identified by GWAS were involved in the direct transcriptional response to thermal stress (72/151 for cold; 59/99 for heat), and overall GWAS candidates were more likely to be differentially expressed than other genes. Further, several GWAS candidates were regulatory genes that may participate in the regulation of stress responses, and gene ontologies related to development and morphogenesis were enriched, suggesting many of these genes influence thermal tolerance through effects on development and physiological status. Overall, our results suggest that thermal tolerance alleles can influence both dynamic plastic responses to thermal stress and preparatory processes that improve thermal resistance. These results also have utility for directly comparing GWAS and transcriptomic approaches for identifying candidate genes associated with thermal tolerance.

Three Quantitative Trait Loci Explain More than 60% of Variation for Chill Coma Recovery Time in a Natural Population of Drosophila ananassae

Ectothermic species such as insects are particularly vulnerable to climatic fluctuations. Nevertheless, many insects that evolved and diversified in the tropics have successfully colonized temperate regions all over the globe. To shed light on the genetic basis of cold tolerance in such species, we conducted a quantitative trait locus (QTL) mapping experiment for chill coma recovery time (CCRT) in Drosophila ananassae, a cosmopolitan species that has expanded its range from tropical to temperate regions. We created a mapping population of recombinant inbred advanced intercross lines (RIAILs) from two founder strains with diverging CCRT phenotypes. The RIAILs were phenotyped for their CCRT and, together with the founder strains, genotyped for polymorphic markers with double-digest restriction site-associated DNA (ddRAD) sequencing. Using a hierarchical mapping approach that combined standard interval mapping and a multiple-QTL model, we mapped three QTL which altogether explained 64% of the phenotypic variance. For two of the identified QTL, we found evidence of epistasis. To narrow down the list of cold tolerance candidate genes, we cross-referenced the QTL intervals with genes that we previously identified as differentially expressed in response to cold in D. ananassae, and with thermotolerance candidate genes of D. melanogaster. Among the 58 differentially expressed genes that were contained within the QTL, GF15058 showed a significant interaction of the CCRT phenotype and gene expression. Further, we identified the orthologs of four D. melanogaster thermotolerance candidate genes, MtnA, klarsicht, CG5246 (D.ana/GF17132) and CG10383 (D.ana/GF14829) as candidates for cold tolerance in D. ananassae.

Cold acclimation triggers major transcriptional changes in Drosophila suzukii

BACKGROUND

Insects have the capacity to adjust their physiological mechanisms during their lifetime to promote cold tolerance and cope with sublethal thermal conditions, a phenomenon referred to as thermal acclimation. The spotted wing drosophila, Drosophila suzukii, is an invasive fruit pest that, like many other species, enhances its thermotolerance in response to thermal acclimation. However, little is known about the underlying mechanisms of this plastic response. Here, we promoted flies' cold tolerance by gradually increasing acclimation duration (i.e. pre-exposure from 2 h to 9 days at 10 °C), and then compared transcriptomic responses of cold hardy versus cold susceptible phenotypes using RNA sequencing.

RESULTS

Cold tolerance of D. suzukii increased with acclimation duration; the longer the acclimation, the higher the cold tolerance. Cold-tolerant flies that were acclimated for 9 days were selected for transcriptomic analyses. RNA sequencing revealed a total of 2908 differentially expressed genes: 1583 were up- and 1325 were downregulated in cold acclimated flies. Functional annotation revealed many enriched GO-terms among which ionic transport across membranes and signaling were highly represented in acclimated flies. Neuronal activity and carbohydrate metabolism were also enriched GO-terms in acclimated flies. Results also revealed many GO-terms related to oogenesis which were underrepresented in acclimated flies.

CONCLUSIONS

Involvement of a large cluster of genes related to ion transport in cold acclimated flies suggests adjustments in the capacity to maintain ion and water homeostasis. These processes are key mechanisms underlying cold tolerance in insects. Down regulation of genes related to oogenesis in cold acclimated females likely reflects that females were conditioned at 10 °C, a temperature that prevents oogenesis. Overall, these results help to understand the molecular underpinnings of cold tolerance acquisition in D. suzukii. These data are of importance considering that the invasive success of D. suzukii in diverse climatic regions relates to its high thermal plasticity.