Parallelism and Epistasis in Skeletal Evolution Identified through Use of Phylogenomic Mapping Strategies

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

WGCCRR: a web-based tool for genome-wide screening of convergent indels and substitutions of amino acids

Summary

Genome-wide analyses of proteincoding gene sequences are being employed to examine the genetic basis of adaptive evolution in many organismal groups. Previous studies have revealed that convergent/parallel adaptive evolution may be caused by convergent/parallel amino acid changes. Similarly, detailed analysis of lineage-specific amino acid changes has shown correlations with certain lineage-specific traits. However, experimental validation remains the ultimate measure of causality. With the increasing availability of genomic data, a streamlined tool for such analyses would facilitate and expedite the screening of genetic loci that hold potential for adaptive evolution, while alleviating the bioinformatic burden for experimental biologists. In this study, we present a user-friendly web-based tool called WGCCRR (Whole Genome Comparative Coding Region Read) designed to screen both convergent/parallel and lineage-specific amino acid changes on a genome-wide scale. Our tool allows users to replicate previous analyses with just a few clicks, and the exported results are straightforward to interpret. In addition, we have also included amino acid indels that are usually neglected in previous work. Our website provides an efficient platform for screening candidate loci for downstream experimental tests.

Availability and Implementation

The tool is available at: https://fishevo.xmu.edu.cn/.

Transcriptomic and proteomic strategies to reveal the mechanism of Gymnocypris przewalskii scale development

BACKGROUND

Fish scales are typical products of biomineralization and play an important role in the adaptation of fish to their environment. The Gymnocypris przewalskii scales are highly specialized, with scales embedded in only specific parts of the dermis, such as the areas around the anal fin and branchiostegite, making G. przewalskii an ideal material for biomineralization research. In this study, we aimed to unveil genes and pathways controlling scale formation through an integrated analysis of both transcriptome and proteome, of which G. przewalskii tissues of the dorsal skin (no scales) and the rump side skin (with scales) were sequenced. The sequencing results were further combined with cellular experiments to clarify the relationship between genes and signaling pathways.

RESULTS

The results indicated the following: (1) a total of 4,904 differentially expressed genes were screened out, including 3,294 upregulated genes and 1,610 downregulated genes (with a filtering threshold of |log2Fold-Change|> 1 and p-adjust < 0.05). The identified differentially expressed genes contained family members such as FGF, EDAR, Wnt10, and bmp. (2) A total of 535 differentially expressed proteins (DEPs) were filtered out from the proteome, with 204 DEPs downregulated and 331 DEPs upregulated (with a filtering threshold of |Fold-Change|> 1.5 and p < 0.05). (3) Integrated analyses of transcriptome and proteome revealed that emefp1, col1a1, col6a2, col16a1, krt8, and krt18 were important genes contributing to scale development and that PI3K-AKT was the most important signaling pathway involved. (4) With the use of the constructed G. przewalskii fibroblast cell line, emefp1, col1a1, col6a2, col16a1, krt8, and krt18 were confirmed to be positively regulated by the PI3K-AKT signaling pathway.

CONCLUSION

This study provides experimental evidence for PI3K-AKT controlled scale development in G. przewalskii and would benefit further study on stress adaptation, scale biomineralization, and the development of skin appendages.

An active traveling wave of Eda/NF-κB signaling controls the timing and hexagonal pattern of skin appendages in zebrafish

Periodic patterns drive the formation of a variety of tissues, including skin appendages such as feathers and scales. Skin appendages serve important and diverse functions across vertebrates, yet the mechanisms that regulate their patterning are not fully understood. Here, we have used live imaging to investigate dynamic signals regulating the ontogeny of zebrafish scales. Scales are bony skin appendages that develop sequentially along the anterior-posterior and dorsal-ventral axes to cover the fish in a hexagonal array. We have found that scale development requires cell-cell communication and is coordinated through an active wave mechanism. Using a live transcriptional reporter, we show that a wave of Eda/NF-κB activity precedes scale initiation and is required for scale formation. Experiments decoupling the propagation of the wave from dermal placode formation and osteoblast differentiation demonstrate that the Eda/NF-κB activity wavefront controls the timing of the sequential patterning of scales. Moreover, this decoupling resulted in defects in scale size and significant deviations in the hexagonal patterning of scales. Thus, our results demonstrate that a biochemical traveling wave coordinates scale initiation and proper hexagonal patterning across the fish body.

