Dealing With Variable Drug Exposure Due to Variable Hepatic Metabolism: A Proof-of-Concept Application of Liquid Biopsy in Renal Impairment.

Precision dosing strategies require accounting for between-patient variability in pharmacokinetics (PK), affecting drug exposure, and in pharmacodynamics (PD), affecting response achieved at the same drug concentration at the site of action. Although liquid biopsy for assessing different levels of molecular drug targets has yet to be established, individual characterization of drug elimination pathways using liquid biopsy has recently been demonstrated. The feasibility of applying this approach in conjunction with modeling tools to guide individual dosing remains unexplored. In this study, we aimed to individualize physiologically-based pharmacokinetic (PBPK) models based on liquid biopsy measurements in plasma from 25 donors with different grades of renal function who were previously administered oral midazolam as part of a microdose cocktail. Virtual twin models were constructed based on demographics, renal function, and hepatic expression of relevant pharmacokinetic pathways projected from liquid biopsy output. Simulated exposure (AUC) to midazolam was in agreement with observed data (AFE = 1.38, AAFE = 1.78). Simulated AUC variability with three dosing approaches indicated higher variability with uniform dosing (14-fold) and stratified dosing (13-fold) compared with individualized dosing informed by liquid biopsy (fivefold). Further, exosome screening revealed mRNA expression of 532 targets relevant to drug metabolism and disposition (169 enzymes and 361 transporters). Data related to these targets can be used to further individualize PBPK models for pathways relevant to PK of other drugs. This study provides additional verification of liquid biopsy-informed PBPK modeling approaches, necessary to advance strategies that seek to achieve precise dosing from the start of treatment.

Hypoxia Is Associated with Increased Immune Infiltrates and Both Anti-Tumour and Immune Suppressive Signalling in Muscle-Invasive Bladder Cancer.

Hypoxia and a suppressive tumour microenvironment (TME) are both independent negative prognostic factors for muscle-invasive bladder cancer (MIBC) that contribute to treatment resistance. Hypoxia has been shown to induce an immune suppressive TME by recruiting myeloid cells that inhibit anti-tumour T cell responses. Recent transcriptomic analyses show hypoxia increases suppressive and anti-tumour immune signalling and infiltrates in bladder cancer. This study sought to investigate the relationship between hypoxia-inducible factor (HIF)-1 and -2, hypoxia, and immune signalling and infiltrates in MIBC. ChIP-seq was performed to identify HIF1α, HIF2α, and HIF1ß binding in the genome of the MIBC cell line T24 cultured in 1% and 0.1% oxygen for 24 h. Microarray data from four MIBC cell lines (T24, J82, UMUC3, and HT1376) cultured under 1%, 0.2%, and 0.1% oxygen for 24 h were used. Differences in the immune contexture between high- and low-hypoxia tumours were investigated using in silico analyses of two bladder cancer cohorts (BCON and TCGA) filtered to only include MIBC cases. GO and GSEA were used with the R packages "limma" and "fgsea". Immune deconvolution was performed using ImSig and TIMER algorithms. RStudio was used for all analyses. Under hypoxia, HIF1α and HIF2α bound to ~11.5-13.5% and ~4.5-7.5% of immune-related genes, respectively (1-0.1% O2). HIF1α and HIF2α both bound to genes associated with T cell activation and differentiation signalling pathways. HIF1α and HIF2α had distinct roles in immune-related signalling. HIF1 was associated with interferon production specifically, whilst HIF2 was associated with generic cytokine signalling as well as humoral and toll-like receptor immune responses. Neutrophil and myeloid cell signalling was enriched under hypoxia, alongside hallmark pathways associated with Tregs and macrophages. High-hypoxia MIBC tumours had increased expression of both suppressive and anti-tumour immune gene signatures and were associated with increased immune infiltrates. Overall, hypoxia is associated with increased inflammation for both suppressive and anti-tumour-related immune signalling and immune infiltrates, as seen in vitro and in situ using MIBC patient tumours.

Long non-coding RNA HIF1A-As2 and MYC form a double-positive feedback loop to promote cell proliferation and metastasis in KRAS-driven non-small cell lung cancer.

