Search

Abstract Breast cancer remains a highly prevalent disease with considerable inter- and intra-tumoral heterogeneity complicating prognostication and treatment decisions. The utilization and depth of genomic, transcriptomic and proteomic data for cancer has exploded over recent times and the addition of spatial context to this information, by understanding the correlating morphologic and spatial patterns of cells in tissue samples, has created an exciting frontier of research, histo-genomics. At the same time, deep learning (DL), a class of machine learning algorithms employing artificial neural networks, has rapidly progressed in the last decade with a confluence of technical developments - including the advent of modern graphic processing units (GPU), allowing efficient implementation of increasingly complex architectures at scale; advances in the theoretical and practical design of network architectures; and access to larger datasets for training - all leading to sweeping advances in image classification and object detection. In this review, we examine recent developments in the application of DL in breast cancer histology with particular emphasis of those producing biologic insights or novel biomarkers, spanning the extraction of genomic information to the use of stroma to predict cancer recurrence, with the aim of suggesting avenues for further advancing this exciting field.

Despite unprecedented efforts, our therapeutic arsenal against SARS-CoV-2 remains limited. The conserved macrodomain 1 (Mac1) in NSP3 is an enzyme exhibiting ADP-ribosylhydrolase activity and a possible drug target. To determine the role of Mac1 catalytic activity in viral replication, we generated recombinant viruses and replicons encoding a catalytically inactive NSP3 Mac1 domain by mutating a critical asparagine in the active site. While substitution to alanine (N40A) reduced catalytic activity by ~10-fold, mutations to aspartic acid (N40D) reduced activity by ~100-fold relative to wild-type. Importantly, the N40A mutation rendered Mac1 unstable in vitro and lowered expression levels in bacterial and mammalian cells. When incorporated into SARS-CoV-2 molecular clones, the N40D mutant only modestly affected viral fitness in immortalized cell lines, but reduced viral replication in human airway organoids by 10-fold. In mice, the N40D mutant replicated at >1000-fold lower levels compared to the wild-type virus while inducing a robust interferon response; all animals infected with the mutant virus survived infection. Our data validate the critical role of SARS-CoV-2 NSP3 Mac1 catalytic activity in viral replication and as a promising therapeutic target to develop antivirals.

It is hypothesized that there are a number of tumor specific driver genes for metastatic prostate cancer. Here, the authors perform genome-wide CRISPRi screens and integrate these data with metastatic prostate cancer functional and clinical genomics data to show that KIF4A and WDR62 drive aggressive prostate cancer phenotypes.

Circulating DNA detected in plasma can be used for diagnostic purposes. Here, the authors show that the 5-hydroxymethyl cytosine biomarker from plasma-derived cell free DNA can be used to detect early stage pancreatic cancer.

Tracking tumour evolution in a patient via circulating tumour DNA (ctDNA) is complicated due to the unknown mix of fragmented alleles from different cancer lesions. Here, the authors make use of a rapid autopsy program to demonstrate how representative ctDNA profiling is of metastasis, as well as presenting methylation profiling method to track evolutionary change.

In breast cancer the contribution of different genetic variants to disease heritability is complex and not fully understood. Here, the authors present a network-based analysis in 84,567 patients studying ~7.3 million variants, identifying gene modules associated with breast cancer survival.

Cohesin complex mediates cohesion of sister chromatids and DNA loop formation. Here the authors show another role of cohesin in replication fork progression, suggesting a potential therapeutic strategy for cohesin-mutant cancers.

The mechanisms of PARP inhibitor (PARPi) resistance are poorly understood. Here the authors employ a CRISPR mutagenesis approach to identify PARP1 mutants causing PARPi resistance and find that PARP1 mutations are tolerated in BRCA1 mutated cells, suggesting alternative resistance mechanisms.

