Evolution and development of facial bone morphology in threespine sticklebacks

Computational anatomy and geometric shape analysis enables analysis of complex craniofacial phenotypes in zebrafish

Due to the complexity of fish skulls, previous attempts to classify craniofacial phenotypes have relied on qualitative features or sparce 2D landmarks. In this work we aim to identify previously unknown 3D craniofacial phenotypes with a semiautomated pipeline in adult zebrafish mutants. We first estimate a synthetic 'normative' zebrafish template using MicroCT scans from a sample pool of wild-type animals using the Advanced Normalization Tools (ANTs). We apply a computational anatomy (CA) approach to quantify the phenotype of zebrafish with disruptions in bmp1a, a gene implicated in later skeletal development and whose human ortholog when disrupted is associated with Osteogenesis Imperfecta. Compared to controls, the bmp1a fish have larger otoliths, larger normalized centroid sizes, and exhibit shape differences concentrated around the operculum, anterior frontal, and posterior parietal bones. Moreover, bmp1a fish differ in the degree of asymmetry. Our CA approach offers a potential pipeline for high-throughput screening of complex fish craniofacial shape to discover novel phenotypes for which traditional landmarks are too sparce to detect. The current pipeline successfully identifies areas of variation in zebrafish mutants, which are an important model system for testing genome to phenome relationships in the study of development, evolution, and human diseases. This article has an associated First Person interview with the first author of the paper.

Multivariate analysis of morphology, behaviour, growth and developmental timing in hybrids brings new insights into the divergence of sympatric Arctic charr morphs

BACKGROUND

Studying the development of fitness related traits in hybrids from populations diverging in sympatry is a fundamental approach to understand the processes of speciation. However, such traits are often affected by covariance structures that complicate the comprehension of these processes, especially because the interactive relationships between traits of different nature (e.g. morphology, behaviour, life-history) remain largely unknown in this context. In a common garden setup, we conducted an extensive examination of a large suit of traits putatively involved in the divergence of two morphs of Arctic charr (Salvelinus alpinus), and investigated the consequences of potential patterns of trait covariance on the phenotype of their hybrids. These traits were measured along ontogeny and involved growth, yolk sac resorption, developmental timing (hatching and the onset of exogeneous feeding), head morphology and feeding behaviour.

RESULTS

Growth trajectories provided the strongest signal of phenotypic divergence between the two charr. Strikingly, the first-generation hybrids did not show intermediate nor delayed growth but were similar to the smallest morph, suggesting parental biases in the inheritance of growth patterns. However, we did not observe extensive multivariate trait differences between the two morphs and their hybrids. Growth was linked to head morphology (suggesting that morphological variations in early juveniles relate to simple allometric effects) but this was the only strong signal of covariance observed between all the measured traits. Furthermore, we did not report evidence for differences in overall phenotypic variance between morphs, nor for enhanced phenotypic variability in their hybrids.

CONCLUSION

Our study shed light on the multivariate aspect of development in a context of adaptive divergence. The lack of evidence for the integration of most traits into a single covariance structure suggested that phenotypic constraints may not always favour nor impede divergence toward ecological niches differing in numerous physical and ecological variables, as observed in the respective habitats of the two charr. Likewise, the role of hybridization as a disruptive agent of trait covariance may not necessarily be significant in the evolution of populations undergoing resource polymorphism.

Allometric relationships in morphological traits associated with foraging, swimming ability, and predator defense reveal adaptations toward brackish and freshwater environments in the threespine stickleback

Freshwater colonization by threespine stickleback has led to divergence in morphology between ancestral marine and derived freshwater populations, making them ideal for studying natural selection on phenotypes. In an open brackish-freshwater system, we previously discovered two genetically distinct stickleback populations that also differ in geometric shape: one mainly found in the brackish water lagoon and one throughout the freshwater system. As shape and size are not perfectly correlated, the aim of this study was to identify the morphological trait(s) that separated the populations in geometric shape. We measured 23 phenotypes likely to be important for foraging, swimming capacity, and defense against predation. The lateral plate morphs in freshwater displayed few significant changes in trait sizes, but the low plated expressed feeding traits more associated with benthic habitats. When comparing the completely plated genetically assigned populations, the freshwater, the hybrids, the migrants and the lagoon fish, many of the linear traits had different slopes and intercepts in trait-size regressions, precluding our ability to directly compare all traits simultaneously, which most likely results from low variation in body length for the lagoon and migrant population. We found the lagoon stickleback population to be more specialized toward the littoral zone, displaying benthic traits such as large, deep bodies with smaller eyes compared to the freshwater completely plated morph. Further, the lagoon and migrant fish had an overall higher body coverage of lateral plates compared to freshwater fish, and the dorsal and pelvic spines were longer. Evolutionary constraints due to allometric scaling relationships could explain the observed, overall restricted, differences in morphology between the sticklebacks in this study, as most traits have diversified in common allometric trajectories. The observed differences in foraging and antipredation traits between the fish with a lagoon and freshwater genetic signature are likely a result of genetic or plastic adaptations toward brackish and freshwater environments.

Early life disruption of the microbiota affects organ development and cytokine gene expression in threespine stickleback

The microbiota that inhabits vertebrates exerts strong effects on host physiology and can be crucial to the development of a normal phenotype. This includes development of the immune system, somatic growth and maintenance, and morphogenesis. However, the genetic background of the host can also affect these life history traits. To this end, we investigated the effects of the microbiota on growth, development, and immune gene expression on two populations of threespine stickleback (Gasterosteus aculeatus), one anadromous and one freshwater. We tested the hypotheses that microbial colonization and the genetic background of the host would affect survival, cytokine gene expression, growth, and development. We raised in vitro crosses of stickleback larvae with and without conventional microbiota. We then exposed all these treatments to Vibrio anguillarum, a potential fish pathogen, in a full factorial design. We found stickleback raised without conventional microbiota had smaller swim bladders relative to those raised with conventional microbiota. Stickleback raised with conventional microbiota exhibited small increases in cytokine gene expression. We found no differences in growth or survival regardless of treatment. These results are consistent with other investigations that show microbiota disruption, in early life, can alter host organ and tissue development and immune responses.

