A Novel Ephemero- and a New CHeRI Orbivirus Isolated from a Dead Farmed White-Tailed Deer (Odocoileus virginianus) in Florida, USA

Abstract

A novel ephemeral fever rhabdovirus and a CHeRI orbivirus of a previously unidentified genetic lineage were isolated in mosquito cell line C6/36 cells as co-infecting agents from the spleen tissue of a dead farmed white-tailed deer (WTD; Odocoileus virginianus) in Florida. We designated the ephemeral fever rhabdovirus as Hardee County ephemerovirus 1, strain CHeRI ephemerovirus 1. The genetic sequences of the CHeRI orbivirus isolated in this work differ significantly from those of three previously described CHeRI orbivirus lineages. We designated this new virus as CHeRI orbivirus 4, strain CHeRI orbivirus 4-1. Whereas it remains unknown whether one, both, or none of the viruses contributed to the pathology, gross observations revealed that the dead WTD had severely congested and hemorrhagic lungs, and that its heart, kidneys, and spleen were also congested.

1. Introduction

The population of farmed deer in the United States has increased by 35% in the last five years, and in 2022, it was determined to comprise approximately 287,000 animals [1]. Florida has an estimated 140 deer farms and 11,000 farmed deer, consisting mostly of white-tailed deer (WTD; Odocoileus virginianus) [1]. The University of Florida Cervidae Health Research Initiative (CHeRI) collaborates with WTD-farming stakeholders within the state of Florida to increase the health and sustainable production of captive cervids, as well as the native wildlife and ecosystems in which they live [2]. As elsewhere, Florida’s farmed deer are susceptible to various serious illnesses that are caused by vector-borne pathogens. Orbiviruses such as bluetongue virus (BTV) and epizootic hemorrhagic disease virus (EHDV) can cause high morbidity and mortality in farmed and wild populations of WTD in North America [3,4]. Scientists affiliated with CHeRI recently discovered CHeRI orbivirus and found Big Cypress, Mobuck and Yunnan orbiviruses in dead farmed WTD, suggesting that these viruses may also be significant pathogens of this cervid species [5,6]. Whereas BTV and EHDV are arthropod-borne viruses that are transmitted by biting midges of the genus Culicoides [4,7], the vector of CHeRI orbiviruses has not yet been identified.

The family Rhabdoviridae is made up of four subfamilies and consists of at least 434 negative-sense single-stranded RNA viruses [8]. Viruses from this family can infect a broad range of hosts, including mammals, birds, reptiles, and fish, and have also been detected in arthropods [8]. Rhabdoviruses are important pathogens of humans, livestock, fish, and agricultural crops, and this group of viruses includes rabies virus and vesicular stomatitis virus, which are zoonotic agents [8]. Ephemerovirus is a genus of viruses from the subfamily Alpharhabdovirinae and is made up of a large group of arthropod-borne rhabdoviruses [8]. Rhabdoviridae genomes include five genes (N, P, M, G, and L) encoding the structural proteins and multiple additional long open reading frames (ORFs) between the G gene and L gene [8,9]. The best characterized ephemerovirus is bovine ephemeral fever virus (BEFV). It is the causative agent of bovine ephemeral fever (BEF), which is an economically important disease of cattle and water buffalo in Africa, Australia, and Asia. Cattle and water buffalo with BEF experience a sudden onset of fever, lameness, anorexia, and ruminal stasis, followed by a sustained drop in milk production [8]. Mortality rates associated with BEF are usually low (1–2%) but have been reported to exceed 20% [8,10,11,12]. Ephemeroviruses had been isolated exclusively from cattle, pigs, and hematophagous insects until now [13,14]. We have designated this new virus as Hardee County ephemerovirus 1, strain CHeRI ephemerovirus.

CHeRI partners discovered CHeRI orbiviruses that group into three distinct lineages (CHeRI orbivirus types 1–3) from specimens taken from dead WTD [5], wherein it was observed that specimens containing these viruses were commonly co-infected with other viruses [5,6]. We now report the discovery of a fourth CHeRI orbivirus lineage, which we have designated CHeRI orbivirus 4, strain CHeRI orbivirus 4-1.

The discovery of the two novel viruses in this work adds to our understanding of the pathogens causing potentially deadly illnesses in WTD. Moreover, our findings exemplify the need for continued pathogen surveillance to identify viruses previously not described in farmed WTD populations and, in all likelihood, in wild WTD.

2. Materials and Methods

2.1. Animal History and Specimens Collected for Diagnostic Tests

On 7 July 2023, a farmed WTD (animal ID: OV1850) from Hardee County, Florida, USA, was found dead. Mucoid feces was observed near the animal and presumed to have been shed by it. Alopecia was noted around its eyes and behind its right ear, possibly as a result of a postmortem fire-ant attack. A field necropsy was performed on the same day the animal was discovered dead, revealing severely congested and hemorrhagic organs. Aliquots of heart, kidney, liver, lung, and spleen tissues (CT, KT, HT, LT, and ST, respectively) were collected and stored at −80 °C for virology tests at a later time.

