CAP-RNAseq: an integrated pipeline for functional annotation and prioritization of co-expression clusters.

Cluster analysis is one of the most widely used exploratory methods for visualization and grouping of gene expression patterns across multiple samples or treatment groups. Although several existing online tools can annotate clusters with functional terms, there is no all-in-one webserver to effectively prioritize genes/clusters using gene essentiality as well as congruency of mRNA-protein expression. Hence, we developed CAP-RNAseq that makes possible (1) upload and clustering of bulk RNA-seq data followed by identification, annotation and network visualization of all or selected clusters; and (2) prioritization using DepMap gene essentiality and/or dependency scores as well as the degree of correlation between mRNA and protein levels of genes within an expression cluster. In addition, CAP-RNAseq has an integrated primer design tool for the prioritized genes. Herein, we showed using comparisons with the existing tools and multiple case studies that CAP-RNAseq can uniquely aid in the discovery of co-expression clusters enriched with essential genes and prioritization of novel biomarker genes that exhibit high correlations between their mRNA and protein expression levels. CAP-RNAseq is applicable to RNA-seq data from different contexts including cancer and available at http://konulabapps.bilkent.edu.tr:3838/CAPRNAseq/ and the docker image is downloadable from https://hub.docker.com/r/konulab/caprnaseq.

Long-Term Acetylcholinesterase Depletion Alters the Levels of Key Synaptic Proteins while Maintaining Neuronal Markers in the Aging Zebrafish (Danio rerio) Brain.

INTRODUCTION: Interventions targeting cholinergic neurotransmission like acetylcholinesterase (AChE) inhibition distinguish potential mechanisms to delay age-related impairments and attenuate deficits related to neurodegenerative diseases. However, the chronic effects of these interventions are not well described. METHODS: In the current study, global levels of cholinergic, cellular, synaptic, and inflammation-mediating proteins were assessed within the context of aging and chronic reduction of AChE activity. Long-term depletion of AChE activity was induced by using a mutant zebrafish line, and they were compared with the wildtype group at young and old ages. RESULTS: Results demonstrated that AChE activity was lower in both young and old mutants, and this decrease coincided with a reduction in ACh content. Additionally, an overall age-related reduction in AChE activity and the AChE/ACh ratio was observed, and this decline was more prominent in wildtype groups. The levels of an immature neuronal marker were upregulated in mutants, while a glial marker showed an overall reduction. Mutants had preserved levels of inhibitory and presynaptic elements with aging, whereas glutamate receptor subunit levels declined. CONCLUSION: Long-term AChE activity depletion induces synaptic and cellular alterations. These data provide further insights into molecular targets and adaptive responses following the long-term reduction of AChE activity that was also targeted pharmacologically to treat neurodegenerative diseases in human subjects.

Doxorubicin induces prolonged DNA damage signal in cells overexpressing DEK isoform-2.

DEK has a short isoform (DEK isoform-2; DEK2) that lacks amino acid residues between 49-82. The full-length DEK (DEK isoform-1; DEK1) is ubiquitously expressed and plays a role in different cellular processes but whether DEK2 is involved in these processes remains elusive. We stably overexpressed DEK2 in human bone marrow stromal cell line HS-27A, in which endogenous DEKs were intact or suppressed via short hairpin RNA (sh-RNA). We have found that contrary to ectopic DEK1, DEK2 locates in the nucleus and nucleolus, causes persistent γH2AX signal upon doxorubicin treatment, and couldn't functionally compensate for the loss of DEK1. In addition, DEK2 overexpressing cells were more sensitive to doxorubicin than DEK1-cells. Expressions of DEK1 and DEK2 in cell lines and primary tumors exhibit tissue specificity. DEK1 is upregulated in cancers of the colon, liver, and lung compared to normal tissues while both DEK1 and DEK2 are downregulated in subsets of kidney, prostate, and thyroid carcinomas. Interestingly, only DEK2 was downregulated in a subset of breast tumors suggesting that DEK2 can be modulated differently than DEK1 in specific cancers. In summary, our findings show distinct expression patterns and subcellular location and suggest non-overlapping functions between the two DEK isoforms.

