Multimodal single-cell analyses reveal mechanisms of perianal fistula in diverse patients with Crohn's disease.

BACKGROUND: Crohn's disease complicated by perianal fistulae is more prevalent and severe in patients of African ancestry. METHODS: We profiled single cells from diverse patients with Crohn's disease with perianal fistula from colorectal mucosa and fistulous tracts. Immunofluorescence was performed to validate predicted cell-cell interactions. Unstimulated monocytes were chronically cultured in diverse cohorts. A subset was analyzed by single-nucleus RNA + ATAC sequencing. FINDINGS: Fistulous tract cells from complete proctectomies demonstrated enrichment of myeloid cells compared to paired rectal tissues. Ligand-receptor analysis highlights myeloid-stromal cross-talk and cellular senescence, with cellular co-localization validated by immunofluorescence. Chitinase-3 like-protein-1 (CHI3L1) is a top upregulated gene in stromal cells from fistulae expressing both destructive and fibrotic gene signatures. Monocyte cultures from patients of African ancestry and controls demonstrated differences in CHI3L1 and oncostatin M (OSM) expression upon differentiation compared to individuals of European ancestry. Activating protein-1 footprints are present in ATAC-seq peaks in stress response genes, including CHI3L1 and OSM; genome-wide chromatin accessibility including JUN footprints was observed, consistent with reported mechanisms of inflammatory memory. Regulon analyses confirm known cell-specific transcription factor regulation and implicate novel ones in fibroblast subsets. All pseudo-bulked clusters demonstrate enrichment of genetic loci, establishing multicellular contributions. In the most significant African American Crohn's genetic locus, upstream of prostaglandin E receptor 4, lymphoid-predominant ATAC-seq peaks were observed, with predicted RORC footprints. CONCLUSIONS: Population differences in myeloid-stromal cross-talk implicate fibrotic and destructive fibroblasts, senescence, epigenetic memory, and cell-specific enhancers in perianal fistula pathogenesis. The transcriptomic and epigenetic data provided here may guide optimization of promising mesenchymal stem cell therapies for perianal fistula. FUNDING: This work was supported by grants U01DK062422, U24DK062429, and R01DK123758.

Identifying high-impact variants and genes in exomes of Ashkenazi Jewish inflammatory bowel disease patients.

Inflammatory bowel disease (IBD) is a group of chronic digestive tract inflammatory conditions whose genetic etiology is still poorly understood. The incidence of IBD is particularly high among Ashkenazi Jews. Here, we identify 8 novel and plausible IBD-causing genes from the exomes of 4453 genetically identified Ashkenazi Jewish IBD cases (1734) and controls (2719). Various biological pathway analyses are performed, along with bulk and single-cell RNA sequencing, to demonstrate the likely physiological relatedness of the novel genes to IBD. Importantly, we demonstrate that the rare and high impact genetic architecture of Ashkenazi Jewish adult IBD displays significant overlap with very early onset-IBD genetics. Moreover, by performing biobank phenome-wide analyses, we find that IBD genes have pleiotropic effects that involve other immune responses. Finally, we show that polygenic risk score analyses based on genome-wide high impact variants have high power to predict IBD susceptibility.

NOX1 is essential for TNFα-induced intestinal epithelial ROS secretion and inhibits M cell signatures.

