Search

A critical question in biology is how new traits evolve, but studying this in wild animals remains challenging. Here, we probe the genetic basis of trait gain in sea robin fish, which have evolved specialized leg-like appendages for locomotion and digging along the ocean floor. We use genome sequencing, transcriptional profiling, and interspecific hybrid analysis to explore the molecular and developmental basis of leg formation. We identified the ancient, conserved transcription factor tbx3a as a major determinant of sensory leg development. Genome editing confirms that tbx3a is required for normal leg formation in sea robins, and for formation of enlarged central nervous system lobes, sensory papillae, and adult digging behavior. Our study establishes sea robins as a model organism for studying the evolution of major trait gain and illustrates how ancient developmental control genes can underlie novel organ formation., Competing Interests: Declaration of interests The authors declare that they have no competing interests., (Copyright © 2024 The Author(s). Published by Elsevier Inc. All rights reserved.)

How do animals evolve new traits? Sea robins are fish that possess specialized leg-like appendages used to "walk" along the sea floor. Here, we show that legs are bona fide sense organs that localize buried prey. Legs are covered in sensory papillae that receive dense innervation from touch-sensitive neurons, express non-canonical epithelial taste receptors, and mediate chemical sensitivity that drives predatory digging behavior. A combination of developmental analyses, crosses between species with and without papillae, and interspecies comparisons of sea robins from around the world demonstrate that papillae represent a key evolutionary innovation associated with behavioral niche expansion on the sea floor. These discoveries provide unique insight into how molecular-, cellular-, and tissue-scale adaptations integrate to produce novel organismic traits and behavior., Competing Interests: Declaration of interests The authors declare that they have no competing interests., (Copyright © 2024 Elsevier Inc. All rights reserved.)

Despite deep evolutionary conservation, recombination rates vary greatly across the genome and among individuals, sexes and populations. Yet the impact of this variation on adaptively diverging populations is not well understood. Here we characterized fine-scale recombination landscapes in an adaptively divergent pair of marine and freshwater populations of threespine stickleback from River Tyne, Scotland. Through whole-genome sequencing of large nuclear families, we identified the genomic locations of almost 50,000 crossovers and built recombination maps for marine, freshwater and hybrid individuals at a resolution of 3.8 kb. We used these maps to quantify the factors driving variation in recombination rates. We found strong heterochiasmy between sexes but also differences in recombination rates among ecotypes. Hybrids showed evidence of significant recombination suppression in overall map length and in individual loci. Recombination rates were lower not only within individual marine-freshwater-adaptive loci, but also between loci on the same chromosome, suggesting selection on linked gene 'cassettes'. Through temporal sampling along a natural hybrid zone, we found that recombinants showed traits associated with reduced fitness. Our results support predictions that divergence in cis-acting recombination modifiers, whose functions are disrupted in hybrids, may play an important role in maintaining differences among adaptively diverging populations., (© 2024. The Author(s).)

Although promoters and their enhancers are frequently contained within a topologically associating domain (TAD), some developmentally important genes have their promoter and enhancers within different TADs. Hypotheses about molecular mechanisms enabling cross-TAD interactions remain to be assessed. To test these hypotheses, we used optical reconstruction of chromatin architecture to characterize the conformations of the Pitx1 locus on single chromosomes in developing mouse limbs. Our data support a model in which neighboring boundaries are stacked as a result of loop extrusion, bringing boundary-proximal cis-elements into contact. This stacking interaction also contributes to the appearance of architectural stripes in the population average maps. Through molecular dynamics simulations, we found that increasing boundary strengths facilitates the formation of the stacked boundary conformation, counter-intuitively facilitating border bypass. This work provides a revised view of the TAD borders' function, both facilitating and preventing cis-regulatory interactions, and introduces a framework to distinguish border-crossing from border-respecting enhancer-promoter pairs., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

Fins are major functional appendages of fish that have been repeatedly modified in different lineages. To search for genomic changes underlying natural fin diversity, we compared the genomes of 36 percomorph fish species that span over 100 million years of evolution and either have complete or reduced pelvic and caudal fins. We identify 1,614 genomic regions that are well-conserved in fin-complete species but missing from multiple fin-reduced lineages. Recurrent deletions of conserved sequences in wild fin-reduced species are enriched for functions related to appendage development, suggesting that convergent fin reduction at the organismal level is associated with repeated genomic deletions near fin-appendage development genes. We used sequencing and functional enhancer assays to confirm that PelA, a Pitx1 enhancer previously linked to recurrent pelvic loss in sticklebacks, has also been independently deleted and may have contributed to the fin morphology in distantly related pelvic-reduced species. We also identify a novel enhancer that is conserved in the majority of percomorphs, drives caudal fin expression in transgenic stickleback, is missing in tetraodontiform, syngnathid, and synbranchid species with caudal fin reduction, and alters caudal fin development when targeted by genome editing. Our study illustrates a broadly applicable strategy for mapping phenotypes to genotypes across a tree of vertebrate species and highlights notable new examples of regulatory genomic hotspots that have been used to evolve recurrent phenotypes across 100 million years of fish evolution., Competing Interests: Conflict of interest statement. H.I.C., Y.T., G.B., and D.M.K. have no conflicts of interest to declare., (© The Author(s) 2023. Published by Oxford University Press on behalf of Society for Molecular Biology and Evolution.)

Allard et al. provide an overview of sea robins, a group of benthic fish that have evolved leg-like appendages that they use to walk on the sea floor and find prey., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Sophorolipids are biosurfactants derived from the nonpathogenic yeasts such as Starmerella bombicola with potential efficacy in anticancer applications. Simple and cost-effective synthesis of these drugs makes them a promising alternative to traditional chemotherapeutics, pending their success in preliminary drug-screening. Drug-screening typically utilizes 2D cell monolayers due to their simplicity and ease of high-throughput assessment. However, 2D assays fail to capture the complexity and 3D context of the tumor microenvironment and have consequently been implicated in the high percentage of drugs investigated in vitro that later fail in clinical trials. Herein, we screened two sophorolipid candidates and a clinically-used chemotherapeutic, doxorubicin, on in vitro breast cancer models ranging from 2D monolayers to 3D spheroids, employing optical coherence tomography to confirm these morphologies. We calculated corresponding IC 50 values for these drugs and found one of the sophorolipids to have comparable toxicities to the chemotherapeutic control. Our findings show increased drug resistance associated with model dimensionality, such that all drugs tested showed that 3D spheroids exhibited higher IC 50 values than their 2D counterparts. These findings demonstrate promising preliminary data to support the use of sophorolipids as a more affordable alternative to traditional clinical interventions and demonstrate the importance of 3D tumor models in assessing drug response.