Functional analysis of co-expression networks of zebrafish ace2 reveals enrichment of pathways associated with development and disease.

Human Angiotensin I Converting Enzyme 2 (ACE2) plays an essential role in blood pressure regulation and SARS-CoV-2 entry. ACE2 has a highly conserved, one-to-one ortholog (ace2) in zebrafish, which is an important model for human diseases. However, the zebrafish ace2 expression profile has not yet been studied during early development, between genders, across different genotypes, or in disease. Moreover, a network-based meta-analysis for the extraction of functionally enriched pathways associated with differential ace2 expression is lacking in the literature. Herein, we first identified significant development-, tissue-, genotype-, and gender-specific modulations in ace2 expression via meta-analysis of zebrafish Affymetrix transcriptomics datasets (ndatasets = 107); and the correlation analysis of ace2 meta-differential expression profile revealed distinct positively and negatively correlated local functionally enriched gene networks. Moreover, we demonstrated that ace2 expression was significantly modulated under different physiological and pathological conditions related to development, tissue, gender, diet, infection, and inflammation using additional RNA-seq datasets. Our findings implicate a novel translational role for zebrafish ace2 in organ differentiation and pathologies observed in the intestines and liver.

Global miRNA expression of bone marrow mesenchymal stem/stromal cells derived from Fanconi anemia patients.

Fanconi anemia (FA) is a rare genetic disorder characterized by genomic instability, developmental defects, and bone marrow (BM) failure. Hematopoietic stem cells (HSCs) in BM interact with the mesenchymal stem/stromal cells (MSCs); and this partly sustains the tissue homeostasis. MicroRNAs (miRNAs) can play a critical role during these interactions possibly via paracrine mechanisms. This is the first study addressing the miRNA profile of FA BM-MSCs obtained before and after BM transplantation (preBMT and postBMT, respectively). Non-coding RNA expression profiling and quality control analyses were performed in Donors (n = 13), FA preBMT (n = 11), and FA postBMT (n = 6) BM-MSCs using GeneChip miRNA 2.0 Array. Six Donor-FA preBMT pairs were used to identify a differentially expressed miRNA expression signature containing 50 miRNAs, which exhibited a strong correlation with the signature obtained from unpaired samples. Five miRNAs (hsa-miR-146a-5p, hsa-miR-148b-3p, hsa-miR-187-3p, hsa-miR-196b-5p, and hsa-miR-25-3p) significantly downregulated in both the paired and unpaired analyses were used to generate the BM-MSCs' miRNA-BM mononuclear mRNA networks upon integration of a public dataset (GSE16334; studying Donor versus FA samples). Functionally enriched KEGG pathways included cellular senescence, miRNAs, and pathways in cancer. Here, we showed that hsa-miR-146a-5p and hsa-miR-874-3p were rescued upon BMT (n = 3 triplets). The decrease in miR-146a-5p was also validated using RT-qPCR and emerged as a strong candidate as a modulator of BM mRNAs in FA patients.

Discovery of Cancer-Specific and Independent Prognostic Gene Subsets of the Slit-Robo Family Using TCGA-PANCAN Datasets.

The Slit-Robo family of axon guidance molecules works in concert, playing important roles in organ development and cancer. Expressions of individual Slit-Robo genes have been used in calculating univariable hazard ratios (HRuni) for predicting cancer prognosis in the literature. However, Slit-Robo members do not act independently; hence, hazard ratios from multivariable Cox regression (HRmulti) on the whole gene set can further lead to identification of cancer-specific, novel, and independent prognostic gene pairs or modules. Herein, we obtained mRNA expressions of the Slit-Robo family consisting of four Robos (ROBO1/2/3/4) and three Slits (SLIT1/2/3), along with four types of survival outcome across cancers found in the Cancer Genome Atlas (TCGA). We used cluster heat maps to visualize closely associated pairs/modules of prognostic genes across 33 different cancers. We found a smaller number of significant genes in HRmulti than in HRuni, suggesting that the former analysis was less redundant. High ROBO4 expression emerged as relatively protective within the family, in both types of HR analyses. Multivariable Cox regression, on the other hand, revealed significantly more HR signatures containing Slit-Robo pairs acting in opposing directions than those containing Slit-Slit or Robo-Robo pairs for disease-specific survival. Furthermore, we discovered, through the online app SmulTCan's lasso regression, Slit-Robo gene subsets that significantly differentiated between high- versus low-risk prognosis patient groups, particularly for renal cancers and low-grade glioma. The statistical pipeline reported herein can help test independent and significant pairs/modules within a codependent gene family for cancer prognostication, and thus should also prove useful in personalized/precision medicine research.

