Population Size and Site Fidelity of Fundulus heteroclitus in a Macrotidal Saltmarsh Creek

Mitochondrial Ecophysiology: Assessing the Evolutionary Forces That Shape Mitochondrial Variation

The mitonuclear species concept hypothesizes that incompatibilities between interacting gene products of the nuclear and mitochondrial genomes are a major factor establishing and maintaining species boundaries. However, most of the data available to test this concept come from studies of genetic variation in mitochondrial DNA, and clines in the mitochondrial genome across contact zones can be produced by a variety of forces. Here, we show that using a combination of population genomic analyses of the nuclear and mitochondrial genomes and studies of mitochondrial function can provide insight into the relative roles of neutral processes, adaptive evolution, and mitonuclear incompatibility in establishing and maintaining mitochondrial clines, using Atlantic killifish (Fundulus heteroclitus) as a case study. There is strong evidence for a role of secondary contact following the last glaciation in shaping a steep mitochondrial cline across a contact zone between northern and southern subspecies of killifish, but there is also evidence for a role of adaptive evolution in driving differentiation between the subspecies in a variety of traits from the level of the whole organism to the level of mitochondrial function. In addition, studies are beginning to address the potential for mitonuclear incompatibilities in admixed populations. However, population genomic studies have failed to detect evidence for a strong and pervasive influence of mitonuclear incompatibilities, and we suggest that polygenic selection may be responsible for the complex patterns observed. This case study demonstrates that multiple forces can act together in shaping mitochondrial clines, and illustrates the challenge of disentangling their relative roles.

The landscape of extreme genomic variation in the highly adaptable Atlantic killifish

Understanding and predicting the fate of populations in changing environments require knowledge about the mechanisms that support phenotypic plasticity and the adaptive value and evolutionary fate of genetic variation within populations. Atlantic killifish (Fundulus heteroclitus) exhibit extensive phenotypic plasticity that supports large population sizes in highly fluctuating estuarine environments. Populations have also evolved diverse local adaptations. To yield insights into the genomic variation that supports their adaptability, we sequenced a reference genome and 48 additional whole genomes from a wild population. Evolution of genes associated with cell cycle regulation and apoptosis is accelerated along the killifish lineage, which is likely tied to adaptations for life in highly variable estuarine environments. Genome-wide standing genetic variation, including nucleotide diversity and copy number variation, is extremely high. The highest diversity genes are those associated with immune function and olfaction, whereas genes under greatest evolutionary constraint are those associated with neurological, developmental, and cytoskeletal functions. Reduced genetic variation is detected for tight junction proteins, which in killifish regulate paracellular permeability that supports their extreme physiological flexibility. Low-diversity genes engage in more regulatory interactions than high-diversity genes, consistent with the influence of pleiotropic constraint on molecular evolution. High genetic variation is crucial for continued persistence of species given the pace of contemporary environmental change. Killifish populations harbor among the highest levels of nucleotide diversity yet reported for a vertebrate species, and thus may serve as a useful model system for studying evolutionary potential in variable and changing environments.

Evolved resistance to PCB- and PAH-induced cardiac teratogenesis, and reduced CYP1A activity in Gulf killifish (Fundulus grandis) populations from the Houston Ship Channel, Texas

The Houston Ship Channel (HSC), connecting Houston, Texas to Galveston Bay and ultimately the Gulf of Mexico, is heavily industrialized and includes several areas that have historically been identified as containing significant levels of mercury, dioxins, furans, polychlorinated biphenyls (PCBs), and polycyclic aromatic hydrocarbons (PAHs). Gulf killifish, Fundulus grandis, inhabit this entire estuarine system, including the most contaminated areas. F. grandis is the sister species of the well-established estuarine model organism Fundulus heteroclitus, for which heritable resistance to both PCB and PAH toxicity has been documented in several populations. F. grandis collected from two Superfund sites on the HSC and from a reference population were used to establish breeding colonies. F1 embryos from HSC populations were approximately 1000-fold more resistant to PCB126- and 2-5-fold more resistant to coal tar-induced cardiovascular teratogenesis, relative to embryos from the reference population. Reciprocal crosses between reference and contaminated populations exhibit an intermediate level of resistance, confirming that observed protection is genetic and biparentally inherited. Ethoxyresorufin-O-deethylase (EROD) data confirm a reduction in basal and induced cytochrome P4501A (CYP1A) activity in resistant populations of F. grandis. This result is consistent with responses previously described for resistant populations of F. heteroclitus, specifically a recalcitrant aryl hydrocarbon receptor (AHR) pathway. The decreased levels of cardiovascular teratogenesis, and decrease in CYP1A inducibility in response to PCB126 and a PAH mixture, suggest that HSC F. grandis populations have adapted to chronic contaminants exposures via a mechanism similar to that previously described for F. heteroclitus. To the best of our knowledge, this is the first documentation of evolved pollution resistance in F. grandis. Additionally, the mechanistic similarities between the population adaptation observed in this study and previous work in F. heteroclitus suggest that genetic variation predating the evolutionary divergence of these two species may best explain the apparent rapid parallel evolution of pollution resistance in genetically and geographically distinct species and populations.

