Evolution of novel sensory organs in fish with legs

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.

Ancient developmental genes underlie evolutionary novelties in walking fish

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.

Fine-scale contemporary recombination variation and its fitness consequences in adaptively diverging stickleback fish

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.

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.

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.

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.

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.

On-Demand Radial Electrodeposition of Alginate Tubular Structures

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.

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.

A novel enhancer near the Pitx1 gene influences development and evolution of pelvic appendages in vertebrates

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.