Breast cancer is a highly heterogeneous disease that can be classified into multiple subtypes based on the tumor transcriptome. Most of the subtyping schemes used in clinics today are derived from analyses of microarray data from thousands of different tumors together with clinical data for the patients from which the tumors were isolated. However, RNA sequencing is gradually replacing microarrays as the preferred transcriptomics platform, and although transcript abundances measured by the two different technologies are largely compatible, subtyping methods developed for probe-based microarray data are incompatible with RNA sequencing as input data. Here we present an RNA sequencing data processing pipeline, which relies on the mapping of sequencing reads to the probe set target sequences instead of the human reference genome, thereby enabling probe-based subtyping of breast cancer tumor tissue using sequencing-based transcriptomics. By analyzing 66 breast cancer tumors for which gene expression was measured using both microarrays and RNA sequencing, we show that RNA sequencing data can be directly compared to microarray data using our pipeline. Additionally, we demonstrate that the established subtyping method CITBCMST (Guedj et al., 2012), which relies on a 375 probe set-signature to classify samples into the six subtypes basL, lumA, lumB, lumC, mApo, and normL, can be applied without further modifications. This pipeline enables a seamless transition to sequencing-based transcriptomics for future clinical purposes.
Bibliographical noteThis is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.
- Gene expression
- RNA sequencing
- Breast cancer
- Molecular subtyping