Densely interconnected transcriptional circuits control cell states in human hematopoiesis.

Though many individual transcription factors are known to regulate hematopoietic differentiation, major aspects of the global architecture of hematopoiesis remain unknown. Here, we profiled gene expression in 38 distinct purified populations of human hematopoietic cells and used probabilistic models of gene expression and analysis of cis-elements in gene promoters to decipher the general organization of their regulatory circuitry. We identified modules of highly coexpressed genes, some of which are restricted to a single lineage but most of which are expressed at variable levels across multiple lineages. We found densely interconnected cis-regulatory circuits and a large number of transcription factors that are differentially expressed across hematopoietic states. These findings suggest a more complex regulatory system for hematopoiesis than previously assumed.

Derivation of pre-X inactivation human embryonic stem cells under physiological oxygen concentrations.

The presence of two active X chromosomes (XaXa) is a hallmark of the ground state of pluripotency specific to murine embryonic stem cells (ESCs). Human ESCs (hESCs) invariably exhibit signs of X chromosome inactivation (XCI) and are considered developmentally more advanced than their murine counterparts. We describe the establishment of XaXa hESCs derived under physiological oxygen concentrations. Using these cell lines, we demonstrate that (1) differentiation of hESCs induces random XCI in a manner similar to murine ESCs, (2) chronic exposure to atmospheric oxygen is sufficient to induce irreversible XCI with minor changes of the transcriptome, (3) the Xa exhibits heavy methylation of the XIST promoter region, and (4) XCI is associated with demethylation and transcriptional activation of XIST along with H3K27-me3 deposition across the Xi. These findings indicate that the human blastocyst contains pre-X-inactivation cells and that this state is preserved in vitro through culture under physiological oxygen.

Understanding variants of unknown significance and classification of genomic alterations.

Despite recent efforts to issue clinical guidelines outlining strategies to define the pathogenicity of genomic variants, there is currently no standardized framework for which to make these assertions. This review does not present a step-by-step methodology, but rather takes a holistic approach to discuss many aspects which should be taken into consideration when determining variant pathogenicity. Categorization should be curated to reflect relevant findings within the scope of the specific medical context. Functional characterization should evaluate all available information, including results from literature reviews, different classes of genomic data repositories, and applicable computational predictive algorithms. This article further proposes a multidimensional view to infer pathogenic status from many genomic measurements across multiple axes. Notably, tumor suppressors and oncogenes exhibit fundamentally different biology which helps refine the importance of effects on splicing, mutation interactions, copy number thresholds, rearrangement annotations, germline status, and genome-wide signatures. Understanding these relevant datapoints with thoughtful perspective could aid in the reclassification of variants of unknown significance (VUS), which are ambiguously understood and currently have uncertain clinical implications. Ongoing assessments of VUS examining these relevant biological axes could lead to more accurate classification of variant pathogenicity interpretation in diagnostic oncology.

Author Correction: Enhancer decommissioning by LSD1 during embryonic stem cell differentiation.

In this Letter, the western blot for LSD1 in the right panel of Fig. 2b ('TCP +') was inadvertently duplicated from the tubulin blot immediately below. The actual tubulin western blot shows the same result, with no significant change to the levels of tubulin (see Fig. 1 of this Amendment). In addition, the western blots for LSD1 and HDAC1 of Fig. 3b and c have been corrected to include vertical black lines to delineate the juxtaposition of lanes that were non-adjacent in the original blotting experiment (see Fig. 2 of this Amendment). Supplementary Figs. 4a, 6b and 9b have also been corrected to delineate non-adjacent lanes with vertical black lines (see Supplementary Information of this Amendment). The complete raw data images from these western blotting experiments can also be found in the Supplementary Information of this Amendment. The original Letter has not been corrected.

Concordance Between Genomic Alterations Detected by Tumor and Germline Sequencing: Results from a Tertiary Care Academic Center Molecular Tumor Board.

OBJECTIVE: The majority of tumor sequencing currently performed on cancer patients does not include a matched normal control, and in cases where germline testing is performed, it is usually run independently of tumor testing. The rates of concordance between variants identified via germline and tumor testing in this context are poorly understood. We compared tumor and germline sequencing results in patients with breast, ovarian, pancreatic, and prostate cancer who were found to harbor alterations in genes associated with homologous recombination deficiency (HRD) and increased hereditary cancer risk. We then evaluated the potential for a computational somatic-germline-zygosity (SGZ) modeling algorithm to predict germline status based on tumor-only comprehensive genomic profiling (CGP) results. METHODS: A retrospective chart review was performed using an academic cancer center's databases of somatic and germline sequencing tests, and concordance between tumor and germline results was assessed. SGZ modeling from tumor-only CGP was compared to germline results to assess this method's accuracy in determining germline mutation status. RESULTS: A total of 115 patients with 146 total alterations were identified. Concordance rates between somatic and germline alterations ranged from 0% to 85.7% depending on the gene and variant classification. After correcting for differences in variant classification and filtering practices, SGZ modeling was found to have 97.2% sensitivity and 90.3% specificity for the prediction of somatic versus germline origin. CONCLUSIONS: Mutations in HRD genes identified by tumor-only sequencing are frequently germline. Providers should be aware that technical differences related to assay design, variant filtering, and variant classification can contribute to discordance between tumor-only and germline sequencing test results. In addition, SGZ modeling had high predictive power to distinguish between mutations of somatic and germline origin without the need for a matched normal control, and could potentially be considered to inform clinical decision-making.

