Disambiguate: An open-source application for disambiguating two species in next generation sequencing data from grafted samples

Representation of genomic intratumor heterogeneity in multi-region non-small cell lung cancer patient-derived xenograft models

Patient-derived xenograft (PDX) models are widely used in cancer research. To investigate the genomic fidelity of non-small cell lung cancer PDX models, we established 48 PDX models from 22 patients enrolled in the TRACERx study. Multi-region tumor sampling increased successful PDX engraftment and most models were histologically similar to their parent tumor. Whole-exome sequencing enabled comparison of tumors and PDX models and we provide an adapted mouse reference genome for improved removal of NOD scid gamma (NSG) mouse-derived reads from sequencing data. PDX model establishment caused a genomic bottleneck, with models often representing a single tumor subclone. While distinct tumor subclones were represented in independent models from the same tumor, individual PDX models did not fully recapitulate intratumor heterogeneity. On-going genomic evolution in mice contributed modestly to the genomic distance between tumors and PDX models. Our study highlights the importance of considering primary tumor heterogeneity when using PDX models and emphasizes the benefit of comprehensive tumor sampling.

Alarmin S100A8 imparts chemoresistance of esophageal cancer by reprogramming cancer-associated fibroblasts

Chemotherapy remains the first-line treatment for advanced esophageal cancer. However, durable benefits are achieved by only a limited subset of individuals due to the elusive chemoresistance. Here, we utilize patient-derived xenografts (PDXs) from esophageal squamous-cell carcinoma to investigate chemoresistance mechanisms in preclinical settings. We observe that activated cancer-associated fibroblasts (CAFs) are enriched in the tumor microenvironment of PDXs resistant to chemotherapy. Mechanistically, we reveal that cancer-cell-derived S100A8 triggers the intracellular RhoA-ROCK-MLC2-MRTF-A pathway by binding to the CD147 receptor of CAFs, inducing CAF polarization and leading to chemoresistance. Therapeutically, we demonstrate that blocking the S100A8-CD147 pathway can improve chemotherapy efficiency. Prognostically, we found the S100A8 levels in peripheral blood can serve as an indicator of chemotherapy responsiveness. Collectively, our study offers a comprehensive understanding of the molecular mechanisms underlying chemoresistance in esophageal cancer and highlights the potential value of S100A8 in the clinical management of esophageal cancer.

Combining the Tyrosine Kinase Inhibitor Cabozantinib and the mTORC1/2 Inhibitor Sapanisertib Blocks ERK Pathway Activity and Suppresses Tumor Growth in Renal Cell Carcinoma

Current treatment approaches for renal cell carcinoma (RCC) face challenges in achieving durable tumor responses due to tumor heterogeneity and drug resistance. Combination therapies that leverage tumor molecular profiles could offer an avenue for enhancing treatment efficacy and addressing the limitations of current therapies. To identify effective strategies for treating RCC, we selected ten drugs guided by tumor biology to test in six RCC patient-derived xenograft (PDX) models. The multitargeted tyrosine kinase inhibitor (TKI) cabozantinib and mTORC1/2 inhibitor sapanisertib emerged as the most effective drugs, particularly when combined. The combination demonstrated favorable tolerability and inhibited tumor growth or induced tumor regression in all models, including two from patients who experienced treatment failure with FDA-approved TKI and immunotherapy combinations. In cabozantinib-treated samples, imaging analysis revealed a significant reduction in vascular density, and single-nucleus RNA sequencing (snRNA-seq) analysis indicated a decreased proportion of endothelial cells in the tumors. SnRNA-seq data further identified a tumor subpopulation enriched with cell-cycle activity that exhibited heightened sensitivity to the cabozantinib and sapanisertib combination. Conversely, activation of the epithelial-mesenchymal transition pathway, detected at the protein level, was associated with drug resistance in residual tumors following combination treatment. The combination effectively restrained ERK phosphorylation and reduced expression of ERK downstream transcription factors and their target genes implicated in cell-cycle control and apoptosis. This study highlights the potential of the cabozantinib plus sapanisertib combination as a promising treatment approach for patients with RCC, particularly those whose tumors progressed on immune checkpoint inhibitors and other TKIs.

SIGNIFICANCE

The molecular-guided therapeutic strategy of combining cabozantinib and sapanisertib restrains ERK activity to effectively suppress growth of renal cell carcinomas, including those unresponsive to immune checkpoint inhibitors.

Genomic profiling of subcutaneous patient-derived xenografts reveals immune constraints on tumor evolution in childhood solid cancer

Subcutaneous patient-derived xenografts (PDXs) are an important tool for childhood cancer research. Here, we describe a resource of 68 early passage PDXs established from 65 pediatric solid tumor patients. Through genomic profiling of paired PDXs and patient tumors (PTs), we observe low mutational similarity in about 30% of the PT/PDX pairs. Clonal analysis in these pairs show an aggressive PT minor subclone seeds the major clone in the PDX. We show evidence that this subclone is more immunogenic and is likely suppressed by immune responses in the PT. These results suggest interplay between intratumoral heterogeneity and antitumor immunity may underlie the genetic disparity between PTs and PDXs. We further show that PDXs generally recapitulate PTs in copy number and transcriptomic profiles. Finally, we report a gene fusion LRPAP1-PDGFRA. In summary, we report a childhood cancer PDX resource and our study highlights the role of immune constraints on tumor evolution.

Preclinical evidence for the effective use of TL-895, a highly selective and potent second-generation BTK inhibitor, for the treatment of B-cell malignancies

TL-895 (formerly known as M7583) is a potent, highly selective, adenosine triphosphate (ATP)-competitive, second-generation, irreversible inhibitor of Bruton's tyrosine kinase (BTK). We characterized its biochemical and cellular effects in in vitro and in vivo models. TL-895 was evaluated preclinically for potency against BTK using IC50 concentration-response curves; selectivity using a 270-kinase panel; BTK phosphorylation in Ramos Burkitt's lymphoma cells by ProteinSimple Wes analysis of one study; anti-proliferative effects in primary chronic lymphocytic leukemia (CLL) blasts; cell viability effects in diffuse large B-cell lymphoma (DLBCL) and mantle-cell lymphoma (MCL) cell lines; effects on antibody-dependent cell-mediated cytotoxicity (ADCC) from Daudi cells and chromium-51 release from human tumor cell lines; and efficacy in vivo using four MCL xenograft model and 21 DLBCL patient-derived xenograft (PDX) models (subtypes: 9 ABC, 11 GCB, 1 Unclassified). TL-895 was active against recombinant BTK (average IC50 1.5 nM) and inhibited only three additional kinases with IC50 within tenfold of BTK activity. TL-895 inhibited BTK auto-phosphorylation at the Y223 phosphorylation site (IC50 1-10 nM). TL-895 inhibited the proliferation of primary CLL blasts in vitro and inhibited growth in a subset of activated DLBCL and MCL cell lines. TL-895 inhibited the ADCC mechanism of therapeutic antibodies only at supra-clinical exposure levels. TL-895 significantly inhibited tumor growth in the Mino MCL xenograft model and in 5/21 DLBCL PDX models relative to vehicle controls. These findings demonstrate the potency of TL-895 for BTK and its efficacy in models of B-cell lymphoma despite its refined selectivity.

BET inhibition targets ABC-DLBCL constitutive B-cell receptor signaling through PAX5

