Analysis of the structure of a natural alternating d(TA)n sequence in yeast chromatin

Regulatory role of Non-canonical DNA Polymorphisms in human genome and their relevance in Cancer

DNA has the ability to form polymorphic structures like canonical duplex DNA and non-canonical triplex DNA, Cruciform, Z-DNA, G-quadruplex (G4), i-motifs, and hairpin structures. The alteration in the form of DNA polymorphism in the response to environmental changes influences the gene expression. Non-canonical structures are engaged in various biological functions, including chromatin epigenetic and gene expression regulation via transcription and translation, as well as DNA repair and recombination. The presence of non-canonical structures in the regulatory region of the gene alters the gene expression and affects the cellular machinery. Formation of non-canonical structure in the regulatory site of cancer-related genes either inhibits or dysregulate the gene function and promote tumour formation. In the current article, we review the influence of non-canonical structure on the regulatory mechanisms in human genome. Moreover, we have also discussed the relevance of non-canonical structures in cancer and provided information on the drugs used for their treatment by targeting these structures.

STRs: Ancient Architectures of the Genome beyond the Sequence

Short tandem repeats (STRs) are commonly defined as short runs of repetitive nucleotides, consisting of tandemly repeating 2-6- bp motif units, which are ubiquitously distributed throughout genomes. Functional STRs are polymorphic in the population, and their variations influence gene expression, which subsequently may result in pathogenic phenotypes. To understand STR phenotypic effects and their functional roles, we describe four different mutational mechanisms including the unequal crossing-over model, gene conversion, retrotransposition mechanism and replication slippage. Due to the multi-allelic nature, small length, abundance, high variability, codominant inheritance, nearly neutral evolution, extensive genome coverage and simple assaying of STRs, these markers are widely used in various types of biological research, including population genetics studies, genome mapping, molecular epidemiology, paternity analysis and gene flow studies. In this review, we focus on the current knowledge regarding STR genomic distribution, function, mutation and applications.

Non-B DNA: a major contributor to small- and large-scale variation in nucleotide substitution frequencies across the genome

Approximately 13% of the human genome can fold into non-canonical (non-B) DNA structures (e.g. G-quadruplexes, Z-DNA, etc.), which have been implicated in vital cellular processes. Non-B DNA also hinders replication, increasing errors and facilitating mutagenesis, yet its contribution to genome-wide variation in mutation rates remains unexplored. Here, we conducted a comprehensive analysis of nucleotide substitution frequencies at non-B DNA loci within noncoding, non-repetitive genome regions, their ±2 kb flanking regions, and 1-Megabase windows, using human-orangutan divergence and human single-nucleotide polymorphisms. Functional data analysis at single-base resolution demonstrated that substitution frequencies are usually elevated at non-B DNA, with patterns specific to each non-B DNA type. Mirror, direct and inverted repeats have higher substitution frequencies in spacers than in repeat arms, whereas G-quadruplexes, particularly stable ones, have higher substitution frequencies in loops than in stems. Several non-B DNA types also affect substitution frequencies in their flanking regions. Finally, non-B DNA explains more variation than any other predictor in multiple regression models for diversity or divergence at 1-Megabase scale. Thus, non-B DNA substantially contributes to variation in substitution frequencies at small and large scales. Our results highlight the role of non-B DNA in germline mutagenesis with implications to evolution and genetic diseases.

The effects of transcription and recombination on mutational dynamics of short tandem repeats

Short tandem repeats (STR) are ubiquitous components of the genomic architecture of most living organisms. Recent work has highlighted the widespread functional significance of such repeats, particularly around gene regulation, but the mutational processes underlying the evolution of these highly abundant and highly variable sequences are not fully understood. Traditional models assume that strand misalignment during replication is the predominant mechanism, but empirical data suggest the involvement of other processes including recombination and transcription. Despite this evidence, the relative influences of these processes have not previously been tested experimentally on a genome-wide scale. Using deep sequencing, we identify mutations at >200 microsatellites, across 700 generations in replicated populations of two otherwise identical sexual and asexual Saccharomyces cerevisiae strains. Using generalized linear models, we investigate correlates of STR mutability including the nature of the mutation, STR composition and contextual factors including recombination, transcription and replication origins. Sexual capability was not a significant predictor of microsatellite mutability, but, intriguingly, we identify transcription as a significant positive predictor. We also find that STR density is substantially increased in regions neighboring, but not within, recombination hotspots.

