Previous Article | Next Article ![]()
Journal of Virology, December 2005, p. 14781-14792, Vol. 79, No. 23
0022-538X/05/$08.00+0 doi:10.1128/JVI.79.23.14781-14792.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Ryan L. Jensen,
,
Bruce C. Schnepp,
Mary J. Connell,
Richard Shell,
Thomas J. Sferra,
Jeffrey S. Bartlett,
K. Reed Clark, and
Philip R. Johnson*
Center for Gene Therapy, Columbus Children's Research Institute, Columbus Children's Hospital, and Department of Pediatrics, College of Medicine and Public Health, The Ohio State University, Columbus, Ohio 43205
Received 20 June 2005/ Accepted 7 September 2005
|
|
|---|
98% amino acid identity), while the single spleen isolate was intermediate between serotypes 2 and 3. Comparison to the available AAV2 crystal structure revealed that the majority of the amino acid substitutions mapped to surface-exposed hypervariable domains. To further characterize the AAV capsid structure in these samples, we used a novel linear rolling-circle amplification method to amplify episomal AAV DNA and isolate infectious molecular clones from several human tissues. Serotype 2-like viruses were generated from these DNA clones and interestingly, failed to bind to a heparin sulfate column. Inspection of the capsid sequence from these two clones (and the other six AAV2-like isolates) revealed that they lacked arginine residues at positions 585 and 588 of the capsid protein, which are thought to be essential for interaction with the heparin sulfate proteoglycan coreceptor. These data provide a framework with which to explore wild-type AAV persistence in vivo and provide additional tools to further define the biodistribution and form of AAV in human tissues. |
|
|---|
Recently, when viral genomes of unexpected diversity were recovered from a surprising array of tissues and organs, it became clear that AAV infection of primates (human and nonhuman) is much more complex than previously appreciated (11, 13). The genetic diversity observed in these studies suggested that multiple genotypes (many more than the historically accepted five serotypes) of AAV circulate in humans and monkeys. Furthermore, the data implied that whatever the portal of entry, AAV can become widely disseminated following primary infection. The implications of these findings for recombinant AAV gene transfer vectors are unclear, but the immunologic and genetic influence of prior infection on recombinant AAV-mediated gene transfer must be considered.
To extend our understanding of wild-type AAV infection in humans, we set out to characterize AAV genomes directly out of freshly acquired human tissues. Because our assumption was that many (if not all) AAV infections begin in the respiratory tract in children as they are concurrently infected with adenoviruses, we collected tonsils and adenoids from pediatric subjects undergoing surgical excision in an outpatient surgery center. We demonstrated that 7% of these samples contained wild-type AAV DNA. In a follow-on set of experiments, we obtained a range of archived, frozen normal tissues from a repository and showed that 3% of these samples also contained wild-type AAV DNA. Analysis of the cap gene sequences from all samples revealed that most of the isolates were closely related to AAV serotype 2; a single isolate shared significant homology with serotypes 2 and 3. Interestingly, none of the isolates in our study were predicted to bind heparin sulfate, suggesting that this receptor is not necessary for wild-type infection in humans.
|
|
|---|
Adenovirus isolation from tonsils and adenoids. Freshly collected tissue from each subject (n = 101) was minced, mixed, and divided into three aliquots. Two aliquots were frozen for DNA isolation (see below) while the third was further minced and sieved through a cell strainer. The sieved material was adsorbed onto A459 cells cultured in Dulbecco's modified Eagle's medium containing 20% fetal calf serum. Cultures were maintained at 37°C in 5% CO2 for 14 days and then blind passed into fresh A549 cells and monitored for an additional 14 days. Cultures were scored positive or negative solely on the appearance of characteristic cytopathic effects.
DNA isolation from human tissue samples. Freshly thawed tissue (0.2 to 0.5 g) was digested for 15 h in 3 ml of digestion buffer (10 mM Tris, pH 8; 100 mM NaCl; 0.5% sodium dodecyl sulfate; 25 mM EDTA) supplemented with 2 mg/ml proteinase K at 55°C with constant agitation. DNA was extracted twice with phenol-chloroform/isoamyl alcohol (25:24:1) and RNA removed by DNase-free RNase A digestion (20 µg/ml) (QIAGEN Inc.) at 37°C for 30 min. RNase A was removed by two sequential phenol-chloroform-isoamyl alcohol extractions. A final chloroform extraction was performed, followed by DNA ethanol precipitation using 0.3 M sodium acetate (pH 5.2) and 2 volumes of ethanol. The DNA pellet was air-dried and suspended in 1 ml of 10 mM Tris, pH 8.0, and DNA concentration determined by A260.
AAV capsid PCR. Nested PCR was used to amplify a 255-bp conserved region of the AAV cap gene. An initial round of PCR was performed on genomic DNA (100 ng) using a conserved degenerate primer set (CapSS2978: 5'-GGYGCCGAYGGAGTGGGYARTKCC-3' and Cap 18S: 5'-GAWKCCCCARTWGTTGTTRATGAGTC-3'), where Y is C or T, R is A or G, K is G or T, W is A or T, M is A or C, S is C or G, I is inosine. One µl of the first-round PCR served as the template for the second round using nested, degenerate PCR primers (Cap 19S: 5'-GYARTKCCTCRGGWRATTGGCA-3' and CapSS3189: 5'-GATGAGTCKYTGCCAGTCWCGKGG-3'). Reaction components for both rounds were 400 nM of each primer, 400 nM deoxynucleoside triphosphates, 0.5 unit SureStart Taq polymerase (Stratagene), 1X SureStart reaction buffer in a final volume of 25 µl.
PCR cycling conditions were: 1 cycle at 94°C for 12 min; 36 cycles at 94°C for 30 seconds, 52.5°C for 30 seconds, and 72°C for 1 min; followed by a 4-min extension step at 72°C. To confirm PCR amplicon identity, the AAV nested capsid PCR (or Ad PCR products-see below) were resolved on 0.8% agarose gels and in-gel Southern blot hybridization was performed using AAV2 cap or adenovirus hexon sequences as the [
-32P]dCTP-labeled hybridization probes. DNA samples identified as containing AAV sequences were subjected to further PCR to isolate the complete AAV capsid coding region.
Due to the length of the cap gene product, a dual PCR approach was employed whereby the 5' half of the cap gene was amplified as a 1.8-kb PCR product, while the 3' half was amplified as a 1.5-kb PCR product. Two different forward 1.8-kb primers were alternatively utilized to amplify these PCR fragments depending on which primer yielded the greatest amplification efficiency (AAV2-1.8F1: 5'-AACATGTGCGCCGTGATTGACGGG-3' or AAV2-1.8F2, 5'-GACCGGATGTTCAAATTTGAACTC-3'). Similarly, 2 different reverse 1.5-kb primers were alternately employed for amplification (AAVCap3'Rev, 5'-TCGTTTCAGTTGAACTTTGGTCTCTGCG-3' or AAVCap3'RevDeg: 5'-CARWRTTYWACTGAMACGAAT-3'). PCR conditions and primer concentrations were the same as for the 255-bp conserved capsid region. Amplified PCR products were agarose gel purified and cloned into pCR2.1-TOPO vector using TOPO TA cloning kit according to manufacturer's instructions (Invitrogen, Inc.). DNA clones were sequenced using BigDye terminator chemistry and an ABI 727 capillary electrophoresis automatic sequencer (PE Applied BioSystems Inc.) by the Columbus Children's Research Institute Core Sequencing Laboratory.
Adenovirus hexon PCR. For detection of adenovirus sequences, total cellular DNA (100 ng) was subjected to identical PCR conditions as that used for AAV capsid PCR. A primer set targeting a conserved region (300 bp) of the hexon gene was employed (9): hex1, 5'-GCCSCARTGGKCWTACATGCACATC-3', and, hex2, 5'-CAGCACSCCICGRATGTCAAA-3'.
Quantitative PCR. AAV genome copy number in tissue samples was quantified using real-time TaqMan PCR analysis (ABI 7700, PE Applied BioSystems). The primers and probe set were selected following alignment of 255-bp cap DNA sequences (ForCAPSS: 5'-AACGACAACCACTACTTTGGC-3' (50 nM); RevCAPSS: 5'-AAGTGGCAGTGGAATCTGTTG-3' (900 nM); probe, [6-FAM]5'-CTACAGCACCCCCTGGGGGTATTTTGA-3') [6-carboxyfluorescein (FAM)-tetramethylrhodamine] (270 nM). PCR conditions were: 50°C 2 min, 95°C 10 min, 40 cycles of 95°C 15 seconds, and 60°C 1 min using 250 ng of human total cellular DNA in 1X Taqman PCR master mix.
Sequence analyses. The DNA and putative protein sequences were aligned and analyzed using the Clustal W method implemented in MegAlgn software in DNASTAR (DNASTAR, Inc). The phylogenetic relationship of all AAV DNA sequences and corresponding putative protein sequences were carried out using Neighbor-Joining method with Kimura two-parameter model implemented in MEGA2 package (19). Similarity percentages between AAV2 and the new AAV sequences were determined using one-pair alignment according to the Lipman-Pearson method. Recombination analysis was performed by using the Similarity Plot method as implemented in the SimPlot software (available at http://www.med.jhu.edu/deptmed/sray/) (21, 22).
Heparin binding analysis. To assess the ability of different AAV capsids to bind heparin, preparations of infectious AAV (harboring the capsid of interest) were subjected to iodixanol density gradient purification (33) and then applied to a POROS HE-20 heparin column (1.7 ml bed volume) using a Biocad Sprint high-pressure liquid chromatography apparatus as previously described (7). Virus was eluted using a linear salt gradient (0.1 to 1 M NaCl). Flowthrough, wash, and eluate (1 ml fractions) were collected for AAV DNase-resistant particle determination.
AAV serology. Subjects with the diagnosis of cystic fibrosis were recruited from Columbus Children's Hospital's Cystic Fibrosis Clinic and were eligible for the study if blood was to be drawn or a centrally placed catheter accessed for a clinical indication. Consent for participation was obtained from subjects (or their legal guardian if they were under 18 years of age) after protocol approval from the Columbus Children's Hospital Institutional Review Board.
To detect antibodies reactive with AAV2 capsid proteins, an enzyme-linked immunosorbent assay (ELISA) was performed. Polystyrene 96-well plates (Nunc Immuno Plate/Maxisorp Surface, Nalge Nunc International) were coated with 50 ng of viral capsid protein. AAV2 capsids were affinity purified as previously described (7) from lysates of 293 cells infected with a recombinant adenovirus type 5 carrying the AAV type 2 capsid coding sequences driven by a standard human cytomegalovirus promoter. After plates were blocked (1% normal sheep serum in phosphate-buffered saline containing 5% dry skim milk) for 2 h at room temperature, serum samples diluted 1:100 in phosphate-buffered saline/5% dry skim milk were added to the test wells and incubated for 2 h at room temperature. After washing 3 times (phosphate-buffered saline containing 0.05% Tween), horseradish peroxidase-linked sheep anti-human immunoglobulin G (Sigma Chemical Co., St. Louis, MO) was added, and the plates incubated for an additional 2 h at room temperature. After three washes, o-phenylenediamine dihydrochloride (Sigma) was added to each well and developed in the dark. Optical density at 450 nm (OD450) was recorded after 30 min (HTS 7000 Bio Assay Reader; Perkin-Elmer Corp.). The OD450 reported was the difference between capsid and mock-coated wells. To show that the AAV capsid preparation was not contaminated with adenovirus proteins, AAV2 coated wells were also developed with adenovirus type 5 immune rabbit serum, and were shown to be at background levels. The relationship between age and OD value was estimated using the Pearson correlation.
Serum neutralization activity directed toward AAV type 2 was determined as previously described (20). Briefly, dilutions of sera were incubated at 57°C for 15 min to inactivate complement and were then mixed for 1 h at 37°C with 1,000 infectious units of a recombinant AAV2 vector that expressed Escherichia coli ß-galactosidase (rAAV/ß-gal). The samples were applied onto monolayers of the C12 cells (8) and incubated at 37°C for 4 h. Subsequently, the monolayers were infected with adenovirus type 5 at a multiplicity of infection of 20. At 48 h after rAAV/ß-gal transduction, cell monolayers were stained for ß-galactosidase activity using the substrate 5-bromo-4-chloro-3-indolyl-beta-D-galactopyranoside (X-Gal). All serum samples were screened at a dilution of 1:10, and those found to have >90% neutralization (reduction in the number of infectious units compared to controls) were assayed again over a range of twofold dilutions (1:100 to 1:1,600).
GenBank accession numbers. The nucleotide sequences described in this paper have been submitted to GenBank. The accession numbers are AY695370, AY695371, AY695372, AY695373, AY695374, AY695375, AY695376, AY695377, and AY695378.
|
|
|---|
Using standard virus culture techniques, including blind passage of cellular lysates after 2 weeks in culture, we were unable to isolate replicating adenovirus from any of the 101 tonsil and adenoid samples. Because replicating adenoviruses were not isolated, we made no further attempt to identify AAV in these cultures. Instead, we prepared total cellular DNA from the 101 samples and then screened for the presence of AAV DNA sequences by PCR. Using a previously identified region within the AAV cap gene as our target (13), we designed and validated a nested, degenerate primer combination that readily detected the published AAV serotypes (Fig. 1A). Control reactions were performed on human cellular DNA that was spiked with plasmid DNA representing cap genes of AAV serotypes 1 to 5. Using the optimized primer sets, we achieved a sensitivity of 15 copies in a background of 100 ng of total cellular DNA (data not shown).
![]() View larger version (36K): [in a new window] |
FIG. 1. PCR schematic for amplification of the complete AAV capsid coding region. (A) The diagram depicts the relative location of degenerate primers (given in Materials and Methods) used to amplify the AAV cap gene. Initially samples were screened with degenerate nested primers (Cap18S, Cap19S, CapSS3189, CapSS2978) to two conserved regions that flank the HVR3 coding region (gray box). To amplify the complete capsid gene, another set of nested primers were constructed (AAV2-1.8F1, AAV2-1.8F2, AAVCap3'Rev, AAVCap3'RevDeg) that bind to 3' regions of rep and cap and amplify 1.8- and 1.5-kb DNA amplicons. (B) Representative amplification of the 255-bp conserved AAV sequence from human tissue DNA (100 ng) following nested PCR (see Materials and Methods for reaction conditions). Asterisks indicate the samples that are positive for AAV amplification.
|
|
View this table: [in a new window] |
TABLE 1. Summary of AAV and adenovirus sequence detection in human tissuesa
|
|
View this table: [in a new window] |
TABLE 2. AAV sequence relatedness and DNA copy number in pediatric tissues
|
AAV DNA in other tissues. To extend the findings described above, we obtained 74 frozen, archived normal tissues from individuals aged 0 to 30 years. The range of tissues and number of samples analyzed are shown in Table 1. Subjects under 0.5 years of age (n = 49) were analyzed but were considered unlikely to have been exposed to adenovirus and AAV. In addition, for three of the 74 samples, the age of the subjects could not be determined. Thus, there were 22 samples from individuals aged 0.5 to 30 years available for analysis, and 16 of the 22 were from individuals between the ages of 0.5 and 14 years. Using the same PCR scheme described for tonsils and adenoids, we found AAV DNA in 2 of the 16 (12.5%) samples from children aged 0.5 to 14 years. One was lung from a 1-year old subject and the other was spleen from an 8-year-old subject. DNA copy numbers were lower in these two samples than in the tonsils and adenoids (Table 2). None of the other 72 samples contained detectable AAV DNA and none of the 74 samples contained adenovirus DNA (Table 1).
AAV cap and rep gene sequences from human tissues.
DNA sequence analysis of the 255-bp amplified cap gene sequence from the nine positive samples revealed significant homology with the corresponding region of AAV2. To isolate complete AAV capsid genes, we synthesized additional degenerate primers at the 3' ends of the rep and cap genes (Fig. 1A). These primers were combined with the forward and reverse conserved sequence primers to amplify 1.8-kb and 1.5-kb PCR products representing the 5' and 3' halves of the cap gene. Individual PCR products were cloned and sequenced; a minimum of 4 individual clones were analyzed for each 1.8-kb and 1.5-kb PCR product (intraclone variation was
0.1%). After a single contiguous sequence was assembled for each tissue isolate, nucleotide sequence alignments revealed significant identity with AAV2 (96 to 97%) for eight sequences. The spleen sequence (S17) was the most divergent and intermediate between serotypes 2 and 3 (Table 2). The complete rep coding sequences were subsequently determined for six of the nine isolates (T32, T40, T70, T71, T88, and S17) and all shared >99% amino acid identity with AAV2. Three conservative mutations in Rep (for all six isolates) were identified (T183A, V508A, and F619S).
The capsid gene amino acid translation and alignment for all nine clones is shown in Fig. 2. Consistent with the nucleotide sequence analysis, eight of the amino acid sequences shared 98% identity with AAV2. Moreover, the majority of the observed amino acid substitutions (relative to AAV2) found in the seven of the tonsil sequences and the lung sequence were conserved among the individual samples. This suggests that a specific virus isolate was circulating in the local population during the time period (winter 2002 to 2003) of tissue procurement. The majority of the observed amino acid substitutions were located in previously identified hypervariable regions (HVRs) 5 to 7, 9, and 10 (11), all of which were predicted to be exposed on the surface of the virion. Two isolates possessed identical sequences (T41 and T71), and two others (T17 and T32) were nearly identical (two amino acid differences). The deduced phylogenetic relationship among the nine cap gene sequences is depicted in Fig. 3, along with previously identified clade B and C viruses (11-13). As expected, seven of the eight sequences clustered closely with each other within AAV2-like clade B.
![]() View larger version (102K): [in a new window] |
FIG. 2. Predicted VP1 capsid amino acid alignment of AAV2 and novel human AAVs. Diagram shows sequence alignment using the CLUSTAL W program. Black boxes designate amino acid substitutions compared to the AAV2 sequence. The locations of previously identified HVR regions (11) are labeled (HVR 1 to 12), as is an additional region (HVR 2') that possesses several substitutions. Several HRV regions (5 to 7, 9, and 10) are colored to facilitate visualization of these regions onto the known atomic structure of AAV2, while invariant HVRs are labeled with black boxes. The locations of R585S and R588T are starred, and arrows denote the approximate locations of nested primers used to amplify the 255-bp HVR3 fragment.
|
![]() View larger version (17K): [in a new window] |
FIG. 3. Phylogenetic analysis of VP1 capsid nucleotide sequences. A neighbor-joining program with a Kimura two-parameter setting was used to derive phylogenetic distances based on 2,200 bp of VP1 sequence. Recently described AAV clade nomenclature (12) was adopted and organized by vertical brackets. The human isolates identified herein are designated in teal type. Due to space restrictions, only a few representative isolates from clades A, D, and E are shown. Sequence isolates are labeled with reference to the source species (bb, baboon; ch, chimpanzee; cy, cynomolgus macaque; hu, human; rh, rhesus macaque). Clade B sequences possessing R585 and R588 amino acids and predicted to bind HSPG efficiently are labeled in red type. The scale for genetic distance is indicated in the bottom left corner.
|
![]() View larger version (16K): [in a new window] |
FIG. 4. S17 sequence homology comparison with AAV2 and AAV3. Simplot analyses of similarity percentages of S17 VP1 versus AAV 2 (red) and AAV3 (blue) are shown. Data were plotted within a sliding window of 200 bp, centered on the position plotted, with a step size between data points of 20 bp. Positions containing gaps were excluded from the comparison. The bar on the top shows the predicted composition of the S17 capsid gene. The corresponding positions of the HVRs are labeled as magenta boxes (HVR 2' in gray).
|
Infectious AAV preparations representing T70 and T88 were generated from full-length molecular clones derived directly from the tonsil and adenoid tissue. A full description of these and other infectious clones derived from tissues is presented elsewhere (Schnepp et al., in preparation). DNA sequences of the clones were identical to the original cap gene nucleotide sequences generated directly from tissues. When applied to a standard heparin chromatography column (HE20-POROS) under low-salt conditions, 93.5% of the total DNase-resistant particles applied to the column were found in the flowthrough or wash buffers (Table 3). In contrast, prototype AAV2 readily bound to the column, with on average 86% of the DNase-resistant particles eluting at 300 mM NaCl (7). Following chromatography, the identity of each virus isolate was further confirmed by amplifying and sequencing a 600-bp portion of viral DNA that was collected from either the flowthrough or peak fraction. In each case, the predicted sequence was recovered (data not shown).
|
View this table: [in a new window] |
TABLE 3. Heparin sulfate column binding of AAV2 virus preparations
|
![]() View larger version (21K): [in a new window] |
FIG. 5. Surface diagrams of AAV2 trimer atomic models. (A) Electrostatic surface potential of the VP3 AAV2 trimer viewed down the threefold axis (yellow triangle) calculated with GRASP (23) running from negative (red) to positive (blue). Labeled arrows indicate the positions of residues implicated in HSPG binding. (B) Predicted electrostatic surface potential of AAV2 VP3 trimer as a result of R585S and R588T substitutions. Amino acid substitutions were modeled using energy minimization simulations with Quanta (Accelrys, San Diego, CA) prior to generating the electrostatic potential map in GRASP. The surface electrostatic potential scale is the same as depicted in panel A. Highlighted regions denote predicted HSPG coreceptor engagement domains in the VP3 trimer.
|
Sequence variation in the AAV capsid maps to surface exposed regions. Given the recently described atomic ribbon structure of AAV2 (25), we were able to map the observed amino acid substitutions in the 8 AAV2-like sequences onto the capsid surface (Fig. 6A and 6B). The majority of amino acid substitutions were located in areas of the capsid predicted to be surface exposed and previously identified as hypervariable regions (HVRs 5 to 7, 9, and 10) (Fig. 3). In contrast, regions predicted to encode core ß-barrel domains responsible for structural integrity were almost invariant. Similarly, the more divergent S17 capsid amino acid sequence was modeled onto the same atomic structure (Fig. 6C and 6D). The majority of the S17 amino acid substitutions also mapped to surface exposed HVR regions 5 to 7, 9, and 10. Interestingly, several HVR 10 substitutions (purple shaded region) were located near the central pore complex that is thought to directly interact with the viral encapsidation complex.
![]() View larger version (53K): [in a new window] |
FIG. 6. Ribbon diagrams of atomic models of AAV2 VP3 trimers showing the location of predicted amino acid substitutions in the human AAV isolates. (A) Ribbon drawing viewed down the threefold axis of symmetry of the AAV2 VP3 trimer. C backbones for the three VP3 monomers are rendered as teal ribbons. Predicted locations of the observed amino acid substitutions present within the eight AAV2-like sequences are color coded to reflect HVR location (HVR 5 to 7, 9, and 10) within the primary sequence (Fig. 2). White space-filling amino acid substitutions mapped outside the known HVRs. (B) Side view of the predicted location of the observed amino acid substitution demonstrating surface display (right side). (C) Superimposition of observed S17 amino acid substitutions relative to the AAV2 VP3 trimer atomic structure viewed down the threefold axis. (D) Side view of the predicted location of the observed amino acid substitutions in isolate S17 (surface display oriented on right side). Images were generated in NAMD/VMD (UIUC Theoretical Biophysics Group) and rendered using Raster3D.
|
|
View this table: [in a new window] |
TABLE 4. Summary of synonymous and nonsynonymous nucleotide substitutions in rep and cap
|
Seroreactivity to AAV2. To further characterize AAV infection in children, we would have liked to evaluate the AAV2 serologic status of our 175 subjects. Unfortunately, we were unable to collect or acquire serum from these individuals. Instead, we screened a separate cohort of 68 individuals (ages 3 to 39) drawn from a local cystic fibrosis clinic. Sera were assayed both for binding antibodies (ELISA) to the AAV2 capsid and for the ability to neutralize AAV2 in vitro (Fig. 7).
![]() View larger version (13K): [in a new window] |
FIG. 7. Serology of AAV infection as a function of age. OD450 values from a standard ELISA (see Materials and Methods) are plotted versus the age of the subject. Sera were tested at a 1:100 dilution. OD values below 0.2 (thin solid line) were considered negative. The same sera were tested for neutralization activity against AAV2 (see Materials and Methods). Data points that are circled represent samples that had neutralization titers of >1:100.
|
After the age of 14 years, the number of seropositive samples rose dramatically to 76% (32 of 42). The correlation with age was significant (r = 0.464, P < 0.001 by Pearson correlation). Only 24% (10 of 42) of the sera from those older than 14 years possessed significant in vitro neutralizing activity. Considering all ages, only 16% (11 of 68) of the samples had significant (titer > 1:100) in vitro neutralizing activity. There were six samples from those over 14 years that had significant ELISA activity (OD > 0.4) that did not mediate significant in vitro neutralization. Such ELISA reactivity probably represented antibodies that bound to the capsid but did not neutralize the virus.
|
|
|---|
Molecular epidemiology of AAV infection in children. Our discovery of AAV2 DNA in tonsils and adenoids from children was not unexpected. AAV2 appears to be the most prevalent of the human AAV serotypes (3, 4, 17, 31), and importantly, our cohort was temporally (winter) and geographically (central Ohio) restricted. Moreover, the portal of entry for AAV would be expected to parallel that of adenovirus. Thus, the oropharynx would be a fertile field from which to harvest AAV (and adenovirus).
While the DNA sequences recovered from these children were similar to the prototype AAV2 and to each other, two significant observations emerged from the sequence analyses. First, the majority of the observed amino acid substitutions were located in previously defined hypervariable regions (11) predicted to be exposed on the surface of the virion (Fig. 6A and 6B). This observation also extended to the more distantly related S17 (AAV2/3) capsid sequence (Fig. 6C and 6D). These data, considered together with the observed nucleotide substitution patterns, suggest that the HVR domains are structurally flexible and possess the capacity to evolve, perhaps in response to host immune pressure. The fact that rep gene sequences from the same tissues were nearly invariant further supported the idea of functional constraints. The identification of the S17 AAV2/3 chimeric sequence confirmed earlier published observations (12), and suggested another mechanism (intermolecular recombination) involved in AAV evolution.
Second, and perhaps most intriguing, was the discovery that none of the seven AAV2-like sequences derived from tonsils and adenoids were predicted to bind HSPG. This was also true for the AAV2-like sequence recovered from lung tissue. To confirm that other compensatory mutations that would restore HSPG binding in the capsid had not occurred, we generated infectious molecular clones from 2 of the tissue samples and demonstrated they indeed did not bind HSPG (Table 3). To extend this observation, we examined other available AAV2-like sequences (12) and noted that 77% of the 31 currently identified clade B AAV2-like sequences lack R585 and R588. Thus, these data suggest that preponderance of AAV2-like isolates do not bind HSPG, and that HSPG is not required for natural infection in humans. This notion is consistent with the fact that most other serotypes of AAV do not bind HSPG.
Seroepidemiology of AAV2 infection in children. To see if the rates of AAV infection as judged by DNA isolation (7 to 12%) were consistent with serological estimates, we analyzed a set of sera from children in the same geographic locale. Ideally, we would have also analyzed sera from the cohort of children who donated tissue samples, but we were unable to obtain serum from these children. In those children 14 years and younger (same ages as our tonsil and adenoid cohort), the seropositivity rate was roughly 12% and thereby confirmed our estimates derived from AAV DNA detection in tissues.
AAV disseminates beyond the portal of entry. Previously published work has described the presence of AAV DNA in multiple adult human tissues (12). We have now extended these findings to children with the discovery of AAV DNA in tissues within and beyond the probable portal of entry (orophaynx). In children ages 0.5 to 14 years, we found AAV DNA in 2 (lung and spleen) of the 16 nonoropharyngeal tissue samples (12.5%) available for analysis. Although the number of samples analyzed was smaller, the percentage containing AAV DNA compared favorably (12.5% versus 18%) with data from adults (12).
To extend beyond the portal entry, infectious agents gene-rally use one of 3 pathways: direct (contiguous) spread, lymphatic, or bloodstream. In our samples, AAV could have easily reached the lung by contiguous spread from the oropharynx, either with or without adenovirus. To reach the spleen, however, the route of viral spread was almost certainly hematogenous, in the form of free or cell-associated (perhaps leukocytes) virus.
Biology of wild-type AAV infection. The now emerging picture is that AAV infection of humans is more complex than previously appreciated. It has been known for decades that AAV is not associated with disease or pathology, and that infection incidence in humans generally parallels that of adenovirus. Not surprisingly, infectious AAV has been isolated from sites where adenovirus is traditionally recovered, including the gastrointestinal tract (1). More recently, it has been appreciated that AAV DNA can be found in many human tissues, including liver, muscle, lymph nodes, leukocytes, kidney, and cervical tissues (11, 12, 15, 16, 28, 29), and now tonsils and adenoids.
Considered together, these data suggest the following biologic scenario. AAV most likely enters the body through the oropharynx in association with adenovirus. Replication ensues (with adenovirus help) and new AAV particles are formed and released from infected cells in the oropharynx. Secondary rounds of replication in newly infected cells follow, again creating new waves of AAV particles. Such rounds of replication would continue until the host immune system responds and blunts the infectious process. By the time replication is controlled, AAV has had the opportunity to spread to the lungs (contiguous) and through the bloodstream to more distant sites. This scenario does not address the rare AAV5 serotype, which apparently enters the body through the genital tract (1, 14). However, there remains only a single isolate of AAV5, and recent studies of human tissues have failed to find AAV5-like DNA.
While the events envisioned above are entirely plausible, many important questions regarding the in vivo biology of natural AAV infection remain unanswered. For example, while AAV and adenovirus appear to be linked early in the infectious process, there appears to be an unlinking sometime during and following dissemination, allowing AAV DNA (but not adenovirus) to persist in organs and tissues outside the oropharynx. It is formally possible that AAV DNA persists by integrating in target cells on chromosome 19 (AAVS1), but this has not been demonstrated in vivo. Moreover, the target cells for AAV replication and persistence have not been identified, nor have the specific cellular receptors for viral attachment been defined.
It should be remembered that even with widespread dissemination, AAV has not been found to cause disease or pathology. Nonetheless, the effect of prior (or subsequent) wild-type AAV infection on gene transfer mediated by recombinant AAV vectors is unknown, and a more thorough understanding of the natural AAV infectious process is needed.
This work was funded by National Institutes of Health (NIAID/DAIDS) grant 2-PO1-AI56354.
Present address: Children's Hospital of Philadelphia, Philadelphia, PA 19104. ![]()
C.-L.C. and R.L.J. contributed equally to this work. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»