In search of autophagy biomarkers in breast cancer: Receptor status and drug agnostic transcriptional changes during autophagy flux in cell lines

Autophagy drives drug resistance and drug-induced cancer cell cytotoxicity. Targeting the autophagy process could greatly improve chemotherapy outcomes. The discovery of specific inhibitors or activators has been hindered by challenges with reliably measuring autophagy levels in a clinical setting. We investigated drug-induced autophagy in breast cancer cell lines with differing ER/PR/Her2 receptor status by exposing them to known but divergent autophagy inducers each with a unique molecular target, tamoxifen, trastuzumab, bortezomib or rapamycin. Differential gene expression analysis from total RNA extracted during the earliest sign of autophagy flux showed both cell- and drug-specific changes. We analyzed the list of differentially expressed genes to find a common, cell- and drug-agnostic autophagy signature. Twelve mRNAs were significantly modulated by all the drugs and 11 were orthogonally verified with Q-RT-PCR (Klhl24, Hbp1, Crebrf, Ypel2, Fbxo32, Gdf15, Cdc25a, Ddit4, Psat1, Cd22, Ypel3). The drug agnostic mRNA signature was similarly induced by a mitochondrially targeted agent, MitoQ. In-silico analysis on the KM-plotter cancer database showed that the levels of these mRNAs are detectable in human samples and associated with breast cancer prognosis outcomes of Relapse-Free Survival in all patients (RSF), Overall Survival in all patients (OS), and Relapse-Free Survival in ER⁺ Patients (RSF ER⁺). High levels of Klhl24, Hbp1, Crebrf, Ypel2, CD22 and Ypel3 were correlated with better outcomes, whereas lower levels of Gdf15, Cdc25a, Ddit4 and Psat1 were associated with better prognosis in breast cancer patients. This gene signature uncovers candidate autophagy biomarkers that could be tested during preclinical and clinical studies to monitor the autophagy process.

Anthracycline-Induced Cardiotoxicity: Molecular Insights Obtained from Human-Induced Pluripotent Stem Cell-Derived Cardiomyocytes (hiPSC-CMs)

Anthracyclines are a class of chemotherapy drugs that are highly effective for the treatment of human cancers, but their clinical use is limited by associated dose-dependent cardiotoxicity. The precise mechanisms by which individual anthracycline induces cardiotoxicity are not fully understood. Human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) are emerging as a physiologically relevant model to assess drugs cardiotoxicity. Here, we describe an assay platform by coupling hiPSC-CMs and impedance measurement, which allows real-time monitoring of cardiomyocyte cellular index, beating amplitude, and beating rate. Using this approach, we have performed comparative studies on a panel of four anthracycline drugs (doxorubicin, epirubicin, idarubicin, and daunorubicin) which share a high degree of structural similarity but are associated with distinct cardiotoxicity profiles and maximum cumulative dose limits. Notably, results from our hiPSC-CMs impedance model (dose-dependent responses and EC₅₀ values) agree well with the recommended clinical dose limits for these drugs. Using time-lapse imaging and RNAseq, we found that the differences in anthracycline cardiotoxicity are closely linked to extent of cardiomyocyte uptake and magnitude of activation/inhibition of several cellular pathways such as death receptor signaling, ROS production, and dysregulation of calcium signaling. The results provide molecular insights into anthracycline cardiac interactions and offer a novel assay system to more robustly assess potential cardiotoxicity during drug development.

NGS Evaluation of Colorectal Cancer Reveals Interferon Gamma Dependent Expression of Immune Checkpoint Genes and Identification of Novel IFNγ Induced Genes

