Genomic approaches with natural fish populations

Regional reef fish assemblage maps provide baseline biogeography for tropicalization monitoring

The Anthropocene rise in global temperatures is facilitating the expansion of tropical species into historically non-native subtropical locales, including coral reef fish. This redistribution of species, known as tropicalization, has serious consequences for economic development, livelihoods, food security, human health, and culture. Measuring the tropicalization of subtropical reef fish assemblages is difficult due to expansive species ranges, temporal distribution shifts with the movement of isotherms, and many dynamic density-dependent factors affecting occurrence and density. Therefore, in locales where tropical and subtropical species co-occur, detecting tropicalization changes relies on regional analyses of the relative densities and occurrence of species. This study provides a baseline for monitoring reef fish tropicalization by utilizing extensive monitoring data from a pivotal location in southeast Florida along a known transition between tropical and subtropical ecotones to define regional reef fish assemblages and use benthic habitat maps to spatially represent their zoogeography. Assemblages varied significantly by ecoregion, habitat depth, habitat type, and topographic relief. Generally, the southern assemblages had higher occurrences and densities of tropical species, whereas the northern assemblages had a higher occurrence and density of subtropical species. A total of 108 species were exclusive to regions south of the Bahamas Fracture Zone (BFZ) (South Palm Beach, Deerfield, Broward-Miami) and 35 were exclusive to the north (North Palm Beach, Martin), supporting the BFZ as a pivotal location that affects the coastal biogeographic extent of tropical marine species in eastern North America. Future tropicalization of reef fish assemblages are expected to be evident in temporal deviance of percent occurrence and/or relative species densities between baseline assemblages, where the poleward expansion of tropical species is expected to show the homogenization of assemblage regions as adjacent regions become more similar or the regional boundaries expand poleward. Ecoregions, habitat depth, habitat type, and relief should be incorporated into the stratification and analyses of reef fish surveys to statistically determine assemblage differences across the seascape, including those from tropicalization.

Diversity and genetic structure of freshwater shark Wallago attu: an emerging species of commercial interest

Pakistan has natural freshwater resources acting as a hotspot for diverse fish fauna. However, this aquatic fauna is declining at an alarming rate due to over-exploitation, habitat degradation, water pollution, climate change, and certain anthropogenic activities. The freshwater shark, Wallago attu, is a popular edible catfish inhabiting these freshwater ecosystems. Habitat degradation, overfishing, and human activities are heavily impacting the natural population of this species. So, sound knowledge about its population structure is necessary for its proper management in natural waters. The current study involves utilizing two mtDNA markers (COI, Cytb) to assess the genetic structure and differentiation among W. attu populations of Pakistani Rivers. Genetic variability analysis indicated a high haplotype (0.343 ± 0.046-0.870 ± 0.023) and low nucleotide diversity (0.0024 ± 0.012-0.0038 ± 0.018) among single and combined gene sequences, respectively. Overall, River Indus was populated with more diverse fauna of Wallago attu as compared to River Chenab and River Ravi. Population pairwise, Fst values (0.40-0.61) were found to be significantly different (p < 0.01) among three Riverine populations based upon combined gene sequences. The gene flow for the combined gene (COI + Cytb) dataset among three populations was less than 1.0. The transition/transversion bias value R (0.58) was calculated for testing of neutral evolution, and it declared low genetic polymorphism among natural riverine populations of Wallago attu. The current study's findings would be meaningful in planning the management and conservation of this economically important catfish in future.

Linking Pedigree Information to the Gene Expression Phenotype to Understand Differential Family Survival Mechanisms in Highly Fecund Fish: A Case Study in the Larviculture of Pacific Bluefin Tuna

