The main findings of this study were as follows: 1) Arg70Gln was the most common substitution of core aa 70, which was significantly associated with an increased risk of cirrhosis; 2) Arg70Gln was strongly correlated with the HCV subtype 1b; 3) Core aa 70 substitutions were more common in the rs12979860 TT genotype and the rs8099917 non-TT genotypes, and a mild correlation was found between core aa 70 substitution and polymorphisms near the IFNL3 gene; 4) Lower HCV RNA levels were significantly associated with the rs12979860 TT and the rs8099917 GG genotypes as well as HCV subtype 3a; 5) the rs8099917 GG and rs12979860 CC genotypes were observed more frequently in patients infected with HCV subtypes 1b and 3a, respectively.
To the best of our knowledge, HCV subtype 1a has been previously demonstrated as the most common in Iran (
4,
6,
12,
37), and this prevalence was confirmed in this study. Few studies (Isfahan province (
14) and Yazd province (
38), central Iran; and Mashhad, northeast Iran (
7)) have reported HCV subtype 3a as the most common (61.2%, 50.3%, and 40%, respectively). The different age of the studied patients could justify these differences, as it has been shown that HCV subtype 3a is more common in patients younger than 40 years (41.5%) (
12). Studies from other countries have also confirmed an association between a young age and HCV subtype 3a (
39-
42). One study on Iranian HCV-infected thalassemic patients in central Iran also revealed that HCV subtype 1a is the most common genotype (
15). On the other hand, a significantly higher frequency of HCV subtype 1b was identified in hemophilic patients, which was compatible with the findings of a previous study (
18).
The results of the present study regarding core aa 70 substitutions in Iranian HCV-infected patients demonstrated an association of Arg70Gln with an increased risk of cirrhosis; interestingly, previous evidence supports the observations of this study (
43). Based on these findings, Arg70Gln could be implied as a plausible predictor for cirrhosis and severe liver disease, especially when associated with HCV subtype 1b (
44). It is worthwhile to note that from a molecular perspective, a probable role has been suggested for core protein in HCV-related cirrhosis via the up-regulation of vascular endothelial growth factor expression (
45,
46). Interestingly, some rare substitutions at core aa 70 have been observed; 18% of core aa 70 substitutions included histidine, proline, or leucine.
In this study, it was revealed that Arg70Gln was significantly associated with HCV subtype 1b. While most of the previous studies on the therapeutic impacts of HCV subtype 1b have been conducted in Japan, a few other investigations have also reported high frequencies of Arg70Gln in patients infected with HCV subtype 1b in different populations, such as in Venezuela, the United States, and China (
47). The Arg70Gln frequency in HCV subtype 1b was reported within a wide range, from 10% in China to 79% in Venezuela, and is compatible with the reported frequency in this study (54.5%) (
47). Jaspe et al. (
47) did not report any frequencies of Arg70Gln in HCV subtypes 1a and 3a that were higher than 3.7% and 11%, respectively, which were similar to the findings of the present study. An association of Arg70Gln with HCV subtype 1b has been reported along with a higher risk of poor treatment response (
44,
48), lipid accumulation and steatosis, and an increased risk of HCC and insulin resistance (
46). Although a higher rate of HCC is reported in patients with eradicated HCV subtype 1b infection associated with core Arg70Gln (
49), early therapeutic intervention with an appropriate regimen to confront the infection by HCV subtype 1b associated with core Arg70Gln may improve the outcome and prevent long-term sequelae. Future studies are required to clarify these impacts. The high frequency of core aa 70 substitution that was found in hemophilic patients could also be due to the increased frequency of HCV subtype 1b within hemophilic patients and the association between HCV subtype 1b and core aa 70 substitution (
47).
The genetic polymorphisms located near the
IFNL3 gene and core aa 70 substitutions of HCV have been demonstrated to be predictive of the treatment response and prognosis (
2,
19,
20,
25,
50). Based on a recent meta-analysis, the rs12979860 CC and rs8099917 TT genotypes have been associated with a higher probability of sustained virological response (SVR) (
51). These genotypes have been also correlated with higher HCV RNA levels in the absence of treatment (
52-
54). Along the same lines, the rs12979860 TT and rs8099917 GG genotypes were correlated with lower HCV RNA levels in this study of HCV-infected patients. Lower HCV RNA levels in patients with the rs12979860 TT genotype were also previously reported by Haj-Sheykholeslami et al. (
25). This finding could be interpreted in light of the higher expression of interferon alpha and beta genes and also interferon-stimulated genes in patients with the rs12979860 TT genotype (
55-
57).
Concerning the genotype frequencies of polymorphisms near the
IFNL3 gene in Iranian patients, the rs12979860 CC, CT, and TT genotype frequencies were reported to be as high as 38%, 48.8%, and 13.2%, respectively, by Sharafi et al. (
24). Additionally, the frequencies of the rs8099917 TT, GT, and GG genotypes were 58.3%, 37.1%, and 4.6%, respectively (
24). The frequencies of these aforementioned genotypes observed in this study were similar to those in the mentioned report (
24).
It is of particular importance to note that the thalassemic patients in this study showed a higher frequency of the rs12979860 T allele, which has been demonstrated to be associated with a poor treatment response (
2,
19). In addition, the lower HCV RNA level found in thalassemic patients might be interpreted in the context of the higher prevalence of cirrhosis and also the probable role of hepatic iron deposition in hepatocyte metabolism, as well as an additional possible role for the core protein in respect to enhancing hepatocyte iron levels (
58). Further prospective studies on separate naïve and non-naïve subgroups of HCV-infected thalassemic patients will be required to clarify the impacts of rs12979860 polymorphism on the course of CHC in thalassemic patients.
A significant link between core aa 70 substitution and the rs8099917 genotypes has been shown in previous studies, and treatment-resistant substitutions of core aa 70 (glutamine/histidine) were significantly more frequent in patients with treatment-resistant rs8099917 non-TT genotypes than the TT genotype (
45,
59,
60). It has also been demonstrated that the cumulative change rates from arginine to glutamine/histidine in core aa 70 were significantly higher in rs8099917 non-TT genotypes than in TT genotypes (
22). The present study agrees with these findings, as the rs8099917 non-TT genotypes had a higher frequency of core aa 70 substitutions; only a mild correlation was found. Interestingly, another correlation was also identified between the rs12979860 TT genotype and core aa 70 substitution. According to previous studies, the rs12979860 C allele is associated with a higher treatment response (
54,
61-
63). The association of these unfavorable factors may suggest an interaction between core aa 70 substitution and polymorphisms near the
IFNL3 gene as well as a possible role in determining the disease outcome. Future studies are required to determine the existence and clarification of this interaction.
The previous study of Hajia et al. (
11) on Iranian HCV-positive patients showed that HCV subtype 3a had a 30.76% frequency of high HCV RNA levels (≥ 800,000 IU/mL), while 69.22% of patients infected with HCV genotype 1 had elevated HCV RNA levels. Additionally, Chakravarti et al. (
64) reported that HCV genotype 1 had the highest HCV RNA level in comparison to HCV genotypes 2 and 3 in another Asian population in India. These findings agree with the present study in that HCV subtype 3a showed the lowest HCV RNA levels.
Based on recent meta-analyses, rs8099917 TT genotype and T allele contribute to higher rates of SVR (
65) and spontaneous HCV clearance in Caucasians mainly in individuals infected with HCV genotype 1 (
66). This study demonstrated an association between rs8099917 GG genotype and HCV subtype 1b, which is in line with previously mentioned studies and could be due to the aforementioned contributory role of the rs8099917 TT genotype and T allele. The finding of a higher frequency of the rs12979860 CC genotype in HCV subtype 3a-infected patients of this study was also in agreement with the results of a recent meta-analysis study, which reported an association of the rs12979860 CC genotype and C allele with HCV subtype 3a (
67), and also with similar findings of other studies on Iranian populations (
24,
26).
We are witnessing a rapid evolution in the treatment of CHC (
68). Emerging direct-acting antivirals (DAAs) have been taking promising steps towards improving the efficacy of therapeutic regimens, augmenting SVR rates to ≥ 90% (
69), and countervailing the low success rates and severe adverse reactions of previous treatments (
70). Selecting an appropriate therapeutic regimen based on the viral and host characteristics is emerging as the best practice routine for the management of HCV infection, but further studies and randomized clinical trials are still required in Iranian patients.
In this study, we faced some limitations; it was not possible to design a prospective study to assess the links between the molecular characteristics of HCV infection and the disease prognosis. Another limitation was the inevitable impact of possible time-related changes in HCV molecular characteristics, such as the substitution of core aa 70 after current assessment in this study, as the study population included both naïve and treatment-experienced patients and confronted the authors with a significant difference between these two groups regarding core aa 70 substitution and the rs12979860 genotypes, especially in hemophilic patients.
In conclusion, several associations between polymorphisms near the IFNL3 gene, HCV core aa 70 substitution and the HCV genotype were demonstrated. Moreover, it was found that HCV subtype 1a was the most prevalent subtype in Iranian patients with CHC. Arg70Gln was revealed to be the most common substitution and was significantly associated with an increased risk of cirrhosis. Additionally, the host characteristics (rs12979860 TT and rs8099917 GG genotypes) were linked with lower HCV RNA levels. Future studies seem reasonable, important, and essential to clarify any unrecognized and possibly important interactions between polymorphisms near the IFNL3 gene, the HCV genotype and core aa 70 substitution to reveal the future route of HCV infection therapy.