Gary S. Pittman

Discovery and verification of functional single nucleotide polymorphisms in regulatory genomic regions: Current and developing technologies

The most common form of genetic variation, single nucleotide polymorphisms or SNPs, can affect the way an individual responds to the environment and modify disease risk. Although most of the millions of SNPs have little or no effect on gene regulation and protein activity, there are many circumstances where base changes can have deleterious effects. Non-synonymous SNPs that result in amino acid changes in proteins have been studied because of their obvious impact on protein activity. It is well known that SNPs within regulatory regions of the genome can result in disregulation of gene transcription. However, the impact of SNPs located in putative regulatory regions, or rSNPs, is harder to predict for two primary reasons. First, the mechanistic roles of non-coding genomic sequence remain poorly defined. Second, experimental validation of the functional consequences of rSNPs is often slow and laborious. In this review, we summarize traditional and novel methodologies for candidate rSNPs selection, in particular in silico techniques that aid in candidate rSNP selection. Additionally we will discuss molecular biological techniques that assess the impact of rSNPs on binding of regulatory machinery, as well as functional consequences on transcription. Standard techniques such as EMSA and luciferase reporter constructs are still widely used to assess effects of rSNPs on binding and gene transcription; however, these protocols are often bottlenecks in the discovery process. Therefore, we highlight novel and developing high-throughput protocols that promise to aid in shortening the process of rSNP validation. Given the large amount of genomic information generated from a multitude of re-sequencing and genome-wide SNP array efforts, future focus should be to develop validation techniques that will allow greater understanding of the impact these polymorphisms have on human health and disease.

Glutathione S-transferase genotype as a susceptibility factor in smoking-related coronary heart disease.

Cancer studies suggest that the null polymorphisms of glutathione S-transferase M1 or T1 (GSTM1/GSTT1) may affect the ability to detoxify or activate chemicals in cigarette smoke. The potential modification of the association between smoking and coronary heart disease (CHD) by GSTM1 and GSTT1 has not been studied in humans. A case-cohort study was conducted to test the hypotheses that specific genotypes of GSTM1 or GSTT1 affect susceptibility to smoking-related CHD. CHD cases (n=400) accrued during 1987-1993 and a cohort-representative sample (n=924) were selected from a biracial cohort of 15792 middle-aged men and women in four US communities. A significantly higher frequency of GSTM1-0 and a lower frequency of GSTT1-0 were found in whites (GSTM1-0=47.1%, GSTT1-0=16.4%) than in African-Americans (AAs) (GSTM1-0=17.5%, GSTT1-0=25.9%). A smoking-GSTM1-0 interaction for the risk of CHD was statistically significant on an additive scale, with ever-smokers with GSTM1-0 at a approximately 1.5-fold higher risk relative to ever-smokers with GSTM1-1 and a approximately 2-fold higher risk relative to never-smokers with GSTM1-0, after adjustment for other CHD risk factors. The interaction between having smoked >/=20 pack-years and GSTT1-1 was statistically significant on both multiplicative and additive scales. The risk of CHD given both GSTT1-1 and >/=20 pack-years of smoking was approximately three times greater than the risk given exposure to >/=20 pack-years of smoking alone, and approximately four times greater than the risk given exposure to GSTT1-1 alone. The modification of the smoking-CHD association by GSTM1 or GSTT1 suggests that chemicals in cigarette smoke that are substrates for glutathione S-transferases may be involved in the etiology of CHD.

Manganese superoxide dismutase Ala-9Val polymorphism and risk of breast cancer in a population-based case–control study of African Americans and whites

