Pancreatic adenocarcinoma presents as a spectrum of a highly aggressive disease in patients. The basis of this disease heterogeneity has proved difficult to resolve due to poor tumor cellularity and ...extensive genomic instability. To address this, a dataset of whole genomes and transcriptomes was generated from purified epithelium of primary and metastatic tumors. Transcriptome analysis demonstrated that molecular subtypes are a product of a gene expression continuum driven by a mixture of intratumoral subpopulations, which was confirmed by single-cell analysis. Integrated whole-genome analysis uncovered that molecular subtypes are linked to specific copy number aberrations in genes such as mutant KRAS and GATA6. By mapping tumor genetic histories, tetraploidization emerged as a key mutational process behind these events. Taken together, these data support the premise that the constellation of genomic aberrations in the tumor gives rise to the molecular subtype, and that disease heterogeneity is due to ongoing genomic instability during progression.
To determine the impact of basal-like and classical subtypes in advanced pancreatic ductal adenocarcinoma (PDAC) and to explore GATA6 expression as a surrogate biomarker.
Within the COMPASS trial, ...patients proceeding to chemotherapy for advanced PDAC undergo tumor biopsy for RNA-sequencing (RNA-seq). Overall response rate (ORR) and overall survival (OS) were stratified by subtypes and according to chemotherapy received. Correlation of
with the subtypes using gene expression profiling,
hybridization (ISH) was explored.
Between December 2015 and May 2019, 195 patients (95%) had enough tissue for RNA-seq; 39 (20%) were classified as basal-like and 156 (80%) as classical. RECIST response data were available for 157 patients; 29 basal-like and 128 classical where the ORR was 10% versus 33%, respectively (
= 0.02). In patients with basal-like tumors treated with modified FOLFIRINOX (
= 22), the progression rate was 60% compared with 15% in classical PDAC (
= 0.0002). Median OS in the intention-to-treat population (
= 195) was 9.3 months for classical versus 5.9 months for basal-like PDAC (HR, 0.47; 95% confidence interval, 0.32-0.69;
= 0.0001).
expression by RNA-seq highly correlated with the classifier (
< 0.001) and ISH predicted the subtypes with sensitivity of 89% and specificity of 83%. In a multivariate analysis, GATA6 expression was prognostic (
= 0.02). In exploratory analyses, basal-like tumors, could be identified by keratin 5, were more hypoxic and enriched for a T-cell-inflamed gene expression signature.
The basal-like subtype is chemoresistant and can be distinguished from classical PDAC by GATA6 expression.
.
We integrated clinical, genomic, and transcriptomic data from 224 primaries and 95 metastases from 289 patients to characterize progression of pancreatic ductal adenocarcinoma (PDAC). Driver gene ...alterations and mutational and expression-based signatures were preserved, with truncations, inversions, and translocations most conserved. Cell cycle progression (CCP) increased with sequential inactivation of tumor suppressors, yet remained higher in metastases, perhaps driven by cell cycle regulatory gene variants. Half of the cases were hypoxic by expression markers, overlapping with molecular subtypes. Paired tumor heterogeneity showed cancer cell migration by Halstedian progression. Multiple PDACs arising synchronously and metachronously in the same pancreas were actually intra-parenchymal metastases, not independent primary tumors. Established clinical co-variates dominated survival analyses, although CCP and hypoxia may inform clinical practice.
Display omitted
•Higher cell cycle progression in PDAC metastases; increases with driver gene loss•Half of PDACs are hypoxic and are associated with subtypes and treatment response•Paired tumors show molecular conservation and Halstedian progression•Multiple PDACs arising in the same pancreas are intra-parenchymal metastases
Connor et al. molecularly characterize primary and metastatic PDAC and show conserved alterations between primary and metastatic lesions. Clinical features outperform molecular alterations in survival analyses, but cell cycle progression and hypoxia signatures may inform clinical practice.
Homologous recombination deficiency (HRD) in pancreatic ductal adenocarcinoma (PDAC), remains poorly defined beyond germline (g) alterations in BRCA1, BRCA2, and PALB2.
We interrogated whole genome ...sequencing (WGS) data on 391 patients, including 49 carriers of pathogenic variants (PVs) in gBRCA and PALB2. HRD classifiers were applied to the dataset and included (1) the genomic instability score (GIS) used by Myriad’s MyChoice HRD assay; (2) substitution base signature 3 (SBS3); (3) HRDetect; and (4) structural variant (SV) burden. Clinical outcomes and responses to chemotherapy were correlated with HRD status.
Biallelic tumor inactivation of gBRCA or PALB2 was evident in 43 of 49 germline carriers identifying HRD-PDAC. HRDetect (score ≥0.7) predicted gBRCA1/PALB2 deficiency with highest sensitivity (98%) and specificity (100%). HRD genomic tumor classifiers suggested that 7% to 10% of PDACs that do not harbor gBRCA/PALB2 have features of HRD. Of the somatic HRDetecthi cases, 69% were attributed to alterations in BRCA1/2, PALB2, RAD51C/D, and XRCC2, and a tandem duplicator phenotype. TP53 loss was more common in BRCA1- compared with BRCA2-associated HRD-PDAC. HRD status was not prognostic in resected PDAC; however in advanced disease the GIS (P = .02), SBS3 (P = .03), and HRDetect score (P = .005) were predictive of platinum response and superior survival. PVs in gATM (n = 6) or gCHEK2 (n = 2) did not result in HRD-PDAC by any of the classifiers. In 4 patients, BRCA2 reversion mutations associated with platinum resistance.
Germline and parallel somatic profiling of PDAC outperforms germline testing alone in identifying HRD-PDAC. An additional 7% to 10% of patients without gBRCA/PALB2 mutations may benefit from DNA damage response agents.
Pancreatic cancer, a highly aggressive tumour type with uniformly poor prognosis, exemplifies the classically held view of stepwise cancer development. The current model of tumorigenesis, based on ...analyses of precursor lesions, termed pancreatic intraepithelial neoplasm (PanINs) lesions, makes two predictions: first, that pancreatic cancer develops through a particular sequence of genetic alterations (KRAS, followed by CDKN2A, then TP53 and SMAD4); and second, that the evolutionary trajectory of pancreatic cancer progression is gradual because each alteration is acquired independently. A shortcoming of this model is that clonally expanded precursor lesions do not always belong to the tumour lineage, indicating that the evolutionary trajectory of the tumour lineage and precursor lesions can be divergent. This prevailing model of tumorigenesis has contributed to the clinical notion that pancreatic cancer evolves slowly and presents at a late stage. However, the propensity for this disease to rapidly metastasize and the inability to improve patient outcomes, despite efforts aimed at early detection, suggest that pancreatic cancer progression is not gradual. Here, using newly developed informatics tools, we tracked changes in DNA copy number and their associated rearrangements in tumour-enriched genomes and found that pancreatic cancer tumorigenesis is neither gradual nor follows the accepted mutation order. Two-thirds of tumours harbour complex rearrangement patterns associated with mitotic errors, consistent with punctuated equilibrium as the principal evolutionary trajectory. In a subset of cases, the consequence of such errors is the simultaneous, rather than sequential, knockout of canonical preneoplastic genetic drivers that are likely to set-off invasive cancer growth. These findings challenge the current progression model of pancreatic cancer and provide insights into the mutational processes that give rise to these aggressive tumours.
Gene set enrichment (GSE) analysis allows researchers to efficiently extract biological insight from long lists of differentially expressed genes by interrogating them at a systems level. In recent ...years, there has been a proliferation of GSE analysis methods and hence it has become increasingly difficult for researchers to select an optimal GSE tool based on their particular dataset. Moreover, the majority of GSE analysis methods do not allow researchers to simultaneously compare gene set level results between multiple experimental conditions.
The ensemble of genes set enrichment analyses (EGSEA) is a method developed for RNA-sequencing data that combines results from twelve algorithms and calculates collective gene set scores to improve the biological relevance of the highest ranked gene sets. EGSEA's gene set database contains around 25 000 gene sets from sixteen collections. It has multiple visualization capabilities that allow researchers to view gene sets at various levels of granularity. EGSEA has been tested on simulated data and on a number of human and mouse datasets and, based on biologists' feedback, consistently outperforms the individual tools that have been combined. Our evaluation demonstrates the superiority of the ensemble approach for GSE analysis, and its utility to effectively and efficiently extrapolate biological functions and potential involvement in disease processes from lists of differentially regulated genes.
EGSEA is available as an R package at http://www.bioconductor.org/packages/EGSEA/ . The gene sets collections are available in the R package EGSEAdata from http://www.bioconductor.org/packages/EGSEAdata/ .
monther.alhamdoosh@csl.com.au mritchie@wehi.edu.au.
Supplementary data are available at Bioinformatics online.
The two major subtypes of diffuse large B cell lymphoma (DLBCL)--activated B cell-like (ABC) and germinal center B cell-like (GCB)--arise by distinct mechanisms, with ABC selectively acquiring ...mutations that target the B cell receptor (BCR), fostering chronic active BCR signaling. The ABC subtype has a ∼40% cure rate with currently available therapies, which is worse than the rate for GCB DLBCL, and highlights the need for ABC subtype-specific treatment strategies. We hypothesized that ABC, but not GCB, DLBCL tumors would respond to ibrutinib, an inhibitor of BCR signaling. In a phase 1/2 clinical trial that involved 80 subjects with relapsed or refractory DLBCL, ibrutinib produced complete or partial responses in 37% (14/38) of those with ABC DLBCL, but in only 5% (1/20) of subjects with GCB DLBCL (P = 0.0106). ABC tumors with BCR mutations responded to ibrutinib frequently (5/9; 55.5%), especially those with concomitant myeloid differentiation primary response 88 (MYD88) mutations (4/5; 80%), a result that is consistent with in vitro cooperation between the BCR and MYD88 pathways. However, the highest number of responses occurred in ABC tumors that lacked BCR mutations (9/29; 31%), suggesting that oncogenic BCR signaling in ABC does not require BCR mutations and might be initiated by non-genetic mechanisms. These results support the selective development of ibrutinib for the treatment of ABC DLBCL.
Pancreatic ductal adenocarcinoma (PDAC) has the worst prognosis among solid malignancies and improved therapeutic strategies are needed to improve outcomes. Patient-derived xenografts (PDX) and ...patient-derived organoids (PDO) serve as promising tools to identify new drugs with therapeutic potential in PDAC. For these preclinical disease models to be effective, they should both recapitulate the molecular heterogeneity of PDAC and validate patient-specific therapeutic sensitivities. To date however, deep characterization of the molecular heterogeneity of PDAC PDX and PDO models and comparison with matched human tumour remains largely unaddressed at the whole genome level. We conducted a comprehensive assessment of the genetic landscape of 16 whole-genome pairs of tumours and matched PDX, from primary PDAC and liver metastasis, including a unique cohort of 5 'trios' of matched primary tumour, PDX, and PDO. We developed a pipeline to score concordance between PDAC models and their paired human tumours for genomic events, including mutations, structural variations, and copy number variations. Tumour-model comparisons of mutations displayed single-gene concordance across major PDAC driver genes, but relatively poor agreement across the greater mutational load. Genome-wide and chromosome-centric analysis of structural variation (SV) events highlights previously unrecognized concordance across chromosomes that demonstrate clustered SV events. We found that polyploidy presented a major challenge when assessing copy number changes; however, ploidy-corrected copy number states suggest good agreement between donor-model pairs. Collectively, our investigations highlight that while PDXs and PDOs may serve as tractable and transplantable systems for probing the molecular properties of PDAC, these models may best serve selective analyses across different levels of genomic complexity.
Background
Human leptospirosis mainly affects people in close occupational contact with domestic livestock and their products in New Zealand. The disease has an unquantified impact on both human ...health and animal production in the country. This study aimed to estimate the burden of leptospirosis in terms of disability‐adjusted life years (DALYs) and cost associated with loss due to absence from work, treatment of disease, animal production loss and cost of vaccination.
Methods
Previously published studies of abattoir workers farmers, and veterinarians, reporting annual risks of influenza‐like illness attributable to Leptospira infection, were used to estimate the expected number of cases in a year. The cost of lost animal production was based on results of observational studies in beef cattle, sheep and deer conducted in New Zealand.
Results
Expected median annual number of severe and mild cases of human leptospirosis was 2,025 (95% probability interval 95% PI 1,138–3,422). Median annual DALYs were 0.42 (95% PI: 0.06–2.40) per 100,000 people for the entire population, and 15.82 (95% PI: 2.09–90.80) per 100,000 people working in at‐risk occupations (i.e. abattoir workers, farmers and veterinarians). Human infection resulted in a median cost of 4.42 (95% PI: 2.04–8.62) million US dollars (USD) due to absence from work and disease treatment. Median production loss cost in beef cattle, sheep and deer was USD 7.92 (95% PI: 3.75–15.48) million, while median vaccination cost in cattle, (including dairy), sheep and deer was USD 6.15 (95% PI: 5.30–7.03) million. Total annual cost of leptospirosis plus vaccination was USD 18.80 (95% PI: 13.47–27.15) million, equivalent to USD 440,000 (95% PI: 320,000–640,000) per 100,000 people.
Conclusion
This study provides an estimate of the disease burden and cost of leptospirosis in New Zealand that could support occupational health authorities and livestock industries in assessing interventions for this disease.
Systemic inflammatory scores may aid prognostication and patient selection for trials. We compared five scores in advanced pancreatic adenocarcinoma (PDAC).
Unresectable/metastatic PDAC patients ...enrolled in the Comprehensive Molecular Characterisation of Advanced Pancreatic Ductal Adenocarcinoma for Better Treatment Selection trial (NCT02750657) were included. Patients had pre-treatment biopsies for whole genome and RNA sequencing. CD8 immunohistochemistry was available in a subset. The neutrophil-to-lymphocyte ratio (NLR), platelet-to-lymphocyte ratio, Prognostic Nutritional Index, Gustave Roussy Immune Score (GRIm-S), and Memorial Sloan Kettering Prognostic Score (MPS) were calculated. Overall survival (OS) was estimated using Kaplan-Meier methods. Associations between inflammatory scores, clinical/genomic characteristics, and OS were analysed.
We analysed 263 patients. High-risk NLR, GRIm-S and MPS were poorly prognostic. The GRIm-S had the highest predictive ability: median OS 6.4 vs. 10 months for high risk vs. low-risk (P < 0.001); HR 2.26 (P < 0.001). ECOG ≥ 1, the basal-like subtype, and low-HRDetect were additional poor prognostic factors (P < 0.01). Inflammatory scores did not associate with RNA-based classifiers or homologous recombination repair deficiency genotypes. High-risk MPS (P = 0.04) and GRIm-S (P = 0.02) patients had lower median CD8 + tumour-infiltrating lymphocytes.
Inflammatory scores incorporating NLR have prognostic value in advanced PDAC. Understanding immunophenotypes of poor-risk patients and using these scores in trials will advance the field.