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.

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.

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.

A fully-automated method discovers loss of mouse-lethal and human-monogenic disease genes in 58 mammals.

Gene losses provide an insightful route for studying the morphological and physiological adaptations of species, but their discovery is challenging. Existing genome annotation tools focus on annotating intact genes and do not attempt to distinguish nonfunctional genes from genes missing annotation due to sequencing and assembly artifacts. Previous attempts to annotate gene losses have required significant manual curation, which hampers their scalability for the ever-increasing deluge of newly sequenced genomes. Using extreme sequence erosion (amino acid deletions and substitutions) and sister species support as an unambiguous signature of loss, we developed an automated approach for detecting high-confidence gene loss events across a species tree. Our approach relies solely on gene annotation in a single reference genome, raw assemblies for the remaining species to analyze, and the associated phylogenetic tree for all organisms involved. Using human as reference, we discovered over 400 unique human ortholog erosion events across 58 mammals. This includes dozens of clade-specific losses of genes that result in early mouse lethality or are associated with severe human congenital diseases. Our discoveries yield intriguing potential for translational medical genetics and evolutionary biology, and our approach is readily applicable to large-scale genome sequencing efforts across the tree of life.

A functional enrichment test for molecular convergent evolution finds a clear protein-coding signal in echolocating bats and whales.

Distantly related species entering similar biological niches often adapt by evolving similar morphological and physiological characters. How much genomic molecular convergence (particularly of highly constrained coding sequence) contributes to convergent phenotypic evolution, such as echolocation in bats and whales, is a long-standing fundamental question. Like others, we find that convergent amino acid substitutions are not more abundant in echolocating mammals compared to their outgroups. However, we also ask a more informative question about the genomic distribution of convergent substitutions by devising a test to determine which, if any, of more than 4,000 tissue-affecting gene sets is most statistically enriched with convergent substitutions. We find that the gene set most overrepresented (q-value = 2.2e-3) with convergent substitutions in echolocators, affecting 18 genes, regulates development of the cochlear ganglion, a structure with empirically supported relevance to echolocation. Conversely, when comparing to nonecholocating outgroups, no significant gene set enrichment exists. For aquatic and high-altitude mammals, our analysis highlights 15 and 16 genes from the gene sets most affected by molecular convergence which regulate skin and lung physiology, respectively. Importantly, our test requires that the most convergence-enriched set cannot also be enriched for divergent substitutions, such as in the pattern produced by inactivated vision genes in subterranean mammals. Showing a clear role for adaptive protein-coding molecular convergence, we discover nearly 2,600 convergent positions, highlight 77 of them in 3 organs, and provide code to investigate other clades across the tree of life.

An MTF1 binding site disrupted by a homozygous variant in the promoter of ATP7B likely causes Wilson Disease.

Approximately 2% of the human genome accounts for protein-coding genes, yet most known Mendelian disease-causing variants lie in exons or splice sites. Individuals who symptomatically present with monogenic disorders but do not possess function-altering variants in the protein-coding regions of causative genes may harbor variants in the surrounding gene regulatory domains. We present such a case: a male of Afghani descent was clinically diagnosed with Wilson Disease-a disorder of systemic copper buildup-but was found to have no function-altering coding variants in ATP7B (ENST00000242839.4), the typically causative gene. Our analysis revealed the homozygous variant chr13:g.52,586,149T>C (NC_000013.10, hg19) 676 bp into the ATP7B promoter, which disrupts a metal regulatory transcription factor 1 (MTF1) binding site and diminishes expression of ATP7B in response to copper intake, likely resulting in Wilson Disease. Our approach to identify the causative variant can be generalized to systematically discover function-altering non-coding variants underlying disease and motivates evaluation of gene regulatory variants.

Alternative Progenitor Cells Compensate to Rebuild the Coronary Vasculature in Elabela- and Apj-Deficient Hearts.

Organogenesis during embryonic development occurs through the differentiation of progenitor cells. This process is extraordinarily accurate, but the mechanisms ensuring high fidelity are poorly understood. Coronary vessels of the mouse heart derive from at least two progenitor pools, the sinus venosus and endocardium. We find that the ELABELA (ELA)-APJ signaling axis is only required for sinus venosus-derived progenitors. Because they do not depend on ELA-APJ, endocardial progenitors are able to expand and compensate for faulty sinus venosus development in Apj mutants, leading to normal adult heart function. An upregulation of endocardial SOX17 accompanied compensation in Apj mutants, which was also seen in Ccbe1 knockouts, indicating that the endocardium is activated in multiple cases where sinus venosus angiogenesis is stunted. Our data demonstrate that by diversifying their responsivity to growth cues, distinct coronary progenitor pools are able to compensate for each other during coronary development, thereby providing robustness to organ development.

DACH1 stimulates shear stress-guided endothelial cell migration and coronary artery growth through the CXCL12-CXCR4 signaling axis.

