Background and Aims We previously identified small molecules predicted to reverse an ileal gene signature for future Crohn’s Disease (CD) strictures. Here we used a new human intestinal organoid (HIO) model system containing macrophages to test a lead candidate, eicosatetraynoic acid (ETYA). Methods Induced pluripotent stem cell lines (iPSC) were derived from CD patients and differentiated into macrophages and HIOs. Macrophages and macrophage:HIO co-cultures were exposed to lipopolysaccharide (LPS) with and without ETYA pre-treatment. Cytospin and flow cytometry characterized macrophage morphology and activation markers, and RNA sequencing defined the global pattern of macrophage gene expression. TaqMan Low Density Array, Luminex multiplex assay, immunohistologic staining, and sirius red polarized light microscopy were performed to measure macrophage cytokine production and HIO pro-fibrotic gene expression and collagen content. Results iPSC-derived macrophages exhibited morphology similar to primary macrophages and expressed inflammatory macrophage cell surface markers including CD64 and CD68. LPS-stimulated macrophages expressed a global pattern of gene expression enriched in CD ileal inflammatory macrophages and matrisome secreted products, and produced cytokines and chemokines including CCL2, IL1B, and OSM implicated in refractory disease. ETYA suppressed CD64 abundance and pro-fibrotic gene expression pathways in LPS stimulated macrophages. Co-culture of LPS-primed macrophages with HIO led to up-regulation of fibroblast activation genes including ACTA2 and COL1A1 , and an increase in HIO collagen content. ETYA pre-treatment prevented pro-fibrotic effects of LPS-primed macrophages. Conclusions ETYA inhibits pro-fibrotic effects of LPS-primed macrophages upon co-cultured HIO. This model may be used in future untargeted screens for small molecules to treat refractory CD.

Nutrient deprivation induces a reversible cell cycle arrest state termed quiescence, which often accompanies transcriptional silencing and chromatin compaction. Paradoxically, nutrient deprivation is associated with activated fibroblast states in pathological microenvironments in which fibroblasts drive extracellular matrix (ECM) remodeling to alter tissue environments. The relationship between nutrient deprivation and fibroblast activation remains unclear. Here, we report that serum deprivation extensively activates transcription of ECM remodeling genes in cultured fibroblasts, despite the induction of quiescence. Starvation-induced transcriptional activation accompanied large-scale histone acetylation of putative distal enhancers, but not promoters. The starvation-activated putative enhancers were enriched for non-coding genetic risk variants associated with inflammatory bowel disease (IBD), suggesting that the starvation-activated gene regulatory network may contribute to fibroblast activation in IBD. Indeed, the starvation-activated gene PLAU, encoding uPA serine protease for plasminogen and ECM, was upregulated in inflammatory fibroblasts in the intestines of IBD patients. Furthermore, the starvation-activated putative enhancer at PLAU, which harbors an IBD risk variant, gained chromatin accessibility in IBD patient fibroblasts. This study implicates nutrient deprivation in transcriptional activation of ECM remodeling genes in fibroblasts and suggests nutrient deprivation as a potential mechanism for pathological fibroblast activation in IBD.

Inflammatory Bowel Disease ( IBD ) is a chronic and often debilitating autoinflammatory condition, with an increasing incidence in children. Standard-of-care therapies lead to sustained transmural healing and clinical remission in fewer than one-third of patients. For children, TNFα inhibition remains the only FDA-approved biologic therapy, providing an even greater urgency to understanding mechanisms of response. Genome-wide association studies ( GWAS ) have identified 418 independent genetic risk loci contributing to IBD, yet the majority are noncoding and their mechanisms of action are difficult to decipher. If causal, they likely alter transcription factor ( TF ) binding and downstream gene expression in particular cell types and contexts. To bridge this knowledge gap, we built a novel resource: multiome-seq (tandem single-nuclei ( sn )RNA-seq and chromatin accessibility ( snATAC )-seq) of intestinal tissue from pediatric IBD patients, where anti-TNF response was defined by endoscopic healing. From the snATAC-seq data, we generated a first-time atlas of chromatin accessibility (putative regulatory elements) for diverse intestinal cell types in the context of IBD. For cell types/contexts mediating genetic risk, we reasoned that accessible chromatin will co-localize with genetic disease risk loci. We systematically tested for significant co-localization of our chromatin accessibility maps and risk variants for 758 GWAS traits. Globally, genetic risk variants for IBD, autoimmune and inflammatory diseases are enriched in accessible chromatin of immune populations, while other traits (e.g., colorectal cancer, metabolic) are enriched in epithelial and stromal populations. This resource opens new avenues to uncover the complex molecular and cellular mechanisms mediating genetic disease risk.

