Tianlong Zhang
Ying Cui
Siyi Jiang1
Lijun Song
Min Li- 1Department of Critical Care Medicine, The Fourth Affiliated Hospital of School of Medicine, and International School of Medicine, International Institutes of Medicine, Zhejiang University, Yiwu, China
- 2Department of Critical Care Medicine, Jinhua Hospital Affiliated to Zhejiang University, Jinhua, Zhejiang, China
Background: Clinical studies have indicated a comorbidity between sepsis and kidney diseases. Individuals with specific mutations that predispose them to kidney conditions are also at an elevated risk for developing sepsis, and vice versa. This suggests a potential shared genetic etiology that has not been fully elucidated.
Methods: Summary statistics data on exposure and outcomes were obtained from genome-wide association meta-analysis studies. We utilized these data to assess genetic correlations, employing a pleiotropy analysis method under the composite null hypothesis to identify pleiotropic loci. After mapping the loci to their corresponding genes, we conducted pathway analysis using Generalized Gene-Set Analysis of GWAS Data (MAGMA). Additionally, we utilized MAGMA gene-test and eQTL information (whole blood tissue) for further determination of gene involvement. Further investigation involved stratified LD score regression, using diverse immune cell data, to study the enrichment of SNP heritability in kidney-related diseases and sepsis. Furthermore, we employed Mendelian Randomization (MR) analysis to investigate the causality between kidney diseases and sepsis.
Results: In our genetic correlation analysis, we identified significant correlations among BUN, creatinine, UACR, serum urate, kidney stones, and sepsis. The PLACO analysis method identified 24 pleiotropic loci, pinpointing a total of 28 nearby genes. MAGMA gene-set enrichment analysis revealed a total of 50 pathways, and tissue-specific analysis indicated significant enrichment of five pairs of pleiotropic results in kidney tissue. MAGMA gene test and eQTL information (whole blood tissue) identified 33 and 76 pleiotropic genes, respectively. Notably, genes PPP2R3A for BUN, VAMP8 for UACR, DOCK7 for creatinine, and HIBADH for kidney stones were identified as shared risk genes by all three methods. In a series of immune cell-type-specific enrichment analyses of pleiotropy, we identified a total of 37 immune cells. However, MR analysis did not reveal any causal relationships among them.
Conclusions: This study lays the groundwork for shared etiological factors between kidney and sepsis. The confirmed pleiotropic loci, shared pathogenic genes, and enriched pathways and immune cells have enhanced our understanding of the multifaceted relationships among these diseases. This provides insights for early disease intervention and effective treatment, paving the way for further research in this field.
1 Background
Sepsis is a life-threatening disease that endangers human lives and results in multiple organ dysfunction (1). It remains one of the leading causes of morbidity and mortality worldwide, highlighting the critical need to understand its impact on public health (2). In 2017 alone, an estimated 48.9 million cases of sepsis occurred globally, resulting in 11 million deaths, accounting for nearly 20% of all global deaths (3). Due to limited treatment options, identifying the underlying mechanisms of sepsis is essential for developing novel therapeutic strategies.
Sepsis presents with a heterogeneous range of clinical symptoms (4), with renal dysfunction being a frequent complication. Early stages of various kidney diseases often manifest as renal impairment, characterized by changes in blood creatinine, blood urea nitrogen (BUN), urinary albumin-to-creatinine ratio (UACR), uric acid (UA), estimated glomerular filtration rate (eGFR), and other relevant indicators. For instance, elevated blood creatinine levels are frequently observed in sepsis patients, a phenomenon also characteristic of acute kidney injury and chronic kidney disease (5). Similarly, changes in BUN, UACR, and UA have been closely linked to sepsis (6–8). Additionally, studies suggest an association between stone formation and adverse outcomes like sepsis (9). The connections between IgA nephropathy, diabetic nephropathy, renal tumors, and sepsis also require further elucidation. The shared genetic underpinnings of these associations remain unexplored. Therefore, a systematic analysis is crucial to investigate the presence of common pleiotropic risk variants between kidney-related diseases or relevant indicators and sepsis, along with the potential involvement of specific biological pathways and shared risk genes.
The rise of genome-wide association studies (GWAS) in recent years has been instrumental. GWAS have become a powerful tool for identifying a wide range of genetic variations associated with complex diseases and elucidating genes related to the occurrence, progression, and treatment of diseases (10). In our GWAS study, patients with acute kidney injury and chronic kidney disease predominantly exhibited elevated indicators such as creatinine and eGFR. Therefore, we used these kidney-related indicators as proxies for the two conditions. Utilizing large-scale GWAS summary data, we conducted a genome-wide analysis for nine kidney-related diseases or indicators (BUN, UACR, UA, eGFR, creatine, kidney neoplasm, kidney stone, IgA nephropathy, and diabetic nephropathy) and sepsis. We employed various genetic methods, including pleiotropy analysis for significant results, to systematically explore pleiotropic associations at both the gene and pathway levels. This analysis aimed to uncover potential shared genetic etiology. Additionally, we conducted gene set enrichment analysis and immune cell enrichment analysis to discover potential correlations. Finally, Mendelian randomization (MR) analysis was employed to explore causal relationships among these factors.
2 Methods
2.1 Study design
This study investigates the potential genetic associations between kidney disease or relevant kidney function indicators and sepsis. A schematic representation of the study design is provided in Figure 1.
Figure 1 An overview of the study design.
2.2 GWAS data sources
In our analyses, we investigated the genetic associations between the following traits: BUN (N=480,698), UACR (N=547,361), UA (N=110,347), eGFR (N=480,698), creatinine (N=110,051), kidney neoplasm (N=456,348), kidney stone disease (N=456,348), IgA nephropathy (N=477,784), diabetic nephropathy (N=452,280), and sepsis (N=344,644) (11–15). Details regarding the data sources for the GWAS employed in this analysis are presented in Supplementary Table S1. All studies included in these consortia were approved by local research ethics committees and institutional review boards, and written informed consent was obtained from all participants.
2.3 Statistical analyses
2.3.1 Genetic correlation analysis
Linkage Disequilibrium Score Regression (LDSC) was employed to assess shared polygenic structures between traits, with LD scores calculated from European ancestry samples of the 1000 Genomes Project serving as the reference group (1). Stringent quality control measures were implemented for SNPs: (i) exclusion of non-biallelic SNPs and those with ambiguous allele chains; (ii) removal of SNPs lacking rs labels, duplicates, or those not included in the 1000 Genomes Project, with alleles mismatching that dataset, or located within the major histocompatibility complex region (chr6: 28.5–33.5Mb); (v) retention of SNPs with a minor allele frequency (MAF) > 0.01. Given the distinct origins of kidney-related diseases and sepsis, there was limited sample overlap between them. Finally, we preserved pertinent information for each SNP, such as effect size, standard error, effect allele, and P-value, for further analysis.
Stratified LD score regression was applied to immune cell data to study whether specific cell types harbor a greater proportion of the genetic variation (SNP heritability) underlying kidney-related traits and sepsis. We aimed to identify immune cell types with significant enrichment of SNP heritability. Data on 292 immune cell types from the ImmGen consortium (including B cells, gamma delta T cells, alpha beta T cells, innate lymphocytes, myeloid cells, stromal cells, and stem cells) were used (16). Following adjustment for baseline models and all gene sets, significance of SNP heritability enrichment estimates in each tissue and cell type was assessed using p-values from regression coefficients z-scores.
2.3.2 Analysis of pleiotropy under the composite null hypothesis
SNP-Level PLACO leverages summary-level genotype-phenotype association statistics to investigate pleiotropic loci among complex traits (17). We computed the squared Z-scores for each variant, excluding SNPs with extremely high Z2 (>80). Additionally, considering the potential correlation between kidney-related diseases and sepsis, we estimated the correlation matrix of Z-scores. Subsequently, a Horizontal α-level intersection-union test (IUT) method was employed to test the hypothesis of no pleiotropy. The final p-value from the IUT test represents the strongest evidence against the null hypothesis.To explore the shared biological mechanisms of these pleiotropic loci identified by PLACO, we further mapped them to nearby genes. MAGMA analysis (18) was conducted on genes located at or overlapping with the pleiotropic loci, utilizing both PLACO outputs and single-trait GWAS results. A significance level of P< 0.05/Ngenes = 3E-06 was applied. To identify potential functional annotations beyond the immediate risk loci, we utilized functional maps and annotations from the Functional Mapping and Annotation of Genome-Wide Association Studies (FUMA) (19). Mapping and Annotation of Genome-Wide Association Studies (FUMA) (19). Additionally, enrichment analysis of mapped genes was performed using pathways from the Molecular Signatures Database (MSigDB) (20). Additionally, we employed eQTL mapping, which identifies genes whose expression is influenced by genetic variants, to explore associations beyond the risk loci themselves. Tissue data originated from the GTEx V8 eQTL summary data for whole blood tissue (21).
2.3.3 MR analysis
We used PLINK’s clumping program (22) to identify potential causal variants among loci significantly associated with exposure (P< 5×10-8). These loci were treated as instrumental variables (IVs) in downstream analyses. All IVs underwent linkage disequilibrium (LD) clumping (r2 = 0.001; distance = 5,00 kb) to mitigate the influence of correlated SNPs. To ensure the strength of the IVs, we computed the R2 and F-statistic , for each IV, where R2 represents the proportion of variance explained by the genetic instrument, N is the effective sample size of GWAS, and k is the number of SNPs (23).
Mendelian randomization primarily relies on the Inverse Variance Weighting (IVW) approach, but this method hinges on three key assumptions. First, the genetic variants (instrumental variables) must be strongly associated with the exposure of interest. Second, these variants should not be influenced by confounding factors that affect both the exposure and the outcome. Finally, the variants’ effect on the outcome must solely act through the exposure, with no independent pathways. To assess potential violations of these assumptions, we conducted several sensitivity analyses. Firstly, both the IVW and MR-Egger methods employed the Q-test to detect heterogeneity in associations among individual instrumental variables. This test helps identify potential violations of the core assumptions. Secondly, MR-Egger was used to estimate horizontal pleiotropy based on its intercept, ensuring the genetic variants are independent influences on the exposure and outcome (24). To further strengthen the results’ stability and robustness, we performed supplementary analyses using additional MR methods with varying modeling assumptions and advantages, such as the weighted median and weighted mode.
All statistical analyses were conducted using R version 4.2.3. MR analyses were performed using the MendelianRandomization package (25).
3 Results
3.1 Genetic correlation analysis
Genetic correlation analysis revealed significant genetic correlations between sepsis and BUN (rg = 0.233, P = 4×10-4), kidney stone (rg = 0.433, P = 1.85×10-5), creatinine (rg =0.244, P = 1.85×10-7), UA (rg = 0.470, P = 7.12×10-12), and UACR (rg = 0.211, P = 5.64×10-5). To account for multiple testing, we applied Bonferroni correction (0.05/9 = 0.006). Notably, all the aforementioned associations with sepsis remained statistically significant after correction. Figure 2 illustrates the results of the genetic correlation analysis.
Figure 2 The results of the genetic correlation analysis.
3.2 Analysis of PLACO
PLACO pleiotropy analysis was performed on the five trait pairs exhibiting significant correlations (Figure 3). Details of the identified pleiotropic loci are provided in Table 1 and Supplementary Table S3. The analysis collectively identified 24 genomic regions. Notably, the 4q21.1 locus was concurrently associated with both BUN and UACR. Lead SNPs, rs10008637 and rs17319721 (both introns), reside within the SHROOM3 gene. Similarly, the 11p14.1 locus associated with BUN and Creatinine harbors lead SNPs rs685270 and rs963837 (intergenic) near genes RP5–1024C24.1 and DCDC1. The 12q13.3 locus linked to BUN and Serum urate contains lead SNPs rs540730 and rs1106766 (introns) within the R3HDM2 gene. Finally, the 13q14.3 locus associated with Creatinine and UACR harbors lead SNPs rs1327653 and rs1239707 (intergenic) with currently unidentified associated genes.
Figure 3 The Manhattan plots for pleiotropic loci in 5 pairs of results. (A) BUN and Sepsis; (B) Creatinine and Sepsis; (C) Serum urate and Sepsis; (D) UACR and Sepsis; (E) Kidney stone and Sepsis.
Table 1 Pleiotropic loci (p-value< 5E-8) identified among the 5 significant result pairs.
The QQ plot (Supplementary Figure S1) did not exhibit premature deviation from the reference line, indicating no genomic inflation. Essential information for each genomic risk locus is provided in Supplementary Figures S2, S3. Supplementary Figure S3 explores the functional consequences of cis-acting pleiotropic SNPs on genes. In the association analysis involving BUN, UACR, and kidney stones, these pleiotropic SNPs predominantly affect introns. For serum urate, they mainly impact the intergenic regions, while their influence on creatinine is relatively minor.
We further conducted MAGMA gene set enrichment analysis on the pleiotropic results using Gene Ontology terms, presenting the top 10 significantly enriched gene sets in Figure 4 and Supplementary Table S4. Functional enrichment analysis revealed that the most significant set for BUN was “GOBP_MYELOID_CELL_ACTIVATION_INVOLVED_IN_IMMUNE_RESPONSE” (p=4.95×10-6), while for creatinine, it was “GOBP_CARDIAC_VENTRICLE_DEVELOPMENT” (p = 1.42×10-5). For serum urate, the most significant was “GOBP_NEGATIVE_REGULATION_OF_ALPHA_BETA_T_CELL_DIFFERENTIATION” (p=1.25×10-6), and for UACR, the most significant was “GOBP_INNER_EAR_MORPHOGENESIS” (p=1.33×10-5). Finally, for kidney stone, the most significant enriched gene set was “GOMF_ALKALINE_PHOSPHATASE_ACTIVITY” (p=1.87×10-5). Tissue-specific MAGMA analysis revealed significant enrichment of maternal effects in the pleiotropy of the 5 pairs of results in kidney tissue. Additionally, reports of enrichment were observed in other tissues such as bladder, liver, prostate, and thyroid, among others (Supplementary Figure S4). Notably, the MAGMA gene set and tissue-specific analyses in this section were conducted using the complete distribution of SNP p-values.
Figure 4 MAGMA pathway enrichment (top 10). The blue line represents a significance level of 0.05, and the red line indicates results that remain significant after multiple corrections (P< 0.05/15485 = 3.23E-6).
We leveraged positional information of lead SNPs to identify genes associated with pleiotropic risk loci across the five significantly correlated pairs. A total of 28 genes were mapped to these loci, and details are provided in the last column of Table 1 (also available in Supplementary Table S5). We further employed MAGMA gene-set analysis to identify pleiotropic genes with genome-wide significance (Supplementary Table S6, Supplementary Figure S5). This analysis identified 33 pleiotropic genes (listed in Supplementary Table S6, Supplementary Figure S5). The QQ plot for the MAGMA gene-test is presented in Supplementary Figure S6. Additionally, we used eQTL information from whole blood tissue to explore genes potentially regulated by these pleiotropic risk loci. Specific details are provided in Supplementary Table S5.
Several pleiotropic genes displayed significant differential expression in tissues relevant to our study, such as the kidney and pancreas (Supplementary Figure S7, Supplementary Table S7). This finding is further supported by tissue-specific enrichment analysis, which revealed an enrichment of these positionally mapped genes in pancreas and kidney tissues (Figure 5, Supplementary Table S7). Pathway enrichment analysis, based on MAGMA gene testing, is illustrated in Supplementary Figure S8. The most significant enrichment was observed in “8p23.1 copy number variation syndrome”. Cell-type enrichment analysis results are displayed in Supplementary Figure S9, with the most notable enrichment found in “Aizarani liver C4 Epcam pos bile duct cells1”. Finally, the protein-protein interaction (PPI) network analysis in Supplementary Figure S10 revealed a network of six genes that interact with each other.
Figure 5 The enrichment of pleiotropic genes (based on MAGMA gene testing) across different tissues.
To further explore the immunological underpinnings of these associations, we conducted LDSC immune cell-type-specific enrichment analyses for sepsis and the kidney-related traits (Figure 6). We investigated a total of 37 immune cell types in relation to these phenotypes. Among the traits, BUN exhibited associations with 11 immune cell types, with “T.8Mem.Sp.OT1.d106.VSVOva” showing the strongest enrichment (p = 2.87×10-3). Similarly, significant immune cell enrichments were observed for Creatinine (4 identified cells, most significant: “MEChi.GFP+.Adult”, p = 1.20×10-2), Serum urate (7 identified cells, most significant: “NK.49CI+.Sp “, p = 2.05×10-2), UACR (3 identified cells, most significant: “T.DP69+.Th.v2”, p = 1.34×10-2), and Kidney stone (12 identified cells, most significant: “NK.49CI-.Sp”, p = 6.47×10-3). Interestingly, the analysis also revealed shared immune cell enrichments between certain traits. For example, “T.8Mem.Sp.OT1.d106.VSVOva, NK.MCMV7.Sp, NK.CD49b+.Lv, NK.b2m-.Sp” cells were enriched for both BUN and Kidney stone, while “NK.49CI+.Sp” and “NK.49CI-.Sp” were enriched for both Serum urate and Kidney stone. Additionally, “ILC1.CD49b-.Lv” and “Ep.5wk.MEClo.Th” cells showed enrichment for both BUN and Creatinine. Full details of all identified immune cell enrichments can be found in Supplementary Table S9.
Figure 6 The LDSC cell-type-specific heritability enrichment analysis.
3.3 MR analysis
We employed the two-sample Mendelian randomization method to investigate the causal relationship between the kidney-related diseases and sepsis. However, the IVW analysis did not identify a statistically significant causal association between them (Supplementary Table S10). Sensitivity analyses yielded results consistent with the main findings, showing no evidence of directional pleiotropy or heterogeneity (Supplementary Table S11).
4 Discussion
This study investigated the complex genetic relationship between kidney diseases and sepsis. We employed two main approaches: first, LDSC analysis to evaluate their genetic correlation, and second, a pleiotropy approach to identify shared genetic loci and genes potentially influencing both conditions. We then assessed causal associations through a two-sample bidirectional MR method. Our analyses revealed significant genetic correlations between BUN, creatinine, urine albumin, UACR, kidney stones, and sepsis. We identified 24 risk loci across five pairs of significantly correlated traits, including four loci shared by multiple traits. Further integrated analysis using MAGMA gene-test and eQTL data identified PPP2R3A, VAMP8, DOCK7, and HIBADH as shared risk genes. Genetic enrichment analysis revealed 50 potentially relevant pathways and 37 immune cell types associated with both diseases. These findings suggest a shared genetic architecture and potentially overlapping pathogenic mechanisms between kidney diseases and sepsis. This deeper understanding of pleiotropy in sepsis paves the way for future research on disease prevention.
Genetic correlation analysis revealed a significant positive genetic correlation between BUN, creatinine, UACR, serum urate, kidney stones, and sepsis. Previous studies have also reported associations among these factors, and we have summarized the findings of these studies in Supplementary Table S1 of the Supplementary Files. Our findings align with and further validate the results of these aforementioned studies. While existing reports suggest relationships between eGFR, IgA nephropathy, kidney neoplasm, diabetic nephropathy, and sepsis (26–29), our study did not identify significant connections among them. Several factors may explain these null findings: sample size limitations, population structure complexities, disease heterogeneity, and potential confounding or reverse causation inherent to observational studies could all play a role. For instance, diabetic nephropathy often results from diabetes, and the genetic correlation between diabetes and sepsis might be insignificant. Additionally, our kidney neoplasm data included cases of metastasized tumors, potentially introducing inaccuracies. The lack of association with IgA nephropathy could be due to a true lack of correlation or its complex etiology. Interestingly, while eGFR showed no significant genetic correlation with sepsis, serum creatinine levels did. This might be because eGFR is also influenced by factors like age and weight, impacting the results. Therefore, these non-significant relationships were not elaborated in the following discussion. Although MR analysis did not establish causality between kidney diseases and sepsis, our genetic correlation and pleiotropy analyses revealed a shared genetic architecture. This suggests that multiple risk loci contribute to both conditions, indicating a potential pleiotropic effect. Individuals with specific mutations for kidney diseases may also have an increased risk of sepsis, and vice versa. However, this shared genetic architecture does not imply that kidney disease and sepsis are inevitable co-occurrences. Clinicians should remain vigilant for potential sepsis in patients with persistently elevated serum creatinine, BUN, UACR or serum urate, as well as those with kidney stones. Prompt antibiotic use to prevent infections is advisable, rather than relying on advanced antibiotics after sepsis develops, which can contribute to bacterial resistance.
PLACO pleiotropy analysis of the five significant trait pairs identified 24 pleiotropic loci and 28 potentially associated genes. Notably, a shared risk locus (11p14.1) and common genes (RP5–1024C24.1, DCDC1) were found in both BUN and creatinine. While not previously linked to sepsis, this locus has been associated with sex hormone secretion (30). The gene RP5–1024C24.1 lacks current research, while DCDC1 has been implicated in esophageal cancer (31). Similarly, BUN and UACR shared a risk locus (4q21.1) and the gene SHROOM3. This locus has been previously linked to essential tremor (32), while SHROOM3 is known to play a role in kidney function and chronic kidney disease (33). Our findings suggest a potential novel genetic link between SHROOM3 and sepsis. BUN and serum urate shared another risk locus (12q13.3) and the gene R3HDM2. This previously unreported locus for BUN, serum urate, and sepsis is a novel finding for nodular disease (34). R3HDM2, linked to uric acid transport (35), aligns with our investigation and strengthens the potential genetic connection between BUN and sepsis. Creatinine and UACR shared a common risk locus (13q14.3) without identified nearby genes. While currently associated with leukemia (36), there is no existing evidence linking 13q14.3 to creatinine, UACR, or sepsis. These newly discovered loci and genes hold promise as future therapeutic targets for sepsis, warranting further investigation.
MAGMA gene-testing, whole blood tissue eQTL information, and analysis of nearby genes all converged on four genes: PPP2R3A, VAMP8, DOCK7, and HIBADH. A mouse model study suggested PPP2R3A’s role in initiating lung inflammation (37). While VAMP8 has been linked to inflammation in periodontal disease (38). DOCK7 is associated with serum lipid level (39), which are known to play a role in both sepsis and kidney health (40). Finally, HIBADH, previously connected to kidney stones (41), has been documented to be involved in valine metabolism disorders (42),a pathway recently implicated in sepsis (43). The established roles of these genes in inflammation, lipid metabolism, and metabolic disorders, alongside our findings in sepsis and kidney diseases, suggest a potential interplay in the pathogenesis of these conditions. However, further research is required to elucidate these potential connections. While tissue expression analysis revealed significantly higher levels of certain pleiotropic genes, including RAP1GAP, PBXIP1, VAMP8, SUMO1, SPP1, and SCARB2, in kidney, pancreas, and liver tissues, tissue-specific enrichment analysis further confirmed this enrichment for these positionally-mapped genes. Notably, VAMP8 emerged as a common pathogenic gene across all three methods employed. This gene exhibits high expression in the kidneys, but its precise physiological and molecular roles remain unclear (44). Interestingly, research suggests its interaction with apolipoproteins, potentially contributing to chronic kidney disease (45). Our study additionally identified its association with the pancreas, which aligns with previous reports highlighting its role in pancreatitis (46). This potentially explains the frequent clinical observation of pancreatitis progressing to sepsis and subsequent renal dysfunction. Moreover, genes like RAP1GAP, expressed in multiple tissues, have been linked to glomerular damage (47). Similarly, PBXIP1 is associated with renal fibrosis and chronic kidney disease (48), while SPP1 is related to damage in renal tubules (49). These co-expressed genes, along with VAMP8, suggest an intrinsic link between the kidneys and other organs, potentially explaining the multi-organ dysfunction observed in sepsis and the need for diverse clinical indicators for diagnosis. However, further research is crucial to fully elucidate this intricate relationship.
Several gene sets were enriched across the 5 pairs of results. T cells are categorized into αβT cells and γδT cells based on the expression of their T-cell receptors. γδT cells, known to be associated with sepsis (50), have been implicated in this process. While αβT cells are more prevalent in the kidney and linked to acute kidney injury (51), their specific role in sepsis remains unclear. Our findings suggest that serum urate, possibly through its influence on αβT cells, could play a role in sepsis development. Myeloid cells are another enriched pathway identified in our analysis. Established research confirms their involvement in sepsis (52). Consequently, elevated blood urea nitrogen levels may modulate myeloid cell function, potentially influencing susceptibility to sepsis. Interestingly, the link between creatinine and ventricular function emerged from our analysis. Studies in animal models with reduced kidney function (5/6 nephrectomy) observed left ventricular hypertrophy and progressive creatinine elevation (53). Similarly, septic cardiomyopathy with altered ventricular morphology is a common observation in sepsis patients (54). Our findings suggest a potential underlying mechanism where increased creatinine levels might impact ventricular development, influencing susceptibility to sepsis. This hypothesis warrants further investigation. Intriguingly, MAGMA analysis identified the enrichment of ion channel and transport protein pathways. These proteins are abundantly expressed in both the inner ear and the kidneys (55). A study involving Alport mice with glomerular disease and hearing loss (56) demonstrated that dual inhibition of endothelin and angiotensin receptors improved kidney and inner ear lesions. This suggests a link between the kidneys and the inner ear. In a mouse model of sepsis, observations indicated that sepsis can affect the morphology and function of the inner ear. We speculate that there might be a link connection between UACR, inner ear morphology, and sepsis. In a retrospective analysis, alkaline phosphatase was identified as a predictive factor for stone formation (57). Interestingly, research has explored its potential role in treating sepsis-related acute kidney injury (58), suggesting its potential as a novel therapeutic target for inflammation in sepsis-related kidney complications (59). Therefore, our study suggests that kidney stones might alter alkaline phosphatase, thereby influencing sepsis. At present, we have only discussed the most significant pathways in the pathway analysis among the five pairs of results. No overlapping pathways were found among the five pairs of results. Nevertheless, all enriched pathways hold promise as potential therapeutic targets for sepsis, warranting further investigation. While our pathway analysis focused on kidney-related genes, enrichment signals were also observed in other tissues such as bladder, liver, prostate, and thyroid. This suggests a broader influence beyond the kidneys, potentially involving a complex interplay between multiple organs in sepsis. The close association of BUN, creatinine, serum urate, UACR, and kidney stones with renal function strengthens the link between the kidneys and sepsis. Our findings suggest a potential genetic connection, where sepsis can trigger changes in these well-established kidney function indicators. This aligns with common observations in retrospective studies, where alterations in these markers frequently accompany sepsis events (6, 60–62). Furthermore, diseases that frequently alter these indicators may also predispose individuals to sepsis. Clinicians should therefore maintain a heightened awareness for potential sepsis development when encountering changes in these established markers and implement timely preventive measures. The insights gleaned from our results hold promise for guiding the development of effective treatments targeting these identified associations. Pathway enrichment analysis of pleiotropic genes identified by MAGMA gene testing revealed a significant association with “8p23.1 copy number variation syndrome”. Similarly, cell-type enrichment analysis highlighted “Aizarani liver C4 Epcam pos bile duct cells1”. However, no research has explored their specific connection to sepsis. While PPI network analysis identified six interacting genes, only HFE has established links to both chronic kidney disease (63) and sepsis (64). The absence of existing literature on the remaining genes should not preclude the inference of a potential genetic link between the kidney and sepsis based on their interactive network. Further research is thus necessary to elucidate the specific roles of these genes, pathways, and cell types in the context of sepsis.
The most significant cell type, “T.8Mem.Sp.OT1.d106.VSVOva”, observed in BUN and sepsis, also showed significance in kidney stone. This cell type is closely linked to CD8 T cells in the spleen. Intriguingly, CD8 T cells have been independently implicated in sepsis (65), BUN (66), and kidney stones (67), suggesting a potentially crucial role for these cells in the interplay between these conditions. Creatinine and sepsis exhibited the most significant enrichment for the immune cell type, “MEChi.GFP+.Adult”, while “ILC1.CD49b-.Lv” was co-expressed in BUN and creatinine. Current research on these two cell types is limited, and their specific relevance to sepsis remains to be elucidated. Analysis of serum urate and sepsis revealed “NK.49CI+.Sp” as the most significant cell type, which was also enriched in kidney stones. Similarly, kidney stones exhibited the most significant enrichment for “NK.49CI-.Sp”, a cell type co-expressed with serum urate. Additionally, other NK cells identified in kidney stones also showed expression in BUN. These findings collectively suggest a potential role for NK cells in the interplay between these conditions. Supporting this notion, previous studies have established the importance of NK cells in sepsis. David et al. demonstrated that NK cell count can predict mortality rates in severe sepsis (68). Furthermore, single-cell RNA sequencing studies on sepsis highlight the immunomodulatory role of NK cells (69). A connection between uric acid and NK cells has also been documented. A single-cell sequencing study of gout patients identified NK cells as potential contributors to novel therapeutic targets (70), and a mouse model demonstrated that NK cells can ameliorate gout inflammation (71). However, research on the correlation between kidney stones, BUN, and NK cells is sparse. A study on urinary tract stones with infection noted changes in NK cells, but the specific mechanisms were not explored (72). Another study investigating risk factors for death in COVID-19 patients found that elevated urea nitrogen and NK cells correlated with a poor prognosis (73). Finally, the immune cell type “T.DP69+.Th.v2” emerged as the most significant finding in the UACR and sepsis analysis. Currently, no relevant research on this specific cell type is available.
Our study has several limitations that should be acknowledged. First, the participants were predominantly of European descent. Since sepsis and kidney issues are global concerns, generalizability of these findings to other populations is crucial and requires further investigation. Second, due to limitations in acquiring individual-level data and the lack of large, well-characterized cohorts, our study relied heavily on publicly available GWAS databases. This limited the ability to perform subgroup or stratified analyses. Third, the MR analysis did not establish a causal relationship between kidney diseases and sepsis, suggesting a potential comorbid association rather than a direct causal effect. Fourth, some of the newly identified gene loci, pathways, and immune cells lack prior mention in existing literature. In the future, we aim to further explore and validate our findings through additional studies. Validation of these findings necessitates additional, well-designed experiments or clinical research. Despite these limitations, this study offers the best available evidence to date on genetic correlations between kidney diseases and sepsis.
5 Conclusion
By integrating comprehensive GWAS data and applying diverse methodologies, our study identified novel pleiotropic loci, shared pathogenic genes, enriched pathways, and immune cells. This comprehensive analysis sheds light on the complex interplay between kidney disease and sepsis, potentially revealing a shared underlying etiology. These findings offer valuable insights for early intervention and treatment strategies, paving the way for further research in this critical area.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
In our study, we utilized large-scale GWAS datasets rather than individual-level data. The studies included in these consortia obtained approval from local research ethics committees and institutional review boards, with all participants providing written informed consent.
Author contributions
TZ: Conceptualization, Writing – original draft, Writing – review & editing. YC: Data curation, Formal analysis, Writing – original draft. SJ: Data curation, Formal analysis, Writing – review & editing. LJ: Data curation, Formal analysis, Writing – review & editing. LS: Data curation, Formal analysis, Writing – review & editing. LH: Visualization, Writing – review & editing. YL: Visualization, Writing – review & editing. JY: Writing – original draft. ML: Conceptualization, Writing – original draft, Writing – review & editing.
Funding
The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.
Acknowledgments
The authors offer their thanks and appreciation to the UK Biobank, CARDIoGRAMplusC4D consortium, Köttgen A et al. GWAS, FinnGen Biobank.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2024.1396041/full#supplementary-material
References
1. Gotts JE, Matthay MA. Sepsis: pathophysiology and clinical management. Bmj. (2016) 353:i1585. doi: 10.1136/bmj.i1585
2. Cecconi M, Evans L, Levy M, Rhodes A. Sepsis and septic shock. Lancet. (2018) 392:75–87. doi: 10.1016/S0140-6736(18)30696-2
3. Rudd KE, Johnson SC, Agesa KM, Shackelford KA, Tsoi D, Kievlan DR, et al. Global, regional, and national sepsis incidence and mortality, 1990-2017: analysis for the global burden of disease study. Lancet. (2020) 395:200–11. doi: 10.1016/S0140-6736(19)32989-7
4. Wang W, Liu CF. Sepsis heterogeneity. World J Pediatr. (2023) 19:919–27. doi: 10.1007/s12519-023-00689-8
5. Bellomo R, Kellum JA, Ronco C, Wald R, Martensson J, Maiden M, et al. Acute kidney injury in sepsis. Intensive Care Med. (2017) 43:816–28. doi: 10.1007/s00134-017-4755-7
6. Weng J, Hou R, Zhou X, Xu Z, Zhou Z, Wang P, et al. Development and validation of a score to predict mortality in icu patients with sepsis: A multicenter retrospective study. J Transl Med. (2021) 19:322. doi: 10.1186/s12967-021-03005-y
7. Chuang CC, Shiesh SC, Chi CH, Tu YF, Hor LI, Shieh CC, et al. Serum total antioxidant capacity reflects severity of illness in patients with severe sepsis. Crit Care. (2006) 10:R36. doi: 10.1186/cc4826
8. Wang HE, Gamboa C, Warnock DG, Muntner P. Chronic kidney disease and risk of death from infection. Am J Nephrol. (2011) 34:330–6. doi: 10.1159/000330673
9. Dimke H, Winther-Jensen M, Allin KH, Lund L, Jess T. Risk of urolithiasis in patients with inflammatory bowel disease: A nationwide danish cohort study 1977-2018. Clin Gastroenterol Hepatol. (2021) 19:2532–40.e2. doi: 10.1016/j.cgh.2020.09.049
10. Zhang T, Cao Y, Zhao J, Yao J, Liu G. Assessing the causal effect of genetically predicted metabolites and metabolic pathways on stroke. J Transl Med. (2023) 21:822. doi: 10.1186/s12967-023-04677-4
11. Richardson TG, Leyden GM, Wang Q, Bell JA, Elsworth B, Davey Smith G, et al. Characterising metabolomic signatures of lipid-modifying therapies through drug target mendelian randomisation. PloS Biol. (2022) 20:e3001547. doi: 10.1371/journal.pbio.3001547
12. Wuttke M, Li Y, Li M, Sieber KB, Feitosa MF, Gorski M, et al. A catalog of genetic loci associated with kidney function from analyses of a million individuals. Nat Genet. (2019) 51:957–72. doi: 10.1038/s41588-019-0407-x
13. Teumer A, Li Y, Ghasemi S, Prins BP, Wuttke M, Hermle T, et al. Genome-wide association meta-analyses and fine-mapping elucidate pathways influencing albuminuria. Nat Commun. (2019) 10:4130. doi: 10.1038/s41467-019-11576-0
14. Köttgen A, Albrecht E, Teumer A, Vitart V, Krumsiek J, Hundertmark C, et al. Genome-wide association analyses identify 18 new loci associated with serum urate concentrations. Nat Genet. (2013) 45:145–54. doi: 10.1038/ng.2500
15. Kurki MI, Karjalainen J, Palta P, Sipilä TP, Kristiansson K, Donner KM, et al. Finngen provides genetic insights from a well-phenotyped isolated population. Nature. (2023) 613:508–18. doi: 10.1038/s41586-022-05473-8
17. Ray D, Chatterjee N. A powerful method for pleiotropic analysis under composite null hypothesis identifies novel shared loci between type 2 diabetes and prostate cancer. PloS Genet. (2020) 16:e1009218. doi: 10.1371/journal.pgen.1009218
18. de Leeuw CA, Mooij JM, Heskes T, Posthuma D. Magma: generalized gene-set analysis of gwas data. PloS Comput Biol. (2015) 11:e1004219. doi: 10.1371/journal.pcbi.1004219
19. Watanabe K, Taskesen E, van Bochoven A, Posthuma D. Functional mapping and annotation of genetic associations with fuma. Nat Commun. (2017) 8:1826. doi: 10.1038/s41467-017-01261-5
20. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U.S.A. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102
22. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. Plink: A tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. (2007) 81:559–75. doi: 10.1086/519795
23. Burgess S, Thompson SG. Avoiding bias from weak instruments in mendelian randomization studies. Int J Epidemiol. (2011) 40:755–64. doi: 10.1093/ije/dyr036
24. Burgess S, Thompson SG. Interpreting findings from mendelian randomization using the mr-egger method. Eur J Epidemiol. (2017) 32:377–89. doi: 10.1007/s10654-017-0255-x
25. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The mr-base platform supports systematic causal inference across the human phenome. Elife. (2018) 7:20210601. doi: 10.7554/eLife.34408
26. Xu H, Gasparini A, Ishigami J, Mzayen K, Su G, Barany P, et al. Egfr and the risk of community-acquired infections. Clin J Am Soc Nephrol. (2017) 12:1399–408. doi: 10.2215/CJN.00250117
27. Cordeiro Cabral DB, de Sandes-Freitas TV, Medina-Pestana JO, Mastroianni-Kirsztajn G. Clinical features, treatment and prognostic factors of post-transplant immunoglobulin a nephropathy. Ann Transplant. (2018) 23:166–75. doi: 10.12659/AOT.907167
28. Liu Z, Mahale P, Engels EA. Sepsis and risk of cancer among elderly adults in the United States. Clin Infect Dis. (2019) 68:717–24. doi: 10.1093/cid/ciy530
29. Jones A, Holmes J, Stephens M, Geen J, Williams J, Donovan K, et al. Using electronic aki alerts to define the epidemiology of acute kidney injury in renal transplants. J Nephrol. (2021) 34:829–38. doi: 10.1007/s40620-020-00869-z
30. Hayes MG, Urbanek M, Ehrmann DA, Armstrong LL, Lee JY, Sisk R, et al. Genome-wide association of polycystic ovary syndrome implicates alterations in gonadotropin secretion in european ancestry populations. Nat Commun. (2015) 6:7502. doi: 10.1038/ncomms8502
31. Du P, Huang P, Huang X, Li X, Feng Z, Li F, et al. Comprehensive genomic analysis of oesophageal squamous cell carcinoma reveals clinical relevance. Sci Rep. (2017) 7:15324. doi: 10.1038/s41598-017-14909-5
32. Clark LN, Gao Y, Wang GT, Hernandez N, Ashley-Koch A, Jankovic J, et al. Whole genome sequencing identifies candidate genes for familial essential tremor and reveals biological pathways implicated in essential tremor aetiology. EBioMedicine. (2022) 85:104290. doi: 10.1016/j.ebiom.2022.104290
33. Köttgen A, Glazer NL, Dehghan A, Hwang SJ, Katz R, Li M, et al. Multiple loci associated with indices of renal function and chronic kidney disease. Nat Genet. (2009) 41:712–7. doi: 10.1038/ng.377
34. Hofmann S, Fischer A, Nothnagel M, Jacobs G, Schmid B, Wittig M, et al. Genome-wide association analysis reveals 12q13.3-Q14.1 as new risk locus for sarcoidosis. Eur Respir J. (2013) 41:888–900. doi: 10.1183/09031936.00033812
35. Fernández-Torres J, Martínez-Nava GA, Oliviero F, López-Reyes AG, Martínez-Flores K, Garrido-Rodríguez D, et al. Common gene variants interactions related to uric acid transport are associated with knee osteoarthritis susceptibility. Connect Tissue Res. (2019) 60:219–29. doi: 10.1080/03008207.2018.1483359
36. Mian M, Scandurra M, Chigrinova E, Shen Y, Inghirami G, Greiner TC, et al. Clinical and molecular characterization of diffuse large B-cell lymphomas with 13q14.3 deletion. Ann Oncol. (2012) 23:729–35. doi: 10.1093/annonc/mdr289
37. Shi X, Wang J, Zhang X, Yang S, Luo W, Wang S, et al. Grem1/ppp2r3a expression in heterogeneous fibroblasts initiates pulmonary fibrosis. Cell Biosci. (2022) 12:123. doi: 10.1186/s13578-022-00860-0
38. Loos BG, Van Dyke TE. The role of inflammation and genetics in periodontal disease. Periodontol. (2020) 83:26–39. doi: 10.1111/prd.12297
39. Guo T, Yin RX, Huang F, Yao LM, Lin WX, Pan SL. Association between the dock7, pcsk9 and galnt2 gene polymorphisms and serum lipid levels. Sci Rep. (2016) 6:19079. doi: 10.1038/srep19079
40. Stasi A, Franzin R, Fiorentino M, Squiccimarro E, Castellano G, Gesualdo L. Multifaced roles of hdl in sepsis and sars-cov-2 infection: renal implications. Int J Mol Sci. (2021) 22:20210601. doi: 10.3390/ijms22115980
41. Wang L, Zhou Z, Yang Y, Gao P, Lin X, Wu Z. A genetic polymorphism in the wdr72 gene is associated with calcium nephrolithiasis in the chinese han population. Front Genet. (2022) 13:897051. doi: 10.3389/fgene.2022.897051
42. Meyer M, Hollenbeck JC, Reunert J, Seelhöfer A, Rust S, Fobker M, et al. 3-hydroxyisobutyrate dehydrogenase (Hibadh) deficiency-a novel disorder of valine metabolism. J Inherit Metab Dis. (2021) 44:1323–9. doi: 10.1002/jimd.12410
43. Xia L, Hantrakun V, Teparrukkul P, Wongsuvan G, Kaewarpai T, Dulsuk A, et al. Plasma metabolomics reveals distinct biological and diagnostic signatures for melioidosis. Am J Respir Crit Care Med. (2023) 209(3):288–98. doi: 10.1164/rccm.202207-1349OC
44. Wang CC, Ng CP, Shi H, Liew HC, Guo K, Zeng Q, et al. A role for vamp8/endobrevin in surface deployment of the water channel aquaporin 2. Mol Cell Biol. (2010) 30:333–43. doi: 10.1128/MCB.00814-09
45. Lecordier L, Heo P, Graversen JH, Hennig D, Skytthe MK, Cornet d'Elzius A, et al. Apolipoproteins L1 and L3 control mitochondrial membrane dynamics. Cell Rep. (2023) 42:113528. doi: 10.1016/j.celrep.2023.113528
46. Cosen-Binker LI, Binker MG, Wang CC, Hong W, Gaisano HY. Vamp8 is the V-snare that mediates basolateral exocytosis in a mouse model of alcoholic pancreatitis. J Clin Invest. (2008) 118:2535–51. doi: 10.1172/JCI34672
47. Potla U, Ni J, Vadaparampil J, Yang G, Leventhal JS, Campbell KN, et al. Podocyte-specific rap1gap expression contributes to focal segmental glomerulosclerosis-associated glomerular injury. J Clin Invest. (2014) 124:1757–69. doi: 10.1172/JCI67846
48. Khumukcham SS, Manavathi B. Two decades of a protooncogene hpip/pbxip1: uncovering the tale from germ cell to cancer. Biochim Biophys Acta Rev Cancer. (2021) 1876:188576. doi: 10.1016/j.bbcan.2021.188576
49. Tsai YC, Kuo MC, Huang JC, Chang WA, Wu LY, Huang YC, et al. Single-cell transcriptomic profiles in the pathophysiology within the microenvironment of early diabetic kidney disease. Cell Death Dis. (2023) 14:442. doi: 10.1038/s41419-023-05947-1
50. Zhang X, Rocha-Ferreira E, Li T, Vontell R, Jabin D, Hua S, et al. Γδt cells but not αβt cells contribute to sepsis-induced white matter injury and motor abnormalities in mice. J Neuroinflamm. (2017) 14:255. doi: 10.1186/s12974-017-1029-9
51. Newman-Rivera AM, Kurzhagen JT, Rabb H. Tcrαβ+ Cd4-/cd8- "Double negative" T cells in health and disease-implications for the kidney. Kidney Int. (2022) 102:25–37. doi: 10.1016/j.kint.2022.02.035
52. Schrijver IT, Théroude C, Roger T. Myeloid-derived suppressor cells in sepsis. Front Immunol. (2019) 10:327. doi: 10.3389/fimmu.2019.00327
53. Leineweber K, Heinroth-Hoffmann I, Pönicke K, Abraham G, Osten B, Brodde OE. Cardiac beta-adrenoceptor desensitization due to increased beta-adrenoceptor kinase activity in chronic uremia. J Am Soc Nephrol. (2002) 13:117–24. doi: 10.1681/ASN.V131117
54. Beesley SJ, Sorensen J, Walkey AJ, Tonna JE, Lanspa MJ, Hirshberg E, et al. Long-term implications of abnormal left ventricular strain during sepsis. Crit Care Med. (2021) 49:e444–e53. doi: 10.1097/CCM.0000000000004886
55. Lang F, Vallon V, Knipper M, Wangemann P. Functional significance of channels and transporters expressed in the inner ear and kidney. Am J Physiol Cell Physiol. (2007) 293:C1187–208. doi: 10.1152/ajpcell.00024.2007
56. Cosgrove D, Gratton MA, Madison J, Vosik D, Samuelson G, Meehan D, et al. Dual inhibition of the endothelin and angiotensin receptor ameliorates renal and inner ear pathologies in alport mice. J Pathol. (2023) 260:353–64. doi: 10.1002/path.6087
57. Dhayat NA, Schneider L, Popp AW, Lüthi D, Mattmann C, Vogt B, et al. Predictors of bone mineral density in kidney stone formers. Kidney Int Rep. (2022) 7:558–67. doi: 10.1016/j.ekir.2021.12.003
58. Peters E, Masereeuw R, Pickkers P. The potential of alkaline phosphatase as a treatment for sepsis-associated acute kidney injury. Nephron Clin Pract. (2014) 127:144–8. doi: 10.1159/000363256
59. Hümmeke-Oppers F, Hemelaar P, Pickkers P. Innovative drugs to target renal inflammation in sepsis: alkaline phosphatase. Front Pharmacol. (2019) 10:919. doi: 10.3389/fphar.2019.00919
60. Ren Y, Zhang L, Xu F, Han D, Zheng S, Zhang F, et al. Risk factor analysis and nomogram for predicting in-hospital mortality in icu patients with sepsis and lung infection. BMC Pulm Med. (2022) 22:17. doi: 10.1186/s12890-021-01809-8
61. Chen X, Li S, Shi C, Zhang W, Liu Z, Jiang J, et al. Risk factors and predictors of urogenous sepsis after percutaneous nephrolithotomy for idiopathic calcium oxalate nephrolithiasis. Transl Androl Urol. (2023) 12:1002–15. doi: 10.21037/tau-23-219
62. Leelahavanichkul A, Souza AC, Street JM, Hsu V, Tsuji T, Doi K, et al. Comparison of serum creatinine and serum cystatin C as biomarkers to detect sepsis-induced acute kidney injury and to predict mortality in cd-1 mice. Am J Physiol Renal Physiol. (2014) 307:F939–48. doi: 10.1152/ajprenal.00025.2013
63. Dopsaj V, Topić A, Savković M, Milinković N, Novaković I, Ćujić D, et al. Associations of common variants in hfe and tmprss6 genes with hepcidin-25 and iron status parameters in patients with end-stage renal disease. Dis Markers. (2019) 2019:4864370. doi: 10.1155/2019/4864370
64. Wizorek JJ, Turnbull IR, Buchman TG. Iron overload before cecal ligation and puncture increases mortality. Shock. (2003) 20:52–5. doi: 10.1097/01.shk.0000065770.72937.fd
65. Wesche-Soldato DE, Chung CS, Gregory SH, Salazar-Mather TP, Ayala CA, Ayala A. Cd8+ T cells promote inflammation and apoptosis in the liver after sepsis: role of fas-fasl. Am J Pathol. (2007) 171:87–96. doi: 10.2353/ajpath.2007.061099
66. Zhang S, Sun Y, Wang J, Lu Y, Yuan H, Zong Y, et al. Shuyu decoction exhibits anti-fatigue properties via alleviating exercise-induced immune dysfunction. J Ethnopharmacol. (2024) 319:117259. doi: 10.1016/j.jep.2023.117259
67. Yu R, Ge J, Lei Y. Effects of different nursing modes on immune function and renal function in patients with renal calculus undergoing percutaneous nephrolithotomy. Arch Esp Urol. (2023) 76:703–10. doi: 10.56434/j.arch.esp.urol.20237609.86
68. Andaluz-Ojeda D, Iglesias V, Bobillo F, Almansa R, Rico L, Gandía F, et al. Early natural killer cell counts in blood predict mortality in severe sepsis. Crit Care. (2011) 15:R243. doi: 10.1186/cc10501
69. Darden DB, Dong X, Brusko MA, Kelly L, Fenner B, Rincon JC, et al. A novel single cell rna-seq analysis of non-myeloid circulating cells in late sepsis. Front Immunol. (2021) 12:696536. doi: 10.3389/fimmu.2021.696536
70. Wang M, Chen W, Zhang X, Mei L, Wu X, Chen X, et al. Single-cell analysis in blood reveals distinct immune cell profiles in gouty arthritis. J Immunol. (2023) 210:745–52. doi: 10.4049/jimmunol.2200422
71. Wang J, Yang Q, Zhang Q, Yin C, Zhou L, Zhou J, et al. Invariant natural killer T cells ameliorate monosodium urate crystal-induced gouty inflammation in mice. Front Immunol. (2017) 8:1710. doi: 10.3389/fimmu.2017.01710
72. Zhao Y, Zhang Y, Liu J. Regulatory effect of pseudomonas aeruginosa mannose-sensitive hemagglutinin on inflammation and immune function in percutaneous nephrolithotomy patients with upper urinary tract calculi complicated with infection. Front Immunol. (2023) 14:1181688. doi: 10.3389/fimmu.2023.1181688
Keywords: sepsis, kidney, genetic correlation, blood urea nitrogen, creatinine, urinary albumin-to-creatinine ratio, urate, kidney stones
Citation: Zhang T, Cui Y, Jiang S, Jiang L, Song L, Huang L, Li Y, Yao J and Li M (2024) Shared genetic correlations between kidney diseases and sepsis. Front. Endocrinol. 15:1396041. doi: 10.3389/fendo.2024.1396041
Received: 05 March 2024; Accepted: 02 July 2024;
Published: 17 July 2024.
Edited by:
Hongbing Liu, Tulane University, United StatesCopyright © 2024 Zhang, Cui, Jiang, Jiang, Song, Huang, Li, Yao and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Min Li, 2504188@zju.edu.cn; Jiali Yao, yaojiali0706@163.com