Full Genome Sequence Analysis Research Articles

Genomic events contributing to the high prevalence of amantadine-resistant influenza A/H3N2

The prevalence of amantadine-resistant influenza A/H3N2 viruses (belonging to the N-lineage), possessing an S31N mutation in the M2 protein and S193F and D225N substitutions in their HA1 subunit, has significantly increased worldwide since 2005. The aim of this study was to clarify the genomic events contributing to the evolution and continuity of the N-lineage amantadine-resistant viruses. The full genome sequence of A/H3N2 isolates, including both amantadine-resistant and amantadine-sensitive viruses, collected in Japan between 2006 and 2008, was determined and phylogenetically compared with isolates obtained from the database. On the basis of the full genome sequence analysis, the N-lineage could be further divided into three genetically related clades: N1 (A/Wisconsin/67/2005-like amantadine-resistant viruses from years 2005-2007), N2 (amantadine-sensitive viruses from 2007) and N3 (A/Brisbane/10/2007-like amantadine-resistant viruses from 2007 and 2008). The 2006/2007 season showed cocirculation of antigenic variants of amantadine-resistant viruses of clades N1 and N3 in addition to the N2-sensitive viruses. In the 2007/2008 season, the clade N3 amantadine-resistant lineage dominated and replaced other strains. Phylogenetic analysis of each individual segment suggested that N2 and N3 were generated from two independent reassortment events involving clade N1 viruses and pre-N-lineage strains. Our data show that several reassortment events have contributed to the evolution of amantadine-resistant A/H3N2 strains and, consequently, to the successful spread of this lineage. Although amantadine resistance is caused by single amino acid mutations in the M2 protein, genome-wide adjustment involving multiple genes appears to be necessary to obtain efficient replication and transmission of resistant viruses. Such adjustments are attainable through reassortment of segments among different virus lineages.

Identification de virus de la fièvre catarrhale ovine réassortis dans la région Ouest méditerranéenne

Bluetongue virus (BTV) is the species of the genus Orbivirus of the family Reoviridae. BTV can infect most ruminants, caus­ing a severe haemorrhagic disease called bluetongue in sheep. BTV is transmitted among ruminant hosts by certain species of Culicoides (biting midges). The BTV genome is composed of ten linear double-stranded (ds) RNA genome segments, encod­ing seven structural and three non-structural proteins. The seg­mented nature of the genome allows different BTV strains infect­ing the same cell to exchange (reassort) genome segments. Nine BTV serotypes have been detected in Europe since 1998, including strains belonging to both eastern (BTV-1, 9, 16) and western lineages (BTV-1, 2, 4, 6, 8, 11). Live attenuated mono­valent vaccine strains of BTV-2, 4, 8, 9 (western group) and BTV- 16 (eastern group) have also been used in the Mediterranean region, in attempts to minimise virus circulation. The release of these vaccine strains, some of which have persisted in the field (including BTV-2 and 16), has added further genetic diversity, generating an unprecedented mix of field and vaccine strain viruses. These events have provided unique opportunities for genome segment exchange (reassortment) between different BTV strains and topotypes. Indeed, a strain of BTV-2, derived by reassortment between vaccine strains of BTV-2 and 16 has been previously detected in Italy during 2002. Full genome sequence analysis of BTV-2 and BTV-4 isolates (1999-2004) from the Western Mediterranean region helped to identify multiple reassortant viruses, involving the exchange of several different genome segments. Reassortant viruses were identified as containing genes derived from different western field strains, from western field and vaccine strains, and from eastern field and western vaccine strains. The detection of these reassortant BTVs in Europe highlights concerns about the use of live BTV vaccines in the region.

Evolution of Dengue Virus Type 2 during Two Consecutive Outbreaks with an Increase in Severity in Southern Taiwan in 2001–2002

