Carotenoids are commonly ingested plant pigments found in fruits and vegetables that are responsible for red, orange and yellow colouring. Evidence suggests that a high carotenoid intake is associated with lower incidence of several chronic diseases(Reference Michaud, Feskanich and Rimm1–Reference Voutilainen, Nurmi and Mursu3). Because carotenoids are not found in significant concentrations in other foods, but are widely distributed in fruits and vegetables, they are considered to be the best current biomarker for fruit and vegetable intake(4), and levels of carotenoids in blood and tissue are known to correlate with dietary carotenoid intake(Reference Peng, Peng and Lin5). Methods for assessing carotenoid intake typically involve self-report from dietary questionnaires, which can be subject to bias and measurement error(Reference Natarajan, Flatt and Sun6), or biochemical analysis (HPLC) of plasma or tissue samples to estimate carotenoid status objectively(Reference Parker7, Reference Karppi, Nurmi and Olmedilla-Alonso8), which is invasive and expensive.
Growing evidence supports the validity of resonance Raman spectroscopy (RRS) as an alternative to HPLC for the measurement of carotenoids in living tissues, and studies around the world are now using this technology in research settings(Reference Mayne, Cartmel and Scarmo9–Reference Haag, Taskoparan and Darvin14). A novel, non-invasive technique for measuring carotenoid status in the skin using visible light, RRS utilises a small probe with a laser at a blue wavelength (λ = 488 nm) to objectively measure total carotenoid levels in skin(Reference Ermakov, Ermakova and McClane15). Our group has previously shown that RRS is a reliable (reproducible) and valid method for estimating skin carotenoid status, when compared with HPLC analyses of biopsied skin in healthy adults, with correlation coefficients ranging from r 0·7(Reference Mayne, Cartmel and Scarmo9) to 0·9(Reference Ermakov and Gellermann10). We(Reference Ermakov, Sharifzadeh and Ermakova12) and others(Reference Rerksuppaphol and Rerksuppaphol16) have been evaluating skin carotenoid status as a biomarker of carotenoid or fruit and vegetable intake in cross-sectional studies, as well as in studies of carotenoid or fruit and vegetable interventions(Reference Meinke, Darvin and Vollert17–Reference Jahns, Whigham and Johnson20).
With the exception of reference biomarkers of nutritional status (e.g. doubly-labelled water and urinary nitrogen), most other biomarkers of nutritional status can be affected by factors other than dietary intake. For example, plasma 25-hydroxyvitamin D is used as a biomarker of vitamin D status, but circulating levels are known to vary seasonally with UV exposure(Reference Ross, Taylor, Yaktine and Del Valle21–Reference Wang, Jacobs and McCullough23) and also to be affected by factors such as obesity(Reference Brock, Huang and Fraser24) and skin pigmentation(Reference McCullough, Weinstein and Freedman22, Reference Millen, Wactawski-Wende and Pettinger25), which researchers often take into account in studies using this biomarker. Given the interest in using RRS to measure skin carotenoid status as a biomarker, studies examining possible factors that may have an impact on the nutrient biomarker are needed. In our prior publication validating RRS measures against skin biopsies, we reported how mean baseline RRS values varied as a function of selected variables of interest (univariate analysis only)(Reference Mayne, Cartmel and Scarmo9). Because some of those selected variables are associated with each other, the purpose of the present analysis was two-fold: (a) to extend our prior work by performing multivariate analyses to identify significant, independent predictors of carotenoid status and (b) to examine how one RRS measure performed as a marker of status, as compared with multiple measures over time. Identification of factors (beyond carotenoid intake) that correlate with ‘usual’ carotenoid status allows investigators, when appropriate, to take into account those factors when interpreting data that use RRS as a biomarker of diet and or/nutrient status.
Subjects and methods
Study population
The Skin Antioxidant Assessment Study was a prospective investigation to examine the reliability and validity of RRS as a measure of skin carotenoid status in adults(Reference Mayne, Cartmel and Scarmo9). A total of seventy-five subjects participated in the study, which assessed the intra- and interperson variability in carotenoid status, measured by RRS within a body site (palm, inner forearm and outer forearm), across these different body sites and over time. Details of the study recruitment are reported elsewhere(Reference Mayne, Cartmel and Scarmo9). Briefly, participants were healthy adults, aged 21 to 65 years, who resided in the state of Connecticut, in the northeastern USA. Since smokers are known to have lower plasma carotenoid concentrations than non-smokers(Reference Faure, Preziosi and Roussel26, Reference Galan, Viteri and Bertrais27), we recruited both smokers and non-smokers in our sample. The present study was conducted according to the guidelines laid down in the Declaration of Helsinki, and all procedures involving human subjects were approved by the Institutional Review Board at the Yale University School of Medicine. Written informed consent was obtained from all subjects.
Data collection
Subject characteristics
Participation in the study required subjects to complete six visits at the following time points: baseline, 1 week, 2 weeks, 1, 3 and 6 months. The 6-month follow-up period allowed us to capture maximal seasonal variation in dietary intake of fruits/vegetables. At the baseline visit, a questionnaire was administered in-person by a trained interviewer to collect information on basic demographics (sex, age, marital status, education level, usual height and weight, use of tobacco and alcohol, diet, supplement use and self-identified race/ethnicity). Participants were also asked to self-identify skin colour (inner arm), under the guidance of a research assistant, using a colour wheel of skin tone samples provided by a prosthetics manufacturer (Steeper USA©). At the baseline and at each follow-up visit, the season of RRS measurement was recorded and the individuals were asked to estimate the number of cigarettes smoked per week (later collapsed into current smoking status; yes v. no) and the number of hours during the past 3 d their skin had been exposed to the sun without sunscreen.
Resonance Raman spectroscopy
Details of RRS instrumentation have been reported elsewhere(Reference Ermakov and Gellermann10–Reference Ermakov, Sharifzadeh and Ermakova12, Reference Ermakov, Ermakova and Gellermann28). To assess the reproducibility of the instrument, RRS scans were performed at baseline and at follow-up visits on all participants to assess skin carotenoid status in the palm, inner arm and outer arm over a period of 6 months. On each day when the measurements were taken, the instrument was calibrated with an external diamond standard, and laser power was measured to assess the potential drift. A single body site was scanned three times for reliability, and the mean value of the Raman measure from a body site was used in data analysis. Mean calibration-adjusted values were used in this analysis.
Dietary intake data
All participants were required to complete a self-administered FFQ from the Nutrition Assessment Shared Resource at the Fred Hutchinson Cancer Research Center (Seattle, Washington, USA) during the baseline and at the 6-month follow-up visit. This ten-page diet questionnaire has been validated in the Women's Health Initiative Dietary Assessment Study for nutrient content from foods (Pearson r 0·52 for β-carotene estimated from the FFQ and mean of 8 d of intake from 24 h recalls)(Reference Patterson, Kristal and Tinker29). Information was collected on the amount and frequency of usual food and drink consumption over the past 2 months, including questions on seasonal consumption of selected fruits and vegetables. Responses provided an estimate of daily food servings consumed, as well as nutrients associated with those servings for each study participant, including intake of total dietary carotenoids. Participants currently taking nutrient supplements at the baseline visit, such as multivitamins that may contain carotenoids, were asked to bring in bottles and continue taking the supplement for the 6-month study period. If participants were not currently consuming a nutrient supplement at baseline, they were asked to refrain from supplements until the study had concluded. No subjects reported taking single-nutrient carotenoid supplements, so the only source of supplemental carotenoids was multivitamins in the present study. The number of pills consumed per week and the nutrient content of β-carotene, lycopene and lutein in the multivitamins was recorded. Intake of total carotenoids (mg/d) was estimated by summing carotenoids consumed from food and supplements.
Statistical analysis
Mean RRS scores over time (average of six time points), indicative of ‘usual’ carotenoid status, were examined by demographic characteristics. A repeated-measures ANOVA was used to test whether mean RRS scores were significantly different by subject characteristics. To avoid having a stratum with one observation, the subject with race marked ‘other’ was coded as missing for this variable.
Main predictors of skin carotenoid status were examined in adjusted analyses at baseline by multivariate linear regression and over time (skin carotenoid status measured by RRS at six time points over 6 months) by linear mixed effects regression. The linear mixed effects models were constructed with time as a six-level class variable, random intercept and a first-order autoregressive(Reference Michaud, Feskanich and Rimm1) covariance structure for the within-subject errors. The longitudinal models took into account the timing of the RRS measurement (v. an average) and accounted for any other time-dependent variables (e.g. season of RRS measurement, recent sun exposure and current smoking status) as well as time-fixed variables assessed at baseline only (e.g. race/ethnicity, skin tone, BMI and usual diet). Because we were interested in the association of usual diet with RRS (rather than how RRS changed with diet in this non-intervention study), we used the calculated average intake of total carotenoids from the baseline and 6-month diet assessment in the models. There is evidence suggesting that energy adjustment provides a more reliable estimate of dietary intake(Reference Schatzkin, Kipnis and Carroll30). As a result, we adjusted total carotenoid intake for total energy consumed per d. The final multivariate models were based on backwards selection, with some forced predictors (e.g. BMI and current smoking status) selected a priori, based on preliminary data from our previous analysis(Reference Mayne, Cartmel and Scarmo9), as well as other literature indicating that both BMI and current smoking status correlate with plasma carotenoids(Reference Faure, Preziosi and Roussel26, Reference Galan, Viteri and Bertrais27, Reference Wang, Gaziano and Norkus31).
A cross-tabulation was constructed to compare the agreement between subjects’ baseline and usual skin carotenoid status (here, mean skin carotenoid status measured at six time points), examined by quartiles. A simple and weighted κ statistic was estimated to assess agreement. Because the six measures were not taken in equal time increments, we also examined agreement only using the average of the baseline, 3- and 6-month time points as the estimate of ‘usual’. All statistical computations were conducted using SAS version 9.2 (SAS Institute, Inc.).
Results
A full description of the study population has been reported previously(Reference Mayne, Cartmel and Scarmo9). Of the seventy-five subjects recruited for the present study, one participant came in for the baseline scan but did not complete the baseline questionnaire or any subsequent scanning visits, and therefore was excluded from all analyses. A total of seventy-four people were included in the final study sample, with a mean age of 36·6 years (median 33 years) and 62·2 % female. The majority of the population was identified as Caucasian (83·8 %) and non-smoking at baseline (85·1 %). We observed the lowest usual skin carotenoid status (here, the average of six measures without accounting for time) among African American subjects, while Asian subjects had a significantly higher usual status (P< 0·01) compared with Caucasians. Baseline smoking status (P= 0·07) and dark skin tone (P= 0·06) were both marginally significant predictors of lower usual skin carotenoid status in univariate analyses.
We examined the main predictors of skin carotenoid status at baseline (one measure only) and over time (longitudinally, using all six independent measurements to account for time) in multivariate analysis (Table 1). We observed that many of the same factors were independent predictors of skin carotenoid status both at baseline and over time, with adjustment for covariates. After taking into account dietary intake of total carotenoids, which as expected was significantly predictive of higher skin carotenoid status at baseline and over time (P< 0·01 for both), we observed that season of measurement and race were also independently predictive of skin carotenoid status, measured by RRS. In particular, significantly lower skin carotenoid status was observed when RRS measurements were obtained during the summer months (P= 0·05 at baseline; P< 0·01 over time), followed by spring/autumn measurements (not significant at baseline; P< 0·01 over time), as compared with the winter season. Asian race was significantly associated with higher skin carotenoid status at baseline and over time (P< 0·01 for both), while African American race was associated with a non-significant lower status; Caucasian race was the comparison group. Current smoking status was marginally associated with lower carotenoid status at baseline (P= 0·07). Finally, self-reported recent sun exposure (number of hours in the 3 d prior to RRS measurement when skin had been exposed to the sun without sunscreen) emerged as a significant predictor of skin carotenoid status over time (P= 0·01), while it was not a significant predictor when using only one measure of sun exposure taken at the baseline visit (P= 0·96).
Main predictors of skin carotenoid status at baseline and over time*, with adjustment for covariates (Mean values and standard deviations; number of subjects and percentages)
N/A, not applicable; RRS, resonance Raman spectroscopy.
* Skin carotenoid status measured by RRS at six time points over 6-months (baseline, 1-week, 2-week, 1-month, 3-month and 6-month follow-up visits).
† Including carotenoid intake from supplements; adjusted for mean daily energy intake (average of baseline and 6-month estimates).
‡ To make β-coefficients more interpretable, they are scaled to the interquartile range for each variable as follows (total carotenoids per 10 mg/d, recent sun exposure per 2 h/d and BMI per 5 kg/m2).
§ Mean RRS values for total study population at baseline was 1·55 (sd 0·57) counts.
We observed good agreement between skin carotenoid status measured by RRS at baseline and multiple time points. Of the eighteen subjects categorised in the highest quartile at baseline, fourteen remained in the highest group when carotenoid status was estimated as the average of six time points. In all, sixteen of the eighteen subjects categorised as having the lowest carotenoid status at baseline remained in the lowest category when their status was averaged over time. Overall, fifty-six of the total seventy-four subjects remained in the same category of skin carotenoid status, and the weighted κ showed an 80 % agreement (95 % CI 0·72, 0·88) between the measures (Table 2). When we repeated the analysis, restricting the data to baseline, 3- and 6-month follow-up visits (so the usual measure of carotenoid status was less influenced by the earlier measures), the weighted κ showed a 69 % agreement (95 % CI 0·59, 0·80) between the measures (Table 3).
Agreement between baseline and usual* skin carotenoid status quartiles, measured by resonance Raman spectroscopy at six time points (n 74)†
Q, quartile.
* Mean skin carotenoid status in the palm measured at six time points; Q4 is highest quartile (baseline, 1-week, 2-week, 1-month, 3-month and 6-month follow-up visits).
† Simple κ = 0·68 (95 % CI 0·54, 0·81); weighted κ = 0·80 (95 % CI 0·72, 0·88).
Agreement between baseline and usual* skin carotenoid status quartiles, measured by resonance Raman spectroscopy at three time points (n 74)†
Q, quartile.
* Mean skin carotenoid status in the palm measured at three time points (baseline, 3-month and 6-month follow-up visits); Q4 is highest quartile.
† Simple κ = 0·50 (95 % CI 0·35, 0·64); weighted κ = 0·69 (95 % CI 0·59, 0·80).
Discussion
We previously reported on the predictors of skin carotenoid status in univariate analysis at baseline only(Reference Mayne, Cartmel and Scarmo9). Because many of those covariates were correlated, here, we report the results of a multivariate analysis to identify predictors of carotenoid status, with adjustment for other factors that correlate with the RRS biomarker. As expected, many of the predictors identified for usual carotenoid status (e.g. over time) were the same as those identified at one time point, including self-reported intake of total carotenoids and summer season of RRS measurement.
In the present analyses, we used an average of the baseline and 6-month dietary intake, which minimises the ability to detect changes in an individual's diet over time but allows us to capture usual diet. Evidence suggests that there are seasonal effects on fruit and vegetable consumption, with variation in the availability and variety of these foods(Reference Locke, Coronado and Thompson32, Reference Ziegler, Wilcox and Mason33), which is probably reflected in overall carotenoid status in blood and tissue. The half-life of carotenoids in skin is believed to be longer than the half-life in blood (the latter being approximately 2 weeks)(Reference Rock, Swendseid and Jacob34), and our own intervention data suggest that skin reflects carotenoid intake over the past 8+ weeks (with values still increasing at 8 weeks after feeding subjects a high-carotenoid diet)(Reference Jahns, Whigham and Johnson20). Therefore, the lower RRS values that we observed during the summer months (June, July and August) may, in part, reflect dietary intake from the late winter/early spring months (March, April and May), when the availability and variety of fruits and vegetables is more limited in the northeast USA.
In multivariate analysis, self-reported sun exposure (queried at all six time points) emerged as a significant predictor of skin carotenoid status over time, and one that was not significant using baseline-only sun exposure (one time point) in relation to baseline RRS measures. Few studies have examined the effect of UV radiation on carotenoid concentrations in human skin, and most assessed β-carotene only with a variety of measurement techniques, making it difficult to directly compare study results(Reference Goralczyk, Wertz, Britton, Liaaen-Jensen and Pfander35). Experimental studies suggest that skin carotene and retinol concentrations decrease immediately following UV radiation in hairless rodents, but levels can be recovered with carotenoid supplementation(Reference Someya, Totsuka and Murakoshi36). There is also evidence from human studies to suggest that the degradation of skin carotenoids following exposure to UV radiation may be differential, with more lycopene being degraded than β-carotene(Reference Ribaya-Mercado, Garmyn and Gilchrest37), possibly serving as the first line of defence in human skin photoprotection, with lycopene being a predominant carotenoid in human skin(Reference Scarmo, Cartmel and Lin38). Because the present findings suggest that recent UV radiation from sun exposure is correlated with lower carotenoid status at a localised body site, it may be helpful to consider adjusting for sun exposure in future studies using RRS as a biomarker. It is important to highlight that our study subjects were queried on recent sun exposure (only reporting on the 3 d prior to RRS scans); we observed a range of 3-d sun exposure measured at baseline from zero to eighteen total hours, with the average 3-d sun exposure equalling to 3 h. While recent sun exposure was significantly inversely correlated with RRS, the magnitude of the effect (β coefficient) was modest, suggesting a relatively minor overall effect of UV radiation on RRS measures, at least under the conditions of sun exposure that we observed. Because we had no measure of overall sun exposure over the duration of 6 months, our multivariate models could only control for recent exposure, which may or may not be correlated with an individual's usual sun exposure. Thus, the lower levels of skin carotenoids that we observed when RRS measurements were taken in summer could also, in part, be due to greater sun exposure. Future studies of RRS in subjects with extensive sun exposure are needed to better understand the impact of sun exposure on the RRS biomarker of skin carotenoid status.
A significant association was observed between self-identified race and skin carotenoid status. Asian subjects had significantly higher skin carotenoid status, compared with white and African American subjects; however, our sample included only six Asian subjects. This may be a chance finding given the small sample size, but it is plausible that the higher status in Asians reflects their lifestyle patterns or genetic differences. While genetics might also explain some of the low skin carotenoid status in our African American subjects, another possibility is greater melanin content or differences in body composition (adiposity). With our limited sample size, consisting of only three individuals with dark skin tone and only five individuals who identified as African American, we were unable to examine this further with multivariate modelling. However, as previously reported(Reference Mayne, Cartmel and Scarmo9), we examined the ratio of measured RRS values in the outer arm (area of most skin pigmentation) v. the palm (area of least skin pigmentation) among participants who were identified as dark skinned, as compared with similar ratios obtained from subjects with light skin. We observed no evidence that the RRS ratios varied significantly with skin pigmentation, as might be expected if melanin had a filtering effect on the laser(Reference Mayne, Cartmel and Scarmo9). Nevertheless, we measured the palm of the hand in our studies, with one of the reasons being that the dermal melanin is light and less variable among subjects of different races and ethnicities. Additional research in racially diverse populations with larger sample sizes and objective quantitation of skin pigmentation is needed to explain the true effect, if any, of melanin on the RRS measures.
The present results indicate that a single RRS measure of the palm of the hand was a reasonably good indicator of a subjects’ usual skin carotenoid status. Multiple measures did not greatly improve the ability of the biomarker to objectively categorise carotenoid status. In terms of feasibility of use, a single RRS measure can save follow-up efforts in studies of the biomarker and health outcomes. Among those participants for whom one measure of skin carotenoid status misclassified the usual skin carotenoid status (by quartiles), the degree of misclassification did not exceed one quartile in difference. This suggests that a single RRS measure can be used in future studies to identify subjects with low carotenoid status, who may be at risk for subsequent nutrition-related diseases.
A strength of the present study was the prospective design, which provided up to six repeated measurements for each participant over the course of a 6-month follow-up in conditions with marked seasonal variation in climate and availability of fresh fruit and vegetables. A limitation of the present study was the use of self-reported data for estimating sun exposure, as well as the lack of an objective measure of melanin. The primary limitation was the small sample size, consisting primarily of Caucasian non-smoking women.
In summary, results from the present study suggest that a single measure of skin carotenoids by RRS is a reasonably robust indicator of usual carotenoid status. We also identified factors, such as season of measurement and sun exposure, which correlate with skin carotenoid status. These findings contribute to a growing body of research advancing the development of this innovative nutrition-related biomarker.
Acknowledgements
We wish to thank the participants in the present study. This study was supported by grants R01 CA096838 (S. T. M.), sub-award no. A05202 and M03A00158 (I. V. E. and W. G.) and R01 EY11600 (P. S. B.) from the National Institutes of Health, USA; and an unrestricted departmental grant from the Research to Prevent Blindness (New York, NY, USA) (P. S. B.). All authors have made substantial contributions and the final approval of the conceptions, drafting and the final version of this manuscript. The authors' contributions are as follows: S. S. was involved in the recruitment of the study subjects and data collection, statistical analysis and manuscript drafting and review; B. C. contributed to the experimental design, supervision of the human study, data collection and manuscript review; H. L. contributed to the experimental design, advised on the statistical analysis and manuscript review; D. J. L. contributed to the experimental design, clinical support and manuscript review; I. V. E. provided spectroscopy support and assisted with training and data interpretation, and manuscript review; P. S. B. was involved in the experimental design and manuscript review; W. G. participated in the design and construction of the RRS system, oversaw training for the data collection and manuscript review; S. T. M. contributed to the experimental design, acquiring funding and manuscript review. W. G. and P. S. B. hold patents for the methods described in the present paper. The authors report no other competing interests.
