Descriptive and functional characterization of variation in the Fundulus heteroclitus Ldh-B proximal promoter

Genotype - environment correlations in corals from the Great Barrier Reef

Background

Knowledge of genetic markers that are correlated to stress tolerance may improve spatial mapping of reef vulnerability and can inform restoration efforts, including the choice of genotypes for breeding and reseeding. In this manuscript we present two methods for screening transcriptome data for candidate genetic markers in two reef building corals, Acropora millepora and Pocillopora damicornis (types α and β). In A. millepora, Single Nucleotide Polymorphisms (SNPs) were pre-selected by targeting genes believed to be involved in the coral thermal stress responses. In P. damicornis (type α and β), SNPs showing varying allele frequencies between two populations from distinct environments were pre-selected. Allele frequencies at nine, five and eight of the pre-selected SNP loci were correlated against gradients of water clarity and temperature in a large number of populations along the Great Barrier Reef.

Results

A significant correlation between environmental category and SNP allele frequency was detected in up to 55% of the tested loci, which is an exceptional success rate for these types of tests. In P. damicornis, SNP allele frequencies of β-hexosaminidase and Elongation factor 1-α were significantly correlated to temperature in type α and to temperature and/or water clarity respectively in type β. Type α also showed a correlation between water clarity and SNP allele frequency in a gene of unknown function. In A. millepora, allele frequencies at five (β-gamma crystallin, Galaxin, Ubiquitin, Ligand of Numb X2 and Thioredoxin) SNP loci showed significant correlations.

Conclusions

After validation of these candidate loci through laboratory or field assessment of relative stress tolerance of colonies harbouring different alleles, it is anticipated that a proportion of these markers may represent the first coral candidate Quantitative Trait Loci for environmental stress tolerance and provide an important genetic tool that can be incorporated into spatial management decisions and restoration efforts of coral reefs. One pertinent example would be to combine spatial data of tolerant populations with genetic connectivity and thus identify high priority conservation reefs and implement targeted coral husbandry and active restoration efforts that use locally- and stress-adapted genotypes.

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.

Gene expression analysis for the identification of selection and local adaptation in fishes

In recent years, variation in gene expression has been recognized as an important component of environmental adaptation in multiple model species, including a few fish species. There is, however, still little known about the genetic basis of adaptation in gene expression resulting from variation in the aquatic environment (e.g. temperature, salinity and oxygen) and the physiological effect and costs of such differences in gene expression. This review presents and discusses progress and pitfalls of applying gene expression analyses to fishes and suggests simple frameworks to get started with gene expression analysis. It is emphasized that well-planned gene expression studies can serve as an important tool for the identification of selection in local populations of fishes, even for non-traditional model species where limited genomic information is available. Recent studies focusing on gene expression variation among natural fish populations are reviewed, highlighting the latest applications that combine genetic evidence from neutral markers and gene expression data.

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.

Classic approach revitalizes genomics: Complete characterization of a candidate gene for thermal adaptation in two coral reef fishes

Lactate dehydrogenase-B (ldh-b) encodes a metabolic enzyme (LDH-B) which plays an important role in maintaining aerobic performance and in thermal acclimation and/or adaptation of fish. As the first step in understanding the effect this enzyme has on the ability of tropical coral reef fishes to cope with thermal stress, we characterized both coding and non-coding regions of ldh-b in two congeneric perciformes, Plectropomus leopardus and Plectropomus laevis. Ldh-b was 4666 and 4539bp in length in P. leopardus and P. laevis, respectively, with coding regions comprising 1005bp in both species. We report a high level of sequence homology between the coding regions of ldh-b in these two species, with 98.1% identity of nucleotides corresponding to 100% amino acid identity between the deduced protein sequences. Comparison between non-coding (intron) regions of both species revealed the presence of several indels, despite the high level of homology observed (95.9% identity of intron nucleotides). Potential regulatory motifs and elements, including twenty-six simple sequence repeat motifs (mono-, di-, tri- and tetranucleotide) and twenty-three putative microRNA elements are identified within the introns of both species, further supporting recent demonstrations that such short motifs and elements exhibit widespread positioning throughout non-coding regions of the genome. This novel characterization of ldh-b in these two coral reef fishes allows for a wide range of future studies (e.g. analytical comparisons of ldh-b and LDH-B among different fish genera from different thermal environments and habitats).

Comparative characterization of a temperature responsive gene (lactate dehydrogenase-B, ldh-b) in two congeneric tropical fish, Lates calcarifer and Lates niloticus

