Gene Expression and Pathway Activation Biomarkers of Breast Cancer Sensitivity to Taxanes

Taxanes are one of the most widely used classes of breast cancer (BC) therapeutics. Despite the long history of clinical usage, the molecular mechanisms of their action and cancer resistance are still not fully understood. Here we aimed to identify gene expression and molecular pathway activation biomarkers of BC sensitivity to taxane drugs paclitaxel and docetaxel. We used to our knowledge the biggest collection of clinically annotated publicly available literature BC gene expression data (12 datasets, n = 1250) and the experimental clinical BC cohort (n = 12). Seven literature datasets were used for biomarker discovery (n = 914), and the remaining five literature plus one experimental datasets (n = 336) - for the validation. We totally found 34 genes and 29 molecular pathways which could strongly discriminate good and poor responders to taxane treatments. The biomarker genes and pathways were associated with molecular processes related to formation of mitotic spindle and centromeres, and with the spindle assembly mitotic checkpoint. Furthermore, we created gene expression and pathway activation signatures predicting BC response to taxanes. These signatures were tested on the validation BC cohort and demonstrated strong biomarker potential reflected by mean AUC values of 0.76 and 0.77, respectively, which outperforms previously reported analogs. Taken together, these findings can deepen our understanding of mechanism of action of taxanes and potentially improve personalization of treatment in BC.

Pan-cancer experimental characteristic of human transcriptional patterns connected with telomerase reverse transcriptase (TERT) gene expression status

The TERT gene encodes the reverse transcriptase subunit of telomerase and is normally transcriptionally suppressed in differentiated human cells but reactivated in cancers where its expression is frequently associated with poor survival prognosis. Here we experimentally assessed the RNA sequencing expression patterns associated with TERT transcription in 1039 human cancer samples of 27 tumor types. We observed a bimodal distribution of TERT expression where ∼27% of cancer samples did not express TERT and the rest showed a bell-shaped distribution. Expression of TERT strongly correlated with 1443 human genes including 103 encoding transcriptional factor proteins. Comparison of TERT- positive and negative cancers showed the differential activation of 496 genes and 1975 molecular pathways. Therein, 32/38 (84%) of DNA repair pathways were hyperactivated in TERT+ cancers which was also connected with accelerated replication, transcription, translation, and cell cycle progression. In contrast, the level of 40 positive cell cycle regulator proteins and a set of epithelial-to-mesenchymal transition pathways was specific for the TERT- group suggesting different proliferation strategies for both groups of cancer. Our pilot study showed that the TERT+ group had ∼13% of cancers with C228T or C250T mutated TERT promoter. However, the presence of promoter mutations was not associated with greater TERT expression compared with other TERT+ cancers, suggesting parallel mechanisms of its transcriptional activation in cancers. In addition, we detected a decreased expression of L1 retrotransposons in the TERT+ group, and further decreased L1 expression in promoter mutated TERT+ cancers. TERT expression was correlated with 17 genes encoding molecular targets of cancer therapeutics and may relate to differential survival patterns of TERT- positive and negative cancers.

Large-scale assessment of pros and cons of autopsy-derived or tumor-matched tissues as the norms for gene expression analysis in cancers

Normal tissues are essential for studying disease-specific differential gene expression. However, healthy human controls are typically available only in postmortal/autopsy settings. In cancer research, fragments of pathologically normal tissue adjacent to tumor site are frequently used as the controls. However, it is largely underexplored how cancers can systematically influence gene expression of the neighboring tissues. Here we performed a comprehensive pan-cancer comparison of molecular profiles of solid tumor-adjacent and autopsy-derived "healthy" normal tissues. We found a number of systemic molecular differences related to activation of the immune cells, intracellular transport and autophagy, cellular respiration, telomerase activation, p38 signaling, cytoskeleton remodeling, and reorganization of the extracellular matrix. The tumor-adjacent tissues were deficient in apoptotic signaling and negative regulation of cell growth including G2/M cell cycle transition checkpoint. We also detected an extensive rearrangement of the chemical perception network. Molecular targets of 32 and 37 cancer drugs were over- or underexpressed, respectively, in the tumor-adjacent norms. These processes may be driven by molecular events that are correlated between the paired cancer and adjacent normal tissues, that mostly relate to inflammation and regulation of intracellular molecular pathways such as the p38, MAPK, Notch, and IGF1 signaling. However, using a model of macaque postmortal tissues we showed that for the 30 min - 24-hour time frame at 4ºC, an RNA degradation pattern in lung biosamples resulted in an artifact "differential" expression profile for 1140 genes, although no differences could be detected in liver. Thus, such concerns should be addressed in practice.