Genetic divergence outpaces phenotypic evolution among threespine stickleback populations in old freshwater habitats

Species such as threespine stickleback (Gasterosteus aculeatus) that inhabit divergent selective environments and that have diversified on different time scales can be of value for understanding evolutionary processes. Here we synthesize high-resolution genotypic and phenotypic data to explore a largely unstudied distribution of threespine stickleback populations living in oceanic and freshwater habitats along coastal and inland regions of Oregon. Many inland aquatic habitats of Oregon remained unglaciated during the last ice age, meaning that some extant Oregon lake and river stickleback may have descended from freshwater populations established long before more well-studied, post-glacial freshwater populations. To address the degree of congruence between genetic and phenotypic divergence, we directly compared Oregon stickleback to much younger (post-glacial) Alaskan populations. We found phenotypic variation in Oregon stickleback to be primarily partitioned between oceanic and freshwater habitats, as has been documented in other stickleback systems. However, the main axis of genetic divergence was between coastal and inland regions regardless of habitat type. Furthermore, when comparing patterns between Oregon and Alaska we found similar levels of phenotypic divergence, but much greater genetic divergence among Oregon's populations. The Oregon stickleback system therefore appears well suited for future studies linking genotypic and phenotypic change, further extending the utility of this small fish to provide general insights into evolutionary processes.

Developmental tuning of mineralization drives morphological diversity of gill cover bones in sculpins and their relatives

The role of osteoblast placement in skeletal morphological variation is relatively well understood, but alternative developmental mechanisms affecting bone shape remain largely unknown. Specifically, very little attention has been paid to variation in later mineralization stages of intramembranous ossification as a driver of morphological diversity. We discover the occurrence of specific, sometimes large, regions of nonmineralized osteoid within bones that also contain mineralized tissue. We show through a variety of histological, molecular, and tomographic tests that this “extended” osteoid material is most likely nonmineralized bone matrix. This tissue type is a significant determinant of gill cover bone shape in the teleostean suborder Cottoidei. We demonstrate repeated evolution of extended osteoid in Cottoidei through ancestral state reconstruction and test for an association between extended osteoid variation and habitat differences among species. Through measurement of extended osteoid at various stages of gill cover development in species across the phylogeny, we gain insight into possible evolutionary developmental origins of the trait. We conclude that this fine‐tuned developmental regulation of bone matrix mineralization reflects heterochrony at multiple biological levels and is a novel mechanism for the evolution of diversity in skeletal morphology. This research lays the groundwork for a new model in which to study bone mineralization and evolutionary developmental processes, particularly as they may relate to adaptation during a prominent evolutionary radiation of fishes.

Craniofacial skeleton of MEXICAN tetra (Astyanax mexicanus): As a bone disease model

A small fresh water fish, the Mexican tetra (Astyanax mexicanus) is a novel animal model in evolutionary developmental biology. The existence of morphologically distinct surface and cave morphs of this species allows simultaneous comparative analysis of phenotypic changes at different life stages. The cavefish harbors many favorable constructive traits (i.e., large jaws with an increased number of teeth, neuromast cells, enlarged olfactory pits and excess storage of adipose tissues) and regressive traits (i.e., reduced eye structures and pigmentation) which are essential for cave adaptation. A wide spectrum of natural craniofacial morphologies can be observed among the different cave populations. Recently, the Mexican tetra has been identified as a human disease model. The fully sequenced genome along with modern genome editing tools has allowed researchers to generate transgenic and targeted gene knockouts with phenotypes that resemble human pathological conditions. This review will discuss the anatomy of the craniofacial skeleton of A. mexicanus with a focus on morphologically variable facial bones, jaws that house continuously replacing teeth and pharyngeal skeleton. Furthermore, the possible applications of this model animal in identifying human congenital and metabolic skeletal disorders is addressed. Developmental Dynamics 248:153-161, 2019. © 2018 Wiley Periodicals, Inc.

Neural crest and the origin of species-specific pattern

For well over half of the 150 years since the discovery of the neural crest, the special ability of these cells to function as a source of species-specific pattern has been clearly recognized. Initially, this observation arose in association with chimeric transplant experiments among differentially pigmented amphibians, where the neural crest origin for melanocytes had been duly noted. Shortly thereafter, the role of cranial neural crest cells in transmitting species-specific information on size and shape to the pharyngeal arch skeleton as well as in regulating the timing of its differentiation became readily apparent. Since then, what has emerged is a deeper understanding of how the neural crest accomplishes such a presumably difficult mission, and this includes a more complete picture of the molecular and cellular programs whereby neural crest shapes the face of each species. This review covers studies on a broad range of vertebrates and describes neural-crest-mediated mechanisms that endow the craniofacial complex with species-specific pattern. A major focus is on experiments in quail and duck embryos that reveal a hierarchy of cell-autonomous and non-autonomous signaling interactions through which neural crest generates species-specific pattern in the craniofacial integument, skeleton, and musculature. By controlling size and shape throughout the development of these systems, the neural crest underlies the structural and functional integration of the craniofacial complex during evolution.

Study of morphological variation of northern Neotropical Ariidae reveals conservatism despite macrohabitat transitions

BACKGROUND

Morphological convergence triggered by trophic adaptations is a common pattern in adaptive radiations. The study of shape variation in an evolutionary context is usually restricted to well-studied fish models. We take advantage of the recently revised systematics of New World Ariidae and investigate skull shape evolution in six genera of northern Neotropical Ariidae. They constitute a lineage that diversified in the marine habitat but repeatedly adapted to freshwater habitats. 3D geometric morphometrics was applied for the first time in catfish skulls and phylogenetically informed statistical analyses were performed to test for the impact of habitat on skull diversification after habitat transition in this lineage.

RESULTS