Mass spawning in fish culture often brings about a marked variance in family size, which can cause a reduction in effective population sizes in seed production for stock enhancement. This study reports an example of combined pedigree information and gene expression phenotypes to understand differential family survival mechanisms in early stages of Pacific bluefin tuna, Thunnus orientalis, in a mass culture tank. Initially, parentage was determined using the partial mitochondrial DNA control region sequence and 11 microsatellite loci at 1, 10, 15, and 40 days post-hatch (DPH). A dramatic proportional change in the families was observed at around 15 DPH; therefore, transcriptome analysis was conducted for the 15 DPH larvae using a previously developed oligonucleotide microarray. This analysis successfully addressed the family-specific gene expression phenotypes with 5739 differentially expressed genes and highlighted the importance of expression levels of gastric-function-related genes at the developmental stage for subsequent survival. This strategy demonstrated herein can be broadly applicable to species of interest in aquaculture to comprehend the molecular mechanism of parental effects on offspring survival, which will contribute to the optimization of breeding technologies.

Population structure and genetic diversity of hatchery stocks as revealed by combined mtDNA fragment sequences in Indian major carp, Catla catla

Catla catla is the second most important Indian major carp due to high growth rate and acceptance to consumers for food value. It is widely cultured in the Indian subcontinent as monoculture or polyculture. In the present study, genetic diversity among hatchery stocks (total 218 samples of catla) collected from different geographical regions of India was examined using mtDNA fragment sequence of Cyt b (306 bp) and D loop (710 bp). High numbers (57) of population specific haplotypes were observed in the present study. The results revealed significant genetic heterogeneity for the sequence data (F_ST = 0.27546, p < .05). Analysis of molecular variance revealed significant genetic differentiation among different catla populations. The information generated in present study could be useful to develop broad genetic base populations of catla.

Discovery of Genome-Wide Microsatellite Markers in Scombridae: A Pilot Study on Albacore Tuna

Recent developments in sequencing technologies and bioinformatics analysis provide a greater amount of DNA sequencing reads at a low cost. Microsatellites are the markers of choice for a variety of population genetic studies, and high quality markers can be discovered in non-model organisms, such as tuna, with these recent developments. Here, we use a high-throughput method to isolate microsatellite markers in albacore tuna, Thunnus alalunga, based on coupling multiplex enrichment and next-generation sequencing on 454 GS-FLX Titanium pyrosequencing. The crucial minimum number of polymorphic markers to infer evolutionary and ecological processes for this species has been described for the first time. We provide 1670 microsatellite design primer pairs, and technical and molecular genetics selection resulting in 43 polymorphic microsatellite markers. On this panel, we characterized 34 random and selectively neutral markers («neutral») and 9 «non-neutral» markers. The variability of «neutral» markers was screened with 136 individuals of albacore tuna from southwest Indian Ocean (42), northwest Indian Ocean (31), South Africa (31), and southeast Atlantic Ocean (32). Power analysis demonstrated that the panel of genetic markers can be applied in diversity and population genetics studies. Global genetic diversity for albacore was high with a mean number of alleles at 16.94; observed heterozygosity 66% and expected heterozygosity 77%. The number of individuals was insufficient to provide accurate results on differentiation. Of the 9 «non-neutral» markers, 3 were linked to a sequence of known function. The one is located to a sequence having an immunity function (ThuAla-Tcell-01) and the other to a sequence having energy allocation function (ThuAla-Hki-01). These two markers were genotyped on the 136 individuals and presented different diversity levels. ThuAla-Tcell-01 has a high number of alleles (20), heterozygosity (87–90%), and assignment index. ThuAla-Hki-01 has a lower number of alleles (9), low heterozygosity (24–27%), low assignment index and significant inbreeding. Finally, the 34 «neutral» and 3 «non-neutral» microsatellites markers were tested on four economically important Scombridae species—Thunnus albacares, Thunnus thynnus, Thunnus obesus, and Acanthocybium solandri.

Spatial patterns in markers of contaminant exposure, glucose and glycogen metabolism, and immunological response in juvenile winter flounder (Pseudoplueronectes americanus)

