Somatic mutations of esophageal adenocarcinoma: a comparison between Black and White patients.

Esophageal adenocarcinoma is the most common histological subtype of esophageal cancer in Western countries and shows poor prognosis with rapid growth. EAC is characterized by a strong male predominance and racial disparity. EAC is up to fivefold more common among Whites than Blacks, yet Black patients with EAC have poorer survival rates. The racial disparity remains largely unknown, and there is limited knowledge of mutations in EAC regarding racial disparities. We used whole-exome sequencing to show somatic mutation profiles derived from tumor samples from 18 EAC male patients. We identified three molecular subgroups based on the pre-defined esophageal cancer-specific mutational signatures. Group 1 is associated with age and NTHL1 deficiency-related signatures. Group 2 occurs primarily in Black patients and is associated with signatures related to DNA damage from oxidative stress and NTHL1 deficiency-related signatures. Group 3 is associated with defective homologous recombination-based DNA often caused by BRCA mutation in White patients. We observed significantly mutated race related genes (LCE2B in Black, SDR39U1 in White) were (q-value < 0.1). Our findings underscore the possibility of distinct molecular mutation patterns in EAC among different races. Further studies are needed to validate our findings, which could contribute to precision medicine in EAC.

Inter and intra-host diversity of RSV in hematopoietic stem cell transplant adults with normal and delayed viral clearance.

Respiratory syncytial virus (RSV) infection in immunocompromised individuals often leads to prolonged illness, progression to severe lower respiratory tract infection, and even death. How the host immune environment of the hematopoietic stem cell transplant (HCT) adults can affect viral genetic variation during an acute infection is not understood well. In the present study, we performed whole genome sequencing of RSV/A or RSV/B from samples collected longitudinally from HCT adults with normal (<14 days) and delayed (≥14 days) RSV clearance who were enrolled in a ribavirin trial. We determined the inter-host and intra-host genetic variation of RSV and the effect of mutations on putative glycosylation sites. The inter-host variation of RSV is centered in the attachment (G) and fusion (F) glycoprotein genes followed by polymerase (L) and matrix (M) genes. Interestingly, the overall genetic variation was constant between normal and delayed clearance groups for both RSV/A and RSV/B. Intra-host variation primarily occurred in the G gene followed by non-structural protein (NS1) and L genes; however, gain or loss of stop codons and frameshift mutations appeared only in the G gene and only in the delayed viral clearance group. Potential gain or loss of O-linked glycosylation sites in the G gene occurred both in RSV/A and RSV/B isolates. For RSV F gene, loss of N-linked glycosylation site occurred in three RSV/B isolates within an antigenic epitope. Both oral and aerosolized ribavirin did not cause any mutations in the L gene. In summary, prolonged viral shedding and immune deficiency resulted in RSV variation, especially in structural mutations in the G gene, possibly associated with immune evasion. Therefore, sequencing and monitoring of RSV isolates from immunocompromised patients are crucial as they can create escape mutants that can impact the effectiveness of upcoming vaccines and treatments.

Proteogenomic Markers of Chemotherapy Resistance and Response in Triple-Negative Breast Cancer.

Microscaled proteogenomics was deployed to probe the molecular basis for differential response to neoadjuvant carboplatin and docetaxel combination chemotherapy for triple-negative breast cancer (TNBC). Proteomic analyses of pretreatment patient biopsies uniquely revealed metabolic pathways, including oxidative phosphorylation, adipogenesis, and fatty acid metabolism, that were associated with resistance. Both proteomics and transcriptomics revealed that sensitivity was marked by elevation of DNA repair, E2F targets, G2-M checkpoint, interferon-gamma signaling, and immune-checkpoint components. Proteogenomic analyses of somatic copy-number aberrations identified a resistance-associated 19q13.31-33 deletion where LIG1, POLD1, and XRCC1 are located. In orthogonal datasets, LIG1 (DNA ligase I) gene deletion and/or low mRNA expression levels were associated with lack of pathologic complete response, higher chromosomal instability index (CIN), and poor prognosis in TNBC, as well as carboplatin-selective resistance in TNBC preclinical models. Hemizygous loss of LIG1 was also associated with higher CIN and poor prognosis in other cancer types, demonstrating broader clinical implications. SIGNIFICANCE: Proteogenomic analysis of triple-negative breast tumors revealed a complex landscape of chemotherapy response associations, including a 19q13.31-33 somatic deletion encoding genes serving lagging-strand DNA synthesis (LIG1, POLD1, and XRCC1), that correlate with lack of pathologic response, carboplatin-selective resistance, and, in pan-cancer studies, poor prognosis and CIN. This article is highlighted in the In This Issue feature, p. 2483.

