Gene expression, single nucleotide variant and fusion transcript discovery in archival material from breast tumors

Systematic review and feasibility study on pre-analytical factors and genomic analyses on archival formalin-fixed paraffin-embedded breast cancer tissue

Formalin-fixed paraffin-embedded (FFPE) tissue represents a valuable source for translational cancer research. However, the widespread application of various downstream methods remains challenging. Here, we aimed to assess the feasibility of a genomic and gene expression analysis workflow using FFPE breast cancer (BC) tissue. We conducted a systematic literature review for the assessment of concordance between FFPE and fresh-frozen matched tissue samples derived from patients with BC for DNA and RNA downstream applications. The analytical performance of three different nucleic acid extraction kits on FFPE BC clinical samples was compared. We also applied a newly developed targeted DNA Next-Generation Sequencing (NGS) 370-gene panel and the nCounter BC360® platform on simultaneously extracted DNA and RNA, respectively, using FFPE tissue from a phase II clinical trial. Of the 3701 initial search results, 40 articles were included in the systematic review. High degree of concordance was observed in various downstream application platforms. Moreover, the performance of simultaneous DNA/RNA extraction kit was demonstrated with targeted DNA NGS and gene expression profiling. Exclusion of variants below 5% variant allele frequency was essential to overcome FFPE-induced artefacts. Targeted genomic analyses were feasible in simultaneously extracted DNA/RNA from FFPE material, providing insights for their implementation in clinical trials/cohorts.

Implementation of a High-Accuracy Targeted Gene Expression Panel for Clinical Care

This study describes the validation of a clinical RNA expression panel with evaluation of concordance between gene copy gain by a next-generation sequencing (NGS) assay and high gene expression by an RNA expression panel. The RNA Salah Targeted Expression Panel (RNA STEP) was designed with input from oncologists to include 204 genes with utility for clinical trial prescreening and therapy selection. RNA STEP was validated with the nanoString platform using remnant formalin-fixed, paraffin-embedded-derived RNA from 102 patients previously tested with a validated clinical NGS panel. The repeatability, reproducibility, and concordance of RNA STEP results with NGS results were evaluated. RNA STEP demonstrated high repeatability and reproducibility, with excellent correlation (r > 0.97, P < 0.0001) for all comparisons. Comparison of RNA STEP high gene expression (log2 ratio ≥ 2) versus NGS DNA-based gene copy number gain (copies ≥ 5) for 38 mutually covered genes revealed an accuracy of 93.0% with a positive percentage agreement of 69.4% and negative percentage agreement of 93.8%. Moderate correlation was observed between platforms (r = 0.53, P < 0.0001). Concordance between high gene expression and gene copy number gain varied by specific gene, and some genes had higher accuracy between assays. Clinical implementation of RNA STEP provides gene expression data complementary to NGS and offers a tool for prescreening patients for clinical trials.

Application of the NanoString nCounter System as an Alternative Method to Investigate Molecular Mechanisms Involved in Host Plant Responses to Plasmodiophora brassicae

Clubroot, caused by the soilborne pathogen Plasmodiophora brassicae, is an important disease of canola (Brassica napus) and other crucifers. The recent application of RNA sequencing (RNA-seq) technologies to study P. brassicae−host interactions has generated large amounts of gene expression data, improving knowledge of the molecular mechanisms of pathogenesis and host resistance. Quantitative PCR (qPCR) analysis has been widely applied to examine the expression of a limited number of genes and to validate the results of RNA-seq studies, but may not be ideal for analyzing larger suites of target genes or increased sample numbers. Moreover, the need for intermediate steps such as cDNA synthesis may introduce variability that could affect the accuracy of the data generated by qPCR. Here, we report the validation of gene expression data from a previous RNA-seq study of clubroot using the NanoString nCounter System, which achieves efficient gene expression quantification in a fast and simple manner. We first confirm the robustness of the NanoString system by comparing the results with those generated by qPCR and RNA-seq and then discuss the importance of some candidate genes for resistance or susceptibility to P. brassicae in the host. The results show that the expression of genes measured using NanoString have a high correlation with the values obtained using the other two technologies, with R > 0.90 and p < 0.01, and the same expression patterns for most genes. The three methods (qPCR, RNA-seq, and NanoString) were also compared in terms of laboratory procedures, time, and cost. We propose that the NanoString nCounter System is a robust, sensitive, highly reproducible, and simple technology for gene expression analysis. NanoString could become a common alternative to qPCR to validate RNA-seq data or to create panels of genes for use as markers of resistance/susceptibility when plants are challenged with different P. brassicae pathotypes.

Laying the groundwork for the Biobank of Rare Malignant Neoplasms at the service of the Hellenic Network of Precision Medicine on Cancer

Biobanks constitute an integral part of precision medicine. They provide a repository of biospecimens that may be used to elucidate the pathophysiology, support diagnoses, and guide the treatment of diseases. The pilot biobank of rare malignant neoplasms has been established in the context of the Hellenic Network of Precision Medicine on Cancer and aims to enhance future clinical and/or research studies in Greece by collecting, processing, and storing rare malignant neoplasm samples with associated data. The biobank currently comprises 553 samples; 384 samples of hematopoietic and lymphoid tissue malignancies, 72 samples of pediatric brain tumors and 97 samples of malignant skin neoplasms. In this article, sample collections and their individual significance in clinical research are described in detail along with computational methods developed specifically for this project. A concise review of the Greek biobanking landscape is also delineated, in addition to recommended technologies, methodologies and protocols that were integrated during the creation of the biobank. This project is expected to re‑enforce current clinical and research studies, introduce advances in clinical and genetic research and potentially aid in future targeted drug discovery. It is our belief that the future of medical research is entwined with accessible, effective, and ethical biobanking and that our project will facilitate research planning in the '‑omic' era by contributing high‑quality samples along with their associated data.

Therapeutic Intervention Using a Smad7-Based Tat Protein to Treat Radiation-Induced Oral Mucositis

PURPOSE

Recent studies reported therapeutic effects of Smad7 on oral mucositis in mice without compromising radiation therapy-induced cancer cell killing in neighboring oral cancer. This study aims to assess whether a Smad7-based biologic can treat oral mucositis in a clinically relevant setting by establishing an oral mucositis model in dogs and analyzing molecular targets.

METHODS AND MATERIALS

We created a truncated human Smad7 protein fused with the cell-penetrating Tat tag (Tat-PYC-Smad7). We used intensity modulated radiation therapy to induce oral mucositis in dogs and applied Tat-PYC-Smad7 to the oral mucosa in dose-finding studies after intensity modulated radiation therapy. Clinical outcomes were evaluated. Molecular targets were analyzed in biopsies and serum samples.

RESULTS

Tat-PYC-Smad7 treatment significantly shortened the duration of grade 3 oral mucositis based on double-blinded Veterinary Radiation Therapy Oncology Group scores and histopathology evaluations. Topically applied Tat-PYC-Smad7 primarily penetrated epithelial cells and was undetectable in serum. NanoString nCounter Canine IO Panel identified that, compared to the vehicle samples, top molecular changes in Tat-PYC-Smad7 treated samples include reductions in inflammation and cell death and increases in cell growth and DNA repair. Consistently, immunostaining shows that Tat-PYC-Smad7 reduced DNA damage and neutrophil infiltration with attenuated TGF-β and NFκB signaling. Furthermore, IL-1β and TNF-α were lower in Tat-PYC-Smad7 treated mucosa and serum samples compared to those in vehicle controls.

CONCLUSIONS

Topical Tat-PYC-Smad7 application demonstrated therapeutic effects on oral mucositis induced by intensity modulated radiation therapy in dogs. The local effects of Tat-PYC-Smad7 targeted molecules involved in oral mucositis pathogenesis as well as reduced systemic inflammatory cytokines.

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.

Gene Expression Profiling of FFPE Samples: A Titration Test

The gene expression analysis of formalin-fixed paraffin-embedded (FFPE) tissues is often hampered by poor RNA quality, which results from the oxidation, cross-linking and other chemical modifications induced by the inclusion in paraffin. Yet, FFPE samples are a valuable source for molecular studies and can provide great insights into disease progression and prognosis. With the advancement of genomic technologies, new methods have been established that offer reliable and accurate gene expression workflows on samples of poor quality. NanoString is a probe-based technology that allows the direct counting of the mRNA transcripts and can be applied to degraded samples. Here, we have tested 2 RNA extraction methods for FFPE samples, and we have performed a titration experiment to evaluate the impact of RNA degradation and RNA input on the gene expression profiles assessed using the NanoString IO360 panel. We have selected FFPE samples of different DV200 values and assessed them on the nCounter platform with 2 different amounts of input RNA. This study concludes that the nCounter is a robust and reliable platform to assess the gene expression of RNA samples with DV200 > 30%; its robustness and ease of use could be of particular benefit to clinical settings.

Nanostring-Based Identification of the Gene Expression Profile in Trigger Finger Samples

Trigger finger is a common yet vastly understudied fibroproliferative hand pathology, severely affecting patients' quality of life. Consistent trauma due to inadequate positioning within the afflicted finger's tendon/pulley system leads to cellular dysregulation and eventual fibrosis. While the genetic characteristics of the fibrotic tissue in the trigger finger have been studied, the pathways that govern the initiation and propagation of fibrosis are still unknown. The complete gene expression profile of the trigger finger has never been explored. Our study has used the Nanostring nCounter gene expression assay to investigate the molecular signaling involved in trigger finger pathogenesis. We collected samples from patients undergoing trigger finger (n = 4) release surgery and compared the gene expression to carpal tunnel tissue (n = 4). Nanostring nCounter analysis identified 165 genes that were differentially regulated; 145 of these genes were upregulated, whereas 20 genes were downregulated. We found that several collagen genes were significantly upregulated, and a regulatory matrix metalloproteinase (MMP), MMP-3, was downregulated. Bioinformatic analysis revealed that several known signaling pathways were dysregulated, such as the TGF-β1 and Wnt signaling pathways. We also found several novel signaling pathways (e.g., PI3K, MAPK, JAK-STAT, and Notch) differentially regulated in trigger finger. The outcome of our study helps in understanding the molecular signaling pathway involved in the pathogenesis of the trigger finger.

Transcriptome sequencing of archived lymphoma specimens is feasible and clinically relevant using exome capture technology

Formalin-fixed, paraffin-embedded (FFPE) specimens are an underutilized resource in medical research, particularly in the setting of transcriptome sequencing, as RNA from these samples is often degraded. We took advantage of an exome capture-based RNA-sequencing protocol to explore global gene expression in paired fresh-frozen (FF) and FFPE samples from 16 diffuse large B-cell lymphoma (DLBCL) patients. While FFPE samples generated fewer mapped reads compared to their FF counterparts, these reads captured the same library complexity and had a similar number of genes expressed on average. Furthermore, gene expression demonstrated a high correlation when comparing housekeeping genes only or across the entire transcriptome (r = 0.99 for both comparisons). Differences in gene expression were primarily seen in lowly expressed genes and genes with small or large coding sequences. Using cell-of-origin classifiers and clinically relevant gene expression signatures for DLBCL, FF, and FFPE samples from the same biopsy paired nearly perfectly in clustering analysis. This was further confirmed in a validation cohort of 50 FFPE DLBCL samples. In summary, we found the biological differences between tumors to be far greater than artifacts created as a result of degraded RNA. We conclude that exome capture transcriptome sequencing data from archival samples can confidently be used for cell-of-origin classification of DLBCL samples.

Transcriptomic analysis to identify genes associated with selective hippocampal vulnerability in Alzheimer's disease

Selective vulnerability of different brain regions is seen in many neurodegenerative disorders. The hippocampus and cortex are selectively vulnerable in Alzheimer's disease (AD), however the degree of involvement of the different brain regions differs among patients. We classified corticolimbic patterns of neurofibrillary tangles in postmortem tissue to capture extreme and representative phenotypes. We combined bulk RNA sequencing with digital pathology to examine hippocampal vulnerability in AD. We identified hippocampal gene expression changes associated with hippocampal vulnerability and used machine learning to identify genes that were associated with AD neuropathology, including SERPINA5, RYBP, SLC38A2, FEM1B, and PYDC1. Further histologic and biochemical analyses suggested SERPINA5 expression is associated with tau expression in the brain. Our study highlights the importance of embracing heterogeneity of the human brain in disease to identify disease-relevant gene expression.

MethCORR infers gene expression from DNA methylation and allows molecular analysis of ten common cancer types using fresh-frozen and formalin-fixed paraffin-embedded tumor samples

Background

Transcriptional analysis is widely used to study the molecular biology of cancer and hold great biomarker potential for clinical patient stratification. Yet, accurate transcriptional profiling requires RNA of a high quality, which often cannot be retrieved from formalin-fixed, paraffin-embedded (FFPE) tumor tissue that is routinely collected and archived in clinical departments. To overcome this roadblock to clinical testing, we previously developed MethCORR, a method that infers gene expression from DNA methylation data, which is robustly retrieved from FFPE tissue. MethCORR was originally developed for colorectal cancer and with this study, we aim to: (1) extend the MethCORR method to 10 additional cancer types and (2) to illustrate that the inferred gene expression is accurate and clinically informative.

Results

Regression models to infer gene expression information from DNA methylation were developed for ten common cancer types using matched RNA sequencing and DNA methylation profiles (HumanMethylation450 BeadChip) from The Cancer Genome Atlas Project. Robust and accurate gene expression profiles were inferred for all cancer types: on average, the expression of 11,000 genes was modeled with good accuracy and an intra-sample correlation of R² = 0.90 between inferred and measured gene expression was observed. Molecular pathway analysis and transcriptional subtyping were performed for breast, prostate, and lung cancer samples to illustrate the general usability of the inferred gene expression profiles: overall, a high correlation of r = 0.96 (Pearson) in pathway enrichment scores and a 76% correspondence in molecular subtype calls were observed when using measured and inferred gene expression as input. Finally, inferred expression from FFPE tissue correlated better with RNA sequencing data from matched fresh-frozen tissue than did RNA sequencing data from FFPE tissue (P < 0.0001; Wilcoxon rank-sum test).

