Assessing urinary flow rate, creatinine, osmolality and other hydration adjustment methods for urinary biomonitoring using NHANES arsenic, iodine, lead and cadmium data
© The Author(s). 2016
Received: 26 February 2016
Accepted: 30 May 2016
Published: 10 June 2016
There are numerous methods for adjusting measured concentrations of urinary biomarkers for hydration variation. Few studies use objective criteria to quantify the relative performance of these methods. Our aim was to compare the performance of existing methods for adjusting urinary biomarkers for hydration variation.
Creatinine, osmolality, excretion rate (ER), bodyweight adjusted ER (ERBW) and empirical analyte-specific urinary flow rate (UFR) adjustment methods on spot urinary concentrations of lead (Pb), cadmium (Cd), non-arsenobetaine arsenic (AsIMM) and iodine (I) from the US National Health and Nutrition Examination Survey (NHANES) (2009–2010 and 2011–2012) were evaluated. The data were divided into a training dataset (n = 1,723) from which empirical adjustment coefficients were derived and a testing dataset (n = 428) on which quantification of the performance of the adjustment methods was done by calculating, primarily, the correlation of the adjusted parameter with UFR, with lower correlations indicating better performance and, secondarily, the correlation of the adjusted parameters with blood analyte concentrations (Pb and Cd), with higher correlations indicating better performance.
Overall performance across analytes was better for Osmolality and UFR based methods. Excretion rate and ERBW consistently performed worse, often no better than unadjusted concentrations.
Osmolality adjustment of urinary biomonitoring data provides for more robust adjustment than either creatinine based or ER or ERBW methods, the latter two of which tend to overcompensate for UFR. Modified UFR methods perform significantly better than all but osmolality in removing hydration variation, but depend on the accuracy of UFR calculations. Hydration adjustment performance is analyte-specific and further research is needed to establish a robust and consistent framework.
KeywordsBiomonitoring Hydration adjustment Creatinine Osmolality Urinary flow rate NHANES
Urinary biomonitoring is the preferred method of exposure and nutritional assessment for many chemical elements and metabolites given its non-invasiveness, logistical appeal and ease of measurement with modern analytical techniques . This is true for potentially harmful elements like arsenic (As) , essential nutrients like iodine (I)  and drug and organic compounds [4, 5], making the meaningful interpretation of urinary data a requirement with implications for public health, occupational health and forensic applications. The US National Health and Nutrition Examination Survey (NHANES) has proved an invaluable, growing resource of chemical biomonitoring data , but the value of such data depend on their correct interpretation , challenges with which are currently limiting the full potential of urinary chemical biomarkers .
Urinary analyte concentrations are susceptible to variation from factors extending beyond exposure and are categorised  as follows: (i) time of sampling relative to exposure; (ii) inter-individual toxico-kinetic factors and (iii) physiological characteristics of the biomonitoring matrix. While the first two factors should not be ignored, the third, specifically the variation in dilution among spot urine samples, is addressed here.
While collection of 24 hr urine samples is preferred, it is not feasible for large biomonitoring studies due to resource limitations, cumbersome sample nature and volunteer compliance issues . First morning void (FMV) or spot collections are common substitutes, but are limited in that they reflect the hydration status of the individual at the time of collection thus may differ markedly in dilution as a result of differences in urinary flow rate (UFR). Spot/FMV samples are nevertheless widely deemed acceptable provided that the effect of sample dilution is quantified and appropriately adjusted . Several methods for adjusting spot/FMV data are currently employed but there is no consensus on which is the most appropriate.
The most common technique employed is creatinine adjustment, whereby urinary analyte concentrations are ratioed to creatinine concentrations. This method implicitly assumes that urinary creatinine is excreted at a constant rate and varies only as a function of UFR. However, these assumptions are of questionable validity, since creatinine concentrations have been shown to depend upon all of demographic group , protein intake , muscle mass  and malnutrition .
Alternative methods such as specific gravity (SG) or osmolality adjustment are commonly reported. Close agreement has been demonstrated between these methods  but osmolality, as measured by osmometry, has been described as the definitive measure of urinary concentration  despite being previously considered prohibitively expensive . Osmometry is not susceptible to the same interferences as SG, conventionally measured by refractometry, which may be confounded in subjects with, for example, proteinuria and glucosuria. Urinary osmolality was a post-2008 inclusion in NHANES and, while similar factors affecting creatinine excretion were found to be responsible for variation in urinary osmolality, less influence was observed on osmolality than on creatinine in the US population .
Adjusting urinary biomonitoring results according to Eqs. 2 and 3 was recently proposed to directly account for hydration status and address demographic variations in UFR more effectively than creatinine and osmolality adjustments . While UFR, providing accurate measurement of time and volume, is more reflective of hydration status than surrogate measures, such as creatinine and osmolality, its application in Eqs. 2 and 3 directly incorporates hydration bias into results. Given that ERs can still apply to restricted time periods, they are dependent on UFR, i.e. hydration, at that time. Strong positive Spearman’s correlation coefficients (rs) have been reported between ERs of urinary analytes and UFR [20, 21], indicating that the adjustment in Eq. 2 is not theoretically robust.
where a and b (referred to here as Araki’s b value) are analyte-dependent, empirically determined regression coefficients.
where Araki’s b values were derived for a number of analytes using multiple voids from single individuals subjected to water loading and water restrictive conditions. This UFR adjustment was found to be more effective than ER, creatinine and SG adjustment in removing UFR-dependent variation from adjusted urinary analyte concentrations [20, 21]. For datasets, such as NHANES, that do not contain extensive analyte concentration data for multiple voids from single individuals, a previously reported iterative method  may be used to calculate population-level Araki’s b values by optimising appropriate performance criteria.
Suggested criteria for assessing performance of urinary biomonitoring adjustment methods
Correlation between adjusted spot analyte concentrations and UFR.
Weaker correlations indicating good performance.
Correlation between adjusted spot analyte concentration and an independent measure of internal dose, e.g. analyte concentration in blood.
Stronger correlations indicating good performance.
Correlation of spot analyte concentrations with analyte excretion over 24 hr or composite 24 hr concentrations.
Closer agreement/lower variation in spot samples indicating good performance.
Correlation of spot analyte concentration with an independent measure of/proxy for external exposure e.g. drinking water analyte concentration.
Stronger correlations indicating good performance.
This paper aims to compare the performance of urinary biomonitoring hydration adjustment techniques using NHANES (2009–2010, 2011–2012) spot urinary concentrations of selected chemical analytes and, in particular, test whether or not analyte-specific UFR adjusted concentrations provide a more robust adjustment than creatinine, osmolality, ER or ERBW adjustments. Arsenic and iodine were selected as chemicals on which to make these tests as respectively toxic and essential elements for which urinary biomonitoring is widely used. Additionally, Pb and Cd were selected for study, due to the availability of paired urine-blood samples in the NHANES database and the applicability of criterion B to these elements. This provides the opportunity to compare the adjustment performance characteristics of two independent assessment criteria.
Data from the NHANES 2009–2010 and 2011–2012 surveys were acquired from the NHANES website . Volunteer consent information and dataset access can be found online: http://www.cdc.gov/nchs/nhanes.htm. Data on demographic variables; body measurements; standard biochemistry profile; diabetes; kidney conditions; plasma fasting glucose; urinary flow rates; urinary creatinine; urinary osmolality; urinary metals; total and speciation urinary As; urinary I and blood metals were downloaded in SAS (.xpt) format. Data were converted to MS Excel (.xlsx) format using the R programming environment SASxport and xlsx packages [30, 31], before being matched by sequence number (SQN) in MS Access. Volunteers with data present on gender, age, bodyweight, urinary creatinine, urinary osmolality, UFR, As speciation and blood metals, of either non-Hispanic white, non-Hispanic black or Mexican American ethnicity were initially included.
Volunteers with evidence of health conditions that could affect the performance of urinary adjustment calculations were excluded using previously published criteria [10, 19]: urinary albumin-creatinine ratio >30 mg/g creatinine was treated as albuminuria and diabetics were identified by self-reported physician diagnosis or plasma glucose ≥126 mg/dL (≥8 h fasting) or ≥200 mg/dL (<8 h fasting). Chronic kidney disease (CKD) was identified by self-reported physician diagnosis or an estimated glomerular filtration rate (eGFR) <60 mL/min/1.73 m2 using the Chronic Kidney Disease Epidemiology Collaboration (CKD-EPI) equation . For volunteers aged <18, eGFR was estimated using the Bedside Schwartz equation . Finally, volunteers without detectable concentrations of urinary Pb, Cd, Total As and I were excluded to limit the effects of censoring on analyses.
Additional data from a single volunteer consisting of multiple spot Cd concentrations and UFR measurements were reproduced  for observational purposes only.
Detailed analytical methodologies for the analytes investigated in this paper, plus other analytical components of NHANES, are reported online: http://www.cdc.gov/nchs/nhanes/nhanes2011-2012/lab_methods_11_12.htm. Urine samples were provided by volunteers at the NHANES mobile examination center (MEC). Volunteers were asked to fully void the bladder and report the last time of previously doing so. Volumes of the samples provided at the MEC were measured and used, with the previous void times reported by volunteers, to calculate UFR as per Eq. 1. For volunteers with initial urinary volumes below requirement, subsequent voids were collected and composite UFRs were calculated using the total volumes and times covered by all voids. This ensured that laboratory measurements made on pooled samples consisting of multiple voids corresponded to the correct UFRs. Urinary osmolality was measured using freezing-point depression (cryoscopic) osmometry performed with an Osmette II, Model 5005 Automatic Osmometer (Precision Systems Inc.). Urinary creatinine was determined using an enzymatic (creatininase) reaction and a Roche/Hitachi Modular P Chemistry Analyzer. Urinary total I, Pb and Cd and whole blood Pb and Cd were determined using inductively coupled plasma dynamic reaction cell mass spectrometry (ICP-DRC-MS) (PerkinElmer ELAN® 6100 DRCPlus or ELAN® DRC II). Urinary As speciation was performed using high performance liquid chromatography (HPLC) coupled to ICP-DRC-MS. Combined urinary inorganic As and methylated metabolites (AsIMM) was calculated as the sum of arsenous acid (AsIII), arsenic acid (AsV), monomethylarsonic acid (MMA) and dimethylarsonic acid (DMA) species, as this is the routine biomarker of As exposure and does not incorporate non-toxic arsenobetaine. Osmolality and UFR were both determined at the MEC shortly after urine collection. The remaining measurements were made after urine and blood samples had been frozen at -20 °C and shipped to relevant laboratories where they remained frozen until analysis to prevent evaporation and the inter-conversion of As species.
Urinary analyte adjustment calculations
Training dataset: 80 %, reserved for Araki’s b value derivation.
Testing dataset: 20 %, for applying and assessing the performance of adjustment calculations.
The partition of 80:20 % was deemed suitable  and was selected to (i) ensure sufficient training data were available; (ii) retain a testing dataset of a size comparable to that of a biomonitoring study in which these kinds of adjustments may be employed and (iii) to preserve the distribution of demographic and analytical variables between both datasets.
where, for consistency with recent publications , Osmref is the median osmolality (mOsm/kg) of training data volunteers (734 mOsm/kg) and Osmmeas is that measured in the individual specimen. Araki’s b values were extracted using an adaption of a previously published approach  which involved a simple numeric method. Pearson correlation coefficients (rp) for criterion A and criterion B were calculated from the training data for values of Araki’s b from 0 to 1.5 at intervals of 0.01. Araki’s b values that yielded optimum correlations for criterion A (minimizing absolute value of rp) and separately for criterion B (maximizing rp) were determined. Araki’s b values were also derived on demographic subsets of specific age groups and specific genders/ethnicities to identify patterns between groups. Optimum Araki’s b values were used to perform Araki’s modified UFR adjustment (Eq. 5), henceforth referred to as UFRA and UFRB when adjusted using optimum Araki’s b values for criteria A and B, respectively. An R script has been provided in Additional file 1 to allow other groups to derive Araki’s b values and perform hydration adjustments.
Due to the specific application of NHANES data in assessing adjustment methods rather than making inferences of biomonitoring measurements in the US population, sample weights were not incorporated into analyses. Statistical tests (and graphical presentations) were performed using R version 3.0.0 (base package) . Urinary and blood analyte data were positively skewed and, therefore, geometric means (GM) were calculated as opposed to arithmetic means. For the same reason, Pearson correlations of urinary analyte concentrations against UFR and blood analyte concentrations were calculated on natural log (ln) transformed data with significance tests (p-values) and 95 % confidence intervals (CI) using the ‘cor.test’ function. Pearson’s, as opposed to Spearman’s, correlation was selected to prevent the loss of information that occurs when data are reduced to ranks in the process of calculating Spearman’s correlation. It was necessary to test the significance of the difference between correlations of, for example, urinary Pb with blood Pb adjusted by different methods. These correlations are not independent (because of the common variable, blood Pb), so the Williams’s test  was performed using the r.test function in the psych package . Point density contour lines were added to plots using two-dimensional kernel density estimation in the MASS package .
Exploratory analyses – training data
Demographic characteristics and unadjusted analyte geometric means (GM) and ranges for training and testing datasets
Demographic group, n (%)
6–11 years old
12–19 years old
20–39 years old
40–59 years old
>60 years old
Analytical measurement, GM (range)
Urinary AsIMM, μg/L
Urinary I, μg/L
Urinary Pb, μg/L
Urinary Cd, μg/L
Blood Pb, μg/dL
Blood Cd, μg/L
Derivation of Araki’s b values – training data
Araki’s b values derived for Pb, Cd, AsIMM and I in the present study (NHANES 2009–2012 (CDC, 2015) training dataset) compared with previously reported literature values
Araki’s b value
Criterion A (UFR) optimised (present study)
Criterion B (Blood) optimised (present study)
Araki et al. (1986)a,c
Araki et al. (1990)b,c
Comparison of adjustment methods – testing data
Geometric means (GM) and ranges of urinary analytes following adjustment by various methods (testing dataset)
Urinary analyte, GM (range)
Creatinine-adjusted, μg/g creatinine
Osmolality-adjusted, μg/L, 734 mOsm/kg
UFRA, μg/L, UFR 1 mL/min
UFRB, μg/L, UFR 1 mL/min
Pearson correlations for performance Criterion A across the range of adjustment methods investigated for NHANES 2009–2012 (CDC, 2015) (testing dataset)
rp (95 % CI)
-0.33*** (-0.41, -0.24)
-0.25*** (-0.34, -0.26)
-0.37*** (-0.45, -0.28)
-0.39*** (-0.47, -0.31)
0.18*** (0.09, 0.27)
0.18*** (0.09, 0.27)
0.09 (-0.001, 0.19)
-0.001 (-0.10, 0.09)
0.02 (-0.07, 0.12)
-0.09 (-0.18, 0.004)
0.52*** (0.45, 0.59)
0.47*** (0.39, 0.54)
0.70*** (0.65, 0.74)
0.45*** (0.37, 052)
0.43*** (0.35, 0.50)
0.42*** (0.34, 0.50)
0.60*** (0.53, 0.66)
0.35*** (0.26, 0.43)
0.01 (-0.08, 0.11)
-0.01 (-0.10, 0.09)
-0.03 (-0.12, 0.07)
-0.01 (-0.10, 0.09)
0.18*** (0.09, 0.27)
0.22*** (0.13, 0.31)
Pearson correlations for performance Criterion B across the range of adjustment methods investigated for NHANES 2009–2012 (CDC, 2015) (testing dataset)
rp (95 % CI)
0.67 (0.61, 0.72)
0.58 (0.51, 0.64)
0.79 (0.75, 0.82)
0.65 (0.59, 0.70)
0.81 (0.78, 0.84)
0.66 (0.60, 0.71)
0.59 (0.52, 0.64)
0.63 (0.57, 0.68)
0.57 (0.50, 0.63)
0.74 (0.70, 0.78)
0.62 (0.56, 0.67)
0.75 (0.70, 0.79)
0.62 (0.56, 0.68)
Pb: UFRA, Osmolality > Creatinine, UFRB > Unadjusted > ERBW > ER
Cd: UFRA, Osmolality > Creatinine, UFRB > Unadjusted > ERBW > ER
AsIMM: UFRA > Osmolality > Creatinine > Unadjusted > ERBW > ER
I: UFRA > Osmolality > Creatinine > ERBW > Unadjusted > ER
Pb: Osmolality ≥ Creatinine ≥ UFRB ≥ UFRA > ER ≥ Unadjusted ≥ ERBW
Cd: Osmolality ≥ Creatinine ≥ UFRB ≥ UFRA ≥ ER ≥ Unadjusted ≥ ERBW
Irrespective of whether criterion A or criterion B was used to assess adjustment method performance, it was evident that UFRA, UFRB, creatinine and osmolality adjustment methods all provided for a statistically significant improvement relative to unadjusted analyte concentrations. Of these, osmolality adjustment was determined to be the optimal adjustment method except for AsIMM and I for which UFRA showed a marginally better performance. The criterion A-based performance of osmolality and UFRA adjustments were equally good for Pb and Cd. Indeed, osmolality adjustment resulted in a weak (r p = 0.10) significant (p < 0.05) correlation only in the case of osmolality adjusted AsIMM versus UFR. Creatinine adjustment, in contrast, yielded significant positive correlations with UFR for Pb (r p = 0.18), Cd (r p = 0.18) and AsIMM (r p = 0.32). Iodine was an exception, with no significant correlation of creatinine adjusted I concentrations against UFR.
Excretion rate adjustment methods (ER, ERBW) performed worse than any of UFRA, UFRB, creatinine or osmolality adjustments according to both criteria. Furthermore, although ER and ERBW adjustments removed observed negative correlations of unadjusted analyte concentrations with UFR, they mostly resulted in positive correlations of an equal or greater magnitude. Excretion rate and ERBW adjustments thus performed no better than implementing no adjustment at all, with the sole exception of I, for which ERBW adjustment performed marginally better (r p = 0.35 cf. -0.39).
Araki’s b values derived for AsIMM, I, Pb and Cd on specific demographic sub-groups of the present study group (training dataset)
Optimum Araki’s b value
Non-Hispanic white male
Non-Hispanic white female
Non-Hispanic black male
Non-Hispanic black female
Mexican American male
Mexican American female
6–11 years old
12–19 years old
20–39 years old
40–59 years old
>60 years old
Osmolality and UFRA adjustment methods provided the best or near best performance of the adjustment methods tested using criterion A. Osmolality and creatinine adjusted concentrations yielded the strongest correlations using criterion B. In its nature, UFRA is tailored to optimise criterion A performance, however, osmolality and creatinine methods tended to perform better than Araki’s UFR-based adjustment methods with respect to criterion B. This may reflect greater uncertainties in volunteer reported times of initial voids than in uncertainties in objectively measured UFR surrogates such as osmolality or creatinine. The issue of reliance on the accuracy of volunteer reported void times when calculating UFR in NHANES has been raised previously . We note this as a limitation of the present study, in that no efforts were made to refine UFR data or quantify their accuracy.
Alternatively, differences in Araki’s b values for given analytes, between individuals of different ages, genders and ethnicities may partly account for deficiencies in UFR based adjustments. Attempts were made to derive Araki’s b values on demographic subsets of the training dataset. Differences in optimum values were observed between groups but this was possibly an artefact of different group sizes. When age-specific b values were implemented to the adjustment of urinary Pb and Cd, no significant improvements were observed in Criterion B correlations. Future efforts should be made to derive Araki’s b values for multiple individuals, collecting multiple voids at various states of hydration, to more closely represent the relationship between analyte concentrations and UFR, as illustrated in Fig. 1f. As proposed previously , these derivations should be made on specific demographic groups to investigate whether the relationship between analyte concentrations and UFR vary in a characteristic manner, something which was not achievable given the constraints of the data utilised in the present study.
The finding that Araki’s b values calculated here by the UFRA method were generally lower than those calculated previously [20, 21], might give rise to questioning of the validity of the UFRA based values. We therefore identify a potential flaw in the validity of criterion A when using the present dataset. As discussed, it was not possible to derive Araki’s b values in the conventional manner using the NHANES dataset. This requires multiple voids from single volunteers at different hydration states, circumstances under which the ‘true’ relationship between UFR and analyte concentration is observed due to the relatively constant internal dose of a given analyte, for a single individual, over the timescale investigated. NHANES data consist of single voids from multiple individuals with greater inter-individual ranges of internal doses which have the potential to alter the observed slope (b value). Under circumstances where the standard deviation of the distribution of internal dose is considerably bigger than that of the distribution of UFR for the studied population, the calculated Araki’s b value may be positively biased. This in-turn determines the optimum value for criterion A - the b value that describes the slope between UFR and analyte – evident in the agreement between the b values that describe the slopes presented in Fig. 1 and those derived using the numeric method. This is further illustrated by the difference between criteria A and B optimum values for Pb and Cd (Fig. 2a and b), the criterion B optimums could be considered more robust for the NHANES dataset used here as they are independent of UFR. We note, however, that this paper does not purport to suggest specific Araki’s b values for application elsewhere, but is successful in reiterating proof of concept of their necessary implementation to remove hydration variation from spot analyte concentrations and, in doing so (Fig. 2a and b), better reflect internal dose.
Our findings reiterate the analyte-specific nature of hydration adjustment, exemplified by the difference in criteria A correlations between creatinine adjusted AsIMM and I. Creatinine adjusted I concentrations yielded no significant correlation with UFR, whereas creatinine adjusted AsIMM concentrations were more strongly correlated with UFR than I, Pb and Cd. Inspection of Fig. 1c indicates differing relationships with UFR of urinary AsIMM for high-DMA and low-DMA samples – this suggests that the biochemistry of different species of the same element (and different elements) influence Araki’s b values.
As noted previously , for studies outside the NHANES framework or similar population-scale biomonitoring designs, the collection of UFRs may provide additional value providing that accurate recordings of time and volume are obtained. It is recognised that this may not be logistically feasible for all studies and surrogates such as creatinine and osmolality are attractive alternatives. Modifications of these surrogate based adjustments have also been explored [22, 23, 42] using methodologies based on the work of Araki et al. (1986). These approaches, such as modified SG  or modified creatinine  adjustment, were not addressed in the present study and further exploration of such alternatives may prove valuable for studies that are restricted in their ability to directly measure UFR, particularly in low-budget circumstances or developing countries. Similarly, measurements of additional urinary constituents that are indicative of medical conditions, such as glucose, protein, ketones and bilirubin, could have provided additional data exclusion criteria had they been available to us. A comparison of the performance of different adjustment methods between volunteers with and without the presence of such analytes, and the medical conditions that were available as exclusion criteria (e.g. CKD and diabetes) will make for an important matter of further research.
The implications of the findings made in this investigation, and the questions that remain unanswered, have implications for environmental and epidemiological studies using urinary biomonitoring to assess human exposures and investigate the dose-response relationships between environmental chemicals and health end-points. The differences in biomarker levels yielded by different adjustment methods is evident (Table 4). This impacts the interpretation of results when making comparisons to existing guidance or reference values. Furthermore, it impacts the derivation of reference values themselves. For example, a large and much needed body of work has been undertaken to derive biomonitoring equivalents to be used in comparison with various urinary biomarkers . Some derivations have utilised creatinine adjustment, which may have limited their applicability to studies using alternative adjustment methods, or other studies using creatinine adjustment with different demographic structures. This problem also extends to studies that explore relationships between urinary analyte concentrations and health outcomes- a widely used application of NHANES data [44–46]. A robust, standardised framework of urinary hydration adjustment is warranted to ensure the validity and sensitivity of such analyses.
Osmolality consistently performs as the best or near best adjustment method against two performance criteria: minimum correlation of adjusted urinary analyte concentration with UFR; maximum correlation of adjusted urinary analyte concentration with blood analyte concentration.
The method of Araki et al. (1986, 1990) for objectively determining Araki’s b values to adjust urinary analyte concentrations also performs well, but is limited by the requirement for accurately determined UFR data and age/gender/ethnicity specific b values.
Creatinine adjustment methods may be suitable for some analytes (e.g. I) that have similar b values to creatinine, but can otherwise result in significant biases.
ER and ERBW based adjustment methods are shown here to overcompensate for UFR and invariable performed worse than osmolality, creatinine, UFRA, UFRB and often worse than unadjusted concentrations.
Thus, we demonstrate that conventional application of UFR is limiting the full potential of this metric. The under-performance of both ER and ERBW in relation to two independent performance criteria support previous findings [20, 21] that using UFR to calculate excretion rates in the conventional manner propagates inaccurate results. The inclusion of Araki’s b values into adjustment calculations was demonstrated to significantly improve the performance of UFR adjustment for both criteria relative to ER and ERBW adjustments.
The derivation of specific Araki’s b values requires substantial further work by collecting multiple voids from the same individuals. By compiling a range of Araki’s b values for a range of analytes and demographic characteristics, their determining factors can be assessed, enabling the development of more sophisticated framework for urinary biomarker adjustment. Finally, additional constraints on the interpretation of urinary biomarker concentrations need addressing, such as time of sampling relative to exposure , which hydration adjustment cannot overcome.
AsIII, arsenous acid; AsIMM, Inorganic arsenic and Methylated Metabolites; AsV, arsenic acid; BW, bodyweight; CI, confidence interval; CKD, chronic kidney disease; DMA, dimethylarsonic acid; eGFR, estimated glomerular filtration rate; ER, excretion rate; ERBW, excretion rate adjusted for body weight; FMV, first morning void; GM, geometric mean; HPLC, high performance liquid chromatography; ICP-DRC-MS, inductively coupled plasma dynamic reaction cell mass spectrometry; MEC, Mobile Examination Center; MMA, monomethylarsonic acid; NHANES, National Health and Nutrition Examination Survey; rp, Pearson’s correlation coefficient; rs, Spearman’s correlation coefficient; SG, Specific Gravity; UFR, urinary flow rate; UFRA, urinary flow rate adjustment optimised to criterion A; UFRB, urinary flow rate adjustment optimised to criterion B
We gratefully acknowledge Professor Andrew Meharg for the provision of comparison data to compliment our investigation and Dr Simon Chenery and Dr Mark Cave for scientific review and advice. We thank Olivier S. Humphrey and Ahmed A. N. Al Bualy for testing the R script provided in Additional file 1.
Funding was provided by the Natural Environment Research Council (NERC) via a University of Manchester/BUFI (Centre for Environmental Geochemistry) studentship (Contract No. GA/125/017, BUFI Ref: S204.2).
Availability of data and material
The data forming this investigation is from NHANES 2009–2010 and 2011–2012. Data are available for download from the NHANES website: http://www.cdc.gov/nchs/nhanes.htm.
DRSM, MJW and DAP conceived the investigation. DRSM acquired the data and performed statistical analyses. RML coordinated the data management statistical methodologies employed and CJM constructed the R programming scripts to conduct analyses. All authors contributed to and reviewed the manuscript.
DRSM is completing a PhD on human biomonitoring of arsenic exposure with emphasis on the interpretation of urinary biomonitoring data. DAP is a Professor of Environmental Geochemistry at the School of Earth, Atmospheric and Environmental Sciences and Williamson Research Centre for Molecular Environmental Science at the University of Manchester. DAP’s current research interests include human biomonitoring of exposure and subsequent genetic damage from consumption of arsenic in food and drinking water. MJW is the Head of Inorganic Geochemistry at the Centre for Environmental Geochemistry, British Geological Survey and Honorary Associate Professor at the University of Nottingham, School of Biosciences, conducting research on assessing human exposures and deficiencies to harmful elements and essential minerals, respectively. RML is an Environmental Statistician at the British Geological Survey whose interests include sampling design and environmental monitoring. CJM is a geochemist with a background in groundwater chemistry and statistical programming.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
Data reported in this manuscript from human subjects were obtained from NHANES 2009–2010 and 2011–2012, which reports obtaining suitable written informed consent. Additional data were provided for comparison purposes from a previously published study (Meharg et al. 2014, Environ Pollut, 194: 181–187) in which written informed consent was obtained from all volunteers.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Aylward LL, Hays SM, Smolders R, Koch HM, Cocker J, Jones K, et al. Sources of variability in biomarker concentrations. J Toxicol Environ Health B Crit Rev. 2014;17:45–61.View ArticleGoogle Scholar
- Middleton DRS, Watts MJ, Hamilton EM, Ander EL, Close RM, Exley KS, Crabbe H, Leonardi GS, Fletcher T and Polya D.A. 2016. Urinary arsenic profiles reveal substantial exposures to inorganic arsenic from private drinking water supplies in Cornwall, UK, Sci Rep, DOI: SREP25656.Google Scholar
- Watts M, Joy E, Young S, Broadley M, Chilimba A, Gibson R, et al. Iodine source apportionment in the Malawian diet. Sci Rep. 2015;5:15251. doi:https://doi.org/10.1038/srep15251.View ArticleGoogle Scholar
- Cone EJ, Caplan YH, Moser F, Robert T, Shelby MK, Black DL. Normalization of urinary drug concentrations with specific gravity and creatinine. J Anal Toxicol. 2009;33:1–7.View ArticleGoogle Scholar
- Moeller KE, Lee KC, Kissack JC. Urine drug screening: Practical guide for clinicians. In: Proceedings of the Mayo Clinic Proceedings. Amsterdam, Netherlands: 2008;83: 66–76.Google Scholar
- Sobus JR, DeWoskin RS, Tan Y-M, Pleil JD, Phillips MB, George BJ, et al. Uses of NHANES biomarker data for chemical risk assessment: Trends, challenges, and opportunities. Environ Health Perspect. 2015;123:919–27.View ArticleGoogle Scholar
- Barrett JR. Urinary biomarkers as exposure surrogates: Controlling for possible bias. Environ Health Perspect. 2015;123:A97.View ArticleGoogle Scholar
- Yeh HC, Lin YS, Kuo CC, Weidemann D, Weaver V, Fadrowski J, et al. Urine osmolality in the US population: Implications for environmental biomonitoring. Environ Res. 2015;136:482–90.View ArticleGoogle Scholar
- Rivera-Núñez Z, Meliker JR, Linder AM, Nriagu JO. Reliability of spot urine samples in assessing arsenic exposure. Int J Hyg Environ Health. 2010;213:259–64.View ArticleGoogle Scholar
- Barr DB, Wilder LC, Caudill SP, Gonzalez AJ, Needham LL, Pirkle JL. Urinary creatinine concentrations in the US population: Implications for urinary biologic monitoring measurements. Environ Health Perspect. 2005;113:192.View ArticleGoogle Scholar
- Mayersohn M, Conrad KA, Achari R. The influence of a cooked meat meal on creatinine plasma concentration and creatinine clearance. Br J Clin Pharmacol. 1983;15:227–30.View ArticleGoogle Scholar
- Baxmann AC, Ahmed MS, Marques NC, Menon VB, Pereira AB, Kirsztajn GM, et al. Influence of muscle mass and physical activity on serum and urinary creatinine and serum cystatin C. Clin J Am Soc Nephrol. 2008;3:348–54.View ArticleGoogle Scholar
- Nermell B, Lindberg AL, Rahman M, Berglund M, Åke Persson L, El Arifeen S, et al. Urinary arsenic concentration adjustment factors and malnutrition. Environ Res. 2008;106:212–8.View ArticleGoogle Scholar
- Barber T, Wallis G. Correction of urinary mercury concentration by specific gravity, osmolality, and creatinine. J Occup Environ Med. 1986;28:354–9.Google Scholar
- Leech S, Penney M. Correlation of specific gravity and osmolality of urine in neonates and adults. Arch Dis Child. 1987;62:671–3.View ArticleGoogle Scholar
- Dossin O, Germain C, Braun JP. Comparison of the techniques of evaluation of urine dilution/concentration in the dog. J Vet Med A Physiol Pathol Clin Med. 2003;50:322–5.View ArticleGoogle Scholar
- Imran S, Eva G, Christopher S, Flynn E, Henner D. Is specific gravity a good estimate of urine osmolality? J Clin Lab Anal. 2010;24:426–30.View ArticleGoogle Scholar
- Sauvé JF, Lévesque M, Huard M, Drolet D, Lavoué J, Tardif R, et al. Creatinine and specific gravity normalization in biological monitoring of occupational exposures. J Occup Environ Hyg. 2015;12:123–9.View ArticleGoogle Scholar
- Hays SM, Aylward LL, Blount BC. Variation in urinary flow rates according to demographic characteristics and body mass index in NHANES: potential confounding of associations between health outcomes and urinary biomarker concentrations. Environ Health Perspect. 2015;123:293–300.Google Scholar
- Araki S, Murata K, Aono H, Yanagihara S, Niinuma Y, Yamamoto R, et al. Comparison of the effects of urinary flow on adjusted and non‐adjusted excretion of heavy metals and organic substances in ‘healthy’ men. J Appl Toxicol. 1986;6:245–51.View ArticleGoogle Scholar
- Araki S, Sata F, Murata K. Adjustment for urinary flow rate: an improved approach to biological monitoring. Int Arch Occup Environ Health. 1990;62:471–7.View ArticleGoogle Scholar
- Sorahan T, Pang D, Esmen N, Sadhra S. Urinary concentrations of toxic substances: an assessment of alternative approaches to adjusting for specific gravity. J Occup Environ Hyg. 2008;5:721–3.View ArticleGoogle Scholar
- Vij HS, Howell S. Improving the specific gravity adjustment method for assessing urinary concentrations of toxic substances. Am Ind Hyg Assoc J. 1998;59:375–80.View ArticleGoogle Scholar
- CDC (Centers for Disease Control and Prevention). 2013. National Biomonitoring Program Lead Biomonitoring Summary Webpage. Available: http://www.cdc.gov/biomonitoring/Lead_BiomonitoringSummary.html [Accessed 05 Sep 2015].
- Sata F, Araki S, Yokoyama K, Murata K. Adjustment of creatinine-adjusted values in urine to urinary flow rate: a study of eleven heavy metals and organic substances. Int Arch Occup Environ Health. 1995;68:64-68.Google Scholar
- CDC (Centers for Disease Control and Prevention). 2013. National Biomonitoring Program Cadmium Biomonitoring Summary Webpage. Available: http://www.cdc.gov/biomonitoring/Cadmium_BiomonitoringSummary.html [Accessed 05 Sep 2015].
- Akerstrom M, Barregard L, Lundh T, Sallsten G. The relationship between cadmium in kidney and cadmium in urine and blood in an environmentally exposed population. Toxicol Appl Pharmacol. 2013;268:286–93.View ArticleGoogle Scholar
- Middleton D, Watts M, Hamilton E, Fletcher T, Leonardi G, Close R, et al. Prolonged exposure to arsenic in UK private water supplies: Toenail, hair and drinking water concentrations. Env Sci Process Impact. 2016. doi:https://doi.org/10.1039/C6EM00072J.Google Scholar
- CDC (Centers for Disease Control and Prevention). 2015. National Health and Nutrition Examination Survey Homepage. Available: http://www.cdc.gov/nchs/nhanes.htm [Accessed 01 Oct 2015].
- Warnes GR. 2014. Sasxport: Read and write sas xport files. https://cran.r-project.org/web/packages/SASxport/index.html.
- Dragulescu AA. 2014. Xlsx: Read, write, format excel 2007 and excel 97/2000/xp/2003 files. https://cran.r-project.org/web/packages/xlsx/index.html.
- Levey AS, Stevens LA, Schmid CH, Zhang YL, Castro AF, Feldman HI, et al. A new equation to estimate glomerular filtration rate. Ann Intern Med. 2009;150:604–12.View ArticleGoogle Scholar
- Schwartz GJ, Muñoz A, Schneider MF, Mak RH, Kaskel F, Warady BA, et al. New equations to estimate GFR in children with CKD. J Am Soc Nephrol. 2009;20:629–37.View ArticleGoogle Scholar
- Meharg A, Williams P, Deacon C, Norton G, Hossain M, Louhing D, et al. Urinary excretion of arsenic following rice consumption. Environ Pollut. 2014;194:181–7.View ArticleGoogle Scholar
- R Core Team. 2013. R: A language and environment for statistical computing, r foundation for statistical computing, Vienna, Austria, www.R-project.Org.
- Kuhn M. 2015. Caret: Classification and regression training. https://cran.r-project.org/web/packages/caret/index.html.
- Dobbin KK, Simon RM. Optimally splitting cases for training and testing high dimensional classifiers. BMC Med Genomics. 2011;4:1.View ArticleGoogle Scholar
- Levine L, Fahy JP. Evaluation of urinary lead concentrations. 1. The significance of the specific gravity. J Ind Hyg Toxicol. 1945;27:217–23.Google Scholar
- Wilcox RR, Tian T. Comparing dependent correlations. J Gen Psychol. 2008;135:105–12.View ArticleGoogle Scholar
- Revelle W. 2014. Psych: Procedures for psychological, psychometric, and personality research, Northwestern University, Evanston, Illinois, https://cran.r-project.org/web/packages/psych/index.html.
- Venables WNR, Ripley BD. Modern applied statistics with S. 4th ed. New York: Springer; 2002. ISBN 0-387-95457-0.View ArticleGoogle Scholar
- Sata F, Araki S. Adjustment of creatinine-adjusted value to urine flow rate in lead workers. Archives of Environmental Health: An International; 1996.Google Scholar
- Hays S, Becker R, Leung H, Aylward L, Pyatt D. Biomonitoring equivalents: a screening approach for interpreting biomonitoring results from a public health risk perspective. Regul Toxicol Pharmacol. 2007;47:96–109.View ArticleGoogle Scholar
- Navas-Acien A, Silbergeld EK, Pastor-Barriuso R, Guallar E. Arsenic exposure and prevalence of type 2 diabetes in US adults. JAMA. 2008;300:814–22.View ArticleGoogle Scholar
- Tellez-Plaza M, Navas-Acien A, Crainiceanu CM, Guallar E. 2008. Cadmium exposure and hypertension in the 1999-2004 national health and nutrition examination survey (NHANES). Environ Health Perspect: 116:51-56.Google Scholar
- Haddow JE, McClain MR, Palomaki GE, Hollowell JG. Urine iodine measurements, creatinine adjustment, and thyroid deficiency in an adult United States population. J Clin Endocrinol Metab. 2007;92:1019–22.View ArticleGoogle Scholar