Article Text
Abstract
Background and aims A total of 105 patients were identified as accidentally infected with hepatitis C virus genotype 1b (HCV1b) through blood transfusion from a single blood donor. This group provides a unique patient population to study host factors involved in the spontaneous clearance of HCV and disease progression.
Methods Clinical markers, HCV RNA and eight single nucleotide polymorphisms (SNPs) of interleukin-28B (IL-28B) were detected. Exome capture and sequencing were analysed for association with HCV clearance.
Results Among the 85 patients with the positive HCV antibody, 27 cases (31.8%) were HCV RNA negative over a period of 9–12 years. Of the 58 patients with positive HCV RNA, 22.4% developed chronic hepatitis, and 5.2% developed cirrhosis. Age was found to be associated with HCV1b clearance. IL-28 rs10853728 CC showed the trend. By exon sequencing, 39 SNPs were found to be significantly different in spontaneous clearance patients (p<0.001). Two SNPs in the tenascin receptor (TNR), five in the transmembrane protease serine 11A (TMPRSS11A), and one in the serine peptidase inhibitor kunitz type 2 (SPINT2) showed the closest associations (p<10−5).
Conclusions Host genetic analyses on the unique, single source HCV1b-infected patient population has suggested that age and mutations in TNR, TMPRSS11A and SPINT2 genes may be factors associated with HCV clearance.
- HCV
- GENETICS
- CHRONIC HEPATITIS
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/
Statistics from Altmetric.com
Summary box
What is already known about this subject?
▸ Host interleukin-28B (IL-28B) polymorphisms were known to be associated with spontaneous hepatitis C virus (HCV) clearance and also response to treatment. HCV is the other factor contributing to clearance. When both host and viral factors are mixed involving in HCV clearance and disease progression, it is difficult to tell the important factors.
What are the new findings?
▸ This is a the study on a unique group of patients with HCV1b-infection (n=105) accidentally transmitted from a single blood donor infected with genotype 1b in Guizhou province, southwest China. With the sole resource of the virus, the clear-known infected time, the similar ethnicity and environments, it is better to understand the host factors for HCV spontaneous clearance and disease progress.
How might it impact on clinical practice in the foreseeable future?
▸ Add the knowledge of how the host factors may affect HCV clearance.
Introduction
Hepatitis C virus (HCV) infection affects hundreds of millions of people worldwide. It has been reported that about 20% of HCV-infected adults can spontaneously clear the virus, while 30% of patients with chronic infection progress to cirrhosis and hepatocellular carcinoma (HCC).1
Viral and host factors are involved in HCV spontaneous clearance and disease progression. Virus factors include HCV genotypes, quasispecies, viral load and co-infection. Host factors include gender, age at infection or the ageing process, race, alanine aminotransferase (ALT) elevation and genetic factors.2 Recently, interleukin-28B (IL-28B) polymorphisms have been reported to be associated with spontaneous HCV clearance and also response to treatment.3–7
The purpose of the current study was to analyse a group of patients infected with the same HCV genotype 1b (HCV1b) source in order to focus on host parameters that may be involved in resolution or persistence of HCV infection. These patients in the current study are unique for several reasons. First, the sole resource of the HCV1b virus excludes virus genotypic differences. HCV1b is a difficult-to-treat genotype with interferon-based therapy. Second, the known date of infection provides data on the natural history of HCV infection over the course of 9–12 years. Third, the common ethnicity and similar environments of the patients reduce some variables into the analysis of factors involved in HCV spontaneous clearance and disease progression. Lastly, the broad age range of patients is helpful to study the importance of host age.
Materials and methods
Study subjects
All patients had received blood transfusions, from 1998 to 2002, from a single blood donor who was subsequently found to have had HCV1b. All recipients were Chinese from Pingtang, Guizhou province, southwest China. Inclusion criteria were transfusion of blood or blood-products from the identified contaminated batches of the same donor. Patients who died from causes other than HCV-related liver disease, and patients we were unable to contact, were excluded. All patients with positive HCV RNA were tested and found to be genotype 1b.
Patients were identified and blood samples were collected from 2010 to 2011, 9–12 years postinfection (median 10 years).
The study was approved by the ethics committee from Guiyang Medical College and conformed to the ethical guidelines of the Declaration of Helsinki; informed consent had been obtained from each individual included in the study.
Serum HCV antibody and RNA assays
Serum biochemical parameters including ALT levels were measured by routine automated methods according to the manufacturer's instructions. Anti-HCV antibody levels were measured by electrochemiluminescence immunoassay (ECLIA) using Abbott Architect i2000 (ABBOTT, Wiesbaden, Germany) according to the manufacturer's instructions. HCV RNA was detected by the commercial quantitative reverse transcription PCR (RT-PCR; COBAS AMPLICOR, Roche Diagnostic Systems, Indianapolis, USA) according to the manufacturer's instructions. The lower limit of detection was 15 IU/mL. The HCV genotype was determined by Versant HCV Genotype 2.0 (LiPA; Siemens Healthcare Diagnostics, Tarrytown, New York, USA).
Fibroscan detection
Transient elastography was performed using FibroScan (Echosens, France). The examination was performed on the right lobe of the liver through the seventh or eighth intercostal space. The measurement depth was between 25 and 65 mm. As suggested by the manufacturer, only results obtained with 10 valid measurements, with a success-rate of at least 60% and with an IQR ≤30%, were considered reliable, as described previously.8–10
PCR amplification and sequencing of IL-28B polymorphisms
Genomic DNAs were isolated from 0.5 mL whole blood using the QIAamp DNA Mini Kit (Qiagen, Hilden, Germany). IL-28 rs12979860, rs8099917, rs10853728, rs12980275, rs4803219, rs4803223, rs8105790 and rs28416813 were amplified by PCR. The PCR protocol involved initial denaturation at 95°C for 10 min, 35 cycles of denaturation for 30 s at 95°C, annealing of primers for 30 s at 55°C and extension for 40 s at 72°C, followed by final extension at 72°C for 10 min. The amplified DNA fragments were separated on a 2% agarose gel, and purified with the QIAquick gel extraction kit (Qiagen, Hilden, Germany) according to the manufacturer's instructions Nucleotide sequences were determined by Sanger sequencing using the Applied Biosystems Automated 3730 DNA Analyzer.
Exome capture and sequencing assay
Exome capture
Purified genomic DNA samples were randomly fragmented by Covaris, with the size of library fragments mainly distributed between 150 and 200 bp.11–13 Adapters were ligated to both ends and purified by the Agencourt AMPure SPRI beads according to the manufacturer's instructions. Fragments with insert sizes of about 250 bp were excised. Extracted DNA was amplified by ligation-mediated PCR (LM-PCR), purified and hybridised to the SureSelect Biotinylated RNA Library (BAITS) for enrichment. Hybridised fragments were bound to the strepavidin beads, whereas non-hybridised fragments were washed out after 24 h. Captured LM-PCR products were analysed using an Agilent 2100 Bioanalyzer to estimate the magnitude of enrichment. Each captured library was then loaded on a Hiseq2000 platform, to ensure that each sample meets the desired average sequencing depth for high-throughput sequencing. Raw image files were processed by Illumina basecalling Software V.1.7 with default parameters and the sequences of each individual were generated as 90 bp pair-end reads.11
Read mapping
SOAPaligner (V.2.21) was used to align the sequencing reads to the NCBI human genome reference assembly (build 36.3) with a maximum of 3 mismatches and the parameters were set as -a -b -D -o -u -p -2 -m -x -s 40 -l 35 -v 3. To evaluate exome capture efficiency, the proportions of reads mapping to target regions and to their flanking regions (within 200 bp) were calculated for each individual. Reads that aligned to the designed target region (TR±200 bp) were collected for single nucleotide polymorphism (SNP) identification and subsequent analysis.11
Individual genotype calling
On the basis of SOAP alignment results, the software SOAPsnp was used to call genotypes. The following parameters were set: -r 0.0005 -e 0.001 -t -u -2 -i -d -o -M -L 90 -s -T (http://soap.genomics.org.cn/ for details). To obtain an accurate genotype set, the genotype was filtered on the basis of the following criteria: the SNP should be observed in at least one individual with a quality score >20, and in a way that the number of reads containing mutant alleles was larger than the reads containing reference alleles. At the same time, at least 90% of all individuals got a quality score >20 and coverage >4.12 ,13 Our next steps are all based on this genotype set.
Principal component analysis
To check the population stratification, we did a principal component analysis with EIGENSTRAT software. We chose the SNP from 1000 genomes (http://www.1000genomes.org/) to check our samples.
Determination of associations
We checked the genome-wide association results using PLINK software (the http://pngu.mgh.harvard.edu/~purcell/plink/) to perform a standard case/control association analysis using: plink --file mydata --assoc. All the top sites SNPs with statistical significance were also checked for Hardy-Weinberg equilibrium, and sites which had p<0.001 were filtered.12 ,13
Statistical analysis
Clinical statistical analyses were performed using the Statistical Program for Social Sciences (V.11.5; SPSS). Continuous variables such as ALT, HCV RNA and time were presented as medians (range), and categorical variables as frequencies. Categorical variables were tested by χ2 test or Fisher's exact test. p Values <0.05 were considered to be statistically significant.
Results
Patient demographic data
The demographic features of the patients are shown in figure 1 and table 1. A total of 105 receipts were found with the single same donor from medical records of the blood transfusions from 1998 to 2002. Twenty of them with a negative HCV antibody could not be contacted or died from causes other than liver disease. Among the 85 patients with a positive HCV antibody, the male-to-female ratio was 45:40 with a median age of 32 years (9–71 years). Twenty-seven cases (31.8%) were HCV RNA negative, and were considered to have achieved spontaneous clearance of HCV. The other 58 cases (68.2%) were HCV RNA positive with levels ranging from 3.0 log to 7.1 log IU/mL, median 4.8 log IU/mL.
Characteristics of HCV spontaneous clearance
Among the 85 patients with a positive HCV antibody, host factors including the gender, age and IL-28 allele were analysed for an association with spontaneous clearance. It was shown that an age less than 20 years had the most significant association with spontaneous clearance (OR=2.04, 95% CI (1.13 to 3.69), p=0.028), while other ages did not have any significant relation to virus clearance (table 1). Gender had no relation with viral clearance (p=0.462; figure 2).
Eight IL-28B polymorphisms, alleles rs12979860, rs8099917, rs10853728, rs12980275, rs4803219, rs4803223, rs8105790 and rs28416813 were PCR amplified and sequenced. As shown in table 2, IL-28B rs10853728 showed a strong trend, which fell just short of statistical significance (p=0.058). For the most well-known SNP, IL-28B rs12979860, 22 of the 24 (91.7%) patients who had spontaneous clearance had the CC allele, while 31 of the 40 (77.5%) in the non-clearance group had the CC allele (p=0.132). The TT allele was not found in any of the patients (figure 3).
Risk factors for HCV disease progression
Nine to 12 years after transfusion, none of the 27 patients with spontaneous clearance had disease progression by liver function tests, abdominal ultrasonography and fibroscan tests. Among the 58 patients with positive HCV RNA, 9 cases (15.5%) had elevated serum ALT ranging from 41 to 192 IU/mL. A total of 13 cases (22.4%) developed chronic hepatitis with mild to moderate fibrosis as determined by clinical manifestations, fibroscan values higher than 7.1 kPa, and enhanced and coarse echogenicity of the liver by ultrasonography. A total of three patients (5.2%) developed cirrhosis with decreased albumin, fibroscan values higher than 9.5 kPa, and splenomegaly as determined by ultrasound.8–10 Neither decompensated cirrhosis nor HCC was found. According to the infected time and duration, it was estimated that the rate of HCV progression to mild or moderate fibrosis was 2.2% per year, and to cirrhosis was 0.5% per year.
Multivariate regression analysis showed that gender (p=0.393), HCV RNA level (p=0.262) and IL-28B allele frequencies (p=0.565) were not statistically associated with disease progression. An age less than 40 years (OR=0.13, 95% CI (1.13 to 3.69), p=0.020) had a negative association with disease progression.
Whole-exome capture and sequencing
Using exome capture and sequencing, 64 449 SNPs were identified in the sample population. A total of 17 081 coding genes were sequenced with coverage of each individual exome at an average depth of 33.9-fold. On average, about 95% of the target regions were covered by at least one read. More than 86% of the target regions were covered by more than 4 reads (figure 4). The population analysis showed that samples used for the association analysis had no significant population stratification.
The top 20 SNPs which had the closest association are listed in table 3. There were SNPs from 11 exons from IL-28B rs2239818 and rs34842046, two from the tenascin receptor (TNR), one (rs3745948) from the serine peptidase inhibitor kunitz type 2 (SPINT2), one (rs7627615) from the 5-hydroxytryptamine receptor 3 family member E (HTR3E), four (rs1370840, rs11930532, rs28437478 and rs6552134) from the transmembrane protease serine 11A (TMPRSS11A), one (rs1263810) from the sal-like protein 2 (SALL2), three (rs9901726, rs2291604 and rs9900543) from the spermatogenesis associated 22 (SPATA 22), one (rs607332) from the nicotinamide mononucleotide adenyltransferase 2 (NMNAT2), two (novel) from the NCK-associated protein 1 (NCKAP1), one (rs2303225) from the MARVEL domain containing 3 (MARVELD3), three (novel) from the zinc finger protein 491, 440 and 439 (ZNF491, ZNF 440 and ZNF 439), and one (rs2307075) from the carbonate dehydratase II (CA2; p<10−4).
From the function analysis, there were five SNPs of receptors, TNR, HTR3E and MARVELD3, and four SNPs from the transmembrane protease TMPRSS11A, which may affect HCV binding to the receptors and entry into hepatocytes (p<10−4).
Discussion
The differences between HCV genotypes coupled with the large diversity of host factors make an analysis of favourable and morbid outcomes in populations difficult to determine and analyse. Our study focused on a specific group of patients with HCV1b-infection due to a unique sole source reservoir of virus, the clear known time of infection, the similar ethnic and environmental background, and the broad scale of 85 patients including both sexes and almost all ages. With the same HCV source, this information may help us to better study the host factors involved in virus clearance and/or disease progress, in exclusion of HCV viral differences.
There are only two studies comparable to the current study, one from Ireland and one from Germany. Both were reports of single-source infections with HCV1b involving homogeneous women of childbearing age who had received HCV-contaminated anti-D immunoglobulin to prevent Rh isoimmunisation from 1977 to 1979. The Irish group had 376 women with a mean age of 28±6 years, who had been infected for about 17 years.14 The German cohort had 1018 women with a median age of 24 years, who had been followed for 20 years.15 Both studies provided very important information on HCV spontaneous clearance and disease progression in large groups of patients with known durations of infection.16–22 However, owing to the relatively uniform age and gender, these two studies provide less data on how age and gender affect HCV spontaneous clearance and disease progression.
IL-28B has recently been found to be a promising gene marker associated with treatment response and spontaneous clearance.3–7 However, data on the association between IL-28B genetic variants and spontaneous viral clearance in a Chinese study involving 376 HCV-infected paid plasma donors did not show an association with rs12979860. The other four SNPs, rs8099917, rs8105790, rs12980275 and rs10853728, were significantly associated with spontaneous HCV clearance.23 In our study, IL-28 rs10853728 showed a stronger association than the other seven SNPs, but the association fell short of statistical significance (p=0.058). This may be because the number of patients in this cohort was not large enough to provide statistical power to the differences. Also, in the highly prevalent IL-28B, favourable genotype area and, in particular, the Chinese population (>80%) may be different from the Caucasian population (40–50%).3 ,4
TNR rs2239818 and rs34842046 are involved with receptor binding, extracellular matrix organisation and negative regulation of cell adhesion in the cell surface or extracellular region.24 ,25 Transmembrane protease TMPRSS11A rs1370840, rs11930532, rs28437478, rs6552134 and rs977728 are expressed in the normal liver, oesophagus, colon and lung, but downregulated in tumours.26 ,27 Neither TNR nor TMPRSS11A had been reported to be associated with spontaneous clearance of HCV. SPINT2 rs3745948 functions as an endopeptidase inhibitor within the plasma membrane, cytoplasm or extracellular region. Methylated SPINT2 and SRD5A2, combined with AFP and PIVKA-II, have been reported to be the most satisfactory panel to detect HCC in patients with chronically infected HCV.28 ,29 The mechanisms by which these three categories of SNPs result in the associations are not known. However, data from functional analyses suggest that they may affect HCV binding to the receptors and entering hepatocytes.
There were two main limitations to the study. The first was the small number of patients. Though we cannot increase the sample size of this unique group since HCV were accidentally infected, the limited statistical power made us interpret the data more conservatively. Future studies, based on another larger scale cohort of patients with HCV and controls, are ongoing to verify those observations. Since this was a retrospective study, the other limitation was that we could not get the clinical data from the acute HCV phase about 10 years ago. The development of jaundice, other symptoms and laboratory findings during the acute phase could not be evaluated for spontaneous HCV clearance.
In summary, this unique single-source HCV1b-infected patient population allows analysis of HCV1b with spontaneous clearance from a new perspective. Host gene SNPs within Tenascin-R, TMPRSS11A, and SPINT2 and IL-28 most likely play roles in the HCV spontaneous clearance and disease progression.
References
Footnotes
HY, SL and YX are the first coauthors.
The abstract of this paper has been presented in the 22nd Congress of the Asian Pacific Association for the Study of the Liver (APASL 2012, PS04-03), and published in the abstract book.
Contributors HY contributed to patient management and analysis of the risk factors. SL was involved in local patient management. YX conducted local patients’ follow-up and maintained records. RC conducted the IL-28 SNPs analysis. YaS contributed to IL-28 detection. JR conducted the Exome capture and sequencing experiments and analysis. KW conducted SNPs data analysis. XJ contributed to SNPs comparison and analysis. YuS conducted bioinformatics analysis. HaZ was involved in HCV RNA detection. JL conducted patient follow-up and is a consultant. LW was involved in HCV RNA analysis and disease stage confirmation. HuZ was involved in study design and is a consultant. MC contributed to sample management and patient management. JJ contributed to the study design and final approval of the paper.
Funding This work was supported by the Program for National Science and Technology Major Project (2013ZX10002004, 2012ZX10002003), Key Project from the Education Bureau of Beijing (KZ201210025024).
Competing interests None.
Patient consent Obtained.
Ethics approval The study was approved by the Ethics Committee from Guiyang Medical College, Guizhou, China.
Provenance and peer review Not commissioned; externally peer reviewed.
Data sharing statement No additional data are available.