Low Temperature Adaptation Is Not the Opposite Process of High Temperature Adaptation in Terms of Changes in Amino Acid Composition

Genetic response of a perennial grass to warm and wet environments interacts and is associated with trait means as well as plasticity

The potential for rapid evolution is an important mechanism allowing species to adapt to changing climatic conditions. Although such potential has been largely studied in various short-lived organisms, to what extent we can observe similar patterns in long-lived plant species, which often dominate natural systems, is largely unexplored. We explored the potential for rapid evolution in Festuca rubra, a long-lived grass with extensive clonal growth dominating in alpine grasslands. We used a field sowing experiment simulating expected climate change in our model region. Specifically, we exposed seeds from five independent seed sources to novel climatic conditions by shifting them along a natural climatic grid and explored the genetic profiles of established seedlings after 3 years. Data on genetic profiles of plants selected under different novel conditions indicate that different climate shifts select significantly different pools of genotypes from common seed pools. Increasing soil moisture was more important than increasing temperature or the interaction of the two climatic factors in selecting pressure. This can indicate negative genetic interaction in response to the combined effects or that the effects of different climates are interactive rather than additive. The selected alleles were found in genomic regions, likely affecting the function of specific genes or their expression. Many of these were also linked to morphological traits (mainly to trait plasticity), suggesting these changes may have a consequence on plant performance. Overall, these data indicate that even long-lived plant species may experience strong selection by climate, and their populations thus have the potential to rapidly adapt to these novel conditions.

Metagenomic Thermometer

Various microorganisms exist in environments, and each of them has its optimal growth temperature (OGT). The relationship between genomic information and OGT of each species has long been studied, and one such study revealed that OGT of prokaryotes can be accurately predicted based on the fraction of seven amino acids (IVYWREL) among all encoded amino-acid sequences in its genome. Extending this discovery, we developed a 'Metagenomic Thermometer' as a means of predicting environmental temperature based on metagenomic sequences. Temperature prediction of diverse environments using publicly available metagenomic data revealed that the Metagenomic Thermometer can predict environmental temperatures with small temperature changes and little influx of microorganisms from other environments. The accuracy of the Metagenomic Thermometer was also confirmed by a demonstration experiment using an artificial hot water canal. The Metagenomic Thermometer was also applied to human gut metagenomic samples, yielding a reasonably accurate value for human body temperature. The result further suggests that deep body temperature determines the dominant lineage of the gut community. Metagenomic Thermometer provides a new insight into temperature-driven community assembly based on amino-acid composition rather than microbial taxa.

Multi-omics for studying and understanding polar life

Polar ecosystems are experiencing amongst the most rapid rates of regional warming on Earth. Here, we discuss 'omics' approaches to investigate polar biodiversity, including the current state of the art, future perspectives and recommendations. We propose a community road map to generate and more fully exploit multi-omics data from polar organisms. These data are needed for the comprehensive evaluation of polar biodiversity and to reveal how life evolved and adapted to permanently cold environments with extreme seasonality. We argue that concerted action is required to mitigate the impact of warming on polar ecosystems via conservation efforts, to sustainably manage these unique habitats and their ecosystem services, and for the sustainable bioprospecting of novel genes and compounds for societal gain.

The methanogen core and pangenome: conservation and variability across biology's growth temperature extremes

Temperature is a key variable in biological processes. However, a complete understanding of biological temperature adaptation is lacking, in part because of the unique constraints among different evolutionary lineages and physiological groups. Here we compared the genomes of cultivated psychrotolerant and thermotolerant methanogens, which are physiologically related and span growth temperatures from -2.5°C to 122°C. Despite being phylogenetically distributed amongst three phyla in the archaea, the genomic core of cultivated methanogens comprises about one-third of a given genome, while the genome fraction shared by any two organisms decreases with increasing phylogenetic distance between them. Increased methanogenic growth temperature is associated with reduced genome size, and thermotolerant organisms-which are distributed across the archaeal tree-have larger core genome fractions, suggesting that genome size is governed by temperature rather than phylogeny. Thermotolerant methanogens are enriched in metal and other transporters, and psychrotolerant methanogens are enriched in proteins related to structure and motility. Observed amino acid compositional differences between temperature groups include proteome charge, polarity and unfolding entropy. Our results suggest that in the methanogens, shared physiology maintains a large, conserved genomic core even across large phylogenetic distances and biology's temperature extremes.

