Skip to main content

BRIEF RESEARCH REPORT article

Front. Vet. Sci., 12 January 2023
Sec. Veterinary Epidemiology and Economics
This article is part of the Research Topic Sequencing and Phylogenetic Analysis as a Tool in Molecular Epidemiology of Veterinary Infectious Diseases View all 19 articles

Retrospective genomic analysis of the first Lumpy skin disease virus outbreak in China (2019)

\nYu-Rong WeiYu-Rong Wei1Wen-Ge MaWen-Ge Ma1Ping WangPing Wang1Wen WangWen Wang2Xiao-Hui SuXiao-Hui Su2Xue-Yun YangXue-Yun Yang1Xiao-Yun MiXiao-Yun Mi1Jian-Yong WuJian-Yong Wu1Jiong Huang
Jiong Huang1*
  • 1Xinjiang Key Laboratory of Animal Infectious Diseases, Institute of Veterinary Medicine, Xinjiang Academy of Animal Science, Urumqi, China
  • 2Center for Animal Disease Prevention and Control of Xinjiang, Urumqi, China

Lumpy skin disease caused by Lumpy skin disease virus (LSDV) is a severe systemic disease affecting cattle and other ruminants. Lumpy skin disease was first reported in northwest China in August 2019 and has severely threatened the cattle breeding industry in China. However, there have been limited genomic studies of LSDV from the first outbreak and its subsequent epidemics. This study aims to characterize the comparative genomic evolution of the LSDV strain from the first outbreak in China. The etiological agent was isolated in a Madin-Darby bovine kidney cell culture and subsequently identified by PCR and Sanger sequencing of six selected genes. The genome sequence was determined using Illumina sequencing and analyzed through genome alignment and phylogenetic tree. The results showed that all six genes were successfully amplified and genetically clustered into LSDV. The virus presented the highest homology to strain China/GD01/2020, which shared 100% identities among 150 open reading frames (ORFs), and 97.1–99.7% identities among additional 6 ORFs. Bayesian inference tree analysis revealed that the virus shared a common ancestor with LSDV strains from China and Vietnam. The study provides an additional genomic data for LSDV tracking and control in China and neighboring countries.

1. Introduction

Lumpy skin disease (LSD) is a significant transboundary viral disease that affects cattle water buffalos and giraffe, and other ruminants (14). The disease is caused by the Lumpy skin disease virus (LSDV), a member of the family Poxviridae, genus Capripoxvirus (5, 6). The virus appears to be mechanically transmitted by blood-sucking arthropods such as flies, mosquitoes and ticks, and to a lesser extent by direct contact between cattle (7, 8). In addition, LSDV-contaminated milk, blood, nasal secretions, and saliva are alternative sources of infection through feeding or drinking routes (9). The affected animals mainly manifest fevers and nodular lesions, and they produce dramatically less milk and undergo weight loss (10). LSDV can cause a high incidence of 5–45% when introduced into a herd, and the case fatality rate ranges from 0.5 to 7.0% (11, 12). As a result, LSD poses a significant economic threat to the cattle-breeding industry.

LSD was initially described in Zambia in 1929 and identified as a communicable disease in the 1940's (6). In the 1950's, LSDV spread rapidly through Central and Eastern Africa. It then spread out of Africa into the Middle East in the 1980's (13). Since 2015, LSD outbreaks have occurred in countries neighboring China, such as Kazakhstan, and Russia (14, 15). In August 2019, the first outbreak of LSD was reported in the Xinjiang Uygur Autonomous Region, Northwest China, in which borders Kazakhstan (16). The disease was then reported in eight provincial administrative regions (Anhui, Fujian, Hong Kong, Guangdong, Jiangxi, Sichuan, Taiwan, and Zhejiang), resulting in 10 LSDV outbreaks in China (1618). From January 2013 to July 2021, there were 28,442 LSDV outbreaks worldwide, resulting in 326,300 cases and 15,500 deaths (19). The global dissemination of LSDV has resulted in serious risk of this contagious disease affecting the large cattle population (more than 95.6 million cattle in stock) in China.

