Email updates

Keep up to date with the latest news and content from BMC Immunology and BioMed Central.

Open Access Research article

The impact of polymorphisms in STAT6 on treatment outcome in HCV infected Taiwanese Chinese

Yun-Ping Lim12, Yu-An Hsu3, Kun-Hsi Tsai4, Fuu-Jen Tsai5*, Cheng-Yuan Peng67, Wen-Ling Liao8, Dong-Zong Hung12, Ni Tien109, Chien-Yih Lin11* and Lei Wan115*

Author Affiliations

1 Department of Pharmacy, College of Pharmacy, China Medical University, Taichung 40402, Taiwan

2 Department of Emergency, Toxicology Center, China Medical University Hospital, Taichung 40447, Taiwan

3 Department of Life Science, National Tsing Hua University, Hsinchu 30013, Taiwan

4 Department of Emergency Medicine, Chi Mei Hospital, Liouying, Tainan, Taiwan

5 School of Chinese Medicine, China Medical University, No. 91, Hsueh-Shih Road, Taichung 40402, Taiwan

6 Department of Internal Medicine, China Medical University, Taichung 40402, Taiwan

7 Division of Digestive System and Gastroenterology, Department of Internal Medicine, China Medical University Hospital, Taichung 40447, Taiwan

8 Center for Personalized Medicine, China Medical University Hospital, Taichung 40447, Taiwan

9 Department of Laboratory Medicine, China Medical University Hospital, Taichung 40447, Taiwan

10 Department of Veterinary Medicine, National Chung Hsing University, Taichung 40227, Taiwan

11 Department of Biotechnology, Asia University, Taichung 41354, Taiwan

For all author emails, please log on.

BMC Immunology 2013, 14:21  doi:10.1186/1471-2172-14-21

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2172/14/21


Received:8 January 2013
Accepted:25 April 2013
Published:8 May 2013

© 2013 Lim et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Genetic polymorphisms observed in various disease states associated with sensitivity or resistance to specific treatments have been a robust area of investigation for decades, with the potential to allow clinicians to make evidence-based decisions on the appropriate course of treatment. This study aimed to evaluate whether genetic polymorphisms of the signal transducer and activator of transcription 6 gene (STAT6) could be associated with a sustained virological response (SVR) among patients infected with hepatitis C virus genotypes 1 and 2 (HCV-1 and HCV-2) who were treated with peginterferon plus ribavirin (PEG-IFNα-RBV). We analyzed the associations between SVR to PEG-IFNα-RBV therapy and 4 single nucleotide polymorphisms (SNPs) in STAT6. This study included Taiwanese Chinese patients infected with either HCV-1 (n = 265) or HCV-2 (n = 195) in the presence or absence of an SVR. Among the STAT6 SNPs examined, the dosage effect of the A allele and allele frequency in rs1059513 were inversely correlated with SVR in patients infected with HCV-1 (P = 0.0179 and P = 0.0235, respectively). This effect was not observed in patients infected with HCV-2. The GG, GGG, and GGGC STAT6 haplotypes comprising 2, 3, and 4 SNPs (rs1059513, rs703817, rs324015, and rs3024974) were found to be associated with SVR, and their presence may increase the probability of a successful treatment outcome in patients infected with HCV-1 (P = 0.0273, 0.0352, and 0.0368, respectively). Moreover, a multivariate logistic regression model for predicting an SVR revealed that the presence of the GGGC haplotype carriers mutually affected the outcome of PEG-IFNα-RBV treatment. The presence of STAT6 SNPs and the association with SVR demonstrated that STAT6 polymorphisms might influence the therapeutic outcomes of patients infected with HCV-1 under standard-of-care (SOC) treatment.

Keywords:
Hepatitis C virus; Standard of care treatment; Sustained virological response; Signal transducer and activator of transcription 6

Background

Hepatitis C is a common infectious disease caused by hepatitis C virus (HCV). Infection with HCV is the most common cause of chronic liver disease, hepatoma, and liver transplantation [1-3]. Typically, HCV spreads through exposure to infected blood, especially via intravenous drug abuse, which is a major route of transmission. In addition to blood-borne exposure, HCV may also spread through sexual contact and vertical transmission, albeit at lower risk [3]. It is estimated that approximately 170 million people worldwide are infected with HCV [4,5].

