Tangent normalization for somatic copy-number inference in cancer genome analysis

DAGIP: alleviating cell-free DNA sequencing biases with optimal transport

Cell-free DNA (cfDNA) is a rich source of biomarkers for various pathophysiological conditions. Preanalytical variables, such as the library preparation protocol or sequencing platform, are major confounders of cfDNA analysis. We present DAGIP, a novel data correction method that builds on optimal transport theory and deep learning, which explicitly corrects for the effect of such preanalytical variables and can infer technical biases. Our method improves cancer detection and copy number alteration analysis by alleviating the sources of variation that are not of biological origin. It also enhances fragmentomic analysis of cfDNA. DAGIP allows the integration of cohorts from different studies.

Prognostic and therapeutic implications of tumor-restrictive type III collagen in the breast cancer microenvironment

Collagen plays a critical role in regulating breast cancer progression and therapeutic resistance. An improved understanding of both the features and drivers of tumor-permissive and -restrictive collagen matrices are critical to improve prognostication and develop more effective therapeutic strategies. In this study, using a combination of in vitro, in vivo and bioinformatic experiments, we show that type III collagen (Col3) plays a tumor-restrictive role in human breast cancer. We demonstrate that Col3-deficient, human fibroblasts produce tumor-permissive collagen matrices that drive cell proliferation and suppress apoptosis in non-invasive and invasive breast cancer cell lines. In human triple-negative breast cancer biopsy samples, we demonstrate elevated deposition of Col3 relative to type I collagen (Col1) in non-invasive compared to invasive regions. Similarly, bioinformatics analysis of over 1000 breast cancer patient biopsies from The Cancer Genome Atlas BRCA cohort revealed that patients with higher Col3:Col1 bulk tumor expression had improved overall, disease-free, and progression-free survival relative to those with higher Col1:Col3 expression. Using an established 3D culture model, we show that Col3 increases spheroid formation and induces the formation of lumen-like structures that resemble non-neoplastic mammary acini. Finally, our in vivo study shows co-injection of murine breast cancer cells (4T1) with rhCol3-supplemented hydrogels limits tumor growth and decreases pulmonary metastatic burden compared to controls. Taken together, these data collectively support a tumor-suppressive role for Col3 in human breast cancer and suggest that strategies that increase Col3 may provide a safe and effective therapeutic modality to limit recurrence in breast cancer patients.

Conumee 2.0: enhanced copy-number variation analysis from DNA methylation arrays for humans and mice

MOTIVATION

Copy-number variations (CNVs) are common genetic alterations in cancer and their detection may impact tumor classification and therapeutic decisions. However, detection of clinically relevant large and focal CNVs remains challenging when sample material or resources are limited. This has motivated us to create a software tool to infer CNVs from DNA methylation arrays which are often generated as part of clinical routines and in research settings.

RESULTS

We present our R package, conumee 2.0, that combines tangent normalization, an adjustable genomic binning heuristic, and weighted circular binary segmentation to utilize DNA methylation arrays for CNV analysis and mitigate technical biases and batch effects. Segmentation results were validated in a lung squamous cell carcinoma dataset from TCGA (n = 367 samples) by comparison to segmentations derived from genotyping arrays (Pearson's correlation coefficient of 0.91). We further introduce a segmented block bootstrapping approach to detect focal alternations that achieved 60.9% sensitivity and 98.6% specificity for deletions affecting CDKN2A/B (60.0% and 96.9% for RB1, respectively) in a low-grade glioma cohort from TCGA (n = 239 samples). Finally, our tool provides functionality to detect and summarize CNVs across large sample cohorts.

AVAILABILITY AND IMPLEMENTATION

Conumee 2.0 is available under open-source license at: https://github.com/hovestadtlab/conumee2.

Cancer aneuploidies are shaped primarily by effects on tumour fitness

Aneuploidies-whole-chromosome or whole-arm imbalances-are the most prevalent alteration in cancer genomes1,2. However, it is still debated whether their prevalence is due to selection or ease of generation as passenger events1,2. Here we developed a method, BISCUT, that identifies loci subject to fitness advantages or disadvantages by interrogating length distributions of telomere- or centromere-bounded copy-number events. These loci were significantly enriched for known cancer driver genes, including genes not detected through analysis of focal copy-number events, and were often lineage specific. BISCUT identified the helicase-encoding gene WRN as a haploinsufficient tumour-suppressor gene on chromosome 8p, which is supported by several lines of evidence. We also formally quantified the role of selection and mechanical biases in driving aneuploidy, finding that rates of arm-level copy-number alterations are most highly correlated with their effects on cellular fitness1,2. These results provide insight into the driving forces behind aneuploidy and its contribution to tumorigenesis.

Evolutionary history of transformation from chronic lymphocytic leukemia to Richter syndrome

Richter syndrome (RS) arising from chronic lymphocytic leukemia (CLL) exemplifies an aggressive malignancy that develops from an indolent neoplasm. To decipher the genetics underlying this transformation, we computationally deconvoluted admixtures of CLL and RS cells from 52 patients with RS, evaluating paired CLL-RS whole-exome sequencing data. We discovered RS-specific somatic driver mutations (including IRF2BP2, SRSF1, B2M, DNMT3A and CCND3), recurrent copy-number alterations beyond del(9p21)(CDKN2A/B), whole-genome duplication and chromothripsis, which were confirmed in 45 independent RS cases and in an external set of RS whole genomes. Through unsupervised clustering, clonally related RS was largely distinct from diffuse large B cell lymphoma. We distinguished pathways that were dysregulated in RS versus CLL, and detected clonal evolution of transformation at single-cell resolution, identifying intermediate cell states. Our study defines distinct molecular subtypes of RS and highlights cell-free DNA analysis as a potential tool for early diagnosis and monitoring.