Tumour immune characterisation of primary triple-negative breast cancer using automated image quantification of immunohistochemistry-stained immune cells.

The tumour immune microenvironment (TIME) in breast cancer is acknowledged with an increasing role in treatment response and prognosis. With a growing number of immune markers analysed, digital image analysis may facilitate broader TIME understanding, even in single-plex IHC data. To facilitate analyses of the latter an open-source image analysis pipeline, Tissue microarray MArker Quantification (TMArQ), was developed and applied to single-plex stainings for p53, CD3, CD4, CD8, CD20, CD68, FOXP3, and PD-L1 (SP142 antibody) in a 218-patient triple negative breast cancer (TNBC) cohort with complementary pathology scorings, clinicopathological, whole genome sequencing, and RNA-sequencing data. TMArQ's cell counts for analysed immune markers were on par with results from alternative methods and consistent with both estimates from human pathology review, different quantifications and classifications derived from RNA-sequencing as well as known prognostic patterns of immune response in TNBC. The digital cell counts demonstrated how immune markers are coexpressed in the TIME when considering TNBC molecular subtypes and DNA repair deficiency, and how combination of immune status with DNA repair deficiency status can improve the prognostic stratification in chemotherapy treated patients. These results underscore the value and potential of integrating TIME and specific tumour intrinsic alterations/phenotypes for the molecular understanding of TNBC.

FHL2 expression by cancer-associated fibroblasts promotes metastasis and angiogenesis in lung adenocarcinoma.

Cancer-associated fibroblasts (CAFs) contribute to the progression of lung cancer. Four and a half LIM domain protein-2 (FHL2) is a component of focal adhesion structures. We analyzed the function of FHL2 expressed by CAFs in lung adenocarcinoma. Expression of FHL2 in fibroblast subtypes was investigated using database of single-cell RNA-sequencing of lung cancer tissue. The role of FHL2 in the proliferation and migration of CAFs was assessed. The effects of FHL2 knockout on the migration and invasion of human lung adenocarcinoma cells and tube formation of endothelial cells induced by CAF-conditioned medium (CM) were evaluated. The effect of FHL2 knockout in CAFs on metastasis was determined using a murine orthotopic lung cancer model. The prognostic significance of stromal FHL2 was assessed by immunohistochemistry in human adenocarcinoma specimens. FHL2 is highly expressed in myofibroblasts in cancer tissue. TGF-ß1 upregulated FHL2 expression in CAFs and FHL2 knockdown attenuated CAF proliferation. FHL2 knockout reduced CAF induced migration of A110L and H23 human lung adenocarcinoma cell lines, and the induction of tube formation of endothelial cells. FHL2 knockout reduced CAF-induced metastasis of lung adenocarcinomas in an orthotopic model in vivo. The concentration of Osteopontin (OPN) in CM from CAF was downregulated by FHL2 knockout. siRNA silencing and antibody blocking of OPN reduced the pro-migratory effect of CM from CAF on lung cancer cells. In resected lung adenocarcinoma specimens, positive stromal FHL2 expression was significantly associated with higher microvascular density and worse prognosis. In conclusion, FHL2 expression by CAFs enhances the progression of lung adenocarcinoma by promoting angiogenesis and metastasis.

Genomic alterations associated with resistance and circulating tumor DNA dynamics for early detection of progression on CDK4/6 inhibitor in advanced breast cancer.