Poly ADP-ribosylation (PARylation) is a highly dynamic post-translation protein modification, but most methods only detect stable PARylation events. Here the authors develop a split-GFP-based sensor for PARylation detection in live cells and use it to identify a new centrosomal PARylation target.

Summary: We have performed an in-depth single-cell phenotypic characterization of high-grade serous ovarian cancer (HGSOC) by multiparametric mass cytometry (CyTOF). Using a CyTOF antibody panel to interrogate features of HGSOC biology, combined with unsupervised computational analysis, we identified noteworthy cell types co-occurring across the tumors. In addition to a dominant cell subset, each tumor harbored rarer cell phenotypes. One such group co-expressed E-cadherin and vimentin (EV), suggesting their potential role in epithelial mesenchymal transition, which was substantiated by pairwise correlation analyses. Furthermore, tumors from patients with poorer outcome had an increased frequency of another rare cell type that co-expressed vimentin, HE4, and cMyc. These poorer-outcome tumors also populated more cell phenotypes, as quantified by Simpson’s diversity index. Thus, despite the recognized genomic complexity of the disease, the specific cell phenotypes uncovered here offer a focus for therapeutic intervention and disease monitoring.

Abstract Background Plausible biological reasons exist regarding why smoking could affect breast cancer risk, but epidemiological evidence is inconsistent. Methods We used serial questionnaire information from the Generations Study cohort (United Kingdom) to estimate HRs for breast cancer in relation to smoking adjusted for potentially confounding factors, including alcohol intake. Results Among 102,927 women recruited 2003–2013, with an average of 7.7 years of follow-up, 1815 developed invasive breast cancer. The HR (reference group was never smokers) was 1.14 (95% CI 1.03–1.25; P = 0.010) for ever smokers, 1.24 (95% CI 1.08–1.43; P = 0.002) for starting smoking at ages

As a consequence of responding to colleagues who asked about the publication of the original article [1], the authors have determined that the data published in Table 4 of the paper are incorrect.

Mutations in the BAF SWI/SNF complex subunits are frequent in cancers but selective therapeutic approaches are not available yet. Here, the authors demonstrate that defects ofARID1Aand other subunits sensitizes cancer cells to the DNA checkpoint kinase inhibitor ATR in a synthetic lethal manner.

[This corrects the article DOI: 10.1371/journal.pone.0004264.].

One approach to identifying cancer-specific vulnerabilities and therapeutic targets is to profile genetic dependencies in cancer cell lines. Here, we describe data from a series of siRNA screens that identify the kinase genetic dependencies in 117 cancer cell lines from ten cancer types. By integrating the siRNA screen data with molecular profiling data, including exome sequencing data, we show how vulnerabilities/genetic dependencies that are associated with mutations in specific cancer driver genes can be identified. By integrating additional data sets into this analysis, including protein-protein interaction data, we also demonstrate that the genetic dependencies associated with many cancer driver genes form dense connections on functional interaction networks. We demonstrate the utility of this resource by using it to predict the drug sensitivity of genetically or histologically defined subsets of tumor cell lines, including an increased sensitivity of osteosarcoma cell lines to FGFR inhibitors and SMAD4 mutant tumor cells to mitotic inhibitors.

[This corrects the article DOI: 10.1371/journal.pone.0176578.].

[This corrects the article DOI: 10.1371/journal.pone.0149099.].

A previously fit and well man developed acute respiratory failure due to environmental mould exposure from living in damp rental accommodation. Despite aggressive intensive care management he rapidly deteriorated and required respiratory and cardiac Extracorporeal Membrane Oxygenation. We hypothesize that poor domiciliary conditions may make an underestimated contribution to community respiratory disease. These conditions may present as acute and severe illness with non-typical pathogens identified.

