Longitudinal cell-free DNA characterization by low-coverage whole-genome sequencing in patients undergoing high-dose radiotherapy.

BACKGROUND AND PURPOSE: Current radiotherapy guidelines rely heavily on imaging-based monitoring. Liquid biopsy monitoring promises to complement imaging by providing frequent systemic information about the tumor. In particular, cell-free DNA (cfDNA) sequencing offers a tumor-agnostic approach, which lends itself to monitoring heterogeneous cohorts of cancer patients. METHODS: We collected plasma cfDNA from oligometastatic patients (OMD) and head-and-neck cancer patients (SCCHN) at six time points before, during, and after radiotherapy, and compared them to the plasma samples of healthy and polymetastatic volunteers. We performed low-pass (on average 7x) whole-genome sequencing on 93 plasma cfDNA samples and correlated copy number alterations and fragment length distributions to clinical and imaging findings. RESULTS: We observed copy number alterations in 4/7 polymetastatic cancer patients, 1/7 OMD and 1/7 SCCHN patients, these patients' imaging showed progression following radiotherapy. Using unsupervised learning, we identified cancer-specific fragment length features that showed a strong correlation with copy number-based tumor fraction estimates. In 4/4 HPV-positive SCCHN patient samples, we detected viral DNA that enabled the monitoring of very low tumor fraction samples. CONCLUSIONS: Our results indicate that an elevated tumor fraction is associated with tumor aggressiveness and systemic tumor spread. This information may be used to adapt treatment strategies. Further, we show that by detecting specific sequences such as viral DNA, the sensitivity of detecting cancer from cell-free DNA sequencing data can be greatly increased.

Concomitant genomic features stratify prognosis to patients with advanced EGFR mutant lung cancer.

This study aimed to explore the clinical significance of genomics features including tumor mutation burden (TMB) and copy number alteration (CNA) for advanced EGFR mutant lung cancer. We retrospectively identified 1378 patients with advanced EGFR mutant lung cancer and next-generation sequencing tests from three cohorts. Multiple co-occurring genomics alternations occurred in a large proportion (97%) of patients with advanced EGFR mutant lung cancers. Both TMB and CNA were predictive biomarkers for these patients. A joint analysis of TMB and CNA found that patients with high TMB and high CNA showed worse responses to EGFR-TKIs and predicted worse outcomes. TMBhighCNAhigh, as a high-risk genomic feature, showed predictive ability in most of the subgroups based on clinical characteristics. These patients had larger numbers of metastatic sites, and higher rates of EGFR copy number amplification, TP53 mutations, and cell-cycle gene alterations, which showed more potential survival gain from combination treatment. Furthermore, a nomogram based on genomic features and clinical features was developed to distinguish prognosis. Genomic features could stratify prognosis and guide clinical treatment for patients with advanced EGFR mutant lung cancer.

NKX2­1 copy number alterations are associated with oncogenic, immunological and prognostic remodeling in non­small cell lung cancer.

NK2 homeobox 1 (NKX2-1) copy number alterations (CNAs) are frequently observed in lung cancer. However, little is known about the complete landscape of focal alterations in NKX2-1 copy number (CN), their clinical significance and their therapeutic implications in non-small cell lung cancer (NSCLC). The correlations between NKX2-1 expression and EGFR driver mutations and programmed death ligand 1 (PD-L1) co-expression were studied using immunohistochemistry and PCR from the tumors of recruited Filipino patients (n=45). Clinical features of NSCLC with NKX2-1 CNAs were resolved at the tumor and clonal levels using the molecular profiles of patients with lung adenocarcinoma and lung squamous cell carcinoma from The Cancer Genome Atlas (n=1,130), and deconvoluted single-cell RNA-seq data from the Bivona project (n=1,654), respectively. Despite a significant and positive correlation between expression and CN (r=0.264; P<0.001), NKX2-1 CNAs exerted a stronger influence on the combined EGFR and PD-L1 status of NSCLC tumors than expression. NKX2-1 CN gain was prognostic of favorable survival (P=0.018) and a better response to targeted therapy. NKX2-1 CN loss predicted a worse survival (P=0.041). Mutational architecture in the Y-chromosome differentiated the two prognostic groups. There were 19,941 synonymous mutations and 1,408 genome-wide CN perturbations associated with NKX2-1 CNAs. Tumors with NKX2-1 CN gain expressed lymphocyte markers more heterogeneously than those with CN loss. Higher expression of tumor-infiltrating lymphocyte gene signatures in CN gain was prognostic of longer disease-free survival (P=0.005). Tumors with NKX2-1 CN gain had higher B-cell (P<0.001) and total T-cell estimates (P=0.003). NKX2-1 CN loss was associated with immunologically colder tumors due to higher M2 macrophage infiltrates (P=0.011) and higher expression of immune checkpoint proteins, CD274 (P=0.025), VTCN1 (P<0.001) and LGALS9 (P=0.002). In conclusion, NKX2-1 CNAs are associated with tumors that exhibit clinically diverse characteristics, and with unique oncogenic, immunological and prognostic signatures.