Combined CDK4/6 inhibitor (CDK4/6i) and endocrine therapy significantly improves outcome for patients with estrogen receptor-positive (ER+) metastatic breast cancer, but drug resistance and thus disease progression inevitably occur. Herein, we aimed to identify genomic alterations associated with combined CDK4/6i and endocrine therapy resistance, and follow the levels of specific mutations in longitudinal circulating tumor DNA (ctDNA) for early detection of progression. From a cohort of 86 patients with ER+ metastatic breast cancer we performed whole exome sequencing or targeted sequencing of paired tumor (N = 8) or blood samples (N = 5) obtained before initiation of combined CDK4/6i and endocrine therapy and at disease progression. Mutations in oncogenic genes at progression were rare, while amplifications of growth-regulating genes were more frequent. The most frequently acquired alterations observed were PIK3CA and TP53 mutations and PDK1 amplification. Longitudinal ctDNA dynamics of mutant PIK3CA or private mutations revealed increased mutation levels at progression in 8 of 10 patients (80%). Impressively, rising levels of PIK3CA-mutated ctDNA were detected 4-17 months before imaging. Our data add to the growing evidence supporting longitudinal ctDNA analysis for real-time monitoring of CDK4/6i response and early detection of progression in advanced breast cancer. Further, our analysis suggests that amplification of growth-related genes may contribute to combined CDK4/6i and endocrine therapy resistance.

Kataegis in clinical and molecular subgroups of primary breast cancer.

Kataegis is a hypermutation phenomenon characterized by localized clusters of single base pair substitution (SBS) reported in multiple cancer types. Despite a high frequency in breast cancer, large-scale analyses of kataegis patterns and associations with clinicopathological and molecular variables in established breast cancer subgroups are lacking. Therefore, WGS profiled primary breast cancers (n = 791) with associated clinical and molecular data layers, like RNA-sequencing data, were analyzed for kataegis frequency, recurrence, and associations with genomic contexts and functional elements, transcriptional patterns, driver alterations, homologous recombination deficiency (HRD), and prognosis in tumor subgroups defined by ER, PR, and HER2/ERBB2 status. Kataegis frequency was highest in the HER2-positive(p) subgroups, including both ER-negative(n)/positive(p) tumors (ERnHER2p/ERpHER2p). In TNBC, kataegis was neither associated with PAM50 nor TNBC mRNA subtypes nor with distant relapse in chemotherapy-treated patients. In ERpHER2n tumors, kataegis was associated with aggressive characteristics, including PR-negativity, molecular Luminal B subtype, higher mutational burden, higher grade, and expression of proliferation-associated genes. Recurrent kataegis loci frequently targeted regions commonly amplified in ER-positive tumors, while few recurrent loci were observed in TNBC. SBSs in kataegis loci appeared enriched in regions of open chromatin. Kataegis status was not associated with HRD in any subgroup or with distinct transcriptional patterns in unsupervised or supervised analysis. In summary, kataegis is a common hypermutation phenomenon in established breast cancer subgroups, particularly in HER2p subgroups, coinciding with an aggressive tumor phenotype in ERpHER2n disease. In TNBC, the molecular implications and associations of kataegis are less clear, including its prognostic value.

A 14-gene B-cell immune signature in early-stage triple-negative breast cancer (TNBC): a pooled analysis of seven studies.

BACKGROUND: Early-stage triple-negative breast cancer (TNBC) displays clinical and biological diversity. From a biological standpoint, immune infiltration plays a crucial role in TNBC prognosis. Currently, there is a lack of genomic tools aiding in treatment decisions for TNBC. This study aims to assess the effectiveness of a B-cell/immunoglobulin signature (IGG) alone, or in combination with tumor burden, in predicting prognosis and treatment response in patients with TNBC. METHODS: Genomic and clinical data were retrieved from 7 cohorts: SCAN-B (N = 874), BrighTNess (n = 482), CALGB-40603 (n = 389), METABRIC (n = 267), TCGA (n = 118), GSE58812 (n = 107), GSE21653 (n = 67). IGG and a risk score integrating IGG with tumor/nodal staging (IGG-Clin) were assessed for event-free survival (EFS) and overall survival (OS) in each cohort. Random effects model was used to derive pooled effect sizes. Association of IGG with pathological complete response (pCR) was assessed in CALGB-40603 and BrighTNess. Immune significance of IGG was estimated through CIBERSORTx and EcoTyper. FINDINGS: IGG was associated with improved EFS (pooled HR = 0.77, [95% CI = 0.70-0.85], I2 = 18%) and OS (pooled HR = 0.79, [0.73-0.85], I2 = 0%) across cohorts, and was predictive of pCR in CALGB-40603 (OR 1.25, [1.10-1.50]) and BrighTNess (OR 1.57 [1.25-1.98]). IGG-Clin was predictive of recurrence (pooled HR = 2.11, [1.75-2.55], I2 = 0%) and death (pooled HR = 1.99, 95% [0.84-4.73], I2 = 79%) across cohorts. IGG was associated with adaptive immune response at CIBERSORTx and EcoTyper analysis. INTERPRETATION: IGG is linked to improved prognosis and pCR in early-stage TNBC. The integration of IGG alongside tumor and nodal staging holds promise as an approach to identify patients benefitting from intensified or de-intensified treatments. FUNDING: This study received funding from: Associació Beca Marta Santamaria, European Union's Horizon 2020 research and innovation and Marie Sklodowska-Curie Actions programs, Fundación FERO, Fundación CRIS contra el cáncer, Agència de Gestó d'Ajuts Universitaris i de Recerca, Instituto de Salud Carlos III, Fundación Contigo, Asociación Cáncer de Mama Metastásico IV, Breast Cancer Research Foundation, RESCUER, Fundación científica AECC and FSEOM.