Lung cancer is the leading cause of cancer-related deaths worldwide. KRAS is the main oncogenic driver in lung cancer that can be activated by gene mutation or amplification, but whether long non-coding RNAs (lncRNAs) regulate its activation remains unknown. Through gain and loss of function approaches, we identified that lncRNA HIF1A-As2, a KRAS-induced lncRNA, is required for cell proliferation, epithelial-mesenchymal transition (EMT) and tumor propagation in non-small cell lung cancer (NSCLC) in vitro and in vivo. Integrative analysis of HIF1A-As2 transcriptomic profiling reveals that HIF1A-As2 modulates gene expression in trans, particularly regulating transcriptional factor genes including MYC. Mechanistically, HIF1A-As2 epigenetically activates MYC by recruiting DHX9 on MYC promoter, consequently stimulating the transcription of MYC and its target genes. In addition, KRAS promotes HIF1A-As2 expression via the induction of MYC, suggesting HIF1A-As2 and MYC form a double-regulatory loop to strengthen cell proliferation and tumor metastasis in lung cancer. Inhibition of HIF1A-As2 by LNA GapmeR antisense oligonucleotides (ASO) significantly improves sensitization to 10058-F4 (a MYC-specific inhibitor) and cisplatin treatment in PDX and KRASLSLG12D-driven lung tumors, respectively.

Tumour infiltrating B cells discriminate checkpoint blockade-induced responses.

BACKGROUND: Immune cell-driven anti-cancer activity is paramount for effective responses to checkpoint inhibitors (ICB). However, the contribution of the different immune cell subsets in the circulation and within the tumour is poorly understood. MATERIALS AND METHODS: To elucidate the role of the different cell subsets in anti-tumour responses elicited by ICB, we performed single-cell analysis of the transcriptome and surface proteome of paired pre- and early on-treatment metastatic melanoma tumour biopsies and matched peripheral blood mononuclear cell samples. We next compared the survival of metastatic melanoma patients treated with ICB according to the abundance of pre-treatment tumour-infiltrating B cell clonotypes. RESULTS: We identified cell clusters associated with disease control or progression, defined differential expression of biological pathways likely involved in the immune awakening against the tumour and examined how cell-cell communication patterns between the tumour cell subsets change during treatment. Furthermore, we discovered that B cells (immunoglobulin expression and abundance of B cell clonotypes) discriminate the clinical response after ICB and propose that B cells likely contribute to anti-tumour immunity by antigen presentation through major histocompatibility complex molecules. Finally, we demonstrated that the abundance of tumour-infiltrating B cell clonotypes at baseline identifies two distinct risk groups, a finding that we confirmed in an independent cohort. CONCLUSIONS: Our exploratory translational study provides new insights on the mechanistic role of B cells in anti-melanoma immunity during treatment with ICB. Additionally, we support pre-treatment B cell tumour infiltration as a promising prognostic biomarker to be further validated as a tool for clinical risk stratification.

AMPKα loss promotes KRAS-mediated lung tumorigenesis.

AMP-activated protein kinase (AMPK) is a critical sensor of energy status that coordinates cell growth with energy balance. In non-small cell lung cancer (NSCLC) the role of AMPKα is controversial and its contribution to lung carcinogenesis is not well-defined. Furthermore, it remains largely unknown whether long non-coding RNAs (lncRNAs) are involved in the regulation of AMPK-mediated pathways. Here, we found that loss of AMPKα in combination with activation of mutant KRASG12D increased lung tumour burden and reduced survival in KrasLSLG12D/+/AMPKαfl/fl mice. In agreement, functional in vitro studies revealed that AMPKα silencing increased growth and migration of NSCLC cells. In addition, we identified an AMPKα-modulated lncRNA, KIMAT1 (ENSG00000228709), which in turn regulates AMPKα activation by stabilizing the lactate dehydrogenase B (LDHB). Collectively, our study indicates that AMPKα loss promotes KRAS-mediated lung tumorigenesis and proposes a novel KRAS/KIMAT1/LDHB/AMPKα axis that could be exploited for therapeutic purposes.

