Phylogenomic analysis of the Lake Kronotskoe species flock of Dolly Varden charr reveals genetic and developmental signatures of sympatric radiation

Recent adaptive radiations provide evolutionary case studies, which provide the context to parse the relationship between genomic variation and the origins of distinct phenotypes. Sympatric radiations of the charr complex (genus Salvelinus) present a trove for phylogenetics as charrs have repeatedly diversified into multiple morphs with distinct feeding specializations. However, species flocks normally comprise only two to three lineages. Dolly Varden charr inhabiting Lake Kronotske represent the most extensive radiation described for the charr genus, containing at least seven lineages, each with defining morphological and ecological traits. Here, we perform the first genome-wide analysis of this species flock to parse the foundations of adaptive change. Our data support distinct, reproductively isolated lineages with little evidence of hybridization. We also find that specific selection on thyroid signaling and craniofacial genes forms a genomic basis for the radiation. Thyroid hormone is further implicated in subsequent lineage partitioning events. These results delineate a clear genetic basis for the diversification of specialized lineages, and highlight the role of developmental mechanisms in shaping the forms generated during adaptive radiation.

Growth patterns of caudal fin rays are informed by both external signals from the regenerating organ and remembered identity autonomous to the local tissue

Regenerating tissues must remember or interpret their spatial position, using this information to restore original size and patterning. The external skeleton of the zebrafish caudal fin is composed of 18 rays; after any portion of the fin is amputated, position-dependent regenerative growth restores each ray to its original length. We tested for transcriptional differences during regeneration of proximal versus distal tissues and identified 489 genes that differed in proximodistal expression. Thyroid hormone directs multiple aspects of ray patterning along the proximodistal axis, and we identified 364 transcripts showing a proximodistal expression pattern that was dependent on thyroid hormone context. To test what aspects of ray positional identity are directed by extrinsic cues versus remembered identity autonomous to the tissue itself, we transplanted distal portions of rays to proximal environments and evaluated regeneration within the new location. While neighboring proximal tissue showed robust expression of scpp7, a transcript with thyroid-regulated proximal enrichment, regenerating rays originating from transplanted distal tissue showed reduced (distal-like) expression during outgrowth. These distal-to-proximal transplants regenerated far beyond the length of the graft itself, indicating that cues from the proximal environment promoted additional growth. Nonetheless, these transplants initially regenerated at a much slower rate compared to controls, suggesting retained memory of distal identity. This early growth retardation caused rays that originated from transplants to become noticeably shorter than their native neighboring rays. While several aspects of fin ray morphology (bifurcation, segment length) were found to be determined by the environment, regeneration speed and ray length are remembered autonomously by tissues, persisting across multiple rounds of amputation and regeneration.

Thyroid hormone regulates proximodistal patterning in fin rays

Processes that regulate size and patterning along an axis must be highly integrated to generate robust shapes; relative changes in these processes underlie both congenital disease and evolutionary change. Fin length mutants in zebrafish have provided considerable insight into the pathways regulating fin size, yet signals underlying patterning have remained less clear. The bony rays of the fins possess distinct patterning along the proximodistal axis, reflected in the location of ray bifurcations and the lengths of ray segments, which show progressive shortening along the axis. Here, we show that thyroid hormone (TH) regulates aspects of proximodistal patterning of the caudal fin rays, regardless of fin size. TH promotes distal gene expression patterns, coordinating ray bifurcations and segment shortening with skeletal outgrowth along the proximodistal axis. This distalizing role for TH is conserved between development and regeneration, in all fins (paired and medial), and between Danio species as well as distantly related medaka. During regenerative outgrowth, TH acutely induces Shh-mediated skeletal bifurcation. Zebrafish have multiple nuclear TH receptors, and we found that unliganded Thrab-but not Thraa or Thrb-inhibits the formation of distal features. Broadly, these results demonstrate that proximodistal morphology is regulated independently from size-instructive signals. Modulating proximodistal patterning relative to size-either through changes to TH metabolism or other hormone-independent pathways-can shift skeletal patterning in ways that recapitulate aspects of fin ray diversity found in nature.

Static Stability and Swim Bladder Volume in the Bluegill Sunfish (Lepomis macrochirus)