We found that skull shape is conserved throughout phylogeny. A morphospace analysis revealed that freshwater and marine species occupy extreme ends of the first principal component axis and that they exhibit similar Procrustes variances. Yet freshwater species occupy the smallest shape space compared to marine and brackish species (based on partial disparity), and marine and freshwater species have the largest Procrustes distance to each other. We observed a single case of shape convergence as derived from 'C-metrics', which cannot be explained by the occupation of the same habitat.

CONCLUSIONS

Although Ariidae occupy such a broad spectrum of different habitats from sea to freshwater, the morphospace analysis and analyses of shape and co-variation with habitat in a phylogenetic context shows that conservatism dominates skull shape evolution among ariid genera.

Ancient genomic variation underlies repeated ecological adaptation in young stickleback populations

Adaptation in the wild often involves standing genetic variation (SGV), which allows rapid responses to selection on ecological timescales. However, we still know little about how the evolutionary histories and genomic distributions of SGV influence local adaptation in natural populations. Here, we address this knowledge gap using the threespine stickleback fish (Gasterosteus aculeatus) as a model. We extend restriction site-associated DNA sequencing (RAD-seq) to produce phased haplotypes approaching 700 base pairs (bp) in length at each of over 50,000 loci across the stickleback genome. Parallel adaptation in two geographically isolated freshwater pond populations consistently involved fixation of haplotypes that are identical-by-descent. In these same genomic regions, sequence divergence between marine and freshwater stickleback, as measured by d_XY , reaches tenfold higher than background levels and genomic variation is structured into distinct marine and freshwater haplogroups. By combining this dataset with a de novo genome assembly of a related species, the ninespine stickleback (Pungitius pungitius), we find that this habitat-associated divergent variation averages six million years old, nearly twice the genome-wide average. The genomic variation that is involved in recent and rapid local adaptation in stickleback has therefore been evolving throughout the 15-million-year history since the two species lineages split. This long history of genomic divergence has maintained large genomic regions of ancient ancestry that include multiple chromosomal inversions and extensive linked variation. These discoveries of ancient genetic variation spread broadly across the genome in stickleback demonstrate how selection on ecological timescales is a result of genome evolution over geological timescales, and vice versa.

Detecting differential copy number variation between groups of samples

We present a method to detect copy number variants (CNVs) that are differentially present between two groups of sequenced samples. We use a finite-state transducer where the emitted read depth is conditioned on the mappability and GC-content of all reads that occur at a given base position. In this model, the read depth within a region is a mixture of binomials, which in simulations matches the read depth more closely than the often-used negative binomial distribution. The method analyzes all samples simultaneously, preserving uncertainty as to the breakpoints and magnitude of CNVs present in an individual when it identifies CNVs differentially present between the two groups. We apply this method to identify CNVs that are recurrently associated with postglacial adaptation of marine threespine stickleback (Gasterosteus aculeatus) to freshwater. We identify 6664 regions of the stickleback genome, totaling 1.7 Mbp, which show consistent copy number differences between marine and freshwater populations. These deletions and duplications affect both protein-coding genes and cis-regulatory elements, including a noncoding intronic telencephalon enhancer of DCHS1. The functions of the genes near or included within the 6664 CNVs are enriched for immunity and muscle development, as well as head and limb morphology. Although freshwater stickleback have repeatedly evolved from marine populations, we show that freshwater stickleback also act as reservoirs for ancient ancestral sequences that are highly conserved among distantly related teleosts, but largely missing from marine stickleback due to recent selective sweeps in marine populations.

The genetic and molecular architecture of phenotypic diversity in sticklebacks

A major goal of evolutionary biology is to identify the genotypes and phenotypes that underlie adaptation to divergent environments. Stickleback fish, including the threespine stickleback (Gasterosteus aculeatus) and the ninespine stickleback (Pungitius pungitius), have been at the forefront of research to uncover the genetic and molecular architecture that underlies phenotypic diversity and adaptation. A wealth of quantitative trait locus (QTL) mapping studies in sticklebacks have provided insight into long-standing questions about the distribution of effect sizes during adaptation as well as the role of genetic linkage in facilitating adaptation. These QTL mapping studies have also provided a basis for the identification of the genes that underlie phenotypic diversity. These data have revealed that mutations in regulatory elements play an important role in the evolution of phenotypic diversity in sticklebacks. Genetic and molecular studies in sticklebacks have also led to new insights on the genetic basis of repeated evolution and suggest that the same loci are involved about half of the time when the same phenotypes evolve independently. When the same locus is involved, selection on standing variation and repeated mutation of the same genes have both contributed to the evolution of similar phenotypes in independent populations.This article is part of the themed issue 'Evo-devo in the genomics era, and the origins of morphological diversity'.

Divergence of insulin superfamily ligands, receptors and Igf binding proteins in marine versus freshwater stickleback: Evidence of selection in known and novel genes

