Genome Medicine

Patient-derived xenograft models capture genomic heterogeneity in endometrial cancer

Abstract Background Endometrial cancer (EC) is a major gynecological cancer with increasing incidence. It comprises four molecular subtypes with differing etiology, prognoses, and responses to chemotherapy. In the future, clinical trials testing new single agents or combination therapies will be targeted to the molecular subtype most likely to respond. As pre-clinical models that faithfully represent the molecular subtypes of EC are urgently needed, we sought to develop and characterize a panel of novel EC patient-derived xenograft (PDX) models. Methods Here, we report whole exome or whole genome sequencing of 11 PDX models and their matched primary tumor. Analysis of multiple PDX lineages and passages was performed to study tumor heterogeneity across lineages and/or passages. Based on recent reports of frequent defects in the homologous recombination (HR) pathway in EC, we assessed mutational signatures and HR deficiency scores and correlated these with in vivo responses to the PARP inhibitor (PARPi) talazoparib in six PDXs representing the copy number high/p53-mutant and mismatch-repair deficient molecular subtypes of EC. Results PDX models were successfully generated from grade 2/3 tumors, including three uterine carcinosarcomas. The models showed similar histomorphology to the primary tumors and represented all four molecular subtypes of EC, including five mismatch-repair deficient models. The different PDX lineages showed a wide range of inter-tumor and intra-tumor heterogeneity. However, for most PDX models, one arm recapitulated the molecular landscape of the primary tumor without major genomic drift. An in vivo response to talazoparib was detected in four copy number high models. Two models (carcinosarcomas) showed a response consistent with stable disease and two models (one copy number high serous EC and another carcinosarcoma) showed significant tumor growth inhibition, albeit one consistent with progressive disease; however, all lacked the HR deficiency genomic signature. Conclusions EC PDX models represent the four molecular subtypes of disease and can capture intra-tumor heterogeneity of the original primary tumor. PDXs of the copy number high molecular subtype showed sensitivity to PARPi; however, deeper and more durable responses will likely require combination of PARPi with other agents.

Developing an optimal stratification model for colorectal cancer screening and reducing racial disparities in multi-center population-based studies

Abstract Background Early detection of colorectal neoplasms can reduce the colorectal cancer (CRC) burden by timely intervention for high-risk individuals. However, effective risk prediction models are lacking for personalized CRC early screening in East Asian (EAS) population. We aimed to develop, validate, and optimize a comprehensive risk prediction model across all stages of the dynamic adenoma-carcinoma sequence in EAS population. Methods To develop precision risk-stratification and intervention strategies, we developed three trans-ancestry PRSs targeting colorectal neoplasms: (1) using 148 previously identified CRC risk loci (PRS148); (2) SNPs selection from large-scale meta-analysis data by clumping and thresholding (PRS183); (3) PRS-CSx, a Bayesian approach for genome-wide risk prediction (PRSGenomewide). Then, the performance of each PRS was assessed and validated in two independent cross-sectional screening sets, including 4600 patients with advanced colorectal neoplasm, 4495 patients with non-advanced adenoma, and 21,199 normal individuals from the ZJCRC (Zhejiang colorectal cancer set; EAS) and PLCO (the Prostate, Lung, Colorectal, and Ovarian Cancer Screening Trial; European, EUR) studies. The optimal PRS was further incorporated with lifestyle factors to stratify individual risk and ultimately tested in the PLCO and UK Biobank prospective cohorts, totaling 350,013 participants. Results Three trans-ancestry PRSs achieved moderately improved predictive performance in EAS compared to EUR populations. Remarkably, the PRSs effectively facilitated a thorough risk assessment across all stages of the dynamic adenoma-carcinoma sequence. Among these models, PRS183 demonstrated the optimal discriminatory ability in both EAS and EUR validation datasets, particularly for individuals at risk of colorectal neoplasms. Using two large-scale and independent prospective cohorts, we further confirmed a significant dose–response effect of PRS183 on incident colorectal neoplasms. Incorporating PRS183 with lifestyle factors into a comprehensive strategy improves risk stratification and discriminatory accuracy compared to using PRS or lifestyle factors separately. This comprehensive risk-stratified model shows potential in addressing missed diagnoses in screening tests (best NPV = 0.93), while moderately reducing unnecessary screening (best PPV = 0.32). Conclusions Our comprehensive risk-stratified model in population-based CRC screening trials represents a promising advancement in personalized risk assessment, facilitating tailored CRC screening in the EAS population. This approach enhances the transferability of PRSs across ancestries and thereby helps address health disparity. Graphical Abstract