Hydrogel microbeads are engineered spherical microgels widely used for biomedical applications in cell cultures, tissue engineering, and drug delivery. Their mechanical and physical properties (i.e., modulus, porosity, diffusion) heavily influence their utility by affecting encapsulated cellular behavior, biopayload elution kinetics, and stability for longer term cultures. There is a need to quantify these properties to guide microbead design for effective application. However, there are few techniques with the μN-level resolution required to evaluate these relatively small, compliant constructs. To circumvent mechanically testing individual microbeads, researchers often approximate microbead properties by characterizing larger bulk gel analogs of the same material formulation. This approach provides some insight into the hydrogel properties. However, bulk gels possess key structural and mechanical differences compared to their microbead equivalents, which may limit their accuracy and utility as analogs for estimating microbead properties. Herein, we explore how microbead properties are influenced by hydrogel formulation (i.e., alginate concentration, divalent cation crosslinker, and crosslinker concentration), and whether these trends are accurately reflected in bulk gel analogs. To accomplish this, we utilize laser direct-write bioprinting to create 12 × 12 arrays of alginate microbeads and characterize all 144 microbeads in parallel using a commercially available microcompression system. In this way, the compressive load is distributed across a large number of beads, thus amplifying sample signal. Comparing microbead properties to those of their bulk gel analogs, we found that their trends in modulus, porosity, and diffusion with hydrogel formulation are consistent, yet bulk gels exhibit significant discrepancies in their measured values., (Copyright © 2023 by ASME.)

Understanding the mechanisms leading to new traits or additional features in organisms is a fundamental goal of evolutionary biology. We show that HOXDB regulatory changes have been used repeatedly in different fish genera to alter the length and number of the prominent dorsal spines used to classify stickleback species. In Gasterosteus aculeatus (typically 'three-spine sticklebacks'), a variant HOXDB allele is genetically linked to shortening an existing spine and adding an additional spine. In Apeltes quadracus (typically 'four-spine sticklebacks'), a variant HOXDB allele is associated with lengthening a spine and adding an additional spine in natural populations. The variant alleles alter the same non-coding enhancer region in the HOXDB locus but do so by diverse mechanisms, including single-nucleotide polymorphisms, deletions and transposable element insertions. The independent regulatory changes are linked to anterior expansion or contraction of HOXDB expression. We propose that associated changes in spine lengths and numbers are partial identity transformations in a repeating skeletal series that forms major defensive structures in fish. Our findings support the long-standing hypothesis that natural Hox gene variation underlies key patterning changes in wild populations and illustrate how different mutational mechanisms affecting the same region may produce opposite gene expression changes with similar phenotypic outcomes., (© 2022. The Author(s).)

In adaptive radiations, single lineages rapidly diversify by adapting to many new niches. Little is known yet about the genomic mechanisms involved, that is, the source of genetic variation or genomic architecture facilitating or constraining adaptive radiation. Here, we investigate genomic changes associated with repeated invasion of many different freshwater niches by threespine stickleback in the Haida Gwaii archipelago, Canada, by resequencing single genomes from one marine and 28 freshwater populations. We find 89 likely targets of parallel selection in the genome that are enriched for old standing genetic variation. In contrast to theoretical expectations, their genomic architecture is highly dispersed with little clustering. Candidate genes and genotype-environment correlations match the three major environmental axes predation regime, light environment, and ecosystem size. In a niche space with these three dimensions, we find that the more divergent a new niche from the ancestral marine habitat, the more loci show signatures of parallel selection. Our findings suggest that the genomic architecture of parallel adaptation in adaptive radiation depends on the steepness of ecological gradients and the dimensionality of the niche space., (© 2022 The Authors. Evolution published by Wiley Periodicals LLC on behalf of The Society for the Study of Evolution.)

Many key signaling molecules used to build tissues during embryonic development are re-activated at injury sites to stimulate tissue regeneration and repair. Bone morphogenetic proteins provide a classic example, but the mechanisms that lead to reactivation of BMPs following injury are still unknown. Previous studies have mapped a large "injury response element" (IRE) in the mouse Bmp5 gene that drives gene expression following bone fractures and other types of injury. Here we show that the large mouse IRE region is also activated in both zebrafish tail resection and mechanosensory hair cell injury models. Using the ability to test multiple constructs and image temporal and spatial dynamics following injury responses, we have narrowed the original size of the mouse IRE region by over 100 fold and identified a small 142 bp minimal enhancer that is rapidly induced in both mesenchymal and epithelial tissues after injury. These studies identify a small sequence that responds to evolutionarily conserved local signals in wounded tissues and suggest candidate pathways that contribute to BMP reactivation after injury., (Copyright © 2021. Published by Elsevier Inc.)

Complete genome sequencing has identified millions of DNA changes that differ between humans and chimpanzees. Although a subset of these changes likely underlies important phenotypic differences between humans and chimpanzees, it is currently difficult to distinguish causal from incidental changes and to map specific phenotypes to particular genome locations. To facilitate further genetic study of human-chimpanzee divergence, we have generated human and chimpanzee autotetraploids and allotetraploids by fusing induced pluripotent stem cells (iPSCs) of each species. The resulting tetraploid iPSCs can be stably maintained and retain the ability to differentiate along ectoderm, mesoderm, and endoderm lineages. RNA sequencing identifies thousands of genes whose expression differs between humans and chimpanzees when assessed in single-species diploid or autotetraploid iPSCs. Analysis of gene expression patterns in interspecific allotetraploid iPSCs shows that human-chimpanzee expression differences arise from substantial contributions of both cis -acting changes linked to the genes themselves and trans -acting changes elsewhere in the genome. To enable further genetic mapping of species differences, we tested chemical treatments for stimulating genome-wide mitotic recombination between human and chimpanzee chromosomes, and CRISPR methods for inducing species-specific changes on particular chromosomes in allotetraploid cells. We successfully generated derivative cells with nested deletions or interspecific recombination on the X chromosome. These studies confirm an important role for the X chromosome in trans regulation of expression differences between species and illustrate the potential of this system for more detailed cis and trans mapping of the molecular basis of human and chimpanzee evolution., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

