Next Article in Journal
Cooperative Cellular Uptake and Activity of Octaarginine Antisense Peptide Nucleic Acid (PNA) Conjugates
Next Article in Special Issue
Medical Therapy of Patients Contaminated with Radioactive Cesium or Iodine
Previous Article in Journal
The Involvement of the Banana F-Box Protein MaEBF1 in Regulating Chilling-Inhibited Starch Degradation through Interaction with a MaNAC67-Like Protein
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Biomolecules
Volume 9
Issue 10
10.3390/biom9100553
biomolecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Significant Changes in Metabolic Profiles after Intervention with Selenium and Coenzyme Q10 in an Elderly Population

by
Urban Alehagen
1,*,
Peter Johansson
2,3,4,
Jan Aaseth
5,
Jan Alexander
6,
Izabella Surowiec
7,
Katrin Lundstedt-Enkel
7 and
Torbjörn Lundstedt
7
1
Division of Cardiovascular Medicine, Department of Medical and Health Sciences, Linköping University, SE-581 85 Linköping, Sweden
2
Department of Social and Welfare Studies, Linköping University, SE-581 83 Linköping, Sweden
3
Department of Internal Medicine, Linköping University, SE-581 83 Linköping, Sweden
4
Department of Medical and Health Sciences, Linköping University, SE-581 83 Linköping, Sweden
5
Research Department, Innlandet Hospital Trust, N-2380 Brumunddal, Norway
6
Norwegian Institute of Public Health, N-0403 Oslo, Norway
7
AcureOmics AB, Tvistevägen 48, SE-907 36 Umeå, Sweden
*
Author to whom correspondence should be addressed.
Biomolecules 2019, 9(10), 553; https://doi.org/10.3390/biom9100553
Submission received: 27 August 2019 / Revised: 25 September 2019 / Accepted: 26 September 2019 / Published: 30 September 2019
(This article belongs to the Special Issue Toxic and Essential Metals in Human Health and Disease)
Browse Figures
Versions Notes

Abstract

:
Selenium and coenzyme Q10 (SeQ10) are important for normal cellular function. Low selenium intake leads to increased cardiovascular mortality. Intervention with these substances with healthy elderly persons over a period of four years in a double-blind, randomised placebo-controlled prospective study showed reduced cardiovascular mortality, increased cardiac function, and a lower level of NT-proBNP. Therefore, we wanted to evaluate changes in biochemical pathways as a result of the intervention with SeQ10 using metabolic profiling. From a population of 443 healthy elderly individuals that were given 200 µg selenium and 200 mg coenzyme Q10, or placebo daily for four years, we selected nine males on active intervention and nine males on placebo for metabolic profiling in the main study. To confirm the results, two validation studies (study 1 n = 60 males, study 2 n = 37 males) were conducted. Principal component analyses were used on clinical and demographic data to select representative sets of samples for analysis and to divide the samples into batches for analysis. Gas chromatography time-of-flight mass spectrometry-based metabolomics was applied. The metabolite data were evaluated using univariate and multivariate approaches, mainly T-tests and orthogonal projections to latent structures (OPLS) analyses. Out of 95 identified metabolites, 19 were significantly decreased due to the intervention after 18 months of intervention. Significant changes could be seen in the pentose phosphate, the mevalonate, the beta-oxidation and the xanthine oxidase pathways. The intervention also resulted in changes in the urea cycle, and increases in the levels of the precursors to neurotransmitters of the brain. This adds information to previous published results reporting decreased oxidative stress and inflammation. This is the first-time metabolic profiling has been applied to elucidate the mechanisms behind an intervention with SeQ10. The study is small and should be regarded as hypothesis-generating; however, the results are interesting and, therefore, further research in the area is needed. This study was registered at Clinicaltrials.gov, with the identifier NCT01443780.

1. Introduction

Selenium (Se) is a trace element that is essential for all living cells within the human body. There are 25 genes encoding for selenoproteins in humans [1]. Several of the selenoproteins have important roles in anti-oxidative stress defence, and in immune surveillance [2,3]. Some of the most important selenoproteins are the glutathione peroxidases, thioredoxin reductase 1, selenoprotein P, and iodothyronine deiodinases [4]. There is a relation between the Se content in soil, food grown on these soils and dietary Se intake, where both high and low intakes have been reported [5,6]. In Europe there is in general a low concentration of Se in the soil [7,8,9,10,11] and, hence, a low intake, estimated at about 40 µg Se/day [1]. The daily intake of selenium needed for the optimal expression of selenoprotein P in plasma has been debated. Xia et al. demonstrated that an intake of at least 75 µg/day, to achieve a plasma concentration of Se above 90 µg/L was required in normal weight Caucasians [12], whereas Hurst et al. reported that an intake of at least 105 µg/day was required [13]. In Sweden, reports have shown an intake of selenium below the recommended levels [14,15], which appears to increase the cardiovascular risk [16]. The soils in North America have a significantly higher content of selenium than in Europe, with reported serum selenium levels generally above 120 μg/L in the US population [17,18]. In a population of elderly healthy Swedish persons with a mean serum selenium concentration of 67 μg/L, we showed that those with low selenium concentration had a higher cardiovascular risk compared to those with higher selenium concentration, thus confirming previous reports in the literature [16,19].
There is an important relationship between Se and ubiquinone (coenzyme Q10, hereafter abbreviated Q10) in living cells. In order to synthesise selenocysteine-containing enzymes a functional mevalonate pathway is required, and Q10 is a product in this pathway [20]. Additionally, the selenoenzyme TrxR may reduce Q10 to its active form, ubiquinol; hence, the presence of an adequate concentration of selenium is required in the cell for optimal Q10 function.
In addition to being a strong antioxidant and essential for all living cells, Q10 is of particular importance for the crucial ATP-generating steps within the mitochondrial respiratory chain and in oxidative phosphorylation mediated through the pentose phosphate pathway [21]. The endogenous production of Q10 decreases with increasing age, and in the myocard, by the age of 80 years it is reduced to half the level that exists at 20 years [22]. Furthermore, the need for Q10 is even greater during diseases [23,24].
In 2013, we presented the results of four years of intervention using Se and Q10 combined (SeQ10) as a dietary supplement in an elderly healthy population in a randomised prospective double-blind, placebo-controlled study [25]. The main result of the intervention was a significant reduction of cardiovascular mortality by more than 50%, better cardiac function as demonstrated on echocardiography, and a lower plasma concentration of the natriuretic peptide N-terminal fragment of proBNP (NT-proBNP) indicating less myocardial wall tension [25]. Due to the striking positive clinical results, in an effort to shed light on possible mechanisms behind the clinical effects, we used a metabolic profiling approach to pinpoint changes in metabolites and metabolic pathways following SeQ10 administration.
Studies of metabolic patterns, by analysing body fluids or sampled tissue, have evolved into an important routine in several fields of the clinical and pharmacological discipline [26]. Non-targeted metabolite profiling—known as metabolomics—provides a functional molecular picture of the biochemistry connecting both genomics and proteomics to a physiological state, and it also incorporates environmental interactions. Metabolomics is, thus, a mapping of the actual state of the phenotype.
The metabolomic approach in characterising biological systems is frequently applied in several situations; (1) In the diagnostic process, i.e., where the metabolic profiles of healthy subjects, and subjects with a disease are used for classification of new cases. (2) In the identification of perturbed metabolic pathways, i.e., the detection of biomarkers and connecting them to, e.g., disease pathogenesis, and (3) as a tool in disease follow-up and treatment, i.e., following changes of metabolites after changes in lifestyle and/or food and/or pharmacological intervention [27,28,29,30].
In the present work, we used metabolomics to identify the impacts of SeQ10 intervention on the plasma metabolome. As previous results from this intervention study indicated that males and females differed strongly in cardiovascular mortality, and as we could not exclude a hormonal influence, we decided to include only males in this study in order to have a more homogenous group to analyse.

2. Results

The baseline characteristics of the total male study population are presented in Table 1.
As can be seen, all variables were balanced between the two groups included in this evaluation, with the exception that among those receiving a placebo there were more participants in the NYHA functional group II and fewer in group I.

2.1. Main Study