2.2. Preparation of Tissue and Skin Homogenates for Diagnostic Tests

After being thawed, samples of the tissue specimens of approximately 30 mg were excised and placed in a 2.0 mL conical bottom microcentrifuge tube (Benchmark Scientific Inc., Sayreville, NJ, USA) containing 0.75 g of high-density zirconium oxide 2 mm beads and 0.15 g of high-density zirconium oxide 0.1 mm beads (Glen Mills Inc., Clifton, NJ, USA). Thereafter, 600 µL of lysis buffer made up of Buffer RLT (RNeasy Mini Kit, Qiagen, Valencia, CA, USA) and dithiothreitol (Thermo Fisher Scientific, Waltham, MA, USA) was added for a final concentration of 40 µM before the samples were homogenized for two minutes using a Mini-Beadbeater-16 (Model 607, BioSpec Products, Bartlesville, OK, USA).

2.3. Extraction and Purification of RNA from Tissue Homogenates

Total RNA was extracted from the tissue homogenates using an RNeasy Mini Kit, according to the manufacturer’s protocol (Qiagen, Valencia, CA, USA). The RNA was eluted from the RNeasy kit’s silica columns into 100 µL of Ambion^TM RNA Storage Solution (Thermo Fisher Scientific) and stored at −80 °C until further analysis. RNA concentrations were determined using a Nanodrop One (Thermo Fisher, Scientific), then diluted to 50 ng/µL and screened for bovine viral diarrhea virus (BVDV), BTV, and EHDV genomic RNAs using a One Step RT-PCR kit (Applied Biosystems, Waltham, MA, USA) as described in Wernike and colleagues [15]. The RT-qPCR assays were performed using a QuantStudio 5 Real Time PCR system (Applied Biosystems) with the following cycling conditions: reverse transcription at 48 °C for 10 min, inactivation of the reverse transcriptase and initial denaturation at 95 °C for 10 min, followed by 40 cycles of three-step cycling consisting of 95 °C for 15 s, 57 °C for 45 s, and 68 °C for 45 s. Amplification and multicomponent graphs were subsequently assessed: Ct values of 36 and under were considered positive, 37–39 suspect (equivocal), and 40 and above negative. Total RNA purified from the homogenates tested negative during initial RT-qPCR screens for BTV, BVDV, and EHDV (Table 1).

Since gross pathology observations indicated hemorrhagic disease, it was possible that other orbiviruses detected in WTD may have been involved in this case. The RNA samples were therefore tested with a one-step multiplex qPCR test for CHeRI orbiviruses 1–3, Mobuck virus (MOV), Big Cypress orbivirus (BCOV), and Yunnan orbivirus (YUOV) that we developed for the CHeRI program (Table 2). The primers and probes of the two newly developed one-step multiplex qPCRs were designed and optimized according to the TaqMan Multiplex Optimization User Guide for optimum assay efficiency (MAN0010189, MAN0014269, Applied Biosystems) and using PrimerExpress v2.0 with default settings. The primer and probe combinations with the lowest penalty values were selected for the multiplex assay and are summarized in Table 2. For both the CHeRI orbivirus 1–3 and the MOV–BCOV–YUOV multiplex assays, RNA extracts were screened using a TaqPath™ 1-Step Multiplex Master Mix kit (Applied Biosystems) with primer concentrations of 72 µM each, and probes of 20 µM. 16 µL PCR mix and 4 µL RNA, diluted to 12.5 ng/µL, were loaded into a MicroAmp Fast 96-Well Reaction Plate (0.1 mL) (Applied Biosystems) and sealed with a MicroAmp Optical Adhesive Film (Applied Biosystems). The plate was then run on a QuantStudio 5 Real Time PCR system (Applied Biosystems). The cycling conditions for both CHeRI orbivirus 1–3 and the MOV–BCOV–YUOV multiplex assays are as follows: 25 °C for 2 min, 53 °C for 10 min, and 95 °C for 2 min, followed by 40 cycles of 95 °C for 3 s and 60 °C for 30 s. Amplification and multicomponent graphs were assessed as previously described. The samples were tested in duplicate (Table 1).

2.4. Cell Cultures