The characterization of candidate loci is a critical step in obtaining insight into adaptation and acclimation of organisms. In this study of two non-model tropical (to sub-tropical) congeneric perciformes (Lates calcarifer and Lates niloticus) we characterized both coding and non-coding regions of lactate dehydrogenase-B (ldh-b), a locus which exhibits temperature-adaptive differences among temperate and sub-tropical populations of the North American killifish Fundulus heteroclitus. Ldh-b was 5,004 and 3,527 bp in length in L. calcarifer and L. niloticus, respectively, with coding regions comprising 1,005 bp in both species. A high level of sequence homology existed between species for both coding and non-coding regions of ldh-b (> 97% homology), corresponding to a 98.5% amino acid sequence homology. All six known functional sites within the encoded protein sequence (LDH-B) were conserved between the two Lates species. Ten simple sequence repeat (SSR) motifs (mono-, di-, tri- and tetranucleotide) and thirty putative microRNA elements (miRNAs) were identified within introns 1, 2, 5 and 6 of both Lates species. Five single nucleotide polymorphisms (SNPs) were also identified within miRNA containing intron regions. Such SNPs are implicated in several complex human conditions and/or diseases (as demonstrated by extensive genome-wide association studies). This novel characterization serves as a platform to further examine how non-model species may respond to changes in their native temperatures, which are expected to increase by up to 6 degrees C over the next century.

A novel hypoxia-response element in the lactate dehydrogenase-B gene of the killifish Fundulus heteroclitus

Previous studies have suggested that the lactate dehydrogenase-B gene (Ldh-B) of the Atlantic killifish, Fundulus heteroclitus, is a hypoxia-responsive gene. Here, we demonstrate that the F. heteroclitus Ldh-B promoter confers hypoxia-dependence upon reporter gene expression in transiently transfected mammalian (Hep3B) and fish (RTG-2 and RTH-149) cells in culture. Mutation and deletion analyses identified a putative hypoxia-response element (HRE) between 109 and 90 nucleotides upstream of the major start site. This HRE is characterized by the sequence 5'-GATGTG-3' spaced by 8 nucleotides from a perfect inverted repeat, and both sites are necessary for hypoxic induction of reporter gene expression in mammalian and fish cells. This HRE differs from the canonical sequence at one nucleotide position that is invariant among HREs from a wide range of hypoxia-sensitive genes. In fish cells, maximal induction of reporter gene expression driven by this HRE occurred at the lowest oxygen level tested (0.5%), took 48 h to 96 h, and was independent of glucose concentration (between 5.6 and 25 mM). Under all conditions tested, hypoxic induction of gene expression was lower in RTH-149 cells than in RTG-2, suggesting a potential defect in hypoxia signaling in RTH-149 cells. These results demonstrate that the F. heteroclitus Ldh-B promoter contains a novel HRE that is capable of driving reporter gene expression in a sequence-specific and oxygen-, time-, and cell line-dependent manner.

Evolvability of physiological and biochemical traits: evolutionary mechanisms including and beyond single-nucleotide mutation

A longstanding challenge for biologists has been to explain not just how organisms are adapted to diverse environments, but how these adaptations arise. Although natural selection is clearly sufficient to act on heritable variation, is this heritable variation sufficient to yield complex adaptations and how does this variation itself arise? Much prior focus has been on mutation of single nucleotides in genes. This process is common and can have dramatic phenotypes, but could be limited in its ability to culminate in complex adaptations for two kinds of reasons: (i) because natural selection is powerful, it can purge genetic variation, and (ii) evolutionary transition from the absence to the presence of a complex adaptation seemingly requires multiple mutations at the right place and time and in the right sequence, with each intermediate stage having increased overall fitness; this seems highly improbable. Because the networks that organisms comprise are hierarchical and redundant and have modular structure, however, single-nucleotide mutations can have large and tolerable impacts. Diverse mechanisms, collectively evolutionary capacitors, can shield genetic variation from the purgative of selection. These features can enable evolution to proceed via single-nucleotide mutation. Importantly, single-nucleotide mutation usually only modifies existing genes rather than creating new ones, and numerous other mechanisms eclipse single-nucleotide mutation in creating genetic variation. These include gene duplication (both segmental and whole-genome), lateral gene transfer, hybridization, mobile genetic elements and symbiosis. Other processes can scramble and reassemble nucleotide sequence. The mechanisms beyond single-gene mutation offer considerable promise in detailing the evolution of complex physiological and biochemical traits, and have already done so for several morphological traits.

Fundulus as the premier teleost model in environmental biology: opportunities for new insights using genomics