Transcriptomic profiling of tissue environments critical for post-embryonic patterning and morphogenesis of zebrafish skin

Pigment patterns and skin appendages are prominent features of vertebrate skin. In zebrafish, regularly patterned pigment stripes and an array of calcified scales form simultaneously in the skin during post-embryonic development. Understanding the mechanisms that regulate stripe patterning and scale morphogenesis may lead to the discovery of fundamental mechanisms that govern the development of animal form. To learn about cell types and signaling interactions that govern skin patterning and morphogenesis, we generated and analyzed single-cell transcriptomes of skin from wild-type fish as well as fish having genetic or transgenically induced defects in squamation or pigmentation. These data reveal a previously undescribed population of epidermal cells that express transcripts encoding enamel matrix proteins, suggest hormonal control of epithelial-mesenchymal signaling, clarify the signaling network that governs scale papillae development, and identify a critical role for the hypodermis in supporting pigment cell development. Additionally, these comprehensive single-cell transcriptomic data representing skin phenotypes of biomedical relevance should provide a useful resource for accelerating the discovery of mechanisms that govern skin development and homeostasis.

An active traveling wave of Eda/NF-kB signaling controls the timing and hexagonal pattern of skin appendages in zebrafish

Periodic patterns make up a variety of tissues, including skin appendages such as feathers and scales. Skin appendages serve important and diverse functions across vertebrates, yet the mechanisms that regulate their patterning are not fully understood. Here, we have used live imaging to investigate dynamic signals regulating the ontogeny of zebrafish scales. Scales are bony skin appendages which develop sequentially along the anterior-posterior and dorsal-ventral axes to cover the fish in a hexagonal array. We have found that scale development requires cell-cell communication and is coordinated through an active wave mechanism. Using a live transcriptional reporter, we show that a wave of Eda/NF-κB activity precedes scale initiation and is required for scale formation. Experiments decoupling the propagation of the wave from dermal placode formation and osteoblast differentiation demonstrate that the Eda/NF-kB activity wavefront times the sequential patterning of scales. Moreover, this decoupling resulted in defects in scale size and significant deviations in the hexagonal patterning of scales. Thus, our results demonstrate that a biochemical traveling wave coordinates scale initiation and proper hexagonal patterning across the fish body.

Dermal appendage-dependent patterning of zebrafish atoh1a+ Merkel cells

Touch system function requires precise interactions between specialized skin cells and somatosensory axons, as exemplified by the vertebrate mechanosensory Merkel cell-neurite complex. Development and patterning of Merkel cells and associated neurites during skin organogenesis remain poorly understood, partly due to the in utero development of mammalian embryos. Here, we discover Merkel cells in the zebrafish epidermis and identify Atonal homolog 1a (Atoh1a) as a marker of zebrafish Merkel cells. We show that zebrafish Merkel cells derive from basal keratinocytes, express neurosecretory and mechanosensory machinery, extend actin-rich microvilli, and complex with somatosensory axons, all hallmarks of mammalian Merkel cells. Merkel cells populate all major adult skin compartments, with region-specific densities and distribution patterns. In vivo photoconversion reveals that Merkel cells undergo steady loss and replenishment during skin homeostasis. Merkel cells develop concomitant with dermal appendages along the trunk and loss of Ectodysplasin signaling, which prevents dermal appendage formation, reduces Merkel cell density by affecting cell differentiation. By contrast, altering dermal appendage morphology changes the distribution, but not density, of Merkel cells. Overall, our studies provide insights into touch system maturation during skin organogenesis and establish zebrafish as an experimentally accessible in vivo model for the study of Merkel cell biology.

Convergent genomics of longevity in rockfishes highlights the genetics of human life span variation

