Clinical features and laboratory findings of human parvovirus B19 in human immunodeficiency virus-infected patients

Immunocompromised patients may develop severe chronic anaemia when infected by human parvovirus B19 (B19V). However, this is not the case in human immunodeficiency virus (HIV)-infected patients with good adherence to highly active antiretroviral treatment (HAART). In this study, we investigated the clinical evolution of five HIV-infected patients receiving HAART who had B19V infections confirmed by serum polymerase chain reaction. Four of the patients were infected with genotype 1a strains and the remaining patient was infected with a genotype 3b strain. Anaemia was detected in three of the patients, but all patients recovered without requiring immunoglobulin and/or blood transfusions. In all cases, the attending physicians did not suspect the B19V infections. There was no apparent relationship between the infecting genotype and the clinical course. In the HAART era, B19V infections in HIV-positive patients may be limited, subtle or unapparent.     

Characterization of Parvovirus B19 genotype 2 in KU812Ep6 cells

An infectious parvovirus B19 (B19V) genotype 2 variant was identified as a high-titer contaminant in a human plasma donation. Genome analysis revealed a 138-bp insertion within the p6 promoter. The inserted sequence was represented by an additional 30 bp from the end of the inverted terminal repeat adjacent to a 108-bp element found also, in inverted orientation, at the extreme right end of the unique sequence of the genome. However, despite the profound variations in the promoter region, the pattern of gene expression and DNA replication did not differ between genotype 1 and genotype 2 in permissive erythroid KU812Ep6 cells. Capsid proteins of both genotypes differ in their amino acid sequences. However, equivalent kinetics of virus inactivation at 56 degrees C or pH 4 indicated a comparable physicochemical stability of virus capsids. Sera from six individuals infected by B19V genotype 1 were investigated on cross-neutralization of B19V genotype 2 in vitro. Similar neutralization of both B19V genotypes was observed in sera from three individuals, while the sera from three other individuals showed weaker cross-neutralization for genotype 2. In conclusion, the in vitro replication characteristics and physical stability of B19V capsids are very similar between human parvovirus B19 genotypes 1 and 2, and cross-neutralization indicates a close antigenic relation of genotypes 1 and 2.     

Tracking of peptide-specific CD4+ T-cell responses after an acute resolving viral infection: a study of parvovirus B19

The evolution of peptide-specific CD4(+) T-cell responses to acute viral infections of humans is poorly understood. We analyzed the response to parvovirus B19 (B19), a ubiquitous and clinically significant pathogen with a compact and conserved genome. The magnitude and breadth of the CD4(+) T-cell response to the two B19 capsid proteins were investigated using a set of overlapping peptides and gamma interferon-specific enzyme-linked immunospot assays of peripheral blood mononuclear cells (PBMCs) from a cohort of acutely infected individuals who presented with acute arthropathy. These were compared to those for a cohort of B19-specific immunoglobulin M-negative (IgM(-)), IgG(+) remotely infected individuals. Both cohorts of individuals were found to make broad CD4(+) responses. However, while the responses following acute infection were detectable ex vivo, responses in remotely infected individuals were only detected after culture. One epitope (LASEESAFYVLEHSSFQLLG) was consistently targeted by both acutely (10/12) and remotely (6/7) infected individuals. This epitope was DRB1*1501 restricted, and a major histocompatibility complex peptide tetramer stained PBMCs from acutely infected individuals in the range of 0.003 to 0.042% of CD4(+) T cells. Tetramer-positive populations were initially CD62L(lo); unlike the case for B19-specific CD8(+) T-cell responses, however, CD62L was reexpressed at later times, as responses remained stable or declined slowly. This first identification of B19 CD4(+) T-cell epitopes, including a key immunodominant peptide, provides the tools to investigate the breadth, frequency, and functions of cellular responses to this virus in a range of specific clinical settings and gives an important reference point for analysis of peptide-specific CD4(+) T cells during acute and persistent virus infections of humans.     

Parvovirus B19 1A complete genome from a fatal case in Brazil

