Risk factors and predictors of urogenous sepsis after percutaneous nephrolithotomy for idiopathic calcium oxalate nephrolithiasis
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.
Table 1
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).
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).
Table 2
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
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
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
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
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 (). 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
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
- Vasudevan V, Samson P, Smith AD, et al. The genetic framework for development of nephrolithiasis. Asian J Urol 2017;4:18-26. [Crossref] [PubMed]
- 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]
- 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]
- 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]
- 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]
- 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]
- Michel MS, Trojan L, Rassweiler JJ. Complications in percutaneous nephrolithotomy. Eur Urol 2007;51:899-906; discussion 906. [Crossref] [PubMed]
- Wagenlehner FM, Pilatz A, Naber KG, et al. Therapeutic challenges of urosepsis. Eur J Clin Invest 2008;38:45-9. [Crossref] [PubMed]
- 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]
- 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]
- 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]
- 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]
- Karki R, Kanneganti TD. Diverging inflammasome signals in tumorigenesis and potential targeting. Nat Rev Cancer 2019;19:197-214. [Crossref] [PubMed]
- 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]
- Wang B, Yin Q. AIM2 inflammasome activation and regulation: A structural perspective. J Struct Biol 2017;200:279-82. [Crossref] [PubMed]
- Man SM, Kanneganti TD. Regulation of inflammasome activation. Immunol Rev 2015;265:6-21. [Crossref] [PubMed]
- Khan SR, Pearle MS, Robertson WG, et al. Kidney stones. Nat Rev Dis Primers 2016;2:16008. [Crossref] [PubMed]
- 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]
- 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]
- Marien TP, Miller NL. Characteristics of renal papillae in kidney stone formers. Minerva Urol Nefrol 2016;68:496-515. [PubMed]
- 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]
- Rashid AO, Fakhulddin SS. Risk factors for fever and sepsis after percutaneous nephrolithotomy. Asian J Urol 2016;3:82-7. [Crossref] [PubMed]
- 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]
- 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]
- Pan Y, Cai W, Huang J, et al. Pyroptosis in development, inflammation and disease. Front Immunol 2022;13:991044. [Crossref] [PubMed]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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]
- 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)