Author Correction: The evolutionary history of lethal metastatic prostate cancer.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The extracellular heparan sulfatase SULF2 limits myeloid IFNß signaling and Th17 responses in inflammatory arthritis.

Heparan sulfate (HS) proteoglycans are important regulators of cellular responses to soluble mediators such as chemokines, cytokines and growth factors. We profiled changes in expression of genes encoding HS core proteins, biosynthesis enzymes and modifiers during macrophage polarisation, and found that the most highly regulated gene was Sulf2, an extracellular HS 6-O-sulfatase that was markedly downregulated in response to pro-inflammatory stimuli. We then generated Sulf2+/- bone marrow chimeric mice and examined inflammatory responses in antigen-induced arthritis, as a model of rheumatoid arthritis. Resolution of inflammation was impaired in myeloid Sulf2+/- chimeras, with elevated joint swelling and increased abundance of pro-arthritic Th17 cells in synovial tissue. Transcriptomic and in vitro analyses indicated that Sulf2 deficiency increased type I interferon signaling in bone marrow-derived macrophages, leading to elevated expression of the Th17-inducing cytokine IL6. This establishes that dynamic remodeling of HS by Sulf2 limits type I interferon signaling in macrophages, and so protects against Th17-driven pathology.

APOBEC3 mutational signatures are associated with extensive and diverse genomic instability across multiple tumour types.

BACKGROUND: The APOBEC3 (apolipoprotein B mRNA editing enzyme catalytic polypeptide 3) family of cytidine deaminases is responsible for two mutational signatures (SBS2 and SBS13) found in cancer genomes. APOBEC3 enzymes are activated in response to viral infection, and have been associated with increased mutation burden and TP53 mutation. In addition to this, it has been suggested that APOBEC3 activity may be responsible for mutations that do not fall into the classical APOBEC3 signatures (SBS2 and SBS13), through generation of double strand breaks.Previous work has mainly focused on the effects of APOBEC3 within individual tumour types using exome sequencing data. Here, we use whole genome sequencing data from 2451 primary tumours from 39 different tumour types in the Pan-Cancer Analysis of Whole Genomes (PCAWG) data set to investigate the relationship between APOBEC3 and genomic instability (GI). RESULTS AND CONCLUSIONS: We found that the number of classical APOBEC3 signature mutations correlates with increased mutation burden across different tumour types. In addition, the number of APOBEC3 mutations is a significant predictor for six different measures of GI. Two GI measures (INDELs attributed to INDEL signatures ID6 and ID8) strongly suggest the occurrence and error prone repair of double strand breaks, and the relationship between APOBEC3 mutations and GI remains when SNVs attributed to kataegis are excluded.We provide evidence that supports a model of cancer genome evolution in which APOBEC3 acts as a causative factor in the development of diverse and widespread genomic instability through the generation of double strand breaks. This has important implications for treatment approaches for cancers that carry APOBEC3 mutations, and challenges the view that APOBECs only act opportunistically at sites of single stranded DNA.

The architecture of clonal expansions in morphologically normal tissue from cancerous and non-cancerous prostates.

BACKGROUND: Up to 80% of cases of prostate cancer present with multifocal independent tumour lesions leading to the concept of a field effect present in the normal prostate predisposing to cancer development. In the present study we applied Whole Genome DNA Sequencing (WGS) to a group of morphologically normal tissue (n = 51), including benign prostatic hyperplasia (BPH) and non-BPH samples, from men with and men without prostate cancer. We assess whether the observed genetic changes in morphologically normal tissue are linked to the development of cancer in the prostate. RESULTS: Single nucleotide variants (P = 7.0 × 10-03, Wilcoxon rank sum test) and small insertions and deletions (indels, P = 8.7 × 10-06) were significantly higher in morphologically normal samples, including BPH, from men with prostate cancer compared to those without. The presence of subclonal expansions under selective pressure, supported by a high level of mutations, were significantly associated with samples from men with prostate cancer (P = 0.035, Fisher exact test). The clonal cell fraction of normal clones was always higher than the proportion of the prostate estimated as epithelial (P = 5.94 × 10-05, paired Wilcoxon signed rank test) which, along with analysis of primary fibroblasts prepared from BPH specimens, suggests a stromal origin. Constructed phylogenies revealed lineages associated with benign tissue that were completely distinct from adjacent tumour clones, but a common lineage between BPH and non-BPH morphologically normal tissues was often observed. Compared to tumours, normal samples have significantly less single nucleotide variants (P = 3.72 × 10-09, paired Wilcoxon signed rank test), have very few rearrangements and a complete lack of copy number alterations. CONCLUSIONS: Cells within regions of morphologically normal tissue (both BPH and non-BPH) can expand under selective pressure by mechanisms that are distinct from those occurring in adjacent cancer, but that are allied to the presence of cancer. Expansions, which are probably stromal in origin, are characterised by lack of recurrent driver mutations, by almost complete absence of structural variants/copy number alterations, and mutational processes similar to malignant tissue. Our findings have implications for treatment (focal therapy) and early detection approaches.

