Widespread selection and gene flow shape the genomic landscape during a radiation of monkeyflowers.

Speciation genomic studies aim to interpret patterns of genome-wide variation in light of the processes that give rise to new species. However, interpreting the genomic "landscape" of speciation is difficult, because many evolutionary processes can impact levels of variation. Facilitated by the first chromosome-level assembly for the group, we use whole-genome sequencing and simulations to shed light on the processes that have shaped the genomic landscape during a radiation of monkeyflowers. After inferring the phylogenetic relationships among the 9 taxa in this radiation, we show that highly similar diversity (π) and differentiation (FST) landscapes have emerged across the group. Variation in these landscapes was strongly predicted by the local density of functional elements and the recombination rate, suggesting that the landscapes have been shaped by widespread natural selection. Using the varying divergence times between pairs of taxa, we show that the correlations between FST and genome features arose almost immediately after a population split and have become stronger over time. Simulations of genomic landscape evolution suggest that background selection (BGS; i.e., selection against deleterious mutations) alone is too subtle to generate the observed patterns, but scenarios that involve positive selection and genetic incompatibilities are plausible alternative explanations. Finally, tests for introgression among these taxa reveal widespread evidence of heterogeneous selection against gene flow during this radiation. Combined with previous evidence for adaptation in this system, we conclude that the correlation in FST among these taxa informs us about the processes contributing to adaptation and speciation during a rapid radiation.

Evolutionary divergence of a Hoxa2b hindbrain enhancer in syngnathids mimics results of functional assays.

Hoxa2 genes provide critical patterning signals during development, and their regulation and function have been extensively studied. We report a previously uncharacterized significant sequence divergence of a highly conserved hindbrain hoxa2b enhancer element in the family syngnathidae (pipefishes, seahorses, pipehorses, seadragons). We compared the hox cis-regulatory element variation in the Gulf pipefish and two species of seahorse against eight other species of fish, as well as human and mouse. We annotated the hoxa2b enhancer element binding sites across three species of seahorse, four species of pipefish, and one species of ghost pipefish. Finally, we performed in situ hybridization analysis of hoxa2b expression in Gulf pipefish embryos. We found that all syngnathid fish examined share a modified rhombomere 4 hoxa2b enhancer element, despite the fact that this element has been found to be highly conserved across all vertebrates examined previously. Binding element sequence motifs and spacing between binding elements have been modified for the hoxa2b enhancer in several species of pipefish and seahorse, and that the loss of the Prep/Meis binding site and further space shortening happened after ghost pipefish split from the rest of the syngnathid clade. We showed that expression of this gene in rhombomere 4 is lower relative to the surrounding rhombomeres in developing Gulf pipefish embryos, reflecting previously published functional tests for this enhancer. Our findings highlight the benefits of studying highly derived, diverse taxa for understanding of gene regulatory evolution and support the hypothesis that natural mutations can occur in deeply conserved pathways in ways potentially related to phenotypic diversity.

The Caudal Skeleton of the Zebrafish, Danio rerio, from a Phylogenetic Perspective: A Polyural Interpretation of Homologous Structures.

The structure of the caudal skeleton of extant teleost fishes has been interpreted in two different ways. In a diural interpretation, a caudal skeleton is composed of two centra articulated with one to six hypurals. Most subsequent authors have followed this interpretation. In contrast, a polyural interpretation considers the teleost fin to be derived from a fully metameristic ancestral bauplan originally composed of a one-to-one relationship between neural arches, centra (when present), and hypurals. Three different interpretations of the identity and homology of skeletal components of the caudal skeleton of the teleost fish Danio rerio have been proposed, two from a diural perspective and one from a polyural perspective. We examine each caudal skeletal component of Danio rerio from both a developmental and phylogenetic perspective. We propose that a polyural interpretation of structures is consistent with the current interpretation of the basal neopterygian caudal fin for this model organism rather than the older diural interpretation that does not take into account the metamerism observed in caudal structures during development. The polyural interpretation suggests several shared evolutionary innovations of major clades that would remain undiscovered under the older diural naming paradigm and makes the terminology of the parts of the caudal fin of Danio rerio strictly comparable to more basal fishes.

