Search

Supplementary Figures S1-2,4, describe cell viability curves of tested cell lines. Figure S3 and S5 show bioinformatic data for resistance upregulated gene sets (S3) and cisplatin-taxane comparison analysis(S5). Figure S6 shows ECAR and OCR experiments. Figure S7 is a heatmap of combination indices of cisplatin with JQ1. Figure S8 shows apoptosis rates of cell lines. S9-12 describe cell viability dynamics of SOX9 KOs.

Chemoresistance is driven by unique regulatory networks in the genome that are distinct from those necessary for cancer development. Here, we investigate the contribution of enhancer elements to cisplatin resistance in ovarian cancers. Epigenome profiling of multiple cellular models of chemoresistance identified unique sets of distal enhancers, super-enhancers (SE), and their gene targets that coordinate and maintain the transcriptional program of the platinum-resistant state in ovarian cancer. Pharmacologic inhibition of distal enhancers through small-molecule epigenetic inhibitors suppressed the expression of their target genes and restored cisplatin sensitivity in vitro and in vivo. In addition to known drivers of chemoresistance, our findings identified SOX9 as a critical SE-regulated transcription factor that plays a critical role in acquiring and maintaining the chemoresistant state in ovarian cancer. The approach and findings presented here suggest that integrative analysis of epigenome and transcriptional programs could identify targetable key drivers of chemoresistance in cancers.Significance:Integrative genome-wide epigenomic and transcriptomic analyses of platinum-sensitive and -resistant ovarian lines identify key distal regulatory regions and associated master regulator transcription factors that can be targeted by small-molecule epigenetic inhibitors.

Immune checkpoint blockade (ICB) therapy has improved long-term survival for patients with advanced melanoma. However, there is critical need to identify potential biomarkers of response and actionable strategies to improve response rates. Through generation and analysis of 148 chromatin modification maps for 36 melanoma samples from patients treated with anti-PD- 1, we identified significant enrichment of active enhancer states in non-responders at baseline. Analysis of an independent cohort of 20 samples identified a set of 437 enhancers that predicted response to anti-PD-1 therapy (Area Under the Curve of 0.8417). The activated non-responder enhancers marked a group of key regulators of several pathways in melanoma cells (including c- MET, TGFβ, EMT and AKT) that are known to mediate resistance to ICB therapy and several checkpoint receptors in T cells. Epigenetic editing experiments implicated involvement of c-MET enhancers in the modulation of immune response. Finally, inhibition of enhancers and repression of these pathways using bromodomain inhibitors along with anti-PD-1 therapy significantly decreased melanoma tumor burden and increased T-cell infiltration. Together, these findings identify a potential enhancer-based biomarker of resistance to anti-PD-1 and suggest enhancer blockade in combination with ICB as a potential strategy to improve responses.

Inguinal white adipose tissue (iWAT) is essential for the beneficial effects of exercise training on metabolic health. Extracellular matrix (ECM) composition, innervation, and vascularization are all important regulators of iWAT function, yet whether exercise training improves these structural components of iWAT is unknown. Using biochemical, imaging, and multi-omics analyses we find that 11-days of wheel running in male mice causes profound iWAT remodeling including decreased ECM deposition and increased vascularization and innervation. We identify adipose stem cells as the main contributors to training-induced ECM remodeling, determine that training causes a shift from hypertrophic to insulin-sensitive adipocyte subpopulations, show that the PRDM16 transcriptional complex is necessary for iWAT remodeling and beiging, and discover neuronal growth regulator 1 (NEGR1) as a link between PRDM16 and neuritogenesis. Exercise training leads to remarkable adaptations to iWAT structure and cell-type composition that can confer beneficial changes in tissue metabolism.Graphical Abstract

