Risk factors and predictors of urogenous sepsis after percutaneous nephrolithotomy for idiopathic calcium oxalate nephrolithiasis
Original Article

Risk factors and predictors of urogenous sepsis after percutaneous nephrolithotomy for idiopathic calcium oxalate nephrolithiasis

Xinfeng Chen#, Siqi Li#, Chunmei Shi, Wei Zhang, Zhenmin Liu, Jie Jiang, Yong Zhang, Zhan Chen, Bing Zheng, Hua Zhu

Department of Urology Surgery, The Second Affiliated Hospital of Nantong University, Nantong, China

Contributions: (I) Conception and design: X Chen, H Zhu; (II) Administrative support: S Li, B Zheng; (III) Provision of study materials or patients: C Shi, W Zhang; (IV) Collection and assembly of data: Z Liu, J Jiang; (V) Data analysis and interpretation: Y Zhang, Z Chen; (VI) Manuscript writing: All authors; (VII) Final approval of manuscript: All authors.

#These authors contributed to the work equally and should be considered as co-first authors.

Correspondence to: Hua Zhu, MD; Bing Zheng, MD. Department of Urology Surgery, The Second Affiliated Hospital of Nantong University, 666 Shengli Road, Nantong 226001, China. Email: ntzhuhua@126.com; ntzb2008@163.com.

Background: The aim of this study was to use bioinformatics approaches to screen and identify the key genes of idiopathic calcium oxalate nephrolithiasis, and explore its potential molecular mechanism.

Methods: The GSE73680 kidney stone data set was downloaded from the Gene Expression Omnibus (GEO). R software (The R Foundation for Statistical Computing) was used to screen differentially expressed genes. GeneMANIA and Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) databases were used to analyze related genes interacting with crucial genes, and a protein-protein interaction (PPI) network was constructed. The differential genes were then subjected to the Gene Ontology (GO) functional annotation and Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathway enrichment analysis using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) database. The clinical data of 156 patients who received percutaneous nephrolithotomy (PCNL) therapy at our facility between January 2013 and December 2017 were retrospectively analyzed. The various parameters associated with postoperative urogenous sepsis were identified using multivariable logistic regression analysis.

Results: The study discovered one differentially expressed gene was nucleotide-binding oligomerization domain-containing protein 2 (NOD2). GO and KEGG analysis showed that NOD2 might affect the occurrence of idiopathic calcium oxalate kidney stones by affecting inflammation, receptor expression, immune environment, necrosis, apoptosis, and other pathways. The clinical parameter of patients who participated in the study, including preoperative urinary white blood cell (WBC) count, preoperative urinary nitrite, stone diameter, operation time, WBC count, and WBC D values, were statistically different between the systemic inflammatory response syndrome (SIRS) group and the urosepsis group. According to multivariate logistic regression analysis, the preoperative urine nitrite, calculus diameter, blood WBC, and NOD2 expression 3 hours after surgery were all independently associated with the urosepsis development.

Conclusions: Preoperative urinary nitrite positive status, postoperative WBC count ≥2.98×109/L 3 hours after operation, stone diameter >6 cm, and low expression of NOD2 in renal papillary tissue are more likely to cause the urinary source of idiopathic calcium oxalate nephrolithiasis after PCNL urogenous sepsis. These parameters also offer a viable treatment paradigm for the perioperative management of PCNL in treating idiopathic calcium oxalate kidney stones.

Keywords: Percutaneous nephrolithotomy (PCNL); idiopathic calcium oxalate nephrolithiasis; urogenous sepsis; bioinformatics; predictors


Submitted Apr 10, 2023. Accepted for publication May 26, 2023. Published online Jun 26, 2023.

doi: 10.21037/tau-23-219


Highlight box

Key findings

• Preoperative positive urinary nitrite status, 3-hour postoperative white blood cell count ≥2.98×109/L, stone diameter >6 cm, and low expression of nucleotide-binding oligomerization domain-containing protein 2 (NOD2) in renal papilla tissue are more likely to cause urogenous idiopathic calcium oxalate kidney stones after percutaneous nephrolithotomy (PCNL) sepsis.

What is known and what is new?

• Bioinformatics analysis based on gene chip data may effectively scan and evaluate critical genes associated with the onset and progression of diseases due to the extensive use of high-throughput genomics research technologies.

• Our investigation first used bioinformatics techniques to explore the kidney stone gene chip data set in the public gene expression database and to discover new targets and insights for treating and preventing idiopathic calcium oxalate kidney stones.

What is the implication, and what should change now?

• Early judgment and intervention based on the above-mentioned independent factors will improve the effect of treatment for these patients.


Introduction

More than 1–5% of people worldwide have a history of urolithiasis (1). Calcium oxalate stones account for more than 80% of kidney stones and usually occur due to a high intake of oxalic acid or calcium forming calcium oxalate stones in the urine, which can be seen on X-ray (2). The most popular method of treatment for kidney stones is surgery. Open surgery is no longer the first option for the surgical removal of kidney stones due to the advancement and popularization of minimally invasive surgery. The preferred method of treating large, complicated upper urinary tract stones or stones that have not responded to retrograde intrarenal surgery or shock wave lithotripsy is percutaneous nephrolithotomy (PCNL) (3), the safety and efficacy of which has been improved through advancements in technology and practice. However, due to bacteria entering the patient’s bloodstream during lithotripsy, postoperative systemic inflammatory response syndrome (SIRS), with a reported incidence of ranging from 9.8% to 37%, continues to be one of the most frequent complications in PCNL (4,5). SIRS is the first step in the sepsis cascade and is closely related to the development of sepsis (6), and 0.3–4.7% of patients can develop severe life-threatening sepsis (7), while the mortality rate can reach 50–66% (8).

