Molecular-functional studies of adaptive genetic variation in prokaryotes and eukaryotes

Effect of phenol formaldehyde-associated microplastics on soil microbial community, assembly, and functioning

Increasing investigations explore the effects of plastic pollutants on bacterial communities, diversity, and functioning in various ecosystems. However, the impact of microplastics (MPs) on the eukaryotic community, microbial assemblages, and interactions is still limited. Here, we investigated bacterial and micro-eukaryotic communities and functioning in soils with different concentrations of phenol formaldehyde-associated MPs (PF-MPs), and revealed the factors, such as soil properties, microbial community assembly, and interactions between microbes, influencing them. Our results showed that a high concentration (1%) of PF-MPs decreased the microbial interactions and the contribution of deterministic processes to the community assembly of microbes, and consequently changed the communities of bacteria, but not eukaryotes. A significant and negative relationship was determined between N₂O emission rate and functional genes related to nitrification, indicating that the competitive interactions between functional microbes would affect the nitrogen cycling of soil ecosystem. We further found that vegetable biomass weakly decreased in treatments with a higher concentration of PF-MPs and positively related to the diversity of micro-eukaryotic communities and functional diversity of bacterial communities. These results suggest that a high concentration of the PF-MPs would influence crop growth by changing microbial communities, interactions, and eukaryotic and functional diversity. Our findings provide important evidence for agriculture management of phenol formaldehyde and suggest that we must consider their threats to microbial community compositions, diversity, and assemblage in soils due to the accumulation of PF-MPs widely used in the field.

The evolution of multi-gene families and metabolic pathways in the evening primroses (Oenothera: Onagraceae): A comparative transcriptomics approach

The plant genus Oenothera has played an important role in the study of plant evolution of genomes and plant defense and reproduction. Here, we build on the 1kp transcriptomic dataset by creating 44 new transcriptomes and analyzing a total of 63 transcriptomes to present a large-scale comparative study across 29 Oenothera species. Our dataset included 30.4 million reads per individual and 2.3 million transcripts on average. We used this transcriptome resource to examine genome-wide evolutionary patterns and functional diversification by searching for orthologous genes and performed gene family evolution analysis. We found wide heterogeneity in gene family evolution across the genus, with section Oenothera exhibiting the most pronounced evolutionary changes. Overall, more significant gene family expansions occurred than contractions. We also analyzed the molecular evolution of phenolic metabolism by retrieving proteins annotated for phenolic enzymatic complexes. We identified 1,568 phenolic genes arranged into 83 multigene families that varied widely across the genus. All taxa experienced rapid phenolic evolution (fast rate of genomic turnover) involving 33 gene families, which exhibited large expansions, gaining about 2-fold more genes than they lost. Upstream enzymes phenylalanine ammonia-lyase (PAL) and 4-coumaroyl: CoA ligase (4CL) accounted for most of the significant expansions and contractions. Our results suggest that adaptive and neutral evolutionary processes have contributed to Oenothera diversification and rapid gene family evolution.

Ecological Load and Balancing Selection in Circumboreal Barnacles

Acorn barnacle adults experience environmental heterogeneity at various spatial scales of their circumboreal habitat, raising the question of how adaptation to high environmental variability is maintained in the face of strong juvenile dispersal and mortality. Here, we show that 4% of genes in the barnacle genome experience balancing selection across the entire range of the species. Many of these genes harbor mutations maintained across 2 My of evolution between the Pacific and Atlantic oceans. These genes are involved in ion regulation, pain reception, and heat tolerance, functions which are essential in highly variable ecosystems. The data also reveal complex population structure within and between basins, driven by the trans-Arctic interchange and the last glaciation. Divergence between Atlantic and Pacific populations is high, foreshadowing the onset of allopatric speciation, and suggesting that balancing selection is strong enough to maintain functional variation for millions of years in the face of complex demography.

The Variety and Inscrutability of Polar Environments as a Resource of Biotechnologically Relevant Molecules

The application of an ever-increasing number of methodological approaches and tools is positively contributing to the development and yield of bioprospecting procedures. In this context, cold-adapted bacteria from polar environments are becoming more and more intriguing as valuable sources of novel biomolecules, with peculiar properties to be exploited in a number of biotechnological fields. This review aims at highlighting the biotechnological potentialities of bacteria from Arctic and Antarctic habitats, both biotic and abiotic. In addition to cold-enzymes, which have been intensively analysed, relevance is given to recent advances in the search for less investigated biomolecules, such as biosurfactants, exopolysaccharides and antibiotics.

Biological Effects of Single-Nucleotide Polymorphisms in the Drosophila melanogaster Malic Enzyme Locus

A pair of amino acid polymorphisms within the Drosophila melanogaster Malic enzyme (Men) locus presents an interesting case of genetic variation that appears to be under selection. The two alleles at each site are biochemically distinct, but their biological effects are unknown. One polymorphic site is near the active site and the other is buried within the protein. Strikingly, in twelve different populations, the first polymorphism is always found at approximately a 50:50 allelic frequency, whereas the second polymorphism is always found at approximately 90:10. The consistency of the frequencies between populations suggests that the polymorphisms are under selection and it is possible that balancing selection is at play. We used 16 lines of flies to create the nine genotypes needed to quantify both effects of the polymorphic sites and possible genetic background effects, which we found to be widespread. The alleles at each site differ, but in different biochemical characteristics. The first site significantly influences MEN K_m and V_max, whereas the second site affects the K_m and the V_max/K_m ratio (relative activity). Interestingly, the rarest allele is the most biochemically distinct. We also assayed three more distal phenotypes, triglyceride concentration, carbohydrate concentration, and longevity. In all cases, the phenotypes of the heterozygous genotypes are intermediate between those of the respective homozygotes suggesting that if balancing selection is maintaining the observed allele frequencies it is not through non-linear combinations of the biochemical phenotypes.

Sequence-Based Analysis of Thermal Adaptation and Protein Energy Landscapes in an Invasive Blue Mussel (Mytilus galloprovincialis)