Virus isolation was attempted in Vero E6 (Cercopithecus aethiops [African green monkey]) (ATCC, Manassas, VA, USA, Cat. no. ATCC CRL1586) and C6/36 cells (Aedes albopictus [Asian tiger mosquito] CRL Cat/no. 1660) obtained from the American Type Culture Collection (ATCC). The cells were propagated as monolayers in 25 cm² vented tissue culture flasks (25 cm² flask, Corning Inc., Corning, NY, USA) using Advanced Dulbecco’s Modified Eagle’s Medium (aDMEM, Invitrogen Corp. Thermo Fisher Scientific, Waltham, MA, USA) supplemented with 2 mM L-alanyl-L-glutamine (GlutaMAX^TM, Invitrogen Corp.), antibiotics (PSN; 50 μg/mL penicillin, 50 μg/mL streptomycin, 100 μg/mL neomycin [Invitrogen Corp.]), and 10% low-antibody, heat-inactivated, gamma-irradiated fetal bovine serum (FBS, Hyclone, GE, Healthcare Life Sciences, Pittsburgh, PA, USA). Vero E6 cells were incubated at 37 °C, and C6/36 cells at 28 °C, in 5% CO₂ atmospheres within humidified incubators.

2.5. Inoculation of Cell Cultures with Tissue Homogenates

Fifty µL of the tissue homogenates were added to 3 mL of supplemented aDMEM and filtered through a 0.45 µm pore-size syringe-tip filter (Grainger, Lake Forest, IL, USA) to remove contaminating bacteria and fungi. The resulting filtrates were then used to inoculate confluent monolayers of Vero E6 and C6/36 in 25 cm² vented tissue culture flasks (Corning Inc.). Mock-inoculated cells were maintained in parallel with the inoculated flasks. The inoculated cells were monitored for formation of virus-induced cytopathic effects (CPEs) using an inverted microscope with phase-contrast optics, with refeeds of the cells performed every 3 days. Aliquots of the spent cell culture media of cells displaying CPEs were collected and stored at –80 °C for follow-up analyses at a future time.

2.6. Next-Generation Sequencing (NGS)

The spent cell culture media of the C6/36 cells inoculated with ST was chosen for analyses, based on our previous work [5]. After thawing on ice, RNA was extracted from the virions in the spent growth media using a QIAamp Viral RNA Mini Kit (Qiagen, Valencia, CA, USA) according to the manufacturer’s protocol. A cDNA library was generated using a NEBNext Ultra RNA Library Prep kit (Illumina, San Diego, CA, USA) and sequenced on an Illumina NextSeq 1000 sequencer. Cell culture host sequences were removed using Kraken v2.0 [16], with Ae. albopictus genome sequences (GCA_001876365.2) as reference. Thereafter, de novo assembly of the remaining paired-end reads was performed using MEGAHIT v1.1.4 [17]. The assembled contigs were subjected to BLASTX searches against the National Center for Biotechnology Information (NCBI) nonredundant protein database using OmicsBox v1.2.

2.7. Phylogenetic Analyses

Maximum likelihood (ML) phylogenetic trees were constructed to assess the evolutionary relationships of Hardee County ephemerovirus 1 to other ephemeroviruses, and of CHeRI orbivirus 4-1 to other orbiviruses. For ephemeroviruses, the amino acid sequence of the deduced L protein from Hardee County ephemerovirus 1 was aligned with the sequences from 17 other ephemeroviruses, along with a sequence from one tibrovirus (which served as an outgroup), all of which were retrieved from the NCBI GenBank database. The alignments were performed using MAFFT [18], and the ML trees were constructed in IQ-TREE v1.6.12, with 1000 bootstrap replicates performed to test the robustness of the clades [19]. To determine whether Hardee County ephemerovirus 1 meets the International Committee on Taxonomy of Viruses (ICTV) classification criteria for the Ephemerovirus genus, sequence identity matrices for the nucleotide sequences of the L and N genes and the amino acid sequences of their corresponding proteins were generated using the Sequence Demarcation Tool v1.2 [20]. Similarly, ML phylogenetic trees were constructed for CHeRI orbivirus 4-1 using nucleotide and amino acid alignments of the RNA-dependent RNA polymerase (VP1), the major outer capsid protein (VP2), and the innermost sub-core capsid T2 protein (VP3), along with sequences from 33 other orbiviruses retrieved from the NCBI GenBank database. St. Croix River virus was included as an outgroup. Sequence identity matrices were generated for nucleotide sequences of these genes and the amino acid sequences of their corresponding proteins using the Sequence Demarcation Tool v1.2.

3. Results

3.1. Gross Examinations

The lungs of farmed WTD (animal ID: OV1850) were severely congested and hemorrhagic, and its heart, spleen, and kidneys were congested (Figure 1). The trachea was filled with white foam and liquid. However, we found no significant hemorrhage in the oral cavities or the rectum. The spleen appeared normal externally, with the capsule and organ maintaining appropriate size and shape; however, internally, the organ was dominated by congested gelatinous red pulp (Figure 1). Other internal organs appeared normal, such as the liver, stomach, and intestines, but were still collected, as is standard practice within CHeRI. These findings indicate viral infection may take place primarily in the affected organs in comparison to the organs and systems that seem normal in contrast.