Three-spine stickleback (Gasterosteus aculeatus) is a teleost model for understanding genetic, physiological and morphological changes accompanying freshwater (FW) adaptation. There is growing evidence that the insulin superfamily plays important roles in traits involved in marine and FW adaptation. We performed a candidate gene analysis to look for evidence of selection on 33 insulin superfamily ligand-receptor genes and insulin-like growth factor binding proteins (Igfbp's) in stickleback. Using genotype data from 11 marine and 10 FW populations, we calculated the number of SNPs per site in regulatory and intronic regions, the number of synonymous and nonsynonymous mutations in coding regions, Wright's fixation index (Fst), and performed t-tests to identify SNPs with divergent genotype frequencies between marine/FW versus Atlantic/Pacific populations. Next, we analysed genome-wide transcriptome data from eight tissues to assess differential gene expression. Two Igfbp's (Igfbp2a and Igfbp5a) show evidence of divergent adaptation between life-history types, and a cluster of nonsynonymous mutations in Igfbp5a exhibit high Fst in exons apparently alternatively spliced in gill. We find evidence of selection on the relaxin family ligand-receptor gene pair, Insl3-Rxfp2, known to be involved in male spermatogenesis and bone metabolism, and in the 5' regulatory region of Igf2. We also confirmed the gene and coding sequence of two unannotated relaxin family ligands. These analyses underscore the utility of candidate gene studies and indicate directions for further exploration of the function of insulin superfamily genes in FW adaptation.

Host Genotype and Microbiota Contribute Asymmetrically to Transcriptional Variation in the Threespine Stickleback Gut

Recent studies of interactions between hosts and their resident microbes have revealed important ecological and evolutionary consequences that emerge from these complex interspecies relationships, including diseases that occur when the interactions go awry. Given the preponderance of these interactions, we hypothesized that effects of the microbiota on gene expression in the developing gut—an important aspect of host biology—would be pervasive, and that these effects would be both comparable in magnitude to and contingent on effects of the host genetic background. To evaluate the effects of the microbiota, host genotype, and their interaction on gene expression in the gut of a genetically diverse, gnotobiotic host model, the threespine stickleback (Gasterosteus aculeatus), we compared RNA-seq data among 84 larval fish. Surprisingly, we found that stickleback population and family differences explained substantially more gene expression variation than the presence of microbes. Expression levels of 72 genes, however, were affected by our microbiota treatment. These genes, including many associated with innate immunity, comprise a tractable subset of host genetic factors for precise, systems-level study of host–microbe interactions in the future. Importantly, our data also suggest subtle signatures of a statistical interaction between host genotype and the microbiota on expression patterns of genetic pathways associated with innate immunity, coagulation and complement cascades, focal adhesion, cancer, and peroxisomes. These genotype-by-environment interactions may prove to be important leads to the understanding of host genetic mechanisms commonly at the root of sometimes complex molecular relationships between hosts and their resident microbes.

Teleosts as Model Organisms To Understand Host-Microbe Interactions

Host-microbe interactions are influenced by complex host genetics and environment. Studies across animal taxa have aided our understanding of how intestinal microbiota influence vertebrate development, disease, and physiology. However, traditional mammalian studies can be limited by the use of isogenic strains, husbandry constraints that result in small sample sizes and limited statistical power, reliance on indirect characterization of gut microbial communities from fecal samples, and concerns of whether observations in artificial conditions are actually reflective of what occurs in the wild. Fish models are able to overcome many of these limitations. The extensive variation in the physiology, ecology, and natural history of fish enriches studies of the evolution and ecology of host-microbe interactions. They share physiological and immunological features common among vertebrates, including humans, and harbor complex gut microbiota, which allows identification of the mechanisms driving microbial community assembly. Their accelerated life cycles and large clutch sizes and the ease of sampling both internal and external microbial communities make them particularly well suited for robust statistical studies of microbial diversity. Gnotobiotic techniques, genetic manipulation of the microbiota and host, and transparent juveniles enable novel insights into mechanisms underlying development of the digestive tract and disease states. Many diseases involve a complex combination of genes which are difficult to manipulate in homogeneous model organisms. By taking advantage of the natural genetic variation found in wild fish populations, as well as of the availability of powerful genetic tools, future studies should be able to identify conserved genes and pathways that contribute to human genetic diseases characterized by dysbiosis.

Evolution of opercle bone shape along a macrohabitat gradient: species identification using mtDNA and geometric morphometric analyses in neotropical sea catfishes (Ariidae)

Transitions between the marine and freshwater macrohabitat have occurred repeatedly in the evolution of teleost fishes. For example, ariid catfishes have moved from freshwater to marine environments, and vice versa. Opercles, a skeletal feature that has been shown to change during such transitions, were subjected to 2D geometric morphometric analyses in order to investigate evolutionary shape changes during habitat transition in ariid catfishes and to test the influence of habitat on shape changes. A mtDNA marker, which proved useful in previous studies, was used to verify species identities. It greatly improved the assignment of specimens to a species, which are difficult to assign by morphology alone. The application of a mtDNA marker confirmed the occurrence of Notarius biffi in Central America, South of El Salvador. Molecular identification together with principal component analysis (PCA) and further morphological inspection of neurocrania indicated the existence of a cryptic species within Bagre pinnimaculatus. Principal component (PC) scores of individual specimens clustered in morphospace by genus rather than by habitat. Strong phylogenetic structure was detected using a permutation test of PC scores of species means on a phylogenetic tree. Calculation of Pagel's λ suggested that opercle shape evolved according to a Brownian model of evolution. Yet canonical variate analysis (CVA) conducted on the habitat groups showed significant differences in opercle shapes among freshwater and marine species. Overall, opercle shape in tropical American Ariidae appears to be phylogenetically constrained. This verifies the application of opercle shape as a taxonomic tool for species identification in fossil ariid catfishes. At the same time, adaptation to freshwater habitats shows characteristic opercle shape trajectories in ariid catfishes, which might be used to detect habitat preferences in fossils.

Craniofacial shape transition across the house mouse hybrid zone: implications for the genetic architecture and evolution of between-species differences

Craniofacial shape differences between taxa have often been linked to environmental adaptation, e.g., new food sources, or have been studied in the context of domestication. Evidence for the genetic basis of such phenotypic differences to date suggests that between-species as well as between-population variation has an oligogenic basis, i.e., few loci of large effect explain most of the variation. In mice, it has been shown that within-population craniofacial variation has a highly polygenic basis, but there are no data regarding the genetic basis of between-species differences in natural populations. Here, we address this question using a phenotype-focused approach. Using 3D geometric morphometrics, we phenotyped a panel of mice derived from a natural hybrid zone between Mus musculus domesticus and Mus mus musculus and quantify the transition of craniofacial shape along the hybridization gradient. We find a continuous shape transition along the hybridization gradient and unaltered developmental stability associated with hybridization. This suggests that the morphospace between the two subspecies is continuous despite reproductive isolation and strong barriers to gene flow. We show that quantitative changes in overall genome composition generate quantitative changes in craniofacial shape; this supports a highly polygenic basis for between-species craniofacial differences in the house mouse. We discuss our findings in the context of oligogenic versus polygenic models of the genetic architecture of morphological traits.

Dissection and Flat-mounting of the Threespine Stickleback Branchial Skeleton

The posterior pharyngeal segments of the vertebrate head give rise to the branchial skeleton, the primary site of food processing in fish. The morphology of the fish branchial skeleton is matched to a species' diet. Threespine stickleback fish (Gasterosteus aculeatus) have emerged as a model system to study the genetic and developmental basis of evolved differences in a variety of traits. Marine populations of sticklebacks have repeatedly colonized countless new freshwater lakes and creeks. Adaptation to the new diet in these freshwater environments likely underlies a series of craniofacial changes that have evolved repeatedly in independently derived freshwater populations. These include three major patterning changes to the branchial skeleton: reductions in the number and length of gill raker bones, increases in pharyngeal tooth number, and increased branchial bone lengths. Here we describe a detailed protocol to dissect and flat-mount the internal branchial skeleton in threespine stickleback fish. Dissection of the entire three-dimensional branchial skeleton and mounting it flat into a largely two-dimensional prep allows for the easy visualization and quantification of branchial skeleton morphology. This dissection method is inexpensive, fast, relatively easy, and applicable to a wide variety of fish species. In sticklebacks, this efficient method allows the quantification of skeletal morphology in genetic crosses to map genomic regions controlling craniofacial patterning.

Goldfish morphology as a model for evolutionary developmental biology

Morphological variation of the goldfish is known to have been established by artificial selection for ornamental purposes during the domestication process. Chinese texts that date to the Song dynasty contain descriptions of goldfish breeding for ornamental purposes, indicating that the practice originated over one thousand years ago. Such a well-documented goldfish breeding process, combined with the phylogenetic and embryological proximities of this species with zebrafish, would appear to make the morphologically diverse goldfish strains suitable models for evolutionary developmental (evodevo) studies. However, few modern evodevo studies of goldfish have been conducted. In this review, we provide an overview of the historical background of goldfish breeding, and the differences between this teleost and zebrafish from an evolutionary perspective. We also summarize recent progress in the field of molecular developmental genetics, with a particular focus on the twin-tail goldfish morphology. Furthermore, we discuss unanswered questions relating to the evolution of the genome, developmental robustness, and morphologies in the goldfish lineage, with the goal of blazing a path toward an evodevo study paradigm using this teleost species as a new model species. For further resources related to this article, please visit the WIREs website.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.

Innate immune responses to gut microbiota differ between oceanic and freshwater threespine stickleback populations

Animal hosts must co-exist with beneficial microbes while simultaneously being able to mount rapid, non-specific, innate immune responses to pathogenic microbes. How this balance is achieved is not fully understood, and disruption of this relationship can lead to disease. Excessive inflammatory responses to resident microbes are characteristic of certain gastrointestinal pathologies such as inflammatory bowel disease (IBD). The immune dysregulation of IBD has complex genetic underpinnings that cannot be fully recapitulated with single-gene-knockout models. A deeper understanding of the genetic regulation of innate immune responses to resident microbes requires the ability to measure immune responses in the presence and absence of the microbiota using vertebrate models with complex genetic variation. Here, we describe a new gnotobiotic vertebrate model to explore the natural genetic variation that contributes to differences in innate immune responses to microbiota. Threespine stickleback, Gasterosteus aculeatus, has been used to study the developmental genetics of complex traits during the repeated evolution from ancestral oceanic to derived freshwater forms. We established methods to rear germ-free stickleback larvae and gnotobiotic animals monoassociated with single bacterial isolates. We characterized the innate immune response of these fish to resident gut microbes by quantifying the neutrophil cells in conventionally reared monoassociated or germ-free stickleback from both oceanic and freshwater populations grown in a common intermediate salinity environment. We found that oceanic and freshwater fish in the wild and in the laboratory share many intestinal microbial community members. However, oceanic fish mount a strong immune response to residential microbiota, whereas freshwater fish frequently do not. A strong innate immune response was uniformly observed across oceanic families, but this response varied among families of freshwater fish. The gnotobiotic stickleback model that we have developed therefore provides a platform for future studies mapping the natural genetic basis of the variation in immune response to microbes.

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.

The aggregate site frequency spectrum for comparative population genomic inference

Understanding how assemblages of species responded to past climate change is a central goal of comparative phylogeography and comparative population genomics, an endeavour that has increasing potential to integrate with community ecology. New sequencing technology now provides the potential to perform complex demographic inference at unprecedented resolution across assemblages of nonmodel species. To this end, we introduce the aggregate site frequency spectrum (aSFS), an expansion of the site frequency spectrum to use single nucleotide polymorphism (SNP) data sets collected from multiple, co-distributed species for assemblage-level demographic inference. We describe how the aSFS is constructed over an arbitrary number of independent population samples and then demonstrate how the aSFS can differentiate various multispecies demographic histories under a wide range of sampling configurations while allowing effective population sizes and expansion magnitudes to vary independently. We subsequently couple the aSFS with a hierarchical approximate Bayesian computation (hABC) framework to estimate degree of temporal synchronicity in expansion times across taxa, including an empirical demonstration with a data set consisting of five populations of the threespine stickleback (Gasterosteus aculeatus). Corroborating what is generally understood about the recent postglacial origins of these populations, the joint aSFS/hABC analysis strongly suggests that the stickleback data are most consistent with synchronous expansion after the Last Glacial Maximum (posterior probability = 0.99). The aSFS will have general application for multilevel statistical frameworks to test models involving assemblages and/or communities, and as large-scale SNP data from nonmodel species become routine, the aSFS expands the potential for powerful next-generation comparative population genomic inference.

Evolution of opercle shape in cichlid fishes from Lake Tanganyika - adaptive trait interactions in extant and extinct species flocks

Phenotype-environment correlations and the evolution of trait interactions in adaptive radiations have been widely studied to gain insight into the dynamics underpinning rapid species diversification. In this study we explore the phenotype-environment correlation and evolution of operculum shape in cichlid fishes using an outline-based geometric morphometric approach combined with stable isotope indicators of macrohabitat and trophic niche. We then apply our method to a sample of extinct saurichthyid fishes, a highly diverse and near globally distributed group of actinopterygians occurring throughout the Triassic, to assess the utility of extant data to inform our understanding of ecomorphological evolution in extinct species flocks. A series of comparative methods were used to analyze shape data for 54 extant species of cichlids (N = 416), and 6 extinct species of saurichthyids (N = 44). Results provide evidence for a relationship between operculum shape and feeding ecology, a concentration in shape evolution towards present along with evidence for convergence in form, and significant correlation between the major axes of shape change and measures of gut length and body elongation. The operculum is one of few features that can be compared in extant and extinct groups, enabling reconstruction of phenotype-environment interactions and modes of evolutionary diversification in deep time.

Regulation of Jaw Length During Development, Disease, and Evolution

Molecular and cellular mechanisms that control jaw length are becoming better understood. This is significant since the jaws are not only critical for species-specific adaptation and survival, but they are often affected by a variety of size-related anomalies including mandibular hypoplasia, retrognathia, asymmetry, and clefting. This chapter overviews how jaw length is established during the allocation, proliferation, differentiation, and growth of jaw precursor cells, which originate from neural crest mesenchyme (NCM). The focus is mainly on results from experiments transplanting NCM between quail and duck embryos. Quail have short jaws whereas those of duck are relatively long. Quail-duck chimeras reveal that the determinants of jaw length are NCM mediated throughout development and include species-specific differences in jaw progenitor number, differential regulation of various signaling pathways, and the autonomous activation of programs for skeletal matrix deposition and resorption. Such insights help make the goal of devising new therapies for birth defects, diseases, and injuries to the jaw skeleton seem ever more likely.

Integration of Random Forest with population-based outlier analyses provides insight on the genomic basis and evolution of run timing in Chinook salmon (Oncorhynchus tshawytscha)

Anadromous Chinook salmon populations vary in the period of river entry at the initiation of adult freshwater migration, facilitating optimal arrival at natal spawning. Run timing is a polygenic trait that shows evidence of rapid parallel evolution in some lineages, signifying a key role for this phenotype in the ecological divergence between populations. Studying the genetic basis of local adaptation in quantitative traits is often impractical in wild populations. Therefore, we used a novel approach, Random Forest, to detect markers linked to run timing across 14 populations from contrasting environments in the Columbia River and Puget Sound, USA. The approach permits detection of loci of small effect on the phenotype. Divergence between populations at these loci was then examined using both principle component analysis and FST outlier analyses, to determine whether shared genetic changes resulted in similar phenotypes across different lineages. Sequencing of 9107 RAD markers in 414 individuals identified 33 predictor loci explaining 79.2% of trait variance. Discriminant analysis of principal components of the predictors revealed both shared and unique evolutionary pathways in the trait across different lineages, characterized by minor allele frequency changes. However, genome mapping of predictor loci also identified positional overlap with two genomic outlier regions, consistent with selection on loci of large effect. Therefore, the results suggest selective sweeps on few loci and minor changes in loci that were detected by this study. Use of a polygenic framework has provided initial insight into how divergence in a trait has occurred in the wild.

Modular skeletal evolution in sticklebacks is controlled by additive and clustered quantitative trait Loci

Understanding the genetic architecture of evolutionary change remains a long-standing goal in biology. In vertebrates, skeletal evolution has contributed greatly to adaptation in body form and function in response to changing ecological variables like diet and predation. Here we use genome-wide linkage mapping in threespine stickleback fish to investigate the genetic architecture of evolved changes in many armor and trophic traits. We identify >100 quantitative trait loci (QTL) controlling the pattern of serially repeating skeletal elements, including gill rakers, teeth, branchial bones, jaws, median fin spines, and vertebrae. We use this large collection of QTL to address long-standing questions about the anatomical specificity, genetic dominance, and genomic clustering of loci controlling skeletal differences in evolving populations. We find that most QTL (76%) that influence serially repeating skeletal elements have anatomically regional effects. In addition, most QTL (71%) have at least partially additive effects, regardless of whether the QTL controls evolved loss or gain of skeletal elements. Finally, many QTL with high LOD scores cluster on chromosomes 4, 20, and 21. These results identify a modular system that can control highly specific aspects of skeletal form. Because of the general additivity and genomic clustering of major QTL, concerted changes in both protective armor and trophic traits may occur when sticklebacks inherit either marine or freshwater alleles at linked or possible "supergene" regions of the stickleback genome. Further study of these regions will help identify the molecular basis of both modular and coordinated changes in the vertebrate skeleton.

Identification of major and minor QTL for ecologically important morphological traits in three-spined sticklebacks (Gasterosteus aculeatus)

Quantitative trait locus (QTL) mapping studies of Pacific three-spined sticklebacks (Gasterosteus aculeatus) have uncovered several genomic regions controlling variability in different morphological traits, but QTL studies of Atlantic sticklebacks are lacking. We mapped QTL for 40 morphological traits, including body size, body shape, and body armor, in a F₂ full-sib cross between northern European marine and freshwater three-spined sticklebacks. A total of 52 significant QTL were identified at the 5% genome-wide level. One major QTL explaining 74.4% of the total variance in lateral plate number was detected on LG4, whereas several major QTL for centroid size (a proxy for body size), and the lengths of two dorsal spines, pelvic spine, and pelvic girdle were mapped on LG21 with the explained variance ranging from 27.9% to 57.6%. Major QTL for landmark coordinates defining body shape variation also were identified on LG21, with each explaining ≥15% of variance in body shape. Multiple QTL for different traits mapped on LG21 overlapped each other, implying pleiotropy and/or tight linkage. Thus, apart from providing confirmatory data to support conclusions born out of earlier QTL studies of Pacific sticklebacks, this study also describes several novel QTL of both major and smaller effect for ecologically important traits. The finding that many major QTL mapped on LG21 suggests that this linkage group might be a hotspot for genetic determinants of ecologically important morphological traits in three-spined sticklebacks.

Complex craniofacial changes in blind cave-dwelling fish are mediated by genetically symmetric and asymmetric loci

The genetic regulators of regressive craniofacial morphologies are poorly understood. To shed light on this problem, we examined the freshwater fish Astyanax mexicanus, a species with surface-dwelling and multiple independent eyeless cave-dwelling forms. Changes affecting the skull in cavefish include morphological alterations to the intramembranous circumorbital bones encircling the eye. Many of these modifications, however, have evolved separately from eye loss, such as fragmentation of the third suborbital bone. To understand the genetic architecture of these eye-independent craniofacial alterations, we developed and scored 33 phenotypes in the context of an F2 hybrid mapping pedigree bred from Pachón cavefish and surface fish. We discovered several individuals exhibiting dramatic left-right differences in bone formation, such as extensive fragmentation on the right side only. This observation, along with well-known eye size asymmetry in natural cave-dwelling animals, led us to further evaluate left-right genetic differences for the craniofacial complex. We discovered three phenotypes, inclusive of bone fragmentation and fusion, which demonstrated a directional heritable basis only on one side. Interestingly, the overall areas of affected bones were genetically symmetric. Phenotypic effect plots of these novel craniofacial QTL revealed that cave alleles are associated with abnormal conditions such as bony fusion and fragmentation. Moreover, many linked loci overlapped with other cave-associated traits, suggesting regressive craniofacial changes may evolve through linkage or as antagonistic pleiotropic consequences of cave-associated adaptations. These novel findings illuminate significant craniofacial changes accompanying evolution in complete darkness and reveal complex changes to the skull differentially influenced by genetic changes affecting the left and right sides.

Genomic divergence between nine- and three-spined sticklebacks

Background

Comparative genomics approaches help to shed light on evolutionary processes that shape differentiation between lineages. The nine-spined stickleback (Pungitius pungitius) is a closely related species of the ecological ‘supermodel’ three-spined stickleback (Gasterosteus aculeatus). It is an emerging model system for evolutionary biology research but has garnered less attention and lacks extensive genomic resources. To expand on these resources and aid the study of sticklebacks in a phylogenetic framework, we characterized nine-spined stickleback transcriptomes from brain and liver using deep sequencing.

Results

We obtained nearly eight thousand assembled transcripts, of which 3,091 were assigned as putative one-to-one orthologs to genes found in the three-spined stickleback. These sequences were used for evaluating overall differentiation and substitution rates between nine- and three-spined sticklebacks, and to identify genes that are putatively evolving under positive selection. The synonymous substitution rate was estimated to be 7.1 × 10^-9 per site per year between the two species, and a total of 165 genes showed patterns of adaptive evolution in one or both species. A few nine-spined stickleback contigs lacked an obvious ortholog in three-spined sticklebacks but were found to match genes in other fish species, suggesting several gene losses within 13 million years since the divergence of the two stickleback species. We identified 47 SNPs in 25 different genes that differentiate pond and marine ecotypes. We also identified 468 microsatellites that could be further developed as genetic markers in nine-spined sticklebacks.

Conclusion

With deep sequencing of nine-spined stickleback cDNA libraries, our study provides a significant increase in the number of gene sequences and microsatellite markers for this species, and identifies a number of genes showing patterns of adaptive evolution between nine- and three-spined sticklebacks. We also report several candidate genes that might be involved in differential adaptation between marine and freshwater nine-spined sticklebacks. This study provides a valuable resource for future studies aiming to identify candidate genes underlying ecological adaptation in this and other stickleback species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-756) contains supplementary material, which is available to authorized users.

