A genome-scale analysis of the cis-regulatory circuitry underlying sonic hedgehog-mediated patterning of the mammalian limb

Multiple steps characterise ventricular layer attrition to form the ependymal cell lining of the adult mouse spinal cord central canal

The ventricular layer of the spinal cord is remodelled during embryonic development and ultimately forms the ependymal cell lining of the adult central canal, which retains neural stem cell potential. This anatomical transformation involves the process of dorsal collapse; however, accompanying changes in tissue organisation and cell behaviour as well as the precise origin of cells contributing to the central canal are not well understood. Here, we describe sequential localised cell rearrangements which accompany the gradual attrition of the spinal cord ventricular layer during development. This includes local breakdown of the pseudostratified organisation of the dorsal ventricular layer prefiguring dorsal collapse and evidence for a new phenomenon, ventral dissociation, during which the ventral‐most floor plate cells separate from a subset that are retained around the central canal. Using cell proliferation markers and cell‐cycle reporter mice, we further show that following dorsal collapse, ventricular layer attrition involves an overall reduction in cell proliferation, characterised by an intriguing increase in the percentage of cells in G1/S. In contrast, programmed cell death does not contribute to ventricular layer remodelling. By analysing transcript and protein expression patterns associated with key signalling pathways, we provide evidence for a gradual decline in ventral sonic hedgehog activity and an accompanying ventral expansion of initial dorsal bone morphogenetic protein signalling, which comes to dominate the forming the central canal lining. This study identifies multiple steps that may contribute to spinal cord ventricular layer attrition and adds to increasing evidence for the heterogeneous origin of the spinal cord ependymal cell population, which includes cells from the floor plate and the roof plate as well as ventral progenitor domains.

Coordination of Leaf Development Across Developmental Axes

Leaves are initiated as lateral outgrowths from shoot apical meristems throughout the vegetative life of the plant. To achieve proper developmental patterning, cell-type specification and growth must occur in an organized fashion along the proximodistal (base-to-tip), mediolateral (central-to-edge), and adaxial–abaxial (top-bottom) axes of the developing leaf. Early studies of mutants with defects in patterning along multiple leaf axes suggested that patterning must be coordinated across developmental axes. Decades later, we now recognize that a highly complex and interconnected transcriptional network of patterning genes and hormones underlies leaf development. Here, we review the molecular genetic mechanisms by which leaf development is coordinated across leaf axes. Such coordination likely plays an important role in ensuring the reproducible phenotypic outcomes of leaf morphogenesis.

Dmrt factors determine the positional information of cerebral cortical progenitors via differential suppression of homeobox genes

The spatiotemporal identity of neural progenitors and the regional control of neurogenesis are essential for the development of cerebral cortical architecture. Here, we report that mammalian DM domain factors (Dmrt) determine the identity of cerebral cortical progenitors. Among the Dmrt family genes expressed in the developing dorsal telencephalon, Dmrt3 and Dmrta2 show a medial^high/lateral^low expression gradient. Their simultaneous loss confers a ventral identity to dorsal progenitors, resulting in the ectopic expression of Gsx2 and massive production of GABAergic olfactory bulb interneurons in the dorsal telencephalon. Furthermore, double-mutant progenitors in the medial region exhibit upregulated Pax6 and more lateral characteristics. These ventral and lateral shifts in progenitor identity depend on Dmrt gene dosage. We also found that Dmrt factors bind to Gsx2 and Pax6 enhancers to suppress their expression. Our findings thus reveal that the graded expression of Dmrt factors provide positional information for progenitors by differentially repressing downstream genes in the developing cerebral cortex.

Evolution of the avian digital pattern

Variation in digit number has occurred multiple times in the history of archosaur evolution. The five digits of dinosaur limbs were reduced to three in bird forelimbs, and were further reduced in the vestigial forelimbs of the emu. Regulation of digit number has been investigated previously by examining genes involved in anterior-posterior patterning in forelimb buds among emu (Dromaius novaehollandiae), chicken (Gallus gallus) and zebra finch (Taeniopygia guttata). It was described that the expression of posterior genes are conserved among these three birds, whereas expression of anterior genes Gli3 and Alx4 varied significantly. Here we re-examined the expression pattern of Gli3 and Alx4 in the forelimb of emu, chicken and zebra finch. We found that Gli3 is expressed in the anterior region, although its range varied among species, and that the expression pattern of Alx4 in forelimb buds is broadly conserved in a stage-specific manner. We also found that the dynamic expression pattern of the BMP antagonist Gremlin1 (Grem1) in limb buds, which is critical for autopodial expansion, was consistent with the digital pattern of emu, chicken and zebra finch. Furthermore, in emu, variation among individuals was observed in the width of Grem1 expression in forelimb buds, as well as in the adult skeletal pattern. Our results support the view that the signalling system that regulates the dynamic expression of Grem1 in the limb bud contributes substantially to variations in avian digital patterns.

Understanding Pattern Formation in Embryos: Experiment, Theory, and Simulation

Twenty-five years ago, Lewis Wolpert, the eminent developmental biologist, asked the question, "Do We Understand Development?" He concluded that such rapid progress had been made in the preceding two decades that "It is not unreasonable to think that enough will eventually be known to program a computer and simulate some aspects of development." This prediction has been fulfilled, at least partially, with data-driven simulations of several different developmental processes being developed in the intervening years. Nevertheless, the question remains of whether we "understand" development and if simulations are sufficient to provide an explanation of development. While in silico replications and models are undoubtedly an important tool in the investigation and dissection of developmental processes, which complement traditional experimental methods, these need to be supplemented by theory that identifies principles and provides coherent explanations. Here, I use the example of pattern formation in the vertebrate neural tube to illustrate this idea.

Hoxa13 regulates expression of common Hox target genes involved in cartilage development to coordinate the expansion of the autopodal anlage

To elucidate the role of Hox genes in limb cartilage development, we identified the target genes of HOXA11 and HOXA13 by ChIP‐Seq. The ChIP DNA fragment contained evolutionarily conserved sequences and multiple highly conserved HOX binding sites. A substantial portion of the HOXA11 ChIP fragment overlapped with the HOXA13 ChIP fragment indicating that both factors share common targets. Deletion of the target regions neighboring Bmp2 or Tshz2 reduced their expression in the autopod suggesting that they function as the limb bud‐specific enhancers. We identified the Hox downstream genes as exhibiting expression changes in the Hoxa13 knock out (KO) and Hoxd11‐13 deletion double mutant (Hox13 dKO) autopod by Genechip analysis. The Hox downstream genes neighboring the ChIP fragment were defined as the direct targets of Hox. We analyzed the spatial expression pattern of the Hox target genes that encode two different categories of transcription factors during autopod development and Hox13dKO limb bud. (a) Bcl11a, encoding a repressor of cartilage differentiation, was expressed in the E11.5 autopod and was substantially reduced in the Hox13dKO. (b) The transcription factors Aff3, Bnc2, Nfib and Runx1t1 were expressed in the zeugopodal cartilage but not in the autopod due to the repressive or relatively weak transcriptional activity of Hox13 at E11.5. Interestingly, the expression of these genes was later observed in the autopodal cartilage at E12.5. These results indicate that Hox13 transiently suspends the cartilage differentiation in the autopodal anlage via multiple pathways until establishing the paddle‐shaped structure required to generate five digits.

Hedgehog signaling in the airway epithelium of patients with chronic obstructive pulmonary disease

Genome-wide association studies have linked gene variants of the receptor patched homolog 1 (PTCH1) with chronic obstructive pulmonary disease (COPD). However, its biological role in the disease is unclear. Our objective was to determine the expression pattern and biological role of PTCH1 in the lungs of patients with COPD. Airway epithelial-specific PTCH1 protein expression and epithelial morphology were assessed in lung tissues of control and COPD patients. PTCH1 mRNA expression was measured in bronchial epithelial cells obtained from individuals with and without COPD. The effects of PTCH1 siRNA knockdown on epithelial repair and mucous expression were evaluated using human epithelial cell lines. Ptch1^+/− mice were used to assess the effect of decreased PTCH1 on mucous expression and airway epithelial phenotypes. Airway epithelial-specific PTCH1 protein expression was significantly increased in subjects with COPD compared to controls, and its expression was associated with total airway epithelial cell count and thickness. PTCH1 knockdown attenuated wound closure and mucous expression in airway epithelial cell lines. Ptch1^+/− mice had reduced mucous expression compared to wildtype mice following mucous induction. PTCH1 protein is up-regulated in COPD airway epithelium and may upregulate mucous expression. PTCH1 provides a novel target to reduce chronic bronchitis in COPD patients.