Large-scale copy number alterations are enriched for synthetic viability in BRCA1/BRCA2 tumors

AbstractBackgroundPathogenicBRCA1orBRCA2germline mutations contribute to hereditary breast, ovarian, prostate, and pancreatic cancer. Paradoxically, bi-allelic inactivation ofBRCA1orBRCA2(bBRCA1/2) is embryonically lethal and decreases cellular proliferation.The compensatory mechanisms that facilitate oncogenesis in bBRCA1/2 tumors remain unclear.MethodsWe identified recurrent genetic alterations enriched in human bBRCA1/2 tumors and experimentally validated if these improved proliferation in cellular models. We analyzed mutations and copy number alterations (CNAs) in bBRCA1/2 breast and ovarian cancer from the TCGA and ICGC. We used Fisher’s exact test to identify CNAs enriched in bBRCA1/2 tumors compared to control tumors that lacked evidence of homologous recombination deficiency. Genes located in CNA regions enriched in bBRCA1/2 tumors were further screened by gene expression and their effects on proliferation in genome-wide CRISPR/Cas9 screens. A set of candidate genes was functionally validated with in vitro clonogenic survival and functional assays to validate their influence on proliferation in the setting of bBRCA1/2 mutations.ResultsWe found that bBRCA1/2 tumors harbor recurrent large-scale genomic deletions significantly more frequently than histologically matched controls (n = 238 cytobands in breast and ovarian cancers). Within the deleted regions, we identified 277 BRCA1-related genes and 218 BRCA2-related genes that had reduced expression and increased proliferation in bBRCA1/2 but not in wild-type cells in genome-wide CRISPR screens. In vitro validation of 20 candidate genes with clonogenic proliferation assays validated 9 genes, includingRIC8AandATMIN(ATM-Interacting protein). We identified loss ofRIC8A, which occurs frequently in both bBRCA1/2 tumors and is synthetically viable with loss of bothBRCA1andBRCA2. Furthermore, we found that metastatic homologous recombination deficient cancers acquire loss-of-function mutations inRIC8A. Lastly, we identified thatRIC8Adoes not rescue homologous recombination deficiency but may influence mitosis in bBRCA1/2 tumors, potentially leading to increased micronuclei formation.ConclusionsThis study provides a means to solve the tumor suppressor paradox by identifying synthetic viability interactions and causal driver genes affected by large-scale CNAs in human cancers.

Efficacy and safety of novel multiple-chain DAP-CAR-T cells targeting mesothelin in ovarian cancer and mesothelioma: a single-arm, open-label and first-in-human study

Despite remarkable achievements in applying chimeric antigen receptor (CAR)-T cells to treat hematological malignancies, they remain much less effective against solid tumors, facing several challenges affecting their clinical use. We previously showed that multichain DNAX-activating protein (DAP) CAR structures could enhance the safety and efficacy of CAR-T cells when used against solid tumors. In particular, mesothelin (MSLN)-targeted CAR-T cell therapy has therapeutic potential in MSLN-positive solid tumors, including ovarian cancer and mesothelioma. In vitro cell killing assays and xenograft model were utilized to determine the anti-tumor efficacy of MSLN targeting DAP-CAR-T cells and other CAR-T cells. ELISA and flow cytometry analysis were used to assess the cytokine secretion capacity and proliferation ability. Eight patients with MSLN expression were enrolled to evaluate the safety and efficacy of MSLN-DAP CAR-T cell therapy. Single-cell sequencing was performed to explore the dynamics of immune cells in patients during treatment and to identify the transcriptomic signatures associated with efficacy and toxicity. We found that multichain DAP-CAR formed by combining a natural killer cell immunoglobulin-like receptor truncator and DAP12 exhibited better cytotoxicity and tumor-killing capacity than other natural killer cell-activated receptors associated with DAP12, DAP10, or CD3Z. The safety and efficacy of MSLN-DAP CAR-T cell therapy in patients with ovarian cancer and mesothelioma were evaluated in a single-arm, open-label clinical trial (ChiCTR2100046544); two patients achieved partial response, while four patients had a stable disease status. Furthermore, single-cell sequencing analysis indicated that KT032 CAR-T cell infusion could recruit more immune cells and temporarily remodel the TME. Our study highlights the safety and therapeutic efficacy of multiple-chain DAP-CAR-T cell therapy targeting MSLN to treat patients with ovarian cancer and mesothelioma. ChiCTR.org.cn, ChiCTR2100046544 . May 21, 2021.

