SOX11 regulates SWI/SNF complex components as member of the adrenergic neuroblastoma core regulatory circuitry.

The pediatric extra-cranial tumor neuroblastoma displays a low mutational burden while recurrent copy number alterations are present in most high-risk cases. Here, we identify SOX11 as a dependency transcription factor in adrenergic neuroblastoma based on recurrent chromosome 2p focal gains and amplifications, specific expression in the normal sympatho-adrenal lineage and adrenergic neuroblastoma, regulation by multiple adrenergic specific (super-)enhancers and strong dependency on high SOX11 expression in adrenergic neuroblastomas. SOX11 regulated direct targets include genes implicated in epigenetic control, cytoskeleton and neurodevelopment. Most notably, SOX11 controls chromatin regulatory complexes, including 10 SWI/SNF core components among which SMARCC1, SMARCA4/BRG1 and ARID1A. Additionally, the histone deacetylase HDAC2, PRC1 complex component CBX2, chromatin-modifying enzyme KDM1A/LSD1 and pioneer factor c-MYB are regulated by SOX11. Finally, SOX11 is identified as a core transcription factor of the core regulatory circuitry (CRC) in adrenergic high-risk neuroblastoma with a potential role as epigenetic master regulator upstream of the CRC.

MEK inhibition causes BIM stabilization and increased sensitivity to BCL-2 family member inhibitors in RAS-MAPK-mutated neuroblastoma.

Introduction: Mutations affecting the RAS-MAPK pathway occur frequently in relapsed neuroblastoma tumors and are associated with response to MEK inhibition in vitro. However, these inhibitors alone do not lead to tumor regression in vivo, indicating the need for combination therapy. Methods and results: Via high-throughput combination screening, we identified that the MEK inhibitor trametinib can be combined with BCL-2 family member inhibitors, to efficiently inhibit growth of neuroblastoma cell lines with RAS-MAPK mutations. Suppressing the RAS-MAPK pathway with trametinib led to an increase in pro-apoptotic BIM, resulting in more BIM binding to anti-apoptotic BCL-2 family members. By favoring the formation of these complexes, trametinib treatment enhances sensitivity to compounds targeting anti-apoptotic BCL-2 family members. In vitro validation studies confirmed that this sensitizing effect is dependent on an active RAS-MAPK pathway. In vivo combination of trametinib with BCL-2 inhibitors led to tumor inhibition in NRAS-mutant and NF1-deleted xenografts. Conclusion: Together, these results show that combining MEK inhibition with BCL-2 family member inhibition could potentially improve therapeutic outcomes for RAS-MAPK-mutated neuroblastoma patients.

FOXO3A-short is a novel regulator of non-oxidative glucose metabolism associated with human longevity.

Intronic single-nucleotide polymorphisms (SNPs) in FOXO3A are associated with human longevity. Currently, it is unclear how these SNPs alter FOXO3A functionality and human physiology, thereby influencing lifespan. Here, we identify a primate-specific FOXO3A transcriptional isoform, FOXO3A-Short (FOXO3A-S), encoding a major longevity-associated SNP, rs9400239 (C or T), within its 5' untranslated region. The FOXO3A-S mRNA is highly expressed in the skeletal muscle and has very limited expression in other tissues. We find that the rs9400239 variant influences the stability and functionality of the primarily nuclear protein(s) encoded by the FOXO3A-S mRNA. Assessment of the relationship between the FOXO3A-S polymorphism and peripheral glucose clearance during insulin infusion (Rd clamp) in a cohort of Danish twins revealed that longevity T-allele carriers have markedly faster peripheral glucose clearance rates than normal lifespan C-allele carriers. In vitro experiments in human myotube cultures utilizing overexpression of each allele showed that the C-allele represses glycolysis independently of PI3K signaling, while overexpression of the T-allele represses glycolysis only in a PI3K-inactive background. Supporting this finding inducible knockdown of the FOXO3A-S C-allele in cultured myotubes increases the glycolytic rate. We conclude that the rs9400239 polymorphism acts as a molecular switch which changes the identity of the FOXO3A-S-derived protein(s), which in turn alters the relationship between FOXO3A-S and insulin/PI3K signaling and glycolytic flux in the skeletal muscle. This critical difference endows carriers of the FOXO3A-S T-allele with consistently higher insulin-stimulated peripheral glucose clearance rates, which may contribute to their longer and healthier lifespans.

Mesenchymal and adrenergic cell lineage states in neuroblastoma possess distinct immunogenic phenotypes.

