Direct formalin fixation induces widespread transcriptomic effects in archival tissue samples

The effect of Proteinase K treatment on GeoMx digital spatial profiling data quality from formalin-fixed, paraffin-embedded tissue

The emergence of spatial profiling technologies in recent years has accelerated opportunities to profile in detail the molecular attributes of a wide range of tissue pathologies using archival specimens. However, tissue treatment for fixation and storage does not always support generation of high-quality genomic data. The purpose of this study was to investigate the impacts of Proteinase K (ProtK) treatment, as a way to increase target transcript exposure, on downstream sequencing data quality metrics for spatial transcriptomic data using formalin-fixed, paraffin-embedded samples. In a series of four independent assessments using different tissue types (nasal mucosa, tonsil, pancreas), varying concentrations of ProtK (ranging from 0.1 to 1 μg/mL) were used as part of the sample processing workflow to generate transcriptomic data using the Nanostring GeoMx DSP and Illumina NextSeq 2000 platforms. Use of higher concentrations of ProtK was generally found to increase total reads (2-4-fold). However, negative probe counts also tended to be increased (2-12-fold), resulting in reductions in the signal-to-noise ratio (10-70% lower) and the number of genes detected above background (50-80% lower). These effects were not seen in all tissues and impacts of tissue handling and processing, beyond ProtK treatment, on data quality metrics, also require consideration. Regardless, these observations highlight the need for careful consideration of a range of sample processing factors and benefits that may be achieved through the optimisation of sample processing workflows for specific tissues as a way to maximise the generation of quality data using spatial transcriptomic approaches.

LINCATRA: Two-cycle method to amplify RNA for transcriptome analysis from formalin-fixed paraffin-embedded tissue

