Mattias Rantalainen

Reprogramming Tumor-Associated Macrophages by Antibody Targeting Inhibits Cancer Progression and Metastasis

Tumors are composed of multiple cell types besides the tumor cells themselves, including innate immune cells such as macrophages. Tumor-associated macrophages (TAMs) are a heterogeneous population of myeloid cells present in the tumor microenvironment (TME). Here, they contribute to immunosuppression, enabling the establishment and persistence of solid tumors as well as metastatic dissemination. We have found that the pattern recognition scavenger receptor MARCO defines a subtype of suppressive TAMs and is linked to clinical outcome. An anti-MARCO monoclonal antibody was developed, which induces anti-tumor activity in breast and colon carcinoma, as well as in melanoma models through reprogramming TAM populations to a pro-inflammatory phenotype and increasing tumor immunogenicity. This anti-tumor activity is dependent on the inhibitory Fc-receptor, FcγRIIB, and also enhances the efficacy of checkpoint therapy. These results demonstrate that immunotherapies using antibodies designed to modify myeloid cells of the TME represent a promising mode of cancer treatment.

Beta-Poisson model for single-cell RNA-seq data analyses

MOTIVATION Single-cell RNA-sequencing technology allows detection of gene expression at the single-cell level. One typical feature of the data is a bimodality in the cellular distribution even for highly expressed genes, primarily caused by a proportion of non-expressing cells. The standard and the over-dispersed gamma-Poisson models that are commonly used in bulk-cell RNA-sequencing are not able to capture this property. RESULTS We introduce a beta-Poisson mixture model that can capture the bimodality of the single-cell gene expression distribution. We further integrate the model into the generalized linear model framework in order to perform differential expression analyses. The whole analytical procedure is called BPSC. The results from several real single-cell RNA-seq datasets indicate that ∼90% of the transcripts are well characterized by the beta-Poisson model; the model-fit from BPSC is better than the fit of the standard gamma-Poisson model in > 80% of the transcripts. Moreover, in differential expression analyses of simulated and real datasets, BPSC performs well against edgeR, a conventional method widely used in bulk-cell RNA-sequencing data, and against scde and MAST, two recent methods specifically designed for single-cell RNA-seq data. AVAILABILITY AND IMPLEMENTATION An R package BPSC for model fitting and differential expression analyses of single-cell RNA-seq data is available under GPL-3 license at https://github.com/nghiavtr/BPSC CONTACT: yudi.pawitan@ki.se or mattias.rantalainen@ki.se SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Digital image analysis outperforms manual biomarker assessment in breast cancer

In the spectrum of breast cancers, categorization according to the four gene expression-based subtypes ‘Luminal A,’ ‘Luminal B,’ ‘HER2-enriched,’ and ‘Basal-like’ is the method of choice for prognostic and predictive value. As gene expression assays are not yet universally available, routine immunohistochemical stains act as surrogate markers for these subtypes. Thus, congruence of surrogate markers and gene expression tests is of utmost importance. In this study, 3 cohorts of primary breast cancer specimens (total n=436) with up to 28 years of survival data were scored for Ki67, ER, PR, and HER2 status manually and by digital image analysis (DIA). The results were then compared for sensitivity and specificity for the Luminal B subtype, concordance to PAM50 assays in subtype classification and prognostic power. The DIA system used was the Visiopharm Integrator System. DIA outperformed manual scoring in terms of sensitivity and specificity for the Luminal B subtype, widely considered the most challenging distinction in surrogate subclassification, and produced slightly better concordance and Cohen’s κ agreement with PAM50 gene expression assays. Manual biomarker scores and DIA essentially matched each other for Cox regression hazard ratios for all-cause mortality. When the Nottingham combined histologic grade (Elston–Ellis) was used as a prognostic surrogate, stronger Spearman’s rank-order correlations were produced by DIA. Prognostic value of Ki67 scores in terms of likelihood ratio χ2 (LR χ2) was higher for DIA that also added significantly more prognostic information to the manual scores (LR−Δχ2). In conclusion, the system for DIA evaluated here was in most aspects a superior alternative to manual biomarker scoring. It also has the potential to reduce time consumption for pathologists, as many of the steps in the workflow are either automatic or feasible to manage without pathological expertise.

Validation of risk stratification models in acute myeloid leukemia using sequencing-based molecular profiling.

Risk stratification of acute myeloid leukemia (AML) patients needs improvement. Several AML risk classification models based on somatic mutations or gene-expression profiling have been proposed. However, systematic and independent validation of these models is required for future clinical implementation. We performed whole-transcriptome RNA-sequencing and panel-based deep DNA sequencing of 23 genes in 274 intensively treated AML patients (Clinseq-AML). We also utilized the The Cancer Genome Atlas (TCGA)-AML study (N=142) as a second validation cohort. We evaluated six previously proposed molecular-based models for AML risk stratification and two revised risk classification systems combining molecular- and clinical data. Risk groups stratified by five out of six models showed different overall survival in cytogenetic normal-AML patients in the Clinseq-AML cohort (P-value<0.05; concordance index >0.5). Risk classification systems integrating mutational or gene-expression data were found to add prognostic value to the current European Leukemia Net (ELN) risk classification. The prognostic value varied between models and across cohorts, highlighting the importance of independent validation to establish evidence of efficacy and general applicability. All but one model replicated in the Clinseq-AML cohort, indicating the potential for molecular-based AML risk models. Risk classification based on a combination of molecular and clinical data holds promise for improved AML patient stratification in the future.

