Multi-ancestry meta-analysis of tobacco use disorder identifies 461 potential risk genes and reveals associations with multiple health outcomes.

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviours and although strides have been made using genome-wide association studies to identify risk variants, most variants identified have been for nicotine consumption, rather than TUD. Here we leveraged four US biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records) in 653,790 individuals (495,005 European, 114,420 African American and 44,365 Latin American) and data from UK Biobank (ncombined = 898,680). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviours in children and hundreds of medical outcomes, including HIV infection, heart disease and pain. This work furthers our biological understanding of TUD and establishes electronic health records as a source of phenotypic information for studying the genetics of TUD.

Myocardial Recovery in Recent Onset Dilated Cardiomyopathy: Role of CDCP1 and Cardiac Fibrosis.

BACKGROUND: Dilated cardiomyopathy (DCM) is a major cause of heart failure and carries a high mortality rate. Myocardial recovery in DCM-related heart failure patients is highly variable, with some patients having little or no response to standard drug therapy. A genome-wide association study may agnostically identify biomarkers and provide novel insight into the biology of myocardial recovery in DCM. METHODS: A genome-wide association study for change in left ventricular ejection fraction was performed in 686 White subjects with recent-onset DCM who received standard pharmacotherapy. Genome-wide association study signals were subsequently functionally validated and studied in relevant cellular models to understand molecular mechanisms that may have contributed to the change in left ventricular ejection fraction. RESULTS: The genome-wide association study identified a highly suggestive locus that mapped to the 5'-flanking region of the CDCP1 (CUB [complement C1r/C1s, Uegf, and Bmp1] domain containing protein 1) gene (rs6773435; P=7.12×10-7). The variant allele was associated with improved cardiac function and decreased CDCP1 transcription. CDCP1 expression was significantly upregulated in human cardiac fibroblasts (HCFs) in response to the PDGF (platelet-derived growth factor) signaling, and knockdown of CDCP1 significantly repressed HCF proliferation and decreased AKT (protein kinase B) phosphorylation. Transcriptomic profiling after CDCP1 knockdown in HCFs supported the conclusion that CDCP1 regulates HCF proliferation and mitosis. In addition, CDCP1 knockdown in HCFs resulted in significantly decreased expression of soluble ST2 (suppression of tumorigenicity-2), a prognostic biomarker for heart failure and inductor of cardiac fibrosis. CONCLUSIONS: CDCP1 may play an important role in myocardial recovery in recent-onset DCM and mediates its effect primarily by attenuating cardiac fibrosis.

Multi-ancestry meta-analysis of tobacco use disorder prioritizes novel candidate risk genes and reveals associations with numerous health outcomes.

Tobacco use disorder (TUD) is the most prevalent substance use disorder in the world. Genetic factors influence smoking behaviors, and although strides have been made using genome-wide association studies (GWAS) to identify risk variants, the majority of variants identified have been for nicotine consumption, rather than TUD. We leveraged five biobanks to perform a multi-ancestral meta-analysis of TUD (derived via electronic health records, EHR) in 898,680 individuals (739,895 European, 114,420 African American, 44,365 Latin American). We identified 88 independent risk loci; integration with functional genomic tools uncovered 461 potential risk genes, primarily expressed in the brain. TUD was genetically correlated with smoking and psychiatric traits from traditionally ascertained cohorts, externalizing behaviors in children, and hundreds of medical outcomes, including HIV infection, heart disease, and pain. This work furthers our biological understanding of TUD and establishes EHR as a source of phenotypic information for studying the genetics of TUD.

Polygenic risk for prostate cancer: Decreasing relative risk with age but little impact on absolute risk.

