Context dependence of meiotic recombination hotspots in yeast: the relationship between recombination activity of a reporter construct and base composition

Borde and colleagues reported that a reporter plasmid inserted at different genomic locations in Saccharomyces cerevisiae had different levels of meiotic recombination activity. We show that the level of recombination activity is very significantly correlated with the GC content of DNA sequences flanking the insertion.     

Diversity in G + C content at the third position of codons in vertebrate genes and its cause.

Correlation was positive between the G + C content at the codon third position in genes of vertebrates and the G + C content of the genome portion surrounding each gene. Exons of genes with a high G + C% at the codon 3rd position are surrounded by G + C-rich introns and G + C-rich flanking sequences, and those with a low G + C% at the position by A + T-rich introns and flanking sequences. Analysis of G + C content distribution along DNA sequences using a DNA Sequence Data Bank supported the view that the vertebrate genome is a mosaic of regions with clear differences in their G + C content. The biological significance of the variation in G + C content throughout the vertebrate genome is discussed in connection with chromosomal banding.     

Cruciform extrusion in plasmids bearing the replicative intermediate configuration of a poxvirus telomere

The transition from lineform DNA to cruciform DNA (cruciformation) within the cloned telomere sequences of the Leporipoxvirus Shope fibroma virus (SFV) has been studied. The viral telomere sequences have been cloned in recombination-deficient Escherichia coli as a 322 base-pair, imperfect palindromic insert in pUC13. The inverted repeat configuration is equivalent to the arrangement of the telomere structures observed within viral DNA replicative intermediates. A major cruciform structure in the purified recombinant plasmid has been identified and mapped using, as probes, the enzymes AflII, nuclease S1 and bacteriophage T7 endonuclease I. It was extruded from the central axis of the cloned viral inverted repeat and, by unrestricted branch migration, attained a size commensurate with the superhelical density of the plasmid molecule at native superhelical densities. This major cruciform extrusion event was the only detectable duplex DNA perturbation, induced by negative superhelical torsion, in the insert viral sequences. No significant steady-state pool of extruded cruciform was identified in E. coli. However, the identification of a major deletion variant generated even in the recombination-deficient E. coli strain DB1256 (recA recBC sbcB) suggested that the cruciform may be extruded transiently in vivo. The lineform to cruciform transition has been further characterized in vitro using two-dimensional agarose gel electrophoresis. The transition was marked by a high energy of formation (delta Gf = 44 kcal/mol), and an apparently low activation energy that enabled facile transitions at physiological temperatures provided there was sufficient torsional energy. By comparing cruciformation in a series of related bidirectional central axis deletions of the telomeric insert, it has been concluded that the presence of extrahelical bases in the terminal hairpin structures contributes substantially to the high delta Gf value. Also, viral sequences flanking the extruded cruciform were shown to influence the measured delta Gf value. Several general features of poxvirus telomere structure that would be expected to influence the facility of cruciform extrusion are discussed along with the implications of the observed cruciform transition event on the replicative process of poxviruses in vivo.     

Unpaired structures in SCA10 (ATTCT)n.(AGAAT)n repeats

A number of human hereditary diseases have been associated with the instability of DNA repeats in the genome. Recently, spinocerebellar ataxia type 10 has been associated with expansion of the pentanucleotide repeat (ATTCT)(n).(AGAAT)(n) from a normal range of ten to 22 to as many as 4500 copies. The structural properties of this repeat cloned in circular plasmids were studied by a variety of methods. Two-dimensional gel electrophoresis and atomic force microscopy detected local DNA unpairing in supercoiled plasmids. Chemical probing analysis indicated that, at moderate superhelical densities, the (ATTCT)(n).(AGAAT)(n) repeat forms an unpaired region, which further extends into adjacent A+T-rich flanking sequences at higher superhelical densities. The superhelical energy required to initiate duplex unpairing is essentially length-independent from eight to 46 repeats. In plasmids containing five repeats, minimal unpairing of (ATTCT)(5).(AGAAT)(5) occurred while 2D gel analysis and chemical probing indicate greater unpairing in A+T-rich sequences in other regions of the plasmid. The observed experimental results are consistent with a statistical mechanical, computational analysis of these supercoiled plasmids. For plasmids containing 29 repeats, which is just above the normal human size range, flanked by an A+T-rich sequence, atomic force microscopy detected the formation of a locally condensed structure at high superhelical densities. However, even at high superhelical densities, DNA strands within the presumably compact A+T-rich region were accessible to small chemicals and oligonucleotide hybridization. Thus, DNA strands in this "collapsed structure" remain unpaired and accessible for interaction with other molecules. The unpaired DNA structure functioned as an aberrant replication origin, in that it supported complete plasmid replication in a HeLa cell extract. A model is proposed in which unscheduled or aberrant DNA replication is a critical step in the expansion mutation.     

