Between genotype and phenotype: protein chaperones and evolvability

A new lysine biosynthetic enzyme from a bacterial endosymbiont shaped by genetic drift and genome reduction

The effect of population bottlenecks and genome reduction on enzyme function is poorly understood. Candidatus Liberibacter solanacearum is a bacterium with a reduced genome that is transmitted vertically to the egg of an infected psyllid-a population bottleneck that imposes genetic drift and is predicted to affect protein structure and function. Here, we define the function of Ca. L. solanacearum dihydrodipicolinate synthase (CLsoDHDPS), which catalyzes the committed branchpoint reaction in diaminopimelate and lysine biosynthesis. We demonstrate that CLsoDHDPS is expressed in Ca. L. solanacearum and expression is increased ~2-fold in the insect host compared to in planta. CLsoDHDPS has decreased thermal stability and increased aggregation propensity, implying mutations have destabilized the enzyme but are compensated for through elevated chaperone expression and a stabilized oligomeric state. CLsoDHDPS uses a ternary-complex kinetic mechanism, which is to date unique among DHDPS enzymes, has unusually low catalytic ability, but an unusually high substrate affinity. Structural studies demonstrate that the active site is more open, and the structure of CLsoDHDPS with both pyruvate and the substrate analogue succinic-semialdehyde reveals that the product is both structurally and energetically different and therefore evolution has in this case fashioned a new enzyme. Our study suggests the effects of genome reduction and genetic drift on the function of essential enzymes and provides insights on bacteria-host co-evolutionary associations. We propose that bacteria with endosymbiotic lifestyles present a rich vein of interesting enzymes useful for understanding enzyme function and/or informing protein engineering efforts.

Mutational robustness and the role of buffer genes in evolvability

Organisms rely on mutations to fuel adaptive evolution. However, many mutations impose a negative effect on fitness. Cells may have therefore evolved mechanisms that affect the phenotypic effects of mutations, thus conferring mutational robustness. Specifically, so-called buffer genes are hypothesized to interact directly or indirectly with genetic variation and reduce its effect on fitness. Environmental or genetic perturbations can change the interaction between buffer genes and genetic variation, thereby unmasking the genetic variation's phenotypic effects and thus providing a source of variation for natural selection to act on. This review provides an overview of our understanding of mutational robustness and buffer genes, with the chaperone gene HSP90 as a key example. It discusses whether buffer genes merely affect standing variation or also interact with de novo mutations, how mutational robustness could influence evolution, and whether mutational robustness might be an evolved trait or rather a mere side-effect of complex genetic interactions.

Unraveling stress resilience: Insights from adaptations to extreme environments by Astyanax mexicanus cavefish

Extreme environmental conditions have profound impacts on shaping the evolutionary trajectory of organisms. Exposure to these conditions elicits stress responses, that can trigger phenotypic changes in novel directions. The Mexican Tetra, Astyanax mexicanus, is an excellent model for understanding evolutionary mechanisms in response to extreme or new environments. This fish species consists of two morphs; the classical surface-dwelling fish and the blind cave-dwellers that inhabit dark and biodiversity-reduced ecosystems. In this review, we explore the specific stressors present in cave environments and examine the diverse adaptive strategies employed by cave populations to not only survive but thrive as successful colonizers. By analyzing the evolutionary responses of A. mexicanus, we gain valuable insights into the genetic, physiological, and behavioral adaptations that enable organisms to flourish under challenging environmental conditions.

Proteostasis governs differential temperature sensitivity across embryonic cell types

Embryonic development is remarkably robust, but temperature stress can degrade its ability to generate animals with invariant anatomy. Phenotypes associated with environmental stress suggest that some cell types are more sensitive to stress than others, but the basis of this sensitivity is unknown. Here, we characterize hundreds of individual zebrafish embryos under temperature stress using whole-animal single-cell RNA sequencing (RNA-seq) to identify cell types and molecular programs driving phenotypic variability. We find that temperature perturbs the normal proportions and gene expression programs of numerous cell types and also introduces asynchrony in developmental timing. The notochord is particularly sensitive to temperature, which we map to a specialized cell type: sheath cells. These cells accumulate misfolded protein at elevated temperature, leading to a cascading structural failure of the notochord and anatomic defects. Our study demonstrates that whole-animal single-cell RNA-seq can identify mechanisms for developmental robustness and pinpoint cell types that constitute key failure points.

Intraspecific variation among Chinook Salmon populations indicates physiological adaptation to local environmental conditions

Understanding interpopulation variation is important to predicting species responses to climate change. Recent research has revealed interpopulation variation among several species of Pacific salmonids; however, the environmental drivers of population differences remain elusive. We tested for local adaptation and countergradient variation by assessing interpopulation variation among six populations of fall-run Chinook Salmon from the western United States. Juvenile fish were reared at three temperatures (11, 16 and 20°C), and five physiological metrics were measured (routine and maximum metabolic rate, aerobic scope, growth rate and critical thermal maximum). We then tested associations between these physiological metrics and 15 environmental characteristics (e.g. rearing temperature, latitude, migration distance, etc.). Statistical associations between the five physiological metrics and 15 environmental characteristics supported our hypotheses of local adaptation. Notably, latitude was a poor predictor of population physiology. Instead, our results demonstrate that populations from warmer habitats exhibit higher thermal tolerance (i.e. critical thermal maxima), faster growth when warm acclimated and greater aerobic capacity at high temperatures. Additionally, populations with longer migrations exhibit higher metabolic capacity. However, overall metabolic capacity declined with warm acclimation, indicating that future climate change may reduce metabolic capacity, negatively affecting long-migrating populations. Linking physiological traits to environmental characteristics enables flexible, population-specific management of disparate populations in response to local conditions.

Population-based heteropolymer design to mimic protein mixtures

Biological fluids, the most complex blends, have compositions that constantly vary and cannot be molecularly defined1. Despite these uncertainties, proteins fluctuate, fold, function and evolve as programmed2-4. We propose that in addition to the known monomeric sequence requirements, protein sequences encode multi-pair interactions at the segmental level to navigate random encounters5,6; synthetic heteropolymers capable of emulating such interactions can replicate how proteins behave in biological fluids individually and collectively. Here, we extracted the chemical characteristics and sequential arrangement along a protein chain at the segmental level from natural protein libraries and used the information to design heteropolymer ensembles as mixtures of disordered, partially folded and folded proteins. For each heteropolymer ensemble, the level of segmental similarity to that of natural proteins determines its ability to replicate many functions of biological fluids including assisting protein folding during translation, preserving the viability of fetal bovine serum without refrigeration, enhancing the thermal stability of proteins and behaving like synthetic cytosol under biologically relevant conditions. Molecular studies further translated protein sequence information at the segmental level into intermolecular interactions with a defined range, degree of diversity and temporal and spatial availability. This framework provides valuable guiding principles to synthetically realize protein properties, engineer bio/abiotic hybrid materials and, ultimately, realize matter-to-life transformations.