The evolutionary history of lethal metastatic prostate cancer.

Cancers emerge from an ongoing Darwinian evolutionary process, often leading to multiple competing subclones within a single primary tumour. This evolutionary process culminates in the formation of metastases, which is the cause of 90% of cancer-related deaths. However, despite its clinical importance, little is known about the principles governing the dissemination of cancer cells to distant organs. Although the hypothesis that each metastasis originates from a single tumour cell is generally supported, recent studies using mouse models of cancer demonstrated the existence of polyclonal seeding from and interclonal cooperation between multiple subclones. Here we sought definitive evidence for the existence of polyclonal seeding in human malignancy and to establish the clonal relationship among different metastases in the context of androgen-deprived metastatic prostate cancer. Using whole-genome sequencing, we characterized multiple metastases arising from prostate tumours in ten patients. Integrated analyses of subclonal architecture revealed the patterns of metastatic spread in unprecedented detail. Metastasis-to-metastasis spread was found to be common, either through de novo monoclonal seeding of daughter metastases or, in five cases, through the transfer of multiple tumour clones between metastatic sites. Lesions affecting tumour suppressor genes usually occur as single events, whereas mutations in genes involved in androgen receptor signalling commonly involve multiple, convergent events in different metastases. Our results elucidate in detail the complex patterns of metastatic spread and further our understanding of the development of resistance to androgen-deprivation therapy in prostate cancer.

Development of a multivariable risk model integrating urinary cell DNA methylation and cell-free RNA data for the detection of significant prostate cancer.

BACKGROUND: Prostate cancer exhibits severe clinical heterogeneity and there is a critical need for clinically implementable tools able to precisely and noninvasively identify patients that can either be safely removed from treatment pathways or those requiring further follow up. Our objectives were to develop a multivariable risk prediction model through the integration of clinical, urine-derived cell-free messenger RNA (cf-RNA) and urine cell DNA methylation data capable of noninvasively detecting significant prostate cancer in biopsy naïve patients. METHODS: Post-digital rectal examination urine samples previously analyzed separately for both cellular methylation and cf-RNA expression within the Movember GAP1 urine biomarker cohort were selected for a fully integrated analysis (n = 207). A robust feature selection framework, based on bootstrap resampling and permutation, was utilized to find the optimal combination of clinical and urinary markers in a random forest model, deemed ExoMeth. Out-of-bag predictions from ExoMeth were used for diagnostic evaluation in men with a clinical suspicion of prostate cancer (PSA ≥ 4 ng/mL, adverse digital rectal examination, age, or lower urinary tract symptoms). RESULTS: As ExoMeth risk score (range, 0-1) increased, the likelihood of high-grade disease being detected on biopsy was significantly greater (odds ratio = 2.04 per 0.1 ExoMeth increase, 95% confidence interval [CI]: 1.78-2.35). On an initial TRUS biopsy, ExoMeth accurately predicted the presence of Gleason score ≥3 + 4, area under the receiver-operator characteristic curve (AUC) = 0.89 (95% CI: 0.84-0.93) and was additionally capable of detecting any cancer on biopsy, AUC = 0.91 (95% CI: 0.87-0.95). Application of ExoMeth provided a net benefit over current standards of care and has the potential to reduce unnecessary biopsies by 66% when a risk threshold of 0.25 is accepted. CONCLUSION: Integration of urinary biomarkers across multiple assay methods has greater diagnostic ability than either method in isolation, providing superior predictive ability of biopsy outcomes. ExoMeth represents a more holistic view of urinary biomarkers and has the potential to result in substantial changes to how patients suspected of harboring prostate cancer are diagnosed.