Mutational landscape of inflammatory breast cancer.

BACKGROUND: Inflammatory breast cancer (IBC) is the most pro-metastatic form of BC. Better understanding of its enigmatic pathophysiology is crucial. We report here the largest whole-exome sequencing (WES) study of clinical IBC samples. METHODS: We retrospectively applied WES to 54 untreated IBC primary tumor samples and matched normal DNA. The comparator samples were 102 stage-matched non-IBC samples from TCGA. We compared the somatic mutational profiles, spectra and signatures, copy number alterations (CNAs), HRD and heterogeneity scores, and frequencies of actionable genomic alterations (AGAs) between IBCs and non-IBCs. The comparisons were adjusted for the molecular subtypes. RESULTS: The number of somatic mutations, TMB, and mutational spectra were not different between IBCs and non-IBCs, and no gene was differentially mutated or showed differential frequency of CNAs. Among the COSMIC signatures, only the age-related signature was more frequent in non-IBCs than in IBCs. We also identified in IBCs two new mutational signatures not associated with any environmental exposure, one of them having been previously related to HIF pathway activation. Overall, the HRD score was not different between both groups, but was higher in TN IBCs than TN non-IBCs. IBCs were less frequently classified as heterogeneous according to heterogeneity H-index than non-IBCs (21% vs 33%), and clonal mutations were more frequent and subclonal mutations less frequent in IBCs. More than 50% of patients with IBC harbored at least one high-level of evidence (LOE) AGA (OncoKB LOE 1-2, ESCAT LOE I-II), similarly to patients with non-IBC. CONCLUSIONS: We provide the largest mutational landscape of IBC. Only a few subtle differences were identified with non-IBCs. The most clinically relevant one was the higher HRD score in TN IBCs than in TN non-IBCs, whereas the most intriguing one was the smaller intratumor heterogeneity of IBCs.

CoT: a transformer-based method for inferring tumor clonal copy number substructure from scDNA-seq data.

Single-cell DNA sequencing (scDNA-seq) has been an effective means to unscramble intra-tumor heterogeneity, while joint inference of tumor clones and their respective copy number profiles remains a challenging task due to the noisy nature of scDNA-seq data. We introduce a new bioinformatics method called CoT for deciphering clonal copy number substructure. The backbone of CoT is a Copy number Transformer autoencoder that leverages multi-head attention mechanism to explore correlations between different genomic regions, and thus capture global features to create latent embeddings for the cells. CoT makes it convenient to first infer cell subpopulations based on the learned embeddings, and then estimate single-cell copy numbers through joint analysis of read counts data for the cells belonging to the same cluster. This exploitation of clonal substructure information in copy number analysis helps to alleviate the effect of read counts non-uniformity, and yield robust estimations of the tumor copy numbers. Performance evaluation on synthetic and real datasets showcases that CoT outperforms the state of the arts, and is highly useful for deciphering clonal copy number substructure.

POFUT1 and PLAGL2 are characteristic markers of mucinous colorectal cancer associated with MUC2 expression.