SmulTCan: A Shiny application for multivariable survival analysis of TCGA data with gene sets.

BACKGROUND: Survival analysis is widely used in cancer research, and although several methods exist in R, there is the need for a more interactive, flexible, yet comprehensive online tool to analyze gene sets using Cox proportional hazards (CPH) models. The web-based Shiny application (app) SmulTCan extends existing tools to multivariable CPH models of gene sets-as exemplified using the netrins and their receptors (netrins-receptors). It can be used to identify survival gene signatures (GSs) and select the best subsets of input gene, microRNA, methylation level, and copy number variation sets from the Cancer Genome Atlas (TCGA). OBJECTIVES: To create a tool for CPH model building and best subset selection, using survival data from TCGA with input gene expression files from UCSC Xena. Furthermore, we aim to analyze the input TSV file of netrins-receptors in SmulTCan and discuss our findings. METHODS: SmulTCan uses Shiny's reactivity with built-in R functions from packages for CPH model analysis and best subset selection including "survminer", "riskRegression", "rms", "glmnet", and "BeSS". RESULTS: Results from the SmulTCan app with the netrins-receptors gene set indicated unique hazard ratio GSs in certain renal and neural cancers, while the best subsets for this gene set, obtained via the app, could differentiate between prognostic outcomes in these cancers. AVAILABILITY: SmulTCan is available at http://konulabapps.bilkent.edu.tr:3838/SmulTCan/. The input file for netrins-receptors is available in the online version of this paper. TCGA dataset folders containing survival files are available through https://github.com/aozh7/SmulTCan/. SUPPLEMENTARY INFORMATION: The supplementary information (SI) accompanies the online version of this article.

CoVrimer: A tool for aligning SARS-CoV-2 primer sequences and selection of conserved/degenerate primers.

As mutations in SARS-CoV-2 virus accumulate rapidly, novel primers that amplify this virus sensitively and specifically are in demand. We have developed a webserver named CoVrimer by which users can search for and align existing or newly designed conserved/degenerate primer pair sequences against the viral genome and assess the mutation load of both primers and amplicons. CoVrimer uses mutation data obtained from an online platform established by NGDC-CNCB (12 May 2021) to identify genomic regions, either conserved or with low levels of mutations, from which potential primer pairs are designed and provided to the user for filtering based on generalized and SARS-CoV-2 specific parameters. Alignments of primers and probes can be visualized with respect to the reference genome, indicating variant details and the level of conservation. Consequently, CoVrimer is likely to help researchers with the challenges posed by viral evolution and is freely available at http://konulabapps.bilkent.edu.tr:3838/CoVrimer/.

Advanced preclinical models for evaluation of drug-induced liver injury - consensus statement by the European Drug-Induced Liver Injury Network [PRO-EURO-DILI-NET].

