Fragmentomics features of cell-free DNA represent promising non-invasive biomarkers for cancer diagnosis. A lack of systematic evaluation of biases in feature quantification hinders the adoption of such applications. We compare features derived from whole-genome sequencing of ten healthy donors using nine library kits and ten data-processing routes and validated in 1182 plasma samples from published studies. Our results clarify the variations from library preparation and feature quantification methods. We design the Trim Align Pipeline and cfDNAPro R package as unified interfaces for data pre-processing, feature extraction, and visualization to standardize multi-modal feature engineering and integration for machine learning.

Colorectal cancer (CRC) has traditionally been thought to develop through stepwise mutation of the APC tumour suppressor and other driver genes, coupled with expansion of positively selected clones. However, recent publications show that many premalignant lesions comprise multiple clones expressing different mutant APC proteins1-4. Here, by mediating transformation on different mouse backgrounds containing mutations in Kras or other common CRC driver genes, we establish that the presence of diverse priming events in the normal mouse intestinal epithelium can change the transformation and clonal-selection landscape, permitting the fixation of strong driver mutations in Apc and Ctnnb1 that are otherwise lost due to negative selection. These findings, combined with our demonstration of mutational patterns consistent with similar priming events in human CRC, suggest that the order in which driver mutations occur in intestinal epithelium can determine whether clones are positively or negatively selected and can shape subsequent tumour development.

Antituberculosis drugs, mostly developed over 60 years ago, combined with a poorly effective vaccine, have failed to eradicate tuberculosis. More worryingly, multiresistant strains of Mycobacterium tuberculosis (MTB) are constantly emerging. Innovative strategies are thus urgently needed to improve tuberculosis treatment. Recently, host-directed therapy has emerged as a promising strategy to be used in adjunct with existing or future antibiotics, by improving innate immunity or limiting immunopathology. Here, using high-content imaging, we identified novel 1,2,4-oxadiazole-based compounds, which allow human macrophages to control MTB replication. Genome-wide gene expression analysis revealed that these molecules induced zinc remobilization inside cells, resulting in bacterial zinc intoxication. More importantly, we also demonstrated that, upon treatment with these novel compounds, MTB became even more sensitive to antituberculosis drugs, in vitro and in vivo, in a mouse model of tuberculosis. Manipulation of heavy metal homeostasis holds thus great promise to be exploited to develop host-directed therapeutic interventions.

The estrogen receptor-α (ER) drives 75% of breast cancers. On activation, the ER recruits and assembles a 1-2 MDa transcriptionally active complex. These complexes can modulate tumour growth, and understanding the roles of individual proteins within these complexes can help identify new therapeutic targets. Here, we present the discovery of ER and ZMIZ1 within the same multi-protein assembly by quantitative proteomics, and validated by proximity ligation assay. We characterise ZMIZ1 function by demonstrating a significant decrease in the proliferation of ER-positive cancer cell lines. To establish a role for the ER-ZMIZ1 interaction, we measured the transcriptional changes in the estrogen response post-ZMIZ1 knockdown using an RNA-seq time-course over 24 h. Gene set enrichment analysis of the ZMIZ1-knockdown data identified a specific delay in the response of estradiol-induced cell cycle genes. Integration of ENCODE data with our RNA-seq results identified that ER and ZMIZ1 both bind the promoter of E2F2. We therefore propose that ER and ZMIZ1 interact to enable the efficient estrogenic response at subset of cell cycle genes via a novel ZMIZ1-ER-E2F2 signalling axis. Finally, we show that high ZMIZ1 expression is predictive of worse patient outcome, ER and ZMIZ1 are co-expressed in breast cancer patients in TCGA and METABRIC, and the proteins are co-localised within the nuclei of tumour cell in patient biopsies. In conclusion, we establish that ZMIZ1 is a regulator of the estrogenic cell cycle response and provide evidence of the biological importance of the ER-ZMIZ1 interaction in ER-positive patient tumours, supporting potential clinical relevance.