OBJECTIVE: Loss-of-function mutations in genes generating reactive oxygen species (ROS), such as NOX1, are associated with IBD. Mechanisms whereby loss of ROS drive IBD are incompletely defined. DESIGN: ROS measurements and single-cell transcriptomics were performed on colonoids stratified by NOX1 genotype and TNFα stimulation. Clustering of epithelial cells from human UC (inflamed and uninflamed) scRNASeq was performed. Validation of M cell induction was performed by immunohistochemistry using UEA1 (ulex europaeus agglutin-1 lectin) and in vivo with DSS injury. RESULTS: TNFα induces ROS production more in NOX1-WT versus NOX1-deficient murine colonoids under a range of Wnt-mediated and Notch-mediated conditions. scRNASeq from inflamed and uninflamed human colitis versus TNFα stimulated, in vitro colonoids defines substantially shared, induced transcription factors; NOX1-deficient colonoids express substantially lower levels of STAT3 (signal transducer and activator of transcription 3), CEBPD (CCAAT enhancer-binding protein delta), DNMT1 (DNA methyltransferase) and HIF1A (hypoxia-inducible factor) baseline. Subclustering unexpectedly showed marked TNFα-mediated induction of M cells (sentinel cells overlying lymphoid aggregates) in NOX1-deficient colonoids. M cell induction by UEA1 staining is rescued with H2O2 and paraquat, defining extra- and intracellular ROS roles in maintenance of LGR5+ stem cells. DSS injury demonstrated GP2 (glycoprotein-2), basal lymphoplasmacytosis and UEA1 induction in NOX1-deficiency. Principal components analyses of M cell genes and decreased DNMT1 RNA velocity correlate with UC inflammation. CONCLUSIONS: NOX1 deficiency plus TNFα stimulation contribute to colitis through dysregulation of the stem cell niche and altered cell differentiation, enhancing basal lymphoplasmacytosis. Our findings prioritise ROS modulation for future therapies.

Large-scale sequencing identifies multiple genes and rare variants associated with Crohn's disease susceptibility.

Genome-wide association studies (GWASs) have identified hundreds of loci associated with Crohn's disease (CD). However, as with all complex diseases, robust identification of the genes dysregulated by noncoding variants typically driving GWAS discoveries has been challenging. Here, to complement GWASs and better define actionable biological targets, we analyzed sequence data from more than 30,000 patients with CD and 80,000 population controls. We directly implicate ten genes in general onset CD for the first time to our knowledge via association to coding variation, four of which lie within established CD GWAS loci. In nine instances, a single coding variant is significantly associated, and in the tenth, ATG4C, we see additionally a significantly increased burden of very rare coding variants in CD cases. In addition to reiterating the central role of innate and adaptive immune cells as well as autophagy in CD pathogenesis, these newly associated genes highlight the emerging role of mesenchymal cells in the development and maintenance of intestinal inflammation.

Utility of polygenic embryo screening for disease depends on the selection strategy.

Polygenic risk scores (PRSs) have been offered since 2019 to screen in vitro fertilization embryos for genetic liability to adult diseases, despite a lack of comprehensive modeling of expected outcomes. Here we predict, based on the liability threshold model, the expected reduction in complex disease risk following polygenic embryo screening for a single disease. A strong determinant of the potential utility of such screening is the selection strategy, a factor that has not been previously studied. When only embryos with a very high PRS are excluded, the achieved risk reduction is minimal. In contrast, selecting the embryo with the lowest PRS can lead to substantial relative risk reductions, given a sufficient number of viable embryos. We systematically examine the impact of several factors on the utility of screening, including: variance explained by the PRS, number of embryos, disease prevalence, parental PRSs, and parental disease status. We consider both relative and absolute risk reductions, as well as population-averaged and per-couple risk reductions, and also examine the risk of pleiotropic effects. Finally, we confirm our theoretical predictions by simulating 'virtual' couples and offspring based on real genomes from schizophrenia and Crohn's disease case-control studies. We discuss the assumptions and limitations of our model, as well as the potential emerging ethical concerns.

Stratification of risk of progression to colectomy in ulcerative colitis via measured and predicted gene expression.

