key: cord-0018676-2a2ebitr
authors: Lu, Furong; Chen, Hongpeng; Hong, Yuan; Lin, Yao; Liu, Lizhen; Wei, Ning; Wu, Qinyan; Liao, Shuanglin; Yang, Shuai; He, Junbing; Shao, Yiming
title: A gain-of-function NLRP3 3′-UTR polymorphism causes miR-146a-mediated suppression of NLRP3 expression and confers protection against sepsis progression
date: 2021-06-25
journal: Sci Rep
DOI: 10.1038/s41598-021-92547-8
sha: ea0a36b5096fb9d09db89b388fa3cc162931b8eb
doc_id: 18676
cord_uid: 2a2ebitr

Nucleotide-binding domain and leucine-rich repeat (LRR)-containing family protein 3 (NLRP3) regulated the maturation of inflammation-related cytokines by forming NLRP3 inflammasome, which plays pivotal roles in sepsis pathogenesis. In this study, we evaluated the genetic association of NLRP3 polymorphisms with sepsis (640 patients and 769 controls) and characterized the impact of NLRP3 polymorphisms on NLRP3 expression and inflammatory responses. No significant differences were observed in genotype/allelic frequencies of NLRP3 29940G>C between sepsis cases and controls. The G allele was significantly overrepresented in patients with septic shock than those in sepsis subgroup, and the GC/GG genetypes were related to the 28-day mortality of sepsis. Lipopolysaccharide challenge to peripheral blood mononuclear cells showed a significant suppression of NLRP3 mRNA expression and release of IL-1β and TNF-α in CC compared with the GC/GG genotype category. Functional experiments with luciferase reporter vectors containing the NLRP3 3′-UTR with the 29940 G-to-C variation in HUVECs and THP-1 cells showed a potential suppressive effect of miR-146a on NLRP3 transcription in the presence of the C allele. Taken together, these results demonstrated that the 29940 G-to-C mutation within the NLRP3 3′-UTR was a gain-of-function alteration that caused the suppression of NLRP3 expression and downstream inflammatory cytokine production via binding with miR-146a, which ultimately protected patients against susceptibility to sepsis progression and poor clinical outcome.

Interleukin-1 beta IL- 6 Interleukin-6 cDNA Complementary deoxyribonucleic acid PBMC Peripheral blood mononuclear cell HUVEC Human umbilical vein endothelial cell LPS Lipopolysaccharides QUANTO Quantity adjusting option Sepsis is a public health problem that is primarily caused by bacterial and fungal infections that is mediated by serious systematic inflammatory responses and leads to septic shock and multiple organ dysfunction syndrome (MODS) 1 . Even though extensive progress has been made in the development of therapeutic strategies for sepsis, it remains the main cause of mortality in intensive care units because of the extreme heterogeneity in disease outcomes possibly due to host genetic characteristics 2 . It is clear that polymorphisms in genes encoding inflammatory mediators or innate immune effectors such as IL-27, MCP-1 and ADAM17 contribute considerably to the disease course and sepsis outcome [3] [4] [5] . Thus, identifying the responsible variants will help us to elucidate genetic mechanisms that determine sepsis heterogeneity and potentially reveal new therapeutic targets for the treatment of sepsis.

In the initial stage of sepsis pathogenesis, the damage-associated molecular patterns (DAMPs) and pathogenassociated molecular patterns (PAMPs) released from microbial components and damaged tissues may interact with germline-encoded receptors called pattern recognition receptors (PRRs) and trigger chronic inflammatory responses 6 . Among these receptors, nucleotide-binding domain and leucine-rich repeat (LRR)-containing family protein 3 (NLRP3) can recruit and regulate the activation of proinflammatory caspase-1 by forming protein complexes called inflammasomes, ultimately leading to cell pyroptosis and the generation of mature forms of the inflammatory cytokines interleukin (IL)-18, IL-1β, and IL-33 7, 8 . The activation of the NLRP3 inflammasome requires NLRP3 expression, which is stimulated by nuclear factor-kappa B (NF-kB) in response to exogenous and endogenous danger signals such as high glucose and bacterial lipopolysaccharide (LPS) 9, 10 . The NLRP3 inflammasome has been regarded as a regulator of the pathomechanisms of various inflammation-related diseases (e.g., ulcerative colitis and sepsis), as well as a disordered host response to pathogenic infections (e.g., S. pyogenes and Staphylococcus aureus) [11] [12] [13] . Genomic deletion of NLRP3 inhibited sepsis-induced inflammatory responses and protected against sepsis in mice through the production of LXB4 and inhibition of proinflammatory cytokines 14 . In addition, NLRP3 inhibitors or gene knockout in mice resulted in increased neutrophil-dependent bacterial clearance, decreased concentrations of IL-18 and IL-1β and protection against septic lethality 15, 16 . These lines of evidence demonstrated a significant role of NLRP3 in sepsis pathophysiology, and blocking NLRP3 may be a promising therapeutic strategy.