Genome Sequencing in the Parkinson Disease Clinic.

Background and Objectives: Genetic variants affect both Parkinson disease (PD) risk and manifestations. Although genetic information is of potential interest to patients and clinicians, genetic testing is rarely performed during routine PD clinical care. The goal of this study was to examine interest in comprehensive genetic testing among patients with PD and document reactions to possible findings from genome sequencing in 2 academic movement disorder clinics. Methods: In 203 subjects with PD (age = 63 years, 67% male), genome sequencing was performed and filtered using a custom panel, including 49 genes associated with PD, parkinsonism, or related disorders, as well as a 90-variant PD genetic risk score. Based on the results, 231 patients (age = 67 years, 63% male) were surveyed on interest in genetic testing and responses to vignettes covering (1) familial risk of PD (LRRK2); (2) risk of PD dementia (GBA); (3) PD genetic risk score; and (4) secondary, medically actionable variants (BRCA1). Results: Genome sequencing revealed a LRRK2 variant in 3% and a GBA risk variant in 10% of our clinical sample. The genetic risk score was normally distributed, identifying 41 subjects with a high risk of PD. Medically actionable findings were discovered in 2 subjects (1%). In our survey, the majority (82%) responded that they would share a LRRK2 variant with relatives. Most registered unchanged or increased interest in testing when confronted with a potential risk for dementia or medically actionable findings, and most (75%) expressed interest in learning their PD genetic risk score. Discussion: Our results highlight broad interest in comprehensive genetic testing among patients with PD and may facilitate integration of genome sequencing in clinical practice.

Polygenic transcriptome risk scores for COPD and lung function improve cross-ethnic portability of prediction in the NHLBI TOPMed program.

While polygenic risk scores (PRSs) enable early identification of genetic risk for chronic obstructive pulmonary disease (COPD), predictive performance is limited when the discovery and target populations are not well matched. Hypothesizing that the biological mechanisms of disease are shared across ancestry groups, we introduce a PrediXcan-derived polygenic transcriptome risk score (PTRS) to improve cross-ethnic portability of risk prediction. We constructed the PTRS using summary statistics from application of PrediXcan on large-scale GWASs of lung function (forced expiratory volume in 1 s [FEV1] and its ratio to forced vital capacity [FEV1/FVC]) in the UK Biobank. We examined prediction performance and cross-ethnic portability of PTRS through smoking-stratified analyses both on 29,381 multi-ethnic participants from TOPMed population/family-based cohorts and on 11,771 multi-ethnic participants from TOPMed COPD-enriched studies. Analyses were carried out for two dichotomous COPD traits (moderate-to-severe and severe COPD) and two quantitative lung function traits (FEV1 and FEV1/FVC). While the proposed PTRS showed weaker associations with disease than PRS for European ancestry, the PTRS showed stronger association with COPD than PRS for African Americans (e.g., odds ratio [OR] = 1.24 [95% confidence interval [CI]: 1.08-1.43] for PTRS versus 1.10 [0.96-1.26] for PRS among heavy smokers with ≥ 40 pack-years of smoking) for moderate-to-severe COPD. Cross-ethnic portability of the PTRS was significantly higher than the PRS (paired t test p < 2.2 × 10-16 with portability gains ranging from 5% to 28%) for both dichotomous COPD traits and across all smoking strata. Our study demonstrates the value of PTRS for improved cross-ethnic portability compared to PRS in predicting COPD risk.

Transcriptional Reprogramming Differentiates Active from Inactive ESR1 Fusions in Endocrine Therapy-Refractory Metastatic Breast Cancer.

Genomic analysis has recently identified multiple ESR1 gene translocations in estrogen receptor alpha-positive (ERα+) metastatic breast cancer (MBC) that encode chimeric proteins whereby the ESR1 ligand binding domain (LBD) is replaced by C-terminal sequences from many different gene partners. Here we functionally screened 15 ESR1 fusions and identified 10 that promoted estradiol-independent cell growth, motility, invasion, epithelial-to-mesenchymal transition, and resistance to fulvestrant. RNA sequencing identified a gene expression pattern specific to functionally active ESR1 gene fusions that was subsequently reduced to a diagnostic 24-gene signature. This signature was further examined in 20 ERα+ patient-derived xenografts and in 55 ERα+ MBC samples. The 24-gene signature successfully identified cases harboring ESR1 gene fusions and also accurately diagnosed the presence of activating ESR1 LBD point mutations. Therefore, the 24-gene signature represents an efficient approach to screening samples for the presence of diverse somatic ESR1 mutations and translocations that drive endocrine treatment failure in MBC. SIGNIFICANCE: This study identifies a gene signature diagnostic for functional ESR1 fusions that drive poor outcome in advanced breast cancer, which could also help guide precision medicine approaches in patients harboring ESR1 mutations.