Current and emerging sequencing-based tools for precision cancer medicine.

Current precision cancer medicine is dependent on the analyses of a plethora of clinically relevant genomic aberrations. During the last decade, next-generation sequencing (NGS) has gradually replaced most other methods for precision cancer diagnostics, spanning from targeted tumor-informed assays and gene panel sequencing to global whole-genome and whole-transcriptome sequencing analyses. The shift has been impelled by a clinical need to assess an increasing number of genomic alterations with diagnostic, prognostic and predictive impact, including more complex biomarkers (e.g. microsatellite instability, MSI, and homologous recombination deficiency, HRD), driven by the parallel development of novel targeted therapies and enabled by the rapid reduction in sequencing costs. This review focuses on these sequencing-based methods, puts their emergence in a historic perspective, highlights their clinical utility in diagnostics and decision-making in pediatric and adult cancer, as well as raises challenges for their clinical implementation. Finally, the importance of applying sensitive tools for longitudinal monitoring of treatment response and detection of measurable residual disease, as well as future avenues in the rapidly evolving field of sequencing-based methods are discussed.

Disruption of the TP53 locus in osteosarcoma leads to TP53 promoter gene fusions and restoration of parts of the TP53 signalling pathway.

TP53 is the most frequently mutated gene in human cancer. This gene shows not only loss-of-function mutations but also recurrent missense mutations with gain-of-function activity. We have studied the primary bone malignancy osteosarcoma, which harbours one of the most rearranged genomes of all cancers. This is odd since it primarily affects children and adolescents who have not lived the long life thought necessary to accumulate massive numbers of mutations. In osteosarcoma, TP53 is often disrupted by structural variants. Here, we show through combined whole-genome and transcriptome analyses of 148 osteosarcomas that TP53 structural variants commonly result in loss of coding parts of the gene while simultaneously preserving and relocating the promoter region. The transferred TP53 promoter region is fused to genes previously implicated in cancer development. Paradoxically, these erroneously upregulated genes are significantly associated with the TP53 signalling pathway itself. This suggests that while the classical tumour suppressor activities of TP53 are lost, certain parts of the TP53 signalling pathway that are necessary for cancer cell survival and proliferation are retained. In line with this, our data suggest that transposition of the TP53 promoter is an early event that allows for a new normal state of genome-wide rearrangements in osteosarcoma. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Improved detection of clinically relevant fusion transcripts in cancer by machine learning classification.

