Open Access

Nucleotide identity and variability among different Pakistani hepatitis C virus isolates

  • Muhammad Idrees1Email author,
  • Sadia Butt1,
  • Zunaira Awan1,
  • Mahwish Aftab1,
  • Bushra Khubaib1,
  • Irshad-ur Rehman1,
  • Madiha Akram1,
  • Sobia Manzoor1,
  • Haji Akbar1,
  • Shazia Rafiqe1 and
  • Sheikh Riazuddin1
Virology Journal20096:130

https://doi.org/10.1186/1743-422X-6-130

Received: 10 July 2009

Accepted: 24 August 2009

Published: 24 August 2009

Abstract

Background

The variability within the hepatitis C virus (HCV) genome has formed the basis for several genotyping methods and used widely for HCV genotyping worldwide.

Aim

The aim of the present study was to determine percent nucleotide identity and variability in HCV isolates prevalent in different geographical regions of Pakistan.

Methods

Sequencing analysis of the 5'noncoding region (5'-NCR) of 100 HCV RNA-positive patients representing all the four provinces of Pakistan were carried out using ABI PRISM 3100 Genetic Analyzer.

Results

The results showed that type 3 is the predominant genotypes circulating in Pakistan, with an overall prevalence of 50%. Types 1 and 4 viruses were 9% and 6% respectively. The overall nucleotide similarity among different Pakistani isolates was 92.50% ± 0.50%. Pakistani isolates from different areas showed 7.5% ± 0.50% nucleotide variability in 5'NCR region. The percent nucleotide identity (PNI) was 98.11% ± 0.50% within Pakistani type 1 sequences, 98.10% ± 0.60% for type 3 sequences, and 99.80% ± 0.20% for type 4 sequences. The PNI between different genotypes was 93.90% ± 0.20% for type 1 and type 3, 94.80% ± 0.12% for type 1 and type 4, and 94.40% ± 0.22% for type 3 and type 4.

Conclusion

Genotype 3 is the most prevalent HCV genotype in Pakistan. Minimum and maximum percent nucleotide divergences were noted between genotype 1 and 4 and 1 and 3 respectively.

Background

Hepatitis C virus (HCV) belongs to the family Flaviviridae, genus Hepacivirus and is responsible for the second most common cause of viral hepatitis [1]. Presently, nearly 8-10% of Pakistani population [2], 2% of the USA population and 3% people worldwide are HCV carriers [3]. HCV has a positive-sense genome of approximately 9.6 kb and is subject to high rates of mutational changes [4]. Genetic heterogeneity of HCV isolated from different geographical regions was documented and at least six major genotypes with a series of subtypes of HCV have been identified so far [5]. The relative prevalence of these genotypes varies among different geographic regions such as subtypes 1a, 1b, 2a, 2c and 3a account for more than 90% of the HCV infections in North and South America, Europe, Russia, China, Japan, Australia, New Zealand and India [6, 7]. Type 4 is prevalent in Egypt, North Africa, Central Africa, and the Middle East; type 5 has been described in South Africa and type 6 is primarily found in Southeast Asia [8].

HCV variants studies have been made in the neighboring countries of Pakistan including India, Thailand, Vietnam, Indonesia and Burma and it is clear from all theses studies that type 1, type 2, type 3, and type 6 variants are prevalent in these areas [911]. From Pakistan few studies are available on the distribution of various hepatitis C virus genotypes [12, 13] however; none contained information on percent nucleotide identity among different isolates and geographic variation in the prevalence of various HCV genotypes. Therefore; 5'NCR sequence analysis followed by phylogenetic analysis was used for identifying different HCV variants, subtypes and genotypes in chronic HCV patients belonging to different geographical regions of Pakistan.

Methods

Patients and samples

One Hundred serum samples from chronic HCV carriers showing HCV RNA positivity and representing the four different areas of Pakistan such as Punjab (East), North West Frontier Province (NWFP) (North-west), Sindh (South-east) and Balochistan (South-west) were included in the study. The isolates from Punjab (number of isolates [n] = 25); NWFP (n = 25); Sindh (n = 25); or Balochistan (n = 25); are designated as P, N, S, or B, respectively, to identify the origin of the samples. A printed questionnaire was completed by each participant before the blood sample was collected after written informed consent. The study protocol was approved by the Institutional Ethical Committee. The demographic characteristics of the sequenced patients are shown in Table 1.
Table 1

