Genetically Encoded Lizard Color Divergence for Camouflage and Thermoregulation.

Local adaptation is critical in speciation and evolution, yet comprehensive studies on proximate and ultimate causes of local adaptation are generally scarce. Here, we integrated field ecological experiments, genome sequencing, and genetic verification to demonstrate both driving forces and molecular mechanisms governing local adaptation of body coloration in a lizard from the Qinghai-Tibet Plateau. We found dark lizards from the cold meadow population had lower spectrum reflectance but higher melanin contents than light counterparts from the warm dune population. Additionally, the colorations of both dark and light lizards facilitated the camouflage and thermoregulation in their respective microhabitat simultaneously. More importantly, by genome resequencing analysis, we detected a novel mutation in Tyrp1 that underpinned this color adaptation. The allele frequencies at the site of SNP 459# in the gene of Tyrp1 are 22.22% G/C and 77.78% C/C in dark lizards and 100% G/G in light lizards. Model-predicted structure and catalytic activity showed that this mutation increased structure flexibility and catalytic activity in enzyme TYRP1, and thereby facilitated the generation of eumelanin in dark lizards. The function of the mutation in Tyrp1 was further verified by more melanin contents and darker coloration detected in the zebrafish injected with the genotype of Tyrp1 from dark lizards. Therefore, our study demonstrates that a novel mutation of a major melanin-generating gene underpins skin color variation co-selected by camouflage and thermoregulation in a lizard. The resulting strong selection may reinforce adaptive genetic divergence and enable the persistence of adjacent populations with distinct body coloration.

The impacts of thermal heterogeneity across microhabitats on post-settlement selection of intertidal mussels.

Rapid genetic selection is critical for allowing natural populations to adapt to different thermal environments such as those that occur across intertidal microhabitats with high degrees of thermal heterogeneity. To address the question of how thermal regimes influence selection and adaptation in the intertidal black mussel Mytilisepta virgata, we continuously recorded environmental temperatures in both tidal pools and emergent rock microhabitats and then assessed genetic differentiation, gene expression patterns, RNA editing level, and cardiac performance. Our results showed that the subpopulations in the tidal pool and on emergent rocks had different genetic structures and exhibited different physiological and molecular responses to high-temperature stress. These results indicate that environmental heterogeneity across microhabitats is important for driving genetic differentiation and shed light on the importance of post-settlement selection for adaptively modifying the genetic composition and thermal responses of these intertidal mussels.

Exploring the adaptability of the secondary structure of mRNA to temperature in intertidal snails based on SHAPE experiments.

RNA-based thermal regulation is an important strategy for organisms to cope with temperature changes. Inhabiting the intertidal rocky shore, a key interface of the ocean, atmosphere and terrestrial environments, intertidal species have developed variable thermal adaptation mechanisms; however, adaptions at the RNA level remain largely uninvestigated. To examine the relationship between mRNA structural stability and species distribution, in the present study, the secondary structure of cytosolic malate dehydrogenase (cMDH) mRNA of Echinolittorina malaccana, Echinolittorina radiata and Littorina brevicula was determined using selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE), and the change in folding free energy of formation (ΔGfold) was calculated. The results showed that ΔGfold increased as the temperature increased. The difference in ΔGfold (ΔΔGfold) between two specific temperatures (25 versus 0°C, 37 versus 0°C and 57 versus 0°C) differed among the three species, and the ΔΔGfold value of E. malaccana was significantly lower than those of E. radiata and L. brevicula. The number of stems of cMDH mRNA of the snails decreased with increasing temperature, and the breakpoint temperature of E. malaccana was the highest among these. The number of loops was also reduced with increasing temperature, while the length of the loop structure increased accordingly. Consequently, these structural changes can potentially affect the translational efficiency of mRNA. These results imply that there were interspecific differences in the thermal stability of RNA secondary structures in intertidal snails, and these differences may be related to snail distribution.

An integrated, multi-level analysis of thermal effects on intertidal molluscs for understanding species distribution patterns.

