Shallow whole genome sequencing for robust copy number profiling of formalin-fixed paraffin-embedded breast cancers

The reckoning of chromosomal instability: past, present, future

Quantitative measures of CIN are crucial to our understanding of its role in cancer. Technological advances have changed the way CIN is quantified, offering increased accuracy and insight. Here, we review measures of CIN through its rise as a field, discuss considerations for its measurement, and look forward to future quantification of CIN.

Clonal Patterns Between Pouch Neoplasia and Prior Colorectal Neoplasia in Inflammatory Bowel Disease Patients: An Exploratory Cohort Study

Prior colorectal neoplasia is the strongest predictor of pouch neoplasia in inflammatory bowel disease, but the underlying mechanism is unknown. We observed clonality between colorectal and pouch neoplasia in 30% of patients, indicating that most pouch neoplasia develops clonally independent from prior colorectal lesions.

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

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

High levels of chromosomal aberrations negatively associate with benefit to checkpoint inhibition in NSCLC

BACKGROUND

Immune checkpoint inhibitors (ICIs) targeting the programmed cell death 1/programmed death-ligand 1 axis have transformed the management of advanced non-small cell lung cancer (NSCLC). However, many patients do not benefit from this type of treatment, and thus several molecular biomarkers of benefit have been explored. The value of somatic copy number alterations (SCNAs) burden remains elusive.

PATIENTS AND METHODS

We assembled a cohort of 109 patients with NSCLC treated with ICIs and available tumor samples. We performed shallow whole-genome sequencing on 89 patients to determine genome-wide SCNAs and targeted gene expression analysis on 63 patients to study immune infiltration. We analyzed SCNAs burden in different ways (ie, the fraction of the genome altered or number of events) and studied their association with ICIs benefit based on survival analysis. We correlated SCNAs burden and immune infiltration on 35 patients of our cohort and on patients with lung adenocarcinoma from The Cancer Genome Atlas (TCGA).

RESULTS

High SCNAs burden, computed in diverse ways, is negatively associated with ICIs progression-free survival (PFS), with the fraction of the genome altered (FGA) by arm and chromosome events showing the strongest association with PFS (p=0.002) (n=77). Nevertheless, we found differences in SCNAs across some clinicopathological features (sample site origin). A multivariate analysis adjusted for relevant characteristics showed that the FGA of arm and chromosome alterations was strongly associated with PFS (HR=2.21, p=3.3 x 10^-5). Finally, we confirmed that SCNAs burden negatively correlates with tumor immune infiltration (n=35), although this correlation was not found for the males studied. Similar results were observed in the TCGA cohort.

CONCLUSIONS

SCNAs burden is a potential biomarker of benefit to ICIs in patients with NSCLC, although there appear to be some nuances worth consideration. Further studies will be needed to establish its role as a biomarker of benefit to ICIs.

Oral Etoposide and Trastuzumab Use for HER2-Positive Metastatic Breast Cancer: A Retrospective Study from the Institut Curie Hospitals

BACKGROUND

The TOP2A and ERBB2 genes are co-amplified in about 40% of HER2 positive (HER2+) breast cancers. Oral etoposide (VP16), an inhibitor of topoisomerase-II (encoded by TOP2A), has demonstrated clinical activity in metastatic breast cancer (MBC). The benefit of oral VP16 combined with trastuzumab (VP16-T) in HER2+ MBC has not yet been evaluated.

METHODS

Patients treated at the Institut Curie Hospitals with VP16-T for HER2+ MBC were retrieved by an in silico search. Progression-free survival (PFS), overall survival (OS), response rate, prolonged PFS (defined as at least 6 months), clinical benefit, and toxicity were assessed. The co-amplification of ERBB2 and TOP2A was assessed by shallow whole genome sequencing on tumor tissue whenever available.

RESULTS

Forty-three patients received VP16-T after a median number of six prior treatment lines for HER2+ MBC. Median PFS and OS were 2.9 months (95% CI [2.4-4.7]) and 11.3 months (95% CI [8.3-25.0]), respectively. Three patients had a complete response, while 12/40 (30%) experienced clinical benefit. Only three patients stopped treatment for toxicity. Seven (35%) patients displayed a TOP2A/ERBB2 co-amplification. No statistically significant correlation was found between outcome and TOP2A/ERBB2 co-amplification.

CONCLUSION

Our analysis suggests a favorable efficacy and toxicity profile for VP16-T in patients with heavily pretreated HER2+ MBC.

Genomic alterations and evolution of cell clusters in metastatic invasive micropapillary carcinoma of the breast

Invasive micropapillary carcinoma (IMPC) has very high rates of lymphovascular invasion and lymph node metastasis and has been reported in several organs. However, the genomic mechanisms underlying its metastasis are unclear. Here, we perform whole-genome sequencing of tumor cell clusters from primary IMPC and paired axillary lymph node metastases. Cell clusters in multiple lymph node foci arise from a single subclone of the primary tumor. We find evidence that the monoclonal metastatic ancestor in primary IMPC shares high frequency copy-number loss of PRDM16 and IGSF9 and the copy number gain of ALDH2. Immunohistochemistry analysis further shows that low expression of IGSF9 and PRDM16 and high expression of ALDH2 are associated with lymph node metastasis and poor survival of patients with IMPC. We expect these genomic and evolutionary profiles to contribute to the accurate diagnosis of IMPC.

Liquid Biopsy to Detect Minimal Residual Disease: Methodology and Impact

One reason why some patients experience recurrent disease after a curative-intent treatment might be the persistence of residual tumor cells, called minimal residual disease (MRD). MRD cannot be identified by standard radiological exams or clinical evaluation. Tumor-specific alterations found in the blood indirectly diagnose the presence of MRD. Liquid biopsies thus have the potential to detect MRD, allowing, among other things, the detection of circulating tumor DNA (ctDNA), circulating tumor cells (CTC), or tumor-specific microRNA. Although liquid biopsy is increasingly studied, several technical issues still limit its clinical applicability: low sensitivity, poor standardization or reproducibility, and lack of randomized trials demonstrating its clinical benefit. Being able to detect MRD could give clinicians a more comprehensive view of the risk of relapse of their patients and could select patients requiring treatment escalation with the goal of improving cancer survival. In this review, we are discussing the different methodologies used and investigated to detect MRD in solid cancers, their respective potentials and issues, and the clinical impacts that MRD detection will have on the management of cancer patients.

Copy Number Signatures and Clinical Outcomes in Upper Tract Urothelial Carcinoma

Tumor staging of upper tract urothelial carcinomas (UTUCs) is relatively difficult to assert accurately before surgery. Here, we used copy number (CN) signatures as a tool to explore their clinical significance of molecular stratification in UTUC. CN signatures were extracted by non-negative matrix factorization from the whole-genome sequencing (WGS) data of 90 Chinese UTUC primary tumor samples. A validation UTUC cohort (n = 56) and a cohort from urinary cell-free DNA (cfDNA) of urothelial cancer patients (n = 94) and matched primary tumors were also examined. Survival analyses were measured using the Kaplan–Meier, and Cox regression was used for multivariate analysis. Here, we identified six CN signatures (Sig1–6). Patients with a high contribution of Sig6 (Sig6high) were associated with higher microsatellite instability level and papillary architecture and had a favorable outcome. Patients with a low weighted genome integrity index were associated with positive lymph node and showed the worst outcome. Sig6high was identified to be an independently prognostic factor. The predictive significance of CN signature was identified by a validation UTUC cohort. CN signatures retained great concordance between primary tumor and urinary cfDNA. In conclusion, our results reveal that CN signature assessment for risk stratification is feasible and provides a basis for clinical studies that evaluate therapeutic interventions and prognosis.

LiquidCNA: Tracking subclonal evolution from longitudinal liquid biopsies using somatic copy number alterations

