Previous Article | Next Article ![]()
Journal of Virology, March 2003, p. 3495-3504, Vol. 77, No. 6
0022-538X/03/$08.00+0 DOI: 10.1128/JVI.77.6.3495-3504.2003
Copyright © 2003, American Society for Microbiology. All Rights Reserved.
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 4 October 2002/ Accepted 18 December 2002
|
|
|---|
|
|
|---|
Wild-type (wt) AAV integrates primarily (
60%) in a site-specific fashion into human chromosome 19 in cultured cells (27, 29, 45). This site-specific integration appears to be mediated by virally encoded rep gene products through the recognition and binding of specific viral and cellular consensus sequences (18, 56). Importantly, almost nothing is known about wt AAV integration or persistence in vivo. Approximately 80% of adults have antibodies reactive with AAV-2 (5), but it is not known whether latent provirus persists (or in what form) within these individuals (39). To date, only a single study has examined wt AAV infection in primates (22). In this study, nine rhesus macaques were inoculated with wt AAV in the presence or absence of wt adenovirus. Using PCR, AAV DNA was found in several animals, while Southern blot analysis of DNA taken from sites of inoculation showed no positive signals. Site-specific integration (measured by PCR amplification of wt AAV-cellular DNA junctions) was detected in one animal (a very faint signal was seen in two others). Thus, it can be concluded that even when conditions are permissive for in vivo site-specific integration, this unusual biologic phenomenon is not a common event (22).
It is important that most rAAV-2 vectors lack the AAV rep gene and thus are not considered to have the potential to undergo targeted integration. Nonetheless, rAAV vectors have been shown to integrate in vitro in many cell types. However, the vast majority of these studies were performed with vectors encoding drug resistance markers that demonstrated integration of rAAV vector in drug-resistant clones (21, 30, 36, 38, 43, 50, 55). In studies in which rAAV vector persistence has been analyzed in the absence of selection in vitro, it appears that rAAV vectors integrated at a very low frequency. In fact, transformed cells in culture transduced with rAAV vectors (in the absence of selection) have shown a loss of transgene expression over time (2, 34). The loss of transgene expression corresponded to a loss of rAAV vector DNA, suggesting that rAAV vectors persisted predominately as episomes that declined in a dividing cell population. These observations stand in stark contrast with the stable gene expression found following retroviral vector transduction (2, 34).
A somewhat analogous situation appears to exist for plasmid DNA transfected into transformed cells in culture. There is a loss of plasmid DNA copies in cells over time due to the episomal nature of this DNA (14). However, under selective conditions, plasmid DNA has been shown to readily integrate in vitro, predominately in head-to-tail arrays (58). Conversely, several in vivo studies have demonstrated that plasmid DNA persists only at very low levels following intramuscular injection in mouse tissue (31, 35).
Although the basis for rAAV genome persistence remains unclear, three possible events seem plausible: (i) nonhomologous recombination-integration into cellular DNA; (ii) maintenance of extrachromosomal forms; or (iii) a combination of integrated and episomal forms. Several studies have suggested, on the basis of hybridization patterns observed following Southern blot analysis of total cellular DNA from skeletal muscle, liver, and brain, that rAAV vector DNA might be integrated in vivo (7, 8, 15, 33, 59). In these studies, rAAV vector DNA was found associated with high-molecular-weight DNA, although no formal proof of integration was offered (i.e., isolation of vector-cellular junctions or proviruses). Shortly thereafter, based on fluorescent in situ hybridization data, rAAV vector sequences were reported to integrate in the liver (37); however, the interphase nature of the published chromosomal spreads made interpretation of data difficult. Also, rAAV integration was observed in a genetic selection model that identified rare insertional events in which a small population of hepatocytes transduced with a therapeutic rAAV vector were able to repopulate a defective liver (4). In another study, vector-cellular DNA junctions from mouse liver DNA were isolated by PCR (40). However, the frequency of these integrative events within the liver appeared to be very low, as evidenced by the complete loss of detectable rAAV vector in regenerating liver tissue following portal vein infusion (42).
In contrast, episomal forms of vector DNA have been easier to document and characterize in vivo. We previously observed a low-molecular-weight rAAV species that was characterized as a double-stranded episome (8). Additionally, other investigators using Southern blot hybridization have observed monomeric and dimeric rAAV vector forms in different tissues (1, 41, 42, 47, 48, 51). More recently, monomeric and multimeric circularized rAAV vector forms have been directly isolated from mouse skeletal muscle (13).
Considered together, these data do not allow us to clearly define the genetic fate of rAAV vector genomes in transduced tissues from animals. On the one hand, there appears to be evidence that rAAV can integrate, but it is likely a very rare event that remains unquantified. On the other hand, episomal forms of vector DNA appear to be commonly found in transduced tissues, although no attempt has been made to assign a relative frequency to episome formation.
Because of this uncertainty, we sought to quantitatively define the proportion of total rAAV vector genomes that persist in one form or the other in skeletal muscle tissue. To accomplish this goal, we used three distinct and complementary analytical approaches to examine the genetic fate of rAAV vector genomes in transduced muscle from mice. These analyses allowed us to conclude that rAAV vector genomes in muscle persist primarily (>99.5%) as large and small episomal concatemers. We were unable to demonstrate integration events in muscle even by using a sensitive PCR-based assay.
|
|
|---|
Cells and viruses. HeLa, 293, and C2C12 cells were purchased from the American Type Culture Collection (Rockville, Md.). HeLa, 293, and C2C12 cells were maintained in Dulbecco's minimal essential medium supplemented with 10% fetal bovine serum and penicillin-streptomycin. rAAV producer cell lines were maintained in the above-described medium supplemented with 700 µg of G418/ml. Production of rAAV vectors encoding the mouse ß-glucuronidase gene or the chicken ovalbumin gene were generated from producer cell line GUS-g2 or ova-6, respectively (9). Both vectors contain the human cytomegalovirus (HCMV) promoter/enhancer for high-level transgene expression. rAAV was generated by growing 1010 producer cells for 5 days in a Corning Cell Cube and then infecting the cells with adenovirus type 5 at a multiplicity of infection (MOI) of 50. After 72 h, the cells were concentrated by low-speed centrifugation (1,000 x g) and resuspended at a density of 107 cells/ml in 20 mM Tris (pH 8.0)-1 mM MgCl2-150 mM NaCl. Cells were lysed by treatment with 0.5% deoxycholate and Benzonase (35 units/ml). Vector purification was achieved using heparin affinity column chromatography as described previously (6). Using the PE Applied Biosystems Prism 7700 TaqMan sequence detector system, DNase-resistant particle (DRP) titers on both vectors were determined as previously described (6). The HCMV promoter-specific primer-probe sequences used were as follows: forward primer, 5' TGAAATCCCCGTGAGTCAA 3' (300 nM final concentration); reverse primer, 5' CATGGTGATGCGGTTTTGG 3' (200 nM final concentration); and HCMV probe, 5' FAM-CCGCTATCCACGCCCATTGATG-TAMRA 3' (200 nM final concentration). The presence of wt-like AAV was assessed by a double-blind passage assay using 1% of the rAAV stock as previously described (6). Based on assay sensitivity (1 IU), each vector preparation contained <1 IU of wt-like AAV per 1011 rAAV particles (DRP). The presence of replication-competent adenovirus in the rAAV preparations was assayed by passing 1% of the rAAV vector stocks onto 293 cells and scoring for adenovirus cytopathic effect after 7 days of incubation (6). The vector stocks did not possess detectable replication-competent adenovirus (<1 PFU per 1011 rAAV particles). Vector stocks contained <1 endotoxin unit per 50-µl dose (Limulus Ameobecyte Assay; Cape Cod Inc.). Adenovirus type 5 was purified from infected 293 cells by CsCl step and isopycnic gradient centrifugations according to standard protocols (19).
rAAV/eGFP-zeo stable cell lines were generated by transducing C2C12 cells with an rAAV vector encoding an enhanced green fluorescent protein (eGFP)-zeocin fusion gene at an MOI of 104 DRP per cell. After 48 h, the cells were plated into 100-mm-diameter dishes at a density of 1 x 105 cells/dish and selected in the above-described medium supplemented with 800 µg of zeocin/ml. The stable cell line pool was passed 10 times (
40 cell doublings) in selective medium to ensure that extrachromosomal rAAV vector forms were lost via dilution.
Animals. Adult male BALB/c mice were used for all studies described in this paper, and experiments were carried out in accordance with local institutional and National Institutes of Health guidelines. The animals were anesthetized with Avertin (240 mg/kg of body weight, given intraperitoneally by injection intramuscularly into the hind limb). A skin incision (0.5 cm in length) was made over the distal femur, and using a 28-gauge needle, 2.5 x 1011 DRP of rAAV/gus or 1.0 x 1011 DRP of rAAV/ova (50-µl volume) was injected along the longitudinal axis of the quadriceps muscle. Control animals received PBS in an identical manner. The entire quadriceps muscle group was removed and rapidly frozen in liquid nitrogen prior to DNA extraction.
Phage library construction and screening. Total cellular DNA was isolated from mouse skeletal muscle at 12 weeks postinoculation with rAAV/ova (1011 DRP). TaqMan PCR revealed that this muscle DNA sample possessed 36 vector copies/nucleus. Using a QIAEX II extraction kit (Qiagen, Inc.) according to the manufacturer's specifications, 4 µg of DNA was digested with EcoRI and XhoI (XhoI cuts once within the rAAV vector sequence) and separated on a 0.8% agarose gel. DNA ranging in size from 3 to 12 kb was extracted from the gel. The recovered fragments were ligated, using a Uni-ZAP XR vector system (Stratagene), into lambda phage vector arms (predigested with XhoI and EcoRI) and packaged (using a Uni-ZAP XR Gigapack cloning kit according to the manufacturer's instructions) into recombinant phage. An analysis of 50 rescued phagemids revealed that the phage library was 98% recombinant, with an average insert size of 4 kb. Using standard methods (44), a total of 1.5 x 106 phages were screened with an HCMV-specific radiolabeled probe.
B1 PCR.
Total cellular DNA (20 to 40 ng) isolated from mouse skeletal muscle or C2C12 cells was used in all B1 PCRs; this amount of DNA gave the greatest level and size distribution of B1 amplicon amplification. Higher sample DNA amounts appeared to rapidly exhaust the primer and nucleotide pools, resulting in PCR products with a smaller size distribution. The DNA was first digested with PacI (to reduce sample viscosity) and then amplified in 20 mM Tris-HCl (pH 8.75)-10 mM KCl-10 mM (NH4)2SO4-2 mM MgCl2-400 µM each deoxynucleoside triphosphate-5 U of Herculase Hotstart DNA polymerase (Stratagene)-2 µg of B1 primer in a 50-µl reaction mixture. The B1-repeat primer sequence used was 5' AGTTCCAGGACAGCGAGGGCTAYACAGA 3'. The reaction conditions were as follows: denaturation at 95°C for 2 min, followed by 10 cycles of 95°C for 30 s, 64°C for 30 s, and 72°C for 7 min 30 s, with an additional 20 cycles of 95°C for 30 s, 64°C for 30 s, and 72°C for 10 min 30 s. A portion (2.5 µl) of the B1 PCR mixture was then used as the template for real-time PCR (TaqMan) to quantitate vector sequence-specific amplification in the B1 PCR. The TaqMan PCR conditions used were as follows: 50°C for 2 min and 95°C for 10 min, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min in a 25-µl reaction volume. For sample quantification, a plasmid DNA standard curve was generated in triplicate using of 101 to 106 plasmid copies; experimental values were obtained by extrapolation (using system software) from the plasmid DNA standard curve (coefficient of linearity,
0.995).
To determine the initial rAAV vector copy number in the muscle samples, 25 ng of total muscle DNA was used in a standard TaqMan PCR with the HCMV primer-probe combination. The resulting vector copy number was divided by 4,167 nuclei (assuming 6 pg per nuclei) to obtain the number of vector copies per nucleus. The ovalbumin-specific TaqMan PCR primers and probe were as follows: forward, 5' CTGACCATCCATTCCTCTTCTGT 3' (300 nM final concentration); reverse, 5' GAAACACATCTGCCAAAGAAGAGA 3' (300 nM final concentration); and probe, 5' FAM-TCAAGCACATCGCAACCAACGCC-TAMRA-3' (200 nM final concentration). To analyze the integrity of genomic DNA, the following mouse gus-specific primers were used: forward, 5' ACCCTGCGGTTGTGATGTG 3' (50 nM final concentration); reverse, 5' AATATGCGG CGGGTTTCAG 3' (300 nM final concentration); and probe, 5' FAM-TGTGGCCAATGAGCCTTCCTCTGC-TAMRA 3' (200 nM final concentration).
Cloning and expression analysis of rAAV in vivo vector forms. One microgram of total muscle DNA, isolated from an rAAV/ova-injected mouse, was digested with EcoRI (an enzyme which does not cut within the rAAV vector) and treated with Plasmid-Safe DNase as described above. This DNA was used as a template in a PCR mixture containing 20 mM Tris (pH 8.4), 50 mM potassium chloride, 1.5 mM magnesium chloride, 1 mM deoxynucleoside triphosphates, and 2.5 U of Platinum Taq DNA polymerase (Gibco BRL) in the following reaction conditions: 94°C for 2 min for 1 cycle and 94°C for 30 s, 55°C for 30 s, and 72°C for 5 min for 30 cycles. The forward primer used was 5' AACAGCTGCAGATCAAGCCAGAGAGCTCAT 3', and the reverse primer was 5' TTGATCTGCAGCTGTTTGAAAGTTGATAGG 3'. The PCR product was cloned into a TOPO-TA cloning vector (Invitrogen). Positive colonies were identified using a radioactive HCMV probe by a standard colony hybridization method (44). rAAV vector fragments were released by PstI digestion and self-ligated using T4 DNA ligase (New England Biolabs). Using Superfect Transfection reagent (Qiagen), the ligation reaction mixture was directly transfected into 3 x 105 HeLa cells. A cell lysate was generated after 48 h, electrophoresed on a SDS-10% polyacrylamide electrophoresis gel, and subsequently transferred to a polyvinylidene difluoride membrane (Amersham). The membrane was blocked for 10 min at 25°C using TBST (10 mM Tris [pH 8.0], 150 mM NaCl, 0.1% Tween 20) plus 5% nonfat dry milk. The primary antibody was a monoclonal anti-chicken egg albumin (Accurate Chemical) and was hybridized in blocking solution at 1:1,000 for 1 h. The blot was washed three times for 10 min in TBST and incubated for 30 min in TBST plus 2.5% dry milk with a horseradish peroxidase-conjugated anti-mouse secondary antibody (Vector Laboratories) at a 1:10,000 dilution. The blot was washed as above and developed using an ECL Plus detection system (Amersham).
|
|
|---|
PCR-based assay for rAAV integration. To complement and extend the direct cloning approach described above, we developed a PCR-based assay that was designed to detect integrated rAAV vector DNA in mouse tissue. This approach was predicated on the random amplification of mouse genomic DNA that is flanked by Alu-like mouse interspersed repetitive sequences (known as B1 elements) by using a B1-specific PCR primer. Mouse B1 repetitive elements are approximately 135 bp long and are present at 80,000 to 100,000 copies per haploid genome (25). Thus, B1 elements comprise about 1% of the total mouse genome and exist on average every 15 kb. Analysis of B1 PCR products using mouse genomic DNA as the template and a single B1 primer revealed a smear ranging in size from >17 kb down to 200 bp (data not shown). Thus, the size distribution of our optimized B1 PCR reflected the expected size distribution based on B1 repeat frequency.
In this assay, rAAV vector DNA integrated between two B1 elements is amplified and the amplified products (containing rAAV integrants) are then quantified using TaqMan PCR. The assay readout is expressed as fold amplification of the rAAV vector DNA (which occurs in the B1 PCR) over the carryover input rAAV vector sequences in the subsequent TaqMan PCR. Hence, net sequence amplification (over input) of the rAAV target sequence following the B1 PCR is diagnostic for rAAV integration. Importantly, episomal rAAV vector DNA is not amplified using the B1 PCR strategy and no net amplification of episomal rAAV vector DNA is predicted to occur.
For example, consider the situations illustrated in Fig. 1. In the case in which there are only nonintegrated copies of rAAV vector DNA (Fig. 1, left lane), no amplification occurs (episomes not flanked by B1 sequences) and there is no fold amplification over input. In the case in which there are only integrated forms (Fig. 1, center lane), the input signal (100 integrants) is amplified 1,000-fold in the B1 PCR. In the last case, both episomal and integrated forms are present (number in episomal form > number in integrated form). The same degree of amplification occurs, but the fold amplification is decreased because the episomal forms are detected in the TaqMan PCR quantification. However, this level of amplification still allows integration events to be detected.
![]() View larger version (35K): [in a new window] |
FIG. 1. Schematic of the B1 PCR strategy. Shown in this schematic are three potential forms of rAAV vector genomes in the host cell nucleus. As shown in the left panel, episomal vector forms do not amplify in the B1 PCR because B1 sequence elements are not present. As shown in the middle panel, integrated rAAV vector DNA are expected to amplify in the B1 PCR. Finally, as shown in the right panel, when both forms are present only the integrated proviral vector DNA amplifies in the B1 PCR. As demonstrated here, nonamplified episomal vector DNA carries over into the subsequent TaqMan PCR and is detected as input signal.
|
|
View this table: [in a new window] |
TABLE 1. Characterization of B1 PCR amplification
|
Because it is possible that rAAV genomes in vivo are integrated and episomal, we next tested the ability of the B1 PCR assay to identify integrated rAAV vector DNA within a competing episomal vector DNA background. The effect of competing episomal vector DNA on B1 PCR assay sensitivity can be readily illustrated by comparing the fold amplification required to observe net vector amplification in the absence of competing episomal DNA to that required in the presence of competing episomal DNA. For example, to detect 150 integrated vector copies in the absence of competing episomal DNA, we would simply need to amplify the integrated sequences by a factor greater than 4.8-fold (our previously defined cutoff value) to score this sample as positive for vector integration. In fact, at the 150-copy level, we observed average amplification of 451-fold over input vector (Table 1). On the other hand, to detect 150 integrated vector copies in the presence of 15,000 episomal copies, the integrated sequences would need to be amplified by a minimum of 385-fold to see net vector amplification. This is because if 150 integrated vector copies were amplified 385-fold in the B1 PCR assay, after amplification there would be a total of 57,700 integrated vector copies plus 15,000 unamplified episomal copies (72,700 total copies). If 5% of this reaction mixture (0.05 x 72,700 = 3,636) and 5% of the input carryover (0.05 x 15,150 = 757.5) were quantified by TaqMan PCR, the fold amplification would be 3,636/757.5 = 4.8. In practice, at the 150-copy level, we observed an average amplification of 16.3-fold over input vector (Table 1). Thus, the effect of the presence of high levels of episomal vector target is to reduce the net amplification, because the input carryover (denominator) is increased (Fig. 1).
In light of these observations, a series of episomal and integrated rAAV vector spike experiments were performed to define assay sensitivity in the presence of competing episomal vector DNA (Table 1). Each reaction mixture contained 1.5 x 104 plasmid copies with different percentages of integrated rAAV vector DNA. We tested dilutions ranging from 50% of the integrated target (1.5 x 104 plasmid copies + 1.5 x 104 integrated copies) down to 0.01% of the integrated target (1.5 x 104 plasmid copies + 1.5 x 100 integrated copies). The B1 PCR results indicated that we were able to detect significant amplification of 75 integrated rAAV vector copies within an episomal background of 1.5 x 104 plasmid copies (Table 1). Thus, the assay was able to detect rAAV vector integration in as few as 16.7 cell equivalents (assuming 4.5 copies per cell) in a background of 3,333 cells (assuming 6 pg of DNA per cell), even with competing input episomal DNA at levels exceeding 104 vector copies.
To verify that all genomic templates were amplifiable in the B1 PCR, we used the native mouse gus locus as an internal control to monitor B1 PCR efficiency. Based on the known genomic sequence of the gus gene, a portion of the gus locus is flanked by two B1 elements spaced 1.5 kb apart. Significant levels of gus locus amplification were observed in all genomic DNA samples (>700-fold; data not shown) following B1 PCR assay analysis.
Analysis of skeletal muscle for integrated rAAV vector sequences using B1 PCR. Total cellular DNA was isolated at various times (6 to 57 weeks) postinjection from quadriceps tissue of rAAV/ova-injected animals. The average number of vector copies per nucleus for each muscle sample was determined by quantitative PCR (Table 2). Muscle from the earlier time points contained higher vector copies/nucleus than later time points; this was attributable to the diminution of single-stranded input vector genomes over time. As discussed above, the B1 PCR assay was optimized for the use of 1.5 x 104 vector copies in 20 to 40 ng of genomic DNA. Therefore, mouse DNA samples having more than 1.5 x 104 vector copies in 20 ng of genomic DNA were diluted in total cellular DNA isolated from naïve mice to achieve this final copy number. Mouse DNA samples having fewer than 1.5 x 104 vector copies in 20 ng of genomic DNA (mice J1, J2, and I4) were analyzed using 40 ng of DNA so that a total of 1.5 x 104 vector copies could be analyzed.
|
View this table: [in a new window] |
TABLE 2. B1 PCR of transduced mouse muscle
|
Use of a novel DNA exonuclease to characterize the form of vector DNA. Considered together, the failed attempts to identify integrated forms of rAAV vector DNA directly implied that the high-molecular-weight vector DNA which persisted in muscle over long periods of time existed as episomal concatemers. To investigate this possibility, we developed an assay using a novel DNA exonuclease (Plasmid-Safe, henceforth called PS-DNase) to selectively analyze circular DNA. PS-DNase selectively hydrolyzes linear double-stranded DNA and linear and closed-circular single-stranded DNA but not double-stranded, closed circular supercoiled or nicked DNA. This nuclease is routinely used in the course of plasmid DNA purification to selectively remove contaminating bacterial genomic DNA from circular plasmid DNA preparations.
Given these unique properties, we reasoned that PS-DNase might be useful to discriminate between integrated (susceptible to digestion) and episomal (resistant to digestion) rAAV vector DNA that was associated with high-molecular-weight genomic DNA (Fig. 2). When used in combination with specific restriction enzymes, we predicted that episomal rAAV DNA would be protected from degradation by PS-DNase and that this protection (or degradation) would be evident on Southern blots.
![]() View larger version (43K): [in a new window] |
FIG. 2. Schematic of the PS-DNase assay for episomal rAAV genomes. The predicted outcomes of treating rAAV-transduced muscle DNA with selected restriction enzymes in combination with PS-DNase are shown in this diagram. In the case of treatment with an enzyme that cuts once within the genome (Single Cutter), subsequent treatment with PS-DNase leads to degradation of all rAAV genomes, regardless of their forms (episomal or integrated). In the second case, when the total DNA is cut with an enzyme that does not cut within the rAAV genome (Zero Cutter), after treatment with PS-DNase, integrated forms are degraded while episomal forms are protected.
|
![]() View larger version (40K): [in a new window] |
FIG. 3. Specificity of PS-DNase on plasmid DNA substrates. Southern blot analysis was performed on 1 ng of the indicated plasmid DNA templates with (+) or without (-) PS-DNase treatment. Linear plasmid DNA (lanes 1 and 2) was completely hydrolyzed in the presence of PS-DNase (lane 2). As expected, open circular (lanes 3 and 4) and supercoiled DNA (lanes 5 and 6) forms were resistant to PS-DNase treatment (lanes 4 and 6).
|
![]() View larger version (63K): [in a new window] |
FIG. 4. Specificity of PS-DNase on rAAV replicative forms. Southern blot analyses of rAAV vector forms (1 ng) treated with (+) or without (-) PS-DNase are shown. Recombinant AAV single-stranded genomes were isolated from a highly purified rAAV preparation, and the DNA was denatured by boiling (note the presence of single-stranded genomes of both polarities) (lane 1). A similarly treated DNA sample was allowed to anneal (lane 2) and shows a prominent band corresponding to double-stranded unit length vector, in addition to a smear of complex annealed forms. All reannealed forms were susceptible to PS-DNase treatment (lane 3). The results of an analysis of the susceptibility of cellular rAAV replicative forms (Rf DNA) to PS-DNase digestion are shown in lanes 4 and 5. Monomeric (Rfm), dimeric (Rfd), and higher-order (lane 4) forms were completely hydrolyzed by PS-DNase coincubation (lane 5).
|
Total cellular muscle DNA was isolated from two animals at 10 or 15 weeks after rAAV injection with either rAAV/gus or rAAV/ova, respectively. The resulting Southern blot hybridization patterns from the two animals were nearly identical before and after treatment with the indicated enzymes; data from the 15-week animal are shown in Fig. 5. Uncut genomic DNA contained several rAAV vector DNA hybridizing species, including high-molecular-weight forms, a monomeric vector form, and residual input single-stranded vector DNA. When genomic DNA was digested with a restriction enzyme that cut once within the rAAV vector, the high-molecular-weight and monomeric bands were reduced to a prominent band which (based upon size) corresponded to a unit length head-to-tail fragment (Fig. 5, lane 2). This was consistent with the rAAV vector DNA being organized in concatameric head-to-tail arrays. When this same linearized genomic DNA was exposed to PS-DNase, complete loss of signal was observed on both the Southern blot (Fig. 5, lane 3) and the same gel stained with ethidium bromide (Fig. 5, lane 8).
![]() View larger version (101K): [in a new window] |
FIG. 5. Southern blot analysis of rAAV vector forms in skeletal muscle following PS-DNase treatment. Total cellular DNA was isolated from skeletal muscle at 15 weeks following injection with rAAV/ova. In the left panel, a Southern blot hybridized with an HCMV vector-specific probe is shown. In the right panel, the identical gel stained with ethidium bromide is shown. Lanes 1 and 6, uncut DNA; lanes 2 and 7, PstI-digested DNA (Single cut [in the rAAV genome]); lanes 3 and 8, PstI-digested DNA after PS-DNase digestion; lanes 4 and 9, EcoRI-digested DNA (Zero cut [within the rAAV genome]); lanes 5 and 10, EcoRI-digested DNA after PS-DNase digestion. PS-DNase completely degraded linearized rAAV vector DNA (lanes 3 and 8) but did not degrade EcoRI-digested (Zero cut) rAAV vector DNA (lane 5, arrowheads), even though the cellular genomic DNA was completely degraded (lane 10).
|
Episomal rAAV vector genomes are transcriptionally active. Since the PS-DNase and B1 PCR experiments both indicated that the majority of rAAV vector DNA persists in an episomal form, we were interested in whether these forms were transcriptionally active and contributed to the observed transgene expression. Total cellular DNA was isolated from rAAV/ova-inoculated skeletal muscle at 15 weeks postinjection and digested with a no-cut restriction enzyme (EcoRI) and was then treated with PS-DNase to degrade linear genomic fragments. By performing TaqMan analysis with a gus primer-probe set, we verified that >99.9% of the genomic DNA signal had been eliminated. The PS-DNase-resistant rAAV vector DNA (episomal form) was used as a template for circle PCR, in which PCR primers designed to amplify away from each other were positioned at a unique PstI restriction site. The products of the PCR were cloned into a TA cloning vector, and bacterial colonies were screened using a vector-specific hybridization probe. Positive clones ranged in size from 0.3 to 3.2 kb in length, and 15 clones (>3.0 kb) were chosen for further analysis. To determine whether the rescued rAAV sequences were capable of expressing the ovalbumin transgene, the inserts were released from the cloning vector with PstI digestion and self-ligated to recapitulate a circular episomal form. The ligated rAAV vector DNA was directly transfected into HeLa cells, and ovalbumin expression was confirmed by Western blot analysis in 6 of the 15 clones analyzed (4 are shown in Fig. 6). Sequence analysis of the clones revealed that they all possessed rearranged ITR elements, with various degrees of deletion. Deletion of a portion of the ova expression cassette was observed in three clones. Although the forms of the rearrangements and sequence deletions are similar to those of previously published in vivo ITR (15), we cannot formally rule out the possibility that the observed sequence alterations occurred as a result of the circle PCR process or upon passage through bacteria.
![]() View larger version (35K): [in a new window] |
FIG. 6. Gene expression from in vivo rAAV vector forms. A circle PCR approach was used to clone the predicted in vivo monomeric circular vector form. PCR primers, overlapping the unique PstI site at their 3' ends, were used to amplify rAAV vectors from PS-DNase-treated muscle DNA. The resulting PCR products were cloned into a TA-cloning vector for further analysis. To determine whether the cloned vectors were transcriptionally active, they were excised from the TA vector with PstI, self-ligated, and transfected into HeLa cells. Cell lysates were analyzed by Western blotting with an anti-ova-horseradish peroxidase antibody. Four independently rescued clones (lanes 1 to 4) that express chicken ovalbumin protein are shown. Lane U, untransfected cell lysate; lane C, lysate from cells transfected with a DNA expression plasmid encoding chicken ovalbumin.
|
|
|
|---|
The B1 PCR assay is by its nature inefficient. B1 elements comprise about 1% of the total mouse genome and are present on average every 15 kb (25). The actual distances between B1 elements differ, and shorter B1 amplicons are expected to amplify more efficiently than longer amplicons. Hence, in certain instances large B1 amplicons containing integrated vector sequences might not be detected in the assay. In fact, our data directly support this proposition. We detected an rAAV integration event in 1 out of every 196 rAAV/eGFP-zeo-transduced C2C12 cells assayed. However, it is important to note that the C2C12 rAAV-transduced pool allowed us to control for the inefficiency of rAAV detection and define the limits of sensitivity for the assay. The integrated pool was composed of mouse muscle (C2C12) cells that were transduced with an rAAV vector at a high MOI (104 DRP/cell). This potentially mimicked in vivo transduction conditions and, following drug selection, generated a population of cells that had the rAAV vector genome randomly integrated into the C2C12 muscle cell genome. Thus, by subjecting this pool of DNA to the B1 PCR assay, we experimentally defined the limits of detection on this heterogeneous mixture of integrants (relative to B1 amplicon position and length) following random rAAV vector integration into host genomic DNA.
Importantly, none of the mice that had received rAAV injections into muscle had detectable integrated vector, and considering the limit of detection, this indicated that at least 99.5% of the rAAV vector in skeletal muscle was not integrated. Thus, two independent experimental approaches (direct cloning and B1 PCR) failed to detect integrated forms of rAAV vector DNA.
These negative findings compelled us to consider an alternative hypothesis, namely, that rAAV vector DNA was persisting as extrachromosomal high-molecular-weight forms. To characterize these high-molecular-weight molecules, we used a novel DNase (PS-DNase) that digests linear molecules but not double-stranded episomes. Based on the predicted activity of PS-DNase and fragment analysis after restriction enzyme digestion, the high-molecular-weight moieties were shown to be concatameric episomes. In addition to the large episomes, we also characterized an abundant circular, double-stranded monomeric vector form. The circular nature of these molecules was predicted by their resistance to PS-DNase digestion and confirmed by direct cloning and sequencing from tissue. Sequence analyses revealed structures with deleted sequences within the viral ITR, an observation consistent with previous in vivo and in vitro data showing significant rearrangement of ITR elements following rAAV transduction (15, 32, 38, 40, 43, 57, 60). Other investigators have also reported the isolation of smaller circular rAAV vector forms from skeletal muscle which were demonstrated to be capable of directing transgene expression (13, 61). Yang et al. (61) also reported that rescued rAAV circular intermediates isolated from skeletal muscle at 120 days postinfection were able to express the transgene. The biochemistry of rAAV episome formation remains undefined and is to be the subject of future investigations.
Thus, considered together, these data argue that very little (if any) rAAV vector DNA was integrated in transduced mouse muscle and that the vast majority of the vector sequences persist as large and small concatameric episomes. In retrospect, these findings are not surprising, because the experimental system did not predict integration. The rAAV vectors used in these studies did not contain the AAV rep gene, and the mouse host cell DNA did not contain the AAV S1 locus (45). Furthermore, even when conditions are permissive for in vivo site-specific integration, this unusual biologic phenomenon is not a common event (22). However, we cannot assert that rAAV vectors never integrate. In fact, for many common modalities of in vivo exogenous DNA transfer, random integration of introduced DNA (via cellular repair or recombination mechanisms) has been shown to occur. Integration has been documented to occur at low levels for nonintegrating vectors such as naked plasmid DNA and recombinant adenovirus vectors (20, 58)
An important question then arises: what is the potential risk (mutagenic or oncogenic) associated with such low-frequency events? For plasmid DNA gene transfer following intramuscular injection, theoretical calculations based on experimental data have allowed investigators to assert that the integration frequency (mutation rate) is at least 3,000 times below the spontaneous mutation rate (31, 35). When these same sorts of considerations are applied to the present rAAV data, we can derive similar numbers. Based on assay sensitivity, we would have detected a maximum level of 75 integrated vector copies in 20 ng of genomic DNA, which corresponds to a maximum of 3,750 integration events per microgram of DNA (assuming integration of a single vector copy per site). Moreover, assuming that 1 µg of mouse DNA is equivalent to 1.67 x 105 genomes and that there are 7.5 x 104 genes per mouse genome (35), there would then be a maximum of 2.9 x 10-7 integrated rAAV vector copies/gene. Comparison of this mutation rate to the generally accepted spontaneous mutation rate of a single base in the mammalian genome (1 x 10-5 mutations/gene) (11) reveals that the maximum mutation rate induced by rAAV integration in this study is at least 33 times lower than the spontaneous mutation rate. This value is likely an overestimation of the actual rate, because rAAV vectors typically integrate at a single locus in head-to-tail tandem arrays. Thus, assuming a minimum of two to four rAAV vectors integrated into a single site, the maximum mutation rate is further refined to a level at least 66 to 132 times lower than that of the spontaneous mutation rate. Although these numbers are slightly lower than those estimated for plasmid DNA, it should be noted that rAAV vectors uniformly have a higher initial transduction frequency than plasmid DNA, which leads to a much higher residual vector copy number in muscle (7).
Alternatively, rAAV integration can be viewed as an insertional mutagenic event due to possible tumor suppressor gene inactivation. To derive an estimate of the maximum number of rAAV integration events per injected muscle, we can assume that there are 1.6 million nuclei per quadriceps muscle (12, 24), yielding 9.6 µg of genomic DNA. Therefore, our data predict that the maximum number of integration events per injected muscle would be 3.6 x 104 integration events (3,750 integration events per microgram). The probability that integration of an intracellular DNA molecule could inactivate a tumor suppressor gene is estimated to be at most 10-6 to 10-9 (49). Thus, the probability of a single mutagenic event occurring in a tumor suppressor gene (worst-case scenario) in the entire injected muscle mass is at most between 0.0036% and 3.6%. Moreover, since only a small percentage of cells in the body can actually become oncogenic and since carcinogenesis itself is a multistep process, the possibility that rAAV integration occurs at the maximum defined levels determined herein appears to be remote.
In conclusion, our data suggest that rAAV genomes persist in mouse muscle as transcriptionally active large and small concatameric episomes. When considered in the context of potential and theoretical integration events, rAAV vectors show similarities to plasmid DNA in muscle. The nature of rAAV vector genome persistence in tissues other than muscle and species other than mice remains to be defined. It is worth noting that rAAV vectors currently in use in human clinical trials have been very well tolerated, even at high doses given by different routes over extended periods of time (23, 26, 52, 53).
This work was funded by grants from the National Institutes of Health.
|
|
|---|
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»