The human NLRP3 gene is located on chromosome 1q44 and contains 9 coding exons. Recently, increasing evidence has shown that genetic mutations within the key domains of the NLRP3 gene play pivotal roles in the modulation of immune inflammatory responses, ultimately resulting in a genetic predisposition to various acute and chronic inflammatory diseases 17 . Approximately 60 genetic variants within the entire NLRP3 gene have been identified so far, among which rs10754558 (29940G>C), rs35829419 (Q705K), rs10925019, rs4925648 and rs4612666 have been well studied in association with susceptibility to ulcerative colitis, Alzheimer's disease, Crohn's disease and acne vulgaris [18] [19] [20] [21] . The 29940 G>C polymorphism within the 3′-UTR of the NLRP3 gene was significantly associated with NLRP3 transcription levels in peripheral blood mononuclear cells (PBMCs) and conferred susceptibility to ischemic stroke 22 . The gain-of-function variant Q705K within exon 3 of the NLRP3 gene generates hyperactive inflammasomes, along with the maturation and overproduction of IL-1β, leading to inflammatory damage in response to infection 23, 24 . However, little is known about the clinical relationship of the NLRP3 genetic variations with sepsis susceptibility and development, and the underlying mechanisms involved in the regulatory effects of these genetic variants on sepsis-induced inflammatory processes remain to be explored.

In the present study, 640 septic patients and 769 controls were enrolled to explore the clinical relationship between the NLRP3 genetic variant 29940G>C and the susceptibility and progression of sepsis in a Han Chinese population. Mechanistic investigations in vitro were conducted to determine the impact of this functional NLRP3 polymorphism on NLRP3 transcription and sepsis-associated inflammatory responses. Our results provide a new understanding of the genetic effect of the NLRP3 inflammasome on sepsis pathophysiology and new approaches for risk assessment, as well as therapeutic strategies for sepsis patients. Besides, this study was restricted to the Chinese Han population with small sample size. Thus our results may not be generalized to other populations, and need to be further validated in a larger and multiethnic population of subjects. Genotyping of NLRP3 polymorphism. TIANamp Blood DNA kit (Tiangen Biotech, Beijing, China) was applied to prepare genomic DNA from blood leukocytes of sepsis patients and healthy individuals. The iMLDR™ multiple SNP typing kit (Shanghai Tanhao biotechnology co., LTD) was applied to genotype the NLRP3 genetic variant rs10754558. The experimental results collected from the sequencer were analysed by GeneMapper 4.0 (Applied Biosystems, USA). Ten percent of samples were randomly selected for re-examination to verify the accuracy.

Preparation of PBMCs for cell culture and LPS stimulation. A total of 50 healthy controls were randomly selected for isolation of PBMCs. After isolation from the peripheral blood by Lymphoprep™ (Axis-Shield PoCAS, Oslo, Norway), the PBMCs were stored at − 80 °C until used for the in vitro LPS-stimulated experiments. RPMI-1640 medium (Thermo Fisher Scientific) was used for cell culture. Approximately 5 × 10 5 cells with 500 ng/mL LPS were cultured in 12-well plates (n = 25), and treatment with PBS vehicle was set as the control (n = 25). After 8 h of incubation, cells were separated for qRT-PCR analysis and supernatants were separated for cytokine measurements. Cell culture and transfection. Human umbilical vein endothelial cells (HUVECs) and THP-1 cells were obtained from Shanghai Institute of Cell Biology. The cells were incubated in RPMI-1640 medium supplemented with FBS (10%) and penicillin/streptomycin (1%). Lentivirus (LV) packaging of LV3-miR-146a (miR-146a overexpression) with sequence 5′-TGA GAA CTG AAT TCC ATG GGTT-3′, LV3-miR-580-3p (miR-580-3p overexpression) with sequence 5′-TTG AGA ATG ATG AAT CAT TAGG-3′, LV3-miR-589-5p (miR-589-5p overexpression) with sequence 5′-TGA GAA CCA CGT CTG CTC TGAG-3′ and LV3-NC (negative control) with sequence 5′-TTC TCC GAA CGT GTC ACG T-3′ was carried out by GenePharma (Shanghai, China). The cells were infected by LV3-miR-146a, LV3-miR-580-3p or LV3-miR-589-5p to establish stable expression cell strain, and infection with LV3-NC was set as the control cells. After 48 h of infection, qRT-PCR analysis was carried out to determine the miR-146a, miR-580-3p and miR-589-5p expression.

