Skip to main content

Molecular remodeling of adipose tissue is associated with metabolic recovery after weight loss surgery

Abstract

Background

Bariatric surgery is an effective therapy for individuals with severe obesity to achieve sustainable weight loss and to reduce comorbidities. Examining the molecular signature of subcutaneous adipose tissue (SAT) following different types of bariatric surgery may help in gaining further insight into their distinct metabolic impact.

Results

Subjects undergoing biliopancreatic diversion with duodenal switch (BPD-DS) showed a significantly higher percentage of total weight loss than those undergoing gastric bypass or sleeve gastrectomy (RYGB + SG) (41.7 ± 4.6 vs 28.2 ± 6.8%; p = 0.00005). Individuals losing more weight were also significantly more prone to achieve both type 2 diabetes and dyslipidemia remission (OR = 0.75; 95%CI = 0.51–0.91; p = 0.03). Whole transcriptome and methylome profiling showed that bariatric surgery induced a profound molecular remodeling of SAT at 12 months postoperative, mainly through gene down-regulation and hypermethylation. The extent of changes observed was greater following BPD-DS, with 61.1% and 49.8% of up- and down-regulated genes, as well as 85.7% and 70.4% of hyper- and hypomethylated genes being exclusive to this procedure, and mostly associated with a marked decrease of immune and inflammatory responses. Weight loss was strongly associated with genes being simultaneously differentially expressed and methylated in BPD-DS, with the strongest association being observed for GPD1L (r2 = 0.83; p = 1.4 × 10–6).

Conclusions

Present findings point to the greater SAT molecular remodeling following BPD-DS as potentially linked with higher metabolic remission rates. These results will contribute to a better understanding of the metabolic pathways involved in the response to bariatric surgery and will eventually lead to the development of gene targets for the treatment of obesity.

Trial registration ClinicalTrials.gov NCT02390973.

Introduction

According to the World Health Organization, obesity prevalence has almost tripled since 1975, with over 650 million adults worldwide suffering from this condition in 2016 [1]. Obesity is related to many comorbid conditions including cardiovascular disease, type 2 diabetes, non-alcoholic fatty liver disease, obstructive sleep apnea, reproductive dysfunction, musculoskeletal disorders, certain types of cancer, as well as psychosocial consequences [2, 3]. Unfortunately, non-surgical treatment of severe obesity using modalities such as diet, exercise, or even pharmacological treatment have low to moderate success rates, especially when considering medium- to long-term weight maintenance [4]. Bariatric surgery, also known as metabolic surgery, emerged as an effective therapy for individuals with severe obesity to achieve significant sustainable weight loss, as well as to reduce the associated comorbidities [3, 5, 6]. The term metabolic surgery was proposed to acknowledge the physiological changes induced by these approaches, which contribute to a more favorable metabolic profile following the surgery, beyond the traditional view that these surgeries induce beneficial effects through weight loss alone [6,7,8]. Yet, these physiological changes are still not fully elucidated.

Among the different types of bariatric surgery, sleeve gastrectomy (SG) is one the most common and simple from a surgical standpoint, consisting of a surgical removal of about 80 percent of the stomach along the greater curvature, physically restricting food intake. In addition to reducing gastric volume, Roux-en-Y gastric bypass (RYGB) also decreases the efficiency of nutrient absorption in the small intestine by creating a little gastric pouch directly connected to the jejunum, bypassing the duodenum. Biliopancreatic diversion with duodenal switch (BPD-DS) is the most complex of these bariatric procedures, combining gastric restriction induced by SG and a more marked malabsorption than that observed after RYGB. Briefly, a pylorus-preserving SG procedure is combined with a transection of the duodenum near the pylorus creating an alimentary limb that connects the biliary limb near the ileo-cecal valve to create a short common channel where nutrients are absorbed [9].

Epigenetics may help in our understanding of how an individual will respond to bariatric surgery as the latter may be viewed as an environmental factor modifying the epigenome, although certain epigenetic marks may be inheritable [10, 11]. DNA methylation is the most widely investigated epigenetic mechanism, and some studies have predicted weight loss or weight regain after bariatric surgery according to baseline gene methylation [12, 13]. Bariatric surgery also promotes modifications in methylation profiles of individuals with obesity, which are more akin to those who are lean [14, 15]. Some authors have also observed lower overall methylation levels after RYGB in subcutaneous adipose tissue (SAT) [16, 17]. Conversely, more recent studies have observed higher methylation levels at cytosine-phosphate-guanine dinucleotides, or CpG sites, after RYGB and SG procedures in peripheral blood [15], as well as higher global methylation levels in skeletal muscle after RYGB [18]. These discrepancies may partly be explained by tissue-specific DNA methylation [19]. As compared to RYGB and SG, BPD-DS is a surgery that creates a greater nutrient malabsorption, due to a reduced length of the intestinal segment allowed for absorption, especially impacting dietary lipids [20, 21], and it has been proven to be particularly effective among individuals with severe obesity [22, 23]. By contrast, the impact of BPD-DS on the epigenetic profile of SAT is almost completely unknown.

The SAT is much more than a site of storage for excess energy intake, and it is recognized to play a crucial role in energy homeostasis control and inflammation [24]. Although significant modifications occur in the SAT after surgically induced weight loss [25, 26], only few studies have investigated gene expression changes in this tissue [27,28,29,30,31]. Genes previously identified as differentially expressed in SAT following bariatric surgery are involved in lipid and energy metabolism, inflammatory and immunological responses, insulin signaling, cell differentiation, oxidative stress regulation and gene transcription [27, 30]. A recent study [32] has also observed long-term effects of RYGB on gene expression in abdominal SAT, with enriched pathways related to lipid metabolism, fat cell differentiation and immune response. Again, most of the studies examining gene expression changes in SAT have been conducted among individuals with obesity who had undergone RYGB or SG, but none had investigated the impact of BPD-DS on the transcriptomic SAT profile [27,28,29,30,31].

Examining molecular changes taking place following a weight-loss procedure may then help in gaining further insight into its distinct metabolic impact and may eventually aid in targeting the most beneficial intervention. Thus, our objective was to compare methylation and gene expression changes in SAT following three different types of bariatric surgeries, namely BPD-DS, RYGB and SG. We hypothesized that the more metabolically effective BPD-DS leads to greater modifications at the methylation and gene expression level than the extensively studied RYGB and SG bariatric surgeries.

Results

BPD-DS induced a more pronounced weight loss than RYGB + SG

Of the 21 participants, 7 underwent BPD-DS, 5 RYGB and 9 SG (Fig. 1). Altogether, participants from the entire cohort showed a preoperative mean body mass index (BMI) of 44.4 kg/m2 ± 6.1, with a %TWL of 32.7% ± 8.9 at 12 months following the surgery. Due to the similar mean percent total weight loss (%TWL) shown by both RYGB (31.9% ± 6.4) and SG (26.2% ± 6.5) participants, these two groups were combined into the RYGB + SG group and compared to participants undergoing BPD-DS. Table 1 shows detailed information about the anthropometric measurements of BPD-DS and RYGB + SG surgery groups. Characteristics of subjects from the three different surgery procedures, BPD-DS, RYGB and SG, are detailed in Additional file 1: Table S1. When comparing BPD-DS and RYGB + SG surgery groups, we found that BPD-DS participants had significantly higher mean body weight, BMI, waist circumference and fat mass than RYGB + SG participants at baseline (Table 1). The proportion of men and women was not significantly different between surgery groups (p = 0.4) Following the surgery, at 12 months postoperative, BPD-DS participants had significantly higher mean delta BMI and %TWL than RYGB + SG subjects. No significant differences were found according to mean body weight, BMI, percent excess weight loss (%EWL), waist circumference or body composition (fat and lean mass). BPD-DS participants showed significantly lower neck circumference at 12 months postoperative, as compared to RYGB + SG participants (Table 1). No significant differences were found for adipocyte size 12 months following the surgery between BPD-DS and RYGB + SG groups, but participants undergoing SG showed significantly larger adipocytes than those who had BPD-DS and RYGB separately (Additional file 1: Table S1; Figure S1).

Fig. 1
figure 1

Flow diagram of study participants

Table 1 Characteristics of participants

Total remission was significantly higher in BPD-DS

