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.

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 arsenic's 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 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.

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 individual's 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.

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.

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 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 Illumina's 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.

Interpretation of genome-wide infinium methylation data from ligated DNA in formalin-fixed, paraffin-embedded paired tumor and normal tissue.

BACKGROUND: Formalin-fixed, paraffin-embedded (FFPE) samples are a highly desirable resource for epigenetic studies, but there is no suitable platform to assay genome-wide methylation in these widely available resources. Recently, Thirlwell et al. (2010) have reported a modified ligation-based DNA repair protocol to prepare FFPE DNA for the Infinium methylation assay. In this study, we have tested the accuracy of methylation data obtained with this modification by comparing paired fresh-frozen (FF) and FFPE colon tissue (normal and tumor) from colorectal cancer patients. We report locus-specific correlation and concordance of tumor-specific differentially methylated loci (DML), both of which were not previously assessed. METHODS: We used Illumina's Infinium Methylation 27K chip for 12 pairs of FF and 12 pairs of FFPE tissue from tumor and surrounding healthy tissue from the resected colon of the same individual, after repairing the FFPE DNA using Thirlwell's modified protocol. RESULTS: For both tumor and normal tissue, overall correlation of ß values between all loci in paired FF and FFPE was comparable to previous studies. Tissue storage type (FF or FFPE) was found to be the most significant source of variation rather than tissue type (normal or tumor). We found a large number of DML between FF and FFPE DNA. Using ANOVA, we also identified DML in tumor compared to normal tissue in both FF and FFPE samples, and out of the top 50 loci in both groups only 7 were common, indicating poor concordance. Likewise, while looking at the correlation of individual loci between FFPE and FF across the patients, less than 10% of loci showed strong correlation (r ≥ 0.6). Finally, we checked the effect of the ligation-based modification on the Infinium chemistry for SNP genotyping on an independent set of samples, which also showed poor performance. CONCLUSION: Ligation of FFPE DNA prior to the Infinium genome-wide methylation assay may detect a reasonable number of loci, but the numbers of detected loci are much fewer than in FF samples. More importantly, the concordance of DML detected between FF and FFPE DNA is suboptimal, and DML from FFPE tissues should be interpreted with great caution.

A genome-wide DNA methylation study in colorectal carcinoma.

BACKGROUND: We performed a genome-wide scan of 27,578 CpG loci covering 14,475 genes to identify differentially methylated loci (DML) in colorectal carcinoma (CRC). METHODS: We used Illumina's Infinium methylation assay in paired DNA samples extracted from 24 fresh frozen CRC tissues and their corresponding normal colon tissues from 24 consecutive diagnosed patients at a tertiary medical center. RESULTS: We found a total of 627 DML in CRC covering 513 genes, of which 535 are novel DML covering 465 genes. We also validated the Illumina Infinium methylation data for top-ranking genes by non-bisulfite conversion q-PCR-based methyl profiler assay in a subset of the same samples. We also carried out integration of genome-wide copy number and expression microarray along with methylation profiling to see the functional effect of methylation. Gene Set Enrichment Analysis (GSEA) showed that among the major "gene sets" that are hypermethylated in CRC are the sets: "inhibition of adenylate cyclase activity by G-protein signaling", "Rac guanyl-nucleotide exchange factor activity", "regulation of retinoic acid receptor signaling pathway" and "estrogen receptor activity". Two-level nested cross validation showed that DML-based predictive models may offer reasonable sensitivity (around 89%), specificity (around 95%), positive predictive value (around 95%) and negative predictive value (around 89%), suggesting that these markers may have potential clinical application. CONCLUSION: Our genome-wide methylation study in CRC clearly supports most of the previous findings; additionally we found a large number of novel DML in CRC tissue. If confirmed in future studies, these findings may lead to identification of genomic markers for potential clinical application.