Polygenic risk scores (PRSs) for a variety of diseases have recently been shown to have relative risks that depend on age, and genetic relative risks decrease with increasing age. A refined understanding of the age dependency of PRSs for a disease is important for personalized risk predictions and risk stratification. To further evaluate how the PRS relative risk for prostate cancer depends on age, we refined analyses for a validated PRS for prostate cancer by using 64,274 prostate cancer cases and 46,432 controls of diverse ancestry (82.8% European, 9.8% African American, 3.8% Latino, 2.8% Asian, and 0.8% Ghanaian). Our strategy applied a novel weighted proportional hazards model to case-control data to fully utilize age to refine how the relative risk decreased with age. We found significantly greater relative risks for younger men (age 30-55 years) compared with older men (70-88 years) for both relative risk per standard deviation of the PRS and dichotomized according to the upper 90th percentile of the PRS distribution. For the largest European ancestral group that could provide reliable resolution, the log-relative risk decreased approximately linearly from age 50 to age 75. Despite strong evidence of age-dependent genetic relative risk, our results suggest that absolute risk predictions differed little from predictions that assumed a constant relative risk over ages, from short-term to long-term predictions, simplifying implementation of risk discussions into clinical practice.

Antimullerian Hormone as a Serum Biomarker for Risk of Chemotherapy-Induced Amenorrhea.

Antimullerian hormone (AMH) is a promising biomarker for ovarian reserve. In this study, we assessed AMH before and 1 year after initiation of adjuvant chemotherapy on National Surgical Adjuvant Breast and Bowel Project (NSABP)/NRG Oncology B-47 in female participants aged 42 years and younger (median age = 39 years). At baseline, median AMH was 1.2 ng/mL; 13 (4.7%) values were less than 0.1 ng/mL (the threshold for detectable levels, in the perimenopause and menopause range), and 57 values (20.6%) were less than 0.5 ng/mL. At 1 year, 215 (77.6%) were less than 0.1 ng/mL, and 264 (95.3%) were less than 0.5 ng/mL. Postchemotherapy menses were reported by 46.2% of participants. Multivariable logistic regression found that the odds of having postchemotherapy menses increased with younger age, higher body mass index, and higher prechemotherapy AMH, but not by trastuzumab administration or by the choice of chemotherapy (doxorubicin-cyclophosphamide followed by paclitaxel vs docetaxel-cyclophosphamide). We conclude that higher prechemotherapy AMH predicts a lower risk of chemotherapy-induced amenorrhea and that AMH 1 year after chemotherapy initiation is not informative in this setting because it is likely to be very low.

Adult diffuse glioma GWAS by molecular subtype identifies variants in D2HGDH and FAM20C.

BACKGROUND: Twenty-five germline variants are associated with adult diffuse glioma, and some of these variants have been shown to be associated with particular subtypes of glioma. We hypothesized that additional germline variants could be identified if a genome-wide association study (GWAS) were performed by molecular subtype. METHODS: A total of 1320 glioma cases and 1889 controls were used in the discovery set and 799 glioma cases and 808 controls in the validation set. Glioma cases were classified into molecular subtypes based on combinations of isocitrate dehydrogenase (IDH) mutation, telomerase reverse transcriptase (TERT) promoter mutation, and 1p/19q codeletion. Logistic regression was applied to the discovery and validation sets to test for associations of variants with each of the subtypes. A meta-analysis was subsequently performed using a genome-wide P-value threshold of 5 × 10-8. RESULTS: Nine variants in or near D-2-hydroxyglutarate dehydrogenase (D2HGDH) on chromosome 2 were genome-wide significant in IDH-mutated glioma (most significant was rs5839764, meta P = 2.82 × 10-10). Further stratifying by 1p/19q codeletion status, one variant in D2HGDH was genome-wide significant in IDH-mutated non-codeleted glioma (rs1106639, meta P = 4.96 × 10-8). Further stratifying by TERT mutation, one variant near FAM20C (family with sequence similarity 20, member C) on chromosome 7 was genome-wide significant in gliomas that have IDH mutation, TERT mutation, and 1p/19q codeletion (rs111976262, meta P = 9.56 × 10-9). Thirty-six variants in or near GMEB2 on chromosome 20 near regulator of telomere elongation helicase 1 (RTEL1) were genome-wide significant in IDH wild-type glioma (most significant was rs4809313, meta P = 2.60 × 10-10). CONCLUSIONS: Performing a GWAS by molecular subtype identified 2 new regions and a candidate independent region near RTEL1, which were associated with specific glioma molecular subtypes.

Genetic predictors of chemotherapy-related amenorrhea in women with breast cancer.

