Complete-genome phylogenetic approach to varicella-zoster virus evolution: genetic divergence and evidence for recombination

Recent studies of varicella-zoster virus (VZV) DNA sequence variation, involving large numbers of globally distributed clinical isolates, suggest that this virus has diverged into at least three distinct genotypes designated European (E), Japanese (J), and mosaic (M). In the present study, we determined and analyzed the complete genomic sequences of two M VZV strains and compared them to the sequences of three E strains and two J strains retrieved from GenBank (including the Oka vaccine preparation, V-Oka). Except for a few polymorphic tandem repeat regions, the whole genome, representing approximately 125,000 nucleotides, is highly conserved, presenting a genetic similarity between the E and J genotypes of approximately 99.85%. These analyses revealed that VZV strains distinctly segregate into at least four genotypes (E, J, M1, and M2) in phylogenetic trees supported by high bootstrap values. Separate analyses of informative sites revealed that the tree topology was dependent on the region of the VZV genome used to determine the phylogeny; collectively, these results indicate the observed strain variation is likely to have resulted, at least in part, from interstrain recombination. Recombination analyses suggest that strains belonging to the M1 and M2 genotypes are mosaic recombinant strains that originated from ancestral isolates belonging to the E and J genotypes through recombination on multiple occasions. Furthermore, evidence of more recent recombination events between M1 and M2 strains is present in six segments of the VZV genome. As such, interstrain recombination in dually infected cells seems to figure prominently in the evolutionary history of VZV, a feature it has in common with other herpesviruses. In addition, we report here six novel genomic targets located in open reading frames 51 to 58 suitable for genotyping of clinical VZV isolates.     

Frequent recovery and broad genotype 2 diversity characterize hepatitis C virus infection in Ghana,West Africa

Hepatitis C virus (HCV) infection is thought to mostly become chronic and rarely resolves. HCV infection was serologically screened in 4,984 samples from Ghanaian blood donors, and 1.3% prevalence was found. At least 53% of confirmed anti-HCV carriers had no detectable viral RNA and were considered to have cleared the virus and recovered from the infection. Confirmation was authenticated by the presence of antibodies specific to at least two viral antigens, mostly NS3 and E2. Reactivity to HCV core antigens was lower in Ghanaian than United Kingdom blood donors. The minority of chronically infected donors carried a viral load significantly lower than an unselected comparative group of United Kingdom blood donors (2.5 x 10(5) versus 2.9 x 10(6) IU/ml; P = 0.004). HCV genotype 2 was largely predominant (87%). Sequence clustering was similarly broad in the E1/E2 and NS5 regions. The phylogenetic diversity and the incapacity to distinguish subtypes within genotype 2 in our and others' West African strains suggested that West Africa may be the origin of HCV genotype 2. The genetic diversity extended to the identification of strains clearly separated from known subtypes of genotype 2 and genotype 1. One strain appears to be part of a new HCV genotype. HCV infection in Ghana is characterized by a high rate of recovery and the predominance of broadly divergent genotype 2 strains.     

2010年山东省环境污水中脊髓灰质炎病毒的监测及其基因特征分析

本研究在山东省开展了脊髓灰质炎病毒(Poliovirus,PV)的外环境监测,从济南、临沂两地采集污水标本,浓缩处理后进行病毒分离,对分离到的PV采用中和试验进行血清定型,并对其VP1及3D区进行序列测定,分析其基因突变和重组情况。2010年,共采集污水标本32份,PV阳性10份,阳性率31.3%;分离到18株PV(PV1型3株,PV2型9株,PV3型6株),均为疫苗相关株,VP1完整编码区核苷酸变异数在0~4个之间,在3株PV2型病毒和4株PV3型病毒的基因组中发现重组;对VP1区影响神经毒力的减毒位点分析发现,PV1型病毒中有1株在nt 2 749发生突变(A→G),PV2型病毒中有1株在nt2 908发生A→G突变,3株在nt2 909发生U→C突变,6株PV3型病毒全部在nt2 493发生C→U突变。环境污水中可以分离到PV,其基因重组率和主要减毒位点的回复突变率较高,未发现脊灰野毒株和疫苗衍生株脊灰病毒(Vaccine-derived poliovirus,VDPV)。     

