Next Article in Journal
Comprehensive Profiling of EBV Gene Expression and Promoter Methylation Reveals Latency II Viral Infection and Sporadic Abortive Lytic Activation in Peripheral T-Cell Lymphomas
Next Article in Special Issue
The Role of the Leishmania infantum Infected Dogs as a Potential Reservoir Host for Toscana Virus in a Zoonotic Visceral Leishmaniasis Focus of Northern Tunisia
Previous Article in Journal
Temporal Series Analysis of Population Cycle Threshold Counts as a Predictor of Surge in Cases and Hospitalizations during the SARS-CoV-2 Pandemic
Previous Article in Special Issue
Monitoring of Insecticide Resistance Mutations and Pathogen Circulation in Sand Flies from Emilia-Romagna, a Leishmaniasis Endemic Region of Northern Italy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

New Isolation of Ponticelli III Virus (Bunyavirales: Phenuiviridae) in Emilia-Romagna Region, Italy

1
Microbial Biotechnologies, Agrosciences and Environment Laboratory (BioMAgE), Faculty of Sciences Semlalia, Cadi Ayyad University, Marrakesh 40000, Morocco
2
Istituto Zooprofilattico Sperimentale della Lombardia e dell’Emilia-Romagna “Bruno Ubertini”, 42124 Brescia, Italy
3
Higher Institute of Nursing and Technical Health Occupations (ISPITS), Ministry of Health, Rabat 10000, Morocco
*
Author to whom correspondence should be addressed.
Viruses 2023, 15(2), 422; https://doi.org/10.3390/v15020422
Submission received: 13 October 2022 / Revised: 16 January 2023 / Accepted: 29 January 2023 / Published: 2 February 2023
(This article belongs to the Special Issue Sand Fly-Borne Phleboviruses, Volume II)

Abstract

:
The number of newly described sandfly-borne phleboviruses has been steadily growing in recent years. Some phleboviruses are human pathogens, but their health relevance is largely uncharacterized. We aimed to investigate the circulation of these viruses in the Emilia-Romagna region where several have already been described. A total of 482 sandflies were collected in a site in Reggio Emilia in 2019 and 2020. Sandflies collected in 2020 were grouped in 21 pools with a maximum of 25 sandflies per pool, submitted to real time PCR, and isolated in Vero cell culture. Complete genome sequencing showed the isolation of a strain of a Ponticelli III virus. This virus, which belongs to the species Adana phlebovirus, differed in the M segment from the Ponticelli I and Ponticelli II viruses. Analysis performed on the genomic segments of the newly isolated virus compared with other phleboviruses highlighted a strong purifying selection in the L segments, and different substitution saturation, highest in the M segments. Future research should address the ecological processes driving the occurrence of these novel phleboviruses and their possible impact on public health.

1. Introduction

The family Phenuiviridae includes 20 genera, the members of three of which are largely transmitted by arthropods: Phlebovirus, Uukuvirus and Ixovirus [1,2]. The genus Phlebovirus comprises about 70 named viruses widely distributed in Europe, Africa, Central Asia and the Americas. As provided by the International Committee on Taxonomy of Viruses (ICTV), these viruses are currently classified into species according to the threshold of 95% identity in the amino acid sequence of the RNA-dependent RNA polymerase (RdRp) protein [3]. Some of these viruses are medically relevant, such as the Rift Valley virus (RVFV) transmitted by mosquitoes and the Toscana virus (TOSV). Other medically relevant viruses transmitted by sandflies include sandfly fever Naples and Sicilian phleboviruses in the Old World, and Alenquer, Candiru, Chagres, Cocle, Echarate, Maldonado, and Punta Toro phleboviruses in the New World [4].
Phleboviruses are single-stranded negative-sense RNA with a genomic organization consisting of three segments, designated S (small), M (medium), and L (large). The S segment bears two genes encoding, by an ambisense strategy, the nucleoprotein (Nc) and non-structural (Ns) protein; the M segment encodes for enveloping the glycoproteins; and the L segment for the RdRp [5,6]. Phleboviruses are morphologically similar, with a diameter of 80-120 nm, icosahedral symmetry, and an envelope formed by short, spike-like glycoproteins embedded in a lipid bilayer [7,8].
In the Old World, at least 250 million people have been exposed to Phlebovirus infections [9,10]. One of the main examples is TOSV, which is known as an emerging pathogen causing meningitis in the Mediterranean region, for which defined reservoirs have not yet been identified. It is unlikely that humans are the reservoir for TOSV because human viremia duration is too short [11].
Sandfly-borne phleboviruses’ distribution and circulation peaks are closely related to vector presence and activity [12]. Sandflies are broadly distributed in regions bordering the Mediterranean basin. They are often abundant in peri-urban, rural and natural environments, which are concentrated near human and domestic animal habitats in suitable areas [13,14].
Recently, several studies conducted in Mediterranean countries showed an increase in Phlebovirus circulation with newly discovered viruses, many of which were detected in sandflies [15,16]. One example is the Adana virus, found in the eponymous name of Adana province in Turkey. The DNA barcoding characterization of the positive pool showed that the pool from which the virus was isolated contained Phlebotomus (Ph.) tobbi, Ph. perfiliewi, and Ph. papatasi [17]. The genetic and phylogenetic analysis revealed that Adana virus is closely related to Salehabad virus and Arbia virus [17]. Conspecific viruses of Adana virus were successively isolated in Italy (Ponticelli I, Ponticelli II, Ponticelli III viruses) [18], in Croatia, and in Northern Macedonia (Bregalaka virus), however they showed different M segments with respect to Adana virus [19].
In our current study, we report a new isolation of a Ponticelli III virus (PONTV III) strain in Emilia-Romagna at about 80 km from the site of the first report of this virus in the region [18].