Unique and overlapping GLI1 and GLI2 transcriptional targets in neoplastic chondrocytes

Excessive Hedgehog (Hh) signaling in chondrocytes is sufficient to cause formation of enchondroma-like lesions which can progress to chondrosarcoma. To elucidate potential underlying mechanisms, we identified GLI1 and GLI2 target genes in human chondrosarcoma. Using chromatin immunoprecipitation (ChIP) sequencing and microarray data, in silico analyses were conducted to identify and characterize unique and overlapping GLI1 and GLI2 binding regions in neoplastic chondrocytes. After overlaying microarray data from human chondrosarcoma, 204 upregulated and 106 downregulated genes were identified as Hh-responsive Gli binding targets. After overlaying published Gli ChIP-on-chip data from mouse, 48 genes were identified as potential direct downstream targets of Hedgehog signaling with shared GLI binding regions in evolutionarily conserved DNA elements. Among these was BMP2, pointing to potential cross-talk between TGF beta signaling and Hh signaling. Our identification of potential target genes that are unique and common to GLI1 and GLI2 in neoplastic chondrocytes contributes to elucidating potential pathways through which Hh signaling impacts cartilage tumor biology.

p.E95K mutation in Indian hedgehog causing brachydactyly type A1 impairs IHH/Gli1 downstream transcriptional regulation

Background

Brachydactyly type A1 (BDA1, OMIM 112500) is a rare inherited malformation characterized primarily by shortness or absence of middle bones of fingers and toes. It is the first recorded disorder of the autosomal dominant Mendelian trait. Indian hedgehog (IHH) gene is closely associated with BDA1, which was firstly mapped and identified in Chinese families in 2000. Previous studies have demonstrated that BDA1-related mutant IHH proteins affected interactions with its receptors and impaired IHH signaling. However, how the altered signaling pathway affects downstream transcriptional regulation remains unclear.

Results

Based on the mouse C3H10T1/2 cell model for IHH signaling activation, two recombinant human IHH-N proteins, including a wild type protein (WT, amino acid residues 28–202) and a mutant protein (MT, p.E95k), were analyzed. We identified 347, 47 and 4 Gli1 binding sites in the corresponding WT, MT and control group by chromatin immunoprecipitation and the overlapping of these three sets was poor. The putative cis regulated genes in WT group were enriched in sensory perception and G-protein coupled receptor-signaling pathway. On the other hand, putative cis regulated genes were enriched in Runx2-related pathways in MT group. Differentially expressed genes in WT and MT groups indicated that the alteration of mutant IHH signaling involved cell-cell signaling and cellular migration. Cellular assay of migration and proliferation validated that the mutant IHH signaling impaired these two cellular functions.

Conclusions

In this study, we performed integrated genome-wide analyses to characterize differences of IHH/Gli1 downstream regulation between wild type IHH signaling and the E95K mutant signaling. Based on the cell model, our results demonstrated that the E95K mutant signaling altered Gli1-DNA binding pattern, impaired downstream gene expressions, and leaded to weakened cellular proliferation and migration. This study may help to deepen the understanding of pathogenesis of BDA1 and the role of IHH signaling in chondrogenesis.

Electronic supplementary material

The online version of this article (10.1186/s12863-018-0697-5) contains supplementary material, which is available to authorized users.

Phenotype loss is associated with widespread divergence of the gene regulatory landscape in evolution

Detecting the genomic changes underlying phenotypic changes between species is a main goal of evolutionary biology and genomics. Evolutionary theory predicts that changes in cis-regulatory elements are important for morphological changes. We combined genome sequencing, functional genomics and genome-wide comparative analyses to investigate regulatory elements in lineages that lost morphological traits. We first show that limb loss in snakes is associated with widespread divergence of limb regulatory elements. We next show that eye degeneration in subterranean mammals is associated with widespread divergence of eye regulatory elements. In both cases, sequence divergence results in an extensive loss of transcription factor binding sites. Importantly, diverged regulatory elements are associated with genes required for normal limb patterning or normal eye development and function, suggesting that regulatory divergence contributed to the loss of these phenotypes. Together, our results show that genome-wide decay of the phenotype-specific cis-regulatory landscape is a hallmark of lost morphological traits.

Cis-regulatory elements are important factors for morphological changes. Here, the authors show widespread divergence of limb and eye regulatory elements in limb loss in snakes and eye degeneration in subterranean mammals respectively.

Gli3 controls the onset of cortical neurogenesis by regulating the radial glial cell cycle through Cdk6 expression

The cerebral cortex contains an enormous number of neurons, allowing it to perform highly complex neural tasks. Understanding how these neurons develop at the correct time and place and in accurate numbers constitutes a major challenge. Here, we demonstrate a novel role for Gli3, a key regulator of cortical development, in cortical neurogenesis. We show that the onset of neuron formation is delayed in Gli3 conditional mouse mutants. Gene expression profiling and cell cycle measurements indicate that shortening of the G1 and S phases in radial glial cells precedes this delay. Reduced G1 length correlates with an upregulation of the cyclin-dependent kinase gene Cdk6, which is directly regulated by Gli3. Moreover, pharmacological interference with Cdk6 function rescues the delayed neurogenesis in Gli3 mutant embryos. Overall, our data indicate that Gli3 controls the onset of cortical neurogenesis by determining the levels of Cdk6 expression, thereby regulating neuronal output and cortical size.

Ectopic Hedgehog Signaling Causes Cleft Palate and Defective Osteogenesis

Cleft palate is a common birth defect that frequently occurs in human congenital malformations caused by mutations in components of the Sonic Hedgehog (S HH) signaling cascade. Shh is expressed in dynamic, spatiotemporal domains within epithelial rugae and plays a key role in driving epithelial-mesenchymal interactions that are central to development of the secondary palate. However, the gene regulatory networks downstream of Hedgehog (Hh) signaling are incompletely characterized. Here, we show that ectopic Hh signaling in the palatal mesenchyme disrupts oral-nasal patterning of the neural crest cell-derived ectomesenchyme of the palatal shelves, leading to defective palatine bone formation and fully penetrant cleft palate. We show that a series of Fox transcription factors, including the novel direct target Foxl1, function downstream of Hh signaling in the secondary palate. Furthermore, we demonstrate that Wnt/bone morphogenetic protein (BMP) antagonists, in particular Sostdc1, are positively regulated by Hh signaling, concomitant with downregulation of key regulators of osteogenesis and BMP signaling effectors. Our data demonstrate that ectopic Hh-Smo signaling downregulates Wnt/BMP pathways, at least in part by upregulating Sostdc1, resulting in cleft palate and defective osteogenesis.

Syndromes associated with holoprosencephaly

Holoprosencephaly (HPE) is partial or complete failure of the forebrain to divide into hemispheres and can be an isolated finding or associated with a syndrome. Most cases of HPE are associated with a syndrome and roughly 40%-60% of fetuses with HPE have trisomy 13 which is the most common etiology of HPE. Other syndromes associated with HPE include additional aneuploidies like trisomy 18 and single gene disorders such as Smith-Lemli-Opitz syndrome. There are a number of syndromes such as pseudotrisomy 13 which do not have a known molecular etiology; therefore, this review has two parts: syndromes with a molecular diagnosis and syndromes where the etiology is yet to be found. As most HPE is syndromic, this review provides a comprehensive list and description of syndromes associated with HPE that may be used as a differential diagnosis and starting point for evaluating individuals with HPE.

Synergistic co-regulation and competition by a SOX9-GLI-FOXA phasic transcriptional network coordinate chondrocyte differentiation transitions

The growth plate mediates bone growth where SOX9 and GLI factors control chondrocyte proliferation, differentiation and entry into hypertrophy. FOXA factors regulate hypertrophic chondrocyte maturation. How these factors integrate into a Gene Regulatory Network (GRN) controlling these differentiation transitions is incompletely understood. We adopted a genome-wide whole tissue approach to establish a Growth Plate Differential Gene Expression Library (GP-DGEL) for fractionated proliferating, pre-hypertrophic, early and late hypertrophic chondrocytes, as an overarching resource for discovery of pathways and disease candidates. De novo motif discovery revealed the enrichment of SOX9 and GLI binding sites in the genes preferentially expressed in proliferating and prehypertrophic chondrocytes, suggesting the potential cooperation between SOX9 and GLI proteins. We integrated the analyses of the transcriptome, SOX9, GLI1 and GLI3 ChIP-seq datasets, with functional validation by transactivation assays and mouse mutants. We identified new SOX9 targets and showed SOX9-GLI directly and cooperatively regulate many genes such as Trps1, Sox9, Sox5, Sox6, Col2a1, Ptch1, Gli1 and Gli2. Further, FOXA2 competes with SOX9 for the transactivation of target genes. The data support a model of SOX9-GLI-FOXA phasic GRN in chondrocyte development. Together, SOX9-GLI auto-regulate and cooperate to activate and repress genes in proliferating chondrocytes. Upon hypertrophy, FOXA competes with SOX9, and control toward terminal differentiation passes to FOXA, RUNX, AP1 and MEF2 factors.

