Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Strong Associations Exist among Oxidative Stress and Antioxidant Biomarkers in the Circulating, Cellular and Urinary Anatomical Compartments in Guatemalan Children from the Western Highlands

  • María J. Soto-Méndez ,

    clinicanutricion@gmail.com

    Affiliation Center for the Studies of Sensory Impairment, Aging, and Metabolism–CeSSIAM–Guatemala City, Guatemala

  • Concepción M. Aguilera,

    Affiliations Department of Biochemistry and Molecular Biology II, Institute of Nutrition and Food Technology “José Mataix”, Center for Biomedical Research, University of Granada, Granada, Spain, Networking Biomedical Research for Obesity and Nutrition–CIBERobn-, Madrid, Spain

  • María D. Mesa,

    Affiliations Department of Biochemistry and Molecular Biology II, Institute of Nutrition and Food Technology “José Mataix”, Center for Biomedical Research, University of Granada, Granada, Spain, Thematic Networks of Cooperative Research–RETIC–, Carlos III Health Institute–ISCIII–, General Sub-Directorate for Research Assessment and Promotion and the European Regional Development Fund–ERDF–ref. RD12/0026, Madrid, Spain

  • Laura Campaña-Martín,

    Affiliation Department of Biochemistry and Molecular Biology II, Institute of Nutrition and Food Technology “José Mataix”, Center for Biomedical Research, University of Granada, Granada, Spain

  • Victoria Martín-Laguna,

    Affiliation Department of Biochemistry and Molecular Biology II, Institute of Nutrition and Food Technology “José Mataix”, Center for Biomedical Research, University of Granada, Granada, Spain

  • Noel W. Solomons,

    Affiliation Center for the Studies of Sensory Impairment, Aging, and Metabolism–CeSSIAM–Guatemala City, Guatemala

  • Klaus Schümann,

    Affiliation Molecular Nutrition Unit, ZIEL, Research Center for Nutrition and Food Science, Technische Universität München, Freising, Germany

  • Ángel Gil

    Affiliations Department of Biochemistry and Molecular Biology II, Institute of Nutrition and Food Technology “José Mataix”, Center for Biomedical Research, University of Granada, Granada, Spain, Networking Biomedical Research for Obesity and Nutrition–CIBERobn-, Madrid, Spain

Correction

12 Feb 2016: Soto-Méndez MJ, Aguilera CM, Mesa MD, Campaña-Martín L, Martín-Laguna V, et al. (2016) Correction: Strong Associations Exist among Oxidative Stress and Antioxidant Biomarkers in the Circulating, Cellular and Urinary Anatomical Compartments in Guatemalan Children from the Western Highlands. PLOS ONE 11(2): e0149740. https://doi.org/10.1371/journal.pone.0149740 View correction

Abstract

Background

A series of antioxidant enzymes and non-enzymatic compounds act to protect cells from uncontrolled propagation of free radicals. It is poorly understood, though, to what extent and how their interaction is harmonized.

Objectives

To explore associative interactions among a battery of urinary and blood biomarkers of oxidative stress and enzymatic and non-enzymatic markers of the antioxidant defense system in children from low income households.

Methods

For this cross-sectional descriptive study, urine, red cells, and plasma were sampled in 82 preschool children attending three daycare centers in Quetzaltenango Guatemala. The urinary oxidative stress biomarkers studied were F2-isoprostanes and 8-hydroxy-deoxy-guanosine. Red cell enzyme activities measured were: catalase, superoxide dismutase, glutathione peroxidase and glutathione reductase. Circulating non-enzymatic antioxidants selected were: retinol, tocopherols, β-carotene and coenzymes Q9 and Q10.

Results

In a Spearman rank-order correlation hemi-matrix, of 55 paired combinations of the 11 biomarkers, 28 (51%) were significantly correlated among each other (p≤ 0.05), with the strongest association being retinol and tocopherols (r = 0.697, p<0.001), and 4 associations (9%) showed a trend (p> 0.5 to ≤ 0.10). F2-isoprostanes showed the greatest number of cross-associations, having significant interactions with 8 of the 10 remaining biomarkers. Goodness-of-fit modeling improved or maintained the r value for 24 of the significant interactions and for one of the 5 borderline associations. Multiple regression backward stepwise analysis indicated that plasma retinol, β-carotene and coenzyme Q10 were independent predictors of urinary F2-isoprostanes.

Conclusion

Numerous significant associations resulted among biomarkers of oxidation and responders to oxidation. Interesting findings were the apparent patterns of harmonious interactions among the elements of the oxidation-antioxidation systems in this population.

Introduction

The oxidation of substrates in the mitochondria to generate ATP is associated with continuous free radical formation [1,2]. At the same time, free radicals can damage cellular and organelle membranes, the cell nucleus, and protein chains. Free radical action forms isoprostanes from lipids [36], 8-hydroxy-deoxyguanosine from DNA [68], and carbonyl-amino acids from peptide chains [9]. With respect to the origins of the marker of DNA damage, the formation seems to arise in the repair process for the damaged nucleoside of the nucleotide chain [10]. Action of redox enzymes is part of the normal cellular metabolism in living organisms [11,12]. Thus, oxidative mechanisms are used to destroy microbes within macrophages [13,14]. Renal excretion of such split products from lipids and nucleic acids converts them into urinary biomarkers of systemic oxidative stress that help to assess the protective reaction/response by antioxidant mechanisms [15,16].

