d(GA x TC)(n) microsatellite DNA sequences enhance homologous DNA recombination in SV40 minichromosomes

MicrosatNavigator: exploring nonrandom distribution and lineage-specificity of microsatellite repeat motifs on vertebrate sex chromosomes across 186 whole genomes

Microsatellites are short tandem DNA repeats, ubiquitous in genomes. They are believed to be under selection pressure, considering their high distribution and abundance beyond chance or random accumulation. However, limited analysis of microsatellites in single taxonomic groups makes it challenging to understand their evolutionary significance across taxonomic boundaries. Despite abundant genomic information, microsatellites have been studied in limited contexts and within a few species, warranting an unbiased examination of their genome-wide distribution in distinct versus closely related-clades. Large-scale comparisons have revealed relevant trends, especially in vertebrates. Here, "MicrosatNavigator", a new tool that allows quick and reliable investigation of perfect microsatellites in DNA sequences, was developed. This tool can identify microsatellites across the entire genome sequences. Using this tool, microsatellite repeat motifs were identified in the genome sequences of 186 vertebrates. A significant positive correlation was noted between the abundance, density, length, and GC bias of microsatellites and specific lineages. The (AC)n motif is the most prevalent in vertebrate genomes, showing distinct patterns in closely related species. Longer microsatellites were observed on sex chromosomes in birds and mammals but not on autosomes. Microsatellites on sex chromosomes of non-fish vertebrates have the lowest GC content, whereas high-GC microsatellites (≥ 50 M% GC) are preferred in bony and cartilaginous fishes. Thus, similar selective forces and mutational processes may constrain GC-rich microsatellites to different clades. These findings should facilitate investigations into the roles of microsatellites in sex chromosome differentiation and provide candidate microsatellites for functional analysis across the vertebrate evolutionary spectrum.

Homopurine guanine-rich sequences in complex with N-methyl mesoporphyrin IX form parallel G-quadruplex dimers and display a unique symmetry tetrad

DNA can fold into G-quadruplexes (GQs), non-canonical secondary structures formed by π-π stacking of G-tetrads. GQs are important in many biological processes, which makes them promising therapeutic targets. We identified a 42-nucleotide long, purine-only G-rich sequence from human genome, which contains eight G-stretches connected by A and AAAA loops. We divided this sequence into five unique segments, four guanine stretches each, named GA1-5. In order to investigate the role of adenines in GQ structure formation, we performed biophysical and X-ray crystallographic studies of GA1-5 and their complexes with a highly selective GQ ligand, N-methyl mesoporphyrin IX (NMM). Our data indicate that all variants form parallel GQs whose stability depends on the number of flexible AAAA loops. GA1-3 bind NMM with 1:1 stoichiometry. The Ka for GA1 and GA3 is modest, ∼0.3 μM -1, and that for GA2 is significantly higher, ∼1.2 μM -1. NMM stabilizes GA1-3 by 14.6, 13.1, and 7.0 °C, respectively, at 2 equivalents. We determined X-ray crystal structures of GA1-NMM (1.98 Å resolution) and GA3-NMM (2.01 Å). The structures confirm the parallel topology of GQs with all adenines forming loops and display NMM binding at the 3' G-tetrad. Both complexes dimerize through the 5' interface. We observe two novel structural features: 1) a 'symmetry tetrad' at the dimer interface, which is formed by two guanines from each GQ monomer and 2) a NMM dimer in GA1-NMM. Our structural work confirms great flexibility of adenines as structural elements in GQ formation and contributes greatly to our understanding of the structural diversity of GQs and their modes of interaction with small molecule ligands.

Revealing structural peculiarities of homopurine GA repetition stuck by i-motif clip

Non-canonical forms of nucleic acids represent challenging objects for both structure-determination and investigation of their potential role in living systems. In this work, we uncover a structure adopted by GA repetition locked in a parallel homoduplex by an i-motif. A series of DNA oligonucleotides comprising GAGA segment and C3 clip is analyzed by NMR and CD spectroscopies to understand the sequence-structure-stability relationships. We demonstrate how the relative position of the homopurine GAGA segment and the C3 clip as well as single-base mutations (guanine deamination and cytosine methylation) affect base pairing arrangement of purines, i-motif topology and overall stability. We focus on oligonucleotides C3GAGA and methylated GAGAC3 exhibiting the highest stability and structural uniformity which allowed determination of high-resolution structures further analyzed by unbiased molecular dynamics simulation. We describe sequence-specific supramolecular interactions on the junction between homoduplex and i-motif blocks that contribute to the overall stability of the structures. The results show that the distinct structural motifs can not only coexist in the tight neighborhood within the same molecule but even mutually support their formation. Our findings are expected to have general validity and could serve as guides in future structure and stability investigations of nucleic acids.