Elucidating the physiological mechanisms that underlie thermal stress and discovering how species differ in capacities for phenotypic acclimatization and evolutionary adaptation to this stress is critical for understanding current latitudinal and vertical distribution patterns of species and for predicting their future state in a warming world. Such mechanistic analyses require careful choice of study systems (species and temperature-sensitive traits) and design of laboratory experiments that reflect the complexities of in situ conditions. Here, we critically review a wide range of studies of intertidal molluscs that provide mechanistic accounts of thermal effects across all levels of biological organization - behavioural, organismal, organ level, cellular, molecular, and genomic - and show how temperature-sensitive traits govern distribution patterns and capacities for coping with thermal stress. Comparisons of congeners from different thermal habitats are especially effective means for identifying adaptive variation. We employ these mechanistic analyses to illustrate how species differ in the severity of threats posed by rising temperature. Counterintuitively, we show that some of the most heat-tolerant species may be most threatened by increases in temperatures because of their small thermal safety margins and minimal abilities to acclimatize to higher temperatures. We discuss recent molecular biological and genomic studies that provide critical foundations for understanding the types of evolutionary changes in protein structure, RNA secondary structure, genome content, and gene expression capacities that underlie adaptation to temperature. Duplication of stress-related genes, as found in heat-tolerant molluscs, may provide enhanced capacity for coping with higher temperatures. We propose that the anatomical, behavioural, physiological, and genomic diversity found among intertidal molluscs, which commonly are of critical importance and high abundance in these ecosystems, makes this group of animals a highly appropriate study system for addressing questions about the mechanistic determinants of current and future distribution patterns of intertidal organisms.

Thermal adaptation of mRNA secondary structure: stability versus lability.

Macromolecular function commonly involves rapidly reversible alterations in three-dimensional structure (conformation). To allow these essential conformational changes, macromolecules must possess higher order structures that are appropriately balanced between rigidity and flexibility. Because of the low stabilization free energies (marginal stabilities) of macromolecule conformations, temperature changes have strong effects on conformation and, thereby, on function. As is well known for proteins, during evolution, temperature-adaptive changes in sequence foster retention of optimal marginal stability at a species' normal physiological temperatures. Here, we extend this type of analysis to messenger RNAs (mRNAs), a class of macromolecules for which the stability-lability balance has not been elucidated. We employ in silico methods to determine secondary structures and estimate changes in free energy of folding (ΔGfold) for 25 orthologous mRNAs that encode the enzyme cytosolic malate dehydrogenase in marine mollusks with adaptation temperatures spanning an almost 60 °C range. The change in free energy that occurs during formation of the ensemble of mRNA secondary structures is significantly correlated with adaptation temperature: ΔGfold values are all negative and their absolute values increase with adaptation temperature. A principal mechanism underlying these adaptations is a significant increase in synonymous guanine + cytosine substitutions with increasing temperature. These findings open up an avenue of exploration in molecular evolution and raise interesting questions about the interaction between temperature-adaptive changes in mRNA sequence and in the proteins they encode.

Physiological determinants of biogeography: The importance of metabolic depression to heat tolerance.

A quantitative understanding of physiological thermal responses is vital for forecasting species distributional shifts in response to climate change. Many studies have focused on metabolic rate as a global metric for analyzing the sublethal effects of changing environments on physiology. Thermal performance curves (TPCs) have been suggested as a viable analytical framework, but standard TPCs may not fully capture physiological responses, due in part to failure to consider the process of metabolic depression. We derived a model based on the nonlinear regression of biological temperature-dependent rate processes and built a heart rate data set for 26 species of intertidal molluscs distributed from 33°S to ~40°N. We then calculated physiological thermal performance limits with continuous heating using T 1 / 2 H , the temperature at which heart rate is decreased to 50% of the maximal rate, as a more realistic measure of upper thermal limits. Results indicate that heat-induced metabolic depression of cardiac performance is a common adaptive response that allows tolerance of harsh environments. Furthermore, our model accounted for the high inter-individual variability in the shape of cardiac TPCs. We then used these TPCs to calculate physiological thermal safety margins (pTSM), the difference between the maximal operative temperature (95th percentile of field temperatures) and T 1 / 2 H of each individual. Using pTSMs, we developed a physiological species distribution model (pSDM) to forecast future geographic distributions. pSDM results indicate that climate-induced species range shifts are potentially less severe than predicted by a simple correlative SDM. Species with metabolic depression below the optimum temperature will be more thermal resistant at their warm trailing edges. High intraspecific variability further suggests that models based on species-level vulnerability to environmental change may be problematic. This multi-scale, mechanistic understanding that incorporates metabolic depression and inter-individual variability in thermal response enables better predictions about the relationship between thermal stress and species distributions.