Transfer of chimeric plasmids among Salmonella typhimurium strains by P22 transduction

Salmonella typhimurium bacteriophage P22 transduced plasmids having P22 sequences inserted in the vector pBR322 with high frequency. Analysis of the structure of the transducing particle DNA and the transduced plasmids indicates that this plasmid transduction involves two homologous recombination events. In the donor cell, a single recombination between the phage and the homologous sequences on the plasmid inserted the plasmid into the phage chromosome, which was then packaged by headfuls into P22 particles. The transducing particle DNA contained duplications of the region of homology flanking the integrated plasmid vector sequences and lacked some phage genes. When these defective phage genomes containing the inserted plasmid infected a recipient cell, recombination between the duplicated regions regenerated the plasmid. A useful consequence of this sequence of events was that genetic markers in the region of homology were readily transferred from phage to plasmid. Plasmid transduction required homology between the phage and the plasmid, but did not depend on the presence of any specific P22 sequence in the plasmid. When the infecting P22 carried a DNA sequence homologous to the ampicillin resistance region of pBR322, the vector plasmid having no P22 insert could be transduced. P22-mediated transduction is a useful way to transfer chimeric plasmids, since most S. typhimurium strains are poorly transformed by plasmid DNA.     

Sequence and analysis of the 46.6-kb plasmid pA1 from Sphingomonas sp. A1 that corresponds to the typical IncP-1beta plasmid backbone without any accessory gene

Sphingomonas sp. A1 (strain A1) is capable of directly incorporating macromolecules (e.g., alginate) through the specialized import system--"super-channel." Here, we report the complete DNA sequence and genetic organization of plasmid pA1 from strain A1. Nucleotide sequence analysis revealed that pA1 comprises 46,557 bp encoding 49 open reading frames (ORFs) with 65% G+C content and abundant GCCG/CGGC motifs. Many predicted pA1 ORFs showed high similarity to pA81 ORFs; pA81 is supposedly a self-transmissible promiscuous incompatibility (Inc) group P-1beta plasmid. Unlike any reported IncP-1 plasmids, pA1 contains no inserted mobile genetic elements. The genetic organization and predicted pA1 ORFs showed greater similarity to the IncP-1beta plasmid backbone than to the IncP-1alpha plasmid backbone. pA1 contains restriction site-associated repeat sequences typical of the IncP-1beta but absent in the IncP-1alpha and delta subgroups. Thus, the overall pA1 structure corresponds to that of the typical IncP-1beta plasmids. Phylogenetic analysis of the replication-associated proteins suggested that pA1 may have diverged later along with the two IncP-1beta plasmids--pA81 and pB4. The 2.4-kb duplicates of stable inheritance genes klcAB and korC in pA1 possibly resulted from insertion and/or recombination events via the repeat sequences flanking these duplicates.     

Meiotic Gene Conversion and Crossing over between Dispersed Homologous Sequences Occurs Frequently in Saccharomyces cerevisiae

We have examined meiotic recombination between two defined leu2 heteroalleles present at the normal LEU2 locus and in leu2-containing plasmids inserted at four other genomic locations. In diploids where the two leu2 markers were present at allelic locations on parental homologs, the frequency of Leu2+ spores varied 38-fold, in a location-dependent manner. These results indicate that recombination in a genetic interval can be modulated by sequences at least 2.7 kb outside that interval. Leu2+ meiotic segregants were also recovered from diploids where LEU2 was marked with one heteroallele, and the other leu2 heteroallele was inserted at another genomic location. These products of ectopic interactions, between dispersed copies of leu2 sharing only 2.2 kb of homology, were recovered at a frequency comparable to that observed in corresponding allelic crosses. This high frequency of ectopic meiotic recombination was observed in crosses where both recombining partners could potentially pair with sequences at an allelic position. In addition, a significant fraction (22-50%) of these ectopic recombinants were associated with crossing over of flanking sequences.     

Directed insertion of a selectable marker into a circular plasmid of Borrelia burgdorferi.

