Oncogenic chimeric transcription factors drive tumor-specific transcription, processing, and translation of silent genomic regions.

Many cancers are characterized by gene fusions encoding oncogenic chimeric transcription factors (TFs) such as EWS::FLI1 in Ewing sarcoma (EwS). Here, we find that EWS::FLI1 induces the robust expression of a specific set of novel spliced and polyadenylated transcripts within otherwise transcriptionally silent regions of the genome. These neogenes (NGs) are virtually undetectable in large collections of normal tissues or non-EwS tumors and can be silenced by CRISPR interference at regulatory EWS::FLI1-bound microsatellites. Ribosome profiling and proteomics further show that some NGs are translated into highly EwS-specific peptides. More generally, we show that hundreds of NGs can be detected in diverse cancers characterized by chimeric TFs. Altogether, this study identifies the transcription, processing, and translation of novel, specific, highly expressed multi-exonic transcripts from otherwise silent regions of the genome as a new activity of aberrant TFs in cancer.

Targeted long-read sequencing of the Ewing sarcoma 6p25.1 susceptibility locus identifies germline-somatic interactions with EWSR1-FLI1 binding.

Ewing sarcoma (EwS) is a rare bone and soft tissue malignancy driven by chromosomal translocations encoding chimeric transcription factors, such as EWSR1-FLI1, that bind GGAA motifs forming novel enhancers that alter nearby expression. We propose that germline microsatellite variation at the 6p25.1 EwS susceptibility locus could impact downstream gene expression and EwS biology. We performed targeted long-read sequencing of EwS blood DNA to characterize variation and genomic features important for EWSR1-FLI1 binding. We identified 50 microsatellite alleles at 6p25.1 and observed that EwS-affected individuals had longer alleles (>135 bp) with more GGAA repeats. The 6p25.1 GGAA microsatellite showed chromatin features of an EWSR1-FLI1 enhancer and regulated expression of RREB1, a transcription factor associated with RAS/MAPK signaling. RREB1 knockdown reduced proliferation and clonogenic potential and reduced expression of cell cycle and DNA replication genes. Our integrative analysis at 6p25.1 details increased binding of longer GGAA microsatellite alleles with acquired EWSR-FLI1 to promote Ewing sarcomagenesis by RREB1-mediated proliferation.

Study of chromatin remodeling genes implicates SMARCA4 as a putative player in oncogenesis in neuroblastoma.

In neuroblastoma (NB), genetic alterations in chromatin remodeling (CRGs) and epigenetic modifier genes (EMGs) have been described. We sought to determine their frequency and clinical impact. Whole exome (WES)/whole genome sequencing (WGS) data and targeted sequencing (TSCA®) of exonic regions of 33 CRGs/EMGs were analyzed in tumor samples from 283 NB patients, with constitutional material available for 55 patients. The frequency of CRG/EMG variations in NB cases was then compared to the Genome Aggregation Database (gnomAD). The sequencing revealed SNVs/small InDels or focal CNAs of CRGs/EMGs in 20% (56/283) of all cases, occurring at a somatic level in 4 (7.2%), at a germline level in 12 (22%) cases, whereas for the remaining cases, only tumor material could be analyzed. The most frequently altered genes were ATRX (5%), SMARCA4 (2.5%), MLL3 (2.5%) and ARID1B (2.5%). Double events (SNVs/small InDels/CNAs associated with LOH) were observed in SMARCA4 (n = 3), ATRX (n = 1) and PBRM1 (n = 1). Among the 60 variations, 24 (8.4%) targeted domains of functional importance for chromatin remodeling or highly conserved domains but of unknown function. Variations in SMARCA4 and ATRX occurred more frequently in the NB as compared to the gnomAD control cohort (OR = 4.49, 95%CI: 1.63-9.97, p = 0.038; OR 3.44, 95%CI: 1.46-6.91, p = 0.043, respectively). Cases with CRG/EMG variations showed a poorer overall survival compared to cases without variations. Genetic variations of CRGs/EMGs with likely functional impact were observed in 8.4% (24/283) of NB. Our case-control approach suggests a role of SMARCA4 as a player of NB oncogenesis.