To investigate viral determinants and evolution linked to outbreak with increased severity, we examined dengue virus type 2 (DENV-2) sequences from plasma of 31 patients (14 dengue fever, 17 dengue hemorrhagic fever, DHF) continuously during the 2001 and 2002 outbreaks in southern Taiwan, in which both the total cases and proportion of DHF cases increased. Analysis of envelope (E) and full-genome sequences between viruses of the two outbreaks revealed 5 nucleotide changes in E, NS1, NS4A, and NS5 genes. None was identical to those reported in the DENV-2 outbreak in Cuba in 1997, suggesting viral determinants linked to severe outbreak are genotype dependent. Compared with previous reports of lineage turnover years apart, our findings that the 2002 viruses descended from a minor variant of the 2001 viruses in less than 6 months was novel, and may represent a mechanism of evolution of DENV from one outbreak to another.

Identification of New CRF43_02G and CRF25_cpx in Saudi Arabia Based on Full Genome Sequence Analysis of Six HIV Type 1 Isolates

Recently, we reported a high level of HIV-1 strain diversity in patients at the King Faisal Specialist Hospital and Research Center, Jeddah, Saudi Arabia. Based on phylogenetic analysis of gag p24, pol integrase, and env gp41 sequences, subtypes A, B, C, D, and G, and CRF02_AG, as well as unique recombinant forms were identified. Subtype G accounted for 25% of the infections in the Saudi population and this high prevalence was unexpected. Although subtype G is found in west central Africa, pure subtype G strains are uncommon. To further characterize the subtype G infections in Saudi Arabia, six strains that appeared to be pure subtype G were selected for full genome sequencing. Near full-length genomes were obtained using RT-PCR amplification to generate overlapping fragments from viral RNA extracted from plasma. The six strains are not subtype G throughout their entire genome. Four isolates have a recombinant structure composed of CRF02_AG and subtype G and share three identical breakpoints. This recombinant form defines a new CRF designated CRF43_02G. The remaining two isolates are CRF25_cpx, a circulating recombinant form identified in Cameroon composed of subtypes A and G and unclassified segments. Reanalysis of the previously reported Saudi HIV-1 partial genome sequences revealed additional isolates classified as CRF43_02G and CRF25_cpx and one isolate was reclassified to CRF22_01A. Identification of CRF43_02G in Saudi Arabia could indicate a transmission network within the country. Alternatively, the new CRF could have been introduced from an external source where this CRF is not yet recognized.

The VG/GA strain of Newcastle disease virus: mucosal immunity, protection against lethal challenge and molecular analysis

The Villegas–Glisson/University of Georgia (VG/GA) strain of Newcastle disease virus (NDV) isolated from the intestine of healthy turkeys has been proposed to replicate in the respiratory and intestinal tract of chickens. In the present study, the virus distribution, the mucosal and systemic immune response, the efficacy against lethal challenge and the full genome sequence of the VG/GA strain were compared with the La Sota strain of NDV. The VG/GA strain was detected at different time points in the respiratory and intestinal tract of chickens with a preferential tropism for the latter. Both the VG/GA and La Sota strains induced NDV-specific immunoglobulin A (IgA) at the upper respiratory tract. IgA levels were higher in the trachea for the La Sota strain, while they were higher in the bile and intestine for the VG/GA strain. Positive correlation between virus distribution of the viruses and IgA production was observed. Despite the presence of the maternal antibodies in broilers, early vaccination with the VG/GA strain afforded 95% to 100% protection against lethal challenge, equivalent to the protection conferred by the La Sota strain. Full genome sequence analysis classified the VG/GA strain within class II, genotype II viruses, which also include most of the respirotropic vaccine strains. Differences with the La Sota strain at the nucleotide and amino acid levels that may explain the differential phenotype of the VG/GA were observed; however, verification of the significance of those changes is required. Taken together, these results validate field observations on the efficacy of VG/GA vaccination and demonstrated the unique characteristics of the strain.

Comparative full genome analysis revealed E1: A226V shift in 2007 Indian Chikungunya virus isolates

The resurgence of Chikungunya virus (CHIKV) in the form of unprecedented explosive epidemic after a gap of 32 years in India is a point of major public health concern. In 2007 again there was outbreak in Kerala, India, affecting more than 25,000 cases with many reported mortalities. To understand the molecular basis of this high virulence and its implication in large-scale epidemic, a detailed systematic serological, virological and molecular investigation was undertaken with the epidemic samples of Kerala-2007. The comparative analysis of full genome sequence of Chikungunya virus isolate of 2007 with 2006 revealed three unique substitutions in structural and non-structural genes of 2007 isolate [two in E1 region (V14A and A226V) and one in Nsp1 (M184T)]. Our finding further substantiates the association of A226V shift in E1 gene with evolutionary success possibly due to adaptation in the mosquito vector with progression of epidemic, as observed in Reunion Island. This A226V shift which was absent in all 2006 Indian isolates, is found to be present in the four 2007 isolates, analysed in this study. These unique molecular features of the 2007 isolates with the progression of the epidemic from 2005 to 2007 demonstrate their high evolutionary and epidemic potential and thereby suggesting possible implication in higher magnitude and virulence of this outbreak.

Novel Approach for Genotyping Varicella-Zoster Virus Strains from Germany

In this study, we present a novel genotyping scheme to classify German wild-type varicella-zoster virus (VZV) strains and to differentiate them from the Oka vaccine strain (genotype B). This approach is based on analysis of four loci in open reading frames (ORFs) 51 to 58, encompassing a total length of 1,990 bp. The new genotyping scheme produced identical clusters in phylogenetic analyses compared to full-genome sequences from well-characterized VZV strains. Based on genotype A, D, B, and C reference strains, a dichotomous identification key (DIK) was developed and applied for VZV strains obtained from vesicle fluid and liquor samples originating from 42 patients suffering from varicella or zoster between 2003 and 2006. Sequencing of regions in ORFs 51, 52, 53, 56, 57, and 58 identified 18 single-nucleotide polymorphisms (SNPs), including two novel ones, SNP 89727 and SNP 92792 in ORF51 and ORF52, respectively. The DIK as well as phylogenetic analysis by Bayesian inference showed that 14 VZV strains belonged to genotype A, and 28 VZV strains were classified as genotype D. Neither Japanese (vaccine)-like B strains nor recombinant-like C strains were found within the samples from Germany. The novel genotyping scheme and the DIK were demonstrated to be practical and simple and allow the highly efficient replication of phylogenetic patterns in VZV initially derived from full-genome DNA sequence analyses. Therefore, this approach may allow us to draw a more comprehensive picture of wild-type VZV strains circulating in Germany and Central Europe by high-throughput procedures in the future.

Resistance to Alpha/Beta Interferon Is a Determinant of West Nile Virus Replication Fitness and Virulence

The emergence of West Nile virus (WNV) in the Western Hemisphere is marked by the spread of pathogenic lineage I strains, which differ from typically avirulent lineage II strains. To begin to understand the virus-host interactions that may influence the phenotypic properties of divergent lineage I and II viruses, we compared the genetic, pathogenic, and alpha/beta interferon (IFN-alpha/beta)-regulatory properties of a lineage II isolate from Madagascar (MAD78) with those of a new lineage I isolate from Texas (TX02). Full genome sequence analysis revealed that MAD78 clustered, albeit distantly, with other lineage II strains, while TX02 clustered with emergent North American isolates, more specifically with other Texas strains. Compared to TX02, MAD78 replicated at low levels in cultured human cells, was highly sensitive to the antiviral actions of IFN in vitro, and demonstrated a completely avirulent phenotype in wild-type mice. In contrast to TX02 and other pathogenic forms of WNV, MAD78 was defective in its ability to disrupt IFN-induced JAK-STAT signaling, including the activation of Tyk2 and downstream phosphorylation and nuclear translocation of STAT1 and STAT2. However, replication of MAD78 was rescued in cells with a nonfunctional IFN-alpha/beta receptor (IFNAR). Consistent with this finding, the virulence of MAD78 was unmasked upon infection of mice lacking IFNAR. Thus, control of the innate host response and IFN actions is a key feature of WNV pathogenesis and replication fitness.

An Attenuated LC16m8 Smallpox Vaccine: Analysis of Full-Genome Sequence and Induction of Immune Protection