B-cell receptor (BCR) signaling is essential for the diffuse large B-cell lymphoma (DLBCL) subtype that originates from activated B-cells (ABCs). ABC-DLBCL cells are sensitive to Bruton tyrosine kinase intervention. However, patients with relapsed or refractory ABC-DLBCL had overall response rates from 33% to 37% for Bruton tyrosine kinase inhibitors, suggesting the evaluation of combination-based treatment for improved efficacy. We investigated the efficacy and mechanism of the bromodomain and extraterminal motif (BET) inhibitor AZD5153 combined with the Bruton tyrosine kinase inhibitor acalabrutinib in ABC-DLBCL preclinical models. AZD5153 is a bivalent BET inhibitor that simultaneously engages the 2 bromodomains of BRD4. Adding AZD5153 to acalabrutinib demonstrated combination benefits in ABC-DLBCL cell line and patient-derived xenograft models. Differential expression analyses revealed PAX5 transcriptional activity as a novel downstream effector of this drug combination. PAX5 is a transcription factor for BCR signaling genes and may be critical for perpetually active BCR signaling in ABC-DLBCL. Our analyses further indicated significant alterations in BCR, RELB/alternative NF-κB, and toll-like receptor/interferon signaling. Validation of these results mapped a positive-feedback signaling loop regulated by PAX5. We demonstrated that AZD5153 decreased PAX5 expression, whereas acalabrutinib disruption of BCR signaling inhibited PAX5 activation. Furthermore, several interferon levels were decreased by AZD5153 and acalabrutinib in tumors. Adding interferon-beta1 (IFNβ1) to cells treated with acalabrutinib partially rescued PAX5 activation. Our results demonstrate that AZD5153 enhances the efficacy of acalabrutinib through PAX5 and BCR mechanisms that are critical for ABC-DLBCL.

Genotoxic effects of base and prime editing in human hematopoietic stem cells

Base and prime editors (BEs and PEs) may provide more precise genetic engineering than nuclease-based approaches because they bypass the dependence on DNA double-strand breaks. However, little is known about their cellular responses and genotoxicity. Here, we compared state-of-the-art BEs and PEs and Cas9 in human hematopoietic stem and progenitor cells with respect to editing efficiency, cytotoxicity, transcriptomic changes and on-target and genome-wide genotoxicity. BEs and PEs induced detrimental transcriptional responses that reduced editing efficiency and hematopoietic repopulation in xenotransplants and also generated DNA double-strand breaks and genotoxic byproducts, including deletions and translocations, at a lower frequency than Cas9. These effects were strongest for cytidine BEs due to suboptimal inhibition of base excision repair and were mitigated by tailoring delivery timing and editor expression through optimized mRNA design. However, BEs altered the mutational landscape of hematopoietic stem and progenitor cells across the genome by increasing the load and relative proportions of nucleotide variants. These findings raise concerns about the genotoxicity of BEs and PEs and warrant further investigation in view of their clinical application.

Toward Practical Integration of Omic and Imaging Data in Co-Clinical Trials

Co-clinical trials are the concurrent or sequential evaluation of therapeutics in both patients clinically and patient-derived xenografts (PDX) pre-clinically, in a manner designed to match the pharmacokinetics and pharmacodynamics of the agent(s) used. The primary goal is to determine the degree to which PDX cohort responses recapitulate patient cohort responses at the phenotypic and molecular levels, such that pre-clinical and clinical trials can inform one another. A major issue is how to manage, integrate, and analyze the abundance of data generated across both spatial and temporal scales, as well as across species. To address this issue, we are developing MIRACCL (molecular and imaging response analysis of co-clinical trials), a web-based analytical tool. For prototyping, we simulated data for a co-clinical trial in "triple-negative" breast cancer (TNBC) by pairing pre- (T0) and on-treatment (T1) magnetic resonance imaging (MRI) from the I-SPY2 trial, as well as PDX-based T0 and T1 MRI. Baseline (T0) and on-treatment (T1) RNA expression data were also simulated for TNBC and PDX. Image features derived from both datasets were cross-referenced to omic data to evaluate MIRACCL functionality for correlating and displaying MRI-based changes in tumor size, vascularity, and cellularity with changes in mRNA expression as a function of treatment.

Homologous recombination deficiency derived from whole-genome sequencing predicts platinum response in triple-negative breast cancers

The high frequency of homologous recombination deficiency (HRD) is the main rationale of testing platinum-based chemotherapy in triple-negative breast cancer (TNBC), however, the existing methods to identify HRD are controversial and there is a medical need for predictive biomarkers. We assess the in vivo response to platinum agents in 55 patient-derived xenografts (PDX) of TNBC to identify determinants of response. The HRD status, determined from whole genome sequencing, is highly predictive of platinum response. BRCA1 promoter methylation is not associated with response, in part due to residual BRCA1 gene expression and homologous recombination proficiency in different tumours showing mono-allelic methylation. Finally, in 2 cisplatin sensitive tumours we identify mutations in XRCC3 and ORC1 genes that are functionally validated in vitro. In conclusion, our results demonstrate that the genomic HRD is predictive of platinum response in a large cohort of TNBC PDX and identify alterations in XRCC3 and ORC1 genes driving cisplatin response.

Therapeutic implications of transcriptomics in head and neck cancer patient-derived xenografts

There are currently no clinical strategies utilizing tumor gene expression to inform therapeutic selection for patients with head and neck squamous cell carcinoma (HNSCC). One of the challenges in developing predictive biomarkers is the limited characterization of preclinical HNSCC models. Patient-derived xenografts (PDXs) are increasingly recognized as translationally relevant preclinical avatars for human tumors; however, the overall transcriptomic concordance of HNSCC PDXs with primary human HNSCC is understudied, especially in human papillomavirus-associated (HPV+) disease. Here, we characterized 64 HNSCC PDXs (16 HPV+ and 48 HPV-) at the transcriptomic level using RNA-sequencing. The range of human-specific reads per PDX varied from 64.6%-96.5%, with a comparison of the most differentially expressed genes before and after removal of mouse transcripts revealing no significant benefit to filtering out mouse mRNA reads in this cohort. We demonstrate that four previously established HNSCC molecular subtypes found in The Cancer Genome Atlas (TCGA) are also clearly recapitulated in HNSCC PDXs. Unsupervised hierarchical clustering yielded a striking natural division of HNSCC PDXs by HPV status, with C19orf57 (BRME1), a gene previously correlated with positive response to cisplatin in cervical cancer, among the most significantly differentially expressed genes between HPV+ and HPV- PDXs. In vivo experiments demonstrated a possible relationship between increased C19orf57 expression and superior anti-tumor responses of PDXs to cisplatin, which should be investigated further. These findings highlight the value of PDXs as models for HPV+ and HPV- HNSCC, providing a resource for future discovery of predictive biomarkers to guide treatment selection in HNSCC.

Patient-Derived Xenografts and Organoids Recapitulate Castration-Resistant Prostate Cancer with Sustained Androgen Receptor Signaling

Castration-resistant prostate cancer (CRPC) remains an incurable and lethal malignancy. The development of new CRPC treatment strategies is strongly impeded by the scarcity of representative, scalable and transferable preclinical models of advanced, androgen receptor (AR)-driven CRPC. Here, we present contemporary patient-derived xenografts (PDXs) and matching PDX-derived organoids (PDXOs) from CRPC patients who had undergone multiple lines of treatment. These models were comprehensively profiled at the morphologic, genomic (n = 8) and transcriptomic levels (n = 81). All are high-grade adenocarcinomas that exhibit copy number alterations and transcriptomic features representative of CRPC patient cohorts. We identified losses of PTEN and RB1, MYC amplifications, as well as genomic alterations in TP53 and in members of clinically actionable pathways such as AR, PI3K and DNA repair pathways. Importantly, the clinically observed continued reliance of CRPC tumors on AR signaling is preserved across the entire set of models, with AR amplification identified in four PDXs. We demonstrate that PDXs and PDXOs faithfully reflect donor tumors and mimic matching patient drug responses. In particular, our models predicted patient responses to subsequent treatments and captured sensitivities to previously received therapies. Collectively, these PDX-PDXO pairs constitute a reliable new resource for in-depth studies of treatment-induced, AR-driven resistance mechanisms. Moreover, PDXOs can be leveraged for large-scale tumor-specific drug response profiling critical for accelerating therapeutic advances in CRPC.

Pharmacogenomic Profiling of Pediatric Acute Myeloid Leukemia to Identify Therapeutic Vulnerabilities and Inform Functional Precision Medicine

