mRNA-seq whole transcriptome profiling of fresh frozen versus archived fixed tissues

A genetic basis for cancer sex differences revealed in Xp11 translocation renal cell carcinoma

Xp11 translocation renal cell carcinoma (tRCC) is a female-predominant kidney cancer driven by translocations between the TFE3 gene on chromosome Xp11.2 and partner genes located on either chrX or on autosomes. The rearrangement processes that underlie TFE3 fusions, and whether they are linked to the female sex bias of this cancer, are largely unexplored. Moreover, whether oncogenic TFE3 fusions arise from both the active and inactive X chromosomes in females remains unknown. Here we address these questions by haplotype-specific analyses of whole-genome sequences of 29 tRCC samples from 15 patients and by re-analysis of 145 published tRCC whole-exome sequences. We show that TFE3 fusions universally arise as reciprocal translocations with minimal DNA loss or insertion at paired break ends. Strikingly, we observe a near exact 2:1 female:male ratio in TFE3 fusions arising via X:autosomal translocation (but not via X inversion), which accounts for the female predominance of tRCC. This 2:1 ratio is at least partially attributable to oncogenic fusions involving the inactive X chromosome and is accompanied by partial re-activation of silenced chrX genes on the rearranged chromosome. Our results highlight how somatic alterations involving the X chromosome place unique constraints on tumor initiation and exemplify how genetic rearrangements of the sex chromosomes can underlie cancer sex differences.

Diagnostic classification of childhood cancer using multiscale transcriptomics

The causes of pediatric cancers' distinctiveness compared to adult-onset tumors of the same type are not completely clear and not fully explained by their genomes. In this study, we used an optimized multilevel RNA clustering approach to derive molecular definitions for most childhood cancers. Applying this method to 13,313 transcriptomes, we constructed a pediatric cancer atlas to explore age-associated changes. Tumor entities were sometimes unexpectedly grouped due to common lineages, drivers or stemness profiles. Some established entities were divided into subgroups that predicted outcome better than current diagnostic approaches. These definitions account for inter-tumoral and intra-tumoral heterogeneity and have the potential of enabling reproducible, quantifiable diagnostics. As a whole, childhood tumors had more transcriptional diversity than adult tumors, maintaining greater expression flexibility. To apply these insights, we designed an ensemble convolutional neural network classifier. We show that this tool was able to match or clarify the diagnosis for 85% of childhood tumors in a prospective cohort. If further validated, this framework could be extended to derive molecular definitions for all cancer types.

Meta-Analysis of MS-Based Proteomics Studies Indicates Interferon Regulatory Factor 4 and Nucleobindin1 as Potential Prognostic and Drug Resistance Biomarkers in Diffuse Large B Cell Lymphoma

The prognosis of diffuse large B cell lymphoma (DLBCL) is inaccurately predicted using clinical features and immunohistochemistry (IHC) algorithms. Nomination of a panel of molecules as the target for therapy and predicting prognosis in DLBCL is challenging because of the divergences in the results of molecular studies. Mass spectrometry (MS)-based proteomics in the clinic represents an analytical tool with the potential to improve DLBCL diagnosis and prognosis. Previous proteomics studies using MS-based proteomics identified a wide range of proteins. To achieve a consensus, we reviewed MS-based proteomics studies and extracted the most consistently significantly dysregulated proteins. These proteins were then further explored by analyzing data from other omics fields. Among all significantly regulated proteins, interferon regulatory factor 4 (IRF4) was identified as a potential target by proteomics, genomics, and IHC. Moreover, annexinA5 (ANXA5) and nucleobindin1 (NUCB1) were two of the most up-regulated proteins identified in MS studies. Functional enrichment analysis identified the light zone reactions of the germinal center (LZ-GC) together with cytoskeleton locomotion functions as enriched based on consistent, significantly dysregulated proteins. In this study, we suggest IRF4 and NUCB1 proteins as potential biomarkers that deserve further investigation in the field of DLBCL sub-classification and prognosis.

Comparison of whole transcriptome sequencing of fresh, frozen, and formalin-fixed, paraffin-embedded cardiac tissue

