Genome Alteration Print (GAP): a tool to visualize and mine complex cancer genomic profiles obtained by SNP arrays

High level of aneuploidy and recurrent loss of chromosome 11 as relevant features of somatotroph pituitary tumors

BACKGROUND

Somatotroph neuroendocrine pituitary tumors (sPitNET) are a subtype of pituitary tumors that commonly cause acromegaly. Our study aimed to determine the spectrum of DNA copy number abnormalities (CNAs) in sPitNETs and their relevance.

METHODS

A landscape of CNAs in sPitNETs was determined using combined whole-genome approaches involving low-pass whole genome sequencing and SNP microarrays. Fluorescent in situ hybridization (FISH) was used for microscopic validation of CNAs. The tumors were also subjected to transcriptome and DNA methylation analyses with RNAseq and microarrays, respectively.

RESULTS

We observed a wide spectrum of cytogenetic changes ranging from multiple deletions, recurrent chromosome 11 loss, stable genomes, to duplication of the majority of the chromosomes. The identified CNAs were confirmed with FISH. sPitNETs with multiple duplications were characterized by intratumoral heterogeneity in chromosome number variation in individual tumor cells, as determined with FISH. These tumors were separate CNA-related sPitNET subtype in clustering analyses with CNA signature specific for whole genome doubling-related etiology. This subtype encompassed GNAS-wild type, mostly densely granulated tumors with favorable expression level of known prognosis-related genes, notably enriched with POUF1/NR5A1-double positive PitNETs. Chromosomal deletions in sPitNETs are functionally relevant. They occurred in gene-dense DNA regions and were related to genes downregulation and increased DNA methylation in the CpG island and promoter regions in the affected regions. Recurrent loss of chromosome 11 was reflected by lowered MEN1 and AIP. No such unequivocal relevance was found for chromosomal gains. Comparisons of transcriptomes of selected most cytogenetically stable sPitNETs with tumors with recurrent loss of chromosome 11 showed upregulation of processes related to gene dosage compensation mechanism in tumors with deletion. Comparison of stable tumors with those with multiple duplications showed upregulation of processes related to mitotic spindle, DNA repair, and chromatin organization. Both comparisons showed upregulation of the processes related to immune infiltration in cytogenetically stable tumors and deconvolution of DNA methylation data indicated a higher content of specified immune cells and lower tumor purity in these tumors.

CONCLUSIONS

sPitNETs fall into three relevant cytogenetic groups: highly aneuploid tumors characterized by known prognostically favorable features and low aneuploidy tumors including specific subtype with chromosome 11 loss.

Chromothripsis is rare in IDH-mutant gliomas compared to IDH-wild-type glioblastomas whereas whole-genome duplication is equally frequent in both tumor types

Background

Adult-type diffuse gliomas comprise IDH (isocitrate dehydrogenase)-mutant astrocytomas, IDH-mutant 1p/19q-codeleted oligodendrogliomas (ODG), and IDH-wild-type glioblastomas (GBM). GBM displays genome instability, which may result from 2 genetic events leading to massive chromosome alterations: Chromothripsis (CT) and whole-genome duplication (WGD). These events are scarcely described in IDH-mutant gliomas. The better prognosis of the latter may be related to their genome stability compared to GBM.

Methods

Pangenomic profiles of 297 adult diffuse gliomas were analyzed at initial diagnosis using SNP arrays, including 192 GBM and 105 IDH-mutant gliomas (61 astrocytomas and 44 ODG). Tumor ploidy was assessed with Genome Alteration Print and CT events with CTLPScanner and through manual screening. Survival data were compared using the Kaplan-Meier method.

Results

At initial diagnosis, 37 GBM (18.7%) displayed CT versus 5 IDH-mutant gliomas (4.7%; P = .0008), the latter were all high-grade (grade 3 or 4) astrocytomas. WGD was detected at initial diagnosis in 18 GBM (9.3%) and 9 IDH-mutant gliomas (5 astrocytomas and 4 oligodendrogliomas, either low- or high-grade; 8.5%). Neither CT nor WGD was associated with overall survival in GBM or in IDH-mutant gliomas.

Conclusions

CT is less frequent in IDH-mutant gliomas compared to GBM. The absence of CT in ODG and grade 2 astrocytomas might, in part, explain their genome stability and better prognosis, while CT might underlie aggressive biological behavior in some high-grade astrocytomas. WGD is a rare and early event occurring equally in IDH-mutant gliomas and GBM.

A new tumorgraft panel to accelerate precision medicine in prostate cancer

Background

Despite the significant advances in the management of advanced prostate cancer (PCa), metastatic PCa is currently considered incurable. For further investigations in precision treatment, the development of preclinical models representing the complex prostate tumor heterogeneity are mandatory. Accordingly, we aimed to establish a resource of patient-derived xenograft (PDX) models that exemplify each phase of this multistage disease for accurate and rapid evaluation of candidate therapies.

Methods

Fresh tumor samples along with normal corresponding tissues were obtained directly from patients at surgery. To ensure that the established models reproduce the main features of patient's tumor, both PDX tumors at multiple passages and patient's primary tumors, were processed for histological characteristics. STR profile analyses were also performed to confirm patient identity. Finally, the responses of the PDX models to androgen deprivation, PARP inhibitors and chemotherapy were also evaluated.

Results

In this study, we described the development and characterization of 5 new PDX models of PCa. Within this collection, hormone-naïve, androgen-sensitive and castration-resistant (CRPC) primary tumors as well as prostate carcinoma with neuroendocrine differentiation (CRPC-NE) were represented. Interestingly, the comprehensive genomic characterization of the models identified recurrent cancer driver alterations in androgen signaling, DNA repair and PI3K, among others. Results were supported by expression patterns highlighting new potential targets among gene drivers and the metabolic pathway. In addition, in vivo results showed heterogeneity of response to androgen deprivation and chemotherapy, like the responses of patients to these treatments. Importantly, the neuroendocrine model has been shown to be responsive to PARP inhibitor.

Conclusion

We have developed a biobank of 5 PDX models from hormone-naïve, androgen-sensitive to CRPC primary tumors and CRPC-NE. Increased copy-number alterations and accumulation of mutations within cancer driver genes as well as the metabolism shift are consistent with the increased resistance mechanisms to treatment. The pharmacological characterization suggested that the CRPC-NE could benefit from the PARP inhibitor treatment. Given the difficulties in developing such models, this relevant panel of PDX models of PCa will provide the scientific community with an additional resource for the further development of PDAC research.

Homologous recombination deficiency derived from whole-genome sequencing predicts platinum response in triple-negative breast cancers

The high frequency of homologous recombination deficiency (HRD) is the main rationale of testing platinum-based chemotherapy in triple-negative breast cancer (TNBC), however, the existing methods to identify HRD are controversial and there is a medical need for predictive biomarkers. We assess the in vivo response to platinum agents in 55 patient-derived xenografts (PDX) of TNBC to identify determinants of response. The HRD status, determined from whole genome sequencing, is highly predictive of platinum response. BRCA1 promoter methylation is not associated with response, in part due to residual BRCA1 gene expression and homologous recombination proficiency in different tumours showing mono-allelic methylation. Finally, in 2 cisplatin sensitive tumours we identify mutations in XRCC3 and ORC1 genes that are functionally validated in vitro. In conclusion, our results demonstrate that the genomic HRD is predictive of platinum response in a large cohort of TNBC PDX and identify alterations in XRCC3 and ORC1 genes driving cisplatin response.

Deep genomic analysis of malignant peripheral nerve sheath tumor cell lines challenges current malignant peripheral nerve sheath tumor diagnosis

Malignant peripheral nerve sheath tumors (MPNSTs) are soft-tissue sarcomas of the peripheral nervous system that develop either sporadically or in the context of neurofibromatosis type 1 (NF1). MPNST diagnosis can be challenging and treatment outcomes are poor. We present here a resource consisting of the genomic characterization of 9 widely used human MPNST cell lines for their use in translational research. NF1-related cell lines recapitulated primary MPNST copy number profiles, exhibited NF1, CDKN2A, and SUZ12/EED tumor suppressor gene (TSG) inactivation, and presented no gain-of-function mutations. In contrast, sporadic cell lines collectively displayed different TSG inactivation patterns and presented kinase-activating mutations, fusion genes, altered mutational frequencies and COSMIC signatures, and different methylome-based classifications. Cell lines re-classified as melanomas and other sarcomas exhibited a different drug-treatment response. Deep genomic analysis, methylome-based classification, and cell-identity marker expression, challenged the identity of common MPNST cell lines, opening an opportunity to revise MPNST differential diagnosis.

Detection of acquired TERT amplification in addition to predisposing p53 and Rb pathways alterations in EGFR-mutant lung adenocarcinomas transformed into small-cell lung cancers

INTRODUCTION

