Advancing human disease research with fish evolutionary mutant models.

Model organism research is essential to understand disease mechanisms. However, laboratory-induced genetic models can lack genetic variation and often fail to mimic the spectrum of disease severity. Evolutionary mutant models (EMMs) are species with evolved phenotypes that mimic human disease. EMMs complement traditional laboratory models by providing unique avenues to study gene-by-environment interactions, modular mutations in noncoding regions, and their evolved compensations. EMMs have improved our understanding of complex diseases, including cancer, diabetes, and aging, and illuminated mechanisms in many organs. Rapid advancements of sequencing and genome-editing technologies have catapulted the utility of EMMs, particularly in fish. Fish are the most diverse group of vertebrates, exhibiting a kaleidoscope of specialized phenotypes, many that would be pathogenic in humans but are adaptive in the species' specialized habitat. Importantly, evolved compensations can suggest avenues for novel disease therapies. This review summarizes current research using fish EMMs to advance our understanding of human disease.

Leafy and weedy seadragon genomes connect genic and repetitive DNA features to the extravagant biology of syngnathid fishes.

Seadragons are a remarkable lineage of teleost fishes in the family Syngnathidae, renowned for having evolved male pregnancy. Comprising three known species, seadragons are widely recognized and admired for their fantastical body forms and coloration, and their specific habitat requirements have made them flagship representatives for marine conservation and natural history interests. Until recently, a gap has been the lack of significant genomic resources for seadragons. We have produced gene-annotated, chromosome-scale genome models for the leafy and weedy seadragon to advance investigations of evolutionary innovation and elaboration of morphological traits in seadragons as well as their pipefish and seahorse relatives. We identified several interesting features specific to seadragon genomes, including divergent noncoding regions near a developmental gene important for integumentary outgrowth, a high genome-wide density of repetitive DNA, and recent expansions of transposable elements and a vesicular trafficking gene family. Surprisingly, comparative analyses leveraging the seadragon genomes and additional syngnathid and outgroup genomes revealed striking, syngnathid-specific losses in the family of fibroblast growth factors (FGFs), which likely involve reorganization of highly conserved gene regulatory networks in ways that have not previously been documented in natural populations. The resources presented here serve as important tools for future evolutionary studies of developmental processes in syngnathids and hold value for conservation of the extravagant seadragons and their relatives.

Reaching Latinx Communities with Algorithmic Optimization for SARS-CoV-2 Testing Locations.

The COVID-19 pandemic has disproportionately affected communities of color, including Latinx communities. Oregon Saludable: Juntos Podemos (OSJP) is a randomized clinical trial aimed at reducing this disparity by both increasing access to testing for SARS-CoV-2, the virus that causes COVID-19, for Oregon Latinx community members and studying the effectiveness of health and behavioral health interventions on turnout and health outcomes. OSJP established SARS-CoV-2 testing events at sites across Oregon. A critical early question was how to locate these sites to best serve Latinx community members. To propose sites in each participating county, we implemented an algorithmic approach solving a facilities location problem. This algorithm was based on minimizing driving time from Latinx population centers to SARS-CoV-2 testing locations. OSJP staff presented these proposed testing locations to community partners as a starting place for identifying final testing sites. Due to differences in geography, population distributions, and potential site accessibility, the study sites exhibited variation in how well the algorithmic optimization objectives could be satisfied. From this variation, we inferred the effects of the drive time optimization metric on the likelihood of Latinx community members utilizing SARS-CoV-2 testing services. After controlling for potential confounders, we found that minimizing the drive time optimization metric was strongly correlated with increased turnout among Latinx community members. This paper presents the algorithm and data sources used for site proposals and discusses challenges and opportunities for community-based health promotion research when translating algorithm proposals into action across a range of health outcomes.

Genome-wide analysis facilitates estimation of the amount of male contribution in meiotic gynogenetic three-spined stickleback (Gasterosteus aculeatus).

