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.

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.

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.

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.

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.

A study of telomere length, arsenic exposure, and arsenic toxicity in a Bangladeshi cohort.

BACKGROUND: Chronic arsenic exposure is associated with increased risk for arsenical skin lesions, cancer, and other adverse health outcomes. One potential mechanism of arsenic toxicity is telomere dysfunction. However, prior epidemiological studies of arsenic exposure, telomere length (TL), and skin lesion are small and cross-sectional. We investigated the associations between arsenic exposure and TL and between baseline TL and incident skin lesion risk among individuals participating in the Health Effects of Arsenic Longitudinal Study in Bangladesh (2000-2009). METHODS: Quantitative PCR was used to measure the average TL of peripheral blood DNA collected at baseline. The association between baseline arsenic exposure (well water and urine) and TL was estimated in a randomly-selected subcohort (n = 1469). A nested case-control study (466 cases and 464 age- and sex-matched controls) was used to estimate the association between baseline TL and incident skin lesion risk (diagnosed < 8 years after baseline). RESULTS: No association was observed between arsenic exposure (water or urine) and TL. Among incident skin lesion cases and matched controls, we observed higher skin lesion risk among individuals with shorter TL (Ptrend = 1.5 × 10-5) with odds ratios of 2.60, 1.59, and 1.10 for the first (shortest), second, and third TL quartiles compared to the fourth (longest). CONCLUSIONS: Arsenic exposure was not associated with TL among Bangladeshi adults, suggesting that leukocyte TL may not reflect a primary mode of action for arsenic's toxicity. However, short TL was associated with increased skin lesion risk, and may be a biomarker of arsenic susceptibility modifying arsenic's effect on skin lesion risk.

A novel pooled-sample multiplex luminex assay for high-throughput measurement of relative telomere length.

OBJECTIVES: Relative telomere length (RTL) is a potential biomarker of aging and risk for chronic disease. Previously, we developed a probe-based RTL assay on Luminex platform, where probes for Telomere (T) and reference gene (R) for a given DNA sample were tested in a single well. Here, we describe a method of pooling multiple samples in one well to increase the throughput and cost-effectiveness. METHODS: We used four different microbeads for the same T-probe and four different microbeads for the same R-probe. Each pair of probe sets were hybridized to DNA in separate plates and then pooled in a single plate for all the subsequent steps. We used DNA samples from 60 independent individuals and repeated in multiple batches to test the precision. RESULTS: The precision was good to excellent with Intraclass correlation coefficient (ICC) of 0.908 (95% CI 0.856-0.942). More than 67% of the variation in the RTL could be explained by sample-to-sample variation; less than 0.1% variation was due to batch-to-batch variation and 0.3% variation was explained by bead-to-bead variation. We increased the throughput of RTL Luminex assay from 60 to 240 samples per run. The new assay was validated against the original Luminex assay without pooling (r = 0.79, P = 1.44 × 10-15 ). In an independent set of samples (n = 550), the new assay showed a negative correlation of RTL with age (r = -0.41), a result providing external validation for the method. CONCLUSION: We describe a novel high throughput pooled-sample multiplex Luminex assay for RTL with good to excellent precision suitable for large-scale studies.

Limited influence of germline genetic variation on all-cause mortality in women with early onset breast cancer: evidence from gene-based tests, single-marker regression, and whole-genome prediction.

PURPOSE: Women diagnosed with breast cancer have heterogeneous survival outcomes that cannot be fully explained by known prognostic factors, and germline variation is a plausible but unconfirmed risk factor. METHODS: We used three approaches to test the hypothesis that germline variation drives some differences in survival: mortality loci identification, tumor aggressiveness loci identification, and whole-genome prediction. The 2954 study participants were women diagnosed with breast cancer before age 50, with a median follow-up of 15 years who were genotyped on an exome array. We first searched for loci in gene regions that were associated with all-cause mortality. We next searched for loci in gene regions associated with five histopathological characteristics related to tumor aggressiveness. Last, we also predicted 10-year all-cause mortality on a subset of 1903 participants (3,245,343 variants after imputation) using whole-genome prediction methods. RESULTS: No risk loci for mortality or tumor aggressiveness were identified. This null result persisted when restricting to women with estrogen receptor-positive tumors, when examining suggestive loci in an independent study, and when restricting to previously published risk loci. Additionally, the whole-genome prediction model also found no evidence to support an association. CONCLUSION: Despite multiple complementary approaches, our study found no evidence that mortality in women with early onset breast cancer is influenced by germline variation.

