Population proteomics: quantitative variation within and among populations in cardiac protein expression

DIGE Analysis of Fish Tissues

Two-dimensional difference gel electrophoresis (2D-DIGE) appears to be especially useful in quantitative approaches, allowing the co-separation of proteins of control samples and proteins of treated/disease samples on the same gel, eliminating gel-to-gel variability. The principle of 2D-DIGE is to label proteins prior to isoelectric focusing and use three spectrally resolvable fluorescent dyes, allowing the independent labeling of control and experimental samples. This procedure makes it possible to reduce the number of gels in an experiment, allowing the accurate and reproducible quantification of multiple samples. 2D-DIGE has been found to be an excellent methodical tool in several areas of fish research, including environmental pollution and toxicology, the mechanisms of development and disorders, reproduction, nutrition, evolution, and ecology.

DIGE Analysis of Fish Tissues

Two-dimensional difference gel electrophoresis (2D-DIGE) appears to be especially useful in quantitative approaches, allowing the co-separation of proteins of control samples from proteins of treatment/disease samples on the same gel, eliminating gel-to-gel variability. The principle of 2D-DIGE is to label proteins prior to isoelectric focusing and use three spectrally resolvable fluorescent dyes, allowing the independent labeling of control and experimental samples. This procedure makes it possible to reduce the number of gels in an experiment, allowing the accurate and reproducible quantification of multiple samples. 2D-DIGE has been found to be an excellent methodical tool in several areas of fish research, including environmental pollution and toxicology, the mechanisms of development and disorders, reproduction, nutrition, evolution, and ecology.

Acclimation and acute temperature effects on population differences in oxidative phosphorylation

Temperature changes affect metabolism on acute, acclamatory, and evolutionary time scales. To better understand temperature's affect on metabolism at these different time scales, we quantified cardiac oxidative phosphorylation (OxPhos) in three Fundulus taxa acclimated to 12 and 28°C and measured at three acute temperatures (12, 20, and 28°C). The Fundulus taxa (northern Maine and southern Georgia F. heteroclitus, and a sister taxa, F. grandis) were used to identify evolved changes in OxPhos. Cardiac OxPhos metabolism was quantified by measuring six traits: state 3 (ADP and substrate-dependent mitochondrial respiration); E state (uncoupled mitochondrial activity); complex I, II, and IV activities; and LEAK ratio. Acute temperature affected all OxPhos traits. Acclimation only significantly affected state 3 and LEAK ratio. Populations were significantly different for state 3. In addition to direct effects, there were significant interactions between acclimation and population for complex I and between population and acute temperature for state 3. Further analyses suggest that acclimation alters the acute temperature response for state 3, E state, and complexes I and II: at the low acclimation temperature, the acute response was dampened at low assay temperatures, and at the high acclimation temperature, the acute response was dampened at high assay temperatures. Closer examination of the data also suggests that differences in state 3 respiration and complex I activity between populations were greatest between fish acclimated to low temperatures when assayed at high temperatures, suggesting that differences between the populations become more apparent at the edges of their thermal range.

Variability among individuals is generated at the gene expression level

Selection acts on individuals, specifically on their differences. To understand adaptation and responses to change therefore requires knowledge of how variation is generated and distributed across traits. Variation occurs on different biological scales, from genetic through physiological to morphological, yet it is unclear which of these carries the most variability. For example, if individual variation is mainly generated by differences in gene expression, variability should decrease progressively from coding genes to morphological traits, whereas if post-translational and epigenetic effects increase variation, the opposite should occur. To test these predictions, we compared levels of variation among individuals in various measures of gene expression, physiology (including activity), and morphology in two abundant and geographically widespread Antarctic molluscs, the clam Laternula elliptica and the limpet Nacella concinna. Direct comparisons among traits as diverse as heat shock protein QPCR assays, whole transcription profiles, respiration rates, burying rate, shell length, and ash-free dry mass were made possible through the novel application of an established metric, the Wentworth Scale. In principle, this approach could be extended to analyses of populations, communities, or even entire ecosystems. We found consistently greater variation in gene expression than morphology, with physiological measures falling in between. This suggests that variability is generated at the gene expression level. These findings have important implications for refining current biological models and predictions of how biodiversity may respond to climate change.

Statistical analysis of the individual variability of 1D protein profiles as a tool in ecology: an application to parasitoid venom

