Early-life exposure to perfluorinated alkyl substances modulates lipid metabolism in progression to celiac disease

Lisanna Sinisalu, Partho Sen, Samira Salihović, Suvi M. Virtanen, Heikki Hyöty, Jorma Ilonen, Jorma Toppari, Riitta Veijola, Matej Orešič, Mikael Knip, Tuulia Hyötyläinen
doi: https://doi.org/10.1101/2020.04.02.20051359

Abstract

OBJECTIVES Celiac disease (CD) is a systemic immune-mediated disorder with increased frequency in the developed countries over the last decades implicating the potential causal role of various environmental triggers in addition to gluten. Herein, we apply determination of perfluorinated alkyl substances (PFAS) and combine the results with the determination of bile acids (BAs) and molecular lipids, with the aim to elucidate the impact of prenatal exposure on risk of progression to CD in a prospective series of children prior the first exposure to gluten (at birth and at 3 months of age).

METHODS We analyzed PFAS, BAs and lipidomic profiles in 76 plasma samples at birth and at 3 months of age in the Type 1 Diabetes Prediction and Prevention (DIPP) study (n=17 progressors to CD, n=16 healthy controls, HCs).

RESULTS Plasma PFAS levels showed a significant inverse association with the age of CD diagnosis in infants who later progressed to the disease. Associations between BAs and triacylglycerols (TGs) showed different patterns already at birth in CD progressors, indicative of different absorption of lipids in these infants.

DISCUSSION PFAS exposure may modulate lipid and BA metabolism, and the impact is different in the infants who develop CD later in life, in comparison to HCs. The results indicate more efficient uptake of PFAS in such infants. Higher PFAS exposure during prenatal and early life may accelerate the progression to CD in the genetically predisposed children.

WHAT IS KNOWN Several observational studies have implicated a role of early life environmental triggers other than gluten in the development of CD. This is supported by the findings showing dysregulation of lipids already prior to the first introduction of gluten.

WHAT IS NEW HERE We show that prenatal exposure to perfluorinated compounds is associated with changes in the lipid metabolism, most likely through the bile acids, and that a high exposure during prenatal and early life may accelerate the progression to CD in the genetically predisposed children.

TRANSLATIONAL IMPACT Exposure to environmental chemicals may impact the rate of progression to CD and should be assessed as a potential risk factor of CD in larger clinical cohort settings.

Introduction

Celiac disease (CD) is a systemic immune-mediated disorder, which is triggered by gluten and other prolamins in genetically susceptible individuals1. The frequency of CD has increased in the developed countries over the last decades, implicating the potential causal role of various environmental triggers in addition to gluten. Similarly to CD, the rate of other autoimmune diseases such as type 1 diabetes (T1D) has also increased during the last decades2,3. CD and T1D share common predisposing alleles in the class II HLA-region4,5, and approximately 10% of patients with T1D also develop clinical CD6,7, while subjects with CD are at-risk for developing T1D before 20 years of age8. The possible triggers that have been indicated to affect the progression to CD include the composition of the intestinal microbiota, infant feeding, and the use of antibiotics1,9,10.

The hygiene hypothesis, stating that a decrease of the infectious burden is associated with the rise of allergic and autoimmune diseases, has also been proposed in CD because CD has shown to be more common in developed countries11-13. Another possible explanation for varying incidence in different populations implicates the role of different infant feeding patterns (including the amount and timing of gluten introduction) in families with low socioeconomic status12. Socioeconomical status (SES) may also have a broader role in increasing the risk for CD, such as different exposures to environmental pollutants. Higher SOS has in several studies been linked to higher burdens of several persistent organic pollutants (POP), particularly to perfluorinated alkyl substances (PFAS)14,15. Exposure to POPs has, however, been poorly investigated as a risk factor for CD. In T1D, recent studies identified associations between the exposures to environmental chemicals and the increased risk of islet autoimmunity16-19. We have shown that fetal exposure to PFAS increases the risk of T1D in genetically susceptible individuals20, and further, that the metabolic profiles of the infants are affected by the prenatal exposure to PFAS. In particular, the PFAS exposure has a marked impact on bile acid (BA) profiles, which, in turn, were linked to changes in specific phospholipids, e.g., lysophosphatodylcholines (LPCs), phosphatodylcholines (PCs) and sphingomyelins (SMs). Similarly, we have recently identified systematic differences in plasma lipid profiles between children who later progressed to clinical CD during the follow-up, when compared to children who remained healthy5. Importantly, these differences were observed prior to the first introduction of gluten in the diet and before the first signs of CD-associated autoimmunity. Similar results were recently reported in an Italian cohort study21.

In line with the observations of intestinal dysbiosis in CD1, changes in lipid metabolism5,21 as well as the documented link of CD with liver disorders22, it has been shown that circulating BAs are elevated in CD, including in children23,24. BAs not only facilitate the digestion and absorption of dietary lipids in the small intestine, they are also important metabolic regulators involved in the maintenance of lipid and glucose homeostasis25. BAs are produced in the liver, and their homeostasis is maintained through the tightly controlled enterohepatic circulation. Moreover, there is a close interplay between BA and gut microbiota. Gut microbiota is involved in the biotransformation of secondary BAs, while BAs can modulate the microbial composition in the gut due to their antimicrobial activity as well as through regulation of the composition of the intestinal microbiota through farnesoid X receptor (FXR) and Takeda G-protein receptor 5 (TGR-5)26,27. PFAS, on the other hand, have been shown to undergo similar enterohepatic circulation as the BAs, and PFAS can also suppress the BA biosynthesis in the liver28.