3.2. Evidence of Virus Isolation in Cultured Cells

Virus-induced CPEs were observed in C6/36 cells by 10 days post-inoculation (dpi) with LT and ST homogenates. In comparison, the Vero E6 cells were passaged at 15 dpi and maintained until 20 dpi, at which point samples of cell supernatant were obtained and frozen at −80 °C for further analyses; no CPEs were observed in the Vero E6 cells by 20 dpi. The CPEs included darkening of the cell cytoplasm, followed by detachment of dead cells from the growing surface of the cell culture flasks. A representative image is shown in Figure 2.

3.3. Genomic Sequence Analyses

Next-generation sequencing generated 71,362,254 reads, of which 97.6% were classified as host genome sequences that were subsequently removed using Kraken v2.0. The remaining 1,736,566 (2.4%) of paired-end reads were subjected to de novo assembly. The BLASTX analysis recovered a genome of Hardee County ephemerovirus 1, with a length of 14,682 bp and an average coverage of 987 reads per nucleotide. Additionally, a novel orbivirus genome, designated as CHeRI orbivirus 4-1, was identified. This genome consists of 10 segments, with a total length of 18,819 bp and an average coverage of 9688 reads per nucleotide. The genome sequences of both viruses are available at GenBank: (a) Hardee County ephemerovirus 1 (GenBank no. PQ480188) and CHeRI orbivirus 4-1 (PQ471672–PQ471681). A schematic rendition of the ephemerovirus genome is shown in Figure 3. A description of the genomic positions, gene products, and amino acid lengths of the ephemerovirus genome is defined within Table 3.

3.4. Phylogenetic, Amino Acid, and Nucleotide Sequence Analyses of Hardee County Ephemerovirus 1 and CHeRI Orbiviruses 4-1

Maximum likelihood (ML) analysis based on L protein (RdRp) amino acid alignments supports Hardee County ephemerovirus 1 as a sister species to Yata virus, with strong bootstrap support (89%) (Figure 4A). Similarly, the amino acid sequence identity analyses reveal that the L protein Hardee County ephemerovirus 1 and Yata virus share the highest identity (59.2%) among the ephemerovirus proteins that were analyzed (Figure 4B). However, the N protein shares 52.7% identity with that of Yata virus, and its highest identity similarity (58.9%) is with that of Kotonkan virus (Figure 4C). These findings suggest that Hardee County ephemerovirus 1 represents a genetically distinct lineage within the Ephemerovirus genus.

Phylogenetic analysis of CHeRI orbivirus 4-1, based on both nucleotide and amino acid sequences of the RNA-dependent RNA polymerase (VP1) (Figure 5A), major outer capsid protein (VP2) (Figure 5B), and innermost sub-core capsid T2 protein (VP3) (Figure 5C), support its classification as a unique member within the CHeRI orbivirus group. Sequence identity analyses of the nucleotide sequences of each gene and the amino acid sequences of their corresponding proteins (Figure 5D–F) further confirm these phylogenetic relationships. The nucleotide sequence identity (lower half of each matrix) reveals that CHeRI orbivirus 4-1 shares the highest VP1 gene identity (74.8%) with CHeRI orbivirus 3-1, while the amino acid sequence identity (upper half of each matrix) is 86.4%. The VP2 gene exhibits the highest nucleotide identity (57.3%) and amino acid identity (48.5%) with CHeRI orbivirus 2-1. The VP3 (T2) gene follows a similar pattern to VP1, with nucleotide identities of 64.5% to 76% and amino acid identities of 63.9% to 85.3% observed in CHeRI orbivirus 3 variants.

4. Discussion

In this manuscript, we report the discovery of a novel ephemerovirus, Hardee County ephemerovirus 1, and a novel orbivirus, CHeRI orbivirus 4-1. Based on the ICTV species demarcation criteria for ephemerovirus classification, Hardee County ephemerovirus 1 exhibits >15% amino acid divergence in the L protein and >8% divergence in the N protein compared to its closest known ephemeroviruses. These levels of divergence exceed the ICTV thresholds, supporting the classification of Hardee County ephemerovirus 1 as a novel species within the Ephemerovirus genus [8]. According to ICTV criteria for the Orbivirus genus, CHeRI orbivirus 4-1 is classified within this genus because its deduced RdRp (VP1) amino acid sequences share over 30% identity with those of other orbiviruses. Furthermore, ICTV classification relies on T2 protein (VP3) amino acid sequences to delineate distinct orbivirus species, with species demarcated by a threshold of <91% sequence identity. CHeRI orbivirus 4-1 exhibits the highest amino acid identity (85.3%) with CHeRI orbivirus 3-3, which is well below the species threshold. Additionally, the nucleotide identity of the VP2 gene is less than 74%, further distinguishing CHeRI orbivirus 4-1 from related species and meeting the ICTV criteria for a new species [5,21]. Phylogenetic analysis further supports the distinctiveness of CHeRI orbivirus 4-1 by placing it in a separate clade from its closest known relatives, the CHeRI orbivirus 3 variants. Taken together, these findings provide genetic and phylogenetic evidence that CHeRI orbivirus 4-1 represents a novel species within the genus Orbivirus.