Parvovirus B19 (B19V) infects individuals worldwide and is associated with an ample range of pathologies and clinical manifestations. B19V is classified into three distinct genotypes, all identified in Brazil. Here, we report a complete sequence of a B19V genotype 1A that was obtained by high-throughput metagenomic sequencing. This genome provides information that will contribute to the studies on B19V epidemiology and evolution.     

Correlation of RNA binding affinity of avian oncornavirus p19 proteins with the extent of processing of virus genome RNA in cells.

We purified the p19 proteins from the Prague C strain of Rous sarcoma virus, avian myeloblastosis virus, B77 sarcoma virus, myeloblastosis-associated virus-2(0), and PR-E 95-C virus and measured their binding affinities for 60S viral RNA by the nitrocellulose filter binding technique. The apparent association constants of the p19 proteins from Rous sarcoma virus Prague C, avian myeloblastosis virus, and B77 sarcoma virus for homologous and heterologous 60S RNAs were similar (1.5 x 10(11) to 2.6 x 10(11) liters/mol), whereas those of myeloblastosis-associated virus-2(0) and PR-E 95-C virus were 10-fold lower. The sizes and relative amounts of the virus-specific polyadenylic acid-containing RNAs in the cytoplasms of cells infected with Rous sarcoma virus Prague C, myeloblastosis-associated virus-2(0), and PR-E 95-C virus were determined by fractionating the RNAs on agarose gels containing methylmercury hydroxide, transferring them to diazobenzyloxymethyl paper and hybridizing them to a 70-nucleotide complementary DNA probe. In cells infected with Rous sarcoma virus Prague C we detected 3.4 x 10(6)-, 1.9 x 10(6)-, and 1.1 x 10(6)-dalton RNAs, in PR-E 95-C virus-infected cells we detected 3.4 x 10(6)-, 1.9 x 10(6)- and 0.7 x 10(6)-dalton RNAs, and in cells infected with myeloblastosis-associated virus-2(0) we detected 3 x 10(6)- and 1.3 x 10(6)-dalton RNAs. Each of these RNA species contained RNA sequences derived from the 5' terminus of genome-length RNA, as evidenced by hybridization with the 5' 70-nucleotide complementary DNA. The ratios of subgenomic mRNA's to genome-length RNAs in cells infected with myeloblastosis-associated virus-2(0) and PR-E 95-C virus were three- to five-fold higher than the ratio in cells infected with Rous sarcoma virus Prague C. These results suggest that more processing of viral RNA in infected cells is correlated with lower binding affinities of the p19 protein for viral RNA, and they are consistent with the hypothesis that the p19 protein controls processing of viral RNA in cells.     

Recombination Enhances HIV-1 Envelope Diversity by Facilitating the Survival of Latent Genomic Fragments in the Plasma Virus Population

HIV-1 is subject to immune pressure exerted by the host, giving variants that escape the immune response an advantage. Virus released from activated latent cells competes against variants that have continually evolved and adapted to host immune pressure. Nevertheless, there is increasing evidence that virus displaying a signal of latency survives in patient plasma despite having reduced fitness due to long-term immune memory. We investigated the survival of virus with latent envelope genomic fragments by simulating within-host HIV-1 sequence evolution and the cycling of viral lineages in and out of the latent reservoir. Our model incorporates a detailed mutation process including nucleotide substitution, recombination, latent reservoir dynamics, diversifying selection pressure driven by the immune response, and purifying selection pressure asserted by deleterious mutations. We evaluated the ability of our model to capture sequence evolution in vivo by comparing our simulated sequences to HIV-1 envelope sequence data from 16 HIV-infected untreated patients. Empirical sequence divergence and diversity measures were qualitatively and quantitatively similar to those of our simulated HIV-1 populations, suggesting that our model invokes realistic trends of HIV-1 genetic evolution. Moreover, reconstructed phylogenies of simulated and patient HIV-1 populations showed similar topological structures. Our simulation results suggest that recombination is a key mechanism facilitating the persistence of virus with latent envelope genomic fragments in the productively infected cell population. Recombination increased the survival probability of latent virus forms approximately 13-fold. Prevalence of virus with latent fragments in productively infected cells was observed in only 2% of simulations when we ignored recombination, while the proportion increased to 27% of simulations when we allowed recombination. We also found that the selection pressures exerted by different fitness landscapes influenced the shape of phylogenies, diversity trends, and survival of virus with latent genomic fragments. Our model predicts that the persistence of latent genomic fragments from multiple different ancestral origins increases sequence diversity in plasma for reasonable fitness landscapes.     