A novel stratification framework for predicting outcome in patients with prostate cancer.

BACKGROUND: Unsupervised learning methods, such as Hierarchical Cluster Analysis, are commonly used for the analysis of genomic platform data. Unfortunately, such approaches ignore the well-documented heterogeneous composition of prostate cancer samples. Our aim is to use more sophisticated analytical approaches to deconvolute the structure of prostate cancer transcriptome data, providing novel clinically actionable information for this disease. METHODS: We apply an unsupervised model called Latent Process Decomposition (LPD), which can handle heterogeneity within individual cancer samples, to genome-wide expression data from eight prostate cancer clinical series, including 1,785 malignant samples with the clinical endpoints of PSA failure and metastasis. RESULTS: We show that PSA failure is correlated with the level of an expression signature called DESNT (HR = 1.52, 95% CI = [1.36, 1.7], P = 9.0 × 10-14, Cox model), and that patients with a majority DESNT signature have an increased metastatic risk (X2 test, P = 0.0017, and P = 0.0019). In addition, we develop a stratification framework that incorporates DESNT and identifies three novel molecular subtypes of prostate cancer. CONCLUSIONS: These results highlight the importance of using more complex approaches for the analysis of genomic data, may assist drug targeting, and have allowed the construction of a nomogram combining DESNT with other clinical factors for use in clinical management.

A four-group urine risk classifier for predicting outcomes in patients with prostate cancer.

OBJECTIVES: To develop a risk classifier using urine-derived extracellular vesicle (EV)-RNA capable of providing diagnostic information on disease status prior to biopsy, and prognostic information for men on active surveillance (AS). PATIENTS AND METHODS: Post-digital rectal examination urine-derived EV-RNA expression profiles (n = 535, multiple centres) were interrogated with a curated NanoString panel. A LASSO-based continuation ratio model was built to generate four prostate urine risk (PUR) signatures for predicting the probability of normal tissue (PUR-1), D'Amico low-risk (PUR-2), intermediate-risk (PUR-3), and high-risk (PUR-4) prostate cancer. This model was applied to a test cohort (n = 177) for diagnostic evaluation, and to an AS sub-cohort (n = 87) for prognostic evaluation. RESULTS: Each PUR signature was significantly associated with its corresponding clinical category (P < 0.001). PUR-4 status predicted the presence of clinically significant intermediate- or high-risk disease (area under the curve = 0.77, 95% confidence interval [CI] 0.70-0.84). Application of PUR provided a net benefit over current clinical practice. In an AS sub-cohort (n = 87), groups defined by PUR status and proportion of PUR-4 had a significant association with time to progression (interquartile range hazard ratio [HR] 2.86, 95% CI 1.83-4.47; P < 0.001). PUR-4, when used continuously, dichotomized patient groups with differential progression rates of 10% and 60% 5 years after urine collection (HR 8.23, 95% CI 3.26-20.81; P < 0.001). CONCLUSION: Urine-derived EV-RNA can provide diagnostic information on aggressive prostate cancer prior to biopsy, and prognostic information for men on AS. PUR represents a new and versatile biomarker that could result in substantial alterations to current treatment of patients with prostate cancer.

Benzo pyrene-induced DNA adducts and gene expression profiles in target and non-target organs for carcinogenesis in mice.

BACKGROUND: Gene expression changes induced by carcinogens may identify differences in molecular function between target and non-target organs. Target organs for benzo[a]pyrene (BaP) carcinogenicity in mice (lung, spleen and forestomach) and three non-target organs (liver, colon and glandular stomach) were investigated for DNA adducts by 32P-postlabelling, for gene expression changes by cDNA microarray and for miRNA expression changes by miRNA microarray after exposure of animals to BaP. RESULTS: BaP-DNA adduct formation occurred in all six organs at levels that did not distinguish between target and non-target. cDNA microarray analysis showed a variety of genes modulated significantly by BaP in the six organs and the overall gene expression patterns were tissue specific. Gene ontology analysis also revealed that BaP-induced bioactivities were tissue specific; eight genes (Tubb5, Fos, Cdh1, Cyp1a1, Apc, Myc, Ctnnb1 and Cav) showed significant expression difference between three target and three non-target organs. Additionally, several gene expression changes, such as in Trp53 activation and Stat3 activity suggested some similarities in molecular mechanisms in two target organs (lung and spleen), which were not found in the other four organs. Changes in miRNA expression were generally tissue specific, involving, in total, 21/54 miRNAs significantly up- or down-regulated. CONCLUSIONS: Altogether, these findings showed that DNA adduct levels and early gene expression changes did not fully distinguish target from non-target organs. However, mechanisms related to early changes in p53, Stat3 and Wnt/ß-catenin pathways may play roles in defining BaP organotropism.