Longevity is a defining, heritable trait that varies dramatically between species. To resolve the genetic regulation of this trait, we have mined genomic variation in rockfishes, which range in longevity from 11 to over 205 years. Multiple shifts in rockfish longevity have occurred independently and in a short evolutionary time frame, thus empowering convergence analyses. Our analyses reveal a common network of genes under convergent evolution, encompassing established aging regulators such as insulin signaling, yet also identify flavonoid (aryl-hydrocarbon) metabolism as a pathway modulating longevity. The selective pressures on these pathways indicate the ancestral state of rockfishes was long lived and that the changes in short-lived lineages are adaptive. These pathways were also used to explore genome-wide association studies of human longevity, identifying the aryl-hydrocarbon metabolism pathway to be significantly associated with human survival to the 99th percentile. This evolutionary intersection defines and cross-validates a previously unappreciated genetic architecture that associates with the evolution of longevity across vertebrates.

Speciation genomics and the role of depth in the divergence of rockfishes ( Sebastes ) revealed through Pool‐seq analysis of enriched sequences

Speciation in the marine environment is challenged by the wide geographic distribution of many taxa and potential for high rates of gene flow through larval dispersal mechanisms. Depth has recently been proposed as a potential driver of ecological divergence in fishes, and yet it is unclear how adaptation along these gradients' shapes genomic divergence. The genus Sebastes contains numerous species pairs that are depth‐segregated and can provide a better understanding of the mode and tempo of genomic diversification. Here, we present exome data on two species pairs of rockfishes that are depth‐segregated and have different degrees of divergence: S. chlorostictus–S. rosenblatti and S. crocotulus–S. miniatus. We were able to reliably identify “islands of divergence” in the species pair with more recent divergence (S. chlorostictus–S. rosenblatti) and discovered a number of genes associated with neurosensory function, suggesting a role for this pathway in the early speciation process. We also reconstructed demographic histories of divergence and found the best supported model was isolation followed by asymmetric secondary contact for both species pairs. These results suggest past ecological/geographic isolation followed by asymmetric secondary contact of deep to shallow species. Our results provide another example of using rockfish as a model for studying speciation and support the role of depth as an important mechanism for diversification in the marine environment.

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

Seadragons are widely recognized for their derived traits, which include leaf-like appendages and extreme spinal curvature. Efforts to understand the genetic basis of these unique traits and conserve these species and their relatives have been limited by genomic resource gaps. In this paper we present full, annotated genomes of leafy and weedy seadragons, which we use to uncover surprising features of gene family and genome architecture evolution that likely relate to the extravagant phenotypic traits of seadragons and their pipefish and seahorse relatives. These genomes and their analyses are important advances for the study of elaborate vertebrate traits, leveraging this diverse, morphologically exceptional group of fishes.

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

Scale Development-Related Genes Identified by Transcriptome Analysis

Scales, as key structures of fish skin, play an important role in physiological function. The study of fish scale development mechanisms provides a basis for exploring the molecular-level developmental differences between scaled and non-scaled fishes. In this study, alizarin red staining was used to divide the different stages of zebrafish (Danio rerio) scale development. Four developmental stages, namely stage I (~17 dpf, scales have not started to grow), stage II (~33 dpf, the point at which scales start to grow), stage III (~41 dpf, the period in which the scales almost cover the whole body), and stage IV (~3 mpf, scales cover the whole body), were determined and used for subsequent transcriptome analysis. WGCNA (weighted correlation network analysis) and DEG (differentially expressed gene) analysis were used for screening the key genes. Based on the comparison between stage II and stage I, 54 hub-genes were identified by WGCNA analysis. Key genes including the Scpp family (Scpp7, Scpp6, Scpp5, and Scpp8), the Fgf family (Fgfr1b and Fgfr3), Tcf7, Wnt10b, Runx2b, and Il2rb were identified by DEG analysis, which indicated that these genes played important roles in the key nodes of scale development signal pathways. Combined with this analysis, the TGF-β, Wnt/β-catenin, and FGF signaling pathways were suggested to be the most important signal pathways for scales starting to grow. This study laid a foundation for exploring the scale development mechanism of other fishes. The scale development candidate genes identified in the current study will facilitate functional gene identifications in the future.

Expression variations in ectodysplasin-A gene (eda) may contribute to morphological divergence of scales in haplochromine cichlids

BACKGROUND