An important goal of clinical genomics is to be able to estimate the risk of adverse disease outcomes. Between 5% and 10% of individuals with ulcerative colitis (UC) require colectomy within 5 years of diagnosis, but polygenic risk scores (PRSs) utilizing findings from genome-wide association studies (GWASs) are unable to provide meaningful prediction of this adverse status. By contrast, in Crohn disease, gene expression profiling of GWAS-significant genes does provide some stratification of risk of progression to complicated disease in the form of a transcriptional risk score (TRS). Here, we demonstrate that a measured TRS based on bulk rectal gene expression in the PROTECT inception cohort study has a positive predictive value approaching 50% for colectomy. Single-cell profiling demonstrates that the genes are active in multiple diverse cell types from both the epithelial and immune compartments. Expression quantitative trait locus (QTL) analysis identifies genes with differential effects at baseline and week 52 follow-up, but for the most part, differential expression associated with colectomy risk is independent of local genetic regulation. Nevertheless, a predicted polygenic transcriptional risk score (PPTRS) derived by summation of transcriptome-wide association study (TWAS) effects identifies UC-affected individuals at 5-fold elevated risk of colectomy with data from the UK Biobank population cohort studies, independently replicated in an NIDDK-IBDGC dataset. Prediction of gene expression from relatively small transcriptome datasets can thus be used in conjunction with TWASs for stratification of risk of disease complications.

A myeloid-stromal niche and gp130 rescue in NOD2-driven Crohn's disease.

Crohn's disease is a chronic inflammatory intestinal disease that is frequently accompanied by aberrant healing and stricturing complications. Crosstalk between activated myeloid and stromal cells is critical in the pathogenicity of Crohn's disease1,2, and increases in intravasating monocytes are correlated with a lack of response to anti-TNF treatment3. The risk alleles with the highest effect on Crohn's disease are loss-of-function mutations in NOD24,5, which increase the risk of stricturing6. However, the mechanisms that underlie pathogenicity driven by NOD2 mutations and the pathways that might rescue a lack of response to anti-TNF treatment remain largely uncharacterized. Here we use direct ex vivo analyses of patients who carry risk alleles of NOD2 to show that loss of NOD2 leads to dysregulated homeostasis of activated fibroblasts and macrophages. CD14+ peripheral blood mononuclear cells from carriers of NOD2 risk alleles produce cells that express high levels of collagen, and elevation of conserved signatures is observed in nod2-deficient zebrafish models of intestinal injury. The enrichment of STAT3 regulation and gp130 ligands in activated fibroblasts and macrophages suggested that gp130 blockade might rescue the activated program in NOD2-deficient cells. We show that post-treatment induction of the STAT3 pathway is correlated with a lack of response to anti-TNF treatment in patients, and demonstrate in vivo in zebrafish the amelioration of the activated myeloid-stromal niche using the specific gp130 inhibitor bazedoxifene. Our results provide insights into NOD2-driven fibrosis in Crohn's disease, and suggest that gp130 blockade may benefit some patients with Crohn's disease-potentially as a complement to anti-TNF therapy.

Inflamed Ulcerative Colitis Regions Associated With MRGPRX2-Mediated Mast Cell Degranulation and Cell Activation Modules, Defining a New Therapeutic Target.

BACKGROUND & AIMS: Recent literature has implicated a key role for mast cells in murine models of colonic inflammation, but their role in human ulcerative colitis (UC) is not well established. A major advance has been the identification of mrgprb2 (human orthologue, MRGPX2) as mediating IgE-independent mast cell activation. We sought to define mechanisms of mast cell activation and MRGPRX2 in human UC. METHODS: Colon tissues were collected from patients with UC for bulk RNA sequencing and lamina propria cells were isolated for MRGPRX2 activation studies and single-cell RNA sequencing. Genetic association of all protein-altering G-protein coupled receptor single-nucleotide polymorphism was performed in an Ashkenazi Jewish UC case-control cohort. Variants of MRGPRX2 were transfected into Chinese hamster ovary (CHO) and human mast cell (HMC) 1.1 cells to detect genotype-dependent effects on ß-arrestin recruitment, IP-1 accumulation, and phosphorylated extracellular signal-regulated kinase. RESULTS: Mast cell-specific mediators and adrenomedullin (proteolytic precursor of PAMP-12, an MRGPRX2 agonist) are up-regulated in inflamed compared to uninflamed UC. MRGPRX2 stimulation induces carboxypeptidase secretion from inflamed UC. Of all protein-altering GPCR alleles, a unique variant of MRGPRX2, Asn62Ser, was most associated with and was bioinformatically predicted to alter arrestin recruitment. We validated that the UC protective serine allele enhances ß-arrestin recruitment, decreases IP-1, and increases phosphorylated extracellular signal-regulated kinase with MRGPRX2 agonists. Single-cell RNA sequencing defines that adrenomedullin is expressed by activated fibroblasts and epithelial cells and that interferon gamma is a key upstream regulator of mast cell gene expression. CONCLUSION: Inflamed UC regions are distinguished by MRGPRX2-mediated activation of mast cells, with decreased activation observed with a UC-protective genetic variant. These results define cell modules of UC activation and a new therapeutic target.

