Skip to main content

Molecular characterization and phylogenetic analysis of bovine ephemeral fever viruses in Khuzestan province of Iran in 2018 and 2020

Abstract

Background

Bovine ephemeral fever (BEF) is an arthropod-borne viral disease caused by the BEF virus (BEFV). This single-stranded RNA virus that affects cattle and water buffalo is endemic in tropical and subtropical regions including Iran. While BEF is a major disease of cattle in Iran, information regarding its agent, molecular characterization, and circulating viruses are highly limited. The current study aimed to, firstly, determine the genetic and antigenic characteristics of BEFV strains in Khuzestan province in Southwest of Iran in 2018 and 2020 and, secondly, to compare them with strains obtained from other areas.

Results

By phylogenetic analysis based on the Glycoprotein gene, BEFV strains were divided into four clusters of Middle East, East Asia, South Africa, and Australia; in which the 2018 and 2020 Iranian BEFV strains were grouped in the Middle East cluster with the Turkish, Indian, and Israeli strains. Depending on the chronology and geographical area, the outbreaks of Turkey (2020), Iran (2018 and 2020), and India (2018 and 2019) are proposed to be related. These BEFVs had the highest identity matrix and the lowest evolutionary distance among the studied strains. Multiple sequence alignment of G1, G2, and G3 antigenic sites showed that these neutralizing epitopes are highly conserved among the strains of the Middle East cluster; however, the strains previously identified in Iran differed in three amino acids placed in G1 and G2 epitopes.

Conclusion

The findings revealed that BEFVs circulating in the Middle East are closely related phylogenetically and geographically. They also have similar antigenic structures; therefore, developing a vaccine based on these strains can be effective for controlling BEF in the Middle East.

Peer Review reports

Introduction

Bovine ephemeral fever (BEF), also called three-day sickness, is a viral disease that affects cattle and water buffalo. BEF virus (BEFV) is endemic in tropical and subtropical regions of Australia, Africa, and Asia, especially the Middle East countries [1, 2]. The disease is marked by a bi-phasic or multiple-phasic fever, anorexia, muscle stiffness, salivation, depression, lameness, ocular and nasal discharge, cessation of rumination, and constipation [3]. BEF causes considerable economic losses, mainly because of its consequences and imposing international trade restrictions [4, 5]. BEF, caused by an arthropod-borne virus, usually occurs in the summer and early autumn months [3]. Weather conditions plays a unique role in the spread of the BEF. Rainfall, prevailing wind patterns, streams, or other ground water sources have a significant effect on its occurrence, particularly in wet and dry seasons [6].

BEFV belongs to the genus Ephemerovirus, Rhabdoviridae family and has a negative-sense single-stranded RNA (ssRNA) [7]. The length of the viral genome is 14,900 nucleotides and encodes five structural proteins of the matrix, nucleoprotein, phosphoprotein, surface glycoprotein, and RNA-dependent RNA polymerase. Surface glycoprotein (G) is a class I transmembrane glycoprotein that includes the major neutralizing epitopes and can induce protective immunity in cattle [8, 9]. Four distinct antigenic sites (G1-G4) have been determined on the surface of this protein [9,10,11]. The viral capsid binds to the cell surface via G protein and penetrates into the cells through receptor-mediated endocytosis. Therefore, G protein plays an important role in the virus replication process as well as induction of immune response in the host body [6].

BEF disease is prevalent in Middle East countries and is considered an enzootic disease in some areas of this region [4, 12]. The first detailed report of BEF in the Middle East was described in Egypt in 1909, in which the disease spread from Aswan, Nile Valley to Cairo, and spread across the Delta to the coast [11]. Afterward, the disease was reported in the Palestine and Jordan Valley in 1931 [12, 13], followed by Israel, Turkey, Iraq, Saudi Arabia, Kuwait, Yemen, Syria, and Iran [12, 14,15,16]. In Iran, the first report of BEFV isolation date back to the occurrence of the virus in the southern and eastern regions of the country in 1974 [14]. The disease is mostly sporadic and is prevalent in southern and western areas with tropical climates [14]. Three outbreaks are also recorded in 2006 [17] and in autumn of 2012 and 2013 [14] in Iran, before the last outbreak seen in 2020. The disease was appeared in Khuzestan in March as sporadic cases and then gradually appeared in the other provinces and soon reached to at least fourteen provinces (unpublished data). The peak of the disease was observed in July and August. While no precise estimate is available about losses caused by the outbreak; 70% of livestock in the infected areas showed clinical signs of disease (unpublished data).

