Frozen tissue coring and layered histological analysis improves cell type-specific proteogenomic characterization of pancreatic adenocarcinoma

BACKGROUND

Omics characterization of pancreatic adenocarcinoma tissue is complicated by the highly heterogeneous and mixed populations of cells. We evaluate the feasibility and potential benefit of using a coring method to enrich specific regions from bulk tissue and then perform proteogenomic analyses.

METHODS

We used the Biopsy Trifecta Extraction (BioTExt) technique to isolate cores of epithelial-enriched and stroma-enriched tissue from pancreatic tumor and adjacent tissue blocks. Histology was assessed at multiple depths throughout each core. DNA sequencing, RNA sequencing, and proteomics were performed on the cored and bulk tissue samples. Supervised and unsupervised analyses were performed based on integrated molecular and histology data.

RESULTS

Tissue cores had mixed cell composition at varying depths throughout. Average cell type percentages assessed by histology throughout the core were better associated with KRAS variant allele frequencies than standard histology assessment of the cut surface. Clustering based on serial histology data separated the cores into three groups with enrichment of neoplastic epithelium, stroma, and acinar cells, respectively. Using this classification, tumor overexpressed proteins identified in bulk tissue analysis were assigned into epithelial- or stroma-specific categories, which revealed novel epithelial-specific tumor overexpressed proteins.

CONCLUSIONS

Our study demonstrates the feasibility of multi-omics data generation from tissue cores, the necessity of interval H&E stains in serial histology sections, and the utility of coring to improve analysis over bulk tissue data.

Detection of Pancreatic Ductal Adenocarcinoma-Associated Proteins in Serum

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancer types, partly because it is frequently identified at an advanced stage, when surgery is no longer feasible. Therefore, early detection using minimally invasive methods such as blood tests may improve outcomes. However, studies to discover molecular signatures for the early detection of PDAC using blood tests have only been marginally successful. In the current study, a quantitative glycoproteomic approach via data-independent acquisition mass spectrometry was utilized to detect glycoproteins in 29 patient-matched PDAC tissues and sera. A total of 892 N-linked glycopeptides originating from 141 glycoproteins had PDAC-associated changes beyond normal variation. We further evaluated the specificity of these serum-detectable glycoproteins by comparing their abundance in 53 independent PDAC patient sera and 65 cancer-free controls. The PDAC tissue-associated glycoproteins we have identified represent an inventory of serum-detectable PDAC-associated glycoproteins as candidate biomarkers that can be potentially used for the detection of PDAC using blood tests.

Clinically impactful metabolic subtypes of pancreatic ductal adenocarcinoma (PDAC)

Background: Pancreatic ductal adenocarcinoma (PDAC) is a lethal disease characterized by a diverse tumor microenvironment. The heterogeneous cellular composition of PDAC makes it challenging to study molecular features of tumor cells using extracts from bulk tumor. The metabolic features in tumor cells from clinical samples are poorly understood, and their impact on clinical outcomes are unknown. Our objective was to identify the metabolic features in the tumor compartment that are most clinically impactful. Methods: A computational deconvolution approach using the DeMixT algorithm was applied to bulk RNASeq data from The Cancer Genome Atlas to determine the proportion of each gene's expression that was attributable to the tumor compartment. A machine learning algorithm designed to identify features most closely associated with survival outcomes was used to identify the most clinically impactful metabolic genes. Results: Two metabolic subtypes (M1 and M2) were identified, based on the pattern of expression of the 26 most important metabolic genes. The M2 phenotype had a significantly worse survival, which was replicated in three external PDAC cohorts. This PDAC subtype was characterized by net glycogen catabolism, accelerated glycolysis, and increased proliferation and cellular migration. Single cell data demonstrated substantial intercellular heterogeneity in the metabolic features that typified this aggressive phenotype. Conclusion: By focusing on features within the tumor compartment, two novel and clinically impactful metabolic subtypes of PDAC were identified. Our study emphasizes the challenges of defining tumor phenotypes in the face of the significant intratumoral heterogeneity that typifies PDAC. Further studies are required to understand the microenvironmental factors that drive the appearance of the metabolic features characteristic of the aggressive M2 PDAC phenotype.

A clinically useful and biologically informative genomic classifier for papillary thyroid cancer

Clinical management of papillary thyroid cancer depends on estimations of prognosis. Standard care, which relies on prognostication based on clinicopathologic features, is inaccurate. We applied a machine learning algorithm (HighLifeR) to 502 cases annotated by The Cancer Genome Atlas Project to derive an accurate molecular prognostic classifier. Unsupervised analysis of the 82 genes that were most closely associated with recurrence after surgery enabled the identification of three unique molecular subtypes. One subtype had a high recurrence rate, an immunosuppressed microenvironment, and enrichment of the EZH2-HOTAIR pathway. Two other unique molecular subtypes with a lower rate of recurrence were identified, including one subtype with a paucity of BRAFV600E mutations and a high rate of RAS mutations. The genomic risk classifier, in addition to tumor size and lymph node status, enabled effective prognostication that outperformed the American Thyroid Association clinical risk stratification. The genomic classifier we derived can potentially be applied preoperatively to direct clinical decision-making. Distinct biological features of molecular subtypes also have implications regarding sensitivity to radioactive iodine, EZH2 inhibitors, and immune checkpoint inhibitors.

Clinical and metabolomic characterization of Brivanib-Induced hypertension in metastatic colorectal cancer

BACKGROUND

Trials of tyrosine kinase inhibitors (TKI) have not demonstrated dramatic benefits in advanced colorectal cancer (CRC), and this may be a function of poor patient selection. TKI-induced hypertension is reportedly a surrogate marker for treatment benefit for some tumor types. Our objective was to determine whether hypertension was associated with benefit in the context of CRC treatment, and also to gain insight on the pathogenesis of TKI-induced hypertension by monitoring associated changes in the circulating metabolome.

PATIENTS AND METHODS

Clinical data were acquired from clinical trial patients with metastatic CRC randomized to cetuximab ± the TKI brivanib (N = 750). Outcomes were evaluated as a function of treatment-induced hypertension. For metabolomic studies, plasma samples were taken at baseline, as well as at 1, 4, and 12 weeks after treatment initiation. Samples were submitted to gas chromatography-mass spectrometry to identify treatment-related metabolomic changes associated with TKI-induced hypertension, compared to pre-treatment baseline. A model based on changes in metabolite concentrations was generated using orthogonal partial least squares discriminant analysis (OPLS-DA).

RESULTS

In the brivanib treated group, 95 patients had treatment-related hypertension within 12 weeks of initiating treatment. TKI-induced hypertension was not associated with a significantly higher response rate, nor was it associated with improved progression-free or overall survival. In metabolomic studies, 386 metabolites were identified. There were 29 metabolites that changed with treatment and distinguished patients with and without TKI-induced hypertension. The OPLS-DA model for brivanib-induced hypertension was significant and robust (R2 Y score = 0.89, Q2 Y score = 0.70, CV-ANOVA = 2.01 e-7). Notable metabolomic features previously reported in pre-eclampsia and associated with vasoconstriction were found.

CONCLUSION

TKI-induced hypertension was not associated with clinical benefit in metastatic CRC. We have identified changes in the metabolome that are associated with the development of worsening brivanib-induced hypertension that may be useful in future efforts of characterizing this toxicity.

Abstract 2049: Clinically and biologically distinct molecular subtypes of lung squamous cell carcinoma (LUSC)

Background: Lung cancer is the leading cause of cancer death worldwide, with LUSC accounting for a third of cases. Stage I and II LUSC is typically treated with resection, and decisions related to adjuvant chemotherapy are related to recurrence risk, which is a function of stage. However, disease stage does not fully capture tumor biology, which may be a more important determinant of recurrence risk. Our objective was to identify biologically significant molecular subgroups of LUSC that more accurately reflect recurrence risk. Methods: Transcriptomic data for resected stage I/II LUSC were obtained from The Cancer Genome Atlas (TCGA). A training set consisting of patients who underwent resection without adjuvant chemotherapy or radiotherapy was used for discovery (N=161). A proprietary machine learning algorithm (HighLifeR™) was used to identify the genes most associated with disease-free survival (DFS). Molecular subgroups were identified by unsupervised clustering of prognostic genes. Functional differences between molecular subgroups were identified by gene set enrichment analysis (GSEA) and Ingenuity Pathway Analysis (IPA). Results: HighLifeR™ identified 60 highly prognostic genes. Two molecular subgroups were identified with significantly different median DFS: a high-risk group with 59.300 months and hazard ratio of 11.551, and a low-risk group that did not reach 50% DFS (log-rank p = 1.96 x 10-5). These subgroups did not differ in age, sex, race, smoking status, or overall stage (p > 0.05). Importantly, univariate and multivariate Cox regression analysis showed that the molecular subgroups outperformed clinical stage in predicting DFS. The molecular subgroups were biologically distinct. On GSEA, the high-risk group was significantly enriched in genes involved in mitotic spindle assembly (NES=1.94, p<0.001); TGFα signaling and PI3K/AKT/mTOR signaling were modestly enriched. IPA pathway analysis demonstrated significant enrichment of numerous pathways linked to neuronal growth and signaling (eg: CREB signaling, S100 signaling, pathways in myelination and synaptogenesis), as well as TGFα signaling and AMPK signaling (all p<0.001). Conclusions: Our prognostic transcriptomic signature identified two biologically distinct molecular subgroups of LUSC. Molecular subgroup classification was more predictive of DFS in stage I and II LUSC than any other clinical or pathological variable. If further validation confirms this, then this biomarker may form the basis of a diagnostic test that helps inform which patients could be considered for adjuvant chemotherapy.