Inshore winter flounder (Pseudoplueronectes americanus) populations in NY, USA have reached record low numbers in recent years, and recruitment into the fishery appears to be limited by survival of post-settlement juvenile fish. In order to identify cellular pathways associated with site-specific variation in condition and mortality, we examined differential mRNA expression in juvenile winter flounder collected from six different bays across a gradient in human population density and sewage inputs. Illumina sequencing of pooled samples of flounder from contrasting degraded sites and less impacted sites was used to guide our choice of targets for qPCR analysis. 253 transcripts of >100bp were differentially expressed, with 60% showing strong homology to mostly teleost sequences within the NCBI database. Based on these data, transcripts representing nine genes of interest associated with contaminant exposure, immune response and glucose and glycogen metabolism were examined by qPCR in individual flounder from each site. Statistically significant site-specific differences were observed in expression of all but one gene, although patterns in expression were complex with only one (vitellogenin), demonstrating a west to east gradient consistent with known loadings of municipal sewage effluent. Principal components analysis (PCA) identified relationships among the genes evaluated. Our data indicate that juvenile winter flounder are responding to estrogenic chemicals in more urbanized coastal bays, and suggests potential mechanistic links between immune response, contaminant exposure and energy metabolism.

Crater Lake Apoyo revisited--population genetics of an emerging species flock

The polytypic Nicaraguan Midas cichlids (Amphilophus cf. citrinellus) have been established as a model system for studying the mechanisms of speciation and patterns of diversification in allopatry and sympatry. The species assemblage in Crater Lake Apoyo has been accepted as a textbook example for sympatric speciation. Here, we present a first comprehensive data set of population genetic (mtDNA & AFLPs) proxies of species level differentiation for a representative set of individuals of all six endemic Amphilophus species occurring in Crater Lake Apoyo. AFLP genetic differentiation was partitioned into a neutral and non-neutral component based on outlier-loci detection approaches, and patterns of species divergence were explored with Bayesian clustering methods. Substantial levels of admixture between species were detected, indicating different levels of reproductive isolation between the six species. Analysis of neutral genetic variation revealed several A. zaliosus as being introgressed by an unknown contributor, hereby rendering the sympatrically evolving L. Apoyo flock polyphyletic. This is contrasted by the mtDNA analysis delivering a clear monophyly signal with Crater Lake Apoyo private haplotypes characterising all six described species, but also demonstrating different demographic histories as inferred from pairwise mismatch distributions.

Transposable elements in fish chromosomes: a study in the marine cobia species

Rachycentron canadum, a unique representative of the Rachycentridae family, has been the subject of considerable biotechnological interest due to its potential use in marine fish farming. This species has undergone extensive research concerning the location of genes and multigene families on its chromosomes. Although most of the genome of some organisms is composed of repeated DNA sequences, aspects of the origin and dispersion of these elements are still largely unknown. The physical mapping of repetitive sequences on the chromosomes of R. canadum proved to be relevant for evolutionary and applied purposes. Therefore, here, we present the mapping by fluorescence in situ hybridization of the transposable element (TE) Tol2, the non-LTR retrotransposons Rex1 and Rex3, together with the 18S and 5S rRNA genes in the chromosome of this species. The Tol2 TE, belonging to the family of hAT transposons, is homogeneously distributed in the euchromatic regions of the chromosomes but with huge colocalization with the 18S rDNA sites. The hybridization signals for Rex1 and Rex3 revealed a semi-arbitrary distribution pattern, presenting differentiated dispersion in euchromatic and heterochromatic regions. Rex1 elements are associated preferentially in heterochromatic regions, while Rex3 shows a scarce distribution in the euchromatic regions of the chromosomes. The colocalization of TEs with 18S and 5S rDNA revealed complex chromosomal regions of repetitive sequences. In addition, the nonpreferential distribution of Rex1 and Rex3 in all heterochromatic regions, as well as the preferential distribution of the Tol2 transposon associated with 18S rDNA sequences, reveals a distinct pattern of organization of TEs in the genome of this species. A heterogeneous chromosomal colonization of TEs may confer different evolutionary rates to the heterochromatic regions of this species.

