The association of cigarette smoking with DNA methylation and gene expression in human tissue samples.

Cigarette smoking adversely affects many aspects of human health, and epigenetic responses to smoking may reflect mechanisms that mediate or defend against these effects. Prior studies of smoking and DNA methylation (DNAm), typically measured in leukocytes, have identified numerous smoking-associated regions (e.g., AHRR). To identify smoking-associated DNAm features in typically inaccessible tissues, we generated array-based DNAm data for 916 tissue samples from the GTEx (Genotype-Tissue Expression) project representing 9 tissue types (lung, colon, ovary, prostate, blood, breast, testis, kidney, and muscle). We identified 6,350 smoking-associated CpGs in lung tissue (n = 212) and 2,735 in colon tissue (n = 210), most not reported previously. For all 7 other tissue types (sample sizes 38-153), no clear associations were observed (false discovery rate 0.05), but some tissues showed enrichment for smoking-associated CpGs reported previously. For 1,646 loci (in lung) and 22 (in colon), smoking was associated with both DNAm and local gene expression. For loci detected in both lung and colon (e.g., AHRR, CYP1B1, CYP1A1), top CpGs often differed between tissues, but similar clusters of hyper- or hypomethylated CpGs were observed, with hypomethylation at regulatory elements corresponding to increased expression. For lung tissue, 17 hallmark gene sets were enriched for smoking-associated CpGs, including xenobiotic- and cancer-related gene sets. At least four smoking-associated regions in lung were impacted by lung methylation quantitative trait loci (QTLs) that co-localize with genome-wide association study (GWAS) signals for lung function (FEV1/FVC), suggesting epigenetic alterations can mediate the effects of smoking on lung health. Our multi-tissue approach has identified smoking-associated regions in disease-relevant tissues, including effects that are shared across tissue types.

Sequencing-based fine-mapping and in silico functional characterization of the 10q24.32 arsenic metabolism efficiency locus across multiple arsenic-exposed populations.

Inorganic arsenic is highly toxic and carcinogenic to humans. Exposed individuals vary in their ability to metabolize arsenic, and variability in arsenic metabolism efficiency (AME) is associated with risks of arsenic-related toxicities. Inherited genetic variation in the 10q24.32 region, near the arsenic methyltransferase (AS3MT) gene, is associated with urine-based measures of AME in multiple arsenic-exposed populations. To identify potential causal variants in this region, we applied fine mapping approaches to targeted sequencing data generated for exposed individuals from Bangladeshi, American Indian, and European American populations (n = 2,357, 557, and 648 respectively). We identified three independent association signals for Bangladeshis, two for American Indians, and one for European Americans. The size of the confidence sets for each signal varied from 4 to 85 variants. There was one signal shared across all three populations, represented by the same SNP in American Indians and European Americans (rs191177668) and in strong linkage disequilibrium (LD) with a lead SNP in Bangladesh (rs145537350). Beyond this shared signal, differences in LD patterns, minor allele frequency (MAF) (e.g., rs12573221 ~13% in Bangladesh ~0.2% among American Indians), and/or heterogeneity in effect sizes across populations likely contributed to the apparent population specificity of the additional identified signals. One of our potential causal variants influences AS3MT expression and nearby DNA methylation in numerous GTEx tissue types (with rs4919690 as a likely causal variant). Several SNPs in our confidence sets overlap transcription factor binding sites and cis-regulatory elements (from ENCODE). Taken together, our analyses reveal multiple potential causal variants in the 10q24.32 region influencing AME, including a variant shared across populations, and elucidate potential biological mechanisms underlying the impact of genetic variation on AME.

Air quality and cancer risk in the All of Us Research Program.