Adaptive responses to thermal stress in poikilotherms plays an important role in determining competitive ability and species distributions. Amino acid substitutions that affect protein stability and modify the thermal optima of orthologous proteins may be particularly important in this context. Here, we examine a set of 2,770 protein-coding genes to determine if proteins in a highly invasive heat tolerant blue mussel (Mytilus galloprovincialis) contain signals of adaptive increases in protein stability relative to orthologs in a more cold tolerant M. trossulus. Such thermal adaptations might help to explain, mechanistically, the success with which the invasive marine mussel M. galloprovincialis has displaced native species in contact zones in the eastern (California) and western (Japan) Pacific. We tested for stabilizing amino acid substitutions in warm tolerant M. galloprovincialis relative to cold tolerant M. trossulus with a generalized linear model that compares in silico estimates of recent changes in protein stability among closely related congeners. Fixed substitutions in M. galloprovincialis were 3,180.0 calories per mol per substitution more stabilizing at genes with both elevated d_N/d_S ratios and transcriptional responses to heat stress, and 705.8 calories per mol per substitution more stabilizing across all 2,770 loci investigated. Amino acid substitutions concentrated in a small number of genes were more stabilizing in M. galloprovincialis compared with cold tolerant M. trossulus. We also tested for, but did not find, enrichment of a priori GO terms in genes with elevated d_N/d_S ratios in M. galloprovincialis. This might indicate that selection for thermodynamic stability is generic across all lineages, and suggests that the high change in estimated protein stability that we observed in M. galloprovincialis is driven by selection for extra stabilizing substitutions, rather than by higher incidence of selection in a greater number of genes in this lineage. Nonetheless, our finding of more stabilizing amino acid changes in the warm adapted lineage is important because it suggests that adaption for thermal stability has contributed to M. galloprovincialis’ superior tolerance to heat stress, and that pairing tests for positive selection and tests for transcriptional response to heat stress can identify candidates of protein stability adaptation.

RNA thermosensors: how might animals exploit their regulatory potential?

The secondary and tertiary orders of RNA structure are crucial for a suite of RNA-related functions, including regulation of translation, gene expression and RNA turnover. The temperature sensitivity of RNA secondary and tertiary structures is exploited by bacteria to fabricate RNA thermosensing systems that allow a rapid adaptive response to temperature change. RNA thermometers (RNATs) present in non-coding regions of certain mRNAs of pathogenic bacteria enable rapid upregulation of translation of virulence proteins when the temperature of the bacterium rises after entering a mammalian host. Rapid upregulation of translation of bacterial heat-shock proteins likewise is governed in part by RNATs. Turnover of mRNA may be regulated by temperature-sensitive RNA structures. Whereas the roles of temperature-sensitive RNA structures similar to RNATs in Eukarya and Archaea are largely unknown, there would appear to be a potential for all taxa to adaptively regulate their thermal physiology through exploitation of RNA-based thermosensory responses akin to those of bacteria. In animals, these responses might include regulation of translation of stress-induced proteins, alternative splicing of messenger RNA precursors, differential expression of allelic proteins, modulation of activities of small non-coding RNAs, regulation of mRNA turnover and control of RNA editing. New methods for predicting, detecting and experimentally modifying RNA secondary structure offer promising windows into these fascinating aspects of RNA biochemistry. Elucidating whether animals too have exploited the types of RNA thermosensing tools that are used so effectively by bacteria seems likely to provide exciting new insights into the mechanisms of evolutionary adaptation and acclimatization to temperature.

Global study of holistic morphological effectors in the budding yeast Saccharomyces cerevisiae

Background

The size of the phenotypic effect of a gene has been thoroughly investigated in terms of fitness and specific morphological traits in the budding yeast Saccharomyces cerevisiae, but little is known about gross morphological abnormalities.

Results

We identified 1126 holistic morphological effectors that cause severe gross morphological abnormality when deleted, and 2241 specific morphological effectors with weak holistic effects but distinctive effects on yeast morphology. Holistic effectors fell into many gene function categories and acted as network hubs, affecting a large number of morphological traits, interacting with a large number of genes, and facilitating high protein expression. Holistic morphological abnormality was useful for estimating the importance of a gene to morphology. The contribution of gene importance to fitness and morphology could be used to efficiently classify genes into functional groups.

Conclusion

Holistic morphological abnormality can be used as a reproducible and reliable gene feature for high-dimensional morphological phenotyping. It can be used in many functional genomic applications.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4526-z) contains supplementary material, which is available to authorized users.

Hubby and Lewontin on Protein Variation in Natural Populations: When Molecular Genetics Came to the Rescue of Population Genetics

The 1966 GENETICS papers by John Hubby and Richard Lewontin were a landmark in the study of genome-wide levels of variability. They used the technique of gel electrophoresis of enzymes and proteins to study variation in natural populations of Drosophila pseudoobscura, at a set of loci that had been chosen purely for technical convenience, without prior knowledge of their levels of variability. Together with the independent study of human populations by Harry Harris, this seminal study provided the first relatively unbiased picture of the extent of genetic variability in protein sequences within populations, revealing that many genes had surprisingly high levels of diversity. These papers stimulated a large research program that found similarly high electrophoretic variability in many different species and led to statistical tools for interpreting the data in terms of population genetics processes such as genetic drift, balancing and purifying selection, and the effects of selection on linked variants. The current use of whole-genome sequences in studies of variation is the direct descendant of this pioneering work.

The mechanism of transactivation regulation due to polymorphic short tandem repeats (STRs) using IGF1 promoter as a model

Functional short tandem repeats (STR) are polymorphic in the population, and the number of repeats regulates the expression of nearby genes (known as expression STR, eSTR). STR in IGF1 promoter has been extensively studied for its association with IGF1 concentration in blood and various clinical traits and represents an important eSTR. We previously used an in-vitro luciferase reporter model to examine the interaction between STRs and SNPs in IGF1 promoter. Here, we further explored the mechanism how the number of repeats of the STR regulates gene transcription. An inverse correlation between the number of repeats and the extent of transactivation was found in a haplotype consisting of three promoter SNPs (C-STR-T-T). We showed that these adjacent SNPs located outside the STR were required for the STR to function as eSTR. The C allele of rs35767 provides a binding site for CCAAT/enhancer-binding-protein δ (C/EBPD), which is essential for the gradational transactivation property of eSTR and FOXA3 may also be involved. Therefore, we propose a mechanism in which the gradational transactivation by the eSTR is caused by the interaction of one or more transcriptional complexes located outside the STR, rather than by direct binding to a repeat motif of the STR.

New views on the selection acting on genetic polymorphism in central metabolic genes

