Parental inflammatory bowel disease and autism in the offspring: Triangulating the evidence using four complementary study designs.

Importance: Evidence linking parental diagnoses of inflammatory bowel disease (IBD) with offspring autism is inconclusive. Objective: To investigate associations between parental diagnoses of IBD and offspring autism and elucidate their underlying aetiology by conducting four complementary studies. Design, Setting and Participants: (1) Nationwide population-based cohort study using Swedish registers to examine associations between parental IBD diagnoses and autism diagnoses in offspring, (2) Linkage disequilibrium (LD)-score regression to estimate the genetic correlation between the phenotypes. (3) Polygenic risk score (PRS) analyses in the Avon Longitudinal Study of Parents and Children (ALSPAC) to investigate associations between maternal genetic liability to IBD and autism factor mean score in offspring. (4) Two-sample Mendelian randomization (MR) to assess bidirectional causal links between genetic liability to IBD and autism. Results: Observational analyses provided evidence of an association between parental IBD diagnoses and offspring autism diagnosis in mutually adjusted models (maternal: OR= 1.32; 95%CI: 1.25 to 1.40; p<0.001; paternal: OR= 1.09; 95%CI: 1.02 to 1.17; p=0.012; n= 2 324 227, 52.3% male). PRS analyses in ALSPAC indicated associations between maternal PRS for IBD subtypes and a measure of broad autism phenotype, autism factor mean score, in the offspring (UC: b=0.02; 95%CI: 0.003 to 0.05; p= 0.02; R2= 0.06; Crohn's: b= 0.03; 95%CI: 0.01 to 0.05; p= 0.004; R2= 0.06; n= 7357, 50.3% male).MR analyses provided evidence of a potential causal effect of genetic liability for IBD, especially ulcerative colitis, on autism (OR= 1.03; 95%CI: 1.01 to 1.06). There was little evidence to suggest a causal effect of genetic liability to autism on risk of IBD, or a genetic correlation between the two conditions. Conclusions and relevance: Triangulating evidence from a nationwide register-based cohort study, genetic correlation, polygenic risk score analyses and MR, we found evidence of a potentially causal link between parental, particularly maternal, diagnoses and genetic liability to IBD and offspring autism. Perinatal immune system dysregulation, micronutrient malabsorption and anaemia may be implicated.


4
of a potentially causal link between parental, particularly maternal, diagnoses and genetic liability to IBD and offspring autism. Perinatal immune system dysregulation, micronutrient malabsorption and anaemia may be implicated.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 12, 2021. ; https://doi.org/10.1101/2021.06.09.21258393 doi: medRxiv preprint . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 12, 2021. ; https://doi.org/10.1101/2021.06.09.21258393 doi: medRxiv preprint INTRODUCTION Autism spectrum disorder (autism) is a long-term neurodevelopmental condition with a highly variable clinical manifestation 1 . A substantial proportion of individuals with autism diagnoses present with gastrointestinal symptoms 2,3 , leading to interest in potential associations with gastrointestinal conditions, such as inflammatory bowel disease (IBD). In particular, a link between parental IBD and offspring autism has been hypothesised. There is evidence suggesting that autism and IBD might share risk genes 4 , while characteristics of IBD such as immune system dysregulation, micronutrient malabsorption, and anaemia 5-7 may be perinatal risk factors for autism [8][9][10][11] . However, observational studies have provided inconclusive evidence [12][13][14][15] , and the underlying aetiology of any associations is unclear.
Elucidating the associations between parental IBD and offspring autism may offer mechanistic insights into the origins of autism.
We conducted four complementary studies ( Figure 1) to investigate: (1) associations between parental diagnoses of IBD and offspring autism in a nationwide cohort in Sweden; (2) genetic correlation between IBD and autism using genome-wide association study (GWAS) summary statistics; (3) polygenic associations between maternal genetic liability to IBD and offspring autistic traits in a large UK birth cohort; and (4) potential causal effects of genetic liability to IBD on autism and the possibility of reverse causation using bidirectional two-sample Mendelian randomization (MR).
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Study 1: Investigating associations between parental diagnoses of IBD and offspring autism-Swedish cohort study.
We used individual-level data from 'Psychiatry Sweden', a comprehensive national register linkage, to investigate whether parental IBD diagnosis is associated with offspring autism diagnosis.  Figure S1). The Stockholm Regional Ethical Review Committee (DNR 2010/1185-31/5) approved the study.

