Genetic diversity in the env V1-V2 region of proviral quasispecies from long-term controller MHC-typed cynomolgus macaques infected with SHIVSF162P4cy.

Intra-host evolution of human immunodeficiency virus (HIV) and simian immunodeficiency virus (SIV) has been shown by viral RNA analysis in subjects who naturally suppress plasma viremia to low levels, known as controllers. However, little is known about the variability of proviral DNA and the inter-relationships among contained systemic viremia, rate of reservoir reseeding and specific major histocompatibility complex (MHC) genotypes, in controllers. Here, we analysed the proviral DNA quasispecies of the env V1-V2 region, in PBMCs and in anatomical compartments of 13 long-term controller monkeys after 3.2 years of infection with simian/human immunodeficiency virus (SHIV)SF162P4cy. A considerable variation in the genetic diversity of proviral quasispecies was present among animals. Seven monkeys exhibited env V1-V2 proviral populations composed of both clusters of identical ancestral sequences and new variants, whereas the other six monkeys displayed relatively high env V1-V2 genetic diversity with a large proportion of diverse novel sequences. Our results demonstrate that in SHIVSF162P4cy-infected monkeys there exists a disparate pattern of intra-host viral diversity and that reseeding of the proviral reservoir occurs in some animals. Moreover, even though no particular association has been observed between MHC haplotypes and the long-term control of infection, a remarkably similar pattern of intra-host viral diversity and divergence was found within animals carrying the M3 haplotype. This suggests that in animals bearing the same MHC haplotype and infected with the same virus, viral diversity follows a similar pattern with similar outcomes and control of infection.

Klebsiella pneumoniae blaKPC-3 nosocomial epidemic: Bayesian and evolutionary analysis.

K. pneumoniae isolates carrying blaKPC-3 gene were collected to perform Bayesian phylogenetic and selective pressure analysis and to apply homology modeling to the KPC-3 protein. A dataset of 44 blakpc-3 gene sequences from clinical isolates of K. pneumoniae was used for Bayesian phylogenetic, selective pressure analysis and homology modeling. The mean evolutionary rate for blakpc-3 gene was 2.67×10-3 substitution/site/year (95% HPD: 3.4×10-4-5.59×10-3). The root of the Bayesian tree dated back to the year 2011 (95% HPD: 2007-2012). Two main clades (I and II) were identified. The population dynamics analysis showed an exponential growth from 2011 to 2013 and the reaching of a plateau. The phylogeographic reconstruction showed that the root of the tree had a probable common ancestor in the general surgery ward. Selective pressure analysis revealed twelve positively selected sites. Structural analysis of KPC-3 protein predicted that the amino acid mutations are destabilizing for the protein and could alter the substrate specificity. Phylogenetic analysis and homology modeling of blaKPC-3 gene could represent a useful tool to follow KPC spread in nosocomial setting and to evidence amino acid substitutions altering the substrate specificity.

Phylogenesys and homology modeling in Zika virus epidemic: food for thought.

Zika virus (ZIKV) is an emerging Flavivirus that have recently caused an outbreak in Brazil and rapid spread in several countries. In this study, the consequences of ZIKV evolution on protein recognition by the host immune system have been analyzed. Evolutionary analysis was combined with homology modeling and T-B cells epitope predictions. Two separate clades, the African one with the Uganda sequence, as the most probable ancestor, and the second one containing all the most recent sequences from the equatorial belt were identified. Brazilian strains clustered all together and closely related to the French Polynesia isolates. A strong presence of a negatively selected site in the envelope gene (Env) protein was evidenced, suggesting a probable purging of deleterious polymorphisms in functionally important genes. Our results show relative conservancy of ZIKV sequences when envelope and other non-structural proteins (NS3 and NS5) are analyzed by homology modeling. However, some regions within the consensus sequence of NS5 protein and to a lesser extent in the envelope protein, show localized high mutation frequency corresponding to a considerable alteration in protein stability. In terms of viral immune escape, envelope protein is under a higher selective pressure than NS5 and NS3 proteins for HLA class I and II molecules. Moreover, envelope mutations that are not strictly related to T-cell immune responses are mostly located on the surface of the protein in putative B-cell epitopes, suggesting an important contribution of B cells in the immune response as well.

Molecular epidemiology, evolution and phylogeny of Chikungunya virus: An updating review.

