Analysis Of Copy Number Alterations Research Articles

Abstract P4-01-10: Targeted exome sequencing of circulating tumor cells in breast cancer

Abstract Background: Next-generation sequencing (NGS) has provided the means to comprehensively analyze somatic mutations in solid tumors. The major challenge is the accurate detection of somatic variants in low-quality DNA from formalin-fixed paraffin-embedded (FFPE) tissues and circulating tumor cells (CTCs) as white blood cells (WBCs) or other background cell populations confound variant identification. The aim of our study was to evaluate the feasibility of combining whole genome amplification (WGA) with hybridization capture-based, targeted exome sequencing in identifying clinically relevant variants in low-quality and limited input metastatic breast cancer (MBC) samples. Methods: For clinical samples, we collected 7.5 mL of blood in Streck tubes and FFPE biopsies from 5 newly diagnosed MBC patients. The CTCs were isolated with the ANGLE Parsortix using a 10um cassette within 24h of blood draw. WBCs from the same patients were profiled to eliminate germline contaminants. We used spiked cell lines (CTC mimics) as a positive control. MDA-MB-231 cells (n=50) were spiked into 7.5ml healthy donor whole-blood. We used peripheral blood not processed on the Parsortix as a negative control (no spiked cells). We performed exome sequencing on a total of 20 samples (experimental samples: spike-in MDA-MB-231 samples (S1 and S2, n=2), positive control bulk MDA-MB-231 cells (P1 and P2, n=2) and negative control WBCs (WL1 and WL2, n=2); clinical samples: CTCs (n=5), MBC (n=4) and WBCs (n=5)). We used the Repli-G WGA method and the SeqCap EZ MedExome Probes that targets the human exome for ~4600 medically relevant genes, followed by sequencing using the Illumina HiSeq2500. Somatic mutations were called with GATK MuTect2 in clinical samples and Samtools using Partek Flow for experimental samples. The variant annotation was performed using SnpEff. Results: The sequencing quality of all 6 spike-in experimental samples showed that >75% of the reads were overlapping in genomic features. The total number of variants identified in P1 and P2 were 69,954 and 69,848, respectively. The number of variants identified in S1 and S2 was 51,359 and 50,798 respectively. We found a substantial overlap of 67% (47,201/ 69,954) common variants in P1 and S1. Likewise, an overlap of 64% common variants (45,392/69,848) was detected in P2 and S2. We analyzed the pattern of somatic variants in both coding and noncoding regions, resulting in the identification of 8786 variants in MBC (n=4) and 638 variants in CTCs (n=5). With further filtering steps, we detected 1959 coding variants in MBC and 118 variants in CTCs. SnpEff predicted that 74% of these variants in CTCs (88/118) and 14% in MBC (274/1959) had a high impact on protein function. A total of 25 common genes were also identified in CTCs and MBC in which variants were found in protein-coding regions. We matched a list of genes having coding variants in CTCs and MBC with the actionable gene panel of a commercial genomic testing platform (FoundationOne) and a precision oncology database (OncoKB) cancer gene list. The CTCs (3/5 samples) harbored alterations in a total of 11 FoundationOne genes and 9 OncoKB genes. We identified 54 actionable genes that showed overlap with FoundationOne and 104 genes with OncoKB in all 4 MBC samples. The copy number alterations (CNAs) analysis showed a larger number of amplifications in CTCs and deletions in MBC. Conclusion: We established the feasibility of exome sequencing of medically relevant genes in detecting known and unknown variants in a breast cancer cell line spike in experiment. This workflow has further been validated on CTCs and MBC samples by sequencing both low-input material and low-quality FFPE samples. This approach has a high potential for the identification of the complete repertoire of mutations along with CNAs and for monitoring the presence of actionable alterations throughout the course of treatment. Citation Format: Pushpinder K. Bains, Daniel Campo, Tania B. Porras, Janice Lu, Julie E. Lang. Targeted exome sequencing of circulating tumor cells in breast cancer [abstract]. In: Proceedings of the 2019 San Antonio Breast Cancer Symposium; 2019 Dec 10-14; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2020;80(4 Suppl):Abstract nr P4-01-10.

Comprehensive analysis of somatic copy number alterations in clear cell renal cell carcinoma.

Somatic copy number alterations (SCNAs) are important biological characteristics that can identify genome‐wide alterations in renal cell carcinoma (RCC). Recent studies have shown that SCNAs have potential value for determining the prognosis of RCC. We examined SCNAs using the Affymetrix platform to analyze samples from 59 patients with clear cell RCCs (ccRCCs) including first cohort (30 cases) and second cohort (validation cohort, 29 cases). We stratified SCNAs in the ccRCCs using a hierarchical cluster analysis based on SCNA types, including gain, loss of heterozygosity (LOH), copy neutral LOH, mosaic, and mixed types. In this way, the examined two cohorts were categorized into two subgroups (1 and 2). Although the frequency of mixed type was higher in subgroup 1 than in subgroup 2 in the two cohorts, the association did not reach statistical significance. There was a significant difference in the frequency of metachronous metastasis between subgroups 1 and 2 (subgroup 2 > 1). In addition, subgroup 2 was retained in multivariate analysis of both cohorts. We examined whether there were specific alleles differing between subgroups 1 and 2 in both cohorts. We found that there was indeed a statistically significant difference in the 3p mixed types. Among the 3p mixed type, we found that 3p24.3 mixed type was inversely correlated with the presence of metachronous metastasis in ccRCC. The association was also retained in multivariate analysis in second cohort. We suggest that the 3p24.3 mixed type may be a novel marker to predict a favorable prognosis in ccRCC.

Low-Pass Whole-Genome Sequencing as a Method of Determining Copy Number Variations in Uveal Melanoma Tissue Samples

