Next Article in Journal
The Relationship between Acute and Chronic Pancreatitis with Pancreatic Adenocarcinoma: Review
Next Article in Special Issue
Ensemble Models for Tick Vectors: Standard Surveys Compared with Convenience Samples
Previous Article in Journal
Potential Combination Drug Therapy to Prevent Redox Stress and Mitophagy Dysregulation in Retinal Müller Cells under High Glucose Conditions: Implications for Diabetic Retinopathy
Previous Article in Special Issue
Introduction and Tribute to Charlie Calisher
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Diseases
Volume 9
Issue 4
10.3390/diseases9040092
diseases-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic Evaluation of the Genus Coltivirus Indicates Genetic Diversity among Colorado Tick Fever Virus Strains and Demarcation of a New Species

by
Holly R. Hughes
*,
Jason O. Velez
,
Kelly Fitzpatrick
,
Emily H. Davis
,
Brandy J. Russell
,
Amy J. Lambert
,
J. Erin Staples
and
Aaron C. Brault
Arboviral Diseases Branch, Division of Vector-Borne Diseases, Centers for Disease Control and Prevention, Fort Collins, CO 80521, USA
*
Author to whom correspondence should be addressed.
Diseases 2021, 9(4), 92; https://doi.org/10.3390/diseases9040092
Submission received: 24 November 2021 / Revised: 13 December 2021 / Accepted: 15 December 2021 / Published: 17 December 2021
(This article belongs to the Special Issue Recent Studies of Arthropod-, Bat- and Rodent-Borne Viruses: A Theme Issue in Honor of Professor Charles H. Calisher)
Browse Figures
Versions Notes

Abstract

:
The type species of the genus Coltivirus, Colorado tick fever virus (CTFV), was discovered in 1943 and is the most common tick-borne viral infection in the Western US. Despite its long history, very little is known about the molecular diversity of viruses classified within the species Colorado tick fever coltivirus. Previous studies have suggested genetic variants and potential serotypes of CTFV, but limited genetic sequence information is available for CTFV strains. To address this knowledge gap, we report herein the full-length genomes of five strains of CTFV, including Salmon River virus and California hare coltivirus (CTFV-Ca). The sequence from the full-length genome of Salmon River virus identified a high genetic identity to the CTFV prototype strain with >90% amino acid identity in all the segments except segment four, suggesting Salmon River virus is a strain of the species Colorado tick fever coltivirus. Additionally, analysis suggests that segment four has been associated with reassortment in at least one strain. The CTFV-Ca full-length genomic sequence was highly variable from the prototype CTFV in all the segments. The genome of CTFV-Ca was most similar to the Eyach virus, including similar segments six and seven. These data suggest that CTFV-Ca is not a strain of CTFV but a unique species. Additional sequence information of CTFV strains will improve the molecular surveillance tools and provide additional taxonomic resolution to this understudied virus.

1. Introduction

The family Reoviridae is comprised of two subfamilies, Spinareovirdae and Sedoreoviridae [1]. The genus Coltivirus is one of nine genera in the subfamily Spinareroviridae, with two tick-borne viruses, Colorado tick fever virus (CTFV) and Eyach virus, representing the first viruses to comprise the genus [2]. Viruses belonging in the genus Coltivirus have genomes made of 12 segments of double-stranded RNA, encoding for 13 viral proteins (VP) 1–12 [3]. Segment nine of coltiviruses produces two proteins, VP9 and VP9′, via a readthrough of the termination codon [4]. Colorado tick fever virus was isolated in 1943 from a febrile patient and has subsequently been found to be distributed in North America in the Rocky Mountain region between 4000 and 10,000 feet, where the vector, Dermacentor andersoni, is common [5,6]. Additional serotypes of CTFV have been described, including California hare coltivirus (referred to in the literature as CTFV-Ca) and Salmon River virus [7,8]. Eyach virus was isolated in Germany in 1976 from an Ixodes species tick [9] and was subsequently identified in France in 1981 [10]. In contrast to the febrile disease of CTFV, Eyach virus is known to cause neuroinvasive disease in humans. Recent discoveries of Kundal virus in India [11], Tarumizu tick virus in Japan [12], and Tai Forest reovirus in Cote d’Ivoire [13] have expanded the genus to include five viral species. Of these newly identified coltiviruses, Tai Forest reovirus has only been isolated from bats, and a tick vector has not been identified.
CTFV is the most commonly reported tick-borne viral disease in the Western US. Cases occur in the early spring and summer months when ticks are most active. CTFV causes a febrile illness, often with a biphasic fever and leukopenia. Up to 30% of the cases require hospitalization [14]. Life-threatening complications and deaths due to CTFV infections are rare, with fatalities associated with meningoencephalitis or disseminated intravascular coagulation. Other reported complications include meningoencephalitis, hepatitis, pericarditis, myocarditis, atypical pneumonia, disseminated intravascular coagulation, and epididymo-orchitis [14,15]. In general, the number of CTFV cases reported annually has decreased over time [14,16]. It is unclear if this reflects changes in the testing or reporting practices and certain states where the disease is found not having given it a notifiable condition. However, recent increases in case counts were seen in selected states, but it is unclear if the incidence is related to human activity or increased clinical awareness [16]. Although CTFV has been a recognized cause of human illness since the 1940s, very limited genetic information is available for CTFV, with only one complete genome sequence available for the exemplar CTFV isolate, Florio. Herein, we report the full genomes of two strains of CTFV, two strains of Salmon River virus, and one strain of CTFV-Ca. Increasing the available genetic sequences will improve our knowledge of the genetic diversity, support surveillance, and aid in optimizing molecular diagnostic testing.

