Concise Analysis of Single-Stranded DNA of Recombinant Adeno-Associated Virus By Automated Electrophoresis System

Abstract

Recombinant adeno-associated virus (rAAV) is a prominent viral vector currently available for human gene therapy. The diameter of the rAAV capsid is ∼25 nm, and a positive or negative single-stranded DNA is packaged within the vector capsid. In this report, we describe a concise method to examine the extracted rAAV genome using an automated electrophoresis system. The rAAV genome, prepared from vector particles through either heat treatment at 95°C for 10 min or the phenol–chloroform extraction method, was analyzed using an automated electrophoresis system under denaturation conditions. The heat treatment protocol demonstrated a comparable yield with the phenol–chloroform extraction protocol, and the quantified amounts of the rAAV genome obtained using the automated electrophoresis system were consistent with those quantitated by quantitative PCR. Additionally, crude rAAV extractions could also be analyzed by the automated electrophoresis system after DNase I treatment. These results indicated that this simple and quick analysis using automated electrophoresis is highly useful for confirming the purity and integrity of the rAAV genome.

INTRODUCTION

Recombinant adeno-associated virus (rAAV) has been widely used in gene therapy as a vector to transfer genes to patients with genetic disorders.^1,2 rAAV has several advantages, including its ability to infect dividing and nondividing cells and its low potential for causing an inflammatory response, although AAV2 has recently been reported as a suspicious etiological factor causing severe acute hepatitis in children.^3
–5 The wild-type AAV genome is a single-stranded DNA (ssDNA) of ∼4.7 kb encoding Rep, Cap, membrane-associated accessory protein,⁶ and assembly-activating protein⁷ enclosed within two, 145-base, inverted terminal repeats (ITRs).^8,9

rAAV is constructed by replacing viral genes with the gene of interest (GOI) to be expressed in target cells, thus the packaging capacity is limited to ∼5 kb.^1,2,10,11 The clinical dose of rAAV is determined based on the vector genome (vg) titer per mL and it requires the availability of accurate quality control methods. Quantitative real-time PCR (qPCR) and droplet digital PCR (ddPCR) are the primary analytical methods widely used in quantification of rAAV owing to their robustness and wide measurement range.^{12

–19} Transgene encapsidation is not highly efficient in the vector manufacturing process, and the produced rAAV still contains empty and partially loaded components, which are considered to be impurities without therapeutic function.^20,21 Analytical ultracentrifugation (AUC)^22
–24 and transmission electron microscopy (TEM)^25,26 are currently used for analysis of the products. AUC is used to fractionate empty and filled rAAV derivatives and can also be used to determine the relative proportion of partially filled capsids in an AAV batch,²⁷ and TEM could morphologically discriminate empty particles from filled particles.^28,29 Mass photometry has also been reported as a method to characterize capsids based on their mass, distinguishing empty, partially filled, full, and overloaded capsids.^30
–32 Traditionally, verification of the rAAV genome size has been conducted through denaturing agarose gel electrophoresis.^33,34 These methods are time-consuming and have limited precision owing to the need for optimization of the analytical conditions. Recently, alternative approaches such as capillary electrophoresis (CE)^34,35 or electrophoresis-mediated microfluidics³⁶ have been reported.

The aim of our study was to examine whether the TapeStation system can be used for the analysis of rAAV ssDNA under denaturation conditions. It is frequently employed for sample quality control of both DNA and RNA for massive parallel sequencing. RNA was concisely analyzed under denaturation conditions using the RNA ScreenTape^® to analyze size, quantity, and integrity. In this study, we show that the analysis of the single-stranded form of the rAAV genome using the TapeStation system is fast, concise, and relatively robust for examining the purity, quantity, and integrity of the rAAV genome.

MATERIALS AND METHODS

Plasmid, DNA, and rAAV

The following reagents were used for automated electrophoresis: pAAV-ZsGreen1 (Takara Bio Inc., Kusatsu, Japan), a 1-kb DNA ladder (New England Biolabs Ltd., Ipswich, MA, USA), and an RNA ladder (Agilent Technologies, Inc., Santa Clara, CA, USA). To prepare the rAAV-ZsGreen1 DNA, 10 μg of the plasmid was digested with PvuII (New England Biolabs) and DraI (New England Biolabs), purified by phenol–chloroform extraction, and precipitated with ethanol.

Purified AAV8-AAT-FIXp and the crude extracts of rAAV2-CMV-ZsGreen1 prepared from VPC2.0 cells on days 1, 2, and 3 post-transfection were provided by the Manufacturing Technology Association of Biologics. The rAAV samples, 1908_rAAV1-CMV-ZsGreen1, 1909_rAAV2-CMV-ZsGreen1, and 1911_rAAV6-CMV-ZsGreen1, were also provided by the Manufacturing Technology Association of Biologics. For rAAV production, 293T cells (American Type Culture Collection, Manassas, VA, USA) were cultured in Dulbecco's modified Eagle's medium (DMEM) (Sigma-Aldrich Co., St. Louis, MO, USA) supplemented with 10% fetal bovine serum (Biowest, Bradenton, FL, USA). PEIpro (Polyplus Transfection, Vandœuvre-lès-Nancy, USA) was used for transfection with pAAV-ZsGreen1, pHelper (Takara Bio Inc.), and rep-cap plasmids (pRC1, pRC2-mi342, or pRC6; Takara Bio Inc.) in 293T cells. The ratio of PEIpro (Polyplus Transfection) to DNA weight was 1:1 in serum-free DMEM. After 72 h post-transfection, cells were lysed with Triton X-100 buffer (0.5% Triton X-100 and 2 mM MgCl₂ in phosphate-buffered saline [PBS]). The volume of the cell lysate was reduced to 1/8 using tangential flow filtration (Spectrum Inc., Rancho Dominguez, CA, USA) and applied to a HiTrap AVB Sepharose Column (GE Healthcare, Chicago, IL, USA). After washing, rAAV was eluted with 50 mM glycine-HCl (pH 2.7). Further purification of rAAV was performed by CsCl density gradient centrifugation at 148,500 g for 46 h at 21°C. Full rAAV particles were then subjected to dialysis against PBS and stored at −80°C until use.