Elasmoid scales are one of the most common dermal appendages and can be found in almost all species of bony fish differing greatly in their shape. Whilst the genetic underpinnings behind elasmoid scale development have been investigated, not much is known about the mechanisms involved in moulding of scales. To investigate the links between gene expression differences and morphological divergence, we inferred shape variation of scales from two different areas of the body (anterior and posterior) stemming from ten haplochromine cichlid species from different origins (Lake Tanganyika, Lake Malawi, Lake Victoria and riverine). Additionally, we investigated transcriptional differences of a set of genes known to be involved in scale development and morphogenesis in fish.

RESULTS

We found that scales from the anterior and posterior part of the body strongly differ in their overall shape, and a separate look on scales from each body part revealed similar trajectories of shape differences considering the lake origin of single investigated species. Above all, nine as well as 11 out of 16 target genes showed expression differences between the lakes for the anterior and posterior dataset, respectively. Whereas in posterior scales four genes (dlx5, eda, rankl and shh) revealed significant correlations between expression and morphological differentiation, in anterior scales only one gene (eda) showed such a correlation. Furthermore, eda displayed the most significant expression difference between species of Lake Tanganyika and species of the other two younger lakes. Finally, we found genetic differences in downstream regions of eda gene (e.g., in the eda-tnfsf13b inter-genic region) that are associated with observed expression differences. This is reminiscent of a genetic difference in the eda-tnfsf13b inter-genic region which leads to gain or loss of armour plates in stickleback.

CONCLUSION

These findings provide evidence for cross-species transcriptional differences of an important morphogenetic factor, eda, which is involved in formation of ectodermal appendages. These expression differences appeared to be associated with morphological differences observed in the scales of haplochromine cichlids indicating potential role of eda mediated signal in divergent scale morphogenesis in fish.

Advancing human disease research with fish evolutionary mutant models

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

Evolution and development: From the pet shop to the pelagic zone

Flying fish and some of their relatives have evolved the ability to elegantly escape predators by gliding through air. A new study - involving a pet shop zebrafish mutant - offers glimpses into how fins might have been modified to enable this stunt.

Modulation of bioelectric cues in the evolution of flying fishes

Changes to allometry, or the relative proportions of organs and tissues within organisms, is a common means for adaptive character change in evolution. However, little is understood about how relative size is specified during development and shaped during evolution. Here, through a phylogenomic analysis of genome-wide variation in 35 species of flying fishes and relatives, we identify genetic signatures in both coding and regulatory regions underlying the convergent evolution of increased paired fin size and aerial gliding behaviors. To refine our analysis, we intersected convergent phylogenomic signatures with mutants with altered fin size identified in distantly related zebrafish. Through these paired approaches, we identify a surprising role for an L-type amino acid transporter, lat4a, and the potassium channel, kcnh2a, in the regulation of fin proportion. We show that interaction between these genetic loci in zebrafish closely phenocopies the observed fin proportions of flying fishes. The congruence of experimental and phylogenomic findings point to conserved, non-canonical signaling integrating bioelectric cues and amino acid transport in the establishment of relative size in development and evolution.

Footprints in the Sand: Deep Taxonomic Comparisons in Vertebrate Genomics to Unveil the Genetic Programs of Human Longevity

With the modern quality, quantity, and availability of genomic sequencing across species, as well as across the expanse of human populations, we can screen for shared signatures underlying longevity and lifespan. Knowledge of these mechanisms would be medically invaluable in combating aging and age-related diseases. The diversity of longevities across vertebrates is an opportunity to look for patterns of genetic variation that may signal how this life history property is regulated, and ultimately how it can be modulated. Variation in human longevity provides a unique window to look for cases of extreme lifespan within a population, as well as associations across populations for factors that influence capacity to live longer. Current large cohort studies support the use of population level analyses to identify key factors associating with human lifespan. These studies are powerful in concept, but have demonstrated limited ability to resolve signals from background variation. In parallel, the expanding catalog of sequencing and annotation from diverse species, some of which have evolved longevities well past a human lifespan, provides independent cases to look at the genomic signatures of longevity. Recent comparative genomic work has shown promise in finding shared mechanisms associating with longevity among distantly related vertebrate groups. Given the genetic constraints between vertebrates, we posit that a combination of approaches, of parallel meta-analysis of human longevity along with refined analysis of other vertebrate clades having exceptional longevity, will aid in resolving key regulators of enhanced lifespan that have proven to be elusive when analyzed in isolation.

