Abstract
Immune checkpoint inhibitors targeting the PD-1/PD-L1 axis lead to durable clinical responses in subsets of cancer patients across multiple indications, including non-small cell lung cancer ...(NSCLC), urothelial carcinoma (UC) and renal cell carcinoma (RCC). Herein, we complement PD-L1 immunohistochemistry (IHC) and tumor mutation burden (TMB) with RNA-seq in 366 patients to identify unifying and indication-specific molecular profiles that can predict response to checkpoint blockade across these tumor types. Multiple machine learning approaches failed to identify a baseline transcriptional signature highly predictive of response across these indications. Signatures described previously for immune checkpoint inhibitors also failed to validate. At the pathway level, significant heterogeneity is observed between indications, in particular within the PD-L1
+
tumors. mUC and NSCLC are molecularly aligned, with cell cycle and DNA damage repair genes associated with response in PD-L1- tumors. At the gene level, the CDK4/6 inhibitor
CDKN2A
is identified as a significant transcriptional correlate of response, highlighting the association of non-immune pathways to the outcome of checkpoint blockade. This cross-indication analysis reveals molecular heterogeneity between mUC, NSCLC and RCC tumors, suggesting that indication-specific molecular approaches should be prioritized to formulate treatment strategies.
Integrated multi-omics evaluation of 823 tumors from advanced renal cell carcinoma (RCC) patients identifies molecular subsets associated with differential clinical outcomes to angiogenesis blockade ...alone or with a checkpoint inhibitor. Unsupervised transcriptomic analysis reveals seven molecular subsets with distinct angiogenesis, immune, cell-cycle, metabolism, and stromal programs. While sunitinib and atezolizumab + bevacizumab are effective in subsets with high angiogenesis, atezolizumab + bevacizumab improves clinical benefit in tumors with high T-effector and/or cell-cycle transcription. Somatic mutations in PBRM1 and KDM5C associate with high angiogenesis and AMPK/fatty acid oxidation gene expression, while CDKN2A/B and TP53 alterations associate with increased cell-cycle and anabolic metabolism. Sarcomatoid tumors exhibit lower prevalence of PBRM1 mutations and angiogenesis markers, frequent CDKN2A/B alterations, and increased PD-L1 expression. These findings can be applied to molecularly stratify patients, explain improved outcomes of sarcomatoid tumors to checkpoint blockade versus antiangiogenics alone, and develop personalized therapies in RCC and other indications.
Display omitted
•Genomics of 823 RCC tumors, including 134 sarcomatoid tumors, reveals 7 subtypes•Subtype specific angiogenesis, immune, metabolic, stromal, and cell-cycle profiles•Differential prevalence of PBRM1, KDM5C, CDKN2A/2B, and TP53 alterations in subsets•Differential outcomes to VEGF blockade alone or in combination with anti-PD-L1
Motzer et al. perform integrative multi-omics analyses of 823 renal cancer tumors from a randomized clinical trial. A robust molecular classification scheme, based on transcriptional and gene alteration profiles and differential clinical outcomes to VEGF blockade alone or in combination with anti-PD-L1, informs personalized treatment strategies and future therapeutic development in RCC.
The normalization of RNA-seq data is essential for accurate downstream inference, but the assumptions upon which most normalization methods are based are not applicable in the single-cell setting. ...Consequently, applying existing normalization methods to single-cell RNA-seq data introduces artifacts that bias downstream analyses. To address this, we introduce SCnorm for accurate and efficient normalization of single-cell RNA-seq data.
Although elevated plasma interleukin-8 (pIL-8) has been associated with poor outcome to immune checkpoint blockade
, this has not been comprehensively evaluated in large randomized studies. Here we ...analyzed circulating pIL-8 and IL8 gene expression in peripheral blood mononuclear cells and tumors of patients treated with atezolizumab (anti-PD-L1 monoclonal antibody) from multiple randomized trials representing 1,445 patients with metastatic urothelial carcinoma (mUC) and metastatic renal cell carcinoma. High levels of IL-8 in plasma, peripheral blood mononuclear cells and tumors were associated with decreased efficacy of atezolizumab in patients with mUC and metastatic renal cell carcinoma, even in tumors that were classically CD8
T cell inflamed. Low baseline pIL-8 in patients with mUC was associated with increased response to atezolizumab and chemotherapy. Patients with mUC who experienced on-treatment decreases in pIL-8 exhibited improved overall survival when treated with atezolizumab but not with chemotherapy. Single-cell RNA sequencing of the immune compartment showed that IL8 is primarily expressed in circulating and intratumoral myeloid cells and that high IL8 expression is associated with downregulation of the antigen-presentation machinery. Therapies that can reverse the impacts of IL-8-mediated myeloid inflammation will be essential for improving outcomes of patients treated with immune checkpoint inhibitors.
Messenger RNA expression is important in normal development and differentiation, as well as in manifestation of disease. RNA-seq experiments allow for the identification of differentially expressed ...(DE) genes and their corresponding isoforms on a genome-wide scale. However, statistical methods are required to ensure that accurate identifications are made. A number of methods exist for identifying DE genes, but far fewer are available for identifying DE isoforms. When isoform DE is of interest, investigators often apply gene-level (count-based) methods directly to estimates of isoform counts. Doing so is not recommended. In short, estimating isoform expression is relatively straightforward for some groups of isoforms, but more challenging for others. This results in estimation uncertainty that varies across isoform groups. Count-based methods were not designed to accommodate this varying uncertainty, and consequently, application of them for isoform inference results in reduced power for some classes of isoforms and increased false discoveries for others.
Taking advantage of the merits of empirical Bayesian methods, we have developed EBSeq for identifying DE isoforms in an RNA-seq experiment comparing two or more biological conditions. Results demonstrate substantially improved power and performance of EBSeq for identifying DE isoforms. EBSeq also proves to be a robust approach for identifying DE genes.
An R package containing examples and sample datasets is available at http://www.biostat.wisc.edu/kendzior/EBSEQ/.
Supplementary data are available at Bioinformatics online.
Here, we report the derivation of arterial endothelial cells from human pluripotent stem cells that exhibit arterial-specific functions in vitro and in vivo. We combine single-cell RNA sequencing of ...embryonic mouse endothelial cells with an EFNB2-tdTomato/EPHB4-EGFP dual reporter human embryonic stem cell line to identify factors that regulate arterial endothelial cell specification. The resulting xeno-free protocol produces cells with gene expression profiles, oxygen consumption rates, nitric oxide production levels, shear stress responses, and TNFα-induced leukocyte adhesion rates characteristic of arterial endothelial cells. Arterial endothelial cells were robustly generated from multiple human embryonic and induced pluripotent stem cell lines and have potential applications for both disease modeling and regenerative medicine.
High-throughput expression profiling experiments with ordered conditions (e.g. time-course or spatial-course) are becoming more common for studying detailed differentiation processes or spatial ...patterns. Identifying dynamic changes at both the individual gene and whole transcriptome level can provide important insights about genes, pathways, and critical time points.
We present an R package, Trendy, which utilizes segmented regression models to simultaneously characterize each gene's expression pattern and summarize overall dynamic activity in ordered condition experiments. For each gene, Trendy finds the optimal segmented regression model and provides the location and direction of dynamic changes in expression. We demonstrate the utility of Trendy to provide biologically relevant results on both microarray and RNA-sequencing (RNA-seq) datasets.
Trendy is a flexible R package which characterizes gene-specific expression patterns and summarizes changes of global dynamics over ordered conditions. Trendy is freely available on Bioconductor with a full vignette at https://bioconductor.org/packages/release/bioc/html/Trendy.html .
Papillary renal cell carcinoma (PRCC) accounts for about 10 percent of all renal cell carcinomas, and the prognosis is poor for people with advanced disease. Interleukin-20 receptor subunit beta ...(IL20RB) is a single-pass type I membrane protein of the type II cytokine receptor family and is related to the pathogenesis of chronic inflammation and autoimmune diseases, including psoriasis, glaucoma, vitiligo, rheumatoid arthritis, and inflammatory bowel disease. However, little has been reported on IL20RB with respect to cancer, especially in PRCC. Thus, we performed this study to explore its biological characteristics in PRCC. Data from the TCGA database were used to analyze the expression and prognosis of IL20RB. qRT-PCR was used to detect the expression of IL20RB in PRCC cells in vitro. After knockdown of IL20RB with small interfering RNA (siRNA) technology, the proliferation, migration, and invasion of Ketr-3 cells and the expression of related proteins in the epithelial-mesenchymal transition (EMT) pathway were measured with Cell Counting Kit-8 (CCK-8), transwell, and western blot assays. The findings demonstrated that the expression of IL20RB was upregulated in both PRCC tissues and cells and that the high expression of IL20RB led to low overall survival (OS). Furthermore, after knockdown of IL20RB in vitro, the proliferation, migration, and invasion of Ketr-3 cells were reduced, and the expression of related proteins in the EMT pathway declined, suggesting that IL20RB plays a vital role in PRCC through the EMT pathway. These results reveal the biological significance of IL20RB in PRCC and provide new insight for future targeted drugs.
The axolotl (Ambystoma mexicanum) has long been the subject of biological research, primarily owing to its outstanding regenerative capabilities. However, the gene expression programs governing its ...embryonic development are particularly underexplored, especially when compared to other amphibian model species. Therefore, we performed whole transcriptome polyA+ RNA sequencing experiments on 17 stages of embryonic development. As the axolotl genome is unsequenced and its gene annotation is incomplete, we built de novo transcriptome assemblies for each stage and garnered functional annotation by comparing expressed contigs with known genes in other organisms. In evaluating the number of differentially expressed genes over time, we identify three waves of substantial transcriptome upheaval each followed by a period of relative transcriptome stability. The first wave of upheaval is between the one and two cell stage. We show that the number of differentially expressed genes per unit time is higher between the one and two cell stage than it is across the mid-blastula transition (MBT), the period of zygotic genome activation. We use total RNA sequencing to demonstrate that the vast majority of genes with increasing polyA+ signal between the one and two cell stage result from polyadenylation rather than de novo transcription. The first stable phase begins after the two cell stage and continues until the mid-blastula transition, corresponding with the pre-MBT phase of transcriptional quiescence in amphibian development. Following this is a peak of differential gene expression corresponding with the activation of the zygotic genome and a phase of transcriptomic stability from stages 9–11. We observe a third wave of transcriptomic change between stages 11 and 14, followed by a final stable period. The last two stable phases have not been documented in amphibians previously and correspond to times of major morphogenic change in the axolotl embryo: gastrulation and neurulation. These results yield new insights into global gene expression during early stages of amphibian embryogenesis and will help to further develop the axolotl as a model species for developmental and regenerative biology.
•Sequenced transcriptomes of 17 stages of embryonic development in the axolotl.•Peaks of transcriptional dynamism punctuate relatively stable states.•Stable states occur during major morphogenetic events: gastrulation, neurulation.
The salamander has the remarkable ability to regenerate its limb after amputation. Cells at the site of amputation form a blastema and then proliferate and differentiate to regrow the limb. To better ...understand this process, we performed deep RNA sequencing of the blastema over a time course in the axolotl, a species whose genome has not been sequenced. Using a novel comparative approach to analyzing RNA-seq data, we characterized the transcriptional dynamics of the regenerating axolotl limb with respect to the human gene set. This approach involved de novo assembly of axolotl transcripts, RNA-seq transcript quantification without a reference genome, and transformation of abundances from axolotl contigs to human genes. We found a prominent burst in oncogene expression during the first day and blastemal/limb bud genes peaking at 7 to 14 days. In addition, we found that limb patterning genes, SALL genes, and genes involved in angiogenesis, wound healing, defense/immunity, and bone development are enriched during blastema formation and development. Finally, we identified a category of genes with no prior literature support for limb regeneration that are candidates for further evaluation based on their expression pattern during the regenerative process.