Association between genome-wide epigenetic and genetic alterations in breast cancer tissue and response to HER2-targeted therapies in HER2-positive breast cancer patients: new findings and a systematic review

Recent evidence suggests that genetic and epigenetic mechanisms might be associated with acquired resistance to cancer therapies. The aim of this study was to assess the association of genome-wide genetic and epigenetic alterations with the response to anti-HER2 agents in HER2-positive breast cancer patients. PubMed was screened for articles published until March 2021 on observational studies investigating the association of genome-wide genetic and epigenetic alterations, measured in breast cancer tissues or blood, with the response to targeted treatment in HER2-positive breast cancer patients. Sixteen studies were included in the review along with ours, in which we compared the genome-wide DNA methylation pattern in breast tumor tissues of patients who acquired resistance to treatment (case group, n = 6) to that of patients who did not develop resistance (control group, n = 6). Among genes identified as differentially methylated between the breast cancer tissue of cases and controls, one of them, PRKACA, was also reported as differentially expressed in two studies included in the review. Although included studies were heterogeneous in terms of methodology and study population, our review suggests that genes of the PI3K pathway may play an important role in developing resistance to anti-HER2 agents in breast cancer patients. Genome-wide genetic and epigenetic alterations measured in breast cancer tissue or blood might be promising markers of resistance to anti-HER2 agents in HER2-positive breast cancer patients. Further studies are needed to confirm these data.

Human Blood Serum Inhibits Ductal Carcinoma Cells BT474 Growth and Modulates Effect of HER2 Inhibition

Trastuzumab, a HER2-targeted antibody, is widely used for targeted therapy of HER2-positive breast cancer (BC) patients; yet, not all of them respond to this treatment. We investigated here whether trastuzumab activity on the growth of HER2-overexpressing BT474 cells may interfere with human peripheral blood endogenous factors. Among 33 individual BC patient blood samples supplemented to the media, BT474 sensitivity to trastuzumab varied up to 14 times. In the absence of trastuzumab, human peripheral blood serum samples could inhibit growth of BT474, and this effect varied ~10 times for 50 individual samples. In turn, the epidermal growth factor (EGF) suppressed the trastuzumab effect on BT474 cell growth. Trastuzumab treatment increased the proportion of BT474 cells in the G0/G1 phases of cell cycle, while simultaneous addition of EGF decreased it, yet not to the control level. We used RNA sequencing profiling of gene expression to elucidate the molecular mechanisms involved in EGF- and human-sera-mediated attenuation of the trastuzumab effect on BT474 cell growth. Bioinformatic analysis of the molecular profiles suggested that trastuzumab acts similarly to the inhibition of PI3K/Akt/mTOR signaling axis, and the mechanism of EGF suppression of trastuzumab activity may be associated with parallel activation of PKC and transcriptional factors ETV1-ETV5.

Experimentally Deduced Criteria for Detection of Clinically Relevant Fusion 3′ Oncogenes from FFPE Bulk RNA Sequencing Data

