Background: Increasing evdence supports the suggestion that the immune cell infiltration (ICI) patterns play a pivotal role in tumor progression in breast cancer (BRCA). Nonetheless, there has been no comprehensive analysis of the ICI patterns effects on the clinical outcomes and immunotherapy.Methods: Multiomic data for BRCA samples were downloaded from TCGA. ESTIMATE algorithm, ssGSEA method, and CIBERSORT analysis were used to uncover the landscape of the tumor immune microenvironment (TIME). BRCA subtypes based on the ICI pattern were identified by consensus clustering and principal-component analysis was performed to obtain the ICI scores to quantify the ICI patterns in individual tumors. Their prognostic value was validated by the Kaplan-Meier survival curves. Gene set enrichment analysis (GSEA) was applied for functional annotation. Immunophenoscore (IPS) was employed to explore the immunotherapeutic role of the ICI scores. Finally, the mutation data was analyzed by using the “maftools” R package.Results: Three different immune infiltration patterns with a distinct prognosis and biological signature were recognized among 1,198 BRCA samples. The characteristics of TIME under these three patterns were highly consistent with three known immune profiles: immune- excluded, immune-desert, and immune-inflamed phenotypes, respectively. The identification of the ICI patterns within individual tumors based on the ICI score, developed under the ICI-related signature genes, contributed into dissecting biological processes, clinical outcome, immune cells infiltration, immunotherapeutic effect, and genetic variation. High ICI score subtype, characterized with a suppression of immunity, suggested an immune-exhausted phenotype. Abundant effective immune cells were discovered in the low ICI score patients, which corresponded to an immune-activated phenotype and might present an immunotherapeutic advantage. Immunophenoscore was implemented as a surrogate of immunotherapeutic outcome, low-ICI scores samples obtained a significantly higher immunophenoscore. Enrichment of the JAK/STAT and VEGF signal pathways were activated in the ICI low-score subgroup. Finally, the synergistic effect between the ICI score and the tumor mutation burden (TMB) was confirmed.Conclusion: This work comprehensively elucidated that the ICI patterns served as an indispensable player in complexity and diversity of TIME. Quantitative identification of the ICI patterns in individual tumor will contribute into mapping the landscape of TIME further optimizing precision immunotherapy.

Background: Hepatocellular carcinoma (HCC) was the sixth common malignancies characteristic with highly aggressive in the world. It was well established that tumor mutation burden (TMB) act as indicator of immunotherapeutic responsiveness in various tumors. However, the role of TMB in tumor immune microenvironment (TIME) is still obscure. Method: The mutation data was analyzed by employing “maftools” package. Weighted gene co-expression network analysis (WGCNA) was implemented to determine candidate module and significant genes correlated with TMB value. Differential analysis was performed between different level of TMB subgroups employing R package “limma”. Gene ontology (GO) enrichment analysis was implemented with “clusterProfiler”, “enrichplot” and “ggplot2” packages. Then risk score signature was developed by systematical bioinformatics analyses. K-M survival curves and receiver operating characteristic (ROC) plot were further analyzed for prognostic validity. To depict comprehensive context of TIME, XCELL, TIMER, QUANTISEQ, MCPcounter, EPIC, CIBERSORT, and CIBERSORT-ABS algorithm were employed. Additionally, the potential role of risk score on immune checkpoint blockade (ICB) immunotherapy was further explored. The quantitative real-time polymerase chain reaction was performed to detect expression of HTRA3. Results: TMB value was positively correlated with older age, male gender and early T status. A total of 75 intersection genes between TMB-related genes and differentially expressed genes (DEGs) were screened and enriched in extracellular matrix-relevant pathways. Risk score based on three hub genes significantly affected overall survival (OS) time, infiltration of immune cells, and ICB-related hub targets. The prognostic performance of risks score was validated in the external testing group. Risk-clinical nomogram was constructed for clinical application. HTRA3 was demonstrated to be a prognostic factor in HCC in further exploration. Finally, mutation of TP53 was correlated with risk score and do not interfere with risk score-based prognostic prediction. Conclusion: Collectively, a comprehensive analysis of TMB might provide novel insights into mutation-driven mechanism of tumorigenesis further contribute to tailored immunotherapy and prognosis prediction of HCC.

