The brain is a common site of metastatic disease in patients with breast cancer, which has few therapeutic options and dismal outcomes. The purpose of our study was to identify common and rare events that underlie breast cancer brain metastasis. We performed deep genomic profiling, which integrated gene copy number, gene expression and DNA methylation datasets on a collection of breast brain metastases. We identified frequent large chromosomal gains in 1q, 5p, 8q, 11q, and 20q and frequent broad-level deletions involving 8p, 17p, 21p and Xq. Frequently amplified and overexpressed genes included ATAD2, BRAF, DERL1, DNMTRB and NEK2A. The ATM, CRYAB and HSPB2 genes were commonly deleted and underexpressed. Knowledge mining revealed enrichment in cell cycle and G2/M transition pathways, which contained AURKA, AURKB and FOXM1. Using the PAM50 breast cancer intrinsic classifier, Luminal B, Her2+/ER negative, and basal-like tumors were identified as the most commonly represented breast cancer subtypes in our brain metastasis cohort. While overall methylation levels were increased in breast cancer brain metastasis, basal-like brain metastases were associated with significantly lower levels of methylation. Integrating DNA methylation data with gene expression revealed defects in cell migration and adhesion due to hypermethylation and downregulation of PENK, EDN3, and ITGAM. Hypomethylation and upregulation of KRT8 likely affects adhesion and permeability. Genomic and epigenomic profiling of breast brain metastasis has provided insight into the somatic events underlying this disease, which have potential in forming the basis of future therapeutic strategies.

Advanced cholangiocarcinoma continues to harbor a difficult prognosis and therapeutic options have been limited. During the course of a clinical trial of whole genomic sequencing seeking druggable targets, we examined six patients with advanced cholangiocarcinoma. Integrated genome-wide and whole transcriptome sequence analyses were performed on tumors from six patients with advanced, sporadic intrahepatic cholangiocarcinoma (SIC) to identify potential therapeutically actionable events. Among the somatic events captured in our analysis, we uncovered two novel therapeutically relevant genomic contexts that when acted upon, resulted in preliminary evidence of anti-tumor activity. Genome-wide structural analysis of sequence data revealed recurrent translocation events involving the FGFR2 locus in three of six assessed patients. These observations and supporting evidence triggered the use of FGFR inhibitors in these patients. In one example, preliminary anti-tumor activity of pazopanib (in vitro FGFR2 IC50≈350 nM) was noted in a patient with an FGFR2-TACC3 fusion. After progression on pazopanib, the same patient also had stable disease on ponatinib, a pan-FGFR inhibitor (in vitro, FGFR2 IC50≈8 nM). In an independent non-FGFR2 translocation patient, exome and transcriptome analysis revealed an allele specific somatic nonsense mutation (E384X) in ERRFI1, a direct negative regulator of EGFR activation. Rapid and robust disease regression was noted in this ERRFI1 inactivated tumor when treated with erlotinib, an EGFR kinase inhibitor. FGFR2 fusions and ERRFI mutations may represent novel targets in sporadic intrahepatic cholangiocarcinoma and trials should be characterized in larger cohorts of patients with these aberrations.

Genomic analyses of cutaneous melanoma (CM) have yielded biological and therapeutic insights, but understanding of non-ultraviolet (UV)-derived CMs remains limited. Deeper analysis of acral lentiginous melanoma (ALM), a rare sun-shielded melanoma subtype associated with worse survival than CM, is needed to delineate non-UV oncogenic mechanisms. We thus performed comprehensive genomic and transcriptomic analysis of 34 ALM patients. Unlike CM, somatic alterations were dominated by structural variation and absence of UV-derived mutation signatures. Only 38% of patients demonstrated driver mutations. In contrast with CM, we observed copy gains in 15% of patients, and somatic translocations, copy gains, and missense and promoter mutations, or germline events, in 41% of patients. We further show that in vitro TERT inhibition has cytotoxic effects on primary ALM cells. These findings provide insight into the role of in ALM tumorigenesis and reveal preliminary evidence that TERT inhibition represents a potential therapeutic strategy in ALM.

Brain metastasis (BM) can affect ~25% of non-small cell lung cancer (NSCLC) patients during their lifetime. Efforts to characterize patients that will develop BM have been disappointing. microRNAs (miRNAs) regulate the expression of target mRNAs. miRNAs play a role in regulating a variety of targets and, consequently, multiple pathways, which make them a powerful tool for early detection of disease, risk assessment, and prognosis. We investigated miRNAs that may serve as biomarkers to differentiate between NSCLC patients with and without BM. miRNA microarray profiling was performed on samples from clinically matched NSCLC from seven patients with BM (BM+) and six without BM (BM−). Using t-test and further qRT-PCR validation, eight miRNAs were confirmed to be significantly differentially-expressed. Of these, expression of miR-328 and miR-330-3p were able to correctly classify BM+ vs. BM− patients. This classifier was used on a validation cohort (n=15) and it correctly classified 12/15 patients. Gene expression analysis comparing A549 parental and A549 cells stably transfected to over-express miR-328 (A549–328) identified several significantly differentially-expressed genes. PRKCA was one of the genes over-expressed in A549–328 cells. Additionally, A549–328 cells had significantly increased cell migration compared to A549 cells, which was significantly reduced upon PRKCA knockdown. In summary, miR-328 has a role in conferring migratory potential to NSCLC cells working in part through PRKCA and with further corroboration in additional independent cohorts, these miRNAs may be incorporated into clinical treatment decision making to stratify NSCLC patients at higher risk for developing BM.

