TRACTS: A program to map oligopurine.oligopyrimidine and other binary DNA tracts

A program to map the locations and frequencies of DNA tracts composed of only two bases ('Binary DNA') is described. The program, TRACTS (URL http://bioportal.weizmann.ac.il/tracts/tracts.html and/or http://bip.weizmann.ac.il/miwbin/servers/tracts) is of interest because long tracts composed of only two bases are highly over-represented in most genomes. In eukaryotes, oligopurine.oligopyrimidine tracts ('R.Y tracts') are found in the highest excess. In prokaryotes, W tracts predominate (A,T 'rich'). A pre-program, ANEX, parses database annotation files of GenBank and EMBL, to produce a convenient one-line list of every gene (exon, intron) in a genome. The main unit lists and analyzes tracts of the three possible binary pairs (R.Y, K.M and S;W). As an example, the results of R.Y tract mapping of mammalian gene p53 is described.     

Long W tracts are over-represented in the Escherichia coli and Haemophilus influenzae genomes.

The occurrence of DNA tracts of the three binary base combinations: R.Y, K.M and W;S has been mapped in the complete genomes of Haemophilus influenzae and Escherichia coli. A highly significant over-representation of W tracts is observed in both bacteria. The excess of W tracts is particularly striking in the 10% intercoding regions. Subdivision of intercoding regions into divergent (promoting), convergent (terminating) and sequential subregions shows that the excess of W tracts is most concentrated in the promoter regions. A particularly high excess of W tracts is observed in the first 200 bases 5' upstream of coding start sites. The data suggest that W tracts have a role in promoter function. A function as unwinding centers, analogous to the role of R.Y tracts in eukaryotes, is proposed. R.Y and K.M tracts are only modestly over-represented in the two bacteria.     

Genome-wide comparative analysis of putative bidirectional promoters from rice, Arabidopsis and Populus

A bidirectional promoter can regulate the expression of two flanking genes arranged in a divergent manner. Although reports pertaining to bidirectional promoters on a genomic scale exist in mammals, little progress has been made in plants. In the present study, we performed a computational analysis of this unique class of promoters to identify overrepresented cis-regulatory motifs from three sequenced plant genomes: rice (Oryza sativa), Arabidopsis thaliana, and Populus trichocarpa using the Plant Cis-acting Regulatory DNA Elements (PLACE) and PLANT CARE databases. We describe these overrepresented elements and their possible regulatory mechanisms. We also discuss similarities and differences with human bidirectional promoters. Furthermore, we describe in detail a few coexpressed and evolutionarily conserved divergent gene pairs and their bidirectional promoters. This study provides insights into bidirectional promoters in three plant species, thereby laying a foundation for their experimental analysis.     

The frequency of two-base tracts in eukaryotic genomes

The frequency of two-base tracts is surveyed in a wide range of eukaryotic genomes using the special program TRACTS. All three two-base families are surveyed: R.Y (A,G.C,T), K.M (A,C.G,T), and S;W (A.T and G.C). Data for the human β-globin complex, for the tobacco chloroplast, and for 247 nt mammalian promoter regions are presented. All two-base tracts longer than three or four bases are overrepresented to an extent surpassing by far their occurrence in a randomized DNA population in the majority of the genomic regions analyzed; 20–30 long tracts are quite frequent, against the statistical odds. R.Y tracts are found at the largest excess, K.M tract to a slightly lesser extent, while S.W tracts are found at a moderate yet significant excess. The majority of the tracts manifest only a limited extent of tandem repeat structures. The idea that the two base tracts serve as unwinding elements is considered. Preseented at the NATO Advanced Research Workshop onGenome Organization and Evolution, Spetsai, Greece, 16–22 September 1992     

The clustering of CpG islands may constitute an important determinant of the 3D organization of interphase chromosomes

We used the 4C-Seq technique to characterize the genome-wide patterns of spatial contacts of several CpG islands located on chromosome 14 in cultured chicken lymphoid and erythroid cells. We observed a clear tendency for the spatial clustering of CpG islands present on the same and different chromosomes, regardless of the presence or absence of promoters within these CpG islands. Accordingly, we observed preferential spatial contacts between Sp1 binding motifs and other GC-rich genomic elements, including the DNA sequence motifs capable of forming G-quadruplexes. However, an anchor placed in a gene/CpG island-poor area formed spatial contacts with other gene/CpG island-poor areas on chromosome 14 and other chromosomes. These results corroborate the two-compartment model of the spatial organization of interphase chromosomes and suggest that the clustering of CpG islands constitutes an important determinant of the 3D organization of the eukaryotic genome in the cell nucleus. Using the ChIP-Seq technique, we mapped the genome-wide CTCF deposition sites in the chicken lymphoid and erythroid cells that were used for the 4C analysis. We observed a good correlation between the density of CTCF deposition sites and the level of 4C signals for the anchors located in CpG islands but not for an anchor located in a gene desert. It is thus possible that CTCF contributes to the clustering of CpG islands observed in our experiments.     

The 5S rDNA related repetitive sequences in the sex chromosomes of the spiny eel (Mastacembelus aculeatus)

The karyotype of the spiny eel (Mastacembelus aculeatus) has highly evolved heteromorphic sex chromosomes. X and Y chromosomes differ from each other in the distribution of heterochromatin blocks. To characterize the repetitive sequences in these heterochromatic regions, we microdissected the X chromosome, constructed an X chromosome library, amplified the genomic DNA using PCR and isolated a repetitive sequence DNA family by screening the library. All family members were clusters of two simple repetitive monomers, MaSRS1 and MaSRS2. We detected a conserved 5S rDNA gene sequence within monomer MaSRS2; thus, tandem-arranged MaSRS1s and MaSRS2s may co-compose 5S rDNA multigenes and NTSs in M. aculeatus. FISH analysis revealed that MaSRS1 and MaSRS2were the main components of the heterochromatic regions of the X and Y chromosomes. This finding contributes additional data about differentiation of heteromorphic sex chromosomes in lower vertebrates.     

The clustering of CpG islands may constitute an important determinant of the 3D organization of interphase chromosomes

We used the 4C-Seq technique to characterize the genome-wide patterns of spatial contacts of several CpG islands located on chromosome 14 in cultured chicken lymphoid and erythroid cells. We observed a clear tendency for the spatial clustering of CpG islands present on the same and different chromosomes, regardless of the presence or absence of promoters within these CpG islands. Accordingly, we observed preferential spatial contacts between Sp1 binding motifs and other GC-rich genomic elements, including the DNA sequence motifs capable of forming G-quadruplexes. However, an anchor placed in a gene/CpG island-poor area formed spatial contacts with other gene/CpG island-poor areas on chromosome 14 and other chromosomes. These results corroborate the two-compartment model of the spatial organization of interphase chromosomes and suggest that the clustering of CpG islands constitutes an important determinant of the 3D organization of the eukaryotic genome in the cell nucleus. Using the ChIP-Seq technique, we mapped the genome-wide CTCF deposition sites in the chicken lymphoid and erythroid cells that were used for the 4C analysis. We observed a good correlation between the density of CTCF deposition sites and the level of 4C signals for the anchors located in CpG islands but not for an anchor located in a gene desert. It is thus possible that CTCF contributes to the clustering of CpG islands observed in our experiments.     

Generation of Five High-Complexity Painting Probe Libraries from Flow-Sorted Mouse Chromosomes

Mouse metaphase chromosomes were purified by flow sorting from the murine fibroblast cell line Mus spretus clone 5A. We sorted chromosomes that fell into five individual peaks based on the Hoechst 33258/chromomycin A3 DNA histogram: three peaks corresponding to the least amount of DNA and two peaks representing chromosomes with the most DNA content. This is the first example of the successful application of bivariate flow karyotyping to murine chromosome sorting. We then applied primer-directed in vitro DNA amplification using the polymerase chain reaction (PCR) to generate and label larger amounts of chromosome-specific DNA. In situ hybridization showed specific binding of the PCR products to mouse chromosomes Y, 19, 18, 3, and X as well as chromosomes 1 and 2. The combination of chromosome sorting from the M. spretus cell line and PCR proved to be highly valuable for generation of pools of DNA fragments that exhibit specific binding to mouse chromosomes and can be used to identify and delineate mouse metaphase chromosomes.     

Epigenetic control of meiotic recombination in plants

Meiotic recombination is a deeply conserved process within eukaryotes that has a profound effect on patterns of natural genetic variation. During meiosis homologous chromosomes pair and undergo DNA double strand breaks generated by the Spo11 endonuclease. These breaks can be repaired as crossovers that result in reciprocal exchange between chromosomes. The frequency of recombination along chromosomes is highly variable, for example, crossovers are rarely observed in heterochromatin and the centromeric regions. Recent work in plants has shown that crossover hotspots occur in gene promoters and are associated with specific chromatin modifications, including H2 A.Z. Meiotic chromosomes are also organized in loop-base arrays connected to an underlying chromosome axis, which likely interacts with chromatin to organize patterns of recombination.Therefore, epigenetic information exerts a major influence on patterns of meiotic recombination along chromosomes, genetic variation within populations and evolution of plant genomes.     

Meiotic recombination in Drosophila Msh6 mutants yields discontinuous gene conversion tracts

Crossovers (COs) generated through meiotic recombination are important for the correct segregation of homologous chromosomes during meiosis. Several models describing the molecular mechanism of meiotic recombination have been proposed. These models differ in the arrangement of heteroduplex DNA (hDNA) in recombination intermediates. Heterologies in hDNA are usually repaired prior to the recovery of recombination products, thereby obscuring information about the arrangement of hDNA. To examine hDNA in meiotic recombination in Drosophila melanogaster, we sought to block hDNA repair by conducting recombination assays in a mutant defective in mismatch repair (MMR). We generated mutations in the MMR gene Msh6 and analyzed recombination between highly polymorphic homologous chromosomes. We found that hDNA often goes unrepaired during meiotic recombination in an Msh6 mutant, leading to high levels of postmeiotic segregation; however, hDNA and gene conversion tracts are frequently discontinuous, with multiple transitions between gene conversion, restoration, and unrepaired hDNA. We suggest that these discontinuities reflect the activity of a short-patch repair system that operates when canonical MMR is defective.     

Defining the DNA uptake specificity of naturally competent Haemophilus influenzae cells

Some naturally competent bacteria exhibit both a strong preference for DNA fragments containing specific ‘uptake sequences’ and dramatic overrepresentation of these sequences in their genomes. Uptake sequences are often assumed to directly reflect the specificity of the DNA uptake machinery, but the actual specificity has not been well characterized for any bacterium. We produced a detailed analysis of Haemophilus influenzae’s uptake specificity, using Illumina sequencing of degenerate uptake sequences in fragments recovered from competent cells. This identified an uptake motif with the same consensus as the motif overrepresented in the genome, with a 9 bp core (AAGTGCGGT) and two short flanking T-rich tracts. Only four core bases (GCGG) were critical for uptake, suggesting that these make strong specific contacts with the uptake machinery. Other core bases had weaker roles when considered individually, as did the T-tracts, but interaction effects between these were also determinants of uptake. The properties of genomic uptake sequences are also constrained by mutational biases and selective forces acting on USSs with coding and termination functions. Our findings define constraints on gene transfer by natural transformation and suggest how the DNA uptake machinery overcomes the physical constraints imposed by stiff highly charged DNA molecules.     

Towards a Drosophila genome map.

A physical map of the genome of Drosophila melanogaster has been created using 965 yeast artificial chromosome (YAC) clones assigned to locations in the cytogenetic map by in situ hybridization with the polytene salivary gland chromosomes. Clones with insert sizes averaging about 200 kb, totaling 1.7 genome equivalents, have been mapped. More than 80% of the euchromatic genome is included in the mapped clones, and 75% of the euchromatic genome is included in 161 cytological contigs ranging in size up to 2.5 Mb (average size 510 kb). On the other hand, YAC coverage of the one-third of the genome constituting the heterochromatin is incomplete, and clones containing long tracts of highly repetitive simple satellite DNA sequences have not been recovered.     

Analysis of the genomic organisation of a small chromosome of Leishmania braziliensis M2903 reveals two genes encoding GTP-binding proteins, one of which belongs to a new G-protein family and is an antigen

Leishmania braziliensis M2903 contains a highly amplified small chromosome. This work is aimed at resolving its structural organization and determining whether this unusual chromosome contains specific genes encoding proteins with important functions in disease pathology or drug resistance. Our results show that the M2903 250-kb small chromosome contains LD1 sequences and has an inverted repeat structure. The LD1 sequences and two cDNAs (cDNA2 and cDNA53) were mapped on a cosmid contig, and the two cDNAs and the corresponding genomic fragments from the small chromosome were sequenced. The gene encoding cDNA2 predicts a putative GTP-binding protein with homology to other GTP-binding proteins only in the G-1 domain region; however, four other conserved motifs can be recognized. Sequence similarity to cDNA53 is located in at least five chromosomes, and its small chromosome copy is a pseudogene. An open reading frame downstream of the cDNA53 pseudogene predicts another GTP-binding protein that belongs to a new G-protein family with an unusual conserved GTP-binding domain and a newly characterized conserved sequence motif. A portion of this GTP-binding protein gene was studied previously in L. aethiopica as a recombinant antigen that reacts with human antibodies.