Electrophoresis of DNA in the denaturation condition

The 4150 TapeStation System (Agilent Technologies, Inc.) was used for DNA sample electrophoresis. The system uses microfluidic technology to perform CE: nucleic acid molecules migrate through a gel matrix within the tape in an electric field, and the rate of migration is determined by the size and charge of the molecules. For DNA loading, an equal volume of the loading buffer (high-sensitivity [HS] RNA loading buffer) was added, and the mixture was incubated at 75°C for 5 min. Electrophoresis was performed using an HS RNA ScreenTape (Agilent Technologies, Inc.). The quantitative range was 0.5–10 ng/μL and sizing range was 0.1–6.0 kb.

rAAV genome DNA extraction from rAAV samples

For analysis of rAAV genomic DNA extracted from rAAV particles (1908_rAAV1-CMV-ZsGreen1, 1909_rAAV2-CMV-ZsGreen1, and 1911_rAAV6-CMV-ZsGreen1), 1.5 × 10⁹ rAAV capsids (20 μL) were destroyed by phenol–chloroform extraction, and the extracted genomic DNA was purified using AMPure XP beads (Beckman Coulter Inc., Brea, CA, USA). Phenol–chloroform extraction was performed twice with an equal volume of phenol: chloroform: isoamyl alcohol (25:24:1), followed by chloroform extraction twice. To 20 μL of genomic DNA (2.9 ng/μL), 36 μL of AMPure XP beads was added, and the genomic DNA was washed twice with 50 μL of 70% ethanol by magnetic precipitation and eluted with 10 μL of 10 mM Tris-HCl and 1 mM EDTA (pH 8.0).

Three rAAV genomic DNA extraction protocols were performed using 1908_rAAV1-CMV-ZsGreen1: protocol #1, heat treatment of purified rAAV capsids at 95°C for 10 min; protocol #2, phenol–chloroform extraction/purification with AMPure XP beads (Beckman Coulter Inc.), as described above; and protocol #3, heat treatment (95°C for 10 min) and purification with AMPure XP beads.

qPCR and ddPCR

All primers and probes were custom synthesized and high-performance liquid chromatography purified (Eurofins Genomics K.K., Tokyo, Japan). For rAAV-ZsGreen1 DNA, the forward primer 5′-TTCGTGATCACCGGCGAGGGCAT-3′, reverse primer 5′-CCGTACATGAAGGCGGCGGACAA-3′, and probe [FAM]AACCTGTGCGTGGTGGAGGGCGGC[BHQ1] were used for qPCR and ddPCR. For qPCR, pAAV-ZsGreen1 DNA (Takara Bio Inc.) digested with PvuII (New England Biolabs) was used to generate a standard curve. The reaction mixture comprised primers and probes at a final concentration of 0.4 and 0.1 μM for each using the QuantiTect Probe PCR Kit (204343; Qiagen, Hilden, Germany) in a final volume of 20 μL. qPCR was performed using the StepOnePlus system (Applied Biosystems, Waltham, MA, USA). The cycling conditions were as follows: 15 min at 95°C, followed by 40 cycles of a two-step thermal profile comprising 15 s at 94°C and 60 s at 60°C. Each DNA template was tested in triplicate at two dilutions.

For ddPCR, the reaction mixture comprised 10 μL of ddPCR Supermix for probes (Bio-Rad Laboratories, Hercules, CA, USA), 2 μL of primers (final concentration of each primer was 0.9 μM), 1 μL of probe (final concentration 0.25 μM), and 1 μL of the template diluted with TE buffer (10 mM tris(hydroxymethyl)aminomethane-HCl, pH 8.0, 1 mM ethylenediaminetetraacetic acid) containing 0.01% of PF-68 in a final volume of 20 μL. Each DNA template was tested in triplicate at two dilutions. The plates were transferred to a QX200 Automated Droplet Generator (Bio-Rad Laboratories). A 96-well plate containing the generated droplets was transferred to a C1000 Touch Thermal Cycler (Bio-Rad Laboratories). The cycling conditions were as follows: 10 min at 95°C, followed by 40 cycles of a two-step thermal profile comprising 30 s at 94°C and 60 s at 60°C. The plate was then transferred to a QX200 droplet reader (Bio-Rad Laboratories), and data analysis was performed using QuantaSoft software (version 1.7.4.0917; Bio-Rad Laboratories). The threshold separating the negative and positive droplets was manually set just above the cluster of negative droplets or just below the cluster of positive droplets, respectively.