BACKGROUND: Genomic rearrangements in cancer cells can create fusion genes that encode chimeric proteins or alter the expression of coding and non-coding RNAs. In some cancer types, fusions involving specific kinases are used as targets for therapy. Fusion genes can be detected by whole genome sequencing (WGS) and targeted fusion panels, but RNA sequencing (RNA-Seq) has the advantageous capability of broadly detecting expressed fusion transcripts. RESULTS: We developed a pipeline for validation of fusion transcripts identified in RNA-Seq data using matched WGS data from The Cancer Genome Atlas (TCGA) and applied it to 910 tumors from 11 different cancer types. This resulted in 4237 validated gene fusions, 3049 of them with at least one identified genomic breakpoint. Utilizing validated fusions as true positive events, we trained a machine learning classifier to predict true and false positive fusion transcripts from RNA-Seq data. The final precision and recall metrics of the classifier were 0.74 and 0.71, respectively, in an independent dataset of 249 breast tumors. Application of this classifier to all samples with RNA-Seq data from these cancer types vastly extended the number of likely true positive fusion transcripts and identified many potentially targetable kinase fusions. Further analysis of the validated gene fusions suggested that many are created by intrachromosomal amplification events with microhomology-mediated non-homologous end-joining. CONCLUSIONS: A classifier trained on validated fusion events increased the accuracy of fusion transcript identification in samples without WGS data. This allowed the analysis to be extended to all samples with RNA-Seq data, facilitating studies of tumor biology and increasing the number of detected kinase fusions. Machine learning could thus be used in identification of clinically relevant fusion events for targeted therapy. The large dataset of validated gene fusions generated here presents a useful resource for development and evaluation of fusion transcript detection algorithms.

Comparison of SP142 and 22C3 PD-L1 assays in a population-based cohort of triple-negative breast cancer patients in the context of their clinically established scoring algorithms.

BACKGROUND: Immunohistochemical (IHC) PD-L1 expression is commonly employed as predictive biomarker for checkpoint inhibitors in triple-negative breast cancer (TNBC). However, IHC evaluation methods are non-uniform and further studies are needed to optimize clinical utility. METHODS: We compared the concordance, prognostic value and gene expression between PD-L1 IHC expression by SP142 immune cell (IC) score and 22C3 combined positive score (CPS; companion IHC diagnostic assays for atezolizumab and pembrolizumab, respectively) in a population-based cohort of 232 early-stage TNBC patients. RESULTS: The expression rates of PD-L1 for SP142 IC ≥ 1%, 22C3 CPS ≥ 10, 22C3 CPS ≥ 1 and 22C3 IC ≥ 1% were 50.9%, 27.2%, 53.9% and 41.8%, respectively. The analytical concordance (kappa values) between SP142 IC+ and these three different 22C3 scorings were 73.7% (0.48, weak agreement), 81.5% (0.63) and 86.6% (0.73), respectively. The SP142 assay was better at identifying 22C3 positive tumors than the 22C3 assay was at detecting SP142 positive tumors. PD-L1 (CD274) gene expression (mRNA) showed a strong positive association with all two-categorical IHC scorings of the PD-L1 expression, irrespective of antibody and cut-off (Spearman Rho ranged from 0.59 to 0.62; all p-values < 0.001). PD-L1 IHC positivity and abundance of tumor infiltrating lymphocytes were of positive prognostic value in univariable regression analyses in patients treated with (neo)adjuvant chemotherapy, where it was strongest for 22C3 CPS ≥ 10 and distant relapse-free interval (HR = 0.18, p = 0.019). However, PD-L1 status was not independently prognostic when adjusting for abundance of tumor infiltrating lymphocytes in multivariable analyses. CONCLUSION: Our findings support that the SP142 and 22C3 IHC assays, with their respective clinically applied scoring algorithms, are not analytically equivalent where they identify partially non-overlapping subpopulations of TNBC patients and cannot be substituted with one another regarding PD-L1 detection. Trial registration The Swedish Cancerome Analysis Network - Breast (SCAN-B) study, retrospectively registered 2nd Dec 2014 at ClinicalTrials.gov; ID NCT02306096.

Perturbation and stability of PAM50 subtyping in population-based primary invasive breast cancer.