All children born in
Offspring autism was identified in the National Patient Register (NPR) using ICD-9 and ICD-10 codes (eMethods S1). Lifetime history of parental IBD, Crohn's disease (Crohn's) and ulcerative colitis (UC) were identified using ICD-9 and ICD-10 codes in the NPR (eMethods S1).
Using STATA/MP15, we estimated the odds ratios and 95% confidence intervals of the association of mother's and father's diagnosis of IBD (any IBD, Crohn's, or UC) with offspring autism using generalised estimating logistic models with robust standard errors accounting for clustering of multiple children born to the same parents.
Model 1 was unadjusted. Model 2 was adjusted for parental age at delivery 17 , migrant status 18 , education level, family income quintile at birth 19 , parents' history of psychiatric diagnosis prior to the birth of the child and offspring sex, birth year and birth order. Model 3 . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted June 12, 2021. ; https://doi.org/10.1101/2021.06.09.21258393 doi: medRxiv preprint was additionally mutually adjusted for maternal and paternal IBD diagnoses to avoid bias from assortative mating 20 . As a sensitivity analysis, we restricted parental IBD diagnoses to those recorded prior to the birth of the index person and investigated associations with offspring autism. Additionally, we investigated associations between any parental IBD diagnoses and offspring autism with and without intellectual disabilities (ID) separately, since these groups may have distinct genetic and environmental risk factors 21-24 and outcomes 25,26 .

Study 2: Investigating genetic correlations-LD-Score regression
We used LD-score regression to estimate the genetic correlation between genetic liability to autism and IBD, Crohn's and UC.
LD-score regression allows the estimation of the genetic correlation between polygenic traits using GWAS summary statistics by capitalising on patterns of linkage disequilibrium among common genetic variants 27 . We used the latest available GWAS summary data on autism (N cases = 18,381; N controls = 27,969) 28 , IBD (N cases = 25,042; N controls = 34,915) 29 , Crohn's (N cases = 12,194; N controls = 28,072) 29 and UC (N cases = 12,366; N controls = 33,609) 29 . Detailed information on study samples and case definition can be found in the original publications.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted June 12, 2021. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted June 12, 2021. ; Broad autism phenotype-autism factor mean score We used a measure of the broad autism phenotype previously estimated in ALSPAC as the mean score of 7 factors derived from a factor analysis of 93 measures related to autism in ALSPAC 33 . The measure was available in 13,128 children and strongly predictive of the autism diagnosis measured independently via school records, record linkage and parental reports 33 . Other autism trait measures or diagnoses were not used as there were fewer genotyped mothers and children with these measures.

Calculation of Polygenic Risk Scores in ALSPAC and statistical analysis
PRS were calculated using PLINK version 1.9, applying the method described by the Psychiatric Genomics Consortium (PGC) 34 . SNPs with mismatching alleles between the discovery and target dataset were removed. The MHC region was removed (25 Mb -34 Mb), except for one SNP representing the strongest signal within the region. Using ALSPAC data as reference panel, SNPs were clumped with an r 2 of 0.25 and a physical distance threshold of 500 kB. The optimal p-value threshold for PRS is dependent on discovery and target sample sizes, as well as SNP inclusion parameters (e.g., r 2 ) 35,36 . For this reason, we calculated PRS for each participant across 13 p-value thresholds (5e-8 to 0.5), standardised by subtracting the mean and dividing by the standard deviation. We defined PRS corresponding to p-value threshold 0.05 as our primary exposure, based on a previous ALSPAC study 37 . This threshold has been found to have sufficient predictive ability for IBD and its subtypes 38 . We could not directly assess the predictive power and optimal p-value threshold of our PRSs in our target sample as there were few UC (n=12) and Crohn's cases (n=16).
After constructing PRS for IBD, UC and Crohn's in mothers and children, we performed linear regressions using STATA/MP 15 to examine associations with the standardised autism factor mean score in childhood. Analyses were adjusted for child's sex and the first 10 principal components of the ALSPAC genotype data to avoid population stratification bias 35 .
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted June 12, 2021.

Study 4: Investigating bidirectional causal links-Two-sample Mendelian randomisation
We performed two-sample Mendelian randomisation (MR) to assess bidirectional causal links between genetic liability to autism and IBD and its subtypes, and vice versa.
MR can be implemented as an instrumental variable approach, utilising common genetic variants as instruments for exposures of interest, allowing assessment of causal effects and their direction on outcomes. MR relies on the following assumptions : (i) there must be a robust association between the common genetic variants and the exposure, (ii) the variants should operate on the outcome entirely via the exposure, (iii) the variants should not be associated with any confounders of the associations between the exposure and the outcome 39 . In this study, we applied two-sample MR, in which the effect sizes and standard errors of the instruments for the exposure and the outcome were extracted from separate GWASs conducted in independent samples from the same underlying population 40 .

Genetic Instruments
Genetic instruments were extracted from the overlapping set of SNPs between the autism 28 , IBD 29 , UC 29 , and Crohn's 29 GWASs. This ensured that all selected genetic instruments would be present in the outcome GWAS.
GWAS summary data were restricted to a common set of SNPs and then clumped in PLINK 1.90 using the 1000Genomes 30 phase 3 European ancestry reference panel, and an r2= 0.01, within a 10,000 kb window. Among the independent variants, instruments were defined using a genome-wide significance threshold of p≤5×10 -08 . The only exception was autism, as only two independent and genome-wide significant variants were identified. We therefore relaxed the p-value threshold to 5×10 -07 to improve statistical power, as used previously 41 .
Supplementary Figure S1 illustrates the process of instrument definition, and supplementary table S1 contains information on the genetic instruments used.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Sensitivity Analyses to test robustness of causal effect estimates
We assessed the strength of the instruments by estimating the mean F statistic. As a rule of thumb, the IVW is unlikely to suffer from weak instrument bias if mean F>10 44 .
We assessed the consistency of the IVW causal effect estimates using sensitivity analyses, Two-sample MR analyses were performed using the TwoSampleMR R package 43 in R version 3.5.1.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Study 1: Associations between parental IBD diagnoses and offspring autism
Cohort characteristics by exposure to parental IBD diagnoses are presented in supplementary tables S7 and S8.
Maternal IBD diagnosis was associated with offspring autism in crude and adjusted models (Any IBD diagnosis: OR MODEL3 = 1.32; 95% CIs: 1.25 to 1.40; Table 1). Similar results were observed in analyses of maternal UC and Crohn's diagnoses and offspring autism (Table 1), and in analyses restricted to maternal IBD diagnoses prior to the index person's birth (Any IBD diagnosis: OR MODEL3 = 1.20; 95% CIs: 1.09 to 1.32; Supplementary Table S9). The paternal IBD associations with autism were weaker (OR MODEL3 = 1.09; 95% CIs 1.02 to 1.17) than the maternal associations (Table 1). Point estimates for associations of parental IBD diagnoses to autism without ID were higher than those for autism with ID, although confidence intervals overlapped (Supplementary Table S10).

Study 2: Genetic correlation between IBD and autism
We found little evidence of a genetic correlation between genetic liability to autism and IBD, UC, or Crohn's ( Table 2). Heritability scores, chi-squares and intercepts satisfied the conditions to provide reliable LD-score regression estimates 27,48 (Supplementary Table S11).

Study 3: Associations between polygenic risk for IBD, UC, Crohn's and broad autism phenotype in ALSPAC
The study sample derivation and characteristics can be found in Supplementary Figure S4.
Maternal polygenic risk for IBD, UC, Crohn's and offspring broad autism phenotype Maternal polygenic risk for UC and Crohn's was associated with a higher autism factor mean score in the offspring (UC: 95%CIs: 0.01 to 0.05; p= 0.004). Similar results were found across other p-value thresholds (0.5-0.05). The effect size of the association between maternal polygenic risk for IBD and autism factor mean score, was comparable to that of UC and Crohn's, although confidence . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

Study 4: Causal effects of genetic liability to IBD on risk of autism
The mean F statistics of the IBD, UC and Crohn's instruments were 67, 68 and 70, respectively, suggesting adequate strength. There was evidence of a causal effect of genetic liability to UC on risk of autism ( IVW OR= 1.04; 95% CIs: 1.01 to 1.07; p= 0.006). Evidence for the effect of genetic liability to IBD and Crohn's on autism risk was weaker, although the magnitude and direction of the effect estimates was comparable to the UC results ( Table 4).
The magnitude and direction of causal effect estimates was consistent across all sensitivity analyses, and there was no evidence to suggest the influence of horizontal pleiotropy (Supplementary Table S14). Results of analyses with instruments extracted from the autism GWAS excluding ID cases were comparable to our primary effect estimates (Supplementary   Table S15).

Causal effects of genetic liability to autism on risk of IBD
The mean F statistic of the autism instruments was 28, suggesting adequate strength. There was no evidence of a causal effect of genetic liability to autism on risk of IBD, UC or Crohn's ( Table 4). The estimates were consistent across sensitivity analyses, with overlapping confidence intervals, and were unlikely to be influenced by horizontal pleiotropy (Supplementary Table S16). Repeating our analyses with instruments extracted from the autism GWAS excluding all ID cases yielded similar results (Supplementary Table S17).
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 12, 2021. ; https://doi.org/10.1101/2021.06.09.21258393 doi: medRxiv preprint DISCUSSION We used four complementary approaches to investigate the associations between parental diagnoses and genetic liability to IBD and offspring autism. Conducting a nationwide register-based cohort study in Sweden we found evidence of associations between parental diagnoses of IBD and offspring autism. Importantly, the maternal effect sizes were larger than paternal, without overlapping confidence intervals. PRS analyses in the ALPSAC birth cohort suggested associations between maternal genetic liability to IBD and offspring autism, while two-sample MR studies provided evidence of a potentially causal effects of genetic liability to IBD on autism risk. There was little evidence to suggest a genetic correlation between autism and IBD, as indicated by LD-score regression analyses.
The use of different study designs to triangulate our findings is a notable strength of our approach 49 . The Swedish nationwide register-based cohort study of over 2 million parentchild pairs is the largest to date on parental IBD-offspring autism. The ALSPAC cohort is a unique resource to assess polygenic associations through rich intergenerational genetic and phenotypic data. Finally, in the MR analyses we used the largest GWAS data available for all conditions and conducted several sensitivity analyses to test the robustness of our findings.
Considering study limitations, in the Swedish registers the possibility of measurement error in IBD diagnoses cannot be excluded. However, this is likely to be non-differential in relation to our study outcome and would therefore bias our findings towards the null. Secondly, while PRSs were based on large GWAS samples, it was not possible for us to investigate the variance explained by the PRSs in our target sample and the variance explained in the phenotype might be relatively small. Thirdly, in two-sample MR analyses investigating the effects of genetic liability to autism on risk of IBD, we used a relaxed instrument inclusion pvalue threshold, this could potentially result in including weak instruments and therefore bias the causal effect estimates. The F statistic of the autism instruments in our analyses suggested that weak instrument bias is unlikely. Finally, using GWAS data we could only . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 12, 2021. ; 1 6 investigate the possible contribution of common variants acting under an additive model and not any contribution from rare variation which has been found to be implicated in autism 51,52 .
Larger samples and GWAS, as they become available, may help replicate and build on our findings and their potential explanations. Specifically, under liability-threshold models of inheritance, parental genetic liability to IBD could reflect sub-phenotypic manifestations of the condition 55-58 and potentially influence fetal development in utero. These subclinical manifestations could include immune system dysregulation, micronutrient malabsorption and anaemia ,6,7,53,54 that have been found to be implicated in IBD 5-7 and autism 8-11 .
In conclusion, triangulating evidence from a nationwide register-based cohort study, genetic correlation, polygenic risk score analyses and MR, we found evidence suggesting associations between parental, particularly maternal, diagnoses of IBD and offspring autism.
Links between maternal genetic liability to IBD and offspring autism may reflect the influence of the maternal genotype on the prenatal/intrauterine environment. Investigating the mechanisms behind these findings may provide valuable insights into the origins of autism.
. CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 12, 2021. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 12, 2021. is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 12, 2021. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 12, 2021. ; https://doi.org/10.1101/2021.06.09.21258393 doi: medRxiv preprint . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 12, 2021. .

CC-BY 4.0 International license
It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 12, 2021. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 12, 2021. is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity.

(which was not certified by peer review)
The copyright holder for this preprint this version posted June 12, 2021. . CC-BY 4.0 International license It is made available under a is the author/funder, who has granted medRxiv a license to display the preprint in perpetuity. (which was not certified by peer review) The copyright holder for this preprint this version posted June 12, 2021. ; https://doi.org/10.1101/2021.06.09.21258393 doi: medRxiv preprint