Phylogenetic analysis and biological characteristics of an Akabane virus isolated in China

Akabane virus (AKAV) is an arbovirus that can cause miscarriage, premature birth, congenital malformations, and encephalomyelitis in young livestock. It is currently widely prevalent in China. Vero and MDBK cells were utilized to determine the virus titer and growth curve of the AKAV FS202301 strain in this study. Subsequently, the viral solution was intracranially (IC) or intraperitoneally (IP) inoculated into 8-day-old suckling mice, and the pathogenicity was explored by observing clinicopathological changes, and hematoxylin–eosin (HE) staining. Additionally, the S, M, and L segment sequences of the FS202301 strain were analyzed, phylogenetic trees were constructed, and antigenic epitopes were predicted to investigate its genetic variation. The results revealed that the viral titer of strain FS202301 was 10⁶ TCID₅₀/mL, with the number of viral copies peaking 24 h post-infection (hpi). This strain predominantly induced damage to the central nervous system, culminating in the death of suckling mice. It was classified as belonging to the genogroup Ia, exhibiting the highest degree of variation in the M segment, lower degrees of variation in the S and L segments, and no recombination events in any of the genes. The Gc protein encoded by the M segment contains more amino acid mutation sites and predicts a greater number of antigenic epitopes. This study aims to enhance the understanding of AKAV genetic variation in China and to establish a theoretical foundation for the future prevention and control of AKAV epidemics.

1 Introduction

Akabane disease (AKAD) is an insect-borne disease caused by the Akabane virus (AKAV) that leads to abortions, premature births, stillbirths, central nervous system damage in newborn fetuses, as well as congenital arthrogryposis-hydranencephaly syndrome in pregnant sheep and cows, and it is mainly transmitted by midges and mosquitoes (1–4). AKAV is classified within the family Bunyaviridae and the genus Orthobunyavirus, a single-stranded negative-stranded RNA virus characterized by a capsid, and a genome comprising three segments: S, M, and L (5, 6). The S segment is highly conserved and encodes the nucleoprotein (N) as well as the nonstructural protein (NSs); the M segment is highly variable and encodes a precursor polyprotein that is post-translationally cleaved to produce the viral envelope glycoproteins (Gn and Gc) and the nonstructural protein (NSm); and the L segment encodes solely the viral RNA polymerase (L protein) (7). Based on the M segment sequences, AKAV was classified into four distinct groups (I-IV), with genogroup I further subdivided into two subgroups: Ia and Ib (8). Kono et al. (8) demonstrated that AKAV genogroup Ia is more neuronophilic and is capable of causing encephalomyelitis in calves, as well as reproductive disorders in pregnant cows.

In 1959, AKAV was first isolated in Gunma Prefecture, Japan, and was subsequently reported in Australia, Israel, Indonesia, South Korea, and China (9, 10). Although no large-scale outbreak of AKAV has been documented in China to date. However, Wang et al. (11) conducted serological surveys across 24 provinces in China from 2006 to 2015, revealing that cattle and sheep tested seropositive for AKAV antibodies in most regions. This indicates a widespread prevalence of AKAV infection among cattle and sheep in China. In 1992, Liao et al. (12) isolated 11 strains of AKAV from dairy farms in Taiwan, China, where abortion and fetal malformations had been observed in pregnant cows, marking the first report. In 2010–2011, Cao et al. (13) isolated two strains of AKAV (HN10169 and HN10174) from mosquitoes, representing the first documented occurrence of AKAV isolation from field insects in China. Between 2013 and 2016, five novel AKAV strains were isolated from bamboo rats in Guangxi Province, China (14). In 2016, the AKAV strain TJ2016 was isolated from bovine serum in China, and sequence analysis of the S and M segments indicated that it belonged to genogroup II (15). In 2019, Tang et al. (10) successfully isolated the AKAV strain GXDH 01 from goat blood, with phylogenetic analysis indicating that this strain is classified within genogroup Ia. Furthermore, it has been demonstrated that AKAV can infect a variety of animals, including cattle, sheep, pigs, bamboo rats, camels, deer, and horses (10, 14, 16–19), resulting in significant economic losses to the global animal husbandry industry. This study presents an analysis of the genetic evolution and pathogenicity of the AKAV FS202301 strain (GenBank accession no. PQ567126-PQ567128) isolated and preserved in the lab (20), aimed at examining the genetic variation of this strain and providing a theoretical foundation for advancing molecular epidemiology and detection research related to AKAV.

2 Materials and methods

2.1 Virus titer determination