Microsatellite landscape evolutionary dynamics across 450 million years of vertebrate genome evolution

The evolutionary dynamics of simple sequence repeats (SSRs or microsatellites) across the vertebrate tree of life remain largely undocumented and poorly understood. In this study, we analyzed patterns of genomic microsatellite abundance and evolution across 71 vertebrate genomes. The highest abundances of microsatellites exist in the genomes of ray-finned fishes, squamate reptiles, and mammals, while crocodilian, turtle, and avian genomes exhibit reduced microsatellite landscapes. We used comparative methods to infer evolutionary rates of change in microsatellite abundance across vertebrates and to highlight particular lineages that have experienced unusually high or low rates of change in genomic microsatellite abundance. Overall, most variation in microsatellite content, abundance, and evolutionary rate is observed among major lineages of reptiles, yet we found that several deeply divergent clades (i.e., squamate reptiles and mammals) contained relatively similar genomic microsatellite compositions. Archosauromorph reptiles (turtles, crocodilians, and birds) exhibit reduced genomic microsatellite content and the slowest rates of microsatellite evolution, in contrast to squamate reptile genomes that have among the highest rates of microsatellite evolution. Substantial branch-specific shifts in SSR content in primates, monotremes, rodents, snakes, and fish are also evident. Collectively, our results support multiple major shifts in microsatellite genomic landscapes among vertebrates.

Genomic correlates of recombination rate and its variability across eight recombination maps in the western honey bee (Apis mellifera L.)

Background

Meiotic recombination has traditionally been explained based on the structural requirement to stabilize homologous chromosome pairs to ensure their proper meiotic segregation. Competing hypotheses seek to explain the emerging findings of significant heterogeneity in recombination rates within and between genomes, but intraspecific comparisons of genome-wide recombination patterns are rare. The honey bee (Apis mellifera) exhibits the highest rate of genomic recombination among multicellular animals with about five cross-over events per chromatid.

Results

Here, we present a comparative analysis of recombination rates across eight genetic linkage maps of the honey bee genome to investigate which genomic sequence features are correlated with recombination rate and with its variation across the eight data sets, ranging in average marker spacing ranging from 1 Mbp to 120 kbp. Overall, we found that GC content explained best the variation in local recombination rate along chromosomes at the analyzed 100 kbp scale. In contrast, variation among the different maps was correlated to the abundance of microsatellites and several specific tri- and tetra-nucleotides.

Conclusions

The combined evidence from eight medium-scale recombination maps of the honey bee genome suggests that recombination rate variation in this highly recombining genome might be due to the DNA configuration instead of distinct sequence motifs. However, more fine-scale analyses are needed. The empirical basis of eight differing genetic maps allowed for robust conclusions about the correlates of the local recombination rates and enabled the study of the relation between DNA features and variability in local recombination rates, which is particularly relevant in the honey bee genome with its exceptionally high recombination rate.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1281-2) contains supplementary material, which is available to authorized users.

A genome-wide map of adeno-associated virus-mediated human gene targeting

To determine which genomic features promote homologous recombination, we created a genome-wide map of gene targeting sites. We used an adeno-associated virus vector to target identical loci introduced as transcriptionally active retroviral vectors. A comparison of ~2,000 targeted and untargeted sites showed that targeting occurred throughout the human genome and was not influenced by the presence of nearby CpG islands, sequence repeats or DNase I-hypersensitive sites. Targeted sites were preferentially located within transcription units, especially when the target loci were transcribed in the opposite orientation to their surrounding chromosomal genes. We determined the impact of DNA replication by mapping replication forks, which revealed a preference for recombination at target loci transcribed toward an incoming fork. Our results constitute the first genome-wide screen of gene targeting in mammalian cells and demonstrate a strong recombinogenic effect of colliding polymerases.

Genome wide survey and analysis of small repetitive sequences in caulimoviruses