The PARP1 selective inhibitor saruparib (AZD5305) elicits potent and durable antitumor activity in patient-derived BRCA1/2-associated cancer models

Distinct transcriptional programs stratify ovarian cancer cell lines into the five major histological subtypes

Abstract Background Epithelial ovarian cancer (OC) is a heterogenous disease consisting of five major histologically distinct subtypes: high-grade serous (HGSOC), low-grade serous (LGSOC), endometrioid (ENOC), clear cell (CCOC) and mucinous (MOC). Although HGSOC is the most prevalent subtype, representing 70–80% of cases, a 2013 landmark study by Domcke et al. found that the most frequently used OC cell lines are not molecularly representative of this subtype. This raises the question, if not HGSOC, from which subtype do these cell lines derive? Indeed, non-HGSOC subtypes often respond poorly to chemotherapy; therefore, representative models are imperative for developing new targeted therapeutics. Methods Non-negative matrix factorisation (NMF) was applied to transcriptomic data from 44 OC cell lines in the Cancer Cell Line Encyclopedia, assessing the quality of clustering into 2–10 groups. Epithelial OC subtypes were assigned to cell lines optimally clustered into five transcriptionally distinct classes, confirmed by integration with subtype-specific mutations. A transcriptional subtype classifier was then developed by trialling three machine learning algorithms using subtype-specific metagenes defined by NMF. The ability of classifiers to predict subtype was tested using RNA sequencing of a living biobank of patient-derived OC models. Results Application of NMF optimally clustered the 44 cell lines into five transcriptionally distinct groups. Close inspection of orthogonal datasets revealed this five-cluster delineation corresponds to the five major OC subtypes. This NMF-based classification validates the Domcke et al. analysis, in identifying lines most representative of HGSOC, and additionally identifies models representing the four other subtypes. However, NMF of the cell lines into two clusters did not align with the dualistic model of OC and suggests this classification is an oversimplification. Subtype designation of patient-derived models by a random forest transcriptional classifier aligned with prior diagnosis in 76% of unambiguous cases. In cases where there was disagreement, this often indicated potential alternative diagnosis, supported by a review of histological, molecular and clinical features. Conclusions This robust classification informs the selection of the most appropriate models for all five histotypes. Following further refinement on larger training cohorts, the transcriptional classification may represent a useful tool to support the classification of new model systems of OC subtypes.

Genomic landscape of endometrial polyps

Endometrial polyps are common, localized overgrowths of endometrial glands and stroma that protrude into the uterine cavity. These tumor-like lesions can cause symptoms like abnormal uterine bleeding and infertility, and they may undergo malignant transformation. The etiology of endometrial polyps remains largely unknown. Here, we conducted whole-genome sequencing and global gene expression profiling on 23 polyps. Major findings were validated with targeted DNA (Sanger sequencing) and protein (immunohistochemistry) level analyses. Sanger sequencing was also utilized to validate the observed novel alterations in an additional set of 54 polyp samples. The most common alterations were chromosomal rearrangements affecting HMGA1 and HMGA2, identified in 74% (17/23) of the polyps. These rearrangements involved LRMDA, RAD51B, TRAF3IP2, and 7p15.2 as recurrent rearrangement partners. 3'RNA sequencing indicated corresponding overexpression of HMGA1 and HMGA2 as well as a downstream target PLAG1. Elevated protein level expression of HMGA1 and HMGA2 was further shown using immunohistochemistry. In addition to frequent HMGA1 and HMGA2 alterations, we found UBE2A as a novel candidate driver gene with highly specific recurrent mutations. We also identified recurrent low-allelic fraction mutations in well-established cancer genes KRAS, PIK3CA, PIK3R1, and PTEN. Here, we have characterized the genomic landscape of endometrial polyps. We show that chromosomal alterations affecting HMGA1 and HMGA2 are a major underlying cause for polyp development. In addition, we present UBE2A as a novel candidate gene for human tumorigenesis. Our results contribute to a better understanding of endometrial polyp development and pave the way towards the development of targeted, non-invasive treatment options.

