Key driver genes as potential therapeutic targets in renal allograft rejection

LILRB3 genetic variation is associated with kidney transplant failure in African American recipients

African American (AA) kidney transplant recipients exhibit a higher rate of graft loss compared with other racial and ethnic populations, highlighting the need to identify causative factors. Here, in the Genomics of Chronic Allograft Rejection cohort, pretransplant blood RNA sequencing revealed a cluster of four consecutive missense single-nucelotide polymorphisms (SNPs), within the leukocyte immunoglobulin-like receptor B3 (LILRB3) gene, strongly associated with death-censored graft loss. This SNP cluster (named LILRB3-4SNPs) encodes missense mutations at amino acids 617-618 proximal to a SHP1/2 phosphatase-binding immunoreceptor tyrosine-based inhibitory motif. The LILRB3-4SNPs cluster is specifically enriched within AA individuals and exhibited a strong association with death-censored graft loss and estimated glomerular filtration rate decline in the AA participants from multiple transplant cohorts. In two large Biobanks (BioMe and All-of-Us), the LILRB3-4SNPs cluster was associated with the early onset of end-stage renal disease and acted synergistically with the apolipoprotein L1 (APOL1) G1/G2 allele to accelerate disease progression. The SNPs were also linked to multiple immune-related diseases in AA individuals. Last, on multiomics analysis of blood and biopsies, recipients with LILRB3-4SNPs showed enhanced inflammation and monocyte ferroptosis. While larger and prospective studies are needed, our data provide insights on the genetic variation underlying kidney transplant outcomes.

Decoding the hallmarks of allograft dysfunction with a comprehensive pan-organ transcriptomic atlas

The pathogenesis of allograft (dys)function has been increasingly studied using 'omics'-based technologies, but the focus on individual organs has created knowledge gaps that neither unify nor distinguish pathological mechanisms across allografts. Here we present a comprehensive study of human pan-organ allograft dysfunction, analyzing 150 datasets with more than 12,000 samples across four commonly transplanted solid organs (heart, lung, liver and kidney, n = 1,160, 1,241, 1,216 and 8,853 samples, respectively) that we leveraged to explore transcriptomic differences among allograft dysfunction (delayed graft function, acute rejection and fibrosis), tolerance and stable graft function. We identified genes that correlated robustly with allograft dysfunction across heart, lung, liver and kidney transplantation. Furthermore, we developed a transfer learning omics prediction framework that, by borrowing information across organs, demonstrated superior classifications compared to models trained on single organs. These findings were validated using a single-center prospective kidney transplant cohort study (a collective 329 samples across two timepoints), providing insights supporting the potential clinical utility of our approach. Our study establishes the capacity for machine learning models to learn across organs and presents a transcriptomic transplant resource that can be employed to develop pan-organ biomarkers of allograft dysfunction.

Improving liver transplant outcomes with transplant-omics and network biology

PURPOSE OF REVIEW

Molecular omics data is increasingly ubiquitous throughout medicine. In organ transplantation, recent large-scale research efforts are generating the 'transplant-ome' - the entire set of molecular omics data, including the genome, transcriptome, proteome, and metabolome. Importantly, early studies in anesthesiology have demonstrated how perioperative interventions alter molecular profiles in various patient populations. The next step for anesthesiologists and intensivists will be to tailor perioperative care to the transplant-ome of individual liver transplant patients.

RECENT FINDINGS

In liver transplant patients, elements of the transplant-ome predict complications and point to novel interventions. Importantly, molecular profiles of both the donor organ and recipient contribute to this risk, and interventions like normothermic machine perfusion influence these profiles. As we can now measure various omics molecules simultaneously, we can begin to understand how these molecules interact to form molecular networks and emerging technologies offer noninvasive and continuous ways to measure these networks throughout the perioperative period. Molecules that regulate these networks are likely mediators of complications and actionable clinical targets throughout the perioperative period.

SUMMARY

The transplant-ome can be used to tailor perioperative care to the individual liver transplant patient. Monitoring molecular networks continuously and noninvasively would provide new opportunities to optimize perioperative management.

Female Gene Networks Are Expressed in Myofibroblast-Like Smooth Muscle Cells in Vulnerable Atherosclerotic Plaques

