Deep sequencing of multiple regions of glial tumors reveals spatial heterogeneity for mutations in clinically relevant genes

Background

The extent of intratumoral mutational heterogeneity remains unclear in gliomas, the most common primary brain tumors, especially with respect to point mutation. To address this, we applied single molecule molecular inversion probes targeting 33 cancer genes to assay both point mutations and gene amplifications within spatially distinct regions of 14 glial tumors.

Results

We find evidence of regional mutational heterogeneity in multiple tumors, including mutations in TP53 and RB1 in an anaplastic oligodendroglioma and amplifications in PDGFRA and KIT in two glioblastomas (GBMs). Immunohistochemistry confirms heterogeneity of TP53 mutation and PDGFRA amplification. In all, 3 out of 14 glial tumors surveyed have evidence for heterogeneity for clinically relevant mutations.

Conclusions

Our results underscore the need to sample multiple regions in GBM and other glial tumors when devising personalized treatments based on genomic information, and furthermore demonstrate the importance of measuring both point mutation and copy number alteration while investigating genetic heterogeneity within cancer samples.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0530-z) contains supplementary material, which is available to authorized users.     

Reporting Tumor Molecular Heterogeneity in Histopathological Diagnosis

Background

Detection of molecular tumor heterogeneity has become of paramount importance with the advent of targeted therapies. Analysis for detection should be comprehensive, timely and based on routinely available tumor samples.

Aim

To evaluate the diagnostic potential of targeted multigene next-generation sequencing (TM-NGS) in characterizing gastrointestinal cancer molecular heterogeneity.

Methods

35 gastrointestinal tract tumors, five of each intestinal type gastric carcinomas, pancreatic ductal adenocarcinomas, pancreatic intraductal papillary mucinous neoplasms, ampulla of Vater carcinomas, hepatocellular carcinomas, cholangiocarcinomas, pancreatic solid pseudopapillary tumors were assessed for mutations in 46 cancer-associated genes, using Ion Torrent semiconductor-based TM-NGS. One ampulla of Vater carcinoma cell line and one hepatic carcinosarcoma served to assess assay sensitivity. TP53, PIK3CA, KRAS, and BRAF mutations were validated by conventional Sanger sequencing.

Results

TM-NGS yielded overlapping results on matched fresh-frozen and formalin-fixed paraffin-embedded (FFPE) tissues, with a mutation detection limit of 1% for fresh-frozen high molecular weight DNA and 2% for FFPE partially degraded DNA. At least one somatic mutation was observed in all tumors tested; multiple alterations were detected in 20/35 (57%) tumors. Seven cancers displayed significant differences in allelic frequencies for distinct mutations, indicating the presence of intratumor molecular heterogeneity; this was confirmed on selected samples by immunohistochemistry of p53 and Smad4, showing concordance with mutational analysis.

Conclusions

TM-NGS is able to detect and quantitate multiple gene alterations from limited amounts of DNA, moving one step closer to a next-generation histopathologic diagnosis that integrates morphologic, immunophenotypic, and multigene mutational analysis on routinely processed tissues, essential for personalized cancer therapy.     

Novel HRPT2/CDC73 Gene Mutations and Loss of Expression of Parafibromin in Chinese Patients with Clinically Sporadic Parathyroid Carcinomas

Objective

It is widely recognized that the diagnosis of parathyroid carcinoma (PC) is often difficult because of the overlap of characteristics between malignant and benign parathyroid tumors, especially at an early stage. Based on the identification of tumor suppressor gene HRPT2/CDC73 and its association with hereditary and sporadic PC, screening of gene mutations and detection of parafibromin immunoreactivity have been suggested as diagnostic instruments of PC in Whites. There is little information about HRPT2/CDC73 mutations and its corresponding protein expression in patients with sporadic PC in Chinese population, and the long-term follow-up data is scarce.

Methods

