Hepatitis C virus genotype 3A in a population of injecting drug users in Montenegro: Bayesian and evolutionary analysis.

Few reports are available on HCV molecular epidemiology among IDUs in Eastern Europe, and none in Montenegro. The aim of this study was to investigate the HCV genotype distribution in Montenegro among IDUs and to perform Bayesian and evolutionary analysis of the most prevalent HCV genotype circulating in this population. Sixty-four HCV-positive IDUs in Montenegro were enrolled between 2013 and 2014, and the NS5B gene was sequenced. The Bayesian analysis showed that the most prevalent subtype was HCV-3a. Phylogenetic data showed that HCV-3a reached Montenegro in the late 1990s, causing an epidemic that exponentially grew between the 1995 and 2005. In the dated tree, four different entries, from 1990 (clade D), 1994 (clade A) to 1999 (clade B) and 2001 (clade C), were identified. In the NS5B protein model, the amino acids variations were located mainly in the palm domain, which contains most of the conserved structural elements of the active site. This study provides an analysis of the virus transmission pathway and the evolution of HCV genotype 3a among IDUs in Montenegro. These data could represent the basis for further strategies aimed to improve disease management and surveillance program development in high-risk populations.

Hepatitis B virus genotype and subgenotype prevalence and distribution in Montenegro.

The Mediterranean area and the Balkans in particular show the highest level of genetic heterogeneity of HBV in Europe. Data about the circulation of HBV genotypes in Montenegro are lacking. It was studied the prevalence and distribution of HBV genot/subgenotypes in a total of 150 HBV infected patients living in Montenegro. Phylogenetic analysis of 136 successfully amplified P sequences showed a high degree of genetic heterogeneity of HBV in Montenegro. Subgenotype D2 (36.8%) and D3 (32.3%) were the most prevalent, followed by genotype A (subgenotype A2 in all of the cases-19.8%). Eight patients were infected with recombinant strains. HBV-D1 which is the most spread HBV subgenotype in the south-eastern Mediterranean countries, seems to be relatively rare in Montenegro, suggesting a penetration of HBV more probably from North-East or West than from Eastern Mediterranean countries. The relatively different prevalence of D3 and A2 among subjects infected through sexual route, seems to confirm the association of these subgenotypes with different route of transmissions (mainly parenteral for D3 and mainly sexual for A2) even in Montenegro. The low prevalence of D2 among children and its absence in perinatal transmission, suggests that this subgenotype circulated prevalently in the past. If this is due to changes in the most prevalent way of transmission and in the recent different contacts of Montenegro with other European countries, it remains to be established by other larger studies.

Recurrence of another hepatitis B virus escape mutant comes back in a patient infected with HIV and low CD4+ count.

The case of an HIV-positive patient who cleared formerly acute HBV infection and experienced HBV exacerbation before starting antiretroviral therapy (ART) followed by a severe elevation of liver enzymes one month after the introduction of ART is reported. Laboratory tests revealed the presence of a subgenotype A2 carrying at least seven different mutations within the S protein and high lopinavir plasma concentrations. Irrespectively of the many possible mechanisms of HBV reactivation in patients infected with HIV (escape mutants selection, host immune deficiency progression, antiretroviral-related immune reconstitution, glucocorticosteroids, chemotherapy), the message is that HIV patients should be monitored strictly not only by serological and virological assays to recognize early HBV reactivation, but also for drug over-exposure.

Phylogenetic analysis of HCV-4d in Turkey: the curious case of Kayseri province.

