Proteomic Analysis of Signaling Network Regulation in Renal Cell Carcinomas with Differential Hypoxia-Inducible Factor-2α Expression

Background

The loss of von Hippel–Lindau (VHL) protein function leads to highly vascular renal tumors characterized by an aggressive course of disease and refractoriness to chemotherapy and radiotherapy. Loss of VHL in renal tumors also differs from tumors of other organs in that the oncogenic cascade is mediated by an increase in the levels of hypoxia-inducible factor-2α (HIF2α) instead of hypoxia-inducible factor-1α (HIF1α).

Methods and Principal Findings

We used renal carcinoma cell lines that recapitulate the differences between mutant VHL and wild-type VHL genotypes. Utilizing a method relying on extracted peptide intensities as a label-free approach for quantitation by liquid chromatography–mass spectrometry, our proteomics study revealed regulation of key proteins important for cancer cell survival, proliferation and stress-resistance, and implicated differential regulation of signaling networks in VHL-mutant renal cell carcinoma. We also observed upregulation of cellular energy pathway enzymes and the stress-responsive mitochondrial 60-kDa heat shock protein. Finding reliance on glutaminolysis in VHL-mutant renal cell carcinoma was of particular significance, given the generally predominant dependence of tumors on glycolysis. The data have been deposited to the ProteomeXchange with identifier PXD000335.

Conclusions and Significance

Pathway analyses provided corroborative evidence for differential regulation of molecular and cellular functions influencing cancer energetics, metabolism and cell proliferation in renal cell carcinoma with distinct VHL genotype. Collectively, the differentially regulated proteome characterized by this study can potentially guide translational research specifically aimed at effective clinical interventions for advanced VHL-mutant, HIF2α-over-expressing tumors.     

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.     

VHL Type 2B Mutations Retain VBC Complex Form and Function

Background

von Hippel-Lindau disease is characterized by a spectrum of hypervascular tumors, including renal cell carcinoma, hemangioblastoma, and pheochromocytoma, which occur with VHL genotype-specific differences in penetrance. VHL loss causes a failure to regulate the hypoxia inducible factors (HIF-1α and HIF-2α), resulting in accumulation of both factors to high levels. Although HIF dysregulation is critical to VHL disease-associated renal tumorigenesis, increasing evidence points toward gradations of HIF dysregulation contributing to the degree of predisposition to renal cell carcinoma and other manifestations of the disease.

Methodology/Principal Findings

This investigation examined the ability of disease-specific VHL missense mutations to support the assembly of the VBC complex and to promote the ubiquitylation of HIF. Our interaction analysis supported previous observations that VHL Type 2B mutations disrupt the interaction between pVHL and Elongin C but maintain partial regulation of HIF. We additionally demonstrated that Type 2B mutant pVHL forms a remnant VBC complex containing the active members ROC1 and Cullin-2 which retains the ability to ubiquitylate HIF-1α.

Conclusions

Our results suggest that subtypes of VHL mutations support an intermediate level of HIF regulation via a remnant VBC complex. These findings provide a mechanism for the graded HIF dysregulation and genetic predisposition for cancer development in VHL disease.     

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.     

Systematic Study on Genetic and Epimutational Profile of a Cohort of Amsterdam Criteria-Defined Lynch Syndrome in Singapore

Background

Germline defects of mismatch repair (MMR) genes underlie Lynch Syndrome (LS). We aimed to gain comprehensive genetic and epigenetic profiles of LS families in Singapore, which will facilitate efficient molecular diagnosis of LS in Singapore and the region.

Methods

Fifty nine unrelated families were studied. Mutations in exons, splice-site junctions and promoters of five MMR genes were scanned by high resolution melting assay followed by DNA sequencing, large fragment deletions/duplications and promoter methylation in MLH1, MSH2, MSH6 and PMS2 were evaluated by multiplex ligation-dependent probe amplification. Tumor microsatellite instability (MSI) was assessed with five mononucleotide markers and immunohistochemical staining (IHC) was also performed.

Results

Pathogenic defects, all confined to MLH1 and MSH2, were identified in 17 out of 59 (28.8%) families. The mutational spectrum was highly heterogeneous and 28 novel variants were identified. One recurrent mutation in MLH1 (c.793C>T) was also observed. 92.9% sensitivity for indication of germline mutations conferred by IHC surpassed 64.3% sensitivity by MSI. Furthermore, 15.6% patients with MSS tumors harbored pathogenic mutations.

Conclusions

Among major ethnic groups in Singapore, all pathogenic germline defects were confined to MLH1 and MSH2. Caution should be applied when the Amsterdam criteria and consensus microsatellite marker panel recommended in the revised Bethesda guidelines are applied to the local context. We recommend a screening strategy for the local LS by starting with tumor IHC and the hotspot mutation testing at MLH1 c.793C>T followed by comprehensive mutation scanning in MLH1 and MSH2 prior to proceeding to other MMR genes.     

Household Income Is Associated with the p53 Mutation Frequency in Human Breast Tumors

Background

A study from Scotland reported that the p53 mutation frequency in breast tumors is associated with socio-economic deprivation.

Methods

We analyzed the association of the tumor p53 mutational status with tumor characteristics, education, and self-reported annual household income (HI) among 173 breast cancer patients from the greater Baltimore area, United States.

Results

p53 mutational frequency was significantly associated with HI. Patients with < $15,000 HI had the highest p53 mutation frequency (21%), followed by the income group between $15,000 and $60,000 (18%), while those above $60,000 HI had the fewest mutations (5%). When dichotomized at $60,000, 26 out of 135 patients in the low income category had acquired a p53 mutation, while only 2 out of 38 with a high income carried a mutation (P < 0.05). In the adjusted logistic regression analysis with 3 income categories (trend test), the association between HI and p53 mutational status was independent of tumor characteristics, age, race/ethnicity, tobacco smoking and body mass. Further analyses revealed that HI may impact the p53 mutational frequency preferentially in patients who develop an estrogen receptor (ER)-negative disease. Within this group, 42% of the low income patients (< $15,000 HI) carried a mutation, followed by the middle income group (21%), while those above $60,000 HI did not carry mutations (P _trend < 0.05).

Conclusions:

HI is associated with the p53 mutational frequency in patients who develop an ER-negative disease. Furthermore, high income patients may acquire fewer p53 mutations than other patients, suggesting that lifetime exposures associated with socio-economic status may impact breast cancer biology.     

The Specificity of the FOXL2 c.402C>G Somatic Mutation: A Survey of Solid Tumors

Background

A somatic mutation in the FOXL2 gene is reported to be present in almost all (97%; 86/89) morphologically defined, adult-type, granulosa-cell tumors (A-GCTs). This FOXL2 c.402C>G mutation changes a highly conserved cysteine residue to a tryptophan (p.C134W). It was also found in a minority of other ovarian malignant stromal tumors, but not in benign ovarian stromal tumors or unrelated ovarian tumors or breast cancers.

Methodology/Principal Findings

Herein we studied other cancers and cell lines for the presence of this mutation. We screened DNA from 752 tumors of epithelial and mesenchymal origin and 28 ovarian cancer cell lines and 52 other cancer cell lines of varied origin. We found the FOXL2 c.402C>G mutation in an unreported A-GCT case and the A-GCT-derived cell line KGN. All other tumors and cell lines analyzed were mutation negative.

Conclusions/Significance

In addition to proving that the KGN cell line is a useful model to study A-GCTs, these data show that the c.402C>G mutation in FOXL2 is not commonly found in a wide variety of other cancers and therefore it is likely pathognomonic for A-GCTs and closely related tumors.     

The UBC-40 Urothelial Bladder Cancer cell line index: a genomic resource for functional studies

Background

Urothelial bladder cancer is a highly heterogeneous disease. Cancer cell lines are useful tools for its study. This is a comprehensive genomic characterization of 40 urothelial bladder carcinoma (UBC) cell lines including information on origin, mutation status of genes implicated in bladder cancer (FGFR3, PIK3CA, TP53, and RAS), copy number alterations assessed using high density SNP arrays, uniparental disomy (UPD) events, and gene expression.