Among the different mechanisms of acquired resistance to EGFR tyrosine kinase inhibitors (TKIs) reported in EGFR-mutated lung adenocarcinoma (EGFR-LUAD) patients, histological transformation into small cell carcinoma (SCLC) occurs in 3-14% of resistant cases, regardless of the generation of EGFR-TKI. In recent studies, bi-allelic inactivation of TP53 and RB1 has been identified in a vast majority of transformed SCLCs. However, the molecular mechanisms driving this histologic transformation remain largely unknown, mainly due to the rarity of samples.

PATIENTS AND METHODS

Out of an initial cohort of 64 patients, tumor tissues of adequate quality and quantity for whole exome sequencing (WES) analysis were available for nine tumors for six patients: paired pre- and post-SCLC transformation samples for three Patients and post-SCLC transformation samples for three other patients.

RESULTS

Mutational analyses showed concurrent TP53 mutations and Rb pathway alterations in five of the six patients analyzed, confirming their suggested role as predisposing genetic alterations to SCLC transformation. In addition, TERT amplification was detected in four of the six patients and found to be an event acquired during SCLC transformation. Clonal history evolution analyses from the paired LUAD/SCLC samples showed different evolution patterns. In two patients, a large proportion of mutations were present in the most recent common ancestor cell of the initial LUAD and the transformed SCLC clones, whereas in the third patient, few clonal mutations were common between the LUAD and SCLC samples and the ancestor clone that lead to SCLC was present at low frequency in the initial LUAD.

CONCLUSION

Despite varied clinical presentations and clonal history evolution patterns, in addition to p53 and Rb pathways alterations, TERT amplification emerged as another common genetic mechanism of EGFR-LUAD to SCLC transformation in our cohort, and could represent a candidate therapeutic target in this subset of SCLC tumors.

Oligosarcomas, IDH-mutant are distinct and aggressive

Oligodendrogliomas are defined at the molecular level by the presence of an IDH mutation and codeletion of chromosomal arms 1p and 19q. In the past, case reports and small studies described gliomas with sarcomatous features arising from oligodendrogliomas, so called oligosarcomas. Here, we report a series of 24 IDH-mutant oligosarcomas from 23 patients forming a distinct methylation class. The tumors were recurrences from prior oligodendrogliomas or developed de novo. Precursor tumors of 12 oligosarcomas were histologically and molecularly indistinguishable from conventional oligodendrogliomas. Oligosarcoma tumor cells were embedded in a dense network of reticulin fibers, frequently showing p53 accumulation, positivity for SMA and CALD1, loss of OLIG2 and gain of H3K27 trimethylation (H3K27me3) as compared to primary lesions. In 5 oligosarcomas no 1p/19q codeletion was detectable, although it was present in the primary lesions. Copy number neutral LOH was determined as underlying mechanism. Oligosarcomas harbored an increased chromosomal copy number variation load with frequent CDKN2A/B deletions. Proteomic profiling demonstrated oligosarcomas to be highly distinct from conventional CNS WHO grade 3 oligodendrogliomas with consistent evidence for a smooth muscle differentiation. Expression of several tumor suppressors was reduced with NF1 being lost frequently. In contrast, oncogenic YAP1 was aberrantly overexpressed in oligosarcomas. Panel sequencing revealed mutations in NF1 and TP53 along with IDH1/2 and TERT promoter mutations. Survival of patients was significantly poorer for oligosarcomas as first recurrence than for grade 3 oligodendrogliomas as first recurrence. These results establish oligosarcomas as a distinct group of IDH-mutant gliomas differing from conventional oligodendrogliomas on the histologic, epigenetic, proteomic, molecular and clinical level. The diagnosis can be based on the combined presence of (a) sarcomatous histology, (b) IDH-mutation and (c) TERT promoter mutation and/or 1p/19q codeletion, or, in unresolved cases, on its characteristic DNA methylation profile.

Patient-derived xenograft models capture genomic heterogeneity in endometrial cancer

Background

Endometrial cancer (EC) is a major gynecological cancer with increasing incidence. It comprises four molecular subtypes with differing etiology, prognoses, and responses to chemotherapy. In the future, clinical trials testing new single agents or combination therapies will be targeted to the molecular subtype most likely to respond. As pre-clinical models that faithfully represent the molecular subtypes of EC are urgently needed, we sought to develop and characterize a panel of novel EC patient-derived xenograft (PDX) models.

Methods

Here, we report whole exome or whole genome sequencing of 11 PDX models and their matched primary tumor. Analysis of multiple PDX lineages and passages was performed to study tumor heterogeneity across lineages and/or passages. Based on recent reports of frequent defects in the homologous recombination (HR) pathway in EC, we assessed mutational signatures and HR deficiency scores and correlated these with in vivo responses to the PARP inhibitor (PARPi) talazoparib in six PDXs representing the copy number high/p53-mutant and mismatch-repair deficient molecular subtypes of EC.

Results

PDX models were successfully generated from grade 2/3 tumors, including three uterine carcinosarcomas. The models showed similar histomorphology to the primary tumors and represented all four molecular subtypes of EC, including five mismatch-repair deficient models. The different PDX lineages showed a wide range of inter-tumor and intra-tumor heterogeneity. However, for most PDX models, one arm recapitulated the molecular landscape of the primary tumor without major genomic drift. An in vivo response to talazoparib was detected in four copy number high models. Two models (carcinosarcomas) showed a response consistent with stable disease and two models (one copy number high serous EC and another carcinosarcoma) showed significant tumor growth inhibition, albeit one consistent with progressive disease; however, all lacked the HR deficiency genomic signature.

Conclusions

EC PDX models represent the four molecular subtypes of disease and can capture intra-tumor heterogeneity of the original primary tumor. PDXs of the copy number high molecular subtype showed sensitivity to PARPi; however, deeper and more durable responses will likely require combination of PARPi with other agents.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00990-z.

Multilayer spectrum of intratumoral heterogeneity in basal bladder cancer

Basal/squamous (Ba/Sq) subtype represents an intrinsic and robust group in the consensus molecular classification of muscle-invasive bladder cancer (MIBC), with poor outcome and controversial chemosensitivity. We aimed to investigate the spectrum of intratumor heterogeneity (ITH) in the Ba/Sq subtype. First, we validated a 29-gene NanoString CodeSet to predict the Ba/Sq subtype for FFPE samples. We identified heterogeneous Ba/Sq tumors in a series of 331 MIBC FFPE samples using dual GATA3/KRT5/6 immunohistochemistry (IHC). Heterogeneous regions with distinct immunostaining patterns were studied separately for gene expression using the 29-gene CodeSet, for mutations by targeted next-generation sequencing, and for copy number alteration (CNA) by microarray hybridization. Among 83 Ba/Sq tumors identified by GATA3/KRT5/6 dual staining, 19 tumors showed heterogeneity at the IHC level. In one third of the 19 cases, regions from the same tumor were classified in different distinct molecular subtypes. The mutational and CNA profiles confirmed the same clonal origin for IHC heterogeneous regions with possible subclonal evolution. Overall, two patterns of intratumoral heterogeneity (ITH) were observed in Ba/Sq tumors: low ITH (regions with distinct immunostaining, but common molecular subtype and shared CNA) or high ITH (regions with distinct immunostaining, molecular subtype, and CNA). These results showed multilayer heterogeneity in Ba/Sq MIBC. In view of personalized medicine, this heterogeneity adds complexity and should be taken into account for sampling procedures used for diagnosis and treatment choice. © 2021 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Integrated Genomic Analysis Identifies Driver Genes and Cisplatin-Resistant Progenitor Phenotype in Pediatric Liver Cancer

Pediatric liver cancers (PLC) comprise diverse diseases affecting infants, children, and adolescents. Despite overall good prognosis, PLCs display heterogeneous response to chemotherapy. Integrated genomic analysis of 126 pediatric liver tumors showed a continuum of driver mechanisms associated with patient age, including new targetable oncogenes. In 10% of patients with hepatoblastoma, all before three years old, we identified a mosaic premalignant clonal expansion of cells altered at the 11p15.5 locus. Analysis of spatial and longitudinal heterogeneity revealed an important plasticity between "hepatocytic," "liver progenitor," and "mesenchymal" molecular subgroups of hepatoblastoma. We showed that during chemotherapy, "liver progenitor" cells accumulated massive loads of cisplatin-induced mutations with a specific mutational signature, leading to the development of heavily mutated relapses and metastases. Drug screening in PLC cell lines identified promising targets for cisplatin-resistant progenitor cells, validated in mouse xenograft experiments. These data provide new insights into cisplatin resistance mechanisms in PLC and suggest alternative therapies. SIGNIFICANCE: PLCs are deadly when they resist chemotherapy, with limited alternative treatment options. Using a multiomics approach, we identified PLC driver genes and the cellular phenotype at the origin of cisplatin resistance. We validated new treatments targeting these molecular features in cell lines and xenografts.This article is highlighted in the In This Issue feature, p. 2355.

A high-risk retinoblastoma subtype with stemness features, dedifferentiated cone states and neuronal/ganglion cell gene expression