A KRAS-responsive long non-coding RNA controls microRNA processing.

Wild-type KRAS (KRASWT) amplification has been shown to be a secondary means of KRAS activation in cancer and associated with poor survival. Nevertheless, the precise role of KRASWT overexpression in lung cancer progression is largely unexplored. Here, we identify and characterize a KRAS-responsive lncRNA, KIMAT1 (ENSG00000228709) and show that it correlates with KRAS levels both in cell lines and in lung cancer specimens. Mechanistically, KIMAT1 is a MYC target and drives lung tumorigenesis by promoting the processing of oncogenic microRNAs (miRNAs) through DHX9 and NPM1 stabilization while halting the biogenesis of miRNAs with tumor suppressor function via MYC-dependent silencing of p21, a component of the Microprocessor Complex. KIMAT1 knockdown suppresses not only KRAS expression but also KRAS downstream signaling, thereby arresting lung cancer growth in vitro and in vivo. Taken together, this study uncovers a role for KIMAT1 in maintaining a positive feedback loop that sustains KRAS signaling during lung cancer progression and provides a proof of principle that interfering with KIMAT1 could be a strategy to hamper KRAS-induced tumorigenesis.

CKAP2L Promotes Non-Small Cell Lung Cancer Progression through Regulation of Transcription Elongation.

Chromosomal instability (CIN) is a driver of clonal diversification and intratumor heterogeneity, providing genetic diversity that contributes to tumor progression. It is estimated that approximately 80% of solid cancers, including non-small cell lung cancer (NSCLC), exhibit features of CIN, which affects tumor growth and response to therapy. However, the molecular mechanisms connecting CIN to tumor progression are still poorly understood. Through an RNAi screen performed on genes involved in CIN and overexpressed in human lung adenocarcinoma samples, we identified the cytoskeleton-associated protein 2-like (CKAP2L) as a potential oncogene that promotes lung cancer proliferation and growth in vitro and in vivo. Mechanistically, CKAP2L directly interacted with RNA Pol II and regulated transcription elongation of key genes involved in spindle assembly checkpoint, chromosome segregation, cell cycle, and E2F signaling. Furthermore, depletion of CKAP2L increased the sensitivity of NSCLC cells to alvocidib, a pan-CDK inhibitor, leading to a significant reduction of cell proliferation and an increase in cell death. Altogether, these findings shed light on the molecular mechanisms through which CKAP2L, a protein involved in CIN, promotes cancer progression and suggest that its inhibition represents a novel therapeutic strategy in NSCLC. SIGNIFICANCE: These findings demonstrate the oncogenic function of CKAP2L through regulation of transcription elongation and suggest that targeting CKAP2L could enhance therapeutic response in patients with NSCLC.

Antagonistic Inflammatory Phenotypes Dictate Tumor Fate and Response to Immune Checkpoint Blockade.

Inflammation can support or restrain cancer progression and the response to therapy. Here, we searched for primary regulators of cancer-inhibitory inflammation through deep profiling of inflammatory tumor microenvironments (TMEs) linked to immune-dependent control in mice. We found that early intratumoral accumulation of interferon gamma (IFN-γ)-producing natural killer (NK) cells induced a profound remodeling of the TME and unleashed cytotoxic T cell (CTL)-mediated tumor eradication. Mechanistically, tumor-derived prostaglandin E2 (PGE2) acted selectively on EP2 and EP4 receptors on NK cells, hampered the TME switch, and enabled immune evasion. Analysis of patient datasets across human cancers revealed distinct inflammatory TME phenotypes resembling those associated with cancer immune control versus escape in mice. This allowed us to generate a gene-expression signature that integrated opposing inflammatory factors and predicted patient survival and response to immune checkpoint blockade. Our findings identify features of the tumor inflammatory milieu associated with immune control of cancer and establish a strategy to predict immunotherapy outcomes.

Vulnerability of drug-resistant EML4-ALK rearranged lung cancer to transcriptional inhibition.

