Author Correction: Interaction between the microbiome and TP53 in human lung cancer.

Following publication of the original paper [1], the authors submitted a new Additional file 5 to replace the one containing formatting issues. The updated Additional file 5 is published in this correction.

Author Correction: Recurrent PTPRT/JAK2 mutations in lung adenocarcinoma among African Americans.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Recurrent PTPRT/JAK2 mutations in lung adenocarcinoma among African Americans.

Reducing or eliminating persistent disparities in lung cancer incidence and survival has been challenging because our current understanding of lung cancer biology is derived primarily from populations of European descent. Here we show results from a targeted sequencing panel using NCI-MD Case Control Study patient samples and reveal a significantly higher prevalence of PTPRT and JAK2 mutations in lung adenocarcinomas among African Americans compared with European Americans. This increase in mutation frequency was validated with independent WES data from the NCI-MD Case Control Study and TCGA. We find that patients carrying these mutations have a concomitant increase in IL-6/STAT3 signaling and miR-21 expression. Together, these findings suggest the identification of these potentially actionable mutations could have clinical significance for targeted therapy and the enrollment of minority populations in clinical trials.

Rapid and reversible suppression of ALT by DAXX in osteosarcoma cells.

Many tumors maintain chromosome-ends through a telomerase-independent, DNA-templated mechanism called alternative lengthening of telomeres (ALT). While ALT occurs in only a subset of tumors, it is strongly associated with mutations in the genes ATRX and DAXX, which encode components of an H3.3 histone chaperone complex. The role of ATRX and DAXX mutations in potentiating the mechanism of ALT remains incompletely understood. Here we characterize an osteosarcoma cell line, G292, with wild-type ATRX but a unique chromosome translocation resulting in loss of DAXX function. While ATRX and DAXX form a complex in G292, this complex fails to localize to nuclear PML bodies. We demonstrate that introduction of wild type DAXX suppresses the ALT phenotype and restores the localization of ATRX/DAXX to PML bodies. Using an inducible system, we show that ALT-associated PML bodies are disrupted rapidly following DAXX induction and that ALT is again restored following withdrawal of DAXX.

Interaction between the microbiome and TP53 in human lung cancer.

BACKGROUND: Lung cancer is the leading cancer diagnosis worldwide and the number one cause of cancer deaths. Exposure to cigarette smoke, the primary risk factor in lung cancer, reduces epithelial barrier integrity and increases susceptibility to infections. Herein, we hypothesize that somatic mutations together with cigarette smoke generate a dysbiotic microbiota that is associated with lung carcinogenesis. Using lung tissue from 33 controls and 143 cancer cases, we conduct 16S ribosomal RNA (rRNA) bacterial gene sequencing, with RNA-sequencing data from lung cancer cases in The Cancer Genome Atlas serving as the validation cohort. RESULTS: Overall, we demonstrate a lower alpha diversity in normal lung as compared to non-tumor adjacent or tumor tissue. In squamous cell carcinoma specifically, a separate group of taxa are identified, in which Acidovorax is enriched in smokers. Acidovorax temporans is identified within tumor sections by fluorescent in situ hybridization and confirmed by two separate 16S rRNA strategies. Further, these taxa, including Acidovorax, exhibit higher abundance among the subset of squamous cell carcinoma cases with TP53 mutations, an association not seen in adenocarcinomas. CONCLUSIONS: The results of this comprehensive study show both microbiome-gene and microbiome-exposure interactions in squamous cell carcinoma lung cancer tissue. Specifically, tumors harboring TP53 mutations, which can impair epithelial function, have a unique bacterial consortium that is higher in relative abundance in smoking-associated tumors of this type. Given the significant need for clinical diagnostic tools in lung cancer, this study may provide novel biomarkers for early detection.

Characterization of BK Polyomaviruses from Kidney Transplant Recipients Suggests a Role for APOBEC3 in Driving In-Host Virus Evolution.