Sequencing of a central nervous system tumor demonstrates cancer transmission in an organ transplant.

Four organ transplant recipients from an organ donor diagnosed with anaplastic pleomorphic xanthoastrocytoma developed fatal malignancies for which the origin could not be confirmed by standard methods. We identified the somatic mutational profiles of the neoplasms using next-generation sequencing technologies and tracked the relationship between the different samples. The data were consistent with the presence of an aggressive clonal entity in the donor and the subsequent proliferation of descendent tumors in each recipient. Deleterious mutations in BRAF, PIK3CA, SDHC, DDR2, and FANCD2, and a chromosomal deletion spanning the CDKN2A/B genes, were shared between the recipients' lesions. In addition to demonstrating that DNA sequencing tracked a donor/recipient cancer transmission, this study established that the genetic profile of a donor tumor and its potential aggressive phenotype could have been determined before transplantation was considered. As the genetic correlates of tumor invasion and metastases become better known, adding genetic profiling by DNA sequencing to the data considered for transplant safety should be considered.

Molecular Features of Cancers Exhibiting Exceptional Responses to Treatment.

A small fraction of cancer patients with advanced disease survive significantly longer than patients with clinically comparable tumors. Molecular mechanisms for exceptional responses to therapy have been identified by genomic analysis of tumor biopsies from individual patients. Here, we analyzed tumor biopsies from an unbiased cohort of 111 exceptional responder patients using multiple platforms to profile genetic and epigenetic aberrations as well as the tumor microenvironment. Integrative analysis uncovered plausible mechanisms for the therapeutic response in nearly a quarter of the patients. The mechanisms were assigned to four broad categories-DNA damage response, intracellular signaling, immune engagement, and genetic alterations characteristic of favorable prognosis-with many tumors falling into multiple categories. These analyses revealed synthetic lethal relationships that may be exploited therapeutically and rare genetic lesions that favor therapeutic success, while also providing a wealth of testable hypotheses regarding oncogenic mechanisms that may influence the response to cancer therapy.

Genomic Profiling of Childhood Tumor Patient-Derived Xenograft Models to Enable Rational Clinical Trial Design.

Accelerating cures for children with cancer remains an immediate challenge as a result of extensive oncogenic heterogeneity between and within histologies, distinct molecular mechanisms evolving between diagnosis and relapsed disease, and limited therapeutic options. To systematically prioritize and rationally test novel agents in preclinical murine models, researchers within the Pediatric Preclinical Testing Consortium are continuously developing patient-derived xenografts (PDXs)-many of which are refractory to current standard-of-care treatments-from high-risk childhood cancers. Here, we genomically characterize 261 PDX models from 37 unique pediatric cancers; demonstrate faithful recapitulation of histologies and subtypes; and refine our understanding of relapsed disease. In addition, we use expression signatures to classify tumors for TP53 and NF1 pathway inactivation. We anticipate that these data will serve as a resource for pediatric oncology drug development and will guide rational clinical trial design for children with cancer.

Molecular profiling predicts meningioma recurrence and reveals loss of DREAM complex repression in aggressive tumors.

Meningiomas account for one-third of all primary brain tumors. Although typically benign, about 20% of meningiomas are aggressive, and despite the rigor of the current histopathological classification system there remains considerable uncertainty in predicting tumor behavior. Here, we analyzed 160 tumors from all 3 World Health Organization (WHO) grades (I through III) using clinical, gene expression, and sequencing data. Unsupervised clustering analysis identified 3 molecular types (A, B, and C) that reliably predicted recurrence. These groups did not directly correlate with the WHO grading system, which classifies more than half of the tumors in the most aggressive molecular type as benign. Transcriptional and biochemical analyses revealed that aggressive meningiomas involve loss of the repressor function of the DREAM complex, which results in cell-cycle activation; only tumors in this category tend to recur after full resection. These findings should improve our ability to predict recurrence and develop targeted treatments for these clinically challenging tumors.