Chikungunya virus (CHIKV) is a mosquito-transmitted alphavirus belonging to the Togaviridae family, causing a febrile illness associated with severe arthralgia and rash. In this review, we summarized a series of articles published from 2013 to 2016 concerning CHIKV epidemiology, phylogeny, vaccine and therapies, to give an update of our most recent article written in 2014 (Lo Presti et al.,2014). CHIKV infection was first reported in 1952 from Makonde plateaus and since this time caused many outbreaks worldwide, involving the Indian Ocean region, African countries, American continent and Italy. CHIKV infection is still underestimated and it is normally associated with clinical symptoms overlapping with dengue virus, recurring epidemics and mutations within the viral genome. These characteristics promote the geographical spread and the inability to control vector-mediated transmission of the virus. For these reasons, the majority of studies were aimed to describe outbreaks and to enhance knowledge on CHIKV biology, pathogenesis, infection treatment, and prevention. In this review, 16 studies on CHIKV phylogenetic and phylodinamics were considered, during the years 2013-2016. Phylogenetic and phylodinamic analysis are useful tools to investigate how the genealogy of a pathogen population is influenced by pathogen's demographic history, host immunological milieu and environmental/ecological factors. Phylogenetic tools were revealed important to reconstruct the geographic spread of CHIKV during the epidemics wave and to have information on the circulating strains of the virus, that are important for the prediction and control of the epidemics, as well as for vaccines and antiviral drugs development. In conclusion, this updating review can give a critical appraisal of the epidemiology, therapeutic and phylogenesis of CHIKV, reinforcing the need to monitor the geographic spread of virus and vectors.

Procalcitonin, MR-Proadrenomedullin, and Cytokines Measurement in Sepsis Diagnosis: Advantages from Test Combination.

BACKGROUND: Elevated cytokines levels correlate with sepsis severity and mortality but their role in the diagnosis is controversial, whereas Procalcitonin (PCT) has been largely used. Recently, the mid-regional proadrenomedullin (MR-proADM) has been combined with PCT for diagnosis optimization. In this study the combined measurement of PCT, MR-proADM, and cytokines in patients with sepsis was evaluated. METHODS: One hundred and four septic patients and 101 controls were enrolled. Receiver operating characteristic (ROC) analysis and multiple logistic regression were used to evaluate applicant markers for sepsis diagnosis. Markers with best Odds Ratio (OR) were combined, and the posttest probability and a composite score were computed. RESULTS: Based upon ROC curves analysis, PCT, MR-proADM, IL-6, IL-10, TNF-α, and MCP-1 were considered applicant for sepsis diagnosis. Among these PCT, MR-proADM , IL-6, and TNF-α showed the best OR. A better posttest probability was found with the combination of PCT with MR-proADM and PCT with IL-6 or TNF-α compared to the single marker. A composite score of PCT, MR-proADM, and TNF-α showed the best ROC curve in the early diagnosis of sepsis. CONCLUSION: The combination of PCT with other markers should expedite diagnosis and treatment of sepsis optimizing clinical management.

Impact of spatial dispersion, evolution, and selection on Ebola Zaire Virus epidemic waves.

Ebola virus Zaire (EBOV) has reemerged in Africa, emphasizing the global importance of this pathogen. Amidst the response to the current epidemic, several gaps in our knowledge of EBOV evolution are evident. Specifically, uncertainty has been raised regarding the potential emergence of more virulent viral variants through amino acid substitutions. Glycoprotein (GP), an essential component of the EBOV genome, is highly variable and a potential site for the occurrence of advantageous mutations. For this study, we reconstructed the evolutionary history of EBOV by analyzing 65 GP sequences from humans and great apes over diverse locations across epidemic waves between 1976 and 2014. We show that, although patterns of spatial dispersion throughout Africa varied, the evolution of the virus has largely been characterized by neutral genetic drift. Therefore, the radical emergence of more transmissible variants is unlikely, a positive finding, which is increasingly important on the verge of vaccine deployment.

Amino acid mutations in Ebola virus glycoprotein of the 2014 epidemic.