The potential threat of smallpox bioterrorism has made urgent the development of lower-virulence vaccinia virus vaccines. An attenuated LC16m8 (m8) vaccine was developed in 1975 from the Lister strain used in the World Health Organization smallpox eradication program but was not used against endemic smallpox. Today, no vaccines can be tested with variola virus for efficacy in humans, and the mechanisms of immune protection against the major intracellular mature virion (IMV) and minor extracellular enveloped virion (EEV) populations of poxviruses are poorly understood. Here, we determined the full-genome sequences of the m8, parental LC16mO (mO), and grandparental Lister (LO) strains and analyzed their evolutionary relationships. Sequence data and PCR analysis indicated that m8 was a progeny of LO and that m8 preserved almost all of the open reading frames of vaccinia virus except for the disrupted EEV envelope gene B5R. In accordance with this genomic background, m8 induced 100% protection against a highly pathogenic vaccinia WR virus in mice by a single vaccination, despite the lack of anti-B5R and anti-EEV antibodies. The immunogenicity and priming efficacy with the m8 vaccine consisting mainly of IMV were as high as those with the intact-EEV parental mO and grandparental LO vaccines. Thus, mice vaccinated with 10(7) PFU of m8 produced low levels of anti-B5R antibodies after WR challenge, probably because of quick clearance of B5R-expressing WR EEV by strong immunity induced by the vaccination. These results suggest that priming with m8 IMV provides efficient protection despite undetectable levels of immunity against EEV.

Biochemical Clustering of Monomeric GTPases of the Ras Superfamily

To date phylogeny has been used to compare entire families of proteins based on their nucleotide or amino acid sequence. Here we developed a novel analytical platform allowing a systematic comparison of protein families based on their biochemical properties. This approach was validated on the Rho subfamily of GTPases. We used two high throughput methods, referred to as AlphaScreen and FlashPlate, to measure nucleotide binding capacity, exchange, and hydrolysis activities of small monomeric GTPases. These two technologies have the characteristics to be very sensitive and to allow homogenous and high throughput assays. To analyze and integrate the data obtained, we developed an algorithm that allows the classification of GTPases according to their enzymatic activities. Integration and hierarchical clustering of these results revealed unexpected features of the small Rho GTPases when compared with primary sequence-based trees. Hence we propose a novel phylobiochemical classification of the Ras superfamily of GTPases.

Nucleotide mutations associated with hepatitis B e antigen negativity

One hundred and forty four patients with chronic hepatitis B were tested to identify new mutations associated with hepatitis B e antigen (HBeAg) negativity, using a full genome sequence analysis. All the patients were Chinese and had hepatitis B virus infection of genotype C. Patients with none of the pre-core or core promoter mutations were significantly (P < 0.001) less common in the group with anti-HBe (13%) than in the group with HBeAg (56%). The complete nucleotide sequence was determined in four anti-HBe-positive patients who had neither pre-core nor core promoter mutations and in five HBeAg-positive patients who also had neither of these mutations (the groups were matched for age and sex). Six mutations were found to be significantly more common in the former group than in the latter: G529A (3/4 vs. 0/5), C934A (4/4 vs. 1/5), A1053G (4/4 vs. 1/5), G1915T/A (4/4 vs. 0/5), T2005C/A (4/4 vs. 0/5), and C3026T (3/4 vs. 0/5). Three of the six mutations were significantly more common in the four anti-HBe-positive patients who had neither pre-core nor core promoter mutations, compared to 11 HBeAg-positive patients who had pre-core and core promoter mutations, and also compared to 15 anti-HBe-positive patients who had pre-core and core promoter mutations, suggesting further the specificity of these mutations. Of the six mutations, two resulted in amino acid substitution in the polymerase protein, and one is located near the enhancer I region. The results suggest that the six newly discovered mutations are associated with HBeAg negativity.

Molecular basis of adaptive convergence in experimental populations of RNA viruses.

Characterizing the molecular basis of adaptation is one of the most important goals in modern evolutionary genetics. Here, we report a full-genome sequence analysis of 21 independent populations of vesicular stomatitis ribovirus evolved on the same cell type but under different demographic regimes. Each demographic regime differed in the effective viral population size. Evolutionary convergences are widespread both at synonymous and nonsynonymous replacements as well as in an intergenic region. We also found evidence for epistasis among sites of the same and different loci. We explain convergences as the consequence of four factors: (1) environmental homogeneity that supposes an identical challenge for each population, (2) structural constraints within the genome, (3) epistatic interactions among sites that create the observed pattern of covariation, and (4) the phenomenon of clonal interference among competing genotypes carrying different beneficial mutations. Using these convergences, we have been able to estimate the fitness contribution of the identified mutations and epistatic groups. Keeping in mind statistical uncertainties, these estimates suggest that along with several beneficial mutations of major effect, many other mutations got fixed as part of a group of epistatic mutations.

