Inflammatory Markers at Birth and Risk of Early-Onset Inflammatory Bowel Disease
Jonas J. Rudbaek, Nis Borbye‐Lorenzen, Gry Poulsen, Adam Koziol, Kristin Skogstrand, Tine Jess
Abstract
The development of inflammatory bowel disease (IBD) may be preceded by a long preclinical phase characterized by several disease-associated immunologic changes, as observed for other immune-mediated inflammatory diseases.1Rudbaek J.J. et al.Nat Rev Gastroenterol Hepatol. 2023; 21: 86-100Crossref PubMed Scopus (10) Google Scholar As an example, disease-associated changes in systemic inflammatory proteins have been observed up to a decade before a diagnosis of Crohn’s disease (CD) or ulcerative colitis (UC).2Agrawal M. et al.EClinicalMedicine. 2021; 36100884Abstract Full Text Full Text PDF Scopus (58) Google Scholar,3Vestergaard M.V. et al.Cell Rep Med. 2023; 4101263Google Scholar Despite increasing evidence of early life events being important for later risk of disease,4Kim E.S. et al.Gastroenterology. 2021; 160: 1118-1130.e3Abstract Full Text Full Text PDF PubMed Scopus (43) Google Scholar,5Torres J. et al.Gut. 2020; 69: 42-51Crossref PubMed Scopus (123) Google Scholar no studies have, to our knowledge, linked systemic inflammation at time of birth with later risk of IBD. In a unique cohort of neonatal blood spots from individuals who developed IBD during childhood, we aimed to study how systemic inflammatory changes around birth are associated with development of (1) IBD before the age of 19 years and (2) very-early-onset (VEO) IBD diagnosed before the age of 6 years. Our cohort included 464 individuals (CD: n = 243, 46% female; UC: n = 221, 50% female) with pediatric-onset IBD diagnosed at a median age of 11 years (interquartile range, 7–15 years) identified via the Danish National Patient Registry. In the cohort, a total of 35 individuals (3.8% of the cohort) had mothers with IBD, and 26 of these individuals developed IBD later in life. In the VEO subcohort, 4.5% had a mother with IBD, and 75% of these developed IBD later in life. For each IBD case, we identified a non-IBD control individual from the general population based on time of sampling, sex, and gestational age or birthweight (ie, when gestational age was not available). Dried blood spots on case and control individuals (n = 928) were identified in the Danish National Biobank, which stores blood spots collected 4–6 days after birth from Danish individuals since 1982, thus providing an opportunity for investigating disease-associated changes at birth. Nine cytokines were measured in the fg/mL range using a highly sensitive multiplex immunoassay (S-PLEX, MesoScale) (Supplementary Table 1). Six of the 9 cytokines were measurable in at least 20% of individuals (interleukin [IL] 1β, IL-4, IL-6, IL-17A, interferon gamma, and tumor necrosis factor α), and further analyses were based on these. We performed univariate and multivariate logistic regression to assess whether the level of an individual inflammatory marker, used as a continuous variable, was associated with future risk of IBD, CD, or UC. Analyses were repeated for the subsample of patients with VEO-IBD. Results are reported as odds ratios (ORs) with 95% confidence intervals (CIs). The measured concentration in fg/mL was divided by 100 such that the shown OR corresponds to the risk associated with an increase of 0.1 pg/mL (or 100 fg/mL). Correction for multiple testing by false discovery rate (FDR) was performed in all analyses described. P values and Q values (ie, FDR-corrected P values) lower than .05 were considered statistically significant. In the overall cohort of pediatric-onset IBD, we found no statistically significant association between levels of individual neonatal inflammatory proteins and following risk of IBD (Figure 1A). This was also the case when studying CD and UC separately (Figure 1B and C). In a sensitivity analysis, cytokines were analyzed as categorical variables, with each marker divided into 4 groups based on concentration, which showed similar results (Supplementary Figure 1). We then confined our analyses to the 88 patients in the pediatric cohort who were diagnosed with VEO-IBD. Of these, 36 had CD (28% female), 52 had UC (54% female), and their median age at diagnosis was 2 years (range, 0–5 years). We assessed neonatal inflammatory levels in the 88 VEO-IBD patients and their matched control individuals, using cytokine concentrations as a continuous variable and adjusting for gestational age, birth year, mode of delivery (vaginal or cesarean birth), and sex. We observed a statistically significant protective effect of IL-4 on later development of VEO-IBD (OR, 0.43; 95% CI, 0.19–0.83; P = .02). In contrast, increased IL-17A was marginally significantly associated with increased risk of VEO-IBD (OR, 1.17; 95% CI, 1.01–1.42; P = .06) (Figure 1D). The significance was lost after adjustment for multiple testing, likely due to low power. To our knowledge, we present the first study investigating systemic inflammatory markers around birth and later risk of IBD. In a unique cohort of 464 pediatric IBD cases, we measured levels of interferon gamma, tumor necrosis factor α, IL-1β, IL-4, IL-6, and IL-17A a few days after birth. In the overall pediatric patient cohort diagnosed with IBD at a median age of 11 years, neonatal levels of cytokines were not associated with later disease occurrence. However, in the subcohort of 88 patients diagnosed with VEO-IBD before age 6 years, we observed a statistically significant association between increased levels of IL-4 and reduced risk of IBD. Increased IL-17A levels were, contrarily, borderline statistically significantly associated with an increased risk of IBD. IL-4, a cytokine generally associated with the T helper (Th) type 2 immune response, has been shown to be implicated in wound healing and tissue repair.6Gause W.C. et al.Nat Rev Immunol. 2013; 13: 607-614Crossref PubMed Scopus (365) Google Scholar The specific role of the Th2 response in IBD is still elusive7Imam T. et al.Front Immunol. 2018; 9: 1212Crossref PubMed Scopus (174) Google Scholar; however, IL-4–treated human macrophages have been observed to promote epithelial wound repair in vitro, and the transfer of such cells to an acute colitis mouse model reduced disease severity.8Jayme T.S. et al.Sci Adv. 2020; 6: 1-13Crossref Scopus (46) Google Scholar Thus, a decrease in IL-4 in VEO-IBD patients might indicate a reduced potential for homeostatic tissue repair. IL-17A, on the other hand, is a well-known proinflammatory cytokine, generally produced by Th17 and ILC3 cells, which has previously been associated with IBD as part of the IL-23/IL-17 axis and has been observed to be up-regulated in active IBD patients.9Schmitt H. et al.Front Immunol. 2021; 12: 1-12Crossref Scopus (96) Google Scholar,10Fujino S. et al.Gut. 2003; 52: 65-70Crossref PubMed Scopus (1451) Google Scholar The strength of this study is access to a unique population-based resource of neonatal blood spots from the Danish National Biobank, which, in the Danish PREDICT project, is linked to nationwide health information from a health care system that covers the entire population. This enables unique insights into disease development. Limitations include variations in measured concentrations due to varying sample storage time, difference in hematocrit levels, and differences in handling of samples over time. To mitigate this risk, we matched case and control individuals on time of sampling and adjusted for the year of birth, that is, year of sampling. Also, likely due to a lack of statistical power, significance disappeared after adjustment for multiple testing in the VEO-IBD analysis. In conclusion, in a unique sample of neonatal blood spots from 464 individuals who later developed IBD, we measured cytokine concentrations a few days after birth. Whereas neonatal cytokine levels were not associated with pediatric IBD in the overall cohort, we observed decreased levels of IL-4 and increased levels of IL-17A a few days after birth in patients who developed IBD before the age of 6 years. The observation of inflammatory changes already at birth in patients who develop IBD at an early age is novel and may shed light on disease occurrence and pathophysiology. The authors would like to thank Karin Skaarup and Nadja Riis for their assistance with the study. Jonas J. Rudbæk, MSc (Conceptualization: Lead; Formal analysis: Lead; Investigation: Lead; Methodology: Lead; Visualization: Lead; Writing – original draft: Lead; Writing – review & editing: Equal). Nis Borbye-Lorenzen, MSc, PhD (Conceptualization: Equal; Investigation: Supporting; Methodology: Supporting; Supervision: Supporting; Writing – review & editing: Equal). Gry Juul Poulsen, PhD (Formal analysis: Supporting; Methodology: Supporting; Writing – review & editing: Equal). Adam Koziol, PhD (Formal analysis: Supporting; Methodology: Supporting; Writing – review & editing: Equal). Kristin Skogstrand, PhD (Conceptualization: Equal; Investigation: Supporting; Methodology: Supporting; Supervision: Supporting; Writing – review & editing: Equal). Tine Jess, MD, PhD, DMSc (Conceptualization: Equal; Funding acquisition: Lead; Supervision: Lead; Writing – review & editing: Equal). Supplementary Table 1Analytical Parameters for the Performed Multiplex ImmunoassayAnalyteWorking range, fg/mLAssay CV, %In detection range, %aSample with a signal within the working range for the given analyte was considered to be within detection range.LowHighIntra-assayInterassayIL-1β105.0159,0007.313.595.1IL-215.485,8006.016.910.6IL-424.855,30010.810.978.3IL-619.139,50012.615.295.0IL-1037.1148,0007.620.917.4IL-12p70231.0484,0007.012.25.4IL-17A87.4210,00015.417.629.3IFN-γ14.028,8006.414.055.7TNF-α20.753,4006.213.270.7a Sample with a signal within the working range for the given analyte was considered to be within detection range. Open table in a new tab The study cohort was established via the Danish National Patient Registry and the Danish National Biobank. We included 464 individuals with predominantly pediatric-onset IBD, born during 1997–2018, and identified via the Danish National Patient Registry using International Classification of Diseases, 8th Revision codes (56301, 56302, 56308, 56309, 56319, 56904) and International Classification of Diseases, 10th Revision codes (DK50, DK51). To qualify for the cohort, case patients should have either 2 hospital/outpatient visits within 2 years or 2 inpatient contacts and have lived in Denmark in the 2 years before diagnosis. One healthy control individual was matched to each IBD case patient based on time of sampling, sex, and gestational age or birthweight (when information about gestational age was not available). Information about gestational age, birthday, mode of delivery (vaginal or cesarean birth) and maternal IBD status was collected from the Danish National Patient Register, the National Birth Register, and the Central Person Register. For all IBD case and control individuals, dried blood spots (DBSs) were identified in the Danish National Biobank. All DBSs had been sampled 4–6 days after birth as part of the routine screening program of all neonates born in Denmark since 1981 and were stored at –20°C. Two 3.2-mm disks were punched from the DBS into microtiter wells using a DBS Puncher instrument (PerkinElmer), and 130-μL extraction buffer (1× phosphate-buffered saline [Gibco], containing protease inhibitor [Roche cOmplete]) was added. Plates were sealed and incubated for 1 hour at room temperature with shaking (450 revolutions/min). The protein extract solution was transferred to 0.5-mL Matrix tubes (Thermo Scientific) and stored at –80°C until use. Samples were analyzed for the levels of cytokines interferon gamma (IFN-γ), IL-1β, IL-2, IL-4, IL-6, IL-10, IL-12p70, IL-17A, and tumor necrosis factor α (TNF-α) at the fg/mL range using S-PLEX proinflammatory panel 1 kit (no. K15396S, MesoScale Diagnostics) per the supplier’s instructions. All washing steps were performed on the Biotek Model 405 LS plate washer using the settings recommended by the supplier. Controls were made in-house from part of the calibrator solution in one batch and aliquoted into portions appropriate for one plate and stored at –20°C until use. Calibrator standards were made fresh for each plate. The concentrations of analytes were calculated from the calibrator curves of each plate using 4-parameter logistic regression in the Meso Scale Discovery Workbench software. Intra-assay variation was calculated from 16 replicates of high controls (calibrator 4) on the same plate. Interassay variation was calculated from 2 replicates of high and low controls (calibrators 4 and 3) on each plate, 13 plates in total. Lower limits of detection (LLODs) were calculated as 2.5 standard deviations from 26 repeated measurements of the 0 standard. Outliers within the repeated measurements were identified using the Tukey method (also known as the Tukey fence method) and were removed before the calculation of the LLOD. The upper detection limits (ULODs) were defined as the highest calibrator concentrations. Concentrations for measurements below the LLOD were defined as the LLOD for the given analyte divided by 2; similarly, concentrations above the ULOD were defined as the ULOD for the given analyte times 2 (Supplementary Table 1). All statistical analyses were performed in R Studio using the following packages: tidyverse, gt, gtsummary, metafor, patchwork, DescTools, and ggpubr. In all analyses described, P values and Q values (FDR-corrected P values) lower than .05 were deemed to be statistically significant. Differences between CD patients, UC patients, and healthy control individuals in the distribution of sex was estimated using the Pearson chi-square test, and difference in age at diagnosis (in years) and gestational age (in days) were estimated using the Kruskal-Wallis rank sum test. Before statistical modeling, cytokine concentrations in fg/mL were divided by 100, such that the shown calculated ORs correspond to the risk associated with an increase of 0.1 pg/mL (or 100 fg/mL). We then performed both univariate and multivariate logistic regression to estimate the association between cytokine concentrations and later IBD, CD, or UC development. In the multivariate model, all 6 cytokines were included as fixed explanatory variables, and the model was adjusted for gestational age (in days), year of birth (in years), mode of delivery (vaginal or cesarean birth), and sex. ORs and 95% CIs were calculated, and P values were adjusted for multiple testing by FDR correction. Furthermore, as a sensitivity analysis, the multivariate analysis was performed as described, but with cytokines grouped into 4 groups corresponding to increasing levels of cytokine abundance. For sparsely abundant cytokines (see Supplementary Table 1: IL-4, IL-17a, IFN-γ, and TNF-α) group 1 comprised cytokines that were measured below the LLOD. The remaining groups (groups 2–4) were then grouped by percentiles (P1 = 33.33%, P2 = 66.66%). In contrast, highly abundant cytokines (IL-1β and IL-6) were grouped into quartiles based on cytokine concentrations. Exploratory analyses were performed to investigate associations between the measured cytokine levels and VEO cases of IBD. VEO IBD was defined as individuals diagnosed with IBD before the age of 6 years. The subcohort was analyzed by univariate and multivariate logistic regression using cytokine concentrations in fg/mL as a continuous variable. The analyses were performed as described for the full cohort. The present project has been approved by the Regional Ethical Committee of the Capital Region of Denmark and by the Scientific Board of the Danish National Biobank.