Aneuploidy in human cancer: new tools and perspectives.

Chromosome copy number imbalances, otherwise known as aneuploidies, are a common but poorly understood feature of cancer. Here, we describe recent advances in both detecting and manipulating aneuploidies that have greatly advanced our ability to study their role in tumorigenesis. In particular, new clustered regularly interspaced short palindromic repeats (CRISPR)-based techniques have been developed that allow the creation of isogenic cell lines with specific chromosomal changes, thereby facilitating experiments in genetically controlled backgrounds to uncover the consequences of aneuploidy. These approaches provide increasing evidence that aneuploidy is a key driver of cancer development and enable the identification of multiple dosage-sensitive genes encoded on aneuploid chromosomes. Consequently, measuring aneuploidy may inform clinical prognosis, while treatment strategies that target aneuploidy could represent a novel method to counter malignant growth.

labelSeg: segment annotation for tumor copy number alteration profiles.

Somatic copy number alterations (SCNAs) are a predominant type of oncogenomic alterations that affect a large proportion of the genome in the majority of cancer samples. Current technologies allow high-throughput measurement of such copy number aberrations, generating results consisting of frequently large sets of SCNA segments. However, the automated annotation and integration of such data are particularly challenging because the measured signals reflect biased, relative copy number ratios. In this study, we introduce labelSeg, an algorithm designed for rapid and accurate annotation of CNA segments, with the aim of enhancing the interpretation of tumor SCNA profiles. Leveraging density-based clustering and exploiting the length-amplitude relationships of SCNA, our algorithm proficiently identifies distinct relative copy number states from individual segment profiles. Its compatibility with most CNA measurement platforms makes it suitable for large-scale integrative data analysis. We confirmed its performance on both simulated and sample-derived data from The Cancer Genome Atlas reference dataset, and we demonstrated its utility in integrating heterogeneous segment profiles from different data sources and measurement platforms. Our comparative and integrative analysis revealed common SCNA patterns in cancer and protein-coding genes with a strong correlation between SCNA and messenger RNA expression, promoting the investigation into the role of SCNA in cancer development.

Refinement of computational identification of somatic copy number alterations using DNA methylation microarrays illustrated in cancers of unknown primary.

High-throughput genomic technologies are increasingly used in personalized cancer medicine. However, computational tools to maximize the use of scarce tissues combining distinct molecular layers are needed. Here we present a refined strategy, based on the R-package 'conumee', to better predict somatic copy number alterations (SCNA) from deoxyribonucleic acid (DNA) methylation arrays. Our approach, termed hereafter as 'conumee-KCN', improves SCNA prediction by incorporating tumor purity and dynamic thresholding. We trained our algorithm using paired DNA methylation and SNP Array 6.0 data from The Cancer Genome Atlas samples and confirmed its performance in cancer cell lines. Most importantly, the application of our approach in cancers of unknown primary identified amplified potentially actionable targets that were experimentally validated by Fluorescence in situ hybridization and immunostaining, reaching 100% specificity and 93.3% sensitivity.

CDKN2A somatic copy number amplification in normal tissues surrounding gastric carcinoma reduces cancer metastasis risk in droplet digital PCR analysis.

BACKGROUND: The CDKN2A gene is frequently affected by somatic copy number variations (SCNVs, including deletions and amplifications [SCNdel and SCNamp]) in the cancer genome. Using surgical gastric margin tissue samples (SMs) as the diploid reference in SCNV analysis via CDKN2A/P16-specific real-time PCR (P16-Light), we previously reported that the CDKN2A SCNdel was associated with a high risk of metastasis of gastric carcinoma (GC). However, the status of CDKN2A SCNVs in SMs and their clinical significance have not been reported. METHODS: Peripheral white blood cell (WBC) and frozen GC and SM tissue samples were collected from patients (n = 80). Droplet digital PCR (ddPCR) was used to determine the copy number (CN) of the CDKN2A gene in tissue samples using paired WBCs as the diploid reference. RESULTS: A novel P16-ddPCR system was initially established with a minimal proportion (or limit, 10%) of the detection of CDKN2A CN alterations. While CDKN2A SCNamp events were detected in both SMs and GCs, fewer CDKN2A SCNdel events were detected in SMs than in GCs (15.0% vs. 41.3%, P = 4.77E-04). Notably, significantly more SCNamp and fewer SCNdel of the CDKN2A gene were detected in SMs from GC patients without metastasis than in those from patients with lymph node metastasis by P16-ddPCR (P = 0.023). The status of CDKN2A SCNVs in SM samples was significantly associated with overall survival (P = 0.032). No cancer deaths were observed among the 11 patients with CDKN2A SCNamp. CONCLUSION: CDKN2A SCNVs in SMs identified by P16-ddPCR are prevalent and significantly associated with GC metastasis and overall survival.