Abstract BACKGROUND Mutations in the DUOX2 intestinal epithelial cell NADPH oxidase are associated with risk for Crohn’s Disease (CD), and may influence wound healing and the development of strictures. Patients who develop strictures exhibit reduced expression of genes encoding mitochondrial sub-units, and increased expression of genes encoding extra-cellular matrix proteins, in the ileum at diagnosis. We conducted a perturbagen bioinformatics analysis which predicted that the cyclooxygenase and lipoxygenase inhibitor eicosatetraynoic-acid (ETYA) would reverse the ileal gene signature associated with stricture development. In the current study we tested effects of ETYA in a human intestinal organoid model (HIO) system. METHODS HIO were derived from wild type (WT) and DUOX2 mutant (DUOX2var) induced pluripotent stem cells (iPSCs) prepared from CD patients. WT and DUOX2var HIO were examined under basal conditions, and following exposure to ETYA for 72 hours or 12 days. Superoxide production in EPCAM+ epithelial cells and CD90+ stromal cells was measured by flow cytometry. RNA was prepared and expression of mitochondrial and extra-cellular matrix genes linked to strictures in patients was determined using RNA sequencing and a TaqMan Low Density Array (TLDA) card. HIO tissue stiffness was measured using atomic force microscopy (AFM). RESULTS ETYA reduced superoxide production by HIO EPCAM+ epithelial cells only in WT organoids; no effect was observed in DUOX2var HIO. This was specific, as no effect of ETYA upon superoxide production was detected in CD90+ stromal cells. Under basal conditions, RNA sequencing demonstrated that WT HIO expressed core mitochondrial and extra-cellular matrix genes and enriched biologic functions implicated in CD strictures. Expression of extra-cellular matrix and wound healing genes was increased in DUOX2var HIO under basal conditions. ETYA up-regulated expression of the COX5B, POLG2, and SLC25A27 mitochondrial genes associated with lower rates of strictures, while reducing expression of the ACTA2, VIM, and COL1A1 extra-cellular matrix genes associated with higher rates of strictures, independent of DUOX2 genotype. WT and DUOX2var HIO exhibited tissue stiffness comparable to normal human ileum under basal conditions; this was significantly reduced by ETYA exposure only in DUOX2var HIO. CONCLUSIONS ETYA regulates mitochondrial and extra-cellular matrix genes implicated in stricture formation in an HIO model system. DUOX2var HIO exhibit increased extra-cellular matrix gene expression under basal conditions, which is reduced by ETYA exposure in conjunction with a reduction in tissue stiffness. Collectively, these data confirm that HIO provide a relevant model system to study mechanisms regulating stricture formation in CD, including screening of small molecules prioritized by perturbagen bioinformatics analysis.