Retinoblastoma is the most frequent intraocular malignancy in children, originating from a maturing cone precursor in the developing retina. Little is known on the molecular basis underlying the biological and clinical behavior of this cancer. Here, using multi-omics data, we demonstrate the existence of two retinoblastoma subtypes. Subtype 1, of earlier onset, includes most of the heritable forms. It harbors few genetic alterations other than the initiating RB1 inactivation and corresponds to differentiated tumors expressing mature cone markers. By contrast, subtype 2 tumors harbor frequent recurrent genetic alterations including MYCN-amplification. They express markers of less differentiated cone together with neuronal/ganglion cell markers with marked inter- and intra-tumor heterogeneity. The cone dedifferentiation in subtype 2 is associated with stemness features including low immune and interferon response, E2F and MYC/MYCN activation and a higher propensity for metastasis. The recognition of these two subtypes, one maintaining a cone-differentiated state, and the other, more aggressive, associated with cone dedifferentiation and expression of neuronal markers, opens up important biological and clinical perspectives for retinoblastomas.

New insights from Whole Genome Sequencing: BCLAF1 deletion as a structural variant that predisposes cells towards cellular transformation

Cancer arises from a multi‑step cellular transformation process where some mutations may be inherited, while others are acquired during the process of malignant transformation. Aberrations in the BCL2 associated transcription factor 1 (BCLAF1) gene have previously been identified in patients with cancer and the aim of the present study was to identify structural variants (SVs) and the effects of BCLAF1 gene silencing on cell transformation. Whole‑genome sequencing was performed on DNA isolated from tumour biopsies with a histologically confirmed diagnosis of oesophageal squamous cell carcinoma (OSCC). Paired‑end sequencing was performed on the Illumina HiSeq2000, with 300 bp reads. Reads were aligned to the Homo sapiens reference genome (NCBI37) using ELAND and CASAVA software. SVs reported from the alignment were collated with gene loci, using the variant effect predictor of Ensembl. The affected genes were subsequently cross‑checked against the Genetic Association Database for disease and cancer associations. BCLAF1 deletion was identified as a noteworthy SV that could be associated with OSCC. Transient small interfering RNA‑mediated knockdown of BCLAF1 resulted in the altered expression of several downstream genes, including downregulation of the proapoptotic genes Caspase‑3 and BAX and the DNA damage repair genes exonuclease 1, ATR‑interacting protein and transcription regulator protein BACH1. BCLAF1 deficiency also attenuated P53 gene expression. Inhibition of BCLAF1 expression also resulted in increased colony formation. These results provide evidence that the abrogation of BCLAF1 expression results in the dysregulation of several cancer signalling pathways and abnormal cell proliferation.

Copy Number Aberration Analysis to Predict Response to Neoadjuvant Anti-HER2 Therapy: Results from the NeoALTTO Phase III Clinical Trial

PURPOSE

The heterogeneity of response to anti-HER2 agents represents a major challenge in patients with HER2-positive breast cancer. To better understand the sensitivity and resistance to trastuzumab and lapatinib, we investigated the role of copy number aberrations (CNA) in predicting pathologic complete response (pCR) and survival outcomes in the NeoALTTO trial.

EXPERIMENTAL DESIGN

The neoadjuvant phase III NeoALTTO trial enrolled 455 patients with HER2-positive early-stage breast cancer. DNA samples from 269 patients were assessed for genome-wide copy number profiling. Recurrent CNAs were found with GISTIC2.0.

RESULTS

CNA estimates were obtained for 184 patients included in NeoALTTO. Among those, matched transcriptome and whole-exome data were available for 154 and 181 patients, respectively. A significant association between gene copy number and pCR was demonstrated for ERBB2 amplification. Nevertheless, ERBB2 amplification ceased to be predictive once ERBB2 expression level was considered. GISTIC2.0 analysis revealed 159 recurrent CNA regions. Lower copy number levels of the 6q23-24 locus predicted absence of pCR in the whole cohort and in the estrogen receptor-positive subgroup. 6q23-24 deletion was significantly more frequent in TP53 wild-type (WT) compared with TP53-mutated, resulting in copy number levels significantly associated with lack of pCR only in the TP53 WT subgroup. Interestingly, a gene-ontology analysis highlighted several immune processes correlated to 6q23-24 copy number.

CONCLUSIONS

Our analysis identified ERBB2 copy number as well as 6q23-24 CNAs as predictors of response to anti-HER2-based treatment. ERBB2 expression outperformed ERBB2 amplification. The complexity of the 6q23-24 region warrants further investigation.

BRCAness, SLFN11, and RB1 loss predict response to topoisomerase I inhibitors in triple-negative breast cancers

Topoisomerase I (TOP1) inhibitors trap TOP1 cleavage complexes resulting in DNA double-strand breaks (DSBs) during replication, which are repaired by homologous recombination (HR). Triple-negative breast cancer (TNBC) could be eligible for TOP1 inhibitors given the considerable proportion of tumors with a defect in HR-mediated repair (BRCAness). The TOP1 inhibitor irinotecan was tested in 40 patient-derived xenografts (PDXs) of TNBC. BRCAness was determined with a single-nucleotide polymorphism (SNP) assay, and expression of Schlafen family member 11 (SLFN11) and retinoblastoma transcriptional corepressor 1 (RB1) was evaluated by real-time polymerase chain reaction (RT-PCR) and immunohistochemistry analyses. In addition, the combination of irinotecan and the ataxia telangiectasia and Rad3-related protein (ATR) inhibitor VE-822 was tested in SLFN11-negative PDXs, and two clinical non-camptothecin TOP1 inhibitors (LMP400 and LMP776) were tested. Thirty-eight percent of the TNBC models responded to irinotecan. BRCAness combined with high SLFN11 expression and RB1 loss identified highly sensitive tumors, consistent with the notion that deficiencies in cell cycle checkpoints and DNA repair result in high sensitivity to TOP1 inhibitors. Treatment by the ATR inhibitor VE-822 increased sensitivity to irinotecan in SLFN11-negative PDXs and abolished irinotecan-induced phosphorylation of checkpoint kinase 1 (CHK1). LMP400 (indotecan) and LMP776 (indimitecan) showed high antitumor activity in BRCA1-mutated or BRCAness-positive PDXs. Last, low SLFN11 expression was associated with poor survival in 250 patients with TNBC treated with anthracycline-based chemotherapy. In conclusion, a substantial proportion of TNBC respond to irinotecan. BRCAness, high SLFN11 expression, and RB1 loss are highly predictive of response to irinotecan and the clinical indenoisoquinoline TOP1 inhibitors.

Chromosomal translocations inactivating CDKN2A support a single path for malignant peripheral nerve sheath tumor initiation

Malignant peripheral nerve sheath tumors (MPNST) are aggressive soft tissue sarcomas with poor prognosis, developing either sporadically or in persons with neurofibromatosis type 1 (NF1). Loss of CDKN2A/B is an important early event in MPNST progression. However, many reported MPNSTs exhibit partial or no inactivation of CDKN2A/B, raising the question of whether there is more than one molecular path for MPNST initiation. We present here a comprehensive genomic analysis of MPNST cell lines and tumors to explore in depth the status of CDKN2A. After accounting for CDKN2A deletions and point mutations, we uncovered a previously unnoticed high frequency of chromosomal translocations involving CDKN2A in both MPNST cell lines and primary tumors. Most identified translocation breakpoints were validated by PCR amplification and Sanger sequencing. Many breakpoints clustered in an intronic 500 bp hotspot region adjacent to CDKN2A exon 2. We demonstrate the bi-allelic inactivation of CDKN2A in all tumors (n = 15) and cell lines (n = 8) analyzed, supporting a single molecular path for MPNST initiation in both sporadic and NF1-related MPNSTs. This general CDKN2A inactivation in MPNSTs has implications for MPNST diagnostics and treatment. Our findings might be relevant for other tumor types with high frequencies of CDKN2A inactivation.

Screening and detection of chromosomal copy number alterations in the domestic horse using SNP-array genotyping data

Chromosomal abnormalities are a common cause of infertility in horses. However, they are difficult to detect using automated methods. Here, we propose a simple methodology based on single nucleotide polymorphism (SNP)-array data that allows us to detect the main chromosomal abnormalities in horses in a single procedure. As proof of concept, we were able to detect chromosomal abnormalities in 33 out of 268 individuals, including monosomies, chimerisms, and male and female sex-reversions, by analyzing the raw signal intensity produced by an SNP array-based genotyping platform. We also demonstrated that the procedure is not affected by the SNP density of the array employed or by the inbreeding level of the individuals. Finally, the methodology proposed in this study could be performed in an open bioinformatic environment, thus permitting its integration as a flexible screening tool in diagnostic laboratories and genomic breeding programs.

Telomere length is key to hepatocellular carcinoma diversity and telomerase addiction is an actionable therapeutic target

BACKGROUND & AIMS

Telomerase activation is the earliest event in hepatocellular carcinoma (HCC) development. Thus, we aimed to elucidate the role of telomere length maintenance during liver carcinogenesis.

METHODS

