RNA Sequencing in Comparison to Immunohistochemistry for Measuring Cancer Biomarkers in Breast Cancer and Lung Cancer Specimens

Comparative profiling of surgically resected primary tumors and their lymph node metastases in small-cell lung cancer

BACKGROUND

Profiling studies in small-cell lung cancer (SCLC) have mainly focused on primary tumors, omitting the potential molecular changes that might occur during lymphatic metastasis formation. Here, we assessed the molecular discordance between primary SCLCs and corresponding lymph node (LN) metastases in the light of subtype distribution and expression of clinically relevant proteins.

METHODS

Comparative profiling of 32 surgically resected primary SCLCs and their LN metastases was achieved by RNA expression analysis and immunohistochemistry (IHC). In addition to subtype markers (ASCL1, NEUROD1, POU2F3, and YAP1), the expression of nine cancer-specific proteins was evaluated.

RESULTS

The selected clinically relevant molecules showed no significant differences in their RNA expression profile when assessing the primary tumors and their corresponding LN metastases. Nevertheless, IHC analyses revealed significantly higher DLL3 expression in the primary tumors than in the LN metastases (P = 0.008). In contrast, NEUROD1 expression was significantly lower in the primary tumors (versus LN metastases, P < 0.001). No statistically significant difference was found by IHC analysis in the case of other clinically relevant proteins. Concerning SCLC molecular subtypes, a change in subtype distribution was detected in 21 cases. Phenotype switching from neuroendocrine (NE) subtypes toward non-NE lesions and from non-NE landscape toward NE subtypes were both detected.

CONCLUSIONS

Although the molecular landscape of SCLC LN metastases largely resembles that of the tumor of origin, key differences exist in terms of DLL3 and NEUROD1 expression, and in subtype distribution. These diagnostic pitfalls should be considered when establishing the tumors' molecular profile for future clinical trials solely based on LN biopsies.

Prospective trial of the Oncobox platform RNA sequencing bioinformatic analysis for personalized prescription of targeted drugs

Interrogating gene expression in tumor can identify up- and downregulated molecular targets of cancer drugs. Here we report the results of prospective clinical investigation of using RNA sequencing analysis for personalized cancer therapy. Transcriptomic profiles were analyzed using Oncobox platform that identifies altered expression of drug target genes and molecular pathways and builds a personalized rating of targeted therapeutics. Totally, 239 adult solid cancer patients were enrolled: 135 received cancer drug therapy, others received palliative treatment or radiotherapy, or died before therapy started. Oncobox recommended drugs were prescribed in 59 % of the cases receiving therapy. Otherwise, patients received non-targeted therapy or targeted therapy predicted as inefficient by Oncobox (controls). Patients in the Oncobox group were significantly pre-treated compared to controls, but we observed a longer progression-free survival (PFS) trend in the Oncobox group. Furthermore, post-hoc analysis revealed that time between biopsy collection and tumor profiling significantly impacts Oncobox predictive capacity. Excluding patient cases with biopsy obtained more than 7 months before sequencing lead to a significant difference in PFS between Oncobox and control groups with hazard ratio of 0.45 (p-value = 0.023). These results suggest that transcriptomic profiling provides clinically relevant therapeutic match and can improve disease control rate in solid cancers.

High expression level of ERBB2 and efficacy of trastuzumab deruxtecan in desmoplastic small round cell tumour: a monocentric case series report

BACKGROUND

Desmoplastic small round cell tumours (DSRCTs) represent an ultra-rare subtype of soft tissue sarcoma characterized by a recurrent EWSR1::WT1 oncogenic translocation. Considered as an extremely aggressive cancer, the prognosis remains poor with a median overall survival not exceeding 24-36 months and a 5-year survival <10%.

PATIENTS AND METHODS

We analysed ERBB2/human epidermal growth factor receptor 2 (HER2) expression levels in a series of 13 DSRCT patients, using whole-exome RNA sequencing on formalin-fixed paraffin-embedded samples from a local biopathological database. In addition, a retrospective case series describes the clinical outcome of three successive DSRCT patients treated with trastuzumab deruxtecan (T-DXd).

RESULTS