A cell-of-origin epigenetic tracer reveals clinically distinct subtypes of high-grade serous ovarian cancer

Abstract Background High-grade serous ovarian cancer (HGSOC) is a major unmet need in oncology. The remaining uncertainty on its originating tissue has hampered the discovery of molecular oncogenic pathways and the development of effective therapies. Methods We used an approach based on the retention in tumors of a DNA methylation trace (OriPrint) that distinguishes the two putative tissues of origin of HGSOC, the fimbrial (FI) and ovarian surface epithelia (OSE), to stratify HGSOC by several clustering methods, both linear and non-linear. The identified tumor subtypes (FI-like and OSE-like HGSOC) were investigated at the RNAseq level to stratify an in-house cohort of macrodissected HGSOC FFPE samples to derive overall and disease-free survival and identify specific transcriptional alterations of the two tumor subtypes, both by classical differential expression and weighted correlation network analysis. We translated our strategy to published datasets and verified the co-occurrence of previously described molecular classification of HGSOC. We performed cytokine analysis coupled to immune phenotyping to verify alterations in the immune compartment associated with HGSOC. We identified genes that are both differentially expressed and methylated in the two tumor subtypes, concentrating on PAX8 as a bona fide marker of FI-like HGSOC. Results We show that: - OriPrint is a robust DNA methylation tracer that exposes the tissue of origin of HGSOC. - The tissue of origin of HGSOC is the main determinant of DNA methylation variance in HGSOC. - The tissue of origin is a prognostic factor for HGSOC patients. - FI-like and OSE-like HGSOC are endowed with specific transcriptional alterations that impact patients’ prognosis. - OSE-like tumors present a more invasive and immunomodulatory phenotype, compatible with its worse prognostic impact. - Among genes that are differentially expressed and regulated in FI-like and OSE-like HGSOC, PAX8 is a bona fide marker of FI-like tumors. Conclusions Through an integrated approach, our work demonstrates that both FI and OSE are possible origins for human HGSOC, whose derived subtypes are both molecularly and clinically distinct. These results will help define a new roadmap towards rational, subtype-specific therapeutic inroads and improved patients’ care.

A functionally impaired missense variant identified in French Canadian families implicates FANCI as a candidate ovarian cancer-predisposing gene

Abstract Background Familial ovarian cancer (OC) cases not harbouring pathogenic variants in either of the BRCA1 and BRCA2 OC-predisposing genes, which function in homologous recombination (HR) of DNA, could involve pathogenic variants in other DNA repair pathway genes. Methods Whole exome sequencing was used to identify rare variants in HR genes in a BRCA1 and BRCA2 pathogenic variant negative OC family of French Canadian (FC) ancestry, a population exhibiting genetic drift. OC cases and cancer-free individuals from FC and non-FC populations were investigated for carrier frequency of FANCI c.1813C>T; p.L605F, the top-ranking candidate. Gene and protein expression were investigated in cancer cell lines and tissue microarrays, respectively. Results In FC subjects, c.1813C>T was more common in familial (7.1%, 3/42) than sporadic (1.6%, 7/439) OC cases ( P = 0.048). Carriers were detected in 2.5% (74/2950) of cancer-free females though female/male carriers were more likely to have a first-degree relative with OC (121/5249, 2.3%; Spearman correlation = 0.037; P = 0.011), suggesting a role in risk. Many of the cancer-free females had host factors known to reduce risk to OC which could influence cancer risk in this population. There was an increased carrier frequency of FANCI c.1813C>T in BRCA1 and BRCA2 pathogenic variant negative OC families, when including the discovery family, compared to cancer-free females (3/23, 13%; OR = 5.8; 95%CI = 1.7–19; P = 0.005). In non-FC subjects, 10 candidate FANCI variants were identified in 4.1% (21/516) of Australian OC cases negative for pathogenic variants in BRCA1 and BRCA2 , including 10 carriers of FANCI c.1813C>T. Candidate variants were significantly more common in familial OC than in sporadic OC ( P = 0.04). Localization of FANCD2, part of the FANCI-FANCD2 (ID2) binding complex in the Fanconi anaemia (FA) pathway, to sites of induced DNA damage was severely impeded in cells expressing the p.L605F isoform. This isoform was expressed at a reduced level, destabilized by DNA damaging agent treatment in both HeLa and OC cell lines, and exhibited sensitivity to cisplatin but not to a poly (ADP-ribose) polymerase inhibitor. By tissue microarray analyses, FANCI protein was consistently expressed in fallopian tube epithelial cells and only expressed at low-to-moderate levels in 88% (83/94) of OC samples. Conclusions This is the first study to describe candidate OC variants in FANCI , a member of the ID2 complex of the FA DNA repair pathway. Our data suggest that pathogenic FANCI variants may modify OC risk in cancer families.