Background Perturbagen analysis of Crohn’s disease (CD) ileal gene expression data identified small molecules including eicosatetraynoic acid (ETYA), which may exert an antifibrotic effect. We developed a patient-specific human intestinal organoid (HIO) model system to test small molecule regulation of mitochondrial and wound-healing functions implicated in stricturing behavior. Methods HIOs were made from CD induced pluripotent stem cells with and without a loss-of-function haplotype in the DUOX2 gene implicated in ileal homeostasis and characterized under basal conditions and following exposure to butyrate and ETYA using RNA sequencing, flow cytometry, and immunofluorescent and polarized light microscopy. Mitochondrial activity was measured using high-resolution respirometry and tissue stiffness using atomic force microscopy. Results HIOs expressed core mitochondrial and extracellular matrix (ECM) genes and enriched biologic functions implicated in CD ileal strictures; ECM gene expression was suppressed by both butyrate and ETYA, with butyrate also suppressing genes regulating epithelial proliferation. Consistent with this, butyrate, but not ETYA, exerted a profound effect on HIO epithelial mitochondrial function, reactive oxygen species production, and cellular abundance. Butyrate and ETYA suppressed HIO expression of alpha smooth muscle actin expressed by myofibroblasts, type I collagen, and collagen protein abundance. HIOs exhibited tissue stiffness comparable to normal human ileum; this was reduced by chronic ETYA exposure in HIOs carrying the DUOX2 loss-of-function haplotype. Conclusions ETYA regulates ECM genes implicated in strictures and suppresses collagen content and tissue stiffness in an HIO model. HIOs provide a platform to test personalized therapeutics, including small molecules prioritized by perturbagen analysis.

Abstract Introduction The DUOX2 intestinal epithelial NADPH oxidase is upregulated in Crohn’s Disease (CD), and DUOX2 mutations are associated with increased CD risk. Oxidative stress and loss of mitochondrial function disrupt the intestinal barrier promoting inflammatory responses to commensals. The relative impact of DUOX2 mutations and microbial products in this regard is poorly understood. Hypothesis We hypothesized that DUOX2 genetic variation would be associated with differences in cellular reactive oxygen species (ROS) production and mitochondrial function in a Human Intestinal Organoid (HIO) model system. Methods Induced pluripotent stem cell lines derived from pediatric CD patients with and without combined DUOX2 missense mutations(R701Q, P982A, and H678R) were used to generate wild type (WT) and DUOX2mut HIOs. Reactive oxygen species (ROS) production was measured using the two-color ROS-ID® Total ROS/Superoxide detection kit, and the mitochondrial membrane potential (MMP) was measured using JC1 staining by flow cytometry in HIO EpCAM+ epithelial cells and CD90+ stromal cells. Expression of inflammatory and mitochondrial genes which varied with DUOX2 mutation carriage in CD patent ileal biopsies was measured by RT-PCR. HIO mitochondrial complex I and II activity was measured using an Oroboros respirometer. Results Epithelial ROS production was reduced in DUOX2mut HIO under basal conditions; this difference was not observed following pyocyanin stimulation (Fig. 1A). A profound suppression of epithelial ROS production was observed following butyrate treatment. Butyrate did not alter stromal cell ROS production. Under these conditions, induction of ROS by pyocyanin was abrogated in WT, but not DUOX2mut HIO epithelial cells (Fig. 1B). Butyrate increased expression of core genes regulating the mitochondrial respiratory chain and DNA synthesis (COX5B, NDUFA1, POLG2, SLC25A27) and HIF1A implicated in barrier function, independent of genotype (p<0.05). The epithelial and stromal cell mitochondrial membrane potential (MMP) (Fig. 2A), and HIO mitochondrial complex I activity (Fig. 2B), were reduced in DUOX2mut HIO under basal conditions. This was specific, as mitochondrial complex II activity did not vary with DUOX2 genotype. Conclusions We confirmed epithelial effects of DUOX2 genotype and butyrate exposure on ROS production in the HIO model system. Genotype dependent effects on basal ROS production were largely abrogated by the microbial products pyocyanin and butyrate, although butyrate inhibition of pyocyanin induced ROS production was dependent on intact DUOX2 function. Data suggest a previously unanticipated effect of DUOX2 genetic variation on the epithelial and stromal cell MMP and cellular respiration. This may have implications for mechanisms by which DUOX2 regulates barrier function and inflammatory responses to commensals in CD.

