Evidence of prezygotic isolation, but not assortative mating, between locally adapted populations of Fundulus heteroclitus across a salinity gradient.

Selection along environmental gradients can drive reproductive isolation and speciation. Among fishes, salinity is a major factor limiting species distributions, and despite its importance in generating species diversity, speciation events between marine and freshwater are rare. Here, we tested for mechanisms of reproductive isolation between locally adapted freshwater and brackish water-native populations of killifish, Fundulus heteroclitus, from either side of a hybrid zone along a salinity gradient. There was evidence for pre-zygotic endogenous reproductive isolation with reduced fertilization success between crosses of freshwater-native males and brackish water-native females. Exogenous pre-zygotic isolation was also present where females had highest fertilization in their native salinity. We used a replicated mass spawning design to test for mate choice in both brackish and fresh water. After genotyping 187 parents and 2523 offspring at 2347 SNPs across the genome, 85% of offspring were successfully assign to their parents. However, no reinforcing mate choice was observed. These results therefore demonstrate emerging, yet limited, reproductive isolation and incipient speciation across a marine to freshwater salinity gradient and suggest that both endogenous and exogenous mechanisms, but not assortative mating, contribute to divergence.

Evolved differences in energy metabolism and growth dictate the impacts of ocean acidification on abalone aquaculture.

Ocean acidification (OA) poses a major threat to marine ecosystems and shellfish aquaculture. A promising mitigation strategy is the identification and breeding of shellfish varieties exhibiting resilience to acidification stress. We experimentally compared the effects of OA on two populations of red abalone (Haliotis rufescens), a marine mollusc important to fisheries and global aquaculture. Results from our experiments simulating captive aquaculture conditions demonstrated that abalone sourced from a strong upwelling region were tolerant of ongoing OA, whereas a captive-raised population sourced from a region of weaker upwelling exhibited significant mortality and vulnerability to OA. This difference was linked to population-specific variation in the maternal provisioning of lipids to offspring, with a positive correlation between lipid concentrations and survival under OA. This relationship also persisted in experiments on second-generation animals, and larval lipid consumption rates varied among paternal crosses, which is consistent with the presence of genetic variation for physiological traits relevant for OA survival. Across experimental trials, growth rates differed among family lineages, and the highest mortality under OA occurred in the fastest growing crosses. Identifying traits that convey resilience to OA is critical to the continued success of abalone and other shellfish production, and these mitigation efforts should be incorporated into breeding programs for commercial and restoration aquaculture.

A draft reference genome of the red abalone, Haliotis rufescens, for conservation genomics.

Red abalone, Haliotis rufescens, are herbivorous marine gastropods that primarily feed on kelp. They are the largest and longest-lived of abalone species with a range distribution in North America from central Oregon, United States, to Baja California, MEX. Recently, red abalone have been in decline as a consequence of overharvesting, disease, and climate change, resulting in the closure of the commercial fishery in the 1990s and the recreational fishery in 2018. Protecting this ecologically and economically important species requires an understanding of their current population dynamics and connectivity. Here, we present a new red abalone reference genome as part of the California Conservation Genomics Project (CCGP). Following the CCGP genome strategy, we used Pacific Biosciences HiFi long reads and Dovetail Omni-C data to generate a scaffold-level assembly. The assembly comprises 616 scaffolds for a total size of 1.3 Gb, a scaffold N50 of 45.7 Mb, and a BUSCO complete score of 97.3%. This genome represents a significant improvement over a previous assembly and will serve as a powerful tool for investigating seascape genomic diversity, local adaptation to temperature and ocean acidification, and informing management strategies.

Genomic and physiological mechanisms underlying skin plasticity during water to air transition in an amphibious fish.

The terrestrial radiation of vertebrates required changes in skin that resolved the dual demands of maintaining a mechanical and physiological barrier while also facilitating ion and gas transport. Using the amphibious killifish Kryptolebias marmoratus, we found that transcriptional regulation of skin morphogenesis was quickly activated upon air exposure (1 h). Rapid regulation of cell-cell adhesion complexes and pathways that regulate stratum corneum formation was consistent with barrier function and mechanical reinforcement. Unique blood vessel architecture and regulation of angiogenesis likely supported cutaneous respiration. Differences in ionoregulatory transcripts and ionocyte morphology were correlated with differences in salinity acclimation and resilience to air exposure. Evolutionary analyses reinforced the adaptive importance of these mechanisms. We conclude that rapid plasticity of barrier, respiratory and ionoregulatory functions in skin evolved to support the amphibious lifestyle of K. marmoratus; similar processes may have facilitated the terrestrial radiation of other contemporary and ancient fishes.