Microsatellites are known to exhibit ubiquitous presence across all kingdoms of life including viruses. Members of the Caulimoviridae family severely affect growth of vegetable and fruit plants and reduce economic yield in diverse cropping systems worldwide. Here, we analyzed the nature and distribution of both simple and complex microsatellites present in complete genome of 44 species of Caulimoviridae. Our results showed, in all analyzed genomes, genome size and GC content had a weak influence on number, relative abundance and relative density of microsatellites, respectively. For each genome, mono- and dinucleotide repeats were found to be highly predominant and are overrepresented in genome of majority of caulimoviruses. AT/TA and GAA/AAG/AGA was the most abundant di- and trinucleotide repeat motif, respectively. Repeats larger than trinucleotide were rarely found in these genomes. Comparative study of occurrence, abundance and density of microsatellite among available RNA and DNA viral genomes indicated that simple repeats were least abundant in genomes of caulimoviruses. Polymorphic repeats even though rare were observed in the large intergenic region of the genome, indicating strand slippage and/or unequal recombination processes do occur in caulimoviruses. To our knowledge, this is the first analysis of microsatellites occurring in any dsDNA viral genome. Characterization of such variations in repeat sequences would be important in deciphering the origin, mutational processes, and role of repeat sequences in viral genomes.

Distinct mutational behaviors differentiate short tandem repeats from microsatellites in the human genome

A tandem repeat's (TR) propensity to mutate increases with repeat number, and can become very pronounced beyond a critical boundary, transforming it into a microsatellite (MS). However, a clear understanding of the mutational behavior of different TR classes and motifs and related mechanisms is lacking, as is a consensus on the existence of a boundary separating short TRs (STRs) from MSs. This hinders our understanding of MSs' mutational properties and their effective use as genetic markers. Using indel calls for 179 individuals from 1000 Genomes Pilot-1 Project, we determined polymorphism incidence for four major TR classes, and formalized its varying relationship with repeat number using segmented regression. We observed a biphasic regime with a transition from a faster to a slower exponential growth at 9, 5, 4, and 4 repeats for mono-, di-, tri-, and tetranucleotide TRs, respectively. We used an in vitro mutagenesis assay to evaluate the contribution of strand slippage errors to mutability. STRs and MSs differ in their absolute polymorphism levels, but more importantly in their rates of mutability growth. Although strand slippage is a major factor driving mononucleotide polymorphism incidence, dinucleotide polymorphism incidence is greater than that expected due to strand slippage alone, indicating that additional cellular factors might be driving dinucleotide mutability in the human genome. Leveraging on hundreds of human genomes, we present the first comprehensive, genome-wide analysis of TR mutational behavior, encompassing several motif sizes and compositions.

Chromosomal position effects on AAV-mediated gene targeting

The effects of chromosomal position and neighboring genomic elements on gene targeting in human cells remain largely unexplored. To study these, we used a shuttle vector system in which murine leukemia virus (MLV)-based proviral targets present at different chromosomal locations and containing mutations in the neomycin phosphotransferase (neo) gene were corrected by adeno-associated virus (AAV)-mediated gene targeting. Sixteen identical target loci present in HT-1080 human sarcoma cells were all successfully corrected by gene targeting. The gene targeting frequencies varied by as much as 10-fold, and there was a clear bias for correction of one of the targets in clones containing two target sites. The targeting frequency at each site was correlated to the proximity and density of various genomic elements, and we found a significant association of higher targeting frequencies at loci near a subset of dinucleotide microsatellite repeats (r = –0.55, P < 0.05), in particular GT repeats (r = –0.87, P < 0.0001). Additionally, there was a correlation between meiotic recombination rates and targeting frequencies at the target loci (r = 0.52, P < 0.05). There was no correlation between surrounding chromosomal transcription units and targeting frequencies. Our results indicate that certain chromosomal positions are preferred sites for gene targeting in human cells.

Detecting microsatellites within genomes: significant variation among algorithms

Background

Microsatellites are short, tandemly-repeated DNA sequences which are widely distributed among genomes. Their structure, role and evolution can be analyzed based on exhaustive extraction from sequenced genomes. Several dedicated algorithms have been developed for this purpose. Here, we compared the detection efficiency of five of them (TRF, Mreps, Sputnik, STAR, and RepeatMasker).

Results

