hepatitis c virus ( hcv ) infection is still one of the major causes of mortality and morbidity worldwide ( 1 , 2 ) . the world health organization ( who ) has estimated 3% prevalence of hcv equating approximately to 180 million individuals globally ( 3 - 5 ) . the prevalence of hcv infection in iran is estimated to be less than 0.5% ( 6 , 7 ) . it is well established that hepatitis c contributes to the increasing risk of fatal - related diseases including cirrhosis and hepatocellular carcinoma ( hcc ) ( 8 , 9 ) . in spite of numerous progressions in hepatitis c treatment , the high prevalence of hepatitis c in developing countries is still a major concern ( 10 ) . the only one open reading frame ( orf ) of hcv comprises about 9024 base pairs and encodes a polyprotein of about 3000 amino acids . up to now , hcv is classified into seven different genotypes and more than 67 subtypes according to the genetic variability and viral sequences ( 11 ) . the linkage of geographic distribution of hcv genotypes in different populations and specific risk groups to genetic diversity of hcv is obvious ( 12 ) . genotype of hcv is a considerable factor with clinical and epidemiological importance because it determines the rate of response to the hcv therapy , and may help trace the source of infection and clarify the possible modes of transmission ( 13 ) . in terms of response to the hcv treatment protocol with pegylated - interferon and ribavirin ( pegifn / rbv ) , the hcv genotype 1 ( hcv-1 ) and hcv-4 infections are more difficult - to - treat than hcv-2 and hcv-3 infections ( 10 , 14 , 15 ) . however , the development of new pangenotypic hcv treatments in recent years has contributed to the hcv elimination in the world ( 13 ) . the eastern mediterranean regional office ( emro ) is one of the six regional divisions of the who around the world , which serves 22 countries and territories in the middle east , the north africa , the horn of africa , and central asia with a total population of 605 million people . the frequency of hepatitis c infection , as estimated by the who , has revealed that at least 23 million individuals in the emro countries are infected with hcv ( 16 ) . the distribution of hcv genotypes in emro countries has two main patterns : arab countries ( except jordan ) with hcv-4 as the predominant isolate and non - arab countries with predominance of hcv genotypes other than hcv-4 such as iran by dominance of hcv-1 and afghanistan and pakistan by dominance of hcv-3 ( 17 , 18 ) . distribution of hcv genotypes in iran is different from most of the other neighboring and middle eastern countries , as the most frequent hcv genotype in iran is hcv-1a , followed by hcv-3a and -1b , respectively ( 19 ) . interestingly , in turkey , azerbaijan and russia , as neighboring countries in north of iran , the most prevalent genotype is hcv-1b ( 20 ) , in afghanistan and pakistan in east of iran , the most prevalent genotype is hcv-3 ( 21 ) , and in most of the neighboring arab countries in west and south of iran including iraq and saudi arabia , the most prevalent genotype is hcv-4 ( 18 , 22 , 23 ) . it is thought that the genetic diversity patterns of hcv in iran are similar to the pattern observed in north america and somewhat western europe ( 24 ) . the controversial issue regarding the source of hcv in iran in comparison with surrounding countries may be elucidated by phylogenetic analysis . there is heterogeneity in the regions sequenced along the hcv genome , such as core , ns5b , hvr-1 , e2 , and a segment of the ns5a gene associated with interferon sensitivity . guidelines propose to use either the full genome , core / e1 or ns5b sequences of hcv for classification of genotypes / subtypes ( 25 ) . furthermore , genotyping of hcv by nucleotide sequence analysis of ns5b is an effective procedure that allows discrimination of hcv subtypes properly . moreover , ns5b is an appropriate gene region to study the molecular epidemiology of hcv ( 25 ) . phylogenetic is a study for inferring or estimating evolutionary relationships among individuals or groups of organisms e.g. species or populations ( 26 ) . in the area of virology , phylogenetic trees contain a lot of information about the inferred evolutionary relationships between a set of viruses and it can be a potential widely used molecular tool to study rapidly - evolving rna viruses such as hcv . the main objective of the current study was to investigate the genetic relationship among all hcv-1a and -1b sequences derived from iran , emro , middle eastern , and some european and north american countries by applying phylogenetic analysis , and understanding of the source of spread of hcv-1a and -1b in iran . an electronic systematic search of available systematic reviews was conducted on all literature to find relevant studies reporting molecular prevalence and evaluation of hcv-1a and -1b in different iranian patient groups , and also all studies about hcv-1 molecular epidemiology in emro and middle eastern countries . to compare the obtained findings the search was performed on all peer - reviewed journals indexed in pubmed , scopus , and web of science databases . the literature review was carried out using the following key terms :  hepatitis c virus  ,  hcv  ,  genotype  ,  genotype 1a and 1b  ,  molecular sequence data  ,  sequence analysis  ,  phylogeography  ,  phylogenetic analysis  , and in addition to aforementioned search terms , the names of twenty - two emro countries were added to our search as follows : afghanistan , bahrain , djibouti , egypt , iran , iraq , jordan , kuwait , lebanon , libya , morocco , oman , pakistan , palestine , qatar , saudi arabia , somalia , sudan , syria , tunisia , united arab emirates , and yemen . moreover , middle eastern countries that were not in emro including cyprus and turkey were added in the search . furthermore , azerbaijan was added in the search , as a neighboring country of iran with wide commuting between the countries . also , some european and north american countries including france , germany , italy , netherlands , spain , uk , usa and canada , were added to the search strategy . in addition , to find appropriate sequences from unpublished studies which were registered in the genbank database , we searched into genbank using the names of aforementioned emro and middle eastern countries in addition to  hepatitis c virus  and  ns5b  . after all searches were completed , the ns5b sequences of all hcv-1a and -1b were selected and fasta format of the sequences were extracted from the genbank database . all published and unpublished studies with proper data were surveyed according to the following criteria : 1 ) all molecular studies in english among different patient groups with hcv-1a and -1b enrolled from iran , emro , middle eastern , and some aforementioned european and north american countries , 2 ) molecular studies that reported nucleotide sequence accession numbers based on ns5b gene sequences which were registered in the genbank database , 3 ) gene sequences with 243 base pair ( bp ) coverage which were between nucleotides 8319 - 8561 for hcv-1a ( the coverage obtained after full alignment along with hcv-1a reference with accession number m62321 ) , and 232 bp coverage between nucleotides 8315 - 8546 for hcv-1b ( the coverage obtained after full alignment along with hcv-1b reference with accession number u84014 ) . the mentioned frames were obtained after alignment and trimming of the included nucleotide sequences and exclusion of the sequences with less than 200 bp lengths . generally , the short nucleotide sequences or the nucleotide sequences which were not in the proper coverage were removed from the mega file . the exclusion criteria were as follows : 1 ) studies with possible errors and confusing data , 2 ) studies that used hcv genomic regions other than ns5b , 3 ) the hcv-2 to -6 isolates and also , hcv-1 other than hcv-1a and -1b isolates . to confirm the genotype / subtype of the selected isolates , we used ncbi viral genotyping tool ( http://www.ncbi.nlm.nih.gov/projects/genotyping/ formpage.cgi ) and hcv geno2pheno ( http://www.geno2pheno.org ) . for the eligibility criteria , all favorable articles obtained through the search strategy were independently reviewed by three authors ( kh - h , h - sh and a - njtsh ) . if there was any discrepancy between authors , it was resolved by consulting the supervisor of the study ( sma ) . the reviewing and screening processes in this study were based on the prisma guidelines for reporting systematic reviews ( 27 ) . we independently screened the title , abstract , and full - text of papers identified through the database searches . after full - text screening , the following data were extracted from each study : first author s name , publication year , country of origin , date of study , type of patient groups , and genbank accession numbers . the quality of the included studies was assessed using a modified strobe checklist ( 28 ) . in this study , the evolutionary relationships of isolates were inferred using the neighbor - joining method and kimura 2-parameter model . the percentage of replicate trees which were associated with clusters was estimated using bootstrap test ( 500 replicates ) . the results in this study were based on the clusters identified using the neighbor - joining phylogeny according to bootstrap test with a cut - off value 50% for defining the clusters ( values > 50% have been shown next to the branches ) . the trees were drawn to scale , with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree ( 29 ) . pair - wise and multiple alignments of the hcv ns5b sequences were performed using multiple sequence alignment - muscle ( multiple sequence comparison by log - expectation ) through molecular evolutionary genetics analysis software version 7.0 ( mega 7.0 ) ( 30 - 32 ) . to find the coverage and to edit the nucleotide sequences , the downloaded sequences were transferred to clc software ( clc main workbench 5 ) . following the full alignment of the sequences and manually trimming , the phylogenetic trees of suitable sequences were constructed by mega 7.0 software . initially , the trees were drawn with traditional - rectangular branch style . due to large number of data an electronic systematic search of available systematic reviews was conducted on all literature to find relevant studies reporting molecular prevalence and evaluation of hcv-1a and -1b in different iranian patient groups , and also all studies about hcv-1 molecular epidemiology in emro and middle eastern countries . to compare the obtained findings the search was performed on all peer - reviewed journals indexed in pubmed , scopus , and web of science databases . the literature review was carried out using the following key terms :  hepatitis c virus  ,  hcv  ,  genotype  ,  genotype 1a and 1b  ,  molecular sequence data  ,  sequence analysis  ,  phylogeography  ,  phylogenetic analysis  , and in addition to aforementioned search terms , the names of twenty - two emro countries were added to our search as follows : afghanistan , bahrain , djibouti , egypt , iran , iraq , jordan , kuwait , lebanon , libya , morocco , oman , pakistan , palestine , qatar , saudi arabia , somalia , sudan , syria , tunisia , united arab emirates , and yemen . moreover , middle eastern countries that were not in emro including cyprus and turkey were added in the search . furthermore , azerbaijan was added in the search , as a neighboring country of iran with wide commuting between the countries . also , some european and north american countries including france , germany , italy , netherlands , spain , uk , usa and canada , were added to the search strategy . in addition , to find appropriate sequences from unpublished studies which were registered in the genbank database , we searched into genbank using the names of aforementioned emro and middle eastern countries in addition to  hepatitis c virus  and  ns5b . after all searches were completed , the ns5b sequences of all hcv-1a and -1b were selected and fasta format of the sequences were extracted from the genbank database . all published and unpublished studies with proper data were surveyed according to the following criteria : 1 ) all molecular studies in english among different patient groups with hcv-1a and -1b enrolled from iran , emro , middle eastern , and some aforementioned european and north american countries , 2 ) molecular studies that reported nucleotide sequence accession numbers based on ns5b gene sequences which were registered in the genbank database , 3 ) gene sequences with 243 base pair ( bp ) coverage which were between nucleotides 8319 - 8561 for hcv-1a ( the coverage obtained after full alignment along with hcv-1a reference with accession number m62321 ) , and 232 bp coverage between nucleotides 8315 - 8546 for hcv-1b ( the coverage obtained after full alignment along with hcv-1b reference with accession number u84014 ) . the mentioned frames were obtained after alignment and trimming of the included nucleotide sequences and exclusion of the sequences with less than 200 bp lengths . generally , the short nucleotide sequences or the nucleotide sequences which were not in the proper coverage were removed from the mega file . the exclusion criteria were as follows : 1 ) studies with possible errors and confusing data , 2 ) studies that used hcv genomic regions other than ns5b , 3 ) the hcv-2 to -6 isolates and also , hcv-1 other than hcv-1a and -1b isolates . to confirm the genotype / subtype of the selected isolates , we used ncbi viral genotyping tool ( http://www.ncbi.nlm.nih.gov/projects/genotyping/ formpage.cgi ) and hcv geno2pheno ( http://www.geno2pheno.org ) . for the eligibility criteria , all favorable articles obtained through the search strategy were independently reviewed by three authors ( kh - h , h - sh and a - njtsh ) . if there was any discrepancy between authors , it was resolved by consulting the supervisor of the study ( sma ) . the reviewing and screening processes in this study were based on the prisma guidelines for reporting systematic reviews ( 27 ) . we independently screened the title , abstract , and full - text of papers identified through the database searches . after full - text screening , the following data were extracted from each study : first author s name , publication year , country of origin , date of study , type of patient groups , and genbank accession numbers . the quality of the included studies was assessed using a modified strobe checklist ( 28 ) . in this study , the evolutionary relationships of isolates were inferred using the neighbor - joining method and kimura 2-parameter model . the percentage of replicate trees which were associated with clusters was estimated using bootstrap test ( 500 replicates ) . the results in this study were based on the clusters identified using the neighbor - joining phylogeny according to bootstrap test with a cut - off value 50% for defining the clusters ( values > 50% have been shown next to the branches ) . the trees were drawn to scale , with branch lengths in the same units as those of the evolutionary distances used to infer the phylogenetic tree ( 29 ) . pair - wise and multiple alignments of the hcv ns5b sequences were performed using multiple sequence alignment - muscle ( multiple sequence comparison by log - expectation ) through molecular evolutionary genetics analysis software version 7.0 ( mega 7.0 ) ( 30 - 32 ) . to find the coverage and to edit the nucleotide sequences , the downloaded sequences were transferred to clc software ( clc main workbench 5 ) . following the full alignment of the sequences and manually trimming , the phylogenetic trees of suitable sequences were constructed by mega 7.0 software . initially , the trees were drawn with traditional - rectangular branch style . due to large number of data in the present study , hcv ns5b nucleotide sequences from 30 various studies that had been totally conducted on more than 3000 subjects were collected . we studied all available nucleotide sequences reported from iran and other emro and middle eastern countries . out of seven studies from iran , four studies were published in 2004 - 2014 ( 33 - 36 ) and three were unpublished . furthermore , 15 published and unpublished studies from emro and middle eastern countries including afghanistan ( 37 ) , azerbaijan , cyprus ( 38 , 39 ) , egypt , morocco ( 40 ) , pakistan ( 41 ) , and tunisia ( 42 , 43 ) were collected . furthermore , other studies from the us and european countries were investigated , randomly ( 44 - 49 ) . based on figure 1 , after removal of 187 duplicates , 382 irrelevant titles , and 81 papers with irrelevant abstracts , finally , 14 studies were eligible to be assessed in phylogenetic analyses ( 33 - 43 , 50 - 52 ) . in total , 505 sequences of hcv-1a and -1b were obtained from these 14 studies . moreover , after searching the genbank database to find unpublished studies with registered sequences , 1831 sequences were obtained . after removal of 1283 sequences ( including 47 sequences because of non - ns5b gene sequences and 1236 sequences because of hcv genotypes other than hcv-1a and -1b ) , a total of 548 sequences were obtained . therefore , a total of 1053 sequences from both published and unpublished studies were collected . furthermore , 195 different sequences from italy , france , netherlands , spain , the united kingdom , and the united states were randomly included ( 44 - 49 ) . after appropriate sequence collection , all of the 835 sequences were transferred into mega software for full alignment . finally , 161 sequences were removed because they were not in the coverage setting ; thus , a total of 674 sequences were used in the phylogenetic analyses including 442 sequences for hcv-1a and 232 sequences for hcv-1b . although 8 included studies were unpublished ( from iran , azerbaijan , egypt , pakistan , and tunisia ) , the nucleotide sequences of these studies had been registered in the genbank database . the characteristics of unpublished studies with registered sequences in the genbank database are shown in table 2 . abbreviations : idu , intravenous drug user ; na , not available . out of 442 extracted sequences for hcv-1a , 325 ( 73.5% ) sequences were obtained from iranian studies . these sequences were derived from different groups including blood donors , inherited bleeding disorders , thalassemia , intravenous drug users , patients on hemodialysis , and patients without known risk factors . there were sequences from the uk intravenous drug users clustered with some french blood donors and some sequences from iranian patients in the phylogenetic tree . also , the specific sequences from iranian patients clustering with the uk and cypriots intravenous drug users were observed ( figure 2 ) . it can be concluded from the tree that a proportion of the iranian isolates clustered along with each other . a comparison of sequences of hcv-1b is shown in the phylogenetic tree in figure 3 . seventy one ( 30.6% ) sequences of total 232 sequences were isolates from iranian patients . these sequences were isolated from blood donors , inherited bleeding disorders , thalassemia , and patients without known risk factors . based on the phylogenetic analysis , the iranian sequences of hcv-1b had heterogeneous dispersion ( figure 3 ) . some of the iranian sequences clustered with each other and some clustered with european sequences particularly sequences from france , spain and italy . most likely , the hcv-1b isolates from iranian patients may have similarities with the european ones . also , in the phylogenetic tree there are isolates from different geographical regions which clustered together . in the present study , hcv ns5b nucleotide sequences from 30 various studies that had been totally conducted on more than 3000 subjects were collected . we studied all available nucleotide sequences reported from iran and other emro and middle eastern countries . out of seven studies from iran , four studies were published in 2004 - 2014 ( 33 - 36 ) and three were unpublished . furthermore , 15 published and unpublished studies from emro and middle eastern countries including afghanistan ( 37 ) , azerbaijan , cyprus ( 38 , 39 ) , egypt , morocco ( 40 ) , pakistan ( 41 ) , and tunisia ( 42 , 43 ) were collected . furthermore , other studies from the us and european countries were investigated , randomly ( 44 - 49 ) . based on figure 1 , after removal of 187 duplicates , 382 irrelevant titles , and 81 papers with irrelevant abstracts , finally , 14 studies were eligible to be assessed in phylogenetic analyses ( 33 - 43 , 50 - 52 ) . in total , 505 sequences of hcv-1a and -1b were obtained from these 14 studies . moreover , after searching the genbank database to find unpublished studies with registered sequences , 1831 sequences were obtained . after removal of 1283 sequences ( including 47 sequences because of non - ns5b gene sequences and 1236 sequences because of hcv genotypes other than hcv-1a and -1b ) , a total of 548 sequences were obtained . therefore , a total of 1053 sequences from both published and unpublished studies were collected . furthermore , 195 different sequences from italy , france , netherlands , spain , the united kingdom , and the united states were randomly included ( 44 - 49 ) . after appropriate sequence collection , all of the 835 sequences were transferred into mega software for full alignment . finally , 161 sequences were removed because they were not in the coverage setting ; thus , a total of 674 sequences were used in the phylogenetic analyses including 442 sequences for hcv-1a and 232 sequences for hcv-1b . although 8 included studies were unpublished ( from iran , azerbaijan , egypt , pakistan , and tunisia ) , the nucleotide sequences of these studies had been registered in the genbank database . the characteristics of unpublished studies with registered sequences in the genbank database are shown in table 2 . out of 442 extracted sequences for hcv-1a , 325 ( 73.5% ) sequences were obtained from iranian studies . these sequences were derived from different groups including blood donors , inherited bleeding disorders , thalassemia , intravenous drug users , patients on hemodialysis , and patients without known risk factors . there were sequences from the uk intravenous drug users clustered with some french blood donors and some sequences from iranian patients in the phylogenetic tree . also , the specific sequences from iranian patients clustering with the uk and cypriots intravenous drug users were observed ( figure 2 ) . it can be concluded from the tree that a proportion of the iranian isolates clustered along with each other . a comparison of sequences of hcv-1b is shown in the phylogenetic tree in figure 3 . seventy one ( 30.6% ) sequences of total 232 sequences were isolates from iranian patients . these sequences were isolated from blood donors , inherited bleeding disorders , thalassemia , and patients without known risk factors . based on the phylogenetic analysis , the iranian sequences of hcv-1b had heterogeneous dispersion ( figure 3 ) . some of the iranian sequences clustered with each other and some clustered with european sequences particularly sequences from france , spain and italy . most likely , the hcv-1b isolates from iranian patients may have similarities with the european ones . also , in the phylogenetic tree there are isolates from different geographical regions which clustered together . currently , hcv is considered as one of the most important viruses threatening human life . the distribution pattern of hcv genotypes is various among hcv - infected individuals which depends on different status of public health behavior and social risk factors . it is clear to all that the predominance of risk factors for hcv transmission has changed over time , from blood transfusion to intravenous drug use ( 24 ) . the distribution of hcv genotypes and subtypes in iran and other middle eastern and emro countries has a very diverse pattern . altogether , the hcv-1 ( -1a and -1b ) is the predominant genotype in iran , so that more than half of the hcv - infected patients in iran are infected with this genotype ( with a rate of 54% ) ( 12 , 18 ) . the genetic diversity of hcv is due to the unique characteristic of the rna molecule . the genetic variation stems from the error - prone ns5b polymerase . as a result , quasispecies  are produced , almost with a single mutation in each cycle of replication . the production of highly different viruses through the dynamic replication process of hcv will occur with a count of 10 trillion viruses per day ( 53 ) . the results of the detailed phylogenetic analysis by using ns5b sequences of hcv-1a indicated that a proportion of the iranian hcv-1a isolates was in common clades . a proportion of iranian hcv-1a isolates most probably has domestic origin  . in this study , these results showed that the hcv-1a sequences from afghan patients were likely similar to the isolates from iran . this indicates that the hcv-1a sequences from both countries are closely related to each other genetically . this may be because of the fact that the afghanistan land was a part of iranian territory years ago and now after occurring fled wars in afghanistan since 1978 , iran s border gates were opened on afghan refugees and it provided the conditions for the large numbers of afghans to immigrate to iran . perhaps this could be one of the main explanations for the hcv genotypes similarity in afghanistan and iran . in pakistan it is plausible that the most of pakistani patients who were infected with hcv-1a , acquired this infection due to unsafe medical practices during surgeries ( 54 ) . in a previous phylogenetic analysis of hcv-1a in pakistan , the virus was identified with polyphyletic origin , and the sequences were found to be closely related to european strains ( 55 ) . cyprus , morocco , tunisia , and turkey are the countries wherein hcv-1b is predominant . the emro and middle eastern countries take diverse patterns of hcv genotype distribution , dominantly hcv - g4 in arab countries to hcv - g3 and -g1 ( -1a and -1b ) in non - arab countries ( 18 ) . as mentioned earlier , some of the hcv-1b sequences of iranian isolates were similar to counterpart sequences from european isolates including those from france and spain . it seems that this similarity is more likely among iranian patients with inherited congenital bleeding disorder such as hemophilia . the prevalence of hcv among iranian hemophilia patients is high , as it is said that the overall prevalence in these iranian patients group is 40.8% , with a range from 13.3% to 80.5% ( 8 , 56 ) . it is noteworthy that in the 1980s , with the arrival of blood and blood products from france to iran , a large number of patients , particularly hemophiliac patients , were infected with hiv and hcv . previously , it was shown that hcv-1b is more frequently observed in iranian hemophiliac patients than other iranian hcv infected groups ( 57 ) . this may indicate the possible similarity of hcv-1b of iranian and european isolates and suggests an infection through blood products such as clotting factors imported to iran . phylogenetic analysis of hcv genotypes and subtypes is a useful molecular method which helps scientists in every geographic region provide a substantial contribution to monitoring the virus for any purpose including hcv molecular tracing , ancestral studies , performing different genetic assessments , and guidance for any treatment decision . there are some limitations in this study : 1 ) the short ( < 300 bp ) nucleotide sequences which were used in the phylogenetic analyses , 2 ) hcv ns5b sequences were available from limited number of emro and middle eastern countries , and 3 ) we could not establish a proper analysis for assessment of the relationship between hcv risk factors and the phylogeny of hcv ns5b sequences . however , more studies should be conducted in next future to find more genetic relationships between these different sequences from different regions and patient groups including genetic distances for measuring genetic divergence , phylodynamic inference , and evolutionary methods to define circulating strains and molecular clock analysis for understanding the ancestral relationships . in conclusion , the ns5b sequences of different infected - patients were phylogenetically- evolutionarily analyzed for molecular tracing of hcv-1 in iran . the phylogenetic trees of hcv-1a and -1b according to 500 pseudo - replicates indicated many clades and codon positions with ancestral relationships of all data . phylogenetic reconstruction of all sequences of hcv-1 pinpoints phylogenetic dispersion of most of hcv-1b of iranian isolates among other european ones with a considerable diversity ; whereas , most of iranian hcv-1a isolates are genetically defined probably with domestic origin .
contexthepatitis c virus ( hcv ) is characterized by a high degree of genetic heterogeneity and classified into 7 genotypes and different subtypes . it heterogeneously distributed through various risk groups and geographical regions . 
 a well - established phylogenetic relationship can simplify the tracing of hcv hierarchical strata into geographical regions . 
 the current study aimed to find genetic phylogeny of subtypes 1a and 1b of hcv isolates based on ns5b nucleotide sequences in iran and other members of eastern mediterranean regional office of world health organization , as well as other middle eastern countries , with a systematic review of available published and unpublished studies.evidence acquisitionthe phylogenetic analyses were performed based on the nucleotide sequences of ns5b gene of hcv genotype 1 ( hcv-1 ) , which were registered in the genbank database . 
 the literature review was performed in two steps : 1 ) searching studies evaluating the ns5b sequences of hcv-1 , on pubmed , scopus , and web of science , and 2 ) searching sequences of unpublished studies registered in the genbank database.resultsin this study , 442 sequences from hcv-1a and 232 from hcv-1b underwent phylogenetic analysis . 
 phylogenetic analysis of all sequences revealed different clusters in the phylogenetic trees . 
 the results showed that the proportion of hcv-1a and -1b isolates from iranian patients probably originated from domestic sources . 
 moreover , the hcv-1b isolates from iranian patients may have similarities with the european ones.conclusionsin this study , phylogenetic reconstruction of hcv-1 sequences clearly indicated for molecular tracing and ancestral relationships of the hcv genotypes in iran , and showed the likelihood of domestic origin for hcv-1a and various origin for hcv-1b .