Understanding the forces that shape eco-evolutionary patterns often requires linking phenotypes to genotypes, allowing characterization of these patterns at the molecular level. DNA-based markers are less informative in this aim compared to markers associated with gene expression and, more specifically, with protein quantities. The characterization of eco-evolutionary patterns also usually requires the analysis of large sample sizes to accurately estimate interindividual variability. However, the methods used to characterize and compare protein samples are generally expensive and time-consuming, which constrains the size of the produced data sets to few individuals. We present here a method that estimates the interindividual variability of protein quantities based on a global, semi-automatic analysis of 1D electrophoretic profiles, opening the way to rapid analysis and comparison of hundreds of individuals. The main original features of the method are the in silico normalization of sample protein quantities using pictures of electrophoresis gels at different staining levels, as well as a new method of analysis of electrophoretic profiles based on a median profile. We demonstrate that this method can accurately discriminate between species and between geographically distant or close populations, based on interindividual variation in venom protein profiles from three endoparasitoid wasps of two different genera (Psyttalia concolor, Psyttalia lounsburyi and Leptopilina boulardi). Finally, we discuss the experimental designs that would benefit from the use of this method.

Stressor-induced proteome alterations in zebrafish: a meta-analysis of response patterns

Proteomics approaches are being increasingly applied in ecotoxicology on the premise that the identification of specific protein expression changes in response to a particular chemical would allow elucidation of the underlying molecular pathways leading to an adverse effect. This in turn is expected to promote the development of focused testing strategies for specific groups of toxicants. Although both gel-based and gel-free global characterization techniques provide limited proteome coverage, the conclusions regarding the cellular processes affected are still being drawn based on the few changes detected. To investigate how specific the detected responses are, we analyzed a set of studies that characterized proteome alterations induced by various physiological, chemical and biological stressors in zebrafish, a popular model organism. Our analysis highlights several proteins and protein groups, including heat shock and oxidative stress defense proteins, energy metabolism enzymes and cytoskeletal proteins, to be most frequently identified as responding to diverse stressors. In contrast, other potentially more specifically responding protein groups are detected much less frequently. Thus, zebrafish proteome responses to stress reported by different studies appear to depend mostly on the level of stress rather than on the specific stressor itself. This suggests that the most broadly used current proteomics technologies do not provide sufficient proteome coverage to allow in-depth investigation of specific mechanisms of toxicant action. We suggest that the results of any differential proteomics experiment performed with zebrafish should be interpreted keeping in mind the list of the most frequent responders that we have identified. Similar reservations should apply to any other species where proteome responses are analyzed by global proteomics methods. Careful consideration of the reliability and significance of observed changes is necessary in order not to over-interpret the experimental results and to prevent the proliferation of false positive linkages between the chemical and the cellular functions it perturbs. We further discuss the implications of the identified "top lists" of frequently responding proteins and protein families, and suggest further directions for proteomics research in ecotoxicology. Apart from improving the proteome coverage, further research should focus on defining the significance of the observed stress response patterns for organism phenotypes and on searching for common upstream regulators that can be targeted by specific assays.

Ten years of transcriptomics in wild populations: what have we learned about their ecology and evolution?

Molecular ecology has moved beyond the use of a relatively small number of markers, often noncoding, and it is now possible to use whole-genome measures of gene expression with microarrays and RNAseq (i.e. transcriptomics) to capture molecular response to environmental challenges. While transcriptome studies are shedding light on the mechanistic basis of traits as complex as personality or physiological response to catastrophic events, these approaches are still challenging because of the required technical expertise, difficulties with analysis and cost. Still, we found that in the last 10 years, 575 studies used microarrays or RNAseq in ecology. These studies broadly address three questions that reflect the progression of the field: (i) How much variation in gene expression is there and how is it structured? (ii) How do environmental stimuli affect gene expression? (iii) How does gene expression affect phenotype? We discuss technical aspects of RNAseq and microarray technology, and a framework that leverages the advantages of both. Further, we highlight future directions of research, particularly related to moving beyond correlation and the development of additional annotation resources. Measuring gene expression across an array of taxa in ecological settings promises to enrich our understanding of ecology and genome function.

Spatial variation in transcript and protein abundance of Atlantic salmon during feeding migration in the Baltic Sea