Zaire Ebola virus (EBOV) is an enveloped non-segmented negative strand RNA virus of 19 kb in length belonging to the family Filoviridae. The virus was isolated and identified in 1976 during the epidemic of hemorrhagic fever in Zaire. The most recent outbreak of EBOV among humans, was that occurred in the forested areas of south eastern Guinea, that began in February 2014 and is still ongoing. The recent Ebola outbreak, is affecting other countries in West Africa, in addiction to Guinea: Liberia, Nigeria, and Sierra Leone. In this article, a selective pressure analysis and homology modeling based on the G Glycoprotein (GP) sequences retrieved from public databases were used to investigate the genetic diversity and modification of antibody response in the recent outbreak of Ebola Virus. Structural and the evolutionary analysis underline the 2014 epidemic virus being under negative selective pressure does not change with respect to the old epidemic in terms of genome adaptation.

Analysis of the ORFK1 hypervariable regions reveal distinct HHV-8 clustering in Kaposi's sarcoma and non-Kaposi's cases.

BACKGROUND: Classical Kaposi's Sarcoma (cKS) is a rare vascular tumor, which develops in subjects infected with Human Herpesvirus-8 (HHV-8). Beside the host predisposing factors, viral genetic variants might possibly be related to disease development. The aim of this study was to identify HHV-8 variants in patients with cKS or in HHV-8 infected subjects either asymptomatic or with cKS-unrelated cutaneous lymphoproliferative disorders. METHODS: The VR1 and VR2 regions of the ORF K1 sequence were analyzed in samples (peripheral blood and/or lesional tissue) collected between 2000 and 2010 from 27 subjects with HHV-8 infection, established by the presence of anti-HHV-8 antibodies. On the basis of viral genotyping, a phylogenetic analysis and a time-scaled evaluation were performed. RESULTS: Two main clades of HHV-8, corresponding to A and C subtypes, were identified. Moreover, for each subtype, two main clusters were found distinctively associated to cKS or non-cKS subjects. Selective pressure analysis showed twelve sites of the K1 coding gene (VR1 and VR2 regions) under positive selective pressure and one site under negative pressure. CONCLUSION: Thus, present data suggest that HHV-8 genetic variants may influence the susceptibility to cKS in individuals with HHV-8 infection.

Back to the origin of HCV 2c subtype and spreading to the Calabria region (Southern Italy) over the last two centuries: a phylogenetic study.

Circulation of HCV genotype 2 has been described in European Countries where numerous subtypes and unclassified HCV 2 lineages have been reported. In Italy, subtype 1b is the most prevalent, followed by genotype 2. In the present study, phylogeny of HCV 2c was investigated. The phylogeny of HCV 2c isolated from 54 Italian patients in the Calabria region (Southern Italy) was investigated by analyzing a fragment of the NS5B gene. Patients came from 5 metropolitan areas and a small village (Sersale). These areas were geographically dispersed throughout the entire region. A Bayesian coalescent-based framework was used to estimate origin and spreading of HCV 2c in this region. Phylogenetic analysis showed that 28 Italian sequences were intermixed with foreign HCV 2c reference sequences and grouped into 3 major clades: A, B, and C. Nineteen inter-clade sequences were associated uniquely with surgery as risk factor for HCV acquisition. By contrast, a sub-cluster within clade B was associated with blood transfusion. Moreover, sequences from Sersale village grouped in the Italian sub-cluster and were intermixed with 10 sequences from metropolitan areas. The three isolates with the longest branch came from Sersale and belonged to patients who had glass syringes as risk factor. HCV 2c isolates from the Calabria region shared a common ancestor whose origin was traced back to 1889. Our results suggest that, after its introduction - possibly as a result of population movements between Italy and African Countries during Italian colonialism - HCV 2c spread through multiple risk factors, not including intravenous drug use. So, transmission chains followed a pathway different from other European Countries. Although HCV incidence is decreasing, these ways are still ongoing, possibly justifying stability in the relative prevalence of HCV 2c.

The deduced evolution history of Omsk hemorrhagic fever virus.

Omsk hemorrhagic fever (OHF) is a severe disease that emerged in the 1940s in Siberia, Russia. It is caused by the OHF virus (OHFV), belonging to the Flavivirus genus. In wildlife, the principal vector of OHFV is the Dermacentor reticulatus tick. However, humans are mainly infected after contact with an infected muskrat Ondatra zibethicus. The evolutionary history of OHFV is not yet clarified. In an attempt to reconstruct the temporal and spatial phylodynamics of OHFV, a collection of 25 OHFV strains was studied. Maximum likelihood analysis, the Bayesian MCMC approach, and procedures based on continuous-time Markov Chain process, were used for the investigation of OHFV E gene nucleotide sequences. Six statistically supported clusters of OHFV were identified; five of them joined in a main clade A. The first revealed evolutionary event, when OHFV clade A and clade B divided, dated to about 700 years ago. Clusters C, D, and E, within clade A, separated 230 years ago and further evolved during last century. The phylogeographic analysis showed that clade A originated in the Omsk Province, whereas clusters B, C, and E appeared to originate in Kurgan, Novosibirsk, and Omsk Provinces, respectively. In conclusion, OHFV, as a member of the mammalian tick-borne group of flaviviruses, evolved in Western Siberia during the last millennium. When a highly susceptible species, O.zibethicus, was introduced into the region, in the 1930s, OHFV used this species as an amplifying host that lead to numerous fatal epizootics in muskrats and to human OHF outbreaks.