Upwelling-level acidification and pH/pCO2 variability moderate effects of ocean acidification on brain gene expression in the temperate surfperch, Embiotoca jacksoni

Acidification-induced changes in neurological function have been documented in several tropical marine fishes. Here, we investigate whether similar patterns of neurological impacts are observed in a temperate Pacific fish that naturally experiences regular and often large shifts in environmental pH/pCO₂ . In two laboratory experiments, we tested the effect of acidification, as well as pH/pCO₂ variability, on gene expression in the brain tissue of a common temperate kelp forest/estuarine fish, Embiotoca jacksoni. Experiment 1 employed static pH treatments (target pH = 7.85/7.30), while Experiment 2 incorporated two variable treatments that oscillated around corresponding static treatments with the same mean (target pH = 7.85/7.70) in an eight-day cycle (amplitude ± 0.15). We found that patterns of global gene expression differed across pH level treatments. Additionally, we identified differential expression of specific genes and enrichment of specific gene sets (GSEA) in comparisons of static pH treatments and in comparisons of static and variable pH treatments of the same mean pH. Importantly, we found that pH/pCO₂ variability decreased the number of differentially expressed genes detected between high and low pH treatments, and that inter-individual variability in gene expression was greater in variable treatments than static treatments. These results provide important confirmation of neurological impacts of acidification in a temperate fish species and, critically, that natural environmental variability may mediate the impacts of ocean acidification.

Reduced physiological plasticity in a fish adapted to stable temperatures

Significance Plastic individuals can buffer environmental changes, maintaining a stable performance across gradients. Plasticity is therefore thought to be particularly beneficial for the survival of wild populations that experience large environmental fluctuations, such as diel and seasonal temperature changes. Maintaining plasticity is widely assumed to be costly; however, empirical evidence demonstrating this cost is scarce. Here, we predict that if plasticity is costly, it would be readily lost in a stable environment, such as a laboratory. To test this, we measured a diverse range of phenotypic traits, spanning gene expression, physiology, and behavior, in wild and laboratory zebrafish acclimated to 15 temperatures. We show that laboratory fish have lost plasticity in many traits, demonstrating that maintaining plasticity carries a cost.

In vitro evolution of uracil glycosylase towards DnaKJ and GroEL binding evolves different misfolded states

Natural evolution is driven by random mutations that improve fitness. In vitro evolution mimics this process, however, on a short time-scale and is driven by the given bait. Here, we used directed in vitro evolution of a random mutant library of Uracil glycosylase (eUNG) displayed on yeast surface to select for binding to chaperones GroEL, DnaK+DnaJ+ATP (DnaKJ) or E.coli cell extract (CE), using binding to the eUNG inhibitor Ugi as probe for native fold. The CE selected population was further divided to Ugi binders (+U) or non-binders (-U). The aim here was to evaluate the sequence space and physical state of the evolved protein binding the different baits. We found that GroEL, DnaKJ and CE-U select and enrich for mutations causing eUNG to misfold, with the three being enriched in mutations in buried and conserved positions, with a tendency to increase positive charge. Still, each selection had its own trajectory, with GroEL and CE-U selecting mutants highly sensitive to protease cleavage while DnaKJ selected partially structured misfolded species with a tendency to refold, making them less sensitive to proteases. More general, our results show that GroEL has a higher tendency to purge promiscuous misfolded protein mutants from the system, while DnaKJ binds misfolding-prone mutant species that are, upon chaperone release, more likely to natively refold. CE-U shares some of the properties of GroEL- and DnaKJ-selected populations, while harboring also unique properties that can be explained by the presence of additional chaperones in CE, such as Trigger factor, HtpG and ClpB.

Protein folding in vitro and in the cell: From a solitary journey to a team effort

Correct protein folding is essential for the health and function of living organisms. Yet, it is not well understood how unfolded proteins reach their native state and avoid aggregation, especially within the cellular milieu. Some proteins, especially small, single-domain and apparent two-state folders, successfully attain their native state upon dilution from denaturant. Yet, many more proteins undergo misfolding and aggregation during this process, in a concentration-dependent fashion. Once formed, native and aggregated states are often kinetically trapped relative to each other. Hence, the early stages of protein life are absolutely critical for proper kinetic channeling to the folded state and for long-term solubility and function. This review summarizes current knowledge on protein folding/aggregation mechanisms in buffered solution and within the bacterial cell, highlighting early stages. Remarkably, teamwork between nascent chain, ribosome, trigger factor and Hsp70 molecular chaperones enables all proteins to overcome aggregation propensities and reach a long-lived bioactive state.

Non-adaptive complexity and biochemical function

Intricate biochemical structures are usually thought to be useful, because natural selection preserves them from degradation by a constant hail of destructive mutations. Biochemists therefore often deliberately disrupt them to understand how complexity improves protein function or fitness. However, evolutionary theory suggests that even useless complexity that never improved fitness can become completely essential if a simple set of evolutionary conditions is fulfilled. We review evidence that stable protein complexes, protein-chaperone interactions, and complexes consisting of several paralogs all fulfill these conditions. This makes reverse genetics or destructive mutagenesis unsuitable for assigning functions to these kinds of complexity. Instead, we advocate that incorporating evolutionary approaches into biochemistry overcomes this difficulty and allows us to distinguish useless from useful biochemical complexity.

Morphological canalization, integration, and plasticity in response to population density in Abutilon theophrasti: Influences of soil conditions and growth stages

Phenotypic integration and developmental canalization have been hypothesized to constrain the degree of phenotypic plasticity, but little evidence exists, probably due to the lack of studies on the relationships among the three processes, especially for plants under different environments. We conducted a field experiment by subjecting plants of Abutilon theophrasti to three densities, under infertile and fertile soil conditions, and analyzing correlations among canalization, integration, and plasticity in a variety of measured morphological traits after 50 and 70 days, to investigate the relationships among the three variables in response to density and how these responses vary with soil conditions and growth stages. Results showed trait canalization decreased and phenotypic integration and the degree of plasticity (absolute plasticity) in traits increased with density. Phenotypic integration often positively correlated with absolute plasticity, whereas correlations between trait canalization and plasticity were insignificant in most cases, with a few positive ones between canalization and absolute plasticity at low and medium densities. As plants grew, these correlations intensified in infertile soil and attenuated in fertile soil. Our findings suggested the complexity of the relationship between canalization and plasticity: Decreased canalization is more likely to facilitate active plastic responses under more favorable conditions, whereas increased level of integration should mainly be an outcome of plastic responses. Soil conditions and growth stage may affect responses of these correlations to density via modifying plant size, competition strength, and plastic responses in traits. We also predicted that decreased canalization can be advantageous or disadvantageous, and the lack of response to stress may demonstrate a stronger ability of adaptation than passive response, thus should be adaptive plasticity as active response.

Variables Influencing Differences in Sequence Conservation in the Fission Yeast Schizosaccharomyces pombe