The gene expression analysis demonstrated a consistent high RNA expression level of ERBB2 in DSRCT, with elevated levels [>5 log2(transcripts per million + 1)] across all samples of the cohort and the expression level was the highest compared with all other sarcoma subtypes. In addition to these results, T-DXd showed a marked activity in all three DSRCT patients who presented with metastatic disease refractory to previous standard chemotherapy. So far, the treatment has been overall well tolerated and is currently pursued in the three patients (duration of response >3 months for all three), which warrants additional investigation.

CONCLUSIONS

This case series presents a major information, suggesting that HER2 is a therapeutic target in DSRCT and T-DXd might represent a novel therapeutic option. Those results require to be rapidly shared with the scientific community and confirmed in a prospective clinical trial in this context of very poor prognosis disease and urgent unmet need.

Detection of fusion events by RNA sequencing in FFPE versus freshly frozen colorectal cancer tissue samples

Gene fusion events result in chimeric proteins that are frequently found in human cancers. Specific targeted therapies are available for several types of cancer fusions including receptor tyrosine kinase gene moieties. RNA sequencing (RNAseq) can directly be used for detection of gene rearrangements in a single test, along with multiple additional biomarkers. However, tumor biosamples are usually formalin-fixed paraffin-embedded (FFPE) tissue blocks where RNA is heavily degraded, which in theory may result in decreased efficiency of fusion detection. Here, for the first time, we compared the efficacy of gene fusion detection by RNAseq for matched pairs of freshly frozen in RNA stabilizing solution (FF) and FFPE tumor tissue samples obtained from 29 human colorectal cancer patients. We detected no statistically significant difference in the number of chimeric transcripts in FFPE and FF RNAseq profiles. The known fusion KANSL1-ARL17A/B occurred with a high frequency in 69% of the patients. We also detected 93 new fusion genes not mentioned in the literature or listed in the ChimerSeq database. Among them, 11 were found in two or more patients, suggesting their potential role in carcinogenesis. Most of the fusions detected most probably represented read-through, microdeletion or local duplication events. Finally, in one patient, we detected a potentially clinically actionable in-frame fusion of LRRFIP2 and ALK genes not previously described in colorectal cancer with an intact tyrosine kinase domain that can be potentially targeted by ALK inhibitors.

Biomarkers of success of anti-PD-(L)1 immunotherapy for non-small cell lung cancer derived from RNA- and whole-exome sequencing: results of a prospective observational study on a cohort of 85 patients

Background

Immune checkpoint inhibitors (ICIs) treatment have shown high efficacy for about 15 cancer types. However, this therapy is only effective in 20-30% of cancer patients. Thus, the precise biomarkers of ICI response are an urgent need.

Methods

We conducted a prospective observational study of the prognostic potential ofseveral existing and putative biomarkers of response to immunotherapy in acohort of 85 patients with lung cancer (LC) receiving PD-1 or PD-L1 targeted ICIs. Tumor biosamples were obtained prior to ICI treatment and profiled by whole exome and RNA sequencing. The entire 403 putative biomarkers were screened, including tumor mutation burden (TMB) and number of cancer neoantigens, 131 specific HLA alleles, homozygous state of 11 HLA alleles and their superfamilies; four gene mutation biomarkers, expression of 45 immune checkpoint genes and closely related genes, and three previously published diagnostic gene signatures; for the first time, activation levels of 188 molecular pathways containing immune checkpoint genes and activation levels of 19 pathways algorithmically generated using a human interactome model centered around immune checkpoint genes. Treatment outcomes and/or progression-free survival (PFS) times were available for 61 of 85 patients with LC, including 24 patients with adenocarcinoma and 27 patients with squamous cell LC, whose samples were further analyzed. For the rest 24 patients, both treatment outcomes and PFS data could not be collected. Of these, 54 patients were treated with PD1-specific and 7 patients with PD-L1-specific ICIs. We evaluated the potential of biomarkers based on PFS and RECIST treatment response data.

Results

In our sample, 45 biomarkers were statistically significantly associated with PFS and 44 with response to treatment, of which eight were shared. Five of these (CD3G and NCAM1 gene expression levels, and levels of activation of Adrenergic signaling in cardiomyocytes, Growth hormone signaling, and Endothelin molecular pathways) were used in our signature that showed an AUC of 0.73 and HR of 0.27 (p=0.00034) on the experimental dataset. This signature was also reliable (AUC 0.76, 0.87) for the independent publicly available LC datasets GSE207422, GSE126044 annotated with ICI response data and demonstrated same survival trends on independent dataset GSE135222 annotated with PFS data. In both experimental and one independent datasets annotated with samples' histotypes, the signature worked better for lung adenocarcinoma than for squamous cell LC.