In the development of the mammalian growth plate, while several transcription factors are individually well known for their key roles in regulating phases of chondrocyte differentiation, there is little information on how they interact and cooperate with each other. We took an unbiased genome wide approach to identify the transcription factors and signaling pathways that play dominant roles in the chondrocyte differentiation cascade. We developed a searchable library of differentially expressed genes, GP-DGEL, which has fine spatial resolution and global transcriptomic coverage for discovery of processes, pathways and disease candidates. Our work identifies a novel regulatory mechanism that integrates the action of three transcription factors, SOX9, GLI and FOXA. SOX9-GLI auto-regulate and cooperate to activate and repress genes in proliferating chondrocytes. Upon entry into prehypertrophy, FOXA competes with SOX9, and control of hypertrophy passes to FOXA, RUNX, AP1 and MEF2 factors.

A novel role for cilia-dependent sonic hedgehog signaling during submandibular gland development

BACKGROUND

Submandibular glands (SMGs) are specialized epithelial structures which generate saliva necessary for mastication and digestion. Loss of SMGs can lead to inflammation, oral lesions, fungal infections, problems with chewing/swallowing, and tooth decay. Understanding the development of the SMG is important for developing therapeutic options for patients with impaired SMG function. Recent studies have suggested Sonic hedgehog (Shh) signaling in the epithelium plays an integral role in SMG development; however, the mechanism by which Shh influences gland development remains nebulous.

RESULTS

Using the Kif3a^f/f ;Wnt1-Cre ciliopathic mouse model to prevent Shh signal transduction by means of the loss of primary cilia in neural crest cells, we report that mesenchymal Shh activity is necessary for gland development. Furthermore, using a variety of murine transgenic lines with aberrant mesenchymal Shh signal transduction, we determine that loss of Shh activity, by means of loss of the Gli activator, rather than gain of Gli repressor, is sufficient to cause the SMG aplasia. Finally, we determine that loss of the SMG correlates with reduced Neuregulin1 (Nrg1) expression and lack of innervation of the SMG epithelium.

CONCLUSIONS

Together, these data suggest a novel mechanistic role for mesenchymal Shh signaling during SMG development. Developmental Dynamics 247:818-831, 2018. © 2018 Wiley Periodicals, Inc.

Noncanonical hedgehog pathway activation through SRF-MKL1 promotes drug resistance in basal cell carcinomas

Hedgehog pathway-dependent cancers can escape smoothened (SMO) inhibition through canonical pathway mutations, however, 50% of resistant BCCs lack additional variants in hedgehog genes. Here we use multi-dimensional genomics in human and mouse resistant BCCs to identify a non-canonical hedgehog activation pathway driven by the transcription factor, serum response factor (SRF). Active SRF along with its co-activator megakaryoblastic leukemia 1 (MKL1) form a novel protein complex and share chromosomal occupancy with the hedgehog transcription factor GLI1, causing amplification of GLI1 transcriptional activity. We show cytoskeletal activation by Rho and the formin family member Diaphanous (mDia) are required for SRF/MKL-driven GLI1 activation and tumor cell viability. Remarkably, we use nuclear MKL1 staining in mouse and human patient tumors to define drug responsiveness to MKL inhibitors highlighting the therapeutic potential of targeting this pathway. Thus, our studies illuminate for the first time cytoskeletal-driven transcription as a personalized therapeutic target to combat drug resistant malignancies.

Gli3 is a negative regulator of Tas1r3-expressing taste cells

Mouse taste receptor cells survive from 3-24 days, necessitating their regeneration throughout adulthood. In anterior tongue, sonic hedgehog (SHH), released by a subpopulation of basal taste cells, regulates transcription factors Gli2 and Gli3 in stem cells to control taste cell regeneration. Using single-cell RNA-Seq we found that Gli3 is highly expressed in Tas1r3-expressing taste receptor cells and Lgr5+ taste stem cells in posterior tongue. By PCR and immunohistochemistry we found that Gli3 was expressed in taste buds in all taste fields. Conditional knockout mice lacking Gli3 in the posterior tongue (Gli3CKO) had larger taste buds containing more taste cells than did control wild-type (Gli3WT) mice. In comparison to wild-type mice, Gli3CKO mice had more Lgr5+ and Tas1r3+ cells, but fewer type III cells. Similar changes were observed ex vivo in Gli3CKO taste organoids cultured from Lgr5+ taste stem cells. Further, the expression of several taste marker and Gli3 target genes was altered in Gli3CKO mice and/or organoids. Mirroring these changes, Gli3CKO mice had increased lick responses to sweet and umami stimuli, decreased lick responses to bitter and sour taste stimuli, and increased glossopharyngeal taste nerve responses to sweet and bitter compounds. Our results indicate that Gli3 is a suppressor of stem cell proliferation that affects the number and function of mature taste cells, especially Tas1r3+ cells, in adult posterior tongue. Our findings shed light on the role of the Shh pathway in adult taste cell regeneration and may help devise strategies for treating taste distortions from chemotherapy and aging.

Hedgehog signaling drives medulloblastoma growth via CDK6

Medulloblastoma, an aggressive cancer of the cerebellum, is among the most common pediatric brain tumors. Approximately one-third of medulloblastomas are associated with misactivation of the Hedgehog (Hh) pathway. GLI family zinc finger 2 (GLI2) coordinates the Hh transcriptional program; however, the GLI2 targets that promote cancer cell proliferation are unknown. Here, we incorporated a Gli2-EGFP allele into 2 different genetic mouse models of Hh-associated medulloblastoma. Hh signaling induced GLI2 binding to the Cdk6 promoter and activated Cdk6 expression, thereby promoting uncontrolled cell proliferation. Genetic or pharmacological inhibition of CDK6 in mice repressed the growth of Hh-associated medulloblastoma and prolonged survival through inhibition of cell proliferation. In human medulloblastoma, misactivation of Hh signaling was associated with high levels of CDK6, pointing to CDK6 as a direct transcriptional target of the Hh pathway. These results suggest that CDK6 antagonists may be a promising therapeutic approach for Hh-associated medulloblastoma in humans.

John Saunders' ZPA, Sonic hedgehog and digit identity - How does it really all work?

Among John Saunders' many seminal contributions to developmental biology, his discovery of the limb 'zone of polarizing activity' (ZPA) is arguably one of the most memorable and ground-breaking. This discovery introduced the limb as a premier model for understanding developmental patterning and promoted the concept of patterning by a morphogen gradient. In the 50 years since the discovery of the ZPA, Sonic hedgehog (Shh) has been identified as the ZPA factor and the basic components of the signaling pathway and many aspects of its regulation have been elucidated. Although much has also been learned about how it regulates growth, the mechanism by which Shh patterns the limb, how it acts to instruct digit 'identity', nevertheless remains an enigma. This review focuses on what has been learned about Shh function in the limb and the outstanding puzzles that remain to be solved.

Towards quantitative analysis of gene regulation by enhancers

Enhancers are regulatory DNA sequences that can activate transcription over large distances. Recent studies have revealed the widespread role of distant activation in eukaryotic gene regulation and in the development of various human diseases, including cancer. Here we review recent progress in the field, focusing on new experimental and computational approaches that quantify the role of chromatin structure and dynamics during enhancer-promoter interactions in vitro and in vivo.

A SHH-FOXF1-BMP4 signaling axis regulating growth and differentiation of epithelial and mesenchymal tissues in ureter development