The active treatment group was evaluated against the placebo group using the OPLS-DA model assessment at 18 months (T18) where 95 metabolites analysed in plasma were included. This evaluation demonstrated significant differences between the two groups by using seven-fold cross-validation (the OPLS-DA model with one significant component that explained 57% of the variance in Y, R2X = 0.19, R2(cum) = 0.57 and Q2(cum) = 0.18) (Figure 1).
One of the samples from the placebo group had a position in the plot that indicated scores that deviated from those of the rest of the placebo group and, thus, negatively influenced the power of the OPLS-DA model used. After checking with this individual, it was revealed that despite being a participant of the placebo group, the person with sample no. 186 supplemented himself with selenium on his own initiative, which the model including each individual’s metabolite pattern in plasma, could detect and demonstrate (Figure 1).
This individual was thus excluded from all further modelling, resulting in the final OPLS-DA model with one predictive and one orthogonal component including 17 individuals (R2X = 0.13, R2(cum) = 0.91, Q2(cum) = 0.41). The difference in the whole metabolome between the two studied groups was evident when comparing predicted cross-validated Y values for samples from both groups (p < 0.01) (Figure 2).
Using the jack-knifing method in the multivariate analyses and applying a 95% confidence interval, the results showed that 19 metabolites out of the 95 identified and modelled metabolites were significantly decreased in the SeQ10-treated men. A further examination applying a T-test analysis, showed that 13 of these 19 metabolites were significantly lower in the SeQ10 -treated men when using univariate statistics (Table 2).
No metabolites were significantly higher in the SeQ10 men as compared with the controls. The p(corr) loading plot from the OPLS-DA model indicated that the difference could largely be attributed to changes in amino acids, which were at lower levels in the SeQ10-treated group in comparison with the placebo group (Figure 3).
However, a decrease could also be seen in the polyunsaturated fatty acids arachidonic acid, eicosapenatenoic and docosahexaenoic acid, and an increase in the saturated fatty acids (hexadecanoic acid, myristic and stearic acid) and monounsaturated fatty acids hexadecenoic acid and oleic acid.

2.2. Projection of Deceased Participants

Two individuals that participated in the project for at least 18 months, but later died from cardiovascular disease, delivered samples at T0 and T18. The metabolic profiles from these individuals were projected into the principal component analysis (PCA) models to investigate whether the T18 blood samples from these individuals deviated from the other T18 samples. In the score plots (not shown) it could be seen that the metabolic patterns of the samples from the deceased participants did not deviate from the surviving participants. It was, therefore, not possible to use the existing models to predict the fate of individual participants after T18.

2.3. Validation Studies

Both univariate analyses and multivariate modelling were carried out in each batch separately. The results are summarised in Table 3 and Table 4 where the T-tests with corresponding variables’ importance to the projection (VIP) values from the OPLS-DA models are presented, and the values of the combined loading vectors obtained are presented in Table 5. In that table, 62 out of 107 (57.9%) evaluated metabolites showed congruence in change, or no change between the main study and at least one of the validation studies.

2.4. Validation Studies 1 and 2

Whereas no metabolites by univariate analysis were significantly changed in any of the studies, comparison of metabolic profiles in relation to the SeQ10 administration supports the results from the main study at T18. There is also consistency between the two validation studies, where 62% of the metabolites detected in at least two studies showed the same changes (higher or lower metabolite levels) in SeQ10-treated individuals as compared with those receiving a placebo. The differences between the groups were largely attributed to amino acids that were at lower levels in the actively treated group compared with the placebo group (Table 3).

3. Discussion

Our group has recently presented clinical data showing reduced cardiovascular mortality, increased cardiac function according to echocardiographic evaluation, and decreased concentration of the cardiac biomarker NT-proBNP, as a result of intervention with SeQ10 as a dietary supplement given to an elderly healthy population over four years. To acquire knowledge of the possible biochemical mechanisms behind these clinical results, we have applied metabolic profiling to investigate possible metabolic changes due to the intervention.
We have shown an effect on the metabolome in the SeQ10-treated men after 18 months when compared with those receiving the placebo (Figure 1, Figure 2 and Figure 3, Table 2 and Table 3).
For a majority of the changes observed in the main study (57.9%), the same changes could be seen in one or both validation studies, supporting the validity of the results from the main study. Together, the changes in the metabolite patterns after SeQ10 treatment could be attributed to biochemical pathways discussed further in the following paragraphs.
An interesting example of the effect of the intervention could be seen in the level of 1,5-anhydro-D-glucitol, a monosaccharide that is found in many foods. This substance is reabsorbed by the proximal renal tubule in competition with glucose. A high blood concentration of 1,5-anhydro-D-glucito is indicative of a low level of glucose in primary urine, and is therefore regarded as positive [31]. It has also been reported in the literature that low levels of 15-anhydro-D-glycitol are associated with increased cardiovascular mortality or increased coronary artery disease risk [32], even in cases without signs of diabetes [33]. In our study we reported a higher concentration of the substance in those receiving SeQ10, compared with those on placebo, indicating a change in the metabolomic profile that could be a clinically important effect of the intervention.

3.1. Biochemical Pathways Altered by the Intervention with SeQ10

(A) The pentose phosphate pathway, which is an alternative metabolic pathway to that of glycolysis with generation of, e.g., NADPH and 5-carbon sugars (pentoses), was enhanced. In those receiving SeQ10, a decrease in 1,6-anhydro-β-D-glucose, as well as an increase in glycerol-3-phosphate, fructose and inositol-1-phosphate could be seen. This strongly indicates that the equilibria are shifted towards the generation of reducing equivalents in the form of NADPH, i.e., further providing selenoproteins, e.g., thioredoxin reductase, glutathione peroxidases and selenop, that may scavenge reactive oxygen species (ROS) within the cell. Previously, increased oxidative stress with ROS formation has been discussed as one of the main mechanisms in the development of heart disease [34]; therefore, intervention to decrease oxidative stress is of great interest [35]. Our group have previously reported signs of less oxidative stress as a result of the SeQ10 intervention [36].
(B) The mevalonate pathway, in which squalene is an important precursor in the biosynthesis of cholesterol, heme, vitamin K, coenzyme Q and all steroid hormones [37], was influenced by the intervention with selenium and coenzyme Q10. We observed that several phytosterols e.g., campesterol, and sitosterol showed higher levels in the SeQ10-treated participants than in those receiving the placebo. Usually, phytosterols are exogenous compounds produced by plants and found in the diet. Some phytosterols may be added to food (usually margarine), and they compete with cholesterol as regards intestinal absorption and can thereby lower LDL cholesterol [38]. It has been reported that phytosterols reduce LDL levels and increase HDL in the blood and, thus, have a potential place in the treatment of hyperlipidaemia [39,40,41]. A potential source of the phytosterols in the present study could be the yeast in the tablets containing yeast selenium. One might also speculate that the intervention modified the gastro-intestinal flora, resulting in increased intestinal release of phytosterols from the food matrix or uptake of such compounds [42].
(C) The beta-oxidation pathway is increased by the treatment with SeQ10. This is supported by the decrease of unsaturated fatty acids with an even number of double bonds and the increased levels of monounsaturated fatty acids (Figure 3). Thus, fatty acids with multiple double bonds decrease, whereas less unsaturated fatty acids increase as a result of the intervention. In combination with the cardio-protective steroids produced in the mevalonate pathway that could be demonstrated in the analyses in the present study, the altered fatty acid composition indicates that SeQ10 should have a positive effect on individuals with metabolic syndrome and/or type 2 diabetes [43].
(D) The xanthine oxidase pathway, which is a metabolic pathway for uric acid formation, is down-regulated as the levels of uric acid in plasma are decreased in SeQ10 -treated men. Uric acid is a major antioxidant in a hydrophilic environment. Reduced synthesis might indicate a reduced need for antioxidants because ROS are taken care of by selenoproteins and/or the reduced form of Q10 in the supplemented group. This supports our previous results that SeQ10 appeared to have an anti-inflammatory effect [44]. From a clinical point of view, reports indicating an important association between inflammation and vascular disease can be found in the literature. Therefore, this aspect is important [45,46]. An anti-oxidative effect caused by the intervention with SeQ10 and shown by the decrease in uric acid levels is also supported by previously presented findings in the same study population on the concentrations of copeptin and the midregional segment of proadrenomedullin (MR-proADM) as markers for anti-oxidative activity [36].
(E) An interesting observation is that an influence on the urea cycle could be seen, with decreased levels of non-essential amino acids as well as creatinine, but without increased levels of the corresponding ketoacids. It could be inferred that the non-essential amino acids were part of the production of proteins. This suggests that the treatment with SeQ10 favours an anabolic process instead of an oxidative process with catabolism.
Decreased levels of the essential amino acids could also be seen in those who were supplemented with SeQ10, compared with those on placebo. As it is reported in the literature that a high concentration of some essential, branched-chain amino acids (leucine, isoleucine and valine) is associated with diabetes mellitus [47,48], and also with oxidative stress and inflammation [49], decreased levels could be interpreted as positive.

