Hassana Kyari Mangga1,2,3*
Jamilu Abubakar Bala1,4
Krishnan Nair Balakrishnan1Alhaji Modu Bukar1,5Zaharaddeen Lawan6
Auwal Gambo1,7
Faez Firdaus Abdullah Jesse8Mustapha M. Noordin2
Mohd-Lila Mohd-Azmi1*- 1Virology Unit, Department of Pathology and Microbiology, Faculty of Veterinary Medicine, Universiti Putra Malaysia, Serdang, Malaysia
- 2Department of Veterinary Pathology and Microbiology, Faculty of Veterinary Medicine, Universiti Putra Malaysia, Serdang, Malaysia
- 3Department of Microbiology, Faculty of Science, University of Maiduguri, Maiduguri, Nigeria
- 4Department of Medical Laboratory Sciences, Faculty of Allied Health Sciences, Bayero University Kano, Kano, Nigeria
- 5Department of Science Laboratory Technology, Ramat Polytechnic Maiduguri, Maiduguri, Nigeria
- 6Department of Agricultural Technology, College of Agriculture, Hussaini Adamu Federal Polytechnic, Kazaure, Nigeria
- 7Department of Microbiology, Faculty of Science, Usmanu Danfodiyo University Sokoto, Sokoto, Nigeria
- 8Department of Veterinary Clinical Studies, Faculty of Veterinary Medicine, Universiti Putra Malaysia, Serdang, Malaysia
Contagious ecthyma commonly known as Orf is a globally important, highly contagious zoonotic, transboundary disease that affects domestic and wild ruminants. The disease is of great economic significance causing an immense impact on animal health, welfare, productivity, and trade. Detailed analysis of the viral genome is crucial to further elucidate the molecular mechanism of Orf virus (ORFV) pathogenesis. In the present study, a confluent monolayer of lamb testicle cells was infected with the processed scab sample obtained from an infected goat. The presence of the virus was confirmed using polymerase chain reaction and electron microscopy, while its genome was sequenced using next-generation sequencing technology. The genome sequence of Malaysian ORFV strain UPM/HSN-20 was found to contain 132,124 bp with a G + C content of 63.7%. The homology analysis indicates that UPM/HSN-20 has a high level of identity 97.3–99.0% with the other reference ORFV strain. Phylogenetic analysis revealed that ORFV strain UPM/HSN-20 is genetically more closely related to ORFV strain XY and NP from China. The availability of the genome-wide analysis of ORFV UPM/HSN-20 strain from Malaysia will serve as a good platform for further understanding of genetic diversity, ORFV infection, and strategic development for control measures.
Introduction
Orf virus (ORFV), Bovine papular stomatitis virus (BPSV), Pseudocowpox virus (PCPV), and Parapoxvirus of red deer in New Zealand (PVNZ), collectively known as parapoxviruses, are among the most economically significant causal agents of diseases in the livestock industry for decades (Hatcher et al., 2015; Bala et al., 2018a). ORFV is the causative agent of contagious ecthyma, or contagious pustular dermatitis virus was the first virus to be described in this group, the prototype of the genus Parapoxvirus, family Poxviridae. It is an ovoid, enveloped, double-stranded positive-sense DNA virus with a linear genome of 130–140 kb (Fleming et al., 2015; Sahu et al., 2020b). Its genome consists of two regions: the core and variable regions. The core region is conserved across most members of the poxvirus group and encodes essential proteins for structural components, nucleotide biosynthesis, genome replication, transcription, and assembly. On the other hand, the variable regions are found at the 5′ and 3′ terminal ends. The proteins encoded by these terminal genes are the major determinants of host range and viral pathogenesis (Delhon et al., 2004).
The disease has a worldwide distribution and mainly affects sheep and goats. But recently, cross-infection of other various domestic and wild ruminants has been reported in different parts of the world (Spyrou and Valiakos, 2015). The major route of transmission is via skin abrasion. After an incubation period of 8–14 days, the infected animals develop proliferative lesions commonly on the lips, muzzle, ears, eyelids, tongue, nostril, and occasionally the infection spreads to other non-wooly areas legs, feet, and udders (Chi et al., 2013). The disease is self-limiting, and lesions slowly heal within 1–2 months after progression from erythema, vesicle, pustules, papules, and scabs (Peralta et al., 2018). Epidemiological reports indicated morbidity is high, but mortality associated with CE is usually low. However, in severe cases, complications from fly infestation and secondary infection by bacteria and fungi have a significant impact leading to increased death of young animals, decrease in optimum productivity, and restriction of domestic and international trade (Bala et al., 2018b; Lawal et al., 2021). For zoonotic infection, in immunocompetent individuals, human Orf is localized and self-limiting, and disease resolves spontaneously within weeks (Hasheminasab et al., 2016; Kassa, 2021). In immunocompromised patients, complications such as lymphangitis, erythema multiforme, and toxic erythema have been reported (Lederman et al., 2007; Cubells et al., 2016).
In Malaysia, several studies have reported the occurrence of ORFV among sheep and goats from the diverse geographical region at a different time based on seroprevalence and molecular detection based on partial gene amplification (Al-Ajeeli, 1995; Abdullah et al., 2015; Sadiq et al., 2017; Bala et al., 2020). Recently, a high prevalence rate of 22.8% was reported during a serosurvey involving 504 animals in Terengganu (Bala et al., 2019). Another study reported an infection rate of 36.7 and 7.8% among goats and sheep from a small ruminant farm in Selangor (Jesse et al., 2018). To date, most of the published studies focused only on the partial sequence analysis based on conserved genes and seroprevalence studies. Despite the high incidence of the disease, there are no comprehensive molecular studies of the endemic ORFV affecting small ruminant farms in Malaysia. This is the first report on genome-wide sequence analysis of ORFV isolated in Malaysia and it is an effective approach to an enhanced understanding of the virus pathogenesis, epidemiology, and strategic development for control measures.
Materials and Methods
Sample Processing
Scab samples were collected from the tissue of an infected goat with clinical symptoms of Orf from a farm located in Selangor, Malaysia. The sample was processed according to the method described by Chi et al. (2015). One gram of scab was homogenized in 0.01 Phosphate-Buffered Saline (PBS). The suspension 10%w/v obtained was then centrifuged at 3000 rpm for 10 min. The supernatant was carefully aspirated and supplemented with penicillin (100 units/ml) and streptomycin (100 μg/ml) then stored at −80°C.
Virus Propagation
Lamb testis cells obtained from ATCC (OA3.Ts CRL-6546) were maintained in Dulbecco’s Modified Eagle’s Medium (DMEM, GIBCO BRL Grand Island, United States), supplemented with 10% fetal bovine serum (FBS, GIBCO BRL Grand Island, United States), penicillin G (100 units/ml), and streptomycin (100 μg/ml). Confluent monolayer of the cells in 25 cm3 flask was infected with 100 μl of processed virus sample. The flask was placed in a CO2 incubator at 37C for 1 h. Fresh DMEM containing 1% FBS was added to the cells after washing with PBS. Frequently, the cells are observed daily for virus-induced cytopathic effect (CPE). When the CPE is 80%, the virus is harvested by methods described by Azmi and Field (1993) and Yang et al. (2014).
Virus Concentration and Purification
After four consecutive passages of the ORFV isolate using a T25 cm3 flask, the virus was subsequently cultivated for purification. This is conducted by sucrose gradient ultracentrifugation as described by Guo et al. (2003). The viral supernatant obtained after three cycles of freezing and thawing was subjected to concentration using polyethylene glycol 6,000 (PEG 6000; Calbiochem, Darmstadt, Germany). We added 10 ml of 36% sucrose in TE buffer onto an ultracentrifuge tube, then slowly transferred the concentrated virus suspension and centrifuged at 65000 × g for 1 h at 4°C. The collected pellet was further purified using five different concentrations of sodium diatrizoate fractions at 50, 40, 30, 20, and 10% (w/v) in TE buffer. The white opalescent virus band detected was carefully aspirated from the gradient tubes and resuspended in PBS. The purified virus suspension was negatively stained with a drop of 2% phosphor-tungstic acid, followed by examination under a transmission electron microscope (TEM: Hitachi, H7100) at an accelerated voltage of 75 kV at microscopy unit at Institute of Bioscience, Universiti Putra Malaysia.
DNA Extraction
Viral DNA from scab samples, cell cultures passages, and purified virions was extracted using a highly specific viral nucleic acid extraction kit (innuPREP virus DNA kit, Analytik Jena, Germany) according to the manufacturer’s instructions. The purity and concentration of the DNA were determined spectrophotometrically using BioPhotometer™ plus (Eppendorf, Hamburg Germany) at the ratio of 260 –280 nm.
Polymerase Chain Reaction
A set of primer designated as forward, and reverse as previously described by Chi et al. (2013) was designed and synthesized by Integrated DNA Technology (forward: 5′-ATG TGG CCG TTC TCC TCT ATC-3′; Reverse: 5′-TTA ATT TAT TGG CTT GCA G-3′) for the amplification of fragments of major immunogenic envelope gene (B2L gene) of ORFV isolate. The polymerase chain reaction was carried out using Qiagen Kit (Qiagen, United States) following the manufacturer’s protocol. Briefly, PCR was performed containing 122 ng extracted DNA, 10 nmol of each primer, Taq Master Mix, and nuclease-free water. The thermal cycler was programmed based on the following cycling condition: the initial denaturation was at 94°C for 3 min, followed by the second cycle of denaturation at 94°C for 30 s; annealing temperature was 55°C for 30 s; extension at 72°C for 1 min; and the program was set to repeat from cycle 2 for 35 cycles and hold at 10°C. Gel electrophoresis of the PCR products was carried out using 1.25% agarose gel at 80 V for 35 min. The gel was stained with red safe, and DNA bands were viewed using a UV transilluminator (Gel Doc, BioRad).
Genome Sequencing and Assembling
Viral DNA extracted from the purified virions was subjected to library preparation using Nextera DNA Flex Library Prep kit (Illumina Inc., United States) according to the manufacturer’s instruction. The libraries were sequenced using MiSeq technology (Illumina Inc., United States). Raw reads were filtered to remove artificial sequences and a Phred quality score of N30 was employed for quality control in BBDuk software (BBDuk version 36). MEGAHIT software (version 1.2.8) was used for de novo assembly of the paired-end reads to contigs. The contigs were further analyzed using RAST and Prodigal PROKKA software (version 1.14.0) for gene prediction.
Open Reading Frame Analysis and Annotation
All the predicted proteins were then subjected to analysis using online Expasy software and BlastP at the National Center for Biotechnology Information1 to identify the gene family. Thus, the predicted ORFs are arranged, numbered, and named according to the reference ORFV strains. In addition, the percentage of nucleotide sequence among the isolates was determined using Clustal Omega (1). After annotation, to confirm the missing of ORF 001 and ORF 134 from the present isolate, specific primers were used for the amplification of a 928 bp region forward primer (GCCGTGGCCGAGTTGTAG) reverse primer (CTCGGTGACCTGCCTGAC) at the 5′ end and 1,172 bp region forward primer (TGCACCAGCATCTGCTAAAC) reverse primer (ACATTCCAAGCTCTCGTCGT) at the 3′ end of the genome.
Phylogenetic Tree and Distance Matrix
An alignment containing 26 complete genome sequences of parapoxviruses and other closely related Squirrelpox virus and Molluscum Contagiosum virus as an outgroup based on the conserved gene (ORF111) was used for comparison. The Phylogenetic tree was constructed using the maximum likelihood methods and Tamura-Nei model (Tamura and Nei, 1993). The reliability of the tree was assessed by a bootstrap test with 1,000 replicates, constructed with MEGA-X (Kumar et al., 2015). The identity value of the nucleotide sequence among parapoxviruses and other selected poxviruses was calculated using the pairwise distance method using MEGA-X software.
Genetic Diversity and the Host
Another phylogenetic tree was constructed using a gene at the variable region (ORF008) to determine the relationship of the ORFVs included in this study base on host. The Phylogenetic tree was constructed using the maximum likelihood methods and Tamura-Nei model (Tamura and Nei, 1993). The reliability of the tree was assessed by a bootstrap test with 1,000 replicates, constructed with MEGA-X (Kumar et al., 2015). The sequences of the ORF117 encoding granulocyte-macrophage colony-stimulating factor (GM-CSF) inhibitory factor were analyzed to study virus genetic diversity in relation to host using Bioedit bioinformatic software.
Results
Virus Isolation and Identification
The clinical sample collected from a goat suspected of contagious ecthyma infection was processed using standard procedure. A confluent monolayer of LT cells infected with the processed sample displays typical characteristics of ORFV CPE. These morphological changes include ballooning, rounding, and degeneration of cells which begins to appear 3 days post-inoculation. The CPE becomes more pronounced in 6–7 days with an increase in cytoplasmic granulation and 70–80% cell detachment. Purification of the concentrated virus using sodium diatrizoate density gradient (Robinson et al., 1982; Guo et al., 2003) revealed a distinct white opalescent band typical of ORFV appeared between 40 and 30% at approximately 36%. Total DNA extracted from the scab samples, cell culture passages, and purified virions particles were subjected to molecular analysis. PCR amplification using a specific primer targeting the B2L gene produces positive bands at the expected region of approximately 1,137 bp from all the samples. Transmission electron microscopic examination of the purified virion reveals the presence of negatively stained virion particles (Figure 1).
Figure 1. Representative of clinical case of Orf virus (ORFV) infection, Cytopathic effect following ORFV infection of LT cells, and electron microscopy of the purified virus. (A) Goat showing multiple lesions. (B) Uninoculated LT cells. (C) Infected cells with cytopathic effect post-inoculation. (D) Electron micrograph of the purified virus.
Genome Sequence Analysis
Orf virus strain UPM/HSN-20 was successfully sequenced in this study using next-generation sequencing (NGS) on the Illumina platform. The raw reads were processed and mapped to the respective reference genome. De novo assembly of the quality sequence reads revealed a linear contiguous sequence of 132,124 bp with G + C content of 63.7% (Table 1). The assembly demonstrates high levels of identity 97.3–99.0% with the other reference ORFV strains deposited in the GenBank. Whereas, the percentage’s identity with other members of the parapoxviruses ranges from 78.2 to 93.8% when compared using NCBI BLAST. The open reading frames (ORFs) were predicted using RAST and BLASTP. A total of 127 predicted protein best hits with ORFV reference strain; the nucleotide sequence of ORF 002, 005, and 132 has low hit, while two genes 001 at the 5′ end and 134 at 3′ variable ends were confirmed to be deleted using comparative sequence analysis and absence of amplification. The order and numbering of the predicted genes including the newly recognized genes (12.5 and 107.5) were adopted from Mercer et al. (2006).
Table 1. List of reference parapoxviruses used in this study.
Like other poxviruses, the terminal regions of the ORFV genome are bounded by inverted terminal repeats (ITR). The analysis of ITR sequences of 17 ORFV strains revealed sequence variability (Figure 2). The left end (5′) sequence of SJ1, YX, GO, SY17, NA17, GZ18, CL18, and MP contained terminal BamHI site along with telomere resolution motifs (ATTTTTT-N(8)-TAAAT), as indicated with the black box (A), while NP, NZ2, NA1/11, OV-HN3/12, OV-IA82 and IHUMI-1, UPM/HSN-20, TVL, and OV-SA00 lack this motif. The right end (3′) sequence of SJ1, YX, GO, NP, SY17, NA17, GZ18, CL18, MP, and NZ2 contained terminal BamHI site along with telomere resolution motifs (ATTTTTT-N(8)-TAAAT), as indicated with the black box (B), while NA1/11, OV-HN3/12, OV-IA82 and IHUMI-1, UPM/HSN-20, TVL, and OV-SA00 lack this motif.
Figure 2. Analysis of ITRs of 17 Orfvimses (A) the left end (5') sequence alignment (B) the right end (3') sequence alignment. The tenninal BamHI site, telomere resolution motif are indicated with black box (ATTTTTT-N(8)-TAAAT).
Phylogenetic and Distance Analysis
A Phylogenetic tree was constructed for tracing the evolutionary relationship between parapoxviruses and closely related poxviruses using highly conserved chordopoxvirus conserved ORF111 gene. The most closely related strain with isolates in this present study based on the highly conserved gene is XY and NP ORFV strain from China, as shown in Figure 3. The phylogenetic tree also revealed that ORFVs are distantly related to Squirrelpox virus, while the Molluscum Contagiosum virus appears as an outgroup. In the case of the pairwise analysis based on ORF111, the result shows a lesser mutation rate among ORFVs, while the other parapoxviruses has shown greater genetic variation Table 2.
Figure 3. Phylogenetic tree based on highly conserved ORFIII gene of parapoxviruses and other closely related poxviruses. The tree was constructed using the maximum likelihood methods and Tamura-Nei model with 1000 bootstrap replicates in NIEGA-X software.
Table 2. Estimates of Evolutionary Divergence between Orf virus UPM/HSN-20 and other reference isolates.
Genetic Diversity and Host
A phylogenetic tree was constructed to investigate the relationship between viruses and the infected host based on the ORF008 gene. The ORFVs form two major clusters (Figure 4), with the lower major cluster, consisting of eight (8) isolates all originating from goats (SJ1, UPM/HSN-20, OV-SA00, MP, NA17, NP, XY, and GO), while the upper major cluster consists of nine (9) ORFVs, eight (8) isolates originating from sheep (OV-IA82, NZ2, TVL, OV-HN3/12, NA1/11, D1701, SY17, and CL18) and one (1) isolate originating from goat (GZ18). The GZ18 strain forms a sub-cluster with SY17 and CL18 with strains originating from sheep.
Figure 4. Phylogenetic tree based on ORF008 of parapoxviruses. The tree was constructed using the maximum likelihood methods and Tamura-Nei model with 1000 bootstrap replicates in NIEGA-X software.
The comparison of amino acids sequences of the ORF117 gene among 17 ORFV strains highlighted that the following strains such as GZ18, D1701, TVL, OV-IA82, SY17, CL18, OV-HN3/12, NA1/11, and NZ2 contained Arginine at position 97, Arginine at position 151, Asparagine at position 176, and Cysteine at position 178, while other strains consist Glycine, Cysteine, Serine, and Phenylalanine in that position (Figure 5). In addition, Alanine was found at position 148 in GZ18, D1701, TVL, OV-IA82, SY17, CL18, OV-HN3/12, NA1/11, and NZ2, while other strains contained Leucine except NP presented Proline.
Figure 5. Multiple alignment of deduced amino acid sequence (ORF117) of Orf virus UPM/HSN-20 and other reference isolated.
Discussion
Livestock farming plays an important role in the development of the agricultural economy and is essential to the livelihoods of millions of farmers in developing countries (Dror and Allen, 2011; Wodajo et al., 2020). It is a well-known fact that the progress of this agricultural sub-sector is mainly determined by the health status of the animal which, in turn, depends upon the intake of high-quality feed and prevention of occurrence of infectious diseases among herds. In Malaysia, ORFV has been recognized for many years as an important pathogen affecting the growth of the livestock industry (Al-Ajeeli, 1995; Bala et al., 2019, 2020). To limit the impact of the outbreaks and to take efficient steps toward control and eradication of the disease, studies on molecular heterogeneity of ORFVs circulation are of utmost importance.
The infected goat is presented with pyrexia and coalescing papillary, severe ulcerative, and scabby lesions of various sizes on the epidermis of the muzzle and lips (Figure 1). During the outbreak, two animals were severely affected, one died, while the other survived under specialized veterinary care. Other clinical symptoms observed include nasal discharge, painful cough, difficulty in breathing, and slow movement. Previous clinical case studies described that infected Boar goats developed multifocal, severe proliferative dermatitis accompanied by other clinical conditions including chronic pneumonia, arthritis, and lymphadenopathy (Nandi et al., 2011; Sadiq et al., 2017). The clinical signs and symptoms observed in this study are consistent with the characteristics of ORFV infection.
To support the clinical diagnosis, a confluent monolayer of LT cells was infected with the processed sample. Virus propagation and isolation in cell culture are regarded as a “gold standard” for the detection of poxviruses (Grossegesse et al., 2017). Although the method is technical, time-consuming, and does not provide rapid diagnosis of infectious diseases, it is the only reliable method for generating sufficient stock of viruses for further studies. Various primary cells and continuous cell lines including lamb testis cell, Madin-Darby bovine kidney (MDBK), and Madin-Darby ovine kidney (MDOK) have been reported to be suitable for the cultivation of ORFVs (Abdullah et al., 2015; Kumar et al., 2018; Wang G. et al., 2019). Viral replication is followed by morphological and biochemical changes which ultimately lead to cell destruction (Zamri-Saad et al., 1999; Bala et al., 2020). In this study, typical characteristics of ORFV growth in cell culture were observed in LT cells line similar to earlier reports (Hematian et al., 2016; Lawal et al., 2021). These morphological changes commonly induced by ORFV include rounding and enlargement of cells together with swelling known as ballooning which indicates the starting point of CPE and begins to appear as early as 3 days post-inoculation. The continuity of this process leads to advanced CPE where the accumulated swollen and enlarged cells clumps and detaches from the flask.
Cytopathic effect alone is not sufficient for identification of the presence of the virus because several factors can alter the results, an additional test is needed for visualization of virus particles. The introduction of negative staining and the wider availability of electron microscopes have made electron microscopy essential in the differential diagnostic cell culture and clinical samples (Hazelton and Gelderblom, 2003; Shen et al., 2013). Pattern recognition of the size and particle morphology leads to direct and rapid visualization of an infectious agent (Gelderblom and Männel, 2003; Dal Pozzo et al., 2011). Some recent studies have employed the use of transmission electron microscopy for in vivo and in vitro identification of ORFV from clinical suspension by negative staining (Nashiruddullah et al., 2018; Galante et al., 2019). Therefore, this study has utilized this technique to ascertain the presence of the virus following the purification procedure. Large virion particles with typical morphological characteristics of parapoxviruses were observed (Figure 1). The enveloped virions presented as an ovoid shape with a unique spiral crisscross pattern that clearly distinguishes parapoxvirus genera from other poxvirus genera. For example, the TEM method has been useful for differential diagnosis between ORFV and capripoxviruses, all of which cause disease with similar clinical manifestation in sheep and goats (Rajesh and Swathi, 2018).
Next, to elucidate the genomic features of ORFV strain UPM/HSN-20, a high-throughput DNA-sequencing approach was employed, a technique that proved to be efficient in revealing detailed genetic information of an organism. In addition to identifying pathogens faster, NGS technology and bioinformatics can provide new insights into the molecular mechanism of microbial pathogenesis (Balakrishnan et al., 2017; Gwinn et al., 2019). Compared to Sanger sequencing, NGS technology has enabled researchers to study and understand the microbial world from broader and deeper perspectives. Pathogen’s genomics and epidemiology enhance public health efforts in preventing transmission of infectious diseases (Mohd-Azmi et al., 2000; Cao et al., 2017). Sequence analysis of our isolate revealed a linear sequence of 132,124 bp (GenBank accession number MW537048) with G + C content is 63.7%. All these features are consistent with the reference strain. The genome appears to be typical of PPVs, containing a central conserved region and variable region with a genus-specific GC profile (Chi et al., 2015; Zhong et al., 2019). Comparative sequence analysis of ORFV strain UPM/HSN-20 with the reference data available reveals a high degree of identity 97.3–99.0% among ORFVs. In addition, our sequence analysis indicates that PCPVs are genetically more closely related to ORFV in comparison to BPSV. Even though, BPSV and PCPV share the same host, whereas ORFV natural hosts are sheep and goats (Hautaniemi et al., 2010). The NCBI blast analysis revealed that the average nucleotide sequence identity between ORFV and PCPV is approximately 93% in contrast to 85% identity between ORFV and BPSV.
We discovered the deletion of two genes in the genome’s terminal regions in the present study (001 at the left end and 134 at the right end). Chi et al. (2015) previously reported similar deletion in the ORFV strain’s genome terminal region (NP: KP010355). The core part of the genome is well conserved, the terminal regions show significant variability and surprising evidence of rearrangement of the terminal sequences. Monkeypox virus was the first poxvirus to show genome rearrangement between non-homologous sequences (Esposito et al., 1981). Our knowledge of genome rearrangement resulting in duplication, transposition, and deletion is still in its infancy, and relatively few studies have reported sequence deletion among less and highly adapted cell culture ORFV isolates (Fleming et al., 1995; Cottone et al., 1998; McInnes et al., 2001; Chi et al., 2015).
To determine the evolutionary pathways and connections between ORFVs and other groups of closely related poxviruses, a phylogenetic tree was constructed. A molecular phylogenetic tree is a functional tool scientists use to show the evolutionary history of a set of species or groups of organisms (McLennan, 2010). From the literature as well as looking at the genome map, we found that ORFV ORF111 gene, formerly known as A35 gene (Rehm and Roper, 2011) by the vaccinia virus designation, is highly conserved in all mammalian-tropic poxvirus and has no similarity to non-poxvirus proteins, giving no clues as to its function. However, the bioinformatics analysis we conducted reveals that the ORFV ORF111 protein is antigenic based on VaxiJen protein function prediction software. In addition, the wide conservation of the ORF111 gene among all poxviruses with a mammalian host range suggests that this protein has an important role. Since the ORF111 gene is highly conserved, we constructed a phylogenetic tree to understand the relationship between parapoxviruses and other closely related poxviruses having high GC content. The phylogenetic tree constructed (Figure 2) in this study showed ORFV strain formed two separate branches on the tree, with the isolate under study (MW537048) forming a distinct cluster with XY (KP010353) and NP (KP010355) with Chinese isolates, signifying a common ancestor. Indeed, the Malaysian isolate (MW537048) and the two Chinese isolates (KP010353 and KP010355) are so genetically related that they are believed to share a common ancestor during evolution. The phylogenetic tree also revealed that ORFVs are distantly related to Squirrelpox virus, while the Molluscum Contagiosum virus appears as an outgroup. This is consistent with previous reports on evolutionary analysis of ORFVs and other poxviruses (Andreani et al., 2019). Besides the phylogenetic tree, the pairwise distance is one of the most widely used methods for the determination of mutation rates among isolates and/or groups of organisms. As indicated in Table 2, a lesser mutation rate was observed among ORFVs. The isolate in this study shows a closer identity with Chinese isolates. In addition, generally, higher genetic variation is observed between Orf viruses and other parapoxviruses. Mutations are a mandatory event for survival as such the local environment and host are important determining factors for virus growth and survival.
To further understand ORFV genetic diversity based on host predilection, a phylogenetic tree was constructed using the ORF008 gene. The ORF008 gene was among the genes used in previous studies conducted by Chi et al. (2015), and Zhong et al. (2019), which show ORFV strain diversity in relation to originating host. The tree they constructed in their studies shows that Orf viruses that were isolated from Sheep cluster together and those isolated from goats cluster together (Chi et al., 2015; Velazquez-Salinas et al., 2018; Zhong et al., 2019), suggesting that different genes in the genome of ORFV may be associated with the host predilection. While in the present study, we found that one goat isolate forms a cluster with two sheep isolates (GZ18, SY17, and CL18). These three isolates share a significant sequence similarity of 99.5%. However, considering the limited ORFV genome sequenced data available worldwide in the database, more studies are needed to confirm the virus diversity in relation to the host. In addition, ORFV infects not only sheep and goats (Spyrou and Valiakos, 2015), and including other hosts will greatly enhance the understanding of virus genetic diversity in relation to the host.
The genes in the conserved region perform essential functions and some genes such as B2L and F1L encoded by ORF11 and ORF059 have been widely used for rapid clinical diagnosis of ORFV (Yogisharadhya et al., 2017, 2018). While the variable region encodes some proteins that interfere with the early pathogen recognition, antiviral response, host’s cytokines, and cellular signaling, leading to pro-inflammatory responses to promote viral survival and replication. An example of such protein is a granulocyte-macrophage colony-stimulating factor (GM-CSF) encoded by the ORF117 gene. The GM-CSF inhibitory factor inhibits the downstream antiviral signaling of granulocyte-macrophage colony-stimulating factor and interleukin-2 (Wang R. et al., 2019). Interestingly, the amino acid sequence alignment of this biological significance protein reveals some unique amino acid sequence variations among ORFV strains (Figure 5). The effect of this sequence variation among the ORFV strain is yet to be reported. It is a well-known fact that amino acid changes can lead to changes in protein structure or function, which can sometimes affect the pathogenic characteristics of the virus.
Conclusion
The present study has demonstrated the functional merit of both non-molecular and molecular approaches for the identification of ORFV obtained from a clinical specimen. Our study provided an overall genomic mapping of the recent Malaysian ORFV strain UPM-HSN-20. A comprehensive study has found that the putative biological characteristics of the Malaysian isolates were similar to those of other Orf viruses. Multiple sequence and phylogenetic analyses revealed that our isolate in this study is closely related to two ORFV isolates NP and XY from China. The analysis of ITR sequences of 17 ORFV strains revealed significant sequence variability. The study has provided useful information for further comprehensive understanding of ORFV molecular biology, epidemiology, and control of the disease.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found at: https://www.ncbi.nlm.nih.gov/, MW537048.
Author Contributions
M-LM-A, HM, JB, and KB participated in conceptualization of the idea and study design. HM, JB, and KB performed the experiments. HM, AB, ZL, and AG participated in the data analysis, interpretation, and preparation of manuscript draft. FJ and MN supervised the analysis and edited the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by grants from Malaysian Skim Geran Penyelidikan Fundamental (FRGS) FASA1/2019 Universiti Putra Malaysia: FRGS/1/2019/STG03/UPM/01/2.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
Abdullah, A. A., Ismail, M. F. B., Balakrishnan, K. N., Bala, J. A., Hani, H., Abba, Y., et al. (2015). Isolation and phylogenetic analysis of caprine Orf virus in Malaysia. VirusDisease 26, 255–259. doi: 10.1007/s13337-015-0278-4
Al-Ajeeli, K. S. A. (1995). Studies on Malaysian Isolates of Orf Virus. Doctoral dissertation. Universiti Pertanian Malaysia.
Andreani, J., Fongue, J., Khalil, J. Y. B., David, L., Mougari, S., Le Bideau, M., et al. (2019). Human infection with Orf virus and description of its whole genome, France, 2017. Emerg. Infect. Dis. 25, 2197–2204. doi: 10.3201/eid2512.181513
Azmi, M., and Field, H. J. (1993). Interactions between equine herpesvirus type 1 and equine herpesvirus type 4: T cell responses in a murine infection model. J. Gen. Virol. 74, 2339–2345. doi: 10.1099/0022-1317-74-11-2339
Bala, J. A., Balakrishnan, K. N., Abdullah, A. A., Adamu, L., Noorzahari, M. S., May, L. K., et al. (2019). An association of Orf virus infection among sheep and goats with herd health programme in Terengganu state, eastern region of the peninsular Malaysia. BMC Vet. Res. 15, 250–215. doi: 10.1186/s12917-019-1999-1
Bala, J. A., Balakrishnan, K. N., Abdullah, A. A., Kimmy, T., Abba, Y., Bin Mohamed, R., et al. (2018a). Dermatopathology of Orf virus (Malaysian isolates) in mice experimentally inoculated at different sites with and without dexamethasone administration. J. pathogens 2018:9207576. doi: 10.1155/2018/9207576
Bala, J. A., Balakrishnan, K. N., Abdullah, A. A., Mohamed, R., Jesse, F. F. A., Noordin, M. M., et al. (2018b). The re-emerging of orf virus infection: a call for surveillance, vaccination and effective control measures. Microb. Pathog. 120, 55–63. doi: 10.1016/j.micpath.2018.04.057
Bala, J. A., Balakrishnan, K. N., Jesse, F. F. A., Abdullah, A. A., bin Noorzahari, M. S., Ghazali, M. T., et al. (2020). Identification of strain diversity and phylogenetic analysis based on two major essential proteins of Orf viruses isolated from several clinical cases reported in Malaysia. Infect. Genet. Evol. 77:104076. doi: 10.1016/j.meegid.2019.104076
Balakrishnan, K. N., Abdullah, A. A., Bala, J., Abba, Y., Sarah, S. A., Jesse, F. F. A., et al. (2017). Identification and comparison of RCMV ALL 03 open reading frame (ORF) among several different strains of cytomegalovirus worldwide. Infect. Genet. Evol. 54, 81–90. doi: 10.1016/j.meegid.2017.06.020
Cao, Y., Fanning, S., Proos, S., Jordan, K., and Srikumar, S. (2017). A review on the applications of next generation sequencing technologies as applied to food-related microbiome studies. Front. Microbiol. 8:1829. doi: 10.3389/fmicb.2017.01829
Chen, H., Li, W., Kuang, Z., Chen, D., Liao, X., Li, M., et al. (2017). The whole genomic analysis of Orf virus strain HN3/12 isolated from Henan province, Central China. BMC Vet. Res. 13, 1–10. doi: 10.1186/s12917-017-1178-1
Chi, X., Zeng, X., Hao, W., Li, M., Li, W., Huang, X., et al. (2013). Heterogeneity among Orf virus isolates from goats in Fujian Province, Southern China. PLoS One 8:e66958. doi: 10.1371/journal.pone.0066958
Chi, X., Zeng, X., Li, W., Hao, W., Li, M., Huang, X., et al. (2015). Genome analysis of Orf virus isolates from goats in the Fujian Province of southern China. Front. Microbiol. 6:1135. doi: 10.3389/fmicb.2015.01135
Cottone, R., Büttner, M., Bauer, B., Henkel, M., Hettich, E., and Rziha, H. J. (1998). Analysis of genomic rearrangement and subsequent gene deletion of the attenuated Orf virus strain D1701. Virus Res. 56, 53–67. doi: 10.1016/S0168-1702(98)00056-2
Cubells, J. R. E., Braverman, I., Kashgarian, M., and Lazova, R. (2016). A 65-year-old female from Connecticut with orf infection. Dermatopathology 3, 55–60. doi: 10.1159/000447125
Dal Pozzo, F., Martinelle, L., Gallina, L., Mast, J., Sarradin, P., Thiry, E., et al. (2011). Original findings associated with two cases of bovine papular stomatitis. J. Clin. Microbiol. 49, 4397–4400. doi: 10.1128/JCM.05281-11
Delhon, G. U. S. T. A. V. O., Tulman, E. R., Afonso, C. L., Lu, Z., De La Concha-Bermejillo, A., Lehmkuhl, H. D., et al. (2004). Genomes of the parapoxviruses Orf virus and bovine papular stomatitis virus. J. Virol. 78, 168–177. doi: 10.1128/JVI.78.1.168-177.2004
Dror, D. K., and Allen, L. H. (2011). The importance of milk and other animal-source foods for children in low-income countries. Food Nutr. Bull. 32, 227–243. doi: 10.1177/156482651103200307
Esposito, J. J., Cabradilla, C. D., Nakano, J. H., and Obijeski, J. F. (1981). Intragenornic sequence transposition in monkeypox virus. Virology 109, 231–243. doi: 10.1016/0042-6822(81)90495-5
Fleming, S. B., Lyttle, D. J., Sullivan, J. T., Mercer, A. A., and Robinson, A. J. (1995). Genomic analysis of a transposition-deletion variant of Orf virus reveals a 3.3 kbp region of non-essential DNA. J. Gen. Virol. 76, 2969–2978. doi: 10.1099/0022-1317-76-12-2969
Fleming, S. B., Wise, L. M., and Mercer, A. A. (2015). Molecular genetic analysis of Orf virus: a poxvirus that has adapted to skin. Viruses 7, 1505–1539. doi: 10.3390/v7031505
Friederichs, S., Krebs, S., Blum, H., Lang, H., and Büttner, M. (2015). Parapoxvirus (PPV) of red deer reveals subclinical infection and confirms a unique species. J. Gen. Virol. 96, 1446–1462. doi: 10.1099/vir.0.000080
Galante, D., Cafiero, M. A., Raele, D. A., Pugliese, N., Padalino, I., Cavaliere, N., et al. (2019). Identification and characterization of Orf viruses isolated from sheep and goats in southern Italy. Vet. Ital. 55, 347–353. doi: 10.12834/VetIt.1025.5477.2
Gelderblom, H., and Männel, A. (2003). Diagnostic electron microscopy in infectious diseases: update 2003. Microsc. Microanal. 9, 432–433. doi: 10.1017/S1431927603035049
Grossegesse, M., Doellinger, J., Tyshaieva, A., Schaade, L., and Nitsche, A. (2017). Combined proteomics/genomics approach reveals proteomic changes of mature Virions as a novel poxvirus adaptation mechanism. Viruses 9:337. doi: 10.3390/v9110337
Günther, T., Haas, L., Alawi, M., Wohlsein, P., Marks, J., Grundhoff, A., et al. (2017). Recovery of the first full-length genome sequence of a parapoxvirus directly from a clinical sample. Sci. Rep. 7, 1–8. doi: 10.1038/s41598-017-03997-y
Guo, J., Zhang, Z., Edwards, J. F., Ermel, R. W., Taylor, C. Jr., and De La Concha-Bermejillo, A. (2003). Characterization of a north American Orf virus isolated from a goat with persistent, proliferative dermatitis. Virus Res. 93, 169–179. doi: 10.1016/S0168-1702(03)00095-9
Gwinn, M., MacCannell, D., and Armstrong, G. L. (2019). Next generation sequencing of infectious pathogens. JAMA 321, 893–894. doi: 10.1001/jama.2018.21669
Hasheminasab, S. S., Mahmoodi, A., Mahmoodi, P., and Maghsood, H. (2016). Orf virus infection in human ecthyma contagiosum: a report of two cases in the west of Iran. Virus 27, 209–210. doi: 10.1007/s13337-016-0304-1
Hatcher, E. L., Wang, C., and Lefkowitz, E. J. (2015). Genome variability and gene content in chordopoxviruses: dependence on microsatellites. Viruses 7, 2126–2146. doi: 10.3390/v7042126
Hautaniemi, M., Ueda, N., Tuimala, J., Mercer, A. A., Lahdenperä, J., and McInnes, C. J. (2010). The genome of pseudocowpoxvirus: comparison of a reindeer isolates and a reference strain. J. Gen. Virol. 91, 1560–1576. doi: 10.1099/vir.0.018374-0
Hazelton, P. R., and Gelderblom, H. R. (2003). Electron microscopy for rapid diagnosis of emerging infectious agents. Emerg. Infect. Dis. 9, 294–303. doi: 10.3201/eid0903.020327
Heare, D. L., Little, S. V., Weise, D. W., Harris, J. R., Hillhouse, A. E., Konganti, K., et al. (2020). Whole-genome sequence of an Orf virus isolate derived from a cell culture infected with contagious Ecthyma vaccine. Microbiol. Res. Ann. 9:e00720. doi: 10.1128/MRA.00752-20
Hematian, A., Sadeghifard, N., Mohebi, R., Taherikalani, M., Nasrolahi, A., Amraei, M., et al. (2016). Traditional and modern cell culture in virus diagnosis. Osong. Public Health Res. Perspect. 7, 77–82. doi: 10.1016/j.phrp.2015.11.011
Huang, T., Tulman, E. R., Diel, D. G., Khatiwada, S., Sims, W., Edwards, J. F., et al. (2015). Coinfection with multiple strains of bovine papular stomatitis virus. Arch. Virol. 160, 1527–1532. doi: 10.1007/s00705-015-2394-2
Jesse, F. F. A., Latif, S. N. A. A., Abba, Y., Hambali, I. U., Bitrus, A. A., Peter, I. D., et al. (2018). Seroprevalence of orf infection based on IgM antibody detection in sheep and goats from selected small ruminant farms in Malaysia. Comp. Clin. Pathol. 27, 499–503. doi: 10.1007/s00580-017-2619-8
Kassa, T. (2021). A review on human Orf: a neglected viral zoonosis. Res. Rep. Tropic. Med. 12, 153–172. doi: 10.2147/RRTM.S306446
Kumar, S., Stecher, G., Li, M., Knyaz, C., and Tamura, K. (2018). MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 35, 1547–1549. doi: 10.1093/molbev/msy096
Kumar, R., Trivedi, R. N., Bhatt, P., Khan, S. H., Khurana, S. K., Tiwari, R., et al. (2015). Contagious pustular dermatitis (orf disease)–epidemiology, diagnosis, control and public health concerns. Adv. Anim. Vet. Sci. 3, 649–676. doi: 10.14737/journal.aavs/2015/3.12.649.676
Lawal, N. U., Ibrahim, M., Onawala, D. A., Bello, M. B., Aliyu, R. M., Baraya, Y. U. S. A., et al. (2021). Molecular characterization and phylogenetic analysis of Orf virus isolated from goats in Sokoto metropolis, Nigeria. Future Sci. OA 7:FSO700. doi: 10.2144/fsoa-2020-0162
Lederman, E. R., Green, G. M., DeGroot, H. E., Dahl, P., Goldman, E., Greer, P. W., et al. (2007). Progressive Orf virus infection in a patient with lymphoma: successful treatment using imiquimod. Clin. Infect. Dis. 44, e100–e103. doi: 10.1086/517509
Li, W., Hao, W., Peng, Y., Duan, C., Tong, C., Song, D., et al. (2015). Comparative genomic sequence analysis of Chinese Orf virus strain NA1/11 with other parapoxviruses. Arch. Virol. 160, 253–266. doi: 10.1007/s00705-014-2274-1
McGuire, M. J., Johnston, S. A., and Sykes, K. F. (2012). Novel immune-modulator identified by a rapid, functional screen of the parapoxvirus ovis (Orf virus) genome. Proteome Sci. 10:4. doi: 10.1186/1477-5956-10-4
McInnes, C. J., Wood, A. R., Nettleton, P. F., and Gilray, J. A. (2001). 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 22, 141–150. doi: 10.1023/A:1008117127729
McInnes, C. J., Wood, A. R., Thomas, K., Sainsbury, A. W., Gurnell, J., Dein, F. J., et al. (2006). Genomic characterization of a novel poxvirus contributing to the decline of the red squirrel (Sciurus vulgaris) in the UK. J. Gen. Virol. 87, 2115–2125. doi: 10.1099/vir.0.81966-0
McLennan, D. A. (2010). How to read a phylogenetic tree. Evol. Educ. Out. 3, 506–519. doi: 10.1007/s12052-010-0273-6
Mercer, A. A., Ueda, N., Friederichs, S. M., Hofmann, K., Fraser, K. M., Bateman, T., et al. (2006). Comparative analysis of genome sequences of three isolates of Orf virus reveals unexpected sequence variation. Virus Res. 116, 146–158. doi: 10.1016/j.virusres.2005.09.011
Mohd-Azmi, M. L., Ali, A. S., and Kheng, W. K. (2000). DNA fingerprinting of red jungle fowl, village chicken and broilers. Asian-Australasian Journal of Animal Sciences 13, 1040–1043. doi: 10.5713/ajas.2000.1040
Nandi, S., De, U. K., and Chowdhury, S. (2011). Current status of contagious ecthyma or orf disease in goat and sheep—a global perspective. Small Rumin. Res. 96, 73–82. doi: 10.1016/j.smallrumres.2010.11.018
Nashiruddullah, N., Pathak, D. C., Barman, N. N., Ahmed, J. A., Borah, P., Begum, S. S., et al. (2018). In vitro and in vivo assessment of Orf virus (ORFV) by electron microscopy. Veterinarski arhiv 88, 847–861. doi: 10.24099/vet.arhiv.0229
Peralta, A., Robles, C. A., Micheluod, J. F., Rossanigo, C. E., Martinez, A., Carosio, A., et al. (2018). Phylogenetic analysis of Orf viruses from five contagious ecthyma outbreaks in Argentinian goats. Front. Vet. Sci. 5:134. doi: 10.3389/fvets.2018.00134
Rajesh, K., and Swathi, B. (2018). Diagnosis of goat pox with electron microscopy and its management. Pharma Inno. J. 7, 567–569.
Rehm, K. E., and Roper, R. L. (2011). Deletion of the A35 gene from modified Vaccinia virus Ankara increases immunogenicity and isotype switching. Vaccine 29, 3276–3283. doi: 10.1016/j.vaccine.2011.02.023
Robinson, A. J., Ellis, G., and Balassu, T. (1982). The genome of Orf virus: restriction endonuclease analysis of viral DNA isolated from lesions of orf in sheep. Arch. Virol. 71, 43–55. doi: 10.1007/BF01315174
Sadiq, M. A. A., Abba, Y., Jesse, F. F. A., Chung, E. L. T., Bitrus, A. A., Abdullah, A. A., et al. (2017). Severe persistent case of con-tagious ecthyma (Orf) in goats. J. Anim. Health Prod 5, 24–28. doi: 10.14737/journal.jahp/2017/5.1.24.28
Sahu, B., Majee, P., Singh, R., and Nayak, D. (2020a). Recombination may drive the emergence of Orf virus diversity: evidence from the first complete genome of Indian Orf virus and comparative genomic analysis. Arch. Virol. 167, 1571–1576. doi: 10.1007/s00705-022-05443-5
Sahu, B. P., Majee, P., Singh, R. R., Sahoo, A., and Nayak, D. (2020b). Comparative analysis, distribution, and characterization of microsatellites in Orf virus genome. Sci. Rep. 10:13852. doi: 10.1038/s41598-020-70634-6
Shen Ni, L., Allaudin, Z. N. b., Mohd Lila, M. A. b., Othman, A. M. b., and Othman, F. b. (2013). Selective apoptosis induction in MCF-7 cell line by truncated minimal functional region of Apoptin. BMC Cancer 13:488. doi: 10.1186/1471-2407-13-488
Spyrou, V., and Valiakos, G. (2015). Orf virus infection in sheep or goats. Vet. Microbiol. 181, 178–182. doi: 10.1016/j.vetmic.2015.08.010
Tamura, K., and Nei, M. (1993). Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol. Biol. Evol. 10, 512–526. doi: 10.1093/oxfordjournals.molbev.a040023
Velazquez-Salinas, L., Ramirez-Medina, E., Bracht, A. J., Hole, K., Brito, B. P., Gladue, D. P., et al. (2018). Phylodynamics of parapoxvirus genus in Mexico (2007–2011). Infect. Genet. Evol. 65, 12–14. doi: 10.1016/j.meegid.2018.07.005
Wang, G., Wang, Y., Kong, J., Li, Y., Wu, J., Chen, Y., et al. (2019). Comparison of the sensitivity of three cell cultures to ORFV. BMC Vet. Res. 15, 1–8. doi: 10.1186/s12917-018-1760-1
Wang, R., Wang, Y., Liu, F., and Luo, S. (2019). Orf virus: A promising new therapeutic agent. Rev. Med. Virol. 29:e2013. doi: 10.1002/rmv.2013
Wodajo, H. D., Gemeda, B. A., Kinati, W., Mulem, A. A., van Eerdewijk, A., and Wieland, B. (2020). Contribution of small ruminants to food security for Ethiopian smallholder farmers. Small Rumin. Res. 184:106064. doi: 10.1016/j.smallrumres.2020.106064
Yang, H., Meng, Q., Qiao, J., Peng, Y., Xie, K., Liu, Y., et al. (2014). Detection of genetic variations in Orf virus isolates epidemic in Xinjiang China. J. Basic Microbiol. 54, 1273–1278. doi: 10.1002/jobm.201300911
Yogisharadhya, R., Kumar, A., Bhanuprakash, V., and Shivachandra, S. B. (2018). Evaluation of a recombinant major envelope protein (F1L) based indirect-ELISA for sero-diagnosis of orf in sheep and goats. J. Virol. Methods 261, 112–120. doi: 10.1016/j.jviromet.2018.08.015
Yogisharadhya, R., Kumar, A., Ramappa, R., Venkatesan, G., Bhanuprakash, V., and Shivachandra, S. B. (2017). Functional characterization of recombinant major envelope protein (rB2L) of Orf virus. Arch. Virol. 162, 953–962. doi: 10.1007/s00705-016-3178-z
Zamri-Saad, M., Effendy, A. W. M., Israf, D. A., and Azmi, M. L. (1999). Cellular and humoral responses in the respiratory tract of goats following intranasal stimulation using formalin-killed Pasteurella haemolytica A2. Vet. Microbiol. 65, 233–240. doi: 10.1016/s0378-1135(98)00298-3
Zhong, J., Guan, J., Zhou, Y., Cui, S., Wang, Z., Zhou, S., et al. (2019). Genomic characterization of two Orf virus isolates from Jilin province in China. Virus Genes 55, 490–501. doi: 10.1007/s11262-019-01666-y
Keywords: ruminants, cell culture, genome, next-generation sequencing, Orf virus
Citation: Mangga HK, Bala JA, Balakrishnan KN, Bukar AM, Lawan Z, Gambo A, Jesse FFA, Noordin MM and Mohd-Azmi M-L (2022) Genome-Wide Analysis and Molecular Characterization of Orf Virus Strain UPM/HSN-20 Isolated From Goat in Malaysia. Front. Microbiol. 13:877149. doi: 10.3389/fmicb.2022.877149
Edited by:
Erna Geessien Kroon, Federal University of Minas Gerais, BrazilReviewed by:
Andrea Peralta, Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), ArgentinaLauro Velazquez-Salinas, United States Department of Agriculture, United States
Copyright © 2022 Mangga, Bala, Balakrishnan, Bukar, Lawan, Gambo, Jesse, Noordin and Mohd-Azmi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Hassana Kyari Mangga, hkmangga@unimaid.edu.ng; Mohd-Lila Mohd-Azmi, azmi@upm.edu.my