The virus solution of the FS202301 strain, which has been plaque-purified, was diluted to a 10-fold dilution and inoculated into a 96-well plate containing a monolayer of Vero cells, with 100 μL added to each well. For each dilution level, seven duplicate wells were set up, and a negative control was included to ensure the accuracy and reliability of the experiment. Following incubation at 37 °C in a 5% CO₂ environment for 5 to 7 days, the number of wells exhibiting cytopathic effects (CPE) at each dilution was observed using a microscope. This experiment was conducted three times, and the half tissue culture infectious dose (TCID₅₀) of the virus was calculated using the Reed-Muench method (21).

2.2 Determination of growth curves

The virus solution (infection complex MOI = 1) was inoculated into MDBK cells, and harvest viral RNA at 10 time points ranging from 0 to 84 h post-infection (hpi). The quantity of viral copies at each time point was quantified using the RT-qPCR assay. Viral proliferation curves were plotted using GraphPad Prism 9.1.0 (GraphPad, San Diego, CA).

2.3 Genome-wide analysis of strain FS202301

Fifteen reference sequences for AKAV S, M, and L were selected from the GenBank database on the NCBI website (Supplementary Tables 1–3). Nucleotide and amino acid similarity analyses, as well as amino acid sequence alignments, were conducted on the fragment sequences of the FS202301 strain using GraphPad Prism (Version 9.0.0) and the MegaAlign function of DNAStar software (Supplementary Figure 2).

2.4 Phylogenetic analysis

Phylogenetic analyses were performed on 225 selected AKAV S sequences, 108 AKAV M sequences, and 44 AKAV L sequences (Supplementary Tables 4–6). Phylogenetic trees were constructed using the maximum likelihood (ML) method in MEGA software (version 11.0.13, Mega Limited, Auckland, New Zealand). The ML method was repeated 1,000 times with bootstrap, while the remaining parameters were maintained at default values. The resulting phylogenetic trees were then imported into the online beautification editing software (The Interactive Tree of Life)¹ for subsequent sorting and editing.

2.5 Recombination detection

Preliminary identification of potential recombination events in the collected AKAV S, M, and L fragments was conducted using seven algorithms available in RDP software (version 4.0): RDP, GENECONV, BootScan, MaxChi, Chimera, SiScan, and 3Seq. The results of the analysis indicate that strains exhibiting four or more positive (+) results with a significance level of p < 0.05 are classified as recombinant strains.

2.6 Prediction of N, Gn, NSm, and Gc protein antigenic epitopes

Antigenic epitopes of N, Gn, NSm and Gc proteins of AKAV FS202301 strain were predicted using an online prediction site for antigenic epitopes.²

2.7 Pathogenicity study in suckling mice

To evaluate the pathogenicity of the FS202301 strain, 8-day-old suckling mice were divided into three groups (n = 8 per group), two of which were inoculated intracranially (IC) and inoculated intraperitoneally (IP) with 10 or 100 μL of 10⁶ TCID₅₀/mL viral fluids, and the negative control (NC) was injected with phosphate buffer solution (PBS). The mental status of suckling mice was monitored daily following inoculation. Suckling mice in the IC group were dissected at 1, 3, and 5 day post inoculation (dpi), while those in the IP group were dissected at 2, 5, 8, and 11 dpi to assess the pathological conditions of various organs (Euthanasia was performed on suckling mice by cervical dislocation). Following dissection, tissue specimens, including brain, heart, liver, spleen, lung, kidney, and intestine, were collected for pathological analysis. All experimental procedures were conducted as described by Gao et al. (22). Total RNA was extracted from brain, heart, liver, spleen, lung, kidney, and intestine tissues of suckling mice using TRIzol reagent (Sangon Biotech, Songjiang, Shanghai, China). RNA samples were reverse transcribed using Hiscript® III All-in-One RT SuperMix Perfect for qPCR (Vazyme, Nanjing, China). Primer sequences were as follows: forward (F): 5’-TAAGACGCCACAAC CAAGTGT-3′, reverse (R): 5’-CCGAAATGCGATGGAGCGTA-3′. The PCR reaction conditions included pre-denaturation at 95 °C for 2 min, followed by 40 cycles of amplification (95 °C for 10 s, 60 °C for 30 s). The AKAV NSS gene was selected as the target for viral genome copy number detection. This experiment received approval from the Animal Ethics Committee of Foshan University.

3 Result

3.1 Virus titer and determination of growth curve

After 48 h of the virus solution was inoculated with Vero cells, the virus titer of FS202301 strain was measured to be 10⁶ TCID₅₀/mL. By evaluating viral replication levels across distinct time intervals, we demonstrated that viral replication levels of the FS202301 strain exhibited an exponential increase from 0 to 24 hpi in MDBK cells, reaching a peak at 24 h. Subsequently, the viral copy number gradually declined and stabilized after 72 h (Supplementary Figure 1).

3.2 Genome-wide analysis of strain FS202301