Demographic characteristics of patients (N = 100).

Characteristics

Punjab

NWFP*

Sindh

Balochistan

Total (N = 100)

Sex-No. (%)

     

Male

13 (52)

15 (60)

11 (44)

18 (72)

57

Female

12 (48)

10 (40)

14 (56)

7 (28)

43

Age range-years

25-61

21-65

18-55

20-57

21-65

Mean age (Y)

± SD

40 ± 5.0

35 ± 7.0

47 ± 8.0

38 ± 9.8

43 ± 10.4

Socio-economic Status

     

No. (%)

     

Lower class

17 (68)

15 (60)

19 (76)

21 (84)

72

Middle class

08 (32)

10 (40)

06 (24)

04 (16)

28

Educational level No. (%)

     

Middle/above school

21 (84)

18 (72)

14 (56)

06 (24)

59

No/Primary school

04 (16)

07 (28)

11 (44)

19 (76)

41

Mode of contamination

     

No. (%)

     

Known

20 (80)

21 (84)

18 (72)

20 (80)

79

Unknown

05 (20)

04 (16)

07 (28)

05 (20)

21

History of previous Surgeries/dental procedure No. (%)

     

Yes

07(28)

03 (12)

06 (24)

03 (12)

19

No

18 (72)

22 (88)

19 (76)

22 (88)

81

Injected antibiotics/vitamins with used needle No. (%)

     

Yes

12 (48)

06 (24)

04 (16)

04 (16)

53

No

13 (52)

19 (76)

21 (84)

21 (84)

47

Blood transfusion/blood products No. (%)

     

Yes

02 (8)

01 (4)

00 (0)

01 (4)

06

No

23 (92)

24 (96)

25 (100)

24 (96)

94

HCV RNA level

     

< 400,000 IU/mL$

16 (64)

11 (44)

13 (52)

09 (36)

49

>400,000 IU/mL

09 (36)

14 (56)

12 (48)

16 (64)

51

Cirrhosis-No (%)

     

Present

03 (12)

02 (8)

05 (20)

2 (8)

12

Absent

22 (88)

23 (92)

20 (80)

23 (92)

88

Standard deviation

*NWFP, North West Frontier Province

$IU/mL, international units per milliliter

HCV RNA extraction and RT-PCR

HCV RNA was extracted from 100 μl serum sample using Gentra (Puregene, Minneapolis, MN 55441 USA) RNA isolation Kit according to the procedure given in the kit protocol. cDNA was synthesized at 37°C for 50 minutes using 1 μM of outer anti-sense primer and single tube nested PCR was done for 285-bp 5'NCR gene as described previously (Idrees et al. 2008). The PCR products were analyzed on 2% agarose gel.

Sequencing PCR of 5'UTR region

The purified DNA was used as templates for sequencing PCR in the Big-Dye Terminator cycle sequencing ready reaction kit (Applied Biosystems). Samples were analyzed on an automated sequencer (ABI PRISM 3100 genetic analyzer; Applied Biosystems). Products were sequenced from both strands to get consensus sequences. Placed the reaction tubes in thermal cycler (PE 2700, ABI) and set the volume to 20 μl. The samples were preheated at 96°C for one minute and then run 35 cycles with the following parameters: at 96°C for 10 seconds, 50°C for 5 seconds and 60°C for 4 minutes.

Purifying extension sample electrophoresis

The extension products were purified using ethanol precipitation method as described in the manual. Re-hydrated the pellet in 15 μl formamide and mixed well by up/down pipetting. Kept at room temperature for 15 minutes in dark. Heat denatured at 95°C for 5 minutes in thermal cycler and immediately put on ice for 5 minutes. The sequenced samples with BigDye terminators were electrophoresed on ABI PRISM 3100 instrument that is equipped with required modules and dye set/primer files.

Phylogenetic analysis