Here we hypothesized that early life PFAS exposure affects the BA profiles in children and explains the association between the early changes in lipid profiles in those children who later progress to clinical CD. We investigated the impact of early-life exposure to PFAS in children who progressed to clinical CD during their follow-up, and compared the exposures to those observed in children who remained healthy during the follow-up. In addition, we also analyzed BAs in the same samples, and performed integrative analysis of previously published lipidomic profiles5, together with PFAS and BA levels.

Methods

Study design

The current study is a part of the Type 1 Diabetes Prediction and Prevention study in Finland (DIPP), which is an ongoing prospective birth cohort study initiated in 1994. In DIPP, parents of newborn children at the university hospitals of Turku, Tampere and Oulu in Finland are asked for permission to screening for T1D-conferring HLA risk alleles in the umbilical cord blood. Families of children identified as having an increased HLA-conferred risk for T1D were invited to join the follow-up study. This study included children born in Tampere University Hospital between August 1999 and September 2005. During that period, 23,839 children were screened at birth for increased risk of T1D, and 2,642 eligible children were enrolled in the follow-up and had at least two visits to the study center. These children carry the high-risk HLA DQB1*02/*03:02 genotype or the moderate-risk HLA-DQB1*03:02/x genotype (x ≠ DQB1*02, 03:01, or 06:02). More than 1,200 of these children took part in the DIPP-CD study. Because of the enrolment criteria 60% of the children were male. At each visit, the families were interviewed for dietary changes, infections, growth, important family related issues and the children gave a non-fasting venous blood sample. The children were followed for four CD-related antibodies: anti-transglutaminase 2 (anti-TG2), anti-endomysium (EMA), antigliadin (AgA-IgG and AgA-IgA) and anti-reticulin (ARA) antibodies, and for four T1D-associated autoantibodies: islet cell autoantibodies (ICA), autoantibodies against insulin (IAA), tyrosine phosphatase-like protein (IA-2A) and glutamate decarboxylase (GADA). IgA deficiency was excluded. All children participating in this study were of Caucasian origin. None of the mothers had CD. The children in the CD follow-up cohort were annually screened for anti-TG2 antibodies using a commercial kit (Celikey Pharmacia Diagnostics, Freiburg, Germany). If a child’s sample was found positive, all previous and forthcoming samples were analyzed for the entire set of CD-related antibodies. A duodenal biopsy was recommended for all anti-TG2-positive children. If the biopsy was consistent with the ESPGHAN criteria of 1990, a gluten-free diet (GFD) was recommended. None of the children participating in this study had any T1D-associated autoantibodies in any samples during the follow-up. The first exposure to gluten in this study was at the median age of 5 months (Table 1).

View this table:
Table 1.

Demographic and clinical characteristics of the study subjects.

We have randomly selected 17 children with biopsy proven CD (progressors) and a matched control for each progressor, carrying similar risk HLA alleles, born within ±1 month of each other, having had each sample taken within ±1 month of each other and living in the same region of the country throughout the whole follow-up period. Samples at birth (cord blood) and at 3 months of age were analyzed. Altogether, 76 plasma samples from children developing CD and from their HCs were analyzed. The clinical and genetic data of the participants are found in Table 1.

The ethics committee of the Tampere University Hospital approved the study. Written informed consent was obtained from the parents for HLA-screening, autoantibody analysis and intestinal biopsies.

Analysis of PFAS and BAs

The BAs and PFAS were analyzed using the established method as described previously29. In brief, the sample preparation was done with 25 mg Ostro Protein Precipitation and Phospholipid Removal 96-well plate (Waters Corporation, Milford, USA), using 30 μL of serum and a set of BA and PFAS internal standards. Matrix-matched calibration standards were made using newborn bovine serum. Analyses were performed on an Acquity UPLC system coupled to a triple quadrupole mass spectrometer (Waters Corporation, Milford, USA) with an atmospheric electrospray interface operating in negative ion mode. Aliquots of 10 μL of samples were injected into the Acquity UPLC BEH C18 2.1 mm × 100 mm, 1.7 μm column (Waters Corporation). A trap column (PFC Isolator column, Waters Corporation) was installed between the pump and injector and used to retain fluorinated compounds originating from the HPLC system and the mobile phase. The eluent system consisted of (A) 2 mM NH4Ac in 70% MiliQ: 30% methanol and (B) 2 mM NH4Ac in 100% methanol. The gradient was programmed as follows: 0–1 min, 1 % solvent B; 1–13 min, 100 % solvent B; 13-16 min, 100 % solvent B; 16-17 min, 1 % solvent B, flow rate 0.3 mL/min. The total run time for UPLC-MS/MS analysis was 17 minutes, while the total run time for each sample injection was 20 minutes, including the reconditioning of the analytical column. MS analysis was performed in multiple reaction monitoring (MRM) mode and experimental details of the MS/MS method are given in Supplementary Table S1.

View this table:
Supplementary Table S1.

Acquisition parameters including the list or targets compounds ordered by retention time.

Lipidomic analyses

The lipidomic analyses were performed as described previously5. The plasma samples (10 μL) were extracted using a modified version of the previously published Folch procedure30. The samples were analyzed using an ultra-high-performance liquid chromatography quadrupole time-of-flight mass spectrometry method (UHPLC-Q-TOF-MS from Agilent Technologies (Santa Clara, CA, USA). The analysis was carried out on an ACQUITY UPLC® BEH C18 column (2.1 mm × 100 mm, particle size 1.7 μm) by Waters (Milford, USA). Internal standard mixture was used for normalization and lipid-class specific calibration was used for quantitation as previously described. Quality control was performed throughout the dataset by including blanks, pure standard samples, extracted standard samples and control plasma samples. Relative standard deviations (%RSDs) for lipid standards representing each lipid class in the control plasma samples (n=8) and in the pooled serum samples (n = 20) were on average 11.7% (raw variation). The lipid concentrations in the pooled control samples was on average 8.4% and 11.4% in the standard samples. This shows that the method is reliable and repeatable throughout the sample set. MS data processing was performed using open source software MZmine 2.1831.

Statistical analyses

Statistical analysis was performed using R (v3.6.0) statistical programming language32. Previously published lipidomics dataset5 was obtained from MetaboLights with the study identifier MTBLS729. The data was log-transformed and missing values were imputed by half of the row’s minimum. The datasets were auto-scaled prior to multivariate analyses.

To integrate across different data types, we applied sparse generalized correlation discriminant analysis via the DIABLO framework, part of mixOmics package (v6.10.8)33. The method constructs components across different data blocks, by maximizing their covariance with each other and a given response (Y) variable. Heterogenous data such as PFAS, BAs and lipid levels were partitioned into three different blocks and regressed to a binary response variable, i.e., CD progressors (CD) or HCs. Regularized sparse partial least squares discriminant (sPLS-DA)33,34 models were fitted. The optimal number of components that achieve the best performance based on the overall error rate or Balanced Error Rate (BER) were determined. Block sPLSDA models were developed at two time-points, i.e., at birth (cord blood) and at 3 months of age. Moreover, these models were cross-validated35 by 5-fold cross-validation (CV) with (N =100 repeats). The final model performances were assessed by area under the curve (AUC), overall misclassification error rate, and BER generated using ‘perf’,’auroc’,’block.splsda’ functions.

The key predictors/contributors that are jointly associated with the response variable of interest across all the input data matrices were identified by their Variable Importance (VI) score, i.e (absolute weighted loadings). The direction (up or down) of contribution of a particular predictor was determined by estimating the mean relative log fold change in the CD progressors vs. HCs at a particular time-point.

In order to understand which variables influence the outcomes, multivariate correlations were performed by a method described in33. In addition, bivariate correlation was performed using Pearson correlation. Debiased Sparse Partial Correlation algorithm (DSPC) was used for estimating partial correlation networks, visualized by the MetaboAnalyst 436 version with cutoff values of correlations between +/- 0.22 to 0.75. Univariate analysis (unpaired two sample t-test and paired t-test) using the ‘t.test’ function, was deployed to identify mean differences in the concentration of individual PFAS and BAs between CD progressors vs. HCs at birth and at 3 months of age. Libraries/packages such as ‘Heatmap.2’, ‘mixOmics’, ‘boxplot’, ‘beanplot’, ‘gplot’ and ‘ggplot2’ were used for data visualization.

Results

Levels of PFAS and BAs in the infants

Seven PFAS compounds were detected in the samples, namely PFHpA, PFHxS, PFOA, PFNA, PFOS, PFDA and PFUnDA, both at birth and at 3 months of age (Supplementary Table S2). PFOS and PFOA had the highest concentrations and they were detected in the majority of samples, whilst PFDA and PFNA were detected in less than 10% of the samples and were excluded from further analyses. The levels of PFAS were lower in the cord blood (CB) than at the age of 3 months (Figure 1A-C). No significant differences were observed between cases and controls either at birth or at 3 months of age (Figure 1A-C). Interestingly, the levels of the total PFAS were significantly elevated from birth to 3 months of age only in the CD group (fold change (FC) 3 months vs. at birth = 1.72, p = 0.003 Supplementary Table S3) although also controls had an increasing trend (FC = 1.02, p = 0.35). Among the individual PFAS, PFOA showed nearly twofold increase, being significant in both groups (Figure 1A-C) (Supplementary Table S3).

View this table:
Supplementary Table S2.

PFAS concentrations.

View this table:
Supplementary Table S3.

Change of PFAS levels from cord blood to 3 months.

Figure 1.
Figure 1. Bean plots showing distribution of selected PFAS and BAs in CD progressors and HCs.

A-C) Levels (log2 intensities) of selected PFAS in cord blood (0 Mo) and at 3 months (3 Mo) of age in the CD progressors (CD) and HCs. D-F) Levels (log2 intensities) of selected BAs in the cord blood and at 3 months of age. The red dotted line denotes the mean of the population. The black solid lines in the bean plots represent the group mean. The group mean difference between case-control at a particular time-point was tested by unpaired two-sample t test. Mean difference between cord blood and 3 months samples obtain from the same infants were tested by paired t test. CD, HCs and BA denotes CD progressors, healthy controls and bile acids respectively.

In addition to the two primary BAs, we measured five primary/conjugated and secondary BAs (Supplementary Table S2). The levels of secondary BAs were mostly below the detection limits, particularly in the cord blood samples, where the total BA pool consisted mainly of GCA and GCDCA. Similarly, as for PFAS, no significant differences were observed between cases and controls either at birth or at 3 months of age (Figure 1D-F). The BA profiles at birth and at 3 months of age were markedly different in both groups (Figure 1D-F). The total BA pool showed significant, almost sevenfold upregulation (p = 0.0001). Particularly, GHCA and CDCA, which were close the limits of detection at birth, have increased substantially with age.

Associations between PFAS, BA and lipid levels

Partial correlation network analysis between PFAS, BAs and lipids (classes) measured incord blood showed a different pattern of association in HCs (Figure 2A) and CD progressors (Figure 2B). In the HCs, the level of deoxycholic acid (DCA) was inversely asociated with the sex of the infant, whereas no such association was seen in the CD progressors. Interestingly, DCA concentration was inversely associated with the level of PFOA. The CD group showed a strong association between age at diagnosis and PFOA exposure at both time points (Figure 2B, Supplementary Figure 1). In general, there was direct association between PFAS and several lipid classes. Age at diagnosis was also inversely associated with total TGs at birth. The association showed a clearly different pattern in the controls, both at birth and at the later time point (Figure 2).

Supplementary Figure 1.
Supplementary Figure 1.

(A) correlation plots of lipid-BA and lipid-PFAS (3 months), Pearson correlations, significant correlations marked with * (p<0.0001 ****, p<0.001 ***, p<0.05 **, p<0.1*) and (B) Partial correlation network of the previously identified panel of CD-related lipids, BAs and PFAS (at birth and at 3 months of age).

Significant correlation was observed in the CD group between BAs and the lipids. Particularly, GCDCA, CDCA and DCA showed significant positive correlations with TGs of low carbon number and double bond count, while TCDCA and GUDCA showed inverse associations with PCs at 3 months of age. At the same time point, PFOA, PFOS and PFUnDA, and total PFAS showed positive correlations with the TGs containing low carbon number and double bond count and inverse associations with PCs. In the partial correlation network, PFOA and PFHxS showed an inverse association with the age at diagnosis and PFHxS further positive associations with CE(18:2) and PC(34:3) (Figure 4B). PFHpA showed an inverse association with one TG and one PC species and a positive association with GHCA. BAs showed mainly inverse associations with several lipid species.

Figure 2.
Figure 2. Molecular network of PFAS (yellow color), BAs (green color), lipid classes (grey color), and clinical variables (orange color) measured in the cord blood:

A) Healthy controls. B) CD progressors. Red and blue lines denote positive and negative correlations respectively. The width of the lines represent the strength of correlation.

To understand how the individual PFAS and BAs are associated with the molecular lipid species, we performed multivariate33 and bivariate correlation analysis. Bivariate correlation analysis identified different correlation patterns in HCs and in CD progressors (Figure 3). In the control group there were significant associations mainly between the levels of the lipid classes at birth, and only the total concentrations of TGs showed associations with BAs (Figure 3A). In the CD progressor group, we observed significant correlations between the levels of BAs (GUDCA, GHCA, CDCA) and with multiple lipid classes in cord blood (PE, TG, SM) (Figure 3B). At 3 months of age, the patterns were different from those seen at birth. In the control group, specific BAs (CDCA, GUDCA, TCDCA and GCHA) showed a significant association with lipid levels, mainly with cholesterol esters (CEs) and TGs, with variable patterns depending on the fatty acyl saturation level (Figure 3C). Only a weak correlation was observed between the levels of PFAS and LPCs, while overall the PFAS levels were not associated with either lipids or BAs. In CD progressors, the patterns were different from the controls also at 3 months of age (Figure 3D). In this group, PFAS (PFOS and PFUnDA) levels and concentrations of TGs showed significant correlations, particularly with TGs containing saturated FAs. In addition, concentrations of BAs were also associated with those of TGs (DCA, GCA, CDCA and TCDA, total BAs), the directions being dependent on the bile acid species.

Figure 3.
Figure 3. Correlation plots of PFAS, BA and lipid classes

in (A) Controls at birth, (B) CD at birth, (C) controls at 3 months and (D) CD at three months. Red and blue colors denote positive and negative correlations respectively. Significant Pearson correlations are marked with *.

The main difference between the groups was the absence of correlation between the LPCs and the other phospholipid classes (PC, SM, PC) in CD progressors, while there was a strong and significant co-regulation of these lipid classes in the control group. At birth, there were no significantly different individual lipids between the groups5, however, the LPCs showed overall increased levels and most of the PCs decreased levels. LPC/PC ratio, with selected PC species containing PUFAs, showed significant upregulation in the CD group (FC=1.3, p=0.01).

Molecular Predictors of PFAS, BA and Lipids in CD progression

Multivariate correlation analysis showed strong associations (r > 0.8) among PFAS, BAs and molecular lipids such as cholesterol esters (CEs), lysophosphatidylcholines (LPCs), phosphatidylcholines (PCs), phosphatidylethanolamines (PEs), sphingomyelin (SMs) and triacylglycerols (TGs) (Figure 4A). These associations become more complex at 3 months of age (Figure 4D).

Figure 4.
Figure 4. Integrative multiblock component analysis of PFAS, BAs and molecular lipids in the CD progressors vs. HCs.

A, D) Circos plots (2 principal components) showing overall association (full-design matrix) between covariates of PFAS, BAs and molecular lipids in cord blood(CB), and at 3 months of age. Red and blue line denotes positive and negative correlations respectively. Sparse generalized correlation33 (coefficient, r > 0.5) are shown. B-C) Horizontal barplots showing top selected predictors/contributors/discriminators from each PFAS and BA data blocks, where the length of bar corresponds to the absolute loading weights (variable importance) of the feature, that tends to discriminate CD progressors vs. HCs at birth. The direction of regulation of the predictors are shown with light brown and blue colors that corresponds to the CD progressors or HCs respectively; wherever the mean intensity of the predictor is maximum. E-F) Key predictors of PFAS and BAs, that helped to classify CD progressors vs. HCs at 3 months (3 Mo) of age. The model performance on each block of data are given by cumulative (principal components) area under the curves (AUCs).

Next, we selected those lipids that showed significant differences between CD and HCs in infancy in our previously reported study5. These lipids together with the BAs and PFAS were subjected to multiblock (MB) analysis. MB analysis identified PFAS, BA and molecular lipid species as predictors (discriminative features) that helped to classify CD progression vs. HCs. The association between top predictors (high variable importance scores) are shown (Figure 4A, D).

In cord blood, the levels of PFOA and PFOS were identified as the top linear predictors of CD progression (Figure 4B). On the other hand, among the BAs, the levels of primary BAs such as GCA and GHCA are the strongest predictors of CD progression (Figure 3C). Moreover, elevated concentration of PFOA in the CD progressors was positively associated with the levels of GCA and GHCA. This suggests that prenatal exposure to PFOA and/or PFOS might accelerate the risk of CD progression. Presumably, it alters the BA and lipid metabolism. Interestingly, lipids such as CE(18:1) (intermediate of BA synthesis) were identified as predictors of CD progression; the levels of PFOA and PFOS were inversely correlated with those of CEs (Figures 4A and 5A). In addition, many ether-linked PCs such as PC(O-34:3), PC(O-36:3,4) and PC(O-38:4,5,6) were identified as key contributors, that were downregulated in the CD progressors as compared to the HCs (Figures 4A and 5A).

At 3 months of age, concentrations of PFHxS was marked as a top predictor associated with CD progression, along with levels of PFOA and PFOS. Moreover, the level of PFHxS were increase in the CD progressors as compared to HCs (Figures 1A, 4E and 5B). At this age, elevated levels of secondary BAs such as DCA and CDCA were marked as the key predictors for CD progression (Figure 4F). As stated before, the level of PFOA was inversely correlated with that of DCA (Figure 2B), suggesting an interaction between the PFAS and BAs in the CD progressors. Moreover, PFOA was inversely linked to age at diagnosis of CD (Figure 2B and Figure 6). In addition, the MB model identified several PCs, SM(d36:2), PC(O-38:4,5) as key predictors that discriminate the lipid profiles of CD progressors from those seen in the HCs (Figure 5B).

Figure 5.
Figure 5.

Key lipid contributors/predictors, that jointly (together with PFAS and BAs), aided in the separation of CD progressors vs. HCs, at birth and at 3 months of age. The direction of regulation of the predictors are shown with light brown and blue colors that corresponds to the CD progressors or HCs respectively; wherever the mean intensity of the predictor is maximum.

Figure 6.
Figure 6.

Correlation of PFOA concentration at birth and the age of diagnosis. Pearson correlation r=-0.575.

Discussion

We observed dysregulation between intralipid class regulation at birth in those infants that progressed to CD, particularly in LPC/PC ratios. PFOS and PFHxS levels were also positively associated with LPCs in the CD progressors, but not in HCs. However, 3 months later, the ratio was similar in both groups. Increased LPCs have been indicated to be characteristic for CD progression at later time points, both our earlier study as well as in another recent study showed increased LPCs already in those 4-month-old infants who progressed later to CD21. Furthermore, our results implied that the associations between BAs and complex lipids show different patterns already at birth in those infants who later developed CD as compared with healthy controls. Significant correlations were observed between TGs and specific BAs already at birth in the CD progressors. At 3 months of age, both groups showed a significant correlation between the lipids and BAs, although the patterns were different. At this later time point, the CD group showed significant positive correlations between different classes of TGs while the healthy controls showed inverse correlations. Thus, the data agrees well with our earlier study, in which we observed distinct lipid changes, particularly in the TG class, in children who later progressed to CD already prior to the first exposure to dietary gluten. The data suggested that the specific TGs, found elevated in the CD progressors, may be due to a host response to compromised intake of essential lipids in the small intestine, requiring de novo lipogenesis. The current study suggests that the changes in TG metabolism are also related to alteration in BA metabolism, and that early-life exposure to PFAS may contribute to these changes. Currently, there is very limited amount of data on BA and lipid metabolism in infants. The BA pool in the fetus and in newborn infants is unique. More than 90% of fetal BAs represent conjugated forms of the primary BAs, cholic acid (CA) and chenodeoxycholic acid (CDCA), since the intestinal bacteria necessary to transform primary BAs into secondary BAs are thought to be absent in utero, Based on recent findings, the fetal BAs may be very different than later in life37,38. The enterohepatic circulation of BAs in utero is also minimal. Instead, there is a transplacental gradient for BAs in the fetal-to-mother direction, except for secondary and tertiary BAs, which are more abundant in maternal serum39.

PFAS exposure has been shown to impact BA metabolism through several mechanisms. BAs and PFAS have been shown to have a similar enterohepatic circulation25,28, and it has been indicated that 7-alpha-hydroxylase (CYP7A1), which catalyzes the first and rate-limiting step in the classical pathway in the formation of BAs from cholesterol, may be down-regulated by PFAS25,40. This might lead to increased re-uptake of BAs, which would generate negative feedback loops via the farnesyl-X-receptor and subsequently reduce their de novo synthesis25,41. It has also been shown that PFOA inhibits the function of the hepatocyte nuclear factor 4α42, which plays a central role in the regulation of BA metabolism in the liver, and is linked both with the synthesis and conjugation of primary BAs. However, in neonates, the alternative (acidic) pathway is the major pathway for bile acid synthesis and this pathway is governed by mitochondrial sterol 27-hydroxylase (CYP27A1), which can initiate a process independent of CYP7A1. Only after weaning, CYP7A1 is expressed and the classic pathway becomes the major pathway for bile acid synthesis in the liver25. There is currently no data on the impact of PFAS exposure on the regulation of the CYP27A1. On the other hand, the impact of PFAS on the TG metabolism may be modulated through the bile acids, as the bile acid receptor FXR has also a regulatory role in triglyceride metabolism43-45. In humans, recent data suggests that FXR is not activated directly by PFAS43. In the liver, FXR activation, by, e.g., bile acids, would result in downregulation of CYP7A1, which in addition to inhibition of the classical BA synthetic pathway also reduces the expression of several genes mediating free fatty acid synthesis, thereby attenuating de novo lipogenesis25. Thus, activation of FXR modulates free FA oxidation and the clearance of TG to the circulation.

The acylalkyl PCs type of lipids that showed significant association with the CD progression were positively associated with the levels of total BAs and PFHxS, indicating potentially disturbed hepatic synthesis of these lipids, particularly as the intra-lipid regulation between these lipids and particularly TGs was disturbed in CD progressors. These ether lipids are perceived to function as endogenous antioxidants in addition to their structural roles in cell membranes, and emerging studies suggest that they are involved in cell differentiation and signaling pathways as well46. Interestingly, these lipids have been observed to be endogenous antigens to activate invariant natural killer T cells (iNKT)47, which are subset of innate immune cells. Recently, therapeutic potential of iNKT cell antigens against autoimmunity have been implicated48. Thus, the reduced levels of the alkyl ether lipids could suggest a compromised response to oxidative stress. Interestingly, these lipids have previously been found to be associated with disease development also in other autoimmune diseases, such as type 1 diabetes49, thus indicating that they may play an important role in the development of autoimmune disorders. Due to their role in cell membranes, the decrease in ether-linked PCs particularly, PC(O-38:4,5) might also attribute to compromised intestinal permeability in the CD progressors.

Taken together, our results show that PFAS exposure may modulate lipid and BA metabolism, and the impact is different in the infants who develop CD later in life, when compared with healthy controls. Although we did not observe any significant differences in the levels of PFAS exposure in those children that later developed CD, which may be due to the small sample size, we did observe a significant increase in their levels in progressors to CD from the time of birth to 3 months of age, suggesting more efficient uptake of PFAS. Furthermore, age at diagnosis of CD was strongly associated with the PFAS exposure. Our study thus suggests that further investigations of the impact of exposures to environmental chemicals in the development of autoimmune diseases are merited.

Data Availability

The lipidomics data sets and the clinical metadata generated in this study are available from MetaboLights, under accession number MTBLS729.

Conflicts of interest

Guarantor of the article

Tuulia Hyötyläinen, Ph.D.

Specific author contributions

T.H., M.K., M.O.: conception and design of the study; L.S., S.S., S.V., H.H., J.I., J.T., R.V., M.O., M.K., T.H.: generation, collection, assembly, analysis, and/or interpretation of data; T.H. and M.O.: drafting or revision of the manuscript; all authors: approval of the final version of the manuscript.

Financial support

This study was supported by funding from VetenskapsrÅdet (to T.H.; grant no. 2016-05176) and Formas (to T.H. and M.O.; grant no. 2019-00869). The DIPP study was supported by the following grants: JDRF (grants 1-SRA-2016-342-M-R, 1SRA-2019-732-M-B); European Union (grant BMH4-CT98-3314); Novo Nordisk Foundation; Academy of Finland (Decision No 292538 and Centre of Excellence in Molecular Systems Immunology and Physiology Research 2012-2017, Decision No. 250114) and Special Research Funds for University Hospitals in Finland.

Potential competing interests

None to report.

Acknowledgments

We thank the families participating in the DIPP study for making this study possible. We also thank the expert staff of the DIPP study for their excellent work with the participating research families and sample collection.

Footnotes

  • ↵* Shared senior authorship.

References

  1. 1.↵
    Caio G, Volta U, Sapone A, et al. Celiac disease: a comprehensive current review. BMC Medicine. 2019;17(1):142.
    OpenUrl
  2. 2.↵
    Holmes GKT, Muirhead A. Epidemiology of coeliac disease in a single centre in Southern Derbyshire 1958–2014. BMJ Open Gastro. 2017;4(1):e000137.
    OpenUrl
  3. 3.↵
    Assa A, Frenkel-Nir Y, Leibovici-Weissman Ya, et al. Anthropometric measures and prevalence trends in adolescents with coeliac disease: a population based study. Arch Dis Child. 2017;102(2):139–144.
    OpenUrlAbstract/FREE Full Text
  4. 4.↵
    Ilonen J, Kiviniemi M, Lempainen J, et al. Genetic susceptibility to type 1 diabetes in childhood – estimation of HLA class II associated disease risk and class II effect in various phases of islet autoimmunity. 2016;17(S22):8–16.
    OpenUrl
  5. 5.↵
    Sen P, Carlsson C, Virtanen SM, et al. Persistent Alterations in Plasma Lipid Profiles Before Introduction of Gluten in the Diet Associated With Progression to Celiac Disease. Clin Transl Gastroenterol. 2019;10(5):1–10.
    OpenUrl
  6. 6.↵
    Rewers M, Liu E, Simmons J, Redondo MJ, Hoffenberg EJ. Celiac disease associated with type 1 diabetes mellitus. Endocrinol Metab Clin North Am. 2004;33(1):197–214, xi.
    OpenUrlCrossRefPubMedWeb of Science
  7. 7.↵
    Scaramuzza AE, Mantegazza C, Bosetti A, Zuccotti GV. Type 1 diabetes and celiac disease: The effects of gluten free diet on metabolic control. World J Diabetes. 2013;4(4):130–134.
    OpenUrl
  8. 8.↵
    Ludvigsson JF, Ludvigsson J, Ekbom A, Montgomery SM. Celiac disease and risk of subsequent type 1 diabetes: a general population cohort study of children and adolescents. Diabetes Care. 2006;29(11):2483–2488.
    OpenUrlAbstract/FREE Full Text
  9. 9.↵
    Kirchberg FF, Werkstetter KJ, Uhl O, et al. Investigating the early metabolic fingerprint of celiac disease – a prospective approach. Journal of Autoimmunity. 2016;72:95–101.
    OpenUrl
  10. 10.↵
    Lionetti E, Catassi C. The Role of Environmental Factors in the Development of Celiac Disease: What Is New? Diseases. 2015;3(4):282–293.
    OpenUrl
  11. 11.↵
    Whyte L, Kotecha S, Watkins W, Jenkins H. Coeliac disease is more common in children with high socio-economic status. Acta Paediatr. 2014;103(3):5.
    OpenUrl
  12. 12.↵
    Zingone F, West J, Crooks CJ, et al. Socioeconomic variation in the incidence of childhood coeliac disease in the UK. Arch Dis Child. 2015;100(5):466–473.
    OpenUrlAbstract/FREE Full Text
  13. 13.↵
    Olén O, Bihagen E, Rasmussen F, Ludvigsson JF. Socioeconomic position and education in patients with coeliac disease. Digestive and Liver Disease. 2012;44(6):471–476.
    OpenUrlCrossRefPubMed
  14. 14.↵
    Tyrrell J, Melzer D, Henley W, Galloway TS, Osborne NJ. Associations between socioeconomic status and environmental toxicant concentrations in adults in the USA: NHANES 2001–2010. Environment International. 2013;59:328–335.
    OpenUrlCrossRefPubMedWeb of Science
  15. 15.↵
    Buekers J, Colles A, Cornelis C, Morrens B, Govarts E, Schoeters G. Socio-Economic Status and Health: Evaluation of Human Biomonitored Chemical Exposure to Per- and Polyfluorinated Substances across Status. Int J Environ Res Public Health. 2018;15(12):2818.
    OpenUrl
  16. 16.↵
    Conway B, Innes KE, Long D. Perfluoroalkyl substances and beta cell deficient diabetes. Journal of Diabetes and its Complications. 2016;30(6):993–998.
    OpenUrl
  17. 17.
    Nisticò L, Iafusco D, Galderisi A, et al. Emerging Effects of Early Environmental Factors over Genetic Background for Type 1 Diabetes Susceptibility: Evidence from a Nationwide Italian Twin Study. The Journal of Clinical Endocrinology & Metabolism. 2012;97(8):E1483–E1491.
    OpenUrl
  18. 18.
    Predieri B, Iughetti L, Guerranti C, Bruzzi P, Perra G, Focardi SE. High Levels of Perfluorooctane Sulfonate in Children at the Onset of Diabetes. Int J Endocrinol. 2015;2015:234358.
    OpenUrl
  19. 19.↵
    Stene LC, Gale EAM. The prenatal environment and type 1 diabetes. Diabetologia. 2013;56(9):1888–1897.
    OpenUrlCrossRefPubMedWeb of Science
  20. 20.↵
    McGlinchey A, Sinioja T, Lamichhane S, et al. Prenatal Exposure To Environmental Chemicals Modulates Serum Phospholipids In Newborn Infants, Increasing Later Risk Of Type 1 Diabetes. bioRxiv. 2019(588350):588350.
  21. 21.↵
    Auricchio R, Galatola M, Cielo D, et al. A Phospholipid Profile at 4 Months Predicts the Onset of Celiac Disease in at-Risk Infants. Scientific Reports. 2019;9(1):14303.
    OpenUrl
  22. 22.↵
    Hoffmanová I, Sánchez D, Tučková L, Tlaskalová-Hogenová H. Celiac Disease and Liver Disorders: From Putative Pathogenesis to Clinical Implications. Nutrients. 2018;10(7):892.
    OpenUrl
  23. 23.↵
    Vuoristo M, Miettinen TA. The Role of Fat and Bile Acid Malabsorption in Diarrhoea of Coeliac Disease. Scandinavian Journal of Gastroenterology. 1987;22(3):289–294.
    OpenUrlPubMed
  24. 24.↵
    Ejderhamn J, Samuelson K, Strandvik B. Serum Primary Bile Acids in the Course of Celiac Disease in Children. J Pediatr Gastroenterol Nutr. 1992;14(4):443–449.
    OpenUrlPubMed
  25. 25.↵
    Chiang JY. Recent advances in understanding bile acid homeostasis. F1000Research. 2017;6:2029–2029.
    OpenUrl
  26. 26.↵
    Barrasa JI, Olmo N, Lizarbe MA, Turnay J. Bile acids in the colon, from healthy to cytotoxic molecules. Toxicology in Vitro. 2013;27(2):964–977.
    OpenUrlCrossRefPubMed
  27. 27.↵
    Ramírez-Pérez O, Cruz-Ramón V, Chinchilla-López P, Méndez-Sánchez N. The Role of the Gut Microbiota in Bile Acid Metabolism. Annals of Hepatology. 2017;16:S21–S26.
    OpenUrl
  28. 28.↵
    Zhao W, Zitzow JD, Ehresman DJ, et al. Na+/Taurocholate Cotransporting Polypeptide and Apical Sodium-Dependent Bile Acid Transporter Are Involved in the Disposition of Perfluoroalkyl Sulfonates in Humans and Rats. Toxicological sciences : an official journal of the Society of Toxicology. 2015;146(2):363–373.
    OpenUrl
  29. 29.↵
    Salihović S, Dickens AM, Schoultz I, et al. Simultaneous determination of perfluoroalkyl substances and bile acids in human serum using ultra-high-performance liquid chromatography–tandem mass spectrometry. Anal Bioanal Chem 2019.
  30. 30.↵
    Folch J, Lees M, Sloane Stanley GH. A simple method for the isolation and purification of total lipides from animal tissues. J Biol Chem. 1957;226(1):497–509.
    OpenUrlFREE Full Text
  31. 31.↵
    Pluskal T, Castillo S, Villar-Briones A, Oresic M. MZmine 2: modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. BMC Bioinformatics. 2010;11:395.
    OpenUrlCrossRefPubMed
  32. 32.↵
    R Development Core Team. R: A language and environment for statistical computing. In:Vienna: R Foundation for Statistical Computing; 2018: http://www.R-project.org.
  33. 33.↵
    Rohart F, Gautier B, Singh A, Lê Cao K-A. mixOmics: An R package for ‘omics feature selection and multiple data integration. PLoS computational biology. 2017;13(11):e1005752.
    OpenUrl
  34. 34.↵
    Le Cao KA, Boitard S, Besse P. Sparse PLS discriminant analysis: biologically relevant feature selection and graphical displays for multiclass problems. BMC Bioinformatics. 2011;12:253.
    OpenUrlCrossRefPubMed
  35. 35.↵
    Westerhuis JA, Hoefsloot HC, Smit S, et al. Assessment of PLSDA cross validation. Metabolomics. 2008;4(1):81–89.
    OpenUrlCrossRefWeb of Science
  36. 36.
    Chong J, Wishart DS, Xia J. Using MetaboAnalyst 4.0 for Comprehensive and Integrative Metabolomics Data Analysis. Nucleic Acids Research. 2019;68(1):e86.
    OpenUrl
  37. 37.↵
    Jebbink J, Veenboer G, Boussata S, et al. Total bile acids in the maternal and fetal compartment in relation to placental ABCG2 expression in preeclamptic pregnancies complicated by HELLP syndrome. Biochimica et biophysica acta. 2015;1852(1):131–136.
    OpenUrl
  38. 38.↵
    Seki Y, Matsushita M, Kimura A, et al. Maternal and fetal circulation of unusual bile acids: A pilot study. Pediatr Int. 2011;53(6):1028–1033.
    OpenUrlPubMed
  39. 39.↵
    Macias RIR, Marin JJG, Serrano MA. Excretion of biliary compounds during intrauterine life. World journal of gastroenterology. 2009;15(7):817–828.
    OpenUrl
  40. 40.↵
    Beggs KM, McGreal SR, McCarthy A, et al. The role of hepatocyte nuclear factor 4- alpha in perfluorooctanoic acid- and perfluorooctanesulfonic acid-induced hepatocellular dysfunction. Toxicol Appl Pharmacol. 2016;304:18–29.
    OpenUrl
  41. 41.↵
    Pandak WM, Kakiyama G. The acidic pathway of bile acid synthesis: Not just an alternative pathway. Liver Research. 2019;3(2):88–98.
    OpenUrl
  42. 42.↵
    Buhrke T, Kruger E, Pevny S, Rossler M, Bitter K, Lampen A. Perfluorooctanoic acid (PFOA) affects distinct molecular signalling pathways in human primary hepatocytes. Toxicology. 2015;333:53–62.
    OpenUrl
  43. 43.↵
    Behr A-C, Plinsch C, Braeuning A, Buhrke T. Activation of human nuclear receptors by perfluoroalkylated substances (PFAS). Toxicology in Vitro. 2020;62:104700.
    OpenUrlCrossRef
  44. 44.
    Bjork JA, Butenhoff JL, Wallace KB. Multiplicity of nuclear receptor activation by PFOA and PFOS in primary human and rodent hepatocytes. Toxicology. 2011;288(1):8–17.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    Zhang L, Nichols RG, Correll J, et al. Persistent Organic Pollutants Modify Gut Microbiota-Host Metabolic Homeostasis in Mice Through Aryl Hydrocarbon Receptor Activation. Environmental health perspectives. 2015;123(7):679–688.
    OpenUrlCrossRefPubMed
  46. 46.↵
    Dean JM, Lodhi IJ. Structural and functional roles of ether lipids. Protein & Cell. 2018;9(2):196–206.
    OpenUrlCrossRefPubMed
  47. 47.↵
    Facciotti F, Ramanjaneyulu GS, Lepore M, et al. Peroxisome-derived lipids are self antigens that stimulate invariant natural killer T cells in the thymus. Nature Immunology. 2012;13(5):474–480.
    OpenUrlCrossRefPubMed
  48. 48.↵
    Van Kaer L, Wu L. Therapeutic Potential of Invariant Natural Killer T Cells in Autoimmunity. Front Immunol. 2018;9:519–519.
    OpenUrl
  49. 49.↵
    Oresic M, Gopalacharyulu P, Mykkanen J, et al. Cord serum lipidome in prediction of islet autoimmunity and type 1 diabetes. Diabetes. 2013;62(9):3268–3274.
    OpenUrlAbstract/FREE Full Text
Back to top
Posted April 15, 2020.
Download PDF
Data/Code
Email
Early-life exposure to perfluorinated alkyl substances modulates lipid metabolism in progression to celiac disease
Lisanna Sinisalu, Partho Sen, Samira Salihović, Suvi M. Virtanen, Heikki Hyöty, Jorma Ilonen, Jorma Toppari, Riitta Veijola, Matej Orešič, Mikael Knip, Tuulia Hyötyläinen
medRxiv 2020.04.02.20051359; doi: https://doi.org/10.1101/2020.04.02.20051359
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area