Studies of the polymorphism of central metabolic genes as a source of fitness variation in natural populations date back to the discovery of allozymes in the 1960s. The unique features of these genes and their enzymes and our knowledge base greatly facilitates the systems-level study of this group. The expectation that pathway flux control is central to understanding the molecular evolution of genes is discussed, as well as studies that attempt to place gene-specific molecular evolution and polymorphism into a context of pathway and network architecture. There is an increasingly complex picture of the metabolic genes assuming additional roles beyond their textbook anabolic and catabolic reactions. In particular, this review emphasizes the potential role of these genes as part of the energy-sensing machinery. It is underscored that the concentrations of key cellular metabolites are the reflections of cellular energy status and nutritional input. These metabolites are the top-down signaling messengers that set signaling through signaling pathways that are involved in energy economy. I propose that the polymorphisms in central metabolic genes shift metabolite concentrations and in that fashion act as genetic modifiers of the energy-state coupling to the transcriptional networks that affect physiological trade-offs with significant fitness consequences.

Effector Diversification Contributes to Xanthomonas oryzae pv. oryzae Phenotypic Adaptation in a Semi-Isolated Environment

Understanding the processes that shaped contemporary pathogen populations in agricultural landscapes is quite important to define appropriate management strategies and to support crop improvement efforts. Here, we took advantage of an historical record to examine the adaptation pathway of the rice pathogen Xanthomonas oryzae pv. oryzae (Xoo) in a semi-isolated environment represented in the Philippine archipelago. By comparing genomes of key Xoo groups we showed that modern populations derived from three Asian lineages. We also showed that diversification of virulence factors occurred within each lineage, most likely driven by host adaptation, and it was essential to shape contemporary pathogen races. This finding is particularly important because it expands our understanding of pathogen adaptation to modern agriculture.

Evolutionary history of a dispersal-associated locus across sympatric and allopatric divergent populations of a wing-polymorphic beetle across Atlantic Europe

Studying the evolutionary history of trait divergence, in particular those related to dispersal capacity, is of major interest for the process of local adaptation and metapopulation dynamics. Here, we reconstruct the evolution of different alleles at the nuclear-encoded mitochondrial NADP(+)-dependent isocitrate dehydrogenase (mtIdh) locus of the ground beetle Pogonus chalceus that are differentially and repeatedly selected in short- and long-winged populations in response to different hydrological regimes at both allopatric and sympatric scales along the Atlantic European coasts. We sequenced 2788 bp of the mtIdh locus spanning a ~7-kb genome region and compared its variation with that of two supposedly neutral genes. mtIdh sequences show (i) monophyletic clustering of the short-winged associated mtIDH-DE haplotypes within the long-winged associated mtIDH-AB haplotypes, (ii) a more than tenfold lower haplotype diversity associated with the mtIDH-DE alleles compared to the mtIDH-AB alleles and (iii) a high number of fixed nucleotide differences between both mtIDH haplotype clusters. Coalescent simulations suggest that this observed sequence variation in the mtIdh locus is most consistent with a singular origin in a partially isolated subpopulation, followed by a relatively recent spread of the mtIDH-DE allele in short-winged populations along the Atlantic coast. These results demonstrate that even traits associated with decreased dispersal capacity can rapidly spread and that reuse of adaptive alleles plays an important role in the adaptive potential within this sympatric mosaic of P. chalceus populations.

Patterns of variation during adaptation in functionally linked loci

An understanding of the distribution of natural patterns of genetic variation is relevant to such fundamental biological fields as evolution and development. One recent approach to understanding such patterns has been to focus on the constraints that may arise as a function of the network or pathway context in which genes are embedded. Despite theoretical expectations of higher evolutionary constraint for genes encoding upstream versus downstream enzymes in metabolic pathways, empirical results have varied. Here we combine two complementary models from population genetics and enzyme kinetics to explore genetic variation as a function of pathway position when selection acts on whole-pathway flux. We are able to qualitatively reproduce empirically observed patterns of polymorphism and divergence and suggest that expectations should vary depending on the evolutionary trajectory of a population. Upstream genes are initially more polymorphic and diverge faster after an environmental change, while we see the opposite trend as the population approaches its fitness optimum.

Hybrid incompatibility arises in a sequence-based bioenergetic model of transcription factor binding

Postzygotic isolation between incipient species results from the accumulation of incompatibilities that arise as a consequence of genetic divergence. When phenotypes are determined by regulatory interactions, hybrid incompatibility can evolve even as a consequence of parallel adaptation in parental populations because interacting genes can produce the same phenotype through incompatible allelic combinations. We explore the evolutionary conditions that promote and constrain hybrid incompatibility in regulatory networks using a bioenergetic model (combining thermodynamics and kinetics) of transcriptional regulation, considering the bioenergetic basis of molecular interactions between transcription factors (TFs) and their binding sites. The bioenergetic parameters consider the free energy of formation of the bond between the TF and its binding site and the availability of TFs in the intracellular environment. Together these determine fractional occupancy of the TF on the promoter site, the degree of subsequent gene expression and in diploids, and the degree of dominance among allelic interactions. This results in a sigmoid genotype-phenotype map and fitness landscape, with the details of the shape determining the degree of bioenergetic evolutionary constraint on hybrid incompatibility. Using individual-based simulations, we subjected two allopatric populations to parallel directional or stabilizing selection. Misregulation of hybrid gene expression occurred under either type of selection, although it evolved faster under directional selection. Under directional selection, the extent of hybrid incompatibility increased with the slope of the genotype-phenotype map near the derived parental expression level. Under stabilizing selection, hybrid incompatibility arose from compensatory mutations and was greater when the bioenergetic properties of the interaction caused the space of nearly neutral genotypes around the stable expression level to be wide. F2's showed higher hybrid incompatibility than F1's to the extent that the bioenergetic properties favored dominant regulatory interactions. The present model is a mechanistically explicit case of the Bateson-Dobzhansky-Muller model, connecting environmental selective pressure to hybrid incompatibility through the molecular mechanism of regulatory divergence. The bioenergetic parameters that determine expression represent measurable properties of transcriptional regulation, providing a predictive framework for empirical studies of how phenotypic evolution results in epistatic incompatibility at the molecular level in hybrids.

Mapping the fitness landscape of gene expression uncovers the cause of antagonism and sign epistasis between adaptive mutations

How do adapting populations navigate the tensions between the costs of gene expression and the benefits of gene products to optimize the levels of many genes at once? Here we combined independently-arising beneficial mutations that altered enzyme levels in the central metabolism of Methylobacterium extorquens to uncover the fitness landscape defined by gene expression levels. We found strong antagonism and sign epistasis between these beneficial mutations. Mutations with the largest individual benefit interacted the most antagonistically with other mutations, a trend we also uncovered through analyses of datasets from other model systems. However, these beneficial mutations interacted multiplicatively (i.e., no epistasis) at the level of enzyme expression. By generating a model that predicts fitness from enzyme levels we could explain the observed sign epistasis as a result of overshooting the optimum defined by a balance between enzyme catalysis benefits and fitness costs. Knowledge of the phenotypic landscape also illuminated that, although the fitness peak was phenotypically far from the ancestral state, it was not genetically distant. Single beneficial mutations jumped straight toward the global optimum rather than being constrained to change the expression phenotypes in the correlated fashion expected by the genetic architecture. Given that adaptation in nature often results from optimizing gene expression, these conclusions can be widely applicable to other organisms and selective conditions. Poor interactions between individually beneficial alleles affecting gene expression may thus compromise the benefit of sex during adaptation and promote genetic differentiation.