Activating KRAS mutations are found in approximately 20% of human cancers but no RAS-directed therapies are currently available. Here we describe a novel, robust, KRAS synthetic lethal interaction with the cyclin dependent kinase, CDK1. This was discovered using parallel siRNA screens in KRAS mutant and wild type colorectal isogenic tumour cells and subsequently validated in a genetically diverse panel of 26 colorectal and pancreatic tumour cell models. This established that the KRAS/CDK1 synthetic lethality applies in tumour cells with either amino acid position 12 (p.G12V, pG12D, p.G12S) or amino acid position 13 (p.G13D) KRAS mutations and can also be replicated in vivo in a xenograft model using a small molecule CDK1 inhibitor. Mechanistically, CDK1 inhibition caused a reduction in the S-phase fraction of KRAS mutant cells, an effect also characterised by modulation of Rb, a master control of the G1/S checkpoint. Taken together, these observations suggest that the KRAS/CDK1 interaction is a robust synthetic lethal effect worthy of further investigation.

Abstract PARP inhibitors have been proposed as a potential targeted therapy for patients with triple‐negative (ER‐, PR‐, HER2‐negative) breast cancers. However, it is as yet unclear as to whether single agent or combination therapy using PARP inhibitors would be most beneficial. To better understand the mechanisms that determine the response to PARP inhibitors, we investigated whether enzymes involved in metabolism of the PARP substrate, β‐NAD+, might alter the response to a clinical PARP inhibitor. Using an olaparib sensitization screen in a triple‐negative (TN) breast cancer model, we identified nicotinamide phosphoribosyltransferase (NAMPT) as a non‐redundant modifier of olaparib response. NAMPT is a rate‐limiting enzyme involved in the generation of the PARP substrate β‐NAD+ and the suppression of β‐NAD+ levels by NAMPT inhibition most likely explains these observations. Importantly, the combination of a NAMPT small molecule inhibitor, FK866, with olaparib inhibited TN breast tumour growth in vivo to a greater extent than either single agent alone suggesting that assessing NAMPT/PARP inhibitor combinations for the treatment of TN breast cancer may be warranted.

Invasive lobular breast cancer (ILC) accounts for 10-15% of all invasive breast carcinomas. It is generally ER positive (ER+) and often associated with lobular carcinoma in situ (LCIS). Genome-wide association studies have identified more than 70 common polymorphisms that predispose to breast cancer, but these studies included predominantly ductal (IDC) carcinomas. To identify novel common polymorphisms that predispose to ILC and LCIS, we pooled data from 6,023 cases (5,622 ILC, 401 pure LCIS) and 34,271 controls from 36 studies genotyped using the iCOGS chip. Six novel SNPs most strongly associated with ILC/LCIS in the pooled analysis were genotyped in a further 516 lobular cases (482 ILC, 36 LCIS) and 1,467 controls. These analyses identified a lobular-specific SNP at 7q34 (rs11977670, OR (95%CI) for ILC = 1.13 (1.09-1.18), P = 6.0 × 10(-10); P-het for ILC vs IDC ER+ tumors = 1.8 × 10(-4)). Of the 75 known breast cancer polymorphisms that were genotyped, 56 were associated with ILC and 15 with LCIS at P

