Abstract
Equine Infectious Anemia virus (EIAV) is the causative agent of one of the most impacting infectious diseases affecting equids. EIAV is genetically diverse with several lineages circulating globally. To elucidate EIAV global spread patterns, we studied its spatiotemporal dynamics applying Bayesian phylodynamic analyses, using a worldwide compiled dataset composed of unique sequences of the gag gene. In addition, we performed a scoping review of 1.470 publications on EIAV to characterize the spatiotemporal trends in EAIV research. Phylogeographic reconstruction suggested Hungary as the most likely country of origin for current EIAV circulation (root state posterior probability = 0.21), and one of the most important centers of diversification for the disease. Historical EIAV spread was predominantly characterized by long-distance spread across continents. As a result, we found that the American and Asian circulating EIAV (i.e., Mongolia) are more related to European lineages than to other Asian countries, being Europe the continent with the highest EIAV phylogenetic diversity. Our bibliometric analysis showed a continuous increment in the number of publications per year, where the United States and China appeared as the countries with highest EIAV-related scientific production. This study provides a historical geographic mapping of the EIAV linages’ spread patterns. In addition, we identified important asymmetry between the current research effort and the availability of genetic data which restricted our abilities to quantify for example the time frame of virus dispersal.
Introduction
Equine infectious anemia (EIA) is a viral disease affecting equines that entails great economic losses worldwide mostly due to the obligatory elimination or isolation of infected animals (Montelaro et al., 1993; Souza, 2009; USDA, 2017). EIAV is currently one of the most significant infectious diseases of equids (Cook et al., 2013). It is caused by a lentivirus (EIAV) from the Retroviridae family (Leroux et al., 2004) known to be mostly transmitted by vectors from the family Tabanidae (Rojas, 1981; Barros and Foil, 2007). Due to the environmental needs of these vectors, EIAV is more prevalent in warm and wet landscapes (Issel and Foil, 1984; Silva et al., 1999; Barros and Foil, 2007; Parr et al., 2014), but owing to a globalized animal trade and fast global warming, it has spread extensively through the globe in the last decades (Lima and Cintra, 2015). In the United States alone, the economic impact of EIAV has been estimated to be US $38 million, which may be much bigger in countries with higher horse-related industry, such as Brazil (Diniz et al., 2017), known to have the largest population of equids in Latin America (Oliveira et al., 2017).
The field of evolutionary epidemiology has seen the implementation of integrative analyses of enormous multivariate datasets that encompass not only genetic and phenotypic data, but also geographic and ecological information. Besides, within the field of evolutionary epidemiology, molecular phylogenetics is increasingly becoming an accurate and impacting approach owing to its combination of genetic and environmental data to predict and control infectious diseases (Galvani, 2003; Mideo et al., 2008; Lion and Gandon, 2016). As a result, an increasing number of studies have finely quantified the factors and trajectories of lineage evolution over space and time. For example, Chu et al. (2013) revealed the spatial and phylogenetic radiation routes followed by Dengue virus worldwide. Later, Alkhamis et al. (2018) showed the phylodynamic spread patterns of African Swine Fever in Africa and Eurasia, while Jamal and Belsham (2018), Omondi et al. (2019) and Jara et al. (2019), performed similar studies focused on the molecular evolution and epidemiology of Foot-and-Mouth Disease virus. The addition of spatial analyses to these methods has proven to be a remarkably useful tool to identify hot areas where treatment and control activities should be deployed. This combined approach has been widely applied to explain virulence, spread, invasion and species establishment of pathogens worldwide (Magiorkinis et al., 2009; Auguste et al., 2015; Pigeault et al., 2015; Weaver and Forrester, 2015; Dellicour et al., 2016, 2017; Rogalski et al., 2017; Forni et al., 2018; Fountain-Jones et al., 2018). However, these methodologies are still not commonly applied for food animal diseases and therefore, major evolution and global diversification of several infectious disease remain an important bottleneck for better and universally applicable molecular testing (Dorey-Robinson et al., 2019), especially in diseases affecting horse populations, such as EIAV (Antolin, 2008; Struchiner et al., 2008; Dorey-Robinson et al., 2019).
Phylogenetic relationships of EIAV have been previously assessed by the analysis of outbreak cases reported worldwide (Nagarajan and Simard, 2007; Cappelli et al., 2011; Capomaccio et al., 2012), mainly focused in the areas with high horse-related industry, such as Brazil (Bicout et al., 2006; Nogueira et al., 2017; Barzoni et al., 2018; Cursino et al., 2018). However, the evolutionary and spread history of EIAV is still unknown (Dorey-Robinson et al., 2019).
In this study, we investigated the spatiotemporal dynamics of EIAV by using Bayesian phylodynamic analyses of a comprehensive genetic, geographical and temporal dataset regarding EIAV occurrences. In addition, we evaluated the global publication trends on EIAV, based on the assumption that temporal changes in research productivity may identify research priority needs. Results may be used to advance EAIV surveillance and support regional disease control policies related to intervention.
Materials and Methods
Data Collection and Curation
We built a comprehensive genetic database by gathering information of the EIAV gag gene sequence (∼1,400 bp) along with the collection date and location as metadata (GenBank ID numbers and metadata available in Supplementary Material). This gene was chosen due to its extended use in literature for phylogenetic studies since it is one of the most conserved regions of the lentivirus genome (Olsen, 1998; Cappelli et al., 2011; Cook et al., 2013), allowing the comparison and update of previously published information to better understand the evolutionary pathways of this pathogen (Capomaccio et al., 2012). Samples lacking any of the metadata information were discarded for subsequent analyses. This dataset comprises 297 sequences of 17 countries from 1947 to 2015. To determine phylogenetic relationship accurately, and to avoid sampling bias among EIAV sequences, we removed all duplicated sequences that represents the same outbreak. Thus, the cleaned database comprises 51 sequences. Data is available in the Supplementary Table S1.
Discrete Phylogeographical Analysis
Sequences were aligned using Mega X, available at www.megasoftware.net (Kumar et al., 2008). The recombination detection program (RDP) v5.3 (Martin et al., 2015) was used to search for evidence of recombination within our dataset using five different methods (BootScan, Chimaera, MaxChi, RDP, and SiScan). Thus, four recombinant sequences were observed in EIAV dataset.
To determine whether there was a sufficient temporal molecular evolutionary signal of the EIAV sequences used for each serotype phylogeny, we used TempEst v1.5 (Rambaut et al., 2016). To calculate the P-values associated with the phylogenetic signal analysis, we used the approach described by Murray et al. (2016) based on 1,000 random permutations of the sequence sampling dates (Navascués et al., 2010).
Phylogeographic history of EIAV dispersal was recovered from the obtained spatiotemporal phylogeny. The phylogenetic tree was generated by a discrete phylogeography estimation by Bayesian inference through Markov Chain Monte Carlo (MCMC), implemented in BEAST v2.5.0 (Bouckaert et al., 2014), applying a Hasegawa-Kishino-Yano (HKY+G; Hasegawa et al., 1985) substitution model with gamma-distributed rate heterogeneity among sites (Shapiro et al., 2006), to determine the best fitting model of evolution we used PhyML with Smart Model Selection (SMS) (Lefort et al., 2017). By using Nested Sampling Beast package v1.0.4 (Russel et al., 2019) we compared different molecular clock models to find the one that showed the best fit for our data. The marginal likelihood value supported the use of uncorrelated lognormal relaxed molecular clock (Drummond et al., 2006) using a constant population size model. To infer the epidemic demographic histories of EIAV we estimated the effective number of infections through time by using the Bayesian skyline plot approach (Drummond et al., 2005). All analyses were developed for 200 million generations, sampling every 10,000th generation and removing 10% as chain burn-in. The Markov Chain Monte Carlo analysis was investigated using Tracer software v1.7 (Rambaut et al., 2018) to ensure adequate effective sample sizes (ESS) (above 200) which were obtained for all parameters. Final trees were summarized and visualized via Tree Annotator v. 2.3.0 and FigTree 1.4.3 respectively (included in BEAST v2.5.0) (Rambaut and Drummond, 2016; Rambaut, 2017).
To reconstruct the ancestral-state phylogeographic transmission across countries, we used the discrete-trait extension implemented in BEASTv2.5.0 (Rambaut and Drummond, 2016; Rambaut, 2017). In addition, to explore the most important historical dispersal routes for the spread of EIAV across countries, we used a Bayesian stochastic search variable selection (BSSVS) (Lemey et al., 2009). Using BSSVS approach, we identified and eliminated the nonzero rates of change between each pair of discrete traits (countries) based on its Bayes factor value obtained (lower than 3). To perform this analysis, a symmetric rate matrix was assumed (Drummond and Bouckaert, 2015). To infer the intensity of directional transitions (forward and backward) within a matrix of the discrete trait (countries), we used a Markov jumps approach. To interpret the Bayes factors, a value of <3 is not considered significant, BF = 3.1–20 represents positive support, BF = 20.1–150 represents strong support, while >150.1 represents an overwhelming support (Drummond and Bouckaert, 2015). Finally, to visualize the spatiotemporal diffusion of EIAV we used Spatial Phylogenetic Reconstruction of Evolutionary Dynamics using Data-Driven Documents (D3) SPREAD3 software (Bielejec et al., 2016).
Global EIAV Research Trends
We developed an overview of the research done on EIAV through a scoping review by searching for publications with “equine infectious anemia” in the title using Scopus database. We captured information related about authors, author affiliation, journal, year, title, abstract and keywords. Bibliometric analysis was performed in R using the Bibliometrix package v2.2.0 (Aria and Cuccurullo, 2017), obtaining a total of 1.470 publications from 1946 to 2019, from 57 countries. We assigned the publications to each country by searching for the country name within the corresponding author’s affiliation, title, and abstract. Publications with multiple countries listed within the title and abstract were assigned to all countries involved. All country name variants were standardized. Using this information, we mapped the scientific productivity per country. To visualize the average annual trends, we fitted a locally weighted regression (“LOESS”) line (Cleveland, 1979). In addition, we visualized the multinational collaboration network of the 20 most productive countries on EIAV research. Finally, we extracted the total horse population per country (Gilbert et al., 2018) which was used to correlate with scientific production.
Results
We found a total of 297 EIAV gag gene sequences from the 17 countries. After removing all highly recombinant sequences, as well as the similar sequences that represent the same outbreak event, 47 sequences were used for the phylogeographical analysis (Figure 1 and Supplementary Table S1 for more details). The relationship found between root-to-tip divergence and sampling dates (years) supported the use of molecular clock analysis in this study. The obtained results supported a significant temporal signal (r2 = 0.71, p-value <0.05).
FIGURE 1
Phylogeography of Equine Infectious Anemia Virus
The phylogeographic reconstruction of EIAV indicated Hungary as the most likely origin of EIAV [root state posterior probability (RSPP) = 0.21]. On the other hand, Italy and the United States represented the highest source of dispersal events, here represent regions from where a considerable number of lineages spread to other areas (Figure 1A). Bayesian skyline plot (BSP) showed a steady pattern in the genetic diversity (population size) until around 1970, when there was a sharp increase in its effective population size, evidencing a constant decrease around mid-1990s follow another increase in diversity in early 2000 (Figure 1B).
From early times on its evolutionary history, EIAV transmission has been characterized by its spread to distant regions (i.e., to countries not sharing international borders with the country of origin). This long-distance spread occurred in all directions, as seen in the case of Hungary, where the disease spread through Europe and Asia, including very distant countries, such as Russia, Japan, and Thailand. Likewise, we observed the outbreak occurred in United States in 1947, this event gave origin to the spread of EIAV across most of the American countries. On the other hand, Italy was the most likely responsible for the first introduction of the virus into Romania and Cuba (Figure 2A). Another important aspect of the spatiotemporal dynamics of EIAV was the multi-directionality of its spread, expressed as the movement of the disease in both directions between Hungary-United States and United States-Canada. Bayesian stochastic search variable selection (BSSVS-BF) results showed the most significant viral transmission routes for EIAV, where the most intense were represented from Germany to Japan (BF = 16,080.9), followed by the spread from Belgium to the United States (BF = 804.5), from Mongolia to Cuba (BF = 402.2) and from Hungary to Ireland (BF = 201.1) (Figure 2B).
FIGURE 2
Scientific Production on Equine Infectious Anemia Virus
Bibliometric analyses resulted in 1.470 unique publications from 1946 to 2019, from 57 countries. There were countries with a minimum of 1 and maximum of 1.291 publications, the most productive countries were the United States, China and the United Kingdom (Figure 3A).
FIGURE 3
The global collaboration network resulted from the 20 most productive countries, revealed that the United States had the largest collaboration network, followed by the United Kingdom, Germany, and France. Asian countries and Australia had the lowest levels of multinational research within this group, and we did not observe any African countries on this list, which indicates their low levels of international collaboration compared to the rest of continents (Figure 3B). Bibliometric analysis also showed that the research effort over time appeared to be relatively constant from 1946 until mid-1960. From this point on there is a clear increment in the scientific productivity which also vary substantially over the time. However, during the period 1980–1985 we observed a reduction in the number of published works. Overall, the average annual growth rate of publications was 23.3%, and the most productive year was 2013 with 59 articles (Figure 3C). The average annual global trends in EIAV scientific productivity showed a Gaussian shape in the number of publications, with a peak of productivity on 2000, after that, we observe a decreasing tendency that still lasts until today (Figure 3C). We also observed that around mid-1970s, Japan was the country with more publications on EIAV. Later, from 1978, United States started consolidating its position as the country with the highest number of publications in the majority of years, showing its highest research effort around 2002. However, from 2007, the scientific productivity of United States showed a continuous decrease, and was surpassed in 2013 by China, that obtained the first position, and later in 2018, Brazil appeared to be on top of the global research productivity on EIAV (Figure 3C).
Unsurprisingly, the highest scientific productivity was related to the countries with largest horse population (i.e., China, United States, and Brazil), being United States the country that exhibited the highest values in almost all the categories. United States showed the highest scientific production on EIAV (Figure 3A), accompanied by the highest number of EIAV genetic sequences, as well as and the most complex collaboration network (Figures 3A–C). China, which is the country with the largest horse population and the second most productive country in terms of scientific articles on EIAV but had not yet contributed to the pool of available sequences (as needed for phylogeographic analysis, see methods for detailed information). However, countries like Hungary, Italy and Ireland, which have a small horse population and little scientific productivity, showed a big amount of publicly available sequences which contributed to the phylogeographic analysis (Figure 4).
FIGURE 4
Discussion
This study represents an attempt to investigate the global spatiotemporal patterns of EIAV from a multidisciplinary perspective. We combined spatial and evolutionary approaches with a scoping review framework to provide an overview of the historical knowledge of the disease. We reconstructed the evolutionary history of the available lineages of EIAV along with the most likely origin of the analyzed sequences, the most relevant countries that have been recognized as important centers of virus diversification, and the most important spread routes. We quantified the global patterns of the scientific literature.
Phylodynamic Patterns of Equine Infectious Anemia Virus Spread
Equine Infectious Anemia virus is widely distributed throughout the world (Cook et al., 2001). Our phylogeographic results suggested Hungary as most likely center of origin for the currently available sequences. From there, EIAV diversified in several directions, starting with its neighboring European countries (i.e., Slovenia, Bosnia and Herzegovina, Romania and Italy). These findings support what was previously described by Kuhar and Malovrh (2016). Our results showed that the outside-Europe spread occurred first toward the Americas (∼1950), followed by a dispersion through Asian countries (∼2000). Capomaccio et al. (2012) also highlighted the important spread route of EIAV between Europe and the Americas, which could be directly related to the intensive movement of animals between these continents. Likewise, Tigre et al. (2017), found high similarities among Brazilian EIAV isolates and European and North American strains.
Similarly to the historical spread exhibited by other animal diseases, such as Foot and Mouth Disease (Bachanek-Bankowska et al., 2018) or Avian avulavirus (Hicks et al., 2019), EIAV has shown the ability to spread rapidly through long distances, with evident jumps across continents. Spatiotemporal patterns of EIAV visualized through skyline plot showed that its worldwide spread have been accompanied by a constant increase in the genetic diversity that started around 1970 and ceased 20 years later, evidencing an abrupt decrease that could be explained by the intensification of disease control programs in countries with elevated count cases and the development of better ELISA-based EIAV serologic tests (Issel and Cook, 1993). For example, EU is continuously stressing the efforts on developing risk assessments aimed to prevent and reduce the potential introduction of animal diseases into free areas by controlling the movement of live animals across EU countries (Zepeda et al., 2001). Likewise, in United States, the improvements in the surveillance methods have allowed an increase in the number diagnostic tests from 2,777 in 1972 to 993,712 in 1991 (Issel et al., 2013). However, during the last decade, we observed some irregular trends in the genetic diversity of the virus, as a result of numerous outbreaks occurred worldwide, causing dramatic economic consequences (Timoney, 2000; Cook et al., 2013; Herholz et al., 2013). This epidemic pattern has been facilitated by the continuously growing international horse trade (Metcalf, 2001; Dominguez et al., 2015).
Based on our results, it is possible to observe a high genetic diversity of the virus, especially in Europe, where the countries surrounding two of the major centers of diversification (Hungary and Italy) accumulated the highest viral diversity, which supports the findings previously described by Kuhar and Malovrh (2016), and Gaudaire et al. (2018). Our results showed that the spread of EIAV into the Americas occurred due to the arrival of a strain from the United States, which later diversified and spread through Canada and Brazil. These results support a previous study of Tigre et al. (2017) that evidenced strong similarities between the Brazilian strains with the ones circulating in the United States. Moreover, our results also showed that most of the American strains returned to Europe, increasing the similarities between the strains found in these areas. On the other hand, the colonization of Asia seemed to occur exclusively from the arrival of European lineages, which explains the previously reported high similarities between certain strains from both continents in comparison with other lineages from Asia or the Americas (Herholz et al., 2013; Sharav et al., 2017). Finally, while the genetic data used here comes from hotspots areas of EIAV, there are still important gaps of information exemplified in the absence of genetic information in areas such as Africa, where despite it has evidenced the presence of three endemic EIAV serotypes (Herholz et al., 2013; Alnaeem and Hemida, 2019), there is a remarkable lack of genetic information, which could be highly related to an underreporting or lack of laboratory capacities which should be advocated (Barzoni et al., 2018).
Global Trends in Equine Infectious Anemia Virus Research
The bibliometric analysis synthesized the global research trends relevant to map areas that have produce more scientific publications on EIAV. We found that overall, in opposition to Lotka’s law (Coile, 1977), the average number of scientific publications per year showed a steady growth until 2000, followed by a slow reduction maintained until today. However, China and Brazil break this pattern by showing a remarkable increase in their scientific production since around 2005 (Figure 3C). Our results also highlighted the global importance of United States and China as the countries with more publications on EIAV. These results were unsurprising since both of them are known to be within the top countries in scientific production worldwide (Zhou and Leydesdorff, 2006; Nafade et al., 2018). It is well known that variations in scientific productivity between countries can be related to differences in factors such as language, training, development, and funding among others (Ramos et al., 2004), which we propose is one of the main reasons influencing the ability of less-developed countries to carry out the necessary scientific research studies to detect and report EIAV. We observed some clear examples of this in countries like Mexico, Ethiopia, Russia, and Mongolia, which despite being in the top 10 of the countries with higher horse population, showed very low scientific production about EIAV. Likewise, as in previous bibliometric research (Vanni et al., 2014; Ekundayo and Okoh, 2018) we observed asymmetries in the multinational alliances between developing and developed countries that have been historically endemic for EIAV (i.e., African countries).
Overall, our results showed that historically, the United States has been the most productive country which is strongly related to the fact that is one of the countries with highest horse population worldwide, being only surpassed by China (FAO, 2019). Another important aspect to highlight is the complex collaboration network that the United States shows, which has been historically considered as an important factor for the country scientific success. As demonstrated by previous bibliometric analyses that have shown a continuous increment in the multinational scientific collaboration over the past two decades (Cronin et al., 2004; Cronin, 2005; Sonnenwald, 2007), along with a direct link between scientific collaboration and the impact of the produced research (Figg et al., 2006). These differences in collaboration and impact can be observed in a number of Asian countries (i.e., China, Japan, Turkey, and India), which interestingly were observed as some of the less connected countries to the multinational collaboration network and yet were seen to be part of the most productive in relation to EIAV, although the impact of their research tends to be lower.
During the last 5 years, the patterns on the publication rates of the different countries have not followed a similar pattern. We observed that most of the countries, including the United States, have shown a reduction in the number of publications on EIAV, which started on 2003. On the other hand, we observed that China and Brazil increased their scientific productivity with more publications on this topic worldwide, tendency that remained until 2018. This change may be due to the increase on biosafety awareness and the need to understand how to control a disease that directly affects countries’ economies, since they are among the countries with highest equine-related animal trade worldwide. Also shows that countries such as the United States lower their research priorities in the past few years.
Limitations
A relevant limitation encountered during the development of this study was caused by the lack of genetic data and inconsistencies in data availability, which introduced bias to our phylodynamic results (Xia, 2014). This is evidence by the under-representations of some countries with known and long history of being endemic for EIAV form which either no genome sequences were found, or the amount of available information was limited. Similarly, countries such as China, South Africa or Nigeria have published information about the disease, but no genetic data (including date of sample collection) are available. These inconsistencies are one of the major factors impeding the development of more accurate analysis assessing EIAV and may influence the results obtained in our phylogeographic analyses, which might need to be updated as more genetic dataset become available. This lack of information is hindering the development of optimal control strategies to minimize the impact of EIAV and other infectious diseases worldwide. Furthermore, this study also has some limitations related to the bibliometric section. (i) We only used publications from Scopus database, and (ii) the countries analyzed as the locations where research was conducted were based on the affiliation address of the first author, which does not necessarily reflect that information, especially in articles with multinational collaborations.
Final Remarks
This study illuminates the EIAV global evolution and the main disproportion between occurrence of EIAV and presence of genetic information. This study suggested that Hungary was likely the origin EIAV and also acted as main source for global dispersal. The most significant viral transmission routes were from Germany to Japan. Here we identified countries with great importance for EIAV spread, consequently our results can be used to redirect resource allocation such as those dedicated to expanding molecular diagnostic and disease surveillance. Importantly, we highlight the importance of the genetic characterization and detailed report of the virus in future outbreaks. In addition, the incorporation of a bibliometric analysis provided a holistic view of the historical global research efforts on EIAV, ultimately provided evidence of knowledge gaps which can be used for example by local stakeholder while defining disease control programs and prioritize funding for future research.
Statements
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found in the Supplementary Material.
Author contributions
MJ contributed to the design of the project, the phylogeographic analyses, the visualization of the results, and wrote the manuscript. AF-D-D contributed to the data collection, and wrote the manuscript. GM contributed to the design of the project, supervised the methods and wrote the manuscript.
Acknowledgments
We acknowledge the Department of Population Health and Pathobiology- North Carolina State University provided startup funds for GM and MJ. We also acknowledge RISF and CVM-NCSU for funding this project partially.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2020.00127/full#supplementary-material
References
1
AlkhamisM. A.GallardoC.JuradoC.SolerA.AriasM.Sánchez-VizcaínoJ. M. (2018). Phylodynamics and evolutionary epidemiology of African swine fever p72-CVR genes in Eurasia and Africa.PLoS One13:e0192565. 10.1371/journal.pone.0192565
2
AlnaeemA. A.HemidaM. G. (2019). Surveillance of the equine infectious anemia virus in eastern and central Saudi Arabia during 2014-2016.Vet. World12719–723. 10.14202/vetworld.2019.719-723
3
AntolinM. F. (2008). Unpacking β: within-host dynamics and the evolutionary ecology of pathogen transmission.Annu. Rev. Ecol. Evol. Syst.39415–437. 10.1146/annurev.ecolsys.37.091305.110119
4
AriaM.CuccurulloC. (2017). bibliometrix: an R-tool for comprehensive science mapping analysis.J. Inform.11959–975. 10.1016/j.joi.2017.08.007
5
AugusteA. J.LiriaJ.ForresterN. L.GiambalvoD.MoncadaM.LongK. C.et al (2015). Evolutionary and ecological characterization of mayaro virus strains isolated during an outbreak, Venezuela, 2010.Emerg. Infect. Dis.211742–1750. 10.3201/eid2110.141660
6
Bachanek-BankowskaK.Di NardoA.WadsworthJ.MiouletV.PezzoniG.GrazioliS.et al (2018). Reconstructing the evolutionary history of pandemic foot-and-mouth disease viruses: the impact of recombination within the emerging O/ME-SA/Ind-2001 lineage.Sci. Rep.8:14693. 10.1038/s41598-018-32693-8
7
BarrosA. T. M.FoilL. D. (2007). The influence of distance on movement of tabanids (Diptera: Tabanidae) between horses.Vet. Parasitol.144380–384. 10.1016/j.vetpar.2006.09.041
8
BarzoniC. S.NogueiraD. M. P.MarquesG. D.DiehlG. N.PellegriniD.daC. P.et al (2018). Equine infectious anemia in the western region of Rio Grande do Sul, Brazil.Ciênc. Rural48:e20170809. 10.1590/0103-8478cr20170809
9
BicoutD. J.CarvalhoR.Chalvet-MonfrayK.SabatierP. (2006). Distribution of equine infectious anemia in horses in the North of Minas Gerais State, Brazil.J. Vet. Diagn. Investig.18479–482. 10.1177/104063870601800511
10
BielejecF.BaeleG.VranckenB.SuchardM.RambautA.LemeyP. (2016). SpreaD3: interactive visualization of spatiotemporal history and trait evolutionary processes.Mol. Biol. Evol.332167–2169. 10.1093/molbev/msw082
11
BouckaertR.HeledJ.KühnertD.VaughanT.WuC. H.XieD.et al (2014). BEAST 2: a software platform for Bayesian evolutionary analysis.PLoS Comput. Biol.10:e1003537. 10.1371/journal.pcbi.1003537
12
CapomaccioS.WillandZ. A.CookS. J.IsselC. J.SantosE. M.ReisJ. K. P.et al (2012). Detection, molecular characterization and phylogenetic analysis of full-length equine infectious anemia (EIAV) gag genes isolated from Shackleford Banks wild horses.Vet. Microbiol.157320–332. 10.1016/j.vetmic.2012.01.015
13
CappelliK.FelicettiM.MarenzoniM. L.CoppolaG.Verini-suppliziA.ColettiM.et al (2011). Molecular detection, epidemiology, and genetic characterization of novel european field isolates of equine infectious anemia virus.J. Clin. Microbiol.4927–33. 10.1128/JCM.01311-10
14
ChuP.-Y.KeG.-M.ChenP.-C.LiuL.-T.TsaiY.-C.TsaiJ.-J. (2013). spatiotemporal dynamics and epistatic interaction sites in dengue virus type 1: a comprehensive sequence-based analysis.PLoS One8:e74165. 10.1371/journal.pone.0074165
15
ClevelandW. S. (1979). Robust locally weighted regression and smoothing scatterplots.J. Am. Stat. Assoc.74829–836. 10.1080/01621459.1979.10481038
16
CoileR. C. (1977). Lotka’s frequency distribution of scientific productivity.J. Am. Soc. Inf. Sci.28366–370. 10.1002/asi.4630280610
17
CookR. F.LerouxC.IsselC. J. (2013). Equine infectious anemia and equine infectious anemia virus in 2013: a review.Vet. Microbiol.167181–204. 10.1016/j.vetmic.2013.09.031
18
CookS. J.CookR. F.MontelaroR. C.IsselC. J. (2001). Differential responses of Equus caballus and Equus asinus to infection with two pathogenic strains of equine infectious anemia virus.Vet. Microbiol.7993–109. 10.1016/S0378-1135(00)00348-5
19
CroninB. (2005). The Hand of Science: Academic Writing and its Rewards.Lanham, MD: Scarecrow Press.
20
CroninB.ShawD.BarreK. (2004). Visible, less visible, and invisible work: patterns of collaboration in 20th century chemistry.J. Am. Soc. Inf. Sci. Technol.55160–168. 10.1002/asi.10353
21
CursinoA.VilelaA.Franco-LuizA. P. M.de OliveiraJ. G.NogueiraM. F.JúniorJ. P. A.et al (2018). Equine infectious anemia virus in naturally infected horses from the Brazilian Pantanal.Arch. Virol.1632385–2394. 10.1007/s00705-018-3877-8
22
DellicourS.LemeyP.RoseR.FariaN. R.PybusO. G.VieiraL. F. P.et al (2017). Using viral gene sequences to compare and explain the heterogeneous spatial dynamics of virus epidemics.Mol. Biol. Evol.342563–2571.
23
DellicourS.RoseR.PybusO. G. (2016). Explaining the geographic spread of emerging epidemics: a framework for comparing viral phylogenies and environmental landscape data.BMC Bioinformatics17:82. 10.1186/s12859-016-0924-x
24
DinizR. S.FreitasF. J. C.LimaJ. M.SakamotoS. M.LeiteR. C.IsselC. J.et al (2017). Equine infectious anemia prevalence in feral donkeys from Northeast Brazil.Prev. Vet. Med.14030–37. 10.1016/j.prevetmed.2017.02.015
25
Dorey-RobinsonD. L. W.LockerN.SteinbachF.ChoudhuryB. (2019). Molecular characterization of equine infectious anaemia virus strains detected in England in 2010 and 2012.Transbound. Emerg. Dis.662311–2317. 10.1111/tbed.13286
26
DominguezM.MünstermannS.MurrayG.TimoneyP. (2015). Highhealth, high-performance’horses: risk mitigation strategies for OIE-listed diseases.Rev. Sci. Tech.34837–848. 10.20506/rst.34.3.2399
27
DrummondA.RambautA.ShapiroB.PybusO. G. (2005). Bayesian coalescent inference of past population dynamics from molecular sequences.Mol. Biol. Evol.221185–1192. 10.1093/molbev/msi103
28
DrummondA. J.BouckaertR. R. (2015). Bayesian Evolutionary Analysis with BEAST.Cambridge: Cambridge University Press.
29
DrummondA. J.HoS. Y. W.PhillipsM. J.RambautA. (2006). Relaxed phylogenetics and dating with confidence.PLoS Biol.4:e88. 10.1371/journal.pbio.0040088
30
EkundayoT. C.OkohA. I. (2018). A global bibliometric analysis of Plesiomonas-related research (1990 – 2017).PLoS One13:e0207655. 10.1371/journal.pone.0207655
31
FAO (2019). FAOSTAT Database.Rome: FAO.
32
FiggW. D.DunnL.LiewehrD. J.SteinbergS. M.ThurmanP. W.BarrettJ. C.et al (2006). Scientific collaboration results in higher citation rates of published articles.Pharmacother. J. Hum. Pharmacol. Drug Ther.26759–767. 10.1592/phco.26.6.759
33
ForniD.CaglianiR.ClericiM.SironiM. (2018). Origin and dispersal of Hepatitis e virus article.Emerg. Microbes Infect.7:11. 10.1038/s41426-017-0009-6
34
Fountain-JonesN. M.PearseW. D.EscobarL. E.Alba-CasalsA.CarverS.DaviesT. J.et al (2018). Towards an eco-phylogenetic framework for infectious disease ecology.Biol. Rev.93950–970. 10.1111/brv.12380
35
GalvaniA. P. (2003). Epidemiology meets evolutionary ecology.Trends Ecol. Evol.18132–139. 10.1016/S0169-5347(02)00050-2
36
GaudaireD.LecouturierF.PonçonN.MorillandE.LaugierC.ZientaraS.et al (2018). Molecular characterization of equine infectious anaemia virus from a major outbreak in southeastern France.Transbound. Emerg. Dis.65e7–e13. 10.1111/tbed.12657
37
GilbertM.NicolasG.CinardiG.Van BoeckelT. P.VanwambekeS. O.WintG. R. W.et al (2018). Global distribution data for cattle, buffaloes, horses, sheep, goats, pigs, chickens and ducks in 2010.Sci. Data5:180227. 10.1038/sdata.2018.227
38
HasegawaM.KishinoH.YanoT. (1985). Dating of the human-ape splitting by a molecular clock of mitochondrial DNA.J. Mol. Evol.22160–174. 10.1007/BF02101694
39
HerholzC.SchwermerH.FüsselA.PerlerL.BinggeliM.Bürgi TschanD.et al (2013). International horse movements and spread of equine diseases: equine infectious Anaemia and Glanders - two examples.Pferdeheilkd. Equine Med.29445–450. 10.21836/pem20130402
40
HicksJ. T.DimitrovK. M.AfonsoC. L.RameyA. M.BahlJ. (2019). Global phylodynamic analysis of avian paramyxovirus-1 provides evidence of inter-host transmission and intercontinental spatial diffusion.BMC Evol. Biol.19:108. 10.1186/s12862-019-1431-2
41
IsselC. J.CookR. F. (1993). A review of techniques for the serologic diagnosis of equine infectious anemia.J. Vet. Diagn. Invest.5137–141. 10.1177/104063879300500136
42
IsselC. J.FoilL. D. (1984). Studies on equine infectious anemia virus transmission by insects.J. Am. Vet. Med. Assoc.184293–297.
43
IsselC. J.SciclunaM. T.CookS. J.CookR. F.CaprioliA.RicciI.et al (2013). Challenges and proposed solutions for more accurate serological diagnosis of equine infectious anaemia.Vet. Rec.172:210. 10.1136/vr-2012-100735
44
JamalS. M.BelshamG. J. (2018). Molecular epidemiology, evolution and phylogeny of foot-and-mouth disease virus.Infect. Genet. Evol.5984–98. 10.1016/j.meegid.2018.01.020
45
JaraM.Frias-De-DiegoA.DellicourS.BaeleG.MachadoG. (2019). Global phylogeographical patterns in the spread of foot-and-mouth disease virus Running title.bioRxiv [Preprint]. 10.1101/590612
46
KuharU.MalovrhT. (2016). High genetic diversity of equine infectious anaemia virus strains from Slovenia revealed upon phylogenetic analysis of the p15 gag gene region.Equine Vet. J.48205–210. 10.1111/evj.12404
47
KumarS.NeiM.DudleyJ.TamuraK. (2008). MEGA: a biologist-centric software for evolutionary analysis of DNA and protein sequences.Brief. Bioinform.9299–306. 10.1093/bib/bbn017
48
LefortV.LonguevilleJ.-E.GascuelO. (2017). SMS: smart model selection in PhyML.Mol. Biol. Evol.342422–2424. 10.1093/molbev/msx149
49
LemeyP.RambautA.DrummondA. J.SuchardM. A. (2009). Bayesian phylogeography finds its roots.PLoS Comput. Biol.5:e1000520. 10.1371/journal.pcbi.1000520
50
LerouxC.CadoréJ. L.MontelaroR. C. (2004). Equine infectious anemia virus (EIAV): What has HIV’s country cousin got tell us?Vet. Res.35485–512. 10.1051/vetres:2004020
51
LimaR. D. S.CintraA. G. (2015). Revisão do Estudo do Complexo do Agronegócio do Cavalo.Brasília: MAPA.
52
LionS.GandonS. (2016). Spatial evolutionary epidemiology of spreading epidemics.Proc. R. Soc. B Biol. Sci.283:20161170. 10.1098/rspb.2016.1170
53
MagiorkinisG.MagiorkinisE.ParaskevisD.HoS. Y. W.ShapiroB.PybusO. G.et al (2009). The global spread of hepatitis C virus 1a and 1b: a phylodynamic and phylogeographic analysis.PLoS Med.6:e1000198. 10.1371/journal.pmed.1000198
54
MartinD. P.MurrellB.GoldenM.KhoosalA.MuhireB. (2015). RDP4: detection and analysis of recombination patterns in virus genomes.Virus Evol.1:vev003. 10.1093/ve/vev003
55
MetcalfE. S. (2001). The role of international transport of equine semen on disease transmission.Anim. Reprod. Sci.68229–237. 10.1016/S0378-4320(01)00159-2
56
MideoN.AlizonS.DayT. (2008). Linking within- and between-host dynamics in the evolutionary epidemiology of infectious diseases.Trends Ecol. Evol.23511–517. 10.1016/j.tree.2008.05.009
57
MontelaroR. C.BallJ. M.RushlowK. E. (1993). “Equine retroviruses,” in The Retroviridae, ed.LevyJ. A. (New York, NY: Plenum Press), 257–360. 10.1007/978-1-4899-1627-3_5
58
MurrayG. G. R.WangF.HarrisonE. M.PatersonG. K.MatherA. E.HarrisS. R.et al (2016). The effect of genetic structure on molecular dating and tests for temporal signal.Methods Ecol. Evol.780–89. 10.1111/2041-210X.12466
59
NafadeV.NashM.HuddartS.PandeT.GebreselassieN.LienhardtC.et al (2018). A bibliometric analysis of tuberculosis research, 2007–2016.PLoS One13:e0199706. 10.1371/journal.pone.0199706
60
NagarajanM. M.SimardC. (2007). Gag genetic heterogeneity of equine infectious anemia virus (EIAV) in naturally infected horses in Canada.Virus Res.129228–235. 10.1016/j.virusres.2007.07.013
61
NavascuésM.DepaulisF.EmersonB. (2010). Combining contemporary and ancient DNA in population genetic and phylogeographical studies.Mol. Ecol. Resour.10760–772. 10.1111/j.1755-0998.2010.02895.x
62
NogueiraM. F.OliveiraJ. M.SantosC. J. S.PetzoldH. V.AguiarD. M.JulianoR. S.et al (2017). Equine infectious anaemia in equids of Southern Pantanal, Brazil: seroprevalence and evaluation of the adoption of a control programme.Pesq. Vet. Bras.37227–233. 10.1590/s0100-736x2017000300005
63
OliveiraF. G.CookR. F.NavesJ. H. F.OliveiraC. H. S.DinizR. S.FreitasF. J. C.et al (2017). Equine infectious anemia prevalence in feral donkeys from Northeast Brazil.Prev. Vet. Med.14030–37. 10.1016/J.PREVETMED.2017.02.015
64
OlsenJ. C. (1998). Gene transfer vectors derived from equine infectious anemia virus.Gene Ther.51481–1487. 10.1038/sj.gt.3300768
65
OmondiG.AlkhamisM. A.ObandaV.GakuyaF.SangulaA.PauszekS.et al (2019). Phylogeographical and cross-species transmission dynamics of SAT1 and SAT2 foot-and-mouth disease virus in Eastern Africa.Mol. Ecol.282903–2916. 10.1111/mec.15125
66
ParrC. S.WilsonN.LearyP.SchulzK.LansK.WalleyL.et al (2014). The encyclopedia of life v2: providing global access to knowledge about life on earth.Biodivers. Data J.2:e1079. 10.3897/BDJ.2.e1079
67
PigeaultR.VézilierJ.CornetS.ZéléF.NicotA.PerretP.et al (2015). Avian malaria: a new lease of life for an old experimental model to study the evolutionary ecology of Plasmodium.Philos. Trans. R. Soc. B Biol. Sci.370:20140300. 10.1098/rstb.2014.0300
68
RambautA. (2017). FigTree-version 1.4.3, A Graphical Viewer of Phylogenetic Trees. Available online at: https://bioweb.pasteur.fr/packages/pack@FigTree@1.4.3(accessed January 15, 2019).
69
RambautA.DrummondA. (2016). TreeAnnotator v.2.4.3. Available online at: http://beast.community/treeannotator(accessed January 15, 2019).
70
RambautA.LamT.Max CarvalhoL.PybusO. (2016). Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen).Virus Evol.2:vew007. 10.1093/ve/vew007
71
RambautA. R.DrummondA. L. J.DongX.BaeleG.SuchardM. (2018). Posterior summarization in bayesian phylogenetics using tracer 1. 7.Syst. Biol.67901–904. 10.1093/sysbio/syy032
72
RamosJ. M.GutiérrezF.MasíaM.Martín-HidalgoA. (2004). Publication of European Union research on infectious diseases (1991-2001): a bibliometric evaluation.Eur. J. Clin. Microbiol. Infect. Dis.23180–184. 10.1007/s10096-003-1074-4
73
RogalskiM. A.GowlerC. D.ShawC. L.HufbauerR. A.DuffyM. A. (2017). Human drivers of ecological and evolutionary dynamics in emerging and disappearing infectious disease systems.Philos. Trans. R. Soc. B Biol. Sci.372:20160043. 10.1098/rstb.2016.0043
74
RojasW. R. (1981). Anemia infecciosa equina.Monogr. Med. Vet.3:560.
75
RusselP. M.BrewerB. J.KlaereS.BouckaertR. R. (2019). Model selection and parameter inference in phylogenetics using nested sampling.Syst. Biol.68219–233. 10.1093/sysbio/syy050
76
ShapiroB.RambautA.DrummondA. J. (2006). Choosing appropriate substitution models for the phylogenetic analysis of protein-coding sequences.Mol. Biol. Evol.237–9. 10.1093/molbev/msj021
77
SharavT.KonnaiS.OchirkhuuN.TsE. O.MekataH.SakodaY.et al (2017). Detection and molecular characterization of equine infectious anemia virus in Mongolian horses.J. Vet. Med. Sci.791884–1888. 10.1292/jvms.17-0202
78
SilvaR. A. M. S.RamirezL.SouzaS. S.OrtizA. G.PereiraS. R.DávilaA. M. R. (1999). Hematology of natural bovine trypanosomosis in the Brazilian Pantanal and Bolivian wetlands.Vet. Parasitol.8587–93. 10.1016/S0304-4017(99)00081-3
79
SonnenwaldD. H. (2007). Scientific collaboration.Annu. Rev. Inf. Sci. Technol.41643–681. 10.1002/aris.2007.1440410121
80
SouzaJ. G. (2009). Programas Nacionais de Saúde Animal do Brasil.Brasília: MAPA.
81
StruchinerC. J.LuzP. M.CodeçoC. T.MassadE. (2008). The many faces of epidemiology: evolutionary epidemiology.Ciên. Saúde Colet.131743–1752. 10.1590/S1413-81232008000600009
82
TigreD. M.BrandãoC. F. L.de PaulaF. L.ChinaliaF. A.CamposG. S.SardiS. I. (2017). Characterization of isolates of equine infectious anemia virus in Brazil.Arch. Virol.162873–877. 10.1007/s00705-016-3172-5
83
TimoneyP. J. (2000). Factors influencing the international spread of equine diseases.Vet. Clin. North Am. Equine Pract.16537–551. 10.1016/S0749-0739(17)30094-9
84
USDA (2017). Animal and Plant Health Inspection Service: Equine Infectious Anemia (EIA).Washington, DC: USDA.
85
VanniT.Mesa-FriasM.Sanchez-GarciaR.RoeslerR.SchwartsmannG.GoldaniM. Z.et al (2014). International scientific collaboration in HIV and HPV: a network analysis.PLoS One9:e93376. 10.1371/journal.pone.0093376
86
WeaverS. C.ForresterN. L. (2015). Chikungunya: evolutionary history and recent epidemic spread.Antiviral Res.12032–39. 10.1016/j.antiviral.2015.04.016
87
XiaX. (2014). “Phylogenetic bias in the likelihood method caused by missing data coupled with among-site rate variation: an analytical approach,” in Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics), edsBasuM.PanY.WangJ. (Cham: Springer Verlag), 12–23. 10.1007/978-3-319-08171-7_2
88
ZepedaC.SalmanM.RupannerR. (2001). International trade, animal health and veterinary epidemiology: challenges and opportunities.Prev. Vet. Med.48261–271.
89
ZhouP.LeydesdorffL. (2006). The emergence of China as a leading nation in science.Res. Policy3583–104.
Summary
Keywords
equine anemia, infectious diseases, molecular epidemiology, phylodynamics, viral evolution
Citation
Jara M, Frias-De-Diego A and Machado G (2020) Phylogeography of Equine Infectious Anemia Virus. Front. Ecol. Evol. 8:127. doi: 10.3389/fevo.2020.00127
Received
20 December 2019
Accepted
17 April 2020
Published
13 May 2020
Volume
8 - 2020
Edited by
Ajit Chande, Indian Institute of Science Education and Research, India
Reviewed by
Raj Kumar Singh, National Research Centre on Equines and Veterinary Type Culture Collection, India; Eleonora Cella, University of Central Florida, United States
Updates
Copyright
© 2020 Jara, Frias-De-Diego and Machado.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gustavo Machado, gmachad@ncsu.edu
This article was submitted to Phylogenetics, Phylogenomics, and Systematics, a section of the journal Frontiers in Ecology and Evolution
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.