3.2. General Effects of SeQ10 Treatment

In addition to the changes discussed above we also observed that the levels of the precursors to neurotransmitters of the brain, i.e., the aromatic amino acids, were decreased in those receiving SeQ10. This might explain the observed reduction in sense of fatigue in the intervention group reported previously [50] The same has also been reported for Q10 given in monotherapy [51,52].
An interesting observation is that the participants treated with SeQ10 had changes in their metabolome that were opposite to changes in patients suffering from severe systemic lupus erythematosus (SLE) [53]. Therefore, patients with severe SLE might also benefit from intervention with SeQ10, as this could normalise their pathological metabolic profile.
At T18 months we also observed that the SeQ10-treated men had lower levels of salicylic acid in their plasma, which could imply less consumption of salicylic acid-containing painkillers than the placebo group. However, we do not have any anamnestic data to support this interesting observation, even though there are reports in the literature indicating such an effect by coenzyme Q10 [54,55,56].
Another observation was decreased levels of the essential amino acid methionine in the SeQ10-treated group. Methionine is the main precursor for cysteine and glutathione and a central molecule in one-carbon metabolism as a precursor for S-adenosylmethionine (SAM), the principal donor of methyl groups. In mice, selenium modulates one-carbon metabolism and increases hepatic methylation capacity, resulting in increased global DNA methylation and affecting lipogenesis [43,57]. The DNA methylation regulating gene expression is also related to ageing [58]. In a previous publication our group presented data indicating less fibrosis in seven out of eight evaluated biomarkers for fibrosis in the population receiving SeQ10 compared with those on placebo. This could indicate a possible anti-ageing effect on the cardiovascular system, as the ageing process is highly associated with increased fibrosis of the cardiovascular system [59]. Increased expression and synthesis of selenoprotein P following selenium supplementation requires SAM-dependent protein methylation [57].
Taken together, the presented changes in metabolic profiles as seen in the present evaluation could provide important information to support the positive clinical results that have been demonstrated in the previous publications. Therefore, for those with insufficient selenium and coenzyme Q10 intake, the metabolic profiling results could strengthen the argument for supplementation.

4. Material

4.1. Aim

The aim of the present work was to investigate changes in biochemical pathways after intervention with SeQ10 to a healthy elderly male community-living population by metabolic profiling of plasma metabolites.

4.2. Design

A double-blind randomised placebo-controlled study where intervention with a dietary supplement of Se and Q10 combined (SeQ10) in an elderly healthy rural population of 443 individuals in the age range of 70–88 years has been previously reported [25,60]. The first participant was included in January 2003, and the last participant concluded the study in February 2010. The participants received the intervention for 48 months during which time they were re-examined every six months, and throughout the study, all mortality was registered. During the study 221 individuals received active supplementation of 200 μg/day organic Se (SelenoPrecise®, PharmaNord, Denmark), plus 200 mg/day of Q10 (Bio-Quinon®, PharmaNord, Denmark), and 222 individuals received a placebo. At inclusion, all participants underwent clinical examination, a new patient record was obtained, the New York Heart Association (NYHA) functional class was assessed (this classification grades how a patient with heart disease experiences symptoms of tiredness, breathlessness or chest pain; it is graded from I–IV, where IV is symptoms already at rest), an ECG and a Doppler-echocardiography were performed. Informed consent was obtained from each patient. The study was approved by the Regional Ethical Committee in Linköping (Forskningsetikkommmitten, Hälsouniversitetet, SE-581 85 Linköping, Sweden; no. D03-176), and it conforms to the ethical guidelines of the 1975 Declaration of Helsinki. The Medical Product Agency declined to review the study protocol since the study was not considered a trial of a medication for a certain disease, but rather one of food supplement commodities that are commercially available. The study was registered at clinicaltrials.gov. (NCT01443780).

5. Methods

5.1. Echocardiography

Doppler-echocardiography examinations were performed with the participants in the left lateral position. The ejection fraction (EF) readings were categorised into four classes with interclass limits placed at 30%, 40% and 50% [61,62]. Normal systolic function was defined as EF ≥ 50%, while severely impaired systolic function was defined as EF < 30%.

5.2. Blood Samples

Blood samples were collected while the participants were resting in a supine position. Pre-chilled, EDTA vials (Terumo EDTA K-3) were used. The vials were placed on ice and centrifuged for 10 min at 3000× g, +4 °C, and the plasma was immediately transferred to storage at −70 °C until further analysis. No sample was thawed more than once.

5.3. Sample Selection

The available biobank contained samples from 443 individuals, including both males and females, and the male population was chosen. From those, the individuals who delivered blood samples at one (T0), two (T0 + T18 months), or three time points (T0 + T18 months + T48 months) were chosen after taking into account clinical and demographic variables in order to get balanced groups. The samples were divided into a main study as well as two validation studies. In each of the validation studies the samples were divided into analytical batches run on separate days. Each batch contained a balanced mix of treatment and placebo samples, with all time point samples from an individual included in the same batch. As explained below, samples selected for analysis were representative of the whole study and covered the whole expected intervention period, the same was true for batch separation where each batch itself contained representative information found in the whole cohort [63,64]. Therefore, the individuals included in the study were selected to meet those criteria.

5.4. Main Study

In the participant database, clinical and demographic data were used for the selection of participants. From the original dataset, a selection was made by use of multivariate design (MVD), which, in turn, was based on a combination of an experimental design and multivariate data analysis (MVDA) [63,64]. By using MVD we ensured that the selected subjects together spanned as large an area of the available properties domain as possible, thus including as large a natural variation as possible, and at the same time this ensured that the properties’ domain in the placebo group covered the same domain as the group receiving SeQ10. Principal component analysis of the clinical and demographic information that was collected on the inclusion date of each participant was applied. A total of 182 variables representing all available demographic and clinical data were used in the PCA model. Based on two separate overview models of the SeQ10 and the placebo-treated men, a representative selection for an initial metabolomic study was obtained. From the PCA score plot (which plotted the subjects in the first two principal components) 10 individuals were selected from each group (two from each of the four quadrants plus two from the centre, thus spanning the whole properties domain). In addition to these 20 participants, seven others (two SeQ10 and five placebo) were selected representing the participants that had participated up until T18 but had later died of cardiovascular disease (CVD). Those who died a cardiovascular death (n = 2) during the intervention time were not included in the main analyses but were evaluated separately.
When retrieving the samples from the freezer it was noticed that one of the placebo samples was contaminated with hepatitis and this was therefore excluded from analysis. During sample preparation, one of the SeQ10 samples did not derivatise so no metabolite analysis could be performed.
All derivatised samples were analysed by GC-MS, and 315 metabolites were detected. Using metabolite databases, 95 of these were identified and were used in the modelling of metabolite profiles.
Ultimately, we included nine male SeQ10 and nine placebo participants (T18) in the main study, and we treated the samples from the individuals that died from CVD as a separate group.

5.5. Validation Studies

To confirm the results from the main study, it was decided to perform two validation studies; “Validation study 1” and “Validation study 2”, each made up of several batches. Each batch contained a balanced mix of SeQ10 and placebo samples, selected according to the MVD approach (as described for the main study).

5.5.1. Validation Study 1

Metabolic profiles from male individuals were investigated. A total of 60 male individuals (31 SeQ10, 29 placebo) were included, all having three samples each (T0, T18 and T48). One individual (SeQ10) from the first study with two time points was included as well. All samples belonging to Validation study 1 were divided into five analytical batches (subsets).

5.5.2. Validation Study 2

In validation study 2, samples from 37 men (23 SeQ10, 14 placebo) were included.

5.6. Metabolomic Profiling

5.6.1. Metabolite Analysis (Extraction, Quantification and Identification)