Functional Mechanisms of Microsatellite DNA in Eukaryotic Genomes

Microsatellite repeat DNA is best known for its length mutability, which is implicated in several neurological diseases and cancers, and often exploited as a genetic marker. Less well-known is the body of work exploring the widespread and surprisingly diverse functional roles of microsatellites. Recently, emerging evidence includes the finding that normal microsatellite polymorphism contributes substantially to the heritability of human gene expression on a genome-wide scale, calling attention to the task of elucidating the mechanisms involved. At present, these are underexplored, but several themes have emerged. I review evidence demonstrating roles for microsatellites in modulation of transcription factor binding, spacing between promoter elements, enhancers, cytosine methylation, alternative splicing, mRNA stability, selection of transcription start and termination sites, unusual structural conformations, nucleosome positioning and modification, higher order chromatin structure, noncoding RNA, and meiotic recombination hot spots.

Recognition of Local DNA Structures by p53 Protein

p53 plays critical roles in regulating cell cycle, apoptosis, senescence and metabolism and is commonly mutated in human cancer. These roles are achieved by interaction with other proteins, but particularly by interaction with DNA. As a transcription factor, p53 is well known to bind consensus target sequences in linear B-DNA. Recent findings indicate that p53 binds with higher affinity to target sequences that form cruciform DNA structure. Moreover, p53 binds very tightly to non-B DNA structures and local DNA structures are increasingly recognized to influence the activity of wild-type and mutant p53. Apart from cruciform structures, p53 binds to quadruplex DNA, triplex DNA, DNA loops, bulged DNA and hemicatenane DNA. In this review, we describe local DNA structures and summarize information about interactions of p53 with these structural DNA motifs. These recent data provide important insights into the complexity of the p53 pathway and the functional consequences of wild-type and mutant p53 activation in normal and tumor cells.

Cruciform structures are a common DNA feature important for regulating biological processes

DNA cruciforms play an important role in the regulation of natural processes involving DNA. These structures are formed by inverted repeats, and their stability is enhanced by DNA supercoiling. Cruciform structures are fundamentally important for a wide range of biological processes, including replication, regulation of gene expression, nucleosome structure and recombination. They also have been implicated in the evolution and development of diseases including cancer, Werner's syndrome and others.

Cruciform structures are targets for many architectural and regulatory proteins, such as histones H1 and H5, topoisomerase IIβ, HMG proteins, HU, p53, the proto-oncogene protein DEK and others. A number of DNA-binding proteins, such as the HMGB-box family members, Rad54, BRCA1 protein, as well as PARP-1 polymerase, possess weak sequence specific DNA binding yet bind preferentially to cruciform structures. Some of these proteins are, in fact, capable of inducing the formation of cruciform structures upon DNA binding. In this article, we review the protein families that are involved in interacting with and regulating cruciform structures, including (a) the junction-resolving enzymes, (b) DNA repair proteins and transcription factors, (c) proteins involved in replication and (d) chromatin-associated proteins. The prevalence of cruciform structures and their roles in protein interactions, epigenetic regulation and the maintenance of cell homeostasis are also discussed.

Superhelical destabilization in regulatory regions of stress response genes

Stress-induced DNA duplex destabilization (SIDD) analysis exploits the known structural and energetic properties of DNA to predict sites that are susceptible to strand separation under negative superhelical stress. When this approach was used to calculate the SIDD profile of the entire Escherichia coli K12 genome, it was found that strongly destabilized sites occur preferentially in intergenic regions that are either known or inferred to contain promoters, but rarely occur in coding regions. Here, we investigate whether the genes grouped in different functional categories have characteristic SIDD properties in their upstream flanks. We report that strong SIDD sites in the E. coli K12 genome are statistically significantly overrepresented in the upstream regions of genes encoding transcriptional regulators. In particular, the upstream regions of genes that directly respond to physiological and environmental stimuli are more destabilized than are those regions of genes that are not involved in these responses. Moreover, if a pathway is controlled by a transcriptional regulator whose gene has a destabilized 5′ flank, then the genes (operons) in that pathway also usually contain strongly destabilized SIDD sites in their 5′ flanks. We observe this statistically significant association of SIDD sites with upstream regions of genes functioning in transcription in 38 of 43 genomes of free-living bacteria, but in only four of 18 genomes of endosymbionts or obligate parasitic bacteria. These results suggest that strong SIDD sites 5′ to participating genes may be involved in transcriptional responses to environmental changes, which are known to transiently alter superhelicity. We propose that these SIDD sites are active and necessary participants in superhelically mediated regulatory mechanisms governing changes in the global pattern of gene expression in prokaryotes in response to physiological or environmental changes.