Citation Format: Ashar Siddiqui, Cynthia Stretch, Farshad Farshidfar, Oliver F. Bathe. Clinically and biologically distinct molecular subtypes of lung squamous cell carcinoma (LUSC) [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 2049.

Dietary EPA+DHA Mitigate Hepatic Toxicity and Modify the Oxylipin Profile in an Animal Model of Colorectal Cancer Treated with Chemotherapy

Irinotecan (CPT-11) and 5-fluorouracil (5-FU) are commonly used to treat metastatic colorectal cancer, but chemotherapy-associated steatosis/steatohepatitis (CASSH) frequently accompanies their use. The objective of this study was to determine effect of CPT-11+5-FU on liver toxicity, liver oxylipins, and cytokines, and to explore whether these alterations could be modified by dietary eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) in the form of fish oil (EPA+DHA). Tumor-bearing animals were administered CPT-11+5-FU and maintained on a control diet or a diet containing EPA+DHA (2.3 g/100 g). Livers were collected one week after chemotherapy for the analysis of oxylipins, cytokines, and markers of liver pathology (oxidized glutathione, GSSH; 4-hydroxynonenal, 4-HNE, and type-I collagen fiber). Dietary EPA+DHA prevented the chemotherapy-induced increases in liver GSSH (p < 0.011) and 4-HNE (p < 0.006). Compared with the tumor-bearing animals, ten oxylipins were altered (three/ten n-6 oxylipins were elevated while seven/ten n-3 oxylipins were reduced) following chemotherapy. Reductions in the n-3 fatty-acid-derived oxylipins that were evident following chemotherapy were restored by dietary EPA+DHA. Liver TNF-α, IL-6 and IL-10 were elevated (p < 0.05) following chemotherapy; dietary EPA+DHA reduced IL-6 (p = 0.09) and eotaxin (p = 0.007) levels. Chemotherapy-induced liver injury results in distinct alterations in oxylipins and cytokines, and dietary EPA+DHA attenuates these pathophysiological effects.

Abstract B053: Targeting SMURF1 with low-dose proteasome inhibitors in pancreatic cancer organoids

Global initiatives focused on whole genome and transcriptome analysis (WGTA) of metastatic pancreatic ductal adenocarcinoma (PDAC) tumor cohorts are driving progress in understanding the clinical impact of PDAC heterogeneity. This knowledge is key to continue the discovery of subsets of patients who could benefit from biomarker-informed targeted therapy and expand the currently limited PDAC treatment options. To identify clinically actionable subtypes in metastatic PDAC (mPDAC), our team analyzes prospectively collected WGTA data from patients enrolled in the PanGen trial (NCT02869802), in parallel with analysis of a set of patient-derived organoids (PDOs). In a cohort of 69 sequenced tumors, four cases (5.8%) had an amplification of chromosome 7q22, which included copy gains of transcriptional co-factor TRRAP, drug metabolizing cytochrome P450 genes CYP3A4 and CYP3A5, and SMURF1. SMURF1 is a ubiquitin-protein ligase that regulates TGFβ receptor signaling, in part via its interaction with SMAD7, and has been implicated in the epithelial-to-mesenchymal program and tumor invasiveness in PDAC. The proteasome inhibitor (PI) bortezomib has been shown to attenuate SMURF1 levels. We investigated whether SMURF1 signaling axis is associated with sensitivity to PIs in mPDAC, by analyzing the cytotoxic effects of the PIs bortezomib, carfilzomib and ixazomib in PDOs. Eight PDOs were treated with PIs at concentrations between 0.1pM-1mM. The presence of live and dead cells was quantified using the IN Cell Analyzer, and cell toxicity was analyzed using GRtoxic metrics (grcalculator.org). The PIs bind catalytic subunits of the 26S proteasome, primarily proteasome subunit beta-5, encoded by PSMB5. The cytotoxic activity of all three PIs negatively correlated with tumor biopsy PSMB5 expression (p&lt;0.05), suggesting that PDOs established from tumors with low expression of PSMB5 are more sensitive to PI-induced cell death. There was a positive correlation between ixazomib cytotoxicity and SMURF1 expression (p&lt;0.05). PanGen patients with high (&gt;75th percentile) tumor SMURF1 expression showed shorter overall survival (OS) compared to the rest of the cohort (median OS 9 vs. 14 months, respectively; p=0.004). Immunohistochemistry analysis of a tumor tissue microarray comprising 175 resected PDAC cases detected an association of concomitant low SMURF1 and high SMAD7 protein levels with adjuvant therapy response (median OS: 40 vs. 12 months with no adjuvant treatment, p=0.002), suggesting that attenuation of SMURF1 in the presence of SMAD7 may improve chemotherapy response in early disease stages. In summary, we present data indicating increased sensitivity to PIs in a subset of tumors with low PSMB5 and high SMURF1 expression, and highlight the translational utility of the investigation of genomic and clinically annotated pre-clinical models in PDAC predictive biomarker discovery.