Our analysis was first conducted on the human X chromosome, and microsatellite distributions were characterized by microsatellite number, length, and divergence from a pure motif. The algorithms work with user-defined parameters, and we demonstrate that the parameter values chosen can strongly influence microsatellite distributions. The five algorithms were then compared by fixing parameters settings, and the analysis was extended to three other genomes (Saccharomyces cerevisiae, Neurospora crassa and Drosophila melanogaster) spanning a wide range of size and structure. Significant differences for all characteristics of microsatellites were observed among algorithms, but not among genomes, for both perfect and imperfect microsatellites. Striking differences were detected for short microsatellites (below 20 bp), regardless of motif.

Conclusion

Since the algorithm used strongly influences empirical distributions, studies analyzing microsatellite evolution based on a comparison between empirical and theoretical size distributions should therefore be considered with caution. We also discuss why a typological definition of microsatellites limits our capacity to capture their genomic distributions.

Studies of DNA dumbbells VIII. Melting analysis of DNA dumbbells with dinucleotide repeat stem sequences

Melting curves and circular dichroism spectra were measured for a number of DNA dumbbell and linear molecules containing dinucleotide repeat sequences of different lengths. To study effects of different sequences on the melting and spectroscopic properties, six DNA dumbbells whose stems contain the central sequences (AA)(10), (AC)(10), (AG)(10), (AT)(10), (GC)(10), and (GG)(10) were prepared. These represent the minimal set of 10 possible dinucleotide repeats. To study effects of dinucleotide repeat length, dumbbells with the central sequences (AG)(n), n = 5 and 20, were prepared. Control molecules, dumbbells with a random central sequence, (RN)(n), n = 5, 10, and 20, were also prepared. The central sequence of each dumbbell was flanked on both sides by the same 12 base pairs and T(4) end-loops. Melting curves were measured by optical absorbance and differential scanning calorimetry in solvents containing 25, 55, 85, and 115 mM Na(+). CD spectra were collected from 20 to 45 degrees C and [Na(+)] from 25 to 115 mM. The spectral database did not reveal any apparent temperature dependence in the pretransition region. Analysis of the melting thermodynamics evaluated as a function of Na(+) provided a means for quantitatively estimating the counterion release with melting for the different sequences. Results show a very definite sequence dependence, indicating the salt-dependent properties of duplex DNA are also sequence dependent. Linear DNA molecules containing the (AG)(n) and (RN)(n), sequences, n = 5, 10, 20, and 30, were also prepared and studied. The linear DNA molecules had the exact sequences of the dumbbell stems. That is, the central repeat sequence in each linear duplex was flanked on both sides by the same 12-bp sequence. Melting and CD studies were also performed on the linear DNA molecules. Comparison of results obtained for the same sequences in dumbbell and linear molecular environments reveals several interesting features of the interplay between sequence-dependent structural variability, sequence length, and the unconstrained (linear) or constrained (dumbbell) molecular environments.

Advances in mechanisms of genetic instability related to hereditary neurological diseases

Substantial progress has been realized in the past several years in our understanding of the molecular mechanisms responsible for the expansions and deletions (genetic instabilities) of repeating tri-, tetra- and pentanucleotide repeating sequences associated with a number of hereditary neurological diseases. These instabilities occur by replication, recombination and repair processes, probably acting in concert, due to slippage of the DNA complementary strands relative to each other. The biophysical properties of the folded-back repeating sequence strands play a critical role in these instabilities. Non-B DNA structural elements (hairpins and slipped structures, DNA unwinding elements, tetraplexes, triplexes and sticky DNA) are described. The replication mechanisms are influenced by pausing of the replication fork, orientation of the repeat strands, location of the repeat sequences relative to replication origins and the flap endonuclease. Methyl-directed mismatch repair, nucleotide excision repair, and repair of damage caused by mutagens are discussed. Genetic recombination and double-strand break repair advances in Escherichia coli, yeast and mammalian models are reviewed. Furthermore, the newly discovered capacities of certain triplet repeat sequences to cause gross chromosomal rearrangements are discussed.

Genomic context of paralogous recombination hotspots mediating recurrent NF1 region microdeletion