Paraffin-embedded tissues were obtained from 13 patients with PC, 13 patients with parathyroid adenoma (PA) and 7 patients with parathyroid hyperplasia(PH), and 6 normal parathyroid (NP) tissues as controls. Peripheral blood from 11 patients with PC was collected. PCR products using Genomic DNA extracted from tumor tissues or blood as template was sequenced for HRPT2/CDC73 gene. Expression of parafibromin in tumor tissues was evaluated by immunohistochemical analysis.

Results

Six mutations in 6 of 13 patients with PC were identified, with three being novel. Four of them were germ-line mutations. Patients with mutations were susceptible to recurrence of the PC. Complete (8/13, 61.5%) or partial (5/13, 38.5%) loss of parafibromin expression was observed in PC tissues. All of tissue samples from normal parathyroid or benign parathyroid tumors displayed positive immunostaining of parafibromin except one adenoma.

Conclusions

The present study supplies information on the mutations and protein expression of HRPT2/CDC73 gene and phenotypes of parathyroid carcinoma in Chinese population. And the expanded mutation database of this gene may benefit patients in the diagnosis and treatment of this disease.     

The mutational landscape of chromatin regulatory factors across 4,623 tumor samples

Background

Chromatin regulatory factors are emerging as important genes in cancer development and are regarded as interesting candidates for novel targets for cancer treatment. However, we lack a comprehensive understanding of the role of this group of genes in different cancer types.

Results

We have analyzed 4,623 tumor samples from thirteen anatomical sites to determine which chromatin regulatory factors are candidate drivers in these different sites. We identify 34 chromatin regulatory factors that are likely drivers in tumors from at least one site, all with relatively low mutational frequency. We also analyze the relative importance of mutations in this group of genes for the development of tumorigenesis in each site, and in different tumor types from the same site.

Conclusions

We find that, although tumors from all thirteen sites show mutations in likely driver chromatin regulatory factors, these are more prevalent in tumors arising from certain tissues. With the exception of hematopoietic, liver and kidney tumors, as a median, the mutated factors are less than one fifth of all mutated drivers across all sites analyzed. We also show that mutations in two of these genes, MLL and EP300, correlate with broad expression changes across cancer cell lines, thus presenting at least one mechanism through which these mutations could contribute to tumorigenesis in cells of the corresponding tissues.     

Mutational Profile of GNAQ Q209 in Human Tumors

Background

Frequent somatic mutations have recently been identified in the ras-like domain of the heterotrimeric G protein α-subunit (GNAQ) in blue naevi 83%, malignant blue naevi (50%) and ocular melanoma of the uvea (46%). The mutations exclusively affect codon 209 and result in GNAQ constitutive activation which, in turn, acts as a dominant oncogene.

Methodology

To assess if the mutations are present in other tumor types we performed a systematic mutational profile of the GNAQ exon 5 in a panel of 922 neoplasms, including glioblastoma, gastrointestinal stromal tumors (GIST), acute myeloid leukemia (AML), blue naevi, skin melanoma, bladder, breast, colorectal, lung, ovarian, pancreas, and thyroid carcinomas.

Principal Findings

We detected the previously reported mutations in 6/13 (46%) blue naevi. Changes affecting Q209 were not found in any of the other tumors. Our data indicate that the occurrence of GNAQ mutations display a unique pattern being present in a subset of melanocytic tumors but not in malignancies of glial, epithelial and stromal origin analyzed in this study.     

Mitochondrial genome variations in advanced stage endometriosis: a study in South Indian population

Background

Endometriosis is a chronic gynecological benign disease that shares several features similar to malignancy. Mitochondrial DNA (mtDNA) mutations have been reported in all most all types of tumors. However, it is not known as to whether mtDNA mutations are associated with endometriosis.

Methodology

We sequenced the entire mitochondrial genome of analogous ectopic and eutopic endometrial tissues along with blood samples from 32 advanced stage endometriosis patients to analyze the role of somatic and germ-line mtDNA variations in pathogenesis of endometriosis. All ectopic tissues were screened for tumor-specific mtDNA deletions and microsatellite instability (MSI). We also performed mtDNA haplogrouping in 128 patients and 90 controls to identify its possible association with endometriosis risk.

Principal Findings