Epidemiological history and phylogeography of West Nile virus lineage 2.

West Nile virus (WNV) was first isolated in Uganda. In Europe WNV was sporadically detected until 1996, since then the virus has been regularly isolated from birds and mosquitoes and caused several outbreaks in horses and humans. Phylogenetic analysis showed two main different WNV lineages. The lineage 1 is widespread and segregates into different subclades (1a-c). WNV-1a includes numerous strains from Africa, America, and Eurasia. The spatio-temporal history of WNV-1a in Europe was recently described, identifying two main routes of dispersion, one in Eastern and the second in Western Europe. The West Nile lineage 2 (WNV-2) is mainly present in sub-Saharan Africa but has been recently emerged in Eastern and Western European countries. In this study we reconstruct the phylogeny of WNV-2 on a spatio-temporal scale in order to estimate the time of origin and patterns of geographical dispersal of the different isolates, particularly in Europe. Phylogeography findings obtained from E and NS5 gene analyses suggest that there were at least two separate introductions of WNV-2 from the African continent dated back approximately to the year 1999 (Central Europe) and 2000 (Russia), respectively. The epidemiological implications and clinical consequences of lineage 1 and 2 cocirculation deserve further investigations.

Human polyomaviruses identification by logic mining techniques.

BACKGROUND: Differences in genomic sequences are crucial for the classification of viruses into different species. In this work, viral DNA sequences belonging to the human polyomaviruses BKPyV, JCPyV, KIPyV, WUPyV, and MCPyV are analyzed using a logic data mining method in order to identify the nucleotides which are able to distinguish the five different human polyomaviruses. RESULTS: The approach presented in this work is successful as it discovers several logic rules that effectively characterize the different five studied polyomaviruses. The individuated logic rules are able to separate precisely one viral type from the other and to assign an unknown DNA sequence to one of the five analyzed polyomaviruses. CONCLUSIONS: The data mining analysis is performed by considering the complete sequences of the viruses and the sequences of the different gene regions separately, obtaining in both cases extremely high correct recognition rates.

Origin, evolution, and phylogeography of recent epidemic CHIKV strains.

Chikungunya virus (CHIKV) is an arthropod-borne virus of the Alphavirus genus, which is transmitted to humans by Aedes spp. mosquitoes and was firstly identified in Tanzania in the mid 1950s. In this article, the findings of a phylogenetic and phylogeographic analysis of the recent CHIKV pandemic are reported. We estimated time of origin of the ancestral virus, time and place of occurrence of A226V mutation, and the flow of viral strains from an area to the other. The Bayesian phylogenetic and phylogeographic analysis was performed on the whole dataset, which consisted of 195 E1 (envelope 1) CHIKV sequences, and on a subset (D2), including 146 of the 195 previous sequences. Using the relaxed clock model, we estimated a CHIKV E1 mean evolutionary rate (in the whole dataset) of 1.4 × 10(-3)substitution/site/year (95% highest posterior density interval HPD 6.4 × 10(-4)-2.5 × 10(-3)), and of 2.2 × 10(-3) (95% HPD 9.6 × 10(-4)-3.8 × 10(-3)) in the D2 subset, including only the strains involved in the recent Indian Ocean epidemic. The phylogeographical analysis suggested an African origin of CHIKV with a tMRCA of 146 years corresponding to 1863 (95% HPD 1741-1941). Moreover D2 subset most probably originated in Kenya, with a tMRCA corresponding to the year 2002 (95% HPD 2000-2004), then spread following two distinct routes: one throughout the Indian Ocean (Reunion, Comoros) and the other moving from India then scattered in the South East Asia and reached Italy. In conclusion, we reconstructed the geographic spread of CHIKV during the last epidemic wave, which showed an eastward path from Africa to Indian Ocean island to India, and from there to other South East Asian countries. Whether A226V variants followed the same migration path remains undefined, since local independent mutations, followed by fixation due to selective advantage conferred by better adaptation to local vectors of infection, cannot be excluded.

