Transgene repeats in aspen: molecular characterisation suggests simultaneous integration of independent T-DNAs into receptive hotspots in the host genome

Rearrangements of T-DNAs during genetic transformation of plants can result in the insertion of transgenes in the form of repeats into the host genome and frequently lead to loss of transgene expression. To obtain insight into the mechanism of repeat formation we screened 45 transgenic lines of aspen and hybrid aspen transformed with six different gene constructs. The frequency of T-DNA repeat formation among randomly screened transgenic lines was found to be about 21%. In ten transgenic lines direct repeats were detected. An inverted repeat was found in one other transgenic line. Sequencing of the junctions between the T-DNA inserts revealed identical residual right-border repeat sequences at the repeat junctions in all ten transgenic lines that had direct repeats. Formation of "precise" junctions based on short regions of sequence similarity between recombining strands was observed in three transgenic lines transformed with the same plasmid. Additional DNA sequences termed filler DNAs were found to be inserted between the T-DNA repeats at eight junctions where there was no similarity between recombining ends. The length of the filler DNAs varied from 4 to almost 300 bp. Small filler DNAs--a few base pairs long--were in most cases copied from T-DNA near the break points. The large filler sequences of about 300 bp in two transgenic lines were found to be of host plant origin, suggesting that transgene repeat formation occurred as a result of the simultaneous invasion of a receptive site in the host genome by two independent T-DNA strands. On the basis of the results obtained, and in the light of previous reports on T-DNA/plant DNA junctions in aspen and other crop plants, a mechanistic model for transgene rearrangement and filler formation is suggested.     

Sequence organization of the rat genome by electron microscopy.

The size and arrangement of repetitive and inverted repeat (foldback) sequences in rat DNA were studied by visualization of hybrid and heteroduplex structures in the electron microscope. The self-reassociation of repetitive sequence-bearing DNA strands often results in the formation of four-ended "H" structures, whose duplex regions equal the repetitive sequence length and can be measured in the electron microscope. In this way, it was determined that the average size of the class of numerous short repetitive sequences is 0.40 +/- 0.15 kbp. Heteroduplex structures were prepared between long whole DNA single strands and short repeat-sequence-bearing strands. The analysis of these structures confirms that the size of the repetitive sequences in 0.4 kbp on average. Length measurements between adjacent duplexes show that the average spacing between two interspersed repeats is at least 1.5-1.8 kbp. By examining 29.4-kbp single strands after brief renaturation, the size and distribution of foldback sequences were determined. There are 1.9 X 10(5) foldback apirs per rat genome, spaced an average of 9.7 kbp apart according to our measurement. Repetitive, inverted repeat and unique sequences are interspersed with each other in at least half the genome.     

Gene stability in transgenic aspen (Populus). I. Flanking DNA sequences and T-DNA structure