We identified 51 somatic (novel: 31; reported: 20) and 583 germ-line mtDNA variations (novel: 53; reported: 530) in endometriosis patients. The A13603G, a novel missense mutation which leads to a substitution from serine to glycine at the codon 423 of ND5 gene showed 100% incidence in ectopic tissues. Interestingly, eutopic endometrium and peripheral leukocytes of all the patients showed heteroplasmy (A/G; 40–80%) at this locus, while their ectopic endometrium showed homoplasmic mutant allele (G/G). Superimposition of native and mutant structures of ND5 generated by homology modeling revealed no structural differences. Tumor-specific deletions and MSI were not observed in any of the ectopic tissues. Haplogrouping analysis showed a significant association between haplogroup M5 and endometriosis risk (P: 0.00069) after bonferroni correction.

Conclusions

Our findings substantiate the rationale for exploring the mitochondrial genome as a biomarker for the diagnosis of endometriosis.     

High-Resolution Mechanical Imaging of Glioblastoma by Multifrequency Magnetic Resonance Elastography

Objective

To generate high-resolution maps of the viscoelastic properties of human brain parenchyma for presurgical quantitative assessment in glioblastoma (GB).

Methods

Twenty-two GB patients underwent routine presurgical work-up supplemented by additional multifrequency magnetic resonance elastography. Two three-dimensional viscoelastic parameter maps, magnitude |G*|, and phase angle φ of the complex shear modulus were reconstructed by inversion of full wave field data in 2-mm isotropic resolution at seven harmonic drive frequencies ranging from 30 to 60 Hz.

Results

Mechanical brain maps confirmed that GB are composed of stiff and soft compartments, resulting in high intratumor heterogeneity. GB could be easily differentiated from healthy reference tissue by their reduced viscous behavior quantified by φ (0.37±0.08 vs. 0.58±0.07). |G*|, which in solids more relates to the material''s stiffness, was significantly reduced in GB with a mean value of 1.32±0.26 kPa compared to 1.54±0.27 kPa in healthy tissue (P = 0.001). However, some GB (5 of 22) showed increased stiffness.

Conclusion

GB are generally less viscous and softer than healthy brain parenchyma. Unrelated to the morphology-based contrast of standard magnetic resonance imaging, elastography provides an entirely new neuroradiological marker and contrast related to the biomechanical properties of tumors.     

KRAS Genotypic Changes of Circulating Tumor Cells during Treatment of Patients with Metastatic Colorectal Cancer

Introduction

Circulating tumor cells (CTCs) could represent a non-invasive source of cancer cells used for longitudinal monitoring of the tumoral mutation status throughout the course of the disease. The aims of the present study were to investigate the detection of KRAS mutations in CTCs from patients with metastatic colorectal cancer (mCRC) and to compare their mutation status during treatment or disease progression with that of the corresponding primary tumors.

Materials and Methods

Identification of the seven most common KRAS mutations on codons 12 and 13 was performed by Peptide Nucleic Acid (PNA)-based qPCR method. The sensitivity of the assay was determined after isolation of KRAS mutant cancer cells spiked into healthy donors'' blood, using the CellSearch Epithelial Cell kit. Consistent detection of KRAS mutations was achieved in samples containing at least 10 tumor cells/7.5 ml of blood.

Results

The clinical utility of the assay was assessed in 48 blood samples drawn from 31 patients with mCRC. All patients had PIK3CA and BRAF wild type primary tumors and 14 KRAS mutant tumors. CTCs were detected in 65% of specimens obtained from 74% of patients. KRAS mutation analysis in CTC-enriched specimens showed that 45% and 16.7% of patients with mutant and wild type primary tumors, respectively, had detectable mutations in their CTCs. Assessing KRAS mutations in serial blood samples revealed that individual patient''s CTCs exhibited different mutational status of KRAS during treatment.

Conclusions

The current findings support the rationale for using the CTCs as a dynamic source of tumor cells which, by re-evaluating their KRAS mutation status, could predict, perhaps more accurately, the response of mCRC patients to targeted therapy.     

