ACD-containing chaperones reveal the divergent thermo-tolerance in penaeid shrimp

The α-crystallin domain-containing (ACD-containing) gene family, which includes typical small heat shock proteins (sHSPs), is the most ubiquitous and diverse family of putative chaperones in all organisms, including eukaryotes and prokaryotes. In the present study, approximately 54-117 ACD-containing genes were identified in five penaeid shrimp species, yielding a significant expansion in comparison with other crustaceans (generally 6-20 ACD-containing genes). Unlike typical sHSPs, which contain a single ACD domain, the ACD-containing genes of penaeid shrimp contain additional ACD domains (3-7 domains, in general), thus having a larger molecular weight and a more complex 3D structure. As indicated by the RNA-seq and qRT-PCR results, the ACD-containing genes of penaeid shrimp showed a strong response to high temperatures. Furthermore, heterologous expression and citrate synthase assays of three representative ACD-containing genes confirmed that their chaperone activity could enhance the thermo-tolerance of E. coli and prevent the aggregation of substrate proteins at high temperatures. Compared with penaeid shrimp species with a relatively low thermo-tolerance (Fenneropenaeus chinensis and Marsupenaeus japonicus), the species with high thermo-tolerance (Litopenaeus vannamei and Fenneropenaeus indicus) contained more ACD-containing genes due to tandem duplication and exhibited biased expression levels under high temperatures. This can explain the divergent thermo-tolerance of different penaeid shrimp species. In conclusion, the ACD-containing genes in penaeid shrimp could be assigned as new chaperones and contribute to their divergent thermo-tolerance phenotypes and adaptations to the ecological environment.

Loss of mitochondrial performance at high temperatures is correlated with upper thermal tolerance among populations of an intertidal copepod

Environmental temperatures have pervasive effects on the performance and tolerance of ectothermic organisms, and thermal tolerance limits likely play key roles underlying biogeographic ranges and responses to environmental change. Mitochondria are central to metabolic processes in eukaryotic cells, and these metabolic functions are thermally sensitive; however, potential relationships between mitochondrial function, thermal tolerance limits and local thermal adaptation in general remain unresolved. Loss of ATP synthesis capacity at high temperatures has recently been suggested as a mechanistic link between mitochondrial function and upper thermal tolerance limits. Here we use a common-garden experiment with seven locally adapted populations of intertidal copepods (Tigriopus californicus), spanning approximately 21.5° latitude, to assess genetically based variation in the thermal performance curves of maximal ATP synthesis rates in isolated mitochondria. These thermal performance curves displayed substantial variation among populations with higher ATP synthesis rates at lower temperatures (20-25 °C) in northern populations than in southern populations. In contrast, mitochondria from southern populations maintained ATP synthesis rates at higher temperatures than the temperatures that caused loss of ATP synthesis capacity in mitochondria from northern populations. Additionally, there was a tight correlation between the thermal limits of ATP synthesis and previously determined variation in upper thermal tolerance limits among populations. This suggests that mitochondria may play an important role in latitudinal thermal adaptation in T. californicus, and supports the hypothesis that loss of mitochondrial performance at high temperatures is linked to whole-organism thermal tolerance limits in this ectotherm.

Experimental evolution reveals the synergistic genomic mechanisms of adaptation to ocean warming and acidification in a marine copepod

Metazoan adaptation to global change relies on selection of standing genetic variation. Determining the extent to which this variation exists in natural populations, particularly for responses to simultaneous stressors, is essential to make accurate predictions for persistence in future conditions. Here, we identified the genetic variation enabling the copepod Acartia tonsa to adapt to experimental ocean warming, acidification, and combined ocean warming and acidification (OWA) over 25 generations of continual selection. Replicate populations showed a consistent polygenic response to each condition, targeting an array of adaptive mechanisms including cellular homeostasis, development, and stress response. We used a genome-wide covariance approach to partition the allelic changes into three categories: selection, drift and replicate-specific selection, and laboratory adaptation responses. The majority of allele frequency change in warming (57%) and OWA (63%) was driven by shared selection pressures across replicates, but this effect was weaker under acidification alone (20%). OWA and warming shared 37% of their response to selection but OWA and acidification shared just 1%, indicating that warming is the dominant driver of selection in OWA. Despite the dominance of warming, the interaction with acidification was still critical as the OWA selection response was highly synergistic with 47% of the allelic selection response unique from either individual treatment. These results disentangle how genomic targets of selection differ between single and multiple stressors and demonstrate the complexity that nonadditive multiple stressors will contribute to predictions of adaptation to complex environmental shifts caused by global change.