Religion in the Age of Social Distancing: How COVID-19 Presents New Directions for Research.

In this brief note written during a global pandemic, we consider some of the important ways this historical moment is altering the religious landscape, aiming our investigative lens at how religious institutions, congregations, and individuals are affected by the social changes produced by COVID-19. This unprecedented time prompts scholars of religion to reflect on how to strategically approach the study of religion in the time of "social distancing," as well as moving forward. Particularly important considerations include developing heuristic, innovative approaches for revealing ongoing changes to religion, as well as how religion continues to structure social life across a wide range of contexts, from the most intimate and personal to the most public and global. Although our note can only be indicative rather than exhaustive, we do suggest that the initial groundwork for reconsiderations might productively focus on several key analytical themes, including: Epidemiology, Ideology, Religious Practice, Religious Organizations and Institutions, as well as Epistemology and Methodology. In offering these considerations as a starting point, we remain aware (and hopeful) that inventive and unanticipated approaches will also emerge.

Integrative Population and Physiological Genomics Reveals Mechanisms of Adaptation in Killifish.

Adaptive divergence between marine and freshwater (FW) environments is important in generating phyletic diversity within fishes, but the genetic basis of this process remains poorly understood. Genome selection scans can identify adaptive loci, but incomplete knowledge of genotype-phenotype connections makes interpreting their significance difficult. In contrast, association mapping (genome-wide association mapping [GWAS], random forest [RF] analyses) links genotype to phenotype, but offer limited insight into the evolutionary forces shaping variation. Here, we combined GWAS, RF, and selection scans to identify loci important in adaptation to FW environments. We utilized FW-native and brackish water (BW)-native populations of Atlantic killifish (Fundulus heteroclitus) as well as a naturally admixed population between the two. We measured morphology and multiple physiological traits that differ between populations and may contribute to osmotic adaptation (salinity tolerance, hypoxia tolerance, metabolic rate, body shape) and used a reduced representation approach for genome-wide genotyping. Our results show patterns of population divergence in physiological capabilities that are consistent with local adaptation. Population genomic scans between BW-native and FW-native populations identified genomic regions evolving by natural selection, whereas association mapping revealed loci that contribute to variation for each trait. There was substantial overlap in the genomic regions putatively under selection and loci associated with phenotypic traits, particularly for salinity tolerance, suggesting that these regions and genes are important for adaptive divergence between BW and FW environments. Together, these data provide insight into the mechanisms that enable diversification of fishes across osmotic boundaries.

Mitochondria, sex and variation in routine metabolic rate.

Variation in the metabolic costs associated with organismal maintenance may play a key role in determining fitness, and thus these differences among individuals are likely to be subject to natural selection. Although the evolvability of maintenance metabolism depends on its underlying genetic architecture, relatively little is known about the nature of genetic variation that underlies this trait. To address this, we measured variation in routine metabolic rate (MO2routine ), an index of maintenance metabolism, within and among three populations of Atlantic killifish, Fundulus heteroclitus, including a population from a region of genetic admixture between two subspecies. Polygenic association tests among individuals from the admixed population identified 54 single nucleotide polymorphisms (SNPs) that were associated with MO2routine , and these SNPs accounted for 43% of interindividual variation in this trait. However, genetic associations with MO2routine involved different SNPs if females and males were analysed separately, and there was a sex-dependent effect of mitochondrial genotype on variation in routine metabolism. These results imply that there are sex-specific genetic mechanisms, and potential mitonuclear interactions, that underlie variation in MO2routine . Additionally, there was evidence for epistatic interactions between 17% of the possible pairs of trait-associated SNPs, suggesting that epistatic effects on MO2routine are common. These data demonstrate not only that phenotypic variation in this ecologically important trait has a polygenic basis with considerable epistasis among loci, but also that these underlying genetic mechanisms, and particularly the role of mitochondrial genotype, may be sex-specific.

Phenotypic flexibility in respiratory traits is associated with improved aerial respiration in an amphibious fish out of water.