Differences in Cell Proliferation and Craniofacial Phenotype of Closely Related Species in the Pupfish Genus Cyprinodon

Understanding the genetic basis for phenotypic differences is fundamental to the study of macroevolutionary patterns of biological diversity. While technological advances in DNA sequencing have made researching genetic variation in wild taxa routine, fully understanding how these variants affect phenotype requires taking the next step to investigate how genetic changes alter cell and tissue interactions that ultimately produce phenotypes. In this article, we investigate a role for cell proliferation as a developmental source of craniofacial diversity in a radiation of 3 species of Cyprinodon from San Salvador Island, Bahamas. Patterns of cell proliferation in the heads of hatching-age fish differ among species of Cyprinodon, and correlate with differences in allometric growth rate among the jaws of 3 distinct species. Regional patterns of cell proliferation in the head are complex, resulting in an unintuitive result in which lower levels of cell proliferation in the posterior head region are associated with the development of relatively larger jaws in one species. We combine these data with previously published morphological and genomic data to show how studying the mechanisms generating phenotype at the cellular and tissue levels of biological organization can help mechanistically link genomic studies with classic morphological studies.

Fgfr3 mutation disrupts chondrogenesis and bone ossification in zebrafish model mimicking CATSHL syndrome partially via enhanced Wnt/β-catenin signaling

CATSHL syndrome, characterized by camptodactyly, tall stature and hearing loss, is caused by loss-of-function mutations of fibroblast growth factor receptors 3 (FGFR3) gene. Most manifestations of patients with CATSHL syndrome start to develop in the embryonic stage, such as skeletal overgrowth, craniofacial abnormalities, however, the pathogenesis of these phenotypes especially the early maldevelopment remains incompletely understood. Furthermore, there are no effective therapeutic targets for this skeleton dysplasia. Methods: We generated fgfr3 knockout zebrafish by CRISPR/Cas9 technology to study the developmental mechanisms and therapeutic targets of CATSHL syndrome. Several zebrafish transgenic lines labeling osteoblasts and chondrocytes, and live Alizarin red staining were used to analyze the dynamical skeleton development in fgfr3 mutants. Western blotting, whole mount in situ hybridization, Edu labeling based cell proliferation assay and Wnt/β-catenin signaling antagonist were used to explore the potential mechanisms and therapeutic targets. Results: We found that fgfr3 mutant zebrafish, staring from early development stage, showed craniofacial bone malformation with microcephaly and delayed closure of cranial sutures, chondroma-like lesion and abnormal development of auditory sensory organs, partially resembling the clinical manifestations of patients with CATSHL syndrome. Further studies showed that fgfr3 regulates the patterning and shaping of pharyngeal arches and the timely ossification of craniofacial skeleton. The abnormal development of pharyngeal arch cartilage is related to the augmented hypertrophy and disordered arrangement of chondrocytes, while decreased proliferation, differentiation and mineralization of osteoblasts may be involved in the delayed maturation of skull bones. Furthermore, we revealed that deficiency of fgfr3 leads to enhanced IHH signaling and up-regulated canonical Wnt/β-catenin signaling, and pharmacological inhibition of Wnt/β-catenin could partially alleviate the phenotypes of fgfr3 mutants. Conclusions: Our study further reveals some novel phenotypes and underlying developmental mechanism of CATSHL syndrome, which deepens our understanding of the pathogenesis of CATSHL and the role of fgfr3 in skeleton development. Our findings provide evidence that modulation of Wnt/β-catenin activity could be a potential therapy for CATSHL syndrome and related skeleton diseases.

The Genetic Basis of Scale-Loss Phenotype in the Rapid Radiation of Takifugu Fishes

Rapid radiation associated with phenotypic divergence and convergence provides an opportunity to study the genetic mechanisms of evolution. Here we investigate the genus Takifugu that has undergone explosive radiation relatively recently and contains a subset of closely-related species with a scale-loss phenotype. By using observations during development and genetic mapping approaches, we show that the scale-loss phenotype of two Takifugu species, T. pardalis Temminck & Schlegel and T. snyderi Abe, is largely controlled by an overlapping genomic segment (QTL). A search for candidate genes underlying the scale-loss phenotype revealed that the QTL region contains no known genes responsible for the evolution of scale-loss phenotype in other fishes. These results suggest that the genes used for the scale-loss phenotypes in the two Takifugu are likely the same, but the genes used for the similar phenotype in Takifugu and distantly related fishes are not the same. Meanwhile, Fgfrl1, a gene predicted to function in a pathway known to regulate bone/scale development was identified in the QTL region. Since Fgfr1a1, another memebr of the Fgf signaling pathway, has been implicated in scale loss/scale shape in fish distantly related to Takifugu, our results suggest that the convergence of the scale-loss phenotype may be constrained by signaling modules with conserved roles in scale development.