Cell-free DNA (cfDNA) measured via liquid biopsies provides a way for minimally invasive monitoring of tumor evolutionary dynamics during therapy. Here we present liquidCNA, a method to track subclonal evolution from longitudinally collected cfDNA samples sequenced through cost-effective low-pass whole-genome sequencing. LiquidCNA utilizes somatic copy number alteration (SCNA) to simultaneously genotype and quantify the size of the dominant subclone without requiring B-allele frequency information, matched-normal samples, or prior knowledge on the genetic identity of the emerging clone. We demonstrate the accuracy of liquidCNA in synthetically generated sample sets and in vitro mixtures of cancer cell lines. In vivo application in patients with metastatic lung cancer reveals the progressive emergence of a novel tumor subpopulation. LiquidCNA is straightforward to use, is computationally inexpensive, and enables continuous monitoring of subclonal evolution to understand and control-therapy-induced resistance.

Evaluation of tools for identifying large copy number variations from ultra-low-coverage whole-genome sequencing data

BACKGROUND

Detection of copy number variations (CNVs) from high-throughput next-generation whole-genome sequencing (WGS) data has become a widely used research method during the recent years. However, only a little is known about the applicability of the developed algorithms to ultra-low-coverage (0.0005-0.8×) data that is used in various research and clinical applications, such as digital karyotyping and single-cell CNV detection.

RESULT

Here, the performance of six popular read-depth based CNV detection algorithms (BIC-seq2, Canvas, CNVnator, FREEC, HMMcopy, and QDNAseq) was studied using ultra-low-coverage WGS data. Real-world array- and karyotyping kit-based validation were used as a benchmark in the evaluation. Additionally, ultra-low-coverage WGS data was simulated to investigate the ability of the algorithms to identify CNVs in the sex chromosomes and the theoretical minimum coverage at which these tools can accurately function. Our results suggest that while all the methods were able to detect large CNVs, many methods were susceptible to producing false positives when smaller CNVs (< 2 Mbp) were detected. There was also significant variability in their ability to identify CNVs in the sex chromosomes. Overall, BIC-seq2 was found to be the best method in terms of statistical performance. However, its significant drawback was by far the slowest runtime among the methods (> 3 h) compared with FREEC (~ 3 min), which we considered the second-best method.

CONCLUSIONS

Our comparative analysis demonstrates that CNV detection from ultra-low-coverage WGS data can be a highly accurate method for the detection of large copy number variations when their length is in millions of base pairs. These findings facilitate applications that utilize ultra-low-coverage CNV detection.

A practical guide to cancer subclonal reconstruction from DNA sequencing

Subclonal reconstruction from bulk tumor DNA sequencing has become a pillar of cancer evolution studies, providing insight into the clonality and relative ordering of mutations and mutational processes. We provide an outline of the complex computational approaches used for subclonal reconstruction from single and multiple tumor samples. We identify the underlying assumptions and uncertainties in each step and suggest best practices for analysis and quality assessment. This guide provides a pragmatic resource for the growing user community of subclonal reconstruction methods.

Landscape of G-quadruplex DNA structural regions in breast cancer

Response and resistance to anticancer therapies vary due to intertumor and intratumor heterogeneity1. Here, we map differentially enriched G-quadruplex (G4) DNA structure-forming regions (∆G4Rs) in 22 breast cancer patient-derived tumor xenograft (PDTX) models. ∆G4Rs are associated with the promoters of highly amplified genes showing high expression, and with somatic single-nucleotide variants. Differences in ΔG4R landscapes reveal seven transcription factor programs across PDTXs. ∆G4R abundance and locations stratify PDTXs into at least three G4-based subtypes. ∆G4Rs in most PDTXs (14 of 22) were found to associate with more than one breast cancer subtype, which we also call an integrative cluster (IC)2. This suggests the frequent coexistence of multiple breast cancer states within a PDTX model, the majority of which display aggressive triple-negative IC10 gene activity. Short-term cultures of PDTX models with increased ∆G4R levels are more sensitive to small molecules targeting G4 DNA. Thus, G4 landscapes reveal additional IC-related intratumor heterogeneity in PDTX biopsies, improving breast cancer stratification and potentially identifying new treatment strategies.

Artificial intelligence (AI) and big data in cancer and precision oncology