Static stability is a property inherent to every organism. More stable bodies benefit from a lower energy cost associated with maintaining a desired orientation, while less stable bodies can be more maneuverable. The static stability of a fish is determined by the relative locations of its center of mass (COM) and center of buoyancy (COB), which may change with changes in swim bladder volume. We hypothesized, however, that fish would benefit from consistent static stability, and predicted that changes in swim bladder volume would not alter the overall pattern of COM and COB locations. We used micro-computed tomography to estimate the locations of the COM and COB in bluegill sunfish (Lepomis macrochirus). Using this technique, we were able to find a small but significant difference between the location of the COM and COB for a given orientation. We found that the swim bladder can change shape within the body cavity, changing relative locations of the COM and COB. At one extreme, the COB is located 0.441 ± 0.007 BL from the snout and 0.190 ± 0.010 BL from the ventral surface of the pelvic girdle, and that the COM is 0.0030 ± 0.0020 BL posterior and 0.0006 ± 0.0005 BL ventral to the COB, a pattern that causes a nose-up pitching torque. At the other extreme, the COM is anterior and dorsal to the COB, a pattern that causes the opposite torque. These changes in location seems to be caused by changes in the shape and centroid location of the swim bladder within the body: The centroid of the swim bladder is located significantly more posteriorly in fish oriented head-down. The air in the bladder "rises" while heavier tissues "sink," driving a change in tissue distribution and changing the location of the COM relative to the COB. Supporting our hypothesis, we found no correlation between swim bladder volume and the distance between the COM and COB. We conclude that bluegill are statically unstable, requiring them to expend energy constantly to maintain their normal orientation, but that the pitch angle of the body could alter the relative locations of COM and COB, changing their static stability.

Dynamics of the Zebrafish Skeleton in Three Dimensions During Juvenile and Adult Development

Zebrafish are a valuable model for normal vertebrate skeletogenesis and the study of myriad bone disorders. Bones grow, ossify and change shape throughout the zebrafish lifetime, and 3D technologies allow us to examine skeletogenic processes in detail through late developmental stages. To facilitate analysis of shape, orientation and tissue density of skeletal elements throughout ontogeny and adulthood, we generated a high-resolution skeletal reference dataset of wild-type zebrafish development. Using microCT technology, we produced 3D models of the skeletons of individuals ranging from 12 to 25 mm standard length (SL). We analyzed the dynamics of skeletal density and volume as they increase during juvenile and adult growth. Our resource allows anatomical comparisons between meristic units within an individual-e.g., we show that the vertebral canal width increases posteriorly along the spine. Further, structures may be compared between individuals at different body sizes: we highlight the shape changes that the lower jaw undergoes as fish mature from juvenile to adult. We show that even reproductively mature adult zebrafish (17-25 mm SL) continue to undergo substantial changes in skeletal morphology and composition with continued adult growth. We provide a segmented model of the adult skull and a series of interactive 3D PDFs at a range of key stages. These resources allow changes in the skeleton to be assessed quantitatively and qualitatively through late stages of development, and can serve as anatomical references for both research and education.

Thyroid hormone shapes craniofacial bones during postembryonic zebrafish development

Changing the shape of craniofacial bones can profoundly alter ecological function, and understanding how developmental conditions sculpt skeletal phenotypes can provide insight into evolutionary adaptations. Thyroid hormone (TH) stimulates metamorphosis and regulates skeletal morphogenesis across vertebrates. To assess the roles of this hormone in sculpting the craniofacial skeleton of a non-metamorphic vertebrate, we tested zebrafish for developmental periods of TH-induced craniofacial shape change. We analyzed shapes of specific bones that function in prey detection, capture and processing. We quantified these elements from late-larval through adult stages under three developmental TH profiles. Under wild-type conditions, each bone progressively grows allometrically into a mature morphology over the course of postembryonic development. In three of the four bones, TH was required to sculpt an adult shape: hypothyroidism inhibited aspects of shape change, and allowed some components of immature shape to be retained into adulthood. Excess developmental TH stimulated aspects of precocious shape change leading to abnormal morphologies in some bones. Skeletal features with functional importance showed high sensitivities to TH, including the transformator process of the tripus, the mandibular symphysis of the lower jaw, the scutiform lamina of the hyomandibula, and the anterior arm of the pharyngeal jaw. In all, we found that TH is necessary for shaping mature morphology of several essential skeletal elements; this requirement is particularly pronounced during larval development. Altered TH titer leads to abnormal morphologies with likely functional consequences, highlighting the potential of TH and downstream pathways as targets for evolutionary change.