Molecular Evolutionary and Epidemiological Dynamics of Genotypes 1G and 2B of Rubella Virus

Rubella Virus (RV), which causes measles-like rashes in children, puts millions of infants at risk of congenital defects across the globe. Employing phylogenetic approaches to the whole genome sequence data and E1 glycoprotein sequence data, the present study reports the substitution rates and dates of emergence of all thirteen previously described rubella genotypes, and gains important insights into the epidemiological dynamics of two geographically widely distributed genotypes 1G and 2B. The overall nucleotide substitution rate of this non-vector-borne RV is in the order of 10⁻³ substitutions/site/year, which is considerably higher than the substitution rates previously reported for the vector-borne alphaviruses within the same family. Currently circulating strains of RV share a common ancestor that existed within the last 150 years, with 95% Highest Posterior Density values ranging from 1868 to 1926 AD. Viral strains within the respective genotypes began diverging between the year 1930 s and 1980 s. Both genotype 1G and 2B have shown a decline in effective number of infections since 1990 s, a period during which mass immunization programs against RV were adapted across the globe. Although both genotypes showed some extent of spatial genetic structuring, the analyses also depicted an inter-continental viral dispersal. Such a viral dispersal pattern could be related to the migration of infected individuals across the regions coupled with a low coverage of MMR vaccination.     

Recombination in the genesis and evolution of hepatitis B virus genotypes

Hepatitis B virus (HBV) infection is widely distributed in both human and ape populations throughout the world and is a major cause of human morbidity and mortality. HBV variants are currently classified into the human genotypes A to H and species-associated chimpanzee and gibbon/orangutan groups. To examine the role of recombination in the evolution of HBV, large-scale data retrieval and automated phylogenetic analysis (TreeOrder scanning) were carried out on all available published complete genome sequences of HBV. We detected a total of 24 phylogenetically independent potential recombinants (different genotype combinations or distinct breakpoints), eight of which were previously undescribed. Instances of intergenotype recombination were observed in all human and ape HBV variants, including evidence for a novel gibbon/genotype C recombinant among HBV variants from Vietnam. By recording sequence positions in trees generated from sequential fragments across the genome, violations of phylogeny between trees also provided evidence for frequent intragenotype recombination between members of genotypes A, D, F/H, and gibbon variants but not in B, C, or the Asian B/C recombinant group. In many cases, favored positions for both inter- and intragenotype recombination matched positions of phylogenetic reorganization between the human and ape genotypes, such as the end of the surface gene and the core gene, where sequence relationships between genotypes changed in the TreeOrder scan. These findings provide evidence for the occurrence of past, extensive recombination events in the evolutionary history of the currently classified genotypes of HBV and potentially in changes in its global epidemiology and associations with human disease.     

Substitution Model of Sequence Evolution for the Human Immunodeficiency Virus Type 1 Subtype B gp120 Gene over the C2-V5 Region