Pakistani isolates sequenced in the present study were aligned with the representative number of sequences for each major genotype and subtype selected from the GenBank database with the help of the Multalign program. Pairwise comparisons for percent nucleotide homology and evolutionary distance were made. The accession numbers of the prototype genotype sequences used to compare the 5' NC sequences were as follows: 1a, M62321; 1b, D90208; 2a, D00944; 2b, D01221; 2c, D10075; 3a, D14307; 3b, D11443; 3c, D16612; 4a, M84848; 4b, M84845; 4c, M84862; 4d, M84832; 4e, M84828; 4f, M84829; 5a, M84860; and 6a, M84827. The phylogenetic analysis of HCV isolates was performed with MEGA 3.0 software [14]. Jukes-Cantor algorithms were utilized, and phylogenetic trees were constructed by the neighbor-joining method. The reliability of different phylogenetic groupings was evaluated by using the bootstrap-resampling test from the MEGA program (1,000 bootstrap replications).

Results

On the basis of phylogenetic analysis, the 100 Pakistani isolates were classified as follows: 50% type 3, 9% type 1 and 6% type 4. Thirty five isolates still remained untypable (Fig 1). It was not possible to differentiate between type 1b and 1c isolates further into different subtypes as both types clustered together. In the case of the type 3 isolates, there was a clear clustering of isolates into subtypes 3a and 3b but still there were isolates that were not clustering to any of the subtypes and these may be new subtypes. Frequency distributions of HCV genotypes were not similar in all the four regions of the country as can be seen in table 2. In the North-west region 60% of isolates were not typed (Table 2).
Figure 1

Phylogenetic tree of HCV 5'UTR (nt 35 to 319) sequences of 100 HCV isolates. To identify the origins of the samples, the isolates of HCV patients belonged to areas of Punjab, N.W.F.P., Sindh or Balochistan are designated as PP, PN, PS or PB respectively. Sequences for each major subtype were selected from GenBank database for analysis. The accession numbers of the reference sequences are as follows: M67463 (1a), D90208 (1b), AY051292 (1c), AF238485 (2a), D82034 (2b), D10075 54 (2c), AF046866 (3a), D11443 (3b), D16612 (3c), D16620 (3d), D16618 (3e), D16614 (3f), X91421 (3g), Y11604 (4a), M84845 (4b), M84862 (4c), M84832 (4d), M84828 (4e), 84829 (4f), M8486 (5a), and Y12083 (6a).

Table 2

HCV genotypes prevailing in Pakistan based on 5' NCR$ sequence analysis (N = 100).

HCV Type

NWFP*

Punjab

Sindh

Balochistan

Total

1

2 (8%)

01 (4%)

03 (12%)

03 (12%)

9 (9%)

3

7 (28%)

17 (68%)

12 (48%)

14 (56%)

50 (50%)

4

01 (4%)

1 (4%)

3 (12%)

1 (4%)

6 (6%)

Not typed

15 (60%)

6 (24%)

7 (28%)

7 (28%)

35 (35%)

Total

25

25

25

25

100

HCV, hepatitis C virus

$NCR, Noncoding region

*NWFP, North West Frontier Province

The overall nucleotide similarity among these different Pakistani HCV sequenced isolates was 92.50% ± 0.50%. The percent nucleotide identity (PNI) was 98.11% ± 0.50% within Pakistani type 1 sequences, 98.10% ± 0.60% for type 3 sequences, and 99.80% ± 0.20% for type 4 sequences. The PNI between different genotypes was 93.90% ± 0.20% for type 1 and type 3, 94.80% ± 0.12% for type 1 and type 4, and 94.40% ± 0.22% for type 3 and type 4. There was a stretch of hypervariable region from nt: 83 to 171 in the 5'NCR of different HCV isolates. Pakistani isolates from different areas showed 7.5% ± 0.50% nucleotide variability in the sequenced 5'NCR region. The comparatively conserved stretch from nt 172 to 285 showed only 3.30% ± 1.06% variation. Minimum and maximum percent nucleotide divergences were noted between genotype 1 and 4 and 1 and 3. The sequence data of all the 100 sequences were submitted to GeneBank. The Accession Numbers provided for our nucleotide sequences by the GeneBank are from EF173931 to EF174030.

Discussion

