Stability, structure and complexity of yeast chromosome III

The Genomic HyperBrowser: inferential genomics at the sequence level

The immense increase in the generation of genomic scale data poses an unmet analytical challenge, due to a lack of established methodology with the required flexibility and power. We propose a first principled approach to statistical analysis of sequence-level genomic information. We provide a growing collection of generic biological investigations that query pairwise relations between tracks, represented as mathematical objects, along the genome. The Genomic HyperBrowser implements the approach and is available at http://hyperbrowser.uio.no.

A stitch in time: efficient computation of genomic DNA melting bubbles

Background

It is of biological interest to make genome-wide predictions of the locations of DNA melting bubbles using statistical mechanics models. Computationally, this poses the challenge that a generic search through all combinations of bubble starts and ends is quadratic.

Results

An efficient algorithm is described, which shows that the time complexity of the task is O(NlogN) rather than quadratic. The algorithm exploits that bubble lengths may be limited, but without a prior assumption of a maximal bubble length. No approximations, such as windowing, have been introduced to reduce the time complexity. More than just finding the bubbles, the algorithm produces a stitch profile, which is a probabilistic graphical model of bubbles and helical regions. The algorithm applies a probability peak finding method based on a hierarchical analysis of the energy barriers in the Poland-Scheraga model.

Conclusion

Exact and fast computation of genomic stitch profiles is thus feasible. Sequences of several megabases have been computed, only limited by computer memory. Possible applications are the genome-wide comparisons of bubbles with promotors, TSS, viral integration sites, and other melting-related regions.

The human genomic melting map

In a living cell, the antiparallel double-stranded helix of DNA is a dynamically changing structure. The structure relates to interactions between and within the DNA strands, and the array of other macromolecules that constitutes functional chromatin. It is only through its changing conformations that DNA can organize and structure a large number of cellular functions. In particular, DNA must locally uncoil, or melt, and become single-stranded for DNA replication, repair, recombination, and transcription to occur. It has previously been shown that this melting occurs cooperatively, whereby several base pairs act in concert to generate melting bubbles, and in this way constitute a domain that behaves as a unit with respect to local DNA single-strandedness. We have applied a melting map calculation to the complete human genome, which provides information about the propensities of forming local bubbles determined from the whole sequence, and present a first report on its basic features, the extent of cooperativity, and correlations to various physical and biological features of the human genome. Globally, the melting map covaries very strongly with GC content. Most importantly, however, cooperativity of DNA denaturation causes this correlation to be weaker at resolutions fewer than 500 bps. This is also the resolution level at which most structural and biological processes occur, signifying the importance of the informational content inherent in the genomic melting map. The human DNA melting map may be further explored at http://meltmap.uio.no.

Through a genome, darkly: comparative analysis of plant chromosomal DNA

Plant nuclear genomes encompass a wide range of variation in size and nucleotide composition with diverse arrangements of chromosomal segments, repetitive sequences and distribution of genes. Comparative genomic analysis may be undertaken at different levels of organisation, which are reflected in this review, together with a focus on the genetic and functional significance of the observed variation. Patterns of genome organisation have been revealed which reflect the different underlying mechanisms and constraints driving change. Thus comparative issues of genome size, nucleotide sequence composition and genome heterogeneity are provided as a background to understanding the different levels of segmental and repetitive sequence duplication and distribution of genes. The extent of synteny and collinearity revealed by recent genetic and sequence comparisons is discussed, together with a consideration of problems associated with such analyses. The possible origins and mechanisms of variation in genome size and organisation are covered, including the prevalence of duplication at different levels of organisation. The likely genetic, functional and adaptive consequences of replicated loci are discussed with evidence from comparative studies. The scope for comparative analysis of epigenetic plant genome variation is considered. Finally, opportunities for applying comparative genomics to isolating genes and understanding complex crop genomes are addressed.

Genes and the physics of the DNA double-helix

The processing of the genetic information stored in the double-helical DNA implies the separation of the two strands, the physics of which is described by the helix-coil transition model. Is there a relationship between genetic maps and DNA physical stability maps that plot the sequence-specific propensity for the thermal disruption of the double-helix? Here, with appropriate methodological formulations, such maps are derived for a large set of sequences, including complete genomes. The superposition of the two maps leads to a contrasted picture with correlations ranging between two extremes: from almost perfect (with the genes precisely delineated as stable regions) to more or less complete unrelatedness. The simplest explanation for the results is that the observed striking correlations correspond to the relics of a primeval organisation of the genetic message, with the physics of DNA playing a role in the delimitation of coding regions. In order to trace the evolutionary fate of this signal further, a detailed study of the yeast complete genome is performed. In this study, the superposition of the genetic and physical stability maps is examined in the light of information concerning gene duplication. On the basis of this analysis it is concluded that the 'signature' associated with the supposed archaic signal is in the process of being erased, most probably because the underlying feature is no longer under selective pressure. There are many evolutionary implications for the results presented and for their proposed interpretations, notably concerning models of mutational dynamics in relation to erasure processes.