The use of fresh tissue for molecular studies is preferred but often impossible. Instead, frozen or formalin-fixed, paraffin-embedded (FFPE) tissues are widely used and constitute valuable resources for retrospective studies. We assessed the utility of cardiac tissue stored in different ways for gene expression analyses by whole transcriptome sequencing of paired fresh, frozen, and FFPE tissues. RNA extracted from FFPE was highly degraded. Sequencing of RNA from FFPE tissues yielded higher proportions of intronic and intergenic reads compared to RNA from fresh and frozen tissues. The global gene expression profiles varied with the storage conditions, particularly mitochondrial and long non-coding RNAs. However, we observed high correlations among protein-coding transcripts (ρ > 0.94) with the various storage conditions. We did not observe any significant storage effect on the allele-specific gene expression. However, FFPE had statistically significantly (p < 0.05) more discordant variant calls compared to fresh and frozen tissue. In conclusion, we found that frozen and FFPE tissues can be used for reliable gene expression analyses, provided that proper quality control is performed and caution regarding the technical variability is withheld.

Improved Bladder Tumor RNA Isolation from Archived Tissues Using Methylene Blue for Normalization, Multiplex RNA Hybridization, Sequencing and Subtyping

Methylene blue (MB) is a dye used for histology with clinical importance and intercalates into nucleic acids. After MB staining of formalin fixed paraffin embedded (FFPE) muscle invasive bladder cancer (MIBC) and normal urothelium, specific regions could be microdissected. It is not known if MB influences RNA used for gene expression studies. Therefore, we analyzed MIBC using five different RNA isolation methods comparing patient matched FFPE and fresh frozen (FF) tissues pre-stained with or without MB. We demonstrate a positive impact of MB on RNA integrity with FF tissues using real time PCR with no interference of its chemical properties. FFPE tissues showed no improvement of RNA integrity, which we propose is due to formalin induced nucleotide crosslinks. Using direct multiplex RNA hybridization the best genes for normalization of MIBC and control tissues were identified from 34 reference genes. In addition, 5SrRNA and 5.8SrRNA were distinctive reference genes detecting <200 bp fragments important for mRNA analyses. Using these normalized RNAs from MB stained MIBC and applying multiplex RNA hybridization and mRNA sequencing, a minimal gene expression panel precisely identified luminal and basal MIBC tumor subtypes, important for diagnosis, prognosis and chemotherapy response.

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.

Determining clinical course of diffuse large B-cell lymphoma using targeted transcriptome and machine learning algorithms

Multiple studies have demonstrated that diffuse large B-cell lymphoma (DLBCL) can be divided into subgroups based on their biology; however, these biological subgroups overlap clinically. Using machine learning, we developed an approach to stratify patients with DLBCL into four subgroups based on survival characteristics. This approach uses data from the targeted transcriptome to predict these survival subgroups. Using the expression levels of 180 genes, our model reliably predicted the four survival subgroups and was validated using independent groups of patients. Multivariate analysis showed that this patient stratification strategy encompasses various biological characteristics of DLBCL, and only TP53 mutations remained an independent prognostic biomarker. This novel approach for stratifying patients with DLBCL, based on the clinical outcome of rituximab, cyclophosphamide, doxorubicin, vincristine, and prednisone therapy, can be used to identify patients who may not respond well to these types of therapy, but would otherwise benefit from alternative therapy and clinical trials.

Analysis of Non-Relapsed and Relapsed Adult Type Granulosa Cell Tumors Suggests Stable Transcriptomes during Tumor Progression

Adult-type granulosa cell tumor (AGCT) is a rare ovarian malignancy characterized by slow growth and hormonal activity. The prognosis of AGCT is generally favorable, but one-third of patients with low-stage disease experience a late relapse, and over half of them die of AGCT. To identify markers that would distinguish patients at risk for relapse, we performed Lexogen QuantSeq 3′ mRNA sequencing on formalin-fixed paraffin-embedded, archival AGCT tissue samples tested positive for the pathognomonic Forkhead Box L2 (FOXL2) mutation. We compared the transcriptomic profiles of 14 non-relapsed archival primary AGCTs (follow-up time 17–26 years after diagnosis) with 13 relapsed primary AGCTs (follow-up time 1.7–18 years) and eight relapsed tumors (follow-up time 2.8–18.9 years). Non-relapsed and relapsed primary AGCTs had similar transcriptomic profiles. In relapsed tumors three genes were differentially expressed: plasmalemma vesicle associated protein (PLVAP) was upregulated (p = 0.01), whereas argininosuccinate synthase 1 (ASS1) (p = 0.01) and perilipin 4 (PLIN4) (p = 0.02) were downregulated. PLVAP upregulation was validated using tissue microarray RNA in situ hybridization. In our patient cohort with extremely long follow-up, we observed similar gene expression patterns in both primary AGCT groups, suggesting that relapse is not driven by transcriptomic changes. These results reinforce earlier findings that molecular markers do not predict AGCT behavior or risk of relapse.