Studies of the biology of Borrelia burgdorferi and the pathogenesis of Lyme disease are severely limited by the current lack of genetic tools. As an initial step toward facile genetic manipulation of this pathogenic spirochete, we have investigated gene inactivation by allelic exchange using a mutated borrelial gyrB gene that confers resistance to the antibiotic coumermycin A1 as a selectable marker. We have transformed B. burgdorferi by electroporation with a linear fragment of DNA in which this selectable marker was flanked by sequences from a native borrelial 26-kb circular plasmid. We have identified coumermycin A1-resistant transformants in which gyrB had interrupted the targeted site on the 26-kb plasmid via homologous recombination with the flanking sequences. Antibiotic resistance conferred by the mutated gyrB gene on the plasmid is dominant, and transformed spirochetes carrying this plasmid do not contain any unaltered copies of the plasmid. Coumermycin A1 resistance can be transferred to naive B. burgdorferi by transformation with borrelial plasmid DNA from the initial transformants. This work represents the first example of a directed mutation in B. burgdorferi whereby a large segment of heterologous DNA (gyrB) has been inserted via homologous recombination with flanking sequences, thus demonstrating the feasibility of specific gene inactivation by allelic exchange.     

Susceptibility to superhelically driven DNA duplex destabilization: a highly conserved property of yeast replication origins

Strand separation is obligatory for several DNA functions, including replication. However, local DNA properties such as A+T content or thermodynamic stability alone do not determine the susceptibility to this transition in vivo. Rather, superhelical stresses provide long-range coupling among the transition behaviors of all base pairs within a topologically constrained domain. We have developed methods to analyze superhelically induced duplex destabilization (SIDD) in genomic DNA that take into account both this long-range stress-induced coupling and sequence-dependent local thermodynamic stability. Here we apply this approach to examine the SIDD properties of 39 experimentally well-characterized autonomously replicating DNA sequences (ARS elements), which function as replication origins in the yeast Saccharomyces cerevisiae. We find that these ARS elements have a strikingly increased susceptibility to SIDD relative to their surrounding sequences. On average, these ARS elements require 4.78 kcal/mol less free energy to separate than do their immediately surrounding sequences, making them more than 2,000 times easier to open. Statistical analysis shows that the probability of this strong an association between SIDD sites and ARS elements arising by chance is approximately 4 × 10⁻¹⁰. This local enhancement of the propensity to separate to single strands under superhelical stress has obvious implications for origin function. SIDD properties also could be used, in conjunction with other known origin attributes, to identify putative replication origins in yeast, and possibly in other metazoan genomes.     

Nucleotide Composition Bias Affects Amino Acid Content in Proteins Coded by Animal Mitochondria

We show that in animal mitochondria homologous genes that differ in guanine plus cytosine (G + C) content code for proteins differing in amino acid content in a manner that relates to the G + C content of the codons. DNA sequences were analyzed using square plots, a new method that combines graphical visualization and statistical analysis of compositional differences in both DNA and protein. Square plots divide codons into four groups based on first and second position A + T (adenine plus thymine) and G + C content and indicate differences in amino acid content when comparing sequences that differ in G + C content. When sequences are compared using these plots, the amino acid content is shown to correlate with the nucleotide bias of the genes. This amino acid effect is shown in all protein-coding genes in the mitochondrial genome, including cox I, cox II, and cyt b, mitochondrial genes which are commonly used for phylogenetic studies. Furthermore, nucleotide content differences are shown to affect the content of all amino acids with A + T- and G + C-rich codons. We speculate that phylogenetic analysis of genes so affected may tend erroneously to indicate relatedness (or lack thereof) based only on amino acid content. Received: 3 July 1996 / Accepted: 6 November 1996     

Comparison of physical and genetic properties of palindromic DNA sequences.

Some viable palindromic DNA sequences were found to cause an increase in the recovery of genetic recombinants. Although these palindromes contained no Chi sites, their presence in cis caused apparent recA+-dependent recombination to increase severalfold. This biological property did not correlate with the physical properties of the palindromes' extrusion of cruciform structures in vitro. Thus, two unrelated palindromes with similar effects on recombination in both Escherichia coli and Pseudomonas syringae displayed quite different kinetics of cruciform formation. In plasmids of native superhelical density, one palindrome underwent rapid cruciform formation at 55 degrees C, whereas the other did not form detectable cruciforms at any temperature. A shorter palindrome with similarly rapid kinetics of cruciform formation did not affect recombination detectably. The lack of a clear relationship between physical and genetic properties was also demonstrated in the case of longer, inviable palindromes. Here we found that the degree of asymmetry required in vivo to rescue a long palindrome from inviability far exceeded that required to kinetically prohibit cruciform extrusion in vitro.     

