Next Article in Journal
GDF-15 Predicts In-Hospital Mortality of Critically Ill Patients with Acute Kidney Injury Requiring Continuous Renal Replacement Therapy: A Multicenter Prospective Study
Previous Article in Journal
Clinical Relevance of Immersive Virtual Reality in the Assessment and Treatment of Addictive Disorders: A Systematic Review and Future Perspective
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Changes in the Proteome Profile of People Achieving Remission of Type 2 Diabetes after Bariatric Surgery

1
The School of Medicine and Manchester Academic Health Sciences Centre, Manchester University, Manchester M13 9PL, UK
2
Department of Endocrinology, Diabetes and Metabolism, Salford Royal Foundation Trust, Salford M6 8HD, UK
3
Stoller Biomarker Discovery Centre, Division of Cancer Sciences, School of Medical Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester M13 9PL, UK
4
Manchester National Institute for Health Research Biomedical Research Centre, Manchester M13 9WL, UK
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this manuscript.
J. Clin. Med. 2021, 10(16), 3659; https://doi.org/10.3390/jcm10163659
Submission received: 22 July 2021 / Revised: 17 August 2021 / Accepted: 17 August 2021 / Published: 18 August 2021

Abstract

:
Bariatric surgery (BS) results in metabolic pathway recalibration. We have identified potential biomarkers in plasma of people achieving type 2 diabetes mellitus (T2DM) remission after BS. Longitudinal analysis was performed on plasma from 10 individuals following Roux-en-Y gastric bypass (n = 7) or sleeve gastrectomy (n = 3). Sequential window acquisition of all theoretical fragment ion spectra mass spectrometry (SWATH-MS) was done on samples taken at 4 months before (baseline) and 6 and 12 months after BS. Four hundred sixty-seven proteins were quantified by SWATH-MS. Principal component analysis resolved samples from distinct time points after selection of key discriminatory proteins: 25 proteins were differentially expressed between baseline and 6 months post-surgery; 39 proteins between baseline and 12 months. Eight proteins (SHBG, TF, PRG4, APOA4, LRG1, HSPA4, EPHX2 and PGLYRP) were significantly different to baseline at both 6 and 12 months post-surgery. The panel of proteins identified as consistently different included peptides related to insulin sensitivity (SHBG increase), systemic inflammation (TF and HSPA4—both decreased) and lipid metabolism (APOA4 decreased). We found significant changes in the proteome for eight proteins at 6- and 12-months post-BS, and several of these are key components in metabolic and inflammatory pathways. These may represent potential biomarkers of remission of T2DM.

1. Introduction

Obesity is a chronic disease which in many cases requires complex management, while being recognized as the fastest growing problem affecting public health worldwide [1,2,3].
Obesity often has multiple predisposing and precipitating factors and is itself associated with a high mortality rate and with comorbid conditions such as type 2 diabetes mellitus (T2DM), metabolic syndrome, hypertension, dyslipidemia, several cancers [4], premature cell ageing [5], sleep apnea and osteoarthritis [6,7]. Furthermore, obesity is a risk factor for dementia and mild cognitive impairment [8].
Bariatric surgery (BS) has been shown not only to be an effective therapy for weight loss, but also improves a variety of metabolic parameters conferring protection from cardiovascular (CVD) and other diseases [9,10]. Immediate post-operative weight-independent effects combined with weight loss can result in remission of T2DM in up to 80% of patients [10]. Remission is observed also in patients undergoing gastrointestinal surgery similar in design to BS for other reasons such as gastric cancer, which reportedly has better T2DM remission rates [11]. Normalization of glucose metabolism occurs through mechanisms that have not yet been fully determined [11,12]. However, different studies have shown that reductions in waist circumference after BS are associated with a greater probability of T2DM remission [13,14,15].
The study of proteins and their functions has been important in helping investigators to decipher the cellular mechanisms that relate to particular phenotypes [16,17], in addition to accelerating the search for biomarkers for predicting disease pathogenesis and outcomes.
Sequential window acquisition of all theoretical fragment ion spectra mass spectrometry (SWATH-MS) is a data independent acquisition method, which aims to enable larger clinical studies by offering relative quantification across multiple samples [18]. It allows a complete and permanent recording of all fragment ions of the detectable peptide precursors present in a biological sample [19].
In order to have a wider view of the changes that occur after weight loss induced by BS as they relate to remission of T2DM and begin to understand if there are potential early biomarkers for T2DM remission, we set out to identify the changes in the plasma proteome of T2DM patients following BS.

2. Materials and Methods

