Establishment of a Replicon Reporter of the Emerging Tick-Borne Bourbon Virus and Use It for Evaluation of Antivirals

Bourbon virus (BRBV) was first isolated from a patient hospitalized at the University of Kansas Hospital in 2014. Since then, several deaths have been reported to be caused by BRBV infection in the Midwest and Southern United States. BRBV is a tick-borne virus that is widely carried by lone star ticks. It belongs to genus Thogotovirus of the Orthomyxoviridae family. Currently, there are no treatments or vaccines available for BRBV or thogotovirus infection caused diseases. In this study, we reconstituted a replicon reporter system, composed of plasmids expressing the RNA-dependent RNA polymerase (RdRP) complex (PA, PB1, and PB2), nucleocapsid (NP) protein, and a reporter gene flanked by the 3′ and 5′ untranslated region (UTR) of the envelope glycoprotein (GP) genome segment. By using the luciferase reporter, we screened a few small molecule compounds of anti-endonuclease that inhibited the nicking activity by parvovirus B19 (B19V) NS1, as well as FDA-approved drugs targeting the RdRP of influenza virus. Our results demonstrated that myricetin, an anti-B19V NS1 nicking inhibitor, efficiently inhibited the RdRP activity of BRBV and virus replication. The IC50 and EC50 of myricetin are 2.22 and 4.6 μM, respectively, in cells. Myricetin had minimal cytotoxicity in cells, and therefore the therapeutic index of the compound is high. In conclusion, the BRBV replicon system is a useful tool to study viral RNA replication and to develop antivirals, and myricetin may hold promise in treatment of BRBV infected patients.


