1. Introduction
Endocrine disrupting chemicals (EDCs) are a heterogeneous group of exogenous compounds that can restrict with several facets of endogenous hormones and accelerate aging, immune, metabolic or neurobehavioral disorders to threaten quality of life [
1]. Persistent exposure to EDC can disrupt homeostasis in the body and creates oxidative stress that can lead to aging and chronic inflammation [
2]. These characteristics were also found to be significant in the observation of telomere length, which is a measure of aging [
3]. As aging is a complex process influenced by genetic, environmental, and lifestyle factors, changes in endogenous hormone levels are part of this intricate interplay.
Tobacco smoke is an environmental mixture with over 7,000 chemicals including polycyclic aromatic hydrocarbons (PAHs) [
4] and has ability to obstruct and interfere in the function of endocrine system and has been entitled as EDCs [
5]. In addition to oxidative stress or inflammasome activation by tobacco smoking, a chronic activation of aryl hydrocarbon receptor signaling can contribute to premature aging and the development of neoplasms by affecting metabolism, extracellular matrix remodeling, inflammation, pigmentation, DNA repair, and apoptosis [
6]. Therefore, there is growing concern about the potential impact of tobacco smoking as an EDC, while research in this field is ongoing.
New-generation tobacco nicotine products are intended to replace combustible cigarettes to avoid tobacco smell and to consider health benefit. However, the health benefit of the new types of tobacco is still unclear [
7] and only combustible cigarettes showed product use dependent exposure to PAHs [
8], which are regarded as endocrine disruptors [
9]. Therefore, traditional combustible cigarettes can be one kind of tobacco smoking to monitor the proper exposure to PAHs.
In addition, exposure timing and early life environment for tobacco smoking have been considered, because the response to the toxic chemicals in children or youth can be manifested into adulthood via cellular memory, such as epigenetic changes [
10]. Therefore, the tobacco starting age can be the first issue for susceptibility differences. For other susceptibility biomarkers of tobacco smoking, various genome-wide association studies on tobacco smoking behavior have identified multiple loci associated with tobacco addiction [
11,
12] For example, genetic polymorphisms of cytochrome P450 (CYP) 2A6 and nicotinic acetylcholine receptor (nAChRs) are notable for tobacco-addiction [
13,
14]. We previously reported differential gene expression patterns between smokers and nonsmokers in a Korean pilot study [
15].
For exposure biomarkers of tobacco smoking, parent chemicals in tobacco smoke or their specific metabolites have been used [
16]. The most representative biomarkers are cotinine, a metabolite of nicotine; t,t-muconic acid (TTMA), a metabolite of benzene, which are quite attributable to their specificity and sensitivity among the various urinary metabolites of tobacco [
17]. In addition, altered biological responses to tobacco smoke exposure have potential as early diagnostic biomarkers. For example, alteration in genomic and mitochondrial (mt) DNA or oxidative stress parameters, such as malondialdehyde (MDA), can be considered as popular biomarkers as responsive to tobacco smoking [
15].
Therefore, we performed a molecular epidemiological approach with in-depth biological monitoring of combustible cigarette smoking among adolescents and adults in South Korea, where the current smoking rate of adolescents is 6.7% in 2019 [
18]and for adults is 20.6 % in 2020 [
19], with various exposure, susceptibility and response biomarkers to clarify how tobacco smoking differently affects adolescents and adults. This approach may reveal the underlying mechanisms of tobacco smoking as an EDC and provide preventive ways from EDCs.
3. Discussion
A molecular epidemiological approach with in-depth biological monitoring can provide reliable mechanisms and health end points of tobacco smoking as an EDC, including high susceptible population for precision prevention. For example of biomonitoring, tobacco smoking and the levels of cotinine showed association with a pronounced (~50%) reduction in fecundability, resulting in a longer time-to-pregnancy [
20]Because their developmental processes amplify or magnify the toxic responses via epigenetic or cell memory [
21]. From the self-reported cigarette pack-year, the present adolescents consumed over 1/6 volume of that of adults. However, urinary levels of cotinine in youth were 1/3 of adults and urinary MDA or TTMA of youth were similar to the adults. Although cotinine and TTMA are metabolites of tobacco products, MDA is an exposure and response biomarker for tobacco smoking via oxidative stress from tobacco chemicals or their combustion. If we assume that youth’s answers for tobacco consumption were not faked, youth can be highly susceptible to the bio-produced tobacco metabolites with small amounts of tobacco, compared to the adults. Around the onset of puberty, the activities of most of metabolic enzymes begin a gradual decline that continues throughout adolescence and concludes with attainment of adult capacity at the completion of pubertal development [
22]. In addition, a current human liver study showed that the activities of CYP2A6 and CYP2E1, major metabolic enzymes of nicotine and benzene, were somewhat lower in adults (≤69 years) than young people (21-45 years) [
23]. Thus, the same exposure to tobacco smoking can result in high levels of metabolites of tobacco components among adolescents than adults and they can be more bioactive or toxic than the parent chemicals in tobacco.
Particularly, oxidative stress and aging can be a main mechanism and a health endpoint of tobacco smoking, respectively [
24]. Chronic inhalation of cigarette smoke is a prominent cause of chronic obstructive pulmonary disease (COPD) and provides an important source of exogenous oxidants [
25]. In the present study, we found urinary MDA, a biomarker for oxidative stress or reactive oxygen species (ROS), was associated with tobacco smoking as urinary cotinine (
Table 2;
Figure 1). In addition, the growth rate of youth was negatively related to urinary MDA levels (
Figure 2). To confirm the quality of the self-report, we recalled the questionnaire and confirmed the homogeneity of the answer for the growth rate with height and body weight. Thus, the present biological monitoring suggests tobacco reduces growth rate via oxidative stress in youth. As excessive ROS might react with nucleic acid, lipids, carbohydrates, and protein causing inflammation and oxidative stress that are the main causes for the development of various metabolic disorders [
5]. In addition, the exposure to ROS from tobacco smoking can cause a higher rate of mutations in the mitochondrial genome that accumulate over time and reduce the efficiency of mtDNA repair systems [
26,
27]. The present study also shows that youth smokers, who reached to approx. 50 % of the smokers, showed high frequency of ‘SNPs for inconsistent bases between buccal and blood cells (
Table 3). Thus, high susceptibility in youth to tobacco smoking was confirmed with mtDNA alteration and MDA-related growth relay.
Finally, the pathway analyses of exposure-related gene–gene interaction suggest that MDA-related genes,
MGMT, XRCC, and
TP53, were known to interact for carcinogenesis, a degenerative disease (
Figure 3C). EDCs including tobacco smoking can increase an overall risk of ovarian aging, leading to the diminish of ovarian reserve, decline of fertility or fecundity, irregularity of the menstrual cycle and an earlier age at menopause, and/or premature ovarian insufficiency/failure in epidemiological studies [
28,
29]. We previously found tobacco smoking up-regulated aging genes, such as
DEFA4 for hearing loss in adults [
15]. In addition, others reported that tobacco smoking advanced brain aging in middle age women [
30]. Moreover, tobacco smoking-related oxidative stress has been emphasized in the appearance of the clinical manifestation of skin aging [
30,
31]. Hexane-soluble tobacco smoke extract may induce
matrix metalloproteinase-
1 expression in human skin fibroblasts through the activation of the aryl hydrocarbon receptor pathway, which is pathogenetically involved in extrinsic skin aging [
32]. Thus, the present susceptibility/genetic biomarkers support that oxidative stress and aging are ones of the mechanisms or health end points of tobacco smoking as an EDC.
Various kinds of tobacco attract youth and it can be a potential risk for their various future diseases [
18,
33] and as a reaction to that, we also consider effective tobacco control, such as raising the minimum legal age for access to tobacco [
34]. The present study also can be a fresh approach to provide scientific evidences to protect adolescents from tobacco.
In conclusion, the present biological monitoring provides high susceptible population, such as adolescents rather than adults, and reliable genetic factors affecting tobacco exposure. The inferred environment-gene-gene interaction suggests accelerating aging in adolescents by tobacco smoking as an EDC.
4. Materials and Methods
4.1. Subjects and Sampling
We recruited adult volunteers, who visited Eulji University Hospital in Daejeon, South Korea for regular examination, and adolescents from high schools around Geumsan-gun near Daejeon. All subjects provided written informed consents and completed extensive questionnaires including medical and smoking history, dietary patterns, alcohol drinking, environment of residency, and smoking behaviors, such as smoking cessation, cigarettes smoker per day (CPD), duration of smoking, and smoking initiation. In addition, the Fagerstrom test of nicotine dependence (FTND) was used to assess nicotine dependence [
35]. None of these subjects had any history of pulmonary, cardiovascular, endocrine, or gastrointestinal disorders. In addition, we measured carbon monoxide (CO) during exhalation with Micro CO Monitor (On-site Lab, Seoul, Korea).
Peripheral blood samples (10 ml) were collected into evacuated tubes containing sodium heparin as an anticoagulant (BD Vacutainer, Franklin Lakes, NJ, USA). In addition, spot urine specimens, the first voids of urine (40 ml) before breakfast, were collected into 50 ml of conical tubes. Both urine and blood samples were stored at -20℃ until analyses. All study protocols for this study were approved by the Institutional Review Board of Eulji University Hospital (Daejeon, Korea).
To compare mtDNA alteration in buccal cells to that in blood cells, we also collected buccal cells from the subjects with sterile cotton swabs, following our previous method [
15].
4.2. Analyses of Urinary Cotinine
We analyzed urinary cotinine by our previous ion-pair HPLC/UVD method [
15] with minor modifications. In brief, 900 μl of each urine sample was mixed with 100 μl of
80 μM 2-phenylimidazole as an internal standard and 330 μl of 3 M NaOH. The mixture was twice extracted with 3 ml of CH2Cl2 each time. After evaporating CH2Cl2 -extract, we dissolved the residue in 1 ml of water and injected 20 μl of its supernatant fraction to HPLC. The HPLC system consisted of dual Younglin SP930D pumps (Younglin, Seoul, Korea), an MIDAS COOL autosampler (Spark Holland, Emme, The Netherlands), an SPD-10A UV-VIS detector (Shimadzu, Kyoto, Japan), and a TSK gel ODS-80™ column (5 μm, 4.6 mm × 150 mm, Toyo Soda Co., Tokyo, Japan). Analyses were carried out with the following gradient mode: mobile phase A, a mixture of acetonitrile/water (15/85) containing 20 mM KH2PO4 and 3 mM sodium 1-octanesulfonate (pH 4.5); B, methanol; Flow rate, 0.7 ml/min; 0–20 min, ratio of A to B = 100:0; 20–25 min, ratio of A to B = 100:0 to 50:50; 25–30 min, ratio of A to B = 50:50; 30–35 min, ratio of A to B = 50:50 to 100:0; and 35-45 min, ratio of A to B = 100:0. The column was kept at 50 ℃ and the absorbance was observed at wavelength of 254 nm.
4.3. Analyses of Urinary MDA
We quantified urinary MDA as adducts of 2-thiobarbituric acid (TBA, CAS number: 504-17-6) with HPLC/UVD [
36]. TBA-MDA adducts were detected at 532 nm with isocratic mode. The mobile phase was a mixture of 50 mM potassium phosphate buffer (pH 6.8) and methanol (58:42, v/v). Flow rate was set at 0.6 ml/min.
4.4. Analyses of Urinary TTMA
We analyzed urinary TTMA with UPLC-MS/MS, using the method previously described by Gagne et al. [
37] with a minor modification. Briefly, TTMA standard and deuterated internal standard, d
4-TTMA, were obtained from Sigma-Aldrich and CDN Isotopes Inc. (Pointe-Claire, Quebec, Canada), respectively. Standard TTMA solutions (0.025-2.5 ng/ml) were prepared in 50% methanol. In short, 950 μl of 0.1 % formic acid containing 2.38 μg/ml of d
4- TTMA was mixed with 50 μl of TTMA standards or urine samples. After centrifugation at 13,000 rpm for 10 min, 5 μl of each supernatant was analyzed with UPLC-MS/MS. The UPLC-MS/MS system consisted of a Waters Acquity UPLC coupled with a Waters Xevo TQ triple quadrupole mass spectrometer (Beverly, MA), and an Acquity UPLC BEH C18 (1.7 μm, 2.1 mm × 50 mm, Waters). Mobile phases were composed of 0.1% formic acid in methanol (eluant A) and in water (eluant B). UPLC separation was achieved with a gradient from 10 to 95 % of eluant A for 1.25 min. Eluant A composition was then held constant for 0.5 min followed by a 0.5 min equilibration period at 10 % of eluant A. The flow rate was set at 0.5 mL/min and the column temperature was kept at 50ºC. The Xevo TQ was operated in negative mode. The capillary voltage was set at 2.8 kV. The source temperature was at 150 ºC. The desolvation temperature was at 500 ºC. Desolvation flow rate was at 900 L/hr and collision gas flow rate was at 0.15 mL/min. Data were acquired in multiple reaction monitoring (MRM) mode.
Urinary cotinine, MDA, and TTMA were adjusted for creatinine, measured with ion-pair HPLC/UVD method [
15].
4.5. Targeted Genotyping
Genomic DNA of peripheral blood was isolated with a QIAamp DNA Blood Mini kit (Qiagen, Hilden, Germany) according to the manufacturer’s instructions. We measured the purity and concentration of the isolated genomic DNA (gDNA) using a NanoDrop® ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE). For the genotyping, we used gDNA samples with a 260:280 ratio of 1.5 or higher and concentrations of 60 ng/µl or higher.
We selected the 96 target SNPs, based on the SNPs of tobacco smoking-responsive genes, of which expression levels were altered by tobacco smoking in our previous microarray study [
15], i.e.,
ACTG1, DEFA4, VAV3,
FCGR3A, etc., as well as the SNPs related to metabolism, e.g.
, CYP2A6, CYP1B1, CYP2E1, and
NQO1, addiction [
38], e.g.,
CHRNA5/A3/B4, 5-HTTLPR, and
DRD2 [
39]
, risks for lung cancer, e.g.,
CYP1A1, GSTP1, and
MPO1 [
39,
40], DNA repair, e.g.,
ERCC1, MGMT, XRCC1, etc. [
41], and epigenetic modulation, e.g.,
HDAC1 and
MTHFR [
33]
, of tobacco smoking (supplement
Table 1). For genotyping assays, the 96.96 Dynamic Array™ integrated fluidic circuits (Fluidigm Corp., South San Francisco, CA) was used. Prior to genotyping, specific target amplification (STA) was performed for each gDNA to enrich targeted SNP sequences. Briefly, 70 ng of gDNA was mixed with 50 nM of STA primer mixture, 50 nM of locus specific primer mixture, and 2.5 µL of Qiagen 2X Multiplex PCR Master Mix (Qiagen) in a final volume of 5µL. PCR reactions were performed on an Arktik Thermal Cycler (Thermo Scientific, Rockford, IL) with the following cycling conditions: 10 min at 95°C, followed by 14 cycles of a 2-step amplification profile of 15 sec at 95°C and 4 min at 60°C. STA products were 100-fold diluted in DNA suspension buffer and 2.5 µL of each product was combined with 3.0 µL of 2X Maxima
® Probe/ROX qPCR Master Mix (Fermentas, St. Leon-Rot, Germany), 0.3 µL of SNPtype 20 X Sample Loading Reagent (Fluidigm), 0.1 µL of SNPtype Reagent (Fluidigm), and 0.1 µL of nuclease-free water. In parallel, 1 µL of each SNPtype assay was mixed with 2.5 µL of 2 X Assay Loading Reagent (Fluidigm) and 1.5 µL of nuclease-free water and loaded into the Fluidigm 96.96 Dynamic Genotyping Arrays. PCR and image processing were carried out on an EP1 system (Fluidigm). We analyzed data with an automated genotype calling algorithm using Fluidigm SNP Genotyping Software (v3.1.1).
4.6. Analyses of mtDNA Alteration between Blood and Buccal Cells
After extraction of DNAs from buccal cells and blood with DNeasy blood & tissue kit (Qiagen), we prepared a 96 well plate (Bioline, London, UK) to load 50 ul of DNA samples (5 ng/µl =250 ng). From the master plates, 5ng of DNA was used for a conventional PCR to amplify HV2 region. Total reaction volume was 20 µl, containing 5ng DNA and 75 µM of forward and backward primers, F015 5’-CAC CCT ATT AAC CAC TCA CG-3’ and R569 5’-GGT GTC TTT GGG GTT TGG TTG-3’, respectively. PCR conditions for amplification were performed on Primus 96 plus (MWG_Biotech, Huntsville, AL): initial denaturation at 96°C for 5 min; 35 cycles of denaturation at 96°C for 30s; annealing at 56°C for 30 s; elongation at 72°C for 1 min, with MyGenieTM 96 Gradient Thermal Block (Bioneer, Daejeon, Korea). For sequencing target mtDNA, we used BigDye Terminator v3.1 sequencing kit and ABI3730XL (Applied Biosystems, Waltham, MA ).
For data analyses, we used DNASTAR Lasergene SeqMan Pro version 7.1.0. Among the reference sequence of the mtDNA HV2 region (015-560 bp), NCBI (ref|NC_012920.1) corresponding location from100 bp to 322 bp of mtDNA HV2 region was analyzed for detecting discrepancies between blood and oral DNA on (a/g or t/c) polymorphism at the position 263 bp and polyC region (303~314 in reference C7(T)C5).
4.7. Statistical Analyses
We removed imputed SNPs with < 0.05 genotype information content, low call rates with < 0.90, and minor allele frequency (MAF) with < 0.05. The Shapiro-Wilk W test was used to test distributional normality for levels of exposure biomarkers (i.e., urinary cotinine, TTMA and MDA). ANOVA or Wilcoxon rank sum test was used to analyze differences in characteristics among never, ex-, and current smokers. Associations between SNPs and levels of exposure biomarkers were investigated with Wilcoxon Mann-Whitney Test or multi linear regression models after adjusting for age, body mass index (BMI), and alcohol consumption under a dominant model. All statistical analyses were conducted using JMP package v. 4.0.2 (SAS Institute, Cary, NC)