Artificial intelligence (AI) and machine learning have significantly influenced many facets of the healthcare sector. Advancement in technology has paved the way for analysis of big datasets in a cost- and time-effective manner. Clinical oncology and research are reaping the benefits of AI. The burden of cancer is a global phenomenon. Efforts to reduce mortality rates requires early diagnosis for effective therapeutic interventions. However, metastatic and recurrent cancers evolve and acquire drug resistance. It is imperative to detect novel biomarkers that induce drug resistance and identify therapeutic targets to enhance treatment regimes. The introduction of the next generation sequencing (NGS) platforms address these demands, has revolutionised the future of precision oncology. NGS offers several clinical applications that are important for risk predictor, early detection of disease, diagnosis by sequencing and medical imaging, accurate prognosis, biomarker identification and identification of therapeutic targets for novel drug discovery. NGS generates large datasets that demand specialised bioinformatics resources to analyse the data that is relevant and clinically significant. Through these applications of AI, cancer diagnostics and prognostic prediction are enhanced with NGS and medical imaging that delivers high resolution images. Regardless of the improvements in technology, AI has some challenges and limitations, and the clinical application of NGS remains to be validated. By continuing to enhance the progression of innovation and technology, the future of AI and precision oncology show great promise.

A Novel Approach for the Genetic Analysis of Biliary Tract Cancer Specimens Obtained Through Endoscopic Ultrasound-Guided Fine Needle Aspiration Using Targeted Amplicon Sequencing

Biliary tract cancer (BTC) is an aggressive malignant tumor, and biomarker-based clinical trials for this cancer are currently ongoing. Endoscopic ultrasound–guided fine needle aspiration (EUS-FNA) is a safe procedure and enables pathological diagnoses; however, it is uncertain whether a tiny tumor sample of BTC obtained through EUS-FNA can be analyzed for diverse genetic alterations in the development and tolerance of BTC. Thus, we aimed to verify the feasibility of genetic analyses with EUS-FNA samples of BTC.

Ultra-low depth sequencing of plasma cell DNA for the detection of copy number aberrations in multiple myeloma

Cytogenetic abnormalities are powerful prognostic factors in multiple myeloma (MM) and are routinely analyzed by FISH on bone marrow (BM) plasma cells (PC). Although considered the gold standard, FISH experiments can be laborious and expensive. Therefore, array-CGH (aCGH) has been introduced as an alternative approach for detecting copy number aberrations (CNA), reducing the number of FISH experiments per case and yielding genome-wide information. Currently, next generation sequencing (NGS) technologies offer new perspectives for the diagnostic workup of malignant disorders. In this study, we examined ultra-low depth whole genome sequencing (LDS) as a valid alternative for aCGH for the detection of CNA in BM PC in MM. To this end, BM aspirates obtained in a diagnostic setting from 20 MM cases were analyzed. CD138+ cell-sorted samples were subjected to FISH analysis. DNA was extracted for subsequent aCGH and LDS analysis. CNA were detected by aCGH and LDS in all but one case. Importantly, all CNA identified by parallel first generation aCGH analysis were also detected by LDS, along with six additional CNA in five cases. One of these additional aberrations was in a region of prognostic importance in MM and was confirmed using FISH. However, risk stratification in these particular cases was unaffected. Thus, a perfectly concordant prognostication between array-CGH and LDS was observed. This validates LDS as a novel and cost-efficient tool for the detection of CNA in MM.

Detection of Copy Number Alterations by Shallow Whole-Genome Sequencing of Formalin-Fixed, Paraffin-Embedded Tumor Tissue

CONTEXT.—

In routine clinical practice, tumor tissue is stored in formalin-fixed, paraffin-embedded blocks. However, the use of formalin-fixed, paraffin-embedded tissue for genome analysis is challenged by poorer DNA quality and quantity. Although several studies have reported genome-wide massive parallel sequencing applied on formalin-fixed, paraffin-embedded samples for mutation analysis, copy number analysis is not yet commonly performed.

OBJECTIVE.—

