Lin Tong

Genome-Wide Association Study Identifies Chromosome 10q24.32 Variants Associated with Arsenic Metabolism and Toxicity Phenotypes in Bangladesh

Arsenic contamination of drinking water is a major public health issue in many countries, increasing risk for a wide array of diseases, including cancer. There is inter-individual variation in arsenic metabolism efficiency and susceptibility to arsenic toxicity; however, the basis of this variation is not well understood. Here, we have performed the first genome-wide association study (GWAS) of arsenic-related metabolism and toxicity phenotypes to improve our understanding of the mechanisms by which arsenic affects health. Using data on urinary arsenic metabolite concentrations and approximately 300,000 genome-wide single nucleotide polymorphisms (SNPs) for 1,313 arsenic-exposed Bangladeshi individuals, we identified genome-wide significant association signals (P<5×10−8) for percentages of both monomethylarsonic acid (MMA) and dimethylarsinic acid (DMA) near the AS3MT gene (arsenite methyltransferase; 10q24.32), with five genetic variants showing independent associations. In a follow-up analysis of 1,085 individuals with arsenic-induced premalignant skin lesions (the classical sign of arsenic toxicity) and 1,794 controls, we show that one of these five variants (rs9527) is also associated with skin lesion risk (P = 0.0005). Using a subset of individuals with prospectively measured arsenic (n = 769), we show that rs9527 interacts with arsenic to influence incident skin lesion risk (P = 0.01). Expression quantitative trait locus (eQTL) analyses of genome-wide expression data from 950 individuals lymphocyte RNA suggest that several of our lead SNPs represent cis-eQTLs for AS3MT (P = 10−12) and neighboring gene C10orf32 (P = 10−44), which are involved in C10orf32-AS3MT read-through transcription. This is the largest and most comprehensive genomic investigation of arsenic metabolism and toxicity to date, the only GWAS of any arsenic-related trait, and the first study to implicate 10q24.32 variants in both arsenic metabolism and arsenical skin lesion risk. The observed patterns of associations suggest that MMA% and DMA% have distinct genetic determinants and support the hypothesis that DMA is the less toxic of these two methylated arsenic species. These results have potential translational implications for the prevention and treatment of arsenic-associated toxicities worldwide.

Arsenic metabolism efficiency has a causal role in arsenic toxicity: Mendelian randomization and gene-environment interaction

BACKGROUND Arsenic exposure through drinking water is a serious global health issue. Observational studies suggest that individuals who metabolize arsenic efficiently are at lower risk for toxicities such as arsenical skin lesions. Using two single nucleotide polymorphisms(SNPs) in the 10q24.32 region (near AS3MT) that show independent associations with metabolism efficiency, Mendelian randomization can be used to assess whether the association between metabolism efficiency and skin lesions is likely to be causal. METHODS Using data on 2060 arsenic-exposed Bangladeshi individuals, we estimated associations for two 10q24.32 SNPs with relative concentrations of three urinary arsenic species (representing metabolism efficiency): inorganic arsenic (iAs), monomethylarsonic acid(MMA) and dimethylarsinic acid (DMA). SNP-based predictions of iAs%, MMA% and DMA% were tested for association with skin lesion status among 2483 cases and 2857 controls. RESULTS Causal odds ratios for skin lesions were 0.90 (95% confidence interval[CI]: 0.87, 0.95), 1.19 (CI: 1.10, 1.28) and 1.23 (CI: 1.12, 1.36)for a one standard deviation increase in DMA%, MMA% and iAs%,respectively. We demonstrated genotype-arsenic interaction, with metabolism-related variants showing stronger associations with skin lesion risk among individuals with high arsenic exposure (synergy index: 1.37; CI: 1.11, 1.62). CONCLUSIONS We provide strong evidence for a causal relationship between arsenic metabolism efficiency and skin lesion risk. Mendelian randomization can be used to assess the causal role of arsenic exposure and metabolism in a wide array of health conditions.exposure and metabolism in a wide array of health conditions.Developing interventions that increase arsenic metabolism efficiency are likely to reduce the impact of arsenic exposure on health.