2. Materials and Methods

All protocols followed manufacturers’ recommendations unless otherwise noted.

2.1. Viruses and RNA

Viruses sequenced in the study were provided by the Arboviral Diseases Branch, Division of Vector-Borne Diseases, Centers for Disease Control and Prevention’s Arbovirus Reference Collection (ARC). Isolate details are described in Table 1. RNA was extracted from frozen cell culture stocks (vero cell culture supernatant or suckling mouse brain homogenate, Table 1) using the QIAmp Viral RNA mini kit (Qiagen, Germantown, MD, USA). RNA was treated with 50% v/v dimethyl sulfoxide and incubated 1.5 h at 65 °C to denature double stranded RNA [17]. RNA was purified using an RNeasy cleanup kit (Qiagen).

2.2. Genome Sequencing and Assembly

The cDNA was generated from extracted RNA using random primers provided in the Ovation RNA-seq system V2 kit (Tecan, Redwood City, CA, USA). Sequencing libraries were prepared from the cDNA using the Ion Xpress plus fragment library kit (Thermo Fisher, Waltham, MD, USA), barcoded with Ion Xpress barcodes (Thermo Fisher), and quantified using the Ion library TaqMan Quantitation kit (Thermo Fisher). Sequencing templates were prepared using the Ion One Touch 2 system and Ion Hi-Q View OT2 kits (Thermo Fisher). Whole genome sequencing was performed on the Ion Torrent Personal Genomics Machine system using the Ion Hi-Q view sequencing kit and 318 v2 chips (Thermo Fisher), sequencing two libraries per chip.
Sequences (fastq) with quality phred Q > 20 were uploaded into CLC Genomic workbench v20 (Qiagen). Genomes were assembled via de novo assembly with a bubble size of the average read length (150 bp). The longest contigs with deepest coverage were submitted for BLAST analysis through NCBI. Consensus sequences for each segment were extracted from CLC genomic workbench and used for reference guided assembly in SeqMan NGen v17 (DNASTAR, Madison, WI, USA). Whole genomes were obtained for each virus with a minimum of 20× coverage, and a minimum average of 8.9 read depth per nucleotide. Putative open-reading frames were identified using EditSeq (DNASTAR). Terminal ends were determined for all sequences using 5′ RACE System 2.0 kits (Thermo Fisher) and 3′ poly(A) tailing kit (Ambion, Austin, TX, USA), followed by 3′ RACE kit (Thermo Fisher) using gene specific primers (IDT, Coralville, IA, USA). RACE DNA templates were cloned with TOPO-TA sequencing kits (Thermo Fisher) and sequenced by capillary sequencing on an ABI 3130 instrument (Thermo Fisher) using the primers provided in the TOPO-TA kits. All sequences generated in this study have been deposited into GenBank (Table 1).

2.3. Molecular and Phylogenetic Analysis

