Pancreatic cancer is a highly lethal malignancy with few effective therapies. We performed exome sequencing and copy number analysis to define genomic aberrations in a prospectively accrued clinical cohort (n = 142) of early (stage I and II) sporadic pancreatic ductal adenocarcinoma. Detailed analysis of 99 informative tumours identified substantial heterogeneity with 2,016 non-silent mutations and 1,628 copy-number variations. We define 16 significantly mutated genes, reaffirming known mutations (KRAS, TP53, CDKN2A, SMAD4, MLL3, TGFBR2, ARID1A and SF3B1), and uncover novel mutated genes including additional genes involved in chromatin modification (EPC1 and ARID2), DNA damage repair (ATM) and other mechanisms (ZIM2, MAP2K4, NALCN, SLC16A4 and MAGEA6). Integrative analysis with in vitro functional data and animal models provided supportive evidence for potential roles for these genetic aberrations in carcinogenesis. Pathway-based analysis of recurrently mutated genes recapitulated clustering in core signalling pathways in pancreatic ductal adenocarcinoma, and identified new mutated genes in each pathway. We also identified frequent and diverse somatic aberrations in genes described traditionally as embryonic regulators of axon guidance, particularly SLIT/ROBO signalling, which was also evident in murine Sleeping Beauty transposon-mediated somatic mutagenesis models of pancreatic cancer, providing further supportive evidence for the potential involvement of axon guidance genes in pancreatic carcinogenesis.

SUMMARY We performed integrated genomic, transcriptomic and proteomic profiling of 150 pancreatic ductal adenocarcinoma (PDAC) specimens, including samples with characteristic low neoplastic cellularity. Deep whole-exome sequencing revealed recurrent somatic mutations in KRAS, TP53, CDKN2A, SMAD4, RNF43, ARID1A, TGFβR2, GNAS, RREB1 and PBRM1. KRAS wild-type tumors harbored alterations in other oncogenic drivers, including GNAS, BRAF, CTNNB1 and additional RAS pathway genes. A subset of tumors harbored multiple KRAS mutations, with some showing evidence of biallelic mutations. Protein profiling identified a favorable prognosis subset with low epithelial-mesenchymal transition and high MTOR pathway scores. Associations of non-coding RNAs with tumor-specific mRNA subtypes were also identified. Our integrated multi-platform analysis reveals a complex molecular landscape of PDAC and provides a roadmap for precision medicine.

No abstract

Cancer cells adapt metabolically to proliferate under nutrient limitation. Here we used combined transcriptional-metabolomic network analysis to identify metabolic pathways that support glucose-independent tumor cell proliferation. We found that glucose deprivation stimulated re-wiring of the tricarboxylic acid (TCA) cycle and early steps of gluconeogenesis to promote glucose-independent cell proliferation. Glucose limitation promoted the production of phosphoenolpyruvate (PEP) from glutamine via the activity of mitochondrial PEP-carboxykinase (PCK2). Under these conditions, glutamine-derived PEP was used to fuel biosynthetic pathways normally sustained by glucose, including serine and purine biosynthesis. PCK2 expression was required to maintain tumor cell proliferation under limited-glucose conditions in vitro and tumor growth in vivo. Elevated PCK2 expression is observed in several human tumor types and enriched in tumor tissue from non-small-cell lung cancer (NSCLC) patients. Our results define a role for PCK2 in cancer cell metabolic reprogramming that promotes glucose-independent cell growth and metabolic stress resistance in human tumors.

Whole-genome sequencing of patient DNA can facilitate diagnosis of a disease, but its potential for guiding treatment has been under-realized. We interrogated the complete genome sequences of a 14-year-old fraternal twin pair diagnosed with dopa (3,4-dihydroxyphenylalanine)–responsive dystonia (DRD; Mendelian Inheritance in Man #128230). DRD is a genetically heterogeneous and clinically complex movement disorder that is usually treated with l-dopa, a precursor of the neurotransmitter dopamine. Whole-genome sequencing identified compound heterozygous mutations in the SPR gene encoding sepiapterin reductase. Disruption of SPR causes a decrease in tetrahydrobiopterin, a cofactor required for the hydroxylase enzymes that synthesize the neurotransmitters dopamine and serotonin. Supplementation of l-dopa therapy with 5-hydroxytryptophan, a serotonin precursor, resulted in clinical improvements in both twins.

Genitopatellar syndrome (GPS) is a skeletal dysplasia with cerebral and genital anomalies for which the molecular basis has not yet been determined. By exome sequencing, we found de novo heterozygous truncating mutations in KAT6B (lysine acetyltransferase 6B, formerly known as MYST4 and MORF) in three subjects; then by Sanger sequencing of KAT6B, we found similar mutations in three additional subjects. The mutant transcripts do not undergo nonsense-mediated decay in cells from subjects with GPS. In addition, human pathological analyses and mouse expression studies point to systemic roles of KAT6B in controlling organismal growth and development. Myst4 (the mouse orthologous gene) is expressed in mouse tissues corresponding to those affected by GPS. Phenotypic differences and similarities between GPS, the Say-Barber-Biesecker variant of Ohdo syndrome (caused by different mutations of KAT6B), and Rubinstein-Taybi syndrome (caused by mutations in other histone acetyltransferases) are discussed. Together, the data support an epigenetic dysregulation of the limb, brain, and genital developmental programs.

The Warburg effect is a tumorigenic metabolic adaptation process characterized by augmented aerobic glycolysis, which enhances cellular bioenergetics. In normal cells, energy homeostasis is controlled by AMPK; however, its role in cancer is not understood, as both AMPK-dependent tumor-promoting and -inhibiting functions were reported. Upon stress, energy levels are maintained by increased mitochondrial biogenesis and glycolysis, controlled by transcriptional coactivator PGC-1α and HIF, respectively. In normoxia, AMPK induces PGC-1α, but how HIF is activated is unclear. Germline mutations in the gene encoding the tumor suppressor folliculin (FLCN) lead to Birt-Hogg-Dubé (BHD) syndrome, which is associated with an increased cancer risk. FLCN was identified as an AMPK binding partner, and we evaluated its role with respect to AMPKdependent energy functions. We revealed that loss of FLCN constitutively activates AMPK, resulting in PGC-1α-mediated mitochondrial biogenesis and increased ROS production. ROS induced HIF transcriptional activity and drove Warburg metabolic reprogramming, coupling AMPK-dependent mitochondrial biogenesis to HIF-dependent metabolic changes. This reprogramming stimulated cellular bioenergetics and conferred a HIF-dependent tumorigenic advantage in FLCN-negative cancer cells. Moreover, this pathway is conserved in a BHD-derived tumor. These results indicate that FLCN inhibits tumorigenesis by preventing AMPK-dependent HIF activation and the subsequent Warburg metabolic transformation.

No abstract