The differentiated cell types of the epithelial and mesenchymal tissue compartments of the mature ureter of the mouse arise in a precise temporal and spatial sequence from uncommitted precursor cells of the distal ureteric bud epithelium and its surrounding mesenchyme. Previous genetic efforts identified a member of the Hedgehog (HH) family of secreted proteins, Sonic hedgehog (SHH) as a crucial epithelial signal for growth and differentiation of the ureteric mesenchyme. Here, we used conditional loss- and gain-of-function experiments of the unique HH signal transducer Smoothened (SMO) to further characterize the cellular functions and unravel the effector genes of HH signaling in ureter development. We showed that HH signaling is not only required for proliferation and SMC differentiation of cells of the inner mesenchymal region but also for survival of cells of the outer mesenchymal region, and for epithelial proliferation and differentiation. We identified the Forkhead transcription factor gene Foxf1 as a target of HH signaling in the ureteric mesenchyme. Expression of a repressor version of FOXF1 in this tissue completely recapitulated the mesenchymal and epithelial proliferation and differentiation defects associated with loss of HH signaling while re-expression of a wildtype version of FOXF1 in the inner mesenchymal layer restored these cellular programs when HH signaling was inhibited. We further showed that expression of Bmp4 in the ureteric mesenchyme depends on HH signaling and Foxf1, and that exogenous BMP4 rescued cell proliferation and epithelial differentiation in ureters with abrogated HH signaling or FOXF1 function. We conclude that SHH uses a FOXF1-BMP4 module to coordinate the cellular programs for ureter elongation and differentiation, and suggest that deregulation of this signaling axis occurs in human congenital anomalies of the kidney and urinary tract (CAKUT).

The mammalian ureter is a simple tube with a specialized multi-layered epithelium, the urothelium, and a surrounding coat of fibroblasts and peristaltically active smooth muscle cells. Besides its important function in urinary drainage, the ureter represents a simple model system to study epithelial and mesenchymal tissue interactions in organ development. The differentiated cell types of the ureter coordinately arise from precursor cells of the distal ureteric bud and its surrounding mesenchyme. How their survival, growth and differentiation is regulated and coordinated within and between the epithelial and mesenchymal tissue compartments is largely unknown. Previous work identified Sonic hedgehog (SHH) as a crucial epithelial signal for growth and differentiation of the ureteric mesenchyme, but the entirety of the cellular functions and the molecular mediators of its mesenchymal signaling pathway have remained obscure. Here we showed that epithelial SHH acts in a paracrine fashion onto the ureteric mesenchyme to activate a FOXF1-BMP4 regulatory module that directs growth and differentiation of both ureteric tissue compartments. HH signaling additionally acts in outer mesenchymal cells as a survival factor. Thus, SHH is an epithelial signal that coordinates various cellular programs in early ureter development.

Neuropilin-1 promotes Hedgehog signaling through a novel cytoplasmic motif

Hedgehog (HH) signaling critically regulates embryonic and postnatal development as well as adult tissue homeostasis, and its perturbation can lead to developmental disorders, birth defects, and cancers. Neuropilins (NRPs), which have well-defined roles in Semaphorin and VEGF signaling, positively regulate HH pathway function, although their mechanism of action in HH signaling remains unclear. Here, using luciferase-based reporter assays, we provide evidence that NRP1 regulates HH signaling specifically at the level of GLI transcriptional activator function. Moreover, we show that NRP1 localization to the primary cilium, a key platform for HH signal transduction, does not correlate with HH signal promotion. Rather, a structure-function analysis suggests that the NRP1 cytoplasmic and transmembrane domains are necessary and sufficient to regulate HH pathway activity. Furthermore, we identify a previously uncharacterized, 12-amino acid region within the NRP1 cytoplasmic domain that mediates HH signal promotion. Overall, our results provide mechanistic insight into NRP1 function within and potentially beyond the HH signaling pathway. These insights have implications for the development of novel modulators of HH-driven developmental disorders and diseases.

The befuddling nature of mouse lemur hands and feet at Bezà Mahafaly, SW Madagascar

The reddish-gray mouse lemur (Microcebus griseorufus) possesses striking phenotypic and behavioral variation. This project investigates differences in autopod proportions in neighboring populations of M. griseorufus from the Special Reserve at Bezà Mahafaly in southwest Madagascar. One population resides in an environment generally preferred by M. griseorufus-a spiny forest with large-trunked trees, vertically-oriented supports, and more open ground, while the other resides in a gallery forest with abundant small, often horizontal peripheral branches in high canopy. We demonstrate significant interpopulation differences in autopod morphophology despite no evidence of divergence in mitochondrial cytochrome b. We test two hypotheses regarding ultimate causation. The first, based on the Fine Branch Arborealism Hypothesis (FBAH), holds that autopod differences are related to different locomotor practices in the two environments, and the second, based on the Narrow Niche Hypothesis (NNH), holds that the observed differences reflect a relaxation (from ancestral to descendant conditions) of selective pressure for terrestrial locomotion and/or use of large, vertical supports combined with positive selection for locomoting in peripheral branch settings. Our data conform well to FBAH expectations and show some support for the NNH. Individuals from the gallery forest possess disproportionally long posterior digits that facilitate locomotion on small, flexible canopy supports while individuals from the spiny forest possess shorter posterior digits and a longer pollex/hallux that increase functional grasping diameter for large vertical supports and facilitate efficient ground locomotion. Focal individual data confirm differences in how often individuals descend to the ground and use vertical supports. We further show that predispersal juveniles, like adults, possess autopod morphologies suited to their natal forest. We explore two proximate mechanisms that could generate these cheiridial differences. The first posits an in vivo plastic response to different locomotor behaviors, the second posits differences that manifest in early development.

Analysis of transcription factors expressed at the anterior mouse limb bud

Limb bud patterning, outgrowth, and differentiation are precisely regulated in a spatio-temporal manner through integrated networks of transcription factors, signaling molecules, and downstream genes. However, the exact mechanisms that orchestrate morphogenesis of the limb remain to be elucidated. Previously, we have established EMBRYS, a whole-mount in situ hybridization database of transcription factors. Based on the findings from EMBRYS, we focused our expression pattern analysis on a selection of transcription factor genes that exhibit spatially localized and temporally dynamic expression patterns with respect to the anterior-posterior axis in the E9.5–E11.5 limb bud. Among these genes, Irx3 showed a posteriorly expanded expression domain in Shh^-/- limb buds and an anteriorly reduced expression domain in Gli3^-/- limb buds, suggesting their importance in anterior-posterior patterning. To assess the stepwise EMBRYS-based screening system for anterior regulators, we generated Irx3 transgenic mice in which Irx3 was expressed in the entire limb mesenchyme under the Prrx1 regulatory element. The Irx3 gain-of-function model displayed complex phenotypes in the autopods, including digit loss, radial flexion, and fusion of the metacarpal bones, suggesting that Irx3 may contribute to the regulation of limb patterning, especially in the autopods. Our results demonstrate that gene expression analysis based on EMBRYS could contribute to the identification of genes that play a role in patterning of the limb mesenchyme.

Cilia-dependent GLI processing in neural crest cells is required for tongue development

Ciliopathies are a class of diseases caused by the loss of a ubiquitous, microtubule-based organelle called a primary cilium. Ciliopathies commonly result in defective development of the craniofacial complex, causing midfacial defects, craniosynostosis, micrognathia and aglossia. Herein, we explored how the conditional loss of primary cilia on neural crest cells (Kif3a^f/f;Wnt1-Cre) generated aglossia. On a cellular level, our data revealed that aglossia in Kif3a^f/f;Wnt1-Cre embryos was due to a loss of mesoderm-derived muscle precursors migrating into and surviving in the tongue anlage. To determine the molecular basis for this phenotype, we performed RNA-seq, in situ hybridization, qPCR and Western blot analyses. We found that transduction of the Sonic hedgehog (Shh) pathway, rather than other pathways previously implicated in tongue development, was aberrant in Kif3a^f/f;Wnt1-Cre embryos. Despite increased production of full-length GLI2 and GLI3 isoforms, previously identified GLI targets important for mandibular and glossal development (Foxf1, Foxf2, Foxd1 and Foxd2) were transcriptionally downregulated in Kif3a^f/f;Wnt1-Cre embryos. Genetic removal of GLI activator (GLIA) isoforms in neural crest cells recapitulated the aglossia phenotype and downregulated Fox gene expression. Genetic addition of GLIA isoforms in neural crest cells partially rescued the aglossia phenotype and Fox gene expression in Kif3a^f/f;Wnt1-Cre embryos. Together, our data suggested that glossal development requires primary cilia-dependent GLIA activity in neural crest cells. Furthermore, these data, in conjunction with our previous work, suggested prominence specific roles for GLI isoforms; with development of the frontonasal prominence relying heavily on the repressor isoform and the development of the mandibular prominence/tongue relying heavily on the activator isoform.

Hedgehog Signaling: From Basic Biology to Cancer Therapy