RESULTS

Denaturation of double-stranded DNA into ssDNA

We first established the analytical conditions for ssDNA with a 1-kb DNA ladder containing 10 DNA fragments (0.5, 1.5, 1.0, 2.0, 3.0, 5.0, 6.0, 7.0, 8.0, and 10.0 kbp) using the TapeStation system. We analyzed double-stranded DNA (dsDNA) fragments using RNA ScreenTape for ssDNA analysis (Fig. 1). According to the manufacturer's instructions for RNA denaturation, we added a 0.5 volume of RNA sample buffer to the RNA sample and subjected the mixture to heat treatment at 72°C for 3 min. This denaturation condition was sufficient to linearize the RNA ladder (Fig. 1A, lane 1), but did not provide clear resolution for dsDNA with and without heat treatment (Fig. 1A, lanes 2 and 3). These results indicated that addition of a 0.5 volume of RNA sample buffer to the dsDNA mixture (1-kb ladder) was not sufficient for denaturation of dsDNA into a single-stranded form, even with heat treatment suitable for RNA denaturation. Next, we added an equal volume of RNA sample buffer to the 1-kb DNA ladder and compared the electropherograms with and without heat treatment at 75°C for 5 min (Fig. 1A, lanes 4 and 5). It should be noted that addition of an equal volume of the denaturation solution, along with heat treatment, provided a clear resolution of the 10 ssDNA fragments, and the mobility and fluorescence intensity of the DNA ladder changed compared with the bands without heat treatment (Fig. 1A). These results indicated that the dsDNA fragments were successfully denatured into a single-stranded form. The decrease in fluorescence intensity of ssDNA indicates the loss of intercalated fluorescent molecules between the complementary strands. In addition, the mobility of ssDNA fragments was consistent with the size of RNA (Fig. 1B). Taken together, RNA ScreenTape can also be used for dsDNA, although the denaturation condition for RNA (addition of 0.5 volume of HS RNA sample buffer with heat treatment at 73°C for 3 min) instructed by the manufacturer was not sufficient for complete denaturation of dsDNA (Fig. 1, lanes 3 and 5). For complete denaturation of dsDNA, addition of an equal volume of the sample buffer with incubation at 75°C for 5 min was required, and the denaturation conditions were valid for dsDNA up to 10 kb in length. Longer heat treatment, until 15 min, did not improve dsDNA denaturation (Fig. 1B).

Figure 1.

Analysis of DNA ladder by TapeStation with an HS RNA ScreenTape. (A) Analysis of a 1-kb DNA ladder with and without heat treatment. A 0.5 volume of HS RNA sample buffer was added to the HS RNA ladder, and the mixture was loaded after heat treatment at 72°C for 3 min, as per the manufacturer's protocol (lane 1). A 0.5 volume of HS RNA sample buffer was added to a 1-kb DNA ladder (12.5 ng/μL), and the mixture was loaded without heat treatment (lane 2). A 0.5 volume of HS RNA sample buffer was added to a 1-kb DNA ladder (12.5 ng/μL), and the mixture was loaded after heat treatment at 72°C for 3 min (lane 3). An equal volume of HS RNA sample buffer was added to a 1-kb DNA ladder (12.5 ng/μL), and the mixture was loaded without heat treatment (lane 4). An equal volume of HS RNA sample buffer was added to a 1-kb DNA ladder (12.5 ng/μL), and the mixture was loaded after heat treatment at 75°C for 5 min (lane 5). (B) An equal volume of HS RNA sample buffer was added to a 1-kb DNA ladder (12.5 ng/μL), and the mixture was loaded with heat treatment at 75°C for 0 min (lane 1), 3 min (lane 2), 5 min (lane 3), 7 min (lane 4), 10 min (lane 5), and 15 min (lane 6) using the HS RNA ScreenTape. Denaturation of a 1-kb dsDNA ladder after addition of HS RNA sample buffer; 0.5 volume of HS RNA sample buffer was added to the HS RNA ladder and the mixture was loaded after heat treatment at 72°C for 3 min, as per the manufacturer's protocol (lane 7). HS, high-sensitivity.

Denaturation of rAAV dsDNA and ITRs

Next, we confirmed the denaturation of the rAAV-ZsGreen1 dsDNA fragment generated from pAAV-ZsGreen1 digested with PvuII and DraI. As shown in Fig. 1, an equal volume of RNA sample buffer was added to the digested DNA fragments, and the mixture was incubated at 75°C for 5 min. The DNA mixture contained a 2,577-bp rAAV-ZsGreen1 fragment as well as 1,194-, 935-, 692-, 192-, 107-, and 19-bp fragments from the pUC portion (Fig. 2A). The electropherogram revealed that the mobility of rAAV-ZsGreen1 corresponded to a length of 2.6 kb based on the denatured 1-kb DNA ladder and RNA ladder (Fig. 2B). Quantification of rAAV-ZsGreen1 ssDNA using RNA ScreenTape was consistent with the amount of DNA applied after denaturation (Fig. 2C). These results indicate that the autoelectrophoresis system can be used to quantify rAAV ssDNA.

Figure 2.