To our knowledge, this is the first time that an ephemerovirus has been identified in a WTD. Whereas some ephemeroviruses are known causes of severe illness in bovids, we are unaware of reports implicating these viruses as cervid pathogens. However, there are similarities in the gross pathology presentations in deer affected by CHeRI orbiviruses 1-3 [5,6] with those observed in the animal that was the subject of this report. Some of the animals infected with CHeRI orbiviruses in previous cases also presented with foam within the respiratory system, a spleen with a gelatinous interior, and normal GI tracts [5]. Of interest, all reported detections of CHeRI orbiviruses (including this one) were identified as co-infections. Yata virus, the ephemerovirus most related to Hardee County ephemerovirus 1 (Figure 4), has been isolated twice from Ma. uniformis mosquitos but has not been associated with the infection of vertebrates [22]. An ephemerovirus that is closely related to the Hardee County ephemerovirus 1 is New Kent County virus (Ephemerovirus kent) (Figure 4), which was identified in deer ticks (Ixodes scapularis) through metagenomic sequencing [8,15]. This suggests the possibility that deer ticks are the vector of Hardee County ephemerovirus 1.

Our future goals are to understand the role of the CHeRI orbivirus 4-1 and Hardee County ephemerovirus 1 in the disease(s) of deer. The discovery of these two viruses has important implications for our research through the CHeRI program, but even more so for WTD farmers in Florida. Florida WTD farmers rely on accurate and up-to-date findings regarding biological threats facing their WTD herds to inform their vaccination, insemination, and pest management systems. As this animal was infected with two viruses, it is not clear whether one or a combination of the two was the cause of death. More work must be conducted to determine whether Hardee County ephemerovirus 1 has an etiological role similar to BEFV, potentially placing other farming infrastructure in Florida at risk. It is likely that the vector of this virus is an arthropod, but more work will need to be completed to establish whether hematophagous insects (such as biting midges or mosquitos) or other arthropods such as ticks are competent vectors of the virus.

Foremost among the next steps into understanding the significance of these two viruses in deer pathology are to isolate the viruses from each other from the co-infected C6/36 cells and to establish specific RT-qPCR tests for them. Following that, we must find mammalian cells as well as arthropod cells other than C6/36 cells wherein the two viruses can complete their life cycles. It may be necessary to use primary deer, biting midge, and tick cells. The ability to propagate these viruses in suitable cells in vitro will allow us to establish preliminary growth kinetics studies. An interesting question now is whether these viruses are emerging or re-emerging viruses, and insights can be attained by testing suitable archived deer specimens.

Author Contributions

Conceptualization: E.D., P.H.O.V., J.A.L., S.M.W. and K.S.; methodology: E.D., P.H.O.V., Z.W., T.C.S.R., P.K., J.M.C.K., J.A.L., S.M.W., K.S. and T.C.S.R.; formal analysis: E.D., P.H.O.V., K.W., J.M.C.K., P.K., T.C.S.R., J.M.C.K., J.A.L., K.S. and T.B.W.; investigation: E.D., P.H.O.V., K.W., Z.W., A.R., MU, P.K., J.M.C.K., S.M.W., K.S. and J.A.L.; resources: E.D., A.R., M.U., P.H.O.V., K.W., P.K., J.M.C.K., S.M.W., K.S., J.A.L. and T.B.W.; writing—original draft: E.D., K.W., P.K., J.M.C.K., S.M.W., K.S. and J.A.L.; writing—review and editing: E.D., P.H.O.V., K.W., P.K., J.M.C.K., J.A.L., S.M.W., K.S., Z.W. and T.B.W.; visualization: E.D., P.H.O.V., K.W., P.K., J.M.C.K., J.A.L., S.M.W. and K.S.; supervision: J.M.C.K., J.A.L., S.M.W. and K.S.; project administration: J.M.C.K., J.A.L., S.M.W. and K.S.; funding acquisition: J.M.C.K., J.A.L., S.M.W. and K.S. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

This work was approved by the University of Florida Institutional Animal Care and Use Committee (IACUC), with protocol numbers 201609390 and 201909390.

Informed Consent Statement

Not applicable.

Data Availability Statement

The complete gene coding sequences for all 10 segments of the CHeRI orbivirus genomes can be found under GenBank accession Nos. PQ471672–PQ4716881 and for the rhabdovirus under GenBank accession number PQ480188.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:

WTD	White-tailed deer
CHeRI	Cervidae Health Research Initiative
BTV	Bluetongue virus
EHDV	Epizootic hemorrhagic disease virus
BEFV	Bovine ephemeral fever virus
MOV	Mobuck virus
BCOV	Big Cypress orbivirus
YUOV	Yunnan orbivirus
aDMEM	Advanced Dulbecco’s Modified Eagle’s Medium
PSN	Penicillin–streptomycin–neomycin
FBS	Fetal bovine serum
CPE	Cytopathic effects
DPI	Days post inoculation
NCBI	National Center for Biotechnology Information
ML	Maximum likelihood

Footnotes

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Figures and Tables

Figure 1 Gross observations of farmed white-tailed deer naturally infected with Hardee County ephemerovirus 1 and CHeRI orbivirus-4. (A) Heart, external view. (B) Lung, internal view. (C) Spleen, external and internal views.

Figure 2 Cytopathic effects in C6/36 cells inoculated with ST homogenate from animal OV1850. (A) Mock-inoculated C6/36 cells, 10 dpi. (B) C6/36 cells inoculated with ST homogenate, 10 dpi. Original images were taken at 400× magnification.

Figure 3 Hardee County ephemerovirus 1. Schematic representation of the genome of strain CHeRI ephemerovirus 1 shown in reverse (positive-sense) polarity. The approximate positions of the structural ORFs are identified as N, P, M, G, and RdRp, which encode for the nucleo-, phospho-, matrix-, glyco-, and RNA-dependent RNA polymerase proteins. The non-structural ORFs are labeled are labeled as GNS, α1, α2, β, γ, and δ. The GNS ORF encodes a transmembrane glycoprotein, the α1 ORF encodes a viroporin, whereas the roles/functions of the α2, β, γ, and δ proteins are not known.

Figure 4 Phylogenetic and sequence identity analyses of Hardee County ephemerovirus 1. (A) Maximum likelihood phylogenetic tree based on the amino acid sequence of the L protein of Hardee County ephemerovirus 1 (highlighted in red) and related ephemeroviruses. Bootstrap values are indicated at each node, and the branch lengths represent the number of inferred substitutions, as shown by the scale. Ekpoma virus 2 was used as an outgroup. (B,C) Amino acid sequence identity matrices for the L protein (B) and N protein (C). The matrices display the percentage sequence identities between Hardee County ephemerovirus 1 (highlighted in red) and related ephemeroviruses.

Figure 5 Phylogenetic and sequence identity analyses of CHeRI orbivirus 4-1. (A–C) Maximum likelihood phylogenetic trees based on the amino acid sequences of the RNA-dependent RNA polymerase (VP1) (A), major outer capsid protein (VP2) (B), and innermost sub-core capsid T2 protein (VP3) (C) of CHeRI orbivirus 4-1 (highlighted in red) and related orbiviruses. Bootstrap values are shown at each node, and branch lengths represent the number of inferred substitutions, as indicated by the scale bar. (D–F) Sequence identity matrices for the RNA-dependent RNA polymerase (VP1) (D), major outer capsid protein (VP2) (E), and innermost sub-core capsid T2 protein (VP3) (F) between CHeRI orbivirus 4-1 (highlighted in red) and related orbiviruses. Amino acid sequence identities are presented in the upper half of each matrix, while the corresponding nucleotide sequence identities are shown in the lower half.

Table 1

Results of preliminary RT-qPCR screening tests.

Virus Tested	Ct Values
BVDV	Ct > 40
BTV	Ct > 40
EHDV	Ct > 40
CHeRI orbiviruses 1–3	Ct > 40
YUOV	Ct > 40
MOV	Ct > 40
BCOV	Ct > 40

Table 2

Sequence of primers and probes used in this study to identify orbiviruses 1–3.

Species	Target	Oligonucleotide ID	Sequence (5′ to 3′)
CHeRI orbv 1	VP7	CHeRI orbv1-F	CGCGCATTCTCCGTTAGATC
		CHeRI orbv1-R	GGTAATCCCGCTTGTCCATTT
		CHeRI orbv1-P	FAM- TGAGGTTCCACAAGGAG-MGB
CHeRI orbv 2	NS3	CHeRI orbv2-F	ATGCATGCAGCCATTTCAAC
		CHeRI orbv2-R	GGGACAACAGCTAGCATTGGA
		CHeRI orbv2-P	VIC-ATGAAACGTGAACAGACG-MGB
CHeRI orbv 3	NS1	CHeRI orbv3-F	CGGAACGCTCGTGCATAGAT
		CHeRI orbv3-R	TGCCTGTCAAACTCGCAAGT
		CHeRI orbv3-P	JUN-TCTACGTGGAGAGGGAA-QSY
YUOV	VP4	YUOV-F	TGGAGAGAGGACAAAGCAAATG
		YUOV-R	TGCTGCGCGCACATCT
		YUOV-P	VIC-ACCATGGCATCGAGTTG-MGB
