HATCHet2: clone- and haplotype-specific copy number inference from bulk tumor sequencing data.

Bulk DNA sequencing of multiple samples from the same tumor is becoming common, yet most methods to infer copy-number aberrations (CNAs) from this data analyze individual samples independently. We introduce HATCHet2, an algorithm to identify haplotype- and clone-specific CNAs simultaneously from multiple bulk samples. HATCHet2 extends the earlier HATCHet method by improving identification of focal CNAs and introducing a novel statistic, the minor haplotype B-allele frequency (mhBAF), that enables identification of mirrored-subclonal CNAs. We demonstrate HATCHet2's improved accuracy using simulations and a single-cell sequencing dataset. HATCHet2 analysis of 10 prostate cancer patients reveals previously unreported mirrored-subclonal CNAs affecting cancer genes.

Accurate detection of copy number aberrations in FFPE samples using the mFAST-SeqS approach.

BACKGROUND: Shallow whole genome sequencing (Shallow-seq) is used to determine the copy number aberrations (CNA) in tissue samples and circulating tumor DNA. However, costs of NGS and challenges of small biopsies ask for an alternative to the untargeted NGS approaches. The mFAST-SeqS approach, relying on LINE-1 repeat amplification, showed a good correlation with Shallow-seq to detect CNA in blood samples. In the present study, we evaluated whether mFAST-SeqS is suitable to assess CNA in small formalin-fixed paraffin-embedded (FFPE) tissue specimens, using vulva and anal HPV-related lesions. METHODS: Seventy-two FFPE samples, including 36 control samples (19 vulva;17 anal) for threshold setting and 36 samples (24 vulva; 12 anal) for clinical evaluation, were analyzed by mFAST-SeqS. CNA in vulva and anal lesions were determined by calculating genome-wide and chromosome arm-specific z-scores in comparison with the respective control samples. Sixteen samples were also analyzed with the conventional Shallow-seq approach. RESULTS: Genome-wide z-scores increased with the severity of disease, with highest values being found in cancers. In vulva samples median and inter quartile ranges [IQR] were 1[0-2] in normal tissues (n = 4), 3[1-7] in premalignant lesions (n = 9) and 21[13-48] in cancers (n = 10). In anal samples, median [IQR] were 0[0-1] in normal tissues (n = 4), 14[6-38] in premalignant lesions (n = 4) and 18[9-31] in cancers (n = 4). At threshold 4, all controls were CNA negative, while 8/13 premalignant lesions and 12/14 cancers were CNA positive. CNA captured by mFAST-SeqS were mostly also found by Shallow-seq. CONCLUSION: mFAST-SeqS is easy to perform, requires less DNA and less sequencing reads reducing costs, thereby providing a good alternative for Shallow-seq to determine CNA in small FFPE samples.

NestedBD: Bayesian inference of phylogenetic trees from single-cell copy number profiles under a birth-death model.

Copy number aberrations (CNAs) are ubiquitous in many types of cancer. Inferring CNAs from cancer genomic data could help shed light on the initiation, progression, and potential treatment of cancer. While such data have traditionally been available via "bulk sequencing," the more recently introduced techniques for single-cell DNA sequencing (scDNAseq) provide the type of data that makes CNA inference possible at the single-cell resolution. We introduce a new birth-death evolutionary model of CNAs and a Bayesian method, NestedBD, for the inference of evolutionary trees (topologies and branch lengths with relative mutation rates) from single-cell data. We evaluated NestedBD's performance using simulated data sets, benchmarking its accuracy against traditional phylogenetic tools as well as state-of-the-art methods. The results show that NestedBD infers more accurate topologies and branch lengths, and that the birth-death model can improve the accuracy of copy number estimation. And when applied to biological data sets, NestedBD infers plausible evolutionary histories of two colorectal cancer samples. NestedBD is available at https://github.com/Androstane/NestedBD .

Evolution of chromosome-arm aberrations in breast cancer through genetic network rewiring.

