Unbiased kidney-centric molecular categorization of chronic kidney disease as a step towards precision medicine.

Current classification of chronic kidney disease (CKD) into stages using indirect systemic measures (estimated glomerular filtration rate (eGFR) and albuminuria) is agnostic to the heterogeneity of underlying molecular processes in the kidney thereby limiting precision medicine approaches. To generate a novel CKD categorization that directly reflects within kidney disease drivers we analyzed publicly available transcriptomic data from kidney biopsy tissue. A Self-Organizing Maps unsupervised artificial neural network machine-learning algorithm was used to stratify a total of 369 patients with CKD and 46 living kidney donors as healthy controls. Unbiased stratification of the discovery cohort resulted in identification of four novel molecular categories of disease termed CKD-Blue, CKD-Gold, CKD-Olive, CKD-Plum that were replicated in independent CKD and diabetic kidney disease datasets and can be further tested on any external data at kidneyclass.org. Each molecular category spanned across CKD stages and histopathological diagnoses and represented transcriptional activation of distinct biological pathways. Disease progression rates were highly significantly different between the molecular categories. CKD-Gold displayed rapid progression, with significant eGFR-adjusted Cox regression hazard ratio of 5.6 [1.01-31.3] for kidney failure and hazard ratio of 4.7 [1.3-16.5] for composite of kidney failure or a 40% or more eGFR decline. Urine proteomics revealed distinct patterns between the molecular categories, and a 25-protein signature was identified to distinguish CKD-Gold from other molecular categories. Thus, patient stratification based on kidney tissue omics offers a gateway to non-invasive biomarker-driven categorization and the potential for future clinical implementation, as a key step towards precision medicine in CKD.

Small Gene Networks Delineate Immune Cell States and Characterize Immunotherapy Response in Melanoma.

Single-cell technologies have elucidated mechanisms responsible for immune checkpoint inhibitor (ICI) response, but are not amenable to a clinical diagnostic setting. In contrast, bulk RNA sequencing (RNA-seq) is now routine for research and clinical applications. Our workflow uses transcription factor (TF)-directed coexpression networks (regulons) inferred from single-cell RNA-seq data to deconvolute immune functional states from bulk RNA-seq data. Regulons preserve the phenotypic variation in CD45+ immune cells from metastatic melanoma samples (n = 19, discovery dataset) treated with ICIs, despite reducing dimensionality by >100-fold. Four cell states, termed exhausted T cells, monocyte lineage cells, memory T cells, and B cells were associated with therapy response, and were characterized by differentially active and cell state-specific regulons. Clustering of bulk RNA-seq melanoma samples from four independent studies (n = 209, validation dataset) according to regulon-inferred scores identified four groups with significantly different response outcomes (P < 0.001). An intercellular link was established between exhausted T cells and monocyte lineage cells, whereby their cell numbers were correlated, and exhausted T cells predicted prognosis as a function of monocyte lineage cell number. The ligand-receptor expression analysis suggested that monocyte lineage cells drive exhausted T cells into terminal exhaustion through programs that regulate antigen presentation, chronic inflammation, and negative costimulation. Together, our results demonstrate how regulon-based characterization of cell states provide robust and functionally informative markers that can deconvolve bulk RNA-seq data to identify ICI responders.

Distinct mechanisms govern populations of myeloid-derived suppressor cells in chronic viral infection and cancer.

Myeloid-derived suppressor cells (MDSCs) are major negative regulators of immune responses in cancer and chronic infections. It remains unclear if regulation of MDSC activity in different conditions is controlled by similar mechanisms. We compared MDSCs in mice with cancer and lymphocytic choriomeningitis virus (LCMV) infection. Chronic LCMV infection caused the development of monocytic MDSCs (M-MDSCs) but did not induce polymorphonuclear MDSCs (PMN-MDSCs). In contrast, both MDSC populations were present in cancer models. An acquisition of immune-suppressive activity by PMN-MDSCs in cancer was controlled by IRE1α and ATF6 pathways of the endoplasmic reticulum (ER) stress response. Abrogation of PMN-MDSC activity by blockade of the ER stress response resulted in an increase in tumor-specific immune response and reduced tumor progression. In contrast, the ER stress response was dispensable for suppressive activity of M-MDSCs in cancer and LCMV infection. Acquisition of immune-suppressive activity by M-MDSCs in spleens was mediated by IFN-γ signaling. However, it was dispensable for suppressive activity of M-MDSCs in tumor tissues. Suppressive activity of M-MDSCs in tumors was retained due to the effect of IL-6 present at high concentrations in the tumor site. These results demonstrate disease- and population-specific mechanisms of MDSC accumulation and the need for targeting different pathways to achieve inactivation of these cells.

A Framework for Multi-Omic Prediction of Treatment Response to Biologic Therapy for Psoriasis.

Biologic therapies have shown high efficacy in psoriasis, but individual response varies and is poorly understood. To inform biomarker discovery in the Psoriasis Stratification to Optimise Relevant Therapy (i.e., PSORT) study, we evaluated a comprehensive array of omics platforms across three time points and multiple tissues in a pilot investigation of 10 patients with severe psoriasis, treated with the tumor necrosis factor (TNF) inhibitor, etanercept. We used RNA sequencing to analyze mRNA and small RNA transcriptome in blood, lesional and nonlesional skin, and the SOMAscan platform to investigate the serum proteome. Using an integrative systems biology approach, we identified signals of treatment response in genes and pathways associated with TNF signaling, psoriasis pathology, and the major histocompatibility complex region. We found association between clinical response and TNF-regulated genes in blood and skin. Using a combination of differential expression testing, upstream regulator analysis, clustering techniques, and predictive modeling, we show that baseline samples are indicative of patient response to biologic therapies, including signals in blood, which have traditionally been considered unreliable for inference in dermatology. In conclusion, our pilot study provides both an analytical framework and empirical basis to estimate power for larger studies, specifically the ongoing PSORT study, which we show as powered for biomarker discovery and patient stratification.

Haplosaurus computes protein haplotypes for use in precision drug design.

Selecting the most appropriate protein sequences is critical for precision drug design. Here we describe Haplosaurus, a bioinformatic tool for computation of protein haplotypes. Haplosaurus computes protein haplotypes from pre-existing chromosomally-phased genomic variation data. Integration into the Ensembl resource provides rapid and detailed protein haplotypes retrieval. Using Haplosaurus, we build a database of unique protein haplotypes from the 1000 Genomes dataset reflecting real-world protein sequence variability and their prevalence. For one in seven genes, their most common protein haplotype differs from the reference sequence and a similar number differs on their most common haplotype between human populations. Three case studies show how knowledge of the range of commonly encountered protein forms predicted in populations leads to insights into therapeutic efficacy. Haplosaurus and its associated database is expected to find broad applications in many disciplines using protein sequences and particularly impactful for therapeutics design.