The pace and outcome of a series of adaptive steps in an evolving lineage depends upon how well different beneficial mutations stack on top of each other. We found that independent beneficial mutations that affected gene expression for a metabolic pathway did not work well together, and were often jointly deleterious. The most beneficial mutations interacted the most poorly with others, which was a trend we found common in other biological systems. Through generating a model that accounted for enzymatic benefits and expression costs, we uncovered that this antagonism was caused by a phenotype to fitness mapping that had an intermediate peak. This allowed us to predict the fitness effect of double mutants and to uncover that the single winning mutations tended to move straight to the peak in a single step. These findings demonstrate the importance of considering the phenotypic changes that cause nonlinear interactions between mutations upon fitness, and thus influence how populations evolve.

Adherence to Bürgi-Dunitz stereochemical principles requires significant structural rearrangements in Schiff-base formation: insights from transaldolase complexes

The Bürgi-Dunitz angle (αBD) describes the trajectory of approach of a nucleophile to an electrophile. The adoption of a stereoelectronically favorable αBD can necessitate significant reactive-group repositioning over the course of bond formation. In the context of enzyme catalysis, interactions with the protein constrain substrate rotation, which could necessitate structural transformations during bond formation. To probe this theoretical framework vis-à-vis biocatalysis, Schiff-base formation was analysed in Francisella tularensis transaldolase (TAL). Crystal structures of wild-type and Lys→Met mutant TAL in covalent and noncovalent complexes with fructose 6-phosphate and sedoheptulose 7-phosphate clarify the mechanism of catalysis and reveal that substrate keto moieties undergo significant conformational changes during Schiff-base formation. Structural changes compelled by the trajectory considerations discussed here bear relevance to bond formation in a variety of constrained enzymic/engineered systems and can inform the design of covalent therapeutics.

Nature's inordinate fondness for metabolic enzymes: why metabolic enzyme loci are so frequently targets of selection

Metabolic enzyme loci were some of the first genes accessible for molecular evolution and ecology research. New technologies now make the whole genome, transcriptome or proteome readily accessible, allowing unbiased scans for loci exhibiting significant differences in allele frequency or expression level and associated with phenotypes and/or responses to natural selection. With surprising frequency and in many cases in proportions greater than chance relative to other genes, glycolysis and TCA cycle enzyme loci appear among the genes with significant associations in these studies. Hence, there is an ongoing need to understand the basis for fitness effects of metabolic enzyme polymorphisms. Allele-specific effects on the binding affinity and catalytic rate of individual enzymes are well known, but often of uncertain significance because metabolic control theory and in vivo studies indicate that many individual metabolic enzymes do not affect pathway flux rate. I review research, so far little used in evolutionary biology, showing that metabolic enzyme substrates affect signalling pathways that regulate cell and organismal biology, and that these enzymes have moonlighting functions. To date there is little knowledge of how alleles in natural populations affect these phenotypes. I discuss an example in which alleles of a TCA enzyme locus associate with differences in a signalling pathway and development, organismal performance, and ecological dynamics. Ultimately, understanding how metabolic enzyme polymorphisms map to phenotypes and fitness remains a compelling and ongoing need for gaining robust knowledge of ecological and evolutionary processes.

Evolution of tryptophan biosynthetic pathway in microbial genomes: a comparative genetic study

Biosynthetic pathway evolution needs to consider the evolution of a group of genes that code for enzymes catalysing the multiple chemical reaction steps leading to the final end product. Tryptophan biosynthetic pathway has five chemical reaction steps that are highly conserved in diverse microbial genomes, though the genes of the pathway enzymes show considerable variations in arrangements, operon structure (gene fusion and splitting) and regulation. We use a combined bioinformatic and statistical analyses approach to address the question if the pathway genes from different microbial genomes, belonging to a wide range of groups, show similar evolutionary relationships within and between them. Our analyses involved detailed study of gene organization (fusion/splitting events), base composition, relative synonymous codon usage pattern of the genes, gene expressivity, amino acid usage, etc. to assess inter- and intra-genic variations, between and within the pathway genes, in diverse group of microorganisms. We describe these genetic and genomic variations in the tryptophan pathway genes in different microorganisms to show the similarities across organisms, and compare the same genes across different organisms to find the possible variability arising possibly due to horizontal gene transfers. Such studies form the basis for moving from single gene evolution to pathway evolutionary studies that are important steps towards understanding the systems biology of intracellular pathways.

Positive selection in glycolysis among Australasian stick insects

BACKGROUND

The glycolytic pathway is central to cellular energy production. Selection on individual enzymes within glycolysis, particularly phosphoglucose isomerase (Pgi), has been associated with metabolic performance in numerous organisms. Nonetheless, how whole energy-producing pathways evolve to allow organisms to thrive in different environments and adopt new lifestyles remains little explored. The Lanceocercata radiation of Australasian stick insects includes transitions from tropical to temperate climates, lowland to alpine habitats, and winged to wingless forms. This permits a broad investigation to determine which steps within glycolysis and what sites within enzymes are the targets of positive selection. To address these questions we obtained transcript sequences from seven core glycolysis enzymes, including two Pgi paralogues, from 29 Lanceocercata species.

RESULTS

Using maximum likelihood methods a signature of positive selection was inferred in two core glycolysis enzymes. Pgi and Glyceraldehyde 3-phosphate dehydrogenase (Gaphd) genes both encode enzymes linking glycolysis to the pentose phosphate pathway. Positive selection among Pgi paralogues and orthologues predominately targets amino acids with residues exposed to the protein's surface, where changes in physical properties may alter enzyme performance.

CONCLUSION

Our results suggest that, for Lancerocercata stick insects, adaptation to new stressful lifestyles requires a balance between maintaining cellular energy production, efficiently exploiting different energy storage pools and compensating for stress-induced oxidative damage.

Specific-gene studies of evolutionary mechanisms in an age of genome-wide surveying