Conclusion

The high reliability of our signature to predict response and PFS after ICI treatment was demonstrated using experimental and 3 independent datasets. Additionally, annotated molecular profiles obtained in this study were made publicly accessible.

RNA sequencing of formalin fixed paraffin-embedded heart tissue provides transcriptomic information about chemotherapy-induced cardiotoxicity

Gene expression of formalin-fixed paraffin-embedded (FFPE) tissue may serve for molecular studies on cardiovascular diseases. Chemotherapeutics, such as doxorubicin (DOX) may cause heart injury, but the mechanisms of these side effects of DOX are not well understood. This study aimed to investigate whether DOX-induced gene expression in archival FFPE heart tissue in experimental rats would correlate with the gene expression in fresh-frozen heart tissue by applying RNA sequencing technology. The results showed RNA from FFPE samples was degraded, resulting in a lower number of uniquely mapped reads. However, DOX-induced differentially expressed genes in FFPE were related to molecular mechanisms of DOX-induced cardiotoxicity, such as inflammation, calcium binding, endothelial dysfunction, senescence, and cardiac hypertrophy signaling. Our data suggest that, despite the limitations, RNA sequencing of archival FFPE heart tissue supports utilizing FFPE tissues from retrospective studies on cardiovascular disorders, including DOX-induced cardiotoxicity.

Sample-Wise and Gene-Wise Comparisons Confirm a Greater Similarity of RNA and Protein Expression Data at the Level of Molecular Pathways and Suggest an Approach for the Data Quality Check in High-Throughput Expression Databases

Identification of genes and molecular pathways with congruent profiles in the proteomic and transcriptomic datasets may result in the discovery of promising transcriptomic biomarkers that would be more relevant to phenotypic changes. In this study, we conducted comparative analysis of 943 paired RNA and proteomic profiles obtained for the same samples of seven human cancer types from The Cancer Genome Atlas (TCGA) and NCI Clinical Proteomic Tumor Analysis Consortium (CPTAC) [two major open human cancer proteomic and transcriptomic databases] that included 15,112 protein-coding genes and 1611 molecular pathways. Overall, our findings demonstrated statistically significant improvement of the congruence between RNA and proteomic profiles when performing analysis at the level of molecular pathways rather than at the level of individual gene products. Transition to the molecular pathway level of data analysis increased the correlation to 0.19-0.57 (Pearson) and 0.14-057 (Spearman), or 2-3-fold for some cancer types. Evaluating the gain of the correlation upon transition to the data analysis the pathway level can be used to refine the omics data by identifying outliers that can be excluded from the comparison of RNA and proteomic profiles. We suggest using sample- and gene-wise correlations for individual genes and molecular pathways as a measure of quality of RNA/protein paired molecular data. We also provide a database of human genes, molecular pathways, and samples related to the correlation between RNA and protein products to facilitate an exploration of new cancer transcriptomic biomarkers and molecular mechanisms at different levels of human gene expression.

Building flexible and robust analysis frameworks for molecular subtyping of cancers

Molecular subtyping is essential to infer tumor aggressiveness and predict prognosis. In practice, tumor profiling requires in-depth knowledge of bioinformatics tools involved in the processing and analysis of the generated data. Additionally, data incompatibility (e.g., microarray versus RNA sequencing data) and technical and uncharacterized biological variance between training and test data can pose challenges in classifying individual samples. In this article, we provide a roadmap for implementing bioinformatics frameworks for molecular profiling of human cancers in a clinical diagnostic setting. We describe a framework for integrating several methods for quality control, normalization, batch correction, classification and reporting, and develop a use case of the framework in breast cancer.

Human Blood Serum Counteracts EGFR/HER2-Targeted Drug Lapatinib Impact on Squamous Carcinoma SK-BR-3 Cell Growth and Gene Expression