BCOV	NS1	BCOV-F	CTAAGACACGTGGGCGATTGT
		BCOV-R	ACGCATTCCGCGTCTGAA
		BCOV-P	FAM-AGTGTGGCATGCTGG-MGB
MOV	VP7	MOV-F	AGCAAACAGAGATCCACGAG
		MOV-R	GTCTACTGTACTAGGCGTTTGTG
		MOV-P	JUN-TGATTGGTCTCTCGGTTGTGGG-QSY

Table 3

Genomic and deduced protein features of Hardee County ephemerovirus 1 (GenBank PQ480188.1).

Gene	Length (Ribonucleotides)	Genomic Positions	Gene Product	Amino Acid Length
Leader sequence	39	1–39	N/A a	N/A
N	1308	40–1347	Nucleoprotein	435
P	870	1372–2241	Phosphoprotein	289
M	672	2267–2938	Matrix	223
G	1815	2979–4793	Glycoprotein	604
GNS	1584	4820–6403	Nonstructural glycoprotein	527
Alpha1	324	6419–6742	Alpha1 protein	107
Alpha2	420	6642–7061	Alpha2 protein	139
Beta	441	7090–7530	Beta protein	146
Gamma	309	7578–7895	Gamma protein	102
Delta	309	7903–8214	Delta protein	102
L	6360	8255–14,614	Polymerase	2119
Trailer sequence	68	14,615–14,682	N/A	N/A

^a N/A; not applicable.

References

1. United States Department of AgricultureNational Agricultural Statistics Service. 2022 Census of Agriculture; Summary and State Data; National Agricultural Statistics Service: Washington, DC, USA, 2024; 24.

2. CHeRI Cervidae Health Research Initiative. Available online: https://wec.ifas.ufl.edu/cheri/ (accessed on 15 November 2024).

3. Clarke, L.L.; Mead, D.G.; Ruder, M.G.; Howerth, E.W.; Stallknecht, D. North American Arboviruses and White-Tailed Deer (Odocoileus virginianus): Associated Diseases and Role in Transmission. Vector-Borne Zoonotic Dis.; 2022; 22, pp. 425-442. [DOI: https://dx.doi.org/10.1089/vbz.2022.0005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35867036]

4. Rivera, N.A.; Varga, C.; Ruder, M.G.; Dorak, S.J.; Roca, A.L.; Novakofski, J.E.; Mateus-Pinilla, N.E. Bluetongue and Epizootic Hemorrhagic Disease in the United States of America at the Wildlife–Livestock Interface. Pathogens; 2021; 10, 915. [DOI: https://dx.doi.org/10.3390/pathogens10080915] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/34451380]

5. Ahasan, M.S.; Subramaniam, K.; Campos Krauer, J.M.; Sayler, K.A.; Loeb, J.C.; Goodfriend, O.F.; Barber, H.M.; Stephenson, C.J.; Popov, V.L.; Charrel, R.N. . Three New Orbivirus Species Isolated from Farmed White-Tailed Deer (Odocoileus virginianus) in the United States. Viruses; 2019; 12, 13. [DOI: https://dx.doi.org/10.3390/v12010013] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31861885]

6. Rodrigues, T.C.S.; Lednicky, J.A.; Loeb, J.C.; Campos Krauer, J.M.; Wisely, S.M.; Waltzek, T.B.; Subramaniam, K. Genome Sequence of a CHeRI Orbivirus 3 Strain Isolated from a Dead White-Tailed Deer (Odocoileus virginianus) in Florida, USA. Microbiol. Resour. Announc.; 2020; 9, e00523-20. [DOI: https://dx.doi.org/10.1128/MRA.00523-20] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32586868]

7. McGregor, B.L.; Erram, D.; Alto, B.W.; Lednicky, J.A.; Wisely, S.M.; Burkett-Cadena, N.D. Vector Competence of Florida Culicoides Insignis (Diptera: Ceratopogonidae) for Epizootic Hemorrhagic Disease Virus Serotype-2. Viruses; 2021; 13, 410. [DOI: https://dx.doi.org/10.3390/v13030410] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33807536]

8. Walker, P.J.; Freitas-Astúa, J.; Bejerman, N.; Blasdell, K.R.; Breyta, R.; Dietzgen, R.G.; Fooks, A.R.; Kondo, H.; Kurath, G.; Kuzmin, I.V. . ICTV Virus Taxonomy Profile: Rhabdoviridae 2022: This Article Is Part of the ICTV Virus Taxonomy Profiles Collection. J. Gen. Virol.; 2022; 103, [DOI: https://dx.doi.org/10.1099/jgv.0.001689] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35723908]

9. Walker, P.J.; Firth, C.; Widen, S.G.; Blasdell, K.R.; Guzman, H.; Wood, T.G.; Paradkar, P.N.; Holmes, E.C.; Tesh, R.B.; Vasilakis, N. Evolution of Genome Size and Complexity in the Rhabdoviridae. PLoS Pathog; 2015; 11, e1004664. [DOI: https://dx.doi.org/10.1371/journal.ppat.1004664] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25679389]

10. Hsieh, Y.-C.; Chen, S.-H.; Chou, C.-C.; Ting, L.-J.; Itakura, C.; Wang, F.-I. Bovine Ephemeral Fever in Taiwan (2001–2002). J. Vet. Med. Sci.; 2005; 67, pp. 411-416. [DOI: https://dx.doi.org/10.1292/jvms.67.411] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15876792]

11. Tonbak, S.; Berber, E.; Yoruk, M.D.; Azkur, A.K.; Pestil, Z.; Bulut, H. A Large-Scale Outbreak of Bovine Ephemeral Fever in Turkey, 2012. J. Vet. Med. Sci.; 2013; 75, pp. 1511-1514. [DOI: https://dx.doi.org/10.1292/jvms.13-0085]

12. Zheng, F.; Qiu, C. Phylogenetic Relationships of the Glycoprotein Gene of Bovine Ephemeral Fever Virus Isolated from Mainland China, Taiwan, Japan, Turkey, Israel and Australia. Virol. J.; 2012; 9, 268. [DOI: https://dx.doi.org/10.1186/1743-422X-9-268] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23150932]

13. Walker, P.J.; Klement, E. Epidemiology and Control of Bovine Ephemeral Fever. Vet. Res.; 2015; 46, 124. [DOI: https://dx.doi.org/10.1186/s13567-015-0262-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26511615]

14. Wu, Q.; Yang, Z.; Lu, Z.; Mi, S.; Feng, Y.; He, B.; Zhu, G.; Gong, W.; Tu, C. Identification of Two Novel Ephemeroviruses in Pigs Infected by Classical Swine Fever Virus. Infect. Genet. Evol.; 2022; 100, 105273. [DOI: https://dx.doi.org/10.1016/j.meegid.2022.105273] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/35321840]

15. Wernike, K.; Hoffmann, B.; Beer, M. Simultaneous Detection of Five Notifiable Viral Diseases of Cattle by Single-Tube Multiplex Real-Time RT-PCR. J. Virol. Methods; 2015; 217, pp. 28-35. [DOI: https://dx.doi.org/10.1016/j.jviromet.2015.02.023] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25746154]

16. Wood, D.E.; Salzberg, S.L. Kraken: Ultrafast Metagenomic Sequence Classification Using Exact Alignments. Genome Biol.; 2014; 15, R46. [DOI: https://dx.doi.org/10.1186/gb-2014-15-3-r46] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24580807]

17. Li, D.; Liu, C.-M.; Luo, R.; Sadakane, K.; Lam, T.-W. MEGAHIT: An Ultra-Fast Single-Node Solution for Large and Complex Metagenomics Assembly via Succinct de Bruijn Graph. Bioinformatics; 2015; 31, pp. 1674-1676. [DOI: https://dx.doi.org/10.1093/bioinformatics/btv033] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25609793]

18. Katoh, K.; Rozewicki, J.; Yamada, K.D. MAFFT Online Service: Multiple Sequence Alignment, Interactive Sequence Choice and Visualization. Brief. Bioinform.; 2019; 20, pp. 1160-1166. [DOI: https://dx.doi.org/10.1093/bib/bbx108] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28968734]

19. Nguyen, L.-T.; Schmidt, H.A.; von Haeseler, A.; Minh, B.Q. IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies. Mol. Biol. Evol.; 2015; 32, pp. 268-274. [DOI: https://dx.doi.org/10.1093/molbev/msu300] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25371430]

20. Muhire, B.M.; Varsani, A.; Martin, D.P. SDT: A Virus Classification Tool Based on Pairwise Sequence Alignment and Identity Calculation. PLoS ONE; 2014; 9, e108277. [DOI: https://dx.doi.org/10.1371/journal.pone.0108277] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25259891]

21. Perinet, L.C.; Mutebi, J.-P.; Powers, A.M.; Lutwama, J.J.; Mossel, E.C. Yata Virus (Family Rhabdoviridae, Genus Ephemerovirus) Isolation from Mosquitoes from Uganda, the First Reported Isolation since 1969. Diseases; 2023; 11, 21. [DOI: https://dx.doi.org/10.3390/diseases11010021] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36810535]

22. Matthijnssens, J.; Attoui, H.; Bányai, K.; Brussaard, C.P.D.; Danthi, P.; del Vas, M.; Dermody, T.S.; Duncan, R.; Fang, Q.; Johne, R. . ICTV Virus Taxonomy Profile: Sedoreoviridae 2022. J. Gen. Virol.; 2022; 103, 001782. [DOI: https://dx.doi.org/10.1099/jgv.0.001782] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/36215107]