BF integrase genes of HIV-1 circulating in São Paulo, Brazil, with a recurrent recombination region

Although some studies have shown diversity in HIV integrase (IN) genes, none has focused particularly on the gene evolving in epidemics in the context of recombination. The IN gene in 157 HIV-1 integrase inhibitor-na?ve patients from the S?o Paulo State, Brazil, were sequenced tallying 128 of subtype B (23 of which were found in non-B genomes), 17 of subtype F (8 of which were found in recombinant genomes), 11 integrases were BF recombinants, and 1 from subtype C. Crucially, we found that 4 BF recombinant viruses shared a recurrent recombination breakpoint region between positions 4900 and 4924 (relative to the HXB2) that includes 2 gRNA loops, where the RT may stutter. Since these recombinants had independent phylogenetic origin, we argue that these results suggest a possible recombination hotspot not observed so far in BF CRF in particular, or in any other HIV-1 CRF in general. Additionally, 40% of the drug-na?ve and 45% of the drug-treated patients had at least 1 raltegravir (RAL) or elvitegravir (EVG) resistance-associated amino acid change, but no major resistance mutations were found, in line with other studies. Importantly, V151I was the most common minor resistance mutation among B, F, and BF IN genes. Most codon sites of the IN genes had higher rates of synonymous substitutions (dS) indicative of a strong negative selection. Nevertheless, several codon sites mainly in the subtype B were found under positive selection. Consequently, we observed a higher genetic diversity in the B portions of the mosaics, possibly due to the more recent introduction of subtype F on top of an ongoing subtype B epidemics and a fast spread of subtype F alleles among the B population.     

Molecular epidemiology of the foot-and-mouth disease virus outbreak in the United Kingdom in 2001

The objective of this study was to quantify the extent to which the genetic diversity of foot-and-mouth disease virus (FMDV) arising over the course of infection of an individual animal becomes fixed, is transmitted to other animals, and thereby accumulates over the course of an outbreak. Complete consensus sequences of 23 genomes (each of 8,200 nucleotides) of FMDV were recovered directly from epithelium tissue acquired from 21 farms infected over a nearly 7-month period during the 2001 FMDV outbreak in the United Kingdom. An analysis of these consensus sequences revealed very few apparently ambiguous sites but clear evidence of 197 nucleotide substitutions at 191 different sites. We estimated the rate of nucleotide substitution to be 2.26 x 10(-5) per site per day (95% confidence interval [CI], 1.75 x 10(-5) to 2.80 x 10(-5)) and nucleotide substitutions to accrue in the consensus sequence at an average rate of 1.5 substitutions per farm infection. This is a sufficiently high rate showing that detailed histories of the transmission pathways can be reliably reconstructed. Coalescent methods indicated that the date at which FMDV first infected livestock in the United Kingdom was 7 February 2001 (95% CI, 20 January to 19 February 2001), which was identical to estimates obtained on the basis of purely clinical evidence. Nucleotide changes appeared to have occurred evenly across the genome, and within the open reading frame, the ratio of nonsynonymous-to-synonymous change was 0.09. The ability to recover particular transmission pathways of acutely acting RNA pathogens from genetic data will help resolve uncertainties about how virus is spread and could help in the control of future epidemics.     

Genetic heterogeneity of hypervariable region 1 of the hepatitis C virus (HCV) genome and sensitivity of HCV to alpha interferon therapy