A strong foundation of basic and applied research documents that the estuarine fish Fundulus heteroclitus and related species are unique laboratory and field models for understanding how individuals and populations interact with their environment. In this paper we summarize an extensive body of work examining the adaptive responses of Fundulus species to environmental conditions, and describe how this research has contributed importantly to our understanding of physiology, gene regulation, toxicology, and ecological and evolutionary genetics of teleosts and other vertebrates. These explorations have reached a critical juncture at which advancement is hindered by the lack of genomic resources for these species. We suggest that a more complete genomics toolbox for F. heteroclitus and related species will permit researchers to exploit the power of this model organism to rapidly advance our understanding of fundamental biological and pathological mechanisms among vertebrates, as well as ecological strategies and evolutionary processes common to all living organisms.

Frequent gain and loss of functional transcription factor binding sites

Cis-regulatory sequences are not always conserved across species. Divergence within cis-regulatory sequences may result from the evolution of species-specific patterns of gene expression or the flexible nature of the cis-regulatory code. The identification of functional divergence in cis-regulatory sequences is therefore important for both understanding the role of gene regulation in evolution and annotating regulatory elements. We have developed an evolutionary model to detect the loss of constraint on individual transcription factor binding sites (TFBSs). We find that a significant fraction of functionally constrained binding sites have been lost in a lineage-specific manner among three closely related yeast species. Binding site loss has previously been explained by turnover, where the concurrent gain and loss of a binding site maintains gene regulation. We estimate that nearly half of all loss events cannot be explained by binding site turnover. Recreating the mutations that led to binding site loss confirms that these sequence changes affect gene expression in some cases. We also estimate that there is a high rate of binding site gain, as more than half of experimentally identified S. cerevisiae binding sites are not conserved across species. The frequent gain and loss of TFBSs implies that cis-regulatory sequences are labile and, in the absence of turnover, may contribute to species-specific patterns of gene expression.

Research in the field of molecular evolution is focused on understanding the genetic basis of functional differences between species. Protein coding sequences have traditionally been the focus of these studies, as the genetic code enables a detailed study of the strength of selection acting on amino acid sequences. However, from the earliest cross-species sequence comparisons, it was clear that protein sequences among closely related species are too similar to explain the observed phenotypic diversity. This led to the hypothesis that the evolution of gene regulation has played a key role in generating diversity between species. The availability of numerous complete genome sequences has made it possible to begin testing this hypothesis. In this work, the authors use an evolutionary model to identify functional divergence within transcription factor binding sites, the core functional elements involved in gene regulation. Applying this model to the baker's yeast, Saccharomyces cerevisiae, and its three closest relatives, the authors find that a substantial fraction of the ancestral binding sites have been lost in a species-specific manner. In some cases the loss of the binding site creates gene expression differences that may be indicative of species-specific changes in gene regulation. This work provides a useful computational framework that will allow further study of the conservation of cis-regulatory sequences and their role in molecular evolution.

Aryl hydrocarbon receptor polymorphisms and dioxin resistance in Atlantic killifish (Fundulus heteroclitus)

The aryl hydrocarbon receptor (AHR) gene encodes a ligand-activated transcription factor through which planar halogenated aromatic hydrocarbons (HAHs) such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) as well as polynuclear aromatic hydrocarbons (PAHs) cause altered gene expression and toxicity. To understand the role of AHR genetic variability in differential sensitivity to HAHs and PAHs, we are currently studying a population of the teleost Fundulus heteroclitus (Atlantic killifish) that has evolved genetic resistance to the toxic and biochemical effects of these compounds. Here, we report that the killifish AHR1 locus is highly polymorphic and that the frequencies of the major allele types differ between dioxin-sensitive and dioxin-resistant populations. Twenty-five single nucleotide polymorphisms (SNPs), nine of which are non-synonymous, were identified in the AHR1 coding sequence. Seven identified alleles were assigned to three groups, designated AHR1*1, AHR1*2 and AHR1*3. AHR1*1 alleles were under-represented in a population of dioxin- and polychlorinated biphenyl (PCB)-resistant fish from a PCB-contaminated Superfund site (New Bedford Harbor, Massachusetts, USA) compared to dioxin-sensitive fish from a less contaminated reference site (Scorton Creek, Massachusetts, USA). To determine the possible role of these AHR1 variants in differential HAH sensitivity, we expressed representative variant proteins from the two most divergent allelic groups (AHR1*1 and AHR1*3) by in-vitro transcription and translation and assessed their functional properties. AHR1*1A and AHR1*3A proteins displayed similar binding capacities and affinities for [H]TCDD. In transient transfection assays using mammalian cells, AHR1*1A and AHR1*3A exhibited similar abilities to support TCDD-dependent transactivation of a luciferase reporter gene under control of AHR-responsive enhancer elements. We discuss the possibility of other functional differences in AHR1 variants or their interaction with other killifish loci (AHR2, AHRR) that may contribute to differences in dioxin sensitivity.