Amphibious fishes have evolved multiple adaptive strategies for respiring out of water, but there has been less focus on reversible plasticity. We tested the hypothesis that when amphibious fishes leave water, enhanced respiratory performance on land is the result of rapid functional phenotypic flexibility of respiratory traits. We acclimated four isogenic strains of Kryptolebias marmoratus to air for 0, 1, 3 or 7 days. We compared respiratory performance out of water with traits linked to the O2 cascade. Aerial O2 consumption rate was measured over a step-wise decrease in O2 levels. There were significant differences between strains, but time out of water had the largest impact on measured parameters. Kryptolebiasmarmoratus had improved respiratory performance [lower aerial critical oxygen tension (Pcrit), higher regulation index (RI)] after only 1 day of air exposure, and these changes were strongly associated with the change in hematocrit and dorsal cutaneous angiogenesis. Additionally, we found that 1 h of air exposure induced the expression of four angiogenesis-associated genes - vegfa, angpt2, pecam-1 and efna1 - in the skin. After 7 days in air, respiratory traits were not significantly linked to the variation in either aerial Pcrit or RI. Overall, our data indicate that there are two phases involved in the enhancement of aerial respiration: an initial rapid response (1 day) and a delayed response (7 days). We found evidence for the hypothesis that respiratory performance on land in amphibious fishes is the result of rapid flexibility in both O2 uptake and O2 carrying capacity.

Contaminant Exposure Linked to Cellular and Endocrine Biomarkers in Southern California Bottlenose Dolphins.

Cetaceans in the Southern California Bight (SCB) are exposed to high levels of halogenated organic contaminants (HOCs), which have previously been linked to impaired reproductive health and immune responses. We used a combination of molecular tools to examine the potential physiological impacts of HOC exposure in two bottlenose dolphin ( Tursiops truncatus) ecotypes in the SCB. We quantified 25 HOCs in the blubber of 22 biopsies collected from males between 2012 and 2016. We then analyzed genome-wide gene expression in skin using RNA-sequencing and measured blubber testosterone to compare HOC exposure with cellular and endocrine biomarkers. We found high levels of HOCs in both ecotypes with significantly higher total polychlorinated biphenyls (PCBs), polybrominated diphenyl ethers (PBDEs), tris(4-chlorophenyl)methanol (TCPMOH), and chlordane-related compounds in the coastal ecotype versus the offshore ecotype. We found evidence of PBDE bioaccumulation in both ecotypes, however, the pattern of bioaccumulation or endocrine disruption for other HOCs was different between the ecotypes, suggesting potential endocrine disruption in the coastal ecotype. We also observed correlations between HOCs and gene coexpression networks enriched for xenobiotic metabolism, hormone metabolism, and immune response that could indicate cellular effects from HOC exposure. By integrating measurements of HOC load with both transcriptome profiling and endocrine biomarkers, our approach provides insight into HOC exposure and potential impacts on wild cetacean health in southern California.

Phylogenomic analysis of Fundulidae (Teleostei: Cyprinodotiformes) using RNA-sequencing data.

Fishes of the New World cyprinodontiform family Fundulidae display a wide variety of tolerance to environmental conditions, making them a valuable model system for comparative, evolutionary, and environmental studies. Despite numerous attempts to resolve the phylogenetic relationships of family Fundulidae, the basal structure of the phylogeny remains unresolved. The lack of a robust and fully resolved phylogeny for family Fundulidae and its most speciose genus Fundulus is an impediment to future research. This study utilized novel RNA-sequencing data for phylogenetic inference among16 members of Fundulidae to better refine the basal nodes of the family and confront long-standing questions regarding (1) the monophyletic status of genus Fundulus, and validity of the Lucania and recently synonymized Adinia genera; (2) the relationship of the west coast endemic Fundulus parvipinnis and F. lima to other Fundulus species; and (3) the validity of subgeneric classifications. In addition, previously published nuclear gene sequences for 32 Fundulidae species were re-analyzed in combination with novel RNA-sequencing data. Maximum likelihood and Bayesian analyses generated identical phylogenies with strong statistical support at nearly all nodes, demonstrating the utility of RNA-sequencing data in constructing robust phylogenies not achievable by previous methods. While many past hypothesized evolutionary relationships for Fundulidae were reinforced, several alternative relationships are hypothesized at basal nodes resulting in a re-analysis of the deeper structure of family Fundulidae. These results reveal family Fundulidae as a paraphyletic grouping of members of genus Fundulus and Lucania and supports the previous synonymy of genus Adinia with genus Fundulus.