Competition for finite resources as coordination mechanism for morphogenesis: An evolutionary algorithm study of digital embryogeny

The standard view of embryogenesis is one of cooperation driven by the cells' shared genetics and evolutionary interests. However, numerous examples from developmental biology and agriculture reveal a surprising amount of competition among body cells, tissues, and organs for both metabolic and informational resources. To explain the existence of such competition we had hypothesized that evolution uses limiting "reservoirs" of resource molecules as a communication medium - a global scratchpad, to enable tissues across the body to coordinate growth. Here, we test this hypothesis via an evolutionary simulation of embryogeny in silico. Genomes encode state transition rules for cells, such as proliferation, differentiation, and resource use, enabling virtual embryos to develop a specific large-scale morphology. An evolutionary algorithm operates over the genomes, with fitness defined as a function of specific morphological requirements for the final embryo shape. We found that not only does such an algorithm rapidly discover rules for cellular behavior that reliably make embryos with specific anatomical properties, but that it discovers the strategy of using finite resources to coordinate development. Given the option of using finite or infinite reservoirs (which determine cells' ability to carry out specific actions), evolution preferentially uses finite reservoirs, which results in higher fitness and increased consistency (without needing direct selection for morphological invariance). We report aspects of anatomical, physiological/transcriptional, and genomic analysis of evolved virtual embryos that help understand how evolution can use competition among genetically identical subunits within a multicellular body to coordinate reliable, complex morphogenesis. Our results suggest that under some conditions, composite multi-scale systems will promote conflict and artificial scarcity for their components.

Too Hot to Eat: Wild and Lab-Bred Lymnaea stagnalis Differ in Feeding Response Following Repeated Heat Exposure

AbstractAcute extreme heat events are increasing in frequency and intensity. Understanding their effects on ectothermic organisms' homeostasis is both important and urgent. In this study we found that the exposure to an acute heat shock (30 °C for 1 hour) repeated for a seven-day period severely suppressed the feeding behavior of laboratory-inbred (W-strain) Lymnaea stagnalis, whereas the first-generation offspring of freshly collected wild (F₁ D-strain) snails raised and maintained under similar laboratory conditions did not show any alterations. The W-strain snails might have inadvertently been selected against heat tolerance since they were first brought into the laboratory many (∼70) years ago. We also posit that the F₁ D-strain snails do not perceive the heat shock as a sufficient stressor to alter their feeding response because their parental populations in wild environments have repeatedly experienced temperature fluctuations, thus becoming more tolerant and resilient to heat. The different responses exhibited by two strains of the same species highlight the importance of selecting the most appropriate strain for addressing questions about the impacts of global warming on organisms' physiology and behavior.

Spatial variation in high temperature-regulated gene expression predicts evolution of plasticity with climate change in the scarlet monkeyflower

A major way that organisms can adapt to changing environmental conditions is by evolving increased or decreased phenotypic plasticity. In the face of current global warming, more attention is being paid to the role of plasticity in maintaining fitness as abiotic conditions change over time. However, given that temporal data can be challenging to acquire, a major question is whether evolution in plasticity across space can predict adaptive plasticity across time. In growth chambers simulating two thermal regimes, we generated transcriptome data for western North American scarlet monkeyflowers (Mimulus cardinalis) collected from different latitudes and years (2010 and 2017) to test hypotheses about how plasticity in gene expression is responding to increases in temperature, and if this pattern is consistent across time and space. Supporting the genetic compensation hypothesis, individuals whose progenitors were collected from the warmer-origin northern 2017 descendant cohort showed lower thermal plasticity in gene expression than their cooler-origin northern 2010 ancestors. This was largely due to a change in response at the warmer (40°C) rather than cooler (20°C) treatment. A similar pattern of reduced plasticity, largely due to a change in response at 40°C, was also found for the cooler-origin northern versus the warmer-origin southern population from 2017. Our results demonstrate that reduced phenotypic plasticity can evolve with warming and that spatial and temporal changes in plasticity predict one another.