Conserved domains were identified in each segment using CD-search through NCBI (https://www.ncbi.nlm.nih.gov/Structure/cdd/wrpsb.cgi (accessed on 29 January 2021)). Nucleotide sequences comprising the open-reading frames for each protein were codon aligned using the Clustal W function of Mega 7 [18]. Percent sequence identity was calculated in Mega 7 using the p-distance model with G + I frequencies and complete deletion of missing data. Substitution models for each segment were determined using the model fit function of Mega 7. Bayesian inference was completed using BEAST v 1.8.4 [19] executed through the CIPRES Scientific Gateway (www.phylo.org (accessed on 1 February 2021)) [20]. A lognormal relaxed clock and a coalescent constant tree prior was used with an MCMC of 50 million generations and 10% burn-in. Priors were tested using the generalized stepping-stone sampling method [21]. Convergence of parameters was verified using TRACER v1.5 [22]. Maximum clade credibility trees were generated using TreeAnnotator [19] and FigTree V 1.4.3 showing the posterior probability of each branch. Maximum likelihood trees generated in Mega 7 were used to verify Bayesian tree topologies.
Reassortment was evaluated by concatenating the open-reading frames of 12 segments representing the complete genomes. Only reassortment events associated with a full segment of the concatenated genome that had p < 0.05 using the Bonferroni adjustment and detected by three or more models using Recombination Detection Program (RDP) v 4.9 [23] were considered significant [24].

2.4. Serological Evaluation

Positive control hyperimmune acidic fluids were diluted 1:5 in a total of 120 µL, followed by a serial 2-fold dilution with BA-1 media. Dilutions were then challenged with 100 plaque forming units of virus in a total of 60 µL, and incubated 1 hr at 37 °C. Following incubation, 100 µL of the virus with control mixture was added to a confluent Vero Cell monolayer on a 6-well plate (Corning, Corning, NY, USA) and incubated at 37 °C for 1 hr. Following incubation, 3 mL of yeast extract-lactalbumin (Ye-lah) overlay media with 0.5% (w/v) Agarose (Lonza, Basel, Switzerland) and 3% (v/v) bicarbonate (Gibco, Waltham, MD, USA) were added to each well, and plates were incubated at 37 °C. 3 mL of secondary overlay of Ye-lah media with 0.5% (w/v) agarose, 3% (v/v) bicarbonate, and 2% (v/v) neutral red (Sigma-Aldrich, St. Louis, MO, USA) was applied after an appropriate number of days to allow for differences in viral growth: 5 days for Eyach virus, 4 days for Salmon River, 4 days for CTVF-Ca, and 2 days for CTFV. Plaque counts were then noted on the following day. Titers were determined by taking the reciprocal of the highest serum dilution that reduced plaque counts by 90%, as indicated by the back titration of each virus.

3. Results

3.1. Phylogenetic Construction of the Genus Coltivirus

Phylogenies were constructed using Bayesian inference. The phylogenetic tree for segment one, the putative polymerase, places the CTFV and Salmon River viruses in a well-supported monophyletic clade, while CTFV-Ca is in a clade with Eyach virus, with the Kundal and Tarimuzu Tick viruses forming a third clade (Figure 1a). Similar phylogenetic trees were inferred for all the segments (Supplementary Figure S1), except segments four and six. Trees for segments four (Figure 1b) and six (Supplementary Figure S1) place the CTFV isolate R124719 basal to the node for the CTFV clade. Recombination detection supports the placement of segment four of R124719. Taken together, these data suggest segment four of isolate R124719 is the product of a reassortment event between CTFV and an unknown donor.

3.2. Segment Analysis

Each segment was evaluated for size, open reading frames, and evolutionary distances. All the CTFV and Salmon River viruses sequenced in this study have consistent segment and ORF sizes (Table 2). The California hare coltivirus and Eyach viruses differ in the segment length and ORF size in segments two and seven, having different lengths compared to CTFV, Salmon River virus, and each other. The ORF of segment two is 3804-nt in length for CTFV and Salmon river virus. However, segment two of the CTFV-Ca and Eyach viruses is 3858 and 3828-nt in length, respectively. The ORF in segment seven for CTFV and Salmon River virus is 1983-nt, while the ORF of segment six for CTFV-Ca is 2097-nt, and for Eyach virus is 2100-nt in length. Conversely, the segment six ORF of CTFV and Salmon River virus is 2094-nt in length, while the ORF of segment seven in CTFV-Ca and Eyach is also 2094-nt in length. Additionally, conserved domains were identified in segment seven of CTFV-Ca, indicating the presence of nucleotide binding proteins, such as chromosome segregation proteins in amino acids 311–575. These similar regions have previously been identified in Eyach segment seven and CTFV segment six [25]. These data suggest that CTFV-Ca has a similar genetic structure to Eyach virus. Interestingly, the ORF of segment 12 was found to be 558-nt in length in CTFV, Salmon River virus, and CTFV-Ca but 555-nt in Eyach virus. Although the ORF of segment 12 in CTFV-Ca is similar in size to CTFV, phylogenetic inference placed CTFV-Ca and Eyach in a well-supported clade (Supplementary Figure S1). Of note, segment four in CTFV isolate R124719 had a longer ORF than the viruses evaluated in this study. Surprisingly, the CTFV sequenced in this study differed in size between the reference CTFV, isolate Florio, in segments two and eleven. When comparing the reference sequence (RefSeq) from CTFV–Florio–RefSeq for segment two (NC_004182) to CTFV–Florio and CTFV R124719 sequenced in this study, the reference sequence contains one nucleotide insertion at position 850, one nucleotide deletion at position 886, and one nucleotide deletion at position 3631. All three identified discrepancies were attributed to a thymine deletion or insertion in the homopolymer regions. Additionally, the CTFV–Florio–RefSeq sequence displays a 4-nt deletion beginning at position 699 for segment 11 (NC_004191) compared to the CTFV strains sequenced in this study.
Both CTFV-Ca and Salmon River viruses have the canonical terminal ends attributed to coltiviruses. The most commonly observed terminal end sequence in both viruses was 5′-GACATTTG…TTGCAGTC-3′. Although subtle differences between the segments were observed with CTFV-Ca, the overall 3′ terminal end sequences were determined to be 5′-DTGCAGTC-3′.

3.3. Evolutionary Distances

Percent identities comparing nucleotides (Figure 2a) and amino acids (Figure 2b) were calculated for CTFV, Salmon River, CTFV-Ca, and Eyach viruses (Supplementary Table S1). Overall, percent identities between CTFV and Salmon River viruses are highest and ranged from 98% nucleotide identity in segment 12 to 75% in segment four. California hare coltivirus and Eyach displayed the lowest identity to CTFV and Salmon River virus, as low as 59% nucleotide identity in segment seven. Overall, amino acid percent identities were increased compared to nucleotide identities in nearly all the segments, except CTFV-Ca and Eyach exhibited lower amino acid than nucleotide identities in segments six (<59%), seven (<58%), and twelve (<71%) compared to the CTFV and Salmon River viruses. Interestingly, CTFV isolate R124719 had reduced identities in segment four compared to CTFV Florio and the Salmon River viruses. This low identity supports the segment four phylogenetic tree placing CTFV R124719 on a node basal to similar CTFV isolates and the observation of a longer ORF than other coltiviruses evaluated.

3.4. Serological Evaluation

Plaque reduction neutralization testing was completed for CTFV, Salmon River virus, and CTFV-Ca to understand the serological relationships (Table 3). CTFV Florio was neutralized efficiently by mouse hyperimmune acidic fluid raised against CTFV Florio (PRNT90 of 1:5,120) and was also cross-neutralized by Salmon River virus antiserum, although at a dilution >4-fold less than homologous neutralization (PRNT90 1:320). Salmon River virus was neutralized equally by antiserum generated against CTFV Florio and Salmon River virus (PRNT90 1:640). Salmon River virus was also cross-neutralized by antiserum produced against CTFV-Ca, but at a very low end-point titration (PRNT90 1:20). Moreover, CTFV-Ca was only neutralized by homologous antiserum (PRNT90 1:160) and CTFV-Ca antiserum did not cross-neutralize CTFV Florio or Salmon River virus. Eyach virus was cross-neutralized at very low dilutions by antiserum generated against Salmon River and CTFV Florio (PRNT90 1:10) but was not neutralized by antiserum from CTFV-Ca.

4. Discussion

The genus Coltivirus has expanded in recent years to include newly described viruses as new species (Kundal virus, Tai Forest reovirus, and Tarumizu Tick virus); however, very little full-length sequence information has been available for the known viruses in the CTFV species complex, specifically CTFV, Salmon River virus, and CTFV-Ca. This study was initiated to improve the current knowledge of CTFV diversity and expand the known database of sequence information. We were able to generate full-length sequence information for two isolates of CTFV, two isolates of Salmon River virus, and the only known isolate of CTFV-Ca. Collectively, the data presented herein highlight a previously underappreciated diversity of CTFV and suggest the placement of Salmon River virus as a strain of CTFV, while CTFV-Ca likely represents a unique coltivirus species.
California hare coltivirus was isolated in California in 1976 from a hare (Lepus californicus) and given the designation S6-14-03 [26]. Previous studies have determined partial sequence information for this virus [25], and, here, we describe the first complete genome sequences for all 12 segments of CTFV-Ca. Phylogenetic inference places CTFV-Ca consistently in a monophyletic clade with Eyach virus for all the segments. These tree topologies are supported by evolutionary distances as CTFV-Ca is more closely related to Eyach than CTFV. California hare coltivirus also had a similar genetic structure to Eyach virus as segment seven of these viruses have similar conserved functional domains and are most related to CTFV segment six [25]. The phylogenetic inference and genetic architecture of CTFV-Ca being more similar to Eyach virus seemingly supports the theory of the evolution of European coltiviruses from North American coltivirus following the migration of host species from North America through Asia [25,27]. In our serological analysis, CTFV-Ca does not display any cross-neutralization with any coltivirus tested, which is similar to previous serological analyses [26]. Collectively, these data suggest that CTFV-Ca is not a strain of CTFV but a unique species within the genus. This conclusion is not unique and supports previous serological studies suggesting CTFV-Ca is more closely related to Eyach but is a distinct species from both CTFV and Eyach [8]. However, applying taxonomic species demarcation >55%, >57%, >60% amino acid identity in VP6, VP7, and VP12, respectively, CTFV-Ca is considered a strain of both CTFV and Eyach virus. Given this newly described genetic information, the species demarcation criteria may need to be revised for the genus Coltivirus.
Salmon River virus was isolated from an ill person in Idaho in 1990. Presently, Salmon River virus is an unclassified coltivirus [1]. In our analysis, Salmon River virus formed well supported monophyletic clades with reference strains of CTFV and a recent Salmon River virus isolate. This analysis is supported by the high percent identities described between the CTFV isolates and the two Salmon River virus isolates. Moreover, some cross-neutralization was observed between the two viruses, with strong neutralization of the Salmon River virus by CTFV antiserum and 4-fold less neutralization of the CTFV by Salmon River virus antiserum. While previous studies have described Salmon River virus as serologically distinct from CTFV, genetic relatedness was demonstrated [7]. The full-length sequences of Salmon River virus generated in this study confirm these previous observations and display a high degree of sequence identity between Salmon River and CTFV [7]. Using the current species demarcation criteria for the genus Coltivirus, Salmon River virus could be considered a strain of the species Colorado tick fever Coltivirus.
The highest genetic diversity between the CTFV strains was observed within segment four. The nucleotide diversity in segment four is supported by the unique phylogenetic placement of the CTFV strains compared to other segments, as well as RDP analysis, suggesting this segment was subject to reassortment. Our observation of segment reassortment is supportive of previous studies suggesting the reassortment of multiple CTFV segments in the field [28,29]. A low sequence identity in segment four was found between CTFV isolates as well as other coltiviruses, with the nucleotide percent identity averaging 77%. The low level of nucleotide identity described in this study is supported by previous studies using RNA hybridization techniques, which also determined segment four to exhibit high genetic variation between CTFV isolates [30]. Although the sequence identity to segment four was low, the amino acid identities were higher, suggesting, although the RNA displays variation, the product of the gene is under strong purifying selection pressure. The high amino acid identity is likely due to the prediction that the product of segment four is a capsid protein and responsible for essential membrane binding [31]. Overall, the percent identities across all the segments with the CTFV strains were found to be <90%, with the exception of segments nine and twelve with very high identity, suggesting that these segments may be good targets for the molecular detection of both CTFV and Salmon River viruses [7].
Of unique interest, the sequences derived in this study display sequence variation from the CTFV strains previously sequenced and deposited into public databases. The sequences derived from CTFV Florio prototype virus obtained from the Arbovirus Reference Collection differed in segments two and eleven from the published reference sequences [3]. Prior to this study, only one whole genome reference was available for CTFV. While a previous study has partially sequenced segment eleven from multiple CTFV isolates [32], the amplified region published did not overlap with the four nucleotide deletion of segment eleven data observed in this study. Additionally, the published partial sequences of CTFV-Ca isolate S6-14-03 [25] match more closely to the CTFV and Salmon River viruses sequenced in this study than to the CTFV-Ca sequences presented here. The sequences of CTFV-Ca derived here show a higher degree of relatedness to Eyach virus than to the CTFV reference sequences and CTFV sequences derived in this study, which agrees with the serological and molecular data suggesting CTFV-Ca is distinct from CTFV [7,8]. In addition, the genetic stability of coltivirus strains when subjected to passaging remains unknown and could also explain the observed diversity. Overall, these data highlight the importance of the independent evaluation of reference strain sequences and depositing virus isolates in multiple reference collections to preserve historical value [33,34]. Additionally, these data highlight the importance of the periodic sequencing of newer strains and strains from other locations to ensure the validity of molecular targets for clinical diagnostic testing.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/diseases9040092/s1, Figure S1: Phylogenetic inference of the genus Coltivirus, Table S1: Segment percent identities.

Author Contributions

Conceptualization, H.R.H., A.J.L., J.E.S., and A.C.B.; methodology, H.R.H., A.J.L., J.E.S., and A.C.B.; data generation and validation, H.R.H., J.O.V., K.F., B.J.R., and E.H.D.; formal analysis, H.R.H., E.H.D., and A.J.L.; resources, A.C.B., A.J.L., and B.J.R.; writing—original draft preparation, H.R.H.; writing—review and editing, J.O.V., E.H.D., K.F., B.J.R., A.J.L., J.E.S., and A.C.B.; All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Sequences generated in this study have been deposited into GenBank (Table 1).

Acknowledgments

The findings and conclusions in this report are those of the authors and do not necessarily represent the official position of the Centers for Disease Control and Prevention. Mention of trade names is for specific information only and does not imply endorsement by the Centers for Disease Control and Prevention.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. King, A.M.; Lefkowitz, E.; Adams, M.J.; Carstens, E.B. Virus Taxonomy: Ninth Report of the International Committee on Taxonomy of Viruses; Elsevier: London, UK, 2011; Volume 9. [Google Scholar]
  2. Attoui, H.; Billoir, F.; Biagini, P.; de Micco, P.; de Lamballerie, X. Complete sequence determination and genetic analysis of Banna virus and Kadipiro virus: Proposal for assignment to a new genus (Seadornavirus) within the family Reoviridae. J. Gen. Virol. 2000, 81 Pt 6, 1507–1515. [Google Scholar] [CrossRef] [PubMed]
  3. Attoui, H.; Billoir, F.; Biagini, P.; Cantaloube, J.F.; de Chesse, R.; de Micco, P.; de Lamballerie, X. Sequence determination and analysis of the full-length genome of colorado tick fever virus, the type species of genus Coltivirus (Family Reoviridae). Biochem. Biophys. Res. Commun. 2000, 273, 1121–1125. [Google Scholar] [CrossRef] [PubMed]
  4. Jaafar, F.M.; Attoui, H.; de Micco, P.; de Lamballerie, X. Termination and read-through proteins encoded by genome segment 9 of Colorado tick fever virus. J. Gen. Virol. 2004, 85 Pt 8, 2237–2244. [Google Scholar] [CrossRef]
  5. Florio, L.; Miller, M.S.; Mugrage, E.R. Colorado tick fever; isolation of the virus from Dermacentar andersoni in Nature and a laboratory study of the transmission of the virus in the tick. J. Immunol. 1950, 64, 257–263. [Google Scholar] [PubMed]
  6. Florio, L.; Stewart, M.O.; Mugrage, E.R. The etiology of Colorado tick fever. J. Exp. Med. 1946, 83, 1–10. [Google Scholar] [CrossRef] [PubMed]
  7. Johnson, A.J.; Karabatsos, N.; Lanciotti, R.S. Detection of Colorado tick fever virus by using reverse transcriptase PCR and application of the technique in laboratory diagnosis. J. Clin. Microbiol. 1997, 35, 1203–1208. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Karabatsos, N.; Poland, J.D.; Emmons, R.W.; Mathews, J.H.; Calisher, C.H.; Wolff, K.L. Antigenic variants of Colorado tick fever virus. J. Gen. Virol. 1987, 68 Pt 5, 1463–1469. [Google Scholar] [CrossRef] [PubMed]
  9. Rehse-Kupper, B.; Casals, J.; Rehse, E.; Ackermann, R. Eyach--an arthropod-borne virus related to Colorado tick fever virus in the Federal Republic of Germany. Acta Virol. 1976, 20, 339–342. [Google Scholar] [PubMed]
  10. Chastel, C.; Main, A.J.; Couatarmanac’h, A.; Le Lay, G.; Knudson, D.L.; Quillien, M.C.; Beaucournu, J.C. Isolation of Eyach virus (Reoviridae, Colorado tick fever group) from Ixodes ricinus and I. ventalloi ticks in France. Arch. Virol. 1984, 82, 161–171. [Google Scholar] [CrossRef] [PubMed]
  11. Yadav, P.D.; Whitmer, S.L.M.; Sarkale, P.; Fei Fan Ng, T.; Goldsmith, C.S.; Nyayanit, D.A.; Esona, M.D.; Shrivastava-Ranjan, P.; Lakra, R.; Pardeshi, P.; et al. Characterization of Novel Reoviruses Wad Medani Virus (Orbivirus) and Kundal Virus (Coltivirus) Collected from Hyalomma anatolicum Ticks in India during Surveillance for Crimean Congo Hemorrhagic Fever. J. Virol. 2019, 93, e00106-19. [Google Scholar] [CrossRef] [Green Version]
  12. Fujita, R.; Ejiri, H.; Lim, C.K.; Noda, S.; Yamauchi, T.; Watanabe, M.; Kobayashi, D.; Takayama-Ito, M.; Murota, K.; Posadas-Herrera, G.; et al. Isolation and characterization of Tarumizu tick virus: A new coltivirus from Haemaphysalis flava ticks in Japan. Virus Res. 2017, 242, 131–140. [Google Scholar] [CrossRef] [PubMed]
  13. Weiss, S.; Dabrowski, P.W.; Kurth, A.; Leendertz, S.A.J.; Leendertz, F.H. A novel Coltivirus-related virus isolated from free-tailed bats from Cote d’Ivoire is able to infect human cells in vitro. Virol. J. 2017, 14, 181. [Google Scholar] [CrossRef]
  14. Yendell, S.J.; Fischer, M.; Staples, J.E. Colorado tick fever in the United States, 2002–2012. Vector-Borne Zoonotic Dis. 2015, 15, 311–316. [Google Scholar] [CrossRef] [PubMed]
  15. Stapels, J.E.; Fischer, M. Coltiviruses (Colorado Tick Fever). In Principles and Practice of Pediatric Infectious Diseases, 5th ed.; Long, S.S., Prober, G.C., Fischer, M., Eds.; Elsevier: Philadelphia, PA, USA, 2018; pp. 1119–1121. [Google Scholar]
  16. McDonald, E.; George, D.; Rekant, S.; Curren, E.; DeBess, E.; Hedberg, K.; Lutz, J.; Faith, J.; Kaisner, H.; Fawcett, R.; et al. Notes from the Field: Investigation of Colorado Tick Fever Virus Disease Cases—Oregon, 2018. Morb. Mortal. Wkly. Rep. 2019, 68, 289–290. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Wilcox, A.H.; Delwart, E.; Diaz-Munoz, S.L. Next-generation sequencing of dsRNA is greatly improved by treatment with the inexpensive denaturing reagent DMSO. Microb. Genom. 2019, 5, e000315. [Google Scholar] [CrossRef]
  18. Kumar, S.; Stecher, G.; Tamura, K. MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Drummond, A.J.; Suchard, M.A.; Xie, D.; Rambaut, A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol. Biol Evol. 2012, 29, 1969–1973. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Miller, M.A.; Pfeiffer, W.; Schwartz, T. Creating the CIPRES Science Gateway for Inference of Large Phylogenetic Trees. In Proceedings of the Gateway Computing Environments Workshop (GCE), New Orleans, LA, USA, 14 November 2010; pp. 1–8. [Google Scholar]
  21. Baele, G.; Lemey, P.; Suchard, M.A. Genealogical Working Distributions for Bayesian Model Testing with Phylogenetic Uncertainty. Syst. Biol. 2016, 65, 250–264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Drummond, A.J.; Rambaut, A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol. Biol. 2007, 7, 214. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Martin, D.P.; Murrell, B.; Golden, M.; Khoosal, A.; Muhire, B. RDP4: Detection and analysis of recombination patterns in virus genomes. Virus Evol. 2015, 1, vev003. [Google Scholar] [CrossRef] [Green Version]
  24. Maljkovic Berry, I.; Melendrez, M.C.; Li, T.; Hawksworth, A.W.; Brice, G.T.; Blair, P.J.; Halsey, E.S.; Williams, M.; Fernandez, S.; Yoon, I.K.; et al. Frequency of influenza H3N2 intra-subtype reassortment: Attributes and implications of reassortant spread. BMC Biol. 2016, 14, 117. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Attoui, H.; Mohd Jaafar, F.; Biagini, P.; Cantaloube, J.F.; de Micco, P.; Murphy, F.A.; de Lamballerie, X. Genus Coltivirus (family Reoviridae): Genomic and morphologic characterization of Old World and New World viruses. Arch. Virol. 2002, 147, 533–561. [Google Scholar] [CrossRef] [PubMed]
  26. Lane, R.S.; Emmons, R.W.; Devlin, V.; Dondero, D.V.; Nelson, B.C. Survey for evidence of Colorado tick fever virus outside of the known endemic area in California. Am. J. Trop. Med. Hyg. 1982, 31, 837–843. [Google Scholar] [CrossRef]
  27. Chastel, C. Erve and Eyach: Two viruses isolated in France, neuropathogenic for man and widely distributed in Western Europe. Bull. Acad. Natl. Med. 1998, 182, 801–809. [Google Scholar] [PubMed]
  28. Brown, S.E.; Miller, B.R.; McLean, R.G.; Knudson, D.L. Co-circulation of multiple Colorado tick fever virus genotypes. Am. J. Trop. Med. Hyg. 1989, 40, 94–101. [Google Scholar] [CrossRef]
  29. Williamson, B.N.; Fischer, R.J.; Lopez, J.E.; Ebihara, H.; Schwan, T.G. Prevalence and Strains of Colorado Tick Fever Virus in Rocky Mountain Wood Ticks in the Bitterroot Valley, Montana. Vector-Borne Zoonotic Dis. 2019, 19, 694–702. [Google Scholar] [CrossRef] [PubMed]
  30. Bodkin, D.K.; Knudson, D.L. Genetic relatedness of Colorado tick fever virus isolates by RNA-RNA blot hybridization. J. Gen. Virol. 1987, 68 Pt 4, 1199–1204. [Google Scholar] [CrossRef]
  31. Suzuki, N.; Supyani, S.; Maruyama, K.; Hillman, B.I. Complete genome sequence of Mycoreovirus-1/Cp9B21, a member of a novel genus within the family Reoviridae, isolated from the chestnut blight fungus Cryphonectria parasitica. J. Gen. Virol. 2004, 85 Pt 11, 3437–3448. [Google Scholar] [CrossRef]
  32. Attoui, H.; Charrel, R.N.; Billoir, F.; Cantaloube, J.F.; de Micco, P.; de Lamballerie, X. Comparative sequence analysis of American, European and Asian isolates of viruses in the genus Coltivirus. J. Gen. Virol. 1998, 79 Pt 10, 2481–2489. [Google Scholar] [CrossRef] [Green Version]
  33. Forshey, B.M.; Castillo, R.M.; Hang, J. Group C orthobunyavirus genomic sequences require validation. J. Virol. 2014, 88, 3052–3053. [Google Scholar] [CrossRef] [Green Version]
  34. Nunes, M.R.; Travassos da Rosa, A.P.; Weaver, S.C.; Tesh, R.B.; Vasconcelos, P.F. Reply to “Group C orthobunyavirus genomic sequences require validation”. J. Virol. 2014, 88, 3054. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Diseases 09 00092 g001 550
Figure 1. Bayesian Maximum Clade Credibility trees of Colorado tick fever virus strains and select members of the genus Coltivirus. Nucleotide coding sequences depicting phylogeny of (a) segment 1 and (b) segment 4. Viruses sequenced in this study are labeled with an asterisk. Taxa are labeled with virus name, strain, and GenBank accession number. Branches are labeled with the posterior probabilities, and scale bar depicts nucleotide substitutions per site. Colorado tick fever virus (CTFV).
Figure 1. Bayesian Maximum Clade Credibility trees of Colorado tick fever virus strains and select members of the genus Coltivirus. Nucleotide coding sequences depicting phylogeny of (a) segment 1 and (b) segment 4. Viruses sequenced in this study are labeled with an asterisk. Taxa are labeled with virus name, strain, and GenBank accession number. Branches are labeled with the posterior probabilities, and scale bar depicts nucleotide substitutions per site. Colorado tick fever virus (CTFV).
Diseases 09 00092 g001
Diseases 09 00092 g002 550
Figure 2. Comparison of genetic identities between Colorado tick fever virus strains and Eyach virus. Segment open-reading frames were codon aligned and percent identities were calculated for (a) nucleotide and (b) amino acids. Graphed points indicate median and range. CTF = Colorado tick fever.
Figure 2. Comparison of genetic identities between Colorado tick fever virus strains and Eyach virus. Segment open-reading frames were codon aligned and percent identities were calculated for (a) nucleotide and (b) amino acids. Graphed points indicate median and range. CTF = Colorado tick fever.
Diseases 09 00092 g002
Table
Table 1. Colorado tick fever viruses sequenced in this study obtained from the Arbovirus Reference Collection.
Table 1. Colorado tick fever viruses sequenced in this study obtained from the Arbovirus Reference Collection.
VirusIsolate DesignationDateLocationPassageGenBank No.
Colorado tick feverFlorio1943ColoradoSM14MZ484422-33
Colorado tick feverR1247192017WyomingVero 1MZ484434-45
Salmon River720896-11990IdahoSM3MW67046-57
Salmon RiverR1305442019OregonVero 1MZ484446-57
California hare coltivirusS6-14-031976CaliforniaSM5MW396980-91
SM = suckling mouse brain.
Table
Table 2. Segment comparison of Colorado tick fever viruses sequenced in this study and reference coltiviruses.
Table 2. Segment comparison of Colorado tick fever viruses sequenced in this study and reference coltiviruses.
Virus Nucleotide Open Reading Frame Size (Segment Size)
SegmentCTFV-Florio-RefSeq 1CTFV-Florio-ARC 2CTFV R124719Salmon RiverSalmon River R130544California Hare ColtivirusEyach
14308 (4350)4308 (4350)4308 (4350)4308 (4349)4308 (4349)4308 (4349)4308 (4349)
23630 4 (3909)3804 (3909)3804 (3909)3804 (3908)3804 (3908)3858 (3964)3828 (3934)
33549 (3586)3549 (3586)3549 (3586)3549 (3585)3549 (3585)3549 (3585)3549 (3585)
43084 (3157)3084 (3156)3087 (3159)3084 (3156)3084 (3156)3084 (3156)3084 (3156)
52256 (2432)2256 (2432)2256 (2432)2256 (2431)2256 (2431)2256 (2419)2256 (2398)
6 32094 (2141)2094 (2141)2094 (2141)2094 (2140)2094 (2140)2097 (2175)2100 (2178)
72055 (2133)2055 (2133)2055 (2133)2055 (2133)2055 (2132)2094 (2139)2094 (2139)
81983 (2029)1983 (2029)1983 (2029)1983 (2028)1983 (2028)1983 (2029)1983 (2028)
91809 (1884)1809 (1884)1809 (1884)1809 (1884)1809 (1884)1809 (1884)1809 (1884)
101818 (1880)1818 (1880)1818 (1880)1818 (1880)1818 (1880)1818 (1879)1818 (1879)
11750 (998)927 (1002)927 (1002)927 (1002)927 (1002)927 (1002)927 (1002)
12558 (675)558 (675)558 (676)558 (675)558 (675)558 (678)555 (678)
1 Colorado tick fever virus, Florio strain, sequences obtained from the reference sequence database (NC_004180-90). 2 Colorado tick fever virus, Florio strain, obtained from the Arbovirus Reference Collection (ARC). 3 Bold segments highlight the similarity between segments 6 of CTFV and 7 of Eyach and California hare coltivirus. 4 Underlined values represent differences in size between reference strains and viruses sequenced in this study. CTFV = Colorado tick fever virus.
Table
Table 3. Serological evaluation of Colorado tick fever virus strains.
Table 3. Serological evaluation of Colorado tick fever virus strains.
Plaque Reduction Neutralization Test 90%
Hyperimmune Fluid AntibodySalmon River Virus 1CTFV-FlorioCalifornia Hare ColtivirusEyach Virus 3
Salmon River6402320<1010
CTFV Florio6405120<1010
California hare coltivirus20<10160<10
1 Bold text shows virus neutralization titers with homologous antibody. 2 Reciprocal endpoint neutralization titers. <10 is negative neutralization. 3 Homologous hyperimmune fluid was unavailable for Eyach virus. CTFV = Colorado tick fever virus.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hughes, H.R.; Velez, J.O.; Fitzpatrick, K.; Davis, E.H.; Russell, B.J.; Lambert, A.J.; Staples, J.E.; Brault, A.C. Genomic Evaluation of the Genus Coltivirus Indicates Genetic Diversity among Colorado Tick Fever Virus Strains and Demarcation of a New Species. Diseases 2021, 9, 92. https://doi.org/10.3390/diseases9040092

AMA Style

Hughes HR, Velez JO, Fitzpatrick K, Davis EH, Russell BJ, Lambert AJ, Staples JE, Brault AC. Genomic Evaluation of the Genus Coltivirus Indicates Genetic Diversity among Colorado Tick Fever Virus Strains and Demarcation of a New Species. Diseases. 2021; 9(4):92. https://doi.org/10.3390/diseases9040092

Chicago/Turabian Style

Hughes, Holly R., Jason O. Velez, Kelly Fitzpatrick, Emily H. Davis, Brandy J. Russell, Amy J. Lambert, J. Erin Staples, and Aaron C. Brault. 2021. "Genomic Evaluation of the Genus Coltivirus Indicates Genetic Diversity among Colorado Tick Fever Virus Strains and Demarcation of a New Species" Diseases 9, no. 4: 92. https://doi.org/10.3390/diseases9040092

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop