Abstract

Chronic inflammation is a pivotal contributor to the liver damage mediated by hepatitis C virus (HCV). The NOD-like receptor, pyrin domain-containing 3 (NLRP3) inflammasome is activated by HCV in both hepatocytes and Kupffer cells. The aim of our study was to investigate the association of nine single-nucleotide polymorphisms in four inflammasome genes (NLRP3, CARD8, IL-1β, and IL-18) with the susceptibility to HCV infection and outcome of interferon treatment in 201 Egyptian chronic hepatitis C patients and 95 healthy controls. The genotyping was conducted using TaqMan predesigned SNP assay. In the comparative analysis, the CC genotype of the NLRP3 rs1539019 was found to be associated with the lower risk to chronic HCV infection (OR: 0.33, 95% CI: 0.17-0.62). This association was also found for the CA genotype and the A allele of the NLRP3 rs35829419 (OR: 0.18 and 0.22, respectively), in addition to the GG genotype and G allele of IL-18 rs1946518 (OR: 0.55 and 0.61, respectively). In contrast, the AA genotype of the IL-1β rs1143629 was significantly more frequent in HCV patients (OR: 1.7, 95% CI: 1-2.86). Notably, the frequency of the AA genotype of NLRP3 rs1539019 was significantly higher in patients with lack of response (NR) to the interferon treatment (OR: 1.95, 95% CI: 1-3.7). A similar association was found for both the CC genotype and C allele of the NLRP3 rs35829419 (OR: 2.78 and 2.73, respectively) and for the TT genotype and T allele of CARD8 rs2043211 (OR: 2.64 and 1.54, respectively). Yet, the IL-1β (rs1143629, rs1143634) and IL-18 (rs187238, rs1946518) polymorphisms did not show any significant association with response to interferon treatment. In conclusion, this study reports, for the first time, the association of genetic variations in NLRP3 with hepatitis C susceptibility and response to treatment in Egyptian patients. However, further large-scale studies are recommended to confirm our findings.

1. Introduction

Hepatitis C virus (HCV) is a global public health problem with more than 185 million infections worldwide [1]. The highest prevalence of HCV was reported in Egypt; approximately 12% of the Egyptian population are HCV positives, and 10% are chronically infected [2].

Hepatitis C virus often causes persistent infection in humans, which leads to chronic hepatitis in up to 60–80% of patients and can progress to liver fibrosis, cirrhosis, and eventually hepatocellular carcinoma [3].

A growing body of evidence has suggested that proinflammatory cytokines (IL-1β and IL-18) play a substantial role in both acute and chronic hepatitis C [4, 5]. This finding was supported by previous studies that reported elevated IL-1β (interleukin-1β) levels in HCV-related liver diseases than other forms of liver damage [6].

In acute infection, these cytokines have an antiviral effect [79] and persuade the expression of interferon-stimulated genes (ISGs) to increase the effectiveness of interferon actions [1012].

However, IL-1β modulates various cellular processes in chronic HCV infection, including induction of COX-2, nitric oxide, and TNF-α [13]. The persistent production of IL-1β by hepatic macrophages in chronic HCV patients would then help to recruit immune cells to the liver causing inflammation which also affects the response to interferon therapy [14, 15].

Interleukin-1β is synthesized as pro-IL-1β that is converted into the active form by the inflammasome, which is a multiprotein complex consisting mainly of a Nod-like receptor (NLR), adaptor molecule, and procaspase-1. The recruitment of the inflammasome results in the activation of caspase-1, which cleaves pro-IL-1β and pro-IL-18 into their mature form [16].

Recent studies reported production of mature IL-1β through the NLRP3 inflammasome assembly in hepatic macrophages and hepatocytes infected with HCV [5, 17]. Furthermore, the activity of this inflammasome is regulated by the caspase recruitment domain-containing protein 8 (CARD8) [18].

The combination of pegylated interferon alpha (PEG-IFN) and ribavirin (RBV) is used for treatment of chronic HCV in Egypt. The main goal of HCV therapy is to produce a sustained virological response (SVR), which is defined as the absence of HCV RNA in the serum at the end of treatment and 6 months later [19]. Unfortunately, the rate of SVR is around 50% in HCV genotype 4-infected patients [20]. Of note, genotype 4 is the major HCV genotype (>90%) among the Egyptian patients [21, 22].