Gynogenetic embryos - those inheriting only maternal DNA - can be experimentally created by fertilizing eggs with radiation-treated sperm containing inactivated paternal chromosomes. Diploidy in the zygotes can be maintained through prevention of the second meiosis or restored by preventing the first mitosis after the maternal chromosome complement has been replicated. These gynogenetic organisms are useful in many fields including aquaculture, evolutionary biology and genomics. Although gynogenetic organisms have been created in numerous species, the completeness of uni-parental inheritance has often been assumed rather than thoroughly quantified across the genome. Instead, when tests of uni-parental inheritance occur, they typically rely on well-studied genetically determined phenotypes that represent a very small sub-set of the genome. Only assessing small genomic regions for paternal inheritance leaves the question of whether some paternal contributions to offspring might still have occurred. In this study, the authors quantify the efficacy of creating gynogenetic diploid three-spined stickleback fish (Gasterosteus aculeatus). To this end, the authors mirrored previous assessments of paternal contribution using well-studied genetically determined phenotypes including sex and genetically dominant morphological traits but expanded on previous studies using dense restriction site-associated DNA sequencing (RAD-seq) markers in parents and offspring to assess paternal inheritance genome-wide. In the gynogenetic diploids, the authors found no male genotypes underlying their phenotypes of interest - sex and dominant phenotypic traits. Using genome-wide assessments of paternal contribution, nevertheless, the authors found evidence of a small, yet potentially important, amount of paternally "leaked" genetic material. The application of this genome-wide approach identifies the need for more widespread assessment of paternal contributions to gynogenetic animals and promises benefits for many aspects of aquaculture, evolutionary biology and genomics.

Effects of a Health Education Intervention for COVID-19 Prevention in Latinx Communities: A Cluster-Randomized Controlled Trial.

To promote COVID-19 preventive attitudes and behaviors among Latinx individuals, researchers and community partners implemented a culturally tailored health education intervention across 12 Oregon counties from February 2021 through April 2022. We did not identify any significant intervention effects on preventive attitudes and behaviors but did observe significant decreases in psychological distress. Although Latinx individuals' preventive attitudes and behaviors were not associated with the health education intervention, findings suggest the intervention has value in promoting their well-being (ClinicalTrials.gov Identifier: NCT04793464). (Am J Public Health. 2022;112(S9):S923-S927. https://doi.org/10.2105/AJPH.2022.307129).

Evolutionary divergence of a Hoxa2b hindbrain enhancer in syngnathids mimics results of functional assays.

Hoxa2 genes provide critical patterning signals during development, and their regulation and function have been extensively studied. We report a previously uncharacterized significant sequence divergence of a highly conserved hindbrain hoxa2b enhancer element in the family syngnathidae (pipefishes, seahorses, pipehorses, seadragons). We compared the hox cis-regulatory element variation in the Gulf pipefish and two species of seahorse against eight other species of fish, as well as human and mouse. We annotated the hoxa2b enhancer element binding sites across three species of seahorse, four species of pipefish, and one species of ghost pipefish. Finally, we performed in situ hybridization analysis of hoxa2b expression in Gulf pipefish embryos. We found that all syngnathid fish examined share a modified rhombomere 4 hoxa2b enhancer element, despite the fact that this element has been found to be highly conserved across all vertebrates examined previously. Binding element sequence motifs and spacing between binding elements have been modified for the hoxa2b enhancer in several species of pipefish and seahorse, and that the loss of the Prep/Meis binding site and further space shortening happened after ghost pipefish split from the rest of the syngnathid clade. We showed that expression of this gene in rhombomere 4 is lower relative to the surrounding rhombomeres in developing Gulf pipefish embryos, reflecting previously published functional tests for this enhancer. Our findings highlight the benefits of studying highly derived, diverse taxa for understanding of gene regulatory evolution and support the hypothesis that natural mutations can occur in deeply conserved pathways in ways potentially related to phenotypic diversity.

Genome-wide single nucleotide polymorphism markers reveal population structure and dispersal direction of an expanding nuisance algal bloom species.

Species invasion and range expansion are currently under scrutiny due to increasing anthropogenic impact on the natural environment. This is also true for harmful algal blooms, which have been reported to have increased in frequency. However, this research is challenging due to the ephemeral nature, small size and mostly low concentrations of microalgae in the environment. One such species is the nuisance microalga Gonyostomum semen (Raphidophyceae), which has increased in occurrence in northern Europe in recent decades. The question of whether the species has expanded its habitat range or if it was already present in the lakes but was too rare to be detected remains unanswered. The aim of the present study was to determine the genetic structure and dispersal pathways of G. semen using RAD (restriction-site-associated DNA) tag sequencing. For G. semen, which has a huge genome (32 Gbp), we faced particular challenges, but were nevertheless able to recover over 1000 single nucleotide polymorphisms at high coverage. Our data revealed a distinct population genetic structure, demonstrating a divide of western and eastern populations that probably represent different lineages. Despite significant genetic differentiation among lakes, we found only limited isolation-by-distance. While we had expected a pattern of recent expansion northwards, the data demonstrated gene flow from the northeast/east towards the southwest/west. This genetic signature suggests that the observed gene flow may be due to dispersal by autumn migratory birds, which act as dispersal vectors of resistant resting propagules that form at the end of the G. semen blooms.