Drugs targeting receptor tyrosine kinase (RTK) oncogenic fusion proteins demonstrate impressive anti-cancer activities. The fusion presence in the cancer is the respective drug prescription biomarker, but their identification is challenging as both the breakpoint and the exact fusion partners are unknown. RNAseq offers the advantage of finding both fusion parts by screening sequencing reads. Paraffin (FFPE) tissue blocks are the most common way of storing cancer biomaterials in biobanks. However, finding RTK fusions in FFPE samples is challenging as RNA fragments are short and their artifact ligation may appear in sequencing libraries. Here, we annotated RNAseq reads of 764 experimental FFPE solid cancer samples, 96 leukemia samples, and 2 cell lines, and identified 36 putative clinically relevant RTK fusions with junctions corresponding to exon borders of the fusion partners. Where possible, putative fusions were validated by RT-PCR (confirmed for 10/25 fusions tested). For the confirmed 3′RTK fusions, we observed the following distinguishing features. Both moieties were in-frame, and the tyrosine kinase domain was preserved. RTK exon coverage by RNAseq reads upstream of the junction site were lower than downstream. Finally, most of the true fusions were present by more than one RNAseq read. This provides the basis for automatic annotation of 3′RTK fusions using FFPE RNAseq profiles.

Next-Generation Grade and Survival Expression Biomarkers of Human Gliomas Based on Algorithmically Reconstructed Molecular Pathways

In gliomas, expression of certain marker genes is strongly associated with survival and tumor type and often exceeds histological assessments. Using a human interactome model, we algorithmically reconstructed 7494 new-type molecular pathways that are centered each on an individual protein. Each single-gene expression and gene-centric pathway activation was tested as a survival and tumor grade biomarker in gliomas and their diagnostic subgroups (IDH mutant or wild type, IDH mutant with 1p/19q co-deletion, MGMT promoter methylated or unmethylated), including the three major molecular subtypes of glioblastoma (proneural, mesenchymal, classical). We used three datasets from The Cancer Genome Atlas and the Chinese Glioma Genome Atlas, which in total include 527 glioblastoma and 1097 low grade glioma profiles. We identified 2724 such gene and 2418 pathway survival biomarkers out of total 17,717 genes and 7494 pathways analyzed. We then assessed tumor grade and molecular subtype biomarkers and with the threshold of AUC > 0.7 identified 1322/982 gene biomarkers and 472/537 pathway biomarkers. This suggests roughly two times greater efficacy of the reconstructed pathway approach compared to gene biomarkers. Thus, we conclude that activation levels of algorithmically reconstructed gene-centric pathways are a potent class of new-generation diagnostic and prognostic biomarkers for gliomas.

Transcriptomic Portraits and Molecular Pathway Activation Features of Adult Spinal Intramedullary Astrocytomas

In this study, we report 31 spinal intramedullary astrocytoma (SIA) RNA sequencing (RNA-seq) profiles for 25 adult patients with documented clinical annotations. To our knowledge, this is the first clinically annotated RNA-seq dataset of spinal astrocytomas derived from the intradural intramedullary compartment. We compared these tumor profiles with the previous healthy central nervous system (CNS) RNA-seq data for spinal cord and brain and identified SIA-specific gene sets and molecular pathways. Our findings suggest a trend for SIA-upregulated pathways governing interactions with the immune cells and downregulated pathways for the neuronal functioning in the context of normal CNS activity. In two patient tumor biosamples, we identified diagnostic KIAA1549-BRAF fusion oncogenes, and we also found 16 new SIA-associated fusion transcripts. In addition, we bioinformatically simulated activities of targeted cancer drugs in SIA samples and predicted that several tyrosine kinase inhibitory drugs and thalidomide analogs could be potentially effective as second-line treatment agents to aid in the prevention of SIA recurrence and progression.

Gene Expression-Based Signature Can Predict Sorafenib Response in Kidney Cancer