Telomere length was measured in the tumor and non-tumor liver tissues of 1,502 patients (978 with HCC) and integrated with TERT alterations and expression, as well as clinical and molecular (analyzed by genome, exome, targeted and/or RNA-sequencing) features of HCC. The preclinical efficacy of anti-TERT antisense oligonucleotides (ASO) was assessed in vitro in 26 cell lines and in vivo in a xenograft mouse model.

RESULTS

Aging, liver fibrosis, male sex and excessive alcohol consumption were independent determinants of liver telomere attrition. HCC that developed in livers with long telomeres frequently had wild-type TERT with progenitor features and BAP1 mutations. In contrast, HCC that developed on livers with short telomeres were enriched in the non-proliferative HCC class and frequently had somatic TERT promoter mutations. In HCCs, telomere length is stabilized in a narrow biological range around 5.7 kb, similar to non-tumor livers, by various mechanisms that activate TERT expression. Long telomeres are characteristic of very aggressive HCCs, associated with the G3 transcriptomic subclass, TP53 alterations and poor prognosis. In HCC cell lines, TERT silencing with ASO was efficient in highly proliferative and poorly differentiated cells. Treatment for 3 to 16 weeks induced cell proliferation arrest in 12 cell lines through telomere shortening, DNA damage and activation of apoptosis. The therapeutic effect was also obtained in a xenograft mouse model.

CONCLUSIONS

Telomere maintenance in HCC carcinogenesis is diverse, and is associated with tumor progression and aggressiveness. The efficacy of anti-TERT ASO treatment in cell lines revealed the oncogenic addiction to TERT in HCC, providing a preclinical rationale for anti-TERT ASO treatment in HCC clinical trials.

LAY SUMMARY

Telomeres are repeated DNA sequences that protect chromosomes and naturally shorten in most adult cells because of the inactivation of the TERT gene, coding for the telomerase enzyme. Here we show that telomere attrition in the liver, modulated by aging, sex, fibrosis and alcohol, associates with specific clinical and molecular features of hepatocellular carcinoma, the most frequent primary liver cancer. We also show that liver cancer is dependent on TERT reactivation and telomere maintenance, which could be targeted through a novel therapeutic approach called antisense oligonucleotides.

Characterizing allele- and haplotype-specific copy numbers in single cells with CHISEL

Single-cell barcoding technologies enable genome sequencing of thousands of individual cells in parallel, but with extremely low sequencing coverage (<0.05×) per cell. While the total copy number of large multi-megabase segments can be derived from such data, important allele-specific mutations-such as copy-neutral loss of heterozygosity (LOH) in cancer-are missed. We introduce copy-number haplotype inference in single cells using evolutionary links (CHISEL), a method to infer allele- and haplotype-specific copy numbers in single cells and subpopulations of cells by aggregating sparse signal across hundreds or thousands of individual cells. We applied CHISEL to ten single-cell sequencing datasets of ~2,000 cells from two patients with breast cancer. We identified extensive allele-specific copy-number aberrations (CNAs) in these samples, including copy-neutral LOHs, whole-genome duplications (WGDs) and mirrored-subclonal CNAs. These allele-specific CNAs affect genomic regions containing well-known breast-cancer genes. We also refined the reconstruction of tumor evolution, timing allele-specific CNAs before and after WGDs, identifying low-frequency subpopulations distinguished by unique CNAs and uncovering evidence of convergent evolution.

RNA-sequencing of IDH-wild-type glioblastoma with chromothripsis identifies novel gene fusions with potential oncogenic properties

Glioblastoma (GBM) is the most frequent and most aggressive form of glioma. It is characterized by marked genomic instability, which suggests that chromothripsis (CT) might be involved in GBM initiation. Recently, CT has emerged as an alternative mechanism of cancer development, involving massive chromosome rearrangements in a one-step catastrophic event. The aim of the study was to detect CT in GBM and identify novel gene fusions in CT regions. One hundred and seventy IDH-wild-type GBM were screened for CT patterns using whole-genome single nucleotide polymorphism (SNP) arrays. RNA sequencing was performed in 52 GBM with CT features to identify gene fusions within CT regions. Forty tumors (40/52, 77%) harbored at least one gene fusion within CT regions. We identified 120 candidate gene fusions, 30 of which with potential oncogenic activities. We validated 11 gene fusions, which involved the most recurrent fusion partners (EGFR, SEPT14, VOPP1 and CPM), by RT-PCR and Sanger sequencing. The occurrence of CT points to underlying gene fusions in IDH-wild-type GBM. CT provides exciting new research avenues in this highly aggressive cancer.

Genetic alterations of SUGP1 mimic mutant-SF3B1 splice pattern in lung adenocarcinoma and other cancers

Genes involved in 3′-splice site recognition during mRNA splicing constitute an emerging class of oncogenes. SF3B1 is the most frequently mutated splicing factor in cancer, and SF3B1 mutants corrupt branchpoint recognition leading to usage of cryptic 3′-splice sites and subsequent aberrant junctions. For a comprehensive determination of alterations leading to this splicing pattern, we performed a pan-TCGA screening for SF3B1-specific aberrant acceptor usage. While the most of aberrant 3′-splice patterns were explained by SF3B1 mutations, we also detected nine SF3B1 wild-type tumors (including five lung adenocarcinomas). Genomic profile analysis of these tumors identified somatic mutations combined with loss-of-heterozygosity in the splicing factor SUGP1 in five of these cases. Modeling of SUGP1 loss and mutations in cell lines showed that both alterations induced mutant-SF3B1-like aberrant splicing. Our study provides definitive evidence that genetic alterations of SUGP1 genocopy SF3B1 mutations in lung adenocarcinoma and other cancers.

PLK1 inhibition exhibits strong anti-tumoral activity in CCND1-driven breast cancer metastases with acquired palbociclib resistance

A significant proportion of patients with oestrogen receptor (ER) positive breast cancers (BC) develop resistance to endocrine treatments (ET) and relapse with metastatic disease. Here we perform whole exome sequencing and gene expression analysis of matched primary breast tumours and bone metastasis-derived patient-derived xenografts (PDX). Transcriptomic analyses reveal enrichment of the G2/M checkpoint and up-regulation of Polo-like kinase 1 (PLK1) in PDX. PLK1 inhibition results in tumour shrinkage in highly proliferating CCND1-driven PDX, including different RB-positive PDX with acquired palbociclib resistance. Mechanistic studies in endocrine resistant cell lines, suggest an ER-independent function of PLK1 in regulating cell proliferation. Finally, in two independent clinical cohorts of ER positive BC, we find a strong association between high expression of PLK1 and a shorter metastases-free survival and poor response to anastrozole. In conclusion, our findings support clinical development of PLK1 inhibitors in patients with advanced CCND1-driven BC, including patients progressing on palbociclib treatment.

Intratumor heterogeneity of prognostic DNA-based molecular markers in adrenocortical carcinoma

BACKGROUND

The prognosis of adrenocortical carcinoma (ACC) is heterogeneous. Genomic studies have identified ACC subgroups characterized by specific molecular alterations, including features measured at DNA level (somatic mutations, chromosome alterations, DNA methylation), which are closely associated with outcome. The aim of this study was to evaluate intratumor heterogeneity of prognostic molecular markers at the DNA level.

METHODS

Two different tissue samples (primary tumor, local recurrence or metastasis) were analyzed in 26 patients who underwent surgery for primary or recurrent ACC. DNA-related biomarkers with prognostic role were investigated in frozen and paraffin-embedded samples. Somatic mutations of p53/Rb and Wnt/β-catenin pathways were assessed using next-generation sequencing (n = 26), chromosome alteration profiles were determined using SNP arrays (n = 14) and methylation profiles were determined using four-gene bisulfite pyrosequencing (n = 12).

RESULTS

Somatic mutations for ZNRF3, TP53, CTNN1B and CDKN2A were found in 7, 6, 6 and 4 patients, respectively, with intratumor heterogeneity in 8/26 patients (31%). Chromosome alteration profiles were 'Noisy' (numerous and anarchic alterations) in 8/14 and 'Chromosomal' (extended patterns of loss of heterozygosity) in 5/14 of the study samples. For these profiles, no intratumor heterogeneity was observed. Methylation profiles were hypermethylated in 5/12 and non-hypermethylated in 7/12 of the study samples. Intratumor heterogeneity of methylation profiles was observed in 2/12 patients (17%).

CONCLUSIONS

Intratumor heterogeneity impacts DNA-related molecular markers. While somatic mutation can differ, prognostic DNA methylation and chromosome alteration profile seem rather stable and might be more robust for the prognostic assessment.

Whole-genome sequencing of human malignant mesothelioma tumours and cell lines