Despite the expanding portfolio of targeted therapies for adults with acute myeloid leukemia (AML), direct implementation in children is challenging due to inherent differences in underlying genetics. Here we established the pharmacologic profile of pediatric AML by screening myeloblast sensitivity to approved and investigational agents, revealing candidates of immediate clinical relevance. Drug responses ex vivo correlated with patient characteristics, exhibited age-specific alterations, and concorded with activities in xenograft models. Integration with genomic data uncovered new gene-drug associations, suggesting actionable therapeutic vulnerabilities. Transcriptome profiling further identified gene-expression signatures associated with on- and off-target drug responses. We also demonstrated the feasibility of drug screening-guided treatment for children with high-risk AML, with two evaluable cases achieving remission. Collectively, this study offers a high-dimensional gene-drug clinical data set that could be leveraged to research the unique biology of pediatric AML and sets the stage for realizing functional precision medicine for the clinical management of the disease.

SIGNIFICANCE

We conducted integrated drug and genomic profiling of patient biopsies to build the functional genomic landscape of pediatric AML. Age-specific differences in drug response and new gene-drug interactions were identified. The feasibility of functional precision medicine-guided management of children with high-risk AML was successfully demonstrated in two evaluable clinical cases. This article is highlighted in the In This Issue feature, p. 476.

Genomic and Molecular Signatures of Successful Patient-Derived Xenografts for Oral Cavity Squamous Cell Carcinoma

Background

Oral cavity squamous cell carcinoma (OSCC) is an aggressive malignant tumor with high recurrence and poor prognosis in the advanced stage. Patient-derived xenografts (PDXs) serve as powerful preclinical platforms for drug testing and precision medicine for cancer therapy. We assess which molecular signatures affect tumor engraftment ability and tumor growth rate in OSCC PDXs.

Methods

Treatment-naïve OSCC primary tumors were collected for PDX models establishment. Comprehensive genomic analysis, including whole-exome sequencing and RNA-seq, was performed on case-matched tumors and PDXs. Regulatory genes/pathways were analyzed to clarify which molecular signatures affect tumor engraftment ability and the tumor growth rate in OSCC PDXs.

Results

Perineural invasion was found as an important pathological feature related to engraftment ability. Tumor microenvironment with enriched hypoxia, PI3K-Akt, and epithelial–mesenchymal transition pathways and decreased inflammatory responses had high engraftment ability and tumor growth rates in OSCC PDXs. High matrix metalloproteinase-1 (MMP1) expression was found that have a great graft advantage in xenografts and is associated with pooled disease-free survival in cancer patients.

Conclusion

This study provides a panel with detailed genomic characteristics of OSCC PDXs, enabling preclinical studies on personalized therapy options for oral cancer. MMP1 could serve as a biomarker for predicting successful xenografts in OSCC patients.

HSF2 cooperates with HSF1 to drive a transcriptional program critical for the malignant state

Heat shock factor 1 (HSF1) is well known for its role in the heat shock response (HSR), where it drives a transcriptional program comprising heat shock protein (HSP) genes, and in tumorigenesis, where it drives a program comprising HSPs and many noncanonical target genes that support malignancy. Here, we find that HSF2, an HSF1 paralog with no substantial role in the HSR, physically and functionally interacts with HSF1 across diverse types of cancer. HSF1 and HSF2 have notably similar chromatin occupancy and regulate a common set of genes that include both HSPs and noncanonical transcriptional targets with roles critical in supporting malignancy. Loss of either HSF1 or HSF2 results in a dysregulated response to nutrient stresses in vitro and reduced tumor progression in cancer cell line xenografts. Together, these findings establish HSF2 as a critical cofactor of HSF1 in driving a cancer cell transcriptional program to support the anabolic malignant state.

An In Vivo CRISPR Screen Identifies Stepwise Genetic Dependencies of Metastatic Progression

Blood-borne metastasis of breast cancer involves a series of tightly regulated sequential steps, including the growth of a primary tumor lesion, intravasation of circulating tumor cells (CTC), and adaptation in various distant metastatic sites. The genes orchestrating each of these steps are poorly understood in physiologically relevant contexts, owing to the rarity of experimental models that faithfully recapitulate the biology, growth kinetics, and tropism of human breast cancer. Here, we conducted an in vivo loss-of-function CRISPR screen in newly derived CTC xenografts, unique in their ability to spontaneously mirror the human disease, and identified specific genetic dependencies for each step of the metastatic process. Validation experiments revealed sensitivities to inhibitors that are already available, such as PLK1 inhibitors, to prevent CTC intravasation. Together, these findings present a new tool to reclassify driver genes involved in the spread of human cancer, providing insights into the biology of metastasis and paving the way to test targeted treatment approaches.

SIGNIFICANCE

A loss-of-function CRISPR screen in human CTC-derived xenografts identifies genes critical for individual steps of the metastatic cascade, suggesting novel drivers and treatment opportunities for metastatic breast cancers.

Continued Androgen Signalling Inhibition improves Cabazitaxel Efficacy in Prostate Cancer

BACKGROUND

The androgen receptor (AR) pathway is a key driver of neoplastic behaviour in the different stages of metastatic prostate cancer (mPCa). Targeting the AR therefore remains the cornerstone for mPCa treatment. We have previously reported that activation of AR signalling affects taxane chemo-sensitivity in preclinical models of castration resistant PCa (CRPC). Here, we explored the anti-tumour efficacy of the AR targeted inhibitor enzalutamide combined with cabazitaxel.

METHODS

We used the AR positive CRPC model PC346C-DCC-K to assess the in vitro and in vivo activity of combining enzalutamide with cabazitaxel. Subsequent validation studies were performed using an enzalutamide resistant VCaP model. To investigate the impact of AR signalling on cabazitaxel activity we used quantitative live-cell imaging of tubulin stabilization and apoptosis related nuclear fragmentation.

FINDINGS

Enzalutamide strongly amplified cabazitaxel anti-tumour activity in the patient-derived xenograft models PC346C-DCC-K (median time to humane endpoint 77 versus 48 days, P<0.0001) and VCaP-Enza-B (median time to humane endpoint 80 versus 53 days, P<0.001). Although enzalutamide treatment by itself was ineffective in reducing tumour growth, it significantly suppressed AR signalling in PC346C-DCC-K tumours as shown by AR target gene expression. The addition of enzalutamide enhanced cabazitaxel induced apoptosis as shown by live-cell imaging (P<0.001).

INTERPRETATION

Our study demonstrates that cabazitaxel efficacy can be improved by simultaneous blocking of AR signalling by enzalutamide, even if AR targeted treatment no longer affects tumour growth. These findings support clinical studies that combine AR targeted inhibitors with cabazitaxel in CRPC.

Intussusceptive Angiogenesis in Human Metastatic Malignant Melanoma

Angiogenesis supplies oxygen and nutrients to growing tumors. Inhibiting angiogenesis may stop tumor growth, but vascular endothelial growth factor inhibitors have limited effect in most tumors. This limited effect may be explained by an additional, less vascular endothelial growth factor-driven form of angiogenesis known as intussusceptive angiogenesis. The importance of intussusceptive angiogenesis in human tumors is not known. Epifluorescence and confocal microscopy was used to visualize intravascular pillars, the hallmark structure of intussusceptive angiogenesis, in tumors. Human malignant melanoma metastases, patient-derived melanoma xenografts in mice (PDX), and genetically engineered v-raf murine sarcoma viral oncogene homolog B1 (BRAF)-induced, phosphatase and TENsin homolog deleted on chromosome 10 (PTEN)-deficient (BPT) mice (BrafCA/+Ptenf/fTyr-Cre+/0-mice) were analyzed for pillars. Gene expression in human melanoma metastases and PDXs was analyzed by RNA sequencing. Matrix metalloproteinase 9 (MMP9) protein expression and T-cell and macrophage infiltration in tumor sections were determined with multiplex immunostaining. Intravascular pillars were detected in human metastases but rarely in PDXs and not in BPT mice. The expression of MMP9 mRNA was higher in human metastases compared with PDXs. High expression of MMP9 protein as well as infiltration of macrophages and T-cells were detected in proximity to intravascular pillars. MMP inhibition blocked formation of pillars, but not tubes or tip cells, in vitro. In conclusion, intussusceptive angiogenesis may contribute to the growth of human melanoma metastases. MMP inhibition blocked pillar formation in vitro and should be further investigated as a potential anti-angiogenic drug target in metastatic melanoma.