Morphology and genomic hallmarks of breast tumours developed by ATM deleterious variant carriers.

BACKGROUND: The ataxia telangiectasia mutated (ATM) gene is a moderate-risk breast cancer susceptibility gene; germline loss-of-function variants are found in up to 3% of hereditary breast and ovarian cancer (HBOC) families who undergo genetic testing. So far, no clear histopathological and molecular features of breast tumours occurring in ATM deleterious variant carriers have been described, but identification of an ATM-associated tumour signature may help in patient management. METHODS: To characterise hallmarks of ATM-associated tumours, we performed systematic pathology review of tumours from 21 participants from ataxia-telangiectasia families and 18 participants from HBOC families, as well as copy number profiling on a subset of 23 tumours. Morphology of ATM-associated tumours was compared with that of 599 patients with no BRCA1 and BRCA2 mutations from a hospital-based series, as well as with data from The Cancer Genome Atlas. Absolute copy number and loss of heterozygosity (LOH) profiles were obtained from the OncoScan SNP array. In addition, we performed whole-genome sequencing on four tumours from ATM loss-of-function variant carriers with available frozen material. RESULTS: We found that ATM-associated tumours belong mostly to the luminal B subtype, are tetraploid and show LOH at the ATM locus at 11q22-23. Unlike tumours in which BRCA1 or BRCA2 is inactivated, tumours arising in ATM deleterious variant carriers are not associated with increased large-scale genomic instability as measured by the large-scale state transitions signature. Losses at 13q14.11-q14.3, 17p13.2-p12, 21p11.2-p11.1 and 22q11.23 were observed. Somatic alterations at these loci may therefore represent biomarkers for ATM testing and harbour driver mutations in potentially 'druggable' genes that would allow patients to be directed towards tailored therapeutic strategies. CONCLUSIONS: Although ATM is involved in the DNA damage response, ATM-associated tumours are distinct from BRCA1-associated tumours in terms of morphological characteristics and genomic alterations, and they are also distinguishable from sporadic breast tumours, thus opening up the possibility to identify ATM variant carriers outside the ataxia-telangiectasia disorder and direct them towards effective cancer risk management and therapeutic strategies.

Germline BAP1 mutations predispose to renal cell carcinomas.

The genetic cause of some familial nonsyndromic renal cell carcinomas (RCC) defined by at least two affected first-degree relatives is unknown. By combining whole-exome sequencing and tumor profiling in a family prone to cases of RCC, we identified a germline BAP1 mutation c.277A>G (p.Thr93Ala) as the probable genetic basis of RCC predisposition. This mutation segregated with all four RCC-affected relatives. Furthermore, BAP1 was found to be inactivated in RCC-affected individuals from this family. No BAP1 mutations were identified in 32 familial cases presenting with only RCC. We then screened for germline BAP1 deleterious mutations in familial aggregations of cancers within the spectrum of the recently described BAP1-associated tumor predisposition syndrome, including uveal melanoma, malignant pleural mesothelioma, and cutaneous melanoma. Among the 11 families that included individuals identified as carrying germline deleterious BAP1 mutations, 6 families presented with 9 RCC-affected individuals, demonstrating a significantly increased risk for RCC. This strongly argues that RCC belongs to the BAP1 syndrome and that BAP1 is a RCC-predisposition gene.

Circulating tumor DNA and circulating tumor cells in metastatic triple negative breast cancer patients.

Circulating tumor DNA (ctDNA) is a new circulating tumor biomarker which might be used as a prognostic biomarker in a way similar to circulating tumor cells (CTCs). Here, we used the high prevalence of TP53 mutations in triple negative breast cancer (TNBC) to compare ctDNA and CTC detection rates and prognostic value in metastatic TNBC patients. Forty patients were enrolled before starting a new line of treatment. TP53 mutations were characterized in archived tumor tissues and in plasma DNA using two next generation sequencing (NGS) platforms in parallel. Archived tumor tissue was sequenced successfully for 31/40 patients. TP53 mutations were found in 26/31 (84%) of tumor samples. The same mutation was detected in the matched plasma of 21/26 (81%) patients with an additional mutation found only in the plasma for one patient. Mutated allele fractions ranged from 2 to 70% (median 5%). The observed correlation between the two NGS approaches (R(2) = 0.903) suggested that ctDNA levels data were quantitative. Among the 27 patients with TP53 mutations, CTC count was ≥1 in 19 patients (70%) and ≥5 in 14 patients (52%). ctDNA levels had no prognostic impact on time to progression (TTP) or overall survival (OS), whereas CTC numbers were correlated with OS (p = 0.04) and marginally with TTP (p = 0.06). Performance status and elevated LDH also had significant prognostic impact. Here, absence of prognostic impact of baseline ctDNA level suggests that mechanisms of ctDNA release in metastatic TNBC may involve, beyond tumor burden, biological features that do not dramatically affect patient outcome.