The standard-of-care (SOC) treatment for chronic hepatitis C includes a combination of pegylated interferon alpha (PEG-IFNα) and ribavirin (RBV) [6,7]. The effectiveness of 2 types of peginterferon, alfa-2a and alfa-2b, both differing in their pharmacokinetic and chemical properties, has been determined [8,9]. Peginterferon alfa-2a plus ribavirin is more effective in treating chronic hepatitis C than interferon alfa-2b plus ribavirin or peginterferon alfa-2a alone [8]. In general, successful treatment is defined as the persistent absence of serum HCV RNA for 6 months or longer post therapy, known as a sustained virological response (SVR) [10].

IFNα exert multifunctional biological activities and consists of highly homologous proteins mediated by binding to the receptor chains, IFNαR1 or IFNαR2 [11]. Binding of IFNα to its receptor leads to its activation via two receptor-associated tyrosine kinases (Tyk2) and Janus kinases 1 and 2 (JAK kinases 1 and 2), which in turn triggers the activation of cytoplasmic factors, namely, the STAT (signal transducers and activators of transcription) family of proteins [11]. STAT6 is a cytosolic protein, the tyrosine phosphorylation of which is catalyzed by JAK kinases 1 and 2 [11]. Subsequently, phosphorylated STAT6 proteins dimerize via their SH2 domains and translocate into the nucleus, where they bind to response elements of target genes and regulate cytokine production [11,12]. STAT6 induces two forms of cytokines, IL-4 and IL-13. Formation of STAT homo/heterodimer in response to IFNα activation may initiate transcription of several genes, including genes encoding the transcription factor IFN regulatory factor 1 (IRF-1) and IFN-stimulated gene factor 3 (ISGF3). The ISGF3 complex may trigger several subsets of genes termed IFN-stimulated response elements (ISREs) by binding to different parts of cis-acting elements [12]. The JAK/STAT signaling pathway is activated when the cytokine interleukins 4 (IL-4), IL-12, and/or type I interferon (IFN) bind to their respective receptors. Since STAT6 protein is essential for responsiveness to IL-4 and IL-12, STAT6-deficient mice do not produce significant numbers of type 2 CD4+ T cells, thereby blocking the differentiation of these cells [13,14]. Furthermore, multiple studies have demonstrated that STAT6 knockout mice are resistant to tumor growth and spontaneous metastatic disease [15]. Therefore, IFN-mediated activation of STAT6 occurs via production of secreted IL-1Ra during inflammatory responses [16]. STAT6 activation has been reported to play an important role in the regulation of genes responsible for their antitumor or antiviral effects.

Previous studies have demonstrated that the SVR rates for patients infected with HCV genotype 1 (HCV-1) are significantly lower than the rates for those who are infected with HCV genotype 2 (HCV-2). In order to validate the potentially useful prognostic biomarkers that are able to predict SVR during early treatment, previous studies have used real-time PCR to estimate the expression profiles of 68 messenger RNAs isolated from HCV-1 infected patients, as well as their correlation with SVR [17]. Investigators found that STAT6 expression levels are highly correlated with SVR. It has been suggested that, in addition to inducing important molecular pathways, the expression of several early genes during HCV therapy may also influence the capacity of residual SVR.

In the present study, we investigated the association between single nucleotide polymorphisms (SNPs) in the STAT6 gene and their susceptibility to SVR in Chinese patients in Taiwan receiving PEG-IFNα-RBV treatment. Our results support the hypothesis that STAT6 is a potential candidate gene for predicting therapeutic outcomes of HCV-1 infected patients.

Methods

Patients

In the present study, 265 patients infected with HCV-1 and 195 patients infected with HCV-2 at the China Medical University Hospital, Taichung, Taiwan were enrolled and actively observed. All participating subjects provided informed consent, and the study protocol was approved by the chairman of the Ethics Committee of the China Medical University Hospital in accordance with the guidelines of the Declaration of Helsinki. Diagnosis of HCV infection was based on persistent elevation of serum transaminase levels for at least 6 months and serum anti-HCV positivity, coupled with the periodic detection of serum HCV RNA. Patients positive for hepatitis B surface antigen, antibodies, and human immunodeficiency virus 1 and 2 were excluded in this study. Patients received weekly injections of PEG-IFNα (1.5 μg/kg body weight) plus body weight–adapted doses of RBV (800 mg/day for < 60 kg; 1000 mg/day for 60–75 kg; 1200 mg/day for > 75 kg) by oral administration for 48 weeks (HCV-1) or 24 weeks (HCV-2).

SNP selection