The molecular tools of genomics have great power to reveal patterns of genetic difference within or among species, but must be complemented by the mechanistic study of the genetic variants found if these variants' evolutionary meaning is to be well understood. Central to this purpose is knowledge of the organisms' genotype-phenotype-environment interactions, which embody biological adaptation and constraint and thus drive natural selection. The history of this approach is briefly reviewed. Strategies embracing the complementarity of genomics and specific-gene studies in evolution are considered. Implementation of these strategies, and examples showing their feasibility and power, are discussed. Initial generalizations emphasize: (1) reproducibility of adaptive mechanisms; (2) evolutionary co-importance of variation in protein sequences and expression; (3) refinement of rudimentary molecular functions as an origin of evolutionary innovations; (4) identification of specific-gene mechanisms as underpinnings of genomic or quantitative genetic variation; and (5) multiple forms of adaptive or constraining epistasis among genes. Progress along these lines will advance understanding of evolution and support its use in addressing urgent medical and environmental applications.

A genetic polymorphism evolving in parallel in two cell compartments and in two clades

BACKGROUND

The enzyme phosphoenolpyruvate carboxykinase, PEPCK, occurs in its guanosine-nucleotide-using form in animals and a few prokaryotes. We study its natural genetic variation in Colias (Lepidoptera, Pieridae). PEPCK offers a route, alternative to pyruvate kinase, for carbon skeletons to move between cytosolic glycolysis and mitochondrial Krebs cycle reactions.

RESULTS

PEPCK is expressed in both cytosol and mitochondrion, but differently in diverse animal clades. In vertebrates and independently in Drosophila, compartment-specific paralogous genes occur. In a contrasting expression strategy, compartment-specific PEPCKs of Colias and of the silkmoth, Bombyx, differ only in their first, 5', exons; these are alternatively spliced onto a common series of following exons. In two Colias species from distinct clades, PEPCK sequence is highly variable at nonsynonymous and synonymous sites, mainly in its common exons. Three major amino acid polymorphisms, Gly 335 ↔ Ser, Asp 503 ↔ Glu, and Ile 629 ↔ Val occur in both species, and in the first two cases are similar in frequency between species. Homology-based structural modelling shows that the variants can alter hydrogen bonding, salt bridging, or van der Waals interactions of amino acid side chains, locally or at one another's sites which are distant in PEPCK's structure, and thus may affect its enzyme function. We ask, using coalescent simulations, if these polymorphisms' cross-species similarities are compatible with neutral evolution by genetic drift, but find the probability of this null hypothesis is 0.001 ≤ P ≤ 0.006 under differing scenarios.

CONCLUSION

Our results make the null hypothesis of neutrality of these PEPCK polymorphisms quite unlikely, but support an alternative hypothesis that they are maintained by natural selection in parallel in the two species. This alternative can now be justifiably tested further via studies of PEPCK genotypes' effects on function, organismal performance, and fitness. This case emphasizes the importance, for evolutionary insight, of studying gene-specific mechanisms affected by natural genetic variation as an essential complement to surveys of such variation.

Evolution of flux control in the glucosinolate pathway in Arabidopsis thaliana

Network characteristics of biochemical pathways are believed to influence the rate of evolutionary change in constituent enzymes. One characteristic that may affect rate heterogeneity is control of the amount of product produced by a biochemical pathway or flux control. In particular, theoretical analyses suggest that adaptive substitutions should be concentrated in the enzyme(s) that exert the greatest control over flux. Although a handful of studies have found a correlation between position in a pathway and evolutionary rate, these investigations have not examined the relationship between evolutionary rate and flux control. Given that genes with greater control will experience stronger selection and that the probability of fixation is proportional to the selective advantage, we ask the following: 1) do upstream enzymes have majority flux control, 2) do enzymes with majority flux control accumulate adaptive substitutions, and 3) are upstream enzymes under higher selective constraint? First, by perturbing the enzymes in the aliphatic glucosinolate pathway in Arabidopsis thaliana with gene insertion lines, we show that flux control is focused in the first enzyme in the pathway. Next, by analyzing several sequence signatures of selection, we also show that this enzyme is the only one in the pathway that shows convincing evidence of selection. Our results support the hypothesis that natural selection preferentially acts on enzymes with high flux control.

The relationship between the hierarchical position of proteins in the human signal transduction network and their rate of evolution

Background

Proteins evolve at disparate rates, as a result of the action of different types and strengths of evolutionary forces. An open question in evolutionary biology is what factors are responsible for this variability. In general, proteins whose function has a great impact on organisms’ fitness are expected to evolve under stronger selective pressures. In biosynthetic pathways, upstream genes usually evolve under higher levels of selective constraint than those acting at the downstream part, as a result of their higher hierarchical position. Similar observations have been made in transcriptional regulatory networks, whose upstream elements appear to be more essential and subject to selection. Less well understood is, however, how selective pressures distribute along signal transduction pathways.

Results

Here, I combine comparative genomics and directed protein interaction data to study the distribution of evolutionary forces across the human signal transduction network. Surprisingly, no evidence was found for higher levels of selective constraint at the upstream network genes (those occupying more hierarchical positions). On the contrary, purifying selection was found to act more strongly on genes acting at the downstream part of the network, which seems to be due to downstream genes being more highly and broadly expressed, performing certain functions and, in particular, encoding proteins that are more highly connected in the protein–protein interaction network. When the effect of these confounding factors is discounted, upstream and downstream genes evolve at similar rates. The trends found in the overall signaling network are exemplified by analysis of the distribution of purifying selection along the mammalian Ras signaling pathway, showing that upstream and downstream genes evolve at similar rates.

Conclusions

These results indicate that the upstream/downstream position of proteins in the signal transduction network has, in general, no direct effect on their rates of evolution, suggesting that upstream and downstream genes are similarly important for the function of the network. This implies that natural selection differently distributes across signal transduction networks and across biosynthetic and transcriptional regulatory networks, which might reflect fundamental differences in their function and organization.

The evolution of genes in branched metabolic pathways

Simulation models of the evolution of genes in a branched metabolic pathway subject to stabilizing selection on flux are described and analyzed. The models are based either on metabolic control theory (MCT), with the assumption that enzymes are far from saturation, or on Michaelis-Menten kinetics, which allows for saturation and near saturation. Several predictions emerge from the models: (1) flux control evolves to be concentrated at pathway branch points, including the first enzyme in the pathway. (2) When flux is far from its optimum, adaptive substitutions occur disproportionately often in branching enzymes. (3) When flux is near its optimum, adaptive substitutions occur disproportionately often in nonbranching enzymes. (4) Slightly deleterious substitutions occur disproportionately often in nonbranching enzymes. (5) In terms of both flux control and patterns of substitution, pathway branches are similar to those predicted for linear pathways. These predictions provide null hypotheses for empirical examination of the evolution of genes in metabolic pathways.