A Novel HRD Signature Is Predictive of FOLFIRINOX Benefit in Metastatic Pancreatic Cancer.

BACKGROUND: Pancreatic cancer (PC) represents an aggressive disease with median overall survival (OS) of less than 1 year in the front-line setting. FOLFIRINOX and gemcitabine and paclitaxel (GP) are standard of care options for these patients; however, optimal selection of therapy is challenging. METHODS: Comprehensive genomic profiling was performed on 8358 PC patients. Outcomes were available for 1149 metastatic PC patients treated with 1L FOLFIRINOX or GP. A scar-based measure of HRD was called using a machine learning-based algorithm incorporating copy number and indel features. RESULTS: A scar-based HRD signature (HRDsig) was identified in 9% of patients. HRDsig significantly co-occurred with biallelic alterations in BRCA1/2, PALB2, BARD1, and RAD51C/D, but encompassed a larger population than that defined by BRCA1/BRCA2/PALB2 (9% vs. 6%). HRDsig was predictive of 1L FOLFIRNOX chemotherapy benefit with doubled OS relative to gemcitabine and paclitaxel (GP) (rwOS aHR 0.37 [0.22-0.62]), including 25% of the population with long-term (2 year+) survival in a real-world cohort of patients. Less benefit from FOLFIRINOX was observed in the HRDsig(-) population. Predictive value was seen in both the BRCA1/2/PALB2 mutant and wildtype populations, suggesting additional value to mutational profiling. CONCLUSION: A scar-based HRD biomarker was identified in a significant fraction of PC patients and is predictive of FOLFIRINOX benefit. Incorporating a biomarker like HRDsig could identify the right patients for platinum chemotherapy and potentially reduce FOLFIRINOX use by over 40%, minimizing toxicities with similar survival outcomes. Confirmatory studies should be performed.

Connecting microRNA genes to the core transcriptional regulatory circuitry of embryonic stem cells.

MicroRNAs (miRNAs) are crucial for normal embryonic stem (ES) cell self-renewal and cellular differentiation, but how miRNA gene expression is controlled by the key transcriptional regulators of ES cells has not been established. We describe here the transcriptional regulatory circuitry of ES cells that incorporates protein-coding and miRNA genes based on high-resolution ChIP-seq data, systematic identification of miRNA promoters, and quantitative sequencing of short transcripts in multiple cell types. We find that the key ES cell transcription factors are associated with promoters for miRNAs that are preferentially expressed in ES cells and with promoters for a set of silent miRNA genes. This silent set of miRNA genes is co-occupied by Polycomb group proteins in ES cells and shows tissue-specific expression in differentiated cells. These data reveal how key ES cell transcription factors promote the ES cell miRNA expression program and integrate miRNAs into the regulatory circuitry controlling ES cell identity.

An ErbB2 splice variant lacking exon 16 drives lung carcinoma.

Lung cancer causes more deaths annually than any other malignancy. A subset of non-small cell lung cancer (NSCLC) is driven by amplification and overexpression or activating mutation of the receptor tyrosine kinase (RTK) ERBB2 In some contexts, notably breast cancer, alternative splicing of ERBB2 causes skipping of exon 16, leading to the expression of an oncogenic ERBB2 isoform (ERBB2ΔEx16) that forms constitutively active homodimers. However, the broader implications of ERBB2 alternative splicing in human cancers have not been explored. Here, we have used genomic and transcriptomic analysis to identify elevated ERBB2ΔEx16 expression in a subset of NSCLC cases, as well as splicing site mutations facilitating exon 16 skipping and deletions of exon 16 in a subset of these lung tumors and in a number of other carcinomas. Supporting the potential of ERBB2ΔEx16 as a lung cancer driver, its expression transformed immortalized lung epithelial cells while a transgenic model featuring inducible ERBB2ΔEx16 specifically in the lung epithelium rapidly developed lung adenocarcinomas following transgene induction. Collectively, these observations indicate that ERBB2ΔEx16 is a lung cancer oncogene with potential clinical importance for a proportion of patients.