Structure and sequence conservation of a putative hypoxia response element in the lactate dehydrogenase-B gene of Fundulus

Many aquatic habitats are characterized by periodic or sustained episodes of low oxygen concentration, or hypoxia, and organisms that survive in these habitats do so by utilizing a suite of behavioral, physiological and biochemical adjustments to low oxygen (1-3). In the killifish Fundulus heteroclitus, one response to prolonged exposure to hypoxia is an increase in the activity of lactate dehydrogenase-B (LDH-B), the terminal enzyme of anaerobic glycolysis, in liver tissue (4). An increase in glycolytic enzyme activity also occurs in mammalian cells during hypoxia, a process due, in part, to increased rates of gene transcription mediated by the hypoxia-inducible transcription factor, HIF-1 (5). Given that a homolog of HIF-1 has been identified in fish (6), we hypothesized that HIF might be involved in the observed up-regulation of LDH-B in F. heteroclitus. Herein, we describe the presence of DNA elements in intron 2 of the Ldh-B gene from F. heteroclitus that resemble hypoxia response elements (HRE) describedfor mammalian genes (7-10). Specifically, over a region of approximately 50 base pairs we identified two consensus HIF-1 binding sites, as well as DNA elements that may bind other transcription factors (e.g., cyclic AMP response elements; CRE). We found that these sites were perfectly conserved among geographically diverse populations of F. heteroclitus, as well as being highly conserved among multiple species in the genus Fundulus. The spacing, orientation, and sequence conservation of these putative regulatory elements suggest that they may be functionally involved in the hypoxic regulation of Ldh-B in these fish.

Environmental adaptations as windows on molecular evolution

Changes in gene regulation may play an important role in adaptive evolution, particularly during adaptation to a changing environment. However, little is known about the molecular mechanisms underlying adaptively significant variation in gene regulation. To address this question, we are using environmental adaptations in populations of a fish, Fundulus heteroclitus as a window into the molecular evolution of gene regulation. F. heteroclitus are found along the East Coast of North America, with populations distributed along a steep thermal gradient. At the extremes of the species range, populations have undergone local adaptation to their habitat temperatures. A variety of genes differ in their regulation between these populations. We have determined the mechanism responsible for changes in lactate dehydrogenase-B (Ldh-B) gene regulation. A limited number of mutations in the regulatory sequence of this gene result in changes in its expression. Both the phenotypic (increased LDH activity) and genotypic (changes in Ldh-B regulatory sequences) differences between populations have been shown to be affected by natural selection, rather than genetic drift. Therefore, even a small number of mutations within important regulatory sequences can provide evolutionarily significant variation and have an impact on environmental adaptation.

Adaptive variation in lactate dehydrogenase-B gene expression: role of a stress-responsive regulatory element

Although changes in gene regulation may play an important role in adaptive evolution, there have been few attempts to investigate the molecular mechanisms responsible for adaptively significant variation in gene expression. Here we describe the mechanism underlying an adaptive difference in the expression of the lactate dehydrogenase-B gene (Ldh-B) between northern and southern populations of the fish Fundulus heteroclitus. Ldh-B regulatory sequences from northern and southern individuals, coupled to a luciferase reporter gene, were introduced into the livers of live fish. Deletion studies indicated that sequence changes between 400 and 500 bp upstream of the transcription start site resulted in a 2-fold difference in reporter gene transcription. These sequence changes can account for the previously observed 2-fold difference in Ldh-B transcription between populations. Variation in transcription factors did not play an important role. Sequences within the functionally important region resemble a mammary tumor virus glucocorticoid responsive element (MTV-GRE) in southern alleles, whereas northern alleles differ from the consensus by 1 bp. To test the hypothesis that this element is involved in the variation between populations of F. heteroclitus, we exposed transiently transgenic fish containing Ldh-B regulatory sequence/reporter gene constructs to handling stress or injected cortisol. Both treatments increased reporter gene transcription driven by southern alleles but not northern alleles, as expected if an MTV-GRE sequence were involved. This finding suggests that sequence variation in a GRE is the cause of the adaptive differences in Ldh-B gene expression between populations and demonstrates that small changes in gene regulatory sequences can have important evolutionary consequences.