To evaluate the use of formalin-fixed, paraffin-embedded tissue for copy number alteration detection using shallow whole-genome sequencing, more generally referred to as copy number variation sequencing.

DESIGN.—

We selected samples from 21 patients, covering a range of different tumor entities. The performance of copy number detection was compared across 3 setups: array comparative genomic hybridization in combination with fresh material; copy number variation sequencing on fresh material; and copy number variation sequencing on formalin-fixed, paraffin-embedded material.

RESULTS.—

Very similar copy number profiles between paired samples were obtained. Although formalin-fixed, paraffin-embedded profiles often displayed more noise, detected copy numbers seemed equally reliable if the tumor fraction was at least 20%.

CONCLUSIONS.—

Copy number variation sequencing of formalin-fixed, paraffin-embedded material represents a trustworthy method. It is very likely that copy number variation sequencing of routinely obtained biopsy material will become important for individual patient care and research. Moreover, the basic technology needed for copy number variation sequencing is present in most molecular diagnostics laboratories.

WisecondorX: improved copy number detection for routine shallow whole-genome sequencing

Shallow whole-genome sequencing to infer copy number alterations (CNAs) in the human genome is rapidly becoming the method par excellence for routine diagnostic use. Numerous tools exist to deduce aberrations from massive parallel sequencing data, yet most are optimized for research and often fail to redeem paramount needs in a clinical setting. Optimally, a read depth-based analytical software should be able to deal with single-end and low-coverage data-this to make sequencing costs feasible. Other important factors include runtime, applicability to a variety of analyses and overall performance. We compared the most important aspect, being normalization, across six different CNA tools, selected for their assumed ability to satisfy the latter needs. In conclusion, WISECONDOR, which uses a within-sample normalization technique, undoubtedly produced the best results concerning variance, distributional assumptions and basic ability to detect true variations. Nonetheless, as is the case with every tool, WISECONDOR has limitations, which arise through its exclusiveness for non-invasive prenatal testing. Therefore, this work presents WisecondorX in addition, an improved WISECONDOR that enables its use for varying types of applications. WisecondorX is freely available at https://github.com/CenterForMedicalGeneticsGhent/WisecondorX.

Circulating cell-free DNA for non-invasive cancer management

Cell-free DNA (cfDNA) was first identified in human plasma in 1948 and is thought to be released from cells throughout the body into the circulatory system. In cancer, a portion of the cfDNA originates from tumour cells, referred to as circulating-tumour DNA (ctDNA), and can contain mutations corresponding to the patient's tumour, for instance specific TP53 alleles. Profiling of cfDNA has recently become an area of increasing clinical relevance in oncology, in particular due to advances in the sensitivity of molecular biology techniques and development of next generation sequencing technologies, as this allows tumour mutations to be identified and tracked non-invasively. This has opened up new possibilities for monitoring tumour evolution and acquisition of resistance, as well as for guiding treatment decisions when tumour biopsy tissue is insufficient or unavailable. In this review, we will discuss the biology of cell-free nucleic acids, methods of analysis, and the potential clinical uses of these techniques, as well as the on-going clinical development of ctDNA assays.

Multiple Facets of Marine Invertebrate Conservation Genomics

Conservation genomics aims to preserve the viability of populations and the biodiversity of living organisms. Invertebrate organisms represent 95% of animal biodiversity; however, few genomic resources currently exist for the group. The subset of marine invertebrates includes the most ancient metazoan lineages and possesses codes for unique gene products and possible keys to adaptation. The benefits of supporting invertebrate conservation genomics research (e.g., likely discovery of novel genes, protein regulatory mechanisms, genomic innovations, and transposable elements) outweigh the various hurdles (rare, small, or polymorphic starting materials). Here we review best conservation genomics practices in the laboratory and in silico when applied to marine invertebrates and also showcase unique features in several case studies of acroporid corals, crown-of-thorns starfish, apple snails, and abalone. Marine conservation genomics should also address how diversity can lead to unique marine innovations, the impact of deleterious variation, and how genomic monitoring and profiling could positively affect broader conservation goals (e.g., value of baseline data for in situ/ex situ genomic stocks).