BK polyomavirus (BKV) frequently causes nephropathy (BKVN) in kidney transplant recipients (KTRs). BKV has also been implicated in the etiology of bladder and kidney cancers. We characterized BKV variants from two KTRs who developed BKVN followed by renal carcinoma. Both patients showed a swarm of BKV sequence variants encoding non-silent mutations in surface loops of the viral major capsid protein. The temporal appearance and disappearance of these mutations highlights the intra-patient evolution of BKV. Some of the observed mutations conferred resistance to antibody-mediated neutralization. The mutations also modified the spectrum of receptor glycans engaged by BKV during host cell entry. Intriguingly, all observed mutations were consistent with DNA damage caused by antiviral APOBEC3 cytosine deaminases. Moreover, APOBEC3 expression was evident upon immunohistochemical analysis of renal biopsies from KTRs. These results provide a snapshot of in-host BKV evolution and suggest that APOBEC3 may drive BKV mutagenesis in vivo.

Somatic mutations in murine models of leukemia and lymphoma: Disease specificity and clinical relevance.

Malignant transformation is a multistep process that is dictated by the acquisition of multiple genomic aberrations that provide growth and survival advantage. During the post genomic era, high throughput genomic sequencing has advanced exponentially, leading to identification of countless cancer associated mutations with potential for targeted therapy. Mouse models of cancer serve as excellent tools to examine the functionality of gene mutations and their contribution to the malignant process. However, it remains unclear whether the genetic events that occur during transformation are similar in mice and humans. To address that, we chose several transgenic mouse models of hematopoietic malignancies and identified acquired mutations in these mice by means of targeted re-sequencing of known cancer-associated genes as well as whole exome sequencing. We found that mutations that are typically found in acute myeloid leukemia or T cell acute lymphoblastic leukemia patients are also common in mouse models of the respective disease. Moreover, we found that the most frequent mutations found in a mouse model of lymphoma occur in a set of epigenetic modifier genes, implicating this pathway in the generation of lymphoma. These results demonstrate that genetically engineered mouse models (GEMM) mimic the genetic evolution of human cancer and serve as excellent platforms for target discovery and validation.

Progenitor B-1 B-cell acute lymphoblastic leukemia is associated with collaborative mutations in 3 critical pathways.

B-1 and B-2 lymphocytes are derived from distinct developmental pathways and represent layered arms of the innate and adaptive immune systems, respectively. In contrast to a majority of murine B-cell malignancies, which stain positive with the B220 antibody, we discovered a novel form of B-cell leukemia in NUP98-PHF23 (NP23) transgenic mice. The immunophenotype (Lin- B220- CD19+ AA4.1+) was identical to that of progenitor (pro) B-1 cells, and VH gene usage was skewed toward 3' V regions, similar to murine fetal liver B cells. Moreover, the gene expression profile of these leukemias was most similar to that of fetal liver pro-B fraction BC, a known source of B-1 B cells, further supporting a pro-B-1 origin of these leukemias. The NP23 pro-B-1 acute lymphoblastic leukemias (ALLs) acquired spontaneous mutations in both Bcor and Janus kinase (Jak) pathway (Jak1/2/3 and Stat5a) genes, supporting a hypothesis that mutations in 3 critical pathways (stem-cell self-renewal, B-cell differentiation, and cytokine signaling) collaborate to induce B-cell precursor (BCP) ALL. Finally, the thymic stromal lymphopoietin (Tslp) cytokine is required for murine B-1 development, and chromosomal rearrangements resulting in overexpression of the TSLP receptor (CRLF2) are present in some patients with high-risk BCP-ALL (referred to as CRLF2r ALL). Gene expression profiles of NP23 pro-B-1 ALL were more similar to that of CRLF2r ALL than non-CRLF2r ALL, and analysis of VH gene usage from patients with CRLF2r ALL demonstrated preferential usage of VH regions used by human B-1 B cells, leading to the suggestion that this subset of patients with BCP-ALL has a malignancy of B-1, rather than B-2, B-cell origin.

Point Mutations in Exon 1B of APC Reveal Gastric Adenocarcinoma and Proximal Polyposis of the Stomach as a Familial Adenomatous Polyposis Variant.

