Free access
Research Article
15 June 2004

Clades of Adeno-Associated Viruses Are Widely Disseminated in Human Tissues


The potential for using Adeno-associated virus (AAV) as a vector for human gene therapy has stimulated interest in the Dependovirus genus. Serologic data suggest that AAV infections are prevalent in humans, although analyses of viruses and viral sequences from clinical samples are extremely limited. Molecular techniques were used in this study to successfully detect endogenous AAV sequences in 18% of all human tissues screened, with the liver and bone marrow being the most predominant sites. Sequence characterization of rescued AAV DNAs indicated a diverse array of molecular forms which segregate into clades whose members share functional and serologic similarities. One of the most predominant human clades is a hybrid of two previously described AAV serotypes, while another clade was found in humans and several species of nonhuman primates, suggesting a cross-species transmission of this virus. These data provide important information regarding the biology of parvoviruses in humans and their use as gene therapy vectors.
Adeno-associated virus (AAV) is a member of the genus Dependovirus, which lies within the Parvoviridae family (17). An interest in this family of viruses has been stimulated because of their potential use as gene transfer vectors (14).
Little is known about the biology of AAV infections, although a significant proportion of humans and nonhuman primates have antibodies in their blood that react to some of the six existing serotypes of AAV (5, 7). This suggests that primates are hosts for infection with AAV, although the clinical sequelae of these infections have yet to be identified.
The study of AAV has been limited to the previously described six serotypes, of which five were isolated as contaminants in laboratory preparations of adenoviruses (1, 3, 16). Our lack of understanding of AAV clinical infections has complicated the search for clinical isolates of the virus. Members of our laboratory recently described a strategy for evaluating latent or persistent AAV genomes from tissues of asymptomatic nonhuman primates through the use of PCR. These studies led to the discovery of two novel AAV serotypes, called AAV7 and AAV8, that have improved properties as vectors for gene therapy (10). In nonhuman primates, AAV sequences were quite prevalent and heterogenous (9).
The goal of this study was to determine if latent AAVs exist in humans, and if so, to characterize their structural, serologic, and functional properties.


Collection of primate tissues.

Our sources of nonhuman primate tissues were described previously (9). Human tissues were collected under two independent IRB protocols approved by the Institutional Review Board of the University of Pennsylvania from either surgical procedures, postmortem examinations, or organ donors through two major national human tissue providers, the Cooperative Human Tissue Network and the National Disease Research Interchange. The human tissues used for this study were comprised of 18 different tissue types that included the colon, liver, lung, spleen, kidney, brain, small bowel, bone marrow, heart, lymph nodes, skeletal muscle, ovary, pancreas, stomach, esophagus, cervix, testes, and prostate. The tissue samples came from a diverse group of individuals of different genders, races (Caucasian, African American, Asian, and Hispanic), and ages (23 to 83 years). Among the 259 samples analyzed from 250 individuals, approximately 28% of the tissues were associated with pathology.

Detection and isolation of AAV sequences.

Total cellular DNAs were extracted from human and nonhuman primate tissues as described previously (10). The molecular prevalence and tissue distribution of AAVs in humans were determined by either signature or full-length cap PCRs using primers and conditions that were similar to those used for nonhuman primate analyses (9, 10). The same PCR cloning strategy used for the isolation and characterization of an expanded family of AAVs in nonhuman primates was deployed for the isolation of AAVs from selected human tissues (9, 10). A total of 67 capsid clones isolated from human tissues were characterized (hu.1 to hu.67). From nonhuman primate tissues, 86 cap clones were sequenced, among which 70 clones were from rhesus macaques, 6 clones were from cynomolgus macaques, 3 clones were from pigtailed macaques, 2 clones were from a baboon, and 5 clones were from a chimpanzee.

Computational analysis of primate AAV sequences.