Tolerance traits related to climate change resilience are independent and polygenic.

The resilience of organisms to climate change through adaptive evolution is dependent on the extent of genetically based variation in key phenotypic traits and the nature of genetic associations between them. For aquatic animals, upper thermal tolerance and hypoxia tolerance are likely to be a important determinants of sensitivity to climate change. To determine the genetic basis of these traits and to detect associations between them, we compared naturally occurring populations of two subspecies of Atlantic killifish, Fundulus heteroclitus, that differ in both thermal and hypoxia tolerance. Multilocus association mapping demonstrated that 47 and 35 single nucleotide polymorphisms (SNPs) explained 43.4% and 51.9% of variation in thermal and hypoxia tolerance, respectively, suggesting that genetic mechanisms underlie a substantial proportion of variation in each trait. However, no explanatory SNPs were shared between traits, and upper thermal tolerance varied approximately linearly with latitude, whereas hypoxia tolerance exhibited a steep phenotypic break across the contact zone between the subspecies. These results suggest that upper thermal tolerance and hypoxia tolerance are neither phenotypically correlated nor genetically associated, and thus that rates of adaptive change in these traits can be independently fine-tuned by natural selection. This modularity of important traits can underpin the evolvability of organisms to complex future environmental change.

Comprehensive Screening Links Halogenated Organic Compounds with Testosterone Levels in Male Delphinus delphis from the Southern California Bight.

While environmental pollutants have been associated with changes in endocrine health in cetaceans, efforts to link contaminant exposure with hormones have largely been limited to a list of known, targeted contaminants, overlooking minimally characterized or unknown compounds of emerging concern. To address this gap, we analyzed a suite of potential endocrine disrupting halogenated organic compounds (HOCs) in blubber from 16 male short-beaked common dolphins ( Delphinus delphis) with known maturity status collected from fishery bycatch in the Southern California Bight. We employed a suspect screening mass spectrometry-based method to investigate a wide range of HOCs that were previously observed in cetaceans from the same region. Potential endocrine effects were assessed through the measurement of blubber testosterone. We detected 167 HOCs, including 81 with known anthropogenic sources, 49 of unknown origin, and 37 with known natural sources. The sum of 11 anthropogenic and 4 unknown HOC classes were negatively correlated with blubber testosterone. Evidence suggests that elevated anthropogenic HOC load contributes to impaired testosterone production in mature male D. delphis. The application of this integrative analytical approach to cetacean contaminant analysis allows for inference of the biological consequences of accumulation of HOCs and prioritization of compounds for future environmental toxicology research.

Local adaptation to osmotic environment in killifish, Fundulus heteroclitus, is supported by divergence in swimming performance but not by differences in excess post-exercise oxygen consumption or aerobic scope.

Regulation of internal ion homeostasis is essential for fishes inhabiting environments where salinities differ from their internal concentrations. It is hypothesized that selection will reduce energetic costs of osmoregulation in a population's native osmotic habitat, producing patterns of local adaptation. Killifish, Fundulus heteroclitus, occupy estuarine habitats where salinities range from fresh to seawater. Populations inhabiting an environmental salinity gradient differ in physiological traits associated with acclimation to acute salinity stress, consistent with local adaptation. Similarly, metabolic rates differ in populations adapted to different temperatures, but have not been studied in regard to salinity. We investigated evidence for local adaptation between populations of killifish native to fresh and brackish water habitats. Aerobic scope (the difference between minimum and maximum metabolic rates), excess post-exercise oxygen consumption, and swimming performance (time and distance to reach exhaustion) were used as proxies for fitness in fresh and brackish water treatments. Swimming performance results supported local adaptation; fish native to brackish water habitats performed significantly better than freshwater-native fish at high salinity while low salinity performance was similar between populations. However, results from metabolic measures did not support this conclusion; both populations showed an increase in resting metabolic rate and a decrease of aerobic scope in fresh water. Similarly, excess post-exercise oxygen consumption was higher for both populations in fresh than in brackish water. While swimming results suggest that environmentally dependent performance differences may be a result of selection in divergent osmotic environments, the differences between populations are not coupled with divergence in metabolic performance.

Natural selection canalizes expression variation of environmentally induced plasticity-enabling genes.

