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.

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.

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.

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.

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.

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.

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.

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.

Reconstruction of the Evolutionary Dynamics of A(H3N2) Influenza Viruses Circulating in Italy from 2004 to 2012.

BACKGROUND: Influenza A viruses are characterised by their rapid evolution, and the appearance of point mutations in the viral hemagglutinin (HA) domain causes seasonal epidemics. The A(H3N2) virus has higher mutation rate than the A(H1N1) virus. The aim of this study was to reconstruct the evolutionary dynamics of the A(H3N2) viruses circulating in Italy between 2004 and 2012 in the light of the forces driving viral evolution. METHODS: Phylodinamic analyses were made using a Bayesian method, and codon-specific positive selection acting on the HA coding sequence was evaluated. RESULTS: Global and local phylogenetic analyses showed that the Italian strains collected between 2004 and 2012 grouped into five significant Italian clades that included viral sequences circulating in different epidemic seasons. The time of the most recent common ancestor (tMRCA) of the tree root was between May and December 2003. The tMRCA estimates of the major clades suggest that the origin of a new viral strain precedes the effective circulation of the strain in the Italian population by 6-31 months, thus supporting a central role of global migration in seeding the epidemics in Italy. The study of selection pressure showed that four codons were under positive selection, three of which were located in antigenic sites. Analysis of population dynamics showed the alternation of periods of exponential growth followed by a decrease in the effective number of infections corresponding to epidemic and inter-epidemic seasons. CONCLUSIONS: Our analyses suggest that a complex interaction between the immune status of the population, migrations, and a few selective sweeps drive the influenza A(H3N2) virus evolution. Our findings suggest the possibility of the year-round survival of local strains even in temperate zones, a hypothesis that warrants further investigation.

Reliable timescale inference of HBV genotype A origin and phylodynamics.

The worldwide distributed Hepatitis B virus (HBV) genotype A is classified into three subgenotypes, and one quasi-subgenotype. The majority of HBV-A subgenotypes are widespread in Africa and in ethnic groups that have relatively recently emigrated from African countries, whereas HBV-A2 is highly prevalent among subjects at high risk for sexual exposure to HBV in north-western Europe and the USA. The aim of this study was to reconstruct the origin and dispersion of HBV-A subgenotypes on a reliable timescale using short-term calibration based on heterochronous sampling for HBV-A2, and long-term calibration based on historical data for the other subgenotypes. To this aim, we analysed 113 newly characterised HBV-A isolates with 247 reference sequences retrieved from a public database. The phylodynamic reconstruction was performed by a Bayesian framework. The common ancestor of the currently circulating A subgenotypes was placed in west-central Africa a mean 1057 years ago. The genotype diverged into two main clades at the beginning of the 13th century: one including all of the west-central African quasi-subgenotypes and the other corresponding to subgenotype A1, originating in east Africa and further segregating into two main subclades: an "African" and a "cosmopolitan" clade. It is likely that the slave trade was the main source the spread of cosmopolitan HBV-A1, which was exported to Asia in the 17th century as a result of Arab or Portuguese trade, and to Latin America in the 18th centuries through the trans-Atlantic slave trade. The origin of the currently circulating A2 strains dates back to the first decades of the 20th century, and the evolutionary demography analysis suggests an exponential growth of infections, between 1970s and the mid-1990s. In conclusion, the very different epidemiological and evolutionary histories of HBV-A subgenotypes justify the use of different calibration approaches to reconstruct their reciprocal phylodynamics.

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.

Extended genetic diversity of bovine viral diarrhea virus and frequency of genotypes and subtypes in cattle in Italy between 1995 and 2013.

Genetic typing of bovine viral diarrhea virus (BVDV) has distinguished BVDV-1 and BVDV-2 species and an emerging putative third species (HoBi-like virus), recently detected in southern Italy, signaling the occurrence of natural infection in Europe. Recognizing the need to update the data on BVDV genetic variability in Italy for mounting local and European alerts, a wide collection of 5' UTR sequences (n = 371) was selected to identify the frequency of genotypes and subtypes at the herd level. BVDV-1 had the highest frequency, followed by sporadic BVDV-2. No novel HoBi-like viruses were identified. Four distribution patterns of BVDV-1 subtypes were observed: highly prevalent subtypes with a wide temporal-spatial distribution (1b and 1e), low prevalent subtypes with a widespread geographic distribution (1a, 1d, 1g, 1h, and 1k) or a restricted geographic distribution (1f), and sporadic subtypes detected only in single herds (1c, 1j, and 1l). BVDV-1c, k, and l are reported for the first time in Italy. A unique genetic variant was detected in the majority of herds, but cocirculation of genetic variants was also observed. Northern Italy ranked first for BVDV introduction, prevalence, and dispersion. Nevertheless, the presence of sporadic variants in other restricted areas suggests the risk of different routes of BVDV introduction.

Enigmatic origin of hepatitis B virus: an ancient travelling companion or a recent encounter?

Hepatitis B virus (HBV) is the leading cause of liver disease and infects an estimated 240 million people worldwide. It is characterised by a high degree of genetic heterogeneity because of the use of a reverse transcriptase during viral replication. The ten genotypes (A-J) that have been described so far further segregate into a number of subgenotypes which have distinct ethno-geographic distribution. Genotypes A and D are ubiquitous and the most prevalent genotypes in Europe (mainly represented by subgenotypes D1-3 and A2); genotypes B and C are restricted to eastern Asia and Oceania; genotype E to central and western Africa; and genotypes H and F (classified into 4 subgenotypes) to Latin America and Alaska. This review summarises the data obtained by studying the global phylodynamics and phylogeography of HBV genotypes, particularly those concerning the origin and dispersion histories of genotypes A, D, E and F and their subgenotypes. The lack of any consensus concerning the HBV substitution rate and the conflicting data obtained using different calibration approaches make the time of origin and divergence of the various genotypes and subgenotypes largely uncertain. It is hypothesised that HBV evolutionary rates are time dependent, and that the changes depend on the main transmission routes of the genotypes and the dynamics of the infected populations.

Phylogeny and population dynamics of respiratory syncytial virus (Rsv) A and B.

Respiratory syncytial virus (RSV) is a major cause of lower respiratory tract infections in infants and young children. RSV is characterised by high variability, especially in the G glycoprotein, which may play a significant role in RSV pathogenicity by allowing immune evasion. To reconstruct the origin and phylodynamic history of RSV, we evaluated the genetic diversity and evolutionary dynamics of RSV A and RSV B isolated from children under 3 years old infected in Italy from 2006 to 2012. Phylogenetic analysis revealed that most of the RSV A sequences clustered with the NA1 genotype, and RSV B sequences were included in the Buenos Aires genotype. The mean evolutionary rates for RSV A and RSV B were estimated to be 2.1 × 10(-3) substitutions (subs)/site/year and 3.03 × 10(-3) subs/site/year, respectively. The time of most recent common ancestor for the tree root went back to the 1940s (95% highest posterior density-HPD: 1927-1951) for RSV A and the 1950s (95%HPD: 1951-1960) for RSV B. The RSV A Bayesian skyline plot (BSP) showed a decrease in transmission events ending in about 2005, when a sharp growth restored the original viral population size. RSV B BSP showed a similar trend. Site-specific selection analysis identified 10 codons under positive selection in RSV A sequences and only one site in RSV B sequences. Although RSV remains difficult to control due to its antigenic diversity, it is important to monitor changes in its coding sequences, to permit the identification of future epidemic strains and to implement vaccine and therapy strategies.