Building flexible and robust analysis frameworks for molecular subtyping of cancers.

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

Dasatinib overcomes glucocorticoid resistance in B-cell acute lymphoblastic leukemia.

Resistance to glucocorticoids (GC) is associated with an increased risk of relapse in B-cell progenitor acute lymphoblastic leukemia (BCP-ALL). Performing transcriptomic and single-cell proteomic studies in healthy B-cell progenitors, we herein identify coordination between the glucocorticoid receptor pathway with B-cell developmental pathways. Healthy pro-B cells most highly express the glucocorticoid receptor, and this developmental expression is conserved in primary BCP-ALL cells from patients at diagnosis and relapse. In-vitro and in vivo glucocorticoid treatment of primary BCP-ALL cells demonstrate that the interplay between B-cell development and the glucocorticoid pathways is crucial for GC resistance in leukemic cells. Gene set enrichment analysis in BCP-ALL cell lines surviving GC treatment show enrichment of B cell receptor signaling pathways. In addition, primary BCP-ALL cells surviving GC treatment in vitro and in vivo demonstrate a late pre-B cell phenotype with activation of PI3K/mTOR and CREB signaling. Dasatinib, a multi-kinase inhibitor, most effectively targets this active signaling in GC-resistant cells, and when combined with glucocorticoids, results in increased cell death in vitro and decreased leukemic burden and prolonged survival in an in vivo xenograft model. Targeting the active signaling through the addition of dasatinib may represent a therapeutic approach to overcome GC resistance in BCP-ALL.

cyCombine allows for robust integration of single-cell cytometry datasets within and across technologies.

Combining single-cell cytometry datasets increases the analytical flexibility and the statistical power of data analyses. However, in many cases the full potential of co-analyses is not reached due to technical variance between data from different experimental batches. Here, we present cyCombine, a method to robustly integrate cytometry data from different batches, experiments, or even different experimental techniques, such as CITE-seq, flow cytometry, and mass cytometry. We demonstrate that cyCombine maintains the biological variance and the structure of the data, while minimizing the technical variance between datasets. cyCombine does not require technical replicates across datasets, and computation time scales linearly with the number of cells, allowing for integration of massive datasets. Robust, accurate, and scalable integration of cytometry data enables integration of multiple datasets for primary data analyses and the validation of results using public datasets.

A Comparison of Tools for Copy-Number Variation Detection in Germline Whole Exome and Whole Genome Sequencing Data.

Copy-number variations (CNVs) have important clinical implications for several diseases and cancers. Relevant CNVs are hard to detect because common structural variations define large parts of the human genome. CNV calling from short-read sequencing would allow single protocol full genomic profiling. We reviewed 50 popular CNV calling tools and included 11 tools for benchmarking in a reference cohort encompassing 39 whole genome sequencing (WGS) samples paired current clinical standard-SNP-array based CNV calling. Additionally, for nine samples we also performed whole exome sequencing (WES), to address the effect of sequencing protocol on CNV calling. Furthermore, we included Gold Standard reference sample NA12878, and tested 12 samples with CNVs confirmed by multiplex ligation-dependent probe amplification (MLPA). Tool performance varied greatly in the number of called CNVs and bias for CNV lengths. Some tools had near-perfect recall of CNVs from arrays for some samples, but poor precision. Several tools had better performance for NA12878, which could be a result of overfitting. We suggest combining the best tools also based on different methodologies: GATK gCNV, Lumpy, DELLY, and cn.MOPS. Reducing the total number of called variants could potentially be assisted by the use of background panels for filtering of frequently called variants.

Clinical implications of intrinsic molecular subtypes of breast cancer for sentinel node status.

Axillary lymph node status is an important prognostic factor for breast cancer patients and sentinel lymph node biopsy (SLNB) is a less invasive surgical proxy. We examined if consecutively derived molecular subtypes from primary breast cancers provide additional predictive value for SLNB status. 1556 patients with a breast cancer > 10 mm underwent primary surgical procedure including SLNB and tumor specimens were assigned with a transcriptomics-based molecular subtype. 1020 patients had a negative sentinel node (SN) and 536 a positive. A significant association between tumor size and SN status (p < 0.0001) was found across all samples, but no association between size and SN status (p = 0.14) was found for BasL tumors. A BasL subtype was a predictor of an SN-negative status (p = 0.001, OR 0.58, 95% CI 0.38;0.90) and among the BasL, postmenopausal status was a predictor for SN-negative status (p = 0.01). Overall survival was significantly lower (p = 0.02) in patients with BasL tumors and a positive SN. Interestingly, we identified a significant correlation between hormone receptor activity and SN status within the BasL subtype. Taken together, molecular subtypes and hormone receptor activity of breast cancers add predictive value for SLNB status.

Using microarray-based subtyping methods for breast cancer in the era of high-throughput RNA sequencing.

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 (RNA-Seq) 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-Seq as input data. Here, we present an RNA-Seq 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-Seq, we show that RNA-Seq data can be directly compared to microarray data using our pipeline. Additionally, we demonstrate that the established subtyping method CITBCMST (Guedj et al., ), 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.