Tumor Phosphatidylinositol-3-Kinase Signaling and Development of Metastatic Disease in Locally Advanced Rectal Cancer

Background

Recognizing EGFR as key orchestrator of the metastatic process in colorectal cancer, but also the substantial heterogeneity of responses to anti-EGFR therapy, we examined the pattern of composite tumor kinase activities governed by EGFR-mediated signaling that might be implicated in development of metastatic disease.

Patients and Methods

Point mutations in KRAS, BRAF, and PIK3CA and ERBB2 amplification were determined in primary tumors from 63 patients with locally advanced rectal cancer scheduled for radical treatment. Using peptide arrays with tyrosine kinase substrates, ex vivo phosphopeptide profiles were generated from the same baseline tumor samples and correlated to metastasis-free survival.

Results

Unsupervised clustering analysis of the resulting phosphorylation of 102 array substrates defined two tumor classes, both consisting of cases with and without KRAS/BRAF mutations. The smaller cluster group of patients, with tumors generating high ex vivo phosphorylation of phosphatidylinositol-3-kinase-related substrates, had a particularly aggressive disease course, with almost a half of patients developing metastatic disease within one year of follow-up.

Conclusion

High phosphatidylinositol-3-kinase-mediated signaling activity of the primary tumor, rather than KRAS/BRAF mutation status, was identified as a hallmark of poor metastasis-free survival in patients with locally advanced rectal cancer undergoing radical treatment of the pelvic cavity.     

Characterization of novel StAR (steroidogenic acute regulatory protein) mutations causing non-classic lipoid adrenal hyperplasia

Context

Steroidogenic acute regulatory protein (StAR) is crucial for transport of cholesterol to mitochondria where biosynthesis of steroids is initiated. Loss of StAR function causes lipoid congenital adrenal hyperplasia (LCAH).

Objective

StAR gene mutations causing partial loss of function manifest atypical and may be mistaken as familial glucocorticoid deficiency. Only a few mutations have been reported.

Design

To report clinical, biochemical, genetic, protein structure and functional data on two novel StAR mutations, and to compare them with published literature.

Setting

Collaboration between the University Children''s Hospital Bern, Switzerland, and the CIBERER, Hospital Vall d''Hebron, Autonomous University, Barcelona, Spain.

Patients

Two subjects of a non-consanguineous Caucasian family were studied. The 46,XX phenotypic normal female was diagnosed with adrenal insufficiency at the age of 10 months, had normal pubertal development and still has no signs of hypergonodatropic hypogonadism at 32 years of age. Her 46,XY brother was born with normal male external genitalia and was diagnosed with adrenal insufficiency at 14 months. Puberty was normal and no signs of hypergonadotropic hypogonadism are present at 29 years of age.

Results

StAR gene analysis revealed two novel compound heterozygote mutations T44HfsX3 and G221S. T44HfsX3 is a loss-of-function StAR mutation. G221S retains partial activity (∼30%) and is therefore responsible for a milder, non-classic phenotype. G221S is located in the cholesterol binding pocket and seems to alter binding/release of cholesterol.

Conclusions

StAR mutations located in the cholesterol binding pocket (V187M, R188C, R192C, G221D/S) seem to cause non-classic lipoid CAH. Accuracy of genotype-phenotype prediction by in vitro testing may vary with the assays employed.     

Deletion of Chromosomes 13q and 14q Is a Common Feature of Tumors with BRCA2 Mutations

Introduction

Germline BRCA1 or BRCA2 mutations account for 20–30% of familial clustering of breast cancer. The main indication for BRCA2 screening is currently the family history but the yield of mutations identified in patients selected this way is low.

Methods

To develop more efficient approaches to screening we have compared the gene expression and genomic profiles of BRCA2-mutant breast tumors with those of breast tumors lacking BRCA1 or BRCA2 mutations.

Results

