Transcriptome- and proteome-wide association studies nominate determinants of kidney function and damage

Integrative analysis of transcriptome and proteome wide association studies prioritized functional genes for obesity

BACKGROUND

Genome-wide association studies have identified dozens of genomic loci for obesity. However, functional genes and their detailed genetic mechanisms underlying these loci are mainly unknown. In this study, we conducted an integrative study to prioritize plausibly functional genes by combining information from genome-, transcriptome- and proteome-wide association analyses.

METHODS

We first conducted proteome-wide association analyses and transcriptome-wide association analyses for the six obesity-related traits. We then performed colocalization analysis on the identified loci shared between the proteome- and transcriptome-association analyses. Finally, we validated the identified genes with other plasma/blood reference panels. The highlighted genes were assessed for expression of other tissues, single-cell and tissue specificity, and druggability.

RESULTS

We prioritized 4 high-confidence genes (FASN, ICAM1, PDCD6IP, and YWHAB) by proteome-wide association studies, transcriptome-wide association studies, and colocalization analyses, which consistently influenced the variation of obesity traits at both mRNA and protein levels. These 4 genes were successfully validated using other plasma/blood reference panels. These 4 genes shared regulatory structures in obesity-related tissues. Single-cell and tissue-specific analyses showed that FASN and ICAM1 were explicitly expressed in metabolism- and immunity-related tissues and cells. Furthermore, FASN and ICAM1 had been developed as drug targets.

CONCLUSION

Our study provided novel promising protein targets for further mechanistic and therapeutic studies of obesity.

Multi-omic characterization of air pollution effects: Applications of AirSigOmniTWP Hub

Air pollution, particularly fine particulate matter and gaseous pollutants including NO2 and NOx, presents significant public health challenges. While the harmful effects of these pollutants are well-documented, the molecular mechanisms underlying their impact on health remain incompletely understood. In this study, we utilized genome-wide association study (GWAS) data from the UK Biobank, expression quantitative trait loci (eQTL) data from the Genotype-Tissue Expression (GTEx) project, and protein quantitative trait loci (pQTL) data from the Atherosclerosis Risk in Communities (ARIC) study to conduct comprehensive analyses using the Unified Test for Molecular Signatures (UTMOST), Transcriptome-wide Association Studies (TWAS), and Proteome-wide Association Studies (PWAS). To integrate and synthesize these analyses, we developed the AirSigOmniTWP Hub, a specialized platform designed to consolidate and interpret the results from UTMOST, TWAS, and PWAS. TWAS analysis identified a significant association between PM10 exposure and the gene INO80E in females (P = 4.37×10⁻⁵, FDR = 0.0383), suggesting a potential role in chromatin remodeling. PWAS analysis revealed a significant association between NOx exposure and the gene PIP in females (P = 2.28×10⁻⁵, FDR = 0.0299), implicating its involvement in inflammatory pathways. Additionally, UTMOST analyses uncovered significant associations between various pollutants and genes including NCOA4P3 and SPATS2L with PM2.5 exposure, indicating potential mechanisms related to transcriptional regulation and gene-environment interactions.

Brief research report: WGCNA-driven identification of histone modification genes as potential biomarkers in AQP4-Associated optic neuritis

Introduction

Neuromyelitis Optica spectrum disorder (NMOSD) is an autoimmune disease characterized by anti-aquaporin-4 (AQP4) auto-antibodies. The discovery of antibodies AQP4 and myelin oligodendrocyte glycoprotein (MOG) has expanded our understanding of the pathogenesis of neuromyelitis optica. However, the molecular mechanisms underlying the disease, particularly AQP4-associated optic neuritis (AQP4-ON), remain to be fully elucidated.

Methods

In this study, we utilized Weighted Gene Co-expression Network Analysis (WGCNA) to investigate the transcriptomic profiles of peripheral blood samples from patients with AQP4-ON and MOG-positive optic neuritis (MOG-ON), compared to healthy controls.

Results

WGCNA revealed a brown module (ME brown) strongly associated with AQP4-ON, which correlated positively with post-onset visual acuity decline. A total of 132 critical genes were identified, mainly involved in histone modification and microtubule dynamics. Notably, genes HDAC4, HDAC7, KDM6A, and KDM5C demonstrated high AUC values in ROC analysis, indicating their potential as biomarkers for AQP4-ON.

Conclusion

Our findings provide novel insights into the molecular signature of AQP4-ON and highlight the potential of systems biology approaches in identifying biomarkers for NMOSD. The identified histone modification genes warrant further investigation for their role in disease pathogenesis and as therapeutic targets.

Association of Integrated Proteomic and Metabolomic Modules with Risk of Kidney Disease Progression

Key Points

Integrated analysis of proteome and metabolome identifies modules associated with CKD progression and kidney failure. Ephrin transmembrane proteins and podocyte-expressed CRIM1 and NPNT emerged as central components and warrant experimental and clinical investigation.