The microbial diversity of Polar environments is a fertile ground for bioprospecting

The term bioprospecting has been adopted for systematic searches in nature for new bioactive compounds, genes, proteins, microorganisms and other products with potential for commercial use. Much effort has been focused on microorganisms able to thrive under harsh conditions, including the Polar environments. Both the lipid and protein cellular building blocks of Polar microorganisms are shaped by their adaptation to the permanently low temperatures. In addition, strongly differing environments, such as permafrost, glaciers and sea ice, have contributed to additional functional diversity. Emerging massive-parallel sequencing technologies have revealed the existence of a huge, hitherto unseen diversity of low-abundance phylotypes--the rare biosphere--even in the Polar environments. This realization has further strengthened the need to employ cultivation-independent approaches, including metagenomics and single-cell genomic sequencing, to get comprehensive access to the genetic diversity of microbial communities for bioprospecting purposes. In this review, we present an updated snapshot of recent findings on the molecular basis for adaptation to the cold and the phylogenetic diversities of different Polar environments. Novel approaches in bioprospecting are presented and we conclude by showing recent bioprospecting outcomes in terms of new molecules patented or applied by some biotech companies.

Experimental approaches to evaluate the contributions of candidate protein-coding mutations to phenotypic evolution

Identifying mechanisms of molecular adaptation can provide important insights into the process of phenotypic evolution, but it can be exceedingly difficult to quantify the phenotypic effects of specific mutational changes. To verify the adaptive significance of genetically based changes in protein function, it is necessary to document functional differences between the products of derived and wild-type alleles and to demonstrate that such differences impinge on higher-level physiological processes (and ultimately, fitness). In the case of metabolic enzymes, this requires documenting in vivo differences in reaction rate that give rise to differences in flux through the pathway in which the enzymes function. These measured differences in pathway flux should then give rise to differences in cellular or systemic physiology that affect fitness-related variation in whole-organism performance. Efforts to establish these causal connections between genotype, phenotype, and fitness require experiments that carefully control for environmental variation and background genetic variation. Here, we discuss experimental approaches to evaluate the contributions of amino-acid mutations to adaptive phenotypic change. We discuss conceptual and methodological issues associated with in vitro and in vivo studies of protein function, and the evolutionary insights that can be gleaned from such studies. We also discuss the importance of isolating the effects of individual mutations to distinguish between positively selected substitutions that directly contribute to improvements in protein function versus positively selected, compensatory substitutions that mitigate negative pleiotropic effects of antecedent changes.

The biochemical basis of life history adaptation: molecular and enzymological causes of NADP(+)-isocitrate dehydrogenase activity differences between morphs of Gryllus firmus that differ in lipid biosynthesis and life history

Although whole-organism aspects of life-history physiology are well studied and molecular information (e.g., transcript abundance) on life-history variation is accumulating rapidly, much less information is available on the biochemical (enzymological) basis of life-history adaptation. The present study investigated the biochemical and molecular causes of specific activity differences of the lipogenic enzyme, NADP(+)-isocitrate dehydrogenase, between genetic lines of the wing-polymorphic cricket, Gryllus firmus, which differ in lipid biosynthesis and life history. With one exception, variation among 21 Nadp(+)-Idh genomic sequences, which spanned the entire coding sequence of the gene, was restricted to a few synonymous substitutions within and among replicate flight-capable or flightless lines. No NADP(+)-IDH electromorph variation was observed among individuals within or among lines as determined by polyacrylamide gel electrophoresis. Nor did any NADP(+)-IDH kinetic or stability parameter, such as K(M) for substrate or cofactor, k(cat), or thermal denaturation, differ between flight-capable and flightless lines. By contrast, line differences in NADP(+)-IDH-specific activity strongly covaried with transcript abundance and enzyme protein concentration. These results demonstrate that NADP(+)-IDH-specific activity differences between artificially selected lines of G. firmus are due primarily, if not exclusively, to genetic variation in regulators of NADP(+)-IDH gene expression, with no observed contribution from altered catalytic efficiency of the enzyme due to changes in amino acid sequence or posttranslational modification. Kinetic analyses indicate that in vitro differences in enzyme-specific activity between flight-capable and flightless lines likely occur in vivo. This study constitutes the most comprehensive analysis to date of the biochemical and molecular causes of naturally occurring genetic variation in enzyme activity that covaries strongly with life history.

Molecular population genetics and selection in the glycolytic pathway

In this review, I discuss the evidence for differential natural selection acting across enzymes in the glycolytic pathway in Drosophila. Across the genome, genes evolve at very different rates and possess markedly varying levels of molecular polymorphism, codon bias and expression variation. Discovering the underlying causes of this variation has been a challenge in evolutionary biology. It has been proposed that both the intrinsic properties of enzymes and their pathway position have direct effects on their molecular evolution, and with the genomic era the study of adaptation has been taken to the level of pathways and networks of genes and their products. Of special interest have been the energy-producing pathways. Using both population genetic and experimental approaches, our laboratory has been engaged in a study of molecular variation across the glycolytic pathway in Drosophila melanogaster and its close relatives. We have observed a pervasive pattern in which genes at the top of the pathway, especially around the intersection at glucose 6-phosphate, show evidence for both contemporary selection, in the form of latitudinal allele clines, and inter-specific selection, in the form of elevated levels of amino acid substitutions between species. To further explore this question, future work will require corroboration in other species, expansion into tangential pathways, and experimental work to better characterize metabolic control through the pathway and to examine the pleiotropic effects of these genes on other traits and fitness components.

Integrating evolutionary and functional approaches to infer adaptation at specific loci

Inferences about adaptation at specific loci are often exclusively based on the static analysis of DNA sequence variation. Ideally,population-genetic evidence for positive selection serves as a stepping-off point for experimental studies to elucidate the functional significance of the putatively adaptive variation. We argue that inferences about adaptation at specific loci are best achieved by integrating the indirect, retrospective insights provided by population-genetic analyses with the more direct, mechanistic insights provided by functional experiments. Integrative studies of adaptive genetic variation may sometimes be motivated by experimental insights into molecular function, which then provide the impetus to perform population genetic tests to evaluate whether the functional variation is of adaptive significance. In other cases, studies may be initiated by genome scans of DNA variation to identify candidate loci for recent adaptation. Results of such analyses can then motivate experimental efforts to test whether the identified candidate loci do in fact contribute to functional variation in some fitness-related phenotype. Functional studies can provide corroborative evidence for positive selection at particular loci, and can potentially reveal specific molecular mechanisms of adaptation.