OBJECTIVE: To study how genetics may play a role in determining risk of chemotherapy-related amenorrhea (CRA) in young women with breast cancer. DESIGN: Genome-wide association study. SETTING: Not applicable. PATIENT(S): Premenopausal women ≤45 years of age enrolled in one of these three trials were included if they had at least one menstrual case report form after chemotherapy ended and if they were of European ancestry. Forms during and up to 3 months after receipt of GnRH agonist were excluded. INTERVENTION(S): None. MAIN OUTCOME MEASURE(S): The association of single-nucleotide polymorphisms with post-chemotherapy menstruation adjusted for trial and arm, age, tamoxifen use, and nodal status. RESULT(S): The median age of the 1,168 women was 41 years (range 19-45). Among these, 457 (39%) never resumed menses after chemotherapy. Older age, tamoxifen use, and node-negative disease were associated with increased risk of CRA. Adjusting for these, rs147451859, in an intron of PPCDC (phosphopantothenoylcysteine decarboxylase), and rs17587029, located 5' upstream of RPS20P11 (ribosomal protein S20 pseudogene 11), were associated with post-chemotherapy menstruation. CONCLUSION(S): Genetic variation may contribute to risk of CRA. Better prediction of who will experience CRA may inform reproductive and treatment decision making in young women with cancer.

Soluble Neprilysin in the General Population: Clinical Determinants and Its Relationship to Cardiovascular Disease.

Background Neprilysin is a metalloprotease involved in proteolysis of numerous peptides, including natriuretic peptides, and is of prognostic and therapeutic importance in heart failure with reduced ejection fraction. No studies have investigated circulating neprilysin in the community, its clinical correlates, or its relationship to cardiovascular disease in the general population. Methods and Results Plasma neprilysin was measured in 1536 participants from Olmsted County, Minnesota, using a commercially available sandwich ELISA assay. Clinical and echocardiographic correlates and subsequent outcomes were determined. Soluble neprilysin is non-normally distributed in the community (median: 3.9 ng/mL; interquartile range: 1.0-43.0 ng/mL). There was no relationship between plasma neprilysin and age (Spearman correlation: -0.04, P=0.16); body mass index (Spearman correlation: -0.04, P=0.16); glomerular filtration rate (Spearman correlation: -0.007, P=0.8); or A-, B-, or C-type natriuretic peptides (Spearman correlation: 0.03, P=0.22; -0.001, P=0.96; 0.01, P=0.67, respectively). Among tertiles of neprilysin, the lowest tertile group had the highest prevalence of smokers (P<0.001), hypertension (P=0.04), dyslipidemia (P=0.03), and diastolic dysfunction (P=0.02). Soluble neprilysin was not prospectively associated with death or heart failure over a median of 10.7 years. Conclusions In a large community-based cohort, for the first time, we described the distribution of circulating neprilysin in the general community. We observed that neprilysin does not correlate with natriuretic peptide levels and is not independently associated with adverse outcomes. The novel associations observed between low soluble neprilysin levels and an adverse cardiometabolic and smoking profile requires further investigation.

Anastrozole Aromatase Inhibitor Plasma Drug Concentration Genome-Wide Association Study: Functional Epistatic Interaction Between SLC38A7 and ALPPL2.

Anastrozole is a widely prescribed aromatase inhibitor for the therapy of estrogen receptor positive (ER+) breast cancer. We performed a genome-wide association study (GWAS) for plasma anastrozole concentrations in 687 postmenopausal women with ER+ breast cancer. The top single-nucleotide polymorphism (SNP) signal mapped across SLC38A7 (rs11648166, P = 2.3E-08), which we showed to encode an anastrozole influx transporter. The second most significant signal (rs28845026, P = 5.4E-08) mapped near ALPPL2 and displayed epistasis with the SLC38A7 signal. Both of these SNPs were cis expression quantitative trait loci (eQTL)s for these genes, and patients homozygous for variant genotypes for both SNPs had the highest drug concentrations, the highest SLC38A7 expression, and the lowest ALPPL2 expression. In summary, our GWAS identified a novel gene encoding an anastrozole transporter, SLC38A7, as well as epistatic interaction between SNPs in that gene and SNPs near ALPPL2 that influenced both the expression of the transporter and anastrozole plasma concentrations.

SLCO1B1 polymorphisms and plasma estrone conjugates in postmenopausal women with ER+ breast cancer: genome-wide association studies of the estrone pathway.