DNA in vivo experiences regulated amounts of untwisting stress. If sufficiently large, these stresses can destabilize the double helix at specific locations. These sites then become favored locations for strand separations. Gene expression and DNA replication, the two major jobs of DNA, both require the strands of the duplex to be separated. Thus, events that affect the ease of strand separation can regulate the initiation of these processes. Stress-induced DNA duplex destabilization (SIDD) has been implicated in mechanisms regulating several biological processes, including the initiation of gene expression and replication. We have developed computational methods that accurately predict the locations and extents of destabilization within genomic DNA sequences that occur in response to specified stress levels. Here, we report that the easily destabilized sites we find in the Escherichia coli K12 genome are statistically significantly overrepresented in the upstream regions of genes encoding proteins that regulate transcription. In particular, the regions upstream of genes that directly respond to physiological and environmental stimuli are more destabilized than are those regions of genes that are not involved in these responses. These results suggest that strong SIDD sites upstream of participating genes may be involved in transcriptional responses to environmental changes.

Characterization of Z-DNA as a nucleosome-boundary element in yeast Saccharomyces cerevisiae

In this article, the effect of a d(CG) DNA dinucleotide repeat sequence on RNA polymerase II transcription is examined in yeast Saccharomyces cerevisiae. Our previous report shows that a d(CG)n dinucleotide repeat sequence located proximally upstream of the TATA box enhances transcription from a minimal CYC1 promoter in a manner that depends on its surrounding negative supercoiling. Here, we demonstrate that the d(CG)9 repeat sequence stimulates gene activity by forming a Z-DNA secondary structure. Furthermore, the extent of transcriptional enhancement by Z-DNA is promoter-specific and determined by its separation distance relative to the TATA box. The stimulatory effect exerted by promoter proximal Z-DNA is not affected by helical phasing relative to the TATA box, suggesting that Z-DNA effects transcription without interacting with the general transcription machinery by looping-out the intervening DNA. A nucleosome-scanning assay reveals that the d(CG)9 repeat sequence in the Z conformation blocks nucleosome formation, and it is found in the linker DNA with two flanking nucleosomes. This result suggests that Z-DNA formation proximally upstream of a promoter is sufficient to demarcate the boundaries of its neighboring nucleosomes, which produces transcriptionally favorable locations for the TATA box near the nucleosomal DNA-entry site and at dyad positions on the nucleosome. These findings suggest that Z-DNA formation in chromatin is a part of the "genomic code" for nucleosome positioning in vivo.

The analysis of stress-induced duplex destabilization in long genomic DNA sequences

We present a method for calculating predicted locations and extents of stress-induced DNA duplex destabilization (SIDD) as functions of base sequence and stress level in long DNA molecules. The base pair denaturation energies are assigned individually, so the influences of near neighbors, methylated bases, adducts, or lesions can be included. Sample calculations indicate that copolymeric energetics give results that are close to those derived when full near-neighbor energetics are used; small but potentially informative differences occur only in the calculated SIDD properties of moderately destabilized regions. The method presented here for analyzing long sequences calculates the destabilization properties within windows of fixed length N, with successive windows displaced by an offset distance d(o). The final values of the relevant destabilization parameters for each base pair are calculated as weighted averages of the values computed for each window in which that base pair appears. This approach implicitly assumes that the strength of the direct coupling between remote base pairs that is induced by the imposed stress attenuates with their separation distance. This strategy enables calculations of the destabilization properties of DNA sequences of any length, up to and including complete chromosomes. We illustrate its utility by calculating the destabilization properties of the entire E. coli genomic DNA sequence. A preliminary analysis of the results shows that promoters are associated with SIDD regions in a highly statistically significant manner, suggesting that SIDD attributes may prove useful in the computational prediction of promoter locations in prokaryotes.