Hepatitis C virus (HCV) populations persist in vivo as a mixture of heterogeneous viruses called quasispecies. The relationship between the genetic heterogeneity of these variants and their responses to antiviral treatment remains to be elucidated. We have studied 26 virus strains to determine the influence of hypervariable region 1 (HVR-1) of the HCV genome on the effectiveness of alpha interferon (IFN-alpha) therapy. Following PCR amplification, we cloned and sequenced HVR-1. Pretreatment serum samples from 13 individuals with chronic hepatitis C whose virus was subsequently eradicated by treatment were compared with samples from 13 nonresponders matched according to the major factors known to influence the response, i.e., sex, genotype, and pretreatment serum HCV RNA concentration. The degree of virus variation was assessed by analyzing 20 clones per sample and by calculating nucleotide sequence entropy (complexity) and genetic distances (diversity). Types of mutational changes were also determined by calculating nonsynonymous substitutions per nonsynonymous site (K(a)) and synonymous substitutions per synonymous site (K(s)). The paired-comparison analysis of the nucleotide sequence entropy and genetic distance showed no statistical differences between responders and nonresponders. By contrast, nonsynonymous substitutions were more frequent than synonymous substitutions (P 相似文献   

BF Integrase Genes of HIV-1 Circulating in S?o Paulo,Brazil, with a Recurrent Recombination Region

Although some studies have shown diversity in HIV integrase (IN) genes, none has focused particularly on the gene evolving in epidemics in the context of recombination. The IN gene in 157 HIV-1 integrase inhibitor-naïve patients from the São Paulo State, Brazil, were sequenced tallying 128 of subtype B (23 of which were found in non-B genomes), 17 of subtype F (8 of which were found in recombinant genomes), 11 integrases were BF recombinants, and 1 from subtype C. Crucially, we found that 4 BF recombinant viruses shared a recurrent recombination breakpoint region between positions 4900 and 4924 (relative to the HXB2) that includes 2 gRNA loops, where the RT may stutter. Since these recombinants had independent phylogenetic origin, we argue that these results suggest a possible recombination hotspot not observed so far in BF CRF in particular, or in any other HIV-1 CRF in general. Additionally, 40% of the drug-naïve and 45% of the drug-treated patients had at least 1 raltegravir (RAL) or elvitegravir (EVG) resistance-associated amino acid change, but no major resistance mutations were found, in line with other studies. Importantly, V151I was the most common minor resistance mutation among B, F, and BF IN genes. Most codon sites of the IN genes had higher rates of synonymous substitutions (dS) indicative of a strong negative selection. Nevertheless, several codon sites mainly in the subtype B were found under positive selection. Consequently, we observed a higher genetic diversity in the B portions of the mosaics, possibly due to the more recent introduction of subtype F on top of an ongoing subtype B epidemics and a fast spread of subtype F alleles among the B population.     

Evolutionary dynamics of local pandemic H1N1/2009 influenza virus lineages revealed by whole-genome analysis

Virus gene sequencing and phylogenetics can be used to study the epidemiological dynamics of rapidly evolving viruses. With complete genome data, it becomes possible to identify and trace individual transmission chains of viruses such as influenza virus during the course of an epidemic. Here we sequenced 153 pandemic influenza H1N1/09 virus genomes from United Kingdom isolates from the first (127 isolates) and second (26 isolates) waves of the 2009 pandemic and used their sequences, dates of isolation, and geographical locations to infer the genetic epidemiology of the epidemic in the United Kingdom. We demonstrate that the epidemic in the United Kingdom was composed of many cocirculating lineages, among which at least 13 were exclusively or predominantly United Kingdom clusters. The estimated divergence times of two of the clusters predate the detection of pandemic H1N1/09 virus in the United Kingdom, suggesting that the pandemic H1N1/09 virus was already circulating in the United Kingdom before the first clinical case. Crucially, three clusters contain isolates from the second wave of infections in the United Kingdom, two of which represent chains of transmission that appear to have persisted within the United Kingdom between the first and second waves. This demonstrates that whole-genome analysis can track in fine detail the behavior of individual influenza virus lineages during the course of a single epidemic or pandemic.     

Recombination in Tula Hantavirus Evolution: Analysis of Genetic Lineages from Slovakia

To examine the evolution of Tula hantavirus (TUL), carried by the European common vole (Microtus arvalis and M. rossiaemeridionalis), we have analyzed genetic variants from Slovakia, the country where the virus is endemic. Phylogenetic analysis (PHYLIP) based on either partial (nucleotides [nt] 441 to 898) or complete N-protein-encoding sequences divided Slovakian TUL variants into two main lineages: (i) strains from eastern Slovakia, which clustered with Russian strains, and (ii) strains from western Slovakia situated closer to those from the Czech Republic. We found genetic diversity of 19% between the two groups and 4% within the western Slovakian TUL strains. Phylogenetic analysis of the 3′ noncoding region (3′-NCR), however, placed the eastern Slovakian strains closer to those from western Slovakia and the Czech Republic, with a greater distance to the Russian strains, suggesting a recombinant nature of the S segment in the eastern Slovakian TUL lineage. A bootscan search of the S-segment sequences of TUL strains revealed at least two recombination points in the S sequences of eastern Slovakian TUL strains (nt 400 to 415 and around 1200) which agreed well with the pattern of amino acid substitutions in the N protein and deletions/insertions in the 3′-NCR of the S segment. These data suggest that homologous recombination events occurred in the evolution of hantaviruses.     

Lack of Genetic Differentiation between Two Geographically Diverse Samples of Candida albicans Isolated from Patients Infected with Human Immunodeficiency Virus

The patterns of genetic variation of samples of Candida albicans isolated from patients infected with human immunodeficiency virus in Durham, N.C., and Vitória, Brazil, were compared. Genotypes for 126 strains were obtained at 16 polymorphic restriction sites distributed on nine PCR fragments. The results indicated evidence of clonality both within and between these two geographically diverse samples. The samples are genetically very similar, with little evidence of genetic differentiation.     

Phylogenetic and Genome-Wide Deep-Sequencing Analyses of Canine Parvovirus Reveal Co-Infection with Field Variants and Emergence of a Recent Recombinant Strain

Canine parvovirus (CPV), a fast-evolving single-stranded DNA virus, comprises three antigenic variants (2a, 2b, and 2c) with different frequencies and genetic variability among countries. The contribution of co-infection and recombination to the genetic variability of CPV is far from being fully elucidated. Here we took advantage of a natural CPV population, recently formed by the convergence of divergent CPV-2c and CPV-2a strains, to study co-infection and recombination. Complete sequences of the viral coding region of CPV-2a and CPV-2c strains from 40 samples were generated and analyzed using phylogenetic tools. Two samples showed co-infection and were further analyzed by deep sequencing. The sequence profile of one of the samples revealed the presence of CPV-2c and CPV-2a strains that differed at 29 nucleotides. The other sample included a minor CPV-2a strain (13.3% of the viral population) and a major recombinant strain (86.7%). The recombinant strain arose from inter-genotypic recombination between CPV-2c and CPV-2a strains within the VP1/VP2 gene boundary. Our findings highlight the importance of deep-sequencing analysis to provide a better understanding of CPV molecular diversity.     

Superinfection prevents recombination of the alphaherpesvirus bovine herpesvirus 1

Homologous recombination between strains of the same alphaherpesvirus species occurs frequently both in vitro and in vivo. This process has been described between strains of herpes simplex virus type 1, herpes simplex virus type 2, pseudorabies virus, feline herpesvirus 1, varicella-zoster virus, and bovine herpesvirus 1 (BoHV-1). In vivo, the rise of recombinant viruses can be modulated by different factors, such as the dose of the inoculated viruses, the distance between inoculation sites, the time interval between inoculation of the first and the second virus, and the genes in which the mutations are located. The effect of the time interval between infections with two distinguishable BoHV-1 on recombination was studied in three ways: (i) recombination at the level of progeny viruses, (ii) interference induced by the first virus infection on β-galactosidase gene expression of a superinfecting virus, and (iii) recombination at the level of concatemeric DNA. A time interval of 2 to 8 h between two successive infections allows the establishment of a barrier, which reduces or prevents any successful superinfection needed to generate recombinant viruses. The dramatic effect of the time interval on the rise of recombinant viruses is particularly important for the risk assessment of recombination between glycoprotein E-negative marker vaccine and field strains that could threaten BoHV-1 control and eradication programs.     

Geo-spatial hotspots of hemorrhagic fever with renal syndrome and genetic characterization of Seoul variants in Beijing, China

Background

Hemorrhagic fever with renal syndrome (HFRS) is highly endemic in mainland China, and has extended from rural areas to cities recently. Beijing metropolis is a novel affected region, where the HFRS incidence seems to be diverse from place to place.

Methodology/Principal Findings

The spatial scan analysis based on geographical information system (GIS) identified three geo-spatial “hotspots” of HFRS in Beijing when the passive surveillance data from 2004 to 2006 were used. The Relative Risk (RR) of the three “hotspots” was 5.45, 3.57 and 3.30, respectively. The Phylogenetic analysis based on entire coding region sequence of S segment and partial L segment sequence of Seoul virus (SEOV) revealed that the SEOV strains circulating in Beijing could be classified into at least three lineages regardless of their host origins. Two potential recombination events that happened in lineage #1 were detected and supported by comparative phylogenetic analysis. The SEOV strains in different lineages and strains with distinct special amino acid substitutions for N protein were partially associated with different spatial clustered areas of HFRS.

Conclusion/Significance

Hotspots of HFRS were found in Beijing, a novel endemic region, where intervention should be enhanced. Our data suggested that the genetic variation and recombination of SEOV strains was related to the high risk areas of HFRS, which merited further investigation.     

Divergence and recombination of clinical herpes simplex virus type 2 isolates

Herpes simplex virus type 2 (HSV-2) infects the genital mucosa and is one of the most common sexually transmitted viruses. Here we sequenced a segment comprising 3.5% of the HSV-2 genome, including genes coding for glycoproteins G, I, and E, from 27 clinical isolates from Tanzania, 10 isolates from Norway, and 10 isolates from Sweden. The sequence variation was low compared to that described for clinical HSV-1 isolates, with an overall similarity of 99.6% between the two most distant HSV-2 isolates. Phylogenetic analysis revealed a divergence into at least two genogroups arbitrarily designated A and B, supported by high bootstrap values and evolutionarily separated at the root. Genogroup A contained isolates collected in Tanzania, and genogroup B contained isolates collected in Tanzania and Scandinavia, implying that the genetic variability of HSV-2 is higher in Tanzania than in Scandinavia. Recombination network analysis and bootscan analysis revealed a complex pattern of phylogenetically conflicting informative sites in the sequence alignments. These signals were present in synonymous and nonsynonymous sites in all three genes and were not accumulated in specific regions, observations arguing against positive selection. Since the PHI test applied solely to synonymous sites revealed a high statistical probability of recombination, we suggest as a novel finding that homologous recombination is, as reported earlier for HSV-1 and varicella-zoster virus, a prominent feature in the evolution of HSV-2.     

Characterization, Genetic Diversity, and Evolutionary Link of Cucumber mosaic virus Strain New Delhi from India

The genome of Cucumber mosaic virus New Delhi strain (CMV-ND) from India, obtained from tomato, was completely sequenced and compared with full genome sequences of 14 known CMV strains from subgroups I and II, for their genetic diversity. Sequence analysis suggests CMV-ND shares maximum sequence identity at the nucleotide level with a CMV strain from Taiwan. Among all 15 strains of CMV, the encoded protein 2b is least conserved, whereas the coat protein (CP) is most conserved. Sequence identity values and phylogram results indicate that CMV-ND belongs to subgroup I. Based on the recombination detection program result, it appears that CMV is prone to recombination, and different RNA components of CMV-ND have evolved differently. Recombinational analysis of all 15 CMV strains detected maximum recombination breakpoints in RNA2; CP showed the least recombination sites.     

Evidence of recombination in intrapatient populations of hepatitis C virus

Hepatitis C virus (HCV) is a major cause of liver disease worldwide and a potential cause of substantial morbidity and mortality in the future. HCV is characterized by a high level of genetic heterogeneity. Although homologous recombination has been demonstrated in many members of the family Flaviviridae, to which HCV belongs, there are only a few studies reporting recombination on natural populations of HCV, suggesting that these events are rare in vivo. Furthermore, these few studies have focused on recombination between different HCV genotypes/subtypes but there are no reports on the extent of intra-genotype or intra-subtype recombination between viral strains infecting the same patient. Given the important implications of recombination for RNA virus evolution, our aim in this study has been to assess the existence and eventually the frequency of intragenic recombination on HCV. For this, we retrospectively have analyzed two regions of the HCV genome (NS5A and E1-E2) in samples from two different groups: (i) patients infected only with HCV (either treated with interferon plus ribavirin or treatment naïve), and (ii) HCV-HIV co-infected patients (with and without treatment against HIV). The complete data set comprised 17712 sequences from 136 serum samples derived from 111 patients. Recombination analyses were performed using 6 different methods implemented in the program RDP3. Recombination events were considered when detected by at least 3 of the 6 methods used and were identified in 10.7% of the amplified samples, distributed throughout all the groups described and the two genomic regions studied. The resulting recombination events were further verified by detailed phylogenetic analyses. The complete experimental procedure was applied to an artificial mixture of relatively closely viral populations and the ensuing analyses failed to reveal artifactual recombination. From these results we conclude that recombination should be considered as a potentially relevant mechanism generating genetic variation in HCV and with important implications for the treatment of this infection.     

Evolution of the core genome of Pseudomonas syringae, a highly clonal, endemic plant pathogen

Pseudomonas syringae is a common foliar bacterium responsible for many important plant diseases. We studied the population structure and dynamics of the core genome of P. syringae via multilocus sequencing typing (MLST) of 60 strains, representing 21 pathovars and 2 nonpathogens, isolated from a variety of plant hosts. Seven housekeeping genes, dispersed around the P. syringae genome, were sequenced to obtain 400 to 500 nucleotides per gene. Forty unique sequence types were identified, with most strains falling into one of four major clades. Phylogenetic and maximum-likelihood analyses revealed a remarkable degree of congruence among the seven genes, indicating a common evolutionary history for the seven loci. MLST and population genetic analyses also found a very low level of recombination. Overall, mutation was found to be approximately four times more likely than recombination to change any single nucleotide. A skyline plot was used to study the demographic history of P. syringae. The species was found to have maintained a constant population size over time. Strains were also found to remain genetically homogeneous over many years, and when isolated from sites as widespread as the United States and Japan. An analysis of molecular variance found that host association explains only a small proportion of the total genetic variation in the sample. These analyses reveal that with respect to the core genome, P. syringae is a highly clonal and stable species that is endemic within plant populations, yet the genetic variation seen in these genes only weakly predicts host association.     

Molecular characteristics and evolutionary analysis of a very virulent infectious bursal disease virus

Infectious bursal disease virus(IBDV) poses a significant threat to the poultry industry. Viral protein 2(VP2), the major structural protein of IBDV, has been subjected to frequent mutations that have imparted tremendous genetic diversity to the virus. To determine how amino acid mutations may affect the virulence of IBDV, we built a structural model of VP2 of a very virulent strain of IBDV identified in China, vv IBDV Gx, and performed a molecular dynamics simulation of the interaction between virulence sites. The study showed that the amino acid substitutions that distinguish vv IBDV from attenuated IBDV(H253Q and T284A) favor a hydrophobic and flexible conformation of ?-barrel loops in VP2, which could promote interactions between the virus and potential IBDV-specific receptors. Population sequence analysis revealed that the IBDV strains prevalent in East Asia show a significant signal of positive selection at virulence sites 253 and 284. In addition, a signal of co-evolution between sites 253 and 284 was identified. These results suggest that changes in the virulence of IBDV may result from both the interaction and the co-evolution of multiple amino acid substitutions at virulence sites.     

High-frequency RNA recombination of murine coronaviruses.

The RNA genome of coronaviruses consists of a single species of nonsegmented RNA. In this communication, we demonstrate that the RNA genomes of different strains of murine coronaviruses recombine during mixed infection at a very high frequency. Susceptible cells were coinfected with a temperature-sensitive mutant of one strain of mouse hepatitis virus (MHV) and a wild-type virus of a different strain. Of 21 randomly isolated viruses released from the coinfected cells at the nonpermissive temperature, 2 were recombinants which differed in the site of recombination. After three serial passages of the original virus pool derived from the mixed infection, the majority of the progeny viruses were recombinants. These recombinant viruses represented at least five different recombination sites between the two parental MHV strains. Such a high-frequency recombination between nonsegmented RNA genomes of MHV suggests that segmented RNA intermediates might be generated during MHV replication. We propose that the RNA replication of MHV proceeds in a discontinuous and nonprocessive manner, thus generating free segmented RNA intermediates, which could be used in RNA recombination via a copy-choice mechanism.