Nucleotide and amino acid similarity analyses were conducted on each fragment of the FS202301 strain using GraphPad Prism (Version 9.0.0) and the MegaAlign function of the DNAStar software (Supplementary Figure 2). The results indicated that the S segment of strain FS202301 exhibited the highest nucleotide similarity (99.6%) with strain HN10169-2010 and the lowest nucleotide similarity (81.9%) with strain AN_9398–2018; the highest amino acid homology (99.6%) was observed with strains DHL10M110–2010, GD18240-China-2018 and HN10169-2010, while the lowest amino acid similarity (88.9%) was with AN_9398–2018 strain. The M segment of strain FS202301 exhibited the highest nucleotide and amino acid similarity with the HN10174-2010 strain, at 99 and 98.9%, respectively. Conversely, it displayed the lowest nucleotide and amino acid similarity to strain MP496-1972, at 70.2 and 74.3%. In contrast, the L segment of strain FS202301 demonstrated the highest nucleotide and amino acid similarity to strain HN10174-2010 (99.4 and 99.8%), the lowest nucleotide similarity to strain CS0016-1975 (86.2%), and the lowest amino acid similarity to strain R7949-1968 (95.9%).

The segments of strain FS202301 were analyzed and compared utilizing the MegAlign function of DNAStar software (Table 1). The results indicated that the S segment exhibited a high level of conservation, with only a single mutation identified at amino acid 221 (T211A). In contrast, the M segment harbored 40 amino acid mutations, with the Gc protein exhibiting the highest number of mutation sites, followed by the NSm protein, while the Gn protein exhibited the fewest mutation sites. Additionally, seven amino acid mutations were identified in the RdRp protein encoded by the L segment.

Table 1

Table 1. Mutation sites for S, M and L segments.

3.3 Phylogenetic tree and recombination analysis

To investigate the genetic variation of the various segments of the FS202301 strain, AKAV S, M, and L segments (225 S segments, 108 M segments, and 44 L segments) from different countries and years were selected from NCBI for genetic evolutionary analysis alongside the FS202301 strain, resulting in the construction of phylogenetic trees (Figure 1). The result revealed that the FS202301 strain is classified within genogroup Ia. In the phylogenetic trees based on the S and L segments, the FS202301 strain exhibited the closest genetic relationship to the GXLCH02-2016 strain. Conversely, in the phylogenetic tree constructed from the M segments, the FS202301 strain was found to have the closest genetic distance to the HN10714-2010 strain.

Figure 1

Figure 1. The construction of phylogenetic trees were performed using the ML method in MEGA software with 1000 bootstrap replicates. (a) Construct a phylogenetic tree based on the S segment. (b) Construct a phylogenetic tree based on the M segment. (c) Construct a phylogenetic tree based on the L segment.

To evaluate the presence of recombination events within the segmented sequences of the isolated FS202301 strain, RDP software (version 4.0) was utilized to analyze the multiple alignment of 226 AKAV S segments, 109 AKAV M segments, and 45 AKAV L segment sequences. The results of this analysis indicated that no recombination events were detected.

3.4 Prediction of N, Gn, NSm, and Gc protein antigenic epitopes

Antigen epitope prediction was conducted for the N, Gn, NSm, and Gc proteins of FS202301 strain. The results indicated that 9 antigenic epitopes were predicted on the N protein, 5 on the Gn protein, and 7 on the NSm protein, while approximately 21 antigenic epitopes were identified on the Gc protein (Supplementary Figure 3). The Gc protein exhibited the highest potential for antigenic epitopes.

3.5 Pathogenicity study in suckling mice

Following IC inoculation, six suckling mice exhibited symptoms such as anorexia, lethargy, ataxia, and moribund state by 3 dpi, with the remaining mice were in a state of moribund by 5 dpi. After 6-day post-IP inoculation, the suckling mice gradually exhibited symptoms such as depression, ataxia, incontinence, and anorexia. By 7–8 dpi, six suckling mice were near death, one was in a state of moribund by 11 dpi, and by 13 dpi, the remaining suckling mice were all in a moribund state (Figure 2b). Moribund suckling mice were euthanized and dissected. Gross pathological examination revealed that both IC and IP groups exhibited primary lesions localized to the brain, characterized by cerebral edema, pallor of brain tissue, effacement of cortical sulci, and smoothing of cerebral gyri (Figure 2a).

Figure 2

Figure 2. Pathogenicity assessment of FS202301 isolate in suckling mice: mortality, viral load dynamics, and tissue tropism. (a) Postmortem examination of suckling mice; (b) Survival curves of suckling mice; (c) Viral load dynamics in various organs of suckling mice at different time points. The data are shown as means ± SD (error bars). Asterisks (*) indicate a significant difference between groups (ns, p < 0.12; *, p < 0.033; **, p < 0.002; ***, p < 0.001).