Malat-1-PRC2-EZH1 interaction supports adaptive oxidative stress dependent epigenome remodeling in skeletal myotubes

PRC2-mediated epigenetic function involves the interaction with long non-coding RNAs (lncRNAs). Although the identity of some of these RNAs has been elucidated in the context of developmental programs, their counterparts in postmitotic adult tissue homeostasis remain uncharacterized. To this aim, we used terminally differentiated postmitotic skeletal muscle cells in which oxidative stress induces the dynamic activation of PRC2-Ezh1 through Embryonic Ectoderm Develpment (EED) shuttling to the nucleus. We identify lncRNA Malat-1 as a necessary partner for PRC2-Ezh1-dependent response to oxidative stress. We show that in this pathway, PRC2-EZH1 dynamic assembly, and in turn stress induced skeletal muscle targeted genes repression, depends specifically on Malat-1. Our study reports about PRC2-RNA interactions in the physiological context of adaptive oxidative stress response and identifies the first lncRNA involved in PRC2-Ezh1 function.

Comprehensive characterization of 536 patient-derived xenograft models prioritizes candidatesfor targeted treatment

Development of candidate cancer treatments is a resource-intensive process, with the research community continuing to investigate options beyond static genomic characterization. Toward this goal, we have established the genomic landscapes of 536 patient-derived xenograft (PDX) models across 25 cancer types, together with mutation, copy number, fusion, transcriptomic profiles, and NCI-MATCH arms. Compared with human tumors, PDXs typically have higher purity and fit to investigate dynamic driver events and molecular properties via multiple time points from same case PDXs. Here, we report on dynamic genomic landscapes and pharmacogenomic associations, including associations between activating oncogenic events and drugs, correlations between whole-genome duplications and subclone events, and the potential PDX models for NCI-MATCH trials. Lastly, we provide a web portal having comprehensive pan-cancer PDX genomic profiles and source code to facilitate identification of more druggable events and further insights into PDXs' recapitulation of human tumors.

Mechanisms of responsiveness to and resistance against trabectedin in murine models of human myxoid liposarcoma

Myxoid liposarcoma (MLPS) is a rare soft-tissue sarcoma characterised by the expression of FUS-DDIT3 chimera. Trabectedin has shown significant clinical anti-tumour activity against MLPS. To characterise the molecular mechanism of trabectedin sensitivity and of resistance against it, we integrated genomic and transcriptomic data from treated mice bearing ML017 or ML017/ET, two patient-derived MLPS xenograft models, sensitive to and resistant against trabectedin, respectively. Longitudinal RNA-Seq analysis of ML017 showed that trabectedin acts mainly as a transcriptional regulator: 15 days after the third dose trabectedin modulates the transcription of 4883 genes involved in processes that sustain adipocyte differentiation. No such differences were observed in ML017/ET. Genomic analysis showed that prolonged treatment causes losses in 4p15.2, 4p16.3 and 17q21.3 cytobands leading to acquired-resistance against the drug. The results dissect the complex mechanism of action of trabectedin and provide the basis for novel combinatorial approaches for the treatment of MLPS that could overcome drug-resistance.

Genome sequence, transcriptome, and annotation of rodent malaria parasite Plasmodium yoelii nigeriensis N67

Background

Rodent malaria parasites are important models for studying host-malaria parasite interactions such as host immune response, mechanisms of parasite evasion of host killing, and vaccine development. One of the rodent malaria parasites is Plasmodium yoelii, and multiple P. yoelii strains or subspecies that cause different disease phenotypes have been widely employed in various studies. The genomes and transcriptomes of several P. yoelii strains have been analyzed and annotated, including the lethal strains of P. y. yoelii YM (or 17XL) and non-lethal strains of P. y. yoelii 17XNL/17X. Genomic DNA sequences and cDNA reads from another subspecies P. y. nigeriensis N67 have been reported for studies of genetic polymorphisms and parasite response to drugs, but its genome has not been assembled and annotated.

Results

We performed genome sequencing of the N67 parasite using the PacBio long-read sequencing technology, de novo assembled its genome and transcriptome, and predicted 5383 genes with high overall annotation quality. Comparison of the annotated genome of the N67 parasite with those of YM and 17X parasites revealed a set of genes with N67-specific orthology, expansion of gene families, particularly the homologs of the Plasmodium chabaudi erythrocyte membrane antigen, large numbers of SNPs and indels, and proteins predicted to interact with host immune responses based on their functional domains.

Conclusions

The genomes of N67 and 17X parasites are highly diverse, having approximately one polymorphic site per 50 base pairs of DNA. The annotated N67 genome and transcriptome provide searchable databases for fast retrieval of genes and proteins, which will greatly facilitate our efforts in studying the parasite biology and gene function and in developing effective control measures against malaria.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07555-9.

Chromatin conformation capture (Hi-C) sequencing of patient-derived xenografts: analysis guidelines

BACKGROUND

Sequencing of patient-derived xenograft (PDX) mouse models allows investigation of the molecular mechanisms of human tumor samples engrafted in a mouse host. Thus, both human and mouse genetic material is sequenced. Several methods have been developed to remove mouse sequencing reads from RNA-seq or exome sequencing PDX data and improve the downstream signal. However, for more recent chromatin conformation capture technologies (Hi-C), the effect of mouse reads remains undefined.

RESULTS

We evaluated the effect of mouse read removal on the quality of Hi-C data using in silico created PDX Hi-C data with 10% and 30% mouse reads. Additionally, we generated 2 experimental PDX Hi-C datasets using different library preparation strategies. We evaluated 3 alignment strategies (Direct, Xenome, Combined) and 3 pipelines (Juicer, HiC-Pro, HiCExplorer) on Hi-C data quality.

CONCLUSIONS

Removal of mouse reads had little-to-no effect on data quality as compared with the results obtained with the Direct alignment strategy. Juicer extracted more valid chromatin interactions for Hi-C matrices, regardless of the mouse read removal strategy. However, the pipeline effect was minimal, while the library preparation strategy had the largest effect on all quality metrics. Together, our study presents comprehensive guidelines on PDX Hi-C data processing.

Fast lightweight accurate xenograft sorting

Motivation

With an increasing number of patient-derived xenograft (PDX) models being created and subsequently sequenced to study tumor heterogeneity and to guide therapy decisions, there is a similarly increasing need for methods to separate reads originating from the graft (human) tumor and reads originating from the host species’ (mouse) surrounding tissue. Two kinds of methods are in use: On the one hand, alignment-based tools require that reads are mapped and aligned (by an external mapper/aligner) to the host and graft genomes separately first; the tool itself then processes the resulting alignments and quality metrics (typically BAM files) to assign each read or read pair. On the other hand, alignment-free tools work directly on the raw read data (typically FASTQ files). Recent studies compare different approaches and tools, with varying results.

Results

We show that alignment-free methods for xenograft sorting are superior concerning CPU time usage and equivalent in accuracy. We improve upon the state of the art sorting by presenting a fast lightweight approach based on three-way bucketed quotiented Cuckoo hashing. Our hash table requires memory comparable to an FM index typically used for read alignment and less than other alignment-free approaches. It allows extremely fast lookups and uses less CPU time than other alignment-free methods and alignment-based methods at similar accuracy. Several engineering steps (e.g., shortcuts for unsuccessful lookups, software prefetching) improve the performance even further.

Availability

Our software xengsort is available under the MIT license at http://gitlab.com/genomeinformatics/xengsort. It is written in numba-compiled Python and comes with sample Snakemake workflows for hash table construction and dataset processing.

The Promise of Patient-Derived Preclinical Models to Accelerate the Implementation of Personalised Medicine for Children with Neuroblastoma