Natural variation in enzyme activity of the African cichlid Pseudocrenilabrus multicolor victoriae

This study describes the metabolic capacities of the African cichlid Pseudocrenilabrus multicolor victoriae from four sites in Uganda, East Africa. Fish were captured during the dry season, from two aquatic systems in different regions (Lake Nabugabo and Mpanga River). Within the Lake Nabugabo region, individuals were sampled from Lake Kayanja (normoxic) and Lwamunda Swamp (hypoxic); within the Mpanga River system, individuals were sampled from Bunoga and Kahunge (characterized by seasonal variation in dissolved oxygen (D.O.)). Enzyme activity levels of pyruvate kinase, lactate dehydrogenase, citrate synthase, and cytochrome C oxidase were measured in four tissues: white skeletal muscle, heart, brain, and liver. Two additional enzymes were measured in the liver, malate dehydrogenase and fructose 1,6-bisphosphatase. Regional differences between enzyme activities in most tissues were evident; however, little variation was observed between two sites within a region despite differences in D.O. In general, P. multicolor from the Mpanga River system displayed greater anaerobic enzyme activity in white skeletal muscle, lower gluconeogenic enzyme activity in the liver, and an overall higher enzyme activity in the heart and brain tissues than fish from the Nabugabo region. The latter may reflect a long-term adaptation to low-oxygen conditions at the metapopulation level in the Nabugabo region.

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.

Variation in gene expression along a salinity gradient in wild populations of the euryhaline black-chinned tilapia Sarotherodon melanotheron

This study evaluated variation in expression of 11 genes within and among six wild populations of the black-chinned tilapia Sarotherodon melanotheron distributed along a salinity gradient from 0 to 100. Previous laboratory studies had shown that expression of these genes was sensitive to water salinity; the current study confirmed that a number of them also varied in expression in wild populations along the salinity gradient. Principal component analysis (PCA) first distinguished two, not mutually exclusive, sets of genes: trade-off genes that were highly expressed at one or other extreme of the salinity gradient and stress genes that were up-regulated at the two salinity extremes (i.e. a U-shaped expression pattern). The PCA clearly partitioned the populations into three groups based on their gene expression patterns and their position along the salinity gradient: a freshwater (GL; 0) population, four brackish and seawater (GB, HB, SM, SF; ranging from 20 to 50) populations and a hypersaline (SK, 100) population. Individual variation in gene expression was significantly greater within the populations at the extreme compared to intermediate salinities. These results reveal phenotypically plastic regulation of gene expression in S. melanotheron, and greater osmoregulatory and plasticity costs at extreme salinities, where fitness-related traits are known to be altered.

Karyotyping and cytogenetic mapping of Atlantic cod (Gadus morhua Linnaeus, 1758)

The Atlantic cod (Gadus morhua) is an important natural resource for northern societies and is now also considered to be a promising candidate for aquaculture. In recent years, much effort has been directed towards the development of genomic tools, and genome initiatives for Atlantic cod have been established. Despite the growing attention devoted to the Atlantic cod genomics, basic aspects of its genome structure and organization remain unknown. Thus, the present work aims to study cytogenetic features of the Atlantic cod as a contribution to the knowledge of this species' genome. The Atlantic cod displays a diploid number of 46 chromosomes, with a karyotypic formula 16 m/sm + 30 st/t. Conventional karyotyping was improved by chromosomal mapping of two classes of repetitive sequences. 18S rDNA clusters were assigned to pairs 2 and 4; small amounts of 18S rDNA clusters were occasionally detected on pair 5. These findings could not be related to the geographical origin of the specimens, but were consistent with the variability of these repeated genes in fish in general. 5S ribosomal gene clusters, apparently corresponding to a single 5S rDNA class, were detected on twelve chromosomes (pairs 11, 12, 14, 17, 20 and 21). The present update of the existing but meagre information on the karyotype of Atlantic cod, plus the first physical mapping of repetitive genes in this species herein, opens the way for an integrated approach that combines genetic and physical mapping with the assembly of the genome of this commercially important species.

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.