SOX2 dosage sustains tumor-promoting inflammation to drive disease aggressiveness by modulating the FOSL2/IL6 axis.

BACKGROUND: Inflammation is undoubtedly a hallmark of cancer development. Its maintenance within tumors and the consequences on disease aggressiveness are insufficiently understood. METHODS: Data of 27 tumor entities (about 5000 samples) were downloaded from the TCGA and GEO databases. Multi-omic analyses were performed on these and in-house data to investigate molecular determinants of tumor aggressiveness. Using molecular loss-of-function data, the mechanistic underpinnings of inflammation-induced tumor aggressiveness were addressed. Patient specimens and in vivo disease models were subsequently used to validate findings. RESULTS: There was significant association between somatic copy number alterations (sCNAs) and tumor aggressiveness. SOX2 amplification was the most important feature among novel and known aggressiveness-associated alterations. Mechanistically, SOX2 regulates a group of genes, in particular the AP1 transcription factor FOSL2, to sustain pro-inflammatory signaling pathways, such as IL6-JAK-STAT3, TNFA and IL17. FOSL2 was found overexpressed in tumor sections of specifically aggressive cancers. In consequence, prolonged inflammation induces immunosuppression and activates cytidine deamination and thus DNA damage as evidenced by related mutational signatures in aggressive tumors. The DNA damage affects tumor suppressor genes such as TP53, which is the most mutated gene in aggressive tumors compared to less aggressive ones (38% vs 14%), thereby releasing cell cycle control. These results were confirmed by analyzing tissues from various tumor types and in vivo studies. CONCLUSION: Our data demonstrate the implication of SOX2 in promoting DNA damage and genome instability by sustaining inflammation via FOSL2/IL6, resulting in tumor aggressiveness.

Somatic Copy Number Alterations in Colorectal Cancer Lead to a Differentially Expressed ceRNA Network (ceRNet).

Colorectal cancer (CRC) represents the second deadliest malignancy worldwide. Around 75% of CRC patients exhibit high levels of chromosome instability that result in the accumulation of somatic copy number alterations. These alterations are associated with the amplification of oncogenes and deletion of tumor-ppressor genes and contribute to the tumoral phenotype in different malignancies. Even though this relationship is well known, much remains to be investigated regarding the effect of said alterations in long non-coding RNAs (lncRNAs) and, in turn, the impact these alterations have on the tumor phenotype. The present study aimed to evaluate the role of differentially expressed lncRNAs coded in regions with copy number alterations in colorectal cancer patient samples. We downloaded RNA-seq files of the Colorectal Adenocarcinoma Project from the The Cancer Genome Atlas (TCGA) repository (285 sequenced tumor tissues and 41 non-tumor tissues), evaluated differential expression, and mapped them over genome sequencing data with regions presenting copy number alterations. We obtained 78 differentially expressed (LFC > 1|< -1, padj < 0.05) lncRNAs, 410 miRNAs, and 5028 mRNAs and constructed a competing endogenous RNA (ceRNA) network, predicting significant lncRNA-miRNA-mRNA interactions. Said network consisted of 30 lncRNAs, 19 miRNAs, and 77 mRNAs. To understand the role that our ceRNA network played, we performed KEGG and GO analysis and found several oncogenic and anti-oncogenic processes enriched by the molecular players in our network. Finally, to evaluate the clinical relevance of the lncRNA expression, we performed survival analysis and found that C5orf64, HOTAIR, and RRN3P3 correlated with overall patient survival. Our results showed that lncRNAs coded in regions affected by SCNAs form a complex gene regulatory network in CCR.

CNA2Subpathway: identification of dysregulated subpathway driven by copy number alterations in cancer.