Wholistic approach - transcriptomic analysis and beyond using archival material for molecular diagnosis

Many neoplasms remain unclassified after histopathological examination, which requires further molecular analysis. To this regard, mesenchymal neoplasms are particularly challenging due to the combination of their rarity and the large number of subtypes, and many entities still lack robust diagnostic hallmarks. RNA transcriptomic profiles have proven to be a reliable basis for the classification of previously unclassified tumors and notably for mesenchymal neoplasms. Using exome-based RNA capture sequencing on more than 5000 samples of archival material (FFPE), the combination of expression profiles analyzes (including several clustering methods), fusion genes, and small nucleotide variations has been developed at the Centre Léon Bérard (CLB) in Lyon for the molecular diagnosis of challenging neoplasms and the discovery of new entities. The molecular basis of the technique, the protocol, and the bioinformatics algorithms used are described herein, as well as its advantages and limitations.

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.

Application of the 3' mRNA-Seq using unique molecular identifiers in highly degraded RNA derived from formalin-fixed, paraffin-embedded tissue

Background

Archival formalin-fixed, paraffin-embedded (FFPE) tissue samples with clinical and histological data are a singularly valuable resource for developing new molecular biomarkers. However, transcriptome analysis remains challenging with standard mRNA-seq methods as FFPE derived-RNA samples are often highly modified and fragmented. The recently developed 3′ mRNA-seq method sequences the 3′ region of mRNA using unique molecular identifiers (UMI), thus generating gene expression data with minimal PCR bias. In this study, we evaluated the performance of 3′ mRNA-Seq using Lexogen QuantSeq 3′ mRNA-Seq Library Prep Kit FWD with UMI, comparing with TruSeq Stranded mRNA-Seq and RNA Exome Capture kit. The fresh-frozen (FF) and FFPE tissues yielded nucleotide sizes range from 13 to > 70% of DV200 values; input amounts ranged from 1 ng to 100 ng for validation.

Results

The total mapped reads of QuantSeq 3′ mRNA-Seq to the reference genome ranged from 99 to 74% across all samples. After PCR bias correction, 3 to 56% of total sequenced reads were retained. QuantSeq 3′ mRNA-Seq data showed highly reproducible data across replicates in Universal Human Reference RNA (UHR, R > 0.94) at input amounts from 1 ng to 100 ng, and FF and FFPE paired samples (R = 0.92) at 10 ng. Severely degraded FFPE RNA with ≤30% of DV200 value showed good concordance (R > 0.87) with 100 ng input. A moderate correlation was observed when directly comparing QuantSeq 3′ mRNA-Seq data with TruSeq Stranded mRNA-Seq (R = 0.78) and RNA Exome Capture data (R > 0.67).

Conclusion

In this study, QuantSeq 3′ mRNA-Seq with PCR bias correction using UMI is shown to be a suitable method for gene quantification in both FF and FFPE RNAs. 3′ mRNA-Seq with UMI may be applied to severely degraded RNA from FFPE tissues generating high-quality sequencing data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08068-1.

Downregulated luteolytic pathways in the transcriptome of early pregnancy bovine corpus luteum are mimicked by interferon-tau in vitro

Background

Maintenance of the corpus luteum (CL) beyond the time of luteolysis is essential for establishing pregnancy. Identifying the distinct features of early pregnancy CL remains unresolved, hence we analyzed here the transcriptome of CL on day 18 pregnant (P) and non-pregnant (NP) cows using RNA-Seq. CL of P cows expressed ISGs, verifying exposure to the pregnancy recognition signal, interferon-tau (IFNT), whereas the CL of NP cows had elevated luteal progesterone levels, implying that luteolysis had not yet commenced.

Results

The DEGs, IPA, and metascape canonical pathways, along with GSEA analysis, differed markedly in the CL of P cows from those of NP cows, at the same day of the cycle. Both metascape and IPA identified similar significantly enriched pathways such as interferon alpha/beta, sonic hedgehog pathway, TNFA, EDN1, TGFB1, and PDGF. However, type-1 interferon and sonic hedgehog pathways were positively enriched whereas most of the enriched pathways were downregulated in the P compared to NP samples. Thirty-four % of these pathways are known to be elevated by PGF2A during luteolysis. Notably, selective DEGs in luteinized granulosa cells were modulated by IFNT in vitro in a similar manner to their regulation in the CL of P cows.

Conclusion