Evolution of stickleback in 50 years on earthquake-uplifted islands.

How rapidly can animal populations in the wild evolve when faced with sudden environmental shifts? Uplift during the 1964 Great Alaska Earthquake abruptly created freshwater ponds on multiple islands in Prince William Sound and the Gulf of Alaska. In the short time since the earthquake, the phenotypes of resident freshwater threespine stickleback fish on at least three of these islands have changed dramatically from their oceanic ancestors. To test the hypothesis that these freshwater populations were derived from oceanic ancestors only 50 y ago, we generated over 130,000 single-nucleotide polymorphism genotypes from more than 1,000 individuals using restriction site-associated DNA sequencing (RAD-seq). Population genomic analyses of these data support the hypothesis of recent and repeated, independent colonization of freshwater habitats by oceanic ancestors. We find evidence of recurrent gene flow between oceanic and freshwater ecotypes where they co-occur. Our data implicate natural selection in phenotypic diversification and support the hypothesis that the metapopulation organization of this species helps maintain a large pool of genetic variation that can be redeployed rapidly when oceanic stickleback colonize freshwater environments. We find that the freshwater populations, despite population genetic analyses clearly supporting their young age, have diverged phenotypically from oceanic ancestors to nearly the same extent as populations that were likely founded thousands of years ago. Our results support the intriguing hypothesis that most stickleback evolution in fresh water occurs within the first few decades after invasion of a novel environment.

Naturalistic Experimental Designs as Tools for Understanding the Role of Genes and the Environment in Prevention Research.

Before genetic approaches were applied in experimental studies with human populations, they were used by animal and plant breeders to observe, and experimentally manipulate, the role of genes and environment on specific phenotypic or behavioral outcomes. For obvious ethical reasons, the same level of experimental control is not possible in human populations. Nonetheless, there are natural experimental designs in human populations that can serve as logical extensions of the rigorous quantitative genetic experimental designs used by animal and plant researchers. Applying concepts such as cross-fostering and common garden rearing approaches from the life science discipline, we describe human designs that can serve as naturalistic proxies for the controlled quantitative genetic experiments facilitated in life sciences research. We present the prevention relevance of three such human designs: (1) children adopted at birth by parents to whom they are not genetically related (common garden approach); (2) sibling designs where one sibling is reared from birth with unrelated adoptive parents and the other sibling is reared from birth by the biological mother of the sibling pair (cross-fostering approach); and (3) in vitro fertilization designs, including egg donation, sperm donation, embryo donation, and surrogacy (prenatal cross-fostering approach). Each of these designs allows for differentiation of the effects of the prenatal and/or postnatal rearing environment from effects of genes shared between parent and child in naturalistic ways that can inform prevention efforts. Example findings from each design type are provided and conclusions drawn about the relevance of naturalistic genetic designs to prevention science.

Developmental timing differences underlie armor loss across threespine stickleback populations.

Comparing ontogenetic patterns within a well-described evolutionary context aids in inferring mechanisms of change, including heterochronies or deletion of developmental pathways. Because selection acts on phenotypes throughout ontogeny, any within-taxon developmental variation has implications for evolvability. We compare ontogenetic order and timing of locomotion and defensive traits in three populations of threespine stickleback that have evolutionarily divergent adult forms. This analysis adds to the growing understanding of developmental genetic mechanisms of adaptive change in this evolutionary model species by delineating when chondrogenesis and osteogenesis in two derived populations begin to deviate from the developmental pattern in their immediate ancestors. We found that differences in adult defensive morphologies arise through abolished or delayed initiation of these traits rather than via an overall heterochronic shift, that intra-population ontogenetic variation is increased for some derived traits, and that altered armor developmental timing differentiates the derived populations from each other despite parallels in adult lateral plate armor phenotypes. We found that changes in ossified elements of the pelvic armor are linked to delayed and incomplete development of an early-forming pelvic cartilage, and that this disruption likely presages the variable pelvic vestiges documented in many derived populations.