2. Materials and Methods

2.1. Sandf Fly Trapping and Identification

Samples were collected in June 2019 and September 2020 in Casalgrande locality (44°33′40.3″ N, 10°43′27.8″ E), Reggio Emilia province. The site was located in a hilly area of the Emilia-Romagna region at 257 m altitude (Table 1). Sandfly sampling was performed overnight using two types of traps: carbon-dioxide-baited traps and CDC light traps. In 2019, sandflies were morphologically identified with an optical microscope after clarification, while in 2020, the collected sandflies were stored in dry ice and used directly for virus isolation in cell culture with a maximum of 25 specimens per pool.

2.2. Virus Isolation on Cell Cultures

Female specimens pools were homogenized in minimal essential medium, supplemented with penicillin and streptomycin, using a sterilized glass potter and clarified by centrifugation at 3000× g for 15 min. Samples were inoculated in a confluent monolayer of VERO cells (African green monkey kidney cells) at passage 170 (cell culture bio bank of IZSLER, code BSCL86), incubated at 37 °C with 5% CO2 and observed daily for 7 days to observe the development of cytopathic effects (CPE) [18]. To exclude the presence of two or more viral strains in culture, selected isolates were duplicated by plaque purification in two rounds of serial dilutions in 96-well plates (maximum dilution of 10−8). One of the obtained duplicates was recultured on VERO cells and subjected to sequencing.

2.3. Virus Detection and Sequencing

Viral RNA was extracted from supernatant cultures that displayed CPE using TRIzol LS Reagent (Invitrogen, Carlsbad, CA, USA); cDNA synthesis was performed using random hexamers (Roche Diagnostics, Mannheim, Germany) and SuperScript™ II reverse transcriptase (Invitrogen, Carlsbad, CA, USA). Complementary DNA was submitted to PCR for the detection of phleboviruses (Pan-Phlebo PCR) using one reverse and two forward primers and targeting a 370-nucleotide region of the S segment [20]. The obtained amplicons were sequenced, and the sequences were utilized for virus identification by BLAST analysis against the GenBank (GB) database (https://blast.ncbi.nlm.nih.gov/Blast.cgi, accessed on 28 January 2023).
The selected viruses isolated from the cell culture supernatant were sequenced on a MiSeq Instrument (Illumina Inc., San Diego, CA, USA). The preparation DNA libraries for whole-genome sequencing and viral RNA amplification were performed following the Kampmann et al. (2011) protocol [21].

2.4. Phylogenetic Analysis

The complete genome of the segments was aligned with the homologous sequences of other phleboviruses available in GB. The nucleotide alignments were obtained using the translated amino acidic sequence with the Mega11 software [22] (muscle algorithm with default parameters) and refined manually. Obtained alignments were used to infer the maximum likelihood trees by IQ tree software, with the selection of the best fit model performed by the software and branch supports inferred with ultrafast bootstrap analysis [23]. Software RDP5 [24] was used to check possible recombination events between different viruses. High variability was recorded in the alignments, so to minimize the rate of false positive result we used the methods MAXCHI, GENECONV and CHIMAERA, which showed most optimal performance in this situation. Only potential recombination events detected by all three methods were considered as suitable [25]. Aligned nucleotide sequences of the different segments of the PONTV were analyzed with DAMBE software to obtain a saturation plot in which transitions and transversions were plotted against GTR evolutionary distance. This allowed us to evaluate the genetic saturation of the different segments [26]. The number of transitions and transversions versus divergence offered a visual display of substitution saturation. To evaluate per-site selection in the same nucleotide alignment, the FEL (fixed effects likelihood) method of data monkey (using a p value of 0.05) was used, FEL infers nonsynoymous (dN) and synonymous (dS) substitution rates using a maximum-likelihood (ML) approach for a given coding alignment [27]. Saturation plot and per-site selection of the two ambisense genes of the S segment were conducted by concatenating the two sequences in the same reading frame. The p-distances (the complement number of identity) of amino acid and nucleotide sequence) were computed with MEGA X software [28]. Graphs were visualized by the ggplot2 package in R.