Sufficient blood flow to tissues relies on arterial blood vessels, but the mechanisms regulating their development are poorly understood. Many arteries, including coronary arteries of the heart, form through remodeling of an immature vascular plexus in a process triggered and shaped by blood flow. However, little is known about how cues from fluid shear stress are translated into responses that pattern artery development. Here, we show that mice lacking endothelial Dach1 had small coronary arteries, decreased endothelial cell polarization, and reduced expression of the chemokine Cxcl12 Under shear stress in culture, Dach1 overexpression stimulated endothelial cell polarization and migration against flow, which was reversed upon CXCL12/CXCR4 inhibition. In vivo, DACH1 was expressed during early arteriogenesis but was down in mature arteries. Mature artery-type shear stress (high, uniform laminar) specifically down-regulated DACH1, while the remodeling artery-type flow (low, variable) maintained DACH1 expression. Together, our data support a model in which DACH1 stimulates coronary artery growth by activating Cxcl12 expression and endothelial cell migration against blood flow into developing arteries. This activity is suppressed once arteries reach a mature morphology and acquire high, laminar flow that down-regulates DACH1. Thus, we identified a mechanism by which blood flow quality balances artery growth and maturation.

Pericytes are progenitors for coronary artery smooth muscle.

Epicardial cells on the heart's surface give rise to coronary artery smooth muscle cells (caSMCs) located deep in the myocardium. However, the differentiation steps between epicardial cells and caSMCs are unknown as are the final maturation signals at coronary arteries. Here, we use clonal analysis and lineage tracing to show that caSMCs derive from pericytes, mural cells associated with microvessels, and that these cells are present in adults. During development following the onset of blood flow, pericytes at arterial remodeling sites upregulate Notch3 while endothelial cells express Jagged-1. Deletion of Notch3 disrupts caSMC differentiation. Our data support a model wherein epicardial-derived pericytes populate the entire coronary microvasculature, but differentiate into caSMCs at arterial remodeling zones in response to Notch signaling. Our data are the first demonstration that pericytes are progenitors for smooth muscle, and their presence in adult hearts reveals a new potential cell type for targeting during cardiovascular disease.

The sinus venosus contributes to coronary vasculature through VEGFC-stimulated angiogenesis.

Identifying coronary artery progenitors and their developmental pathways could inspire novel regenerative treatments for heart disease. Multiple sources of coronary vessels have been proposed, including the sinus venosus (SV), endocardium and proepicardium, but their relative contributions to the coronary circulation and the molecular mechanisms regulating their development are poorly understood. We created an ApjCreER mouse line as a lineage-tracing tool to map SV-derived vessels onto the heart and compared the resulting lineage pattern with endocardial and proepicardial contributions to the coronary circulation. The data showed a striking compartmentalization to coronary development. ApjCreER-traced vessels contributed to a large number of arteries, capillaries and veins on the dorsal and lateral sides of the heart. By contrast, untraced vessels predominated in the midline of the ventral aspect and ventricular septum, which are vessel populations primarily derived from the endocardium. The proepicardium gave rise to a smaller fraction of vessels spaced relatively uniformly throughout the ventricular walls. Dorsal (SV-derived) and ventral (endocardial-derived) coronary vessels developed in response to different growth signals. The absence of VEGFC, which is expressed in the epicardium, dramatically inhibited dorsal and lateral coronary growth but left vessels on the ventral side unaffected. We propose that complementary SV-derived and endocardial-derived migratory routes unite to form the coronary vasculature and that the former requires VEGFC, revealing its role as a tissue-specific mediator of blood endothelial development.

VEGF-C and aortic cardiomyocytes guide coronary artery stem development.

Coronary arteries (CAs) stem from the aorta at 2 highly stereotyped locations, deviations from which can cause myocardial ischemia and death. CA stems form during embryogenesis when peritruncal blood vessels encircle the cardiac outflow tract and invade the aorta, but the underlying patterning mechanisms are poorly understood. Here, using murine models, we demonstrated that VEGF-C-deficient hearts have severely hypoplastic peritruncal vessels, resulting in delayed and abnormally positioned CA stems. We observed that VEGF-C is widely expressed in the outflow tract, while cardiomyocytes develop specifically within the aorta at stem sites where they surround maturing CAs in both mouse and human hearts. Mice heterozygous for islet 1 (Isl1) exhibited decreased aortic cardiomyocytes and abnormally low CA stems. In hearts with outflow tract rotation defects, misplaced stems were associated with shifted aortic cardiomyocytes, and myocardium induced ectopic connections with the pulmonary artery in culture. These data support a model in which CA stem development first requires VEGF-C to stimulate vessel growth around the outflow tract. Then, aortic cardiomyocytes facilitate interactions between peritruncal vessels and the aorta. Derangement of either step can lead to mispatterned CA stems. Studying this niche for cardiomyocyte development, and its relationship with CAs, has the potential to identify methods for stimulating vascular regrowth as a treatment for cardiovascular disease.