Pleural mesothelioma is a cancer of serosal surfaces caused by environmental exposure to asbestos. Clinical outcome remains poor and while trials of new treatments are ongoing it remains an understudied cancer. Mesothelioma cell lines can readily be grown from primary tumour and from tumour cells shed into pleural effusion with the latter representing a particularly valuable source of DNA in clinical settings, procurable without the need for additional invasive procedures. However, it is not well understood how accurately patient-derived cultured tumour cells represent the molecular characteristics of their primary tumour. We used whole-genome sequencing of primary tumour and matched cultured cells to comprehensively characterize mutations and structural alterations. Most cases had complex rearranged genomes with evidence of chromoanagenesis and rearrangements reminiscent of chromoplexy. Many of the identified driver mutations were structural, indicating that mesothelioma is often caused by structural alterations and catastrophic genomic events, rather than point mutations. Because the majority of genomic changes detected in tumours were also displayed by the genomes of cultured tumour cells, we conclude that low-passage cultured tumour cells are generally suitable for molecular characterization of mesothelioma and may be particularly useful where tissue samples with high tumour cell content are not available. However, the subclonal compositions of the cell lines did not fully recapitulate the subclonal diversity of the primary tumours. Furthermore, longitudinal acquisition of major alterations in subclonal cell populations was observed after long-term passaging. These two factors define limitations of tumour-derived cell lines as genomic substrate for clinical purposes.

BAP1 mutations define a homogeneous subgroup of hepatocellular carcinoma with fibrolamellar-like features and activated PKA

BACKGROUND & AIMS

DNAJB1-PRKACA fusion is a specific driver event in fibrolamellar carcinoma (FLC), a rare subtype of hepatocellular carcinoma (HCC) that occurs in adolescents and young adults. In older patients, molecular determinants of HCC with mixed histological features of HCC and FLC (mixed-FLC/HCC) remain to be discovered.

METHODS

A series of 151 liver tumors including 126 HCC, 15 FLC, and 10 mixed-FLC/HCC were analyzed by RNAseq and whole-genome- or whole-exome sequencing. Western blots were performed to validate genomic discoveries. Results were validated using the TCGA database.

RESULTS

Most of the mixed-FLC/HCC RNAseq clustered in a robust subgroup of 17 tumors, which all had mutations or translocations inactivating BAP1, the gene encoding BRCA1-associated protein-1. Like FLC, BAP1-HCC were significantly enriched in females, patients with a lack of chronic liver disease, and fibrotic tumors compared to non-BAP1 HCC. However, patients were older and had a poorer prognosis than those with FLC. BAP1 tumors were immune hot, showed progenitor features and did not show DNAJB1-PRKACA fusion, while almost none of these tumors had mutations in CTNNB1, TP53 and TERT promoter. In contrast, 80% of the BAP1 tumors showed a chromosome gain of PRKACA at 19p13, combined with a loss of PRKAR2A (coding for the inhibitory regulatory subunit of PKA) at 3p21, leading to a high PRKACA/PRKAR2A ratio at the mRNA and protein levels.

CONCLUSION

We have characterized a subgroup of BAP1-driven HCC with fibrolamellar-like features and a dysregulation of the PKA pathway, which could be at the root of the clinical and histological similarities between BAP1 tumors and DNAJB1-PRKACA FLCs.

LAY SUMMARY

Herein, we have defined a homogeneous subgroup of hepatocellular carcinomas in which the BAP1 gene is inactivated. This leads to the development of cancers with features similar to those of fibrolamellar carcinoma. These tumors more frequently develop in females without chronic liver disease or cirrhosis. The presence of PKA activation and T cell infiltrates suggest that these tumors could be treated with PKA inhibitors or immunomodulators.

Genetic characterization of a unique neuroendocrine transdifferentiation prostate circulating tumor cell-derived eXplant model

Transformation of castration-resistant prostate cancer (CRPC) into an aggressive neuroendocrine disease (CRPC-NE) represents a major clinical challenge and experimental models are lacking. A CTC-derived eXplant (CDX) and a CDX-derived cell line are established using circulating tumor cells (CTCs) obtained by diagnostic leukapheresis from a CRPC patient resistant to enzalutamide. The CDX and the derived-cell line conserve 16% of primary tumor (PT) and 56% of CTC mutations, as well as 83% of PT copy-number aberrations including clonal TMPRSS2-ERG fusion and NKX3.1 loss. Both harbor an androgen receptor-null neuroendocrine phenotype, TP53, PTEN and RB1 loss. While PTEN and RB1 loss are acquired in CTCs, evolutionary analysis suggest that a PT subclone harboring TP53 loss is the driver of the metastatic event leading to the CDX. This CDX model provides insights on the sequential acquisition of key drivers of neuroendocrine transdifferentiation and offers a unique tool for effective drug screening in CRPC-NE management.

Clinical Interest of Combining Transcriptomic and Genomic Signatures in High-Grade Serous Ovarian Cancer

High-grade serous ovarian cancer is one of the deadliest gynecological malignancies and remains a clinical challenge. There is a critical need to effectively define patient stratification in a clinical setting. In this study, we address this question and determine the optimal number of molecular subgroups for ovarian cancer patients. By studying several independent patient cohorts, we observed that classifying high-grade serous ovarian tumors into four molecular subgroups using a transcriptomic-based approach did not reproducibly predict patient survival. In contrast, classifying these tumors into only two molecular subgroups, fibrosis and non-fibrosis, could reliably inform on patient survival. In addition, we found complementarity between transcriptomic data and the genomic signature for homologous recombination deficiency (HRD) that helped in defining prognosis of ovarian cancer patients. We also established that the transcriptomic and genomic signatures underlined independent biological processes and defined four different risk populations. Thus, combining genomic and transcriptomic information appears as the most appropriate stratification method to reliably subgroup high-grade serous ovarian cancer patients. This method can easily be transferred into the clinical setting.

CNApp, a tool for the quantification of copy number alterations and integrative analysis revealing clinical implications

Somatic copy number alterations (CNAs) are a hallmark of cancer, but their role in tumorigenesis and clinical relevance remain largely unclear. Here, we developed CNApp, a web-based tool that allows a comprehensive exploration of CNAs by using purity-corrected segmented data from multiple genomic platforms. CNApp generates genome-wide profiles, computes CNA scores for broad, focal and global CNA burdens, and uses machine learning-based predictions to classify samples. We applied CNApp to the TCGA pan-cancer dataset of 10,635 genomes showing that CNAs classify cancer types according to their tissue-of-origin, and that each cancer type shows specific ranges of broad and focal CNA scores. Moreover, CNApp reproduces recurrent CNAs in hepatocellular carcinoma and predicts colon cancer molecular subtypes and microsatellite instability based on broad CNA scores and discrete genomic imbalances. In summary, CNApp facilitates CNA-driven research by providing a unique framework to identify relevant clinical implications. CNApp is hosted at https://tools.idibaps.org/CNApp/.

Pangenomic Classification of Pituitary Neuroendocrine Tumors

Pituitary neuroendocrine tumors (PitNETs) are common, with five main histological subtypes: lactotroph, somatotroph, and thyrotroph (POU1F1/PIT1 lineage); corticotroph (TBX19/TPIT lineage); and gonadotroph (NR5A1/SF1 lineage). We report a comprehensive pangenomic classification of PitNETs. PitNETs from POU1F1/PIT1 lineage showed an epigenetic signature of diffuse DNA hypomethylation, with transposable elements expression and chromosomal instability (except for GNAS-mutated somatotrophs). In TPIT lineage, corticotrophs were divided into three classes: the USP8-mutated with overt secretion, the USP8-wild-type with increased invasiveness and increased epithelial-mesenchymal transition, and the large silent tumors with gonadotroph transdifferentiation. Unexpected expression of gonadotroph markers was also found in GNAS-wild-type somatotrophs (SF1 expression), challenging the current definition of SF1/gonadotroph lineage. This classification improves our understanding and affects the clinical stratification of patients with PitNETs.

A systematic evaluation of copy number alterations detection methods on real SNP array and deep sequencing data

BACKGROUND

The Copy Number Alterations (CNAs) are discovered to be tightly associated with cancers, so accurately detecting them is one of the most important tasks in the cancer genomics. A series of CNAs detection methods have been proposed and new ones are still being developed. Due to the complexity of CNAs in cancers, no CNAs detection method has been accepted as the gold standard caller. Several evaluation works have made attempts to reveal typical CNAs detection methods' performance. Limited by the scale of evaluation data, these different comparison works don't reach a consensus and the researchers are still confused on how to choose one proper CNAs caller for their analysis. Therefore, it needs a more comprehensive evaluation of typical CNAs detection methods' performance.

RESULTS

In this work, we use a large-scale real dataset from CAGEKID consortium to evaluate total 12 typical CNAs detection methods. These methods are most widely used in cancer researches and always used as benchmark for the newly proposed CNAs detection methods. This large-scale dataset comprises of SNP array data on 94 samples and the whole genome sequencing data on 10 samples. Evaluations are comprehensively implemented in current scenarios of CNAs detection, which include that detect CNAs on SNP array data, on sequencing data with tumor and normal matched samples and on sequencing data with single tumor sample. Three SNP based methods are firstly ranked. Subsequently, the best SNP based method's results are used as benchmark to compare six matched samples based methods and three single tumor sample based methods in terms of the preprocessing, recall rate, Jaccard index and segmentation characteristics.

CONCLUSIONS