Apart from the anti-GD2 antibody, immunotherapy for neuroblastoma has had limited success due to immune evasion mechanisms, coupled with an incomplete understanding of predictors of response. Here, from bulk and single-cell transcriptomic analyses, we identify a subset of neuroblastomas enriched for transcripts associated with immune activation and inhibition and show that these are predominantly characterized by gene expression signatures of the mesenchymal lineage state. By contrast, tumors expressing adrenergic lineage signatures are less immunogenic. The inherent presence or induction of the mesenchymal state through transcriptional reprogramming or therapy resistance is accompanied by innate and adaptive immune gene activation through epigenetic remodeling. Mesenchymal lineage cells promote T cell infiltration by secreting inflammatory cytokines, are efficiently targeted by cytotoxic T and natural killer cells and respond to immune checkpoint blockade. Together, we demonstrate that distinct immunogenic phenotypes define the divergent lineage states of neuroblastoma and highlight the immunogenic potential of the mesenchymal lineage.

DHODH is an independent prognostic marker and potent therapeutic target in neuroblastoma.

Despite intensive therapy, children with high-risk neuroblastoma are at risk of treatment failure. We applied a multiomic system approach to evaluate metabolic vulnerabilities in human neuroblastoma. We combined metabolomics, CRISPR screening, and transcriptomic data across more than 700 solid tumor cell lines and identified dihydroorotate dehydrogenase (DHODH), a critical enzyme in pyrimidine synthesis, as a potential treatment target. Of note, DHODH inhibition is currently under clinical investigation in patients with hematologic malignancies. In neuroblastoma, DHODH expression was identified as an independent risk factor for aggressive disease, and high DHODH levels correlated to worse overall and event-free survival. A subset of tumors with the highest DHODH expression was associated with a dismal prognosis, with a 5-year survival of less than 10%. In xenograft and transgenic neuroblastoma mouse models treated with the DHODH inhibitor brequinar, tumor growth was dramatically reduced, and survival was extended. Furthermore, brequinar treatment was shown to reduce the expression of MYC targets in 3 neuroblastoma models in vivo. A combination of brequinar and temozolomide was curative in the majority of transgenic TH-MYCN neuroblastoma mice, indicating a highly active clinical combination therapy. Overall, DHODH inhibition combined with temozolomide has therapeutic potential in neuroblastoma, and we propose this combination for clinical testing.

Mesenchymal-Type Neuroblastoma Cells Escape ALK Inhibitors.

Cancer therapy frequently fails due to the emergence of resistance. Many tumors include phenotypically immature tumor cells, which have been implicated in therapy resistance. Neuroblastoma cells can adopt a lineage-committed adrenergic (ADRN) or an immature mesenchymal (MES) state. They differ in epigenetic landscape and transcription factors, and MES cells are more resistant to chemotherapy. Here we analyzed the response of MES cells to targeted drugs. Activating anaplastic lymphoma kinase (ALK) mutations are frequently found in neuroblastoma and ALK inhibitors (ALKi) are in clinical trials. ALKi treatment of ADRN neuroblastoma cells with a tumor-driving ALK mutation induced cell death. Conversely, MES cells did not express either mutant or wild-type ALK and were resistant to ALKi, and MES cells formed tumors that progressed under ALKi therapy. In assessing the role of MES cells in relapse development, TRAIL was identified to specifically induce apoptosis in MES cells and to suppress MES tumor growth. Addition of TRAIL to ALKi treatment of neuroblastoma xenografts delayed relapses in a subset of the animals, suggesting a role for MES cells in relapse formation. While ADRN cells resembled normal embryonal neuroblasts, MES cells resembled immature precursor cells, which also lacked ALK expression. Resistance to targeted drugs can therefore be an intrinsic property of immature cancer cells based on their resemblance to developmental precursors. SIGNIFICANCE: In neuroblastoma, mesenchymal tumor cells lack expression of the tumor-driving ALK oncogene and are resistant to ALKi, but dual treatment with ALKi and mesenchymal cell-targeting TRAIL delays tumor relapse.

Defects in 8-oxo-guanine repair pathway cause high frequency of C > A substitutions in neuroblastoma.

