Frequency and spectrum of genomic integration of recombinant adeno-associated virus serotype 8 vector in neonatal mouse liver

Neonatal injection of recombinant adeno-associated virus serotype 8 (rAAV8) vectors results in widespread transduction in multiple organs and therefore holds promise in neonatal gene therapy. On the other hand, insertional mutagenesis causing liver cancer has been implicated in rAAV-mediated neonatal gene transfer. Here, to better understand rAAV integration in neonatal livers, we investigated the frequency and spectrum of genomic integration of rAAV8 vectors in the liver following intraperitoneal injection of 2.0 × 10¹¹ vector genomes at birth. This dose was sufficient to transduce a majority of hepatocytes in the neonatal period. In the first approach, we injected mice with a β-galactosidase-expressing vector at birth and quantified rAAV integration events by taking advantage of liver regeneration in a chronic hepatitis animal model and following partial hepatectomy. In the second approach, we performed a new, quantitative rAAV vector genome rescue assay by which we identified rAAV integration sites and quantified integrations. As a result, we find that at least ~0.05% of hepatocytes contained rAAV integration, while the average copy number of integrated double-stranded vector genome per cell in the liver was ~0.2, suggesting concatemer integration. Twenty-three of 34 integrations (68%) occurred in genes, but none of them were near the mir-341 locus, the common rAAV integration site found in mouse hepatocellular carcinoma. Thus, rAAV8 vector integration occurs preferentially in genes at a frequency of 1 in approximately 10³ hepatocytes when a majority of hepatocytes are once transduced in the neonatal period. Further studies are warranted to elucidate the relationship between vector dose and integration frequency or spectrum.     

Structure and distribution of inverted repeats (palindromes)

The size and distribution of renatured inverted repeats (palindromes) in Mus musculus DNA were examined by electron microscopy (EM). The majority (85%) of palindromes were found to be clustered in about one half of the DNA strands. The rest of the DNA strands were seen with a solitary looped structure. — The unlooped palindromes constituted 53% of all palindromes and were always clustered. There was a significant reduction in the number of unlooped palindromes in comparison to D. melanogaster DNA (Biezunski, 1981) and as a result the palindrome clusters were smaller and contained 2–8 palindromes [4–16 inverted repeats (ir)] per DNA strand. The looped palindromes had a wide and regular distribution with spacing lengths similar to those found in D. melanogaster DNA, and showed some periodicity. The average spacing between centers of all palindromes (inside a cluster) was 4.325 kb, and between centers of looped palindromes 8.544 kb. — The lengths of the ir of unlooped and looped palindromes were grouped (similar to D. melanogaster DNA) in one size-class with a range of 30–240 bp and an average length of 130 bp. Longer ir were also observed and the average length of ir in unlooped palindromes was 186 bp, in looped 588 bp, and the total average length was 375 bp. — It was calculated that there are about 224,000–320,000 palindromes (ir pairs) in the mouse genome, with the spacing between centers of all palindromes about 13-9 kb in length. — In high molecular weight mouse DNA, complex looped structures composed of rows of 5–8 looped palindromes one on top the other, formed by renaturation of multiple ir, were observed. It is suggested, that clustered repetitive sequences, in direct and inverted orientation, might be of one family and homologous to one another, and be able to reassociate, in vitro and in vivo, into structures of different forms, which could function as binding sites for various regulatory proteins during mouse development.     

Structure and distribution of inverted repeats (palindromes)

The size and distribution of renatured inverted repeats (palindromes) in D. melanogaster DNA were studied by electron microscopy (EM). The results of these studies differ from the previously published observations regarding the number, distribution and the size of inverted repeats (ir) present in DNA. -1. In contrast to the previous published observation almost all (96%) of the ir were found in crowded clusters. The DNA strands with clustered palindromes contained 2-21 palindromes (4-42 ir), with an average of 7.25 palindromes (14.5 ir) per strand. No correlation could be found between the length of the DNA strands and the number of ir per strand. -2, Also contrary to some previously published results, most (80%) of the ir formed on renaturation unlooped palindromes and these were always clustered. Looped palindromes (hairpins, formed by renaturation of ir separated by a non-homologous sequence long enough to be seen in EM as single-stranded loop) were found 1-2 per DNA strand, as part of clusters or as solitary palindromes in a DNA strand. The average spacing length (inside clusters) between centers of all palindromes was 2.349 kb, and between centers of looped palindromes 7.6 kb. - 3. The length of the ir was found to be smaller than documented in most of the previously published results. The majority, 80-90%, of the ir found in the unlooped and looped palindromes, respectively, belonged to one main-size class with a range of 30-210 bp and an average length of 100 bp, but longer ir were also observed. The average length of the ir in unlooped palindromes was 124 bp, in looped 244 bp, and the total average was 148 bp - 4. It was calculated that there are about 30,000 palindromes (60,000 ir) in the D, melanogaster genome, of which about 24,000 are unlooped and 6,000 looped, with the spacing between centers of all palindromes averaging about 4.4 kb in length.     