Exogenous iodide ameliorates perchlorate-induced thyroid phenotypes in threespine stickleback.

Perchlorate is a ubiquitous environmental contaminant that has widespread endocrine disrupting effects in vertebrates, including threespine stickleback (Gasterosteus aculeatus). The target of perchlorate is thyroid tissue where it induces changes in the organization, activation, and morphology of thyroid follicles and surrounding tissues. To test the hypothesis that some phenotypes of perchlorate toxicity are not mediated by thyroid hormone, we chronically exposed stickleback beginning at fertilization to perchlorate (10, 30, 100ppm) or control water with and without supplementation of either iodide or thyroxine (T4). Stickleback were sampled across a one-year timespan to identify potential differences in responses to treatment combinations before and after sexual maturation. We found that most thyroid histomorphological phenotypes induced by perchlorate (follicle proliferation, reduced follicle area (adults only), colloid depletion, thyrocyte hypertrophy (subadults only)) were significantly ameliorated by exogenous iodide supplementation. In contrast, treatment with exogenous T4 did not correct any of the thyroid-specific histopathologies induced by perchlorate. Whole-body thyroid hormone concentrations were not significantly affected by perchlorate exposure; however, supplementation with iodide and T4 significantly increased T4 concentrations. This study also revealed an increased erythrocyte area in the thyroid region of perchlorate-exposed adults, while lipid droplet number increased in perchlorate-exposed subadults. Increased erythrocyte area was ameliorated by both iodide and T4, while neither supplement was able to correct lipid droplet number. Our finding on lipid droplets indicates that exposure to perchlorate in early development may have obesogenic effects.

Effects of crossovers between homeologs on inheritance and population genomics in polyploid-derived salmonid fishes.

A whole genome duplication occurred in the ancestor of all salmonid fishes some 50-100 million years ago. Early inheritance studies with allozymes indicated that loci in the salmonid genome are inherited disomically in females. However, some pairs of duplicated loci showed patterns of inheritance in males indicating pairing and recombination between homeologous chromosomes. Nearly 20% of loci in the salmonid genome are duplicated and share the same alleles (isoloci), apparently due to homeologous recombination. Half-tetrad analysis revealed that isoloci tend to be telomeric. These results suggested that residual tetrasomic inheritance of isoloci results from homeologous recombination near chromosome ends and that continued disomic inheritance resulted from homologous pairing of centromeric regions. Many current genetic maps of salmonids are based on single nucleotide polymorphisms and microsatellites that are no longer duplicated. Therefore, long sections of chromosomes on these maps are poorly represented, especially telomeric regions. In addition, preferential multivalent pairing of homeologs from the same species in F1 hybrids results in an excess of nonparental gametes (so-called pseudolinkage). We consider how not including duplicated loci has affected our understanding of population and evolutionary genetics of salmonids, and we discuss how incorporating these loci will benefit our understanding of population genomics.

Perchlorate disrupts embryonic androgen synthesis and reproductive development in threespine stickleback without changing whole-body levels of thyroid hormone.

Perchlorate, an environmental contaminant, disrupts normal functioning of the thyroid. We previously showed that perchlorate disrupts behavior and gonad development, and induces external morphological changes in a vertebrate model organism, the threespine stickleback. Whether perchlorate alters these phenotypes via a thyroid-mediated mechanism, and the extent to which the effects depend on dose, are unknown. To address these questions, we chronically exposed stickleback to control conditions and to three concentrations of perchlorate (10, 30 and 100ppm) at various developmental stages from fertilization to reproductive maturity. Adults chronically exposed to perchlorate had increased numbers of thyroid follicles and decreased numbers of thyrocytes. Surprisingly, T4 and T3 levels in larval, juvenile, and adult whole fish chronically exposed to perchlorate did not differ from controls, except at the lowest perchlorate dose, suggesting a non-monotonic dose response curve. We found no detectable abnormalities in external phenotype at any dose of perchlorate, indicating that the increased number of thyroid follicles compensated for the disruptive effects of these doses. In contrast to external morphology, gonadal development was altered substantially, with the highest dose of perchlorate causing the largest effects. Perchlorate increased the number both of early stage ovarian follicles in females and of advanced spermatogenic stages in males. Perchlorate also disrupted embryonic androgen levels. We conclude that chronic perchlorate exposure may not result in lasting adult gross morphological changes but can produce lasting modifications to gonads when compensation of T3 and T4 levels occurs by thyroid follicle hyperplasia. Perchlorate may therefore affect vertebrate development via both thyroidal and non-thyroidal mechanisms.