Applications of Urinary Extracellular Vesicles in the Diagnosis and Active Surveillance of Prostate Cancer.

Prostate cancer is the most common non-cutaneous cancer among men in the UK, causing significant health and economic burdens. Diagnosis and risk prognostication can be challenging due to the genetic and clinical heterogeneity of prostate cancer as well as uncertainties in our knowledge of the underlying biology and natural history of disease development. Urinary extracellular vesicles (EVs) are microscopic, lipid bilayer defined particles released by cells that carry a variety of molecular cargoes including nucleic acids, proteins and other molecules. Urine is a plentiful source of prostate-derived EVs. In this narrative review, we summarise the evidence on the function of urinary EVs and their applications in the evolving field of prostate cancer diagnostics and active surveillance. EVs are implicated in the development of all hallmarks of prostate cancer, and this knowledge has been applied to the development of multiple diagnostic tests, which are largely based on RNA and miRNA. Common gene probes included in multi-probe tests include PCA3 and ERG, and the miRNAs miR-21 and miR-141. The next decade will likely bring further improvements in the diagnostic accuracy of biomarkers as well as insights into molecular biological mechanisms of action that can be translated into opportunities in precision uro-oncology.

Cancer invasion and anaerobic bacteria: new insights into mechanisms.

There is growing evidence that altered microbiota abundance of a range of specific anaerobic bacteria are associated with cancer, including Peptoniphilus spp., Porphyromonas spp., Fusobacterium spp., Fenollaria spp., Prevotella spp., Sneathia spp., Veillonella spp. and Anaerococcus spp. linked to multiple cancer types. In this review we explore these pathogenic associations. The mechanisms by which bacteria are known or predicted to interact with human cells are reviewed and we present an overview of the interlinked mechanisms and hypotheses of how multiple intracellular anaerobic bacterial pathogens may act together to cause host cell and tissue microenvironment changes associated with carcinogenesis and cancer cell invasion. These include combined effects on changes in cell signalling, DNA damage, cellular metabolism and immune evasion. Strategies for early detection and eradication of anaerobic cancer-associated bacterial pathogens that may prevent cancer progression are proposed.

Genomic evolution shapes prostate cancer disease type.

The development of cancer is an evolutionary process involving the sequential acquisition of genetic alterations that disrupt normal biological processes, enabling tumor cells to rapidly proliferate and eventually invade and metastasize to other tissues. We investigated the genomic evolution of prostate cancer through the application of three separate classification methods, each designed to investigate a different aspect of tumor evolution. Integrating the results revealed the existence of two distinct types of prostate cancer that arise from divergent evolutionary trajectories, designated as the Canonical and Alternative evolutionary disease types. We therefore propose the evotype model for prostate cancer evolution wherein Alternative-evotype tumors diverge from those of the Canonical-evotype through the stochastic accumulation of genetic alterations associated with disruptions to androgen receptor DNA binding. Our model unifies many previous molecular observations, providing a powerful new framework to investigate prostate cancer disease progression.

Caution regarding the specificities of pan-cancer microbial structure.

Results published in an article by Poore et al. (Nature. 2020;579:567-574) suggested that machine learning models can almost perfectly distinguish between tumour types based on their microbial composition using machine learning models. Whilst we believe that there is the potential for microbial composition to be used in this manner, we have concerns with the paper that make us question the certainty of the conclusions drawn. We believe there are issues in the areas of the contribution of contamination, handling of batch effects, false positive classifications and limitations in the machine learning approaches used. This makes it difficult to identify whether the authors have identified true biological signal and how robust these models would be in use as clinical biomarkers. We commend Poore et al. on their approach to open data and reproducibility that has enabled this analysis. We hope that this discourse assists the future development of machine learning models and hypothesis generation in microbiome research.