Human papillomaviruses 53 and 66: clinical aspects and genetic analysis.

Variants of HPV types may have different oncogenic potential. While HPV 16 and 18 variants have been extensively studied, little is known on the less frequent high-risk types such as HPV 53 and 66. Here, we analyzed the genetic variability of HPV 53 and 66 by sequencing the E6, E7, L1 genes and the Long Control Region (LCR) sequences of HPV 53 and HPV 66 from infected women. Fisher's exact-test was performed to correlate viral variants with cervical lesions. Higher-order interactions among identified mutations were analyzed by co-variation and cluster analyses. Antigenic-index alterations following L1 mutations were predicted by Jameson-Wolf algorithm. In HPV53, novel variants were identified in L1 (N=9) and E6 (N=1) genes. The novel L1 mutation P432L was statistically associated with L-SIL lesions (P=0.04) and its development reduced the L1 predicted antigenicity (up to -2.3 for Glu433). HPV 53 E6 and L1 sequences clustered phylogenetically into two main clades. In HPV 66, novel polymorphisms were identified in L1 (N=4) and E6 (N=4) genes. The L1 protein mutations S405P and D458N were exclusively found in patients with L-SILs. Seven E7 variants and 10 LCR variants were for the first time analyzed. Novel HPV 53 and 66 variants were identified in this study. Some of these mutations were significantly associated with L-SIL lesions and affected the antigenic index of the L1 protein with possible interesting implications in vaccine design.

Genetic variability and circulation pattern of human metapneumovirus isolated in Italy over five epidemic seasons.

Human metapneumovirus (hMPV) is an important aetiological agent of respiratory tract infection (RTI) in infants. Based on genetic differences, hMPVs are separated into two main groups, A and B, and two subgroups 1 and 2. To better understand the genetic diversity of hMPV, we analyzed 435 bp fragments of the F gene from 49 isolates obtained from a paediatric clinical centre in Northern Italy from 2005 to 2009. The phylogenetic analysis showed that our hMPV sequences clustered into five main clades (A1, A2a, A2b, B1, B2) statistically supported. Most of the strains belong to A2 (49%) and to B1 (28.5%) lineages. The intermixing in the phylogenetic tree showed no seasonal distribution between samples collected over a 5 year period. Mean genetic distances showed that clade A2 was more heterogeneous than others. In our F glycoprotein dataset we observed only two positively selected sites with an w value of 1.408 and 1.429, respectively, and 27 negatively selected sites. The two positively selected sites could be considered evolutionary "hot spots" because they were under positive selection and/or relaxed selective constraints. The abundant negatively selected sites reflect a high degree of conservation of F protein, probably necessary for viral infection.

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.

Identification of Merkel cell polyomavirus in the lower respiratory tract of Italian patients.

Merkel cell polyomavirus (MCPyV) has been found to be integrated monoclonally in a rare skin cancer named Merkel cell carcinoma. More recently, MCPyV has been detected in the upper respiratory tract of pediatric and adult patients. However, the mode of transmission and pathogenic role of MCPyV in the respiratory system has not been determined. In this study, MCPyV was sought in the lower respiratory tract of adult patients admitted to the hospital. MCPyV DNA was detected in 15 (17.24%) out of 87 lower respiratory tract samples. Most of the patients with MCPyV were over 50 years old. Nucleotide sequence of the t-antigen of MCPyV identified in respiratory secretions showed a homology to those found in Merkel cell carcinoma. In addition, phylogenetic analysis undertaken on the t-antigen sequences of Italian isolates and other MCPyVs identified in healthy and cancer tissues showed that all these isolates belonged to the same clade. Selective pressure analysis for the t-antigen revealed the presence of five sites under positive selection (omega = 4.3), with a posterior probabilities above 0.99. The alpha parameter of the gamma distribution was 0.01, showing that this distribution has a characteristic L-shape and suggesting a strong nucleotide substitution rate heterogeneity across sites. This study shows that MCPyV can infect the lower respiratory tract, but further investigations are needed to define its pathogenic role in respiratory diseases. J. Med. Virol. 82:505-509, 2010. (c) 2010 Wiley-Liss, Inc.