DNA repair and the contribution to chemotherapy resistance

Patterns of genomic instability in > 2000 patients with ovarian cancer across six clinical trials evaluating olaparib

Abstract Background The introduction of poly(ADP-ribose) polymerase (PARP) inhibitors represented a paradigm shift in the treatment of ovarian cancer. Genomic data from patients with high-grade ovarian cancer in six phase II/III trials involving the PARP inhibitor olaparib were analyzed to better understand patterns and potential causes of genomic instability. Patients and methods Homologous recombination deficiency (HRD) was assessed in 2147 tumor samples from SOLO1, PAOLA-1, Study 19, SOLO2, OPINION, and LIGHT using next-generation sequencing technology. Genomic instability scores (GIS) were assessed in BRCA1 and/or BRCA2 (BRCA)-mutated (BRCAm), non-BRCA homologous recombination repair-mutated (non-BRCA HRRm), and non-HRRm tumors. Results BRCAm was identified in 1021/2147 (47.6%) tumors. BRCAm tumors had significantly higher GIS than non-BRCAm tumors (P < 0.001) and high biallelic loss (815/838; 97.3%) regardless of germline (658/672; 97.9%) or somatic (101/108; 93.5%) BRCAm status. In non-BRCA HRRm tumors (n = 121) a similar proportion were HRD-positive (GIS ≥ 42: 55/121; 45.5%) relative to HRD-negative (GIS < 42: 52/121; 43.0%). GIS was highly variable in non-BRCA HRRm (median 42 [interquartile range (IQR) 29–58]) and non-HRRm (n = 1005; median 32 [IQR 20–55]) tumors. Gene mutations with high GIS included HRR genes BRIP1 (median 46 [IQR 41–58]), RAD51C (median 58 [IQR 48–66]), RAD51D (median 62 [IQR 54–69]), and PALB2 (median 64 [IQR 58–74]), and non-HRR genes NF1 (median 49 [IQR 25–60]) and RB1 (median 55 [IQR 30–71]). CCNE1-amplified and PIK3CA-mutated tumors had low GIS (CCNE1-amplified: median 24 [IQR 18–29]; PIK3CA-mutated: median 32 [IQR 14–52]) and were predominantly non-BRCAm. Conclusions These analyses provide valuable insight into patterns of genomic instability and potential drivers of HRD, besides BRCAm, in ovarian cancer and will help guide future research into the potential clinical effectiveness of anti-cancer treatments in ovarian cancer, including PARP inhibitors as well as other precision oncology agents. Trial registration The SOLO1 trial was registered at ClinicalTrials.gov (NCT01844986) on April 30, 2013; the PAOLA-1 trial was registered at ClinicalTrials.gov (NCT02477644) on June 18, 2015 (retrospectively registered); Study 19 was registered at ClinicalTrials.gov (NCT00753545) on September 12, 2008 (retrospectively registered); the SOLO2 trial was registered at ClinicalTrials.gov (NCT01874353) on June 7, 2013; the OPINION trial was registered at ClinicalTrials.gov (NCT03402841) on January 3, 2018; the LIGHT trial was registered at ClinicalTrials.gov (NCT02983799) on November 4, 2016.

Replication-associated mechanisms contribute to an increased CpG > TpG mutation burden in mismatch repair-deficient cancers