IntroductionA polymorphism in the manganese superoxide dismutase (MnSOD) gene, Ala-9Val, has been examined in association with breast cancer risk in several epidemiologic studies. Results suggest that the Ala allele increases the risk of breast cancer and modifies the effects of environmental exposures that produce oxidative damage to DNA.MethodsWe examined the role of the MnSOD Ala-9Val polymorphism in a population-based case–control study of invasive and in situ breast cancer in North Carolina. Genotypes were evaluated for 2025 cases (760 African Americans and 1265 whites) and for 1812 controls (677 African Americans and 1135 whites).ResultsThe odds ratio for MnSOD Ala/Ala versus any MnSOD Val genotypes was not elevated in African Americans (odds ratio = 0.9, 95% confidence interval = 0.7–1.2) or in whites (odds ratio = 1.0, 95% confidence interval = 0.8–1.2). Greater than additive joint effects were observed for the Ala/Ala genotype and smoking, radiation to the chest, and occupational exposure to ionizing radiation. Antagonism was observed between the Ala/Ala genotype and the use of nonsteroidal anti-inflammatory drugs.ConclusionsThe MnSOD genotype may contribute to an increased risk of breast cancer in the presence of specific environmental exposures. These results provide further evidence for the importance of reactive oxygen species and of oxidative DNA damage in the etiology of breast cancer.

Genetic Variation and Antioxidant Response Gene Expression in the Bronchial Airway Epithelium of Smokers at Risk for Lung Cancer

Prior microarray studies of smokers at high risk for lung cancer have demonstrated that heterogeneity in bronchial airway epithelial cell gene expression response to smoking can serve as an early diagnostic biomarker for lung cancer. As a first step in applying functional genomic analysis to population studies, we have examined the relationship between gene expression variation and genetic variation in a central molecular pathway (NRF2-mediated antioxidant response) associated with smoking exposure and lung cancer. We assessed global gene expression in histologically normal airway epithelial cells obtained at bronchoscopy from smokers who developed lung cancer (SC, n = 20), smokers without lung cancer (SNC, n = 24), and never smokers (NS, n = 8). Functional enrichment analysis showed that the NRF2-mediated, antioxidant response element (ARE)-regulated genes, were significantly lower in SC, when compared with expression levels in SNC. Importantly, we found that the expression of MAFG (a binding partner of NRF2) was correlated with the expression of ARE genes, suggesting MAFG levels may limit target gene induction. Bioinformatically we identified single nucleotide polymorphisms (SNPs) in putative ARE genes and to test the impact of genetic variation, we genotyped these putative regulatory SNPs and other tag SNPs in selected NRF2 pathway genes. Sequencing MAFG locus, we identified 30 novel SNPs and two were associated with either gene expression or lung cancer status among smokers. This work demonstrates an analysis approach that integrates bioinformatics pathway and transcription factor binding site analysis with genotype, gene expression and disease status to identify SNPs that may be associated with individual differences in gene expression and/or cancer status in smokers. These polymorphisms might ultimately contribute to lung cancer risk via their effect on the airway gene expression response to tobacco-smoke exposure.

Probing the functional impact of sequence variation on p53-DNA interactions using a novel microsphere assay for protein-DNA binding with human cell extracts.

The p53 tumor suppressor regulates its target genes through sequence-specific binding to DNA response elements (REs). Although numerous p53 REs are established, the thousands more identified by bioinformatics are not easily subjected to comparative functional evaluation. To examine the relationship between RE sequence variation—including polymorphisms—and p53 binding, we have developed a multiplex format microsphere assay of protein-DNA binding (MAPD) for p53 in nuclear extracts. Using MAPD we measured sequence-specific p53 binding of doxorubicin-activated or transiently expressed p53 to REs from established p53 target genes and p53 consensus REs. To assess the sensitivity and scalability of the assay, we tested 16 variants of the p21 target sequence and a 62-multiplex set of single nucleotide (nt) variants of the p53 consensus sequence and found many changes in p53 binding that are not captured by current computational binding models. A group of eight single nucleotide polymorphisms (SNPs) was examined and binding profiles closely matched transactivation capability tested in luciferase constructs. The in vitro binding characteristics of p53 in nuclear extracts recapitulated the cellular in vivo transactivation capabilities for eight well-established human REs measured by luciferase assay. Using a set of 26 bona fide REs, we observed distinct binding patterns characteristic of transiently expressed wild type and mutant p53s. This microsphere assay system utilizes biologically meaningful cell extracts in a multiplexed, quantitative, in vitro format that provides a powerful experimental tool for elucidating the functional impact of sequence polymorphism and protein variation on protein/DNA binding in transcriptional networks.

DNA Methylation of the Aryl Hydrocarbon Receptor Repressor Associations with Cigarette Smoking and Subclinical Atherosclerosis