Drug-induced liver injury (DILI) is a major cause of acute liver failure (ALF) and one of the leading indications for liver transplantation in Western societies. Given the wide use of both prescribed and over the counter drugs, DILI has become a major health issue for which there is a pressing need to find novel and effective therapies. Although significant progress has been made in understanding the molecular mechanisms underlying DILI, our incomplete knowledge of its pathogenesis and inability to predict DILI is largely due to both discordance between human and animal DILI in preclinical drug development and a lack of models that faithfully recapitulate complex pathophysiological features of human DILI. This is exemplified by the hepatotoxicity of acetaminophen (APAP) overdose, a major cause of ALF because of its extensive worldwide use as an analgesic. Despite intensive efforts utilising current animal and in vitro models, the mechanisms involved in the hepatotoxicity of APAP are still not fully understood. In this expert Consensus Statement, which is endorsed by the European Drug-Induced Liver Injury Network, we aim to facilitate and outline clinically impactful discoveries by detailing the requirements for more realistic human-based systems to assess hepatotoxicity and guide future drug safety testing. We present novel insights and discuss major players in APAP pathophysiology, and describe emerging in vitro and in vivo pre-clinical models, as well as advanced imaging and in silico technologies, which may improve prediction of clinical outcomes of DILI.

Transcriptome profiles associated with selenium-deficiency-dependent oxidative stress identify potential diagnostic and therapeutic targets in liver cancer cells.

Hepatocellular carcinoma (HCC) is one of the most common cancer types with high mortality rates and displays increased resistance to various stress conditions such as oxidative stress. Conventional therapies have low efficacies due to resistance and off-target effects in HCC. Here we aimed to analyze oxidative stress-related gene expression profiles of HCC cells and identify genes that could be crucial for novel diagnostic and therapeutic strategies. To identify important genes that cause resistance to reactive oxygen species (ROS), a model of oxidative stress upon selenium (Se) deficiency was utilized. The results of transcriptome-wide gene expression data were analyzed in which the differentially expressed genes (DEGs) were identified between HCC cell lines that are either resistant or sensitive to Se-deficiency-dependent oxidative stress. These DEGs were further investigated for their importance in oxidative stress resistance by network analysis methods, and 27 genes were defined to have key roles; 16 of which were previously shown to have impact on liver cancer patient survival. These genes might have Se-deficiency-dependent roles in hepatocarcinogenesis and could be further exploited for their potentials as novel targets for diagnostic and therapeutic approaches.

Increased SGK1 activity potentiates mineralocorticoid/NaCl-induced kidney injury.

Serum and glucocorticoid-regulated kinase 1 (SGK1) stimulates aldosterone-dependent renal Na+ reabsorption and modulates blood pressure. In addition, genetic ablation or pharmacological inhibition of SGK1 limits the development of kidney inflammation and fibrosis in response to excess mineralocorticoid signaling. In this work, we tested the hypothesis that a systemic increase in SGK1 activity would potentiate mineralocorticoid/salt-induced hypertension and kidney injury. To that end, we used a transgenic mouse model with increased SGK1 activity. Mineralocorticoid/salt-induced hypertension and kidney damage was induced by unilateral nephrectomy and treatment with deoxycorticosterone acetate and NaCl in the drinking water for 6 wk. Our results show that although SGK1 activation did not induce significantly higher blood pressure, it produced a mild increase in glomerular filtration rate, increased albuminuria, and exacerbated glomerular hypertrophy and fibrosis. Transcriptomic analysis showed that extracellular matrix- and immune response-related terms were enriched in the downregulated and upregulated genes, respectively, in transgenic mice. In conclusion, we propose that systemically increased SGK1 activity is a risk factor for the development of mineralocorticoid-dependent kidney injury in the context of low renal mass and independently of blood pressure.NEW & NOTEWORTHY Increased activity of the protein kinase serum and glucocorticoid-regulated kinase 1 may be a risk factor for accelerated renal damage. Serum and glucocorticoid-regulated kinase 1 expression could be a marker for the rapid progression toward chronic kidney disease and a potential therapeutic target to slow down the process.

CHRNA5 belongs to the secondary estrogen signaling network exhibiting prognostic significance in breast cancer.