Recombination between paralogs that flank the NF1 gene at 17q11.2 typically results in a 1.5-Mb microdeletion that includes NF1 and at least 13 other genes. We show that the principal sequences responsible are two 51-kb blocks with 97.5% sequence identity (NF1REP-P1-51 and NF1REP-M-51). These blocks belong to a complex group of paralogs with three components on 17q11.2 and another on 19p13.13. Breakpoint sequencing of deleted chromosomes from multiple patients revealed two paralogous recombination hot spots within the 51-kb blocks. Lack of sequence similarity between these sites failed to suggest or corroborate any putative cis-acting recombinogenic motifs. However, the NF1 REPs showed relatively high alignment mismatch between recombining paralogs, and we note that the NF1REP hot spots were regions of good alignment bordered by relatively large alignment gaps. Statistical tests for gene conversion detected a single significant tract of perfect match between the NF1REPs that was 700 bp long and coincided with PRS2, the predominant recombination hot spot. Tracts of perfect match occurring by chance may contribute to breakpoint localization, but our result suggests that perfect tracts at recombination hot spots may be a result of gene conversion at sites at which preferential pairing occurs for other, as-yet-unknown reasons.

Structure-dependent Recombination Hot Spot Activity of GAA·TTC Sequences from Intron 1 of the Friedreich's Ataxia Gene

The recombinational properties of long GAA.TTC repeating sequences were analyzed in Escherichia coli to gain further insights into the molecular mechanisms of the genetic instability of this tract as possibly related to the etiology of Friedreich's ataxia. Intramolecular and intermolecular recombination studies showed that the frequency of recombination between the GAA.TTC tracts was as much as 15 times higher than the non-repeating control sequences. Homologous, intramolecular recombination between GAA.TTC tracts and GAAGGA.TCCTTC repeats also occurred with a very high frequency (approximately 0.8%). Biochemical analyses of the recombination products demonstrated the expansions and deletions of the GAA.TTC repeats. These results, together with our previous studies on the CTG.CAG sequences, suggest that the recombinational hot spot characteristics may be a common feature of all triplet repeat sequences. Unexpectedly, we found that the recombination properties of the GAA.TTC tracts were unique, compared with CTG.CAG repeats, because they depended on the DNA secondary structure polymorphism. Increasing the length of the GAA.TTC repeats decreased the intramolecular recombination frequency between these tracts. Also, a correlation was found between the propensity of the GAA.TTC tracts to adopt the sticky DNA conformation and the inhibition of intramolecular recombination. The use of novobiocin to modulate the intracellular DNA topology, i.e. the lowering of the negative superhelical density, repressed the formation of the sticky DNA structure, thereby restoring the expected positive correlation between the length of the GAA.TTC tracts and the frequency of intramolecular recombination. Hence, our results demonstrate that sticky DNA exists and functions in E. coli.

Circular dichroism spectroscopy of conformers of (guanine + adenine) repeat strands of DNA

(Guanine+adenine) strands of DNA are known to associate into guanine tetraplexes, homodimerize into parallel or antiparallel duplexes, and fold into a cooperatively melting single strand resembling the protein alpha helix. Using CD spectroscopy and other methods, we studied how this conformational polymorphism depended on the primary structure of DNA. The study showed that d(GGGA)(5) and d(GGA)(7) associated into homoduplexes at low salt or in the presence of LiCl but were prone to guanine tetraplex formation, especially in the presence of KCl. In addition, they yielded essentially the same CD spectrum in the presence of ethanol as observed with the ordered single strand of d(GA)(10). Strands of d(GA)(10), d(GGAA)(5), d(GAA)(7), and d(GAAA)(5) associated into homoduplexes in both LiCl and KCl solutions, but not into guanine tetraplexes. d(GAAA)(5) and d(GAA)(7) further failed to form the single-stranded conformer in aqueous ethanol. Adenine protonation, however, stabilized the single-stranded conformer even in these adenine-rich fragments. The ordered single strands, homoduplexes as well as the guanine tetraplexes, all provided strikingly similar CD spectra, indicating that all of the conformers shared similar base stacking geometries. The increasing adenine content only decreased the conformer thermostability.

Long CTG.CAG repeats from myotonic dystrophy are preferred sites for intermolecular recombination