Association between genome-wide copy number variation and arsenic-induced skin lesions: a prospective study.

BACKGROUND: Exposure to arsenic in drinking water is a global health problem and arsenic-induced skin lesions are hallmark of chronic arsenic toxicity. We and others have reported germline genetic variations as risk factors for such skin lesions. The role of copy number variation (CNV) in the germline DNA in this regard is unknown. METHODS: From a large prospectively followed-up cohort, exposed to arsenic, we randomly selected 2171 subjects without arsenic-induced skin lesions at enrollment and genotyped their whole blood DNA samples on Illumina Cyto12v2.1 SNP chips to generate DNA copy number. Participants were followed up every 2 years for a total of 8 years, especially for the development of skin lesions. In Cox regression models, each CNV segment was used as a predictor, accounting for other potential covariates, for incidence of skin lesions. RESULT: The presence of genomic deletion(s) in a number of genes (OR5J2, GOLGA6L7P, APBA2, GALNTL5, VN1R31P, PHKG1P2, SGCZ, ZNF658) and lincRNA genes (RP11-76I14.1, CTC-535 M15.2, RP11-73B2.2) were associated with higher risk [HR between 1.67 (CI 1.3-2.1) and 2.15 (CI 1.5-2.9) for different CNVs] for development of skin lesions independent of gender, age, and arsenic exposure. Some deletions had stronger effect in a specific gender (ZNF658 in males, SGCZ in females) and some had stronger effect in higher arsenic exposure (lincRNA CTD-3179P9.1) suggesting a possible gene-environment interaction. CONCLUSION: This first genome-wide CNV study in a prospectively followed-up large cohort, exposed to arsenic, suggests that DNA deletion in several genes and lincRNA genes may predispose an individual to a higher risk of development of arsenic-induced skin lesions.

Mediation analysis demonstrates that trans-eQTLs are often explained by cis-mediation: a genome-wide analysis among 1,800 South Asians.

A large fraction of human genes are regulated by genetic variation near the transcribed sequence (cis-eQTL, expression quantitative trait locus), and many cis-eQTLs have implications for human disease. Less is known regarding the effects of genetic variation on expression of distant genes (trans-eQTLs) and their biological mechanisms. In this work, we use genome-wide data on SNPs and array-based expression measures from mononuclear cells obtained from a population-based cohort of 1,799 Bangladeshi individuals to characterize cis- and trans-eQTLs and determine if observed trans-eQTL associations are mediated by expression of transcripts in cis with the SNPs showing trans-association, using Sobel tests of mediation. We observed 434 independent trans-eQTL associations at a false-discovery rate of 0.05, and 189 of these trans-eQTLs were also cis-eQTLs (enrichment P<0.0001). Among these 189 trans-eQTL associations, 39 were significantly attenuated after adjusting for a cis-mediator based on Sobel P<10-5. We attempted to replicate 21 of these mediation signals in two European cohorts, and while only 7 trans-eQTL associations were present in one or both cohorts, 6 showed evidence of cis-mediation. Analyses of simulated data show that complete mediation will be observed as partial mediation in the presence of mediator measurement error or imperfect LD between measured and causal variants. Our data demonstrates that trans-associations can become significantly stronger or switch directions after adjusting for a potential mediator. Using simulated data, we demonstrate that this phenomenon is expected in the presence of strong cis-trans confounding and when the measured cis-transcript is correlated with the true (unmeasured) mediator. In conclusion, by applying mediation analysis to eQTL data, we show that a substantial fraction of observed trans-eQTL associations can be explained by cis-mediation. Future studies should focus on understanding the mechanisms underlying widespread cis-mediation and their relevance to disease biology, as well as using mediation analysis to improve eQTL discovery.