Summary
Objective There are few prospective data exploring the combined effect of genetic and lifestyle factors on the incidence of inflammatory bowel disease (IBD), and whether genetic risk of IBD may be mitigated by lifestyle remains unclear.
Design We conducted a prospective cohort study based on the UK Biobank to examine the associations across genetic risk, modifiable lifestyle factors, and risk of Crohn’s disease (CD) and ulcerative colitis (UC). Genetic susceptibility to CD and UC was estimated by polygenic risk scores using common genetic variants identified by genome-wide association studies and was further categorized into high, intermediate, and low genetic risk categories. Weighted unhealthy lifestyle scores were constructed based on disease-related lifestyle factors, including ever smoking, unhealthy diet, physical inactivity, obesity, and abnormal sleep duration, and were categorized into ‘favorable’, ‘intermediate’, and ‘unfavorable’ categories. The Cox proportional hazard regression model was used to estimate the HRs and 95% CIs for their associations and accumulative risk for developing CD and UC were estimated for each risk group.
Results During a median follow-up of 12.0 years, 707 CD and 1576 UC cases were diagnosed. Genetic risk and unhealthy lifestyle categories were monotonically associated with CD and UC risk and there was no multiplicative interaction between them. Compared with participants with low genetic risk, the hazard ratios (HRs) of CD and UC were 2.24 (95% confidence interval [CI] 1.75-2.86) and 2.15 (95% CI 1.82-2.53) for those with high genetic risk, respectively. The HRs of CD and UC for individuals in unfavorable lifestyle category were 1.94 (95% CI 1.61-2.33) and 1.98 (95% CI 1.73-2.27), respectively, compared with those in favorable category. When considering genetic risk and lifestyle jointly, the HRs of individuals with high genetic risk but a favorable lifestyle (2.33, 95% CI 1.58-3.44 for CD, and 2.05, 95% CI, 1.58-2.66 for UC) were reduced nearly by half comparing to those with high genetic risk but an unfavorable lifestyle (4.40, 95% CI, 2.91-6.66 for CD and 4.44, 95% CI, 3.34-5.91 for UC).
Conclusion Genetic and lifestyle factors were independently associated with susceptibility to CD and UC. Participants at high genetic risk could reduce nearly 50% risk of CD and UC by adherence to a favorable lifestyle.
Funding XL: the Natural Science Fund for Distinguished Young Scholars of Zhejiang Province (LR22H260001). XYW: National Natural Science Foundation of China (81970494) and Key Project of Research and Development Plan of Hunan Province(2019SK2041); SCL: the Swedish Heart-Lung Foundation (Hjärt-Lungfonden, 20210351), the Swedish Research Council (Vetenskapsrådet, 2019-00977), and the Swedish Cancer Society (Cancerfonden); ET: CRUK Career Development Fellowship (C31250/A22804); KFD: Project of the regional diagnosis and treatment center of the Health Planning Committee (No. JBZX-201903).
What is already known on this topic Previous studies have identified a number of genetic variants and several modifiable risk factors for IBD. However, there is a lack of studies on the combined effects of genetic and lifestyle factors on IBD risk.
What this study adds In this prospective cohort study, genetic risk and modifiable lifestyle factors were independently associated with the risk of incident Crohn’s disease and ulcerative colitis. Participants at high genetic risk could reduce nearly 50% of their genetic susceptibility to Crohn’s disease and ulcerative colitis by adherence to a favorable lifestyle.
How this study might affect research, practice or policy Promoting a healthy lifestyle is an effective strategy to lower the incidence of these diseases, especially among those with high-risk genetic background.
Introduction
Inflammatory bowel disease (IBD), which includes two main subtypes (Crohn’s disease, CD, and ulcerative colitis, UC), is a global health problem with substantial disease burden, especially in the industrialized countries.1,2 Although onset in childhood and early adulthood is well recognized, epidemiological studies now highlight the increasing incidence and prevalence of IBD onset in middle age, or later life. Compared to IBD in children or adolescents, the etiology of adult-onset IBD is believed to be more multifactorial, with genetic and environmental factors playing important roles in its development.3,4
Genome-wide association studies have identified over a hundred of risk loci, such as TYK2, IL2RA, and IL23R, to be associated with IBD.3,5-9 Although a single genetic variant accounts for only a small fraction of the genetic variability of IBD, polygenic risk scores combining multiple risk loci can be used as an indicator to identify individuals at higher genetic susceptibility to IBD.10 Compared with rare genetic mutations with larger effect (e.g., NOD2), polygenic risk scores can identify a larger fraction of population at comparable or greater disease risk, which poses opportunities for clinical utility. However, it is largely unknown whether a PRS of IBD can identify individuals at high genetic risk for potential personalized prevention via adoption of healthy lifestyles in later life.
Observational studies have identified several potentially modifiable risk factors in relation to IBD, including cigarette smoking, unhealthy diet, physical inactivity, obesity, and abnormal sleep duration.4,11 The associations of these lifestyle factors with the risk of CD and UC appear complex. For instance, active smoking has been reported to be protective against UC but risky for CD.11 The association between alcohol drinking and IBD is inconclusive and remains elusive.12 Sleep is one of the common lifestyle factors in maintaining physical and psychological health, however, whether sleep behavior is associated with IBD risk has been scarcely investigated.13 Comprehensive appraisal of the associations between these modifiable lifestyle factors and IBD risk will deepen the understanding of the etiology of IBD and provide clues for IBD prevention.
So far, there is lack of comprehensive investigation on the combined effect of genetic and lifestyle factors on the development of IBD and its subtypes. Herein, we conducted a prospective cohort study based on the UK Biobank to examine the associations across genetic risk, modifiable lifestyle factors, and risk of IBD, to test whether there is any multiplicative interaction between genetic risk and lifestyle factors, and to figure out to what extent the genetic risk of IBD may be mitigated by adherence to healthy lifestyle choices.
Methods
Study population
This cohort study is based on data collected from the UK Biobank including approximately 500,000 participants recruited across the United Kingdom between 2006 and 2010.14 Individuals of non-European ancestry (due to limited numbers and to minimize population structure bias), without genetic information, or with baseline IBD diagnosis, new onset IBD within 1-year follow-up, or unclear IBD diagnosis were excluded, leaving 453,492 individuals (Figure 1).
Flowchart of study population selection in the UK Biobank study. CD, Crohn’s disease; IBD, inflammatory bowel disease; UC, ulcerative colitis.
Genetic risk profiling
We applied two strategies to estimate the genetic susceptibility to CD and UC for UK Biobank population. We firstly constructed a polygenetic risk score (PRS), using the common genetic variants that were identified to be strongly associated with CD and UC (P<5×10−8) from a genome-wide association meta-analysis of up to 86,640 individuals of European ancestry, including 5,956 CD and 6,968 UC cases.3 After removing genetic variants in linkage disequilibrium (r2>0.001), 51 and 30 independent single nucleotide polymorphisms (SNPs) were used to calculate the PRS of CD and UC, respectively. Polygenic risk scores were constructed for each participant by summing up the number of risk-increasing alleles for each SNP weighted by effect size on genetic liability to CD or UC (Supplementary Table 1). Because using merely the genome-wide significant SNPs may omit some correlated informative signals that are independently associated with CD and UC, we additionally constructed genomic risk scores by including all SNPs at suggestive significance level (P<1×10−5) reported by the GWAS. Genomic risk scores were calculated by using the LDpred2.15 Either polygenic risk score or genomic risk score with better stratification ability was taken froward to proxy the genetic susceptibility of CD and UC, and was further used to categorize the low (the lowest quintile), intermediate (quintiles 2 to 4), and high (highest quintile) genetic risk groups.
Modifiable lifestyle factors
Six common lifestyle factors, including cigarette smoking, diet, alcohol consumption, physical inactivity, BMI, and sleep duration, were examined for their associations with CD and UC risk respectively. These lifestyle factors were chosen based on pre-existing evidence on their associations with either CD or UC, as reported by a recent umbrella review and cohort studies.11,16 Detailed information on definitions of common lifestyle factors is displayed in Supplementary Method and Supplementary Table 2. Healthy lifestyle scores (HLS) were constructed based on aforementioned lifestyle factors. Individuals were assigned 1 point for each of lifestyle behaviors if they were classified into the healthy group. A higher lifestyle score indicates higher adherence to healthy lifestyle. The unweighted lifestyle score was categorized as favorable (4 or 5 healthy lifestyle factors), intermediate (3 healthy lifestyle factors), and unfavorable (0-2 healthy lifestyle factor) lifestyles. We then constructed a weighted standardized healthy lifestyle score based on the β coefficient of each lifestyle factor in the Cox proportional hazards model adjusted for age, age-square, sex, education, Townsend deprivation index, and first 20 principal components of ancestry using the formula and then the weighted standardized healthy lifestyle score was categorized into unfavorable (the lowest quintile), intermediate (quintiles 2-4), and favorable (the highest quintile) groups.
Cases ascertainment and follow-up
Diagnostic information was obtained from the primary care data and hospital inpatient records. Incident CD and UC cases were ascertained by a primary or secondary diagnosis defined by corresponding International Classification of Diseases codes (ICD-9: 555, 556; ICD-10: K50, K51). Participants were followed up from the baseline (2006-2010) until the date of first diagnosis of IBD, date of death, date of loss to follow-up, the last date of hospital admission (i.e., HES and SMR: 31 March 2021, and PEDW: 28 February 2018), whichever came first. Disease locations were obtained from diagnosis records for subgroup analyses.17
Statistical analysis
We used the Cox proportional hazard regression model to examine the associations of genetic risk categories, lifestyle categories, and genetic risk and lifestyle combined categories (9 categories with high genetic risk and unfavorable lifestyle as reference) with risk of incident CD and UC. The model was adjusted for age, age-square, sex, education, Townsend deprivation index, Charlson comorbidities index and first 20 principal components of ancestry. The interactions between lifestyle factors and polygenic risk scores were also examined using a multiplicative interaction model. The proportionality of hazards assumption was assessed using the Schoenfeld residuals method and found to be satisfied (P>0.15). To examine the consistency of the association in subpopulations, we conducted stratification analyses by age (≥60 and <60 years), sex (female and male), education attainment (≥college/university and < college/university), and the tertiles of Townsend deprivation index (from low to high, T1-3). We also stratified the analysis on the associations of the healthy lifestyle categories with CD and UC risk by genetic risk. Sensitivity analyses and subgroup analyses by considering disease locations and age of diagnosis of UC and CD (to retrospectively include prevalent cases and stratify analysis by age of onset) were also performed to thoroughly examine their complex associations (Supplementary Method). The cumulative incidence of CD and UC by categories of genetic risk and lifestyle scores were obtained using the cumulative incidence function of competing risk regression.18 All tests were two-sided and the association with the P value <0.05 was deemed significant. All analyses were performed using R software, version 3.6.3.
Results
Table 1 shows baseline characteristics of included participants by incident disease status. Over a median follow-up of 12.0 years (interquartile range, 11.2-12.7 years), 707 CD and 1576 UC cases were diagnosed. The mean age of diagnosis was 64 (range: 43-82 for CD and 65 (range: 43-82) for UC.
Characteristics of participants by incident Crohn’s disease and ulcerative colitis status.
For genetic susceptibility, both polygenic risk score and genomic risk score showed significantly associations with the risk of CD and UC (Supplementary Table 4). Comparing to the PRS, the genomic risk score showed no further improvement on the stratification of genetic risk groups (Table 2, Supplementary Table 4); therefore, only the PRS was used in the joint analysis. The PRS was normally distributed (Supplementary Figure 1) and showed no associations with lifestyle factors with the exception for an association between polygenic risk score of CD and smoking status (Supplementary Table 5). Risk of incident CD and UC increased across genetic risk categories (low to high) in a linear fashion (Table 2). Compared with participants with low genetic risk, the hazard ratios (HRs) of CD and UC were 2.24 (95% confidence interval, CI: 1.75-2.86, P<0.001) and 2.15 (95% CI: 1.82-2.53, P<0.001), for those with high genetic risk, respectively. The associations remained significant after additional adjustment for lifestyle factors. The same pattern of associations was observed in the analysis using the genetic risk quintiles instead of categories (Supplementary Table 4).
Risk of incident Crohn’s disease and ulcerative colitis according to genetic risk and lifestyle categories.
The associations of individual lifestyle factors with the risk of CD and UC are presented in Table 3. Less healthy behavior was in general associated with increased risk of CD and UC, compared with those meeting healthy lifestyle guidelines (the reference category) for each component of healthy lifestyles although not all risk estimates were statistically significant. Exceptions were noted for alcohol drinking, which was neither associated with CD nor associated with UC (Supplementary Table 6). Given that there was no well-established evidence supporting their associations from previous evidence either,12 we therefore excluded alcohol consumption from the construction of healthy lifestyle scores. For UC, the association with obesity was not statistically significant (P=0.29) and irregular physical activity was marginally associated with an increased risk of UC (P=0.068); nevertheless, given that obesity and physical activity were well-established lifestyle factors related to IBD based on previous evidence,19,20 we decided to include these variables for the construction of healthy lifestyle scores. The associations of smoking status and frequency with the risk of CD and UC at different disease locations and age of diagnosis are shown with details in Supplementary Tables 7-11. Briefly, both previous and current smoking were consistently associated with an increased risk of CD and UC among older populations (Supplementary Table 11). We therefore simplified the smoking exposure into ever versus never smoking behaviors as one of the components of a healthier lifestyle.
Risk of incident Crohn’s disease and ulcerative colitis with each lifestyle factor
Having a healthier lifestyle score was significantly related to a reduced risk of CD and UC in a dose-response manner (P for trend < 0.001) (Table 2). The HRs of CD and UC for individuals in the unfavorable category were 1.94 (95% CI: 1.61-2.33; P<0.001) and 1.98 (95% CI: 1.73-2.27; P<0.001), respectively, compared with those in the favorable category. The associations did not change in the sensitivity analysis with further adjustment for genetic risk (Table 2), in the analysis using the number of healthy lifestyle factors instead of categories (Supplementary Table 12) and in the analysis using the unweighted healthy lifestyle score (Supplementary Table 13). The cumulative incidence rate of CD and UC during the follow-up was higher in the group with an unfavorable lifestyle compared to the group with a favorable lifestyle (Figure 2).
Cumulative incidence plot of the risk of incident Crohn’s disease (A) and ulcerative colitis (B) by polygenic risk score categories and healthy lifestyle categories.
In the analysis of joint categories for genetic risk and healthy lifestyle, the HR of CD and UC showed a linear increase with increasing genetic risk and decreasing healthy lifestyle score (Figure 3). Compared with individuals with low genetic risk and favorable lifestyle, the HRs of CD and UC for those with high genetic risk and unfavorable lifestyle were 4.40 (95% CI: 2.91-6.66; P<0.001) and 4.44 (95% CI: 3.34-5.91; P<0.001), respectively. We observed no significant difference in the HRs of CD or UC between the group of high-genetic risk but having a favorable lifestyle (HR=2.33, 95%CI: 1.58-3.44 for CD and HR=2.05, 95%CI: 1.58-2.66 for UC) and the group of low-genetic risk but having an unfavorable lifestyle (HR=2.32, 95%CI: 1.44-3.74 for CD and HR=1.77, 95%CI: 1.23-2.56 for UC). The analysis on the associations of healthy lifestyle categories with incident CD and UC risk in groups defined by genetic risk confirmed that the unfavorable lifestyle was associated with higher risk of CD and UC across all genetic groups (Table 4). Specifically, in individuals with low genetic risk, the HRs of CD and UC were 2.32 (95% CI, 1.42-3.81) and 1.70 (95% CI, 1.17-2.47) for participants with an unfavorable lifestyle compared to those with a favorable lifestyle. We did not detect any multiplicative interaction between the genetic risk and the weighted healthy lifestyle score (p=0.85 for CD and p=0.87 for UC, Supplementary Table 14). The observed associations remained statistically significant in a series of sensitivity analyses (Supplementary Tables 15-20).
Risk of incident Crohn’s disease and ulcerative colitis according to lifestyle categories within each genetic risk category.
Risk of incident Crohn’s disease and ulcerative colitis by joint categorization for genetic risk and healthy lifestyle score. CI, confidence interval; HR, hazard ratio.
We calculated the cumulative risk of CD and UC over 12 years for each group defined jointly by genetic risk and healthy lifestyle scores (Figure 4). Compared to those with low genetic risk and favorable lifestyle (accumulative risk: 0.08% for CD, 0.18% for UC), individuals with high genetic risk and unfavorable lifestyle had 4.87 times higher accumulative risk of CD (equivalent to an excess risk of 0.31% due to high genetic susceptibility and unfavorable lifestyle together) and 5.28 times higher accumulative risk of UC (equivalent to an excess risk of 0.77% due to high genetic susceptibility and unfavorable lifestyle together).
Cumulative incidence plot of the risk of incident Crohn’s disease (A) and ulcerative colitis (B) by joint categorization for genetic risk and healthy lifestyle. CI, confidence interval.
Discussion
We conducted a cohort study using data from the UK Biobank to investigate the associations across genetic susceptibility, modifiable lifestyle factors, and risk of CD and UC. We found that a polygenic risk score and modifiable lifestyle factors were independently associated with risk of incident CD and UC. No multiplicative interaction was observed between the polygenic risk score and lifestyle scores. Individuals with high genetic risk and an unfavorable lifestyle had elevated risk of CD and UC compared to their counterparts with low genetic risk and a favorable lifestyle. The cumulative incidence of CD and UC over 12 years was approximately five times higher in those with highest genetic risk and unfavorable lifestyle compared to those with the lowest genetic risk and favorable lifestyle.
Even though no studies have been conducted to examine the effect of overall healthy lifestyle on CD or UC risk, a high adherence to a healthy lifestyle was associated with reduced mortality in CD or UC patients from three large cohort studies.21 For individual lifestyle factors, there are epidemiological studies assessing the associations of each lifestyle factor with CD and UC risk. Current smoking was found to be positively associated with risk of CD, but inversely associated with risk of UC in a meta-analysis of 9 and 13 studies, respectively.22 However, most included studies are cross-sectional studies where residual confounding and reverse causality could not be eliminated. In a subsequent prospective cohort study using data from Nurses’ Health Study (NHS) and Nurses’ Health Study II, the risk of CD was observed to increase in former and current smokers after an over 18 years’ follow-up period.23 A positive association between smoking and UC risk was observed in some subsequent studies.24,25 Smoking has been also associated with progression of CD, 26 but not with that of UC.27 In this study, we performed a series of sensitivity/subgroup analyses by considering disease locations and age of diagnosis to thoroughly examine their complex associations with smoking. Our study verified that both previous or current smoking were consistently associated with increased risk of CD regardless of disease locations and age of onset. However, the association between smoking and UC appears complex, in which smoking was associated with reduced risk for early-onset UC (<=20 years) but increased risk for later-onset UC (>40 years). Our previous study showed that smoking habit influences the age at diagnosis and changes in disease extent in UC 28. The mechanisms of the observed inverse association are not clear, but reverse causality is again an important point to consider. Since the initial age structure of UK biobank cohort is older, we therefore concentrated on older onset cases and considered smoking as a risk factor of CD and UC.
Evidence on the association between obesity and CD is inconsistent with a positive association in a meta-analysis of 5 cohort studies20 but a null finding in another meta-analysis29. In a recent prospective analysis of 5 cohorts, obesity defined by BMI was associated with an increased risk of older-onset CD but not UC.30 A recent Mendelian randomization study found that genetically predicted higher BMI and body fat percentage was associated with an increased risk of CD, but that genetically predicted higher BMI was associated with a lower risk of UC.31 Regular physical activity has been associated with a lower risk of CD, but not UC.19,32 Even though a high adherence to Mediterranean diet has associated with a low risk of CD 33 and abnormal sleep duration has been associated with a high risk of UC,34 there are few corresponding prior studies that examined these associations jointly. Our analysis using data from UK Biobank verified the associations of smoking, physical inactivity, unhealthy diet, and abnormal sleep duration with CD risk, with the exception of alcohol consumption. As for UC, this cohort study found positive associations between all aforementioned lifestyle factors with risk of UC, but no significant association for alcohol drinking either. The null association of alcohol drinking with the risk of CD and UC reported by the present study is consistent with evidence from a recent cohort study.12
To our knowledge, no previous studies have examined the association of a combination of healthy lifestyle and multiple genetic factors with risk of incident CD and UC. However, the interaction effects between genes and environmental factors, such as smoking and certain dietary nutrients, on CD and UC have been assessed. A study including 19,735 IBD cases (10,856 CD cases and 8879 UC cases) of known smoking status found that two variants in HLA and NOD2 gene regions interacted with smoking in influencing CD risk and smoking modified the disease risk of some variants in opposite directions for CD vs UC,35 which indicated that the effects of smoking on IBD risk depend on genetic variants. Nevertheless, an increased risk of CD and a decreased risk of UC were found in smokers in twin or sibling studies where the cases and controls shared genetic risk for the disease 36,37. The interaction effects were also observed for certain dietary nutrients, such as dietary fatty acids, potassium, and iron intake; however, these findings are far from being established to determine the gene-environment interaction on CD and UC risk.38 Our study found independent associations of genetic risk and healthy lifestyle with IBD risk, but no overall interaction for their joint effects on IBD risk. Among individuals with high genetic risk, those with unfavorable lifestyle had double the risk of CD or UC compared to those with favorable lifestyle. This finding has important clinical implications by indicating that promoting a healthy lifestyle is an effective strategy to lower incidence of these diseases, even among those with high-risk genetic background.
The strengths of this study include the joint analysis of the genetic and lifestyle factors to gain a comprehensive understanding on the risk of CD and UC, in which polygenetic risk scores and healthy lifestyle scores were constructed to examine their associations with the disease risk in a large prospective cohort of UK Biobank participants.3 We made efforts to account for additional genetic susceptibility that were not captured by genome-wide significant SNPs, while the more sophisticated genomic risk score showed no superior capacity comparing to the simple PRS in risk stratification. Limitations of the present study should also be acknowledged when interpreting the findings. Firstly, adherence to a healthy lifestyle might change in the follow-up and influence the association estimation. Nonetheless, the bias caused by the change should be non-differential and therefore attenuate the estimates in a conservative way due to the prospective nature of the design. Secondly, even though important known confounders were adjusted in the models, there is possible residual confounding. Thus, the causality of the association for lifestyle factors cannot be exclusively determined. Thirdly, information on lifestyle factors was partially collected via a self-administrative questionnaire survey, and thus the probability of misclassifications of lifestyle factors and recall bias cannot be fully excluded. This non-differential bias is likely to drive these findings toward to null. Additionally, misclassification of outcome possibly caused by cases undocumented in medical records could attenuate the effect estimates. Fourthly, given that the current analyses were confined to individuals of European ancestry and concentrated more on older population, our findings may not be generalizable to other populations of different ethnicities and/ or ages.
In summary, both high genetic risk and an unfavorable lifestyle were associated with increased risk of CD and UC among adults without IBD. An unfavorable lifestyle was associated with higher risk of CD and UC in individuals regardless of genetic strata. Participants at high genetic risk could reduce nearly 50% risk of CD and UC by adherence to a favorable lifestyle.
Data Availability
Researchers can request the data we used from the UK Biobank (www.ukbiobank.ac.uk/).
Additional information
Contributors
XL, JC and XYW conceptualized the project. YHS, SY and XJC performed the analyses and wrote the first draft. JS, LLY, LJW, XXK, TH, XK, GH, MD, KFD, SCL, JS, ET and ELG helped with the review and editing. XL, JC and XYW revised the final version. XL, JC and XYW contributed equally to this work and should be considered co-corresponding authors. All authors critically reviewed the manuscript for important intellectual content. XL is the study guarantor. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.
No competing interests
“All authors have completed the ICMJE uniform disclosure form at http://www.icmje.org/disclosure-of-interest/ and declare: no support from any organization for the submitted work; no financial relationships with any organizations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.”
Ethical approval
The manuscript didn’t previously publish or be under consideration for publication elsewhere. And the ethical approval was granted for the UK Biobank by the North West-Haydock Research Ethics Committee (REC reference: 16/NW/0274).
Data sharing
Researchers can request the data we used from the UK Biobank (www.ukbiobank.ac.uk/).
Transparency
The lead author (XL) affirms that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned have been explained.
Dissemination to participants and related patient and public communities
The results of the research will be disseminated to the public through broadcasts, popular science articles, and newspapers.
Acknowledgement
This research was conducted using the UK Biobank study under Application Number 66354. We want to thank all UK Biobank participants and the management team for their participation and assistance.