The fitness and reproductive output of fishes can be affected by environmental disturbances. In this study, transcriptomics and label-free proteomics were combined to investigate Atlantic salmon (Salmo salar) sampled from three different field locations within the Baltic Sea (Baltic Main Basin (BMB), Gulf of Finland (GoF), and Bothnian Sea (BS)) during marine migration. The expression of several stress related mRNAs and proteins of xenobiotic metabolism, oxidative stress, DNA damage, and cell death were increased in salmon from GoF compared to salmon from BMB or BS. Respiratory electron chain and ATP synthesis related gene ontology-categories were upregulated in GoF salmon, whereas those associated with RNA processing and synthesis, translation, and protein folding decreased. Differences were seen also in metabolism and immune function related gene expression. Comparisons of the transcriptomic and proteomic profiles between salmon from GoF and salmon from BMB or BS suggest environmental stressors, especially exposure to contaminants, as a main explanation for differences. Salmon feeding in GoF are thus “disturbed by hazardous substances”. The results may also be applied in evaluating the conditions of pelagic ecosystems in the different parts of Baltic Sea.

Proteome variance differences within populations of European whitefish (Coregonus lavaretus) originating from contrasting salinity environments

Variation in gene expression is an important component of the phenotypic differences observed in nature. Gene expression variance across biological groups and environmental conditions has been studied extensively and has revealed specific genes and molecular mechanisms of interest. However, little is known regarding the importance of within-population gene expression variation to environmental adaptation. To address this issue, we quantified the proteomes of individuals of European whitefish (Coregonus lavaretus) from populations that have previously been shown to have adapted during early development to freshwater and brackishwater salinity environments. Using MS-based label-free proteomics, we studied 955 proteins in eight hatch-stage fish embryos from each population that had been reared in either freshwater or brackishwater salinity conditions. By comparing the levels of within-population protein expression variance over individuals and per protein between populations, we found that fish embryos from the population less affected by salinity level had also markedly higher levels of expression variance. Gene Ontologies and molecular pathways associated with osmoregulation showed the most significant difference of within-population proteome variance between populations. Several new candidate genes for salinity adaptation were identified, emphasising the added value of combining assessments of within-population gene expression variation with standard gene expression analysis practices for better understanding the mechanisms of environmental adaptation.

BIOLOGICAL SIGNIFICANCE

We demonstrate the benefits of studying within-population gene expression variance together with more typical methods of gene expression profiling. Proteome variance differences within European whitefish populations originating from different salinity environments allowed us to identify several new candidate genes for salinity adaptation in teleost fish and generate many further hypotheses to be tested. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Assessing gene network stability and individual variability in the fathead minnow (Pimephales promelas) transcriptome

Transcriptomics is increasingly used to assess biological responses to environmental stimuli and stressors such as aquatic pollutants. However, fundamental studies characterizing individual variability in mRNA levels are lacking, which currently limits the use of transcriptomics in environmental monitoring assessments. To address individual variability in transcript abundance, we performed a meta-analysis on 231 microarrays that were conducted in the fathead minnow (FHM), a widely used toxicological model. The mean variability for gene probes was ranked from most to least variable based upon the coefficient of variation. Transcripts that were the most variable in individual tissues included NADH dehydrogenase flavoprotein 1, GTPase IMAP family member 7-like and v-set domain-containing T-cell activation inhibitor 1-like while genes encoding ribosomal proteins (rpl24 and rpl36), basic transcription factor 3, and nascent polypeptide-associated complex alpha subunit were the least variable in individuals across a range of microarray experiments. Gene networks that showed high variability (based upon the variation in expression of individual members within the network) included cell proliferation, metabolism (steroid, lipids, and glucose), cell adhesion, vascularization, and regeneration while those that showed low variability (more stability) included mRNA and rRNA processing, regulation of translational fidelity, RNA splicing, and ribosome biogenesis. Real-time PCR was conducted on a subset of genes for comparison of variability collected from the microarrays. There was a significant positive relationship between the two methods when measuring individual variability, suggesting that variability detected in microarray data can be used to guide decisions on sample sizes for measuring transcripts in real-time PCR experiments. A power analysis revealed that measuring estrogen receptor ba (esrba) requires fewer biological replicates than that of estrogen receptor bb (esrbb) in the gonad and samples sizes required to detect a 50% change for reproductive-related transcripts is between 12 and 20. Characterizing individual variability at the molecular level will prove necessary as efforts are made toward integrating molecular tools into environmental risk assessments.

RNA-Seq reveals complex genetic response to Deepwater Horizon oil release in Fundulus grandis

Background