BACKGROUND: Estrone (E1), the major circulating estrogen in postmenopausal women, promotes estrogen-receptor positive (ER+) breast tumor growth and proliferation. Two major reactions contribute to E1 plasma concentrations, aromatase (CYP19A1) catalyzed E1 synthesis from androstenedione and steroid sulfatase (STS) catalyzed hydrolysis of estrone conjugates (E1Cs). E1Cs have been associated with breast cancer risk and may contribute to tumor progression since STS is expressed in breast cancer where its activity exceeds that of aromatase. METHODS: We performed genome-wide association studies (GWAS) to identify SNPs associated with variation in plasma concentrations of E1Cs, E1, and androstenedione in 774 postmenopausal women with resected early-stage ER+ breast cancer. Hormone concentrations were measured prior to aromatase inhibitor therapy. RESULTS: Multiple SNPs in SLCO1B1, a gene encoding a hepatic influx transporter, displayed genome-wide significant associations with E1C plasma concentrations and with the E1C/E1 ratio. The top SNP for E1C concentrations, rs4149056 (p = 3.74E-11), was a missense variant that results in reduced transporter activity. Patients homozygous for the variant allele had significantly higher average E1C plasma concentrations than did other patients. Furthermore, three other SLCO1B1 SNPs, not in LD with rs4149056, were associated with both E1C concentrations and the E1C/E1 ratio and were cis-eQTLs for SLCO1B3. GWAS signals of suggestive significance were also observed for E1, androstenedione, and the E1/androstenedione ratio. CONCLUSION: These results suggest a mechanism for genetic variation in E1C plasma concentrations as well as possible SNP biomarkers to identify ER+ breast cancer patients for whom STS inhibitors might be of clinical value.

Metformin pharmacogenomics: a genome-wide association study to identify genetic and epigenetic biomarkers involved in metformin anticancer response using human lymphoblastoid cell lines.

Metformin is currently considered as a promising anticancer agent in addition to its anti-diabetic effect. To better individualize metformin therapy and explore novel molecular mechanisms in cancer treatment, we conducted a pharmacogenomic study using 266 lymphoblastoid cell lines (LCLs). Metformin cytotoxicity assay was performed using the MTS assay. Genome-wide association (GWA) analyses were performed in LCLs using 1.3 million SNPs, 485k DNA methylation probes, 54k mRNA expression probe sets, and metformin cytotoxicity (IC50s). Top candidate genes were functionally validated using siRNA screening, followed by MTS assay in breast cancer cell lines. Further study of one top candidate, STUB1, was performed to elucidate the mechanisms by which STUB1 might contribute to metformin action. GWA analyses in LCLs identified 198 mRNA expression probe sets, 12 SNP loci, and 5 DNA methylation loci associated with metformin IC50 with P-values <10−4 or <10−5. Integrated SNP/methylation loci-expression-IC50 analyses found 3 SNP loci or 5 DNA methylation loci associated with metformin IC50 through trans-regulation of expression of 11 or 26 genes with P-value <10−4. Functional validation of top 61 candidate genes in 4 IPA networks indicated down regulation of 14 genes significantly altered metformin sensitivity in two breast cancer cell lines. Mechanistic studies revealed that the E3 ubiquitin ligase, STUB1, could influence metformin response by facilitating proteasome-mediated degradation of cyclin A. GWAS using a genomic data-enriched LCL model system, together with functional and mechanistic studies using cancer cell lines, help us to identify novel genetic and epigenetic biomarkers involved in metformin anticancer response.

Aromatase inhibitor-associated bone fractures: a case-cohort GWAS and functional genomics.