Quantification of 2.6-kb rAAV-ZsGreen1 DNA fragment generated from pAAV-ZsGreen1. (A) The structure of pAAV-ZsGreen1 and its digestion sites for PvuII and DraI. (B) An equal volume of HS RNA sample buffer was added to a 1-kb DNA ladder (10 ng/μL), and the mixture was loaded after heat treatment at 72°C for 5 min (lane 1). An equal volume of HS RNA sample buffer was added to PvuII/DraI-digested pAAV-ZsGreen1 at the following concentrations: 0.11 ng/μL (lane 2), 0.23 ng/μL (lane 3), 0.45 ng/μL (lane 4), 0.68 ng/μL (lane 5), 0.90 ng/μL (lane 6), and 1.10 ng/μL (lane 7); 0.5 volume of HS RNA sample buffer was added to the HS RNA ladder, and the mixture was loaded after heat treatment at 72°C for 3 min (lane 8). (C) Detected concentration of the 2.6-kb rAAV-ZsGreen1 fragment by TapeStation using an RNA screen strip was plotted against the loaded concentration determined by Qubit. The RNA sample buffer contained a 25-base standard oligonucleotide (0.7 ng/μL) for quantification. AAV, adeno-associated virus; ITR, inverted terminal repeat.

The rAAV genome contained a GC-rich and palindromic ITR at both ends, which formed a T-shaped secondary structure (Fig. 3A). To address whether the denaturation conditions for dsDNA were also effective in linearizing the ITR secondary structure, we analyzed three synthesized oligonucleotides, ITR-oligo_1–50, ITR-oligo_1–100, and ITR-oligo_1–150, using the TapeStation system after addition of an equal volume of RNA sample buffer with or without heat treatment at 75°C for 5 min (Fig. 3B). Without heat treatment, the intensities of the 50-, 100-, and 150-base bands were faint, but became clearer after heat treatment, indicating that partial linearization of the ITR-oligos occurred without treatment and was successfully completed after heat treatment. The oligonucleotides exhibited size-dependent mobility compared with that of the ssDNA ladder. These results confirmed that the denaturation conditions were also effective for ITR structures.

Figure 3.

Analysis of synthesized ITR oligonucleotide. (A) Sequences of the synthesized ITR-oligo_1–50, ITR-oligo_1–100, and ITR-oligo_1–150. (B) ITR-oligos were analyzed after addition of an equal amount of HS RNA sample buffer into oligo (10 ng/μL) with or without heat treatment at 75°C for 5 min. ITR-oligo_1–50 without (lane 2) and with heat treatment (lane 3). ITR-oligo_1–100 without (lane 4) and with (lane 5) heat. ITR-oligo_1–150 without (lane 6) and with (lane 7) heat. An equal volume of HS RNA sample buffer was added to a 50-kb DNA ladder (10 ng/μL) with heat treatment at 75°C for 5 min (lane 1); 0.5 volume of HS RNA sample buffer was added to the HS RNA ladder, and the mixture was heated at 72°C for 3 min (lane 8).

Analysis of rAAV DNA extracted by phenol–chloroform

To apply this assay to rAAV, we prepared ssDNA from three purified rAAV samples, 1908_rAAV1-CMV-ZsGreen1, 1909_rAAV2-CMV-ZsGreen1, and 1911_rAAV6-CMV-ZsGreen1, using the phenol–chloroform extraction method. The extracted ssDNA was analyzed using a TapeStation system with RNA ScreenTape after the denaturation treatment described above (Fig. 4A, lanes 2, 3, and 4). The intact rAAV genome packaged in the capsids was 2.6 kb in size. The extracted rAAV DNA revealed clear resolution without other bands, indicating the absence of a truncated partial genome and DNA impurities. The concentrations of rAAV ssDNA were 3.15 ng/μL for 1908_rAAV1-CMV-ZsGreen1, 3.19 ng/μL for 1909_rAAV2-CMV-ZsGreen1, and 3.19 ng/μL for 1911_rAAV6-CMV-ZsGreen1. The same samples were also quantitated by qPCR, and the concentration of rAAV was 3.09 ng/μL for 1908_rAAV1-CMV-ZsGreen1, 3.08 ng/μL for 1909_rAAV2-CMV-ZsGreen1, and 4.10 ng/μL for 1911_rAAV6-CMV-ZsGreen1 (Fig. 4B). The concentration of rAAV DNA quantified by the TapeStation system was almost consistent with that determined by qPCR, although the effective concentration range of ssDNA measurements was not as broad as that of qPCR.

Figure 4.