Plasma samples were extracted, derivatised, and analysed using GC–TOF–MS as previously described [65]. Masses were acquired in the mass range 50–800 mass/charge numbers of ions (m/z) at a rate of 30 spectra/s. The acceleration voltage was turned on after a solvent delay of 150 s. An alkane series (C10–C32) was run together with all samples. Non-processed MS files from GC/TOFMS analysis were exported in NetCDF format to MATLAB software 7.3 (Mathworks, Natick, MA, USA), where all data pre-treatment procedures, such as baseline correction, chromatogram alignment, time-window setting and multivariate curve resolution were performed using custom scripts [66].
For the identification of metabolites, The National Institute of Standards and Technology (NIST) MS Search software (version 2.0, Gaithersburg, MD, USA) was used. With this we compared the mass spectra of all detected compounds with spectra in the NIST library, the in-house mass spectra library established by Umeå Plant Science Centre and Swedish Metabolomics Centre and the mass spectra library maintained by the Max Planck Institute, Golm, Germany (http://csbdb.mpimp-golm.mpg.de/csbdb/gmd/gmd.html). A retention index comparison was performed, with a retention index deviation < ±10 (in addition to a high spectral match, above 800).

5.6.2. Statistical Analysis

Univariate Statistics

Univariate analysis was performed on each metabolite separately, with p values calculated by performing two-sample, unequal variance, two-tailed t-tests.

Multivariate Data Analysis

Metabolite levels reflect biochemical processes, and as few processes in an organism are truly independent, metabolite levels are also highly correlated. The MVDA approaches are best suited to analyses of highly correlated data and allow the identification of the following: small changes in the variable patterns between different groups; internal relations/correlations between samples and between metabolites; and how samples and metabolites are connected to each other. These differences can be missed by using a univariate approach alone [67,68,69].

OPLS-DA Models

The orthogonal partial least-squares discriminant analysis (OPLS-DA) models were calculated between the two intervention groups (SeQ10 and placebo) to identify the metabolites that showed different levels in these groups.

Data Normalisation

Data were normalised using 11 internal standards (eluting over the whole chromatographic time range). To obtain accurate peak areas for the internal standards, two unique masses for each compound were specified and the samples were reprocessed using an in-house MATLAB-based script. Principal component analyses were performed on the peak areas of internal standards with Unit Variance Normalisation and the first component (t1) score value for each sample was used to normalise the resolved data. This was done by dividing the peak areas in each sample by the corresponding score value [70].
All multivariate modelling was performed using the software SIMCA version 14 (Umetrics, Umeå, Sweden). The PCA approach was used for a metabolomics data overview [71], and the OPLS-DA approach was used to elucidate the differences concerning metabolomics between the two groups of subjects (SeQ10 and placebo) [72,73]. Unit Variance scaling was used for both types of models [74]. The significant differences between the SeQ10 group and the placebo-treated individuals were identified by using confidence intervals calculation by use of jack-knifing [75]. A potential model significance was found by means of cross-validation [75]. Furthermore, cross-validated predicted Y values from the cross-validation procedure were used to visualise the ability of the OPLS-DA model to separate individuals receiving SeQ10 from those receiving a placebo.
The VIP values were calculated (SIMCA v.14) as VIP scores. This was done in order to estimate the importance of each variable in the projection used in the OPLS-DA model. A variable with a VIP score of about 0.8 or greater was considered significant in the model.
In order to obtain an overview model of the different batches, the p(corr)-loadings for each of the SeQ10 treated vs. placebo OPLS-DA models were combined into an X-matrix which was used to create a PCA model where no scaling was applied. The p [1] loading profile from this PCA model was used as a combined metabolic profile representing the metabolic effects of the SeQ10 intake in validation studies 1 and 2, and compared with the p(corr) vector from the main study.
Small differences in the number of significant compounds when detected by univariate and multivariate approaches were expected, since both approaches are based on different testing principles; the former (t-tests) is a statistical hypothesis test, whereas the multivariate testing, using jack-knifing, is a resampling technique applied to the multivariate space. These differences were, thus, included in the evaluation of the results.

6. Conclusions

There was a clear difference in the metabolic profile in plasma from men receiving SeQ10 as compared with men receiving placebo. The major differences can be seen in the pentose phosphate pathway, in the mevalonate pathway, and in the beta-oxidation pathway. An interesting observation was downregulation of plasma methionine, indicating increased methylating activity, as also previously observed experimentally. It can also be reported that xanthine oxidase was down-regulated. Moreover, there were effects on the urea cycle where the non-essential amino acids were down-regulated. All these effects, together with signs of less inflammation and less oxidative stress observed in the previous studies of this group, indicate a positive health effect by the supplementation with selenium and coenzyme Q10 compared with those on placebo.
However, as this is the first intervention with the SeQ10 combination in an elderly population low in Se, the results should be regarded as hypothesis-generating, and more research is needed to confirm the positive health impact of the SeQ10 combination in elderly persons.

Author Contributions

U.A. and P.J. conceived and designed the research project. U.A. and P.J. conducted the research. T.L. and I.S. provided the essential reagents, and performed the analyses. U.A., T.L. and K.L.-E. analysed data and performed the statistical analyses. U.A., P.J., J.A. (Jan Aaseth), J.A. (Jan Alexander), T.L., I.S. and K.L.-E. wrote the paper. U.A. had the final responsibility for the final content. All the authors have read and approved the final manuscript.

Funding

The analysis costs were partially supported by grants from Pharma Nord Aps, Denmark, the County Council of Östergötland, Linköping University. The funding organisations had no role in the design, management, analysis, or interpretation of the data, nor in the preparation, review or approval of the manuscript. No economic compensation was distributed.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Rayman, M.P. Selenium and human health. Lancet 2012, 379, 1256–1268. [Google Scholar] [CrossRef]
  2. Selenius, M.; Rundlöf, A.K.; Olm, E.; Fernandes, A.P.; Björnstedt, M. Selenium and the Selenoprotein Thioredoxin Reductase in the Prevention, Treatment and Diagnostics of Cancer. Antioxid. Redox Signal. 2010, 12, 867–880. [Google Scholar] [CrossRef] [PubMed]
  3. Huang, Z.; Rose, A.H.; Hoffmann, P.R. The Role of Selenium in Inflammation and Immunity: From Molecular Mechanisms to Therapeutic Opportunities. Antioxid. Redox Signal. 2012, 16, 705–743. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Fairweather-Tait, S.J.; Bao, Y.; Broadley, M.R.; Collings, R.; Ford, D.; Hesketh, J.E.; Hurst, R. Selenium in Human Health and Disease. Antioxid. Redox Signal. 2011, 14, 1337–1383. [Google Scholar] [CrossRef]
  5. Yang, G.; Zhou, R. Further observations on the human maximum safe dietary selenium intake in a seleniferous area of China. J. Trace Elem. Electrolytes Health Dis. 1994, 8, 159–165. [Google Scholar]
  6. Cheng, Y.Y.; Qian, P.C. The effect of selenium-fortified table salt in the prevention of Keshan disease on a population of 1.05 million. Biomed. Environ. Sci. 1990, 3, 422–428. [Google Scholar]
  7. Akbaraly, N.T.; Arnaud, J.; Hininger-Favier, I.; Gourlet, V.; Roussel, A.-M.; Berr, C. Selenium and Mortality in the Elderly: Results from the EVA Study. Clin. Chem. 2005, 51, 2117–2123. [Google Scholar] [CrossRef] [Green Version]
  8. Adame, E.M.; Florea, D.; Pérez, L.S.; López, J.M.; López-González, B.; De La Cruz, A.P.; Del Pozo, E.P. Deficient selenium status of a healthy adult Spanish population. Nutrición Hospitalaria 2012, 27, 524–528. [Google Scholar]
  9. Letsiou, S.; Nomikos, T.; Panagiotakos, D.; Pergantis, S.A.; Fragopoulou, E.; Antonopoulou, S.; Pitsavos, C.; Stefanadis, C. Serum total selenium status in Greek adults and its relation to age. The ATTICA study cohort. Biol. Trace Elem. Res. 2009, 128, 8–17. [Google Scholar] [CrossRef]
  10. Benes, B.; Spevácková, V.; Smíd, J.; Batáriová, A.; Cejchanová, M.; Zítková, L. Effects of age, BMI, smoking and contraception on levels of Cu, Se and Zn in the blood of the population in the Czech Republic. Central Eur. J. Public Health 2005, 13, 202–207. [Google Scholar]
  11. Reszka, E.; Wieczorek, E.; Jablonska, E.; Janasik, B.; Fendler, W.; Wasowicz, W. Association between plasma selenium level and NRF2 target genes expression in humans. J. Trace Elem. Med. Biol. 2015, 30, 102–106. [Google Scholar] [CrossRef] [PubMed]
  12. Xia, Y.; Hill, K.E.; Li, P.; Xu, J.; Zhou, D.; Motley, A.K.; Wang, L.; Byrne, D.W.; Burk, R.F. Optimization of selenoprotein P and other plasma selenium biomarkers for the assessment of the selenium nutritional requirement: A placebo-controlled, double-blind study of selenomethionine supplementation in selenium-deficient Chinese subjects. Am. J. Clin. Nutr. 2010, 92, 525–531. [Google Scholar] [CrossRef] [PubMed]
  13. Hurst, R.; Armah, C.N.; Dainty, J.R.; Hart, D.J.; Teucher, B.; Goldson, A.J.; Broadley, M.R.; Motley, A.K.; Fairweather-Tait, S.J. Establishing optimal selenium status: Results of a randomized, double-blind, placebo-controlled trial. Am. J. Clin. Nutr. 2010, 91, 923–931. [Google Scholar] [CrossRef] [PubMed]
  14. Gao, H.; Hägg, S.; Sjögren, P.; Lambert, P.C.; Ingelsson, E.; Van Dam, R.M. Serum selenium in relation to measures of glucose metabolism and incidence of Type 2 diabetes in an older Swedish population. Diabet. Med. 2014, 31, 787–793. [Google Scholar] [CrossRef] [PubMed]
  15. Harris, H.R.; Bergkvist, L.; Wolk, A. Selenium intake and breast cancer mortality in a cohort of Swedish women. Breast Cancer Res. Treat. 2012, 134, 1269–1277. [Google Scholar] [CrossRef] [PubMed]
  16. Alehagen, U.; Johansson, P.; Björnstedt, M.; Rosén, A.; Post, C.; Aaseth, J. Relatively high mortality risk in elderly Swedish subjects with low selenium status. Eur. J. Clin. Nutr. 2016, 70, 91–96. [Google Scholar] [CrossRef] [PubMed]
  17. Kafai, M.R.; Ganji, V. Sex, age, geographical location, smoking, and alcohol consumption influence serum selenium concentrations in the USA: Third national health and nutrition examination survey, 1988–1994. J. Trace Elements Med. Boil. 2003, 17, 13–18. [Google Scholar] [CrossRef]
  18. Bleys, J.; Navas-Acien, A.; Laclaustra, M.; Pastor-Barriuso, R.; Menke, A.; Ordovas, J.; Stranges, S.; Guallar, E. Serum Selenium and Peripheral Arterial Disease: Results From the National Health and Nutrition Examination Survey, 2003–2004. Am. J. Epidemiol. 2009, 169, 996–1003. [Google Scholar] [CrossRef] [PubMed]
  19. Alehagen, U.; Alexander, J.; Aaseth, J. Supplementation with Selenium and Coenzyme Q10 Reduces Cardiovascular Mortality in Elderly with Low Selenium Status. A Secondary Analysis of a Randomised Clinical Trial. PLoS ONE 2016, 11, e0157541. [Google Scholar] [CrossRef]
  20. Xia, L.; Nordman, T.; Olsson, J.M.; Damdimopoulos, A.; Björkhem-Bergman, L.; Nalvarte, I.; Eriksson, L.C.; Arnér, E.S.; Spyrou, G.; Björnstedt, M. The mammalian cytosolic selenoenzyme thioredoxin reductase reduces ubiquinone. A novel mechanism for defense against oxidative stress. J. Biol. Chem. 2003, 278, 2141–2146. [Google Scholar] [CrossRef]
  21. Acosta, M.J.; Fonseca, L.V.; Desbats, M.A.; Cerqua, C.; Zordan, R.; Trevisson, E.; Salviati, L. Coenzyme Q biosynthesis in health and disease. Biochim. Biophys. Acta 2016, 1857, 1079–1085. [Google Scholar] [CrossRef] [PubMed]
  22. Kalén, A.; Appelkvist, E.-L.; Dallner, G. Age-related changes in the lipid compositions of rat and human tissues. Lipids 1989, 24, 579–584. [Google Scholar] [CrossRef] [PubMed]
  23. Folkers, K.; Moesgaard, S.; Morita, M. A one year bioavailability study of coenzyme Q10 with 3 months withdrawal period. Mol. Asp. Med. 1994, 15, s281–s285. [Google Scholar] [CrossRef]
  24. Molyneux, S.L.; Florkowski, C.M.; George, P.M.; Pilbrow, A.P.; Frampton, C.M.; Lever, M.; Richards, A.M. Coenzyme Q10: An independent predictor of mortality in chronic heart failure. J. Am. Coll. Cardiol. 2008, 52, 1435–1441. [Google Scholar] [CrossRef]
  25. Alehagen, U.; Johansson, P.; Björnstedt, M.; Rosén, A.; Dahlström, U. Cardiovascular mortality and N-terminal-proBNP reduced after combined selenium and coenzyme Q10 supplementation: A 5-year prospective randomized double-blind placebo-controlled trial among elderly Swedish citizens. Int. J. Cardiol. 2013, 167, 1860–1866. [Google Scholar] [CrossRef] [Green Version]
  26. Madsen, R.K.; Lundstedt, T.; Gabrielsson, J.; Sennbro, C.-J.; Alenius, G.-M.; Moritz, T.; Rantapää-Dahlqvist, S.; Trygg, J. Diagnostic properties of metabolic perturbations in rheumatoid arthritis. Arthritis Res. Ther. 2011, 13, R19. [Google Scholar] [CrossRef]
  27. Madsen, R.; Lundstedt, T.; Trygg, J. Chemometrics in metabolomics—A review in human disease diagnosis. Anal. Chim. Acta 2010, 659, 23–33. [Google Scholar] [CrossRef]
  28. Nordstrom, A.; Lewensohn, R. Metabolomics: Moving to the clinic. J. Neuroimmune Pharmacol. 2010, 5, 4–17. [Google Scholar] [CrossRef]
  29. Vinayavekhin, N.; Homan, E.A.; Saghatelian, A. Exploring disease through metabolomics. ACS Chem. Biol. 2010, 5, 91–103. [Google Scholar] [CrossRef]
  30. Scholz, S.; Fischer, S.; Gündel, U.; Küster, E.; Luckenbach, T.; Voelker, D. The zebrafish embryo model in environmental risk assessment—applications beyond acute toxicity testing. Environ. Sci. Pollut. Res. 2008, 15, 394–404. [Google Scholar] [CrossRef]
  31. Yamanouchi, T.; Minoda, S.; Yabuuchi, M.; Akanuma, Y.; Akanuma, H.; Miyashita, H.; Akaoka, I. Plasma 1,5-anhydro-D-glucitol as new clinical marker of glycemic control in NIDDM patients. Diabetes 1989, 38, 723–729. [Google Scholar] [CrossRef] [PubMed]
  32. Fujiwara, T.; Yoshida, M.; Yamada, H.; Tsukui, T.; Nakamura, T.; Sakakura, K.; Wada, H.; Arao, K.; Katayama, T.; Funayama, H.; et al. Lower 1,5-anhydroglucitol is associated with denovo coronary artery disease in patients at high cardiovascular risk. Heart Vessels 2015, 30, 469–476. [Google Scholar] [CrossRef]
  33. Ouchi, S.; Shimada, K.; Miyazaki, T.; Takahashi, S.; Sugita, Y.; Shimizu, M.; Murata, A.; Kadoguchi, T.; Kato, T.; Aikawa, T.; et al. Low 1,5-anhydroglucitol levels are associated with long-term cardiac mortality in acute coronary syndrome patients with hemoglobin A1c levels less than 7.0. Cardiovasc. Diabetol. 2017, 16, 151. [Google Scholar] [CrossRef] [PubMed]
  34. Münzel, T.; Camici, G.G.; Maack, C.; Bonetti, N.R.; Fuster, V.; Kovacic, J.C. Impact of Oxidative Stress on the Heart and Vasculature: Part 2 of a 3-Part Series. J. Am. Coll. Cardiol. 2017, 70, 212–229. [Google Scholar] [CrossRef] [PubMed]
  35. Siti, H.N.; Kamisah, Y.; Kamsiah, J. The role of oxidative stress, antioxidants and vascular inflammation in cardiovascular disease (a review). Vasc. Pharmacol. 2015, 71, 40–56. [Google Scholar] [CrossRef] [PubMed]
  36. Alehagen, U.; Aaseth, J.; Johansson, P. Less increase of copeptin and MR-proADM due to intervention with selenium and coenzyme Q10 combined: Results from a 4-year prospective randomized double-blind placebo-controlled trial among elderly Swedish citizens. BioFactors 2015, 41, 443–452. [Google Scholar] [CrossRef]
  37. Holstein, S.A.; Hohl, R.J. Isoprenoids: Remarkable diversity of form and function. Lipids 2004, 39, 293–309. [Google Scholar] [CrossRef]
  38. Shaghaghi, M.A.; AbuMweis, S.S.; Jones, P.J. Cholesterol-Lowering Efficacy of Plant Sterols/Stanols Provided in Capsule and Tablet Formats: Results of a Systematic Review and Meta-Analysis. J. Acad. Nutr. Diet. 2013, 113, 1494–1503. [Google Scholar] [CrossRef]
  39. Marangoni, F.; Poli, A. Phytosterols and cardiovascular health. Pharmacol. Res. 2010, 61, 193–199. [Google Scholar] [CrossRef]
  40. Dumolt, J.H.; Rideout, T.C. The Lipid-lowering Effects and Associated Mechanisms of Dietary Phytosterol Supplementation. Curr. Pharm. Des. 2017, 23, 5077–5085. [Google Scholar] [CrossRef]
  41. Cabral, C.E.; Klein, M.R.S.T. Phytosterols in the Treatment of Hypercholesterolemia and Prevention of Cardiovascular Diseases. Arq. Bras. Cardiol. 2017, 109, 475–482. [Google Scholar] [CrossRef] [PubMed]
  42. Shandilya, U.K.; Sharma, A. Functional Foods and Their Benefints: An Overview. J. Nutr. Health Food Eng. 2017, 7, 247. [Google Scholar]
  43. Speckmann, B.; Schulz, S.; Hiller, F.; Hesse, D.; Schumacher, F.; Kleuser, B.; Geisel, J.; Obeid, R.; Grune, T.; Kipp, A.P. Selenium increases hepatic DNA methylation and modulates one-carbon metabolism in the liver of mice. J. Nutr. Biochem. 2017, 48, 112–119. [Google Scholar] [CrossRef] [PubMed]
  44. Alehagen, U.; Lindahl, T.L.; Aaseth, J.; Svensson, E.; Johansson, P. Levels of sP-selectin and hs-CRP Decrease with Dietary Intervention with Selenium and Coenzyme Q10 Combined: A Secondary Analysis of a Randomized Clinical Trial. PLoS ONE 2015, 10, e0137680. [Google Scholar] [CrossRef] [PubMed]
  45. Walker, K.A.; Power, M.C.; Hoogeveen, R.C.; Folsom, A.R.; Ballantyne, C.M.; Knopman, D.S.; Windham, B.G.; Selvin, E.; Jack, C.R.; Gottesman, R.F. Midlife Systemic Inflammation, Late-Life White Matter Integrity, and Cerebral Small Vessel Disease: The Atherosclerosis Risk in Communities Study. Stroke 2017, 48, 3196–3202. [Google Scholar] [CrossRef] [PubMed]
  46. Harrington, R.A. Targeting Inflammation in Coronary Artery Disease. N. Engl. J. Med. 2017, 377, 1197–1198. [Google Scholar] [CrossRef] [PubMed]
  47. Würtz, P.; Soininen, P.; Kangas, A.J.; Rönnemaa, T.; Lehtimäki, T.; Kähönen, M.; Viikari, J.S.; Raitakari, O.T.; Ala-Korpela, M. Branched-Chain and Aromatic Amino Acids Are Predictors of Insulin Resistance in Young Adults. Diabetes Care 2013, 36, 648–655. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Tillin, T.; Hughes, A.D.; Wang, Q.; Würtz, P.; Ala-Korpela, M.; Sattar, N.; Forouhi, N.G.; Godsland, I.F.; Eastwood, S.V.; McKeigue, P.M.; et al. Diabetes risk and amino acid profiles: Cross-sectional and prospective analyses of ethnicity, amino acids and diabetes in a South Asian and European cohort from the SABRE (Southall And Brent REvisited) Study. Diabetologia 2015, 58, 968–979. [Google Scholar] [CrossRef]
  49. Zhenyukh, O.; Civantos, E.; Ruiz-Ortega, M.; Sánchez, M.S.; Vázquez, C.; Peiró, C.; Egido, J.; Mas, S. High concentration of branched-chain amino acids promotes oxidative stress, inflammation and migration of human peripheral blood mononuclear cells via mTORC1 activation. Free Radic. Boil. Med. 2017, 104, 165–177. [Google Scholar] [CrossRef] [PubMed]
  50. Johansson, P.; Dahlstrom, O.; Dahlstrom, U.; Alehagen, U. Improved health-related quality of life, and more days out of hospital with supplementation with selenium and coenzyme Q10 combined. Results from a double blind, placebo-controlled prospective study. J. Nutr. Health Aging 2015, 19, 870–877. [Google Scholar] [CrossRef] [PubMed]
  51. Iwase, S.; Kawaguchi, T.; Yotsumoto, D.; Doi, T.; Miyara, K.; Odagiri, H.; Kitamura, K.; Ariyoshi, K.; Miyaji, T.; Ishiki, H.; et al. Efficacy and safety of an amino acid jelly containing coenzyme Q10 and L-carnitine in controlling fatigue in breast cancer patients receiving chemotherapy: A multi-institutional, randomized, exploratory trial (JORTC-CAM01). Support. Care Cancer 2016, 24, 637–646. [Google Scholar] [CrossRef] [PubMed]
  52. Sanoobar, M.; Dehghan, P.; Khalili, M.; Azimi, A.; Seifar, F. Coenzyme Q10 as a treatment for fatigue and depression in multiple sclerosis patients: A double blind randomized clinical trial. Nutr. Neurosci. 2016, 19, 138–143. [Google Scholar] [CrossRef] [PubMed]
  53. Bengtsson, A.A.; Trygg, J.; Wuttge, D.M.; Sturfelt, G.; Theander, E.; Donten, M.; Moritz, T.; Sennbro, C.-J.; Torell, F.; Lood, C.; et al. Metabolic Profiling of Systemic Lupus Erythematosus and Comparison with Primary Sjögren’s Syndrome and Systemic Sclerosis. PLoS ONE 2016, 11, 0159384. [Google Scholar] [CrossRef] [PubMed]
  54. Di Pierro, F.; Rossi, A.; Consensi, A.; Giacomelli, C.; Bazzichi, L. Role for a water-soluble form of CoQ10 in female subjects affected by fibromyalgia. A preliminary study. Clin. Exp. Rheumatol. 2017, 35 (Suppl. 105), 20–27. [Google Scholar]
  55. Zhang, Y.P.; Eber, A.; Yuan, Y.; Yang, Z.; Rodriguez, Y.; Levitt, R.C.; Takacs, P.; Candiotti, K.A. Prophylactic and Antinociceptive Effects of Coenzyme Q10 on Diabetic Neuropathic Pain in a Mouse Model of Type 1 Diabetes. Anesthesiology 2013, 118, 945–954. [Google Scholar] [CrossRef] [Green Version]
  56. Zhang, Y.P.; Mei, S.; Yang, J.; Rodriguez, Y.; Candiotti, K.A. Acute Hypoglycemia Induces Painful Neuropathy and the Treatment of Coenzyme Q10. J. Diabetes Res. 2016, 2016, 1–7. [Google Scholar] [CrossRef] [PubMed]
  57. Jackson, M.I.; Cao, J.; Zeng, H.; Uthus, E.; Combs, G.F. S-Adenosylmethionine-dependent Protein Methylation Is Required for Expression of Selenoprotein P and Gluconeogenic Enzymes in HepG2 Human Hepatocytes. J. Boil. Chem. 2012, 287, 36455–36464. [Google Scholar] [CrossRef] [Green Version]
  58. Jung, M.; Pfeifer, G.P. Aging and DNA methylation. BMC Biol. 2015, 13, 7. [Google Scholar] [CrossRef]
  59. Alehagen, U.; Aaseth, J.; Alexander, J.; Svensson, E.; Johansson, P.; Larsson, A. Less fibrosis in elderly subjects supplemented with selenium and coenzyme Q10-A mechanism behind reduced cardiovascular mortality? BioFactors 2017, 44, 137–147. [Google Scholar] [CrossRef]
  60. Johansson, P.; Dahlström, Ö; Dahlström, U.; Alehagen, U. Effect of selenium and Q10 on the cardiac biomarker NT-proBNP. Scand. Cardiovasc. J. 2013, 47, 281–288. [Google Scholar] [CrossRef]
  61. Jensen-Urstad, K.; Bouvier, F.; Höjer, J.; Ruiz, H.; Hulting, J.; Samad, B.; Thorstrand, C.; Jensen-Urstad, M. Comparison of Different Echocardiographic Methods With Radionuclide Imaging for Measuring Left Ventricular Ejection Fraction During Acute Myocardial Infarction Treated by Thrombolytic Therapy. Am. J. Cardiol. 1998, 81, 538–544. [Google Scholar] [CrossRef]
  62. van Royen, N.; Jaffe, C.C.; Krumholz, H.M.; Johnson, K.M.; Lynch, P.J.; Natale, D.; Atkinson, P.; Deman, P.; Wackers, F.J. Comparison and reproducibility of visual echocardiographic and quantitative radionuclide left ventricular ejection fractions. Am. J. Cardiol. 1996, 77, 843–850. [Google Scholar] [CrossRef]
  63. Surowiec, I.; Johansson, E.; Torell, F.; Idborg, H.; Gunnarsson, I.; Svenungsson, E.; Jakobsson, P.-J.; Trygg, J. Multivariate strategy for the sample selection and integration of multi-batch data in metabolomics. Metabolomics 2017, 13, 114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Surowiec, I.; Vikström, L.; Hector, G.; Johansson, E.; Vikström, C.; Trygg, J. Generalized Subset Designs in Analytical Chemistry. Anal. Chem. 2017, 89, 6491–6497. [Google Scholar] [CrossRef] [PubMed]
  65. A, J.; Trygg, J.; Gullberg, J.; Johansson, A.I.; Jonsson, P.; Antti, H.; Marklund, S.L.; Moritz, T. Extraction and GC/MS analysis of the human blood plasma metabolome. Anal. Chem. 2005, 77, 8086–8094. [Google Scholar] [CrossRef] [PubMed]
  66. Jonsson, P.; Johansson, A.I.; Gullberg, J.; Trygg, J.; A, J.; Grung, B.; Marklund, S.; Sjöström, M.; Antti, H.; Moritz, T. High-Throughput Data Analysis for Detecting and Identifying Differences between Samples in GC/MS-Based Metabolomic Analyses. Anal. Chem. 2005, 77, 5635–5642. [Google Scholar] [CrossRef] [PubMed]
  67. Trygg, J.; Wold, S. PLS regressionon waveletcompressed NIR spectra. Chemom. Intell. Lab. Syst. 1998, 42, 209–220. [Google Scholar] [CrossRef]
  68. Wold, S.; Trygg, J.; Berglund, A.; Antti, H. Some recent developments in PLS modeling. Chemom. Intell. Lab. Syst. 2001, 58, 131–150. [Google Scholar] [CrossRef]
  69. Wold, S. Pattern-Recognition by Means of Disjoint Principal Component Analysis. Pattern Recognit. 1976, 8, 127–139. [Google Scholar] [CrossRef]
  70. Redestig, H.; Fukushima, A.; Stenlund, H.; Moritz, T.; Arita, M.; Saito, K.; Kusano, M. Compensation for Systematic Cross-Contribution Improves Normalization of Mass Spectrometry Based Metabolomics Data. Anal. Chem. 2009, 81, 7974–7980. [Google Scholar] [CrossRef]
  71. Wold, S.; Esbensen, K.; Geladi, P. Principal Component Analysis. Chemom. Intell. Lab. Syst. 1987, 2, 37–52. [Google Scholar] [CrossRef]
  72. Trygg, J.; Wold, S. Orthogonal projections to latent structures (O-PLS). J. Chemom. 2002, 16, 119–128. [Google Scholar] [CrossRef]
  73. Bylesjo, M.; Rantalainen, M.; Cloarec, O.; Nicholson, J.K.; Holmes, E.; Trygg, J. OPLS discriminant analysis: Combiningthe strengths of PLS-DA and SIMCA classification. J. Chemom. 2006, 20, 341–351. [Google Scholar] [CrossRef]
  74. Wiklund, S.; Johansson, E.; Sjöström, L.; Mellerowicz, E.J.; Edlund, U.; Shockcor, J.P.; Gottfries, J.; Moritz, T.; Trygg, J. Visualization of GC/TOF-MS-Based Metabolomics Data for Identification of Biochemically Interesting Compounds Using OPLS Class Models. Anal. Chem. 2008, 80, 115–122. [Google Scholar] [CrossRef] [PubMed]
  75. Efron, B.; Gong, G. A Leisurely Look at the Bootstrap, the Jackknife, and Cross-Validation. Am. Stat. 1983, 37, 36–48. [Google Scholar]
Biomolecules 09 00553 g001 550
Figure 1. OPLS-DA score plot based on each individual’s metabolic profile. (R2X = 0.19, R2Y = 0.57, Q2Y = 0.18). Note: Round circles indicate the men that received the placebo (Pl, green; n = 9), and the squares are men treated with SeQ10 (Ki, blue; n = 9). Note: The metabolite pattern was determined in plasma at the time point of 18 months (T18). Note: One of the individuals receiving placebo took SeQ10 by self-administration.
Figure 1. OPLS-DA score plot based on each individual’s metabolic profile. (R2X = 0.19, R2Y = 0.57, Q2Y = 0.18). Note: Round circles indicate the men that received the placebo (Pl, green; n = 9), and the squares are men treated with SeQ10 (Ki, blue; n = 9). Note: The metabolite pattern was determined in plasma at the time point of 18 months (T18). Note: One of the individuals receiving placebo took SeQ10 by self-administration.
Biomolecules 09 00553 g001
Biomolecules 09 00553 g002 550
Figure 2. Predicted cross-validated Y values (from the OPLS-DA model). Note: Actual values for controls are denoted 0 and for SeQ10-treated are denoted 1, showing that the two groups’ metabolic patterns are significantly different (p = 0.006).
Figure 2. Predicted cross-validated Y values (from the OPLS-DA model). Note: Actual values for controls are denoted 0 and for SeQ10-treated are denoted 1, showing that the two groups’ metabolic patterns are significantly different (p = 0.006).
Biomolecules 09 00553 g002
Biomolecules 09 00553 g003 550
Figure 3. P(corr) loading vector from the OPLS-DA model with the 95 metabolites identified in plasma of men (n = 9) receiving a daily addition of SeQ10 for 18 months as compared with placebo-treated men (controls; n = 8). Note: Bars pointing to the left from mean (0) indicate metabolites with lower levels in SeQ10-treated men. Amino acids are marked in light blue.
Figure 3. P(corr) loading vector from the OPLS-DA model with the 95 metabolites identified in plasma of men (n = 9) receiving a daily addition of SeQ10 for 18 months as compared with placebo-treated men (controls; n = 8). Note: Bars pointing to the left from mean (0) indicate metabolites with lower levels in SeQ10-treated men. Amino acids are marked in light blue.
Biomolecules 09 00553 g003
Table
Table 1. Baseline characteristics of the male study population receiving an intervention of dietary supplementation of selenium and coenzyme Q10 combined over four years.
Table 1. Baseline characteristics of the male study population receiving an intervention of dietary supplementation of selenium and coenzyme Q10 combined over four years.
ActivePlacebop Value
N5544
Age (years)76.1 (2.9)76.3 (2.9)
History
Smokers (present) n (%)5 (7.3)6 (13.6)0.47
Diabetes n (%)11 (20.0)8 (18.2)0.82
Hypertension n (%)34 (61.8)31 (70.5)0.91
IHD n (%)12 (21.8)10 (22.7)0.91
NYHA class I n (%)37 (67.3)22 (50.0)0.08
NYHA class II n (%)12 (21.8)18 (40.9)0.04
NYHA class III n (%)6 (10.9)4 (9.1)0.77
NYHA class IV n (%)00
Unclassified NYHA n (%)00
Medications
ACEI or ARBn (%)11 (20.0)12 (27.3)0.39
Beta-blockers n (%)20 (36.4)17 (38.6)0.82
Digitalis n (%)3 (5.5)0 (0)n/a
Diuretics n (%)13 (23.6)13 (29.5)0.51
Statins n (%)13 (23.6)7 (15.9)0.34
Examinations
EF < 40% n (%)6 (10.9)2 (4.5)0.25
Note: ACEI: ACE- inhibitors; ARB; Angiotension receptor blockers; EF: Ejection fraction; IHD; Ischaemic heart disease; NYHA: New York Heart Association functional class; p values indicate the Chi-square test was used for discrete variables.
Table
Table 2. Metabolites that showed significantly lower levels in plasma of men receiving a daily addition of SeQ10 for 18 months as compared with the controls, the placebo-treated men.
Table 2. Metabolites that showed significantly lower levels in plasma of men receiving a daily addition of SeQ10 for 18 months as compared with the controls, the placebo-treated men.
Metabolitep-Value VIP (OPLS-DA)
Decreased in SeQ10 group
Alanine0.0011.79
Alanine-3-cyano0.0161.62
Arachidonic acidNS1.30
Aspartic acid0.0171.65
Butyric acid, 2-aminoNS1.16
Cysteine0.0481.05
FucoseNS1.21
Glutamic acid0.0221.60
Glycine0.0101.50
Isoleucine0.0311.49
LeucineNS1.45
Methionine0.0301.51
Myo-inositol0.0221.42
Ornithine-1,5-lactam0.0131.29
PhenylalanineNS1.25
Proline-4-hydroxy0.0261.39
Salicylic acidNS1.34
Serine0.0321.44
Threonine0.0491.36
Note: No metabolite showed significantly higher levels in SeQ10 men; Note: VIP: Variables’ Importance to the Projection. Note: T-tests were used in the statistical evaluation. Note: NS: Not significant.
Table
Table 3. Compounds that were found significant according to T-test in at least one batch when SeQ10 and controls at 18 months were compared, with the corresponding p-value and VIP value from this batch.
Table 3. Compounds that were found significant according to T-test in at least one batch when SeQ10 and controls at 18 months were compared, with the corresponding p-value and VIP value from this batch.
Compoundp Value
(95% Confidence Level)		VIP
Decreased in the SeQ10 group
Alanine0.0261.94
Glycine0.0251.95
Isoleucine0.0161.93
Leucine0.0311.60
Lysine0.0251.37
Tryptophan0.0250.95
Increased
1-monohexadecanoylglycerol0.0041.79
Aspartic acid0.0072.19
Hexadecanoic acid0.0131.71
Hexadecenoic acid0.0391.63
Lauric acid0.0081.51
Myristic acid0.0261.44
Oleic acid0.0331.72
Note: p-values obtained from t-tests. Note: VIP values were obtained from the OPLS-DA model.
Table
Table 4. Compounds that were found significant according to t-test in at least one batch when SeQ10 and controls for men at 18 months were compared, with the corresponding p-value and VIP value from this batch.
Table 4. Compounds that were found significant according to t-test in at least one batch when SeQ10 and controls for men at 18 months were compared, with the corresponding p-value and VIP value from this batch.
Compoundp Value
(95% Confidence Level)		VIP
Decreased in the SeQ10 group
Arabinose0.0131.83
Ribose0.0141.76
Sucrose0.0341.83
Xylitol0.0251.41
Increased
Fructose0.0221.96
1,5-anhydroglucitol0.0161.93
Note: p values obtained from t-tests; VIP values were obtained from the OPLS-DA model.
Table
Table 5. Metabolites identified in the main study and validation studies 1 and 2 together with the combined loading vectors from these studies.
Table 5. Metabolites identified in the main study and validation studies 1 and 2 together with the combined loading vectors from these studies.
CompoundMain StudyValidation 1Validation 2
1,5-anhydro-d-glucitolN↑↑
1-dodecanoyl-sn-glycero-3-phosphocholineNN↑
1-MonohexadecanoylglycerolN↑↑
1-Palmitoyl-sn-glycero-3-phosphocholineN↓N
2-aminobutyric acidN↓↓
3-hydroxybutyric acidN↑↑
Alanine↓↓↓
Alanine, 3-cyano↓NN
AllothreonineN↓N
Alpha-ketoglutarate↓NN
Aminomalonic acidNN↓
Arabinose↓↓↓
Arabitol/ribitol↓NN
Arachidonic acid↓NN
Arginine↓↓↓
Asparagine, DL-↓↓↓
Aspartic acid, DL↓↓↑
Benzyl alcohol↓NN
Beta-Alanine↓NN
Butanoic acid, 2,4-dihydroxy-↑NN
Butyric acid, 2-amino↓NN
Campesterol↑↑↑
Cholesterol↑NN
Citric acid↓↑↑
Creatinine↓↓↑
Cysteine↓↓↓
Cystine↓↑↓
Docosanoic acid↑↑N
Docosahexaenoic acid, 4,7,10,13,16,19-(Z,Z,Z,Z,Z,Z)↓N↓
Eicosanoic acid, n-↓NN
Eicosapentaenoic acid↓NN
Ethanolamine↓↑↑
Fructose↑↑↑
Fucose↓NN
Glucose↑NN
Glucose, 1,6-anhydro, beta↓↓N
Glutamic acid↓↓↓
Glutamine↓↓↑
Glyceric acid↓↓↓
Glycerol↑NN
Glycerol-2-phosphate↓NN
Glycerol-3-phosphate↑↓↓
Glycine↓↓↑
Glycolic acidNN↑
Heptadecanoic acid, n-↑↓N
Heptanoic acid↑NN
Hexadecanoic acid↑↑↑
Hexadecenoic acid↑↑↑
Hippuric acid↓NN
Histidine↓↓↓
Indole-3-acetic acid↓NN
Inosine↑NN
Inositol, myo↓↑↓
Inositol-1-phosphate↑↓↑
Isoleucine↓↓↓
Itaconic acid↓NN
Lactic acid,↓↓↑
Lauric acid (dodedecanoic acid)↑↑↑
Leucine↓↓↓
Linoleic acid↑↑↑
Lysine↓↓↓
Malic acid↓↓↑
Malonic acid, 2-amino↑NN
Maltose↓↓↑
MannitolNN↓
Mannose↑NN
Methionine↓NN
Monomethylphosphate↑↑↑
Myristic acid↑↑↑
Nonanoic acid,n↓↑↓
Oleic acid↑↑↑
Ornithine↓↓↓
Ornithine-1,5-lactam↓NN
Oxalic acidNN↓
Phenylalanine↓↓↓
Phosphoric acid↑NN
Pipecolic acid↓NN
Proline↓↓↓
Proline, 4-hydroxy↓NN
Protocatechuic acid↓NN
PseudouridineN↓↑
Putrescine↓NN
Pyridine, 3-hydroxy↓NN
Pyroglutamic acid↓↓↑
Pyruvic acid↑↓↓
Ribose↓↓↓
Salicylic acid↓↓↓
Serine↓↓↓
Sitosterol↑↓N
Sorbitol↑NN
Stearic acid↑↑↑
Succinic acid↑N↑
Sucrose↑↑↑
Taurine↓↓↓
Threitol↑NN
Threonic acid↓↑↓
Threonine↓↓↓
Tocopherol, alpha-↑↑↓
Tocopherol, gamma-↓↓↑
Tryptophan↓↓↓
Tyrosine↓↓↓
Urea↓N↓
Uric acid↓↓↑
Valine↓↓↓
Xylitol↑↓↓
Xylose↓↓↑
Xylulose↓N↑
Note: Arrows pointing down are metabolites with lower levels in the SeQ10-treated men compared with the placebo group, and arrows pointing up are those metabolites that showed higher levels in the SeQ10 treated men. “N” stands for no change as a result of the intervention.

Share and Cite

MDPI and ACS Style

Alehagen, U.; Johansson, P.; Aaseth, J.; Alexander, J.; Surowiec, I.; Lundstedt-Enkel, K.; Lundstedt, T. Significant Changes in Metabolic Profiles after Intervention with Selenium and Coenzyme Q10 in an Elderly Population. Biomolecules 2019, 9, 553. https://doi.org/10.3390/biom9100553

AMA Style

Alehagen U, Johansson P, Aaseth J, Alexander J, Surowiec I, Lundstedt-Enkel K, Lundstedt T. Significant Changes in Metabolic Profiles after Intervention with Selenium and Coenzyme Q10 in an Elderly Population. Biomolecules. 2019; 9(10):553. https://doi.org/10.3390/biom9100553

Chicago/Turabian Style

Alehagen, Urban, Peter Johansson, Jan Aaseth, Jan Alexander, Izabella Surowiec, Katrin Lundstedt-Enkel, and Torbjörn Lundstedt. 2019. "Significant Changes in Metabolic Profiles after Intervention with Selenium and Coenzyme Q10 in an Elderly Population" Biomolecules 9, no. 10: 553. https://doi.org/10.3390/biom9100553

APA Style

Alehagen, U., Johansson, P., Aaseth, J., Alexander, J., Surowiec, I., Lundstedt-Enkel, K., & Lundstedt, T. (2019). Significant Changes in Metabolic Profiles after Intervention with Selenium and Coenzyme Q10 in an Elderly Population. Biomolecules, 9(10), 553. https://doi.org/10.3390/biom9100553

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