RT-qPCR targeting the NSs gene was conducted on samples collected from the brain, heart, liver, spleen, lungs, kidneys, and intestines of suckling mice (Figure 2c). In the IC group, AKAV was detected exclusively in the brain and liver at 1 dpi, while viral presence was observed in all examined tissues (brain, heart, liver, spleen, lungs, kidneys, and intestines) at 3 and 5 dpi. Quantitative analysis revealed that the viral load in the brain was significantly higher than in other tissues, followed by the heart. In contrast, the IP group showed nearly undetectable viral levels at 2 dpi, followed by elevated viral loads in the intestines and heart at 5 dpi. By 8 and 11 dpi, the brain exhibited significantly higher viral loads compared to other tissues, with the heart showing the second highest levels, while remaining tissues maintained low viral loads. Brain tissue was selected as the representative organ for AKAV detection. The quantitative analysis demonstrated distinct viral replication patterns between the two groups: in the IC group, brain viral loads peaked at 3 dpi before gradually declining, whereas in the IP group, brain viral loads showed a rapid increase starting at 5 dpi, with the rate of increase slowing by 8 dpi (Figure 2c).

Histopathological examination through HE staining revealed distinct histopathological damage in the brains of suckling mice, which varied depending on the inoculation route (Figure 3). In the NC group, neither the brain nor the heart displayed notable pathological abnormalities. In the IC group, histopathological examination revealed distinct alterations in brain tissue architecture. Within the hippocampal formation, neurons in the CA2 region and dentate gyrus (DG) exhibited morphological changes, characterized by cellular shrinkage (indicated by dark blue arrows), increased basophilia, and loss of distinct nuclear-cytoplasmic boundaries. Similarly, a subset of neurons in the hypothalamus and brainstem demonstrated comparable pathological features. Cardiac tissue showed mild hydropic degeneration of cardiomyocytes (marked by red arrows), swollen cells, sparse and lightly stained cytoplasm and interstitium, including connective tissue and blood vessels, without obvious necrosis or other abnormalities. In the IP group, histopathological analysis revealed distinct neuropathological alterations. Focal neuronal necrosis (indicated by black arrows) was occasionally observed in the hippocampal CA1 region, thalamus, hypothalamus, and brainstem, characterized by nuclear fragmentation and pyknosis. The hippocampal CA2 region exhibited sporadic vacuolar degeneration in pyramidal cells (marked by green arrows), manifesting as rounded cytoplasmic vacuoles. Notably, no significant pathological changes, such as cardiomyocyte degeneration or necrosis, were observed in the myocardial tissue.

Figure 3

Figure 3. Histopathological analysis of FS202301 strain inoculated in suckling mice.

4 Discussion

AKAV is an arbovirus primarily transmitted through Culicoides and mosquitoes (23–25), widely prevalent in tropical and temperate climate regions such as Africa, Asia, the Middle East, and Australia (22), posing a significant threat to livestock farming, particularly in cattle and sheep populations. A serological survey on Chinese cattle and sheep revealed that AKAV is widely present in China, with a significantly higher incidence observed in cattle compared to sheep and goats (22). Currently, genogroup Ia strains are predominantly identified in Japan and China, while genogroup Ib strains are primarily located in Japan and Israel; genogroup II strains are mainly distributed in Japan, genogroup III strains are chiefly observed in Australia, and genogroup IV strains are primarily found in South Korea (15, 26). Notably, genogroup Ia strains exhibit the highest prevalence, resulting in significant economic losses to the farming industry.

In this study, nucleotide and amino acid similarity analyses of the S, M, and L segment sequences of the FS202301 strain compared to the reference strain demonstrated that the nucleotide similarity of the S segment ranged from 81.9 to 99.6%, while the amino acid similarity ranged from 88.9 to 99.6%. The nucleotide similarity of the M segment ranged from 70.2 to 99%, and the amino acid similarity ranged from 74.3 to 98.9%. For the L segment, the nucleotide similarity was between 86.2 and 99.4%, whereas the amino acid similarity was between 95.9 and 99.8%. These results indicate that the nucleotide and amino acid differences in the M segment are the most pronounced, reflecting a relatively high degree of variability. In contrast, the nucleotide and amino acid similarities between the S and L segments show smaller differences, suggesting a reduced level of variation. The M segment is recognized as the genomic segment exhibiting the greatest variability among Bunyaviruses, whereas the S segment demonstrates a lower degree of variation. The RdRp encoded by the L segment, which is essential for viral genome replication and transcription, has a conserved structure (27–30). This indicates that the results of this study are consistent with existing research findings.