In Turkey, genotype 1, especially type 1b virus, causes approximately 90% of these infections, while types 2, 3, and 4 exist, albeit in low prevalences and are due to relatively recent and limited introductions. Two recent reports from Kayseri, a relatively large city in Central Anatolia, indicated unusually high prevalence for type 4 infections in the province reaching a 35% among patients admitted to hospitals for treatment of chronic hepatitis C. In this study, the origin, the demographic history, and the dynamic of the epidemic of unusual HCV genotype 4d in Turkey by using Bayesian coalescent-based method were investigated. A gene flow migration approach was also used to describe the synchronous geographical dispersal and genetic diversification of this unusual genotype in Kayseri province. The Turkish clade had a tMRCA of 44 years corresponding to the year 1967 and seems to have a different origin being completely segregated from the European one. Gene flow migration analysis indicated that Kayseri province appeared to be the epicenter of HCV-4d epidemic, exporting the infections. The demographic history of HCV-4d showed that the epidemic started in 1970s year then following a slow exponential growth until 1980s. The Turkish monophyletic clade suggests a segregate circulation of the epidemic in this region mostly due to unsafe parenteral medical procedures (with drug addiction playing a relatively negligible role).

Evolutionary dynamics of HBV-D1 genotype epidemic in Turkey.

Hepatitis B virus (HBV), is the leading cause of liver diseases infecting an estimated 240 million persons worldwide. The HBV prevalence rates are variables between different countries, with an high level of endemicity in the south-eastern part of Europe. Seven main HBV-D subgenotypes have been described until now (D1-D7). Turkey, seems to have played an important role in the penetration of HBV-D1 in the Mediterranean area. The importance of Turkey in the European epidemiology of HBV is also suggested by the observation that the highest spread of HBV infection in the Continent are reported in Turkey with Romania, Bulgaria, Greece, Albania and some southern regions of Italy. In this paper the molecular epidemiology and the epidemiological history of HBV-D in Turkey was studied, by characterizing 34 new Turkish isolates and performing a phylogeographic reconstruction. By using a phylodynamic and phylogeographic Bayesian approach, the analysis suggested that HBV-D1 originated in Turkey about in the early 1940s. The large prevalence of D1 in comparison to the other subgenotypes in Turkey confirms the importance of this Country as epidemiological reservoir of HBV-D1 dispersion. The phylogeny suggests that after each initial introduction of the virus in a specific population, separate transmission clusters have been evolving along independent phylogenetic lineages. Better characterization and continuous monitoring of such groups are going to be crucial to understand in detail the epidemiology of HBV-D1 subgenotype in Turkey and to assess the efficacy of prevention, vaccination and therapy in controlling the epidemic.

From Clinical Specimen to Whole Genome Sequencing of A(H3N2) Influenza Viruses: A Fast and Reliable High-Throughput Protocol.

(1) Background: Over the last few years, there has been growing interest in the whole genome sequencing (WGS) of rapidly mutating pathogens, such as influenza viruses (IVs), which has led us to carry out in-depth studies on viral evolution in both research and diagnostic settings. We aimed at describing and determining the validity of a WGS protocol that can obtain the complete genome sequence of A(H3N2) IVs directly from clinical specimens. (2) Methods: RNA was extracted from 80 A(H3N2)-positive respiratory specimens. A one-step RT-PCR assay, based on the use of a single set of specific primers, was used to retro-transcribe and amplify the entire IV type A genome in a single reaction, thus avoiding additional enrichment approaches and host genome removal treatments. Purified DNA was quantified; genomic libraries were prepared and sequenced by using Illumina MiSeq platform. The obtained reads were evaluated for sequence quality and read-pair length. (3) Results: All of the study specimens were successfully amplified, and the purified DNA concentration proved to be suitable for NGS (at least 0.2 ng/µL). An acceptable coverage depth for all eight genes of influenza A(H3N2) virus was obtained for 90% (72/80) of the clinical samples with viral loads >105 genome copies/mL. The mean depth of sequencing ranged from 105 to 200 reads per position, with the majority of the mean depth values being above 103 reads per position. The total turnaround time per set of 20 samples was four working days, including sequence analysis. (4) Conclusions: This fast and reliable high-throughput sequencing protocol should be used for influenza surveillance and outbreak investigation.

Time and Mode of Epidemic HCV-2 Subtypes Spreading in Europe: Phylodynamics in Italy and Albania.