Neuroblastomas are childhood tumors with frequent fatal relapses after induction treatment, which is related to tumor evolution with additional genomic events. Our whole-genome sequencing data analysis revealed a high frequency of somatic cytosine > adenine (C > A) substitutions in primary neuroblastoma tumors, which was associated with poor survival. We showed that increased levels of C > A substitutions correlate with copy number loss (CNL) of OGG1 or MUTYH Both genes encode DNA glycosylases that recognize 8-oxo-guanine (8-oxoG) lesions as a first step of 8-oxoG repair. Tumor organoid models with CNL of OGG1 or MUTYH show increased 8-oxoG levels compared to wild-type cells. We used CRISPR-Cas9 genome editing to create knockout clones of MUTYH and OGG1 in neuroblastoma cells. Whole-genome sequencing of single-cell OGG1 and MUTYH knockout clones identified an increased accumulation of C > A substitutions. Mutational signature analysis of these OGG1 and MUTYH knockout clones revealed enrichment for C > A signatures 18 and 36, respectively. Clustering analysis showed that the knockout clones group together with tumors containing OGG1 or MUTYH CNL. In conclusion, we demonstrate that defects in 8-oxoG repair cause accumulation of C > A substitutions in neuroblastoma, which contributes to mutagenesis and tumor evolution.

A G316A Polymorphism in the Ornithine Decarboxylase Gene Promoter Modulates MYCN-Driven Childhood Neuroblastoma.

Ornithine decarboxylase (ODC1), a critical regulatory enzyme in polyamine biosynthesis, is a direct transcriptional target of MYCN, amplification of which is a powerful marker of aggressive neuroblastoma. A single nucleotide polymorphism (SNP), G316A, within the first intron of ODC1, results in genotypes wildtype GG, and variants AG/AA. CRISPR-cas9 technology was used to investigate the effects of AG clones from wildtype MYCN-amplified SK-N-BE(2)-C cells and the effect of the SNP on MYCN binding, and promoter activity was investigated using EMSA and luciferase assays. AG clones exhibited decreased ODC1 expression, growth rates, and histone acetylation and increased sensitivity to ODC1 inhibition. MYCN was a stronger transcriptional regulator of the ODC1 promoter containing the G allele, and preferentially bound the G allele over the A. Two neuroblastoma cohorts were used to investigate the clinical impact of the SNP. In the study cohort, the minor AA genotype was associated with improved survival, while poor prognosis was associated with the GG genotype and AG/GG genotypes in MYCN-amplified and non-amplified patients, respectively. These effects were lost in the GWAS cohort. We have demonstrated that the ODC1 G316A polymorphism has functional significance in neuroblastoma and is subject to allele-specific regulation by the MYCN oncoprotein.

Super enhancers define regulatory subtypes and cell identity in neuroblastoma.

Half of the children diagnosed with neuroblastoma (NB) have high-risk disease, disproportionately contributing to overall childhood cancer-related deaths. In addition to recurrent gene mutations, there is increasing evidence supporting the role of epigenetic deregulation in disease pathogenesis. Yet, comprehensive cis-regulatory network descriptions from NB are lacking. Here, using genome-wide H3K27ac profiles across 60 NBs, covering the different clinical and molecular subtypes, we identified four major super-enhancer-driven epigenetic subtypes and their underlying master regulatory networks. Three of these subtypes recapitulated known clinical groups; namely, MYCN-amplified, MYCN non-amplified high-risk and MYCN non-amplified low-risk NBs. The fourth subtype, exhibiting mesenchymal characteristics, shared cellular identity with multipotent Schwann cell precursors, was induced by RAS activation and was enriched in relapsed disease. Notably, CCND1, an essential gene in NB, was regulated by both mesenchymal and adrenergic regulatory networks converging on distinct super-enhancer modules. Overall, this study reveals subtype-specific super-enhancer regulation in NBs.

A cancer drug atlas enables synergistic targeting of independent drug vulnerabilities.

Personalized cancer treatments using combinations of drugs with a synergistic effect is attractive but proves to be highly challenging. Here we present an approach to uncover the efficacy of drug combinations based on the analysis of mono-drug effects. For this we used dose-response data from pharmacogenomic encyclopedias and represent these as a drug atlas. The drug atlas represents the relations between drug effects and allows to identify independent processes for which the tumor might be particularly vulnerable when attacked by two drugs. Our approach enables the prediction of combination-therapy which can be linked to tumor-driving mutations. By using this strategy, we can uncover potential effective drug combinations on a pan-cancer scale. Predicted synergies are provided and have been validated in glioblastoma, breast cancer, melanoma and leukemia mouse-models, resulting in therapeutic synergy in 75% of the tested models. This indicates that we can accurately predict effective drug combinations with translational value.