Which variables determine the constraints on gene sequence evolution is one of the most central questions in molecular evolution. In the fission yeast Schizosaccharomyces pombe, an important model organism, the variables influencing the rate of sequence evolution have yet to be determined. Previous studies in other single celled organisms have generally found gene expression levels to be most significant, with numerous other variables such as gene length and functional importance identified as having a smaller impact. Using publicly available data, we used partial least squares regression, principal components regression, and partial correlations to determine the variables most strongly associated with sequence evolution constraints. We identify centrality in the protein–protein interactions network, amino acid composition, and cellular location as the most important determinants of sequence conservation. However, each factor only explains a small amount of variance, and there are numerous variables having a significant or heterogeneous influence. Our models explain more than half of the variance in dN, raising the possibility that future refined models could quantify the role of stochastics in evolutionary rate variation.

The Adaptive Potential of the Middle Domain of Yeast Hsp90

The distribution of fitness effects (DFEs) of new mutations across different environments quantifies the potential for adaptation in a given environment and its cost in others. So far, results regarding the cost of adaptation across environments have been mixed, and most studies have sampled random mutations across different genes. Here, we quantify systematically how costs of adaptation vary along a large stretch of protein sequence by studying the distribution of fitness effects of the same ≈2,300 amino-acid changing mutations obtained from deep mutational scanning of 119 amino acids in the middle domain of the heat shock protein Hsp90 in five environments. This region is known to be important for client binding, stabilization of the Hsp90 dimer, stabilization of the N-terminal-Middle and Middle-C-terminal interdomains, and regulation of ATPase–chaperone activity. Interestingly, we find that fitness correlates well across diverse stressful environments, with the exception of one environment, diamide. Consistent with this result, we find little cost of adaptation; on average only one in seven beneficial mutations is deleterious in another environment. We identify a hotspot of beneficial mutations in a region of the protein that is located within an allosteric center. The identified protein regions that are enriched in beneficial, deleterious, and costly mutations coincide with residues that are involved in the stabilization of Hsp90 interdomains and stabilization of client-binding interfaces, or residues that are involved in ATPase–chaperone activity of Hsp90. Thus, our study yields information regarding the role and adaptive potential of a protein sequence that complements and extends known structural information.

Host's guardian protein counters degenerative symbiont evolution

When laying eggs, plataspid stinkbugs deposit small packets called “symbiont capsules.” Newborn stinkbugs suck the capsules to acquire a bacterial mutualist. Without the symbiont, the babies cannot grow and die. Due to the long-lasting symbiosis, the symbiont has experienced genome reduction and become uncultivable. Within the capsules, however, the symbiont can survive for over a week outside the host. Why and how? Here, we uncover a molecular secret implemented in the symbiont capsules. Mother stinkbugs massively produce a special intestinal secretion protein, PMDP, at the expense of their own survival. PMDP embeds the fragile symbiont cells and protects them within the capsules. The host-provisioned molecule for sustaining symbiosis may be utilized for cultivation and/or preservation of fastidious microorganisms.

Microbial symbioses significantly contribute to diverse organisms, where long-lasting associations tend to result in symbiont genome erosion, uncultivability, extinction, and replacement. How such inherently deteriorating symbiosis can be harnessed to stable partnership is of general evolutionary interest. Here, we report the discovery of a host protein essential for sustaining symbiosis. Plataspid stinkbugs obligatorily host an uncultivable and genome-reduced gut symbiont, Ishikawaella. Upon oviposition, females deposit “capsules” for symbiont delivery to offspring. Within the capsules, the fragile symbiotic bacteria survive the harsh conditions outside the host until acquired by newborn nymphs to establish vertical transmission. We identified a single protein dominating the capsule content, which is massively secreted by female-specific intestinal organs, embedding the symbiont cells, and packaged into the capsules. Knockdown of the protein resulted in symbiont degeneration, arrested capsule production, symbiont transmission failure, and retarded nymphal growth, unveiling its essential function for ensuring symbiont survival and vertical transmission. The protein originated from a lineage of odorant-binding protein-like multigene family, shedding light on the origin of evolutionary novelty regarding symbiosis. Experimental suppression of capsule production extended the female’s lifespan, uncovering a substantial cost for maintaining symbiosis. In addition to the host’s guardian protein, the symbiont’s molecular chaperone, GroEL, was overproduced in the capsules, highlighting that the symbiont’s eroding functionality is compensated for by stabilizer molecules of host and symbiont origins. Our finding provides insight into how intimate host–symbiont associations can be maintained over evolutionary time despite the symbiont’s potential vulnerability to degeneration and malfunctioning.

Connexin 43 contributes to phenotypic robustness of the mouse skull

BACKGROUND

We compared skull shape and variation among genetically modified mice that exhibit different levels of connexin43 (Cx43) channel function, to determine whether Cx43 contributes to craniofacial phenotypic robustness. Specifically, we used two heterozygous mutant mouse models (G60S/+ and I130T/+) that, when compared to their wildtype counterparts, have an ~80% and ~50% reduction in Cx43 function, respectively.

RESULTS

Both mutant strains showed significant differences in skull shape compared to wildtype littermates and while these differences were more severe in the G60S/+ mouse, shape differences were localized to similar regions of the skull in both mutants. However, increased skull shape variation was observed in G60S/+ mutants only. Additionally, covariation of skull structures was disrupted in the G60S/+ mutants only, indicating that while a 50% reduction in Cx43 function is sufficient to cause a shift in mean skull shape, the threshold for Cx43 function for disrupting craniofacial phenotypic robustness is lower.

CONCLUSIONS

Collectively, our results indicate Cx43 can contribute to phenotypic robustness of the skull through a nonlinear relationship between Cx43 gap junctional function and phenotypic outcomes.

Linking hsp90's role as an evolutionary capacitator to the development of cancer

This paper links Heat Shock Protein Hsp90 as an evolutionary capacitator with the development of cancer. Hsp90 stabilises proteins associated with cancer in a number of ways. Canalisation allows for the accumulation of malignant mutations in the genome, and selection of beneficial phenotypes when cancer cells are stressed, allowing oncogenic development and progression. Hsp90 may allow for mutational 'big bangs' that can trigger primary malignant transformation. Hsp90 buffers catastrophic mutations in the oncogenome to prevent protein degradation and cellular apoptosis. Hsp90 was found to prevent the degradation of mutated p53, encouraging uncontrolled proliferation of cancer cells. Hsp90 buffering of mutations in response to cytotoxic therapy can lead to expression of beneficial phenotypes when Hsp90 is supressed and development of drug resistance. Trials with Hsp90 inhibitors have shown some success as an adjunctive therapy in preventing cancer progression, development of drug resistance, and even re-sensitisation to therapy after chemoresistance has developed.

Evaluating the correlation between genome-wide diversity and the release of plastic phenotypic variation in experimental translocations to novel natural environments