A 439-bp mRNA in the 3′UTR of NLRP3 containing rs10754558 G-to-C mutant was submitted to the RNAfold Web Server (http:// rna. tbi. univie. ac. at/ cgi-bin/ RNAfo ld. cgi) setting default parameters to produce positional entropies, centroid structures and Potential minimum free energy (MFE) structures. Next, a fragment with 313 bp length in the 3′UTR region of NLRP3 containing rs10754558 G-to-C mutant was inserted downstream to luciferase reporter gene in pmirGLO luciferase reporter vectors. Lipofectamine 3000 (Invitrogen, MD, USA) was used for cell transfection, and dualluciferase assay kit (Beyotime, Shanghai, China) was used to measure the activities of luciferase in a SpectraMax L Chemiluminescent reader (Molecular Devices, USA).

Enzyme linked immunosorbent assay (ELISA). The IL-1β, IL-6, and TNF-α concentrations in supernatants separated from studied subjects plasma or cell culture medium were detected by each specific ELISA kit (Boster Biological Technology Co., Ltd., Wuhan, China). Infinite F50 microplate reader (Tecan, Switzerland) was used to detect the absorbance at 450 nm.

Data analysis and statistics. GraphPad Prism 6.0 software (GraphPad Software Inc., San Diego, CA, USA) was applied to data analysis. Fisher's exact test or Chi-squared test was performed to calculate the geno-

The association of NLRP3 29940G>C polymorphism with sepsis susceptibility and progression. The genotypes and allelic distributions of the NLRP3 29940G>C polymorphism were consistent with the Hardy-Weinberg equilibrium in sepsis patients and healthy controls (All p > 0.05; Supplementary Table S1 ). As shown in Table 2 , no significant differences in the genotype/allelic frequencies of the 29940G>C polymorphism were found between sepsis cases and healthy individuals (all p > 0.05). When the cases were separated www.nature.com/scientificreports/ into sepsis and septic shock, the data showed that patients with the GG/GC genotypes were associated with a significantly increased risk of progression from sepsis to septic shock (GG/GC vs. CC: p = 0.005, OR = 1.608, 95% CI 1.152-2.245, Table 3 ). Similarly, G allele carriers had an increased risk of sepsis progression (G allele vs. C allele: p = 0.003, OR = 1.402, 95% CI 1.120-1.754). These significant differences were still observed after the Benjamini-Hochberg correction.

The association of the NLRP3 29940G>C polymorphism with the mortality of septic patients. We further evaluated the effect of NLRP3 29940G>C polymorphism on the 28-day mortality of septic patients. As presented in Table 4 , the frequencies of the G allele and GG/GC genotypes in 28-day nonsurviving sepsis patients were significantly overrepresented compared to those in the surviving sepsis patients (GG/GC vs. CC: p = 0.008, OR = 1.700, 95% CI 1.146-2.522; G allele vs. C allele: p = 0.010, OR = 1.395, 95% CI 1.083-1.795). These differences still exhibited statistically significance after the Benjamini-Hochberg multipletesting corrections (GG/GC vs. CC: p = 0.032; G allele vs. C allele: p = 0.020). Furthermore, Kaplan-Meier survival analysis showed that 28-day survival exhibited much worse in septic patients with the GG/GC genotypes than in septic patients carrying the CC genotype (log-rank = 7.518, p = 0.006; Fig. 1 ).

Effect of the NLRP3 29940G>C polymorphism on NLRP3 expression and cytokine production. A total of 50 septic patients and 50 healthy controls were randomly selected to evaluate the potential association between the NLRP3 29940G>C polymorphism and the expression of NLRP3 and related inflamma- www.nature.com/scientificreports/ tory cytokines. As presented in Fig. 2 , sepsis patients exhibited significantly higher expression and production of NLRP3 and the proinflammatory cytokines (TNF-α, IL-6 and IL-1β) than healthy controls, and these levels were significantly increased with sepsis progression. When sepsis cases were separated into subgroups by genotype, carriers of the GG/GC genotypes at the NLRP3 29940G>C site exhibited significantly higher NLRP3 mRNA and TNF-α levels than carriers with the CC genotype. No significant differences in NLRP3 and cytokine levels were found in healthy controls with different genotypes.

To further validate whether the NLRP3 29940G>C polymorphism was likely influence NLRP3 expression at the transcriptional level, qRT-PCR analysis was conducted in LPS-stimulated PBMCs from 50 healthy controls to demonstrate the connection between genotype and phenotype. Our results showed that the NLRP3 29940G>C variant was closely associated with the responsiveness of PBMCs to LPS stimulation. LPS stimulation induced significantly higher NLRP3 expression and proinflammatory cytokine (IL-1β and TNF-α) production in PBMCs www.nature.com/scientificreports/ with the 29940GG/GC genotypes compared with those with the CC genotype, while these differences were not observed without LPS stimulation (Fig. 3) . However, no significant difference in IL-6 production was observed between PBMCs with different 29940G>C genotypes.

Influence of the NLRP3 3′-UTR polymorphism 29940G>C on NLRP3 mRNA stability. Considering the location of the 29940G>C polymorphism within the 3′-UTR of the NLRP3 gene, we further examined the effect of the sepsis-associated G-to-C mutation at the rs10754558 site on NLRP3 mRNA stability. As shown in Fig. 4 , NLRP3 mRNA secondary structure prediction analysis with the RNAfold Web Server indicated that the MFE value in the 3′-UTR of NLRP3 with the G-to-C mutation at the rs10754558 site had a notable variation from ∆G = − 141.40 to − 137.80 kcal/mol. The centroid secondary structure MFE value in this predicted mRNA secondary structure with the G allele of the 29940G>C polymorphism (∆G = − 80.96 kcal/mol) was also observed to be lower than that with the C allele (∆G = − 54.05 kcal/mol). This G-to-C mutation in the internal base pairing contributed to markedly different bubble and stem-loop structures, ultimately leading to a possible decrease in NLRP3 mRNA stability. The human NLRP3 gene containing nine coding exons is located on chromosome 1q44 (247416156-247449108), and the 29940G>C polymorphism within 3′-UTR of NLRP3 gene was positioned at chromosome 1 (247448734) according to the GRCh 38. p12 Primary Assembly (Fig. 5A ). Next, a luciferase assay with pmirGLO luciferase reporter vectors containing the G-to-C mutation at the rs10754558 site (Fig. 5B ) was conducted in HUVECs and THP-1 cells to gain further insights into the functional impact of this NLRP3 3′-UTR polymorphism on NLRP3 expression. The pmirGLO luciferase reporter vector with the NLRP3 29940G allele showed significantly higher activities than the vector with the C allele (Fig. 5C,D) . www.nature.com/scientificreports/

ing miR-146a-5p binding. As the 29940G>C polymorphism was located in the 3′-UTR of the NLRP3 gene, we focused on its effect on gene interactions with miRNAs. Bioinformatics analyses predicted by miR-NASNP-v3 (http:// bioin fo. life. hust. edu. cn/ miRNA SNP/) and TargetScan (http:// www. targe tscan. org/) revealed a target binding site of miR-146a-5p enabled by the 29940 G-to-C mutation within the 3′-UTR of the NLRP3 gene ( Fig. 6) 27 . We next detected the miR-146a expression in 50 septic patients and 50 healthy controls. Our results showed that sepsis patients exhibited significantly lower expression level of miR-146a compared to the healthy controls (Fig. 6A) . No significant differences in miR-146a expression were observed in patients with different types of sepsis or different genotypes of rs10754558 polymorphism (Fig. 6B,C) . Similarly, in vitro LPS stimulation induced decreased expression of miR-146a in the PBMCs isolated from 50 healthy subjects, while differential levels of miR-146a were not detected in PBMCs with different genotypes of NLRP3 polymorphism (Fig. 6D,E) .

To validate the potential role of miR-146a-5p in the attenuation of NLRP3 expression mediated by the NLRP3 29940 G-to-C mutation, HUVECs and THP-1 cells were infected with LV3-miR-146a (miR-146a overexpression vector) and LV3-NC (negative control) to establish a stably overexpressing cell line (Fig. 6F-H) . Next, pmirGLO-NLRP3-3′-UTR-29940G and pmirGLO-NLRP3-3′-UTR-29940C were transfected into the cells for 48 h. As shown in Fig. 6I ,J, miR-146a-5p significantly decreased the luciferase reporter activity of the pmirGLO-NLRP3-3′-UTR-29940C vector, while no obvious variations were found in the pmirGLO-NLRP3-3′-UTR-29940G vector. With regards to the other miRNAs predicted by miRNASNP-v3 and TargetScan, neither miR-580-3p nor miR-589-5p overexpression exerted a significant effect on the activities of the luciferase reporter vectors containing the NLRP3 3′-UTR with 29940 G-to-C variation in THP-1 cells (Supplementary Fig. S1 ). These results indicated that NLRP3-29940G>C is a gain-of-function alteration that causes the suppression of NLRP3 expression by altering the binding of miR-146a-5p to the 3′-UTR of NLRP3. However, no significant difference in NLRP3 expression was observed between THP-1 cells infected with LV3-NC and LV3-miR-146a (Fig. 6K).   p36.3 p36.2 p36.1 p35 p34.2 p33 p32.3 p31. 

C C C C C SalI B D rs10754558 G rs10754558 C C C C C rs10754558

It has been demonstrated that inappropriate systemic inflammatory responses contribute to the poor prognosis of sepsis patients 28 . A number of publications have indicated that a significant portion of the variability in predisposition to the occurrence and progression of sepsis was potentially attributed to genetic factors, particularly certain functional genetic variants 29 . Increasing evidence from our previous studies and others has revealed an important effect of genetic variations in at least 40 genes encoding inflammation-related factors on interindividual differences in the immune inflammatory response and susceptibility to sepsis outcomes [30] [31] [32] . In this study, we explored the clinical relevance of NLRP3 genetic variation, rs10754558 (29940G>C), to sepsis susceptibility and progression in a Han Chinese population (640 septic patients and 769 controls) for the first time. Our data indicated that patients with the 29940 GG/GC genotypes exhibited significantly increased risk of progression from sepsis to septic shock, and the 28-day survival in patients with sepsis carrying the GG/GC genotypes exhibited much worse than in the CC genotype carriers, which indicated that the 29940G>C may act The miR-146a expression in septic patients (n = 50) and healthy controls (n = 50) with different rs10754558 polymorphism; (D, E) PBMCs isolated from 50 healthy subjects with different rs10754558 genotypes were stimulated with 500 ng/mL of LPS (n = 25) or PBS (n = 25) in vitro, and then miR-146a expression was detected; (F) Bioinformatics predicted that the rs10754558 G-to-C mutation was located at the binding site between 3′UTR of NLRP3 and miR-146a-5p; (G, H) HUVECs and THP-1 cells were infected by LV3-miR-146a (miR-146a overexpression) or LV3-NC, and the infection efficiency of LV3-miR-146a was determined by qRT-PCR analysis; (I, J) The effect of miR-146a on the transcriptional activities of rs10754558 G>C in HUVECs and THP-1 cells were determined by dual-luciferase report assays; (K) The effect of miR-146a on NLRP3 expression in THP-1 cells was determined by qRT-PCR analysis. Bar graphs indicate the mean ± SEM for a minimum of three independent experiments. *p < 0.05, **p < 0.01, ***p < 0.001. www.nature.com/scientificreports/ as a risk factor for the development of sepsis and poor clinical outcome. Furthermore, functional experiments suggested that this sepsis-associated NLRP3-29940G>C polymorphism might be a gain-of-function alteration that suppresses NLRP3 expression and downstream cytokine production by altering the binding of miR-146a-5p to the 3′-UTR of NLRP3. NOD-like receptors (NLRs) are pattern recognition receptors that recognize various pathogens and damageassociated molecular patterns and are important in the context of sepsis 33 . Among these receptors, NLRP3 (also known as NALP3 and cryopyrin) forms a cytosolic signaling complex termed the NLRP3 inflammasome, which is an important molecular platform that controls the activation of proinflammatory caspase-1 and the maturation of the potent inflammatory mediators IL-18 and IL-1β, which play pivotal roles in sepsis pathogenesis 9 . The expression of NLRP3 inflammasome mRNA is significantly upregulated in patients with sepsis and various murine models of sepsis and is closely associated with the disease course and the clinical prognosis of sepsis patients 34, 35 . Consistently, our data showed that sepsis patients exhibited significantly higher NLRP3 gene expression and protein production than healthy volunteers, and these levels were also increased with disease severity. These results verified the pivotal involvement of NLRP3 in the clinical course of sepsis as a proinflammatory regulator, as well as an indicator of disease severity.

Several studies have paid close attention to in-depth genetic analysis of the NLRP3 gene in various inflammation-related diseases [18] [19] [20] [21] . Herrmann et al. discovered that the rs7525979 C-to-T mutation in the NLRP3coding region conferred protection against susceptibility to Parkinson's disease by impacting NLRP3 translation efficiency 36 . The 29940 G>C genetic variation in the 3′-UTR of NLRP3 attenuated the stability and expression of NLRP3 mRNA in a gain-of-function manner, which resulted in susceptibility to several inflammation-related disorders, such as asthma, ischemic stroke and rheumatoid arthritis 22, 37, 38 . Other studies have demonstrated that the Q705K genetic variation in exon 3 of the NLRP3 gene is also a gain-of-function mutation that generates overactive NLRP3 inflammasomes and overproduction of IL-18, IL-1β and TNF-α, which are closely associated with atopic dermatitis, ulcerative colitis and abdominal aortic aneurysms 39, 40 . However, little is known about the clinical relevance of the gain-of-function variant 29940 G>C to the susceptibility and progression of sepsis.

The initial results in this study showed no statistically significant differences in the genotype/allelic frequencies of the NLRP3 29940G>C variation between sepsis patients and healthy volunteers, suggesting that there was no effect on the onset of sepsis. Novel findings from additional stratification analyses showed that the frequencies of the G allele and the GG/GC genotypes in the septic shock subgroup were significantly overrepresented compared to those in the sepsis subgroup. Carriers of 29940 GG/GC exhibited significantly increased NLRP3 mRNA expression compared to that of the CC genotype carriers, indicating a role of the 29940 G>C polymorphism in the susceptibility to sepsis progression by altering NLRP3 gene transcription. As the mRNA secondary structure is an important feature of the functional efficiency of cis-acting factors in the 3′-UTR, we examined the impact of the 29940 G-to-C mutation in the NLRP3 3′-UTR on NLRP3 mRNA stability. RNA secondary structure prediction indicated that the G-to-C mutation in the internal base pairing contributed to markedly different bubble and stem-loop structures and an increased MFE, ultimately leading to a possible decrease in NLRP3 mRNA stability. Next, a luciferase reporter assay was conducted to evaluate the functional impact of the 29940 G-to-C mutation on NLRP3 expression and showed that the allele-specific luciferase vector with the sepsis risk-associated G allele exhibited significantly higher activities than the C allele, which corroborated the findings of several previous studies 22, 37 . To further confirm these results, we performed in vitro LPS stimulation experiments to verify the impact of the 29940 G-to-C mutation on NLRP3 mRNA expression. As predicted, LPS-induced NLRP3 gene expression was significantly increased in PBMCs with the sepsis-associated risk GG/ GC genotypes compared with the CC genotype. However, no significant differences in NLRP3 expression or cytokine production were observed in PBMCs with different genotypes in the absence of LPS stimulation. These results suggested that the hereditary effects of the NLRP3 29940 G-to-C mutation made a real difference in the response to exogenous and endogenous danger signals such as bacterial LPS, which predisposed certain sepsis patients to disease progression from sepsis to sepsis shock rather than the occurrence of sepsis.

MiRNAs are small noncoding RNAs that bind to specific mRNA molecules and cause mRNA degradation or suppression of target gene expression, and they play pivotal roles in the modulation of multiple cellular processes, such as proliferation, differentiation, inflammation and apoptosis 41 . Increasing evidence has demonstrated the important impact of specific genetic variations within the 3′-UTR on the binding of miRNAs to cause the gain or loss of regulation by certain miRNAs, ultimately resulting in susceptibility to multiple diseases 42, 43 . According to bioinformatics analyses predicted by miRNASNP-v3 (http:// bioin fo. life. hust. edu. cn/ miRNA SNP/) and a previous study 27, 38 , the 29940 G-to-C mutation in the 3′-UTR of NLRP3 might establish a new functional binding site for miR-146a-5p, which contributes to potential miR-146-5p recruitment and binding. Research so far has demonstrated an anti-inflammatory effect of miR-146a-5p on the regulation of multiple signaling pathways involved in the modulation of various inflammation-related cytokines 44 . Our results showed that miR-146a-5p expression was significantly downregulated under a sepsis or in vitro LPS stimulation, in accordance with the findings of our previous study 45 . However, no significant differences in miR-146a expression were observed in patients with different types of sepsis or different genotypes of NLRP3 polymorphism. The luciferase reporter experiment showed a significantly suppressive effect of miR-146a on NLRP3 transcription in the presence of the mutant-type C allele in HUVEC and THP-1 cells, while these differences were not observed in the presence of the wild-type G allele. These results further validated that the 29940 G>C polymorphism in the 3′-UTR of the NLRP3 gene was a gain-of-function alteration causing the suppression of NLRP3 expression by miR-146a-5p, which conferred protection against progression from sepsis to septic shock. However, we failed to detect the inhibitory effect of miR-146a on NLRP3 expression in a monocytic THP-1 cell line, which might be attributed to the wild-type G allele at rs10754558 site in THP-1 cells. Our further study would explore this gain-of-function mutation in NLRP3 gene by constructing the NLRP3 gene point mutation cell line via CRISPR/Cas9 technology. www.nature.com/scientificreports/ Upon sensing PAMPs and DAMPs, NLRP3 assembles a group of high-molecular-weight protein complexes that can activate caspase 1-mediated release of bioactive proinflammatory cytokines, as well as nonclassical secretion of various cytosolic proteins, causing excessive inflammatory responses and pyroptosis 7 . NLRP3 activation contributes to the cascade release of inflammatory cytokines, including IL-33, IL-18, IL-1β and HMGB1, and prostaglandins and leukotrienes are also released by the NLRP3/Caspase-1 axis 46, 47 . NLRP3 inhibitors or genetic deficiency have been shown to decrease the production of various inflammatory cytokines, including TNF-α, IL-6, IL-1β, MIP-2 and lipid mediators, and provide protection against septic lethality [14] [15] [16] . Consequently, we measured the production of IL-6, TNF-α and IL-1β in sepsis patients and healthy controls and in PBMCs under in vitro LPS stimulation to investigate the effect of the NLRP3 29940 G>C variation on the production of these related cytokines. It has been reported that the NLRP3 inflammasome is significantly activated due to a gain-of-function alteration in rs35829419 in NLRP3, causing excessive IL-1β, IL-18 and TNF-α production and exaggerated inflammatory responses in a caspase-1-dependent manner 40 . Our results showed a close association between TNF-α production but not IL-1β or IL-6 with the NLRP3 29940 G>C variation in sepsis patients, which might be due to the complex inflammatory system in the host. Furthermore, PBMCs with the sepsis-associated risk allele 29940G in NLRP3 exhibited significantly enhanced expression of NLRP3 and release of IL-1β and TNF-α compared to those with the C allele under LPS stimulation. These results further confirmed that the 29940 G-to-C mutation in the NLRP3 3′-UTR suppressed NLRP3 expression and downstream inflammatory cytokine production by binding with miR-146a in a gain-of-function manner, which ultimately protected patients against susceptibility to sepsis progression and exhibited potentially important diagnostic and therapeutic value (Fig. 7) .

Several limitations in the present study should be mentioned. The undiscovered underlying confounding factors, including pre-existing conditions in sepsis patients, could have influenced those observed differences in our results, even though exclusion and inclusion criteria were used to obtain more homogeneity. Furthermore, this study was restricted to the Chinese Han population with small sample size. According to the international HapMap project, the frequency of G allele at rs10754558 site in Han Chinese descent (0.41) is similar to those in several populations, such as Finnish in Finland (0.42) and British in England and Scotland (0.41), while it differed from those in other populations, such as Americans of African Ancestry (0.24) and African Carribbeans in Barbados (0.31). Therefore, one should exercise caution when generalizing the findings to other ethnic groups, and these findings need to be further validated in a larger and multiethnic population of subjects.

The data in this study indicated that the 29940 G>C polymorphism in the 3′-UTR of NLRP3 conferred susceptibility to the development of sepsis and a poor clinical outcome. This 29940 G-to-C mutation was demonstrated to be a functional polymorphism that cause the suppression of NLRP3 transcription and downstream inflammatory Figure 7 . The molecular mechanism of NLRP3 3′-UTR polymorphism-mediated miR-146a-5p modulation on NLRP3 expression in the progression of sepsis. The 29940 G-to-C mutation within the NLRP3 3′-UTR was a gain-of-function alteration that caused the suppression of NLRP3 expression and downstream inflammatory cytokine production via binding with miR-146a-5p, which ultimately protected patients against susceptibility to sepsis progression. www.nature.com/scientificreports/ cytokine production by binding with miR-146a-5p in a gain-of-function manner, which ultimately conferred protection against progression from sepsis to septic shock. These results may provide a new understanding of the genetic effect of the NLRP3 inflammasome on sepsis pathophysiology and new approaches for risk assessments and therapeutic strategies for sepsis patients.

Getting sepsis therapy right

The pathogenesis of sepsis and potential therapeutic targets

The interleukin-27-964A>G polymorphism enhances sepsis-induced inflammatory responses and confers susceptibility to the development of sepsis

Association study of MCP-1 promoter polymorphisms with the susceptibility and progression of sepsis

Association study between promoter polymorphisms of ADAM17 and progression of sepsis

The immune system's role in sepsis progression, resolution, and long-term outcome

Metabolic regulation of NLRP3

Mechanism and regulation of NLRP3 inflammasome activation

The NLRP3 inflammasome and its role in sepsis development

High glucose and lipopolysaccharide prime NLRP3 inflammasome via ROS/TXNIP pathway in mesangial cells

Targeting the NLRP3 inflammasome in inflammatory diseases

Activation of the Nlrp3 inflammasome by Streptococcus pyogenes requires streptolysin O and NF-kappa B activation but proceeds independently of TLR signaling and P2X7 receptor

Orchestration of human macrophage NLRP3 inflammasome activation by Staphylococcus aureus extracellular vesicles

NLRP3 inflammasome deficiency protects against microbial sepsis via increased lipoxin B4 synthesis

Deletion of Nlrp3 augments survival during polymicrobial sepsis by decreasing autophagy and enhancing phagocytosis

Critical role for NALP3/CIAS1/Cryopyrin in innate and adaptive immunity through its regulation of caspase-1

Gene polymorphisms in the NALP3 inflammasome are associated with interleukin-1 production and severe inflammation: Relation to common inflammatory diseases?

Association of NLRP3 single nucleotide polymorphisms with ulcerative colitis: A case-control study

NLRP3 polymorphisms are associated with late-onset Alzheimer's disease in Han Chinese

The genetics of NOD-like receptors in Crohn's disease

Genetic association between the NLRP3 gene and acne vulgaris in a Chinese population

A polymorphism within the 3′UTR of NLRP3 is associated with susceptibility for ischemic stroke in Chinese Population

Interaction of the inflammasome genes CARD8 and NLRP3 in abdominal aortic aneurysms

Polimorphisms in inflammasome genes are involved in the predisposition to systemic lupus erythematosus

The third international consensus definitions for sepsis and septic shock (Sepsis-3)

A functional polymorphism-mediated disruption of EGR1/ADAM10 pathway confers the risk of sepsis progression

miRNASNP-v3: A comprehensive database for SNPs and disease-related variations in miRNAs and miRNA targets

Sepsis and septic shock

Host genetic variants in sepsis risk: A field synopsis and meta-analysis

Association of TNF-alpha-308 (G>A) (rs1800629) gene polymorphism with adverse outcomes of sepsis in critically ill patients

The promoter polymorphisms of receptor for advanced glycation end products were associated with the susceptibility and progression of sepsis

Association between genetic polymorphisms in the autophagy-related 5 gene promoter and the risk of sepsis

Danger of frustrated sensors: Role of Toll-like receptors and NOD-like receptors in aseptic and septic inflammations around total hip replacements

Inflammasome gene profile is modulated in septic patients, with a greater magnitude in non-survivors

P2X7 receptor induces mitochondrial failure in monocytes and compromises NLRP3 inflammasome activation during sepsis

NLRP3 expression in mesencephalic neurons and characterization of a rare NLRP3 polymorphism associated with decreased risk of Parkinson's disease. NPJ Parkinson's disease 4, 24

Associations of functional NLRP3 polymorphisms with susceptibility to food-induced anaphylaxis and aspirininduced asthma

Polymorphisms and expression of inflammasome genes are associated with the development and severity of rheumatoid arthritis in Brazilian patients

The Q705K polymorphism in NLRP3 is a gain-of-function alteration leading to excessive interleukin-1beta and IL-18 production

NLRP3 rs35829419 polymorphism is associated with increased susceptibility to multiple diseases in humans

Role of sepsis modulated circulating microRNAs. Ejifcc

MicroRNA binding site polymorphism in inflammatory genes associated with colorectal cancer: Literature review and bioinformatics analysis

The relevance of SNPs at 3′UTR region of CASP7 and miR-371b-5p associated diseases: A computational analysis

A comprehensive review on miR-146a molecular mechanisms in a wide spectrum of immune and non-immune inflammatory diseases

The functional polymorphisms of miR-146a are associated with susceptibility to severe sepsis in the Chinese population

Bacterial RNA and small antiviral compounds activate caspase-1 through cryopyrin/Nalp3

NLRP3 (NALP3, Cryopyrin) facilitates in vivo caspase-1 activation, necrosis, and HMGB1 release via inflammasome-dependent and -independent pathways

Writing-Original Draft Preparation

All data generated or analysed during this study are included in this published article. 

This work was supported by the National Natural Science Foundation of China (81772048, 82072151), the Natural Science Foundation of Guangdong Province (2020A1515010109; 2019A1515010933), the Fund of Southern Marine Science and Engineering Guangdong Laboratory (Zhanjiang) (ZJW-2019-007), the Novel Coronavirus Infection Prevention technology and product project of Zhanjiang City (2020E01007) and the Science and Technology Project of Jieyang City (xgfy021).

The authors declare no competing interests.

The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598-021-92547-8.Correspondence and requests for materials should be addressed to J.H. or Y.S.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.