A NOTCH feed-forward loop drives reprogramming from adrenergic to mesenchymal state in neuroblastoma.

Transition between differentiation states in development occurs swift but the mechanisms leading to epigenetic and transcriptional reprogramming are poorly understood. The pediatric cancer neuroblastoma includes adrenergic (ADRN) and mesenchymal (MES) tumor cell types, which differ in phenotype, super-enhancers (SEs) and core regulatory circuitries. These cell types can spontaneously interconvert, but the mechanism remains largely unknown. Here, we unravel how a NOTCH3 intracellular domain reprogrammed the ADRN transcriptional landscape towards a MES state. A transcriptional feed-forward circuitry of NOTCH-family transcription factors amplifies the NOTCH signaling levels, explaining the swift transition between two semi-stable cellular states. This transition induces genome-wide remodeling of the H3K27ac landscape and a switch from ADRN SEs to MES SEs. Once established, the NOTCH feed-forward loop maintains the induced MES state. In vivo reprogramming of ADRN cells shows that MES and ADRN cells are equally oncogenic. Our results elucidate a swift transdifferentiation between two semi-stable epigenetic cellular states.

Recruitment of BRCA1 limits MYCN-driven accumulation of stalled RNA polymerase.

MYC is an oncogenic transcription factor that binds globally to active promoters and promotes transcriptional elongation by RNA polymerase II (RNAPII)1,2. Deregulated expression of the paralogous protein MYCN drives the development of neuronal and neuroendocrine tumours and is often associated with a particularly poor prognosis3. Here we show that, similar to MYC, activation of MYCN in human neuroblastoma cells induces escape of RNAPII from promoters. If the release of RNAPII from transcriptional pause sites (pause release) fails, MYCN recruits BRCA1 to promoter-proximal regions. Recruitment of BRCA1 prevents MYCN-dependent accumulation of stalled RNAPII and enhances transcriptional activation by MYCN. Mechanistically, BRCA1 stabilizes mRNA decapping complexes and enables MYCN to suppress R-loop formation in promoter-proximal regions. Recruitment of BRCA1 requires the ubiquitin-specific protease USP11, which binds specifically to MYCN when MYCN is dephosphorylated at Thr58. USP11, BRCA1 and MYCN stabilize each other on chromatin, preventing proteasomal turnover of MYCN. Because BRCA1 is highly expressed in neuronal progenitor cells during early development4 and MYC is less efficient than MYCN in recruiting BRCA1, our findings indicate that a cell-lineage-specific stress response enables MYCN-driven tumours to cope with deregulated RNAPII function.

Mesenchymal Neuroblastoma Cells Are Undetected by Current mRNA Marker Panels: The Development of a Specific Neuroblastoma Mesenchymal Minimal Residual Disease Panel.

Patients with neuroblastoma in molecular remission remain at considerable risk for disease recurrence. Studies have found that neuroblastoma tissue contains adrenergic (ADRN) and mesenchymal (MES) cells; the latter express low levels of commonly used markers for minimal residual disease (MRD). We identified MES-specific MRD markers and studied the dynamics of these markers during treatment. PATIENTS AND METHODS: Microarray data were used to identify genes differentially expressed between ADRN and MES cell lines. Candidate genes were then studied using real-time quantitative polymerase chain reaction in cell lines and control bone marrow and peripheral blood samples. After selecting a panel of markers, serial bone marrow, peripheral blood, and peripheral blood stem cell samples were obtained from patients with high-risk neuroblastoma and tested for marker expression; survival analyses were also performed. RESULTS: PRRX1, POSTN, and FMO3 mRNAs were used as a panel for specifically detecting MES mRNA in patient samples. MES mRNA was detected only rarely in peripheral blood; moreover, the presence of MES mRNA in peripheral blood stem cell samples was associated with low event-free survival and overall survival. Of note, during treatment, serial bone marrow samples obtained from 29 patients revealed a difference in dynamics between MES mRNA markers and ADRN mRNA markers. Furthermore, MES mRNA was detected in a higher percentage of patients with recurrent disease than in those who remained disease free (53% v 32%, respectively; P = .03). CONCLUSION: We propose that the markers POSTN and PRRX1, in combination with FMO3, be used for real-time quantitative polymerase chain reaction-based detection of MES neuroblastoma mRNA in patient samples because these markers have a unique pattern during treatment and are more prevalent in patients with poor outcome. Together with existing markers of MRD, these new markers should be investigated further in large prospective studies.

ALK positively regulates MYCN activity through repression of HBP1 expression.