A Genome-wide Association Study of Early-Onset Breast Cancer Identifies PFKM as a Novel Breast Cancer Gene and Supports a Common Genetic Spectrum for Breast Cancer at Any Age

Early-onset breast cancer (EOBC) causes substantial loss of life and productivity, creating a major burden among women worldwide. We analyzed 1,265,548 Hapmap3 single-nucleotide polymorphisms (SNP) among a discovery set of 3,523 EOBC incident cases and 2,702 population control women ages ≤ 51 years. The SNPs with smallest P values were examined in a replication set of 3,470 EOBC cases and 5,475 control women. We also tested EOBC association with 19,684 genes by annotating each gene with putative functional SNPs, and then combining their P values to obtain a gene-based P value. We examined the gene with smallest P value for replication in 1,145 breast cancer cases and 1,142 control women. The combined discovery and replication sets identified 72 new SNPs associated with EOBC (P < 4 × 10−8) located in six genomic regions previously reported to contain SNPs associated largely with later-onset breast cancer (LOBC). SNP rs2229882 and 10 other SNPs on chromosome 5q11.2 remained associated (P < 6 × 10−4) after adjustment for the strongest published SNPs in the region. Thirty-two of the 82 currently known LOBC SNPs were associated with EOBC (P < 0.05). Low power is likely responsible for the remaining 50 unassociated known LOBC SNPs. The gene-based analysis identified an association between breast cancer and the phosphofructokinase-muscle (PFKM) gene on chromosome 12q13.11 that met the genome-wide gene-based threshold of 2.5 × 10−6. In conclusion, EOBC and LOBC seem to have similar genetic etiologies; the 5q11.2 region may contain multiple distinct breast cancer loci; and the PFKM gene region is worthy of further investigation. These findings should enhance our understanding of the etiology of breast cancer. Cancer Epidemiol Biomarkers Prev; 23(4); 658–69. ©2014 AACR.

Gene-specific differential DNA methylation and chronic arsenic exposure in an epigenome-wide association study of adults in Bangladesh.

Background: Inorganic arsenic is one of the most common naturally occurring contaminants found in the environment. Arsenic is associated with a number of health outcomes, with epigenetic modification suggested as a potential mechanism of toxicity. Objective: Among a sample of 400 adult participants, we evaluated the association between arsenic exposure, as measured by blood and urinary total arsenic concentrations, and epigenome-wide white blood cell DNA methylation. Methods: We used linear regression models to examine the associations between arsenic exposure and methylation at each CpG site, adjusted for sex, age, and batch. Differentially methylated loci were subsequently examined in relation to corresponding gene expression for functional evidence of gene regulation. Results: In adjusted analyses, we observed four differentially methylated CpG sites with urinary total arsenic concentration and three differentially methylated CpG sites with blood arsenic concentration, based on the Bonferroni-corrected significance threshold of p < 1 × 10–7. Methylation of PLA2G2C (probe cg04605617) was the most significantly associated locus in relation to both urinary (p = 3.40 × 10–11) and blood arsenic concentrations (p = 1.48 × 10–11). Three additional novel methylation loci—SQSTM1 (cg01225779), SLC4A4 (cg06121226), and IGH (cg13651690)—were also significantly associated with arsenic exposure. Further, there was evidence of methylation-related gene regulation based on gene expression for a subset of differentially methylated loci. Conclusions: We observed significant associations between arsenic exposure and gene-specific differential white blood cell DNA methylation, suggesting that epigenetic modifications may be an important pathway underlying arsenic toxicity. The specific differentially methylated loci identified may inform potential pathways for future interventions. Citation: Argos M, Chen L, Jasmine F, Tong L, Pierce BL, Roy S, Paul-Brutus R, Gamble MV, Harper KN, Parvez F, Rahman M, Rakibuz-Zaman M, Slavkovich V, Baron JA, Graziano JH, Kibriya MG, Ahsan H. 2015. Gene-specific differential DNA methylation and chronic arsenic exposure in an epigenome-wide association study of adults in Bangladesh. Environ Health Perspect 123:64–71; http://dx.doi.org/10.1289/ehp.1307884

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.