Lapatinib is a targeted therapeutic inhibiting HER2 and EGFR proteins. It is used for the therapy of HER2-positive breast cancer, although not all the patients respond to it. Using human blood serum samples from 14 female donors (separately taken or combined), we found that human blood serum dramatically abolishes the lapatinib-mediated inhibition of growth of the human breast squamous carcinoma SK-BR-3 cell line. This antagonism between lapatinib and human serum was associated with cancelation of the drug induced G1/S cell cycle transition arrest. RNA sequencing revealed 308 differentially expressed genes in the presence of lapatinib. Remarkably, when combined with lapatinib, human blood serum showed the capacity of restoring both the rate of cell growth, and the expression of 96.1% of the genes expression of which were altered by the lapatinib treatment alone. Co-administration of EGF with lapatinib also restores the cell growth and cancels alteration of expression of 95.8% of the genes specific to lapatinib treatment of SK-BR-3 cells. Differential gene expression analysis also showed that in the presence of human serum or EGF, lapatinib was unable to inhibit the Toll-Like Receptor signaling pathway and alter expression of genes linked to the Gene Ontology term of Focal adhesion.

PD-L1 Expression by RNA-Sequencing in Non-Small Cell Lung Cancer: Concordance with Immunohistochemistry and Associations with Pembrolizumab Treatment Outcomes

Programmed cell death ligand (PD-L1) expression by immunohistochemistry (IHC) lacks sensitivity for pembrolizumab immunotherapy selection in non-small cell lung cancer (NSCLC), particularly for tumors with low expression. We retrospectively evaluated transcriptomic PD-L1 by mRNA next-generation sequencing (RNA-seq). In an unselected NSCLC patient cohort (n = 3168) tested during standard care (2017-2021), PD-L1 IHC and RNA-seq demonstrated moderate concordance, with 80% agreement overall. Most discordant cases were either low or negative for PD-L1 expression by IHC but high by RNA-seq. RNA-seq accurately discriminated PD-L1 IHC high from low tumors by receiver operator curve (ROC) analysis but could not distinguish PD-L1 IHC low from negative tumors. In a separate pembrolizumab monotherapy cohort (n = 102), NSCLC tumors classified as PD-L1 high versus not high by RNA-seq had significantly improved response, progression-free survival, and overall survival as an individual measure and in combination with IHC high or low status. PD-L1 IHC status (high or low) trended toward but had no significant associations with improved outcomes. Conventional PD-L1 IHC testing has inherent limitations, making it an imperfect reference standard for evaluating novel testing technologies. RNA-seq offers an objective PD-L1 measure that could represent a complementary method to IHC to improve NSCLC patient selection for immunotherapy.

Human Blood Serum Can Diminish EGFR-Targeted Inhibition of Squamous Carcinoma Cell Growth through Reactivation of MAPK and EGFR Pathways

Regardless of the presence or absence of specific diagnostic mutations, many cancer patients fail to respond to EGFR-targeted therapeutics, and a personalized approach is needed to identify putative (non)responders. We found previously that human peripheral blood and EGF can modulate the activities of EGFR-specific drugs on inhibiting clonogenity in model EGFR-positive A431 squamous carcinoma cells. Here, we report that human serum can dramatically abolish the cell growth rate inhibition by EGFR-specific drugs cetuximab and erlotinib. We show that this phenomenon is linked with derepression of drug-induced G1S cell cycle transition arrest. Furthermore, A431 cell growth inhibition by cetuximab, erlotinib, and EGF correlates with a decreased activity of ERK1/2 proteins. In turn, the EGF- and human serum-mediated rescue of drug-treated A431 cells restores ERK1/2 activity in functional tests. RNA sequencing revealed 1271 and 1566 differentially expressed genes (DEGs) in the presence of cetuximab and erlotinib, respectively. Erlotinib- and cetuximab-specific DEGs significantly overlapped. Interestingly, the expression of 100% and 75% of these DEGs restores to the no-drug level when EGF or a mixed human serum sample, respectively, is added along with cetuximab. In the case of erlotinib, EGF and human serum restore the expression of 39% and 83% of DEGs, respectively. We further assessed differential molecular pathway activation levels and propose that EGF/human serum-mediated A431 resistance to EGFR drugs can be largely explained by reactivation of the MAPK signaling cascade.

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.

Depiction of neuroendocrine features associated with immunotherapy response using a novel one-class predictor in lung adenocarcinoma

BACKGROUND

Tumours with no evidence of neuroendocrine transformation histologically but harbouring neuroendocrine features are collectively referred to as non-small cell lung cancer (NSCLC) with neuroendocrine differentiation (NED). Investigating the mechanisms underlying NED is conducive to designing appropriate treatment options for NSCLC patients.

METHODS