Embryonic development of goldfish (Carassius auratus): a model for the study of evolutionary change in developmental mechanisms by artificial selection

BACKGROUND

Highly divergent morphology among the different goldfish strains (Carassius auratus) may make it a suitable model for investigating how artificial selection has altered developmental mechanisms. Here we describe the embryological development of the common goldfish (the single fin Wakin), which retains the ancestral morphology of this species.

RESULTS

We divided goldfish embryonic development into seven periods consisting of 34 stages, using previously reported developmental indices of zebrafish and goldfish. Although several differences were identified in terms of their yolk size, epiboly process, pigmentation patterns, and development rate, our results indicate that the embryonic features of these two teleost species are highly similar in their overall morphology from the zygote to hatching stage.

CONCLUSIONS

These results provide an opportunity for further study of the evolutionary relationship between domestication and development, through applying well-established zebrafish molecular biological resources to goldfish embryos.

Ecomorphological disparity in an adaptive radiation: opercular bone shape and stable isotopes in Antarctic icefishes

To assess how ecological and morphological disparity is interrelated in the adaptive radiation of Antarctic notothenioid fish we used patterns of opercle bone evolution as a model to quantify shape disparity, phylogenetic patterns of shape evolution, and ecological correlates in the form of stable isotope values. Using a sample of 25 species including representatives from four major notothenioid clades, we show that opercle shape disparity is higher in the modern fauna than would be expected under the neutral evolution Brownian motion model. Phylogenetic comparative methods indicate that opercle shape data best fit a model of directional selection (Ornstein–Uhlenbeck) and are least supported by the “early burst” model of adaptive radiation. The main evolutionary axis of opercle shape change reflects movement from a broad and more symmetrically tapered opercle to one that narrows along the distal margin, but with only slight shape change on the proximal margin. We find a trend in opercle shape change along the benthic–pelagic axis, underlining the importance of this axis for diversification in the notothenioid radiation. A major impetus for the study of adaptive radiations is to uncover generalized patterns among different groups, and the evolutionary patterns in opercle shape among notothenioids are similar to those found among other adaptive radiations (three-spined sticklebacks) promoting the utility of this approach for assessing ecomorphological interactions on a broad scale.

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.