Digits in a dish: An in vitro system to assess the molecular genetics of hand/foot development at single-cell resolution.

In vitro models allow for the study of developmental processes outside of the embryo. To gain access to the cells mediating digit and joint development, we identified a unique property of undifferentiated mesenchyme isolated from the distal early autopod to autonomously re-assemble forming multiple autopod structures including: digits, interdigital tissues, joints, muscles and tendons. Single-cell transcriptomic analysis of these developing structures revealed distinct cell clusters that express canonical markers of distal limb development including: Col2a1, Col10a1, and Sp7 (phalanx formation), Thbs2 and Col1a1 (perichondrium), Gdf5, Wnt5a, and Jun (joint interzone), Aldh1a2 and Msx1 (interdigital tissues), Myod1 (muscle progenitors), Prg4 (articular perichondrium/articular cartilage), and Scx and Tnmd (tenocytes/tendons). Analysis of the gene expression patterns for these signature genes indicates that developmental timing and tissue-specific localization were also recapitulated in a manner similar to the initiation and maturation of the developing murine autopod. Finally, the in vitro digit system also recapitulates congenital malformations associated with genetic mutations as in vitro cultures of Hoxa13 mutant mesenchyme produced defects present in Hoxa13 mutant autopods including digit fusions, reduced phalangeal segment numbers, and poor mesenchymal condensation. These findings demonstrate the robustness of the in vitro digit system to recapitulate digit and joint development. As an in vitro model of murine digit and joint development, this innovative system will provide access to the developing limb tissues facilitating studies to discern how digit and articular joint formation is initiated and how undifferentiated mesenchyme is patterned to establish individual digit morphologies. The in vitro digit system also provides a platform to rapidly evaluate treatments aimed at stimulating the repair or regeneration of mammalian digits impacted by congenital malformation, injury, or disease.

Single-Cell Identification of Melanoma Biomarkers in Circulating Tumor Cells.

The current standard for investigating tumors is surgical biopsy, which is costly, invasive, and difficult to perform serially. As an adjunct, circulating tumor cells (CTCs)-cells that have broken away from the primary tumor or metastatic sites-can be obtained from a blood draw and offer the potential for obtaining serial genetic information and serving as biomarkers. Here, we detail the potential for melanoma CTCs to serve as biomarkers and discuss a clinically viable methodology for single-cell CTC isolation and analysis that overcomes previous limitations. We explore the use of melanoma CTC biomarkers by isolating and performing single-cell RNA sequencing on CTCs from melanoma patients. We then compared transcriptional profiles of single melanoma CTCs against A375 cells and peripheral blood mononuclear cells to identify unique genes differentially regulated in circulating melanoma tumor cells. The information that can be obtained via analysis of these CTCs has significant potential in disease tracking.

Genetic analysis of multiple primary melanomas arising within the boundaries of congenital nevi depigmentosa.

Here, we present a rare case of a patient who developed multiple primary melanomas within the boundaries of two nevi depigmentosa. The melanomas were excised, and as a preventive measure, the remainder of the nevi depigmentosa were removed. We performed whole-exome sequencing on excised tissue from the nevus depigmentosus, adjacent normal skin, and saliva to explain this intriguing phenomenon. We also performed a GeneTrails Comprehensive Solid Tumor Panel analysis on one of the melanoma tissues. Genetic analysis revealed germline MC1R V92M and TYR R402Q polymorphisms and a MET E168D germline mutation that may have increased the risk of melanoma development. This genetic predisposition, combined with a patient-reported history of substantial sun exposure and sunburns, which were more severe within the boundaries of the nevi depigmentosa due to the lack of photoprotective melanin, produced numerous somatic mutations in the melanocytes of the nevi depigmentosa. Fitting with this paradigm for melanoma development in chronically sun-damaged skin, the patient's melanomas harbored somatic mutations in CDKN2A (splice site), NF1, and ATRX and had a tumor mutation burden in the 90-95th percentile for melanoma.