Patient-derived preclinical models are now a core component of cancer research and have the ability to drastically improve the predictive power of preclinical therapeutic studies. However, their development and maintenance can be challenging, time consuming, and expensive. For neuroblastoma, a developmental malignancy of the neural crest, it is possible to establish patient-derived models as xenografts in mice and zebrafish, and as spheroids and organoids in vitro. These varied approaches have contributed to comprehensive packages of preclinical evidence in support of new therapeutics for neuroblastoma. We discuss here the ethical and technical considerations for the creation of patient-derived models of neuroblastoma and how their use can be optimized for the study of tumour evolution and preclinical therapies. We also discuss how neuroblastoma patient-derived models might become avatars for personalised medicine for children with this devastating disease.

Chromosome 8 gain is associated with high-grade transformation in MPNST

One of the most common malignancies affecting adults with Neurofibromatosis type 1 (NF1) is the malignant peripheral nerve sheath tumor (MPNST), an aggressive and often fatal sarcoma that commonly arises from benign plexiform neurofibromas. Despite advances in our understanding of MPNST pathobiology, there are few effective therapeutic options, and no investigational agents have proven successful in clinical trials. To further understand the genomic heterogeneity of MPNST, and to generate a preclinical platform that encompasses this heterogeneity, we developed a collection of NF1-MPNST patient-derived xenografts (PDX). These PDX were compared with the primary tumors from which they were derived using copy number analysis, whole exome sequencing, and RNA sequencing. We identified chromosome 8 gain as a recurrent genomic event in MPNST and validated its occurrence by FISH in the PDX and parental tumors, in a validation cohort, and by single-cell sequencing in the PDX. Finally, we show that chromosome 8 gain is associated with inferior overall survival in soft-tissue sarcomas. These data suggest that chromosome 8 gain is a critical event in MPNST pathogenesis and may account for the aggressive nature and poor outcomes in this sarcoma subtype.

Profiling of conditionally reprogrammed cell lines for in vitro chemotherapy response prediction of pancreatic cancer

BACKGROUND

The establishment of patient-derived models for pancreatic ductal adenocarcinoma (PDAC) using conventional methods has been fraught with low success rate, mainly because of the small number of tumour cells and dense fibrotic stroma. Here, we sought to establish patient-derived model of PDAC and perform genetic analysis with responses to anticancer drug by using the conditionally reprogrammed cell (CRC) methodology.

METHODS

We performed in vitro and in vivo tumourigenicity assays and analysed histological characteristics by immunostaining. We investigated genetic profiles including mutation patterns and copy number variations using targeted deep sequencing and copy-number analyses. We assessed the responses of cultured CRCs to the available clinical anticancer drugs based on patient responsiveness.

FINDINGS

We established a total of 28 CRCs from patients. Of the 28 samples, 27 showed KRAS mutations in codon 12/13 or codon 61. We found that somatic mutations were shared in the primary-CRC pairs and shared mutations included key oncogenic mutations such as KRAS (9 pairs), TP53 (8 pairs), and SMAD4 (3 pairs). Overall, CRCs preserved the genetic characteristics of primary tumours with high concordance, with additional confirmation of low-AF NPM1 mutation in CRC (35 shared mutations out of 36 total, concordance rate=97.2%). CRCs of the responder group were more sensitive to anticancer agents than those of the non-responder group (P < 0.001).

INTERPRETATION

These results show that a pancreatic cancer cell line model can be efficiently established using the CRC methodology, to better support a personalized therapeutic approach for pancreatic cancer patients.

FUNDING

2014R1A1A1006272, HI19C0642-060019, 2019R1A2C2008050, 2020R1A2C209958611, and 2020M3E5E204028211.

Proteomic Resistance Biomarkers for PI3K Inhibitor in Triple Negative Breast Cancer Patient-Derived Xenograft Models

Simple Summary

The objective of this study is to identify potential proteomic biomarkers in triple negative breast cancer (TNBC) that associate with response to PI3K inhibitors which are in clinical trials. We tested a panel of TNBC patient-derived xenograft (PDX) models for their tumor growth response to a pan-PI3K inhibitor, BKM120. Proteomic analyses by reverse phase protein array (RPPA) of 182 markers were performed on baseline and post short-term treatment PDX samples, to correlate with tumor growth response. We identified several baseline and treatment induced proteomic biomarkers in association with resistance. These results provide important insights for the development of PI3K inhibitors in TNBC.

Abstract

PI3K pathway activation is frequently observed in triple negative breast cancer (TNBC). However, single agent PI3K inhibitors have shown limited anti-tumor activity. To investigate biomarkers of response and resistance mechanisms, we tested 17 TNBC patient-derived xenograft (PDX) models representing diverse genomic backgrounds and varying degrees of PI3K pathway signaling activities for their tumor growth response to the pan-PI3K inhibitor, BKM120. Baseline and post-treatment PDX tumors were subjected to reverse phase protein array (RPPA) to identify protein markers associated with tumor growth response. While BKM120 consistently reduced PI3K pathway activity, as demonstrated by reduced levels of phosphorylated AKT, percentage tumor growth inhibition (%TGI) ranged from 35% in the least sensitive to 84% in the most sensitive model. Several biomarkers showed significant association with resistance, including elevated baseline levels of growth factor receptors (EGFR, pHER3 Y1197), PI3Kp85 regulatory subunit, anti-apoptotic protein BclXL, EMT (Vimentin, MMP9, IntegrinaV), NFKB pathway (IkappaB, RANKL), and intracellular signaling molecules including Caveolin, CBP, and KLF4, as well as treatment-induced increases in the levels of phosphorylated forms of Aurora kinases. Interestingly, increased AKT phosphorylation or PTEN loss at baseline were not significantly correlated to %TGI. These results provide important insights into biomarker development for PI3K inhibitors in TNBC.

A Preclinical Trial and Molecularly Annotated Patient Cohort Identify Predictive Biomarkers in Homologous Recombination–deficient Pancreatic Cancer

Purpose:

Pancreatic ductal adenocarcinoma (PDAC) arising in patients with a germline BRCA1 or BRCA2 (gBRCA) mutation may be sensitive to platinum and PARP inhibitors (PARPi). However, treatment stratification based on gBRCA mutational status alone is associated with heterogeneous responses.

Experimental Design:

We performed a seven-arm preclinical trial consisting of 471 mice, representing 12 unique PDAC patient-derived xenografts, of which nine were gBRCA mutated. From 179 patients whose PDAC was whole-genome and transcriptome sequenced, we identified 21 cases with homologous recombination deficiency (HRD), and investigated prognostic biomarkers.

Results:

We found that biallelic inactivation of BRCA1/BRCA2 is associated with genomic hallmarks of HRD and required for cisplatin and talazoparib (PARPi) sensitivity. However, HRD genomic hallmarks persisted in xenografts despite the emergence of therapy resistance, indicating the presence of a genomic scar. We identified tumor polyploidy and a low Ki67 index as predictors of poor cisplatin and talazoparib response. In patients with HRD PDAC, tumor polyploidy and a basal-like transcriptomic subtype were independent predictors of shorter survival. To facilitate clinical assignment of transcriptomic subtype, we developed a novel pragmatic two-marker assay (GATA6:KRT17).

Conclusions:

In summary, we propose a predictive and prognostic model of gBRCA-mutated PDAC on the basis of HRD genomic hallmarks, Ki67 index, tumor ploidy, and transcriptomic subtype.

Neoplastic-Stromal Cell Cross-talk Regulates Matrisome Expression in Pancreatic Cancer