Salinity- and population-dependent genome regulatory response during osmotic acclimation in the killifish (Fundulus heteroclitus) gill

The killifish Fundulus heteroclitus is abundant in osmotically dynamic estuaries and it can quickly adjust to extremes in environmental salinity. We performed a comparative osmotic challenge experiment to track the transcriptomic and physiological responses to two salinities throughout a time course of acclimation, and to explore the genome regulatory mechanisms that enable extreme osmotic acclimation. One southern and one northern coastal population, known to differ in their tolerance to hypo-osmotic exposure, were used as our comparative model. Both populations could maintain osmotic homeostasis when transferred from 32 to 0.4 p.p.t., but diverged in their compensatory abilities when challenged down to 0.1 p.p.t., in parallel with divergent transformation of gill morphology. Genes involved in cell volume regulation, nucleosome maintenance, ion transport, energetics, mitochondrion function, transcriptional regulation and apoptosis showed population- and salinity-dependent patterns of expression during acclimation. Network analysis confirmed the role of cytokine and kinase signaling pathways in coordinating the genome regulatory response to osmotic challenge, and also posited the importance of signaling coordinated through the transcription factor HNF-4α. These genome responses support hypotheses of which regulatory mechanisms are particularly relevant for enabling extreme physiological flexibility.

Common mechanism underlies repeated evolution of extreme pollution tolerance

Human alterations to the environment can exert strong evolutionary pressures, yet contemporary adaptation to human-mediated stressors is rarely documented in wildlife populations. A common-garden experimental design was coupled with comparative transcriptomics to discover evolved mechanisms enabling three populations of killifish resident in urban estuaries to survive normally lethal pollution exposure during development, and to test whether mechanisms are unique or common across populations. We show that killifish populations from these polluted sites have independently converged on a common adaptive mechanism, despite variation in contaminant profiles among sites. These populations are united by a similarly profound desensitization of aryl-hydrocarbon receptor-mediated transcriptional activation, which is associated with extreme tolerance to the lethal effects of toxic dioxin-like pollutants. The rapid, repeated, heritable and convergent nature of evolved tolerance suggests that ancestral killifish populations harboured genotypes that enabled adaptation to twentieth-century industrial pollutants.

Embryonic gene expression among pollutant resistant and sensitive Fundulus heteroclitus populations

Changes in gene expression, coupled with biochemical, physiological, and behavioral alterations, play a critical role in adaptation to environmental stress. Our goal was to explore ways natural populations may have adapted to local, polluted environments. We took advantage of natural populations of Fundulus heteroclitus, one of the few studied fish species in North America that has established resistant populations in highly contaminated urban estuaries. We analyzed morphology, physiology, and gene expression of developing F. heteroclitus embryos during late organogenesis (stage 31); these embryos were from both resistant and sensitive populations and were raised in a common, unpolluted environment. While cardiac heart rates show significant differences between embryos of parents from clean and heavily contaminated Superfund sites, time-to-stage, embryo morphology, and gene expression profile analyses do not differ significantly between untreated embryos from resistant and sensitive populations. Further evaluation that includes tissue-specific approaches in gene expression analysis and larger sample sizes may be necessary to highlight important phenotypes associated with mechanisms of sensitivity and resistance among natural F. heteroclitus embryo populations. Alternatively, population differences may be masked by developmental canalization, and biologically important differences between sensitive and resistant embryos may only manifest with exposure (e.g., be dependent on gene by environment interactions).

Relative influences of historical and contemporary forces shaping the distribution of genetic variation in the Atlantic killifish, Fundulus heteroclitus