3. Results

A total of 482 sandflies were collected (73 in 2019 and 409 in 2020) in Casalgrande locality, Emilia-Romagna region, Italy. Only some specimens collected in 2019 were identified using morphological characteristics; of these, 72 were Ph. perfiliewi and 1 was Ph. perniciosus (Table 1). The 135 collected sandflies in 2020 were submitted to viral isolation. Of the five 25-specimen pools and one 10-specimen pool of sandflies, only one pool showed a CPE on VERO cells (Table 1). CPE observed on VERO cells, from day 3 post-infection, was marked by diffuse degeneration of the cell monolayer, with progressive rounding and subsequent detachment of the cells which appeared numerous and fluctuating in the culture medium. Cell culture supernatants with CPE were subjected to Pan-Phlebo PCR, producing a sequence ascribable to a PONTV (Figure 1). This sample was subjected to complete genome sequencing.
According to the previously obtained sequences, and particularly the M sequence, the new isolate was classified as PONTV III, strain 270315. Obtained sequences were deposited in GB by the accession numbers OP293791 for the L segment, OP293792 for the M segment, and OP293793 for the S segment. These sequences showed high identity for the PONTV II 238134-4/2016 (MG911981) in the L segment (99.13%), PONTV III 238134-3/2016 (MG911977) in the M segment (99.46%), and PONTV III 195684-2/2016 (MG911986) for the S segment (98.33%).
Sequences of the three segments were used with homologous sequences of sandfly- borne phleboviruses from the Mediterranean and Middle East that were available in GB for phylogenetic analysis. Two alignments were obtained with L and M segments (Figure 1) and one with the Ns and Nc (Figure S1). The evolutionary model selected for inferring ML trees, according to the lowest Bayesian information criterion, was the transition model with empirical base frequencies and free rate heterogeneity (TIM2+F+R4) for the L alignment, the general time reversible model with empirical base frequencies and free rate heterogeneity (GTR+F+R4) for the M alignment, the transition model with empirical base frequencies and gamma rate heterogeneity (TIM2+F+G4) for the Nc alignment, and the TPM2u with empirical base frequencies and gamma rate heterogeneity (TPM2u+F+G4) for the Ns alignment. In the L, Ns and Nc trees, the Ponticelli viruses were in a well-supported clade, in the M tree the three viruses were in separate branches (Figure 1 and Figure S1).
None of the putative recombination events were detected by all three methods selected by RDP5 software in the alignments obtained by sequences of L and M segments.
The phylogenetic analysis obtained of L segments (Figure 1A) and the p-distance (Table 2) confirmed that our PONTV III virus, as for other PONTVs, belonged to Adana phlebovirus species, according to the 95% identity threshold. The phylogenetic tree obtained for M segments (Figure 1B) showed some topologic incongruences with respect to the L, Nc and Ns segment trees (Figure 1A and Figure S1), most likely due to reassortment events. Interestingly, among the viruses of the Salehabad complex employed in the analysis, the Bregalaka virus had an M segment strictly related to PONTV II (Table 2, Figure 1B). The overall mean nucleotide p-distance between the different Ponticelli viruses was 0.017 for the S segments, 0.020 for the L segments, and rose to 0.21 for the different M segments, while the p-value between the same type of M segments ranged between 0.015 (PONTV I, PONTV III) and 0.012 (PONTV II).
The FEL analysis highlighted a pervasive purifying selection in all sites: the 83% non-invariant sites under purifying selection in the L segment (1388 on 1673 non invariant sites tested), 63% in the S segment (293 on 464 non invariant sites tested), and 73% in the M segment (947 on 1295 non invariant sites tested). The FEL analysis identified only one of the sites under diversifying positive selection in segment M. Plotting nucleotide p-distance against aminoacid distance, we observed that the different segments seemed subject to a differential dynamic of evolution, with an increasing amino acidic and nucleotide diversity in M segments, while L segments showed less aminoacid diversity for the same amount of nucleotide diversity to the other segments (Figure 2). The analysis performed on the three PONTV demonstrated that purifying selection was more intense in the L segment, as demonstrated by the FEL analysis and the lowest aminoacid diversity at equal nucleotide diversity of the L segment, respectively, to M and S segments.
The saturation plot showed a higher relative number of transversions in the M segment, which exceeded that of transitions with the increase in GTR distance. This situation strongly suggested a substation saturation in the M segments (Figure 3).