Although BEF is known to occur as a single serotype worldwide, some antigenic variations are detected at virus epitope sites [18]. Detection of these variations would be useful for identifying the epidemiology of BEFV, genetic evolution, and design broad-spectrum vaccines [18]. To determine the genetic and antigenic variations of BEFV strains in Iran in 2018 and 2020, analyzed the nucleotide and deduced amino acid sequences of complete G gene of BEFV isolates and compared the sequences with those of the GenBank database.

Results

Detection of BEFV

RT-nested PCR showed that all suspicious samples were positive for BEFV and a specific product of 420 bp was found. To determine the genetic characteristics of BEFVs circulating in Iran during 2018 and 2020, the glycoprotein encoding gene of 8 BEFV strains (4 strains related to 2018 and 4 strains related to 2020) were sequenced. Because the G gene sequences of the prevalent strains in each year were the same, only the sequence of one strain in each year was submitted to the GenBank. The accession numbers of MZ51169 and MZ51168 were acquired for 2018 and 2020 Iranian BEFV strains, respectively.

Identity and evolutionary distance analysis

G gene nucleotide and deduced amino acid sequences of Iranian BEFV strains were compared with other BEFV sequences in the GenBank using multiple sequence alignment. The IR-2018 and IR-2020 strains showed the lowest identity with South African strains (86.8–87.1% nucleotide identity and 93–93.5% amino-acid identity). The IR-2018 strain shared the most nucleotide identity with the IR-2020 strain (98.4%), and Indian strains (98.2%), respectively, while concerning amino acid sequence, it showed the highest identity (99%) with the 2020 Turkish, 2006 Israeli and IR-2020 strains. The G gene sequence of the IR-2020 strain showed the highest nucleotide identity with strains identified in Turkey in 2020 (99.4–99.6% identity) as well as Indian strains (99.1% identity). Two BEFV strains of 2020 Iran and Turkey were quite similar in terms of amino acid sequence and were more than 99% similar to the Indian strains. Two investigated Iranian strains had 98.4 and 99% nucleotide and amino acid sequence identity, respectively (Table 1 and Supplementary Table 1).

Table 1 Distance estimation (left, bottom) and Identity Matrix (right, top) of BEFVs from different countries according to G gene

The evolutionary distance analysis of Iranian strains showed similar results with identity. These strains had the highest and lowest evolutionary distances with the South African and the Middle East strains, respectively (Table 1).

Phylogenetic analysis

The phylogenetic tree was constructed with the G ectodomain encoding sequences of the viruses. As shown in Fig. 1, the BEFV strains formed four major clusters of East Asian, Australian, South Africa, and Middle East strains. The IR-2018 and IR-2020 BEFV strains were grouped in the Middle East cluster with the Turkish, Indian, and Israeli strains. The viruses related to 2020 outbreaks in Turkey and Iran were placed on the same branch, while BEFV strains of Indian-2018 and 2019, Turkish-2020, and Iranian-2018 and 2020 formed a sub-cluster.

Fig. 1
figure 1

The phylogenetic tree of G gene ectodomain sequences constructed by the Maximum-likelihood using the MegaX program. BEFV strains were grouped into four clusters of I: East Asia, II: Middle East, III: Australia, and IV: South Africa. Sequencing viruses in this study are marked with a black square and previously studied viruses in Iran are marked with a black circle

Analysis of antigenic sites