A siRNA screen identifies RAD21, EIF3H, CHRAC1 and TANC2 as driver genes within the 8q23, 8q24.3 and 17q23 amplicons in breast cancer with effects on cell growth, survival and transformation.

RNA interference has boosted the field of functional genomics, by making it possible to carry out 'loss-of-function' screens in cultured cells. Here, we performed a small interfering RNA screening, in three breast cancer cell lines, for 101 candidate driver genes overexpressed in amplified breast tumors and belonging to eight amplicons on chromosomes 8q and 17q, investigating their role in cell survival/proliferation. This screening identified eight driver genes that were amplified, overexpressed and critical for breast tumor cell proliferation or survival. They included the well-described oncogenic driver genes for the 17q12 amplicon, ERBB2 and GRB7. Four of six other candidate driver genes-RAD21 and EIF3H, both on chromosome 8q23, CHRAC1 on chromosome 8q24.3 and TANC2 on chromosome 17q23-were confirmed to be driver genes regulating the proliferation/survival of clonogenic breast cancer cells presenting an amplification of the corresponding region. Indeed, knockdown of the expression of these genes decreased cell viability, through both cell cycle arrest and apoptosis induction, and inhibited the formation of colonies in anchorage-independent conditions, in soft agar. Strategies for inhibiting the expression of these genes or the function of the proteins they encode are therefore of potential value for the treatment of breast cancers presenting amplifications of the corresponding genomic region.

Polarity gene alterations in pure invasive micropapillary carcinomas of the breast.

INTRODUCTION: Pure invasive micropapillary carcinoma (IMPC) is a special type of breast carcinoma characterised by clusters of cells presenting polarity abnormalities. The biological alterations underlying this pattern remain unknown. METHODS: Pangenomic analysis (n=39), TP53 (n=43) and PIK3CA (n=41) sequencing in a series of IMPCs were performed. A subset of cases was also analysed with whole-exome sequencing (n=4) and RNA sequencing (n=6). Copy number variation profiles were compared with those of oestrogen receptors and grade-matched invasive ductal carcinomas (IDCs) of no special type. RESULTS: Unsupervised analysis of genomic data distinguished two IMPC subsets: one (Sawtooth/8/16) exhibited a significant increase in 16p gains (71%), and the other (Firestorm/Amplifier) was characterised by a high frequency of 8q (35%), 17q (20% to 46%) and 20q (23% to 30%) amplifications and 17p loss (74%). TP53 mutations (10%) were more frequently identified in the amplifier subset, and PIK3CA mutations (4%) were detected in both subsets. Compared to IDC, IMPC exhibited specific loss of the 6q16-q22 region (45%), which is associated with downregulation of FOXO3 and SEC63 gene expression. SEC63 and FOXO3 missense mutations were identified in one case each (2%). Whole-exome sequencing combined with RNA sequencing of IMPC allowed us to identify somatic mutations in genes involved in polarity, DNAH9 and FMN2 (8% and 2%, respectively) or ciliogenesis, BBS12 and BBS9 (2% each) or genes coding for endoplasmic reticulum protein, HSP90B1 and SPTLC3 (2% each) and cytoskeleton, UBR4 and PTPN21 (2% each), regardless of the genomic subset. The intracellular biological function of the mutated genes identified by gene ontology analysis suggests a driving role in the clinicopathological characteristics of IMPC. CONCLUSION: In our comprehensive molecular analysis of IMPC, we identified numerous genomic alterations without any recurrent fusion genes. Recurrent somatic mutations of genes participating in cellular polarity and shape suggest that they, together with other biological alterations (such as epigenetic modifications and stromal alterations), could contribute to the morphological pattern of IMPC. Though none of the individual abnormalities demonstrated specificity for IMPC, whether their combination in IMPC may have a cumulative effect that drives the abnormal polarity of IMPC needs to be examined further with in vitro experiments.