The genome of LSDV is a linear double-stranded 145–152 kb DNA molecule that contains 150–156 predicted open reading frames (ORFs). The first complete genome sequence was determined in 2001 from primary lamb testicle cells of the Neethling type strain 2,490, which contained 156 annotated genes (20). Of these encoded genes, the G protein-coupled chemokine receptor (GPRC) and the RNA polymerase 30 kDa subunit (RP030) were recognized as markers for differentiating the poxviruses at the family and genus levels (21, 22). In addition, the LSDV contains 90 core genes conserved in all chordopoxviruses, and they have been used for phylogenetic analyses (17). Since 2001, approximately 40 LSDV genome sequences of different origins have been sequenced and deposited in GenBank, but comparative genomic data and evolutionary studies are still limited.

The genomes of Chinese strains from Guangdong and Hong Kong in South China have been sequenced (17, 23), and previous studies of these circulating Chinese LSDV strains showed close genetic relationship with the LSDV/Russia/Saratov/2017 (accession no. MH646674.1) strain or Neethling vaccine strain based on single or multiple genes (11). Thus, the genetic relationships between the LSDV isolate present in the first outbreak in China and those of subsequent epidemics in China and other countries still needs to be clarified. In this study, we employed next-generation sequencing to obtain the complete genome sequence of the LSDV/China/XJ01/2019 strain isolated from the only cow that died during the first LSDV outbreak in China. We also performed a detailed genomic comparison of LSDV/China/XJ01/2019 and related genomic sequences reported before and after this outbreak. This study provides insight into the spread of LSDV during the epidemic in China.

2. Materials and methods

2.1. Specimens

In August 2019, the first ever LSD outbreak occurred in Ili Kazakh Autonomous Prefecture, Xinjiang Uygur Autonomous Region, Northwest China, which shares an approximately 50-km border with Kazakhstan (Supplementary Figure S1). During the outbreak, a Holstein cow was found to have died from LSD. After an examination and dissection of the dead dairy cow, a skin nodule sample was collected. It was then transported and shipped to the laboratory under cold conditions and immediately stored at −80°C for further testing. The owner of the animal was informed about the purpose and process of this study. The farmer agreed to allow the skin nodule samples to be collected from his died cow.

2.2. DNA extraction and amplification