Association of CD274 (PD-L1) Copy Number Changes with Immune Checkpoint Inhibitor Clinical Benefit in Non-Squamous Non-Small Cell Lung Cancer.

BACKGROUND: We sought to characterize response to immune checkpoint inhibitor (ICI) in non-squamous non-small cell lung cancer (NSCLC) across various CD274 copy number gain and loss thresholds and identify an optimal cutoff. MATERIALS AND METHODS: A de-identified nationwide (US) real-world clinico-genomic database was leveraged to study 621 non-squamous NSCLC patients treated with ICI. All patients received second-line ICI monotherapy and underwent comprehensive genomic profiling as part of routine clinical care. Overall survival (OS) from start of ICI, for CD274 copy number gain and loss cohorts across varying copy number thresholds, were assessed. RESULTS: Among the 621 patients, patients with a CD274 CN greater than or equal to specimen ploidy +2 (N = 29) had a significantly higher median (m) OS when compared with the rest of the cohort (N = 592; 16.1 [8.9-37.3] vs 8.6 [7.1-10.9] months, hazard ratio (HR) = 0.6 [0.4-1.0], P-value = .05). Patients with a CD274 copy number less than specimen ploidy (N = 299) trended toward a lower mOS when compared to the rest of the cohort (N = 322; 7.5 [5.9-11.3] vs 9.6 [7.9-12.8] months, HR = 0.9 [0.7-1.1], P-value = .3). CONCLUSION: This work shows that CD274 copy number gains at varying thresholds predict different response to ICI blockade in non-squamous NSCLC. Considering these data, prospective clinical trials should further validate these findings, specifically in the context of PD-L1 IHC test results.

Prediction and characterization of diffuse large B-cell lymphoma cell-of-origin subtypes using targeted sequencing.

The aim of the present study was to determine cell of origin (COO) from a platform using a DNA-based method, COO DNA classifier (COODC). A targeted exome-sequencing platform that applies the mutational profile of a sample was used to classify COO subtype. Two major mutational signatures associated with COO were identified: Catalogue of Somatic Mutations in Cancer (COSMIC) signature 23 enriched in activated B cell (ABC) and COSMIC signature 3, which suggested increased frequency in germinal center B cell (GCB). Differential mutation signatures linked oncogenesis to mutational processes during B-cell activation, confirming the putative origin of GCB and ABC subtypes. Integrating COO with comprehensive genomic profiling enabled identification of features associated with COO and demonstrated the feasibility of determining COO without RNA.

Lay abstract To determine subtypes of diffuse large B-cell lymphoma (DLBCL), a cancer with poor survival, we aimed to identify DLBCL subtypes using DNA mutation-based tools. A targeted gene-sequencing platform, which measures the number and types of DNA mutations in a sample, was used to categorize DLBCL subtypes. Two major patterns of mutations associated with subtypes were identified: Catalogue of Somatic Mutations in Cancer (COSMIC) signature 23 and COSMIC signature 3. Differences in how the subtypes developed suggest a link between tumor developments and B cells being activated normally, confirming where the DLBCL subtypes came from. Combining this information with comprehensive genomic profiling, which determines all of the genes that a person has, allowed identification of features that are associated with DLBCL subtypes and showed that a DLBCL subtype can be determined without using RNA.

Tumor Mutational Burden as a Predictive Biomarker for Response to Immune Checkpoint Inhibitors: A Review of Current Evidence.

Treatment with immune checkpoint inhibitors (ICPIs) extends survival in a proportion of patients across multiple cancers. Tumor mutational burden (TMB)-the number of somatic mutations per DNA megabase (Mb)-has emerged as a proxy for neoantigen burden that is an independent biomarker associated with ICPI outcomes. Based on findings from recent studies, TMB can be reliably estimated using validated algorithms from next-generation sequencing assays that interrogate a sufficiently large subset of the exome as an alternative to whole-exome sequencing. Biological processes contributing to elevated TMB can result from exposure to cigarette smoke and ultraviolet radiation, from deleterious mutations in mismatch repair leading to microsatellite instability, or from mutations in the DNA repair machinery. A variety of clinical studies have shown that patients with higher TMB experience longer survival and greater response rates following treatment with ICPIs compared with those who have lower TMB levels; this includes a prospective randomized clinical trial that found a TMB threshold of ≥10 mutations per Mb to be predictive of longer progression-free survival in patients with non-small cell lung cancer. Multiple trials are underway to validate the predictive values of TMB across cancer types and in patients treated with other immunotherapies. Here we review the rationale, algorithm development methodology, and existing clinical data supporting the use of TMB as a predictive biomarker for treatment with ICPIs. We discuss emerging roles for TMB and its potential future value for stratifying patients according to their likelihood of ICPI treatment response. IMPLICATIONS FOR PRACTICE: Tumor mutational burden (TMB) is a newly established independent predictor of immune checkpoint inhibitor (ICPI) treatment outcome across multiple tumor types. Certain next-generation sequencing-based techniques allow TMB to be reliably estimated from a subset of the exome without the use of whole-exome sequencing, thus facilitating the adoption of TMB assessment in community oncology settings. Analyses of multiple clinical trials across several cancer types have demonstrated that TMB stratifies patients who are receiving ICPIs by response rate and survival. TMB, alongside other genomic biomarkers, may provide complementary information in selecting patients for ICPI-based therapies.