Selection of representative STAT6 SNPs was based on SNP genotype information, downloaded in December 2008 from the HapMap Chinese Han in Beijing (CHB) + JPT population. HapMap genotypes were further analyzed via Haploview software (version 4.2; Broad Institute, Cambridge, MA, USA). Tag SNPs were selected by using the Tagger function according to the following criteria: (1) a minor allele frequency (MAF) in the HapMap CHB + JPT population of > 0.10; and (2) a genotyping score of ≥ 0.6 (Illumina, Inc., San Diego, CA), as recommended by the manufacturer, to achieve a successful genotyping rate. Four SNPs in the STAT6 gene met the above criteria and were selected, including rs1059513 (S1; A/G at 3′ UTR), rs703817 (S2; A/G at 3′ UTR), rs324015 (S3; A/G at 3′ UTR), and rs3024974 (S4; C/T at boundary of intron 17).

HCV genotyping and RNA measurements

Genotyping of HCV according to the classification of Simmonds et al.[18] was performed by a reverse hybridization assay (INNO LiPA HCV-II; Innogenetics, Ghent, Belgium). The virological response was assessed by qualitative HCV RNA assays with a lower sensitivity of 30–50 IU/ml (HCV Amplicor™ 2.0, Roche Diagnostics, Branchburg, NJ, USA). According to the qualitative HCV RNA results (reported in copies/mL), patients were defined as (1) sustained virological responders (HCV RNA undetectable at week 24 posttreatment) and (2) non-sustained virological responder (HCV RNA detected at week 24 posttreatment). Thus, the subjects were classified as either SVR (+) or SVR (−), respectively.

Genomic DNA extraction and genotyping

Genomic DNA was isolated from the peripheral blood of all participants by using a commercial kit (Genomic DNA kit; QIAGEN, Valencia, CA, USA) according to the manufacturer’s instructions. Four STAT6 polymorphisms were detected by an allele-specific extension method and ligation assay kit (Illumina, San Diego, CA, USA), in accordance with the manufacturer’s instructions.

Statistical analysis