Loss-of-function mutations in the tumour suppressor APC are an initial step in intestinal tumorigenesis1,2. APC-mutant intestinal stem cells outcompete their wild-type neighbours through the secretion of Wnt antagonists, which accelerates the fixation and subsequent rapid clonal expansion of mutants3-5. Reports of polyclonal intestinal tumours in human patients and mouse models appear at odds with this process6,7. Here we combine multicolour lineage tracing with chemical mutagenesis in mice to show that a large proportion of intestinal tumours have a multiancestral origin. Polyclonal tumours retain a structure comprising subclones with distinct Apc mutations and transcriptional states, driven predominantly by differences in KRAS and MYC signalling. These pathway-level changes are accompanied by profound differences in cancer stem cell phenotypes. Of note, these findings are confirmed by introducing an oncogenic Kras mutation that results in predominantly monoclonal tumour formation. Further, polyclonal tumours have accelerated growth dynamics, suggesting a link between polyclonality and tumour progression. Together, these findings demonstrate the role of interclonal interactions in promoting tumorigenesis through non-cell autonomous pathways that are dependent on the differential activation of oncogenic pathways between clones.

Anti-tuberculosis drugs, mostly developed over 60 years ago, combined with a poorly effective vaccine, have failed to eradicate tuberculosis. More worryingly, multi-resistant strains of Mycobacterium tuberculosis are constantly emerging. Innovative strategies are thus urgently needed to improve tuberculosis treatment. Recently, host-directed therapy has emerged as a promising strategy to be used in adjunct with existing or future antibiotics, by improving innate immunity or limiting immunopathology. Here, using high content imaging, we identified novel 1,2,4-oxadiazole-based compounds, that allow human macrophages to control MTB replication. Genome-wide gene expression analysis revealed that these molecules induced zinc remobilization inside cells, resulting in bacterial zinc intoxication. More importantly, we also demonstrated that, upon treatment with these novel compounds, M. tuberculosis became even more sensitive to anti-tuberculosis drugs, in vitro and in vivo , in a mouse model of tuberculosis. Manipulation of heavy metal homeostasis holds thus great promise to be exploited to develop host-directed therapeutic interventions.

High-grade serous ovarian carcinoma (HGSOC) is characterised by poor outcome and extreme chromosome instability (CIN). Therapies targeting centrosome amplification (CA), a key mediator of chromosome missegregation, may have significant clinical utility in HGSOC. However, the prevalence of CA in HGSOC, its relationship to genomic biomarkers of CIN and its potential impact on therapeutic response have not been defined. Using high-throughput multi-regional microscopy on 287 clinical HGSOC tissues and 73 cell lines models, here we show that CA through centriole overduplication is a highly recurrent and heterogeneous feature of HGSOC and strongly associated with CIN and genome subclonality. Cell-based studies showed that high-prevalence CA is phenocopied in ovarian cancer cell lines, and that high CA is associated with increased multi-treatment resistance; most notably to paclitaxel, the commonest treatment used in HGSOC. CA in HGSOC may therefore present a potential driver of tumour evolution and a powerful biomarker for response to standard-of-care treatment.

Abstract NF-κB driven cellular immunity is essential for both pro- and anti-inflammatory responses to microbes, which makes it one of the most frequently targeted pathways by bacteria during pathogenesis. How NF-κB tunes the epithelial response to Streptococcus pneumoniae across the spectrum of commensal to pathogenic outcomes is not fully understood. In this study, we compare a commensal-like 6B ST90 strain to an invasive TIGR4 isolate and demonstrate that TIGR4 both blunts and antagonizes NF-κB activation. We identified, through comparative mass spectrometry of the p65 interactome, that the 6B ST90 isolate drives a non-canonical NF-κB RelB cascade, whereas TIGR4 induces p65 degradation though aggrephagy. Mechanistically, TIGR4 challenge triggers a novel interaction of COMMD2 with p65 and p62. This complex mediates export of p65 for degradation and is necessary for repressing host inflammatory transcription. With these results, we reveal for the first time a new bacterial pathogenesis mechanism to repress host inflammatory response though COMMD2 mediated turnover of p65, and present a paradigm for diverging NF-κB responses to pneumococcus.