Hepatitis C virus (HCV) genotype 2 causes about 10% of global infections and has the most variable circulation profile in Europe. The history of "endemic" HCV-2 subtypes has been satisfactorily reconstructed, instead there is little information about the recent spread of the "epidemic" subtypes, including HCV-2c. To investigate the origin and dispersion pathways of HCV-2c, 245 newly characterized Italian and Albanian HCV-2 NS5B sequences were aligned with 247 publicly available sequences and included in phylogeographic and phylodynamic analyses using the Bayesian framework. Our findings show that HCV-2c was the most prevalent subtype in Italy and Albania. The phylogeographic analysis suggested an African origin of HCV-2c before it reached Italy about in the 1940s. Phylodynamic analysis revealed an exponential increase in the effective number of infections and Re in Italy between the 1940s and 1960s, and in Albania between the 1990s and the early 2000s. It seems very likely that HCV-2c reached Italy from Africa at the time of the second Italian colonization but did not reach Albania until the period of dramatic migration to Italy in the 1990s. This study contributes to reconstructing the history of the spread of epidemic HCV-2 subtypes to Europe.

Genetic variability of the small t antigen of the novel KI, WU and MC polyomaviruses.

Recently, three novel human polyomaviruses were discovered: KIPyV, WUPyV and MCPyV. To study the genetic variability of these viruses, an evolutionary analysis of the large T antigen, small t antigen, VP1, VP2 and VP3 genes was carried out. Of the genes analyzed, only the small t antigen of KIPyV and MCPyV was found to be under statistically supported selective pressure. In KIPyV, 31 selected sites were found under positive selection (ω(2) = 1.54), with posterior probabilities above 0.75. Of these sites, 17 are located in the J domain and two of these are located in the HPDGG box. Thirteen of these sites were identified as being under positive selection at the 99% confidence level. Moreover, seven codons were under significant purifying selection, particularly in the J domain. In MCPyV, we confirmed the five sites already found to be under positive selection (posterior probabilities 0.99). Selective pressure analysis may provide useful information on viral evolution.

Origin and transmission of Feline coronavirus type I in domestic cats from Northern Italy: a phylogeographic approach.

Feline coronavirus (FCoV) is responsible, along with an inadequate immune response of the host, for Feline infectious peritonitis (FIP), one of the most frequent and deadly infectious feline disease worldwide. This study analyzed the genetic characteristics of the spike (S) gene of 33 FCoVs circulating in Northern Italy between 2011 and 2015 in cats with or without FIP. In order to reconstruct the most probable places of origin and dispersion of FCoV among Italian cats, a phylogeographic approach was performed based on 106 FCoV S gene partial sequences from cats, including the 33 novel Italian sequences and 73 retrieved from public databases. Only FCoV type I was found in the Italian cats. The estimated mean evolutionary rate of FCoV was 2.4 × 10-2 subs/site/year (95% HPD: 1.3-3.7 × 10-2), confirming the high genetic variability in the circulating strains. All the isolates clustered in a unique highly significant clade that likely originated from USA between the 1950s and the 1970s, confirming the first descriptions of the disease in American cats. Our results suggest that from USA the virus likely entered Germany and thereafter spread to other European countries. Phylogeography showed that sequences segregated mainly by geographical origin. In the 2010s Italian sequences clustered in different subclades, confirming that different strains cocirculate in Italy. Further studies on archival samples and other genetic regions of FCoV are suggested in order to confirm the present results and to reconstruct a more in-depth detailed virus dispersion pattern for the definition of possible control measures.

Identification and Genetic Characterization of a Novel Respirovirus in Alpine Chamois (rupicapra rupicapra rupicapra).