Analysis of specific somatic copy number alterations (SCNAs) using multiplex ligation-dependent probe amplification (MLPA) is used routinely as a prognostic test for uveal melanoma (UM). This technique requires relatively large amounts of input DNA, unattainable from many small fine-needle aspirate biopsy specimens. Herein, we compared the use of MLPA with whole-genome amplification (WGA) combined with low-pass whole-genome sequencing (LP-WGS) for detection of SCNA profiles in UM biopsy specimens. DNA was extracted from 21 formalin-fixed, paraffin-embedded UM samples and SCNAs were assessed using MLPA and WGA followed by LP-WGS. Cohen's κ was used to assess the concordance of copy number calls of each individual chromosome arm for each patient. MLPA and WGA/LP-WGS detection of SCNAs in chromosomes 1p, 3, 6, and 8 were compared and found to be highly concordant with a Cohen's κ of 0.856 (bias-corrected and accelerated 95% CI, 0.770-0.934). Only 13 of 147 (8.8%) chromosomal arms investigated resulted in discordant calls, predominantly SCNAs detected by WGA/LP-WGS but not MLPA. These results indicate that LP-WGS might be a suitable alternative or adjunct to MLPA for the detection of SCNAs associated with prognosis of UM, for cases with limiting tissue or DNA yields.

Genome-wide somatic copy number alteration analysis and database construction for cervical cancer.