Biological pathways reflect the key cellular mechanisms that dictate disease states, drug response and altered cellular function. The local areas of pathways are defined as subpathways (SPs), whose dysfunction has been reported to be associated with the occurrence and development of cancer. With the development of high-throughput sequencing technology, identifying dysfunctional SPs by using multi-omics data has become possible. Moreover, the SPs are not isolated in the biological system but interact with each other. Here, we propose a network-based calculated method, CNA2Subpathway, to identify dysfunctional SPs is driven by somatic copy number alterations (CNAs) in cancer through integrating pathway topology information, multi-omics data and SP crosstalk. This provides a novel way of SP analysis by using the SP interactions in the system biological level. Using data sets from breast cancer and head and neck cancer, we validate the effectiveness of CNA2Subpathway in identifying cancer-relevant SPs driven by the somatic CNAs, which are also shown to be associated with cancer immune and prognosis of patients. We further compare our results with five pathway or SP analysis methods based on CNA and gene expression data without considering SP crosstalk. With these analyses, we show that CNA2Subpathway could help to uncover dysfunctional SPs underlying cancer via the use of SP crosstalk. CNA2Subpathway is developed as an R-based tool, which is freely available on GitHub (https://github.com/hanjunwei-lab/CNA2Subpathway).

Simultaneous Copy Number Alteration and Single-Nucleotide Variation Analysis in Matched Aqueous Humor and Tumor Samples in Children with Retinoblastoma.

Retinoblastoma (RB) is a childhood cancer that forms in the developing retina of young children; this tumor cannot be biopsied due to the risk of provoking extraocular tumor spread, which dramatically alters the treatment and survival of the patient. Recently, aqueous humor (AH), the clear fluid in the anterior chamber of the eye, has been developed as an organ-specific liquid biopsy for investigation of in vivo tumor-derived information found in the cell-free DNA (cfDNA) of the biofluid. However, identifying somatic genomic alterations, including both somatic copy number alterations (SCNAs) and single nucleotide variations (SNVs) of the RB1 gene, typically requires either: (1) two distinct experimental protocols-low-pass whole genome sequencing for SCNAs and targeted sequencing for SNVs-or (2) expensive deep whole genome or exome sequencing. To save time and cost, we applied a one-step targeted sequencing method to identify both SCNAs and RB1 SNVs in children with RB. High concordance (median = 96.2%) was observed in comparing SCNA calls derived from targeted sequencing to the traditional low-pass whole genome sequencing method. We further applied this method to investigate the degree of concordance of genomic alterations between paired tumor and AH samples from 11 RB eyes. We found 11/11 AH samples (100%) had SCNAs, and 10 of them (90.1%) with recurrent RB-SCNAs, while only nine out of 11 tumor samples (81.8%) had positive RB-SCNA signatures in both low-pass and targeted methods. Eight out of the nine (88.9%) detected SNVs were shared between AH and tumor samples. Ultimately, 11/11 cases have somatic alterations identified, including nine RB1 SNVs and 10 recurrent RB-SCNAs with four focal RB1 deletions and one MYCN gain. The results presented show the feasibility of utilizing one sequencing approach to obtain SCNA and targeted SNV data to capture a broad genomic scope of RB disease, which may ultimately expedite clinical intervention and be less expensive than other methods.

Low-cost and clinically applicable copy number profiling using repeat DNA.

BACKGROUND: Somatic copy number alterations (SCNAs) are an important class of genomic alteration in cancer. They are frequently observed in cancer samples, with studies showing that, on average, SCNAs affect 34% of a cancer cell's genome. Furthermore, SCNAs have been shown to be major drivers of tumour development and have been associated with response to therapy and prognosis. Large-scale cancer genome studies suggest that tumours are driven by somatic copy number alterations (SCNAs) or single-nucleotide variants (SNVs). Despite the frequency of SCNAs and their clinical relevance, the use of genomics assays in the clinic is biased towards targeted gene panels, which identify SNVs but provide limited scope to detect SCNAs throughout the genome. There is a need for a comparably low-cost and simple method for high-resolution SCNA profiling. RESULTS: We present conliga, a fully probabilistic method that infers SCNA profiles from a low-cost, simple, and clinically-relevant assay (FAST-SeqS). When applied to 11 high-purity oesophageal adenocarcinoma samples, we obtain good agreement (Spearman's rank correlation coefficient, rs=0.94) between conliga's inferred SCNA profiles using FAST-SeqS data (approximately £14 per sample) and those inferred by ASCAT using high-coverage WGS (gold-standard). We find that conliga outperforms CNVkit (rs=0.89), also applied to FAST-SeqS data, and is comparable to QDNAseq (rs=0.96) applied to low-coverage WGS, which is approximately four-fold more expensive, more laborious and less clinically-relevant. By performing an in silico dilution series experiment, we find that conliga is particularly suited to detecting SCNAs in low tumour purity samples. At two million reads per sample, conliga is able to detect SCNAs in all nine samples at 3% tumour purity and as low as 0.5% purity in one sample. Crucially, we show that conliga's hidden state information can be used to decide when a sample is abnormal or normal, whereas CNVkit and QDNAseq cannot provide this critical information. CONCLUSIONS: We show that conliga provides high-resolution SCNA profiles using a convenient, low-cost assay. We believe conliga makes FAST-SeqS a more clinically valuable assay as well as a useful research tool, enabling inexpensive and fast copy number profiling of pre-malignant and cancer samples.