Chemistry-First Approach for Nomination of Personalized Treatment in Lung Cancer.

Diversity in the genetic lesions that cause cancer is extreme. In consequence, a pressing challenge is the development of drugs that target patient-specific disease mechanisms. To address this challenge, we employed a chemistry-first discovery paradigm for de novo identification of druggable targets linked to robust patient selection hypotheses. In particular, a 200,000 compound diversity-oriented chemical library was profiled across a heavily annotated test-bed of >100 cellular models representative of the diverse and characteristic somatic lesions for lung cancer. This approach led to the delineation of 171 chemical-genetic associations, shedding light on the targetability of mechanistic vulnerabilities corresponding to a range of oncogenotypes present in patient populations lacking effective therapy. Chemically addressable addictions to ciliogenesis in TTC21B mutants and GLUT8-dependent serine biosynthesis in KRAS/KEAP1 double mutants are prominent examples. These observations indicate a wealth of actionable opportunities within the complex molecular etiology of cancer.

Novel CYP2A6 diplotypes identified through next-generation sequencing are associated with in-vitro and in-vivo nicotine metabolism.

OBJECTIVES: Smoking patterns and cessation rates vary widely across smokers and can be influenced by variation in rates of nicotine metabolism [i.e. cytochrome P450 2A6 (CYP2A6), enzyme activity]. There is high heritability of CYP2A6-mediated nicotine metabolism (60-80%) owing to known and unidentified genetic variation in the CYP2A6 gene. We aimed to identify and characterize additional genetic variants at the CYP2A6 gene locus. METHODS: A new CYP2A6-specific sequencing method was used to investigate genetic variation in CYP2A6. Novel variants were characterized in a White human liver bank that has been extensively phenotyped for CYP2A6. Linkage and haplotype structure for the novel single nucleotide polymorphisms (SNPs) were assessed. The association between novel five-SNP diplotypes and nicotine metabolism rate was investigated. RESULTS: Seven high-frequency (minor allele frequencies ≥6%) noncoding SNPs were identified as important contributors to CYP2A6 phenotypes in a White human liver bank (rs57837628, rs7260629, rs7259706, rs150298687 (also denoted rs4803381), rs56113850, rs28399453, and rs8192733), accounting for two times more variation in in-vitro CYP2A6 activity relative to the four established functional CYP2A6 variants that are frequently tested in Whites (CYP2A6*2, *4, *9, and *12). Two pairs of novel SNPs were in high linkage disequilibrium, allowing us to establish five-SNP diplotypes that were associated with CYP2A6 enzyme activity (rate of nicotine metabolism) in-vitro in the liver bank and in-vivo among smokers. CONCLUSION: The novel five-SNP diplotype may be useful to incorporate into CYP2A6 genotype models for personalized prediction of nicotine metabolism rate, cessation success, and response to pharmacotherapies.

The house spider genome reveals an ancient whole-genome duplication during arachnid evolution.

BACKGROUND: The duplication of genes can occur through various mechanisms and is thought to make a major contribution to the evolutionary diversification of organisms. There is increasing evidence for a large-scale duplication of genes in some chelicerate lineages including two rounds of whole genome duplication (WGD) in horseshoe crabs. To investigate this further, we sequenced and analyzed the genome of the common house spider Parasteatoda tepidariorum. RESULTS: We found pervasive duplication of both coding and non-coding genes in this spider, including two clusters of Hox genes. Analysis of synteny conservation across the P. tepidariorum genome suggests that there has been an ancient WGD in spiders. Comparison with the genomes of other chelicerates, including that of the newly sequenced bark scorpion Centruroides sculpturatus, suggests that this event occurred in the common ancestor of spiders and scorpions, and is probably independent of the WGDs in horseshoe crabs. Furthermore, characterization of the sequence and expression of the Hox paralogs in P. tepidariorum suggests that many have been subject to neo-functionalization and/or sub-functionalization since their duplication. CONCLUSIONS: Our results reveal that spiders and scorpions are likely the descendants of a polyploid ancestor that lived more than 450 MYA. Given the extensive morphological diversity and ecological adaptations found among these animals, rivaling those of vertebrates, our study of the ancient WGD event in Arachnopulmonata provides a new comparative platform to explore common and divergent evolutionary outcomes of polyploidization events across eukaryotes.

Recurrent internal tandem duplications of BCOR in clear cell sarcoma of the kidney.