The basal breast cancer subtype is enriched for triple-negative breast cancer (TNBC) and displays consistent large chromosomal deletions. Here, we characterize evolution and maintenance of chromosome 4p (chr4p) loss in basal breast cancer. Analysis of The Cancer Genome Atlas data shows recurrent deletion of chr4p in basal breast cancer. Phylogenetic analysis of a panel of 23 primary tumor/patient-derived xenograft basal breast cancers reveals early evolution of chr4p deletion. Mechanistically we show that chr4p loss is associated with enhanced proliferation. Gene function studies identify an unknown gene, C4orf19, within chr4p, which suppresses proliferation when overexpressed-a member of the PDCD10-GCKIII kinase module we name PGCKA1. Genome-wide pooled overexpression screens using a barcoded library of human open reading frames identify chromosomal regions, including chr4p, that suppress proliferation when overexpressed in a context-dependent manner, implicating network interactions. Together, these results shed light on the early emergence of complex aneuploid karyotypes involving chr4p and adaptive landscapes shaping breast cancer genomes.

Application of a Multiplex Ligation-Dependent Probe Amplification-Based Next-Generation Sequencing Approach for the Detection of Pathogenesis of Duchenne Muscular Dystrophy and Spinal Muscular Atrophy Caused by Copy Number Aberrations.

Both Duchenne muscular dystrophy (DMD; OMIM no. 310200) and spinal muscular atrophy (SMA; OMIM no. 253300/253550/253400/271150) are genetic disorders characterized by progressive muscle degeneration and weakness. Genetic copy number aberrations in the pathogenetic genes DMD and SMN1 lead to alterations in functional proteins, resulting in DMD and SMA, respectively. Multiplex ligation-dependent probe amplification (MLPA) has become a standard method for the detection of common copy number aberrations (CNAs), including DMD and SMN1 deletions, both of which are associated with poor clinical outcomes. However, traditional MLPA assays only accommodate a maximum of 60 MLPA probes per test. To increase the number of targeted sequences in one assay, an MLPA-based next-generation sequencing (NGS) assay has been developed that is based on the standard MLPA procedure, allows high-throughput screening for a large number of fragments and samples by integrating additional indices for detection, and can be analyzed on all Illumina NGS platforms.

A Gradual Transition Toward Anaplasia in Wilms Tumor Through Tolerance to Genetic Damage.

Patients with Wilms tumor (WT) in general have excellent survival, but the prognosis of patients belonging to the subgroup of WT with diffuse anaplasia (DA) is poor due to frequent resistance to chemotherapy. We hypothesized that DA WT cells might undergo changes, such as acquiring a persistent tolerance to DNA damage and copy number aberrations (CNAs), which could eventually lead to their resistance to chemotherapy treatment. Tissue sections from chemotherapy-treated DA WTs (n = 12) were compared with chemotherapy-treated nonanaplastic WTs (n = 15) in a tissue microarray system, enabling analysis of 769 tumor regions. All regions were scored for anaplastic features and immunohistochemistry was used to quantify p53 expression, proliferation index (Ki67), and DNA double-strand breaks (γH2AX). CNAs were assessed by array-based genotyping and TP53 mutations using targeted sequencing. Proliferation index and the frequency of DNA double-strand breaks (γH2AX dot expression) increased with higher anaplasia scores. Almost all (95.6%) areas with full-scale anaplasia had TP53 mutations or loss of heterozygosity, along with an increased amount of CNAs. Interestingly, areas with wild-type TP53 with loss of heterozygosity and only one feature of anaplasia (anaplasia score 1) also had significantly higher proliferation indices, more DNA double-strand breaks, and more CNAs than regions without any anaplastic features (score 0); such areas may be preanaplastic cell populations under selective pressure for TP53 mutations. In conclusion, we suggest that chemoresistance of DA WTs may be partly explained by a high proliferative capability of anaplastic cells, which also have a high burden of double-stranded DNA breaks and CNAs, and that there is a gradual emergence of anaplasia in WT.

Transcriptomic subtyping of malignant peripheral nerve sheath tumours highlights immune signatures, genomic profiles, patient survival and therapeutic targets.

BACKGROUND: Malignant peripheral nerve sheath tumour (MPNST) is an aggressive orphan disease commonly affecting adolescents or young adults. Current knowledge of molecular tumour biology has been insufficient for development of rational treatment strategies. We aimed to discover molecular subtypes of potential clinical relevance. METHODS: Fresh frozen samples of MPNSTs (n = 94) and benign neurofibromas (n = 28) from 115 patients in a European multicentre study were analysed by DNA copy number and/or transcriptomic profiling. Unsupervised transcriptomic subtyping was performed and the subtypes characterized for genomic aberrations, clinicopathological associations and patient survival. FINDINGS: MPNSTs were classified into two transcriptomic subtypes defined primarily by immune signatures and proliferative processes. "Immune active" MPNSTs (44%) had sustained immune signals relative to neurofibromas, were more frequently low-grade (P = 0.01) and had favourable prognostic associations in a multivariable model of disease-specific survival with clinicopathological factors (hazard ratio 0.25, P = 0.003). "Immune deficient" MPNSTs were more aggressive and characterized by proliferative signatures, high genomic complexity, aberrant TP53 and PRC2 loss, as well as high relative expression of several potential actionable targets (EGFR, ERBB2, EZH2, KIF11, PLK1, RRM2). Integrated gene-wise analyses suggested a DNA copy number-basis for proliferative transcriptomic signatures in particular, and the tumour copy number burden further stratified the transcriptomic subtypes according to patient prognosis (P < 0.01). INTERPRETATION: Approximately half of MPNSTs belong to an "immune deficient" transcriptomic subtype associated with an aggressive disease course, PRC2 loss and expression of several potential therapeutic targets, providing a rationale for molecularly-guided intervention trials. FUNDING: Research grants from non-profit organizations, as stated in the Acknowledgements.

Mutations of BRCA1, BRCA2, and PALB2 Genes in Breast Tumor Tissue: Relationship with the Effectiveness of Neoadjuvant Chemotherapy and Disease Prognosis.

It has been shown that the loss of function of the BRCA1, BRCA2, and PALB2 genes due to a number of hereditary mutations or chromosomal aberrations can affect the effectiveness of chemotherapy treatment and disease prognosis in patients with various types of cancer, and in particular in breast cancer. Thus, the aim of the work was to evaluate the predictive and prognostic potential of DNA copy number aberrations and mutations in the BRCA1, BRCA2, and PALB2 genes in breast tumors. MATERIALS AND METHODS: The study included 66 patients with breast cancer. DNA copy number aberrations (CNA) were assessed by high-density CytoScanHD™ Array micro matrix analysis. Gene mutations were assessed by sequencing on the MiSeq™ Sequencing System using the Accel-Amplicon BRCA1, BRCA2, and PALB2 Panel. RESULTS: It has been established that the presence of a normal copy number of PALB2 is associated with a lack of response to chemotherapy in Taxotere-containing treatment regimens (p = 0.05). In addition, the presence of a PALB2 deletion is associated with 100% metastatic survival rates (log-rank test p = 0.04). As a result of sequencing, 25 mutations were found in the BRCA1 gene, 42 mutations in BRCA2, and 27 mutations in the PALB2 gene. The effect of mutations on the effectiveness of treatment is controversial, but an effect on the survival of patients with breast cancer has been shown. So, in the presence of pathogenic mutations in the BRCA2 gene, 100% metastatic survival is observed (log-rank test p = 0.05), as well as in the elimination of PALB2 mutations during treatment (log-rank test p = 0.07). CONCLUSION: Currently, there is little data on the effect of chromosomal aberrations and mutations in the BRCA1/2 and PALB2 genes on the effectiveness of treatment and prognosis of the disease. At the same time, the study of these genes has great potential for testing focused on a personalized approach to the treatment of patients with breast cancer.

HATCHet2: clone- and haplotype-specific copy number inference from bulk tumor sequencing data.