The host genetic variations play a pivotal role in susceptibility to HCV infection and fate of viral infection. Several polymorphisms within cytokine genes such as IL-18, IFN, IL-10, and IL-28 have been reported to be associated with HCV susceptibility and fate of the disease [23, 24]. Despite that NLRP3 gene polymorphisms have been found to be associated with susceptibility to HIV and human T-lymphotropic virus 1 infection [25, 26], their association with HCV infection is still ambiguous. So, the objective of the current study was to investigate the association of nine single-nucleotide polymorphisms in four inflammasome-related genes (NLRP3, IL-1β, IL-18, and CARD8) with susceptibility to HCV infection and outcome.

2. Material and Methods

2.1. Participants

This study included 201 chronic HCV patients recruited from the National Hepatology and Tropical Medicine Research Institute (NHTMRI), Cairo, Egypt. In addition, 95 healthy subjects were recruited as sex- and age-matched controls. The study was carried out from 2014 to 2016.

The study had been approved by the Research Ethical Committee (REC) of National Hepatology and Tropical Medicine Research Institute, serial number 11-2014, date 13 July 2014. This committee works according to the Declaration of Helsinki for human subjects (2008). Moreover, an informed consent was obtained from all the participants. The study has been conducted in accordance with national protocols for the treatment of chronic hepatitis C, which has been issued by the National Committee for Control of Viral Hepatitis, Ministry of Health in Egypt.

The selection criteria for chronic HCV patients were the positive HCV Ab and positive HCV viremia at the time of diagnosis. Liver biopsy was conducted for every patient to assess the degree of fibrosis and for exclusion of other concurrent liver diseases. The exclusion criteria were other hepatic viral infections, schistosomiasis, autoimmune diseases, hematological disorders, alcoholic liver disease, obesity-induced liver disease, renal disease, poorly controlled diabetes, and immunologically mediated disease (ulcerative colitis, Crohn’s disease, systemic lupus erythematous, scleroderma, severe psoriasis, and rheumatic arthritis). Additionally, /day and intravenous or inhaled drug abuse were also excluded. For healthy controls, all were negative for HCV Ab, HIV, HBV, and no other concurrent diseases.

Chronic hepatitis C patients received 48 weeks of treatment with combination therapy, according to the protocol used in NHTMRI, which is composed of IFN administrated subcutaneously in a dose of 1.5 μg/kg once per week (dosage modified in case of cytopenia or symptomatic adverse events during treatment) and 1000–1200 mg of oral ribavirin per day. After discontinuation of treatment, they were followed for another 24 weeks. To determine the response to combination therapy, sera were tested for HCV RNA at the end of therapy and after 6 months of therapy. The chronic HCV patients were classified into two groups: the responder group included 86 subjects which achieved SVR defined as clearance of HCV RNA from serum at the end of treatment and for 24 weeks after the end of treatment.

The second group consisted of 115 “nonresponders/relapsers.” The nonresponder could not achieve end-of-treatment virological response (ETVR). Relapsers could successfully achieve ETVR but fail to achieve SVR (reappearance of HCV RNA).

A complete medical history and biochemical parameters were obtained for all the patients (height and weight to estimate BMI, age, ALT, AST, AFP, total bilirubin, ALP, albumin, Hb, TLC, the degree of liver fibrosis, platelet count, and creatinine); all parameters were initial before the start of treatment. Additionally, the viral load at baseline, week 12, week 24, week 48 (the end of the treatment), and 6 months after the end of treatment were also obtained.

2.2. Blood Sample and DNA Isolation

A venous blood sample of 2 mL was drawn from every participant. The blood sample was collected in sterile anticoagulant tubes. Genomic DNA was extracted from EDTA whole blood using PureLink® Genomic DNA Mini Kit (Invitrogen, USA).

2.3. Genetic Polymorphism Detection

Nine SNPs, three within the NLRP3 gene (rs1539019, rs35829419, and rs10754558), two within IL-1β (rs1143629, rs1143634), two within IL-18 (rs187238, rs1946518), and two within the CARD8 gene (rs6509365, rs2043211), were genotyped in this study. The genotyping was performed using TaqMan® Predesigned SNP Genotyping Assay (40x) (Catalog Number: 4351379); Invitrogen® by life technologies®; Applied Biosystems; USA. For the master mix, we used TaqMan® Universal Master Mix II, no UNG (2x) (Catalog Number: 4440043); Invitrogen® by life technologies®; Applied Biosystem; USA.

The PCR for detection of SNP was carried out according to the manufacturer’s instructions using Rotor-Gene®, Qiagen, Real-Time PCR, USA. Allelic discrimination was performed using Rotor–Gene Q Software version: 2.0.2 (Build 4).

2.4. Statistical Analysis