Vertebrates have repeatedly modified skeletal structures to adapt to their environments. The threespine stickleback is an excellent system for studying skeletal modifications, as different wild populations have either increased or decreased the lengths of their prominent dorsal and pelvic spines in different freshwater environments. Here we identify a regulatory locus that has a major morphological effect on the length of stickleback dorsal and pelvic spines, which we term Maser (major spine enhancer). Maser maps in a closely linked supergene complex that controls multiple armor, feeding, and behavioral traits on chromosome IV. Natural alleles in Maser are differentiated between marine and freshwater sticklebacks; however, alleles found among freshwater populations are also differentiated, with distinct alleles found in short- and long-spined freshwater populations. The distinct freshwater alleles either increase or decrease expression of the bone growth inhibitor gene Stanniocalcin2a in developing spines, providing a simple genetic mechanism for either increasing or decreasing spine lengths in natural populations. Genomic surveys suggest many recurrently differentiated loci in sticklebacks are similarly specialized into three or more distinct alleles, providing multiple ancient standing variants in particular genes that may contribute to a range of phenotypes in different environments., Competing Interests: The authors declare no competing interest., (Copyright © 2021 the Author(s). Published by PNAS.)

Similar forms often evolve repeatedly in nature, raising long-standing questions about the underlying mechanisms. Here, we use repeated evolution in stickleback to identify a large set of genomic loci that change recurrently during colonization of freshwater habitats by marine fish. The same loci used repeatedly in extant populations also show rapid allele frequency changes when new freshwater populations are experimentally established from marine ancestors. Marked genotypic and phenotypic changes arise within 5 years, facilitated by standing genetic variation and linkage between adaptive regions. Both the speed and location of changes can be predicted using empirical observations of recurrence in natural populations or fundamental genomic features like allelic age, recombination rates, density of divergent loci, and overlap with mapped traits. A composite model trained on these stickleback features can also predict the location of key evolutionary loci in Darwin's finches, suggesting that similar features are important for evolution across diverse taxa., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).)

Mutations of small effect underlie most adaptation to new environments, but beneficial variants with large fitness effects are expected to contribute under certain conditions. Genes and genomic regions having large effects on phenotypic differences between populations are known from numerous taxa, but fitness effect sizes have rarely been estimated. We mapped fitness over a generation in an F2 intercross between a marine and a lake stickleback population introduced to a freshwater pond. A quantitative trait locus map of the number of surviving offspring per F2 female detected a single, large-effect locus near Ectodysplasin ( Eda ), a gene having an ancient freshwater allele causing reduced bony armor and other changes. F2 females homozygous for the freshwater allele had twice the number of surviving offspring as homozygotes for the marine allele, producing a large selection coefficient, s = 0.50 ± 0.09 SE. Correspondingly, the frequency of the freshwater allele increased from 0.50 in F2 mothers to 0.58 in surviving offspring. We compare these results to allele frequency changes at the Eda gene in an Alaskan lake population colonized by marine stickleback in the 1980s. The frequency of the freshwater Eda allele rose steadily over multiple generations and reached 95% within 20 y, yielding a similar estimate of selection, s = 0.49 ± 0.05, but a different degree of dominance. These findings are consistent with other studies suggesting strong selection on this gene (and/or linked genes) in fresh water. Selection on ancient genetic variants carried by colonizing ancestors is likely to increase the prevalence of large-effect fitness variants in adaptive evolution., Competing Interests: The authors declare no competing interest.

Three-dimensional (3D) tissue-engineered in vitro models, particularly multicellular spheroids and organoids, have become important tools to explore disease progression and guide the development of novel therapeutic strategies. These avascular constructs are particularly powerful in oncological research due to their ability to mimic several key aspects of in vivo tumors, such as 3D structure and pathophysiologic gradients. Advancement of spheroid models requires characterization of critical features (i.e., size, shape, cellular density, and viability) during model development, and in response to treatment. However, evaluation of these characteristics longitudinally, quantitatively and non-invasively remains a challenge. Herein, Optical Coherence Tomography (OCT) is used as a label-free tool to assess 3D morphologies and cellular densities of tumor spheroids generated via the liquid overlay technique. We utilize this quantitative tool to assess Matrigel's influence on spheroid morphologic development, finding that the absence of Matrigel produces flattened, disk-like aggregates rather than 3D spheroids with physiologically-relevant features. Furthermore, this technology is adapted to quantify cell number within tumor spheroids, and to discern between live and dead cells, to non-destructively provide valuable information on tissue/construct viability, as well as a proof-of-concept for longitudinal drug efficacy studies. Together, these findings demonstrate OCT as a promising noninvasive, quantitative, label-free, longitudinal and cell-based method that can assess development and drug response in 3D cellular aggregates at a mesoscopic scale., Competing Interests: Declaration of Competing Interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2020. Published by Elsevier Ltd.)

Background: Heteromorphic sex chromosomes have evolved repeatedly across diverse species. Suppression of recombination between X and Y chromosomes leads to degeneration of the Y chromosome. The progression of degeneration is not well understood, as complete sequence assemblies of heteromorphic Y chromosomes have only been generated across a handful of taxa with highly degenerate sex chromosomes. Here, we describe the assembly of the threespine stickleback (Gasterosteus aculeatus) Y chromosome, which is less than 26 million years old and at an intermediate stage of degeneration. Our previous work identified that the non-recombining region between the X and the Y spans approximately 17.5 Mb on the X chromosome., Results: We combine long-read sequencing with a Hi-C-based proximity guided assembly to generate a 15.87 Mb assembly of the Y chromosome. Our assembly is concordant with cytogenetic maps and Sanger sequences of over 90 Y chromosome BAC clones. We find three evolutionary strata on the Y chromosome, consistent with the three inversions identified by our previous cytogenetic analyses. The threespine stickleback Y shows convergence with more degenerate sex chromosomes in the retention of haploinsufficient genes and the accumulation of genes with testis-biased expression, many of which are recent duplicates. However, we find no evidence for large amplicons identified in other sex chromosome systems. We also report an excellent candidate for the master sex-determination gene: a translocated copy of Amh (Amhy)., Conclusions: Together, our work shows that the evolutionary forces shaping sex chromosomes can cause relatively rapid changes in the overall genetic architecture of Y chromosomes.