Multi-region DNA sequencing of primary tumors and metastases from individual patients helps identify somatic aberrations driving cancer development. However, most methods to infer copy-number aberrations (CNAs) analyze individual samples. We introduce HATCHet2 to identify haplotype- and clone-specific CNAs simultaneously from multiple bulk samples. HATCHet2 introduces a novel statistic, the mirrored haplotype B-allele frequency (mhBAF), to identify mirrored-subclonal CNAs having different numbers of copies of parental haplotypes in different tumor clones. HATCHet2 also has high accuracy in identifying focal CNAs and extends the earlier HATCHet method in several directions. We demonstrate HATCHet2's improved accuracy using simulations and a single-cell sequencing dataset. HATCHet2 analysis of 50 prostate cancer samples from 10 patients reveals previously-unreported mirrored-subclonal CNAs affecting cancer genes.

Dissecting Molecular Heterogeneity of Circulating Tumor Cells (CTCs) from Metastatic Breast Cancer Patients through Copy Number Aberration (CNA) and Single Nucleotide Variant (SNV) Single Cell Analysis.

Circulating tumor cells' (CTCs) heterogeneity contributes to counteract their introduction in clinical practice. Through single-cell sequencing we aim at exploring CTC heterogeneity in metastatic breast cancer (MBC) patients. Single CTCs were isolated using DEPArray NxT. After whole genome amplification, libraries were prepared for copy number aberration (CNA) and single nucleotide variant (SNV) analysis and sequenced using Ion GeneStudio S5 and Illumina MiSeq, respectively. CTCs demonstrate distinctive mutational signatures but retain molecular traces of their common origin. CNA profiling identifies frequent aberrations involving critical genes in pathogenesis: gains of 1q (CCND1) and 11q (WNT3A), loss of 22q (CHEK2). The longitudinal single-CTC analysis allows tracking of clonal selection and the emergence of resistance-associated aberrations, such as gain of a region in 12q (CDK4). A group composed of CTCs from different patients sharing common traits emerges. Further analyses identify losses of 15q and enrichment of terms associated with pseudopodium formation as frequent and exclusive events. CTCs from MBC patients are heterogeneous, especially concerning their mutational status. The single-cell analysis allows the identification of aberrations associated with resistance, and is a candidate tool to better address treatment strategy. The translational significance of the group populated by similar CTCs should be elucidated.

TAaCGH Suite for Detecting Cancer-Specific Copy Number Changes Using Topological Signatures.

Copy number changes play an important role in the development of cancer and are commonly associated with changes in gene expression. Persistence curves, such as Betti curves, have been used to detect copy number changes; however, it is known these curves are unstable with respect to small perturbations in the data. We address the stability of lifespan and Betti curves by providing bounds on the distance between persistence curves of Vietoris-Rips filtrations built on data and slightly perturbed data in terms of the bottleneck distance. Next, we perform simulations to compare the predictive ability of Betti curves, lifespan curves (conditionally stable) and stable persistent landscapes to detect copy number aberrations. We use these methods to identify significant chromosome regions associated with the four major molecular subtypes of breast cancer: Luminal A, Luminal B, Basal and HER2 positive. Identified segments are then used as predictor variables to build machine learning models which classify patients as one of the four subtypes. We find that no single persistence curve outperforms the others and instead suggest a complementary approach using a suite of persistence curves. In this study, we identified new cytobands associated with three of the subtypes: 1q21.1-q25.2, 2p23.2-p16.3, 23q26.2-q28 with the Basal subtype, 8p22-p11.1 with Luminal B and 2q12.1-q21.1 and 5p14.3-p12 with Luminal A. These segments are validated by the TCGA BRCA cohort dataset except for those found for Luminal A.

Longitudinal monitoring of disease burden and response using ctDNA from dried blood spots in xenograft models.

Whole-genome sequencing (WGS) of circulating tumour DNA (ctDNA) is now a clinically important biomarker for predicting therapy response, disease burden and disease progression. However, the translation of ctDNA monitoring into vital preclinical PDX models has not been possible owing to low circulating blood volumes in small rodents. Here, we describe the longitudinal detection and monitoring of ctDNA from minute volumes of blood in PDX mice. We developed a xenograft Tumour Fraction (xTF) metric using shallow WGS of dried blood spots (DBS), and demonstrate its application to quantify disease burden, monitor treatment response and predict disease outcome in a preclinical study of PDX mice. Further, we show how our DBS-based ctDNA assay can be used to detect gene-specific copy number changes and examine the copy number landscape over time. Use of sequential DBS ctDNA assays could transform future trial designs in both mice and patients by enabling increased sampling and molecular monitoring.