Determinants and Consequences of Arsenic Metabolism Efficiency among 4,794 Individuals: Demographics, Lifestyle, Genetics, and Toxicity

Background: Exposure to inorganic arsenic (iAs), a class I carcinogen, affects several hundred million people worldwide. Once absorbed, iAs is converted to monomethylated (MMA) and then dimethylated forms (DMA), with methylation facilitating urinary excretion. The abundance of each species in urine relative to their sum (iAs%, MMA%, and DMA%) varies across individuals, reflecting differences in arsenic metabolism capacity. Methods: The association of arsenic metabolism phenotypes with participant characteristics and arsenical skin lesions was characterized among 4,794 participants in the Health Effects of Arsenic Longitudinal Study (Araihazar, Bangladesh). Metabolism phenotypes include those obtained from principal component (PC) analysis of arsenic species. Results: Two independent PCs were identified: PC1 appears to represent capacity to produce DMA (second methylation step), and PC2 appears to represent capacity to convert iAs to MMA (first methylation step). PC1 was positively associated (P <0.05) with age, female sex, and BMI, while negatively associated with smoking, arsenic exposure, education, and land ownership. PC2 was positively associated with age and education but negatively associated with female sex and BMI. PC2 was positively associated with skin lesion status, while PC1 was not. 10q24.32/AS3MT region polymorphisms were strongly associated with PC1, but not PC2. Patterns of association for most variables were similar for PC1 and DMA%, and for PC2 and MMA% with the exception of arsenic exposure and SNP associations. Conclusions: Two distinct arsenic metabolism phenotypes show unique associations with age, sex, BMI, 10q24.32 polymorphisms, and skin lesions. Impact: This work enhances our understanding of arsenic metabolism kinetics and toxicity risk profiles. Cancer Epidemiol Biomarkers Prev; 25(2); 381–90. ©2015 AACR.

Genome-wide association study of smoking behaviours among Bangladeshi adults

Background The high prevalence of tobacco use in some developing nations, including Bangladesh, poses several public health challenges for these populations. Smoking behaviour is determined by genetic and environmental factors; however, the genetic determinants of smoking behaviour have not been previously examined in a Bangladeshi or South Asian population. We performed a genome-wide association study (GWAS) of tobacco smoking behaviour among a population-based sample of 5354 (2035 ever smokers and 3319 never smokers) men and women in Bangladesh. Methods Genome-wide association analyses were conducted for smoking initiation (ever vs never smokers), smoking quantity (cigarettes per day), age of smoking initiation, and smoking cessation (former vs current smokers). Sex-stratified associations were performed for smoking initiation. Results We observed associations for smoking initiation in the SLC39A11 region at 17q21.31 (rs2567519, p=1.33×10−7) among men and in the SLCO3A1 region at 15q26 (rs12912184, p=9.32×10−8) among women. Conclusions These findings suggest possible underlying mechanisms related to solute carrier transporter genes, which transport neurotransmitters, nutrients, heavy metals and other substrates into cells, for smoking initiation in a South Asian population in a sex-specific pattern. Genetic markers could have potential translational implications for the prevention or treatment of tobacco use and addiction in South Asian populations and warrant further exploration.

Arsenic exposure, telomere length, and expression of telomere-related genes among Bangladeshi individuals

BACKGROUND Inorganic arsenic is a carcinogen whose mode of action may involve telomere dysfunction. Recent epidemiological studies suggest that chronic arsenic exposure is associated with longer telomeres and altered expression of telomere-related genes in peripheral blood. In this study, we evaluated the association of urinary arsenic concentration with expression of telomere-related genes and telomere length in Bangladeshi individuals with a wide range of arsenic exposure through naturally contaminated drinking water. METHODS We used linear regression models to estimate associations between urinary arsenic and array-based expression measures for 69 telomere related genes using mononuclear cell RNA samples from 1799 individuals. Association between arsenic exposure and a qPCR-based telomere length measure was assessed among 167 individuals. RESULTS Urinary arsenic was positively associated with expression of WRN, and negatively associated with TERF2, DKC1, TERF2IP and OBFC1 (all P<0.00035, Bonferroni-corrected threshold). We detected interaction between urinary arsenic and arsenic metabolism efficiency in relation to expression of WRN (P for interaction =0.00008). In addition, we observed that very high arsenic exposure was associated with longer telomeres compared to very low exposure (P=0.02). DISCUSSION Our findings suggest that arsenics carcinogenic mode of action may involve alteration of telomere maintenance and/or telomere damage. This study extends our knowledge regarding the effect of arsenic on telomere length and expression of telomere-related genes.