Homologous recombination was shown to enable the expansion of CTG.CAG repeat sequences. Other prior investigations revealed the involvement of replication and DNA repair in these genetic instabilities. Here we used a genetic assay to measure the frequency of homologous intermolecular recombination between two CTG.CAG tracts. When compared with non-repeating sequences of similar lengths, long (CTG.CAG)(n) repeats apparently recombine with an approximately 60-fold higher frequency. Sequence polymorphisms that interrupt the homogeneity of the CTG.CAG repeat tracts reduce the apparent recombination frequency as compared with the pure uninterrupted repeats. The orientation of the repeats relative to the origin of replication strongly influenced the apparent frequency of recombination. This suggests the involvement of DNA replication in the recombination process of triplet repeats. We propose that DNA polymerases stall within the CTG.CAG repeat tracts causing nicks or double-strand breaks that stimulate homologous recombination. The recombination process is RecA-dependent.

Structural heterogeneity of pyrimidine/purine-biased DNA sequence analyzed by atomic force microscopy

We report here the direct evidence for the formation of alternative DNA structures in a plasmid DNA, termed pTIR10, containing a 0.23-kb pyrimidine/purine-biased (Pyr/Pur) stretch isolated from the rat genome. Long Pyr/Pur sequences are abundant in eukaryotic genomes, and they may modulate the biological activity of genes and genomes via formation of various types of triplex-related structures. The plasmid DNA in sodium acetate buffer (pH 4.35) was deposited on APS-modified mica, and after drying it was imaged with an atomic force microscope in air. Various types of thick protrusions have been observed on pTIR10 DNA. Structural parameters (width and height) of DNA molecules suggest that the alternative structures observed here are variations on the theme of an intramolecular triplex. The biological relevance of the structural features within Pyr/Pur stretches is discussed.

Interaction in vitro of type III intermediate filament proteins with triplex DNA

As previously shown, type III intermediate filaments (IFs) select from a mixture of linear mouse genomic DNA fragments mobile and repetitive, recombinogenic sequences that have also been identified in SDS-stable crosslinkage products of vimentin and DNA isolated from intact fibroblasts. Because these sequences also included homopurine.homopyrimidine (Pu.Py) tracts known to adopt triple-helical conformation under superhelical tension, and because IF proteins are single-stranded (ss) and supercoiled DNA-binding proteins, it was of interest whether they have a particular affinity for triplex DNA. To substantiate this, IF-selected DNA fragments harboring a (Pu.Py) segment and synthetic d(GA)(n) microsatellites were inserted into a vector plasmid and the constructs analyzed for their capacity to interact with IF proteins. Band shift assays revealed a substantially higher affinity of the IF proteins for the insert-containing plasmids than for the empty vector, with an activity decreasing in the order of vimentin > glial fibrillary acidic protein > desmin. In addition, footprint analyses performed with S1 nuclease, KMnO(4), and OsO(4)/bipyridine showed that the (Pu.Py) inserts had adopted triplex conformation under the superhelical strain of the plasmids, and that the IF proteins protected the triple-helical insert sequences from nucleolytic cleavage and chemical modification. All these activities were largely reduced in extent when analyzed on linearized plasmid DNAs. Because intramolecular triplexes (H-DNA) expose single-stranded loops, and the prokaryotic ssDNA-binding proteins g5p and g32p also protected at least the Pu-strand of the (Pu.Py) inserts from nucleolytic degradation, it seemed likely that the IF proteins take advantage of their ssDNA-binding activity in interacting with H-DNA. However, in contrast to g5p and E. coli SSB, they produced no clear band shifts with single-stranded d(GA)(20) and d(TC)(20), so that the interactions rather appear to occur via the duplex-triplex and triplex-loop junctions of H-DNA. On the other hand, the IF proteins, and also g32p, promoted the formation of intermolecular triplexes from the duplex d[A(GA)(20).(TC)(20)T] and d(GA)(20) and d(TC)(20) single strands, with preference of the Py (Pu.Py) triplex motif, substantiating an affinity of the proteins for the triplex structure as such. This triplex-stabilizing effect of IF proteins also applies to the H-DNA of (Pu.Py) insert-containing plasmids, as demonstrated by the preservation of intramolecular triplex-vimentin complexes upon linearization of their constituent supercoiled DNAs, in contrast to poor complex formation from free, linearized plasmid DNA and vimentin. Considering that (Pu.Py) sequences are found near MAR/replication origins, in upstream enhancer and promoter regions of genes, and in recombination hot spots, these results might point to roles of IF proteins in DNA replication, transcription, recombination, and repair.