No codon insertion, deletion, or frame shift was found in the glycoprotein sequences of IR-2018 and IR-2020 BEFV strains. All changes have resulted from single nucleotide substitutions. Multiple sequence alignments of G1, G2, and G3 antigenic sites showed that these neutralizing epitopes are highly conserved among the strains of the Middle East cluster, while it was found that the three substitutions at positions 218 (K to R), 223 (E to P), and 503 (K to T in IR-2013 strain and E to T in IR-2012 strain) in Iranian BEFV strains previously identified compared to present strains. Positions 218 and 503 are putative N-linked glycosylation sites. If the entire sequence of the G gene is considered, there are at least 22 amino acid substitutions between the strains previously identified in Iran and those present (Fig. 2). The Japanese YHL strain differed at 15 positions from IR-2018 and IR-2020 strains in which two substitutions of aa 170 (N to T) and aa 503 (K to T) were located in G2 and G1 antigenic sites, respectively.

Fig. 2
figure 2

Alignment of deduced amino acid sequences of the neutralizing epitopes G1 (487–503), G2 (168–189), G3a (49–63), G3b (215–231) and G3c (262–271) on the G protein of different bovine ephemeral fever virus strains. Accession number, country and year of each strains were indicated on the left side, respectively. Putative glycosylation sites are boxed

Discussion

Although BEF is one of the most important diseases of cattle in Iran [14], little is available about its agent and characteristics of circulating viruses. Therefore, in the present study, we characterized the prevalent BEFV strains in Khuzestan province southwest of Iran in 2018 and 2020.

In 2018, the incidence of BEF was sporadic and no report of the disease is available in other regions of the country. Our surveillance program of the BEF showed the absence of the disease in 2019. However, many animals suspected of contracting the three-day sickness were reported in Khuzestan province and other parts of the country in 2020. The RT-nested PCR results confirmed the occurrence of the disease in Khuzestan province as well as other regions of the country. The BEF disease epizootic has been reported in Turkey simultaneous with an outbreak in Iran [19]. The phylogenetic analysis, identity matrix, and distance evolution confirmed the high genetic closeness of prevalent strains in Turkey and Iran in 2020.

According to the phylogenetic analysis of the G gene ectodomain region, the BEFV strains were categorized into four clusters of East Asia, Australia, Middle East, and South Africa. However, previous studies have grouped the BEF strains into the three clusters of East Asia, Middle East, and Australia [20,21,22,23,24]. According to our results and the study of Omar et al [25], the BEFV strains of South Africa are distinct from those of other regions of the world and fall into the new cluster of South Africa. These strains showed the most distance evolution and the lowest similarity with the identified strains in the present study.

The phylogenetic tree revealed that the IR-2018 and IR-2020 strains are in the Middle East cluster close to Turkish, Indian, and Israeli BEFV strains. However, the viruses related to the outbreaks of 2012 and 2013 in Iran were grouped in East Asia cluster close to the strains from Taiwan, isolated in 2013, and Turkey, isolated in 2012. Most of the strains identified in Turkey belong to the Middle East group; however, in 2012, in some regions of Turkey, some identified BEF viruses were more similar to those placed in the East Asian group [16, 21, 26]. The BEF viruses that were identified in different regions of Iran during 2012 and 2013 were highly similar to those from Turkey [14]. Phylogenetic, identity matrix, and distance evolution analyzes revealed that the BEF viruses previously identified in Iran have a genetic distance from the strains of the present study and may have different sources. It seems that the virus circulating in 2012 and 2013 probably has expanded to the Middle East from East Asia [16, 26], while circulating strains of 2018 and 2020 were domestic strains of the Middle East. Studies of G gene nucleotide and deduced amino acid sequences from 2000 to 2020 indicate the circulation of a domestic strain in the Middle East. However, in 2012 and 2013, another strain that seemed to be exotic was prevalent in Turkey and Iran, but no evidence of the virus has been available since. It remains to be determined whether the structure of the G protein and its epitopes are effective in the compatibility of the virus with competence vectors [27].

This study demonstrated that BEFVs circulating in the Middle East are closely related phylogenetically and geographically. Pyasi et al. reported the outbreaks of BEFV during 2018 and 2019 in India [28]. Their molecular studies have shown that these strains were evolutionarily very close to those of Turkey and Israel, so they strongly suggested that these strains reached India from Israel and Turkey while maintaining their genomic sequence [28]. The findings of the present study also confirm this view. These viruses probably traversed across the countries like Turkey and Israel reaching Iran, and then India. Iran and India do not have a common border, so it is unclear exactly how the virus has been transmitted from Iran to India. The similarity between BEFVs from Iran’s and Turkey’s 2020-outbreaks (99.4–99.6% nucleotide similarity) was higher than that of IR-2018 and IR-2020 BEFVs (98.4% nucleotide similarity). The simultaneous outbreaks of BEFV in Iran and Turkey and high nucleotide and amino acid similarity indicate the common source of these viruses.