Dissecting Trait Variation across Species Barriers

Dissecting the basis of naturally occurring trait variation is one of the central goals of modern genetics. For eukaryotes, classic methods for this purpose rely on screens of recombinants from matings between distinct parents. These tools cannot be used in studies of species that cannot mate to form recombinant progeny in the first place. However, new approaches are coming online to shuffle the genomes of otherwise incompatible species. With them, geneticists can elucidate how evolution built a new trait, even if it happened millions of years ago, in a lineage that is now reproductively cutoff from its closest relatives.

Historical contingency shapes adaptive radiation in Antarctic fishes

Adaptive radiation illustrates links between ecological opportunity, natural selection and the generation of biodiversity. Central to adaptive radiation is the association between a diversifying lineage and the evolution of phenotypic variation that facilitates the use of new environments or resources. However, is not clear whether adaptive evolution or historical contingency is more important for the origin of key phenotypic traits in adaptive radiation. Here we use targeted sequencing of >250,000 loci across 46 species to examine hypotheses concerning the origin and diversification of key traits in the adaptive radiation of Antarctic notothenioid fishes. Contrary to expectations of adaptive evolution, we show that notothenioids experienced a punctuated burst of genomic diversification and evolved key skeletal modifications before the onset of polar conditions in the Southern Ocean. We show that diversifying selection in pathways associated with human skeletal dysplasias facilitates ecologically important variation in buoyancy among Antarctic notothenioid species, and demonstrate the sufficiency of altered trip11, col1a2 and col1a1a function in zebrafish (Danio rerio) to phenocopy skeletal reduction in Antarctic notothenioids. Rather than adaptation being driven by the cooling of the Antarctic, our results highlight the role of historical contingency in shaping the adaptive radiation of notothenioids. Understanding the historical and environmental context for the origin of key traits in adaptive radiations extends beyond reconstructing events that result in evolutionary innovation, as it also provides a context in forecasting the effects of climate change on the stability and evolvability of natural populations.

In Toto Imaging of Dynamic Osteoblast Behaviors in Regenerating Skeletal Bone

Osteoblasts are matrix-depositing cells that can divide and heal bone injuries. Their deep-tissue location and the slow progression of bone regeneration challenge attempts to capture osteoblast behaviors in live tissue at high spatiotemporal resolution. Here, we have developed an imaging platform to monitor and quantify individual and collective behaviors of osteoblasts in adult zebrafish scales, skeletal body armor discs that regenerate rapidly after loss. Using a panel of transgenic lines that visualize and manipulate osteoblasts, we find that a founder pool of osteoblasts emerges through de novo differentiation within one day of scale plucking. These osteoblasts undergo division events that are largely uniform in frequency and orientation to establish a primordium. Osteoblast proliferation dynamics diversify across the primordium by two days after injury, with cell divisions focused near, and with orientations parallel to, the scale periphery, occurring coincident with dynamic localization of fgf20a gene expression. In posterior scale regions, cell elongation events initiate in areas soon occupied by mineralized grooves called radii, beginning approximately 2 days post injury, with patterned osteoblast death events accompanying maturation of these radii. By imaging at single-cell resolution, we detail acquisition of spatiotemporally distinct cell division, motility, and death dynamics within a founder osteoblast pool as bone regenerates.

Wnt/β-catenin regulates an ancient signaling network during zebrafish scale development