The Respirovirus genus, family Paramamixoviridae, includes respiratory viral pathogens. Here we report the identification and genetic characterization of a respirovirus in an Alpine chamois showing interstitial pneumonia associated with catarrhal bronchopneumonia. The full-genome characterization of this respirovirus, named ChamoisRV/IT2014, revealed low similarities to caprine respirovirus (77.1%), bovine respirovirus (74.5%) and human respirovirus (72.0%). The phylogenetic analyses based on the full-length genome sequence of the novel isolate and reference respirovirus strains showed that ChamoisRV/IT2014 clustered with caprine respirovirus but formed a separate branch. The phylogenetic tree topology of complete large protein amino acid sequences, representing the current species demarcation criterion for Respirovirus genus, showed a 0.05 branch length of ChamoisRV/IT2014 sequence between the nearest node and the tip of the branch, suggesting that this virus belongs to a novel species. This new isolate in a new host species raises several questions to be addressed on the epidemiological role of chamois and the risks of cross-transmission between wild ruminants and livestock.

In-depth phylogenetic analysis of the hemagglutinin gene of influenza A(H3N2) viruses circulating during the 2016-2017 season revealed egg-adaptive mutations of vaccine strains.

Objectives: The authors conducted an in-depth phylogenetic analysis of the hemagglutinin (HA) gene of A(H3N2) identified during the 2016-2017 season to compare the circulating strains with both the egg-derived and cell-derived 2016-2017 candidate vaccine virus (CVV).Methods: 162 HA1 sequences of H3N2 viruses identified during the 2016-2017 season were phylogenetically analyzed and compared with egg-/cell-derived CVV. The predicted vaccine efficacy (pVE) of egg-/cell-derived CVV against field strains was quantified by Pepitope model.Results: All H3N2 belonged to 3C.2a genetic clade, most (80.2%) to 3C.2a1 sub-clade. Several H3N2 variants circulated in the 2016-2017 season. HA sequences of field H3N2 strains displayed greater identity with cell-derived CVV than with egg-derived CVV. The amino acid substitutions in positions 160 and 194 identified in egg-derived CVV affected the pVE, which was lower for egg-derived CVV than for cell-derived CVV.Conclusions: The results suggested that reduced vaccine effectiveness observed in 2016-2017 season was probably due to changes in the HA of vaccine strains acquired upon adaptation in eggs, which are the basis of most manufacturing systems currently used globally. Egg-free vaccine manufacturing systems would be advisable to improve the effectiveness of influenza vaccines.

Time-scaled phylogeography of complete Zika virus genomes using discrete and continuous space diffusion models.

Zika virus (ZIKV), a vector-borne infectious agent that has recently been associated with neurological diseases and congenital microcephaly, was first reported in the Western hemisphere in early 2015. A number of authors have reconstructed its epidemiological history using advanced phylogenetic approaches, and the majority of Zika phylogeography studies have used discrete diffusion models. Continuous space diffusion models make it possible to infer the possible origin of the virus in real space by reconstructing its ancestral location on the basis of geographical coordinates deduced from the latitude and longitude of the sampling locations. We analysed all the ZIKV complete genome isolates whose sampling times and localities were available in public databases at the time the study began, using a Bayesian approach for discrete and continuous phylogeographic reconstruction. The discrete phylogeographic analysis suggested that ZIKV emerged to become endemic/epidemic in the first decade of the 1900s in the Ugandan rainforests, and then reached Western Africa and Asia between the 1930s and 1950s. After a long period of about 40 years, it spread to the Pacific islands and reached Brazil from French Polynesia. Continuous phylogeography of the American epidemic showed that the virus entered in north-eastern Brazil in late 2012 and started to spread in early 2013 from two high probability regions: one corresponding to the entire north-east Brazil and the second surrounding the city of Rio de Janeiro, in a mainly northwesterly direction to Central America, the north-western countries of south America and the Caribbean islands. Our data suggest its cryptic circulation in both French Polynesia and Brazil, thus raising questions about the mechanisms underlying its undetected persistence in the absence of a known animal reservoir, and underline the importance of continuous diffusion models in making more reliable phylogeographic reconstructions of emerging viruses.

High-resolution diffusion pattern of human infections by Salmonella enterica serovar Napoli in Northern Italy explained through phylogeography.