We identified a group of 66 genes showing differential expression in our training set of 7 BRCA2-mutant tumors and in an independent validation set of 19 BRCA2-mutant tumors. The differentially expressed genes include a prominent cluster of genes from chromosomes 13 and 14 whose expression is reduced. Gene set enrichment analysis confirmed that genes in specific bands on 13q and 14q showed significantly reduced expression, suggesting that the affected bands may be preferentially deleted in BRCA2-mutant tumors. Genomic profiling showed that the BRCA2-mutant tumors indeed harbor deletions on chromosomes 13q and 14q. To exploit this information we have created a simple fluorescence in situ hybridization (FISH) test and shown that it detects tumors with deletions on chromosomes 13q and 14q.

Conclusion

Together with previous reports, this establishes that deletions on chromosomes 13q and 14q are a hallmark of BRCA2-mutant tumors. We propose that FISH to detect these deletions would be an efficient and cost-effective first screening step to identify potential BRCA2-mutation carriers among breast cancer patients without a family history of breast cancer.     

High-Throughput Mutation Profiling of Primary and Metastatic Endometrial Cancers Identifies KRAS,FGFR2 and PIK3CA to Be Frequently Mutated

Background

Despite being the most common pelvic gynecologic malignancy in industrialized countries, no targeted therapies are available for patients with metastatic endometrial carcinoma. In order to improve treatment, underlying molecular characteristics of primary and metastatic disease must be explored.

Methodology/Principal Findings

We utilized the mass spectrometric-based mutation detection technology OncoMap to define the types and frequency of point somatic mutations in endometrial cancer. 67 primary tumors, 15 metastases corresponding to 7 of the included primary tumors and 11 endometrial cancer cell lines were screened for point mutations in 28 known oncogenes. We found that 27 (40.3%) of 67 primary tumors harbored one or more mutations with no increase in metastatic lesions. FGFR2, KRAS and PIK3CA were consistently the most frequently mutated genes in primary tumors, metastatic lesions and cell lines.

Conclusions/Significance

Our results emphasize the potential for targeting FGFR2, KRAS and PIK3CA mutations in endometrial cancer for development of novel therapeutic strategies.     

Clonal Analysis in Recurrent Astrocytic, Oligoastrocytic and Oligodendroglial Tumors Implicates IDH1- Mutation as Common Tumor Initiating Event

Background

To investigate the dynamics of inter- and intratumoral molecular alterations during tumor progression in recurrent gliomas.

Methodology/Principal Findings

To address intertumoral heterogeneity we investigated non- microdissected tumor tissue of 106 gliomas representing 51 recurrent tumors. To address intratumoral heterogeneity a set of 16 gliomas representing 7 tumor pairs with at least one recurrence, and 4 single mixed gliomas were investigated by microdissection of distinct oligodendroglial and astrocytic tumor components. All tumors and tumor components were analyzed for allelic loss of 1p/19q (LOH 1p/19q), for TP53- mutations and for R132 mutations in the IDH1 gene. The investigation of non- microdissected tumor tissue revealed clonality in 75% (38/51). Aberrant molecular alterations upon recurrence were detected in 25% (13/51). 64% (9/14) of these were novel and associated with tumor progression. Loss of previously detected alterations was observed in 36% (5/14). One tumor pair (1/14; 7%) was significant for both. Intratumoral clonality was detected in 57% (4/7) of the microdissected tumor pairs and in 75% (3/4) of single microdissected tumors. 43% (3/7) of tumor pairs and one single tumor (25%) revealed intratumoral heterogeneity. While intratumoral heterogeneity affected both the TP53- mutational status and the LOH1p/19q status, all tumors with intratumoral heterogeneity shared the R132 IDH1- mutation as a common feature in both their microdissected components.

Conclusions/Significance

The majority of recurrent gliomas are of monoclonal origin. However, the detection of divertive tumor cell clones in morphological distinct tumor components sharing IDH1- mutations as early event may provide insight into the tumorigenesis of true mixed gliomas.     

Towards accurate characterization of clonal heterogeneity based on structural variation

Background

Recent advances in deep digital sequencing have unveiled an unprecedented degree of clonal heterogeneity within a single tumor DNA sample. Resolving such heterogeneity depends on accurate estimation of fractions of alleles that harbor somatic mutations. Unlike substitutions or small indels, structural variants such as deletions, duplications, inversions and translocations involve segments of DNAs and are potentially more accurate for allele fraction estimations. However, no systematic method exists that can support such analysis.