Sequencing-based breast cancer diagnostics as an alternative to routine biomarkers

Sequencing-based breast cancer diagnostics have the potential to replace routine biomarkers and provide molecular characterization that enable personalized precision medicine. Here we investigate the concordance between sequencing-based and routine diagnostic biomarkers and to what extent tumor sequencing contributes clinically actionable information. We applied DNA- and RNA-sequencing to characterize tumors from 307 breast cancer patients with replication in up to 739 patients. We developed models to predict status of routine biomarkers (ER, HER2,Ki-67, histological grade) from sequencing data. Non-routine biomarkers, including mutations in BRCA1, BRCA2 and ERBB2(HER2), and additional clinically actionable somatic alterations were also investigated. Concordance with routine diagnostic biomarkers was high for ER status (AUC = 0.95;AUC(replication) = 0.97) and HER2 status (AUC = 0.97;AUC(replication) = 0.92). The transcriptomic grade model enabled classification of histological grade 1 and histological grade 3 tumors with high accuracy (AUC = 0.98;AUC(replication) = 0.94). Clinically actionable mutations in BRCA1, BRCA2 and ERBB2(HER2) were detected in 5.5% of patients, while 53% had genomic alterations matching ongoing or concluded breast cancer studies. Sequencing-based molecular profiling can be applied as an alternative to histopathology to determine ER and HER2 status, in addition to providing improved tumor grading and clinically actionable mutations and molecular subtypes. Our results suggest that sequencing-based breast cancer diagnostics in a near future can replace routine biomarkers.

Determining breast cancer histological grade from RNA-sequencing data.

BackgroundThe histologic grade (HG) of breast cancer is an established prognostic factor. The grade is usually reported on a scale ranging from 1 to 3, where grade 3 tumours are the most aggressive. However, grade 2 is associated with an intermediate risk of recurrence, and carries limited information for clinical decision-making. Patients classified as grade 2 are at risk of both under- and over-treatment.MethodsRNA-sequencing analysis was conducted in a cohort of 275 women diagnosed with invasive breast cancer. Multivariate prediction models were developed to classify tumours into high and low transcriptomic grade (TG) based on gene- and isoform-level expression data from RNA-sequencing. HG2 tumours were reclassified according to the prediction model and a recurrence-free survival analysis was performed by the multivariate Cox proportional hazards regression model to assess to what extent the TG model could be used to stratify patients. The prediction model was validated in N=487 breast cancer cases from the The Cancer Genome Atlas (TCGA) data set. Differentially expressed genes and isoforms associated with HGs were analysed using linear models.ResultsThe classification of grade 1 and grade 3 tumours based on RNA-sequencing data achieved high accuracy (area under the receiver operating characteristic curve = 0.97). The association between recurrence-free survival rate and HGs was confirmed in the study population (hazard ratio of grade 3 versus 1 was 2.62 with 95 % confidence interval = 1.04–6.61). The TG model enabled us to reclassify grade 2 tumours as high TG and low TG gene or isoform grade. The risk of recurrence in the high TG group of grade 2 tumours was higher than in low TG group (hazard ratio = 2.43, 95 % confidence interval = 1.13–5.20). We found 8200 genes and 13,809 isoforms that were differentially expressed between HG1 and HG3 breast cancer tumours.ConclusionsGene- and isoform-level expression data from RNA-sequencing could be utilised to differentiate HG1 and HG3 tumours with high accuracy. We identified a large number of novel genes and isoforms associated with HG. Grade 2 tumours could be reclassified as high and low TG, which has the potential to reduce over- and under-treatment if implemented clinically.

Application of single-cell sequencing in human cancer

Abstract Precision medicine is emerging as a cornerstone of future cancer care with the objective of providing targeted therapies based on the molecular phenotype of each individual patient. Traditional bulk-level molecular phenotyping of tumours leads to significant information loss, as the molecular profile represents an average phenotype over large numbers of cells, while cancer is a disease with inherent intra-tumour heterogeneity at the cellular level caused by several factors, including clonal evolution, tissue hierarchies, rare cells and dynamic cell states. Single-cell sequencing provides means to characterize heterogeneity in a large population of cells and opens up opportunity to determine key molecular properties that influence clinical outcomes, including prognosis and probability of treatment response. Single-cell sequencing methods are now reliable enough to be used in many research laboratories, and we are starting to see applications of these technologies for characterization of human primary cancer cells. In this review, we provide an overview of studies that have applied single-cell sequencing to characterize human cancers at the single-cell level, and we discuss some of the current challenges in the field.

