Background & aims

Crohn's disease is a relapsing and remitting inflammatory disorder with a variable clinical course. Although most patients present with an inflammatory phenotype (B1), approximately 20% of patients rapidly progress to complicated disease, which includes stricturing (B2), within 5 years. We analyzed DNA methylation patterns in blood samples of pediatric patients with Crohn's disease at diagnosis and later time points to identify changes that associate with and might contribute to disease development and progression.

Methods

We obtained blood samples from 164 pediatric patients (1-17 years old) with Crohn's disease (B1 or B2) who participated in a North American study and were followed for 5 years. Participants without intestinal inflammation or symptoms served as controls (n = 74). DNA methylation patterns were analyzed in samples collected at time of diagnosis and 1-3 years later at approximately 850,000 sites. We used genetic association and the concept of Mendelian randomization to identify changes in DNA methylation patterns that might contribute to the development of or result from Crohn's disease.

Results

We identified 1189 5'-cytosine-phosphate-guanosine-3' (CpG) sites that were differentially methylated between patients with Crohn's disease (at diagnosis) and controls. Methylation changes at these sites correlated with plasma levels of C-reactive protein. A comparison of methylation profiles of DNA collected at diagnosis of Crohn's disease vs during the follow-up period showed that, during treatment, alterations identified in methylation profiles at the time of diagnosis of Crohn's disease more closely resembled patterns observed in controls, irrespective of disease progression to B2. We identified methylation changes at 3 CpG sites that might contribute to the development of Crohn's disease. Most CpG methylation changes associated with Crohn's disease disappeared with treatment of inflammation and might be a result of Crohn's disease.

Conclusions

Methylation patterns observed in blood samples from patients with Crohn's disease accompany acute inflammation; with treatment, these change to resemble methylation patterns observed in patients without intestinal inflammation. These findings indicate that Crohn's disease-associated patterns of DNA methylation observed in blood samples are a result of the inflammatory features of the disease and are less likely to contribute to disease development or progression.

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.

Molecular mechanisms driving disease course and response to therapy in ulcerative colitis (UC) are not well understood. Here, we use RNAseq to define pre-treatment rectal gene expression, and fecal microbiota profiles, in 206 pediatric UC patients receiving standardised therapy. We validate our key findings in adult and paediatric UC cohorts of 408 participants. We observe a marked suppression of mitochondrial genes and function across cohorts in active UC, and that increasing disease severity is notable for enrichment of adenoma/adenocarcinoma and innate immune genes. A subset of severity genes improves prediction of corticosteroid-induced remission in the discovery cohort; this gene signature is also associated with response to anti-TNFα and anti-α₄β₇ integrin in adults. The severity and therapeutic response gene signatures were in turn associated with shifts in microbes previously implicated in mucosal homeostasis. Our data provide insights into UC pathogenesis, and may prioritise future therapies for nonresponders to current approaches.

BACKGROUND:Lack of evidence-based outcomes data leads to uncertainty in developing treatment regimens in children who are newly diagnosed with ulcerative colitis. We hypothesised that pretreatment clinical, transcriptomic, and microbial factors predict disease course. METHODS:In this inception cohort study, we recruited paediatric patients aged 4-17 years with newly diagnosed ulcerative colitis from 29 centres in the USA and Canada. Patients initially received standardised mesalazine or corticosteroids, with pre-established criteria for escalation to immunomodulators (ie, thiopurines) or anti-tumor necrosis factor-α (TNFα) therapy. We used RNA sequencing to define rectal gene expression before treatment, and 16S sequencing to characterise rectal and faecal microbiota. The primary outcome was week 52 corticosteroid-free remission with no therapy beyond mesalazine. We assessed factors associated with the primary outcome using logistic regression models of the per-protocol population. This study is registered with ClinicalTrials.gov, number NCT01536535. FINDINGS:Between July 10, 2012, and April 21, 2015, of 467 patients recruited, 428 started medical therapy, of whom 400 (93%) were evaluable at 52 weeks and 386 (90%) completed the study period with no protocol violations. 150 (38%) of 400 participants achieved week 52 corticosteroid-free remission, of whom 147 (98%) were taking mesalazine and three (2%) were taking no medication. 74 (19%) of 400 were escalated to immunomodulators alone, 123 (31%) anti-TNFα therapy, and 25 (6%) colectomy. Low baseline clinical severity, high baseline haemoglobin, and week 4 clinical remission were associated with achieving week 52 corticosteroid-free remission (n=386, logistic model area under the curve [AUC] 0·70, 95% CI 0·65-0·75; specificity 77%, 95% CI 71-82). Baseline severity and remission by week 4 were validated in an independent cohort of 274 paediatric patients with newly diagnosed ulcerative colitis. After adjusting for clinical predictors, an antimicrobial peptide gene signature (odds ratio [OR] 0·57, 95% CI 0·39-0·81; p=0·002) and abundance of Ruminococcaceae (OR 1·43, 1·02-2·00; p=0·04), and Sutterella (OR 0·81, 0·65-1·00; p=0·05) were independently associated with week 52 corticosteroid-free remission. INTERPRETATION:Our findings support the utility of initial clinical activity and treatment response by 4 weeks to predict week 52 corticosteroid-free remission with mesalazine alone in children who are newly diagnosed with ulcerative colitis. The development of personalised clinical and biological signatures holds the promise of informing ulcerative colitis therapeutic decisions. FUNDING:US National Institutes of Health.