Hepatitis C Virus Evolutionary Patterns Studied Through Analysis of Full-Genome Sequences

The evolutionary patterns of hepatitis C virus (HCV), including the best-fitting nucleotide substitution model and the molecular clock hypothesis, were investigated by analyzing full-genome sequences available in the HCV database. The likelihood ratio test allowed us to discriminate among different evolutionary hypotheses. The phylogeny of the six major HCV types was accurately inferred, and the final tree was rooted by reconstructing the hypothetical HCV common ancestor with the maximum likelihood method. The presence of phylogenetic noise and the relative nucleotide substitution rates in the different HCV genes were also examined. These results offer a general guideline for the future of HCV phylogenetic analysis and also provide important insights on HCV origin and evolution.

Tempo and mode of human and simian T-lymphotropic virus (HTLV/STLV) evolution revealed by analyses of full-genome sequences.

We investigated the tempo and mode of evolution of the primate T-lymphotropic viruses (PTLVs). Several different models of nucleotide substitution were tested on a general phylogenetic tree obtained using the 20 full-genome HTLV/STLV sequences available. The likelihood ratio test showed that the Tamura and Nei model with discrete gamma-distributed rates among sites is the best-fitting substitution model. The heterogeneity of nucleotide substitution rates along the PTLV genome was further investigated for different genes and at different codon positions (cdp's). Tests of rate constancy showed that different PTLV lineages evolve at different rates when first and second cdp's are considered, but the molecular-clock hypothesis holds for some PTLV lineages when the third cdp is used. Negative selection was evident throughout the genome. However, in the gp46 region, a small fragment subjected to positive selection was identified using a Monte Carlo simulation based on a likelihood method. Employing correlations of the virus divergence times with anthropologically documented migrations of their host, a possible timescale was estimated for each important node of the PTLV tree. The obtained results on these slow-evolving viruses could be used to fill gaps in the historical records of some of the host species. In particular, the HTLV-I/STLV-I history might suggest a simian migration from Asia to Africa not much earlier than 19,500-60,000 years ago.

Full-Genome Sequence Analyses of Hepatitis B Virus (HBV) Strains Recovered from Chimpanzees Infected in the Wild: Implications for an Origin of HBV

Hepatitis B virus (HBV) belongs to the genus Orthohepadnavirus of the family Hepadnaviridae. Having been found in various animals (duck, heron, woodchuck, ground squirrel, and primates), hepadnaviruses must have undergone a long history of evolution and may comprise more members than currently recognized. Chimpanzees may also have their own hepadnavirus, even if it might be very close to HBV. We analyzed HBV-like sequences from three chimpanzees (Pan troglodytes) that were most likely infected during their life in Africa in the wild. Two chimpanzees (Ch256 and Ch258) possessed a viral genome of 3182 nt in length with a 33-nt deletion in the preS1 region, which could not be classified into any of the six genotypes (A–F) of human HBV but was very homologous to a previously reported isolate from a London Zoo chimpanzee. Phylogenetically distinct from the HBV-like sequences from gibbons, orangutans, and a gorilla so far reported, the Ch256 and Ch258 isolates would represent an indigenous chimpanzee HBV (tentatively ChHBV). A third chimpanzee (Ch195) had a 3212-nt genome, classifiable into the genotype E of HBV. Because HBV-E has been found mostly in Africans, Ch195 may have been infected from a human source in Africa. However, an inverse scenario is also possible: a spread of HBV-E might have occurred from chimpanzees to humans a long time ago in Africa. Analysis of the arginine-rich C-terminal region of the core protein, which is well conserved among mammalian hepadnaviruses, indicated that HBV-E/F and nonhuman primate hepadnaviruses are much closer than HBV-A/B/C/D to the hepadnaviruses of woodchuck and ground squirrel. Our results support an “ex-nonhuman primate” hypothesis for the origin of HBV.

Full-genome sequence analysis of clinically representative hepatitis C virus isolates

Full-genome sequence analysis of clinically representative hepatitis C virus isolates