- Open Access
Phylogenetic analysis of eight sudanese camel contagious ecthyma viruses based on B2L gene sequence
Virology Journalvolume 12, Article number: 124 (2015)
Camel contagious ecthyma (CCE) is an important viral disease of camelids caused by a poxvirus of the genus parapoxvirus (PPV) of the family Poxviridae. The disease has been reported in west and east of the Sudan causing economical losses. However, the PPVs that cause the disease in camels of the Sudan have not yet subjected to genetic characterization. At present, the PPV that cause CCE cannot be properly classified because only few isolates that have been genetically analyzed.
Methods and results
PCR was used to amplify the B2L gene of the PPV directly from clinical specimens collected from dromedary camels affected with contagious ecthyma in the Sudan between 1993 and 2013. PCR products were sequenced and subjected to genetic analysis. The results provided evidence for close relationships and genetic variation of the camel PPV (CPPV) represented by the circulation of both Pseudocowpox virus (PCPV) and Orf virus (ORFV) strains among dromedary camels in the Sudan. Based on the B2L gene sequence the available CPPV isolates can be divided into two genetic clades or lineages; the Asian lineage represented by isolates from Saudi Arabia, Bahrain and India and the African lineage comprising isolates from the Sudan.
The camel parapoxvirus is genetically diverse involving predominantly viruses close to PCPV in addition to ORFVs, and can be divided into two genetically distant lineages. Based on sequences of the B2L gene it is not possible to suggest that the viruses that cause CCE form a monophylogenetic group or species within the PPV phylogeny.
Camel contagious ecthyma (CCE), also known as Auzdik, Orf in camels and Pustular dermatitis is a contagious skin disease of camelids caused by a pox virus of the genus parapoxvirus (PPV), subfamily Chordopoxvirinae of the family Poxviridae. The disease has a worldwide distribution and has been reported in Mongolia , Kenya , Kazakhstan and Turkmenistan , Somalia , Sudan [5–7], Libya , Saudi Arabia [9, 10], Bahrain  and India . The disease is endemic in affected areas with variations in intensity of infection, morbidity and mortality rates and tends to occur annually in the rainy season, as the situation in Sudan and Saudi Arabia, affecting mostly young animals. The age group at risk are those less than one year of age  including cases reported in month-old camel calves in Saudi Arabia . The marked seasonality associated with the rainy season may be due to the optimum condition for the survival and perpetuation of the virus and skin abrasions caused by browsing thorny trees . In most cases, the disease caused no mortality, but when camel calves are severely affected the pox-lesions interferes with the calves’ ability to suckle or graze and extends to eyelids leading to blindness, particularly at the Savanah belt in the Sudan, leading to mortality rates that can reach 9 % . Clinically, the pox-lesions first appeared on the lips of affected animals as small papules that progressively developed into scabs on the lips, muzzle, nares and eyelids culminating into fissured crusts on the lips. Swelling of the head and sometimes the neck have been observed in the field. The lesion is proliferative and highly vascularized and may extend into gum, palate and tongue. The major factors associated with increased likelihood of CCE occurrence are season of the year, camel age, camel movements and location and their association with thorny trees .
Camel contagious ecthyma is a sparsely studied disease and the causative virus have only recently genetically characterized. According to Abubakr et al  sequence homologies and phylogenetic analysis of the major envelope gene (B2L) of 2 CCE skin scabs (two from Bahrain and one from Saudi Arabia) showed that these viruses are closely related to the pseudocowpox virus (PCPV) species of the PPV genus of the family poxviridae. In addition, Nagarajan et al  showed that CCE in Indian dromedary camels is also caused by a PCPV.
Currently, there are four established species within the PPV genus of the family Poxviridae; Orf virus (ORFV), the type species of the genus, which causes disease mainly in sheep and goats, PCPV and bovine papular stomatitis virus (BPSV) both infect cattle, in addition to Parapoxvirus of red deer in New Zealand (PVNZ), which has only been isolated from red deer in New Zealand. BPSV and PCPV affect mainly cattle, but differ from ORFVs in the site of the pox lesion, as BPSV is restricted to the muzzle and PCPV the teat. Tentative species include the camel PPV (CPPV), reindeer parapox virus, musk ox and seal parapoxvirus [13–15]. At present, the CPPV that cause CCE cannot be properly classified because only few isolates that have been genetically analyzed.
Camel contagious ecthyma is believed to be endemic in the Sudan since long time, but was first diagnosed in the disease outbreaks west of the country in 1986 . Later on, the disease was reported in eastern Sudan [6, 7]. However, the PPVs that cause the disease in dromedary camels of the Sudan have not yet subjected to genetic characterization.
In the present study, we investigated lesion-derived DNAs from eight skin scabs collected from dromedary camels in the Sudan presented with CCE lesions and showed for the first time, that the CPPV does not form a monophyletic clade within the PPV phylogeny and the disease is caused by viruses from more than one PPV species.
In this investigation, we screened archive scab sample collected between 1993 and 2013 in addition to the new scab collection from Sudan in August 2013. The 2013 outbreak of the disease involved 17 camel herds around Showak area of eastern Sudan, causing morbidity and mortality rates of 20 and 1.5 %, respectively. Pustular dermatitis lesion characteristic for CCE were observed (Fig. 1). Out of 10 scab materials, eight gave positive specific amplification in real-time PCR and gel-based PCR (data not shown). Table 1 shows the origin of positive samples used in this study. To characterize the DNA fragments generated by PCR, the PCR products were sequenced and subjected to genetic analysis.
The B2L gene fragment from the eight DNAs has a high G + C content of 62.6 to 63.9 % (Table 2). The seven isolates (SD-V4, V8, V13, V20, K1, K2 and K3) shared 96–98 % homology with viruses in the species PCPV while one isolate (V34-SD) shared 96–99 % homology with viruses in the ORFV species of the genus PPV (Table 2). Seven CPPV Sudanese isolates shared 96–98 % nucleotide sequence identity with virus strains in PCPV species, while one isolate (V34-SD) shared only 94–95 % identity with PCPV, but 96–99 % identity with ORFV.
Results of multiple sequence alignment and per cent identity matrix based on 160 amino acids sequences (Table 3) revealed that the Sudanese CPPVs from 1993 to 2013 examined in the present study, shared variable amino acid sequence identity ranged from 95.04 to 100 % (intra-group similarity). Five isolates (V4-SD, V20-SD, K1-SD, K1-SD and K3-SD) possessed 100 % identity, isolate V13-SD has a less homology of 95.04 to 99.69 % with the rest Sudanese isolates. On the other hand, isolate V34-SD, which has the closest identity to ORFV species, shared a relatively lower per cent identity of 95.04–96.21 with the rest seven isolates (Table 3). As shown in Table 3, the previously published five CPPV isolates Judhpur-IN and Cam/09 from India, BH-1 and BH-3 from Bahrain and SA-98 from Saudi Arabia shared intra group similarity of 98.94 to 100 % homology and a 95.04 to 99.69 % homology with the eight Sudanese isolates analysed in the present study. On the other hand, isolate CE41, which was classified as BPSV shared only 81.21 to 82.55 % homology with all the CPPVs.
A phylogenetic tree was created for the amino acids sequences of the B2L gene of 40 PPV strains, including the eight strains described in the present study using the Neighbor-Joining method. The result revealed that the Sudanese CPPVs clustered with other PPVs published earlier and that the analyzed PPVs displayed distinct clusters representing BPSV and PVNZ species (Fig. 2). Members of the other two PPV species, namely ORFV and PCPV are distributed between the two main branches shown in the phylogenetic tree. Seven of the eight Sudanese CPPV isolates sequenced in the present study clustered together closer to PCPV isolates from Finland and Germany. One Sudanese isolate (V34-SD) clustered in a different sub-branch that contains ORFV isolates from Finland and Germany and a PCPV isolate from Finland. On the other hand, the five previously published sequences representing CPPV isolates from India, Saudi Arabia and Bahrain clustered at the other branch of the phylogenetic tree together with PCPV isolates from Finland and ORFV isolates from Asian countries (India, Korea, China and Taiwan). Isolate CE41 isolated from a dromedary camel in the Sudan clustered with the BPSV species of PPV.
The obtained nucleotide sequences were submitted to GenBank and assigned the accession numbers KR231664, KR231665, KR231666, KR231667, KR231668, KR231669, KR231670 and KR231671 for isolates V4-SD, V8-Sd, V13-Sd, V20-SD, V34-SD, K1-SD, K2-SD and K3-SD, respectively.
A preferred target gene to generate PCR amplified DNA fragments for sequence analysis and comparison of PPV DNA is the open reading frame (ORF) 011 (B2L gene), the PPV homologue of the vaccinia virus Copenhagen (VACV) gene F13L, which encodes the major envelope antigen p37K .
In the present study, a genetic analysis of eight Sudanese isolates of the CPPV was performed. BLAST results unveiled that seven isolates are genetically nearly related to PCPV and one is closer to ORFV. The B2L gene fragment from the eight DNAs has a high G + C content of 62.6 to 63.9 % in agreement with that of the whole PPV genome .
The findings of the present study provide evidence for a close relationship of seven Sudanese strains despite the different geographic location and year of collection. Interestingly, five isolates collected at the same area of eastern Sudan, but during a time span between 1993 and 2013 have 100 % amino acid homology pointing to a circulation of genetically stable virus and a stable genetic makeup and a wide distribution in the country.
The dominant CPPV circulating in the Sudan is found to be genetically close to PCPV species of the PPV and a genetic variability in this country, of at least two PPV species: PCPV (n = 7) and ORFV (n = 1). Three isolates, which were collected over 20 years, possessed 100 % genetic homology.
Search of the literature showed that five isolates have been antecedently classified as PCPV [11, 12]. The results of this study ascent the total number of CPPV that belong to PCPV species to 12. However, the study also identified one isolate that genetically closer to ORFV than PCPV. Besides, we also corroborated previous classification of isolate CE 41 as BPSV . Taken together, the present study shows heterogeneity among the CPPVs and suggest that clearly, these viruses do not form a monophyletic species within PPV phylogeny, at least based on the B2L gene analysis. A search of the literature revealed the presence of only one CPPV strain (Cam/09) that has been classified as ORFV (accession # GU460370.1) as it has 99 % nucleotide homology with an Indian ORFV from sheep . Here we confirm previous research work and present grounds that CCE in camels can be caused by ORFV. This is what has been anticipated before any genetic analysis was made, as lesions produced by the causative virus is characterized by a progressive nodule-vesicle-pustule-scab formation found mainly on the lips that resembles lesions induced by ORFV in sheep and goats, but not PCPV or BPSV . The findings also might points to a transmission from sheep or goats, as it is common for camels, sheep, and goats to share the same pasture in the Sudan. Contrary to such speculation, Abu Elzein et al , in Saudi Arabia failed to experimentally infect sheep and goats with an isolate of CPPV. Unfortunately, the PPV circulating in sheep and goats in the Sudan has not yet genetically analyzed. Therefore, molecular analysis of ORFV circulating in small ruminants in the same area and experimental infection of sheep and goats with a CPPV genetically closer to ORFV is expected to rule out such possibility.
Previous reports from Saudi Arabia  and India  implicated PCPV as the cause of the disease in camels. Of note, Abubakr et al  on the bases of their first connection of the CPPV to PCPV postulated a relatively recent introduction of the virus to camels, possibly from humans, that have been exposed to PCPV-infected ruminants. Yet, no PCPV infection in human or animals has been reported in the Sudan or Arabian Peninsula, but only one case of ORFV infection of a man in United Arab Emirates . Of interest is to determine the reservoir host that maintain the virus in nature during the inter-epizootic period. A serological survey involving domestic and wild species of animals in the affected area of eastern Sudan may identify such reservoir host and contribute to its control.
On the other hand, a nucleotide sequence of the B2L gene from a cell culture- adapted CPPV (CE41) was found to be genetically closely related to BPSV . This genetic placement might be caused by genomic rearrangements known to occur when PPVs are passaged in cell culture [24–26]. Nevertheless, taking into consideration the data presented in this report, it can be postulated that contagious ecthyma in dromedary camels is caused by viruses from different species of PPV. This situation is not uncommon as some PPVs obtained from wild Japanese serows are genetically closer to ORFV, and some are closer to BPSV . However, a possibility still exists that a separate virus species of PPV genetically related to both PCPV and ORFV based on the sequence of the B2L gene, but differ in other genomic domains causes this disease. To this end whole genome sequencing and comparison with known PPV species is a requisite and anticipated to bring out a better approximation to reality.
At present 11 out of the 14 published CPPV sequences are grouped along with PCPV, two are ORFV related and one belonged to the BPSV species of PPV indicating genetic variability of the CPPV. Further investigations should examine other genes of the virus and perform complete genome sequencing in order to better describe and classify the CPPV.
According to previous phylogenetic analysis based on the B2L gene of PPV, PCPV and ORFV strains are usually characterized by a certain level of variability, while PVNZ and BPSV strains may possess 100 % identity . In line with that, the phylogenetic tree generated in this study for 40 PPV sequences showed distinct clustering of BPSV and PVNZ strains, with ORFV and PCPV strains bundling in more than one subgroup mostly with a link to the geographic location of respective disease outbreaks. Similarly Oem et al  showed that some PPVs of a particular classification clustered with different species in the phylogenetic tree. Therefore, the present study supports the contention by Oem et al  that phylogenetic analysis using partial B2L gene sequences may not be suited owing to the recent increase in the publication of more sequences of PPV from different countries and hosts. Besides, it should also be kept in mind that infections with multiple PPV strains of a single species of PPV can occur in one animal as recently described .
CCE is a sparsely studied disease of camels and presently does not have to be reported to the World Organization of Animal Health (OIE) despite several publications reporting wide distribution of the disease. The disease is economically important and a recent publication from Iran  described severe cases of CCE that caused a mortality rate of 6 %, therefore, this important camel disease has to be looked at seriously.
We genetically analyzed based on sequences of the B2L gene eight strains of CPPV collected from two areas of the Sudan between 1993 and 2013. The prevailing virus is genetically close to PCPV species of the PPV and a genetic variability in this country, of at least two PPV species exists: PCPV (n = 7) and ORFV (n = 1). Based on sequences of the B2L gene it is not possible to suggest that the viruses that cause CCE form a monophylogenetic group or species within the PPV phylogeny.
Origin of specimens
This study was approved by the Institutional Review Board of the Veterinary Research Institute, Sudan. This study was approved by the Institutional Review Board of the Veterinary Research Institute, Sudan Skin scab specimens (n = 6) were collected from dromedary camels (camelus dromedarius) showing symptoms of contagious ecthyma in east and west Sudan during field epidemiological investigations between 1993 and 2005. Diagnosis of the disease was made on the bases of clinical evidence and electron microscopy and later by PCR . Details on disease symptoms and epidemiology were published . Additional four specimens were collected in 2013. Data on the clinical picture, morbidity and mortality rates resulted in the 2013 disease outbreaks in Showak area of eastern Sudan were collected. Table 1 shows some data on positive specimens analyzed in the present study.
Tissue homogenization and DNA extraction
A 20 % suspension (weight/volume) was made of the scab material in tris-EDTA (TE) buffer (pH 7.4), freeze-thawed at −30 °C, mechanically homogenized using a mechanical Homogenizer (TissueRuptor, Qiagen, Hilden, Germany) and centrifuged at 1500 × g for 10 min at 4 °C. Total viral DNA was extracted from supernatants using a QiaAmp DNA Mini Kit (Qiagen, Hilden, Germany) according to manufacturer’s instructions.
To check its PPV identity each extracted DNA was first tested using the quantitative real-time PCR assay described by Nitsche et al . The amplification was carried out on Mx3000P qPCR system (Agilent Technologies, USA) using Brilliant II QPCR Master Mix (Agilent Technologies, USA). Second, with DNA of known PPV identity, a PCR protocol was used to obtain DNA fragment for sequencing. This PCR assay targets the commonly used B2L gene (open reading frame 011) which is homologue of vaccinia virus Copenhagen (VACV) gene F13L encoding the major envelop antigen p37K . PCR amplification was done using forward primer (5′-TTAATTTATTGGCTTGCAGAACTCCGAGCGC-3′) and reverse primer (5′-ATGTGGCCGTTCTCCTCCATC-3′)  that amplify 1200 bp DNA sequence. PCR amplification of the envelope gene was performed using the following thermal profiles initial denaturation at 94 °C for 10 min, followed by 35 cycles of denaturation at 94 °C for 1 min, annealing at 55 °C for 30 s, extension at 72 °C for 1.5 min, and final extension at 72 °C for 10 and DNA was amplified with a thermal cycler TP3000 (Biometra, Germany). The PCR products were then checked in agarose gel, purified using QIAquick PCR Purification Kit (Qiagen, Hilden, Germany) and then sent for sequencing. Sequencing was completed using the BigDye® Terminator v3.1 cycle sequencing kit chemistry and each primer pairs. Nucleotide positions were confirmed based on two independent sequencing reactions in both directions.
BLAST and phylogenetic analysis
The terminal unreliable nucleotides of the obtained B2L gene sequences were first trimmed-off automatically and both forward and reverse sequences aligned in the Geneious® 8.1.4 package. The biologically correct sequences of the eight CPPV isolates were subjected to basic local alignment search was compared with nucleotide sequence in the GenBank database using the online BLASTN program on the NCBI website .
Phylogenetic tree based on amino acid sequences was constructed for 40 PPV strains (Table 4), including sequences from the four PPV species; ORFV (n = 9), PCPV (n = 10), BPSV (n = 6) and PVNZ (n = 3) and eight strains from this study, using the Neighbor-Joining method in MEGA6 . The amino acid substitution model used was the Johen-Tylor-Thornton model . The significance of all deduced phylogenetic trees was verified by bootstrap analysis of 1000 replicates.
DashtserenTs SF, Varejka F, Khokho A. Camel contagious ecthyma (Pustular dermatitis). Acta Virol. 1984;28:122–7.
Munz E, Schillinger D, Reimann M, Mahnel H. Electron microscopical diagnosis of ecthyma contagiosum in camels (camelus dromedarius). First report of the disease in Kenya. J Vet Medicine B. 1986;33:73–7.
Buchnev KN, Tulebaev SZ, Sansyzbaev AB. Infectious diseases of camels in the USSR. Revue Sci Tech Off Int Epi. 1987;6(2):487–95.
Moallin ASM, Zessin KH. Outbreak of camel contagious ecthyma in central Somalia. Trop Anim Health Prod. 1988;20:185–6.
Ali AO, Kheir AM, Abdulamir H, Barri MES. Camel contagious ecthyma in Sudan. A case report. Rev Elev Med Vet Pays Trop. 1991;44:143–5.
Khalafalla AI, Mohamed MEH. Epizootiology of camel contagious ecthyma in eastern Sudan. Revue Elev Vet Pays Trop. 1997;50:99–103.
Khalafalla AI. Camel contagious ecthyma: risks in young calves. Revue Elev Vet Pays Trop. 2000;53(2):173–6.
Azwai SM, Carter SD, Woldehiwet Z. Immune responses of the camel (Camelus dromedarius) to contagious ecthyma (Orf) virus infection. Vet Microbiol. 1995;47(1 –2):119–31.
Abu Elzein EME, Housawi FMT, Al-Afaleq AI, Ramadan RO, Gameel AA, Al-Gundi O. Clinicopathological response of dromedary camels and sheep to cross-experimental infection with two virulent orf viruses originating from camels and sheep. J Camel Pract Res. 2004;11(No. 1):15–9.
Housawi FM, Abu Elzein EME, Gameel A, Mustafa M, Al-Afaleq A, Gilary J, et al. Severe Auzdyk infection in one-month old camel calves (camelus dromedarius). Vet Archives. 2004;74:467–74.
Abubakr MI, Abu-Elzein EME, Housawi FM, Abdelrahman AO, Fadlallah ME, Nayel MN, et al. Pseudo cowpox virus: The etiological agent of contagious ecthyma (Auzdyk) in camels (Camelus dromedarius) in the Arabian Peninsula. Vector Borne Zoonotic Dis. 2007;7(2):257–60.
Nagarajan G, Ghorui SK, Kumar S, Pathak KM. Complete nucleotide sequence of the envelope gene of pseudocowpox virus isolates from Indian dromedaries (Camelus dromedarius). Arch Virol. 2010;155(10):1725–8.
Mercer A, Fleming S, Robinson A, Nettleton P, Reid H. Molecular genetic analyses of parapoxviruses pathogenic for humans. Arch Virol Suppl. 1997;13:25–34.
Becher P, Koenig M, Mueller G, Siebert U, Thiel H-J. Characterization of sealpox virus, a separate member of the parapoxviruses. Arch Virol. 2002;147:1133–40.
Vikøren T, Lillehaug A, Akerstedt J, Bretten T, Haugum M, Tryland M. A severe outbreak of contagious ecthyma (orf) in a free-ranging musk ox (Ovibosm oschatus) population in Norway. Vet Microbiol. 2008;127:10–20.
Saitou N, Nei M. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol. 1987;4:406–25.
Zuckerkandl E, Pauling L. Evolutionary divergence and convergence in proteins. In: Bryson V, Vogel HJ, editors. Evolving Genes and Proteins. New York: Academic; 1965. p. 97–166.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Sullivan JT, Mercer AA, Fleming SB, Robinson AJ. Identification and characterization of an orf virus homologue of the vaccinia virus gene encoding the major envelope protein p37K. Virology. 1994;202:471–8.
Delhon G, Tulman ER, Afonso CL, Lu Z, De la Concha-Bermejillo A, Lehmkuhl HD, et al. Genomes of the parapoxviruses ORF virus and bovine papular stomatitis virus. J Virol. 2004;78:168–77.
Dal Pozzo F, Martinelle L, Gallina L, Mast J, Sarradin P, Thir E, et al. Original findings associated with two cases of bovine papular stomatitis. J Clin Microbiol. 2011;49:4397–400.
Li H, Zhu X, Zheng Y, Wang S, Liu Z, Dou Y, et al. Phylogenetic analysis of two Chinese orf virus isolates based on sequences of B2L and VIR genes. Arch Virol. 2013;158(7):1477–85.
Al-Salam S, Nowotny N, Sohail MR, Kolodziejek J, Berger TG. Ecthyma contagiosum (orf)—report of a human case from the United Arab Emirates and review of the literature. J Cutan Pathol. 2008;35(6):603–7.
Cottone R, Büttner M, Bauer B, Henkel M, Hettich E, Rziha HJ. Analysis of genomic rearrangement and subsequent gene deletion of the attenuated Orf virus strain D1701. Virus Res. 1998;56:53–67.
Fleming SB, Lyttle DJ, Sullivan JT, Mercer AA, Robinson AJ. Genomic analysis of a transposition-deletion variant of orf virus reveals a 3.3 kbp region of non-essential DNA. J Gen Virol. 1995;76:2969–78.
McInnes CJ, Wood AR, Nettleton PE, Gilray JA. Genomic comparison of an avirulent strain of Orf virus with that of a virulent wild type isolate reveals that the Orf virus G2L gene is non-essential for replication. Virus Genes. 2001;22:141–50.
Inoshima Y, Murakami K, Yokoyama T, Sentsui H. Genetic heterogeneity among parapoxviruses isolated from sheep, cattle and Japanese serows (Capricornis crispus). J Gen Virol. 2001;82:1215–20.
Oem JK, Lee EY, Lee KK, Kim SH, Lee MH, Hyun BH. Bovine papular stomatitis virus (BPSV) infections in Korean native cattle. J Vet Med Sci. 2013;75(5):675–8.
Huang T, Tulman ER, Diel DG, Khatiwada S, Sims W, Edwards JF, et al. Coinfection with multiple strains of bovine papular stomatitis virus. Arch Virol. 2015;160(6):1527–32. doi:10.1007/s00705-015-2394-2.
Mombeni EG, Mombeini MG, Varshovi HR, Khalaj M, Kenarkohi M, Goudarzi M, et al. Outbreak of contagious ecthyma in camels (Camelus dromedarius and Camelus bactrianus) in Southwest Iran. Rev Elev Med Vet Pays Trop. 2013;66(4):113–5.
Khalafalla AI, Büttner M, Rziha HJ. Polymerase chain reaction (PCR) for rapid diagnosis and differentiation of para-and orthopox virus infections in camels. In: H.P.S.Makkar and G.J.Viljoen (eds). Applications of Gene-Based Technologies for Improving Animal Production and Health in Developing Countries. FAO/IAEA Publications, SPRINGER; 2005: p.335–42.
Nitsche A, Büttner M, Wilhelm S, Pauli G, Meyer H. Real-time PCR detection of Parapoxvirus DNA. Clin Chem. 2005;52:316–9.
Guo J, Rasmussen J, Wünschmann A, de La Concha-Bermejillo A. Genetic characterization of orf viruses isolated fromvarious ruminant species of a zoo. Vet Microbiol. 2004;99(2):81–92.
Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.
Tikkanen MK. Recent isolates of parapoxvirus of Finnish reindeer (Rangifer tarandus tarandus) are closely related to bovine pseudocowpox virus. J Gen Virol. 2004;85:1413–8.
Hosamani M, Bhanuprakash V, Scagliarini A, Singh RK. Comparative sequence analysis of major envelope protein gene (B2L) of Indian orf viruses isolated from sheep and goats. Vet Microbiol. 2006;116:317–24.
Chan KW, Hsu WL, Wang CY, Yang CH, Lin FY, Chulakasian S, et al. Differential diagnosis of orf viruses by a single-step PCR. J Virol Methods. 2009;160:85–9.
Friederichs S, Krebs S, Blum H, Wolf E, Lang H, von Buttlar H, et al. Comparative and retrospective molecular analysis of Parapoxvirus (PPV) isolates. Virus Res. 2014;181:11–21.
Bora DP, Barman NN, Das SK, Bhanuprakash V, Yogisharadhya R, Venkatesan G, et al. Identification and phylogenetic analysis of orf viruses isolated from outbreaks in goats of Assam, a northeastern state of India. Virus Genes. 2012;45(1):98–104.
Zhang K, Liu Y, Kong H, Shang Y, Liu X. Comparison and phylogenetic analysis based on the B2L gene of orf virus from goats and sheep in China during 2009-2011. Arch Virol. 2014;159(6):1475–9.
Lojkic I, Cac Z, Beck A, Bedekovic T, Cvetnic Z, Sostaric B. Phylogenetic analysis of Croatian orf viruses isolated from sheep and goats. Virol J. 2010;7(1):314.
Hautaniemi M, Vaccari F, Scacliarini A, Laaksonen S, Huovilainen A, McInnes CJ. Analysis of deletion within the reindeer pseudocowpoxvirus genome. Virus Res. 2011;160:326–32.
Gassmann U, Wyler R, Wittek R. Analysis of parapoxvirus genomes. Arch Virol. 1985;83:17–31.
Scagliarini A, Vaccari F, Turrin F, Bianchi A, Cordioli P, Lavazza A. Para-poxvirus infections of Red Deer, Italy. Emerg Infect Dis. 2011;17:684–7.
Jones DT, Taylor WR, Thornton JM. The rapid generation of mutation data matrices from protein sequences. Comput Appl Biosci. 1992;8:275–82.
This research was partially supported by a research grant no 140105 from the Deanship of Scientific Research, King Faisal University, Kingdom of Saudi Arabia. We thank Prof. Marzook M Al-Eknah, Director Camel Research Center for the continuous support.
The authors declare they have no competing interests.
AIK and YHA initiated and participated in the planning of the study and collected the field specimens. IME, KAA and AIA participated in the planning of the study and in the lab work. AIK wrote the first draft of the manuscript. All authors contributed to data analysis and overseen manuscript revision.