Conclusions

In all cancers investigated, MethCORR enabled DNA methylation-based transcriptional analysis, thus enabling future analysis of cancer in situations where high-quality DNA, but not RNA, is available. Here, we provide the framework and resources for MethCORR modeling of ten common cancer types, thereby widely expanding the possibilities for transcriptional studies of archival FFPE material.

3' MACE RNA-sequencing allows for transcriptome profiling in human tissue samples after long-term storage

This study aims to compare the potential of standard RNA-sequencing (RNA-Seq) and 3' massive analysis of c-DNA ends (MACE) RNA-sequencing for the analysis of fresh tissue and describes transcriptome profiling of formalin-fixed paraffin-embedded (FFPE) archival human samples by MACE. To compare MACE to standard RNA-Seq on fresh tissue, four healthy conjunctiva from four subjects were collected during vitreoretinal surgery, halved and immediately transferred to RNA lysis buffer without prior fixation and then processed for either standard RNA-Seq or MACE RNA-Seq analysis. To assess the impact of FFPE preparation on MACE, a third part was fixed in formalin and processed for paraffin embedding, and its transcriptional profile was compared with the unfixed specimens analyzed by MACE. To investigate the impact of FFPE storage time on MACE results, 24 FFPE-treated conjunctival samples from 24 patients were analyzed as well. Nineteen thousand six hundred fifty-nine transcribed genes were detected by both MACE and standard RNA-Seq on fresh tissue, while 3251 and 2213 transcripts were identified explicitly by MACE or RNA-Seq, respectively. Standard RNA-Seq tended to yield longer detected transcripts more often than MACE technology despite normalization, indicating that the MACE technology is less susceptible to a length bias. FFPE processing revealed negligible effects on MACE sequencing results. Several quality-control measurements showed that long-term storage in paraffin did not decrease the diversity of MACE libraries. We noted a nonlinear relation between storage time and the number of raw reads with an accelerated decrease within the first 1000 days in paraffin, while the numbers remained relatively stable in older samples. Interestingly, the number of transcribed genes detected was independent on FFPE storage time. RNA of sufficient quality and quantity can be extracted from FFPE samples to obtain comprehensive transcriptome profiling using MACE technology. We thus present MACE as a novel opportunity for utilizing FFPE samples stored in histological archives.

MethCORR modelling of methylomes from formalin-fixed paraffin-embedded tissue enables characterization and prognostication of colorectal cancer

Transcriptional characterization and classification has potential to resolve the inter-tumor heterogeneity of colorectal cancer and improve patient management. Yet, robust transcriptional profiling is difficult using formalin-fixed, paraffin-embedded (FFPE) samples, which complicates testing in clinical and archival material. We present MethCORR, an approach that allows uniform molecular characterization and classification of fresh-frozen and FFPE samples. MethCORR identifies genome-wide correlations between RNA expression and DNA methylation in fresh-frozen samples. This information is used to infer gene expression information in FFPE samples from their methylation profiles. MethCORR is here applied to methylation profiles from 877 fresh-frozen/FFPE samples and comparative analysis identifies the same two subtypes in four independent cohorts. Furthermore, subtype-specific prognostic biomarkers that better predicts relapse-free survival (HR = 2.66, 95%CI [1.67-4.22], P value < 0.001 (log-rank test)) than UICC tumor, node, metastasis (TNM) staging and microsatellite instability status are identified and validated using DNA methylation-specific PCR. The MethCORR approach is general, and may be similarly successful for other cancer types.

RNA-seq from archival FFPE breast cancer samples: molecular pathway fidelity and novel discovery

BACKGROUND

Formalin-fixed, paraffin-embedded (FFPE) tissues for RNA-seq have advantages over fresh frozen tissue including abundance and availability, connection to rich clinical data, and association with patient outcomes. However, FFPE-derived RNA is highly degraded and chemically modified, which impacts its utility as a faithful source for biological inquiry.

METHODS

True archival FFPE breast cancer cases (n = 58), stored at room temperature for 2-23 years, were utilized to identify key steps in tissue selection, RNA isolation, and library choice. Gene expression fidelity was evaluated by comparing FFPE data to public data obtained from fresh tissues, and by employing single-gene, gene set and transcription network-based regulon analyses.

RESULTS

We report a single 10 μm section of breast tissue yields sufficient RNA for RNA-seq, and a relationship between RNA quality and block age that was not linear. We find single-gene analysis is limiting with FFPE tissues, while targeted gene set approaches effectively distinguish ER+ from ER- breast cancers. Novel utilization of regulon analysis identified the transcription factor KDM4B to associate with ER+ disease, with KDM4B regulon activity and gene expression having prognostic significance in an independent cohort of ER+ cases.

CONCLUSION

Our results, which outline a robust FFPE-RNA-seq pipeline for broad use, support utilizing FFPE tissues to address key questions in the breast cancer field, including the delineation between indolent and life-threatening disease, biological stratification and molecular mechanisms of treatment resistance.

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.

An expanded landscape of human long noncoding RNA

Long noncoding RNAs (lncRNAs) are emerging as key regulators of multiple essential biological processes involved in physiology and pathology. By analyzing the largest compendium of 14,166 samples from normal and tumor tissues, we significantly expand the landscape of human long noncoding RNA with a high-quality atlas: RefLnc (Reference catalog of LncRNA). Powered by comprehensive annotation across multiple sources, RefLnc helps to pinpoint 275 novel intergenic lncRNAs correlated with sex, age or race as well as 369 novel ones associated with patient survival, clinical stage, tumor metastasis or recurrence. Integrated in a user-friendly online portal, the expanded catalog of human lncRNAs provides a valuable resource for investigating lncRNA function in both human biology and cancer development.

Pathway activity profiling of growth factor receptor network and stemness pathways differentiates metaplastic breast cancer histological subtypes

BACKGROUND

Gene expression profiling of rare cancers has proven challenging due to limited access to patient materials and requirement of intact, non-degraded RNA for next-generation sequencing. We customized a gene expression panel compatible with degraded RNA from formalin-fixed, paraffin-embedded (FFPE) patient cancer samples and investigated its utility in pathway activity profiling in patients with metaplastic breast cancer (MpBC).

METHODS

Activity of various biological pathways was profiled in samples from nineteen patients with MpBC and 8 patients with invasive ductal carcinoma with triple negative breast cancer (TNBC) phenotype using a custom gene expression-based assay of 345 genes.

RESULTS