The sample size was calculated using OpenEpi info (http://www.openepi.com) with the following parameters: confidence interval , , statistical power , and minor allele frequency (MAF) estimated using Haploview bioinformatics software version 4.2, . Relative excess heterozygosity (REH) was calculated to inspect the harmony of genotype frequencies within the control groups with the Hardy-Weinberg equilibrium (HWE) directly according to Wellek et al. [27]. In addition, values from the standard exact HWE lack-of-fit test using 100,000 permutations () were also studied. HWE analysis was performed using R 3.4.0 [28] with the “Hardy Weinberg version 1.5.8” package [29].

Risk assessment statistical analysis was conducted by SPSS version 21 (SPSS, Chicago, IL, USA). Continuous variables of demographic and anthropometric characteristics were compared using Student’s -test. The genotype and allele frequencies were calculated using the gene counting method. Univariate comparisons of categorical variables were performed with the chi–square test. Allele frequencies and genotype distributions between the studied groups were compared by the chi-square test and Fisher’s exact test when expected were found for any cell. Unadjusted (crude or exposure) odds ratio and 95% confidence interval (CI) were estimated. The nominal level of statistical significance for all analyses was .

The pairwise linkage disequilibrium between SNPs located on the same chromosomes (, ) and haplotype proportions were estimated using Haploview bioinformatics software version 4.2. The correlation between the linkage disequilibrium haplotypes and HCV susceptibility were estimated using the SHEsis program (http://analysis.bio-x.cn/myAnalysis.php).

3. Results

Nine single-nucleotide polymorphisms in four inflammasome genes were genotyped in 201 chronic HCV patients and 95 healthy controls. The general characteristics of patients and controls are displayed in Table 1. The patients and the controls were age-, sex-, and BMI-matched.

The SNP characteristics are listed in Table 2. No significant deviation from the Hardy-Weinberg equilibrium was found in any of the studied SNPs in the healthy control group, indicating that the genotypes were appropriately assigned.

The genotypes and allele frequencies of all the studied SNPs are shown in Table 3. The HCV patients revealed a significantly higher AC genotype frequency for the NLRP3 rs1539019 compared to the healthy control group (OR 2.1, 95% CI: 1.3-3.5). On the other hand, the genotype CC was associated with significantly lower risk of chronic HCV infection (OR 0.33, 95% CI: 0.17-0.62).

Additionally, the frequencies of the CC genotype and C allele of the NLRP3 rs35829419 were significantly higher in the HCV patients compared to the healthy control group (OR 6.05, 95% CI: 3.4-10.66; OR 4.59, 95% CI: 2.79-7.55, respectively). Nevertheless, genotype CA and allele A were associated with decreased risk of HCV positivity. For NLRP3 rs10754558, no significant differences in genotypes or allele frequencies were found between the HCV patients and the healthy control group (Table 3).

For IL-1β rs1143629, the AA genotype was significantly more frequent in HCV patients (). However, the IL-1β rs1143634 polymorphism revealed no association with HCV infection.

In contrast to the IL-18 rs187238, which did not show any association with HCV infection, genotype GG and allele G in the rs1946518 of the same gene were significantly more frequent in the healthy control indicating a protective effect against infection (OR: 0.55 and 0.61, respectively). However, the T allele of this SNP was associated with increased risk of HCV positivity (OR: 1.63).

On the other hand, the SNPs in the adaptor molecule CARD8 (rs6509365 and rs2043211) did not reveal any significant association with HCV infection (Table 3).

The pairwise linkage disequilibrium analysis for the studied SNPs within NLRP3 is shown in Figure 1, and those within IL-1β, IL-18, and CARD8 are presented in Figure 2. It has been disclosed that the studied SNPs within NLRP3, IL-1β, and IL-18 were in linkage equilibrium. Furthermore, rs6509365 and rs2043211 within CARD8 were found to be in linkage disequilibrium ( and ). This means that the latter SNPs were inherited as haplotypes. Further investigation reveals that none of these haplotypes associated with HCV susceptibility as shown in Figure 2(d).

As the patients recruited for this study already finished their interferon therapy protocols, we decided to test the association between the polymorphisms in the selected SNPs and the response to the interferon treatment. Hereafter, the HCV patients were classified based on their responses to interferon therapy into nonresponders (NR) and sustained virological responders (SVR) as shown in Table 4.

The patients with SVR showed a significantly lower AA genotype frequency of the NLRP3 rs1539019 compared to the nonresponder group (OR: 1.95). However, the frequencies of the CC genotype and C allele of the NLRP3 rs35829419 were significantly higher in the nonresponder group compared with the SVR group suggesting their association with the increased risk of lack of response to interferon treatment (OR: 2.78 and 2.73, respectively). On the other hand, the CA genotype and the A allele were associated with sustained response to IFN therapy (OR: 0.39, 0.37, respectively). No significant differences were found in the genotypes or alleles frequencies in the NLRP3 rs10754558 between the SVR and NR groups. For the IL-1β (rs1143629, rs1143634) and IL-18 (rs187238, rs1946518) polymorphisms, no significant association was found in terms of the response to IFN treatment.

When considering the CARD8 polymorphisms, the frequencies of the TT genotype and T allele of CARD8 rs2043211 were significantly higher in the NR group (OR: 2.64 and 1.54, respectively) making them associated with risk of failure of IFN treatment. However, the A allele of the same SNP was associated with more response to IFN treatment (OR: 0.65).

Upon classification of HCV patients based on the degree of liver fibrosis, we cannot find any significant difference in either genotypes or alleles frequencies of any of the studied SNPs (Table 5). This indicates lack of association between the studied inflammasome genes’ polymorphisms and the severity of fibrosis.

4. Discussion

The inflammation is known to be a major player in the progression of the HCV-associated liver diseases. Recently, many studies have shown that the NLRP3 inflammasome is implicated in recognizing HCV and eliciting the subsequent innate immune and inflammatory responses [5, 17, 30]. The genetic factors might affect the susceptibility, severity, and chronicity of HCV infection [31]. To the best of our knowledge, this is the first study investigating the association of polymorphisms in NLRP3 and CARD8 genes with the risk of HCV infection and its response to IFN therapy. The targeted SNPs in our study were chosen based on the previous publications [25, 26].

We studied three SNPs in the NLRP3 gene, and we found a statistically significant association between rs1539019 and rs35829419 polymorphisms and HCV infection. Our data demonstrated that the NLRP3 rs35829419 C>A polymorphism was associated with increased risk to HCV infection. Mounting efforts have been devoted for investigating the effect of the genetic variation in NLRP3 on the susceptibility to various diseases. Several studies demonstrated the association of NLRP3 rs35829419 C>A polymorphism with the increased risk of colorectal cancer, HIV-1 infection, rheumatoid arthritis, leprosy, and atopic dermatitis [3236].

Of note, a recent study found no correlation between the rs35829419 in NLRP3 and the renal damage caused by HCV-related cryoglobulinemia [37]. However, all the subjects recruited in this study were HCV patients (with or without renal damage) and no healthy control group was included. Thus, they did not test the association between the rs35829419 polymorphism and chronic HCV infection.

Interestingly, in our results the A allele of rs35829419 was more prevalent in the healthy control and the SVR groups indicating a protective effect of this allele against HCV infection. Recently, studies revealed that the NLRP3 Q705K (rs35829419) polymorphism is a gain-of-function change which results in more favourable innate immune response and this may provide a partial explanation for the protective effect of the A allele [38].

The association between the NLRP3 rs1539019 polymorphism and HCV infection is in the same line with a study demonstrating a strong association between the NLRP3 rs1539019 polymorphism and pneumoconiosis in Chinese coal workers [39]. Additionally, Dehghan et al. reported that rs1539019 is associated with the circulation fibrinogen levels and consequently the cardiovascular disease risk [40]. On the other hand, a study on a Japanese population reported no statistically significant association between the rs1539019 polymorphism and the risk of essential hypertension [41]. Notably, the NLRP3 rs1539019 is an intronic polymorphism whose function is not fully understood. However, several studies have reported that many transcription factors bind to intronic sites that may play a role in regulating the gene expression.

Interleukin-1 beta is a proinflammatory cytokine which is secreted by monocytes and macrophages and play a pivotal role in the innate immune response to infection [42]. There is mounting evidence that HCV patients have higher plasma IL-1β levels compared to healthy subjects emphasizing the implication of IL-1β in HCV pathogenesis [17]. However, the expression of IL-1β (mRNA levels) in the liver biopsies of chronic hepatitis C patients was found to be comparable to the controls [43]. In our study, the IL-1β rs1143629 revealed a significant association with HCV infection. In agreement with our findings, rs1143629 has been reported to be associated with malaria [44] and HIV-1 susceptibility [33]. Being an intronic SNP whose function is not fully understood does not exclude the importance of rs1143629 in influencing the IL-1β gene expression.

On the other hand, the IL-1β rs1143634, which was shown recently to be associated with advanced fibrosis in nonalcoholic steatohepatitis in Caucasians, did not reveal any association with the fibrosis degree in HCV patients in our study. rs1143634 is a synonymous SNP (F105F) in exon 5 of IL-1β whose function is not completely elucidated; yet, it was found to be associated with several diseases including HIV-1 [33] and periodontitis [45, 46]. However, no significant difference in rs1143634 genotyping or allele frequencies was detected in our study between HCV patients and healthy controls.

Interestingly, our results revealed a significant association between IL-18 rs1946518 and the susceptibility to HCV infection, but no significant association was observed for the IL-18 rs187238 polymorphism. Notably, both SNPs reside in the promoter region of the IL-18 gene and affect its expression and activity by disruption of transcription factor-binding sites [47, 48].

The proinflammatory cytokine IL-18 was reported to be involved in the HCV pathogenesis [49, 50]. In many studies, the IL-18 serum levels of HCV patients were significantly correlated with degree of liver injury and response to antiviral treatment [51, 52].

Our results showed that the T allele of IL-18 rs1946518 was associated with susceptibility to HCV; however, it was not associated with either response to treatment or fibrosis degree. Despite that there are many published studies which investigated the effect of IL-18 polymorphisms and the HCV outcome, the published data are very contradictory. A study conducted on Tunisian HCV patients reported a significant association between IL-18 polymorphism (rs1946518) and the severity of HCV [53]. A similar finding was reported in a study conducted on Indian patients [54]. However, a recent study stated that IL-18 rs1946518 polymorphism could not predict the outcome of the response of a combined interferon/ribavirin treatment in Egyptian HCV patients [55], which comes with agreement with our findings. The different findings among the studies could be explained by the differences in HCV genotypes and the ethnic groups.

The adaptor protein CARD8 is a crucial regulator of IL-1β production via inhibition of caspase-1 and/or inhibition of NFκB-mediated gene expression [56, 57]. To the best of our knowledge, our study is the first study investigating the effect of genetic variations in CARD8 on HCV pathogenesis and progression. Our results revealed no association between the rs6509365 or rs2043211 polymorphism with the susceptibility to HCV or the severity of the disease (fibrosis degree). However, only CARD8 rs2043211 was associated with the response to interferon therapy in our study. Previous studies have reported that CARD8 rs2043211 was found to be associated with susceptibility to gout [58], nodular melanoma [59], and decreased risk of ankylosing spondylitis [60]. Nonetheless, CARD8 rs2043211 could not predict the occurrence or progression of coronary artery diseases [61] or the susceptibility and outcome of meningitis [62]. The CARD8 rs2043211 (C10X) polymorphism leads to a nonfunctional truncated CARD8 and subsequently more active inflammasome and prolonged inflammation [38, 63].

We are aware that our study has some limitations including the small number of patients enrolled from only one population (Egyptian) which makes it hard to generalize the conclusions of our study. So, a consortium multicenters’ study with patients from different countries all over the world should be considered in order to draw a general conclusion on the correlation of the genetic variations in inflammasome genes and susceptibility to HCV. Also, this study did not include the spontaneously resolved patients (anti-HCV positive, negative serum viremia), which is another limitation. However, most of Egyptian patients with spontaneous clearance of HCV skip the proper diagnosis, and this made the recruitment of this group almost impossible. Moreover, we are also aware that some of the HCV treatment protocols have been shifted to the directly acting antiviral agents (DAA), which achieve a better sustained virologic response. Unfortunately, when we conducted our study, the DAA treatment protocol was just starting in Egypt with a very limited number of patients which makes it very hard to recruit them in our study. However, in the future we are planning to investigate the genetic variation in NLRP3 in correlation with the DAA response in HCV patients.

In conclusion, the present study reported, for the first time, that the genetic polymorphisms in NLRP3 (rs1539019 and rs35829419) are associated with hepatitis C susceptibility and response to treatment in Egyptian patients. Thus, our work highlights the importance of the genetic variations of innate immune factors for the predisposition to HCV infection and outcome. However, further large-scale studies are recommended to confirm our findings.

Abbreviations

NLRP3:NOD-like receptor, pyrin domain-containing 3
IL-1β:Interleukin-1β
ISGs:Interferon-stimulated genes
PEG-IFN:Pegylated interferon alpha
CARD8:Caspase recruitment domain-containing protein 8
SVR:Sustained virological response
HWE:Hardy-Weinberg equilibrium
SNP:Single-nucleotide polymorphism.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no competing interests to disclose.

Authors’ Contributions

Shady Z. K. Estfanous and Sahar A. Ali have equal contribution to the work.

Acknowledgments

The authors would like to express their gratitude to all the patients and healthy controls accepted to participate in the study. We would like to thank Dr. Sameera Azzat for the critical review of the manuscript. The study is supported by the Research Support and Technology Development Grant, project ID: 12636, from the Science and Technology Development Fund (STDF), Egypt.