Cervical cancer is a common gynecological malignancy with high incidence and mortality. Somatic copy number alterations (CNAs) play an important role in identifying tumor suppressor genes and oncogenes and are a useful diagnostic indicator for many cancer types. However, the genomic landscape of CNAs in cervical cancer has not yet been comprehensively characterized. In the present study, we collected 974 cervical cancer samples from different data sources. All samples were analyzed by genomic arrays to obtain high-resolution CNAs. Focal genomic regions with CNA events and potential cancer driver genes were identified by GISTIC2.0. Meanwhile, we constructed a comprehensive cervical cancer database by PHP and self-written Perl and R scripts. In total, 54 recurrent regions of amplification and deletion were detected. Frequently altered tumor suppressor genes were found in these regions, including PIK3CA, ERBB2, EP300 and FBXW7. CNA hotspots and related enriched functional categories were also identified. The incidence of chromothripsis in cervical cancer was estimated to be 6.06%, and the chromosome pulverization hotspot regions were detected. Based on the curated data, we developed CNAdbCC (http://cailab.labshare.cn/CNAdbCC/), a comprehensive database for copy number alterations in cervical cancer. We provide a user-friendly Web interface for data mining and visualization. It is the most comprehensive public database devoted exclusively to genomic alterations in cervical cancer. These results extend our molecular understanding of cervical cancer. The database will enable researchers to explore specific CNA patterns in this lethal cancer and facilitate the discovery of therapeutic candidates.

Correlation of HER2, MDM2, c-MYC, c-MET, and TP53 Copy Number Alterations in Circulating Tumor Cells with Tissue in Gastric Cancer Patients: A Pilot Study

Background:The analysis of the gene copy number alterations in tumor samples are increasingly used for diagnostic and prognostic purposes in patients with GC. However, these procedures are not always applicable due to their invasive nature. In this study, we have analyzed the copy number alterations of five genes (HER2, MDM2, c-MYC, c-MET, and TP53) with a fixed relevance for GC in the CTC of GC patients, and, accordingly, as a potential approach, evaluated their usage to complete primary tumor biopsy. Methods:We analyzed the status of the copy number alterations of the selected genes in CTC and matched biopsy tissues from 37 GC patients using FISH. Results: HER2 amplification was observed in 2 (5.41%) samples. HER2 gene status in CTC showed a strong agreement with its status in 36 out of 37 patients’ matched tissue samples (correlation: 97.29%; Kappa: 0.65; p < 0.001). MDM2 amplification was found only in 1 (2.70%) sample; however, the amplification of this gene was not detectable in the CTC isolated from this patient. c-MYC amplification was observed in 3 (8.11%) samples, and the status of its amplification in the CTC indicated a complete agreement with its status in the matched tissue samples (correlation: 100%; Kappa: 1.0). Conclusion:Our work suggests that the amplification of HER2 and c-MYC is in concordance with the CTC and achieved biopsies, and, consequently, CTC may act as a non-invasive alternative for recording the amplification of these genes among GC patients.

Comprehensive Bioinformatics Analysis of Key Methyltransferases and Demethylases for Histone Lysines in Hepatocellular Carcinoma.

Background & Aims:Methylation of lysines on histones, controlled by various methyltransferases and demethylases, is an important component of epigenetic modifications, and abnormal regulation of such enzymes serves as common events in hepatocellular carcinoma. We determined to identify important methyltransferases and demethylases that might regulate the development of hepatocellular carcinoma by bioinformatics.Methods:The Oncomine and UALCAN databases were used to retrieve mRNA expression levels of histone lysine methyltransferases and demethylases in hepatocellular carcinoma. Data analyses of genetic alterations, mainly mutations and copy number alterations, were performed on the cBioportal platform. Protein-protein interactions were established in the STRING database.Results:mRNA expression of 8 genes correlated with clinical staging and grading, whereas 4 genes indicated a role in the prognosis, all co-expressed with SEDB1 and WHSC1. Genetically, 12 genes showing an alteration rate higher than 5% were identified, and only 3 were indicative of prognosis. Copy number gains in ASH1L, SETDB1, and KDM5B might partially contribute to the upregulation of their mRNA expression. The close relationship of mutations in MLL2/MLL3 with driver gene mutations in hepatocellular carcinoma provided a rationale for further investigation.Conclusions:We identified 11 methyltransferases and demethylases for major histone lysines that might be promising research targets in the pathogenesis, development, and prediction of prognosis in hepatocellular carcinoma using bioinformatics.

TERT promoter mutation associated with multifocal phenotype and poor prognosis in patients with IDH wild-type glioblastoma.

BackgroundAlthough mutations in the promoter region of the telomerase reverse transcriptase (TERTp) gene are the most common alterations in glioblastoma (GBM), their clinical significance remains unclear. Therefore, we investigated the impact of TERTp status on patient outcome and clinicopathological features in patients with GBM over a long period of follow-up.MethodsWe retrospectively analyzed 153 cases of GBM. Six patients with isocitrate dehydrogenase 1 (IDH1) or H3F3A gene mutations were excluded from this study. Among the 147 cases of IDH wild-type GBM, 92 (62.6%) had the TERTp mutation. Clinical, immunohistochemical, and genetic factors (BRAF, TP53 gene mutation, CD133, ATRX expression, O6-methylguanine-DNA methyltransferase [MGMT] promoter methylation) and copy number alterations (CNAs) were investigated.ResultsGBM patients with the TERTp mutation were older at first diagnosis versus those with TERTp wild type (66.0 vs. 60.0 years, respectively, P = .034), and had shorter progression-free survival (7 vs. 10 months, respectively, P = .015) and overall survival (16 vs. 24 months, respectively, P = .017). Notably, magnetic resonance imaging performed showed that TERTp-mutant GBM was strongly associated with multifocal/distant lesions (P = .004). According to the CNA analysis, TERTp mutations were positively correlated with EGFR amp/gain, CDKN2A deletion, and PTEN deletion; however, these mutations were negatively correlated with PDGFR amp/gain, CDK4 gain, and TP53 deletion.Conclusions TERTp mutations were strongly correlated with multifocal/distant lesions and poor prognosis in patients with IDH wild-type GBM. Less aggressive GBM with TERTp wild type may be a distinct clinical and molecular subtype of IDH wild-type GBM.

A unique mutator phenotype reveals complementary oncogenic lesions leading to acute leukemia.

Mice homozygous for a hypomorphic allele of DNA replication factor minichromosome maintenance protein 2 (designated Mcm2cre/cre) develop precursor T cell lymphoblastic leukemia/lymphoma (pre-T LBL) with 4-32 small interstitial deletions per tumor. Mice that express a NUP98-HOXD13 (NHD13) transgene develop multiple types of leukemia, including myeloid and T and B lymphocyte. All Mcm2cre/cre NHD13+ mice develop pre-T LBL, and 26% develop an unrelated, concurrent B cell precursor acute lymphoblastic leukemia (BCP-ALL). Copy number alteration (CNA) analysis demonstrated that pre-T LBLs were characterized by homozygous deletions of Pten and Tcf3 and partial deletions of Notch1 leading to Notch1 activation. In contrast, BCP-ALLs were characterized by recurrent deletions involving Pax5 and Ptpn1 and copy number gain of Abl1 and Nup214 resulting in a Nup214-Abl1 fusion. We present a model in which Mcm2 deficiency leads to replicative stress, DNA double strand breaks (DSBs), and resultant CNAs due to errors in DNA DSB repair. CNAs that involve critical oncogenic pathways are then selected in vivo as malignant lymphoblasts because of a fitness advantage. Some CNAs, such as those involving Abl1 and Notch1, represent attractive targets for therapy.

Molecular Characterization of Locally Relapsed Head and Neck Cancer after Concomitant Chemoradiotherapy.

To investigate the pathobiological origin of local relapse after chemoradiotherapy, we studied genetic relationships of primary tumors (PT) and local relapses (LR) of patients treated with chemoradiotherapy. First, low-coverage whole genome sequencing was performed on DNA from 44 biopsies of resected head and neck squamous cell carcinoma (HNSCC) specimens (median 3 biopsies/tumor) to assess suitability of copy number alterations (CNAs) as biomarker for genetic relationships. CNAs were compared within and between tumors and an algorithm was developed to assess genetic relationships with consideration of intratumor heterogeneity. Next, this CNA-based algorithm was combined with target enrichment sequencing of genes frequently mutated in HNSCC to assess the genetic relationships of paired tumors and LRs of patients treated with chemoradiotherapy. Genetic relationship analysis using CNAs could accurately (96%) predict tumor biopsy pairs as patient-matched or independent. However, subsequent CNA analysis of PTs and LRs after chemoradiotherapy suggested genetic relationships in only 20% of cases, and absence in 80%. Target enrichment sequencing for mutations confirmed absence of any genetic relationship in half of the paired PTs and LRs. There are minor variations in CNA profiles within different areas of HNSCC tumors and many between independent tumors, suggesting that CNA profiles could be exploited as a marker of genetic relationship. Using CNA profiling and mutational analysis of cancer driver genes, relapses after chemoradiotherapy appear to be partially genetically related to the corresponding PTs, but seem often genetically unrelated. This remarkable observation warrants further studies and will impact therapeutic innovations and prognostic modeling when using index tumor characteristics.

Euroclonality-NGS DNA Capture Panel for Integrated Analysis of IG/TR Rearrangements, Translocations, Copy Number and Sequence Variation in Lymphoproliferative Disorders

Introduction: Current diagnostic standards for lymphoproliferative disorders include detection of clonal immunoglobulin (IG) and/or T cell receptor (TR) rearrangements, translocations, copy number alterations (CNA) and somatic mutations. These analyses frequently require a series of separate tests such as clonality PCR, fluorescence in situ hybridisation and/or immunohistochemistry, MLPA or SNParrays and sequencing. The EuroClonality-NGS DNA capture (EuroClonality-NDC) panel, developed by the EuroClonality-NGS Working Group, was designed to characterise all these alterations by capturing variable, diversity and joining IG and TR genes along with additional clinically relevant genes for CNA and mutation analysis. Methods: Well characterised B and T cell lines (n=14) representing a diverse repertoire of IG/TR rearrangements were used as a proficiency assessment to ensure 7 testing EuroClonality centres achieved optimal sequencing performance using the EuroClonality-NDC optimised and standardised protocol. A set of 56 IG/TR rearrangements across the 14 cell lines were compiled based on detection by Sanger, amplicon-NGS and capture-NGS sequencing technologies. For clinical validation of the NGS panel, clinical samples representing both B and T cell malignancies (n=280), with ≥ 5% tumour infiltration were collected from 10 European laboratories, with 88 (31%) being formalin fixed paraffin-embedded samples. Samples were distributed to the 7 centres for library preparation, hybridisation with the EuroClonality-NDC panel and sequencing on a NextSeq 500, using the EuroClonality-NDC standard protocol. Sequencing data were analysed using a customised version of ARResT/Interrogate, with independent review of the results by 2 centres. All cases exhibiting discordance between the benchmark and capture NGS results were submitted to an internal review committee comprising members of all participating centres. Results: All 7 testing centres detected all 56 rearrangements of the proficiency assessment and continued through to the validation phase. A total of 10/280 (3.5%) samples were removed from the validation analysis due to NGS failures (n=1), tumour infiltration &amp;lt; 5% (n=7), and sample misidentification (n=2). The EuroClonality-NDC panel detected B cell clonality (i.e. detection of at least one clonal rearrangement at IGH, IGK or IGL loci) in 189/197 (96%) B cell malignancies. Seven of the 8 discordant cases were post-germinal centre malignancies exhibiting Ig somatic hypermutation. The EuroClonality-NDC panel detected T cell clonality (i.e. detection of at least one clonal rearrangement at TRA, TRB, TRD or TRG loci) in 70/73 (96%) T cell malignancies. In all 3 discordant cases analysis of benchmark PCR data was not able to detect clonality at any TR loci. Next, we examined whether the EuroClonality-NDC panel could detect clonality at each of the individual loci, resulting in sensitivity values of 95% or higher for all IG/TR loci, with the exception of those where limited benchmark data were available, i.e. IGL (n=3) and TRA (n=7). The specificity of the panel was assessed on benign reactive lesions (n=21) that did not contain clonal IG/TR rearrangements based on BIOMED-2/EuroClonality PCR results; no clonality was observed by EuroClonality-NDC in any of the 21 cases. Limit of detection (LOD) assessment to detect IG/TR rearrangements was performed using cell line blends with each of the 7 centres receiving blended cell lines diluted to 10%, 5.0%, 2.5% and 1.25%. Across all 7 centres the overall detection rate was 100%, 94.1%, 76.5% and 32.4% respectively, giving an overall LOD of 5%. Sufficient data were available in 239 samples for the analysis of translocations. The correct translocation was detected in 137 out of 145 cases, resulting in a sensitivity of 95%. Table 1 shows how translocations identified by the EuroClonality-NDC protocol were restricted to disease subtypes known to harbour those types of translocations. Analysis of CNA and somatic mutations in all samples is underway and will be presented at the meeting. Conclusions: The EuroClonality-NDC panel, with an optimised laboratory protocol and bioinformatics pipeline, detects IG and TR rearrangements and translocations with high sensitivity and specificity with a LOD ≤ 5% and provides a single end-to-end workflow for the simultaneous detection of IG/TR rearrangements, translocations, CNA and sequence variants. Table. Disclosures Stamatopoulos: Janssen: Honoraria, Research Funding; Abbvie: Honoraria, Research Funding. Klapper:Roche, Takeda, Amgen, Regeneron: Honoraria, Research Funding. Ferrero:Gilead: Speakers Bureau; Janssen: Consultancy, Membership on an entity's Board of Directors or advisory committees, Speakers Bureau; EUSA Pharma: Membership on an entity's Board of Directors or advisory committees; Servier: Speakers Bureau. van den Brand:Gilead: Speakers Bureau. Groenen:Gilead: Speakers Bureau. Brüggemann:Incyte: Membership on an entity's Board of Directors or advisory committees; Amgen: Membership on an entity's Board of Directors or advisory committees; Roche: Consultancy. Langerak:Gilead: Research Funding, Speakers Bureau; F. Hoffmann-La Roche Ltd: Research Funding; Genentech, Inc.: Research Funding; Janssen: Speakers Bureau. Gonzalez:Roche: Honoraria, Research Funding; AstraZeneca: Consultancy, Honoraria, Research Funding, Speakers Bureau.

The Genomic and Transcriptomic Landscape of Double-Refractory Multiple Myeloma

In Multiple myeloma (MM) no treatment has a curative potential and even complete response to novel agents such as proteasome inhibitors (PIs) and immunomodulatory agents (IMiDs) are followed by relapse over time. Next generation sequencing (NGS) has showed how MM at diagnosis is defined by several somatic mutations, but only few drivers, even fewer "druggable" mutations, and many found at a subclonal level. At relapse, targeted studies have shown occasional mutations in drug target genes but the genomic and transcriptomic determinants of chemoresistance in MM remains elusive. We selected 42 MM patients refractory to both lenalidomide and PIs. Whole exome sequencing was performed in 40 of them, and RNAseq in 27. Clinical annotation was available for all patients. Standard analysis pipelines where applied to analyze mutations, copy number alterations (CNAs), mutational signatures, gene expression and expressed mutations. Patients received a median of 3 lines of treatment, with median overall survival of 14.6 months from sampling. We found a median of 77.5 mutations per patient, which is more than what reported at diagnosis (Bolli et al, Nature Communications 2014;5:2997). 100% of samples showed evidence of subclonality, and 37% of them exhibited a higher number of subclonal than clonal variants. Therefore, even at this advanced stage the MM genome is evolving and is composed of different subclones that may display different chemosensitivity. The mutational landscape was also different. TP53 mutations were the second most common after KRAS (20% and 17.5%, respectively). Interestingly TP53 mutations all clustered in patients receiving bortezomib as the last line of treatment. Only 2 patients showed a CRBN mutation, both subclonal. Combining mutations and CNA analysis, the TP53 pathway was the most frequently inactivated (45% of patients). Altogether, mutations or deletions of genes in the CRBN E3 ubiquitin ligase complex were found in 32.5% of patients, while proteasomal subunit genes were infrequently hit. Refractory cases were also uniquely characterized by a novel signature linked to exposure to alkylating agents, whose activity was more pronounced after high-dose melphalan suggesting a mutagenic effect of the drug on residual cells at the time of transplant. Whether this has any pathogenetic role on the disease course remains to be elucidated. RNAseq analysis did not show any influence of treatment or mutational data on the clustering of samples, which was mainly influenced by karyotypic events. The main cluster was composed by non-hyperdiploid patients with both amp(1q) and del(13): these showed CCND2 and MCL1 upregulation, the latter representing a marker of venetoclax resistance and novel target of experimental treatments. Only 26.3% of mutations were expressed, and this correlated with the clonality level of the mutation. However, most mutations in driver genes were expressed, with the notable exception of those causing nonsense mediated decay. Overall, classical high-risk features or CRBN pathway mutations were found in 65% of the cohort. However, only amp(1q) predicted survival in our cohort. The lack of prognostic value of high-risk lesions is likely explained by a higher prevalence of such features in double-refractory stages. Our data suggest that gene mutation is not a preferred mode of evolution of drug resistance in MM. Chemoresistance of the bulk tumor population is likely attained though differential, yet converging evolution of different subclones that are overall highly variable from patient to patient and within the same patient. Disclosures Kastritis: Prothena: Honoraria; Genesis: Honoraria; Amgen: Honoraria, Research Funding; Janssen: Honoraria, Research Funding; Takeda: Honoraria; Pfizer: Honoraria. Dimopoulos:Sanofi Oncology: Research Funding. Cavo:bms: Honoraria, Membership on an entity's Board of Directors or advisory committees, Speakers Bureau; janssen: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Other: travel accommodations, Speakers Bureau; sanofi: Honoraria, Membership on an entity's Board of Directors or advisory committees, Speakers Bureau; novartis: Honoraria; takeda: Honoraria, Membership on an entity's Board of Directors or advisory committees, Speakers Bureau; celgene: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Other: travel accommodations, Speakers Bureau; AbbVie: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees; amgen: Consultancy, Honoraria, Membership on an entity's Board of Directors or advisory committees, Speakers Bureau. Corradini:Janssen: Honoraria, Other: Travel Costs; Takeda: Honoraria, Other: Travel Costs; Jazz Pharmaceutics: Honoraria; Gilead: Honoraria, Other: Travel Costs; Daiichi Sankyo: Honoraria; Celgene: Honoraria, Other: Travel Costs; Amgen: Honoraria; AbbVie: Consultancy, Honoraria, Other: Travel Costs; KiowaKirin: Honoraria; Kite: Honoraria; BMS: Other: Travel Costs; Sanofi: Honoraria; Servier: Honoraria; Roche: Honoraria; Novartis: Honoraria, Other: Travel Costs. Bolli:Celgene: Honoraria; Novartis: Honoraria; Gilead: Other: travel expenses.

Novel Molecular Pathogenesis and Therapeutic Target in Acute Erythroid Leukemia

Acute erythroid leukemia (AEL) is a unique subtype of acute myeloid leukemia (AML) characterized by erythroid predominance and dysplasia, which is morphologically classified into two subtypes: pure erythroid and myeloid/erythroid leukemias (PEL and MEL). Clinically, AEL exhibits aggressive and rapid clinical course, which is considered to be linked to frequent TP53 mutations and complex karyotypes found in this subtype. Although multiple sequencing studies were conducted to explain such phenotypic variety and to search promising therapeutic targets, neither the diagnostic discrimination of PEL and MEL subtypes nor the establishment of standard therapeutic strategy has fully been successful. To clarify genetic characteristics and to identify novel therapeutic targets, we comprehensively characterized the mutation and copy number alteration (CNA) profiles and assessed differential impacts of the frequent driver genetic abnormalities on disease phenotypes and clinical outcomes. We initially analyzed 105 cases with AEL (PEL; n=10, MEL; n=69, and unspecified; n=26). Among these, 9 cases were serially analyzed at 2-6 time points. All samples were obtained according to protocols approved by the ethics board of each participating institution. As control, 907 cases with non-erythroid AML (NEL) distinguished by the low fraction (&amp;lt;50%) of erythroid cells in bone marrow were enrolled. Publicly available information on clinical and genetic findings was collected from previous publications, including The Cancer Genome Atlas and BeatAML datasets. In total, 1,012 cases with AML were enrolled in this study. The genetic profile in PEL and MEL were also compared with that of NEL. Finally, we recognized biological significance of the molecular defects by functional analysis of xenograft model, and showed novel therapeutic possibilities targeting pathogenic pathways. Most frequently observed were mutations in TP53 (39%), STAG2 (20%), and NPM1 (15%), and KMT2A-partial tandem duplication (PTD) (18%) in AEL. On the basis of mutational profiles, consensus clustering divided AEL into 4 subgroups. Group A (n=41, 39%) was defined by TP53 mutations into which all PELs were categorized. In CNA analysis, amplification in 9p, 19p, and 21q were significantly more enriched in PEL than NEL. Intriguingly, amplification of 9p, 19p, and 21q lesions commonly included JAK2, EPOR, and ERG loci, respectively. These genes were affected by amplification or mutations more frequently in AEL than NEL (Figure 1). Group B (n=19, 18%) was enriched for mutations in STAG2, which was significantly more frequent in MEL than in PEL (P&amp;lt;0.05). Compared to STAG2-mutated cases in NEL, those in MEL showed significantly higher frequency of KMT2A-PTD. Group C (n=16, 15%) was defined by NPM1 mutations frequently accompanied by PTPN11 mutations (50% in Group C). This combination was also distinct in MEL and was significantly more frequent than in NPM1-mutated NEL cohort (P&amp;lt;0.05). In Group D (n=29, 28%) without any TP53, STAG2, or NPM1 mutations (triple-negative), RUNX1, ASXL1, BCOR, PHF6, and U2AF1 mutations were enriched compared to triple-negative NEL. In addition, novel loss-of-function USP9X mutations were exclusively identified in Group D AEL (P&amp;lt;0.05, compared to triple-negative NEL), which has been implicated in erythroid leukemogenesis through the TGF beta pathway. Serial sample sequencing analysis revealed that JAK2 mutation positive myeloproliferative neoplasms evolved into PEL. Clinically, patients in Group A showed shorter overall survival time compared to those in any other group. To further search for a novel therapeutic target to the most severe type of AEL, JAK2 inhibitor was administered to xenograft models (from 2 cases) having JAK2 and EPOR amplifications and TP53 mutations. Engrafted cells showed in vitro activation of STAT5 and proliferative capacity enhanced with erythropoietin. In vivo analysis showed that inhibition of JAK/STAT pathway significantly suppressed cell growth and prolonged overall survival. In summary, our findings suggest that AEL is classified into 4 major subgroups having unique genetic and clinical features. PEL and MEL are genetically distinct subtypes, which was also highlighted by comparison to NEL. Frequent involvement of EPOR/JAK/STAT pathway suggests therapeutic indication of JAK inhibition for AEL especially in the most aggressive type of AEL cases with TP53 mutations. Disclosures Yoda: Chordia Therapeutics Inc.: Research Funding. Miyazaki:Kyowa-Kirin: Honoraria; Dainippon-Sumitomo: Honoraria; Nippon-Shinyaku: Honoraria; Chugai: Research Funding; Otsuka: Honoraria; Novartis: Honoraria. Takaori-Kondo:Pfizer: Honoraria; Chugai: Research Funding; Janssen: Honoraria; Kyowa Kirin: Research Funding; Takeda: Research Funding; Ono: Research Funding; Bristol-Myers Squibb: Honoraria, Research Funding; Celgene: Honoraria, Research Funding; Novartis: Honoraria. Nakagawa:Sumitomo Dainippon Pharma Co., Ltd.: Research Funding. Usuki:Astellas Pharma Inc: Research Funding, Speakers Bureau; Daiichi Sankyo Co., Ltd.: Research Funding, Speakers Bureau. Heuser:Bayer Pharma AG, Berlin: Research Funding; Synimmune: Research Funding. Maciejewski:Novartis: Consultancy; Alexion: Consultancy. Ogawa:Qiagen Corporation: Patents &amp; Royalties; RegCell Corporation: Equity Ownership; Kan Research Laboratory, Inc.: Consultancy; ChordiaTherapeutics, Inc.: Consultancy, Equity Ownership; Asahi Genomics: Equity Ownership; Dainippon-Sumitomo Pharmaceutical, Inc.: Research Funding.

Integrated Analysis of Copy-Number Alterations and Gene Mutations in 2,000 Patients with Myeloid Neoplasms

Integrated Analysis of Copy-Number Alterations and Gene Mutations in 2,000 Patients with Myeloid Neoplasms

PDTM-16. PEDIATRIC HIGH GRADE GLIOMA RESOURCES FROM THE CHILDREN’S BRAIN TUMOR TISSUE CONSORTIUM (CBTTC) AND PEDIATRIC BRAIN TUMOR ATLAS (PBTA)

Abstract BACKGROUND Pediatric high grade glioma (pHGG) remains a fatal disease. Access to richly annotated biospecimens and patient derived tumor models will accelerate pHGG research and support translation of research discoveries. This work describes the pediatric high grade glioma set of the Children’s Brain Tumor Tissue Consortium (CBTTC) from the first release of the Pediatric Brain Tumor Atlas (PBTA). METHODS pHGG tumors with associated clinical data and imaging were prospectively collected through the CBTTC and analyzed as the Pediatric Brain Tumor Atlas (PBTA) with processed genomic data deposited into PedcBioPortal for broad access and visualization. Matched tumor was cultured to create high grade glioma cell lines analyzed by targeted and WGS and RNA-seq. A tissue microarray (TMA) of primary pHGG tumors was also created. RESULTS The pHGG set includes 87 collection events (73 patients, 60% at diagnosis, median age of 9 yrs, 55% female, 46% hemispheric). Operative reports, pathology reports and histology images are available for nearly all cases. Pre- and post-operative MRI images and reports are also available for a subset. Tumor WGS/RNAseq is available for 70 subjects. Analysis of somatic mutations and copy number alterations of known glioma genes were of expected distribution (36% H3.3, 47% TP53, 24% ATRX and 7% BRAFV600E variants). In our panel of pHGG, six patients (8 tumors) harbored germline mismatch repair mutations with tumor hyper-mutation. A pHGG TMA (n=77), includes 36 patient tumors with matched sequencing. At least one established glioma cell line was generated from 23 patients (32%). Unique reagents include those derived from a H3.3 G34R glioma and from tumors with mismatch repair deficiency. CONCLUSION The CBTTC and PBTA have created an openly available integrated resource of over 2,000 tumors, including a rich set of pHGG primary tumors, corresponding cell lines and archival fixed tissue to advance translational research for pHGG.

Comprehensive Genomic Review of TCGA Head and Neck Squamous Cell Carcinomas (HNSCC).

The aim of this present study was to comprehensively describe somatic DNA alterations and transcriptional alterations in the last extension of the HNSCC subsets in TCGA, encompassing a total of 528 tumours. In order to achieve this goal, transcriptional analysis, functional enrichment assays, survival analysis, somatic copy number alteration analysis and somatic alteration analysis were carried out. A total of 3491 deregulated genes were found in HNSCC patients, and the functional analysis carried out determined that tissue development and cell differentiation were the most relevant signalling pathways in upregulated and downregulated genes, respectively. Somatic copy number alteration analysis showed a “top five” altered HNSCC genes: CDKN2A (deleted in 32.03% of patients), CDKN2B (deleted in 28.34% of patients), PPFIA1 (amplified in 26.02% of patients), FADD (amplified in 25.63% of patients) and ANO1 (amplified in 25.44% of patients). Somatic mutations analysis revealed TP53 mutation in 72% of the tumour samples followed by TTN (39%), FAT1 (23%) and MUC16 (19%). Another interesting result is the mutual exclusivity pattern that was discovered between the TP53 and PIK3CA mutations, and the co-occurrence of CDKN2A with the TP53 and FAT1 alterations. On analysis to relate differential expression genes and somatic copy number alterations, some genes were overexpressed and amplified, for example, FOXL2, but other deleted genes also showed overexpression, such as CDKN2A. Survival analysis revealed that overexpression of some oncogenes, such as EGFR, CDK6 or CDK4 were associated with poorer prognosis tumours. These new findings help us to develop new therapies and programs for the prevention of HNSCC.

Identification of driver genes regulating immune cell infiltration in cervical cancer by multiple omics integration

Cervical cancer (CC) is one of the most common cancers in women. However, copy number alteration (CNA)-driven dysregulated genes and their functions in CC are still rarely investigated. In the present study, we conducted integrative analysis of CNA and gene expression data from The Cancer Genome Atlas (TCGA) cervical cancer to identify dysregulated genes triggered by CNAs. The integration of copy number status and RNA expression revealed 763 amplified and 1,391 deleted genes significantly dysregulated by the CNAs (P-value < 1e-8). Among these CNA genes, five driver genes, including PI3KCA, PI3KCB, DVL3, WWTR1, and ERBB2, exhibited a strong association with immune cell infiltration, suggesting that the pathways that they participate in may be involved in regulating immune cell infiltration. Moreover, we also observed that the genes of immunotherapeutic targets were abundantly expressed in the wild-type samples, suggesting that immunotherapy based on these immunotherapeutic targets may be applied to wild-type samples. In addition, the two CNA driver genes, DVL3 and ERBB2, might be sensitive and resistant biomarkers for examining the tumor’s response to chemoradiotherapy, respectively. Particularly, the expression of ERBB2 was also observed to be higher in responders of chemotherapy than non-responders. Furthermore, a subset of CNA genes was identified to predict the prognosis of cervical cancer. In summary, our systematic data analysis of these CNA genes not only improved our understanding of the veiled mechanism behind immune cell infiltration, but also provided the potential clinical application of these CNA genes in cervical cancer.

VCF2CNA: A tool for efficiently detecting copy-number alterations in VCF genotype data and tumor purity

VCF2CNA is a tool (Linux commandline or web-interface) for copy-number alteration (CNA) analysis and tumor purity estimation of paired tumor-normal VCF variant file formats. It operates on whole genome and whole exome datasets. To benchmark its performance, we applied it to 46 adult glioblastoma and 146 pediatric neuroblastoma samples sequenced by Illumina and Complete Genomics (CGI) platforms respectively. VCF2CNA was highly consistent with a state-of-the-art algorithm using raw sequencing data (mean F1-score = 0.994) in high-quality whole genome glioblastoma samples and was robust to uneven coverage introduced by library artifacts. In the whole genome neuroblastoma set, VCF2CNA identified MYCN high-level amplifications in 31 of 32 clinically validated samples compared to 15 found by CGI’s HMM-based CNA model. Moreover, VCF2CNA achieved highly consistent CNA profiles between WGS and WXS platforms (mean F1 score 0.97 on a set of 15 rhabdomyosarcoma samples). In addition, VCF2CNA provides accurate tumor purity estimates for samples with sufficient CNAs. These results suggest that VCF2CNA is an accurate, efficient and platform-independent tool for CNA and tumor purity analyses without accessing raw sequence data.

Abstract 909: Single-cell transcriptomic profiling of non-small cell lung cancer uncovers inter- and intracell population structure across TCGA lung adenocarcinoma and lung squamous cancer subtypes

Abstract Lung adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC) are the most prevalent types of non-all cell lung cancer (NSCLC), a leading cause of cancer death worldwide. In this study, we analyzed the transcriptomes of ~45,000 single cells (scRNA) from 13 NSCLC patients, including 5 LUAD cases which were collected and profiled at our institution. To correlate genomes and transcriptomes we performed Whole-Genome Sequencing (WGS) on 3 of these 5 LUAD cases. By comparing tumor tissue with matched adjacent non-malignant lung tissue we are able to confidently distinguish 13 cell-type specific clusters that unambiguously match previously characterized lineages. We developed algorithms for the identification of malignant cells derived from tumor tissue through scRNA analysis of copy number alterations and single nucleotide variants (SNV). Joint analysis of WGS and scRNA confirmed an enrichment of tobacco-associated SNVs among malignant cells of the tumor. Stromal cell types demonstrated consistent expression patterns across cases, while malignant cells demonstrated both inter- and intra-tumoral heterogeneity in their expression of signatures related to GPCR signaling, 3’ UTR mediated translational regulation, and cell-cell junction organization. In particular, one case displayed a unique pattern of intra-tumoral heterogeneity, as a subset of malignant cells robustly express a marker of pulmonary neuroendocrine cells, CGRP. Employing immunohistochemistry, the spatial organization of these malignant cells is revealed to be mutually exclusive within the tumor microenvironment and overlapping in expression of clinical markers of small-cell lung cancer. Finally, we deconvolved bulk TCGA LUAD and LUSC gene expression samples and analyzed the relationship between cell type specific gene expression in cell types of the lung and passenger mutation topographies. Our results provide insight into the molecular and clinical correlates of deconvolved NSCLC transcriptomes and provide a novel methodology with which to explore genomic variation at a single cell resolution. Furthermore, our dataset provides a resource for illuminating cancer-cell transcriptional changes and revealing key molecular drivers of tumor-stromal interactions in lung cancer. Citation Format: Kofi E. Gyan, Aditya Deshpande, Shaham Beg, Huasong Tian, Joel Rosiene, Marlon Stoeckius, Peter Smibert, Davide Risso, Juan Miguel Mosquera, Marcin Imielinski. Single-cell transcriptomic profiling of non-small cell lung cancer uncovers inter- and intracell population structure across TCGA lung adenocarcinoma and lung squamous cancer subtypes [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 909.

Abstract 1390: Molecular characterization of circulating tumor cells in cholangiocarcinoma patients: A new tool for treatment management

Abstract BACKGROUND: Cholangiocarcinoma (CCA) is a highly fatal disease mainly treated with standard chemotherapy, albeit with limited efficacy. New therapeutic options are greatly needed, but the use of targeted treatments is often prevented by the impossibility to obtain tissue biopsies for molecular characterization. Here, we propose the use of circulating tumor cells (CTCs) as an alternative source of tumor material to perform molecular characterization for the identification of novel therapeutic targets. MATERIALS AND METHODS: Blood samples (10 ml) from patients with advanced CCA were processed for CTC isolation as follows: -CTC enrichment with Parsortix -identification and single-cell recovery of epithelial CTCs (expressing epithelial markers) and non-conventional CTCs (lacking epithelial and leukocyte markers) using the DEPArray -whole genome amplification and quality check using Ampli1 kit and Ampli1 QC kit -mutational profiling using Ion AmpliSeq Cancer HotSpot Panel v2 and AmpliSeq somatic pipeline for variant calling -copy number alteration (CNA) analysis using Ampli1 LowPass kit, plus unsupervised clustering and frequency alteration analyses. RESULTS: We analyzed 88 single CTCs isolated from 38 blood samples longitudinally collected from 23 patients (12 with intrahepatic, 9 with extrahepatic CCA and 2 with gallbladder cancer). CNA profiles showed a high level of both inter- and intra-patient heterogeneity, with each CTC displaying a unique profile. Intra-patient heterogeneity was further confirmed by clustering analysis as, in most cases, CTCs from the same patient clustered independently. CTC clustering was also not affected by sampling time (before/during chemotherapy), nor by the anatomical location of primary tumor. Conversely, we observed an enrichment of CTCs derived from patients non-responding to therapy (showing a PD according to RECIST criteria) in 2 of the 4 identified clusters (p=0.00041). By pairwise comparison of CNAs among clusters, we identified 2 regions more frequently altered in one cluster enriched for CTCs from non-responders: 10q22.2 and 3p11.1. The latter encodes, among others, for EPHA3, a targetable gene whose involvement in chemoresistance will be investigated by in vitro studies. Mutational profiling of 19 CTCs (from 6 patients) also confirmed the high intra-patient heterogeneity with most mutations being present in only 1 CTC. This limits the applicability of this approach in patients with few CTCs. Nonetheless, in 1 patient presenting 9 CTCs, we identified 1 mutation in KIT shared by 7/9 CTCs, indicating it as a possible treatment target for this patient. CONCLUSIONS: Our results support the possibility of using CTC molecular characterization to identify both resistance mechanisms and patient-specific targets, thus opening the way for a shift in treatment management of CCA towards an innovative and personalized therapy. Citation Format: Carolina Reduzzi, Marta Vismara, Marco Silvestri, Monica Niger, Rosita Motta, Giorgia Peverelli, Patrizia Miodini, Luigi Celio, Filippo De Braud, Maria G. Daidone, Vera Cappelletti. Molecular characterization of circulating tumor cells in cholangiocarcinoma patients: A new tool for treatment management [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 1390.

Abstract 1670: A hidden markov modeling approach for identifying tumor subclones in next-generation sequencing studies

Abstract Allele-specific copy number alteration (ASCNA) analysis is for identifying copy number abnormalities in tumor cells. Unlike normal cells, tumor cells are heterogeneous with a combination of dominant and minor subclones with distinct copy number profiles. Estimating the clonal proportion and identifying main and subclone genotypes across the genome is important for understanding tumor progression. Several ASCNA tools have recently been developed, but they have been limited to the identification of subclone regions, and not the genotype of the subclones. In this paper, we propose subHMM, a hidden Markov model-based approach that estimates both subclone region as well as region-specific subclone genotype and clonal proportion. We specify a hidden state variable representing the conglomeration of clonal genotypes and subclone status. We propose a two-step algorithm for parameter estimation, where in the first step, a standard hidden Markov model with this conglomerated state space is fit. Then, in the second step, region-specific estimates of the clonal proportions are obtained by maximizing region-specific pseudo-likelihoods. We apply subHMM to a renal cell carcinoma dataset from The Cancer Genome Atlas. In addition, we conduct simulation studies that show the good performance of the proposed approach. Citation Format: Bin Zhu, Hyoyoung Choo-Wosoba, Paul Albert. A hidden markov modeling approach for identifying tumor subclones in next-generation sequencing studies [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 1670.