Salmonella enterica serovar Napoli (serovar Napoli) is an emerging cause of human salmonellosis in Northern Italy. No specific reservoirs of serovar Napoli have been identified in Italy, so far. However, the environment, especially surface waters, has been hypothesized as an important source of infection based on the observation that genotypically different clusters of serovar Napoli are detected in different geographical macro-areas. To further support the hypothesis of a spatially-restricted pattern of serovar Napoli diffusion, a spatial segregation of serovar Napoli lineages should be observed also at smaller geographical scale. However, classical genotyping techniques used for Salmonella, such as pulsed-field gel electrophoresis (PFGE), did not possess enough discriminatory power to highlight spatial clustering of serovar Napoli within the macro-areas. To this purpose, we performed phylogeographical analyses based on genome-wide single nucleotide polymorphisms to test whether spatio-temporal evolution patterns of serovar Napoli in Northern Italy could be recognized with high geographical resolution, i.e. at local level. Specifically, we analyzed the local spread of the main PFGE clonal group, responsible for more than 60% of human infections in the study area, that did not show any geographical differentiation by PFGE within Northern Italy, i.e. the macro-area considered in the study. Both discrete and continuous phylogeography highlighted the existence of two main geographically-restricted clades: a Southern clade corresponding to the Po Valley and a Northern clade corresponding to the Pre-Alps area. Furthermore, the phylogeographical analyses suggested that the most probable site of origin of the clone was in an area of the Po Valley at the confluence of the Po and Ticino rivers, one of the most important Italian wetlands. These findings provide further support to the hypothesis that environmental transmission may play an important role in the ecology of serovar Napoli.

Highlighting priority areas for bovine viral diarrhea control in Italy: A phylogeographic approach.

The prevalence and genetic diversity of bovine viral diarrhea virus (BVDV) in a geographic area are largely influenced by live animal trade and management practices. Despite control and eradication programs currently underway in several European countries, the risk of BVDV spread within and among countries is still present. BVDV-1 is the predominant type circulating in European cattle population. In this study, a phylogeographic analysis was applied to the BVDV-1 highest prevalent subtypes in Italy to reconstruct the origin and spatial-temporal distribution and to trace main viral flows between different locations to highlight priority areas for BVDV control. A comprehensive dataset of BVDV-1b (n = 173) and 1e (n = 172) 5' UTR sequences was analysed, including both novel and published sequences from Italy and from European countries bordering and/or with commercial cattle flows with Italy. A common phylogeographic pattern was observed for BVDV-1b and 1e subtypes: interspersion from multiple Italian areas and European countries was widespread until the end of the last century, whereas significant local clusters were observed starting from 2000. These findings support a continuous viral flow among different areas over long time scales with no evidence of significant geographical structure, while local transmission networks are limited to more recent years. Northern Italy has been confirmed as the area of origin of the main clades of both BVDV subtypes at national level, acting both as a crucial area for introduction and a maintenance source for other areas. Piedmont, Central and Southern Italian regions contributed to limited geographical distribution and local BVDV-1b and 1e persistence. On the whole, priority control measures for BVDV-1b and 1e in Italy should be focused on: i) implementation of BVDV systematic control in all Northern Italian regions to break the viral flow from larger to smaller animal populations; and ii) breaking the dynamics of infections in regions with self-maintenance of BVDV by voluntary control programs.

Phylogeny of human T cell lymphotropic virus type 1 in Peru: a high degree of evolutionary relatedness with South African isolates.

We investigated the prevalence and molecular epidemiology of human T cell lymphotropic virus type 1 in Peruvian HIV-1-positive subjects, and found a 10.1% prevalence in a consecutive series of 318 HIV-1-positive patients living in Lima. Phylogenetic analysis of the long terminal repeat of 10 patient isolates showed that all of them belonged to the HTLV-1aA (Transcontinental) subgroup. Although the majority of the Peruvian sequences included in the analysis formed a clade with other Latin American sequences, the isolates of three patients clustered significantly with South African strains. These data show a high prevalence of HTLV-1 infection in HIV-1-positive subjects living in Lima and confirm the presence in Latin America of HTLV-1 strains probably arising from South Africa.