Whole-genome sequencing as an alternative to analyze copy number abnormalities in acute myeloid leukemia and myelodysplastic syndrome.

Copy number aberrations (CNA) are the core determinants for diagnosis, risk stratification and prognosis in acute myeloid leukemia (AML) and myelodysplastic syndrome (MDS). In this study, a shallow whole-genome sequencing-based assay, LeukoPrint, was utilized to depict genomic CNA profiles from the bone marrow of 137 newly diagnosed AML/MDS patients. It demonstrated 98.1% concordance of CNA profiles with cytogenetics and/or fluorescence in situ hybridization (FISH). It is advantageous in detecting CNAs of short segments (1 Mb) and from samples with low leukemic cell content, more accurate for describing complex karyotypes and less confounded by subjective bias. LeukoPrint improved the overall diagnostic yield by redefining the risk categories for 16 patients by presenting new information. In summary, LeukoPrint provided an automated, convenient, and cost-effective approach to describe genomic CNA profiles. It brought greater diagnostic yield and risk stratification information by incorporating into the routine cytogenetics based on the CNA-related criteria of standard ELN/IPSS-R guidelines.

DNA Copy Number Aberrations and Expression of ABC Transporter Genes in Breast Tumour: Correlation with the Effect of Neoadjuvant Chemotherapy and Prognosis of the Disease.

One of the important reasons for the ineffectiveness of chemotherapy in breast cancer (BC) is considered to be the formation of a multidrug resistance phenotype in tumour cells, which is caused by the expression of energy-dependent ABC transporters. The aim of this work was to assess chromosomal aberrations and the level of transcripts of all 49 known ABC transporter genes in breast tumours. MATERIALS AND METHODS: The study included 129 patients with breast cancer. A microarray study of all tumour samples was carried out on microchips. RESULTS: This study established that the presence of a deletion in genes ABCB1, ABCB4, ABCB8, ABCC7, ABCC11, ABCC12, ABCF2, and ABCG4 is associated with an objective response to treatment (p ≤ 0.05). A decrease in the expression of genes was associated with a good response to chemotherapy, whereas an increase in expression caused the progression and stabilization of the tumour. Analysis of metastatic-free survival rates showed that the presence of ABCB1/4 and ABCC1/6 deletions was associated with 100% survival (log-rank test p = 0.01 and p = 0.03). CONCLUSIONS: The study showed that the aberrant state of ABC transporter genes, as well as a decrease in the expression of these genes, is a predictor of the effectiveness of therapeutic treatment and a potential prognostic marker of metastatic survival.

Recent Advances in Pathology: the 2022 Annual Review Issue of The Journal of Pathology.

The 2022 Annual Review Issue of The Journal of Pathology, Recent Advances in Pathology, contains 15 invited reviews on research areas of growing importance in pathology. This year, the articles include those that focus on digital pathology, employing modern imaging techniques and software to enable improved diagnostic and research applications to study human diseases. This subject area includes the ability to identify specific genetic alterations through the morphological changes they induce, as well as integrating digital and computational pathology with 'omics technologies. Other reviews in this issue include an updated evaluation of mutational patterns (mutation signatures) in cancer, the applications of lineage tracing in human tissues, and single cell sequencing technologies to uncover tumour evolution and tumour heterogeneity. The tissue microenvironment is covered in reviews specifically dealing with proteolytic control of epidermal differentiation, cancer-associated fibroblasts, field cancerisation, and host factors that determine tumour immunity. All of the reviews contained in this issue are the work of invited experts selected to discuss the considerable recent progress in their respective fields and are freely available online (https://onlinelibrary.wiley.com/journal/10969896). © 2022 The Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Single-Cell Molecular Characterization to Partition the Human Glioblastoma Tumor Microenvironment Genetic Background.