Results

Based on gene mutation patterns and genomic changes we identify lines representative of the FGFR3-driven tumor pathway and of the TP53/RB tumor suppressor-driven pathway. High-density array copy number analysis identified significant focal gains (1q32, 5p13.1-12, 7q11, and 7q33) and losses (i.e. 6p22.1) in regions altered in tumors but not previously described as affected in bladder cell lines. We also identify new evidence for frequent regions of UPD, often coinciding with regions reported to be lost in tumors. Previously undescribed chromosome X losses found in UBC lines also point to potential tumor suppressor genes. Cell lines representative of the FGFR3-driven pathway showed a lower number of UPD events.

Conclusions

Overall, there is a predominance of more aggressive tumor subtypes among the cell lines. We provide a cell line classification that establishes their relatedness to the major molecularly-defined bladder tumor subtypes. The compiled information should serve as a useful reference to the bladder cancer research community and should help to select cell lines appropriate for the functional analysis of bladder cancer genes, for example those being identified through massive parallel sequencing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1450-3) contains supplementary material, which is available to authorized users.     

Microsatellite abundance across the Anthozoa and Hydrozoa in the phylum Cnidaria

Background

Microsatellite loci have high mutation rates and thus are indicative of mutational processes within the genome. By concentrating on the symbiotic and aposymbiotic cnidarians, we investigated if microsatellite abundances follow a phylogenetic or ecological pattern. Individuals from eight species were shotgun sequenced using 454 GS-FLX Titanium technology. Sequences from the three available cnidarian genomes (Nematostella vectensis, Hydra magnipapillata and Acropora digitifera) were added to the analysis for a total of eleven species representing two classes, three subclasses and eight orders within the phylum Cnidaria.

Results

Trinucleotide and tetranucleotide repeats were the most abundant motifs, followed by hexa- and dinucleotides. Pentanucleotides were the least abundant motif in the data set. Hierarchical clustering and log likelihood ratio tests revealed a weak relationship between phylogeny and microsatellite content. Further, comparisons between cnidaria harboring intracellular dinoflagellates and those that do not, show microsatellite coverage is higher in the latter group.

Conclusions

Our results support previous studies that found tri- and tetranucleotides to be the most abundant motifs in invertebrates. Differences in microsatellite coverage and composition between symbiotic and non-symbiotic cnidaria suggest the presence/absence of dinoflagellates might place restrictions on the host genome.

Electronic supplementary material

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

A Pathogenic Mosaic TP53 Mutation in Two Germ Layers Detected by Next Generation Sequencing

Background

Li-Fraumeni syndrome is caused by germline TP53 mutations and is clinically characterized by a predisposition to a range of cancers, most commonly sarcoma, brain tumours and leukemia. Pathogenic mosaic TP53 mutations have only rarely been described.

Methods and Findings

We describe a 2 years old child presenting with three separate cancers over a 6 month period; two soft tissue mesenchymal tumors and an aggressive metastatic neuroblastoma. As conventional testing of blood DNA by Sanger sequencing for mutations in TP53, ALK, and SDH was negative, whole exome sequencing of the blood DNA of the patient and both parents was performed to screen more widely for cancer predisposing mutations. In the patient''s but not the parents'' DNA we found a c.743 G>A, p.Arg248Gln (CCDS11118.1) TP53 mutation in 3–20% of sequencing reads, a level that would not generally be detectable by Sanger sequencing. Homozygosity for this mutation was detected in all tumor samples analyzed, and germline mosaicism was demonstrated by analysis of the child''s newborn blood spot DNA. The occurrence of separate tumors derived from different germ layers suggests that this de novo mutation occurred early in embryogenesis, prior to gastrulation.

Conclusion

The case demonstrates pathogenic mosaicim, detected by next generation deep sequencing, that arose in the early stages of embryogenesis.     

The genetic heterogeneity and mutational burden of engineered melanomas in zebrafish models

Background

Melanoma is the most deadly form of skin cancer. Expression of oncogenic BRAF or NRAS, which are frequently mutated in human melanomas, promote the formation of nevi but are not sufficient for tumorigenesis. Even with germline mutated p53, these engineered melanomas present with variable onset and pathology, implicating additional somatic mutations in a multi-hit tumorigenic process.

