- Research article
- Open Access
- Open Peer Review
Molecular characterization of Hepatitis C virus 3a in Peshawar
BMC Infectious Diseases volume 16, Article number: 163 (2016)
The purpose of this study was to explore molecular epidemiology of HCV genotype 3a in Peshawar based on sequencing and phylogenetic analysis of Core region of HCV genome.
Chronically infected Hepatitis C virus infected patients enrolled under the Prime Minister Hepatitis C control program at three Tertiary care units of Peshawar [Khyber Teaching Hospital Peshawar, Lady Reading Hospital Peshawar, Hayat Abad Medical Complex Peshawar] were included in this cross sectional observational study. Qualitative detection of HCV and HCV genotyping was carried out by a modified reverse transcription-polymerase chain reaction (RT-PCR) and type specific genotyping assay. The Core gene of HCV genotype 3a was amplified, cloned and sequenced. The sequences obtained were used for phylogenetic analysis using MEGA 6 software.
Among the 422 (82.75 %) PCR positive samples, 192 (45.5 %) were identified as having HCV genotype 3a infection. HCV Core gene sequencing was carried out randomly for the characterization of HCV 3a. Nucleotide sequence analysis of the obtained viral genomic sequences based on partial HCV 3a Core gene sequences with reference sequences from different countries showed that our sequences clustered with some local and regional sequences with high bootstrap values.
HCV 3a is highly prevalent in Peshawar, Pakistan and its phylogenetics based on Core gene sequences indicate the prevalence of different lineages of HCV 3a in Peshawar which may have consequences for disease management strategies causing more economic pressure on the impoverished population due to possible antiviral resistance.
Hepatitis C virus (HCV), first identified in 1989 is the primary cause of chronic liver diseases including cirrhosis (60–85 %) and hepatocellular carcinoma , with an annual mortality rate of 3.5-5.5 million due to complications of end stage liver diseases . An estimated 185 million people (3 % of the world’s population) are infected worldwide with relatively high prevalence rates in developing countries [3, 4]. In Pakistan, approximately 10–17 million people are chronically infected with HCV [5, 6] with an overall prevalence rate of 5 % . Frequency of HCV infection differs among the four provinces with an estimated 1.1 million people infected in KPK [8, 9].
HCV is an enveloped, positive sense ssRNA genome  in the family Flaviviridae consisting of an oral reading frame (ORF) approximately of 9.6 kb in length . HCV genome shows substantial nucleotide sequence variability in both the structural and nonstructural coding regions, with different isolates of HCV showing as much as 30 % nucleotide sequence divergence over the entire genome which is sufficient to alter the antigenic and biological characteristics of seven major genotypes of HCV . These genotypes have distinct geographical distributions. Although HCV genotypes 1, 2 and 3 appear to have a worldwide distribution, their relative prevalence varies from one geographic area to another . Identification of HCV genotype is extremely important as different genotypes are relevant to the epidemiology and clinical management of chronic HCV infection and is the strongest predictive parameter for sustained virological response .
The primary means to identify and classify new genotypes is either by phylogenetic analysis of sequences or by measures of pairwise sequence similarity . Although it is generally true that longer sequences are more informative for classification, it is usually possible to identify genotypes by sequence comparison of relatively short sub genomic regions of HCV . While region encoding the 5′ non-coding region (5′ NCR) is too conserved for HCV subtypes identification, sequencing and phylogenetic analysis of nucleotide sequences amplified in the region of the genome encoding the Core protein is frequently used for classification of HCV into genotypes and subtypes [17–19]. HCV Core gene has sufficient genetic diversity and can produce topologically identical trees to those obtained upon analysis of complete genome sequences .
Molecular epidemiological studies previously conducted in Peshawar, Pakistan have indicated that the most prevalent HCV subtype is 3a accounting for 70–90 % of HCV infections based solely on type specific PCR based genotyping methods developed 15 to 20 years ago [21, 22]. At present, there is little information about the genetic history and evolution of HCV genotype 3a in Peshawar which needs to be investigated for several reasons including increasing antiviral resistance in the case of HCV 3a in Peshawar, KP province and its evolutionary relationship with other regional or global isolates. As the prevalence of HCV 3a has been reported to be very high among the general population of Peshawar, therefore, we embarked on the current study to investigate it at molecular level by Type-specific assay, sequencing of the Core gene and phylogenetic analysis in order to have a clear epidemiological picture of the prevalent HCV 3a isolates.
Chronically infected Hepatitis C virus infected patients enrolled under the Prime minister Hepatitis C control program at three Tertiary care units [Khyber Teaching Hospital Peshawar, Lady Reading Hospital Peshawar, Hayat Abad Medical Complex Peshawar] willingly participated in this observational study. A non-probability convenient sampling technique was used to collect data. Written consent was obtained from all the eligible study participants for participation and publication of the data and the Institutional Ethics Committee (Khyber Medical University Ethics Board) approved the study. Blood samples were collected from the patients in sterile vacutainers and sera extraction was carried out at the Institute of Biotechnology and Genetic Engineering Peshawar.
Total RNA was extracted from each sample using FavorGen RNA isolation Kit (Favorgen Biotech corp, Taiwan, CAT No FAVNK 001) according to the manufacturer’s instructions. The isolated RNA was stored either stored at-80C or immediately used for RT PCR.
Reverse Transcription PCR followed by Type-specific PCR for identification of HCV subtypes was carried out essentially according to  with modifications in the Universal and Type-specific primers according to the latest HCV sequencing data (Table 1).
Amplification of HCV 3a Core gene
A representative number of samples (58/192), which turned out positive for HCV 3a by the Type-specific assay described earlier, were used for a second round of RNA extraction, RT PCR and subsequent regular nested PCR for the amplification of Core gene using gene-specific primer .
Cloning and sequencing of the Core gene
The Core gene products were eluted from agarose gel using Pure Link™ Quick Gel Extraction Kit (https://www.thermofisher.com/pk/en/home.html). The eluted products were cloned in pGEMT-Easy vector and transformed in DH5-alfa strain of E.coli (Promega). Sequencing PCR of the Core gene was carried out using vector-specific sense and antisense primers in a sequencing PCR (Big Dye Deoxy Terminator method). Bidirectional sequencing run was performed using ABI PRISM 310 DNA sequencer [Applied Biosystems].
In case of each isolate, three independent clones were sequenced. The respective sequences were aligned and consensus sequences were derived. The sequences were blasted using NCBI Blast tool and genotype authentication was further validated with the help of an online software HCV Genotyper (http://www.bioafrica.net/rega-genotype/html/subtypinghcv.html). The confirmed and curetted sequences representing the diversity of 58 sequenced isolates were submitted to GenBank
HCV Core gene sequences from various geographical regions of the world were retrieved from the online sequence repositories and Phylogenetic tree was constructed by using Mega6 .
SPSS version 20 was used for data analysis. The qualitative variables were described using percentage and the quantitative variables were described using the mean, median and standard deviation.
For identification of HCV genotype 3a infection among the actively infected HCV samples, RT PCR followed by Type specific nested PCR based genotyping assay was carried out. Genotyping analysis revealed that out of the 422 (82.75 %) samples processed. Out of 422 actively infected patients, 257 (61 %) were females and 165 (39 %) were male) HCV genotype 3a was identified as the most abundant HCV subtype 192 (45.50 %). Mean age of patients included in the study was 38 ± 11.8 (Mean ± SD. For authentication of HCV genotype 3a, Core gene amplification and sequencing was carried out by randomly selecting a representative number (58/192) of samples and it was established that the patients had HCV 3a infection. These results were also validated by an online sub-typing tool (http://www.bioafrica.net/rega-genotype/html/subtypinghcv.html). Detailed sequence and phylogenetic analysis of the obtained viral genomic sequences was carried out. Phylogenetic tree was developed with Maximum Likelihood algorithm (1000 bootstrap replicates). To analyze HCV genotype 3a of Peshawar origin, the available subtype 3a sequences of the Core gene were retrieved from NCBI and aligned with the respective subtype identified in this study. Pair wise nucleotide and deduced amino acid sequence comparison of our sequences with 36 reported sequences from different countries was performed which indicated that some of the HCV 3a isolates of Peshawar clustered together with each other and with isolates from Pakistan (KP869047, KP869048, KP869049, KP869055, KP869056, KP869057, KP869058), while others clustered closer to European, Indian, Iranian or Chinese HCV sequences with high bootstrap support (Fig. 1). Two of the representative isolates (KP869053 and KP869054) were clustered with the European isolates. Isolate KP869050 and KP869051 shared homology with both the Indian, Iranian and Chinese isolates. Moreover isolate KP869052 was found closely related to those from Iran and Pakistan.
Accumulation of nucleotide substitutions in the HCV genome results in diversification and evolution into seven major genotypes and a series of subtypes . The classification of HCV by viral genotype is not only important for appropriate treatment regimen and assessment of global viral evolution but their epidemiological data can reveal transmission activity and migration movement of infected individuals from the endemic area [27–29]. These viral types and subtypes show differing distribution in different geographic regions which have provided investigators with an epidemiologic marker that can be used to trace the source of HCV infection in a given population . Sequencing and phylogenetic analysis of HCV Core gene nucleotide sequences has earlier been used for identification and classification of HCV isolates into genotypes and subtypes . HCV genotype 3a is the most abundant in Pakistan and it has earlier been documented that HCV genotype 3a is relatively more responsive to therapy as compared to genotype 1 or 4 [28, 29, 31]. However, in Pakistan resistance is being observed in case of relatively responsive and easy to treat HCV genotype 3a .
The present study is the first report describing the genotypic and evolutionary analysis of HCV 3a isolated from chronic HCV infected patients of Peshawar, based on a modified genotyping assay and subsequent sequencing and phylogenetic analysis. Investigation of the different circulating genotypes and their evolution is not only crucial for epidemiological and clinical analysis but might be helpful for the improvement of diagnostic tests and treatment regimens . HCV infection is highly prevalent in Pakistan with an overall prevalence rate of 5 % . The estimated HCV prevalence in KP province is 1.1 % . In this study, out of total 510 seroactive patients investigated, 82.75 % were actively infected with HCV having high proportion of HCV genotype 3a (45.50 %) which is in agreement with most recent studies in Peshawar and Pakistan reporting HCV 3a to be the most prevalent genotype [35, 36], however the percentage prevalence is lower than the previous reports claiming much higher prevalence of HCV genotype 3a (60–74 %) [22, 37]. One possible explanation for this changing HCV genotype 3a landscape could be the change in epidemiological pattern over times as a result of people migration  or the inadequate sensitivity of old genotyping assays to correctly type the circulating viral strains due to substantial genetic heterogeneity inherent to RNA viruses . Moreover the emerging resistance to interferon therapy experienced in case of HCV genotype 3a  might possibly be due to the less prevalent genotypes other than HCV 3a or variants of HCV 3a evolved overtimes which can cause substantial economic and health burden over the infected population. Internationally accepted guidelines for the treatment of hepatitis C are rarely followed in KPK and people undergo successive therapies once the initial response of antiviral drugs is negative . Due to lack of awareness, the practice mentioned above is not only causing economic and health related losses but is also contributing towards the evolution of more resistant HCV 3a types.
Molecular evolutionary analysis of the obtained viral genome sequences revealed that some HCV 3a isolates of Peshawar clustered closer to local isolates (Fig. 1) indicating the previous existence of similar types in other parts of Pakistan, which may have spread to Peshawar, KPK province via various transmission routes. Some of the HCV 3a isolates grouped with European, south Asian, Iranian and Chinese HCV 3a isolates (Fig. 1). Peshawar city is located in the northwestern region of Pakistan. It shares international border with the Afghanistan which has been home to conflicts for the past 40 years resulting in migration of various ethnic groups into and out of Afghanistan. These migrations have changed epidemiological patterns of various pathogens overtimes including HCV. Phylogenetic analysis indicate that foreign presence in Afghanistan and migration of Afghanis to Peshawar might have contributed towards the spread of isolates which are genetically closer to European, Indian, Chinese and Iranian isolates as HCV genotype 3a is also highly prevalent in neighboring countries including China , India  and Iran .
This study has some limitations. Due to convenient sampling, selection bias might have occurred. Moreover we have characterized HCV 3a entirely on the bases of partial core gene sequences. To elucidate the epidemiology of emerging HCV 3a in Peshawar and to further improve the accuracy of diagnostic assays and treatment regimens, there is a need to analyze complete coding sequences of more diverse regions of HCV genome on a comparatively large scale.
This study concludes that HCV 3a is highly prevalent in Peshawar, Pakistan and its phylogenetics based on Core gene sequences indicate the prevalence of different lineages of HCV 3a in Peshawar which may have consequences for disease management strategies causing more economic pressure on the impoverished population due to possible antiviral resistance.
Availability of data and materials
The datasets supporting the conclusions of this article are available in the GenBank [National Center for Biotechnology Information] repository under the following accession numbers [KP869047, KP869048, KP869049, KP869050, KP869051, KP869052, KP869053, KP869054, KP869055, KP869056, KP869057, KP869058]
hepatitis C virus
Petta S, Craxi A. Current and future HCV therapy: do we still need other anti-HCV drugs? Liver Int. 2015;35(1):4–10.
Kamal-Yanni M. Hepatitis C, drug affordability. Lancet Glob Health. 2015;3(2):e73–4.
Gower E, Estes C, Blach S, Razavi-Shearer K, Razavi H. Global epidemiology and genotype distribution of the hepatitis C virus infection. J Hepatol. 2014;61(1 Suppl):S45–57.
Hill A, Khoo S, Fortunak J, Simmons B, Ford N. Minimum costs for producing hepatitis C direct-acting antivirals for use in large-scale treatment access programs in developing countries. Clin Infect Dis. 2014;58(7):928–36.
Ali A, Nisar M, Ahmad H, Saif N, Idrees M, Bajwa MA. Determination of HCV genotypes and viral loads in chronic HCV infected patients of Hazara Pakistan. Virol J. 2011;8:466.
Athar MA, Xu Y, Xie X, Xu Z, Ahmad V, Hayder Z, et al. Rapid detection of HCV genotyping 1a, 1b, 2a, 3a, 3b and 6a in a single reaction using two-melting temperature codes by a real-time PCR-based assay. J Virol Methods. 2015;222:85–90.
Nouroz FS, S Mujtaba G, Noreen S. An overview on hepatitis C virus genotypes, Egypt. J Med Hum Genet. 2015;16(4):291–8.
Umar M, Bilal M, Hepatitis C. A Mega Menace: A Pakistani Perspective. J Pioneer Med Sci. 2012;2(2):68–72.
Salam A, Khan B, Sarwar N, Masoom G. Seroprevalnce of HBsAg and anti HCV in general population of Pakistan during 2012 and 2013. Prof Med J. 2014;21(4):766.
Quer J, Gregori J, Rodriguez-Frias F, Buti M, Madejon A, Perez-del-Pulgar S, et al. High-resolution hepatitis C virus subtyping using NS5B deep sequencing and phylogeny, an alternative to current methods. J Clin Microbiol. 2015;53(1):219–26.
Tuplin A, Struthers M, Cook J, Bentley K, Evans DJ. Inhibition of HCV translation by disrupting the structure and interactions of the viral CRE and 3′ X-tail. Nucleic Acids Res. 2015;43(5):2914–26.
Echeverria N, Moratorio G, Cristina J, Moreno P. Hepatitis C virus genetic variability and evolution. World J Hepatol. 2015;7(6):831–45.
Rasheed A, Ullah S, Naeem S, Zubair M, Ahmad W, Hussain Z. Occurrence of HCV genotypes in different age groups of patients from Lahore, Pakistan. Adv Life Sci. 2014;1(2):89–95.
Messina JP, Humphreys I, Flaxman A, Brown A, Cooke GS, Pybus OG, et al. Global distribution and prevalence of hepatitis C virus genotypes. Hepatology. 2015;61(1):77–87.
Simmonds P, Mellor J, Sakuldamrongpanich T, Nuchaprayoon C, Tanprasert S, Holmes EC, et al. Evolutionary analysis of variants of hepatitis C virus found in South-East Asia: comparison with classifications based upon sequence similarity. J Gen Virol. 1996;77(12):3013–24.
Marascio N, Torti C, Liberto M, Foca A. Update on different aspects of HCV variability: focus on NS5B polymerase. BMC Infect Dis. 2014;14(5):S1.
Li HC, Lo SY. Hepatitis C virus: Virology, diagnosis and treatment. World J Hepatol. 2015;7(10):1377–89.
Saludes V, Gonzalez V, Planas R, Matas L, Ausina V, Martro E. Tools for the diagnosis of hepatitis C virus infection and hepatic fibrosis staging. World J Gastroenterol. 2014;20(13):3431–42.
Simmonds P, Bukh J, Combet C, Deleage G, Enomoto N, Feinstone S, et al. Consensus proposals for a unified system of nomenclature of hepatitis C virus genotypes. Hepatology. 2005;42(4):962–73.
Simmonds P, Holmes EC, Cha TA, Chan SW, McOmish F, Irvine B, et al. 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(11):2391–9.
Waqas M, Parveen F, Waqar M, Ali A, Khan MA, Wasim M, et al. Distribution Pattern of Various Genotypes of HCV Circulating in District Mardan, Khyber Pakhtunkhwa, Pakistan. Advan Biol Res. 2015;9(2):69–74.
Khan S, Attaullah S, Ayaz S, Niaz Khan S, Shams S, Ali I, et al. Molecular epidemiology of hcv among health care workers of khyber pakhtunkhwa. Virol J. 2011;8:105.
Ohno O, Mizokami M, Wu RR, Saleh MG, Ohba K, Orito E, et al. New hepatitis C virus (HCV) genotyping system that allows for identification of HCV genotypes 1a, 1b, 2a, 2b, 3a, 3b, 4, 5a, and 6a. J Clin Microbiol. 1997;35(1):201–7.
Kanwal S, Mahmood T. Hepatitis C viral heterogeneity based on core gene and an attempt to design small interfering RNA against strains resistant to interferon in rawalpindi, pakistan. Hepat Mon. 2012;12(6):398–407.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30(12):2725–9.
Smith DB, Bukh J, Kuiken C, Muerhoff AS, Rice CM, Stapleton JT, et al. Expanded classification of hepatitis C virus into 7 genotypes and 67 subtypes: updated criteria and genotype assignment web resource. Hepatology. 2014;59(1):318–27.
Sharafi H, Alavian S, Keshvari M. Efficacy of 24-Week Pegylated Interferon Alpha and Ribavirin Combination Therapy in Highly Selected Patients Infected With Hepatitis C Virus Genotype 1. Hepat Mon. 2015;15(1), e24955.
Humphreys IS, Von Delft A, Brown A, Hibbert L, Collier JD, Foster GR, et al. HCV genotype-3a T cell immunity: specificity, function and impact of therapy. Gut. 2012;61(11):1589–99.
Wasitthankasem R, Vongpunsawad S, Siripon N, Suya C, Chulothok P, Chaiear K, et al. Genotypic distribution of hepatitis C virus in Thailand and Southeast Asia. PLoS One. 2015;10(5), e0126764.
Zein NN. Clinical significance of hepatitis C virus genotypes. Clin Microbiol Rev. 2000;13(2):223–35.
Abro AH, Al-Dabal L, Younis NJ. Distribution of hepatitis C virus genotypes in Dubai, United Arab Emirates. J Pak Med Assoc. 2010;60(12):987–90.
Ali S, Ali I, Azam S, Ahmad B. Frequency distribution of HCV genotypes among chronic hepatitis C patients of Khyber Pakhtunkhwa. Virol J. 2011;8:193.
Khan A, Tanaka Y, Azam Z, Abbas Z, Kurbanov F, Saleem U, et al. Epidemic spread of hepatitis C virus genotype 3a and relation to high incidence of hepatocellular carcinoma in Pakistan. J Med Virol. 2009;81(7):1189–97.
Afridi SQ, Ali MM, Awan F, Zahid MN, Afridi IQ, Afridi SQ, et al. Molecular epidemiology and viral load of HCV in different regions of Punjab, Pakistan. Virol J. 2014;11:24.
Afridi SQ, Zahid MN, Shabbir MZ, Hussain Z, Mukhtar N, Tipu MY, et al. Prevalence of HCV genotypes in district Mardan. Virol J. 2013;10:90.
Ali A, Ahmed H, Idrees M. Molecular epidemiology of Hepatitis C virus genotypes in Khyber Pakhtoonkhaw of Pakistan. Virol J. 2010;7:203.
Khan TM, Mehr MT, Ullah H, Khan H, Iman NU. Frequency of hepatitis C genotypes in the North of Pakistan. Gomal J Med Sci. 2014;12(2):106–9.
Yan Z, Fan K, Wang Y, Fan Y, Tan Z, Deng G. Changing pattern of clinical epidemiology on hepatitis C virus infection in southwest china. Hepat Mon. 2012;12(3):196–204.
Chakravarti A, Ashraf A, Malik S. A study of changing trends of prevalence and genotypic distribution of hepatitis C virus among high risk groups in North India. Indian J Med Microbiol. 2013;31(4):354–9.
Jahanbakhsh Sefidi F, Keyvani H, Monavari SH, Alavian SM, Fakhim S, Bokharaei-Salim F. Distribution of hepatitis C virus genotypes in Iranian chronic infected patients. Hepat Mon. 2013;13(1), e7991.
Authors acknowledge support of the Directorate of Science and Technology KP, for the study [Project No. Dirtt: S&T/KP/PRS/FINANCIAL ASSISTANCE/2011-12].
None of the authors have any competing interests in the manuscript.
IA and JA designed the study. AG, IA, FZ and NZ carried out the experimental work. JA, IA and IAK helped with data analysis, writing the MS and reviewing it thoroughly. All authors read and approved the final MS.
IA is a Fulbrighter working at CIIT Islamabad as Associate Professor and Leader of the Virology group. JA is a professor and director of the Institute of Basic Medical Sciences, Khyber Medical University Peshawar. His major interest is microbiology. AG is a PhD scholar while FZ and NZ are graduate students who have worked in IA’s lab. IAK is a professor at the University of Agriculture Peshawar and his major interest is entomology and microbiology.