INTRODUCTION: The NIH All of Us Research Program has enrolled over 544,000 participants across the US with unprecedented racial/ethnic diversity, offering opportunities to investigate myriad exposures and diseases. This paper aims to investigate the association between PM2.5 exposure and cancer risks. MATERIALS AND METHODS: This work was performed on data from 409,876 All of Us Research Program participants using the All of Us Researcher Workbench. Cancer case ascertainment was performed using data from electronic health records and the self-reported Personal Medical History questionnaire. PM2.5 exposure was retrieved from NASA's Earth Observing System Data and Information Center and assigned using participants' 3-digit zip code prefixes. Multivariate logistic regression was used to estimate the odds ratio (OR) and 95% confidence interval (CI). Generalized additive models (GAMs) were used to investigate non-linear relationships. RESULTS: A total of 33,387 participants and 46,176 prevalent cancer cases were ascertained from participant EHR data, while 20,297 cases were ascertained from self-reported survey data from 18,133 participants; 9,502 cancer cases were captured in both the EHR and survey data. Average PM2.5 level from 2007 to 2016 was 8.90 µg/m3 (min 2.56, max 15.05). In analysis of cancer cases from EHR, an increased odds for breast cancer (OR 1.17, 95% CI 1.09-1.25), endometrial cancer (OR 1.33, 95% CI 1.09-1.62) and ovarian cancer (OR 1.20, 95% CI 1.01-1.42) in the 4th quartile of exposure compared to the 1st. In GAM, higher PM2.5 concentration was associated with increased odds for blood cancer, bone cancer, brain cancer, breast cancer, colon and rectum cancer, endocrine system cancer, lung cancer, pancreatic cancer, prostate cancer, and thyroid cancer. CONCLUSIONS: We found evidence of an association of PM2.5 with breast, ovarian, and endometrial cancers. There is little to no prior evidence in the literature on the impact of PM2.5 on risk of these cancers, warranting further investigation.

Duration-sensitive association between air pollution exposure and changes in cardiometabolic biomarkers: Evidence from a predominantly African American cohort.

BACKGROUND: Ambient fine particulate matter (PM2.5) exposure has been related to cardiometabolic diseases, but the underlying biological pathways remain unclear at the population level. OBJECTIVE: To investigate the effect of PM2.5 exposure on changes in multiple cardiometabolic biomarkers across different exposure durations. METHOD: Data from a prospective cohort study were analyzed. Ten cardiometabolic biomarkers were measured, including ghrelin, resistin, leptin, C-peptide, creatine kinase myocardial band (CK-MB), monocyte chemoattractant protein-1 (MCP-1), tumor necrosis factor alpha (TNF-alpha), N-terminal pro B-type natriuretic peptide (NT-proBNP), troponin, and interleukin-6 (IL-6). PM2.5 levels across exposure durations from 1 to 36 months were assessed. Mixed effect model was used to estimate changes in biomarker levels against 1 µg/m3 increase in PM2.5 level across different exposure durations. RESULTS: Totally, 641 participants were included. The average PM2.5 exposure level was 9 µg/m3. PM2.5 exposure was inversely associated with ghrelin, and positively associated with all other biomarkers. The magnitudes of these associations were duration-sensitive and exhibited a U-shaped or inverted-U-shaped trend. For example, the association of resistin were ß = 0.05 (95% CI: 0.00, 0.09) for 1-month duration, strengthened to ß = 0.27 (95% CI: 0.14, 0.41) for 13-month duration, and weakened to ß = 0.12 (95% CI: -0.03, 0.26) for 24-month duration. Similar patterns were observed for other biomarkers except for CK-MB, of which the association direction switched from negative to positive as the duration increased. Resistin, leptin, MCP-1, TNF-alpha, and troponin had a sensitive exposure duration of nearly 12 months. Ghrelin and C-peptide were more sensitive to longer-term exposure (>18 months), while NT-proBNP and IL-6 were more sensitive to shorter-term exposure (<6 months). CONCLUSION: PM2.5 exposure was associated with elevated levels in cardiometabolic biomarkers related to insulin resistance, inflammation, and heart injury. The magnitudes of these associations depended on the exposure duration. The most sensitive exposure durations of different biomarkers varied.

Long-term exposure to ambient air pollution and measures of central hemodynamics and arterial stiffness among multiethnic Chicago residents.