Gastric adenocarcinoma and proximal polyposis of the stomach (GAPPS) is an autosomal-dominant cancer-predisposition syndrome with a significant risk of gastric, but not colorectal, adenocarcinoma. We mapped the gene to 5q22 and found loss of the wild-type allele on 5q in fundic gland polyps from affected individuals. Whole-exome and -genome sequencing failed to find causal mutations but, through Sanger sequencing, we identified point mutations in APC promoter 1B that co-segregated with disease in all six families. The mutations reduced binding of the YY1 transcription factor and impaired activity of the APC promoter 1B in luciferase assays. Analysis of blood and saliva from carriers showed allelic imbalance of APC, suggesting that these mutations lead to decreased allele-specific expression in vivo. Similar mutations in APC promoter 1B occur in rare families with familial adenomatous polyposis (FAP). Promoter 1A is methylated in GAPPS and sporadic FGPs and in normal stomach, which suggests that 1B transcripts are more important than 1A in gastric mucosa. This might explain why all known GAPPS-affected families carry promoter 1B point mutations but only rare FAP-affected families carry similar mutations, the colonic cells usually being protected by the expression of the 1A isoform. Gastric polyposis and cancer have been previously described in some FAP-affected individuals with large deletions around promoter 1B. Our finding that GAPPS is caused by point mutations in the same promoter suggests that families with mutations affecting the promoter 1B are at risk of gastric adenocarcinoma, regardless of whether or not colorectal polyps are present.

Whole Genome Sequencing of Newly Established Pancreatic Cancer Lines Identifies Novel Somatic Mutation (c.2587G>A) in Axon Guidance Receptor Plexin A1 as Enhancer of Proliferation and Invasion.

The genetic profile of human pancreatic cancers harbors considerable heterogeneity, which suggests a possible explanation for the pronounced inefficacy of single therapies in this disease. This observation has led to a belief that custom therapies based on individual tumor profiles are necessary to more effectively treat pancreatic cancer. It has recently been discovered that axon guidance genes are affected by somatic structural variants in up to 25% of human pancreatic cancers. Thus far, however, some of these mutations have only been correlated to survival probability and no function has been assigned to these observed axon guidance gene mutations in pancreatic cancer. In this study we established three novel pancreatic cancer cell lines and performed whole genome sequencing to discover novel mutations in axon guidance genes that may contribute to the cancer phenotype of these cells. We discovered, among other novel somatic variants in axon guidance pathway genes, a novel mutation in the PLXNA1 receptor (c.2587G>A) in newly established cell line SB.06 that mediates oncogenic cues of increased invasion and proliferation in SB.06 cells and increased invasion in 293T cells upon stimulation with the receptor's natural ligand semaphorin 3A compared to wild type PLXNA1 cells. Mutant PLXNA1 signaling was associated with increased Rho-GTPase and p42/p44 MAPK signaling activity and cytoskeletal expansion, but not changes in E-cadherin, vimentin, or metalloproteinase 9 expression levels. Pharmacologic inhibition of the Rho-GTPase family member CDC42 selectively abrogated PLXNA1 c.2587G>A-mediated increased invasion. These findings provide in-vitro confirmation that somatic mutations in axon guidance genes can provide oncogenic gain-of-function signals and may contribute to pancreatic cancer progression.

The CDX1-microRNA-215 axis regulates colorectal cancer stem cell differentiation.

The transcription factor caudal-type homeobox 1 (CDX1) is a key regulator of differentiation in the normal colon and in colorectal cancer (CRC). CDX1 activates the expression of enterocyte genes, but it is not clear how the concomitant silencing of stem cell genes is achieved. MicroRNAs (miRNAs) are important mediators of gene repression and have been implicated in tumor suppression and carcinogenesis, but the roles of miRNAs in differentiation, particularly in CRC, remain poorly understood. Here, we identified microRNA-215 (miR-215) as a direct transcriptional target of CDX1 by using high-throughput small RNA sequencing to profile miRNA expression in two pairs of CRC cell lines: CDX1-low HCT116 and HCT116 with stable CDX1 overexpression, and CDX1-high LS174T and LS174T with stable CDX1 knockdown. Validation of candidate miRNAs identified by RNA-seq in a larger cell-line panel revealed miR-215 to be most significantly correlated with CDX1 expression. Quantitative ChIP-PCR and promoter luciferase assays confirmed that CDX1 directly activates miR-215 transcription. miR-215 expression is depleted in FACS-enriched cancer stem cells compared with unsorted samples. Overexpression of miR-215 in poorly differentiated cell lines causes a decrease in clonogenicity, whereas miR-215 knockdown increases clonogenicity and impairs differentiation in CDX1-high cell lines. We identified the genome-wide targets of miR-215 and found that miR-215 mediates the repression of cell cycle and stemness genes downstream of CDX1. In particular, the miR-215 target gene BMI1 has been shown to promote stemness and self-renewal and to vary inversely with CDX1. Our work situates miR-215 as a link between CDX1 expression and BMI1 repression that governs differentiation in CRC.