Genomic Analysis Made It Possible to Identify Gene-Driver Alterations Covering the Time Window between Diagnosis of Neuroblastoma 4S and the Progression to Stage 4.

Neuroblastoma (NB) is a tumor of the developing sympathetic nervous system. Despite recent advances in understanding the complexity of NB, the mechanisms that determine its regression or progression are still largely unknown. Stage 4S NB is characterized by a favorable course of disease and often by spontaneous regression, while progression to true stage 4 is a very rare event. Here, we focused on genomic analysis of an NB case that progressed from stage 4S to stage 4 with a very poor outcome. Array-comparative genomic hybridization (a-CGH) on tumor-tissue DNA, and whole-exome sequencing (WES) on exosomes DNA derived from plasma collected at the onset and at the tumor progression, pointed out relevant genetic changes that can explain this clinical worsening. The combination of a-CGH and WES data allowed for the identification iof somatic copy number aberrations and single-nucleotide variants in genes known to be responsible for aggressive NB. KLRB1, MAPK3 and FANCA genes, which were lost at the time of progression, were studied for their possible role in this event by analyzing in silico the impact of their expression on the outcome of 786 NB patients.

Somatic DNA copy number alterations in non-dysplastic Barrett's esophagus.

BACKGROUND: The purpose of this study was to analyze non-dysplastic Barrett's esophagus (NDBE) biopsy tissue and compare the rate of somatic DNA copy number alterations (CNAs) in patients who subsequently progressed to high-grade dysplasia (HGD) or esophageal adenocarcinoma (EAC) to those patients who did not. METHODS: A retrospectively collected database of Barrett's esophagus (BE) patients spanning a 16-year period was queried. Patients who progressed from NDBE to HGD or EAC were identified and compared to patients who did not. Initial biopsy specimens were microdissected and extracted DNA underwent Multiplex Ligation-dependent Probe Amplification (MLPA) for CNAs. Comparisons between progressors and non-progressors were made with Fisher's exact and two-sample t tests. Logistic regression assessed factors associated with progression. RESULTS: Of the 2459 patients in the BE database, 36 patients progressed from NDBE to either HGD or EAC. There were eight progressors who had biopsy specimens with adequate DNA for analysis. The progressor and non-progressor cohort had similar demographic information and medical history. The progressor group trended towards being older at diagnosis (72 ± 10 vs. 64 ± 13 years, p = 0.097) and fewer progressors reported reflux symptoms (50 vs. 94.7%, p < 0.001). Progressor specimens had more overall CNAs (75% vs. 33.6%, p = 0.026). On univariable analysis, there was an association between progression and absence of GERD symptoms (OR 16.54 [3.42-80.03], p = 0.001), any CNA (OR 5.10 [1.18-23.30], p = 0.035), and CNA in GATA6 or ERBB2 (OR 6.72 [1.18-38.22], p = 0.032). CONCLUSIONS: Patients who progressed from NDBE to HGD or EAC were older at first diagnosis of BE and fewer of the progressors reported symptoms of reflux when compared to non-progressors. Progression was associated with the presence of any CNA and specific CNAs in GATA6 or ERBB2.

Variability in retinoblastoma genome stability is driven by age and not heritability.