All participants included in the study had type 2 diabetes before the surgery, with a mean duration of 8.4 ± 8.3 years since diagnostic, based on the Canadian Diabetes Association guidelines [33] (fasting plasma glucose ≥ 7.0 mmol/L, or glycated hemoglobin (HbA1c) ≥ 6.5%, or 2 h glycemia in a 75 g oral glucose tolerance test (OGTT) ≥ 11.1 mmol/L, or random glycemia ≥ 11.1 mmol/L). Type 2 diabetes was treated with either oral hypoglycemic drugs (n = 16, 76%), insulin (n = 4, 19%), or diet alone (n = 1, 7%). Most of participants also presented other comorbidities such as hypertension (n = 18, 86%), dyslipidemia (n = 18, 86%), or cardiovascular disease (n = 5, 24%). At 12 months postoperative, 13 participants (62%) were in complete type 2 diabetes remission, 4 were (19%) in partial remission and 4 (19%) had improvements in their glycemic control. For dyslipidemia, one participant from BPD-DS and two from RYGB + SG had no preoperative data, leaving 12 (67%) participants with complete remission, and 6 (33%) with improved, unchanged, or abnormal values. Following BPD-DS, 100% of participants achieved complete remission of type 2 diabetes (7 out of 7) and dyslipidemia (6 out of 6), while this percentage was 50% following RYGB + SG for type 2 diabetes (7 out of 14) (Fisher’s p = 0.05) and dyslipidemia remission (6 out of 12) (Fisher’s p = 0.05). Consequently, total remission rate was significantly different between the two surgery groups, with 100% of participants showing total remission following BPD-DS (6 out of 6), as compared to only 25% (3 out of 12) following RYGB + SG (Fisher's p = 0.009).

Extensive transcriptomic remodeling occurred in the SAT of BPD-DS participants

To examine the impact of bariatric surgery on SAT gene expression, we first compared SAT transcriptomic profiles between baseline and 12 months postoperative for each surgery group. First, we found that most of differentially expressed genes (DEGs) were exclusive to the BPD-DS group. Concretely, 713 up- and 943 down-regulated DEGs, representing 61.1% and 49.8% from the total DEGs, respectively (Fig. 2A). Conversely, only a few were exclusive to the RYGB + SG surgery group, with 170 (14.6%) up- and 169 (8.9%) down-regulated DEGs (Fig. 2B). Thus, there were approximately four times more up-regulated and six times more down-regulated DEGs exclusive to BPD-DS than to RYGB + SG. Regarding DEGs common to both surgery groups, there were three times more down-regulated than up-regulated genes, that is 782 (41.3%) down- versus 283 (24.3%) up-regulated DEGs (Fig. 2A–B). Interestingly, among common DEGs, the mean fold change of down-regulated DEGs was more than 50% greater in the BPD-DS group (1.89 log2 FC) than in the RYGB + SG group (1.26 log2 FC) (Fig. 2C), and approximately 35% greater for up-regulated DEGs (1.14 versus 0.85 log2 FC) (Fig. 2D).

Fig. 2
figure 2

Gene expression changes in subcutaneous adipose tissue were more pronounced following BPD-DS. Panel A shows differentially expressed genes (DEGs) exclusive for BPD-DS (green dots) and common to both surgery groups (red dots). Panel B shows DEGs exclusive for RYGB-DS (blue dots) and common to both surgery groups (purple dots). DEGs were considered significant when false discovery rate (FDR)-corrected p-value < 0.05 and fold change (FC) > 1.5. Panels C and D show density plots of FC distribution among down-regulated and up-regulated DEGs, respectively. Green and blue colors stand for DEGs exclusive for BPD-DS and RYGB + SG, respectively. Red and purple colors stand for DEGs common to both surgery groups but showing the specific FC distribution for each BPD-DS and RYGB + SG surgery, respectively. Dotted lines stand for mean FC values for each surgery group

Most of the differentially methylated genes were exclusive to BPD-DS

Differentially methylated loci containing at least one differentially methylated CpG site were examined in SAT at baseline and 12 months following the surgery for each group. As for the DEGs, most of the differentially methylated genes (DMGs) were exclusive to BPD-DS, with 8 094 DMGs being significantly hypermethylated and 6 369 hypomethylated, representing 85.7% and 70.4% of the total DMGs (Fig. 3A). Only a few DMGs were exclusive to the RYGB + SG group, with 114 (1.2%) hypermethylated and 769 (8.5%) hypomethylated DMGs (Fig. 3B). Furthermore, the mean fold change of hypermethylated DMGs common to both surgery groups was 50% greater following BPD-DS (1.26 log2 FC), as compared to RYGB + SG (0.68 log2 FC) (Fig. 3C), as well as for hypomethylated DMGs, although to a lesser extent (Fig. 3D).

Fig. 3
figure 3

Most of gene methylation changes in subcutaneous adipose tissue occur following BPD-DS. Panel A shows differentially methylated genes (DMGs) exclusive for BPD-DS (green dots) and common to both surgery groups (red dots). Panel B shows DMGs exclusive for RYGB-DS (blue dots) and common to both surgery groups (purple dots). DMGs were defined as loci with at least one differentially methylated CpG site (false discovery rate (FDR)-corrected p-value < 0.05 and fold change (FC) > 1.5. Panels C and D show density plots of FC distribution among hypermethylated and hypomethylated DMGs, respectively. Green and blue colors stand for DEGs exclusive for BPD-DS and RYGB + SG, respectively. Red and purple colors stand for DMGs common to both surgery groups but showing the specific FC distribution for each BPD-DS and RYGB + SG surgery, respectively. Dotted lines stand for mean FC values for each surgery

Immune-related pathways were differentially enriched and down-regulated in the SAT of BPD-DS participants

Pathway enrichment analysis was used to gain further insights into the biological processes significantly enriched with both DEGs and DMGs. A total of 15 and 225 pathways were found to be significantly enriched with up- and down-regulated DEGs common to both surgeries, whereas 27 and 122 were significantly enriched with up- and down-regulated DEGs exclusive to the BPD-DS group. In contrast, we did not find any pathway significantly enriched with up- or down-regulated DEGs exclusive to RYGB + SG. For illustrative purposes, most significantly enriched pathways exclusive to BPD-DS are shown in Fig. 4A and compared to pathways enriched with common and RYGB + SG exclusive DEGs. First, it is worth noting that most of up-regulated pathways exclusive to BPD-DS were also enriched with DEGs common to both surgery groups and were linked to the establishment of protein localization, RNA catabolic processes and protein translation, among others. On the other hand, down-regulated pathways were much more numerous with both BPD-DS and common DEGs, as already mentioned, and almost all the top enriched pathways were shared between these two groups (Fig. 4A). Down-regulated pathways were mostly related to immune and inflammatory biological processes, such as neutrophil-mediated immunity, response to bacterial molecules, leukocyte chemotaxis or cytokine production.

Fig. 4
figure 4

Immune-related pathways were markedly down-regulated following bariatric surgery. Left panel shows top Gene Ontology-Biological Process (GO-BP) terms significantly enriched with up-regulated (red blocks, up) and down-regulated (blue blocks, down) differentially expressed genes (DEGs). Right panel shows top GO-BP terms significantly enriched with hypermethylated (red blocks, up) and hypomethylated (blue blocks, down) differentially methylated genes (DMGs). Each column represents pathways enriched with DEGs specific to BPD-DS, RYGB + SG or common to both surgery groups. Pathways were considered significantly enriched when composed with at least 20 DEGs or DMGs and with FDR-adjusted p-value < 0.05

It is worth highlighting that, among the top metabolic pathways significantly enriched with hypermethylated DMGs exclusive to BPD-DS, most of them were related to immunity and inflammation, as previously shown for down-regulated DEGs (Fig. 4B). Similarly, most of these pathways were also found to be significantly enriched with hypermethylated DMGs common to both surgery groups. Again, no pathways were significantly enriched with hypermethylated DMGs exclusive to the RYGB + SG group. On the other hand, extracellular matrix and structure organization, actin filament organization or cell-substrate adhesion were among the top biological processes significantly enriched with hypomethylated DMGs exclusive to BPD-DS, with only half of them being also enriched with DMGs common to both surgeries (Fig. 4B). Only the axonogenesis pathway was significantly enriched in both BPD-DS and RYGB + SG surgery groups.

Most of genes that were both differentially expressed and methylated were found in the SAT of BPD-DS participants

We assessed whether gene expression changes were taking place in the same genes in which epigenetic modifications were also observed. Globally, almost 70% of the DEGs exclusive to BPD-DS surgery group were also differentially methylated. Concretely, most of down-regulated (796, 84.4%; Fig. 5A) and almost half of up-regulated DEGs (347, 48.6%; Fig. 5B) exclusive to the BPD-DS surgery group were also identified as DMGs. On opposite, only a few down-regulated (4, 2.4%; Fig. 5A) and up-regulated DEGs (7, 4.1%; Fig. 5B) exclusive for RYGB + SG surgery group were also identified as DMGs. Finally, most of the genes simultaneously identified as DEGs and DMGs common to both surgery groups (126, 16.1%) were down-regulated DEGs, as compared to only 20 (7.1%) up-regulated DEGs.

Fig. 5
figure 5

Genes being simultaneously differentially expressed and methylated largely belonged to the BPD-DS surgery group. A and B show respectively the proportion of differentially expressed genes (DEGs) down- and up-regulated that are simultaneously identified as differentially methylated genes (DMGs). The proportion of DEGs common to both surgery groups (COMMON), as well as exclusive to BPD-DS (BPD) and RYGB + SG (GAS) is shown in the inner ring. The proportion of hypermethylated and hypomethylated DMGs is shown in the outer ring. C and D show respectively the proportion of hypermethylated and hypomethylated CpG sites located within body or promoter regions of genes being simultaneously DEGs and DMGs. The proportion of genes being simultaneously DEGs and DMGs is shown in the inner ring, and the proportion of CpG sites for each gene location and surgery group is shown in the outer ring

We further investigated whether differentially methylated CpG sites associated to genes being both DEGs and DMGs were located within critical regions for gene transcription or within intergenic regions. As illustrated in Fig. 5C, an important proportion of differentially hypermethylated CpG sites (37.3%) were located within promoter regions, and an even greater proportion within gene bodies (62.7%). Approximately the same proportions were observed for hypomethylated CpG sites (Fig. 5D). Interestingly, none of the CpG sites was located within intergenic regions, suggesting an actual relationship between differential methylation and gene expression changes in genes identified as both DEGs and DMGs.

Total weight loss was strongly associated with genes being simultaneously differentially expressed and methylated in BPD-DS

According to the results showing higher remission rates following BPD-DS, we also found that individuals losing more weight were also those more prone to achieve total remission (OR = 0.75; 95%CI = 0.51–0.91; p = 0.03) (Fig. 6A). Similarly, %TWL was also significantly associated with type 2 diabetes remission (OR = 0.78; 95%CI = 0.62–0.99; p = 0.04) and showed a trend for association with dyslipidemia remission (OR = 0.69; 95%CI = 0.46–1.03; p = 0.07). We then tested whether %TWL was associated with gene expression changes in SAT from baseline to 12 months postoperative. Within DEGs common to both surgeries, a total of 24 up-regulated (8%) and 38 down-regulated DEGs (5%) were significantly associated with %TWL. Among common DEGs, the strongest positive associations were observed for the GJC3 (r2 = 0.84) and APOE (r2 = 0.82) genes, and a strong inverse correlation was found for CD248 (r2 = 0.83) (Fig. 6D). %TWL was also significantly associated with 75 up-regulated (11%) and 49 down-regulated DEGs (5%) exclusive to BPD-DS, with GPD1L showing the strongest positive association (r2 = 0.83) (Fig. 6D). None of DEGs exclusive to RYGB + SG was significantly associated with %TWL. Most of down-regulated (33, 67%) and almost half of up-regulated DEGs (36, 48%) exclusive to BPD-DS and showing a significant association with %TWL were also identified as DMGs.

Fig. 6
figure 6

Differentially expressed and methylated genes were associated with weight loss, adipocyte size and neck circumference. A to C show the predicted probability (red dots from 0 to 1), obtained by binomial logistic regression, of each participant to have a complete (0) or a partial remission (1), based on %TWL, %adipocyte size and %neck circumference. OR is the odds ratio with 95% confidence intervals (CI) and P is the p value for the linear trend of association. Gray and blue dots refer BPD-DS and RYGB + SG, respectively. D to F respectively show associations between differentially expressed genes (DEGs) in each surgery group with the percentage of total weight loss (%TWL), adipocyte size change (%Adipocyte) and neck circumference change (%Neck). Green, blue and red dots respectively stand for associations at non-adjusted p < 0.05 with DEGs exclusive to BPD-DS, RYGB + SG or common to both surgery groups. Grey dots represent not significant associations. Dot size is proportional to the magnitude (r2) of the association. Results are from multivariate linear regression models adjusted for sex, age and BMI

On the other hand, adipocyte size change did not show a significant impact on total remission rates (OR = 0.95; 95%CI = 0.86–1.02; p = 0.19) (Fig. 6B), nor on type 2 diabetes or dyslipidemia remission, and was only significantly associated with four DEGs, all of them being down-regulated DEGs exclusive to BPD-DS: INSYN1 (r2 = 0.68), SRXN1 (r2 = 0.65), CORO1C (r2 = 0.65) and SNPH (r2 = 0.64) (Fig. 6B). Interestingly, all of them were also significantly associated with %TWL, and three were also identified as DMGs (INSYN1, CORO1C and SNPH).

Finally, neck circumference change was also found to be a significant predictor of total remission (OR = 0.71; 95%CI = 0.47–0.92; p = 0.04) (Fig. 6C) and dyslipidemia remission (OR = 0.74; 95%CI = 0.55–0.99; p = 0.04), and showed a trend for association with type 2 diabetes remission (OR = 0.80; 95%CI = 0.63–1.01; p = 0.06). Accordingly, 110 up-regulated and 6 down-regulated DEGs exclusive to BPD-DS were significantly associated to this variable, as well as 23 up-regulated and 9 down-regulated DEGs common to both surgery groups, but none were exclusive to RYGB + SG. A total of 47 DEGs significantly associated with neck circumference change were also associated with %TWL.

Discussion

Following bariatric surgery in any of the procedures analyzed, most of the participants lost 20% or more of their initial weight, which has been suggested as a threshold for intervention success [34, 35]. Also, global transcriptomic and methylomic findings highlighted that, independently of the type of procedure, bariatric surgery induced a profound molecular remodeling in the SAT of patients with severe obesity, mainly through gene down-regulation and hypermethylation. However, weight loss at 12 months postoperative was far more important in participants who underwent BPD-DS, as compared to those undergoing RYGB + SG, which may partly explain the more extensive transcriptomic and methylomic modifications observed in the SAT of BPD-DS participants. It should be noted that less than 35% of all DEGs and only 17% of DMGs were common to both BPD-DS and RYGB + SG surgery groups. Moreover, the extent of gene down-regulation and hypermethylation among common DEGs and DMGs was 50% greater following BPD-DS than after RYGB + SG. Even more striking was that more than 50% of DEGs and almost 80% of DMGs were exclusive to BPD-DS, suggesting more extensive transcriptomic and methylomic modifications in SAT following this surgery, which also represents 70% of genes being simultaneously identified as DEGs and DMGs. Interestingly, an important proportion of differentially methylated CpG sites was located within promoter regions, through which gene expression can be regulated. Part of the methylation profile may be acquired during embryogenesis and is thought to remain stable across the lifespan [36]. However, some methylation marks may also be modified from the exposition to environmental factors such as diet, exercise, obesity, ageing or bariatric surgery [10, 11, 36].

Greater remission rates in type 2 diabetes following BPD-DS in comparison to RYGB and SG procedures have been previously reported [23, 37, 38]. These beneficial effects on type 2 diabetes are maintained over time with slightly more than 90% of the patients who discontinued diabetic therapy 10 years after undergoing BPD-DS surgery [21]. Besides the impressively favorable impact on plasma glucose and insulin levels, other beneficial metabolic shifts have been reported to be more persistent in long-term follow-ups with BPD-DS compared to RYGB, such as improved lipid profile and blood pressure lowering [21, 39, 40]. Whether these greater modifications are due to a more important and sustained weight loss or to more profound metabolic modifications after BPD-DS is still unknown. In this regard, it is worth highlighting that other factors than surgery type itself, such as sex and initial BMI, have been revealed to have a critical impact on weight loss and health outcomes following a bariatric surgery, with heavier men usually having the worst prognosis [41,42,43]. We also previously observed that men are overrepresented among subjects with a reduced weight loss response after surgery and, more importantly, that initial BMI is the best predictor of weight loss [44]. Concretely, we reported that the probability to show reduced weight loss following bariatric surgery significantly increases with initial BMI and mostly in men. In the present study, we tried to minimize a potential sex bias by keeping similar proportions of men and women among the different procedures, as well as by excluding for the analysis all the transcripts and CpG sites located on sexual chromosomes. Similarly, we used %TWL instead of %EWL as a measure of body weight loss, in an effort to reduce the impact of initial BMI, as recently reported [45]. Also, in order to take into account baseline differences between participants, differential gene expression analysis was performed by using a paired design, and linear as well as logistic models were adjusted for age, sex and initial BMI. Herein, having these considerations into account, we have established a potential link between a greater weight loss reduction with a more extensive transcriptomic and methylomic remodeling in SAT, which ultimately may contribute to these metabolic improvements. Similarly, gene expression changes in SAT were strongly associated with the reduction of neck circumference following bariatric surgery, which was previously suggested to be a reliable predictor of remission rates [46].

In animal models, metabolic changes seem to be mostly attributable to the malabsorptive effect of BPD-DS [47]. Vink et al. [48] have observed that a very-low caloric diet compared to a low caloric diet with similar weight loss induced greater gene expression modifications in pathways related to mitochondrial function, adipogenesis, immunity and inflammation. Thus, it is possible that the more pronounced effects ascertained in this study in BPD-DS compared to RYGB + SG were also partly attributable to a decrease in lipid absorption [49]. Moreover, the jejunum maybe crucial in regulating insulin sensitivity and is completely bypassed in BPD-DS [49]. Among mouse models, a high-fat compared to a low-fat isocaloric diet led to greater modifications in gene expression and methylation [50]. However, mice in the high-fat diet group were also significantly heavier following the diet than mice in the low-fat group [50]. It is not known whether diet or weight gain was responsible to induce these greater changes. In this regard, different bariatric procedures may also lead to distinct changes in diet, but actual data do not allow to test whether these changes are also impacting the present results.

In our study, many significantly enriched biologic processes related to protein translation and ribosomal activity were observed for up-regulated DEGs exclusive to BPD-DS and common to both surgery groups. Enrichments in similar pathways have been previously reported after RYGB [32] and diet-induced weight loss [48]. Genes involved in protein translation are also differently expressed among metabolically unhealthy individuals with obesity [51]. Protein translation may be regulated through an enhanced insulin sensitivity following bariatric surgery or dietary-induced weight loss [52]. For instance, insulin activates eukaryotic initiation and elongation factors, and increases the cellular content of ribosomal proteins [53]. Herein, many genes encoding ribosomal proteins (RPS and RPL genes) were up-regulated. In addition, genes encoding eukaryotic translation initiation factor 4E binding proteins (EIF4EBP1, EIF4EBP2 and EIF4EBP3 genes) and eukaryotic translation initiation factors (EIF genes) were found among the most up-regulated DEGs. Although these genes all encode proteins involved in translation, they may also be linked to adiposity, adipogenesis or glucose homeostasis. In fact, some EIF genes have been reported to correlate with genes encoding insulin receptor (INSR) and insulin receptor substrate-1 (IRS-1) [54], which were both also significantly up-regulated in the present study.

More than a simple reservoir of energy surplus, SAT has important endocrine and paracrine functions which regulate many biological processes [55]. Overweight and obesity may lead to SAT dysfunction including several perturbations, as an impaired expandability and adipocyte hypertrophy, altered innate and adaptative immune functions, changes in the secretion of pro- and anti-inflammatory peptides and eventually fibrosis [56, 57]. Transcriptomic modifications in pathways related to immunity and/or inflammation following weight loss, induced either by bariatric surgery or by diet, have been previously observed [27, 28]. Here, almost all the biological processes found to be significantly enriched with down-regulated DEGs were related to immune-related functions. Interestingly, most of such pathways were also significantly enriched with hypermethylated DMGs, pointing to a marked rearrangement of inflammation and immune molecular processes in the SAT of study participants. Although this effect was observed independently of the bariatric procedure analyzed, it was particularly pronounced following BPD-DS. By contrast, the fact that none of these pathways was significantly enriched with DEGs or DMGs exclusive to RYGB + SG emphasizes the procedure-specific nature of many of these molecular changes. Among significantly enriched immune-related pathways, biological processes such as T cell activation, leukocyte cell–cell adhesion, neutrophil activation and ERK1/2 cascades were significantly enriched. It has been observed that individuals with obesity have an increased quantity of T lymphocytes in their SAT [58]. Moreover, T regulatory lymphocytes have been shown to be reduced following bariatric surgery [59]. Neutrophils are an essential component of the innate immune response. Following RYGB, Poitou et al. [60] identified several DEGs which were related to neutrophil-mediated inflammation. In their study [60], DEGs related to neutrophils function or activity included S100A8, S100A9, S100A12, CD300E, VNN2, FPR2 and APOBEC3A, which were all also differentially expressed in both surgery groups of the present study, but around twice as much in BPD-DS than in RYGB + SG. Interestingly, Kerr et al. [32] have observed that there were continuously down-regulated genes 5 years following RYGB, despite significant weight regain occurring from 2 to 5 years postoperative. The authors observed that these continuously down-regulated genes were involved in biological processes such as cytokine production, cell chemotaxis and neutrophil activation [32], suggesting that these gene modifications might not be linked exclusively to weight variations. Moreover, these long-term changes in gene expression may be sustained through epigenetic mechanisms.

Two genes encoding for cytokines were among DEGs with the greatest change extent for both surgery groups, CSF3, which encodes for colony stimulating factor 3, a cytokine that has been reported to be elevated among individuals with obesity [61], and IL6, which is a well-known cytokine involved in inflammation. Following weight loss induced either by diet or bariatric surgery, a down-regulation in gene expression was previously noticed for genes such as CCL2 [28], NFKB1 [32], NLRP3 [32], HIF1A [27], CLEC7A [27] and IL4R [27]. These genes have also been found to be significantly down-regulated herein. Among them, HIF1A may indirectly activate NLRP3, which encodes for NLRP3 inflammasomes, contributing to the inflammatory responses via IL-1β activation, which is down-regulated to a greater extent in BPD-DS in this study, and could be of key importance in the development of type 2 diabetes [62].

Biological processes significantly enriched with hypomethylated DMGs were mostly related to extracellular structure and matrix organization, actin filament organization, cell-substrate and matrix adhesion, as well as cell-substrate junction and assembly, among others. These changes were again more pronounced following BPD-DS. Kelehmainen et al.[63] reported a down-regulation of DEGs involved in extracellular matrix following weight loss. The excessive accumulation of extracellular matrix components associated with obesity can lead to adipose tissue fibrosis which contributes to the dysfunction of adipocytes [64]. Moreover, higher SAT fibrosis may lessen the weight loss response following RYGB [64]. In the present study, these pathways were hypomethylated but not up-regulated. Thus, the functional impact of this hypomethylation remains unknown. It is possible that changes in gene expression were transient and no longer present at 12 months postoperative, since they potentially occurred earlier following the bariatric surgery, as previously shown in skeletal muscle [31].

Among significantly down-regulated DEGs common to both surgery groups, a strong inverse association with %TWL was observed for CD248, a gene which encodes for tumor endothelial marker 1/endosialin, a transmembrane glycoprotein known to be expressed in proliferating tissues, especially during embryogenesis, tumor growth and inflammatory lesions [65]. More recently, Petrus et al. [66] have demonstrated that CD248 is up-regulated in the SAT of individuals with obesity and insulin resistance and is potentially involved in the response to hypoxia. They reported that both CCL2 and IL6, respectively involved in extracellular matrix remodeling and inflammation, were correlated positively with CD248 gene expression [66]. Among up-regulated DEGs common to both surgery groups, GJC3 was associated with %TWL and it has been reported to be down-regulated in obesity [67]. On the other hand, GPD1L was the up-regulated DEG exclusive to BPD-DS most strongly associated with %TWL. In a long-term follow-up study, GPD1L was reported to be regulated by weight loss and regain after RYGB [32]. Furthermore, GPD1L was recently identified as potentially playing a causal role in obesity and insulin resistance [68]. During weight loss and weight maintenance induced by a low caloric diet, GPD1L was found to be up-regulated, while being down-regulated during weight gain induced by a high-fat diet [68]. It is worth noting that in the present study, most of DEGs exclusive to BPD-DS and associated with %TWL were also identified as DMGs, while none of the DEGs exclusive to RYGB + SG group were significantly associated to %TWL. Moreover, a total of four down-regulated DEGs exclusive to BPD-DS also showed an association with adipocyte size change. Among them, CORO1C, a gene recently identified to be up-regulated in the SAT of individuals with obesity [69], was found to be closely linked to %TWL and was also identified as hypomethylated. From a broader perspective, methylomic changes observed in this study were mostly exclusive to BPD-DS, which points to an epigenetic-mediated mechanism by which gene expression changes in SAT may occur in a greater extent in patients undergoing this type of surgery.

Present findings thus provide evidence that BPD-DS induce larger methylomic and transcriptomic modifications than RYGB + SG, which may be partly explained by greater weight loss and malabsorption created by this surgical approach. However, it is also possible that BPD-DS participants, who had higher BMI and waist circumference before surgery, started with a more deteriorated metabolic profile than participants who underwent RYGB or SG, which ultimately led to the more extensive transcriptomic and methylomic modifications observed. Results shown herein were obtained at 12 months postoperative, and it is possible that participants may still be losing weight or not being weight stable, which could affect transcriptome and methylome profiles. However, it has been reported that weight loss is as its nadir around 12 to 18 months following either RYGB [70] or BPD-DS [71].

Conclusions

To our knowledge, this is the first study examining the impact of bariatric surgery on SAT transcriptomic and methylomic profiles by using two high throughput technologies, RNA sequencing and genome-wide DNA methylation analysis. Our findings provide a novel overview of the transcriptomic and methylomic changes taking place 12 months following a bariatric surgery, concretely for BPD-DS compared to RYGB and SG. These results also confirm those obtained in previous transcriptomic studies following RYGB and SG. For instance, many of the enriched biological pathways found herein are shared with those previously found but have been observed to be of greater magnitude following BPD-DS. Globally, enriched biological processes in SAT following BPD-DS pointed to a strong decrease in immune and inflammatory responses and to an increase in protein translation, as well as to a shift towards modifications in other components of SAT, such as extracellular structure and actin filaments. These results will contribute to a better understanding of the metabolic pathways involved in the response to bariatric surgery and will eventually lead to the development of potential gene targets for the treatment of obesity and its related complications.

Methods

Study population

A total of 32 subjects with severe obesity, defined as BMI greater than or equal to 35 kg/m2, and aged between 18 and 60 years old, were recruited from the bariatric surgery clinic of the Institut universitaire de cardiologie et de pneumologie de Québec (IUCPQ). Recruitment occurred from September 2015 to November 2017. Exclusion criteria included pregnancy or desired pregnancy during the study; previous esophageal, digestive or bariatric surgery; abnormal bowel habits including irritable bowel syndrome, unexplained intermittent vomiting, severe abdominal pain, as well as chronic diarrhea or constipation in the last 60 days; history of gastric or duodenal ulcers; hypoalbuminemia (< 35 g/L); history of renal, hepatic, cardiac or pulmonary severe disease; evidence of psychiatric problems that may affect the capacity to understand the project and comply with the medical, surgical and/or behavioral recommendations; history of drug use or alcohol abuse in the last 12 months and during the study, as well as history of gastrointestinal inflammatory diseases. Of the initial 32 participants enrolled in the study, analyses were finally conducted on the 21 subjects for which SAT biopsies were successfully performed during bariatric surgery and 12 months later, including 10 men and 11 women. The clinical trial REMISSION is registered at Clinicaltrials.gov (NCT02390973).

Short-term prospective study protocol

Clinical exams including fasting biochemistry and anthropometric measurements were performed preoperatively and 12 months following the bariatric surgery. Participants underwent either BPD-DS, RYGB or SG according to National Institutes of Health (NIH) consensus for gastrointestinal surgery criteria [72] and based on the surgeon–patient’s choice at the bariatric surgery clinic of IUCPQ. After the surgery, participants followed a standardized postoperative protocol including feeding and a supplementation with vitamins and minerals. A detailed description of the BDP-DS procedure is given elsewhere [73].

Anthropometric measurements

Height and body weight were measured preoperatively and 12 months following bariatric surgery, and BMI was calculated as the weight in kilograms divided by the height in square meters. As recommended by the American Society for Metabolic and Bariatric Surgery (ASMBS) for reporting weight loss outcomes [8], we present the following information in the results section: mean initial BMI, change in BMI, %TWL and %EWL. %TWL was calculated as follows: [(Initial Weight)—(Postoperative Weight)] / [(Initial Weight)] * 100. %EWL was calculated as follows: [(Initial Weight)—(Postoperative Weight)] / [(Initial Weight)—(Ideal Weight)] * 100. Ideal weight is defined as the weight corresponding to a BMI of 23 kg/m2. A recent systematic review[75]investigated weight loss outcomes of RYGB and SG concluded that %TWL should be preferred over %EWL to minimize baseline BMI influence [45]. In this view, %TWL was used in the present study as the main weight loss outcome. Neck circumference, recently reported as a reliable predictor for the success of bariatric surgery [46], was also measured preoperatively and 12 months following the surgery.

Remission measurements

Overnight fasting blood samples were collected in the morning of each visit. Briefly, cholesterol and triglyceride levels were measured in plasma and lipoprotein fractions with a Technicon RA analyzer (Bayer, Etobicoke, ON, Canada) using enzymatic methods. Dyslipidemia remission was qualified according to plasma levels of low- (LDL) and high-density lipoproteins (HDL), total cholesterol and triglycerides [8]. Glucose was measured using the glucose oxidase method and insulin was quantified by radioimmunoassay (Linco Research, St. Charles, MO, US). The homeostasis model assessment of insulin resistance (HOMA-IR) index was calculated using the following formula: fasting insulin (μU/mL) * fasting glucose (mmol/L)/22.5. Diabetes remission was defined as suggested by the ASMBS [8] (HbA1c < 6.0% and fasting glycemia < 7.0 mmol/L in the absence of anti-diabetic pharmacological treatment), partial remission (HbA1c 6%-6.4% and fasting glycemia 5.6 mmol/L-6.9 mmol/L), improvement (statistical reduction in HbA1c and fasting glycemia not meeting criteria for remission or decrease in antidiabetic medications requirement). For comparative purposes, intermediate dyslipidemia and type 2 diabetes remission rates were grouped into two larger groups: partial and complete remission. A novel parameter called total remission was defined herein as the complete remission of both type 2 diabetes and dyslipidemia.

Adipose tissue sampling

Samples of SAT were collected at the site of the surgical incision in the lower abdomen. Immediately following surgical removal, fresh adipose tissue samples were carried to the laboratory where a portion of each sample was flash frozen in liquid nitrogen and stored at − 80 °C for further RNA and DNA extraction. Another portion of SAT sample was digested and used for adipocyte isolation and cell sizing. Briefly, tissue samples were digested with collagenase in Krebs–Ringer-Henseleit buffer for 45 min at 37 °C according to a modified version of the Rodbell method, as previously described [74]. Cell suspensions were filtered through nylon mesh and washed with Krebs–Ringer-Henseleit buffer. To determine adipocyte diameter, pictures of 250 cells were taken with a light microscope and analyzed with Scion Image software [74].

RNA sequencing

Total RNA from SAT biopsies obtained before and 12 months following the surgery was extracted using a RNeasy Lipid Tissue Mini Kit (Qiagen, Mississauga, ON, Canada) following the manufacturer’s instructions and treated with DNase (Qiagen) to avoid DNA contamination. RNA integrity was evaluated using the Agilent 2100 Bioanalyzer system (Agilent, Santa Clara, CA, US). RNA sequencing was performed at the McGill University and Génome Québec Innovation Centre (MUGQIC). Library preparation was carried out using the Illumina NEB stranded mRNA library preparation kit (Illumina, San Diego, CA, US) and sequencing was performed on the Illumina NovaSeq 6000 S4 platform (Illumina) using 100 bp paired-end reads. Raw reads were first trimmed at 50 bases and at a Phred quality score of 30 using Trim Galore! (v0.6.4) [75], a wrapper tool around Cutadapt (v3.2) [76] and FastQC (v0.11.9) [77]⁠. Read quantification was performed using kallisto (v0.46.2) [78] with 100 bootstraps. Reads were aligned to the GRCh38 human reference transcriptome and transcripts located on sexual chromosomes were excluded for further analyses. The obtained transcript counts were used to infer gene-level abundance estimates with the tximport (v1.20.0) R package [79]. Gene expression was then normalized, and lowly expressed genes were filtered out with the filterByExpr function in edgeR (v3.34.0) [80], leaving a total of 18 862 genes for further analyses. Differential gene expression analysis was performed between pre-surgical and 12-month follow-up levels in edgeR using a paired design, which can be viewed as a generalization of a paired t-test. DEGs were considered at false discovery rate (FDR)-corrected p-value < 0.05 and fold change (FC) > 1.5.

Genome-wide methylation analysis

Genomic DNA of the 21 study participants was extracted from 200 mg of SAT biopsy samples obtained before and 12 months following surgery using the DNeasy Blood & Tissue kit (Qiagen). Following quantification of DNA using both NanoDrop Spectrophotometer (Thermo Scientific, Wilmington, DE, US) and PicoGreen DNA methods, DNA (1 μg) was bisulfite converted. Quantitative genome-wide methylation analysis was conducted using the EPIC platform (Illumina), interrogating over 850 000 CpG sites at single-nucleotide resolution. Methylation arrays were processed at the MUGQIC according to manufacturer’s instructions. All probes with low detection p-values (< 0.05) were removed, as well as those located in sex chromosomes. Polymorphic and cross-reactive probes were also excluded, leaving a total of 774 177 probes for further differential methylation analyses. Methylation data was normalized with the quantile method, as previously described [81], using the minfi (v1.38.0) R package [82]⁠. Methylation levels (beta values, β) were estimated as the ratio of signal intensity of the methylated alleles to the sum of methylated and unmethylated intensity signals. The β values varied from 0 (no methylation) to 1 (100% methylation). Differentially methylated CpG sites were considered at FDR-corrected p-value < 0.05 and FC > 1.5. (DMGs were defined as loci with at least one differentially methylated mapped CpG site.

Pathway enrichment analysis

The functional significance of DEGs and DMGs was explored by pathway enrichment analysis using the clusterProfiler v3.16.0 R package [83]. The clusterProfiler package implements statistical methods to analyze functional profiles of genes and gene clusters and produces FDR-adjusted p-values to show significantly enriched pathways. The Gene Ontology Biological Processes (GO-BP) database was used for functional enrichment analysis. Pathway enrichment analysis was performed with DEGs and DMGs common to both surgery groups analyzed, as well as with those exclusive for BPD-DS and RYGB + SG. Pathways were considered significantly enriched at FDR-adjusted p-value < 0.05 and composed with at least 20 DEGs or DMGs.

Statistics

Clinical data were checked for normality with the Kolmogorov–Smirnov test and two-group comparisons were tested with two-tailed Student’s t test for independent samples. Fisher exact test was used to compare remission success rates between surgery groups, as well as the proportion of men and women. Multivariate linear models adjusted by sex, age and BMI were implemented to test for associations between gene expression and gene methylation levels with %TWL, adipocyte size change and neck circumference change. Linear associations were considered significant when FDR-corrected p-value < 0.05. Binomial logistic regression was used to predict the probability of total remission success set as a dichotomous variable. %TWL, adipocyte size change and neck circumference change were tested as independent predictors, with age, sex and BMI set as covariables. All the statistical analysis were implemented in R (v4.1.0) [84]⁠.

Availability of data and materials

The datasets generated during the current study are not publicly available due to privacy and confidentiality reasons but are available from the corresponding author on reasonable request.

References

  1. World Health Organization. Obesity and overweight. 2021. Acessed from: https://www.who.int/news-room/fact-sheets/detail/obesity-and-overweight

  2. Center for disease control and prevention. Adults obesity facts. 2021. Acessed from: https://www.cdc.gov/obesity/data/adult.html

  3. Wolfe BM, Kvach E, Eckel RH. Treatment of Obesity. Circ Res. 2016;11:1844.

    Article  CAS  Google Scholar 

  4. Wing RR, Phelan S. Long-term weight loss maintenance. Am J Clin Nutr. 2005;82:222s.

    CAS  PubMed  Article  Google Scholar 

  5. O’Brien PE, Hindle A, Brennan L, Skinner S, Burton P, Smith A, et al. Long-term outcomes after bariatric surgery: a systematic review and meta-analysis of weight loss at 10 or more years for all bariatric procedures and a single-centre review of 20-year outcomes after adjustable gastric banding. Obes Surg. 2019;29(1):3.

    PubMed  Article  Google Scholar 

  6. Vidal J, Corcelles R, Jiménez A, Flores L, Lacy AM. Metabolic and bariatric surgery for obesity. Gastroenterology. 2017;152:1780.

    PubMed  Article  Google Scholar 

  7. Buchwald H, Buchwald JN. Metabolic (bariatric and nonbariatric) surgery for type 2 diabetes: a personal perspective review. Diabetes Care. 2019;42:331.

    PubMed  Article  Google Scholar 

  8. Brethauer SA, Kim J, el Chaar M, Papasavas P, Eisenberg D, Rogers A, et al. Standardized outcomes reporting in metabolic and bariatric surgery. Surg Obes Relat Dis. 2015;11:489.

    PubMed  Article  Google Scholar 

  9. Biertho L, Lebel S, Marceau S, Hould F-SS, Julien F, Biron S. Biliopancreatic diversion with duodenal switch. Surg Clin North Am. 2016;96:815–26.

    PubMed  Article  Google Scholar 

  10. Metere A, Graves CE. Factors influencing epigenetic mechanisms: is there a role for bariatric surgery? High Throughput. 2020;9(1):6.

    PubMed Central  Article  Google Scholar 

  11. Samblas M, Milagro FI, Martínez A. DNA methylation markers in obesity, metabolic syndrome, and weight loss. Epigenetics. 2019;14:421.

    PubMed  PubMed Central  Article  Google Scholar 

  12. Nicoletti CF, Pinhel MS, Noronha NY, Jácome A, Crujeiras AB, Nonino CB. Association of MFSD3 promoter methylation level and weight regain after gastric bypass: assessment for 3 y after surgery. Nutrition. 2020;70:110499.

    CAS  PubMed  Article  Google Scholar 

  13. Assem S, Abdelbaki TN, Mohy-El Dine SH, Ketat AF, Abdelmonsif DA. SERPINE-1 gene methylation and protein as molecular predictors of laparoscopic sleeve gastrectomy outcome. Obes Surg. 2020;30(7):2620.

    PubMed  Article  Google Scholar 

  14. Barres R, Kirchner H, Rasmussen M, Yan J, Kantor FR, Krook A, et al. Weight loss after gastric bypass surgery in human obesity remodels promoter methylation. Cell Rep. 2013. https://doi.org/10.1016/j.celrep.2013.05.017.

    Article  PubMed  Google Scholar 

  15. Fraszczyk E, Luijten M, Spijkerman AMW, Snieder H, Wackers PFK, Bloks VW, et al. The effects of bariatric surgery on clinical profile, DNA methylation, and ageing in severely obese patients. Clin Epigenet. 2020. https://doi.org/10.1186/s13148-019-0790-2.

    Article  Google Scholar 

  16. Dahlman I, Sinha I, Gao H, Brodin D, Thorell A, Rydén M, et al. The fat cell epigenetic signature in post-obese women is characterized by global hypomethylation and differential DNA methylation of adipogenesis genes. Int J Obes. 2015. https://doi.org/10.1038/ijo.2015.31.

    Article  Google Scholar 

  17. Benton MC, Johnstone A, Eccles D, Harmon B, Hayes MT, Lea RA, et al. An analysis of DNA methylation in human adipose tissue reveals differential modification of obesity genes before and after gastric bypass and weight loss. Genome Biol. 2015. https://doi.org/10.1186/s13059-014-0569-x.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Garcia LA, Day SE, Coletta RL, Campos B, Benjamin TR, de Filippis E, et al. Weight loss after roux-en-y gastric bypass surgery reveals skeletal muscle DNA methylation changes. Clin Epigent. 2021. https://doi.org/10.1186/s13148-021-01086-6.

    Article  Google Scholar 

  19. Ghosh S, Yates AJ, Frühwald MC, Miecznikowski JC, Plass C, Smiraglia D. Tissue specific DNA methylation of CpG islands in normal human adult somatic tissues distinguishes neural from non-neural tissues. Epigenetics. 2010;5(6):527.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  20. Biertho L, Simon-Hould F, Marceau S, Lebel S, Lescelleur O, Biron S. Current Outcomes of Laparoscopic Duodenal Switch. Ann Surg Innov Res. 2016. https://doi.org/10.1186/s13022-016-0024-7.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Marceau P, Biron S, Marceau S, Hould F-S, Lebel S, Lescelleur O, et al. Long-term metabolic outcomes 5 to 20 years after biliopancreatic diversion. Obes Surg. 2015;25(9):1584.

    PubMed  Article  Google Scholar 

  22. Bolckmans R, Himpens J. Long-term (>10 Yrs) outcome of the laparoscopic biliopancreatic diversion with duodenal switch. Ann Surg. 2016;264(6):1029.

    PubMed  Article  Google Scholar 

  23. Skogar ML, Sundbom M. Duodenal switch is superior to gastric bypass in patients with super obesity when evaluated with the bariatric analysis and reporting outcome system (BAROS). Obes Surg. 2017;27(9):2308.

    PubMed  PubMed Central  Article  Google Scholar 

  24. Longo M, Zatterale F, Naderi J, Parrillo L, Formisano P, Raciti GA, et al. Adipose tissue dysfunction as determinant of obesity-associated metabolic complications. Int J Mol Sci. 2019;20(9):3258.

    Article  CAS  Google Scholar 

  25. Adami GF, Carbone F, Montecucco F, Camerini G, Cordera R. Adipose tissue composition in obesity and after bariatric surgery. Obes Surg. 2019;29(9):3030.

    PubMed  Article  Google Scholar 

  26. Grenier-Larouche T, Carreau AM, Geloën A, Frisch F, Biertho L, Marceau S, et al. Fatty acid metabolic remodeling during type 2 diabetes remission after bariatric surgery. Diabetes. 2017;66:2743–55.

    CAS  PubMed  Article  Google Scholar 

  27. de Pinhel MA, Noronha SNY, Nicoletti CF, de Oliveira BAP, Cortes-Oliveira C, Pinhanelli VC, et al. Changes in global transcriptional profiling of women following obesity surgery bypass. Obes Surg. 2018;28:176.

    PubMed  Article  Google Scholar 

  28. Beisani M, Pappa S, Moreno P, Martínez E, Tarascó J, Granada ML, et al. Laparoscopic sleeve gastrectomy induces molecular changes in peripheral white blood cells. Clin Nutr. 2020;39:592.

    PubMed  Article  Google Scholar 

  29. ElGendy K, Malcomson FC, Bradburn DM, Mathers JC. Effects of bariatric surgery on DNA methylation in adults: a systematic review and meta-analysis. Surg Obes Relat Dis. 2020;16:128.

    PubMed  Article  Google Scholar 

  30. Pinhel MAS, Noronha NY, Nicoletti CF, Pereira VA, de Oliveira BA, Cortes-Oliveira C, et al. Changes in DNA methylation and gene expression of insulin and obesity-related gene PIK3R1 after roux-en-Y gastric bypass. Int J Mol Sci. 2020;21:4476.

    CAS  PubMed Central  Article  Google Scholar 

  31. Gancheva S, Ouni M, Jelenik T, Koliaki C, Szendroedi J, Toledo FGS, et al. Dynamic changes of muscle insulin sensitivity after metabolic surgery. Nat Commun. 2019. https://doi.org/10.1038/s41467-019-12081-0.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Kerr AG, Andersson DP, Rydén M, Arner P, Dahlman I. Long-term changes in adipose tissue gene expression following bariatric surgery. J Int Med. 2020. https://doi.org/10.1111/joim.13066.

    Article  Google Scholar 

  33. Punthakee Z, Goldenberg R, Katz P. Definition, classification and diagnosis of diabetes, prediabetes and metabolic syndrome. Can J Diabetes. 2018;42:10–5.

    Article  Google Scholar 

  34. Grover BT, Morell MC, Kothari SN, Borgert AJ, Kallies KJ, Baker MT. Defining weight loss after bariatric surgery: a call for standardization. Obes Surg. 2019. https://doi.org/10.1007/s11695-019-04022-z.

    Article  PubMed  Google Scholar 

  35. Corcelles R, Boules M, Froylich D, Hag A, Daigle CR, Aminian A, et al. Total Weight Loss as the Outcome Measure of Choice After Roux-en-Y Gastric Bypass. Obes Surg. 2016. https://doi.org/10.1007/s11695-015-2022-y.

    Article  PubMed  Google Scholar 

  36. van Dijk SJ, Molloy PL, Varinli H, Morrison JL, Muhlhausler BS. Epigenetics and human obesity. Int J Obes. 2015;39:85.

    Article  CAS  Google Scholar 

  37. Mingrone G, Panunzi S, de Gaetano A, Guidone C, Iaconelli A, Nanni G, et al. Bariatric–metabolic surgery versus conventional medical treatment in obese patients with type 2 diabetes: 5 year follow-up of an open-label, single-centre, randomised controlled trial. Lancet. 2015;386:964.

    PubMed  Article  Google Scholar 

  38. Skroubis G, Kouri N, Mead N, Kalfarentzos F. Long-term results of a prospective comparison of roux-en-Y gastric bypass versus a variant of biliopancreatic diversion in a non-superobese population (BMI 35–50 kg/m2). Obes Surg. 2014;24:197.

    PubMed  Article  Google Scholar 

  39. Sjöström CD, Peltonen M, Sjöström L. Blood pressure and pulse pressure during long-term weight loss in the obese: the swedish obese subjects (SOS) intervention study. Obes Res. 2001;9:188.

    PubMed  Article  Google Scholar 

  40. Edholm D, Svensson F, Näslund I, Karlsson FA, Rask E, Sundbom M. Long-term results 11 years after primary gastric bypass in 384 patients. Surg Obes Relat Dis. 2013;9:708.

    PubMed  Article  Google Scholar 

  41. MacHado MB, Velasco IT, Scalabrini-Neto A. Gastric bypass and cardiac autonomic activity: influence of gender and age. Obes Surg. 2009;19:332–8.

    PubMed  Article  Google Scholar 

  42. Roberts K, Duffy A, Kaufman J, Burrell M, Dziura J, Bell R. Size matters: gastric pouch size correlates with weight loss after laparoscopic Roux-en-Y gastric bypass. Surg Endosc. 2007;21:1397–402.

    CAS  PubMed  Article  Google Scholar 

  43. Melton GB, Steele KE, Schweitzer MA, Lidor AO, Magnuson TH. Suboptimal weight loss after gastric bypass surgery: correlation of demographics, comorbidities, and insurance status with outcomes. J Gastrointest Surg. 2008;12:250–5.

    PubMed  Article  Google Scholar 

  44. de Toro-Martín J, Guénard F, Tchernof A, Pérusse L, Marceau S, Vohl M-C. Polygenic risk score for predicting weight loss after bariatric surgery. JCI Insight. 2018. https://doi.org/10.1172/jci.insight.122011.

    Article  PubMed  PubMed Central  Google Scholar 

  45. van Rijswijk A-S, van Olst N, Schats W, van der Peet DL, van de Laar AW. What is weight loss after bariatric surgery expressed in percentage total weight Loss (%TWL)? a systematic Review. Obes Surg. 2021;31:3833.

    PubMed  Article  Google Scholar 

  46. Bioletto F, Pellegrini M, D’Eusebio C, Boschetti S, Rahimi F, de Francesco A, et al. Development and validation of a scoring system for pre-surgical and early post-surgical prediction of bariatric surgery unsuccess at 2 years. Sci Rep. 2021;11:1–10.

    Article  Google Scholar 

  47. Baraboi E-D, Li W, Labbé SM, Roy M-C, Samson P, Hould F-S, et al. Metabolic changes induced by the biliopancreatic diversion in diet-induced obesity in male rats: the contributions of sleeve gastrectomy and duodenal switch. Endocrinology. 2015;156:13106.

    Article  CAS  Google Scholar 

  48. Vink RG, Roumans NJ, Fazelzadeh P, Tareen SHK, Boekschoten M, v, van Baak MA, et al. Adipose tissue gene expression is differentially regulated with different rates of weight loss in overweight and obese humans. Int J Obes. 2017;41:309.

    CAS  Article  Google Scholar 

  49. Castagneto-Gissey L, Angelini G, Casella-Mariolo JR, Marini P, Mingrone G, Casella G. The jejunum is the key factor in insulin resistance. Surg Obes Relat Dis. 2020;16:509.

    PubMed  Article  Google Scholar 

  50. Keleher MR, Zaidi R, Hicks L, Shah S, Xing X, Li D, et al. A high-fat diet alters genome-wide DNA methylation and gene expression in SM/J mice. BMC Geno. 2018. https://doi.org/10.1186/s12864-018-5327-0.

    Article  Google Scholar 

  51. Gaye A, Doumatey AP, Davis SK, Rotimi CN, Gibbons GH. Whole-genome transcriptomic insights into protective molecular mechanisms in metabolically healthy obese African Americans. Geno Med. 2018. https://doi.org/10.1038/s41525-018-0043-x.

    Article  Google Scholar 

  52. Yoshino M, Kayser BD, Yoshino J, Stein RI, Reeds D, Eagon JC, et al. Effects of diet versus gastric bypass on metabolic function in diabetes. New Engl J Med. 2020;383:721.

    CAS  PubMed  Article  Google Scholar 

  53. Proud CG. Regulation of protein synthesis by insulin. Biochem Soci Trans. 2006;34:213.

    CAS  Article  Google Scholar 

  54. MacLaren R, Cui W, Simard S, Cianflone K. Influence of obesity and insulin sensitivity on insulin signaling genes in human omental and subcutaneous adipose tissue. J Lipid Res. 2008;49:308.

    CAS  PubMed  Article  Google Scholar 

  55. Scheja L, Heeren J. The endocrine function of adipose tissues in health and cardiometabolic disease. Nat Rev Endocrinol. 2019;15:507.

    CAS  PubMed  Article  Google Scholar 

  56. Zatterale F, Longo M, Naderi J, Raciti GA, Desiderio A, Miele C, et al. Chronic adipose tissue inflammation linking obesity to insulin resistance and type 2 diabetes. Frontiers in physiology. 2020;1607.

  57. DeBari MK, Abbott RD. Adipose tissue fibrosis: mechanisms, models, and importance. Int J Mol Sci. 2020;21:6030.

    CAS  PubMed Central  Article  Google Scholar 

  58. Duffaut C, Zakaroff-Girard A, Bourlier V, Decaunes P, Maumus M, Chiotasso P, et al. Interplay between human adipocytes and T lymphocytes in obesity. Arteriosclerosis Thrombo Vasc Biol. 2009. https://doi.org/10.1161/ATVBAHA.109.192583.

    Article  Google Scholar 

  59. Villarreal-Calderón JR, Cuéllar RX, Ramos-González MR, Rubio-Infante N, Castillo EC, Elizondo-Montemayor L, et al. Interplay between the adaptive immune system and insulin resistance in weight loss induced by bariatric surgery. Oxid Med Cell Longev. 2019. https://doi.org/10.1155/2019/3940739.

    Article  PubMed  PubMed Central  Google Scholar 

  60. Poitou C, Perret C, Mathieu F, Truong V, Blum Y, Durand H, et al. Bariatric surgery induces disruption in inflammatory signaling pathways mediated by immune cells in adipose tissue: A RNA-Seq study. PLOS ONE. 2015;10:e0125718.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  61. Cancello R, Henegar C, Viguerie N, Taleb S, Poitou C, Rouault C, et al. Reduction of macrophage infiltration and chemoattractant gene expression changes in white adipose tissue of morbidly obese subjects after surgery-induced weight loss. Diabetes. 2005;54:2277.

    CAS  PubMed  Article  Google Scholar 

  62. Ding S, Xu S, Ma Y, Liu G, Jang H, Fang J. Modulatory mechanisms of the NLRP3 inflammasomes in diabetes. Biomolecules. 2019. https://doi.org/10.3390/biom9120850.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Kolehmainen M, Salopuro T, Schwab US, Kekäläinen J, Kallio P, Laaksonen DE, et al. Weight reduction modulates expression of genes involved in extracellular matrix and cell death: the GENOBIN study. Int J Obes. 2008;32(1):292–303.

    CAS  Article  Google Scholar 

  64. Abdennour M, Reggio S, le Naour G, Liu Y, Poitou C, Aron-Wisnewsky J, et al. Association of adipose tissue and liver fibrosis with tissue stiffness in morbid obesity: links with diabetes and bmi loss after gastric bypass. J Clin Endocrinol Metab. 2014;99(3):898–907.

    CAS  PubMed  Article  Google Scholar 

  65. Valdez Y, Maia M, M. Conway E. CD248: reviewing its role in health and disease. Current Drug Targets. 2012;13:432.

    CAS  PubMed  Article  Google Scholar 

  66. Petrus P, Fernandez TL, Kwon MM, Huang JL, Lei V, Safikhan NS, et al. Specific loss of adipocyte CD248 improves metabolic health via reduced white adipose tissue hypoxia, fibrosis and inflammation. EBioMedicine. 2019;44:489.

    PubMed  PubMed Central  Article  Google Scholar 

  67. Lu Z, Meng L, Sun Z, Shi X, Shao W, Zheng Y, et al. Differentially expressed genes and enriched signaling pathways in the adipose tissue of obese people. Front Genet. 2021. https://doi.org/10.3389/fgene.2021.620740.

    Article  PubMed  PubMed Central  Google Scholar 

  68. He H, Sun D, Zeng Y, Wang R, Zhu W, Cao S, et al. A Systems genetics approach identified GPD1L and its molecular mechanism for obesity in human adipose tissue. Sci Rep. 2017. https://doi.org/10.1038/s41598-017-01517-6.

    Article  PubMed  PubMed Central  Google Scholar 

  69. Joshi H, Vastrad B, Joshi N, Vastrad C, Tengli A, Kotturshetti I. Identification of key pathways and genes in obesity using bioinformatics analysis and molecular docking studies. Front Endocrinol. 2021. https://doi.org/10.3389/fendo.2021.628907.

    Article  Google Scholar 

  70. Courcoulas AP, King WC, Belle SH, Berk P, Flum DR, Garcia L, et al. 7-Year weight trajectories and health outcomes in the longitudinal assessment of bariatric surgery (LABS) study. JAMA Surg. 2018;153:427.

    PubMed  Article  Google Scholar 

  71. Topart P, Becouarn G, Delarue J. Weight loss and nutritional outcomes 10 years after biliopancreatic diversion with duodenal switch. Obes Surg. 2017;27(11):1645–50.

    PubMed  Article  Google Scholar 

  72. Gastrointestinal surgery for severe obesity: National institutes of health consensus development conference statement. Am J Clini Nutr. 1992;55.

  73. Marceau P, Biron S, Hould F-S, Lebel S, Marceau S, Lescelleur O, et al. Duodenal switch longterm results. Obes Surg. 2007;17(11):1421–30.

    PubMed  Article  Google Scholar 

  74. Grenier-Larouche T, Galinier A, Casteilla L, Carpentier AC, Tchernof A. Omental adipocyte hypertrophy relates to coenzyme Q10 redox state and lipid peroxidation in obese women. J Lipid Res. 2015;56:1985.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  75. Krueger F. Trim Galore. 2019. http://www.bioinformatics.babraham.ac.uk/projects/trim_galore/.

  76. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J EMBnet Stichting. 2011;17:10–2.

    Article  Google Scholar 

  77. Andrews S. FastQC: a quality control tool for high throughput sequence data. 2019. http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.

  78. Bray NL, Pimentel H, Melsted P, Pachter L. Near-optimal probabilistic RNA-seq quantification. Nat Biotechnol Nat Publ Gr. 2016;34:525–7.

    CAS  Article  Google Scholar 

  79. Soneson C, Love MI, Robinson MD. Differential analyses for RNA-seq: transcript-level estimates improve gene-level inferences. F1000Res. Faculty of 1000 Ltd; 2015;4.

  80. Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics Oxf Univ Press. 2009;26:139–40.

    Article  CAS  Google Scholar 

  81. Touleimat N, Tost J. Complete pipeline for Infinium® Human Methylation 450K BeadChip data processing using subset quantile normalization for accurate DNA methylation estimation. Epigenomics. 2012. https://doi.org/10.2217/epi1221.

    Article  PubMed  Google Scholar 

  82. Fortin J-P, Triche TJ, Hansen KD. Preprocessing, normalization and integration of the illumina humanmethylationEPIC array with minfi. Bioinformatics. 2016;33:btw691.

    Article  CAS  Google Scholar 

  83. Yu G, Wang LG, Han Y, ClusterProfiler He QY. An R package for comparing biological themes among gene clusters. OMICS J Integr Biol. 2012;16:284–7.

    CAS  Article  Google Scholar 

  84. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2021. https://www.R-project.org/.

Download references

Acknowledgements

We would like to acknowledge the contribution of surgeons, nurses, and the medical team of the bariatric surgery program at IUCPQ. The coinvestigators and collaborators of the REMISSION study are (alphabetical order): Bégin C, Biertho L, Bouvier M, Biron S, Cani P, Carpentier A, Dagher A, Dubé F, Fergusson A, Fulton S, Hould FS, Julien F, Kieffer T, Laferrère B, Lafortune A, Lebel S, Lescelleur O, Levy E, Marette A, Marceau S, Michaud A, Picard F, Poirier P, Richard D, Schertzer J, Tchernof A, Vohl MC.

Funding

This study is supported by a Team grant from the Canadian Institutes of Health Research (CIHR) on bariatric care (TB2-138776) and an Investigator-initiated study grant from Johnson & Johnson Medical Companies (Grant ETH-14-610). Funding sources for the trial had no role in the design, conduct or management of the study, in data collection, analysis or interpretation of data, or in the preparation of the present manuscript and decision to publish. MCV is Tier 1 Canada Research Chair in Genomics Applied to Nutrition and Metabolic Health.

Author information

Authors and Affiliations

Authors

Contributions

AT designed the experiments and coordinated the project. LB, OL and SL collected the data and followed the patients; MN collected the data and reviewed the article; ABM and JTM performed the statistical analysis and drafted the manuscript. DR, MCV and AT reviewed the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Marie-Claude Vohl.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the ethics committee of IUCPQ. All participants provided written informed consent prior to participation in accordance with the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

AT and LB received research funding from Johnson & Johnson for the present study in conjunction with a team grant from the Canadian Institutes of Health Research. They also receive funding from Medtronic and GI Windows for studies unrelated to the present article. AT received consulting fees from Novo Nordisk, Eli Lilly and Bausch Health. The remaining authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Table S1. Participants characteristics. Figure S1. Subcutaneous adipose tissue morphology before and after the bariatric surgery.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Bouchard-Mercier, A., de Toro-Martín, J., Nadeau, M. et al. Molecular remodeling of adipose tissue is associated with metabolic recovery after weight loss surgery. J Transl Med 20, 283 (2022). https://doi.org/10.1186/s12967-022-03485-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12967-022-03485-6

Keywords

  • Whole genome
  • Transcriptomic
  • Methylomic
  • Obesity
  • Bariatric surgery
  • Remission
  • Type 2 diabetes
  • Dyslipidemia
  • Subcutaneous adipose tissue