Recurrent epimutation of SDHC in gastrointestinal stromal tumors.

Succinate dehydrogenase (SDH) is a conserved effector of cellular metabolism and energy production, and loss of SDH function is a driver mechanism in several cancers. SDH-deficient gastrointestinal stromal tumors (dSDH GISTs) collectively manifest similar phenotypes, including hypermethylated epigenomic signatures, tendency to occur in pediatric patients, and lack of KIT/PDGFRA mutations. dSDH GISTs often harbor deleterious mutations in SDH subunit genes (SDHA, SDHB, SDHC, and SDHD, termed SDHx), but some are SDHx wild type (WT). To further elucidate mechanisms of SDH deactivation in SDHx-WT GIST, we performed targeted exome sequencing on 59 dSDH GISTs to identify 43 SDHx-mutant and 16 SDHx-WT cases. Genome-wide DNA methylation and expression profiling exposed SDHC promoter-specific CpG island hypermethylation and gene silencing in SDHx-WT dSDH GISTs [15 of 16 cases (94%)]. Six of 15 SDHC-epimutant GISTs occurred in the setting of the multitumor syndrome Carney triad. We observed neither SDHB promoter hypermethylation nor large deletions on chromosome 1q in any SDHx-WT cases. Deep genome sequencing of a 130-kbp (kilo-base pair) window around SDHC revealed no recognizable sequence anomalies in SDHC-epimutant tumors. More than 2000 benign and tumor reference tissues, including stem cells and malignancies with a hypermethylator epigenotype, exhibit solely a non-epimutant SDHC promoter. Mosaic constitutional SDHC promoter hypermethylation in blood and saliva from patients with SDHC-epimutant GIST implicates a postzygotic mechanism in the establishment and maintenance of SDHC epimutation. The discovery of SDHC epimutation provides a unifying explanation for the pathogenesis of dSDH GIST, whereby loss of SDH function stems from either SDHx mutation or SDHC epimutation.

A specific missense mutation in GTF2I occurs at high frequency in thymic epithelial tumors.

We analyzed 28 thymic epithelial tumors (TETs) using next-generation sequencing and identified a missense mutation (chromosome 7 c.74146970T>A) in GTF2I at high frequency in type A thymomas, a relatively indolent subtype. In a series of 274 TETs, we detected the GTF2I mutation in 82% of type A and 74% of type AB thymomas but rarely in the aggressive subtypes, where recurrent mutations of known cancer genes have been identified. Therefore, GTF2I mutation correlated with better survival. GTF2I ß and δ isoforms were expressed in TETs, and both mutant isoforms were able to stimulate cell proliferation in vitro. Thymic carcinomas carried a higher number of mutations than thymomas (average of 43.5 and 18.4, respectively). Notably, we identified recurrent mutations of known cancer genes, including TP53, CYLD, CDKN2A, BAP1 and PBRM1, in thymic carcinomas. These findings will complement the diagnostic assessment of these tumors and also facilitate development of a molecular classification and assessment of prognosis and treatment strategies.

Characterization of fibroblast growth factor receptor 1 in small-cell lung cancer.