PAM50 gene expression subtypes represent a cornerstone in the molecular classification of breast cancer and are included in risk prediction models to guide therapy. We aimed to illustrate the impact of included genes and biological processes on subtyping while considering a tumor's underlying clinical subgroup defined by ER, PR, and HER2 status. To do this we used a population-representative and clinically annotated early-stage breast tumor cohort of 6233 samples profiled by RNA sequencing and applied a perturbation strategy of excluding co-expressed genes (gene sets). We demonstrate how PAM50 nearest-centroid classification depends on biological processes present across, but also within, ER/PR/HER2 subgroups and PAM50 subtypes themselves. Our analysis highlights several key aspects of PAM50 classification. Firstly, we demonstrate the tight connection between a tumor's nearest and second-nearest PAM50 centroid. Additionally, we show that the second-best subtype is associated with overall survival in ER-positive, HER2-negative, and node-negative disease. We also note that ERBB2 expression has little impact on PAM50 classification in HER2-positive disease regardless of ER status and that the Basal subtype is highly stable in contrast to the Normal subtype. Improved consciousness of the commonly used PAM50 subtyping scheme will aid in our understanding and interpretation of breast tumors that have seemingly conflicting PAM50 classification when compared to clinical biomarkers. Finally, our study adds further support in challenging the common misconception that PAM50 subtypes are distinct classes by illustrating that PAM50 subtypes in tumors represent a continuum with prognostic implications.

Diagnostic gastrointestinal markers in primary lung cancer and pulmonary metastases.

Histopathological diagnosis of pulmonary tumors is essential for treatment decisions. The distinction between primary lung adenocarcinoma and pulmonary metastasis from the gastrointestinal (GI) tract may be difficult. Therefore, we compared the diagnostic value of several immunohistochemical markers in pulmonary tumors. Tissue microarrays from 629 resected primary lung cancers and 422 resected pulmonary epithelial metastases from various sites (whereof 275 colorectal cancer) were investigated for the immunohistochemical expression of CDH17, GPA33, MUC2, MUC6, SATB2, and SMAD4, for comparison with CDX2, CK20, CK7, and TTF-1. The most sensitive markers for GI origin were GPA33 (positive in 98%, 60%, and 100% of pulmonary metastases from colorectal cancer, pancreatic cancer, and other GI adenocarcinomas, respectively), CDX2 (99/40/100%), and CDH17 (99/0/100%). In comparison, SATB2 and CK20 showed higher specificity, with expression in 5% and 10% of mucinous primary lung adenocarcinomas and both in 0% of TTF-1-negative non-mucinous primary lung adenocarcinomas (25-50% and 5-16%, respectively, for GPA33/CDX2/CDH17). MUC2 was negative in all primary lung cancers, but positive only in less than half of pulmonary metastases from mucinous adenocarcinomas from other organs. Combining six GI markers did not perfectly separate primary lung cancers from pulmonary metastases including subgroups such as mucinous adenocarcinomas or CK7-positive GI tract metastases. This comprehensive comparison suggests that CDH17, GPA33, and SATB2 may be used as equivalent alternatives to CDX2 and CK20. However, no single or combination of markers can categorically distinguish primary lung cancers from metastatic GI tract cancer.

Molecular characteristics of breast tumors in patients screened for germline predisposition from a population-based observational study.