Biotechnologies for the management of genetic resources for food and agriculture

In recent years, the land area under agriculture has declined as also has the rate of growth in agricultural productivity while the demand for food continues to escalate. The world population now stands at 7 billion and is expected to reach 9 billion in 2045. A broad range of agricultural genetic diversity needs to be available and utilized in order to feed this growing population. Climate change is an added threat to biodiversity that will significantly impact genetic resources for food and agriculture (GRFA) and food production. There is no simple, all-encompassing solution to the challenges of increasing productivity while conserving genetic diversity. Sustainable management of GRFA requires a multipronged approach, and as outlined in the paper, biotechnologies can provide powerful tools for the management of GRFA. These tools vary in complexity from those that are relatively simple to those that are more sophisticated. Further, advances in biotechnologies are occurring at a rapid pace and provide novel opportunities for more effective and efficient management of GRFA. Biotechnology applications must be integrated with ongoing conventional breeding and development programs in order to succeed. Additionally, the generation, adaptation, and adoption of biotechnologies require a consistent level of financial and human resources and appropriate policies need to be in place. These issues were also recognized by Member States at the FAO international technical conference on Agricultural Biotechnologies for Developing Countries (ABDC-10), which took place in March 2010 in Mexico. At the end of the conference, the Member States reached a number of key conclusions, agreeing, inter alia, that developing countries should significantly increase sustained investments in capacity building and the development and use of biotechnologies to maintain the natural resource base; that effective and enabling national biotechnology policies and science-based regulatory frameworks can facilitate the development and appropriate use of biotechnologies in developing countries; and that FAO and other relevant international organizations and donors should significantly increase their efforts to support the strengthening of national capacities in the development and appropriate use of pro-poor agricultural biotechnologies.

BiMFG: BIOINFORMATICS TOOLS FOR MARINE AND FRESHWATER SPECIES

Biomolecule sequences and structures of land, air and water species are determined rapidly and the data entries are unevenly distributed for different organisms. It frequently leads to the BLAST results of homologous search containing undesirable entries from organisms living in different environments. To reduce irrelevant searching results, a separate database for comparative genomics is urgently required. A comprehensive bioinformatics tool set and an integrated database, named Bioinformatics tools for Marine and Freshwater Genomics (BiMFG), are constructed for comparative analyses among model species and underwater species. Novel matching techniques based on conserved motifs and/or secondary structure elements are designed for efficiently and effectively retrieving and aligning remote sequences through cross-species comparisons. It is especially helpful when sequences under analysis possess low similarities and unresolved structural information. In addition, the system provides core techniques of multiple sequence alignment, multiple second structure profile alignment and iteratively refined multiple structural alignments for biodiversity analysis and verification in marine and freshwater biology. The BiMFG web server is freely available for use at .

Evolutionary and functional analyses of cytochrome P4501A promoter polymorphisms in natural populations

The functional importance of variable, transcriptional regulatory sequences within and among natural populations is largely unexplored. We analysed the cytochrome P4501A (CYP1A) promoter in three populations of the minnow, Fundulus heteroclitus, because two SNPs in the promoter and first intron of CYP1A are under selection in populations adapted to pollutants. To define the importance of these SNPs, 1630 bp of the CYP1A promoter and first intron and exon were sequenced in eight individuals from three populations: a population from a polluted environment resistant to some aromatic pollutants and two flanking reference populations. CYP1A is induced by many aromatic pollutants, but in populations adapted to pollutants, CYP1A has been shown to be refractory to induction. We were interested in understanding whether variation in the CYP1A promoter explains mechanism(s) of adaptation to these aromatic pollutants. The CYP1A promoter was extremely variable (an average of 9.3% of the promoter nucleotides varied among all populations) and exhibited no fixed differences between populations. As CYP1A is poorly inducible in adapted fish, we hypothesized that CYP1A promoter regions might vary functionally between populations. Unexpectedly, in vitro analysis showed significantly greater transcription from CYP1A promoters found in the population from the polluted environment relative to promoters found in both reference populations. Thus, despite extensive variation among populations and lack of fixed differences between populations, individuals from a polluted environment have significantly enhanced promoter activity. These data demonstrate that intraspecific variation, which provides the raw material for natural selection to act on, can occur while maintaining promoter function.