Abstract Background Single base substitution (SBS) mutations, particularly C > T and T > C, are increased owing to unrepaired DNA replication errors in mismatch repair-deficient (MMRd) cancers. Excess CpG > TpG mutations have been reported in MMRd cancers defective in mismatch detection (dMutSα), but not in mismatch correction (dMutLα). Somatic CpG > TpG mutations conventionally result from unrepaired spontaneous deamination of 5’-methylcytosine throughout the cell cycle, causing T:G mismatches and signature SBS1. It has been proposed that MutSα detects those mismatches, prior to error correction by base excision repair (BER). However, other evidence appears inconsistent with that hypothesis: for example, MutSα is specifically expressed in S/G2 phases of the cell cycle, and defects in replicative DNA polymerase proofreading specifically cause excess CpG > TpG mutations in signature SBS10b. Methods We analysed mutation spectra and COSMIC mutation signatures in whole-genome sequencing data from 1803 colorectal cancers (164 dMutLα, 20 dMutSα) and 596 endometrial cancers (103 dMutLα, 9 dMutSα) from the UK 100,000 Genomes Project. We mapped each C > T mutation to its genomic features, including normal DNA methylation state, replication timing, transcription strand, and replication strand, to investigate the mechanism(s) by which these mutations arise. Results We confirmed that dMutSα tumours specifically had higher CpG > TpG burdens than dMutLα tumours. We could fully reconstitute the observed dMutSα CpG > TpG mutation spectrum by adding CpG > TpG mutations in proportion to their SBS1 activity to the dMutLα spectrum. However, other evidence indicated that the SBS1 excess in dMutSα cancers did not come from 5’-methylcytosine deamination alone: non-CpG C > T mutations were also increased in dMutSα cancers; and, in contrast to tumours deficient in BER, CpG > TpG mutations were biased to the leading DNA replication strand, at similar levels in dMutSα and dMutLα cancers, suggesting an origin in DNA replication. Other substitution mutations usually corrected by BER were not increased in dMutSα tumours. Conclusions There is a CpG > TpG and SBS1 excess specific to dMutSα MMRd tumours, consistent with previous reports, and we find a general increase in somatic C > T mutations. Contrary to some other studies, the similar leading replication strand bias in both dMutSα and dMutLα tumours indicates that at least some of the excess CpG > TpG mutations arise via DNA replication errors, and not primarily via the replication-independent deamination of 5’-methylcytosine.

Joint inference of mutational signatures from indels and single-nucleotide substitutions reveals prognostic impact of DNA repair deficiencies

Mutational signatures are increasingly used to understand the mechanisms causing cancer. However, their important applications in predicting prognosis and stratifying patients for therapy are hampered by inaccurate inference of the various featureless, dense trinucleotide mutational spectra, which are often confounded with one another. One of them is the homologous recombination deficiency (HRd)-associated signature SBS3, relevant because of its association with prognosis in ovarian and breast cancer and because of its potential as a biomarker for synthetic lethality therapies. Here, we highlight strong benefits of a multimodal approach for mutational signature extraction, applied on top of standard bioinformatic pipelines. By jointly operating on single-base substitution (SBS) and indel (ID) spectra, this method enables accurate identification of various DNA repair deficiency signatures and patient survival prediction. Across four different cohorts of whole-genome sequenced high-grade serous ovarian cancers (HGSOC), the multimodal SBS + ID approach correctly distinguished the commonly confused signatures SBS3, SBS5, SBS8, SBS39, and SBS40. Importantly, we robustly identified two different multimodal SBS3 signatures, m-SBS3a and m-SBS3b, with distinct patterns in the indel spectrum. Multimodal SBS3b signature was strongly predictive of longer survival in ovarian cancer patients, replicating across four cohorts, with effect sizes greatly exceeding other genetic markers. Our m-SBS3 also predicted survival in platinum-treated patients with various cancer types, and moreover, the SBS + ID joint inference was successfully applied to mismatch repair-deficient colorectal cancer and immunotherapy response, supporting a general utility of the multimodal mutational signatures approach. Overall, combining SBS and ID mutations improves detection of HR deficiency-associated signatures and reveals distinct SBS3 subtypes with prognostic value. This multimodal approach outperforms existing markers and is readily applicable to therapy stratification.

Functional pre-therapeutic evaluation by genome editing of variants of uncertain significance of essential tumor suppressor genes