BACKGROUND

Women presenting with coronary artery disease more often present with fibrous atherosclerotic plaques, which are currently understudied. Phenotypically modulated smooth muscle cells (SMCs) contribute to atherosclerosis in women. How these phenotypically modulated SMCs shape female versus male plaques is unknown.

METHODS

Gene regulatory networks were created using RNAseq gene expression data from human carotid atherosclerotic plaques. The networks were prioritized based on sex bias, relevance for smooth muscle biology, and coronary artery disease genetic enrichment. Network expression was linked to histologically determined plaque phenotypes. In addition, their expression in plaque cell types was studied at single-cell resolution using single-cell RNAseq. Finally, their relevance for disease progression was studied in female and male Apoe-/- mice fed a Western diet for 18 and 30 weeks.

RESULTS

Here, we identify multiple sex-stratified gene regulatory networks from human carotid atherosclerotic plaques. Prioritization of the female networks identified 2 main SMC gene regulatory networks in late-stage atherosclerosis. Single-cell RNA sequencing mapped these female networks to 2 SMC phenotypes: a phenotypically modulated myofibroblast-like SMC network and a contractile SMC network. The myofibroblast-like network was mostly expressed in plaques that were vulnerable in women. Finally, the mice ortholog of key driver gene MFGE8 (milk fat globule EGF and factor V/VIII domain containing) showed retained expression in advanced plaques from female mice but was downregulated in male mice during atherosclerosis progression.

CONCLUSIONS

Female atherosclerosis is characterized by gene regulatory networks that are active in fibrous vulnerable plaques rich in myofibroblast-like SMCs.

Functional Genome Analysis for Immune Cells Provides Clues for Stratification of Systemic Lupus Erythematosus

Systemic lupus erythematosus (SLE) is caused by a combination of genetic and environmental factors. Recently, analysis of a functional genome database of genetic polymorphisms and transcriptomic data from various immune cell subsets revealed the importance of the oxidative phosphorylation (OXPHOS) pathway in the pathogenesis of SLE. In particular, activation of the OXPHOS pathway is persistent in inactive SLE, and this activation is associated with organ damage. The finding that hydroxychloroquine (HCQ), which improves the prognosis of SLE, targets toll-like receptor (TLR) signaling upstream of OXPHOS suggests the clinical importance of this pathway. IRF5 and SLC15A4, which are regulated by polymorphisms associated with SLE susceptibility, are functionally associated with OXPHOS as well as blood interferon activity and metabolome. Future analyses of OXPHOS-associated disease-susceptibility polymorphisms, gene expression, and protein function may be useful for risk stratification of SLE.

Characteristic Genes and Immune Infiltration Analysis for Acute Rejection after Kidney Transplantation

Background

Renal transplantation can significantly improve the survival rate and quality of life of patients with end-stage renal disease, but the probability of acute rejection (AR) in adult renal transplant recipients is still approximately 12.2%. Machine learning (ML) is superior to traditional statistical methods in various clinical scenarios. However, the current AR model is constructed only through simple difference analysis or a single queue, which cannot guarantee the accuracy of prediction. Therefore, this study identified and validated new gene sets that contribute to the early prediction of AR and the prognosis prediction of patients after renal transplantation by constructing a more accurate AR gene signature through ML technology.

Methods

Based on the Gene Expression Omnibus (GEO) database and multiple bioinformatic analyses, we identified differentially expressed genes (DEGs) and built a gene signature via LASSO regression and SVM analysis. Immune cell infiltration and immunocyte association analyses were also conducted. Furthermore, we investigated the relationship between AR genes and graft survival status.

Results

Twenty-four DEGs were identified. A 5 gene signature (CPA6, EFNA1, HBM, THEM5, and ZNF683) were obtained by LASSO analysis and SVM analysis, which had a satisfied ability to differentiate AR and NAR in the training cohort, internal validation cohort and external validation cohort. Additionally, ZNF683 was associated with graft survival.

Conclusion

A 5 gene signature, particularly ZNF683, provided insight into a precise therapeutic schedule and clinical applications for AR patients.