ALK mutations occur in 10% of primary neuroblastomas and represent a major target for precision treatment. In combination with MYCN amplification, ALK mutations infer an ultra-high-risk phenotype resulting in very poor patient prognosis. To open up opportunities for future precision drugging, a deeper understanding of the molecular consequences of constitutive ALK signaling and its relationship to MYCN activity in this aggressive pediatric tumor entity will be essential. We show that mutant ALK downregulates the 'HMG-box transcription factor 1' (HBP1) through the PI3K-AKT-FOXO3a signaling axis. HBP1 inhibits both the transcriptional activating and repressing activity of MYCN, the latter being mediated through PRC2 activity. HBP1 itself is under negative control of MYCN through miR-17~92. Combined targeting of HBP1 by PI3K antagonists and MYCN signaling by BET- or HDAC-inhibitors blocks MYCN activity and significantly reduces tumor growth, suggesting a novel targeted therapy option for high-risk neuroblastoma.

TBX2 is a neuroblastoma core regulatory circuitry component enhancing MYCN/FOXM1 reactivation of DREAM targets.

Chromosome 17q gains are almost invariably present in high-risk neuroblastoma cases. Here, we perform an integrative epigenomics search for dosage-sensitive transcription factors on 17q marked by H3K27ac defined super-enhancers and identify TBX2 as top candidate gene. We show that TBX2 is a constituent of the recently established core regulatory circuitry in neuroblastoma with features of a cell identity transcription factor, driving proliferation through activation of p21-DREAM repressed FOXM1 target genes. Combined MYCN/TBX2 knockdown enforces cell growth arrest suggesting that TBX2 enhances MYCN sustained activation of FOXM1 targets. Targeting transcriptional addiction by combined CDK7 and BET bromodomain inhibition shows synergistic effects on cell viability with strong repressive effects on CRC gene expression and p53 pathway response as well as several genes implicated in transcriptional regulation. In conclusion, we provide insight into the role of the TBX2 CRC gene in transcriptional dependency of neuroblastoma cells warranting clinical trials using BET and CDK7 inhibitors.

RAS-MAPK Pathway-Driven Tumor Progression Is Associated with Loss of CIC and Other Genomic Aberrations in Neuroblastoma.

Mutations affecting the RAS-MAPK pathway frequently occur in relapsed neuroblastoma tumors, which suggests that activation of this pathway is associated with a more aggressive phenotype. To explore this hypothesis, we generated several model systems to define a neuroblastoma RAS-MAPK pathway signature. Activation of this pathway in primary tumors indeed correlated with poor survival and was associated with known activating mutations in ALK and other RAS-MAPK pathway genes. Integrative analysis showed that mutations in PHOX2B, CIC, and DMD were also associated with an activated RAS-MAPK pathway. Mutation of PHOX2B and deletion of CIC in neuroblastoma cell lines induced activation of the RAS-MAPK pathway. This activation was independent of phosphorylated ERK in CIC knockout systems. Furthermore, deletion of CIC caused a significant increase in tumor growth in vivo These results show that the RAS-MAPK pathway is involved in tumor progression and establish CIC as a powerful tumor suppressor that functions downstream of this pathway in neuroblastoma.Significance: This work identifies CIC as a powerful tumor suppressor affecting the RAS-MAPK pathway in neuroblastoma and reinforces the importance of mutation-driven activation of this pathway in cancer. Cancer Res; 78(21); 6297-307. ©2018 AACR.

Genomic Amplifications and Distal 6q Loss: Novel Markers for Poor Survival in High-risk Neuroblastoma Patients.

Background: Neuroblastoma is characterized by substantial clinical heterogeneity. Despite intensive treatment, the survival rates of high-risk neuroblastoma patients are still disappointingly low. Somatic chromosomal copy number aberrations have been shown to be associated with patient outcome, particularly in low- and intermediate-risk neuroblastoma patients. To improve outcome prediction in high-risk neuroblastoma, we aimed to design a prognostic classification method based on copy number aberrations. Methods: In an international collaboration, normalized high-resolution DNA copy number data (arrayCGH and SNP arrays) from 556 high-risk neuroblastomas obtained at diagnosis were collected from nine collaborative groups and segmented using the same method. We applied logistic and Cox proportional hazard regression to identify genomic aberrations associated with poor outcome. Results: In this study, we identified two types of copy number aberrations that are associated with extremely poor outcome. Distal 6q losses were detected in 5.9% of patients and were associated with a 10-year survival probability of only 3.4% (95% confidence interval [CI] = 0.5% to 23.3%, two-sided P = .002). Amplifications of regions not encompassing the MYCN locus were detected in 18.1% of patients and were associated with a 10-year survival probability of only 5.8% (95% CI = 1.5% to 22.2%, two-sided P < .001). Conclusions: Using a unique large copy number data set of high-risk neuroblastoma cases, we identified a small subset of high-risk neuroblastoma patients with extremely low survival probability that might be eligible for inclusion in clinical trials of new therapeutics. The amplicons may also nominate alternative treatments that target the amplified genes.