Gender, age, body mass index (BMI), and viral load differences between SVR (+) and SVR (−) groups were estimated using the Mann–Whitney U test. The differences between genotypes and each of the above mentioned parameters were estimated using the Kruskal–Wallis test. The association of each SNP with SVR was assessed using the chi-square (χ2) test or Fisher’s exact test. Genotype and allele frequencies in SVR (+) and SVR (−) patients were compared, and odds ratios (ORs) with a 95% confidence interval (CI) were estimated by applying an unconditional logistic regression model. Haplotype analysis using a sliding window mode and Hardy–Weinberg equilibrium (HWE) analysis were performed using the PLINK program (version 1.07; http://pngu.mgh.harvard.edu/purcell/plink/ webcite) [19]. The distribution of each haplotype frequency was identified using Phase version 2.1 (University of Chicago, Chicago, IL, USA), a computational tool based on Bayesian methods [20]. Multivariate logistic regression analysis was used to predict the independent variables that correlated with a therapeutic response. An OR at a 95% CI was employed for this analysis. Statistical tests were performed using SPSS software (version 20.0 for Windows; SPSS Inc., Chicago, IL, USA). Statistical significance was defined as P < 0.05 and was two-sided. Screening and construction of linkage disequilibrium (LD) plots were performed using Haploview software (version 4.2).

Results

Patient characteristics

In this study, 460 patients infected with HCV-1 (n = 265) and HCV-2 (n = 195) were successfully enrolled. Patients infected with HCV-2 achieved SVR at higher rates than patients infected with HCV-1 (91.3% vs. 61.5%). The basic characteristics and clinical information of all patients infected with HCV-1 and HCV-2 are displayed in Table 1. The information is further stratified by HCV genotype. We observed several significant differences between SVR (+) and SVR (−) groups, including their age (starting with date of entry), body mass index (BMI), platelets and viral load in HCV-1–infected patients (P = 0.002, P = 0.026, P = 0.002, and P = 0.001, respectively). Some effects were not observed significantly in HCV-2 infected populations. However, statistically significant differences in the aforementioned variables were not observed among the 4 SNPs.

Table 1. Characteristics of the HCV genotype 1 and 2 infected patients receiving PEG-IFNα-RBV therapy

The chromosome position, Hardy–Weinberg equilibrium (HWE), and MAF in individuals bearing the aforementioned 4 STAT6 SNPs genotypes examined in this study are displayed in Table 2. The SNPs were in accordance with HWE (P > 0.05).

Table 2. Four single nucleotide polymorphisms in the STAT6 gene in 265 HCV-1 and 195 HCV-2 infected patients receiving PEG-IFNα-RBV therapy with or without a SVR in a Chinese population in Taiwan

Association of tagging SNPs of STAT6 with therapeutic response and SVR

The genotype frequencies and allelic dose-dependent association of each SNP with responsiveness to PEG-IFNα-RBV therapy are displayed in Table 3. In genotype association tests, only genotypes from rs1059513 were significantly associated with SVR in HCV-1 infected patients (P = 0.0179). When the A/A and A/G genotypes were compared, a decreased crude OR of 0.40 (95% CI, 0.18–0.87; P = 0.0179) for SVR (+) versus SVR (−) groups was observed. However, no genotype was significantly associated with SVR in the HCV-2 infected population.

Table 3. Genotype frequencies of STAT6 single nucleotide polymorphisms (SNPs) in HCV-1 and HCV-2 infected patients receiving PEG-IFNα-RBV therapy with and without a SVR in a Chinese population in Taiwan

The distributions of the aforementioned 4 STAT6 SNPs allele frequencies between SVR (+) and SVR (−) groups in patients infected with HCV-1 and HCV-2 are displayed in Table 4. Compared with SVR (+) individuals, patients harboring the A allele of SNP rs1059513 were at an increased risk for non-responsiveness to SOC treatment (P = 0.0235; OR = 0.42, 95% CI, 0.20–0.91), with a significance observed exclusively in patients infected with HCV-1. Thus, HCV-1 infected individuals with STAT6 rs1059513 A/A genotype and carrying the A allele may be at increased risk for non-responsiveness to PEG-IFNα-RBV treatment.

Table 4. Allele frequencies of STAT6 single nucleotide polymorphisms in HCV-1 and HCV-2 infected patients receiving PEG-IFNα-RBV therapy with and without a SVR in a Chinese population in Taiwan

Frequencies of the STAT6 haplotypes and the association with SVR

We performed a comparative analysis of haplotype frequencies and treatment responses by adjustment for gender as a covariate by using the PLINK program. The results of the overall global test and haplotype frequencies are displayed in Table 5. The analysis was introduced by the sliding window mode to examine the potential sizes of all haplotypes (number SNPs per haplotype). As a result, there were 10 sliding windows, 1 of which (the omnibus test) was significantly associated with SVR in HCV-1 infected patients (P < 0.05).

Table 5. Details of sex-adjusted haplotype frequency analysis for 2-SNP, 3-SNP, and 4-SNP windows showing the most significant results among all possible sliding windows in HCV-1 infected patients

In HCV-1 infected populations possessing haplotype GG, the window S1-S2 (composed of rs1059513 and rs703817), gave the most significant P value (P = 0.0273; OR = 2.31). Haplotype-specific analyses showed that the GGG haplotype (S1-S2-S3, composed of rs1059513, rs703817, and rs324015) may increase the rate of SVR (P = 0.0352; OR = 2.24) in the SVR (+) groups compared with the SVR (−) groups. The window S1-S2-S3-S4 with the GGGC haplotype (composed of rs1059513, rs703817, rs324015, and rs3024974) was significantly associated with a higher rate of SVR in HCV-1 infected patients (OR = 2.22; P = 0.0368). Furthermore, the results demonstrated that HCV-1 infected patients with therapeutic response bearing the haplotypes GG, GGG, and GGGC, appeared more frequently in SVR (+) groups than in SVR (−) groups. However, the haplotypes did not play a significant role in HCV-2 infected patients (data not shown). Thus, in the HCV-1 infected individuals, haplotype-specific analyses showed that the aforementioned haplotypes are associated with an increase success rate of SVR and may play a role in the therapeutic outcomes of PEG-IFNα-RBV treatment in HCV-1 infected patients. Table 6 shows that 9.5% of 163 SVR (+) HCV-1 infected patients possessed the GGGC haplotype. The GGGC haplotype carriers achieved a significantly higher SVR rate compared with non-GGGC carriers (95% CI, 1.06–4.89, OR = 2.28; P = 0.0306). These findings confirmed the results of our analyses of genotypes, allele frequencies, and haplotypes.

Table 6. Distribution of haplotype frequencies from SNPs of STAT6 gene in 265 HCV-1 infected populations and their associations with SVR

We performed LD analysis, and the results were viewed using Haploview software (version 4.2). Our results indicated that a low degree of pairwise LD among these SNPs in all HCV-1 infected populations in the presence or absence of SVR. The graphical summary of LD (r2 values) among tested 4 STAT6 SNPs at different loci is displayed in Figure 1. Four tag SNPs were selected and designated in sequential order.

thumbnailFigure 1. Linkage disequilibrium plot of the STAT6 gene for the analyzed SNPs. The figure shows the output of Haploview (version 4.2) linkage disequilibrium plot where each square (with r2 values written within the box correspond to r2 values × 100 as linkage disequilibrium measure range) represents a pairwise linkage disequilibrium relationship between the two SNPs. Darkest colored squares indicate high linkage disequilibrium (r2 = 1); mid colored squares indicate r2 values between 0 and 1; lightest colored squares indicate low linkage disequilibrium (r2 = 0). Figure depicts the linkage disequilibrium pattern in (A) all HCV-1 infected patients; (B) HCV-1 infected PEG-IFNα-RBV therapy with SVR (+) patients; (C) HCV-1 infected PEG-IFNα-RBV therapy with SVR (−) patients.

Multivariate logistic regression analysis

We subsequently determined whether independent factors affected the outcome of PEG-IFNα-RBV therapy. To this end, we performed multivariate logistic regression analyses with respect to GGGC/non-GGGC haplotypes, age, BMI, and HCV RNA levels (low/high). This analysis revealed that the GGGC haplotype carrier acts as an independent predictor of the virological response outcome of PEG-IFNα-RBV therapy (P = 0.043; OR = 2.32, 95% CI, 1.03–5.23) (Table 7). Furthermore, age and HCV viral load may also act together to influence the outcome of PEG-IFNα-RBV treatment (P = 0.012; OR = 0.97, 95% CI, 0.94–1.00; and P = 0.050; OR = 8.11, 95% CI, 0.98–67.32, respectively).

Table 7. Predictive factors associated independently with the SVR to PEG-IFNα-RBV therapy in HCV-1 infected populations by multivariate logistic regression analysis

Discussion

The present study examined the genetic association of HCV infection under PEG-IFNα-RBV therapy on a responsive gene, STAT6, and responsiveness of treatment. We confirmed that this gene, located on chromosome 12, is a novel candidate susceptibility gene. The results showed that rs1059513 at the 3′ UTR of STAT6 may be associated with SVR rates under SOC treatment. Therapeutic outcome association studies showed that the presence of the GGGC haplotype was significantly associated with a higher rate of SVR (OR = 1.06–4.89) in HCV-1 infected but not HCV-2 infected populations. By using multivariate logistic regression analyses, we observed that the GGGC haplotype, age of entry, and viral load may act interactively to influence the therapeutic outcomes. This result suggests that the SVR rates of PEG-IFNα-RBV treatment could, at least partially, depend on STAT6 activation.

Thus far, the standard treatment regimen for chronic hepatitis C infection has been a combination of PEG-IFNα and RBV [8]. About 50% of treatment-naïve populations will typically achieve SVR by using this regimen [8]. The SVR rate is markedly higher in HCV-2/3 infected patients than in HCV-1 infected patients (~ 85% versus 40%) [8]. In order to identify the biomarkers that distinguish the profiles of HCV therapeutic responsiveness, Younossi and colleagues [17] collected RNA transcripts from peripheral blood mononuclear cells (PBMCs) of 68 HCV-infected patients under PEG-IFNα-RBV treatment. After preparing expression profiles of patients under this treatment, the investigators observed that genetic expression of STAT6 and cytokine signaling 1 are putative predictive markers that correlate with SVR [17]. In support of this finding, we also observed a correlation between STAT6 genetic polymorphisms and SVR.

During viral infections, innate immunity remains the most important first line of defense in the body. Viral nucleic acids are detected by endosomal Toll-like receptors (TLRs) and cytoplasmic RIG-I–like receptors (RLRs), which subsequently induce the production of type 1 IFN to inhibit viral replication [22]. Consequently, many cytokines (including type 1 IFNs) trigger the IFN-stimulated gene expression cascade through the canonical JAK-STAT signaling pathway [23]. Previous studies have shown that STAT6 activation is ubiquitously detected during viral infections, concluding that STAT6 may play a fundamental role in the defense mechanism against viral infections [24]. Viral load, age of entry into the study, and BMI were not significantly different among the STAT6 genotypes of our enrolled populations. Thus, it is unlikely that specific STAT6 genotypes predispose individuals to infection with HCV or contribute to the process of spontaneous viral elimination.

Studies using a microarray platform to examine the differential gene expression profiles of IL-4–stimulated B cells between wild-type and Stat6−/− mice were recently performed [25]. Seventy known genes were found to be differentially expressed with significant discrepancies between the 2 mouse groups, which were all negatively or positively regulated by STAT6 via a palindromic consensus sequence TTC(N)2-4GAA [26]. In addition, STAT6 remains a critical determinant in the polarization and differentiation process of T-helper type 2 (Th2) cells in a Stat6 knockout mouse model [27]. The expression levels of Th2-related cytokines (including IL-4, IL-5, and IL-13) were diminished in this mouse model [27]. Moreover, STAT6 may trigger a Th1 response by enhancing the production of IL-12 through the inhibition of IL-10 production in dendritic cells [28]. Therefore, STAT6 is also a critical factor in the homeostasis of inflammatory and hypersensitivity immune responses other than those associated with viral immunity.

Several reports have provided strong evidence that patients infected with HCV-1 will have approximately a 50% probability (in Caucasians) or 80% probability (in African Americans) of a poor response toward PEG-IFNα-RBV treatment [8]. Although viral clearance rates have been strongly associated with clinical features of HCV, for example: gender, age < 40 years, low HCV RNA levels pre-treatment, in the absence of liver cirrhosis, and HCV genotypes 2/3 [29]. Table 7 lists 11 patients with low viremia. Fried et al.[8] concluded that HCV RNA level is one of the main factors affecting treatment outcome; therefore, we performed multivariate analysis to check whether HCV RNA level is an independent factor of treatment outcome. We first determined the HCV RNA baseline without categorizing the patients; however, in this analysis, P = 0.182 and OR = 1.00. This result led us to conclude that it is not the best model for prediction of treatment outcome. Therefore, we categorized the HCV RNA level according to the description of Fujimoto et al.[30], with low HCV RNA level <100 KIU/ml by Amplicore-monitor assay and performed statistical analysis, and the results for this analysis are more obviously (P = 0.050, OR = 8.11).

Many investigators continue to monitor the host genetic factors that may relate to clinical outcomes to provide custom-made therapy for HCV infection. In our study, the overall achievement response rates of SVR were less than 62% and 92% in the HCV-1 and HCV-2 populations, respectively. Thus, a reliable prediction of a subject as a potential non-viral responder (NVR) at the start of treatment would allow the clinician to make appropriate, evidence-based decisions on the compatibility of future treatment, thereby reducing side effects and/or the cost of treatment.

The contributions of host genetic factors to HCV-2/3 clearance are relatively small compared with HCV-1 clearance, because the former is more likely to be eliminated by SOC therapy [8]. In the present study, we observed an association between the STAT6 haplotype GGGC and the rate of SVR, although its statistical significance in the HCV-1 infected populations based on estimation by multivariate logistic regression was marginal. The frequency of the GGGC allele in our population was small, whether the frequency of this haplotype is significant in other ethnic groups remains unknown. We propose that this STAT6 haplotype GGGC might play a role in the outcome of PEG-IFNα-RBV therapy (adjusted r2 = 11.2%). The linkage among the STAT6 polymorphisms and SVR should be confirmed in future studies with larger sample sizes. Moreover, the significance of the genetic effects of STAT6 on other racial and ethnic groups remains to be elucidated.

The mechanisms of action of the aforementioned 4 SNPs located within the 3′ UTR and intronic regions of STAT6 gene are hampered by the lack of published functional studies on these polymorphisms and warrants further investigation at the molecular level.

Conclusion

This study provides evidence that the STAT6 genetic polymorphisms are related to the susceptibility of HCV-infected patients, especially HCV-1 infected patients to SVR. Furthermore, an outcome measure of haplotype GGGC in the clinical response was identified. Thus, STAT6 may play an important role in a broad range of antiviral activities under PEG-IFNα-RBV treatment.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

YPL and LW designed and carried out most of the study. YPL and YAH wrote the manuscript and YPL performed data analysis. KHT, CYP, CYC, WLL, DZH and NT participated in clinical data and information collection. FJT, CYL, and LW conceived and supervised the project and reviewed the manuscript. All authors contributed to and approved the final manuscript by providing constructive suggestions.

Acknowledgments

This study was supported by National Science Council, Executive Yuan, Taiwan, R.O.C. (NSC101-2320-B-039-038 and NSC101-2320-B-039-007-MY3), China Medical University Hospital, Taichung, Taiwan (DMR-99-093), China Medical University, Taichung, Taiwan (CMU98-ASIA-04), and in part by Taiwan Department of Health Clinical Trial and Research Center of Excellence (DOH102-TD-B-111-004).

References

  1. Alter MJ: Epidemiology of hepatitis C virus infection.

    World J Gastroenterol 2007, 13:2436-2441. PubMed Abstract | Publisher Full Text OpenURL

  2. Bowden DS, Berzsenyi MD: Chronic hepatitis C virus infection: genotyping and its clinical role.

    Future Microbiol 2006, 1:103-112. PubMed Abstract | Publisher Full Text OpenURL

  3. Bellentani S, Miglioli L, Bedogni G, Croce LS, Tiribelli C: Epidemiology of hepatitis C virus infection.

    Minerva Gastroenterol Dietol 2005, 51:15-29. PubMed Abstract OpenURL

  4. Shepard CW, Finelli L, Alter MJ: Global epidemiology of hepatitis C virus infection.

    Lancet Infect Dis 2005, 5:558-567. PubMed Abstract | Publisher Full Text OpenURL

  5. Pellicano R, Mladenova I, Dimitrova SM, Bruno CM, Sciacca C, Rizzetto : The epidemiology of hepatitis C virus infection. An update for clinicians.

    Minerva Gastroenterol Dietol 2004, 50:1-7. PubMed Abstract OpenURL

  6. Ascione A, De Luca M, Tartaglione MT, Lampasi F, Di Costanzo GG, Lanza AG, Picciotto FP, Marino-Marsilia G, Fontanella L, Leandro G: Peginterferon alfa-2a plus ribavirin is more effective than peginterferon alfa-2b plus ribavirin for treating chronic hepatitis C virus infection.

    Gastroenterology 2010, 138:116-122. PubMed Abstract | Publisher Full Text OpenURL

  7. Bosques-Padilla F, Trejo-Estrada R, Campollo-Rivas O, Cortez-Hernandez C, Dehesa-Violante M, Maldonado-Garza H, Perez-Gomez R, Cabrera-Valdespino A: Peginterferon alfa-2a plus ribavirin for treating chronic hepatitis C virus infection: analysis of Mexican patients included in a multicenter international clinical trial.

    Ann Hepatol 2003, 2:135-139. PubMed Abstract OpenURL

  8. Fried MW, Shiffman ML, Reddy KR, Smith C, Marinos G, Gonçales FL Jr, Häussinger D, Diago M, Carosi G, Dhumeaux D, Craxi A, Lin A, Hoffman J, Yu J: Peginterferon alfa-2a plus ribavirin for chronic hepatitis C virus infection.

    N Engl J Med 2002, 347:975-982. PubMed Abstract | Publisher Full Text OpenURL

  9. Manns MP, McHutchison JG, Gordon SC, Rustgi VK, Shiffman M, Reindollar R, Goodman ZD, Koury K, Ling M, Albrecht JK: Peginterferon alfa-2b plus ribavirin compared with interferon alfa-2b plus ribavirin for initial treatment of chronic hepatitis C: a randomised trial.

    Lancet 2001, 358:958-965. PubMed Abstract | Publisher Full Text OpenURL

  10. Sokal EM, Bourgois A, Stephenne X, Silveira T, Porta G, Gardovska D, Fischler B, Kelly D: Peginterferon alfa-2a plus ribavirin for chronic hepatitis C virus infection in children and adolescents.

    J Hepatol 2010, 52:827-831. PubMed Abstract | Publisher Full Text OpenURL

  11. Wurster AL, Tanaka T, Grusby MJ: The biology of Stat4 and Stat6.

    Oncogene 2000, 19:2577-2584. PubMed Abstract | Publisher Full Text OpenURL

  12. Dickensheets HL, Venkataraman C, Schindler U, Donnelly RP: Interferons inhibit activation of STAT6 by interleukin 4 in human monocytes by inducing SOCS-1 gene expression.

    Proc Natl Acad Sci USA 1999, 96:10800-10805. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  13. Galka E, Thompson JL, Zhang WJ, Poritz LS, Koltun WA: Stat6 (null phenotype) human lymphocytes exhibit increased apoptosis.

    J Surg Res 2004, 122:14-20. PubMed Abstract | Publisher Full Text OpenURL

  14. Zhang S, Bernard D, Khan WI, Kaplan MH, Bramson JL, Wan Y: CD4+ T-cell-mediated anti-tumor immunity can be uncoupled from autoimmunity via the STAT4/STAT6 signaling axis.

    Eur J Immunol 2009, 39:1252-1259. PubMed Abstract | Publisher Full Text OpenURL

  15. Ostrand-Rosenberg S, Sinha P, Clements V, Dissanayake SI, Miller S, Davis C, Danna E: Signal transducer and activator of transcription 6 (Stat6) and CD1: inhibitors of immunosurveillance against primary tumors and metastatic disease.

    Cancer Immunol Immunother 2004, 53:86-91. PubMed Abstract | Publisher Full Text OpenURL

  16. Wan L, Lin CW, Lin YJ, Sheu JJC, Chen BH, Liao CC, Tsai Y, Lin WY, Lai CH, Tsai FJ: Type I IFN induced IL1-Ra expression in hepatocytes is mediated by activating STAT6 through the formation of STAT2: STAT6 heterodimer.

    J Cell Mol Med 2008, 12:876-888. PubMed Abstract | Publisher Full Text OpenURL

  17. Younossi ZM, Baranova A, Afendy A, Collantes R, Stepanova M, Manyam G, Bakshi A, Sigua CL, Chan JP, Iverson AA, Santini CD, Chang SY: Early gene expression profiles of patients with chronic hepatitis C treated with pegylated interferon-alfa and ribavirin.

    Hepatology 2009, 49:763-774. PubMed Abstract | Publisher Full Text OpenURL

  18. Simmonds P, Holmes EC, Cha TA, Chan SW, McOmish F, Irvine B, Beall E, Yap PL, Kolberg J, Urdea MS: Classification of hepatitis C virus into six major genotypes and a series of subtypes by phylogenetic analysis of the NS-5 region.

    J Gen Virol 1993, 74:2391-2399. PubMed Abstract | Publisher Full Text OpenURL

  19. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, Maller J, Sklar P, de Bakker PI, Daly MJ, Sham PC: PLINK: a tool set for whole-genome association and population-based linkage analyses.

    Am J Hum Genet 2007, 81:559-575. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. Barrett JC, Fry B, Maller J, Daly MJ: Haploview: analysis and visualization of LD and haplotype maps.

    Bioinformatics 2005, 21:263-265. PubMed Abstract | Publisher Full Text OpenURL

  21. The French METAVIR Cooperative Study Group: Intraobserver and interobserver variations in liver biopsy interpretation in patients with chronic hepatitis C.

    Hepatology 1994, 20:15-20. PubMed Abstract OpenURL

  22. Kato H, Sato S, Yoneyama M, Yamamoto M, Uematsu S, Matsui K, Tsujimura T, Takeda K, Fujita T, Takeuchi O, Akira S: Cell type-specific involvement of RIG-I in antiviral response.

    Immunity 2005, 23:19-28. PubMed Abstract | Publisher Full Text OpenURL

  23. Gale M Jr, Foy EM: Evasion of intracellular host defence by hepatitis C virus.

    Nature 2005, 436:939-945. PubMed Abstract | Publisher Full Text OpenURL

  24. Chen H, Sun H, You F, Sun W, Zhou X, Chen L, Yang J, Wang Y, Tang H, Guan Y, Xia W, Gu J, Ishikawa H, Gutman D, Barber G, Qin Z, Jiang Z: Activation of STAT6 by STING is critical for antiviral innate immunity.

    Cell 2011, 147:436-446. PubMed Abstract | Publisher Full Text OpenURL

  25. Schroder AJ, Pavlidis P, Arimura A, Capece D, Rothman PB: Cutting edge: STAT6 serves as a positive and negative regulator of gene expression in IL-4-stimulated B lymphocytes.

    J Immunol 2002, 168:996-1000. PubMed Abstract | Publisher Full Text OpenURL

  26. Schindler U, Wu P, Rothe M, Brasseur M, McKnight SL: Components of a Stat recognition code: evidence for two layers of molecular selectivity.

    Immunity 1995, 2:689-697. PubMed Abstract | Publisher Full Text OpenURL

  27. Kaplan MH, Schindler U, Smiley ST, Grusby MJ: Stat6 is required for mediating responses to IL-4 and for development of Th2 cells.

    Immunity 1996, 4:313-319. PubMed Abstract | Publisher Full Text OpenURL

  28. Yao Y, Li W, Kaplan MH, Chang CH: Interleukin (IL)-4 inhibits IL-10 to promote IL-12 production by dendritic cells.

    J Exp Med 2005, 201:1899-1903. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Gao B, Hong F, Radaeva S: Host factors and failure of interferon-alpha treatment in hepatitis C virus.

    Hepatology 2004, 39:880-890. PubMed Abstract | Publisher Full Text OpenURL

  30. Fujimoto Y, Ochi H, Maekawa T, Abe H, Hayes CN, Kumada H, Nakamura Y, Chayama K: A single nucleotide polymorphism in activated Cdc42 associated tyrosine kinase 1 influences the interferon therapy in hepatitis C patients.

    J Hepatol 2011, 54:629-639. PubMed Abstract | Publisher Full Text OpenURL