OBJECTIVES: To examine whether long-term air pollution exposure is associated with central hemodynamic and brachial artery stiffness parameters. METHODS: We assessed central hemodynamic parameters including central blood pressure, cardiac parameters, systemic vascular compliance and resistance, and brachial artery stiffness measures [including brachial artery distensibility (BAD), compliance (BAC), and resistance (BAR)] using waveform analysis of the arterial pressure signals obtained from a standard cuff sphygmomanometer (DynaPulse2000A, San Diego, CA). The long-term exposures to particles with an aerodynamic diameter < 2.5 µm (PM2.5) and nitrogen dioxide (NO2) for the 3-year periods prior to enrollment were estimated at residential addresses using fine-scale intra-urban spatiotemporal models. Linear mixed models adjusted for potential confounders were used to examine associations between air pollution exposures and health outcomes. RESULTS: The cross-sectional study included 2,387 Chicago residents (76% African Americans) enrolled in the ChicagO Multiethnic Prevention And Surveillance Study (COMPASS) during 2013-2018 with validated address information, PM2.5 or NO2, key covariates, and hemodynamics measurements. We observed long-term concentrations of PM2.5 and NO2 to be positively associated with central systolic, pulse pressure and BAR, and negatively associated with BAD, and BAC after adjusting for relevant covariates. A 1-µg/m3 increment in preceding 3-year exposures to PM2.5 was associated with 1.8 mmHg higher central systolic (95% CI: 0.98, 4.16), 1.0 mmHg higher central pulse pressure (95% CI: 0.42, 2.87), a 0.56%mmHg lower BAD (95% CI: -0.81, -0.30), and a 0.009 mL/mmHg lower BAC (95% CI: -0.01, -0.01). CONCLUSION: This population-based study provides evidence that long-term exposures to PM2.5 and NO2 is related to central BP and arterial stiffness parameters, especially among African Americans.

Association of Microsatellite Instability and Gene Expression Profile in Colorectal Carcinoma and Potential Implications for Therapy.

Background and Objective: In sporadic colorectal carcinomas (CRC), microsatellite instability (MSI) pathways play important roles. Previously, we showed differences in DNA methylation patterns in microsatellite stable (MSS) colorectal carcinomas and MSI-CRC. In the current study, we explore the similarities and differences in gene expression profiles in MSS and MSI at the gene level and at the pathway level to better understand CRC pathogenesis and/or the potential for therapeutic opportunities. Material and Methods: Seventy-one CRC patients (MSI = 18, MSS = 53) were studied. Paired tumor and adjacent normal tissues were used for genome-wide gene expression assays. Result: At the gene level, we compared the list of differentially expressed genes (fold change (FC) ≥ 3 and FDR < 0.05) in tumor tissues compared to corresponding normal tissue in CRC patients with MSI tumors (190 genes) and MSS tumors (129 genes). Of these, 107 genes overlapped. The list of genes that were differentially expressed in MSI tumors only showed enrichment predominantly in two broad categories of pathways-(a) Inflammation-related pathways including the interleukin-17 (IL-17) signaling pathway, tumor necrosis factor (TNF) signaling pathway, chemokine signaling, nuclear factor kappa B (NFκB) signaling, and cytokine-cytokine interactions, and (b) metabolism-related pathways, including retinol metabolism, steroid hormone biosynthesis, drug metabolism, pentose and glucoronate interconversions, and ascorbate and aldarate metabolism. The genes in inflammation-related pathways were up-regulated whereas genes in metabolism-related pathways were down-regulated in MSI tumor tissue. Pathway-level analysis also revealed similar results confirming the gene enrichment findings. For example, the 150 genes involved in the IL-17 signaling pathway were on average up-regulated by 1.19 fold (CI 1.16-1.21) in MSI compared to 1.14 fold (CI 1.13-1.16) in MSS patients (interaction p = 0.0009). Conclusions: We document an association between MSI status and differential gene expression that broadens our understanding of CRC pathogenesis. Furthermore, targeting one or more of these dysregulated pathways could provide the basis for improved therapies for MSI and MSS CRC.

The Impact of Neighborhood Disadvantage on Asthma Prevalence in a Predominantly African-American, Chicago-Based Cohort.

This study aimed to investigate the joint effect of neighborhood disadvantages on asthma prevalence and evaluate whether individual-level variables protect residents against neighborhood disadvantages. Data from the Chicago Multiethnic Prevention and Surveillance Study (from 2013-2020) were analyzed. Eight neighborhood characteristics were measured using the Chicago Health Atlas, including neighborhood unsafety, limited access to healthy food, neighborhood alienation, severe rent burden, vacant housing, single-parent household, neighborhood poverty, and unemployment. A structured questionnaire measured asthma diagnosis (childhood or adulthood) and individual-level variables including sex, age, income, education, and race. Weighted quantile sum (WQS) regression was used to evaluate the impact of neighborhood disadvantages. Stratified analysis was performed by income and education. A total of 6,592 participants (mean age = 53.5 (standard deviation, 11.1) years) were included. Most of the study population were non-Hispanic Black (82.5%) and reported an annual household income less than $15,000 (53%). Asthma prevalence was 23.6%. The WQS index, which represents the overall neighborhood disadvantages, was associated with asthma prevalence (odds ratio = 1.14, 95% confidence interval: 1.07, 1.22) when adjusted for individual-level confounders. Neighborhood poverty contributed 40.8% to the overall impact, followed by vacant housing (23.1%) and neighborhood alienation (22.9%). When stratified by individual-level income or education, no difference was observed for the association between WQS index and asthma prevalence.