Cross-Cohort Analysis Identifies a TEAD4-MYCN Positive Feedback Loop as the Core Regulatory Element of High-Risk Neuroblastoma.

High-risk neuroblastomas show a paucity of recurrent somatic mutations at diagnosis. As a result, the molecular basis for this aggressive phenotype remains elusive. Recent progress in regulatory network analysis helped us elucidate disease-driving mechanisms downstream of genomic alterations, including recurrent chromosomal alterations. Our analysis identified three molecular subtypes of high-risk neuroblastomas, consistent with chromosomal alterations, and identified subtype-specific master regulator proteins that were conserved across independent cohorts. A 10-protein transcriptional module-centered around a TEAD4-MYCN positive feedback loop-emerged as the regulatory driver of the high-risk subtype associated with MYCN amplification. Silencing of either gene collapsed MYCN-amplified (MYCNAmp) neuroblastoma transcriptional hallmarks and abrogated viability in vitro and in vivo Consistently, TEAD4 emerged as a robust prognostic marker of poor survival, with activity independent of the canonical Hippo pathway transcriptional coactivators YAP and TAZ. These results suggest novel therapeutic strategies for the large subset of MYCN-deregulated neuroblastomas.Significance: Despite progress in understanding of neuroblastoma genetics, little progress has been made toward personalized treatment. Here, we present a framework to determine the downstream effectors of the genetic alterations sustaining neuroblastoma subtypes, which can be easily extended to other tumor types. We show the critical effect of disrupting a 10-protein module centered around a YAP/TAZ-independent TEAD4-MYCN positive feedback loop in MYCNAmp neuroblastomas, nominating TEAD4 as a novel candidate for therapeutic intervention. Cancer Discov; 8(5); 582-99. ©2018 AACR.This article is highlighted in the In This Issue feature, p. 517.

Indian Hedgehog Suppresses a Stromal Cell-Driven Intestinal Immune Response.

BACKGROUND & AIMS: Upon intestinal epithelial damage a complex wound healing response is initiated to restore epithelial integrity and defend against pathogenic invasion. Epithelium-derived Indian Hedgehog (Ihh) functions as a critical sensor in this process. Signaling occurs in a paracrine manner because the receptor for Ihh is expressed only in the mesenchyme, but the exact Hedgehog target cell has remained elusive. The aim of this study was to elucidate further the nature of this target cell in the context of intestinal inflammation. METHODS: Hedgehog activity was modulated genetically in both cell type-specific and body-wide models and the resulting animals were analyzed for gene expression profiles and sensitivity for dextran sodium sulfate (DSS) colitis. To characterize the Hedgehog target cell, Gli1-CreERT2-Rosa26-ZsGreen animals were generated, which express ZsGreen in all Hedgehog-responsive cells. These cells were characterized using flow cytometry and immunofluorescence. RESULTS: Loss of Indian Hedgehog from the intestinal epithelium resulted in a rapid increase in expression of inflammation-related genes, accompanied by increased influx of immune cells. Animals with epithelium-specific deletion of Ihh or lacking the Hedgehog receptor Smoothened from Hedgehog target cells were more sensitive to DSS colitis. In contrast, specific deletion of Smoothened in the myeloid compartment did not alter the response to DSS. This suggests that Hedgehog signaling does not repress intestinal immunity through an effect on myeloid cells. Indeed, we found that Hedgehog-responsive cells expressed gp38, smooth muscle actin, and desmin, indicating a fibroblastic nature. Ihh signaling inhibited expression of C-X-C motif chemokine ligand 12 (CXCL12) in fibroblasts in vitro and in vivo, thereby impairing the recruitment of immune cells. CONCLUSIONS: We show that epithelium-derived Indian Hedgehog signals exclusively to fibroblasts in the intestine. Loss of Ihh leads to a rapid immune response with up-regulation of fibroblast-derived CXCL12, and migration of immune cells into the lamina propria.