Sorafenib is a tyrosine kinase inhibitory drug with multiple molecular specificities that is approved for clinical use in second-line treatments of metastatic and advanced renal cell carcinomas (RCCs). However, only 10–40% of RCC patients respond on sorafenib-containing therapies, and personalization of its prescription may help in finding an adequate balance of clinical efficiency, cost-effectiveness, and side effects. We investigated whether expression levels of known molecular targets of sorafenib in RCC can serve as prognostic biomarker of treatment response. We used Illumina microarrays to profile RNA expression in pre-treatment formalin-fixed paraffin-embedded (FFPE) samples of 22 metastatic or advanced RCC cases with known responses on next-line sorafenib monotherapy. Among them, nine patients showed partial response (PR), three patients—stable disease (SD), and 10 patients—progressive disease (PD) according to Response Evaluation Criteria In Solid Tumors (RECIST) criteria. We then classified PR + SD patients as “responders” and PD patients as “poor responders”. We found that gene signature including eight sorafenib target genes was congruent with the drug response characteristics and enabled high-quality separation of the responders and poor responders [area under a receiver operating characteristic curve (AUC) 0.89]. We validated these findings on another set of 13 experimental annotated FFPE RCC samples (for 2 PR, 1 SD, and 10 PD patients) that were profiled by RNA sequencing and observed AUC 0.97 for 8-gene signature as the response classifier. We further validated these results in a series of qRT-PCR experiments on the third experimental set of 12 annotated RCC biosamples (for 4 PR, 3 SD, and 5 PD patients), where 8-gene signature showed AUC 0.83.

Experimental and Meta-Analytic Validation of RNA Sequencing Signatures for Predicting Status of Microsatellite Instability

Microsatellite instability (MSI) is an important diagnostic and prognostic cancer biomarker. In colorectal, cervical, ovarian, and gastric cancers, it can guide the prescription of chemotherapy and immunotherapy. In laboratory diagnostics of susceptible tumors, MSI is routinely detected by the size of marker polymerase chain reaction products encompassing frequent microsatellite expansion regions. Alternatively, MSI status is screened indirectly by immunohistochemical interrogation of microsatellite binding proteins. RNA sequencing (RNAseq) profiling is an emerging source of data for a wide spectrum of cancer biomarkers. Recently, three RNAseq-based gene signatures were deduced for establishing MSI status in tumor samples. They had 25, 15, and 14 gene products with only one common gene. However, they were developed and tested on the incomplete literature of The Cancer Genome Atlas (TCGA) sampling and never validated experimentally on independent RNAseq samples. In this study, we, for the first time, systematically validated these three RNAseq MSI signatures on the literature colorectal cancer (CRC) (n = 619), endometrial carcinoma (n = 533), gastric cancer (n = 380), uterine carcinosarcoma (n = 55), and esophageal cancer (n = 83) samples and on the set of experimental CRC RNAseq samples (n = 23) for tumors with known MSI status. We found that all three signatures performed well with area under the curve (AUC) ranges of 0.94-1 for the experimental CRCs and 0.94-1 for the TCGA CRC, esophageal cancer, and uterine carcinosarcoma samples. However, for the TCGA endometrial carcinoma and gastric cancer samples, only two signatures were effective with AUC 0.91-0.97, whereas the third signature showed a significantly lower AUC of 0.69-0.88. Software for calculating these MSI signatures using RNAseq data is included.

RNA Sequencing for Personalized Treatment of Metastatic Leiomyosarcoma: Case Report