Our survey thoroughly reveals 12 typical methods' superiority and inferiority. We explain why methods show specific characteristics from a methodological standpoint. Finally, we present the guiding principle for choosing one proper CNAs detection method under specific conditions. Some unsolved problems and expectations are also addressed for upcoming CNAs detection methods.

Integrated molecular characterization of chondrosarcoma reveals critical determinants of disease progression

Chondrosarcomas are primary cancers of cartilaginous tissue with highly contrasting prognoses. These tumors are defined by recurrent mutations in the IDH genes and other genetic alterations including inactivation of CDKN2A and COL2A1; however, these have no clinical value. Here we use multi-omics molecular profiles from a series of cartilage tumors and find an mRNA classification that identifies two subtypes of chondrosarcomas defined by a balance in tumor differentiation and cell cycle activation. The microRNA classification reveals the importance of the loss of expression of the 14q32 locus in defining the level of malignancy. Finally, DNA methylation is associated with IDH mutations. We can use the multi-omics classifications to predict outcome. We propose an mRNA-only classifier to reproduce the integrated multi-omics classification, and its application to relapsed tumor samples shows the progressive nature of the classification. Thus, it may be possible to use mRNA-based signatures to detect patients with high-risk chondrosarcomas.

LobSig is a multigene predictor of outcome in invasive lobular carcinoma

Invasive lobular carcinoma (ILC) is the most common special type of breast cancer, and is characterized by functional loss of E-cadherin, resulting in cellular adhesion defects. ILC typically present as estrogen receptor positive, grade 2 breast cancers, with a good short-term prognosis. Several large-scale molecular profiling studies have now dissected the unique genomics of ILC. We have undertaken an integrative analysis of gene expression and DNA copy number to identify novel drivers and prognostic biomarkers, using in-house (n = 25), METABRIC (n = 125) and TCGA (n = 146) samples. Using in silico integrative analyses, a 194-gene set was derived that is highly prognostic in ILC (P = 1.20 × 10-5)-we named this metagene 'LobSig'. Assessing a 10-year follow-up period, LobSig outperformed the Nottingham Prognostic Index, PAM50 risk-of-recurrence (Prosigna), OncotypeDx, and Genomic Grade Index (MapQuantDx) in a stepwise, multivariate Cox proportional hazards model, particularly in grade 2 ILC cases (χ 2, P = 9.0 × 10-6), which are difficult to prognosticate clinically. Importantly, LobSig status predicted outcome with 94.6% accuracy amongst cases classified as 'moderate-risk' according to Nottingham Prognostic Index in the METABRIC cohort. Network analysis identified few candidate pathways, though genesets related to proliferation were identified, and a LobSig-high phenotype was associated with the TCGA proliferative subtype (χ 2, P < 8.86 × 10-4). ILC with a poor outcome as predicted by LobSig were enriched with mutations in ERBB2, ERBB3, TP53, AKT1 and ROS1. LobSig has the potential to be a clinically relevant prognostic signature and warrants further development.

CEP55 is a determinant of cell fate during perturbed mitosis in breast cancer

The centrosomal protein, CEP55, is a key regulator of cytokinesis, and its overexpression is linked to genomic instability, a hallmark of cancer. However, the mechanism by which it mediates genomic instability remains elusive. Here, we showed that CEP55 overexpression/knockdown impacts survival of aneuploid cells. Loss of CEP55 sensitizes breast cancer cells to anti‐mitotic agents through premature CDK1/cyclin B activation and CDK1 caspase‐dependent mitotic cell death. Further, we showed that CEP55 is a downstream effector of the MEK1/2‐MYC axis. Blocking MEK1/2‐PLK1 signaling therefore reduced outgrowth of basal‐like syngeneic and human breast tumors in in vivo models. In conclusion, high CEP55 levels dictate cell fate during perturbed mitosis. Forced mitotic cell death by blocking MEK1/2‐PLK1 represents a potential therapeutic strategy for MYC‐CEP55‐dependent basal‐like, triple‐negative breast cancers.

Recommendations for the clinical interpretation and reporting of copy number gains using gene panel NGS analysis in routine diagnostics

Next-generation sequencing (NGS) panel analysis on DNA from formalin-fixed paraffin-embedded (FFPE) tissue is increasingly used to also identify actionable copy number gains (gene amplifications) in addition to sequence variants. While guidelines for the reporting of sequence variants are available, guidance with respect to reporting copy number gains from gene-panel NGS data is limited. Here, we report on Dutch consensus recommendations obtained in the context of the national Predictive Analysis for THerapy (PATH) project, which aims to optimize and harmonize routine diagnostics in molecular pathology. We briefly discuss two common approaches to detect gene copy number gains from NGS data, i.e., the relative coverage and B-allele frequencies. In addition, we provide recommendations for reporting gene copy gains for clinical purposes. In addition to general QC metrics associated with NGS in routine diagnostics, it is recommended to include clinically relevant quantitative parameters of copy number gains in the clinical report, such as (i) relative coverage and estimated copy numbers in neoplastic cells, (ii) statistical scores to show significance (e.g., z-scores), and (iii) the sensitivity of the assay and restrictions of NGS-based detection of copy number gains. Collectively, this information can guide clinical and analytical decisions such as the reliable detection of high-level gene amplifications and the requirement for additional in situ assays in case of borderline results or limited sensitivity.

Layer-based visualization and biomedical information exploration of multi-channel large histological data

BACKGROUND AND OBJECTIVE

Modern microscopes can acquire multi-channel large histological data from tissues of human beings or animals, which contain rich biomedical information for disease diagnosis and biological feature analysis. However, due to the large size, fuzzy tissue structure, and complicated multiple elements integrated in the image color space, it is still a challenge for current software systems to effectively calculate histological data, show the inner tissue structures and unveil hidden biomedical information. Therefore, we developed new algorithms and a software platform to address this issue.

METHODS

This paper presents a multi-channel biomedical data computing and visualization system that can efficiently process large 3D histological images acquired from high-resolution microscopes. A novelty of our system is that it can dynamically display a volume of interest and extract tissue information using a layer-based data navigation scheme. During the data exploring process, the actual resolution of the loaded data can be dynamically determined and updated, and data rendering is synchronized in four display windows at each data layer, where 2D textures are extracted from the imaging volume and mapped onto the displayed clipping planes in 3D space.

RESULTS

To test the efficiency and scalability of this system, we performed extensive evaluations using several different hardware systems and large histological color datasets acquired from a CryoViz 3D digital system. The experimental results demonstrated that our system can deliver interactive data navigation speed and display detailed imaging information in real time, which is beyond the capability of commonly available biomedical data exploration software platforms.

CONCLUSION

Taking advantage of both CPU (central processing unit) main memory and GPU (graphics processing unit) graphics memory, the presented software platform can efficiently compute, process and visualize very large biomedical data and enhance data information. The performance of this system can satisfactorily address the challenges of navigating and interrogating volumetric multi-spectral large histological image at multiple resolution levels.

Whole genome duplication is an early event leading to aneuploidy in IDH-wild type glioblastoma

Glioblastoma, the most frequent and lethal form of glioma, displays chromosome instability and recurrent somatic copy number alterations (SCNA). Chromothripsis and whole genome duplication (WGD) have been recently identified in cancer. In the present study, we analyzed SCNA and determine the ploidy pattern in 123 IDH-wild-type glioblastomas, using SNP array data. WGD and chromothripsis events were validated using, respectively, FISH and CTLPScanner. WGD was detected in 11.4% glioblastomas (14/123) and was associated with TP53 mutation (p = 0.0068). It was an early event occurring after the recurrent SCNA observed in diffuse high-grade gliomas. Glioblastomas with WGD were more aneuploid compared to glioblastomas without WGD (p < 0.0001). Chromothripsis occurred in 29.3% glioblastomas (36/123) and mostly affected chromosomes 7, 9 and 12, with amplification of oncogenes (EGFR, MDM2/CDK4), and homozygous deletion of tumor suppressor genes (CDKN2A). There was a significant association between chromothripsis and gene rearrangement at a given locus. WGD is an early genetic event significantly associated to TP53 mutation and leading to chromosome instability and aneuploidy in IDH-wild-type glioblastoma. Chromothripsis recurrently targets oncogenes and tumor suppressor genes that are key players in gliomagenesis and tumor progression. The occurrence of chromothripsis points to underlying gene rearrangements (including gene fusions), potential therapeutic targets in glioblastoma.

Breast cancer diagnosed during pregnancy is associated with enrichment of non-silent mutations, mismatch repair deficiency signature and mucin mutations

Breast cancer diagnosed during pregnancy (BCP) is a rare and highly challenging disease. To investigate the impact of pregnancy on the biology of breast cancer, we conducted a comparative analysis of a cohort of BCP patients and non-pregnant control patients by integrating gene expression, copy number alterations and whole genome sequencing data. We showed that BCP exhibit unique molecular characteristics including an enrichment of non-silent mutations, a higher frequency of mutations in mucin gene family and an enrichment of mismatch repair deficiency mutational signature. This provides important insights into the biology of BCP and suggests that these features may be implicated in promoting tumor progression during pregnancy. In addition, it provides an unprecedented resource for further understanding the biology of breast cancer in young women and how pregnancy could modulate tumor biology.