4. Discussion

Detection of phleboviruses in hosts is made difficult by the overall diversity and genetic complexity of clinically relevant strains, their essentially undefined clinical associations, and the resulting lack of diagnostic tools [4]. A paradigmatic example is the Toscana virus (TOSV), the etiological agent of neuroinvasive diseases including meningitis and meningoencephalitis. This virus was first isolated in 1971 by Verani et al. (1982) [29] and only later recognized as an agent of meningitis in humans [30]. In addition to TOSV, Sicilian phlebovirus (SFSV), and Naples phlebovirus (SFNV) also have pathogenic potential for humans, which is generally mild fever [31]. Different species of vertebrates are presumed to be hosts for these viruses, but convincing evidence is lacking. Antibodies to TOSV and SFSV have been found in horses, cats, dogs, sheep, pigs, cattle, and goats [9,31]. TOSV has been isolated from brains of Kuhl bats (Pipistrellus kuhlii) in regions of Italy where several TOSV strains are circulating. Recently, positive sera of TOSV and SFSV measured using a seroneutralization test were detected in European bat colonies (free-tailed bats (Tadarida teniotis) and serovar bats (Eptesicus serotinus)) located in the Aragon and the Balearic Islands regions [10]. The possible infection of PONTV II was recorded by way of seroneuralization in Italy [30]. Moreover, Bregalaka virus (BREV) neutralizing antibodies were detected in 3.3% of human sera in North Macedonia, demonstrating that BREV frequently infects humans [10]. Interestingly, PONTV II and BREV were strictly related, having an identity of 96% in the amino acid sequence of the M segment. The clues of possible infectivity of these strictly related viruses for humans deserve more investigation.
In phleboviruses, the same L and S segments could be associated with different M segments. The presence of viruses differing in M segments was already recorded in two group of phleboviruses beyond the PONTV, Massilia, Arrábida and Granada viruses [32,33,34,35] in Europe and in the Candiru phlebovirus [36] in South America. This condition is recorded also in viruses of the Orthobunyavirus genus, in which reassortant viruses are common, but with the relevant difference that the donor virus of the M segment is often known [37]. The great genetic variability of phleboviruses is linked to the absence of the proofreading activity of the RdRp and counterbalanced by a strong and pervasive purifying selection able to preserve the amino acid sequence, which results more stable than the nucleotide sequence from which is translated, above all in the L segment.
Results obtained by the analysis of the PONTV indicated that the three segments were differentially saturated, with the M segment more variable and the L segment more preserved. More work is needed to understand if different saturations of the segments could be caused by a different level of purifying selection and consequent extinction of deleterious genotype, or if other genetic mechanisms are potentially creating this different evolutionary dynamic.
In the M segment, transversions occur more frequently than transitions, indicating a saturation of substitutions. The characterization of the mechanism at the base of this condition was hampered by the variability in M segment. While it is difficult to reconstruct recombination events that could have happened in the saturated condition of the M segment, we failed to detect a recombination signal between sequences of phleboviruses in this study. Even if it is not possible to completely exclude recombination events, this result argues in favor of the rarity of these events, as already highlighted in negative senses RNA viruses [38,39]. This suggests that a mechanism other than recombination can act to improve variability of the capsid genes. Reassortment is likely one of the more relevant mechanisms in producing the diversity between M segments of phleboviruses.
The M segment, which encodes the envelope proteins, was notably more saturated. We hypothesize that this increased variability can favor the escape of the virus from the hosts’ immunity and allow its occupation of different ecological niches. This appears to be confirmed by the different serological reaction in humans and animals of the three PONTV bearing different M segments [32,40].

5. Conclusions

In this study, PONTV III was isolated from sandflies collected in the field. The genomic structure, nucleotide and amino acid sequence homology, and molecular genetics evolution analyses of the virus confirmed that it is a virus of the Adana phlebovirus species. This is a new isolation of PONTV III from sandflies. Further studies should examine the host and epidemiological range of the virus in the field, especially whether the virus can infect and be pathogenic for humans, livestock, and other animals.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/v15020422/s1. Figure S1: Nucleotide and amino acid uncorrected p-distance between L (lower right) and M (upper left) segments of viruses in the “Salehabad complex”. According to the ICTV rule, virus with an AA diversity of L segment < 0.05 must be considered the same species (in orange in the table).

Author Contributions

