Genomic Landscapes of Pancreatic Neoplasia
Article information
Abstract
Pancreatic cancer is a deadly disease with a dismal prognosis. However, recent advances in sequencing and bioinformatic technology have led to the systematic characterization of the genomes of all major tumor types in the pancreas. This characterization has revealed the unique genomic landscape of each tumor type. This knowledge will pave the way for improved diagnostic and therapeutic approaches to pancreatic tumors that take advantage of the genetic alterations in these neoplasms.
Decades of cancer research have shown that cancer is a genetic disease caused by the accumulation of somatic, and in some cases inherited, mutations in oncogenes and tumor suppressor genes. Pancreatic neoplasms are some of the best characterized at the genomic level. Some of the key genetic drivers in pancreatic neoplasia have been known for years, but the introduction of high throughput sequencing has led to a more complete understanding of pancreatic cancer genomes. The exomes of all of the major tumor types in the pancreas have been analyzed by whole exome sequencing. This systematic genomic characterization has identified unique genetic signatures in each morphologically defined pancreatic tumor type, demonstrating that each type of neoplasm is driven by a distinct set of genetic alterations. A huge opportunity now exists to translate this new knowledge to improve patient care.
PANCREATIC DUCTAL ADENOCARCINOMA AND ITS VARIANTS
Pancreatic ductal adenocarcinoma (PDAC) is the most common pancreatic malignancy, and with a 5-year survival of only 6%, it is one of the deadliest of all human cancers. The KRAS gene, which encodes a small GTPase that mediates downstream signaling from growth factor receptors, is the most commonly mutated oncogene in PDAC [1-3]. Missense mutations in KRAS cluster in specific hotspots (most commonly codon 12), consistent with its role as an oncogene—KRAS mutations occur in >90% of PDACs. In addition to hotspot mutations in the KRAS oncogene, three tumor suppressor genes are frequently mutated in PDAC. The tumor suppressor gene P16/CDKN2A, which encodes a critical cell cycle regulator, is inactivated in >90% of PDACs by several mechanisms, including intragenic mutation coupled with loss of heterozygosity, homozygous deletion, and promoter methylation [4]. Mutations in TP53, a key component of the cellular stress response, are also common in PDAC, reported in approximately 75% of PDACs, most commonly by small intragenic mutation coupled with loss of heterozygosity [5,6]. Mutations in TP53 often result in strong diffuse nuclear expression of p53 protein which can be detected by immunohistochemistry (Fig. 1A). Somatic inactivation of SMAD4, a tumor suppressor gene that codes for a component of the transforming growth factor beta signaling pathway, occurs in approximately 55% of PDACs, usually by homozygous deletion or intragenic mutation coupled with loss of heterozygosity [7,8]. These mutations can also be detected by immunohistochemistry, as they cause loss of Smad4 protein expression (Fig. 1B). Mutations in SMAD4 have potential clinical implications, as these mutations are associated with worse prognosis and widespread metastases (rather than local disease) [9,10]. Although these four genes (often referred to as the four “mountains”) are the most commonly mutated genes in pancreatic cancer (Table 1), only a minority of patients (37% in one study) have mutations in all four genes, highlighting the genetic heterogeneity of the disease [11].
Mutations in these four driver genes were well described before the introduction of high throughput sequencing. However, multiple studies have now examined the whole exomes and whole genomes of large numbers of PDACs, and these studies have deepened our understanding of the pancreatic cancer in many ways [6,12]. First, these studies confirmed the importance of the four key driver genes (KRAS, P16/CDKN2A, TP53, and SMAD4) as the most frequently altered genes in PDAC. These studies also identified numerous other genes that are less commonly somatically mutated in PDACs—the average number of nonsynonymous genetic alterations ranged from 26 to 63 in the two key studies [6,12]. Several of these genes have known roles in tumorigenesis, including MLL3, TGFBR2, ATM, and ARID1A, and thus are likely to be drivers in pancreatic cancer in spite of their low mutation rate. However, many of the other infrequently mutated genes have no known role in cancer, and thus it is not possible to separate out driver mutations (which have a functional effect on tumorigenesis) from passenger mutations (which have no functional effect but instead accumulate randomly through repeated rounds of cell division within a tumor). Overall, these data highlight the heterogeneity at the gene level among different PDACs from different patients. However, although the individual genes altered are markedly heterogeneous, PDACs share many similarities when mutated genes are not considered at the individual gene level, but instead are considered as components of their larger biological pathways: there are several core processes and pathways that are genetically altered in the majority of carcinomas sequenced. These commonly altered pathways, which include KRAS signaling, DNA damage control, and cell adhesion, represent shared features of pancreatic tumorigenesis. In addition to these pathways with clear links to processes critical to tumor formation, one study also identified frequent mutations in genes in the axon guidance pathway as promising potential drivers in PDAC [12].
Somatic mutations have also been used to understand metastasis and model the time course of pancreatic tumorigenesis. Studies of somatic mutations in matched primary tumors and metastases did not identify any genetic alterations that were specific to metastasis [13]. Instead, the vast majority of the mutations identified in metastases were present in a subclonal population of the primary tumor. Modeling based on these data suggests a time period of approximately 15 years between the occurrence of the initiating mutation in PDAC and the acquisition of metastatic ability [13]. These data are encouraging, as they suggest a broad time window for early detection of PDAC while it is still curable by surgery [14].
Although the contribution of protein-altering somatic mutations is well documented in PDAC, it is likely that other genetic and epigenetic alterations also play a role in pancreatic tumorigenesis. This idea is supported by the observation that copy number alterations and promoter methylation affect known drivers of PDAC [15]. Many studies have identified copy number gains and losses as well as complex karyotypes in PDAC, and several studies have examined differential expression of micro-RNAs [16-21]. However, because of difficulties in determining the target genes and functional effects of these types of alterations, it is much more challenging to separate out driver and passenger alterations when analyzing large chromosomal alterations as well as differences in gene and microRNA expression.
There are several uncommon variants of PDAC—most of these are genetically similar to PDAC with respect to mutations in key driver genes (KRAS, P16/CDKN2A, TP53, and SMAD4). However, some variants also have unique genetic alterations. For example, in addition to mutations in previously known drivers, adenosquamous carcinoma has frequent somatic mutations in UPF1, which encodes a crucial component of the RNA degradation pathway of nonsense mediated decay [22,23]. The mutations in UPF1 cause aberrant splicing of the mRNA, eliminating essential protein domains and potentially conferring dominant negative activity [22]. Colloid carcinomas, which are characterized by large pools of stromal mucin and are associated with intestinal-type intraductal papillary mucinous neoplasms (IPMNs, see below), have somatic mutations in GNAS (which are not typically seen in PDAC), as well as KRAS and TP53 [24]. Medullary carcinomas exhibit a high prevalence of microsatellite instability [25,26], and undifferentiated carcinomas show frequent loss of E-cadherin protein expression [27] (Fig. 1C). Although no somatic mutations in CDH1 have been reported in undifferentiated carcinomas, promoter methylation has been reported and may explain the loss of E-cadherin protein expression [27].
FAMILIAL PANCREATIC CANCER GENES
A number of genes have been identified that increase the risk of developing PDAC when altered in the germline (Table 2). These genes are important to recognize for three reasons. First, the risk of developing pancreatic cancer can be quantified when the gene is known; this knowledge can help guide patient care and research screening. For example, germline mutations in PRSS1 cause chronic pancreatitis and greatly increase the risk of pancreatic cancer [28]. Some of these individuals, particularly those with a nonfunctioning pancreas from years of chronic pancreatitis, choose prophylactic pancreatectomy to reduce their risk of developing PDAC [29]. Second, the risk of developing extra-pancreatic malignancies can be quantified when the gene is known. For example, germline P16/CDN2KA mutations increase the risk of both PDAC and melanoma and lives can be saved by screening carriers of a germline P16/CDN2KA mutation and their biological relatives for melanoma [30,31]. Third, some of the germline mutations result in changes in the PDACs that are potentially therapeutically targetable. For example, germline BRCA2 mutations increase the risk of ovarian, breast and pancreatic cancer, and the pancreatic cancers that arise in patients with a germline BRCA2 mutation may be particularly sensitive to poly (ADP-ribose) polymerase (PARP) inhibitors and to DNA cross-linking agents such as cisplatinum [32,33].
PRECURSORS TO PANCREATIC DUCTAL ADENOCARCINOMA: PANCREATIC INTRAEPITHELIAL NEOPLASIAS
PDACs arise from histologically well-defined precursor lesions. The majority of PDACs arise from pancreatic intraepithelial neoplasia (PanIN), which are microscopic intraductal lesions and too small to be detected using currently available imaging technologies. However, a significant minority of PDACs arise from cystic precursors that can be identified with currently available imaging technologies—these cystic precursor lesions include IPMNs and mucinous cystic neoplasms (MCNs). Although some alterations are shared among these precursor lesions, there are genetic features that distinguish PanINs, IPMNs, and MCNs.
PanINs, which are categorized based on the morphological grade of dysplasia as low-grade (PanIN-1 and PanIN-2) and highgrade (PanIN-3), sequentially acquire the driver gene mutations that characterize PDAC. KRAS mutations are an early event in PanIN formation, occurring in >90% of even the lowest grade PanINs [34]. P16/CDKN2A mutations also occur in PanIN-1, though far less frequently than KRAS mutations, and the prevelance of P16/CDKN2A mutations increase with PanIN grade [34]. However, loss of Smad4 expression and TP53 alterations are late events. Smad4 loss has been reported only in high-grade PanIN (PanIN-3) and invasive carcinoma—the same pattern has been reported for TP53 alterations [35,36]. In addition to mutations in these driver genes, telomere shortening is a very early event in PanINs, occurring in approximately 90% of PanIN-1s [37]. More extensive genetic characterization of PanINs is complicated by several factors. First, these lesions are small and thus only yield enough DNA for specific targeted genetic analyses. Second, as these lesions are only visible microscopically, it is very difficult to prospectively identify them and harvest fresh tissue; thus, the vast majority of PanIN lesions are only available from formalinfixed paraffin-embedded tissue. The study of high-grade PanIN (PanIN-3) is particularly difficult. Because PDAC can invade into and then spread along the duct system, it can be impossible to distinguish true PanIN-3 from the intraductal spread of invasive carcinoma. In spite of these difficulties, one study has reported whole exome sequencing of PDAC and associated PanINs [38]. This study found that approximately two-thirds of the somatic mutations were shared between PanINs and associated invasive carcinoma, while 10% were unique to the invasive carcinoma and 25% were unique to the PanINs. This study confirms the genetic relationship between PanINs and invasive carcinoma but also suggests additional genetic evolution with clinical progression of the disease.
PRECURSORS TO PANCREATIC DUCTAL ADENOCARCINOMA: IPMNS AND MCNS
Cystic pancreatic cancer precursors, IPMNs and MCNs, can also give rise to PDAC. Both are cystic lesions but have distinct clinical and morphological differences. IPMNs are more common in the head of the pancreas and involve the pancreatic duct system. MCNs are more common in the tail of the pancreas, occur almost exclusively in women, do not involve the duct system, and have a characteristic ovarian-type stroma. Like PanINs, IPMNs and MCNs are also classified by their grade of dysplasia (low, intermediate, or high). In addition, IPMNs are also classified by the direction of differentiation of the neoplastic epithelium into gastric, intestinal, pancreatobiliary, and oncocytic subtypes.
IPMNs share mutations in the four key drivers of invasive PDAC (Table 1). KRAS mutations are prevalent in IPMNs and are present in up to 80% of tumors, occurring even in IPMNs with low-grade dysplasia [24,39]. Loss of p16 expression occurs in both noninvasive IPMNs and in IPMN associated invasive carcinomas, but it is much more common in invasive carcinomas (100% of invasive carcinomas vs 10% of noninvasive IPMNs in one study) [40]. Similarly, Smad4 loss is rare in noninvasive IPMNs but occurs in about one-third of IPMN-associated invasive carcinomas, identifying this as a late alteration in IPMN associated tumorigenesis [40,41]. TP53 mutations and protein overexpression are also limited to areas of high-grade dysplasia or invasive carcinoma in IPMNs [42-45].
In addition to these drivers shared with conventional PDAC, IPMNs also have alterations in genes not typically targeted in PDAC (Table 1). Approximately 60% of IPMNs have somatic mutations in the oncogenic hotspot of GNAS, which encodes a protein that couples transmembrane receptors to their downstream signaling proteins [24]. These mutations occur in all grades of dysplasia and are most prevalent in intestinal-type IPMNs [39,46]. Other genetic alterations that underlie the morphological differences between the IPMN subtypes remain unknown. In addition, inactivating mutations in RNF43, a gene coding for an ubiquitin ligase, occur in approximately 60% of IPMNs [47]. The mutations in GNAS and RNF43 are striking in their prevalence as well as mutation types. The mutations in GNAS cluster at a single oncogenic hotspot codon, unequivocally establishing this gene as an oncogene. In contrast, the mutations in RNF43 are striking enriched for nonsense and other inactivating mutations, indicating that this is a tumor suppressor gene. Other potential driver genes that are mutated in IPMN but not PDAC include PIK3CA (10% of IPMNs) and STK11/LKB1 (5% of IPMNs) [48,49]. Whole exome sequencing of IPMNs revealed an average of 26 nonsynonymous somatic mutations per tumor, far fewer than in PDAC [47].
Like IPMNs, MCNs share mutations in the four key drivers of PDAC—KRAS mutation is an early event, while tumor suppressor gene mutations (TP53, SMAD4, P16/CDKN2A) occur in MCNs with high-grade dysplasia or associated invasive carcinomas (Table 1). However, unlike IPMNs, GNAS mutations have not been identified in MCNs. Inactivating RNF43 mutations do occur in MCNs, but at a lower prevalence than in IPMNs [47]. Whole exome sequencing revealed an average of 16 nonsynonymous somatic mutations per IPMN [47].
OTHER CYSTIC NEOPLASMS
In addition to the cystic precursors to pancreatic cancer (IPMNs and MCNs), serous cystadenomas (SCAs) are clinically, pathologically, and genetically distinct from IPMNs and MCNs. SCAs are important to recognize because they are benign neoplasms with negligible risk of malignant behavior—these neoplasms are only resected if they cause symptoms due to mass effect or if there is clinical/radiologic concern for an IPMN. They are not precursors to invasive cancer. Whole exome sequencing of SCAs revealed an average of only 10 nonsynonymous somatic mutations per tumor, far fewer than in PDAC or its cystic precursors [47]. Only VHL was identified as a definite tumor suppressor gene in this tumor type, with somatic mutation in 50% and loss of heterozygosity in 90% of tumors (Table 1) [47]. SCAs occur as part of the clinical syndrome in von Hippel–Lindau syndrome, which results in clear cell neoplasms in various organs and is caused by germline mutations in VHL on chromosome 3—as such, the occurrence of somatic VHL mutations in sporadic SCAs is not surprising. Taken together, these data support the assertion that VHL is the major tumor suppressor responsible for the formation of both familial and sporadic SCAs. In addition, recent data show that vascular endothelial growth factor, a downstream target of VHL, is a marker of SCAs in pancreatic cyst fluid [50]. No mutations have been identified in PDAC or IPMN/MCN driver genes in SCAs. In addition to frequent loss of heterozygosity on chromosome 3p (the VHL locus), allelic loss of chromosome 10q has also been reported in 50% of cases and losses in several other chromosomes have been reported in multiple SCAs; however, no target genes for these losses has been identified [51].
SOLID-PSEUDOPAPILLARY NEOPLASMS
Solid-pseudopapillary neoplasms (SPNs) are uncommon pancreatic neoplasms with low malignant potential—although surgical resection is curative for most patients, metastases can rarely occur [52,53]. Activating somatic mutations in CTNNB1 occur in almost all cases, leading to abnormal nuclear accumulation of β-catenin which can be detected with immunohistochemistry (Fig. 1D) [54,55]. Importantly, this nuclear accumulation of β-catenin on immunohistochemical stains can be used to support the histological diagnosis of SPN. Whole exome sequencing of SPNs revealed shockingly few somatic mutations, an average of only three nonsynonymous somatic mutations per tumor [47]. Some tumors in this sequencing study had only one somatic mutation, but all SPNs had a somatic mutation in CTNNB1. This demonstrates a unique feature of SPNs: a remarkably low number of somatic mutations, lower than any tumor type described to date, including pediatric tumors. Additional studies are required to determine whether tumorigenesis in this tumor type is driven by other types of alterations, such as large chromosomal translocations or epigenetic events.
PANCREATIC NEUROENDOCRINE TUMORS
Pancreatic neuroendocrine tumors (PanNETs), also known as islet cell tumors, are a distinct pancreatic neoplasm with neuroendocrine differentiation, as demonstrated by morphology or immunohistochemistry. They are classified as “syndromic” if they express hormones that result in a clinical syndrome, and as “familial” if they arise in a patient with a genetic predisposition. PanNETs are graded based on their proliferation rate, as assessed by mitotic count or Ki67 labeling index [56]. Tumors with >20 mitoses per high power field or Ki67 index >20% are considered neuroendocrine carcinomas, and neuroendocrine carcinomas can be further subdivided into large cell and small cell carcinomas.
Whole exome sequencing of sporadic low-grade PanNETs identified an average of only 16 nonsynonymous somatic mutations per tumor [57]. Several driver genes unique to this tumor type were identified (Table 1). Somatic mutations in the MEN1 gene were identified in 45% of sporadic PanNETs [57,58]. Considering PanNETs are a major clinical manifestation of multiple endocrine neoplasia type 1 syndrome caused by germline mutation in MEN1, the presence of somatic mutations in this gene in sporadic PanNETs is not surprising. In addition, mutually exclusive somatic mutations in the chromatin remodeling genes ATRX and DAXX were identified in 45% of sporadic PanNETs [57]. These genes, which were enriched for inactivating mutations in PanNETs, function as part of a complex that is important for telomere maintenance—inactivation of ATRX and DAXX is associated with the telomerase independent telomere maintenance mechanism known as alternative lengthening of telomeres [59]. Studies in patients with MEN1 syndrome identified loss of ATRX and DAXX expression only in large PanNETs, suggesting that alterations in these genes are late events in pancreatic neuroendocrine tumorigenesis [60]. Somatic mutations in genes in the mammalian target of rapamycin cell signaling pathway (including PIK3CA, PTEN, and TSC2) occur in approximately 15% of sporadic PanNETs [57]. These alterations carry clinical significance, as drugs targeting this pathway have been developed for clinical use, and one such drug (Everolimus) has shown promise in PanNETs [61]. PanNETs lack frequent genetic alterations in major driver genes of PDAC—rare TP53 mutations have been reported, but no mutations in KRAS or P16/CDKN2A have been reported [57,62]. However, P16/CDKN2A promoter hypermethylation has been reported to be frequent in gastrinomas [62]. Although SMAD4 mutations were reported in one small study of nonfunctional PanNETs, this finding was not replicated in the whole exome sequencing study of PanNETs [57,63].
Intriguingly, high-grade neuroendocrine carcinomas are genetically distinct from low-grade PanNETs. Neuroendocrine carcinomas retain expression of ATRX and DAXX and instead have alterations of P53 and RB1 [64].
NEOPLASMS WITH ACINAR DIFFERENTIATION
Two rare neoplasms of the pancreas exhibit acinar differentiation: acinar cell carincoma and pancreatoblastoma. While acinar cell carcinoma typically occurs in adults, pancreatoblastoma is more common in children and can occur in patients with Beck-with-Wiedemann syndrome, a disorder associated with imprinting dysregulation of chromosome 11p leading to overgrowth of various organs and predisposition to embryonal tumors [65].
Acinar cell carcinomas are characterized by striking genomic instability at both the base pair and chromosomal level. In a recent whole exome sequencing study, each carcinoma had an average of 119 nonsynonymous somatic mutations, a larger number than any other pancreatic neoplasm [66]. Approximately 10% of acinar cell carincomas had microsatellite instability and a very large number of somatic mutations. Chromosomal instability was demonstrated by high fractional allelic losses in acinar cell carcinomas—sites of frequent loss included chromosome 11p (which was previously reported to be lost in a large proportion of acinar cell carcinomas) as well as tumor suppressor loci on chromosome 17p (TP53) and 18q (SMAD4) [67]. There was also a striking diversity in the genes altered by small somatic mutations in the whole exome sequencing—no single gene was altered in more than 30% of the acinar cell carcinomas sequenced. Mutations were identified in genes known to be drivers of PDAC (SMAD4 and TP53), cystic neoplasms (GNAS and RNF43), and PanNET (MEN1), as well as genes known to be drivers in other extrapancreatic tumor types. The mutations identified in APC and CTNNB1 confirmed previously reported findings of the importance of the Wnt signaling pathway in acinar cell carcinomas [67]. Other intriguing driver gene mutations identified in small subsets of acinar cell carcinomas included JAK1, BRAF, RB1, PTEN, ARID1A, MLL3, and BAP1. Importantly, these data suggest that therapies targeting mutations in JAK1 and BRAF could show promise in patients with acinar cell carcinoma. Many other genes were also mutated in 10%–20% of acinar cell carcinomas in the whole exome sequencing study, but further functional studies are required to determine the role of these mutations (if any) in tumorigenesis. Recent genetic studies have also identified frequent genomic rearrangements involving RAF genes (BRAF and RAF1) in almost one-fourth of acinar cell carcinomas [68]. Functional studies of the most frequently formed fusion gene SND1-BRAF show that it activates the mitogen-activated protein kinase pathway and confers sensitivity to MEK inhibition, pointing to a promising targeted therapy for a subset of patients with acinar cell carcinoma [68]. This study also confirmed the observation in the whole exome sequencing data of frequent mutations in DNA repair genes, suggesting potential utility for platinum-based chemotherapy or PARP inhibitors.
Loss of chromosome 11p is the most frequent genetic alteration in pancreatoblastomas, occurring in >80% of tumors [69]. Like acinar cell carcinomas, pancreatoblastomas also have mutations in Wnt pathway genes: activating CTNNB1 mutations are more common than inactivating APC mutations in pancreatoblastoma, although both can occur [69]. Pancreatoblastomas lack frequent mutations in the key PDAC driver genes. Although at least a subset of pancreatoblastomas contain large chromosomal alterations, these alterations are difficult to interpret since no target genes in these alterations have been identified [70-73]. Recent whole exome sequencing studies included two pancreatoblastomas, and although the number of tumors sequenced is small, a few observations can be made [66]. First, these tumors contained far fewer mutations than the acinar cell carcinomas—17 and 18 nonsynonymous somatic mutations per pancreatoblastoma. Second, both pancreatoblastomas had somatic mutations in CTNNB1, supporting previous studies demonstrating the importance of Wnt pathway alterations in this tumor type.
CONCLUSION
Although the genetic alterations underlying pancreatic neoplasms have been studied for decades, the recent development of high throughput sequencing has enabled systematic characterization of the genomes of all the major tumor types in the pancreas. These studies have confirmed previously identified drivers of pancreatic neoplasia as well as identified previously unknown genes that likely play a crucial role in tumorigenesis. Each type of pancreatic neoplasm has a unique genetic profile, and many potential therapeutic targets have been identified. As we enter into the era of genomic medicine, knowledge of the molecular alterations in pancreatic neoplasms will become a critical component of clinical care, as these alterations will likely form the basis for early detection strategies and targeted therapeutic approaches. This has the potential to lead to great improvements in the lives of patients with pancreatic neoplasms, as diagnostic and therapeutic approaches will be targeted to an individual patient’s tumor.
Notes
Conflicts of Interest
Dr. Hruban receives royalty payments from Myriad Genetics for the PALB2 invention.
Acknowledgements
Supported by the Kaya Tuncer Career Development Award in GI Cancer Prevention; the AGA-Bernard Lee Schwartz Research Scholar Award in Pancreatic Cancer; Sigma Beta Sorority; Joseph C. Monastra Fundation; Michael Rolfe Pancreatic Cancer Foundation; Sol Goldman Pancreatic Cancer Research Center; NIH grant CA62924.