Genetic variations, such as single nucleotide polymorphisms (SNPs) in microRNAs (miRNA) or in the miRNA binding sites may affect the miRNA dependent gene expression regulation, which has been implicated in various cancers, including breast cancer, and may alter individual susceptibility to cancer. We investigated associations between miRNA related SNPs and breast cancer risk. First we evaluated 2,196 SNPs in a case-control study combining nine genome wide association studies (GWAS). Second, we further investigated 42 SNPs with suggestive evidence for association using 41,785 cases and 41,880 controls from 41 studies included in the Breast Cancer Association Consortium (BCAC). Combining the GWAS and BCAC data within a meta-analysis, we estimated main effects on breast cancer risk as well as risks for estrogen receptor (ER) and age defined subgroups. Five miRNA binding site SNPs associated significantly with breast cancer risk: rs1045494 (odds ratio (OR) 0.92; 95% confidence interval (CI): 0.88-0.96), rs1052532 (OR 0.97; 95% CI: 0.95-0.99), rs10719 (OR 0.97; 95% CI: 0.94-0.99), rs4687554 (OR 0.97; 95% CI: 0.95-0.99, and rs3134615 (OR 1.03; 95% CI: 1.01-1.05) located in the 3' UTR of CASP8, HDDC3, DROSHA, MUSTN1, and MYCL1, respectively. DROSHA belongs to miRNA machinery genes and has a central role in initial miRNA processing. The remaining genes are involved in different molecular functions, including apoptosis and gene expression regulation. Further studies are warranted to elucidate whether the miRNA binding site SNPs are the causative variants for the observed risk effects.

Genetic perturbation screens have the potential to dissect a wide range of cellular phenotypes. Such screens have historically been difficult in diploid mammalian cells. The recent derivation of haploid embryonic stem cells provides an opportunity to cause loss of function mutants with a random mutagen in a mammalian cell with a normal genetic background. We describe an approach to genetic screens that exploits the highly active piggyBac transposon in haploid mammalian cells. As an example of haploid transposon (HTP) screening, we apply this approach to identifying determinants of cancer drug toxicity and resistance. In a screen for 6-thioguanine resistance we recovered components of the DNA mismatch repair pathway, a known requirement for toxicity. In a further screen for resistance to the clinical poly(ADP-ribose) polymerase (PARP) inhibitor olaparib we recovered multiple Parp1 mutants. Our results show that olaparib toxicity to normal cells is mediated predominantly via Parp1, and suggest that the clinical side effects of olaparib may be on target. The transposon mutant libraries are stable and can be readily reused to screen other drugs. The screening protocol described has several advantages over other methods such as RNA interference: it is rapid and low cost, and mutations can be easily reverted to establish causality.

The 6q25.1 locus was first identified via a genome-wide association study (GWAS) in Chinese women and marked by single nucleotide polymorphism (SNP) rs2046210, approximately 180 Kb upstream of ESR1. There have been conflicting reports about the association of this locus with breast cancer in Europeans, and a GWAS in Europeans identified a different SNP, tagged here by rs12662670. We examined the associations of both SNPs in up to 61,689 cases and 58,822 controls from forty-four studies collaborating in the Breast Cancer Association Consortium, of which four studies were of Asian and 39 of European descent. Logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI). Case-only analyses were used to compare SNP effects in Estrogen Receptor positive (ER+) versus negative (ER-) tumours. Models including both SNPs were fitted to investigate whether the SNP effects were independent. Both SNPs are significantly associated with breast cancer risk in both ethnic groups. Per-allele ORs are higher in Asian than in European studies [rs2046210: OR (A/G) = 1.36 (95% CI 1.26-1.48), p = 7.6 × 10(-14) in Asians and 1.09 (95% CI 1.07-1.11), p = 6.8 × 10(-18) in Europeans. rs12662670: OR (G/T) = 1.29 (95% CI 1.19-1.41), p = 1.2 × 10(-9) in Asians and 1.12 (95% CI 1.08-1.17), p = 3.8 × 10(-9) in Europeans]. SNP rs2046210 is associated with a significantly greater risk of ER- than ER+ tumours in Europeans [OR (ER-) = 1.20 (95% CI 1.15-1.25), p = 1.8 × 10(-17) versus OR (ER+) = 1.07 (95% CI 1.04-1.1), p = 1.3 × 10(-7), p(heterogeneity) = 5.1 × 10(-6)]. In these Asian studies, by contrast, there is no clear evidence of a differential association by tumour receptor status. Each SNP is associated with risk after adjustment for the other SNP. These results suggest the presence of two variants at 6q25.1 each independently associated with breast cancer risk in Asians and in Europeans. Of these two, the one tagged by rs2046210 is associated with a greater risk of ER- tumours.