Stress-induced DNA duplex destabilization (SIDD) in the E. coli genome: SIDD sites are closely associated with promoters

We present the first analysis of stress-induced DNA duplex destabilization (SIDD) in a complete chromosome, the Escherichia coli K12 genome. We used a newly developed method to calculate the locations and extents of stress-induced destabilization to single-base resolution at superhelix density sigma = -0.06. We find that SIDD sites in this genome show a statistically highly significant tendency to avoid coding regions. And among intergenic regions, those that either contain documented promoters or occur between divergently transcribing coding regions, and hence may be inferred to contain promoters, are associated with strong SIDD sites in a statistically highly significant manner. Intergenic regions located between convergently transcribing genes, which are inferred not to contain promoters, are not significantly enriched for destabilized sites. Statistical analysis shows that a strongly destabilized intergenic region has an 80% chance of containing a promoter, whereas an intergenic region that does not contain a strong SIDD site has only a 24% chance. We describe how these observations may illuminate specific mechanisms of regulation, and assist in the computational identification of promoter locations in prokaryotes.

Cytosine deamination plays a primary role in the evolution of mammalian isochores

DNA melting is rate-limiting for cytosine deamination, from which we infer that the rate of cytosine deamination should decline twofold for each 10% increase in GC content. Analysis of human DNA sequence data confirms that this is the case for 5-methylcytosine. Several lines of evidence further confirm that it is also the case for unmethylated cytosine and that cytosine deamination causes the majority of all C-->T and G-->A transitions in mammals. Thus, cytosine deamination and DNA base composition each affect the other, forming a positive feedback loop that facilitates divergent genetic drift to high or low GC content. Because a 10 degrees C increase in temperature in vitro increases the rate of cytosine deamination 5. 7-fold, cytosine deamination must be highly dependent on body temperature, which is consistent with the dramatic differences between the isochores of warm-blooded versus cold-blooded vertebrates. Because this process involves both DNA melting and positive feedback, it would be expected to spread progressively (in evolutionary time) down the length of the chromosome, which is consistent with the large size of isochores in modern mammals.

Statistical analysis of yeast genomic downstream sequences reveals putative polyadenylation signals

The study of a few genes has permitted the identification of three elements that constitute a yeast polyadenyl-ation signal: the efficiency element (EE), the positioning element and the actual site for cleavage and poly-adenyl-ation. In this paper we perform an analysis of oligonucleotide composition on the sequences located downstream of the stop codon of all yeast genes. Several oligonucleotide families appear over-represented with a high significance (referred to herein as 'words'). The family with the highest over-representation includes the oligonucleotides shown experimentally to play a role as EEs. The word with the highest score is TATATA, followed, among others, by a series of single-nucleotide variants (TATGTA, TACATA, TAAATA.) and one-letter shifts (ATATAT). A position analysis reveals that those words have a high preference to be in 3' flanks of yeast genes and there they have a very uneven distribution, with a marked peak around 35 bp after the stop codon. Of the predicted ORFs, 85% show one or more of those sequences. Similar results were obtained using a data set of EST sequences. Other clusters of over-represented words are also detected, namely T- and A-rich signals. Using these results and previously known data we propose a general model for the 3' trailers of yeast mRNAs.

An inverse correlation between stress resistance and stuck fermentations in wine yeasts. A molecular study

During alcoholic fermentations yeast cells are subjected to several stress conditions and, therefore, yeasts have developed molecular mechanisms in order to resist this adverse situation. The mechanisms involved in stress response have been studied in Saccharomyces cerevisiae laboratory strains. However a better understanding of these mechanisms in wine yeasts could open the possibility to improve the fermentation process. In this work an analysis of the stress response in three wine yeasts has been carried out by studying the expression of several representative genes under several stress conditions which occur during fermentation. We propose a simplified method to study how these stress conditions affect the viability of yeast cells. Using this approach an inverse correlation between stress-resistance and stuck fermentations has been found. We also have preliminary data about the use of the HSP12 gene as a molecular marker for stress-resistance in wine yeasts.