Comparison of G1, G2, and G3 antigenic sites showed that these neutralizing epitopes are highly conserved among the Middle East strains; however, the strains previously identified in Iran differed in three amino acids placed in G1 and G2 epitopes, which two of them were putative attachment site of oligosaccharides. Japanese YHL strain (vaccine strain used in Iran) differed from the Middle East BEFVs, especially IR-2018 and IR-2020 strains in G1 and G2 epitopes. The amino acid substituted in the G1 epitope is a putative N-linked glycosylation site; therefore, it may affect conformation and recognition of the neutralizing epitopes. Consistent with the differences in antigenic sites, Iranian and YHL BEFVs were placed as two distinct clusters. Due to these differences in antigenic structure, the use of heterologous strains as a vaccine may not provide full protective immunity [29]. According to the similarity of circulating strains in the Middle East, the use of these strains to develop vaccines can play an important role in preventing this disease.

Conclusion

Comparison of the nucleotide and deduced amino acid sequences of the BEFV G gene related to different countries showed that these viruses are geographically divided into four clusters: Middle East, East Asia, South Africa, and Australia. Based on the chronology and geographical area, the outbreaks of Turkey (2020), Iran (2018 and 2020), and India (2018 and 2019) are proposed to be related, which also suggests that these strains can be used to develop a BEFV vaccine.

Materials and methods

Ethics approval

Written consent was obtained from all farm owners before entering the study. Sampling was performed by a specialist veterinarian to minimize pain and injury to the animal. In addition, the current study is approved by the Ethics Committee of Shahid Chamran University of Ahvaz. All methods were performed in accordance with the relevant guidelines.

Collecting samples

This study is conducted in Khuzestan province, southwest of Iran, with a tropical climate. Heparinized whole blood samples were collected from cattle with clinical symptoms existing in affected farms. The most important clinical symptoms in these animals were persistent fever, muscle stiffness, constipation, ocular and nasal discharge, and spontaneous recovery after 3 days. A total of 50 and 40 samples were taken from suspicious animals in 2018 and 2020, respectively. Blood samples were centrifuged at 3000 rpm for 15 min to separate buffy coats. The buffy coat layer was transferred to a new tube and then rinsed three times.

RNA extraction and RT-nested PCR

RNA was extracted from buffy coats using a viral RNA extraction kit (Bioneer; South Korea) according to the manufacturer’s protocol. Then, extracted RNA samples were reverse-transcribed using a kit (Yekta Tajhiz; Iran) to obtain cDNA. Screening of the positive BEFV samples was performed based on the amplification of the partial G gene via a nested PCR [30, 31]. Primer sequences are indicated in Table 2.

Table 2 List of primes used for detection and sequencing of glycoprotein (G) gene

Sequencing the full-length G gene

BEFV G gene was amplified as described by Hsieh et al with some modification [32]. A nested PCR was designed according to primers introduced by Hsieh et al. In the first run, G1F and G4R primers were used to amplify an 1872 length fragment with the following protocol: 94 ̊C for 2 min, 25 cycles of 94 ̊C for 50 s, 50 ̊C for 50 s, 72 ̊C for 75 s, and a final extension of 72 ̊C for 5 min. In the second run, the PCR product was used as a template and G1-G4 fragments were amplified. PCR products were directly sequenced or subcloned into the pTG19-T vector (Vivantis; Malaysia) using standard techniques and then sequenced again. Sequencing was performed by Bioneer Company (South Korea) with the same PCR primers in two directions. Sequences were trimmed with BioEdit software version 7.0.4.1 (mbio, Inc., North Carolina, USA). Obtained sequences were submitted to the GenBank and are available under the accession numbers MZ51169 and MZ51168.

Phylogenetic analysis