Analysis of rAAV genome pretreated using the phenol–chloroform extraction method. (A) rAAV DNA was prepared from 1908_rAAV1-CMV-ZsGreen1, 1909_rAAV2-CMV-ZsGreen1, or 1911_rAAV6-CMV-ZsGreen1 by the phenol–chloroform extraction method, and the extracted rAAV DNA (10 ng/μL) was mixed with an equal volume of HS RNA sample buffer and treated at 75°C for 5 min. The mixtures were analyzed by TapeStation using HS RNA ScreenTape. Lane 1, 1-kb DNA ladder (10 ng/μL); lane 2, 1908_rAAV1-CMV-ZsGreen1; lane 3, 1909_rAAV2-CMV-ZsGreen1; lane 4, 1911_rAAV6-CMV-ZsGreen1; and lane 5, 0.5 volume of HS RNA sample buffer was added to the HS RNA ladder and the mixture was heated at 72°C for 3 min. (B) Concentration of 1908_rAAV1-CMV-ZsGreen1, 1909_rAAV2-CMV-ZsGreen1, and 1911_rAAV6-CMV-ZsGreen1 determined by qPCR and TapeStation using the HS RNA ScreenTape. (C) rAAV2-CMV-ZsGreen1 extracted from the CVP2.0 cells on days 1, 2, and 3 post-transfection was treated with DNase I. The rAAV2-CMV-ZsGreen1 DNA sample was analyzed by TapeStation using the HS RNA ScreenTape. Lane 1, 1-kb DNA ladder (10 ng/μL); lane 2, day 1 sample (0.5 × 10⁷ vg/μL); lane 3, day 2 sample (1 × 10⁸ vg/μL); and lane 4, day 3 sample (0.5 × 10⁸ vg/μL). The arrow indicates 2.6-kb rAAV2-CMV-ZsGreen1. (D) rAAV8-AAT-FIXp (1 × 10⁸ vg/μL) DNA was prepared by the phenol–chloroform extraction method, and the extracted rAAV DNA (10 ng/μL) was mixed with an equal volume of HS RNA sample buffer and treated at 75°C for 5 min. The mixtures were analyzed by TapeStation using the HS RNA ScreenTape. Lane 1, 1-kb DNA ladder (10 ng/μL); lane 2, 3 ng/μL 1908_rAAV1-CMV-ZsGreen1; and lane 3, 0.17 ng/μL, lane 4, 0.67 ng/μL, and lane 5, 1.7 ng/μL AAV8-AAT-FIXp. The arrow indicates 3.2-kb AAV8-AAT-FIXp DNA. qPCR, quantitative real-time PCR; rAAV, recombinant adeno-associated virus.

The crude extracts of rAAV2-CMV-ZsGreen1 prepared from VPC2.0 cells transfected with the triple plasmids on days 1, 2, and 3 (Fig. 4C, lanes 2, 3, and 4) post-transfection were analyzed by TapeStation following DNase I treatment. rAAV2-CMV-ZsGreen1 DNA (2.6 kb) could not be detected on day 1, but the 2.6-kb viral genome in the crude extracts could be detected on days 2 and 3. These results indicated that this analytical system was also useful to monitor rAAV genome DNA containing DNA impurities. Factor IX-associated hemophilia, known as hemophilia B or Christmas disease, results from a deficiency or dysfunction of clotting factor IX. Sustained therapeutic factor IX levels were attained after liver-directed gene therapy with AAV vectors in patients with hemophilia B.^37,38 We analyzed rAAV DNA carrying the FIX gene (AAV8-AAT-FIXp) as a GOI encapsulated in an AAV8 capsid (Fig. 4D).

Analysis of rAAV DNA extracted by heat treatment

Next, we examined whether the simple DNA extraction protocol was valid for the TapeStation assay. It has been reported that heat treatment of rAAV particles could release the encapsidated rAAV genome from the packaged capsids.³⁹ For extraction of the rAAV genome, vector particles were heated at 95°C for 10 min with or without purification using AMPure XP beads. The extracted rAAV genomes were then analyzed using the TapeStation system (Fig. 5A, B, lanes 1 and 3). For comparison, a sample prepared by the phenol–chloroform extraction/AMPure XP bead purification method was analyzed using the TapeStation system (Fig. 5B, lane 2). These three samples revealed similar yields of genomic DNA quantified using the TapeStation system, which were consistent with the results obtained by ddPCR (Fig. 5C). These results indicated that heat treatment at 95°C for 10 min in PBS is a simple and useful pretreatment of rAAV genome extraction and suitable for TapeStation analysis.

Figure 5.

rAAV DNA preparation and measurement on RNA ScreenTape. (A) Three DNA extraction methods for 1908_rAAV1-CMV-ZsGreen1: 1, heat treatment; 2, phenol–chloroform extraction/AMPure XP purification; and 3, heat treatment/AMPure XP purification. (B) Electrophoresis of 1908_rAAV1-CMV-ZsGreen1 DNA extracted by three methods; 0.5 volume of HS RNA sample buffer was added to a 1-kb DNA ladder (10 ng/μL) by heat treatment (method 1, lane 2); phenol–chloroform extraction/AMPure XP purification (method 2, lane 3); and heat treatment/AMPure XP purification (method 3, lane 4). (C) Quantitation of 1908_rAAV1-CMV-ZsGreen1 DNA by ddPCR and TapeStation with three pretreatment protocols. ddPCR, droplet digital PCR; PBS, phosphate-buffered saline; SD, standard deviation.

Serially diluted rAAV1-CMV-ZsGreen1 was analyzed to quantify the rAAV genome using the TapeStation system (Fig. 6A). The amount of rAAV1-CMV-ZsGreen1 DNA quantified using the TapeStation system was plotted against results obtained using ddPCR (Fig. 6B, C). Quantification was available in the range of 1 to 20 ng/μL with a linear regression (R ² = 0.99), but TapeStation measurement slightly overestimated the values (slope 1.1) compared with ddPCR. These results indicated that PBS did not affect the quantification of ssDNA using TapeStation. However, as per the manufacturer's instructions, attention should be paid to the buffer composition when RNA samples are analyzed using the TapeStation system.

Figure 6.