Genomic approaches with natural fish populations from polluted environments

Transcriptomics and population genomics are two complementary genomic approaches that can be used to gain insight into pollutant effects in natural populations. Transcriptomics identify altered gene expression pathways, and population genomics approaches more directly target the causative genomic polymorphisms. Neither approach is restricted to a predetermined set of genes or loci. Instead, both approaches allow a broad overview of genomic processes. Transcriptomics and population genomic approaches have been used to explore genomic responses in populations of fish from polluted environments and have identified sets of candidate genes and loci that appear biologically important in response to pollution. Often differences in gene expression or loci between polluted and reference populations are not conserved among polluted populations, suggesting a biological complexity that we do not yet fully understand. As genomic approaches become less expensive with the advent of new sequencing and genotyping technologies, they will be more widely used in complementary studies. However, although these genomic approaches are immensely powerful for identifying candidate genes and loci, the challenge of determining biological mechanisms that link genotypes and phenotypes remains.

Ecologically and evolutionarily important SNPs identified in natural populations

Evolution by natural selection acts on natural populations amidst migration, gene-by-environmental interactions, constraints, and tradeoffs, which affect the rate and frequency of adaptive change. We asked how many and how rapidly loci change in populations subject to severe, recent environmental changes. To address these questions, we used genomic approaches to identify randomly selected single nucleotide polymorphisms (SNPs) with evolutionarily significant patterns in three natural populations of Fundulus heteroclitus that inhabit and have adapted to highly polluted Superfund sites. Three statistical tests identified 1.4-2.5% of SNPs that were significantly different from the neutral model in each polluted population. These nonneutral patterns in populations adapted to highly polluted environments suggest that these loci or closely linked loci are evolving by natural selection. One SNP identified in all polluted populations using all tests is in the gene for the xenobiotic metabolizing enzyme, cytochrome P4501A (CYP1A), which has been identified previously as being refractory to induction in the three highly polluted populations. Extrapolating across the genome, these data suggest that rapid evolutionary change in natural populations can involve hundreds of loci, a few of which will be shared in independent events.

Functional genomics of physiological plasticity and local adaptation in killifish

Evolutionary solutions to the physiological challenges of life in highly variable habitats can span the continuum from evolution of a cosmopolitan plastic phenotype to the evolution of locally adapted phenotypes. Killifish (Fundulus sp.) have evolved both highly plastic and locally adapted phenotypes within different selective contexts, providing a comparative system in which to explore the genomic underpinnings of physiological plasticity and adaptive variation. Importantly, extensive variation exists among populations and species for tolerance to a variety of stressors, and we exploit this variation in comparative studies to yield insights into the genomic basis of evolved phenotypic variation. Notably, species of Fundulus occupy the continuum of osmotic habitats from freshwater to marine and populations within Fundulus heteroclitus span far greater variation in pollution tolerance than across all species of fish. Here, we explore how transcriptome regulation underpins extreme physiological plasticity on osmotic shock and how genomic and transcriptomic variation is associated with locally evolved pollution tolerance. We show that F. heteroclitus quickly acclimate to extreme osmotic shock by mounting a dramatic rapid transcriptomic response including an early crisis control phase followed by a tissue remodeling phase involving many regulatory pathways. We also show that convergent evolution of locally adapted pollution tolerance involves complex patterns of gene expression and genome sequence variation, which is confounded with body-weight dependence for some genes. Similarly, exploiting the natural phenotypic variation associated with other established and emerging model organisms is likely to greatly accelerate the pace of discovery of the genomic basis of phenotypic variation.