Boundary stacking interactions enable cross-TAD enhancer-promoter communication during limb development.

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.

The genetic basis of novel trait gain in walking fish.

A major goal in biology is to understand how organisms evolve novel traits. Multiple studies have identified genes contributing to regressive evolution, the loss of structures that existed in a recent ancestor. However, fewer examples exist for genes underlying constructive evolution, the gain of novel structures and capabilities in lineages that previously lacked them. Sea robins are fish that have evolved enlarged pectoral fins, six mobile locomotory fin rays (legs) and six novel macroscopic lobes in the central nervous system (CNS) that innervate the corresponding legs. Here, we establish successful husbandry and use a combination of transcriptomics, CRISPR-Cas9 editing, and behavioral assays to identify key transcription factors that are required for leg formation and function in sea robins. We also generate hybrids between two sea robin species with distinct leg morphologies and use allele-specific expression analysis and gene editing to explore the genetic basis of species-specific trait diversity, including a novel sensory gain of function. Collectively, our study establishes sea robins as a new model for studying the genetic basis of novel organ formation, and demonstrates a crucial role for the conserved limb gene tbx3a in the evolution of chemosensory legs in walking fish.

Whole-genome Comparisons Identify Repeated Regulatory Changes Underlying Convergent Appendage Evolution in Diverse Fish Lineages.

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.

Sea robins.

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.

Sophorolipid Candidates Demonstrate Cytotoxic Efficacy against 2D and 3D Breast Cancer Models.

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 IC50 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 IC50 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.

Whole-genome comparisons identify repeated regulatory changes underlying convergent appendage evolution in diverse fish lineages.

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 wild fish species that 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 (CONDELs) 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, s yngnathid, and synbranchid species with caudal fin reduction, and which alters caudal fin development when targeted by genome editing. Our study illustrates a general 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 during 100 million years of fish evolution.

Viscoelastic Properties of Bioprinted Alginate Microbeads Compared to Their Bulk Hydrogel Analogs.

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.

Evolution of stickleback spines through independent cis-regulatory changes at HOXDB.

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.

Genomic changes underlying repeated niche shifts in an adaptive radiation.

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.

Characterization of mouse Bmp5 regulatory injury element in zebrafish wound models.

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.

Genetic studies of human-chimpanzee divergence using stem cell fusions.

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.

Longer or shorter spines: Reciprocal trait evolution in stickleback via triallelic regulatory changes in Stanniocalcin2a.

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.

Predicting future from past: The genomic basis of recurrent and rapid stickleback evolution.

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.

Fitness maps to a large-effect locus in introduced stickleback populations.

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.

Non-Destructive Tumor Aggregate Morphology and Viability Quantification at Cellular Resolution, During Development and in Response to Drug.

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.

Assembly of the threespine stickleback Y chromosome reveals convergent signatures of sex chromosome evolution.

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.

Laser-based 3D bioprinting for spatial and size control of tumor spheroids and embryoid bodies.

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.

Multi-modal characterization of polymeric gels to determine the influence of testing method on observed elastic modulus.

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.

DNA fragility in the parallel evolution of pelvic reduction in stickleback fish.

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.

Efficient CRISPR-Cas9 editing of major evolutionary loci in sticklebacks.

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.