Full-Length Transcriptome Comparison Provides Novel Insights into the Molecular Basis of Adaptation to Different Ecological Niches of the Deep-Sea Hydrothermal Vent in Alvinocaridid Shrimps

The deep-sea hydrothermal vent ecosystem is one of the extreme chemoautotrophic environments. Shinkaicaris leurokolos Kikuchi and Hashimoto, 2000, and Alvinocaris longirostris Kikuchi and Ohta, 1995, are typically co-distributed and closely related alvinocaridid shrimps in hydrothermal vent areas with different ecological niches, providing an excellent model for studying the adaptive evolution mechanism of animals in the extreme deep-sea hydrothermal vent environment. The shrimp S. leurokolos lives in close proximity to the chimney vent discharging high-temperature fluid, while A. longirostris inhabits the peripheral areas of hydrothermal vents. In this study, full-length transcriptomes of S. leurokolos and A. longirostris were generated using a combination of single-molecule real-time (SMRT) and Illumina RNA-seq technology. Expression analyses of the transcriptomes showed that among the top 30% of highly expressed genes of each species, more genes related to sulfide and heavy metal metabolism (sulfide: quinone oxidoreductase, SQR; persulfide dioxygenase, ETHE1; thiosulfate sulfurtransferase, TST, and ferritin, FRI) were specifically highly expressed in S. leurokolos, while genes involved in maintaining epibiotic bacteria or pathogen resistance (beta-1,3-glucan-binding protein, BGBP; endochitinase, CHIT; acidic mammalian chitinase, CHIA, and anti-lipopolysaccharide factors, ALPS) were highly expressed in A. longirostris. Gene family expansion analysis revealed that genes related to anti-oxidant metabolism (cytosolic manganese superoxide dismutase, SODM; glutathione S-transferase, GST, and glutathione peroxidase, GPX) and heat stress (heat shock cognate 70 kDa protein, HSP70 and heat shock 70 kDa protein cognate 4, HSP7D) underwent significant expansion in S. leurokolos, while CHIA and CHIT involved in pathogen resistance significantly expanded in A. longirostris. Finally, 66 positively selected genes (PSGs) were identified in the vent shrimp S. leurokolos. Most of the PSGs were involved in DNA repair, antioxidation, immune defense, and heat stress response, suggesting their function in the adaptive evolution of species inhabiting the extreme vent microhabitat. This study provides abundant genetic resources for deep-sea invertebrates, and is expected to lay the foundation for deep decipherment of the adaptive evolution mechanism of shrimps in a deep-sea chemosynthetic ecosystem based on further whole-genome comparison.

Mitogenomics provides new insights into the phylogenetic relationships and evolutionary history of deep-sea sea stars (Asteroidea)

