Studies of DNA methylation (DNAm) in solid human tissues are relatively scarce; tissue-specific characterization of DNAm is needed to understand its role in gene regulation and its relevance to ...complex traits. We generated array-based DNAm profiles for 987 human samples from the Genotype-Tissue Expression (GTEx) project, representing 9 tissue types and 424 subjects. We characterized methylome and transcriptome correlations (eQTMs), genetic regulation in cis (mQTLs and eQTLs) across tissues and e/mQTLs links to complex traits. We identified mQTLs for 286,152 CpG sites, many of which (>5%) show tissue specificity, and mQTL colocalizations with 2,254 distinct GWAS hits across 83 traits. For 91% of these loci, a candidate gene link was identified by integration of functional maps, including eQTMs, and/or eQTL colocalization, but only 33% of loci involved an eQTL and mQTL present in the same tissue type. With this DNAm-focused integrative analysis, we contribute to the understanding of molecular regulatory mechanisms in human tissues and their impact on complex traits.
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.
The gut microbiome likely plays a role in the etiology of multiple health conditions, especially those affecting the gastrointestinal tract. Little consensus exists as to the best, standard methods ...to collect fecal samples for future microbiome analysis. We evaluated three distinct populations (N = 132 participants) using 16S rRNA gene amplicon sequencing data to investigate the reproducibility, stability, and accuracy of microbial profiles in fecal samples collected and stored via fecal occult blood test (FOBT) or Flinders Technology Associates (FTA) cards, fecal immunochemical tests (FIT) tubes, 70% and 95% ethanol, RNAlater, or with no solution. For each collection method, based on relative abundance of select phyla and genera, two alpha diversity metrics, and four beta diversity metrics, we calculated intraclass correlation coefficients (ICCs) to estimate reproducibility and stability, and Spearman correlation coefficients (SCCs) to estimate accuracy of the fecal microbial profile. Comparing duplicate samples, reproducibility ICCs for all collection methods were excellent (ICCs ≥75%). After 4-7 days at ambient temperature, ICCs for microbial profile stability were excellent (≥75%) for most collection methods, except those collected via no-solution and 70% ethanol. SCCs comparing each collection method to immediately-frozen no-solution samples ranged from fair to excellent for most methods; however, accuracy of genus-level relative abundances differed by collection method. Our findings, taken together with previous studies and feasibility considerations, indicated that FOBT/FTA cards, FIT tubes, 95% ethanol, and RNAlater are excellent choices for fecal sample collection methods in future microbiome studies. Furthermore, establishing standard collection methods across studies is highly desirable.
Kirsten Rat Sarcoma (KRAS) is the most commonly mutated oncogene in colorectal carcinoma (CRC). We have previously reported the interactions between microsatellite instability (MSI), DNA promoter ...methylation, and gene expression. In this study, we looked for associations between KRAS mutation, gene expression, and methylation that may help with precision medicine. Genome-wide gene expression and DNA methylation were done in paired CRC tumor and surrounding healthy tissues. The results suggested that (a) the magnitude of dysregulation of many major gene pathways in CRC was significantly greater in patients with the KRAS mutation, (b) the up- and down-regulation of these dysregulated gene pathways could be correlated with the corresponding hypo- and hyper-methylation, and (c) the up-regulation of CDKN2A was more pronounced in tumors with the KRAS mutation. A recent cell line study showed that there were higher CDKN2A levels in 5-FU-resistant CRC cells and that these could be down-regulated by Villosol. Our findings suggest the possibility of a better response to anti-CDKN2A therapy with Villosol in KRAS-mutant CRC. Also, the more marked up-regulation of genes in the proteasome pathway in CRC tissue, especially with the KRAS mutation and MSI, may suggest a potential role of a proteasome inhibitor (bortezomib, carfilzomib, or ixazomib) in selected CRC patients if necessary.
Inherited genetic variation affects local gene expression and DNA methylation in humans. Most expression quantitative trait loci (cis-eQTLs) occur at the same genomic location as a methylation QTL ...(cis-meQTL), suggesting a common causal variant and shared mechanism. Using DNA and RNA from peripheral blood of Bangladeshi individuals, here we use co-localization methods to identify eQTL-meQTL pairs likely to share a causal variant. We use partial correlation and mediation analyses to identify >400 of these pairs showing evidence of a causal relationship between expression and methylation (i.e., shared mechanism) with many additional pairs we are underpowered to detect. These co-localized pairs are enriched for SNPs showing opposite associations with expression and methylation, although many SNPs affect multiple CpGs in opposite directions. This work demonstrates the pervasiveness of co-regulated expression and methylation in the human genome. Applying this approach to other types of molecular QTLs can enhance our understanding of regulatory mechanisms.
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.
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.
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 (
) among the LOCRC patients may suggest a better response of
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.
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.
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.
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).
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.
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.
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.