Background

Proteins and metabolites play crucial roles in various biological functions and are frequently interconnected through enzymatic or transport processes.

Methods

We present an integrated analysis of 4091 proteins and 630 metabolites in the Chronic Renal Insufficiency Cohort study (N =1708; average follow-up for kidney failure, 9.5 years, with 537 events). Proteins and metabolites were integrated using an unsupervised clustering method, and we assessed associations between clusters and CKD progression and kidney failure using Cox proportional hazards models. Analyses were adjusted for demographics and risk factors, including the eGFR and urine protein–creatinine ratio. Associations were identified in a discovery sample (random two thirds, n =1139) and then evaluated in a replication sample (one third, n =569).

Results

We identified 139 modules of correlated proteins and metabolites, which were represented by their principal components. Modules and principal component loadings were projected onto the replication sample, which demonstrated a consistent network structure. Two modules, representing a total of 236 proteins and 82 metabolites, were robustly associated with both CKD progression and kidney failure in both discovery and validation samples. Using gene set enrichment, several transmembrane-related terms were identified as overrepresented in these modules. Transmembrane–ephrin receptor activity displayed the largest odds (odds ratio=13.2, P value = 5.5×10−5). A module containing CRIM1 and NPNT expressed in podocytes demonstrated particularly strong associations with kidney failure (P value = 2.6×10−5).

Conclusions

This study demonstrates that integration of the proteome and metabolome can identify functions of pathophysiologic importance in kidney disease.

Identification of novel therapeutic targets for chronic kidney disease and kidney function by integrating multi-omics proteome with transcriptome

BACKGROUND

Chronic kidney disease (CKD) is a progressive disease for which there is no effective cure. We aimed to identify potential drug targets for CKD and kidney function by integrating plasma proteome and transcriptome.

METHODS

We designed a comprehensive analysis pipeline involving two-sample Mendelian randomization (MR) (for proteins), summary-based MR (SMR) (for mRNA), and colocalization (for coding genes) to identify potential multi-omics biomarkers for CKD and combined the protein-protein interaction, Gene Ontology (GO), and single-cell annotation to explore the potential biological roles. The outcomes included CKD, extensive kidney function phenotypes, and different CKD clinical types (IgA nephropathy, chronic glomerulonephritis, chronic tubulointerstitial nephritis, membranous nephropathy, nephrotic syndrome, and diabetic nephropathy).

RESULTS

Leveraging pQTLs of 3032 proteins from 3 large-scale GWASs and corresponding blood- and tissue-specific eQTLs, we identified 32 proteins associated with CKD, which were validated across diverse CKD datasets, kidney function indicators, and clinical types. Notably, 12 proteins with prior MR support, including fibroblast growth factor 5 (FGF5), isopentenyl-diphosphate delta-isomerase 2 (IDI2), inhibin beta C chain (INHBC), butyrophilin subfamily 3 member A2 (BTN3A2), BTN3A3, uromodulin (UMOD), complement component 4A (C4a), C4b, centrosomal protein of 170 kDa (CEP170), serologically defined colon cancer antigen 8 (SDCCAG8), MHC class I polypeptide-related sequence B (MICB), and liver-expressed antimicrobial peptide 2 (LEAP2), were confirmed. To our knowledge, 20 novel causal proteins have not been previously reported. Five novel proteins, namely, GCKR (OR 1.17, 95% CI 1.10-1.24), IGFBP-5 (OR 0.43, 95% CI 0.29-0.62), sRAGE (OR 1.14, 95% CI 1.07-1.22), GNPTG (OR 0.90, 95% CI 0.86-0.95), and YOD1 (OR 1.39, 95% CI 1.18-1.64,) passed the MR, SMR, and colocalization analysis. The other 15 proteins were also candidate targets (GATM, AIF1L, DQA2, PFKFB2, NFATC1, activin AC, Apo A-IV, MFAP4, DJC10, C2CD2L, TCEA2, HLA-E, PLD3, AIF1, and GMPR1). These proteins interact with each other, and their coding genes were mainly enrichment in immunity-related pathways or presented specificity across tissues, kidney-related tissue cells, and kidney single cells.

CONCLUSIONS

Our integrated analysis of plasma proteome and transcriptome data identifies 32 potential therapeutic targets for CKD, kidney function, and specific CKD clinical types, offering potential targets for the development of novel immunotherapies, combination therapies, or targeted interventions.

Plasma proteome analysis implicates novel proteins as potential therapeutic targets for chronic kidney disease: A proteome-wide association study