PURPOSE: Cholinergic signals can be important modulators of cellular signaling in cancer. We recently have shown that knockdown of nicotinic acetylcholine receptor subunit alpha 5, CHRNA5, diminishes the proliferative potential of breast cancer cells. However, modulation of CHRNA5 expression in the context of estrogen signaling and its prognostic implications in breast cancer remained unexplored. METHODS: Meta-analyses of large breast cancer microarray cohorts were used to evaluate the association of CHRNA5 expression with estrogen (E2) treatment, estrogen receptor (ER) status and patient prognosis. The results were validated through RT-qPCR analyses of multiple E2 treated cell lines, CHRNA5 depleted MCF7 cells and across a breast cancer patient cDNA panel. We also calculated a predicted secondary (PS) score representing direct/indirect induction of gene expression by E2 based on a public dataset (GSE8597). Co-expression analysis was performed using a weighted gene co-expression network analysis (WGCNA) pipeline. Multiple other publicly available datasets such as CCLE, COSMIC and TCGA were also analyzed. RESULTS: Herein we found that CHRNA5 expression was induced by E2 in a dose- and time-dependent manner in breast cancer cell lines. ER- breast tumors exhibited higher CHRNA5 expression levels than ER+ tumors. Independent meta-analysis for survival outcome revealed that higher CHRNA5 expression was associated with a worse prognosis in untreated breast cancer patients. Furthermore, CHRNA5 and its co-expressed gene network emerged as secondarily induced targets of E2 stimulation. These targets were largely downregulated by exposure to CHRNA5 siRNA in MCF7 cells while the response of primary ESR1 targets was dependent on the direction of the PS-score. Moreover, primary and secondary target genes were uncoupled and clustered distinctly based on multiple public datasets. CONCLUSION: Our findings strongly associate increased expression of CHRNA5 and its co-expression network with secondary E2 signaling and a worse prognosis in breast cancer.

Zebrafish Xenotransplantation Models for Studying Gene Function and Drug Treatment in Hepatocellular Carcinoma.

INTRODUCTION: Zebrafish is a promising model organism for human disease including hepatocellular cancer (HCC). Recently, zebrafish has emerged also as a host for xenograft studies of liver cancer cell lines and patient derived tumors of HCC. Zebrafish embryos enable drug screening and gene function studies of xenografted cells via ease of microinjection and visualization of tumor growth and metastasis. OBJECTIVES: In this review, we aimed to overview zebrafish HCC and liver cancer xenotransplantation studies focusing on 'gene functional analysis' and 'drug/chemical screening'. METHODS: Herein, a comprehensive literature search was performed for liver and HCC xenografts in zebrafish on PubMed using different key words and filters for molecular modifications or drug exposure. RESULTS: Our literature search revealed around 250 studies which were filtered and summarized in a table (Table 1) revealing comprehensive collection of experimental and technical details on microinjection, injected cell lines, molecular modifications of injected cells, types and doses of drug treatments as well as biological assessments. CONCLUSION: This review provides a platform for HCC and liver xenografts and highlights studies performed to understand gene functionality and drug efficacy in vivo in zebrafish.

ZenoFishDb v1.1: A Database for Xenotransplantation Studies in Zebrafish.

