Douglas L. Crawford

Variation in gene expression within and among natural populations

Evolution may depend more strongly on variation in gene expression than on differences between variant forms of proteins. Regions of DNA that affect gene expression are highly variable, containing 0.6% polymorphic sites. These naturally occurring polymorphic nucleotides can alter in vivo transcription rates. Thus, one might expect substantial variation in gene expression between individuals. But the natural variation in mRNA expression for a large number of genes has not been measured. Here we report microarray studies addressing the variation in gene expression within and between natural populations of teleost fish of the genus Fundulus. We observed statistically significant differences in expression between individuals within the same population for approximately 18% of 907 genes. Expression typically differed by a factor of 1.5, and often more than 2.0. Differences between populations increased the variation. Much of the variation between populations was a positive function of the variation within populations and thus is most parsimoniously described as random. Some genes showed unexpected patterns of expression—changes unrelated to evolutionary distance. These data suggest that substantial natural variation exists in gene expression and that this quantitative variation is important in evolution.

Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing

We present a de novo assembly of a eukaryote transcriptome using 454 pyrosequencing data. The Glanville fritillary butterfly (Melitaea cinxia; Lepidoptera: Nymphalidae) is a prominent species in population biology but had no previous genomic data. Sequencing runs using two normalized complementary DNA collections from a genetically diverse pool of larvae, pupae, and adults yielded 608 053 expressed sequence tags (mean length = 110 nucleotides), which assembled into 48 354 contigs (sets of overlapping DNA segments) and 59 943 singletons. blast comparisons confirmed the accuracy of the sequencing and assembly, and indicated the presence of c. 9000 unique genes, along with > 6000 additional microarray‐confirmed unannotated contigs. Average depth of coverage was 6.5‐fold for the longest 4800 contigs (348–2849 bp in length), sufficient for detecting large numbers of single nucleotide polymorphisms. Oligonucleotide microarray probes designed from the assembled sequences showed highly repeatable hybridization intensity and revealed biological differences among individuals. We conclude that 454 sequencing, when performed to provide sufficient coverage depth, allows de novo transcriptome assembly and a fast, cost‐effective, and reliable method for development of functional genomic tools for nonmodel species. This development narrows the gap between approaches based on model organisms with rich genetic resources vs. species that are most tractable for ecological and evolutionary studies.

Variation within and among species in gene expression: raw material for evolution

Heritable variation in regulatory or coding regions is the raw material for evolutionary processes. The advent of microarrays has recently promoted examination of the extent of variation in gene expression within and among taxa and examination of the evolutionary processes affecting variation. This review examines these issues. We find: (i) microarray‐based measures of gene expression are precise given appropriate experimental design; (ii) there is large inter‐individual variation, which is composed of a minor nongenetic component and a large heritable component; (iii) variation among populations and species appears to be affected primarily by neutral drift and stabilizing selection, and to a lesser degree by directional selection; and (iv) neutral evolutionary divergence in gene expression becomes nonlinear with greater divergence times due to functional constraint. Evolutionary analyses of gene expression reviewed here provide unique insights into partitioning of regulatory variation in nature. However, common limitations of these studies include the tendency to assume a linear relationship between expression divergence and species divergence, and failure to test explicit hypotheses that involve the ecological context of evolutionary divergence.

Variation in tissue-specific gene expression among natural populations

BackgroundVariation in gene expression is extensive among tissues, individuals, strains, populations and species. The interactions among these sources of variation are relevant for physiological studies such as disease or toxic stress; for example, it is common for pathologies such as cancer, heart failure and metabolic disease to be associated with changes in tissue-specific gene expression or changes in metabolic gene expression. But how conserved these differences are among outbred individuals and among populations has not been well documented. To address this we examined the expression of a selected suite of 192 metabolic genes in brain, heart and liver in three populations of the teleost fish Fundulus heteroclitus using a highly replicated experimental design.ResultsHalf of the genes (48%) were differentially expressed among individuals within a population-tissue group and 76% were differentially expressed among tissues. Differences among tissues reflected well established tissue-specific metabolic requirements, suggesting that these measures of gene expression accurately reflect changes in proteins and their phenotypic effects. Remarkably, only a small subset (31%) of tissue-specific differences was consistent in all three populations.ConclusionsThese data indicate that many tissue-specific differences in gene expression are unique to one population and thus are unlikely to contribute to fundamental differences between tissue types. We suggest that those subsets of treatment-specific gene expression patterns that are conserved between taxa are most likely to be functionally related to the physiological state in question.

Fish as models for environmental genomics

Fish offer important advantages for defining the organism–environment interface and responses to natural or anthropogenic stressors. Genomic approaches using fish promise increased investigative power, and have already provided insights into the mechanisms that underlie short-term and long-term environmental adaptations. The range of fish species for which genomic resources are available is increasing, but will require significant further expansion for the optimal application of fish environmental genomics.

Natural variation in cardiac metabolism and gene expression in Fundulus heteroclitus