AbstractBackgroundTargeted therapies in oncology are promising but variants of uncertain significance (VUS) limit their use for clinical management and necessitate functional testing in vitro. UsingBRCA1andBRCA2variants, which have consequences on PARP inhibitor sensitivity, andPOLEvariants, potential biomarkers of immunotherapy response, we developed a rapid functional assay based on CRISPR-Cas9 genome editing to determine the functional consequences of these variants having potentially direct implications on patients’ access to targeted therapies.MethodsWe first evaluated the functional impact of 26BRCA1and 7BRCA2variants by editing and comparing NGS results between the variant of interest and a silent control variant. Ten of these variants had already been classified as benign or pathogenic and were used as controls. Finally, we extended this method to the characterization ofPOLEVUS.ResultsFor the 23 variants that were unclassified or for which conflicting interpretations had been reported, 15 were classified as functionally normal and 6 as functionally abnormal. Another two variants were found to have intermediate consequences, both with potential impacts on splicing. We then compared these scores to the patients’ responses to PARP inhibitors when possible. Finally, to prove the application of our method to the classification of variants from other tumor suppressor genes, we exemplified with threePOLEVUS. Among them, two were classified with an intermediate functional impact and one was functionally abnormal. Eventually, fourPOLEvariants previously classified in databases were also evaluated. However, we found evidence of a discordance with the classification, variant p.Leu424Val being found here functionally normal.ConclusionsOur new rapid functional assay can be used to characterize the functional implication ofBRCA1andBRCA2variants, giving patients whose variants were evaluated as functionally abnormal access to PARP inhibitor treatment. Retrospective analysis of patients’ responses to PARP inhibitors, where accessible, was consistent with our functional score evaluation and confirmed the accuracy of our protocol. This method could potentially be extended to the classification of VUS from all essential tumor suppressor genes and can be performed within a timeframe compatible with clinical applications, thereby having a direct theranostic impact.

High-grade serous tubo-ovarian cancer refined with single-cell RNA sequencing: specific cell subtypes influence survival and determine molecular subtype classification

Abstract Background High-grade serous tubo-ovarian cancer (HGSTOC) is characterised by extensive inter- and intratumour heterogeneity, resulting in persistent therapeutic resistance and poor disease outcome. Molecular subtype classification based on bulk RNA sequencing facilitates a more accurate characterisation of this heterogeneity, but the lack of strong prognostic or predictive correlations with these subtypes currently hinders their clinical implementation. Stromal admixture profoundly affects the prognostic impact of the molecular subtypes, but the contribution of stromal cells to each subtype has poorly been characterised. Increasing the transcriptomic resolution of the molecular subtypes based on single-cell RNA sequencing (scRNA-seq) may provide insights in the prognostic and predictive relevance of these subtypes. Methods We performed scRNA-seq of 18,403 cells unbiasedly collected from 7 treatment-naive HGSTOC tumours. For each phenotypic cluster of tumour or stromal cells, we identified specific transcriptomic markers. We explored which phenotypic clusters correlated with overall survival based on expression of these transcriptomic markers in microarray data of 1467 tumours. By evaluating molecular subtype signatures in single cells, we assessed to what extent a phenotypic cluster of tumour or stromal cells contributes to each molecular subtype. Results We identified 11 cancer and 32 stromal cell phenotypes in HGSTOC tumours. Of these, the relative frequency of myofibroblasts, TGF-β-driven cancer-associated fibroblasts, mesothelial cells and lymphatic endothelial cells predicted poor outcome, while plasma cells correlated with more favourable outcome. Moreover, we identified a clear cell-like transcriptomic signature in cancer cells, which correlated with worse overall survival in HGSTOC patients. Stromal cell phenotypes differed substantially between molecular subtypes. For instance, the mesenchymal, immunoreactive and differentiated signatures were characterised by specific fibroblast, immune cell and myofibroblast/mesothelial cell phenotypes, respectively. Cell phenotypes correlating with poor outcome were enriched in molecular subtypes associated with poor outcome. Conclusions We used scRNA-seq to identify stromal cell phenotypes predicting overall survival in HGSTOC patients. These stromal features explain the association of the molecular subtypes with outcome but also the latter’s weakness of clinical implementation. Stratifying patients based on marker genes specific for these phenotypes represents a promising approach to predict prognosis or response to therapy.

The vaginal microbiota and innate immunity after local excisional treatment for cervical intraepithelial neoplasia