A subset of lung adenocarcinomas is driven by the EML4-ALK translocation. Even though ALK inhibitors in the clinic lead to excellent initial responses, acquired resistance to these inhibitors due to on-target mutations or parallel pathway alterations is a major clinical challenge. Exploring these mechanisms of resistance, we found that EML4-ALK cells parental or resistant to crizotinib, ceritinib or alectinib are remarkably sensitive to inhibition of CDK7/12 with THZ1 and CDK9 with alvocidib or dinaciclib. These compounds robustly induce apoptosis through transcriptional inhibition and downregulation of anti-apoptotic genes. Importantly, alvocidib reduced tumour progression in xenograft mouse models. In summary, our study takes advantage of the transcriptional addiction hypothesis to propose a new treatment strategy for a subset of patients with acquired resistance to first-, second- and third-generation ALK inhibitors.

Profiling of Circulating Free DNA Using Targeted and Genome-wide Sequencing in Patients with SCLC.

INTRODUCTION: SCLC accounts for approximately 250,000 deaths worldwide each year. Acquisition of adequate tumor biopsy samples is challenging, and liquid biopsies present an alternative option for patient stratification and response monitoring. METHODS: We applied whole genome next-generation sequencing to circulating free DNA (cfDNA) from 39 patients with limited-stage (LS) SCLC and 30 patients with extensive-stage SCLC to establish genome-wide copy number aberrations and also performed targeted mutation analysis of 110 SCLC associated genes. Quantitative metrics were calculated for copy number aberrations, including percent genome amplified (PGA [the percentage of genomic regions amplified]), Z-score (a measure of standard deviation), and Moran's I (a measure of spatial autocorrelation). In addition CellSearch, an epitope-dependent enrichment platform, was used to enumerate circulating tumor cells (CTCs) from a parallel blood sample. RESULTS: Genome-wide and targeted cfDNA sequencing data identified tumor-related changes in 94% of patients with LS SCLC and 100% of patients with extensive-stage SCLC. Parallel analysis of CTCs based on at least 1 CTC/7.5 mL of blood increased tumor detection frequencies to 95% for LS SCLC. Both CTC counts and cfDNA readouts correlated with disease stage and overall survival. CONCLUSIONS: We demonstrate that a simple cfDNA genome-wide copy number approach provides an effective means of monitoring patients through treatment and show that targeted cfDNA sequencing identifies potential therapeutic targets in more than 50% of patients. We are now incorporating this approach into additional studies and trials of targeted therapies.

Analysis of circulating cell-free DNA identifies KRAS copy number gain and mutation as a novel prognostic marker in Pancreatic cancer.

Serial biopsy of pancreatic ductal adenocarcinoma (PDAC), to chart tumour evolution presents a significant challenge. We examined the utility of circulating free DNA (cfDNA) as a minimally invasive approach across a cohort of 55 treatment-naïve patients with PDAC; 31 with metastatic and 24 with locally advanced disease. Somatic mutations in cfDNA were detected using next generation sequencing in 15/24 (62.5%) and 27/31 (87%) of patients with locally advanced and metastatic disease, respectively. Copy number changes were detected in cfDNA of 10 patients of whom 7 exhibited gain of chromosome 12p harbouring KRAS as well as a canonical KRAS codon 12 mutation. In multivariable Cox Regression analysis, we show for the first time that patients with KRAS copy number gain and KRAS mutation have significantly worse outcomes, suggesting that this may be linked to PDAC progression. The simple cfDNA assay we describe will enable determination of the presence of KRAS copy number gain and KRAS mutations in larger studies and clinical trials.

Signaling pathway screening platforms are an efficient approach to identify therapeutic targets in cancers that lack known driver mutations: a case report for a cancer of unknown primary origin.

