1. Introduction
In developed countries, type 2 diabetes mellitus (T2DM) represents a major public health problem mainly by its relationship with different cardiovascular diseases (CVD) [
1]. Long-term hyperglycemia results in the synthesis of a number of molecules like advanced glycation end products (AGEs), advanced oxidation protein products (AOPPs), and low-density lipoprotein susceptibility to oxidation (oxLDL). Besides chronic hyperglycemia, other factors such as insulin resistance, dyslipidemia and states of inflammation and oxidation are related to vascular injury in diabetes through several underlying processes [
2]. However, only some diabetic patients develop cardiovascular disease while others do not develop these complications despite having the same common risk factors. The etiologic pathway linking impaired glucose tolerance and cardiovascular disease remains to be clarified. Many factors, including genetic components may be involved and not all are well established. Currently, multiple areas of research are open to explain this complex phenomenon [
3] and the precise role of the different disturbed metabolic pathways is not well established. In this context, the identification of new molecules that take part in the development of these vascular complications in T2DM may be of great importance for improving the outcomes of this population or to design new therapeutic targets. Thus, the underlying factors that may lead to the development of CVD require further in-depth research. The proteomic analysis is a hypothesis-free approach integrating genetic and epigenetic influences by examining the protein expression profiles and is not limited by the above knowledge.
Human serum is rich in potential biomarkers reflecting the pathophysiological state of the human body and related to the impaired metabolic pathways present in different disorders. In this way, it is able to facilitate the early detection of many diseases. However, finding biomarkers in serum is difficult due to the interference of the main proteins that provide a significant background covering low abundance proteins [
4].
To overcome this problem, serum depletion is essential for removing the major proteins and to improve the sensitivity of detection methods for these minor proteins.
The main objective of this cross-sectional study was to compare the serum proteome of the T2DM patients with and without CVD and non-diabetic controls to find protein differences in the protein profile expression of the study groups. Candidate proteins were identified by MALDI TOF/TOF or LC-MS/MS in order to identify essential proteins related to severe cardiovascular disease in patients with T2DM.
Here we report that proteins involved in inflammatory and redox state, such as retinol binding protein (RBP4) and glutathione peroxidase 3 (GPx-3) have increased levels in the serum of T2DM patients with severe CVD compared to those without CVD and non-diabetic controls. These differentially expressed proteins could be potential markers used for diagnosis or prognosis of CVD.
3. Discussion
Diabetes mellitus currently represents a serious health problem since affects a large proportion of the population [
5]. The patients with type 2 diabetes mellitus have increased risk for developing many complications mainly cardiovascular events depicting the major cause of mortality of this disease [
6,
7]. The probability of developing vascular complications in this population depends on a number of traditional factors and others that are not fully known. The knowledge of new factors involved in these disorders may facilitate the identification of the subjects in whom an intensive approach on cardiovascular risk factors should be established before irreversible damage occurs.
The present study was conducted to analyze the differences in the serum proteome profile of the T2DM patients with and without severe CVD in order to identify the proteins which could be related to higher cardiovascular risk in this population. Since the majority of patients with prolonged duration of T2DM usually show some sign of subclinical vascular disease, we select patients with established cardiovascular disease to identify the proteomic profile associated with more severe and advanced cardiovascular disease.
Our proteomic analysis by using 2D-DIGE and LC-MS/MS revealed changes in proteins related to the inflammatory and redox state. Over recent years, the inflammatory state is becoming more important as a factor involved in the cardiovascular disorders. We found increased RBP4 concentrations in the T2DM patients, being significantly higher in those with CVD compared to those without CVD and the control subjects. RBP4 is a protein, member of the lipocalin family, involved in diverse functions like sensory transduction and carrier of retinol. In plasma, the RBP4-retinol complex interacts with transthyretin, preventing its loss by filtration through the kidney glomeruli. RBP4 is generated by mature adipocytes [
8] and activated macrophages [
9]. The involvement of RBP4 in insulin resistance has been described in some studies [
8,
10,
11]. Clinical studies have shown that exercise, or other weight loss interventions, could improve the insulin sensitivity by reducing RBP4 serum levels [
12,
13]. Furthermore, a very recent study determines the mechanism by which RBP4 induces insulin resistance through activation of the antigen presenting cells (APC) of the adipose tissue. This activation induces higher levels of proinflammatory cytokines resulting in an inflammation of the adipose tissue that leads to insulin resistance [
14].
However, the association of RBP4 with cardiovascular events and their potential utility as a biomarker of CVD is not well known. In our study we found an increased expression of RBP4 in diabetic patients with cardiovascular disease in regard to the others groups. These findings are consistent with two preliminary studies which show a relationship between RBP4 levels and common cardiovascular risk factors in the elderly population [
15,
16]. Furthermore, Farjo
et al. [
17] have recently shown that RBP4 promotes the expression of proinflammatory molecules such as vascular cell adhesion molecule 1 (VCAM-1), intercellular adhesion molecule 1 (ICAM-1),
E-selectin and interleukin 6 (IL-6) through activation of NF-κB and NADPH oxidase in an
in vitro study. This effect is RBP4-dose dependent and induces the adhesion of leukocytes to the endothelium leading to endothelial inflammation. Our results are in agreement with these recent data suggesting that the elevation of RBP4 could be one of the main causes involved in the development and the progression of CVD associated to T2DM.
Another factor involved in the inflammatory state and vascular lesions is the oxidative stress due to the increase of reactive oxygen species (ROS) [
18]. To offset the toxic effects of free radicals produced during vascular lesions over cells, a cascade of redox reactions between nitric oxide (NO), a neurotransmitter with vasodilatory function, and ROS takes place. This process increases lipid peroxidation [
19] causing cell damage of the vascular system and aggravating CVD in these patients. Since oxidized lipoproteins constitute one of the most important factors involved in the development of atherosclerosis [
20], the endogenous defense systems are essential in the prevention of these complications. According to this, our results revealed an increase of GPx-3 in the serum of diabetes patients with CVD compared to those without CVD and controls. GPx-3 belongs to the family of glutathione peroxidases (GPx’s), a group of enzymes that reduce the oxidative stress decreasing vascular injury. GPx-3 is a homotetrameric protein present in high density lipoproteins (HDL) particles and secreted to the plasma that protects cells and enzymes catalyzing the reduction of hydrogen peroxide, lipid peroxides and organic hydroperoxides. It contributes to maintaining the vascular bioavailability of NO. We postulate that in patients with CVD an increase of GPx-3 levels takes place as a protective mechanism to reduce lipid peroxidation and to maintain the NO levels. There is some evidence supporting this potential protective effect. Regarding this, some studies link increased levels of GPx-3 with vascular complications [
21], and a deficiency or a reduction in GPx-3 are related to an increase in the predisposition to developing a thrombotic disorder [
22]. Thus, levels of GPx-3 could act as a predictor of cardiovascular risk.
Our results also show a significant decrease in TTR levels in the serum of the T2DM patients presenting CVD as well as a decrease in the chain A of C3 complement levels in the T2DM group compared to the other groups. Although several studies support these findings [
23,
24], our results about C3 and TTR are not conclusive since these proteins are depleted by the spin cartridges. Thus, results obtained by 2D-DIGE referents to C3 complement and TTR may not reflect the pathophysiological serum levels of patients and controls. So, future studies are needed to clarify the role of these proteins in the development of T2DM and CVD.
The strengths of this study lie in the novel evaluation of the differences in the serum proteome between T2DM patients with and without severe CVD and the identification of the proteins probably involved in the pathogenesis of CVD associated to T2DM and in the strict selection of study subjects. This allows the minimization of confounding factors resulting in a high reproducibility of the results in all replicates. However, our study has certain limitations. A case-control study does not allow us to establish a cause-effect relationship. Moreover, a biochemical validation of the protein candidates by ELISA in a larger T2DM study population will be helpful in order to confirm these preliminary results.
In conclusion, our study shows that proteins involved in the inflammatory state and oxidative stress are associated with cardiovascular complications in patients with type 2 diabetes. The usefulness of these proteins as biomarkers of cardiovascular risk or potential therapeutic targets of vascular disease will require further study.
4. Subjects and Methods
4.1. Study Population
Our cross-sectional study included 18 males divided according to three study groups: (i) T2DM patients with CVD; (ii) T2DM patients without CVD; (iii) Non-diabetic subjects as controls. To achieve 85.0% statistical power to detect differences among groups it is necessary to include six experimental units per group, resulting in a total sample size of 18 subjects in the experimental study.
Diagnosis of diabetes was according to the American Diabetes Association criteria (2005). From January 2006 to December 2007 we consecutively recruited patients who had been referred to our outpatient clinic from primary care centers for treatment of diabetes. The inclusion criteria for patients with CVD were coronary heart disease (previous myocardial infarction or coronary revascularization surgery), cerebrovascular disease (ischemic stroke), or ischemic peripheral arterial disease.
The control group consisted of non-diabetic subjects consecutively recruited from the general community in the same period of time and matched with diabetic patients for demographic characteristics.
All were Caucasians, ambulatory and had neither renal, hepatic, gastrointestinal nor thyroid diseases. All T2DM patients were receiving medications for diabetes, including metformin, sulfonylureas, insulin or a combination of these drugs.
Samples of venous blood were obtained from patients in the Endocrinology Unit of University Hospital San Cecilio in Granada. Samples of venous blood in vacutainer tubes containing no anticoagulant were taken in the morning after fasting overnight. Samples were incubated at room temperature for 30–45 min and centrifuged for 15 min at 7500 rpm. Then, the supernatant was carefully aspirated by a serologic pipette, aliquoted into 50 µL cryovials and stored at −80 °C until the examination.
The study was conducted with the approval of the ethical committee of the San Cecilio University Hospital and conformed to the relevant ethical guidelines for human and animal research (Project ID: PI 0514-2012. Research Ethics Committee of Granada Center (CEI-Granada) at 26 November 2012). Written informed consent was obtained from all subjects.
4.2. Depletion, Concentration and Cleaning of Serum Samples
A chromatographic removal of the 14 most abundant serum proteins (albumin, IgG, antitrypsin, IgA, transferrin, haptoglobin, fibrinogen, alpha 2-macroglobulin, alpha 1-acid glycoprotein, IgM, apolipoprotein AI, apolipoprotein AII, complement C3 and transthyretin) was performed in 50 µL of whole serum using Multiple Affinity Removal Spin Cartridges (Agilent, Santa Clara, CA, USA) according to the user manual. The depleted fractions were concentrated by a 5 kDa Mw cut off spin concentrators (Agilent) at 4000 rpm, 60 min at 10 °C. The ReadyPrep 2D cleanup kit (BioRad, Hercules, CA, USA) was used to reduce streaking, background staining, and other gel artifacts that may contaminate the 2D/IEF samples. The resultant pellet was dispersed into the adequate volume of 2D gel electrophoresis buffer (8 M Urea, 4% CHAPS) to solubilize protein samples. The samples were separated in aliquots of 11.5 µL for DIGE assay and 0.5 µL for protein determination, pH adjustment and 1D electrophoresis.
4.3. 2D-DIGE
4.3.1. CyDye Fluors Preparation and Sample Labeling
The CyDye DIGE fluors (Kit labeling minimal CyDye DIGE fluor 1 × 5 nmol (GE Healthcare, Freiburg, Germany)) were dissolved in dimethylformamide up to 400 pmol/mL of working solution.
Sodium hydroxide at 0.1 M diluted 1:10 was added to each sample to adjust pH at 8.5–9 for the labeling reaction. Randomization of samples across gels removes any bias from the experiments. Fifty µg of protein from each sample was minimally labeled on ice for 30 min with 1 µL of corresponding fluorescent dye working solution according to the Ettan DIGE System user manual (GE Healthcare). The internal standard (IS) was generated by combining equal amounts of each sample and labeled with Cy2 fluorescent dye. The reaction was stopped by addition and incubation of 10 mM lysine for 10 min on ice in the dark (1 µL Lys per 50 µg protein). To prepare the samples for each gel, two samples were mixed with 50 µg of IS according to the user manual of 2D-DIGE GE Healthcare (
Table 3). Twenty mM DTT and 0.5 mM PMSF were added to the mix. Finally, a total of nine 2D gels were run.
Table 3.
Sample arrangement on each gel and staining with each corresponding fluorophore. A: T2DM + CVD group; B: T2DM − CVD group; IS: Internal standard.
Table 3.
Sample arrangement on each gel and staining with each corresponding fluorophore. A: T2DM + CVD group; B: T2DM − CVD group; IS: Internal standard.
Gel Number | Cy2 | Cy3 | Cy5 |
---|
1 | IS | Control 1 | Sample B4 |
2 | IS | Sample A1 | Control 4 |
3 | IS | Sample B1 | Sample A4 |
4 | IS | Control 2 | Sample B5 |
5 | IS | Sample A2 | Control 5 |
6 | IS | Sample B2 | Sample A5 |
7 | IS | Control 3 | Sample B6 |
8 | IS | Sample A3 | Control 6 |
9 | IS | Sample B3 | Sample A6 |
4.3.2. 2D Gel Electrophoresis
IPG (0.5%) buffer was added to destreak solution (GE Healthcare) after 30 min to temper the solution at RT. The adequate volume was added to the sample and placed into the strip holders. The IPG strips (pH 3–10, 24 cm, GE Healthcare) were located over the sample, covered with mineral oil and focused on an IPGphor III (GE Healthcare) following these steps: 1 h at 500 V, 1 h at 1000 V, 3 h at 8000 V, 8000 V until 20,000 V/h and 500 V as maintenance step. The run was monitorized by IPGphor III software (GE Healthcare) ensuring that isoelectric focusing (IEF) process was successful. The strips were twice equilibrated for 15 min in the adequate buffer containing 75 mM Tris buffer (pH 8.8), 6 M Urea, 30% glycerol, 2% SDS and Bromophenol Blue. The first equilibration was with 1% DTT, and the second one with 2.5% Iodoacetamide in the equilibration buffer. The second dimension was carried out using 12% polyacrylamide gels at 13 mA per gel for 1 h in the first step and 17 mA per gel for 5.50 h in the second step. The proteins were visualized by Typhoon 9400 fluorescence scanner (GE Healthcare). The image analysis was performed by the specific software DeCyder 7.0 (GE Healthcare), that performs an automatic analysis and a co-detection of the fluorescent gels images by background subtraction, normalization, inter-gel matching and quantization. By using the statistical analysis by Decyder 7.0 (GE Healthcare) the sum of the pixel values within a spot minus background for each experimental sample is compared directly to the internal standard. Thereby, the protein abundance for each spot is expressed as a normalized ratio relative to the internal standard. This analysis compares the average ratio and variation within each group to the average ratio and variation in the other groups to see if any change between the groups is significant. Thus, it is possible to compare the protein abundance for a protein of interest in different samples. The total proteins were detected by Sypro staining (BioRad).
4.3.3. Identification of the Candidate Biomarker by MALDI-TOF/TOF
The analysis by Decyder 7.0 performs the matching of the gels and creates a pick list after the scan. Selected protein spots are automatically picked from 2D representative gel by Ettan Spot Picker (GE Healthcare) using the pick list generated. The gel plugs are transferred into microplates for further digestion prior to identification.
The gel pieces were reduced with DTT for 60 min at room temperature and cysteines were carbamidomethylated with iodoacetamide for 30 min. The proteins were digested with trypsin (Gold, MS Grade, Promega, Madison, WI, USA) for 8 h at 37 °C. The resulting peptide mixture was spotted on a MALDI plate with CHCA and analyzed using the 4800 MALDI TOF/TOF mass spectrometer (AB-Sciex, Madrid, Spain). The MALDI-TOF/TOF spectra were interpreted by database search MASCOT (Matrix Science, Boston, MA, USA) with a significance threshold of the MOWSE score of p < 0.05. The database used for identification was NCBI restricted to mammalian proteins with the following parameters: peptide mass tolerance 100 ppm, fragment tolerance 0.2 Da, enzyme set as trypsin and allowance up to one missed cleavage, variable modification of methionine oxidation (+16 Da), fixed modification of cysteine carbamidomethylation (+57 Da). Only the proteins with the better identification parameters and the best score during the search in databases have been included in the list of results. Other identified proteins with lower score were identified (Data not shown).
4.3.4. LC-MS/MS Analysis
The tryptic extracts were reanalyzed by high resolution LC-MS/MS in data-dependent mode. The MS system used was an Orbitrap XL (ThermoFisher, Waltham, MA, USA) equipped with a microESI ion source (Proxeon). The tryptic extracts were diluted to 20 µL with 5% methanol and 1% formic acid. Then, extracts were loaded into a chromatographic system consisting of a C18 preconcentration cartridge (Agilent, Santa Clara, CA, USA) connected to a 15 cm long, 100 µm i.d. C18 column (Nikkyo Technos Co., Tokyo, Japan). The separation was done at 0.4 µL/min in a 30-min acetonitrile gradient from 3% to 40% (solvent A: 0.1% formic acid, solvent B: acetonitrile 0.1% formic acid). The HPLC system was composed of a 1200 capillary nano pump (Agilent), a binary pump, a thermostated micro injector and a micro switch valve. The Orbitrap XL was operated in the positive ion mode with a spray voltage of 2 kV. The scan range for full scans was m/z 400–2000. The spectrometric analysis was performed in a data-dependent mode, acquiring a full scan followed by 10 MS/MS scans of the 10 most intense signals detected in the MS scan. An exclusion time of 30 s was included to avoid repetitive MS/MS analysis of the dominant MS signals. LC-MS/MS spectra were searched using SEQUEST (Proteome Discoverer v1.3, ThermoFisher) with the following parameters: peptide mass tolerance 10 ppm, fragment tolerance 0.8 Da, enzyme set as trypsin and allowance up to two missed cleavages, dynamic modification of methionine oxidation (+16 Da), fixed modification of cysteine carbamidomethylation (+57 Da).
The database used for searching included human proteins (Uniprot taxonomy_9606). The peptide identifications were filtered for 0.1% FDR and only proteins identified with two or more peptides were considered. Only the proteins with the better identification parameters and the best score during the search in databases have been included in the list of results. Other proteins with lower score were identified (Data not shown).
4.3.5. Western Blot Analysis
A biochemical validation of protein candidates was performed by western blot analysis using standard techniques. In order to get a better approximation to the physiologic state of the individuals we used whole serum samples. We performed the division of 15 samples from each study group into 3 groups of 5 individuals each one, giving a total of 9 groups (3 groups of Control, 3 groups of T2DM + CVD and 3 groups of T2DM − CVD). The 9 groups were compared in three independent Western blots. Equal volumes of the pooled samples were separated on a 1D SDS-PAGE using 15% polyacrylamide gels during 60 min at 150 V. The separated proteins were transferred to PVDF membrane on Trans-Blot Turb Transfer System (BioRad), following the transfer protocol for 1.5 mm gels. The membranes were blocked in phosphate-buffered saline (PBS) supplemented with 5% (w/v) non-fat dry milk (BioRad) for 2 h at RT. Afterwards the membranes were washed with PBS 0.05% Tween 20 (PBST) during 10 min and subsequently they were incubated overnight with polyclonal human anti-goat primary antibodies: Anti-serum retinol binding protein (RBP4) ((SAB2500875), Anti-Transthyretin (TTR) (SAB25001931), Anti-glutathione peroxidase 3 (GPx-3) (SAB2500485) Sigma-Aldrich, Madrid, Spain); Anti-C3 complement ((A304) Quidel, San Diego, CA, USA), at 4 °C in PBS containing 0.1% albumin (PBSA). After three washes with PBST, the membranes were incubated with corresponding secondary horseradish peroxidase-conjugated antibody (A5420, Sigma-Aldrich) diluted 1:20,000 in PBSA for 1 h at RT. The proteins were detected by chemiluminiscence using Clarity Western ECL Substrate (BioRad) and the ChemiDoc gel imaging system (BioRad). The bands were quantified in three independent experiments by Quantity One software (BioRad), normalized to value 1 corresponding to the control and expressed as intensity difference relative to the control group (ratio disease group/control group).
4.4. Statistical Analysis
The data for continuous variables are presented as mean ± standard deviation (SD). The data for categorical variables are presented as numbers and/or percentages. Kolmogorov-Smirnov test was used to test the normal distribution of continuous variables. Comparisons of categorical variables among groups were performed using the Chi-square test. Comparison of continuous variables between case-control was performed using the unpaired Student’s t test, or the Mann-Whitney test. The statistical significance was set at p < 0.05 (two-tailed) and performed with the specific software, SPSS 18.0 (Chicago, IL, USA).
For statistical analysis by Decyder 7.0, the individual gels were processed in Differential In-gel Analysis module (DIA) to spot detection and first in-gel analysis for the images from each single gel removing background and artefacts and normalizing the gel images. The Biological Variation Analysis module (BVA) was used to analyze DIA spot maps of all gels in one DIGE experiment, to merge gel spot maps, to normalize intensities and to obtain statistical differences between groups. The analysis of variance (ANOVA) was applied to matched spots. The data was filtered to retain spots with ANOVA p-values of 0.05 or less.
To get a power of 85% to find differences in testing the null hypothesis “The proportions of the three groups are equal using a Chi-square test for three independent samples, considering that the significance level is 5% and assuming that the expected proportions of the three groups are specified by the researcher (85%, 10%, 10%)”, it is necessary to include 18 experimental units distributed equally among the 3 groups (6 units per group).