Colorectal mucinous adenocarcinoma (MAC) is one of the most lethal histological types of colorectal cancer, and its mechanism of development is not well understood. In this study, we aimed to clarify the molecular characteristics of MAC via in silico analysis using The Cancer Genome Atlas database. The expression of genes on chromosome 20q (Chr20q) was negatively associated with the expression of MUC2, which is a key molecule that can be used to distinguish between MAC and nonmucinous adenocarcinoma (NMAC). This was consistent with a significant difference in copy number alteration of Chr20q between the two histological types. We further identified 475 differentially expressed genes (DEGs) between MAC and NMAC, and some of the Chr20q genes among the DEGs are considered to be pivotal genes used to define MAC. Both in vitro and in vivo analysis showed that simultaneous knockdown of POFUT1 and PLAGL2, both of which are located on Chr20q, promoted MUC2 expression. Moreover, these genes were highly expressed in NMAC but not in MAC according to the results of immunohistological studies using human samples. In conclusion, POFUT1 and PLAGL2 are considered to be important for defining MAC, and these genes are associated with MUC2 expression.

A genome-wide study of gastric intramucosal neoplasia based on somatic copy number alterations, gene mutations, and mRNA expression patterns.

We performed comprehensive analyses of somatic copy number alterations (SCNAs) and gene expression profiles of gastric intramucosal neoplasia (IMN) using array-based methods in 97 intestinal-type IMNs, including 39 low-grade dysplasias (LGDs), 37 high-grade dysplasias (HGDs), and 26 intramucosal carcinomas (IMCs) with stromal invasion of the lamina propria to identify the molecular mechanism of IMN. In addition, we examined gene mutations using gene panel analyses. We used cluster analyses for exclusion of arbitrariness to identify SCNA patterns and expression profiles. IMNs were classified into two distinct subgroups (subgroups 1 and 2) based on SCNA patterns. Subgroup 1 showed a genomic stable pattern due to the low frequency of SCNAs, whereas subgroup 2 exhibited a chromosomal instability pattern due to the high frequencies of SCNAs and TP53 mutations. Interestingly, although the frequencies of LGD and HGD were significantly higher in subgroup 1 than in subgroup 2, IMC was commonly found in both types. Although the expression profiles of specific mRNAs could be used to categorise subgroups 1 and 2, no clinicopathological findings correlated with either subgroup. We examined signalling pathways specific to subgroups 1 and 2 to identify the association of each subgroup with signalling pathways based on gene ontology tree visualisation: subgroups 1 and 2 were associated with haem metabolism and chromosomal instability, respectively. These findings reveal a comprehensive genomic landscape that highlights the molecular complexity of IMNs and provide a road map to facilitate our understanding of gastric IMNs.

Genomic alterations in ovarian endometriosis and subsequently diagnosed ovarian carcinoma.