Common and Rare Variant Prediction and Penetrance of IBD in a Large, Multi-ethnic, Health System-based Biobank Cohort.

BACKGROUND AND AIMS: Polygenic risk scores (PRS) may soon be used to predict inflammatory bowel disease (IBD) risk in prevention efforts. We leveraged exome-sequence and single nucleotide polymorphism (SNP) array data from 29,358 individuals in the multiethnic, randomly ascertained health system-based BioMe biobank to define effects of common and rare IBD variants on disease prediction and pathophysiology. METHODS: PRS were calculated from European, African American, and Ashkenazi Jewish (AJ) reference case-control studies, and a meta-GWAS run using all three association datasets. PRS were then combined using regression to assess which combination of scores best predicted IBD status in European, AJ, Hispanic, and African American cohorts in BioMe. Additionally, rare variants were assessed in genes associated with very early-onset IBD (VEO-IBD), by estimating genetic penetrance in each BioMe population. RESULTS: Combining risk scores based on association data from distinct ancestral populations improved IBD prediction for every population in BioMe and significantly improved prediction among European ancestry UK Biobank individuals. Lower predictive power for non-Europeans was observed, reflecting in part substantially lower African IBD case-control reference sizes. We replicated associations for two VEO-IBD genes, ADAM17 and LRBA, with high dominant model penetrance in BioMe. Autosomal recessive LRBA risk alleles are associated with severe, early-onset autoimmunity; we show that heterozygous carriage of an African-predominant LRBA protein-altering allele is associated with significantly decreased LRBA and CTLA-4 expression with T-cell activation. CONCLUSIONS: Greater genetic diversity in African populations improves prediction across populations, and generalizes some VEO-IBD genes. Increasing African American IBD case-collections should be prioritized to reduce health disparities and enhance pathophysiological insight.

Single-Cell Analysis of Crohn's Disease Lesions Identifies a Pathogenic Cellular Module Associated with Resistance to Anti-TNF Therapy.

Clinical benefits of cytokine blockade in ileal Crohn's disease (iCD) are limited to a subset of patients. Here, we applied single-cell technologies to iCD lesions to address whether cellular heterogeneity contributes to treatment resistance. We found that a subset of patients expressed a unique cellular module in inflamed tissues that consisted of IgG plasma cells, inflammatory mononuclear phagocytes, activated T cells, and stromal cells, which we named the GIMATS module. Analysis of ligand-receptor interaction pairs identified a distinct network connectivity that likely drives the GIMATS module. Strikingly, the GIMATS module was also present in a subset of patients in four independent iCD cohorts (n = 441), and its presence at diagnosis correlated with failure to achieve durable corticosteroid-free remission upon anti-TNF therapy. These results emphasize the limitations of current diagnostic assays and the potential for single-cell mapping tools to identify novel biomarkers of treatment response and tailored therapeutic opportunities.

Prioritizing Crohn's disease genes by integrating association signals with gene expression implicates monocyte subsets.