Uterine leiomyosarcoma (UL) is a rare malignant tumor that develops from the uterine smooth muscle tissue. Due to the low frequency and lack of sufficient data from clinical trials there is currently no effective treatment that is routinely accepted for UL. Here we report a case of a 65-years-old female patient with metastatic UL, who progressed on ifosfamide and doxorubicin therapy and developed severe hypertensive crisis after administration of second line pazopanib, which lead to treatment termination. Rapid progression of the tumor stressed the need for the alternative treatment options. We performed RNA sequencing and whole exome sequencing profiling of the patient's biopsy and applied Oncobox bioinformatic algorithm to prioritize targeted therapeutics. No clinically relevant mutations associated with drug efficiencies were found, but the Oncobox transcriptome analysis predicted regorafenib as the most effective targeted treatment option. Regorafenib administration resulted in a complete metabolic response which lasted for 10 months. In addition, RNA sequencing analysis revealed a novel cancer fusion transcript of YWHAE gene with fusion partner JAZF1. Several chimeric transcripts for YWHAE and JAZF1 genes were previously found in uterine neoplasms and some of them were associated with tumor prognosis. However, their combination was detected in this study for the first time. Taken together, these findings evidence that RNA sequencing may complement analysis of clinically relevant mutations and enhance management of oncological patients by suggesting putative treatment options.

Large-Scale Transcriptomics-Driven Approach Revealed Overexpression of CRNDE as a Poor Survival Prognosis Biomarker in Glioblastoma

Glioblastoma is the most common and malignant brain malignancy worldwide, with a 10-year survival of only 0.7%. Aggressive multimodal treatment is not enough to increase life expectancy and provide good quality of life for glioblastoma patients. In addition, despite decades of research, there are no established biomarkers for early disease diagnosis and monitoring of patient response to treatment. High throughput sequencing technologies allow for the identification of unique molecules from large clinically annotated datasets. Thus, the aim of our study was to identify significant molecular changes between short- and long-term glioblastoma survivors by transcriptome RNA sequencing profiling, followed by differential pathway-activation-level analysis. We used data from the publicly available repositories The Cancer Genome Atlas (TCGA; number of annotated cases = 135) and Chinese Glioma Genome Atlas (CGGA; number of annotated cases = 218), and experimental clinically annotated glioblastoma tissue samples from the Institute of Pathology, Faculty of Medicine in Ljubljana corresponding to 2-58 months overall survival (n = 16). We found one differential gene for long noncoding RNA CRNDE whose overexpression showed correlation to poor patient OS. Moreover, we identified overlapping sets of congruently regulated differential genes involved in cell growth, division, and migration, structure and dynamics of extracellular matrix, DNA methylation, and regulation through noncoding RNAs. Gene ontology analysis can provide additional information about the function of protein- and nonprotein-coding genes of interest and the processes in which they are involved. In the future, this can shape the design of more targeted therapeutic approaches.

Effect of the Expression of ELOVL5 and IGFBP6 Genes on the Metastatic Potential of Breast Cancer Cells

Breast cancer (BC) is the leading cause of death from malignant neoplasms among women worldwide, and metastatic BC presents the biggest problems for treatment. Previously, it was shown that lower expression of ELOVL5 and IGFBP6 genes is associated with a higher risk of the formation of distant metastases in BC. In this work, we studied the change in phenotypical traits, as well as in the transcriptomic and proteomic profiles of BC cells as a result of the stable knockdown of ELOVL5 and IGFBP6 genes. The knockdown of ELOVL5 and IGFBP6 genes was found to lead to a strong increase in the expression of the matrix metalloproteinase (MMP) MMP1. These results were in good agreement with the correlation analysis of gene expression in tumor samples from patients and were additionally confirmed by zymography. The knockdown of ELOVL5 and IGFBP6 genes was also discovered to change the expression of a group of genes involved in the formation of intercellular contacts. In particular, the expression of the CDH11 gene was markedly reduced, which also complies with the correlation analysis. The spheroid formation assay showed that intercellular adhesion decreased as a result of the knockdown of the ELOVL5 and IGFBP6 genes. Thus, the obtained data indicate that malignant breast tumors with reduced expression of the ELOVL5 and IGFBP6 genes can metastasize with a higher probability due to a more efficient invasion of tumor cells.

Expression of Human Endogenous Retroviruses in Systemic Lupus Erythematosus: Multiomic Integration With Gene Expression