Results

To decipher the genetics of these melanomas, we sequence the protein coding exons of 53 primary melanomas generated from several BRAF^V600E or NRAS^Q61K driven transgenic zebrafish lines. We find that engineered zebrafish melanomas show an overall low mutation burden, which has a strong, inverse association with the number of initiating germline drivers. Although tumors reveal distinct mutation spectrums, they show mostly C > T transitions without UV light exposure, and enrichment of mutations in melanogenesis, p53 and MAPK signaling. Importantly, a recurrent amplification occurring with pre-configured drivers BRAF^V600E and p53^-/- suggests a novel path of BRAF cooperativity through the protein kinase A pathway.

Conclusion

This is the first analysis of a melanoma mutational landscape in the absence of UV light, where tumors manifest with remarkably low mutation burden and high heterogeneity. Genotype specific amplification of protein kinase A in cooperation with BRAF and p53 mutation suggests the involvement of melanogenesis in these tumors. This work is important for defining the spectrum of events in BRAF or NRAS driven melanoma in the absence of UV light, and for informed exploitation of models such as transgenic zebrafish to better understand mechanisms leading to human melanoma formation.     

Comparative sequencing analysis reveals high genomic concordance between matched primary and metastatic colorectal cancer lesions

Background

Colorectal cancer is the second leading cause of cancer death in the United States, with over 50,000 deaths estimated in 2014. Molecular profiling for somatic mutations that predict absence of response to anti-EGFR therapy has become standard practice in the treatment of metastatic colorectal cancer; however, the quantity and type of tissue available for testing is frequently limited. Further, the degree to which the primary tumor is a faithful representation of metastatic disease has been questioned. As next-generation sequencing technology becomes more widely available for clinical use and additional molecularly targeted agents are considered as treatment options in colorectal cancer, it is important to characterize the extent of tumor heterogeneity between primary and metastatic tumors.

Results

We performed deep coverage, targeted next-generation sequencing of 230 key cancer-associated genes for 69 matched primary and metastatic tumors and normal tissue. Mutation profiles were 100% concordant for KRAS, NRAS, and BRAF, and were highly concordant for recurrent alterations in colorectal cancer. Additionally, whole genome sequencing of four patient trios did not reveal any additional site-specific targetable alterations.

Conclusions

Colorectal cancer primary tumors and metastases exhibit high genomic concordance. As current clinical practices in colorectal cancer revolve around KRAS, NRAS, and BRAF mutation status, diagnostic sequencing of either primary or metastatic tissue as available is acceptable for most patients. Additionally, consistency between targeted sequencing and whole genome sequencing results suggests that targeted sequencing may be a suitable strategy for clinical diagnostic applications.

Electronic supplementary material

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

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.     

Bimodal signatures of germline methylation are linked with gene expression plasticity in the coral Acropora millepora

Background

In invertebrates, genes belonging to dynamically regulated functional categories appear to be less methylated than “housekeeping” genes, suggesting that DNA methylation may modulate gene expression plasticity. To date, however, experimental evidence to support this hypothesis across different natural habitats has been lacking.

Results

Gene expression profiles were generated from 30 pairs of genetically identical fragments of coral Acropora millepora reciprocally transplanted between distinct natural habitats for 3 months. Gene expression was analyzed in the context of normalized CpG content, a well-established signature of historical germline DNA methylation. Genes with weak methylation signatures were more likely to demonstrate differential expression based on both transplant environment and population of origin than genes with strong methylation signatures. Moreover, the magnitude of expression differences due to environment and population were greater for genes with weak methylation signatures.

Conclusions

Our results support a connection between differential germline methylation and gene expression flexibility across environments and populations. Studies of phylogenetically basal invertebrates such as corals will further elucidate the fundamental functional aspects of gene body methylation in Metazoa.

Electronic supplementary material

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

Genetic control of renal tumorigenesis by the mouse Rtm1 locus

Background