The yeast FBP1 poly(A) signal functions in both orientations and overlaps with a gene promoter

This report provides an analysis of a region of chromosome XII in which the FBP1 and YLR376c genes transcribe in the same direction. Our investigation indicates that the Saccharomyces cerevisiae FBP1 gene contains strong signals for polyadenylation and transcription termination in both orientations in vivo . A (TA)14 element plays a major role in directing polyadenylation in both orientations. While this region has four nonoverlapping copies of a TATATA hexanucleotide, which is a very potent polyadenylation efficiency element in yeast, it alone is not sufficient for full activation in the reverse orientation of a cluster of downstream poly(A) sites, and an additional upstream sequence is required. The putative RNA hairpin formed from the (TA)14 element is not involved in 3'-end formation. Surprisingly, deletion of the entire (TA)14 stretch affects transcription termination in the reverse orientation, in contrast to our previous results with the forward orientation, indicating that the transcription termination element operating in the reverse orientation has very different sequence requirements. Promoter elements for the YLR376c gene overlap with the signal for FBP1 3'-end formation. To our knowledge, this is the first time that overlapping of both types of regulatory signals has been found in two adjacent yeast genes.

Induction of microsatellite instability by oxidative DNA damage

Instability of repetitive sequences, both in intronic sequences and within coding regions, has been demonstrated to be a hallmark of genomic instability in human cancer. Understanding how these mutational events arise may provide an opportunity for prevention or early intervention in cancer development. To study the source of this instability, we have identified a region of the beta-lactamase gene that is tolerant to the insertion of fragments of exogenous DNA as large as 1,614 bp with minimal loss of enzyme activity, as determined by antibiotic resistance. Fragments inserted out-of-frame render Escherichia coli sensitive to antibiotic, and compensatory frameshift mutations that restore the reading frame of beta-lactamase can be selected on the basis of antibiotic resistance. We have utilized this site to insert a synthetic microsatellite sequence within the beta-lactamase gene and selected for mutations yielding frameshifts. This assay provides for detection of one frameshift mutation in a background of 10(6) wild-type sequences. Mismatch repair deficiency increased the observed frameshift frequency approximately 300-fold. Exposure of plasmid containing microsatellite sequences to hydrogen peroxide resulted in frameshift mutations that were localized exclusively to the microsatellite sequences, whereas DNA damage by UV or N-methyl-N'-nitro-N-nitrosoguanidine did not result in enhanced mutagenesis. We postulate that in tumor cells, endogenous production of oxygen free radicals may be a major factor in promoting instability of microsatellite sequences. This beta-lactamase assay may provide a sensitive methodology for the detection and quantitation of mutations associated with the development of cancer.

Stress-induced duplex DNA destabilization in scaffold/matrix attachment regions

S/MARs are DNA elements 300 to several thousand base-pairs long, which are operationally defined by their affinity for the nuclear scaffold or matrix. S/MARs occur exclusively in eukaryotic genomes, where they mediate several functions. Because S/MARs do not have a clearcut consensus sequence, the characteristics that define their activity are thought to be structural. Ubiquitous S/MAR binding proteins have been identified, but to date no unique binding sequence or structural motif has been found. Here we show by computational analysis that S/MARs conform to a specific design whose essential attribute is the presence of stress-induced base-unpairing regions (BURs). Stress-induced destabilization (SIDD) profiles are calculated using a previously developed statistical mechanical procedure in which the superhelical deformation is partitioned between strand separation, twisting within denatured regions, and residual superhelicity. The results of these calculations show that BURs exhibit a succession of evenly spaced destabilized sites that would render part or all of the S/MAR sequence single stranded at sufficient superhelicity. These analyses are performed for a range of sequenced S/MAR elements from the borders of eukaryotic gene domains, from centromeres, and from positions where S/MARs are known to support the action of an enhancer. The results reported here are in excellent agreement with earlier in vitro chemical reactivity studies. This approach demonstrates the potential for computational analysis to predict the points of division of the eukaryotic genome into functional units (domains), and also to locate certain cis-regulatory sequences.