Phylogenetic analyses frequently rely on models of sequence evolution that detail nucleotide substitution rates, nucleotide frequencies, and site-to-site rate heterogeneity. These models can influence hypothesis testing and can affect the accuracy of phylogenetic inferences. Maximum likelihood methods of simultaneously constructing phylogenetic tree topologies and estimating model parameters are computationally intensive, and are not feasible for sample sizes of 25 or greater using personal computers. Techniques that initially construct a tree topology and then use this non-maximized topology to estimate ML substitution rates, however, can quickly arrive at a model of sequence evolution. The accuracy of this two-step estimation technique was tested using simulated data sets with known model parameters. The results showed that for a star-like topology, as is often seen in human immunodeficiency virus type 1 (HIV-1) subtype B sequences, a random starting topology could produce nucleotide substitution rates that were not statistically different than the true rates. Samples were isolated from 100 HIV-1 subtype B infected individuals from the United States and a 620 nt region of the env gene was sequenced for each sample. The sequence data were used to obtain a substitution model of sequence evolution specific for HIV-1 subtype B env by estimating nucleotide substitution rates and the site-to-site heterogeneity in 100 individuals from the United States. The method of estimating the model should provide users of large data sets with a way to quickly compute a model of sequence evolution, while the nucleotide substitution model we identified should prove useful in the phylogenetic analysis of HIV-1 subtype B env sequences. Received: 4 October 2000 / Accepted: 1 March 2001     

Heterogeneity and evolution rates of delta virus RNA sequences.

To investigate the geographical divergence of delta virus RNA sequences, 868 nucleotides (nt), including the delta antigen-coding region, were determined in isolates from two Japanese patients, M and S, by polymerase chain reaction and direct sequencing and compared with three previously reported nucleotide sequences. The sequence obtained for hepatitis delta virus RNA from patient M was approximately 92% identical to sequences previously obtained for two other strains of hepatitis delta virus, whereas the sequence of hepatitis delta virus RNA obtained from patient S was approximately 81% identical to the previously sequenced strains. This suggests that delta agent in Japan has a heterogeneous origin and the delta virus RNA sequence from Japanese patient S is the most divergent delta virus isolate yet analyzed. To study the evolution rate of delta virus RNA, viral isolates obtained 3 and 4 years apart from each of two patients were also sequenced. It was estimated that the substitution rate of viral RNA was 0.57 x 10(-3) nt per site per year in patient M and 0.64 x 10(-3) nt per site per year in patient S for the delta antigen gene.     

Equine infectious anemia virus genomic evolution in progressor and nonprogressor ponies

A primary mechanism of lentivirus persistence is the ability of these viruses to evolve in response to biological and immunological selective pressures with a remarkable array of genetic and antigenic variations that constitute a perpetual natural experiment in genetic engineering. A widely accepted paradigm of lentivirus evolution is that the rate of genetic variation is correlated directly with the levels of virus replication: the greater the viral replication, the more opportunities that exist for genetic modifications and selection of viral variants. To test this hypothesis directly, we examined the patterns of equine infectious anemia virus (EIAV) envelope variation during a 2.5-year period in experimentally infected ponies that differed markedly in clinical progression and in steady-state levels of viral replication as indicated by plasma virus genomic RNA assays. The results of these comprehensive studies revealed for the first time similar extents of envelope gp90 variation in persistently infected ponies regardless of the number of disease cycles (one to six) and viremia during chronic disease. The extent of envelope variation was also independent of the apparent steady-state levels of virus replication during long-term asymptomatic infection, varying from undetectable to 10(5) genomic RNA copies per ml of plasma. In addition, the data confirmed the evolution of distinct virus populations (genomic quasispecies) associated with sequential febrile episodes during acute and chronic EIA and demonstrated for the first time ongoing envelope variation during long-term asymptomatic infections. Finally, comparison of the rates of evolution of the previously defined EIAV gp90 variable domains demonstrated distinct differences in the rates of nucleotide and amino acid sequence variation, presumably reflecting differences in the ability of different envelope domains to respond to immune or other biological selection pressures. Thus, these data suggest that EIAV variation can be associated predominantly with ongoing low levels of virus replication and selection in target tissues, even in the absence of substantial levels of plasma viremia, and that envelope variation continues during all stages of persistent infection as the virus successfully avoids clearance by host defense mechanisms.     

The E1B 19-kilodalton protein is not essential for transformation of rodent cells in vitro by adenovirus type 5.