Molecular Differences between Screen-Detected and Interval Breast Cancers Are Largely Explained by PAM50 Subtypes

Purpose: Interval breast cancer is of clinical interest, as it exhibits an aggressive phenotype and evades detection by screening mammography. A comprehensive picture of somatic changes that drive tumors to become symptomatic in the screening interval can improve understanding of the biology underlying these aggressive tumors. Experimental Design: Initiated in April 2013, Clinical Sequencing of Cancer in Sweden (Clinseq) is a scientific and clinical platform for the genomic profiling of cancer. The breast cancer pilot study consisted of women diagnosed with breast cancer between 2001 and 2012 in the Stockholm/Gotland regions. A subset of 307 breast tumors was successfully sequenced, of which 113 were screen-detected and 60 were interval cancers. We applied targeted deep sequencing of cancer-related genes; low-pass, whole-genome sequencing; and RNA sequencing technology to characterize somatic differences in the genomic and transcriptomic architecture by interval cancer status. Mammographic density and PAM50 molecular subtypes were considered. Results: In the univariate analyses, TP53, PPP1R3A, and KMT2B were significantly more frequently mutated in interval cancers than in screen-detected cancers. Acquired somatic copy number aberrations with a frequency difference of at least 15% between the two groups included gains in 17q23-q25.3 and losses in 16q24.2. Gene expression analysis identified 447 significantly differentially expressed genes, of which 120 were replicated in an independent microarray dataset. After adjusting for PAM50, most differences were no longer significant. Conclusions: Molecular differences by interval cancer status were observed, but they were largely explained by PAM50 subtypes. This work offers new insights into the biological differences between the two tumor groups. Clin Cancer Res; 23(10); 2584–92. ©2016 AACR.

An Integrated Bioinformatics Approach for Identifying Genetic Markers that Predict Cerebrospinal Fluid Biomarker p-tau181/Aβ1-42 Ratio in ApoE4-Negative Mild Cognitive Impairment Patients

Alzheimers disease (AD) is the most common form of dementia, with no disease-modifying treatment yet available. Early detection of patients at risk of developing AD is of central importance. Blood-based genetic signatures can serve as early detection and as population-based screening tools. In this study, we aimed to identify genetic markers and gene signatures associated with cerebrospinal fluid (CSF) biomarkers levels of t-tau, p-tau181, and with the two ratios t-tau/Aβ1-42 and p-tau181/Aβ1-42 in the context of progression from mild cognitive impairment (MCI) to AD, and to identify a panel of genetic markers that can predict CSF biomarker p-tau181/Aβ1-42 ratio with consideration of APOE ε4 stratification. We analyzed genome-wide the Alzheimers Disease Neuroimaging Initiative dataset with up to 48 months follow-up. In the first part of the analysis, the main effect of single nucleotide polymorphisms (SNPs) under an additive genetic model was assessed for each of the four CSF biomarkers. In the second part of the analysis, we performed an integrated analysis of genome-wide association study results with pathway enrichment analysis, predictive modeling and network analysis in the subgroup of ApoE4-negative subjects. We identified a panel of five SNPs, rs6766238, rs1143960, rs1249963, rs11975968, and rs4836493, that are predictive for p-tau181/Aβ1-42 ratio (high/low) with a sensitivity of 66% and a specificity of 70% (AUC 0.74). These results suggest that a panel of SNPs is a potential prognostic biomarker in ApoE4-negative MCI patients.

Prognostic value of Ki67 analysed by cytology or histology in primary breast cancer

Aims The accuracy of biomarker assessment in breast pathology is vital for therapy decisions. The therapy predictive and prognostic biomarkers oestrogen receptor (ER), progesterone receptor, HER2 and Ki67 may act as surrogates to gene expression profiling of breast cancer. The aims of this study were to investigate the concordance of consecutive biomarker assessment by immunocytochemistry on preoperative fine-needle aspiration cytology versus immunohistochemistry (IHC) on the corresponding resected breast tumours. Further, to investigate the concordance with molecular subtype and correlation to stage and outcome. Methods Two retrospective cohorts comprising 385 breast tumours with clinicopathological data including gene expression-based subtype and up to 10-year overall survival data were evaluated. Results In both cohorts, we identified a substantial variation in Ki67 index between cytology and histology and a switch between low and high proliferation within the same tumour in 121/360 cases. ER evaluations were discordant in only 1.5% of the tumours. From cohort 2, gene expression data with PAM50 subtype were used to correlate surrogate subtypes. IHC-based surrogate classification could identify the correct molecular subtype in 60% and 64% of patients by cytology (n=63) and surgical resections (n=73), respectively. Furthermore, high Ki67 in surgical resections but not in cytology was associated with poor overall survival and higher probability for axillary lymph node metastasis. Conclusions This study shows considerable differences in the prognostic value of Ki67 but not ER in breast cancer depending on the diagnostic method. Furthermore, our findings show that both methods are insufficient in predicting true molecular subtypes.