Many organisms survive fluctuating and extreme environmental conditions by manifesting multiple distinct phenotypes during adulthood by means of developmental processes that enable phenotypic plasticity. We report on the discovery of putative plasticity-enabling genes that are involved in transforming the gill of the euryhaline teleost fish, Fundulus heteroclitus, from its freshwater to its seawater gill-type, a process that alters both morphology and function. Gene expression that normally enables osmotic plasticity is inhibited by arsenic. Gene sets defined by antagonistic interactions between arsenic and salinity show reduced transcriptional variation among individual fish, suggesting unusually accurate and precise regulatory control of these genes, consistent with the hypothesis that they participate in a canalized developmental response. We observe that natural selection acts to preserve canalized gene expression in populations of killifish that are most tolerant to abrupt salinity change and that these populations show the least variability in their transcription of genes enabling plasticity of the gill. We found that genes participating in this highly canalized and conserved plasticity-enabling response had significantly fewer and less complex associations with transcriptional regulators than genes that respond only to arsenic or salinity. Collectively these findings, which are drawn from the relationships between environmental challenge, plasticity, and canalization among populations, suggest that the selective processes that facilitate phenotypic plasticity do so by targeting the regulatory networks that gives rise to the response. These findings also provide a generalized, conceptual framework of how genes might interact with the environment and evolve toward the development of plastic traits.

Reciprocal osmotic challenges reveal mechanisms of divergence in phenotypic plasticity in the killifish Fundulus heteroclitus.

The killifish Fundulus heteroclitus is an estuarine species with broad physiological plasticity, enabling acclimation to diverse stressors. Previous work suggests that freshwater populations expanded their physiology to accommodate low salinity environments; however, it is unknown whether this compromises their tolerance to high salinity. We used a comparative approach to investigate the mechanisms of a derived freshwater phenotype and the fate of an ancestral euryhaline phenotype after invasion of a freshwater environment. We compared physiological and transcriptomic responses to high- and low-salinity stress in fresh and brackish water populations and found an enhanced plasticity to low salinity in the freshwater population coupled with a reduced ability to acclimate to high salinity. Transcriptomic data identified genes with a conserved common response, a conserved salinity-dependent response and responses associated with population divergence. Conserved common acclimation responses revealed stress responses and alterations in cell-cycle regulation as important mechanisms in the general osmotic response. Salinity-specific responses included the regulation of genes involved in ion transport, intracellular calcium, energetic processes and cellular remodeling. Genes diverged between populations were primarily those showing salinity-specific expression and included those regulating polyamine homeostasis and the cell cycle. Additionally, when populations were matched with their native salinity, expression patterns were consistent with the concept of 'transcriptomic resilience', suggesting local adaptation. These findings provide insight into the fate of a plastic phenotype after a shift in environmental salinity and help to reveal mechanisms allowing for euryhalinity.

Genomic and physiological footprint of the Deepwater Horizon oil spill on resident marsh fishes.

The biological consequences of the Deepwater Horizon oil spill are unknown, especially for resident organisms. Here, we report results from a field study tracking the effects of contaminating oil across space and time in resident killifish during the first 4 mo of the spill event. Remote sensing and analytical chemistry identified exposures, which were linked to effects in fish characterized by genome expression and associated gill immunohistochemistry, despite very low concentrations of hydrocarbons remaining in water and tissues. Divergence in genome expression coincides with contaminating oil and is consistent with genome responses that are predictive of exposure to hydrocarbon-like chemicals and indicative of physiological and reproductive impairment. Oil-contaminated waters are also associated with aberrant protein expression in gill tissues of larval and adult fish. These data suggest that heavily weathered crude oil from the spill imparts significant biological impacts in sensitive Louisiana marshes, some of which remain for over 2 mo following initial exposures.

Small groups, contexts, and civic engagement: A multilevel analysis of United States Congregational Life Survey data.

Prior research suggests that church-goers are more civically engaged than their non-church-going counterparts. Little is known, however, about how the popular phenomenon of small groups factors into this equation. In the present study, we examine relationships between small group participation at individual and congregation levels and civic engagement. Using multilevel modeling and national data on congregations and individuals from the U.S. Congregational Life Study (n=82,044), we find that: (1) individual-level small group involvement is associated with four measures of civic engagement; (2) congregation-level small group participation is associated with both lower and higher civic engagement in the case of two outcomes; and (3) in the case of three civic outcomes, congregation-level small group participation moderates individual-level small group involvement such that small group members' civic activity more closely resembles the lower civic engagement of small group nonparticipants. In the case of one civic outcome, at high levels of overall small group participation, small group members' civic engagement drops below that of small group nonparticipants. Explanations for these findings, including a "crowding out" effect, are examined including their complex implications for debates regarding small groups, religious involvement, and civic engagement.