Citation Format: Andrew Metcalfe, Joanna M. Karasinska, James T. Topham, Steve E. Kalloger, Hassan Ali, Dawn Ashforth, Marco A. Marra, Janessa Laskin, Patricia A. Tang, Rachel Goodwin, Oliver F. Bathe, Daniel J. Renouf, David F. Schaeffer. Targeting SMURF1 with low-dose proteasome inhibitors in pancreatic cancer organoids [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr B053.

Abstract B066: Elucidating the role of insulin receptor isoform expression in metastatic pancreatic ductal adenocarcinoma

Introduction: While type 2 diabetes mellitus (T2DM) has been shown to be both a risk factor and consequence of pancreatic ductal adenocarcinoma (PDAC), the bidirectional relationship between T2DM and PDAC remains to be fully understood. Due to spatial proximity of insulin-producing beta cells, PDAC tumor cells are presumably exposed to a highly concentrated insulin microenvironment, which may contribute to the dysregulation of the insulin receptor gene INSR. There are two isoforms of insulin receptor: IR-A and IR-B. IR-A activity has been associated with oncogenic function and is upregulated in several cancer types. This leads to questions regarding the role of INSR and more specifically, how the ratio of the two isoforms may potentiate tumor aggressiveness in PDAC. Results: Using data from whole-genome and RNA sequencing of metastatic PDAC (mPDAC) from both the PanGen (n=70; NCT01855477) and Personalized Oncogenomics (n=22; NCT02155621) trials, we identified somatic copy loss of INSR in 39/92 (42%) patient tumors, with the majority of such cases (35/39; 90%) showing heterozygous copy loss. There was no statistical association between INSR copy status and PDAC transcriptomic subtypes. INSR expression was lower in basal-like versus classical subtype tumors (p=1.1e-4), and INSR expression was further attenuated in basal-like tumors with heterozygous copy loss of INSR (p=0.0041). The ratio of IR-A:IR-B expression was heterogenous across samples, and heightened IR-A:IR-B ratio was significantly (p&lt;0.05) associated with expression of genes linked to the PI3K-Akt signaling axis (CCDC88A, THEM4) and glucose metabolism (HK1, G6PC, PKLR) in basal-like tumors. Interestingly,HK1 and IGF1R were significantly upregulated, independent of IR-A:IR-B ratio, among basal-like (p 9.1e-9 and p=1.4e-5, respectively) compared to classical subtype tumors. Conclusion: These data indicate that heightened IR-A:IR-B ratio is associated with the expression of key pathways converging on PI3K signalling and glucose metabolism in mPDAC, while highlighting that several IR-A:IR-B-associated gene correlations are unique to basal-like tumors. Such findings warrant further investigation of relative INSR isoform expression in mPDAC and pre-clinical models to elucidate the role of the oncogenic isoform, IR-A, in driving tumor aggressiveness in basal-like subtype tumors.

Citation Format: Lan V. Tao, James T. Topham, Joanna M. Karasinska, Erica S. Tsang, Andrew Metcalfe, Hassan Ali, Dawn Ashforth, Rachel Goodwin, Patricia A. Tang, Oliver F. Bathe, Janessa Laskin, Marco Marra, Daniel J. Renouf, David F. Schaeffer. Elucidating the role of insulin receptor isoform expression in metastatic pancreatic ductal adenocarcinoma [abstract]. In: Proceedings of the AACR Special Conference on Pancreatic Cancer; 2022 Sep 13-16; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2022;82(22 Suppl):Abstract nr B066.

Neoplastic cell enrichment of tumor tissues using coring and laser microdissection for proteomic and genomic analyses of pancreatic ductal adenocarcinoma

BACKGROUND

The identification of differentially expressed tumor-associated proteins and genomic alterations driving neoplasia is critical in the development of clinical assays to detect cancers and forms the foundation for understanding cancer biology. One of the challenges in the analysis of pancreatic ductal adenocarcinoma (PDAC) is the low neoplastic cellularity and heterogeneous composition of bulk tumors. To enrich neoplastic cells from bulk tumor tissue, coring, and laser microdissection (LMD) sampling techniques have been employed. In this study, we assessed the protein and KRAS mutation changes associated with samples obtained by these enrichment techniques and evaluated the fraction of neoplastic cells in PDAC for proteomic and genomic analyses.

METHODS

Three fresh frozen PDAC tumors and their tumor-matched normal adjacent tissues (NATs) were obtained from three sampling techniques using bulk, coring, and LMD; and analyzed by TMT-based quantitative proteomics. The protein profiles and characterizations of differentially expressed proteins in three sampling groups were determined. These three PDACs and samples of five additional PDACs obtained by the same three sampling techniques were also subjected to genomic analysis to characterize KRAS mutations.

RESULTS

The neoplastic cellularity of eight PDACs ranged from less than 10% to over 80% based on morphological review. Distinctive proteomic patterns and abundances of certain tumor-associated proteins were revealed when comparing the tumors and NATs by different sampling techniques. Coring and bulk tissues had comparable proteome profiles, while LMD samples had the most distinct proteome composition compared to bulk tissues. Further genomic analysis of bulk, cored, or LMD samples demonstrated that KRAS mutations were significantly enriched in LMD samples while coring was less effective in enriching for KRAS mutations when bulk tissues contained a relatively low neoplastic cellularity.

CONCLUSIONS

In addition to bulk tissues, samples from LMD and coring techniques can be used for proteogenomic studies. The greatest enrichment of neoplastic cellularity is obtained with the LMD technique.

Integrative analysis of KRAS wildtype metastatic pancreatic ductal adenocarcinoma reveals mutation and expression-based similarities to cholangiocarcinoma

Oncogenic KRAS mutations are absent in approximately 10% of patients with metastatic pancreatic ductal adenocarcinoma (mPDAC) and may represent a subgroup of mPDAC with therapeutic options beyond standard-of-care cytotoxic chemotherapy. While distinct gene fusions have been implicated in KRAS wildtype mPDAC, information regarding other types of mutations remain limited, and gene expression patterns associated with KRAS wildtype mPDAC have not been reported. Here, we leverage sequencing data from the PanGen trial to perform comprehensive characterization of the molecular landscape of KRAS wildtype mPDAC and reveal increased frequency of chr1q amplification encompassing transcription factors PROX1 and NR5A2. By leveraging data from colorectal adenocarcinoma and cholangiocarcinoma samples, we highlight similarities between cholangiocarcinoma and KRAS wildtype mPDAC involving both mutation and expression-based signatures and validate these findings using an independent dataset. These data further establish KRAS wildtype mPDAC as a unique molecular entity, with therapeutic opportunities extending beyond gene fusion events.

KRAS wildtype metastatic pancreatic ductal adenocarcinoma (mPDAC) could represent a distinct molecular entity from other PDACs. Here, the authors analyse KRAS wildtype mPDAC tumours using genomics and transcriptomics and find molecular similarities with cholangiocarcinomas.

Integrative analysis of KRAS-wildtype pancreatic ductal adenocarcinoma reveals unique similarities to extrahepatic cholangiocarcinoma

587

Background: Oncogenic driver mutations in KRAS represent a hallmark genomic event in approximately 90% of pancreatic adenocarcinoma (PDAC). For the remaining 10% of patients with KRAS wildtype (wt) PDAC, distinct driver mutations have been described, but their transcriptional landscape has not been reported. Here, we leverage sequencing data from the PanGen trial to provide a comprehensive characterization of advanced KRASwt PDAC. Methods: 63 patients with advanced PDAC received whole genome and transcriptome sequencing prior to treatment for metastatic disease as part of the PanGen trial (NCT02869802). Clinical features, somatic mutation data and gene expression patterns were compared between KRASwt and mutant groups. PDAC samples were contrasted with 77 other metastatic carcinoma (colorectal and cholangiocarcinoma) samples from the Personalized OncoGenomics trial (NCT0215562). KRAS wt-associated genes were further investigated using 3 additional PDAC cohorts (COMPASS NCT02750657, TCGA, and ICGC). Results: 9 of 63 (14%) samples were KRASwt, with an earlier median age at diagnosis (51.4 vs. 60.9 years; p=0.03). Clinical features, including diabetes, family history of malignancy, and location of primary tumor, were comparable. CA 19-9 at baseline was lower in the KRASwt group, with median 58 vs. 4900 U/mL in the KRAS-mutant group ( p=0.03). Patients with KRASwt PDAC showed increased overall survival in univariable ( p=0.0024) and multivariable ( p=0.0089) analyses. 6 of 9 (67%) KRASwt tumors had fusions involving NRG1 (n = 3), FGFR2 (n = 1), BRAF (n = 1) or NTRK2 (n = 1), while known actionable fusions were not observed in KRAS mutant patients. KRASwt tumors showed increased expression of genes associated with cholangiocytes and grouped with cholangiocarcinoma samples in unsupervised clustering analysis. Validation using three independent PDAC cohorts revealed a core set of 70 KRAS wt-associated genes that converge on keratinization, ion transport, and hormone metabolism pathways. Conclusions: Patients with KRASwt PDAC show potentially targetable molecular traits with actionable fusions. We also highlight novel mutation and expression-based similarities between KRASwt PDAC and cholangiocarcinoma samples. Recurrent dysregulation of genes involved in cellular structure and metastasis provide impetus for further investigation into the developmental trajectory and potential therapeutic vulnerabilities of KRASwt PDAC.

Preoperative Single-dose Methylprednisolone Prevents Surgical Site Infections After Major Liver Resection: A Randomized Controlled Trial

OBJECTIVE

The primary aim of this study was to evaluate the efficacy of a single preoperative dose of methylprednisolone for preventing postoperative complications after major liver resections.

SUMMARY BACKGROUND DATA

Hepatic resections are associated with a significant acute systemic inflammatory response. This effect subsequently correlates with postoperative morbidity, mortality, and length of recovery. Multiple small trials have proposed that the administration of glucocorticoids may modulate this effect.

METHODS

This study was a parallel, dual-arm, double-blind randomized controlled trial. Adult patients undergoing elective major hepatic resection (≥3 segments) at a quaternary care institution were included (2013-2019). Patients were randomly assigned to receive a single preoperative 500 mg dose of methylprednisolone versus placebo. The main outcome measure was postoperative complications after liver resection, within 90 days of the index operation. Standard statistical methodology was employed (P < 0.05 = significant).

RESULTS

A total of 151 patients who underwent a major hepatic resection were randomized (mean age = 62.8 years; 57% male; body-mass-index = 27.9). No significant differences were identified between the intervention and control groups (age, sex, body-mass-index, preoperative comorbidities, hepatic function, ASA class, portal vein embolization rate) (P > 0.05). Underlying hepatic diagnoses included colorectal liver metastases (69%), hepatocellular carcinoma (18%), noncolorectal liver metastases (7%), and intrahepatic cholangiocarcinoma (6%). There was a significant reduction in the overall incidence of postoperative complications in the methylprednisolone group (31.2% vs 47.3%; P = 0.042). Patients in the glucocorticoid group also displayed less frequent organ space surgical site infections (6.5% vs 17.6%; P = 0.036), as well as a shorter length of hospital stay (8.9 vs 12.5 days; P = 0.015). Postoperative serum bilirubin and prothrombin timeinternational normalized ratio (PT-INR) levels were also lower in the steroid group (P = 0.03 and 0.04, respectively). Multivariate analysis did not identify any additional significant modifying factor relationships (estimated blood loss, duration of surgery, hepatic vascular occlusion (rate or duration), portal vein embolization, drain use, etc) (P > 0.05).

CONCLUSIONS

A single preoperative dose of methylprednisolone significantly reduces the length of hospital stay, postoperative serum bilirubin, and PT-INR, as well as infectious and overall complications following major hepatectomy.

Early detection of treatment futility in patients with metastatic colorectal cancer

PURPOSE

Chemotherapy options for treating CRC have rapidly expanded in recent years, and few have predictive biomarkers. Oncologists are challenged with evidence-based selection of treatments, and response is evaluated retrospectively based on serial imaging beginning after 2-3 months. As a result, cumulative toxicities may appear in patients who will not benefit. Early recognition of non-benefit would reduce cumulative toxicities. Our objective was to determine treatment-related changes in the circulating metabolome corresponding to treatment futility.

METHODS

Metabolomic studies were performed on serial plasma samples from patients with CRC in a randomized controlled trial of cetuximab vs. cetuximab + brivanib (N = 188). GC-MS quantified named 94 metabolites and concentrations were evaluated at baseline, Weeks 1, 4 and 12 after treatment initiation. In a discovery cohort (N = 68), a model distinguishing changes in metabolites associated with radiographic disease progression and response was generated using OPLS-DA. A cohort of 120 patients was used for validation of the model.

RESULTS

By one week after treatment, a stable model of 21 metabolites could distinguish between progression and partial response (R2Y = 0.859; Q2Y = 0.605; P = 5e-4). In the validation cohort, patients with the biomarker had a significantly shorter OS (P < 0.0001). In a separate cohort of patients with HCC on axitinib, appearance of the biomarker also signified a shorter PFS (1.7 months vs. 9.2 months, P = 0.001).

CONCLUSION

We have identified changes in the metabolome that appear within 1 week of starting treatment associated with treatment futility. The novel approach described is applicable to future efforts in developing a biomarker for early assessment of treatment efficacy.

Abstract IA-003: Proteogenomic characterizations of pancreatic ductal adenocarcinoma

Pancreatic cancer is one of the deadliest cancers and the five-year survival rate is less than 10%. Pancreatic ductal adenocarcinoma (PDAC) represents more than 90% of all pancreatic malignancies, and is responsible for the majority of pancreatic cancer-related deaths. Towards understanding the underlying molecular alterations that drive PDAC oncogenesis and identify therapeutic targets for personalized treatments, we comprehensively characterized 140 pancreatic cancers and 67 normal adjacent tissues. To ensure robust, downstream analyses, tumor neoplastic cellularity was assessed via multiple, orthogonal strategies using molecular features, and verified via pathological estimation of tumor cellularity based on histological review to select tumors with sufficient tumor cellularity. We also included the analysis of 9 normal pancreatic ductal tissues. Proteomic, phosphoproteomic, and glycoproteomic analyses were used to characterize proteins and their modifications. In addition, whole genome sequencing, whole exome sequencing, methylation, RNA-seq, and miRNA-seq were performed on the same tissues to facilitate an integrated proteogenomic analysis and determine the impact of genomic alterations on protein expression, signaling pathways, and post-translational modifications. These characterizations revealed functional impacts of genomic and epigenomic alterations on proteins and protein modifications, delineated PDAC cell microenvironment compositions and the immune signatures for immunotherapy, also uncovered putative kinase inhibitors that could be tested for therapy. This integrated proteogenomic characterization of PDAC will serve as a valuable resource for the community, paving the way for early detection and identification of novel therapeutic targets.

Citation Format: Liwei Cao, Chen Huang, Daniel Cui Zhou, Yingwei Hu, Mamie Lih, Sara R. Savage, Karsten Krug, David J. Clark, Michael Schnaubelt, Lijun Chen, Felipe da Veiga Leprevost, Rodrigo Vargas Eguez, Alexey I. Nesvizhskii, D.R. Mani, Gilbert S. Omenn, Emily S. Boja, Mehdi Mesri, Ana I. Robles, Henry Rodriguez, Oliver F. Bathe, Daniel W. Chan, Ralph H. Hruban, Li Ding, Bing Zhang, Hui Zhang, Clinical Proteomic Tumor Analysis Consortium. Proteogenomic characterizations of pancreatic ductal adenocarcinoma [abstract]. In: Proceedings of the AACR Virtual Special Conference on Pancreatic Cancer; 2021 Sep 29-30. Philadelphia (PA): AACR; Cancer Res 2021;81(22 Suppl):Abstract nr IA-003.

Report from the 21st Annual Western Canadian Gastrointestinal Cancer Consensus Conference; Calgary, Alberta; 20-21 September 2019

The 21st annual Western Canadian Gastrointestinal Cancer Consensus Conference (WCGCCC) was held in Calgary, Alberta, 20-21 September 2019. The WCGCCC is an interactive multi-disciplinary conference attended by health care professionals from across Western Canada (British Columbia, Alberta, Saskatchewan, and Manitoba) involved in the care of patients with gastrointestinal cancer. Surgical, medical, and radiation oncologists, pathologists, radiologists, and allied health care professionals such as dietitians and nurses participated in presentation and discussion sessions to develop the recommendations presented here. This consensus statement addresses current issues in the management of hepato-pancreato-biliary (HPB) cancers.

Proteogenomic characterization of pancreatic ductal adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is a highly aggressive cancer with poor patient survival. Toward understanding the underlying molecular alterations that drive PDAC oncogenesis, we conducted comprehensive proteogenomic analysis of 140 pancreatic cancers, 67 normal adjacent tissues, and 9 normal pancreatic ductal tissues. Proteomic, phosphoproteomic, and glycoproteomic analyses were used to characterize proteins and their modifications. In addition, whole-genome sequencing, whole-exome sequencing, methylation, RNA sequencing (RNA-seq), and microRNA sequencing (miRNA-seq) were performed on the same tissues to facilitate an integrated proteogenomic analysis and determine the impact of genomic alterations on protein expression, signaling pathways, and post-translational modifications. To ensure robust downstream analyses, tumor neoplastic cellularity was assessed via multiple orthogonal strategies using molecular features and verified via pathological estimation of tumor cellularity based on histological review. This integrated proteogenomic characterization of PDAC will serve as a valuable resource for the community, paving the way for early detection and identification of novel therapeutic targets.

Evaluation of the Implementation of Multiple Enhanced Recovery After Surgery Pathways Across a Provincial Health Care System in Alberta, Canada

IMPORTANCE

Engaging multidisciplinary care teams in surgical practice is important for the improvement of surgical outcomes.

OBJECTIVE

To evaluate the association of multiple Enhanced Recovery After Surgery (ERAS) pathways with ERAS guideline adherence and outcomes.

DESIGN, SETTING, AND PARTICIPANTS

This quality improvement study compared a pre-ERAS cohort (2013-2017) with a post-ERAS cohort (2014-2018). All patients were from Alberta Health Services in Alberta, Canada, and had available ERAS and up to 1-year postsurgery administrative data. Data collected included age, sex, body mass index, tobacco and alcohol use, diabetes, comorbidity index, and surgical characteristics. Data analysis was performed from May 7, 2020, to February 1, 2021.

INTERVENTIONS

Implementation of 5 ERAS pathways (colorectal, liver, pancreas, gynecologic oncology, and radical cystectomy) across 9 sites.

MAIN OUTCOMES AND MEASURES

Adherence to ERAS guidelines was measured by the percentage of patients whose care met the common ERAS pathway care element criteria. Surgical procedures were grouped by complexity; complications were classified by severity. Outcome measures for the pre-post-ERAS cohorts included length of stay (LOS), readmission, complications, and mortality.

RESULTS

A total of 7757 patients participated in the study, including 984 in the pre-ERAS cohort (median [interquartile range] age, 62 [53-71] years; 526 [53.5%] female) and 6773 in the post-ERAS cohort (median [interquartile range] age, 62 [53-71] years; 3470 [51.2%] male). In the total cohort, care-element adherence improved from 52% to 76% (P < .001), no significant differences were found in serious complications (from 6.2% to 4.9%; P = .08) or 30-day mortality (from 0.71% to 0.93%; P = .50), 1-year mortality decreased from 7.1% to 4.6% (P < .001), mean (SD) LOS decreased from 9.4 (7.0) to 7.8 (5.0) days (P < .001), and 30-day readmission rates were unchanged (from 13.4% to 11.7%; P = .12). After adjustment for patient characteristics, the LOS mean difference decreased 0.71 days (95% CI, -1.13 to -0.29 days; P < .001), with no significant differences in adjusted 30-day readmission (-3.5%; 95% CI, -22.7% to 20.4%; P = .75), serious complications (1.3%; 95% CI, -26.2% to 39.0%; P = .94), or mortality (30-day mortality: 42% [95% CI, -35.4% to 212.3%]; P = .38; 1-year mortality: 8% [95% CI, -20.5% to 46.8%]; P = .62). The adjusted 1-year readmission rate was -15.6% (95% CI, -27.7% to -1.5%; P = .03) in favor of ERAS, and readmission LOS was shorter by 1.7 days (95% CI, -3.3 to -0.1 days; P = .04).

CONCLUSIONS AND RELEVANCE

The results of this quality improvement study suggest that implementation of ERAS across multiple pathways may improve health care practitioner adherence to ERAS guidelines, LOS, and readmission rates at a system level.

Abstract 17: Proteogenomic characterization of pancreatic ductal adenocarcinoma

Pancreatic Ductal Adenocarcinoma (PDAC) is a highly aggressive cancer with 5-year survival rate less than 10%. Although numerous efforts have been made to characterize the functional driver genes of PDAC, the impacts of genomic alterations on protein modifications and molecular mechanism of this cancer type are still not well defined. To this end, Investigators from The Clinical Proteomic Tumor Analysis Consortium (CPTAC) conducted the first comprehensive characterization of 140 pancreatic cancers, 67 normal adjacent tissues, and 9 normal pancreatic ductal tissues using whole genome sequencing, whole exome sequencing, methylation, RNA-seq, miRNA-seq, proteomics, phosphoproteomics, and glycoproteomics. To address the inherent low neoplastic cellularity of pancreatic cancer, multiple orthogonal strategies using molecular features and histology imaging were deployed to deconvolute the cellularity and identify 105 tumors with sufficient neoplastic purity for downstream analyses. Proteomics, phosphoproteomics, and glycoproteomics analyses by mass spectrometry identified and quantified 11,662 proteins, 51,469 phosphosites, and 34,024 glycopeptides, respectively. Genomic data revealed that 97% of the cancers harbored KRAS alterations in our cohort. In addition, we found that somatic mutations and copy number alterations could impact gene and protein expression, as well as protein modifications such as phosphorylation and glycosylation. Glycoproteomic analyses uncovered tumor-associated alterations in glycoprotein biosynthesis and glycoproteins as potential targets for diagnosis or therapeutic intervention. Over-expressed kinase substrates and their corresponding kinases were identified with integration of global proteomic and phosphoproteomic measurements, thus identifying potential therapeutic targets. The results of our molecular and cellular subtyping revealed a small group of immune-hot subtype tumors that could benefit from immunotherapy, as well as the underlying mechanisms associated with major groups of immune-cold subtypes, including endothelial cell remodeling, glycolysis, and cell junction dysregulation. Non-negative matrix factorization (NMF)-based proteogenomics subtyping revealed two clusters with strong prognostic relevance. Finally, we identified proteins and glycoproteins overexpressed in early stage pancreatic cancer that may serve as candidates for early detection. This comprehensive proteogenomic characterization of PDAC provides a valuable resource for uncovering the molecular mechanisms of this cancer and thus paves the way for discovery of novel early detection and therapeutic targets.

Citation Format: Liwei Cao, Chen Huang, Daniel Cui Zhou, Oliver F. Bathe, Daniel W. Chan, Ralph H. Hruban, Li Ding, Bing Zhang, Hui Zhang, Clinical Proteomic Tumor Analysis Consortium Investigators. Proteogenomic characterization of pancreatic ductal adenocarcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 17.

Prognostic value of early changes in CT-measured body composition in patients receiving chemotherapy for unresectable pancreatic cancer

OBJECTIVES

Skeletal muscle mass is a prognostic factor in pancreatic ductal adenocarcinoma (PDAC). However, it remains unclear whether changes in body composition provide an incremental prognostic value to established risk factors, especially the Response Evaluation Criteria in Solid Tumors version 1.1 (RECISTv1.1). The aim of this study was to determine the prognostic value of CT-quantified body composition changes in patients with unresectable PDAC starting chemotherapy.

METHODS

We retrospectively evaluated 105 patients with unresectable (locally advanced or metastatic) PDAC treated with FOLFIRINOX (n = 64) or gemcitabine-based (n = 41) first-line chemotherapy within a multicenter prospective trial. Changes (Δ) in skeletal muscle index (SMI), subcutaneous (SATI), and visceral adipose tissue index (VATI) between pre-chemotherapy and first follow-up CT were assessed. Cox regression models and covariate-adjusted survival curves were used to identify predictors of overall survival (OS).

RESULTS

At multivariable analysis, adjusting for RECISTv1.1-response at first follow-up, ΔSMI was prognostic for OS with a hazard ratio (HR) of 1.2 (95% CI: 1.08-1.33, p = 0.001). No significant association with OS was observed for ΔSATI (HR: 1, 95% CI: 0.97-1.04, p = 0.88) and ΔVATI (HR: 1.01, 95% CI: 0.99-1.04, p = 0.33). At an optimal cutoff of 2.8 cm²/m² per 30 days, the median survival of patients with high versus low ΔSMI was 143 versus 233 days (p < 0.001).

CONCLUSIONS

Patients with a lower rate of skeletal muscle loss at first follow-up demonstrated improved survival for unresectable PDAC, regardless of their RECISTv1.1-category. Assessing ΔSMI at the first follow-up CT may be useful for prognostication, in addition to routine radiological assessment.

KEY POINTS

• In patients with unresectable pancreatic ductal adenocarcinoma, change of skeletal muscle index (ΔSMI) in the early phase of chemotherapy is prognostic for overall survival, even after adjusting for Response Evaluation Criteria in Solid Tumors version 1.1 (RECISTv1.1) assessment at first follow-up. • Changes in adipose tissue compartments at first follow-up demonstrated no significant association with overall survival. • Integrating ΔSMI into routine radiological assessment may improve prognostic stratification and impact treatment decision-making at the first follow-up.

Germline genetic contribution to the immune landscape of cancer

Understanding the contribution of the host's genetic background to cancer immunity may lead to improved stratification for immunotherapy and to the identification of novel therapeutic targets. We investigated the effect of common and rare germline variants on 139 well-defined immune traits in ∼9000 cancer patients enrolled in TCGA. High heritability was observed for estimates of NK cell and T cell subset infiltration and for interferon signaling. Common variants of IFIH1, TMEM173 (STING1), and TMEM108 were associated with differential interferon signaling and variants mapping to RBL1 correlated with T cell subset abundance. Pathogenic or likely pathogenic variants in BRCA1 and in genes involved in telomere stabilization and Wnt-β-catenin also acted as immune modulators. Our findings provide evidence for the impact of germline genetics on the composition and functional orientation of the tumor immune microenvironment. The curated datasets, variants, and genes identified provide a resource toward further understanding of tumor-immune interactions.

A Population-Based Retrospective Study of Biliary Tract Cancers in Alberta, Canada

Background: Biliary tract cancers (BTC) are uncommon malignancies and are underrepresented in the literature. Methods: We performed a retrospective population-based review of adult patients with biopsy-confirmed BTC in Alberta from 2000 to 2015. Demographic data, risk factors, symptoms, treatment, and staging data were collected and analyzed. Survival analyses were completed. Results: A total of 1604 patients were included in our study, of which 766 (47.8%) were male. The median age at diagnosis was 68 (range 19–99). There were 374 (23.3%) patients with resectable tumors at diagnosis versus 597 (37.2%) with unresectable tumors. Of the patients, 380 (21.5%) received chemotherapy (CT) and 81 (5.0%) underwent radiation therapy. There was a clear trend with worsening stage and performance status associated with shorter median overall survival (OS). Ampulla of Vater tumors had the best median OS (25.69 months), while intrahepatic bile duct cancers had the worst (5.78 months). First-line palliative CT regimens included gemcitabine+cisplatin (OS 14.98 months (mo), n = 212), single agent gemcitabine (OS 12.42 mo, n = 22), capecitabine (OS 8.12 mo, n = 8), and capecitabine+gemcitabine (OS 6.93 mo, n = 13). Patients with advanced or metastatic disease who received first-line gemcitabine+cisplatin had a median OS of 11.8 months (n = 119). Conclusion: BTCs have poor survival. Worse outcomes occur in higher stage and poorer Eastern Cooperative Oncology Group (ECOG) performance status patients across all tumor subtypes. Tumor resectability at diagnosis was associated with better OS. Our study supports the use of gemcitabine+cisplatin as a combination first-line palliative CT, as patients treated in Alberta have a comparable OS to that reported in the ABC-02 phase III study.

Cancer-Associated Muscle Wasting-Candidate Mechanisms and Molecular Pathways

Excessive muscle loss is commonly observed in cancer patients and its association with poor prognosis has been well-established. Cancer-associated sarcopenia differs from age-related wasting in that it is not responsive to nutritional intervention and exercise. This is related to its unique pathogenesis, a result of diverse and interconnected mechanisms including inflammation, disordered metabolism, proteolysis and autophagy. There is a growing body of evidence that suggests that the tumor is the driver of muscle wasting by its elaboration of mediators that influence each of these pro-sarcopenic pathways. In this review, evidence for these tumor-derived factors and putative mechanisms for inducing muscle wasting will be reviewed. Potential targets for future research and therapeutic interventions will also be reviewed.

Optimal treatments for hepato-pancreato-biliary trauma in severely injured patients: a narrative scoping review

SUMMARY

Hepato-pancreato-biliary (HPB) injuries can be extremely challenging to manage. This scoping review (8438 citations) offers a number of recommendations. If diagnosis and therapy are rapid, patients with major hepatic injuries who present in physiologic extremis have high survival rates despite prolonged hospital stays. Nonoperative management of major liver injuries, as diagnosed using computed tomography, is typically successful. Adjuncts (e.g., angioembolization, laparoscopic washouts, biliary stents) are essential in managing high-grade injuries. Injury to the extrahepatic biliary tree is rare. Cholecystectomy is indicated for all gallbladder trauma. Full-thickness common bile duct injuries require a hepaticojejunostomy, although damage control remains closed suction drainage. Injuries to the pancreatic head often involve concurrent trauma to regional vasculature. Damage control necessitates drainage after stopping hemorrhage. Injury to the left pancreas commonly requires a distal pancreatectomy. Outcomes for high-grade pancreatic and liver injuries are improved by involving an HPB team. Complications are multidisciplinary and should be managed without delay.

Delving into Early-onset Pancreatic Ductal Adenocarcinoma: How Does Age Fit In?

PURPOSE

With the rising incidence of early-onset pancreatic cancer (EOPC), molecular characteristics that distinguish early-onset pancreatic ductal adenocarcinoma (PDAC) tumors from those arising at a later age are not well understood.

EXPERIMENTAL DESIGN

We performed bioinformatic analysis of genomic and transcriptomic data generated from 269 advanced (metastatic or locally advanced) and 277 resectable PDAC tumor samples. Patient samples were stratified into EOPC (age of onset ≤55 years; n = 117), intermediate (age of onset 55-70 years; n = 264), and average (age of onset ≥70 years; n = 165) groups. Frequency of somatic mutations affecting genes commonly implicated in PDAC, as well as gene expression patterns, were compared between EOPC and all other groups.

RESULTS

EOPC tumors showed significantly lower frequency of somatic single-nucleotide variant (SNV)/insertions/deletions (indel) in CDKN2A (P = 0.0017), and were more likely to achieve biallelic mutation of CDKN2A through homozygous copy loss as opposed to heterozygous copy loss coupled with a loss-of-function SNV/indel mutation, the latter of which was more common for tumors with later ages of onset (P = 1.5e-4). Transcription factor forkhead box protein C2 (FOXC2) was significantly upregulated in EOPC tumors (P = 0.032). Genes significantly correlated with FOXC2 in PDAC samples were enriched for gene sets related to epithelial-to-mesenchymal transition (EMT) and included VIM (P = 1.8e-8), CDH11 (P = 6.5e-5), and CDH2 (P = 2.4e-2).

CONCLUSIONS

Our comprehensive analysis of sequencing data generated from a large cohort of PDAC patient samples highlights a distinctive pattern of biallelic CDKN2A mutation in EOPC tumors. Increased expression of FOXC2 in EOPC, with the correlation between FOXC2 and EMT pathways, represents novel molecular characteristics of EOPC.See related commentary by Lou, p. 8.

Oncogenic HBV variants and integration are present in hepatic and lymphoid cells derived from chronic HBV patients

The hepatitis B virus (HBV) is a major cause of hepatocellular carcinoma (HCC), partly driven by viral integration and specific oncogenic HBV variants. However, the biological significance of HBV genomes within lymphoid cells (i.e., peripheral blood mononuclear cells, PBMCs) is unclear. Here, we collected available plasma, PBMC, liver, and tumor from 52 chronic HBV (CHB) carriers: 32 with HCC, 19 without HCC, and one with dendritic cell sarcoma, DCS. Using highly sensitive sequencing techniques, next generation sequencing, and AluPCR, we demonstrate that viral genomes (i.e., HBV DNA, RNA, and cccDNA), oncogenic variants, and HBV-host integration are often found in all sample types collected from 52 patients (including lymphoid cells and a DCS tumor). Viral integration was recurrently identified (n = 90 such hits) in genes associated with oncogenic consequences in lymphoid and liver cells. Further, HBV genomes increased in PBMCs derived from 7 additional (treated or untreated) CHB carriers after extracellular mitogen stimulation. Our study shows novel HBV molecular data and replication not only liver, but also within 63.8% of lymphoid cells analysed (including a representative lymphoid cell malignancy), that was enhanced in ex vivo stimulated PBMC.

Lipid is heterogeneously distributed in muscle and associates with low radiodensity in cancer patients

BACKGROUND

Low muscle radiodensity is associated with mortality in a variety of cancer types. Biochemical and morphological correlates are unknown. We aimed to evaluate triglyceride (TG) content and location as a function of computed tomography (CT)-derived measures of skeletal muscle radiodensity in cancer patients.

METHODS

Rectus abdominis (RA) biopsies were collected during cancer surgery from 75 patients diagnosed with cancer. Thin-layer chromatography and gas chromatography were used for quantification of TG content of the muscle. Axial CT images of lumbar vertebra were used to measure muscle radiodensity. Oil Red O staining was used to determine the location of neutral lipids in frozen muscle sections.

RESULTS

There was wide variation in RA radiodensity in repeated measures (CV% ranged from 3 to 55% based on 10 serial images) as well as within one slice (CV% ranged from 6 to 61% based on 10 subregions). RA radiodensity and total lumbar muscle radiodensity were inversely associated with TG content of RA (r = -0.396, P < 0.001, and r = -0.355, P = 0.002, respectively). Of the total percentage area of muscle staining positive for neutral lipid, 54 ± 17% was present as extramyocellular lipids (range 23.5-77.8%) and 46 ± 17% (range 22.2-76.5%) present as intramyocellular lipid droplets.

CONCLUSIONS

Repeated measures revealed wide variation in radiodensity of RA muscle, both vertically and horizontally. Low muscle radiodensity reflects high level of TG in patients with cancer. Non-uniform distribution of intramyocellular and extramyocellular lipids was evident using light microscopy. These results warrant investigation of mechanisms resulting in lipid deposition in muscles of cancer patients.

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale1-3. Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter4; identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation5,6; analyses timings and patterns of tumour evolution7; describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity8,9; and evaluates a range of more-specialized features of cancer genomes8,10-18.

A combination of the immunohistochemical markers CK7 and SATB2 is highly sensitive and specific for distinguishing primary ovarian mucinous tumors from colorectal and appendiceal metastases

Primary ovarian mucinous tumors can be difficult to distinguish from metastatic gastrointestinal neoplasms by histology alone. The expected immunoprofile of a suspected metastatic lower gastrointestinal tumor is CK7-/CK20+/CDX2+/PAX8-. This study assesses the addition of a novel marker SATB2, to improve the diagnostic algorithm. A test cohort included 155 ovarian mucinous tumors (105 carcinomas and 50 borderline tumors) and 230 primary lower gastrointestinal neoplasms (123 colorectal adenocarcinomas and 107 appendiceal neoplasms). All cases were assessed for SATB2, PAX8 CK7, CK20, and CDX2 expression on tissue microarrays. Expression was scored in a 3-tier system as absent, focal (1-50% of tumor cells) and diffuse ( >50% of tumor cells) and then categorized into either absent/present or nondiffuse/diffuse. SATB2 and PAX8 expression was further evaluated in ovarian tumors from an international cohort of 2876 patients (expansion cohort, including 159 mucinous carcinomas and 46 borderline mucinous tumors). The highest accuracy of an individual marker in distinguishing lower gastrointestinal from ovarian mucinous tumors was CK7 (91.7%, nondiffuse/diffuse cut-off) followed by SATB2 (88.8%, present/absent cut-off). The most effective combination was CK7 and SATB2 with accuracy of 95.3% using the 3-tier interpretation, absent/focal/diffuse. This combination outperformed the standard clinical set of CK7, CK20 and CDX2 (87.5%). Re-evaluation of outlier cases confirmed ovarian origin for all but one case. The accuracy of SATB2 was confirmed in the expansion cohort (91.5%). SATB2 expression was also detected in 15% of ovarian endometrioid carcinoma but less than 5% of other ovarian histotypes. A simple two marker combination of CK7 and SATB2 can distinguish lower gastrointestinal from ovarian primary mucinous tumors with greater than 95% accuracy. PAX8 and CDX2 have value as second-line markers. The utility of CK20 in this setting is low and this warrants replacement of this marker with SATB2 in clinical practice.

Clinical and biological characterization of skeletal muscle tissue biopsies of surgical cancer patients

BACKGROUND

Researchers increasingly use intraoperative muscle biopsy to investigate mechanisms of skeletal muscle atrophy in patients with cancer. Muscles have been assessed for morphological, cellular, and biochemical features. The aim of this study was to conduct a state-of-the-science review of this literature and, secondly, to evaluate clinical and biological variation in biopsies of rectus abdominis (RA) muscle from a cohort of patients with malignancies.

METHODS

Literature was searched for reports on muscle biopsies from patients with a cancer diagnosis. Quality of reports and risk of bias were assessed. Data abstracted included patient characteristics and diagnoses, sample size, tissue collection and biobanking procedures, and results. A cohort of cancer patients (n = 190, 88% gastrointestinal malignancies), who underwent open abdominal surgery as part of their clinical care, consented to RA biopsy from the site of incision. Computed tomography (CT) scans were used to quantify total abdominal muscle and RA cross-sectional areas and radiodensity. Biopsies were assessed for muscle fibre area (μm² ), fibre types, myosin heavy chain isoforms, and expression of genes selected for their involvement in catabolic pathways of muscle.

RESULTS

Muscle biopsy occurred in 59 studies (total N = 1585 participants). RA was biopsied intraoperatively in 40 studies (67%), followed by quadriceps (26%; percutaneous biopsy) and other muscles (7%). Cancer site and stage, % of male participants, and age were highly variable between studies. Details regarding patient medical history and biopsy procedures were frequently absent. Lack of description of the population(s) sampled and low sample size contributed to low quality and risk of bias. Weight-losing cases were compared with weight stable cancer or healthy controls without considering a measure of muscle mass in 21 out of 44 studies. In the cohort of patients providing biopsy for this study, 78% of patients had preoperative CT scans and a high proportion (64%) met published criteria for sarcopenia. Fibre type distribution in RA was type I (46% ± 13), hybrid type I/IIA (1% ± 1), type IIA (36% ± 10), hybrid type IIA/D (15% ± 14), and type IID (2% ± 5). Sexual dimorphism was prominent in RA CT cross-sectional area, mean fibre cross-sectional area, and in expression of genes associated with muscle growth, apoptosis, and inflammation (P < 0.05). Medical history revealed multiple co-morbid conditions and medications.

CONCLUSIONS

Continued collaboration between researchers and cancer surgeons enables a more complete understanding of mechanisms of cancer-associated muscle atrophy. Standardization of biobanking practices, tissue manipulation, patient characterization, and classification will enhance the consistency, reliability, and comparability of future studies.

Immunohistochemical phenotyping of T cells, granulocytes, and phagocytes in the muscle of cancer patients: association with radiologically defined muscle mass and gene expression

BACKGROUND

Inflammation is a recognized contributor to muscle wasting. Research in injury and myopathy suggests that interactions between the skeletal muscle and immune cells confer a pro-inflammatory environment that influences muscle loss through several mechanisms; however, this has not been explored in the cancer setting. This study investigated the local immune environment of the muscle by identifying the phenotype of immune cell populations in the muscle and their relationship to muscle mass in cancer patients.

METHODS

Intraoperative muscle biopsies were collected from cancer patients (n = 30, 91% gastrointestinal malignancies). Muscle mass was assessed histologically (muscle fiber cross-sectional area, CSA; μm²) and radiologically (lumbar skeletal muscle index, SMI; cm²/m² by computed tomography, CT). T cells (CD4 and CD8) and granulocytes/phagocytes (CD11b, CD14, and CD15) were assessed by immunohistochemistry. Microarray analysis was conducted in the muscle of a second cancer patient cohort.

RESULTS

T cells (CD3+), granulocytes/phagocytes (CD11b+), and CD3-CD4+ cells were identified. Muscle fiber CSA (μm²) was positively correlated (Spearman's r = > 0.45; p = < 0.05) with the total number of T cells, CD4, and CD8 T cells and granulocytes/phagocytes. In addition, patients with the smallest SMI exhibited fewer CD8 T cells within their muscle. Consistent with this, further exploration with gene correlation analyses suggests that the presence of CD8 T cells is negatively associated (Pearson's r = ≥ 0.5; p = <0.0001) with key genes within muscle catabolic pathways for signaling (ACVR2B), ubiquitin proteasome (FOXO4, TRIM63, FBXO32, MUL1, UBC, UBB, UBE2L3), and apoptosis/autophagy (CASP8, BECN1, ATG13, SIVA1).

CONCLUSION

The skeletal muscle immune environment of cancer patients is comprised of immune cell populations from the adaptive and innate immunity. Correlations of T cells, granulocyte/phagocytes, and CD3-CD4+ cells with muscle mass measurements indicate a positive relationship between immune cell numbers and muscle mass status in cancer patients. Further exploration with gene correlation analyses suggests that the presence of CD8 T cells is negatively correlated with components of muscle catabolism.

Molecular Characterization and Clinical Relevance of Metabolic Expression Subtypes in Human Cancers

Metabolic reprogramming provides critical information for clinical oncology. Using molecular data of 9,125 patient samples from The Cancer Genome Atlas, we identified tumor subtypes in 33 cancer types based on mRNA expression patterns of seven major metabolic processes and assessed their clinical relevance. Our metabolic expression subtypes correlated extensively with clinical outcome: subtypes with upregulated carbohydrate, nucleotide, and vitamin/cofactor metabolism most consistently correlated with worse prognosis, whereas subtypes with upregulated lipid metabolism showed the opposite. Metabolic subtypes correlated with diverse somatic drivers but exhibited effects convergent on cancer hallmark pathways and were modulated by highly recurrent master regulators across cancer types. As a proof-of-concept example, we demonstrated that knockdown of SNAI1 or RUNX1—master regulators of carbohydrate metabolic subtypes—modulates metabolic activity and drug sensitivity. Our study provides a system-level view of metabolic heterogeneity within and across cancer types and identifies pathway cross-talk, suggesting related prognostic, therapeutic, and predictive utility.

A strategy for early detection of response to chemotherapy drugs based on treatment-related changes in the metabolome

We describe a biomarker-based approach to delivering chemotherapy that entails monitoring treatment changes in the circulating metabolome that reflect efficacy. In-vitro, multiple tumor cell lines were exposed to numerous chemotherapeutics. Supernatants were collected at baseline and 72 hours post treatment. MTT assays were used to quantify growth inhibition. Clinical samples were derived from a phase II clinical trial of second-line axitinib in patients with advanced hepatocellular carcinoma. Sera were collected at baseline and 2–4 weeks after treatment initiation. Response to therapy was estimated by CT scan at 8 weeks. Samples were analyzed by gas chromatography-mass spectrometry to identify metabolomic changes associated with response. In vitro, we found drug-specific and generalizable patterns of change in the extracellular metabolome accompany growth inhibition. A cell death signature was also identified. This approach was also applied to clinical samples. While the in vitro signatures were detectable in vivo, a more robust signal was identified clinically that appeared within 4 weeks of administering drug that distinguished individuals with a treatment response. These changes were extinguished as tumor growth resumed. Serial monitoring of the metabolome during chemotherapy is a means to follow treatment efficacy and emergence of resistance, informing the oncologist whether to modify treatment.

Acute kidney injury following resection of hepatocellular carcinoma: prognostic value of the acute kidney injury network criteria

Background

Acute kidney injury (AKI) is associated with increased morbidity and mortality after liver resection. Patients with hepatocellular carcinoma (HCC) have a higher risk of AKI owing to the underlying association between hepatic and renal dysfunction. Use of the Acute Kidney Injury Network (AKIN) diagnostic criteria is recommended for patients with cirrhosis, but remains poorly studied following liver resection. We compared the prognostic value of the AKIN creatinine and urine output criteria in terms of postoperative outcomes following liver resection for HCC.

Methods

All patients who underwent a liver resection for HCC from January 2010 to June 2016 were included. We used AKIN urine output and creatinine criteria to assess for AKI within 48 hours of surgery.

Results

Eighty liver resections were performed during the study period. Cirrhosis was confirmed in 80%. Median hospital stay was 9 (interquartile range 7–12) days, and 30-day mortality was 2.5%. The incidence of AKI was higher based on the urine output than on the creatinine criterion (53.8% v. 20%), and was associated with prolonged hospitalization and 30-day postoperative mortality when defined by serum creatinine (hospital stay: 11.2 v. 20.1 d, p = 0.01; mortality: 12.5% v. 0%, p < 0.01), but not urine output (hospital stay: 15.6 v. 10 d, p = 0.05; mortality: 2.3% v. 2.7%, p > 0.99).

Conclusion

The urine output criterion resulted in an overestimation of AKI and compromised the prognostic value of AKIN criteria. Revision may be required to account for the exacerbated physiologic postoperative reduction in urine output in patients with HCC.

Host phenotype is associated with reduced survival independent of tumour biology in patients with colorectal liver metastases

BACKGROUND

Most prognostic scoring systems for colorectal liver metastases (CRLMs) account for factors related to tumour biology. Little is known about the effects of the host phenotype to the tumour. Our objective was to delineate the relationship of systemic inflammation and body composition features [i.e. low skeletal muscle mass (sarcopenia) and low visceral adipose tissue (VAT)], two well-described host phenotypes in cancer.

METHODS

Clinical data and pre-operative blood samples were collected from 99 patients who underwent resection of CRLM. Pre-operative computed tomography scans were available for 97 patients; body composition was analysed at the L3 level, stratified for sex and age. Clinicopathological variables, serum C-reactive protein (CRP), and various body composition variables were evaluated. Overall survival was evaluated as a function of these same variables in multivariate Cox regression analysis.

RESULTS

Skeletal muscle was significantly correlated with VAT (r = 0.46, P < 0.001). Of patients with sarcopenia, 35 (65%) also had low VAT. C-reactive protein was elevated (≥5 mg/mL) in 42 patients (43.3%). Elevated CRP was more common in patients with sarcopenia (73.8% vs. 51.1%, P = 0.029). The most significant prognostic factors were the coincidence of elevated CRP and adverse body composition features (sarcopenia and/or low VAT; hazard ratio 4.3, 95% confidence interval 1.5-13.0, P = 0.008), as well as Fong clinical prognostic score (hazard ratio 2.9, 95% confidence interval 1.5-5.5, P = 0.002).

CONCLUSIONS

Body composition in patients with CRLM is not directly linked to the presence of systemic inflammation. However, when systemic inflammation coincides with sarcopenia and/or low VAT, prognosis is adversely affected, independent of the Fong clinical prognostic score.

Effect of proximity to specialty care on outcomes for biliary cancers: a population-based retrospective cohort study

BACKGROUND

The management of biliary cancers is complex and requires a multidisciplinary approach. Because it is unknown how access to specialty care affects resource use and survival in patients with biliary cancer, we conducted a population-based study to understand the needs of these patients and the relation of geography to care delivery and clinical outcomes for biliary cancer in Alberta.

METHODS

All patients with biliary cancer diagnosed in Alberta from Sept. 1, 2001, to Dec. 31, 2015 were included in this population-based retrospective cohort study. Data were extracted from administrative databases and the 2011 Canadian census. Driving time and types of medical services were tracked throughout the patients' clinical course. We categorized proximity to specialty care according to driving time to the nearest specialist. The primary outcome was overall survival. We conducted Cox proportional hazard regression to evaluate the effects of driving time on overall survival and multivariate logistic regression to evaluate the effect of driving time on treatment types and stage at diagnosis.

RESULTS

We identified 1610 patients with biliary cancer; they accounted for 117 381 medical encounters. Patients living 120 minutes or more from the nearest hepatobiliary surgeon and from the nearest cancer centre had significantly decreased survival (hazard ratio [and 95% confidence interval (CI)] 1.27 [1.17-1.37]) and 1.27 [1.14-1.41], respectively). Location of residence was not associated with advanced stage or probability of undergoing surgery or a biliary drainage procedure. Patients who lived 120 minutes or more from a cancer centre were less likely than those who lived less than 120 minutes away to receive chemotherapy (odds ratio 0.51, 95% CI 0.29-0.88). Subgroup analysis showed that the effect of travel time was especially pronounced among those who received only best supportive care and those who had biliary drains.