Bipolar disorder is a common, heritable mental illness characterized by recurrent episodes of mania and depression. Despite considerable effort to elucidate the genetic underpinnings of bipolar disorder, causative genetic risk factors remain elusive. We conducted a comprehensive genomic analysis of bipolar disorder in a large Old Order Amish pedigree. Microsatellite genotypes and high-density SNP-array genotypes of 388 family members were combined with whole genome sequence data for 50 of these subjects, comprising 18 parent-child trios. This study design permitted evaluation of candidate variants within the context of haplotype structure by resolving the phase in sequenced parent-child trios and by imputation of variants into multiple unsequenced siblings. Non-parametric and parametric linkage analysis of the entire pedigree as well as on smaller clusters of families identified several nominally significant linkage peaks, each of which included dozens of predicted deleterious variants. Close inspection of exonic and regulatory variants in genes under the linkage peaks using family-based association tests revealed additional credible candidate genes for functional studies and further replication in population-based cohorts. However, despite the in-depth genomic characterization of this unique, large and multigenerational pedigree from a genetic isolate, there was no convergence of evidence implicating a particular set of risk loci or common pathways. The striking haplotype and locus heterogeneity we observed has profound implications for the design of studies of bipolar and other related disorders.

This clinical trial evaluated whether whole exome sequencing (WES) and RNA sequencing (RNAseq) of paired normal and tumor tissues could be incorporated into a personalized treatment plan for newly diagnosed patients (<25 years of age) with diffuse intrinsic pontine glioma (DIPG). Additionally, whole genome sequencing (WGS) was compared to WES to determine if WGS would further inform treatment decisions, and whether circulating tumor DNA (ctDNA) could detect the H3K27M mutation to allow assessment of therapy response. Patients were selected across three Pacific Pediatric Neuro‐Oncology Consortium member institutions between September 2014 and January 2016. WES and RNAseq were performed at diagnosis and recurrence when possible in a CLIA‐certified laboratory. Patient‐derived cell line development was attempted for each subject. Collection of blood for ctDNA was done prior to treatment and with each MRI. A specialized tumor board generated a treatment recommendation including up to four FDA‐approved agents based upon the genomic alterations detected. A treatment plan was successfully issued within 21 business days from tissue collection for all 15 subjects, with 14 of the 15 subjects fulfilling the feasibility criteria. WGS results did not significantly deviate from WES‐based therapy recommendations; however, WGS data provided further insight into tumor evolution and fidelity of patient‐derived cell models. Detection of the H3F3A or HIST1H3B K27M (H3K27M) mutation using ctDNA was successful in 92% of H3K27M mutant cases. A personalized treatment recommendation for DIPG can be rendered within a multicenter setting using comprehensive next‐generation sequencing technology in a clinically relevant timeframe.

Large-scale multiplexed identification of somatic alterations in cancer has become feasible with next generation sequencing (NGS). However, calibration of NGS somatic analysis tools has been hampered by a lack of tumor/normal reference standards. We thus performed paired PCR-free whole genome sequencing of a matched metastatic melanoma cell line (COLO829) and normal across three lineages and across separate institutions, with independent library preparations, sequencing, and analysis. We generated mean mapped coverages of 99X for COLO829 and 103X for the paired normal across three institutions. Results were combined with previously generated data allowing for comparison to a fourth lineage on earlier NGS technology. Aggregate variant detection led to the identification of consensus variants, including key events that represent hallmark mutation types including amplified BRAF V600E, a CDK2NA small deletion, a 12 kb PTEN deletion, and a dinucleotide TERT promoter substitution. Overall, common events include >35,000 point mutations, 446 small insertion/deletions, and >6,000 genes affected by copy number changes. We present this reference to the community as an initial standard for enabling quantitative evaluation of somatic mutation pipelines across institutions.

Glioblastoma is an aggressive and molecularly heterogeneous cancer with few effective treatment options. We hypothesized that next-generation sequencing can be used to guide treatment recommendations within a clinically acceptable time frame following surgery for patients with recurrent glioblastoma. We conducted a prospective genomics-informed feasibility trial in adults with recurrent and progressive glioblastoma. Following surgical resection, genome-wide tumor/normal exome sequencing and tumor RNA sequencing were performed to identify molecular targets for potential matched therapy. A multidisciplinary molecular tumor board issued treatment recommendations based on the genomic results, blood-brain barrier penetration of the indicated therapies, drug-drug interactions, and drug safety profiles. Feasibility of generating genomics-informed treatment recommendations within 35 days of surgery was assessed. Of the 20 patients enrolled in the study, 16 patients had sufficient tumor tissue for analysis. Exome sequencing was completed for all patients, and RNA sequencing was completed for 14 patients. Treatment recommendations were provided within the study's feasibility time frame for 15 of 16 (94%) patients. Seven patients received treatment based on the tumor board recommendations. Two patients reached 12-month progression-free survival, both adhering to treatments based on the molecular profiling results. One patient remained on treatment and progression free 21 months after surgery, 3 times longer than the patient's previous time to progression. Analysis of matched nonenhancing tissue from 12 patients revealed overlapping as well as novel putatively actionable genomic alterations. Use of genome-wide molecular profiling is feasible and can be informative for guiding real-time, central nervous system-penetrant, genomics-informed treatment recommendations for patients with recurrent glioblastoma. .