Quantification of packaged rAAV after heat treatment. (A) Procedure for sample prep. (B) An equal volume of HS RNA sample buffer was added to the 1-kb DNA ladder (10 ng/μL), and the mixture was loaded after heat treatment at 72°C for 5 min (lane 1). An equal volume of HS RNA sample buffer was added to 1–20 ng/μL 1908_rAAV1-CMV-ZsGreen1 (lanes 2–6). (C) 1908_rAAV1-CMV-ZsGreen1 concentration quantitated by TapeStation using RNA screen strip was plotted against those quantified by ddPCR.

DISCUSSION

In this study, we demonstrated the feasibility of using an automated electrophoresis device to determine the size, concentration, and integrity of ssDNA extracted from rAAV capsids. One critical quality attribute is the integrity of rAAV ssDNA packaged in the capsids, and it requires appropriate analytical techniques for quality control. The analysis of the size of single-stranded rAAV DNA packaged in capsids used to be performed using alkaline agarose gel electrophoresis.³³ Current analytical methods for rAAV genome integrity include high-throughput sequencing,^40
–42 two-dimensional digital PCR,⁴³ AUC,^22
–24 and capillary gel electrophoresis.^40,44 However, each of them lacks simplicity, speed, and user-friendliness. In this work, we demonstrate a concise and rapid analysis to characterize the integrity of rAAV genomes. The TapeStation system and its reagents for RNA analysis were used for the rAAV ssDNA analysis. In the electrophoretic analysis, formamide, formaldehyde, and glyoxal are commonly used for RNA.^45
–47 Formamide is ineffective against rigid secondary structures, whereas glyoxal exhibits a relatively strong denaturation effect on secondary structures. Although the specific denaturation reagents used in the system were not disclosed by the manufacturer, the denaturation conditions described in this report addressed concerns regarding the analysis of ssDNA with GC-rich secondary structures such as ITRs. The method proposed in this report offers a rapid and high-throughput alternative for confirming the purity and integrity of the rAAV genome, without the need for training in the use of specialized equipment.

In this system, SYBR Gold was used for detection of ssDNA. SYBR Gold stain is a proprietary, unsymmetrical cyanine dye that exhibits >1,000-fold fluorescence enhancement upon binding to nucleic acids and has a high quantum yield upon binding to ssDNA or to RNA. SYBR Gold stain is more sensitive than ethidium bromide, SYBR Green I stain, and SYBR Green II stain for detecting ssDNA and RNA.⁴⁸ In the TapeStation assay, ssDNA revealed a 1.18-fold higher mobility than that of RNA, and the fluorescence intensity of ssDNA was proportionally increased to the amount of applied ssDNA (R ² > 0.99) (data not shown).

The TapeStation system is frequently installed in laboratories and supports massive parallel sequencing for quality control of sequencing samples. Depending on the number of samples analyzed, analysis of 15 or fewer samples can be completed within 1 h, including the denaturation treatment. For the analysis, 1–25 ng (7 × 10⁹–1.8 × 10¹⁰ vg) of rAAV was sufficient to evaluate the ssDNA amount (Fig. 6C) when the rAAV genome size was 2.6 kb, although the concentration range was not broad compared with qPCR.

CONCLUSIONS

The use of an automated electrophoresis system, the TapeStation, along with appropriate denaturation conditions, enables efficient and accurate evaluation of rAAV ssDNA, providing valuable insights into the purity and integrity of the rAAV genome. This method will represent a valuable contribution to the field, streamlining the rAAV analysis process and enhancing the accessibility of quality control.

Footnotes

AUTHORs' CONTRIBUTIONS

Y.Y., K.Y., and K.U. designed the study and Y.Y. performed most experiments with assistance from K.H., N.H., and K.M.K. All the authors have read and approved the final manuscript.

AUTHOR DISCLOSURE

The authors declared no competing interests.

FUNDING INFORMATION

This work was funded by Japan Agency for Medical Research and Development (AMED) under grant numbers JP18ae0201001 and JP18ae0201002.

References

Naso

, Tomkowicz

, Perry

, 3rd, et al. Adeno-associated virus (AAV) as a vector for gene therapy. BioDrugs, 2017; 31(4):317–334; doi: 10.1007/s40259-017-0234-5.

Wang

, Tai

PWL

, Gao

. Adeno-associated virus vector as a platform for gene therapy delivery. Nat Rev Drug Discov, 2019; 18(5):358–378; doi: 10.1038/s41573-019-0012-9.

Barnes

, Draper

, Kurian

, et al. Analysis of AAV-extracted DNA by charge detection mass spectrometry reveals genome truncations. Anal Chem, 2023; 95(9):4310–4316; doi: 10.1021/acs.analchem.2c04234.

Servellita

, Gonzalez

, Lamson

, et al. Adeno-associated virus type 2 in US children with acute severe hepatitis. Nature, 2023; 7961: 574–580; doi: 10.1038/s41586-023-05949-1.

Morfopoulou

, Buddle

, Montaguth

OET

, et al. Genomic investigations of unexplained acute hepatitis in children. Nature, 2023; 617(7961):564–573; doi: 10.1038/s41586-023-06003-w.

Ogden

, Kelsic

, Sinai