INTERPRETATION

Geography and accessibility to specialty care affected survival in patients with biliary cancer. Further study is required to understand how patients with biliary drains and those receiving best supportive care are affected by proximity to specialty care. This will aid in the identification of strategies to provide improved care for this subgroup who are particularly affected by geography.

A Framework for Development of Useful Metabolomic Biomarkers and Their Effective Knowledge Translation

Despite the significant advantages of metabolomic biomarkers, no diagnostic tests based on metabolomics have been introduced to clinical use. There are many reasons for this, centered around substantial obstacles in developing clinically useful metabolomic biomarkers. Most significant is the need for interdisciplinary teams with expertise in metabolomics, analysis of complex clinical and metabolomic data, and clinical care. Importantly, the clinical need must precede biomarker discovery, and the experimental design for discovery and validation must reflect the purpose of the biomarker. Standard operating procedures for procuring and handling samples must be developed from the beginning, to ensure experimental integrity. Assay design is another challenge, as there is not much precedent informing this. Another obstacle is that it is not yet clear how to protect any intellectual property related to metabolomic biomarkers. Viewing a metabolomic biomarker as a natural phenomenon would inhibit patent protection and potentially stifle commercial interest. However, demonstrating that a metabolomic biomarker is actually a derivative of a natural phenomenon that requires innovation would enhance investment in this field. Finally, effective knowledge translation strategies must be implemented, which will require engagement with end users (clinicians and lab physicians), patient advocate groups, policy makers, and payer organizations. Addressing each of these issues comprises the framework for introducing a metabolomic biomarker to practice.

