- Open Access
Centromeric cohesion failure invokes a conserved choreography of chromosomal mis-segregations in pancreatic neuroendocrine tumor
Genome Medicine volume 12, Article number: 38 (2020)
Pancreatic neuroendocrine tumors (PANETs) are rare, slow growing cancers that often present with local and distant metastasis upon detection. PANETS contain distinct karyotypes, epigenetic dysregulation, and recurrent mutations in MEN1, ATRX, and DAXX (MAD+); however, the molecular basis of disease progression remains uncharacterized.
We evaluated associations between aneuploidy and the MAD+ mutational state of 532 PANETs from 11 published genomic studies and 19 new cases using a combination of exome, targeted panel, shallow WGS, or RNA-seq. We mapped the molecular timing of MAD+ PANET progression using cellular fractions corrected for inferred tumor content.
In 287 PANETs with mutational data, MAD+ tumors always exhibited a highly recurrent signature of loss of heterozygosity (LOH) and copy-number alterations affecting 11 chromosomes, typically followed by genome doubling upon metastasis. These LOH chromosomes substantially overlap with those that undergo non-random mis-segregation due to ectopic CENP-A localization to flanking centromeric regions in DAXX-depleted cell lines. Using expression data from 122 PANETs, we found decreased gene expression in the regions immediately adjacent to the centromere in MAD+ PANETs. Using 43 PANETs from AACR GENIE, we inferred this signature to be preceded by mutations in MEN1, ATRX, and DAXX. We conducted a meta-analysis on 226 PANETs from 8 CGH studies to show an association of this signature with metastatic incidence. Our study shows that MAD+ tumors are a genetically diverse and aggressive subtype of PANETs that display extensive chromosomal loss after MAD+ mutation, which is followed by genome doubling.
We propose an evolutionary model for a subset of aggressive PANETs that is initiated by mutation of MEN1, ATRX, and DAXX, resulting in defects in centromere cohesion from ectopic CENP-A deposition that leads to selective loss of chromosomes and the LOH phenotype seen in late-stage metastatic PANETs. These insights aid in disease risk stratification and nominate potential therapeutic vulnerabilities to treat this disease.
Pancreatic neuroendocrine tumors (PANETs) are rare neuroendocrine malignancies largely derived from pancreatic islet alpha-  and beta-cells . Approximately half of all PANETs are non-functional, defined as the absence of hormone secretion, and thus resulting in asymptomatic progression and late detection that typically co-occurs with liver metastasis . These tumors are characterized by mutations in chromatin modifiers MEN1, ATRX, and DAXX (MAD) (in 46, 18, and 31% of tumors, respectively) [3,4,5,6,7] and typically dichotomize into a genome that is either highly aneuploid or largely diploid with few copy-number variants (CNVs) [8,9,10,11,12,13,14,15,16,17]. Scarpa et al. defined a subtype of PANETs with a recurrent pattern of whole chromosomal loss (RPCL) in chromosomes 1, 2, 3, 6, 8, 10, 11, 15, 16, and 22 , while Stumpf et al. defined recurrent gains in the complementary set of chromosomes  suggesting a link via whole-genome duplication mechanisms or a technical difference in data normalization. The RPCL subtype is enriched for MAD mutations as well as an alternative lengthening of the telomere (ALT) phenotype  which indicates a potential functional link between the two. A recent characterization of PANETs highlights the role of epigenetic modifications into distinct subtype of this disease .
DAXX co-immunoprecipitates with both menin and ATRX via its C-terminal  and N-terminal regions, respectively . The menin-DAXX complex assembles on DNA where DAXX is unoccupied by histone variant H3.3/H4 and functions to enhance marks of H3K9me3 at the promoter of membrane metallo-endopeptidase (MME), a colorectal cancer oncogene [20,21,22]. The DAXX-ATRX complex participates in a functionally distinct pathway, catalyzing replication-independent deposition of the histone variant H3.3 at telomeric and pericentric heterochromatin regions [23,24,25]. Directly associated with H3.3 is the H3 variant, CENP-A , a histone protein that is responsible for assembling kinetochore proteins and dependent upon DAXX- [18, 25, 27, 28] and ATRX-mediated [24, 29,30,31,32] histone modifications for its endogenous localization [33, 34]. Dysregulation of DAXX induces mis-localization of CENP-A, resulting in chromosomal instability, neocentromere formation, and micronuclei formation, a common result of premature sister chromatid separation .
In our study of 532 PANETs, we sought to understand the pathogenesis of PANETs by examining the relationship between MAD mutations, chromosomal instability, cohesion, and CENP-A localization. We found that MAD mutations (MAD+) in PANET tumors were strongly predictive of a highly conserved pattern of loss of heterozygosity (LOH) and copy-number (CN) alterations across select chromosomes, typically followed by genome doubling in late-stage disease or metastatic disease. These patterns of chromosome mis-segregation are likely to stem from mis-localization of CENP-A in DAXX-deficient cells, resulting in merotelic attachments and premature sister chromatid separation via cohesion fatigue. Herein, we show that an aggressive subtype of PANETs follows a conserved progression of molecular events that originates from non-random chromosome mis-segregation and may suggest potential therapeutic targets to disrupt this choreography.
Our whole-exome sequencing (WES) cohort originated from 4 patients enrolled in the NET-SEQ study (ClinicalTrials.gov, NCT02586844) at the Princess Margaret Cancer Centre. Of the 7 patients registered in this study, 4 had histological or cytological diagnosis as well-differentiated pancreatic neuroendocrine tumors (PANETs) to be used for exploratory analysis. Our shallow whole-genome (sWGS) cohort was comprised of 15 NET samples provided by the Ontario Tumour Bank. Three sample types were processed: buffy coat blood cells, formalin-fixed paraffin-embedded (FFPE) tissues at time of diagnosis, and fresh-frozen core needle biopsies.
We sequenced DNA from the WES cohort to target a depth of 250× coverage in tumors and 50× coverage in normals. We also generated RNA sequencing (RNA-seq) libraries from these cases, which we sequenced using ~ 80 million reads. We sequenced DNA from the sWGS cohort to 0.34× mean coverage. Sequence data were aligned to the human reference genome sequence build hg19. Variant detection in exome data was performed using MuTect  and HaplotypeCaller , while copy-number profiles were called using VarScan2  and Sequenza . Loss of heterozygosity data was inferred from both DNA and RNA data by determining purity-adjusted allelic fractions. Gene-wise transcript abundances were quantified using the Cufflinks suite of tools . Pseudo allele-specific copy-number profiles were estimated from sWGS data using 500 kB bins tiled across the genome to count the number of reads and the number of heterozygous variants in each bin. To validate these copy-number calls, we paired this analysis with fluorescence in situ hybridization on complementary FFPE tissues.
CENP gene expression analysis
We analyzed 148 PANET gene expression profiles obtained from published microarray datasets: 99 generated by Sadanandam et al.  (GSE73338) and 49 from Chan et al.  (GSE117851). We compared expression patterns to a set of normal pancreatic islet cells from 57 non-diabetic and 20 diabetic donors (GSE41762) generated by Tang et al. . To approximate whether the CN signature was retained in MAD+ PANETs from these datasets, we first separated samples based on whether they carried MAD mutations and computed the z-score for gene expression against the MAD− PANETs on a per gene basis. Genes were mapped back to the human genome assembly hg19, and a loess regression with a 50% smoothing span was fitted to these values.
To calculate whether genes near the centromeres in LOH chromosomes are lower expressed than the rest of the chromosomal arm, we took the aforementioned gene expression z-scores and calculated the arm-level gene expression percentile and fractional distance to the centromere. Chromosomes were stratified into LOH and heterozygous chromosomes, and a loess regression was fit to the gene-level z-scores. To estimate regions of the chromosome arm that were repressed or elevated relative to the rest of the arm, we used an arm-level empirical cumulative density function to estimate the percentile of each gene.
Detection of monoallelic expression
To detect monoallelic expression of genes, we called all SNPs from RNA-seq data using HaplotypeCaller . We tested each gene containing 2 or more SNPs for MAE using a weighted t test. The allelic fractions of all SNPs in a gene, weighted by the number of reads supporting that SNP call, were compared to all SNP allelic fractions across the entire sample. By bootstrapping this calculation 1000 times per gene, we obtained the average z-statistic for each gene and compared it to a null distribution created using a similar test where the gene set is replaced with randomly selected SNPs.
Detection of parental skewing
SNPs from WES data of NET-001 tumor, matching blood DNA, and maternal DNA were estimated using HaplotypeCaller. All SNPs were divided into groups based on chromosomes and were then discretized into either homozygous (AF ≥ 0.8 or AF ≤ 0.2) or heterozygous (AF > 0.2 and AF < 0.8). We only focused on SNPs that were homozygous in the maternal DNA and heterozygous in the NET-001 germline DNA. For each LOH chromosome, we calculated the fraction of SNPs that were homozygous and either matched the maternal SNPs or did not (paternal), or were heterozygous.
CENP-A ChIP analysis
WIG files for the Nechemia dataset were downloaded from GEO:GSE111381 , while BigWIG files for the Nye dataset were downloaded from GEO:GSE120230 . Peaks were assigned to cytobands based on the hg19 reference genome. CENP-A peaks were summarized across a reference “merged peaks” representation defined by Nye et al. using two metrics: the max peak height for each merged peak or reads per kilobase of peaks per million mapped reads (RPKM).
For the Nye dataset, overlapping peaks between DAXX and control groups were compared using a t-statistic. To test for an elevated number of peaks in each cytoband, we calculated the Kolmorogov-Smirnov D-statistic by comparing the peaks found only in that cytoband against peaks found across the entire genome.
Alternative lengthening of telomere
Telomere lengths for all sWGS data were analyzed using Telomerecat . Samples were split between PANETs and GINETs, and a one-sided t test was done on the estimated telomere lengths. A one-sided F test was also conducted to calculate for difference in variance.
Meta-analysis of published datasets
Whole-exome sequencing of the BON-1 and QGP-1 PANET cell line from Vandamme and colleagues  was re-analyzed, and LOH segments were called based on allelic fractions (European Nucleotide Archive study ID: PRJEB8223). Copy-number profiles derived from CGH microarray data were obtained from data tables described in six publications [8,9,10,11,12,13,14,15,16] and transcribed into genomic coordinates (Additional file 1) by mapping to cytobands using the UCSC Table Browser hg19 cytoBandIdeo file (http://hgdownload.cse.ucsc.edu/goldenPath/hg19/database/cytoBandIdeo.txt.gz). Each copy ratio segment was assigned a value corresponding to the copy-status. Jaccard index values were calculated to measure the asymmetric binary concordance between any two copy-number profiles.
Molecular timing in project GENIE
Copy-number profiles and mutational data of PANETs from AACR’s project GENIE (v1.0.1) were downloaded from Sage Synapse (https://www.synapse.org/; synapse IDs: syn7851250, syn7851253, and syn7851246). In total, 43 PANET samples had both copy-number information and mutational information. The molecular timing of these samples was determined by estimating the tumor purity required for every possible copy-number profile to generate the observed tumor purity for all somatic mutations. The simplest copy-number profile that fits the constraints of pathologist purity ± 0.15 and copy-number constraints imposed by the relative copy-states of somatic mutations was used to infer molecular timing of the disease.
Additional methods and detailed version and parameter information are available in the Additional file 2.
Mutations in MEN1, ATRX, and DAXX are characteristic of chromosome mis-segregation errors
To characterize the molecular profiles of PANETs (Additional file 3: Fig. S1), we generated an exome and whole-genome dataset totaling 19 samples. The exome cohort consisted of whole-exome sequencing (WES) paired with RNA-seq across 4 patients, 6 samples (4 metastatic samples, 2 of which are complemented with their diagnostic tissue) (Additional file 4: Table S1a). Our whole-genome cohort consisted of 13 PANETs analyzed using shallow (0.3×) whole-genome sequencing (sWGS) and a 21-gene panel targeted deep sequencing, paired with FISH of 4 centromeric probes across 5 of the 13 samples (Additional file 4: Table S1b). Moreover, we included 165 PANET samples from version 5.0 of the publicly available AACR GENIE dataset  (Additional file 4: Table S1c).
By stratifying our cohort based on MAD status (MAD+ n = 18/19), we discovered a highly recurrent copy-number and LOH pattern that overlapped the RPCL pattern described by Scarpa et al.  in our exome and whole-genome cohorts (Fig. 1a, b). Three cases in our whole-genome cohort did not exhibit this LOH pattern; this pattern is undetected in NET-105 due to low tumor purity obscuring the signal (purity = 0.25, Additional file 4: Table S1b), while NET-130 and NET-131 did not show LOH in chromosomes 15, 16, 21, and 22. Although the inferred ploidy between samples differed, almost every case demonstrated LOH for chromosomes 1, 2, 3, 6, 8, 10, 11, 16, 21, and 22 and retained heterozygosity for chromosomes 4, 5, 7, 9, 12, 13, 14, 17, 19, and 20 (Fig 1d). Chromosomes 15 and 18 showed no consistent pattern of variation with heterozygosity and LOH occurring in equal proportions. Moreover, the LOH chromosomes were largely copy-neutral (i.e., diploid) while the heterozygous chromosomes showed copy-gain (Fig. 1e), mimicking the pattern observed by Stumpf et al. [8,9,10]. We validated the CN and LOH regions identified in our WES samples using the Affymetrix SNP 6.0 array (Additional file 3: Fig. S2a) and allelic skewing in RNA sequencing (Additional file 3: Fig. S2b). Additionally, we validated the CN inference in our 13 sWGS samples using fluorescent in situ hybridization of centromeric probes targeting LOH chromosomes 3 and 10, and zygosity-intact 7 and 17 (Additional file 4: Table S2a). To rule out the possibility of germline LOH, we confirmed that all patients had a diploid heterozygous genome in their germline DNA (Additional file 3: Fig. S3a).
Extended validation of CN and LOH signature
Next, we sought to expand our validation through inclusion of a larger, clinically derived cohort made available through the AACR GENIE consortium . While genome-wide zygosity calls were not available for the GENIE cohort, we were able to stratify the CN profiles of 165 PANETs samples into MAD+ (n = 99) and MAD− (n = 66) subgroups. Consistent with our genome-wide cohort, the targeted clinical panel sequencing data recapitulated the same pattern of losses and gains (Fig. 1c) and a near perfect overlap of MAD+ karyotypes. Chromosomes 8, 9, 15, and 21 had more than one prominent copy-states, suggesting more variable copy-number alterations of these chromosomes (Fig. 1f). PANETs in the GENIE cohort without MAD mutations were largely diploid with fewer recurrent gains and losses compared to those with MAD mutations (Fig. 1g).
We defined a CN signature by utilizing all available CN aberrations (CNA) data to calculate the co-occurence of copy-number states between chromosomes. By taking the copy-number states with the highest propensity (Fig. 1f), we flagged aberrations that are synchronous in their presentation from those that are random independent events (Additional file 3: Fig. S4a). Furthermore, we identified loss of chromosome 7 (Additional file 3: Fig. S4b) and gain of chromosomes 1, 6, and 16 (Additional file 3: Fig. S4d) as aberrations strongly antagonistic of our copy-number signature. Since the aberrations are largely chromosomal in size, we hypothesized that these CNA were likely a result of mis-segregation errors from merotelic events resulting in lagging chromosomes.
To evaluate whether PANET model systems accurately recapitulate these well-defined molecular signatures, we inferred genome-wide zygosity using publicly available WES data from two metastatic PANET cancer cell lines, BON-1 and QGP-1 . While both cell lines exhibited a high degree of aneuploidy and LOH, neither were MAD+ and the affected chromosomes differed dramatically between cell lines and when compared to the MAD+ PANETs in our meta-analysis (Additional file 3: Fig. S5). These results are in agreement with those of Boora et al. , suggesting that BON-1 and QGP-1 are genetically distinct from clinical samples of MAD− and MAD+ PANETs and should be used with caution in understanding PANET cancer biology or for pharmacological screening.
The MAD phenotype is associated with alternative lengthening of telomeres
To test whether MAD+ PANETs in our cohort exhibited ALT phenotype as reported by Jiao et al. [45, 48], we compared the overall length of telomeres between 13 PANETs and 10 MAD− GINETs (gastrointestinal neuroendocrine tumors) using sWGS. We observed longer telomere lengths in PANETs (p = 0.031; one-sided t test) as well as greater variation across samples (p < 0.001; one-sided F test) suggestive of an ALT phenotype (Additional file 3: Fig. S6). In the sWGS cohort, NET-129 lacked any MAD mutations and displayed telomere length concordant with the average length of GINET telomeres. The presence of ALT might suggest disruption of H3.3 incorporation, which we hypothesize is linked to the chromosomal mis-segregation pattern observed.
Mis-segregation errors are associated with DAXX-linked cohesion fatigue
The MAD+ CN signature may be a result of selective pressures favoring the retention and loss of certain chromosomes, or merely a stochastic event that always leads to the same karyotype. We first hypothesized that retention of chromosomes may be a consequence of selective pressures due to monoallelic expression (MAE). By examining 36 PANETs (23 MAD+, 13 MAD−) from our study and Chan et al.  that exhibited the LOH signature (Additional file 3: Fig. S7), we only found 11 genes that exhibited evidence of MAE scattered across 7 of the 8 chromosomes that were always heterozygous, leaving chromosome 14 unexplained (Additional file 3: Fig. S8a). Our second hypothesis was that there is a genetic predisposition to losing chromosomes from one parent. For NET-001, we flagged heterozygous single-nucleotide polymorphisms (SNPs) in germline DNA that were observed to be homozygous in the matched tumor and germline DNA collected from the patient’s mother. We determined that only 8 of the 11 chromosomes exhibiting LOH were of maternal origin (Chr1, 6, 8, 11, 13, 15, 21, and 22) and 3 were paternal (Chr2, 3, and 18) (Additional file 3: Fig. S8b), hence showing no statistical significance for favoring parental origin (p = 0.23, binomial test). Overall, we show that there are no robust selective factors driving our copy-number signature suggesting that another mechanism may be at play.
A recent publication by Worrall et al.  details a non-random chromosome mis-segregation pattern in SW403 colorectal cancer cells similar to our own (Fig. 2a; r = 0.50, point-biserial correlation) that is due to cohesion fatigue, resulting in premature sister chromatid separation and lagging chromosome formation. To explore the possibility of centromere dysfunction as the underlying mechanism driving formation of the MAD+ CN signature, we integrated our genomic data with two ChIP-seq datasets characterizing CENP-A binding locations. Nechemia-Arbely et al. defined centromeric CENP-A loads in HeLa cells throughout the cell cycle under endogenous and ectopic CENP-A expression  (Additional file 3: Fig. S9). Nye et al. characterized non-centromeric CENP-A binding in DAXX-depleted and DAXX-intact SW480 colorectal cancer cells  (Fig. 2b, c).
In HeLa cells, chromosomal CENP-A levels in the centromere were more similar between endogenous and ectopic expression conditions than between cell cycle phase, with chromosomes 2 and 9 containing the highest level of deposition and chromosomes 14, 19, and 21 the lowest (Additional file 3: Fig. S9). Meanwhile, in SW480 cells, there was an increase of ectopic CENP-A occupancy from 397 to 1124 kb when DAXX was depleted, the majority of this coverage occurring in chromosomes 8, 19, and 21 while losing coverage on chromosomes 1, 16, and 17 (Fig. 2d).
Using the single-cell sequencing (SCS) data from Worrall et al., we tested whether mis-segregation of specific chromosomes were related to centromeric CENP-A level and centromere size . We did not observe any correlation between CENP-A levels in centromeric regions and mis-segregation rates (rEndogenous_g1,g2 = − 0.005, − 0.09, rElevated_g1,g2,RC = − 0.11, 0.01, − 0.41), nor with the coverage or number of CENP-A binding sites acquired in DAXX-depleted conditions (rcoverage = 0.08, rcount = 0.12). However, there was a significant correlation between the frequency of mis-segregation per chromosome and the size of centromeres as well as flanking cytobands (hg19: rCEN = 0.58, p < 0.01; rflank = 0.51, p = 0.01; Fig. 2e). The features that most correlated with chromosomal mis-segregation rates were the mean CENP-A levels across all ectopic locations (rDAXX = 0.51, p = 0.02; rControl = 0.43, p = 0.06) and levels in regions close to the centromeric regions under both DAXX-depleted and control conditions (rDAXX = 0.56, p = 0.02; rControl = 0.58, p = 0.02; Fig. 2f). While the majority of ectopic CENP-A peaks localized to the flanking regions of centromere (Fig. 2b), we found that newly acquired CENP-A peaks in DAXX-depletion conditions were primarily localized to these flanking regions (Fig. 2c). These results suggest that DAXX-deficient tumors may mis-localize CENP-A to ectopic sites that flank the centromere, which could possibly seed the formation of neocentromeres and favor merotelic attachments of select chromosomes.
LOH chromosomes in MAD+ PANETs exhibit a gene-repressive environment directly adjacent to the centromere
It has been shown that proximity to chromocenters and pericentromeric regions results in gene repression . Therefore, to test whether PANETs have increased CENP-A loading in a similar fashion to DAXX-depleted colorectal cells, we tested whether there is a corresponding decrease of gene expression in regions proximal to the centromere in the LOH chromosomes. We examined the gene expression profiles of 122 PANETs with known MAD mutational status from two datasets: Sadanandam et al.  (28 MAD+, 47 MAD−) and Chan et al.  (30 MAD+, 17 MAD−). We first verified that the MAD+ samples recapitulate our previously defined copy-number signature (Fig. 3a, b). By normalizing the expression scores of each gene for all genes on the chromosome arms, we sought to identify whether there were regions on the chromosome that had lower or higher expression. We observed that there was a region, between 0 and 0.025 fractional distance of the centromere to the end of the chromosome arms, which had decreased gene expression in LOH chromosomes relative to heterozygous chromosomes (Fig. 3c, d). While this observation was noted in both datasets, it was noticeably absent in a dataset composed of 77 normal pancreatic islet cells  (Fig. 3c). The minor discrepancies between PANET datasets may reflect the inherent noise in RNA-seq data, the stratification of LOH and heterozygous chromosomes without genome data to validate, or the simplifying assumption that chromosomal arms only have single copy-state. As seen in the colorectal cell lines, the increased deposition of CENP-A in the pericentromeric region due to DAXX depletion may be linked to a corresponding decrease of gene expression in this region, suggesting an unseen mechanism for chromosomal mis-segregation pattern in PANETs.
Mutational events in MAD genes precede chromosomal mis-segregation
Given the progression of events proposed by our mis-segregation model, we assessed whether MAD mutations arose prior to LOH events. Hence, we developed a molecular timing analysis for the initial release of AACR GENIE dataset (Additional file 3: Fig. S1) . We obtained allelic fractions for clinical panel sequencing data from 43 mixed primary/metastasis samples (29 MAD+, 14 MAD−) that contain both copy-number and somatic mutation data in the GENIE v1.0 data freeze. Of the MAD+ population, 26/29 samples co-occur with the CN signature versus only 1/14 MAD− samples (Additional file 3: Fig. S10).
We next estimated the allele-specific CN profile of the GENIE cohort using the observed allelic fractions, CN log2 ratios, and pathologist-estimated tumor purities (± 0.15) (Supplementary Data). Of the MAD+ GENIE PANETs, 6/29 samples with low (< 30%) tumor cellularity were excluded from the analysis. The remaining 23 PANETs showed a strong tendency to adopt a CN model with cancer cell fraction of MAD mutations at 1.0, reinforcing the hypothesis that these mutations occur prior to LOH and genome doubling events (Fig. 4). We observed a significant enrichment of MEN1 and DAXX mutations prior to LOH and genome doubling events (Bonferroni adjusted p values: MEN1 = 0.00029, DAXX = 0.00011, binomial test) when using a cutoff of 0.85 cancer cell fraction. ATRX mutations reached significant enrichment at a cancer cell fraction cutoff of 0.63, which is expected as LOH on the X chromosome occurs infrequently in PANETs and may not always require LOH as a second hit due to X-inactivation. We observed that 35/39 MAD+ PANET samples follow a molecular timing model of MAD mutations prior to LOH (4/4 exome cohort, 10/12 whole-genome cohort, 21/23 GENIE cohort) (p = 3.4 × 10−7, binomial test). Collectively, our results provide evidence that acquisition of MEN1 and DAXX/ATRX mutations is an early event that leads to a genome-wide LOH event, likely through centromere fatigue from merotelic attachments.
Meta-analysis of copy-number profiles informs the molecular progression towards late-stage PANETs
PANET CN profiles have been extensively reported in the literature, each with their own reported chromosomes of significance. Scarpa et al.  examined 102 clinically sporadic PANETS and identified 4 CN groups: (G1) CN loss affecting > 50% of the genome, (G2) a diploid genome with little to no LOH or CN loss, (G3) CN gains affecting ~ 100% of the genome, and (G4) a mix of CN-LOH and CN gains. The absolute copy-number profiles of our exome cohort revealed 3 of the 4 groups described by Scarpa et al. The pancreatic diagnostic sample for patient NET-003 displays the G1 signature while the liver-metastasis sample displayed the G4 CN-LOH signature. Similarly, NET-009 presented with a similar transition of the G4 CN-LOH signature towards a G3 whole-genome gain (Fig. 1a), suggesting a mechanism of whole-genome duplications underlying PANET progression.
To further validate these groupings across independent cohorts, we analyzed previously published CGH datasets for the same signatures of absolute loss of LOH chromosomes (G1) or gain of retained chromosomes (G4). Due to the inability of CGH to detect CN-LOH and whole-genome gains, we anticipated tumors with the G3 profile to appear similar to G4 profiles defined by no aberrations in the LOH chromosomes with gains of the retained chromosomes. To compare the CGH copy-number data with our current study, we performed a meta-analysis of 226 NETs from 8 previous reports (Supplementary Data) [8,9,10,11,12,13,14,15,16]. By clustering the absolute copy-number profiles of our NETs and published datasets (Online methods), we demonstrated that tumors were divided into 5 clusters characterized by high and low fractions of genome-wide aneuploidy (Additional file 3: Fig. S11a). PANETs in our exome and whole-genome cohorts were mostly represented in cluster 1 which best represented G4 PANETs. Cluster 5 displayed loss of LOH chromosomes, suggesting that they best represent G1 PANETs. Cluster 4 was composed of 9 samples but contained karyotype that is reminiscent of profiles from the GENIE cohort, suggesting that this may reflect a normalization error rather than true biology. The remaining clusters 3 and 5 were largely diploid with few recurrent CN aberrations such as chromosome 11 loss.
As PANETs with increased chromosomal instability are characteristic of late-stage and more aggressive PANETs , we next sought to validate whether PANETs with the CN signature were in fact more aggressive. We separated samples based on copy-number profiles with high-chromosomal instability (high-CI) or low-chromosomal instability (low-CI) (Additional file 3: Fig. S11b). PANETs with high-CI were more likely to be metastatic (OR 4.35, 95% CI [1.99, 9.52]; p = 0.00; Cochrane’s Q) (Fig. 5), and the majority of the high-CI NETs were those found in clusters 1, 4, and 5 (Additional file 3: Fig. S11a); all tumors follow the proposed model of PANET CN progression. Thus, loss of the LOH chromosomes appears to be an initial step towards metastasis, reflecting a transient state due to the small sample size, resulting in a more stable and metastatic genome-doubled PANET (Fig. 6).
In our study, we integrate large publicly available datasets of PANETs to show a remarkably conserved MEN1- and DAXX/ATRX-driven metastatic disease progression. Across 306 PANETs with inferable copy-number profiles and MAD mutational status, we observed a well-defined pattern of LOH affecting select chromosomes following somatic mutations of MEN1 and DAXX or ATRX. By leveraging expression profiles for 122 of the 306 PANETs [1, 41], we show that this CN signature may be linked to perturbation of core kinetochore processes which would induce chromosomal mis-segregations. A recent publication from Worrall et al. suggests that there may be order in the timing of chromosomal mis-segregation stemming from merotelic attachments and cohesion fatigue . Strikingly, their mis-segregation fractions from SCS partially overlapped our LOH signature, suggesting a potential mechanism that we investigated using large CENP-A ChIP-seq datasets [35, 43]. Our results support the hypothesis that depletion of DAXX is associated with increased chromosome-specific ectopic CENP-A deposition, resulting in decreased gene expression , that correlates with mis-segregation frequencies. Finally, using the AACR GENIE dataset , we show that MEN1 and DAXX/ATRX mutational events all preceded the onset of chromosomal instability in clinical samples, subsequently resulting in LOH and whole-genome duplication to propagate chromosomal stability and increase tumor aggressiveness in 226 PANETs [8,9,10,11,12,13,14,15,16]. Overall, we define the molecular progression mechanisms for an aggressive subtype of PANETs which is also the first known observations to support the non-random chromosome mis-segregation theory  in primary clinical tumor specimens.
We observed that non-random mis-segregation of chromosomes in the SW480 colorectal cancer cell line  largely overlaps the LOH chromosomes in PANETs. This chromosome-specific overlap suggests that merotelic attachment and lagging chromosome formation is the common underlying mechanism guiding patterns of mis-segregation. Unequal chromosome and centromere sizes as well as levels of CENP-A can predispose chromosomes to merotelic attachment . The “placeholder theory” posits that H3.3 acts as a placeholder at centromeric domains during S phase, to be replaced by CENP-A during late G1 . H3.3 deposition at tandem repeat sites in centromeric and pericentromeric regions  is largely guided by the DAXX-ATRX complex. Disruption of H3.3 deposition is evident due to the presence of the alternative lengthening of telomere phenotype exclusively seen in the MAD+ PANETs [4, 24, 51, 52]. Using a dataset produced by Nye et al. that illustrates mis-localization of ectopic CENP-A in SW480 cells under DAXX-depleted conditions , we calculated a significant correlation with increased CENP-A deposition in regions flanking the centromere and mis-segregation fraction per chromosome. To translate these findings to PANETs, we observed a region immediately adjacent to the centromeres in LOH chromosomes in two published PANET datasets that had decreased gene expression relative to the rest of the chromosome arm. Thus, we propose that promiscuous CENP-A deposition to flanking centromeric regions is a factor that induces merotelic attachments, lagging chromosome formation, and the mis-segregation pattern proposed by Worrall et al.
The functional effects of aneuploidy are generally detrimental to cellular proliferation [53, 54] and can induce aneuploidy-associated stresses . Loss of an entire chromosome can have drastic effects, resulting in slowing of cell growth  but may allow advantages through loss of tumor suppressor genes (TSG) . Additionally, duplication of the remaining chromosomal region following a CN loss could harbor advantageous alterations, allowing cells to overcome the negative growth effects of chromosomal loss  or enhance for homozygous expression of preceding oncogenic mutations . A study by Taylor et al. illustrates this point in lung epithelial cells where they used CRISPR-Cas9 to induce loss of chromosome 3p. After several passages of slow growth, the cells acquired whole-chromosome duplication to overcome the negative growth effects incurred from the aberration . PANETs are well characterized as slow growing neoplasms that are clinically detectable only when they have metastasized ; the initial steps of disease progression described as near-global LOH could be an underlying mechanism for this slow growth. Acquired whole-genome duplication would be the tumor’s way to alleviate the negative growth effects of MAD-induced LOH, resulting in a more aggressive tumor that harbors loss of tumor suppressor genes but is not confounded by aneuploidy-associated stresses.
While the sample size of our in-house PANET cohort is small, we were able to leverage publicly available datasets allowing us to create a unifying model of disease progression to explain the remarkable consistency between karyotypes. However, due to the nature of this meta-analysis, we were unable to confidently validate our epigenetic dysregulation hypothesis due to the absence of publicly available ChIP-seq of H3.3 and CENP-A data in MAD+ and MAD− PANETs. Instead, we provide preliminary results leveraging work in SW480 and HeLa cells paired with evidence of a repressive gene environment proximal to centromeres in PANETs to provide compelling evidence to pursue further in vivo validation of this disease progression. Furthermore, we acknowledge that there are minor variations in copy-number and LOH profiles in the literature [16, 60, 61], but we hypothesize that these differences are a reflection of stochastic chromosomal instability events, mutational profiles allowing for more aneuploidy tolerance (e.g., TP53 mutation), or synthetic lethality which may alter which chromosomal losses are tolerated.
In our study, we observed a conserved trend of MEN1-, ATRX-, and DAXX-induced chromosome mis-segregation, leading to the characteristic karyotype of aggressive PANETs. These findings pave the way for functional validation studies to recreate the molecular progression of PANETs in model systems. Understanding the molecular basis of disease progression towards a more metastatic state has several benefits, specifically for risk stratification, treatment design for intermediate progression stages, and possibly even prophylactic treatment in at-risk individuals. As loss of heterozygosity appears as the molecular trigger for metastases, genome duplication acts as a mechanism to stabilize the genome.
Availability of data and materials
All somatic mutations and copy-number aberrations from the exome and whole-genome cohort can be found in Additional files 5 and 6, while raw sequencing data is controlled access at EGAS00001004239 (https://www.ebi.ac.uk/ega/studies/EGAS00001004239).
All CGH analyzed during this study are included in Additional file 1 as well as in the supplementary files of the original publications [8,9,10,11,12,13,14,15,16].
Publicly available gene expression profiles from PANETs were obtained from GEO datasets GSE73338 (doi: 10.1158/2159-8290.CD-15-0068)  and GSE117851 (doi: 10.1038/s41467-018-06498-2) . Gene expression profiles from normal pancreatic islet cells were obtained from the GEO dataset GSE41762 (doi: 10.1126/scitranslmed.3009934) .
Data for CENP-A ChIP-seq data were obtained from GEO datasets GSE111381 (doi: https://doi.org/10.1101/428557)  and GSE120230 (doi: https://doi.org/10.1371/journal.pone.0205948) .
Whole-exome sequencing data from PANET cell lines were obtained from European Nucleotide Archive study ID: PRJEB8223 (doi: https://doi.org/10.1530/JME-14-0304) .
The datasets generated during and/or analyzed during the current study are available in the AACR GENIE repository, http://www.cbioportal.org/genie/ (doi: 10.1158/2159-8290.CD-17-0151) .
American Association for Cancer Research
Alternative lengthening telomere
Core kinetochore genes
Copy-neutral loss of heterozygosity
Fluorescence in situ hybridization
Genomics Evidence Neoplasia Information Exchange
Gastrointestinal neuroendocrine tumor
Genotype-Tissue Expression project
Loss of heterozygosity
Mutant MEN1 and ATRX/DAXX
Wild-type MEN1 and ATRX/DAXX
Pancreatic neuroendocrine tumor
Recurrent pattern of whole chromosomal loss
Reads per kilobase of peaks per million mapped reads
Shallow whole-genome sequencing
Total copy number
tumor suppressor gene
Chan CS, Laddha SV, Lewis PW, Koletsky MS, Robzyk K, Da Silva E, et al. ATRX, DAXX or MEN1 mutant pancreatic neuroendocrine tumors are a distinct alpha-cell signature subgroup. Nat Commun. 2018;9(1):4158. https://doi.org/10.1038/s41467-018-06498-2.
Cejas P, Drier Y, Dreijerink KMA, Brosens LAA, Deshpande V, Epstein CB, et al. Enhancer signatures stratify and predict outcomes of non-functional pancreatic neuroendocrine tumors. Nat Med. 2019;25(8):1260–5.
Yao JC, Hassan M, Phan A, Dagohoy C, Leary C, Mares JE, et al. One hundred years after “carcinoid”: epidemiology of and prognostic factors for neuroendocrine tumors in 35,825 cases in the United States. J Clin Oncol. 2008;26(18):3063–72.
Jiao Y, Shi C, Edil BH, de Wilde RF, Klimstra DS, Maitra A, et al. DAXX/ATRX, MEN1, and mTOR pathway genes are frequently altered in pancreatic neuroendocrine tumors. Science. 2011;331(6021):1199–203.
Francis JM, Kiezun A, Ramos AH, Serra S, Pedamallu CS, Qian ZR, et al. Somatic mutation of CDKN1B in small intestine neuroendocrine tumors. Nat Genet. 2013;45(12):1483–6.
Banck MS, Kanwar R, Kulkarni AA, Boora GK, Metge F, Kipp BR, et al. The genomic landscape of small intestine neuroendocrine tumors. J Clin Invest. 2013;123(6):2502–8.
Cunningham JL, Díaz de Ståhl T, Sjöblom T, Westin G, Dumanski JP, Janson ET. Common pathogenetic mechanism involving human chromosome 18 in familial and sporadic ileal carcinoid tumors. Genes Chromosomes Cancer. 2011;50(2):82–94.
Stumpf E, Aalto Y, Höög A, Kjellman M, Otonkoski T, Knuutila S, Andersson LC. Chromosomal alterations in human pancreatic endocrine tumors. Genes Chromosomes Cancer. 2000;29(1):83–7. https://doi.org/10.1002/1098-2264(2000)9999:9999<::AID-GCC1011>3.0.CO;2-Z.
Terris B, Meddeb M, Marchio A, Danglot G, Fléjou JF, Belghiti J, et al. Comparative genomic hybridization analysis of sporadic neuroendocrine tumors of the digestive system. Genes Chromosomes Cancer. 1998;22(1):50–6. https://doi.org/10.1002/(SICI)1098-2264(199805)22:1<50::AID-GCC7>3.0.CO;2-6.
Tönnies H, Toliat MR, Ramel C, Pape UF, Neitzel H, Berger W, Wiedenmann B. Analysis of sporadic neuroendocrine tumours of the enteropancreatic system by comparative genomic hybridisation. Gut. 2001;48(4):536–41. https://doi.org/10.1136/gut.48.4.536.
Kulke MH, Freed E, Chiang DY, Philips J, Zahrieh D, Glickman JN, Shivdasani RA. High-resolution analysis of genetic alterations in small bowel carcinoid tumors reveals areas of recurrent amplification and loss. Genes Chromosomes Cancer. 2008;47(7):591–603.
Speel EJ, Scheidweiler AF, Zhao J, Matter C, Saremaslani P, Roth J, et al. Genetic evidence for early divergence of small functioning and nonfunctioning endocrine pancreatic tumors: gain of 9Q34 is an early event in insulinomas. Cancer Res. 2001;61(13):5186–92 PMID:11431358.
Zhao J, Moch H, Scheidweiler AF, Baer A, Schäffer AA, Speel EJ, et al. Genomic imbalances in the progression of endocrine pancreatic tumors. Genes Chromosomes Cancer. 2001;32(4):364–72.
Floridia G, Grilli G, Salvatore M, Pescucci C, Moore PS, Scarpa A, Taruscio D. Chromosomal alterations detected by comparative genomic hybridization in nonfunctioning endocrine pancreatic tumors. Cancer Genet Cytogenet. 2005;156(1):23–30. https://doi.org/10.1016/j.cancergencyto.2004.04.015.
Haugvik S-P, Vodák D, Haugom L, Hovig E, Gladhaug IP, Heim S, Micci F. Transcriptomic profiling of tumor aggressiveness in sporadic nonfunctioning pancreatic neuroendocrine neoplasms. Pancreas. 2016;45(8):1196–203.
Nagano Y, Kim DH, Zhang L, White JA, Yao JC, Hamilton SR, Rashid A. Allelic alterations in pancreatic endocrine tumors identified by genome-wide single nucleotide polymorphism analysis. Endocr Relat Cancer. 2007;14(2):483–92.
Scarpa A, Chang DK, Nones K, Corbo V, Patch A-M, Bailey P, et al. Whole-genome landscape of pancreatic neuroendocrine tumours. Nature. 2017;543(7643):65–71.
Feng Z, Wang L, Sun Y, Jiang Z, Domsic J, An C, et al. Menin and Daxx interact to suppress neuroendocrine tumors through epigenetic control of the membrane metallo-endopeptidase. Cancer Res. 2017;77(2):401–11.
Wang X, Zhao Y, Zhang J, Chen Y. Structural basis for DAXX interaction with ATRX. Protein Cell. 2017;8(10):767–71.
Drané P, Ouararhni K, Depaux A, Shuaib M, Hamiche A. The death-associated protein DAXX is a novel histone chaperone involved in the replication-independent deposition of H3.3. Genes Dev. 2010;24(12):1253–65.
Haase G, Gavert N, Brabletz T, Ben-Ze’ev A. A point mutation in the extracellular domain of L1 blocks its capacity to confer metastasis in colon cancer cells via CD10. Oncogene. 2017;36(11):1597–606.
Lewis PW, Elsaesser SJ, Noh K-M, Stadler SC, David Allis C. Daxx is an H3.3-specific histone chaperone and cooperates with ATRX in replication-independent chromatin assembly at telomeres. Proc Natl Acad Sci U S A. 2010;107(32):14075–80.
Watson LA, Goldberg H, Bérubé NG. Emerging roles of ATRX in cancer. Epigenomics. 2015;7(8):1365–78.
Morozov VM, Gavrilova EV, Ogryzko VV, Ishov AM. Dualistic function of Daxx at centromeric and pericentromeric heterochromatin in normal and stress conditions. Nucleus. 2012;3(3):276–85.
Dunleavy EM, Almouzni G, Karpen GH. H3. 3 is deposited at centromeres in S phase as a placeholder for newly assembled CENP-A in G1 phase. Nucleus. 2011;2(2):146–57.
Pluta AF, Earnshaw WC, Goldberg IG. Interphase-specific association of intrinsic centromere protein CENP-C with HDaxx, a death domain-binding protein implicated in Fas-mediated cell death. J Cell Sci. 1998;111(Pt 14):2029–41.
Goldberg AD, Banaszynski LA, Noh K-M, Lewis PW, Elsaesser SJ, Stadler S, et al. Distinct factors control histone variant H3.3 localization at specific genomic regions. Cell. 2010;140(5):678–91.
He Q, Kim H, Huang R, Lu W, Tang M, Shi F, et al. The Daxx/Atrx complex protects tandem repetitive elements during DNA hypomethylation by promoting H3K9 trimethylation. Cell Stem Cell. 2015;17(3):273–86.
Stimpson KM, Sullivan BA. Review of histone H3K4 methylation keeps centromeres open for business. EMBO J. 2011;30(2):233–4.
Baumann C, Viveiros MM, De La Fuente R. Loss of maternal ATRX results in centromere instability and aneuploidy in the mammalian oocyte and pre-implantation embryo. PLoS Genet. 2010;6(9):e1001137.
Ritchie K, Seah C, Moulin J, Isaac C, Dick F, Bérubé NG. Loss of ATRX leads to chromosome cohesion and congression defects. J Cell Biol. 2008;180(2):315–24.
Bergmann JH, Rodríguez MG, Martins NMC, Kimura H, Kelly DA, Masumoto H, et al. Epigenetic engineering shows H3K4me2 is required for HJURP targeting and CENP-A assembly on a synthetic human kinetochore. EMBO J. 2011;30(2):328–40.
González-Barrios R, Soto-Reyes E, Herrera LA. Assembling pieces of the centromere epigenetics puzzle. Epigenetics. 2012;7(1):3–13.
Shrestha RL, Ahn GS, Staples MI, Sathyan KM, Karpova TS, Foltz DR, Basrai MA. Mislocalization of centromeric histone H3 variant CENP-A contributes to chromosomal instability (CIN) in human cells. Oncotarget. 2017;8(29):46781–800.
Nye J, Sturgill D, Athwal R, Dalal Y. HJURP antagonizes CENP-A mislocalization driven by the H3.3 chaperones HIRA and DAXX. PloS One. 2018;13(10):e0205948. https://doi.org/10.1371/journal.pone.0205948.
Cibulskis K, Lawrence MS, Carter SL, Sivachenko A, Jaffe D, Sougnez C, et al. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat Biotechnol. 2013;31(3):213–9.
McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20(9):1297–303.
Koboldt DC, Zhang Q, Larson DE, Shen D, McLellan MD, Lin L, et al. VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 2012;22(3):568–76.
Favero F, Joshi T, Marquard AM, Birkbak NJ, Krzystanek M, Li Q, et al. Sequenza: allele-specific copy number and mutation profiles from tumor sequencing data. Ann Oncol. 2015;26(1):64–70.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.
Sadanandam A, Wullschleger S, Lyssiotis CA, Grötzinger C, Barbi S, Bersani S, et al. A cross-species analysis in pancreatic neuroendocrine tumors reveals molecular subtypes with distinctive clinical, metastatic, developmental, and metabolic characteristics. Cancer Discov. 2015;5(12):1296–313. https://doi.org/10.1158/2159-8290.CD-15-0068.
Tang Y, Axelsson AS, Spégel P, Andersson LE, Mulder H, Groop LC, et al. Genotype-based treatment of type 2 diabetes with an α2A-adrenergic receptor antagonist. Sci Transl Med. 2014;6(257):257ra139. https://doi.org/10.1126/scitranslmed.3009934.
Nechemia-Arbely Y, Miga KH, Shoshani O, Aslanian A, McMahon MA, Lee AY, et al. DNA replication-mediated error correction of ectopic CENP-A deposition maintains centromere identity. bioRxiv. 2018:428557. https://doi.org/10.1101/428557.
Farmery JHR, Smith ML, NIHR BioResource - Rare Diseases, Lynch AG. Telomerecat: a ploidy-agnostic method for estimating telomere length from whole genome sequencing data. Sci Rep. 2018;8(1):1300.
Vandamme T, Peeters M, Dogan F, Pauwels P, Van Assche E, Beyens M, et al. Whole-exome characterization of pancreatic neuroendocrine tumor cell lines BON-1 and QGP-1. J Mol Endocrinol. 2015;54(2):137–47. https://doi.org/10.1530/JME-14-0304.
AACR Project GENIE Consortium. AACR project GENIE: powering precision medicine through an international consortium. Cancer Discov. 2017;7(8):818–31. https://doi.org/10.1158/2159-8290.CD-17-0151.
Boora GK, Kanwar R, Kulkarni AA, Pleticha J, Ames M, Schroth G, et al. Exome-level comparison of primary well-differentiated neuroendocrine tumors and their cell lines. Cancer Genet. 2015;208(7–8):374–81.
Marinoni I, Kurrer AS, Vassella E, Dettmer M, Rudolph T, Banz V, et al. Loss of DAXX and ATRX are associated with chromosome instability and reduced survival of patients with pancreatic neuroendocrine tumors. Gastroenterology. 2014;146(2):453–460.e5.
Worrall JT, Tamura N, Mazzagatti A, Shaikh N, van Lingen T, Bakker B, et al. Non-random mis-segregation of human chromosomes. Cell Rep. 2018;23(11):3366–80.
Wijchers PJ, Geeven G, Eyres M, Bergsma AJ, Janssen M, Verstegen M, et al. Characterization and dynamics of pericentromere-associated domains in mice. Genome Res. 2015;25(7):958.
de Wilde RF, Heaphy CM, Maitra A, Meeker AK, Edil BH, Wolfgang CL, et al. Loss of ATRX or DAXX expression and concomitant acquisition of the alternative lengthening of telomeres phenotype are late events in a small subset of MEN-1 syndrome pancreatic neuroendocrine tumors. Modern Pathol. 2012;25(7):1033–9.
Heaphy CM, de Wilde RF, Jiao Y, Klein AP, Edil BH, Shi C, et al. Altered telomeres in tumors with ATRX and DAXX mutations. Science. 2011;333(6041):425.
Torres EM, Williams BR, Amon A. Aneuploidy: cells losing their balance. Genetics. 2008;179(2):737–46.
Thorburn RR, Gonzalez C, Brar GA, Christen S, Carlile TM, Ingolia NT, et al. Aneuploid yeast strains exhibit defects in cell growth and passage through START. Mol Biol Cell. 2013;24(9):1274–89.
Tang Y-C, Amon A. Gene copy-number alterations: a cost-benefit analysis. Cell. 2013;152(3):394–405.
Davoli T, Xu AW, Mengwasser KE, Sack LM, Yoon JC, Park PJ, Elledge SJ. Cumulative haploinsufficiency and triplosensitivity drive aneuploidy patterns and shape the cancer genome. Cell. 2013;155(4):948–62.
Taylor AM, Shih J, Ha G, Gao GF, Zhang X, Berger AC, et al. Genomic and functional approaches to understanding cancer aneuploidy. Cancer Cell. 2018;33(4):676–689.e3.
O’Keefe C, McDevitt MA, Maciejewski JP. Copy neutral loss of heterozygosity: a novel chromosomal lesion in myeloid malignancies. Blood. 2010;115(14):2731–9.
Cloyd JM, Poultsides GA. Non-functional neuroendocrine tumors of the pancreas: advances in diagnosis and management. World J Gastroenterol: WJG. 2015;21(32):9512–25.
Chung DC, Brown SB, Graeme-Cook F, Tillotson LG, Warshaw AL, Jensen RT, Arnold A. Localization of putative tumor suppressor loci by genome-wide allelotyping in human pancreatic endocrine tumors. Cancer Res. 1998;58(16):3706–11.
Rigaud G, Missiaglia E, Moore PS, Zamboni G, Falconi M, Talamini G, et al. High resolution allelotype of nonfunctional pancreatic endocrine tumors: identification of two molecular subgroups with clinical implications. Cancer Res. 2001;61(1):285–92.
We thank the staff of the Princess Margaret Genomics Centre (www.pmgenomics.ca, Neil Winegarden, Julissa Tsao, Nick Khuu, and Gurbaksh Basi) and the Bioinformatics and High-Performance Computing Core (Carl Virtanen, Zhibin Lu, and Natalie Stickle) for their expertise in generating the sequencing and microarray data used in this study. We thank Dr. Suzanne Kamel-Reid for helping to transfer the biological material to the Pugh lab in the initial steps of this study. We also thank Michael F. Berger of the Memorial Sloan Kettering Cancer Center for facilitating pathology estimates for samples within the AACR GENIE cohort. Biological materials were provided by the Ontario Tumour Bank, which is funded by the Ontario Institute for Cancer Research. This study was conducted with the support of the Ontario Institute for Cancer Research (OICR, PanCuRx Translational Research Initiative).
This project was funded by grants from the Cancer Research Society and the Carcinoid NeuroEndocrine Tumour Society Canada (#19341, TJP) and the Princess Margaret Cancer Centre Neuroendocrine Tumour Research Fund (LLS). TJP is supported by the Canada Research Chairs program. Infrastructure support was provided by the Princess Margaret Cancer Foundation; Canada Foundation for Innovation, Leaders Opportunity Fund, CFI 340 #32383; and Ontario Ministry of Research and Innovation, Ontario Research Fund Small Infrastructure Program (TJP). RQ is supported by a Medical Biophysics Excellence Ontario Student Opportunity Trust Fund Award from the Princess Margaret Cancer Centre, the Province of Ontario, and the University of Toronto. AS is supported by a Conquer Cancer Foundation ASCO Young Investigator Award.
Ethics approval and consent to participate
The University Health Network Research Ethics Board (UHN REB) approved this study (protocol number 13-6061-C) that includes obtaining informed consent from participants to collect, store, use, and share biological samples and data for the study. The UHN REB operates in compliance with the Tri-Council Policy Statement, ICH Guideline for Good Clinical Practice E6(R1), Ontario Personal Health Information Protection Act (2004), Part C Division 5 of the Food and Drug Regulations, and Part 4 of the Natural Health Products Regulations and the Medical Devices Regulations of Health Canada. This research conformed to the Declaration of Helsinki.
Consent for publication
Participants consented to sharing of de-identified data, having their name, address, and telephone stripped from any reports, publications, or presentations originating from this study. All participant information was anonymized.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1.
Details of CGH datasets, listing chromosome gains and losses as reported in the paper of origin.
Additional file 2.
Supplementary methods for the manuscript.
Additional file 3.
Supplementary figures for the manuscript.
Additional file 4.
Supplementary tables for the manuscript.
Additional file 5.
Somatic mutation calls from the exome and whole-genome cohort analyzed in this study.
Additional file 6.
Cellular fraction estimates of SNVs from 43 samples from the AACR GENIE v1.0 cohort as well as the exome and whole-genome datasets generated in this study.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Quevedo, R., Spreafico, A., Bruce, J. et al. Centromeric cohesion failure invokes a conserved choreography of chromosomal mis-segregations in pancreatic neuroendocrine tumor. Genome Med 12, 38 (2020). https://doi.org/10.1186/s13073-020-00730-9
- Exome sequencing
- Pancreatic neuroendocrine tumors
- Molecular timing
- Gene expression profiling
- Whole-genome sequencing
- Molecular cytogenetics
- Loss of heterozygosity
- Genetic instability
- Publicly available data