Results

In this paper, we present a novel maximum-likelihood method that estimates allele fractions of structural variants integratively from various forms of alignment signals. We develop a tool, BreakDown, to estimate the allele fractions of most structural variants including medium size (from 1 kilobase to 1 megabase) deletions and duplications, and balanced inversions and translocations.

Conclusions

Evaluation based on both simulated and real data indicates that our method systematically enables structural variants for clonal heterogeneity analysis and can greatly enhance the characterization of genomically instable tumors.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2105-15-299) contains supplementary material, which is available to authorized users.     

Assessment of EGFR Mutations in Circulating Tumor Cell Preparations from NSCLC Patients by Next Generation Sequencing: Toward a Real-Time Liquid Biopsy for Treatment

Introduction

Assessment of EGFR mutation in non-small cell lung cancer (NSCLC) patients is mandatory for optimization of pharmacologic treatment. In this respect, mutation analysis of circulating tumor cells (CTCs) may be desirable since they may provide real-time information on patient''s disease status.

Experimental Design

Blood samples were collected from 37 patients enrolled in the TRIGGER study, a prospective phase II multi-center trial of erlotinib treatment in advanced NSCLC patients with activating EGFR mutations in tumor tissue. 10 CTC preparations from breast cancer patients without EGFR mutations in their primary tumors and 12 blood samples from healthy subjects were analyzed as negative controls. CTC preparations, obtained by the Veridex CellSearch System, were subjected to ultra-deep next generation sequencing (NGS) on the Roche 454 GS junior platform.

Results

CTCs fulfilling all Veridex criteria were present in 41% of the patients examined, ranging in number between 1 and 29. In addition to validated CTCs, potential neoplastic elements were seen in 33 cases. These included cells not fulfilling all Veridex criteria (also known as “suspicious objects”) found in 5 (13%) of 37 cases, and isolated or clustered large naked nuclei with irregular shape observed in 33 (89%) cases. EGFR mutations were identified by NGS in CTC preparations of 31 (84%) patients, corresponding to those present in matching tumor tissue. Twenty-five (96%) of 26 deletions at exon 19 and 6 (55%) of 11 mutations at exon 21 were detectable (P = 0.005). In 4 (13%) cases, multiple EGFR mutations, suggesting CTC heterogeneity, were documented. No mutations were found in control samples.

Conclusions

We report for the first time that the CellSearch System coupled with NGS is a very sensitive and specific diagnostic tool for EGFR mutation analysis in CTC preparations with potential clinical impact.     

Development of synchronous VHL syndrome tumors reveals contingencies and constraints to tumor evolution

Background

Genomic analysis of multi-focal renal cell carcinomas from an individual with a germline VHL mutation offers a unique opportunity to study tumor evolution.

Results

We perform whole exome sequencing on four clear cell renal cell carcinomas removed from both kidneys of a patient with a germline VHL mutation. We report that tumors arising in this context are clonally independent and harbour distinct secondary events exemplified by loss of chromosome 3p, despite an identical genetic background and tissue microenvironment. We propose that divergent mutational and copy number anomalies are contingent upon the nature of 3p loss of heterozygosity occurring early in tumorigenesis. However, despite distinct 3p events, genomic, proteomic and immunohistochemical analyses reveal evidence for convergence upon the PI3K-AKT-mTOR signaling pathway. Four germline tumors in this young patient, and in a second, older patient with VHL syndrome demonstrate minimal intra-tumor heterogeneity and mutational burden, and evaluable tumors appear to follow a linear evolutionary route, compared to tumors from patients with sporadic clear cell renal cell carcinoma.

Conclusions

In tumors developing from a germline VHL mutation, the evolutionary principles of contingency and convergence in tumor development are complementary. In this small set of patients with early stage VHL-associated tumors, there is reduced mutation burden and limited evidence of intra-tumor heterogeneity.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0433-z) contains supplementary material, which is available to authorized users.