In the present study, we integrated multiple lung cancer datasets to identify neuroendocrine features using a one-class logistic regression (OCLR) machine learning algorithm trained on small cell lung cancer (SCLC) cells, a pulmonary neuroendocrine cell type, based on the transcriptome of NSCLC and named the NED index (NEDI). Single-sample gene set enrichment analysis, pathway enrichment analysis, ESTIMATE algorithm analysis, and unsupervised subclass mapping (SubMap) were performed to assess the altered pathways and immune characteristics of lung cancer samples with different NEDI values.

RESULTS

We developed and validated a novel one-class predictor based on the expression values of 13,279 mRNAs to quantitatively evaluate neuroendocrine features in NSCLC. We observed that a higher NEDI correlated with better prognosis in patients with LUAD. In addition, we observed that a higher NEDI was significantly associated with reduced immune cell infiltration and immune effector molecule expression. Furthermore, we found that etoposide-based chemotherapy might be more effective in the treatment of LUAD with high NEDI values. Moreover, we noted that tumours with low NEDI values had better responses to immunotherapy than those with high NEDI values.

CONCLUSIONS

Our findings improve the understanding of NED and provide a useful strategy for applying NEDI-based risk stratification to guide decision-making in the treatment of LUAD.

Genomic and proteomic biomarker landscape in Clinical Trials

The use of molecular biomarkers to support disease diagnosis, monitor its progression, and guide drug treatment has gained traction in the last decades. While only a dozen biomarkers have been approved for their exploitation in the clinic by the FDA, many more are evaluated in the context of translational research and clinical trials. Furthermore, the information on which biomarkers are measured, for which purpose, and in relation to which conditions are not readily accessible: biomarkers used in clinical studies available through resources such as ClinicalTrials.gov are described as free text, posing significant challenges in finding, analyzing, and processing them by both humans and machines. We present a text mining strategy to identify proteomic and genomic biomarkers used in clinical trials and classify them according to the methodologies by which they are measured. We find more than 3000 biomarkers used in the context of 2600 diseases. By analyzing this dataset, we uncover patterns of use of biomarkers across therapeutic areas over time, including the biomarker type and their specificity. These data are made available at the Clinical Biomarker App at https://www.disgenet.org/biomarkers/, a new portal that enables the exploration of biomarkers extracted from the clinical studies available at ClinicalTrials.gov and enriched with information from the scientific literature. The App features several metrics that assess the specificity of the biomarkers, facilitating their selection and prioritization. Overall, the Clinical Biomarker App is a valuable and timely resource about clinical biomarkers, to accelerate biomarker discovery, development, and application.

Clinically relevant fusion oncogenes: detection and practical implications

Mechanistically, chimeric genes result from DNA rearrangements and include parts of preexisting normal genes combined at the genomic junction site. Some rearranged genes encode pathological proteins with altered molecular functions. Those which can aberrantly promote carcinogenesis are called fusion oncogenes. Their formation is not a rare event in human cancers, and many of them were documented in numerous study reports and in specific databases. They may have various molecular peculiarities like increased stability of an oncogenic part, self-activation of tyrosine kinase receptor moiety, and altered transcriptional regulation activities. Currently, tens of low molecular mass inhibitors are approved in cancers as the drugs targeting receptor tyrosine kinase (RTK) oncogenic fusion proteins, that is, including ALK, ABL, EGFR, FGFR1-3, NTRK1-3, MET, RET, ROS1 moieties. Therein, the presence of the respective RTK fusion in the cancer genome is the diagnostic biomarker for drug prescription. However, identification of such fusion oncogenes is challenging as the breakpoint may arise in multiple sites within the gene, and the exact fusion partner is generally unknown. There is no gold standard method for RTK fusion detection, and many alternative experimental techniques are employed nowadays to solve this issue. Among them, RNA-seq-based methods offer an advantage of unbiased high-throughput analysis of only transcribed RTK fusion genes, and of simultaneous finding both fusion partners in a single RNA-seq read. Here we focus on current knowledge of biology and clinical aspects of RTK fusion genes, related databases, and laboratory detection methods.

Immunohistochemical expression of Drosha is reduced in eutopic and ectopic endometrium of women with adenomyosis