Ten individuals (2 males and 8 females) were selected from the subset of those who had achieved remission of T2DM following Roux-en-Y gastric bypass (n = 7) or sleeve gastrectomy (n = 3) BS as participants in a larger prospective study (please see Adam et al. for further details of methodology [20]). All were of Caucasian ethnicity. Patients were selected on the basis they had achieved remission of T2DM at 12 months. Remission was defined as reduction of glycosylated haemoglobin (HbA1c) below 42 mmol/mol with cessation of all anti-diabetic therapy at 12 months [21]. SWATH-MS proteomics was carried out on 29 plasma samples with the following sample sub-groups: baseline samples, pre-bariatric surgery (taken at 4 months prior to BS (n = 10); post-bariatric surgery samples, 6 months follow up (n = 10); post-bariatric surgery samples, 12 months follow up (n = 9).

2.1. Sample Preparation for SWATH-MS

Plasma samples (10 µL) were immunodepleted using top 12 abundant protein depletion spin columns from Pierce Biotechnology, UK. Depleted plasma was buffer exchanged and concentrated using Amicon centrifugal filters. Total protein concentration was determined by microplate format (bicinchoninic acid assay). Each sample was reduced, alkylated and digested prior to lyophilisation. More details available in McGurk et al., 2020 [22].

2.2. Liquid Chromatography Mass Spectrometry (LC-MS) Instrument Analyses

Samples were analysed by SWATH-MS with a micro-flow LC-MS system comprising an Eksigent nanoLC 400 autosampler and an Eksigent nanoLC 425 pump coupled to a SCIEX 6600 Triple-TOF mass spectrometer with a DuoSpray Ion Source. The system was configured for a trap-elute elute analysis in which sample was injected from the autosampler (8 °C) onto a trap column (YMC- Triart C18; length: 5 mm; ID: 0.5 mm; particle size: 3 µm; pore size: 120 Å) with loading buffer mobile phase (10 µL min-1, 3 min, 2% acetonitrile, 0.1% formic acid) then eluted through an analytical column (YMC-Triart C18; length: 150 mm; ID: 0.3 mm; particle size: 3 µm; pore size: 120 Å; 30 °C) with the required analytical gradient into the mass spectrometer source. The system was controlled by Analyst software v1.7.1 and Eksigent control software v4.2 (SCIEX, Nieuwerkerk aan den IJssel, The Netherlands).
Peptides were eluted with an analytical gradient (pick-up 10 µL, 5 µL/min, 68 min) and analysed using a mass spectrometer method with a total duty cycle of 2.59 s comprising a TOF MS1 scan (m/z 400–1250, 250 ms) followed by 100 SWATH-MS scans (m/z 100–1500, 20 ms) with variable m/z isolation widths, collision energy and collision energy spread. SWATH maps generated via microflow LC on line to mass spectrometry for 68 min (as described in reference [23]). SWATH-MS data were searched using openSWATH (Version 2.0.0) against the Stoller human serum spectral library. Peptide matches were scored using pyProphet (Version 0.18.3) and search results were aligned using the feature alignment script from MSproteomicstools. Downstream analysis was performed in R using the bioconductor (release 3.5) packages SWATH2Stats and MSstats.

2.3. Determination of Glycosylated Haemoglobin, Glucose and Insulin

HbA1c was measured on an Hb9210 Premier autoanalyser (boranate affinity and high-performance liquid chromatography (Menarini Diagnostics, Wokingham, Berkshire, UK). Glucose and insulin were measured using Abcam ELISA kits (Abcam, Cambridge, UK) and using the HOMA2 calculator (https://www.dtu.ox.ac.uk/homacalculator/ accessed on 2 May 2020) were used to calculate the homeostatic model assessment of insulin resistance (HOMA-IR) and beta cell function (HOMA-B) [24].

2.4. Ethics

Informed consent was obtained from each participant before recruitment. This research adhered to the tenets of the Declaration of Helsinki. All participants provided written informed consent. The study was approved by the Greater Manchester Research Ethics Committee (REC No:11/NW/0731, IRAS ID: 85208)

3. Statistical Analysis

Statistical analysis of anthropometric and biochemical measurements was carried out on SPSS for Mac (Version 23.0, IBM Corporation, New York, NY, USA). Normality was determined by using the Shapiro-Wilk test and by visualising the histogram and normal Q-Q plot. To assess within and between group differences we used one-way analysis of variance for parametric variables and Friedman’s test for non-parametric variables. A significant p value was considered to be <0.05 (post hoc—Tukey). Statistical analysis was performed using all proteins quantified at both time points in 3 or more individuals; 6 months vs. baseline = 314 proteins; 12 months vs. baseline = 317 proteins; 6 months vs. 12 months = 309 proteins.
Principal component analysis (PCA) and heatmaps with hierarchical clustering (Euclidean distance) were plotted to assess separation between time points. Differential expression analysis was performed using the LIMMA package in R (version 3.46.0) on protein fold changes calculated for each of the 10 individuals at 6- or 12-months post-surgery relative to pre-surgery samples, taking a p-value ≤ 0.05 as significant.

4. Results

Longitudinal sampling is of value in biomarker analyses. Thus, the anthropometric and biochemical measurements were taken at baseline, 6- and 12-months as shown in Table 1 and SWATH-MS data described below. Two men and 8 women with a mean age of 53 years were included in the analysis. Seven individuals underwent Roux-en-Y gastric bypass and three underwent sleeve gastrectomy.
BS resulted in a significant reduction in body mass index (BMI) after 6- and 12-months (p < 0.0001), systolic and diastolic blood pressure after 12 months (p < 0.05) but not 6 months, HbA1c after 6- and 12-months (p < 0.0001). Fasting plasma glucose declined significantly at both 6 and 12 months (p < 0.05). HOMA2-IR declined at 6 and 12 months but did not reach significance (p = 0.57) and HOMA2-B (%) increased at 6 and 12 months (p = 0.10) but this also did not reach significance (p = 0.14) (Table 1).
All 10 individuals underwent remission of T2DM achieving HbA1c less than 42 off all diabetes medications after 12 months. Following SWATH-MS analysis, principal component analysis (PCA) revealed very little separation between the different time points when including all 467 quantified proteins (Figure 1). Using differential expression analysis, 25 proteins were significantly different between pre-surgery and 6 months post-surgery samples (Table 2, Figure 2A). These include three apolipoproteins (APOC3, APOM and APOA4), sex hormone binding globulin (SHBG), serotransferrin (TF) and angiotensinogen (AGT). Hierarchal clustering with log2 protein abundances scaled across all pre-surgery and 6 months post-surgery samples shows that, with just the 25 significant proteins, it is possible to discriminate 6 months post-surgery from pre-surgery samples (Figure 2B). PCA using log2 protein abundance ratios for each of the 25 discriminatory proteins across all 10 individuals, again shows that the significant proteins largely separate pre-surgery from 6-month post-surgery samples, with the 12-month post-surgery samples appearing more similar to the 6-month post-surgery samples (Figure 2C). Gene ontology (GO) term enrichment analysis of the 25 discriminatory proteins showed significant enrichment for proteins involved in 55 main biological processes, among them and most important to this context: response to stimulus, leukocyte mediated immunity, plasma lipoprotein remodelling and protein-lipid complex remodelling pathways. All pathways displayed had a Benjamini adjusted p-value ≤ 0.05 and a minimum of 3 proteins contributing to the enrichment. The discriminatory proteins with the biggest fold enrichment (FE) changes were APOA4, APOC3, APOM and AGT (FE = 96.35; p = 0.002).
When comparing pre-surgery samples to those collected 12 months post-surgery there was a greater degree of difference, with 39 proteins significantly different (Table 3, Figure 3A). These include two apolipoproteins (APOA1 and APOA4), SHBG, fibronectin (FN1) and haptoglobulin (HP). With those 39 significant proteins, hierarchical clustering of scaled log2 protein abundances showed discrimination between pre-surgery and 12 months post-surgery samples (Figure 3B). PCA with log2 protein abundance ratios for the 39 discriminatory proteins shows that a greater degree of separation was achieved at 12 months post-surgery, with samples taken 6 months post-surgery falling in between the pre-surgery and 12-month samples based on these 39 discriminatory proteins (Figure 3C). GO term enrichment analysis of the 39 discriminatory proteins showed enrichment for proteins involved in 86 significant biological processes, most of them involving healing, blood coagulation and immune activation pathways. All pathways displayed had a Benjamini adjusted p-value ≤ 0.05 and a minimum of 3 proteins contributing to the enrichment. The proteins showing the biggest FE were: Antithrombin-III (SERPINC1), Plasma protease C1 inhibitor (SERPING1), Alpha-2-macroglobulin (A2M) and FN1 (FE = 69.34; p = 0.0003).
The relative quantification information acquired using SWATH-MS for 467 proteins in a longitudinal fashion has a degree of power as temporal analyses in humans provide greater surety of biomarker status. Relative to pre-surgery samples, eight proteins were significantly different at both 6- and 12-months post-surgery (Figure 4). Of those eight proteins, three showed increased expression after BS (FC: fold change baseline vs. 12 months), p value); sex hormone binding globulin (SHBG) (1.95, p < 0.01), leucine-rich alpha-2-glycoprotein (LRG1) (0.59, p < 0.05) and N-acetylmuramoyl-L-alanine amidase (PGLYRP2) (0.43, p < 0.05), whilst the remaining five showed decreased expression; serotranferrin (TF) (−0.78, p < 0.01), proteoglycan 4 (PRG4) (−0.78, p < 0.05), Apolipoprotein A4 (APOA4) (−1.38, p < 0.05), heat shock protein 4 (HSPA4) (−0.38, p < 0.05), bifunctional epoxide hydrolase 2 (EPHX2) (Table 4). The greatest fold change was seen for SHBG approaching a two-fold elevation after BS. The correspondent protein names and gene symbol for each code used to build the figures (volcano plots and heatmaps) can be found in the supplementary material.

5. Discussion

Using the technique of SWATH-MS to generate proteomic maps, we have shown it is possible to separate samples collected pre- and post bariatric surgery on the basis of the plasma proteome with data from a relatively small sample size. We found a significant change in plasma protein levels for a number of metabolically relevant proteins from pre-BS to 6- and 12-months post-surgery with eight proteins showing change at both 6- and 12-months vs. pre-surgery levels. The relevant proteins common to both post-surgery periods analysed are TF, APOA4, HSPA4, LRG1, PGLYRP2, SHBG, EPHX2 and PRG4.
Greater resolution between baseline and post-bariatric surgery samples was noted when proteins with significant fold changes were solely included in the analysis, as a means of discovering potential biomarkers; Figure 2 shows the PCA of the significant proteins between baseline and 12 months and Figure 3 shows the same between baseline and 6 months. The clearest separation occurs between baseline and 12 months which coincides with the onset of the time period at which often the greatest weight loss is observed with BS [25]. We previously described proteomic changes in relation to a lifestyle change intervention in non-diabetic hyperglycaemia and also showed that the baseline levels of certain proteins such as insulin-like growth factor (IGF)-II and vitamin D binding protein were predictive of more or less weight loss with that lifestyle change intervention [26].
Previous work employing SWATH-MS to study biomarkers in patients with impaired glucose tolerance who lost weight in response to a diet and exercise programme and in whom HbA1c resolved identified the following major differentiating proteins: Insulin-like growth factor 2 (IGF-II), Retinol binding protein 4 (RBP4), Fetuin-A (FetA) -Zinc-α2- glycoprotein (ZA2G), Visfatin (NAMPT), fatty acid synthase (FAS) and vitamin D binding protein (VDR) [26]. The disparity between this study and ours suggests a differential process and a potentially specific set of biomarkers for remission of T2DM in patients undergoing BS compared to those achieving weight loss through other structured means. No difference in circulating IGF binding proteins was picked up by the SWATH-MS analysis in the present study. HOMA2-IR as a measure of insulin resistance decreased and HOMA2-B as measure of pancreatic beta cell function increased numerically; this did not reach statistical significance.
Previously Varela-Rodríguez et al. showed, using proteomic analysis, that BS remodels subcutaneous adipose tissue function to influence specific molecular mechanisms with lower inflammation, increased uptake of glucose, higher insulin sensitivity, higher de novo lipogenesis, increased mitochondrial function and decreased adipocyte size [27]. Our results support these findings as the cluster of relevant proteins showing a change after 6 months of BS are involved in biological process such as protein-lipid remodelling, leukocyte mediate immunity and response to stimulus pathways; and after 12 months, proteins involved in healing, blood coagulation and immune activation.
The panel of eight proteins identified as discriminating between the groups, and common at both 6- and 12-months post-surgery, included proteins known to be related to insulin sensitivity and glucose intolerance (TF, SHBG and PRG4) [28,29,30], relate to systemic inflammation (HSPA4 and APOA4) [31,32] ,also those involved in the immune response (LRG1 and PGLYRP2) [33,34,35,36] and cholesterol function (EPHX2) [37,38]. We have highlighted the importance of each protein in the context of obesity, BS and T2DM below.
LRG1 is a circulating protein first discovered in 1977 [39] whose function has remained elusive for many years. Some evidence has suggested that it may function, in part, as an angiogenesis factor [40]. It has been observed as a marker of inflammation in the sputum of asthmatic patients [41], in the serum of patients with ulcerative colitis [42] and autoimmune disease respectively [43]. Shirai et al. showed that it is often up-regulated during the acute phase response [44]. Studies have also found LRG1 to be higher in patients with T2DM and peripheral arterial disease [45]. Notwithstanding previous studies, our results demonstrate a clear increase in LRG1 expression both 6- and 12-months after surgery. This is a similar finding to that of Roriguez-Rivera et al. [46] one year after surgery and Pek et al. who prospectively studied 231 morbidly obese patients undergoing BS and found LRG1 levels to significantly increase just 1.5 months after surgery [47]. These results suggest that inflammation is probably not the sole inducer of LRG1 expression and further research into the role of this glycoprotein is required.
TF levels declined after BS surgery in our cohort, which may reflect the reduced inflammatory state manifested in tandem with weight loss after BS. TF is a part of the innate immune system and acts as a marker of inflammation in patients with T2DM [48]. Furthermore, pre-operative deficiencies of iron and other trace elements are prevalent amongst individuals with obesity undergoing surgery and are often replaced post-surgery [49], which offers an explanation as to the maintained significant decline in this protein.
HSPA4 is largely thought to function as a cytosolic chaperone involved in facilitating protein folding, degradation, complex assembly and translocation. Significantly higher levels of HSPA4 have been reported in T2DM of longer duration compared to newly diagnosed [50]. Nakhiayani et al. showed that HSPA4 was associated with the inhibition of nitric oxide production in individuals with T2DM [51], confirmed that it correlated with CRP [51] whilst Morteza et al. found that the odds ratio for its predictive value for micro-albuminuria in T2DM was highly significant [52]. Garamyolgyi et al. showed that HSPA4 levels correlated well with HbA1c in gestational diabetes [53]. Our work demonstrates a significant reduction in HSPA4 which correlated with percentage change in IL6 between baseline and 12 months (r = 0.87, p < 0.001), this is in keeping with previous work by Styger et al. who demonstrated reduction in both serum and liver HSPA4 in a murine model after BS [54] and also is in agreement with the narrative that HSPA4 appears to be induced by states of chronic inflammation [55]. Some have suggested targeting these proteins as therapy for T2DM [56].
In our cohort, APOA4 was noted to significantly and consistently reduce at both time points after BS which may appear intuitive since this lipid-binding protein is primarily synthesized in the small intestine, which is bypassed during this type of surgery [57]. In animal models circulating Apoa4 appear to confer some protection against diabetes [58,59] and atherosclerosis. However, its role in humans has not fully been delineated. Previous studies conflict in their results surrounding this lipoprotein with some reporting no change after BS [60] whilst others have increases in APOA4 levels post-surgery [46,61]. Rao et al. reported that BS resulted in a decrease in HOMA-IR in patients after RYGB and higher baseline APOA4 levels correlated with the decrease in HOMA-IR. Another recent study, utilizing proteomic methods, analysed a large number of obese individuals who were subject to low calorie diet finding that, amongst others, APO A4 levels declined after the eight week intervention [62]. We speculate that APOA4 may be a physiologic compensation to insulin resistance that resolves after BS. In the circulation, APOA4 is present on chylomicron remnants, high-density lipoproteins, and in lipid-free form. APOA4 is involved in a myriad of physiological processes such as lipid absorption and metabolism, anti-atherosclerosis, platelet aggregation and thrombosis, glucose homeostasis and food intake. APOA4 deficiency is associated with atherosclerosis and diabetes, which renders it as a potential therapeutic target for treatment of these diseases [63].
PGLRYP2 [64] is an enzyme that breaks down glycopeptides and can break down the cell walls of bacteria [65] giving it a principal role in fighting bacterial infection. It is expressed constitutively from hepatocytes [66] and proteomic studies have previously suggested a possible utility in diagnosing sepsis in critically ill patients [67]. Expression can be induced in skin keratinocytes after exposure to certain bacteria [68] as well as in intestinal and oral epithelial cells [69]. Our study demonstrates a significant increase in fold change of this protein at both 6- and 12-months post-surgery with correlations with percentage change in BMI at 6 months (p < 0.05) but not 12 months. We speculate that this change may result from the resolution of non-alcoholic hepatic steatosis (NAHS) in our patients. Previous work has shown that the PGLRYP2 gene is hypermethylated in patients with NAHS [70] and immune-competent mice who overexpress PGLRYP2 display significantly enhanced anti-cancer immune responses against hepatocellular carcinoma (HCC) [35]. Indeed, resolution of both NAHS and lower incidence of HCC are recorded as one of the beneficial effects of BS [71,72].
SHBG is produced and secreted by the liver and acts to regulate the bioavailability of sex steroids. [73]. Low serum SHBG concentrations are associated with metabolic syndrome [74], T2DM [75] and increased risk for cardiovascular problems. [76] In obesity, SHBG levels are low and associated with pro-inflammatory cytokines and hepatic steatosis [74]. Our results showed an increase in SHBG levels after BS consistent with results found previously in men and women after BS [77,78].
EPHX2 is a soluble enzyme that hydrolyses epoxyeicoasatrienoic acids (EETs) to their inactive diols [79] EETS play important roles in vasodilation [80], lipid metabolism [81] improve insulin sensitivity [82] and have both anti-inflammatory [83] and analgesic properties [84]. Gene deletion and pharmacological inhibition of EPHX2 exhibit an increase in insulin sensitivity in a T2DM rodent model [85]. Furthermore, its variants are likely to be associated with neuropathy in T2DM patients [86]. EPHX2 circulating levels are elevated in obese individuals and its expression is attenuated after physical activity in both adipose tissue and peripheral blood mononuclear cells [87]. Higher EPHX2 levels are significantly associated with risk of incident CVD [88]. We observed a consistent reduction of EPHX2 at both 6- and 12- months post-BS. The inhibition of EPHX2 has been suggested as previous therapeutic target to tackle CVD and this requires further study [89].
PRG4 was first identified as a component of the extracellular matrix in synovial fluid implicated in shear force reduction in cartilage [90]. However, recent evidence suggests that it can act as a “buffer” by effecting downstream signalling during an inflammatory response [91]. Furthermore, chronic inflammatory states cause post-transcriptional modifications resulting in dysfunctional PRG4 upsetting homeostasis [92]. PRG4 levels declined significantly at both time points which is a consistent finding by previous groups studying the effect of weight loss on the proteome [93,94] and has previously been found to correlate with insulin sensitivity [94]. Knockout PRG4 mice exhibit better glucose handling when fed a high-fat diet and are protected from hepatic steatosis and white adipose tissue inflammation [30]. The regulation and roles of PRG4 in obesity are largely unknown and require further work to be delineated.

6. Strengths and Limitations

We are showing here preliminary data from a hypothesis-driven study. Our study utilises data on only 10 people, mostly female, which is a limitation. Nevertheless, we have described significant changes in the proteome in relation to remission of T2DM. We have identified significant differences in the plasma levels of a number of metabolically relevant proteins after BS in relation to remission of T2DM using a technique that has high reproducibility and sensitivity, in three time points (baseline, 6- and 12-months post-surgery). In this respect, it offers a potential set of biomarkers for remission that require further investigation, which is now planned. An additional limitation is that we have only here described changes in protein levels in plasma rather than in adipose or other tissues. However, biomarkers are generally measured in biofluids. Finally, our study did not compare remitters with non-remitters, which would add further value in identifying biomarkers relevant to T2DM remission. Further analysis using different and larger cohorts are necessary to replicate our findings.

7. Conclusions

Using SWATH-MS we found significant changes in the proteome for eight proteins from pre-BS to 6/12 and 12/12 post-BS. This panel of proteins identified as consistently different, included peptides related to insulin sensitivity (SHBG increase), systemic inflammation (TF and HSPA4—both decreased) and lipid metabolism (APOA4 decreased). Several of these are key components in metabolic, immune system and inflammatory pathways. Therefore, these protein changes merit further exploration as potential marker signatures in patients who undergo remission of T2DM.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/jcm10163659/s1.

Author Contributions

Conceptualization, A.H.H. and J.M.G.; methodology, I.B.-J., A.E.C., B.G. and A.D.W.; software, A.E.C.; validation, I.B.-J. and A.E.C. and A.D.W.; formal analysis, Z.I. and H.A.F.; investigation, Z.I. and H.A.F.; resources, R.P.D., A.H.H., H.S. and A.S.; data curation, R.P.D., A.H.H., H.S., D.H. and A.S.; writing—original draft preparation, Z.I. and H.A.F.; writing—review and editing, Z.I., H.A.F., D.H. and A.H.H.; visualization, A.H.H.; supervision, A.H.H., R.P.D. and H.S.; project administration, A.H.H., R.P.D., A.S. and H.S.; funding acquisition, A.H.H. and J.M.G. All authors have read and agreed to the published version of the manuscript.

Funding

SWATH-MS was supported with equipment grants from Bloodwise (19007) and Medical Research Council (MR/M008959/1). This work was supported by the CRUK Manchester Centre award (C5759/A25254). ADW is supported by the NIHR Manchester Biomedical Research Centre.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board (or Ethics Committee) of Salford Royal NHS Foundation Trust (NCARC NCA003 04/2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

Any requests for data extracts will be considered by Adrian Heald.

Acknowledgments

We acknowledge support from National Institute for Health Research (NIHR) Manchester Biomedical Research Centre and the patients who took part in this study.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

A2MAlpha-2-macroglobulin
AGTAngiotensinogen
APOA1Apolipoprotein A1
APOA4Apolipoprotein A4
APOC3Apolipoprotein C3
APOMApolipoprotein M
BMIBody Mass Index
BSBariatric Surgery
CRPC—reactive protein
CVDCardiovascular disease
EETsEpoxyeicoasatrienoic acids
EPHX2Bifunctional epoxide hydrolase 2
FEFold Enrichment
FN1Fibronectin
GOGene Ontology
HbA1cglycosylated haemoglobin
HCCHepatocellular carcinoma
HOMA-BHomeostatic model assessment of beta cell function
HOMA-IRHomeostatic model assessment of insulin resistance
HPHaptoglobulin
HSPA4Heat shock 70 kDa protein 4
IGF-IIInsulin-like growth factor II
LC-MSLiquid Chromatography Mass spectrometry
LRG1Leucine-rich alpha-2-glycoprotein
NAHSNon-Alcoholic Hepatic Steatosis
PCAPrincipal component analysis
PGLYRP2N-acetylmuramoyl-L-alanine amidase
PRG4Proteoglycan 4
SERPINC1Antithrombin-III
SERPING1Plasma protease C1 inhibitor
SHBGSex Hormone Binding Globulin
SWATH-MSSequential window acquisition of all theoretical fragment ion spectra Mass Spectrometry
T2DMType 2 Diabetes Mellitus
TFSerotransferrin

References

  1. Jastreboff, A.M.; Kotz, C.M.; Kahan, S.; Kelly, A.S.; Heymsfield, S.B. Obesity as a disease: The obesity society 2018 position statement. Obesity 2019, 27, 7–9. [Google Scholar] [CrossRef] [Green Version]
  2. Bray, G.A.; Heisel, W.E.; Afshin, A.; Jensen, M.D.; Dietz, W.H.; Long, M.; Kushner, R.F.; Daniels, S.R.; Wadden, T.A.; Tsai, A.G. The science of obesity management: An endocrine society scientific statement. Endocr. Rev. 2018, 39, 79–132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Heymsfield, S.B.; Wadden, T.A. Mechanisms, pathophysiology, and management of obesity. N. Engl. J. Med. 2017, 376, 254–266. [Google Scholar] [CrossRef]
  4. Lin, C.; Cai, X.; Yang, W.; Lv, F.; Nie, L.; Ji, L. The Body Weight Alteration and Incidence of Neoplasm in Patients with Type 2 Diabetes: A Meta-Analysis of Randomized Controlled Trials. Front. Endocrinol. 2020, 11, 1001. [Google Scholar] [CrossRef] [PubMed]
  5. Ogrodnik, M.; Zhu, Y.; Langhi, L.G.; Tchkonia, T.; Krüger, P.; Fielder, E.; Victorelli, S.; Ruswhandi, R.A.; Giorgadze, N.; Pirtskhalava, T. Obesity-induced cellular senescence drives anxiety and impairs neurogenesis. Cell Metab. 2019, 29, 1061–1077.e8. [Google Scholar] [CrossRef] [Green Version]
  6. Flegal, K.M.; Graubard, B.I.; Williamson, D.F.; Gail, M.H. Cause-specific excess deaths associated with underweight, overweight, and obesity. Obstet. Gynecol. Surv. 2008, 63, 157–159. [Google Scholar] [CrossRef]
  7. Apovian, C.M. Obesity: Definition, comorbidities, causes, and burden. Am. J. Manag. Care 2016, 22, s176–s185. [Google Scholar] [PubMed]
  8. Dye, L.; Boyle, N.B.; Champ, C.; Lawton, C. The relationship between obesity and cognitive health and decline. Proc. Nutr. Soc. 2017, 76, 443–454. [Google Scholar] [CrossRef] [Green Version]
  9. Ammori, B.J.; Skarulis, M.C.; Soran, H.; Syed, A.A.; Eledrisi, M.; Malik, R.A. Medical and surgical management of obesity and diabetes: What’s new? Diabet. Med. 2020, 37, 203–210. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Iqbal, Z.; Adam, S.; Ho, J.H.; Syed, A.A.; Ammori, B.J.; Malik, R.A.; Soran, H. Metabolic and cardiovascular outcomes of bariatric surgery. Curr. Opin. Lipidol. 2020, 31, 246–256. [Google Scholar] [CrossRef]
  11. Peng, D.; Cheng, Y.X.; Zhang, W. Does Roux-en-Y Construction Really Bring Benefit of Type 2 Diabetes Mellitus Remission after Gastrectomy in Patients with Gastric Cancer? A Systematic Review and Meta-Analysis. Diabetes Ther. 2020, 11, 2863–2872. [Google Scholar] [CrossRef]
  12. Frikke-Schmidt, H.; O’Rourke, R.W.; Lumeng, C.N.; Sandoval, D.A.; Seeley, R.J. Does bariatric surgery improve adipose tissue function? Obes. Rev. 2016, 17, 795–809. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Purnell, J.Q.; Selzer, F.; Wahed, A.S.; Pender, J.; Pories, W.; Pomp, A.; Dakin, G.; Mitchell, J.; Garcia, L.; Staten, M.A. Type 2 diabetes remission rates after laparoscopic gastric bypass and gastric banding: Results of the longitudinal assessment of bariatric surgery study. Diabetes Care 2016, 39, 1101–1107. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Kadera, B.E.; Lum, K.; Grant, J.; Pryor, A.D.; Portenier, D.D.; DeMaria, E.J. Remission of type 2 diabetes after Roux-en-Y gastric bypass is associated with greater weight loss. Surg. Obes. Relat. Dis. 2009, 5, 305–309. [Google Scholar] [CrossRef] [PubMed]
  15. Hamza, N.; Abbas, M.H.; Darwish, A.; Shafeek, Z.; New, J.; Ammori, B.J. Predictors of remission of type 2 diabetes mellitus after laparoscopic gastric banding and bypass. Surg. Obes. Relat. Dis. 2011, 7, 691–696. [Google Scholar] [CrossRef]
  16. Aebersold, R.; Mann, M. Mass-spectrometric exploration of proteome structure and function. Nature 2016, 537, 347–355. [Google Scholar] [CrossRef]
  17. Tanner, S.D.; Baranov, V.I.; Ornatsky, O.I.; Bandura, D.R.; George, T.C. An introduction to mass cytometry: Fundamentals and applications. Cancer Immunol. Immunother. 2013, 62, 955–965. [Google Scholar] [CrossRef]
  18. Zhu, X.; Chen, Y.; Subramanian, R. Comparison of information-dependent acquisition, SWATH, and MSAll techniques in metabolite identification study employing ultrahigh-performance liquid chromatography–quadrupole time-of-flight mass spectrometry. Anal. Chem. 2014, 86, 1202–1209. [Google Scholar] [CrossRef]
  19. Gillet, L.C.; Navarro, P.; Tate, S.; Röst, H.; Selevsek, N.; Reiter, L.; Bonner, R.; Aebersold, R. Targeted data extraction of the MS/MS spectra generated by data-independent acquisition: A new concept for consistent and accurate proteome analysis. Mol. Cell. Proteom. 2012, 11. [Google Scholar] [CrossRef] [Green Version]
  20. Adam, S.; Azmi, S.; Ho, J.H.; Liu, Y.; Ferdousi, M.; Siahmansur, T.; Kalteniece, A.; Marshall, A.; Dhage, S.S.; Iqbal, Z. Improvements in Diabetic Neuropathy and Nephropathy after Bariatric Surgery: A Prospective Cohort Study. Obes. Surg. 2020, 31, 554–563. [Google Scholar] [CrossRef]
  21. Brethauer, S.A.; Kim, J.; El Chaar, M.; Papasavas, P.; Eisenberg, D.; Rogers, A.; Ballem, N.; Kligman, M.; Kothari, S. Standardized Outcomes Reporting in Metabolic and Bariatric Surgery. Obes. Surg. 2015, 25, 587–606. [Google Scholar] [CrossRef] [PubMed]
  22. McGurk, K.A.; Dagliati, A.; Chiasserini, D.; Lee, D.; Plant, D.; Baricevic-Jones, I.; Kelsall, J.; Eineman, R.; Reed, R.; Geary, B. The use of missing values in proteomic data-independent acquisition mass spectrometry to enable disease activity discrimination. Bioinformatics 2020, 36, 2217–2223. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Walker, D.E.; Yang, F.; Carver, J.; Joe, K.; Michels, D.A.; Yu, X.C. A modular and adaptive mass spectrometry-based platform for support of bioprocess development toward optimal host cell protein clearance. mAbs 2017, 9, 654–663. [Google Scholar] [CrossRef] [PubMed]
  24. Matthews, D.; Hosker, J.; Rudenski, A.; Naylor, B.; Treacher, D.; Turner, R. Homeostasis model assessment: Insulin resistance and β-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia 1985, 28, 412–419. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Magro, D.O.; Geloneze, B.; Delfini, R.; Pareja, B.C.; Callejas, F.; Pareja, J.C. Long-term weight regain after gastric bypass: A 5-year prospective study. Obes. Surg. 2008, 18, 648–651. [Google Scholar] [CrossRef] [PubMed]
  26. Malipatil, N.; Fachim, H.A.; Siddals, K.; Geary, B.; Wark, G.; Porter, N.; Anderson, S.; Donn, R.; Harvie, M.; Whetton, A.D. Data independent acquisition mass spectrometry can identify circulating proteins that predict future weight loss with a diet and exercise programme. J. Clin. Med. 2019, 8, 141. [Google Scholar] [CrossRef] [Green Version]
  27. Varela-Rodríguez, B.M.; Juiz-Valiña, P.; Varela, L.; Outeiriño-Blanco, E.; Bravo, S.B.; García-Brao, M.J.; Mena, E.; Noguera, J.F.; Valero-Gasalla, J.; Cordido, F. Beneficial effects of bariatric surgery-induced by weight loss on the proteome of abdominal subcutaneous adipose tissue. J. Clin. Med. 2020, 9, 213. [Google Scholar] [CrossRef] [Green Version]
  28. Liu, Y.; Wang, H.Y.; Shao, J.; Zaro, J.L.; Shen, W.C. Enhanced insulin receptor interaction by a bifunctional insulin-transferrin fusion protein: An approach to overcome insulin resistance. Sci. Rep. 2020, 10, 7724. [Google Scholar] [CrossRef]
  29. Deswal, R.; Yadav, A.; Dang, A.S. Sex hormone binding globulin-an important biomarker for predicting PCOS risk: A systematic review and meta-analysis. Syst. Biol. Reprod. Med. 2018, 64, 12–24. [Google Scholar] [CrossRef]
  30. Nahon, J.E.; Hoekstra, M.; van Harmelen, V.; Rensen, P.C.; Willems van Dijk, K.; Kooijman, S.; Van Eck, M. Proteoglycan 4 deficiency protects against glucose intolerance and fatty liver disease in diet-induced obese mice. Biochim. Biophys. Acta Mol. Basis Dis. 2019, 1865, 494–501. [Google Scholar] [CrossRef]
  31. Goodman, S.C.; Letra, A.; Dorn, S.; Araujo-Pires, A.C.; Vieira, A.E.; de Souza, L.C.; Yadlapati, M.; Garlet, G.P.; Silva, R.M. Expression of heat shock proteins in periapical granulomas. J. Endod. 2014, 40, 830–836. [Google Scholar] [CrossRef]
  32. Doumatey, A.P.; Zhou, J.; Zhou, M.; Prieto, D.; Rotimi, C.N.; Adeyemo, A. Proinflammatory and lipid biomarkers mediate metabolically healthy obesity: A proteomics study. Obesity 2016, 24, 1257–1265. [Google Scholar] [CrossRef]
  33. Druhan, L.J.; Lance, A.; Li, S.; Price, A.E.; Emerson, J.T.; Baxter, S.A.; Gerber, J.M.; Avalos, B.R. Leucine rich α-2 glycoprotein: A novel neutrophil granule protein and modulator of myelopoiesis. PLoS ONE 2017, 12, e0170261. [Google Scholar] [CrossRef]
  34. Hao, L.; Xie, H.; Zhang, B.; Chen, D.; Wang, S.; Zhang, H.; He, S. LRG1 downregulation in allergic airway disorders and its expression in peripheral blood and tissue cells. J. Transl. Med. 2016, 14, 202. [Google Scholar] [CrossRef] [Green Version]
  35. Yang, Z.; Feng, J.; Xiao, L.; Chen, X.; Yao, Y.; Li, Y.; Tang, Y.; Zhang, S.; Lu, M.; Qian, Y. Tumor-Derived Peptidoglycan Recognition Protein 2 Predicts Survival and Antitumor Immune Responses in Hepatocellular Carcinoma. Hepatology 2020, 71, 1626–1642. [Google Scholar] [CrossRef] [Green Version]
  36. Dziarski, R.; Gupta, D. The peptidoglycan recognition proteins (PGRPs). Genome Biol. 2006, 7, 232. [Google Scholar] [CrossRef] [Green Version]
  37. Shen, L.; Peng, H.; Peng, R.; Fan, Q.; Zhao, S.; Xu, D.; Morisseau, C.; Chiamvimonvat, N.; Hammock, B.D. Inhibition of soluble epoxide hydrolase in mice promotes reverse cholesterol transport and regression of atherosclerosis. Atherosclerosis 2015, 239, 557–565. [Google Scholar] [CrossRef] [Green Version]
  38. Lien, C.C.; Chen, C.H.; Lee, Y.M.; Guo, B.C.; Cheng, L.C.; Pan, C.C.; Shyue, S.K.; Lee, T.S. The phosphatase activity of soluble epoxide hydrolase regulates ATP-binding cassette transporter-A1-dependent cholesterol efflux. J. Cell. Mol. Med. 2019, 23, 6611–6621. [Google Scholar] [CrossRef] [PubMed]
  39. Haupt, H.; Baudner, S. Isolation and characterization of an unknown, leucine-rich 3.1-S-alpha2-glycoprotein from human serum (author’s transl). Hoppe-Seyler’s Z. Physiol. Chem. 1977, 358, 639–646. [Google Scholar] [CrossRef]
  40. Wang, X.; Abraham, S.; McKenzie, J.A.; Jeffs, N.; Swire, M.; Tripathi, V.B.; Luhmann, U.F.; Lange, C.A.; Zhai, Z.; Arthur, H.M. LRG1 promotes angiogenesis by modulating endothelial TGF-β signalling. Nature 2013, 499, 306–311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Honda, H.; Fujimoto, M.; Miyamoto, S.; Ishikawa, N.; Serada, S.; Hattori, N.; Nomura, S.; Kohno, N.; Yokoyama, A.; Naka, T. Sputum leucine-rich alpha-2 glycoprotein as a marker of airway inflammation in asthma. PLoS ONE 2016, 11, e0162672. [Google Scholar] [CrossRef] [PubMed]
  42. Serada, S.; Fujimoto, M.; Terabe, F.; Iijima, H.; Shinzaki, S.; Matsuzaki, S.; Ohkawara, T.; Nezu, R.; Nakajima, S.; Kobayashi, T. Serum leucine-rich alpha-2 glycoprotein is a disease activity biomarker in ulcerative colitis. Inflamm. Bowel Dis. 2012, 18, 2169–2179. [Google Scholar] [CrossRef]
  43. Serada, S.; Fujimoto, M.; Ogata, A.; Terabe, F.; Hirano, T.; Iijima, H.; Shinzaki, S.; Nishikawa, T.; Ohkawara, T.; Iwahori, K. iTRAQ-based proteomic identification of leucine-rich α-2 glycoprotein as a novel inflammatory biomarker in autoimmune diseases. Ann. Rheum. Dis. 2010, 69, 770–774. [Google Scholar] [CrossRef] [PubMed]
  44. Shirai, R.; Hirano, F.; Ohkura, N.; Ikeda, K.; Inoue, S. Up-regulation of the expression of leucine-rich α2-glycoprotein in hepatocytes by the mediators of acute-phase response. Biochem. Biophys. Res. Commun. 2009, 382, 776–779. [Google Scholar] [CrossRef]
  45. Pek, S.L.; Tavintharan, S.; Wang, X.; Lim, S.C.; Woon, K.; Yeoh, L.Y.; Ng, X.; Liu, J.; Sum, C.F. Elevation of a novel angiogenic factor, leucine-rich-α2-glycoprotein (LRG1), is associated with arterial stiffness, endothelial dysfunction, and peripheral arterial disease in patients with type 2 diabetes. J. Clin. Endocrinol. Metab. 2015, 100, 1586–1593. [Google Scholar] [CrossRef] [Green Version]
  46. Rodríguez-Rivera, C.; Pérez-García, C.; Muñoz-Rodríguez, J.R.; Vicente-Rodríguez, M.; Polo, F.; Ford, R.M.; Segura, E.; León, A.; Salas, E.; Sáenz-Mateos, L. Proteomic identification of biomarkers associated with eating control and bariatric surgery outcomes in patients with morbid obesity. World J. Surg. 2019, 43, 744–750. [Google Scholar] [CrossRef]
  47. Pek, S.L.; Cheng, A.K.; Lin, M.X.; Wong, M.S.; Chan, E.Z.; Moh, A.M.; Sum, C.F.; Lim, S.C.; Tavintharan, S. Association of circulating proinflammatory marker, leucine-rich-α2-glycoprotein (LRG1), following metabolic/bariatric surgery. Diabetes Metab. Res. Rev. 2018, 34, e3029. [Google Scholar] [CrossRef]
  48. Golizeh, M.; Lee, K.; Ilchenko, S.; Ösme, A.; Bena, J.; Sadygov, R.G.; Kashyap, S.R.; Kasumov, T. Increased serotransferrin and ceruloplasmin turnover in diet-controlled patients with type 2 diabetes. Free Radic. Biol. Med. 2017, 113, 461–469. [Google Scholar] [CrossRef]
  49. Shipton, M.J.; Johal, N.J.; Dutta, N.; Slater, C.; Iqbal, Z.; Ahmed, B.; Ammori, B.J.; Senapati, S.; Akhtar, K.; Summers, L.K. Haemoglobin and Hematinic Status before and after Bariatric Surgery over 4 years of Follow-Up. Obes. Surg. 2020, 31, 682–693. [Google Scholar] [CrossRef]
  50. Nakhjavani, M.; Morteza, A.; Khajeali, L.; Esteghamati, A.; Khalilzadeh, O.; Asgarani, F.; Outeiro, T.F. Increased serum HSP70 levels are associated with the duration of diabetes. Cell Stress Chaperones 2010, 15, 959–964. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Nakhjavani, M.; Morteza, A.; Asgarani, F.; Khalilzadeh, O.; Ghazizadeh, Z.; Bathaie, S.Z.; Esteghamati, A. The dual behavior of heat shock protein 70 and asymmetric dimethylarginine in relation to serum CRP levels in type 2 diabetes. Gene 2012, 498, 107–111. [Google Scholar] [CrossRef]
  52. Morteza, A.; Nakhjavani, M.; Larry, M.; Nargesi, A.A.; Esteghamati, A. Heat shock protein 70 and albuminuria in patients with type 2 diabetes: A matched case control study. Cell Stress Chaperones 2013, 18, 815–819. [Google Scholar] [CrossRef] [Green Version]
  53. Garamvölgyi, Z.; Prohászka, Z.; Rigó, J.; Kecskeméti, A.; Molvarec, A. Increased circulating heat shock protein 70 (HSPA1A) levels in gestational diabetes mellitus: A pilot study. Cell Stress Chaperones 2015, 20, 575–581. [Google Scholar] [CrossRef] [Green Version]
  54. Stygar, D.; Skrzep-Poloczek, B.; Romuk, E.; Chełmecka, E.; Poloczek, J.; Sawczyn, T.; Maciarz, J.; Kukla, M.; Karcz, K.W.; Jochem, J. The influence of high-fat, high-sugar diet and bariatric surgery on HSP70 and HSP90 plasma and liver concentrations in diet-induced obese rats. Cell Stress Chaperones 2019, 24, 427–439. [Google Scholar] [CrossRef] [Green Version]
  55. Van Eden, W.; Van der Zee, R.; Prakken, B. Heat-shock proteins induce T-cell regulation of chronic inflammation. Nat. Rev. Immunol. 2005, 5, 318–330. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Mulyani, W.R.; Sanjiwani, M.I.; Sandra, I.; Lestari, A.A.; Wihandani, D.M.; Suastika, K.; Saraswati, M.R.; Bhargah, A.; Manuaba, I.B. Chaperone-based therapeutic target innovation: Heat shock protein 70 (HSP70) for Type 2 diabetes mellitus. Diabetes Metab. Syndr. Obes. Targets Ther. 2020, 13, 559. [Google Scholar] [CrossRef] [Green Version]
  57. Qu, J.; Ko, C.W.; Tso, P.; Bhargava, A. Apolipoprotein A-IV: A multifunctional protein involved in protection against atherosclerosis and diabetes. Cells 2019, 8, 319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Ordovas, J.M.; Cassidy, D.; Civeira, F.; Bisgaier, C.; Schaefer, E. Familial apolipoprotein AI, C-III, and A-IV deficiency and premature atherosclerosis due to deletion of a gene complex on chromosome 11. J. Biol. Chem. 1989, 264, 16339–16342. [Google Scholar] [CrossRef]
  59. Cohen, R.D.; Castellani, L.W.; Qiao, J.H.; Van Lenten, B.J.; Lusis, A.J.; Reue, K. Reduced aortic lesions and elevated high density lipoprotein levels in transgenic mice overexpressing mouse apolipoprotein A-IV. J. Clin. Investig. 1997, 99, 1906–1916. [Google Scholar] [CrossRef]
  60. Rao, R.; Roche, A.; Febres, G.; Bessler, M.; Tso, P.; Korner, J. Circulating Apolipoprotein A-IV presurgical levels are associated with improvement in insulin sensitivity after Roux-en-Y gastric bypass surgery. Surg. Obes. Relat. Dis. 2017, 13, 468–473. [Google Scholar] [CrossRef] [Green Version]
  61. Culnan, D.M.; Cooney, R.N.; Stanley, B.; Lynch, C.J. Apolipoprotein A-IV, a putative satiety/antiatherogenic factor, rises after gastric bypass. Obesity 2009, 17, 46–52. [Google Scholar] [CrossRef] [Green Version]
  62. Bruderer, R.; Muntel, J.; Müller, S.; Bernhardt, O.M.; Gandhi, T.; Cominetti, O.; Macron, C.; Carayol, J.; Rinner, O.; Astrup, A. Analysis of 1508 plasma samples by capillary-flow data-independent acquisition profiles proteomics of weight loss and maintenance. Mol. Cell. Proteom. 2019, 18, 1242–1254. [Google Scholar] [CrossRef] [Green Version]
  63. Peng, J.; Li, X.P. Apolipoprotein A-IV: A potential therapeutic target for atherosclerosis. Prostaglandins Other Lipid Mediat. 2018, 139, 87–92. [Google Scholar] [CrossRef]
  64. Wang, Z.M.; Li, X.; Cocklin, R.R.; Wang, M.; Wang, M.; Fukase, K.; Inamura, S.; Kusumoto, S.; Gupta, D.; Dziarski, R. Human peptidoglycan recognition protein-L is an N-acetylmuramoyl-L-alanine amidase. J. Biol. Chem. 2003, 278, 49044–49052. [Google Scholar] [CrossRef] [Green Version]
  65. Herbold, D.R.; Glaser, L. Interaction of N-acetylmuramic acid L-alanine amidase with cell wall polymers. J. Biol. Chem. 1975, 250, 7231–7238. [Google Scholar] [CrossRef]
  66. Zhang, Y.; Van Der Fits, L.; Voerman, J.S.; Melief, M.J.; Laman, J.D.; Wang, M.; Wang, H.; Wang, M.; Li, X.; Walls, C.D. Identification of serum N-acetylmuramoyl-l-alanine amidase as liver peptidoglycan recognition protein 2. Biochim. Biophys. Acta BBA -Proteins Proteom. 2005, 1752, 34–46. [Google Scholar] [CrossRef] [PubMed]
  67. Pinheiro da Silva, F.; Cataldi, T.R.; de Lima, T.M.; Starzynski, P.N.; Barbeiro, H.V.; Labate, M.T.; CéMachado, M.C.; de Souza, H.P.; Labate, C.A. Proteomic profiling identifies N-acetylmuramoyl-l-alanine amidase as a novel biomarker of sepsis. Biomark. Med. 2016, 10, 1225–1229. [Google Scholar] [CrossRef] [PubMed]
  68. Wang, H.; Gupta, D.; Li, X.; Dziarski, R. Peptidoglycan recognition protein 2 (N-acetylmuramoyl-L-Ala amidase) is induced in keratinocytes by bacteria through the p38 kinase pathway. Infect. Immun. 2005, 73, 7216–7225. [Google Scholar] [CrossRef] [Green Version]
  69. Uehara, A.; Sugawara, Y.; Kurata, S.; Fujimoto, Y.; Fukase, K.; Kusumoto, S.; Satta, Y.; Sasano, T.; Sugawara, S.; Takada, H. Chemically synthesized pathogen-associated molecular patterns increase the expression of peptidoglycan recognition proteins via toll-like receptors, NOD1 and NOD2 in human oral epithelial cells. Cell. Microbiol. 2005, 7, 675–686. [Google Scholar] [CrossRef]
  70. Wettenhall, J.M.; Smyth, G.K. limmaGUI: A graphical user interface for linear modeling of microarray data. Bioinformatics 2004, 20, 3705–3706. [Google Scholar] [CrossRef] [Green Version]
  71. Njei, B.; McCarty, T.R.; Sharma, P.; Lange, A.; Najafian, N.; Ngu, J.N.; Ngomba, V.E.; Echouffo-Tcheugui, J.B. Bariatric surgery and hepatocellular carcinoma: A propensity score-matched analysis. Obes. Surg. 2018, 28, 3880–3889. [Google Scholar] [CrossRef]
  72. Lassailly, G.; Caiazzo, R.; Buob, D.; Pigeyre, M.; Verkindt, H.; Labreuche, J.; Raverdy, V.; Leteurtre, E.; Dharancy, S.; Louvet, A. Bariatric surgery reduces features of nonalcoholic steatohepatitis in morbidly obese patients. Gastroenterology 2015, 149, 379–388. [Google Scholar] [CrossRef] [Green Version]
  73. Simó, R.; Sáez-López, C.; Barbosa-Desongles, A.; Hernández, C.; Selva, D.M. Novel insights in SHBG regulation and clinical implications. Trends Endocrinol. Metab. 2015, 26, 376–383. [Google Scholar] [CrossRef] [PubMed]
  74. Kalme, T.; Seppälä, M.; Qiao, Q.; Koistinen, R.; Nissinen, A.; Harrela, M.; Loukovaara, M.; Leinonen, P.; Tuomilehto, J. Sex Hormone-Binding Globulin and Insulin-Like Growth Factor-Binding Protein-1 as Indicators of Metabolic Syndrome, Cardiovascular Risk, and Mortality in Elderly Men. J. Clin. Endocrinol. Metab. 2005, 90, 1550–1556. [Google Scholar] [CrossRef] [PubMed]
  75. Ding, E.L.; Song, Y.; Manson, J.E.; Hunter, D.J.; Lee, C.C.; Rifai, N.; Buring, J.E.; Gaziano, J.M.; Liu, S. Sex hormone–binding globulin and risk of type 2 diabetes in women and men. N. Engl. J. Med. 2009, 361, 1152–1163. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Sutton-Tyrrell, K.; Wildman, R.; Matthews, K.; Chae, C.; Lasley, B.; Brockwell, S.; Pasternak, R.; Lloyd-Jones, D.; Sowers, M.; Torrens, J. Sex-hormone-binding globulin and the free androgen index are related to cardiovascular risk factors in multiethnic premenopausal and perimenopausal women enrolled in the Study of Women Across the Nation (SWAN). Circulation 2005, 111, 1242–1249. [Google Scholar] [CrossRef] [Green Version]
  77. Jansen, A.; Berg, J.P.; Klungsøyr, O.; Müller, M.H.; Lyche, J.L.; Aaseth, J.O. The Influence of Persistent Organic Pollutants on Thyroidal, Reproductive and Adrenal Hormones after Bariatric Surgery. Obes. Surg. 2020, 30, 1368–1378. [Google Scholar] [CrossRef] [PubMed]
  78. Sarwer, D.B.; Spitzer, J.C.; Wadden, T.A.; Rosen, R.C.; Mitchell, J.E.; Lancaster, K.; Courcoulas, A.; Gourash, W.; Christian, N.J. Sexual functioning and sex hormones in men who underwent bariatric surgery. Surg. Obes. Relat. Dis. 2015, 11, 643–651. [Google Scholar] [CrossRef]
  79. Larsen, B.T.; Campbell, W.B.; Gutterman, D.D. Beyond vasodilatation: Non-vasomotor roles of epoxyeicosatrienoic acids in the cardiovascular system. Trends Pharmacol. Sci. 2007, 28, 32–38. [Google Scholar] [CrossRef]
  80. Yu, Z.; Xu, F.; Huse, L.M.; Morisseau, C.; Draper, A.J.; Newman, J.W.; Parker, C.; Graham, L.; Engler, M.M.; Hammock, B.D. Soluble epoxide hydrolase regulates hydrolysis of vasoactive epoxyeicosatrienoic acids. Circ. Res. 2000, 87, 992–998. [Google Scholar] [CrossRef] [Green Version]
  81. Newman, J.W.; Morisseau, C.; Hammock, B.D. Epoxide hydrolases: Their roles and interactions with lipid metabolism. Prog. Lipid Res. 2005, 44, 1–51. [Google Scholar] [CrossRef]
  82. Mustafa, S.; Sharma, V.; McNeill, J.H. Insulin resistance and endothelial dysfunction: Are epoxyeicosatrienoic acids the link? Exp. Clin. Cardiol. 2009, 14, e41–e50. [Google Scholar] [PubMed]
  83. Schmelzer, K.R.; Kubala, L.; Newman, J.W.; Kim, I.H.; Eiserich, J.P.; Hammock, B.D. Soluble epoxide hydrolase is a therapeutic target for acute inflammation. Proc. Natl. Acad. Sci. USA 2005, 102, 9772–9777. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Inceoglu, B.; Wagner, K.; Schebb, N.H.; Morisseau, C.; Jinks, S.L.; Ulu, A.; Hegedus, C.; Rose, T.; Brosnan, R.; Hammock, B.D. Analgesia mediated by soluble epoxide hydrolase inhibitors is dependent on cAMP. Proc. Natl. Acad. Sci. USA 2011, 108, 5093–5097. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  85. Luther, J.M.; Brown, N.J. Epoxyeicosatrienoic acids and glucose homeostasis in mice and men. Prostaglandins Other Lipid Mediat. 2016, 125, 2–7. [Google Scholar] [CrossRef] [Green Version]
  86. Habieb, M.S.; Dawood, A.A.; Emara, M.M.; Elhelbawy, M.G.; Elhelbawy, N.G. The Human Genetic Variants CYP2J2 rs2280275 and EPHX2 rs751141 and Risk of Diabetic Nephropathy in Egyptian Type 2 Diabetic Patients. Appl. Clin. Genet. 2020, 13, 165. [Google Scholar] [CrossRef] [PubMed]
  87. Khadir, A.; Kavalakatt, S.; Madhu, D.; Cherian, P.; Al-Mulla, F.; Abubaker, J.; Tiss, A. Soluble epoxide hydrolase 2 expression is elevated in obese humans and decreased by physical activity. Int. J. Mol. Sci. 2020, 21, 2056. [Google Scholar] [CrossRef] [Green Version]
  88. Lee, C.R.; North, K.E.; Bray, M.S.; Fornage, M.; Seubert, J.M.; Newman, J.W.; Hammock, B.D.; Couper, D.J.; Heiss, G.; Zeldin, D.C. Genetic variation in soluble epoxide hydrolase ( EPHX2 ) and risk of coronary heart disease: The Atherosclerosis Risk in Communities (ARIC) study. Hum. Mol. Genet. 2006, 15, 1640–1649. [Google Scholar] [CrossRef]
  89. Imig, J.D.; Hammock, B.D. Soluble epoxide hydrolase as a therapeutic target for cardiovascular diseases. Nat. Rev. Drug Discov. 2009, 8, 794–805. [Google Scholar] [CrossRef] [Green Version]
  90. Swann, D.A.; Slayter, H.S.; Silver, F.H. The molecular structure of lubricating glycoprotein-I, the boundary lubricant for articular cartilage. J. Biol. Chem. 1981, 256, 5921–5925. [Google Scholar] [CrossRef]
  91. Iozzo, R.V.; Schaefer, L. Proteoglycan form and function: A comprehensive nomenclature of proteoglycans. Matrix Biol. 2015, 42, 11–55. [Google Scholar] [CrossRef]
  92. Das, N.; Schmidt, T.A.; Krawetz, R.J.; Dufour, A. Proteoglycan 4: From Mere Lubricant to Regulator of Tissue Homeostasis and Inflammation: Does proteoglycan 4 have the ability to buffer the inflammatory response? Bioessays 2019, 41, 1800166. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  93. Geyer, P.E.; Wewer Albrechtsen, N.J.; Tyanova, S.; Grassl, N.; Iepsen, E.W.; Lundgren, J.; Madsbad, S.; Holst, J.J.; Torekov, S.S.; Mann, M. Proteomics reveals the effects of sustained weight loss on the human plasma proteome. Mol. Syst. Biol. 2016, 12, 901. [Google Scholar] [CrossRef] [PubMed]
  94. Oller Moreno, S.; Cominetti, O.; Núñez Galindo, A.; Irincheeva, I.; Corthesy, J.; Astrup, A.; Saris, W.H.; Hager, J.; Kussmann, M.; Dayon, L. The differential plasma proteome of obese and overweight individuals undergoing a nutritional weight loss and maintenance intervention. Proteom. Clin. Appl. 2018, 12, 1600150. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Principal component analysis of all identified proteins. Pre-surgery samples are coloured in blue, 6-month post-surgery samples in green and 12-month post-surgery samples in red.
Figure 1. Principal component analysis of all identified proteins. Pre-surgery samples are coloured in blue, 6-month post-surgery samples in green and 12-month post-surgery samples in red.
Jcm 10 03659 g001
Figure 2. Protein abundance changes 6 months post-surgery. (A) Volcano plot: green points are significant and have a fold change >1, red points are significant but have a fold change <1 and black points are not significant. Protein abundance changes 6 months post-surgery. (B) Hierarchical clustering and heatmap of all 25 proteins significantly different 6 months post-surgery relative to pre-surgery samples. Hierarchical clustering was performed using the Euclidean distance. Row scaling was performed on the log2 abundance for each protein by the subtraction of the mean from each feature and then dividing by the standard deviation. Protein abundance changes 6 months post-surgery (C) Principal component analysis with all 25 proteins significantly different 6 months post-surgery relative to pre-surgery samples.
Figure 2. Protein abundance changes 6 months post-surgery. (A) Volcano plot: green points are significant and have a fold change >1, red points are significant but have a fold change <1 and black points are not significant. Protein abundance changes 6 months post-surgery. (B) Hierarchical clustering and heatmap of all 25 proteins significantly different 6 months post-surgery relative to pre-surgery samples. Hierarchical clustering was performed using the Euclidean distance. Row scaling was performed on the log2 abundance for each protein by the subtraction of the mean from each feature and then dividing by the standard deviation. Protein abundance changes 6 months post-surgery (C) Principal component analysis with all 25 proteins significantly different 6 months post-surgery relative to pre-surgery samples.
Jcm 10 03659 g002aJcm 10 03659 g002b
Figure 3. Protein abundance changes 12 months post-surgery. (A) Volcano plot: green points are significant and have a fold change >1, red points are significant but have a fold change <1 and black points are not significant. Protein abundance changes 12 months post-surgery. (B) Hierarchical clustering and heatmap of all 39 proteins significantly different 12 months post-surgery relative to pre-surgery samples. Hierarchical clustering was performed using the Euclidean distance. Row scaling was performed on the log2 abundance for each protein by the subtraction of the mean from each feature and then dividing by the standard deviation. Protein abundance changes 12 months post-surgery. (C) Principal component analysis with all 39 proteins significantly different 12 months post-surgery relative to pre-surgery samples.
Figure 3. Protein abundance changes 12 months post-surgery. (A) Volcano plot: green points are significant and have a fold change >1, red points are significant but have a fold change <1 and black points are not significant. Protein abundance changes 12 months post-surgery. (B) Hierarchical clustering and heatmap of all 39 proteins significantly different 12 months post-surgery relative to pre-surgery samples. Hierarchical clustering was performed using the Euclidean distance. Row scaling was performed on the log2 abundance for each protein by the subtraction of the mean from each feature and then dividing by the standard deviation. Protein abundance changes 12 months post-surgery. (C) Principal component analysis with all 39 proteins significantly different 12 months post-surgery relative to pre-surgery samples.
Jcm 10 03659 g003aJcm 10 03659 g003b
Figure 4. (AH): Changes in protein expression between baseline, 6 months and 12 months for the 8 differentially expressed proteins at both 6 months and 12 months post-surgery.
Figure 4. (AH): Changes in protein expression between baseline, 6 months and 12 months for the 8 differentially expressed proteins at both 6 months and 12 months post-surgery.
Jcm 10 03659 g004
Table 1. Anthropometric and clinical measurements taken from patients before and after bariatric surgery.
Table 1. Anthropometric and clinical measurements taken from patients before and after bariatric surgery.
Baseline6 Months12 Monthsp-Value
Number of samples101010
Age (years)53 CI (46–60)
Sex2 Males, 8 Females
Anthropometric and Clinical *
Height (cm)Male: 177 Female: 163 cm
BMI (kg/m2)54.9 (10.3)41.8 (6.96) *38.6 (5.08) * p < 0.0001
Systolic BP135 (22.6)123 (16.7)116 (14.4) *p < 0.05
Diastolic BP79.2 (17.3)67.9 (14.1)67.9 (11) *p < 0.05
Biochemical
HbA1c (mmol/mol)56.0 (14.1)40.0 (5.4) *38.2 (6.25) * p < 0.0001
Fasting plasma glucose7.9 (3.5)5.8 (2.0)4.8 (0.70)p < 0.05
HOMA2-IR0.57 (0.52)0.489 (0.6)0.34 (0.24)p = 0.57
HOMA2-B (%)28.9 (18.2)45.5 (32.7)48.9 (16.7)p = 0.14
* Significant compares to baseline; significant compared to 6 months; data are presented as mean (standard deviation); CI: confidence interval.
Table 2. Proteins showing a significant difference in expression between baseline and six months post-surgery.
Table 2. Proteins showing a significant difference in expression between baseline and six months post-surgery.
Protein Names (Gene Symbol)Fold Change (Log2)p Value
Serotransferrin (TF)−1.068<0.0001
Beta-Ala-His dipeptidase (CNDP1) −0.9520.006
Complement C3 (C3)−0.5550.008
Proteoglycan 4 (PRG4)−0.8980.011
Cholinesterase (BCHE)−0.5810.012
Sex hormone-binding globulin (SHBG) 1.4850.013
Apolipoprotein M (APOM)0.6890.018
Protein disulfide-isomerase A3 (PDIA3)1.5270.021
Bifunctional epoxide hydrolase 2 (EPHX2) −0.4120.022
Apolipoprotein A-IV (APOA4)−0.7130.023
N-acetylmuramoyl-L-alanine amidase (PGLYRP)0.3910.024
Angiotensinogen (AGT) −1.0360.025
T-complex protein 1 subunit theta (CCT8) −1.2350.026
Heat shock 70 kDa protein 4 (HSPA4)−0.4970.026
Leucine-rich alpha-2-glycoprotein (LRG1)0.5320.032
Pigment epithelium-derived factor (SERPINF1)−0.3910.032
C4b-binding protein alpha chain (C4BPA) −0.3930.032
Exportin-1 (XPO1)−0.7460.033
Immunoglobulin heavy constant mu (IGHM)−0.9780.034
Alpha-1-acid glycoprotein 2 (ORM2) −0.7420.036
Sulfhydryl oxidase 1 (QSOX1) −0.7750.037
Apolipoprotein C-III (APOC3)−1.1710.039
Basement membrane-specific heparan sulfate proteoglycan core protein (HSPG2)−0.4350.041
Heparin cofactor 2 (SERPIND1)−0.4150.042
Aminopeptidase (ANPEP)−0.3950.048
Table 3. Proteins showing a significant difference in expression between baseline and 12 months post-surgery.
Table 3. Proteins showing a significant difference in expression between baseline and 12 months post-surgery.
Protein Names (Gene Symbol)Fold Change (log2)p Value
Serum amyloid p-component (SAP) −0.9360.001
Apolipoprotein A-I (APOA1)1.1040.001
Serotransferrin (TF)−0.7760.002
Inter-alpha-trypsin inhibitor heavy chain H3 (ITIH3)0.7990.002
Fibronectin (FN1) −1.7240.003
Sex hormone-binding globulin (SHBG) 1.9510.003
Neural cell adhesion molecule L1-like protein (L1CAM)0.6630.007
Haptoglobin (HP)−1.3830.007
Apolipoprotein A-IV (APOA4)−1.3810.008
Antithrombin-III (SERPINC1) 0.5410.008
Peptidase inhibitor 16 (PI16)0.7420.009
Galectin-3-binding protein (LGALS3BP) −0.5900.009
Laminin subunit gamma-1 (LAMC1)−1.5060.012
Proteoglycan 4 (PRG4) −0.7790.015
Hemoglobin subunit delta (HBD) −1.0880.020
N-acetylmuramoyl-L-alanine amidase (PGLYRP2)0.4380.021
Dermcidin (DCD)0.9050.021
Heat shock 70 kDa protein 4 (HSPA4)−0.3890.022
Actin-related protein 2/3 complex subunit 1B (ARPC1B)−0.4030.023
Bifunctional epoxide hydrolase 2 (EPHX2) −0.4730.025
Hyaluronan-binding protein 2 (HABP2) 0.3580.025
Galactokinase (GALK1)0.8990.027
Alpha-2-macroglobulin (A2M)0.8010.027
Leucine-rich alpha-2-glycoprotein (LRG1)0.5910.027
Peroxiredoxin-6 (PRDX6) 0.6270.028
Thyroxine-binding globulin (SERPINA7)0.3650.030
Mitogen-activated protein 3 kinase 5 (MAP3K5)0.6320.031
Alpha-1B-glycoprotein (A1BG)0.4620.033
Complement factor H (CFH)−0.4150.036
Protein AMBP (AMBP) 0.3520.037
L-lactate dehydrogenase A chain (LDHA)−1.5380.037
Endoplasmin (HSP90B1)−0.4170.038
Angiotensin-converting enzyme (ACE)−0.4050.039
Plasma protease C1 inhibitor (SERPING1)0.9390.039
Desmin (DES)−0.9520.043
Vitamin D-binding protein (GC)0.3810.046
Secreted phosphoprotein 24 (SPP2) 0.4750.048
T-complex protein 1 subunit gamma (CCT3)−0.7220.049
Haemoglobin subunit beta (HBB)−1.1280.049
Table 4. Proteins showing a significant difference in expression at both 6- and 12-months post-surgery. Log fold changes are shown between baseline and 12 months.
Table 4. Proteins showing a significant difference in expression at both 6- and 12-months post-surgery. Log fold changes are shown between baseline and 12 months.
Protein Names (Gene Symbol)Fold Change (log2)p Value
Serotransferrin (TF)−0.7760.002
Apolipoprotein A-IV (APOA4)−1.3810.008
Heat shock 70 kDa protein 4 (HSPA4)−0.3890.022
Leucine-rich alpha-2-glycoprotein (LRG1)0.5910.027
N-acetylmuramoyl-L-alanine amidase (PGLYRP2)0.4380.021
Sex hormone-binding globulin (SHBG) 1.9510.003
Bifunctional epoxide hydrolase 2 (EPHX2) −0.4730.025
Proteoglycan 4 (PRG4)−0.7790.015
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Iqbal, Z.; Fachim, H.A.; Gibson, J.M.; Baricevic-Jones, I.; Campbell, A.E.; Geary, B.; Donn, R.P.; Hamarashid, D.; Syed, A.; Whetton, A.D.; et al. Changes in the Proteome Profile of People Achieving Remission of Type 2 Diabetes after Bariatric Surgery. J. Clin. Med. 2021, 10, 3659. https://doi.org/10.3390/jcm10163659

AMA Style

Iqbal Z, Fachim HA, Gibson JM, Baricevic-Jones I, Campbell AE, Geary B, Donn RP, Hamarashid D, Syed A, Whetton AD, et al. Changes in the Proteome Profile of People Achieving Remission of Type 2 Diabetes after Bariatric Surgery. Journal of Clinical Medicine. 2021; 10(16):3659. https://doi.org/10.3390/jcm10163659

Chicago/Turabian Style

Iqbal, Zohaib, Helene A. Fachim, J. Martin Gibson, Ivona Baricevic-Jones, Amy E. Campbell, Bethany Geary, Rachelle P. Donn, Dashne Hamarashid, Akheel Syed, Anthony D. Whetton, and et al. 2021. "Changes in the Proteome Profile of People Achieving Remission of Type 2 Diabetes after Bariatric Surgery" Journal of Clinical Medicine 10, no. 16: 3659. https://doi.org/10.3390/jcm10163659

APA Style

Iqbal, Z., Fachim, H. A., Gibson, J. M., Baricevic-Jones, I., Campbell, A. E., Geary, B., Donn, R. P., Hamarashid, D., Syed, A., Whetton, A. D., Soran, H., & Heald, A. H. (2021). Changes in the Proteome Profile of People Achieving Remission of Type 2 Diabetes after Bariatric Surgery. Journal of Clinical Medicine, 10(16), 3659. https://doi.org/10.3390/jcm10163659

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop