Integrative genomic analysis reveals cancer-associated mutations at diagnosis of CML in patients with high-risk disease.

Genomic events associated with poor outcome in chronic myeloid leukemia (CML) are poorly understood. We performed whole-exome sequencing, copy-number variation, and/or RNA sequencing for 65 patients to discover mutations at diagnosis and blast crisis (BC). Forty-six patients with chronic-phase disease with the extremes of outcome were studied at diagnosis. Cancer gene variants were detected in 15 (56%) of 27 patients with subsequent BC or poor outcome and in 3 (16%) of 19 optimal responders (P = .007). Frequently mutated genes at diagnosis were ASXL1, IKZF1, and RUNX1 The methyltransferase SETD1B was a novel recurrently mutated gene. A novel class of variant associated with the Philadelphia (Ph) translocation was detected at diagnosis in 11 (24%) of 46 patients comprising fusions and/or rearrangement of genes on the translocated chromosomes, with evidence of fragmentation, inversion, and imperfect sequence reassembly. These were more frequent at diagnosis in patients with poor outcome: 9 (33%) of 27 vs 2 (11%) of 19 optimal responders (P = .07). Thirty-nine patients were tested at BC, and all had cancer gene variants, including ABL1 kinase domain mutations in 58%. However, ABL1 mutations cooccurred with other mutated cancer genes in 89% of cases, and these predated ABL1 mutations in 62% of evaluable patients. Gene fusions not associated with the Ph translocation occurred in 42% of patients at BC and commonly involved fusion partners that were known cancer genes (78%). Genomic analysis revealed numerous relevant variants at diagnosis in patients with poor outcome and all patients at BC. Future refined biomarker testing of specific variants will likely provide prognostic information to facilitate a risk-adapted therapeutic approach.

HENMT1 and piRNA Stability Are Required for Adult Male Germ Cell Transposon Repression and to Define the Spermatogenic Program in the Mouse.

piRNAs are critical for transposable element (TE) repression and germ cell survival during the early phases of spermatogenesis, however, their role in adult germ cells and the relative importance of piRNA methylation is poorly defined in mammals. Using a mouse model of HEN methyltransferase 1 (HENMT1) loss-of-function, RNA-Seq and a range of RNA assays we show that HENMT1 is required for the 2' O-methylation of mammalian piRNAs. HENMT1 loss leads to piRNA instability, reduced piRNA bulk and length, and ultimately male sterility characterized by a germ cell arrest at the elongating germ cell phase of spermatogenesis. HENMT1 loss-of-function, and the concomitant loss of piRNAs, resulted in TE de-repression in adult meiotic and haploid germ cells, and the precocious, and selective, expression of many haploid-transcripts in meiotic cells. Precocious expression was associated with a more active chromatin state in meiotic cells, elevated levels of DNA damage and a catastrophic deregulation of the haploid germ cell gene expression. Collectively these results define a critical role for HENMT1 and piRNAs in the maintenance of TE repression in adult germ cells and setting the spermatogenic program.

Mutations in the mammalian target of rapamycin pathway regulators NPRL2 and NPRL3 cause focal epilepsy.

OBJECTIVE: Focal epilepsies are the most common form observed and have not generally been considered to be genetic in origin. Recently, we identified mutations in DEPDC5 as a cause of familial focal epilepsy. In this study, we investigated whether mutations in the mammalian target of rapamycin (mTOR) regulators, NPRL2 and NPRL3, also contribute to cases of focal epilepsy. METHODS: We used targeted capture and next-generation sequencing to analyze 404 unrelated probands with focal epilepsy. We performed exome sequencing on two families with multiple members affected with focal epilepsy and linkage analysis on one of these. RESULTS: In our cohort of 404 unrelated focal epilepsy patients, we identified five mutations in NPRL2 and five in NPRL3. Exome sequencing analysis of two families with focal epilepsy identified NPRL2 and NPRL3 as the top candidate-causative genes. Some patients had focal epilepsy associated with brain malformations. We also identified 18 new mutations in DEPDC5. INTERPRETATION: We have identified NPRL2 and NPRL3 as two new focal epilepsy genes that also play a role in the mTOR-signaling pathway. Our findings show that mutations in GATOR1 complex genes are the most significant cause of familial focal epilepsy identified to date, including cases with brain malformations. It is possible that deregulation of cellular growth control plays a more important role in epilepsy than is currently recognized.

Discovery of a previously unrecognized microdeletion syndrome of 16p11.2-p12.2.

We have identified a recurrent de novo pericentromeric deletion in 16p11.2-p12.2 in four individuals with developmental disabilities by microarray-based comparative genomic hybridization analysis. The identification of common clinical features in these four individuals along with the characterization of complex segmental duplications flanking the deletion regions suggests that nonallelic homologous recombination mediated these rearrangements and that deletions in 16p11.2-p12.2 constitute a previously undescribed syndrome.

Detection of clinically relevant copy number variants with whole-exome sequencing.

Copy number variation (CNV) is a common source of genetic variation that has been implicated in many genomic disorders. This has resulted in the widespread application of genomic microarrays as a first-tier diagnostic tool for CNV detection. More recently, whole-exome sequencing (WES) has been proven successful for the detection of clinically relevant point mutations and small insertion-deletions exome wide. We evaluate the utility of short-read WES (SOLiD 5500xl) to detect clinically relevant CNVs in DNA from 10 patients with intellectual disability and compare these results to data from two independent high-resolution microarrays. Eleven of the 12 clinically relevant CNVs were detected via read-depth analysis of WES data; a heterozygous single-exon deletion remained undetected by all algorithms evaluated. Although the detection power of WES for small CNVs currently does not match that of high-resolution microarray platforms, we show that the majority (88%) of rare coding CNVs containing three or more exons are successfully identified by WES. These results show that the CNV detection resolution of WES is comparable to that of medium-resolution genomic microarrays commonly used as clinical assays. The combined detection of point mutations, indels, and CNVs makes WES a very attractive first-tier diagnostic test for genetically heterogeneous disorders.

Acquired copy number alterations in adult acute myeloid leukemia genomes.

Cytogenetic analysis of acute myeloid leukemia (AML) cells has accelerated the identification of genes important for AML pathogenesis. To complement cytogenetic studies and to identify genes altered in AML genomes, we performed genome-wide copy number analysis with paired normal and tumor DNA obtained from 86 adult patients with de novo AML using 1.85 million feature SNP arrays. Acquired copy number alterations (CNAs) were confirmed using an ultra-dense array comparative genomic hybridization platform. A total of 201 somatic CNAs were found in the 86 AML genomes (mean, 2.34 CNAs per genome), with French-American-British system M6 and M7 genomes containing the most changes (10-29 CNAs per genome). Twenty-four percent of AML patients with normal cytogenetics had CNA, whereas 40% of patients with an abnormal karyotype had additional CNA detected by SNP array, and several CNA regions were recurrent. The mRNA expression levels of 57 genes were significantly altered in 27 of 50 recurrent CNA regions <5 megabases in size. A total of 8 uniparental disomy (UPD) segments were identified in the 86 genomes; 6 of 8 UPD calls occurred in samples with a normal karyotype. Collectively, 34 of 86 AML genomes (40%) contained alterations not found with cytogenetics, and 98% of these regions contained genes. Of 86 genomes, 43 (50%) had no CNA or UPD at this level of resolution. In this study of 86 adult AML genomes, the use of an unbiased high-resolution genomic screen identified many genes not previously implicated in AML that may be relevant for pathogenesis, along with many known oncogenes and tumor suppressor genes.

Towards Routine Implementation of Liquid Biopsies in Cancer Management: It Is Always Too Early, until Suddenly It Is Too Late.

Blood-based liquid biopsies are considered a new and promising diagnostic and monitoring tool for cancer. As liquid biopsies only require a blood draw, they are non-invasive, potentially more rapid and assumed to be a less costly alternative to genomic analysis of tissue biopsies. A multi-disciplinary workshop (n = 98 registrations) was organized to discuss routine implementation of liquid biopsies in cancer management. Real-time polls were used to engage with experts' about the current evidence of clinical utility and the barriers to implementation of liquid biopsies. Clinical, laboratory and health economics presentations were given to illustrate the opportunities and current levels of evidence, followed by three moderated break-out sessions to discuss applications. The workshop concluded that tumor-informed assays using next-generation sequencing (NGS) or PCR-based genotyping assays will most likely provide better clinical utility than tumor-agnostic assays, yet at a higher cost. For routine application, it will be essential to determine clinical utility, to define the minimum quality standards and performance of testing platforms and to ensure their use is integrated into current clinical workflows including how they complement tissue biopsies and imaging. Early health economic models may help identifying the most viable application of liquid biopsies. Alternative funding models for the translation of complex molecular diagnostics, such as liquid biopsies, may also be explored if clinical utility has been demonstrated and when their use is recommended in multi-disciplinary consensus guidelines.