Understanding how patterning influences cell behaviors to generate three dimensional morphologies is a central goal of developmental biology. Additionally, comparing these regulatory mechanisms among morphologically diverse tissues allows for rigorous testing of evolutionary hypotheses. Zebrafish skin is endowed with a coat of precisely patterned bony scales. We use in-toto live imaging during scale development and manipulations of cell signaling activity to elucidate core features of scale patterning and morphogenesis. These analyses show that scale development requires the concerted activity of Wnt/β-catenin, Ectodysplasin (Eda) and Fibroblast growth factor (Fgf) signaling. This regulatory module coordinates Hedgehog (HH) dependent collective cell migration during epidermal invagination, a cell behavior not previously implicated in skin appendage morphogenesis. Our analyses demonstrate the utility of zebrafish scale development as a tractable system in which to elucidate mechanisms of developmental patterning and morphogenesis, and suggest a single, ancient origin of skin appendage patterning mechanisms in vertebrates.

Genetic analyses in Lake Malawi cichlids identify new roles for Fgf signaling in scale shape variation

Elasmoid scales are the most common epithelial appendage among vertebrates, however an understanding of the genetic mechanisms that underlie variation in scale shape is lacking. Using an F2 mapping cross between morphologically distinct cichlid species, we identified >40 QTL for scale shape at different body positions. We show that while certain regions of the genome regulate variation in multiple scales, most are specific to scales at distinct positions. This suggests a degree of regional modularity in scale development. We also identified a single QTL for variation in scale shape disparity across the body. Finally, we screened a QTL hotspot for candidate loci, and identified the Fgf receptor fgfr1b as a prime target. Quantitative rtPCR and small molecule manipulation support a role for Fgf signaling in shaping cichlid scales. While Fgfs have previously been implicated in scale loss, these data reveal new roles for the pathway in scale shape variation.

Cyclin-dependent kinase 21 is a novel regulator of proliferation and meiosis in the male germline of zebrafish

Germ cell differentiation and maintenance relies on complex regulation of mitotic and meiotic progression. Cyclin-dependent kinases (CDKs) and their activating cyclin partners are known to have specialized roles in regulating cell cycle progression across tissues, including germ cells. Very little is known about CDK/cyclin function in zebrafish or the regulation of germ cell maintenance and differentiation. In a forward genetic screen for gonadogenesis defects in zebrafish, a mutation disrupting cdk21 (cyclin-dependent kinase 21) was identified, which caused gonad hypoplasia, reduced fertility and failure of female sex specification. The cdk21 gene is unique to fishes, though the encoded protein is related to the D-cyclin partners Cdk4 and Cdk6, which are known G1 cell cycle regulators. In the testis, cdk21 mutant germ cells exhibited cell cycle defects such as diminished proliferation, prolonged meiosis and delayed sperm differentiation. Furthermore, cdk21 mutants failed to maintain germ cells following breeding. Based on these findings, we propose that cdk21 regulates spermatogonial proliferation, progression through meiosis and germline stem cell activation in the testis. In addition, we investigated cdk4 and cdk6 in zebrafish development and found that each has distinct expression patterns in the gonads. Mutant analysis demonstrated that cdk6 was necessary for viability beyond larval stages. In contrast, cdk4 mutants were viable but were all male with low breeding success and sperm overabundance. Our analysis demonstrated that zebrafish harbor three genes of the cdk4/6 family, cdk4, cdk6 and cdk21, with cdk21 having an essential role in germ cell development in the testis.

Looking Beyond the Genes: The Interplay Between Signaling Pathways and Mechanics in the Shaping and Diversification of Epithelial Tissues

The core of Evo-Devo lies in the intuition that the way tissues grow during embryonic development, the way they sustain their structure and function throughout lifetime, and the way they evolve are closely linked. Epithelial tissues are ubiquitous in metazoans, covering the gut and internal branched organs, as well as the skin and its derivatives (ie, teeth). Here, we discuss in vitro, in vivo, and in silico studies on epithelial tissues to illustrate the conserved, dynamical, and complex aspects of their development. We then explore the implications of the dynamical and nonlinear nature of development on the evolution of their size and shape at the phenotypic and genetic levels. In rare cases, when the interplay between signaling and mechanics is well understood at the cell level, it is becoming clear that the structure of development leads to covariation of characters, an integration which in turn provides some predictable structure to evolutionary changes. We suggest that such nonlinear systems are prone to genetic drift, cryptic genetic variation, and context-dependent mutational effects. We argue that experimental and theoretical studies at the cell level are critical to our understanding of the phenotypic and genetic evolution of epithelial tissues, including carcinomas.