The release of oil resulting from the blowout of the Deepwater Horizon (DH) drilling platform was one of the largest in history discharging more than 189 million gallons of oil and subject to widespread application of oil dispersants. This event impacted a wide range of ecological habitats with a complex mix of pollutants whose biological impact is still not yet fully understood. To better understand the effects on a vertebrate genome, we studied gene expression in the salt marsh minnow Fundulus grandis, which is local to the northern coast of the Gulf of Mexico and is a sister species of the ecotoxicological model Fundulus heteroclitus. To assess genomic changes, we quantified mRNA expression using high throughput sequencing technologies (RNA-Seq) in F. grandis populations in the marshes and estuaries impacted by DH oil release. This application of RNA-Seq to a non-model, wild, and ecologically significant organism is an important evaluation of the technology to quickly assess similar events in the future.

Results

Our de novo assembly of RNA-Seq data produced a large set of sequences which included many duplicates and fragments. In many cases several of these could be associated with a common reference sequence using blast to query a reference database. This reduced the set of significant genes to 1,070 down-regulated and 1,251 up-regulated genes. These genes indicate a broad and complex genomic response to DH oil exposure including the expected AHR-mediated response and CYP genes. In addition a response to hypoxic conditions and an immune response are also indicated. Several genes in the choriogenin family were down-regulated in the exposed group; a response that is consistent with AH exposure. These analyses are in agreement with oligonucleotide-based microarray analyses, and describe only a subset of significant genes with aberrant regulation in the exposed set.

Conclusion

RNA-Seq may be successfully applied to feral and extremely polymorphic organisms that do not have an underlying genome sequence assembly to address timely environmental problems. Additionally, the observed changes in a large set of transcript expression levels are indicative of a complex response to the varied petroleum components to which the fish were exposed.

The relationship between phenotypic and environmental variation: do physiological responses reduce interindividual differences?

What is the effect of a variable environment on phenotypic variation? Does the physiological response to a new environment increase or decrease the differences among individuals? We provide a speculative hypothesis suggesting that the induction of a physiological response to environmental change minimizes phenotypic differences among individuals in outbred genetically variable populations. Although this suggestion runs counter to the general idea that environmental variation induces phenotypic variation, we provide evidence that this is not always the case. One explanation for this counterintuitive hypothesis is that in a variable environment, the physiological mechanism that maintains homeostasis changes the concentrations of active transcription factors (TFs). This change in TFs reduces the effectiveness of nucleotide polymorphisms in TF binding sites and thus reduces the variation among individuals in mRNA expression and in the phenotypes affected by these mRNA transcripts. Thus, there are fewer differences among individuals in a variable environment compared with the variation observed in a constant environment. Our conjecture is that the physiological mechanisms that maintain homeostasis in response to environmental variation canalize phenotypic variation. If our hypothesis is correct, then the physiological canalization of gene expression in a variable environment hides genetic variation and thereby reduces the evolutionary costs of polymorphism. This hypothesis provides a new perspective on the mechanisms by which high levels of genetic variation can persist in real-world populations.

Comparative genomics in ecological physiology: toward a more nuanced understanding of acclimation and adaptation

Organisms that live in variable environments must adjust their physiology to compensate for environmental change. Modern functional genomics technologies offer global top-down discovery-based tools for identifying and exploring the mechanistic basis by which organisms respond physiologically to a detected change in the environment. Given that populations and species from different niches may exhibit different acclimation abilities, comparative genomic approaches may offer more nuanced understanding of acclimation responses, and provide insight into the mechanistic and genomic basis of variable acclimation. The physiological genomics literature is large and growing, as is the comparative evolutionary genomics literature. Yet, expansion of physiological genomics experiments to exploit taxonomic variation remains relatively undeveloped. Here, recent advances in the emerging field of comparative physiological genomics are considered, including examples of plants, bees and fish, and opportunities for further development are outlined particularly in the context of climate change research. Elements of robust experimental design are discussed with emphasis on the phylogenetic comparative approach. Understanding how acclimation ability is partitioned among populations and species in nature, and knowledge of the relevant genes and mechanisms, will be important for characterizing and predicting the ecological and evolutionary consequences of human-accelerated environmental change.

Proteomics in evolutionary ecology: linking the genotype with the phenotype

