Genome annotation is the analysis of genome sequence of a particular species which includes all the possible analysis of DNA sequence that can be done by computational means. Raw DNA sequence produced by the genome-sequencing projects is taken for the analysis. Analysis and interpretation is carried out which is necessary to extract its biological significance and organize the information into the context of our understanding of biological processes. By genome annotation, an individual gene and its protein (or RNA) product can be predicted by in depth analysis of several gene features, such as ORFs length, % (G + C) content, promoter, polyA tail, homology, etc. The focal point of each such record is the validity of a ORF as a potential functional gene. Moreover, analysis may also include a brief description of the evidence for the assigned or proposed function. Several genome annotation tools are freely available for further analysis of genome sequence. Users can use several individual tools for each task or can use some integrated tools that do many task simultaneously. GLIMMER [1] , GeneMark [2] , ORF finder [3] and FrameD [4] are well-used tools for predicting Open Reading Frame (ORF). BLAST is the most used homology tool for homology prediction but there are many other options are available, such as MPsrch, PSI-BLAST and WU-BLAST [5] . Artemis, Apollo, JBROWSE, etc. are used for Genome browsers for curation and further annotation [6] . There are also many Genome annotation pipelines available, such as PASA and MAKER, in which several analysis tools are integrated [6] . Users are required to select the tools or pipeline according to their needs that serve best for their purpose

Contagious pustular dermatitis; ORF is a common epitheliotropic viral disease of sheep, goats and wild ruminants and is characterized by the formation of papules, nodules, or vesicles caused by the ORF virus [7] . Humans always posses a high risk due to zoonotic characteristics of this disease and often humans can contract this disorder through direct contact with infected animals by the fomites that carry the ORF virus [8] . Purulent-appearing papule is the major symbol causes locally and generally no systemic symptoms is obseved [3] . Normally it infects the finger, hand, arm, face and even the penis [6] - [9] . ORF virus (ORFV) is an oval and enveloped virus containing dsDNA genome within the genus Parapoxvirus, family Poxviridae [6] . The genus also includes pseudocowpox virus (PCPV) and bovine papular stomatitis virus (BPSV) in cattle and parapoxvirus (PPV) of red deer in New Zealand [10] .

Mechanisms involved in ORFV virulence are not well studied [10] . Several putative virulence genes have been identified, such as the vascular endothelial growth factor homologue, an interleukin (IL-10) homologue, a double-stranded RNA-binding protein, a factor inhibiting the cytokines granulocyte-macrophage colony-stimu- lating factor and IL-2 [11] - [16] . The whole genome of ORF virus strain OV-IA82 has been sequenced by the Plum Island Animal Disease Center, USA in 2004 [17] . The genome has high 63% GC content with 89% coding region which codes 130 proteins and length is 139,962 bp [17] . The complete genomic sequence available for OV-IA82 enables to deduce the different complex area including other molecular aspects of pathogenesis of this virus in details [18] . In this study, we tried to identify unique genes in ORFV genome and characterize them by Insilco process.