A major goal of population genetics research is to identify the relative influences of historical and contemporary processes that serve to structure genetic variation. Most population genetic models assume that populations exist in a state of migration-drift equilibrium. However, in the past this assumption has rarely been verified, and is likely rarely achieved in natural populations. We assessed the equilibrium status at both local and regional scales of the Atlantic killifish, Fundulus heteroclitus. This species is a model organism for the study of adaptive clinal variation, but has also experienced a complicated history of range expansion and secondary contact following allopatric divergence, potentially obscuring the influence of contemporary evolutionary processes. Presumptively neutral genetic markers (microsatellites) demonstrated zones of secondary intergradation among coastal populations centred around northern New Jersey and the Chesapeake Bay region. Analysis of genetic variation indicated isolation by distance among some populations and provided supporting evidence that the Delaware Bay, but not the Chesapeake Bay, has acted as a barrier to dispersal among coastal populations. Bayesian estimates indicated large effective population sizes and low migration rates, and were in good agreement with empirically derived estimates of population and neighbourhood size from mark-recapture studies. These data indicate that populations are not in migration-drift equilibrium at a regional scale, and suggest that contributing factors include large population size combined with relatively low migration rates. These conditions should be considered when interpreting the evolutionary significance of the distribution of genetic variation among F. heteroclitus populations.

Genetic diversity and structure of an estuarine fish (Fundulus heteroclitus) indigenous to sites associated with a highly contaminated urban harbor

Intense selection on isolated populations can cause loss of genetic diversity, which if persistent, reduces adaptive potential and increases extinction probability. Phenotypic evidence of inherited tolerance suggests that polychlorinated biphenyls (PCBs), have acted as strong selective agents on populations of a non-migratory fish, Fundulus heteroclitus, indigenous to heavily contaminated sites. To evaluate population genetic structure and test for effects of intense, multi-generational PCB contamination on genetic diversity, we used AFLP analysis on fish collected from six sites along the east coast of North America that varied widely in PCB contamination. The sites included a heavily contaminated urban harbor (New Bedford, MA), an adjacent moderately contaminated sub-estuary (Buzzards Bay, MA), and an uncontaminated estuary 60 km away (Narragansett, RI). AFLP markers distinguished populations at moderate and small scales, suggesting genetic differentiation at distances of 2 km or less. Genetic diversity did not differ across the study sites. Genome-wide diversity may have been preserved because of large effective population sizes and/or because the mechanism for genetic adaptation to these contaminants affected only a small number of loci. Alternatively, loss in diversity may have been restored with moderate levels of migration and relatively short generation time for this species.

Microsatellite analysis of the phylogeography, Pleistocene history and secondary contact hypotheses for the killifish, Fundulus heteroclitus

The mummichog, Fundulus heteroclitus, exhibits extensive latitudinal clinal variation in a number of physiological and biochemical traits, coupled with phylogeographical patterns at mitochondrial and nuclear DNA loci that suggest a complicated history of spatially variable selection and secondary intergradation. This species continues to serve as a model for understanding local and regional adaptation to variable environments. Resolving the influences of historical processes on the distribution of genetic variation within and among extant populations of F. heteroclitus is crucial to a better understanding of how populations evolve in the context of contemporary environments. In this study, we analysed geographical patterns of genetic variation at eight microsatellite loci among 15 populations of F. heteroclitus distributed throughout the North American range of the species from Nova Scotia to Georgia. Genetic variation in Northern populations was lower than in Southern populations and was strongly correlated with latitude throughout the species range. The most common Northern alleles at all eight loci exhibited concordant latitudinal clinal patterns, and the existence of an abrupt transition zone in allele frequencies between Northern and Southern populations was similar to that observed for mitochondrial DNA and allozyme loci. A significant pattern of isolation by distance was observed both within and between northern and southern regions. This pattern was unexpected, particularly for northern populations, given the recent colonization history of post-Pleistocene habitats, and was inconsistent with either a recent northward population expansion or a geographically restricted northern Pleistocene refugium. The data provided no evidence for recent population bottlenecks, and estimates of historical effective population sizes suggest that post-Pleistocene populations have been large throughout the species distribution. These results suggest that F. heteroclitus was broadly distributed throughout most of its current range during the last glacial event and that the abrupt transition in allele frequencies that separate Northern and Southern populations may reflect regional disequilibrium conditions associated with the post-Pleistocene colonization history of habitats in that region.