From all sequence contigs, VP1 open reading frames were analyzed. VP1 protein sequences were aligned with ClustalX (22) and an in-frame DNA alignment was produced with the BioEdit (11) software package. Phylogenies were inferred with the MEGA v2.1 and TreePuzzle packages. Neighbor-joining, maximum parsimony (18), and maximum likelihood (21) algorithms were used to confirm similar clusterings of sequences in monophylic groups. Clades were then defined from a neighbor-joining phylogenetic tree of all protein sequences. The amino acid distances were estimated by making use of Poisson correction. The statistical robustness of the analysis was estimated by bootstrapping with 1,000 replicates. Sequences were considered monophylic when they had a connecting node within a genetic distance of 0.05. A group of sequences originating from three or more sources was considered a clade. Homoplasy was screened for by implementation of the split decomposition algorithm (2). Splits that were picked up in this manner were then further analyzed for recombination by use of the Bootscan algorithm in the Simplot software (20). A sliding window of 400 nucleotides (nt) (10 nt/step) was used to obtain 100 bootstrap replicate neighbor-joining trees. Subsequently, split decomposition and neighbor-joining phylogenies were inferred from the putative recombination fragments. A significant improvement of bootstrap values, a reduction of splits, and a regrouping of the hybrid sequences with their parents were considered the criteria for recombination.

Evaluation of primate AAVs as gene transfer vectors.

Primate AAV vectors were produced by a previously described transencapsidation method (10). Vectors expressing enhanced green fluorescent protein (EGFP) were used to examine their in vitro transduction efficiencies in 84-31 cells and to study their serological properties. For in vivo studies, human α-antitrypsin (A1AT) was selected as a sensitive and quantitative reporter gene for the vectors and was expressed under the control of a cytomegalovirus-enhanced chicken β-actin promoter. Four to 6-week-old NCR nude mice were treated with novel AAV vectors at a dose of 1011 genome copies per animal through intraportal, intratracheal, and intramuscular injections for liver-, lung-, and muscle-directed gene transfer, respectively. Serum samples were collected at different time points after gene transfer and A1AT concentrations were determined by an enzyme-linked immunosorbent assay.

Nucleotide sequence accession numbers.

New AAV sequences from the PCR clones presented in this paper have been submitted to GenBank. The accession numbers are AY530553 to AY530629 .