BACKGROUND: Pathogenic germline variants (PGVs) in certain genes are linked to higher lifetime risk of developing breast cancer and can influence preventive surgery decisions and therapy choices. Public health programs offer genetic screening based on criteria designed to assess personal risk and identify individuals more likely to carry PGVs, dividing patients into screened and non-screened groups. How tumor biology and clinicopathological characteristics differ between these groups is understudied and could guide refinement of screening criteria. METHODS: Six thousand six hundred sixty breast cancer patients diagnosed in South Sweden during 2010-2018 were included with available clinicopathological and RNA sequencing data, 900 (13.5%) of which had genes screened for PGVs through routine clinical screening programs. We compared characteristics of screened patients and tumors to non-screened patients, as well as between screened patients with (n = 124) and without (n = 776) PGVs. RESULTS: Broadly, breast tumors in screened patients showed features of a more aggressive disease. However, few differences related to tumor biology or patient outcome remained significant after stratification by clinical subgroups or PAM50 subtypes. Triple-negative breast cancer (TNBC), the subgroup most enriched for PGVs, showed the most differences between screening subpopulations (e.g., higher tumor proliferation in screened cases). Significant differences in PGV prevalence were found between clinical subgroups/molecular subtypes, e.g., TNBC cases were enriched for BRCA1 PGVs. In general, clinicopathological differences between screened and non-screened patients mimicked those between patients with and without PGVs, e.g., younger age at diagnosis for positive cases. However, differences in tumor biology/microenvironment such as immune cell composition were additionally seen within PGV carriers/non-carriers in ER + /HER2 - cases, but not between screening subpopulations in this subgroup. CONCLUSIONS: Characterization of molecular tumor features in patients clinically screened and not screened for PGVs represents a relevant read-out of guideline criteria. The general lack of molecular differences between screened/non-screened patients after stratification by relevant breast cancer subsets questions the ability to improve the identification of screening candidates based on currently used patient and tumor characteristics, pointing us towards universal screening. Nevertheless, while that is not attained, molecular differences identified between PGV carriers/non-carriers suggest the possibility of further refining patient selection within certain patient subsets using RNA-seq through, e.g., gene signatures. TRIAL REGISTRATION: The Sweden Cancerome Analysis Network - Breast (SCAN-B) was prospectively registered at ClinicalTrials.gov under the identifier NCT02306096.

Molecular characteristics of lung adenocarcinoma with respect to patient age at diagnosis.

Lung cancer is primarily a disease of the elderly, with a median age at diagnosis around 70 years. In our study we sought to address the question of whether and how clinical characteristics, molecular alterations and molecular phenotypes differ between patient populations with early-stage lung adenocarcinoma (AC) with respect to age at diagnosis. Patients were stratified based on age at diagnosis into five systematic age bins (<50, 50-60, 60-70, 70-80 and ≥80 years). To assess clinicopathological variables on a population-based level, we accessed data from the national quality registry for lung cancer in Sweden. In parallel, we used compiled datasets from public cohorts to investigate focal and genome-wide DNA alterations, epigenetic alterations, immune composition and transcriptional patterns in relation to age at diagnosis. Gender, stage, WHO performance and likelihood of receiving chemotherapy as adjuvant treatment were linked to age at diagnosis. Associations between younger patient age and likelihood of harboring certain driver mutations (eg, in EGFR and ALK) were confirmed. We also found an association between age at diagnosis and certain mutational signatures. However, age did not seem to drive transcriptional, copy number, or epigenetic variation in the tumors. Based on our findings, age at diagnosis alone does not appear to provide an additional layer of biological complexity above that of proposed genetic and transcriptional phenotypes of AC.

A trans-omics assessment of gene-gene interaction in early-stage NSCLC.

Epigenome-wide gene-gene (G × G) interactions associated with non-small-cell lung cancer (NSCLC) survival may provide insights into molecular mechanisms and therapeutic targets. Hence, we proposed a three-step analytic strategy to identify significant and robust G × G interactions that are relevant to NSCLC survival. In the first step, among 49 billion pairs of DNA methylation probes, we identified 175 775 G × G interactions with PBonferroni ≤ 0.05 in the discovery phase of epigenomic analysis; among them, 15 534 were confirmed with P ≤ 0.05 in the validation phase. In the second step, we further performed a functional validation for these G × G interactions at the gene expression level by way of a two-phase (discovery and validation) transcriptomic analysis, and confirmed 25 significant G × G interactions enriched in the 6p21.33 and 6p22.1 regions. In the third step, we identified two G × G interactions using the trans-omics analysis, which had significant (P ≤ 0.05) epigenetic cis-regulation of transcription and robust G × G interactions at both the epigenetic and transcriptional levels. These interactions were cg14391855 × cg23937960 (ßinteraction  = 0.018, P = 1.87 × 10-12 ), which mapped to RELA × HLA-G (ßinteraction  = 0.218, P = 8.82 × 10-11 ) and cg08872738 × cg27077312 (ßinteraction  = -0.010, P = 1.16 × 10-11 ), which mapped to TUBA1B × TOMM40 (ßinteraction =-0.250, P = 3.83 × 10-10 ). A trans-omics mediation analysis revealed that 20.3% of epigenetic effects on NSCLC survival were significantly (P = 0.034) mediated through transcriptional expression. These statistically significant trans-omics G × G interactions can also discriminate patients with high risk of mortality. In summary, we identified two G × G interactions at both the epigenetic and transcriptional levels, and our findings may provide potential clues for precision treatment of NSCLC.