Background—Tobacco smoke contains numerous agonists of the aryl hydrocarbon receptor (AhR) pathway, and activation of the AhR pathway was shown to promote atherosclerosis in mice. Intriguingly, cigarette smoking is most strongly and robustly associated with DNA modifications to an AhR pathway gene, the AhR repressor (AHRR). We hypothesized that altered AHRR methylation in monocytes, a cell type sensitive to cigarette smoking and involved in atherogenesis, may be a part of the biological link between cigarette smoking and atherosclerosis. Methods and Results—DNA methylation profiles of AHRR in monocytes (542 CpG sites±150 kb of AHRR, using Illumina 450K array) were integrated with smoking habits and ultrasound-measured carotid plaque scores from 1256 participants of the Multi-Ethnic Study of Atherosclerosis (MESA). Methylation of cg05575921 significantly associated (P=6.1×10−134) with smoking status (current versus never). Novel associations between cg05575921 methylation and carotid plaque scores (P=3.1×10−10) were identified, which remained significant in current and former smokers even after adjusting for self-reported smoking habits, urinary cotinine, and well-known cardiovascular disease risk factors. This association replicated in an independent cohort using hepatic DNA (n=141). Functionally, cg05575921 was located in a predicted gene expression regulatory element (enhancer) and had methylation correlated with AHRR mRNA profiles (P=1.4×10−17) obtained from RNA sequencing conducted on a subset (n=373) of the samples. Conclusions—These findings suggest that AHRR methylation may be functionally related to AHRR expression in monocytes and represents a potential biomarker of subclinical atherosclerosis in smokers.

Variation in genes relevant to aromatic hydrocarbon metabolism and the risk of adult brain tumors