Pancreatic ductal adenocarcinoma (PDAC) is characterized by a highly desmoplastic reaction, warranting intense cancer-stroma communication. In this study, we interrogated the contribution of the BET family of chromatin adaptors to the cross-talk between PDAC cells and the tumor stroma. Short-term treatment of orthotopic xenograft tumors with CPI203, a small-molecule inhibitor of BET proteins, resulted in broad changes in the expression of genes encoding components of the extracellular matrix (matrisome) in both cancer and stromal cells. Remarkably, more than half of matrisome genes were expressed by cancer cells. In vitro cocultures of PDAC cells and cancer-associated fibroblasts (CAF) demonstrated that matrisome expression was regulated by BET-dependent cancer-CAF cross-talk. Disrupting this cross-talk in vivo resulted in diminished growth of orthotopic patient-derived xenograft tumors, reduced proliferation of cancer cells, and changes in collagen structure consistent with that of patients who experienced better survival. Examination of matrisome gene expression in publicly available data sets of 573 PDAC tumors identified a 65-gene signature that was able to distinguish long- and short-term PDAC survivors. Importantly, the expression of genes predictive of short-term survival was diminished in the cancer cells of orthotopic xenograft tumors of mice treated with CPI203. Taken together, these results demonstrate that inhibiting the activity BET proteins results in transcriptional and structural differences in the matrisome are associated with better patient survival. IMPLICATIONS: These studies highlight the biological relevance of the matrisome program in PDAC and suggest targeting of epigenetically driven tumor-stroma cross-talk as a potential therapeutic avenue.

BI-3406, a Potent and Selective SOS1-KRAS Interaction Inhibitor, Is Effective in KRAS-Driven Cancers through Combined MEK Inhibition

KRAS is the most frequently mutated driver of pancreatic, colorectal, and non-small cell lung cancers. Direct KRAS blockade has proved challenging, and inhibition of a key downstream effector pathway, the RAF-MEK-ERK cascade, has shown limited success because of activation of feedback networks that keep the pathway in check. We hypothesized that inhibiting SOS1, a KRAS activator and important feedback node, represents an effective approach to treat KRAS-driven cancers. We report the discovery of a highly potent, selective, and orally bioavailable small-molecule SOS1 inhibitor, BI-3406, that binds to the catalytic domain of SOS1, thereby preventing the interaction with KRAS. BI-3406 reduces formation of GTP-loaded RAS and limits cellular proliferation of a broad range of KRAS-driven cancers. Importantly, BI-3406 attenuates feedback reactivation induced by MEK inhibitors and thereby enhances sensitivity of KRAS-dependent cancers to MEK inhibition. Combined SOS1 and MEK inhibition represents a novel and effective therapeutic concept to address KRAS-driven tumors. SIGNIFICANCE: To date, there are no effective targeted pan-KRAS therapies. In-depth characterization of BI-3406 activity and identification of MEK inhibitors as effective combination partners provide an attractive therapeutic concept for the majority of KRAS-mutant cancers, including those fueled by the most prevalent mutant KRAS oncoproteins, G12D, G12V, G12C, and G13D.See related commentary by Zhao et al., p. 17.This article is highlighted in the In This Issue feature, p. 1.

Squamous trans-differentiation of pancreatic cancer cells promotes stromal inflammation

A highly aggressive subset of pancreatic ductal adenocarcinomas undergo trans-differentiation into the squamous lineage during disease progression. Here, we investigated whether squamous trans-differentiation of human and mouse pancreatic cancer cells can influence the phenotype of non-neoplastic cells in the tumor microenvironment. Conditioned media experiments revealed that squamous pancreatic cancer cells secrete factors that recruit neutrophils and convert pancreatic stellate cells into cancer-associated fibroblasts (CAFs) that express inflammatory cytokines at high levels. We use gain- and loss-of-function approaches to show that squamous-subtype pancreatic tumor models become enriched with neutrophils and inflammatory CAFs in a p63-dependent manner. These effects occur, at least in part, through p63-mediated activation of enhancers at pro-inflammatory cytokine loci, which includes IL1A and CXCL1 as key targets. Taken together, our findings reveal enhanced tissue inflammation as a consequence of squamous trans-differentiation in pancreatic cancer, thus highlighting an instructive role of tumor cell lineage in reprogramming the stromal microenvironment.

Molecular profiling of driver events in metastatic uveal melanoma

Metastatic uveal melanoma is less well understood than its primary counterpart, has a distinct biology compared to skin melanoma, and lacks effective treatments. Here we genomically profile metastatic tumors and infiltrating lymphocytes. BAP1 alterations are overrepresented and found in 29/32 of cases. Reintroducing a functional BAP1 allele into a deficient patient-derived cell line, reveals a broad shift towards a transcriptomic subtype previously associated with better prognosis of the primary disease. One outlier tumor has a high mutational burden associated with UV-damage. CDKN2A deletions also occur, which are rarely present in primaries. A focused knockdown screen is used to investigate overexpressed genes associated withcopy number gains. Tumor-infiltrating lymphocytes are in several cases found tumor-reactive, but expression of the immune checkpoint receptors TIM-3, TIGIT and LAG3 is also abundant. This study represents the largest whole-genome analysis of uveal melanoma to date, and presents an updated view of the metastatic disease.

Cystic renal-epithelial derived induced pluripotent stem cells from polycystic kidney disease patients

Autosomal‐dominant polycystic kidney disease (ADPKD) is the most common inherited kidney disease, leading to kidney failure in most patients. In approximately 85% of cases, the disease is caused by mutations in PKD1. How dysregulation of PKD1 leads to cyst formation on a molecular level is unknown. Induced pluripotent stem cells (iPSCs) are a powerful tool for in vitro modeling of genetic disorders. Here, we established ADPKD patient‐specific iPSCs to study the function of PKD1 in kidney development and cyst formation in vitro. Somatic mutations are proposed to be the initiating event of cyst formation, and therefore, iPSCs were derived from cystic renal epithelial cells rather than fibroblasts. Mutation analysis of the ADPKD iPSCs revealed germline mutations in PKD1 but no additional somatic mutations in PKD1/PKD2. Although several somatic mutations in other genes implicated in ADPKD were identified in cystic renal epithelial cells, only few of these mutations were present in iPSCs, indicating a heterogeneous mutational landscape, and possibly in vitro cell selection before and during the reprogramming process. Whole‐genome DNA methylation analysis indicated that iPSCs derived from renal epithelial cells maintain a kidney‐specific DNA methylation memory. In addition, comparison of PKD1+/− and control iPSCs revealed differences in DNA methylation associated with the disease history. In conclusion, we generated and characterized iPSCs derived from cystic and healthy control renal epithelial cells, which can be used for in vitro modeling of kidney development in general and cystogenesis in particular.

Bacteria Boost Mammalian Host NAD Metabolism by Engaging the Deamidated Biosynthesis Pathway

Nicotinamide adenine dinucleotide (NAD), a cofactor for hundreds of metabolic reactions in all cell types, plays an essential role in metabolism, DNA repair, and aging. However, how NAD metabolism is impacted by the environment remains unclear. Here, we report an unexpected trans-kingdom cooperation between bacteria and mammalian cells wherein bacteria contribute to host NAD biosynthesis. Bacteria confer resistance to inhibitors of NAMPT, the rate-limiting enzyme in the amidated NAD salvage pathway, in cancer cells and xenograft tumors. Mechanistically, a microbial nicotinamidase (PncA) that converts nicotinamide to nicotinic acid, a precursor in the alternative deamidated NAD salvage pathway, is necessary and sufficient for this protective effect. Using stable isotope tracing and microbiota-depleted mice, we demonstrate that this bacteria-mediated deamidation contributes substantially to the NAD-boosting effect of oral nicotinamide and nicotinamide riboside supplementation in several tissues. Collectively, our findings reveal an important role of bacteria-enabled deamidated pathway in host NAD metabolism.

Novel Breast Cancer Brain Metastasis Patient-Derived Orthotopic Xenograft Model for Preclinical Studies

The vast majority of mortality in breast cancer results from distant metastasis. Brain metastases occur in as many as 30% of patients with advanced breast cancer, and the 1-year survival rate of these patients is around 20%. Pre-clinical animal models that reliably reflect the biology of breast cancer brain metastasis are needed to develop and test new treatments for this deadly condition. The patient-derived xenograft (PDX) model maintains many features of a donor tumor, such as intra-tumor heterogeneity, and permits the testing of individualized treatments. However, the establishment of orthotopic PDXs of brain metastasis is procedurally difficult. We have developed a method for generating such PDXs with high tumor engraftment and growth rates. Here, we describe this method and identify variables that affect its outcomes. We also compare the brain-orthotopic PDXs with ectopic PDXs grown in mammary pads of mice, and show that the responsiveness of PDXs to chemotherapeutic reagents can be dramatically affected by the site that they are in.