Available global nucleotide sequences of the BEFV G gene were retrieved from the NCBI Genbank database. All retrieved sequences from different countries and those of Iran generated in this study were aligned with ClustalW [33]. A phylogenetic tree was constructed based on G ectodomain encoding sequences using Molecular Evolutionary Genetics Analysis the MegaX software package (10.0.4) [34]. Tree construction was performed using Maximum-likelihood (ML) method, General Time Reversible (GTR) Model as the nucleotide substitution model [35], and gamma-distributed 4 (G4) according to a recent report [27]. Sub-tree Pruning & Re-grafting (SPR) branch swapping were applied to the analysis. The reliability of the branching was evaluated by the bootstrap analysis with 1000 replicates.

Distance estimation and sequence identity matrix

The pairwise distance was estimated using MegaX software. The maximum composite likelihood was used as a substitution model. A discrete Gamma distribution (G) was used to model evolutionary rate differences among sites. The number of bootstrap replication was considered 1000. The sequence identity matrix was calculated using BioEdit software.

Amino acid sequence analysis

Indicated nucleotide sequences of G encoding genes were translated to amino acid sequences using Expasy found at the expasy.org web server. Then, antigenic sites, including G1 (487aa–503aa), G2 (168aa–189aa), G3a (49–63), G3b (215aa–231aa), and G3c (262aa–271aa) [36] were retrieved and aligned with BioEdit software. These sequences were also compared with the Japanese Yamaguchi (YHL) vaccine strain, which is also used in Iran. In addition, the amino acid changes in putative N-linked glycosylation regions were also determined [18].

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. Sequence data of this project have been deposited in the GenBank of the National Center for Biotechnology Information (NCBI) under the accession number MZ51169 and MZ51168.

Abbreviations

BEF:

Bovine ephemeral fever

BEFV:

BEF virus

ssRNA:

single-stranded RNA

ML:

Maximum-likelihood

GTR:

General Time Reversible

G4:

Gamma-distributed 4

SPR:

Sub-tree Pruning & Re-grafting

References

  1. Mlingo TA, Nthangeni BM, Mokoena NB. Genome sequence of bovine ephemeral fever virus vaccine strain of South African origin. Vet Med Sci. 2021;7(5):1611–15.

  2. Jiang H, Hou P, He H, Wang H. Cell apoptosis regulated by interaction between viral gene alpha 3 and host heterogeneous nuclear ribonucleoprotein K facilitates bovine ephemeral fever virus replication. Vet Microbiol. 2020;240:108510.

    Article  CAS  Google Scholar 

  3. Stokes JE, Darpel KE, Gubbins S, Carpenter S, de Marco MMF, Hernández-Triana LM, et al. Investigation of bovine ephemeral fever virus transmission by putative dipteran vectors under experimental conditions. Parasit Vectors. 2020;13(1):1–11.

    Article  Google Scholar 

  4. Mirzaie K, Bahonar A, Mehrabadi MF, Hajilu G, Yaghoubi M. Determinants of bovine ephemeral fever outbreak during 2013, in Qazvin Province. Iran Asian Pac J Trop Dis. 2017;7:744–7.

    Article  Google Scholar 

  5. Hayama Y, Moriguchi S, Yanase T, Suzuki M, Niwa T, Ikemiyagi K, et al. Epidemiological analysis of bovine ephemeral fever in 2012–2013 in the subtropical islands of Japan. BMC Vet Res. 2016;12(1):1–13.

    Article  Google Scholar 

  6. Walker P. Bovine ephemeral fever in Australia and the world. Curr Top Microbiol Immunol. 2005;292:57–80.

    CAS  PubMed  Google Scholar 

  7. Kun J, Rongrong J, Xiangbin W, Yan Z, Yiping D, Gang L, et al. Genetic characterization of bovine ephemeral fever virus in southern China, 2013–2017. Virus Genes. 2020;56(3):390–5.

    Article  Google Scholar 

  8. Zaghawa AA, Housawi F, Al-Naeem A, Elsify A, Hegazy YM. Bovine ephemeral fever epidemics in kingdom Saudi Arabia: clinical, epidemiological and molecular investigation. J Infect Dev Ctries. 2017;11(11):854–60.

    Article  CAS  Google Scholar 

  9. Ting LJ, Lee MS, Lee SH, Tsai HJ, Lee F. Relationships of bovine ephemeral fever epizootics to population immunity and virus variation. Vet Microbiol. 2014;173(3–4):241–8.

    Article  CAS  Google Scholar 

  10. Johal J, Gresty K, Kongsuwan K, Walker PJ. Antigenic characterization of bovine ephemeral fever rhabdovirus G and GNS glycoproteins expressed from recombinant baculoviruses. Arch Virol. 2008;153(9):1657–65.

    Article  CAS  Google Scholar 

  11. Walker PJ, Klement E. Epidemiology and control of bovine ephemeral fever. Vet Res. 2015;46:124.

    Article  Google Scholar 

  12. Yeruham I, Van Ham M, Stram Y, Friedgut O, Yadin H, Mumcuoglu KY, et al. Epidemiological investigation of bovine ephemeral fever outbreaks in Israel. In: Veterinary medicine international 2010; 2010.

    Google Scholar 

  13. Lee F. Bovine ephemeral fever in Asia: recent status and research gaps. Viruses. 2019;11(5):412.

    Article  Google Scholar 

  14. Bakhshesh M, Abdollahi D. Bovine ephemeral fever in Iran: diagnosis, isolation and molecular characterization. J Arthropod Borne Dis. 2015;9(2):195–203.

    PubMed  PubMed Central  Google Scholar 

  15. Abu Elzein EM, Gameel AA, al Afaleq AI, al Gundi O, al Bashier AM, Zeedan A, et al. Observations on the recent epizootic of bovine ephemeral fever in Saudi Arabia. Rev Sci Tech. 1999;18(3):672–80.

    Article  CAS  Google Scholar 

  16. Alkan F, Albayrak H, Timurkan MO, Ozan E, Coskun N. Assessment of the molecular epidemiology of bovine ephemeral fever in Turkey. Veterinarski arhiv. 2017;87(6):665–75.

    Article  Google Scholar 

  17. Bazargani TT, Movassaghi AR, Bahonar AR, Hassan EB, Hemmatzadeh F, Khedmati K. Identification of bovine ephemeral fever (BEF) outbreak in a dairy farm in Varamin. Iran Iranian J Vet Sci Technol. 2013;5(1):35–44.

    Google Scholar 

  18. Kongsuwan K, Cybinski DH, Cooper J, Walker PJ. Location of neutralizing epitopes on the G protein of bovine ephemeral fever rhabdovirus. J Gen Virol. 1998;79(11):2573–81.

    Article  CAS  Google Scholar 

  19. Karayel-Hacioglu I, Yelken SD, Vezir Y, Unal N, Alkan F. Isolation and genetic characterization of bovine ephemeral fever virus from epidemic-2020 in Turkey. Trop Anim Health Prod. 2021;53(2):1–7.

    Article  Google Scholar 

  20. Abayli H, Tonbak S, Azkur AK, Bulut H. Complete genome analysis of highly pathogenic bovine ephemeral fever virus isolated in Turkey in 2012. Arch Virol. 2017;162(10):3233–8.

    Article  CAS  Google Scholar 

  21. Tonbak S, Berber E, Yoruk MD, Azkur AK, Pestil Z, Bulut H. A large-scale outbreak of bovine ephemeral fever in Turkey, 2012. J Vet Med Sci. 2013;75(11):1511–4.

    Article  Google Scholar 

  22. Kato T, Aizawa M, Takayoshi K, Kokuba T, Yanase T, Shirafuji H, et al. Phylogenetic relationships of the G gene sequence of bovine ephemeral fever virus isolated in Japan. Taiwan Aust Vet Microbiol. 2009;137(3–4):217–23.

    Article  CAS  Google Scholar 

  23. Lin G, Qiu C. Phylogenetic relationships of the partial G gene sequence of bovine ephemeral fever virus isolated from mainland China, Taiwan, Japan, Australia, Turkey and Israel. J Anim Vet Adv. 2012;11(17):3217–22.

    Article  Google Scholar 

  24. Zheng F, Qiu C. Phylogenetic relationships of the glycoprotein gene of bovine ephemeral fever virus isolated from mainland China, Taiwan, Japan, Turkey. Isr Aust Virol J. 2012;9(1):1–8.

    Google Scholar 

  25. Omar R, Van Schalkwyk A, Carulei O, Heath L, Douglass N, Williamson A-L. South African bovine ephemeral fever virus glycoprotein sequences are phylogenetically distinct from those from the rest of the world. Arch Virol. 2020;165(5):1207–10.

    Article  CAS  Google Scholar 

  26. Oğuzoğlu T, Ertürk A, Çizmeci Ş, Koç B, Akça Y. A report on bovine ephemeral fever virus in Turkey: antigenic variations of different strains of EFV in the 1985 and 2012 outbreaks using partial glycoprotein gene sequences. Transbound Emerg Dis. 2015;62(5):e66–70.

    Article  Google Scholar 

  27. Trinidad L, Blasdell KR, Joubert DA, Davis SS, Melville L, Kirkland PD, et al. Evolution of bovine ephemeral fever virus in the Australian episystem. J Virol. 2014;88(3):1525–35.

    Article  Google Scholar 

  28. Pyasi S, Sahu BP, Sahoo P, Dubey PK, Sahoo N, Byrareddy SN, et al. Identification and phylogenetic characterization of bovine ephemeral fever virus (BEFV) of middle eastern lineage associated with 2018–2019 outbreaks in India. Transbound Emerg Dis. 2020;67(5):2226–32.

    CAS  Google Scholar 

  29. Almasi S, Bakhshesh M. Antigenic variation of bovine ephemeral fever viruses isolated in Iran, 2012–2013. Virus Genes. 2019;55(5):654–9.

    Article  CAS  Google Scholar 

  30. Wang FI, Hsu AM, Huang KJ. Bovine ephemeral fever in Taiwan. J Vet Diagn Investig. 2001;13(6):462–7.

    Article  CAS  Google Scholar 

  31. Zheng F-y, Lin G-z, Qiu C-q, Yuan K-z, Song J-y. Expression and antigenic characterization of the epitope-G1 of the bovine ephemeral fever virus glycoprotein in Pichia pastoris. Virol Sin. 2007;22(5):347–52.

    Article  CAS  Google Scholar 

  32. Hsieh Y-C, Wang S-Y, Lee Y-F, Chen S-H, Mak PO, Chu C-Y. DNA sequence analysis of glycoprotein G gene of bovine ephemeral fever virus and development of a double oil emulsion vaccine against bovine ephemeral fever. J Vet Med Sci. 2006;68(6):543–8.

    Article  CAS  Google Scholar 

  33. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997;25(24):4876–82.

    Article  CAS  Google Scholar 

  34. Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35(6):1547.

    Article  CAS  Google Scholar 

  35. Nei M, Kumar S. Molecular evolution and phylogenetics: Oxford university press; 2000.

    Google Scholar 

  36. Walker PJ, Wang Y, Cowley JA, McWilliam SM, Prehaud CJ. Structural and antigenic analysis of the nucleoprotein of bovine ephemeral fever rhabdovirus. J Gen Virol. 1994;75(8):1889–99.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

The authors are very thankful to Shahid Chamran University of Ahvaz for the facilities to accomplish the present research project (Grant Number: SCU.SB98.658; Research Project No: 929).

Funding

No funding.

Author information

Authors and Affiliations

Authors

Contributions

All authors contributed to the design of the experiment. SER, FM and KhM carried out the experiments. Data analysis was performed by SER. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Seyedeh Elham Rezatofighi.

Ethics declarations

Ethics approval and consent to participate

Written consent was obtained from all farm owners before entering the study. Sampling was performed by a specialist veterinarian to minimize pain and injury to the animal. In addition, the current study is approved by the Ethics Committee of Shahid Chamran University of Ahvaz. All methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Supplementary Table 1.

Amino acid identity of BEFVs from different countries according to glycoprotein.

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

Rezatofighi, S.E., Mirzadeh, K. & Mahmoodi, F. Molecular characterization and phylogenetic analysis of bovine ephemeral fever viruses in Khuzestan province of Iran in 2018 and 2020. BMC Vet Res 18, 19 (2022). https://doi.org/10.1186/s12917-021-03119-x

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12917-021-03119-x

Keywords