A metabolome-wide case-control study of african american breast cancer patients.

BACKGROUND: Breast cancer survivors face long-term sequelae compared to the general population, suggesting altered metabolic profiles after breast cancer. We used metabolomics approaches to investigate the metabolic differences between breast cancer patients and women in the general population, aiming to elaborate metabolic changes among breast cancer patients and identify potential targets for clinical interventions to mitigate long-term sequelae. METHODS: Serum samples were retrieved from 125 breast cancer cases recruited from the Chicago Multiethnic Epidemiologic Breast Cancer Cohort (ChiMEC), and 125 healthy controls selected from Chicago Multiethnic Prevention and Surveillance Study (COMPASS). We used liquid chromatography-high resolution mass spectrometry to obtain untargeted metabolic profiles and partial least squares discriminant analysis (PLS-DA) combined with fold change to select metabolic features associated with breast cancer. Pathway analyses were conducted using Mummichog to identify differentially enriched metabolic pathways among cancer patients. As potential confounders we included age, marital status, tobacco smoking, alcohol drinking, type 2 diabetes, and area deprivation index in our model. Random effects of residence for intercept was also included in the model. We further conducted subgroup analysis by treatment timing (chemotherapy/radiotherapy/surgery), lymph node status, and cancer stages. RESULTS: The entire study participants were African American. The average ages were 57.1 for cases and 58.0 for controls. We extracted 15,829 features in total, among which 507 features were eventually selected by our criteria. Pathway enrichment analysis of these 507 features identified three differentially enriched metabolic pathways related to prostaglandin, leukotriene, and glycerophospholipid. The three pathways demonstrated inconsistent patterns. Metabolic features in the prostaglandin and leukotriene pathways exhibited increased abundances among cancer patients. In contrast, metabolic intensity in the glycerolphospholipid pathway was deregulated among cancer patients. Subgroup analysis yielded consistent results. However, changes in these pathways were strengthened when only using cases with positive lymph nodes, and attenuated when only using cases with stage I disease. CONCLUSION: Breast cancer in African American women is associated with increase in serum metabolites involved in prostaglandin and leukotriene pathways, but with decrease in serum metabolites in glycerolphospholipid pathway. Positive lymph nodes and advanced cancer stage may strengthen changes in these pathways.

Arsenic exposure from drinking water and endothelial dysfunction in Bangladeshi adolescents.

Cardiovascular disease (CVD) is the leading cause of morbidity and mortality worldwide, with ∼80% of CVD-related deaths occurring in low- and middle-income countries. Growing evidence suggests that chronic arsenic exposure may contribute to CVD through its effect on endothelial function in adults. However, few studies have examined the influence of arsenic exposure on cardiovascular health in children and adolescents. To examine arsenic's relation to preclinical markers of endothelial dysfunction, we enrolled 200 adolescent children (ages 15-19 years; median 17) of adult participants in the Health Effects of Arsenic Longitudinal Study (HEALS), in Araihazar, Bangladesh. Participants' arsenic exposure was determined by recall of lifetime well usage for drinking water. As part of HEALS, wells were color-coded to indicate arsenic level (<10 µg/L, 10-50 µg/L, >50 µg/L). Endothelial function was measured by recording fingertip arterial pulsatile volume change and reactive hyperemia index (RHI) score, an independent CVD risk factor, was calculated from these measurements. In linear regression models adjusted for participant's sex, age, education, maternal education, land ownership and body weight, individuals who reported always drinking water from wells with >50 µg/L arsenic had a 11.75% lower level of RHI (95% CI: -21.26, -1.09, p = 0.03), as compared to participants who drank exclusively from wells with ≤50 µg/L arsenic. Sex-stratified analyses suggest that these associations were stronger in female participants. As compared to individuals who drank exclusively from wells with ≤50 µg/L arsenic, the use of wells with >50 µg/L arsenic was associated with 14.36% lower RHI (95% CI: -25.69, -1.29, p = 0.03) in females, as compared to 5.35% lower RHI (95% CI: -22.28, 15.37, p = 0.58) in males for the same comparison. Our results suggest that chronic arsenic exposure may be related to endothelial dysfunction in adolescents, especially among females. Further work is needed to confirm these findings and examine whether these changes may increase risk of later adverse cardiovascular health events.