The Hedgehog (HH) signaling pathway was discovered originally as a key pathway in embryonic patterning and development. Since its discovery, it has become increasingly clear that the HH pathway also plays important roles in a multitude of cancers. Therefore, HH signaling has emerged as a therapeutic target of interest for cancer therapy. In this review, we provide a brief overview of HH signaling and the key molecular players involved and offer an up-to-date summary of our current knowledge of endogenous and exogenous small molecules that modulate HH signaling. We discuss experiences and lessons learned from the decades-long efforts toward the development of cancer therapies targeting the HH pathway. Challenges to develop next-generation cancer therapies are highlighted.

Integration of Shh and Fgf signaling in controlling Hox gene expression in cultured limb cells

During embryonic development, fields of progenitor cells form complex structures through dynamic interactions with external signaling molecules. How complex signaling inputs are integrated to yield appropriate gene expression responses is poorly understood. In the early limb bud, for instance, Sonic hedgehog (Shh) is expressed in the distal posterior mesenchyme, where it acts as a mediator of anterior to posterior (AP) patterning, whereas fibroblast growth factor 8 (Fgf8) is produced by the apical ectodermal ridge (AER) at the distal tip of the limb bud to direct outgrowth along the proximal to distal (PD) axis. Here we use cultured limb mesenchyme cells to assess the response of the target Hoxd genes to these two factors. We find that they act synergistically and that both factors are required to activate Hoxd13 in limb mesenchymal cells. However, the analysis of the enhancer landscapes flanking the HoxD cluster reveals that the bimodal regulatory switch observed in vivo is only partially achieved under these in vitro conditions, suggesting an additional requirement for other factors.

IFT56 regulates vertebrate developmental patterning by maintaining IFTB complex integrity and ciliary microtubule architecture

Cilia are key regulators of animal development and depend on intraflagellar transport (IFT) proteins for their formation and function, yet the roles of individual IFT proteins remain unclear. We examined the Ift56^hop mouse mutant and reveal novel insight into the function of IFT56, a poorly understood IFTB protein. Ift56^hop mice have normal cilia distribution but display defective cilia structure, including abnormal positioning and number of ciliary microtubule doublets. We show that Ift56^hop cilia are unable to accumulate Gli proteins efficiently, resulting in developmental patterning defects in Shh signaling-dependent tissues such as the limb and neural tube. Strikingly, core IFTB proteins are unable to accumulate normally within Ift56^hop cilia, including IFT88, IFT81 and IFT27, which are crucial for key processes such as tubulin transport and Shh signaling. IFT56 is required specifically for the IFTB complex, as IFTA components and proteins that rely on IFTA function are unaffected in Ift56^hop cilia. These studies define a distinct and novel role for IFT56 in IFTB complex integrity that is crucial for cilia structure and function and, ultimately, animal development.

Sonic Hedgehog Signaling in Limb Development

The gene encoding the secreted protein Sonic hedgehog (Shh) is expressed in the polarizing region (or zone of polarizing activity), a small group of mesenchyme cells at the posterior margin of the vertebrate limb bud. Detailed analyses have revealed that Shh has the properties of the long sought after polarizing region morphogen that specifies positional values across the antero-posterior axis (e.g., thumb to little finger axis) of the limb. Shh has also been shown to control the width of the limb bud by stimulating mesenchyme cell proliferation and by regulating the antero-posterior length of the apical ectodermal ridge, the signaling region required for limb bud outgrowth and the laying down of structures along the proximo-distal axis (e.g., shoulder to digits axis) of the limb. It has been shown that Shh signaling can specify antero-posterior positional values in limb buds in both a concentration- (paracrine) and time-dependent (autocrine) fashion. Currently there are several models for how Shh specifies positional values over time in the limb buds of chick and mouse embryos and how this is integrated with growth. Extensive work has elucidated downstream transcriptional targets of Shh signaling. Nevertheless, it remains unclear how antero-posterior positional values are encoded and then interpreted to give the particular structure appropriate to that position, for example, the type of digit. A distant cis-regulatory enhancer controls limb-bud-specific expression of Shh and the discovery of increasing numbers of interacting transcription factors indicate complex spatiotemporal regulation. Altered Shh signaling is implicated in clinical conditions with congenital limb defects and in the evolution of the morphological diversity of vertebrate limbs.

Limb development: a paradigm of gene regulation

The limb is a commonly used model system for developmental biology. Given the need for precise control of complex signalling pathways to achieve proper patterning, the limb is also becoming a model system for gene regulation studies. Recent developments in genomic technologies have enabled the genome-wide identification of regulatory elements that control limb development, yielding insights into the determination of limb morphology and forelimb versus hindlimb identity. The modulation of regulatory interactions - for example, through the modification of regulatory sequences or chromatin architecture - can lead to morphological evolution, acquired regeneration capacity or limb malformations in diverse species, including humans.

Suppressor of Fused Chaperones Gli Proteins To Generate Transcriptional Responses to Sonic Hedgehog Signaling

Cellular responses to the graded Sonic Hedgehog (Shh) morphogenic signal are orchestrated by three Gli genes that give rise to both transcription activators and repressors. An essential downstream regulator of the pathway, encoded by the tumor suppressor gene Suppressor of fused (Sufu), plays critical roles in the production, trafficking, and function of Gli proteins, but the mechanism remains controversial. Here, we show that Sufu is upregulated in active Shh responding tissues and accompanies Gli activators translocating into and Gli repressors out of the nucleus. Trafficking of Sufu to the primary cilium, potentiated by Gli activators but not repressors, was found to be coupled to its nuclear import. We have identified a nuclear export signal (NES) motif of Sufu in juxtaposition to the protein kinase A (PKA) and glycogen synthase kinase 3 (GSK3) dual phosphorylation sites and show that Sufu binds the chromatin with both Gli1 and Gli3. Close comparison of neural tube development among individual Ptch1-/-, Sufu-/-, and Ptch1-/-; Sufu-/- double mutant embryos indicates that Sufu is critical for the maximal activation of Shh signaling essential to the specification of the most-ventral neurons. These data define Sufu as a novel class of molecular chaperone required for every aspect of Gli regulation and function.

Understanding Mechanisms of GLI-Mediated Transcription during Craniofacial Development and Disease Using the Ciliopathic Mutant, talpid2

The primary cilium is a ubiquitous, microtubule-based organelle that cells utilize to transduce molecular signals. Ciliopathies are a group of diseases that are caused by a disruption in the structure or function of the primary cilium. Over 30% of all ciliopathies are primarily defined by their craniofacial phenotypes, which typically include midfacial defects, cleft lip/palate, micrognathia, aglossia, and craniosynostosis. The frequency and severity of craniofacial phenotypes in ciliopathies emphasizes the importance of the cilium during development of the craniofacial complex. Molecularly, many ciliopathic mutants, including the avian talpid² (ta²), report pathologically high levels of full-length GLI3 (GLI3_FL), which can go on to function as an activator (GLI_A), and reduced production of truncated GLI3 (GLI3_T), which can go on to function as a repressor (GLI_R). These observations suggest that the craniofacial phenotypes of ciliary mutants like ta² are caused either by excessive activity of the GLI_A or reduced activity of GLI_R. To decipher between these two scenarios, we examined GLI3 occupation at the regulatory regions of target genes and subsequent target gene expression. Using in silico strategies we identified consensus GLI binding regions (GBRs) in the avian genome and confirmed GLI3 binding to the regulatory regions of its targets by chromatin immunoprecipitation (ChIP). In ta² mutants, there was a strikingly low number of GLI3 target genes that had significantly increased expression in facial prominences compared to the control embryo and GLI3 occupancy at GBRs associated with target genes was largely reduced. In vitro DNA binding assays, further supported ChIP results, indicated that the excessive GLI3_FL generated in ta² mutants did not bind to GBRs. In light of these results, we explored the possibility of GLI co-regulator proteins playing a role in regulatory mechanism of GLI-mediated transcription. Taken together our studies suggest that craniofacial ciliopathic phenotypes are produced via reduced GLI_T production, allowing for target gene transcription to be mediated by the combinatorial code of GLI co-regulators.

Explaining the disease phenotype of intergenic SNP through predicted long range regulation