The X-linked BCL-6 co-repressor (BCOR) gene encodes a key constituent of a variant polycomb repressive complex (PRC) that is mutated or translocated in human cancers. Here we report on the identification of somatic internal tandem duplications (ITDs) clustering in the C terminus of BCOR in 23 of 27 (85%) pediatric clear cell sarcomas of the kidney (CCSK) from two independent cohorts. We profile CCSK tumours using a combination of whole-exome, transcriptome and targeted sequencing. Identical ITD mutations are found in primary and relapsed tumour pairs but not in adjacent normal kidney or blood. Mutant BCOR transcripts and proteins are markedly upregulated in ITD-positive tumours. Transcriptome analysis of ITD-positive CCSKs reveals enrichment for PRC2-regulated genes and similarity to undifferentiated sarcomas harbouring BCOR-CCNB3 fusions. The discovery of recurrent BCOR ITDs defines a major oncogenic event in this childhood sarcoma with significant implications for diagnostic and therapeutic approaches to this tumour.

Genomic profiling of Sézary syndrome identifies alterations of key T cell signaling and differentiation genes.

Sézary syndrome is a rare leukemic form of cutaneous T cell lymphoma characterized by generalized redness, scaling, itching and increased numbers of circulating atypical T lymphocytes. It is rarely curable, with poor prognosis. Here we present a multiplatform genomic analysis of 37 patients with Sézary syndrome that implicates dysregulation of cell cycle checkpoint and T cell signaling. Frequent somatic alterations were identified in TP53, CARD11, CCR4, PLCG1, CDKN2A, ARID1A, RPS6KA1 and ZEB1. Activating CCR4 and CARD11 mutations were detected in nearly one-third of patients. ZEB1, encoding a transcription repressor essential for T cell differentiation, was deleted in over one-half of patients. IL32 and IL2RG were overexpressed in nearly all cases. Our results demonstrate profound disruption of key signaling pathways in Sézary syndrome and suggest potential targets for new therapies.

Rise and fall of subclones from diagnosis to relapse in pediatric B-acute lymphoblastic leukaemia.

There is incomplete understanding of genetic heterogeneity and clonal evolution during cancer progression. Here we use deep whole-exome sequencing to describe the clonal architecture and evolution of 20 pediatric B-acute lymphoblastic leukaemias from diagnosis to relapse. We show that clonal diversity is comparable at diagnosis and relapse and clonal survival from diagnosis to relapse is not associated with mutation burden. Six pathways were frequently mutated, with NT5C2, CREBBP, WHSC1, TP53, USH2A, NRAS and IKZF1 mutations enriched at relapse. Half of the leukaemias had multiple subclonal mutations in a pathway or gene at diagnosis, but mostly with only one, usually minor clone, surviving therapy to acquire additional mutations and become the relapse founder clone. Relapse-specific mutations in NT5C2 were found in nine cases, with mutations in four cases being in descendants of the relapse founder clone. These results provide important insights into the genetic basis of treatment failure in ALL and have implications for the early detection of mutations driving relapse.

BCOR-CCNB3 fusions are frequent in undifferentiated sarcomas of male children.

The BCOR-CCNB3 fusion gene, resulting from a chromosome X paracentric inversion, was recently described in translocation-negative 'Ewing-like' sarcomas arising in bone and soft tissue. Genetic subclassification of undifferentiated unclassified sarcomas may potentially offer markers for reproducible diagnosis and substrates for therapy. Using whole transcriptome paired-end RNA sequencing (RNA-seq) we unexpectedly identified BCOR-CCNB3 fusion transcripts in an undifferentiated spindle-cell sarcoma. RNA-seq results were confirmed through direct RT-PCR of tumor RNA and cloning of the genomic breakpoints from tumor DNA. Five additional undifferentiated sarcomas with BCOR-CCNB3 fusions were identified in a series of 42 pediatric and adult unclassified sarcomas. Genomic breakpoint analysis demonstrated unique breakpoint locations in each case at the DNA level even though the resulting fusion mRNA was identical in all cases. All patients with BCOR-CCNB3 sarcoma were males diagnosed in mid childhood (7-13 years of age). Tumors were equally distributed between axial and extra-axial locations. Five of the six tumors were soft-tissue lesions with either predominant spindle-cell morphology or spindle-cell areas interspersed with ovoid to round cells. CCNB3 immunohistochemistry showed strong nuclear positivity in five tumors before oncologic therapy, but was patchy to negative in post-treatment tumor samples. An RT-PCR assay developed to detect the fusion transcript in archival formalin-fixed tissue was positive in all six cases, with high sensitivity and specificity in both pre- and post-treated samples. This study adds to recent reports on the clinicopathologic spectrum of BCOR-CCNB3 fusion-positive sarcomas, a newly emerging entity within the undifferentiated unclassified sarcoma category and describes a simple RT-PCR assay that in conjunction with CCNB3 immunohistochemistry can be useful in diagnosing these tumors.