A missense variant in FTCD is associated with arsenic metabolism and toxicity phenotypes in Bangladesh.

Inorganic arsenic (iAs) is a carcinogen, and exposure to iAs via food and water is a global public health problem. iAs-contaminated drinking water alone affects >100 million people worldwide, including ~50 million in Bangladesh. Once absorbed into the blood stream, most iAs is converted to mono-methylated (MMA) and then di-methylated (DMA) forms, facilitating excretion in urine. Arsenic metabolism efficiency varies among individuals, in part due to genetic variation near AS3MT (arsenite methyltransferase; 10q24.32). To identify additional arsenic metabolism loci, we measured protein-coding variants across the human exome for 1,660 Bangladeshi individuals participating in the Health Effects of Arsenic Longitudinal Study (HEALS). Among the 19,992 coding variants analyzed exome-wide, the minor allele (A) of rs61735836 (p.Val101Met) in exon 3 of FTCD (formiminotransferase cyclodeaminase) was associated with increased urinary iAs% (P = 8x10-13), increased MMA% (P = 2x10-16) and decreased DMA% (P = 6x10-23). Among 2,401 individuals with arsenic-induced skin lesions (an indicator of arsenic toxicity and cancer risk) and 2,472 controls, carrying the low-efficiency A allele (frequency = 7%) was associated with increased skin lesion risk (odds ratio = 1.35; P = 1x10-5). rs61735836 is in weak linkage disequilibrium with all nearby variants. The high-efficiency/major allele (G/Valine) is human-specific and eliminates a start codon at the first 5´-proximal Kozak sequence in FTCD, suggesting selection against an alternative translation start site. FTCD is critical for catabolism of histidine, a process that generates one-carbon units that can enter the one-carbon/folate cycle, which provides methyl groups for arsenic metabolism. In our study population, FTCD and AS3MT SNPs together explain ~10% of the variation in DMA% and support a causal effect of arsenic metabolism efficiency on arsenic toxicity (i.e., skin lesions). In summary, this work identifies a coding variant in FTCD associated with arsenic metabolism efficiency, providing new evidence supporting the established link between one-carbon/folate metabolism and arsenic toxicity.

Correction: A missense variant in FTCD is associated with arsenic metabolism and toxicity phenotypes in Bangladesh.

[This corrects the article DOI: 10.1371/journal.pgen.1007984.].

A Transcriptome and Methylome Study Comparing Tissues of Early and Late Onset Colorectal Carcinoma.

There is an increase in the incidence of early onset colorectal carcinoma (EOCRC). To better understand if there is any difference in molecular pathogenesis of EOCRC and late onset colorectal carcinoma (LOCRC), we compared the clinical, histological, transcriptome, and methylome profile of paired CRC and healthy colonic tissue from 67 EOCRC and 98 LOCRC patients. The frequency of stage 3 CRC, lymph node involvement, lymphovascular invasion, and perineural invasion was higher in the EOCRC group. Many of the cancer related pathways were differentially expressed in CRC tissue in both EOCRC and LOCRC patients. However, the magnitude of differential expression for some groups of genes, such as DNA damage repair genes and replication stress genes, were significantly less pronounced in the EOCRC group, suggesting less efficient DNA damage repair to be associated with EOCRC. A more marked methylation of "growth factor receptor" genes in LOCRC correlated with a more pronounced down-regulation of those genes in that group. From a therapeutic point of view, more over-expression of fatty acid synthase (FASN) among the LOCRC patients may suggest a better response of FASN targeted therapy in that group. The age of onset of CRC did not appear to modify the response of cis-platin or certain immune checkpoint inhibitors. We found some differences in the molecular pathogenesis in EOCRC and LOCRC that may have some biological and therapeutic significance.

Clinical Factors Affecting Bevacizumab Efficacy With and Without Conventional Chemotherapy in Metastatic Colon Cancer.