Precision medicine aims to tailor cancer therapies to target specific tumor-promoting aberrations. For tumors that lack actionable drivers, which occurs frequently in the clinic, extensive molecular characterization and pre-clinical drug efficacy studies will be required. A cell line maintained at low passage and a patient- derived xenograft model (PDX) were generated using a fresh biopsy from a patient with a poorly-differentiated neuroendocrine tumor of unknown primary origin. Next-generation sequencing, high throughput signaling network analysis, and drug efficacy trials were then conducted to identify actionable targets for therapeutic intervention. No actionable mutations were identified after whole exome sequencing of the patient's DNA. However, whole genome sequencing revealed amplification of the 3q and 5p chromosomal arms, that include the PIK3CA and RICTOR genes, respectively. We then conducted pathway analysis, which revealed activation of the AKT pathway. Based on this analysis, efficacy of PIK3CA and AKT inhibitors were evaluated in the tumor biopsy-derived cell culture and PDX, and response to the AKT inhibitor AZD5363 was observed both in vitro and in vivo indicating the patient would benefit from targeted therapies directed against the serine/threonine kinase AKT. In conclusion, our study demonstrates that high throughput signaling pathway analysis will significantly aid in identifying actionable alterations in rare tumors and guide patient stratification into early-phase clinical trials.

Enhancer Activation by Pharmacologic Displacement of LSD1 from GFI1 Induces Differentiation in Acute Myeloid Leukemia.

Pharmacologic inhibition of LSD1 promotes blast cell differentiation in acute myeloid leukemia (AML) with MLL translocations. The assumption has been that differentiation is induced through blockade of LSD1's histone demethylase activity. However, we observed that rapid, extensive, drug-induced changes in transcription occurred without genome-wide accumulation of the histone modifications targeted for demethylation by LSD1 at sites of LSD1 binding and that a demethylase-defective mutant rescued LSD1 knockdown AML cells as efficiently as wild-type protein. Rather, LSD1 inhibitors disrupt the interaction of LSD1 and RCOR1 with the SNAG-domain transcription repressor GFI1, which is bound to a discrete set of enhancers located close to transcription factor genes that regulate myeloid differentiation. Physical separation of LSD1/RCOR1 from GFI1 is required for drug-induced differentiation. The consequent inactivation of GFI1 leads to increased enhancer histone acetylation within hours, which directly correlates with the upregulation of nearby subordinate genes.

KRAS induces lung tumorigenesis through microRNAs modulation.

Oncogenic KRAS induces tumor onset and development by modulating gene expression via different molecular mechanisms. MicroRNAs (miRNAs) are small non-coding RNAs that have been established as main players in tumorigenesis. By overexpressing wild type or mutant KRAS (KRASG12D) and using inducible human and mouse cell lines, we analyzed KRAS-regulated microRNAs in non-small-cell lung cancer (NSCLC). We show that miR-30c and miR-21 are significantly upregulated by both KRAS isoforms and induce drug resistance and enhance cell migration/invasion via inhibiting crucial tumor suppressor genes, such as NF1, RASA1, BID, and RASSF8. MiR-30c and miR-21 levels were significantly elevated in tumors from patients that underwent surgical resection of early stages NSCLC compared to normal lung and in plasma from the same patients. Systemic delivery of LNA-anti-miR-21 in combination with cisplatin in vivo completely suppressed the development of lung tumors in a mouse model of lung cancer. Mechanistically, we demonstrated that ELK1 is responsible for miR-30c and miR-21 transcriptional activation by direct binding to the miRNA proximal promoter regions. In summary, our study defines that miR-30c and miR-21 may be valid biomarkers for early NSCLC detection and their silencing could be beneficial for therapeutic applications.

PDGFR-modulated miR-23b cluster and miR-125a-5p suppress lung tumorigenesis by targeting multiple components of KRAS and NF-kB pathways.