Reconstructing the recent West Nile virus lineage 2 epidemic in Europe and Italy using discrete and continuous phylogeography.

West Nile virus lineage 2 (WNV-2) was mainly confined to sub-Saharan Africa until the early 2000s, when it was identified for the first time in Central Europe causing outbreaks of human and animal infection. The aim of this study was to reconstruct the origin and dispersion of WNV-2 in Central Europe and Italy on a phylodynamic and phylogeographical basis. To this aim, discrete and continuous space phylogeographical models were applied to a total of 33 newly characterised full-length viral genomes obtained from mosquitoes, birds and humans in Northern Italy in the years 2013-2015 aligned with 64 complete sequences isolated mainly in Europe. The European isolates segregated into two highly significant clades: a small one including three sequences and a large clade including the majority of isolates obtained in Central Europe since 2004. Discrete phylogeographical analysis showed that the most probable location of the root of the largest European clade was in Hungary a mean 12.78 years ago. The European clade bifurcated into two highly supported subclades: one including most of the Central/East European isolates and the other encompassing all of the isolates obtained in Greece. The continuous space phylogeographical analysis of the Italian clade showed that WNV-2 entered Italy in about 2008, probably by crossing the Adriatic sea and reaching a central area of the Po Valley. The epidemic then spread simultaneously eastward, to reach the region of the Po delta in 2013, and westward to the border area between Lombardy and Piedmont in 2014; later, the western strain changed direction southward, and reached the central area of the Po valley once again in 2015. Over a period of about seven years, the virus spread all over an area of northern Italy by following the Po river and its main tributaries.

Distribution of Marburg virus in Africa: An evolutionary approach.

The aim of this study was to investigate the origin and geographical dispersion of Marburg virus, the first member of the Filoviridae family to be discovered. Seventy-three complete genome sequences of Marburg virus isolated from animals and humans were retrieved from public databases and analysed using a Bayesian phylogeographical framework. The phylogenetic tree of the Marburg virus data set showed two significant evolutionary lineages: Ravn virus (RAVV) and Marburg virus (MARV). MARV divided into two main clades; clade A included isolates from Uganda (five from the European epidemic in 1967), Kenya (1980) and Angola (from the epidemic of 2004-2005); clade B included most of the isolates obtained during the 1999-2000 epidemic in the Democratic Republic of the Congo (DRC) and a group of Ugandan isolates obtained in 2007-2009. The estimated mean evolutionary rate of the whole genome was 3.3×10(-4) substitutions/site/year (credibility interval 2.0-4.8). The MARV strain had a mean root time of the most recent common ancestor of 177.9years ago (YA) (95% highest posterior density 87-284), thus indicating that it probably originated in the mid-XIX century, whereas the RAVV strain had a later origin dating back to a mean 33.8 YA. The most probable location of the MARV ancestor was Uganda (state posterior probability, spp=0.41), whereas that of the RAVV ancestor was Kenya (spp=0.71). There were significant migration rates from Uganda to the DRC (Bayes Factor, BF=42.0) and in the opposite direction (BF=5.7). Our data suggest that Uganda may have been the cradle of Marburg virus in Africa.

Spatial and Temporal Phylogeny of Border Disease Virus in Pyrenean Chamois (Rupicapra p. pyrenaica).