Phenotypic reaction norms are often shaped and constrained by selection and are important for allowing organisms to respond to environmental change. However, selection cannot constrain reaction norms for environmental conditions that populations have not experienced. Consequently, cryptic neutral genetic variation for the reaction norm can accumulate such that a release of phenotypic variation occurs upon exposure to novel14 conditions. Most genomic diversity behaves as if functionally neutral. Therefore, genome-wide diversity metrics may correlate with levels of cryptic genetic variation and, as a result, exhibit a positive relationship with a release of phenotypic variation in novel environments. To test this hypothesis, we conducted translocations of juvenile brook trout (Salvelinus fontinalis) from 12 populations to novel uninhabited ponds that represented a gradient of environmental conditions. We assessed reaction norms for morphological traits (body size and four morphometric relative warps) across pond environmental gradients and evaluated the effect of genome-wide heterozygosity on phenotypic variability. All traits displayed plastic reaction norms. Overall, we found some evidence that a release of phenotypic variation consistent with cryptic genetic variation can occur in novel environmental conditions. However, the extent to which this release correlated with average genome-wide diversity was limited to only one of five traits examined. Our results suggest a limited link between genomic diversity26 and the accumulation of cryptic genetic variation in reaction norms. Similarly, reaction norms were constrained for many of the morphological traits examined. Past conditions may have constrained reaction norms in the putatively novel environments despite significant deviations from contemporary source population habitat. Additionally, as a generalist colonizing species brook trout may exhibit plastic phenotypes across a wide range of environmental conditions.

How to Build a Biological Machine Using Engineering Materials and Methods

We present work in 3D printing electric motors from basic materials as the key to building a self-replicating machine to colonise the Moon. First, we explore the nature of the biological realm to ascertain its essence, particularly in relation to the origin of life when the inanimate became animate. We take an expansive view of this to ascertain parallels between the biological and the manufactured worlds. Life must have emerged from the available raw material on Earth and, similarly, a self-replicating machine must exploit and leverage the available resources on the Moon. We then examine these lessons to explore the construction of a self-replicating machine using a universal constructor. It is through the universal constructor that the actuator emerges as critical. We propose that 3D printing constitutes an analogue of the biological ribosome and that 3D printing may constitute a universal construction mechanism. Following a description of our progress in 3D printing motors, we suggest that this engineering effort can inform biology, that motors are a key facet of living organisms and illustrate the importance of motors in biology viewed from the perspective of engineering (in the Feynman spirit of “what I cannot create, I cannot understand”).

Highly parallel lab evolution reveals that epistasis can curb the evolution of antibiotic resistance

Genetic perturbations that affect bacterial resistance to antibiotics have been characterized genome-wide, but how do such perturbations interact with subsequent evolutionary adaptation to the drug? Here, we show that strong epistasis between resistance mutations and systematically identified genes can be exploited to control spontaneous resistance evolution. We evolved hundreds of Escherichia coli K-12 mutant populations in parallel, using a robotic platform that tightly controls population size and selection pressure. We find a global diminishing-returns epistasis pattern: strains that are initially more sensitive generally undergo larger resistance gains. However, some gene deletion strains deviate from this general trend and curtail the evolvability of resistance, including deletions of genes for membrane transport, LPS biosynthesis, and chaperones. Deletions of efflux pump genes force evolution on inferior mutational paths, not explored in the wild type, and some of these essentially block resistance evolution. This effect is due to strong negative epistasis with resistance mutations. The identified genes and cellular functions provide potential targets for development of adjuvants that may block spontaneous resistance evolution when combined with antibiotics.

The antibiotic resistance crisis calls for new ways of restricting the ability of bacteria to evolve resistance. Here, Lukačišinová et al. perform highly controlled evolution experiments in E. coli strains to identify genetic perturbations that strongly limit the evolution of antibiotic resistance through epistasis.

Exploration of space to achieve scientific breakthroughs

Living organisms adapt to changing environments using their amazing flexibility to remodel themselves by a process called evolution. Environmental stress causes selective pressure and is associated with genetic and phenotypic shifts for better modifications, maintenance, and functioning of organismal systems. The natural evolution process can be used in complement to rational strain engineering for the development of desired traits or phenotypes as well as for the production of novel biomaterials through the imposition of one or more selective pressures. Space provides a unique environment of stressors (e.g., weightlessness and high radiation) that organisms have never experienced on Earth. Cells in the outer space reorganize and develop or activate a range of molecular responses that lead to changes in cellular properties. Exposure of cells to the outer space will lead to the development of novel variants more efficiently than on Earth. For instance, natural crop varieties can be generated with higher nutrition value, yield, and improved features, such as resistance against high and low temperatures, salt stress, and microbial and pest attacks. The review summarizes the literature on the parameters of outer space that affect the growth and behavior of cells and organisms as well as complex colloidal systems. We illustrate an understanding of gravity-related basic biological mechanisms and enlighten the possibility to explore the outer space environment for application-oriented aspects. This will stimulate biological research in the pursuit of innovative approaches for the future of agriculture and health on Earth.

Distinct metabolic states of a cell guide alternate fates of mutational buffering through altered proteostasis

Metabolic changes alter the cellular milieu; can this also change intracellular protein folding? Since proteostasis can modulate mutational buffering, if change in metabolism has the ability to change protein folding, arguably, it should also alter mutational buffering. Here we find that altered cellular metabolic states in E. coli buffer distinct mutations on model proteins. Buffered-mutants have folding problems in vivo and are differently chaperoned in different metabolic states. Notably, this assistance is dependent upon the metabolites and not on the increase in canonical chaperone machineries. Being able to reconstitute the folding assistance afforded by metabolites in vitro, we propose that changes in metabolite concentrations have the potential to alter protein folding capacity. Collectively, we unravel that the metabolite pools are bona fide members of proteostasis and aid in mutational buffering. Given the plasticity in cellular metabolism, we posit that metabolic alterations may play an important role in cellular proteostasis.

Changes in osmotic homeostasis alter metabolites and therefore chemical milieu of the cells. Here, the authors show that altering metabolites in E. coli also change the cellular capacity for buffering mutations that impair protein folding and influences proteostasis irrespective of molecular chaperones

Predicting the viability of beta-lactamase: How folding and binding free energies correlate with beta-lactamase fitness

One of the long-standing holy grails of molecular evolution has been the ability to predict an organism's fitness directly from its genotype. With such predictive abilities in hand, researchers would be able to more accurately forecast how organisms will evolve and how proteins with novel functions could be engineered, leading to revolutionary advances in medicine and biotechnology. In this work, we assemble the largest reported set of experimental TEM-1 β-lactamase folding free energies and use this data in conjunction with previously acquired fitness data and computational free energy predictions to determine how much of the fitness of β-lactamase can be directly predicted by thermodynamic folding and binding free energies. We focus upon β-lactamase because of its long history as a model enzyme and its central role in antibiotic resistance. Based upon a set of 21 β-lactamase single and double mutants expressly designed to influence protein folding, we first demonstrate that modeling software designed to compute folding free energies such as FoldX and PyRosetta can meaningfully, although not perfectly, predict the experimental folding free energies of single mutants. Interestingly, while these techniques also yield sensible double mutant free energies, we show that they do so for the wrong physical reasons. We then go on to assess how well both experimental and computational folding free energies explain single mutant fitness. We find that folding free energies account for, at most, 24% of the variance in β-lactamase fitness values according to linear models and, somewhat surprisingly, complementing folding free energies with computationally-predicted binding free energies of residues near the active site only increases the folding-only figure by a few percent. This strongly suggests that the majority of β-lactamase's fitness is controlled by factors other than free energies. Overall, our results shed a bright light on to what extent the community is justified in using thermodynamic measures to infer protein fitness as well as how applicable modern computational techniques for predicting free energies will be to the large data sets of multiply-mutated proteins forthcoming.