Retinoblastoma (RB) is a childhood intraocular cancer initiated by biallelic inactivation of the RB tumor suppressor gene (RB1-/- ). RB can be hereditary (germline RB1 pathogenic allele is present) or non-hereditary. Somatic copy number alterations (SCNAs) contribute to subsequent tumorigenesis. Previous studies of only enucleated RB eyes have reported associations between heritability status and the prevalence of SCNAs. Herein, we use an aqueous humor (AH) liquid biopsy to investigate RB genomic profiles in the context of germline RB1 status, age, and International Intraocular Retinoblastoma Classification (IIRC) clinical grouping for both enucleated and salvaged eyes. Between 2014 and 2019, AH was sampled from a total of 54 eyes of 50 patients. Germline RB1 status was determined from clinical blood testing, and cell-free DNA from AH was analyzed for SCNAs. Of the 50 patients, 23 (46.0%; 27 eyes) had hereditary RB, and 27 (54.0%, 27 eyes) had non-hereditary RB. Median age at diagnosis was comparable between hereditary (13 ± 10 months) and non-hereditary (13 ± 8 months) eyes (P = 0.818). There was no significant difference in the prevalence or number of SCNAs based on (1) hereditary status (P > 0.56) or (2) IIRC grouping (P > 0.47). There was, however, a significant correlation between patient age at diagnosis, and (1) number of total SCNAs (r[52] = 0.672, P < 0.00001) and (2) number of highly-recurrent RB SCNAs (r[52] = 0.616, P < 0.00001). This evidence does not support the theory that specific molecular or genomic subtypes exist between hereditary and non-hereditary RB; rather, the prevalence of genomic alterations in RB eyes is strongly related to patient age at diagnosis.

Incorporating genomic signatures into surgical and medical decision-making for elderly glioblastoma patients.

Glioblastoma (GBM) is the most common type of malignant primary brain tumor in adults. It is a uniformly fatal disease (median overall survival 16 months) even with aggressive resection and an adjuvant temozolomide-based chemoradiation regimen. Age remains an independent risk factor for a poor prognosis. Several factors contribute to the dismal outcomes in the elderly population with GBM, including poor baseline health status, differences in underlying genomic alterations, and variability in the surgical and medical management of this subpopulation. The latter arises from a lack of adequate representation of elderly patients in clinical trials, resulting in limited data on the response of this subpopulation to standard treatment. Results from retrospective and some prospective studies have indicated that resection of only contrast-enhancing lesions and administration of hypofractionated radiotherapy in combination with temozolomide are effective strategies for optimizing survival while maintaining baseline quality of life in elderly GBM patients; however, survival remains dismal relative to that in a younger cohort. Here, the authors present historical context for the current strategies used for the multimodal management (surgical and medical) of elderly patients with GBM. Furthermore, they provide insights into elderly GBM patient-specific genomic signatures such as isocitrate dehydrogenase 1/2 (IDH1/2) wildtype status, telomerase reverse transcriptase promoter (TERTp) mutations, and somatic copy number alterations including CDK4/MDM2 coamplification, which are becoming better understood and could be utilized in a clinical trial design and patient stratification to guide the development of more effective adjuvant therapies specifically for elderly GBM patients.

Interaction analysis between germline susceptibility loci and somatic alterations in lung cancer.

Emerging evidence indicates that germline variations may interact with somatic events in carcinogenesis. However, the germline-somatic interaction in lung cancer remains largely unknown. We investigated whether lung cancer driver genes (CDGs) were more likely to locate within cancer susceptibility regions. Pathway analysis was performed to identify common pathways underlying CDGs and cancer susceptibility genes (CSGs). Next, we analyzed the associations between lung cancer risk SNPs and somatic alterations, including mutations and copy number alterations, in the level of genes, pathways, and overall burden of alterations. Enrichment analysis showed that lung CDGs are more likely to locate within cancer susceptibility regions (p = 8.40 × 10-3 ). Both of lung CSGs and CDGs showed significant enrichment in pathways such as cell cycle and p53 signaling pathway. Gene-based analysis showed that rs36600 (22q12.2) was associated with somatic mutations within ARID1A (OR = 2.45, 95%CI: 1.47-4.08, p = 5.78 × 10-4 ). Pathway-based analysis of somatic truncation mutations identified rs2395185 and rs3817963 at 6p22.1 was associated with cell cycle pathway (OR = 1.56, p = 3.61 × 10-4 for rs2395185; OR = 1.58, p = 4.15 × 10-4 for rs3817963), and rs3817963 was also associated with MAPK signaling pathway (OR = 1.54, p = 8.58 × 10-4 ). Further analysis associated rs2395185 at 6p22.1 (HLA class II genes) with increased APOBEC3A expression (p = 9.50 × 10-3 ) and elevated APOBEC mutagenesis (p = 3.58 × 10-3 ). These results indicate germline-somatic interactions in lung tumorigenesis, and help to uncover the molecular mechanisms underlying lung cancer risk SNPs.