Border disease virus (BDV) affects a wide range of ruminants worldwide, mainly domestic sheep and goat. Since 2001 several outbreaks of disease associated to BDV infection have been described in Pyrenean chamois (Rupicapra pyrenaica pyrenaica) in Spain, France and Andorra. In order to reconstruct the most probable places of origin and pathways of dispersion of BDV among Pyrenean chamois, a phylogenetic analysis of 95 BDV 5'untranslated sequences has been performed on chamois and domestic ungulates, including novel sequences and retrieved from public databases, using a Bayesian Markov Chain Monte Carlo method. Discrete and continuous space phylogeography have been applied on chamois sequences dataset, using centroid positions and latitude and longitude coordinates of the animals, respectively. The estimated mean evolutionary rate of BDV sequences was 2.9×10-3 subs/site/year (95% HPD: 1.5-4.6×10-3). All the Pyrenean chamois isolates clustered in a unique highly significant clade, that originated from BDV-4a ovine clade. The introduction from sheep (dated back to the early 90s) generated a founder effect on the chamois population and the most probable place of origin of Pyrenean chamois BDV was estimated at coordinates 42.42 N and 1.9 E. The pathways of virus dispersion showed two main routes: the first started on the early 90s of the past century with a westward direction and the second arise in Central Pyrenees. The virus spread westward for more than 125 km and southward for about 50km and the estimated epidemic diffusion rate was about 13.1 km/year (95% HPD 5.2-21.4 km/year). The strong spatial structure, with strains from a single locality segregating together in homogeneous groups, and the significant pathways of viral dispersion among the areas, allowed to reconstruct both events of infection in a single area and of migrations, occurring between neighboring areas.

Dating the origin and dispersal of Human Papillomavirus type 16 on the basis of ancestral human migrations.

A major limitation when reconstructing the origin and evolution of HPV-16 is the lack of reliable substitution rate estimates for the viral genes. On the basis of the hypothesis of human HPV-16 co-divergence, we estimated a mean evolutionary rate of 1.47×10(-7) (95% HPD=0.64-2.47×10(-7)) subs/site/year for the viral LCR region. The results of a Bayesian phylogeographical analysis suggest that the currently circulating HPV-16 most probably originated in Africa about 110 thousand years ago (Kya), before giving rise to four known geographical lineages: the Asian/European lineage, which most probably originated in Asia a mean 38 Kya, and the Asian/American and two African lineages, which probably respectively originated about 33 and 27 Kya. These data closely reflect current hypotheses concerning modern human expansion based on studies of mitochondrial DNA phylogeny. The correlation between ancient human migration and the present HPV phylogeny may be explained by the co-existence of modes of transmission other than sexual transmission.

Phylogeography, phylodynamics and transmission chains of bovine viral diarrhea virus subtype 1f in Northern Italy.

Bovine viral diarrhea virus (BVDV) type 1 in Italy is characterized by high genetic diversity, with at least 20 subtypes. Subtype 1f is endemic in a restricted geographic area, meaning that it has local distribution. We investigated the population dynamics of BVDV-1f in Northern Italy and characterized the transmission chains of a subset of samples from Piedmont and Aosta Valley regions. A total of 51 samples from 1966 to 2013 were considered and 5' UTR sequences were used for phylogeography. A subset of 12 samples was selected for Npro gene sequencing and further characterization of the transmission chains using both molecular and epidemiological data. Phylogeography estimated the root of BVDV-1f tree in Veneto in 1965. Four significant subclades included sequences clustering by region: Lombardy (n=3), Lombardy and Emilia-Romagna (n=7), Piedmont (n=17), Piedmont and Aosta Valley (n=21). The Piedmont-only subclade has a ladder-like branching structure, while the Piedmont and Aosta Valley subclade has a nearly complete binary structure. In the subset, the outbreak reconstruction identified one sample from Piedmont as the most probable source of infection for the Aosta Valley cases. An ad hoc questionnaire submitted to public veterinarians revealed connections between sampled and non-sampled farms by means of trades, exhibitions and markets. According to the phylogeography, BVDV-1f moved westward, entering from Veneto, and spreading to Lombardy and Emilia-Romagna in the early 1990s, and finally to Piedmont and Aosta Valley in the first decade of 2000s. Both phylogeographic analyses on the whole dataset and on the selection of Npro dataset pointed out that subtype 1f entered Aosta Valley from Piedmont. The integration of molecular and epidemiological data revealed connections between farms, and such approach should be considered in any control plan. In Aosta Valley, the study showed that BVDV1f can be controlled only monitoring the introduction of cattle from Piedmont region.