Biomarkers in Breast Cancer: An Integrated Analysis of Comprehensive Genomic Profiling and PD-L1 Immunohistochemistry Biomarkers in 312 Patients with Breast Cancer.

BACKGROUND: We examined the current biomarker landscape in breast cancer when programmed death-ligand 1 (PD-L1) testing is integrated with comprehensive genomic profiling (CGP). MATERIAL AND METHODS: We analyzed data from samples of 312 consecutive patients with breast carcinoma tested with both CGP and PD-L1 (SP142) immunohistochemistry (IHC) during routine clinical care. These samples were stratified into hormone receptor positive (HR+)/human epidermal growth factor receptor negative (HER2-; n = 159), HER2-positive (n = 32), and triple-negative breast cancer (TNBC) cohorts (n = 121). RESULTS: We found that in the TNBC cohort, 43% (52/121) were immunocyte PD-L1-positive, and in the HR+/HER2- cohort, 30% (48/159) had PIK3CA companion diagnostics mutations, and hence were potentially eligible for atezolizumab plus nab-paclitaxel or alpelisib plus fulvestrant, respectively. Of the remaining 212 patients, 10.4% (22/212) had a BRCA1/2 mutation, which, if confirmed by germline testing, would allow olaparib plus talazoparib therapy. Of the remaining 190 patients, 169 (88.9%) were positive for another therapy-associated marker or a marker that would potentially qualify the patient for a clinical trial. In addition, we examined the relationship between immunocyte PD-L1 positivity and different tumor mutation burden (TMB) cutoffs and found that when a TMB cutoff of ≥9 mutations per Mb was applied (cutoff determined based on prior publication), 11.6% (14/121) patients were TMB ≥9 mutations/Mb and of these, TMB ≥9 mutations per Mb, 71.4% (10/14) were also positive for PD-L1 IHC. CONCLUSION: Our integrated PD-L1 and CGP methodology identified 32% of the tested patients as potentially eligible for at least one of the two new Food and Drug Administration approved therapies, atezolizumab or alpelisib, and an additional 61.2% (191/312) had other biomarker-guided potential therapeutic options. IMPLICATIONS FOR PRACTICE: This integrated programmed death-ligand 1 immunohistochemistry and comprehensive genomic profiling methodology identified 32% of the tested patients as eligible for at least one of the two new Food and Drug Administration-approved therapies, atezolizumab or alpelisib, and an additional 61.2% (191/312) had other biomarker-guided potential therapeutic options. These findings suggest new research opportunities to evaluate the predictive utility of other commonly seen PIK3CA mutations in hormone receptor-positive breast cancers and to standardize tumor mutation burden cutoffs to evaluate its potentially predictive role in triple-negative breast cancer.

The Pan-Cancer Landscape of Coamplification of the Tyrosine Kinases KIT, KDR, and PDGFRA.