Nature versus nurture in heat stress induced learning between inbred and outbred populations of Lymnaea stagnalis

Changing environmental conditions often lead to microevolution of traits that are adaptive under the current selection pressure. Currently, one of the major selection pressures is the rise in temperatures globally that has a severe impact on the behavioral ecology of animals. However, the role of thermal stress on neuronal plasticity and memory formation is not well understood. Thermal tolerance and sensitivity to heat stress show variation across populations of the same species experiencing different thermal regimes. We used two populations of the pond snail Lymnaea stagnalis: one lab-bred W-snails and the other wild Delta snails to test heat shock induced learning and memory formation for the Garcia effect learning paradigm. In Garcia effect, a single pairing of a heat stressor (30 °C for 1h) with a novel taste results in a taste-specific negative hedonic shift lasting 24h as long-term memory (LTM) in lab bred W-snails. In this study we used a repeated heat stress procedure to test for increased or decreased sensitivity to the heat before testing for the Garcia effect. We found that lab-bred W-snails show increased sensitivity to heat stress after repeated heat exposure for 7days, leading to enhanced LTM for Garcia effect with only 15min of heat exposure instead of standard 1h. Surprisingly, the freshly collected wild snails do not show Garcia effect. Additionally, F1 generation of wild snails raised and maintained under laboratory conditions still retain their heat stress tolerance similar to their parents and do not show a Garcia effect under standard learning paradigm or even after repeated heat stressor. Thus, we found a differential effect of heat stress on memory formation in wild and lab bred snails. Most interestingly we also show that local environmental (temperature) conditions for one generation is not enough to alter thermal sensitivity in a wild population of L. stagnalis.

Land use and life history constrain adaptive genetic variation and reduce the capacity for climate change adaptation in turtles

BACKGROUND

Rapid anthropogenic climate change will require species to adapt to shifting environmental conditions, with successful adaptation dependent upon current patterns of genetic variation. While landscape genomic approaches allow for exploration of local adaptation in non-model systems, most landscape genomics studies of adaptive capacity are limited to exploratory identification of potentially important functional genes, often without a priori expectations as to the gene functions that may be most important for climate change responses. In this study, we integrated targeted sequencing of genes of known function and genotyping of single-nucleotide polymorphisms to examine spatial, environmental, and species-specific patterns of potential local adaptation in two co-occuring turtle species: the Blanding's turtle (Emydoidea blandingii) and the snapping turtle (Chelydra serpentina).

RESULTS

We documented divergent patterns of spatial clustering between neutral and putatively adaptive genetic variation in both species. Environmental associations varied among gene regions and between species, with stronger environmental associations detected for genes involved in stress response and for the more specialized Blanding's turtle. Land cover appeared to be more important than climate in shaping spatial variation in functional genes, indicating that human landscape alterations may affect adaptive capacity important for climate change responses.

CONCLUSIONS

Our study provides evidence that responses to climate change will be contingent on species-specific adaptive capacity and past history of exposure to human land cover change.

Allele-Specific Expression and Evolution of Gene Regulation Underlying Acute Heat Stress Response and Local Adaptation in the Copepod Tigriopus californicus