Medullary Breast Carcinoma, a Triple-Negative Breast Cancer Associated with BCLG Overexpression

Medullary breast carcinoma (MBC) is a rare subtype of triple-negative breast cancer with specific genomic features within the spectrum of basal-like carcinoma (BLC). In this study of 19 MBCs and 36 non-MBC BLCs, we refined the transcriptomic and genomic knowledge about this entity. Unsupervised and supervised analysis of transcriptomic profiles confirmed that MBC clearly differs from non-MBC BLC, with 92 genes overexpressed and 154 genes underexpressed in MBC compared with non-MBC BLC. Immunity-related pathways are the most differentially represented pathways in MBC compared with non-MBC BLC. The proapoptotic gene BCLG (official name BCL2L14) is by far the most intensely overexpressed gene in MBC. A quantitative RT-PCR validation study conducted in 526 breast tumors corresponding to all molecular subtypes documented the specificity of BCLG overexpression in MBC, which was confirmed at the protein level by immunohistochemistry. We also found that most MBCs belong to the immunomodulatory triple-negative breast cancer subtype. Using pan-genomic analysis, it was found that MBC harbors more losses of heterozygosity than non-MBC BLC. These observations corroborate the notion that MBC remains a distinct entity that could benefit from specific treatment strategies (such as deescalation or targeted therapy) adapted to this rare tumor type.

Application of the cghRA framework to the genomic characterization of Diffuse Large B-Cell Lymphoma

Motivation

Although sequencing-based technologies are becoming the new reference in genome analysis, comparative genomic hybridization arrays (aCGH) still constitute a simple and reliable approach for copy number analysis. The most powerful algorithms to analyze such data have been freely provided by the scientific community for many years, but combining them is a complex scripting task.

Results

The cghRA framework combines a user-friendly graphical interface and a powerful object-oriented command-line interface to handle a full aCGH analysis, as is illustrated in an original series of 107 Diffuse Large B-Cell Lymphomas. New algorithms for copy-number calling, polymorphism detection and minimal common region prioritization were also developed and validated. While their performances will only be demonstrated with aCGH, these algorithms could actually prove useful to any copy-number analysis, whatever the technique used.

Availability and implementation

R package and source for Linux, MS Windows and MacOS are freely available at http://bioinformatics.ovsa.fr/cghRA.

Contact

mareschal@ovsa.fr or fabrice.jardin@chb.unicancer.fr.

Supplementary information

Supplementary data are available at Bioinformatics online.

Dual origin of relapses in retinoic-acid resistant acute promyelocytic leukemia

Retinoic acid (RA) and arsenic target the t(15;17)(q24;q21) PML/RARA driver of acute promyelocytic leukemia (APL), their combination now curing over 95% patients. We report exome sequencing of 64 matched samples collected from patients at initial diagnosis, during remission, and following relapse after historical combined RA-chemotherapy treatments. A first subgroup presents a high incidence of additional oncogenic mutations disrupting key epigenetic or transcriptional regulators (primarily WT1) or activating MAPK signaling at diagnosis. Relapses retain these cooperating oncogenes and exhibit additional oncogenic alterations and/or mutations impeding therapy response (RARA, NT5C2). The second group primarily exhibits FLT3 activation at diagnosis, which is lost upon relapse together with most other passenger mutations, implying that these relapses derive from ancestral pre-leukemic PML/RARA-expressing cells that survived RA/chemotherapy. Accordingly, clonogenic activity of PML/RARA-immortalized progenitors ex vivo is only transiently affected by RA, but selectively abrogated by arsenic. Our studies stress the role of cooperating oncogenes in direct relapses and suggest that targeting pre-leukemic cells by arsenic contributes to its clinical efficacy.

Historical acute promyelocytic leukemia patients treated with retinoic acid and chemotherapy sometimes did relapse. Here the authors performed exome sequencing on 64 patient's samples from diagnosis/relapse/remission and show relapse associates either with cooperating oncogenes at diagnosis, or with unexpected persistence of ancestral pre-leukemic clones.

Distinct epigenetic landscapes underlie the pathobiology of pancreatic cancer subtypes

Recent studies have offered ample insight into genome-wide expression patterns to define pancreatic ductal adenocarcinoma (PDAC) subtypes, although there remains a lack of knowledge regarding the underlying epigenomics of PDAC. Here we perform multi-parametric integrative analyses of chromatin immunoprecipitation-sequencing (ChIP-seq) on multiple histone modifications, RNA-sequencing (RNA-seq), and DNA methylation to define epigenomic landscapes for PDAC subtypes, which can predict their relative aggressiveness and survival. Moreover, we describe the state of promoters, enhancers, super-enhancers, euchromatic, and heterochromatic regions for each subtype. Further analyses indicate that the distinct epigenomic landscapes are regulated by different membrane-to-nucleus pathways. Inactivation of a basal-specific super-enhancer associated pathway reveals the existence of plasticity between subtypes. Thus, our study provides new insight into the epigenetic landscapes associated with the heterogeneity of PDAC, thereby increasing our mechanistic understanding of this disease, as well as offering potential new markers and therapeutic targets.

Pancreatic ductal adenocarcinoma (PDAC) is a complex disease and its underlying epigenomic heterogeneity is not fully understood. Here, the authors utilize patient-derived PDAC xenografts to define the epigenomic landscape of PDAC, highlighting chromatin states linked to differing disease aggressiveness and survival.

Morphology and genomic hallmarks of breast tumours developed by ATM deleterious variant carriers

Background

The ataxia telangiectasia mutated (ATM) gene is a moderate-risk breast cancer susceptibility gene; germline loss-of-function variants are found in up to 3% of hereditary breast and ovarian cancer (HBOC) families who undergo genetic testing. So far, no clear histopathological and molecular features of breast tumours occurring in ATM deleterious variant carriers have been described, but identification of an ATM-associated tumour signature may help in patient management.

Methods

To characterise hallmarks of ATM-associated tumours, we performed systematic pathology review of tumours from 21 participants from ataxia-telangiectasia families and 18 participants from HBOC families, as well as copy number profiling on a subset of 23 tumours. Morphology of ATM-associated tumours was compared with that of 599 patients with no BRCA1 and BRCA2 mutations from a hospital-based series, as well as with data from The Cancer Genome Atlas. Absolute copy number and loss of heterozygosity (LOH) profiles were obtained from the OncoScan SNP array. In addition, we performed whole-genome sequencing on four tumours from ATM loss-of-function variant carriers with available frozen material.

Results

We found that ATM-associated tumours belong mostly to the luminal B subtype, are tetraploid and show LOH at the ATM locus at 11q22–23. Unlike tumours in which BRCA1 or BRCA2 is inactivated, tumours arising in ATM deleterious variant carriers are not associated with increased large-scale genomic instability as measured by the large-scale state transitions signature. Losses at 13q14.11-q14.3, 17p13.2-p12, 21p11.2-p11.1 and 22q11.23 were observed. Somatic alterations at these loci may therefore represent biomarkers for ATM testing and harbour driver mutations in potentially ‘druggable’ genes that would allow patients to be directed towards tailored therapeutic strategies.

Conclusions

Although ATM is involved in the DNA damage response, ATM-associated tumours are distinct from BRCA1-associated tumours in terms of morphological characteristics and genomic alterations, and they are also distinguishable from sporadic breast tumours, thus opening up the possibility to identify ATM variant carriers outside the ataxia-telangiectasia disorder and direct them towards effective cancer risk management and therapeutic strategies.

Electronic supplementary material

The online version of this article (10.1186/s13058-018-0951-9) contains supplementary material, which is available to authorized users.

Genomic hotspots but few recurrent fusion genes in breast cancer

The advent of next generation sequencing technologies has boosted the interest in exploring the role of fusion genes in the development and progression of solid tumors. In breast cancer, most of the detected gene fusions seem to be "passenger" events while the presence of recurrent and driver fusions is still under study. We performed RNA sequencing in 55 well-characterized breast cancer samples and 10 adjacent normal breast tissues, complemented by an analysis of SNP array data. We explored the presence of fusion genes and defined their association with breast cancer subtypes, clinical-pathologic characteristics and copy number aberrations. Overall, 370 fusions were detected across the majority of the samples. HER2+ samples had significantly more fusions than triple negative and luminal subtypes. The number of fusions was correlated with histological grade, Ki67 and tumor size. Clusters of fusion genes were observed across the genome and a significant correlation of fusions with copy number aberrations and more specifically amplifications was also revealed. Despite the large number of fusion events, only a few were recurrent, while recurrent individual genes forming fusions with different partners were also detected including the estrogen receptor 1 gene in the previously detected ESR1-CCDC170 fusion. Overall we detected novel gene fusion events while we confirmed previously reported fusions. Genomic hotspots of fusion genes, differences between subtypes and small number of recurrent fusions are the most relevant characteristics of these events in breast cancer. Further investigation is necessary to comprehend the biological significance of these fusions.