Chronic kidney disease (CKD) is prevalent globally with limited therapeutic drugs available. To systemically identify novel proteins involved in the pathogenesis of CKD and possible therapeutic targets, we integrated human plasma proteomes with the genome-wide association studies (GWASs) of CKD, estimated glomerular filtration rate (eGFR) and blood urea nitrogen (BUN) to perform proteome-wide association study (PWAS), Mendelian Randomization and Bayesian colocalization analyses. The single-cell RNA sequencing data of healthy human and mouse kidneys were analyzed to explore the cell-type specificity of identified genes. Functional enrichment analysis was conducted to investigate the involved signaling pathways. The PWAS identified 22 plasma proteins significantly associated with CKD. Of them, the significant associations of three proteins (INHBC, LMAN2, and SNUPN) were replicated in the GWASs of eGFR, and BUN. Mendelian Randomization analyses showed that INHBC and SNUPN were causally associated with CKD, eGFR, and BUN. The Bayesian colocalization analysis identified shared causal variants for INHBC in CKD, eGFR, and BUN (all PP4 > 0.75). The single-cell RNA sequencing revealed that the INHBC gene was sparsely scattered within the kidney cells. This proteomic study revealed that INHBC, LMAN2, and SNUPN may be involved in the pathogenesis of CKD, which represent novel therapeutic targets and warrant further exploration in future research.

Genetic imputation of kidney transcriptome, proteome and multi-omics illuminates new blood pressure and hypertension targets

Genetic mechanisms of blood pressure (BP) regulation remain poorly defined. Using kidney-specific epigenomic annotations and 3D genome information we generated and validated gene expression prediction models for the purpose of transcriptome-wide association studies in 700 human kidneys. We identified 889 kidney genes associated with BP of which 399 were prioritised as contributors to BP regulation. Imputation of kidney proteome and microRNAome uncovered 97 renal proteins and 11 miRNAs associated with BP. Integration with plasma proteomics and metabolomics illuminated circulating levels of myo-inositol, 4-guanidinobutanoate and angiotensinogen as downstream effectors of several kidney BP genes (SLC5A11, AGMAT, AGT, respectively). We showed that genetically determined reduction in renal expression may mimic the effects of rare loss-of-function variants on kidney mRNA/protein and lead to an increase in BP (e.g., ENPEP). We demonstrated a strong correlation (r = 0.81) in expression of protein-coding genes between cells harvested from urine and the kidney highlighting a diagnostic potential of urinary cell transcriptomics. We uncovered adenylyl cyclase activators as a repurposing opportunity for hypertension and illustrated examples of BP-elevating effects of anticancer drugs (e.g. tubulin polymerisation inhibitors). Collectively, our studies provide new biological insights into genetic regulation of BP with potential to drive clinical translation in hypertension.

Discovery of genomic and transcriptomic pleiotropy between kidney function and soluble receptor for advanced glycation end-products using correlated meta-analyses: The Long Life Family Study (LLFS)

Patients with chronic kidney disease (CKD) have increased oxidative stress and chronic inflammation, which may escalate the production of advanced glycation end-products (AGE). High soluble receptor for AGE (sRAGE) and low estimated glomerular filtration rate (eGFR) levels are associated with CKD and aging. We evaluated whether eGFR calculated from creatinine and cystatin C share pleiotropic genetic factors with sRAGE. We employed whole-genome sequencing and correlated meta-analyses on combined genomewide association study (GWAS) p -values in 4,182 individuals (age range: 24-110) from the Long Life Family Study (LLFS). We also conducted transcriptome-wide association studies (TWAS) on whole blood in a subset of 1,209 individuals. We identified 59 pleiotropic GWAS loci ( p <5×10 -8 ) and 17 TWAS genes (Bonferroni- p <2.73×10 -6 ) for eGFR traits and sRAGE. TWAS genes, LSP1 and MIR23AHG , were associated with eGFR and sRAGE located within GWAS loci, lncRNA- KCNQ1OT1 and CACNA1A/CCDC130 , respectively. GWAS variants were eQTLs in the kidney glomeruli and tubules, and GWAS genes predicted kidney carcinoma. TWAS genes harbored eQTLs in the kidney, predicted kidney carcinoma, and connected enhancer-promoter variants with kidney function-related phenotypes at p <5×10 -8 . Additionally, higher allele frequencies of protective variants for eGFR traits were detected in LLFS than in ALFA-Europeans and TOPMed, suggesting better kidney function in healthy-aging LLFS than in general populations. Integrating genomic annotation and transcriptional gene activity revealed the enrichment of genetic elements in kidney function and kidney diseases. The identified pleiotropic loci and gene expressions for eGFR and sRAGE suggest their underlying shared genetic effects and highlight their roles in kidney- and aging-related signaling pathways.

Advancing proteomics in nephrology: unraveling causal pathways and therapeutic targets

Proteomics has illuminated disease pathophysiology, unearthed novel biomarkers, and bolstered risk assessment strategies. In nephrology, observational analyses unveil biomarkers associated with adverse outcomes, whereas genetics offer insights into causal pathways. Mendelian randomization offers the potential to link the two, uncovering causal relationships between biomarkers and kidney function. Lanktree et al. demonstrate Mendelian randomization's utility in identifying additional proteins affecting kidney function and kidney disease progression.