Geographic variation in environmental temperature can select for local adaptation among conspecific populations. Divergence in gene expression across the transcriptome is a key mechanism for evolution of local thermal adaptation in many systems, yet the genetic mechanisms underlying this regulatory evolution remain poorly understood. Here we examine gene expression in 2 locally adapted Tigriopus californicus populations (heat tolerant San Diego, SD, and less tolerant Santa Cruz, SC) and their F1 hybrids during acute heat stress response. Allele-specific expression (ASE) in F1 hybrids was used to determine cis-regulatory divergence. We found that the number of genes showing significant allelic imbalance increased under heat stress compared to unstressed controls. This suggests that there is significant population divergence in cis-regulatory elements underlying heat stress response. Specifically, the number of genes showing an excess of transcripts from the more thermal tolerant (SD) population increased with heat stress while that number of genes with an SC excess was similar in both treatments. Inheritance patterns of gene expression also revealed that genes displaying SD-dominant expression phenotypes increase in number in response to heat stress; that is, across loci, gene expression in F1's following heat stress showed more similarity to SD than SC, a pattern that was absent in the control treatment. The observed patterns of ASE and inheritance of gene expression provide insight into the complex processes underlying local adaptation and thermal stress response.

Heat tolerance and thermal preference of the copepod Tigriopus californicus are insensitive to ecologically relevant dissolved oxygen levels

Shifting climate patterns may impose novel combinations of abiotic conditions on animals, yet understanding of the present-day interactive effects of multiple stressors remains under-developed. We tested the oxygen and capacity limited thermal tolerance (OCLTT) hypothesis and quantified environmental preference of the copepod Tigriopus californicus, which inhabits rocky-shore splashpools where diel fluctuations of temperature and dissolved oxygen (DO) are substantial. Egg-mass bearing females were exposed to a 5 h heat ramp to peak temperatures of 34.1–38.0 °C crossed with each of four oxygen levels: 22, 30, 100 and 250% saturation (4.7–5.3, 5.3–6.4, 21.2–21.3, and 50.7–53.3 kPa). Survival decreased at higher temperatures but was independent of DO. The behavioral preference of females was quantified in seven combinations of gradients of both temperature (11–37 °C) and oxygen saturation (17–206% or 3.6–43.6 kPa). Females avoided high temperatures regardless of DO levels. This pattern was more pronounced when low DO coincided with high temperature. In uniform temperature treatments, the distribution shifted toward high DO levels, especially in uniform high temperature, confirming that Tigriopus can sense environmental pO₂. These results question the ecological relevance of OCLTT for Tigriopus and raise the possibility of microhabitat selection being used within splashpool environments to avoid physiologically stressful combinations of conditions.

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

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

Adaptation to climate change through genetic accommodation and assimilation of plastic phenotypes

Theory suggests that evolutionary changes in phenotypic plasticity could either hinder or facilitate evolutionary rescue in a changing climate. Nevertheless, the actual role of evolving plasticity in the responses of natural populations to climate change remains unresolved. Direct observations of evolutionary change in nature are rare, making it difficult to assess the relative contributions of changes in trait means versus changes in plasticity to climate change responses. To address this gap, this review explores several proxies that can be used to understand evolving plasticity in the context of climate change, including space for time substitutions, experimental evolution and tests for genomic divergence at environmentally responsive loci. Comparisons among populations indicate a prominent role for divergence in environmentally responsive traits in local adaptation to climatic gradients. Moreover, genomic comparisons among such populations have identified pervasive divergence in the regulatory regions of environmentally responsive loci. Taken together, these lines of evidence suggest that divergence in plasticity plays a prominent role in adaptation to climatic gradients over space, indicating that evolving plasticity is also likely to play a key role in adaptive responses to climate change through time. This suggests that genetic variation in plastic responses to the environment (G × E) might be an important predictor of species' vulnerabilities to climate-driven decline or extinction. This article is part of the theme issue 'The role of plasticity in phenotypic adaptation to rapid environmental change'.

Consequences of HSF knockdown on gene expression during the heat shock response in Tigriopus californicus