3D multicellular aggregates, and more advanced organotypic systems, have become central tools in recent years to study a wide variety of complex biological processes. Most notably, these model systems have become mainstream within oncology (multicellular tumor spheroids) and regenerative medicine (embryoid bodies) research. However, the biological behavior of these in vitro tissue surrogates is extremely sensitive to their aggregate size and geometry. Indeed, both of these geometrical parameters are key in producing pathophysiological gradients responsible for cellular and structural heterogeneity, replicating in vivo observations. Moreover, the fabrication techniques most widely used for producing these models lack the ability to accurately control cellular spatial location, an essential component for regulating homotypic and heterotypic cell signaling. Herein, we report on a 3D bioprinting technique, laser direct-write (LDW), that enables precise control of both spatial patterning and size of cell-encapsulating microbeads. The generated cell-laden beads are further processed into core-shelled structures, allowing for the growth and formation of self-contained, self-aggregating cells (e.g., breast cancer cells, embryonic stem cells). Within these structures we demonstrate our ability to produce multicellular tumor spheroids (MCTSs) and embryoid bodies (EBs) with well-controlled overall size and shape, that can be designed on demand. Furthermore, we investigated the impact of aggregate size on the uptake of a commonly employed ligand for receptor-mediated drug delivery, Transferrin, indicating that larger tumor spheroids exhibit greater spatial heterogeneity in ligand uptake. Taken together, these findings establish LDW as a versatile biomanufacturing platform for bioprinting and patterning core-shelled structures to generate size-controlled 3D multicellular aggregates. STATEMENT OF SIGNIFICANCE: Multicellular 3D aggregates are powerful in vitro models used to study a wide variety of complex biological processes, particularly within oncology and regenerative medicine. These tissue surrogates are fabricated using environments that encourage cellular self-assembly. However, specific applications require control of aggregate size and position to recapitulate key in vivo parameters (e.g., pathophysiological gradients and homotypic/heterotypic cell signaling). Herein, we demonstrate the ability to create and spatially pattern size-controlled embryoid bodies and tumor spheroids, using laser-based 3D bioprinting. Furthermore, we investigated the effect of tumor spheroid size on internalization of Transferrin, a common ligand for targeted therapy, finding greater spatial heterogeneity in our large aggregates. Overall, this technique offers incredible promise and flexibility for fabricating idealized 3D in vitro models., (Copyright © 2019 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.)

We present an electrodeposition technique for fabricating tubular alginate structures. In this technique, two electrodes (anode and cathode) are suspended in a solution of alginate and insoluble calcium carbonate particles, and the application of an electrical potential produces a localized pH change at the anode surface causing suspended divalent cations to become soluble and cross-link the alginate. We robustly characterize how the fabrication parameters influence the rate of radial deposition on the anode, including deposition time, applied voltage, alginate concentration, type of divalent cation and concentration, and anode diameter. Furthermore, we produce gels with a range of tailorable features, including mechanical properties, dimensions (thick-ness and lumen size), customizable tubular geometries, and radial compositional heterogeneity., Competing Interests: The authors declare no competing financial interest.

Demand for materials that mechanically replicate native tissue has driven development and characterization of various new biomaterials. However, a consequence of materials and characterization technique diversity is a lack of consensus within the field, with no clear way to compare values measured via different modalities. This likely contributes to the difficulty in replicating findings across the research community; recent evidence suggests that different modalities do not yield the same mechanical measurements within a material, and direct comparisons cannot be made across different testing platforms. Herein, we examine whether "material properties" are characterization modality-specific by analyzing the elastic moduli determined by five typical biomaterial mechanical characterization techniques: unconfined-compression, tensiometry, rheometry, and micro-indentation at the macroscopic level, and microscopically using nanoindentation. These analyses were performed in two different polymeric gels frequently used for biological applications, polydimethylsiloxane (PDMS) and agarose. Each was fabricated to span a range of moduli, from physiologic to supraphysiologic values. All five techniques identified the same overall trend within each material group, supporting their ability to appreciate relative moduli differences. However, significant differences were found across modalities, illustrating a difference in absolute moduli values, and thereby precluding direct comparison of measurements from different characterization modalities. These observed differences may depend on material compliance, viscoelasticity, and microstructure. While determining the underlying mechanism(s) of these differences was beyond the scope of this work, these results demonstrate how each modality affects the measured moduli of the same material, and the sensitivity of each modality to changes in sample material composition., (Copyright © 2019 Elsevier Ltd. All rights reserved.)