Comparative analysis of methods for identifying somatic copy number alterations from deep sequencing data.

Somatic copy-number alterations (SCNAs) are an important type of structural variation affecting tumor pathogenesis. Accurate detection of genomic regions with SCNAs is crucial for cancer genomics as these regions contain likely drivers of cancer development. Deep sequencing technology provides single-nucleotide resolution genomic data and is considered one of the best measurement technologies to detect SCNAs. Although several algorithms have been developed to detect SCNAs from whole-genome and whole-exome sequencing data, their relative performance has not been studied. Here, we have compared ten SCNA detection algorithms in both simulated and primary tumor deep sequencing data. In addition, we have evaluated the applicability of exome sequencing data for SCNA detection. Our results show that (i) clear differences exist in sensitivity and specificity between the algorithms, (ii) SCNA detection algorithms are able to identify most of the complex chromosomal alterations and (iii) exome sequencing data are suitable for SCNA detection.

Characterization of genetic aberrations in a single case of metastatic thymic adenocarcinoma.

BACKGROUND: Thymic adenocarcinoma is an extremely rare subtype of thymic epithelial tumors. Due to its rarity, there is currently no sequencing approach for thymic adenocarcinoma. METHODS: We performed whole exome and transcriptome sequencing on a case of thymic adenocarcinoma and performed subsequent validation using Sanger sequencing. RESULTS: The case of thymic adenocarcinoma showed aggressive behaviors with systemic bone metastases. We identified a high incidence of genetic aberrations, which included somatic mutations in RNASEL, PEG10, TNFSF15, TP53, TGFB2, and FAT1. Copy number analysis revealed a complex chromosomal rearrangement of chromosome 8, which resulted in gene fusion between MCM4 and SNTB1 and dramatic amplification of MYC and NDRG1. Focal deletion was detected at human leukocyte antigen (HLA) class II alleles, which was previously observed in thymic epithelial tumors. We further investigated fusion transcripts using RNA-seq data and found an intergenic splicing event between the CTBS and GNG5 transcript. Finally, enrichment analysis using all the variants represented the immune system dysfunction in thymic adenocarcinoma. CONCLUSION: Thymic adenocarcinoma shows highly malignant characteristics with alterations in several cancer-related genes.

Multi-region and single-cell sequencing reveal variable genomic heterogeneity in rectal cancer.

BACKGROUND: Colorectal cancer is a heterogeneous group of malignancies with complex molecular subtypes. While colon cancer has been widely investigated, studies on rectal cancer are very limited. Here, we performed multi-region whole-exome sequencing and single-cell whole-genome sequencing to examine the genomic intratumor heterogeneity (ITH) of rectal tumors. METHODS: We sequenced nine tumor regions and 88 single cells from two rectal cancer patients with tumors of the same molecular classification and characterized their mutation profiles and somatic copy number alterations (SCNAs) at the multi-region and the single-cell levels. RESULTS: A variable extent of genomic heterogeneity was observed between the two patients, and the degree of ITH increased when analyzed on the single-cell level. We found that major SCNAs were early events in cancer development and inherited steadily. Single-cell sequencing revealed mutations and SCNAs which were hidden in bulk sequencing. In summary, we studied the ITH of rectal cancer at regional and single-cell resolution and demonstrated that variable heterogeneity existed in two patients. The mutational scenarios and SCNA profiles of two patients with treatment naïve from the same molecular subtype are quite different. CONCLUSIONS: Our results suggest each tumor possesses its own architecture, which may result in different diagnosis, prognosis, and drug responses. Remarkable ITH exists in the two patients we have studied, providing a preliminary impression of ITH in rectal cancer.