MpBC samples of mesenchymal (chondroid and/or osteoid) histology demonstrated increased SNAI1 and BCL2L11 pathway activity compared to samples with non-mesenchymal histology. Additionally, late cornified envelope and keratinization genes were downregulated in MpBC compared to TNBC, and epithelial-to-mesenchymal transition (EMT) and collagen genes were upregulated in MpBC. Patients with high activity of an invasiveness gene expression signature, as well as high expression of the mesenchymal marker and extracellular matrix glycoprotein gene SPARC, experienced worse outcomes than those with low invasiveness activity and low SPARC expression.

CONCLUSIONS

This study demonstrates the utility of gene expression profiling of metaplastic breast cancer FFPE samples with a custom counts-based assay. Gene expression patterns identified by this assay suggest that, although often histologically triple negative, patients with MpBC have distinct pathway activation compared to patients with invasive ductal TNBC. Incorporation of targeted therapies may lead to improved outcome for MpBC patients, especially in those patients expressing increased activity of invasiveness pathways.

MYB-NFIB gene fusions identified in archival adenoid cystic carcinoma tissue employing NanoString analysis: an exploratory study

BACKGROUND

Adenoid cystic carcinoma (ACC) is a slow growing salivary gland malignancy that is molecularly characterized by t(6:9)(q22-23;p23-24) translocations which predominantly result in MYB-NFIB gene fusions in nearly half of tumours. Detection of MYB-NFIB transcripts is typically performed with fresh ACC tissue using conventional RT-PCR fragment analysis or FISH techniques, which are prone to failure when only archival formalin fixed paraffin embedded (FFPE) tissue is available. The purpose of this pilot study was to evaluate the utility of NanoString probe technology for the detection of MYB-NFIB transcripts in archival ACC tissue.

METHODS

A NanoString probeset panel was designed targeting the junctions of three currently annotated MYB-NFIB fusion genes as well as 5'/3' MYB probesets designed to detect MYB gene expression imbalance. RNA isolated from twenty-five archival ACC specimens was profiled and analyzed. RT-qPCR and sequencing were performed to confirm NanoString results. MYB protein expression was analyzed by immunohistochemistry.

RESULTS

Of the 25 samples analyzed, 11/25 (44%) expressed a high degree of MYB 5'/3' imbalance and five of these samples were positive for at least one specific MYB-NFIB variant in our panel. MYB-NFIB variant detection on NanoString analysis was confirmed by direct cDNA sequencing. No clinical correlations were found to be associated with MYB fusion status.

CONCLUSION

We conclude that the application of NanoString digital probe counting technology is well suited for the detection and quantification of MYB-NFIB fusion transcripts in archival ACC specimens.

Robustness of RNA sequencing on older formalin-fixed paraffin-embedded tissue from high-grade ovarian serous adenocarcinomas

Formalin-fixed paraffin-embedded (FFPE) tissues are among the most widely available clinical specimens. Their potential utility as a source of RNA for transcriptome studies would greatly enhance population-based cancer studies. Although preliminary studies suggest FFPE tissue may be used for RNA sequencing, the effect of storage time on these specimens needs to be determined. We conducted this study to determine whether RNA in archived FFPE high-grade ovarian serous adenocarcinomas from Surveillance, Epidemiology and End Results (SEER) registries was present in sufficient quantity and quality for RNA-Seq analysis. FFPE tissues, stored from 7 to 32 years, were obtained from three SEER sites. RNA was extracted, quantified, quality assessed, and subjected to RNA-Seq (a whole transcriptome sequencing technology). FFPE specimens stored for longer periods of time had poorer RNA sample quality as indicated by negative correlations between specimen storage time and fragment distribution values (DV). In addition, sample contamination was a common issue among the RNA, with 41 of 67 samples having 5% to 48% bacterial contamination. However, regardless of specimen storage time and bacterial contamination, 60% of the samples yielded data that enabled gene expression quantification, identifying more than 10,000 genes, with the correlations among most biological replicates above 0.7. This study demonstrates that FFPE high-grade ovarian serous adenocarcinomas specimens stored in repositories for up to 32 years and under varying storage conditions are a promising source of RNA for RNA-Seq. We also describe certain caveats to be considered when designing RNA-Seq studies using archived FFPE tissues.

Deleterious effects of formalin-fixation and delays to fixation on RNA and miRNA-Seq profiles

The National Cancer Institute conducted the Biospecimen Pre-analytical Variables (BPV) study to determine the effects of formalin fixation and delay to fixation (DTF) on the analysis of nucleic acids. By performing whole transcriptome sequencing and small RNA profiling on matched snap-frozen and FFPE specimens exposed to different delays to fixation, this study aimed to determine acceptable delays to fixation and proper workflow for accurate and reliable Next-Generation Sequencing (NGS) analysis of FFPE specimens. In comparison to snap-freezing, formalin fixation changed the relative proportions of intronic/exonic/untranslated RNA captured by RNA-seq for most genes. The effects of DTF on NGS analysis were negligible. In 80% of specimens, a subset of RNAs was found to differ between snap-frozen and FFPE specimens in a consistent manner across tissue groups; this subset was unaffected in the remaining 20% of specimens. In contrast, miRNA expression was generally stable across various formalin fixation protocols, but displayed increased variability following a 12 h delay to fixation.

Chromoanasynthesis is a common mechanism that leads to ERBB2 amplifications in a cohort of early stage HER2+ breast cancer samples

BACKGROUND

HER2 positive (HER2+) breast cancers involve chromosomal structural alterations that act as oncogenic driver events.

METHODS

We interrogated the genomic structure of 18 clinically-defined HER2+ breast tumors through integrated analysis of whole genome and transcriptome sequencing, coupled with clinical information.

RESULTS

ERBB2 overexpression in 15 of these tumors was associated with ERBB2 amplification due to chromoanasynthesis with six of them containing single events and the other nine exhibiting multiple events. Two of the more complex cases had adverse clinical outcomes. Chromosomes 8 was commonly involved in the same chromoanasynthesis with 17. In ten cases where chromosome 8 was involved we observed NRG1 fusions (two cases), NRG1 amplification (one case), FGFR1 amplification and ADAM32 or ADAM5 fusions. ERBB3 over-expression was associated with NRG1 fusions and EGFR and ERBB3 expressions were anti-correlated. Of the remaining three cases, one had a small duplication fully encompassing ERBB2 and was accompanied with a pathogenic mutation.

CONCLUSION

Chromoanasynthesis involving chromosome 17 can lead to ERBB2 amplifications in HER2+ breast cancer. However, additional large genomic alterations contribute to a high level of genomic complexity, generating the hypothesis that worse outcome could be associated with multiple chromoanasynthetic events.

Systematic evaluation of RNA quality, microarray data reliability and pathway analysis in fresh, fresh frozen and formalin-fixed paraffin-embedded tissue samples