PURPOSE: Bevacizumab (BZ) combined with first line chemotherapy (CC) has shown good clinical outcomes in metastatic colorectal cancer (mCRC). Overall survival (OS) and/or progression free survival in mCRC patients receiving BZ with or without 5FU-based CC is thought to be affected by clinical and morphological factor(s). PATIENTS AND METHODS: We reviewed retrospective medical records of all consecutive mCRC patients treated with BZ with or without CC at tertiary care center between 2003 and 2009 out of which149 patients (m = 77, f = 72) were eligible. RESULTS: Our study population had a mean age at diagnosis of 63.5 years (SD = 11) with median follow-up period of 19.4 months. On initial radiological evaluation following BZ therapy, 56 patients (m = 31, f = 25) had complete or partial response categorized as "early responders." Remaining patients (m = 46, f = 47) who were either stable or showed progressive disease were categorized as "non-responders." Fifty percent among early responders and 60% among non-responders [relative risk (RR) 0.67 (95% confidence interval (CI), 0.43-1.06)] demonstrated disease progression on follow up. There was a slightly better OS among early responders compared to non-responders (median 21.5 months days versus 16.8 months, P = 0.07). Cox regression analysis suggested male sex (RR 0.65, 95% CI, 0.43-0.98), hematochezia (RR 0.63, 95% CI, 0.4-0.98), resectable primary tumor (RR 0.42, 95% CI, 0.24-0.72) and resectable metastatic mass (RR 0.32, 95% CI, 0.14-0.74) were found to be associated with longer OS. Abdominal pain (RR 1.76, 95% CI, 1.1-2.8), accompanying diabetes (RR 1.76, 95% CI, 1.09-2.85), and unexplained weight loss (RR 2.73, 95% CI, 1.73-4.29) were associated with poor OS. CONCLUSIONS: Better OS among mCRC patients with resectable primary and metastatic tumors was seen. This is the first study to demonstrate slightly better outcome in males and negative influence of diabetes on outcome in mCRC treated with BZ.

The Association Between Smoking and Gut Microbiome in Bangladesh.

INTRODUCTION: Epidemiological studies that investigate alterations in the gut microbial composition associated with smoking are lacking. This study examined the composition of the gut microbiome in smokers compared with nonsmokers. AIMS AND METHODS: Stool samples were collected in a cross-sectional study of 249 participants selected from the Health Effects of Arsenic Longitudinal Study in Bangladesh. Microbial DNA was extracted from the fecal samples and sequenced by 16S rRNA gene sequencing. The associations of smoking status and intensity of smoking with the relative abundance or the absence and presence of individual bacterial taxon from phylum to genus levels were examined. RESULTS: The relative abundance of bacterial taxa along the Erysipelotrichi-to-Catenibacterium lineage was significantly higher in current smokers compared to never-smokers. The odds ratio comparing the mean relative abundance in current smokers with that in never-smokers was 1.91 (95% confidence interval = 1.36-2.69) for the genus Catenibacterium and 1.89 (95% confidence interval = 1.39-2.56) for the family Erysipelotrichaceae, the order Erysipelotrichale, and the class Erysipelotrichi (false discovery rate-adjusted p values = .0008-.01). A dose-response association was observed for each of these bacterial taxa. The presence of Alphaproteobacteria was significantly greater comparing current with never-smokers (odds ratio = 4.85, false discovery rate-adjusted p values = .04). CONCLUSIONS: Our data in a Bangladeshi population are consistent with evidence of an association between smoking status and dosage with change in the gut bacterial composition. IMPLICATIONS: This study for the first time examined the relationship between smoking and the gut microbiome composition. The data suggest that smoking status may play an important role in the composition of the gut microbiome, especially among individuals with higher levels of tobacco exposure.

The effect of age on DNA methylation in whole blood among Bangladeshi men and women.