Background: While CTLA-4 and PD-1 immune-checkpoint blockade (ICB) have revolutionized melanoma treatment, these treatments have different mechanisms, response rates, and toxicities. Differential predictors of response are poorly characterized. Combination ICB has the highest response rates but also much higher toxicity, and which patients to treat with single agent vs. combination ICB is controversial. In this study, we set out to understand the differential biology of anti-CTLA-4 and anti-PD-1 by integrating transcriptomic and clinical characterization from large cohorts of patients with metastatic melanoma treated with anti-CTLA-4, anti-PD-1, or combination blockade. Methods: We performed a standardized meta-analysis of all available cohorts of ICB-treated (anti-PD-1 alone, anti-CTLA-4 alone, and combination) melanoma patients with sequencing of pretreatment tumor and clinical annotations (n=572; 268 evaluated post quality-control filtering, cutaneous). Raw data in the form of FASTQs were obtained for all cohorts, and data were reprocessed using standardized pipelines. Immune infiltrate was characterized using single-sample gene set enrichment analysis (ssGSEA) scores from 18 immune cell type signatures, and GSEA was performed on differentially expressed genes between responders and non-responders for each treatment cohort accounting for differences in immune infiltrate. To dissect the cell types and drivers of nominated signatures, we analyzed scRNAseq from a separate cohort of pretreatment metastatic melanoma. Results: We found that immune infiltrate was strongly predictive of response to combination blockade (OR=6.16, p Conclusion: High rates of response to combination in patients with high levels of immune infiltrate suggests that immune high patients may differentially benefit from combination blockade. Identifying features of non-responding patients stratified by immune infiltrate can lead to new biological insights into ICB. Citation Format: Amy Y. Huang, Natalie Vokes, Cora Ricker, Tyler Aprati, Emily Robitschek, Jiekun Yang, Li-Lun Ho, Kyriakitsa Galani, Kelly P. Burke, Giuseppe Tarantino, Jiajia Chen, Arlene H. Sharpe, Eliezer M. Van Allen, Manolis Kellis, Genevieve M. Boland, David Liu. Multiomic meta-analysis of differential response to PD-1 and CTLA-4 blockade in metastatic melanoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3271.

Immune checkpoint blockade (ICB) therapy has improved long-term survival for patients with advanced melanoma. However, there is critical need to identify potential biomarkers of response and actionable strategies to improve response rates. Through generation and analysis of over 200 chromatin modification maps for ICB-treated melanoma patient samples, melanoma cells and T cells, we identified significant enrichment of active enhancer states in non-responders at baseline. Enhancer mapping by bulk ChIP-Seq or single cell ATAC-Seq methods in two independent cohorts of ICB-treated melanoma samples identified an enhancer signature that predicted response to anti-PD-1 therapy. The activated non-responder enhancers marked a group of key regulators of several pathways in melanoma cells (including c-MET, TGFβ, EMT and AKT) that are known to mediate resistance to ICB therapy. In addition, several checkpoint receptors were alternatively activated by aberrant enhancers in T cells. Unbiased CRISPRi screening in melanoma cells and T cells identified novel functional enhancers, such as one upstream of c-MET, that mediate ICB response or T cell mediated killing. Finally, inhibition of enhancers and repression of these pathways using bromodomain inhibitors along with anti-PD-1 therapy significantly decreased melanoma tumor burden and increased T-cell infiltration. Epigenomic experiments identified a signature of 107 genes that could be used as pharmacodynamic marker for BET inhibitor response. Together, these findings identify enhancer upregulation as a key mechanism of adaptive resistance to ICB response, a potential enhancer-based biomarker of resistance to anti-PD-1 and enhancer blockade in combination with ICB as a potential strategy to improve responses. Citation Format: Kunal Rai, Mayinuer Maitituoheti, Alvin Shi, Ming Tang, Li-Lun Ho, Firas Youssef Kreidieh, Christopher Terranova, Kyriakitsa Galani, Emily Z. Keung, Caitlin A. Creasy, Manrong Wu, Jiajia Chen, Nana Chen, Anand K. Singh, Apoorvi Chaudhri, Nazanin E. Anvar, Giuseppe Tarrantino, Jiekun Yang, Sharmistha Sarkar, Shan Jiang, Jared Malke, Lauren Haydu, Elizabeth Burton, Michael A. Davies, Jeffrey E. Gershenwald, Patrick Hwu, Alexander Lazar, Jaime H. Cheah, Christian K. Soule, Stuart S. Levine, Chantale Bernatchez, Srinivas V. Saladi, David Liu, Hussein Tawbi, Jennifer Wargo, Genevieve M. Boland, Manolis Kellis. Functional roles of enhancers in immune microenvironment & immunotherapy response. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 6425.

Pancreatic ductal adenocarcinoma (PDAC) remains one of the most challenging cancer to treat. Due to the asymptomatic nature of the disease and ineffective drug treatment modalities, the survival rate of PDAC patients remains one of the lowest. The recurrent genetic alterations in PDAC are yet to be targeted; therefore, identifying effective therapeutic combinations is desperately needed. Here, we performed anin vivoCRISPR screening in a clinically relevant patient-derived xenograft (PDX) model system to identify synergistic drug combinations for PDAC treatment. Our approach revealed protein arginine methyltransferase gene 5 (PRMT5) as a promising druggable candidate whose inhibition creates synergistic vulnerability of PDAC cells to gemcitabine. Genetic and pharmacological inhibition results indicate that of PRMT5 depletion results in synergistic cytotoxicity with Gem due to depleted replication protein A (RPA) levels and an impaired non-homology end joining (NHEJ) DNA repair. Thus, the novel combination creates conditional lethality through the accumulation of excessive DNA damage and cell death, bothin vitroandin vivo. The findings demonstrate that unbiased genetic screenings combined with a clinically relevant model system is an effective approach in identifying synthetic lethal drug combinations for cancer treatment.STATEMENT of SIGNIFICANCEIdentify synergistic drug combinations for PDAC is a significant unmet need. Through CRISPR screening, we discovered and validated that PRMT5 depletion creates synergistic vulnerability of PDAC cells to gemcitabine. Mechanistically, the combination impairs DNA repair, synergistic accumulation of DNA damage and cell deathin vitroandin vivo.

Regular physical exercise has long been recognized to reverse the effects of diet-induced obesity, but the molecular mechanisms mediating these multi-tissue beneficial effects remain uncharacterized. Here, we address this challenge by studying the opposing effects of exercise training and high-fat diet at single-cell, deconvolution and tissue-level resolutions across 3 metabolic tissues. We profile scRNA-seq in 204,883 cells, grouped into 53 distinct cell subtypes/states in 22 major cell types, from subcuta-neous and visceral white adipose tissue (WAT), and skeletal muscle (SkM) in mice with diet and exercise training interventions. With a great number of mesenchymal stem cells (MSCs) profiled, we compared depot-specific adipose stem cell (ASC) states, and defined 7 distinct fibro-adipogenic progenitor (FAP) states in SkM including discovering and validating a novel CD140+/CD34+/SCA1-FAP population. Exercise- and obesity-regulated proportion, transcriptional and cell-cell interaction changes were most strongly pronounced in and centered around ASCs, FAPs, macrophages and T-cells. These changes reflected thermogenesis-vs-lipogenesis and hyperplasia-vs-hypertrophy shifts, clustered in pathways including extracellular matrix remodeling and circadian rhythm, and implicated complex single- and multi-tissue communication including training-associated shift of a cytokine from binding to its decoy receptor on ASCs to true receptor on M2 macrophages in vWAT. Overall, our work provides new insights on the metabolic protective effects of exercise training, uncovers a previously-underappreciated role of MSCs in mediating tissue-specific and multi-tissue effects, and serves as a model for multitissue single-cell analyses in physiologically complex and multifactorial traits exemplified by obesity and exercise training.

Pancreatic ductal adenocarcinoma (PDA) cells reprogram their transcriptional and metabolic programs to survive the nutrient-poor tumor microenvironment. Through in vivo CRISPR screening, we discovered islet-2 (ISL2) as a candidate tumor suppressor that modulates aggressive PDA growth. Notably, ISL2, a nuclear and chromatin-associated transcription factor, is epigenetically silenced in PDA tumors and high promoter DNA methylation or its reduced expression correlates with poor patient survival. The exogenous ISL2 expression or CRISPR-mediated upregulation of the endogenous loci reduces cell proliferation. Mechanistically, ISL2 regulates the expression of metabolic genes, and its depletion increases oxidative phosphorylation (OXPHOS). As such, ISL2-depleted human PDA cells are sensitive to the inhibitors of mitochondrial complex I in vitro and in vivo. Spatial transcriptomic analysis shows heterogeneous intratumoral ISL2 expression, which correlates with the expression of critical metabolic genes. These findings nominate ISL2 as a putative tumor suppressor whose inactivation leads to increased mitochondrial metabolism that may be exploitable therapeutically.

Exercise training is critical for the prevention and treatment of obesity, but its underlying mechanisms remain incompletely understood given the challenge of profiling heterogeneous effects across multiple tissues and cell types. Here, we address this challenge and opposing effects of exercise and high-fat diet (HFD)-induced obesity at single-cell resolution in subcutaneous and visceral white adipose tissue and skeletal muscle in mice with diet and exercise training interventions. We identify a prominent role of mesenchymal stem cells (MSCs) in obesity and exercise-induced tissue adaptation. Among the pathways regulated by exercise and HFD in MSCs across the three tissues, extracellular matrix remodeling and circadian rhythm are the most prominent. Inferred cell-cell interactions implicate within- and multi-tissue crosstalk centered around MSCs. Overall, our work reveals the intricacies and diversity of multi-tissue molecular responses to exercise and obesity and uncovers a previously underappreciated role of MSCs in tissue-specific and multi-tissue beneficial effects of exercise.

Pancreatic ductal adenocarcinoma (PDAC) remains one of the deadliest cancers. Uncovering mechanisms responsible for the heterogeneous clinical features of this disease is an essential step toward developing improved and more specific therapeutic approaches. Here, we sought to identify transcriptional regulators of aggressive PDAC growth through in vivo CRISPR screening of epigenetic and transcription factors in an orthotopic model. We identified the ISL LIM homeobox 2 (ISL2) gene as a tumor suppressor whose depletion enhances the proliferation of human PDAC cells in vitro and in vivo and cooperates with activated KRAS to initiate PDAC in a murine model. Conversely, the upregulation of ISL2 expression through CRISPR-mediated locus-specific epigenetic editing results in reduced cell proliferation. Importantly, ISL2 is epigenetically silenced through DNA methylation in ~60% of PDAC tumors, which correlates with poor patient outcome. Functional studies showed that ISL2 loss rewires metabolic gene expression, and consequently potentiates oxidative phosphorylation while reducing glycolysis. This metabolic shift creates selective vulnerability to small molecule inhibitors of mitochondrial respiration and fatty acid oxidation. Collectively, these findings reveal ISL2 as a novel tumor suppressor whose inactivation drives metabolic reprogramming in an aggressive PDAC subset and point to potential therapeutic vulnerabilities in these tumors.

Predicting the response and identifying additional targets that will improve the efficacy of chemotherapy is a major goal in cancer research. Through large-scale in vivo and in vitro CRISPR knockout screens in pancreatic ductal adenocarcinoma cells, we identified genes whose genetic deletion or pharmacologic inhibition synergistically increase the cytotoxicity of MEK signaling inhibitors. Furthermore, we show that CRISPR viability scores combined with basal gene expression levels could model global cellular responses to the drug treatment. We develop drug response evaluation by in vivo CRISPR screening (DREBIC) method and validated its efficacy using large-scale experimental data from independent experiments. Comparative analyses demonstrate that DREBIC predicts drug response in cancer cells from a wide range of tissues with high accuracy and identifies therapeutic vulnerabilities of cancer-causing mutations to MEK inhibitors in various cancer types., Predicting the response to chemotherapy is a major goal of cancer research. Here the authors use CRISPR knockout screens in pancreatic ductal adenocarcinoma cells to identify deletions synergistic with MEK inhibitors.

The Aerhada Pb-Zn-Ag deposit is located in the western segment of the Great Hinggan Range Ag-Pb-Zn-Cu-Mo-Au-Fe metallogenic belt in NE China. Orebodies occur mainly as vein type and are hosted by sandstone and siliceous slate. Three stages of primary mineralization, including an early arsenopyrite-pyrite-quartz, a middle polymetallic and silver sulfides-quartz and a late sphalerite-pyrite-calcite-fluorite are recognized. Four types of fluid inclusions have been identified in the ore-bearing quartz and fluorite veins, i.e., liquid-rich, gas-rich, three-phase CO 2 aqueous inclusions, and pure gas or liquid aqueous inclusions. Microthermometric studies on fluid inclusions reveal that homogenization temperatures from early to late stages range from 253° to 430 °C, 195° to 394 °C and 133° to 207 °C, respectively. Fluid salinities range from 2.9 to 14.0 wt.% NaCl equiv. The vapor composition of the ore fluid is dominated by H 2 O, CO 2 and CH 4 , with minor proportions of N 2 . The fluid δ 18 O H2O and δ D H2O values vary from +1.6 to +9.3‰ and −122 to −56‰, respectively, and reflect a magmatic fluid and a meteoric fluid dominant hydrothermal system for the early and late stages of mineralization, respectively. The calculated δ 34 S H2S values of hydrothermal fluids in equilibrium with sulfides range from +5.2 to +7.1‰, suggesting a mixed source for sulfur, i.e., the local magmatic and sedimentary rocks. The Pb isotope compositions of sulfides are similar to those of the local magmatic and sedimentary rocks, implying that lead and possibly silver relate to these sources. The noble gas isotope compositions of fluid inclusions hosted in ore minerals suggest that the ore-forming fluids were dominantly derived from a deep mantle source. Fluid mixing and dilution are inferred as the dominant mechanisms for ore deposition. The Aerhada Pb-Zn-Ag deposit can be classified as a medium to low temperature hydrothermal vein type deposit.

Chemoresistance is driven by unique regulatory networks in the genome that are distinct from those necessary for cancer development. Here, we investigate the contribution of enhancer elements to cisplatin resistance in ovarian cancers. Epigenome profiling of multiple cellular models of chemoresistance identified unique sets of distal enhancers, super-enhancers (SE), and their gene targets that coordinate and maintain the transcriptional program of the platinum-resistant state in ovarian cancer. Pharmacologic inhibition of distal enhancers through small-molecule epigenetic inhibitors suppressed the expression of their target genes and restored cisplatin sensitivity in vitro and in vivo. In addition to known drivers of chemoresistance, our findings identified SOX9 as a critical SE-regulated transcription factor that plays a critical role in acquiring and maintaining the chemoresistant state in ovarian cancer. The approach and findings presented here suggest that integrative analysis of epigenome and transcriptional programs could identify targetable key drivers of chemoresistance in cancers. Significance: Integrative genome-wide epigenomic and transcriptomic analyses of platinum-sensitive and -resistant ovarian lines identify key distal regulatory regions and associated master regulator transcription factors that can be targeted by small-molecule epigenetic inhibitors.

Whole-genome sequencing efforts of tumors and normal tissues have identified numerous genetic mutations, both somatic and germline, that do not overlap with coding genomic sequences. Attributing a functional role to these noncoding mutations and characterizing them using experimental methods has been more challenging compared with coding mutations. In this review, we provide a brief introduction to the world of noncoding mutations. We discuss recent progress in identifying noncoding mutations and the analytic and experimental approaches utilized to interpret their functional roles. We also highlight the potential mechanisms by which a noncoding mutation may exert its effect and discuss future challenges and opportunities.

Background The mutational processes underlying non-coding cancer mutations and their biological significance in tumor evolution are poorly understood. To get better insights into the biological mechanisms of mutational processes in breast cancer, we integrate whole-genome level somatic mutations from breast cancer patients with chromatin states and transcription factor binding events. Results We discover that a large fraction of non-coding somatic mutations in estrogen receptor (ER)-positive breast cancers are confined to ER binding sites. Notably, the highly mutated estrogen receptor binding sites are associated with more frequent chromatin loop contacts and the associated distal genes are expressed at higher level. To elucidate the functional significance of these non-coding mutations, we focus on two of the recurrently mutated estrogen receptor binding sites. Our bioinformatics and biochemical analysis suggest loss of DNA-protein interactions due to the recurrent mutations. Through CRISPR interference, we find that the recurrently mutated regulatory element at the LRRC3C-GSDMA locus impacts the expression of multiple distal genes. Using a CRISPR base editor, we show that the recurrent C→T conversion at the ZNF143 locus results in decreased TF binding, increased chromatin loop formation, and increased expression of multiple distal genes. This single point mutation mediates reduced response to estradiol-induced cell proliferation but increased resistance to tamoxifen-induced growth inhibition. Conclusions Our data suggest that ER binding is associated with localized accumulation of somatic mutations, some of which affect chromatin architecture, distal gene expression, and cellular phenotypes in ER-positive breast cancer. Electronic supplementary material The online version of this article (10.1186/s13059-018-1572-4) contains supplementary material, which is available to authorized users.

Review history. (DOCX 58 kb)

Figure S1. ERBS shared by more patients contain more mutations after controlling for ER binding intensity. Figure S2. Comparable levels of somatic mutation burden are observed at ERBS overlapping promoter, intronic, and intergenic regions. Figure S3. ERBS are protected from somatic insertions and deletions, and the protective effect is correlated with ER binding intensity. Figure S4. Somatic mutations enriched at ERBS show the APOBEC mutational signature. Figure S5. ERBS with more mutations make more frequent chromatin interactions independent of ER binding intensity. Figure S6. ERBS associated with poor/met outcome contain more somatic mutations. Figure S7. Amplification plots for WT and mutant MCF-7 clones based on the qPCR colony screening strategy depicted in Figure 4d. Figure S8. Somatic mutation reduces ZBTB7A binding. ChIP-qPCR analysis shows ZBTB7A enrichment at the mutation site in MCF-7 WT cells and a mutant clone. Figure S9. High expression of TMEM41B, IPO7 and WEE1 is associated with poor survival for breast cancer patients. Figure S10. Somatic mutation burden at ERBS is higher when blood instead of adjacent to tumor breast tissue is used as â normalâ in the mutation calling process. (PDF 2338 kb)

CRISPR-Cas9-induced DNA damage may have deleterious effects at high-copy-number genomic regions. Here, we use CRISPR base editors to knock out genes by changing single nucleotides to create stop codons. We show that the CRISPR-STOP method is an efficient and less deleterious alternative to wild-type Cas9 for gene-knockout studies. Early stop codons can be introduced in ∼17,000 human genes. CRISPR-STOP-mediated targeted screening demonstrates comparable efficiency to WT Cas9, which indicates the suitability of our approach for genome-wide functional screenings.

As part of a larger study on acculturation information seeking (AIS), this poster presents preliminary findings of a survey of 188 international students in the United States. The study surveyed the types of information needs related to academic and everyday life, those that were challenging to fulfill, and the sources used for AIS. Using the Big Five Inventory, the study also tested the relationship between AIS and personality. Findings show that individuals with different personality traits viewed different needs (e.g., health, legal information) as challenging. Similarly, respondents with certain traits favored certain sources that were seldom used by others (e.g., those with high openness use social Q&A sites more frequently). Recommendations for better sharing of acculturation information are discussed.

Experimental approaches to genetic studies of complex traits evolve with technological advances. How do discoveries using different approaches advance our knowledge of the genetic architecture underlying complex diseases/traits? Do most of the findings of newer techniques, such as genome-wide association study (GWAS), provide more information than older ones, for example, genome-wide linkage study? In this review, we address these issues by developing a nicotine dependence (ND) genetic susceptibility map based on the results obtained by the approaches commonly used in recent years, namely, genome-wide linkage, candidate gene association, GWAS and targeted sequencing. Converging and diverging results from these empirical approaches have elucidated a preliminary genetic architecture of this intractable psychiatric disorder and yielded new hypotheses on ND etiology. The insights we obtained by putting together results from diverse approaches can be applied to other complex diseases/traits. In sum, developing a genetic susceptibility map and keeping it updated are effective ways to keep track of what we know about a disease/trait and what the next steps may be with new approaches.

Significant association of CHRNB3 variants with nicotine dependence in multiple ethnic populations

Previous studies have implicated genes encoding the 5-HT3AB receptors (HTR3A and HTR3B) and the serotonin transporter (SLC6A4), both independently and interactively, in alcohol (AD), cocaine (CD), and nicotine dependence (ND). However, whether these genetic effects also exist in subjects with comorbidities remains largely unknown. We used 1,136 African-American (AA) and 2,428 European-American (EA) subjects from the Study of Addiction: Genetics and Environment (SAGE) to determine associations between 88 genotyped or imputed variants within HTR3A, HTR3B, and SLC6A4 and three types of addictions, which were measured by DSM-IV diagnoses of AD, CD, and ND and the Fagerstrom Test for Nicotine Dependence (FTND), an independent measure of ND commonly used in tobacco research. Individual SNP-based association analysis revealed a significant association of rs2066713 in SLC6A4 with FTND in AA (β = −1.39; P = 1.6E − 04). Haplotype-based association analysis found one major haplotype formed by SNPs rs3891484 and rs3758987 in HTR3B that was significantly associated with AD in the AA sample, and another major haplotype T–T-G, formed by SNPs rs7118530, rs12221649, and rs2085421 in HTR3A, which showed significant association with FTND in the EA sample. Considering the biologic roles of the three genes and their functional relations, we used the GPU-based Generalized Multifactor Dimensionality Reduction (GMDR-GPU) program to test SNP-by-SNP interactions within the three genes and discovered two- to five-variant models that have significant impacts on AD, CD, ND, or FTND. Interestingly, most of the SNPs included in the genetic interaction model(s) for each addictive phenotype are either overlapped or in high linkage disequilibrium for both AA and EA samples, suggesting these detected variants in HTR3A, HTR3B, and SLC6A4 are interactively contributing to etiology of the three addictive phenotypes examined in this study.