Evolution generates a remarkable breadth of living forms, but many traits evolve repeatedly, by mechanisms that are still poorly understood. A classic example of repeated evolution is the loss of pelvic hindfins in stickleback fish ( Gasterosteus aculeatus ). Repeated pelvic loss maps to recurrent deletions of a pelvic enhancer of the Pitx1 gene. Here, we identify molecular features contributing to these recurrent deletions. Pitx1 enhancer sequences form alternative DNA structures in vitro and increase double-strand breaks and deletions in vivo. Enhancer mutability depends on DNA replication direction and is caused by TG-dinucleotide repeats. Modeling shows that elevated mutation rates can influence evolution under demographic conditions relevant for sticklebacks and humans. DNA fragility may thus help explain why the same loci are often used repeatedly during parallel adaptive evolution., (Copyright © 2019 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

Background: Stickleback fish are widely used to study the genetic and ecological basis of phenotypic evolution. Although several major loci have now been identified that contribute to evolutionary differences between wild populations, further study of the phenotypes associated with particular genes and mutations has been limited by the difficulty of generating targeted mutations at precise locations in the stickleback genome., Approach and Aims: We compared different methods of expressing single-guide RNAs (sgRNAs) and Cas9 activity in fertilized stickleback eggs. We used an easily scored pigmentation gene ( SLC24A5 ) to screen for molecular lesions, phenotypic effects, and possible germline transmission of newly induced alleles. We then used the optimized CRISPR methods to target two major evolutionary loci in sticklebacks, KITLG and EDA . We hypothesized that coding region mutations in the KITLG gene would alter body pigmentation and possibly sex determination, and that mutations in the EDA gene would disrupt the formation of most armor plates, fin rays, spines, teeth, and gill rakers., Results: Targeted deletions were successfully induced at each target locus by co-injecting one-cell stage stickleback embryos with either Cas9 mRNA or Cas9 protein, together with sgRNAs designed to protein-coding exons. Founder animals were typically mosaic for multiple mutations, which they transmitted through the germline at overall rates of 21 to 100%. We found that the copy of KITLG on the X chromosome ( KITLGX ) has diverged from the KITLG on the Y chromosome ( KITLGY ). Predicted loss-of-function mutations in the KITLGX gene dramatically altered pigmentation in both external skin and internal organ, but the same was not true for KITLGY mutations. Predicted loss-of-function mutations in either the KITLGX or KITLGY genes did not lead to sex reversal or prevent fertility. Homozygous loss-of-function mutations in the EDA gene led to complete loss of armor plates, severe reduction or loss of most soft rays in the dorsal, anal, and caudal fins, and severe reductions in tooth and gill raker number. In contrast, long dorsal and pelvic spines remained intact in EDA mutant animals, suggesting that common co-segregation of plate loss and spine reduction in wild populations is unlikely to be due to pleiotropic effects of EDA mutations., Conclusion: CRISPR-Cas9 approaches can be used to induce germline mutations in key evolutionary loci in sticklebacks. Targeted coding region mutations confirm an important role for KITLG and EDA in skin pigmentation and armor plate reduction, respectively. They also provide new information about the functions of these genes in other body structures.

Vertebrate pelvic reduction is a classic example of repeated evolution. Recurrent loss of pelvic appendages in sticklebacks has previously been linked to natural mutations in a pelvic enhancer that maps upstream of Pitx1 . The sequence of this upstream PelA enhancer is not conserved to mammals, so we have surveyed a large region surrounding the mouse Pitx1 gene for other possible hind limb control sequences. Here we identify a new pelvic enhancer, PelB , that maps downstream rather than upstream of Pitx1. PelB drives expression in the posterior portion of the developing hind limb, and deleting the sequence from mice alters the size of several hind limb structures. PelB sequences are broadly conserved from fish to mammals. A wild stickleback population lacking the pelvis has an insertion/deletion mutation that disrupts the structure and function of PelB , suggesting that changes in this ancient enhancer contribute to evolutionary modification of pelvic appendages in nature., Competing Interests: AT, TC, CG, YC, CI, DM, DK No competing interests declared, (© 2018, Thompson et al.)

Bipolar disorder (BD) and schizophrenia (SCZ) are highly heritable diseases that affect more than 3% of individuals worldwide. Genome-wide association studies have strongly and repeatedly linked risk for both of these neuropsychiatric diseases to a 100 kb interval in the third intron of the human calcium channel gene CACNA1C. However, the causative mutation is not yet known. We have identified a human-specific tandem repeat in this region that is composed of 30 bp units, often repeated hundreds of times. This large tandem repeat is unstable using standard polymerase chain reaction and bacterial cloning techniques, which may have resulted in its incorrect size in the human reference genome. The large 30-mer repeat region is polymorphic in both size and sequence in human populations. Particular sequence variants of the 30-mer are associated with risk status at several flanking single-nucleotide polymorphisms in the third intron of CACNA1C that have previously been linked to BD and SCZ. The tandem repeat arrays function as enhancers that increase reporter gene expression in a human neural progenitor cell line. Different human arrays vary in the magnitude of enhancer activity, and the 30-mer arrays associated with increased psychiatric disease risk status have decreased enhancer activity. Changes in the structure and sequence of these arrays likely contribute to changes in CACNA1C function during human evolution and may modulate neuropsychiatric disease risk in modern human populations., (Copyright © 2018. Published by Elsevier Inc.)

A substantial part of biodiversity is thought to have arisen from adaptive radiations in which one lineage rapidly diversified into multiple lineages specialized to many different niches. However, selection and drift reduce genetic variation during adaptation to new niches and may thus prevent or slow down further niche shifts. We tested whether rapid adaptation is still possible from a highly derived ecotype in the adaptive radiation of threespine stickleback on the Haida Gwaii archipelago, Western Canada. In a 19-year selection experiment, we let giant sticklebacks from a large blackwater lake evolve in a small clearwater pond without vertebrate predators. A total of 56 whole genomes from the experiment and 26 natural populations revealed that adaptive genomic change was rapid in many small genomic regions and encompassed 75% of the change between 12,000-year-old ecotypes. Genomic change was as fast as phenotypic change in defence and trophic morphology, and both were largely parallel between the short-term selection experiment and long-term natural adaptive radiation. Our results show that functionally relevant standing genetic variation can persist in derived radiation members, allowing adaptive radiations to unfold very rapidly.

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., (© 2018 Lowe et al.; Published by Cold Spring Harbor Laboratory Press.)

Background: Dorsal spine reduction in threespine sticklebacks (Gasterosteus aculeatus) is a classic example of recurrent skeletal evolution in nature. Sticklebacks in marine environments typically have long spines that form part of their skeletal armor. Many derived freshwater populations have evolved shorter spines. Changes in spine length are controlled in part by a quantitative trait locus (QTL) previously mapped to chromosome 4, but the causative gene and mutations underlying the repeated evolution of this interesting skeletal trait have not been identified., Results: Refined mapping of the spine length QTL shows that it lies near the MSX2A transcription factor gene. MSX2A is expressed in developing spines. In F1 marine × freshwater fish, the marine allele is preferentially expressed. Differences in expression can be attributed to splicing regulation. Due to the use of an alternative 5 ' splice site within the first exon, the freshwater allele produces greater amounts of a shortened, non-functional transcript and makes less of the full-length transcript. Sequence changes in the MSX2A region are shared by many freshwater fish, suggesting that repeated evolution occurs by reuse of a spine-reduction variant. To demonstrate the effect of full-length MSX2A on spine length, we produced transgenic freshwater fish expressing a copy of marine MSX2A. The spines of the transgenic fish were significantly longer on average than those of their non-transgenic siblings, partially reversing the reduced spine lengths that have evolved in freshwater populations., Conclusions: MSX2A is a major gene underlying dorsal spine reduction in freshwater sticklebacks. The gene is linked to a separate gene controlling bony plate loss, helping explain the concerted effects of chromosome 4 on multiple armor-reduction traits. The nature of the molecular changes provides an interesting example of morphological evolution occurring not through a simple amino acid change, nor through a change only in gene expression levels, but through a change in the ratio of splice products encoding both normal and truncated proteins.

Approximately 15 genes have been directly associated with skin pigmentation variation in humans, leading to its characterization as a relatively simple trait. However, by assembling a global survey of quantitative skin pigmentation phenotypes, we demonstrate that pigmentation is more complex than previously assumed, with genetic architecture varying by latitude. We investigate polygenicity in the KhoeSan populations indigenous to southern Africa who have considerably lighter skin than equatorial Africans. We demonstrate that skin pigmentation is highly heritable, but known pigmentation loci explain only a small fraction of the variance. Rather, baseline skin pigmentation is a complex, polygenic trait in the KhoeSan. Despite this, we identify canonical and non-canonical skin pigmentation loci, including near SLC24A5, TYRP1, SMARCA2/VLDLR, and SNX13, using a genome-wide association approach complemented by targeted resequencing. By considering diverse, under-studied African populations, we show how the architecture of skin pigmentation can vary across humans subject to different local evolutionary pressures., (Copyright © 2017 Elsevier Inc. All rights reserved.)

Ecological speciation with gene flow is widespread in nature [1], but it presents a conundrum: how are associations between traits under divergent natural selection and traits that contribute to assortative mating maintained? Theoretical models suggest that genetic mechanisms inhibiting free recombination between loci underlying these two types of traits (hereafter, "genetic coupling") can facilitate speciation [2-4]. Here, we perform a direct test for genetic coupling by mapping both divergent traits and female mate choice in a classic model of ecological speciation: sympatric benthic and limnetic threespine stickleback (Gasterosteus aculeatus). By measuring mate choice in F2 hybrid females, we allowed for recombination between loci underlying assortative mating and those under divergent ecological selection. In semi-natural mating arenas in which females had access to both benthic and limnetic males, we found that F2 females mated with males similar to themselves in body size and shape. In addition, we found two quantitative trait loci (QTLs) associated with female mate choice that also predicted female morphology along the benthic-limnetic trait axis. Furthermore, a polygenic genetic model that explains adaptation to contrasting benthic and limnetic feeding niches [5] also predicted F2 female mate choice. Together, these results provide empirical evidence that genetic coupling of assortative mating with traits under divergent ecological selection helps maintain species in the face of gene flow, despite a polygenic basis for adaptation to divergent environments., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

The intestinal epithelium serves critical physiologic functions that are shared among all vertebrates. However, it is unknown how the transcriptional regulatory mechanisms underlying these functions have changed over the course of vertebrate evolution. We generated genome-wide mRNA and accessible chromatin data from adult intestinal epithelial cells (IECs) in zebrafish, stickleback, mouse, and human species to determine if conserved IEC functions are achieved through common transcriptional regulation. We found evidence for substantial common regulation and conservation of gene expression regionally along the length of the intestine from fish to mammals and identified a core set of genes comprising a vertebrate IEC signature. We also identified transcriptional start sites and other putative regulatory regions that are differentially accessible in IECs in all 4 species. Although these sites rarely showed sequence conservation from fish to mammals, surprisingly, they drove highly conserved IEC expression in a zebrafish reporter assay. Common putative transcription factor binding sites (TFBS) found at these sites in multiple species indicate that sequence conservation alone is insufficient to identify much of the functionally conserved IEC regulatory information. Among the rare, highly sequence-conserved, IEC-specific regulatory regions, we discovered an ancient enhancer upstream from her6/HES1 that is active in a distinct population of Notch-positive cells in the intestinal epithelium. Together, these results show how combining accessible chromatin and mRNA datasets with TFBS prediction and in vivo reporter assays can reveal tissue-specific regulatory information conserved across 420 million years of vertebrate evolution. We define an IEC transcriptional regulatory network that is shared between fish and mammals and establish an experimental platform for studying how evolutionarily distilled regulatory information commonly controls IEC development and physiology.

Variants in GDF5 are associated with human arthritis and decreased height, but the causal mutations are still unknown. We surveyed the Gdf5 locus for regulatory regions in transgenic mice and fine-mapped separate enhancers controlling expression in joints versus growing ends of long bones. A large downstream regulatory region contains a novel growth enhancer (GROW1), which is required for normal Gdf5 expression at ends of developing bones and for normal bone lengths in vivo. Human GROW1 contains a common base-pair change that decreases enhancer activity and colocalizes with peaks of positive selection in humans. The derived allele is rare in Africa but common in Eurasia and is found in Neandertals and Denisovans. Our results suggest that an ancient regulatory variant in GROW1 has been repeatedly selected in northern environments and that past selection on growth phenotypes explains the high frequency of a GDF5 haplotype that also increases arthritis susceptibility in many human populations.

Repeated adaptation to a new environment often leads to convergent phenotypic changes whose underlying genetic mechanisms are rarely known. Here, we study adaptation of color vision in threespine stickleback during the repeated postglacial colonization of clearwater and blackwater lakes in the Haida Gwaii archipelago. We use whole genomes from 16 clearwater and 12 blackwater populations, and a selection experiment, in which stickleback were transplanted from a blackwater lake into an uninhabited clearwater pond and resampled after 19 y to test for selection on cone opsin genes. Patterns of haplotype homozygosity, genetic diversity, site frequency spectra, and allele-frequency change support a selective sweep centered on the adjacent blue- and red-light sensitive opsins SWS2 and LWS. The haplotype under selection carries seven amino acid changes in SWS2, including two changes known to cause a red-shift in light absorption, and is favored in blackwater lakes but disfavored in the clearwater habitat of the transplant population. Remarkably, the same red-shifting amino acid changes occurred after the duplication of SWS2 198 million years ago, in the ancestor of most spiny-rayed fish. Two distantly related fish species, bluefin killifish and black bream, express these old paralogs divergently in black- and clearwater habitats, while sticklebacks lost one paralog. Our study thus shows that convergent adaptation to the same environment can involve the same genetic changes on very different evolutionary time scales by reevolving lost mutations and reusing them repeatedly from standing genetic variation.

Synovial joints are crucial for support and locomotion in vertebrates, and are the frequent site of serious skeletal defects and degenerative diseases in humans. Growth and differentiation factor 5 (Gdf5) is one of the earliest markers of joint formation, is required for normal joint development in both mice and humans, and has been genetically linked to risk of common osteoarthritis in Eurasian populations. Here, we systematically survey the mouse Gdf5 gene for regulatory elements controlling expression in synovial joints. We identify separate regions of the locus that control expression in axial tissues, in proximal versus distal joints in the limbs, and in remarkably specific sub-sets of composite joints like the elbow. Predicted transcription factor binding sites within Gdf5 regulatory enhancers are required for expression in particular joints. The multiple enhancers that control Gdf5 expression in different joints are distributed over a hundred kilobases of DNA, including regions both upstream and downstream of Gdf5 coding exons. Functional rescue tests in mice confirm that the large flanking regions are required to restore normal joint formation and patterning. Orthologs of these enhancers are located throughout the large genomic region previously associated with common osteoarthritis risk in humans. The large array of modular enhancers for Gdf5 provide a new foundation for studying the spatial specificity of joint patterning in vertebrates, as well as new candidates for regulatory regions that may also influence osteoarthritis risk in human populations., Competing Interests: The authors have declared that no competing interests exist.

Microcapsules are shelled 3D microenvironments, with a liquid core. These core-shelled structures enable cell-cell contact, cellular proliferation and aggregation within the capsule, and can be utilized for controlled release of encapsulated contents. Traditional microcapsule fabrication methods provide limited control of capsule size, and are unable to control capsule placement. To overcome these limitations, we demonstrate size and spatial control of poly-l-lysine and chitosan microcapsules, using laser direct-write (LDW) printing, and subsequent processing, of alginate microbeads. Additionally, microbeads were used as volume pixels (voxels) to form continuous 3D hydrogel structures, which were processed like capsules, to form custom shelled aqueous-core 3D structures of prescribed geometry; such as strands, rings, and bifurcations. Heterogeneous structures were also created with controlled initial locations of different cell types, to demonstrate the ability to prescribe cell signaling (heterotypic and homotypic) in co-culture conditions. Herein, we demonstrate LDW's ability to fabricate intricate 3D structures, essentially with "printed macroporosity," and to precisely control structural composition by bottom-up fabrication in a bead-by-bead manner. The structural and compositional control afforded by this process enables the creation of a wide range of new constructs, with many potential applications in tissue engineering and regenerative medicine. Biotechnol. Bioeng. 2016;113: 2264-2274. © 2016 Wiley Periodicals, Inc., (© 2016 Wiley Periodicals, Inc.)

Changes in bone size and shape are defining features of many vertebrates. Here we use genetic crosses and comparative genomics to identify specific regulatory DNA alterations controlling skeletal evolution. Armor bone-size differences in sticklebacks map to a major effect locus overlapping BMP family member GDF6. Freshwater fish express more GDF6 due in part to a transposon insertion, and transgenic overexpression of GDF6 phenocopies evolutionary changes in armor-plate size. The human GDF6 locus also has undergone distinctive regulatory evolution, including complete loss of an enhancer that is otherwise highly conserved between chimps and other mammals. Functional tests show that the ancestral enhancer drives expression in hindlimbs but not forelimbs, in locations that have been specifically modified during the human transition to bipedalism. Both gain and loss of regulatory elements can localize BMP changes to specific anatomical locations, providing a flexible regulatory basis for evolving species-specific changes in skeletal form., (Copyright © 2016 Elsevier Inc. All rights reserved.)

How predictable is the genetic basis of phenotypic adaptation? Answering this question begins by estimating the repeatability of adaptation at the genetic level. Here, we provide a comprehensive estimate of the repeatability of the genetic basis of adaptive phenotypic evolution in a natural system. We used quantitative trait locus (QTL) mapping to discover genomic regions controlling a large number of morphological traits that have diverged in parallel between pairs of threespine stickleback (Gasterosteus aculeatus species complex) in Paxton and Priest lakes, British Columbia. We found that nearly half of QTL affected the same traits in the same direction in both species pairs. Another 40% influenced a parallel phenotypic trait in one lake but not the other. The remaining 10% of QTL had phenotypic effects in opposite directions in the two species pairs. Similarity in the proportional contributions of all QTL to parallel trait differences was about 0.4. Surprisingly, QTL reuse was unrelated to phenotypic effect size. Our results indicate that repeated use of the same genomic regions is a pervasive feature of parallel phenotypic adaptation, at least in sticklebacks. Identifying the causes of this pattern would aid prediction of the genetic basis of phenotypic evolution., (Copyright © 2015 by the Genetics Society of America.)

Bone morphogenetic proteins (BMPs) are key signaling molecules required for normal development of bones and other tissues. Previous studies have shown that null mutations in the mouse Bmp5 gene alter the size, shape and number of multiple bone and cartilage structures during development. Bmp5 mutations also delay healing of rib fractures in adult mutants, suggesting that the same signals used to pattern embryonic bone and cartilage are also reused during skeletal regeneration and repair. Despite intense interest in BMPs as agents for stimulating bone formation in clinical applications, little is known about the regulatory elements that control developmental or injury-induced BMP expression. To compare the DNA sequences that activate gene expression during embryonic bone formation and following acute injuries in adult animals, we assayed regions surrounding the Bmp5 gene for their ability to stimulate lacZ reporter gene expression in transgenic mice. Multiple genomic fragments, distributed across the Bmp5 locus, collectively coordinate expression in discrete anatomic domains during normal development, including in embryonic ribs. In contrast, a distinct regulatory region activated expression following rib fracture in adult animals. The same injury control region triggered gene expression in mesenchymal cells following tibia fracture, in migrating keratinocytes following dorsal skin wounding, and in regenerating epithelial cells following lung injury. The Bmp5 gene thus contains an "injury response" control region that is distinct from embryonic enhancers, and that is activated by multiple types of injury in adult animals., (Copyright © 2015 Elsevier Inc. All rights reserved.)

Armor plate changes in sticklebacks are a classic example of repeated adaptive evolution. Previous studies identified ectodysplasin (EDA) gene as the major locus controlling recurrent plate loss in freshwater fish, though the causative DNA alterations were not known. Here we show that freshwater EDA alleles have cis-acting regulatory changes that reduce expression in developing plates and spines. An identical T → G base pair change is found in EDA enhancers of divergent low-plated fish. Recreation of the T → G change in a marine enhancer strongly reduces expression in posterior armor plates. Bead implantation and cell culture experiments show that Wnt signaling strongly activates the marine EDA enhancer, and the freshwater T → G change reduces Wnt responsiveness. Thus parallel evolution of low-plated sticklebacks has occurred through a shared DNA regulatory change, which reduces the sensitivity of an EDA enhancer to Wnt signaling, and alters expression in developing armor plates while preserving expression in other tissues.

Most vertebrate organs are composed of epithelium surrounded by support and stromal tissues formed from mesenchyme cells, which are not generally thought to form organized progenitor pools. Here, we use clonal cell labeling with multicolor reporters to characterize individual mesenchymal progenitors in the developing mouse lung. We observe a diversity of mesenchymal progenitor populations with different locations, movements, and lineage boundaries. Airway smooth muscle (ASM) progenitors map exclusively to mesenchyme ahead of budding airways. Progenitors recruited from these tip pools differentiate into ASM around airway stalks; flanking stalk mesenchyme can be induced to form an ASM niche by a lateral bud or by an airway tip plus focal Wnt signal. Thus, mesenchymal progenitors can be organized into localized and carefully controlled domains that rival epithelial progenitor niches in regulatory sophistication., (Copyright © 2014, American Association for the Advancement of Science.)

Phosphate concentration is tightly regulated at the cellular and organismal levels. The first metazoan phosphate exporter, XPR1, was recently identified, but its in vivo function remains unknown. In a genetic screen, we identified a mutation in a zebrafish ortholog of human XPR1, xpr1b. xpr1b mutants lack microglia, the specialized macrophages that reside in the brain, and also displayed an osteopetrotic phenotype characteristic of defects in osteoclast function. Transgenic expression studies indicated that xpr1b acts autonomously in developing macrophages. xpr1b mutants display no gross developmental defects that may arise from phosphate imbalance. We constructed a targeted mutation of xpr1a, a duplicate of xpr1b in the zebrafish genome, to determine whether Xpr1a and Xpr1b have redundant functions. Single mutants for xpr1a were viable, and double mutants for xpr1b;xpr1a were similar to xpr1b single mutants. Our genetic analysis reveals a specific role for the phosphate exporter Xpr1 in the differentiation of tissue macrophages., (Copyright © 2014 The Authors. Published by Elsevier Inc. All rights reserved.)

Developmental genetic studies of evolved differences in morphology have led to the hypothesis that cis-regulatory changes often underlie morphological evolution. However, because most of these studies focus on evolved loss of traits, the genetic architecture and possible association with cis-regulatory changes of gain traits are less understood. Here we show that a derived benthic freshwater stickleback population has evolved an approximate twofold gain in ventral pharyngeal tooth number compared with their ancestral marine counterparts. Comparing laboratory-reared developmental time courses of a low-toothed marine population and this high-toothed benthic population reveals that increases in tooth number and tooth plate area and decreases in tooth spacing arise at late juvenile stages. Genome-wide linkage mapping identifies largely separate sets of quantitative trait loci affecting different aspects of dental patterning. One large-effect quantitative trait locus controlling tooth number fine-maps to a genomic region containing an excellent candidate gene, Bone morphogenetic protein 6 (Bmp6). Stickleback Bmp6 is expressed in developing teeth, and no coding changes are found between the high- and low-toothed populations. However, quantitative allele-specific expression assays of Bmp6 in developing teeth in F1 hybrids show that cis-regulatory changes have elevated the relative expression level of the freshwater benthic Bmp6 allele at late, but not early, stages of stickleback development. Collectively, our data support a model where a late-acting cis-regulatory up-regulation of Bmp6 expression underlies a significant increase in tooth number in derived benthic sticklebacks.

Ecological differences often evolve early in speciation as divergent natural selection drives adaptation to distinct ecological niches, leading ultimately to reproductive isolation. Although this process is a major generator of biodiversity, its genetic basis is still poorly understood. Here we investigate the genetic architecture of niche differentiation in a sympatric species pair of threespine stickleback fish by mapping the environment-dependent effects of phenotypic traits on hybrid feeding and performance under semi-natural conditions. We show that multiple, unlinked loci act largely additively to determine position along the major niche axis separating these recently diverged species. We also find that functional mismatch between phenotypic traits reduces the growth of some stickleback hybrids beyond that expected from an intermediate phenotype, suggesting a role for epistasis between the underlying genes. This functional mismatch might lead to hybrid incompatibilities that are analogous to those underlying intrinsic reproductive isolation but depend on the ecological context.

Hair color differences are among the most obvious examples of phenotypic variation in humans. Although genome-wide association studies (GWAS) have implicated multiple loci in human pigment variation, the causative base-pair changes are still largely unknown. Here we dissect a regulatory region of the KITLG gene (encoding KIT ligand) that is significantly associated with common blond hair color in northern Europeans. Functional tests demonstrate that the region contains a regulatory enhancer that drives expression in developing hair follicles. This enhancer contains a common SNP (rs12821256) that alters a binding site for the lymphoid enhancer-binding factor 1 (LEF1) transcription factor, reducing LEF1 responsiveness and enhancer activity in cultured human keratinocytes. Mice carrying ancestral or derived variants of the human KITLG enhancer exhibit significant differences in hair pigmentation, confirming that altered regulation of an essential growth factor contributes to the classic blond hair phenotype found in northern Europeans.

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.