HCV is an RNA virus is with a high rate of genetic mutation and extensive genetic heterogeneity of HCV exists in infected individuals as a result HCV isolates are found as either a group of isolates with very closely related genomes quasispecies, or distinct groups genetically called genotypes. It is believed that the different HCV variants are relevant to epidemiological questions, vaccine development, clinical management, therapeutic decisions and strategies. Due to this vital importance of HCV variants, the present study was carried out to identifying different HCV genotypes from Pakistan in particular to find out variability among HCV isolates of the same and different genotypes. In the present study we were able to successfully sequence and classify an excellent percent of specimens. Several findings emerged from this study. The first finding is the observation that the direct sequencing of amplification products provides more detailed sequence information and could be useful in the detection of new viral types and subtypes. Further, it is clear from the results of the present study that direct sequencing of the 5'UTR fragment allows good discrimination among the HCV major types. Due to the high degree of conservation found within 5' NCR this approach is not able to completely differentiate between all subtypes.

It is further clear from the findings of the present study that in Pakistan, HCV genotypes show differing distributions in different geographic regions. HCV genotypes 1, 3 and 4 have been detected with genotype 3 being most frequently detected. Although genotype 4 is found almost exclusively in Middle East and western countries [15] this genotype is uncommon in our country. Unexpectedly genotype 4 was seen very rare in Balochistan that is attached to Iran in the South-west where genotype 4 is the second major type existing in that area [16]. Another important finding is the observation of the absence of genotype 2 in all the four different regions of the country though not surprising as from neighbor countries like India and Iran genotype 2 is reported very rare [7, 16].

Next important finding of the present study is the isolation of many type 3 variants from Pakistan. The occurrence of many variants is not surprising because such type of variants have also been reported from neighboring countries particularly from India. The possibility of identifying more and more variants cannot be ruled out in the present situation of high prevalence of hepatitis C in this country. For this purpose, a study representing larger numbers of isolates from all provinces and community is required to generate countrywide data on HCV genotyping and variants.

Conclusion

We conclude that (i) multiple HCV genotypes are prevalent in Pakistan with genotype 3a as the predominant HCV genotype circulating in Pakistan, (ii) 5'NCR sequence analysis is sufficient for the routine genotyping of isolates in clinical settings; however, sequencing is very expensive and needs special laboratory settings, expertise and this method is unable to detect more than one genotype if present in the patient, (iii) Minimum and maximum percent nucleotide divergences were noted between genotype 1 and 4 and 1 and 3 respectively.

Abbreviations

HCV: 

hepatitis C virus

NCR: 

noncoding region

PNI: 

percent nucleotide identity

NWFP: 

North West frontier province

ABI: 

Applied Biosystem Inc.

RT-PCR: 

reverse transcriptase polymerase chain reaction

cDNA: 

complimentary DNA.

Declarations

Acknowledgements

This study was partially supported by Ministry of Science & Technology, Government of Pakistan. We thank all the subjects for their cooperation in the study.

Authors’ Affiliations

(1)
National Centre of Excellence in Molecular Biology, University of the Punjab Lahore