Formalin-fixed paraffin-embedded (FFPE) tissues are valuable resources commonly used in pathology. However, formalin fixation modifies nucleic acids challenging the isolation of high-quality RNA for genetic profiling. Here, we assessed feasibility and reliability of microarray studies analysing transcriptome data from fresh, fresh-frozen (FF) and FFPE tissues. We show that reproducible microarray data can be generated from only 2 ng FFPE-derived RNA. For RNA quality assessment, fragment size distribution (DV200) and qPCR proved most suitable. During RNA isolation, extending tissue lysis time to 10 hours reduced high-molecular-weight species, while additional incubation at 70 °C markedly increased RNA yields. Since FF- and FFPE-derived microarrays constitute different data entities, we used indirect measures to investigate gene signal variation and relative gene expression. Whole-genome analyses revealed high concordance rates, while reviewing on single-genes basis showed higher data variation in FFPE than FF arrays. Using an experimental model, gene set enrichment analysis (GSEA) of FFPE-derived microarrays and fresh tissue-derived RNA-Seq datasets yielded similarly affected pathways confirming the applicability of FFPE tissue in global gene expression analysis. Our study provides a workflow comprising RNA isolation, quality assessment and microarray profiling using minimal RNA input, thus enabling hypothesis-generating pathway analyses from limited amounts of precious, pathologically significant FFPE tissues.

Robust detection of immune transcripts in FFPE samples using targeted RNA sequencing

Current criteria for identifying cancer patients suitable for immunotherapy with immune checkpoint blockers (ICBs) are subjective and prone to misinterpretation, as they mainly rely on the visual assessment of CD274 (best known as PD-L1) expression levels by immunohistochemistry (IHC). To address this issue, we developed a RNA sequencing (RNAseq)-based approach that specifically measures the abundance of immune transcripts in formalin-fixed paraffin embedded (FFPE) specimens. Besides exhibiting superior sensitivity as compared to whole transcriptome RNAseq, our assay requires little starting material, implying that it is compatible with RNA degradation normally caused by formalin. Here, we demonstrate that a targeted RNAseq panel reliably profiles mRNA expression levels in FFPE samples from a cohort of ovarian carcinoma patients. The expression profile of immune transcripts as measured by targeted RNAseq in FFPE versus freshly frozen (FF) samples from the same tumor was highly concordant, in spite of the RNA quality issues associated with formalin fixation. Moreover, the results of targeted RNAseq on FFPE specimens exhibited a robust correlation with mRNA expression levels as measured on the same samples by quantitative RT-PCR, as well as with protein abundance as determined by IHC. These findings demonstrate that RNAseq profiling on archival FFPE tissues can be used reliably in studies assessing the efficacy of cancer immunotherapy.

Robust transcriptional tumor signatures applicable to both formalin-fixed paraffin-embedded and fresh-frozen samples

Formalin-fixed paraffin-embedded (FFPE) samples represent a valuable resource for clinical researches. However, FFPE samples are usually considered an unreliable source for gene expression analysis due to the partial RNA degradation. In this study, through comparing gene expression profiles between FFPE samples and paired fresh-frozen (FF) samples for three cancer types, we firstly showed that expression measurements of thousands of genes had at least two-fold change in FFPE samples compared with paired FF samples. Therefore, for a transcriptional signature based on risk scores summarized from the expression levels of the signature genes, the risk score thresholds trained from FFPE (or FF) samples could not be applied to FF (or FFPE) samples. On the other hand, we found that more than 90% of the relative expression orderings (REOs) of gene pairs in the FF samples were maintained in their paired FFPE samples and largely unaffected by the storage time. The result suggested that the REOs of gene pairs were highly robust against partial RNA degradation in FFPE samples. Finally, as a case study, we developed a REOs-based signature to distinguish liver cirrhosis from hepatocellular carcinoma (HCC) using FFPE samples. The signature was validated in four datasets of FFPE samples and eight datasets of FF samples. In conclusion, the valuable FFPE samples can be fully exploited to identify REOs-based diagnostic and prognostic signatures which could be robustly applicable to both FF samples and FFPE samples with degraded RNA.

Accurate RNA Sequencing From Formalin-Fixed Cancer Tissue To Represent High-Quality Transcriptome From Frozen Tissue

PURPOSE

Accurate transcriptional sequencing (RNA-seq) from formalin-fixation and paraffin-embedding (FFPE) tumor samples presents an important challenge for translational research and diagnostic development. In addition, there are now several different protocols to prepare a sequencing library from total RNA. We evaluated the accuracy of RNA-seq data generated from FFPE samples in terms of expression profiling.

METHODS

We designed a biospecimen study to directly compare gene expression results from different protocols to prepare libraries for RNA-seq from human breast cancer tissues, with randomization to fresh-frozen (FF) or FFPE conditions. The protocols were compared using multiple computational methods to assess alignment of reads to reference genome, and the uniformity and continuity of coverage; as well as the variance and correlation, of overall gene expression and patterns of measuring coding sequence, phenotypic patterns of gene expression, and measurements from representative multigene signatures.

RESULTS

The principal determinant of variance in gene expression was use of exon capture probes, followed by the conditions of preservation (FF versus FFPE), and phenotypic differences between breast cancers. One protocol, with RNase H-based rRNA depletion, exhibited least variability of gene expression measurements, strongest correlation between FF and FFPE samples, and was generally representative of the transcriptome from standard FF RNA-seq protocols.

CONCLUSION

Method of RNA-seq library preparation from FFPE samples had marked effect on the accuracy of gene expression measurement compared to matched FF samples. Nevertheless, some protocols produced highly concordant expression data from FFPE RNA-seq data, compared to RNA-seq results from matched frozen samples.

Identification of CHEK1, SLC26A4, c-KIT, TPO and TG as new biomarkers for human follicular thyroid carcinoma

The search for preoperative biomarkers for thyroid malignancies, in particular for follicular thyroid carcinoma (FTC) diagnostics, is of utmost clinical importance. We thus aimed at screening for potential biomarker candidates for FTC. To evaluate dynamic alterations in molecular patterns as a function of thyroid malignancy progression, a comparative analysis was conducted in clinically distinct subgroups of FTC and poorly differentiated thyroid carcinoma (PDTC) nodules. NanoString analysis of FFPE samples was performed in 22 follicular adenomas, 56 FTC and 25 PDTC nodules, including oncocytic and non-oncocytic subgroups. The expression levels of CHEK1, c-KIT, SLC26A4, TG and TPO were significantly altered in all types of thyroid carcinomas. Based on collective changes of these biomarkers which correlating among each other, a predictive score has been established, allowing for discrimination between benign and FTC samples with high sensitivity and specificity. Additional transcripts related to thyroid function, cell cycle, circadian clock, and apoptosis regulation were altered in the more aggressive oncocytic subgroups only, with expression levels correlating with disease progression. Distinct molecular patterns were observed for oncocytic and non-oncocytic FTCs and PDTCs. A predictive score correlation coefficient based on collective alterations of identified here biomarkers might help to improve the preoperative diagnosis of FTC nodules.

Cell-of-origin of diffuse large B-cell lymphomas determined by the Lymph2Cx assay: better prognostic indicator than Hans algorithm