Clonal competition within complex evolutionary hierarchies shapes AML over time

Clonal heterogeneity and evolution has major implications for disease progression and relapse in acute myeloid leukemia (AML). To model clonal dynamics in vivo, we serially transplanted 23 AML cases to immunodeficient mice and followed clonal composition for up to 15 months by whole-exome sequencing of 84 xenografts across two generations. We demonstrate vast changes in clonality that both progress and reverse over time, and define five patterns of clonal dynamics: Monoclonal, Stable, Loss, Expansion and Burst. We also show that subclonal expansion in vivo correlates with a more adverse prognosis. Furthermore, clonal expansion enabled detection of very rare clones with AML driver mutations that were undetectable by sequencing at diagnosis, demonstrating that the vast majority of AML cases harbor multiple clones already at diagnosis. Finally, the rise and fall of related clones enabled deconstruction of the complex evolutionary hierarchies of the clones that compete to shape AML over time.

Clonal evolution and heterogeneity has strong implications for treatment response in acute myeloid leukemia. Here, the authors use patient derived in vivo modelling to highlight the complex clonal and evolutionary dynamics underpinning acute myeloid leukemia progression.

Impact of mouse contamination in genomic profiling of patient-derived models and best practice for robust analysis

BACKGROUND

Patient-derived xenograft and cell line models are popular models for clinical cancer research. However, the inevitable inclusion of a mouse genome in a patient-derived model is a remaining concern in the analysis. Although multiple tools and filtering strategies have been developed to account for this, research has yet to demonstrate the exact impact of the mouse genome and the optimal use of these tools and filtering strategies in an analysis pipeline.

RESULTS

We construct a benchmark dataset of 5 liver tissues from 3 mouse strains using human whole-exome sequencing kit. Next-generation sequencing reads from mouse tissues are mappable to 49% of the human genome and 409 cancer genes. In total, 1,207,556 mouse-specific alleles are aligned to the human genome reference, including 467,232 (38.7%) alleles with high sensitivity to contamination, which are pervasive causes of false cancer mutations in public databases and are signatures for predicting global contamination. Next, we assess the performance of 8 filtering methods in terms of mouse read filtration and reduction of mouse-specific alleles. All filtering tools generally perform well, although differences in algorithm strictness and efficiency of mouse allele removal are observed. Therefore, we develop a best practice pipeline that contains the estimation of contamination level, mouse read filtration, and variant filtration.

CONCLUSIONS

The inclusion of mouse cells in patient-derived models hinders genomic analysis and should be addressed carefully. Our suggested guidelines improve the robustness and maximize the utility of genomic analysis of these models.

HaTSPiL: A modular pipeline for high-throughput sequencing data analysis

Background

Next generation sequencing methods are widely adopted for a large amount of scientific purposes, from pure research to health-related studies. The decreasing costs per analysis led to big amounts of generated data and to the subsequent improvement of software for the respective analyses. As a consequence, many approaches have been developed to chain different software in order to obtain reliable and reproducible workflows. However, the large range of applications for NGS approaches entails the challenge to manage many different workflows without losing reliability.

Methods

We here present a high-throughput sequencing pipeline (HaTSPiL), a Python-powered CLI tool designed to handle different approaches for data analysis with a high level of reliability. The software relies on the barcoding of filenames using a human readable naming convention that contains any information regarding the sample needed by the software to automatically choose different workflows and parameters. HaTSPiL is highly modular and customisable, allowing the users to extend its features for any specific need.

Conclusions

HaTSPiL is licensed as Free Software under the MIT license and it is available at https://github.com/dodomorandi/hatspil.

Genomic data analysis workflows for tumors from patient-derived xenografts (PDXs): challenges and guidelines

BACKGROUND

Patient-derived xenograft (PDX) models are in vivo models of human cancer that have been used for translational cancer research and therapy selection for individual patients. The Jackson Laboratory (JAX) PDX resource comprises 455 models originating from 34 different primary sites (as of 05/08/2019). The models undergo rigorous quality control and are genomically characterized to identify somatic mutations, copy number alterations, and transcriptional profiles. Bioinformatics workflows for analyzing genomic data obtained from human tumors engrafted in a mouse host (i.e., Patient-Derived Xenografts; PDXs) must address challenges such as discriminating between mouse and human sequence reads and accurately identifying somatic mutations and copy number alterations when paired non-tumor DNA from the patient is not available for comparison.

RESULTS

We report here data analysis workflows and guidelines that address these challenges and achieve reliable identification of somatic mutations, copy number alterations, and transcriptomic profiles of tumors from PDX models that lack genomic data from paired non-tumor tissue for comparison. Our workflows incorporate commonly used software and public databases but are tailored to address the specific challenges of PDX genomics data analysis through parameter tuning and customized data filters and result in improved accuracy for the detection of somatic alterations in PDX models. We also report a gene expression-based classifier that can identify EBV-transformed tumors. We validated our analytical approaches using data simulations and demonstrated the overall concordance of the genomic properties of xenograft tumors with data from primary human tumors in The Cancer Genome Atlas (TCGA).

CONCLUSIONS

The analysis workflows that we have developed to accurately predict somatic profiles of tumors from PDX models that lack normal tissue for comparison enable the identification of the key oncogenic genomic and expression signatures to support model selection and/or biomarker development in therapeutic studies. A reference implementation of our analysis recommendations is available at https://github.com/TheJacksonLaboratory/PDX-Analysis-Workflows .

A RAD51 assay feasible in routine tumor samples calls PARP inhibitor response beyond BRCA mutation

Poly(ADP-ribose) polymerase (PARP) inhibitors (PARPi) are effective in cancers with defective homologous recombination DNA repair (HRR), including BRCA1/2-related cancers. A test to identify additional HRR-deficient tumors will help to extend their use in new indications. We evaluated the activity of the PARPi olaparib in patient-derived tumor xenografts (PDXs) from breast cancer (BC) patients and investigated mechanisms of sensitivity through exome sequencing, BRCA1 promoter methylation analysis, and immunostaining of HRR proteins, including RAD51 nuclear foci. In an independent BC PDX panel, the predictive capacity of the RAD51 score and the homologous recombination deficiency (HRD) score were compared. To examine the clinical feasibility of the RAD51 assay, we scored archival breast tumor samples, including PALB2-related hereditary cancers. The RAD51 score was highly discriminative of PARPi sensitivity versus PARPi resistance in BC PDXs and outperformed the genomic test. In clinical samples, all PALB2-related tumors were classified as HRR-deficient by the RAD51 score. The functional biomarker RAD51 enables the identification of PARPi-sensitive BC and broadens the population who may benefit from this therapy beyond BRCA1/2-related cancers.

BMX-Mediated Regulation of Multiple Tyrosine Kinases Contributes to Castration Resistance in Prostate Cancer

Prostate cancer responds to therapies that suppress androgen receptor (AR) activity (androgen deprivation therapy, ADT) but invariably progresses to castration-resistant prostate cancer (CRPC). The Tec family nonreceptor tyrosine kinase BMX is activated downstream of PI3K and has been implicated in regulation of multiple pathways and in the development of cancers including prostate cancer. However, its precise mechanisms of action, and particularly its endogenous substrates, remain to be established. Here, we demonstrate that BMX expression in prostate cancer is suppressed directly by AR via binding to the BMX gene and that BMX expression is subsequently rapidly increased in response to ADT. BMX contributed to CRPC development in cell line and xenograft models by positively regulating the activities of multiple receptor tyrosine kinases through phosphorylation of a phosphotyrosine-tyrosine (pYY) motif in their activation loop, generating pYpY that is required for full kinase activity. To assess BMX activity in vivo, we generated a BMX substrate-specific antibody (anti-pYpY) and found that its reactivity correlated with BMX expression in clinical samples, supporting pYY as an in vivo substrate. Inhibition of BMX with ibrutinib (developed as an inhibitor of the related Tec kinase BTK) or another BMX inhibitor BMX-IN-1 markedly enhanced the response to castration in a prostate cancer xenograft model. These data indicate that increased BMX in response to ADT contributes to enhanced tyrosine kinase signaling and the subsequent emergence of CRPC, and that combination therapies targeting AR and BMX may be effective in a subset of patients.Significance: The tyrosine kinase BMX is negatively regulated by androgen and contributes to castration-resistant prostate cancer by enhancing the phosphorylation and activation of multiple receptor tyrosine kinases following ADT. Cancer Res; 78(18); 5203-15. ©2018 AACR.