This study unraveled the unique transcriptomic signature of the IFNT-exposed, early pregnancy CL, highlighting the abundance of downregulated pathways known to be otherwise induced during luteolysis. These and IFNT-regulated in vitro pregnancy-specific DEGs suggest that IFNT contributes to the characteristics and maintenance of early pregnancy CL.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07747-3.

Comparison of Lesional Juvenile Myositis and Lupus Skin Reveals Overlapping Yet Unique Disease Pathophysiology

OBJECTIVE

Skin inflammation heralds systemic disease in juvenile myositis, yet we lack an understanding of pathogenic mechanisms driving skin inflammation in this disease. We undertook this study to define cutaneous gene expression signatures in juvenile myositis and identify key genes and pathways that differentiate skin disease in juvenile myositis from childhood-onset systemic lupus erythematosus (SLE).

METHODS

We used formalin-fixed paraffin-embedded skin biopsy samples from 15 patients with juvenile myositis (9 lesional, 6 nonlesional), 5 patients with childhood-onset SLE, and 8 controls to perform transcriptomic analysis and identify significantly differentially expressed genes (DEGs; q ≤ 5%) between patient groups. We used Ingenuity Pathway Analysis (IPA) to highlight enriched biologic pathways and validated DEGs by immunohistochemistry and quantitative real-time polymerase chain reaction.

RESULTS

Comparison of lesional juvenile myositis to control samples revealed 221 DEGs, with the majority of up-regulated genes representing interferon (IFN)-stimulated genes. CXCL10, CXCL9, and IFI44L represented the top 3 DEGs (fold change 23.2, 13.3, and 13.0, respectively; q < 0.0001). IPA revealed IFN signaling as the top canonical pathway. When compared to childhood-onset SLE, lesional juvenile myositis skin shared a similar gene expression pattern, with only 28 unique DEGs, including FBLN2, CHKA, and SLURP1. Notably, patients with juvenile myositis who were positive for nuclear matrix protein 2 (NXP-2) autoantibodies exhibited the strongest IFN signature and also demonstrated the most extensive Mx-1 immunostaining, both in keratinocytes and perivascular regions.

CONCLUSION

Lesional juvenile myositis skin demonstrates a striking IFN signature similar to that previously reported in juvenile myositis muscle and peripheral blood. Further investigation into the association of a higher IFN score with NXP-2 autoantibodies may provide insight into disease endotypes and pathogenesis.

Applications of single-cell and bulk RNA sequencing in onco-immunology

The rising interest for precise characterization of the tumour immune contexture has recently brought forward the high potential of RNA sequencing (RNA-seq) in identifying molecular mechanisms engaged in the response to immunotherapy. In this review, we provide an overview of the major principles of single-cell and conventional (bulk) RNA-seq applied to onco-immunology. We describe standard preprocessing and statistical analyses of data obtained from such techniques and highlight some computational challenges relative to the sequencing of individual cells. We notably provide examples of gene expression analyses such as differential expression analysis, dimensionality reduction, clustering and enrichment analysis. Additionally, we used public data sets to exemplify how deconvolution algorithms can identify and quantify multiple immune subpopulations from either bulk or single-cell RNA-seq. We give examples of machine and deep learning models used to predict patient outcomes and treatment effect from high-dimensional data. Finally, we balance the strengths and weaknesses of single-cell and bulk RNA-seq regarding their applications in the clinic.

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.

Two novel qualitative transcriptional signatures robustly applicable to non-research-oriented colorectal cancer samples with low-quality RNA

Currently, due to the low quality of RNA caused by degradation or low abundance, the accuracy of gene expression measurements by transcriptome sequencing (RNA‐seq) is very challenging for non‐research‐oriented clinical samples, majority of which are preserved in hospitals or tissue banks worldwide with complete pathological information and follow‐up data. Molecular signatures consisting of several genes are rarely applied to such samples. To utilize these resources effectively, 45 stage II non‐research‐oriented samples which were formalin‐fixed paraffin‐embedded (FFPE) colorectal carcinoma samples (CRC) using RNA‐seq have been analysed. Our results showed that although gene expression measurements were significantly affected, most cancer features, based on the relative expression orderings (REOs) of gene pairs, were well preserved. We then developed two REO‐based signatures, which consisted of 136 gene pairs for early diagnosis of CRC, and 4500 gene pairs for predicting post‐surgery relapse risk of stage II and III CRC. The performance of our signatures, which included hundreds or thousands of gene pairs, was more robust for non‐research‐oriented clinical samples, compared to that of two published concise REO‐based signatures. In conclusion, REO‐based signatures with relatively more gene pairs could be robustly applied to non‐research‐oriented CRC samples.