The objective of this study was to evaluate the immunohistochemical expression of Dicer, Drosha, and Exportin-5 in the eutopic and ectopic endometrium of women with adenomyosis. Twenty-two paired ectopic and eutopic endometrium from women with adenomyosis and 10 eutopic endometrium samples from control women undergoing hysterectomy were included in the study. Paraffin-embedded tissue blocks were cut and stained for immunohistochemistry. The percentage of epithelial cells positively marked was identified digitally after an automated slide scanning process. Mann-Whitney test or Wilcoxon signed-rank test was performed for independent and paired groups, respectively. A lower expression of Drosha was observed in the eutopic endometrium of women with adenomyosis than in the eutopic endometrium of women without the disease (69.9±3.4% vs 85.2±2.9%, respectively) (P=0.016; 95%CI: 3.4 to 27.4%). We also detected lower Drosha expression in the ectopic endometrium of women with adenomyosis than in the eutopic endometrium of the same women (59.6±3.2% vs 69.9±3.4%, respectively) (P=0.004; 95%CI: 2.3 to 16.7%). Additionally, we observed a correlation between Drosha expression in the ectopic and paired eutopic endometrium (P=0.034, rho=0.454). No significant difference in Dicer or Exportin expression was observed. Predominant pattern of cytoplasmic staining for the anti-Drosha antibody and both a nuclear and cytoplasmic pattern for the anti-Exportin antibody were observed. Drosha expression was significantly lower in the endometrium of women with adenomyosis compared to the eutopic endometrium of asymptomatic women without the disease. Furthermore, its expression was lower in the ectopic endometrium but correlated to the paired eutopic endometrium.

Quality control recommendations for RNASeq using FFPE samples based on pre-sequencing lab metrics and post-sequencing bioinformatics metrics

Background

Formalin-fixed, paraffin-embedded (FFPE) tissues have many advantages for identification of risk biomarkers, including wide availability and potential for extended follow-up endpoints. However, RNA derived from archival FFPE samples has limited quality. Here we identified parameters that determine which FFPE samples have the potential for successful RNA extraction, library preparation, and generation of usable RNAseq data.

Methods

We optimized library preparation protocols designed for use with FFPE samples using seven FFPE and Fresh Frozen replicate pairs, and tested optimized protocols using a study set of 130 FFPE biopsies from women with benign breast disease. Metrics from RNA extraction and preparation procedures were collected and compared with bioinformatics sequencing summary statistics. Finally, a decision tree model was built to learn the relationship between pre-sequencing lab metrics and qc pass/fail status as determined by bioinformatics metrics.

Results