Somatic and germline activating mutations of the ALK kinase receptor in neuroblastoma.

Neuroblastoma, a tumour derived from the peripheral sympathetic nervous system, is one of the most frequent solid tumours in childhood. It usually occurs sporadically but familial cases are observed, with a subset of cases occurring in association with congenital malformations of the neural crest being linked to germline mutations of the PHOX2B gene. Here we conducted genome-wide comparative genomic hybridization analysis on a large series of neuroblastomas. Copy number increase at the locus encoding the anaplastic lymphoma kinase (ALK) tyrosine kinase receptor was observed recurrently. One particularly informative case presented a high-level gene amplification that was strictly limited to ALK, indicating that this gene may contribute on its own to neuroblastoma development. Through subsequent direct sequencing of cell lines and primary tumour DNAs we identified somatic mutations of the ALK kinase domain that mainly clustered in two hotspots. Germline mutations were observed in two neuroblastoma families, indicating that ALK is a neuroblastoma predisposition gene. Mutated ALK proteins were overexpressed, hyperphosphorylated and showed constitutive kinase activity. The knockdown of ALK expression in ALK-mutated cells, but also in cell lines overexpressing a wild-type ALK, led to a marked decrease of cell proliferation. Altogether, these data identify ALK as a critical player in neuroblastoma development that may hence represent a very attractive therapeutic target in this disease that is still frequently fatal with current treatments.

Sequential analysis of cfDNA reveals clonal evolution in neuroblastoma patients receiving ALK targeted therapy.

PURPOSE: The study of cell free DNA (cfDNA) enables sequential analysis of tumor cell-specific genetic alterations in neuroblastoma patients. EXPERIMENTAL DESIGN: Eighteen patients with relapsing neuroblastoma having received Lorlatinib, a 3rd generation ALK inhibitor, were identified (SACHA national registry and/or in the institution). cfDNA was analyzed at relapse for 9 patients, and sequentially for 5 patients (blood/bone marrow plasma) by performing WGS library construction followed by ALK-targeted ddPCR of the hotspot mutations (F1174L, R1275Q, I1170N) (variant allele fraction (VAF) detection limit 0.1%) and WES to evaluate disease burden and clonal evolution, following comparison with tumor/germline WES. RESULTS: Overall response rate to Lorlatinib was 33% (CI 13-59%), with response observed in 6/10 cases without versus 0/8 cases with MYCN amplification (MNA). ALK VAFs correlated with the overall clinical disease status, with a VAF<0.1% in clinical remission, versus higher VAFs (>30%) at progression. Importantly, sequential ALK ddPCR detected relapse earlier than clinical imaging. cfDNA WES revealed new SNVs, not seen in the primary tumor, in all instances of disease progression after Lorlatinib treatment, indicating clonal evolution, including alterations in genes linked to tumor aggressivity (TP53) or novel targets (EGFR). Gene pathway analysis revealed an enrichment for genes targeting cell differentiation in emerging clones, and cell adhesion in persistent clones. Evidence of clonal hematopoiesis could be observed in follow-up samples. CONCLUSION: We demonstrate the clinical utility of combining ALK cfDNA ddPCR for disease monitoring and cfDNA WES for the study of clonal evolution and resistance mechanisms in neuroblastoma patients receiving ALK targeted therapy.

Oxidative phosphorylation is a metabolic vulnerability of endocrine therapy and palbociclib resistant metastatic breast cancers.