Tumor purity adjusted beta values improve biological interpretability of high-dimensional DNA methylation data.

A common issue affecting DNA methylation analysis in tumor tissue is the presence of a substantial amount of non-tumor methylation signal derived from the surrounding microenvironment. Although approaches for quantifying and correcting for the infiltration component have been proposed previously, we believe these have not fully addressed the issue in a comprehensive and universally applicable way. We present a multi-population framework for adjusting DNA methylation beta values on the Illumina 450/850K platform using generic purity estimates to account for non-tumor signal. Our approach also provides an indirect estimate of the aggregate methylation state of the surrounding normal tissue. Using whole exome sequencing derived purity estimates and Illumina 450K methylation array data generated by The Cancer Genome Atlas project (TCGA), we provide a demonstration of this framework in breast cancer illustrating the effect of beta correction on the aggregate methylation beta value distribution, clustering accuracy, and global methylation profiles.

RNA sequencing-based single sample predictors of molecular subtype and risk of recurrence for clinical assessment of early-stage breast cancer.

Multigene assays for molecular subtypes and biomarkers can aid management of early invasive breast cancer. Using RNA-sequencing we aimed to develop single-sample predictor (SSP) models for clinical markers, subtypes, and risk of recurrence (ROR). A cohort of 7743 patients was divided into training and test set. We trained SSPs for subtypes and ROR assigned by nearest-centroid (NC) methods and SSPs for biomarkers from histopathology. Classifications were compared with Prosigna in two external cohorts (ABiM, n = 100 and OSLO2-EMIT0, n = 103). Prognostic value was assessed using distant recurrence-free interval. Agreement between SSP and NC for PAM50 (five subtypes) was high (85%, Kappa = 0.78) for Subtype (four subtypes) very high (90%, Kappa = 0.84) and for ROR risk category high (84%, Kappa = 0.75, weighted Kappa = 0.90). Prognostic value was assessed as equivalent and clinically relevant. Agreement with histopathology was very high or high for receptor status, while moderate for Ki67 status and poor for Nottingham histological grade. SSP and Prosigna concordance was high for subtype (OSLO-EMIT0 83%, Kappa = 0.73 and ABiM 80%, Kappa = 0.72) and moderate and high for ROR risk category (68 and 84%, Kappa = 0.50 and 0.70, weighted Kappa = 0.70 and 0.78). Pooled concordance for emulated treatment recommendation dichotomized for chemotherapy was high (85%, Kappa = 0.66). Retrospective evaluation suggested that SSP application could change chemotherapy recommendations for up to 17% of postmenopausal ER+/HER2-/N0 patients with balanced escalation and de-escalation. Results suggest that NC and SSP models are interchangeable on a group-level and nearly so on a patient level and that SSP models can be derived to closely match clinical tests.

SRIQ clustering: A fusion of Random Forest, QT clustering, and KNN concepts.