The deep sea (> 200 m) is considered as the largest and most remote biome, which characterized by low temperatures, low oxygen level, scarce food, constant darkness, and high hydrostatic pressure. The sea stars (class Asteroidea) are ecologically important and diverse echinoderms in all of the world’s oceans, occurring from the intertidal to the abyssal zone (to about 6000 m). To date, the phylogeny of the sea stars and the relationships of deep-sea and shallow water groups have not yet been fully resolved. Here, we recovered five mitochondrial genomes of deep-sea asteroids. The A+T content of the mtDNA in deep-sea asteroids were significantly higher than that of the shallow-water groups. The gene orders of the five new mitogenomes were identical to that of other asteroids. The phylogenetic analysis showed that the orders Valvatida, Paxillosida, Forcipulatida are paraphyletic. Velatida was the sister order of all the others and then the cladeValvatida-Spinulosida-Paxillosida-Notomyotida versus Forcipulatida-Brisingida. Deep-sea asteroids were nested in different lineages, instead of a well-supported clade. The tropical Western Pacific was suggested as the original area of asteroids, and the temperate water was initially colonized with asteroids by the migration events from the tropical and cold water. The time-calibrated phylogeny showed that Asteroidea originated during Devonian-Carboniferous boundary and the major lineages of Asteroidea originated during Permian–Triassic boundary. The divergence between the deep-sea and shallow-water asteroids coincided approximately with the Triassic-Jurassic extinction. Total 29 positively selected sites were detected in fifteen mitochondrial genes of five deep-sea lineages, implying a link between deep-sea adaption and mitochondrial molecular biology in asteroids.

CCA-addition in the cold: Structural characterization of the psychrophilic CCA-adding enzyme from the permafrost bacterium Planococcus halocryophilus

•

A high-resolution structure of a psychrophilic RNA polymerase contributes to our knowledge of cold adaptation.

•

While catalytic core motifs are conserved, at least one shows cold adaptation.

•

Loss of helix-capping increases structural flexibility in a catalytic core motif.

•

Overall reduction of alpha-helical elements appears as a strategy for cold adaptation.

CCA-adding enzymes are highly specific RNA polymerases that add and maintain the sequence C-C-A at tRNA 3‘-ends. Recently, we could reveal that cold adaptation of such a polymerase is not only achieved at the expense of enzyme stability, but also at the cost of polymerization fidelity. Enzymes from psychrophilic organisms usually show an increased structural flexibility to enable catalysis at low temperatures. Here, polymerases face a dilemma, as there is a discrepancy between the need for a tightly controlled flexibility during polymerization and an increased flexibility as strategy for cold adaptation. Based on structural and biochemical analyses, we contribute to clarify the cold adaptation strategy of the psychrophilic CCA-adding enzyme from Planococcus halocryophilus, a gram-positive bacterium thriving in the arctic permafrost at low temperatures down to −15 °C. A comparison with the closely related enzyme from the thermophilic bacterium Geobacillus stearothermophilus reveals several features of cold adaptation - a significantly reduced amount of alpha-helical elements in the C-terminal tRNA-binding region and a structural adaptation in one of the highly conserved catalytic core motifs located in the N-terminal catalytic core of the enzyme.

Novel perspectives of environmental proteomics

The connection among genome expression, proteome alteration, metabolism regulation and phenotype change under environmental stresses is very vague. It is a tough task for the traditional research approaches to reveal the related scientific mechanisms of the above connection at molecular and systematic levels. Proteomics approach is an insightful tool for revealing the biological functions, metabolic networks and functional protein interaction networks of cells and organisms under stresses at the systematic level. The purpose of this review is to provide an insightful guideline on how to set up a proteomic investigation for revealing biomolecule mechanisms, protein biomarkers and metabolism networks related to stress response, pollutant recognition, transport and biodegradation, and providing an insightful high-throughput approach for screening functional enzymes and effective microbes based on bioinformatics and functional verification method. Furthermore, the toxicity evaluation of pollutants and byproducts by proteomics approaches provides a scientific insight for early diagnosis of ecological risk and determination of the effectiveness of pollutant treatment techniques.

Genomic analysis of Asian honeybee populations in China reveals evolutionary relationships and adaptation to abiotic stress

