Background Non-alcoholic fatty liver disease (NAFLD) encompasses a wide spectrum of liver pathologies. However, no medical treatment has been approved for the treatment of NAFLD. In our previous study, we found that PKLR could be a potential target for treatment of NALFD. Here, we investigated the effect of PKLR in in vivo model and performed drug repositioning to identify a drug candidate for treatment of NAFLD. Methods Tissue samples from liver, muscle, white adipose and heart were obtained from control and PKLR knock-out mice fed with chow and high sucrose diets. Lipidomics as well as transcriptomics analyses were conducted using these tissue samples. In addition, a computational drug repositioning analysis was performed and drug candidates were identified. The drug candidates were both tested in in vitro and in vivo models to evaluate their toxicity and efficacy. Findings The Pklr KO reversed the increased hepatic triglyceride level in mice fed with high sucrose diet and partly recovered the transcriptomic changes in the liver as well as in other three tissues. Both liver and white adipose tissues exhibited dysregulated circadian transcriptomic profiles, and these dysregulations were reversed by hepatic knockout of Pklr. In addition, 10 small molecule drug candidates were identified as potential inhibitor of PKLR using our drug repositioning pipeline, and two of them significantly inhibited both the PKLR expression and triglyceride level in in vitro model. Finally, the two selected small molecule drugs were evaluated in in vivo rat models and we found that these drugs attenuate the hepatic steatosis without side effect on other tissues. Interpretation In conclusion, our study provided biological insights about the critical role of PKLR in NAFLD progression and proposed a treatment strategy for NAFLD patients, which has been validated in preclinical studies. 

Simple Summary Hepatocellular carcinoma (HCC) is the most common malignancy of liver cancer. However, treatment of HCC is still severely limited due to limitation of drug therapy. We aimed to screen more possible target genes and candidate drugs for HCC, exploring the possibility of drug treatments from systems biological perspective. We identified ten candidate target genes, which are hub genes in HCC co-expression networks, which also possess significant prognostic value in two independent HCC cohorts. The rationality of these target genes was well demonstrated through variety analyses of patient expression profiles. We then screened candidate drugs for target genes and finally identified withaferin-a and mitoxantrone as the candidate drug for HCC treatment. The drug effectiveness was validated in in vitro model and computational analysis, providing more evidence for our drug repositioning method and results. Hepatocellular carcinoma (HCC) is a malignant liver cancer that continues to increase deaths worldwide owing to limited therapies and treatments. Computational drug repurposing is a promising strategy to discover potential indications of existing drugs. In this study, we present a systematic drug repositioning method based on comprehensive integration of molecular signatures in liver cancer tissue and cell lines. First, we identify robust prognostic genes and two gene co-expression modules enriched in unfavorable prognostic genes based on two independent HCC cohorts, which showed great consistency in functional and network topology. Then, we screen 10 genes as potential target genes for HCC on the bias of network topology analysis in these two modules. Further, we perform a drug repositioning method by integrating the shRNA and drug perturbation of liver cancer cell lines and identifying potential drugs for every target gene. Finally, we evaluate the effects of the candidate drugs through an in vitro model and observe that two identified drugs inhibited the protein levels of their corresponding target genes and cell migration, also showing great binding affinity in protein docking analysis. Our study demonstrates the usefulness and efficiency of network-based drug repositioning approach to discover potential drugs for cancer treatment and precision medicine approach.

Introduction: Macrophages and T cells play crucial roles in liver physiology, but their functional diversity in hepatocellular carcinoma (HCC) remains largely unknown. Methods: Two bulk RNA-sequencing (RNA-seq) cohorts for HCC were analyzed using gene co-expression network analysis. Key gene modules and networks were mapped to single-cell RNA-sequencing (scRNA-seq) data of HCC. Cell type fraction of bulk RNA-seq data was estimated by deconvolution approach using single-cell RNA-sequencing data as a reference. Survival analysis was carried out to estimate the prognosis of different immune cell types in bulk RNA-seq cohorts. Cell-cell interaction analysis was performed to identify potential links between immune cell types in HCC. Results: In this study, we analyzed RNA-seq data from two large-scale HCC cohorts, revealing a major and consensus gene co-expression cluster with significant implications for immunosuppression. Notably, these genes exhibited higher enrichment in liver macrophages than T cells, as confirmed by scRNA-seq data from HCC patients. Integrative analysis of bulk and single-cell RNA-seq data pinpointed SPP1+ macrophages as an unfavorable cell type, while VCAN+ macrophages, C1QA+ macrophages, and CD8+ T cells were associated with a more favorable prognosis for HCC patients. Subsequent scRNA-seq investigations and in vitro experiments elucidated that SPP1, predominantly secreted by SPP1+ macrophages, inhibits CD8+ T cell proliferation. Finally, targeting SPP1 in tumor-associated macrophages through inhibition led to a shift towards a favorable phenotype. Discussion: This study underpins the potential of SPP1 as a translational target in immunotherapy for HCC.

New technologies enable single-cell transcriptome analysis, mapping genome-wide expression across the human body. Here, we present an extended analysis of protein-coding genes in all major human tissues and organs, combining single-cell and bulk transcriptomics. To enhance transcriptome depth, 31 tissues were analyzed using a pooling method, identifying 557 unique cell clusters, manually annotated by marker gene expression. Genes were classified by body-wide expression and validated through antibody-based profiling. All results are available in the updated open-access Single Cell Type section of the Human Protein Atlas (www.proteinatlas.org) for genome-wide exploration of genes, proteins, and their spatial distribution in cells.