In NSCLC alterations in PDGF receptors are markers of worst prognosis and efficient targeting of these receptors is yet to be achieved. In this study, we explored PDGFR-regulated microRNAs demonstrating that miR-23b cluster and miR-125a-5p are downregulated by increased expression of PDGFR-α or PDGFR-ß in NSCLC cells. Mechanistically, the expression of these microRNAs is positively regulated by p53 and negatively modulated by NF-kB p65. Forced expression of miR-23b cluster or miR-125a-5p enhanced drug sensitivity and suppressed invasiveness of NSCLC cells by silencing several genes involved in oncogenic KRAS and NF-kB pathways, including SOS1, GRB2, IQGAP1, RALA, RAF-1, IKKß, AKT2, ERK2 and KRAS itself. Of note, an inverse correlation between miR-23b cluster, miR-125a-5p and respective target genes was also found in vivo in a large dataset of lung adenocarcinoma samples. Furthermore, in vivo delivery of miR-23b cluster or miR-125a-5p significantly repressed tumour growth in a highly aggressive NSCLC circulating tumour cell (CTC) patient derived explant (CDX) mouse model. In conclusion, our finding sheds light on the PDGFR signaling and endorses the possibility to employ miR-23b cluster and miR-125a-5p as therapeutic tools to silence simultaneously a range of redundant pathways and main effectors of tumorigenesis in NSCLC.

ProbFold: a probabilistic method for integration of probing data in RNA secondary structure prediction.

MOTIVATION: Recently, new RNA secondary structure probing techniques have been developed, including Next Generation Sequencing based methods capable of probing transcriptome-wide. These techniques hold great promise for improving structure prediction accuracy. However, each new data type comes with its own signal properties and biases, which may even be experiment specific. There is therefore a growing need for RNA structure prediction methods that can be automatically trained on new data types and readily extended to integrate and fully exploit multiple types of data. RESULTS: Here, we develop and explore a modular probabilistic approach for integrating probing data in RNA structure prediction. It can be automatically trained given a set of known structures with probing data. The approach is demonstrated on SHAPE datasets, where we evaluate and selectively model specific correlations. The approach often makes superior use of the probing data signal compared to other methods. We illustrate the use of ProbFold on multiple data types using both simulations and a small set of structures with both SHAPE, DMS and CMCT data. Technically, the approach combines stochastic context-free grammars (SCFGs) with probabilistic graphical models. This approach allows rapid adaptation and integration of new probing data types. AVAILABILITY AND IMPLEMENTATION: ProbFold is implemented in C ++. Models are specified using simple textual formats. Data reformatting is done using separate C ++ programs. Source code, statically compiled binaries for x86 Linux machines, C ++ programs, example datasets and a tutorial is available from http://moma.ki.au.dk/prj/probfold/ CONTACT: : jakob.skou@clin.au.dk SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Derivative scores from site accessibility and ranking of miRNA target predictions.

In the present study, we define derivative scoring functions from PITA and STarMir predictions. The scoring functions are evaluated for up to five selected miRNAs with a relatively large number of validated targets reported by TarBase and miRecords. The average ranking of validated targets returned by PITA and STarMir is compared to the average ranking produced by the new derivatives scores. We obtain an average improvement of 13.6% (STD∼5.7%) relative to the average ranking of validated targets produced by PITA and STarMir.

Ranking of microRNA target prediction scores by Pareto front analysis.

Over the past ten years, a variety of microRNA target prediction methods has been developed, and many of the methods are constantly improved and adapted to recent insights into miRNA-mRNA interactions. In a typical scenario, different methods return different rankings of putative targets, even if the ranking is reduced to selected mRNAs that are related to a specific disease or cell type. For the experimental validation it is then difficult to decide in which order to process the predicted miRNA-mRNA bindings, since each validation is a laborious task and therefore only a limited number of mRNAs can be analysed. We propose a new ranking scheme that combines ranked predictions from several methods and - unlike standard thresholding methods - utilises the concept of Pareto fronts as defined in multi-objective optimisation. In the present study, we attempt a proof of concept by applying the new ranking scheme to hsa-miR-21, hsa-miR-125b, and hsa-miR-373 and prediction scores supplied by PITA and RNAhybrid. The scores are interpreted as a two-objective optimisation problem, and the elements of the Pareto front are ranked by the STarMir score with a subsequent re-calculation of the Pareto front after removal of the top-ranked mRNA from the basic set of prediction scores. The method is evaluated on validated targets of the three miRNA, and the ranking is compared to scores from DIANA-microT and TargetScan. We observed that the new ranking method performs well and consistent, and the first validated targets are elements of Pareto fronts at a relatively early stage of the recurrent procedure, which encourages further research towards a higher-dimensional analysis of Pareto fronts.