Gene expression profiling together with unsupervised analysis methods, typically clustering methods, has been used extensively in cancer research to unravel, e.g., new molecular subtypes that hold promise of disease refinement that may ultimately benefit patients. However, many of the commonly used methods require a prespecified number of clusters to extract and frequently require some type of feature pre-selection, e.g. variance filtering. This introduces subjectivity to the process of cluster discovery and the definition of putative novel tumor subtypes. Here, we introduce SRIQ, a novel unsupervised clustering method that could circumvent some of the issues in commonly used unsupervised analysis methods. SRIQ incorporates concepts from random forest machine learning as well as quality threshold- and k-nearest neighbor clustering. It is implemented as a Java and Python pipeline including data pre-processing, differential expression analysis, and pathway analysis. Using 434 lung adenocarcinomas profiled by RNA sequencing, we demonstrate the technical reproducibility of SRIQ and benchmark its performance compared to the commonly used consensus clustering method. Based on differential gene expression analysis and auxiliary molecular data we show that SRIQ can define new tumor subsets that appear biologically relevant and consistent compared and that these new subgroups seem to refine existing transcriptional subtypes that were defined using consensus clustering. Together, this provides support that SRIQ may be a useful new tool for unsupervised analysis of gene expression data from human malignancies.

Feasibility of EBUS-TBNA for histopathological and molecular diagnostics of NSCLC-A retrospective single-center experience.

Endobronchial ultrasound-guided transbronchial needle aspiration (EBUS-TBNA) is a minimally invasive bronchoscopic procedure, well established as a diagnostic modality of first choice for diagnosis and staging of non-small cell lung cancer (NSCLC). The therapeutic decisions for advanced NSCLC require comprehensive profiling of actionable mutations, which is currently considered to be an essential part of the diagnostic process. The purpose of this study was to evaluate the utility of EBUS-TBNA cytology specimen for histological subtyping, molecular profiling of NSCLC by massive parallel sequencing (MPS), as well as for PD-L1 analysis. A retrospective review of 806 EBUS bronchoscopies was performed, resulting in a cohort of 132 consecutive patients with EBUS-TBNA specimens showing NSCLC cells in lymph nodes. Data on patient demographics, radiology features of the suspected tumor and mediastinal engagement, lymph nodes sampled, the histopathological subtype of NSCLC, and performed molecular analysis were collected. The EBUS-TBNA specimen proved sufficient for subtyping NSCLC in 83% and analysis of treatment predictive biomarkers in 77% (MPS in 53%). The adequacy of the EBUS-TBNA specimen was 69% for EGFR gene mutation analysis, 49% for analysis of ALK rearrangement, 36% for ROS1 rearrangement, and 33% for analysis of PD-L1. The findings of our study confirm that EBUS-TBNA cytology aspirate is appropriate for diagnosis and subtyping of NSCLC and largely also for treatment predictive molecular testing, although more data is needed on the utility of EBUS cytology specimen for MPS and PD-L1 analysis.

Epigenome-wide three-way interaction study identifies a complex pattern between TRIM27, KIAA0226, and smoking associated with overall survival of early-stage NSCLC.

The interaction between DNA methylation of tripartite motif containing 27 (cg05293407TRIM27 ) and smoking has previously been identified to reveal histologically heterogeneous effects of TRIM27 DNA methylation on early-stage non-small-cell lung cancer (NSCLC) survival. However, to understand the complex mechanisms underlying NSCLC progression, we searched three-way interactions. A two-phase study was adopted to identify three-way interactions in the form of pack-year of smoking (number of cigarettes smoked per day × number of years smoked) × cg05293407TRIM27 × epigenome-wide DNA methylation CpG probe. Two CpG probes were identified with FDR-q ≤ 0.05 in the discovery phase and P ≤ 0.05 in the validation phase: cg00060500KIAA0226 and cg17479956EXT2 . Compared to a prediction model with only clinical information, the model added 42 significant three-way interactions using a looser criterion (discovery: FDR-q ≤ 0.10, validation: P ≤ 0.05) had substantially improved the area under the receiver operating characteristic curve (AUC) of the prognostic prediction model for both 3-year and 5-year survival. Our research identified the complex interaction effects among multiple environment and epigenetic factors, and provided therapeutic target for NSCLC patients.