Genome-wide association studies have identified ~170 loci associated with Crohn's disease (CD) and defining which genes drive these association signals is a major challenge. The primary aim of this study was to define which CD locus genes are most likely to be disease related. We developed a gene prioritization regression model (GPRM) by integrating complementary mRNA expression datasets, including bulk RNA-Seq from the terminal ileum of 302 newly diagnosed, untreated CD patients and controls, and in stimulated monocytes. Transcriptome-wide association and co-expression network analyses were performed on the ileal RNA-Seq datasets, identifying 40 genome-wide significant genes. Co-expression network analysis identified a single gene module, which was substantially enriched for CD locus genes and most highly expressed in monocytes. By including expression-based and epigenetic information, we refined likely CD genes to 2.5 prioritized genes per locus from an average of 7.8 total genes. We validated our model structure using cross-validation and our prioritization results by protein-association network analyses, which demonstrated significantly higher CD gene interactions for prioritized compared with non-prioritized genes. Although individual datasets cannot convey all of the information relevant to a disease, combining data from multiple relevant expression-based datasets improves prediction of disease genes and helps to further understanding of disease pathogenesis.

Functional variants in the LRRK2 gene confer shared effects on risk for Crohn's disease and Parkinson's disease.

Crohn's disease (CD), a form of inflammatory bowel disease, has a higher prevalence in Ashkenazi Jewish than in non-Jewish European populations. To define the role of nonsynonymous mutations, we performed exome sequencing of Ashkenazi Jewish patients with CD, followed by array-based genotyping and association analysis in 2066 CD cases and 3633 healthy controls. We detected association signals in the LRRK2 gene that conferred risk for CD (N2081D variant, P = 9.5 × 10-10) or protection from CD (N551K variant, tagging R1398H-associated haplotype, P = 3.3 × 10-8). These variants affected CD age of onset, disease location, LRRK2 activity, and autophagy. Bayesian network analysis of CD patient intestinal tissue further implicated LRRK2 in CD pathogenesis. Analysis of the extended LRRK2 locus in 24,570 CD cases, patients with Parkinson's disease (PD), and healthy controls revealed extensive pleiotropy, with shared genetic effects between CD and PD in both Ashkenazi Jewish and non-Jewish cohorts. The LRRK2 N2081D CD risk allele is located in the same kinase domain as G2019S, a mutation that is the major genetic cause of familial and sporadic PD. Like the G2019S mutation, the N2081D variant was associated with increased kinase activity, whereas neither N551K nor R1398H variants on the protective haplotype altered kinase activity. We also confirmed that R1398H, but not N551K, increased guanosine triphosphate binding and hydrolyzing enzyme (GTPase) activity, thereby deactivating LRRK2. The presence of shared LRRK2 alleles in CD and PD provides refined insight into disease mechanisms and may have major implications for the treatment of these two seemingly unrelated diseases.

Different tissue phagocytes sample apoptotic cells to direct distinct homeostasis programs.

Recognition and removal of apoptotic cells by professional phagocytes, including dendritic cells and macrophages, preserves immune self-tolerance and prevents chronic inflammation and autoimmune pathologies. The diverse array of phagocytes that reside within different tissues, combined with the necessarily prompt nature of apoptotic cell clearance, makes it difficult to study this process in situ. The full spectrum of functions executed by tissue-resident phagocytes in response to homeostatic apoptosis, therefore, remains unclear. Here we show that mouse apoptotic intestinal epithelial cells (IECs), which undergo continuous renewal to maintain optimal barrier and absorptive functions, are not merely extruded to maintain homeostatic cell numbers, but are also sampled by a single subset of dendritic cells and two macrophage subsets within a well-characterized network of phagocytes in the small intestinal lamina propria. Characterization of the transcriptome within each subset before and after in situ sampling of apoptotic IECs revealed gene expression signatures unique to each phagocyte, including macrophage-specific lipid metabolism and amino acid catabolism, and a dendritic-cell-specific program of regulatory CD4+ T-cell activation. A common 'suppression of inflammation' signature was noted, although the specific genes and pathways involved varied amongst dendritic cells and macrophages, reflecting specialized functions. Apoptotic IECs were trafficked to mesenteric lymph nodes exclusively by the dendritic cell subset and served as critical determinants for the induction of tolerogenic regulatory CD4+ T-cell differentiation. Several of the genes that were differentially expressed by phagocytes bearing apoptotic IECs overlapped with susceptibility genes for inflammatory bowel disease. Collectively, these findings provide new insights into the consequences of apoptotic cell sampling, advance our understanding of how homeostasis is maintained within the mucosa and set the stage for development of novel therapeutics to alleviate chronic inflammatory diseases such as inflammatory bowel disease.