The geographic and biological diversity of China has resulted in the differential adaptation of the eastern honeybee, Apis cerana, to these varied habitats. A. cerana were collected from 14 locations in China. Their genomes were sequenced, and nucleotide polymorphisms were identified at more than 9 million sites. Both STRUCTURE and principal component analysis placed the bees into seven groups. Phylogenomic analysis groups the honeybees into many of the same clusters with high bootstrap values (91%-100%). Populations from Tibet and South Yunnan are sister taxa and together represent the earliest diverging lineage included in this study. We propose that the evolutionary origin of A. cerana in China was in the southern region of Yunnan Province and expanded from there into the southeastern regions and into the northeastern mountain regions. The Cold-Temperate West Sichuan Plateau and Tropical Diannan populations were compared to identify genes under adaptive selection in these two habitats. Pathway enrichment analysis showing genes under selection, including the Hippo signaling pathway, GABAergic pathway, and trehalose-phosphate synthase, indicates that most genes under selection pressure are involved in the process of signal transduction and energy metabolism. qRT-PCR analysis reveals that one gene under selection, the AcVIAAT gene, involved in the GABAergic pathway, is responding to cold temperature stress. Through homologous recombination, we show that the AcVIAAT gene is able to replace the CNAG_01904 gene in the fungus Cryptococcus neoformans and that it makes the fungus less sensitive to conditions of oxidative stress and variations in temperature. Our results contribute to our understanding of the evolutionary origin of A. cerana in China and the molecular basis of environmental adaptation.

Predicting the optimal growth temperatures of prokaryotes using only genome derived features

MOTIVATION

Optimal growth temperature is a fundamental characteristic of all living organisms. Knowledge of this temperature is central to the study of a prokaryote, the thermal stability and temperature dependent activity of its genes, and the bioprospecting of its genome for thermally adapted proteins. While high throughput sequencing methods have dramatically increased the availability of genomic information, the growth temperatures of the source organisms are often unknown. This limits the study and technological application of these species and their genomes. Here, we present a novel method for the prediction of growth temperatures of prokaryotes using only genomic sequences.

RESULTS

By applying the reverse ecology principle that an organism's genome includes identifiable adaptations to its native environment, we can predict a species' optimal growth temperature with an accuracy of 5.17°C root-mean-square error and a coefficient of determination of 0.835. The accuracy can be further improved for specific taxonomic clades or by excluding psychrophiles. This method provides a valuable tool for the rapid calculation of organism growth temperature when only the genome sequence is known.

AVAILABILITY AND IMPLEMENTATION

Source code, genomes analyzed and features calculated are available at: https://github.com/DavidBSauer/OGT_prediction.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome Sequencing of Mesonia algae K4-1 Reveals Its Adaptation to the Arctic Ocean

The special ecological environment of the Arctic has brought about a large number of salt-tolerant and psychrotolerant microorganisms. We isolated two culturable bacterial strains of the genus Mesonia; one from the Arctic ocean, Mesonia algae K4-1, and one from the tropical sea, Mesonia sp. HuA40. Our genome analysis and phenotypic experiments indicated that Mesonia algae K4-1 is a moderately halophilic and psychrophilic bacterium. Mesonia algae K4-1 can tolerate 3–14% NaCl and grow at a wide range of temperatures from 4 to 50°C. Mesonia sp. HuA40 is a mesophilic bacterium that can only grow with 3–9% NaCl. In addition, the salt adaptation strategy of Mesonia algae K4-1 accumulates organic osmolytes in the cell. RNA helicases, glutathione and organic compatible solutes may play important roles in maintaining the metabolism and physiological function of Mesonia algae K4-1 under cold stress. Moreover, the ability of Mesonia algae K4-1 to adapt to an oligotrophic marine environment is likely due to the synthesis of a large number of extracellular polysaccharides and the secretion of various families of extracellular proteases. This study systematically analyzed the relationship between genomic differentiation and environmental factors of the Mesonia genus and revealed the possible adaptation mechanism of Mesonia algae K4-1 in the extreme Arctic marine environment at the genomic level.

The complete mitochondrial genome of the largest amphipod, Alicella gigantea: Insight into its phylogenetic relationships and deep sea adaptive characters