M.D. (Mohamed Daoudi), M.C., A.B., S.B. and M.D. (Michele Dottori) conceived and designed the study; M.D. (Mohamed Daoudi), M.C., P.B., D.L., G.R., K.M., A.P. and S.B. performed the experiments and analyzed the data; M.D. (Mohamed Daoudi), S.B. and M.C. wrote the paper; M.D. (Mohamed Daoudi), M.C., S.B., P.B., M.D. (Michele Dottori), A.B., A.P., K.M. and D.L. were involved in the interpretation of the results and critical reading of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Emilia-Romagna Regional Plan for the Surveillance and Control of Arboviruses (2021) and by funds “Ricerca Corrente 2022” (IZSLER 08/22 RC) from the Italian Ministry of Health.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors warmly thank Charles H. Calisher for revising the manuscript. We would also like to acknowledge Jenna M. Lacey for the English review.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Plyusnin, A.; Beaty, B.J.; Elliott, R.M.; Goldback, R.; Kormelink, R.; Lundkvist, A. Family Bunyaviridae. In Virus taxonomy: Ninth Report of the International Committee on Taxonomy of Viruses; King, A.M.Q., Adams, M.J., Carstens, E.B., Lefkowitz, E.J., Eds.; Elsevier: London, UK, 2011; pp. 725–741. [Google Scholar]
  2. Horne, K.M.; VanLandingham, D.L. Bunyavirus-Vector Interactions. Viruses 2014, 6, 4373–4397. [Google Scholar] [CrossRef] [PubMed]
  3. Kuhn, J.H.; Adkins, S.; Alioto, D.; Alkhovsky, S.V.; Amarasinghe, G.K.; Anthony, S.J.; Avšič-Županc, T.; Ayllón, M.A.; Bahl, J.; Balkema-Buschmann, A.; et al. 2020 taxonomic update for phylum Negarnaviricota (Riboviria: Orthornavirae), including the large orders Bunyavirales and Mononegavirales. Arch. Virol. 2020, 165, 3023–3072. [Google Scholar] [CrossRef] [PubMed]
  4. Lambert, A.; Hughes, H. Clinically Important Phleboviruses and Their Detection in Human Samples. Viruses 2021, 13, 1500. [Google Scholar] [CrossRef] [PubMed]
  5. Ruigrok, R.W.; Crépin, T.; Kolakofsky, D. Nucleoproteins and nucleocapsids of negative-strand RNA viruses. Curr. Opin. Microbiol. 2011, 14, 504–510. [Google Scholar] [CrossRef] [PubMed]
  6. Raymond, D.D.; Piper, M.E.; Gerrard, S.R.; Smith, J.L. Structure of the Rift Valley fever virus nucleocapsid protein reveals another architecture for RNA encapsidation. Proc. Natl. Acad. Sci. USA 2010, 107, 11769–11774. [Google Scholar] [CrossRef]
  7. Freiberg, A.N.; Sherman, M.B.; Morais, M.C.; Holbrook, M.R.; Watowich, S.J. Three-Dimensional Organization of Rift Valley Fever Virus Revealed by Cryoelectron Tomography. J. Virol. 2008, 82, 10341–10348. [Google Scholar] [CrossRef]
  8. Bouloy, M. Molecular biology of phleboviruses. In Bunyaviridae, Molecular and Cellular Biology; Plyusnin, A., Elliot, R.M., Eds.; Caister Academic Press: Norfolk, UK, 2011; pp. 95–128. [Google Scholar]
  9. Moriconi, M.; Rugna, G.; Calzolari, M.; Bellini, R.; Albieri, A.; Angelini, P.; Cagarelli, R.; Landini, M.P.; Charrel, R.N.; Varani, S. Phlebotomine sand fly–borne pathogens in the Mediterranean Basin: Human leishmaniasis and phlebovirus infections. PLOS Neglected Trop. Dis. 2017, 11, e0005660. [Google Scholar] [CrossRef]
  10. Ayhan, N.; López-Roig, M.; Monastiri, A.; Charrel, R.; Serra-Cobo, J. Seroprevalence of Toscana Virus and Sandfly Fever Sicilian Virus in European Bat Colonies Measured Using a Neutralization Test. Viruses 2021, 13, 88. [Google Scholar] [CrossRef]
  11. Depaquit, J.; Grandadam, M.; Fouque, F.; Andry, P.; Peyrefitte, C. Arthropod-borne viruses transmitted by Phlebotomine sandflies in Europe: A review. Eurosurveillance 2010, 15, 19507. [Google Scholar] [CrossRef]
  12. Elliott, R.M.; Brennan, B. Emerging phleboviruses. Curr. Opin. Virol. 2014, 5, 50–57. [Google Scholar] [CrossRef] [Green Version]
  13. Killick-Kendrick, R. The biology and control of Phlebotomine sandflies. Clin. Dermatol. 1999, 17, 279–289. [Google Scholar] [CrossRef]
  14. Daoudi, M.; Calzolari, M.; Boussaa, S.; Bonilauri, P.; Torri, D.; Romeo, G.; Lelli, D.; Lavazza, A.; Hafidi, M.; Dottori, M.; et al. Identification of Toscana virus in natural population of sand flies (Diptera: Psychodidae) from Moroccan leishmaniasis foci. J. Infect. Public Health 2022, 15, 406–411. [Google Scholar] [CrossRef]
  15. Papa, A.; Velo, E.; Bino, S. A novel phlebovirus in Albanian sandflies. Clin. Microbiol. Infect. 2011, 17, 585–587. [Google Scholar] [CrossRef]
  16. Anagnostou, V.; Pardalos, G.; Athanasiou-Metaxa, M.; Papa, A. Novel Phlebovirus in Febrile Child, Greece. Emerg. Infect. Dis. 2011, 17, 940–941. [Google Scholar] [CrossRef]
  17. Alkan, C.; Allal-Ikhlef, A.; Alwassouf, S.; Baklouti, A.; Piorkowski, G.; de Lamballerie, X.; Izri, A.; Charrel, R. Virus isolation, genetic characterization and seroprevalence of Toscana virus in Algeria. Clin. Microbiol. Infect. 2015, 21, 1040.e1–1040.e9. [Google Scholar] [CrossRef]
  18. Calzolari, M.; Chiapponi, C.; Bellini, R.; Bonilauri, P.; Lelli, D.; Moreno, A.; Barbieri, I.; Pongolini, S.; Lavazza, A.; Dottori, M. Isolation of three novel reassortant phleboviruses, Ponticelli I, II, III, and of Toscana virus from field-collected sand flies in Italy. Parasites Vectors 2018, 11, 84. [Google Scholar] [CrossRef]
  19. Ayhan, N.; Alten, B.; Ivovic, V.; Cvetkovikj, A.; Stefanovska, J.; Martinkovic, F.; Piorkowski, G.; Moureau, G.; Gould, E.A.; Pettersson, J.H.-O.; et al. Field surveys in Croatia and North Macedonia reveal two novel phleboviruses circulating in sandflies. J. Gen. Virol. 2021, 102, 001674. [Google Scholar] [CrossRef]
  20. Lambert, A.J.; Lanciotti, R.S. Consensus Amplification and Novel Multiplex Sequencing Method for S Segment Species Identification of 47 Viruses of the Orthobunyavirus, Phlebovirus, and Nairovirus Genera of the Family Bunyaviridae. J. Clin. Microbiol. 2009, 47, 2398–2404. [Google Scholar] [CrossRef]
  21. Kampmann, M.-L.; Fordyce, S.L.; Ávila-Arcos, M.C.; Rasmussen, M.; Willerslev, E.; Nielsen, L.P.; Gilbert, M.T.P. A simple method for the parallel deep sequencing of full influenza A genomes. J. Virol. Methods 2011, 178, 243–248. [Google Scholar] [CrossRef]
  22. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 2021, 38, 3022–3027. [Google Scholar] [CrossRef]
  23. 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, 268–274. [Google Scholar] [CrossRef] [PubMed]
  24. Martin, D.P.; Varsani, A.; Roumagnac, P.; Botha, G.; Maslamoney, S.; Schwab, T.; Kelz, Z.; Kumar, V.; Ben Murrell, B. RDP5: A computer program for analyzing recombination in, and removing signals of recombination from, nucleotide sequence datasets. Virus Evol. 2020, 7, veaa087. [Google Scholar] [CrossRef]
  25. Posada, D.; Crandall, K.A. Evaluation of methods for detecting recombination from DNA sequences: Computer simulations. Proc. Natl. Acad. Sci. USA 2001, 98, 13757–13762. [Google Scholar] [CrossRef]
  26. Xia, X.; Xie, Z. DAMBE: Software package for data analysis in molecular biology and evolution. J. Hered. 2001, 92, 371–373. [Google Scholar] [CrossRef] [PubMed]
  27. Pond, S.L.K.; Poon, A.F.Y.; Velazquez, R.; Weaver, S.; Hepler, N.L.; Ben Murrell, B.; Shank, S.D.; Magalis, B.R.; Bouvier, D.; Nekrutenko, A.; et al. HyPhy 2.5—A Customizable Platform for Evolutionary Hypothesis Testing Using Phylogenies. Mol. Biol. Evol. 2019, 37, 295–299. [Google Scholar] [CrossRef] [PubMed]
  28. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [PubMed]
  29. Verani, P.; Ciufolini, M.G.; Nicoletti, L.; Balducci, M.; Sabatinelli, G.; Coluzzi, M.; Paci, P.; Amaducci, L. Ecological and epidemiological studies of Toscana virus, an arbovirus isolated from Phlebotomus. Ann. Ist. Super Sanita 1982, 18, 397–399. [Google Scholar]
  30. Calisher, C.; Weinberg, A.; Muth, D.; Lazuick, J. Toscana virus infection in United States citizen returning from Italy. Lancet 1987, 329, 165–166. [Google Scholar] [CrossRef]
  31. Alkan, C.; Bichaud, L.; de Lamballerie, X.; Alten, B.; Gould, E.A.; Charrel, R.N. Sand fly-borne phleboviruses of Eurasia and Africa: Epidemiology, genetic diversity, geographic range, control measures. Antivir. Res. 2013, 100, 54–74. [Google Scholar] [CrossRef]
  32. Percivalle, E.; Cassaniti, I.; Calzolari, M.; Lelli, D.; Baldanti, F. Thirteen Years of Phleboviruses Circulation in Lombardy, a Northern Italy Region. Viruses 2021, 13, 209. [Google Scholar] [CrossRef]
  33. Charrel, R.N.; Moureau, G.; Temmam, S.; Izri, A.; Marty, P.; Parola, P.; da Rosa, A.T.; Tesh, R.B.; de Lamballerie, X. Massilia Virus, A Novel Phlebovirus (Bunyaviridae) Isolated from Sandflies in the Mediterranean. Vector-Borne Zoonotic Dis. 2009, 9, 519–530. [Google Scholar] [CrossRef] [Green Version]
  34. Collao, X.; De Ory, F.; Navarro, J.M.; Lipkin, I.W.; Sanbonmatsu-Gamez, S.L.; Palacios, G.; Tenorio, A.; Molina, R.; Sánchez-Seco, M.P.; Pérez-Ruiz, M.; et al. Granada Virus: A Natural Phlebovirus Reassortant of the Sandfly Fever Naples Serocomplex with Low Seroprevalence in Humans. Am. J. Trop. Med. Hyg. 2010, 83, 760–765. [Google Scholar] [CrossRef]
  35. Amaro, F.; Hanke, D.; Zé-Zé, L.; Alves, M.J.; Becker, S.C.; Höper, D. Genetic characterization of Arrabida virus, a novel phlebovirus isolated in South Portugal. Virus Res. 2016, 214, 19–25. [Google Scholar] [CrossRef]
  36. Amaro, F.; Zé-Zé, L.; Lourenço, J.; Giovanetti, M.; Becker, S.; Alves, M. Phylogenetic Analysis of Massilia phlebovirus in Portugal. Viruses 2021, 13, 1412. [Google Scholar] [CrossRef]
  37. Briese, T.; Calisher, C.H.; Higgs, S. Viruses of the family Bunyaviridae: Are all available isolates reassortants? Virology 2013, 446, 207–216. [Google Scholar] [CrossRef]
  38. Palacios, G.; Tesh, R.; da Rosa, A.T.; Savji, N.; Sze, W.; Jain, K.; Serge, R.; Guzman, H.; Guevara, C.; Nunes, M.R.T.; et al. Characterization of the Candiru Antigenic Complex (Bunyaviridae: Phlebovirus), a Highly Diverse and Reassorting Group of Viruses Affecting Humans in Tropical America. J. Virol. 2011, 85, 3811–3820. [Google Scholar] [CrossRef]
  39. Han, G.-Z.; Worobey, M. Homologous Recombination in Negative Sense RNA Viruses. Viruses 2011, 3, 1358–1373. [Google Scholar] [CrossRef]
  40. Lelli, D.; Scanferla, V.; Moreno, A.; Sozzi, E.; Ravaioli, V.; Renzi, M.; Tosi, G.; Dottori, M.; Lavazza, A.; Calzolari, M. Serological Evidence of Phleboviruses in Domestic Animals on the Pre-Apennine Hills (Northern Italy). Viruses 2021, 13, 1577. [Google Scholar] [CrossRef]
