The biological importance of measuring individual variation

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

Population analysis of gene expression is typically achieved by quantifying levels of mRNA; however, gene expression is also a function of protein translation and turnover. Therefore, a complete understanding of population variation in gene expression requires quantitative knowledge of protein expression within and among natural populations. We used two-dimensional fluorescence difference gel electrophoresis (2D-DIGE) to quantitatively compare expression of heart ventricle proteins among 18 individuals in three populations of the teleost fish Fundulus. Among populations, expressions between orthologous proteins and mRNAs were generally positively correlated. Additionally, similar to the pattern of cardiac mRNA expression for the same populations, we found considerable variation in protein expression both within and among populations: Of 408 protein features in 2D gels, 34% are significantly different (P < 0.01) among individuals within a population, 9% differ between populations, and 12% have a pattern of expression that suggests they have evolved by natural selection. Although similar to mRNA expression, the frequency of significant differences among populations is larger for proteins. Similar to mRNA expressions, expressions of most proteins are correlated to the expressions of many other proteins. However, the correlations among proteins are more extensive than the correlation for similar RNAs. These correlations suggest a greater coordinate regulation of protein than mRNA expression. The larger frequency of significant differences among populations and the greater frequency of correlated expression among proteins versus among RNAs suggest that the molecular mechanisms affecting protein expression enhance the differences among populations, and these regulatory steps could be a source of variation for adaptation.

The influence of gill and liver metabolism on the predicted bioconcentration of three pharmaceuticals in fish

The potential for xenobiotic compounds to bioconcentrate is typically expressed through the bioconcentration factor (BCF), which has gained increased regulatory significance over the past decade. Due to the expense of in vivo bioconcentration studies and the growing regulatory need to assess bioconcentration potential, BCF is often calculated via single-compartment models, using K(OW) as the primary input. Recent efforts to refine BCF models have focused on physiological factors, including the ability of the organism to eliminate the compound through metabolic transformation. This study looks at the ability of in vitro biotransformation assays using S9 fractions to provide an indication of metabolic potential. Given the importance of the fish gill and liver in metabolic transformation, the metabolic loss of ibuprofen, norethindrone and propranolol was measured using rainbow trout (Oncorhynchus mykiss) and channel catfish (Ictalurus punctatus) gill and liver S9 fractions. Metabolic transformation rates (k(M)) were calculated and integrated into a refined BCF model. A significant difference was noted between BCF solely based on K(OW) and BCF including k(M). These studies indicate that the inclusion of k(M) in BCF models can bring predicted bioconcentration estimates closer to in vivo values.

Comparative transcriptomics implicates mechanisms of evolved pollution tolerance in a killifish population

Wild populations of the killifish Fundulus heteroclitus resident in heavily contaminated North American Atlantic coast estuaries have recently and independently evolved dramatic, heritable, and adaptive pollution tolerance. We compared physiological and transcriptome responses to embryonic polychlorinated biphenyl (PCB) exposures between one tolerant population and a nearby sensitive population to gain insight into genomic, physiological and biochemical mechanisms of evolved tolerance in killifish, which are currently unknown. The PCB exposure concentrations at which developmental toxicity emerged, the range of developmental abnormalities exhibited, and global as well as specific gene expression patterns were profoundly different between populations. In the sensitive population, PCB exposures produced dramatic, dose-dependent toxic effects, concurrent with the alterations in the expression of many genes. For example, PCB-mediated cardiovascular system failure was associated with the altered expression of cardiomyocyte genes, consistent with sarcomere mis-assembly. In contrast, genome-wide expression was comparatively refractory to PCB induction in the tolerant population. Tolerance was associated with the global blockade of the aryl hydrocarbon receptor (AHR) signalling pathway, the key mediator of PCB toxicity, in contrast to the strong dose-dependent up-regulation of AHR pathway elements observed in the sensitive population. Altered regulation of signalling pathways that cross-talk with AHR was implicated as one candidate mechanism for the adaptive AHR signalling repression and the pollution tolerance that it affords. In addition to revealing mechanisms of PCB toxicity and tolerance, this study demonstrates the value of comparative transcriptomics to explore molecular mechanisms of stress response and evolved adaptive differences among wild populations.

Genomic approaches with natural fish populations

Natural populations v. inbred stocks provide a much richer resource for identifying the effects of nucleotide substitutions because natural populations have greater polymorphism. Additionally, natural populations offer an advantage over most common research organisms because they are subject to natural selection, and analyses of these adaptations can be used to identify biologically important changes. Among fishes, these analyses are enhanced by having a wide diversity of species (>28 000 species, more than any other group of vertebrates) living in a huge range of environments (from below freezing to > 46 degrees C, in fresh water to salinities >40 ppt.). Moreover, fishes exhibit many different life-history and reproductive strategies and have many different phenotypes and social structures. Although fishes provide numerous advantages over other vertebrate models, there is still a dearth of available genomic tools for fishes. Fishes make up approximately half of all known vertebrate species, yet <0.2% of fish species have significant genomic resources. Nonetheless, genomic approaches with fishes have provided some of the first measures of individual variation in gene expression and insights into environmental and ecological adaptations. Thus, genomic approaches with natural fish populations have the potential to revolutionize fundamental studies of diverse fish species that offer myriad ecological and evolutionary questions.

Estimating the potential for adaptation of corals to climate warming

The persistence of tropical coral reefs is threatened by rapidly increasing climate warming, causing a functional breakdown of the obligate symbiosis between corals and their algal photosymbionts (Symbiodinium) through a process known as coral bleaching. Yet the potential of the coral-algal symbiosis to genetically adapt in an evolutionary sense to warming oceans is unknown. Using a quantitative genetics approach, we estimated the proportion of the variance in thermal tolerance traits that has a genetic basis (i.e. heritability) as a proxy for their adaptive potential in the widespread Indo-Pacific reef-building coral Acropora millepora. We chose two physiologically different populations that associate respectively with one thermo-tolerant (Symbiodinium clade D) and one less tolerant symbiont type (Symbiodinium C2). In both symbiont types, pulse amplitude modulated (PAM) fluorometry and high performance liquid chromatography (HPLC) analysis revealed significant heritabilities for traits related to both photosynthesis and photoprotective pigment profile. However, quantitative real-time polymerase chain reaction (qRT-PCR) assays showed a lack of heritability in both coral host populations for their own expression of fundamental stress genes. Coral colony growth, contributed to by both symbiotic partners, displayed heritability. High heritabilities for functional key traits of algal symbionts, along with their short clonal generation time and high population sizes allow for their rapid thermal adaptation. However, the low overall heritability of coral host traits, along with the corals' long generation time, raise concern about the timely adaptation of the coral-algal symbiosis in the face of continued rapid climate warming.

Cross-species comparison of fluoxetine metabolism with fish liver microsomes

Pharmaceuticals discharged in municipal wastewater are of emerging concern because of their potential for inducing biological effects in aquatic organisms. Selective serotonin reuptake inhibitors (SSRIs), pharmaceuticals prescribed to treat chronic depression, have been detected in receiving and wastewaters. Fluoxetine is a highly prescribed model SSRI used to assess impacts of antidepressants on aquatic organisms. In this study, in vitro hepaticfluoxetine metabolism was determined in several model fish species: rainbow trout, goldfish, zebrafish and killifish. Incubation of fluoxetine with hepatic microsomes from trout pre-treated with carbamazepine showed a time-dependant loss of fluoxetine, concomitant with an increase in norfluoxetine, the major mammalian demethylated metabolite. However, fluoxetine was not well metabolized in reactions with hepatic microsomes from untreated fish. Fluoxetine loss was greater than norfluoxetine production, indicating that norfluoxetine is not the predominant fluoxetine biotransformation product in fish. Furthermore, norfluoxetine was often undetected, possibly indicating that fluoxetine demethylation is a minor metabolic pathway in fish. Inter-species differences in fluoxetine metabolism were not evident because of high intra-species variability, although killifish appeared to have the highest hepatic metabolic capacity for fluoxetine. Fluoxetine metabolism in mammals is catalyzed by cytochrome P450 (CYP) enzymes. Trout were exposed to knownCYP inducers, carbamazepine and 3-methylcholanthrene, to assess potential induction of hepatic fluoxetine metabolism. Microsomes from 3-methylcholanthrene treated fish did not induce detectable changes in fluoxetine concentrations in vitro, indicating that fish CYP1s are not involved in fluoxetine metabolism; the CYPs involved are still unclear. Identification of metabolites other than norfluoxetine warrants further investigation.

Transcriptomic variation in a coral reveals pathways of clonal organisation

A microarray study was undertaken to examine the potential for clonal gene expression variation in a branching reef building coral, Acropora millepora. The role of small-scale gradients in light and water flow was examined by comparing gene expression levels between branch elevation (tip and base) and position (centre and edge) of replicate coral colonies (n=3). Analyses of variance revealed that almost 60% of variation in gene expression was present between colonies and 34 genes were considered differentially expressed between colonies (minimum P=6.5×10(-4)). These genes are associated with energy metabolism, protein biosynthesis and cell-cell recognition representing either genotypic variation in gene expression or the effects of specific environmental conditions that affect patterns of energy acquisition, growth and pathogen encounters. Less variation was present between central and peripheral branches (7%) and only a single gene was deemed differentially expressed (P=1.493×10(-3)). The function of this gene, a phosphatidylserine decarboxylase, suggests different growth patterns between branch positions within colonies and is consistent with the usual higher growth rates on the perimeter of corymbose-like branching coral colonies such as A. millepora. Four genes were differentially expressed between the tip and base of branches (P=3.239×10(-4)) and were associated with lysosome lipase activity and fluorescence, suggesting that branch tips may encounter higher pathogen loads or levels of mechanical stress and require greater levels of photo-protection associated with higher water flow and light levels. This study therefore confirms transcriptomic variation in response to small-scale environmental gradients consistent with differential resource allocation in clonal coral colonies.

The effect of genetic and environmental variation on metabolic gene expression

What is the relationship between genetic or environmental variation and the variation in messenger RNA (mRNA) expression? To address this, microarrays were used to examine the effect of genetic and environmental variation on cardiac mRNA expression for metabolic genes in three groups of Fundulus heteroclitus: (i) individuals sampled in the field (field), (ii) field individuals acclimated for 6 months to laboratory conditions (acclimated), or (iii) individuals bred for 10 successive generations in a laboratory environment (G10). The G10 individuals have significantly less genetic variation than individuals obtained in the field and had a significantly lower variation in mRNA expression across all genes in comparison to the other two groups (P = 0.001). When examining the gene specific variation, 22 genes had variation in expression that was significantly different among groups with lower variation in G10 individuals than in acclimated individuals. Additionally, there were fewer genes with significant differences in expression among G10 individuals vs. either acclimated or field individuals: 66 genes have statistically different levels of expression vs. 107 or 97 for acclimated or field groups. Based on the permutation of the data, these differences in the number of genes with significant differences among individuals within a group are unlikely to occur by chance (P < 0.01). Surprisingly, variation in mRNA expression in field individuals is lower than in acclimated individuals. Relative to the variation among individual within a group, few genes have significant differences in expression among groups (seven, 2.3%) and none of these are different between acclimated and field individuals. The results support the concept that genetic variation affects variation in mRNA expression and also suggests that temporal environmental variation associated with estuarine environments does not increase the variation among individuals or add to the differences among groups.

Microarray analysis reveals transcriptional plasticity in the reef building coral Acropora millepora

We investigated variation in transcript abundance in the scleractinian coral, Acropora millepora, within and between populations characteristically exposed to different turbidity regimes and hence different levels of light and suspended particulate matter. We examined phenotypic plasticity by comparing levels of gene expression between source populations and following 10 days of acclimatization to a laboratory environment. Analyses of variance revealed that 0.05% of genes were differentially expressed between source populations, 1.32% following translocation into a common laboratory and 0.07% in the interaction (source population-dependent responses to translocation). Functional analyses identified an over-representation of differentially expressed genes associated with metabolism and fluorescence categories (primarily downregulated), and environmental information processing (primarily upregulated) following translocation to a lower light and turbidity environment. Such metabolic downregulation may indicate nonoxidative stress, hibernation or caloric restriction associated with the changed environmental conditions. Green fluorescent protein-related genes were the most differentially expressed and were exclusively downregulated; however, green fluorescent protein levels remained unchanged following translocation. Photophysiological responses of corals from both locations were characterized by a decline when introduced to the common laboratory environment but remained healthy (F(v)/F(m) > 0.6). Declines in total lipid content following translocation were the greatest for inshore corals, suggesting that turbid water corals have a strong reliance on heterotrophic feeding.

Variation in yolk precursor receptor mRNA expression is a key determinant of reproductive phenotype in the zebra finch (Taeniopygia guttata)

The vitellogenin/very low density lipoprotein receptor (VTG/VLDL-R), a 95 kDa protein that belongs to the low density lipoprotein receptor gene family,mediates the uptake of yolk precursors by developing follicles during oocyte growth. However, the extent to which variation in VTG/VLDL-R expression plays a role in determining inter-individual variation in reproductive phenotype(e.g. follicle or egg size) is not known. Here we show that the mRNA sequence of the zebra finch (Taeniopygia guttata) VTG/VLDL-R shows a high degree of sequence identity (92%) with chicken VTG/VLDL-R mRNA. Using quantitative real-time PCR we measured transcriptional expression of VTG/VLDL-R mRNA in various tissues, and for different stages of oocyte growth,in individual female zebra finches. VTG/VLDL-R mRNA was expressed at high levels in vitellogenic oocytes and in skeletal muscle, and was also detectable in liver, but these tissues expressed different splice variants: the short-form LR8–in oocytes and liver, and the LR8+ form in skeletal muscle. There was significant temporal variation in VTG/VLDL-R expression during follicle growth, with highest levels in ovary and a gradual decrease from pre-F3 to F1 vitellogenic follicles. Variation in ovary mRNA expression was correlated with inter-individual variation in clutch size and laying interval. Furthermore, variation in F3 follicle VTG/VLDL-R mRNA expression was correlated with inter-individual variation in egg mass and F1 follicle mass,suggesting that VTG/VLDL receptor mRNA expression is a key determinant of inter-individual variation in reproductive phenotype.

Screening for coping style increases the power of gene expression studies

Background

Individuals of many vertebrate species show different stress coping styles and these have a striking influence on how gene expression shifts in response to a variety of challenges.

Principal Findings

This is clearly illustrated by a study in which common carp displaying behavioural predictors of different coping styles (characterised by a proactive, adrenaline-based or a reactive, cortisol-based response) were subjected to inflammatory challenge and specific gene transcripts measured in individual brains. Proactive and reactive fish differed in baseline gene expression and also showed diametrically opposite responses to the challenge for 80% of the genes investigated.

Significance

Incorporating coping style as an explanatory variable can account for some the unexplained variation that is common in gene expression studies, can uncover important effects that would otherwise have passed unnoticed and greatly enhances the interpretive value of gene expression data.

Technical analysis of cDNA microarrays

BACKGROUND

There is extensive variation in gene expression among individuals within and between populations. Accurate measures of the variation in mRNA expression using microarrays can be confounded by technical variation, which includes variation in RNA isolation procedures, day of hybridization and methods used to amplify and dye label RNA for hybridization.

METHODOLOGY/PRINCIPAL FINDINGS

In this manuscript we analyze the relationship between the amount of mRNA and the fluorescent signal from the microarray hybridizations demonstrating that for a wide-range of mRNA concentrations the fluorescent signal is a linear function of the amount of mRNA. Additionally, the separate isolation, labeling or hybridization of RNA does not add significant amounts of variation in microarray measures of gene expression. However, single or double rounds of amplification for labeling do have small but significant affects on 10% of genes, but this source of technical variation is easy to avoid. To examine both technical and stochastic biological variation, mRNA expression was measured from the same five individuals over a six-week time course.

CONCLUSION

There were few, if any, meaningful differences in gene expression among time points. Thus, microarray measures using standard laboratory procedures can be precise and quantitative and are not subject to significant random biological noise.

The environmental contribution to gene expression profiles

Microarray analysis provides a bridge between the molecular genetic analysis of model organisms in laboratory settings and studies of physiology, development, and adaptation in the wild. By sampling species across a range of environments, it is possible to gain a broad picture of the genomic response to environmental perturbation. Incorporating estimates of genetic relationships into study designs will facilitate genomic analysis of environmental plasticity by aiding the identification of major regulatory loci in natural populations.

Production of different phenotypes from the same genotype in the same environment by developmental variation

The phenotype of an organism is determined by the genes, the environment and stochastic developmental events. Although recognized as a basic biological principle influencing life history, susceptibility to diseases, and probably evolution, developmental variation (DV) has been only poorly investigated due to the lack of a suitable model organism. This obstacle could be overcome by using the recently detected, robust and highly fecund parthenogenetic marbled crayfish as an experimental animal. Batch-mates of this clonal crayfish, which were shown to be isogenic by analysis of nuclear microsatellite loci, exhibited surprisingly broad ranges of variation in coloration, growth, life-span, reproduction, behaviour and number of sense organs, even when reared under identical conditions. Maximal variation was observed for the marmorated coloration, the pattern of which was unique in each of the several hundred individuals examined. Variation among identically raised batch-mates was also found with respect to fluctuating asymmetry, a traditional indicator of the epigenetic part of the phenotype, and global DNA methylation, an overall molecular marker of an animal's epigenetic state. Developmental variation was produced in all life stages, probably by reaction-diffusion-like patterning mechanisms in early development and non-linear, self-reinforcing circuitries involving behaviour and metabolism in later stages. Our data indicate that, despite being raised in the same environment, individual genotypes can map to numerous phenotypes via DV, thus generating variability among clone-mates and individuality in a parthenogenetic species. Our results further show that DV, an apparently ubiquitous phenomenon in animals and plants, can introduce components of randomness into life histories, modifying individual fitness and population dynamics. Possible perspectives of DV for evolutionary biology are discussed.

Merging methods in molecular and ecological genetics to study the adaptation of plants to anthropogenic metal-polluted sites: implications for phytoremediation

Metallophyte species that occur naturally on metal-enriched soils represent major biological resources for the improvement of phytoremediation, a benign and cost-effective technology that uses plants to clean up anthropogenic metal-polluted soils. Within the last decade, molecular genetic studies carried out on several model organisms (including Arabidopsis halleri) have considerably enhanced our understanding of metal tolerance and hyperaccumulation in plants, but the identification of the genes of interest for phytoremediation purposes remains a challenge. To meet this challenge, we propose to combine '-omics' with molecular ecology methods. Using A. halleri, we confronted molecular genetic results with: (i) within-species polymorphism and large-scale population differentiation for zinc tolerance; (ii) the demographical context (e.g. migration pattern) of the species for zinc tolerance evolution; (iii) the Quantitative Trait Loci (QTL) analysis of the genetic architecture for zinc tolerance; and (iv) the fine-scale dissection of identified QTL regions, to discuss more precisely the nature of the genes potentially involved in the adaptation to zinc-polluted soils.