Convergent evolution as a generator of phenotypic diversity in threespine stickleback

Convergent evolution, in which populations produce similar phenotypes in response to similar selection pressure, is strong evidence for the role of natural selection in shaping biological diversity. In some cases, closely related populations can produce functionally similar but phenotypically divergence forms in response to selection. Functional convergence with morphological divergence has been observed in laboratory selection experiments and computer simulations, but while potentially common, is rarely recognized in nature. Here, we present data from the North Pacific threespine stickleback radiation showing that ecologically and functionally similar, but morphologically divergent phenotypes rapidly evolved when an ancestral population colonized freshwater benthic habitats in parallel. In addition, we show that in this system, functional convergence substantially increases morphospace occupation relative to ancestral phenotypes, which suggests that convergent evolution may, paradoxically, be an important and previously underappreciated source of morphological diversity.

Developmental dissociation in morphological evolution of the stickleback opercle

Oceanic threespine sticklebacks have repeatedly and independently evolved new morphologies upon invasions of freshwater habitats. A consistent derived feature of the freshwater form across populations and geography is a shape change of the opercle, a large early developing facial bone. We show that the principal multivariate axis describing opercle shape development from the young larva to the full adult stage of oceanic fish matches the principal axis of evolutionary change associated with relocation from the oceanic to freshwater habitat. The opercle phenotype of freshwater adults closely resembles the phenotype of the bone in juveniles. Thus, evolution to the freshwater condition is in large part by truncation of development; the freshwater fish do not achieve the full ancestral adult bone shape. Additionally, the derived state includes dissociated ontogenetic changes. Dissociability may reflect an underlying modular pattern of opercle development, and facilitate flexibility of morphological evolution.