GC-biased gene conversion and selection affect GC content in the Oryza genus (rice)

Base composition varies among and within eukaryote genomes. Although mutational bias and selection have initially been invoked, more recently GC-biased gene conversion (gBGC) has been proposed to play a central role in shaping nucleotide landscapes, especially in yeast, mammals, and birds. gBGC is a kind of meiotic drive in favor of G and C alleles, associated with recombination. Previous studies have also suggested that gBGC could be at work in grass genomes. However, these studies were carried on third codon positions that can undergo selection on codon usage. As most preferred codons end in G or C in grasses, gBGC and selection can be confounded. Here we investigated further the forces that might drive GC content evolution in the rice genus using both coding and noncoding sequences. We found that recombination rates correlate positively with equilibrium GC content and that selfing species (Oryza sativa and O. glaberrima) have significantly lower equilibrium GC content compared with more outcrossing species. As recombination is less efficient in selfing species, these results suggest that recombination drives GC content. We also detected a positive relationship between expression levels and GC content in third codon positions, suggesting that selection favors codons ending with G or C bases. However, the correlation between GC content and recombination cannot be explained by selection on codon usage alone as it was also observed in noncoding positions. Finally, analyses of polymorphism data ruled out the hypothesis that genomic variation in GC content is due to mutational processes. Our results suggest that both gBGC and selection on codon usage affect GC content in the Oryza genus and likely in other grass species.     

Non-B right-handed DNA conformations of homopurine.homopyrimidine sequences in the murine immunoglobulin C alpha switch region

The switch region of IgA immunoglobulin in mice cloned into a recombinant plasmid contains a supercoil-dependent S1 nuclease hypersensitive site, indicative of a non-B-DNA secondary structure. This site maps to the (AGGAG)28 direct repeat (DR2) of the alpha switch region and appears at a negative superhelical density of greater than 0.02. Studies with P1 nuclease and bromoacetaldehyde indicate that this structure is also present at neutral pH. S1 nuclease sensitivity is retained for the shorter repeat (AGGAG)6GA in a recombinant plasmid but is not seen for the repeat (CTGAG)6, corresponding to the DR1 repeat of the alpha switch region, or in a sequence corresponding to a portion of the consensus sequence which contains a short stretch of alternating pyrine-pyrimidine residues. Fine mapping of the (AGGAG)6GA and flanking sequences with dimethyl sulfate, bromoacetaldehyde, osmium tetroxide, and diethyl pyrocarbonate reveals an asymmetric pattern of modification dependent on both pH and supercoiling. Two-dimensional gel electrophoresis at low pH shows the relaxation of 3 superhelical turns on formation of this structure by the (AGGAG)6GA repeat. These results are most consistent with the formation of an intramolecular triple-strand.     

The RecE recombination pathway mediates recombination between partially homologous DNA sequences: structural analysis of recombination products

Escherichia coli generalized recombination, utilizing the RecA RecB recombination pathway, requires large stretches (70-200 bp) of complete DNA sequence homology. In contrast, we have found that the RecE pathway can promote recombination between DNA with only short stretches of homology. A plasmid containing 10 partially homologous direct repeats was linearized by digestion with specific restriction enzymes. After transformation, a RecE+ (sbcA) host was able to circularize the plasmid by recombination between partially homologous direct repeat sequences. Recombination occurred in regions of as little as 6 bp of perfect homology. Recombination was enhanced in the regions adjacent to restriction sites used to linearize the plasmid, consistent with a role of double-strand breaks in promoting recombination. A mechanism is proposed in which the 5' exonuclease, ExoVIII, produces 3' single-stranded ends from the linearized plasmid. These pair with other sequences of partial homology. Partial homologies in the sequences flanking the actual join serve to stabilize this recombination intermediate. Recombination is completed by a process of "copy and join." This recombination mechanism requires less homology to stabilize intermediates than the degree of homology needed for mechanisms involving strand invasion. Its role in nature may be to increase genomic diversity, for example, by enhancing recombination between bacteriophages and regions of the bacterial chromosome.     

Silencing Homologous RNA Recombination Hot Spots with GC-Rich Sequences in Brome Mosaic Virus