Figure 1. Mid-point rooted ML phylogenetic tree obtained by sequences of L segments (A), and M segments (B), of the viruses isolated from sandflies in the Mediterranean and Middle East. Virus isolated in this study is shown as red. Bootstrap values are shown near the nodes.
Figure 1. Mid-point rooted ML phylogenetic tree obtained by sequences of L segments (A), and M segments (B), of the viruses isolated from sandflies in the Mediterranean and Middle East. Virus isolated in this study is shown as red. Bootstrap values are shown near the nodes.
Viruses 15 00422 g001
Figure 2. Nucleotide and amino acid uncorrected p-distance between viruses of the Salehabad complex.
Figure 2. Nucleotide and amino acid uncorrected p-distance between viruses of the Salehabad complex.
Viruses 15 00422 g002
Figure 3. Saturation plots showing transitions and transversions versus GTR distance of the three genomic sequences of the Ponticelli, I Ponticelli II, and Ponticelli III viruses.
Figure 3. Saturation plots showing transitions and transversions versus GTR distance of the three genomic sequences of the Ponticelli, I Ponticelli II, and Ponticelli III viruses.
Viruses 15 00422 g003
Table 1. Sandflies species composition in Casalgrande locality.
Table 1. Sandflies species composition in Casalgrande locality.
18 June 201930 July 201913 August 201927 August 201910 September 201924 September 20194 September 2020Total
Phlebotomus perfiliewi24015510 72
Phlebotomus perniciosus1 1
Phlebotomus spp.015453 670135 *409
Total3194685770135482
* Isolation on cell culture, 6 pools tested, 1 positive.
Table 2. Nucleotidic and amino acidic p distance between L (lower right) and M segments (upper left). According to the ICTV rule, virus with an aminoacid diversity of L segment < 0.05 can be considered the same species.
Table 2. Nucleotidic and amino acidic p distance between L (lower right) and M segments (upper left). According to the ICTV rule, virus with an aminoacid diversity of L segment < 0.05 can be considered the same species.
VirusAbb.ADAVPONTVIPONTVIIPONTVIIIBREVSALVZABAVARBVMVV
Adana virus 195ADAV 0.34–0.330.34–0.310.35–0.320.35–0.310.36–0.340.34–0.320.31–0.270.35–0.32
Ponticelli I virus 220116–3/2013PONTV I0.19–0.04 0.3–0.230.29–0.220.3–0.240.32–0.260.35–0.320.37–0.360.33–0.27
Ponticelli II virus 181135–8/2013PONTV II0.19–0.040.02–0 0.3–0.240.14–0.040.31–0.230.34–0.310.36–0.340.32–0.26
Ponticelli III virus 270315/2020PONTV III0.19–0.040.02–00.01–0 0.3–0.240.31–0.260.34–0.310.37–0.360.32–0.28
Bregalaka virus M31BREV0.19–0.040.12–0.020.12–0.010.12–0.02 0.31–0.230.34–0.310.36–0.340.32–0.27
Salehabad virus I–81SALV0.25–0.140.26–0.140.25–0.140.26–0.140.25–0.15 0.36–0.320.38–0.380.33–0.29
Zaba virus C48ZABAV0.25–0.150.26–0.150.26–0.150.26–0.150.26–0.150.2–0.05 0.36–0.350.33–0.31
Arbia virus ISS PHL18ARBV0.26–0.150.26–0.150.26–0.150.26–0.150.26–0.150.23–0.10.23–0.11 0.36–0.35
Medjerda Valley virus T131MVV0.26–0.140.25–0.150.25–0.140.25–0.140.25–0.140.23–0.10.23–0.10.19–0.04
Alcube virus S20ACBV0.27–0.150.26–0.160.26–0.160.26–0.160.26–0.160.24–0.110.24–0.110.24–0.110.23–0.1
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.

