Skip to main content
Download PDF

Molecular characterization of recombinant LSDV isolates from 2022 outbreak in Indonesia through phylogenetic networks and whole-genome SNP-based analysis

BMC Genomics volume 25, Article number: 240 (2024) Cite this article

Abstract

Lumpy skin disease (LSD) is a transboundary viral disease of cattle and water buffaloes caused by the LSD virus, leading to high morbidity, low mortality, and a significant economic impact. Initially endemic to Africa only, LSD has spread to the Middle East, Europe, and Asia in the past decade. The most effective control strategy for LSD is the vaccination of cattle with live-attenuated LSDV vaccines. Consequently, the emergence of two groups of LSDV strains in Asian countries, one closely related to the ancient Kenyan LSDV isolates and the second made of recombinant viruses with a backbone of Neethling-vaccine and field isolates, emphasized the need for constant molecular surveillance. This current study investigated the first outbreak of LSD in Indonesia in 2022. Molecular characterization of the isolate circulating in the country based on selected LSDV-marker genes: RPO30, GPCR, EEV glycoprotein gene, and B22R, as well as whole genome analysis using several analytical tools, indicated the Indonesia LSDV isolate as a recombinant of LSDV_Neethling_vaccine_LW_1959 and LSDV_NI-2490. The analysis clustered the Indonesia_LSDV with the previously reported LSDV recombinants circulating in East and Southeast Asia, but different from the recombinant viruses in Russia and the field isolates in South-Asian countries. Additionally, this study has demonstrated alternative accurate ways of LSDV whole genome analysis and clustering of isolates, including the recombinants, instead of whole-genome phylogenetic tree analysis. These data will strengthen our understanding of the pathogens’ origin, the extent of their spread, and determination of suitable control measures required.

Peer Review reports

Introduction

Lumpy skin disease virus (LSDV) is a member of the Capripoxvirus (CaPV) genus that has a linear double-stranded DNA genome [1]. The virus causes lumpy skin disease (LSD) in cattle and water buffaloes [2]. The CaPV genus comprises two other species, the sheeppox virus (SPPV) and goatpox virus (GTPV), which are antigenically and genetically closely related to LSDV [1, 3]. LSD infections are associated with high morbidity (up to 100%) and low mortality (< 10%), with a significant economic impact due to losses related to animal production, restrictions to trade, and treatment/vaccination costs, and therefore a significant World Organisation for Animal Health (WOAH) notifiable disease [2, 4]. Clinically sick animals present symptoms such as fever, enlarged lymph nodes, nasal discharge, and nodular skin lesions [5, 6].

LSDV is transmitted by blood-feeding vectors such as stable flies, mosquitoes, and ticks [7,8,9]. This may have facilitated its spread from the endemic region in Africa to the Middle East, Europe, and South Asia in the past ten years [10]. Other modes of transmission, such as direct contact, seminal, and intrauterine transmission, have also been demonstrated experimentally [11,12,13]. LSD control measures include cattle movement/trade restriction, vector control, disinfection, or culling of the animals showing clinical signs [2, 14], with the most effective control strategy being vaccination using live attenuated vaccines, such as those derived from the Neethling strain [15,16,17].

Understanding the molecular and genetic epidemiology of CaPVs is a prerequisite for better management and control of LSD, given the antigenic similarity of CaPVs that challenges their distinction by serological tests. Among several molecular-based techniques for genetically distinguishing CaPVs include, the sequence analysis of genes encoding the RNA polymerase 30 kDa subunit (RPO30), the G-protein-coupled receptor (GPCR), the CaPV homolog of the variola virus B22R, and the EEV glycoprotein that possess specific signatures that can differentiate the CaPVs [18,19,20,21]. These genes could also help distinguish live attenuated vaccines from field isolates and the recombinant LSDV strains. However, a whole genome-based approach is preferable as it enables proper recombination analysis which has been done for the recently reported recombinants [22, 23]. Nevertheless, phylogenetic tree construction using recombinants’ whole genome is not suitable/accurate since the analysis process considers mainly the mutations and speciation events and does not account for recombination between the sequences or ancestor sequences [24]. Therefore, alternative tools such as restricting the analysis to parts of the virus genome, or phylogenetic networks and single nucleotide polymorphism (SNP) genotyping, which consider events like recombination/hybridization have been shown as more appropriate for whole genome analysis of recombinants [25,26,27].