Transcriptional Differences between Canine Cutaneous Epitheliotropic Lymphoma and Immune-Mediated Dermatoses

Canine cutaneous epitheliotropic T-cell lymphoma (CETL) and immune-mediated T-cell predominant dermatoses (IMD) share several clinical and histopathological features, but differ substantially in prognosis. The discrimination of ambiguous cases may be challenging, as diagnostic tests are limited and may prove equivocal. This study aimed to investigate transcriptional differences between CETL and IMD, as a basis for further research on discriminating diagnostic biomarkers. We performed 100bp single-end sequencing on RNA extracted from formalin-fixed and paraffin-embedded skin biopsies from dogs with CETL and IMD, respectively. DESeq2 was used for principal component analysis (PCA) and differential gene expression analysis. Genes with significantly different expression were analyzed for enriched pathways using two different tools. The expression of selected genes and their proteins was validated by RT-qPCR and immunohistochemistry. PCA demonstrated the distinct gene expression profiles of CETL and IMD. In total, 503 genes were upregulated, while 4986 were downregulated in CETL compared to IMD. RT-qPCR confirmed the sequencing results for 5/6 selected genes tested, while the protein expression detected by immunohistochemistry was not entirely consistent. Our study revealed transcriptional differences between canine CETL and IMD, with similarities to human cutaneous lymphoma. Differentially expressed genes are potential discriminatory markers, but require further validation on larger sample collections.

LncRNA as Cancer Biomarkers

Although the great majority of cancers share a defined group of hallmarks that is responsible for the uncontrolled growth of particular cell types, it is today clear that under the name of cancer we refer to hundreds of different diseases. Furthermore, each of these diseases has an intrinsic variability due to the genetic background in which it develops. The ability to correctly identify these diseases is urgently needed, because each of them may require a specific therapeutic treatment for successful cure. Cancer biomarkers can be extremely valuable tools for efficient diagnosis and prognosis of cancers. In order to succeed in distinguishing between cancer types and progression-associated genetic backgrounds, cancer biomarkers need to have a strong specificity for a particular disease condition. With the development of novel sequencing technologies, it became clear that the set of genes transcribed from human cells is not limited to genes that code for proteins. On the contrary, our cells contain thousands of RNA without any protein-coding potential. The observation that these transcripts have a much higher cell/tissue specificity of expression in comparison to protein-coding genes makes them a potentially very valuable source of novel cancer biomarkers.

Molecular intrinsic versus clinical subtyping in breast cancer: A comprehensive review

Breast cancer is a heterogeneous disease manifesting diversity at the molecular, histological and clinical level. The development of breast cancer classification was centered on informing clinical decisions. The current approach to the classification of breast cancer, which categorizes this disease into clinical subtypes based on the detection of estrogen receptor, progesterone receptor, human epidermal growth factor receptor 2, and proliferation marker Ki67, is not ideal. This is manifested as a heterogeneity of therapeutic responses and outcomes within the clinical subtypes. The newer classification model, based on gene expression profiling (intrinsic subtyping) informs about transcriptional responses downstream from IHC single markers, revealing deeper appreciation for the disease heterogeneity and capturing tumor biology in a more comprehensive way than an expression of a single protein or gene alone. While accumulating evidences suggest that intrinsic subtypes provide clinically relevant information beyond clinical surrogates, it is imperative to establish whether the current conventional immunohistochemistry-based clinical subtyping approach could be improved by gene expression profiling and if this approach has a potential to translate into clinical practice.

Integration of whole-exome and anchored PCR-based next generation sequencing significantly increases detection of actionable alterations in precision oncology

BACKGROUND

Frequency of clinically relevant mutations in solid tumors by targeted and whole-exome sequencing is ∼30%. Transcriptome analysis complements detection of actionable gene fusions in advanced cancer patients. Goal of this study was to determine the added value of anchored multiplex PCR (AMP)-based next-generation sequencing (NGS) assay to identify further potential drug targets, when coupled with whole-exome sequencing (WES).

METHODS

Selected series of fifty-six samples from 55 patients enrolled in our precision medicine study were interrogated by WES and AMP-based NGS. RNA-seq was performed in 19 cases. Clinically relevant and actionable alterations detected by three methods were integrated and analyzed.

RESULTS