Microevolution of intermediary metabolism: evolutionary genetics meets metabolic biochemistry

During the past decade, microevolution of intermediary metabolism has become an important new research focus at the interface between metabolic biochemistry and evolutionary genetics. Increasing recognition of the importance of integrative studies in evolutionary analysis, the rising interest in ‘evolutionary systems biology’, and the development of various ‘omics’ technologies have all contributed significantly to this developing interface. The present review primarily focuses on five prominent areas of recent research on pathway microevolution: lipid metabolism and life-history evolution; the electron transport system, hybrid breakdown and speciation; glycolysis, alcohol metabolism and population adaptation in Drosophila; chemostat selection in microorganisms; and anthocyanin pigment biosynthesis and flower color evolution. Some of these studies have provided a new perspective on important evolutionary topics that have not been investigated extensively from a biochemical perspective (hybrid breakdown, parallel evolution). Other studies have provided new data that augment previous biochemical information, resulting in a deeper understanding of evolutionary mechanisms (allozymes and biochemical adaptation to climate, life-history evolution, flower pigments and the genetics of adaptation). Finally, other studies have provided new insights into how the function or position of an enzyme in a pathway influences its evolutionary dynamics, in addition to providing powerful experimental models for investigations of network evolution. Microevolutionary studies of metabolic pathways will undoubtedly become increasingly important in the future because of the central importance of intermediary metabolism in organismal fitness, the wealth of biochemical data being provided by various omics technologies, and the increasing influence of integrative and systems perspectives in biology.

Comparative genomics of the vertebrate insulin/TOR signal transduction pathway: a network-level analysis of selective pressures

Complexity of biological function relies on large networks of interacting molecules. However, the evolutionary properties of these networks are not fully understood. It has been shown that selective pressures depend on the position of genes in the network. We have previously shown that in the Drosophila insulin/target of rapamycin (TOR) signal transduction pathway there is a correlation between the pathway position and the strength of purifying selection, with the downstream genes being most constrained. In this study, we investigated the evolutionary dynamics of this well-characterized pathway in vertebrates. More specifically, we determined the impact of natural selection on the evolution of 72 genes of this pathway. We found that in vertebrates there is a similar gradient of selective constraint in the insulin/TOR pathway to that found in Drosophila. This feature is neither the result of a polarity in the impact of positive selection nor of a series of factors affecting selective constraint levels (gene expression level and breadth, codon bias, protein length, and connectivity). We also found that pathway genes encoding physically interacting proteins tend to evolve under similar selective constraints. The results indicate that the architecture of the vertebrate insulin/TOR pathway constrains the molecular evolution of its components. Therefore, the polarity detected in Drosophila is neither specific nor incidental of this genus. Hence, although the underlying biological mechanisms remain unclear, these may be similar in both vertebrates and Drosophila.

Genotype by temperature interactions in the metabolic rate of the Glanville fritillary butterfly

Metabolic rate is a highly plastic trait. Here I examine factors that influence the metabolic rate of the Glanville fritillary butterfly (Melitaea cinxia) in pupae and resting and flying adults. Body mass and temperature had consistent positive effects on metabolic rate in pupae and resting adults but not in flying adults. There was also a consistent nonlinear effect of the time of the day, which was strongest in pupae and weakest in flying adults. Flight metabolic rate was strongly affected by an interaction between the phosphoglucose isomerase (Pgi) genotype and temperature. Over a broad range of measurement temperatures, heterozygous individuals at a single nucleotide polymorphism (SNP) in Pgi had higher peak metabolic rate in flight, but at high temperatures homozygous individuals performed better. The two genotypes did not differ in resting metabolic rate, suggesting that the heterozygotes do not pay an additional energetic cost for their higher flight capacity. Mass-independent resting and flight metabolic rates were at best weakly correlated at the individual level, and therefore, unlike in many vertebrates, resting metabolic rate does not serve as a useful surrogate of the metabolic capacity of this butterfly.

The physiology of climate change: how potentials for acclimatization and genetic adaptation will determine 'winners' and 'losers'

Physiological studies can help predict effects of climate change through determining which species currently live closest to their upper thermal tolerance limits, which physiological systems set these limits, and how species differ in acclimatization capacities for modifying their thermal tolerances. Reductionist studies at the molecular level can contribute to this analysis by revealing how much change in sequence is needed to adapt proteins to warmer temperatures--thus providing insights into potential rates of adaptive evolution--and determining how the contents of genomes--protein-coding genes and gene regulatory mechanisms--influence capacities for adapting to acute and long-term increases in temperature. Studies of congeneric invertebrates from thermally stressful rocky intertidal habitats have shown that warm-adapted congeners are most susceptible to local extinctions because their acute upper thermal limits (LT(50) values) lie near current thermal maxima and their abilities to increase thermal tolerance through acclimation are limited. Collapse of cardiac function may underlie acute and longer-term thermal limits. Local extinctions from heat death may be offset by in-migration of genetically warm-adapted conspecifics from mid-latitude 'hot spots', where midday low tides in summer select for heat tolerance. A single amino acid replacement is sufficient to adapt a protein to a new thermal range. More challenging to adaptive evolution are lesions in genomes of stenotherms like Antarctic marine ectotherms, which have lost protein-coding genes and gene regulatory mechanisms needed for coping with rising temperature. These extreme stenotherms, along with warm-adapted eurytherms living near their thermal limits, may be the major 'losers' from climate change.

The evolution of control and distribution of adaptive mutations in a metabolic pathway

In an attempt to understand whether it should be expected that some genes tend to be used disproportionately often by natural selection, we investigated two related phenomena: the evolution of flux control among enzymes in a metabolic pathway and properties of adaptive substitutions in pathway enzymes. These two phenomena are related by the principle that adaptive substitutions should occur more frequently in enzymes with greater flux control. Predicting which enzymes will be preferentially involved in adaptive evolution thus requires an evolutionary theory of flux control. We investigated the evolution of enzyme control in metabolic pathways with two models of enzyme kinetics: metabolic control theory (MCT) and Michaelis-Menten saturation kinetics (SK). Our models generate two main predictions for pathways in which reactions are moderately to highly irreversible: (1) flux control will evolve to be highly unequal among enzymes in a pathway and (2) upstream enzymes evolve a greater control coefficient then those downstream. This results in upstream enzymes fixing the majority of beneficial mutations during adaptive evolution. Once the population has reached high fitness, the trend is reversed, with the majority of neutral/slightly deleterious mutations occurring in downstream enzymes. These patterns are the result of three factors (the first of these is unique to the MCT simulations while the other two seem to be general properties of the metabolic pathways): (1) the majority of randomly selected, starting combinations of enzyme kinetic rates generate pathways that possess greater control for the upstream enzymes compared to downstream enzymes; (2) selection against large pools of intermediate substrates tends to prevent majority control by downstream enzymes; and (3) equivalent mutations in enzyme kinetic rates have the greatest effect on flux for enzymes with high levels of flux control, and these enzymes will accumulate adaptive substitutions, strengthening their control. Prediction 1 is well supported by available data on control coefficients. Data for evaluating prediction 2 are sparse but not inconsistent with this prediction.