PURPOSE: Amplifications of receptor tyrosine kinases (RTKS) are therapeutic targets in multiple tumor types (e.g. HER2 in breast cancer), and amplification of the chromosome 4 segment harboring the three RTKs KIT, PDGFRA, and KDR (4q12amp) may be similarly targetable. The presence of 4q12amp has been sporadically reported in small tumor specific series but a large-scale analysis is lacking. We assess the pan-cancer landscape of 4q12amp and provide early clinical support for the feasibility of targeting this amplicon. EXPERIMENTAL DESIGN: Tumor specimens from 132,872 patients with advanced cancer were assayed with hybrid capture based comprehensive genomic profiling which assays 186-315 genes for all classes of genomic alterations, including amplifications. Baseline demographic data were abstracted, and presence of 4q12amp was defined as 6 or more copies of KIT/KDR/PDGFRA. Concurrent alterations and treatment outcomes with matched therapies were explored in a subset of cases. RESULTS: Overall 0.65% of cases harbored 4q12amp at a median copy number of 10 (range 6-344). Among cancers with >100 cases in this series, glioblastomas, angiosarcomas, and osteosarcomas were enriched for 4q12amp at 4.7%, 4.8%, and 6.4%, respectively (all p < 0.001), giving an overall sarcoma (n = 6,885) incidence of 1.9%. Among 99 pulmonary adenocarcinoma cases harboring 4q12amp, 50 (50%) lacked any other known driver of NSLCC. Four index cases plus a previously reported case on treatment with empirical TKIs monotherapy had stable disease on average exceeding 20 months. CONCLUSION: We define 4q12amp as a significant event across the pan-cancer landscape, comparable to known pan-cancer targets such as NTRK and microsatellite instability, with notable enrichment in several cancers such as osteosarcoma where standard treatment is limited. The responses to available TKIs observed in index cases strongly suggest 4q12amp is a druggable oncogenic target across cancers that warrants a focused drug development strategy. IMPLICATIONS FOR PRACTICE: Coamplification of the receptor tyrosine kinases (rtks) KIT/KDR/PDGFRA (4q12amp) is present broadly across cancers (0.65%), with enrichment in osteosarcoma and gliomas. Evidence for this amplicon having an oncogenic role is the mutual exclusivity of 4q12amp to other known drivers in 50% of pulmonary adenocarcinoma cases. Furthermore, preliminary clinical evidence for driver status comes from four index cases of patients empirically treated with commercially available tyrosine kinase inhibitors with activity against KIT/KDR/PDGFRA who had stable disease for 20 months on average. The sum of these lines of evidence suggests further clinical and preclinical investigation of 4q12amp is warranted as the possible basis for a pan-cancer drug development strategy.

Characterization of Clinical Cases of Malignant PEComa via Comprehensive Genomic Profiling of DNA and RNA.

PURPOSE: Perivascular epithelioid cell tumor (PEComa) is a rare mesenchymal soft tissue neoplasm often linked to mTOR pathway activation via TSC2 mutation. We analyzed a series of 31 consecutive metastatic PEComa (mPEComa) cases using a combined DNA/RNA hybrid capture-based comprehensive genomic profiling (CGP) assay to assess the genomic landscape of mPEComa. PATIENTS AND METHODS: Formalin-fixed, paraffin-embedded (FFPE) blocks or slides were obtained from tumors from 31 unique patients with mPEC-oma. DNA and RNA were extracted and CGP was performed on 405 genes using a targeted next-generation sequencing (NGS) assay in a CLIA-certified lab. RESULTS: All cases had locally advanced or metastatic disease, and 58% of patients were female with a median age of 50 years (range 8-76), and 17 and 14 specimens were from primary and metastatic sites, respectively. One hundred genomic alterations were identified in the cohort, with an average of 3.2 genomic alterations/case including alterations in TSC2 32.3% of cases (10), TSC1 9.6% (3), TFE3 16.1% (5, all fusions), and folliculin (FLCN) 6.4% (2), with all occurring in mutually exclusive fashion. Of TSC2 mutant cases, 70% had biallelic inactivation of this locus, as were 100% of TSC1 mutant cases. Two TSC1/2 wildtype cases harbored truncating mutations in FLCN, both of which were under LOH. Five TFE3 fusion cases were identified including the novel 5' fusion partner ZC3H4. CONCLUSIONS: We describe for the first time mPEComa cases with FLCN mutations under LOH, further characterizing dysregulation of the mTOR pathway as a unifying theme in mPEC-oma. Cumulatively, we demonstrate the feasibility and potential utility of segregating mPEComa by TSC, TFE3, and FLCN status via CGP in clinical care.

Pan-Cancer Analysis of CDK12 Loss-of-Function Alterations and Their Association with the Focal Tandem-Duplicator Phenotype.