According to recent research, Randall’s plaque, the accumulation of calcium phosphate crystals on the renal papilla's surface, is the precursor to stone development (9). Therefore, identifying the key genes and pathways based on Randall calcium plaque is crucial for understanding how idiopathic calcium oxalate kidney stones develop. Bioinformatics analysis based on gene chip data may effectively scan and evaluate critical genes associated with the onset and progression of diseases due to the extensive use of high-throughput genomics research technologies (10). To screen the essential genes and prospective molecular mechanisms in the formation of idiopathic calcium oxalate kidney stones, our investigation first used bioinformatics techniques to explore the kidney stone gene chip data set in the public gene expression database, with the aim to identify new targets and gain insights into treating and preventing idiopathic calcium oxalate kidney stones. Then, to better inform the development of practical settings, we further investigated the association between the key genes we screened and the occurrence of postoperative urogenous sepsis. Then, based on the clinical data of patients with idiopathic calcium oxalate renal stones treated with PCNL in our hospital, we further explore the correlation between the key genes we selected and the occurrence of postoperative urinary sepsis, providing better reference opinions for clinical work. We present this article in accordance with the STREGA reporting checklist (available at https://tau.amegroups.com/article/view/10.21037/tau-23-219/rc).


Methods

Data collection and processing

The National Center for Biotechnology Information's gene expression database, the Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/) provided the kidney stone-related gene expression data set GSE73680 (11). The data set is based on the GPL17077 platform. It includes 23 cases of renal papilla tissue from individuals who had idiopathic calcium oxalate kidney stones and 6 cases with normal human renal papilla tissue. In combination with R software (The R Foundation for Statistical Computing), the platform annotation data were used to translate the gene names from the gene expression profile probe data.

Differentially expressed gene screening