Diffuse large B-cell lymphomas (DLBCLs) are clinically heterogeneous and need a biomarker that can predict the outcome of treatments accurately. To assess the prognostic significance of the cell-of-origin type for DLBCLs, we applied the Lymph2Cx assay using a NanoString gene expression platform on formalin-fixed paraffin wax-embedded pretreatment tissues obtained from 82 patients with de novo DLBCL, not otherwise specified. All patients were treated with rituximab plus cyclophosphamide, doxorubicin, vincristine, and prednisone (R-CHOP) as the first line of chemotherapy. Based on the expression levels of Bcl-6, CD10, and MUM-1 measured by immunohistochemistry, cases were subdivided into germinal center B-cell (GCB) and non-GCB types according to the Hans algorithm. NanoString assay was performed on 82 cases. The Lymph2Cx assay successfully classified 82 cases into three categories: activated B-cell (ABC), GCB, and unclassified types. The concordance rate between the Lymph2Cx assay and the Hans algorithm was 73.6%. The Lymph2Cx-defined ABC type had significantly poorer outcomes compared with the GCB type (5-year overall survival, GCB vs. ABC, 96.6% vs. 77.1%, P = 0.020; 5-year disease-free survival, GCB vs. ABC, 96.6% vs. 79.2%, P = 0.018). In contrast, no significant differences were observed in survival between the two patient subgroups with DLBCL types classified by the Hans algorithm. The Lymph2Cx assay is a robust, reliable method for predicting the outcome of patients with DLBCL treated with R-CHOP chemotherapy.

Genome-wide primary transcriptome analysis of H2-producing archaeon Thermococcus onnurineus NA1

In spite of their pivotal roles in transcriptional and post-transcriptional processes, the regulatory elements of archaeal genomes are not yet fully understood. Here, we determine the primary transcriptome of the H2-producing archaeon Thermococcus onnurineus NA1. We identified 1,082 purine-rich transcription initiation sites along with well-conserved TATA box, A-rich B recognition element (BRE), and promoter proximal element (PPE) motif in promoter regions, a high pyrimidine nucleotide content (T/C) at the -1 position, and Shine-Dalgarno (SD) motifs (GGDGRD) in 5' untranslated regions (5' UTRs). Along with differential transcript levels, 117 leaderless genes and 86 non-coding RNAs (ncRNAs) were identified, representing diverse cellular functions and potential regulatory functions under the different growth conditions. Interestingly, we observed low GC content in ncRNAs for RNA-based regulation via unstructured forms or interaction with other cellular components. Further comparative analysis of T. onnurineus upstream regulatory sequences with those of closely related archaeal genomes demonstrated that transcription of orthologous genes are initiated by highly conserved promoter sequences, however their upstream sequences for transcriptional and translational regulation are largely diverse. These results provide the genetic information of T. onnurineus for its future application in metabolic engineering.

The Utilization of Formalin Fixed-Paraffin-Embedded Specimens in High Throughput Genomic Studies

High throughput genomic assays empower us to study the entire human genome in short time with reasonable cost. Formalin fixed-paraffin-embedded (FFPE) tissue processing remains the most economical approach for longitudinal tissue specimen storage. Therefore, the ability to apply high throughput genomic applications to FFPE specimens can expand clinical assays and discovery. Many studies have measured the accuracy and repeatability of data generated from FFPE specimens using high throughput genomic assays. Together, these studies demonstrate feasibility and provide crucial guidance for future studies using FFPE specimens. Here, we summarize the findings of these studies and discuss the limitations of high throughput data generated from FFPE specimens across several platforms that include microarray, high throughput sequencing, and NanoString.

A Comparison of RNA-Seq Results from Paired Formalin-Fixed Paraffin-Embedded and Fresh-Frozen Glioblastoma Tissue Samples

The molecular classification of glioblastoma (GBM) based on gene expression might better explain outcome and response to treatment than clinical factors. Whole transcriptome sequencing using next-generation sequencing platforms is rapidly becoming accepted as a tool for measuring gene expression for both research and clinical use. Fresh frozen (FF) tissue specimens of GBM are difficult to obtain since tumor tissue obtained at surgery is often scarce and necrotic and diagnosis is prioritized over freezing. After diagnosis, leftover tissue is usually stored as formalin-fixed paraffin-embedded (FFPE) tissue. However, RNA from FFPE tissues is usually degraded, which could hamper gene expression analysis. We compared RNA-Seq data obtained from matched pairs of FF and FFPE GBM specimens. Only three FFPE out of eleven FFPE-FF matched samples yielded informative results. Several quality-control measurements showed that RNA from FFPE samples was highly degraded but maintained transcriptomic similarities to RNA from FF samples. Certain issues regarding mutation analysis and subtype prediction were detected. Nevertheless, our results suggest that RNA-Seq of FFPE GBM specimens provides reliable gene expression data that can be used in molecular studies of GBM if the RNA is sufficiently preserved.

Recurrent hormone-binding domain truncated ESR1 amplifications in primary endometrial cancers suggest their implication in hormone independent growth

The estrogen receptor alpha (ERα) is highly expressed in both endometrial and breast cancers, and represents the most prevalent therapeutic target in breast cancer. However, anti-estrogen therapy has not been shown to be effective in endometrial cancer. Recently it has been shown that hormone-binding domain alterations of ERα in breast cancer contribute to acquired resistance to anti-estrogen therapy. In analyses of genomic data from The Cancer Genome Atlas (TCGA), we observe that endometrial carcinomas manifest recurrent ESR1 gene amplifications that truncate the hormone-binding domain encoding region of ESR1 and are associated with reduced mRNA expression of exons encoding the hormone-binding domain. These findings support a role for hormone-binding alterations of ERα in primary endometrial cancer, with potentially important therapeutic implications.

Assessment of Tumor Heterogeneity, as Evidenced by Gene Expression Profiles, Pathway Activation, and Gene Copy Number, in Patients with Multifocal Invasive Lobular Breast Tumors

BACKGROUND

Invasive lobular carcinoma (ILC) comprises approximately ~10-20% of breast cancers. In general, multifocal/multicentric (MF/MC) breast cancer has been associated with an increased rate of regional lymph node metastases. Tumor heterogeneity between foci represents a largely unstudied source of genomic variation in those rare patients with MF/MC ILC.

METHODS

We characterized gene expression and copy number in 2 or more foci from 11 patients with MF/MC ILC (all ER+, HER2-) and adjacent normal tissue. RNA and DNA were extracted from 3x1.5 mm cores from all foci. Gene expression (730 genes) and copy number (80 genes) were measured using Nanostring PanCancer and Cancer CNV panels. Linear mixed models were employed to compare expression in tumor versus normal samples from the same patient, and to assess heterogeneity (variability) in expression among multiple ILC within an individual.

RESULTS

