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.

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.

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.

GNAS, GNAQ, and GNA11 alterations in patients with diverse cancers.

BACKGROUND: Advances in deep sequencing technology have uncovered a widespread, protumorigenic role of guanine nucleotide-binding (G protein) α (GNA) subunits, particularly GNA subunits Gs (GNAS), Gq (GNAQ), and G11 (GNA11) (GNA*), in a diverse collection of malignancies. The objectives of the current study were: 1) to determine GNA* aberration status in a cohort of 1348 patients with cancer and 2) to examine tumor mutational burden, overall survival rates, and treatment outcomes in patients with GNA*-positive tumors versus those with tumors that had wild-type GNA*. METHODS: For each patient, clinical and genomic data were collected from medical records. Next-generation sequencing was performed for each patient (range, 182-236 genes). RESULTS: Aberrations of GNA* genes were identified in a subset of patients who had 8 of the 12 cancer types examined, and a significant association was observed for appendiceal cancer and ocular melanoma (P < .0001 for both; multivariate analysis). Overall, 4.1% of the cancer population was affected. GNA* abnormalities were associated with higher numbers of co-alterations in univariate (but not multivariate) analysis and were most commonly accompanied by Aurora kinase A (AURKA), Cbl proto-oncogene (CBL), and LYN proto-oncogene (LYN) co-alterations (all P < .0001; multivariate analysis). GNA* alterations were correlated with a trend toward lower median overall survival (P = .085). The median tumor mutational burden was 4 mutations per megabase in both GNA*-altered and GNA* wild-type tumors. For this limited sample of GNA*-positive patients, longer survival was not correlated with any specific treatment regimens. CONCLUSIONS: In the current sample, the genes GNAS, GNAQ, and GNA11 were widely altered across cancer types, and these alterations often were accompanied by specific genomic abnormalities in AURKA, CBL, and LYN. Therefore, targeting GNA* alterations may require drugs that address the GNA* signal and important co-alterations. Cancer 2018;00:000-000. © 2018 American Cancer Society.

Integrated genomic DNA/RNA profiling of hematologic malignancies in the clinical setting.

The spectrum of somatic alterations in hematologic malignancies includes substitutions, insertions/deletions (indels), copy number alterations (CNAs), and a wide range of gene fusions; no current clinically available single assay captures the different types of alterations. We developed a novel next-generation sequencing-based assay to identify all classes of genomic alterations using archived formalin-fixed paraffin-embedded blood and bone marrow samples with high accuracy in a clinically relevant time frame, which is performed in our Clinical Laboratory Improvement Amendments-certified College of American Pathologists-accredited laboratory. Targeted capture of DNA/RNA and next-generation sequencing reliably identifies substitutions, indels, CNAs, and gene fusions, with similar accuracy to lower-throughput assays that focus on specific genes and types of genomic alterations. Profiling of 3696 samples identified recurrent somatic alterations that impact diagnosis, prognosis, and therapy selection. This comprehensive genomic profiling approach has proved effective in detecting all types of genomic alterations, including fusion transcripts, which increases the ability to identify clinically relevant genomic alterations with therapeutic relevance.

Enhancer decommissioning by LSD1 during embryonic stem cell differentiation.

Transcription factors and chromatin modifiers are important in the programming and reprogramming of cellular states during development. Transcription factors bind to enhancer elements and recruit coactivators and chromatin-modifying enzymes to facilitate transcription initiation. During differentiation a subset of these enhancers must be silenced, but the mechanisms underlying enhancer silencing are poorly understood. Here we show that the histone demethylase lysine-specific demethylase 1 (LSD1; ref. 5), which demethylates histone H3 on Lys 4 or Lys 9 (H3K4/K9), is essential in decommissioning enhancers during the differentiation of mouse embryonic stem cells (ESCs). LSD1 occupies enhancers of active genes that are critical for control of the state of ESCs. However, LSD1 is not essential for the maintenance of ESC identity. Instead, ESCs lacking LSD1 activity fail to differentiate fully, and ESC-specific enhancers fail to undergo the histone demethylation events associated with differentiation. At active enhancers, LSD1 is a component of the NuRD (nucleosome remodelling and histone deacetylase) complex, which contains additional subunits that are necessary for ESC differentiation. We propose that the LSD1-NuRD complex decommissions enhancers of the pluripotency program during differentiation, which is essential for the complete shutdown of the ESC gene expression program and the transition to new cell states.

Ronin/Hcf-1 binds to a hyperconserved enhancer element and regulates genes involved in the growth of embryonic stem cells.

Self-renewing embryonic stem (ES) cells have an exceptional need for timely biomass production, yet the transcriptional control mechanisms responsible for meeting this requirement are largely unknown. We report here that Ronin (Thap11), which is essential for the self-renewal of ES cells, binds with its transcriptional coregulator, Hcf-1, to a highly conserved enhancer element that previously lacked a recognized binding factor. The subset of genes bound by Ronin/Hcf-1 function primarily in transcription initiation, mRNA splicing, and cell metabolism; genes involved in cell signaling and cell development are conspicuously underrepresented in this target gene repertoire. Although Ronin/Hcf-1 represses the expression of some target genes, its activity at promoter sites more often leads to the up-regulation of genes essential to protein biosynthesis and energy production. We propose that Ronin/Hcf-1 controls a genetic program that contributes to the unimpeded growth of ES cells.