Resistance to endocrine treatments and CDK4/6 inhibitors is considered a near-inevitability in most patients with estrogen receptor positive breast cancers (ER + BC). By genomic and metabolomics analyses of patients' tumours, metastasis-derived patient-derived xenografts (PDX) and isogenic cell lines we demonstrate that a fraction of metastatic ER + BC is highly reliant on oxidative phosphorylation (OXPHOS). Treatment by the OXPHOS inhibitor IACS-010759 strongly inhibits tumour growth in multiple endocrine and palbociclib resistant PDX. Mutations in the PIK3CA/AKT1 genes are significantly associated with response to IACS-010759. At the metabolic level, in vivo response to IACS-010759 is associated with decreased levels of metabolites of the glutathione, glycogen and pentose phosphate pathways in treated tumours. In vitro, endocrine and palbociclib resistant cells show increased OXPHOS dependency and increased ROS levels upon IACS-010759 treatment. Finally, in ER + BC patients, high expression of OXPHOS associated genes predict poor prognosis. In conclusion, these results identify OXPHOS as a promising target for treatment resistant ER + BC patients.

Reversible transitions between noradrenergic and mesenchymal tumor identities define cell plasticity in neuroblastoma.

Noradrenergic and mesenchymal identities have been characterized in neuroblastoma cell lines according to their epigenetic landscapes and core regulatory circuitries. However, their relationship and relative contribution in patient tumors remain poorly defined. We now document spontaneous and reversible plasticity between the two identities, associated with epigenetic reprogramming, in several neuroblastoma models. Interestingly, xenografts with cells from each identity eventually harbor a noradrenergic phenotype suggesting that the microenvironment provides a powerful pressure towards this phenotype. Accordingly, such a noradrenergic cell identity is systematically observed in single-cell RNA-seq of 18 tumor biopsies and 15 PDX models. Yet, a subpopulation of these noradrenergic tumor cells presents with mesenchymal features that are shared with plasticity models, indicating that the plasticity described in these models has relevance in neuroblastoma patients. This work therefore emphasizes that intrinsic plasticity properties of neuroblastoma cells are dependent upon external cues of the environment to drive cell identity.

BET and CDK Inhibition Reveal Differences in the Proliferation Control of Sympathetic Ganglion Neuroblasts and Adrenal Chromaffin Cells.

Neuroblastoma arising from the adrenal differ from ganglionic neuroblastoma both genetically and clinically, with adrenal tumors being associated with a more severe prognosis. The different tumor properties may be linked to specific tumor founder cells in adrenal and sympathetic ganglia. To address this question, we first set up cultures of mouse sympathetic neuroblasts and adrenal chromaffin cells. These cultures were then treated with various proliferation inhibitors to identify lineage-specific responses. We show that neuroblast and chromaffin cell proliferation was affected by WNT, ALK, IGF1, and PRC2/EZH2 signaling inhibitors to a similar extent. However, differential effects were observed in response to bromodomain and extraterminal (BET) protein inhibitors (JQ1, GSK1324726A) and to the CDK-7 inhibitor THZ1, with BET inhibitors preferentially affecting chromaffin cells, and THZ1 preferentially affecting neuroblasts. The differential dependence of chromaffin cells and neuroblasts on BET and CDK signaling may indicate different mechanisms during tumor initiation in sympathetic ganglia and adrenal.

Single-cell transcriptomics reveals shared immunosuppressive landscapes of mouse and human neuroblastoma.

BACKGROUND: High-risk neuroblastoma is a pediatric cancer with still a dismal prognosis, despite multimodal and intensive therapies. Tumor microenvironment represents a key component of the tumor ecosystem the complexity of which has to be accurately understood to define selective targeting opportunities, including immune-based therapies. METHODS: We combined various approaches including single-cell transcriptomics to dissect the tumor microenvironment of both a transgenic mouse neuroblastoma model and a cohort of 10 biopsies from neuroblastoma patients, either at diagnosis or at relapse. Features of related cells were validated by multicolor flow cytometry and functional assays. RESULTS: We show that the immune microenvironment of MYCN-driven mouse neuroblastoma is characterized by a low content of T cells, several phenotypes of macrophages and a population of cells expressing signatures of myeloid-derived suppressor cells (MDSCs) that are molecularly distinct from the various macrophage subsets. We document two cancer-associated fibroblasts (CAFs) subsets, one of which corresponding to CAF-S1, known to have immunosuppressive functions. Our data unravel a complex content in myeloid cells in patient tumors and further document a striking correspondence of the microenvironment populations between both mouse and human tumors. We show that mouse intratumor T cells exhibit increased expression of inhibitory receptors at the protein level. Consistently, T cells from patients are characterized by features of exhaustion, expressing inhibitory receptors and showing low expression of effector cytokines. We further functionally demonstrate that MDSCs isolated from mouse neuroblastoma have immunosuppressive properties, impairing the proliferation of T lymphocytes. CONCLUSIONS: Our study demonstrates that neuroblastoma tumors have an immunocompromised microenvironment characterized by dysfunctional T cells and accumulation of immunosuppressive cells. Our work provides a new and precious data resource to better understand the neuroblastoma ecosystem and suggest novel therapeutic strategies, targeting both tumor cells and components of the microenvironment.