35 and 34 genes were upregulated (FC>2) and down-regulated (FC<0.5) respectively in ILC tumor relative to adjacent normal tissue, q<0.05. 9/34 down-regulated genes (FIGF, RELN, PROM1, SFRP1, MMP7, NTRK2, LAMB3, SPRY2, KIT) had changes larger than CDH1, a hallmark of ILC. Copy number changes in these patients were relatively few but consistent across foci within each patient. Amplification of three genes (CCND1, FADD, ORAOV1) at 11q13.3 was present in 2/11 patients in both foci. We observed significant evidence of within-patient between-foci variability (heterogeneity) in gene expression for 466 genes (p<0.05 with FDR 8%), including CDH1, FIGF, RELN, SFRP1, MMP7, NTRK2, LAMB3, SPRY2 and KIT.

CONCLUSIONS

There was substantial variation in gene expression between ILC foci within patients, including known markers of ILC, suggesting an additional level of complexity that should be addressed.

Cross-platform normalization of microarray and RNA-seq data for machine learning applications

Large, publicly available gene expression datasets are often analyzed with the aid of machine learning algorithms. Although RNA-seq is increasingly the technology of choice, a wealth of expression data already exist in the form of microarray data. If machine learning models built from legacy data can be applied to RNA-seq data, larger, more diverse training datasets can be created and validation can be performed on newly generated data. We developed Training Distribution Matching (TDM), which transforms RNA-seq data for use with models constructed from legacy platforms. We evaluated TDM, as well as quantile normalization, nonparanormal transformation, and a simple log 2 transformation, on both simulated and biological datasets of gene expression. Our evaluation included both supervised and unsupervised machine learning approaches. We found that TDM exhibited consistently strong performance across settings and that quantile normalization also performed well in many circumstances. We also provide a TDM package for the R programming language.

Evaluating Robustness and Sensitivity of the NanoString Technologies nCounter Platform to Enable Multiplexed Gene Expression Analysis of Clinical Samples

Analysis of clinical trial specimens such as formalin-fixed paraffin-embedded (FFPE) tissue for molecular mechanisms of disease progression or drug response is often challenging and limited to a few markers at a time. This has led to the increasing importance of highly multiplexed assays that enable profiling of many biomarkers within a single assay. Methods for gene expression analysis have undergone major advances in biomedical research, but obtaining a robust dataset from low-quality RNA samples, such as those isolated from FFPE tissue, remains a challenge. Here, we provide a detailed evaluation of the NanoString Technologies nCounter platform, which provides a direct digital readout of up to 800 mRNA targets simultaneously. We tested this system by examining a broad set of human clinical tissues for a range of technical variables, including sensitivity and limit of detection to varying RNA quantity and quality, reagent performance over time, variability between instruments, the impact of the number of fields of view sampled, and differences between probe sequence locations and overlapping genes across CodeSets. This study demonstrates that Nanostring offers several key advantages, including sensitivity, reproducibility, technical robustness, and utility for clinical application.

RNA-seq-based transcriptome profiling reveals differential gene expression in the lungs of Sprague-Dawley rats during early-phase acute hypobaric hypoxia

Individuals subjected to hypobaric hypoxia at high altitudes may exhibit differential physiological responses in terms of susceptibility and tolerance to the development of hypoxia-related disorders. We studied early-phase gene expression in the lungs of Sprague-Dawley rats exhibiting such differential physiological responses after exposure to acute hypobaric hypoxia for 1 h at a simulated altitude of 9144 m. RNA-seq transcriptome profiling of lung tissues revealed differential gene expression in tolerant and susceptible groups, subsequently validated by qRT-PCR for ten selected differentially expressed genes. The gene expression pattern indicated hypometabolism and negative regulation of vasoconstriction in all groups except susceptible rats, coupled with altered MAPK, p53 and JAK-STAT signaling. Upregulation of early-phase response genes including Dusp1 (dual specificity phosphatase), Cdkn1a (cyclin-dependent kinase inhibitor 1A), Txnip (thioredoxin-interacting protein), Rgs1 (regulator of G-protein signaling 1) and Rgs2 (regulator of G-protein signaling 2) in susceptible rats indicated a progression toward growth arrest and apoptosis. Enhanced expression of cell adhesion molecules, wound healing and repair bioprocesses was observed in tolerant males. Upregulated Kcnj15 (potassium inwardly rectifying channel subfamily j membrane 15) and Vsig4 (V-set and Ig domain containing 4) variants in tolerant females suggested adaptation to hypoxia possibly by fluid reabsorption to avoid edematous conditions and suppression of T cell proliferation to avoid acute lung inflammation. Our study might help in understanding the molecular-physiological mechanisms associated with progressive damage in the lung tissues of susceptible and tissue-protective measures in tolerant rats during acute hypobaric hypoxia.

Identification of new biomarkers for human papillary thyroid carcinoma employing NanoString analysis

We previously reported an upregulation of the clock transcript BMAL1, correlating with TIMP1 expression in fresh-frozen samples from papillary thyroid carcinoma (PTC). Since frozen postoperative biopsy samples are difficult to obtain, we aimed to validate the application of high-precision NanoString analysis for formalin-fixed paraffin-embedded (FFPE) thyroid nodule samples and to screen for potential biomarkers associated with PTC. No significant differences were detected between fresh-frozen and FFPE samples. NanoString analysis of 51 transcripts in 17 PTC and 17 benign nodule samples obtained from different donors and in 24 pairs of benign and PTC nodules, obtained from the same donor (multinodular goiters), confirmed significant alterations in the levels of BMAL1, c-MET, c-KIT, TIMP1, and other transcripts. Moreover, we identified for the first time alterations in CHEK1 and BCL2 levels in PTC. A predictive score was established for each sample, based on the combined expression levels of BMAL1, CHEK1, c-MET, c-KIT and TIMP1. In combination with BRAF mutation analysis, this predictive score closely correlated with the clinicopathological characteristics of the analyzed thyroid nodules. Our study identified new thyroid transcripts with altered levels in PTC using the NanoString approach. A predictive score correlation coefficient might contribute to improve the preoperative diagnosis of thyroid nodules.

Comparing Platforms for Messenger RNA Expression Profiling of Archival Formalin-Fixed, Paraffin-Embedded Tissues

Archival formalin-fixed, paraffin-embedded (FFPE) tissue specimens represent a readily available but largely untapped resource for gene expression profiling-based biomarker discovery. Several technologies have been proposed to cope with the bias from RNA cross-linking and degradation associated with archival specimens to generate data comparable with RNA from fresh-frozen materials. Direct comparison studies of these RNA expression platforms remain rare. We compared two commercially available platforms for RNA expression profiling of archival FFPE specimens from clinical studies of prostate and ovarian cancer: the Affymetrix Human Gene 1.0ST Array following whole-transcriptome amplification using the NuGen WT-Ovation FFPE System V2, and the NanoString nCounter without amplification. For each assay, we profiled 7 prostate and 11 ovarian cancer specimens, with a block age of 4 to 21 years. Both platforms produced gene expression profiles with high sensitivity and reproducibility through technical repeats from FFPE materials. Sensitivity and reproducibility remained high across block age within each cohort. A strong concordance was shown for the transcript expression values for genes detected by both platforms. We showed the biological validity of specific gene signatures generated by both platforms for both cohorts. Our study supports the feasibility of gene expression profiling and large-scale signature validation on archival prostate and ovarian tumor specimens using commercial platforms. These approaches have the potential to aid precision medicine with biomarker discovery and validation.