AMP-based NGS detected 48 fusions in 31 samples (55.4%); 31.25% (15/48) were classified as targetable based on published literature. WES revealed 29 samples (51.8%) harbored targetable alterations. TMB-high and MSI-high status were observed in 12.7% and 1.8% of cases. RNA-seq from 19 samples identified 8 targetable fusions (42.1%), also captured by AMP-based NGS. When number of actionable fusions detected by AMP-based NGS were added to WES targetable alterations, 66.1% of samples had potential drug targets. When both WES and RNA-seq were analyzed, 57.8% of samples had targetable alterations.

CONCLUSIONS

This study highlights importance of an integrative genomic approach for precision oncology, including use of different NGS platforms with complementary features. Integrating RNA data (whole transcriptome or AMP-based NGS) significantly enhances detection of potential targets in cancer patients. In absence of fresh frozen tissue, AMP-based NGS is a robust method to detect actionable fusions using low-input RNA from archival tissue.

Human Endogenous Retrovirus Expression Is Upregulated in the Breast Cancer Microenvironment of HIV Infected Women: A Pilot Study

In people living with HIV (PLWH), chronic inflammation can lead to cancer initiation and progression, besides driving a dysregulated and diminished immune responsiveness. HIV infection also leads to increased transcription of Human Endogenous Retroviruses (HERVs), which could increase an inflammatory environment and create a tumor growth suppressive environment with high expression of pro-inflammatory cytokines. In order to determine the impact of HIV infection to HERV expression on the breast cancer microenvironment, we sequenced total RNA from formalin-fixed paraffin-embedded (FFPE) breast cancer samples of women HIV-negative and HIV-positive for transcriptome and retrotranscriptome analyses. We performed RNA extraction from FFPE samples, library preparation and total RNA sequencing (RNA-seq). The RNA-seq analysis shows 185 differentially expressed genes: 181 host genes (178 upregulated and three downregulated) and four upregulated HERV transcripts in HIV-positive samples. We also explored the impact of HERV expression in its neighboring breast cancer development genes (BRCA1, CCND1, NBS1/NBN, RAD50, KRAS, PI3K/PIK3CA) and in long non-coding RNA expression (AC060780.1, also known as RP11-242D8.1). We found a significant positive association of HERV expression with RAD50 and with AC060780.1, which suggest a possible role of HERV in regulating breast cancer genes from PLWH with breast cancer. In addition, we found immune system, extracellular matrix organization and metabolic signaling genes upregulated in HIV-positive breast cancer. In conclusion, our findings provide evidence of transcriptional and retrotranscriptional changes in breast cancer from PLWH compared to non-HIV breast cancer, including dysregulation of HERVs, suggesting an indirect effect of the virus on the breast cancer microenvironment.

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