Thousands of disease-associated SNPs (daSNPs) are located in intergenic regions (IGR), making it difficult to understand their association with disease phenotypes. Recent analysis found that non-coding daSNPs were frequently located in or approximate to regulatory elements, inspiring us to try to explain the disease phenotypes of IGR daSNPs through nearby regulatory sequences. Hence, after locating the nearest distal regulatory element (DRE) to a given IGR daSNP, we applied a computational method named INTREPID to predict the target genes regulated by the DRE, and then investigated their functional relevance to the IGR daSNP's disease phenotypes. 36.8% of all IGR daSNP-disease phenotype associations investigated were possibly explainable through the predicted target genes, which were enriched with, were functionally relevant to, or consisted of the corresponding disease genes. This proportion could be further increased to 60.5% if the LD SNPs of daSNPs were also considered. Furthermore, the predicted SNP-target gene pairs were enriched with known eQTL/mQTL SNP-gene relationships. Overall, it's likely that IGR daSNPs may contribute to disease phenotypes by interfering with the regulatory function of their nearby DREs and causing abnormal expression of disease genes.

Transcription factor PRO1 targets genes encoding conserved components of fungal developmental signaling pathways

The filamentous fungus Sordaria macrospora is a model system to study multicellular development during fruiting body formation. Previously, we demonstrated that this major process in the sexual life cycle is controlled by the Zn(II)₂ Cys₆ zinc cluster transcription factor PRO1. Here, we further investigated the genome-wide regulatory network controlled by PRO1 by employing chromatin immunoprecipitation combined with next-generation sequencing (ChIP-seq) to identify binding sites for PRO1. We identified several target regions that occur in the promoter regions of genes encoding components of diverse signaling pathways. Furthermore, we identified a conserved DNA-binding motif that is bound specifically by PRO1 in vitro. In addition, PRO1 controls in vivo the expression of a DsRed reporter gene under the control of the esdC target gene promoter. Our ChIP-seq data suggest that PRO1 also controls target genes previously shown to be involved in regulating the pathways controlling cell wall integrity, NADPH oxidase and pheromone signaling. Our data point to PRO1 acting as a master regulator of genes for signaling components that comprise a developmental cascade controlling fruiting body formation.

A unique stylopod patterning mechanism by Shox2-controlled osteogenesis

Vertebrate appendage patterning is programmed by Hox-TALE factor-bound regulatory elements. However, it remains unclear which cell lineages are commissioned by Hox-TALE factors to generate regional specific patterns and whether other Hox-TALE co-factors exist. In this study, we investigated the transcriptional mechanisms controlled by the Shox2 transcriptional regulator in limb patterning. Harnessing an osteogenic lineage-specific Shox2 inactivation approach we show that despite widespread Shox2 expression in multiple cell lineages, lack of the stylopod observed upon Shox2 deficiency is a specific result of Shox2 loss of function in the osteogenic lineage. ChIP-Seq revealed robust interaction of Shox2 with cis-regulatory enhancers clustering around skeletogenic genes that are also bound by Hox-TALE factors, supporting a lineage autonomous function of Shox2 in osteogenic lineage fate determination and skeleton patterning. Pbx ChIP-Seq further allowed the genome-wide identification of cis-regulatory modules exhibiting co-occupancy of Pbx, Meis and Shox2 transcriptional regulators. Integrative analysis of ChIP-Seq and RNA-Seq data and transgenic enhancer assays indicate that Shox2 patterns the stylopod as a repressor via interaction with enhancers active in the proximal limb mesenchyme and antagonizes the repressive function of TALE factors in osteogenesis.

Genome-Wide Chromatin Immunoprecipitation Sequencing Analysis of the Penicillium chrysogenum Velvet Protein PcVelA Identifies Methyltransferase PcLlmA as a Novel Downstream Regulator of Fungal Development

Penicillium chrysogenum is the sole industrial producer of the β-lactam antibiotic penicillin, which is the most commonly used drug for treating bacterial infections. In P. chrysogenum and other filamentous fungi, secondary metabolism and morphogenesis are controlled by the highly conserved multisubunit velvet complex. Here we present the first chromatin immunoprecipitation next-generation sequencing (ChIP-seq) analysis of a fungal velvet protein, providing experimental evidence that a velvet homologue in P. chrysogenum (PcVelA) acts as a direct transcriptional regulator at the DNA level in addition to functioning as a regulator at the protein level in P. chrysogenum, which was previously described. We identified many target genes that are related to processes known to be dependent on PcVelA, e.g., secondary metabolism as well as asexual and sexual development. We also identified seven PcVelA target genes that encode putative methyltransferases. Yeast two-hybrid and bimolecular fluorescence complementation analyses showed that one of the putative methyltransferases, PcLlmA, directly interacts with PcVelA. Furthermore, functional characterization of PcLlmA demonstrated that this protein is involved in the regulation of conidiosporogenesis, pellet formation, and hyphal morphology, all traits with major biotechnological relevance. IMPORTANCE Filamentous fungi are of major interest for biotechnological and pharmaceutical applications. This is due mainly to their ability to produce a wide variety of secondary metabolites, many of which are relevant as antibiotics. One of the most prominent examples is penicillin, a β-lactam antibiotic that is produced on the industrial scale by fermentation of P. chrysogenum. In recent years, the multisubunit protein complex velvet has been identified as one of the key regulators of fungal secondary metabolism and development. However, until recently, only a little has been known about how velvet mediates regulation at the molecular level. To address this issue, we performed ChIP-seq (chromatin immunoprecipitation in combination with next-generation sequencing) on and follow-up analysis of PcVelA, the core component of the velvet complex in P. chrysogenum. We demonstrate direct involvement of velvet in transcriptional control and present the putative methyltransferase PcLlmA as a new downstream factor and interaction partner of PcVelA.

Gene regulation at a distance: From remote enhancers to 3D regulatory ensembles

Large-scale identification of elements associated with gene expression revealed that many of them are located extremely far from gene transcriptional start sites. We review here the growing evidence that show that distal cis-acting elements provide key instructions to genes, as genetic variation affecting them is growingly identified as an importance source of phenotypic diversity and disease. We discuss the different mechanisms that allow these elements to exert their regulatory functions, in a robust and specific manner, despite the large genomic distances separating them from their target genes. We particularly focus on the role of the structural organization of the genome in guiding such regulatory interactions.

T-box3 is a ciliary protein and regulates stability of the Gli3 transcription factor to control digit number

Crucial roles for T-box3 in development are evident by severe limb malformations and other birth defects caused by T-box3 mutations in humans. Mechanisms whereby T-box3 regulates limb development are poorly understood. We discovered requirements for T-box at multiple stages of mouse limb development and distinct molecular functions in different tissue compartments. Early loss of T-box3 disrupts limb initiation, causing limb defects that phenocopy Sonic Hedgehog (Shh) mutants. Later ablation of T-box3 in posterior limb mesenchyme causes digit loss. In contrast, loss of anterior T-box3 results in preaxial polydactyly, as seen with dysfunction of primary cilia or Gli3-repressor. Remarkably, T-box3 is present in primary cilia where it colocalizes with Gli3. T-box3 interacts with Kif7 and is required for normal stoichiometry and function of a Kif7/Sufu complex that regulates Gli3 stability and processing. Thus, T-box3 controls digit number upstream of Shh-dependent (posterior mesenchyme) and Shh-independent, cilium-based (anterior mesenchyme) Hedgehog pathway function.

DOI:http://dx.doi.org/10.7554/eLife.07897.001

Mutations in the gene that encodes a protein called T-box3 cause serious birth defects, including deformities of the hands and feet, via poorly understood mechanisms. Several other proteins are also important for ensuring that limbs develop correctly. These include the Sonic Hedgehog protein, which controls a signaling pathway that determines whether a protein called Gli3 is converted into its “repressor” form. The hair-like structures called primary cilia that sit on the surface of animal cells also contain Gli3, and processes within these structures control the production of the Gli3-repressor.

Emechebe, Kumar et al. have now studied genetically engineered mice in which the production of the T-box3 protein was stopped at different stages of mouse development. This revealed that turning off T-box3 production early in development causes many parts of the limb not to form. This type of defect appears to be the same as that seen in mice that lack the Sonic Hedgehog protein.

If the production of T-box3 is turned off later in mouse development in the rear portion of the developing limb, the limb starts to develop but doesn’t develop enough rear toes. When T-box3 production is turned off in the front portion of the developing limbs, mice are born with too many front toes. This latter problem mimics the effects seen in mice that are unable to produce Gli3-repressor or that have defective primary cilia.

Further investigation unexpectedly revealed that T-box3 is found in primary cilia and localizes to the same regions of the cilia as the Gli3-repressor. Furthermore, T-box3 also interacts with a protein complex that controls the stability of Gli3 and processes it into the Gli3-repressor form.

In the future, it will be important to determine how T-box3 controls the stability of Gli3 and whether that process occurs in the primary cilia or in other parts of the cell where T-box3 and Gli3 coexist, such as the nucleus. This could help us understand how T-box3 and Sonic Hedgehog signaling contribute to other aspects of development and to certain types of cancer.