Das methodische Potenzial der neuen Sequenziertechnologien jenseits der Mutationssuche

In diesem Beitrag wird eine Reihe wichtiger Anwendungen der neuen Sequenziertechnologien bzw. des Next Generation Sequencing (NGS) vorgestellt. An ausgewählten Beispielen werden für jede Methode die Anwendungsmöglichkeiten in der humangenetischen Forschung dargestellt, jeweils das prinzipielle Vorgehen beschrieben und mögliche Quellen für ausführliche Arbeitsanweisungen vorgestellt. Die beschriebenen Techniken umfassen im Einzelnen: RNA-Sequenzierung mittels NGS („RNA-Seq“), Chromatinimmunpräzipitation in Kombination mit NGS („ChIP-Seq“), „ribosome profiling“, Immunpräzipitation methylierter DNA-Segmente in Kombination mit NGS („methylated DNA immunoprecipitation“ bzw. „MeDIP-Seq“) und die HiC-Technik, eine Weiterentwicklung der Chromosome-Conformation-Capture(3c)-Methode.

MAP-RSeq: Mayo Analysis Pipeline for RNA sequencing

Background

Although the costs of next generation sequencing technology have decreased over the past years, there is still a lack of simple-to-use applications, for a comprehensive analysis of RNA sequencing data. There is no one-stop shop for transcriptomic genomics. We have developed MAP-RSeq, a comprehensive computational workflow that can be used for obtaining genomic features from transcriptomic sequencing data, for any genome.

Results

For optimization of tools and parameters, MAP-RSeq was validated using both simulated and real datasets. MAP-RSeq workflow consists of six major modules such as alignment of reads, quality assessment of reads, gene expression assessment and exon read counting, identification of expressed single nucleotide variants (SNVs), detection of fusion transcripts, summarization of transcriptomics data and final report. This workflow is available for Human transcriptome analysis and can be easily adapted and used for other genomes. Several clinical and research projects at the Mayo Clinic have applied the MAP-RSeq workflow for RNA-Seq studies. The results from MAP-RSeq have thus far enabled clinicians and researchers to understand the transcriptomic landscape of diseases for better diagnosis and treatment of patients.

Conclusions

Our software provides gene counts, exon counts, fusion candidates, expressed single nucleotide variants, mapping statistics, visualizations, and a detailed research data report for RNA-Seq. The workflow can be executed on a standalone virtual machine or on a parallel Sun Grid Engine cluster. The software can be downloaded from http://bioinformaticstools.mayo.edu/research/maprseq/.

Comparison of RNA-Seq by poly (A) capture, ribosomal RNA depletion, and DNA microarray for expression profiling

BACKGROUND

RNA sequencing (RNA-Seq) is often used for transcriptome profiling as well as the identification of novel transcripts and alternative splicing events. Typically, RNA-Seq libraries are prepared from total RNA using poly(A) enrichment of the mRNA (mRNA-Seq) to remove ribosomal RNA (rRNA), however, this method fails to capture non-poly(A) transcripts or partially degraded mRNAs. Hence, a mRNA-Seq protocol will not be compatible for use with RNAs coming from Formalin-Fixed and Paraffin-Embedded (FFPE) samples.

RESULTS

To address the desire to perform RNA-Seq on FFPE materials, we evaluated two different library preparation protocols that could be compatible for use with small RNA fragments. We obtained paired Fresh Frozen (FF) and FFPE RNAs from multiple tumors and subjected these to different gene expression profiling methods. We tested 11 human breast tumor samples using: (a) FF RNAs by microarray, mRNA-Seq, Ribo-Zero-Seq and DSN-Seq (Duplex-Specific Nuclease) and (b) FFPE RNAs by Ribo-Zero-Seq and DSN-Seq. We also performed these different RNA-Seq protocols using 10 TCGA tumors as a validation set.The data from paired RNA samples showed high concordance in transcript quantification across all protocols and between FF and FFPE RNAs. In both FF and FFPE, Ribo-Zero-Seq removed rRNA with comparable efficiency as mRNA-Seq, and it provided an equivalent or less biased coverage on gene 3' ends. Compared to mRNA-Seq where 69% of bases were mapped to the transcriptome, DSN-Seq and Ribo-Zero-Seq contained significantly fewer reads mapping to the transcriptome (20-30%); in these RNA-Seq protocols, many if not most reads mapped to intronic regions. Approximately 14 million reads in mRNA-Seq and 45-65 million reads in Ribo-Zero-Seq or DSN-Seq were required to achieve the same gene detection levels as a standard Agilent DNA microarray.

CONCLUSIONS

Our results demonstrate that compared to mRNA-Seq and microarrays, Ribo-Zero-Seq provides equivalent rRNA removal efficiency, coverage uniformity, genome-based mapped reads, and consistently high quality quantification of transcripts. Moreover, Ribo-Zero-Seq and DSN-Seq have consistent transcript quantification using FFPE RNAs, suggesting that RNA-Seq can be used with FFPE-derived RNAs for gene expression profiling.

Next-generation sequencing of RNA and DNA isolated from paired fresh-frozen and formalin-fixed paraffin-embedded samples of human cancer and normal tissue

Formalin-fixed, paraffin-embedded (FFPE) tissues are an invaluable resource for clinical research. However, nucleic acids extracted from FFPE tissues are fragmented and chemically modified making them challenging to use in molecular studies. We analysed 23 fresh-frozen (FF), 35 FFPE and 38 paired FF/FFPE specimens, representing six different human tissue types (bladder, prostate and colon carcinoma; liver and colon normal tissue; reactive tonsil) in order to examine the potential use of FFPE samples in next-generation sequencing (NGS) based retrospective and prospective clinical studies. Two methods for DNA and three methods for RNA extraction from FFPE tissues were compared and were found to affect nucleic acid quantity and quality. DNA and RNA from selected FFPE and paired FF/FFPE specimens were used for exome and transcriptome analysis. Preparations of DNA Exome-Seq libraries was more challenging (29.5% success) than that of RNA-Seq libraries, presumably because of modifications to FFPE tissue-derived DNA. Libraries could still be prepared from RNA isolated from two-decade old FFPE tissues. Data were analysed using the CLC Bio Genomics Workbench and revealed systematic differences between FF and FFPE tissue-derived nucleic acid libraries. In spite of this, pairwise analysis of DNA Exome-Seq data showed concordance for 70–80% of variants in FF and FFPE samples stored for fewer than three years. RNA-Seq data showed high correlation of expression profiles in FF/FFPE pairs (Pearson Correlations of 0.90 +/- 0.05), irrespective of storage time (up to 244 months) and tissue type. A common set of 1,494 genes was identified with expression profiles that were significantly different between paired FF and FFPE samples irrespective of tissue type. Our results are promising and suggest that NGS can be used to study FFPE specimens in both prospective and retrospective archive-based studies in which FF specimens are not available.