Genome-Wide Association Studies and Heritability Estimates of Body Mass Index Related Phenotypes in Bangladeshi Adults

Many health outcomes are influenced by a persons body mass index, as well as by the trajectory of body mass index through a lifetime. Although previous research has established that body mass index related traits are influenced by genetics, the relationship between these traits and genetics has not been well characterized in people of South Asian ancestry. To begin to characterize this relationship, we analyzed the association between common genetic variation and five phenotypes related to body mass index in a population-based sample of 5,354 Bangladeshi adults. We discovered a significant association between SNV rs347313 (intron of NOS1AP) and change in body mass index in women over two years. In a linear mixed-model, the G allele was associated with an increase of 0.25 kg/m2 in body mass index over two years (p-value of 2.3·10−8). We also estimated the heritability of these phenotypes from our genotype data. We found significant estimates of heritability for all of the body mass index-related phenotypes. Our study evaluated the genetic determinants of body mass index related phenotypes for the first time in South Asians. The results suggest that these phenotypes are heritable and some of this heritability is driven by variation that differs from those previously reported. We also provide evidence that the genetic etiology of body mass index related traits may differ by ancestry, sex, and environment, and consequently that these factors should be considered when assessing the genetic determinants of the risk of body mass index-related disease.

The Genetic Architecture of Arsenic Metabolism Efficiency:A SNP-Based Heritability Study of Bangladeshi Adults.

Background Consumption of arsenic-contaminated drinking water adversely affects health. There is interindividual variation in arsenic metabolism efficiency, partially due to genetic variation in the arsenic methyltransferase (AS3MT) gene region. Objectives The goal of this study was to assess the overall contribution of genetic factors to variation in arsenic metabolism efficiency, as measured by the relative concentration of dimethylarsinic acid (DMA%) in urine. Methods Using data on genome-wide single nucleotide polymorphisms (SNPs) and urinary DMA% for 2,053 arsenic-exposed Bangladeshi individuals, we employed various SNP-based approaches for heritability estimation and polygenic modeling. Results Using data on all participants, the percent variance explained (PVE) for DMA% by all measured and imputed SNPs was 16% (p = 0.08), which was reduced to 5% (p = 0.34) after adjusting for AS3MT SNPs. Using information on close relatives only, the PVE was 63% (p = 0.0002), but decreased to 41% (p = 0.01) after adjusting for AS3MT SNPs. Regional heritability analysis confirmed 10q24.32 (AS3MT) as a major arsenic metabolism locus (PVE = 7%, p = 4.4 × 10–10), but revealed no additional regions. We observed a moderate association between a polygenic score reflecting elevated DMA% (composed of thousands of non-AS3MT SNPs) and reduced skin lesion risk in an independent sample (p < 0.05). We observed no associations for SNPs reported in prior candidate gene studies of arsenic metabolism. Conclusions Our results suggest that there are common variants outside of the AS3MT region that influence arsenic metabolism in Bangladeshi individuals, but the effects of these variants are very weak compared with variants near AS3MT. The high heritability estimates observed using family-based heritability approaches suggest substantial effects for rare variants and/or unmeasured environmental factors. Citation Gao J, Tong L, Argos M, Scannell Bryan M, Ahmed A, Rakibuz-Zaman M, Kibriya MG, Jasmine F, Slavkovich V, Graziano JH, Ahsan H, Pierce BL. 2015. The genetic architecture of arsenic metabolism efficiency: a SNP-based heritability study of Bangladeshi adults. Environ Health Perspect 123:985–992; http://dx.doi.org/10.1289/ehp.1408909