Escape to Ferality: The Endoferal Origin of Weedy Rice from Crop Rice through De-Domestication.

Domestication is the hallmark of evolution and civilization and harnesses biodiversity through selection for specific traits. In regions where domesticated lines are grown near wild relatives, congeneric sources of aggressive weedy genotypes cause major economic losses. Thus, the origins of weedy genotypes where no congeneric species occur raise questions regarding management effectiveness and evolutionary mechanisms responsible for weedy population success. Since eradication in the 1970s, California growers avoided weedy rice through continuous flood culture and zero-tolerance guidelines, preventing the import, presence, and movement of weedy seeds. In 2003, after decades of no reported presence in California, a weedy rice population was confirmed in dry-seeded fields. Our objectives were to identify the origins and establishment of this population and pinpoint possible phenotypes involved. We show that California weedy rice is derived from a different genetic source among a broad range of AA genome Oryzas and is most recently diverged from O. sativa temperate japonica cultivated in California. In contrast, other weedy rice ecotypes in North America (Southern US) originate from weedy genotypes from China near wild Oryza, and are derived through existing crop-wild relative crosses. Analyses of morphological data show that California weedy rice subgroups have phenotypes like medium-grain or gourmet cultivars, but have colored pericarp, seed shattering, and awns like wild relatives, suggesting that reversion to non-domestic or wild-like traits can occur following domestication, despite apparent fixation of domestication alleles. Additionally, these results indicate that preventive methods focused on incoming weed sources through contamination may miss burgeoning weedy genotypes that rapidly adapt, establish, and proliferate. Investigating the common and unique evolutionary mechanisms underlying global weed origins and subsequent interactions with crop relatives sheds light on how weeds evolve and addresses broader questions regarding the stability of selection during domestication and crop improvement.

Improved integrative framework combining association data with gene expression features to prioritize Crohn's disease genes.

Genome-wide association studies in Crohn's disease (CD) have identified 140 genome-wide significant loci. However, identification of genes driving association signals remains challenging. Furthermore, genome-wide significant thresholds limit false positives at the expense of decreased sensitivity. In this study, we explored gene features contributing to CD pathogenicity, including gene-based association data from CD and autoimmune (AI) diseases, as well as gene expression features (eQTLs, epigenetic markers of expression and intestinal gene expression data). We developed an integrative model based on a CD reference gene set. This integrative approach outperformed gene-based association signals alone in identifying CD-related genes based on statistical validation, gene ontology enrichment, differential expression between M1 and M2 macrophages and a validation using genes causing monogenic forms of inflammatory bowel disease as a reference. Besides gene-level CD association P-values, association with AI diseases was the strongest predictor, highlighting generalized mechanisms of inflammation, and the interferon-γ pathway particularly. Within the 140 high-confidence CD regions, 598 of 1328 genes had low prioritization scores, highlighting genes unlikely to contribute to CD pathogenesis. For select regions, comparably high integrative model scores were observed for multiple genes. This is particularly evident for regions having extensive linkage disequilibrium such as the IBD5 locus. Our analyses provide a standardized reference for prioritizing potential CD-related genes, in regions with both highly significant and nominally significant gene-level association P-values. Our integrative model may be particularly valuable in prioritizing rare, potentially private, missense variants for which genome-wide evidence for association may be unattainable.