Genetic architecture of variation in the lateral line sensory system of threespine sticklebacks

Vertebrate sensory systems have evolved remarkable diversity, but little is known about the underlying genetic mechanisms. The lateral line sensory system of aquatic vertebrates is a promising model for genetic investigations of sensory evolution because there is extensive variation within and between species, and this variation is easily quantified. In the present study, we compare the lateral line sensory system of threespine sticklebacks (Gasterosteus aculeatus) from an ancestral marine and a derived benthic lake population. We show that lab-raised individuals from these populations display differences in sensory neuromast number, neuromast patterning, and groove morphology. Using genetic linkage mapping, we identify regions of the genome that influence different aspects of lateral line morphology. Distinct loci independently affect neuromast number on different body regions, suggesting that a modular genetic structure underlies the evolution of peripheral receptor number in this sensory system. Pleiotropy and/or tight linkage are also important, as we identify a region on linkage group 21 that affects multiple aspects of lateral line morphology. Finally, we detect epistasis between a locus on linkage group 4 and a locus on linkage group 21; interactions between these loci contribute to variation in neuromast pattern. Our results reveal a complex genetic architecture underlying the evolution of the stickleback lateral line sensory system. This study further uncovers a genetic relationship between sensory morphology and non-neural traits (bony lateral plates), creating an opportunity to investigate morphological constraints on sensory evolution in a vertebrate model system.

