Recently, several Insect Specific flaviviruses have been sequenced, characterized, and identified in arthropods [17](Guarido et al., 2021). In this study we aimed to investigate the presence of ISF in mosquitoes and to analyze the natural diversity and phylogenetic relationships among them in the north-central region of Argentina.
The implementation of a screening by means of Sánchez Seco PCR for Flaviviruses [26](Sánchez Seco et al., 2005) made possible to identify the activity in the north-central region of Argentina of three ISF species: CFA, CxFV and KRV. Although only 10 pools belonging to CxFV were able to be amplified by Vázquez PCR technique [30] (Vazquez et al. 2012), our results indicate that the coincidence between both short and large fragments is 100%, making this screening technique PCR an appropriate tool for surveillance and identification of circulating flaviviruses. The use of NS5 gene as target for molecular characterization and identification of flaviviruses is wide used bringing useful information regarding phylogenetic relationships [13, 26, 29, 30 ](Sánchez-Seco et al., 2005; Vázquez et al., 2012; Fang et al., 2018; Talavera et al., 2018). These sequences strongly suggest the presence of other putative viruses not yet isolated, as in the case of OcFV, isolated after the detection of field sequences [5, 6](Calzolari et al., 2012; 2016).
Insect-specific flaviviruses is a highly divergent group within the genus Flavivirus that shares a common ancestor with all other members, including the disease-causing ones. A possible explanation of the origin of ISFVs, accounts for the existence of a co-divergence between vectors and viruses and this hypothesis is supported by the pattern of virus-vector association observed in this work (Fig. 2). In agreement with other studies [1, 19](Huanyu et al., 2012; Bittar et al., 2016), there is a clear divergence between Clade 1: ISFs associated with mosquitoes of the genus Aedes / Stegomyia (CFA, KRV, AeFV and SOcFV) and Clade 2: ISFVs associated with Culex (CLBV, QBV, CxFV, SCxFV). Despite the hypothesis that ISFV underwent multiple introductions with frequent vector changes and occasional genetic recombinations, particularly in CFAV, however, there are no studies to support that CxFV was generated by codivergence between these viruses and their vectors [1, 7, 9](Cook et al. 2009; Bittar et al., 2016; Chatterjee et al., 2021). As expected, based on previous evidence, CFAV (Clade 1) were found in mosquitoes of the genus Mansonia and Ochlerotatus which, according to the phylogeny of mosquitoes, are more related to the genus Aedes / Stegomyia than to Culex [10](Cook et al. 2012).
CxFV primarily infects globally distributed mosquito species of the genus Culex, which are vectors for pathogenic flaviviruses like WNV, SLEV, and Japanese encephalitis virus (JEV) [1](Bittar et al., 2016). According to other studies [14, 25](Fernandez et al., 2018; Miranda et al., 2018), CxFV has a greater diversity of reservoirs and its phylogeny could be influenced by these. It is the first time that CxFV has been detected in the mosquito species Haemagogus spegazinii, incorporating a greater diversity of reservoirs to this viral species. In Fig. 2, it is shown that this CxFV strain grouped within subgroup 1. The phylogenetic tree shows two clades: the first clade is associated with Asian and USA strains, where Cx. pipiens is the main reservoir. However, other mosquito species, such as Culex tritaeniorhynchus and Anopheles sinensis [22](Liang et al 2015), have been reported infected with CxFV belonging to the first clade. These isolates were classified in genotype I (Asian/USA). The second phylogenetic group of CxFV have been reported mostly in the Americas and Africa, with Culex. quinquefasciatus as the main reservoir [1, 15, 22](Goenaga et al., 2014; Bittar et al., 2016; Liang et al., 2015). The secuences detected in this study are more phylogenetically related to this group, described as genotype II, (Africa/Caribbean/Latin America; Fig. 2). Our findings agree with this phylogeographic separation demonstrated in other studies [15, 22](Goenaga et al., 2014; Liang et al., 2015). Our sequences clustered with previous detections made in Argentina (Buenos Aires, Córdoba, Chaco and Santa Fe provinces) [15](Goenaga et al., 2014). Moreover, all these Argentinian sequences are closed to that detected in Mexico and Africa. A wide circulation of CxFV was observed in the Central-North region of Argentina (Fig. 1). CxFV seems to have been introduced several times in the New World in association with the species of Culex, which could explain the distribution of these viruses detected in this work [3, 10, 14] (Bolling et al., 2011; Cook et al ., 2012; Fernandez et al., 2018).
Although viruses in the genus Flavivirus share complex antigenic relationships, they can be divided into four phylogenetic/ecological groups. These divisions largely reflect the selective constraints imposed on the viruses by the vertebrate hosts, the invertebrate vectors, and the associated ecologies. Some authors, describes the evolution and possible origins of individual flaviviruses, correlating ecological and epidemiological characteristics with their phylogenies and geographic dispersal [6, 12, 16](Gould et al., 2003; Calzolari et al., 2016; de Oliveira Ribeiro et al., 2020). Many of the phylogenetic lineages that define viral species diverged relatively recently, and the subsequent dispersal and epidemiology of these viruses have therefore been significantly influenced by increasing human population densities and activities such as recreation, urbanization, land reclamation, transportation, and deforestation [18](Guzman et al., 2018).
Certain authors suggest that by not replicating in vertebrate cells, the main route of transmission and maintenance of these viruses would be the vertical and venereal route [3, 8, 23](Cook et al. 2006, Lutomiah et al. 2007; Bollling et al. 2012). However, in this study this virus was detected for the first time in another genus (Haemagogus spegazzini), raising questions about its mode of transmission. There would be other forms of transmission not discovered to date.
The temporal analysis of CxFV, KRV and CFA in the city of Córdoba gives evidence of an erratic maintenance with the presence of presence-absence events and with a permanence of the genetic identity of the circulating CxFV over time in the city of Córdoba. This is compatible with the seasonal fluctuations detected in the Culex mosquito populations for the temperate regions of our country. The explanation for the absence of this virus during certain periods could be due to the scarcity of mosquitoes collected / analyzed.
In conclusion, ISFVs, more specifically CxFV, KRV and CFAV, circulate in several mosquitoes species in the Central Region of Argentina. The presence of these viruses in mosquitoes could play an important role from the public health perspective, because it has been shown that previous CxFV infection can increase or block the infection of the mosquito by other pathogenic flaviviruses. Information on ISFV is very limited so far and, therefore, not much can be concluded about their frequency, distribution, and host range like in this work.