Individual variation in gene expression is important for evolutionary adaptation and susceptibility to diseases and pathologies. In this study, we address the functional importance of this variation by comparing cardiac metabolism to patterns of mRNA expression using microarrays. There is extensive variation in both cardiac metabolism and the expression of metabolic genes among individuals of the teleost fish Fundulus heteroclitus from natural outbred populations raised in a common environment: metabolism differed among individuals by a factor of more than 2, and expression levels of 94% of genes were significantly different (P < 0.01) between individuals in a population. This unexpectedly high variation in metabolic gene expression explains much of the variation in metabolism, suggesting that it is biologically relevant. The patterns of gene expression that are most important in explaining cardiac metabolism differ between groups of individuals. Apparently, the variation in metabolism seems to be related to different patterns of gene expression in the different groups of individuals. The magnitude of differences in gene expression in these groups is not important; large changes in expression have no greater predictive value than small changes. These data suggest that variation in physiological performance is related to the subtle variation in gene expression and that this relationship differs among individuals.

The genomic landscape of rapid repeated evolutionary adaptation to toxic pollution in wild fish

Mapping genetic adaptations to pollution Many organisms have evolved tolerance to natural and human-generated toxins. Reid et al. performed a genomic analysis of killifish, geographically separate and independent populations of which have adapted recently to severe pollution (see the Perspective by Tobler and Culumber). Sequencing multiple sensitive and resistant populations revealed signals of selective sweeps for variants that may confer tolerance to toxins, some of which were shared between resistant populations. Thus, high genetic diversity in killifish seems to allow selection to act on existing variation, driving rapid adaptation to selective forces such as pollution. Science, this issue p. 1305; see also p. 1232 Genetic diversity in Atlantic killifish has allowed for convergent evolution of pollution tolerance. Atlantic killifish populations have rapidly adapted to normally lethal levels of pollution in four urban estuaries. Through analysis of 384 whole killifish genome sequences and comparative transcriptomics in four pairs of sensitive and tolerant populations, we identify the aryl hydrocarbon receptor–based signaling pathway as a shared target of selection. This suggests evolutionary constraint on adaptive solutions to complex toxicant mixtures at each site. However, distinct molecular variants apparently contribute to adaptive pathway modification among tolerant populations. Selection also targets other toxicity-mediating genes and genes of connected signaling pathways; this indicates complex tolerance phenotypes and potentially compensatory adaptations. Molecular changes are consistent with selection on standing genetic variation. In killifish, high nucleotide diversity has likely been a crucial substrate for selective sweeps to propel rapid adaptation.

The biological importance of measuring individual variation

SUMMARY Functional genomics research using Fundulus heteroclitus has focused on variation among individuals because of the evolutionary importance and value of Fundulus in explaining the human condition (why individual humans are different and are affected differently by stress, disease and drugs). Among different populations and species of Fundulus, there are evolutionarily adaptive differences in gene expression. This natural variation in gene expression seems to affect cardiac metabolism because up to 81% of the variation in glucose utilization observed in isolated heart ventricles is related to specific patterns of gene expression. The surprising result from this research is that among different groups of individuals, the expression of mRNA from different metabolic pathways explains substrate-specific metabolism. For example, variation in oxidative phosphorylation mRNAs explains glucose metabolism for one group of individuals but expression of glucose metabolism genes explains this metabolism in a different group of individuals. This variation among individuals has important implications for studies using inbred strains: conclusions based on one individual or one strain will not necessarily reflect a generalized conclusion for a population or species. Finally, there are surprisingly strong positive and negative correlations among metabolic genes, both within and between pathways. These data suggest that measures of mRNA expression are meaningful, yet there is a complexity in how gene expression is related to physiological processes.

RNA-Seq reveals complex genetic response to deepwater horizon oil release in Fundulus grandis

BackgroundThe 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.ResultsOur 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.ConclusionRNA-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.

Phylogenetic Analysis of Thermal Acclimation of the Glycolytic Enzymes in the Genus Fundulus

Physiological acclimation that alters enzyme activity can compensate for the effect of temperature on function and may be achieved by altering enzyme concentration. This study uses phylogenetic analyses to investigate the evolutionary history of and to test several hypotheses about acclimation responses among all the glycolytic enzymes. These hypotheses are that (1) acclimation increases enzyme concentration at lower temperatures to compensate for reduced activity; (2) equilibrium enzymes tend to show acclimation responses; and (3) acclimation responses are more common in species whose populations experience either large temporal or geographical temperature variations. Using maximal activities as indices of enzyme concentration, the presence of acclimation responses in all the glycolytic enzymes in the heart ventricle was determined for five species in the teleost genus Fundulus. Three of these species are distributed along the steep thermal cline of the North American Atlantic coast, and thus these species experience both seasonal and geographical variation in temperature. The other two species are found in the Gulf of Mexico and experience seasonal variation similar to the Atlantic species but no geographical variation in temperature. Two Atlantic coast species, Fundulus heteroclitus and Fundulus majalis, have unique derived acclimation responses. No derived acclimation responses occur in the Gulf species. A conserved response in hexokinase was observed within one subgenus comprising both Atlantic and Gulf species. In F. heteroclitus, enolase responded to acclimation, and in F. majalis, aldolase, triphosphate isomerase, and lactate dehydrogenase had acclimation responses. These enzymes are equilibrium enzymes, and the concentrations of all of them increase at lower temperatures, which would compensate for the effect of temperature on enzyme activity. The compensatory changes all occur in the Atlantic species and may be a mechanism for species to expand their ranges. These data suggest that physiological acclimation is evolutionarily labile.