J Arthropod-Borne Dis, March 2016, 10(1): 12–23 S Chinikar et al.: Assessment of Recombination … 12 Original Article Assessment of Recombination in the S-segment Genome of Crimean- Congo Hemorrhagic Fever Virus in Iran *Sadegh Chinikar 1, Nariman Shah-Hosseini 1, Saeid Bouzari 2, Mohammad Ali Shokrgozar 3, Ehsan Mostafavi 4, Tahmineh Jalali 1, Sahar Khakifirouz 1, Martin H Groschup 5, Matthias Niedrig 6 1Arboviruses and Viral Hemorrhagic Fevers Laboratory (National Reference Lab), Pasteur Institute of Iran, Tehran, Iran 2Department of Molecular Biology, Pasteur Institute of Iran, Tehran, Iran 3National Cell Bank of Iran, Pasteur Institute of Iran, Tehran, Iran 4Department of Epidemiology, Pasteur Institute of Iran, Tehran, Iran 5Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, Greifswald-Insel Riems, Germany 6Robert Koch Institute, Berlin, Germany (Received 21 May 2014; accepted 29 June 2014) Abstract Background: Crimean-Congo Hemorrhagic Fever Virus (CCHFV) belongs to genus Nairovirus and family Bunyaviridae. The main aim of this study was to investigate the extent of recombination in S-segment genome of CCHFV in Iran. Methods: Samples were isolated from Iranian patients and those available in GenBank, and analyzed by phyloge- netic and bootscan methods. Results: Through comparison of the phylogenetic trees based on full length sequences and partial fragments in the S- segment genome of CCHFV, genetic switch was evident, due to recombination event. Moreover, evidence of multi- ple recombination events was detected in query isolates when bootscan analysis was used by SimPlot software. Conclusion: Switch of different genomic regions between different strains by recombination could contribute to CCHFV diversification and evolution. The occurrence of recombination in CCHFV has a critical impact on epidemi- ological investigations and vaccine design. Keywords: Recombination, Phylogenetic, Diversity, Crimean Congo Hemorrhagic Fever, Iran Introduction Crimean Congo Hemorrhagic Fever (CCHF) is caused CCHF virus (CCHFV), which is a real threat for human life. CCHFV is a tick-borne virus in the family Bunyaviridae, genus Nairovirus (Kinsella et al. 2004). Endemic regions of CCHF have been reported in Africa, the Middle East, Eastern Europe and Western Asia where the vector and/ or reservoir ticks of Hyalomma spp. are distributed (Elevli et al. 2010, Mehr- avaran et al. 2013, Champour et al. 2014a). In the last 10 years, Turkey, Bosnia, and Iran have been reported the most frequent out- breaks of CCHF worldwide (Chinikar et al. 2012a). Clinical symptoms of CCHF include headache, high fever, back pain, joint pain, stomach pain and vomiting (Swanepoel et al. 1989). An infected tick remains infected through- out its life and transmits the infection to large vertebrates. Livestock play a significant role in virus amplification because the ani- mals become viremic for seven days (Cham- pour et al. 2014b). *Corresponding author: Dr Sadegh Chinikar, E-mail: sadeghchinikar@yahoo.com http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 12–23 S Chinikar et al.: Assessment of Recombination … 13 The viral genome is segmented, compris- ing three circular single strands of negative sense RNA (Flick et al. 2003). The small (S) segment (approximately 1.7 kb) has a single open reading frame encoding the nucleocapsid (N) protein. In contrast, the medium (M) seg- ment (approximately 5 kb) encodes a large polyprotein, which is processed into the two surface glycoproteins (Gn and Gc) and sev- eral non-structural proteins. The large (L) segment encodes a single L protein of ap- proximately 460 kDa, which is probably the viral polymerase (Nuttall et al. 2001). To date, different phylogenetic analysis studies on CCHF virus S segment have re- vealed that strains from different regions of the world cluster into several distinct phylo- genetic groups. The genetic diversity of CCHFV can be explained by mutation, reassortment and recombination (Chinikar et al. 2004, Deyde et al. 2006). Although mutations are among the most important reasons for increasing RNA virus genome diversity, RNA reassortment, mix- ing of the genetic material of a species into new combinations, of segmented viruses is also a critical mechanism for rapid novel vi- rus creation. In this regard, one recent paper reported convincing evidence for genome reassortment in CCHFV (Hewson et al. 2004). Concerning recombination in CCHFV, a preliminary study reported the possibility of recombination in CCHFV, without provid- ing clear evidence (Chare et al. 2003). The results of other studies depicted evidence of potential recombination among S segment in CCHFV (Deyde et al. 2006). The only solid evidence of recombination in S-segment of CCHFV was eventually demonstrated by Lukashev (2005). There have been several articles on the molecular epidemiology of CCHFV employ- ing a phylogenetic approach, but only few studies discussed the reasons behind the observed genetic diversity (Chinikar et al. 2010, Chinikar et al. 2004). Given the poten- tial importance of recombination, it is im- portant to determine its frequency for CCHFV. According to the results of Lukashev study, we decided to focus our recombination anal- ysis only on the S-segment by analyzing the extent of recombination events in sequences from samples which were isolated from Ira- nian patients. Materials and Methods Viral extraction, RT-PCR and sequencing Initially, viral RNA was extracted from 140 μl of serum using a QIAamp RNA Mini kit, according to of the manufacturer's in- structions (QI Agen GmbH, Hilden, Ger- many) (Yashina et al. 2003). Specific primers for amplification of full-length CCHFV ge- nome (S-segment) were designed by CLC main workbench software version 5.0. PCR products were amplified using one-step RT- PCR, according to Rodriguez et al. (Rodri- guez et al. 1997). The amplified products were visualized by ethidium bromide agarose gel staining (Yadav et al. 2012). The PCR products were then sequenced using Big Dye Terminator V3.1 Cycle sequencing Kit with Modified Sanger Sequencing Method by ABI Genetic Analyzer 3130 (Chinikar et al. 2013a). Sequence alignments and phylogeny analysis To compare the topology of phylogenetic tree in different regions of S-segment ge- nome of CCHFV, full-length sequences of 6 CCHFV genome from Iranian patients in this study (Khorasan-e-Razavi 72, Kerman 43, Zahedan 19, Tehran 65, Isfahan 78 and Gilan 69), in addition to sixteen full length se- quences available from GenBank at www. ncbi.nih.gov used as role model and com- pared with phylogenetic trees based on re- gions of S-segment genome including nu- cleotides 1–200, 600–800 and 1000–1200 (Table 1). All sequences were aligned with the CLUSTAL W software program version http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 12–23 S Chinikar et al.: Assessment of Recombination … 14 1.6. Phylogenetic tree was generated by the Neighbor-Joining (NJ) with Kimura 2-pa- rameter distance using Mega 5 software. Bootstrap confidence limits were based on 1000 replicates. Similarity plotting and Bootscaning SimPlot software calculates and plots the percent identity of the query sequence to a panel of reference sequences and identifies the potential recombination among analyzed sequences in a sliding window, which is moved across the alignment in steps. To ob- tain evidence of recombination, a similarity plotting based on isolated sequences in this study (Khorasan-e-Razavi 72, Kerman 43, Zahedan 19, Tehran 65, Isfahan 78 and Gilan 69) was performed by SimPlot version 3.5.1 software. A boot-scanning analysis in SimPlot version 3.5.1 was undertaken to determine potential recombination between query genes and CCHFV sequences from other parts of the world. Both Simplot and bootscanning analysis were performed with Kimura (2- parameter), and window and step sizes of 200 and 10 bp (Lole et al. 1999, Worobey and Holmes 1999). Results The phylogenetic tree based on full length sequence created using the NJ algorithm showed the same topology and strongly sup- ported the same phylogenetic groups as the MJ tree (data not shown). The phylogenetic tree based on NJ demonstrated that Tehran 65, Kerman 43 and KhRazavi 72 sequences grouped within clade-IV (Asia-1). The Zahedan 19 sequence fell in clade-IV (Asia-2) and showed 82 % similarity to the isolated strain NIV112143 (JN572089). Meanwhile, Gilan 69 and Isfahan 78 sequences fell within clade-V (Europe) (Fig. 1A). The phylogenetic trees based on region of nucleotides 1–200 showed a completely dif- ferent topology in comparison to full length sequence. Firstly, ex-Afghanistan and KhRazavi 72 strains formed a separate group together. Besides, Baghdad-12 constructed an out group and did not fall in clade-IV (Asia-1). Secondly, in contrast to full length sequence of Oman strain which located in clade-IV (Asia-1), the initial region of Oman strain (nt 1–200) formed an out group in clade-IV (Asia-2). In clade V, Turkey strain (200310849) was an out group with 99 % bootstrap values (Fig. 1B). Given the phylogenetic tree of S-segment fragment including nucleotides 600–800, out groups were different. In this region, ex-Af- ghanistan and KhRazavi 72 strains formed an out group within clade-IV (Asia-2), while full length sequence of these isolates belong to clade-IV (Asia-1). In contrast to Baghdad- 12 strain which was remained as out group, Oman isolate fell in its original clade-IV (Asia-1) (Fig. 1C). In phylogenetic tree based on the region from 1000–1200 nt, the out group was China (79121) and Dubai (616). In this region, isolate Gilan69 was an out group within clade V (Fig. 1D). Evidence for recombination By examining the points at which the similarities between query (Zahedan 19, Kerman 43, Tehran 65, Gilan 69, KhRazavi 72 and Isfahan 78) and reference sequences increased or decreased, we could tentatively identify recombination breakpoints along the S-segment of CCHFV (results for similarity plotting are not shown). Figures show a boot-scan analysis on a window size of 200 nucleotides and moving in steps of 10 nucleotides along the align- ment. Boot-scanning analysis of Zahedan 19 provided evidence for viral recombination. Thus, we speculate that the recently isolated Zahedan19 suggesting a mosaic pattern comprising CCHFV isolates Dubai (Dubai 616) and India (NIV 112143) (Fig. 2A). Moreover, evidence of recombination was http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 12–23 S Chinikar et al.: Assessment of Recombination … 15 observed in Kerman 43 isolate. In this re- gard, recombination was occurred between isolated CCHFV strains from Pakistan, Matin and SR3 (Fig. 2B). Boot-scanning analysis of Tehran65 showed different cross-over site between Pakistani CCHFV strains Matin and SR3. Moreover, two possible cross-over sites at nucleotide positions 600 and 630 were observed be- tween Matin strain and Baghdad 12 strain, while Tehran 65 was considered as query sequence (Fig. 2C). Analysis for Gilan 69 was performed to confirm mosaicism and determine recombi- nation breakpoints. In this regard, bootscanning of the Gilan 69 sequence with two Russian sequences (Kashmanov and Drosdov) and two sequences from Turkey (200310849 and Kelkit 06) revealed the presence of several points of crossover (Fig. 2D). The location of the recombination event in isolate Khorasan-e-Razavi 72 was deter- mined by using the Simplot program and bootscanning analysis. Isolate Khorasan-e- Razavi 72 was compared with four repre- sentative isolates of CCHFV (SR3, Oman, Baghdad 12 and ex-Afghanistan) on the full genome. Comparison of these data indicates that some fragments of the query CCHFV genome (Khorasan-e-Razavi 72) are prone to be replaced by representative isolates, as a result of recombination (Fig. 2E). The result of recombination analysis, which is depicted in Fig. 2F, shows that Is- fahan 78 is a mosaic between Kosovo (1917), Russia (Drosdov) and Turkey (200310849) sequences, and subsequently evidence of re- combination was detected. http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 12–23 S Chinikar et al.: Assessment of Recombination … 16 Fig. 1. Continued… http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 12–23 S Chinikar et al.: Assessment of Recombination … 17 Fig. 1. A) Phylogenetic relationships between Iranian isolated sequences: Zahedan 19, Kerman 43, Tehran 65, Gilan 69, Khorasan-e-Razavi 72 and Isfahan78, and respective representative CCHFV from GenBank based on full length sequence. B) Phylogenetic tree based on nucleotides 1–200. C) Phylogenetic tree based on nucleotides 600–800. D) Phylogenetic tree based on nucleotides 1000–1200. The virus strain and the geographic origin are given for each isolate. The sequences obtained from this study are shown by asterisk symbol. The numbers above the branches indicate the bootstrap values in percentages (of 100 replicates) Table 1. Crimean Congo Hemorrhagic Fever virus isolates used in this study, with associated countries of origin, collection date and GenBank accession numbers Code of virus isolate Location Year of isolation Clade No/ Name (S segment) GenBank accession no S segment Iran-Gilan69 Iran (Northwest) 2012 V/EUR KJ027521 Iran-Isfahan78 Iran (Central) 2013 V/EUR KJ027522 Iran-Kerman43 Iran (Southeast) 2013 IV/ASI-1 KJ196326 Iran-KhRazavi72 Iran (Northeast) 2012 IV/ASI-1 KJ485700 Iran-Tehran65 Iran (North) 2011 IV/ASI-1 KJ566219 Iran-Zahedan19 Iran (Southeast) 2012 IV/ASI-2 KJ676542 NIV 112143 India 2011 IV/ASI-2 JN572089 SCT ex Afghanistan Afghanistan 2012 IV/ASI-1 JX908640 Baghadad12 Iraq 1979 IV/ASI-1 AJ538196 China79121 China 1979 IV/ASI-2 AF358784 Kosovo 1917 Kosovo 2009 V/EUR JN173797 Matin Pakistan 1976 IV/ASI-1 AF527810 Fig. 1. Continued… http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 12–23 S Chinikar et al.: Assessment of Recombination … 18 PakistanSR3 Pakistan 2000 IV/ASI-1 AJ538198 Drosdov Russia 1967 V/EUR DQ211643 Kashmanov Russia 1967 V/EUR DQ211644 ArD15786 Senegal 1972 I/W.AFR-1 DQ211640 Hodzha Uzbekistan 1967 IV/ASI-2 AY223475 AP92 Greece 1975 VI/GREECE DQ211638 Oman Oman 1997 IV/ASI-1 DQ211645 Dubai 616 Dubai 1979 IV/ASI-2 JN108025 Kelkit06 Turkey 2006 V/EUR GQ337053 200310849 Turkey 2003 V/EUR DQ211649 Table 1. Continued… Position Window: 200 bp, Step: 10 bp, GapStrip: On, Reps: 100, Kimura (2-parameter), T/t: 2.0, Neighbor-Joining Position Window: 200 bp, Step: 10 bp, GapStrip: On, Reps: 100, Kimura (2-parameter), T/t: 2.0, Neighbor-Joining BootScan-Query: Iran (Kerman 43)* BootScan-Query: Iran (Zahedan 19)* % o f P er m u te d T re es % o f P er m u te d T re es http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 12–23 S Chinikar et al.: Assessment of Recombination … 19 Fig. 2. Continued… Position Window: 200 bp, Step: 10 bp, GapStrip: On, Reps: 100, Kimura (2-parameter), T/t: 2.0, Neighbor-Joining Position Window: 200 bp, Step: 10 bp, GapStrip: On, Reps: 100, Kimura (2-parameter), T/t: 2.0, Neighbor-Joining BootScan-Query: Iran (Gilan 69)* BootScan-Query: Iran (Tehran 65)* % o f P er m u te d T re es % o f P er m u te d T re es http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 12–23 S Chinikar et al.: Assessment of Recombination … 20 Fig. 2. The location of the recombination event in Iranian isolates as query: Zahedan 19(A), Kerman 43(B), Tehran 65(C), Gilan 69(D), Khorasan-e-Razavi 72(E) and Isfahan 78(F) was determined by using the bootscanning analysis. Bootscaning was conducted with Simplot 3.5.1 software. The window size was 200 bp, with a step size 10 bp Fig. 2. Continued… BootScan-Query: Iran (KhRazavi 72)* BootScan-Query: Iran (KhRazavi 78)* Position Window: 200 bp, Step: 10 bp, GapStrip: On, Reps: 100, Kimura (2-parameter), T/t: 2.0, Neighbor-Joining Position Window: 200 bp, Step: 10 bp, GapStrip: On, Reps: 100, Kimura (2-parameter), T/t: 2.0, Neighbor-Joining % o f P er m u te d T re es % o f P er m u te d T re es http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 12–23 S Chinikar et al.: Assessment of Recombination … 21 Discussion Given that Major Histocompatibility Com- plex (MHC) can appear virus particles with molecular weight more than 10 kDa on the sur- face of host cells and encoding nucleopro- teins by S-segment has approximately 53 kDa weight, these appeared proteins recognize as the predominant antigen inducing a high immune response (Garcia et al. 2006). Thus, recombination in S-segment may provide new antigenic epitopes to the virus provided that recombination occurred in the coding regions of nucleoproteins and subsequently may deceive immune response (Alcami and Koszinowski 2000, Weber and Elliott 2002). Recombination is being recognized in- creasingly as a potentially important reason of creating and shaping genetic diversity in RNA viruses (Worobey and Holmes 1999). This aspect of genomic analysis is rarely in- vestigated for CCHFV and there is only one previously solid documented record of re- combination in CCHFV (Lukashev 2005). The small (S) segment of CCHFV (ap- proximately 1.7 kb) has a single open read- ing frame encoding the nucleocapsid (N) pro- tein. This protein plays a role in encapsidating the viral RNA to form ribonucleoprotein com- plexes (RNP). The N protein is also involved in a range of interactions with other mole- cules including viral RNA, viral polymerase, other viral proteins, host proteins as well as forming multimers with themselves (Han and Rayner 2011). Through comparison of the phylogenetic trees based on full length sequences and par- tial fragments in the S-segment genome of CCHFV, genetic switch was evident, due to recombination event. The phylogenetic re- sults obtained in this study are in accordance with Lukashev, who reported that phyloge- netic grouping according to short fragments of S-segment genome are not usually reliable and the possibility of recombination should be taken into consideration (Lukashev 2005). To obtain more solid evidence of recom- bination in this study, Bootscan analysis was used by SimPlot software. As a result, evi- dence of multiple recombination events was detected in query isolates. By ignoring the presence of recombination in S-segment ge- nome of CCHFV, phylogenetic data based on partial sequences cannot be reliable due to misinterpretation. Thus, this study strong- ly recommends using full sequence or at least long length sequences for future mo- lecular epidemiological studies to obtain more precise phylogenetic trees. Conclusion This report provides evidence of recom- bination in the S-segment of CCHFV. Switch of different genomic regions between different strains by recombination could contribute to CCHFV diversification and evo- lution. Thus, recombination may impact on investigations of virus taxonomy, phyloge- netic investigations and vaccine design. This is of utmost important to countries where multiple CCHFV variants are circulating. Acknowledgements This paper is based on the results of corresponding author's PhD thesis. We thank all members of the Arboviruses and Viral Hemorrhagic Fevers Laboratory (National Reference Lab), Pasteur Institute of Iran, and also Keyhan Azadmanesh, and Ms Niknam and Mohsen Chiani for technical assistance. This research has been financially supported by Pasteur Institute of Iran. The authors declare that there is no conflict of interests. References Alcami A, Koszinowski UH (2000) Viral mech- http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 12–23 S Chinikar et al.: Assessment of Recombination … 22 anisms of immune evasion. Immunolo- gy Today. 21: 447–455. Champour M, Chinikar S, Mohammadi G, Razmi G, Shah-Hosseini N, Khakifirouz S, Mostafavi E, Jalali T (2014a) Molecular epidemiology of Crimean-Congo hemorrhagic fever virus detected from ticks of one humped camels (Camelus dromedarius) pop- ulation in northeastern Iran. J Parasit Dis. 1–6. Doi: 10.1007/s12639-014- 0458-y Champour M, Mohammadi GR, Chinikar S, Razmi G, Shah-Hosseini N, Khakifirouz S, Mostafavi E, Jalali T (2014b) Seroepidemiology of Crimean- Congo hemorrhagic fever virus in one- humped camels (Camelus dro- medarius) population in northeast of Iran. J Vector Borne Dis. 51: 62–65. Chare ER, Gould EA, Holmes EC (2003) Phylogenetic analysis reveals a low rate of homologous recombination in negative-sense RNA viruses. J Gen Virol. 84: 2691–2703. Chinikar S, Shah-Hosseini N, Bouzari S, Jalali T, Shokrgozar MA, Mostafavi E (2013a) New circulating genomic variant of Crimean-Congo hemorrhag- ic fever virus in Iran. Archives Virol. 158: 1085–1088. Chinikar S, Shayesteh M, Khakifirouz S, Jalali T, Rasi Varaie F S, Rafigh M, Mostafavi E, Shah-Hosseini N (2013b) Nosocomial infection of Crimean- Congo haemorrhagic fever in eastern Iran: case report. Travel Med Infect Dis. 11: 252–255. Chinikar S, Shah-Hosseini N, Khakifirouz S, Rafigh M, Hasanzehi (2012a) Se- rological and molecular evaluation of Crimean-Congo Haemorrhagic Fever in Iranian probable patients. Int J Infect Dis. 16: 250. Chinikar S, Mojataba Ghiasi S, Moradi M, Goya M, Reza Shirzadi M, Zeinali M, Mostafavi E, Pourahmad M, Haeri A (2010) Phylogenetic analysis in a recent controlled outbreak of Crimean-Congo haemorrhagic fever in the south of Iran, December 2008. Euro Surveill, 15: 47. Chinikar S, Persson S M, Johansson M, Bladh L, Goya M, Houshmand B, Mirazimi A, Plyusnin A, Lundkvist A, Nilsson M (2004) Genetic analysis of crimean-congo hemorrhagic fever vi- rus in Iran. J Med Virol. 73: 404–411. Deyde VM, Khristova ML, Rollin PE, Ksiazek TG, Nichol ST (2006) Cri- mean-Congo hemorrhagic fever virus genomics and global diversity. J Virol. 80: 8834–8842. Elevli M, Ozkul AA, Civilibal M, Midilli K, Gargili A, Duru NS (2010) A newly identified Crimean-Congo hemorrhagic fever virus strain in Turkey. Int J Infect Dis. 14: 213–216. Flick R, Flick K, Feldmann H, Elgh F (2003) Reverse genetics for Crimean-Congo hemorrhagic fever virus. J Virol. 77: 5997–6006. Garcia S, Chinikar S, Coudrier D, Billecocq A, Hooshmand B, Crance J, Garin D, Bouloy M (2006) Evaluation of a Crimean-Congo hemorrhagic fever vi- rus recombinant antigen expressed by Semliki Forest suicide virus for IgM and IgG antibody detection in human and animal sera collected in Iran. J Clin Virol. 35: 154–159. Han N, Rayner S (2011) Epidemiology and mutational analysis of global strains of Crimean-Congo haemorrhagic fever virus. Virol Sin. 26: 229–244. Hewson R, Gmyl A, Gmyl L, Smirnova S E, Karganova G, Jamil B, Hasan R, Chamberlain J, Clegg C (2004) Evi- dence of segment reassortment in Crimean-Congo haemorrhagic fever virus. J Gen Virol. 85: 3059–3070. Kinsella E, Martin SG, Grolla A, Czub M, Feldmann H, Flick R (2004) Sequence http://jad.tums.ac.ir Published Online: June 27, 2015 J Arthropod-Borne Dis, March 2016, 10(1): 12–23 S Chinikar et al.: Assessment of Recombination … 23 determination of the Crimean-Congo hemorrhagic fever virus L segment. Virol. 321: 23–28. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC (1999) Full-length human immuno- deficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 73: 152–160. Lukashev AN (2005) Evidence for recom- bination in Crimean-Congo hemorrhagic fever virus. J Gen Virol. 86: 2333– 2338. Nuttall P A (2001) Crimean congo hemor- rhagic Fever. In: Service M W (ed.) Encyclopedia of arthropod-transmitted infections. UK: CAB International Wallingford Mehravaran A, Moradi M, Telmadarraiy Z, Mostafavi E, Moradi AR, Khaki-firouz S, Shah-Hosseini N, Varaie FSR, Jalali T, Hekmat S (2013) Molecular detec- tion of Crimean-Congo haemorrhagic fever (CCHF) virus in ticks from south- eastern Iran. Ticks and Tick-borne Dis. 4: 35–38. Swanepoel R, Gill D, Shepherd A, Leman P, Mynhardt J, Harvey S (1989) The clinical pathology of Crimean-Congo hemorrhagic fever. Review of Infec- tious Diseases, 11, S794-S800. Rodriguez LL, Maupin GO, Ksiazek TG, Rollin PE, Khan AS, Schwartz TF, Lofts RS, Smith JF, Noor AM, Peters CJ (1997) Molecular investigation of a multisource outbreak of Crimean- Congo hemorrhagic fever in the United Arab Emirates. Am J Trop Med Hyg. 57: 512–518. Weber F, Elliott RM (2002) Antigenic drift, antigenic shift and interferon antago- nists: how bunyaviruses counteract the immune system. Virus research. 88: 129–136. Worobey M, Holmes EC (1999) Evolution- ary aspects of recombination in RNA viruses. J Gen Virol. 80: 2535–2543. Yadav PD, Cherian SS, Zawar D, Kokate P, Gunjikar R, Jadhav S, Mishra AC, Mourya DT (2012) Genetic charac- terization and molecular clock anal- yses of the Crimean-Congo Hemor- rhagic Fever virus from human and ticks in India, 2010–2011. Infect Genet Evol. 14: 223–231. Yashina L, Vyshemirskii O, Seregin S, Petrova I, Samokhvalov E, Lvov D, Gutorov V, Kuzina I, Tyunnikov G, Tang YW (2003) Genetic analysis of Crimean-Congo hemorrhagic fever virus in Russia. J Clin Microbiol. 41: 860–862. http://jad.tums.ac.ir Published Online: June 27, 2015