Notochord vacuoles absorb compressive bone growth during zebrafish spine formation

The vertebral column or spine assembles around the notochord rod which contains a core made of large vacuolated cells. Each vacuolated cell possesses a single fluid-filled vacuole, and loss or fragmentation of these vacuoles in zebrafish leads to spine kinking. Here, we identified a mutation in the kinase gene dstyk that causes fragmentation of notochord vacuoles and a severe congenital scoliosis-like phenotype in zebrafish. Live imaging revealed that Dstyk regulates fusion of membranes with the vacuole. We find that localized disruption of notochord vacuoles causes vertebral malformation and curving of the spine axis at those sites. Accordingly, in dstyk mutants the spine curves increasingly over time as vertebral bone formation compresses the notochord asymmetrically, causing vertebral malformations and kinking of the axis. Together, our data show that notochord vacuoles function as a hydrostatic scaffold that guides symmetrical growth of vertebrae and spine formation.

Thyroid hormone modulation during zebrafish development recapitulates evolved diversity in danionin jaw protrusion mechanics

Protrusile jaws are a highly useful innovation that has been linked to extensive diversification in fish feeding ecology. Jaw protrusion can enhance the performance of multiple functions, such as suction production and capturing elusive prey. Identifying the developmental factors that alter protrusion ability will improve our understanding of fish diversification. In the zebrafish protrusion arises postmetamorphosis. Fish metamorphosis typically includes significant changes in trophic morphology, accompanies a shift in feeding niche and coincides with increased thyroid hormone production. We tested whether thyroid hormone affects the development of zebrafish feeding mechanics. We found that it affected all developmental stages examined, but that effects were most pronounced after metamorphosis. Thyroid hormone levels affected the development of jaw morphology, feeding mechanics, shape variation, and cranial ossification. Adult zebrafish utilize protrusile jaws, but an absence of thyroid hormone impaired development of the premaxillary bone, which is critical to jaw protrusion. Premaxillae from early juvenile zebrafish and hypothyroid adult zebrafish resemble those from adults in the genera Danionella, Devario, and Microdevario that show little to no jaw protrusion. Our findings suggest that evolutionary changes in how the developing skulls of danionin minnows respond to thyroid hormone may have promoted diversification into different feeding niches.

MicroCT-based phenomics in the zebrafish skeleton reveals virtues of deep phenotyping in a distributed organ system

Phenomics, which ideally involves in-depth phenotyping at the whole-organism scale, may enhance our functional understanding of genetic variation. Here, we demonstrate methods to profile hundreds of phenotypic measures comprised of morphological and densitometric traits at a large number of sites within the axial skeleton of adult zebrafish. We show the potential for vertebral patterns to confer heightened sensitivity, with similar specificity, in discriminating mutant populations compared to analyzing individual vertebrae in isolation. We identify phenotypes associated with human brittle bone disease and thyroid stimulating hormone receptor hyperactivity. Finally, we develop allometric models and show their potential to aid in the discrimination of mutant phenotypes masked by alterations in growth. Our studies demonstrate virtues of deep phenotyping in a spatially distributed organ system. Analyzing phenotypic patterns may increase productivity in genetic screens, and facilitate the study of genetic variants associated with smaller effect sizes, such as those that underlie complex diseases.

Working with zebrafish at postembryonic stages

As the processes of embryogenesis become increasingly well understood, there is growing interest in the development that occurs at later, postembryonic stages. Postembryonic development holds tremendous potential for discoveries of both fundamental and translational importance. Zebrafish, which are small, rapidly and externally developing, and which boast a wealth of genetic resources, are an outstanding model of vertebrate postembryonic development. Nonetheless, there are specific challenges posed by working with zebrafish at these stages, and this chapter is meant to serve as a primer for those working with larval and juvenile zebrafish. Since accurate staging is critical for high-quality results and experimental reproducibility, we outline best practices for reporting postembryonic developmental progress. Emphasizing the importance of accurate staging, we present new data showing that rates of growth and size-stage relationships can differ even between wild-type strains. Finally, since rapid and uniform development is particularly critical when working at postembryonic stages, we briefly describe methods that we use to achieve high rates of growth and developmental uniformity through postembryonic stages in both wild-type and growth-compromised zebrafish.

