ABSTRACT
Background The nicotine metabolite ratio and nicotine equivalents are measures of metabolism rate and intake. Genome-wide prediction of these nicotine biomarkers will extend biomarker studies to cohorts without measured biomarkers and enable tobacco-related behavioral and exposure research.
Methods We screened genetic variants genome-wide using marginal scans and applied statistical learning algorithms on top-ranked genetic variants and age, ethnicity and sex, and cigarettes per day (CPD) (in additional modeling) to build prediction models for the urinary nicotine metabolite ratio (uNMR) and creatinine-standardized total nicotine equivalents (TNE) in 2,239 current cigarette smokers in five ethnic groups. We predicted these nicotine biomarkers using model ensembles, and evaluated external validity using behavioral outcomes in 1,864 treatment-seeking smokers in two ethnic groups.
Results The genomic regions with the most selected and trained variants for measured biomarkers were chr19q13.2 (uNMR, without and with CPD) and chr15q25.1 and chr10q25.3 (TNE, without and with CPD). We observed ensemble correlations between measured and predicted biomarker values for the uNMR and TNE without (with CPD) of 0.67 (0.68), and 0.65 (0.72) in the training sample. We observed inconsistency in penalized regression models of TNE (with CPD) with fewer variants at chr15q25.1 selected and trained. In treatment-seeking smokers, predicted uNMR (without CPD) was significantly associated with CPD, and predicted TNE (without CPD) with CPD, Time-To-First-Cigarette, and Fagerström total score.
Conclusions Nicotine metabolites, genome-wide data and statistical learning approaches develop novel robust predictive models for urinary nicotine biomarkers in multiple ethnic groups. Predicted biomarker associations help define genetically-influenced components of nicotine dependence.
IMPLICATIONS We demonstrate development of robust models and multiethnic prediction of the urinary nicotine metabolite ratio and total nicotine equivalents using statistical and machine learning approaches. Trained variants in models for both biomarkers include top-ranked variants in multiethnic genome-wide studies of smoking behavior, nicotine metabolites and related disease. Association of the two predicted nicotine biomarkers with Fagerstr□m Test for Nicotine Dependence items support models of nicotine biomarkers as predictors of physical dependence and nicotine exposure. Predicted nicotine biomarkers may facilitate tobacco-related disease and treatment research in samples with genomic data and limited nicotine metabolite or tobacco exposure data.
Competing Interest Statement
AWB is an employee of Oregon Research Institute and Oregon Community and Evaluation Services, and serves as a Scientific Advisor and Consultant to BioRealm, LLC. CME is a co-owner and the Principal Biostatistician for BioRealm, LLC. HAT has served as PI on NIH-supported studies for smoking cessation in which the medication was donated by the manufacturer (e.g., Pfizer, varenicline). CSM, SM, LLM, DOM, SEM, YMP and SLP have no conflicts of interest to report. JWB is an employee and an owner of BioRealm, LLC. BioRealm, LLC offers services related to the Smokescreen Genotyping Array and analysis of nicotine biomarkers.
Funding Statement
FUNDING This work was supported by the National Institute on Alcohol Abuse and Alcoholism (R44 AA027675 to AWB, CSM, SM, CME, SLP, and JWB) and by the National Cancer Institute (R01 CA232516 to HAT, and U01 CA164973 and P01 CA138338 to LLM, DOS, SEM, YMP and SLP). The sponsors had no role in the analysis of data, writing of the report, or in the decision to submit the paper for publication. ACKNOWLEDGMENTS The authors thank participants of the Multiethnic Cohort Study and of the University of Wisconsin cessation trials. The authors acknowledge the contribution of data from Genetic Architecture of Smoking and Smoking Cessation accessed through dbGAP (phs000404.v1.p1). Support for genotyping, which was performed at the Center for Inherited Disease Research (CIDR), was provided by 1 X01 HG005274-01. CIDR is fully funded through a federal contract from the National Institutes of Health to The Johns Hopkins University, contract number HHSN268200782096C. Assistance with genotype cleaning, as well as with general study coordination, was provided by the Gene Environment Association Studies (GENEVA) Coordinating Center (U01 HG004446). Funding support for collection of the University of Wisconsin Transdisciplinary Tobacco Use Research Center cessation trials was provided by P50 DA019706 and P50 CA084724.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
Written informed consent was obtained from all participants. The research described herein received approvals from the Institutional Review Boards of BioRealm, the Oregon Research Institute, the University of Hawaii and the NIH Joint Addiction, Aging, and Mental Health Data Access Committee.
All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Data Availability
The Multiethnic Cohort data is available upon application to the Multiethnic Cohort (https://www.uhcancercenter.org/for-researchers/mec-data-sharing). The University of Wisconsin data is available upon application to the Database of Genotypes and Phenotypes (https://www.ncbi.nlm.nih.gov/projects/gap/cgi-bin/study.cgi?study_id=phs000404.v1.p1).