Quality Control in Metagenomics Data.

Experiments involving metagenomics data are become increasingly commonplace. Processing such data requires a unique set of considerations. Quality control of metagenomics data is critical to extracting pertinent insights. In this chapter, we outline some considerations in terms of study design and other confounding factors that can often only be realized at the point of data analysis.In this chapter, we outline some basic principles of quality control in metagenomics, including overall reproducibility and some good practices to follow. The general quality control of sequencing data is then outlined, and we introduce ways to process this data by using bash scripts and developing pipelines in Snakemake (Python).A significant part of quality control in metagenomics is in analyzing the data to ensure you can spot relationships between variables and to identify when they might be confounded. This chapter provides a walkthrough of analyzing some microbiome data (in the R statistical language) and demonstrates a few days to identify overall differences and similarities in microbiome data. The chapter is concluded by discussing remarks about considering taxonomic results in the context of the study and interrogating sequence alignments using the command line.

Gene-Transcript Expression in Urine Supernatant and Urine Cell-Sediment Are Different but Equally Useful for Detecting Prostate Cancer.

There is considerable interest in urine as a non-invasive liquid biopsy to detect prostate cancer (PCa). PCa-specific transcripts such as the TMPRSS2:ERG fusion gene can be found in both urine extracellular vesicles (EVs) and urine cell-sediment (Cell) but the relative usefulness of these and other genes in each fraction in PCa detection has not been fully elucidated. Urine samples from 76 men (PCa n = 40, non-cancer n = 36) were analysed by NanoString for 154 PCa-associated genes-probes, 11 tissue-specific, and six housekeeping. Comparison to qRT-PCR data for four genes (PCA3, OR51E2, FOLH1, and RPLP2) was strong (r = 0.51-0.95, Spearman p < 0.00001). Comparing EV to Cells, differential gene expression analysis found 57 gene-probes significantly more highly expressed in 100 ng of amplified cDNA products from the EV fraction, and 26 in Cells (p < 0.05; edgeR). Expression levels of prostate-specific genes (KLK2, KLK3) measured were ~20× higher in EVs, while PTPRC (white-blood Cells) was ~1000× higher in Cells. Boruta analysis identified 11 gene-probes as useful in detecting PCa: two were useful in both fractions (PCA3, HOXC6), five in EVs alone (GJB1, RPS10, TMPRSS2:ERG, ERG_Exons_4-5, HPN) and four from Cell (ERG_Exons_6-7, OR51E2, SPINK1, IMPDH2), suggesting that it is beneficial to fractionate whole urine prior to analysis. The five housekeeping genes were not significantly differentially expressed between PCa and non-cancer samples. Expression signatures from Cell, EV and combined data did not show evidence for one fraction providing superior information over the other.

Whole exome sequencing study identifies candidate loss of function variants and locus heterogeneity in familial cholesteatoma.

Cholesteatoma is a rare progressive disease of the middle ear. Most cases are sporadic, but some patients report a positive family history. Identifying functionally important gene variants associated with this disease has the potential to uncover the molecular basis of cholesteatoma pathology with implications for disease prevention, surveillance, or management. We performed an observational WES study of 21 individuals treated for cholesteatoma who were recruited from ten multiply affected families. These family studies were complemented with gene-level mutational burden analysis. We also applied functional enrichment analyses to identify shared properties and pathways for candidate genes and their products. Filtered data collected from pairs and trios of participants within the ten families revealed 398 rare, loss of function (LOF) variants co-segregating with cholesteatoma in 389 genes. We identified six genes DENND2C, DNAH7, NBEAL1, NEB, PRRC2C, and SHC2, for which we found LOF variants in two or more families. The parallel gene-level analysis of mutation burden identified a significant mutation burden for the genes in the DNAH gene family, which encode products involved in ciliary structure. Functional enrichment analyses identified common pathways for the candidate genes which included GTPase regulator activity, calcium ion binding, and degradation of the extracellular matrix. The number of candidate genes identified and the locus heterogeneity that we describe within and between multiply affected families suggest that the genetic architecture for familial cholesteatoma is complex.

Major data analysis errors invalidate cancer microbiome findings.