The R package “limma” (https://www.r-project.org/)was used to find the differentially expressed genes in the processed GSE73680 data set, with P<0.01 and |log2 fold change (FC)|>1 as the screening conditions. With the “ggplot2” program in R, a volcano map of the differentially expressed genes was constructed.

Examining the Search Tool for the Retrieval of Interacting Genes/Proteins database and the GeneMANIA database for genes that interact with key genes

For physiological and biochemical processes, pathways, and protein–DNA, genetic, and protein-protein interactions (PPI), we used the GeneMANIA database (http://www.genemania.org) and the Search Tool for the Retrieval of Interacting Genes/Proteins (STRING) database (https://string-db.org/). Following this, we constructed a PPI network by identifying the linked genes associated with key genes from the phenotype screening, protein domain, and protein expression perspectives.

Kyoto Encyclopedia of Genes and Genomes signaling pathway enrichment and Gene Ontology functional annotation analysis

For the differentially expressed genes, analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathway enrichment and Gene Ontology (GO) functional annotation was carried out through the online biological information annotation library, Database for Annotation, Visualization, and Integrated Discovery (DAVID) 6.8, with P<0.05 as the screening threshold.

Clinical information

Statistical analysis was conducted on the clinical data of 421 patients who underwent initial PCNL therapy at the Second Affiliated Hospital of Nantong University, between January 2013 and December 2017. Patients with end-stage renal disease, tumors, blood diseases, immunosuppressive disorders, horseshoe kidneys, and preoperative heartbeats of >90 beats per minute were excluded from the study. Only 135 patients experienced SIRS, and the other 21 patients developed urosepsis. As per the 2015 urinary system infection recommendations set by the European Association of Urology (EAU), we defined SIRS as urinary sepsis with severe urinary tract infection that causes pathogenic microorganisms in urine to enter the blood system through blood vessels, causing bacteremia, sepsis, and other adverse conditions (12). The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The study was approved by ethics board of the Second Affiliated Hospital of Nantong University (No. 2021KT158) and informed consent was taken from all the patients.

Preoperative examination and treatment

Preoperative intravenous urography (IVU), CT, ultrasonography, and kidney, ureter, and bladder (KUB) X-ray were all routinely conducted. According to standard procedure, patients with positive urine cultures receive the appropriate antibiotic treatment for 3–7 days based on culture antibiotic test findings, but urine cultures were not be repeated. Then, 30 minutes before surgery, we administered broad-spectrum prophylactic antibiotics.

Recorded metrics

Clinical characteristics including postoperative white blood cell (WBC) count, stone composition, operation time, stone location, stone size, nitrite status, urinary WBC in preoperative urinalysis, gender, and age were examined. The postoperative WBC was assessed immediately following the operation or within 3 hours. Operative time was measured from the start of the renal biopsy to the end of the nephrostomy tube insertion. Meanwhile, the components of the stones were analyzed. However, the composition of many stones can be mixed, and so we recorded the main components. Meanwhile, renal papilla tissue biopsy was performed during the operation, and the expression of the key genes we screened was judged using the immunohistochemical method.

Treatment

Participants were seated in the prone position for all surgeries while under general anesthesia. After establishing the best possible renal access, the opening was expanded with an 8–18 F fascial dilator and a corresponding Amplatz sheath. Then, we applied a pneumatic lithotripter, an 8/9.8 F semirigid ureteroscope, or Hol: YAG laser lithotripsy. Stone tweezers are used to remove larger fragments, and a pulse perfusion pump was used to flush out tiny debris. After the procedure ended, a silicone rubber nephrostomy tube and a double-J stent were inserted.

Statistical analysis

The statistical evaluation was conducted using SPSS software (v.20.0, IBM, Armonk, NY, USA). Continuous variables were compared using the independent samples t test and are presented as the mean and ± SD. The χ2 test was used to compare categorical variables. To identify independent characteristics linked with the occurrence of postoperative urosepsis, multivariate logistic regression analysis was employed. The aforesaid online database automatically conducted the statistical analyses of the subsequent investigations. An acceptable P value or log-rank P value of <0.05 was used to determine statistical significance.


Results

The differentially expressed genes between the patients with idiopathic calcium oxalate nephrolithiasis and normal human renal papilla tissue

The GSE73680 data set selected in this study contains 6 cases of normal human renal papilla tissue and 26 cases of patients with calcium oxalate stones. With P<0.01 and |log2FC|>1 as screening conditions, R software demonstrated that 7,799 genes were differentially expressed in the renal papilla of participants with idiopathic calcium oxalate kidney stones compared to normal people. Moreover, 2,282 genes were upregulated, and 1,930 were downregulated (see the volcano plot, principal component analysis plot, heat map, and sample-normalized box plot in Figure 1A). Compared with those in healthy people, downregulated genes may have more clinical application feasibility and more research value, so we only selected genes with downregulated expression among common differential genes for further analysis. Considering that the pyroptosis mechanism is one of the research hotspots in the current research on the mechanism of kidney stones, we took the intersection of the downregulated differentially expressed genes and all pyroptotic genes to construct a Venn diagram (Figure 1B). From this, we finally obtained 1 differentially expressed gene, nucleotide-binding oligomerization domain-containing protein 2 (NOD2). The definitions of pyroptosis-related genes vary significantly across the different articles, and thus we selected the 33 most widely used genes from the related literature (13-16), the detailed summary of which is shown in Table 1.

Figure 1 Genes that were differently expressed in idiopathic calcium oxalate nephrolithiasis were screened. (A) The volcano map, principal component analysis map, sample-normalized box plot, and the heat map of differentially expressed genes in the GSE73680 data set. Blue points indicate downregulation, red points indicate upregulation, and black points represent |log FC| <1. (B) Venn diagram of the downregulated genes and all pyroptotic genes in the GSE73680 data set. FC, fold change.

Table 1

All pyroptosis-related genes

Gene Common name
AIM2 Absent in melanoma2
CASP1 Cysteine-aspartic acid protease-1
CASP3 Cysteine-aspartic acid protease-3
CASP4 Cysteine-aspartic acid protease-4
CASP5 Cysteine-aspartic acid protease-5
CASP6 Cysteine-aspartic acid protease-6
CASP8 Cysteine-aspartic acid protease-8
CASP9 Cysteine-aspartic acid protease-9
ELANE Elastase, neutrophil expressed
GPX4 Glutathione peroxidase-4
GSDMA Gasdermin A
GSDMB Gasdermin B
GSDMC Gasdermin C
GSDMD Gasdermin D
GSDME Gasdermin E
IL-18 Interleukin-18
IL-1β Interleukin-1β
IL-6 Interleukin-6
NLRC4 NLR family CARD domain containing 4
NLRP1 NLR family pyrin domain containing 1
NLRP2 NLR family pyrin domain containing 2
NLRP3 NLR family pyrin domain containing 3
NLRP6 NLR family pyrin domain containing 6
NLRP7 NLR family pyrin domain containing 7
NOD1 Nucleotide-binding oligomerization domain containing 1
NOD2 Nucleotide-binding oligomerization domain containing2
PJVK Pejvakin/deafness, autosomal recessive 59
PLCG1 Phospholipase C gamma1
PRKACA Protein kinase camp-activated catalytic subunit alpha
PYCARD PYD and CARD domains containing
SCAF11 SR-related CTD-associated factor 11
TIRAP TIR domain-containing adaptor protein
TNF Tumor necrosis factor

Identification of NOD2 interacting genes and proteins

Using GeneMania, we created a gene-gene interaction network of NOD2 and adjusted the neighboring genes (Figure 2A). The network of NOD2 PPIs was established using the STRING database (Figure 2B).

Figure 2 Identification of the NOD2 interacting genes and proteins. (A) Gene-gene interaction network of NOD2 in the GeneMania database. (B) PPI network of NOD2 in the STRING database. NOD2, nucleotide binding oligomerization domain 2; PPI, protein-protein interaction.

KEGG and GO enrichment analysis

Using the DAVID online database, the enrichment analysis of genes interacting with NOD2 in the PPI network diagram was carried out with the background of Homo sapiens, and enrichment information on GO was obtained. The term “biological process” refers to a group of processes that primarily controls the NLRP3 inflammasome complex, necrotic cell death, cytokine-mediated signaling pathways, cytoplasmic pattern recognition receptor signaling pathways, and I-κB kinase/nuclear factor (NF)-κB signaling, among others (Figure 3A and Table 2). For the term “cellular composition”, interacting genes were mainly enriched in the serine/threonine protein kinase complex, plasma membrane receptor complex, CD40 receptor complex, phagophore, inflammasome complex, and ubiquitin junction enzyme complexes, among others (Figure 3B and Table 3). For the term, "molecular function", the interacting genes were mainly linked to the caspase recruitment domain, ubiquitin-binding, cysteine-type endopeptidase, enzyme regulation activity, apoptosis process, protein serine/threonine amino acid kinase activity, cytokine receptor binding, and protein tyrosine kinase activity (Figure 3C and Table 4). The enrichment analysis of KEGG signaling pathways revealed that the interacting genes were mainly related to T helper 1 (Th1) and Th2 cell differentiation, PI3K-Akt signaling pathway, programmed death-ligand 1 (PD-L1) expression in cancer and programmed cell death protein 1 (PD-1) checkpoint pathway, the NOD (Nucleotide Binding Oligomerization Domain Containing)-like receptor signaling pathway, MAPK(mitogen-activated protein kinase) signaling pathway, tumor necrosis factor (TNF) signaling pathway, and T-cell receptor signaling pathway, among others (Figure 3D and Table 5).

Figure 3 The outcomes of the GO and KEGG enrichment analysis. (A) Results of BP enrichment, (B) CC enrichment, (C) MF enrichment, and (D) KEGG signaling pathway enrichment analysis. GO, Gene Ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; BP, biological process; CC, cell component; MF, molecular function.

Table 2

Results of GO biological process analysis

ID Description Gene ratio Bg-Ratio P value Adjusted P value Q value
GO:0043123 Positive regulation of I-κb kinase/NF-κb signaling 23/49 183/18670 7.11e-34 7.98e-31 4.09e-31
GO:0002753 Cytoplasmic pattern recognition receptor signaling pathway 17/49 62/18670 3.93e-31 8.51e-29 4.36e-29
GO:0002758 Innate immune response-activating signal transduction 24/49 298/18670 1.30e-30 2.50e-28 1.28e-28
GO:0031098 Stress-activated protein kinase signaling cascade 22/49 315/18670 1.57e-26 2.26e-24 1.16e-24
GO:0051090 DNA-binding transcription factor activity regulation 22/49 432/18670 1.71e-23 1.84e-21 9.46e-22
GO:0001959 Cytokine-mediated signaling pathway regulation 16/49 177/18670 5.45e-21 4.29e-19 2.20e-19
GO:0010803 TNF-mediated signaling pathway regulation 12/49 61/18670 3.94e-20 2.72e-18 1.40e-18
GO:0031663 Lipopolysaccharide-mediated signaling pathway 8/49 58/18670 2.14e-12 4.95e-11 2.54e-11
GO:0051222 Protein transport positive regulation 13/49 418/18670 3.73e-11 6.72e-10 3.45e-10
GO:0032642 Chemokine production regulation 8/49 82/18670 3.80e-11 6.78e-10 3.48e-10
GO:0038095 Fc-epsilon receptor signaling pathway 8/49 169/18670 1.26e-08 1.43e-07 7.33e-08
GO:0006914 Autophagy 11/49 496/18670 4.87e-08 4.87e-07 2.50e-07
GO:0042088 Immune response of T-helper type 1 5/49 43/18670 9.02e-08 8.67e-07 4.44e-07
GO:0060544 Necroptotic process regulation 4/49 26/18670 6.00e-07 5.14e-06 2.64e-06
GO:0010939 Necrotic cell death regulation 4/49 38/18670 2.89e-06 2.21e-05 1.13e-05
GO:1900225 NLRP3 inflammasome complex assembly regulation 2/49 13/18670 5.17e-04 0.002 0.001

P<0.05 indicates a statistically significant difference. GO, Gene Ontology.

Table 3

Results of GO cell component analysis

ID Description Gene ratio Bg-Ratio P value Adjusted P value Q value
GO:0035631     CD40 receptor complex 5/49 11/19717 3.51e-11 4.70e-09 3.18e-09
GO:0045335     Phagocytic vesicle 7/49 132/19717 3.49e-08 1.81e-06 1.22e-06
GO:0098802     Plasma membrane receptor complex 9/49 295/19717 4.04e-08 1.81e-06 1.22e-06
GO:0030139     Endocytic vesicle 8/49 303/19717 7.40e-07 2.48e-05 1.67e-05
GO:0030670     Phagocytic vesicle membrane 5/49 76/19717 1.24e-06 3.33e-05 2.25e-05
GO:0005776     Autophagosome 5/49 93/19717 3.39e-06 6.49e-05 4.38e-05
GO:0098562     Cytoplasmic side of membrane 6/49 178/19717 5.04e-06 8.45e-05 5.71e-05
GO:0098589     Membrane region 7/49 328/19717 1.56e-05 1.74e-04 1.18e-04
GO:0044445     Cytosolic part 6/49 247/19717 3.24e-05 3.34e-04 2.26e-04
GO:0001891     Phagocytic cup 3/49 27/19717 4.05e-05 3.87e-04 2.62e-04
GO:0030666     Endocytic vesicle membrane 5/49 167/19717 5.79e-05 5.17e-04 3.50e-04
GO:1902554     Serine/threonine protein kinase complex 4/49 88/19717 6.73e-05 5.36e-04 3.62e-04
GO:0000407     Phagophore assembly site 3/49 32/19717 6.80e-05 5.36e-04 3.62e-04
GO:1902911     Protein kinase complex 4/49 109/19717 1.55e-04 0.001 7.77e-04
GO:0010008     Endosome membrane 7/49 479/19717 1.70e-04 0.001 8.09e-04
GO:0061702     Inflammasome complex 2/49 14/19717 5.40e-04 0.004 0.002
GO:0000421     Autophagosome membrane 2/49 36/19717 0.004 0.023 0.016
GO:0061695     Transferase complex, transferring phosphorus-containing groups 4/49 259/19717 0.004 0.024 0.016
GO:0000151     Ubiquitin ligase complex 4/49 282/19717 0.005 0.031 0.021

P<0.05 indicates a statistically significant difference. GO, Gene Ontology.

Table 4

Results of GO molecular function analysis

ID Description Gene ratio Bg-Ratio P value Adjusted P value Q value
GO:0050700    CARD domain binding 7/49 16/17697 8.95e-15 1.34e-12 7.63e-13
GO:0043130    Ubiquitin binding 6/49 76/17697 6.20e-08 3.90e-06 2.22e-06
GO:0004708    MAP kinase kinase activity 4/49 16/17697 9.21e-08 3.90e-06 2.22e-06
GO:0043028    Cysteine-type endopeptidase regulator activity involved in apoptotic process 5/49 42/17697 1.04e-07 3.90e-06 2.22e-06
GO:0032182    Ubiquitin-like protein binding 6/49 96/17697 2.52e-07 7.56e-06 4.30e-06
GO:0097110    Scaffold protein binding 5/49 59/17697 5.90e-07 1.48e-05 8.39e-06
GO:0005126    Cytokine receptor binding 8/49 286/17697 1.07e-06 2.30e-05 1.30e-05
GO:0019787    Ubiquitin-like protein transferase activity 9/49 407/17697 1.50e-06 2.81e-05 1.60e-05
GO:0004674    Protein serine/threonine kinase activity 9/49 439/17697 2.79e-06 4.65e-05 2.65e-05
GO:0004842    Ubiquitin-protein transferase activity 8/49 382/17697 9.11e-06 1.24e-04 7.06e-05
GO:0008656    Cysteine-type endopeptidase activator activity involved in apoptotic process 3/49 16/17697 1.09e-05 1.26e-04 7.14e-05
GO:0042834    Peptidoglycan binding 3/49 17/17697 1.32e-05 1.32e-04 7.50e-05
GO:0016505    Peptidase activator activity involved in apoptotic process 3/49 19/17697 1.87e-05 1.76e-04 9.98e-05
GO:0004713    Protein tyrosine kinase activity 5/49 134/17697 3.37e-05 2.41e-04 1.37e-04
GO:0005164    Tumor necrosis factor receptor binding 3/49 31/17697 8.49e-05 5.79e-04 3.29e-04

P<0.05 indicates a statistically significant difference. GO, Gene Ontology.

Table 5

Results of Kyoto Encyclopedia of Genes and Genomes analysis

ID Description Gene ratio Bg-Ratio P value Adjusted P value Q value
hsa04621     NOD-like receptor signaling pathway 28/47 181/8076 3.69e-35 4.51e-33 2.10e-33
hsa04668     TNF signaling pathway 16/47 112/8076 6.54e-19 1.99e-17 9.29e-18
hsa04064     NF-kappa B signaling pathway 15/47 104/8076 8.40e-18 2.05e-16 9.54e-17
hsa05235     PD-L1 expression and PD-1 checkpoint pathway in cancer 9/47 89/8076 1.54e-09 9.37e-09 4.37e-09
hsa04010     MAPK signaling pathway 13/47 294/8076 7.00e-09 3.71e-08 1.73e-08
hsa04623     Cytosolic DNA-sensing pathway 7/47 63/8076 6.15e-08 2.89e-07 1.34e-07
hsa04625     C-type lectin receptor signaling pathway 8/47 104/8076 1.20e-07 5.21e-07 2.43e-07
hsa04660     T-cell receptor signaling pathway 6/47 104/8076 2.76e-05 9.62e-05 4.48e-05
hsa04659     Th17 cell differentiation 6/47 107/8076 3.25e-05 1.10e-04 5.12e-05
hsa04662     B-cell receptor signaling pathway 3/47 82/8076 0.012 0.025 0.012
hsa04151     PI3K-Akt signaling pathway 6/47 354/8076 0.016 0.033 0.015
hsa04658     Th1 and Th2 cell differentiation 3/47 92/8076 0.016 0.033 0.015

P<0.05 indicates a statistically significant difference.

General information about the enrolled patients

This retrospective study ultimately included 156 patients treated for PCNL. Of these, 135 patients only developed SIRS, while 21 developed urosepsis postoperatively. Preoperative urinary nitrite, calculus diameter, operating duration, WBC count, and WBC D values all showed statistically significant differences between the 2 groups (P<0.05; Table 6).

Table 6

General information of the enrolled patients

Variables SIRS group Urosepsis group P value
Gender
   Male 62 12 0.338
   Female 73 9
Age (years)
   <60 109 14 0.142
   ≥60 26 7
Preoperative white blood cells in urine
   Negative 54 3 0.023
   Positive 81 18
Preoperative urinary nitrite
   Negative 91 2 <0.001
   Positive 44 19
Preoperative urine culture
   Negative 86 11 0.320
   Positive 49 10
History of hypertension
   No 111 18 0.694
   Yes 24 3
History of diabetes
   No 120 20 0.372
   Yes 15 1
Stone location
   Left 62 9 0.793
   Right 73 12
Stone diameter (cm)
   <6 73 5 0.010
   ≥6 62 16
Time of operation (min)
   <120 110 9 <0.001
   ≥120 25 12
3-h postoperative WBC count (×109/L)
   <2.98 10 19 <0.001
   ≥2.98 125 2
WBC D value (×109/L) 3.19 ± 5.47 −3.77 ± 1.59 <0.001

WBC, white blood cell. D value of WBC = postoperative white blood cell count – preoperative white blood cell count; The results of WBC D value were expressed as (x¯±s). SIRS, systemic inflammatory response syndrome.

Univariate and multivariate analysis of urogenous sepsis risk variables following PCNL

In renal papillary tissue samples, we assessed the expression of NOD2 using immunohistochemistry. The immunohistochemical data were graded on the basis of the degree of staining in the cells and the percentage of positive cells. They were scored according to staining intensity and the percentage of positively stained cells. Staining intensity was scored as follows: 0 = negative or no staining, 1 = weakly positive, 2 = moderately positive, and 3 = strongly positive. Meanwhile, the percentage of staining was scored as follows: 0= negative, 1=≤25%, 2=25–50%, 3=51–75%, and 4≥75%. The result of multiplying the 2 scores represented the final score for each specimen. The samples with a score lower than 7 were ultimately classified as NOD2 low expression after the arithmetic mean of these scores was computed.

The Cox regression analysis was conducted on the factors that affected the occurrence of urogenous sepsis after PCNL. It was ultimately found that preoperative urinary nitrite, stone diameter, 3-hour postoperative blood WBC count, and NOD2 expression were independent prognostic factors that affected the occurrence of urogenous sepsis. Preoperative positive urinary nitrite, diameter, 3-hour postoperative WBC ≥2.98×109/L, stone diameter >6 cm, and low expression of NOD2 in renal papilla tissue were more likely to cause urogenous idiopathic calcium oxalate kidney stones after PCNL sepsis (Table 7).

Table 7

The univariate and multivariate analyses of variables associated with urogenous sepsis following PCNL

Variables Univariate analysis Multivariate analysis
HR 95% CI P value     HR   95% CI P value
Gender 1.052 0.104–3.751 0.554
Age 2.517 0.933–6.172 0.275
Preoperative urine WBC 5.325 1.034–16.271 0.013 2.751 0.769-6.488 0.212
Preoperative urinary nitrite 8.157 3.071–20.613 0.000 4.571   1.049–13.854 0.031
Preoperative urine culture 1.997 0.532–4.066 0.402
History of high blood pressure 0.771 0.210–2.827 0.695
Diabetes history 1.021 0.057–3.338 0.569
Stone composition 1.354 0.345–3.712 0.441
Stone location 1.008 0.097–3.086 0.578
Operation time 2.854 1.030–5.732 0.233
Diameter of stone 9.078 3.306–23.372 0.000 5.125   1.825–13.516 0.020
3-h postoperative blood WBC 7.592 2.901–16.041 0.000 4.351   1.012–12.459 0.039
Expression of NOD2 in renal papillary tissue 10.351 1.831–9.886 0.000 4.593   1.246-10.822 0.025

HR, hazard ratio; CI, confidence interval; NOD2, nucleotide binding oligomerization domain 2; WBC, white blood cell.


Discussion

Kidney stone development is a complicated process. According to the traditional view, kidney stone formation involves several biochemical functions, such as crystal nucleation, urine supersaturation, crystal growth, and aggregation (17). Earlier research presented several models to describe how kidney stones arise, such as vascular calcification and medullary–papillary complex mineralization (18). However, there are certain limitations to these models. The renal papillary Randall calcium spot model is currently a relatively accepted theory. In this model, supersaturated calcium salt crystals are deposited on the top of the renal medulla renal papillae, and the crystals gradually grow and aggregate, eventually forming urinary calculi. With the advancement of laparoscopic technology, research has revealed that Randall's calcium plaques, which frequently have stones attached, are present in the renal papilla of most patients with calcium-containing stones (19). Renal papillary biopsy studies also suggest that Randall's calcium plaques are closely related to urinary calcium levels, calcium oxalate stones, and kidney stones secondary to some metabolic diseases (20). Therefore, bioinformatics methods to study the key genes and signaling pathways of kidney stone formation based on Randall calcium plaques of renal papillae are of great value.

One study showed that urogenous sepsis is a dangerous complication after PCNL (21). The occurrence of urogenous sepsis not only aggravates the primary disease and reduces the effect of PCNL surgery but also prolongs the treatment time and increases the treatment cost and the risk of death, which has become one of the issues of high clinical concern. According to one survey, every 1-hour delay in treating urogenous sepsis increases the risk of death by 10% (22). However, the early clinical symptoms of most patients with urogenous sepsis lack specificity and only manifest as arrhythmia and body temperature changes, making it difficult to diagnose the disease accurately. Therefore, developing an early predictive index with high specificity and sensitivity to diagnose the condition of patients with urosepsis is an urgent clinical need. In order to determine the final factors that influence the development of urogenous sepsis after PCNL and to identify specific clinical values, our study combined the analysis of key genes that promote the occurrence of idiopathic calcium oxalate kidney stones with the clinical data analysis of our hospital.

In this work, the differentially expressed genes for 2,282 upregulated and 1,930 downregulated genes were examined in the renal papilla gene expression matrix of patients with idiopathic calcium oxalate nephrolithiasis, which was retrieved from the GEO database. A variety of factors are known to influence the formation and occurrence of kidney stones: for example, excessive autophagy mediated by endoplasmic reticulum stress effectively protects kidney function and impedes apoptosis and kidney stone formation (23), while protein-mediated renal oxidative damage and inflammation aggravate calcium oxalate stone formation (24); these mechanisms mainly involve inflammation and other factors. Pyroptosis is a newly recognized type of planned cell death in which cells grow until the cell membrane ruptures, spilling the cell's contents and inciting robust inflammatory and immunological responses. In-depth research on pyroptosis, which is closely involved in the development and occurrence of diseases will help us better comprehend its function in the incidence, progression, and outcome of associated diseases and elicit new ideas for clinical treatment and prevention (25,26). Moreover, pyroptosis and inflammation are also complementary in many diseases. Indeed, the mechanism of pyroptosis has emerged as one of the research hotspots in the current research on the mechanism of kidney stones. Recently, there have also been relevant studies, such as miRNA-141-3p mediated by NLRP3 in renal tubular epithelial cell injury induced by calcium oxalate (27); Long chain non coding RNA 00339 (LINC00339) promotes renal tubular epithelial cell apoptosis by regulating the miR-22-3p/NLRP3 axis in calcium oxalate induced renal stones (28). Therefore, we finally constructed a Venn diagram by intersecting the downregulated differentially expressed genes with all known pyroptotic genes and eventually obtained a single differentially expressed gene, NOD2.

Then, we constructed the NOD2-interacting gene network using the GeneMania and STRING databases and performed GO and KEGG analysis on these interacting genes. The results showed that NOD2 might affect inflammation, receptor expression, immune environment, necrosis, apoptosis, and other processes to influence the occurrence of idiopathic calcium oxalate kidney stones.

Next, we retrospectively analyzed 156 patients treated for PCNL in our hospital. Of these, 135 patients only developed SIRS, while 21 developed urosepsis postoperatively. There were statistically significant variations between the 2 groups in terms of preoperative urine WBC count, preoperative urinary nitrite, stone diameter, operation duration, WBC count, and WBC D values. According to Ferry et al. (29), urine sediment microscopy and the WBC nitrite test demonstrated the best diagnostic accuracy for lower urinary tract infections. According to their study, nitrite could improve the diagnostic efficacy for urinary tract infections, which is in agreement with our findings. Preoperative urine nitrite testing revealed that the urosepsis group had a higher positive rate (90.5% vs. 32.6%; P<0.001) than did the SIRS group, and is associated with urinary sepsis as an independent risk factor

Cox regression analysis was conducted on the factors that affect the occurrence of urogenous sepsis after PCNL which ultimately revealed that preoperative urinary nitrite, stone diameter, 3-hour postoperative blood WBC count, and NOD2 expression are the independent prognostic factors that affect the occurrence of urogenous sepsis. WBCs are the primary line of protection against external bacterial infection. In the early stages of urogenic sepsis, WBC will significantly decrease. This is due to the release of endotoxin into the bloodstream and the subsequent production of inflammatory factors that activate the systemic inflammatory response. Typically, this reaction occurs intraoperatively or within the first 2 hours after surgery (30). The WBC count then increases and peaks on the first postoperative day. WBC counts that do not increase are indicative of a poor prognosis. These observations are consistent with our findings, specifically, that urosepsis was independently linked with individuals whose postoperative WBC count was less than 2.98×109/L. A study by Wu et al. (31) reported that a WBC count of 2.85×109/L or below within 2 hours after surgery might predict urinary tract syndrome following endoscopic lithotripsy of the upper urinary tract. Therefore, our research results have certain accuracy and credibility. Although immunohistochemical detection of NOD2 expression is an invasive procedure, this study is more convenient and cost-effective compared with the cost of testing such as single-cell sequencing, so it has good clinical value to help early detection and intervention of urogenic sepsis after PCNL.

Our study does have several limitations, including a limited sample size and its retrospective design. In the future, we will conduct relevant prospective and/or randomized studies to further evaluate the accuracy of these results and provide more valuable information for clinical management.


Conclusions

If the patient has the following preoperative urinary nitrite positivity and blood white blood cells ≥2.98 at 3 hours after surgery × Idiopathic calcium oxalate kidney stones with a stone diameter greater than 6cm and low expression of NOD2 in renal papilla tissue at 109/L are more prone to urinary sepsis after PCNL. The first 6 hours following a medical assessment of urosepsis are regarded as the “golden hour” for interventional therapy (32). Therefore, early judgment and intervention based on the above-mentioned independent factors will improve the effect of treatment for these patients.


Acknowledgments

Funding: This study was supported by the Nantong Science and Education Strong Health project (No. Key 57, Young 68), the Research Project of Nantong Municipal Health Commission (No. MB2021025), the Nantong Science and Technology Project (No. MS22022085), the Nantong Science and Technology Planning Project (No. MSZ2022101), the Jiangsu Provincial Cadre Health Research Project (No. BJ21010), and the Science and Technology Project of Jiangsu Provincial Health Commission (No. LKM2022059).


Footnote

Reporting Checklist: The authors have completed the STREGA reporting checklist. Available at https://tau.amegroups.com/article/view/10.21037/tau-23-219/rc

Data Sharing Statement: Available at https://tau.amegroups.com/article/view/10.21037/tau-23-219/dss

Peer Review File: Available at https://tau.amegroups.com/article/view/10.21037/tau-23-219/prf

Conflicts of Interest: All authors have completed the ICMJE uniform disclosure form (available at https://tau.amegroups.com/article/view/10.21037/tau-23-219/coif). The authors have no conflicts of interest to declare.

Ethical Statement: The authors are accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved. The study was conducted in accordance with the Declaration of Helsinki (as revised in 2013). The study was approved by ethics board of the Second Affiliated Hospital of Nantong University (No. 2021KT158) and informed consent was taken from all the patients.

Open Access Statement: This is an Open Access article distributed in accordance with the Creative Commons Attribution-NonCommercial-NoDerivs 4.0 International License (CC BY-NC-ND 4.0), which permits the non-commercial replication and distribution of the article with the strict proviso that no changes or edits are made and the original work is properly cited (including links to both the formal publication through the relevant DOI and the license). See: https://creativecommons.org/licenses/by-nc-nd/4.0/.


References

  1. Vasudevan V, Samson P, Smith AD, et al. The genetic framework for development of nephrolithiasis. Asian J Urol 2017;4:18-26. [Crossref] [PubMed]
  2. Wang C, Dong X, Yin X, et al. Impact of intestinal flora on calcium oxalate stones. Zhong Nan Da Xue Xue Bao Yi Xue Ban 2021;46:1285-9. [PubMed]
  3. Palazzo S, Colamonico O, Forte S, et al. Experience of percutaneous access under ultrasound guidance in renal transplant patients with allograft lithiasis. Arch Ital Urol Androl 2016;88:337-40. [Crossref] [PubMed]
  4. El-Nahas AR, Shokeir AA, El-Assmy AM, et al. Post-percutaneous nephrolithotomy extensive hemorrhage: a study of risk factors. J Urol 2007;177:576-9. [Crossref] [PubMed]
  5. Singh P, Yadav S, Singh A, et al. Systemic Inflammatory Response Syndrome Following Percutaneous Nephrolithotomy: Assessment of Risk Factors and Their Impact on Patient Outcomes. Urol Int 2016;96:207-11. [Crossref] [PubMed]
  6. Wu W, Zhang D, Jin T, et al. Progress in the study of biomarkers for early prediction of systemic inflammatory response syndrome after percutaneous nephrolithotomy. Front Immunol 2023;14:1142346. [Crossref] [PubMed]
  7. Michel MS, Trojan L, Rassweiler JJ. Complications in percutaneous nephrolithotomy. Eur Urol 2007;51:899-906; discussion 906. [Crossref] [PubMed]
  8. Wagenlehner FM, Pilatz A, Naber KG, et al. Therapeutic challenges of urosepsis. Eur J Clin Invest 2008;38:45-9. [Crossref] [PubMed]
  9. Verrier C, Bazin D, Huguet L, et al. Topography, Composition and Structure of Incipient Randall Plaque at the Nanoscale Level. J Urol 2016;196:1566-74. [Crossref] [PubMed]
  10. Sedlazeck FJ, Lee H, Darby CA, et al. Piercing the dark matter: bioinformatics of long-range sequencing and mapping. Nat Rev Genet 2018;19:329-46. [Crossref] [PubMed]
  11. Taguchi K, Hamamoto S, Okada A, et al. Genome-Wide Gene Expression Profiling of Randall's Plaques in Calcium Oxalate Stone Formers. J Am Soc Nephrol 2017;28:333-47. [Crossref] [PubMed]
  12. Vaggers S, Puri P, Wagenlehner F, et al. A Content Analysis of Mobile Phone Applications for the Diagnosis, Treatment, and Prevention of Urinary Tract Infections, and Their Compliance with European Association of Urology Guidelines on Urological Infections. Eur Urol Focus 2021;7:198-204. [Crossref] [PubMed]
  13. Karki R, Kanneganti TD. Diverging inflammasome signals in tumorigenesis and potential targeting. Nat Rev Cancer 2019;19:197-214. [Crossref] [PubMed]
  14. Xia X, Wang X, Cheng Z, et al. The role of pyroptosis in cancer: pro-cancer or pro-"host"? Cell Death Dis 2019;10:650. [Crossref] [PubMed]
  15. Wang B, Yin Q. AIM2 inflammasome activation and regulation: A structural perspective. J Struct Biol 2017;200:279-82. [Crossref] [PubMed]
  16. Man SM, Kanneganti TD. Regulation of inflammasome activation. Immunol Rev 2015;265:6-21. [Crossref] [PubMed]
  17. Khan SR, Pearle MS, Robertson WG, et al. Kidney stones. Nat Rev Dis Primers 2016;2:16008. [Crossref] [PubMed]
  18. Hsi RS, Ramaswamy K, Ho SP, et al. The origins of urinary stone disease: upstream mineral formations initiate downstream Randall's plaque. BJU Int 2017;119:177-84. [Crossref] [PubMed]
  19. Pless MS, Williams JC Jr, Andreassen KH, et al. Endoscopic observations as a tool to define underlying pathology in kidney stone formers. World J Urol 2019;37:2207-15. [Crossref] [PubMed]
  20. Marien TP, Miller NL. Characteristics of renal papillae in kidney stone formers. Minerva Urol Nefrol 2016;68:496-515. [PubMed]
  21. Béïque L, Witherspoon L, Zvonar R, et al. Duration of Antibiotic Therapy in Sepsis Secondary to Urinary Stones: A Retrospective Observational Study. Can J Hosp Pharm 2019;72:331-3. [Crossref] [PubMed]
  22. Rashid AO, Fakhulddin SS. Risk factors for fever and sepsis after percutaneous nephrolithotomy. Asian J Urol 2016;3:82-7. [Crossref] [PubMed]
  23. Sun Y, Kang J, Tao Z, et al. Effect of endoplasmic reticulum stress-mediated excessive autophagy on apoptosis and formation of kidney stones. Life Sci 2020;244:117232. [Crossref] [PubMed]
  24. An L, Wu W, Li S, et al. Escherichia coli Aggravates Calcium Oxalate Stone Formation via PPK1/Flagellin-Mediated Renal Oxidative Injury and Inflammation. Oxid Med Cell Longev 2021;2021:9949697. [Crossref] [PubMed]
  25. Pan Y, Cai W, Huang J, et al. Pyroptosis in development, inflammation and disease. Front Immunol 2022;13:991044. [Crossref] [PubMed]
  26. Wen S, Deng F, Li L, et al. VX-765 ameliorates renal injury and fibrosis in diabetes by regulating caspase-1-mediated pyroptosis and inflammation. J Diabetes Investig 2022;13:22-33. [Crossref] [PubMed]
  27. Gan XG, Wang ZH, Xu HT. Mechanism of miRNA-141-3p in Calcium Oxalate-Induced Renal Tubular Epithelial Cell Injury via NLRP3-Mediated Pyroptosis. Kidney Blood Press Res 2022;47:300-8. [Crossref] [PubMed]
  28. Song Z, Zhang Y, Gong B, et al. Long noncoding RNA LINC00339 promotes renal tubular epithelial pyroptosis by regulating the miR-22-3p/NLRP3 axis in calcium oxalate-induced kidney stone. J Cell Biochem 2019;120:10452-10462. [Crossref] [PubMed]
  29. Ferry SA, E, Holm S, Ferry BM, et al. High Diagnostic Accuracy of Nitrite Test Paired with Urine Sediment can Reduce Unnecessary Antibiotic Therapy. Open Microbiol J 2015;9:150-9. [Crossref] [PubMed]
  30. Bozkurt IH, Aydogdu O, Yonguc T, et al. Predictive Value of Leukocytosis for Infectious Complications After Percutaneous Nephrolithotomy. Urology 2015;86:25-9. [Crossref] [PubMed]
  31. Wu H, Zhu S, Yu S, et al. Early drastic decrease in white blood count can predict uroseptic shock induced by upper urinary tract endoscopic lithotripsy: a translational study. J Urol 2015;193:2116-22. [Crossref] [PubMed]
  32. Kang MJ, Shin TG, Jo IJ, et al. Factors influencing compliance with early resuscitation bundle in the management of severe sepsis and septic shock. Shock 2012;38:474-9. [Crossref] [PubMed]

(English Language Editor: J. Gray)

Cite this article as: Chen X, Li S, Shi C, Zhang W, Liu Z, Jiang J, Zhang Y, Chen Z, Zheng B, Zhu H. Risk factors and predictors of urogenous sepsis after percutaneous nephrolithotomy for idiopathic calcium oxalate nephrolithiasis. Transl Androl Urol 2023;12(6):1002-1015. doi: 10.21037/tau-23-219

Download Citation