RNA sequencing is considered the gold standard for high-throughput profiling of gene expression at the transcriptional level. Its increasing importance in cancer research and molecular diagnostics is reflected in the growing number of its mentions in scientific literature and clinical trial reports. However, the use of different reagents and protocols for RNA sequencing often produces incompatible results. Recently, we published the Oncobox Atlas of RNA sequencing profiles for normal human tissues obtained from healthy donors killed in road accidents. This is a database of molecular profiles obtained using uniform protocol and reagents settings that can be broadly used in biomedicine for data normalization in pathology, including cancer. Here, we publish new original 39 breast cancer (BC) and 19 lung cancer (LC) RNA sequencing profiles obtained for formalin-fixed paraffin-embedded (FFPE) tissue samples, fully compatible with the Oncobox Atlas. We performed the first correlation study of RNA sequencing and immunohistochemistry-measured expression profiles for the clinically actionable biomarker genes in FFPE cancer tissue samples. We demonstrated high (Spearman's rho 0.65-0.798) and statistically significant (p < 0.00004) correlations between the RNA sequencing (Oncobox protocol) and immunohistochemical measurements for HER2/ERBB2, ER/ESR1 and PGR genes in BC, and for PDL1 gene in LC; AUC: 0.963 for HER2, 0.921 for ESR1, 0.912 for PGR, and 0.922 for PDL1. To our knowledge, this is the first validation that total RNA sequencing of archived FFPE materials provides a reliable estimation of marker protein levels. These results show that in the future, RNA sequencing can complement immunohistochemistry for reliable measurements of the expression biomarkers in FFPE cancer samples.

Biomarker potential of repetitive-element transcriptome in lung cancer

Since repetitive elements (REs) account for nearly 53% of the human genome, profiling its transcription after an oncogenic change might help in the search for new biomarkers. Lung cancer was selected as target since it is the most frequent cause of cancer death. A bioinformatic workflow based on well-established bioinformatic tools (such as RepEnrich, RepBase, SAMTools, edgeR and DESeq2) has been developed to identify differentially expressed RNAs from REs. It was trained and tested with public RNA-seq data from matched sequencing of tumour and healthy lung tissues from the same patient to reveal differential expression within the RE transcriptome. Healthy lung tissues express a specific set of REs whose expression, after an oncogenic process, is strictly and specifically changed. Discrete sets of differentially expressed REs were found for lung adenocarcinoma, for small-cell lung cancer, and for both cancers. Differential expression affects more HERV-than LINE-derived REs and seems biased towards down-regulation in cancer cells. REs behaving consistently in all patients were tested in a different patient cohort to validate the proposed biomarkers. Down-regulation of AluYg6 and LTR18B was confirmed as potential lung cancer biomarkers, while up-regulation of HERVK11D-Int is specific for lung adenocarcinoma and up-regulation of UCON88 is specific for small cell lung cancer. Hence, the study of RE transcriptome might be considered another research target in cancer, making REs a promising source of lung cancer biomarkers.

The impact of RNA extraction method on accurate RNA sequencing from formalin-fixed paraffin-embedded tissues

Background

Utilization of RNA sequencing methods to measure gene expression from archival formalin-fixed paraffin-embedded (FFPE) tumor samples in translational research and clinical trials requires reliable interpretation of the impact of pre-analytical variables on the data obtained, particularly the methods used to preserve samples and to purify RNA.

Methods

Matched tissue samples from 12 breast cancers were fresh frozen (FF) and preserved in RNAlater or fixed in formalin and processed as FFPE tissue. Total RNA was extracted and purified from FF samples using the Qiagen RNeasy kit, and in duplicate from FFPE tissue sections using three different kits (Norgen, Qiagen and Roche). All RNA samples underwent whole transcriptome RNA sequencing (wtRNAseq) and targeted RNA sequencing for 31 transcripts included in a signature of sensitivity to endocrine therapy. We assessed the effect of RNA extraction kit on the reliability of gene expression levels using linear mixed-effects model analysis, concordance correlation coefficient (CCC) and differential analysis. All protein-coding genes in the wtRNAseq and three gene expression signatures for breast cancer were assessed for concordance.

Results

Despite variable quality of the RNA extracted from FFPE samples by different kits, all had similar concordance of overall gene expression from wtRNAseq between matched FF and FFPE samples (median CCC 0.63–0.66) and between technical replicates (median expression difference 0.13–0.22). More than half of genes were differentially expressed between FF and FFPE, but with low fold change (median |LFC| 0.31–0.34). Two out of three breast cancer signatures studied were highly robust in all samples using any kit, whereas the third signature was similarly discordant irrespective of the kit used. The targeted RNAseq assay was concordant between FFPE and FF samples using any of the kits (CCC 0.91–0.96).

Conclusions

The selection of kit to purify RNA from FFPE did not influence the overall quality of results from wtRNAseq, thus variable reproducibility of gene signatures probably relates to the reliability of individual gene selected and possibly to the algorithm. Targeted RNAseq showed promising performance for clinical deployment of quantitative assays in breast cancer from FFPE samples, although numerical scores were not identical to those from wtRNAseq and would require calibration.

A comparative analysis of RNA sequencing methods with ribosome RNA depletion for degraded and low-input total RNA from formalin-fixed and paraffin-embedded samples

Background

Formalin-fixed and paraffin-embedded (FFPE) blocks held in clinical laboratories are an invaluable resource for clinical research, especially in the era of personalized medicine. It is important to accurately quantitate gene expression with degraded and small amounts of total RNA from FFPE materials.

Results

High concordance in transcript quantifications were shown between FF and FFPE samples using the same kit. The gene expression using the TaKaRa kit showed a difference with other kits, which may be due to the different principle of rRNA depletion or the amount of input total RNA. For seriously degraded RNA from FFPE samples, libraries could be constructed with as low as 50 ng of total RNA, although there was residual rRNA in the libraries. Data analysis with HISAT demonstrated that the unique mapping ratio, percentage of exons in unique mapping reads and number of detected genes decreased along with the decreasing quality of input RNA.

Conclusions

The method of RNA library construction with rRNA depletion can be used for clinical FFPE samples. For degraded and low-input RNA samples, it is still possible to obtain repeatable RNA expression profiling but with a low unique mapping ratio and high residual rRNA.

Genome-wide quantification of ADAR adenosine-to-inosine RNA editing activity

Adenosine-to-inosine (A-to-I) RNA editing by the adenosine deaminase that acts on RNA (ADAR) enzymes is a common RNA modification, preventing false activation of the innate immune system by endogenous double-stranded RNAs. Methods for quantification of ADAR activity are sought after, due to an increasing interest in the role of ADARs in cancer and autoimmune disorders, as well as attempts to harness the ADAR enzymes for RNA engineering. Here, we present the Alu editing index (AEI), a robust and simple-to-use computational tool devised for this purpose. We describe its properties and demonstrate its superiority to current quantification methods of ADAR activity. The AEI is used to map global editing across a large dataset of healthy human samples and identify putative regulators of ADAR, as well as previously unknown factors affecting the observed Alu editing levels. These should be taken into account in future comparative studies of ADAR activity. The AEI tool is available at https://github.com/a2iEditing/RNAEditingIndexer.

RNA sequencing for research and diagnostics in clinical oncology

Molecular diagnostics is becoming one of the major drivers of personalized oncology. With hundreds of different approved anticancer drugs and regimens of their administration, selecting the proper treatment for a patient is at least nontrivial task. This is especially sound for the cases of recurrent and metastatic cancers where the standard lines of therapy failed. Recent trials demonstrated that mutation assays have a strong limitation in personalized selection of therapeutics, consequently, most of the drugs cannot be ranked and only a small percentage of patients can benefit from the screening. Other approaches are, therefore, needed to address a problem of finding proper targeted therapies. The analysis of RNA expression (transcriptomic) profiles presents a reasonable solution because transcriptomics stands a few steps closer to tumor phenotype than the genome analysis. Several recent studies pioneered using transcriptomics for practical oncology and showed truly encouraging clinical results. The possibility of directly measuring of expression levels of molecular drugs' targets and profiling activation of the relevant molecular pathways enables personalized prioritizing for all types of molecular-targeted therapies. RNA sequencing is the most robust tool for the high throughput quantitative transcriptomics. Its use, potentials, and limitations for the clinical oncology will be reviewed here along with the technical aspects such as optimal types of biosamples, RNA sequencing profile normalization, quality controls and several levels of data analysis.

Prognostic modeling of the immune-centric transcriptome reveals interleukin signaling candidates contributing to differential patient outcomes

Immunotherapy is a promising advancement in the treatment of non-small-cell lung carcinoma (NSCLC), although much of how lung tumors interact with the immune system in the natural course of disease remains unknown. We investigated the impact of the expression of immune-centric genes and pathways in tumors on patient survival to reveal novel candidates for immunotherapeutic research. Tumor transcriptomes and detailed clinical characteristics were obtained from patients with NSCLC who were participants of either the Inflammation, Health and Lung Epidemiology (INHALE) (discovery, N = 280) or The Cancer Genome Atlas (TCGA) Lung (replication, N = 1026) studies. Expressions of 2253 genes derived from 48 major immune pathways were assessed for association with patient prognosis using a multivariable Cox model and pathway effects were assessed with an in-house implementation of the Gene Set Enrichment Analysis (GSEA) algorithm. Prognosis-guided gene and pathway analysis of immune-centric expression in tumors revealed significant survival enrichments across both cohorts. The 'Interleukin Signaling' pathway, containing 430 genes, was found to be statistically and significantly enriched with prognostic signal in both the INHALE (P = 0.008) and TCGA (P = 0.039) datasets. Subsequent leading-edge analysis identified a subset of genes (N = 23) shared between both cohorts, driving the pathway enrichment. Cumulative expression of this leading-edge gene signature was a strong predictor of patient survival [discovery: hazard ratio (HR) = 1.59, P = 3.0 × 10-8; replication: HR = 1.29, P = 7.4 × 10-7]. These data demonstrate the impact of immune-centric expression on patient outcomes in NSCLC. Furthermore, prognostic gene effects were localized to discrete immune pathways, of which Interleukin Signaling had the greatest impact on overall survival and the subset of genes driving these effects have promise for future therapeutic intervention.

DegNorm: normalization of generalized transcript degradation improves accuracy in RNA-seq analysis

RNA degradation affects RNA-seq quality when profiling transcriptional activities in cells. Here, we show that transcript degradation is both gene- and sample-specific and is a common and significant factor that may bias the results in RNA-seq analysis. Most existing global normalization approaches are ineffective to correct for degradation bias. We propose a novel pipeline named DegNorm to adjust the read counts for transcript degradation heterogeneity on a gene-by-gene basis while simultaneously controlling for the sequencing depth. The robust and effective performance of this method is demonstrated in an extensive set of simulated and real RNA-seq data.