BACKGROUND: Glioblastoma (GB) is a devastating primary brain malignancy. The recurrence of GB is inevitable despite the standard treatment of surgery, chemotherapy, and radiation, and the median survival is limited to around 15 months. The barriers to treatment include the complex interactions among the different cellular components inhabiting the tumor microenvironment. The complex heterogeneous nature of GB cells is helped by the local inflammatory tumor microenvironment, which mostly induces tumor aggressiveness and drug resistance. METHODS: By using fluorescent multiple labeling and a DEPArray cell separator, we recovered several single cells or groups of single cells from populations of different origins from IDH-WT GB samples. From each GB sample, we collected astrocytes-like (GFAP+), microglia-like (IBA1+), stem-like cells (CD133+), and endothelial-like cells (CD105+) and performed Copy Number Aberration (CNA) analysis with a low sequencing depth. The same tumors were subjected to a bulk CNA analysis. RESULTS: The tumor partition in its single components allowed single-cell molecular subtyping which revealed new aspects of the GB altered genetic background. CONCLUSIONS: Nowadays, single-cell approaches are leading to a new understanding of GB physiology and disease. Moreover, single-cell CNAs resource will permit new insights into genome heterogeneity, mutational processes, and clonal evolution in malignant tissues.

Leveraging single-cell sequencing to unravel intratumour heterogeneity and tumour evolution in human cancers.

Intratumour heterogeneity (ITH) and tumour evolution are well-documented phenomena in human cancers. While the advent of next-generation sequencing technologies has facilitated the large-scale capture of genomic data, the field of single-cell genomics is nascent but rapidly advancing and generating many new insights into the complex molecular mechanisms of tumour biology. In this review, we provide an overview of current single-cell DNA sequencing technologies, exploring how recent methodological advancements have enumerated new insights into ITH and tumour evolution. Areas highlighted include the potential power of single-cell genome sequencing studies to explore evolutionary dynamics contributing to tumourigenesis through to progression, metastasis, and therapy resistance. We also explore the use of in situ sequencing technologies to study ITH in a spatial context, as well as examining the use of single-cell genomics to perform lineage tracing in both normal and malignant tissues. Finally, we consider the use of multimodal single-cell sequencing technologies. Taken together, it is hoped that these many facets of single-cell genome sequencing will improve our understanding of tumourigenesis, progression, and lethality in cancer, leading to the development of novel therapies. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Predictive and Prognostic Significance of mRNA Expression and DNA Copies Aberrations of ERCC1, RRM1, TOP1, TOP2A, TUBB3, TYMS, and GSTP1 Genes in Patients with Breast Cancer.

Increasingly, many researchers are focusing on the sensitivity in breast tumors (BC) to certain chemotherapy drugs and have personalized their research based on the assessment of this sensitivity. One such personalized approach is to assess the chemotherapy's gene expression, as well as aberrations in the number of DNA copies-deletions and amplifications with the ability to have a significant effect on the gene's activity. Thus, the aim of this work was to study the predictive and prognostic significance of the expression and chromosomal aberrations of eight chemosensitivity genes in breast cancer patients. MATERIAL AND METHODS: The study involved 97 patients with luminal B breast cancer IIB-IIIB stages. DNA and RNA were isolated from samples of tumor tissue before and after treatment. Microarray analysis was performed for all samples on high-density microarrays (DNA chips) of Affymetrix (USA) CytoScanTM HD Array and Clariom™ S Assay, human. Detection of expression level of seven chemosensitivity genes-RRM1, ERCC1, TOP1, TOP2a, TUBB3, TYMS, and GSTP1-was performed using PCR real-time (RT-qPCR). RESULTS: The expression of the RRM1 (AC scheme), TOP2α, TYMS, and TUBB3 genes in patients with an objective response to treatment (complete and partial regression) is higher than in patients with stabilization and progression (p < 0.05). According to our results, the presence of a high level of GSTP1 in a tumor biopsy is associated with the low efficiency of the NAC CP scheme (p = 0.05). The presence of RRM1 deletion is associated with complete and partial regression, as for the TOP1 and TUBB3 genes (p < 0.05). Higher rates of metastatic survival are associated with a high level of expression and amplification of the GSTP1 gene (log-rank test p = 0.02 and p = 0.05). CONCLUSION: Thus, a complex assessment of the chemotherapy's gene expression is important not only for understanding the heterogeneity and molecular biology of breast cancer but also to obtain a more accurate disease prognosis.