Bone fractures are a major consequence of osteoporosis. There is a direct relationship between serum estrogen concentrations and osteoporosis risk. Aromatase inhibitors (AIs) greatly decrease serum estrogen levels in postmenopausal women, and increased incidence of fractures is a side effect of AI therapy. We performed a discovery case-cohort genome-wide association study (GWAS) using samples from 1071 patients, 231 cases and 840 controls, enrolled in the MA.27 breast cancer AI trial to identify genetic factors involved in AI-related fractures, followed by functional genomic validation. Association analyses identified 20 GWAS single nucleotide polymorphism (SNP) signals with P < 5E-06. After removal of signals in gene deserts and those composed entirely of imputed SNPs, we applied a functional validation "decision cascade" that resulted in validation of the CTSZ-SLMO2-ATP5E, TRAM2-TMEM14A, and MAP4K4 genes. These genes all displayed estradiol (E2)-dependent induction in human fetal osteoblasts transfected with estrogen receptor-α, and their knockdown altered the expression of known osteoporosis-related genes. These same genes also displayed SNP-dependent variation in E2 induction that paralleled the SNP-dependent induction of known osteoporosis genes, such as osteoprotegerin. In summary, our case-cohort GWAS identified SNPs in or near CTSZ-SLMO2-ATP5E, TRAM2-TMEM14A, and MAP4K4 that were associated with risk for bone fracture in estrogen receptor-positive breast cancer patients treated with AIs. These genes displayed E2-dependent induction, their knockdown altered the expression of genes related to osteoporosis, and they displayed SNP genotype-dependent variation in E2 induction. These observations may lead to the identification of novel mechanisms associated with fracture risk in postmenopausal women treated with AIs.

Thiopurine pharmacogenomics: association of SNPs with clinical response and functional validation of candidate genes.

AIM: We investigated candidate genes associated with thiopurine metabolism and clinical response in childhood acute lymphoblastic leukemia. MATERIALS & METHODS: We performed genome-wide SNP association studies of 6-thioguanine and 6-mercaptopurine cytotoxicity using lymphoblastoid cell lines. We then genotyped the top SNPs associated with lymphoblastoid cell line cytotoxicity, together with tagSNPs for genes in the 'thiopurine pathway' (686 total SNPs), in DNA from 589 Caucasian UK ALL97 patients. Functional validation studies were performed by siRNA knockdown in cancer cell lines. RESULTS: SNPs in the thiopurine pathway genes ABCC4, ABCC5, IMPDH1, ITPA, SLC28A3 and XDH, and SNPs located within or near ATP6AP2, FRMD4B, GNG2, KCNMA1 and NME1, were associated with clinical response and measures of thiopurine metabolism. Functional validation showed shifts in cytotoxicity for these genes. CONCLUSION: The clinical response to thiopurines may be regulated by variation in known thiopurine pathway genes and additional novel genes outside of the thiopurine pathway.

Discovery of genetic biomarkers contributing to variation in drug response of cytidine analogues using human lymphoblastoid cell lines.

BACKGROUND: Two cytidine analogues, gemcitabine and cytosine arabinoside (AraC), are widely used in the treatment of a variety of cancers with a large individual variation in response. To identify potential genetic biomarkers associated with response to these two drugs, we used a human lymphoblastoid cell line (LCL) model system with extensive genomic data, including 1.3 million SNPs and 54,000 basal expression probesets to perform genome-wide association studies (GWAS) with gemcitabine and AraC IC50 values. RESULTS: We identified 11 and 27 SNP loci significantly associated with gemcitabine and AraC IC50 values, respectively. Eleven candidate genes were functionally validated using siRNA knockdown approach in multiple cancer cell lines. We also characterized the potential mechanisms of genes by determining their influence on the activity of 10 cancer-related signaling pathways using reporter gene assays. Most SNPs regulated gene expression in a trans manner, except 7 SNPs in the PIGB gene that were significantly associated with both the expression of PIGB and gemcitabine cytotoxicity. CONCLUSION: These results suggest that genetic variation might contribute to drug response via either cis- or trans- regulation of gene expression. GWAS analysis followed by functional pharmacogenomics studies might help identify novel biomarkers contributing to variation in response to these two drugs and enhance our understanding of underlying mechanisms of drug action.

Integrative gene set analysis: application to platinum pharmacogenomics.