Admixture mapping of male nuptial colour and body shape in a recently formed hybrid population of threespine stickleback

Despite recent progress, we still know relatively little about the genetic architecture that underlies adaptation to divergent environments. Determining whether the genetic architecture of phenotypic adaptation follows any predictable patterns requires data from a wide variety of species. However, in many organisms, genetic studies are hindered by the inability to perform genetic crosses in the laboratory or by long generation times. Admixture mapping is an approach that circumvents these issues by taking advantage of hybridization that occurs between populations or species in the wild. Here, we demonstrate the utility of admixture mapping in a naturally occurring hybrid population of threespine sticklebacks (Gasterosteus aculeatus) from Enos Lake, British Columbia. Until recently, this lake contained two species of sticklebacks adapted to divergent habitats within the lake. This benthic-limnetic species pair diverged in a number of phenotypes, including male nuptial coloration and body shape, which were previously shown to contribute to reproductive isolation between them. However, recent ecological disturbance has contributed to extensive hybridization between the species, and there is now a single, admixed population within Enos Lake. We collected over 500 males from Enos Lake and found that most had intermediate nuptial colour and body shape. By genotyping males with nuptial colour at the two extremes of the phenotypic distribution, we identified seven genomic regions on three chromosomes associated with divergence in male nuptial colour. These genomic regions are also associated with variation in body shape, suggesting that tight linkage and/or pleiotropy facilitated adaptation to divergent environments in benthic-limnetic species pairs.

Extensive linkage disequilibrium and parallel adaptive divergence across threespine stickleback genomes

Population genomic studies are beginning to provide a more comprehensive view of dynamic genome-scale processes in evolution. Patterns of genomic architecture, such as genomic islands of increased divergence, may be important for adaptive population differentiation and speciation. We used next-generation sequencing data to examine the patterns of local and long-distance linkage disequilibrium (LD) across oceanic and freshwater populations of threespine stickleback, a useful model for studies of evolution and speciation. We looked for associations between LD and signatures of divergent selection, and assessed the role of recombination rate variation in generating LD patterns. As predicted under the traditional biogeographic model of unidirectional gene flow from ancestral oceanic to derived freshwater stickleback populations, we found extensive local and long-distance LD in fresh water. Surprisingly, oceanic populations showed similar patterns of elevated LD, notably between large genomic regions previously implicated in adaptation to fresh water. These results support an alternative biogeographic model for the stickleback radiation, one of a metapopulation with appreciable bi-directional gene flow combined with strong divergent selection between oceanic and freshwater populations. As predicted by theory, these processes can maintain LD within and among genomic islands of divergence. These findings suggest that the genomic architecture in oceanic stickleback populations may provide a mechanism for the rapid re-assembly and evolution of multi-locus genotypes in newly colonized freshwater habitats, and may help explain genetic mapping of parallel phenotypic variation to similar loci across independent freshwater populations.