Abnormal expression of the transcriptional regulator and hedgehog (Hh) signaling pathway effector Gli3 is known to trigger congenital disease, most frequently affecting the central nervous system (CNS) and the limbs. Accurate delineation of the genomic cis-regulatory landscape controlling Gli3 transcription during embryonic development is critical for the interpretation of noncoding variants associated with congenital defects. Here, we employed a comparative genomic analysis on fish species with a slow rate of molecular evolution to identify seven previously unknown conserved noncoding elements (CNEs) in Gli3 intronic intervals (CNE15-21). Transgenic assays in zebrafish revealed that most of these elements drive activities in Gli3 expressing tissues, predominantly the fins, CNS, and the heart. Intersection of these CNEs with human disease associated SNPs identified CNE15 as a putative mammalian craniofacial enhancer, with conserved activity in vertebrates and potentially affected by mutation associated with human craniofacial morphology. Finally, comparative functional dissection of an appendage-specific CNE conserved in slowly evolving fish (elephant shark), but not in teleost (CNE14/hs1586) indicates co-option of limb specificity from other tissues prior to the divergence of amniotes and lobe-finned fish. These results uncover a novel subset of intronic Gli3 enhancers that arose in the common ancestor of gnathostomes and whose sequence components were likely gradually modified in other species during the process of evolutionary diversification.

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.

ABSTRACT

Gene deserts are defined as genomic regions devoid of protein coding genes and spanning more than 500 kilobases, collectively encompassing about 25% of the human genome. Approximately 30% of all gene deserts are enriched for conserved elements with cis -regulatory signatures. These are located predominantly near developmental transcription factors (TFs) but despite predicted critical functions, the transcriptional contributions and biological necessity of most gene deserts remain elusive. Here, we explore the cis -regulatory impact of a gene desert flanking the Shox2 gene, a TF indispensable for proximal limb, craniofacial and cardiac pacemaker development. Using a functional genomics approach in mouse embryos we identify the gene desert as a hub for numerous Shox2 -overlapping enhancers arranged in a globular chromatin domain with tissue-specific features. In accordance, using endogenous CRISPR deletion, we demonstrate that the gene desert interval is essential for Shox2 transcriptional control in developing limbs, craniofacial compartments, and the heart. Phenotypically, gene desert ablation leads to pacemaker-related embryonic lethality due to Shox2 depletion in the cardiac sinus venosus. We show that this role is partially mediated through a distal gene desert enhancer, providing evidence for intra-gene desert regulatory robustness. Finally, we uncover a multi-layered functional role of the gene desert by revealing an additional requirement for stylopod morphogenesis, mediated through an array of proximal limb enhancers (PLEs). In summary, our study establishes the Shox2 gene desert as a fundamental genomic unit that controls pleiotropic gene expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.