aCNViewer: Comprehensive genome-wide visualization of absolute copy number and copy neutral variations

Motivation

Copy number variations (CNV) include net gains or losses of part or whole chromosomal regions. They differ from copy neutral loss of heterozygosity (cn-LOH) events which do not induce any net change in the copy number and are often associated with uniparental disomy. These phenomena have long been reported to be associated with diseases and particularly in cancer. Losses/gains of genomic regions are often correlated with lower/higher gene expression. On the other hand, loss of heterozygosity (LOH) and cn-LOH are common events in cancer and may be associated with the loss of a functional tumor suppressor gene. Therefore, identifying recurrent CNV and cn-LOH events can be important as they may highlight common biological components and give insights into the development or mechanisms of a disease. However, no currently available tools allow a comprehensive whole-genome visualization of recurrent CNVs and cn-LOH in groups of samples providing absolute quantification of the aberrations leading to the loss of potentially important information.

Results

To overcome these limitations, we developed aCNViewer (Absolute CNV Viewer), a visualization tool for absolute CNVs and cn-LOH across a group of samples. aCNViewer proposes three graphical representations: dendrograms, bi-dimensional heatmaps showing chromosomal regions sharing similar abnormality patterns, and quantitative stacked histograms facilitating the identification of recurrent absolute CNVs and cn-LOH. We illustrated aCNViewer using publically available hepatocellular carcinomas (HCCs) Affymetrix SNP Array data (Fig 1A). Regions 1q and 8q present a similar percentage of total gains but significantly different copy number gain categories (p-value of 0.0103 with a Fisher exact test), validated by another cohort of HCCs (p-value of 5.6e-7) (Fig 2B).

Availability and implementation

aCNViewer is implemented in python and R and is available with a GNU GPLv3 license on GitHub https://github.com/FJD-CEPH/aCNViewer and Docker https://hub.docker.com/r/fjdceph/acnviewer/.

Contact

aCNViewer@cephb.fr

Copy number profiles of paired primary and metastatic colorectal cancers

Liver metastasis is the major cause of death following a diagnosis of colorectal cancer (CRC). In this study, we compared the copy number profiles of paired primary and liver metastatic CRC to better understand how the genomic structure of primary CRC differs from the metastasis. Paired primary and metastatic tumors from 16 patients and their adjacent normal tissue samples were analyzed using single nucleotide polymorphism arrays. Genome-wide chromosomal copy number alterations were assessed, with particular attention to 188 genes known to be somatically altered in CRC and 24 genes that are clinically actionable in CRC. These data were analyzed with respect to the timing of primary and metastatic tissue resection and with exposure to chemotherapy. The genomic differences between the tumor and paired metastases revealed an average copy number discordance of 22.0%. The pairs of tumor samples collected prior to treatment revealed significantly higher copy number differences compared to post-therapy liver metastases (P = 0.014). Loss of heterozygosity acquired in liver metastases was significantly higher in previously treated liver metastasis samples compared to treatment naive liver metastasis samples (P = 0.003). Amplification of the clinically actionable genes ERBB2, FGFR1, PIK3CA or CDK8 was observed in the metastatic tissue of 4 patients but not in the paired primary CRC. These examples highlight the intra-patient genomic discrepancies that can occur between metastases and the primary tumors from which they arose. We propose that precision medicine strategies may therefore identify different actionable targets in metastatic tissue, compared to primary tumors, due to substantial genomic differences.

Pan-cancer analysis of homozygous deletions in primary tumours uncovers rare tumour suppressors

Homozygous deletions are rare in cancers and often target tumour suppressor genes. Here, we build a compendium of 2218 primary tumours across 12 human cancer types and systematically screen for homozygous deletions, aiming to identify rare tumour suppressors. Our analysis defines 96 genomic regions recurrently targeted by homozygous deletions. These recurrent homozygous deletions occur either over tumour suppressors or over fragile sites, regions of increased genomic instability. We construct a statistical model that separates fragile sites from regions showing signatures of positive selection for homozygous deletions and identify candidate tumour suppressors within those regions. We find 16 established tumour suppressors and propose 27 candidate tumour suppressors. Several of these genes (including MGMT, RAD17, and USP44) show prior evidence of a tumour suppressive function. Other candidate tumour suppressors, such as MAFTRR, KIAA1551, and IGF2BP2, are novel. Our study demonstrates how rare tumour suppressors can be identified through copy number meta-analysis.

Homozygous deletions are rare in cancers and often target tumour suppressor genes. Here, the authors conduct pan-cancer analyses and apply statistical modelling to identify 27 candidate tumour suppressors, including MAFTRR, KIAA1551, and IGF2BP2.

Gene expression profiling of patient-derived pancreatic cancer xenografts predicts sensitivity to the BET bromodomain inhibitor JQ1: implications for individualized medicine efforts

c-MYC controls more than 15% of genes responsible for proliferation, differentiation, and cellular metabolism in pancreatic as well as other cancers making this transcription factor a prime target for treating patients. The transcriptome of 55 patient-derived xenografts show that 30% of them share an exacerbated expression profile of MYC transcriptional targets (MYC-high). This cohort is characterized by a high level of Ki67 staining, a lower differentiation state, and a shorter survival time compared to the MYC-low subgroup. To define classifier expression signature, we selected a group of 10 MYC target transcripts which expression is increased in the MYC-high group and six transcripts increased in the MYC-low group. We validated the ability of these markers panel to identify MYC-high patient-derived xenografts from both: discovery and validation cohorts as well as primary cell cultures from the same patients. We then showed that cells from MYC-high patients are more sensitive to JQ1 treatment compared to MYC-low cells, in monolayer, 3D cultured spheroids and in vivo xenografted tumors, due to cell cycle arrest followed by apoptosis. Therefore, these results provide new markers and potentially novel therapeutic modalities for distinct subgroups of pancreatic tumors and may find application to the future management of these patients within the setting of individualized medicine clinics.

Combined tumor genomic profiling and exome sequencing in a breast cancer family implicates ATM in tumorigenesis: A proof of principle study

Familial breast cancers (BCs) account for 10%-20% of all diagnosed BCs, yet only 20% of such tumors arise in the context of a germline mutation in known tumor suppressor genes such as BRCA1 or BRCA2. The vast genetic heterogeneity which characterizes non BRCA1 and non BRCA2 (or BRCAx) families makes grouped studies impossible to perform. Next generation sequencing techniques, however, allow individual families to be studied to identify rare and or private mutations but the high number of genetic variants identified need to be sorted using pathogenicity or recurrence criteria. An additional sorting criterion may be represented by the identification of candidate regions defined by tumor genomic rearrangements. Indeed, comparative genomic hybridization (CGH) using single nucleotide polymorphism (SNP) arrays allows the detection of conserved ancestral haplotypes within recurrent regions of loss of heterozygosity, common to several familial tumors, which can highlight genomic loci harboring a germline mutation in cancer predisposition genes. The combination of both exome sequencing and SNP array-CGH for a series of familial BC revealed a germline ATM mutation associated with a loss of the wild-type allele in two BC from a BRCAx family. The analysis of additional breast tumors from ten BC families in which a germline ATM mutation had been identified revealed a high frequency of wild-type allele loss. This result argues strongly in favor of the involvement of ATM in these tumors as a tumor suppressor gene and confirms that germline ATM mutations are involved in a subset of familial BC.

Clonal composition of human ovarian cancer based on copy number analysis reveals a reciprocal relation with oncogenic mutation status

Intratumoral heterogeneity of cancer cells remains largely unexplored. Here we investigated the composition of ovarian cancer and its biological relevance. A whole-genome single nucleotide polymorphism array was applied to detect the clonal composition of 24 formalin-fixed, paraffin-embedded samples of human ovarian cancer. Genome-wide segmentation data consisting of the log2 ratio (log2R) and B allele frequency (BAF) were used to calculate an estimate of the clonal composition number (CC number) for each tumor. Somatic mutation profiles of cancer-related genes were also determined for the same 24 samples by next-generation sequencing. The CC number was estimated successfully for 23 of the 24 cancer samples. The mean ± SD value for the CC number was 1.7 ± 1.1 (range of 0-4). A somatic mutation in at least one gene was identified in 22 of the 24 ovarian cancer samples, with the mutations including those in the oncogenes KRAS (29.2%), PIK3CA (12.5%), BRAF (8.3%), FGFR2 (4.2%), and JAK2 (4.2%) as well as those in the tumor suppressor genes TP53 (54.2%), FBXW7 (8.3%), PTEN (4.2%), and RB1 (4.2%). Tumors with one or more oncogenic mutations had a significantly lower CC number than did those without such a mutation (1.0 ± 0.8 versus 2.3 ± 0.9, P = 0.0027), suggesting that cancers with driver oncogene mutations are less heterogeneous than those with other mutations. Our results thus reveal a reciprocal relation between oncogenic mutation status and clonal composition in ovarian cancer using the established method for the estimation of the CC number.