A prognostic DNA signature for T1T2 node-negative breast cancer patients.

Predicting evolution of small node-negative breast carcinoma is a real challenge in clinical practice. The aim of this study was to search whether qualitative or quantitative DNA changes may help to predict metastasis of small node-negative breast carcinoma. Small invasive ductal carcinomas without axillary lymph node involvement (T1T2N0) from 168 patients with either good (111 patients with no event at 5 years after diagnosis) or poor (57 patients with early metastasis) outcome were analyzed with comparative genomic hybridization (CGH) array. A CGH classifier, identifying low- and high-risk groups of metastatic recurrence, was established in a training set of 78 patients, then validated, and compared with clinicopathological parameters in a distinct set of 90 patients. The genomic status of regions located on 2p22.2, 3p23, and 8q21-24 and the number of segmental alterations were defined in the training set to classify tumors into low- or high-risk groups. In the validation set, in addition to estrogen receptors and grade, this CGH classifier provided significant prognostic information in multivariate analysis (odds ratio, 3.34; 95% confidence interval 1.01-11.02; P = 4.78 × 10(-2), Wald test). This study shows that tumor DNA contains important prognostic information that may help to predict metastasis in T1T2N0 tumors of the breast.

Definition of Biologically Distinct Groups of Conjunctival Melanomas According to Etiological Factors and Implications for Precision Medicine.