The study of the proteome (proteomics), which includes the dynamics of protein expression, regulation, interactions and its function, has played a less prominent role in evolutionary and ecological investigations in comparison with the study of the genome and transcriptome. There are, however, a number of arguments suggesting that this situation should change. First, the proteome is closer to the phenotype than the genome or the transcriptome, and as such may be more directly responsive to natural selection, and thus closely linked to adaptation. Second, there is evidence of a low correlation between protein and transcript expression levels across genes in many different organisms. Finally, there have been some recent important technological improvements in proteomics methods that make them feasible, practical and useful to address a wide range of evolutionary questions even in nonmodel organisms. The different proteomic methods, their limitations and problems when interpreting empirical data are described and discussed. In addition, the proteomic literature pertaining to evolutionary ecology is reviewed with examples, and potential applications of proteomics in a variety of evolutionary contexts are outlined. New proteomic research trends such as the study of posttranslational modifications and protein-protein interactions, as well as the combined use of the different -omics approaches, are discussed in relation to the development of a more functional and integrated perspective, needed for achieving a more comprehensive knowledge of evolutionary change.

Long-term survival of hydrated resting eggs from Brachionus plicatilis

Background

Several organisms display dormancy and developmental arrest at embryonic stages. Long-term survival in the dormant form is usually associated with desiccation, orthodox plant seeds and Artemia cysts being well documented examples. Several aquatic invertebrates display dormancy during embryonic development and survive for tens or even hundreds of years in a hydrated form, raising the question of whether survival in the non-desiccated form of embryonic development depends on pathways similar to those occurring in desiccation tolerant forms.

Methodology/Principal Findings

To address this question, Illumina short read sequencing was used to generate transcription profiles from the resting and amictic eggs of an aquatic invertebrate, the rotifer, Brachionus plicatilis. These two types of egg have very different life histories, with the dormant or diapausing resting eggs, the result of the sexual cycle and amictic eggs, the non-dormant products of the asexual cycle. Significant transcriptional differences were found between the two types of egg, with amictic eggs rich in genes involved in the morphological development into a juvenile rotifer. In contrast, representatives of classical “stress” proteins: a small heat shock protein, ferritin and Late Embryogenesis Abundant (LEA) proteins were identified in resting eggs. More importantly however, was the identification of transcripts for messenger ribonucleoprotein particles which stabilise RNA. These inhibit translation and provide a valuable source of useful RNAs which can be rapidly activated on the exit from dormancy. Apoptotic genes were also present. Although apoptosis is inconsistent with maintenance of prolonged dormancy, an altered apoptotic pathway has been proposed for Artemia, and this may be the case with the rotifer.

Conclusions

These data represent the first transcriptional profiling of molecular processes associated with dormancy in a non-desiccated form and indicate important similarities in the molecular pathways activated in resting eggs compared with desiccated dormant forms, specifically plant seeds and Artemia.

The effect of short-term hypoxic exposure on metabolic gene expression

The long-term effect of hypoxia is to decrease both the production and use of ATP and thus decrease the reliance on mitochondrial oxidative energy production. Yet, recent studies include more immediate affects of hypoxia on gene expression and these data suggest the maintenance of mitochondrial function. To better understand the short-term physiological response to hypoxia, we quantified metabolic mRNA expression in the heart ventricles and livers of the teleost fish Fundulus grandis exposed to partial oxygen pressure of 2.8 kPa (-13.5% air saturation).Twenty-eight individuals from a single population were exposed to hypoxia for 0, 4, 8, 12, 24, 48, and 96  hr. Liver and cardiac tissues were sampled from the same individuals at 0-48  hr. At 96  hr, only cardiac tissue was assayed. Gene expression was significantly different (ANOVA, P < 0.05) for 17 of 226 metabolic genes (7.5%) in cardiac tissue and for 20 of 256 (7.8%) metabolic genes in hepatic tissue. For the two tissues examined in this study, the maximum response occurred at different times. For cardiac tissue, using Dunnett's post hoc test, most of these significant differences occurred at 96  hr of exposure. For liver, all but one significant difference occurred at 4  hr. Surprisingly, too many (relative to random expectations) of the genes with significant increase in mRNA are involved in the oxidative phosphorylation pathway: 44% of the significant genes at 96  hr in the heart and 33% of the significant genes at 4  hr in the liver are involved in the oxidative phosphorylation pathway. These data indicate that there are tissue-specific differences in the timing of the response to hypoxia, yet both cardiac and hepatic tissues have increases in mRNA that code for enzyme in the oxidative phosphorylation pathway. If these changes in mRNA produce a similar change in protein, then these data suggest that the initial response to hypoxia involves an increase in the oxidative pathway potentially as a mechanism to maintain ATP production.