Integrative genomics has the potential to uncover relevant loci, as clinical outcome and response to chemotherapies are most likely not due to a single gene (or data type) but rather a complex relationship involving genetic variation, mRNA, DNA methylation, and copy number variation. In addition to this complexity, many complex phenotypes are thought to be controlled by the interplay of multiple genes within the same molecular pathway or gene set (GS). To address these two challenges, we propose an integrative gene set analysis approach and apply this strategy to a cisplatin (CDDP) pharmacogenomics study involving lymphoblastoid cell lines for which genome-wide SNP and mRNA expression data was collected. Application of the integrative GS analysis implicated the role of the RNA binding and cytoskeletal part GSs. The genes LMNB1 and CENPF, within the cytoskeletal part GS, were functionally validated with siRNA knockdown experiments, where the knockdown of LMNB1 and CENPF resulted in CDDP resistance in multiple cancer cell lines. This study demonstrates the utility of an integrative GS analysis strategy for detecting novel genes associated with response to cancer therapies, moving closer to tailored therapy decisions for cancer patients.

Identification of genetic variants or genes that are associated with Homoharringtonine (HHT) response through a genome-wide association study in human lymphoblastoid cell lines (LCLs).

Homoharringtonine (HHT) has been widely used in China to treat patients with acute and chronic myeloid leukemia for decades. Since response to HHT varies among patients, our study aimed to identify biomarkers that might influence the response to HHT using a panel of various human lymphoblastoid cell lines (LCLs). Genome-wide association (GWA) analysis using single nucleotide polymorphism (SNP) and mRNA expression data was assessed for association with cytotoxicity to HHT in LCLs. Integrated analysis among SNPs, expression, AUC value was also performed to help select candidate genes for further functional characterization. Functional validation of candidate genes was performed using leukemia cell lines (U937, K562). Candidate genes were knocked down using specific siRNA and its response to HHT was assessed using MTS assay. We found that 15 expression probes were associated with HHT AUC with P < 10(-4), and 96 individual probe sets with P < 10(-3). Eighteen SNPs were associated with HHT AUC with P < 10(-5) and 281 SNPs with P < 10(-4). The integrated analysis identified 4 unique SNPs that were associated with both expression and AUC. Functional validation using siRNA knockdown in leukemia cell lines showed that knocking down CCDC88A, CTBP2, SOCS4 genes in U937 and K562 cells significantly altered HHT cytotoxicity. In summary, this study performed with LCLs can help to identify novel biomarker that might contribute to variation in response to HHT therapy.

Genome-wide association study for biomarker identification of Rapamycin and Everolimus using a lymphoblastoid cell line system.

The mammalian target of rapamycin (mTOR) inhibitors, a set of promising potential anti-cancer agents, has shown response variability among individuals. This study aimed to identify novel biomarkers and mechanisms that might influence the response to Rapamycin and Everolimus. Genome-wide association (GWA) analyses involving single nucleotide polymorphisms (SNPs), mRNA, and microRNAs microarray data were assessed for association with area under the cytotoxicity dose response curve (AUC) of two mTOR inhibitors in 272 human lymphoblastoid cell lines (LCLs). Integrated analysis among SNPs, expression data, microRNA data and AUC values were also performed to help select candidate genes for further functional characterization. Functional validation of candidate genes using siRNA screening in multiple cell lines followed by MTS assays for the two mTOR inhibitors were performed. We found that 16 expression probe sets (genes) that overlapped between the two drugs were associated with AUC values of two mTOR inhibitors. One hundred and twenty seven and one hundred SNPs had P < 10(-4), while 8 and 10 SNPs had P < 10(-5) with Rapamycin and Everolimus AUC, respectively. Functional studies indicated that 13 genes significantly altered cell sensitivity to either one or both drugs in at least one cell line. Additionally, one microRNA, miR-10a, was significantly associated with AUC values for both drugs and was shown to repress expression of genes that were associated with AUC and desensitize cells to both drugs. In summary, this study identified genes and a microRNA that might contribute to response to mTOR inhibitors.

Selective estrogen receptor modulators and pharmacogenomic variation in ZNF423 regulation of BRCA1 expression: individualized breast cancer prevention.

The selective estrogen receptor modulators (SERM) tamoxifen and raloxifene can reduce the occurrence of breast cancer in high-risk women by 50%, but this U.S. Food and Drug Administration-approved prevention therapy is not often used. We attempted to identify genetic factors that contribute to variation in SERM breast cancer prevention, using DNA from the NSABP P-1 and P-2 breast cancer prevention trials. An initial discovery genome-wide association study identified common single-nucleotide polymorphisms (SNP) in or near the ZNF423 and CTSO genes that were associated with breast cancer risk during SERM therapy. We then showed that both ZNF423 and CTSO participated in the estrogen-dependent induction of BRCA1 expression, in both cases with SNP-dependent variation in induction. ZNF423 appeared to be an estrogen-inducible BRCA1 transcription factor. The OR for differences in breast cancer risk during SERM therapy for subjects homozygous for both protective or both risk alleles for ZNF423 and CTSO was 5.71.