Aberrant Splicing of SDHC in Families With Unexplained Succinate Dehydrogenase-Deficient Paragangliomas.

CONTEXT: Germline mutations in the succinate dehydrogenase genes (SDHA/B/C/D, SDHAF2-collectively, "SDHx") have been implicated in paraganglioma (PGL), renal cell carcinoma (RCC), gastrointestinal stromal tumor (GIST), and pituitary adenoma (PA). Negative SDHB tumor staining is indicative of SDH-deficient tumors, usually reflecting an underlying germline SDHx mutation. However, approximately 20% of individuals with SDH-deficient tumors lack an identifiable germline SDHx mutation. METHODS: We performed whole-exome sequencing (WES) of germline and tumor DNA followed by Sanger sequencing validation, transcriptome analysis, metabolomic studies, and haplotype analysis in 2 Italian-Australian families with SDH-deficient PGLs and various neoplasms, including RCC, GIST, and PA. RESULTS: Germline WES revealed a novel SDHC intronic variant, which had been missed during previous routine testing, in 4 affected siblings of the index family. Transcriptome analysis demonstrated aberrant SDHC splicing, with the retained intronic segment introducing a premature stop codon. WES of available tumors in this family showed chromosome 1 deletion with loss of wild-type SDHC in a PGL and a somatic gain-of-function KIT mutation in a GIST. The SDHC intronic variant identified was subsequently detected in the second family, with haplotype analysis indicating a founder effect. CONCLUSIONS: This is the deepest intronic variant to be reported among the SDHx genes. Intronic variants beyond the limits of standard gene sequencing analysis should be considered in patients with SDH-deficient tumors but negative genetic test results.

The mutational burden of therapy-related myeloid neoplasms is similar to primary myelodysplastic syndrome but has a distinctive distribution.

Therapy-related myeloid neoplasms (T-MN) are poorly characterized secondary hematological malignancies following chemotherapy/radiotherapy exposure. We compared the clinical and mutational characteristics of T-MN (n = 129) and primary myelodysplastic syndrome (P-MDS, n = 108) patients. Although the somatic mutation frequency was similar between T-MN and P-MDS patients (93% in both groups), the pattern was distinct. TP53 mutations were more frequent in T-MN (29.5 vs. 7%), while spliceosomal complex mutations were more common in P-MDS (56.5 vs. 25.6%). In contrast to P-MDS, the ring sideroblasts (RS) phenotype was not associated with better survival in T-MN, most probably due to genetic association with TP53 mutations. SF3B1 was mutated in 96% of P-MDS with ≥15% RS, but in only 32% T-MN. TP53 mutations were detected in 92% T-MN with ≥15% RS and SF3B1 wild-type cases. Interestingly, T-MN and P-MDS patients with "Very low" or "Low" Revised International Prognostic Scoring System (IPSS-R) showed similar biological and clinical characteristics. In a Cox regression analysis, TP53 mutation was a poor prognostic factor in T-MN, independent of IPSS-R cytogenetics, disease-modifying therapy, and NRAS mutation. Our data have direct implications for T-MN management and provide evidence that, in addition to conventional disease parameters, mutational analysis should be incorporated in T-MN risk stratification.

DNA Barcoding Reveals Habitual Clonal Dominance of Myeloma Plasma Cells in the Bone Marrow Microenvironment.

Multiple myeloma (MM) is a hematological malignancy resulting from the uncontrolled proliferation of antibody-producing plasma cells in the bone marrow. At diagnosis, independent plasma cell tumors are found throughout the skeleton. The recirculation of mutant plasma cells from the initial lesion and their recolonization of distant marrow sites are thought to occur by a process similar to solid tumor metastasis. However, the efficiency of this bone marrow homing process and the proportion of disseminated cells that actively divide and contribute to new tumor growth in MM are both unknown. We used the C57BL/KaLwRij mouse model of myeloma, lentiviral-mediated DNA barcoding of 5TGM1 myeloma cells, and next-generation sequencing to investigate the relative efficiency of plasma cell migration to, and growth within, the bone marrow. This approach revealed three major findings: firstly, establishment of metastasis within the bone marrow was extremely inefficient, with approximately 0.01% of circulating myeloma cells becoming resident long term in the bone marrow of each long bone; secondly, the individual cells of each metastasis exhibited marked differences in their proliferative fates, with the majority of final tumor burden within a bone being attributable to the progeny of between 1 and 8 cells; and, thirdly, the proliferative fate of individual clonal plasma cells differed at each bone marrow site in which the cells "landed." These findings suggest that individual myeloma plasma cells are subjected to vastly different selection pressures within the bone marrow microenvironment, highlighting the importance of niche-driven factors, which determine the disease course and outcome.

Metastatic squamous cell carcinoma cells that overexpress c-Met exhibit enhanced angiogenesis factor expression, scattering and metastasis in response to hepatocyte growth factor.

We previously performed gene expression profiling in a multistep squamous cell carcinoma (SCC) progression model, and identified growth-regulated oncogene-1 (Gro-1/KC) as a factor that contributes to enhanced angiogenesis, tumorigenesis and metastasis. In the present study, we explored molecular pathways coactivated with Gro-1/KC, and identified a transcript that encodes c-Met, the receptor for hepatocyte growth factor/scatter factor (HGF). Northern, Western blot, and immunohistochemical analyses confirm that expression of c-Met mRNA and protein is increased with SCC progression. In vitro, HGF preferentially promoted scattering in the metastatic LY-1 and LY-2 lines, and enhanced angiogenesis factors Gro-1/KC and vascular endothelial growth factor (VEGF) production by all tumor cell lines. In vivo, tumor growth and lung metastasis were promoted by transfection and overexpression of HGF cDNA in metastatic LY-1 cells. Our data indicate that metastatic SCC cells that overexpress c-Met exhibit angiogenesis factor expression and enhanced scattering in response to HGF in vitro, and tumorigenesis and metastasis in response to HGF in the tumor microenvironment in vivo.

Three microarray platforms: an analysis of their concordance in profiling gene expression.

BACKGROUND: Microarrays for the analysis of gene expression are of three different types: short oligonucleotide (25-30 base), long oligonucleotide (50-80 base), and cDNA (highly variable in length). The short oligonucleotide and cDNA arrays have been the mainstay of expression analysis to date, but long oligonucleotide platforms are gaining in popularity and will probably replace cDNA arrays. As part of a validation study for the long oligonucleotide arrays, we compared and contrasted expression profiles from the three formats, testing RNA from six different cell lines against a universal reference standard. RESULTS: The three platforms had 6430 genes in common. In general, correlation of gene expression levels across the platforms was good when defined by concordance in the direction of expression difference (upregulation or downregulation), scatter plot analysis, principal component analysis, cell line correlation or quantitative RT-PCR. The overall correlations (r values) between platforms were in the range 0.7 to 0.8, as determined by analysis of scatter plots. When concordance was measured for expression ratios significant at p-values of <0.05 and at expression threshold levels of 1.5 and 2-fold, the agreement among the platforms was very high, ranging from 93% to 100%. CONCLUSION: Our results indicate that the long oligonucleotide platform is highly suitable for expression analysis and compares favorably with the cDNA and short oligonucleotide varieties. All three platforms can give similar and reproducible results if the criterion is the direction of change in gene expression and minimal emphasis is placed on the magnitude of change.

Multiple-laboratory comparison of microarray platforms.

Microarray technology is a powerful tool for measuring RNA expression for thousands of genes at once. Various studies have been published comparing competing platforms with mixed results: some find agreement, others do not. As the number of researchers starting to use microarrays and the number of cross-platform meta-analysis studies rapidly increases, appropriate platform assessments become more important. Here we present results from a comparison study that offers important improvements over those previously described in the literature. In particular, we noticed that none of the previously published papers consider differences between labs. For this study, a consortium of ten laboratories from the Washington, DC-Baltimore, USA, area was formed to compare data obtained from three widely used platforms using identical RNA samples. We used appropriate statistical analysis to demonstrate that there are relatively large differences in data obtained in labs using the same platform, but that the results from the best-performing labs agree rather well.