Isolation of Recombinant Adeno-Associated Virus Vector-Cellular DNA Junctions from Mouse Liver

Recombinant adeno-associated virus (rAAV) vectors allow for sustained expression of transgene products from mouse liver following a single portal vein administration. Here a rAAV vector expressing human coagulation factor F.IX (hF.IX), AAV-EF1alpha-F.IX (hF.IX expression was controlled by the human elongation factor 1alpha [EF1alpha] enhancer-promoter) was injected into mice via the portal vein or tail vein, or directly into the liver parenchyma, and the forms of rAAV vector DNA extracted from the liver were analyzed. Southern blot analyses suggested that rAAV vector integrated into the host genome, forming mainly head-to-tail concatemers with occasional deletions of the inverted terminal repeats (ITRs) and their flanking sequences. To further confirm vector integration, we developed a shuttle vector system and isolated and sequenced rAAV vector-cellular DNA junctions from transduced mouse livers. Analysis of 18 junctions revealed various rearrangements, including ITR deletions and amplifications of the vector and cellular DNA sequences. The breakpoints of the vector were mostly located within the ITRs, and cellular DNA sequences were recombined with the vector genome in a nonhomologous manner. Two rAAV-targeted DNA sequences were identified as the mouse rRNA gene and the alpha1 collagen gene. These observations serve as direct evidence of rAAV integration into the host genome of mouse liver and allow us to begin to elucidate the mechanisms involved in rAAV integration into tissues in vivo.     

A limited number of transducible hepatocytes restricts a wide-range linear vector dose response in recombinant adeno-associated virus-mediated liver transduction

Recombinant adeno-associated virus (rAAV) vectors are promising vehicles for achieving stable liver transduction in vivo. However, the mechanisms of liver transduction are not fully understood, and furthermore, the relationships between rAAV dose and levels of transgene expression, total number of hepatocytes transduced, and proportion of integrated vector genomes have not been well established. To begin to elucidate the liver transduction dose response with rAAV vectors, we injected mice with two different human factor IX or Escherichia coli lacZ-expressing AAV serotype 2-based vectors at doses ranging between 4.0 x 10(8) and 1.1 x 10(13) vector genomes (vg)/mouse, in three- to sixfold increments. A 2-log-range linear dose-response curve of transgene expression was obtained from 3.7 x 10(9) to 3.0 x 10(11) vg/mouse. Vector doses above 3.0 x 10(11) vg/mouse resulted in disproportionately smaller increases in both the number of transduced hepatocytes and levels of transgene expression, followed by saturation at doses above 1.8 x 10(12) vg/mouse. In contrast, a linear increase in the number of vector genomes per hepatocyte was observed up to 1.8 x 10(12) vg/mouse concomitantly with enhanced vector genome concatemerization, while the proportion of integrated vector genomes was independent of the vector dose. Thus, the mechanisms that restrict a wide-range linear dose response at high doses likely involve decreased functionality of vector genomes and restriction of transduction to fewer than 10% of total hepatocytes. Such information may be useful to determine appropriate vector doses for in vivo administration and provides further insights into the mechanisms of rAAV transduction in the liver.     

Increased common fragile site expression, cell proliferation defects, and apoptosis following conditional inactivation of mouse Hus1 in primary cultured cells

