Printer Friendly

Genetic and non-genetic predictors of LINE-1 methylation in leukocyte DNA.

DNA methylation plays a fundamental role in regulation of gene expression, genomic imprinting, X-chromosome inactivation, and repression of transposable elements (Jones and Liang 2009). Aberrant DNA methylation has been associated with various cancers and with developmental, autoimmune, and other chronic diseases (Robertson 2005). Global DNA methylation can be directly quantified by measuring 5-methylcytosine content of the genome, or it can be estimated based on methylation of repetitive sequences such as Alu elements or long interspersed nuclear element 1 (LINE-1) (Yang et al. 2004). Age, sex, smoking, and arsenic and lead exposures have been associated with DNA methylation, but findings have been inconsistent among studies (Breitling et al. 2011; El-Maarri et al. 2007; Fraga et al. 2005; Fuke et al. 2004; Terry et al. 2011). The folate and methionine-dependent one-carbon metabolism pathway could modulate DNA methylation by altering the level of 5-adenosylmethionine (SAM), the principal source of methyl groups (Ulrey et al. 2005). Genetic variants might also influence the methylation of CpG loci locally, or might have a global influence on methylation throughout the genome. For example, a single nucleotide polymorphism (SNP) in TRPC3 (transient receptor potential cation channel, subfamily C, member 3)-isoform 2 has been reported to regulate the methylation status of its own promoter (Martin-Trujillo et al. 2011), and variants of the methylene-tetrahydrofolate reductase gene (MTHFR) have been associated with global DNA hypomethylation (Castro et al. 2004; Friso et al. 2002). However, although the determinants of global and site-specific methylation are widely assumed to be likely contributors to health and disease, they are poorly defined at this time.

Assessing the impact of both genetic and non-genetic factors on global DNA methylation may improve our understanding of the molecular pathogenesis of many common diseases. Therefore, we investigated associations of global DNA methylation in LINE-1 from bisulfite-modified granulocyte DNA with genetic variants and personal, demographic, lifestyle, and environmental characteristics.

Methods

Study population. The study population, design, and data collection have been previously described (Garcia-Closas M et al. 2005). Briefly, participating individuals were controls from the Spanish Bladder Cancer/EPICURO study who were admitted to hospitals in five regions of Spain for a range of conditions including hernia, fractures, and other non-cancer diseases, and were 20-81 years of age. We collected demographic and exposure information at the hospitals using computer-assisted personal interviews. Of a total of 1,271 controls who agreed to participate in the study and were interviewed, 1,056 provided blood for DNA extraction. We excluded 23 subjects because of inadequate or poor quality DNA (n = 15) or missing smoking status data (n = 8). Three subjects with missing data on smoking status were included because they had data on other variables including age, sex, region, and body mass index (BMI). To ensure homogeneity, we also excluded one non-Caucasian individual, leaving 925 individuals with granulocyte DNA for bisulfite modification and pyrosequencing. Pyrosequencing failed in 33 individuals; thus, the final study population for the present analysis included 892 participants. We obtained written informed consent from all participants, and the study was approved by the local Spanish institutional review boards and U.S. National Cancer Institute.

Quantification of LINE-1 methylation. We extracted granulocyte DNA using standard methods (Garcia-Closas M et al. 2005). We carried out bisulfite conversion of DNA using the EZ-96 DNA Methylation-Gold[TM] kit (Zymo Research, Irvine, CA, USA) according to the manufacturer's recommendations. We carried out polymerase chain reaction (PCR) amplification of bisulfite-modified DNA using a set of forward and reverse primers reported previously (Estecio et al. 2007). To quantify the methylation level of each of the first four CpG sites next to the pyrosequencing primer, we performed sequencing of the PCR product by pyrosequencing, using the PyroMark[TM] Q24 System (QIAGEN, Valencia, CA, USA) as recommended by the manufacturer. The first four were the CpGs from which we could obtain methylation values of all samples. We extracted the methylation level at each CpG site using the PyroMark[TM] application software, version 2.0.6 (QIAGEN), and we expressed the value as the percentage of methylated cytosines over the sum of methylated and unmethylated cytosines. We used the average methylation level of the first four LINE-1 CpG sites as a surrogate marker of the global DNA methylation level. To determine whether changes in blood cell populations affect LINE-1 methylation levels, we analyzed LINE-1 methylation in independently purified granulocyte and lymphocyte samples. The results showed no significant difference in LINE-1 methylation between granulocyte and lymphocyte samples, thus suggesting that variation in the distribution of peripheral blood cell populations among participants would not contribute to variation in global DNA methylation (data not shown). For the present study, we used DNA extracted from granulocytes to quantify DNA methylation. As a quality control measure, we measured LINE-1 methylation in 129 randomly selected duplicate samples and the within-sample coefficient of variation (CV) was 4.0%. In the analysis, we used the average of the duplicates for those samples.

Nutritional assessment. We estimated the usual intake of vitamins B1, B2, B3, B6, and B12, folate, protein, alcohol, fruit, and vegetables over the 5 years before interview using a validated food frequency questionnaire of 127 items (Garcia-Closas R et al. 2007). Micronutrients and macronutrients included in the present analysis have been suggested as important cofactors and methyl donors in one-carbon metabolism (Stover 2009). We calculated nutrient density variables by dividing the total estimated mass of daily food consumed by the total estimated daily energy intake (micrograms per day per kilocalorie).

Trace elements. We collected toenail clippings to estimate chronic exposure to trace elements. Sample collection and experimental methods used to measure trace elements level have been reported previously (Amaral et al. 2012). Briefly, after cleaning and washing the toenails to remove external contaminants, we quantified elements at the Trace Element Analysis Core (Dartmouth College, Hanover, NH, USA), using inductively coupled plasma-mass spectrometry (Hopkins et al. 2004). We digested the samples with Optima[R] Nitric Acid (Fisher Scientific, St. Louis, MO, USA) at 105[degrees]C followed by addition of hydrogen peroxide and further heating the dilution with deionized water. We recorded gravimetrically all sample preparation steps. As a quality control, each batch of analyses included six standard reference material samples with known trace element content (GBW 07601, powdered human hair; China National Analysis Center for Iron and Steel, Beijing, China) and six analytic blanks, along with the study samples. In total, we determined concentrations (micrograms per gram) of 12 trace elements (aluminum, arsenic, cadmium, chromium, copper, iron, lead, manganese, nickel, selenium, vanadium, and zinc).