Samples that failed bioinformatics qc tended to have low median sample-wise correlation within the cohort (Spearman correlation < 0.75), low number of reads mapped to gene regions (< 25 million), or low number of detectable genes (11,400 # of detected genes with TPM > 4). The median RNA concentration and pre-capture library Qubit values for qc failed samples were 18.9 ng/ul and 2.08 ng/ul respectively, which were significantly lower than those of qc pass samples (40.8 ng/ul and 5.82 ng/ul). We built a decision tree model based on input RNA concentration, input library qubit values, and achieved an F score of 0.848 in predicting QC status (pass/fail) of FFPE samples.

Conclusions

We provide a bioinformatics quality control recommendation for FFPE samples from breast tissue by evaluating bioinformatic and sample metrics. Our results suggest a minimum concentration of 25 ng/ul FFPE-extracted RNA for library preparation and 1.7 ng/ul pre-capture library output to achieve adequate RNA-seq data for downstream bioinformatics analysis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-022-01355-0.

Prognostic value of tumor immune biomarkers in biopsies from patients with refractory solid cancers

PD-L1 and tumor-infiltrating lymphocytes play a key role in the immune escape of cancer, although their prognostic value remains unknown in patients with refractory solid cancer compared to other known prognostic estimation methods. In this ancillary study, we assessed the prognostic value of previously-defined prognostic scores (such as the Royal Marsden Hospital (RMH) score) and of PD-L1, CD3, CD8 and FOXP3 expressions based on immunohistochemistry (IHC) and RNA sequencing (RNAseq) of tumor samples from patients included in the personalized-medicine MOSCATO-02 trial. We collected biopsies with successful IHC analysis from 266 patients treated between April 2016 and September 2017, among whom 170 (63.9%) also had a matched RNAseq. We used a Random Forest model to identify the best prognostic factor, and a Lasso-penalized Cox model to validate the findings. We found that the RMH score was the strongest prognostic factor, with high scores associated with a higher risk of death (Hazard Ratio (HR)=1.29; CI95%[1.19-1.21]). The PD-L1 expression score obtained from IHC analyses was the second-best performing predictor, with the 1+ score (low expression) linked to a lower risk of death (HR=0.564; CI95%[0.539-0.580]). Other tested variables, including primary tumor type and subsequent treatments received following biopsy, were not found significantly linked to prognosis. We found modest correlation between IHC and RNAseq expressions of immune genes, but RNAseq related better to prognosis. Overall, our study supports the use of the RMH score and the assessment of PD-L1 expression in IHC to estimate prognosis in patients with advanced cancer.

OncoboxPD: human 51 672 molecular pathways database with tools for activity calculating and visualization

OncoboxPD (Oncobox pathway databank) available at https://open.oncobox.com is the collection of 51 672 uniformly processed human molecular pathways. Superposition of all pathways formed interactome graph of protein–protein interactions and metabolic reactions containing 361 654 interactions and 64 095 molecular participants. Pathways are uniformly classified by biological processes, and each pathway node is algorithmically functionally annotated by specific activator/repressor role. This enables online calculation of statistically supported pathway activation levels (PALs) with the built-in bioinformatic tool using custom RNA/protein expression profiles. Each pathway can be visualized as static or dynamic graph, where vertices are molecules participating in a pathway and edges are interactions or reactions between them. Differentially expressed nodes in a pathway can be visualized in two-color mode with user-defined color scale. For every comparison, OncoboxPD also generates a graph summarizing top up- and downregulated pathways.

CDK6 Immunophenotype Implicates Potential Therapeutic Application of CDK4/6 Inhibitors in Urothelial Carcinoma

Background

Infiltrating bladder urothelial carcinoma is the most common bladder malignancy with limited therapeutic options and poor prognosis. Identifying new therapeutic targets or strategies has important clinical significance. The data from public sources indicate poor prognosis in urothelial carcinoma cases with high CDK6 mRNA levels. Furthermore, studies have shown that CDK6 expression is elevated in urothelial carcinoma tissue compared to the surrounding urothelium, thus presenting a case for performing CDK4/6 inhibitor targeted research in urothelial carcinoma. However, a phase II trial showed that CDK4/6 inhibitors are not effective for advanced urothelial carcinoma, suggesting that case screening is important for targeted therapy.

Objective

Immunohistochemistry (IHC) is simple and easy to perform and can be used to screen urothelial carcinoma cases with high CDK6 expression in clinical practice. The aim of this study was to determine the CDK6 expression threshold for positive cases.

Methods

We evaluated the correlation between the H-score of CDK6 protein expression and survival or CDK6 mRNA level using RNA sequencing. The effects of different CDK4/6 inhibitors were tested on bladder carcinoma cell lines with different CDK6 expression levels.

Results

The H-score, which predicts poor prognosis and reflects a high CDK6 mRNA level, was determined as the selection criterion for positive cases. Furthermore, we found that urothelial carcinoma cell lines with higher CDK6 expression levels displayed greater sensitivity to CDK4/6 inhibitors than cells with lower expression levels.

Conclusions

IHC staining for CDK6 protein in urothelial carcinoma is proposed as a promising screening platform for CDK4/6 inhibitor targeted therapy.

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.

Obtaining Reliable RT-qPCR Results in Molecular Diagnostics—MIQE Goals and Pitfalls for Transcriptional Biomarker Discovery

In this review, we discuss the development pipeline for transcriptional biomarkers in molecular diagnostics and stress the importance of a reliable gene transcript quantification strategy. Hence, a further focus is put on the MIQE guidelines and how to adapt them for biomarker discovery, from signature validation up to routine diagnostic applications. First, the advantages and pitfalls of the holistic RNA sequencing for biomarker development will be described to establish a candidate biomarker signature. Sequentially, the RT-qPCR confirmation process will be discussed to validate the discovered biomarker signature. Examples for the successful application of RT-qPCR as a fast and reproducible quantification method in routinemolecular diagnostics are provided. Based on the MIQE guidelines, the importance of “key steps” in RT-qPCR is accurately described, e.g., reverse transcription, proper reference gene selection and, finally, the application of automated RT-qPCR data analysis software. In conclusion, RT-qPCR proves to be a valuable tool in the establishment of a disease-specific transcriptional biomarker signature and will have a great future in molecular diagnostics or personalized medicine.

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.

Evaluation of PD-L1 and tumor infiltrating lymphocytes in paired pretreatment biopsies and post neoadjuvant chemotherapy surgical specimens of breast carcinoma

Herein it was evaluated the impact of PD-L1 immunohistochemical expression and stromal tumor-infiltrating lymphocyte (sTIL) counts in pretreatment needle core biopsy on response to neoadjuvant chemotherapy (NACT) for patients with breast carcinomas (BC). In 127 paired pre- and post-NACT BC specimens, immunohistochemical expression of PD-L1 was evaluated in stroma and in neoplastic cells. In the same samples sTILs were semi-quantified in tumor stroma. Post-NACT specimens were histologically rated as having residual cancer burden (RCB of any degree), or with complete pathological response (pCR). PD-L1 expression and higher sTIL counts were associated with histological grade 3 BC. PD-L1 expression was also associated with the non-luminal-HER2+ and triple negative immunohistochemical profiles of BC. Pathological complete response was associated with histological grade 3 tumors, and with the non-luminal-HER2+ and triple negative profiles. Additionally, our results support an association between PD-L1 expression and pCR to NACT. It was also observed that there is a trend to reduction of sTIL counts in the post-NACT specimens of patients with pCR. Of note, PD-L1 was expressed in half of the hormone receptor positive cases, a finding that might expand the potential use of immune checkpoint inhibitors for BC patients.

RNA Sequencing Data for FFPE Tumor Blocks Can Be Used for Robust Estimation of Tumor Mutation Burden in Individual Biosamples

Tumor mutation burden (TMB) is a well-known efficacy predictor for checkpoint inhibitor immunotherapies. Currently, TMB assessment relies on DNA sequencing data. Gene expression profiling by RNA sequencing (RNAseq) is another type of analysis that can inform clinical decision-making and including TMB estimation may strongly benefit this approach, especially for the formalin-fixed, paraffin-embedded (FFPE) tissue samples. Here, we for the first time compared TMB levels deduced from whole exome sequencing (WES) and RNAseq profiles of the same FFPE biosamples in single-sample mode. We took TCGA project data with mean sequencing depth 23 million gene-mapped reads (MGMRs) and found 0.46 (Pearson)–0.59 (Spearman) correlation with standard mutation calling pipelines. This was converted into low (<10) and high (>10) TMB per megabase classifier with area under the curve (AUC) 0.757, and application of machine learning increased AUC till 0.854. We then compared 73 experimental pairs of WES and RNAseq profiles with lower (mean 11 MGMRs) and higher (mean 68 MGMRs) RNA sequencing depths. For higher depth, we observed ~1 AUC for the high/low TMB classifier and 0.85 (Pearson)–0.95 (Spearman) correlation with standard mutation calling pipelines. For the lower depth, the AUC was below the high-quality threshold of 0.7. Thus, we conclude that using RNA sequencing of tumor materials from FFPE blocks with enough coverage can afford for high-quality discrimination of tumors with high and low TMB levels in a single-sample mode.

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.

ACE2 & TMPRSS2 Expressions in Head & Neck Tissues: A Systematic Review

To review the data regarding the expression of angiotensin converting enzyme-2 (ACE2) and transmembrane protease serine-2 (TMPRSS2) in head and neck tissue. Scopus, Cochrane Library, Medrxiv, Google Scholar and PubMED/MEDLINE were searched by four independent investigators for studies investigating ACE2 or TMPRSS2 expressions in head and neck tissues. The following outcomes were considered: sample origin (animal versus human); detection method; anatomical location and cell types. PRISMA checklist and modified population, intervention, comparison, outcome, timing and setting (PICOTS) framework were used to perform the review. Of the 24 identified studies, 17 met our inclusion criteria. Thirteen studies were conducted during the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) pandemic. ACE2 and TMPRSS2 were expressed in oral, pharyngeal, sinusonasal human mucosa. The following cell types expressed ACE2: basal, apical, goblet, minor salivary, and endothelial cells. TMPRSS2 was found in goblet and apical respiratory cells. ACE2 and TMPRSS2 were found in the olfactory region, especially in sustentacular non-neural and neural stem cells. Animal studies suggested that ACE2 expression may vary regarding age. There was an important heterogeneity between studies in the methods used to detect ACE2 and TMPRSS2, leading to a potential identification bias. The SARS-CoV-2 receptors, ACE2 and TMPRSS2, are both expressed in many head and neck tissues, enabling the viral entry into the host organism.

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.