Targeted disruption of the mouse Hus1 cell cycle checkpoint gene results in embryonic lethality and proliferative arrest in cultured cells. To investigate the essential functions of Hus1, we developed a system for the regulated inactivation of mouse Hus1 in primary fibroblasts. Inactivation of a loxP site-flanked conditional Hus1 allele by using a cre-expressing adenovirus resulted in reduced cell doubling, cell cycle alterations, and increased apoptosis. These phenotypes were associated with a significantly increased frequency of gross chromosomal abnormalities and an S-phase-specific accumulation of phosphorylated histone H2AX, an indicator of double-stranded DNA breaks. To determine whether these chromosomal abnormalities occurred randomly or at specific genomic regions, we assessed the stability of common fragile sites, chromosomal loci that are prone to breakage in cells undergoing replication stress. Hus1 was found to be essential for fragile site stability, because spontaneous chromosomal abnormalities occurred preferentially at common fragile sites upon conditional Hus1 inactivation. Although p53 levels increased after Hus1 loss, deletion of p53 failed to rescue the cell-doubling defect or increased apoptosis in conditional Hus1 knockout cells. In summary, we propose that Hus1 loss leads to chromosomal instability during DNA replication, triggering increased apoptosis and impaired proliferation through p53-independent mechanisms.     

Vectors with restriction-site banks. I. pJRD158, a 3903-bp plasmid containing 28 unique cloning sites

A DNA fragment has been constructed that contains many unique cloning sites not present in currently used Escherichia coli plasmid cloning vehicles. Insertion of this fragment into a modified version of pBR322 results in an AmpRTetR vector (pJRD158) of 3903 bp containing 28 unique cloning sites, four "almost unique" cloning sites, and eight unassigned unique 6-bp palindromes. The plasmid has the additional advantages of very high copy number and altered incompatibility. The latter permits it to be stably maintained in the same host as pBR322.     

Parameters determining the efficiency of gene targeting in the moss Physcomitrella patens

In the moss Physcomitrella patens, transforming DNA containing homologous sequences integrates predominantly by homologous recombination with its genomic target. A systematic investigation of the parameters that determine gene targeting efficiency shows a direct relationship between homology length and targeting frequency for replacement vectors (a selectable marker flanked by homologous DNA). Overall homology of only 1 kb is sufficient to achieve a 50% yield of targeted transformants. Targeting may occur through homologous recombination in one arm, accompanied by non-homologous end-joining by the other arm of the vector, or by allele replacement following two homologous recombination events. Allele replacement frequency depends on the symmetry of the targeting vector, being proportional to the length of the shorter arm. Allele replacement may involve insertion of multiple copies of the transforming DNA, accompanied by ectopic insertions at non-homologous sites. Single-copy and single insertions at targeted loci (targeted gene replacements, ‘TGR’) occur with a frequency of 7–20% of all transformants when the minimum requirements for allele replacement are met. Homologous recombination in Physcomitrella is substantially more efficient than in any multicellular eukaryote, recommending it as the outstanding model for the study of homologous recombination in plants.     

DNA instability at chromosomal fragile sites in cancer

Human chromosomal fragile sites are specific genomic regions which exhibit gaps or breaks on metaphase chromosomes following conditions of partial replication stress. Fragile sites often coincide with genes that are frequently rearranged or deleted in human cancers, with over half of cancer-specific translocations containing breakpoints within fragile sites. But until recently, little direct evidence existed linking fragile site breakage to the formation of cancer-causing chromosomal aberrations. Studies have revealed that DNA breakage at fragile sites can induce formation of RET/PTC rearrangements, and deletions within the FHIT gene, resembling those observed in human tumors. These findings demonstrate the important role of fragile sites in cancer development, suggesting that a better understanding of the molecular basis of fragile site instability is crucial to insights in carcinogenesis. It is hypothesized that under conditions of replication stress, stable secondary structures form at fragile sites and stall replication fork progress, ultimately resulting in DNA breaks. A recent study examining an FRA16B fragment confirmed the formation of secondary structure and DNA polymerase stalling within this sequence in vitro, as well as reduced replication efficiency and increased instability in human cells. Polymerase stalling during synthesis of FRA16D has also been demonstrated. The ATR DNA damage checkpoint pathway plays a critical role in maintaining stability at fragile sites. Recent findings have confirmed binding of the ATR protein to three regions of FRA3B under conditions of mild replication stress. This review will discuss recent advances made in understanding the role and mechanism of fragile sites in cancer development.     