SLAMF8 Participates in Acute Renal Transplant Rejection via TLR4 Pathway on Pro-Inflammatory Macrophages

Background

Acute rejection (AR) in kidney transplantation is an established risk factor that reduces the survival rate of allografts. Despite standard immunosuppression, molecules with regulatory control in the immune pathway of AR can be used as important targets for therapeutic operations to prevent rejection.

Methods

We downloaded the microarray data of 15 AR patients and 37 non-acute rejection (NAR) patients from Gene Expression Omnibus (GEO). Gene network was constructed, and genes were classified into different modules using weighted gene co-expression network analysis (WGCNA). Kyoto Encyclopedia of Genes and Genomes (KEGG) and Cytoscape were applied for the hub genes in the most related module to AR. Different cell types were explored by xCell online database and single-cell RNA sequencing. We also validated the SLAMF8 and TLR4 levels in Raw264.7 and human kidney tissues of TCMR.

Results

A total of 1,561 differentially expressed genes were filtered. WGCNA was constructed, and genes were classified into 12 modules. Among them, the green module was most closely associated with AR. These genes were significantly enriched in 20 pathway terms, such as cytokine–cytokine receptor interaction, chemokine signaling pathway, and other important regulatory processes. Intersection with GS > 0.4, MM > 0.9, the top 10 MCC values and DEGs in the green module, and six hub genes (DOCK2, NCKAP1L, IL2RG, SLAMF8, CD180, and PTPRE) were identified. Their expression levels were all confirmed to be significantly elevated in AR patients in GEO, Nephroseq, and quantitative real-time PCR (qRT-PCR). Single-cell RNA sequencing showed that AR patient had a higher percentage of native T, CD1C+_B DC, NKT, NK, and monocytes in peripheral blood mononuclear cells (PBMCs). Xcell enrichment scores of 20 cell types were significantly different (p<0.01), mostly immune cells, such as B cells, CD4+ Tem, CD8+ T cells, CD8+ Tcm, macrophages, M1, and monocytes. GSEA suggests that highly expressed six hub genes are correlated with allograft rejection, interferon γ response, interferon α response, and inflammatory response. In addition, SLAMF8 is highly expressed in human kidney tissues of TCMR and in M1 phenotype macrophages of Raw264.7 cell line WGCNA accompanied by high expression of TLR4.

Conclusion

This study demonstrates six hub genes and functionally enriched pathways related to AR. SLAMF8 is involved in the M1 macrophages via TLR4, which contributed to AR process.

Guidelines for bioinformatics of single-cell sequencing data analysis in Alzheimer's disease: review, recommendation, implementation and application

Alzheimer's disease (AD) is the most common form of dementia, characterized by progressive cognitive impairment and neurodegeneration. Extensive clinical and genomic studies have revealed biomarkers, risk factors, pathways, and targets of AD in the past decade. However, the exact molecular basis of AD development and progression remains elusive. The emerging single-cell sequencing technology can potentially provide cell-level insights into the disease. Here we systematically review the state-of-the-art bioinformatics approaches to analyze single-cell sequencing data and their applications to AD in 14 major directions, including 1) quality control and normalization, 2) dimension reduction and feature extraction, 3) cell clustering analysis, 4) cell type inference and annotation, 5) differential expression, 6) trajectory inference, 7) copy number variation analysis, 8) integration of single-cell multi-omics, 9) epigenomic analysis, 10) gene network inference, 11) prioritization of cell subpopulations, 12) integrative analysis of human and mouse sc-RNA-seq data, 13) spatial transcriptomics, and 14) comparison of single cell AD mouse model studies and single cell human AD studies. We also address challenges in using human postmortem and mouse tissues and outline future developments in single cell sequencing data analysis. Importantly, we have implemented our recommended workflow for each major analytic direction and applied them to a large single nucleus RNA-sequencing (snRNA-seq) dataset in AD. Key analytic results are reported while the scripts and the data are shared with the research community through  GitHub. In summary, this comprehensive review provides insights into various approaches to analyze single cell sequencing data and offers specific guidelines for study design and a variety of analytic directions. The review and the accompanied software tools will serve as a valuable resource for studying cellular and molecular mechanisms of AD, other diseases, or biological systems at the single cell level.