Abstract Background Vaginal microbiota (VMB) composition is altered in women with cervical intra-epithelial neoplasia (CIN) compared to healthy controls and is associated with disease progression. However, the impact of CIN excision on the VMB and innate immunity is not known. This observational study aims to explore the impact of CIN excision on the VMB, antimicrobial peptides (AMP) and proinflammatory cytokines. Methods We sampled 103 non-pregnant, premenopausal women at the time of excisional treatment for CIN and at their 6-month follow-up visit. A further 39 untreated controls with normal cytology were also sampled. We used metataxonomics to group vaginal swab samples into community state types (CSTs) and ELISA to quantify cytokine and AMP levels in matched vaginal secretions. Analyses were performed to compare the bacterial composition and immune analyte levels before and after CIN excision and in healthy controls. Results Women with CIN had significantly higher rates of Lactobacillus species depletion pre-treatment compared to healthy controls (CST IV 21/103, 20% vs 1/39, 3%, p = 0.0081). Excision did not change the VMB composition, with CST IV remaining significantly more prevalent after excision compared to untreated, healthy controls (CST IV 19/103, 20% vs 1/39, 3%, p = 0.0142). Prevotella bivia and Sneathia amnii were significantly higher in samples before treatment compared to untreated controls, and Prevotella bivia remained significantly higher amongst the treated, with less Lactobacillus crispatus compared to untreated controls. IL-1β and IL-8 remained significantly elevated pre- (p < 0.0001 and p = 0.0014, respectively) and post-treatment (p < 0.0001 and p = 0.0035, respectively) compared to untreated controls. Levels of human beta-defensin-1 and secretory leukocyte protease inhibitor were both significantly reduced following CIN excision (p < 0.0001); however, their levels remained lower than controls post-treatment. Conclusions Women with CIN have an increased prevalence of Lactobacillus sp. depletion, high-diversity VMB composition, and higher levels of proinflammatory cytokines and AMPs compared to normal controls. Surgical excision of the disease reduces levels of vaginal AMPs but does not alter VMB composition or cytokine levels. These findings suggest that women with CIN have an inherent predisposition to a high-diversity proinflammatory environment that is not corrected by disease excision. The failure to re-establish a Lactobacillus-enriched CST may explain why women remain at high risk of pre-invasive and invasive disease recurrence.

The WID-CIN test identifies women with, and at risk of, cervical intraepithelial neoplasia grade 3 and invasive cervical cancer

AbstractBackgroundCervical screening is transitioning from primary cytology to primary human papillomavirus (HPV) testing. HPV testing is highly sensitive but there is currently no high-specificity triage method for colposcopy referral to detect cervical intraepithelial neoplasia grade 3 or above (CIN3+) in women positive for high-risk (hr) HPV subtypes. An objective, automatable test that could accurately perform triage, independently of sample heterogeneity and age, is urgently required.MethodsWe analyzed DNA methylation at ~850,000 CpG sites across the genome in a total of 1254 cervical liquid-based cytology (LBC) samples from cases of screen-detected histologically verified CIN1-3+ (98% hrHPV-positive) and population-based control women free from any cervical disease (100% hrHPV-positive). Samples were provided by a state-of-the-art population-based cohort biobank and consisted of (i) a discovery set of 170 CIN3+ cases and 202 hrHPV-positive/cytology-negative controls; (ii) a diagnostic validation set of 87 CIN3+, 90 CIN2, 166 CIN1, and 111 hrHPV-positive/cytology-negative controls; and (iii) a predictive validation set of 428 cytology-negative samples (418 hrHPV-positive) of which 210 were diagnosed with CIN3+ in the upcoming 1–4 years and 218 remained disease-free.ResultsWe developed the WID-CIN (Women’s cancer risk IDentification-Cervical Intraepithelial Neoplasia) test, a DNA methylation signature consisting of 5000 CpG sites. The receiver operating characteristic area under the curve (AUC) in the independent diagnostic validation set was 0.92 (95% CI 0.88–0.96). At 75% specificity (≤CIN1), the overall sensitivity to detect CIN3+ is 89.7% (83.3–96.1) in all and 92.7% (85.9–99.6) and 65.6% (49.2–82.1) in women aged ≥30 and <30. In hrHPV-positive/cytology-negative samples in the predictive validation set, the WID-CIN detected 54.8% (48.0–61.5) cases developing 1–4 years after sample donation in all ages or 56.9% (47.6–66.2) and 53.5% (43.7–63.2) in ≥30 and <30-year-old women, at a specificity of 75%.ConclusionsThe WID-CIN test identifies the vast majority of hrHPV-positive women with current CIN3+ lesions. In the absence of cytologic abnormalities, a positive WID-CIN test result is likely to indicate a significantly increased risk of developing CIN3+ in the near future.