- Research note
- Open Access
Associations of two-pore domain potassium channels and triple negative breast cancer subtype in The Cancer Genome Atlas: systematic evaluation of gene expression and methylation
BMC Research Notes volume 10, Article number: 475 (2017)
It is unclear whether 2-pore domain potassium channels are novel molecular markers with differential expression related to biologically aggressive triple-negative type breast tumors. Our objective was to systematically evaluate associations of 2-pore domain potassium channel gene expression and DNA methylation with triple-negative subtype in The Cancer Genome Atlas invasive breast cancer dataset. Methylation and expression data for all fifteen 2-pore domain potassium family genes were examined for 1040 women, and associations with triple-negative subtype (vs. luminal A) were evaluated using age/race adjusted generalized-linear models, with Bonferroni-corrected significance thresholds. Subtype associated CpG loci were evaluated for functionality related to expression using Spearman’s correlation.
Overexpression of KCNK5, KCNK9 and KCNK12, and underexpression of KCNK6 and KCNK15, were significantly associated with triple-negative subtype (Bonferroni-corrected p < 0.0033). A total of 195 (114 hypomethylated and 81 hypermethylated) CpG loci were found to be significantly associated with triple-negative subtype (Bonferroni-corrected p < 8.22 × 10−8). Significantly negatively correlated expression patterns that were differentially observed in triple-negative vs. luminal A subtype were demonstrated for: KCNK2 (gene body: cg04923840, cg13916421), KCNK5 (gene body: cg05255811, cg18705155, cg09130674, cg21388745, cg00859574) and KCNK9 (TSS1500: cg21415530, cg12175729; KCNK9/TRAPPC9 intergenic region: cg17336929, cg25900813, cg03919980). CpG loci listed for KCNK5 and KCNK9 all showed relative hypomethylation for probability of triple-negative vs. luminal A subtype. Triple-negative subtype was associated with distinct 2-pore domain potassium channel expression patterns. Both KCNK5 and KCNK9 overexpression appeared to be functionally related to CpG loci hypomethylation.
Two-pore domain potassium (K2p) channels enable background leak of potassium (K+) ions and the K2p-family has 15 members [1, 2]. K2p-channels are important for baseline cellular activity at rest including membrane-potential, calcium homeostasis and cell-volume regulation . Evidence from laboratory studies supports the hypothesis that alterations in expression/function of K2p-channels may play a role in cancer development/progression [3,4,5]. The role of K2p-channels in breast cancer (BC) is emerging and recent reviews suggest potential clinical utility [4, 6, 7]. Williams et al. examined K2p-channel expression in a microarray database study, and reported that all but five members showed altered expression in BC .
Preliminary reports suggest that some K2p-channels may be novel molecular markers with differential expression related to biologically-aggressive triple-negative (TN) type BC. In estrogen-receptor-alpha positive (ER+) cell lines, KCNK5 expression appears to be under regulatory control of ER-alpha, and estrogen-response-elements (ERE) are found in the enhancer region of KCNK5 . Clarke et al. suggest that upregulation of KCNK5 is associated with poor outcome for TNBC-related basal-like subtype . Conway et al. suggest that KCNK4 may be differentially methylated according to race, with relatively higher median DNA-methylation observed for non-Hispanic (nH) black vs. nH-white BC . Other studies similarly suggest that differential methylation of KCNK9 may be related to nH-black race and TNBC [11,12,13,14]. Mu et al. first described the KCNK9 gene on 8q24.3 as a potential proto-oncogene; gene-amplification (10%) and protein-overexpression (44%) were detected in BC, but not in normal tissue-controls . DNA copy-number loss has been observed for KCNK12 in a small clinical study of tubular BC .
K2p-channels also hold promise for innovative oncologic therapeutic/prevention strategies [7, 17, 18]. Current molecular-epidemiologic data characterizing K2p-channels in clinical BC are limited. Our study goal was to systematically evaluate associations between K2p-gene-family mRNA-expression and DNA-methylation with TNBC in The-Cancer-Genome-Atlas (TCGA), which contains genomic information from DNA-methylation and mRNA-expression arrays for almost 1100 human BC.
Materials and methods
Sample and procedure
Publicly available BC data was downloaded/aggregated (KAD) from TCGA-Data-Matrix and cBioPortal-for-Cancer-Genomics (provisional dataset) web-sites (February 2016) [19, 20]. TCGA molecular BC dataset was generated by TCGA-Research-Network and details have been previously published [21, 22]. The assembled dataset of invasive BC (using TCGA sample-code suffix −01 for primary solid-tumors) included information on: (1) clinical-pathological factors: age/stage at diagnosis, menopausal status, lymph-node status, estrogen/progesterone receptor (ER/PR) status, human-epidermal-growth-factor-receptor-2 (HER2) status, and race/ethnicity; and (2) K2p-gene-markers including: (a) Level 3 DNA-methylation probes (Illumina-Infinium-Human-Methylation-450K-BeadChip beta-values (450K array), n = 767) ; and (b) Level 3 RNA-Seq gene-expression data (Illumina-HiSeq-platform, gene level RNA-Seq by Expectation Maximization (RSEM)-normalized and log2-transformed values, n = 959) . Cytosine–guanine (GpG) methylation loci within 25 kb from either end of the genes of interest were included for examination (UCSC-Genome-Browser , hg37); 724 GpG loci were included, but was reduced to 608 after exclusion of probes with null-reads. Analytic dataset consisted of 1040 women with information on stage, and all women with available data were included in analyses. The study was approved by UIC IRB.
Breast cancer subtype
The primary study-endpoint was a binary variable of combined TN-status, developed by aggregating data from immunohistochemical (IHC) and intrinsic subtype assays (PAM50) with a referent of luminal-A status (ER/PR+/HER2− status). For sensitivity analyses, data on Intrinsic Subtype classification was obtained from supplementary publication data based on PAM50 assay , and data on Integrative Cluster (IC) classification based on copy-number-change was obtained from Cancer-Research-UK-Cambridge-Institute (Oscar Rueda and Carlos Caldas) [26,27,28]. Secondary endpoints in sensitivity analyses included: ER/PR− status (IHC; referent ER/PR+); basal-subtype (PAM50; ER/PR−/HER2−, and cytokeratin 5/6+ and/or HER1+; referent luminal-A); and IC10-type (referent IC3). Gene-expression patterns were also compared for RSEM and Agilent arrays.
For categorical variables, Chi square test-of-association was used to examine significance of relationships. For continuous variables, t-test was used to compare means as appropriate to distribution characteristics; and for non-normally distributed data, nonparametric Wilcoxon rank-sum test was used.
Tertiles were used for categorization of gene-expression, except where limited by small sample-signals, then binary mRNA-overexpression was defined as a z-score > +1 standard deviation (SD) . Association of all K2p gene-expression and DNA-methylation with TN-subtype was evaluated using age (years) and race (nH-white, nH-black) adjusted generalized-linear regression-models (glm) with binomial-distribution and log-link to estimate prevalence risk-ratios (RRs) and their corresponding 95% confidence-intervals (CIs). The significance threshold for gene-expression and DNA-methylation was p < 0.05 and Bonferroni-corrected p-values were used where appropriate (Bonferroni-adjusted threshold for expression was p < 0.0033; and for methylation p < 8.22 × 10−8).
Methylation glm results were sorted on p-values (smallest to largest) and all loci associated with TN-subtype were selected for reporting and further analysis; loci were then evaluated for functionality related to expression using Spearman’s correlation. Methylation was deemed significantly correlated with expression if the Spearman’s Rho-value was either ≤−0.2 or ≥+0.2 and associated with a p-value <0.05 .
Characteristics of the study sample are presented in Table 1. Overall, most of the sample was age >50 years (73%), postmenopausal (76%), nH-white (86%), stage-II (58%), and ER/PR+ (66%) and HER2− (89%); 631 patients had data on TN (n = 159) vs. luminal-A (n = 472) status.
All K2p gene-expression was non-normally distributed (Additional file 1: Table S1). Table 2 demonstrates that overexpression of KCNK5/KCNK9/KCNK12, and underexpression of KCNK6/KCNK15, were associated with TN-subtype in age/race adjusted models (all p < 0.0033). In fully-adjusted models, KCNK5/KCNK9/KCNK12/KCNK15 remained associated with TN-subtype.
Gene-expression sensitivity analyses (Additional file 2: Table S2) show similar K2p RSEM-expression patterns for ER/PR-, basal and IC10 type tumors; however, overexpression of KCNK7 was also seen with basal-subtype, and expression of KCNK9 and KCNK5 was not associated with basal and IC10 types, respectively. Comparative analyses using the Agilent expression-array (n = 590) (Additional file 2: Table S2) shows similar expression patterns for KCNK5/KCNK6/KCNK12/KCNK15, but additionally shows consistent underexpression of KCNK4 across subtype, and that KCNK9 was not associated with any subtype for this array.
Of 608 CpG loci, 195 (114 hypomethylated and 81 hypermethylated) loci were found to be associated with TN-subtype (p < 8.22 × 10−8, Additional file 3: Table S3). Significantly-associated negatively-correlated expression patterns differentially observed in TN vs. luminal-A subtype (Table 3; and Additional file 4: Table S4) were demonstrated for: KCNK2 (gene-body: cg04923840, cg13916421), KCNK5 (gene-body: cg05255811, cg18705155, cg09130674, cg21388745, cg00859574) and KCNK9 (TSS1500: cg21415530, cg12175729; KCNK9/TRAPPC9 intergenic: cg17336929, cg25900813, cg03919980). Loci listed for KCNK5/KCNK9 all showed relative hypomethylation for probability of TN vs. luminal-A subtype.
K2p RSEM-expression analyses show that specific genes are associated with TN-subtype (overexpression of KCNK5/KCNK9/KCNK12, and underexpression of KCNK6/KCNK15). Sensitivity analyses demonstrate consistency; similar RSEM-expression patterns of association were largely evident for ER/PR-, basal and IC10 related subtypes, although KCNK9 and KCNK5 overexpression failed to meet significance criteria for basal and IC10 types, respectively. Comparative analyses with the Agilent-array supports RSEM findings for KCNK5/KCNK6/KCNK12/KCNK15.
Regarding specificity for association of K2p-expression with tumor vs. normal sample-type, we used the MEXPRESS web tool for visualizing expression, DNA methylation and clinical TCGA data , and found that KCNK9 (p = 2.58 × 10−9) and KCNK12 (p = 4.63 × 10−10) overexpression appeared to be associated with tumor type, while KCNK5 overexpression appeared to have marginal association (p = 0.0757). KCNK6 and KCNK15 overexpression (associated in our study with luminal-A subtype) also appeared to be associated with tumor sample-type (both p < 2.2 × 10−16).
Specific methylation loci on KCNK2/KCNK5/KCNK9 were found to be significantly-associated with negatively-correlated gene-expression and differentially observed in TN vs. luminal-A subtype, which suggests potential functional significance. All select loci on KCNK5/KCNK9 had negative delta-beta values indicating relative hypomethylation in TN-subtype. The overall levels of negative expression Rho-values for KCNK5 and KCNK9 were more than twice larger for TN than luminal-A subtype. Hence it is plausible that KCNK5/KCNK9 overexpression related to TN-subtype may be functionally related to specific tumor CpG loci hypomethylation.
Prior literature suggests KCNK5 may be under control of ER-alpha signaling and upregulation is associated with worse prognosis for basal-type tumors [8, 9]. Although we found association of KCNK5/KCNK9 overexpression and TN-subtype, we were unable to relate to prognosis. Further, 17-beta-estradiol is known to induce KCNK5 gene-expression via EREs found in the enhancer region of the KCNK5 gene . Hence the observed pattern in our study could be interpreted as evidence of loss of regulatory control of ER-alpha signaling pathways. An alternate explanation could be that KCNK5 gene-expression in TNBC is incompletely understood. Prior studies have implicated KCNK9 as a proto-oncogene in BC and suggest importance in TNBC [12,13,14,15, 31]. KCNK9 is a maternally-imprinted gene with monoallelic-expression predominantly in brain tissue, but expression has been observed in breast tissue, and both are of ectodermal origin . It is theorized that KCNK9 overexpression may occur due to acquired relative hypomethylation and subsequent functionally biallelic-expression which may be equivalent to duplication of an active allele . Our findings regarding the functional relationship of KCNK9 methylation/expression may be consistent with this theory, as we observed increased KCNK9-expression together with relative hypomethylation at functional loci. We also observed a similar relationship for KCNK5. There is scant prior data on KCNK12 expression in BC and Williams et al. suggested no alteration in their study ; we found KCNK12 overexpression to be associated with TN subtype. Williams et al. described overexpression of KCNK6/KCNK15; we found that these were specifically overexpressed in luminal-A subtype.
Our findings are consistent with prior reports suggesting that some K2p-channels may be novel molecular-markers with differential expression related to biologically-aggressive TNBC, and may hold promise for innovative oncologic therapeutic/prevention strategies, as experience gained with pharmacological manipulation of K+ channels in other pathologies, might facilitate their use as molecular-targets in precision-medicine . Wallace et al. suggest that K+ channel activity demonstrates promise as a pharmacologic target in BC and may represent a mechanism through which phytoestrogens act (e.g. Genistein) . A recent study showed that inhibition of the intermediate-conductance calcium-activated K+ channel KCNN3 with specific blockers including the antifungal clotrimazole, suppressed cell proliferation, migration and epithelial–mesenchymal transition in TNBC cells, and illustrates the potential of K+ channels as novel targets . Sun et al. demonstrated that a monoclonal antibody (Y4) against KCNK9 extracellular-domain effectively inhibits growth of human lung cancer xenografts and BC metastasis in mice, and suggest that antibody-based KCNK9 targeting is a promising therapeutic strategy in KCNK9 overexpressing malignancies . Skaar et al. have been researching whether loss of KCNK9 imprint-control and/or differential methylation could be adapted in novel approaches for clinical diagnosis and targeted therapy for TNBC [11,12,13,14].
In summary, our results suggest that several K2p-channel methylation/expression markers are related to TNBC. KCNK5 and KCNK9 show overexpression and relative hypomethylation of specific CpG loci, and have negative methylation/gene-expression correlations, which are features that support functional biology. Study strengths include: TCGA is a large, well characterized BC resource; we used the target tumor-tissue of interest, not peripheral blood; and our findings are consistent with, and build on, other reports from TCGA. Our findings are considered preliminary but may lead to further exploration of K2p-channels as potential molecular-targets in precision-medicine, as pharmacological manipulation of K+ channels is currently feasible.
TCGA is not a population-based study or clinical trial and may be subject to selection bias. There was missing data on CpG loci with null-reads which may have been excluded due to quality-control protocols. Missing data may cause loss of analytic precision and power and could potentially bias results . Our analysis used specific CpG probes present on the 450K array, but a more comprehensive fine-mapping sequencing approach could facilitate interrogation of additional loci. Potential modeling inadequacies may also result from failure to include molecular-markers not available to us, such as protein-expression, non-coding RNAs and expression-quantitative-trait-loci. K2p copy-number-change may also be important and this analysis is planned as a future study. Regarding survival analysis, specific treatment data is unavailable, and since follow-up time and number of survival events are limited, we were unable to examine prognosis, as our study is largely under-powered to detect differences between groups. Regarding sensitivity analyses, overlap between the RSEM and Agilent expression-arrays was only 499 cases, and as such sampling differences could account for observed variation in expression profiling. These limitations make the present study results preliminary in nature and our findings require validation in additional studies.
2-pore domain potassium
The Cancer Genome Atlas
human epidermal growth factor receptor-2
University of Illinois at Chicago
Institutional Review Board
Patel AJ, Lazdunski M. The 2P-domain K+ channels: role in apoptosis and tumorigenesis. Pflugers Arch. 2004;448(3):261–73.
Enyedi P, Czirjak G. Molecular background of leak K+ currents: two-pore domain potassium channels. Physiol Rev. 2010;90(2):559–605.
Williams S, Bateman A, O’Kelly I. Altered expression of two-pore domain potassium (K2P) channels in cancer. PLoS ONE. 2013;8(10):e74589.
Comes N, Serrano-Albarras A, Capera J, Serrano-Novillo C, Condom E, Ramon YCS, Ferreres JC, Felipe A. Involvement of potassium channels in the progression of cancer to a more malignant phenotype. Biochem Biophys Acta. 2015;1848((10 Pt B)):2477–92.
Kondratskyi A, Kondratska K, Skryma R, Prevarskaya N. Ion channels in the regulation of apoptosis. Biochem Biophys Acta. 2015;1848(10 Pt B):2532–46.
Lastraioli E, Iorio J, Arcangeli A. Ion channel expression as promising cancer biomarker. Biochem Biophys Acta. 2015;1848(10 Pt B):2685–702.
Pardo LA, Stuhmer W. The roles of K(+) channels in cancer. Nat Rev Cancer. 2014;14(1):39–48.
Alvarez-Baron CP, Jonsson P, Thomas C, Dryer SE, Williams C. The two-pore domain potassium channel KCNK5: induction by estrogen receptor alpha and role in proliferation of breast cancer cells. Mol Endocrinol (Baltimore, Md). 2011;25(8):1326–36.
Clarke C, Madden SF, Doolan P, Aherne ST, Joyce H, O’Driscoll L, Gallagher WM, Hennessy BT, Moriarty M, Crown J, et al. Correlating transcriptional networks to breast cancer survival: a large-scale coexpression analysis. Carcinogenesis. 2013;34(10):2300–8.
Conway K, Edmiston SN, Tse CK, Bryant C, Kuan PF, Hair BY, Parrish EA, May R, Swift-Scanlan T. Racial variation in breast tumor promoter methylation in the carolina breast cancer study. Cancer Epidemiol Biomark Prev. 2015;24(6):921–30.
Luedi PP, Dietrich FS, Weidman JR, Bosko JM, Jirtle RL, Hartemink AJ. Computational and experimental identification of novel human imprinted genes. Genome Res. 2007;17(12):1723–30.
Skaar D, Gould M, Jirtle R, Seewaldt V. Altered imprinted gene DMR regulatory methylation in breast cancer. Cancer Epidemiol Biomark Prev. 2014;23(11 Supplement):02–3 (Abstract CN02-03).
Desai S, Skaar D, Dietze E, Ambrose A, Jirtle R, Seewaldt V. Abstract C65: loss of imprinting of KCNK9 in African American women with triple-negative breast cancer. Cancer Epidemiol Biomark Prev. 2014;23(11 Supplement):C65.
Skaar DA, Seewaldt VL, Gould MN, Jirtle RL. Abstract 4043: loss of imprinting of KCNK9 in breast cancer. Can Res. 2012;72(8 Supplement):4043.
Mu D, Chen L, Zhang X, See LH, Koch CM, Yen C, Tong JJ, Spiegel L, Nguyen KC, Servoss A, et al. Genomic amplification and oncogenic properties of the KCNK9 potassium channel gene. Cancer Cell. 2003;3(3):297–302.
Riener MO, Nikolopoulos E, Herr A, Wild PJ, Hausmann M, Wiech T, Orlowska-Volk M, Lassmann S, Walch A, Werner M. Microarray comparative genomic hybridization analysis of tubular breast carcinoma shows recurrent loss of the CDH13 locus on 16q. Hum Pathol. 2008;39(11):1621–9.
Zhang P, Yang X, Yin Q, Yi J, Shen W, Zhao L, Zhu Z, Liu J. Inhibition of SK4 potassium channels suppresses cell proliferation, migration and the epithelial–mesenchymal transition in triple-negative breast cancer cells. PLoS ONE. 2016;11(4):e0154471.
Wallace JL, Gow IF, Warnock M. The life and death of breast cancer cells: proposing a role for the effects of phytoestrogens on potassium channels. J Membr Biol. 2011;242(2):53–67.
Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, et al. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2(5):401–4.
Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, et al. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6(269):p11.
Cancer-Genome-Atlas-Network. Comprehensive molecular portraits of human breast tumours. Nature. 2012;490(7418):61–70.
Tomczak K, Czerwinska P, Wiznerowicz M. The cancer genome atlas (TCGA): an immeasurable source of knowledge. Contemp oncol (Poznan, Poland). 2015;19(1a):A68–77.
Illumina. Infinium HumanMethylation450K BeadChip. 2017; Array Support. https://support.illumina.com/array/array_kits/infinium_humanmethylation450_beadchip_kit.html. Accessed 09 Aug 2017.
Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011;12(1):323.
Kent WJ, Sugnet CW, Furey TS, Roskin KM, Pringle TH, Zahler AM, Haussler D. The human genome browser at UCSC. Genome Res. 2002;12(6):996–1006.
Ali H, Rueda OM, Chin SF, Curtis C, Dunning MJ, Aparicio S, Caldas C. Genome-driven integrated classification of breast cancer validated in over 7500 samples. Genome Biol. 2014;15(8):431.
Curtis C, Shah SP, Chin SF, Turashvili G, Rueda OM, Dunning MJ, Speed D, Lynch AG, Samarajiwa S, Yuan Y, et al. The genomic and transcriptomic architecture of 2000 breast tumours reveals novel subgroups. Nature. 2012;486(7403):346–52.
Dawson SJ, Rueda OM, Aparicio S, Caldas C. A new genome-driven integrated classification of breast cancer and its implications. EMBO J. 2013;32(5):617–28.
Tell RW, Horvath CM. Bioinformatic analysis reveals a pattern of STAT3-associated gene expression specific to basal-like breast cancers in human tumors. Proc Natl Acad Sci. 2014;111(35):12787–92.
Koch A, De Meyer T, Jeschke J, Van Criekinge W. MEXPRESS: visualizing expression, DNA methylation and clinical TCGA data. BMC Genom. 2015;16:636.
Pei L, Wiser O, Slavin A, Mu D, Powers S, Jan LY, Hoey T. Oncogenic potential of TASK3 (Kcnk9) depends on K+ channel function. Proc Natl Acad Sci USA. 2003;100(13):7803–7.
Sun H, Luo L, Lal B, Ma X, Chen L, Hann CL, Fulton AM, Leahy DJ, Laterra J, Li M. A monoclonal antibody against KCNK9K(+) channel extracellular domain inhibits tumour growth and metastasis. Nat commun. 2016;7:10339.
Sterne JA, White IR, Carlin JB, Spratt M, Royston P, Kenward MG, Wood AM, Carpenter JR. Multiple imputation for missing data in epidemiological and clinical research: potential and pitfalls. BMJ. 2009;338:b2393.
KAD drafted the manuscript and performed the statistical analysis. KAD, WZ, LS and MA contributed to the conception/design of the work and the analysis/interpretation of study data. All authors read and approved the final manuscript.
KAD is Assistant Professor, Epidemiology at University of Wisconsin-Milwaukee Joseph J. Zilber School of Public Health. In part, his research focuses on improving understanding of the molecular underpinnings of biologically aggressive breast tumor types and related molecular disparities, that are more often seen in non-Hispanic black women with breast cancer.
The results published here are in whole based upon data generated by the TCGA Research Network: (http://cancergenome.nih.gov/). The authors also recognize the valuable contributions of the participants, specimen donors and research groups who developed the TCGA breast cancer dataset resource.
The authors would also like to recognize Oscar Rueda and Carlos Caldas from University of Cambridge, Cancer Research UK Cambridge Institute, who shared data on Integrative Cluster classification.
Prior scientific presentations and abstracts:
The ASCO Annual Meeting 2015 (in part).
The AACR Science of Cancer Health Disparities in Racial/Ethnic Minorities and the Medically Underserved 2015 and 2016 (in part).
The authors declare that they have no competing interests.
Availability of data and materials
The original datasets analyzed during the current study are publicly available from the TCGA repository at: https://cghub.ucsc.edu/. Specific datasets used during the current study are also available from the corresponding author on reasonable request.
Consent for publication
Ethics approval and consent to participate
The study was approved by the University of Illinois at Chicago institutional review board (Study Number: 20150198-88043-1). Since the study analyzes publicly available non-identifiable data from TCGA and is technically ‘In Silico’, UIC IRB determined that this research did NOT meet the definition of human subject research. As such informed consent was not required for this study.
The Komen Foundation: Grant # KG111385.
The Cancer Foundation for Minority and Underserved Populations.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Dookeran, K.A., Zhang, W., Stayner, L. et al. Associations of two-pore domain potassium channels and triple negative breast cancer subtype in The Cancer Genome Atlas: systematic evaluation of gene expression and methylation. BMC Res Notes 10, 475 (2017). https://doi.org/10.1186/s13104-017-2777-4
- Potassium channels
- Breast cancer