BACKGROUND: CDK12 loss-of-function (LOF) genomic alterations are associated with focal tandem duplications (FTDs) in ovarian and prostate cancers. Because these FTDs may produce fusion-induced neoantigens (FINAs), CDK12 alteration is a candidate biomarker for immune checkpoint inhibitor sensitivity. Here we determine the prevalence of CDK12-LOF alterations and their association with FTDs across diverse tumor types. MATERIALS AND METHODS: A total of 142,133 tumor samples comprising 379 cancer types were sequenced (August 2014 to April 2018) by hybrid capture-based comprehensive genomic profiling (Foundation Medicine, Cambridge, MA) as part of routine clinical care. Results were analyzed for base substitutions, short insertions/deletions, rearrangements, and copy number alterations. CDK12-LOF genomic alterations were assessed for zygosity status and association with FTDs/focal copy number gain. RESULTS: CDK12 genomic alterations were detected in 1.1% of all cases, most frequently in prostate cancer (5.6%), but were also observed at >1% frequency in 11 cancer types. Across multiple cancer types, including prostate, gastric/esophageal, ovarian, breast, and endometrial cancer, the number of FTDs was significantly increased in CDK12-LOF versus CDK12 wild-type cases. Notably, CDK12-LOF was not consistently associated with a homologous recombination deficiency genomic signature. Quantitative assessment of CDK12-associated FTDs by measurement of single copy number gains identified novel likely deleterious CDK12 kinase-domain mutations in prostate and ovarian cancers. CONCLUSION: Detection of CDK12-LOF genomic alterations and their association with FTDs in a diverse spectrum of malignancies suggests that immunotherapy approaches targeting FINAs derived from CDK12-associated FTDs may be a broadly applicable strategy that could be explored across cancer types in a tumor-agnostic manner. IMPLICATIONS FOR PRACTICE: CDK12 inactivation in ovarian and prostate cancer results in the generation of focal tandem duplications, which can cause fusion-induced neoantigens. In prostate cancer, CDK12 alterations have demonstrated promise as a potential predictive biomarker for response to immune checkpoint blockade. This study evaluated genomic profiling data from >142,000 tumors to determine the prevalence of CDK12 loss-of-function genomic alterations across tumor types and demonstrated that CDK12 alterations are associated with the tandem-duplicator phenotype in cancer types other than ovarian and prostate cancer. The association of CDK12 alterations with focal tandem duplications across broad cancer types suggests that CDK12 inactivation warrants further investigation as a pan-cancer biomarker for immunotherapy benefit.

Analysis of DNA Damage Response Gene Alterations and Tumor Mutational Burden Across 17,486 Tubular Gastrointestinal Carcinomas: Implications for Therapy.

BACKGROUND: Alterations in the DNA damage response (DDR) pathway confer sensitivity to certain chemotherapies, radiation, and other DNA damage repair targeted therapies. BRCA1/2 are the most well-studied DDR genes, but recurrent alterations are described in other DDR pathway members across cancers. Deleterious DDR alterations may sensitize tumor cells to poly (ADP-ribose) polymerase inhibition, but there are also increasing data suggesting that there may also be synergy with immune checkpoint inhibitors. The relevance of DDR defects in gastrointestinal (GI) cancers is understudied. We sought to characterize DDR-defective GI malignancies and to explore genomic context and tumor mutational burden (TMB) to provide a platform for future rational investigations. MATERIALS AND METHODS: Tumor samples from 17,486 unique patients with advanced colorectal, gastroesophageal, or small bowel carcinomas were assayed using hybrid-capture-based comprehensive genomic profiling including sequencing of 10 predefined DDR genes: ARID1A, ATM, ATR, BRCA1, BRCA2, CDK12, CHEK1, CHEK2, PALB2, and RAD51. TMB (mutations per megabase [mut/Mb]) was calculated from up to 1.14 Mb of sequenced DNA. Clinicopathologic features were extracted and descriptive statistics were used to explore genomic relationships among identified subgroups. RESULTS: DDR alterations were found in 17% of cases: gastric adenocarcinoma 475/1,750 (27%), small bowel adenocarcinoma 148/666 (22%), esophageal adenocarcinoma 467/2,501 (19%), and colorectal cancer 1,824/12,569 (15%). ARID1A (9.2%) and ATM (4.7%) were the most commonly altered DDR genes in this series, followed by BRCA2 (2.3%), BRCA1 (1.1%), CHEK2 (1.0%), ATR (0.8%), CDK12 (0.7%), PALB2 (0.6%), CHEK1 (0.1%) and RAD51 (0.1%). More than one DDR gene alteration was found in 24% of cases. High microsatellite instability (MSI-H) and high TMB (TMB-H, ≥20 mut/Mb) were found in 19% and 21% of DDR-altered cases, respectively. Of DDR-altered/TMB-H cases, 87% were also MSI-H. However, even in the microsatellite stable (MSS)/DDR-wild-type (WT) versus MSS/DDR-altered, TMB-high was seen more frequently (0.4% vs. 3.3%, P < .00001.) Median TMB was 5.4 mut/Mb in the MSS/DDR-altered subset versus 3.8 mut/Mb in the MSS/DDR-WT subset (P ≤ .00001), and ATR alterations were enriched in the MSS/TMB-high cases. CONCLUSION: This is the largest study to examine selected DDR defects in tubular GI cancers and confirms that DDR defects are relatively common and that there is an association between the selected DDR defects and a high TMB in more than 20% of cases. Microsatellite stable DDR-defective tumors with elevated TMB warrant further exploration. IMPLICATIONS FOR PRACTICE: Deleterious DNA damage response (DDR) alterations may sensitize tumor cells to poly (ADP-ribose) polymerase inhibition, but also potentially to immune checkpoint inhibitors, owing to accumulation of mutations in DDR-defective tumors. The relevance of DDR defects in gastrointestinal (GI) cancers is understudied. This article characterizes DDR-defective GI malignancies and explores genomic context and tumor mutational burden to provide a platform for future rational investigations.