Knockout of USP19 Deubiquitinating Enzyme Prevents Muscle Wasting by Modulating Insulin and Glucocorticoid Signaling

Muscle atrophy arises because of many chronic illnesses, as well as from prolonged glucocorticoid treatment and nutrient deprivation. We previously demonstrated that the USP19 deubiquitinating enzyme plays an important role in chronic glucocorticoid- and denervation-induced muscle wasting. However, the mechanisms by which USP19 exerts its effects remain unknown. To explore this further, we fasted mice for 48 hours to try to identify early differences in the response of wild-type and USP19 knockout (KO) mice that could yield insights into the mechanisms of USP19 action. USP19 KO mice manifested less myofiber atrophy in response to fasting due to increased rates of protein synthesis. Insulin signaling was enhanced in the KO mice, as revealed by lower circulating insulin levels, increased insulin-stimulated glucose disposal and phosphorylation of Akt and S6K in muscle, and improved overall glucose tolerance. Glucocorticoid signaling, which is essential in many conditions of atrophy, was decreased in KO muscle, as revealed by decreased expression of glucocorticoid receptor (GR) target genes upon both fasting and glucocorticoid treatment. This decreased GR signaling was associated with lower GR protein levels in the USP19 KO muscle. Restoring the GR levels in USP19-deficient muscle was sufficient to abolish the protection from myofiber atrophy. Expression of GR target genes also correlated with that of USP19 in human muscle samples. Thus, USP19 modulates GR levels and in so doing may modulate both insulin and glucocorticoid signaling, two critical pathways that control protein turnover in muscle and overall glucose homeostasis.