TSPYL5 SNPs: association with plasma estradiol concentrations and aromatase expression.

We performed a discovery genome-wide association study to identify genetic factors associated with variation in plasma estradiol (E2) concentrations using DNA from 772 postmenopausal women with estrogen receptor (ER)-positive breast cancer prior to the initiation of aromatase inhibitor therapy. Association analyses showed that the single nucleotide polymorphisms (SNP) (rs1864729) with the lowest P value (P = 3.49E-08), mapped to chromosome 8 near TSPYL5. We also identified 17 imputed SNPs in or near TSPYL5 with P values < 5E-08, one of which, rs2583506, created a functional estrogen response element. We then used a panel of lymphoblastoid cell lines (LCLs) stably transfected with ERα with known genome-wide SNP genotypes to demonstrate that TSPYL5 expression increased after E2 exposure of cells heterozygous for variant TSPYL5 SNP genotypes, but not in those homozygous for wild-type alleles. TSPYL5 knockdown decreased, and overexpression increased aromatase (CYP19A1) expression in MCF-7 cells, LCLs, and adipocytes through the skin/adipose (I.4) promoter. Chromatin immunoprecipitation assay showed that TSPYL5 bound to the CYP19A1 I.4 promoter. A putative TSPYL5 binding motif was identified in 43 genes, and TSPYL5 appeared to function as a transcription factor for most of those genes. In summary, genome-wide significant SNPs in TSPYL5 were associated with elevated plasma E2 in postmenopausal breast cancer patients. SNP rs2583506 created a functional estrogen response element, and LCLs with variant SNP genotypes displayed increased E2-dependent TSPYL5 expression. TSPYL5 induced CYP19A1 expression and that of many other genes. These studies have revealed a novel mechanism for regulating aromatase expression and plasma E2 concentrations in postmenopausal women with ER(+) breast cancer.

Genetic association with overall survival of taxane-treated lung cancer patients - a genome-wide association study in human lymphoblastoid cell lines followed by a clinical association study.

BACKGROUND: Taxane is one of the first line treatments of lung cancer. In order to identify novel single nucleotide polymorphisms (SNPs) that might contribute to taxane response, we performed a genome-wide association study (GWAS) for two taxanes, paclitaxel and docetaxel, using 276 lymphoblastoid cell lines (LCLs), followed by genotyping of top candidate SNPs in 874 lung cancer patient samples treated with paclitaxel. METHODS: GWAS was performed using 1.3 million SNPs and taxane cytotoxicity IC50 values for 276 LCLs. The association of selected SNPs with overall survival in 76 small or 798 non-small cell lung cancer (SCLC, NSCLC) patients were analyzed by Cox regression model, followed by integrated SNP-microRNA-expression association analysis in LCLs and siRNA screening of candidate genes in SCLC (H196) and NSCLC (A549) cell lines. RESULTS: 147 and 180 SNPs were associated with paclitaxel or docetaxel IC50s with p-values <10-4 in the LCLs, respectively. Genotyping of 153 candidate SNPs in 874 lung cancer patient samples identified 8 SNPs (p-value < 0.05) associated with either SCLC or NSCLC patient overall survival. Knockdown of PIP4K2A, CCT5, CMBL, EXO1, KMO and OPN3, genes within 200 kb up-/downstream of the 3 SNPs that were associated with SCLC overall survival (rs1778335, rs2662411 and rs7519667), significantly desensitized H196 to paclitaxel. SNPs rs2662411 and rs1778335 were associated with mRNA expression of CMBL or PIP4K2A through microRNA (miRNA) hsa-miR-584 or hsa-miR-1468. CONCLUSIONS: GWAS in an LCL model system, joined with clinical translational and functional studies, might help us identify genetic variations associated with overall survival of lung cancer patients treated paclitaxel.