The stability of transgenes in the genome of transformed plants depends strongly on their correct physical integration into the host genome as well as on flanking target DNA sequences. For long-lived species like trees, however, no information is available so far concerning inactivation or loss of transgenes due to gene silencing or somatic genome rearrangement events. In this study, four independently transformed 35S-rolC transgenic hybrid aspen plants (Populus tremula L.?×?tremuloides Michx.), each harbouring one copy of the transgene, were investigated during continuous growth in the greenhouse. In one of these transgenic lines (Esch5:35S-rolC-##1) individuals frequently show phenotypic reversions, while in the remaining three lines (Esch5:35S-rolC-#3, -#5, -#16) the gene was essentially stable. Molecular analysis including PCR, Southern and Northern assays clearly showed that the transgene had been lost in the revertant tissue of the unstable line. Sequencing of T-DNA right and left borders, and flanking DNA regions, in all four transgenic aspen lines revealed no differences either in the type of flanking DNA (G-C to A-T ratio) or with respect to the presence of enhancers or MAR (matrix associated repeats)-like structures. Primers located within the left and right flanking regions in the three stable lines could be used to recover the target sites from the untransformed plants. This was not possible, however, with the unstable line, indicating that at least one flanking sequence does not derive from the plant target DNA but is of unknown origin. PCR using other primer pairs, and inverse PCR analysis, revealed an additional truncated T-DNA copy of 1050 nucleotides adjacent to the left border of the complete copy in this line. Sequencing of this truncated T-DNA revealed that it represented an inverted copy of part of the right half of the original construct. This special feature would allow the inverted repeat to pair with right border sequences of the complete copy. This would explain the frequently observed reversion resulting in transgene loss as due to intrachromosomal base-pairing leading to double-stranded loops of single-stranded DNA during mitotic cell divisions.     

Gene stability in transgenic aspen ( Populus ). I. Flanking DNA sequences and T-DNA structure

The stability of transgenes in the genome of transformed plants depends strongly on their correct physical integration into the host genome as well as on flanking target DNA sequences. For long-lived species like trees, however, no information is available so far concerning inactivation or loss of transgenes due to gene silencing or somatic genome rearrangement events. In this study, four independently transformed 35S-rolC transgenic hybrid aspen plants (Populus tremula L. × tremuloides Michx.), each harbouring one copy of the transgene, were investigated during continuous growth in the greenhouse. In one of these transgenic lines (Esch5:35S-rolC-##1) individuals frequently show phenotypic reversions, while in the remaining three lines (Esch5:35S-rolC-#3, -#5, -#16) the gene was essentially stable. Molecular analysis including PCR, Southern and Northern assays clearly showed that the transgene had been lost in the revertant tissue of the unstable line. Sequencing of T-DNA right and left borders, and flanking DNA regions, in all four transgenic aspen lines revealed no differences either in the type of flanking DNA (G-C to A-T ratio) or with respect to the presence of enhancers or MAR (matrix associated repeats)-like structures. Primers located within the left and right flanking regions in the three stable lines could be used to recover the target sites from the untransformed plants. This was not possible, however, with the unstable line, indicating that at least one flanking sequence does not derive from the plant target DNA but is of unknown origin. PCR using other primer pairs, and inverse PCR analysis, revealed an additional truncated T-DNA copy of 1050 nucleotides adjacent to the left border of the complete copy in this line. Sequencing of this truncated T-DNA revealed that it represented an inverted copy of part of the right half of the original construct. This special feature would allow the inverted repeat to pair with right border sequences of the complete copy. This would explain the frequently observed reversion resulting in transgene loss as due to intrachromosomal base-pairing leading to double-stranded loops of single-stranded DNA during mitotic cell divisions. Received: 9 June 1998 / Accepted: 6 October 1998     

The structures of integration sites in transgenic rice

Extensive genomic sequencing and sequence motif analysis have been conducted over the integration sites of two transgenic rice plants, #478 and #559, carrying the luciferase gene and/or hygromycin phosphotransferase gene. The transgenes reside in a region with inverted structure and a large duplication of rice genome over 2 kb. Integration was found at the AT-rich region and/or at the repetitive sequence region, including a SAR-like structure, retrotransposon and telomere repeats. The presence of a patch of sequence homology between plasmid and target DNA, and a small region of duplication involving the target DNA around the recombination site, implicated illegitimate recombination in the process of gene integration. Massive rearrangement of genomic DNA including deletion or translocation was also observed at the integration site and the flanking region of the transgene. The recognition sites of DNA topoisomerases I or II were observed in the rearranged sequences. Since only three junctions of transgenic rice were implicated in the illegitimate recombination and extensive rearrangement of the rice genome, rice protoplasts may be active in this process.     

Tandemly repeated transgenes of the human minisatellite MS32 (D1S8), with novel mouse gamma satellite integration.

The human hypervariable minisatellite MS32 has a well characterised internal repeat unit array and high mutation rates have been observed at this locus. Analysis of MS32 mutants has shown that male germline mutations are polarised to one end of the array and frequently involve complex gene conversion-like events, suggesting that tandem repeat instability may be modulated by cis-acting sequences flanking the locus. In order to investigate the processes affecting MS32 mutation rate and mechanism, we have created transgenic mice harbouring an MS32 allele. Here we describe the organisation of eight transgenic insertions. Analysis of these transgenic loci by MVR-PCR and structural analysis of the junctions between mouse flanking DNA and the transgenic loci has shed light on mechanisms of integration and rearrangement of the tandem repeated transgenes. Sequence analysis of the mouse DNA flanking these transgenes has shown that 5 of the 8 insertions have integrated into mouse gamma satellite repeated sequence. This suggests a non-random integration of the MS32 transgene construct into the mouse genome.     

Instability of long inverted repeats within mouse transgenes.

Various sequences in the mammalian genomes are unstable. One class of sequence arrangement is long inverted repeats, which are known to be unstable in bacteria and yeast. While in mammals some evidence suggests that short inverted repeats (<10 bp long) may show instability, nothing is known about the stability of long inverted repeats. Here we describe two unrelated multicopy transgenes in the mouse (loci 109 and OX1-5), each of which contains a long inverted repeat that shows substantial mitotic instability. This instability also occurs in the germline so that mutant transgenes appear within pedigrees at a high frequency. The mutation processes acting at these two inverted repeats are complex and can involve insertion or deletion, and can result in stabilization of the transgene. At transgene 109 mutational events range from very small rearrangements at the centre of the inverted repeat to complete transgene deletion. In addition we show that the rates of mutation at the inverted repeat of transgene OX1-5 can vary between the male and female germlines and between inbred strains of mice, suggesting the possibility of a genetic analysis to identify loci that modulate inverted repeat instability.     

A 21 kilobase-pair deletion/addition difference in the inverted repeat sequence of chloroplast DNA from Chlamydomonas eugametos and C. moewusii

Summary Our recent physical mapping of chloroplast DNA (cpDNA) from Chlamydomonas moewusii, a unicellular green alga which is interfertile with Chlamydomonas eugametos, has revealed a two-fold size difference between the inverted repeat sequences of these algae. With a size of 42 kbp, the inverted repeat of C. moewusii is the largest yet identified in any chloroplast genome. Here we have compared the arrangement of conserved sequences within the two algal inverted repeats by hybridizing cloned restriction fragments representing over 90% of these repeats to Southern blots of cpDNA digests from the two algae. We found that the size difference between the two algal inverted repeats is due to the presence of an extra DNA segment of 21 kilobase pairs (kbp) in C. moewusii. Except for this sequence, the C. moewusii inverted repeat is highly homologous to the entire C. eugametos repeat and the arrangement of conserved sequences in the two repeats is identical. Southern hybridizations with specific gene probes revealed that the conserved sequences include the rDNA region and the genes coding for the large subunit of ribulose 1,5 bisphosphate carboxylase-oxygenase (rbcL) and for the 32 kilodalton thylakoid membrane protein (psbA). With respect to the conserved sequences, the extra 21 kbp DNA segment of C. moewusii lies in the region of psbA, most probably slightly downstream from this gene.     

Twelve 43-base-pair repeats map in a cis-acting region essential for partition of plasmid mini-F.

The nucleotide sequence of the DNA region involved in partitioning of plasmid mini-F has been determined. The sequence consists of 12 direct tandemly arranged repeats of 43 base pairs (the two flanking repeats, 43 plus 1 base pairs) with extensive homology to each other. Each repeat contains an additional inverted repeat of 7 base pairs.     

Long-range structural effects in supercoiled DNA: statistical thermodynamics reveals a correlation between calculated cooperative melting and contextual influence on cruciform extrusion

We have previously described [K. M. Sullivan and D. M. J. Lilley (1986) Cell 47, 817-827] a set of sequences, called C-type inducing sequences, which cause cruciform extrusion by adjacent inverted repeats to occur by an abnormal kinetic pathway involving a large denatured region of DNA. In this paper we apply statistical thermodynamic DNA helix melting theory to these sequences. We find a marked correlation between the ability of sequences to confer C-type cruciform character experimentally and their calculated propensity to undergo cooperative melting, and no exceptions have been found. The correlations are both qualitative and quantitative. Thus the ColE1 flanking sequences behave as single melting units, while the DNA of the S-type plasmid pIRbke8 exhibits no propensity to melt in the region of the bke cruciform. The results of the calculations are also fully consistent with the following experimental observations: 1. the ability of the isolated colL and colR fragments of the ColE1 flanking sequences, as well as the short sequence col30, to confer C-type character; 2. C-type induction by an A + T rich Drosophila sequence; 3. low-temperature cruciform extrusion by an (AT)34 sequence; 4. the effect of changing sequences at a site 90 base pairs (bp) removed from the inverted repeat; 5. the effects of systematic deletion of the colL sequence; and 6. the effects of insertion of various sequences in between the colL sequence and the xke inverted repeat. These studies show that telestability effects on thermal denaturation as predicted from equilibrium helix melting theory of linear DNA molecules may explain all the features that are revealed by studying the extrusion of cruciforms in circular DNA molecules subjected to superhelical stress.     

Inverted repeat sequences in the drosophila genome

The properties of inverted repeat (foldback) sequences in Drosophila melanogaster DNA have been studied by HAP chromatography and electron microscope methods. Electron microscope observations show that there is a broad distribution of lengths of the duplex regions of the inverted repeats from very short to greater than 15 kb, with number and weight average values of 1.35 kb and 5.0 kb respectively. About 20% of the inverted repeats are separated by a single-strand spacer with lengths too short to observe, but the other 80% have spacers, P, with lengths ranging from 0.5 kb to greater than 30 kb. The number average and weight average spacer lengths for the total sample are 2.7 kb and 6.1 kb. With respect to the lengths of the spacers, P, between inverted repeats, the Drosophila genome differs from that of most organisms which have been studied where the spacers P are mostly too short to be measured. EM and HAP studies suggest that the average center-to-center spacing between sets of inverted repeats is 40–80 kb. The HAP studies show that there is a broad range of thermal stabilities for the duplexes formed by reassociation of inverted repeat sequences. Kinetic analysis shows that all of the frequency components of the Drosophila genome are present in the inverted repeats, the loops P, and the flanking sequences. There is a somewhat larger proportion of middle repetitive DNA in those inverted repeat duplexes which are resistant to digestion by Mung Bean Endonuclease I. These enzyme resistant duplexes comprise about 3% of the entire genome. It is estimated that there are approximately 2000–4000 inverted repeat pairs in the entire genome.     

Sequence analysis of the junctions between a large inverted repeat and single-copy regions in tobacco chloroplast DNA

Summary Tobacco chloroplast DNA contains a large inverted repeat sequence of 26 kilobase pairs (kbp). The inverted repeat is separated by 20 kbp small single-copy and 90 kbp large single-copy regions. We have cloned four DNA fragments containing each junction between the inverted repeat and the single-copy regions. The sequence analysis revealed the exact edges of the inverted repeat. A putative coding region for a ribosomal protein CS19 was found 4 base pairs (bp) away from the inverted repeat on the left margin of the large single-copy region. A sequence AGGAG, which is complementary to the 3 terminal sequence of tobacco chloroplast 16S rRNA, was found within the inverted repeat. A tRNA^His gene was found 5 bp away from the inverted repeat on the right-hand margin of the large single-copy region.     

Sequence, organization, and evolution of the A+T region of Drosophila melanogaster mitochondrial DNA

The long (4.6-kb) A+T region of Drosophila melanogaster mitochondrial DNA has been cloned and sequenced. The A+T region is organized in two large arrays of tandemly repeated DNA sequence elements, with nonrepetitive intervening and flanking sequences comprising only 22% of its length. The first repeat array consists of five repeats of 338-373 bp. The second consists of four intact 464-bp repeats and a fifth partial repeat of 137 bp. Three DNA sequence elements are found to be highly conserved in D. melanogaster and in several Drosophila species with short A+T regions. These include a 300-bp DNA sequence element that overlaps the DNA replication origin and two thymidylate stretches identified on opposite DNA strands. We conclude that the length heterogeneity observed in the A+T regulatory region in mitochondrial DNAs from the genus Drosophila results from the expansion (and contraction) of the number of repeated DNA sequence elements. We also propose that the 300-bp conserved DNA sequence element, in conjunction with another primary sequence determinant, perhaps the adjacent thymidylate stretch, functions in the regulation of mitochondrial DNA replication.      

Genomic organization of the gene that encodes the precursor to EGF-related peptides, exogastrula-inducing peptides, of the sea urchin Anthocidaris crassispina

Exogastrula-inducing peptides (EGIPs) were identified in embryos of the sea urchin Anthocidaris crassispina as polypeptides with structural similarity to epidermal growth factor (EGF) that severely affect gastrulation of sea urchin embryos to induce exogastrulation. Here we have obtained genomic clones for the EGIP precursor gene (EGIP) and determined its genomic organization. The EGIP gene spans the length of 9 kb in the genome and is composed of seven exons and six introns. Each of the four EGF motifs in the precursor protein is encoded by a single exon, and all the exon boundaries are in phase 1, suggesting that EGIP have been generated during evolution by duplication of an exon encoding a single ancient EGIP sequence. The 5'-flanking sequence of EGIP from -4372 to +194 revealed the presence of multiple repeat sequences including direct and inverted repeats as well as two clusters of GGGG/CCCC elements. The function of the upstream flanking region of EGIP was examined by introducing the gene constructs into embryos in which different regions from the flanking DNA were placed upstream to the GFP reporter gene. Systematic deletion of the upstream DNA revealed the presence of potent enhancer activity between -372 and -210.     

Sleeping Beauty transposition in the mouse genome is associated with changes in DNA methylation at the site of insertion

The Sleeping Beauty (SB) transposon (Tn) system is a nonviral gene delivery tool that has widespread application for transfer of therapeutic genes into the mammalian genome. To determine its utility as a gene delivery system, it was important to assess the epigenetic modifications associated with SB insertion into the genome and potential inactivation of the transgene. This study investigated the DNA methylation pattern of an SB Tn as well as the flanking genomic region at insertion sites in the mouse genome. The ubiquitous ROSA26 promoter and an initial part of the eGFP coding sequence in the SB Tn exhibited high levels of CpG methylation in transgenic mouse lines, irrespective of the chromosomal loci of the insertion sites. In contrast, no detectable CpG methylation in the endogenous mouse ROSA26 counterpart was observed in the same animals. Furthermore, significant hypomethylation was detected in neighboring chromosomal sequences of two unique SB Tn insertions compared to wild-type patterns. Taken together, these results suggest that SB Tn insertions into the mouse genome can be discriminated by DNA methylation machinery from an identical endogenous DNA sequence and can profoundly alter the DNA methylation status of the transgene cargo as well as flanking host genomic regions.     

Chromosomal mapping and zygosity check of transgenes based on flanking genome sequences determined by genomic walking.

Transgenes can affect transgenic mice via transgene expression or via the so-called positional effect. DNA sequences can be localized in chromosomes using recently established mouse genomic databases. In this study, we describe a chromosomal mapping method that uses the genomic walking technique to analyze genomic sequences that flank transgenes, in combination with mouse genome database searches. Genomic DNA was collected from two transgenic mouse lines harboring pCAGGS-based transgenes, and adaptor-ligated, enzyme restricted genomic libraries for each mouse line were constructed. Flanking sequences were determined by sequencing amplicons obtained by PCR amplification of genomic libraries with transgene-specific and adaptor primers. The insertion positions of the transgenes were located by BLAST searches of the Ensembl genome database using the flanking sequences of the transgenes, and the transgenes of the two transgenic mouse lines were mapped onto chromosomes 11 and 3. In addition, flanking sequence information was used to construct flanking primers for a zygosity check. The zygosity (homozygous transgenic, hemizygous transgenic and non-transgenic) of animals could be identified by differential band formation in PCR analyses with the flanking primers. These methods should prove useful for genetic quality control of transgenic animals, even though the mode of transgene integration and the specificity of flanking sequences needs to be taken into account.     

Stability of an inverted repeat in a human fibrosarcoma cell.

Deletions and rearrangements of DNA sequences within the genome of human cells result in mutations associated with human disease. We have developed a selection system involving a neo gene containing a DNA sequence inserted into the NcoI site that can be used to quantitatively assay deletion of this sequence from the chromosome. The spontaneous deletion from the neo gene of a 122 bp inverted repeat occurred at a rate of 2.1 x 10(-8) to <3.1 x 10(-9) revertants/cell/generation in three different cell lines. Deletion of the 122 bp inverted repeat occurred between 6 bp flanking direct repeats. Spontaneous deletion of a 122 bp non-palindromic DNA sequence flanked by direct repeats was not observed, indicating a rate of deletion of <3.1 x 10(-9) revertants/cell/generation. This result demonstrates that a 122 bp inverted repeat can exhibit a low level of instability in some locations in the chromosome of a human cell line.     

A new family of widely propagated MB1-repeats in the human genome

A new family of repeats--i.e. MB1 repeats family--the number of copies of which per a human genome constitutes a few hundreds of thousands of copies has been revealed in a human gemone by computer analysis of a noncanonical similarity of nucleic acid sequences. The numbers of that family of repeats have also been revealed in the genomes of mouse and rat, they have been identified as mirror--reflected copies--in purines and pyrimidines--of B1 repeats in the genome of mouse and the Alu repeats in the human genome. The MB1 repeats tend to remain most similar at a length of 70 b.p. They are not flanked by short repeats, neither contain poly(A) region at the 3' end, by which they differ from the repeats of the SINE family. It has been assumed that the member of the Alu repeats family and the MB1 repeats family can form a so called H-form of DNA. The mirror-reflected repeat family could have been formed by replication of parallel DNA strands.