Shieldin complex promotes DNA end-joining and counters homologous recombination in BRCA1-null cells

BRCA1 deficiencies cause breast, ovarian, prostate and other cancers, and render tumours hypersensitive to poly(ADP-ribose) polymerase (PARP) inhibitors. To understand the resistance mechanisms, we conducted whole-genome CRISPR-Cas9 synthetic-viability/resistance screens in BRCA1-deficient breast cancer cells treated with PARP inhibitors. We identified two previously uncharacterized proteins, C20orf196 and FAM35A, whose inactivation confers strong PARP-inhibitor resistance. Mechanistically, we show that C20orf196 and FAM35A form a complex, 'Shieldin' (SHLD1/2), with FAM35A interacting with single-stranded DNA through its C-terminal oligonucleotide/oligosaccharide-binding fold region. We establish that Shieldin acts as the downstream effector of 53BP1/RIF1/MAD2L2 to promote DNA double-strand break (DSB) end-joining by restricting DSB resection and to counteract homologous recombination by antagonizing BRCA2/RAD51 loading in BRCA1-deficient cells. Notably, Shieldin inactivation further sensitizes BRCA1-deficient cells to cisplatin, suggesting how defining the SHLD1/2 status of BRCA1-deficient tumours might aid patient stratification and yield new treatment opportunities. Highlighting this potential, we document reduced SHLD1/2 expression in human breast cancers displaying intrinsic or acquired PARP-inhibitor resistance.

FGFR1 underlies obesity-associated progression of estrogen receptor-positive breast cancer after estrogen deprivation

Obesity increases breast cancer mortality by promoting resistance to therapy. Here, we identified regulatory pathways in estrogen receptor-positive (ER-positive) tumors that were shared between patients with obesity and those with resistance to neoadjuvant aromatase inhibition. Among these was fibroblast growth factor receptor 1 (FGFR1), a known mediator of endocrine therapy resistance. In a preclinical model with patient-derived ER-positive tumors, diet-induced obesity promoted a similar gene expression signature and sustained the growth of FGFR1-overexpressing tumors after estrogen deprivation. Tumor FGFR1 phosphorylation was elevated with obesity and predicted a shorter disease-free and disease-specific survival for patients treated with tamoxifen. In both human and mouse mammary adipose tissue, FGF1 ligand expression was associated with metabolic dysfunction, weight gain, and adipocyte hypertrophy, implicating the impaired response to a positive energy balance in growth factor production within the tumor niche. In conjunction with these studies, we describe a potentially novel graft-competent model that can be used with patient-derived tissue to elucidate factors specific to extrinsic (host) and intrinsic (tumor) tissue that are critical for obesity-associated tumor promotion. Taken together, we demonstrate that obesity and excess energy establish a tumor environment with features of endocrine therapy resistance and identify a role for ligand-dependent FGFR1 signaling in obesity-associated breast cancer progression.

A comparison of next-generation sequencing analysis methods for cancer xenograft samples

The application of next-generation sequencing (NGS) technology in cancer is influenced by the quality and purity of tissue samples. This issue is especially critical for patient-derived xenograft (PDX) models, which have proven to be by far the best preclinical tool for investigating human tumor biology, because the sensitivity and specificity of NGS analysis in xenograft samples would be compromised by the contamination of mouse DNA and RNA. This definitely affects downstream analyses by causing inaccurate mutation calling and gene expression estimates. The reliability of NGS data analysis for cancer xenograft samples is therefore highly dependent on whether the sequencing reads derived from the xenograft could be distinguished from those originated from the host. That is, each sequence read needs to be accurately assigned to its original species. Here, we review currently available methodologies in this field, including Xenome, Disambiguate, bamcmp and pdxBlacklist, and provide guidelines for users.

Computational approach to discriminate human and mouse sequences in patient-derived tumour xenografts

BACKGROUND

Patient-Derived Tumour Xenografts (PDTXs) have emerged as the pre-clinical models that best represent clinical tumour diversity and intra-tumour heterogeneity. The molecular characterization of PDTXs using High-Throughput Sequencing (HTS) is essential; however, the presence of mouse stroma is challenging for HTS data analysis. Indeed, the high homology between the two genomes results in a proportion of mouse reads being mapped as human.

RESULTS

In this study we generated Whole Exome Sequencing (WES), Reduced Representation Bisulfite Sequencing (RRBS) and RNA sequencing (RNA-seq) data from samples with known mixtures of mouse and human DNA or RNA and from a cohort of human breast cancers and their derived PDTXs. We show that using an In silico Combined human-mouse Reference Genome (ICRG) for alignment discriminates between human and mouse reads with up to 99.9% accuracy and decreases the number of false positive somatic mutations caused by misalignment by >99.9%. We also derived a model to estimate the human DNA content in independent PDTX samples. For RNA-seq and RRBS data analysis, the use of the ICRG allows dissecting computationally the transcriptome and methylome of human tumour cells and mouse stroma. In a direct comparison with previously reported approaches, our method showed similar or higher accuracy while requiring significantly less computing time.

CONCLUSIONS

The computational pipeline we describe here is a valuable tool for the molecular analysis of PDTXs as well as any other mixture of DNA or RNA species.

Development of Novel Patient-Derived Xenografts from Breast Cancer Brain Metastases

Brain metastases are an increasing burden among breast cancer patients, particularly for those with HER2+ and triple negative (TN) subtypes. Mechanistic insight into the pathophysiology of brain metastases and preclinical validation of therapies has relied almost exclusively on intracardiac injection of brain-homing cells derived from highly aggressive TN MDA-MB-231 and HER2+ BT474 breast cancer cell lines. Yet, these well characterized models are far from representing the tumor heterogeneity observed clinically and, due to their fast progression in vivo, their suitability to validate therapies for established brain metastasis remains limited. The goal of this study was to develop and characterize novel human brain metastasis breast cancer patient-derived xenografts (BM-PDXs) to study the biology of brain metastasis and to serve as tools for testing novel therapeutic approaches. We obtained freshly resected brain metastases from consenting donors with breast cancer. Tissue was immediately implanted in the mammary fat pad of female immunocompromised mice and expanded as BM-PDXs. Brain metastases from 3/4 (75%) TN, 1/1 (100%) estrogen receptor positive (ER+), and 5/9 (55.5%) HER2+ clinical subtypes were established as transplantable BM-PDXs. To facilitate tracking of metastatic dissemination using BM-PDXs, we labeled PDX-dissociated cells with EGFP-luciferase followed by reimplantation in mice, and generated a BM-derived cell line (F2-7). Immunohistologic analyses demonstrated that parental and labeled BM-PDXs retained expression of critical clinical markers such as ER, progesterone receptor, epidermal growth factor receptor, HER2, and the basal cell marker cytokeratin 5. Similarly, RNA sequencing analysis showed clustering of parental, labeled BM-PDXs and their corresponding cell line derivative. Intracardiac injection of dissociated cells from BM-E22-1, resulted in magnetic resonance imaging-detectable macrometastases in 4/8 (50%) and micrometastases (8/8) (100%) mice, suggesting that BM-PDXs remain capable of colonizing the brain at high frequencies. Brain metastases developed 8-12 weeks after ic injection, located to the brain parenchyma, grew around blood vessels, and elicited astroglia activation characteristic of breast cancer brain metastasis. These novel BM-PDXs represent heterogeneous and clinically relevant models to study mechanisms of brain metastatic colonization, with the added benefit of a slower progression rate that makes them suitable for preclinical testing of drugs in therapeutic settings.