Alicella gigantea (Alicelloidae) is a scavenger with the largest body size among amphipods. It is a participant in the foodweb of deepsea ecosystem and distributed with vast bathymetric and geographic ranges. In this study, the mitochondrial genome of A. gigantea was completely assembled and characterized. The complete sequence has a total length of 16,851 bp, comprising the usual eukaryotic components, with 13 protein-coding genes (PCGs), 2 ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and 2 noncoding control regions (CRs). The gene rearrangement and reverse nucleotide strand bias of its mitochondrial genome are similar to those observed in the deepsea amphipod Eurythenes maldoror (Eurytheneidae), but different from the characters of Halice sp. MT-2017 (Dexaminoidea), an inhabitant of a deeper environment. Phylogenetic analysis indicates that A. gigantea occupies the basal branch of deepsea species-E. maldoror and Hirondellea gigas. This phylogeny supports the hypothesis that the evolution of hadal amphipods has undergone a transition from the abyssal depth. Compared to 41 available shallow water equivalents, the four accessible mitochondrial genomes from the deep sea, including the one produced in this study, show significantly fewer charged amino acids in the 13 PCGs, which suggests an adaption to the deepsea environment.

Characterization of the mitochondrial genome of an ancient amphipod Halice sp. MT-2017 (Pardaliscidae) from 10,908 m in the Mariana Trench

Small amphipods (Halice sp. MT-2017) with body length <1 cm were collected from the Challenger Deep (~10,920 m below sea level). The divergence time of their lineage was approximately 109 Mya, making this group ancient compared to others under study. The mitochondrial genome of Halice sp. shared the usual gene components of metazoans, comprising 13 protein coding genes (PCGs), 22 transfer RNAs (tRNAs), and 2 ribosomal RNAs (rRNAs). The arrangement of these genes, however, differed greatly from that of other amphipods. Of the 15 genes that were rearranged with respect to the pancrustacean gene pattern, 12 genes (2 PCGs, 2 rRNAs, and 8 tRNAs) were both translocated and strand-reversed. In contrast, the mitochondrial genomes in other amphipods never show so many reordered genes, and in most instances, only tRNAs were involved in strand-reversion-coupled translocation. Other characteristics, including reversed strand nucleotide composition bias, relatively higher composition of non-polar amino acids, and lower evolutionary rate, were also identified. Interestingly, the latter two features were shared with another hadal amphipod, Hirondellea gigas, suggesting their possible associations with the adaptation to deep-sea extreme habitats. Overall, our data provided a useful resource for future studies on the evolutionary and adaptive mechanisms of hadal faunas.

Adaptation of Proteins to the Cold in Antarctic Fish: A Role for Methionine?

The evolution of antifreeze glycoproteins has enabled notothenioid fish to flourish in the freezing waters of the Southern Ocean. Whereas successful at the biodiversity level to life in the cold, paradoxically at the cellular level these stenothermal animals have problems producing, folding, and degrading proteins at their ambient temperatures of -1.86 °C. In this first multi-species transcriptome comparison of the amino acid composition of notothenioid proteins with temperate teleost proteins, we show that, unlike psychrophilic bacteria, Antarctic fish provide little evidence for the mass alteration of protein amino acid composition to enhance protein folding and reduce protein denaturation in the cold. The exception was the significant overrepresentation of positions where leucine in temperate fish proteins was replaced by methionine in the notothenioid orthologues. We hypothesize that these extra methionines have been preferentially assimilated into the genome to act as redox sensors in the highly oxygenated waters of the Southern Ocean. This redox hypothesis is supported by analyses of notothenioids showing enrichment of genes associated with responses to environmental stress, particularly reactive oxygen species. So overall, although notothenioid fish show cold-associated problems with protein homeostasis, they may have modified only a selected number of biochemical pathways to work efficiently below 0 °C. Even a slight warming of the Southern Ocean might disrupt the critical functions of this handful of key pathways with considerable impacts for the functioning of this ecosystem in the future.

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.

Roles of Active-Site Aromatic Residues in Cold Adaptation of Sphingomonas glacialis Esterase EstSP1