Genes involved in phase I and phase II regulation of aromatic hydrocarbon-induced effects exhibit sequence variability that may mediate the risk of adult brain tumors. We evaluated associations between gene variants in CYP1A1, CYP1B1, GSTM3, EPHX1, and NQO1 and adult brain tumor incidence. Cases were patients with glioma (n = 489), meningioma (n = 197), or acoustic neuroma (n = 96) diagnosed from 1994 to 1998 at three U.S. hospitals. Controls were 799 patients admitted to the same hospitals for nonmalignant conditions. DNA was extracted from blood samples collected from 1277 subjects, and genotyping was conducted for CYP1A1 I462V, CYP1B1 V432L, EPHX1 Y113H, GSTM3 *A/*B (intron 6 deletion), and NQO1 P187S. The CYP1B1 V432L homozygous variant was associated with decreased risk of meningioma (odds ratio [OR] = 0.6; 95% CI, 0.3-1.0) but not the other tumor types. The GSTM3 *B/*B genotype was associated with increased risk of glioma (OR = 2.3; 95% CI, 1.0-5.2) and meningioma (OR = 3.6; 95% CI, 1.3-9.8). Increased risks associated with GSTM3 *B/*B were observed in younger subjects (age < 50) and older subjects (age > or = 50), in men and women, and within each study site. The magnitude of association for GSTM3 with glioma and meningioma was greater among ever-smokers than among those who had never smoked. None of the other genotypes showed consistent associations with any tumor type. The association with the GSTM3 *B allele, while intriguing, requires replication, and additional research is needed to clarify the function of the GSTM3 alleles studied here.

Human single-nucleotide polymorphisms alter p53 sequence-specific binding at gene regulatory elements

p53 coordinates the expression of an intricate network of genes in response to stress signals. Sequence-specific DNA binding is essential for p53-mediated tumor suppression. We evaluated the impact of single-nucleotide polymorphisms (SNPs) in p53 response elements (p53RE) on DNA binding and gene expression in response to DNA damage. Using a bioinformatics approach based on incorporating p53 binding strength into a position weight matrix, we selected 32 SNPs in putative and validated p53REs. The microsphere assay for protein–DNA binding (MAPD) and allele-specific expression analysis was employed to assess the impact of SNPs on p53-DNA binding and gene expression, respectively. Comparing activated p53 binding in nuclear extracts from doxorubicin- or ionizing radiation (IR)-treated human cells, we observed little difference in binding profiles. Significant p53 binding was observed for most polymorphic REs and several displayed binding comparable to the p21 RE. SNP alleles predicted to lower p53 binding indeed reduced binding in 25 of the 32 sequences. Chromatin immunoprecipitation-sequencing in lymphoblastoid cells confirmed p53 binding to seven polymorphic p53 REs in response to doxorubicin. In addition, five polymorphisms were associated with altered gene expression following doxorubicin treatment. Our findings demonstrate an effective strategy to identify and evaluate SNPs that may alter p53-mediated stress responses.

Distinct epigenetic effects of tobacco smoking in whole blood and among leukocyte subtypes

Tobacco smoke exposure dramatically alters DNA methylation in blood cells and may mediate smoking-associated complex diseases through effects on immune cell function. However, knowledge of smoking effects in specific leukocyte subtypes is limited. To better characterize smoking–associated methylation changes in whole blood and leukocyte subtypes, we used Illumina 450K arrays and Reduced Representation Bisulfite Sequencing (RRBS) to assess genome-wide DNA methylation. Differential methylation analysis in whole blood DNA from 172 smokers and 81 nonsmokers revealed 738 CpGs, including 616 previously unreported CpGs, genome-wide significantly associated with current smoking (p <1.2x10-7, Bonferroni correction). Several CpGs (MTSS1, NKX6-2, BTG2) were associated with smoking duration among heavy smokers (>22 cigarettes/day, n = 86) which might relate to long-term heavy-smoking pathology. In purified leukocyte subtypes from an independent group of 20 smokers and 14 nonsmokers we further examined methylation and gene expression for selected genes among CD14+ monocytes, CD15+ granulocytes, CD19+ B cells, and CD2+ T cells. In 10 smokers and 10 nonsmokers we used RRBS to fine map differential methylation in CD4+ T cells, CD8+ T cells, CD14+, CD15+, CD19+, and CD56+ natural killer cells. Distinct cell-type differences in smoking-associated methylation and gene expression were identified. AHRR (cg05575921), ALPPL2 (cg21566642), GFI1 (cg09935388), IER3 (cg06126421) and F2RL3 (cg03636183) showed a distinct pattern of significant smoking-associated methylation differences across cell types: granulocytes> monocytes>> B cells. In contrast GPR15 (cg19859270) was highly significant in T and B cells and ITGAL (cg09099830) significant only in T cells. Numerous other CpGs displayed distinctive cell-type responses to tobacco smoke exposure that were not apparent in whole blood DNA. Assessing the overlap between these CpG sites and differential methylated regions (DMRs) with RRBS in 6 cell types, we confirmed cell-type specificity in the context of DMRs. We identified new CpGs associated with current smoking, pack-years, duration, and revealed unique profiles of smoking-associated DNA methylation and gene expression among immune cell types, providing potential clues to hematopoietic lineage-specific effects in disease etiology.

Abstract 3647: Dose-dependent alteration of CpG methylation in AHRR and GFI1 in mononuclear cell DNA of smokers.

About 45 million (∼19%) of U.S. adults smoke, 46 million are former smokers and >100 million are exposed to environmental tobacco smoke. Tobacco use is associated with many types of cancer, heart disease, stroke, and respiratory disease. While hundreds of tobacco smoke constituents cause DNA damage, many adverse outcomes are not related to DNA damage and an emerging hypothesis is that exposure-induced epigenetic effects may mediate many of these outcomes. Using epigenome-wide association, Joubert et al (Env Health Perspectives, 2012) observed Bonferroni-corrected statistically significant associations (480,000 tests, p Citation Format: Xuting Wang, Gary S. Pittman, Dan Su, Kelly N. Adamski, Michelle R. Campbell, Bonnie R. Joubert, Zhiqing Y. Huang, Cathrine Hoyo, Susan K. Murphy, Stephanie A. London, Douglas A. Bell. Dose-dependent alteration of CpG methylation in AHRR and GFI1 in mononuclear cell DNA of smokers. [abstract]. In: Proceedings of the 104th Annual Meeting of the American Association for Cancer Research; 2013 Apr 6-10; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2013;73(8 Suppl):Abstract nr 3647. doi:10.1158/1538-7445.AM2013-3647