Human tissues were obtained from a variety of sources, and DNAs were evaluated for endogenous AAV sequences by PCRs with oligonucleotides specific to homologous regions of the cap gene. Figure 1 summarizes a portion of the results of this screen of 259 human samples of 18 different tissue types derived from 250 individuals. The data were compared to similar studies using an expanded pool of tissues from rhesus monkeys, cynomolgus and pigtailed macaques, baboons, and chimpanzees. The prevalence of AAV sequences in human and nonhuman primate tissues was similar (19 and 18%, respectively). Livers and spleens were the predominant sites of AAV infection in both human and nonhuman primates, although endogenous sequences were also frequently found in colons and bone marrow from humans and lymph nodes from nonhuman primates. Quantitative PCR studies indicated that endogenous AAV in human tissues is present in low quantities and unlikely to be present as a result of germ line transmission (data not shown).
To better understand the origins and consequences of endogenous AAVs in humans, we attempted to recover and fully sequence full-length cap structures from human tissues; isolates from nonhuman primates, in addition to what was previously described, were also included. A total of 108 new and unique isolates (from 55 human and 53 nonhuman primates) were identified (clones from the same individual with fewer than four amino acid differences were deemed redundant and eliminated from the analysis).
This pool of primate AAV cap sequences was analyzed for phylogenetic relationships by using a variety of computational approaches. Sequences were aligned with the ClustalX1.81 program and phylogenies were assessed by the neighbor-joining, maximum parsimony, and maximum likelihood algorithms (18, 19, 21). Each method yielded a similar clustering of sequences. The phylogeny of AAV was further evaluated for evidence of recombination through a sequential analysis of split decomposition (2) and bootscanning (20). Split decomposition analysis depicts parallel events in a set of sequences with a tree-like network rather than a bifurcating tree. The bootscanning algorithm then further verifies these putative recombination events by visualizing the mosaic structure of a given sequence. A number of different cap sequences amplified from eight different human subjects showed phylogenetic relationships to AAV2 (5′) and AAV3 (3′) around a common breakpoint at position 1400 of the cap DNA sequence, consistent with recombination and the formation of a hybrid virus (Fig. 2). This is the general region of the cap gene in which recombination was detected in isolates from a mesenteric lymph node of a rhesus macaque (9).
The phylogenetic analyses were repeated, excluding the clones that were positively identified as hybrids. In this analysis, goose and avian AAVs were included as outgroups (6). Figure 3 summarizes a neighbor-joining tree; similar relationships were obtained by maximum parsimony and maximum likelihood analyses. These analyses demonstrated 11 phylogenetic groups, which are summarized in Table 1. When a group contained nonredundant but phylogenetically similar members from three or more sources, it was called a clade; otherwise, it was called a clone or set of clones. The clades were defined as clades A to F, as shown in Table 1. The five categories of clones are as follows: the previously described AAV3 and AAV5 clones that are clearly distinct from one another but were not detected in our screen; a group of three similar clones from a rhesus macaque that are closely related to AAV4; two similar clones, represented by rh.8, from different rhesus macaques that are not related to previously described AAV serotypes; and a single unique clone from chimpanzees, called ch.5. The previously described AAV1 and AAV6 clones are members of a single clade for which four isolates were recovered from three humans. The clades representing AAV2 and the AAV2-AAV3 hybrid are the most abundant of those found in humans (22 isolates from 12 individuals for AAV2 and 17 isolates from 8 individuals for the AAV2-AAV3 hybrid). A clade containing AAV7 is unique to rhesus and cynomolgus macaques, with 15 members being isolated from 10 different animals. The clade containing AAV8 is interesting because it is found in both human and nonhuman primates: 9 isolates were recovered from 7 humans and 21 isolates were obtained from 9 different nonhuman primates, including rhesus macaques, a baboon, and a pigtail monkey.
The last clade was derived from isolates from three humans and did not contain a previously described serotype. Polyclonal antisera were generated against a representative member of this clade, and a comprehensive study of serologic cross-reactivity between the previously described serotypes was performed (Table 2). The isolate from this clade was serologically distinct from the other known serotypes and therefore the clade was called the AAV9 clade.
Additional experiments were performed to evaluate the relationship of phylogenetic relatedness to function, as measured by serologic activity and tropism. Polyclonal antisera generated against the nine known serotypes were used to evaluate cross-neutralization (Table 2). For the purposes of discussion, we defined a new serotype as one for which the neutralization titer by heterologous sera was at least 16-fold less than the neutralization titer against the homologous vector in reciprocal titrations. These data confirmed the phylogenetic groupings of the different clones and clades except for an unanticipated serological reactivity of the structurally distinct AAV5 and AAV1 serotypes (the ratios of heterologous to homologous titers were 1/4 and 1/8, respectively, in reciprocal titrations). It should be noted that the previously described AAV1 and AAV6 serotypes do not segregate by either their phylogeny (Fig. 3) or their serology (the ratios of heterologous to homologous titers were 1/2 and 1/4, respectively, in reciprocal titrations).
The biological tropisms of AAVs were studied by generating transencapsidated vectors in which recombinant AAV2 genomes expressing either GFP or the secreted reporter gene A1AT were packaged with capsids derived from the various clones or clades. The vectors were evaluated for their transduction efficiency in vitro, based on GFP transduction, and their transduction efficiency in vivo in the liver, muscle or lung. Vectors expressing EGFP were used to examine their in vitro transduction efficiencies in 84-31 cells and to study their serological properties. For in vivo studies, human A1AT was selected as a sensitive and quantitative reporter gene for the vectors and was expressed under the control of the cytomegalovirus-enhanced chicken β-actin promoter. Four to 6-week-old NCR nude mice were treated with novel AAV vectors at a dose of 1011 genome copies per animal through intraportal, intratracheal, and intramuscular injections for liver-, lung-, and muscle-directed gene transfers, respectively. Serum samples were collected at different time points after the gene transfer, and A1AT concentrations were determined by an enzyme-linked immunosorbent assay. A representative set of assay results is shown in Fig. 4.
In order to compare unique profiles of transduction, we developed a grading system to characterize the relative transduction efficiency of each in vitro and in vivo model (from 0 [lowest] to 3 [highest]). The cumulative functional difference between two vectors with capsids A and B is the sum of the absolute values of differences between the individual assays as follows: cumulative functional difference = (in vitro A − in vitro B) + (liver A − liver B) + (lung A − lung B) + (muscle A − muscle B). Smaller cumulative functional differences indicate similar profiles with regard to transduction efficiency. Table 3 summarizes the cumulative functional difference scores as well as the % differences in VP1 amino acid sequences in pairwise comparisons.
Unique profiles of biological activity, in terms of the efficiency of gene transfer, were demonstrated for the different clones and clades of AAVs, with substantial concordance between members of a set of clones or a clade (data not shown). This suggests that biological pressures drive the evolution of AAVs.
Our studies point out a number of issues that are relevant to the study of parvoviruses in humans. The prevalence of endogenous AAV sequences in a wide array of human tissues suggests that natural infections with this group of viruses are quite common. The wide tissue distribution of viral sequences and their frequent detection in the liver, spleen, and gut suggest that transmission may occur via the gastrointestinal tract and that viremia may be a feature of the infection. Some earlier reports also documented the detection of AAV sequences in the human female genital tract and suggested that sexual contact could be another route of transmission (8, 23). However, the clinical consequences of infection with AAV have yet to be delineated.
An inspection of the topology of the phylogenetic analysis revealed insight into the relationship between the evolution of the virus and its host restriction. The entire genus Dependovirus appears to be derived from avian AAVs, consistent with the work of Lukashov and Goudsmit (15). After the emergence of AAV4 and AAV5, the family diverged into two monophylic groups (Fig. 3), with one containing clades that are specific to humans (clades A, B, and C) and the other comprised of a mixture of clades that were isolated exclusively from humans (clade F), exclusively from nonhuman primates (clade D), or from both human and nonhuman primates (clade E).
The presence of latent AAVs that are widely disseminated throughout human and nonhuman primates and their apparent predisposition to recombine and to cross species barriers raise important issues. This combination of events has the potential to lead to the emergence of new infectious agents with modified virulence. Assessments of this potential are confounded by the fact that the clinical sequelae of AAV infections in primates have yet to be defined. In addition, the high prevalence of AAV sequences in the liver may contribute to dissemination of the virus in the human population in the setting of allogeneic and xenogeneic liver transplantation. Finally, the finding of endogenous AAVs in humans has implications for the use of AAV for human gene therapy. The fact that wild-type AAV is so prevalent in primates without ever being associated with a malignancy suggests that it is not particularly oncogenic. In fact, the expression of AAV rep genes has been shown to suppress transformation (12). A potential complication of AAV gene therapy, however, could be recombination between the vector and endogenous genomes. This could lead to swapping of the inverted terminal repeats, rearrangement in the transgene cassettes, a loss of regulatory elements for regulated gene expression, and other effects.
FIG. 1.
FIG. 1. Tissue distribution of AAV sequences in primate tissues. Total cellular DNAs were extracted from 10 major tissues, including the blood, from nonhuman primates (top) and humans (bottom) and were subjected to PCR analysis for AAV sequences. The x axis indicates the tissues analyzed and the y axis represents the frequencies of detection of AAV sequences. The total number of samples examined for each tissue of human or nonhuman primate origin is listed at the top of each bar.
FIG. 2.
FIG. 2. Recombination analysis of AAV2-AAV3 hybrid clade representative hu.2. (Top) Bootscanning analysis of hu.2 versus representatives of all other clades showing a distinct phylogenetic relatedness between the 5′ (A) and the 3′ (B) segments (indicated below) of the entire VP1 capsid sequence. (Bottom) Split decomposition analysis of AAV2, AAV3, and the AAV2-AAV3 hybrid group representative hu.2, with AAV8 used as an outgroup. The left panel shows the analysis of the entire VP1 sequence, while on the right the conflicting phylogenies are resolved for segments A and B. The hamming distances are indicated on the splits. All branches were supported by a minimum of 85% bootstrap values (n = 1,000).
FIG. 3.
FIG. 3. Neighbor-joining phylogenies of the VP1 protein sequence of primate AAVs. Major nodes with bootstrap values of <75 are indicated with an “X.” A goose parvovirus and an avian AAV (6) were used as the outgroup. Clades are indicated by name and by vertical lines to the right of the taxa from which they are made. The nomenclature for the taxa is either the serotype name or a reference to the species source (hu, human; rh, rhesus macaque; cy, cynomolgus macaque; bb, baboon; pi, pigtailed macaque; ch, chimpanzee), followed by a number indicating the order in which they were sequenced. Clade C was identified and positively determined to have originated through the recombination of known clades. The AAV2-AAV3 hybrid clade originated after one recombination event, and its unrooted neighbor-joining phylogeny is shown.
FIG. 4.
FIG. 4. Evaluation of gene transfer efficiency of novel primate AAV-based vectors in vitro and in vivo. AAV vectors were pseudotyped with capsids as indicated and measured as percentages of GFP+ cells. Analyses of in vitro transduction with GFP vectors and in vivo transduction with A1AT vectors (measured as serum A1AT levels, in micrograms per milliliter) were performed as described in Materials and Methods. For in vivo comparisons, the serum A1AT level is indicated at the top of each bar, representing the gene transfer efficiency of each AAV vector.
TABLE 1. Classification of number of sources per clade or clonea
Clade or cloneRepresentative membersNo. of sequences for speciesb     Reference for published sequence
  HumansRhesus macaquesCynomolgus macaquesBaboonsChimpanzeesPigtail macaques 
    AAAV1 (AAV6)3 (4)     13
    BAAV212 (22)     1
    CAAV2-AAV3 hybrid8 (17)      
    DAAV7 5 (10)5 (5)   10
    EAAV87 (9)7 (16) 1 (2) 1 (3)10
    FAAV93 (3)      
    AAV3       4
    AAV4  1 (3)    16
    AAV5       3
    ch.5     1 (1)  
    rh.8  2 (2)     
