Genetic Diversity of Toscana Virus

Distribution of Toscana virus (TOSV) is evolving with climate change, and pathogenicity may be higher in nonexposed populations outside areas of current prevalence (Mediterranean Basin). To characterize genetic diversity of TOSV, we determined the coding sequences of isolates from Spain and France. TOSV is more diverse than other well-studied phleboviruses (e.g.,Rift Valley fever virus).

Distribution of Toscana virus (TOSV) is evolving with climate change, and pathogenicity may be higher in nonexposed populations outside areas of current prevalence (Mediterranean Basin). To characterize genetic diversity of TOSV, we determined the coding sequences of isolates from Spain and France. TOSV is more diverse than other well-studied phleboviruses (e.g.,Rift Valley fever virus).
T oscana virus (TOSV) belongs to the family Bunyaviridae and genus Phlebovirus. According to the Eighth Report of the International Committee on Taxonomy of Viruses, phleboviruses are classifi ed into 9 serocomplexes (1), 1 of which includes Rift Valley fever virus (RVFV), a deadly pathogen for cattle and humans. TOSV belongs to the sandfl y virus serotype Naples serocomplex. Phleboviruses have segmented RNA genomes comprised of 3 units: L (large), M (medium), and S (small). Based on the analysis of the G N glycoprotein (M segment), 4 lineages of TOSV have been proposed (2). Moreover, phylogenetic analysis has demonstrated that TOSV isolates from Spain differ from the original isolates from Italy, TOSV strain ISS Phl. 3 (3).
TOSV is widespread in the Mediterranean Basin, and evidence of human infection has been found in Italy, France, Spain, Portugal, Cyprus, and Turkey. The main clinical manifestation is neurologic dysfunction (4). Although the virus is an important pathogen (4), little genomic information is available for TOSV. We believe that obtaining genome sequence information for viruses with poor representation in public databases is an urgent task for the virologic community, especially for programs of virus surveillance and study of emerging pathogens. The increas-ing popularity of nucleic acid-based methods of detection underscores the need for a deeper knowledge of sequence variability of wild-type strains to ensure the sensitivity and specifi city of the diagnostic assays. We therefore characterized the genetic diversity of TOSV by determining the coding sequences of isolates from Spain and France.

The Study
Viruses were recovered from Vero E6 (ATCC CRL-1586) cell cultures of cerebrospinal fl uid samples (CSF) from 2 TOSV-infected patients and from a pool of sandfl ies. The fi rst patient was a 30-year-old man living in a rural area of Granada, Spain. He had been admitted to hospital with headache, nausea, and fever; he was found to have meningeal and CSF pleocytosis. TOSV ESH 62100 was detected in an acute-phase CSF sample by reverse transcription-PCR (RT-PCR) (5). The second patient was a 57-year-old woman from Augbane (Marseille area), France, who had fever, malaise, photophobia, neck rigidity, and vomiting. TOSV-specifi c immunoglobulin M was detected in an acute-phase serum sample; a convalescent-phase sample indicated seroconversion. Her CSF sample was positive for TOSV by PCR (6) and was the source for the isolation of TOSV H/IMTSSA. TOSV EsPhGR40 was isolated from a pool of female sandfl ies captured in the metropolitan area of Granada (3).
Virus sequences were obtained by using conserved primers (primer sequences available on request). Gen-Bank accession numbers of the new TOSV sequences were EF120631 and FJ153279-FJ153286. To assess the potential for establishing a simple cutoff for classifi cation of TOSV genotypes (e.g., similar to programs created for mumps) (7), we used pairwise sequence comparison to compare these sequences with all published phleboviral sequences. Sequence analysis of the M segment was sufficient to enable determination of genotypes ( Figure 1). Two different clusters were clearly distinguished. We propose to name those clusters genotypes A and B. The previously described TOSV lineages that include isolates from Italy all clustered in the proposed TOSV A genotype. Comparision with RVFV, for which 7 lineages (A to G) have been described (8), showed higher divergence in L and M segments and nonstructural (NS) and nucleoprotein (N) genes between the proposed TOSV genotypes than that found among RVFV lineages (Table).
Pairwise comparison demonstrated the following: 1) all TOSV S segments were highly conserved; 2) the L segment demonstrated less conservation than the N gene at the nucleotide or deduced amino acid levels; 3) TOSV M segments were the most divergent; and 4) variation in M segments was higher than that among RVFV strains but less pronounced than that within the group of sandfl y fever Sicilian viruses or sandfl y fever Naples viruses (Table). With regard to the conserved S segments, the NS gene was more variable than the N gene, and despite a nucleotide pairwise difference of 12.1% in the N coding region, the amino acid sequences were completely conserved (100% identity); Phylogenetic analysis of all phlebovirus sequences was performed by using the maximum-likelihood method available in PAUP* under tree bisection and reconnection branch swapping; the best-fi t model of nucleotide substitu-tion in each instance was determined by using MODELT-EST (9) (parameter values available upon request) ( Figure  2). Time to most recent common ancestor (TMRCA) was estimated for the M segment by using the data from 32 dated samples collected over 35 years and the Bayesian Markov Chain Monte Carlo approach (BEAST package; 10). We applied a relaxed molecular clock with an uncorrelated exponential distribution of rates, a general time reversible + I + F4 model of nucleotide substitution and logistic population growth. Only the M segment had enough sequence representation to perform the analysis. Nonetheless, because glycoproteins are present on the surface of virions, they are the proteins most exposed to the selective pressures of the host.
The estimated mean rate of evolutionary change was 9.1 × 10 -5 substitutions/site/year (95% highest probability density [HPD] = 2.5 × 10 -4 to 2.7 × 10 -6 substitutions/site/ year). Under this rate the mean TMRCA was estimated to be 3,265 years, although with wide variance (178-11,000 years). Studies of RVFV, a closely related phlebovirus (11) have estimated the mean RVFV M segment rate to be 2.42 × 10 -4 (95% HPD = 1.8 × 10 -4 to 3.0 × 10 -4 substitutions/site/year) (8). TOSV M segment TMRCA is higher than RVFV M segment TMRCA (117.3 years [variance 95-143 years]) (8). Although the wide variance of the TMRCA calculations might be affected by the scarce sequence information available, the differences observed might also be related to their biological differences. Both are phleboviruses transmitted by the bite of arthropods that could also be implied as amplifi er hosts; however, because TOSV is believed not to be amplifi ed on its  mammalian hosts, it seems that the RVFV cycle includes mammals. The high conservation of RVFV genome sequences has been interpreted as indicating that the overall tolerance for mutation within the RVFV genome is low or that the viruses in the group have a relatively recent common ancestor (12). More studies are needed to understand the signifi cance of the high degree of purifying selection observed in TOSV.

Conclusions
Our main goal, to improve the knowledge of sequence information for this neglected genus, was achieved by the   (13); however, the cocirculation of both genotypes in Portugal is confusing because sequences obtained independently in 2 laboratories from strains allegedly obtained from the same patient clustered in different clades (DQ479890 and AY129737) (14,15). Geographic differences in genotype distribution may relate to differences in vector distribution. The distribution of phlebotomines is evolving with climate change and has implications for the epidemiology of vector-borne infectious diseases. Preexisting immunity almost certainly plays a role in limiting illness in the Mediterranean Basin; pathogenicity may be higher in naive populations outside areas of current TOSV prevalence. Thus, TOSV might be considered a threat for human health, and research and surveillance programs should aim to prevent its spread to new areas.