Phylogenetic analysis of the S, M, and L segments revealed that the FS202301 strain is classified within the genogroup Ia. Among them, phylogenetic analysis of the S segment indicated that FS202301 is genetically closest to the GXLCH02 strain, isolated in 2016. Conversely, phylogenetic analyses based on the M and L segments demonstrated that FS202301 shares greater genetic proximity with the HN10174 strain, isolated in 2010. Previous studies have demonstrated that AKAV genogroup Ia is capable of causing encephalomyelitis in calves, exhibiting a neurotoxicity that is greater than that of genogroup II strains (31). Despite the high seroprevalence of AKAV detected in China, infected animals often present with mild symptoms or remain cryptically infected, with only transient fever and mild neurological disorders observed during the peak breeding season of Culicoides and mosquitoes (22). The FS202301 strain has been classified as belonging to genogroup Ia. It was hypothesized that, similar to other strains within genogroup Ia, it may induce encephalomyelitis in infected animals. The present study examined the pathogenicity of FS202301 strain in 8-day-old suckling mice, and all of the mice exhibited severe neurological symptoms following infection, consistent with previously reported findings. Comparative analysis revealed that IC inoculation induced more rapid disease progression and greater severity compared to IP inoculation, likely attributable to the direct viral penetration of the blood–brain barrier, resulting in immediate neural tissue invasion and subsequent neurological damage. Quantitative viral load assessment across multiple organs at various time points demonstrated that the brain consistently maintained the highest viral titers, which showed progressive increase throughout the infection course. In addition, due to the different modes of virus inoculation, AKAV was detected in all organs of suckling mice in the IP group at different times, whereas in the IC group, AKAV was detected only in the brain at 2 dpi, and the rest of the organs at 3 and 5 dpi, which indicated that the AKAV was systemically infected; and the viral loads of AKAV in the IC group were higher than those of the IP group in all organs. The results of HE staining reveal that both the IC and IP groups exhibited pathological damage in the brain, further confirming the brain as a primary target of AKAV infection. These findings provide compelling evidence that the brain serves as the primary target for AKAV infection and replication.

In this study, a comparative analysis of the amino acid sequences of the segments of the FS202301 strain revealed that the M segment contained the most amino acid mutation sites, whereas the S segment exhibited the fewest. In orthobunyaviruses, the glycoproteins Gc and Gn, encoded by the M segment, assemble into spikes on the viral particles and play crucial roles in viral attachment, cell fusion, and the induction of host immune responses (32). The Gc protein exhibits high immunogenicity and serves as the principal target for neutralization, capable of inducing the production of neutralizing antibodies and facilitating attachment to mammalian cells; the Gn protein is involved in viral attachment to insects and other arthropods, whereas the NSm protein may be implicated in viral assembly and morphogenesis (33–35). The majority of the mutated amino acid sites in the M segment of the FS202301 strain are situated within the Gc proteins, and modifications at these sites may have significant implications for their immunogenicity. Ishihara et al. (36) constructed multiple NSm deletion mutants using reverse genetics and assessed their pathogenicity, identifying several regions that are critical for viral infectivity (amino acids 323–460, 369–447, 323–367, and 447–460). They found that the NSm protein influences AKAV replication by functioning as a virulence factor. In this study, the majority of amino acid mutations identified in the FS202301 NSm protein were located within the region encompassing amino acids 323–460 (T323I, Q364R, H386D, I391T, I409V, V460I, and N496S), hypothesized that mutations at these sites may significantly influence its pathogenicity.

The AKAV N protein is the most abundant protein in the virion and infected cells, which induces the production of complement and triggers complement fixation antibodies after infection of the animal host, and has become widely utilized in virus or antibody detection (37–39). Additionally, there are at least five antigenic regions on the Gc protein, which is the main neutralizing protein of the virus (40). In this study, the antigenic epitopes of the FS202301 strain were predicted, and it was discovered that the Gc protein exhibited the highest number of predicted antigenic epitopes, followed by the N protein. These findings provide a theoretical foundation for the development of monoclonal antibodies targeting the AKAV Gc and N proteins, as well as for future innovations in vaccine and detection methods.

In this study, the viral titer and growth curve of the FS202301 strain were assessed. The results indicated that the viral titer was measured at 10⁶ TCID₅₀/mL, with the viral copy number peaking at 24 h post-infection, followed by a gradual decline and stabilization after 72 h. Pathogenicity tests conducted on 8-day-old suckling mice revealed that the FS202301 strain primarily induced neurological damage. Furthermore, phylogenetic analyses of the S, M, and L segments of strain FS202301 indicated that it is classified within genogroup Ia and exhibits a close genetic relationship with strains GXLCH02 and HN10174. The strain did not undergo recombination. Through comparative analysis of amino acid sequences and the prediction of antigenic epitopes, we identified a greater number of amino acid mutation sites and antigenic epitopes on the Gc protein. These findings provide a theoretical basis for comprehending the potential genetic variations of AKAV and for the development of monoclonal antibodies in China in the future.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary material.

