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.9x10(-3) subs/site/year (95% HPD: 1.5-4.6x10(-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

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

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

Sofosbuvir Selects for Drug-Resistant Amino Acid Variants in the Zika Virus RNA-Dependent RNA-Polymerase Complex In Vitro

The nucleotide analog sofosbuvir, licensed for the treatment of hepatitis C, recently revealed activity against the Zika virus (ZIKV) in vitro and in animal models. However, the ZIKV genetic barrier to sofosbuvir has not yet been characterized. In this study, in vitro selection experiments were performed in infected human hepatoma cell lines. Increasing drug pressure significantly delayed viral breakthrough (p = 0.029). A double mutant in the NS5 gene (V360L/V607I) emerged in 3 independent experiments at 40–80 µM sofosbuvir resulting in a 3.9 ± 0.9-fold half- maximal inhibitory concentration (IC50) shift with respect to the wild type (WT) virus. A triple mutant (C269Y/V360L/V607I), detected in one experiment at 80 µM, conferred a 6.8-fold IC50 shift with respect to the WT. Molecular dynamics simulations confirmed that the double mutant V360L/V607I impacts the binding mode of sofosbuvir, supporting its role in sofosbuvir resistance. Due to the distance from the catalytic site and to the lack of reliable structural data, the contribution of C269Y was not investigated in silico. By a combination of sequence analysis, phenotypic susceptibility testing, and molecular modeling, we characterized a double ZIKV NS5 mutant with decreased sofosbuvir susceptibility. These data add important information to the profile of sofosbuvir as a possible lead for anti-ZIKV drug development

Sofosbuvir Selects for Drug-Resistant Amino Acid Variants in the Zika Virus RNA-Dependent RNA-Polymerase Complex In Vitro

The nucleotide analog sofosbuvir, licensed for the treatment of hepatitis C, recently revealed activity against the Zika virus (ZIKV) in vitro and in animal models. However, the ZIKV genetic barrier to sofosbuvir has not yet been characterized. In this study, in vitro selection experiments were performed in infected human hepatoma cell lines. Increasing drug pressure significantly delayed viral breakthrough (p = 0.029). A double mutant in the NS5 gene (V360L/V607I) emerged in 3 independent experiments at 40–80 µM sofosbuvir resulting in a 3.9 ± 0.9-fold half- maximal inhibitory concentration (IC50) shift with respect to the wild type (WT) virus. A triple mutant (C269Y/V360L/V607I), detected in one experiment at 80 µM, conferred a 6.8-fold IC50 shift with respect to the WT. Molecular dynamics simulations confirmed that the double mutant V360L/V607I impacts the binding mode of sofosbuvir, supporting its role in sofosbuvir resistance. Due to the distance from the catalytic site and to the lack of reliable structural data, the contribution of C269Y was not investigated in silico. By a combination of sequence analysis, phenotypic susceptibility testing, and molecular modeling, we characterized a double ZIKV NS5 mutant with decreased sofosbuvir susceptibility. These data add important information to the profile of sofosbuvir as a possible lead for anti-ZIKV drug development

Phylogenetic analysis of human immunodeficiency virus type 2 group B

Context: Human immunodeficiency virus type 2 (HIV-2) infections are mainly restricted to West Africa; however, in the recent years, the prevalence of HIV-2 is a growing concern in some European countries and the Southwestern region of India. Despite the presence of different HIV-2 groups, only A and B Groups have established human-to-human transmission chains. Aims: This work aimed to evaluate the phylogeographic inference of HIV-2 Group B worldwide to estimate their data of origin and the population dynamics. Materials and Methods: The evolutionary rates, the demographic history for HIV-2 Group B dataset, and the phylogeographic analysis were estimated using a Bayesian approach. The viral gene flow analysis was used to count viral gene out/in flow among different locations. Results: The root of the Bayesian maximum clade credibility tree of HIV-2 Group B dated back to 1957. The demographic history of HIV-2 Group B showed that the epidemic remained constant up to 1970 when started an exponential growth. From 1985 to early 2000s, the epidemic reached a plateau, and then it was characterized by two bottlenecks and a new plateau at the end of 2000s. Phylogeographic reconstruction showed that the most probable location for the root of the tree was Ghana. Regarding the viral gene flow of HIV-2 Group B, the only observed viral gene flow was from Africa to France, Belgium, and Luxembourg. Conclusions: The study gives insights into the origin, history, and phylogeography of HIV-2 Group B epidemic. The growing number of infections of HIV-2 worldwide indicates the need for strengthening surveillance

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

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