AbstractAnalysis of circulating tumor DNA (ctDNA) to monitor cancer dynamics and detect minimal residual disease has been an area of increasing interest. Multiple methods have been proposed but few studies have compared the performance of different approaches. Here, we compare detection of ctDNA in serial plasma samples from patients with breast cancer using different tumor‐informed and tumor‐naïve assays designed to detect structural variants (SVs), single nucleotide variants (SNVs), and/or somatic copy‐number aberrations, by multiplex PCR, hybrid capture, and different depths of whole‐genome sequencing. Our results demonstrate that the ctDNA dynamics and allele fractions (AFs) were highly concordant when analyzing the same patient samples using different assays. Tumor‐informed assays showed the highest sensitivity for detection of ctDNA at low concentrations. Hybrid capture sequencing targeting between 1,347 and 7,491 tumor‐identified mutations at high depth was the most sensitive assay, detecting ctDNA down to an AF of 0.00024% (2.4 parts per million, ppm). Multiplex PCR targeting 21–47 tumor‐identified SVs per patient detected ctDNA down to 0.00047% AF (4.7 ppm) and has potential as a clinical assay.

The drivers of recurrence and resistance in ovarian high grade serous carcinoma remain unclear. We investigate the acquisition of resistance by collecting tumour biopsies from a cohort of 276 women with relapsed ovarian high grade serous carcinoma in the BriTROC-1 study. Panel sequencing shows close concordance between diagnosis and relapse, with only four discordant cases. There is also very strong concordance in copy number between diagnosis and relapse, with no significant difference in purity, ploidy or focal somatic copy number alterations, even when stratified by platinum sensitivity or prior chemotherapy lines. Copy number signatures are strongly correlated with immune cell infiltration, whilst diagnosis samples from patients with primary platinum resistance have increased rates of CCNE1 and KRAS amplification and copy number signature 1 exposure. Our data show that the ovarian high grade serous carcinoma genome is remarkably stable between diagnosis and relapse and acquired chemotherapy resistance does not select for common copy number drivers. 'Treatment resistance is common in ovarian high grade serous carcinoma, often leading to relapse. Here, the authors leverage shallow whole genome and panel sequencing of 276 patients with available diagnostic and relapse samples and show high concordance of copy number and mutation status.

The volume of DNA in a sequencing experiment is often amplified by PCR, leading to the possibility that the same original DNA fragment will be sequenced twice-a "PCR duplicate". Sometimes indistinguishable from these are multiple sequences arising from identical but independent molecules, which can lead to an over-estimation of the PCR duplicate proportion. The PCR duplicate proportion, and other measures derived from it, are important statistics for quality assurance, experimental design, and interpretation of sequencing experiments. Here, we provide a full likelihood basis for a combinatorial approach using heterozygous SNPs as implemented in our R package and demonstrate the efficacy of the approach. We also discuss the association with DNA copy number and demonstrate the impact on a question of inferring mitochondrial DNA copy number that has recently been a feature of several high-profile cancer studies. This is explored through a simulation study.

Dataset containing pre-processed data files required to replicate analysis performed in the publication "The genomic landscape of recurrent ovarian high grade serous carcinoma: the BriTROC-1 study" (Smith & Bradley et al. 2023).

Oesophageal adenocarcinoma (OAC) provides an ideal case study to characterize large-scale rearrangements. Using whole genome short-read sequencing of 383 cases, for which 214 had matched whole transcriptomes, we observed structural variations (SV) with a predominance of deletions, tandem duplications and inter-chromosome junctions that could be identified as LINE-1 mobile element (ME) insertions. Complex clusters of rearrangements resembling breakage-fusion-bridge cycles or extrachromosomal circular DNA accounted for 22% of complex SVs affecting known oncogenes. Counting SV events affecting known driver genes substantially increased the recurrence rates of these drivers. After excluding fragile sites, we identified 51 candidate new drivers in genomic regions disrupted by SVs, including ETV5, KAT6B and CLTC. RUNX1 was the most recurrently altered gene (24%), with many deletions inactivating the RUNT domain but preserved the reading frame, suggesting an altered protein product. These findings underscore the importance of identification of SV events in OAC with implications for targeted therapies.