There are certain kinds of genome annotation tools available for analysis the whole genome, some of which Artemis and Artemis comparison Tools (ACT) has been used extensively for Open Reading Frame (ORF) visualization, editing, determine GC plot as well as retrieve desired nucleotide length within the genome for further analysis [19] . The GLIMMER has been chosen for predicting the genes or ORF in the genome of ORFV based up on its accuracy and flexibility of changing stop and start codons according to the requirement [20] [21] . Potential protein coding ORFs were identified by the following criteria: ORF size larger than 60 amino acid (aa), presence of potential transcriptional start and stop sites, a high GLIMMER score and homology to other known Parapoxvirus or cellular ORFs [21] . To find the similarity and homology, sequences of the GLIMMER pedicted protein, coding ORF were compared with the protein databases such as SwissProt, Trembl, UniProt, etc. In this study, among the tools which are available to find the similarity compared to protein databases, the Blast module Blastall, which supports all five Blast programs (blastp, blastn, blastx, tblastn and tblastx) has been chosen for finding the unique genes. Promoter, poly A signal and CpG island were analyzed for each unique genes in order to determine their potentiality as genes. For promoter prediction upstream ~350 bases were subjected within Neural Network Promoter Prediction Tool developed by University of Berkly, USA with a cut off value 0.8 [22] . PolyA signals were determined by polyADQ which is a poly (A) signal search engine developed by Cold Spring Harbor Laboratory, USA [22] [23] . The software CpGIE was downloaded via the website: [24] . The following cutoff values were used to determine the CpG island in a given genomic sequence: ≥200 nt, G + C content 50%, and an observed: expected CpG ratio 0.6 [25] [26] . To determine the Trans Membrane (TM) domain and peptide signal of the unique genes TMHMM and SignalP 3.0 software was used respectively [27] . The Artemis and Artemis Comparison Tool (ACT), written by Kim Rutherford, Sanger Institute, UK [19] were used for genome analysis and the pair wise comparisons between OV- IA82, VACV and BPSV genomes. For phylogenetics analysis various B2L genes of ORFV were selected on the basis of host specificity such as human, sheep and goats [28] [29] . In order to determine distant relationship BPSV and psedocowpox B2L gene also selected. All the genes sequences were downloaded from gene bank and stored for phylogenetic tree construction as listed in Table 1. The nucleotide sequences of diverse ORf viruses and others were aligned using the Bioedit program and Mega 5 software [22] . One thousand bootstrap replicates were subjected to nucleotide sequence distance and neighbor-joining methods, and the consensus phylogenetic tree was drawn.

In this study we tried to identify unique genes in ORFV on the basis of the availability of complete genome sequence of OV-IA82 by Insilco process in order to determine the novel vaccine strain. Recent advances in bioinformatics enable us to further analysis of the ORFV genome and the results are given below.

ORF virus shares specific genomic features with other poxviruses as VACV and BPSV, in terms of genome organization and gene content. Comparative genome sequences with other two closely related viruses VACV and BPSV here provide a comparative view of Parapoxvirus (PPV) genomics and basic knowledge of viral functions associated with virus replication and manipulation of cellular responses. Based on comparative genomic analysis, the genomes of BPSV and ORFV differ significantly which may be responsible for differences in host range. Modern genome analysis tools, such as Artemis Comparison Tools (ACT), promptus to understand the PPV biology which will permit the engineering of novel vaccine viruses and expression vectors with enhanced efficacy and greater versatility [30] . Nevertheless, we have identified four unique genes with potential importance as virulence factors and thus they could be vaccine candidates in the future. These genes should not correspond to horizontal gene transfer and their characteristic features may be consequences of the specific evolutionary cycle that shapes the ORFV gene repertoires in the context of their parasitic lifestyles [31] . ORF 039 with signal peptides and transmembrane domains may be directly toxic or confer association with the host. Therefore, further focus can be placed on subset of this for functional analysis. The function, subcellular location, average of

hydrophobicity and protein regions that share a significant degree of sequence similarity with known protein family can be detected by using computational approach. Primarily 2D-PAGE might be used for membrane protein analysis for ORF 039 [32] [33] . Besides, Combination of nano liquid chromatography (NanoLC) and mass spectrometry (MS) can be used to detect the transmembrane domian of ORF 039 [33] . ORF 039 can be used for Insilco protein homology modeling. Identification of certain protein ligand of this protein may open a new door for the development of antibody for ORFV and also for potential drug target. Insilco drug target analysis may be carried out to find out the potentiality of the genes for further analysis.

Similar Insilco studies by the following different approach have been carried out to identify potential genes for therapeutic targets [34] [35] . Unique genes or proteins that are involved in a certain pathway or having certain characteristics like outer membrane protein, presence of unique protein ligand family, etc. are always an ideal candidate for drug target [35] - [37] . Insilco identification of target genes and prediction of drug candidate is a well-established methodology in drug discovery. Screening out the target genes and their corresponding drug by Insilco approach makes the drug discovery procedure more robust, quick and economically feasible [36] .

The development of novel vaccine strain will control contagious pustular dermatitis in a small ruminant, potential source of leather, meat and milk which have a significant role in the economy of a country [13] . Therefore, control of such ORF related disease will have a significant role in economic development of a country. This study will be very much useful for further study of the evolution of the ORFV that may provide an encouragement for the development of new diagnostic tools and medicines.

The authors wish to thank Ministry of Education (17/10, M-15/2007/226 and HEQEP, CP-W1-3413), Bangladesh for partially funding this research work.