, et al. Comprehensive AAV capsid fitness landscape reveals a viral gene and enables machine-guided design. Science, 2019; 366(6469):1139–1143; doi: 10.1126/science.aaw2900.

Sonntag

, Schmidt

, Kleinschmidt

. A viral assembly factor promotes AAV2 capsid formation in the nucleolus. Proc Natl Acad Sci U S A, 2010; 107(22):10220–10225; doi: 10.1073/pnas.1001673107.

Large

, Silveria

, Zane

, et al. Adeno-associated virus (AAV) gene delivery: Dissecting molecular interactions upon cell entry. Viruses, 2021; 13(7):1336; doi: 10.3390/v13071336.

Kuzmin

, Shutova

, Johnston

, et al. The clinical landscape for AAV gene therapies. Nat Rev Drug Discov, 2021; 20(3):173–174; doi: 10.1038/d41573-021-00017-7.

10.

, Yang

, Colosi

. Effect of genome size on AAV vector packaging. Mol Ther, 2010; 18(1):80–86; doi: 10.1038/mt.2009.255.

11.

, Samulski

. Engineering adeno-associated virus vectors for gene therapy. Nat Rev Genet, 2020; 21(4):255–272; doi: 10.1038/s41576-019-0205-4.

12.

Lock

, McGorray

, Auricchio

, et al. Characterization of a recombinant adeno-associated virus type 2 Reference Standard Material. Hum Gene Ther, 2010; 21(10):1273–1285; doi: 10.1089/hum.2009.223.

13.

Fagone

, Wright

, Nathwani

, et al. Systemic errors in quantitative polymerase chain reaction titration of self-complementary adeno-associated viral vectors and improved alternative methods. Hum Gene Ther Methods, 2012; 23(1):1–7; doi: 10.1089/hgtb.2011.104.

14.

Ayuso

, Blouin

, Lock

, et al. Manufacturing and characterization of a recombinant adeno-associated virus type 8 reference standard material. Hum Gene Ther, 2014; 25(11):977–987; doi: 10.1089/hum.2014.057.

15.

Werling

, Satkunanathan

, Thorpe

, et al. Systematic comparison and validation of quantitative real-time PCR methods for the quantitation of adeno-associated viral products. Hum Gene Ther Methods, 2015; 26(3):82–92; doi: 10.1089/hgtb.2015.013.

16.

D'Costa

, Blouin

, Broucque

, et al. Practical utilization of recombinant AAV vector reference standards: focus on vector genomes titration by free ITR qPCR. Mol Ther Methods Clin Dev, 2016; 5:16019; doi: 10.1038/mtm.2016.19.

17.

Sanmiguel

, Gao

, Vandenberghe

. Quantitative and digital droplet-based AAV genome titration. Methods Mol Biol, 2019; 1950:51–83; doi: 10.1007/978-1-4939-9139-6_4.

18.

Dobnik

, Kogovšek

, Jakomin

, et al. Accurate quantification and characterization of adeno-associated viral vectors. Front Microbiol, 2019; 10:1570; doi: 10.3389/fmicb.2019.01570.

19.

Wang

, Menon

, Shen

, et al. A qPCR method for AAV genome titer with ddPCR-level of accuracy and precision. Mol Ther Methods Clin Dev, 2020; 19:341–346; doi: 10.1016/j.omtm.2020.09.017.

20.

Penaud-Budloo

, François

, Clément

, et al. Pharmacology of recombinant adeno-associated virus production. Mol Ther Methods Clin Dev, 2018; 8:166–180; doi: 10.1016/j.omtm.2018.01.002.

21.

O'Connor

, Lutomski

, Jarrold

, et al. Lot-to-lot variation in adeno-associated virus serotype 9 (AAV9) preparations. Hum Gene Ther Methods, 2019; 30(6):214–225; doi: 10.1089/hgtb.2019.105.

22.

Burnham

, Nass

, Kong

, et al. Analytical ultracentrifugation as an approach to characterize recombinant adeno-associated viral vectors. Hum Gene Ther Methods, 2015; 26(6):228–242; doi: 10.1089/hgtb.2015.048.

23.

Nass

, Mattingly

, Woodcock

, et al. Universal method for the purification of recombinant AAV vectors of differing serotypes. Mol Ther Methods Clin Dev, 2018; 9:33–46; doi: 10.1016/j.omtm.2017.12.004.

24.

Maruno

, Ishii

, Torisu

, et al. Size distribution analysis of the adeno-associated virus vector by the c(s) analysis of band sedimentation analytical ultracentrifugation with multiwavelength detection. J Pharm Sci, 2023; 112(4):937–946; doi: 10.1016/j.xphs.2022.10.023.

25.

Chen

. Comparative observation of the recombinant adeno-associated virus 2 using transmission electron microscopy and atomic force microscopy. Microsc Microanal, 2007; 13(5):384–389; doi: 10.1017/s1431927607070808.

26.

Zhou

. Seeing engineered loops in a gene delivery vehicle by cryoEM. Structure, 2012; 20(8):1286–1288; doi: 10.1016/j.str.2012.07.008.

27.

Saleun

, Mas

, Le Roy

, et al. Analytical ultracentrifugation sedimentation velocity for the characterization of recombinant adeno-associated virus vectors sub-populations. Eur Biophys J, 2023; 52(4–5):367–377; doi: 10.1007/s00249-023-01650-3.

28.

, Chen