Pancreatic ductal adenocarcinoma (PDAC) is the most frequent and aggressive pancreatic tumor characterized by high metastatic risk and special tumor microenvironment. To comprehensively delineate the complex intra-tumoral heterogeneity and the underlying mechanism during metastatic lesions malignant progression, single-cell RNA sequencing (scRNA-seq) was employed. PCA and TSNE were used for dimensionality reduction analysis and cell clustering. Find All Markers function was used to calculate differential genes in each cluster, and Do Heatmap function was used to plot the distribution of differential genes in each cluster. GSVA was employed to assign pathway activity estimates to individual cells. Lineage trajectory progression was inferred by monocle. CNV status was inferred to compare the heterogeneity among patients and subtypes by infercnv. Ligand-receptor interactions were identified by CellPhoneDB, and regulons network of cells was analyzed by SCENIC. Through RNA-sequencing of 6236 individual cells from 5 liver metastatic PDAC lesions, 10 major cell clusters are identified by using unbiased clustering analysis of expression profiling and well-known cell markers. Cells with high CNV level were considered as malignant cells and pathway analyses were carried out to highlight intratumor heterogeneity in PDAC. Pseudotime trajectory analysis revealed that components of multiple tumor-related pathways and transcription factors (TFs) were differentially expressed along PDAC progression. The complex cellular communication suggested potential immunotherapeutic targets in PDAC. Regulon network identified multiple candidates for promising cell-specific transcriptional factors. Finally, metastatic-related genes expression levels and signaling pathways were validated in bulk RNA Sequencing data. This study contributed a comprehensive single-cell transcriptome atlas and contributed into novel insight of intratumor heterogeneity and molecular mechanism in metastatic PDAC.

Background: Liver cancer stem cells, characterized by self-renewal and initiating cancer cells, were decisive drivers of progression and therapeutic resistance in hepatocellular carcinoma (HCC). However, a comprehensive understanding of HCC stemness has not been identified.Methods: RNA sequencing information, corresponding clinical annotation, and mutation data of HCC were downloaded from The Cancer Genome Atlas-LIHC project. Two stemness indices, mRNA expression-based stemness index (mRNAsi) and epigenetically regulated-mRNAsi, were used to comprehensively analyze HCC stemness. Estimation of Stromal and Immune Cells in Malignant Tumors using Expression Data and single-sample gene-set enrichment analysis algorithm were performed to characterize the context of tumor immune microenvironment (TIME). Next, differentially expressed gene (DEG) analysis and weighted gene co-expression network analysis (WGCNA) were employed to identify significant mRNAsi-related modules with hub genes. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology enrichment pathways were analyzed to functionally annotate these key genes. The least absolute shrinkage and selection operator (LASSO) Cox regression analysis was performed to establish a prognostic signature. Kaplan–Meier survival curves and receiver operating characteristic (ROC) analysis were applied for prognostic value validation. Seven algorithms (XCELL, TIMER, QUANTISEQ, MCPcounter, EPIC, CIBERSORT, and CIBERSORT-ABS) were utilized to draw the landscape of TIME. Finally, the mutation data were analyzed by employing “maftools” package.Results: mRNAsi was significantly elevated in HCC samples. mRNAsi escalated as tumor grade increased, with poor prognosis presenting the higher stemness index. The stemness-related (greenyellow) modules with 175 hub genes were screened based on DEGs and WGCNA. A prognostic signature was established using LASSO analysis of prognostic hub genes to classify samples into two risk subgroups, which exhibited good prognostic performance. Additionally, prognostic risk-clinical nomogram was drawn to estimate risk quantitatively. Moreover, risk score was significantly associated with contexture of TIME and immunotherapeutic targets. Finally, potential interaction between risk score with tumor mutation burden (TMB) was elucidated.Conclusion: This work comprehensively elucidated that stemness characteristics served as a crucial player in clinical outcome, complexity of TIME, and immunotherapeutic prediction from both mRNAsi and mRNA level. Quantitative identification of stemness characteristics in individual tumor will contribute into predicting clinical outcome, mapping landscape of TIME further optimizing precision immunotherapy.

Objective: To investigate the efficacy and safety of febuxostat on renal function in CKD stage 3 diabetic nephropathy patients. Methods: Patients in our hospital with chronic kidney disease (CKD) stage 3 diabetic nephropathy (DN) complicated by high serum uric acid (360 µmol/L) were recruited. Patients were then divided into treatment group and control group according to the random number table method. All the patients received low purine diet, renin-angiotensin-aldosterone system (RAAS) inhibitors, and adequate routine hypoglycemic treatment. Febuxostat was employed only in the treatment group. The levels of blood uric acid (sUA), serum creatinine (Scr), cystatin C (cys-c), eGFR, 24-hour urine protein quantification, albuminuria, and creatinine ratio (ACR) were evaluated in all patients before and after treatment at 4, 8, 12, and 24 week. Results: No difference was found before treatment between the two groups. After treatment at 4, 8, 12, and 24 week, the levels of sUA, SCr, cys-c, and eGFR between the two groups were significant different (P<0.05). There was no difference in 24-hour urine protein quantification, albuminuria, and creatinine ratio between two groups before treatment, and significant differences were observed after treatment. Fifty percent of patients from the treatment group achieved the treatment goal with 20 mg febuxostat at 4 weeks. Tubular markers were also decreased with the treatment. Conclusions: Febuxostat can reduce uric acid and improve renal function effectively in patients with CKD stage 3 diabetic nephropathy, while being well tolerated. However, the conclusion is still uncertain due to the short term of the study.