The newly constructed adenovirus type 5 mutant in1 carries a single AT base pair insertion immediately after nucleotide position 1715 in the E1B gene sequence which destroys the proximal AUG normally present in E1B messages and prevents production of intact E1B 19-kDa protein in infected cells. We have used in1, variants of in1 containing mutant alleles of viral genes known to enhance transformation frequency, and adenovirus type 5 mutant dl337 (S. Pilder, J. Logan, and T. Shenk, J. Virol. 52:664-671, 1984), in which the sequence between nucleotides 1770 and 1916 within the 19-kDa reading frame is deleted, to test the generally accepted hypothesis that this E1B protein is essential for the transformation of rodent cells and maintenance of the transformed phenotype. We find that these mutants transform rat embryo cells, rat kidney and mouse kidney primary cells, and cells of the 3Y1 rat line with decreased frequencies only when virus is added to these various cells at high input multiplicities of infection. In contrast, when lower doses of virus are used, the mutants transform with wild-type frequencies. Cells infected with higher doses of mutant virus show increased levels of DNA degradation and cell killing compared with those of cells infected with the same levels of wild-type virus, and these effects most likely contribute to the decreased transformation frequencies observed. On the basis of these results and the results of phenotypic analyses of numerous transformants, we propose that the E1B 19-kDa protein is not required for induction and/or maintenance of transformed-cell characteristics in rodent cells infected with adenovirus type 5.     

Phylogenetic evidence for the rapid evolution of human B19 erythrovirus

Human B19 erythrovirus is a ubiquitous viral pathogen, commonly infecting individuals before adulthood. As with all autonomous parvoviruses, its small single-stranded DNA genome is replicated with host cell machinery. While the mechanism of parvovirus genome replication has been studied in detail, the rate at which B19 virus evolves is unknown. By inferring the phylogenetic history and evolutionary dynamics of temporally sampled B19 sequences, we observed a surprisingly high rate of evolutionary change, at approximately 10(-4) nucleotide substitutions per site per year. This rate is more typical of RNA viruses and suggests that high mutation rates are characteristic of the Parvoviridae.     

Evolution of a persistent aphthovirus in cytolytic infections: partial reversion of phenotypic traits accompanied by genetic diversification.

Foot-and-mouth disease virus (FMDV) shows a dual potential to be cytolytic or to establish persistent infections in cell culture. FMDV R100, a virus rescued after 100 passages of carrier BHK-21 cells persistently infected with FMDV clone C-S8c1, showed multiple genetic and phenotypic alterations relative to the parental clone C-S8c1. Several FMDV R100 populations have been subjected to 100 serial cytolytic infections in BHK-21 cells, and the reversion of phenotypic and genetic alterations has been analyzed. An extreme temperature sensitivity of R100 reverted totally or partially in some passage series but not in others. The small-plaque morphology reverted to normal size in all cases. The hypervirulence for BHK-21 cells did not revert, and even showed an increase, upon cytolytic passage. Most of the mutations that had been fixed in the R100 genome during persistence did not revert in the course of cytolytic passages, but the extended polyribocytidylate tract of R100 (about 460 residues, versus 290 in C-S8c1) decreased dramatically in length, to the range of 220 to 260 residues in all passage series examined. In passages involving very large viral populations, a variant with two amino acid substitutions (L-144-->V and A-145-->P) next to the highly conserved Arg-Gly-Asp (RGD motif; positions 141 to 143) within the G-H loop of capsid protein VP1 became dominant. A clonal analysis allowed isolation of a mutant with the single replacement A-145-->P. Viral production and growth competition experiments showed the two variants to have a fitness very close to that of the parental virus. The results provide evidence that the repertoire of variants that could potentially become dominant in viral quasispecies may be influenced by the population size of the evolving virus. The net results of a series of persistent-infection passages followed by a series of cytolytic passages was progressive genomic diversification despite reversion or stasis of phenotypic traits. Implications for the evolution of RNA viruses are discussed.     

Molecular evolution of hepatitis B virus over 25 years