Evolutionary patterns of chromosomal instability and mismatch repair deficiency in proximal and distal colorectal cancer.

AIM: Colorectal carcinomas (CRCs) progress through heterogeneous pathways. The aim of this study was to analyse whether or not the cytogenetic evolution of CRC is linked to tumour site, level of chromosomal imbalance and metastasis. METHOD: A set of therapy-naïve pT3 CRCs comprising 26 proximal and 49 distal pT3 CRCs was studied by combining immunohistochemistry of mismatch repair (MMR) proteins, microsatellite analyses and molecular karyotyping as well as clinical parameters. RESULTS: A MMR deficient/microsatellite-unstable (dMMR/MSI-H) status was associated with location of the primary tumour proximal to the splenic flexure, and dMMR/MSI-H tumours presented with significantly lower levels of chromosomal imbalances compared with MMR proficient/microsatellite-stable (pMMR/MSS) tumours. Oncogenetic tree modelling suggested two evolutionary clusters characterized by dMMR/MSI-H and chromosomal instability (CIN), respectively, for both proximal and distal CRCs. In CIN cases, +13q, -18q and +20q were predicted as preferentially early events, and -1p, -4 -and -5q as late events. Separate oncogenetic tree models of proximal and distal cases indicated similar early events independent of tumour site. However, in cases with high CIN defined by more than 10 copy number aberrations, loss of 17p occurred earlier in cytogenetic evolution than in cases showing low to moderate CIN. Differences in the oncogenetic trees were observed for CRCs with lymph node and distant metastasis. Loss of 8p was modelled as an early event in node-positive CRC, while +7p and +8q comprised early events in CRC with distant metastasis. CONCLUSION: CRCs characterized by CIN follow multiple, interconnected genetic pathways in line with the basic 'Vogelgram' concept proposed for the progression of CRC that places the accumulation of genetic changes at centre of tumour evolution. However, the timing of specific genetic events may favour metastatic potential.

The feasibility of using liquid biopsies as a complementary assay for copy number aberration profiling in routinely collected paediatric cancer patient samples.

BACKGROUND: Paediatric tumours are often characterised by the presence of recurrent DNA copy number alterations (CNAs). These DNA copy number profiles, obtained from a tissue biopsy, can aid in the correct prognostic classification and therapeutic stratification of several paediatric cancer entities (e.g. MYCN amplification in neuroblastoma) and are part of the routine diagnostic practice. Liquid biopsies (LQBs) offer a potentially safer alternative for such invasive tumour tissue biopsies and can provide deeper insight into tumour heterogeneity. PROCEDURE: The robustness and reliability of LQB CNA analyses was evaluated. We performed retrospective CNA profiling using shallow whole-genome sequencing (sWGS) on paired plasma circulating cell-free DNA (cfDNA) and tissue DNA samples from routinely collected samples from paediatric patients (n = 128) representing different tumour entities, including osteosarcoma, Ewing sarcoma, rhabdomyosarcoma, Wilms tumour, brain tumours and neuroblastoma. RESULTS: Overall, we observed a good concordance between CNAs in tissue DNA and cfDNA. The main cause of CNA discordance was found to be low cfDNA sample quality (i.e. the ratio of cfDNA (<700 bp) and high molecular weight DNA (>700 bp)). Furthermore, CNAs were observed that were present in cfDNA and not in tissue DNA, or vice-versa. In neuroblastoma samples, no false-positives or false-negatives were identified for the detection of the prognostic marker MYCN amplification. CONCLUSION: In future prospective studies, CNA analysis on LQBs that are of sufficient quality can serve as a complementary assay for CNA analysis on tissue biopsies, as either cfDNA or tissue DNA can contain CNAs that cannot be identified in the other biomaterial.