BACKGROUND: It is well-known that methylation changes occur as humans age, however, understanding how age-related changes in DNA methylation vary by sex is lacking. In this study, we characterize the effect of age on DNA methylation in a sex-specific manner and determine if these effects vary by genomic context. We used the Illumina HumanMethylation 450 K array and DNA derived from whole blood for 400 adult participants (189 males and 211 females) from Bangladesh to identify age-associated CpG sites and regions and characterize the location of these age-associated sites with respect to CpG islands (vs. shore, shelf, or open sea) and gene regions (vs. intergenic). We conducted a genome-wide search for age-associated CpG sites (among 423,604 sites) using a reference-free approach to adjust for cell type composition (the R package RefFreeEWAS) and performed an independent replication analysis of age-associated CpGs. RESULTS: The number of age-associated CpGs (p < 5 x 10- 8) were 986 among men and 3479 among women of which 2027(63.8%) and 572 (64.1%) replicated (using Bonferroni adjusted p < 1.2 × 10- 5). For both sexes, age-associated CpG sites were more likely to be hyper-methylated with increasing age (compared to hypo-methylated) and were enriched in CpG islands and promoter regions compared with other locations and all CpGs on the array. Although we observed strong correlation between chronological age and previously-developed epigenetic age models (r ≈ 0.8), among our top (based on lowest p-value) age-associated CpG sites only 12 for males and 44 for females are included in these prediction models, and the median chronological age compared to predicted age was 44 vs. 51.7 in males and 45 vs. 52.1 in females. CONCLUSIONS: Our results describe genome-wide features of age-related changes in DNA methylation. The observed associations between age and methylation were generally consistent for both sexes, although the associations tended to be stronger among women. Our population may have unique age-related methylation changes that are not captured in the established methylation-based age prediction model we used, which was developed to be non-tissue-specific.

The contribution of parent-to-offspring transmission of telomeres to the heritability of telomere length in humans.

Leukocyte telomere length (LTL) is a heritable trait with two potential sources of heritability (h2): inherited variation in non-telomeric regions (e.g., SNPs that influence telomere maintenance) and variability in the lengths of telomeres in gametes that produce offspring zygotes (i.e., "direct" inheritance). Prior studies of LTL h2 have not attempted to disentangle these two sources. Here, we use a novel approach for detecting the direct inheritance of telomeres by studying the association between identity-by-descent (IBD) sharing at chromosome ends and phenotypic similarity in LTL. We measured genome-wide SNPs and LTL for a sample of 5069 Bangladeshi adults with substantial relatedness. For each of the 6318 relative pairs identified, we used SNPs near the telomeres to estimate the number of chromosome ends shared IBD, a proxy for the number of telomeres shared IBD (Tshared). We then estimated the association between Tshared and the squared pairwise difference in LTL ((ΔLTL)2) within various classes of relatives (siblings, avuncular, cousins, and distant), adjusting for overall genetic relatedness (ϕ). The association between Tshared and (ΔLTL)2 was inverse among all relative pair types. In a meta-analysis including all relative pairs (ϕ > 0.05), the association between Tshared and (ΔLTL)2 (P = 0.01) was stronger than the association between ϕ and (ΔLTL)2 (P = 0.43). Our results provide strong evidence that telomere length (TL) in parental germ cells impacts TL in offspring cells and contributes to LTL h2 despite telomere "reprogramming" during embryonic development. Applying our method to larger studies will enable robust estimation of LTL h2 attributable to direct transmission of telomeres.

Lack of Association between Opioid-Receptor Genotypes and Smoking Cessation Outcomes in a Randomized, Controlled Naltrexone Trial.

AIMS: The present study examined how variation in mu- (OPRM1), kappa- (OPRK), and delta- (OPRD) opioid receptor genes may influence the efficacy of naltrexone in the context of a smoking cessation trial. METHODS: The study's primary objective was to examine the association of the Asn40Asp OPRM1 single nucleotide polymorphism (SNP) with naltrexone's effects on smoking quit rate, weight gain, and heavy drinking behavior during a double-blind, randomized clinical trial in 280 adult DSM-IV nicotine-dependent participants. The secondary goal of the study was to examine the relationship of 20 additional SNPs of OPRM1, OPRK, and OPRD with the aforementioned outcomes. RESULTS: Results indicated a null association between any opioid-receptor gene SNP and naltrexone's effects on smoking quit rate, weight gain, and heavy drinking behavior in this sample of nicotine dependent participants. CONCLUSIONS: In sum, these results do not suggest that genetic variation in opioid-receptors is related to treatment responses to naltrexone in a smoking cessation trial.

Genome-wide association study of telomere length among South Asians identifies a second RTEL1 association signal.