Chromosomal instability is a major challenge to patient stratification and targeted drug development for high-grade serous ovarian carcinoma (HGSOC). Here we show that somatic copy number alterations (SCNAs) in frequently amplified HGSOC cancer genes significantly correlate with gene expression and methylation status. We identify five prevalent clonal driver SCNAs (chromosomal amplifications encompassing MYC, PIK3CA, CCNE1, KRAS and TERT) from multi-regional HGSOC data and reason that their strong selection should prioritise them as key biomarkers for targeted therapies. We use primary HGSOC spheroid models to test interactions between in vitro targeted therapy and SCNAs. MYC chromosomal copy number is associated with in-vitro and clinical response to paclitaxel and in-vitro response to mTORC1/2 inhibition. Activation of the mTOR survival pathway in the context of MYC-amplified HGSOC is statistically associated with increased prevalence of SCNAs in genes from the PI3K pathway. Co-occurrence of amplifications in MYC and genes from the PI3K pathway is independently observed in squamous lung cancer and triple negative breast cancer. In this work, we show that identifying co-occurrence of clonal driver SCNA genes could be used to tailor therapeutics for precision medicine.

PIK3CA mutations are the most common in breast cancer, particularly in the estrogen receptor-positive cohort, but the benefit of PI3K inhibitors has had limited success compared with approaches targeting other less common mutations. We found a frequent allelic expression imbalance between the missense mutant and wild-type PIK3CA alleles in breast tumors from the METABRIC (70.2%) and the TCGA (60.1%) projects. When considering the mechanisms controlling allelic expression, 27.7% and 11.8% of tumors showed imbalance due to regulatory variants in cis, in the two studies respectively. Furthermore, preferential expression of the mutant allele due to cis-regulatory variation is associated with poor prognosis in the METABRIC tumors ( P = 0.031). Interestingly, ER−, PR−, and HER2+ tumors showed significant preferential expression of the mutated allele in both datasets. Our work provides compelling evidence to support the clinical utility of PIK3CA allelic expression in breast cancer in identifying patients of poorer prognosis, and those with low expression of the mutated allele, who will unlikely benefit from PI3K inhibitors. Furthermore, our work proposes a model of differential regulation of a critical cancer-promoting gene in breast cancer.

For many intracellular bacterial pathogens manipulating host cell survival is essential for maintaining their replicative niche, and is a common strategy used to promote infection. The bacterial pathogen Listeria monocytogenes is well known to hijack host machinery for its own benefit, such as targeting the host histone H3 for modification by SIRT2. However, by what means this modification benefits infection, as well as the molecular players involved, were unknown. Here we show that SIRT2 activity supports Listeria intracellular survival by maintaining genome integrity and host cell viability. This protective effect is dependent on H3K18 deacetylation, which safeguards the host genome by counteracting infection-induced DNA damage. Mechanistically, infection causes SIRT2 to interact with the nucleic acid binding protein TDP-43 and localise to genomic R-loops, where H3K18 deacetylation occurs. This work highlights novel functions of TDP-43 and R-loops during bacterial infection and identifies the mechanism through which L. monocytogenes co-opts SIRT2 to allow efficient infection.