References

  1. Leiveven J: Pegasys/RBV Improves Fibrosis in Responders, relapsers & Nonresponders with Advanced Fibrosis. 55th Annual Meeting of the American Association for the Study of Liver Disease: 2004, October 29 - November 2: Boston, MA, USAGoogle Scholar
  2. Idrees M, Lal A, Naseem M, Khalid M: High prevalence of hepatitis C virus infection in the largest province of Pakistan. J Dig Dis 2008, 9: 96-104. 10.1111/j.1751-2980.2008.00329.xView ArticleGoogle Scholar
  3. Artini M, Natoli C, Tinari N, Costanzo A, Marinelli R, Balsano C, Porcari P, Angelucci D, D'Egidio M, Levrero M, Iacobelli S: Elevated serum levels of 90K/MAC-2 BP predict unresponsiveness to alpha-interferon therapy in chronic HCV hepatitis patients. J Hepatol 1996, 25: 212-217. 10.1016/S0168-8278(96)80076-6View ArticlePubMedGoogle Scholar
  4. Liew M, Erali M, Page S, Hillyard D, Wittwer C: Hepatitis C Genotyping by Denaturing High-Performance Liquid Chromatography. J Clinical Microbiol 2004,42(1):158-163. 10.1128/JCM.42.1.158-163.2004View ArticleGoogle Scholar
  5. Bukh J, Miller RH, Purcell RH: Genetic heterogeneity of hepatitis C virus: quasispecies and genotypes. Semin Liver Dis 1995, 15: 41-63. 10.1055/s-2007-1007262View ArticlePubMedGoogle Scholar
  6. Maertens G, Stuyver L: Genotypes and genetic variation hepatitis. In The molecular medicine of viral hepatitis. Edited by: Harrison TJ, Zuckerman A. John Wiley and Sons, Chichester, England; 1997:183-233.Google Scholar
  7. Chowdhury A, Santra A, Chaudhuri S, Dhali GK, Chaudhuri S, Maity SG, Naik TN, Bhattacharya SK, Mazumder DN: Hepatitis c virus infection in the general population: a community-based study in west bengal, india. Hepatology 2003,37(4):802-9. 10.1053/jhep.2003.50157View ArticlePubMedGoogle Scholar
  8. Chamberlain RW, Adams N, Saeed AA, Simmonds P, Elliot RM: Complete nucleotide sequence of a type 4 hepatitis C virus variant, the predominantgenotype in the Middle East. J Gen Virol 1997, 78: 1341-1347.View ArticlePubMedGoogle Scholar
  9. Tokita H, Okamoto H, Tsuda F, Song P, Nakata S, Chosa T, Iizuka H, Mishiro S, Miyakawa Y, Mayumi M: Hepatitis C virus variants from Vietnam are classifiable into the seventh, eighth, and ninth major genetic groups. Proc Natl Acad Sci USA 1994, 91: 11022-11026. 10.1073/pnas.91.23.11022PubMed CentralView ArticlePubMedGoogle Scholar
  10. Okamoto H, Tokita H, Sakamoto M, Horikita M, Kojima H, Mishiro S: Characterization of the genomic sequence of (or 3a) hepatitis C virus isolates and PCR primers for specific detection. J Gen Virol 1993, 74: 2385-2390. 10.1099/0022-1317-74-11-2385View ArticlePubMedGoogle Scholar
  11. Hotta H, Handajani R, Lusida MI, Soemarto W, Doi H, Miyajima H, Homma M: Subtype analysis of hepatitis C virus in Indonesia the basis of NS5b region sequences. J Clin Microbiol 1994, 32: 3049-3051.PubMed CentralPubMedGoogle Scholar
  12. Idrees M, Riazuddin S: Frequency distribution of hepatitis C virus genotypes in different geographical regions of Pakistan and their possible routes of transmission. BMC Infect Dis 2008, 8: 69. 10.1186/1471-2334-8-69PubMed CentralView ArticlePubMedGoogle Scholar
  13. Shah HA, Jafri W, Malik I, Prescott L, Simmonds P: Hepatitis C virus (HCV) genotypes and chronic liver disease in Pakistan. Gastroenterol. Hepatology 1997, 12: 758-761.Google Scholar
  14. Kumar S, Tamura K, Jakobsen IB, Nei M: MEGA2: molecular evolutionary genetics analysis software. Bioinformatics 2001, 17: 1244-1245. 10.1093/bioinformatics/17.12.1244View ArticlePubMedGoogle Scholar
  15. Mellor J, Holmes EC, Jarvis LM, Yap PL, Simmonds P: Investigation of the pattern of hepatitis C virus sequence diversity in different geographical regions: implications for virus classification. J Gen Virol 1995, 76: 2493-2507. 10.1099/0022-1317-76-10-2493View ArticlePubMedGoogle Scholar
  16. Samimi-Rad K, Nategh R, Malekzadeh R, Norder H, Magnius L: Molecular epidemiology of hepatitis C virus in Iran as reflected by phylogenetic analysis of the NS5B region. J Med Virol 2004, 74: 246-252. 10.1002/jmv.20170View ArticlePubMedGoogle Scholar

Copyright

© Idrees et al; licensee BioMed Central Ltd. 2009

This article is published under license to 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.

Advertisement