Heritability of and strong single gene (Pgi) effects on life-history traits in the Glanville fritillary butterfly

We estimated broad-sense heritabilities (H(2)) of 13 female and seven male life-history traits of the Glanville fritillary butterfly (Melitaea cinxia) under semi-natural conditions in a large outdoor population cage. The analysis was based on full-sib families collected as young larvae in the field and reared under common garden conditions. We found significant genetic variance in female lifespan, fecundity, number of matings and host-plant preference as well as in male body mass and mobility. Apart from host-plant preference, female traits that were more strongly correlated with lifetime reproductive success (LRS; measured as total number of eggs laid) had higher H(2). LRS itself exhibited significant heritability. Host-plant preference had very high H(2), consistent with a previously reported genetically determined geographical cline in host-plant preference in the study area. Lifespan and egg hatching rate were significantly associated with a SNP in the coding region of the Pgi gene, for which there is previous evidence for balancing selection. Selection on Pgi, which furthermore shows spatial and temporal variation, may maintain genetic variance in fitness-related life-history traits. In contrast, we found no strong evidence for life-history trade-offs.

Linking genotypes to phenotypes and fitness: how mechanistic biology can inform molecular ecology

The accessibility of new genomic resources, high-throughput molecular technologies and analytical approaches such as genome scans have made finding genes contributing to fitness variation in natural populations an increasingly feasible task. Once candidate genes are identified, we argue that it is necessary to take a mechanistic approach and work up through the levels of biological organization to fully understand the impacts of genetic variation at these candidate genes. We demonstrate how this approach provides testable hypotheses about the causal links among levels of biological organization, and assists in designing relevant experiments to test the effects of genetic variation on phenotype, whole-organism performance capabilities and fitness. We review some of the research programs that have incorporated mechanistic approaches when examining naturally occurring genetic and phenotypic variation and use these examples to highlight the value of developing a comprehensive understanding of the relationship between genotype and fitness. We give suggestions to guide future research aimed at uncovering and understanding the genetic basis of adaptation and argue that further integration of mechanistic approaches will help molecular ecologists better understand the evolution of natural populations.

Transcriptional regulation of metabolism associated with the increased desiccation resistance of the cactophilic Drosophila mojavensis

In Drosophila, adaptation to xeric environments presents many challenges, greatest among them the maintenance of water balance. Drosophila mojavensis, a cactophilic species from the deserts of North America, is one of the most desiccation resistant in the genus, surviving low humidity primarily by reducing its metabolic rate. Genetic control of reduced metabolic rate, however, has yet to be elucidated. We utilized the recently sequenced genome of D. mojavensis to create an oligonucleotide microarray to pursue the identities of the genes involved in metabolic regulation during desiccation. We observed large differences in gene expression between male and female D. mojavensis as well as both quantitative and qualitative sex differences in their ability to survive xeric conditions. As expected, genes associated with metabolic regulation and carbohydrate metabolism were differentially regulated between stress treatments. Most importantly, we identified four points in central metabolism (Glyceraldehyde 3-phosphate dehydrogenase, transaldolase, alcohol dehydrogenase, and phosphoenolpyruvate carboxykinase) that indicate the potential mechanisms controlling metabolic rate reduction associated with desiccation resistance. Furthermore, a large number of genes associated with vision pathways also were differentially expressed between stress treatments, especially in females, that may underlie the initial detection of stressful environments and trigger subsequent metabolic changes.

Molecular identification of differentially regulated genes in the hydrothermal-vent species Bathymodiolus thermophilus and Paralvinella pandorae in response to temperature

BACKGROUND

Hydrothermal vents and cold seeps represent oases of life in the deep-sea environment, but are also characterized by challenging physical and chemical conditions. The effect of temperature fluctuations on vent organisms in their habitat has not been well explored, in particular at a molecular level, most gene expression studies being conducted on coastal marine species. In order to better understand the response of hydrothermal organisms to different temperature regimes, differentially expressed genes (obtained by a subtractive suppression hybridization approach) were identified in the mussel Bathymodiolus thermophilus and the annelid Paralvinella pandorae irlandei to characterize the physiological processes involved when animals are subjected to long term exposure (2 days) at two contrasting temperatures (10 degrees versus 20 degrees C), while maintained at in situ pressures. To avoid a potential effect of pressure, the experimental animals were initially thermally acclimated for 24 hours in a pressurized vessel.

RESULTS

For each species, we produced two subtractive cDNA libraries (forward and reverse) from sets of deep-sea mussels and annelids exposed together to a thermal challenge under pressure. RNA extracted from the gills, adductor muscle, mantle and foot tissue were used for B. thermophilus. For the annelid model, whole animals (small individuals) were used. For each of the four libraries, we sequenced 200 clones, resulting in 78 and 83 unique sequences in mussels and annelids (about 20% of the sequencing effort), respectively, with only half of them corresponding to known genes. Real-time PCR was used to validate differentially expressed genes identified in the corresponding libraries. Strong expression variations have been observed for some specific genes such as the intracellular hemoglobin, the nidogen protein, and Rab7 in P. pandorae, and the SPARC protein, cyclophilin, foot protein and adhesive plaque protein in B. thermophilus.

CONCLUSION

Our results indicate that mussels and worms are not responding in the same way to temperature variations. While the results obtained for the mussel B. thermophilus seem to indicate a metabolic depression (strong decrease in the level of mRNA expression of numerous genes) when temperature increased, the annelid P. pandorae mainly displayed a strong regulation of the mRNA encoding subunits and linkers of respiratory pigments and some proteins involved in membrane structure. In both cases, these regulations seem to be partly due to a possible cellular oxidative stress induced by the simulated thermal environment (10 degrees C to 20 degrees C). This work will serve as a starting point for studying the transcriptomic response of hydrothermal mussels and annelids in future experiments in response to thermal stress at various conditions of duration and temperature challenge.