PTEN (Phosphatase and tensin homolog) is a tumour suppressor gene commonly defective in human cancer, and is thus a potentially important therapeutic target. Targeting tumour suppressor loss-of-function is possible by exploiting the genetic concept of synthetic lethality (SL). By combining the use of isogenic models of PTEN deficiency with high-throughput RNA interference (RNAi) screening, we have identified Nemo-Like Kinase (NLK) inhibition as being synthetically lethal with PTEN deficiency. This SL is likely mediated by the transcription factor FOXO1 (Forkhead box O1), an NLK substrate, as the selectivity of NLK gene silencing for PTEN deficient cells can be reversed by FOXO1 knockdown. In addition, we provide evidence that PTEN defective cells targeted by NLK gene depletion undergo senescence, suggesting that NLK function is critical for the continued proliferation of PTEN deficient cells. Taken together, these data provide new insight into the potential of targeting of NLK to treat a range of tumourigenic conditions characterised by PTEN deficiency.

Next generation DNA sequencing (NGS) technologies have revolutionized the pace at which whole genome and exome sequences can be generated. However, despite these advances, many of the methods for targeted resequencing, such as the generation of high-depth exome sequences, are somewhat limited by the relatively large amounts of starting DNA that are normally required. In the case of tumour analysis this is particularly pertinent as many tumour biopsies often return submicrogram quantities of DNA, especially when tumours are microdissected prior to analysis. Here, we present a method for exome capture and resequencing using as little as 50 ng of starting DNA. The sequencing libraries generated by this minimal starting amount (MSA-Cap) method generate datasets that are comparable to standard amount (SA) whole exome libraries that use three micrograms of starting DNA. This method, which can be performed in most laboratories using commonly available reagents, has the potential to enhance large scale profiling efforts such as the resequencing of tumour exomes.