LSD was first reported in Asia in 2019 in India, Bangladesh, and China and subsequently spread to other Asian countries such as Nepal, Myanmar, Thailand, and Vietnam [28]. Phylogenetic analysis of the Asian LSDV isolates has demonstrated that the isolate circulating in South-Asian countries resembles the ancient LSDV isolate from Kenya [29,30,31], while the isolates in Southeast and Eastern Asia are recombinant viruses with Neethling vaccine strain and field isolates backbones [23, 32,33,34,35]. In Indonesia, the initial suspected clinical cases of LSD were reported in November 2020 cattle and buffaloes from Kampar and buffaloes from Kuansing districts, Riau Province. In 2021, more suspicions were reported in cattle from Padang Pariaman district in North Sumatera (January) and in a buffalo from Muara Bungo district (August). However, none of the suspected LSD cases from 2020 to 2021 could be confirmed by either serological (Capripox ELISA) or molecular tests (Polymerase chain reaction (PCR). In February 2022, the disease spread to cattle from Indragiri Hulu district, where samples were collected and confirmed by PCR at the Research Centre for Veterinary Science (RCVS), National Research and Innovation Agency, enabling the notification to WOAH on February 2022 [36]. This study aims to provide clinical evidence of the first outbreak of LSD in Indonesia in 2022 and the molecular characterization of the LSDV isolate circulating in the country based on selected LSDV-marker genes and whole genome analysis. Due to the recent emergence of LSDV recombinants that make it challenging (less accurate) to cluster LSDV isolates using whole genome phylogenetic trees, the study will also explore some alternative approaches based on existing tools for better accuracy in LSDV whole genome-based clustering.

Results

LSD outbreak investigation and clinical signs

The suspected LSD infected animals presented nodular lesions on different body parts, nasal discharge, and fever (Fig. 1). The morbidity and mortality rates recorded in the Indragiri Hulu district by February 2022 were 28.04% (30/107) and 1.87% (2/107), respectively.

Fig. 1
figure 1

Clinical signs of LSD suspected cases in Indonesia; nodules on cattle from Indragiri Hulu district

Full size image

Molecular detection and diagnosis of LSDV

CaPV genome was detected in ten (83.33%) out of the twelve samples collected from 2022 suspected cattle cases in Indragiri Hulu district, Indonesia, by RT-PCR using the Bowden et al. protocol [6]. Further analysis of the positive samples, using HRM assay, confirmed that the suspected cases were infected with LSDV (Additional file 1: Table s1).

Targeted gene sequence and phylogenetic analysis

For five samples, the targeted fragments were successfully amplified and sequenced, and the assembled contigs of 606, 1146, 327, and 832 bp for the RPO30, GPCR, EEV, and B22R genes, respectively, were deposited in the GenBank with accession numbers OR246967 to OR246986.

The multiple sequence alignments showed that all five Indonesia LSDV samples were 100% identical based on RPO30, GPCR, EEV, and B22R gene sequences. Phylogenetic analysis of the complete RPO30 gene sequences showed that Indonesia sequences were grouped in LSDV Cluster I together with LSDV Neethling-derived vaccines and other LSDV isolates from East and Southeast Asia, such as China, Taiwan, Vietnam, and Thailand (Fig. 2A). However, the GPCR gene sequences grouped the Indonesia sequences in LSDV Cluster II with the common LSDV field isolates, the historical LSDVs from Kenya, and LSDV KSGP-0240 derived vaccine, yet with the other East and Southeast Asia isolates (Fig. 2B).

Fig. 2
figure 2

Neighbour-joining tree based on the complete (A) RPO30 and (B) GPCR gene sequences of CaPVs with LSDVs from Indonesia (in red) visualized on iTOL. The evolutionary distances were computed using the Maximum Composite Likelihood method with 1000 bootstrap replicates on MEGA X

Full size image

The multiple sequence alignment of the partial EEV glycoprotein gene showed that Indonesian LSDV isolates are identical in this region of the genome to LSDV isolates from China, Vietnam, Taiwan, and LSDV Neethling-derived vaccines that have a 27-nucleotide deletion (175–201) (Fig. 3A). The B22R sequence alignment showed that part of the Indonesian sequence was identical to the historical LSDVs from Kenya (353 bp) and the other to the LSDV Neethling-derived vaccine (479 bp). However, the B22R of the Indonesian LSDVs differed from LSDV Neethling and LSDV KSGP-0240 derived vaccines due to one nucleotide InDel at positions 102 and 745, respectively (Fig. 3B).

Fig. 3
figure 3

Multiple sequence alignment of (A) the partial sequences of the EEV glycoprotein gene showing the 27-nucleotide deletion highlighted in the block found in Indonesia isolates and (B) partial sequences of the B22R gene showing the nucleotide insertion (in blocks) in LSDV_Neethling and LSDV_KSGPO-240 vaccines that are absent in Indonesia isolates (in red). The dots indicate the identical nucleotides in the alignment

Full size image

LSDV whole genome SNPs and recombination analysis

The assembled LSDV genomes were 151,010 bp long for each of the two Indonesian isolates; LSDV_Indonesia_2022_S1 and LSDV_Indonesia_2022_S4, with a mean coverage of 2,856 ± 814 and 637 ± 174 respectively. These genomes comprised 158 predicted ORFs each (Additional file 1: Table s2). Sixty-eight (68) of the predicted ORFs were 100% homologous to both LSDV_NI-2490 (LSDV Cluster II) and LSDV_Neethling_vaccine_LW_1959 (LSDV Cluster I) reference strains, while 26 ORFs were only 100% identical to LSDV_NI-2490 and 55 ORFs to LSDV_Neethling_vaccine_LW_1959. The remaining 9 ORFs were 99.5-99.9% homologous to LSDV_NI-2490 and LSDV_Neethling_vaccine_LW_1959 reference strains (Additional file 1: Table s3). NCBI nucleotide blast search based on the whole genome revealed that over the aligned genome fragment, the Indonesia LSDV genome shared 100% nucleotide identity to LSDVs from Thailand (OM033705.1), Vietnam (MZ577076.1) and China (OM803092.1), and 99.98% and 99.94% to LSDVs from Taiwan (OL752713.2) and HongKong (MW732649.1) respectively. In addition, the Indonesian LSDVs shared 99.56% and 99.25% nucleotide identity to LSDV_Neethling_vaccine_LW_1959 and LSDV_NI-2490, respectively. The shared identities between the Indonesian LSDV genomes and the recombinant LSDVs from Russia were 99.41% to Saratov_2017 (MH646674_1), 99.51% to Udmurtiya_2019 (MT134042_1), 99.54% to Saratov_2019 (OM530217_1) and 99.64% to Tyumen_2019 (OL542833_1).

SNP analysis of the aligned LSDV genomes relative to the LSDV_NI-2490 (AF325528.1) reference genome revealed 2145 SNPs distributed across the genome. The heatmap based on the SNP index showed seven clusters of LSDV genotypes; NI-like (related to the ancient NI-2490 LSDV), NW-like (commonly circulating LSDV field isolates), Neethling-like (related to Neethling_vaccine) and four recombinant clusters, recombinant_1-like (Saratov_Russia_isolates), recombinant_2-like (Udmurtiya_Russia_isolate), recombinant_3-like (East and Southeast Asia isolates) and recombinant_4-like (Tyumen_Russia_isolate) (Fig. 4A). The Indonesian LSDVs clustered with the recombinant_3 like LSDVs, together with LSDV isolates from Thailand, Vietnam, Taiwan, and China. The seven clusters of LSDVs based on the SNPs profile were also differentiable by the PCA scatterplot and the neighbor-joining tree (Fig. 4B). Further SNP analysis of the recombinant_3-like LSDVs relative to the LSDV_Indonesia_2022_S1, detected 133 SNPs across the 15 aligned genomes. While there was one SNP between the two Indonesian isolates, the heatmap revealed 2–5 SNPs between the Indonesia LSDVs and other East and Southeast Asian isolates such as Thailand_OM033705_1, China_OM373209_1/OM984485_1, Vietnam_MZ577075_1, and Taiwan_OL752713_2. The highest number of SNPs were detected in Isolates from China, MW355944_1_LSDV/GD01/2020 (11 SNPs) and OP508345_1_LSDV_China/Xinjiang/Cattle/2019 (23 SNPs) and HongKong MW732649_1_LSDV/HongKong/2020 (76 SNPs). Other notable LSDVs in recombinant_3-like cluster were OM793602_1_LSDV_Russia_Tomsk_2020 and OM793603_1_LSDV_Russia_Khabarovsk_2020 from Russia with 11 and 24 SNPs respectively (Fig. 4C).

Fig. 4
figure 4

Analysis of SNPs extracted from LSDV genome alignments. (A) SNP distribution heatmap, (B) PCA-scatter plot, and neighbor-joining tree based on the PCA scores of the SNP data of all LSDV complete genomes, (C) SNP distribution heatmap of recombinant_3-like LSDV genomes

Full size image

The recombination analysis of the LSDV genomes on RDP4 revealed 11 potential recombination events in the Indonesia LSDV genome, with LSDV_Neethling_Indoevaccine_LW_1959 (AF409138.1) and LSDV_NI-2490 (AF325528.1) as the parental donors. Similar recombination events are also present in other LSDV isolates from China, Thailand, Hong Kong, Vietnam, Taiwan, and Tomsk, Russia (Table 1). Some of the ORFs affected by the recombination events include the G-protein coupled receptor gene (LSDV011), the DNA polymerase gene (LSDV039), the DNA ligase-like protein gene (LSDV133), and the variola virus B22R-like protein gene (LSDV134). These four ORFs are among the nine proteins annotated as only 99.5-99.9% homologous to both LSDV_NI-2490 and LSDV_Neethling_vaccine_LW_1959 reference genomes (Additional file 1: Table s3). The most notable ORF was LSDV087, whose sequence analysis showed a nucleotide insertion at position 599 in all Neethling vaccine strains which caused a frameshift and truncated the ORF to 603nt instead of 763nts. No recombination was observed within this ORF, and all known LSDV recombinants were 100% identical to the field isolates in this region of the genome. Further, in addition to the 8 SNPs observed between the LSDV NI-like field isolates and Neethling-like field isolates, the nucleotide alignment of ORF087 showed 14 SNPs between the LSDV NI-like field isolates and Neethling derived vaccines; and 7 SNPs between Neethling-like field isolates and Neethling derived vaccines (Additional file 2: Figure s1).

Table 1 Predicted potential recombination events using different detection methods available in the RDP4 program
Full size table

The phylogenetic analysis of sequence fragments located on the 5’ and 3’ ends of breakpoint events supported the recombination events. For example, analysis of the sequences located on the 5’ end of the breakpoint event 10, clustered the Indonesian LSDVs with LSDV_NI-2490 (Cluster II), while the sequences on the 3’ end of the breakpoint clustered them with LSDV_Neethling_vaccine_LW_1959 (Cluster I) (Additional file 3: Figure s2). The recombination of Indonesian LSDVs was also confirmed using SimPlot, followed by phylogenetic analysis of the sequences located on the 5’ and 3’ ends of the breakpoints (Additional file 4: Figure s3). Similarly, the phylogenetic network analysis of the aligned LSDV genomes using PopART software indicated that the Indonesia LSDVs are likely recombinants of NI-like and Neethling-like LSDVs placed on the opposite ends of the network. The network confirmed that the LSDVs from Indonesia cluster together with other isolates from East and Southeast Asia (recomb_3-like), which have been demonstrated as LSDV recombinants and are different from the known LSDV recombinants from Russia (recomb_1-like, recomb_2-like, and recomb_4-like) (Fig. 5).

Fig. 5
figure 5

Median-joining network inferred from LSDV whole genome sequences using the PopART program. The network shows the Indonesia LSDVs (in red) in cluster Recomb_3-like, between the NI-like and Neethling-like LSDV clusters. The number of mutations between each genome is labeled

Full size image

Discussion

This study has presented typical clinical signs of LSD including skin nodules, nasal discharge, and high fever in suspected cattle cases in 2022 outbreak in Indragiri Hulu district, Indonesia, with morbidity and mortality rates of 28.04% and 1.87%, respectively. It has been reported that the severity of LSD varies based on the immune status of the animal, the breed, the age, and the virus strain, leading to an average mortality rate of < 10% and variable morbidity rate from 5 to 45% [2, 4]. However, the morbidity and mortality rates recorded in the LSD outbreak in Indonesia were consistent with other LSD outbreaks in Asian countries [28, 37, 38].

The initial molecular characterization of the Indonesian LSDV isolates based on targeted CaPV-marker genes indicated possible virus recombination [18,19,20,21]. First, the nucleotide sequence analysis of the four targeted genes showed that the Indonesian LSDV isolates were 100% identical to the recently demonstrated LSDV recombinants from East and Southeast Asia: China, Thailand, Taiwan, Vietnam, and HongKong [23, 32,33,34,35]. Second, on one hand, the complete RPO30 gene and the partial EEV gene sequence analysis showed that Indonesian LSDV sequences have 99.67% and 100% sequence similarity, respectively, to the LSDV Neethling-derived vaccines. On the other hand, the GPCR gene sequence similarity was highest (99.13%) to the common LSDV field strains, while the B22R sequence was 99.40% similar to the Neethling vaccine strain and 99.28% to the NI-2490 LSDV strain.

Further analysis of the Indonesian LSDV whole genome sequences confirmed its closest homology to the recombinant LSDVs from East and Southeast Asia and that it was 99.56% similar to Neethling_LSDV_vaccine and 99.25% to the LSDV_NI-2490 strain. Possible recombination of the Indonesian LSDV genomes was assessed and verified using different tools (RDP4, SimPlot, DAPCA, and PopART network) that have previously been applied in recombination analysis [39, 40]. The PopART network and the PCA scatter plot showed seven clusters of LSDVs and grouped Indonesia isolates with other reported East and Southeast Asia recombinant LSDVs [41]. The RDP4 and SimPlot indicated that Indonesia LSDVs are most likely recombinants of Neethling_LSDV_vaccine and LSDV_NI-2490 isolates with 11 recombination events. Notably, the recombinant events detected in this study have been reported in LSDV genomes from Thailand and China [32, 34]. Although the parental donors for all the reported LSDV recombinants are Neethling vaccine and Kenyan NI-2490-like strains, the Indonesia LSDV recombination was similar to other recombinants from East and Southeast Asia but different from those found in Russia that have 27 (LSDV_ Russia/Saratov/2017), 24 (LSDV_Udmurtiya_Russia_2019) and 25 (Tyumen_2019, OL542833_1) recombination events [42,43,44]. The difference in recombination patterns of LSDVs was demonstrated in this study by the SNP heatmap and the PopART network analysis that showed the phylogenetic relatedness of the recombinants to their parental donors.

In addition to the previously reported genes employed in combination in this study (GPCR, RPO30, EEV, and B22R), the LSDV whole genome analysis suggested that LSDV087 could be an additional DIVA target through PCR and sequencing, that can complement the existing CaPV RT-PCR DIVA methods, due to the unique signatures or SNPs found at DNA level in LSDV Neethling derived vaccines compared to the field strains including the LSDV recombinants [18,19,20,21]. In addition to the intergenic recombination (exchange of full-length gene sequence) observed in Indonesian LSDV genomes and other reported recombinant LSDVs, intragenic recombination (exchange of gene sequence fragments) was also detected in nine genes [45]. Recombination in viruses is common in the field following a co-infection of viruses in a single cell of a host which is influenced by the selection pressure of a virus population and contributes to virus evolution through genetic shift [45, 46]. The intragenic recombinations like those detected in the Indonesia LSDV isolates are related to the regulation of gene expression, hence altering the pathogenic properties of the virus [45, 46]. Virus recombination in the field is likely following vaccination with live attenuated vaccine in a population where wild-type virus is already circulating or contamination in cell culture during vaccine production [47]. No LSDV vaccination program was reported in Indonesia before the LSD outbreak; therefore, it is possible that the virus was introduced into the country from the neighboring countries. This is likely considering that the Indonesian LSDV isolate is nearly identical to the isolates circulating in the neighboring countries in East and Southeast Asia, suggesting a common origin of the virus in the region. A recent study has shown that the emergence of LSDV recombinants in Asia is most likely a result of a spillover from animals vaccinated with the Lumpivax vaccine that appears to be a mixture of Neethling-like LSDV vaccine strain, KSGP-0240 vaccine strain and several LSDV recombinants that are nearly identical to those circulating in Asia [41]. For instance, the initial outbreak of the recombinant_3-like LSDVs was reported in Qapqal Xibe county, Xinjiang province in Northwest China, which is located ∼ 20 km from the border of Kazakhstan, where a mass vaccination campaign with Lumpivax vaccine was launched [48]. Thus, it is likely that in addition to the poor-quality vaccine spillover, cattle movement, trade and insect vectors may have contributed to the spread of the LSDV recombinants into the neighboring countries in the region. Moreover, It would be interesting to further investigate whether the changes in proteins for viral virulence and host range, such as G-protein coupled receptor (LSDV011), IFN-alpha beta receptor glycoprotein (LSDV135), IL-10 (LSDV005) and IL-1 receptor (LSDV006), that were found in LSDV recombinants circulating in Indonesia and other Southeast Asian countries cause higher infectivity and virulence, thereby contributing to the virus’s rapid spread as suggested in earlier studies with LSDV recombinants in Russia [49, 50].

Conclusion

This study has reiterated the importance of applying several analytical tools for accurate LSDV whole genome analysis, especially with datasets including recombinants. Based on these tools, the study has shown that the Indonesia LSDV isolate is a recombinant of Neethling_vaccine_LW_1959 and LSDV_NI-2490 strains and is closely related to the previously reported LSDV recombinants circulating in East and Southeast Asia, suggesting a common origin/source of the LSDV isolate causing the outbreak and spreading in the region. Most importantly, the analysis in this study has identified a unique ORF, LSDV087, as a potential DIVA target for LSDVs through sequencing. Therefore, this study demonstrates the significance of accurate molecular characterization of LSDV for disease surveillance and understanding of its origin, as well as vaccine quality control before animal vaccination.

Materials and methods

Sample collection and processing

Six skin nodules and six nasal swab samples were collected from suspected cattle cases in Indragiri Hulu district in 2022 and sent to the RCVS for further investigation. DNA was extracted using the Qiagen AllPrep DNA/RNA extraction kit and eluted in 50 Î¼l elution buffer.

PCR and differential diagnosis

A real-time PCR (RT-PCR), previously described by Bowden et al. [6], was used to detect the CaPV genome in all twelve collected samples. The PCR reactions were prepared in a 20 Î¼l reaction volume containing 18 Î¼l of the PCR master mix with primers, a probe, and two μl of template DNA. The positive samples were characterized further using the high-resolution melting (HRM) assay capable of differentiating CaPVs based on the melting temperature (Tm) of the PCR amplicon [51].

Amplification and sequencing of selected CaPV genes

Four selected CaPV-marker genes (RPO30, GPCR, EEV, and B22R) that can differentiate LSDV vaccine strains from field isolates were amplified and sequenced for the CaPV-positive samples following the previously described PCR conditions [52]. The PCR products were analyzed on a 2% agarose gel in 1X TAE buffer and visualized on a Gel Documentation System (Bio-Rad). The PCR products were purified and sent for Sanger sequencing from both directions using forward and reverse primers at LGC Genomics (Germany).

Whole genome sequencing of LSDV using Ion S5 technology

Approximately 100 ng of DNA was enzymatically fragmented into 200 bp lengths using Ion shear Plus reagents. The fragmented DNA was ligated to adapters and barcodes using the Ion Xpressâ„¢ Plus Fragment Library Kit and the Ion Xpress barcode adapters (Thermo Fisher Scientific, USA) following the manufacturer’s instructions. Size selection was performed using Pippin Prep (Sage Science, Inc., USA). The libraries were further amplified for eight cycles using Platinumâ„¢ PCR SuperMix high fidelity and library amplification primer mix provided with the Ion Plus Fragment Library Kit. The amplified barcoded libraries at a concentration of 100 pM were pooled in equal volumes and loaded onto the Ion Chefâ„¢ Instrument (Thermo Fisher Scientific, USA) for automated template preparation and chip loading using Ion 540â„¢ Kit-Chef (Thermo Fisher Scientific, USA). Using the Ion Chefâ„¢ Instrument, the pooled libraries were clonally amplified on the Ion spheres (ISPs) by emulsion PCR, followed by automated loading template-enriched ISPs onto an Ion 540 chip. Sequencing was performed on an Ion S5â„¢ Next generation sequencing system (Thermo Fisher Scientific, USA) with 500 flows to generate 200 bp reads.

Sequence and phylogenetic analysis of CaPV-marker genes

The obtained sequences were assembled using Vector NTI software (Invitrogen) v11.5. For each targeted gene, the sequences were analyzed with those of other LSDVs, SPPVs, and GTPVs retrieved from GenBank. Multiple sequence alignments were performed on MEGA X using the Muscle algorithm and the codon option. Neighbor-joining trees for the complete RPO30 and GPCR genes were also built on MEGA X, and the evolutionary distances were computed using the Maximum Composite Likelihood method with 1000 bootstrap replicates [53]. The phylogenetic trees with the associated metadata were visualized on the Interactive Tree of Life (ITOL) tool [54]. The multiple sequence alignments, analysis, and visualization of the partial EEV glycoprotein and B22R genes were conducted using BioEdit v7.2.5.

Whole genome reconstruction and annotation of the Indonesian LSDV isolate

The run was pre-processed using the torrent suite software to remove adapters. Then, the raw reads were quality-filtered using fastq-mcf v1.04.676 (ea-utils), and their quality was assessed with FastQC (v. 011.5). A Phred score cut-off of 20 was used, and the reads with lengths between 50 and 250 bp were selected. Denovo assemblies were performed on a subset of reads using SPAdes (v3.11.1), Unicycler (v0.5.0), and Megahit (v1.2.9). Reference-guided assembly was performed using bowtie2 (2.3.4.1). After mapping the cleaned raw reads against the reference sequence, the Mpileup files were generated using SAMtools (v1.11). Variant calling was performed using the consensus caller method of BCFtools (v1.9). The consensus sequences were produced with vcfutils.pl (VCFtools v0.1.16) and seqtk (v1.3.106) using a Phred score of 20 and compared to the de novo assemblies after alignment with Mafft (v7.453). The assemblies were also visualized using IGV, and mapping statistics were visualized with Qualimap (v.2.2.1).

Given that the Indonesia LSDVs were grouped in LSDV Cluster I or II based on different LSDV-marker genes analysis, the open reading frames (ORFs) of the assembled Indonesia LSDV genomes were predicted with GATU using reference genomes from LSDV Cluster I; LSDV_Neethling_vaccine_LW_1959 (AF409138.1) and Cluster II; LSDV_NI-2490 (AF325528.1) [55]. The whole genome sequences of Indonesia LSDV isolates, LSDV_Indonesia_2022_S1, and LSDV_Indonesia_2022_S4, were deposited in the GenBank database under the accession numbers OR232413 and OR232414, respectively.

LSDV whole genome SNP and recombination analysis

The assembled Indonesian LSDV genomes were aligned with 46 complete LSDV genome sequences retrieved from GenBank using MAFFT. The SNPs from the aligned LSDV dataset relative to the LSDV_NI-2490 (AF325528.1) genome were analyzed using the adegenet package in R. First, the alignment file was converted into a genlight object using the function fasta2genlight and the SNPs were extracted [56]. The SNPs distribution and density of the LSDV genotypes were visualized by heatmap. The discriminant analysis of Principal Component Analysis (DAPCA) of the extracted SNPs was also performed to identify the different clusters of the LSDV genotypes, which were plotted as a PCA scatterplot and neighbor-joining tree. Further SNP-based analysis was performed on all LSDV genomes clustering together with the Indonesia LSDV isolates. Three different tools were applied to assess evidence of recombination in the Indonesia LSDV isolate. The potential recombination events of the Indonesia LSDV isolates were first evaluated using the RDP4 software, which includes seven methods (RDP, GENECONV, BootScan, MaxChi, Chimera, SiScan, and 3Seq) [57]. The recombination was then analyzed using the SimPlot 3.5.1 software with a window size of 1000 bp and a step of 20 bp. Given that recombinant genome sequences are not suitable for whole genome phylogenetic analysis [25, 27], several phylogenetic trees were constructed based on the fragment sequences before and after the breakpoints predicted by RDP4 and SimPlot software on MEGA X, using the neighbor-joining method with 1000 bootstraps and visualized on ITOL. Lastly, the recombination was explored by conducting a phylogenetic network analysis of the aligned LSDV genomes on the PopART program, using the median-joining method with epsilon set to zero [40].

Data availability

The DNA sequences generated and used in the analysis for this study are available in the Gen-Bank under the accession numbers OR246967 to OR246986 for gene fragments and OR232413 and OR232414 for the whole genomes.

Abbreviations

BCF:

Binary variant call format

IGV:

Integrative genomics viewer

LSDV:

Lumpy skin disease virus

PCA:

Principal Component Analysis

SAM:

Sequence Alignment/Map

SNP:

Single nucleotide polymorphism

VCF:

Variant call format

References

  1. Buller RM, Arif BM, Black DN, Dumbell KR, Esposito JJ, Lefkowitz EJ, et al. Family poxviridae. Virus taxonomy: classification and nomenclature of viruses. Eighth report of the International Committee on Taxonomy of viruses. San Diego, CA, USA: Elsevier Academic; 2005. pp. 117–33.

    Google Scholar 

  2. Tuppurainen ESM, Venter EH, Shisler JL, Gari G, Mekonnen GA, Juleff N, et al. Review: Capripoxvirus diseases: current status and opportunities for control. Transbound Emerg Dis. 2017;64:729–45.

    Article  CAS  PubMed  Google Scholar 

  3. Tulman ER, Afonso CL, Lu Z, Zsak L, Sur J-H, Sandybaev NT et al. The genomes of Sheeppox and Goatpox viruses. J Virol. 2002;76.

  4. Babiuk S, Bowden TR, Boyle DB, Wallace DB, Kitching RP. Capripoxviruses: an emerging worldwide threat to sheep, goats and cattle. Transbound Emerg Dis. 2008;55:263–72.

    Article  CAS  PubMed  Google Scholar 

  5. Brenner J, Bellaiche M, Gross E, Elad D, Oved Z, Haimovitz M et al. Appearance of skin lesions in cattle populations vaccinated against lumpy skin disease: statutory challenge. Vaccine. 2009;27.

  6. Bowden TR, Babiuk SL, Parkyn GR, Copps JS, Boyle DB. Capripoxvirus tissue tropism and shedding: a quantitative study in experimentally infected sheep and goats. Virology. 2008;371:380–93.

    Article  CAS  PubMed  Google Scholar 

  7. Tuppurainen ESM, Lubinga JC, Stoltsz WH, Troskie M, Carpenter ST, Coetzer JAW et al. Evidence of vertical transmission of lumpy skin disease virus in Rhipicephalus decoloratus ticks. Ticks Tick Borne Dis. 2013;4.

  8. Sohier C, Haegeman A, Mostin L, De Leeuw I, Campe W, De Van A et al. Experimental evidence of mechanical lumpy skin disease virus transmission by Stomoxys calcitrans biting flies and Haematopota spp. Horseflies Sci Rep. 2019;9.

  9. Chihota CM, Rennie LF, Kitching RP, Mellor PS. Mechanical transmission of lumpy skin disease virus by Aedes aegypti (Diptera: Culicidae). Epidemiol Infect. 2001;126.

  10. Tuppurainen ESM, Oura CAL, Review. Lumpy skin disease: an emerging threat to Europe, the Middle East and Asia. Transbound Emerg Dis. 2012;59:40–8.

    Article  CAS  PubMed  Google Scholar 

  11. Annandale CH, Holm DE, Ebersohn K, Venter EH. Seminal transmission of lumpy skin disease virus in heifers. Transbound Emerg Dis. 2014;61.

  12. Rouby S, Aboulsoud E. Evidence of intrauterine transmission of lumpy skin disease virus. Vet J. 2016;209.

  13. Aleksandr K, Olga B, David WB, Pavel P, Yana P, Svetlana K et al. Non-vector-borne transmission of lumpy skin disease virus. Sci Rep. 2020;10.

  14. Lafar S, Zro K, Ennaji MM. Capripoxvirus diseases: current updates and developed strategies for control. In: Emerging and reemerging viral pathogens: volume 1: fundamental and basic virology aspects of human, animal and plant pathogens. 2019.

  15. Wolff J, Moritz T, Schlottau K, Hoffmann D, Beer M, Hoffmann B. Development of a safe and highly efficient inactivated vaccine candidate against lumpy skin disease virus. Vaccines. 2021;9.

  16. Tuppurainen E, Dietze K, Wolff J, Bergmann H, Beltran-alcrudo D, Fahrion A et al. Review: vaccines and vaccination against lumpy skin disease. Vaccines. 2021;9.

  17. Haegeman A, De Leeuw I, Saduakassova M, Van Campe W, Aerts L, Philips W et al. The importance of quality control of lsdv live attenuated vaccines for its safe application in the field. Vaccines. 2021;9.

  18. Lamien CE, Le Goff C, Silber R, Wallace DB, Gulyaz V, Tuppurainen E, et al. Use of the Capripoxvirus homologue of Vaccinia virus 30 kDa RNA polymerase subunit (RPO30) gene as a novel diagnostic and genotyping target: development of a classical PCR method to differentiate Goat Poxvirus from Sheep Poxvirus. Vet Microbiol. 2011;149:30–9.

    Article  CAS  PubMed  Google Scholar 

  19. Le Goff C, Lamien CE, Fakhfakh E, Chadeyras A, Aba-Adulugba E, Libeau G, et al. Capripoxvirus G-protein-coupled chemokine receptor: a host-range gene suitable for virus animal origin discrimination. J Gen Virol. 2009;90:1967–77.

    Article  PubMed  Google Scholar 

  20. Mishra B, Mondal P, Patel CL, Zafir I, Gangwar R, Singh N, et al. VARV B22R homologue as phylogenetic marker gene for Capripoxvirus classification and divergence time dating. Virus Genes. 2019;55:51–9.

    Article  CAS  PubMed  Google Scholar 

  21. Menasherow S, Erster O, Rubinstein-Giuni M, Kovtunenko A, Eyngor E, Gelman B, et al. A high-resolution melting (HRM) assay for the differentiation between Israeli field and Neethling vaccine lumpy skin disease viruses. J Virol Methods. 2016;232:12–5.

    Article  CAS  PubMed  Google Scholar 

  22. Ma J, Yuan Y, Shao J, Sun M, He W, Chen J et al. Genomic characterization of lumpy skin disease virus in southern China. Transbound Emerg Dis. 2022;69.

  23. Flannery J, Shih B, Haga IR, Ashby M, Corla A, King S et al. A novel strain of lumpy skin disease virus causes clinical disease in cattle in Hong Kong. Transbound Emerg Dis. 2022;69.

  24. Schierup MH, Hein J. Consequences of recombination on traditional phylogenetic analysis. Genetics. 2000;156.

  25. Huson DH, Bryant D. Application of phylogenetic networks in evolutionary studies. Mol Biol Evol. 2006;23.

  26. Rasmussen MD, Hubisz MJ, Gronau I, Siepel A. Genome-wide inference of ancestral recombination graphs. PLoS Genet. 2014;10.

  27. Wohl S, Schaffner SF, Sabeti PC. Genomic analysis of viral outbreaks. Annual Rev Virol. 2016;3.

  28. Azeem S, Sharma B, Shabir S, Akbar H, Venter E. Lumpy skin disease is expanding its geographic range: a challenge for Asian livestock management and food security. Vet J. 2022;279.

  29. Badhy SC, Chowdhury MGA, Settypalli TBK, Cattoli G, Lamien CE, Fakir MAU et al. Molecular characterization of lumpy skin disease virus (LSDV) emerged in Bangladesh reveals unique genetic features compared to contemporary field strains. BMC Vet Res. 2021;17.

  30. Koirala P, Meki IK, Maharjan M, Settypalli BK, Manandhar S, Yadav SK et al. Molecular characterization of the 2020 outbreak of lumpy skin disease in Nepal. Microorganisms. 2022;10.

  31. Sudhakar SB, Mishra N, Kalaiyarasu S, Jhade SK, Singh VP. Genetic and phylogenetic analysis of lumpy skin disease viruses (LSDV) isolated from the first and subsequent field outbreaks in India during 2019 reveals close proximity with unique signatures of historical Kenyan NI-2490/Kenya/KSGP-like field strains. Transbound Emerg Dis. 2022;69.

  32. Li L, Wang Z, Qi C, Liu S, Gong M, Li J et al. Genetic analysis of genome sequence characteristics of two lumpy skin disease viruses isolated from China. BMC Vet Res. 2022;18.

  33. Huang C-W, Ting L-J, Liu Y-P, Lin Y-J, Lee F, Chiou C-J. Complete coding sequence of lumpy skin disease virus isolated from Kinmen Island, Taiwan, in 2020. Microbiol Resour Announc. 2022;11.

  34. Suwankitwat N, Songkasupa T, Boonpornprasert P, Sripipattanakul P, Theerawatanasirikul S, Deemagarn T et al. Rapid spread and genetic characterisation of a recently emerged recombinant lumpy skin disease virus in Thailand. Vet Sci. 2022;9.

  35. Mathijs E, Vandenbussche F, Nguyen L, Aerts L, Nguyen T, De Leeuw I et al. Coding-complete sequences of recombinant lumpy skin disease viruses collected in 2020 from four outbreaks in Northern Vietnam. Microbiol Resour Announc. 2021;10.

  36. WOAH. Indonesia - lumpy skin disease virus (Inf. with) - immediate notification. Retrieved from https://wahis.woah.org/#/in-review/4340. 2022.

  37. Arjkumpa O, Suwannaboon M, Boonrawd M, Punyawan I, Laobannu P, Yantaphan S et al. First emergence of lumpy skin disease in cattle in Thailand, 2021. Transbound Emerg Dis. 2021;68.

  38. Roche X, Rozstalnyy A, TagoPacheco D, Kamata A, Pittiglio C, Alcrudo DB et al. Introduction and spread of lumpy skin disease in South, East and Southeast Asia: qualitative risk assessment and management. FAO Anim Prod Heal. 2020;183.

  39. Salminen M, Martin D. Detecting and characterizing individual recombination events. In: The Phylogenetic Handbook. 2012.

  40. Leigh JW, Bryant D. POPART: full-feature software for haplotype network construction. Methods Ecol Evol. 2015;6.

  41. Vandenbussche F, Mathijs E, Philips W, Saduakassova M, De Leeuw I, Sultanov A et al. Recombinant LSDV strains in Asia: vaccine spillover or natural emergence? Viruses. 2022;14.

  42. Kononov A, Byadovskaya O, Kononova S, Yashin R, Zinyakov N, Mischenko V et al. Detection of vaccine-like strains of lumpy skin disease virus in outbreaks in Russia in 2017. Arch Virol. 2019;164.

  43. Krotova A, Byadovskaya O, Shumilova I, Zinyakov N, van Schalkwyk A, Sprygin A. Molecular characterization of a novel recombinant lumpy skin disease virus isolated during an outbreak in Tyumen, Russia, in 2019. Transbound Emerg Dis. 2022;69.

  44. Sprygin A, Van Schalkwyk A, Shumilova I, Nesterov A, Kononova S, Prutnikov P, et al. Full-length genome characterization of a novel recombinant vaccine-like lumpy skin disease virus strain detected during the climatic winter in Russia, 2019. Arch Virol. 2020;165:2675–7.

    Article  CAS  PubMed  Google Scholar 

  45. Pérez-Losada M, Arenas M, Galán JC, Palero F, González-Candelas F. Recombination in viruses: mechanisms, methods of study, and evolutionary consequences. Infection, Genetics and Evolution. 2015;30.

  46. Gershon PD, Kitching RP, Hammond JM, Black DN. Poxvirus genetic recombination during natural virus transmission. J Gen Virol. 1989;70.

  47. Gundlach BR, Lewis MG, Sopper S, Schnell T, Sodroski J, Stahl-Hennig C et al. Evidence for recombination of live, attenuated immunodeficiency virus vaccine with challenge virus to a more virulent strain. J Virol. 2000;74.

  48. Lu G, Xie J, Luo J, Shao R, Jia K, Li S. Lumpy skin disease outbreaks in China, since 3 August 2019. Transbound Emerg Dis. 2021;68:216–9.

    Article  PubMed  Google Scholar 

  49. Kononova S, Kononov A, Shumilova I, Byadovskaya O, Nesterov A, Prutnikov P, et al. A lumpy skin disease virus which underwent a recombination event demonstrates more aggressive growth in primary cells and cattle than the classical field isolate. Transbound Emerg Dis. 2021;68:1377–83.

    Article  CAS  PubMed  Google Scholar 

  50. Nesterov A, Mazloum A, Byadovskaya O, Shumilova I, Van Schalkwyk A, Krotova A et al. Experimentally controlled study indicates that the naturally occurring recombinant vaccine-like lumpy skin disease strain Udmurtiya/2019, detected during freezing winter in northern latitudes, is transmitted via indirect contact. Front Vet Sci. 2022;9.

  51. Gelaye E, Mach L, Kolodziejek J, Grabherr R, Loitsch A, Achenbach JE et al. A novel HRM assay for the simultaneous detection and differentiation of eight poxviruses of medical and veterinary importance. Sci Rep. 2017;7.

  52. Chibssa TR, Sombo M, Lichoti JK, Adam TIB, Liu Y, Elraouf YA et al. Molecular analysis of east African lumpy skin disease viruses reveals a mixed isolate with features of both vaccine and field isolates. Microorganisms. 2021;9.

  53. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms| Molecular Biology and Evolution| Oxford Academic. Mol Biol Evol. 2018;35.

  54. Letunic I, Bork P. Interactive tree of life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res. 2021;49:W293–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Tcherepanov V, Ehlers A, Upton C. Genome annotation transfer utility (GATU): rapid annotation of viral genomes using a closely related reference genome. BMC Genomics. 2006;7.

  56. Jombart T, Ahmed I. Adegenet 1.3-1: new tools for the analysis of genome-wide SNP data. Bioinformatics. 2011;27.

  57. Martin DP, Murrell B, Golden M, Khoosal A, Muhire B. RDP4: detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015;1.

Download references

Acknowledgements

The authors thank all the staff from the Indonesian Research Center for Veterinary Science, Research Organization for Health, National Research and Innovation Agency for collecting and providing the samples for this study.

Funding

This research was supported by the VETLAB network initiative of the Joint FAO/IAEA Centre of Nuclear Techniques in Food and Agriculture, funded through the Peaceful Uses Initiatives (PUI) by Japan and the United States of America.

Author information

Authors and Affiliations

  1. Research Center for Veterinary Science, Research Organization for Health, National Research and Innovation Agency, West Java, Indonesia

    Indrawati Sendow, Ni Luh Putu Indi Dharmayanti, Heri Hoerudin, Atik Ratnawati, Harimurti Nuradji, Muharam Saepulloh, Rahmat Setya Adji & Nuha Fairusya

  2. Animal Production and Health Laboratory, Joint FAO/IAEA Centre of Nuclear Techniques in Food and Agriculture, Department of Nuclear Sciences and Applications, International Atomic Energy Agency, Wagramer Strasse 5, A-1400, Vienna, P.O. Box 100, Austria

    Irene Kasindi Meki, Tirumala Bharani K. Settypalli, Hatem Ouled Ahmed, Giovanni Cattoli & Charles Euloge Lamien

  3. Regional Livestock Services, Riau Province, Indonesia

    Faralinda Sari

  4. Disease Investigator Center Bukit Tinggi, West Sumatera, Indonesia

    Katamtama Anindita

Authors
  1. Indrawati Sendow
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Irene Kasindi Meki
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Ni Luh Putu Indi Dharmayanti
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Heri Hoerudin
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Atik Ratnawati
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Tirumala Bharani K. Settypalli
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Hatem Ouled Ahmed
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Harimurti Nuradji
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Muharam Saepulloh
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Rahmat Setya Adji
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Nuha Fairusya
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Faralinda Sari
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Katamtama Anindita
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Giovanni Cattoli
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. Charles Euloge Lamien
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

Experimental design and conceptualization: IS, NLPID, and CEL; Sample collection: HN, MS, RSA, NF, FS, and KA; Investigation: AR, HH, TBKS and HOA; Data analysis: CEL and IKM; Writing—original draft: IKM; Supervision: IS, CEL and GC; Resources: IS, CEL and GC; All authors reviewed and edited the manuscript.

Corresponding author

Correspondence to Irene Kasindi Meki.

Ethics declarations

Ethical approval and consent to participate

The authors confirm that the ethical policies of the journal, as stated on the journal’s ethical guidelines have been adhered to. The samples used in this study were collected in animal farms after approval by all owners and submitted to the Research Centre for Veterinary Science, National Research and Innovation Agency, in Indonesia for LSDV diagnosis. Ethical approval was not required. This study was conducted in accordance with ARRIVE guidelines.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sendow, I., Meki, I.K., Dharmayanti, N.L.P.I. et al. Molecular characterization of recombinant LSDV isolates from 2022 outbreak in Indonesia through phylogenetic networks and whole-genome SNP-based analysis. BMC Genomics 25, 240 (2024). https://doi.org/10.1186/s12864-024-10169-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-024-10169-6

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us