Whole transcriptome analysis (WTA) using RNA extracted from Formalin Fixed Paraffin Embedded (FFPE) tissue is an invaluable tool to understand the molecular pathology of disease. RNA extracted from FFPE tissue is either degraded and/or in very low quantities hampering gene expression analysis. Earlier studies described protocols applied for cellular RNA using poly-A primer-based linear amplification. The current study describes a method, LINCATRA (LINear amplifiCAtion of RNA for whole TRAnscriptome analysis). It employs random nonamer primer based method which can amplify short, fragmented RNA with high fidelity from as low as 5 ng to obtain enough material for WTA. The two-cycle method significantly amplified RNA at ∼1000 folds (p < 0.0001) improving the mean read lengths (p < 0.05) in WTA. Overall, increased mean read length positively correlated with on-target reads (Pearson's r = 0.71, p < 0.0001) in both amplified and unamplified RNA-seq analysis. Gene expression analysis compared between unamplified and amplified group displayed substantial overlap of the differentially expressed genes (DEGs) (log2 fold change cut-off < -2 and >2, p < 0.05) identified between lung cancer and asthma cohorts validating the method developed. This method is applicable in clinical molecular pathology field for both diagnostics and elucidation of disease mechanisms.

Assessment of Gene Set Enrichment Analysis using curated RNA-seq-based benchmarks

Pathway enrichment analysis is a ubiquitous computational biology method to interpret a list of genes (typically derived from the association of large-scale omics data with phenotypes of interest) in terms of higher-level, predefined gene sets that share biological function, chromosomal location, or other common features. Among many tools developed so far, Gene Set Enrichment Analysis (GSEA) stands out as one of the pioneering and most widely used methods. Although originally developed for microarray data, GSEA is nowadays extensively utilized for RNA-seq data analysis. Here, we quantitatively assessed the performance of a variety of GSEA modalities and provide guidance in the practical use of GSEA in RNA-seq experiments. We leveraged harmonized RNA-seq datasets available from The Cancer Genome Atlas (TCGA) in combination with large, curated pathway collections from the Molecular Signatures Database to obtain cancer-type-specific target pathway lists across multiple cancer types. We carried out a detailed analysis of GSEA performance using both gene-set and phenotype permutations combined with four different choices for the Kolmogorov-Smirnov enrichment statistic. Based on our benchmarks, we conclude that the classic/unweighted gene-set permutation approach offered comparable or better sensitivity-vs-specificity tradeoffs across cancer types compared with other, more complex and computationally intensive permutation methods. Finally, we analyzed other large cohorts for thyroid cancer and hepatocellular carcinoma. We utilized a new consensus metric, the Enrichment Evidence Score (EES), which showed a remarkable agreement between pathways identified in TCGA and those from other sources, despite differences in cancer etiology. This finding suggests an EES-based strategy to identify a core set of pathways that may be complemented by an expanded set of pathways for downstream exploratory analysis. This work fills the existing gap in current guidelines and benchmarks for the use of GSEA with RNA-seq data and provides a framework to enable detailed benchmarking of other RNA-seq-based pathway analysis tools.

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.

Transcriptomic profiling of nuclei from paraformaldehyde-fixed and formalin-fixed paraffin-embedded brain tissues

BACKGROUND

Paraformaldehyde (PFA) fixation is necessary for histochemical staining, and formalin-fixed and paraffin-embedded (FFPE) tissue archives are the largest repository of clinically annotated specimens. Single-cell gene expression workflows have recently been developed for PFA-fixed and FFPE tissue specimens. However, for tissues where intact cells are hard to recover, including tissues containing highly interconnected neurons, single-nuclear transcriptomics is beneficial. Moreover, since RNA is very unstable, the effects of standard pathological practice on the transcriptome of samples obtained from such archived specimens like FFPE samples are largely anecdotal.

RESULTS

We evaluated the effects of polyformaldehyde (PFA) fixation and paraffin-embedding on transcriptional profiles of the mouse hippocampus obtained by RNA sequencing (RNA-seq). The transcriptomic signatures of nuclei isolated from fresh PFA-fixed and fresh FFPE tissues were comparable to those of cryopreserved samples. However, more differentially expressed genes were obtained for brains after PFA fixation for more than 3 days than in fresh PFA-fixed samples, especially genes involved in spliceosome and synaptic-related pathways. Importantly, the real cell states were destroyed, with oligodendrocyte precursor cells depleted in the 1day fixed hippocampus. After fixation for 3 days, the proportions of neuronal cells and oligodendrocytes decreased and microglia increased; however, relative frequencies remained constant for longer fixation durations. The storage time of FFPE samples had a negligible effect on the cell composition.

SIGNIFICANCE

This represents the first work to investigate the effects of fixation and storage time of brains on its nuclear transcriptome signatures in detail. The fixation time had more influences on the nuclear transcriptomic profiles than FFPE retention time, and the cliff-like effects appeared to occur over a fixed period of 1-3 days. These findings are expected to guide sample preparation for single-nucleus RNA-seq of FFPE samples, particularly in transcriptomic studies focused on brain diseases.

Tumor biology and immune infiltration define primary liver cancer subsets linked to overall survival after immunotherapy

Primary liver cancer is a rising cause of cancer deaths in the US. Although immunotherapy with immune checkpoint inhibitors induces a potent response in a subset of patients, response rates vary among individuals. Predicting which patients will respond to immune checkpoint inhibitors is of great interest in the field. In a retrospective arm of the National Cancer Institute Cancers of the Liver: Accelerating Research of Immunotherapy by a Transdisciplinary Network (NCI-CLARITY) study, we use archived formalin-fixed, paraffin-embedded samples to profile the transcriptome and genomic alterations among 86 hepatocellular carcinoma and cholangiocarcinoma patients prior to and following immune checkpoint inhibitor treatment. Using supervised and unsupervised approaches, we identify stable molecular subtypes linked to overall survival and distinguished by two axes of aggressive tumor biology and microenvironmental features. Moreover, molecular responses to immune checkpoint inhibitor treatment differ between subtypes. Thus, patients with heterogeneous liver cancer may be stratified by molecular status indicative of treatment response to immune checkpoint inhibitors.

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.

Expression of HLA and Autoimmune Pathway Genes in Liver Biopsies of Young Subjects With Autoimmune Hepatitis Type 1

OBJECTIVES

To test the hypothesis that autoimmune hepatitis (AIH type I) in young subjects is due to genetic differences in proinflammatory genes responding to viral triggers in patients and controls.

METHODS

Intrahepatic gene expression was compared between AIH type I (n = 24, age 9-30 years) patients (hereafter referred to as the AIH group) and controls (n = 21, age 4-25 years). RNA sequencing was performed on complementary DNA (cDNA) libraries made from total RNA extracted from formalin-fixed paraffin-embedded (FFPE) liver biopsy samples. Gene expression levels were quantified, and differentially expressed genes were functionally analyzed. Pathway analysis was performed using the databases Kyoto Encyclopedia of Genes and Genomes (KEGG) and PANTHER. The remaining sequences were mapped to the RefSeq complete set of viral genomes.

RESULTS

Differential gene analysis identified 181 genes that were significantly differentially expressed (136 upregulated in the AIH group). Autoimmune pathway genes such as CD19 and CD20 which are important in B cell regulation and maturation as well as, CD8 and LY9 , which are T-cell related, were upregulated in our AIH group. Genes implicated in AIH pathogenesis including CXCL10 , which is thought to be associated with AIH severity and progression, complement genes ( C1QA, C1QB , and C1QC ), and human leucocyte antigen ( HLA ) genes ( HLA-DRB1, HLA-DRA, HLA-B , and HLA-C ) were upregulated in samples from the AIH group. Specific viral etiologies were not found.

CONCLUSIONS

Unbiased next-generation sequencing and differential gene expression analysis of the AIH group has not only added support for the role of B cells in the pathogenesis and treatment of AIH but also has introduced potential new therapeutic targets: CXCL10 (anti- CXCL10 ) and several complement system-related genes.

Translation complex stabilization on messenger RNA and footprint profiling to study the RNA responses and dynamics of protein biosynthesis in the cells

During protein biosynthesis, ribosomes bind to messenger (m)RNA, locate its protein-coding information, and translate the nucleotide triplets sequentially as codons into the corresponding sequence of amino acids, forming proteins. Non-coding mRNA features, such as 5' and 3' untranslated regions (UTRs), start sites or stop codons of different efficiency, stretches of slower or faster code and nascent polypeptide interactions can alter the translation rates transcript-wise. Most of the homeostatic and signal response pathways of the cells converge on individual mRNA control, as well as alter the global translation output. Among the multitude of approaches to study translational control, one of the most powerful is to infer the locations of translational complexes on mRNA based on the mRNA fragments protected by these complexes from endonucleolytic hydrolysis, or footprints. Translation complex profiling by high-throughput sequencing of the footprints allows to quantify the transcript-wise, as well as global, alterations of translation, and uncover the underlying control mechanisms by attributing footprint locations and sizes to different configurations of the translational complexes. The accuracy of all footprint profiling approaches critically depends on the fidelity of footprint generation and many methods have emerged to preserve certain or multiple configurations of the translational complexes, often in challenging biological material. In this review, a systematic summary of approaches to stabilize translational complexes on mRNA for footprinting is presented and major findings are discussed. Future directions of translation footprint profiling are outlined, focusing on the fidelity and accuracy of inference of the native in vivo translation complex distribution on mRNA.

DNA Methylation Profiling: An Emerging Paradigm for Cancer Diagnosis

Histomorphology has been a mainstay of cancer diagnosis in anatomic pathology for many years. DNA methylation profiling is an additional emerging tool that will serve as an adjunct to increase accuracy of pathological diagnosis. Genome-wide interrogation of DNA methylation signatures, in conjunction with machine learning methods, has allowed for the creation of clinical-grade classifiers, most prominently in central nervous system and soft tissue tumors. Tumor DNA methylation profiling has led to the identification of new entities and the consolidation of morphologically disparate cancers into biologically coherent entities, and it will progressively become mainstream in the future. In addition, DNA methylation patterns in circulating tumor DNA hold great promise for minimally invasive cancer detection and classification. Despite practical challenges that accompany any new technology, methylation profiling is here to stay and will become increasingly utilized as a cancer diagnostic tool across a range of tumor types. Expected final online publication date for the Annual Review of Pathology: Mechanisms of Disease, Volume 17 is January 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.