Determining the longitudinal molecular evolution of hepatitis B virus (HBV) is difficult due to HBV's genomic complexity and the need to study paired samples collected over long periods of time. In this study, serial samples were collected from eight hepatitis B virus e antigen-negative asymptomatic carriers of HBV genotype B in 1979 and 2004, thus providing a 25-year period to document the long-term molecular evolution of HBV. The rate, nature, and distribution of mutations that emerged over 25 years were determined by phylogenetic and linear regression analysis of full-length HBV genome sequences. Nucleotide hypervariability was observed within the polymerase and pre-S/S overlap region and within the core gene. The calculated mean number of nucleotide substitutions/site/year (7.9 x 10(-5)) was slightly higher than published estimates (1.5 x 10(-5) to 5 x 10(-5)). Nucleotide changes in the quasispecies population did not significantly alter the molecular evolutionary rate, based on linear regression analysis of evolutionary distances among serial clone pre-S region sequences. Therefore, the directly amplified or dominant sequence was sufficient to estimate the putative molecular evolutionary rate for these long-term serial samples. On average, the ratio of synonymous (dS) to nonsynonymous (dN) substitutions was highest for the polymerase-coding region and lowest for the core-coding region. The low dS/dN ratios observed within the core suggest that selection occurs within this gene region, possibly as an immune evasion strategy. The results of this study suggest that HBV sequence divergence may occur more rapidly than previously estimated, in a host immune phase-dependent manner.     

Dynamics of deletion genotypes in an experimental insect virus population

Defective viruses, that are deficient in certain essential genes, are maintained in the population by trans-complementation, exploiting the gene products of complete genotypes in co-infected cells. This process becomes prevalent only when cells are frequently infected by several virus particles, and only then will the fitness of defective viruses be subjected to frequency-dependent selection. Deletion variants that are not infectious per os are present in a multicapsid nucleopolyhedrovirus (SfMNPV, Baculoviridae) that infects the fall army worm, Spodoptera frugiperda. These variants enhance the pathogenicity and, therefore, the likelihood of transmission of the virus when co-infecting cells with complete genotypes, resulting in occlusion bodies (OBs) that may contain both genotypes co-occluded. Mixtures of complete (B) and defective (C) variants in ratios of 90% B+10% C, 50% B+50% C and 10% B+90% C were used to inoculate by injection S. frugiperda larvae. Viral OBs extracted from diseased insects were subjected to four or five successive rounds of per os infection. Following successive passages, genotype frequencies in all three experimental populations converged to a single equilibrium frequency comprising approximately 20% of deletion genotype C and approximately 80% of complete genotype B. This mirrors the relative proportions of deletion (22%) and complete (78%) genotypes observed in the wild-type SfMNPV population. The pathogenicity of experimental populations at the final passage was not significantly different from that of the wild-type isolate. In contrast, OBs of all genotype mixtures were significantly more pathogenic than OBs of genotype B alone. A population genetics model, in which virus populations were assigned linear frequency-dependent transmissibility values, was in remarkably close agreement to empirical data. Clearly, non-infectious deletion variants can profoundly affect the likelihood of transmission and the genetic structure and stability of virus populations.     

The intricate interplay between RNA viruses and NF-κB

RNA viruses have rapidly evolving genomes which often allow cross-species transmission and frequently generate new virus variants with altered pathogenic properties. Therefore infections by RNA viruses are a major threat to human health. The infected host cell detects trace amounts of viral RNA and the last years have revealed common principles in the biochemical mechanisms leading to signal amplification that is required for mounting of a powerful antiviral response. Components of the RNA sensing and signaling machinery such as RIG-I-like proteins, MAVS and the inflammasome inducibly form large oligomers or even fibers that exhibit hallmarks of prions. Following a nucleation event triggered by detection of viral RNA, these energetically favorable and irreversible polymerization events trigger signaling cascades leading to the induction of antiviral and inflammatory responses, mediated by interferon and NF-κB pathways. Viruses have evolved sophisticated strategies to manipulate these host cell signaling pathways in order to ensure their replication. We will discuss at the examples of influenza and HTLV-1 viruses how a fascinating diversity of biochemical mechanisms is employed by viral proteins to control the NF-κB pathway at all levels.