The species origins of the six AAV clades and five individual AAV clones (or sets of clones) are represented.
Numbers in parentheses indicate the total numbers of nonredundant sequences isolated.
TABLE 2. Serologic evaluation of novel AAV vectorsa
Vector used for rabbit immunizationTiter for vector pseudotype        
AAV2/11/163,840No NABNo NABNo NAB1/40,9601/40,9601/40No NABNo NAB
AAV2/21/801/81,9201/5,1201/20No NAB1/801/401/40No NAB
AAV2/31/1,2801/2,5601/40,9601/201/401/2,5601/1,2801/1,280No NAB
AAV2/41/20No NABNo NAB1/1,2801/40No NABNo NABNo NAB1/40
AAV2/51/20,480No NAB1/80No NAB1/163,8401/5,1201/40No NABNo NAB
AAV2/61/81,920No NAB1/6401/401/401/327,6801/40No NAB1/40
AAV2/71/1,2801/6401/1,2801/20No NAB1/1,2801/163,8401/5,1201/80
AAV2/81/201/1,2801/1,280No NAB1/20No NAB1/6401/327,6801/2,560
AAV2/9No NABNo NABNo NABNo NABNo NABNo NAB1/201/6401/20,480
Rabbit polyclonal antibodies against AAV serotypes 1 to 9 were generated. To test neutralizing and cross-neutralizing reactivities of antisera to each serotype of the AAVCMVEGFP vector, we preincubated the vectors with homologous and heterologous antisera and then incubated them with indicator 84-31 cells for 48 to 72 h. The transduction of 84-31 cells by AAVCMVEGFP vectors was assessed under a UV microscope. Neutralization titers were defined as described previously (10). NAB, neutralizing antibody.
TABLE 3. Functional and structural differences of primate AAVsa
VirusCumulative difference in gene transfer scores or % cap amino acid differences (in bold)         
To examine differences in structure and function among primate AAVs, we selected one representative member from each primate AAV clade for comparison. For instance, the data obtained from AAV1 were used to represent clade A, followed by AAV2 for clade B, rh.34 for the AAV4 clade, AAV7 for clade D, AAV8 for clade E, and AAV hu.14 for clade F. AAV5, AAVch.5, and AAVrh.8 stand as single AAV serotypes or genotypes for the comparison. For functional analysis, the in vitro transduction of different AAVCMVEGFP vectors was measured in 84-31 cells that were seeded in a 96-well plate and infected with pseudotyped AAVCMVEGFP vectors at a multiplicity of infection of 104 GC per cell. The relative EGFP transduction efficiency was scored as 0, 1, 2, or 3, corresponding to 0 to 10%, 10 to 30%, 30 to 70%, and 70 to 100% green cells, as estimated by using a UV microscope at 48 h postinfection. The in vivo potency of different AAVCBhA1AT vectors was evaluated in NCR nude mice for liver-, lung-, and muscle-directed gene transfer at a dose of 1011 GC per animal and was scored as 0, 1, 2, or 3 relative to different serum A1AT levels at day 28 after gene transfer, depending on the route of vector administration (for the liver, 0 = A1AT of <400 ng/ml, 1 = A1AT of 400 to 1,000 ng/ml, 2 = A1AT of 1,000 to 10,000 ng/ml, and 3 = A1AT of >10,000 ng/ml; for the lung, 0 = A1AT of <200 ng/ml, 1 = A1AT of 200 to 1,000 ng/ml, 2 = A1AT of 1,000 to 10,000 ng/ml, and 3 = A1AT of >10,000 ng/ml; for muscle, 0 = A1AT of <100 ng/ml, 1 = A1AT of 100 to 1,000 ng/ml, 2 = A1AT of 1,000 to 10,000 ng/ml, and 3 = A1AT of >10,000 ng/ml). Cumulative differences of the gene transfer scores were calculated as described in the text. ND, not determined. The percentages of difference in cap structure were determined by dividing the number of amino acid differences after a pairwise deletion of gaps by 750, the length of the VP1 protein sequence alignment.