DOI:http://dx.doi.org/10.7554/eLife.07897.002

Anteroposterior Limb Skeletal Patterning Requires the Bifunctional Action of SWI/SNF Chromatin Remodeling Complex in Hedgehog Pathway

Graded Sonic hedgehog (Shh) signaling governs vertebrate limb skeletal patterning along the anteroposterior (AP) axis by regulating the activity of bifunctional Gli transcriptional regulators. The genetic networks involved in this patterning are well defined, however, the epigenetic control of the process by chromatin remodelers remains unknown. Here, we report that the SWI/SNF chromatin remodeling complex is essential for Shh-driven limb AP patterning. Specific inactivation of Srg3/mBaf155, a core subunit of the remodeling complex, in developing limb buds hampered the transcriptional upregulation of Shh/Gli target genes, including the Shh receptor Ptch1 and its downstream effector Gli1 in the posterior limb bud. In addition, Srg3 deficiency induced ectopic activation of the Hedgehog (Hh) pathway in the anterior mesenchyme, resulting in loss of progressive asymmetry. These defects in the Hh pathway accompanied aberrant BMP activity and disruption of chondrogenic differentiation in zeugopod and autopod primordia. Notably, our data revealed that dual control of the Hh pathway by the SWI/SNF complex is essential for spatiotemporal transcriptional regulation of the BMP antagonist Gremlin1, which affects the onset of chondrogenesis. This study uncovers the bifunctional role of the SWI/SNF complex in the Hh pathway to determine the fate of AP skeletal progenitors.

Anteroposterior (AP) limb skeletal patterning is directed by morphogen Sonic hedgehog (Shh) signaling. Modulation of Shh responsiveness and repression of Shh pathway activity in distinct limb bud regions are essential for proper limb skeletal formation. Although the genetic networks involved in these processes have been identified, epigenetic control by chromatin remodeler remains unknown. We have unraveled the function of the SWI/SNF chromatin remodeling complex in Shh signaling during limb patterning. The complex activates the responses of the posterior limb progenitors to Shh, however, it represses the signaling in the anterior limb progenitors. Here we provide genetic evidence for the dual requirement of the SWI/SNF complex in Shh signaling to pattern AP limb skeletal elements.

Hedgehog inhibition causes complete loss of limb outgrowth and transformation of digit identity in Xenopus tropicalis

The study of the tetrapod limb has contributed greatly to our understanding of developmental pathways and how changes to these pathways affect the evolution of morphology. Most of our understanding of tetrapod limb development comes from research on amniotes, with far less known about mechanisms of limb development in amphibians. To better understand the mechanisms of limb development in anuran amphibians, we used cyclopamine to inhibit Hedgehog signaling at various stages of development in the western clawed frog, Xenopus tropicalis, and observed resulting morphologies. We also analyzed gene expression changes resulting from similar experiments in Xenopus laevis. Inhibition of Hedgehog signaling in X. tropicalis results in limb abnormalities including reduced digit number, missing skeletal elements, and complete absence of limbs. In addition, posterior digits assume an anterior identity by developing claws that are usually only found on anterior digits, confirming Sonic hedgehog's role in digit identity determination. Thus, Sonic hedgehog appears to play mechanistically separable roles in digit number specification and digit identity specification as in other studied tetrapods. The complete limb loss observed in response to reduced Hedgehog signaling in X. tropicalis, however, is striking, as this functional role for Hedgehog signaling has not been found in any other tetrapod. This changed mechanism may represent a substantial developmental constraint to digit number evolution in frogs. J. Exp. Zool. (Mol. Dev. Evol.) 9999B:XX-XX, 2016. © 2016 Wiley Periodicals, Inc.

Cyclopamine tartrate, an inhibitor of Hedgehog signaling, strongly interferes with mitochondrial function and suppresses aerobic respiration in lung cancer cells

BACKGROUND

Aberrant Hedgehog (Hh) signaling is associated with the development of many cancers including prostate cancer, gastrointestinal cancer, lung cancer, pancreatic cancer, ovarian cancer, and basal cell carcinoma. The Hh signaling pathway has been one of the most intensely investigated targets for cancer therapy, and a number of compounds inhibiting Hh signaling are being tested clinically for treating many cancers. Lung cancer causes more deaths than the next three most common cancers (colon, breast, and prostate) combined. Cyclopamine was the first compound found to inhibit Hh signaling and has been invaluable for understanding the function of Hh signaling in development and cancer. To find novel strategies for combating lung cancer, we decided to characterize the effect of cyclopamine tartrate (CycT), an improved analogue of cyclopamine, on lung cancer cells and its mechanism of action.

METHODS

The effect of CycT on oxygen consumption and proliferation of non-small-cell lung cancer (NSCLC) cell lines was quantified by using an Oxygraph system and live cell counting, respectively. Apoptosis was detected by using Annexin V and Propidium Iodide staining. CycT's impact on ROS generation, mitochondrial membrane potential, and mitochondrial morphology in NSCLC cells was monitored by using fluorometry and fluorescent microscopy. Western blotting and fluorescent microscopy were used to detect the levels and localization of Hh signaling targets, mitochondrial fission protein Drp1, and heme-related proteins in various NSCLC cells.

RESULTS

Our findings identified a novel function of CycT, as well as another Hh inhibitor SANT1, to disrupt mitochondrial function and aerobic respiration. Our results showed that CycT, like glutamine depletion, caused a substantial decrease in oxygen consumption in a number of NSCLC cell lines, suppressed NSCLC cell proliferation, and induced apoptosis. Further, we found that CycT increased ROS generation, mitochondrial membrane hyperpolarization, and mitochondrial fragmentation, thereby disrupting mitochondrial function in NSCLC cells.

CONCLUSIONS

Together, our work demonstrates that CycT, and likely other Hh signaling inhibitors, can interrupt NSCLC cell function by promoting mitochondrial fission and fragmentation, mitochondrial membrane hyperpolarization, and ROS generation, thereby diminishing mitochondrial respiration, suppressing cell proliferation, and causing apoptosis. Our work provides novel mechanistic insights into the action of Hh inhibitors in cancer cells.

Computational prediction and experimental validation of novel Hedgehog-responsive enhancers linked to genes of the Hedgehog pathway

BACKGROUND

The Hedgehog (Hh) signaling pathway, acting through three homologous transcription factors (GLI1, GLI2, GLI3) in vertebrates, plays multiple roles in embryonic organ development and adult tissue homeostasis. At the level of the genome, GLI factors bind to specific motifs in enhancers, some of which are hundreds of kilobases removed from the gene promoter. These enhancers integrate the Hh signal in a context-specific manner to control the spatiotemporal pattern of target gene expression. Importantly, a number of genes that encode Hh pathway molecules are themselves targets of Hh signaling, allowing pathway regulation by an intricate balance of feed-back activation and inhibition. However, surprisingly few of the critical enhancer elements that control these pathway target genes have been identified despite the fact that such elements are central determinants of Hh signaling activity. Recently, ChIP studies have been carried out in multiple tissue contexts using mouse models carrying FLAG-tagged GLI proteins (GLI(FLAG)). Using these datasets, we tested whether a meta-analysis of GLI binding sites, coupled with a machine learning approach, could reveal genomic features that could be used to empirically identify Hh-regulated enhancers linked to loci of the Hh signaling pathway.

RESULTS

A meta-analysis of four existing GLI(FLAG) datasets revealed a library of GLI binding motifs that was substantially more restricted than the potential sites predicted by previous in vitro binding studies. A machine learning method (kmer-SVM) was then applied to these datasets and enriched k-mers were identified that, when applied to the mouse genome, predicted as many as 37,000 potential Hh enhancers. For functional analysis, we selected nine regions which were annotated to putative Hh pathway molecules and found that seven exhibited GLI-dependent activity, indicating that they are directly regulated by Hh signaling (78% success rate).

CONCLUSIONS

The results suggest that Hh enhancer regions share common sequence features. The kmer-SVM machine learning approach identifies those features and can successfully predict functional Hh regulatory regions in genomic DNA surrounding Hh pathway molecules and likely, other Hh targets. Additionally, the library of enriched GLI binding motifs that we have identified may allow improved identification of functional GLI binding sites.

The Pluripotency Factor NANOG Binds to GLI Proteins and Represses Hedgehog-mediated Transcription