Integration of hepatitis B virus: analysis of unoccupied sites.

Hepatitis B virus (HBV) sequences integrated in the PLC/PRF/5 cell line (Alexander cells), which was derived from a human primary liver carcinoma, were previously extensively studied. Here we describe the analysis of the unoccupied sites of two linearly integrated forms of HBV DNA, AL-14 and AL-26, that were characterized previously. No major cellular DNA rearrangements were seen at the integration sites except for small deletions of host sequences: 2 kilobases of DNA in AL-14 and 17 base pairs (bp) in AL-26. The unoccupied site of AL-26 was found to be missing 182 bp, which previously mapped next to the right end of the integration sites of several independent clones. These were believed to be of cellular origin, but we show here that these 182 bp are in fact from unusual HBV sequences. Surprisingly, a region of this newly detected HBV DNA sequence is more homologous to that of woodchuck HBV DNA. Our analysis shows that the normal counterparts of both AL-14 and AL-26 contain minisatellite-like repetitive sequences. Based on the data presented here and our previous finding of HBV DNA integration at satellite sequences, we propose that genomic simple repetitive sequences are hot spots for HBV DNA integration.     

Site-specific genomic integration produces therapeutic Factor IX levels in mice

We used the integrase from phage phiC31 to integrate the human Factor IX (hFIX) gene permanently into specific sites in the mouse genome. A plasmid containing attB and an expression cassette for hFIX was delivered to the livers of mice by using high-pressure tail vein injection. When an integrase expression plasmid was co-injected, hFIX serum levels increased more than tenfold to approximately 4 microg/ml, similar to normal FIX levels, and remained stable throughout the more than eight months of the experiment. hFIX levels persisted after partial hepatectomy, suggesting genomic integration of the vector. Site-specific integration was proven by characterizing and quantifying genomic integration in the liver at the DNA level. Integration was documented at two pseudo-attP sites, native sequences with partial identity to attP, with one site highly predominant. This study demonstrates in vivo gene transfer in an animal by site-specific genomic integration.     

Improved measurement of dibenzo[a,l]pyrene-induced abasic sites by the aldehyde-reactive probe assay

Dibenzo[a,l]pyrene (DB[a,l]P) induces abundant amounts of depurinating adducts that spontaneously dissociate to form abasic sites in DNA. However, several previous studies that used the aldehyde-reactive probe (ARP) assay, could not verify abasic site formation by DB[a,l]P. Therefore, we examined whether a modification of the ARP assay would allow greater quantification of abasic sites. A previous study indicated that the abasic site quantification is improved by letting ARP trap the nascent abasic sites in cells, before extracting DNA for the assay. To test whether the addition of ARP to the DB[a,l]P-DNA adduct-forming reaction would improve abasic site quantification, we treated calf thymus DNA (0.625 mg/mL) with DB[a,l]P (80 microM) and 3-methylcholanthrene-treated rat liver microsomes with or without ARP (3 mM). The inclusion of ARP in the adduct-forming reaction resulted in significantly greater detection of abasic sites (62 lesions/10(6) bp versus 3.7 lesions/10(6) bp). DB[a,l]P also induces DNA strand breaks. The strand breaks may occur at abasic sites and by other mechanisms, such as oxidative damage. ARP/O-methoxyamine-abasic site conjugates are refractory to strand breakage, however, ARP or O-methoxyamine (3-10 mM) could only partially protect DB[a,l]P-induced DNA degradation, presumably by protecting the abasic sites, but not the other strand breaks. These results suggest that if DNA strand breakages occur at the abasic sites or at bases flanking them, and the fragments are lost during DNA extraction, abasic site estimation could be compromised. To obtain an independent line of evidence for abasic site formation in DB[a,l]P-treated cells, mouse Mbeta16 fibroblasts were treated with DB[a,l]P and O-methoxyamine. O-Methoxyamine is known to potentiate cytotoxicity of abasic site-inducing chemicals by forming abasic site conjugates, which partially inhibits their repair. O-Methoxyamine was found to increase DB[a,l]P cytotoxicity in these cells, supporting the idea that DB[a,l]P formed abasic sites. In summary, the inclusion of ARP in the DB[a,l]P-DNA adduct-forming reaction traps and protects the nascent abasic sites, allowing an improved quantification of abasic sites.     