Ethics statement

The animal study was approved by the Animal Ethics Committee of Foshan University. The study was conducted in accordance with the local legislation and institutional requirements.

Author contributions

GL: Formal analysis, Methodology, Writing – review & editing, Investigation, Writing – original draft. SC: Investigation, Writing – review & editing, Writing – original draft, Formal analysis, Methodology. XL: Methodology, Visualization, Writing – review & editing. HL: Investigation, Writing – review & editing, Methodology. LQ: Writing – review & editing, Methodology, Investigation. DZ: Methodology, Investigation, Writing – review & editing. JH: Software, Writing – review & editing, Validation, Formal analysis. PL: Formal analysis, Software, Writing – review & editing, Validation. YuZ: Formal analysis, Writing – review & editing, Software, Validation. YaZ: Writing – review & editing, Formal analysis, Software, Validation. MZ: Writing – review & editing, Validation. HG: Validation, Writing – review & editing. WT: Funding acquisition, Resources, Writing – review & editing, Data curation. XZ: Writing – review & editing, Data curation, Funding acquisition, Resources. BS: Writing – review & editing, Funding acquisition, Resources, Data curation. WQ: Resources, Funding acquisition, Writing – review & editing, Data curation. MZ: Funding acquisition, Supervision, Writing – review & editing, Project administration. KZ: Project administration, Writing – review & editing, Funding acquisition, Supervision.

Funding

The author(s) declare that financial support was received for the research and/or publication of this article. This research was supported by National Key Research and Development Program Projects (2023YFD1801302, 2023YFD1801301); Guangdong Basic and Applied Basic Research Foundation (2024A1515140152); Special Funds for Rural Revitalization Strategy of Guangdong Provincial Department of Agriculture and Rural Affairs in 2024 (2024CXTD15); Open Foundation of Guangdong Provincial Key Laboratory of Animal Molecular Design and Precise Breeding(2024B03); Base and Talent Program of Science and Technology Plan in Tibet Autonomous Region (XZ202401JD0012); Foshan University High Level Talents Research Start-up Project (CGZ07465/CGZ07001).

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.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

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.2025.1691476/full#supplementary-material

Footnotes

1. ^https://itol.embl.de

2. ^http://tools.immuneepitope.org/bcell/

References

1. Charles, JA. Akabane virus. Vet Clin North Am Food Anim Pract. (1994) 10:525–46.

Google Scholar

2. Inaba, Y, Kurogi, H, and Omori, T. Letter: Akabane disease: epizootic abortion, premature birth, stillbirth and congenital arthrogryposis-hydranencephaly in cattle, sheep and goats caused by Akabane virus. Aust Vet J. (1975) 51:584–5.

Google Scholar

3. Kirkland, PD. Akabane and bovine ephemeral fever virus infections. Vet Clin North Am Food Anim Pract. (2002) 18:501–14. doi: 10.1016/s0749-0720(02)00026-9

PubMed Abstract | Crossref Full Text | Google Scholar

4. St George, TD, Cybinski, D, and Paull, NI. The isolation of Akabane virus from a normal bull. Aust Vet J. (1977) 53:249.

Google Scholar

5. Kuhn, JH, Adkins, S, Alkhovsky, SV, Avšič-Županc, T, Ayllón, MA, Bahl, J, et al. 2022 taxonomic update of phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders Bunyavirales and Mononegavirales. Arch Virol. (2022) 167:2857–906. doi: 10.1007/s00705-022-05546-z

PubMed Abstract | Crossref Full Text | Google Scholar

6. Pattnaik, AK, and Abraham, G. Identification of four complementary RNA species in Akabane virus-infected cells. J Virol. (1983) 47:452–62.

Google Scholar

7. Kobayashi, T, Yanase, T, Yamakawa, M, Kato, T, Yoshida, K, and Tsuda, T. Genetic diversity and reassortments among Akabane virus field isolates. Virus Res. (2007) 130:162–71. doi: 10.1016/j.virusres.2007.06.007

PubMed Abstract | Crossref Full Text | Google Scholar

8. Kono, R, Hirata, M, Kaji, M, Goto, Y, Ikeda, S, Yanase, T, et al. Bovine epizootic encephalomyelitis caused by Akabane virus in southern Japan. BMC Vet Res. (2008) 4:20. doi: 10.1186/1746-6148-4-20

PubMed Abstract | Crossref Full Text | Google Scholar

9. Oya, A, Okuno, T, Ogata, T, Kobayashii,, and Matsuyama, T. Akabane, a new arbor virus isolated in Japan. Jpn J Med Sci Biol. (1961) 14:101–8.

Google Scholar

10. Tang, HB, Ren, P, Qin, S, Lin, J, Bai, A, Qin, S, et al. Isolation, genetic analysis of the first Akabane virus from goat in China. J Vet Med Sci. (2019) 81:1445–9. doi: 10.1292/jvms.18-0602

PubMed Abstract | Crossref Full Text | Google Scholar

11. Wang, J, Blasdell, KR, Yin, H, and Walker, PJ. A large-scale serological survey of Akabane virus infection in cattle, yak, sheep and goats in China. Vet Microbiol. (2017) 207:7–12. doi: 10.1016/j.vetmic.2017.05.014

PubMed Abstract | Crossref Full Text | Google Scholar

12. Liao, YK, Lu, YS, Goto, Y, and Inaba, Y. The isolation of Akabane virus (Iriki strain) from calves in Taiwan. J Basic Microbiol. (1996) 36:33–9.

Google Scholar

13. Cao, Y, Fu, S, Song, S, Cai, L, Zhang, H, Gao, L, et al. Isolation and genome phylogenetic analysis of arthropod-borne viruses, including Akabane virus, from mosquitoes collected in Hunan Province, China. Vector Borne Zoonotic Dis. (2019) 19:62–72. doi: 10.1089/vbz.2018.2267

PubMed Abstract | Crossref Full Text | Google Scholar

14. Tang, HB, Chen, F, Rao, G, Bai, A, Jiang, J, Du, Y, et al. Characterization of Akabane virus from domestic bamboo rat. Southern China Vet Microbiol. (2017) 207:280–5. doi: 10.1016/j.vetmic.2017.06.018

PubMed Abstract | Crossref Full Text | Google Scholar

15. Chen, D, Wang, D, Wei, F, Kong, Y, Deng, J, Lin, X, et al. Characterization and reverse genetic establishment of cattle derived Akabane virus in China. BMC Vet Res. (2021) 17:349. doi: 10.1186/s12917-021-03054-x

PubMed Abstract | Crossref Full Text | Google Scholar

16. Al-Busaidy, S, Hamblin, C, and Taylor, WP. Neutralising antibodies to Akabane virus in free-living wild animals in Africa. Trop Anim Health Prod. (1987) 19:197–202.

Google Scholar

17. Davies, FG, and Jessett, DM. A study of the host range and distribution of antibody to Akabane virus (genus bunyavirus, family Bunyaviridae) in Kenya. J Hyg (Lond). (1985) 95:191–6.

Google Scholar

18. Huang, CC, Huang, TS, Deng, MC, Jong, MH, and Lin, SY. Natural infections of pigs with akabane virus. Vet Microbiol. (2003) 94:1–11. doi: 10.1016/s0378-1135(03)00062-2

PubMed Abstract | Crossref Full Text | Google Scholar

19. Konno, S, Moriwaki, M, and Nakagawa, M. Akabane disease in cattle: congenital abnormalities caused by viral infection. Spontaneous disease. Vet Pathol. (1982) 19:246–66.

Google Scholar

20. Zhang, M, Deng, M, Zhao, S, Zhao, D, Zheng, Y, Qin, L, et al. Isolation and identification of an AKAV strain in dairy cattle in China. Front Vet Sci. (2025) 12:1574667. doi: 10.3389/fvets.2025.1574667

PubMed Abstract | Crossref Full Text | Google Scholar

21. Smither, SJ, Lear-Rooney, C, Biggins, J, Pettitt, J, Lever, MS, and Olinger, GG Jr. Comparison of the plaque assay and 50% tissue culture infectious dose assay as methods for measuring filovirus infectivity. J Virol Methods. (2013) 193:565–71. doi: 10.1016/j.jviromet.2013.05.015

PubMed Abstract | Crossref Full Text | Google Scholar

22. Gao, HF, Wang, JP, Yang, ZX, Xie, JR, He, YW, Hong, QH, et al. Genetic and pathogenic characterisation of a virulent Akabane virus isolated from goats in Yunnan, China. J Vet Res. (2022) 66:35–42. doi: 10.2478/jvetres-2022-0007

PubMed Abstract | Crossref Full Text | Google Scholar

23. Allingham, PG, and Standfast, HA. An investigation of transovarial transmission of Akabane virus in Culicoides brevitarsis. Aust Vet J. (1990) 67:273–4.

Google Scholar

24. Jennings, M, and Mellor, PS. Culicoides: biological vectors of Akabane virus. Vet Microbiol. (1989) 21:125–31.

Google Scholar

25. Yanase, T, Kato, T, Kubo, T, Yoshida, K, Ohashi, S, Yamakawa, M, et al. Isolation of bovine arboviruses from Culicoides biting midges (Diptera: Ceratopogonidae) in southern Japan: 1985–2002. J Med Entomol. (2005) 42:63–7. doi: 10.1093/jmedent/42.1.63

PubMed Abstract | Crossref Full Text | Google Scholar

26. Taylor, WP, and Mellor, PS. The distribution of Akabane virus in the Middle East. Epidemiol Infect. (1994) 113:175–85.

Google Scholar

27. Akashi, H, Kaku, Y, Kong, XG, and Pang, H. Sequence determination and phylogenetic analysis of the Akabane bunyavirus S RNA genome segment. J Gen Virol. (1997) 78:2847–51.

Google Scholar

28. Goller, KV, Höper, D, Schirrmeier, H, Mettenleiter, TC, and Beer, M. Schmallenberg virus as possible ancestor of Shamonda virus. Emerg Infect Dis. (2012) 18:1644–6. doi: 10.3201/eid1810.120835

PubMed Abstract | Crossref Full Text | Google Scholar

29. Poch, O, Sauvaget, I, Delarue, M, and Tordo, N. Identification of four conserved motifs among the RNA-dependent polymerase encoding elements. EMBO J. (1989) 8:3867–74.

Google Scholar

30. Purnomo Edi, S, Ibrahim, A, Sukoco, R, Bunali, L, Taguchi, M, Kato, T, et al. Molecular characterization of an Akabane virus isolate from West Java. Indonesia J Vet Med Sci. (2017) 79:774–9. doi: 10.1292/jvms.17-0009

PubMed Abstract | Crossref Full Text | Google Scholar

31. Oem, JK, Yoon, HJ, Kim, HR, Roh, IS, Lee, KH, Lee, OS, et al. Genetic and pathogenic characterization of Akabane viruses isolated from cattle with encephalomyelitis in Korea. Vet Microbiol. (2012) 158:259–66. doi: 10.1016/j.vetmic.2012.02.017

PubMed Abstract | Crossref Full Text | Google Scholar

32. Wernike, K, Aebischer, A, Sick, F, Szillat, KP, and Beer, M. Differentiation of antibodies against selected Simbu serogroup viruses by a glycoprotein Gc-based triplex ELISA. Vet Sci. (2021) 8:12. doi: 10.3390/vetsci8010012

PubMed Abstract | Crossref Full Text | Google Scholar

33. Ludwig, GV, Israel, BA, Christensen, BM, Yuill, TM, and Schultz, KT. Role of La Crosse virus glycoproteins in attachment of virus to host cells. Virology. (1991) 181:564–71.

Google Scholar

34. Roman-Sosa, G, Brocchi, E, Schirrmeier, H, Wernike, K, Schelp, C, and Beer, M. Analysis of the humoral immune response against the envelope glycoprotein Gc of Schmallenberg virus reveals a domain located at the amino terminus targeted by mAbs with neutralizing activity. J Gen Virol. (2016) 97:571–80. doi: 10.1099/jgv.0.000377

PubMed Abstract | Crossref Full Text | Google Scholar

35. Shi, X, Kohl, A, Leonard, VH, Li, P, McLees, A, and Elliott, RM. Requirement of the N-terminal region of orthobunyavirus nonstructural protein NSm for virus assembly and morphogenesis. J Virol. (2006) 80:8089–99. doi: 10.1128/JVI.00579-06

PubMed Abstract | Crossref Full Text | Google Scholar

36. Ishihara, Y, Shioda, C, Bangphoomi, N, Sugiura, K, Saeki, K, Tsuda, S, et al. Akabane virus nonstructural protein NSm regulates viral growth and pathogenicity in a mouse model. J Vet Med Sci. (2016) 78:1391–7. doi: 10.1292/jvms.16-0140

PubMed Abstract | Crossref Full Text | Google Scholar

37. Elliott, RM. Orthobunyaviruses: recent genetic and structural insights. Nat Rev Microbiol. (2014) 12:673–85. doi: 10.1038/nrmicro3332

PubMed Abstract | Crossref Full Text | Google Scholar

38. Saeed, MF, Li, L, Wang, H, Weaver, SC, and Barrett, ADT. Phylogeny of the Simbu serogroup of the genus Bunyavirus. J Gen Virol. (2001) 82:2173–81. doi: 10.1099/0022-1317-82-9-2173

PubMed Abstract | Crossref Full Text | Google Scholar

39. Wang, J, Chen, D, Wei, F, Deng, J, Su, J, Lin, X, et al. Generation of stable cell lines expressing Akabane virus N protein and insight into its function in viral replication. Pathogens. (2023) 12:1058. doi: 10.3390/pathogens12081058

PubMed Abstract | Crossref Full Text | Google Scholar

40. Ogawa, Y, Eguchi, M, and Shimoji, Y. Two Akabane virus glycoprotein Gc domains induce neutralizing antibodies in mice. J Vet Med Sci. (2022) 84:538–42. doi: 10.1292/jvms.21-0670

PubMed Abstract | Crossref Full Text | Google Scholar