Evolutionary genomics of environmental pollution.

Chemical toxins have been a persistent source of evolutionary challenges throughout the history of life, and deep within the genomic storehouse of evolutionary history lay ancient adaptations to diverse chemical poisons. However, the rate of change of contemporary environments mediated by human-introduced pollutants is rapidly screening this storehouse and severely testing the adaptive potential of many species. In this chapter, we briefly review the deep history of evolutionary adaptation to environmental toxins, and then proceed to describe the attributes of stressors and populations that may facilitate contemporary adaptation to pollutants introduced by humans. We highlight that phenotypes derived to enable persistence in polluted habitats may be multi-dimensional, requiring global genome-scale tools and approaches to uncover their mechanistic basis, and include examples of recent progress in the field. The modern tools of genomics offer promise for discovering how pollutants interact with genomes on physiological timescales, and also for discovering what genomic attributes of populations may enable resistance to pollutants over evolutionary timescales. Through integration of these sophisticated genomics tools and approaches with an understanding of the deep historical forces that shaped current populations, a more mature understanding of the mechanistic basis of contemporary ecological-evolutionary dynamics should emerge.

Genomic mechanisms of evolved physiological plasticity in killifish distributed along an environmental salinity gradient.

Adaptive variation tends to emerge clinally along environmental gradients or discretely among habitats with limited connectivity. However, in Atlantic killifish (Fundulus heteroclitus), a population genetic discontinuity appears in the absence of obvious barriers to gene flow along parallel salinity clines and coincides with a physiologically stressful salinity. We show that populations resident on either side of this discontinuity differ in their abilities to compensate for osmotic shock and illustrate the physiological and functional genomic basis of population variation in hypoosmotic tolerance. A population native to a freshwater habitat, upstream of the genetic discontinuity, exhibits tolerance to extreme hypoosmotic challenge, whereas populations native to brackish or marine habitats downstream of the discontinuity lose osmotic homeostasis more severely and take longer to recover. Comparative transcriptomics reveals a core transcriptional response associated with acute and acclimatory responses to hypoosmotic shock and posits unique mechanisms that enable extreme osmotic tolerance. Of the genes that vary in expression among populations, those that are putatively involved in physiological acclimation are more likely to exhibit nonneutral patterns of divergence between freshwater and brackish populations. It is not the well-known effectors of osmotic acclimation, but rather the lesser-known immediate-early responses, that appear important in contributing to population differences.

Independently evolved pollution resistance in four killifish populations is largely explained by few variants of large effect.

The genetic architecture of phenotypic traits can affect the mode and tempo of trait evolution. Human-altered environments can impose strong natural selection, where successful evolutionary adaptation requires swift and large phenotypic shifts. In these scenarios, theory predicts that adaptation is due to a few adaptive variants of large effect, but empirical studies that have revealed the genetic architecture of rapidly evolved phenotypes are rare, especially for populations inhabiting polluted environments. Fundulus killifish have repeatedly evolved adaptive resistance to extreme pollution in urban estuaries. Prior studies, including genome scans for signatures of natural selection, have revealed some of the genes and pathways important for evolved pollution resistance, and provide context for the genotype-phenotype association studies reported here. We created multiple quantitative trait locus (QTL) mapping families using progenitors from four different resistant populations, and using RAD-seq genetically mapped variation in sensitivity (developmental perturbations) following embryonic exposure to a model toxicant PCB-126. We found that one to two large-effect QTL loci accounted for resistance to PCB-mediated developmental toxicity. QTLs harbored candidate genes that govern the regulation of aryl hydrocarbon receptor (AHR) signaling. One QTL locus was shared across all populations and another was shared across three populations. One QTL locus showed strong signatures of recent natural selection in the corresponding wild population but another QTL locus did not. Some candidate genes for PCB resistance inferred from genome scans in wild populations were identified as QTL, but some key candidate genes were not. We conclude that rapidly evolved resistance to the developmental defects normally caused by PCB-126 is governed by few genes of large effect. However, other aspects of resistance beyond developmental phenotypes may be governed by additional loci, such that comprehensive resistance to PCB-126, and to the mixtures of chemicals that distinguish urban estuaries more broadly, may be more genetically complex.