Sarcopenia and myosteatosis are accompanied by distinct biological profiles in patients with pancreatic and periampullary adenocarcinomas

Introduction

Pancreatic and periampullary adenocarcinomas are associated with abnormal body composition visible on CT scans, including low muscle mass (sarcopenia) and low muscle radiodensity due to fat infiltration in muscle (myosteatosis). The biological and clinical correlates to these features are poorly understood.

Methods

Clinical characteristics and outcomes were studied in 123 patients who underwent pancreaticoduodenectomy for pancreatic or non-pancreatic periampullary adenocarcinoma and who had available preoperative CT scans. In a subgroup of patients with pancreatic cancer (n = 29), rectus abdominus muscle mRNA expression was determined by cDNA microarray and in another subgroup (n = 29) ¹H-NMR spectroscopy and gas chromatography-mass spectrometry were used to characterize the serum metabolome.

Results

Muscle mass and radiodensity were not significantly correlated. Distinct groups were identified: sarcopenia (40.7%), myosteatosis (25.2%), both (11.4%). Fat distribution differed in these groups; sarcopenia associated with lower subcutaneous adipose tissue (P<0.0001) and myosteatosis associated with greater visceral adipose tissue (P<0.0001). Sarcopenia, myosteatosis and their combined presence associated with shorter survival, Log Rank P = 0.005, P = 0.06, and P = 0.002, respectively. In muscle, transcriptomic analysis suggested increased inflammation and decreased growth in sarcopenia and disrupted oxidative phosphorylation and lipid accumulation in myosteatosis. In the circulating metabolome, metabolites consistent with muscle catabolism associated with sarcopenia. Metabolites consistent with disordered carbohydrate metabolism were identified in both sarcopenia and myosteatosis.