IMPORTANCE: Recent reports showing that human cancers have a distinctive microbiome have led to a flurry of papers describing microbial signatures of different cancer types. Many of these reports are based on flawed data that, upon re-analysis, completely overturns the original findings. The re-analysis conducted here shows that most of the microbes originally reported as associated with cancer were not present at all in the samples. The original report of a cancer microbiome and more than a dozen follow-up studies are, therefore, likely to be invalid.

Major data analysis errors invalidate cancer microbiome findings.

We re-analyzed the data from a recent large-scale study that reported strong correlations between microbial organisms and 33 different cancer types, and that created machine learning predictors with near-perfect accuracy at distinguishing among cancers. We found at least two fundamental flaws in the reported data and in the methods: (1) errors in the genome database and the associated computational methods led to millions of false positive findings of bacterial reads across all samples, largely because most of the sequences identified as bacteria were instead human; and (2) errors in transformation of the raw data created an artificial signature, even for microbes with no reads detected, tagging each tumor type with a distinct signal that the machine learning programs then used to create an apparently accurate classifier. Each of these problems invalidates the results, leading to the conclusion that the microbiome-based classifiers for identifying cancer presented in the study are entirely wrong. These flaws have subsequently affected more than a dozen additional published studies that used the same data and whose results are likely invalid as well.

Combining Molecular Subtypes with Multivariable Clinical Models Has the Potential to Improve Prediction of Treatment Outcomes in Prostate Cancer at Diagnosis.

Clinical management of prostate cancer is challenging because of its highly variable natural history and so there is a need for improved predictors of outcome in non-metastatic men at the time of diagnosis. In this study we calculated the model score from the leading clinical multivariable model, PREDICT prostate, and the poor prognosis DESNT molecular subtype, in a combined expression and clinical dataset that were taken from malignant tissue at prostatectomy (n = 359). Both PREDICT score (p < 0.0001, IQR HR = 1.59) and DESNT score (p < 0.0001, IQR HR = 2.08) were significant predictors for time to biochemical recurrence. A joint model combining the continuous PREDICT and DESNT score (p < 0.0001, IQR HR = 1.53 and 1.79, respectively) produced a significantly improved predictor than either model alone (p < 0.001). An increased probability of mortality after diagnosis, as estimated by PREDICT, was characterised by upregulation of cell-cycle related pathways and the downregulation of metabolism and cholesterol biosynthesis. The DESNT molecular subtype has distinct biological characteristics to those associated with the PREDICT model. We conclude that the inclusion of biological information alongside current clinical prognostic tools has the potential to improve the ability to choose the optimal treatment pathway for a patient.

A Model to Detect Significant Prostate Cancer Integrating Urinary Peptide and Extracellular Vesicle RNA Data.

There is a clinical need to improve assessment of biopsy-naïve patients for the presence of clinically significant prostate cancer (PCa). In this study, we investigated whether the robust integration of expression data from urinary extracellular vesicle RNA (EV-RNA) with urine proteomic metabolites can accurately predict PCa biopsy outcome. Urine samples collected within the Movember GAP1 Urine Biomarker study (n = 192) were analysed by both mass spectrometry-based urine-proteomics and NanoString gene-expression analysis (167 gene-probes). Cross-validated LASSO penalised regression and Random Forests identified a combination of clinical and urinary biomarkers for predictive modelling of significant disease (Gleason Score (Gs) ≥ 3 + 4). Four predictive models were developed: 'MassSpec' (CE-MS proteomics), 'EV-RNA', and 'SoC' (standard of care) clinical data models, alongside a fully integrated omics-model, deemed 'ExoSpec'. ExoSpec (incorporating four gene transcripts, six peptides, and two clinical variables) is the best model for predicting Gs ≥ 3 + 4 at initial biopsy (AUC = 0.83, 95% CI: 0.77−0.88) and is superior to a standard of care (SoC) model utilising clinical data alone (AUC = 0.71, p < 0.001, 1000 resamples). As the ExoSpec Risk Score increases, the likelihood of higher-grade PCa on biopsy is significantly greater (OR = 2.8, 95% CI: 2.1−3.7). The decision curve analyses reveals that ExoSpec provides a net benefit over SoC and could reduce unnecessary biopsies by 30%.