Collateral fitness effects of mutations

The distribution of fitness effects of mutation plays a central role in constraining protein evolution. The underlying mechanisms by which mutations lead to fitness effects are typically attributed to changes in protein specific activity or abundance. Here, we reveal the importance of a mutation's collateral fitness effects, which we define as effects that do not derive from changes in the protein's ability to perform its physiological function. We comprehensively measured the collateral fitness effects of missense mutations in the Escherichia coli TEM-1 β-lactamase antibiotic resistance gene using growth competition experiments in the absence of antibiotic. At least 42% of missense mutations in TEM-1 were deleterious, indicating that for some proteins collateral fitness effects occur as frequently as effects on protein activity and abundance. Deleterious mutations caused improper posttranslational processing, incorrect disulfide-bond formation, protein aggregation, changes in gene expression, and pleiotropic effects on cell phenotype. Deleterious collateral fitness effects occurred more frequently in TEM-1 than deleterious effects on antibiotic resistance in environments with low concentrations of the antibiotic. The surprising prevalence of deleterious collateral fitness effects suggests they may play a role in constraining protein evolution, particularly for highly expressed proteins, for proteins under intermittent selection for their physiological function, and for proteins whose contribution to fitness is buffered against deleterious effects on protein activity and protein abundance.

Chaperone client proteins evolve slower than non-client proteins

Chaperones are important molecular machinery that assists proteins to attain their native three-dimensional structure crucial for function. Earlier studies using experimental evolution showed that chaperones impose a relaxation of sequence constraints on their "client" proteins, which may lead to the fixation of slightly deleterious mutations on the latter. However, we hypothesized that such a phenomenon might be harmful to the organism in a natural physiological condition. In this study, we investigated the evolutionary rates of chaperone client and non-client proteins in five model organisms from both prokaryotic and eukaryotic lineages. Our study reveals a slower evolutionary rate of chaperone client proteins in all five organisms. Additionally, the slower folding rate and lower aggregation propensity of chaperone client proteins reveal that the chaperone may play an essential role in rescuing the slightly disadvantageous effects due to random mutations and subsequent protein misfolding. However, the fixation of such mutations is less likely to be selected in the natural population.

Molecular Chaperones Accelerate the Evolution of Their Protein Clients in Yeast

Protein stability is a major constraint on protein evolution. Molecular chaperones, also known as heat-shock proteins, can relax this constraint and promote protein evolution by diminishing the deleterious effect of mutations on protein stability and folding. This effect, however, has only been stablished for a few chaperones. Here, we use a comprehensive chaperone–protein interaction network to study the effect of all yeast chaperones on the evolution of their protein substrates, that is, their clients. In particular, we analyze how yeast chaperones affect the evolutionary rates of their clients at two very different evolutionary time scales. We first study the effect of chaperone-mediated folding on protein evolution over the evolutionary divergence of Saccharomyces cerevisiae and S. paradoxus. We then test whether yeast chaperones have left a similar signature on the patterns of standing genetic variation found in modern wild and domesticated strains of S. cerevisiae. We find that genes encoding chaperone clients have diverged faster than genes encoding non-client proteins when controlling for their number of protein–protein interactions. We also find that genes encoding client proteins have accumulated more intraspecific genetic diversity than those encoding non-client proteins. In a number of multivariate analyses, controlling by other well-known factors that affect protein evolution, we find that chaperone dependence explains the largest fraction of the observed variance in the rate of evolution at both evolutionary time scales. Chaperones affecting rates of protein evolution mostly belong to two major chaperone families: Hsp70s and Hsp90s. Our analyses show that protein chaperones, by virtue of their ability to buffer destabilizing mutations and their role in modulating protein genotype–phenotype maps, have a considerable accelerating effect on protein evolution.

Hypothesis: protein and RNA attributes are continuously optimized over time

Background

Little is known why proteins and RNAs exhibit half-lives varying over several magnitudes. Despite many efforts, a conclusive link between half-lives and gene function could not be established suggesting that other determinants may influence these molecular attributes.

Results

Here, I find that with increasing gene age there is a gradual and significant increase of protein and RNA half-lives, protein structure, and other molecular attributes that tend to affect protein abundance. These observations are accommodated in a hypothesis which posits that new genes at ‘birth’ are not optimized and thus their products exhibit low half-lives and less structure but continuous mutagenesis eventually improves these attributes. Thus, the protein and RNA products of the oldest genes obtained their high degrees of stability and structure only after billions of years while the products of younger genes had less time to be optimized and are therefore less stable and structured. Because more stable proteins with lower turnover require less transcription to maintain the same level of abundance, reduced transcription-associated mutagenesis (TAM) would fixate the changes by increasing gene conservation.

Conclusions

Consequently, the currently observed diversity of molecular attributes is a snapshot of gene products being at different stages along their temporal path of optimization.

Snails in the sun: Strategies of terrestrial gastropods to cope with hot and dry conditions

Terrestrial gastropods do not only inhabit humid and cool environments but also habitat in which hot and dry conditions prevail. Snail species that are able to cope with such climatic conditions are thus expected to having developed multifaceted strategies and mechanisms to ensure their survival and reproduction under heat and desiccation stress. This review paper aims to provide an integrative overview of the numerous adaptation strategies terrestrial snails have evolved to persist in hot and dry environments as well as their mutual interconnections and feedbacks, but also to outline research gaps and questions that remained unanswered. We extracted relevant information from more than 140 publications in order to show how biochemical, cellular, physiological, morphological, ecological, thermodynamic, and evolutionary parameters contribute to provide an overall picture of this classical example in stress ecology. These mechanisms range from behavioral and metabolic adaptations, including estivation, to the induction of chaperones and antioxidant enzymes, mucocyte and digestive gland cell responses and the modification and frequency of morphological features, particularly shell pigmentation. In this context, thermodynamic constraints call for processes of complex adaptation at varying levels of biological organization that are mutually interwoven. We were able to assemble extensive, mostly narrowly focused information from the literature into a web of network parameters, showing that future work on this subject requires multicausal thinking to account for the complexity of relationships involved in snails' adaptation to insolation, heat, and drought.

Candidate genes expression profiling during wilting in chickpea caused by Fusarium oxysporum f. sp. ciceris race 5