Cell polarization occurs along a single axis that is generally determined by a spatial cue. Cells of the budding yeast Saccharomyces cerevisiae select a site for polarized growth in a specific pattern depending on cell type. Haploid a and α cells bud in the axial budding pattern, which depends on a transient marker and requires proteins Bud3, Bud4, Axl1 and Axl2. Here we report that Bud4 functions as a platform that mediates the ordered assembly of the axial landmark at the division site during M and early G1. While Bud4 associates with Bud3 in all cell types and in the absence of Axl1 or Axl2, Bud4 interacts with Axl1 and Axl2 mainly in haploid cells and only in the presence of all other components of the landmark. Bud4 can bind to GTP or GDP, and a GTP binding-defective Bud4 fails to interact with Axl1 in vitro. The same bud4 mutation leads to mis-localization of Axl1 and disrupts the axial budding pattern, indicating that GTP binding to Bud4 is important for its role in bud-site selection. We also show the cell-type-specific association of the axial landmark with Bud5, a GDP/GTP exchange factor for Rsr1. Despite their expression in all cell types, Bud4 and Axl2 associate with Bud5 specifically in haploid cells and in the presence of Axl1, whose expression is limited to a and α cells. Together, our findings suggest that Bud4 plays a critical role in the assembly of the axial landmark and its link to the Rsr1 GTPase module.

In the budding yeast Saccharomyces cerevisiae, selection of the bud site determines the axis of polarized cell growth and eventual oriented cell division. Bud sites are selected in specific patterns depending on cell type. These patterns appear to depend on distinct types of marker proteins in the cell cortex; in particular, the bipolar budding of diploid cells depends on persistent landmarks at the birth-scar-distal and -proximal poles that involve the proteins Bud8p and Bud9p, respectively. Rax1p and Rax2p also appear to function specifically in bipolar budding, and we report here a further characterization of these proteins and of their interactions with Bud8p and Bud9p. Rax1p and Rax2p both appear to be integral membrane proteins. Although commonly used programs predict different topologies for Rax2p, glycosylation studies indicate that it has a type I orientation, with its long N-terminal domain in the extracytoplasmic space. Analysis of rax1 and rax2 mutant budding patterns indicates that both proteins are involved in selecting bud sites at both the distal and proximal poles of daughter cells as well as near previously used division sites on mother cells. Consistent with this, GFP-tagged Rax1p and Rax2p were both observed at the distal pole as well as at the division site on both mother and daughter cells; localization to the division sites was persistent through multiple cell cycles. Localization of Rax1p and Rax2p was interdependent, and biochemical studies showed that these proteins could be copurified from yeast. Bud8p and Bud9p could also be copurified with Rax1p, and localization studies provided further evidence of interactions. Localization of Rax1p and Rax2p to the bud tip and distal pole depended on Bud8p, and normal localization of Bud8p was partially dependent on Rax1p and Rax2p. Although localization of Rax1p and Rax2p to the division site did not appear to depend on Bud9p, normal localization of Bud9p appeared largely or entirely dependent on Rax1p and Rax2p. Taken together, the results indicate that Rax1p and Rax2p interact closely with each other and with Bud8p and Bud9p in the establishment and/or maintenance of the cortical landmarks for bipolar budding.

Polarized cell growth requires the coupling of a defined spatial site on the cell cortex to the apparatus that directs the establishment of cell polarity. In the budding yeast Saccharomyces cerevisiae, the Ras-family GTPase Rsr1p/Bud1p and its regulators select the proper site for bud emergence on the cell cortex. The Rho-family GTPase Cdc42p and its associated proteins then establish an axis of polarized growth by triggering an asymmetric organization of the actin cytoskeleton and secretory apparatus at the selected bud site. We explored whether a direct linkage exists between the Rsr1p/Bud1p and Cdc42p GTPases. Here we show specific genetic interactions between RSR1/BUD1 and particular cdc42 mutants defective in polarity establishment. We also show that Cdc42p coimmunoprecipitated with Rsr1p/Bud1p from yeast extracts. In vitro studies indicated a direct interaction between Rsr1p/Bud1p and Cdc42p, which was enhanced by Cdc24p, a guanine nucleotide exchange factor for Cdc42p. Our findings suggest that Cdc42p interacts directly with Rsr1p/Bud1p in vivo, providing a novel mechanism by which direct contact between a Ras-family GTPase and a Rho-family GTPase links the selection of a growth site to polarity establishment.