The antioxidation system in the human organism is complex and geared to the suppression of uncontrolled oxidation while permitting essential and beneficial oxidation reactions [17]. A series of compounds, including dietary nutrients, have free radical quenching properties or work in conjunction with enzymatic antioxidation reactions or both [18]. Retinol and its derivatives, ascorbate, and tocopherols, as well as diverse carotenoids and co-enzymes Q [19], are among the antioxidant nutrients. Their circulating concentrations roughly reflect the corresponding nutritional status of an individual.

Finally, a series of enzymes protect the cell from oxidation reactions involving molecular or atomic oxygen [20]. These include superoxide dismutase (SOD), catalase (CAT), glutathione peroxidase (GPX) and glutathione reductase (GSR), which are the primary intercellular antioxidant enzymes that detoxify oxygen-containing free radicals produced during normal aerobic respiration [2123]. The zinc and copper-dependent SOD transforms the O2-.into H2O2 and O2 [21], whereas CAT and GSR catalyze the reduction of H2O2 into H2O [22,24]. GPX works in concert with GSR in a system that neutralizes hydrogen peroxide [23].

In a research project entitled “Study on the normative state and inter- and intra-individual variation in growth, hematology, hydration, and markers of oxidation, infection and inflammation in pre-school children with a similar dietary intake”, we collected concurrent data in a series of two urinary biomarkers of oxidative stress, four enzymatic biomarkers of the antioxidant defense system (ADS) in red blood cells (rbc), and five circulating antioxidant nutrients from preschool children in a governmental system of daycare centers. The degree, to which the interactions of biomarkers reflect biological harmonization is the central query of our companion article on the interaction of immunological biomarkers in this population sample [25], which is included as supplemental material (S1 Article). We use parallel analytical approaches of cross-associations with linear and non-linear correlations and backward elimination multiple-regressions to quantify the magnitude of interactions. We present here the result of an exploration of the mutual interplay and interactions of biomarkers of oxidative stress and the ADS in the context of variable environmental and genetic situations within the stabilizing influence of a common institutional dietary offering.

Materials and Methods

Study Design

This study was a descriptive, cross-sectional field research on the variation and association among oxidative stress and antioxidant biomarkers.

Setting and Subjects

The setting for this research project was the Western Highlands of Quetzaltenango, located 136 miles from Guatemala City at 8717 ft above the sea level [26]. Attendees at three daycare centers’ from the SOSEP (Secretaría de Obras Sociales de la Esposa del Presidente {Secretariat of the Beneficial Works of the First Lady}) system were assigned to participate in this study. Center A in La Esperanza was the semi-urban setting, located 2 miles away from downtown Quetzaltenango; Center B in La Puerta del Llano was located in a marginal-urban setting in the outskirts of the city of Quetzaltenango; and Center C was a rural setting situated 15 miles away from Quetzaltenango in La Estancia, San Martín Chile Verde. Almost all children attending the centers were of Maya indigenous ethnicity. However, some of the living habits, pastimes and physical characteristics varied between centers because of the infrastructure and the environment of each location.

Recruitment and Enrolment

Inclusion Criteria.

To be enrolled in the study children had to be attending one of the selected SOSEP Centers, be aged 2 to 7 years, and maintaining an attendance of at least 80% during the 8 weeks of the fieldwork. Furthermore, subjects had to be apparently healthy and with no dietary restrictions related to the foods offered by the 8-week rotating menu from the SOSEP system.

Exclusion Criteria.

Children who did not adhere to the urine collection routine, who refused to participate in the study or whose parents did not sign the consent form were excluded.

Ethical Considerations.

The study protocol was approved by The Human Subjects Committee of the Center for the Studies of Sensory Impairment, Aging and Metabolism (CeSSIAM) and was registered at clinicaltrials.gov as NCT02203890. A parent or guardian signed the written consent form. The SOSEP’s director for the Quetzaltenango area had authorized the study. The diet offered to the children was complemented when required, in order to provide all food items on the menu. STROBE statement for this article is included as a supporting information file (S1 STROBE Checklist)

Collection, Handling and Storage of Biological Samples

We collected 24-h urine samples and a 5 mL blood sample during the last of the 8 study weeks. We started urine collection at each daycare center when a child arrived (between 7:00 and 8:00 a.m.) with SOSEP personnel assisting for the collection, using BD Vacutainer® No.364999 plastic 24-h collection container (Becton-Dickinson, New Jersey, USA). After training, parents continued the collection at home. Urine collection was finished at the center, 24 h after initiation. The collection process was repeated if there was a suspicion that it was incomplete.

The collection was taken to the laboratory, where the entire urine sample was agitated in order to obtain homogeneity before aliquoting. One aliquot was stored at -80°C for 37 to 46 weeks before being sent to Granada, Spain, on dry-ice for measurement of oxidative biomarkers. The other aliquot was stored at the same temperature in Guatemala for other measurements that were held in the country.