Chickpea production may be seriously threatened by Fusarium wilt, a disease caused by the soil-borne fungus Fusarium oxysporum f. sp. ciceris. F. oxysporum race 5 is the most important race in the Mediterranean basin. Recently, the region responsible for resistance race 5 has been delimited within a region on chromosome 2 that spans 820 kb. To gain a better understanding of this genomic region, we used a transcriptomic approach based on quantitative real-time PCR to analyze the expression profiles of 22 selected candidate genes. We used a pair of near-isogenic lines (NILs) differing in their sensitivity to Fusarium race 5 (resistant vs susceptible) to monitor the transcriptional changes over a time-course experiment (24, 48, and 72 hours post inoculation, hpi). Qualitative differences occurred during the timing of regulation. A cluster of 12 genes were induced by the resistant NIL at 24 hpi, whereas a second cluster contained 9 genes induced by the susceptible NIL at 48 hpi. Their possible functions in the molecular defence of chickpea is discussed. Our study provides new insight into the molecular defence against Fusarium race 5 and demonstrates that development of NILs is a rich resource to facilitate the detection of candidate genes. The new genes regulated here may be useful against other Fusarium races.

Filaments and phenotypes: cellular roles and orphan effects associated with mutations in cytoplasmic intermediate filament proteins

Cytoplasmic intermediate filaments (IFs) surround the nucleus and are often anchored at membrane sites to form effectively transcellular networks. Mutations in IF proteins (IFps) have revealed mechanical roles in epidermis, muscle, liver, and neurons. At the same time, there have been phenotypic surprises, illustrated by the ability to generate viable and fertile mice null for a number of IFp-encoding genes, including vimentin. Yet in humans, the vimentin ( VIM) gene displays a high probability of intolerance to loss-of-function mutations, indicating an essential role. A number of subtle and not so subtle IF-associated phenotypes have been identified, often linked to mechanical or metabolic stresses, some of which have been found to be ameliorated by the over-expression of molecular chaperones, suggesting that such phenotypes arise from what might be termed "orphan" effects as opposed to the absence of the IF network per se, an idea originally suggested by Toivola et al. and Pekny and Lane.

The causes of evolvability and their evolution

Evolvability is the ability of a biological system to produce phenotypic variation that is both heritable and adaptive. It has long been the subject of anecdotal observations and theoretical work. In recent years, however, the molecular causes of evolvability have been an increasing focus of experimental work. Here, we review recent experimental progress in areas as different as the evolution of drug resistance in cancer cells and the rewiring of transcriptional regulation circuits in vertebrates. This research reveals the importance of three major themes: multiple genetic and non-genetic mechanisms to generate phenotypic diversity, robustness in genetic systems, and adaptive landscape topography. We also discuss the mounting evidence that evolvability can evolve and the question of whether it evolves adaptively.

Gene regulatory divergence between locally adapted ecotypes in their native habitats

Local adaptation is a key driver of ecological specialization and the formation of new species. Despite its importance, the evolution of gene regulatory divergence among locally adapted populations is poorly understood, especially how that divergence manifests in nature. Here, we evaluate gene expression divergence and allele-specific gene expression responses for locally adapted coastal perennial and inland annual accessions of the yellow monkeyflower, Mimulus guttatus, in a field reciprocal transplant experiment. Overall, 6765 (73%) of surveyed genes were differentially expressed between coastal and inland habitats, while 7213 (77%) were differentially expressed between the coastal perennial and inland annual accessions. Cis-regulatory variation was pervasive, affecting 79% (5532) of differentially expressed genes. We detected trans effects for 52% (3611) of differentially expressed genes. Expression plasticity of alleles across habitats (G × E interactions) appears to be relatively common (affecting 18% of transcripts) and could minimize fitness trade-offs at loci that contribute to local adaptation. We also found evidence that at least one chromosomal inversion may act as supergene by holding together haplotypes of differentially expressed genes, but this pattern depends on habitat context. Our results highlight multiple key patterns regarding the relationship between gene expression and the evolution of locally adapted populations.

Protein evolution speed depends on its stability and abundance and on chaperone concentrations

Some biological evolution is slow (millions of years), and some is fast (months to years). The speed at which a protein evolves depends on how stable a protein’s folded structure is, how well it avoids aggregation, and how well-chaperoned it is. What are the mechanisms? We compute fitness landscapes by combining a model of protein-folding equilibria with sequence-change dynamics. We find that adapting to a new environment is fastest for proteins that are least stably folded, because those sit on steep downhill parts of fitness potentials. The modeling shows that cells should adapt to warmer environments faster than to colder ones, explains why increasing a protein’s abundance slows cell evolution, and explains how chaperones accelerate evolution by mitigating this effect.

Proteins evolve at different rates. What drives the speed of protein sequence changes? Two main factors are a protein’s folding stability and aggregation propensity. By combining the hydrophobic–polar (HP) model with the Zwanzig–Szabo–Bagchi rate theory, we find that: (i) Adaptation is strongly accelerated by selection pressure, explaining the broad variation from days to thousands of years over which organisms adapt to new environments. (ii) The proteins that adapt fastest are those that are not very stably folded, because their fitness landscapes are steepest. And because heating destabilizes folded proteins, we predict that cells should adapt faster when put into warmer rather than cooler environments. (iii) Increasing protein abundance slows down evolution (the substitution rate of the sequence) because a typical protein is not perfectly fit, so increasing its number of copies reduces the cell’s fitness. (iv) However, chaperones can mitigate this abundance effect and accelerate evolution (also called evolutionary capacitance) by effectively enhancing protein stability. This model explains key observations about protein evolution rates.

The developmental-genetics of canalization

Canalization, or robustness to genetic or environmental perturbations, is fundamental to complex organisms. While there is strong evidence for canalization as an evolved property that varies among genotypes, the developmental and genetic mechanisms that produce this phenomenon are very poorly understood. For evolutionary biology, understanding how canalization arises is important because, by modulating the phenotypic variation that arises in response to genetic differences, canalization is a determinant of evolvability. For genetics of disease in humans and for economically important traits in agriculture, this subject is important because canalization is a potentially significant cause of missing heritability that confounds genomic prediction of phenotypes. We review the major lines of thought on the developmental-genetic basis for canalization. These fall into two groups. One proposes specific evolved molecular mechanisms while the other deals with robustness or canalization as a more general feature of development. These explanations for canalization are not mutually exclusive and they overlap in several ways. General explanations for canalization are more likely to involve emergent features of development than specific molecular mechanisms. Disentangling these explanations is also complicated by differences in perspectives between genetics and developmental biology. Understanding canalization at a mechanistic level will require conceptual and methodological approaches that integrate quantitative genetics and developmental biology.

Effects of early-life malnutrition on neurodevelopment and neuropsychiatric disorders and the potential mechanisms

Lines of evidence have demonstrated that early-life malnutrition is highly correlated with neurodevelopment and adulthood neuropsychiatric disorders, while some findings are conflicting with each other. In addition, the biological mechanisms are less investigated. We systematically reviewed the evidence linking early-life nutrition status with neurodevelopment and clinical observations in human and animal models. We summarized the effects of special nutritious on neuropsychiatric disorders and explored the underlying potential mechanisms. The further understanding of the biological regulation of early-life nutritional status on neurodevelopment might shed light on precision nutrition at an integrative systems biology framework.