Developmental timing of sodium perchlorate exposure alters angiogenesis, thyroid follicle proliferation and sexual maturation in stickleback.

Perchlorate, a common aquatic contaminant, is well known to disrupt homeostasis of the hypothalamus-pituitary-thyroid axis. This study utilizes the threespine stickleback (Gasterosteus aculeatus) fish to determine if perchlorate exposure during certain windows of development has morphological effects on thyroid and gonads. Fish were moved from untreated water to perchlorate-contaminated water (30 and 100mg/L) starting at 0, 3, 7, 14, 21, 42, 154 and 305 days post fertilization until approximately one year old. A reciprocal treatment (fish in contaminated water switched to untreated water) was conducted on the same schedule. Perchlorate exposure increased angiogenesis and follicle proliferation in thyroid tissue, delayed gonadal maturity, and skewed sex ratios toward males; effects depended on concentration and timing of exposure. This study demonstrates that perchlorate exposure beginning during the first 42 days of development has profound effects on stickleback reproductive and thyroid tissues, and by implication can impact population dynamics. Long-term exposure studies that assess contaminant effects at various stages of development provide novel information to characterize risk to aquatic organisms, to facilitate management of resources, and to determine sensitive developmental windows for further study of underlying mechanisms.

Perchlorate exposure does not modulate temporal variation of whole-body thyroid and androgen hormone content in threespine stickleback.

Previously we showed that exposure of threespine stickleback (Gasterosteus aculeatus) to the endocrine disruptor perchlorate results in pronounced structural changes in thyroid and gonad, while surprisingly, whole-body thyroid hormone concentrations remain unaffected. To test for hormone titer variations on a finer scale, we evaluated the interactive effects of time (diel and reproductive season) and perchlorate exposure on whole-body contents of triiodothyronine (T3), thyroxine (T4), and 11-ketotestosterone (11-KT) in captive stickleback. Adult stickleback were exposed to 100ppm perchlorate or control water and sampled at 4-h intervals across the 24-hday and at one time-point (1100h) weekly across the reproductive season (May-July). Neither whole-body T3 nor T4 concentration significantly differed across the day in control or perchlorate treated stickleback. Across the reproductive season, whole-body T3 concentration remained stable while T4 significantly increased. However, neither hormone concentration was significantly affected by perchlorate, verifying our previous studies. The concentration of whole-body 11-KT, a major fish androgen, displayed significant diel variation and also steadily declined across the reproductive season in untreated males; perchlorate exposure did not influence the concentration of 11-KT in either diel or reproductive season schedules. Diel and reproductive season variations in 11-KT content in male stickleback are likely related to reproductive physiology and behavior. The observed increase in T4 content across the reproductive season may be reflective of increased energy investment in reproduction near the end of the life cycle.

Developmental timing of perchlorate exposure alters threespine stickleback dermal bone.

Adequate levels of thyroid hormone are critical during development and metamorphosis, and for maintaining metabolic homeostasis. Perchlorate, a common contaminant of water sources, inhibits thyroid function in vertebrates. We utilized threespine stickleback (Gasterosteus aculeatus) to determine if timing of perchlorate exposure during development impacts adult dermal skeletal phenotypes. Fish were exposed to water contaminated with perchlorate (30mg/L or 100mg/L) beginning at 0, 3, 7, 14, 21, 42, 154 or 305days post fertilization until sexual maturity at 1year of age. A reciprocal treatment moved stickleback from contaminated to clean water on the same schedule providing for different stages of initial exposure and different treatment durations. Perchlorate exposure caused concentration-dependent significant differences in growth for some bony traits. Continuous exposure initiated within the first 21days post fertilization had the greatest effects on skeletal traits. Exposure to perchlorate at this early stage can result in small traits or abnormal skeletal morphology of adult fish which could affect predator avoidance and survival.

The transgenerational effects of heat stress in the nematode Caenorhabditis remanei are negative and rapidly eliminated under direct selection for increased stress resistance in larvae.