The genetic basis of susceptibility to renal tumorigenesis has not yet been established in mouse strains. Mouse lines derived by bidirectional phenotypic selection on the basis of their maximal (AIRmax) or minimal (AIRmin) acute inflammatory responsiveness differ widely in susceptibility to spontaneous and urethane-induced renal tumorigenesis. To map the functional loci modulating renal tumor susceptibility in these mice, we carried out a genome-wide genetic linkage study, using SNP arrays, in an (AIRmax x AIRmin)F2 intercross population treated with a single urethane dose at 1 week of age and phenotyped for renal tumors at 35 weeks of age.

Results

AIRmax mice did not develop renal tumors spontaneously nor in response to urethane, whereas in AIRmin mice renal tumors formed spontaneously (in 52% of animals) and after urethane induction (89%). The tumors had a papillary morphology and were positive for alpha-methylacyl-CoA racemase and negative for CD10. By analysis of 879 informative SNPs in 662 mice, we mapped a single quantitative trait locus modulating the incidence of renal tumors in the (AIRmax x AIRmin)F2 intercross population. This locus, which we named Renal tumor modifier QTL 1 (Rtm1), mapped to chromosome 17 at 23.4 Mb (LOD score = 15.8), with SNPs rs3696835 and rs3719497 flanking the LOD score peak. The A allele of rs3719497 from AIRmin mice was associated with a 2.5-fold increased odds ratio for renal tumor development. The LOD score peak included the Tuberous sclerosis 2 (Tsc2) gene which has already been implicated in kidney disease: loss of function by germline retroviral insertion is associated with spontaneous renal tumorigenesis in the Eker rat, and heterozygous-null Tsc2^(+/-) mice develop renal cystadenomas.

Conclusions

We mapped Rtm1 as a single major locus modulating renal tumorigenesis in a murine intercross population. Thus, the AIR mouse lines can be considered a new genetic model for studying the role of germline and somatic molecular alterations in kidney neoplastic disease.     

Detection and Clinical Significance of Intratumoral EGFR Mutational Heterogeneity in Chinese Patients with Advanced Non-Small Cell Lung Cancer

Purpose

This study evaluated occurrence and potential clinical significance of intratumoral EGFR mutational heterogeneity in Chinese patients with non-small cell lung cancer (NSCLC).

Materials and Methods

Eighty-five stage IIIa-IV NSCLC patients who had undergone palliative surgical resection were included in this study. Of these, 45 patients carried EGFR mutations (group-M) and 40 patients were wild-type (group-W). Each tumor sample was microdissected to yield 28–34 tumor foci and Intratumoral EGFR mutation were determined using Denaturing High Performance Liquid Chromatography (DHPLC) and Amplification Refractory Mutation System (ARMS). EGFR copy numbers were measured using fluorescence in situ hybridization (FISH).

Results

Microdissection yielded 1,431 tumor foci from EGFR mutant patients (group-M) and 1,238 foci from wild-type patients (group-W). The EGFR mutant frequencies in group-M were 80.6% (1,154/1,431) and 87.1% (1,247/1,431) using DHPLC and ARMS, respectively. A combination of EGFR-mutated and wild-type cells was detected in 32.9% (28/85) of samples by DHPLC and 28.2% (24/85) by ARMS, supporting the occurrence of intratumoral heterogeneity. Thirty-one patients (36.5%) were identified as EGFR FISH-positive. Patients harboring intratumoral mutational heterogeneity possessed lower EGFR copy numbers than those tumors contained mutant cells alone (16.7% vs. 71.0%, P<0.05). Among 26 patients who had received EGFR-TKIs, the mean EGFR mutation content was higher in patients showing partial response (86.1%) or stable disease (48.7%) compared with patients experiencing progressive disease (6.0%) (P = 0.001). There also showed relationship between progression-free survival (PFS) and different content of EGFR mutation groups (pure wild type EGFR, EGFR mutation with heterogeneity and pure mutated EGFR) (P = 0.001).

Conclusion

Approximately 30% of patients presented intratumoral EGFR mutational heterogeneity, accompanying with relatively low EGFR copy number. EGFR mutant content was correlated with the response and prognosis of EGFR-TKIs.