Genomic copy number alterations of primary and secondary metastasizing pleomorphic adenomas.

AIMS: Metastasizing pleomorphic adenoma (MPA) is a rare tumour, and its mechanism of metastasis still is unknown. To date, there has been no study on MPA genomics. We analysed primary and secondary MPAs with array comparative genomic hybridization to identify somatic copy number alterations and affected genes. METHODS AND RESULTS: Tumour DNA samples from primary (parotid salivary gland) and secondary (scalp skin) MPAs were subjected to array comparative genomic hybridization investigation, and the data were analysed with NEXUS COPY NUMBER DISCOVERY. The primary MPA showed copy number losses affecting 3p22.2p14.3 and 19p13.3p123, and a complex pattern of four different deletions at chromosome 6. The 3p deletion encompassed several genes: CTNNB1, SETD2, BAP1, and PBRM1, among others. The secondary MPA showed a genomic profile similar to that of the primary MPA, with acquisition of additional copy number changes affecting 9p24.3p13.1 (loss), 19q11q13.43 (gain), and 22q11.1q13.33 (gain). CONCLUSION: Our findings indicated a clonal origin of the secondary MPA, as both tumours shared a common profile of genomic copy number alterations. Furthermore, we were able to detect in the primary tumour a specific pattern of copy number alterations that could explain the metastasizing characteristic, whereas the secondary MPA showed a more unbalanced genome.

Subclonal somatic copy number alterations emerge and dominate in recurrent osteosarcoma.

Multiple large-scale tumor genomic profiling efforts have been undertaken in osteosarcoma, however, little is known about the spatial and temporal intratumor heterogeneity and how it may drive treatment resistance. We performed whole-genome sequencing of 37 tumor samples from eight patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. We identified subclonal copy number alterations in all but one patient. We observed that in five patients, a subclonal copy number clone from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clone in 6 out of 7 patients with more than one clone. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy number clones. Our study sheds light on intratumor heterogeneity and the potential drivers of treatment resistance in osteosarcoma.

Genomic analyses of germline and somatic variation in high-grade serous ovarian cancer.

BACKGROUND: High-grade serous ovarian cancers (HGSCs) display a high degree of complex genetic alterations. In this study, we identified germline and somatic genetic alterations in HGSC and their association with relapse-free and overall survival. Using a targeted capture of 557 genes involved in DNA damage response and PI3K/AKT/mTOR pathways, we conducted next-generation sequencing of DNA from matched blood and tumor tissue from 71 HGSC participants. In addition, we performed the OncoScan assay on tumor DNA from 61 participants to examine somatic copy number alterations (SCNA). RESULTS: Approximately one-third of tumors had loss-of-function (LOF) germline (18/71, 25.4%) or somatic (7/71, 9.9%) variants in the DNA homologous recombination repair pathway genes BRCA1, BRCA2, CHEK2, MRE11A, BLM, and PALB2. LOF germline variants also were identified in other Fanconi anemia genes and in MAPK and PI3K/AKT/mTOR pathway genes. Most tumors harbored somatic TP53 variants (65/71, 91.5%). Using the OncoScan assay on tumor DNA from 61 participants, we identified focal homozygous deletions in BRCA1, BRCA2, MAP2K4, PTEN, RB1, SLX4, STK11, CREBBP, and NF1. In total, 38% (27/71) of HGSC patients harbored pathogenic variants in DNA homologous recombination repair genes. For patients with multiple tissues from the primary debulking or from multiple surgeries, the somatic mutations were maintained with few newly acquired point mutations suggesting that tumor evolution was not through somatic mutations. There was a significant association of LOF variants in homologous recombination repair pathway genes and high-amplitude somatic copy number alterations. Using GISTIC analysis, we identified NOTCH3, ZNF536, and PIK3R2 in these regions that were significantly associated with an increase in cancer recurrence and a reduction in overall survival. CONCLUSIONS: From 71 patients with HGCS, we performed targeted germline and tumor sequencing and provided a comprehensive analysis of these 557 genes. We identified germline and somatic genetic alterations including somatic copy number alterations and analyzed their associations with relapse-free and overall survival. This single-site long-term follow-up study provides additional information on genetic alterations related to occurrence and outcome of HGSC. Our findings suggest that targeted treatments based on both variant and SCNA profile potentially could improve relapse-free and overall survival.