Genomic Profiling of Parathyroid Carcinoma Reveals Genomic Alterations Suggesting Benefit from Therapy.

BACKGROUND: Parathyroid carcinoma (PC) is a rare endocrine malignancy that can cause life-threatening hypercalcemia. We queried whether comprehensive genomic profiling (CGP) of PC might identify genomic alterations (GAs), which would suggest benefit from rationally matched therapeutics. METHODS: We performed hybrid-capture-based CGP to identify GAs and tumor mutational burden (TMB) in tumors from patients with this malignancy. RESULTS: There were 85 total GAs in 16 cases (5.3 GAs per case), and the median TMB was 1.7 mutations per megabase (m/Mb), with three cases having >20 m/Mb (18.7%). The genes most frequently harboring GA were CDC73 (38%), TP53 (38%), and MEN1 (31%). All MEN1-mutated cases also had loss of heterozygosity at that locus, but in contrast all CDC73-mutated cases retained heterozygosity. GAs suggesting potential benefit from matched targeted therapy were identified in 11 patients (69%) and most frequently found in PTEN (25%), NF1 (12.5%), KDR (12.5%), PIK3CA (12.5%), and TSC2 (12.5%). A patient whose tumor harbored KDR T668 K and who was treated with cabozantinib experienced a > 50% drop in parathyroid hormone level and radiographic partial response of 5.4 months with duration limited by toxicity. CONCLUSION: CGP identified GAs in PC that suggest benefit from targeted therapy, as supported by an index case of response to a matched tyrosine kinase inhibitor. Moreover, the unexpectedly high frequency of high TMB (>20 m/Mb) suggests a subset of PC may benefit from immune checkpoint inhibitors. IMPLICATIONS FOR PRACTICE: Parathyroid carcinoma (PC) is a rare endocrine malignancy that can cause life-threatening hypercalcemia. However, its molecular characteristics remain unclear, with few systemic therapeutic options available for this tumor. Hybrid-capture-based comprehensive genomic profiling of 16 primary cancers demonstrated presence of potentially actionable genomic alterations, including PTEN, NF1, KDR, PIK3CA, and TSC2, and a subset of hypermutated cancers with more than 20 mutations per megabase, the latter of which could benefit from immune checkpoint inhibitor therapy. A case benefiting from rationally matched targeted therapy for activating KDR mutation is also presented. These findings should be further investigated for their therapeutic potential.

PARP-1 activity (PAR) determines the sensitivity of cervical cancer to olaparib.

OBJECTIVES: Cervical cancer (CC) remains a major health problem worldwide. Poly (adenosine diphosphate [ADP]-ribose) polymerase (PARP) inhibitors (PARPi) have emerged as a promising class of chemotherapeutics in ovarian cancer. We explored the preclinical in vitro and in vivo activity of olaparib against multiple primary whole exome sequenced (WES) CC cells lines and xenografts. METHODS: Olaparib cell-cycle, apoptosis, homologous-recombination-deficiency (HRD), PARP trapping and cytotoxicity activity was evaluated against 9 primary CC cell lines in vitro. PARP and PAR expression were analyzed by Western blot assays. Finally, olaparib in vivo antitumor activity was tested against CC xenografts. RESULTS: While none of the cell lines demonstrated HRD, three out of 9 (33.3%) primary CC cell lines showed strong PARylation activity and demonstrated high sensitivity to olaparib in vitro treatment (cutoff IC50 values < 2 µM, p = 0.0012). Olaparib suppressed CC cell growth through cell cycle arrest in the G2/M phase and caused apoptosis (p < 0.0001). Olaparib activity in CC involved both PARP enzyme inhibition and trapping. In vivo, olaparib significantly impaired CC xenografts tumor growth (p = 0.0017) and increased overall animal survival (p = 0.008). CONCLUSIONS: A subset of CC primary cell lines is highly responsive to olaparib treatment in vitro and in vivo. High level of PARylation correlated with olaparib preclinical activity and may represent a useful biomarker for the identification of CC patients benefitting the most from PARPi.

A computational approach to distinguish somatic vs. germline origin of genomic alterations from deep sequencing of cancer specimens without a matched normal.

