Calibration of Toenail Metal Concentrations for Sample Mass Heterogeneity and Between-Batch Variability: The COMET Approach

Abstract

Background: Toenails are promising biomarkers of long-term metal exposure in epidemiological studies, but their accuracy may be compromised by systematic and random errors associated with heterogeneous toenail sample masses, as well as by substantial variability across laboratory batches. Objectives: We propose a novel modeling approach to calibrate toenail metal concentrations for the heterogeneity in sample masses and the variability between batches. Methods: We developed a heteroscedastic spline mixed model relating sample mass and laboratory batch with measured concentrations, allowing for an average bias in measurements over all batches as a smooth function of sample mass, random variation in mass-related biases across batches, and mass-related heterogeneity in within-batch error variance. The model allowed partitioning the total variance of measured concentrations into the extraneous variances (due to different sample masses and laboratory batches) and the intrinsic variance (resulting from distinct metal exposures). We derived calibrated metal concentrations from the model by removing both sources of extraneous variation and estimating the predicted concentrations had all toenail samples been analyzed in a single batch and of the same mass. We provide the R script COMET (COrrected METals) to fit the proposed model, extract variance components, and calibrate metal concentrations. Results: In a multicase-control study in Spain (MCC-Spain) with toenail determinations for 16 metals in 4,473 incident cases of five common cancers and 3,450 population controls, sample mass and batch accounted for 26%-60% of the total variance of measured concentrations for most metals. In comparison with calibrated concentrations, odds ratios for measured concentrations were biased by > 10% toward or away from the null in one-quarter of the estimated metal-cancer associations. Discussion: The proposed model allows correcting toenail metal concentrations for sample mass heterogeneity and between-batch variability and could be applied to other biological specimens of heterogeneous size, distinct laboratory techniques, and different study designs.