The virus was cultured in Madin-Darby bovine kidney cell, and after three generations of culturing, the viral genome was extracted by using a QIAamp DNA Mini Kit (QIAGEN). The RPO30 and GPRC genes were amplified using the primers listed in Supplementary Table S1. Four additional primers (containing deletion or insertion sequences compared with Goatpox and Sheeppox viruses) were designed to amplify Ankyrin repeat protein (LSDV152), Interleukin-1 receptor-like protein (LSDV013), Putative alpha amanitin-sensitive protein (LSDV009) and Putative late transcription factor (LSDV076) genes that can be used to differentiate the genus Capripoxvirus referred to as the strain, LSDV/Russia/Saratov/2017 (accession no. MH646674.1) (Supplementary Table S1). These six genes were purified using Quick Gel columns (QIAGEN), and then ligated into pMD19-T (Tiangen Biotech) and transformed into Escherichia coli DH5α competent cells (Tiangen Biotech). Quintuplicate positive clones were extracted for Sanger sequencing (Sangon Biotech). Nucleotide sequences of the above six genes were downloaded from the National Center for Biotechnology Information, USA (Supplementary Table S2). Phylogenetic relatedness analyses were carried out using the MEGA 11 (https://megasoftware.net/) using the Maximum Likelihood method and the best fitting DNA model with 1,000 bootstrap replicates (24).

2.3. Genome sequencing and analysis

The LSDV genomic DNA was used for library construction and next-generation sequencing (Novagene). The raw reads were processed using a standard in-house pipeline (Novagene) to remove adapters, host sequences, chimeras, short reads, and low-quality reads. The clean data were further mapped to LSDV/Russia/Saratov/2017 (accession no. MH646674.1) by using Geneious prime software (25). The resulting contigs were assembled into a whole-genome sequence, which was then mapped onto a reference genome, resulting in a draft genome sequence. The clean reads were mapped to the gaps between the draft genome sequence and LSDV/Russia/Saratov/2017 (accession no. MH646674.1) using the medium sensitivity/Fast mode and iterated up to five times using Geneious prime 2020.0.3 software (25). The mapped sequences were used to generate consensus sequences to obtain the primary genomic sequence, which was then manually checked. Genome annotation were performed using GATU software (26) with the 20L81_Bang-Thanh VNM 20 (accession no. MZ577076.1) and Kubash KAZ 16 (accession no. MN642592.1) genomes as reference. The annotations were manually verified and curated using the Ugene software package (27).

The genome sequence of the strain LSDV/China/XJ01/2019 was aligned to a set of reference LSDV sequences retrieved from GenBank using ClustalX 2.1 (http://www.clustal.org/clustal2/). A phylogenetic tree was generated using the alignment and Bayesian approaches in MrBayes v. 3.2.7 (28) to evaluate the relationships between LSDV/China/XJ01/2019 and reference genome sequences in GenBank (Supplementary Table S3). Phylogenetic reconstruction was performed using the GTR evolutionary model including a Γ distribution and two runs of four chains each. The chain convergence was evaluated after 200,000 generations. Results were considered stable if an EES value was >200.

3. Results

3.1. Virus identification

The GPCR and PRO30 genes were successfully amplified from the cell culture, primers for four other genes (Ankyrin repeat protein, Interleukin-1 receptor-like protein, Putative alpha amanitin-sensitive protein, and Putative late transcription factor) were designed to confirm the presence of LSDV, and finally, all six genes were tested as positive. After Sanger sequencing and phylogenetic analysis, all six genes were clustered with LSDV; consequently, the virus was identified as LSDV (Figure 1) and named LSDV/China/XJ01/2019.

FIGURE 1
www.frontiersin.org

Figure 1. Phylogenetic diagram of RPO30 (A), GPCR (B), Putative alpha amanitin-sensitive protein (C), Ankyrin repeat protein gene (D), Interleukin-1 receptor-like protein gene (E) and Putative late transcription factor (F) genes of Lumpy skin disease virus. The phylogenetic tree was constructed using MEGA 11 with the Maximum Likelihood method. The Lumpy skin disease virus, Sheeppox virus, Goatpox virus are labeled with red, blue and green colors, respectively. The sequences obtained from this study are indicated by solid red triangles.

3.2. Genome assembly

The complete genome sequence of the LSDV/China/XJ01/2019 strain was analyzed using an Illumina NovaSeq sequencer (Illumina, USA) generating 150-bp single reads. A total of 7,496,974 (150 × 150 bp) PE150 clean reads were obtained from the Illumina NovaSeq sequencer. Consensus sequences were generated by de novo assembling and mapping to LSDV/Russia/Saratov/2017 (accession no. MH646674.1). The clean reads were subsequently mapped to the consensus sequences and the average coverage was determined to be 728.6 × (Supplementary Figure S1). The viral genome sequence was then submitted and deposited in GenBank under accession no. OM105589.

3.3. Genomic comparisons between LSDV/China/XJ01/2019 and other strains in and around China

Pairwise genome sequence comparison revealed that LSDV/China/XJ01/2019 strain shared the highest similarity (99.9960%) with 20L43_Ly-Quoc/VNM/20, followed by 20L43_Ly-Quoc/VNM/20 (99.9954%), 20L42_Quyet-Thang/VNM/20 (99.9947%), 20L81_Bang-Thanh/VNM/20 (99.9947%) and China/GD01/2020 (99.9893%). The mVista program (http://genome.lbl.gov/vista/mvista/submit.shtml) was used to analyze the genome-wide differences among the strains in and around China. Seven regions of the viral genome were extremely variable, and it contained deletions and mutation mainly in genes LSDV008 (Putative soluble interferon-gamma receptor gene), LSDV011 (G protein-coupled chemokine receptor-like protein gene), LSDV126 (putative EEV glycoprotein gene), LSDV145 (Ankyrin repeat protein gene) and LSDV146 (Phospholipase D-like protein gene). Specifically, complete genome sequences of Chinese strains LSDV/China/XJ01/2019, China/GD01/2020, and LSDV/HongKong/2020 are greatly similar to the genomes of the Vietnamese strains 20L42_Quyet-Thang/VNM/20, 20L43_Ly-Quoc/VNM/20, 20L70_Dinh-To/VNM/20 and 20L81_Bang-Thanh/VNM/20, with only three significantly different regions observed (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. Multiple genomic sequence alignment of the LSDV strains in and around China. Graph showing sequence similarities between LSDV/China/XJ01/2019 and the strains 20L42_Quyet-Thang/VNM/20, 20L43_Ly-Quoc/VNM/20, 20L70_Dinh-To/VNM/20, 20L81_Bang-Thanh/VNM/20, China/GD01/2020, Kubash/KAZ/16, LSDV/China/XJ01/2019, LSDV/HongKong/2020, LSDV/Russia/Dagestan/2015, LSDV/Russia/Saratov/2017 and LSDV/ Russia/Udmurtiya/2019, which are plotted in a sliding 100-bp window. The seven differentiated regions between the genomes of the Chinese strains and reference strain LSDV/Russia/Saratov/2017 are indicated by white dotted boxes.

3.4. Open reading frame comparisons of the LSDV/China/XJ01/2019 with other strains in and around China

The genome sequence was annotated referred to the strain 20L81_Bang-Thanh VNM 20 (accession no. MZ577076.1), a total of 156 open reading frames (ORFs) were identified in the strain, LSDV/China/XJ01/2019. Comparing the ORFs of LSDV/China/XJ01/2019 with the most closely related strain China/GD01/2020 revealed 150 ORFs sharing 100% sequence identities and 6 ORFs sharing 97.1–99.7% sequence identities (Table 1).

TABLE 1
www.frontiersin.org

Table 1. ORFs that are not identical between LSDV/China/XJ01/2019 and other Chinese, Kazakhstan, Russian and Vietnamese Lumpy skin disease virus strains.

A total of 17 variable loci in the genomic sequence, resulting in variations in 7 ORFs and 2 non-coding regions, were found between LSDV/China/XJ01/2019 and the most closely related strain, China/GD01/2020. The changed ORFs represent seven proteins, hypothetical protein (LSDV001), putative myristylated protein (LSDV059), superoxide dismutase-like protein (LSDV131), tyrosine protein kinase-like protein (LSDV143), kelch-like protein (LSDV144), putative ER-localized apoptosis regulator (LSDV154), and hypothetical protein (LSDV155). In total, seven ORFs were missing in the LSDV/China/XJ01/2019 genome compared to the China/GD01/2020 genome. Of these seven missing ORFs, all (LSDV004, LSDV023, LSDV044, LSDV055, LSDV057, LSDV106 and LSDV107) had both start and stop codons in the China/GD01/2020 strain, which may be annotation issues.

The LSDV/China/XJ01/2019 strain shared 153 identical ORFs with the Vietnamese strains (20L42_Quyet-Thang/VNM/20, 20L43_Ly-Quoc/VNM/20, 20L70_Dinh-To/VNM/20 and 20L81_Bang-Thanh/VNM/20) and 3 ORFs had 97.1–99.7% sequence identities. The changed ORFs included three proteins, putative myristylated protein (LSDV059), putative viral membrane protein (LSDV073) and tyrosine protein kinase-like protein (LSDV143). The LSDV/China/XJ01/2019 strain showed no additional ORFs compared to the Vietnamese strains.

A comparison of LSDV/China/XJ01/2019 with the Russian strains (LSDV/Russia/Dagestan/2015, LSDV/Russia/Saratov/2017 and LSDV/Russia/Udmurtiya/2019), identified 75 identical ORFs and 81 ORFs sharing 60.6–99.9% sequence identities. Two ORFs (LSDV114 and LSDV131) in LSDV/Russia/Dagestan/2015 and one ORF (LSDV081) in LSDV/Russia/Udmurtiya/2019 were missing in our genomic sequence.

A comparison of LSDV/China/XJ01/2019 with the Kazakhstani strains (Kubash/KAZ/16, Neethling-RIBSP, and KZ-Kostanay-2018), revealed 53 identical ORFs and 103 ORFs sharing 75.2–99.9% sequence identities. In total, 1 ORF (LSDV023) in Kubash/KAZ/16, 10 ORFs (LSDV001, LSDV002, LSDV003, LSDV004, LSDV023, LSDV027, LSDV135, LSDV148, LSDV155, and LSDV156) in Neethling-RIBSP and 10 ORFs (LSDV001, LSDV002, LSDV003, LSDV004, LSDV114, LSDV131, LSDV153, LSDV154, LSDV155, and LSDV156) in KZ-Kostanay-2018 were missing in the genomic sequence due to lacking of the complete 5′- and 3′-terminal sequences of Neethling-RIBSP and KZ-Kostanay-2018.

3.5. Phylogenetic analyses

LSDV/China/XJ01/2019 together with China/GD01/2020, LSDV/HongKong/2020, 20L42_Quyet-Thang/VNM/20, 20L43_Ly-Quoc/VNM/20, 20L70_Dinh-To/VNM/20, and 20L81_Bang-Thanh/ VNM/20 belonged to the same clade, showing their elevated similarity and comprised a monophyletic group with short tree branches. The phylogenetic tree confirmed that Chinese and Vietnamese strains belong to the same evolutionary lineage compared to other LSDV strains (Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3. Bayesian maximum clade credibility tree of whole-genome sequences of Lumpy skin disease virus. The whole-genome MCC tree was constructed using the BEAST software package (v. 2.4.3) and then visualized by FigTree (v1.4.3). The LSDV genome sequenced in this study are labeled in red.

4. Discussion

LSD is a severe systemic disease that was first reported in China 3 years ago, and it has severely impacted the cattle-breeding industry and brought significant economic losses to the affected areas (11). Before and after it, LSD outbreaks had occurred in several Eurasian countries, including Vietnam, Kazakhstan, Russia, India, and South Korea, and affected cattle, Korean water deer, and giraffe (4, 2936), indicating a wave of highly contagious epidemic. Here we present the genomic characterization of the virulent strain of LSDV, together with its homologous viruses that were subsequently discovered in China and Vietnam. This study provides additional genomic data for LSDV evolution and is crucial for virus tracking and vaccine development.

Genomic comparisons of the LSDV/China/XJ01/2019 strain with those in and around China showed that it showed the highest level of genomic similarity with the Chinese and Vietnamese strains. Phylogenetic analysis and genomic similarity comparisons using whole genome sequences indicated that the strain LSDV/China/XJ01/2019 was most closely related to the strains China/GD01/2020, LSDV/HongKong/2020, 20L42_Quyet-Thang/VNM/20, 20L43_Ly-Quoc/VNM/20, 20L70_Dinh-To/VNM/20 and 20L81_Bang-Thanh/VNM/20, inferred that these LSDV strains might originated from a common ancestor. Genome sequencing of more LSDV strains circulating in East and Southeast Asia may help pinpoint their origins.

A total of 17 variable loci in the genomic sequence caused variations in 7 ORFs and 2 non-coding regions between LSDV/China/XJ01/2019 and the most closely related strain, China/GD01/2020. The changed ORFs encoded seven proteins, which are hypothetical protein (LSDV001), putative myristylated protein (LSDV059), superoxide dismutase-like protein (LSDV131), tyrosine protein kinase-like protein (LSDV143), kelch-like protein (LSDV144), putative ER-localized apoptosis regulator (LSDV154) and hypothetical protein (LSDV155). According to the whole-genome comparison of LSDV, the putative myristylated protein (LSDV059) and tyrosine protein kinase-like protein (LSDV143) were highly variable, indicating that the diversity of these two proteins was probably resulted from adaptive evolutionary pressure. In future, it would be instructive to reveal the association between gene mutations and viral pathogenicity or transmissibility, which would be helpful in understanding LSDV adaptive evolution.

Data on the movements of cattle into the epidemic zone during the first LSD outbreak in China and on LSDV-infected cattle in region bordering China and Kazakhstan are still lacking, indicating the emergence of LSDV in Xinjiang, northwest China may be unexpectedly complex and difficult to trace, as it is still unknown whether the outbreak was imported or localized, highlighting the need for further research. Moreover, the phylogeny of the LSDV strains in Xinjiang, northwest China, and Vietnam, which are thousands of kilometers apart, highlights how little we know about the spread of this lineage and its introduction to these regions.

In summary, we analyzed the LSDV genome from a dairy cow that died during the first LSD outbreak in China and found that the strain LSDV/China/XJ01/2019 strain was genetically close to additional LSDV strains in China and Vietnam. Molecular epidemiological investigations of LSDV in susceptible animals and vectors at national, regional, and global levels are desired to understand the evolution and transmission routes of the latest global LSD epidemic.

Data availability statement

The data presented in the study is deposited in the Genbank, accession number OM105589.

Ethics statement

Ethical review and approval was not required for the animal study because specimens were collected from a dead cow, but it was approved by the Animal Husbandry and Veterinary Bureau of Xinjiang Uygur Autonomous Region. Written informed consent was obtained from the owners for the participation of their animals in this study.

Author contributions

JH conceptualized the study. Y-RW, W-GM, PW, WW, X-HS, X-YY, X-YM, and J-YW carried out the laboratory assays. Y-RW analyzed the data and drafted the paper. All authors approved the final version of the manuscript to be published.

Funding

This work was supported by the National Science and Technology Major Project (Grant No. 2021YFD1800505) and the Science and Technology Major Project of the Xinjiang Uygur Autonomous Region of China (2020A01001).

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.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.1073648/full#supplementary-material

References

1. Das M, Chowdhury M, Akter S, Mondal A, Uddin M, Rahman M, et al. An updated review on lumpy skin disease: a perspective of Southeast Asian countries. J Adv Biotechnol Exp Therap. (2021) 4:322. doi: 10.5455/jabet.2021.d133

CrossRef Full Text | Google Scholar

2. Molini U, Boshoff E, Niel AP, Phillips J, Khaiseb S, Settypalli TBK, et al. Detection of lumpy skin disease virus in an asymptomatic eland (Taurotragus oryx) in Namibia. J Wildl Dis. (2021) 57:708–11. doi: 10.7589/JWD-D-20-00181

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Ko YS, Oh Y, Lee TG, Bae DY, Tark D, Cho HS. Serological and molecular prevalence of lumpy skin disease virus in Korean water deer, native and dairy cattle in Korea. Korean J Vet Serv. (2022) 45:133–7. doi: 10.7853/kjvs.2022.45.2.133

CrossRef Full Text | Google Scholar

4. Tung Duy D, Long Hoang T, Hiep Dinh N, Thuy Thi H, Giang Hoang N, Kien Viet Dung T, et al. Characterization of Lumpy skin disease virus isolated from a C in Vietnam. Transbound Emerg Dis. (2022) 69:E3268–72. doi: 10.1111/tbed.14583

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Vandenbussche F, Mathijs E, Haegeman A, Al-Majali A, Van Borm S, De Clercq K. Complete genome sequence of capripoxvirus strain KSGP 0240 from a commercial live attenuated vaccine. Genome Announc. (2016) 4:e01114–01116. doi: 10.1128/genomeA.01114-16

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Sprygin A, Artyuchova E, Babin Y, Prutnikov P, Kostrova E, Byadovskaya O, et al. Epidemiological characterization of lumpy skin disease outbreaks in Russia in 2016. Transbound Emerg Dis. (2018) 65:1514–21. doi: 10.1111/tbed.12889

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Tuppurainen ESM, Venter EH, Coetzer JW, Bell-Sakyi L. Lumpy skin disease: Attempted propagation in tick cell lines and presence of viral DNA in field ticks collected from naturally-infected cattle. Ticks Tick-Borne Dis. (2015) 6:134–40. doi: 10.1016/j.ttbdis.2014.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Sanz-Bernardo B, Haga IR, Wijesiriwardana N, Basu S, Larner W, Diaz AV, et al. Quantifying and modeling the acquisition and retention of lumpy skin disease virus by hematophagus insects reveals clinically but not subclinically affected cattle are promoters of viral transmission and key targets for control of disease outbreaks. J Virol. (2021) 95:e02239–e02220. doi: 10.1128/JVI.02239-20

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Sprygin A, Pestova Y, Wallace DB, Tuppurainen E, Kononov AV. Transmission of lumpy skin disease virus: a short review. Virus Res. (2019) 269:197637. doi: 10.1016/j.virusres.2019.05.015

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Gupta T, Patial V, Bali D, Angaria S, Sharma M, Chahota R. A review: lumpy skin disease and its emergence in India. Vet Res Commun. (2020) 44:111–8. doi: 10.1007/s11259-020-09780-1

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Liu P, Li J, Chen R, Cheng Z, Shi Y, Baosuyalatu, et al. The first outbreak investigation of lumpy skin disease in China. China Anim Health Inspection. (2020) 37:1–5. doi: 10.3969/j.issn.1005-944X.2020.01.001 (in Chinese)

CrossRef Full Text | Google Scholar

12. Wang Y, Zhao L, Yang J, Shi M, Nie F, Liu S, et al. Analysis of vaccine-like lumpy skin disease virus from flies near the western border of China. Transbound Emerg Dis. (2022) 69:1813–23. doi: 10.1111/tbed.14159

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Morgenstern M, Klement E. The effect of vaccination with live attenuated neethling lumpy skin disease vaccine on milk production and mortality-an analysis of 77 dairy farms in Israel. Vaccines. (2020) 8:324. doi: 10.3390/vaccines8020324

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Salnikov N, Usadov T, Kolcov A, Zhivoderov S, Morgunov Y, Gerasimov V, et al. Identification and characterization of lumpy skin disease virus isolated from cattle in the Republic of North Ossetia-Alania in 2015. Transbound Emerg Dis. (2018) 65:916–20. doi: 10.1111/tbed.12818

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Mathijs E, Vandenbussche F, Saduakassova M, Kabduldanov T, Haegeman A, Aerts L, et al. Complete coding sequence of a lumpy skin disease virus strain isolated during the 2016 outbreak in Kazakhstan. Microbiol Resour Announc. (2020) 9:e01399–19. doi: 10.1128/MRA.01399-19

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Lu G, Xie J, Luo J, Shao R, Jia K, Li S. Lumpy skin disease outbreaks in China, since August 3, 2019. Transbound Emerg Dis. (2021) 68:216–9. doi: 10.1111/tbed.13898

PubMed Abstract | CrossRef Full Text | Google Scholar

17. 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:e336–e343. doi: 10.1111/tbed.14304

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Wu J, Chen Y, Che C, Li B. Epidemiological investigation and analysis of an epidemic of Lumpy skin disease in Sichuan, China. China Anim Health. (2021) 21:20–1. doi: 10.3969/j.issn.1008-4754.2021.01.014 (in Chinese)

CrossRef Full Text | Google Scholar

19. Yang Z, Hou W, Mo Q, Yin J, Yuan D, Zhou R, et al. International and domestic prevalence of Lumpy skin disease. Shandong J Anim Sci Vet Med. (2021) 42:66–71. doi: 10.3969/j.issn.1007-1733.2021.10.020 (in Chinese)

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Tulman ER, Afonso CL, Lu Z, Zsak L, Kutish GF, Rock DL. Genome of lumpy skin disease virus. J Virol. (2001) 75:7122–30. doi: 10.1128/JVI.75.15.7122-7130.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

21. El-Tholoth M, El-Kenawy AA. G-Protein-Coupled Chemokine Receptor Gene in Lumpy Skin Disease Virus Isolates from Cattle and Water Buffalo (Bubalus bubalis) in Egypt. Transbound Emerg Dis. (2016) 63:e288–95. doi: 10.1111/tbed.12344

PubMed Abstract | CrossRef Full Text | Google Scholar

22. 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:1142. doi: 10.3390/microorganisms9061142

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Ma J, Yuan Y, Shao J, Sun M, He W, Chen J, et al. Genomic characterization of lumpy skin disease virus isolated in southeast China. Transbound Emerg Dis. (2022) 69:2788–99. doi: 10.1111/tbed.14432

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Koichiro T, Glen S, Sudhir K. MEGA11: molecular evolutionary genetics analysis version 11. Mol Biol Evol. (2021) 38:3022–7. doi: 10.1093/molbev/msab120

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. (2012) 28:1647–9. doi: 10.1093/bioinformatics/bts199

PubMed Abstract | CrossRef Full Text | Google Scholar

26. 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:150. doi: 10.1186/1471-2164-7-150

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Okonechnikov K, Golosova O, Fursov M, Team U. Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics. (2012) 28:1166–7. doi: 10.1093/bioinformatics/bts091

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Ronquist F, Teslenko M, Van Der Mark P, Ayres DL, Darling A, Hohna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. (2012) 61:539–42. doi: 10.1093/sysbio/sys029

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Acharya KP, Subedi D. First outbreak of lumpy skin disease in Nepal. Transbound Emerg Dis. (2020) 67:2280–1. doi: 10.1111/tbed.13815

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Kumar N, Chander Y, Kumar R, Khandelwal N, Riyesh T, Chaudhary K, et al. Isolation and characterization of lumpy skin disease virus from cattle in India. PLoS ONE. (2021) 16:e0241022. doi: 10.1371/journal.pone.0241022

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Orynbayev MB, Nissanova RK, Khairullin BM, Issimov A, Zakarya KD, Sultankulova KT, et al. Lumpy skin disease in Kazakhstan. Trop Anim Health Prod. (2021) 53:166. doi: 10.1007/s11250-021-02613-6

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Tran HTT, Truong AD, Dang AK, Ly DV, Nguyen CT, Chu NT, et al. Lumpy skin disease outbreaks in vietnam, 2020. Transbound Emerg Dis. (2021) 68:977–80. doi: 10.1111/tbed.14022

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Byadovskaya O, Prutnikov P, Shalina K, Babiuk S, Perevozchikova N, Korennoy F, et al. The changing epidemiology of lumpy skin disease in Russia since the first introduction from 2015 to 2020. Transbound Emerg Dis. (2022) 69:E2551–62. doi: 10.1111/tbed.14599

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Chouhan CS, Parvin S, Ali MY, Sadekuzzaman M, Chowdhury MGA, Ehsan MA, et al. Epidemiology and economic impact of lumpy skin disease of cattle in Mymensingh and Gaibandha districts of Bangladesh. Transbound Emerg Dis. (2022). doi: 10.1111/tbed.14697

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Odonchimeg M, Erdenechimeg D, Tuvshinbayar A, Tsogtgerel M, Bazarragchaa E, Ulaankhuu A, et al. Molecular identification and risk factor analysis of the first Lumpy skin disease outbreak in cattle in Mongolia. J Vet Med Sci. (2022) 84:1244–52. doi: 10.1292/jvms.22-0250

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Sariya L, Paungpin W, Chaiwattanarungruengpaisan S, Thongdee M, Nakthong C, Jitwongwai A, et al. Molecular detection and characterization of lumpy skin disease viruses from outbreaks in Thailand in 2021. Transbound Emerg Dis. (2022) 69:E2145–52. doi: 10.1111/tbed.14552

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: Lumpy skin disease virus, Illumina sequencing, phylogenetic analysis, genome sequence, genome alignment

Citation: Wei Y-R, Ma W-G, Wang P, Wang W, Su X-H, Yang X-Y, Mi X-Y, Wu J-Y and Huang J (2023) Retrospective genomic analysis of the first Lumpy skin disease virus outbreak in China (2019). Front. Vet. Sci. 9:1073648. doi: 10.3389/fvets.2022.1073648

Received: 18 October 2022; Accepted: 13 December 2022;
Published: 12 January 2023.

Edited by:

Moh A. Alkhamis, Kuwait University, Kuwait

Reviewed by:

Mohamed Nayel, University of Sadat City, Egypt
Laxmi Narayan Sarangi, National Dairy Development Board, India

Copyright © 2023 Wei, Ma, Wang, Wang, Su, Yang, Mi, Wu and Huang. 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: Jiong Huang, yes jh124@163.com

Disclaimer: 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.