Conjunctival melanoma (ConjMel) is a potentially deadly ocular melanoma, originating from partially sunlight-exposed mucosa. We explored the mutational landscape of ConjMel and studied the correlation with etiological factors. We collected 47 primary ConjMel samples and performed next-generation sequencing of 400 genes. Hotspot mutations in BRAF, NRAS, HRAS, and KIT were observed in 16 (34%), 5 (11%), 2, and 2 cases, respectively. Patients with BRAF and CDKN2A-mutated ConjMel tended to be younger while the NF1-mutated one tended to be older. The eight tumors arising from nevi were enriched in CTNNB1 mutations (63% vs. 8%; Fisher's exact p-test = 0.001) compared to non-nevi ConjMel and five were devoid of BRAF, RAS, NF1, or KIT mutations, suggesting a specific oncogenic process in these tumors. The two KIT-mutated cases carried SF3B1 mutations and were located on sun-protected mucosa, a genotype shared with genital and anorectal mucosal melanomas. Targetable mutations were observed in ERBB2, IDH1, MET, and MAP2K1 (one occurrence each). Mutational landscape of ConjMel characterizes distinct molecular subtypes with oncogenic drivers common with mucosal and skin melanomas. CTNNB1 mutations were associated with nevus-derived ConjMel. Concomitant KIT/SF3B1 mutations in sun-protected cases suggest a common tumorigenic process with genital and anorectal mucosal melanomas.

STAG2 mutations alter CTCF-anchored loop extrusion, reduce cis-regulatory interactions and EWSR1-FLI1 activity in Ewing sarcoma.

STAG2, a cohesin family gene, is among the most recurrently mutated genes in cancer. STAG2 loss of function (LOF) is associated with aggressive behavior in Ewing sarcoma, a childhood cancer driven by aberrant transcription induced by the EWSR1-FLI1 fusion oncogene. Here, using isogenic Ewing cells, we show that, while STAG2 LOF profoundly changes the transcriptome, it does not significantly impact EWSR1-FLI1, CTCF/cohesin, or acetylated H3K27 DNA binding patterns. In contrast, it strongly alters the anchored dynamic loop extrusion process at boundary CTCF sites and dramatically decreases promoter-enhancer interactions, particularly affecting the expression of genes regulated by EWSR1-FLI1 at GGAA microsatellite neo-enhancers. Down-modulation of cis-mediated EWSR1-FLI1 activity, observed in STAG2-LOF conditions, is associated with enhanced migration and invasion properties of Ewing cells previously observed in EWSR1-FLI1low cells. Our study illuminates a process whereby STAG2-LOF fine-tunes the activity of an oncogenic transcription factor through altered CTCF-anchored loop extrusion and cis-mediated enhancer mechanisms.

Lack of evidence for CDK12 as an ovarian cancer predisposing gene.

CDK12 variants were investigated as a genetic susceptibility to ovarian cancer in a series of 416 unrelated and consecutive patients with ovarian carcinoma and who carry neither germline BRCA1 nor BRCA2 pathogenic variant. The presence of CDK12 variants was searched in germline DNA by massive parallel sequencing on pooled DNAs. The lack of detection of deleterious variants and the observed proportion of missense variants in the series of ovarian carcinoma patients as compared with all human populations strongly suggests that CDK12 is not an ovarian cancer predisposing gene.

Transcriptional Programs Define Intratumoral Heterogeneity of Ewing Sarcoma at Single-Cell Resolution.

EWSR1-FLI1, the chimeric oncogene specific for Ewing sarcoma (EwS), induces a cascade of signaling events leading to cell transformation. However, it remains elusive how genetically homogeneous EwS cells can drive the heterogeneity of transcriptional programs. Here, we combine independent component analysis of single-cell RNA sequencing data from diverse cell types and model systems with time-resolved mapping of EWSR1-FLI1 binding sites and of open chromatin regions to characterize dynamic cellular processes associated with EWSR1-FLI1 activity. We thus define an exquisitely specific and direct enhancer-driven EWSR1-FLI1 program. In EwS tumors, cell proliferation and strong oxidative phosphorylation metabolism are associated with a well-defined range of EWSR1-FLI1 activity. In contrast, a subpopulation of cells from below and above the intermediary EWSR1-FLI1 activity is characterized by increased hypoxia. Overall, our study reveals sources of intratumoral heterogeneity within EwS tumors.

Integrated genomic and transcriptomic analysis of ductal carcinoma in situ of the breast.

PURPOSE: To gain insight into genomic and transcriptomic subtypes of ductal carcinomas in situ of the breast (DCIS). EXPERIMENTAL DESIGN: We did a combined phenotypic and genomic analysis of a series of 57 DCIS integrated with gene expression profile analysis for 26 of the 57 cases. RESULTS: Thirty-two DCIS exhibited a luminal phenotype; 21 were ERBB2 positive, and 4 were ERBB2/estrogen receptor (ER) negative with 1 harboring a bona fide basal-like phenotype. Based on a CGH analysis, genomic types were identified in this series of DCIS with the 1q gain/16q loss combination observed in 3 luminal DCIS, the mixed amplifier pattern including all ERBB2, 12 luminal and 2 ERBB2(-)/ER(-) DCIS, and the complex copy number alteration profile encompassing 14 luminal and 1 ERBB2(-)/ER(-) DCIS. Eight cases (8 of 57; 14%) presented a TP53 mutation, all being amplifiers. Unsupervised analysis of gene expression profiles of 26 of the 57 DCIS showed that luminal and ERBB2-amplified, ER-negative cases clustered separately. We further investigated the effect of high and low copy number changes on gene expression. Strikingly, amplicons but also low copy number changes especially on 1q, 8q, and 16q in DCIS regulated the expression of a subset of genes in a very similar way to that recently described in invasive ductal carcinomas. CONCLUSIONS: These combined approaches show that the molecular heterogeneity of breast ductal carcinomas exists already in in situ lesions and further indicate that DCIS and invasive ductal carcinomas share genomic alterations with a similar effect on gene expression profile.