Multiplex Cre/lox recombination permits selective site-specific DNA targeting to both a natural and an engineered site in the yeast genome.

Variant lox sites having an altered spacer region (heterospecific lox sites) are not proficient for Cre-mediated recombination with the canonical 34 bp loxP site, but can recombine with each other. By placing different heterospecific lox sites at different genomic locations, Cre can catalyze independent DNA recombination events at multiple loci in the same cell without concern that unwanted inter-locus recombination events will be generated. Such heterospecific lox sites also allow Cre to specifically target efficient integration of exogenous DNA to endogenous lox-like sequences that naturally occur in the genome. Specific targeting occurs only with a DNA vector carrying a heterospecific lox site in which the spacer region has been redesigned to match the 'spacer' region of the targeted chromosomal element. Moreover, in cells expressing a catalytically active Cre recombinase, naturally occurring lox-like sequences can exhibit almost 20% mitotic recombination. Thus, in the same cell, heterospecific lox sites can be used independently at multiple loci for integration, for deletion and for enhanced mitotic recombination, thereby increasing the repertoire of genomic manipulations catalyzed by the Cre recombinase.     

Molecular cloning and expression of the functional gene encoding the M2 subunit of mouse ribonucleotide reductase: a new dominant marker gene.

Mammalian ribonucleotide reductase consists of two non-identical subunits, proteins M1 and M2. M2-related DNA sequences are present on mouse chromosomes 4, 7, 12 and 13. However, M2-overproducing mouse cells show amplification of a chromosome 12-specific, single 13 kb HindIII fragment, which probably represents the active gene. We have isolated this fragment from parental mouse cell DNA and used it to clone and characterize the functional M2 gene. The 5770 bp transcribed M2 sequence contains ten exons separated by nine 95-917 bp introns. The 501 bp of 5' flanking DNA is G + C rich and contains TTTAAA and CCAAT sequences as well as potential Sp1 binding sites. The M2-related sequence on chromosome 13, which contains only the last six exons and several internal rearrangements, is a pseudogene. Transfection of BALB/3T3 cells with the M2 gene resulted in stable transformants with a 10-fold reduction in sensitivity to hydroxyurea, compared to control cells. This confirmed that the cloned M2 genomic DNA represents the functional gene and conclusively establishes the link between hydroxyurea resistance and M2 expression in mammalian cells. M2 genomic DNA should be a valuable dominant, selectable marker for identifying and isolating stable co-transformants.     

How common are common fragile sites: variation of aphidicolin-induced chromosomal fragile sites in a population of the deer mouse (Peromyscus maniculatus)

Aphidicolin (APC)-induced chromosomal gaps and breaks were analyzed for ten deer mice (Peromyscus maniculatus) from a natural population. The FSM statistical methodology was used to identify fragile sites as chromosomal loci exhibiting significantly non-random numbers of gaps/breaks in each individual and enabled an assessment of variation in fragile sites among the individuals. The individual deer mice exhibited as few as 7 to as many as 19 of the populational total of 34 sites. Two sites were fragile in all individuals and 13 sites were fragile in single individuals only. Defined by populational frequencies of greater than 50%, high-frequency fragile sites constituted 26% of the populational total. Approximately 35% of the total fragile sites were fragile in 20–40% of the population (low-frequency fragile sites) and about 38% were fragile in single individuals only. Analysis of the data pooled over all individuals identified significantly non-random breakage at 80 sites, 47 of which were not identified as fragile in any single individual. It appears, therefore, that fragile site identifications from pooled data have fostered an inflated estimate of the numbers and frequencies of common fragile sites. Comparison of the fragile site and spontaneous breakage (control) data suggest that APC-induced fragile sites represent regions of chromosomes that experience elevated levels of somatic mutation. Additionally, the occurrence of APC-induced fragile sites at or near the interstitial breakpoints of two pericentric-inversion polymorphisms in this population supports the hypothesis that fragile sites experience an increased rate of meiotic chromosomal mutation and are predisposed to undergo phylogenetic rearrangement. Received: 22 January 1997 / Accepted: 24 February 1997