An experienced phlebotomist collected blood samples using BD Vacutainer® 4 mL tubes with EDTA (No.367861) and Safety-LokTM deposable needles (No.367281). Samples were centrifuged to separate red blood cells from plasma; both were stored in Nalgene® Cryogenic Vials (No.5000-0012) at -80°C until shipment to the Institute of Nutrition and Food Technology, Center of Biomedical Research, University of Granada, Granada, Spain in order to determine antioxidant enzymes activity in red blood cells and antioxidant nutrients in plasma.

Laboratory Assays and Analyses

Plasma retinol, tocopherols, ß-carotene, and coenzymes Q9 and Q10 (Co-Q9 and Co-Q10) were determined with high-performance liquid chromatographic (HPLC) methods using a 100-microliter aliquot of plasma sample previously deproteinized with 1-propanol [27] (CV intra-day 6.1%, CV inter-day 3.8%, LOD 0.012 mg/mL for retinol; CV intra-day 2.7%, CV inter-day 5.6%, LOD 0.062 mg/mL for tocopherols; CV intra-day 6.7%, CV inter-day 13.0%, LOD 0.012 mg/mL for β-carotene; CV intra-day 1.5%, CV inter-day 4.5%, LOD 0.012 mg/mL for Co-Q9; and CV intra-day 3.7%, CV inter-day 9.5%, LOD 0.012 mg/mL for Co-Q10). Plasma levels of retinol, tocopherols, ß-carotene, Co-Q9 and CoQ10 were assayed by high pressure liquid chromatography coupled to mass spectrometry (HPLC-MS), using methanol 0.1% and isocratic formic acid as solvent, with a flow of 0.5 ml/min in a ACQUITY UPLCr BEH C18 50 mm column (internal diameter 2.1 mm and particle size 1.7μm). β-carotene was also determined after extraction with 1-propanol in a HPLC system attached to a multiwavelength ultraviolet detector set at 450 nm. All these compounds were identified by predetermining the retention times of individual standards [28].

Hemoglobin concentration was determined by use of Drabkin´s reagent [29] (Sigma D5941) (CV intra-day 2.4%; CV inter-day 3.1%). The final concentration of Hb was adjusted to 5 mg/mL for antioxidant enzymes analyses. Catalase (CAT) activity was determined using de Aebi [30] (CV intra-day 15.5%, CV inter-day 17.4%) and it is is expressed as nmol/ (L · g Hb). SOD activity was assayed according to the methods of McCord & Fridovich (McCord & Fridovich 1969 [31]), using xanthine and xanthine oxidase to generate superoxide radicals. These radicals oxidize the cytochrome c, generating color measured at 450 nm. The presence of SOD competes with cytochrome c and blocks color generation measured at 450 nm (CV intra-day 5.1%, CV inter-day 7.0%). Data are expressed as U/mg Hb. GR was determined by measuring the rate of reduced nicotinamide adenine dinucleotide phosphate (NADPH) oxidation in the presence of oxidized glutathione (GSSG) according to Carlberg and Mannervik [32] (CV intra-day 8.0%, CV inter-day 10.3%). The results are expressed as U/g Hb. Finally, GPX was analyzed at 340 mn using the procedure developed by Flohé and Günzler [33] with tert-butyl hydroperoxide as substrate (CV intra-day 7.8%, CV inter-day 12.9%). The results are expressed as U/g Hb.