Re-thinking our understanding of immunity: Robustness in the tissue reconstruction system

Robustness, understood as the maintenance of specific functionalities of a given system against internal and external perturbations, is pervasive in today's biology. Yet precise applications of this notion to the immune system have been scarce. Here we show that the concept of robustness sheds light on tissue repair, and particularly on the crucial role the immune system plays in this process. We describe the specific mechanisms, including plasticity and redundancy, by which robustness is achieved in the tissue reconstruction system (TRS). In turn, tissue repair offers a very important test case for assessing the usefulness of the concept of robustness, and identifying different varieties of robustness.

Evolutionary and genetic features of drug targets

In the modern drug discovery pipeline, identification of novel drug targets is a critical step. Despite rapid progress in developing biomedical techniques, it is still a great challenge to find promising new targets from the ample space of human genes. This fact is partially responsible for the situation of "more investments, fewer drugs" in the pharmaceutical industry. A series of recent researches revealed that successfully targeted genes share some common evolutionary and genetic features, which means that the knowledge accumulated in modern evolutionary biology and genetics is very helpful to identify potential drug targets and to find new drugs as well. In this article, we comprehensively summarize the links between human drug targets and genetic diseases and their evolutionary origins, with an attempt to introduce these novel concepts and their medical implications to the biomedical community.

On causal roles and selected effects: our genome is mostly junk

The idea that much of our genome is irrelevant to fitness-is not the product of positive natural selection at the organismal level-remains viable. Claims to the contrary, and specifically that the notion of "junk DNA" should be abandoned, are based on conflating meanings of the word "function". Recent estimates suggest that perhaps 90% of our DNA, though biochemically active, does not contribute to fitness in any sequence-dependent way, and possibly in no way at all. Comparisons to vertebrates with much larger and smaller genomes (the lungfish and the pufferfish) strongly align with such a conclusion, as they have done for the last half-century.

A Functional Link Between Bir1 and the Saccharomyces cerevisiae Ctf19 Kinetochore Complex Revealed Through Quantitative Fitness Analysis

The chromosomal passenger complex (CPC) is a key regulator of eukaryotic cell division, consisting of the protein kinase Aurora B/Ipl1 in association with its activator (INCENP/Sli15) and two additional proteins (Survivin/Bir1 and Borealin/Nbl1). Here, we report a genome-wide genetic interaction screen in Saccharomyces cerevisiae using the bir1-17 mutant, identifying through quantitative fitness analysis deletion mutations that act as enhancers and suppressors. Gene knockouts affecting the Ctf19 kinetochore complex were identified as the strongest enhancers of bir1-17, while mutations affecting the large ribosomal subunit or the mRNA nonsense-mediated decay pathway caused strong phenotypic suppression. Thus, cells lacking a functional Ctf19 complex become highly dependent on Bir1 function and vice versa. The negative genetic interaction profiles of bir1-17 and the cohesin mutant mcd1-1 showed considerable overlap, underlining the strong functional connection between sister chromatid cohesion and chromosome biorientation. Loss of some Ctf19 components, such as Iml3 or Chl4, impacted differentially on bir1-17 compared with mutations affecting other CPC components: despite the synthetic lethality shown by either iml3∆ or chl4∆ in combination with bir1-17, neither gene knockout showed any genetic interaction with either ipl1-321 or sli15-3 Our data therefore imply a specific functional connection between the Ctf19 complex and Bir1 that is not shared with Ipl1.

The Role of DAF-21/Hsp90 in Mouth-Form Plasticity in Pristionchus pacificus

Phenotypic plasticity is increasingly recognized to facilitate adaptive change in plants and animals, including insects, nematodes, and vertebrates. Plasticity can occur as continuous or discrete (polyphenisms) variation. In social insects, for example, in ants, some species have workers of distinct size classes while in other closely related species variation in size may be continuous. Despite the abundance of examples in nature, how discrete morphs are specified remains currently unknown. In theory, polyphenisms might require robustness, whereby the distribution of morphologies would be limited by the same mechanisms that execute buffering from stochastic perturbations, a function attributed to heat-shock proteins of the Hsp90 family. However, this possibility has never been directly tested because plasticity and robustness are considered to represent opposite evolutionary principles. Here, we used a polyphenism of feeding structures in the nematode Pristionchus pacificus to test the relationship between robustness and plasticity using geometric morphometrics of 20 mouth-form landmarks. We show that reducing heat-shock protein activity, which reduces developmental robustness, increases the range of mouth-form morphologies. Specifically, elevated temperature led to a shift within morphospace, pharmacological inhibition of all Hsp90 genes using radicicol treatment increased shape variability in both mouth-forms, and CRISPR/Cas9-induced Ppa-daf-21/Hsp90 knockout had a combined effect. Thus, Hsp90 canalizes the morphologies of plastic traits resulting in discrete polyphenism of mouth-forms.

HSP90 Shapes the Consequences of Human Genetic Variation

HSP90 acts as a protein-folding buffer that shapes the manifestations of genetic variation in model organisms. Whether HSP90 influences the consequences of mutations in humans, potentially modifying the clinical course of genetic diseases, remains unknown. By mining data for >1,500 disease-causing mutants, we found a strong correlation between reduced phenotypic severity and a dominant (HSP90 ≥ HSP70) increase in mutant engagement by HSP90. Examining the cancer predisposition syndrome Fanconi anemia in depth revealed that mutant FANCA proteins engaged predominantly by HSP70 had severely compromised function. In contrast, the function of less severe mutants was preserved by a dominant increase in HSP90 binding. Reducing HSP90's buffering capacity with inhibitors or febrile temperatures destabilized HSP90-buffered mutants, exacerbating FA-related chemosensitivities. Strikingly, a compensatory FANCA somatic mutation from an "experiment of nature" in monozygotic twins both prevented anemia and reduced HSP90 binding. These findings provide one plausible mechanism for the variable expressivity and environmental sensitivity of genetic diseases.

Hsp90 and Physiological Stress Are Linked to Autonomous Transposon Mobility and Heritable Genetic Change in Nematodes

Transposable elements (TEs) have been recognized as potentially powerful drivers of genomic evolutionary change, but factors affecting their mobility and regulation remain poorly understood. Chaperones such as Hsp90 buffer environmental perturbations by regulating protein conformation, but are also part of the PIWI-interacting RNA pathway, which regulates genomic instability arising from mobile TEs in the germline. Stress-induced mutagenesis from TE movement could thus arise from functional trade-offs in the dual roles of Hsp90. We examined the functional constraints of Hsp90 and its role as a regulator of TE mobility by exposing nematodes (Caenorhabditis elegans and Caenorhabditis briggsae) to environmental stress, with and without RNAi-induced silencing of Hsp90. TE excision frequency increased with environmental stress intensity at multiple loci in several strains of each species. These effects were compounded by RNAi-induced knockdown of Hsp90. Mutation frequencies at the unc-22 marker gene in the offspring of animals exposed to environmental stress and Hsp90 RNAi mirrored excision frequency in response to these treatments. Our results support a role for Hsp90 in the suppression of TE mobility, and demonstrate that that the Hsp90 regulatory pathway can be overwhelmed with moderate environmental stress. By compromising genomic stability in germline cells, environmentally induced mutations arising from TE mobility and insertion can have permanent and heritable effects on both the phenotype and the genotype of subsequent generations.