Discussion

Muscle phenotypes differ clinically and biologically. Because these muscle phenotypes are linked to poor survival, it will be imperative to delineate their pathophysiologic mechanisms, including whether they are driven by variable tumor biology or host response.

Differentially expressed alternatively spliced genes in skeletal muscle from cancer patients with cachexia

BACKGROUND

Alternative splicing (AS) is a post-transcriptional gene regulatory mechanism that contributes to proteome diversity. Aberrant splicing mechanisms contribute to various cancers and muscle-related conditions such as Duchenne muscular dystrophy. However, dysregulation of AS in cancer cachexia (CC) remains unexplored. Our objectives were (i) to profile alternatively spliced genes (ASGs) on a genome-wide scale and (ii) to identify differentially expressed alternatively spliced genes (DASGs) associated with CC.

METHODS

Rectus abdominis muscle biopsies obtained from cancer patients were stratified into cachectic cases (n = 21, classified based on International consensus diagnostic framework for CC) and non-cachectic controls (n = 19, weight stable cancer patients). Human transcriptome array 2.0 was used for profiling ASGs using the total RNA isolated from muscle biopsies. Representative DASG signatures were validated using semi-quantitative RT-PCR.

RESULTS

We identified 8960 ASGs, of which 922 DASGs (772 up-regulated and 150 down-regulated) were identified at ≥1.4 fold-change and P < 0.05. Representative DASGs validated by semi-quantitative RT-PCR confirmed the primary findings from the human transcriptome arrays. Identified DASGs were associated with myogenesis, adipogenesis, protein ubiquitination, and inflammation. Up to 10% of the DASGs exhibited cassette exon (exon included or skipped) as a predominant form of AS event. We also observed other forms of AS events such as intron retention, alternate promoters.

CONCLUSIONS

Overall, we have, for the first time, conducted global profiling of muscle tissue to identify DASGs associated with CC. The mechanistic roles of the identified DASGs in CC pathophysiology using model systems is warranted, as well as replication of findings in independent cohorts.

A quantitative multimodal metabolomic assay for colorectal cancer

Background

Early diagnosis of colorectal cancer (CRC) simplifies treatment and improves treatment outcomes. We previously described a diagnostic metabolomic biomarker derived from semi-quantitative gas chromatography-mass spectrometry. Our objective was to determine whether a quantitative assay of additional metabolomic features, including parts of the lipidome could enhance diagnostic power; and whether there was an advantage to deriving a combined diagnostic signature with a broader metabolomic representation.

Methods

The well-characterized Biocrates P150 kit was used to quantify 163 metabolites in patients with CRC (N = 62), adenoma (N = 31), and age- and gender-matched disease-free controls (N = 81). Metabolites included in the analysis included phosphatidylcholines, sphingomyelins, acylcarnitines, and amino acids. Using a training set of 32 CRC and 21 disease-free controls, a multivariate metabolomic orthogonal partial least squares (OPLS) classifier was developed. An independent set of 28 CRC and 20 matched healthy controls was used for validation. Features characterizing 31 colorectal adenomas from their healthy matched controls were also explored, and a multivariate OPLS classifier for colorectal adenoma could be proposed.

Results