Male breast cancer accounts for approximately 1% of all breast cancer. To date, risk factors for male breast cancer are poorly defined, but certain risk factors and genetic features appear common to both male and female breast cancer. Genome-wide association studies (GWAS) have recently identified common single nucleotide polymorphisms (SNPs) that influence female breast cancer risk; 12 of these have been independently replicated. To examine if these variants contribute to male breast cancer risk, we genotyped 433 male breast cancer cases and 1,569 controls. Five SNPs showed a statistically significant association with male breast cancer: rs13387042 (2q35) (odds ratio (OR) = 1.30, p = 7.98×10⁻⁴), rs10941679 (5p12) (OR = 1.26, p = 0.007), rs9383938 (6q25.1) (OR = 1.39, p = 0.004), rs2981579 (FGFR2) (OR = 1.18, p = 0.03), and rs3803662 (TOX3) (OR = 1.48, p = 4.04×10⁻⁶). Comparing the ORs for male breast cancer with the published ORs for female breast cancer, three SNPs--rs13387042 (2q35), rs3803662 (TOX3), and rs6504950 (COX11)--showed significant differences in ORs (p

Capan-1 is a well-characterised BRCA2-deficient human cell line isolated from a liver metastasis of a pancreatic adenocarcinoma. Here we report a genome-wide assessment of structural variations and high-depth exome characterization of single nucleotide variants and small insertion/deletions in Capan-1. To identify potential somatic and tumour-associated variations in the absence of a matched-normal cell line, we devised a novel method based on the analysis of HapMap samples. We demonstrate that Capan-1 has one of the most rearranged genomes sequenced to date. Furthermore, small insertions and deletions are detected more frequently in the context of short sequence repeats than in other genomes. We also identify a number of novel mutations that may represent genetic changes that have contributed to tumour progression. These data provide insight into the genomic effects of loss of BRCA2 function.

Epidemiological studies have shown that one of the strongest risk factors for prostate cancer is a family history of the disease, suggesting that inherited factors play a major role in prostate cancer susceptibility. Germline mutations in BRCA2 predispose to breast and ovarian cancer with its predominant tumour suppressor function thought to be the repair of DNA double-strand breaks. BRCA2 has also been implicated in prostate cancer etiology, but it is unclear the impact that mutations in this gene have on prostate tumourigenesis. Here we have undertaken a genetic analysis in the mouse to determine the role of Brca2 in the adult prostate. We show that deletion of Brca2 specifically in prostate epithelia results in focal hyperplasia and low-grade prostate intraepithelial neoplasia (PIN) in animals over 12 months of age. Simultaneous deletion of Brca2 and the tumour suppressor Trp53 in prostate epithelia gave rise to focal hyperplasia and atypical cells at 6 months, leading to high-grade PIN in animals from 12 months. Epithelial cells in these lesions show an increase in DNA damage and have higher levels of proliferation, but also elevated apoptosis. Castration of Brca2;Trp53 mutant animals led to regression of PIN lesions, but atypical cells persisted that continued to proliferate and express nuclear androgen receptor. This study provides evidence that Brca2 can act as a tumour suppressor in the prostate, and the model we describe should prove useful in the development of new therapeutic approaches.

The majority of new drug approvals for cancer are based on existing therapeutic targets. One approach to the identification of novel targets is to perform high-throughput RNA interference (RNAi) cellular viability screens. We describe a novel approach combining RNAi screening in multiple cell lines with gene expression and genomic profiling to identify novel cancer targets. We performed parallel RNAi screens in multiple cancer cell lines to identify genes that are essential for viability in some cell lines but not others, suggesting that these genes constitute key drivers of cellular survival in specific cancer cells. This approach was verified by the identification of PIK3CA, silencing of which was selectively lethal to the MCF7 cell line, which harbours an activating oncogenic PIK3CA mutation. We combined our functional RNAi approach with gene expression and genomic analysis, allowing the identification of several novel kinases, including WEE1, that are essential for viability only in cell lines that have an elevated level of expression of this kinase. Furthermore, we identified a subset of breast tumours that highly express WEE1 suggesting that WEE1 could be a novel therapeutic target in breast cancer. In conclusion, this strategy represents a novel and effective strategy for the identification of functionally important therapeutic targets in cancer.

The Lkb1 tumour suppressor is a multitasking kinase participating in a range of physiological processes. We have determined the impact of Lkb1 deficiency on intestinal homeostasis, particularly focussing on secretory cell differentiation and development since we observe strong expression of Lkb1 in normal small intestine Paneth and goblet cells. We crossed mice bearing an Lkb1 allele flanked with LoxP sites with those carrying a Cyp1a1-specific inducible Cre recombinase. Lkb1 was efficiently deleted from the epithelial cells of the mouse intestine after intraperitoneal injection of the inducing agent beta-naphthoflavone. Bi-allelic loss of Lkb1 led to the perturbed development of Paneth and goblet cell lineages. These changes were characterised by the lack of Delta ligand expression in Lkb1-deficient secretory cells and a significant increase in the levels of the downstream Notch signalling effector Hes5 but not Hes1. Our data show that Lkb1 is required for the normal differentiation of secretory cell lineages within the intestine, and that Lkb1 deficiency modulates Notch signalling modulation in post-mitotic cells.

Specification of mammary epithelial cell fate occurs during embryogenesis as cells aggregate to form the mammary anlage. Within the embryonic mammary bud, a population of epithelial cells exists that will subsequently proliferate to form a ductal tree filling the stromal compartment, and which can produce milk upon terminal differentiation after birth. Subsequently, these structures can be remodelled and returned to a basal state after weaning before regenerating in future pregnancies. The plasticity of the mammary epithelial cell, and its responsiveness to hormone receptors, facilitates this amazing biological feat, but aberrant signalling may also result in unintended consequences in the form of frequent malignancies. Reflecting this intimate connection, a considerable number of signalling pathways have been implicated in both mammary gland morphogenesis and carcinogenesis.