A key constraint in genomic testing in oncology is that matched normal specimens are not commonly obtained in clinical practice. Thus, while well-characterized genomic alterations do not require normal tissue for interpretation, a significant number of alterations will be unknown in whether they are germline or somatic, in the absence of a matched normal control. We introduce SGZ (somatic-germline-zygosity), a computational method for predicting somatic vs. germline origin and homozygous vs. heterozygous or sub-clonal state of variants identified from deep massively parallel sequencing (MPS) of cancer specimens. The method does not require a patient matched normal control, enabling broad application in clinical research. SGZ predicts the somatic vs. germline status of each alteration identified by modeling the alteration's allele frequency (AF), taking into account the tumor content, tumor ploidy, and the local copy number. Accuracy of the prediction depends on the depth of sequencing and copy number model fit, which are achieved in our clinical assay by sequencing to high depth (>500x) using MPS, covering 394 cancer-related genes and over 3,500 genome-wide single nucleotide polymorphisms (SNPs). Calls are made using a statistic based on read depth and local variability of SNP AF. To validate the method, we first evaluated performance on samples from 30 lung and colon cancer patients, where we sequenced tumors and matched normal tissue. We examined predictions for 17 somatic hotspot mutations and 20 common germline SNPs in 20,182 clinical cancer specimens. To assess the impact of stromal admixture, we examined three cell lines, which were titrated with their matched normal to six levels (10-75%). Overall, predictions were made in 85% of cases, with 95-99% of variants predicted correctly, a significantly superior performance compared to a basic approach based on AF alone. We then applied the SGZ method to the COSMIC database of known somatic variants in cancer and found >50 that are in fact more likely to be germline.

Association of Patient Characteristics and Tumor Genomics With Clinical Outcomes Among Patients With Non-Small Cell Lung Cancer Using a Clinicogenomic Database.

Importance: Data sets linking comprehensive genomic profiling (CGP) to clinical outcomes may accelerate precision medicine. Objective: To assess whether a database that combines EHR-derived clinical data with CGP can identify and extend associations in non-small cell lung cancer (NSCLC). Design, Setting, and Participants: Clinical data from EHRs were linked with CGP results for 28 998 patients from 275 US oncology practices. Among 4064 patients with NSCLC, exploratory associations between tumor genomics and patient characteristics with clinical outcomes were conducted, with data obtained between January 1, 2011, and January 1, 2018. Exposures: Tumor CGP, including presence of a driver alteration (a pathogenic or likely pathogenic alteration in a gene shown to drive tumor growth); tumor mutation burden (TMB), defined as the number of mutations per megabase; and clinical characteristics gathered from EHRs. Main Outcomes and Measures: Overall survival (OS), time receiving therapy, maximal therapy response (as documented by the treating physician in the EHR), and clinical benefit rate (fraction of patients with stable disease, partial response, or complete response) to therapy. Results: Among 4064 patients with NSCLC (median age, 66.0 years; 51.9% female), 3183 (78.3%) had a history of smoking, 3153 (77.6%) had nonsquamous cancer, and 871 (21.4%) had an alteration in EGFR, ALK, or ROS1 (701 [17.2%] with EGFR, 128 [3.1%] with ALK, and 42 [1.0%] with ROS1 alterations). There were 1946 deaths in 7 years. For patients with a driver alteration, improved OS was observed among those treated with (n = 575) vs not treated with (n = 560) targeted therapies (median, 18.6 months [95% CI, 15.2-21.7] vs 11.4 months [95% CI, 9.7-12.5] from advanced diagnosis; P < .001). TMB (in mutations/Mb) was significantly higher among smokers vs nonsmokers (8.7 [IQR, 4.4-14.8] vs 2.6 [IQR, 1.7-5.2]; P < .001) and significantly lower among patients with vs without an alteration in EGFR (3.5 [IQR, 1.76-6.1] vs 7.8 [IQR, 3.5-13.9]; P < .001), ALK (2.1 [IQR, 0.9-4.0] vs 7.0 [IQR, 3.5-13.0]; P < .001), RET (4.6 [IQR, 1.7-8.7] vs 7.0 [IQR, 2.6-13.0]; P = .004), or ROS1 (4.0 [IQR, 1.2-9.6] vs 7.0 [IQR, 2.6-13.0]; P = .03). In patients treated with anti-PD-1/PD-L1 therapies (n = 1290, 31.7%), TMB of 20 or more was significantly associated with improved OS from therapy initiation (16.8 months [95% CI, 11.6-24.9] vs 8.5 months [95% CI, 7.6-9.7]; P < .001), longer time receiving therapy (7.8 months [95% CI, 5.5-11.1] vs 3.3 months [95% CI, 2.8-3.7]; P < .001), and increased clinical benefit rate (80.7% vs 56.7%; P < .001) vs TMB less than 20. Conclusions and Relevance: Among patients with NSCLC included in a longitudinal database of clinical data linked to CGP results from routine care, exploratory analyses replicated previously described associations between clinical and genomic characteristics, between driver mutations and response to targeted therapy, and between TMB and response to immunotherapy. These findings demonstrate the feasibility of creating a clinicogenomic database derived from routine clinical experience and provide support for further research and discovery evaluating this approach in oncology.