To evaluate the expression of immune checkpoint genes, their concordance with expression of IFNγ, and to identify potential novel ICP related genes (ICPRG) in colorectal cancer (CRC), the biological connectivity of six well documented ("classical") ICPs (CTLA4, PD1, PDL1, Tim3, IDO1, and LAG3) with IFNγ and its co-expressed genes was examined by NGS in 79 CRC/healthy colon tissue pairs. Identification of novel IFNγ- induced molecules with potential ICP activity was also sought. In our study, the six classical ICPs were statistically upregulated and correlated with IFNγ, CD8A, CD8B, CD4, and 180 additional immunologically related genes in IFNγ positive (FPKM > 1) tumors. By ICP co-expression analysis, we also identified three IFNγ-induced genes [(IFNγ-inducible lysosomal thiol reductase (IFI30), guanylate binding protein1 (GBP1), and guanylate binding protein 4 (GBP4)] as potential novel ICPRGs. These three genes were upregulated in tumor compared to normal tissues in IFNγ positive tumors, co-expressed with CD8A and had relatively high abundance (average FPKM = 362, 51, and 25, respectively), compared to the abundance of the 5 well-defined ICPs (Tim3, LAG3, PDL1, CTLA4, PD1; average FPKM = 10, 9, 6, 6, and 2, respectively), although IDO1 is expressed at comparably high levels (FPKM = 39). We extended our evaluation by querying the TCGA database which revealed the commonality of IFNγ dependent expression of the three potential ICPRGs in 638 CRCs, 103 skin cutaneous melanomas (SKCM), 1105 breast cancers (BC), 184 esophageal cancers (ESC), 416 stomach cancers (STC), and 501 lung squamous carcinomas (LUSC). In terms of prognosis, based on Pathology Atlas data, correlation of GBP1 and GBP4, but not IFI30, with 5-year survival rate was favorable in CRC, BC, SKCM, and STC. Thus, further studies defining the role of IFI30, GBP1, and GBP4 in CRC are warranted.

Many-to-one binding by intrinsically disordered protein regions

Disordered binding regions (DBRs), which are embedded within intrinsically disordered proteins or regions (IDPs or IDRs), enable IDPs or IDRs to mediate multiple protein-protein interactions. DBR-protein complexes were collected from the Protein Data Bank for which two or more DBRs having different amino acid sequences bind to the same (100% sequence identical) globular protein partner, a type of interaction herein called many-to-one binding. Two distinct binding profiles were identified: independent and overlapping. For the overlapping binding profiles, the distinct DBRs interact by means of almost identical binding sites (herein called "similar"), or the binding sites contain both common and divergent interaction residues (herein called "intersecting"). Further analysis of the sequence and structural differences among these three groups indicate how IDP flexibility allows different segments to adjust to similar, intersecting, and independent binding pockets.

Novel reference genes in colorectal cancer identify a distinct subset of high stage tumors and their associated histologically normal colonic tissues

Background

Reference genes are often interchangeably called housekeeping genes due to 1) the essential cellular functions their proteins provide and 2) their constitutive expression across a range of normal and pathophysiological conditions. However, given the proliferative drive of malignant cells, many reference genes such as beta-actin (ACTB) and glyceraldehyde-3-phosphate-dehydrogenase (GAPDH) which play critical roles in cell membrane organization and glycolysis, may be dysregulated in tumors versus their corresponding normal controls

Methods

Because Next Generation Sequencing (NGS) technology has several advantages over hybridization-based technologies, such as independent detection and quantitation of transcription levels, greater sensitivity, and increased dynamic range, we evaluated colorectal cancers (CRC) and their histologically normal tissue counterparts by NGS to evaluate the expression of 21 “classical” reference genes used as normalization standards for PCR based methods. Seventy-nine paired tissue samples of CRC and their patient matched healthy colonic tissues were subjected to NGS analysis of their mRNAs.

Results

We affirmed that 17 out of 21 classical reference genes had upregulated expression in tumors compared to normal colonic epithelial tissue and dramatically so in some cases. Indeed, tumors were distinguished from normal controls in both unsupervised hierarchical clustering analyses (HCA) and principal component analyses (PCA). We then identified 42 novel potential reference genes with minimal coefficients of variation (CV) across 79 CRC tumor pairs. Though largely consistently expressed across tumors and normal control tissues, a subset of high stage tumors (HSTs) as well as some normal tissue samples (HSNs) located adjacent to these HSTs demonstrated dysregulated expression, thus identifying a subset of tumors with a potentially distinct and aggressive biological profile.

Conclusion

While classical CRC reference genes were found to be differentially expressed between tumors and normal controls, novel reference genes, identified via NGS, were more consistently expressed across malignant and normal colonic tissues. Nonetheless, a subset of HST had profound dysregulation of such genes as did many of the histologically normal tissues adjacent to such HSTs, indicating that the HSTs so distinguished may have unique biological properties and that their histologically normal tissues likely harbor a small population of microscopically undetected but metabolically active tumors.

Electronic supplementary material

The online version of this article (10.1186/s12881-019-0867-y) contains supplementary material, which is available to authorized users.