It has been observed that AU-rich sequences form homologous recombination hot spots in brome mosaic virus (BMV), a tripartite positive-stranded RNA virus of plants (P. D. Nagy and J. J. Bujarski, J. Virol. 71:3799–3810, 1997). To study the effect of GC-rich sequences on the recombination hot spots, we inserted 30-nucleotide-long GC-rich sequences downstream of AU-rich homologous recombination hot spot regions in parental BMV RNAs (RNA2 and RNA3). Although these insertions doubled the length of sequence identity in RNA2 and RNA3, the incidence of homologous RNA2 and RNA3 recombination was reduced markedly. Four different, both highly structured and nonstructured downstream GC-rich sequences had a similar “homologous recombination silencing” effect on the nearby hot spots. The GC-rich sequence-mediated recombination silencing mapped to RNA2, as it was observed when the GC-rich sequence was inserted at downstream locations in both RNA2 and RNA3 or only in the RNA2 component. On the contrary, when the downstream GC-rich sequence was present only in the RNA3 component, it increased the incidence of homologous recombination. In addition, upstream insertions of similar GC-rich sequences increased the incidence of homologous recombination within downstream hot spot regions. Overall, this study reveals the complex nature of homologous recombination in BMV, where sequences flanking the common hot spot regions affect recombination frequency. A replicase-driven template-switching model is presented to explain recombination silencing by GC-rich sequences.     

Flanking sequences affect replication arrest at the Escherichia coli terminator TerB in vivo.

We have analyzed the effect of flanking sequences on Tus-induced replication arrest. pBR322 plasmid derivatives which carry the Escherichia coli replication terminator TerB at different locations were used. Efficiency of the replication arrest was estimated from the plasmid copy number and transformation frequency of tus+ cells. We found that flanking sequences do affect replication arrest efficiency, a weak arrest being correlated with the presence of an AT-rich region which is replicated just before TerB. Some sequences located after the replication terminator can also affect replication termination. We propose that the AT-rich regions might impair binding of the Tus protein to the TerB sequence or facilitate helicase-induced unwinding of DNA and Tus displacement from the TerB site.     

Engineering of homologous recombination hotspots with AU-rich sequences in brome mosaic virus.

Previously, we observed that crossovers sites of RNA recombinants clustered within or close to AU-rich regions during genetic recombination in brome mosaic bromovirus (BMV) (P. D. Nagy and J. J. Bujarski. J. Virol. 70:415-426, 1996). To test whether AU-rich sequences can facilitate homologous recombination, AU-rich sequences were introduced into parental BMV RNAs (RNA2 and RNA3). These insertions created a homologous RNA2-RNA3 recombination hotspot. Two other AU-rich sequences also supported high-frequency homologous recombination if a common sequence with high or average G/C content was present immediately upstream of the AU-rich element. Homologous RNA recombination did not require any additional sequence motifs or RNA structures and was position nonspecific within the 3' noncoding region. These results suggest that nucleotide content (i.e., the presence of common 5' GC-rich or moderately AU-rich and 3' AU-rich regions) is the important factor that determines the sites of homologous recombination. A mechanism that involves replicase switching during synthesis of positive-sense RNA strands is presented to explain the observed results.     

DNA structural polymorphism modulates the kinetics of superhelical DNA cleavage by BamHI restriction endonuclease

A compartmental model developed by Hensley (Hensley, P., Nardone, G., Chirikjian, J.G., and Wastney, M. E., (1990) J. Biol. Chem. 265, 15300-15307) for analysis of the time courses of the cleavage of superhelical DNA substrates by the restriction endonuclease, BamHI, has been used to quantify the effects of changes in temperature, ionic strength, superhelical density, and the DNA substrate on the binding and strand cleavage processes. Studies reported here indicate that changes in topology may be introduced into the DNA substrate solely as a result of the plasmid preparation process and in the absence of covalent bond cleavage and ligation. These changes in topology have qualitatively different effects on the kinetics than those promoted by changes in the superhelical density. The former are removed by briefly warming the DNA prior to assay, suggesting that they are only kinetically stable, while the latter changes are not affected by heating. Increasing the [NaCl] from 0.01 M to 0.1 M increases the overall rate of plasmid cleavage by increasing both the rates of cleavage and enzyme DNA association. To describe the decrease in the overall cleavage rate observed in 0.15 M NaCl, an ionic strength-dependent rate-determining structural transition in the DNA substrate was incorporated into the model. The largest changes in the rate of the cleavage process resulted from changes in the DNA substrate. For the SV40 substrate compared to pBR322, the rate constants describing the two association processes and the first bond cleavage event were increased 6- to 7-fold. The rate of the second bond cleavage process was not affected. These changes may be due to differences in the flanking sequences.