Temperature adaptations of the thermophilic snail Echinolittorina malaccana: insights from metabolomic analysis.

The periwinkle snail Echinolittorina malaccana, for which the upper lethal temperature is near 55°C, is one of the most heat-tolerant eukaryotes known. We conducted a multi-level investigation - including cardiac physiology, enzyme activity, and targeted and untargeted metabolomic analyses - that elucidated a spectrum of adaptations to extreme heat in this organism. All systems examined showed heat intensity-dependent responses. Under moderate heat stress (37-45°C), the snail depressed cardiac activity and entered a state of metabolic depression. The global metabolomic and enzymatic analyses revealed production of metabolites characteristic of oxygen-independent pathways of ATP generation (lactate and succinate) in the depressed metabolic state, which suggests that anaerobic metabolism was the main energy supply pathway under heat stress (37-52°C). The metabolomic analyses also revealed alterations in glycerophospholipid metabolism under extreme heat stress (52°C), which likely reflected adaptive changes to maintain membrane structure. Small-molecular-mass organic osmolytes (glycine betaine, choline and carnitine) showed complex changes in concentration that were consistent with a role of these protein-stabilizing solutes in protection of the proteome under heat stress. This thermophilic species can thus deploy a wide array of adaptive strategies to acclimatize to extremely high temperatures.

Comparing mutagenesis and simulations as tools for identifying functionally important sequence changes for protein thermal adaptation.

Comparative studies of orthologous proteins of species evolved at different temperatures have revealed consistent patterns of temperature-related variation in thermal stabilities of structure and function. However, the precise mechanisms by which interspecific variations in sequence foster these adaptive changes remain largely unknown. Here, we compare orthologs of cytosolic malate dehydrogenase (cMDH) from marine molluscs adapted to temperatures ranging from -1.9 °C (Antarctica) to ∼55 °C (South China coast) and show how amino acid usage in different regions of the enzyme (surface, intermediate depth, and protein core) varies with adaptation temperature. This eukaryotic enzyme follows some but not all of the rules established in comparisons of archaeal and bacterial proteins. To link the effects of specific amino acid substitutions with adaptive variations in enzyme thermal stability, we combined site-directed mutagenesis (SDM) and in vitro protein experimentation with in silico mutagenesis using molecular dynamics simulation (MDS) techniques. SDM and MDS methods generally but not invariably yielded common effects on protein stability. MDS analysis is shown to provide insights into how specific amino acid substitutions affect the conformational flexibilities of mobile regions (MRs) of the enzyme that are essential for binding and catalysis. Whereas these substitutions invariably lie outside of the MRs, they effectively transmit their flexibility-modulating effects to the MRs through linked interactions among surface residues. This discovery illustrates that regions of the protein surface lying outside of the site of catalysis can help establish an enzyme's thermal responses and foster evolutionary adaptation of function.

Structural flexibility and protein adaptation to temperature: Molecular dynamics analysis of malate dehydrogenases of marine molluscs.