Rapidly accumulating literature has proven feasibility of the zebrafish xenograft models in cancer research. Nevertheless, online databases for searching the current zebrafish xenograft literature are in great demand. Herein, we have developed a manually curated database, called ZenoFishDb v1.1 (https://konulab.shinyapps.io/zenofishdb), based on R Shiny platform aiming to provide searchable information on ever increasing collection of zebrafish studies for cancer cell line transplantation and patient-derived xenografts (PDXs). ZenoFishDb v1.1 user interface contains four modules: DataTable, Visualization, PDX Details, and PDX Charts. The DataTable and Visualization pages represent xenograft study details, including injected cell lines, PDX injections, molecular modifications of cell lines, zebrafish strains, as well as technical aspects of the xenotransplantation procedures in table, bar, and/or pie chart formats. The PDX Details module provides comprehensive information on the patient details in table format and can be searched and visualized. Overall, ZenoFishDb v1.1 enables researchers to effectively search, list, and visualize different technical and biological attributes of zebrafish xenotransplantation studies particularly focusing on the new trends that make use of reporters, RNA interference, overexpression, or mutant gene constructs of transplanted cancer cells, stem cells, and PDXs, as well as distinguished host modifications.

Design, synthesis and anticancer/antiestrogenic activities of novel indole-benzimidazoles.

Indole-benzimidazoles have recently gained attention due to their antiproliferative and antiestrogenic effects. However, their structural similarities and molecular mechanisms shared with selective estrogen receptor modulators (SERMs) have not yet been investigated. In this study, we synthesized novel ethylsulfonyl indole-benzimidazole derivatives by substituting the first (R1) and fifth (R2) positions of benzimidazole and indole groups, respectively. Subsequently, we performed 1H NMR, 13C NMR, and Mass spectral and in silico docking analyses, and anticancer activity screening studies of these novel indole-benzimidazoles. The antiproliferative effects of indole-benzimidazoles were found to be more similar between the estrogen (E2) responsive cell lines MCF-7 and HEPG2 in comparison to the Estrogen Receptor negative (ER-) cell line MDA-MB-231. R1:p-fluorobenzyl group members were selected as lead compounds for their potent anticancer effects and moderate structural affinity to ER. Microarray expression profiling and gene enrichment analyses (GSEA) of the selected compounds (R1:p-fluorobenzyl: 48, 49, 50, 51; R1:3,4-difluorobenzyl: 53) helped determine the similarly modulated cellular signaling pathways among derivatives. Moreover, we identified known compounds that have significantly similar gene signatures to that of 51 via queries performed in LINCS database; and further transcriptomics comparisons were made using public GEO datasets (GSE35428, GSE7765, GSE62673). Our results strongly demonstrate that these novel indole-benzimidazoles can modulate ER target gene expression as well as dioxin-mediated aryl hydrocarbon receptor and amino acid deprivation-mediated integrated stress response signaling in a dose-dependent manner.

The Transcription Factor Elf3 Is Essential for a Successful Mesenchymal to Epithelial Transition.

The epithelial to mesenchymal transition (EMT) and the mesenchymal to epithelial transition (MET) are two critical biological processes that are involved in both physiological events such as embryogenesis and development and also pathological events such as tumorigenesis. They present with dramatic changes in cellular morphology and gene expression exhibiting acute changes in E-cadherin expression. Despite the comprehensive understanding of EMT, the regulation of MET is far from being understood. To find novel regulators of MET, we hypothesized that such factors would correlate with Cdh1 expression. Bioinformatics examination of several expression profiles suggested Elf3 as a strong candidate. Depletion of Elf3 at the onset of MET severely impaired the progression to the epithelial state. This MET defect was explained, in part, by the absence of E-cadherin at the plasma membrane. Moreover, during MET, ELF3 interacts with the Grhl3 promoter and activates its expression. Our findings present novel insights into the regulation of MET and reveal ELF3 as an indispensable guardian of the epithelial state. A better understanding of MET will, eventually, lead to better management of metastatic cancers.

Lower connectivity of tumor coexpression networks is not specific to cancer.

Global level network analysis of molecular links is necessary for systems level view of complex diseases like cancer. Using genome-wide expression datasets, we constructed and compared gene co-expression based specific networks of pre-cancerous tumors (adenoma) and cancerous tumors (carcinoma) with paired normal networks to assess for any possible changes in network connectivity. Previously, loss of connectivity was reported as a characteristics of cancer samples. Here, we observed that pre-cancerous conditions also had significantly less connections than paired normal samples. We observed a loss of connectivity trend for colorectal adenoma, aldosterone producing adenoma and uterine leiomyoma. We also showed that the loss of connectivity trend is not specific to positive or negative correlation based networks. Differential hub genes, which were the most highly differentially less connected genes in tumor, were mostly different between different datasets. No common gene list could be defined which underlies the lower connectivity of tumor specific networks. Connectivity of colorectal cancer methylation targets was different from other genes. Extracellular space related terms were enriched in negative correlation based differential hubs and common methylation targets of colorectal carcinoma. Our results indicate a systems level change of lower connectivity as cells transform to not only cancer but also pre-cancerous conditions. This systems level behavior could not be attributed to a group of genes.

PKNOX2 expression and regulation in the bone marrow mesenchymal stem cells of Fanconi anemia patients and healthy donors.

HOX and TALE transcription factors are important regulators of development and homeostasis in determining cellular identity. Deregulation of this process may drive cancer progression. The aim of this study was to investigate the expression of these transcription factors in the bone marrow derived mesenchymal stem cells (BM-MSCs) of Fanconi anemia (FA) patients, which is a cancer-predisposing disease. Expression levels of HOX and TALE genes in BM-MSCs were obtained from FA patients and healthy donors by RT-qPCR and highly conserved expression levels were observed between patient and donor cells, except PKNOX2, which is a member of TALE class. PKNOX2 was significantly downregulated in FA cells compared to donors (P < 0.05). PKNOX2 expression levels did not change with diepoxybutane (DEB), a DNA crosslinking agent, in either donor or FA cells except one patient's with a truncation mutation of FANCA. A difference of PKNOX2 protein level was not obtained between FA patient and donor BM-MSCs by western blot analysis. When human TGF-ß1 (rTGF-ß1) recombinant protein was provided to the cultures, PKNOX2 as well as TGF-ß1 expression increased both in FA and donor BM-MSCs in a dose dependent manner. 5 ng/mL rTGF-ß stimulation had more dominant effect on the gene expression of donor BM-MSCs compared to FA cells. Decreased PKNOX2 expression in FA BM-MSCs may provide new insights into the molecular pathophysiology of the disease and TGF-ß1 levels of the microenvironment may be the cause of PKNOX2 downregulation.

Cholinergic Receptor Nicotinic Alpha 5 (CHRNA5) RNAi is associated with cell cycle inhibition, apoptosis, DNA damage response and drug sensitivity in breast cancer.

Cholinergic Receptor Nicotinic Alpha 5 (CHRNA5) is an important susceptibility locus for nicotine addiction and lung cancer. Depletion of CHRNA5 has been associated with reduced cell viability, increased apoptosis and alterations in cellular motility in different cancers yet not in breast cancer. Herein we first showed the expression of CHRNA5 was variable and positively correlated with the fraction of total genomic alterations in breast cancer cell lines and tumors indicating its potential role in DNA damage response (DDR). Next, we demonstrated that silencing of CHRNA5 expression in MCF7 breast cancer cell line by RNAi affected expression of genes involved in cytoskeleton, TP53 signaling, DNA synthesis and repair, cell cycle, and apoptosis. The transcription profile of CHRNA5 depleted MCF7 cells showed a significant positive correlation with that of A549 lung cancer cell line while exhibiting a negative association with the CHRNA5 co-expression profile obtained from Cancer Cell Line Encylopedia (CCLE). Moreover, it exhibited high similarities with published MCF7 expression profiles obtained from exposure to TP53 inducer nutlin-3a and topoisomerase inhibitors. We then demonstrated that CHRNA5 siRNA treatment reduced cell viability and DNA synthesis indicating G1 arrest while it significantly increased apoptotic sub-G1 cell population. Accordingly, we observed lower levels of phosphorylated RB (Ser807/811) and an increased BAX/BCL2 ratio in RNAi treated MCF7 cells. We also showed that CHRNA5 RNAi transcriptome correlated negatively with DDR relevant gene expression profile in breast cancer gene expression datasets while the coexposure to topoisomerase inhibitors in the presence of CHRNA5 RNAi enhanced chemosensitivity potentially due to reduced DDR. CHRNA5 RNAi consistently lowered total CHEK1 mRNA and protein levels as well as phosphorylated CHEK1 (Ser345) in MCF7 cells. We also detected a significant positive correlation between the expression levels of CHRNA5 and CHEK1 in CCLE, TCGA and METABRIC breast cancer datasets. Our study suggests CHRNA5 RNAi is associated with cell cycle inhibition, apoptosis as well as reduced DDR and increased drug sensitivity in breast cancer yet future studies are warranted since dose- and cell line-specific differences exist in response to CHRNA5 depletion. Gene expression microarray data can be accessed from GEO database under the accession number GSE89333.