Sequence complexity and DNA curvature

A linguistic complexity measure was applied to the complete genomes of HIV-1, Escherichia coli, Bacillus subtilis, Haemophilus influenzae, Mycoplasma genitalium, and to long human and yeast genomic fragments. Complexity values averaged over entire genomic sequences were compared, as were predicted average values of intrinsic DNA curvature. We found that both the most curved and the least complex fragments are located preferentially in non-coding parts of the genome. Analysis of location of the most curved and the simplest regions in bacteria showed that the low-complexity segments are preferentially located in close proximity to the highly curved sequences, which are, in turn, placed from 100 to 200 bases upstream to the start of the nearest coding sequence. We conclude that the parallel analysis of sequence complexity and DNA curvature might provide important information about sequence-structure-function relationship in genomes.

Preferential incision of interstrand crosslinks induced by 8-methoxypsoralen plus UVA in yeast during the cell cycle

Interstrand crosslink (ICL) induction by 8-methoxypsoralen plus UVA and the incision step of the repair have been investigated during the mitotic cell cycle of haploid Saccharomyces cerevisiae. Cells were synchronised by elutriation and events were examined at the level of the MAT alpha and the HML alpha loci in a SIR strain. The DNA sequence of these two loci is identical, but the MAT alpha locus may be replicated earlier in S phase and is transcriptionally active while the HML alpha locus may be replicated later in S phase and is transcriptionally inactive because of Sir repression that creates a heterochromatin-like structure at this locus. ICL were induced to similar extents in both loci during the stages of the cell cycle examined, and these levels were identical to those reported for asynchronous cultures. Preferential incisions occurred for ICL in the MAT alpha locus compared to those in the HML alpha locus, independently of the cell cycle phase studied. The levels of incision were comparable for events in the early G1 phase (eG1), late G1 phase (lG1), early S phase (eS), middle S phase (mS), late S phase (lS) or G2 phase (G2). Thus the preferential incision of ICL observed previously in asynchronous cell culture is maintained throughout the cell cycle and, surprisingly, occurs equally well in G1. Here the opportunities for recombination to further process the incised damaged are substantially limited compared to those in the S and G2 phases.

Distribution of topoisomerase II-mediated cleavage sites and relation to structural and functional landmarks in 830 kb of Drosophila DNA

The pattern of sites for cleavage mediated by topoisomerase II was determined in 830 kb of cloned DNA from the Drosophila X chromosome, with the objectives of comparing it with mapped structural and functional landmarks and examining if the correlations with such landmarks reported in individual loci can be generalized to a region approximately 100 times longer. The relative frequencies of topoisomerase II cleavage sites in 247 restriction fragments from 67 clones were quantified by hybridization with probes prepared from DNA fragments which abutted all cleavage sites in each clone, selected through the covalently bound topoisomerase II subunit; the specificity and quantitative nature of this method were demonstrated using a plasmid DNA model. The 12 restriction fragments with strong nuclear scaffold attachment (SAR) activity, of which seven possess autonomous replication (ARS) activity, show statistically strong coincidence or contiguity ( P </=0.11) with regions of high topoisomerase II cleavage site frequency. These regions show no correlation with repetitive sequence or A/T or C/G content and some extend over >10 kb; their sensitivity is therefore unlikely to be due to alternating purine-pyrimidine repeats or regions of Z conformation, which are preferred motifs. The hypothesis that they possess intrinsic curvature is consistent with the similarity of their length and spacing to regions of predicted curvature in the 315 kb DNA of Saccharomyces cerevisiae chromosome III and with the reported strong binding preference of topoisomerase II for curved DNA. The topoisomerase II cleavage pattern in this DNA further shows that its relationships to functional properties seen in individual loci, especially to MAR/SAR and ARS activity and to the restricted accessibility of DNA to topoisomerase II in vivo, can be generalized to much longer regions of the genome.

Constraints on mutability in a multiallelic gene family

A highly variable family of related DNA sequences was examined in order to determine the effect of local sequence environment on substitution mutation; 29 sequences from the Brassica self-incompatibility gene family, which possess a high level of nonsynonymous mutations, were aligned and grouped according to their similarity and function. The level and distribution of substitution mutations were calculated. A nonrandom distribution of sequence variation was observed along the sequences. The effect of neighbor biases and structural and thermodynamic measures were then compared in the absence of strong codon conservation. Biases were observed in the rates of substitution of the same base pair in different local sequence environments. The effect of the 5' neighbor was such that nucleotide A or C was associated with more mutations than G or T. There were significant interactions of certain dinucleotides with the frequency of mutation. Sequence-dependent measures of helical stability, intrinsic curvature, components of curvature, and stacking interactions were calculated for each sequence. Decreased helical stability was found to be associated with increased mutation. The compound measure of curvature, calculated according to the "wedge" model, showed little association with mutation. However, the components of increased wedge angle and decreased twist both showed an association with increased mutation. A small effect of A-type DNA stacking was found to be associated with mutated bases.