The Molecular Chaperone DnaK Is a Source of Mutational Robustness

Molecular chaperones, also known as heat-shock proteins, refold misfolded proteins and help other proteins reach their native conformation. Thanks to these abilities, some chaperones, such as the Hsp90 protein or the chaperonin GroEL, can buffer the deleterious phenotypic effects of mutations that alter protein structure and function. Hsp70 chaperones use a chaperoning mechanism different from that of Hsp90 and GroEL, and it is not known whether they can also buffer mutations. Here, we show that they can. To this end, we performed a mutation accumulation experiment in Escherichia coli, followed by whole-genome resequencing. Overexpression of the Hsp70 chaperone DnaK helps cells cope with mutational load and completely avoid the extinctions we observe in lineages evolving without chaperone overproduction. Additionally, our sequence data show that DnaK overexpression increases mutational robustness, the tolerance of its clients to nonsynonymous nucleotide substitutions. We also show that this elevated mutational buffering translates into differences in evolutionary rates on intermediate and long evolutionary time scales. Specifically, we studied the evolutionary rates of DnaK clients using the genomes of E. coli, Salmonella enterica, and 83 other gamma-proteobacteria. We find that clients that interact strongly with DnaK evolve faster than weakly interacting clients. Our results imply that all three major chaperone classes can buffer mutations and affect protein evolution. They illustrate how an individual protein like a chaperone can have a disproportionate effect on the evolution of a proteome.

DnaK-Dependent Accelerated Evolutionary Rate in Prokaryotes

Many proteins depend on an interaction with molecular chaperones in order to fold into a functional tertiary structure. Previous studies showed that protein interaction with the GroEL/GroES chaperonine and Hsp90 chaperone can buffer the impact of slightly deleterious mutations in the protein sequence. This capacity of GroEL/GroES to prevent protein misfolding has been shown to accelerate the evolution of its client proteins. Whether other bacterial chaperones have a similar effect on their client proteins is currently unknown. Here, we study the impact of DnaK (Hsp70) chaperone on the evolution of its client proteins. Evolutionary parameters were derived from comparison of the Escherichia coli proteome to 1,808,565 orthologous proteins in 1,149 proteobacterial genomes. Our analysis reveals a significant positive correlation between protein binding frequency with DnaK and evolutionary rate. Proteins with high binding affinity to DnaK evolve on average 4.3-fold faster than proteins in the lowest binding affinity class at the genus resolution. Differences in evolutionary rates of DnaK interactor classes are still significant after adjusting for possible effects caused by protein expression level. Furthermore, we observe an additive effect of DnaK and GroEL chaperones on the evolutionary rates of their common interactors. Finally, we found pronounced similarities in the physicochemical profiles that characterize proteins belonging to DnaK and GroEL interactomes. Our results thus implicate DnaK-mediated folding as a major component in shaping protein evolutionary dynamics in bacteria and supply further evidence for the long-term manifestation of chaperone-mediated folding on genome evolution.

Effects of ambient and preceding temperatures and metabolic genes on flight metabolism in the Glanville fritillary butterfly

Flight is essential for foraging, mate searching and dispersal in many insects, but flight metabolism in ectotherms is strongly constrained by temperature. Thermal conditions vary greatly in natural populations and may hence restrict fitness-related activities. Working on the Glanville fritillary butterfly (Melitaea cinxia), we studied the effects of temperature experienced during the first 2 days of adult life on flight metabolism, genetic associations between flight metabolic rate and variation in candidate metabolic genes, and genotype-temperature interactions. The maximal flight performance was reduced by 17% by 2 days of low ambient temperature (15 °C) prior to the flight trial, mimicking conditions that butterflies commonly encounter in nature. A SNP in phosphoglucose isomerase (Pgi) had a significant association on flight metabolic rate in males and a SNP in triosephosphate isomerase (Tpi) was significantly associated with flight metabolic rate in females. In the Pgi SNP, AC heterozygotes had higher flight metabolic rate than AA homozygotes following low preceding temperature, but the trend was reversed following high preceding temperature, consistent with previous results on genotype-temperature interaction for this SNP. We suggest that these results on 2-day old butterflies reflect thermal effect on the maturation of flight muscles. These results highlight the consequences of variation in thermal conditions on the time scale of days, and they contribute to a better understanding of the complex dynamics of flight metabolism and flight-related activities under conditions that are relevant for natural populations living under variable thermal conditions.

Downregulation of the evolutionary capacitor Hsp90 is mediated by social cues

The relationship between robustness and evolvability is a long-standing question in evolution. Heat shock protein 90 (HSP90), a molecular chaperone, has been identified as a potential capacitor for evolution, since it allows for the accumulation and release of cryptic genetic variation, and also for the regulation of novel genetic variation through transposon activity. However, to date, it is unknown whether Hsp90 expression is regulated upon demand (i.e. when the release of cryptic genetic variation is most needed). Here, we show that Hsp90 has reduced transcription under conditions where the mobilization of genetic variation could be advantageous. We designed a situation that indicates a stressful environment but avoids the direct effects of stress, by placing untreated (focal) red flour beetles, Tribolium castaneum, into groups together with wounded conspecifics, and found a consistent reduction in expression of two Hsp90 genes (Hsp83 and Hsp90) in focal beetles. We moreover observed a social transfer of immunity in this non-eusocial insect: there was increased activity of the phenoloxidase enzyme and downregulation of the immune regulator, imd. Our study poses the exciting question of whether evolvability might be regulated through the use of information derived from the social environment.

Molecular basis of canalization in an ascidian species complex adapted to different thermal conditions

Canalization is a result of intrinsic developmental buffering that ensures phenotypic robustness under genetic variation and environmental perturbation. As a consequence, animal phenotypes are remarkably consistent within a species under a wide range of conditions, a property that seems contradictory to evolutionary change. Study of laboratory model species has uncovered several possible canalization mechanisms, however, we still do not understand how the level of buffering is controlled in natural populations. We exploit wild populations of the marine chordate Ciona intestinalis to show that levels of buffering are maternally inherited. Comparative transcriptomics show expression levels of genes encoding canonical chaperones such as Hsp70 and Hsp90 do not correlate with buffering. However the expression of genes encoding endoplasmic reticulum (ER) chaperones does correlate. We also show that ER chaperone genes are widely conserved amongst animals. Contrary to previous beliefs that expression level of Heat Shock Proteins (HSPs) can be used as a measurement of buffering levels, we propose that ER associated chaperones comprise a cellular basis for canalization. ER chaperones have been neglected by the fields of development, evolution and ecology, but their study will enhance understanding of both our evolutionary past and the impact of global environmental change.