Genotyping. For genotype assays, we extracted DNA from leukocytes as described previously (Garcia-Closas M et al. 2005). We determined genotypes at the Core Genotyping Facility of the Division of Cancer Epidemiology and Genetics, National Cancer Institute, MD, USA. We selected for the analysis a total of 515 SNPs in 24 genes involved in the one-carbon metabolism pathway, including DNA methylation and arsenic metabolism [for a list of the 24 genes evaluated, see Supplemental Material, Table S1 (http://dx.doi.org/10.1289/ehp.1206068)]. We selected these genes because they are critical for the one-carbon metabolism pathway (Lee YL et al. 2009; Ulrey et al. 2005). Previously, we described in detail methods of the genotyping process (Garcia-Closas M et al. 2005, 2007; Rothman et al. 2010). We genotyped SNPs using Illumina Infinium[R] Human1M-Duo, Illumina GoldenGate[R] (Illumina, San Diego, CA, USA) and TaqMan[R] (Applied Biosystems, Foster City, CA, USA) assays (for a complete list of SNPs according to assay, see Supplemental Material, Table S2). In addition, we estimated associations between LINE-1 methylation and glutathione S-transferase mu 1 (GSTM1), glutathione S-transferase theta 1 (GSTT1), and N-acetyltransferase 2 (NAT2) variants because of their relevance to bladder cancer (Cash et al. 2012). These variants were determined as described by Garcia-Closas M et al. (2005). All genotypes included in the study were in Hardy-Weinberg equilibrium in the study population (p > 0.05) (data not shown).

Statistical analysis. The distribution of LINE-1 methylation levels was slightly bimodal and positively skewed [see Supplemental Material, Figure S1 (http://dx.doi.org/10.1289/ ehp.1206068)]. To estimate associations between LINE-1 methylation levels and each of the variables considered, we fitted bivariate robust linear regression models and calculated the corresponding beta coefficients and 95% confidence intervals (CIs). Characteristics analyzed as continuous variables were age, micronutrient intakes, fruit and vegetable intakes, and toenail concentrations of trace elements. Characteristics analyzed as categorical variables were BMI (< 25.0, 25.0-26.99, 27.0-29.99, [greater than or equal to] 30.0), smoking status (non-, occasional, former, current smoker), and tobacco type (non-smoker, blond only, black only, blond and black, unknown).

To identify SNPs for detailed assessment, we first used the Fisher's exact test to screen SNPs that were significantly associated (p < 0.05) with LINE-1 methylation categorized according to tertiles (< 56.7%, 56.7-58.6%, and > 58.6%) according to codominant mode of inheritance. The 22 SNPs identified for further analyses [listed in Supplemental Material, Table S3 (http:// dx.doi.org/10.1289/ehp.1206068)] were subsequently modeled according to all modes of inheritance (additive, codominant, dominant, and recessive). The mode of inheritance that best predicted LINE-1 methylation is reported.

In addition to age and sex, adjusted robust linear regression models for each potential predictor included region, which may be related to diet, micronutrients (Gabriel et al. 2006), and environmental pollution and smoking status, which may be related to trace elements (Moerman and Potts 2011). We also did a sensitivity analysis without adjusting for smoking status to see whether there was a change in the beta estimates. We included all the 515 SNPs in the analysis regardless of linkage disequilibrium. The association between LINE-1 methylation and potential predictors was assessed in the multivariable adjusted model stratifying by sex. Because arsenic and rs7085104 in arsenic (+3 oxidation state) methyltransferase (AS3MT), which is involved in arsenic metabolism, were individually associated with LINE-1 methylation, we assessed the presence of effect modification by including a multiplicative interaction term in a model adjusted for age, sex, region, and smoking status. The Wald test was used to calculate the interaction p-value. We corrected for multiple testing using the Bonferroni method. We conducted a sensitivity analysis, excluding 42 individual participants who had a CV > 4% for LINE-1 methylation in duplicate samples. This exclusion did not result in substantial differences in the beta coefficients; therefore, those participants remained in the analyses (data not shown). All statistical tests were two-sided, and p [less than or equal to] 0.05 was considered significant. We carried out all data analyses using Stata/ SE version 10.1 (StataCorp, College Station, TX, USA).

Results

The characteristics of the 892 participants in the present study and median and mean LINE-1 methylation levels according to each variable of interest are provided in Supplemental Material, Table S4 (http:// dx.doi.org/10.1289/ehp.1206068). The majority of the study participants were male (89%) and regular (former or current) smokers (63.9%), with a median age of 66 years [interquartile range (IQR) = 13]. The mean LINE-1 methylation level was 58.9% (SD = 5.3%) with minimum and maximum value of 37.9% and 85.7%, respectively. Table 1 shows the association between LINE-1 methylation levels and characteristics of study subjects. In the bivariable robust linear regression analysis, only toenail arsenic and nickel concentrations were significantly associated with LINE-1 methylation. However, in multivariable robust linear regression models adjusted for age, sex, region, and smoking status, the levels of LINE-1 methylation were significantly lower in women than in men (adjusted [beta] = -0.7; 95% CI: -1.2, -0.1, p = 0.02) and in smokers of blond tobacco only (adjusted [beta] = -0.7; 95% CI: -1.3, -0.08, p = 0.03) and of both blond and black tobacco (adjusted [beta] = -0.6; 95% CI: -1.1, -0.07, p = 0.03) compared with nonsmokers. Toenail arsenic concentration also was negatively associated with LINE-1 methylation (adjusted [beta] for a 1-[micro]g/g increase = -3.6; 95% CI: -5.9, -1.2, p = 0.003). In contrast, LINE-1 levels were positively associated with 1-[micro]g/g increases in toenail concentrations of iron (adjusted [beta] = 0.002; 95% CI: 0.001, 0.004, p = 0.009) and nickel (adjusted [beta] = 0.02; 95% CI: 0.005, 0.03, p = 0.004). BMI, B vitamins, folate, total protein, alcohol, and fruit and/or vegetable intake were not significantly associated with LINE-1 methylation regardless of adjustment for covariates (Table 1). Results from the multivariable analyses without adjusting for smoking status done as a sensitivity analysis were not different from the associations above (see Supplemental Material, Table S5). If we were to correct for multiple comparisons by Bonferroni method, none of the above would be significant.

Out of the 515 genetic variants assessed, 22 passed a first screening using Fisher's exact test (p [less than or equal to] 0.05 according to codominant mode of inheritance) [see Supplemental Material, Table S3 (http://dx.doi.org/10.1289/ehp.1206068)]. Of these, 7 SNPs in five genes were significantly associated with LINE-1 methylation based on multivariable models adjusted for age, sex, region, and smoking status (Table 2; model-based estimates for the 15 SNPs that were not significantly associated with LINE-1 methylation are reported in Supplemental Material, Table S6.) Significant positive associations were estimated for DNA (cytosine-5-)-methyltransferase 3 alpha (lWMT3A)-rs7581217 (per allele: adjusted [beta] = 0.3; 95% CI: 0.1, 0.6, p = 0.002); trans-cobalamin II (TCAQ)-rs9621049 (recessive: adjusted [beta] = 4.2; 95% CI: 2.8, 5.7, p = 4.7 x [10.sup.-9]), TCA2-rs4820887 (recessive: adjusted P = 4.0; 95% CI: 2.5, 5.6, p = 4.8 x [10.sup.-7]) and TCA2-rs9606756 (recessive: adjusted [beta] = 1.9; 95% CI: 0.5, 3.3, p = 0.008); and AS3MT-rs7085104 (recessive: adjusted [beta] = 0.7; 95% CI: 0.3, 1.2, p = 0.001). In addition, significant associations under the codominant mode of inheritance (based on global p-values) were estimated for solute carrier family 19 (folate transporter), member 1 (SZC75A7)-rs914238 (TC vs. TT, [beta] = 0.5; 95% CI: 0.08, 0.8; CC vs. TT, [beta] = -0.3; 95% CI: -0.7, 0.2; global p = 0.0007) and 5,10-methenyltetrahydrofolate synthetase (MTHFS)-rs1380642 (CT vs. CC, [beta] = 0.3; 95% CI: -0.08, 0.6; TT vs. CC, [beta] = -0.8; 95% CI: -1.6, 0.09; global p = 0.05). After correcting for multiple testing using the Bonferroni method, TCA2-rs9621049 and TCA2-rs4820887 remained significant p < 0.05).

A significant interaction (p = 0.01) was observed between arsenic and AS3MT-rs7085104 on LINE-1 methylation (adjusted [beta] for a 1-[micro]g/g increase in As = -4.1; 95% CI: -6.6, -1.7, p = 0.001 for genotype AA/AG; and adjusted [beta] for a 1-[micro]g/g increase in As = 10.2; 95% CI: -3.2, 23.7, p = 0.1 for genotype GG).

After simultaneously adjusting for age, geographic region, and all factors that were significant predictors of LINE-1 methylation (sex; tobacco type; toenail arsenic, iron, and nickel; and the 5 SNPs noted above), associations with sex, arsenic, nickel, iron, ZWMT3A-rs7581217, TCA2-rs9621049, and MTHFS-rs1380642 remained significant. The association with blond tobacco was nonsignificant although the direction of the point estimate remained unchanged. The association between rs9606756, rs4820887, and LINE-1 methylation become nonsignifcant [see Supplemental Material, Table S7 (http://dx.doi.org/10.1289/ehp.1206068)]. This might be due to reduced sample size in the simultaneously adjusted model due to missing data.

Discussion

In the present study, we used a comprehensive approach to assess associations of genetic and non-genetic factors with LINE-1 methylat ion in a group of participants 20-81 years of age. Lower levels of LINE-1 methylat ion were found among women compared with men, and among smokers of blond tobacco compared with nonsmokers.

In addition, toenail concentrations of arsenic were also negatively associated with LINE-1 methylation. On the other hand, LINE-1 methylation levels were positively associated with toenail concentrations of iron and nickel, and with seven variants in DNMT3A, TCN2, AS3MT, SLC19A1, and MTHFS genes.

Our findings support previous results showing that women have significantly lower levels of LINE-1 methylation (El-Maarri et al. 2007, 2011; Wilhelm et al. 2010; Zhu et al. 2012). DNA methylation is important for X-chromosome inactivation in women (Jones and Liang 2009), and although LINE-1 sequences do not seem to be the major mechanism involved in this process, they may be involved in spreading the X-inactivation signal across the chromosome (Bailey et al. 2000). In support of this, a recent study showed that LINE-1 sequences were hypomethylated in the inactive X-chromosome (Singer et al. 2012). A small study of 33 men and 33 women reported lower levels of blood SAM in women (Poirier et al. 2001). Hormonal factors may also contribute to the difference in methylation levels between sexes. However, a recent in vitro study assessing the roles of estrogen, progesterone, and dihydrotestosterone on DNA methylation in four cell lines found no detectable effect of these hormones on methylation levels at the LINE-1 and Alu repeats (El-Maarri et al. 2011). Further studies are needed to decipher the relationship between sex and LINE-1 methylation.

Because tobacco smoking is an important contributor to disease and is a modifiable behavioral factor, there has been much interest in the relationship between smoking and DNA methylation. Our findings are in line with other studies that reported no association between LINE-1 methylation and smoking status (Terry et al. 2011). In the present study, we found that subjects who smoked blond tobacco had lower levels of global DNA methylation than nonsmokers. An experimental study has shown that cigarette smoke condensates can induce DNA demethylation in repeat elements such as LINE-1 and D4Z4 (Liu et al. 2010). Both black and blond tobacco cause disease although the former is more mutagenic, reflective of the higher levels of N-nitrosamines and aromatic amines in smoke produced by black tobacco (Malaveille et al. 1989). Our findings suggest that the toxic effects of blond tobacco could be mediated by modulating the epigenetic landscape. This may have a public health implication given epigenetic alterations are reversible.

We also provide evidence that arsenic levels were inversely associated with LINE-1 methylation, and that arsenic may have a strong effect on LINE-1 methylation. For every 1-[micro]g/g increase in arsenic there was a 3.6% decrease in DNA methylation level. This inverse association is in agreement with that from a population-based study that used a similar assay to assess LINE-1 methylation levels and toenail concentrations of arsenic (Wilhelm et al. 2010), as well as with several other experimental studies (Reichard and Puga 2010; Ren et al. 2011). The mechanisms through which arsenic exposure influences DNA methylation are not fully understood. Studies in cell lines and mouse models exposed to arsenic for [less than or equal to] 22 and 48 weeks, respectively, have shown that prolonged exposure to sodium arsenite resulted in decreased global DNA methylation, and inhibition of DNA (cytosine-5-)-methyltransferase 1 (DNMT1), DNMT3A, and DNA (cytosine-5-)-methyltransferase 3 beta (DNMT3B) gene expression (Reichard and Puga 2010; Ren et al. 2011). It is likely that through the combined effect of depleting the cellular pool of SAM and inhibiting the activity of DNMTs, both inorganic and organic arsenic may decrease global DNA methylation.

We are not aware of any human studies associating iron and nickel levels and global DNA methylation. In the present study, iron and nickel showed a small but significant positive association with LINE-1 methylation level. Genes involved in hepatocellular carcinoma (HCC) have been found to be hypermethylated in hereditary hemochromatosis, a disease characterized by chronic iron overload that is a risk factor for HCC (Lehmann et al. 2007). Iron, together with 2-oxoglutarate and oxygen, is an essential cofactor for the ten-eleven translocation (TET) family of proteins that hydroxylate 5-methylcytosine to 5-hydroxymethylcytosine and further oxidize to 5-carboxylcytosine and 5-formylcytosine, which have all been suggested to be precursors for both active and passive DNA demethylation (Bhutani et al. 2011). Experimental studies conducted in Chinese hamster cell lines (G12) treated with nickel chloride for up to 3 weeks have shown that nickel chloride treatment leads to both promoter hypermethylation and elevated total genomic DNA methylation (Lee YW et al. 1995, 1998). How nickel induces DNA methylation is not yet understood, but it has been proposed that nickel first induces chromatin condensation followed by de novo methylat ion of heterochromatic DNA (Lee YW et al. 1995).

The three SNPs with the strongest associations with LINE-1 methylation were all in TCN2, including two exonic SNPs that result in missense substitutions (rs9606756 and rs9621049) and one intronic SNP (rs4820887). SNPs rs9621049 and rs4820887 have a linkage disequilibrium [r.sup.2] value of 0.8, implying that the observed effect in LINE-1 methylation may eventually be attributed to either of them. TCN2 encodes for transcobalamin II which binds and transports vitamin B12 into the cell (Regec et al. 1995), suggesting that variations in TCN2 could potentially impair the one-carbon metabolism pathway by altering the cytoplasmic concentration of vitamin B12. TCN2-rs9606756 leads to an I23V substitution located at a NAGNAG tandem acceptor site that is a target of alternative splicing (Hiller et al. 2006). TCN2-rs9621049 leads to a S348F and may also play a role in the availability of vitamin B12 in the cell thereby affecting LINE-1 methylation levels.

Four SNPs in other genes (DNMT3A rs7581217, AS3MT-rs7085104, MTHFS-rs1380642, SLC19A1-rs914238) involved in the one-carbon metabolism were also associated with global DNA methylation in our study population. DNMT3A, a de novo DNA methyltransferase, establishes the patterns of methylation in early embryonic development, along with DNMT3B, and cooperates with DNMT1 to maintain the methylation of repetitive sequences such as LINE-1 and Alu elements (Jones and Liang 2009). Recurrent mutations in DNMT3A have been associated with adult hematologic malignancies (Ley et al. 2010; Yan et al. 2011), and mice lacking Dnmt3a die within the first weeks of postnatal life (Robertson 2005). The product of AS3MT catalyzes the conversion of trivalent arsenic by adding a methyl group to mono-methylarsonic acid and dimethylarsonic acid (Ren et al. 2011); monomethylarsonic acid being the most toxic metabolite (Engstrom et al. 2011). rs7085104, located in the promoter region of AS3MT, has been associated with arsenic metabolism, as evidenced by differences in urinary concentration of arsenic metabolites (Engstrom et al. 2011; Valenzuela et al. 2009). We also observed a significant interaction of this SNP with levels of arsenic on LINE-1 methylation levels. Whereas subjects with at least one copy of the major allele had a 4.1% decrease in methylation level per 1-[micro]g/g increase in arsenic, which is comparable to the overall population (-3.6%), those homozygous for the variant allele had a 10% increase in LINE-1 methylation. These findings support the putative functionality of the association. The product of MTHFS catalyzes the conversion of 5-formyltetrahydrofolate to 5,10-methenyltetrahydrofolate and a genome-wide association study reported an association between a variant in this gene and chronic kidney disease (Kottgen et al. 2008). SLC19A1 is a ubiquitously expressed major transporter of folate and antifolates and regulator of the intracellular concentrations of folate (Matherly et al. 2007). Common variants in this gene have been associated with plasma folate levels, various types of cancer (esophageal, gastric, and acute lymphoblastic leukemia), and altered methotrexate transport and adverse effects of methotrexate (Matherly et al. 2007).

Among the limitations of the study is that the majority of paraticipants were of advanced age (mean = 64 years of age, SD = 10 years) and male. This may explain the lack of association between DNA methylation and age in our study population, in contrast with other studies that included subjects with a broader age range (Fraga et al. 2005). Thus, our findings refer to an adult population of mostly men. Results were consistent with estimates for the population as a whole when stratified by sex, with the exception of nickel, tobacco type, and MTHFS-rs1380642, which become nonsignificant when the point estimates were in the same direction (data not shown). These differences may reflect reduced power to estimate associations among women as a result of the small sample size. The presence of missing data for some of the variables might have resulted in decreased power, but even with the available sample size we were able to reproduce previous results and identify novel predictors of LINE-1 methylation. Furthermore, although the study subjects were recruited from hospitals, none of reasons for hospitalization were significantly associated with LINE-1 methylation.

Strengths of the study include its size and the availability and quality of individual data on demographics, lifestyle, environmental exposures, and genetics. Additionally, we assessed LINE-1 methylation levels, which are considered a good marker of global DNA methylation (Yang et al. 2004), using pyro-sequencing, which gives accurate and reproducible measurements (Estecio et al. 2007; Laird 2010; Tost and Gut 2007). Furthermore, this assessment was made using DNA from granulocytes, thereby avoiding a possible effect of cell blood count in our study.

To the best of our knowledge, this is the first study to identify seven SNPs in association with changes in LINE-1 methylation and to integrate different types of information to assess the determinants of global methylation in blood DNA. Integration of both internal and external exposure data in this study is a step forward in understanding how the exposome modulates DNA methylation patterns.

Conclusions

The present study provides further evidence that DNA methylation levels are influenced by variants in genes involved in the one-carbon metabolism pathway, and exposure to trace elements and tobacco smoke. Given the fact that smoking and some of the genetic variants and trace elements associated with LINE-1 methylation in the present study have also been associated with adverse health outcomes including cancer, our results provide additional insight into the potential mechanism through which these agents participate in the development of those diseases. Furthermore, these factors should be considered as potential confounders in etiologic and interventional studies analyzing the role of DNA methylation in disease. Nevertheless, future studies are required to replicate and extend our findings in different populations.

Appendix 1

Spanish Bladder Cancer/EPICURO Study investigators

Institut Municipal d'Investigacio Medica, Universitat Pompeu Fabra, Barcelona-Coordinating Center (M. Kogevinas, N. Malats, F.X. Real, M. Sala, G. Castano, M. Tora, D. Puente, C. Villanueva, C. Murta-Nascimento, J. Fortuny, E. Lopez, S. Hernandez, R. Jaramillo, G. Vellalta, L. Palencia, F. Fermandez, A. Amoros, A. Alfaro, G. Carretero); Hospital del Mar, Universitat Autonoma de Barcelona, Barcelona (J. Lloreta, S. Serrano, L. Ferrer, A. Gelabert, J. Carles, O. Bielsa, K. Villadiego); Hospital Germans Trias i Pujol, Badalona, Barcelona (L. Cecchini, J.M. Saladie, L. Ibarz); Hospital de Sant Boi, Sant Boi de Llobregat, Barcelona (M. Cespedes); Consorci Hospitalari Parc Taulf, Sabadell (C. Serra, D. Garcfa, J. Pujadas, R. Hernando, A. Cabezuelo, C. Abad, A. Prera, J. Prat); Centre Hospitalari i Cardiologic, Manresa, Barcelona (M. Domenech, J. Badal, J. Malet); Hospital Universitario de Canarias, La Laguna, Tenerife (R. Garcfa-Closas, J. Rodrfguez de Vera, A.I. Martfn); Hospital Universitario Nuestra Senora de la Candelaria, Tenerife (J. Tano, F. Caceres); Hospital General Universitario de Elche, Universidad Miguel Hernandez, Elche, Alicante (A. Carrato, F. Garcfa-Lopez, M. Ull, A. Teruel, E. Andrada, A. Bustos, A. Castillejo, J.L. Soto); Universidad de Oviedo, Oviedo, Asturias (A. Tardon); Hospital San Agustfn, Aviles, Asturias (J.L. Guate, J.M. Lanzas, J. Velasco); Hospital Central Covadonga, Oviedo, Asturias (J.M. Fernandez, J.J. Rodrfguez, A. Herrero); Hospital Central General, Oviedo, Asturias (R. Abascal, C. Manzano, T. Miralles); Hospital de Cabuenes, Gijon, Asturias (M. Rivas, M. Arguelles); Hospital de Jove, Gijon, Asturias (M. Dfaz, J. Sanchez, O. Gonzalez); Hospital de Cruz Roja, Gijon, Asturias (A. Mateos, V. Frade); Hospital Alvarez-Buylla, Mieres, Asturias (P. Muntanola, C. Pravia); Hospital Jarrio, Coana, Asturias (A.M. Huescar, F. Huergo); Hospital Carmen y Severo Ochoa, Cangas, Asturias (J. Mosquera).

REFERENCES

Amaral AF, Porta M, Silverman DT, Milne RL, Kogevinas M, Rothman N, et al. 2012. Pancreatic cancer risk and levels of trace elements. Gut 61:1583-1588.

Bailey JA, Carrel L, Chakravarti A, Eichler EE. 2000. Molecular evidence for a relationship between LINE-1 elements and X chromosome inactivation: the Lyon repeat hypothesis. Proc Natl Acad Sci USA 97:6634-6639.

Bhutani N, Burns DM, Blau HM. 2011. DNA demethylation dynamics. Cell 146:866-872.

Breitling LP, Yang R, Korn B, Burwinkel B, Brenner H. 2011. Tobacco-smoking-related differential DNA methylation: 27K discovery and replication. Am J Hum Genet 88:450-457.

Cash HL, Tao L, Yuan JM, Marsit CJ, Houseman EA, Xiang YB, et al. 2012. LINE-1 hypomethylation is associated with bladder cancer risk among nonsmoking chinese. Int J Cancer 130:1151-1159.

Castro R, Rivera I, Ravasco P, Camilo ME, Jakobs C, Blom HJ, et al. 2004. 5,10-methylenetetrahydrofolate reductase (MTHFR) 677CR T and 1298AR C mutations are associated with DNA hypomethylation. J Med Genet 41:454-458.

El-Maarri O, Becker T, Junen J, Manzoor SS, Diaz-Lacava A, Schwaab R, et al. 2007. Gender specific differences in levels of DNA methylation at selected loci from human total blood: a tendency toward higher methylation levels in males. Hum Genet 122:505-514.

El-Maarri O, Walier M, Behne F, van Glim J, Singer H, Diaz-Lacava A, et al. 2011. Methylation at global LINE-1 repeats in human blood are affected by gender but not by age or natural hormone cycles. PLoS One 6:e16252; doi:10.1371/journal.pone.0016252 [Online 19 January 2011].

Engstrom K, Vahter M, Mlakar SJ, Concha G, Nermell B, Raqib R, et al. 2011. Polymorphisms in arsenic(+III oxidation state) methyltransferase [AS3MT] predict gene expression of AS3MT as well as arsenic metabolism. Environ Health Perspect 119:182-188.

Estecio MRH, Gharibyan V, Shen L, Ibrahim AEK, Doshi K, He R, et al. 2007. LINE-1 hypomethylation in cancer is highly variable and inversely correlated with microsatellite instability. PLoS ONE 2:e399; doi:10.1371/journal.pone.0000399 [Online 2 May 2007].

Fraga MF, Ballestar E, Paz MF, Ropero S, Setien F, Ballestar ML, et al. 2005. Epigenetic differences arise during the lifetime of monozygotic twins. Proc Natl Acad Sci USA 102:10604-10609.

Friso S, Choi SW, Girelli D, Mason JB, Dolnikowski GG, Bagley PJ, et al. 2002. A common mutation in the 5,10-methylenetetrahydrofolate reductase gene affects genomic DNA methylation through an interaction with folate status. Proc Natl Acad Sci USA 99:5606-5611.

Fuke C, Shimabukuro M, Petronis A, Sugimoto J, Oda T, Miura K, et al. 2004. Age related changes in 5-methylcytosine content in human peripheral leukocytes and placentas: an HPLC-based study. Ann Hum Genet 68:196-204.

Gabriel HE, Crott JW, Ghandour H, Dallal GE, Choi SW, Keyes MK, et al. 2006. Chronic cigarette smoking is associated with diminished folate status, altered folate form distribution, and increased genetic damage in the buccal mucosa of healthy adults. Am J Clin Nutr 83:835-841.

Garcia-Closas M, Malats N, Real FX, Yeager M, Welch R, Silverman D, et al. 2007. Large-scale evaluation of candidate genes identifies associations between VEGF polymorphisms and bladder cancer risk. PLoS Genet 3:e29; doi:10.1371/journal.pgen.0030029 [Online 23 February 2007].

Garcia-Closas M, Malats N, Silverman D, Dosemeci M, Kogevinas M, Hein DW, et al. 2005. NAT2 slow acetylation, GSTM1 null genotype, and risk of bladder cancer: results from the Spanish Bladder Cancer study and meta-analyses. Lancet 366:649-659.

Garcia-Closas R, Garcia-Closas M, Kogevinas M, Malats N, Silverman D, Serra C, et al. 2007. Food, nutrient and heterocyclic amine intake and the risk of bladder cancer. Eur J Cancer 43:1731-1740.

Hiller M, Huse K, Szafranski K, Jahn N, Hampe J, Schreiber S, et al. 2006. Single-nucleotide polymorphisms in NAGNAG acceptors are highly predictive for variations of alternative splicing. Am J Hum Genet 78:291-302.

Hopkins WA, Staub BP, Baionno JA, Jackson BP, Roe JH, Ford NB. 2004. Trophic and maternal transfer of selenium in brown house snakes (Lamprophis fuliginosus). Ecotoxicol Environ Saf 58:285-293.

Jones PA, Liang G. 2009. Rethinking how DNA methylation patterns are maintained. Nat Rev Genet 10:805-811.

Kottgen A, Kao WH, Hwang SJ, Boerwinkle E, Yang Q, Levy D, et al. 2008. Genome-wide association study for renal traits in the Framingham Heart and Atherosclerosis Risk in Communities Studies. BMC Med Genet 9:49; doi:10.1186/14712350-9-49 [Online 3 June 2008].

Laird PW. 2010. Principles and challenges of genome-wide DNA methylation analysis. Nat Rev Genet 11:191-203.

Lee YL, Xu X, Wallenstein S, Chen J. 2009. Gene expression profiles of the one-carbon metabolism pathway. J Genet Genomics 36:277-282.

Lee YW, Broday L, Costa M. 1998. Effects of nickel on DNA methyltransferase activity and genomic DNA methylation levels. Mutat Res 415:213-218.

Lee YW, Klein CB, Kargacin B, Salnikow K, Kitahara J, Dowjat K, et al. 1995. Carcinogenic nickel silences gene expression by chromatin condensation and DNA methylation: a new model for epigenetic carcinogens. Mol Cell Biol 15:2547-2557.

Lehmann U, Wingen LU, Brakensiek K, Wedemeyer H, Becker T, Heim A, et al. 2007. Epigenetic defects of hepatocellular carcinoma are already found in non-neoplastic liver cells from patients with hereditary haemochromatosis. Hum Mol Genet 16:1335-1342.

Ley TJ, Ding L, Walter MJ, McLellan MD, Lamprecht T, Larson DE, et al. 2010. DNMT3A mutations in acute myeloid leukemia. N Engl J Med 363:2424-2433.

Liu F, Killian JK, Yang M, Walker RL, Hong JA, Zhang M, et al. 2010. Epigenomic alterations and gene expression profiles in respiratory epithelia exposed to cigarette smoke condensate. Oncogene 29:3650-3664.

Malaveille C, Vineis P, Esteve J, Ohshima H, Brun G, Hautefeuille A, et al. 1989. Levels of mutagens in the urine of smokers of black and blond tobacco correlate with their risk of bladder cancer. Carcinogenesis 10:577-586.

Martin-Trujillo A, Iglesias-Platas I, Coto E, Corral-Juan M, San Nicolas H, Corral J, et al. 2011. Genotype of an individual single nucleotide polymorphism regulates DNA methylation at the TRPC3 alternative promoter. Epigenetics 6:1236-1241.

Matherly LH, Hou Z, Deng Y. 2007. Human reduced folate carrier: translation of basic biology to cancer etiology and therapy. Cancer Metastasis Rev 26:111-128.

Moerman JW, Potts GE. 2011. Analysis of metals leached from smoked cigarette litter. Tob Control 20(suppl 1):i30-i35.

Poirier LA, Wise CK, Delongchamp RR, Sinha R. 2001. Blood determinations of S-adenosylmethionine, S-adenosylhomocysteine, and homocysteine: correlations with diet. Cancer Epidemiol Biomarkers Prev 10:649-655.

Regec A, Quadros EV, Platica O, Rothenberg SP. 1995. The cloning and characterization of the human transcobalamin II gene. Blood 85:2711-2719.

Reichard JF, Puga A. 2010. Effects of arsenic exposure on DNA methylation and epigenetic gene regulation. Epigenomics 2:87-104.

Ren X, McHale CM, Skibola CF, Smith AH, Smith MT, Zhang L. 2011. An emerging role for epigenetic dysregulation in arsenic toxicity and carcinogenesis. Environ Health Perspect 119:11-19. Robertson KD. 2005. DNA methylation and human disease. Nat Rev Genet 6:597-610.

Rothman N, Garcia-Closas M, Chatterjee N, Malats N, Wu X, Figueroa JD, et al. 2010. A multi-stage genome-wide association study of bladder cancer identifies multiple susceptibility loci. Nat Genet 42:978-984.

Singer H, Walier M, Nisgen N, Meesters C, Schreiner F, Woelfle J, et al. 2012. Methylation of L1Hs promoters is lower on the inactive X, has a tendency of being higher on autosomes in smaller genomes and shows inter-individual variability at some loci. Hum Mol Genet 21:219-235.

Stover PJ. 2009. One-carbon metabolism-genome interactions in folate-associated pathologies. J Nutr 139:2402-2405.

Terry MB, Delgado-Cruzata L, Vin-Raviv N, Wu HC, Santella RM. 2011. DNA methylation in white blood cells: association with risk factors in epidemiologic studies. Epigenetics 6:828-837.

Tost J, Gut IG. 2007. DNA methylation analysis by pyrosequencing. Nat Protoc 2:2265-2275.

Ulrey CL, Liu L, Andrews LG, Tollefsbol TO. 2005. The impact of metabolism on DNA methylation. Hum Mol Genet 14(Spec 1):R139-R147.

Valenzuela OL, Drobna Z, Hernandez-Castellanos E, Sanchez-Pena LC, Garcia-Vargas GG, Borja-Aburto VH, et al. 2009. Association of AS3MT polymorphisms and the risk of pre-malignant arsenic skin lesions. Toxicol Appl Pharmacol 239:200-207.

Wilhelm CS, Kelsey KT, Butler R, Plaza S, Gagne L, Zens MS, et al. 2010. Implications of LINE1 methylation for bladder cancer risk in women. Clin Cancer Res 16:1682-1689.

Yan XJ, Xu J, Gu ZH, Pan CM, Lu G, Shen Y, et al. 2011. Exome sequencing identifies somatic mutations of DNA methyl-transferase gene DNMT3A in acute monocytic leukemia. Nat Genet 43:309-315.

Yang AS, Estecio MR, Doshi K, Kondo Y, Tajara EH, Issa JP. 2004. A simple method for estimating global DNA methylation using bisulfite PCR of repetitive DNA elements. Nucleic Acids Res 32:e38; doi:10.1093/nar/gnh032 [Online 18 February 2004].

Zhu ZZ, Hou L, Bollati V, Tarantini L, Marinelli B, Cantone L, et al. 2012. Predictors of global methylation levels in blood DNA of healthy subjects: a combined analysis. Int J Epidemiol 41:126-139.

Salman M. Tajuddin, (1) Andre F. S. Amaral, (1) Agustin F. Fernandez, (2) Sandra Rodriguez-Rodero, (2) Ramon Maria Rodriguez, (2) Lee E. Moore, (3) Adonina Tardon, (4) Alfredo Carrato, (5,6) Montserrat Garcia-Closas, (7) Debra T. Silverman, (3) Brian P. Jackson, (8) Reina Garcia-Closas, (9) Ashley L. Cook, (1) Kenneth P. Cantor, (10) Stephen Chanock, (3) Manolis Kogevinas, (11) Nathaniel Rothman, (3) Francisco X. Real, (12,13) Mario F. Fraga, (2,14) and Nuria Malats (1) for the Spanish Bladder Cancer/EPICURO Study Investigators

(1) Genetic and Molecular Epidemiology Group, Spanish National Cancer Research Centre (CNIO), Madrid, Spain; (2) Cancer Epigenetics Laboratory, Instituto Universitario de Oncologfa del Principado de Asturias (IUOPA), Oviedo, Spain; (3) National Cancer Institute, National Institutes of Health, Department of Health and Human Services, Bethesda, Maryland, USA; (4) Universidad de Oviedo, Oviedo, Spain; (5) Hospital General Universitario de Elche, Elche, Spain; (6) Hospital Ramon y Cajal, Madrid, Spain; (7) Division of Genetics and Epidemiology, Institute of Cancer Research, Sutton, UK; (8) Trace Element Analysis Core, Department of Earth Sciences, Dartmouth College, Hanover, New Hampshire, USA; (9) Unidad de Investigacion, Hospital Universitario de Canarias, La Laguna, Spain; (10) KP Cantor Environmental LLC, Silver Spring, Maryland, USA; (11) Centre for Research in Environmental Epidemiology (CREAL), Barcelona, Spain; (12) Epithelial Carcinogenesis Group, Spanish National Cancer Research Centre (CNIO), Madrid, Spain; (13) Departament de Ciencies Experimentals i de la Salut, Universitat Pompeu Fabra, Barcelona, Spain; (14) Department of Immunology and Oncology, Centro Nacional de Biotecnologfa/CNB-CSIC, Cantoblanco, Madrid, Spain

Address correspondence to N. Malats, Genetic and Molecular Epidemiology Group, Spanish National Cancer Research Centre (CNIO), C/Melchor Fernandez Almagro, 3, 28029 Madrid, Spain. Telephone: 34-912-246-900 (ext. 3330). E-mail: nmalats@cnio.es

* These authors contributed equally to this work.

** Spanish Bladder Cancer/EPICURO study investigators are listed in Appendix 1.

Supplemental Material is available online (http:// dx.doi.org/10.1289/ehp.1206068).

We thank the members of S.M.T.'s thesis committee for comments on the manuscript. We also thank the Spanish Bladder Cancer/EPICURO study participants, technicians, field workers, and coordinators.

This work was partially supported by the Association for International Cancer Research (AICR; grant 09-0780, and a doctoral scholarship awarded to S.M.T.); Fondo de Investigaciones Sanitarias, Instituto de Salud Carlos III, MINECO, Spain (grants 00/0745, PI051436, PI061614, PI09-02102, and G03/174);

Red Tematica de Investigacion Cooperativa en Cancer (grant RD06/0020-RTICC); the U.S. National Institutes of Health (grant RO1-CA089715); a postdoctoral fellowship awarded to A.F.S.A. from the Fundacion Cientifica de la AECC; Fundacio Marato TV3; and The Johns Hopkins University Vredenburg Scholarship awarded to A.L.C.

K.P.C. owns KP Cantor Environmental LLC, Silver Spring, MD, USA. The authors declare they have no actual or potential competing financial interests.

Received 26 September 2012; accepted 1 April 2013.
Table 1. Association between LINE/1 methylation level and individual
characteristics of study subjects in the SBC/EPICURO study.

                                               LINE-1
Variables                            n       methylation
                                            mean (95% CI)

Age (years)                         892          --

Sex
  Men                               792   59.0 (58.6, 59.4)
  Women                             100   58.0 (57.1, 58.8)

Region
  Barcelona                         168   59.3 (58.5, 60.2)
  Valles                            135   58.3 (57.5, 59.0)
  Elche                             73    58.2 (57.1, 59.3)
  Tenerife                          145   58.6 (57.8, 59.4)
  Asturias                          371   59.2 (58.6, 59.8)

BMI (kg/ [m.sup.2])
  < 25.0                            372   58.8 (58.3, 59.3)
  25.0-26.99                        148   58.9 (57.9, 59.9)
  27.0-29.99                        120   59.0 (58.1, 59.8)
  [greater than or                  57    58.8 (57.5, 60.2)
  equal to] 30.0
  Missing data                      195

Smoking status
  Nonsmoker                         255   58.2 (57.7, 58.7)
  Occasional smoker                 66    60.0 (58.4, 61.5)
  Former smoker                     329   58.9 (58.3, 59.5)
  Current smoker                    239   59.4 (58.6, 60.2)
  Missing data                       3

Tobacco type
  Nonsmoker                         255   58.2 (57.7, 58.7)
  Blond only                        99    58.5 (57.6, 59.5)
  Black only                        219   59.6 (58.8, 60.4)
  Blond and black                   154   58.7 (57.8, 59.6)
  Unknown                           97    59.3 (58.2, 60.5)
  Missing data                      68

Controls' diagnosis
  Hernia                            332   58.8 (58.2, 59.4)
  Fracture and trauma               263   59.3 (58.6, 60.0)
  Hydrocele                         122   58.7 (57.8, 59.6)
  Other abdominal surgery           99    58.3 (57.4, 59.2)
  Other diseases                    76    59.2 (57.9, 60.5)
Dietary intake (c)
  Vitamin B1 ([micro]g/day/kcal)    645          --
  Vitamin B2 ([micro]g/day/kcal)    645          --
  Vitamin B3 ([micro]g/day/kcal)    645          --
  Vitamin B6 ([micro]g/day/kcal)    645          --
  Vitamin B12 ([micro]g/day/kcal)   645          --
  Folate ([micro]g/day/kcal)        645          --
  Protein ([micro]g/day/kcal)       645          --
  Alcohol ([micro]g/day/kcal)       645          --
  Fruit (g/day/kcal)                639          --
  Vegetable (g/day/kcal)            640          --
  Fruit and vegetable               639          --
 (g/day/kcal)
Toenail trace elements
([micro]g/g) (d)

  Aluminum                          658          --
  Arsenic                           659          --
  Cadmium                           659          --
  Chromium                          658          --
  Copper                            659          --
  Iron                              657          --
  Lead                              659          --
  Manganese                         659          --
  Nickel                            659          --
  Selenium                          659          --
  Vanadium                          651          --
  Zinc                              659          --

NAT2 phenotype
  Rapid/intermediate acetylator     389   59.0 (58.4, 59.6)
  Slow acetylator                   498   58.8 (58.4, 59.2)
  Missing data                       5

GSTM1 genotype
  (+/+, +/-)                        421   58.9 (58.3, 59.4)
  (-/-)                             462   59.0 (58.5, 59.4)
  Missing data                       9

GSTT1 genotype
  (+/+, +/-)                        688   59.0 (58.6, 59.4)
  (-/-)                             198   58.5 (57.9, 59.2)
  Missing data                       6

                                          Unadjusted

Variables                               [beta] (95% CI)        p-
                                                              Value

Age (years)                          -0.004 (-0.02, 0.01)      0.6

Sex
  Men                                         Ref
  Women                                -0.4 (-0.9, 0.05)      0.08

Region
  Barcelona                                   Ref
  Valles                               -0.2 (-0.7, 0.4)        0.5
  Elche                                -0.4 (-1.1, 0.2)        0.2
  Tenerife                             -0.1 (-0.7, 0.4)        0.7
  Asturias                              0.2 (-0.3, 0.6)        0.5

BMI (kg/ [m.sup.2])
  < 25.0                                      Ref
  25.0-26.99                           -0.2 (-0.6, 0.3)        0.5
  27.0-29.99                            0.1 (-0.3, 0.6)        0.6
  [greater than or                     -0.08 (-0.7, 0.6)       0.8
  equal to] 30.0
  Missing data

Smoking status
  Nonsmoker                                   Ref
  Occasional smoker                     0.4 (-0.2, 1.1)        0.2
  Former smoker                        -0.2 (-0.5, 0.3)        0.6
  Current smoker                       -0.2 (-0.6, 0.2)        0.4
  Missing data

Tobacco type
  Nonsmoker                                   Ref
  Blond only                           -0.5 (-1.0, 0.09)       0.1
  Black only                           0.06 (-0.4, 0.5)        0.8
  Blond and black                      -0.3 (-0.8, 0.2)        0.2
  Unknown                              -0.01 (-0.6, 0.6)       0.9
  Missing data

Controls' diagnosis
  Hernia                                      Ref
  Fracture and trauma                  -0.2 (-0.5, 0.2)        0.4
  Hydrocele                             0.2 (-0.3, 0.7)        0.4
  Other abdominal surgery              -0.2 (-0.7, 0.3)        0.4
  Other diseases                       0.01 (-0.9, 0.6)        0.9
Dietary intake (c)
  Vitamin B1 ([micro]g/day/kcal)        0.5 (-0.6, 1.6)        0.3
  Vitamin B2 ([micro]g/day/kcal)        0.1 (-0.5, 0.8)        0.7
  Vitamin B3 ([micro]g/day/kcal)      0.01 (-0.05, 0.08)       0.7
  Vitamin B6 ([micro]g/day/kcal)        0.6 (-0.2, 1.4)        0.1
  Vitamin B12 ([micro]g/day/kcal)     -0.04 (-0.09, 0.01)      0.1
  Folate ([micro]g/day/kcal)         0.001 (-0.002, 0.004)     0.4
  Protein ([micro]g/day/kcal)         0.01 (-0.01, 0.03)       0.4
  Alcohol ([micro]g/day/kcal)        -0.002 (-0.02, 0.02)      0.8
  Fruit (g/day/kcal)                0.0001 (-0.001, 0.002)     0.9
  Vegetable (g/day/kcal)             0.001 (-0.001, 0.003)     0.3
  Fruit and vegetable               0.0003 (-0.0008, 0.001)    0.6
 (g/day/kcal)
Toenail trace elements
([micro]g/g) (d)

  Aluminum                          -0.003 (-0.008, 0.002)     0.2
  Arsenic                              -2.9 (-5.2, -0.6)      0.02
  Cadmium                              0.08 (-0.4, 0.5)        0.7
  Chromium                             0.06 (-0.01, 0.1)      0.09
  Copper                             -0.002 (-0.06, 0.05)     0.95
  Iron                              -0.002 (-0.006, 0.002)     0.4
  Lead                                -0.05 (-0.1, 0.03)       0.2
  Manganese                           -0.03 (-0.1, 0.09)       0.7
  Nickel                              0.02 (0.006, 0.03)      0.002
  Selenium                              0.1 (-0.8, 1.0)        0.8
  Vanadium                             -0.7 (-2.7, 1.3)        0.5
  Zinc                              -0.002 (-0.004, 0.001)     0.2

NAT2 phenotype
  Rapid/intermediate acetylator               Ref
  Slow acetylator                       0.2 (-0.1, 0.5)        0.3
  Missing data

GSTM1 genotype
  (+/+, +/-)                                  Ref
  (-/-)                                0.04 (-0.3, 0.4)        0.8
  Missing data

GSTT1 genotype
  (+/+, +/-)                                  Ref
  (-/-)                                -0.2 (-0.6, 0.2)        0.3
  Missing data

                                    Adjusted

Variables                            n     [beta] (95% CI) (a,b)

Age (years)                         889    -0.006 (-0.02, 0.01)

Sex
  Men                               789             Ref
  Women                             100      -0.7 (-1.2, -0.1)

Region
  Barcelona                         168             Ref
  Valles                            135      -0.2 (-0.7, 0.3)
  Elche                             73       -0.5 (-1.1, 0.2)
  Tenerife                          144      -0.1 (-0.6, 0.4)
  Asturias                          369       0.1 (-0.3, 0.6)

BMI (kg/ [m.sup.2])
  < 25.0                            370             Ref
  25.0-26.99                        148      -0.2 (-0.6, 0.3)
  27.0-29.99                        120       0.2 (-0.3, 0.7)
  [greater than or                  57       -0.03 (-0.7, 0.6)
  equal to] 30.0
  Missing data

Smoking status
  Nonsmoker                         255             Ref
  Occasional smoker                 66        0.3 (-0.3, 1.0)
  Former smoker                     329      -0.3 (-0.7, 0.1)
  Current smoker                    239      -0.4 (-0.9, 0.07)
  Missing data

Tobacco type
  Nonsmoker                         255             Ref
  Blond only                        99      -0.7 (-1.3, -0.08)
  Black only                        218      -0.2 (-0.6, 0.3)
  Blond and black                   154     -0.6 (-1.1, -0.07)
  Unknown                           97       -0.12 (-0.7, 0.5)
  Missing data

Controls' diagnosis
  Hernia                            330             Ref
  Fracture and trauma               262      -0.02 (-0.4, 0.4)
  Hydrocele                         122       0.2 (-0.3, 0.7)
  Other abdominal surgery           99       -0.09 (-0.6, 0.5)
  Other diseases                    76        0.2 (-0.4, 0.8)
Dietary intake (c)
  Vitamin B1 ([micro]g/day/kcal)    644       0.6 (-0.6, 1.7)
  Vitamin B2 ([micro]g/day/kcal)    644       0.2 (-0.5, 0.8)
  Vitamin B3 ([micro]g/day/kcal)    644     0.02 (-0.05, 0.09)
  Vitamin B6 ([micro]g/day/kcal)    644      0.8 (-0.05, 1.6)
  Vitamin B12 ([micro]g/day/kcal)   644     -0.03 (-0.08, 0.02)
  Folate ([micro]g/day/kcal)        644    0.003 (-0.001, 0.01)
  Protein ([micro]g/day/kcal)       644     0.01 (-0.01, 0.03)
  Alcohol ([micro]g/day/kcal)       644     -0.01 (-0.03, 0.02)
  Fruit (g/day/kcal)                638   0.0001 (-0.001, 0.002)
  Vegetable (g/day/kcal)            639    0.002 (-0.001, 0.004)
  Fruit and vegetable               638   0.0005 (-0.0007, 0.002)
 (g/day/kcal)
Toenail trace elements
([micro]g/g) (d)

  Aluminum                          658   -0.003 (-0.008, 0.002)
  Arsenic                           659      -3.6 (-5.9, -1.2)
  Cadmium                           659       0.1 (-0.3, 0.6)
  Chromium                          659     -0.01 (-0.05, 0.03)
  Copper                            659     -0.01 (-0.07, 0.04)
  Iron                              658    0.002 (0.001, 0.004)
  Lead                              659     -0.06 (-0.1, 0.02)
  Manganese                         659     -0.05 (-0.2, 0.06)
  Nickel                            659     0.02 (0.005, 0.03)
  Selenium                          659       0.2 (-0.7 1.2)
  Vanadium                          651      -0.9 (-2.8, 1.2)
  Zinc                              659   -0.001 (-0.004, 0.002)

NAT2 phenotype
  Rapid/intermediate acetylator     388             Ref
  Slow acetylator                   496       0.2 (-0.1, 0.5)
  Missing data

GSTM1 genotype
  (+/+, +/-)                        419             Ref
  (-/-)                             461      0.01 (-0.3, 0.3)
  Missing data

GSTT1 genotype
  (+/+, +/-)                        685             Ref
  (-/-)                             198      -0.2 (-0.6, 0.2)
  Missing data

Variables                            p-
                                    Value

Age (years)                          0.5

Sex
  Men
  Women                             0.02

Region
  Barcelona
  Valles                             0.5
  Elche                              0.1
  Tenerife                           0.7
  Asturias                           0.6

BMI (kg/ [m.sup.2])
  < 25.0
  25.0-26.99                         0.4
  27.0-29.99                         0.4
  [greater than or                   0.9
  equal to] 30.0
  Missing data

Smoking status
  Nonsmoker
  Occasional smoker                  0.3
  Former smoker                      0.2
  Current smoker                     0.1
  Missing data

Tobacco type
  Nonsmoker
  Blond only                        0.03
  Black only                         0.5
  Blond and black                   0.03
  Unknown                            0.7
  Missing data

Controls' diagnosis
  Hernia
  Fracture and trauma                0.9
  Hydrocele                          0.5
  Other abdominal surgery            0.7
  Other diseases                     0.5
Dietary intake (c)
  Vitamin B1 ([micro]g/day/kcal)     0.3
  Vitamin B2 ([micro]g/day/kcal)     0.6
  Vitamin B3 ([micro]g/day/kcal)     0.5
  Vitamin B6 ([micro]g/day/kcal)    0.07
  Vitamin B12 ([micro]g/day/kcal)    0.3
  Folate ([micro]g/day/kcal)         0.1
  Protein ([micro]g/day/kcal)        0.4
  Alcohol ([micro]g/day/kcal)        0.5
  Fruit (g/day/kcal)                 0.9
  Vegetable (g/day/kcal)             0.1
  Fruit and vegetable                0.4
 (g/day/kcal)
Toenail trace elements
([micro]g/g) (d)

  Aluminum                           0.2
  Arsenic                           0.003
  Cadmium                            0.6
  Chromium                           0.6
  Copper                             0.6
  Iron                              0.009
  Lead                               0.2
  Manganese                          0.4
  Nickel                            0.004
  Selenium                           0.6
  Vanadium                           0.4
  Zinc                               0.4

NAT2 phenotype
  Rapid/intermediate acetylator
  Slow acetylator                    0.2
  Missing data

GSTM1 genotype
  (+/+, +/-)
  (-/-)                              0.9
  Missing data

GSTT1 genotype
  (+/+, +/-)
  (-/-)                              0.4
  Missing data

The exposure contrast for trace elements is 1 pg/g and for dietary
variables is 1 pg/day/kcal.

(a) Adjusted for age, sex, region, and smoking status (non-,
occasional, former, current smoker). Tobacco type's [beta] is not
adjusted for smoking status. (b) The number of observations are
reduced by three because of missing data on smoking status. (c) Data
available for those who completed food frequency questionnaire. (d)
Data available for those who provided toenail for trace element
assessment.

Table 2. Association between LINE-1 methylation levels and SNPs in
genes involved in the one-carbon metabolism pathway.

                       dbSNP [chromosome,
Gene           position in the gene, location (a)]      MAF     n

DNMT3A           rs7581217 [2, intron, 25378448]        0.39   875
AS3MT       rs7085104 [10, flanking 5'UTR, 104618863]   0.38   751
                                                               124
MTHFS (c)   rs1380642 [15, flanking 3'UTR, 77883926]    0.18   585
                                                               258
                                                               32
SLC19A1      rs914238 [21, flanking 5'UTR, 45840089]    0.49   231
  (c)                                                          435
                                                               209
TCN2 (d)         rs9621049 [22, exon, 29343419]         0.11   864
                                                               11
                 rs9606756 [22, exon, 29336860]         0.12   864
                                                               11
                rs4820887 [22, intron, 29346914]        0.10   866
                                                                9

Gene           MOI         Genotype

DNMT3A       Additive    per allele T
AS3MT       Recessive        AA/AG
                              GG
MTHFS (c)   Codominant        CC
                              CT
                              TT
SLC19A1     Codominant        TT
  (c)                         TC
                              CC
TCN2 (d)    Recessive        CC/CT
                              TT
            Recessive        AA/AG
                              GG
            Recessive        GG/GA
                              AA

                         Unadjusted

Gene         [beta] (95% CI)      p-Value

DNMT3A       0.3 (0.1, 0.6)        0.003
AS3MT              Ref
             0.8 (0.3, 1.2)       0.0008
MTHFS (c)          Ref             0.03
            0.3 (-0.05, 0.6)
            -0.8 (-1.6, 0.07)
SLC19A1            Ref            0.0008
  (c)        0.5(0.09, 0.8)
            -0.2 (-0.7, 0.2)
TCN2 (d)           Ref
              4.5(3.1, 5.9)     4.3 x 10-10
                   Ref
             2.2 (0.8, 3.6)        0.003
                   Ref
             4.6 (3.0, 6.2)     9.3 x 10-9

                          Adjusted

Gene        [beta] (b) (95% CI)        p-Value

DNMT3A        0.3 (0.1, 0.6)            0.002
AS3MT               Ref
              0.7 (0.3, 1.2)            0.001
MTHFS (c)           Ref                 0.05
             0.3 (-0.08, 0.6)
             -0.8 (-1.6, 0.09)
SLC19A1             Ref                0.0007
  (c)         0.5 (0.08, 0.8)
             -0.3 (-0.7, 0.2)
TCN2 (d)            Ref
              4.2 (2.8, 5.7)      4.7 x [10.sup.-9]
                    Ref
              1.9 (0.5, 3.3)            0.008
                    Ref
              4.0 (2.5, 5.6)      4.8 x [10.sup.-7]

Abbreviations: MAF, minor allele frequency; MOI, mode of inheritance.

(a) Human Genome Build 36.3 location. (b) Adjusted for age, sex,
region, and smoking status. (c) Global p-value for rs1380642 and
rs914238 was estimated by using a two-degree-of-freedom
likelihood-ratio test. (d) Linkage disequilibrium (rs4820887 vs.
rs9621049) [r.sup.2] = 0.8.
COPYRIGHT 2013 National Institute of Environmental Health Sciences
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2013 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research
Author:Tajuddin, Salman M.; Amaral, Andre F.S.; Fernandez, Agustin F.; Rodriguez-Rodero, Sandra; Rodriguez,
Publication:Environmental Health Perspectives
Article Type:Report
Geographic Code:4EUSP
Date:Jun 1, 2013
Words:8783
Previous Article:Linking geological and health sciences to assess childhood lead poisoning from artisanal gold mining in Nigeria.
Next Article:Sulfated metabolites of polychlorinated biphenyls are high-affinity ligands for the thyroid hormone transport protein transthyretin.
Topics:

Terms of use | Privacy policy | Copyright © 2020 Farlex, Inc. | Feedback | For webmasters