Parents encountering stress environments can influence the phenotype of their offspring in a form of transgenerational phenotypic plasticity that has the potential to be adaptive if offspring are thereby better able to deal with future stressors. Here, we test for the existence of anticipatory parental effects in the heat stress response in the highly polymorphic nematode Caenorhabditis remanei. Rather providing an anticipatory response, parents subject to a prior heat stress actually produce offspring that are less able to survive a severe heat shock. Selection on heat shock resistance within the larvae via experimental evolution leads to a loss of sensitivity (robustness) to environmental variation during both the parental and larval periods. Whole genome transcriptional analysis of both ancestor and selected lines shows that there is weak correspondence between genetic pathways induced via temperature shifts during parental and larval periods. Parental effects can evolve very rapidly via selection acting directly on offspring.

Single Cell RNA Sequencing Provides Clues for the Developmental Genetic Basis of Syngnathidae's Evolutionary Adaptations.

Seahorses, pipefishes, and seadragons are fishes from the family Syngnathidae that have evolved extraordinary traits including male pregnancy, elongated snouts, loss of teeth, and dermal bony armor. The developmental genetic and cellular changes that led to the evolution of these traits are largely unknown. Recent syngnathid genomes revealed suggestive gene content differences and provide the opportunity for detailed genetic analyses. We created a single cell RNA sequencing atlas of Gulf pipefish embryos to understand the developmental basis of four traits: derived head shape, toothlessness, dermal armor, and male pregnancy. We completed marker gene analyses, built genetic networks, and examined spatial expression of select genes. We identified osteochondrogenic mesenchymal cells in the elongating face that express regulatory genes bmp4, sfrp1a, and prdm16. We found no evidence for tooth primordia cells, and we observed re-deployment of osteoblast genetic networks in developing dermal armor. Finally, we found that epidermal cells expressed nutrient processing and environmental sensing genes, potentially relevant for the brooding environment. The examined pipefish evolutionary innovations are composed of recognizable cell types, suggesting derived features originate from changes within existing gene networks. Future work addressing syngnathid gene networks across multiple stages and species is essential for understanding how their novelties evolved.

Stacks: an analysis tool set for population genomics.

Massively parallel short-read sequencing technologies, coupled with powerful software platforms, are enabling investigators to analyse tens of thousands of genetic markers. This wealth of data is rapidly expanding and allowing biological questions to be addressed with unprecedented scope and precision. The sizes of the data sets are now posing significant data processing and analysis challenges. Here we describe an extension of the Stacks software package to efficiently use genotype-by-sequencing data for studies of populations of organisms. Stacks now produces core population genomic summary statistics and SNP-by-SNP statistical tests. These statistics can be analysed across a reference genome using a smoothed sliding window. Stacks also now provides several output formats for several commonly used downstream analysis packages. The expanded population genomics functions in Stacks will make it a useful tool to harness the newest generation of massively parallel genotyping data for ecological and evolutionary genetics.

The population structure and recent colonization history of Oregon threespine stickleback determined using restriction-site associated DNA-sequencing.

Understanding how genetic variation is partitioned across genomes within and among populations is a fundamental problem in ecological and evolutionary genetics. To address this problem, we studied the threespine stickleback fish, which has repeatedly undergone parallel phenotypic and genetic differentiation when oceanic fish have invaded freshwater habitats. While significant evolutionary genetic research has been performed using stickleback from geographic regions that have been deglaciated in the last 20 000 years, less research has focused on freshwater populations that predate the last glacial maximum. We performed restriction-site associated DNA-sequencing (RAD-seq) based population genomic analyses on stickleback from across Oregon, which was not glaciated during the last maximum. We sampled stickleback from coastal, Willamette Basin and central Oregon sites, analysed their genetic diversity using RAD-seq, performed structure analyses, reconstructed their phylogeographic history and tested the hypothesis of recent stickleback introduction into central Oregon, where incidence of this species was only recently documented. Our results showed a clear phylogeographic break between coastal and inland populations, with oceanic populations exhibiting the lowest levels of divergence from one another. Willamette Basin and central Oregon populations formed a clade of closely related populations, a finding consistent with a recent introduction of stickleback into central Oregon. Finally, genome-wide analysis of genetic diversity (π) and correlations of alleles within individuals in subpopulations (FIS) supported a role for introgressive hybridization in coastal populations and a recent expansion in central Oregon. Our results exhibit the power of next-generation sequencing genomic approaches such as RAD-seq to identify both historical population structure and recent colonization history.