INTRODUCTION: There remains a significant therapeutic need for small-cell lung cancer (SCLC). We and others have reported high frequency of copy number gains in cytogenetic bands encoding fibroblast growth factor receptor 1 (FGFR1) in SCLC tumors and cell lines. METHODS: Thirteen SCLC cell lines and 68 SCLC patient tumor samples were studied for FGFR1 amplification. Growth inhibition assays were performed using PD173074, a pan-FGFR inhibitor to determine the correlation between FGFR1 expression and drug sensitivity. RESULTS: We did not detect FGFR1 mutations in SCLC cell lines. Focal amplification of FGFR1 gene was found in five tumor samples (7%), with high-level focal amplification in only one tumor sample (1%). Amplification owing to polysomy of chromosome 8, where FGFR1 locates, was observed in 22 tumor samples (32%). There was no correlation between FGFR1 gene copy number and messenger RNA expression or protein expression in SCLC cells. FGFR inhibitor sensitivity correlated with FGFR1 copy number determined by real-time polymerase chain reaction assay (r= -0.79; p = 0.01). CONCLUSION: FGFR1 gene mutations and focal amplification are rare in SCLC, but polysomy of chromosome 8 is relatively common. FGFR1 copy number gain predicts sensitivity to FGFR inhibition, and FGFR expression correlates inversely with chemosensitivity.

NUP98-PHF23 is a chromatin-modifying oncoprotein that causes a wide array of leukemias sensitive to inhibition of PHD histone reader function.

In this report, we show that expression of a NUP98-PHF23 (NP23) fusion, associated with acute myeloid leukemia (AML) in humans, leads to myeloid, erythroid, T-cell, and B-cell leukemia in mice. The leukemic and preleukemic tissues display a stem cell-like expression signature, including Hoxa, Hoxb, and Meis1 genes. The PHF23 plant homeodomain (PHD) motif is known to bind to H3K4me3 residues, and chromatin immunoprecipitation experiments demonstrated that the NP23 protein binds to chromatin at a specific subset of H3K4me3 sites, including at Hoxa, Hoxb, and Meis1. Treatment of NP23 cells with disulfiram, which inhibits the binding of PHD motifs to H3K4me3, rapidly and selectively killed NP23-expressing myeloblasts; cell death was preceded by decreased expression of Hoxa, Hoxb, and Meis1. Furthermore, AML driven by a related fusion gene, NUP98-JARID1A (NJL), was also sensitive to disulfiram. Thus, the NP23 mouse provides a platform to evaluate compounds that disrupt binding of oncogenic PHD proteins to H3K4me3.

A p21-ZEB1 complex inhibits epithelial-mesenchymal transition through the microRNA 183-96-182 cluster.

The tumor suppressor p21 acts as a cell cycle inhibitor and has also been shown to regulate gene expression by functioning as a transcription corepressor. Here, we identified p21-regulated microRNAs (miRNAs) by sequencing small RNAs from isogenic p21(+/+) and p21(-/-) cells. Three abundant miRNA clusters, miR-200b-200a-429, miR-200c-141, and miR-183-96-182, were downregulated in p21-deficient cells. Consistent with the known function of the miR-200 family and p21 in inhibition of the epithelial-mesenchymal transition (EMT), we observed EMT upon loss of p21 in multiple model systems. To explore a role of the miR-183-96-182 cluster in EMT, we identified its genome-wide targets and found that miR-183 and miR-96 repressed common targets, including SLUG, ZEB1, ITGB1, and KLF4. Reintroduction of miR-200, miR-183, or miR-96 in p21(-/-) cells inhibited EMT, cell migration, and invasion. Conversely, antagonizing miR-200 and miR-183-96-182 cluster miRNAs in p21(+/+) cells increased invasion and elevated the levels of VIM, ZEB1, and SLUG mRNAs. Furthermore, we found that p21 forms a complex with ZEB1 at the miR-183-96-182 cluster promoter to inhibit transcriptional repression of this cluster by ZEB1, suggesting a reciprocal feedback loop.

High prevalence of MAP2K1 mutations in variant and IGHV4-34-expressing hairy-cell leukemias.

To understand the genetic mechanisms driving variant and IGHV4-34-expressing hairy-cell leukemias, we performed whole-exome sequencing of leukemia samples from ten affected individuals, including six with matched normal samples. We identified activating mutations in the MAP2K1 gene (encoding MEK1) in 5 of these 10 samples and in 10 of 21 samples in a validation set (overall frequency of 15/31), suggesting potential new strategies for treating individuals with these diseases.