The Hedgehog (HH) signaling pathway is essential for the maintenance and response of several types of stem cells. To study the transcriptional response of stem cells to HH signaling, we searched for proteins binding to GLI proteins, the transcriptional effectors of the HH pathway in mouse embryonic stem (ES) cells. We found that both GLI3 and GLI1 bind to the pluripotency factor NANOG. The ectopic expression of NANOG inhibits GLI1-mediated transcriptional responses in a dose-dependent fashion. In differentiating ES cells, the presence of NANOG reduces the transcriptional response of cells to HH. Finally, we found thatGli1andNanogare co-expressed in ES cells at high levels. We propose that NANOG acts as a negative feedback component that provides stem cell-specific regulation of the HH pathway.

Identification and Validation of Novel Hedgehog-Responsive Enhancers Predicted by Computational Analysis of Ci/Gli Binding Site Density

The Hedgehog (Hh) signaling pathway directs a multitude of cellular responses during embryogenesis and adult tissue homeostasis. Stimulation of the pathway results in activation of Hh target genes by the transcription factor Ci/Gli, which binds to specific motifs in genomic enhancers. In Drosophila, only a few enhancers (patched, decapentaplegic, wingless, stripe, knot, hairy, orthodenticle) have been shown by in vivo functional assays to depend on direct Ci/Gli regulation. All but one (orthodenticle) contain more than one Ci/Gli site, prompting us to directly test whether homotypic clustering of Ci/Gli binding sites is sufficient to define a Hh-regulated enhancer. We therefore developed a computational algorithm to identify Ci/Gli clusters that are enriched over random expectation, within a given region of the genome. Candidate genomic regions containing Ci/Gli clusters were functionally tested in chicken neural tube electroporation assays and in transgenic flies. Of the 22 Ci/Gli clusters tested, seven novel enhancers (and the previously known patched enhancer) were identified as Hh-responsive and Ci/Gli-dependent in one or both of these assays, including: Cuticular protein 100A (Cpr100A); invected (inv), which encodes an engrailed-related transcription factor expressed at the anterior/posterior wing disc boundary; roadkill (rdx), the fly homolog of vertebrate Spop; the segment polarity gene gooseberry (gsb); and two previously untested regions of the Hh receptor-encoding patched (ptc) gene. We conclude that homotypic Ci/Gli clustering is not sufficient information to ensure Hh-responsiveness; however, it can provide a clue for enhancer recognition within putative Hedgehog target gene loci.

Urogenital development in Pallister-Hall syndrome is disrupted in a cell-lineage-specific manner by constitutive expression of GLI3 repressor

Pallister-Hall syndrome (PHS) is a rare disorder caused by mutations in GLI3 that produce a transcriptional repressor (GLI3R). Individuals with PHS present with a variably penetrant variety of urogenital system malformations, including renal aplasia or hypoplasia, hydroureter, hydronephrosis or a common urogenital sinus. The embryologic mechanisms controlled by GLI3R that result in these pathologic phenotypes are undefined. We demonstrate that germline expression of GLI3R causes renal hypoplasia, associated with decreased nephron number, and hydroureter and hydronephrosis, caused by blind-ending ureters. Mice with obligate GLI3R expression also displayed duplication of the ureters that was caused by aberrant common nephric duct patterning and ureteric stalk outgrowth. These developmental abnormalities are associated with suppressed Hedgehog signaling activity in the cloaca and adjacent vesicular mesenchyme. Mice with conditional expression of GLI3R were utilized to identify lineage-specific effects of GLI3R. In the ureteric bud, GLI3R expression decreased branching morphogenesis. In Six2-positive nephrogenic progenitors, GLI3R decreased progenitor cell proliferation reducing the number of nephrogenic precursor structures. Using mutant mice with Gli3R and Gli3 null alleles, we demonstrate that urogenital system patterning and development is controlled by the levels of GLI3R and not by an absence of full-length GLI3. We conclude that the urogenital system phenotypes observed in PHS are caused by GLI3R-dependent perturbations in nephric duct patterning, renal branching morphogenesis and nephrogenic progenitor self-renewal.

Next generation limb development and evolution: old questions, new perspectives

The molecular analysis of limb bud development in vertebrates continues to fuel our understanding of the gene regulatory networks that orchestrate the patterning, proliferation and differentiation of embryonic progenitor cells. In recent years, systems biology approaches have moved our understanding of the molecular control of limb organogenesis to the next level by incorporating next generation ‘omics’ approaches, analyses of chromatin architecture, enhancer-promoter interactions and gene network simulations based on quantitative datasets into experimental analyses. This Review focuses on the insights these studies have given into the gene regulatory networks that govern limb development and into the fin-to-limb transition and digit reductions that occurred during the evolutionary diversification of tetrapod limbs.

A shift in anterior-posterior positional information underlies the fin-to-limb evolution

The pectoral fins of ancestral fishes had multiple proximal elements connected to their pectoral girdles. During the fin-to-limb transition, anterior proximal elements were lost and only the most posterior one remained as the humerus. Thus, we hypothesised that an evolutionary alteration occurred in the anterior–posterior (AP) patterning system of limb buds. In this study, we examined the pectoral fin development of catshark (Scyliorhinus canicula) and revealed that the AP positional values in fin buds are shifted more posteriorly than mouse limb buds. Furthermore, examination of Gli3 function and regulation shows that catshark fins lack a specific AP patterning mechanism, which restricts its expression to an anterior domain in tetrapods. Finally, experimental perturbation of AP patterning in catshark fin buds results in an expansion of posterior values and loss of anterior skeletal elements. Together, these results suggest that a key genetic event of the fin-to-limb transformation was alteration of the AP patterning network.

DOI:http://dx.doi.org/10.7554/eLife.07048.001

Humans, mice, and other animals with four limbs belong to a group of land-dwelling animals known as the tetrapods. This group of animals evolved from ancient fish and one crucial adaptation to life on land involved the modification of fins to form limbs. The front pair of limbs (the ‘arms’) evolved from the ‘pectoral’ fins of the ancient fish. These fins contain numerous bones that fan out from a set of bones called the pectoral girdle. However, most of the bones nearer the front side (the thumb side in the human limb) were lost in the ancestors of tetrapods as they moved onto land. Only the bone nearest the back remained as the ‘humerus’, which forms the upper part of the limb (i.e., the upper arm of humans).

In the embryos of mice and other animals, the limbs develop from structures called limb buds. For the limb to develop properly, the cells in the limb bud need to receive specific instructions that depend on their position in the bud. A protein called Gli3R provides cells with information about their position along the ‘anterior–posterior’ (or thumb-to-little finger) axis of the bud. This protein regulates several genes that are involved in limb development, and this results in different genes being expressed in cells along the anterior–posterior axis. For example, Alx4 is only expressed in a small area at the anterior end of the bud, while Hand2 expression is found in a large area towards the posterior part.

Gli3R is also found in a fish called the catshark, but it is not clear how it controls the formation of fins. Onimaru et al. show that the pattern of gene expression in the catshark fin bud is different to that of the mouse limb bud. For example, Alx4 is expressed in a larger area of the fin bud that extends further towards the posterior, while Hand2 is only found in a much smaller area at the posterior end of the bud. The experiments also suggest that Gli3R is active in a much larger area of the fin bud than in the limb bud.

Next, Onimaru et al. used a drug on the catshark embryos to increase the activity of another protein that can inhibit Gli3R. The fin buds of these shark had anterior shift in several gene expression domains, and the fins that formed were missing several anterior bones and had only a single bone connected to the pectoral girdle. Onimaru et al.'s findings suggest that during the evolution of the tetrapods, there may have been a shift in the anterior–posterior patterning of the fin bud to form a limb. An important area for future work will be to use genome-wide studies to study the fin/limb buds of other species.

DOI:http://dx.doi.org/10.7554/eLife.07048.002

Hedgehog Signaling in the Maintenance of Cancer Stem Cells

Cancer stem cells (CSCs) represent a rare population of cells with the capacity to self-renew and give rise to heterogeneous cell lineages within a tumour. Whilst the mechanisms underlying the regulation of CSCs are poorly defined, key developmental signaling pathways required for normal stem and progenitor functions have been strongly implicated. Hedgehog (Hh) signaling is an evolutionarily-conserved pathway essential for self-renewal and cell fate determination. Aberrant Hh signaling is associated with the development and progression of various types of cancer and is implicated in multiple aspects of tumourigenesis, including the maintenance of CSCs. Here, we discuss the mounting evidence suggestive of Hh-driven CSCs in the context of haematological malignancies and solid tumours and the novel strategies that hold the potential to block many aspects of the transformation attributed to the CSC phenotype, including chemotherapeutic resistance, relapse and metastasis.