, Argento

, et al. Analytical strategies for quantification of adeno-associated virus empty capsids to support process development. Hum Gene Ther Methods, 2019; 30(4):144–152; doi: 10.1089/hgtb.2019.088.

29.

Gimpel

, Katsikis

, Sha

, et al. Analytical methods for process and product characterization of recombinant adeno-associated virus-based gene therapies. Mol Ther Methods Clin Dev, 2021; 20:740–754; doi: 10.1016/j.omtm.2021.02.010.

30.

Ebberink

, Ruisinger

, Nuebel

, et al. Assessing production variability in empty and filled adeno-associated viruses by single molecule mass analyses. Mol Ther Methods Clin Dev, 2022; 27:491–501; doi: 10.1016/j.omtm.2022.11.003.

31.

Wörner

, Snijder

, Friese

, et al. Assessment of genome packaging in AAVs using Orbitrap-based charge-detection mass spectrometry. Mol Ther Methods Clin Dev, 2022; 24:40–47; doi: 10.1016/j.omtm.2021.11.013.

32.

, Hwang

, Li

, et al. Rapid characterization of adeno-associated virus (AAV) gene therapy vectors by mass photometry. Gene Ther, 2022; 29(12):691–697; doi: 10.1038/s41434-021-00311-4.

33.

Dong

, Nakai

, Xiao

. Characterization of genome integrity for oversized recombinant AAV vector. Mol Ther, 2010; 18(1):87–92; doi: 10.1038/mt.2009.258.

34.

Barnes

, Draper

, Chen

, et al. Quantitative analysis of genome packaging in recombinant AAV vectors by charge detection mass spectrometry. Mol Ther Methods Clin Dev, 2021; 23:87–97; doi: 10.1016/j.omtm.2021.08.002.

35.

Hutanu

, Boelsterli

, Schmidli

, et al. Stronger together: Analytical techniques for recombinant adeno associated virus. Electrophoresis, 2022; 43(9–10):1107–1117; doi: 10.1002/elps.202100302.

36.

Coll De Peña

, Masto

, Atwood

, et al. Electrophoresis-mediated characterization of full and empty adeno-associated virus capsids. ACS Omega, 2022; 7(27):23457–23466; doi: 10.1021/acsomega.2c01813.

37.

Miesbach

, Meijer

, Coppens

, et al. Gene therapy with adeno-associated virus vector 5-human factor IX in adults with hemophilia B. Blood, 2018; 131(9):1022–1031; doi: 10.1182/blood-2017-09-804419.

38.

Nathwani

, Reiss

, Tuddenham

, et al. Long-term safety and efficacy of factor IX gene therapy in hemophilia B. N Engl J Med, 2014; 371(21):1994–2004; doi: 10.1056/NEJMoa1407309.

39.

Horowitz

, Rahman

, Bower

, et al. Biophysical and ultrastructural characterization of adeno-associated virus capsid uncoating and genome release. J Virol, 2013; 87(6):2994–3002; doi: 10.1128/jvi.03017-12.

40.

Tai

PWL

, Xie

, Fong

, et al. Adeno-associated virus genome population sequencing achieves full vector genome resolution and reveals human-vector chimeras. Mol Ther Methods Clin Dev, 2018; 9:130–141; doi: 10.1016/j.omtm.2018.02.002.

41.

Lecomte

, Saleun

, Bolteau

, et al. The SSV-Seq 2.0 PCR-free method improves the sequencing of adeno-associated viral vector genomes containing GC-rich regions and homopolymers. Biotechnol J, 2021; 16(1):e2000016; doi: 10.1002/biot.202000016.

42.

Zanker

, Lázaro-Petri

, Hüser

, et al. Insight and development of advanced recombinant adeno-associated virus analysis tools exploiting single-particle quantification by multidimensional droplet digital PCR. Hum Gene Ther, 2022; 33(17–18):977–989; doi: 10.1089/hum.2021.182.

43.

Furuta-Hanawa

, Yamaguchi

, Uchida

. Two-dimensional droplet digital PCR as a tool for titration and integrity evaluation of recombinant adeno-associated viral vectors. Hum Gene Ther Methods, 2019; 30(4):127–136; doi: 10.1089/hgtb.2019.031.

44.

Luo

, Guttman

. Rapid adeno-associated virus genome integrity analysis by capillary gel electrophoresis. AdV Bopharm Anal, 2020; 33:33–38.

45.

Masek

, Vopalensky

, Suchomelova

, et al. Denaturing RNA electrophoresis in TAE agarose gels. Anal Biochem, 2005; 336(1):46–50; doi: 10.1016/j.ab.2004.09.010.

46.

Rio

. Denaturation and electrophoresis of RNA with glyoxal. Cold Spring Harb Protoc, 2015; 2015(2):223–226; doi: 10.1101/pdb.prot081000.

47.

Rio

. Denaturation and electrophoresis of RNA with formaldehyde. Cold Spring Harb Protoc, 2015; 2015(2):219–222; doi: 10.1101/pdb.prot080994.

48.

Tuma

, Beaudet

, Jin

, et al. Characterization of SYBR Gold nucleic acid gel stain: A dye optimized for use with 300-nm ultraviolet transilluminators. Anal Biochem, 1999; 268(2):278–288; doi: 10.1006/abio.1998.3067.