Thyroid hormone-dependent adult pigment cell lineage and pattern in zebrafish

Pigment patterns are useful for elucidating fundamental mechanisms of pattern formation and how these mechanisms evolve. In zebrafish, several pigment cell classes interact to generate stripes, yet the developmental requirements and origins of these cells remain poorly understood. Using zebrafish and a related species, we identified roles for thyroid hormone (TH) in pigment cell development and patterning, and in postembryonic development more generally. We show that adult pigment cells arise from distinct lineages having distinct requirements for TH and that differential TH dependence can evolve within lineages. Our findings demonstrate critical functions for TH in determining pigment pattern phenotype and highlight the potential for evolutionary diversification at the intersection of developmental and endocrine mechanisms.

Dwarfism and increased adiposity in the gh1 mutant zebrafish vizzini

Somatic growth and adipogenesis are closely associated with the development of obesity in humans. In this study, we identify a zebrafish mutant, vizzini, that exhibits both a severe defect in somatic growth and increased accumulation of adipose tissue. Positional cloning of vizzini revealed a premature stop codon in gh1. Although the effects of GH are largely through igfs in mammals, we found no decrease in the expression of igf transcripts in gh1 mutants during larval development. As development progressed, however, we found overall growth to be progressively retarded and the attainment of specific developmental stages to occur at abnormally small body sizes relative to wild type. Moreover, both subcutaneous (sc) and visceral adipose tissues underwent precocious development in vizzini mutants, and at maturity, the sizes of different fat deposits were greatly expanded relative to wild type. In vivo confocal imaging of sc adipose tissue (SAT) expansion revealed that vizzini mutants exhibit extreme enlargement of adipocyte lipid droplets without a corresponding increase in lipid droplet number. These findings suggest that GH1 signaling restricts SAT hypertrophy in zebrafish. Finally, nutrient deprivation of vizzini mutants revealed that SAT mobilization was greatly diminished during caloric restriction, further implicating GH1 signaling in adipose tissue homeostasis. Overall, the zebrafish gh1 mutant, vizzini, exhibits decreased somatic growth, increased adipose tissue accumulation, and disrupted adipose plasticity after nutrient deprivation and represents a novel model to investigate the in vivo dynamics of vertebrate obesity.

Metamorphosis in teleosts

Teleosts are the largest and most diverse group of vertebrates, and many species undergo morphological, physiological, and behavioral transitions, "metamorphoses," as they progress between morphologically divergent life stages. The larval metamorphosis that generally occurs as teleosts mature from larva to juvenile involves the loss of embryo-specific features, the development of new adult features, major remodeling of different organ systems, and changes in physical proportions and overall phenotype. Yet, in contrast to anuran amphibians, for example, teleost metamorphosis can entail morphological change that is either sudden and profound, or relatively gradual and subtle. Here, we review the definition of metamorphosis in teleosts, the diversity of teleost metamorphic strategies and the transitions they involve, and what is known of their underlying endocrine and genetic bases. We suggest that teleost metamorphosis offers an outstanding opportunity for integrating our understanding of endocrine mechanisms, cellular processes of morphogenesis and differentiation, and the evolution of diverse morphologies and life histories.

Ancient DNA assessment of tiger salamander population in Yellowstone National Park

Recent data indicates that blotched tiger salamanders (Ambystoma tigrinum melanostictum) in northern regions of Yellowstone National Park are declining due to climate-related habitat changes. In this study, we used ancient and modern mitochondrial haplotype diversity to model the effective size of this amphibian population through recent geological time and to assess past responses to climatic changes in the region. Using subfossils collected from a cave in northern Yellowstone, we analyzed >700 base pairs of mitochondrial sequence from 16 samples ranging in age from 100 to 3300 years old and found that all shared an identical haplotype. Although mitochondrial diversity was extremely low within the living population, we still were able to detect geographic subdivision within the local area. Using serial coalescent modelling with Bayesian priors from both modern and ancient genetic data we simulated a range of probable population sizes and mutation rates through time. Our simulations suggest that regional mitochondrial diversity has remained relatively constant even through climatic fluctuations of recent millennia.