Orthologous proteins of species adapted to different temperatures exhibit differences in stability and function that are interpreted to reflect adaptive variation in structural "flexibility." However, quantifying flexibility and comparing flexibility across proteins has remained a challenge. To address this issue, we examined temperature effects on cytosolic malate dehydrogenase (cMDH) orthologs from differently thermally adapted congeners of five genera of marine molluscs whose field body temperatures span a range of ∼60 °C. We describe consistent patterns of convergent evolution in adaptation of function [temperature effects on KM of cofactor (NADH)] and structural stability (rate of heat denaturation of activity). To determine how these differences depend on flexibilities of overall structure and of regions known to be important in binding and catalysis, we performed molecular dynamics simulation (MDS) analyses. MDS analyses revealed a significant negative correlation between adaptation temperature and heat-induced increase of backbone atom movements [root mean square deviation (rmsd) of main-chain atoms]. Root mean square fluctuations (RMSFs) of movement by individual amino acid residues varied across the sequence in a qualitatively similar pattern among orthologs. Regions of sequence involved in ligand binding and catalysis-termed mobile regions 1 and 2 (MR1 and MR2), respectively-showed the largest values for RMSF. Heat-induced changes in RMSF values across the sequence and, importantly, in MR1 and MR2 were greatest in cold-adapted species. MDS methods are shown to provide powerful tools for examining adaptation of enzymes by providing a quantitative index of protein flexibility and identifying sequence regions where adaptive change in flexibility occurs.

Heat-resistant cytosolic malate dehydrogenases (cMDHs) of thermophilic intertidal snails (genus Echinolittorina): protein underpinnings of tolerance to body temperatures reaching 55°C.

Snails of the genus Echinolittorina are among the most heat-tolerant animals; they experience average body temperatures near 41-44°C in summer and withstand temperatures up to at least 55°C. Here, we demonstrate that heat stability of function (indexed by the Michaelis-Menten constant of the cofactor NADH, KMNADH) and structure (indexed by rate of denaturation) of cytosolic malate dehydrogenases (cMDHs) of two congeners (E. malaccana and E. radiata) exceeds values previously found for orthologs of this protein from less thermophilic species. The ortholog of E. malaccana is more heat stable than that of E. radiata, in keeping with the congeners' thermal environments. Only two inter-congener differences in amino acid sequence in these 332 residue proteins were identified. In both cases (positions 48 and 114), a glycine in the E. malaccana ortholog is replaced by a serine in the E. radiata protein. To explore the relationship between structure and function and to characterize how amino acid substitutions alter stability of different regions of the enzyme, we used molecular dynamics simulation methods. These computational methods allow determination of thermal effects on fine-scale movements of protein components, for example, by estimating the root mean square deviation in atom position over time and the root mean square fluctuation for individual residues. The minor changes in amino acid sequence favor temperature-adaptive change in flexibility of regions in and around the active sites. Interspecific differences in effects of temperature on fine-scale protein movements are consistent with the differences in thermal effects on binding and rates of heat denaturation.

Effects of dietary cadmium on growth, antioxidants and bioaccumulation of sea cucumber (Apostichopus japonicus) and influence of dietary vitamin C supplementation.

The effects of dietary cadmium (Cd) supplementation on growth, antioxidant capacity and accumulation of Cd in tissues (body wall, digestive tracts, and respiratory tree) of sea cucumber, Apostichopus japonicus, exposed to sub-chronic concentrations (0, 10, 50, 100, and 500mg Cd/kg dry weight) of Cd were investigated. In addition, the potential protective effects of vitamin C (L-ascorbic acid, AsA) against the effects of Cd on sea cucumbers were investigated. Sea cucumbers were exposed to dietary Cd for 30 days, after which another group of healthy sea cucumbers was supplied diet supplemented with mixed Cd and AsA for another 30 days. Cd exposure for 30 days resulted in increased Cd accumulation in tissues of sea cucumbers with exposure time and concentration. The order of Cd accumulation in organs was digestive tracts>respiratory tree>body wall. On day 30, the body weight gain (BWG) and specific growth rate (SGR) decreased significantly (P<0.05) in the 500mg Cd/kg treatment. Superoxide dismutase (SOD) activity, glutathione peroxidase (GSH-Px) activity and catalase (CAT) activity in the coelomic fluid of sea cucumbers decreased with increasing dietary Cd concentration, but malondialdehyde (MDA) content in the coelomic fluid increased. Providing diet supplemented with Cd and AsA indicated that although sea cucumbers exhibited signs of Cd toxicity, no death occurred in response to 50mg Cd/kg for 30 days. Based on these findings, five treatments were provided: 50mg Cd/kg+0mg AsA/kg, 50mg Cd/kg+ 3000mg AsA/kg, 50mg Cd/kg+ 5000mg AsA/kg, 50mg Cd/kg+10,000mg AsA/kg, and 50mg Cd/kg+15,000mg AsA/kg. The BWG and SGR of sea cucumbers fed the AsA supplemented diet mixed with Cd increased. Additionally, MDA levels in coelomic fluid were negatively correlated with dietary AsA levels, while antioxidant capacities (SOD, GSH-Px and CAT) were positively correlated with dietary AsA levels. Moreover, Cd accumulation in tissues decreased in response to dietary AsA supplementation of treatments. Overall, antioxidant capacity and bioaccumulation in sea cucumber was found to decrease and be induced in response to Cd, but vitamin C mitigated these effects, with 5000mg AsA/kg providing the optimum protection against 50mg/kg Cd.

The effects of dietary lead on growth, bioaccumulation and antioxidant capacity in sea cucumber, Apostichopus japonicus.

Three different diets amended with lead nitrate [Pb(NO3)2] (100, 500 and 1000mg Pb/kg dry weight) and a Pb-free control diet (1.03mg Pb/kg dry weight) were fed to sea cucumber (Apostichopus japonicus) for 30 days. The patterns of Pb accumulation over time were determined in various tissues (body wall, intestine and respiratory tree), as well as growth performance and antioxidant enzymes activities. Pb accumulation in body wall and intestine increased with time in all dietary Pb treatments. When fed the highest Pb diet, the body wall exhibited the greatest Pb burden (16.37mg Pb/kg tissue wet weight), while Pb content in the intestine (2.68mg Pb/kg tissue wet weight) and the respiratory tree (1.78mg Pb/kg tissue wet weight) were lower than Pb content in the body wall by day 30. The body weight gain (BWG), specific growth rate (SGR) and survival rate (SR) had not been affected by 30 days oral administration of Pb supplemented diet. However, the antioxidant enzymes activities [superoxide dismutase (SOD) and glutathione peroxidase (GSH-Px)] of test groups were lower than control group in body wall and malondialdehyde (MDA) concentration in the body wall was opposite after 30 days in sea cucumbers. In summary, this work reports toxic effects in sea cucumber, A. japonicus, after dietary exposure to Pb.

Effects of dietary vitamin C supplementation on lead-treated sea cucumbers, Apostichopus japonicus.

This study aimed to determine the effects of dietary supplementation with vitamin C (l-ascorbic acid, AsA) on lead (Pb) accumulation and toxicity in sea cucumber, Apostichopus japonicus. Three hundred sea cucumbers (10.02±0.02g) fed a basal diet containing 100mg Pb/kg were subjected to 5 levels of l-ascorbate-2-polyphosphate (LAPP) supplementation to achieve 5 AsA treatment levels of 0, 2727, 4630, 9171, 13,893mg AsA/kg. After 30 days, specific growth rate (SGR) and body weight gain (BWG) of supplementation groups were significantly (P<0.05) higher than the control group. Pb contents of the intestine significantly (P<0.05) decreased from 20 days to 30 days in all supplementation groups. Pb contents of the respiratory tree in sea cucumbers in the 2727 and 4630mg AsA/kg treatment groups significantly (P<0.05) decreased from 10 days to 30 days. Superoxide dismutase (SOD) activity, glutathione peroxidase (GSH-Px) activity and total antioxidant capacity (T-AOC) of the sea cucumber body wall in supplementation groups increased after 30 days. Malondialdehyde (MDA) content of the body wall significantly (P<0.05) declined with increasing AsA level. Therefore, it can be demonstrated that vitamin C could mitigate the effects of Pb on sea cucumber and the optimum levels ranged from 2727mg AsA/kg to 4630mg AsA/kg when Pb levels were 100mg/kg.