This work was supported by NIH (NIDDK P30 DK47757 and NHLBI P01 HL59407), the Cystic Fibrosis Foundation, the Juvenile Diabetes Research Foundation and GlaxoSmithKline Pharmaceuticals.
J.M.W. held equity in Targeted Genetics Corp. at the time of this study.


Atchison, R. W., B. C. Casto, and W. M. Hammon. 1965. Adenovirus-associated defective virus particles. Science194:754-756.
Bandelt, H. J., and A. W. Dress. 1992. Split decomposition: a new and useful approach to phylogenetic analysis of distance data. Mol. Phylogenet. Evol.1:242-252.
Bantel-Schaal, U., and H. Zur Hausen. 1984. Characterization of the DNA of a defective human parvovirus isolated from a genital site. Virology134:52-63.
Blacklow, N. R., M. D. Hoggan, and W. P. Rowe. 1967. Isolation of adenovirus-associated viruses from man. Proc. Natl. Acad. Sci. USA58:1410-1415.
Blacklow, N. R., M. D. Hoggan, and W. P. Rowe. 1968. Serologic evidence for human infection with adenovirus-associated viruses. J. Natl. Cancer Inst.40:319-327.
Bossis, I., and J. A. Chiorini. 2003. Cloning of an avian adeno-associated virus (AAAV) and generation of recombinant AAAV particles. J. Virol.77:6799-6810.
Chirmule, N., K. Propert, S. Magosin, Y. Qian, R. Qian, and J. Wilson. 1999. Immune responses to adenovirus and adeno-associated virus in humans. Gene Ther.6:1574-1583.
Friedman-Einat, M., Z. Grossman, F. Mileguir, Z. Smetana, M. Ashkenazi, G. Barkai, N. Varsano, E. Glick, and E. Mendelson. 1997. Detection of adeno-associated virus type 2 sequences in the human genital tract. J. Clin. Microbiol.35:71-78.
Gao, G., M. R. Alvira, S. Somanathan, Y. Lu, L. H. Vandenberghe, J. J. Rux, R. Calcedo, J. Sanmiguel, Z. Abbas, and J. M. Wilson. 2003. Adeno-associated viruses undergo substantial evolution in primates during natural infections. Proc. Natl. Acad. Sci. USA100:6081-6086.
Gao, G. P., M. R. Alvira, L. Wang, R. Calcedo, J. Johnston, and J. M. Wilson. 2002. Novel adeno-associated viruses from rhesus monkeys as vectors for human gene therapy. Proc. Natl. Acad. Sci. USA99:11854-11859.
Hall, T. A. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp. Ser.41:95-98.
Hermonat, P. L. 1989. The adeno-associated virus Rep78 gene inhibits cellular transformation induced by bovine papillomavirus. Virology172:253-261.
Hoggan, M. D., N. R. Blacklow, and W. P. Rowe. 1966. Studies of small DNA viruses found in various adenovirus preparations: physical, biological, and immunological characteristics. Proc. Natl. Acad. Sci. USA55:1467-1474.
Kotin, R. M. 1994. Prospects for the use of adeno-associated virus as a vector for human gene therapy. Hum. Gene Ther.5:793-801.
Lukashov, V. V., and J. Goudsmit. 2001. Evolutionary relationships among parvoviruses: virus-host coevolution among autonomous primate parvoviruses and links between adeno-associated and avian parvoviruses. J. Virol.75:2729-2740.
Mayor, H. D., and J. L. Melnick. 1966. Small deoxyribonucleic acid-containing viruses (picornavirus group). Nature210:331-332.
Muzyczka, N., and K. I. Berns. 2001. Parvoviridae: the viruses and their replication, p. 2327-2359. In D. M. Knipe et al. (ed.), Fields virology, 4th ed. Lippincott Williams & Wilkins, Philadelphia, Pa.
Nei, M., and S. Kumar. 2000. Molecular evolution and phylogenetics. Oxford University Press, New York, N.Y.
Saitou, N., and M. Nei. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Evol.4:406-425.
Salminen, M. O., J. K. Carr, D. S. Burke, and F. E. McCutchan. 1995. Identification of breakpoints in intergenotypic recombinants of HIV type 1 by bootscanning. AIDS Res. Hum. Retrovir.11:1423-1425.
Schmidt, H. A., K. Strimmer, M. Vingron, and A. von Haeseler. 2002. TREE-PUZZLE: maximum likelihood phylogenetic analysis using quartets and parallel computing. Bioinformatics18:502-504.
Thompson, J. D., T. J. Gibson, F. Plewniak, F. Jeanmougin, and D. G. Higgins. 1997. The CLUSTAL X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res.25:4876-4882.
Walz, C. M., T. R. Anisi, J. R. Schlehofer, L. Gissmann, A. Schneider, and M. Muller. 1998. Detection of infectious adeno-associated virus particles in human cervical biopsies. Virology247:97-105.