Although the existence of a cellular heat shock response is nearly universal, its relationship to organismal thermal tolerance is not completely understood. Many of the genes involved are known to be regulated by the highly conserved heat shock transcription factor-1 (HSF-1), yet the regulatory network is not fully characterized. Here, we investigated the role of HSF-1 in gene expression following thermal stress using knockdown of HSF-1 by RNA interference in the intertidal copepod Tigriopus californicus We observed some evidence for decreased transcription of heat shock protein genes following knockdown, supporting the widely acknowledged role of HSF-1 in the heat shock response. However, the majority of differentially expressed genes between the control and HSF-1 knockdown groups were upregulated, suggesting that HSF-1 normally functions to repress their expression. Differential expression observed in genes related to chitin and cuticle formation lends support to previous findings that these processes are highly regulated following heat stress. We performed a genome scan and identified a set of 396 genes associated with canonical heat shock elements. RNA-seq data did not find those genes to be more highly represented in our HSF-1 knockdown treatment, indicating that requirements for binding and interaction of HSF-1 with a given gene are not simply predicted by the presence of HSF-1 binding sites. Further study of the pathways implicated by these results and future comparisons among populations of T. californicus may help us understand the role and importance of HSF-1 in the heat shock response and, more broadly, in organismal thermal tolerance.

Novel signals of adaptive genetic variation in northwestern Atlantic cod revealed by whole-genome sequencing

Selection can create complex patterns of adaptive differentiation among populations in the wild that may be relevant to management. Atlantic cod in the Northwest Atlantic are at a fraction of their historical abundance and a lack of recovery within the Gulf of Maine has created concern regarding the misalignment of fisheries management structures with biological population structure. To address this and investigate genome-wide patterns of variation, we used low-coverage sequencing to perform a region-wide, whole-genome analysis of fine-scale population structure. We sequenced 306 individuals from 20 sampling locations in U.S. and Canadian waters, including the major spawning aggregations in the Gulf of Maine in addition to spawning aggregations from Georges Bank, southern New England, the eastern Scotian Shelf, and St. Pierre Bank. With genotype likelihoods estimated at almost 11 million loci, we found large differences in haplotype frequencies of previously described chromosomal inversions between Canadian and U.S. sampling locations and also among U.S. sampling locations. Our whole-genome resolution also revealed novel outlier peaks, some of which showed significant genetic differentiation among sampling locations. Comparisons between allochronic winter- and spring-spawning populations revealed highly elevated relative (FST ) and absolute (dxy ) genetic differentiation near genes involved in reproduction, particularly genes associated with the brain-pituitary-gonadal axis, which likely control timing of spawning, contributing to prezygotic isolation. We also found genetic differentiation associated with heat shock proteins and other genes of functional relevance, with complex patterns that may point to multifaceted selection pressures and local adaptation among spawning populations. We provide a high-resolution picture of U.S. Atlantic cod population structure, revealing greater complexity than is currently recognized in management. Our genome-scan approach likely underestimates the full suite of adaptive differentiation among sampling locations. Nevertheless, it should inform the revision of stock boundaries to preserve adaptive genetic diversity and evolutionary potential of cod populations.

Variation in developmental temperature alters adulthood plasticity of thermal tolerance in Tigriopus californicus

In response to environmental change, organisms rely on both genetic adaptation and phenotypic plasticity to adjust key traits that are necessary for survival and reproduction. Given the accelerating rate of climate change, plasticity may be particularly important. For organisms in warming aquatic habitats, upper thermal tolerance is likely to be a key trait, and many organisms express plasticity in this trait in response to developmental or adulthood temperatures. Although plasticity at one life stage may influence plasticity at another life stage, relatively little is known about this possibility for thermal tolerance. Here, we used locally adapted populations of the copepod Tigriopus californicus to investigate these potential effects in an intertidal ectotherm. We found that low latitude populations had greater critical thermal maxima (CT_max) than high latitude populations, and variation in developmental temperature altered CT_max plasticity in adults. After development at 25°C, CT_max was plastic in adults, whereas no adulthood plasticity in this trait was observed after development at 20°C. This pattern was identical across four populations, suggesting that local thermal adaptation has not shaped this effect among these populations. Differences in the capacities to maintain ATP synthesis rates and to induce heat shock proteins at high temperatures, two likely mechanisms of local adaptation in this species, were consistent with changes in CT_max owing to phenotypic plasticity, which suggests that there is likely mechanistic overlap between the effects of plasticity and adaptation. Together, these results indicate that developmental effects may have substantial impacts on upper thermal tolerance plasticity in adult ectotherms.