Systemic lupus erythematosus (SLE) is a chronic autoimmune disease characterized by the production of autoantibodies predominantly to nuclear material. Many aspects of disease pathology are mediated by the deposition of nucleic acid containing immune complexes, which also induce the type 1interferon response, a characteristic feature of SLE. Notably, SLE is remarkably heterogeneous, with a variety of organs involved in different individuals, who also show variation in disease severity related to their ancestries. Here, we probed one potential contribution to disease heterogeneity as well as a possible source of immunoreactive nucleic acids by exploring the expression of human endogenous retroviruses (HERVs). We investigated the expression of HERVs in SLE and their potential relationship to SLE features and the expression of biochemical pathways, including the interferon gene signature (IGS). Towards this goal, we analyzed available and new RNA-Seq data from two independent whole blood studies using Telescope. We identified 481 locus specific HERV encoding regions that are differentially expressed between case and control individuals with only 14% overlap of differentially expressed HERVs between these two datasets. We identified significant differences between differentially expressed HERVs and non-differentially expressed HERVs between the two datasets. We also characterized the host differentially expressed genes and tested their association with the differentially expressed HERVs. We found that differentially expressed HERVs were significantly more physically proximal to host differentially expressed genes than non-differentially expressed HERVs. Finally, we capitalized on locus specific resolution of HERV mapping to identify key molecular pathways impacted by differential HERV expression in people with SLE.

Using proteomic and transcriptomic data to assess activation of intracellular molecular pathways

Analysis of molecular pathway activation is the recent instrument that helps to quantize activities of various intracellular signaling, structural, DNA synthesis and repair, and biochemical processes. This may have a deep impact in fundamental research, bioindustry, and medicine. Unlike gene ontology analyses and numerous qualitative methods that can establish whether a pathway is affected in principle, the quantitative approach has the advantage of exactly measuring the extent of a pathway up/downregulation. This results in emergence of a new generation of molecular biomarkers-pathway activation levels, which reflect concentration changes of all measurable pathway components. The input data can be the high-throughput proteomic or transcriptomic profiles, and the output numbers take both positive and negative values and positively reflect overall pathway activation. Due to their nature, the pathway activation levels are more robust biomarkers compared to the individual gene products/protein levels. Here, we review the current knowledge of the quantitative gene expression interrogation methods and their applications for the molecular pathway quantization. We consider enclosed bioinformatic algorithms and their applications for solving real-world problems. Besides a plethora of applications in basic life sciences, the quantitative pathway analysis can improve molecular design and clinical investigations in pharmaceutical industry, can help finding new active biotechnological components and can significantly contribute to the progressive evolution of personalized medicine. In addition to the theoretical principles and concepts, we also propose publicly available software for the use of large-scale protein/RNA expression data to assess the human pathway activation levels.

Algorithmic Annotation of Functional Roles for Components of 3,044 Human Molecular Pathways

Current methods of high-throughput molecular and genomic analyses enabled to reconstruct thousands of human molecular pathways. Knowledge of molecular pathways structure and architecture taken along with the gene expression data can help interrogating the pathway activation levels (PALs) using different bioinformatic algorithms. In turn, the pathway activation profiles can characterize molecular processes, which are differentially regulated and give numeric characteristics of the extent of their activation or inhibition. However, different pathway nodes may have different functions toward overall pathway regulation, and calculation of PAL requires knowledge of molecular function of every node in the pathway in terms of its activator or inhibitory role. Thus, high-throughput annotation of functional roles of pathway nodes is required for the comprehensive analysis of the pathway activation profiles. We proposed an algorithm that identifies functional roles of the pathway components and applied it to annotate 3,044 human molecular pathways extracted from the Biocarta, Reactome, KEGG, Qiagen Pathway Central, NCI, and HumanCYC databases and including 9,022 gene products. The resulting knowledgebase can be applied for the direct calculation of the PALs and establishing large scale profiles of the signaling, metabolic, and DNA repair pathway regulation using high throughput gene expression data. We also provide a bioinformatic tool for PAL data calculations using the current pathway knowledgebase.