Ronald Rahaman

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.

Association Between Arsenic Exposure From Drinking Water and Plasma Levels of Cardiovascular Markers

The authors conducted a cross-sectional study to assess the relation between arsenic exposure from drinking water and plasma levels of markers of systemic inflammation and endothelial dysfunction (matrix metalloproteinase-9, myeloperoxidase, plasminogen activator inhibitor-1, soluble E-selectin, soluble intercellular adhesion molecule-1 (ICAM-1), and soluble vascular adhesion molecule-1 (VCAM-1)) using baseline data from 668 participants (age, >30 years) in the Health Effects of Arsenic Longitudinal Study in Bangladesh (2007-2008). Both well water arsenic and urinary arsenic were positively associated with plasma levels of soluble VCAM-1. For every 1-unit increase in log-transformed well water arsenic (ln μg/L) and urinary arsenic (ln μg/g creatinine), plasma soluble VCAM-1 was 1.02 (95% confidence interval: 1.01, 1.03) and 1.04 (95% confidence interval: 1.01, 1.07) times greater, respectively. There was a significant interaction between arsenic exposure and higher body mass index, such that the increased levels of plasminogen activator inhibitor-1 and soluble VCAM-1 associated with arsenic exposure were stronger among people with higher body mass index. The findings indicate an effect of chronic arsenic exposure from drinking water on vascular inflammation and endothelial dysfunction that could be modified by body mass index and also suggest a potential mechanism underlying the association between arsenic exposure and cardiovascular disease.

A Genome-Wide Study of Cytogenetic Changes in Colorectal Cancer Using SNP Microarrays: Opportunities for Future Personalized Treatment

In colorectal cancer (CRC), chromosomal instability (CIN) is typically studied using comparative-genomic hybridization (CGH) arrays. We studied paired (tumor and surrounding healthy) fresh frozen tissue from 86 CRC patients using Illuminas Infinium-based SNP array. This method allowed us to study CIN in CRC, with simultaneous analysis of copy number (CN) and B-allele frequency (BAF) - a representation of allelic composition. These data helped us to detect mono-allelic and bi-allelic amplifications/deletion, copy neutral loss of heterozygosity, and levels of mosaicism for mixed cell populations, some of which can not be assessed with other methods that do not measure BAF. We identified associations between CN abnormalities and different CRC phenotypes (histological diagnosis, location, tumor grade, stage, MSI and presence of lymph node metastasis). We showed commonalities between regions of CN change observed in CRC and the regions reported in previous studies of other solid cancers (e.g. amplifications of 20q, 13q, 8q, 5p and deletions of 18q, 17p and 8p). From Therapeutic Target Database, we identified relevant drugs, targeted to the genes located in these regions with CN changes, approved or in trials for other cancers and common diseases. These drugs may be considered for future therapeutic trials in CRC, based on personalized cytogenetic diagnosis. We also found many regions, harboring genes, which are not currently targeted by any relevant drugs that may be considered for future drug discovery studies. Our study shows the application of high density SNP arrays for cytogenetic study in CRC and its potential utility for personalized treatment.

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.

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.