Urinary biomarkers 15-isoprostane F2t (F2-Iso) and 8-Hydroxydeoxyguanosine (8-OHdG) were determined using ELISA assay kits (Oxford Biomedical Research, Inc., Catalog # EA84.102606, Michigan, USA and JaICA, Nikken SEIL Co., Ltd, Catalog# IM-KOGHS 040914E, Shizuoka, Japan, respectively).

Data Handling and Statistical Analyses

Data were organized and recorded in an SPSS version 20.0 database in order to run all the statistical analyses. Normality of variables was assessed using the Kolmogrov-Smirnov test. Descriptive statistics are presented as median, 95% CI and range. We ran the Spearman correlation coefficients according to the distribution of the sample. When the Spearman test gave a significant correlation coefficient (p <0.05) or one with a tendency to be significant (p >0.05 to <0.10), we ran goodness-of-fit models of SPSS in order to detect any improvement in the relation. Multiple regression backward stepwise analyses were performed to develop models to predict values of the urinary oxidative biomarkers F2-Iso and 8-OHdG from the antioxidant defense system parameters measured in the present study (SOD, CAT, GPX and GSR, and tocopherols, retinol, β-carotene and coenzymes Q9 and Q10) using the same software. Durbin-Watson statistics were used to assess whether the assumption of independent errors for the variables was tenable. ANOVA testing was done to determine whether the selected model was significantly better at predicting the outcome than using the mean as a “best guess”. The F represents the ratio of the improvement in the prediction as a result of fitting the model relative to the inaccuracy that still exists in the model. Variance inflation factor (VIF) and tolerance statistics were obtained to assess whether there was some co-linearity among the independent variables [34].

Results

Characteristics of the Participants

Of the 87 children enrolled in the study as a whole, binary samples were available variously in from 78 to 82 cases. In this binary sub-sample, 38 participants were girls, and 44 were boys. The ages ranged from 23 to 81 mo, with a mean of 55 ± 16 mo, and a median of 56 mo. Fig 1 illustrates the sex- and age-distribution of the children for each of the three day-care centers, and provides data on the median ages by sex.

thumbnail
Fig 1. Characteristics of the Subjects, by setting and sex.

https://doi.org/10.1371/journal.pone.0146921.g001

Descriptive Statistics of the Biomarkers Concentration

Table 1 shows the median, 95% confidence interval, and minimum and maximum values for each of the 11 biomarkers, as well as normalized values for the urinary oxidative stress biomarkers, adjusted for urinary creatinine.

thumbnail
Table 1. Descriptive Statistics of the Biomarkers Concentration.

https://doi.org/10.1371/journal.pone.0146921.t001

Geographic Selectivity of the Biomarker concentration

Shown in Table 2 are the median values and 95% CI for the 13 biomarker expressions used in the study compared across the three geographical sites of the study. As compared by a three-way Kruskal-Wallis test, there were significant differences by site for all biomarkers, with the superscript “a” and the upward arrow indicating the highest in the row, and the superscript “b” and the downward arrow indicating the lowest. By a priori logic, the higher levels of urinary oxidative stress biomarkers and the lower concentrations of anti-oxidant vitamins were considered unfavorable, and presented in bold font. There is no a priori assignment of anti-oxidant enzyme activities. As seen by the accumulation of bold font, the marginal-urban center “B” was unfavorable for all nine relevant biomarkers. By contrast, the semi-urban center “A” had the favorable direction in eight of the nine indicators.

thumbnail
Table 2. Comparisons of oxidation and anti-oxidation biomarkers concentrations and activities by setting.

https://doi.org/10.1371/journal.pone.0146921.t002

The upward arrow represents the highest value(s) in the cross-row set, whereas the downward arrow represents the lowest value(s). The hyphen represents the intermediary value in a set. Such analysis was not made for antioxidant enzymes, as they do not have a reference value.

The bold font indicates the least favorable median value. Anti-oxidant enzymes have no a priori assignment for favorability.

Hemimatrix of Linear Spearman Correlations for Inter-Biomarker Associations

Fig 2 shows the Spearman correlation coefficients for 55 binary inter-biomarker associations among the 11 variables including 2 urinary oxidative stress biomarkers, 4 rbc antioxidant enzymes and 5 plasma antioxidant compounds. Of these, 28 (51%) met the 5% statistical significance criterion. The strongest significance was found for the association between retinol and tocopherols (r = 0.697, p<0.001) and the weakest between SOD and GSR (r = 0.220, p = 0.048).

thumbnail
Fig 2. Spearman Correlation Coefficient Hemi-Matrix for Inter-Relationships of Biomarkers.

The Spearman rank-order correlation coefficients Hemi-Matrix for mutual, cross inter-relationships of the 11 measured biomarkers are illustrated within the 66 cells of pertinent reference. The dark-shaded cells represent the 11 auto-correlations. The remaining 55 cells illustrate the probability level for the corresponding Spearman r value. The 28 medium-shaded cells have statistically-significant associations, whereas the 27 clear cells have non-significant associations. Corresponding units of concentration or activity are given with the biomarkers along the horizontal axis; n values are given in the vertical axis.

https://doi.org/10.1371/journal.pone.0146921.g002

Comparison of Goodness-of-Fit Correlations for Inter-Biomarker Associations

Table 3 lists the 28 significant binary associations as analyzed by Spearman rank-order correlation and the corresponding goodness-of-fit correlation coefficient with the appropriate curve form. In 19 cases of the series (70%), the r value was improved or remained the same, whereas in 8 cases (30%) the goodness-of-fit r-value was lower than provided by Spearman analysis. The associations showing the greatest improvements for the strength of correlation were 8-OHdG adjusted for Creatinine and retinol (+75%) and SOD and β-carotene (+75%).

thumbnail
Table 3. Comparison of Spearman and Non-Linear Correlation Coefficients in Inter-Biomarker Significant Associations.

https://doi.org/10.1371/journal.pone.0146921.t003

There were 5 correlation coefficients in the probability range of 0.051 to 0.100 including those between 8-OHdG and SOD, CAT and retinol, GSR and retinol, β-carotene and retinol, and 8-OHdG and GSR. Only the latter two associations in this series achieved significance at p<0.05 when applying the goodness-of-fit model.

Multiple Regression Models

Table 4 shows the coefficients of the multiple regression models for the dependent variables urinary F2-isoprostanes and 8-hydroxy-deoxy-guanosine adjusted for creatinine using as independent variables the antioxidant defense system parameters measured in the present study (red blood cells activities of SOD, CAT, GPX and GSR, and plasma concentrations of tocopherols, retinol, β-carotene and coenzyme Q9 and Q10). Plasma retinol and β-carotene were independent predictors of F2-Iso, accounting for about 25% of the total variability, whereas plasma retinol and Co-Q10 were independent predictors of 8-OHdG accounting for 22% of the variability.

thumbnail
Table 4. Coefficients of the multiple regression model for the dependent variables of urinary F2-Iso and 8-OHdG, adjusted for creatinine, using as independent variables the antioxidant defense system parameters measured in the present study (SOD, CAT GPX and GSR, and plasma concentrations of tocopherols, retinol, β-carotene and Co-Q9 and Co-Q10).

https://doi.org/10.1371/journal.pone.0146921.t004

Fig 3 illustrates 6 selected goodness-of-fit association curves showing the different curve-forms that appeared with the transformed regression model. These include three with cubic form and one each with exponential, sigmoid and linear forms as the result of the goodness-of-fit adaptation. Overall, the cubic configuration dominated with 14 (52%) of all transformations performed.

thumbnail
Fig 3. Selected associations between variables with the superimposition of the Goodness-of-Fit curve form and the correlation coefficient “r” for the model’s equation.

https://doi.org/10.1371/journal.pone.0146921.g003

Discussion

The underlying research hypothesis of this field survey was that the status of oxidative stress and other diagnostic domains would show substantial heterogeneity across the sample, despite the narrow age range and a common institutional menu. The findings amply confirm the notion of wide, inter-individual variance. At the second level, however, we observe broad, multiple and harmonious associations among antioxidant micronutrients, antioxidant enzymes and resulting free radical damage to lipid and nuclear cells constituents. Although one can never impute causality nor confirm direction of causality merely from associations, this opens up the suggestion that the control of low-grade oxidative stress, indeed, behaves like an interactive system within the organism. It has been suggested [35] that a number of adverse environmental exposures or chronic nutrient deficiencies or illnesses, common to the unsanitary and impoverished conditions of low-income settings, can induce or exacerbate cellular oxidation. Of special interest is the existence of indoor smoke contamination of the Western Highlands homes, which can provoke oxidative stress [36]. It was in fact, the scheme of common generalized environmental stressors that framed the rationale for our previous publication [25] examining the interaction of inflammatory biomarkers in these same subjects (S1 Article).

Given that the subjects came from residential settings as described above, with high levels of poverty and low levels of hygiene, it would not be surprising to find values of the biomarkers mildly outside of the normative ranges. For the urinary biomarkers of oxidative stress there are no established normal values for children. Using the normative range for F2-Iso, adjusted for creatinine, given as 3.0 ng/mg of creatinine by the manufacturers in the kit (Source: Oxford Medical Research manual), 66 subjects (85%) are above the upper value, when similarly adjusted. Also, all our subjects had higher values than the mean of 0.7 ng/mg of creatinine based on 39 healthy adolescent controls as published earlier, when compared in common units [37]. When comparing our results to normative values 8.4 ng/mg of creatinine for 8-OHdG provided with the used kit (Source: JaICA manual), 70 of our subjects (90%) had values above the adult mean value adjusted for creatinine concentration. Moreover, when compared to the median value of 9.5 ng/mg of creatinine for 14 healthy children aged 2 to 15 y of age [38], 68 of our subjects (87%) were above their median creatinine-adjusted value. Thus, it seems that a mild elevation in both urinary markers of oxidation is the rule in our population.

Unlike the urinary and plasma biomarkers in this study, for which normative ranges can be offered, activity expressions for anti-oxidant enzymes in red cell membrane cannot be treated or interpreted in a parallel manner. The large variation in assays and their expression across the literature precludes assigning diagnostic reference standard criteria, and they can best be interpreted in relative terms, as we have done in the present analysis of findings in Tables 1 and 2.

We measured circulating concentrations of five antioxidant nutrients. With respect to the reference system [39], the normal range of plasma circulating retinol is 30 to 95 μg/dL; in this context 75 of 81 samples (92%) would be considered low. This finding is surprising, as the most recent national nutrition survey [40] found the 6 to 59 mo population of Guatemala to be free of hypovitaminosis A; we suspect that there might be a calibration error in our analyses, but one that would not invalidate the relative rank-order among subjects in any way. The same reference source [39], provides the normative range from 5 to 20 mg/L for circulating levels of tocopherols; in our sample 64 of 81 samples (79%) are considered below the normal range. Variation in higher or lower levels of circulating lipids, which are the carriers of tocopherols, can distort the interpretation of vitamin E status [41,42]. Unfortunately, we do not have measurement of plasma cholesterol or triglycerides for this population. Nonetheless, at best, only about half of the low levels would normalize by adjustment [41], and even 40% prevalence of low tocopherols is an ominous finding.

The circulating levels of carotenoids are so dependent on the consumption of foods of plant origin for which no international references can or should be made. A previous longitudinal study of urban schoolchildren in Guatemala, conducted between February and May, documented a 3-fold increase in total carotene concentration as compared to other seasons, attributed to the mango season [43]. In the present study, blood was sampled from August through October, i.e. during the downslope of the annual ß-carotene cycle from native mangos.

According to Molyneux et al. [44], the normal range for Co-Q10 is 0.41 to 1.55 mg/L. With our values ranging from 0.15 to 0.32 mg/L, all subjects in this study would be considered to have low Co-Q10 status. However, with the reference range for Co-Q9 from 0.007 to 0.037 mg/L [44] and with our subjects’ distribution 0.040 to 0.060 mg/L, our lowest value is above the normative reference range.

After the description of biomarker values in the context of status assessment in Table 1, we examined the status across geographic sites. In terms of the hypothesis of mutual associations, the clustering of favorable and adverse diagnostic implications of the biomarkers, as shown in Table 2, adds another dimension to the confirmation. In this respect, the association is geographic, and those with the greatest evidence of oxidative stress as a group also have the lowest concentrations of the vitamins and vitamin-like constituents in their circulation.

Multiple comparisons to evaluate interactions always run the risk of spuriously significant correlations. In our case, with a 5% probability criterion, we would have expected to find 2 and 3 correlations with a p value of <0.05 by chance alone. With our finding of a total of 28, at least 25 interactions are likely to represent truly significant associations. In addition, the fact that the goodness-of-fit models are also significant and improved in their r value in 86% of the instances is additional evidence for the validity of the Spearman associations encountered to be truly significant, and not spurious.

Although causality cannot be deduced from association alone, assigning directionality in the regression model (Table 3) is illustrative. Of all of the biomarkers, adjusted F2-Iso had the widest range of cross associations, with 7 of the 10 companion biomarkers; only adjusted 8-OHdG, GPX and Co-Q10 did not associate with F2-Iso adjusted for creatinine. Plasma retinol and β-carotene were identified as independent predictors of adjusted F2-Iso after running the backwards multiple regression model. Tocopherols also showed 7 significant cross-associations. The multiple-regression model was also applied for adjusted 8-OHdG; and SOD activity, plasma retinol and Co-Q10 were the resultant independent predictors.

Clearly, excessive oxidation is suggested to be associated with an increase in the urinary biomarkers of oxidative stress F2-Iso and 8-OHdG. Interpretation of the other biomarkers is not so straightforward. In theory, the response to oxidative stress could deplete a vitamin or affect the activities of enzyme biomarkers or produce a compensatory and protective increase. In the present experience, urinary excretion of oxidative biomarkers associates directly with the activity of antioxidant enzymes and inversely with the vitamin concentration. This association suggests that the enzymatic response is compensatory, whereas the antioxidant nutrients may be being consumed by the oxidation process. This interpretation is further confirmed by the fact the relationships than exist between ADS enzymes and non-enzymatic antioxidants are universally inverse.

Strengths and limitations of the study

We acknowledge certain strengths and weaknesses of design and execution of the study. The principal strength is that it deals with childhood, with children being of preschool age. In theory, the settings with similar dietary offering should control one important variable, and narrow overall variance. We were ambitious enough to attempt 24-h urine collection in the design for young children, with relative success in obtaining complete collections. We recognize, however, that we did not obtain full adherence in this effort due to the young age of the subjects; thus reliance on concentration-related variables for urinary analytes had to supersede that on total one-day output. An additional strength is the possibility to relate a large variety of oxidative and antioxidant biomarkers of different classes and collected from different body fluids to each other, although, we have not considered or measured a variety of other endogenous and exogenous small-molecule antioxidants such as oxidized and reduced glutathione, and alpha-lipoic acid, ferritin, uric acid, bilirubin, metallothioneine, L-carnitine and melatonin [45]. In terms of limitations, however, our central thesis of interactions is based on geographical clustering and rank-order correlations and we remain aware that a conclusion for causality can never be derived from associations alone. The use of an ELISA method for urinary 8-OHdG provides a reasonable approach for comparisons made with the same kit, but we realize that the external validity is inferior to methods based on chromatography [46]. Finally, there are some inconsistences between the retinol status observed here and other experiences seen in Guatemala, and the lack of adjustment of tocopherols levels for circulating lipids may be producing an overestimation of deficiency for this nutrient.

Conclusion

In conclusion, the largely common menu offering at the daycare centers during 5 days per week was no guarantee of the nutritional or a physiological uniformity across a sample of low-income preschoolers. Nutritional biomarkers varied widely, with many falling below the cutoff for adequate status. The best available explanation for mildly elevated urinary biomarkers of oxidative stress would be exposure to the constellation of household and community exposures originating in the poor sanitation, hygiene, and air and water quality of the areas of study. Most strikingly, however, are apparent patterns of harmonious interactions among the elements of the oxidation-antioxidation systems within a setting that presents a variety of mild to moderate oxidative stress from the conditions of daily existence.

Supporting Information

S1 STROBE Checklist. STROBE statement for observational studies.

https://doi.org/10.1371/journal.pone.0146921.s001

(DOC)

S1 Article. Article by Soto-Méndez, et al. Associations among Inflammatory Biomarkers in the Circulating, Plasmatic, Salivary and Intraluminal Anatomical Compartments in Apparently Healthy Preschool Children from the Western Highlands of Guatemala.

PLoS ONE DOI:10.1371/journal.pone.0129158 June 15, 2015.

https://doi.org/10.1371/journal.pone.0146921.s002

(PDF)

Acknowledgments

We thank all the persons involved in this project. Raquel Campos, Sheny Romero, Carlos Tánchez, Jeniece Alvey, Victoria Pérez Lima and Pablo Cifuentes assisted on the sample collection alone with SOSEP personnel. We acknowledge the kind assistance of D’Ann Finley in reviewing the use of language in the crafting of the revised manuscript. Mainly, we thank the children who participated and their parents for authorizing their participation and assisting with sample collection. This work will be used in partial fulfillment of the doctoral requirements that will allow the graduate student, Maria Jose Soto-Méndez, to obtain a Ph.D. degree at the University of Granada, Spain.

Author Contributions

Conceived and designed the experiments: NWS KS AG. Performed the experiments: CMA MDM LCM VML. Analyzed the data: MJSM. Wrote the paper: NWS MJSM KS AG. Planned the fieldwork: NWS MJSM.

References

  1. 1. Sinha K, Das J, Pal P, Sil P. Oxidative stress: the mitochondria-dependent and mitochondria-independent pathways of apoptosis. Arch Toxicol 2013;87:1157–80. pmid:23543009
  2. 2. Zorov D, Juhaszova M, Sollott S. Mitochondrial reactive oxygen species (ROS) and ROS-induced ROS release. Physiol Rev 2014;94:909–50. pmid:24987008
  3. 3. Morrow JD, Hill KE, Burk RF, Nammour TM, Badr KF, Roberts LJ. A series of prostaglandin F2-like compounds are produced in vivo in humans by a non-cyclooxygenase, free radical-catalyzed mechanism. Proc Natl Acad Sci U S A 1990;87:9383–7. pmid:2123555
  4. 4. Pratico D. F2-isoprostanes: sensitive and specific non-invasive indices of lipid peroxidation in vivo. Atherosclerosis 1999;147:1–10. pmid:10525118
  5. 5. Roberts LJ, Milne GL. Isoprostanes. J Lipid Res 2009;50 Suppl:S219–23. pmid:18957694
  6. 6. Il’yasova D, Scarbrough P, Spasojevic I. Urinary Biomarkers of Oxidative Status. Clin Chim Acta 2012;413:1446–53. pmid:22683781
  7. 7. Lunec J, Holloway K, Cooke M, Faux S, Griffiths H, Evans M. Urinary 8-OXO-2’-deoxyguanosine: Redox regulation of DNA repair in vivo? Free Radic Biol Med 2002;33:875–85. pmid:12361799
  8. 8. Prabhulkar S, Li C-Z. Assessment of oxidative DNA damage and repair at single cellular level via real-time monitoring of 8-OHdG biomarker. Biosens Bioelectron 2010;26:1743–9. pmid:20863679
  9. 9. Stadtman ER, Levine RL. Free radical-mediated oxidation of free amino acids and amino acid residues in proteins. Amino Acids 2003;25:207–18. pmid:14661084
  10. 10. Evans MD, Saparbaev M, Cooke MS. DNA repair and the origins of urinary oxidized 2'-deoxyribonucleoside. Mutagenesis. 2010 Sep;25(5):433–42. pmid:20522520
  11. 11. Andersen HR, Nielsen JB, Nielsen F, Grandjean P. Antioxidative enzyme activities in human erythrocytes. Clin Chem 1997;43:562–8. pmid:9105255
  12. 12. Léger C, Bertrand P. Direct Electrochemistry of Redox Enzymes as a Tool for Mechanistic Studies. Chem Rev 2008;108:2379–438. pmid:18620368
  13. 13. Watanabe K, Shibuya S, Ozawa Y, Nojiri H, Izuo N, Yokote K, et al. Superoxide dismutase 1 loss disturbs intracellular redox signaling, resulting in global age-related pathological changes. Biomed Res Int 2014;2014:140165. pmid:25276767
  14. 14. Clifford D, Repine J. Hydrogen peroxide mediated killing of bacteria. Mol Cell Biochem 1982;49:143–9. pmid:6298593
  15. 15. Milne G, Dai Q, Roberts L Jr. The isoprostanes-25 years later. Biochim Biophys Acta 2014;pii:S1388–981.
  16. 16. Loft S, Fischer-Nielsen A, Jeding I, Vistisen K, Poulsen H. 8-Hydroxydeoxyguanosine as a urinary biomarker of oxidative DNA damage. J Toxicol Enviromental Heal 1993;40:391–404.
  17. 17. Halliwell B. How to characterize a biological antioxidant. Free Radic Res Commun 1990;9:1–32. pmid:2159941
  18. 18. Monsen E. Dietary reference intakes for the antioxidant nutrients: vitamin C, vitamin E, selenium, and carotenoids. J Am Diet Assoc 2000;100:637–40. pmid:10863565
  19. 19. Zempleni J, Rucker R, McCormick D, Suttie J. Handbook of vitamins. vol. 11. 4th ed. Boca Raton: CRC Press; 1991. https://doi.org/10.1016/0891-5849(91)90176-4
  20. 20. Valko M, Leibfritz D, Moncol J, Cronin MTD, Mazur M, Telser J. Free radicals and antioxidants in normal physiological functions and human disease. Int J Biochem Cell Biol 2007;39:44–84. pmid:16978905
  21. 21. McCord J, Fridovich I. Superoxide dismutase: the first twenty years (1968–1988). Free Radic Biol Med 1988;5:363–9. pmid:2855736
  22. 22. Chelikani P, Fita I, Loewen PC. Diversity of structures and properties among catalases. Cell Mol Life Sci 2004;61:192–208. pmid:14745498
  23. 23. Deponte M. Glutathione catalysis and the reaction mechanisms of glutathione-dependent enzymes. Biochim Biophys Acta 2013;1830:3217–66. pmid:23036594
  24. 24. Łukaszewicz-Hussain A, Moniuszko-Jakoniuk J. Liver Catalase, Glutathione Peroxidase and Reductase Activity, Reduced Glutathione and Hydrogen Peroxide Levels in Acute Intoxication with Chlorfenvinphos, an Organophosphate Insecticide. Pol J Env Stud 2004;13:303–9.
  25. 25. Soto-Méndez MJ, Romero-Abal ME, Aguilera CM, Rico MC, Solomons NW, Schümann K, et al. Associations among Inflammatory Biomarkers in the Circulating, Plasmatic, Salivary and Intraluminal Anatomical Compartments in Apparently Healthy Preschool Children from the Western Highlands of Guatemala. PLoS ONE 10(6): 2015; e0129158. pmid:26075910
  26. 26. Soto-Méndez MJ, Campos-Oliva R, Aguilera C, Solomons N, Schümann K, Gil A. Urinary osmolality of preschool children with a largely common weekday meal offering, from the western highlands of Guatemala. Spanish J Community Nutr 2014;20:13–9.
  27. 27. Battino M, Leone L, Bompadre S. High-performance liquid chromatography-EC assay of mitochondrial coenzyme Q9, coenzyme Q9H2, coenzyme Q10, coenzyme Q10H2, and vitamin E with a simplified on-line solid-phase extraction. Methods Enzymol. 2004;378:156–62. pmid:15038966
  28. 28. Rangel-Huerta OD, Aguilera CM, Martin MV, Soto MJ, Rico MC, Vallejo F, et al. Normal or High Polyphenol Concentration in Orange Juice Affects Antioxidant Activity, Blood Pressure, and Body Weight in Obese or Overweight Adults. J Nutr. 2015 Aug;145(8):1808–16. pmid:26136593
  29. 29. Drabkin Drabkin DL. The standardization of hemoglobin measurement. Am J Med Sci 1948; 215(1): 110 pmid:18919954
  30. 30. Aebi H. Catalase in Vitro. Methods Enzymol 1984;105:121–6. pmid:6727660
  31. 31. McCord JM, Fridovich I. The utility of superoxide dismutase in studying free radical reactions. I. Radicals generated by the interaction of sulfite, dimethyl sulfoxide, and oxygen. J Biol Chem 1969;244:6056–63 pmid:4981789
  32. 32. Carlberg I, Mannervik B. Glthathione Reductase. Methods Enzymol 1985;113:484–90. pmid:3003504
  33. 33. Flohé L, Günzler W. Assays of Glutathione Peroxidase. Methods Enzymol 1984;105:114–21. pmid:6727659
  34. 34. Field A. Discovering statistics using SPSS for Windows. 2nd ed. London: Sage Publications Ltd.; 2005.
  35. 35. Berk M, Williams LJ, Jacka FN, O’Neil A, Pasco JA, Moylan S, et al. So depression is an inflammatory disease, but where does the inflammation come from? BMC Med 2013;11:1–16.
  36. 36. Boy E, Bruce N, Delgado H. Birth weight and exposure to kitchen wood smoke during pregnancy in rural Guatemala. Environ Health Perspect 2002;110:109–14. pmid:11781172
  37. 37. Akohoue SA, Shankar S, Milne GL, Morrow J, Chen KY, Ajayi WU, et al. Energy expenditure, inflammation, and oxidative stress in steady-state adolescents with sickle cell anemia. Pediatr Res 2007;61:233–8. pmid:17237728
  38. 38. Drury J, Jeffers G, Cooke R. Urinary 8-hydroxydeoxyguanosine in infants and children. Free Radic Res 1998;28:423–8. pmid:9684987
  39. 39. McAuley D. GLOBALRPh 2012. http://www.globalrph.com/labs_v.htm (accessed January 07, 2014).
  40. 40. MSPAS. Encuesta Nacional de Salud Materno Infantil 2008 (ENSMI-2008/09). Guatemala: 2010.
  41. 41. Traber MG. Vitamin E Inadequacy in Humans:Causes and Consequences. Adv Nutr 2014;5:503–514. pmid:25469382
  42. 42. Winbauer AN, Pingree SS, Nuttall KL. Evaluating serum alpha-tocopherol (vitamin E) in terms of a lipid ratio. Ann Clin Lab Sci. 1999;29(3):185–91. pmid:10440582
  43. 43. Rosales F, Kang Y, Pfeiffer B, Rau A, Romero-Abal M, Erhardt J, et al. Twice the recommended daily allowance of iron is associated with an increase in plasma alpha-1 antichymotrypsin concentrations in Guatemalan school-aged children. Nutr Res 2004;24:875–87.
  44. 44. Molyneux S, Lever M, Florkowski C, George P. Plasma total coenzyme Q9 (CoQ9) in the New Zealand population: reference interval and biological variation. Clin Chem 2007;53:802–3. pmid:17405953
  45. 45. Rizzo AM, Berselli P, Zava S, Montorfano G, Negroni M, Corsetto P, et al. Endogenous antioxidants and radical scavengers. Adv Exp Med Biol. 2010;698:52–67. pmid:21520703
  46. 46. Barregard L, Møller P, Henriksen T, Mistry V, Koppen G, Rossner P Jr, et al. Human and methodological sources of variability in the measurement of urinary 8-oxo-7,8-dihydro-2'-deoxyguanosine. Antioxid Redox Signal. 2013 Jun 20;18(18):2377–91. pmid:23198723