Background NRG1 gene fusions may be clinically actionable, since cancers carrying the fusion transcripts can be sensitive to tyrosine kinase inhibitors. The NRG1 gene encodes ligands for the HER2(ERBB2)-ERBB3 heterodimeric receptor tyrosine kinase, and the gene fusions are thought to lead to autocrine stimulation of the receptor. The NRG1 fusion expressed in the breast cancer cell line MDA-MB-175 serves as a model example of such fusions, showing the proposed autocrine loop and exceptional drug sensitivity. However, its structure has not been properly characterised, its oncogenic activity has not been fully explained, and there is limited data on such fusions in breast cancer. Methods We analysed genomic rearrangements and transcripts of NRG1 in MDA-MB-175 and a panel of 571 breast cancers. Results We found that the MDA-MB-175 fusion—originally reported as a DOC4(TENM4)-NRG1 fusion, lacking the cytoplasmic tail of NRG1 —is in reality a double fusion, PPP6R3-TENM4-NRG1 , producing multiple transcripts, some of which include the cytoplasmic tail. We hypothesise that many NRG1 fusions may be oncogenic not for lacking the cytoplasmic domain but because they do not encode NRG1’s nuclear-localised form. The fusion in MDA-MB-175 is the result of a very complex genomic rearrangement, which we partially characterised, that creates additional expressed gene fusions, RSF1-TENM4 , TPCN2-RSF1 , and MRPL48-GAB2 . We searched for NRG1 rearrangements in 571 breast cancers subjected to genome sequencing and transcriptome sequencing and found four cases (0.7%) with fusions, WRN-NRG1 , FAM91A1-NRG1 , ARHGEF39-NRG1 , and ZNF704-NRG1 , all splicing into NRG1 at the same exon as in MDA-MB-175. However, the WRN-NRG1 and ARHGEF39-NRG1 fusions were out of frame. We identified rearrangements of NRG1 in many more (8% of) cases that seemed more likely to inactivate than to create activating fusions, or whose outcome could not be predicted because they were complex, or both. This is not surprising because NRG1 can be pro-apoptotic and is in activated in some breast cancers. Conclusions Our results highlight the complexity of rearrangements of NRG1 in breast cancers and confirm that some do not activate but inactivate. Careful interpretation of NRG1 rearrangements will therefore be necessary for appropriate patient management.

Low-coverage or shallow whole genome sequencing (sWGS) approaches can efficiently detect somatic copy number aberrations (SCNAs) at low cost. This is clinically important for many cancers, in particular cancers with severe chromosomal instability (CIN) that frequently lack actionable point mutations and are characterised by poor disease outcome. Absolute copy number (ACN), measured in DNA copies per cancer cell, is required for meaningful comparisons between copy number states, but is challenging to estimate and in practice often requires manual curation. Using a total of 60 cancer cell lines, 148 patient-derived xenograft (PDX) and 142 clinical tissue samples, we evaluate the performance of available tools for obtaining ACN from sWGS. We provide a validated and refined tool called Rascal ( r elative to a bsolute copy number scal ing) that provides improved fitting algorithms and enables interactive visualisation of copy number profiles. These approaches are highly applicable to both pre-clinical and translational research studies on SCNA-driven cancers and provide more robust ACN fits from sWGS data than currently available tools.

Additional file 2: Supplementary Table 1. Match of fusion cDNA AF009227, originally reported by Schaefer et al. [11], to Reference genome. Supplementary Table 2. Exons of PPP6R3, TENM4, NRG1. Supplementary Table 3. Fusions of NRG1 in 571 breast cancers. Supplementary Table 4. All NRG1 rearrangements in first 250 breast cancers.

Glioma‐derived cell‐free DNA (cfDNA) is challenging to detect using liquid biopsy because quantities in body fluids are low. We determined the glioma‐derived DNA fraction in cerebrospinal fluid (CSF), plasma, and urine samples from patients using sequencing of personalized capture panels guided by analysis of matched tumor biopsies. By sequencing cfDNA across thousands of mutations, identified individually in each patient’s tumor, we detected tumor‐derived DNA in the majority of CSF (7/8), plasma (10/12), and urine samples (10/16), with a median tumor fraction of 6.4 × 10−3, 3.1 × 10−5, and 4.7 × 10−5, respectively. We identified a shift in the size distribution of tumor‐derived cfDNA fragments in these body fluids. We further analyzed cfDNA fragment sizes using whole‐genome sequencing, in urine samples from 35 glioma patients, 27 individuals with non‐malignant brain disorders, and 26 healthy individuals. cfDNA in urine of glioma patients was significantly more fragmented compared to urine from patients with non‐malignant brain disorders (P = 1.7 × 10−2) and healthy individuals (P = 5.2 × 10−9). Machine learning models integrating fragment length could differentiate urine samples from glioma patients (AUC = 0.80–0.91) suggesting possibilities for truly non‐invasive cancer detection.