BACKGROUND: Leucocyte telomere length (TL) is a potential biomarker of ageing and risk for age-related disease. Leucocyte TL is heritable and shows substantial differences by race/ethnicity. Recent genome-wide association studies (GWAS) report ~10 loci harbouring SNPs associated with leucocyte TL, but these studies focus primarily on populations of European ancestry. OBJECTIVE: This study aims to enhance our understanding of genetic determinants of TL across populations. METHODS: We performed a GWAS of TL using data on 5075 Bangladeshi adults. We measured TL using one of two technologies (qPCR or a Luminex-based method) and used standardised variables as TL phenotypes. RESULTS: Our results replicate previously reported associations in the TERC and TERT regions (P=2.2×10-8 and P=6.4×10-6, respectively). We observed a novel association signal in the RTEL1 gene (intronic SNP rs2297439; P=2.82×10-7) that is independent of previously reported TL-associated SNPs in this region. The minor allele for rs2297439 is common in South Asian populations (≥0.25) but at lower frequencies in other populations (eg, 0.07 in Northern Europeans). Among the eight other previously reported association signals, all were directionally consistent with our study, but only rs8105767 (ZNF208) was nominally significant (P=0.003). SNP-based heritability estimates were as high as 44% when analysing close relatives but much lower when analysing distant relatives only. CONCLUSIONS: In this first GWAS of TL in a South Asian population, we replicate some, but not all, of the loci reported in prior GWAS of individuals of European ancestry, and we identify a novel second association signal at the RTEL1 locus.

Screening for gene-environment (G×E) interaction using omics data from exposed individuals: an application to gene-arsenic interaction.

Identifying gene-environment interactions is a central challenge in the quest to understand susceptibility to complex, multi-factorial diseases. Developing an understanding of how inter-individual variability in inherited genetic variation alters the effects of environmental exposures will enhance our knowledge of disease mechanisms and improve our ability to predict disease and target interventions to high-risk sub-populations. Limited progress has been made identifying gene-environment interactions in the epidemiological setting using existing statistical approaches for genome-wide searches for interaction. In this paper, we describe a novel two-step approach using omics data to conduct genome-wide searches for gene-environment interactions. Using existing genome-wide SNP data from a large Bangladeshi cohort study specifically designed to assess the effect of arsenic exposure on health, we evaluated gene-arsenic interactions by first conducting genome-wide searches for SNPs that modify the effect of arsenic on molecular phenotypes (gene expression and DNA methylation features). Using this set of SNPs showing evidence of interaction with arsenic in relation to molecular phenotypes, we then tested SNP-arsenic interactions in relation to skin lesions, a hallmark characteristic of arsenic toxicity. With the emergence of additional omics data in the epidemiologic setting, our approach may have the potential to boost power for genome-wide interaction research, enabling the identification of interactions that will enhance our understanding of disease etiology and our ability to develop interventions targeted at susceptible sub-populations.

Occupational pesticide exposure and subclinical hypothyroidism among male pesticide applicators.

OBJECTIVES: Animal studies suggest that exposure to pesticides may alter thyroid function; however, few epidemiologic studies have examined this association. We evaluated the relationship between individual pesticides and thyroid function in 679 men enrolled in a substudy of the Agricultural Health Study, a cohort of licensed pesticide applicators. METHODS: Self-reported lifetime pesticide use was obtained at cohort enrolment (1993-1997). Intensity-weighted lifetime days were computed for 33 pesticides, which adjusts cumulative days of pesticide use for factors that modify exposure (eg, use of personal protective equipment). Thyroid-stimulating hormone (TSH), thyroxine (T4), triiodothyronine (T3) and antithyroid peroxidase (anti-TPO) autoantibodies were measured in serum collected in 2010-2013. We used multivariate logistic regression to estimate ORs and 95% CIs for subclinical hypothyroidism (TSH >4.5 mIU/L) compared with normal TSH (0.4-<4.5 mIU/L) and for anti-TPO positivity. We also examined pesticide associations with TSH, T4 and T3 in multivariate linear regression models. RESULTS: Higher exposure to the insecticide aldrin (third and fourth quartiles of intensity-weighted days vs no exposure) was positively associated with subclinical hypothyroidism (ORQ3=4.15, 95% CI 1.56 to 11.01, ORQ4=4.76, 95% CI 1.53 to 14.82, ptrend <0.01), higher TSH (ptrend=0.01) and lower T4 (ptrend=0.04). Higher exposure to the herbicide pendimethalin was associated with subclinical hypothyroidism (fourth quartile vs no exposure: ORQ4=2.78, 95% CI 1.30 to 5.95, ptrend=0.02), higher TSH (ptrend=0.04) and anti-TPO positivity (ptrend=0.01). The fumigant methyl bromide was inversely associated with TSH (ptrend=0.02) and positively associated with T4 (ptrend=0.01). CONCLUSIONS: Our results suggest that long-term exposure to aldrin, pendimethalin and methyl bromide may alter thyroid function among male pesticide applicators.