Share and Cite

MDPI and ACS Style

Daoudi, M.; Romeo, G.; Marzani, K.; Petrella, A.; Bonilauri, P.; Lelli, D.; Boumezzough, A.; Boussaa, S.; Dottori, M.; Calzolari, M. New Isolation of Ponticelli III Virus (Bunyavirales: Phenuiviridae) in Emilia-Romagna Region, Italy. Viruses 2023, 15, 422. https://doi.org/10.3390/v15020422

AMA Style

Daoudi M, Romeo G, Marzani K, Petrella A, Bonilauri P, Lelli D, Boumezzough A, Boussaa S, Dottori M, Calzolari M. New Isolation of Ponticelli III Virus (Bunyavirales: Phenuiviridae) in Emilia-Romagna Region, Italy. Viruses. 2023; 15(2):422. https://doi.org/10.3390/v15020422

Chicago/Turabian Style

Daoudi, Mohamed, Giuseppe Romeo, Katia Marzani, Angelica Petrella, Paolo Bonilauri, Davide Lelli, Ali Boumezzough, Samia Boussaa, Michele Dottori, and Mattia Calzolari. 2023. "New Isolation of Ponticelli III Virus (Bunyavirales: Phenuiviridae) in Emilia-Romagna Region, Italy" Viruses 15, no. 2: 422. https://doi.org/10.3390/v15020422

APA Style

Daoudi, M., Romeo, G., Marzani, K., Petrella, A., Bonilauri, P., Lelli, D., Boumezzough, A., Boussaa, S., Dottori, M., & Calzolari, M. (2023). New Isolation of Ponticelli III Virus (Bunyavirales: Phenuiviridae) in Emilia-Romagna Region, Italy. Viruses, 15(2), 422. https://doi.org/10.3390/v15020422

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