The aromatic amino acids, Tyr or Trp, which line the active-site walls of esterases, stabilize the catalytic His loop via hydrogen bonding. A Tyr residue is preferred in extremophilic esterases (psychrophilic or hyperthermophilic esterases), whereas a Trp residue is preferred in moderate-temperature esterases. Here, we provide evidence that Tyr and Trp play distinct roles in cold adaptation of the psychrophilic esterase EstSP1 isolated from an Arctic bacterium Sphingomonas glacialis PAMC 26605. Stern-Volmer plots showed that the mutation of Tyr191 to Ala, Phe, Trp, and His resulted in reduced conformational flexibility of the overall protein structure. Interestingly, the Y191W and Y191H mutants showed increased thermal stability at moderate temperatures. All Tyr191 mutants showed reduced catalytic activity relative to wild-type EstSP1. Our results indicate that Tyr with its phenyl hydroxyl group is favored for increased conformational flexibility and high catalytic activity of EstSP1 at low temperatures at the expense of thermal stability. The results of this study suggest that, in the permanently cold Arctic zone, enzyme activity has been selected for psychrophilic enzymes over thermal stability. The results presented herein provide novel insight into the roles of Tyr and Trp residues for temperature adaptation of enzymes that function at low, moderate, and high temperatures.

Cold adaptation of tRNA nucleotidyltransferases: A tradeoff in activity, stability and fidelity

Cold adaptation is an evolutionary process that has dramatic impact on enzymatic activity. Increased flexibility of the protein structure represents the main evolutionary strategy for efficient catalysis and reaction rates in the cold, but is achieved at the expense of structural stability. This results in a significant activity-stability tradeoff, as it was observed for several metabolic enzymes. In polymerases, however, not only reaction rates, but also fidelity plays an important role, as these enzymes have to synthesize copies of DNA and RNA as exact as possible. Here, we investigate the effects of cold adaptation on the highly accurate CCA-adding enzyme, an RNA polymerase that uses an internal amino acid motif within the flexible catalytic core as a template to synthesize the CCA triplet at tRNA 3'-ends. As the relative orientation of these residues determines nucleotide selection, we characterized how cold adaptation impacts template reading and fidelity. In a comparative analysis of closely related psychro-, meso-, and thermophilic enzymes, the cold-adapted polymerase shows a remarkable error rate during CCA synthesis in vitro as well as in vivo. Accordingly, CCA-adding activity at low temperatures is not only achieved at the expense of structural stability, but also results in a reduced polymerization fidelity.

Adaptation and evolution of deep-sea scale worms (Annelida: Polynoidae): insights from transcriptome comparison with a shallow-water species

Polynoid scale worms (Polynoidae, Annelida) invaded deep-sea chemosynthesis-based ecosystems approximately 60 million years ago, but little is known about their genetic adaptation to the extreme deep-sea environment. In this study, we reported the first two transcriptomes of deep-sea polynoids (Branchipolynoe pettiboneae, Lepidonotopodium sp.) and compared them with the transcriptome of a shallow-water polynoid (Harmothoe imbricata). We determined codon and amino acid usage, positive selected genes, highly expressed genes and putative duplicated genes. Transcriptome assembly produced 98,806 to 225,709 contigs in the three species. There were more positively charged amino acids (i.e., histidine and arginine) and less negatively charged amino acids (i.e., aspartic acid and glutamic acid) in the deep-sea species. There were 120 genes showing clear evidence of positive selection. Among the 10% most highly expressed genes, there were more hemoglobin genes with high expression levels in both deep-sea species. The duplicated genes related to DNA recombination and metabolism, and gene expression were only enriched in deep-sea species. Deep-sea scale worms adopted two strategies of adaptation to hypoxia in the chemosynthesis-based habitats (i.e., rapid evolution of tetra-domain hemoglobin in Branchipolynoe or high expression of single-domain hemoglobin in Lepidonotopodium sp.).