Oncogenic ETS fusions deregulate E2F3 target genes in Ewing sarcoma and prostate cancer.

Deregulated E2F transcription factor activity occurs in the vast majority of human tumors and has been solidly implicated in disturbances of cell cycle control, proliferation, and apoptosis. Aberrant E2F regulatory activity is often caused by impairment of control through pRB function, but little is known about the interplay of other oncoproteins with E2F. Here we show that ETS transcription factor fusions resulting from disease driving rearrangements in Ewing sarcoma (ES) and prostate cancer (PC) are one such class of oncoproteins. We performed an integrative study of genome-wide DNA-binding and transcription data in EWSR1/FLI1 expressing ES and TMPRSS2/ERG containing PC cells. Supported by promoter activity and mutation analyses, we demonstrate that a large fraction of E2F3 target genes are synergistically coregulated by these aberrant ETS proteins. We propose that the oncogenic effect of ETS fusion oncoproteins is in part mediated by the disruptive effect of the E2F-ETS interaction on cell cycle control. Additionally, a detailed analysis of the regulatory targets of the characteristic EWSR1/FLI1 fusion in ES identifies two functionally distinct gene sets. While synergistic regulation in concert with E2F in the promoter of target genes has a generally activating effect, EWSR1/FLI1 binding independent of E2F3 is predominantly associated with repressed differentiation genes. Thus, EWSR1/FLI1 appears to promote oncogenesis by simultaneously promoting cell proliferation and perturbing differentiation.

The exomes of the NCI-60 panel: a genomic resource for cancer biology and systems pharmacology.

The NCI-60 cell lines are the most frequently studied human tumor cell lines in cancer research. This panel has generated the most extensive cancer pharmacology database worldwide. In addition, these cell lines have been intensely investigated, providing a unique platform for hypothesis-driven research focused on enhancing our understanding of tumor biology. Here, we report a comprehensive analysis of coding variants in the NCI-60 panel of cell lines identified by whole exome sequencing, providing a list of possible cancer specific variants for the community. Furthermore, we identify pharmacogenomic correlations between specific variants in genes such as TP53, BRAF, ERBBs, and ATAD5 and anticancer agents such as nutlin, vemurafenib, erlotinib, and bleomycin showing one of many ways the data could be used to validate and generate novel hypotheses for further investigation. As new cancer genes are identified through large-scale sequencing studies, the data presented here for the NCI-60 will be an invaluable resource for identifying cell lines with mutations in such genes for hypothesis-driven research. To enhance the utility of the data for the greater research community, the genomic variants are freely available in different formats and from multiple sources including the CellMiner and Ingenuity websites.

Massively-parallel sequencing of genes on a single chromosome: a comparison of solution hybrid selection and flow sorting.

BACKGROUND: Targeted capture, combined with massively-parallel sequencing, is a powerful technique that allows investigation of specific portions of the genome for less cost than whole genome sequencing. Several methods have been developed, and improvements have resulted in commercial products targeting the human or mouse exonic regions (the exome). In some cases it is desirable to custom-target other regions of the genome, either to reduce the amount of sequence that is targeted or to capture regions that are not targeted by commercial kits. It is important to understand the advantages, limitations, and complexity of a given capture method before embarking on a targeted sequencing experiment. RESULTS: We compared two custom targeted capture methods suitable for single chromosome analysis: Solution Hybrid Selection (SHS) and Flow Sorting (FS) of single chromosomes. Both methods can capture targeted material and result in high percentages of genotype identifications across these regions: 59-92% for SHS and 70-79% for FS. FS is amenable to current structural variation detection methods, and variants were detected. Structural variation was also assessed for SHS samples with paired end sequencing, resulting in variant identification. CONCLUSIONS: While both methods can effectively target genomic regions for genotype determination, several considerations make each method appropriate in different circumstances. SHS is well suited for experiments targeting smaller regions in a larger number of samples. FS is well suited when regions of interest cover large regions of a single chromosome. Although whole genome sequencing is becoming less expensive, the sequencing, data storage, and analysis costs make targeted sequencing using SHS or FS a compelling option.