Information & Contributors


Published In

cover image Journal of Virology
Journal of Virology
Volume 78Number 1215 June 2004
Pages: 6381 - 6388
PubMed: 15163731


Received: 8 December 2003
Accepted: 13 February 2004
Published online: 15 June 2004


Request permissions for this article.



Guangping Gao
Gene Therapy Program, Division of Medical Genetics, Department of Medicine, University of Pennsylvania School of Medicine, and The Wistar Institute, Philadelphia, Pennsylvania 19104
Luk H. Vandenberghe
Gene Therapy Program, Rega Institute for Medical Research, Katholieke Universiteit Leuven, B-3000 Louvain, Belgium
Mauricio R. Alvira
Gene Therapy Program, Division of Medical Genetics, Department of Medicine, University of Pennsylvania School of Medicine, and The Wistar Institute, Philadelphia, Pennsylvania 19104
You Lu
Gene Therapy Program, Division of Medical Genetics, Department of Medicine, University of Pennsylvania School of Medicine, and The Wistar Institute, Philadelphia, Pennsylvania 19104
Roberto Calcedo
Gene Therapy Program, Division of Medical Genetics, Department of Medicine, University of Pennsylvania School of Medicine, and The Wistar Institute, Philadelphia, Pennsylvania 19104
Xiangyang Zhou
Gene Therapy Program, Division of Medical Genetics, Department of Medicine, University of Pennsylvania School of Medicine, and The Wistar Institute, Philadelphia, Pennsylvania 19104
James M. Wilson [email protected]
Gene Therapy Program, Division of Medical Genetics, Department of Medicine, University of Pennsylvania School of Medicine, and The Wistar Institute, Philadelphia, Pennsylvania 19104

Metrics & Citations


Note: There is a 3- to 4-day delay in article usage, so article usage will not appear immediately after publication.

Citation counts come from the Crossref Cited by service.


If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.

View Options

Figures and Media






Share the article link

Share with email

Email a colleague

Share on social media

American Society for Microbiology ("ASM") is committed to maintaining your confidence and trust with respect to the information we collect from you on websites owned and operated by ASM ("ASM Web Sites") and other sources. This Privacy Policy sets forth the information we collect about you, how we use this information and the choices you have about how we use such information.
FIND OUT MORE about the privacy policy