STUDY QUESTION: Can the alleged association between ovarian endometriosis and ovarian carcinoma be substantiated by genetic analysis of endometriosis diagnosed prior to the onset of the carcinoma? SUMMARY ANSWER: The data suggest that ovarian carcinoma does not originate from ovarian endometriosis with a cancer-like genetic profile; however, a common precursor is probable. WHAT IS KNOWN ALREADY: Endometriosis has been implicated as a precursor of ovarian carcinoma based on epidemiologic studies and the discovery of common driver mutations in synchronous disease at the time of surgery. Endometrioid ovarian carcinoma and clear cell ovarian carcinoma are the most common endometriosis-associated ovarian carcinomas (EAOCs). STUDY DESIGN, SIZE, DURATION: The pathology biobanks of two university hospitals in Sweden were scrutinized to identify women with surgically removed endometrioma who subsequently developed ovarian carcinoma (1998-2016). Only 45 archival cases with EAOC and previous endometriosis were identified and after a careful pathology review, 25 cases were excluded due to reclassification into non-EAOC (n = 9) or because ovarian endometriosis could not be confirmed (n = 16). Further cases were excluded due to insufficient endometriosis tissue or poor DNA quality in either the endometriosis, carcinoma, or normal tissue (n = 9). Finally 11 cases had satisfactory DNA from all three locations and were eligible for further analysis. PARTICIPANTS/MATERIALS, SETTING, METHODS: Epithelial cells were collected from formalin-fixed and paraffin-embedded (FFPE) sections by laser capture microdissection (endometrioma n = 11) or macrodissection (carcinoma n = 11) and DNA was extracted. Normal tissue from FFPE sections (n = 5) or blood samples collected at cancer diagnosis (n = 6) were used as the germline controls for each included patient. Whole-exome sequencing was performed (n = 33 samples). Somatic variants (single-nucleotide variants, indels, and copy number alterations) were characterized, and mutational signatures and kataegis were assessed. Microsatellite instability and mismatch repair status were confirmed with PCR and immunohistochemistry, respectively. MAIN RESULTS AND THE ROLE OF CHANCE: The median age for endometriosis surgery was 42 years, and 54 years for the subsequent ovarian carcinoma diagnosis. The median time between the endometriosis and ovarian carcinoma was 10 (7-30) years. The data showed that all paired samples harbored one or more shared somatic mutations. Non-silent mutations in cancer-associated genes were frequent in endometriosis; however, the same mutations were never observed in subsequent carcinomas. The degree of clonal dominance, demonstrated by variant allele frequency, showed a positive correlation with the time to cancer diagnosis (Spearman's rho 0.853, P < 0.001). Mutations in genes associated with immune escape were the most conserved between paired samples, and regions harboring these genes were frequently affected by copy number alterations in both sample types. Mutational burdens and mutation signatures suggested faulty DNA repair mechanisms in all cases. LARGE SCALE DATA: Datasets are available in the supplementary tables. LIMITATIONS, REASONS FOR CAUTION: Even though we located several thousands of surgically removed endometriomas between 1998 and 2016, only 45 paired samples were identified and even fewer, 11 cases, were eligible for sequencing. The observed high level of intra- and inter-heterogeneity in both groups (endometrioma and carcinoma) argues for further studies of the alleged genetic association. WIDER IMPLICATIONS OF THE FINDINGS: The observation of shared somatic mutations in all paired samples supports a common cellular origin for ovarian endometriosis and ovarian carcinoma. However, contradicting previous conclusions, our data suggest that cancer-associated mutations in endometriosis years prior to the carcinoma were not directly associated with the malignant transformation. Rather, a resilient ovarian endometriosis may delay tumorigenesis. Furthermore, the data indicate that genetic alterations affecting the immune response are early and significant events. STUDY FUNDING/COMPETING INTEREST(S): The present work has been funded by the Sjöberg Foundation (2021-01145 to K.S.; 2022-01-11:4 to A.S.), Swedish state under the agreement between the Swedish government and the county councils, the ALF-agreement (965552 to K.S.; 40615 to I.H.; 965065 to A.S.), Swedish Cancer Society (21-1848 to K.S.; 21-1684 to I.H.; 22-2080 to A.S.), BioCARE-A Strategic Research Area at Lund University (I.H. and S.W.-F.), Mrs Berta Kamprad's Cancer Foundation (FBKS-2019-28, I.H.), Cancer and Allergy Foundation (10381, I.H.), Region Västra Götaland (A.S.), Sweden's Innovation Agency (2020-04141, A.S.), Swedish Research Council (2021-01008, A.S.), Roche in collaboration with the Swedish Society of Gynecological Oncology (S.W.-F.), Assar Gabrielsson Foundation (FB19-86, C.M.), and the Lena Wäpplings Foundation (C.M.). A.S. declares stock ownership and is also a board member in Tulebovaasta, SiMSen Diagnostics, and Iscaff Pharma. A.S. has also received travel support from EMBL, Precision Medicine Forum, SLAS, and bioMCC. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

MicroRNA copy number alterations in the malignant transformation of pleomorphic adenoma to carcinoma ex pleomorphic adenoma.

OBJECTIVE: This study used array comparative genomic hybridization to assess copy number alterations (CNAs) involving miRNA genes in pleomorphic adenoma (PA), recurrent pleomorphic adenoma (RPA), residual PA, and carcinoma ex pleomorphic adenoma (CXPA). MATERIALS AND METHODS: We analyzed 13 PA, 4 RPA, 29 CXPA, and 14 residual PA using Nexus Copy Number Discovery software. The miRNAs genes affected by CNAs were evaluated based on their expression patterns and subjected to pathway enrichment analysis. RESULTS: Across the groups, we found 216 CNAs affecting 2261 miRNA genes, with 117 in PA, 59 in RPA, 846 in residual PA, and 2555 in CXPA. The chromosome 8 showed higher involvement in altered miRNAs in PAs and CXPA patients. Six miRNA genes were shared among all groups. Additionally, miR-21, miR-455-3p, miR-140, miR-320a, miR-383, miR-598, and miR-486 were prominent CNAs found and is implicated in carcinogenesis of several malignant tumors. These miRNAs regulate critical signaling pathways such as aerobic glycolysis, fatty acid biosynthesis, and cancer-related pathways. CONCLUSION: This study was the first to explore CNAs in miRNA-encoding genes in the PA-CXPA sequence. The findings suggest the involvement of numerous miRNA genes in CXPA development and progression by regulating oncogenic signaling pathways.

Tumor Copy Number Alteration Burden as a Predictor for Resistance to Immune Checkpoint Blockade across Different Cancer Types.

Immune checkpoint blockade (ICB) benefits only a subset of advanced cancer patients, and predictive biomarkers for immunotherapy response are needed. Recently, copy number alteration (CNA) burden has been proposed to predict ICB resistance. We assessed this finding using the publicly accessible data for 1661 ICB-treated patients whose tumors were profiled by MSK-IMPACT, an approved targeted assay in clinical care. We tested the hypothesis that the continuous increase in CNA burden is associated with poor overall survival following ICB. In addition, we hypothesized that the combinatorial biomarkers of tumor mutational burden (TMB) and CNA burden would better stratify patients for immune status and ICB response. Of the 1661 cases, 79% (n = 1307) were treated with anti PD-1/PD-L1 and the remaining 21% (n = 354) with anti CTLA-4 or the combination of both. In a multivariate analysis, increase in CNA burden was associated with poor overall survival [HR = 1.52, 95% CI (1.01-2.30), p = 0.04]. The combination of biomarkers TMB and CNA burden stratified patients into four clinically distinct subsets among which "LowTMB/HighCNA" showed the worst survival (p < 0.0001). The four patient subsets had unique CNA profiles and enriched pathways, which could predict transcriptional and phenotypic effects related to immune signaling and CD8+ T-cell abundance in the tumor microenvironment. CNA burden was associated with poor overall survival in patients receiving ICB and could improve patient stratification when incorporated with TMB. These findings may guide patient selection for immunotherapy or alternative strategies.

Combination of minimal residual disease on day 15 and copy number alterations results in BCR-ABL1-negative pediatric B-ALL: A powerful tool for prediction of induction failure.

The current genomic abnormalities provide prognostic value in pediatric Acute Lymphoblastic Leukemia (ALL). Furthermore, Copy Number Alteration (CNA) has recently been used to improve the genetic risk stratification of patients. This study aimed to evaluate CNA profiles in BCR-ABL1-negative pediatric B-ALL patients and correlate the data with Minimal Residual Disease (MRD) results after induction therapy. We examined 82 bone marrow samples from pediatric BCR-ABL1-negative B-ALL using the MLPA method for the most common CNAs, including IKZF1, CDKN2A/B, PAX5, RB1, BTG1, ETV6, EBF1, JAK2, and PAR1 region. Subsequently, patients were followed-up by multiparameter Flow Cytometry for MRD (MFC-MRD) assessment on days 15 and 33 after induction. Data showed that 58.5 % of patients carried at least one gene deletion, whereas 41.7 % of them carried more than one gene deletion simultaneously. The most frequent gene deletions were CDKN2A/B, ETV6, and IKZF1 (30.5 %, 14.6 %, and 14.6 %, respectively), while the PAR1 region showed predominantly duplication (30.5 %). CDKN2A/B and IKZF1 were related to positive MRD results on day 15 (p = 0.003 and p = 0.007, respectively). The simultaneous presence of more than one deletion was significantly associated with high induction failure (p = 0.001). Also, according to the CNA profile criteria, the CNA with poor risk (CNA-PR) profile was statistically associated with older age and positive MRD results on day 15 (p = 0.014 and p = 0.013, respectively). According to our results, the combined use of CNAs with MRD results on day 15 can predict induction failure and be helpful in ameliorating B-ALL risk stratification and treatment approaches.

Inferring single-cell copy number profiles through cross-cell segmentation of read counts.

BACKGROUND: Copy number alteration (CNA) is one of the major genomic variations that frequently occur in cancers, and accurate inference of CNAs is essential for unmasking intra-tumor heterogeneity (ITH) and tumor evolutionary history. Single-cell DNA sequencing (scDNA-seq) makes it convenient to profile CNAs at single-cell resolution, and thus aids in better characterization of ITH. Despite that several computational methods have been proposed to decipher single-cell CNAs, their performance is limited in either breakpoint detection or copy number estimation due to the high dimensionality and noisy nature of read counts data. RESULTS: By treating breakpoint detection as a process to segment high dimensional read count sequence, we develop a novel method called DeepCNA for cross-cell segmentation of read count sequence and per-cell inference of CNAs. To cope with the difficulty of segmentation, an autoencoder (AE) network is employed in DeepCNA to project the original data into a low-dimensional space, where the breakpoints can be efficiently detected along each latent dimension and further merged to obtain the final breakpoints. Unlike the existing methods that manually calculate certain statistics of read counts to find breakpoints, the AE model makes it convenient to automatically learn the representations. Based on the inferred breakpoints, we employ a mixture model to predict copy numbers of segments for each cell, and leverage expectation-maximization algorithm to efficiently estimate cell ploidy by exploring the most abundant copy number state. Benchmarking results on simulated and real data demonstrate our method is able to accurately infer breakpoints as well as absolute copy numbers and surpasses the existing methods under different test conditions. DeepCNA can be accessed at: https://github.com/zhyu-lab/deepcna . CONCLUSIONS: Profiling single-cell CNAs based on deep learning is becoming a new paradigm of scDNA-seq data analysis, and DeepCNA is an enhancement to the current arsenal of computational methods for investigating cancer genomics.

Isolation and Genomic Analysis of Circulating Tumor Cell Clusters in Cancer Patients.

The role of circulating tumor cell (CTC) clusters in the metastatic dissemination process is gaining increased attention. Besides homotypic clusters, heterotypic clusters that contain tumor cells admixed with normal cells are frequently observed in patients with solid tumors. Current methods used for cluster detection and enumeration do not allow an accurate estimation of the relative fractions of tumor cells. Here we describe a method for estimating tumor fraction of clusters including isolation and collection of single clusters, assessment of copy number alterations of single clusters by low-pass whole genome sequencing, and bioinformatic analysis of sequencing data.

Comprehensive Genomic Analysis of Cemento-Ossifying Fibroma.

Cemento-ossifying fibroma (COF) of the jaws is currently classified as a benign mesenchymal odontogenic tumor, and only targeted approaches have been used to assess its genetic alterations. A minimal proportion of COFs harbor CDC73 somatic mutations, and copy number alterations (CNAs) involving chromosomes 7 and 12 have recently been reported in a small proportion of cases. However, the genetic background of COFs remains obscure. We used a combination of whole-exome sequencing and RNA sequencing to assess somatic mutations, fusion transcripts, and CNAs in a cohort of 12 freshly collected COFs. No recurrent fusions have been identified among the 5 cases successfully analyzed by RNA sequencing, with in-frame fusions being detected in 2 cases (MARS1::GOLT1B and PARG::BMS1 in one case and NCLN::FZR1 and NFIC::SAMD1 in the other case) and no candidate fusions identified for the remaining 3 cases. No recurrent pathogenic mutations were detected in the 11 cases that had undergone whole-exome sequencing. A KRAS p.L19F missense variant was detected in one case, and 2 CDC73 deletions were detected in another case. The other variants were of uncertain significance and included variants in PC, ACTB, DOK6, HACE1, and COL1A2 and previously unreported variants in PTPN14, ATP5F1C, APOBEC1, HDAC5, ATF7IP, PARP2, and ACTR3B. The affected genes do not clearly converge on any signaling pathway. CNAs were detected in 5/11 cases (45%), with copy gains involving chromosome 12 occurring in 3/11 cases (27%). In conclusion, no recurrent fusions or pathogenic variants have been detected in the present COF cohort, with copy gains involving chromosome 12 occurring in 27% of cases.

Genomic alterations driving precancerous to cancerous lesions in esophageal cancer development.

Esophageal squamous cell carcinoma (ESCC) develops through a series of increasingly abnormal precancerous lesions. Previous studies have revealed the striking differences between normal esophageal epithelium and ESCC in copy number alterations (CNAs) and mutations in genes driving clonal expansion. However, due to limited data on early precancerous lesions, the timing of these transitions and which among them are prerequisites for malignant transformation remained unclear. Here, we analyze 1,275 micro-biopsies from normal esophagus, early and late precancerous lesions, and esophageal cancers to decipher the genomic alterations at each stage. We show that the frequency of TP53 biallelic inactivation increases dramatically in early precancerous lesion stage while CNAs and APOBEC mutagenesis substantially increase at late stages. TP53 biallelic loss is the prerequisite for the development of CNAs of genes in cell cycle, DNA repair, and apoptosis pathways, suggesting it might be one of the earliest steps initiating malignant transformation.

Association of CDKN2A/2B deletion with relapse after hematopoietic stem cell transplantation for acute lymphoblastic leukemia.

The most important prognostic factor for Philadelphia chromosome-positive acute lymphoblastic leukemia (Ph+ALL) is minimal residual disease (MRD). Previous studies have reported copy number variants of genes such as IKZF1, CDKN2A/2B, and PAX5. These gene mutations can be analyzed using multiplex ligation-dependent probe amplification (MLPA), which is less costly and easier to perform than large-scale gene mutation analyses. In this study, we performed copy number variant analysis of leukemia cells at the first onset of Ph+ALL in a case series of allogeneic hematopoietic stem cell transplantation (allo-HSCT) using the MLPA method. We analyzed how it influenced allo-HSCT prognosis together with MRD information. CDKN2A/2B copy number variations significantly increased the rate of post-transplant recurrence (P=0.025) and significantly reduced disease-free survival (P=0.015). Additionally, patients with IKZF1 deletions had a significantly higher post-transplant recurrence rate (P=0.042). Although they were positive for pre-transplant MRD, no relapse was observed in patients with wild-type copy number variations in IKZF1 or CDKN2A/2B. CDKN2A/2B copy number variation is a crucial factor that can be confirmed at initial onset as a post-transplant prognostic factor of Ph+ALL.

Classification of HER2-negative breast cancers by ERBB2 copy number alteration status reveals molecular differences associated with chromosome 17 gene aberrations.

Background: Recently, HER2-negative breast cancers have been reclassified by protein expression into 'HER2-low' and 'HER2-zero' subgroups, but the consideration of HER2-low breast cancer as a distinct biological subtype with differing prognoses remains controversial. By contrast, non-neutral ERBB2 copy number alteration (CNA) status is associated with inferior survival outcomes compared to ERBB2 CNA-neutral breast cancer, providing an alternative approach to classification. Methods: Here, we investigated the molecular landscape of non-metastatic HER2-negative BCs in relation to ERBB2 CNA status to elucidate biological differences. Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and The Cancer Genome Atlas (TCGA) TCGA-BRCA datasets (n = 1875) were analyzed. Results: Nearly two-fifths of the cohort harbored ERBB2 CNAs (39.4%), which were significantly enriched within hormone receptor-negative (56.1%) than within hormone receptor-positive BCs (35.5%; p < 0.0001). Globally, CNAs across the genome were significantly higher in ERBB2 non-neutral compared to neutral cohorts (p < 0.0001). Notably, genetic aberrations on chromosome 17 - BRCA1, NF1, TP53, MAP2K4, and NCOR1 - were widespread in the ERBB2 non-neutral cases. While chromosome 17q arm-level alterations were largely in tandem with ERBB2 CNA status, arm-level loss in chromosome 17p was prevalent regardless of ERBB2 gain, amplification, or loss. Differential gene expression analysis demonstrated that pathways involved in the cell cycle, proteasome, and DNA replication were upregulated in ERBB2 non-neutral cases. Conclusion: Classification of HER2-negative BCs according to ERBB2 CNA status reveals differences in the genomic landscape. The implications of concurrent aberrations in other genes on chromosome 17 merit further research in ERBB2 non-neutral BCs.

Functional screening of amplification outlier oncogenes in organoid models of early tumorigenesis.

Somatic copy number gains are pervasive across cancer types, yet their roles in oncogenesis are insufficiently evaluated. This inadequacy is partly due to copy gains spanning large chromosomal regions, obscuring causal loci. Here, we employed organoid modeling to evaluate candidate oncogenic loci identified via integrative computational analysis of extreme copy gains overlapping with extreme expression dysregulation in The Cancer Genome Atlas. Subsets of "outlier" candidates were contextually screened as tissue-specific cDNA lentiviral libraries within cognate esophagus, oral cavity, colon, stomach, pancreas, and lung organoids bearing initial oncogenic mutations. Iterative analysis nominated the kinase DYRK2 at 12q15 as an amplified head and neck squamous carcinoma oncogene in p53-/- oral mucosal organoids. Similarly, FGF3, amplified at 11q13 in 41% of esophageal squamous carcinomas, promoted p53-/- esophageal organoid growth reversible by small molecule and soluble receptor antagonism of FGFRs. Our studies establish organoid-based contextual screening of candidate genomic drivers, enabling functional evaluation during early tumorigenesis.

Identifying the driver miRNAs with somatic copy number alterations driving dysregulated ceRNA networks in cancers.

BACKGROUND: MicroRNAs (miRNAs) play critical roles in cancer initiation and progression, which were critical components to maintain the dynamic balance of competing endogenous RNA (ceRNA) networks. Somatic copy number alterations (SCNAs) in the cancer genome could disturb the transcriptome level of miRNA to deregulate this balance. However, the driving effects of SCNAs of miRNAs were insufficiently understood. METHODS: In this study, we proposed a method to dissect the functional roles of miRNAs under different copy number states and identify driver miRNAs by integrating miRNA SCNAs profile, miRNA-target relationships and expression profiles of miRNA, mRNA and lncRNA. RESULTS: Applying our method to 813 TCGA breast cancer (BRCA) samples, we identified 29 driver miRNAs whose SCNAs significantly and concordantly regulated their own expression levels and further inversely dysregulated expression levels of their targets or disturbed the miRNA-target networks they directly involved. Based on miRNA-target networks, we further constructed dynamic ceRNA networks driven by driver SCNAs of miRNAs and identified three different patterns of SCNA interference in the miRNA-mediated dynamic ceRNA networks. Survival analysis of driver miRNAs showed that high-level amplifications of four driver miRNAs (including has-miR-30d-3p, has-mir-30b-5p, has-miR-30d-5p and has-miR-151a-3p) in 8q24 characterized a new BRCA subtype with poor prognosis and contributed to the dysfunction of cancer-associated hallmarks in a complementary way. The SCNAs of driver miRNAs across different cancer types contributed to the cancer development by dysregulating different components of the same cancer hallmarks, suggesting the cancer specificity of driver miRNA. CONCLUSIONS: These results demonstrate the efficacy of our method in identifying driver miRNAs and elucidating their functional roles driven by endogenous SCNAs, which is useful for interpreting cancer genomes and pathogenic mechanisms.

Cracking the pattern of tumor evolution based on single-cell copy number alterations.

Copy number alterations (CNAs) are a key characteristic of tumor development and progression. The accumulation of various CNAs during tumor development plays a critical role in driving tumor evolution. Heterogeneous clones driven by distinct CNAs have different selective advantages, leading to differential patterns of tumor evolution that are essential for developing effective cancer therapies. Recent advances in single-cell sequencing technology have enabled genome-wide copy number profiling of tumor cell populations at single-cell resolution. This has made it possible to explore the evolutionary patterns of CNAs and accurately discover the mechanisms of intra-tumor heterogeneity. Here, we propose a two-step statistical approach that distinguishes neutral, linear, branching and punctuated evolutionary patterns for a tumor cell population based on single-cell copy number profiles. We assessed our approach using a variety of simulated and real single-cell genomic and transcriptomic datasets, demonstrating its high accuracy and robustness in predicting tumor evolutionary patterns. We applied our approach to single-cell DNA sequencing data from 20 breast cancer patients and observed that punctuated evolution is the dominant evolutionary pattern in breast cancer. Similar conclusions were drawn when applying the approach to single-cell RNA sequencing data obtained from 132 various cancer patients. Moreover, we found that differential immune cell infiltration is associated with specific evolutionary patterns. The source code of our study is available at https://github.com/FangWang-SYSU/PTEM.