Accelerated evolution at chaperone promoters among Antarctic notothenioid fishes

BACKGROUND

Antarctic fishes of the Notothenioidei suborder constitutively upregulate multiple inducible chaperones, a highly derived adaptation that preserves proteostasis in extreme cold, and represent a system for studying the evolution of gene frontloading. We screened for Hsf1-binding sites, as Hsf1 is a master transcription factor of the heat shock response, and highly-conserved non-coding elements within proximal promoters of chaperone genes across 10 Antarctic notothens, 2 subpolar notothens, and 17 perciform fishes. We employed phylogenetic models of molecular evolution to determine whether (i) changes in motifs associated with Hsf1-binding and/or (ii) relaxed purifying selection or exaptation at ancestral cis-regulatory elements coincided with the evolution of chaperone frontloading in Antarctic notothens.

RESULTS

Antarctic notothens exhibited significantly fewer Hsf1-binding sites per bp at chaperone promoters than subpolar notothens and Serranoidei, the most closely-related suborder to Notothenioidei included in this study. 90% of chaperone promoters exhibited accelerated substitution rates among Antarctic notothens relative to other perciformes. The proportion of bases undergoing accelerated evolution (i) was significantly greater in Antarctic notothens than in subpolar notothens and Perciformes in 70% of chaperone genes and (ii) increased among bases that were more conserved among perciformes. Lastly, we detected evidence of relaxed purifying selection and exaptation acting on ancestrally conserved cis-regulatory elements in the Antarctic notothen lineage and its major branches.

CONCLUSION

A large degree of turnover has occurred in Notothenioidei at chaperone promoter regions that are conserved among perciform fishes following adaptation to the cooling of the Southern Ocean. Additionally, derived reductions in Hsf1-binding site frequency suggest cis-regulatory modifications to the classical heat shock response. Of note, turnover events within chaperone promoters were less frequent in the ancestral node of Antarctic notothens relative to younger Antarctic lineages. This suggests that cis-regulatory divergence at chaperone promoters may be greater between Antarctic notothen lineages than between subpolar and Antarctic clades. These findings demonstrate that strong selective forces have acted upon cis-regulatory elements of chaperone genes among Antarctic notothens.

The molecular basis for the neofunctionalization of the juvenile hormone esterase duplication in Drosophila

The Drosophila melanogaster enzymes juvenile hormone esterase (DmJHE) and its duplicate, DmJHEdup, present ideal examples for studying the structural changes involved in the neofunctionalization of enzyme duplicates. DmJHE is a hormone esterase with precise regulation and highly specific activity for its substrate, juvenile hormone. DmJHEdup is an odorant degrading esterase (ODE) responsible for processing various kairomones in antennae. Our phylogenetic analysis shows that the JHE lineage predates the hemi/holometabolan split and that several duplications of JHEs have been templates for the evolution of secreted β-esterases such as ODEs through the course of insect evolution. Our biochemical comparisons further show that DmJHE has sufficient substrate promiscuity and activity against odorant esters for a duplicate to evolve a general ODE function against a range of mid-long chain food esters, as is shown in DmJHEdup. This substrate range complements that of the only other general ODE known in this species, Esterase 6. Homology models of DmJHE and DmJHEdup enabled comparisons between each enzyme and the known structures of a lepidopteran JHE and Esterase 6. Both JHEs showed very similar active sites despite low sequence identity (30%). Both ODEs differed drastically from the JHEs and each other, explaining their complementary substrate ranges. A small number of amino acid changes are identified that may have been involved in the early stages of the neofunctionalization of DmJHEdup. Our results provide key insights into the process of neofunctionalization and the structural changes that can be involved.