The metabolomic profile that distinguished CRC from controls consisted of 48 metabolites (R²Y = 0.83, Q²Y = 0.75, CV-ANOVA p-value < 0.00001). In this quantitative assay, the coefficient of variance for each metabolite was <10%, and this dramatically enhanced the separation of these groups. Independent validation resulted in AUROC of 0.98 (95% CI, 0.93–1.00) and sensitivity and specificity of 93% and 95%. Similarly, we were able to distinguish adenoma from controls (R²Y = 0.30, Q²Y = 0.20, CV-ANOVA p-value = 0.01; internal AUROC = 0.82 (95% CI, 0.72–0.93)). When combined with the previously generated GC-MS signatures for CRC and adenoma, the candidate biomarker performance improved slightly.

Conclusion

The diagnostic power for metabolomic tests for colorectal neoplasia can be improved by utilizing a multimodal approach and combining metabolites from diverse chemical classes. In addition, quantification of metabolites enhances separation of disease-specific metabolomic profiles. Our future efforts will be focused on developing a quantitative assay for the metabolites comprising the optimal diagnostic biomarker.

Electronic supplementary material

The online version of this article (10.1186/s12885-017-3923-z) contains supplementary material, which is available to authorized users.

Monitoring for Response to Antineoplastic Drugs: The Potential of a Metabolomic Approach

For most cancers, chemotherapeutic options are rapidly expanding, providing the oncologist with substantial choices. Therefore, there is a growing need to select the best systemic therapy, for any individual, that effectively halts tumor progression with minimal toxicity. Having the capability to predict benefit and to anticipate toxicity would be ideal, but remains elusive at this time. An alternative approach is an adaptive approach that involves close observation for treatment response and emergence of resistance. Currently, response to systemic therapy is estimated using radiographic tests. Unfortunately, radiographic estimates of response are imperfect and radiographic signs of response can be delayed. This is particularly problematic for targeted agents, as tumor shrinkage is often not apparent with these drugs. As a result, patients are exposed to prolonged courses of toxic drugs that may ultimately be found to be ineffective. A biomarker-based adaptive strategy that involves the serial analysis of the metabolome is attractive. The metabolome changes rapidly with changes in physiology. Changes in the circulating metabolome associated with various antineoplastic agents have been described, but further work will be required to understand what changes signify clinical benefit. We present an investigative approach for the discovery and validation of metabolomic response biomarkers, which consists of serial analysis of the metabolome and linkage of changes in the metabolome to measurable therapeutic benefit. Potential pitfalls in the development of metabolomic biomarkers of response and loss of response are reviewed.

Integrative Genomic Analysis of Cholangiocarcinoma Identifies Distinct IDH-Mutant Molecular Profiles

Cholangiocarcinoma (CCA) is an aggressive malignancy of the bile ducts, with poor prognosis and limited treatment options. Here, we describe the integrated analysis of somatic mutations, RNA expression, copy number, and DNA methylation by The Cancer Genome Atlas, of a set of predominantly intrahepatic CCA cases, and propose a molecular classification scheme. We identified an IDH -mutant enriched subtype with distinct molecular features including low expression of chromatin modifiers, elevated expression of mitochondrial genes, and increased mitochondrial DNA copy number. Leveraging the multi-platform data, we observed that ARID1A exhibited DNA hypermethylation and decreased expression in the IDH -mutant subtype. More broadly, we found that IDH mutations are associated with an expanded histological spectrum of liver tumors with molecular features that stratify with CCA. Our studies reveal insights into the molecular pathogenesis and heterogeneity of cholangiocarcinoma and provide classification information of potential therapeutic significance.

Small RNAome profiling from human skeletal muscle: novel miRNAs and their targets associated with cancer cachexia

BACKGROUND

MicroRNAs (miRs) are small non-coding RNAs that regulate gene (mRNA) expression. Although the pathological role of miRs have been studied in muscle wasting conditions such as myotonic and muscular dystrophy, their roles in cancer cachexia (CC) are still emerging.

OBJECTIVES

The objectives are (i) to profile human skeletal muscle expressed miRs; (ii) to identify differentially expressed (DE) miRs between cachectic and non-cachectic cancer patients; (iii) to identify mRNA targets for the DE miRs to gain mechanistic insights; and (iv) to investigate if miRs show potential prognostic and predictive value.

METHODS

Study subjects were classified based on the international consensus diagnostic criteria for CC. Forty-two cancer patients were included, of which 22 were cachectic cases and 20 were non-cachectic cancer controls. Total RNA isolated from muscle biopsies were subjected to next-generation sequencing.

RESULTS

A total of 777 miRs were profiled, and 82 miRs with read counts of ≥5 in 80% of samples were retained for analysis. We identified eight DE miRs (up-regulated, fold change of ≥1.4 at P < 0.05). A total of 191 potential mRNA targets were identified for the DE miRs using previously described human skeletal muscle mRNA expression data (n = 90), and a majority of them were also confirmed in an independent mRNA transcriptome dataset. Ingenuity pathway analysis identified pathways related to myogenesis and inflammation. qRT-PCR analysis of representative miRs showed similar direction of effect (P < 0.05), as observed in next-generation sequencing. The identified miRs also showed prognostic and predictive value.

CONCLUSIONS

In all, we identified eight novel miRs associated with CC.

Performance of variable selection methods using stability-based selection

BACKGROUND

Variable selection is frequently carried out during the analysis of many types of high-dimensional data, including those in metabolomics. This study compared the predictive performance of four variable selection methods using stability-based selection, a new secondary selection method that is implemented in the R package BioMark. Two of these methods were evaluated using the more well-known false discovery rate (FDR) as well.

RESULTS

Simulation studies varied factors relevant to biological data studies, with results based on the median values of 200 partial area under the receiver operating characteristic curve. There was no single top performing method across all factor settings, but the student t test based on stability selection or with FDR adjustment and the variable importance in projection (VIP) scores from partial least squares regression models obtained using a stability-based approach tended to perform well in most settings. Similar results were found with a real spiked-in metabolomics dataset. Group sample size, group effect size, number of significant variables and correlation structure were the most important factors whereas the percentage of significant variables was the least important.

CONCLUSIONS

Researchers can improve prediction scores for their study data by choosing VIP scores based on stability variable selection over the other approaches when the number of variables is small to modest and by increasing the number of samples even moderately. When the number of variables is high and there is block correlation amongst the significant variables (i.e., true biomarkers), the FDR-adjusted student t test performed best. The R package BioMark is an easy-to-use open-source program for variable selection that had excellent performance characteristics for the purposes of this study.

Integrative Genomic Analysis of Cholangiocarcinoma Identifies Distinct IDH-Mutant Molecular Profiles

Cholangiocarcinoma (CCA) is an aggressive malignancy of the bile ducts, with poor prognosis and limited treatment options. Here, we describe the integrated analysis of somatic mutations, RNA expression, copy number, and DNA methylation by The Cancer Genome Atlas of a set of predominantly intrahepatic CCA cases and propose a molecular classification scheme. We identified an IDH mutant-enriched subtype with distinct molecular features including low expression of chromatin modifiers, elevated expression of mitochondrial genes, and increased mitochondrial DNA copy number. Leveraging the multi-platform data, we observed that ARID1A exhibited DNA hypermethylation and decreased expression in the IDH mutant subtype. More broadly, we found that IDH mutations are associated with an expanded histological spectrum of liver tumors with molecular features that stratify with CCA. Our studies reveal insights into the molecular pathogenesis and heterogeneity of cholangiocarcinoma and provide classification information of potential therapeutic significance.

A population-based study of biliary tract cancers (BTCs) in Alberta, Canada: How do our patients fare?

246

Background: BTCs are poorly studied due to their rarity and heterogeneity. We explored demographics and outcomes of BTC pts over a 15 year period. Methods: All patients (pts) with biopsy-proven BTC (intrahepatic [IC] and extrahepatic [EC] cholangiocarcinomas, gallbladder cancers [GB], and ampulla of vater cancers [AV]) in Alberta were reviewed from January 1, 2000, to December 31, 2015. Demographic, pathologic, and survival data were extracted from electronic charts. Descriptive statistics were utilized. Overall survival (OS) was defined as the time from pathologic diagnosis to death date. Results: A total of 1,718 pts with BTCs were identified. Median age was 68 with 52% being male. The Table demonstrates OS breakdown based on tumour location and stage. Regardless of location of primary tumour, grade impacted survival (median OS in with well differentiated tumours vs undifferentiated tumors 26.6 vs 3.9 months). Pts who received standard of care palliative cisplatin/gemcitabine (Cis/Gem) chemotherapy (n = 123) had a median OS of 15.4 months. Conclusions: Patients with AV and IC have the best and worst prognosis, respectively. Shorter survival is observed with higher stage, grade, and unresectable disease. Pts who received palliative Cis/Gem had better OS than reported in the pivotal phase III trial. Further analysis of prognostic factors will be presented. [Table: see text]

New genetic signatures associated with cancer cachexia as defined by low skeletal muscle index and weight loss

BACKGROUND

Cachexia affects the majority with advanced cancer. Based on current demographic and clinical factors, it is not possible to predict who will develop cachexia or not. Such variation may, in part, be due to genotype. It has recently been proposed to extend the diagnostic criteria for cachexia to include a direct measure of low skeletal muscle index (LSMI) in addition to weight loss (WL). We aimed to explore our panel of candidate single nucleotide polymorphism (SNPs) for association with WL +/- computerized tomography-defined LSMI. We also explored whether the transcription in muscle of identified genes was altered according to such cachexia phenotype METHODS: A retrospective cohort study design was used. Analysis explored associations of candidate SNPs with WL (n = 1276) and WL + LSMI (n = 943). Human muscle transcriptome (n = 134) was analysed using an Agilent platform.

RESULTS

Single nucleotide polymorphisms in the following genes showed association with WL alone: GCKR, LEPR, SELP, ACVR2B, TLR4, FOXO3, IGF1, CPN1, APOE, FOXO1, and GHRL. SNPs in LEPR, ACVR2B, TNF, and ACE were associated with concurrent WL + LSMI. There was concordance between muscle-specific expression for ACVR2B, FOXO1 and 3, LEPR, GCKR, and TLR4 genes and LSMI and/or WL (P < 0.05).

CONCLUSIONS

The rs1799964 in the TNF gene and rs4291 in the ACE gene are new associations when the definition of cachexia is based on a combination of WL and LSMI. These findings focus attention on pro-inflammatory cytokines and the renin-angiotensin system as biomarkers/mediators of muscle wasting in cachexia.