If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
In-vitro Identification of Distinctive Metabolic Signatures of Intact Varicose Vein Tissue via Magic Angle Spinning Nuclear Magnetic Resonance Spectroscopy
Nuclear magnetic resonance (NMR) spectroscopy is an established tool for metabolic profiling of tissues or biofluids with utility in identifying disease biomarkers and changes in enzymatic or gene expression. This pilot study aims to compare the metabolic profiles of intact varicose and non-varicose vein tissue via magic angle spinning (MAS) NMR spectroscopy with a view to promoting the understanding of the pathogenesis of varicose vein formation.
Methods
Varicose vein tissue (n = 8) was collected from patients undergoing varicose veins surgery. Control non-varicose great saphenous vein samples were collected from patients undergoing lower limb amputation (n = 3) and peripheral arterial bypass surgery (n = 5). Intact tissue samples (average weight 10.33 ± 0.8 mg) from each vein segment were analysed using 1D MAS 1H NMR (600 MHz) spectroscopy. For selected vein samples, two-dimensional (2D) NMR experiments were performed. Differences between spectra from varicose and non-varicose tissues were elucidated using a variety of multivariate statistical analyses.
Results
The metabolic profiles of varicose veins samples were clearly differentiated from non-varicose veins samples. Lipid metabolites were present at a higher concentration in the non-varicose veins group whilst creatine, lactate and myo-inositol metabolites were more characteristic of the varicose veins group.
Conclusion
We demonstrate differential metabolic profiles between varicose veins and non-varicose veins. Elucidating the metabolic signature underlying varicose veins can further improve our understanding of the biological mechanisms of disease initiation, progression, and aid in identifying putative therapeutic targets.
Varicose veins affect one-third of adults in western population. Complications of the disease such as leg oedema and ulceration affect the patients' health and lifestyle. Pathogenesis of primary varicose veins is unclear. This novel study has elucidated the metabolic profile of intact varicose vein as compared to non-varicose vein. Understanding the significance of differential metabolic profile in future may untangle the biological mechanism of disease initiation and progression.
Background
Varicose veins affect one-third of adults in the Western world.
Marked morphological differences clearly distinguish varicose veins from non-varicose veins. Varicose veins also exhibit expression of various transcripts and proteins linked to specific genes regulating cell growth and apoptosis, extracellular matrix metabolism and inflammatory processes.
However, transcriptomic (study of individuals' gene expression as RNA) and proteomic (study of translation of RNA into proteins) analyses have not yet allowed elucidation of the cellular and molecular changes described in the wall of varicose veins, or the mechanisms responsible for disease progression.
Metabonomics is defined as “the quantitative measurement of the dynamic multi-parametric metabolic response of living systems to pathophysiological stimuli or genetic modification.”
‘Metabonomics’: Understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological nmr spectroscopic data.
Metabolites are intermediate or end-products of cellular metabolism and hold key information regarding ongoing biological processes. When exposed to radiofrequency pulses in an external magnetic field, magnetically active atomic nuclei such as protons in a biological sample show slight differences in their transition frequency depending upon their immediate electronic and chemical environment. These differences are the basis of the spectral fingerprint containing discrete signals reflecting multiple metabolites and their concentrations.
Various stimuli including disease, diet and drugs affect cellular metabolism and hence the metabolite pools in the extracellular tissues and fluids are affected as a result of this altered metabolism.
High-resolution magic angle spinning (MAS)-nuclear magnetic resonance (NMR) spectroscopy is a non-destructive approach to analyse intact tissue and is found to be valuable in providing information regarding the metabolites' concentration and compartmentalisation within the tissues.
Several studies have demonstrated toxicological and diagnostic applications of high-resolution MAS-NMR spectroscopy of intact tissues including colonic, renal and hepatic tissues.
Metabolic profiling of human colorectal cancer using high-resolution magic angle spinning nuclear magnetic resonance (hr-mas nmr) spectroscopy and gas chromatography mass spectrometry (gc/ms).
The biochemical composition of tissue spectra collected by MAS-NMR is then analysed via mathematical modelling using chemometric software to extract and identify robust information on cellular biological events. This can lead to the elucidation of cellular pathways, provide information on gene or enzyme functions and can aid in the determination of the influence, efficacy and toxicity of targeting therapeutics.
‘Metabonomics’: Understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological nmr spectroscopic data.
Identifying new cellular pathways using metabonomic approaches may address the undisclosed biological events, which have not been picked up by other systems biology approaches, including transcriptomics and proteomics. This preliminary study aims to assess the viability of using a metabonomic approach to investigate the biochemical composition of vein wall tissue and to explore the pathogenesis of varicose veins disease, which may help to improve our understanding of its aetiology and its management.
Methods
Research ethics committee approval was obtained from the Riverside Research Ethics Committee, London, UK. Varicose veins were collected from patients undergoing surgery for lower extremity varicose veins (n = 8). The presence of truncal reflux in the great saphenous vein (GSV) was confirmed and vein marking was undertaken, both using duplex ultrasound prior to surgery. Retrieved GSV varicose samples were also inspected carefully to select for analysis only dilated and tortuous vein segments. Three control non-varicose GSV samples were collected from patients undergoing lower limb amputation and five non-varicose GSV samples from individuals having peripheral arterial bypass surgery. Each vein sample was collected fresh at the time of surgery, snap frozen in liquid nitrogen and stored at −80 °C. The entire circumference of vein tissue was cut using a sterile scalpel and forceps; any peri-vascular fat was gently removed. To assess tissue homogeneity within an individual sample and to establish whether the metabolic profile is relatively uniform across a length of vein tissue, one vein segment from each group (varicose and non-varicose veins control) was divided into three adjacent parts and analysed independently. In total, 10 samples from eight varicose veins patients and 10 samples from eight non-varicose veins control patients were analysed. Due to ethical limitations upon the study design, patients from whom varicose veins tissue samples were obtained were younger and had fewer co-morbidities as compared with those from whom non-varicose veins were collected (Table 1).
Table 1Patient demographic and clinical features of varicose veins and non-varicose veins groups.
Metabolic profiling of human colorectal cancer using high-resolution magic angle spinning nuclear magnetic resonance (hr-mas nmr) spectroscopy and gas chromatography mass spectrometry (gc/ms).
and spun at 5 kHz in high-resolution MAS 600 MHz 1H-NMR spectroscope (Bruker BioSpin, Rheinstetten, Germany).
NMR spectroscopy of vein tissue
Standard one-dimensional spectra were used with suppression of the water signal to adjust the dynamic range to the metabolites of interest. The Carr–Purcell–Meiboom–Gill (CPMG) spin-echo pulse sequence was also applied in a second measurement to reduce effects of broader resonances from high-molecular-weight compounds, therefore enabling sharp resonances from low-molecular-weight compounds to be easily identified. A total spin–spin relaxation time of 120 ms was chosen (n = 300, τ = 200 μs). For each sample, 256 scans were acquired into 64 K data points. Spectral width of δ 20 ppm was used with relaxation delay of 2 s and acquisition time of 2.72 s. Each spectrum was phased, calibrated using the chemical shift of the alanine methyl group at δ1.48 or the β-glucose (1H) doublet at δ 4.644, and baseline was corrected using TOPSPIN 2.0a software (Bruker BioSpin, Rheinstetten, Germany).
Chemometric analysis of NMR spectra
Spectra were imported into MATLAB R2009b (Mathworks™, 2009) using in-house software. The region containing the water resonances (from δ4.68 to δ5.24) was removed from each spectrum and all spectra were subsequently aligned and normalised using probabilistic quotient normalisation.
Data were transferred to SIMCA-P+ 11.5 statistical software (UMETRICS™, Sweden). Principal components analysis (PCA) and orthogonal partial least square regression (OPLS) analysis were employed using unit variance (UV) scaled data to identify the presence of inherent similarities in the data and investigate whether the varicose veins group systematically differed in biochemical composition from the control group. OPLS is a supervised linear regression approach which can model separately the co-varying variation from structured noise (orthogonal variation), while simultaneously maximising the X and Y covariance, that is, the differences between two groups.
An in-house mathematical modelling algorithm in MATLAB (Mathworks, 2009) was employed to identify the strength of correlation of the influential metabolite peaks with a particular class. The colour-coded orthogonal 2 partial least square (O2PLS) regression coefficient plot identifies the strength of correlation of a metabolite to a particular tissue type, therefore highlighting significant differences between varicose and non-varicose veins. We performed receiver operating characteristic (ROC) curve analysis to examine the diagnostic accuracy of partial least-squares driven class discrimination. An area between 0.9 and 1 under the ROC curve represents a strong diagnostic accuracy of the model.
Individual metabolite identities were assigned by using data from other studies examining human tissue metabolic profiles via MAS NMR and using the human metabolome database.
Metabolic profiling of human colorectal cancer using high-resolution magic angle spinning nuclear magnetic resonance (hr-mas nmr) spectroscopy and gas chromatography mass spectrometry (gc/ms).
Heteronuclear 1h-31p statistical total correlation nmr spectroscopy of intact liver for metabolic biomarker assignment: Application to galactosamine-induced hepatotoxicity.
The assignments were confirmed by two-dimensional (2D) NMR experiments on selected varicose veins and non-varicose veins samples, and through the use of the Bruker Amix/SBase database. 2D NMR disperses the data into two spectra showing mainly correlation between adjacent protons (1H–1H Correlation Spectroscopy or COSY) or proton and carbon atoms 1H–13C heteronuclear single quantum coherence (HSQC) spectroscopy, facilitating elucidation of the structure of each metabolite.
Heteronuclear 1h-31p statistical total correlation nmr spectroscopy of intact liver for metabolic biomarker assignment: Application to galactosamine-induced hepatotoxicity.
J-resolved spectra which display the coupling patterns and J-coupling values in a second orthogonal dimension were also used to identify metabolites and interpret overlapping signal patterns.
1H NMR spectra from all samples with relevant varicose (n = 10) and non-varicose veins (n = 10) spectra are shown in Fig. 1. Assigned spectra from each group are shown in Figs. 2 and 3. Spectral similarities within each group and differences between the two groups of veins were evident. The large peaks of triglycerides (unsaturated and saturated) were found in abundance in the non-varicose veins group (Fig. 2). Creatine, myo-inositol, lactate and glutamate were found in relatively higher levels in varicose veins tissues (Fig. 3). These findings were supported by the results of 2D NMR on selected vein tissue.
Figure 1Spectra from 10 non-varicose vein samples and 10 varicose vein samples. Similarities of spectral peaks within the varicose veins group representing homogeneity within the group. In comparison, 10 spectra from non-varicose vein samples are different to the varicose veins spectra. Spectra coloured in blue in each group represent the samples taken from the single vein segment of a patient of that group. There is clearly more intra-segmental versus inter-sample homogeneity.
Analysis of the PCA scores plot of CPMG spectra of varicose and non-varicose veins (Fig. 4a) showed that the variance in the first two components was skewed by two outliers. These outliers were two samples out of a series of three biological replicates from the same vein. Outlier samples contained substantially higher concentrations of triglycerides compared to other samples. However, after exclusion of these two outlying samples, variation was still relatively higher in the non-varicose veins group compared to varicose veins samples (Fig. 4b).
Figure 4Principal components analysis (PCA) scores plot of varicose vein and non-varicose vein CPMG spectra, showing variance between samples and the direction of the largest variance in the data. Each ▪ or represents a vein sample. Percentage of variation in the NMR data explained by first and second principal components were 79% and 8%, respectively, in the presence of all replicates (a), and 53% and 19%, respectively, when two out of three replicates were removed from each group (b). Varicose vein samples are grouped together in the unsupervised PCA analysis. Intra-segmental versus inter-sample homogeneity is also noticeable in the non-varicose veins group.
There was evidence for clustering of a subset of the samples from the varicose veins group and so OPLS analysis was used to identify systematic variation relating to disease class. To avoid potential bias in OPLS discriminatory analyses, only one out of three replicates from the same vein segment in each group was included for the OPLS discriminatory analysis. Although the number of samples in this pilot study was small, the OPLS scores plot (Fig. 5) clearly highlights the separation of the two groups, with tighter clustering of the varicose veins samples compared to non-varicose veins samples.
Figure 5Orthogonal partial least square (OPLS) analysis separating the two classes (varicose and non-varicose veins) based on their metabolic profiles detected using MAS-1H NMR. Probabilistic quotient normalisation and unit variance settings were applied for data analysis. R2Ycum was 0.788 and Q2cum was 0.406 for orthogonal component. Q2 indicates the fraction of variation predicted by the model for cross validation. One non-varicose vein sample is closer to the varicose veins group (spectrum 5 in Fig. 1).
Examination of the O2PLS coefficients plot showed that creatine, myo-inositol and lactate were highly correlated with the varicose veins group (Fig. 6). The ROC curve analysis confirmed that the diagnostic accuracy of our model was excellent with an area under the curve score of 0.92 (Fig. 7).
Figure 6O2PLS co-efficient plot showing 1H NMR CPMG spectral data of veins using MAS technology. Spectra have been processed using probabilistic quotient normalisation and unit variance scaling. Red peaks represent highly correlated metabolites to the respective group and are responsible for the distinction between the two groups, yellow are less variable metabolites, and blue the least correlated. Highly correlated metabolites to the varicose veins group are creatine, lactate and myo-inositol. Statistical total correlated spectroscopy (STOCSY) of metabolic features at δ3.26 confirms the correlation of myo-inositol peaks.
Figure 7Receiver operating characteristic (ROC) curve diagnostic analysis for predictive O2PLS model. Area under the curve score of 0.9218 shows excellent predictive accuracy of our model.
We observed that the metabolic profiles of varicose veins vary significantly from non-varicose veins. The main differentiating biomarkers are creatine, myo-inositol and lactate, which were higher in varicose veins tissue than non-varicose veins tissue. On the other hand, triglycerides were found in abundance in non-varicose veins.
Creatine is a nitrogenous organic acid originating from the diet, and also synthesised in the liver. It is stored in skeletal muscle where it acts as a cellular energy source.
The effects of creatine ethyl ester supplementation combined with heavy resistance training on body composition, muscle performance, and serum and muscle creatine levels.
Intracellular compartmentation, structure and function of creatine kinase isoenzymes in tissues with high and fluctuating energy demands: The ‘phosphocreatine circuit’ for cellular energy homeostasis.
Phosphocreatine is found in a concentration ranging between 20 and 35 mM in skeletal muscle cells and 5–10 mM in other excitable tissue including brain and smooth muscle.
Intracellular compartmentation, structure and function of creatine kinase isoenzymes in tissues with high and fluctuating energy demands: The ‘phosphocreatine circuit’ for cellular energy homeostasis.
In our study, varicose veins samples were observed to be thick walled on inspection. Hypertrophied and hyperplastic smooth muscle mass in the varicose veins wall is therefore responsible for the creatine abundance observed in this group. However, varicose veins can also have smooth muscle atrophy and be thin walled; therefore, inclusion of a large sample including both thin and thick-walled varicose veins, with metabolite extraction for NMR and mass spectrometry analysis, is recommended for future studies.
Myo-inositol, in addition to regulating intracellular osmolality in the brain,
forms a structural component of the phosphatidylinositol and phosphatidylinositides, which function as second messengers. Phosphatidylinositides are involved in lipid signalling and cell signalling via intracellular pathways containing phosphatidylinositol kinases,
Decreased production of collagen type iii in cultured smooth muscle cells from varicose vein patients is due to a degradation by mmps: Possible implication of mmp-3.
NFκB signalling plays a role in the regulation of adhesion molecules, cytokines, chemokines, growth factors and matrix metalloproteinases (MMPs) in vascular cells.
Varicose veins have been found to be associated with upregulation of growth factors, including vascular endothelial growth factor, transforming growth factor-β1 and acidic fibroblast growth factors.
which occurs in a number of states including hypoxia. Importantly, hypoxia and increased activation of hypoxia-inducible pathway have been recently implicated in the pathogenesis of varicose veins.
Increased activity of the anaerobic fraction of LDH has also been reported in varicose veins suggesting increased anaerobic glycolysis in this context.
Association of the metabolites lactate and myo-inositol with such key signalling pathways proposes these pathways for future research.
To our knowledge, this is the first study to characterise the metabolic profile of varicose veins by NMR spectroscopy. From this study, it is apparent that metabonomic techniques are suitable for characterising vein pathology. A limitation of our study is that the control group did not accurately match with the varicose veins group in terms of demographic features, past medical history and medications. Factors such as age, medications and past medical history may influence the metabolic profile of vein tissue. However, to our knowledge, there has been no demonstration of the effect of age on metabolic profile of vein tissue so far. Statins and most of the other cardiovascular medications have a short half-life (less than 24 h)
and neither parent statins nor their metabolites were detected in the NMR profiles. Although this does not evade the possibility that an endogenous metabolic response to chronic use of statins might be detected, in light of the nature of the metabolic changes it seems unlikely that statins would be responsible for any of the observed differences in metabolic profile between the varicose vein and control groups. Each vein segment was washed with deuterium oxide before the experiment. It is challenging to obtain non-varicose vein control tissue from age-matched (often young and otherwise healthy) individuals. Use of inferior epigastric vein retrieved at inguinal hernia surgery may facilitate patient matching. However, this will not reflect the relevant haemodynamic stresses to which a GSV is exposed. We further plan to examine control vein tissue from various anatomical regions, including facial vein at the time of carotid endarterectomy, arm vein during haemodialysis fistula formation, inferior epigastric veins from inguinal hernia repairs and inferior mesenteric vein from bowel resections, to document any metabolic differences between the veins from different anatomical regions and to assign a metabolic profile of non-varicose vein.
This study uses 1H MAS NMR spectroscopy to rapidly obtain a picture of the metabolic profile of vein tissue. Furthermore, metabolites can also be extracted from tissue and analysed by NMR. Correlation of metabolic profiles with histopathological analysis of tissue samples may further advance our understanding of the pathophysiology of varicose veins. These approaches will form part of our ongoing investigative plan, to be applied in larger sample sets, for establishing differential metabolic signature of varicose veins as compared to non-varicose veins.
Conclusion
This preliminary study demonstrates differential metabolic expression in varicose as compared with non-varicose veins. Elucidating the metabolic signature underlying varicose vein disease has the potential to further our understanding of biological mechanisms of disease initiation and progression, identify putative biomarkers and targets for prevention and treatment and merits further explanation.
Acknowledgements
The authors would like to acknowledge financial support from the Waters Corporation, Imperial College Healthcare Trust Fund, Graham-Dixon Charitable Trust, European Venous Forum, and Royal Society of Chemistry, UK. We acknowledge the support of Olaf Beckonert and Kirill Veselkov for their help in chemometric analysis and performing ROC curve analysis.
Ethical Approval
Ethical approval was obtained from the Riverside Research Ethics Committee, RREC 3092.
Funding
Acknowledgement section contains the names of research funding organizations who supported this project.
‘Metabonomics’: Understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological nmr spectroscopic data.
Metabolic profiling of human colorectal cancer using high-resolution magic angle spinning nuclear magnetic resonance (hr-mas nmr) spectroscopy and gas chromatography mass spectrometry (gc/ms).
Heteronuclear 1h-31p statistical total correlation nmr spectroscopy of intact liver for metabolic biomarker assignment: Application to galactosamine-induced hepatotoxicity.
The effects of creatine ethyl ester supplementation combined with heavy resistance training on body composition, muscle performance, and serum and muscle creatine levels.
Intracellular compartmentation, structure and function of creatine kinase isoenzymes in tissues with high and fluctuating energy demands: The ‘phosphocreatine circuit’ for cellular energy homeostasis.
Decreased production of collagen type iii in cultured smooth muscle cells from varicose vein patients is due to a degradation by mmps: Possible implication of mmp-3.
In the study of varicose vein pathophysiology, there has been much attention focussing on the venous epidemiology, symptoms and risk factors,1–3 identification of histologic abnormalities in the venous wall,4 the structural changes that occur in collagen, elastin, laminin, fibronectin and tenascin,5 the differential expression of smooth muscle cell collagen I versus collagen III expression and the post-transcriptional modification of collagen by matrix metalloproteinases (MMPs) that take place affecting the production of collagen III,6,7 the presence of MMPs in varicose veins8–10 and how MMPs affect venous dilation and the mechanisms involving the endothelium and smooth muscle.
To submit a comment for a journal article, please use the space above and note the following:
We will review submitted comments as soon as possible, striving for within two business days.
This forum is intended for constructive dialogue. Comments that are commercial or promotional in nature, pertain to specific medical cases, are not relevant to the article for which they have been submitted, or are otherwise inappropriate will not be posted.
We require that commenters identify themselves with names and affiliations.