INTRODUCTION
Bourbon virus (BRBV) is a member in the genus Thogotovirus of the Orthomyxoviridae family. There are seven genera in the family of Orthomyxoviridae, which are segmented negative-strand RNA viruses (Kawaoka and Palese, 2006), including four types of influenza viruses (influenza virus A, B, C, and D), thogotovirus, quaranjavirus, and isavirus. Many viruses in the Orthomyxoviridae family are important pathogens to humans or animals. The epidemic/pandemic influenza A viruses have caused millions of human deaths in the past century and are still circulating, posing a huge threat to human health and the economy. Thogotoviruses mainly circulate in domestic animals, such as sheep, cattle, and camels, causing neural diseases and abortion. Two thogotoviruses, thogoto virus, and dhori virus, have been reported to infect humans and cause deaths (Butenko et al., 1987;Kosoy et al., 2015). Human antibodies against thogoto virus and dhori virus have been detected in Europe, Asia, and Africa (Filipe et al., 1985;Hubalek and Rudolf, 2012). It has been reported that thogoto virus can cause human infections in a vector-free manner, possibly by an aerosol route (Butenko et al., 1987), highlighting the potential to infect humans in a large population.
Bourbon virus was first isolated from a blood sample of a hospitalized male patient from Bourbon County, Kansas, United States, in the spring of 2014 at the University of Kansas Hospital (Kosoy et al., 2015). The patient died due to a complex syndrome of leukopenia, lymphopenia, thrombocytopenia, hyponatremia, and increased levels of aspartate aminotransferase and alanine aminotransferase. Due to the high level of viremia and unique identification of the virus in the serum taken from the patient, BRBV was believed to be the cause of the illness and the death of the patient. Later, two additional cases of BRBV infection have been reported. In 2015, a patient from Payne County, Oklahoma, United States tested positive for neutralization antibodies to BRBV before fully recovering (Savage et al., 2017). In June 2017, a 58-year-old woman from Missouri died from an infection of BRBV after she had been misdiagnosed for a significant period of time (Bricker et al., 2019). Not surprisingly, a high seroprevalence of BRBVneutralizing antibodies in raccoons (50%) and white-tailed deer (86%) was detected in MO, United States (Jackson et al., 2019). However, the Center for Disease Control and Prevention (CDC), United States, currently does not know if the virus may be found in other areas of the United States, since a seroprevalence of BRBV infection has not been evaluated outside of these epidemic regions. Nevertheless, BRBV is the first species of the genus Thogotovirus of the Orthomyxoviridae family to be identified as a human pathogen in the New World .
Bourbon virus is widely carried by the lone star tick (Amblyomma americanum), a species that is aggressive, feeds on humans, and is widely distributed across the East, Southeast, and Midwest States. BRBV was found in three pools of lone star ticks from retrospective tests in 39,096 ticks from northwestern Missouri, 240 km from Bourbon County, Kansas, United States (Savage et al., 2017). The human BRBV and the strain isolated from tick pools share >99.0% sequence at the amino acid level and 95.0% identity at the RNA sequence level (Savage et al., 2017). BRBV replicates in cell lines derived from the hard ticks Amblyomma, Hyalomma, and Rhipicephalus . Together with the geographic location of the BRBV infection and the geographic distribution of a number of A. americanum ticks (Savage et al., 2017(Savage et al., , 2018, these studies strongly suggest that the lone star tick is a vector of BRBV transmission to humans. Little is known about the biology of BRBV. Due to the high mortality of BRBV infection, specific treatments with antiviral drugs are in demand to save lives. In this study, we established a replicon reporter of BRBV. We then used the replicon reporter to examine anti-influenza drugs that target viral RNA-dependent RNA polymerase (RdRP) and endonuclease inhibitors of human parvovirus B19 (B19V) for inhibition of BRBV RdRP activity and BRBV replication in HEK293T and Vero cells.

Plaque Assays
Cells were seeded in six-well plates at a density of 2 × 10 6 cells and were confluent the following day. We used the cell growth media to serially dilute the virus stock at 10fold. 200 µl of the diluent were added to each well and incubated for 1 h with gently rocking the plate every 20 min.

pcDNA3-Based Viral cDNA Plasmids
Bourbon virus RNA was isolated from a stock of BRBV-KS passaged on Vero cells. cDNA of each viral RNA genome fragment was synthesized using the reverse primer and amplified by PCR. The amplification program started with one cycle at 48 • C for 45 min and one cycle at 94 • C for 2 min. These cycles were followed by 40 cycles at 94 • C for 20 s, 52 • C for 30 s, and 72 • C for 40 s; the program ended with one cycle at 72 • C for 5 min. The PCR products were visualized by agarose gel electrophoresis, extracted, and cloned into pcDNA3 vector (Invitrogen), resulting in pcDNA-(BRBV) PA, PB1, PB2, M, NP, and glycoprotein (GP) plasmids.

pHW2000-GP Plasmid
The plasmid pHW2000 (Hoffmann et al., 2002) contains the human ribosome RNA polymerase I (pol I) promoter and the murine terminator sequence separated by two BsmBI sites. The pol I promoter and terminator elements are flanked by a truncated immediate-early promoter of cytomegalovirus (CMV) and by the bovine growth hormone polyadenylation signal (bGHpA). The pcDNA3-GP was used as a template to amplify the GP cDNA with a primer set that has the BsmBI sites (Table 1).
After digestion of the PCR products with BsmBI, the fragments were cloned into the pHW2000 vector, which resulted in the pHW2000-(BRBV)GP plasmid.

Reporter Plasmids
pHW-CMV-GP was constructed by deletion of the CMV promoter in pHW2000-GP. Then, we constructed pHW-CMV-GP-GFP and pHW-CMV-GP-gLuc plasmids inserting the GFP open reading frame (ORF) through the NdeI and EagI sites of the GP ORF and replacing the GP ORF with the Gaussia luciferase (gLuc) ORF, respectively, into the pHW-CMV-GP plasmid ( Figure 4A). The sequences of all the primers are listed in Table 1. All the constructed plasmids were sequenced to confirm the cloned cDNA sequences at MCLAB (South San Francisco, CA, United States). DNA sequence analyses were performed using SnapGene 4.0 (SnapGene, Chicago, IL, United States).

FIGURE 2 | Comparisons of BRBV-KS infection between Vero and HEK293T cells. (A)
BRBV growth kinetics. Cells were infected with BRBV-KS at an MOI = 10 pfu/cell. Cell culture media (supernatants) were collected daily for 6 days. Viral RNA was extracted from the samples of supernatants and then titrated by a RT-qPCR assay. The data present the means with standard deviations, which are obtained from three independent experiments, with each experiment analyzing samples in duplicate. (B) Western blotting.
∼5 × 10 4 infected cells were lyzed in 1 × Laemmli sample buffer and loaded in SDS-10% polyacrylamide gel for separation. The gel was transferred to a polyvinylidene difluoride membrane, which was blotted with an anti-BRBV M protein rat polyclonal antibody. The M protein in BRBV-infected Vero and HEK293T cells was detected at ∼30 kDa as indicated, but not in mock-infected cells. β-actin was probed as a loading control. Original images can be found in "Supplementary Material".

Plasmid DNA Transfection
HEK293T cells were transfected using the PEImax (Cat# 24765-2, MW 40,000, Polysciences, Inc.) at a ratio of 1:3 of DNA: PEI in 200 µl of opti-MEM (Invitrogen) . The total amounts of plasmid DNA were kept constant (2-3 µg per well of six-well plate) in each group by supplementation with an empty vector.

BRBV Reporter Assays
BRBV GFP Reporter Assay pHW-CMV-GP-GFP was co-transfected with the four pcDNA-PA, PB1, PB2, and NP plasmids into HEK293T cells confluent in six-well plates. After 7 days post-transfection, the GFP signal was observed under a fluorescent microscope (Nikon Ti-S). The transfection of pHW-CMV-GP-GFP with three pcDNA plasmids (PB1, PB2, NP) was set as a negative control.

BRBV Luciferase Reporter Assay
pHW-CMV-GP-gLuc and the four pcDNA-based plasmids were co-transfected into HEK293T cells in 96-well plate. At 3 days post-transfection, 10 µl of the cell culture supernatants were taken and mixed with 50 µl of Working solution or luciferase activity, using the gLuc Flash Assay kit (Pierce TM , #16159). Luminescent signals were detected on a Synergy H1 microplate reader (BioTek, Winooski, VT, United States). Three pcDNA plasmids (PB1, PB2, and NP) and pHW-CMV-GP-gLuc cotransfected cells were set up as a negative control.

Half Cytotoxic Concentration (CC 50 ) Assay
The cytotoxicity of compounds on HEK293T and Vero cells was evaluated using the CytoTox-Glo TM Cytotoxicity Assay kit (Promega, #G9290), according to the manufacturer's instructions. Cells were seeded at 4 × 10 4 cells per well in 96-well plates. After overnight incubation, different concentrations of the compounds diluted in culture media were added to each well. Cells treated with DMSO were used for the negative control. Cell viability was measured after 3 days of incubation.

Half Maximal Effective Concentration (EC 50 ) Assay
Confluent monolayers of Vero and HEK293T cells in 96well plates were inoculated with 10 pfu of BRBV. Different concentrations of the compounds, diluted in medium, were added to the wells. Each concentration of the compounds was tested in duplicate per experiment, and the experiment was repeated three times. At 3 days post-infection, culture supernatants were collected and used to quantify the nuclease digestion-resistant vgc numbers using reverse transcription and quantitative PCR (RT-qPCR).

Western Blotting
Cells were collected by centrifugation and lysed using a 1 ml tuberculin syringes. The cell lysates were loaded, along with 2.5 µl of a pre-stained protein ladder (#P008; GoldBio, FIGURE 4 | Establishment of a BRBV replicon report system. (A) Diagram of the reporter system. In a CMV pol II-deleted pHW-CMV-GP plasmid, a GFP open read frame (ORF) or a gLuc ORF was inserted into the pol I expression cassette to obtain pHW-CMV-GP-GFP and pHW-CMV-GP-gLuc. In cells transfected with pcDNA-PA, PB1, PB2, and NP plasmids, pol l transcribes negative sense viral (reporter) RNA [(-)rRNA], which is used as a template for replication of (-)rRNA and production of reporter mRNA (GP-GFP or GP-gLuc) by the RdRP complex and NP protein as diagrammed. Then reporter mRNA is translated to express GFP or gLuc. Since the CMV pol II promoter was removed, the production of reporter mRNA is solely dependent on the replication of the reporter RNA (the activity of the RdRP complex). "C" at the 5 end of the (+) reporter mRNA represents the 5 m7G cap. (B) Expression of PA, PB1, PB2, and NP supported GFP expression from the GP-GFP gene. Co-transfection of the pHW-CMV-GP-GFP with pcDNA-PA, PB1, PB2, and NP in HEK293T cells revealed green fluorescence under a UV microscope, but the co-transfection with only pcDNA-PB1, PB2, and NP did not. (C) PA, PB1, PB2, and NP expression supported gLuc expression from GP-gLuc gene. Co-transfection of the pHW-CMV-GP-gLuc with pcDNA-PA, PB1, PB2, and NP in HEK293T cells showed a high luciferase activity (green bar), but not with only three pcDNA plasmids or no plasmids (control). All luciferase activities (Luc) were detected by the Gaussia Luciferase Glow Assay Kit (Pierce). Individual values are represented by black solid dots and the columns are the means with standard deviations from eight experiments performed in duplicate. ****P < 0.0001 by one-way Student t-test.
St. Louis, MO, United States), and separated on a 10% SDS-polyacrylamide gel. Proteins were transferred onto a polyvinylidene difluoride membrane (#IPVH00010; Millipore, Bedford, MA, United States). The membrane was blocked and probed with primary and secondary antibodies sequentially. Signals were visualized by enhanced chemiluminescence, and the pre-stained protein ladder was imaged under bright light simultaneously under a Fuji LAS 4000 imaging system.
An anti-BRBV M protein polyclonal antibody was produced in rats by immunizing them with the BRBV M protein, which was expressed and purified in bacteria, following a protocol previously described (Sun et al., 2009). Anti-β-actin (#A5441) antibody was purchased from Sigma (St. Louis, MO, United States).

Statistical Analysis
Inhibitory concentration (IC 50 ), cytotoxic concentrations (CC 50 ), and Effective Concentration (EC 50 ) were determined with the means and standard deviations obtained from three independent experiments performed in duplicate. Calculations of IC 50 , CC 50 , and EC 50 and statistical analysis were done by using GraphPad Prism version 8.0. Error bars represent means and standard deviations (SD), and statistical significance (P value) was determined by using 1-way ANOVA analysis, followed by Tukey-Kramer post-test for comparison of three or more groups and unpaired (Student) t-test for comparison of two groups.

BRBV-KS Infection of HEK293T Cells
We propagated BRBV-KS, obtained from BEI Resources 1 , in Vero cells. At 7 days post-infection, when most of the cells appeared cytopathic, the media were collected and centrifugated. The supernatant was collected, aliquoted, and stored at −80 • C. The virus was titrated on Vero cells with a titer of 1.5 × 10 8 plaque 1 www.beiresources.org/ forming units (pfu)/ml ( Figure 1A). We next developed a RT-qPCR to quantify the vgc numbers in the BRBV stock. We found a high correlation of the vgc numbers with the pfu titrated (R 2 = 0.9875) (Figure 1B). In the following experiments, we used the RT-qPCR to quantify the virus.
We tested human HEK293T cells for infection of BRBV for comparison with infection in Vero cells. Virus growth was determined over the course of infection, and infected cells were collected at 6 days post-infection for Western blotting of matrix (M) protein expression. The results showed that BRBV infected HEK293T cells displayed similar growth kinetics as those in Vero cells (Figure 2A), which had a peak of virus production of approximately (∼) 10 9 vgc/ml in the media, as well as detection of the M protein at a size of ∼30 kDa ( Figure 2B).  Table 1). Each viral gene was reversely transcribed using each forward primer and PCR amplified using respective paired primers ( Table 1). The PCR cDNA fragments were first cloned into pcDNA3 and sequenced. Sequences of GP, PA, and PB2 genes were confirmed as identical to those originally deposited in GenBank; whereas NP, PB1, and M genes have 10, 7, and 1 nucleotide variations, respectively (Figure 3A), which result in amino acid variations in NP and PB1 (Figure 3B).

Cloning of the BRBV cDNA Fragments
We next used the set of pHW primers ( Table 1) that had the BsmBI sites to amplify the GP cDNA and cloned it in a BsmBIdigested pHW2000 bidirectional cloning vector (Hoffmann et al., , on BRBV RdRP activity was quantified using the GP-gLuc reporter assay at different concentrations as indicated. DMSO was used as a positive control of the reporter (Positive), and transfection with only pcDNA-PB1, PB2, and NP was used as a negative control of the reporter (Negative). At 3 days post-transfection, the activity of the expressed Gaussia luciferase was quantified using a Gaussia Luciferase Glow Assay Kit (Pierce). Results from each compound were normalized to the vehicle control in each experiment and were expressed with the means and standard deviations obtained from three experiments performed in duplicate. (C,D) Cytotoxicity assay. Compounds were assayed for cytotoxicity in HEK293T cells at the concentrations indicated using a CytoTox-Glo TM Cytotoxicity Assay kit (Promega). The results are normalized to the mock-treated control (set as 1.0). Values are the means and standard deviations obtained from three experiments, each performed in duplicate.
(E) IC 50 of compound #9 on RdRP activity. Compound #9 was applied to HEK293T cells transfected with the gLuc reporter system of five plasmids at various concentrations as indicated. The polymerase activity was quantified by measuring the luciferase activity. Values are the means with standard deviations and are normalized to the mock group from three experiments performed in duplicate. The half maximal inhibitory concentration (IC 50 ) was calculated with GraphPad Prism software. P values are calculated using one-way and Tukey-Kramer post-test, compared with DMSO or mock control. *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001; and n.s. (P > 0.05) denotes no statistical significance. 2002). The pHW2000-GP clone that had the right direction of BRBV cDNAs was confirmed by sequencing, in which there is the immediate early promoter of CMV, RNA polymerase II (pol II) promoter, a pol I terminator in the 3 UTR, a ribosome RNA pol I promoter and a bGHpA in the 5 UTR.

Establishment of a Viral RNA Dependent RNA Polymerase (RdRP) Activity Assay, a Replicon Reporter System of BRBV
Working on BRBV requires a BSL3 facility, thus, a replicon system of BRBV is important for assessing antivirals in a BSL2 setting and for high throughput screening of antivirals. Based on the pHW-CMV-GP plasmid, we inserted a GFP ORF and a secreted gLuc ORF between the 3 and 5 UTR of the GP segment (Figure 4A), which resulted in plasmids pHW-CMV-GP-GFP and pHW-CMV-gLuc. When pHW-CMV-GP-GFP was co-transfected with pcDNA-PA, PB1, PB2, and NP plasmids in HEK293T cells, at 7 days post-transfection, green fluorescence was clearly observed in cells transfected with pcDNA-PA, PB1, PB2, and NP plasmids but not in cells transfected with all plasmids except pcDNA-PA ( Figure 4B). Similarly, a high luciferase activity was detected in the cells cotransfected with pHW-CMV-GP-gLuc and pcDNA-PA, PB1, PB2, and NP plasmids, but nearly background levels were detected in cells transfected with pHW-CMV-GP-gLuc and pcDNA-PB1, PB2, and NP plasmids ( Figure 4C). Overall, there was a 10-fold increase in luciferase activity by the addition of the pcDNA-PA plasmid, suggesting a function of the RdRP activity from expression of PA, PB1, PB2, and NP.
Collectively, our results demonstrated that a parvoviral endonuclease activity inhibitor, compound #9 (myricetin), inhibits BRBV replication in cells through inhibition of the RdRP activity with a therapeutic index (TI = CC 50 /EC 50 ) of >50.

DISCUSSION
In this study, we established a GFP/luciferase expression-based replicon reporter system for BRBV, and demonstrated it is suitable for screening of antivirals in a BSL2 setting. We identified myricetin (compound #9) as a potent anti-BRBV drug. As BRBV is classified a BSL3 agent, the replicon system of BRBV is important for screening antivirals in a high throughput format in a BSL2 setting. Our replicon reporter system of BRBV contains 4 plasmids, one in cis plasmid (pHW-CMV-gLuc) and three in trans plasmids (pcDNA-PA, PB1, PB2, and NP). Transfection of all them into HEK293T cells results replication of the gLuc gene that is flanked with the 5 and 3 UTR of the BRBV GP genome. It will be ideal to integrate the pol I expression cassette (5 UTR-gLuc-3 UTR) and a CMV promoter expression cassette of PA-PB1-PB2-NP into one plasmid and use it in transfection, or to generate a cell line that hosts the replicon reporter, which is currently under construction.
In contrast to influenza viruses, thogotoviruses are transmitted mainly through tick vectors and thus are also called "tickborne viruses" (Haig et al., 1965;Anderson and Casals, 1973). Phylogenetically, BRBV is most closely related to dhori virus and its subtype, batken virus, which have been known to occur in regions throughout Africa, Asia and Europe (Lvov et al., 1974;Moore et al., 1975;Filipe et al., 1985;Hubalek and Rudolf, 2012). Dhori and batken viruses have been isolated from Hyalomma  (Waterhouse et al., 2018). NP forms a tetramer, as shown, which has four arms and a center to capture viral RNA and pass it to the RdRP. The region of NP carrying the α-helix that includes aa 67-71 (SIVMA) is enlarged from one of the monomers. ticks, and antibodies against dhori virus have been detected in camels, goats, horses, cattle, and humans (Lvov et al., 1974;Moore et al., 1975;Filipe et al., 1985;Hubalek and Rudolf, 2012). Interestingly, batken virus has also been isolated from several mosquito species (Lvov et al., 1974). BRBV replicates at a high level in a variety of invertebrate and vertebrate cell lines, including HEK293T cells, and produces relatively high titers in all mammalian cells that are compatible with the known susceptibility of the human host to BRBV infection and diseases . However, it is possible that the lack of the endogenous expression of RIG-I (Fitzgerald et al., 2017), a key cytosolic RNA sensor, in HEK293T cells, or virus-expressed protein(s) that interrupt the innate immunity response in human cells, facilitates BRBV replication in HEK293 cells or human cells.
The cell cultured BRBV-KS has a series of sequence variations in one α-helix of the NP protein, compared to the originally reported sequences (KP657749) by direct sequencing of the patient's samples ( Figure 3B). We modeled the structure of BRBV NP protein for the amino acids' electric charges and propensities to form an α-helix against the structure of influenza D virus (Donchet et al., 2019) using SWISS-MODEL (Waterhouse et al., 2018) (Figure 8). The modeling indicates that the variations at aa 67-71 did not change the charges of any amino acids, had minimal effects on the polarity of the variant sequences, and the sequence was still able to form an α-helix. Interestingly, the BRBV-STL strain, which was isolated from a fatal BRBV case in St. Louis, MO, United States in 2017 (Bricker et al., 2019), had the identical sequence at aa 67-71 in the NP protein (MK453525) as that of the cell cultured BRBV-KS.
Currently, there is no treatment or vaccine for BRBVinfection caused diseases. Mice (CD-1) were susceptible to BRBV infection as evidenced by seroconversion, but the infection did not cause disease symptoms or death of the animals . Recently, mice lacking IFN-α/β receptor expression (IFNαAR −/− ) have been used as an animal model to evaluate antivirals against BRBV (Bricker et al., 2019;Fuchs et al., 2019). Favipiravir, a pan inhibitor of RdRP (Furuta et al., 2005), was shown to protect IFNαAR −/− mice from lethal BRBV infection. The EC 50 of favipiravir against BRBV infection of Vero was 310 µM (Bricker et al., 2019). In animals, administration of 150 mg/kg of favipiravir twice daily at 3 days post-infection protected all animals from death of BRBV infection. However, it was speculated that favipiravir reached a concentration of 1.28 mM in serum, leading to the argument that it is not practical to use favipiravir to treat BRBV-infected patients. Notably myricetin, that was identified in our study, has a 15-fold lower EC 50 (20 µM) compared with the 310 µM of favipiravir in Vero cells, suggesting that myricetin may be a more potent inhibitor of BRBV infection than favipiravir. Due to the lack of an animal BSL3 facility that can host BRBV-infected IFNαAR −/− mice, we were not able to examine myricetin in animals.
Myricetin, a common plant-derived flavonoid, exhibits a wide range of activities, including antioxidant, anticancer, antidiabetic, and anti-inflammatory activities, as well as antimicrobial activities (Semwal et al., 2016). Myricetin was reported as a strong inhibitor of reverse transcriptase of retroviruses (Ono et al., 1990). It has also been shown to be active in inhibition of nsp13, a helicase of the severe acute respiratory syndrome (SARS) coronavirus helicase that targets the ATPase activity in vitro (IC 50 = 2.71 µM) (Keum and Jeong, 2012;Yu et al., 2012). We found that myricetin does not exert obvious cytotoxicity in Vero and HEK239T cells, as well as in normal breast epithelial cells (Keum and Jeong, 2012;Yu et al., 2012). We previously found that myricetin inhibited over 90% of the endonuclease activity of B19V NS1N (the endonuclease domain of NS1) in vitro at 10 µM (Xu et al., 2018), and thus speculate that myricetin may inhibit the endonuclease activity of BRBV PA, which warrants further investigation.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

AUTHOR CONTRIBUTIONS
JQ, JT, SH, and SG contributed to the conception of the study. SH, KN, XW, and FC performed the experiments. SH, KN, JW, and JQ analyzed the data. JQ, SH, and JT wrote the manuscript. All authors contributed to the article and approved the submitted version.

FUNDING
This study was supported by PHS grant AI144564 from the National Institute of Allergy and Infectious Diseases. The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.