Mutational landscape of aggressive cutaneous squamous cell carcinoma.

PURPOSE: Aggressive cutaneous squamous cell carcinoma (cSCC) is often a disfiguring and lethal disease. Very little is currently known about the mutations that drive aggressive cSCC. EXPERIMENTAL DESIGN: Whole-exome sequencing was performed on 39 cases of aggressive cSCC to identify driver genes and novel therapeutic targets. Significantly, mutated genes were identified with MutSig or complementary methods developed to specifically identify candidate tumor suppressors based upon their inactivating mutation bias. RESULTS: Despite the very high-mutational background caused by UV exposure, 23 candidate drivers were identified, including the well-known cancer-associated genes TP53, CDKN2A, NOTCH1, AJUBA, HRAS, CASP8, FAT1, and KMT2C (MLL3). Three novel candidate tumor suppressors with putative links to cancer or differentiation, NOTCH2, PARD3, and RASA1, were also identified as possible drivers in cSCC. KMT2C mutations were associated with poor outcome and increased bone invasion. CONCLUSIONS: The mutational spectrum of cSCC is similar to that of head and neck squamous cell carcinoma and dominated by tumor-suppressor genes. These results improve the foundation for understanding this disease and should aid in identifying and treating aggressive cSCC.

Exonuclease mutations in DNA polymerase epsilon reveal replication strand specific mutation patterns and human origins of replication.

Tumors with somatic mutations in the proofreading exonuclease domain of DNA polymerase epsilon (POLE-exo*) exhibit a novel mutator phenotype, with markedly elevated TCT→TAT and TCG→TTG mutations and overall mutation frequencies often exceeding 100 mutations/Mb. Here, we identify POLE-exo* tumors in numerous cancers and classify them into two groups, A and B, according to their mutational properties. Group A mutants are found only in POLE, whereas Group B mutants are found in POLE and POLD1 and appear to be nonfunctional. In Group A, cell-free polymerase assays confirm that mutations in the exonuclease domain result in high mutation frequencies with a preference for C→A mutation. We describe the patterns of amino acid substitutions caused by POLE-exo* and compare them to other tumor types. The nucleotide preference of POLE-exo* leads to increased frequencies of recurrent nonsense mutations in key tumor suppressors such as TP53, ATM, and PIK3R1. We further demonstrate that strand-specific mutation patterns arise from some of these POLE-exo* mutants during genome duplication. This is the first direct proof of leading strand-specific replication by human POLE, which has only been demonstrated in yeast so far. Taken together, the extremely high mutation frequency and strand specificity of mutations provide a unique identifier of eukaryotic origins of replication.

Mutually exclusive recurrent somatic mutations in MAP2K1 and BRAF support a central role for ERK activation in LCH pathogenesis.

Langerhans cell histiocytosis (LCH) is a myeloproliferative disorder characterized by lesions composed of pathological CD207(+) dendritic cells with an inflammatory infiltrate. BRAFV600E remains the only recurrent mutation reported in LCH. In order to evaluate the spectrum of somatic mutations in LCH, whole exome sequencing was performed on matched LCH and normal tissue samples obtained from 41 patients. Lesions from other histiocytic disorders, juvenile xanthogranuloma, Erdheim-Chester disease, and Rosai-Dorfman disease were also evaluated. All of the lesions from histiocytic disorders were characterized by an extremely low overall rate of somatic mutations. Notably, 33% (7/21) of LCH cases with wild-type BRAF and none (0/20) with BRAFV600E harbored somatic mutations in MAP2K1 (6 in-frame deletions and 1 missense mutation) that induced extracellular signal-regulated kinase (ERK) phosphorylation in vitro. Single cases of somatic mutations of the mitogen-activated protein kinase (MAPK) pathway genes ARAF and ERBB3 were also detected. The ability of MAPK pathway inhibitors to suppress MAPK kinase and ERK phosphorylation in cell culture and primary tumor models was dependent on the specific LCH mutation. The findings of this study support a model in which ERK activation is a universal end point in LCH arising from pathological activation of upstream signaling proteins.