The sequence (dC-dA)n X (dG-dT)n forms left-handed Z-DNA in negatively supercoiled plasmids

ZSCAN4-binding motif-TGCACAC is conserved and enriched in CA/TG microsatellites in both mouse and human genomes

The Zinc finger and SCAN domain containing 4 (ZSCAN4) protein, expressed transiently in pluripotent stem cells, gametes, and early embryos, extends telomeres, enhances genome stability, and improves karyotypes in mouse embryonic stem (mES) cells. To gain insights into the mechanism of ZSCAN4 function, we identified genome-wide binding sites of endogenous ZSCAN4 protein using ChIP-seq technology in mouse and human ES cells, where the expression of endogenous ZSCAN4 was induced by treating cells with retinoic acids or by overexpressing DUX4. We revealed that both mouse and human ZSCAN4 bind to the TGCACAC motif located in CA/TG microsatellite repeats, which are known to form unstable left-handed duplexes called Z-DNA that can induce double-strand DNA breaks and mutations. These ZSCAN4 binding sites are mostly located in intergenic and intronic regions of the genomes. By generating ZSCAN4 knockout in human ES cells, we showed that ZSCAN4 does not seem to be involved in transcriptional regulation. We also found that ectopic expression of mouse ZSCAN4 enhances the suppression of chromatin at ZSCAN4-binding sites. These results together suggest that some of the ZSCAN4 functions are mediated by binding to the error-prone regions in mouse and human genomes.

Z-form extracellular DNA is a structural component of the bacterial biofilm matrix

Biofilms are community architectures adopted by bacteria inclusive of a self-formed extracellular matrix that protects resident bacteria from diverse environmental stresses and, in many species, incorporates extracellular DNA (eDNA) and DNABII proteins for structural integrity throughout biofilm development. Here, we present evidence that this eDNA-based architecture relies on the rare Z-form. Z-form DNA accumulates as biofilms mature and, through stabilization by the DNABII proteins, confers structural integrity to the biofilm matrix. Indeed, substances known to drive B-DNA into Z-DNA promoted biofilm formation whereas those that drive Z-DNA into B-DNA disrupted extant biofilms. Importantly, we demonstrated that the universal bacterial DNABII family of proteins stabilizes both bacterial- and host-eDNA in the Z-form in situ. A model is proposed that incorporates the role of Z-DNA in biofilm pathogenesis, innate immune response, and immune evasion.

Sequence-dependent cost for Z-form shapes the torsion-driven B-Z transition via close interplay of Z-DNA and DNA bubble

Despite recent genome-wide investigations of functional DNA elements, the mechanistic details about their actions remain elusive. One intriguing possibility is that DNA sequences with special patterns play biological roles, adopting non-B-DNA conformations. Here we investigated dynamics of thymine-guanine (TG) repeats, microsatellite sequences and recurrently found in promoters, as well as cytosine–guanine (CG) repeats, best-known Z-DNA forming sequence, in the aspect of Z-DNA formation. We measured the energy barriers of the B–Z transition with those repeats and discovered the sequence-dependent penalty for Z-DNA generates distinctive thermodynamic and kinetic features in the torque-induced transition. Due to the higher torsional stress required for Z-form in TG repeats, a bubble could be induced more easily, suppressing Z-DNA induction, but facilitate the B–Z interconversion kinetically at the transition midpoint. Thus, the Z-form by TG repeats has advantages as a torsion buffer and bubble selector while the Z-form by CG repeats likely behaves as torsion absorber. Our statistical physics model supports quantitatively the populations of Z-DNA and reveals the pivotal roles of bubbles in state dynamics. All taken together, a quantitative picture for the transition was deduced within the close interplay among bubbles, plectonemes and Z-DNA.

DNA fragility in the parallel evolution of pelvic reduction in stickleback fish

Evolution generates a remarkable breadth of living forms, but many traits evolve repeatedly, by mechanisms that are still poorly understood. A classic example of repeated evolution is the loss of pelvic hindfins in stickleback fish (Gasterosteus aculeatus). Repeated pelvic loss maps to recurrent deletions of a pelvic enhancer of the Pitx1 gene. Here, we identify molecular features contributing to these recurrent deletions. Pitx1 enhancer sequences form alternative DNA structures in vitro and increase double-strand breaks and deletions in vivo. Enhancer mutability depends on DNA replication direction and is caused by TG-dinucleotide repeats. Modeling shows that elevated mutation rates can influence evolution under demographic conditions relevant for sticklebacks and humans. DNA fragility may thus help explain why the same loci are often used repeatedly during parallel adaptive evolution.

Exploring biomolecular energy landscapes

The potential energy landscape perspective provides both a conceptual and a computational framework for predicting, understanding and designing molecular properties. In this Feature Article, we highlight some recent advances that greatly facilitate structure prediction and analysis of global thermodynamics and kinetics in proteins and nucleic acids. The geometry optimisation procedures, on which these calculations are based, can be accelerated significantly using local rigidification of selected degrees of freedom, and through implementations on graphics processing units. Results of progressive local rigidification are first summarised for trpzip1, including a systematic analysis of the heat capacity and rearrangement rates. Benchmarks for all the essential optimisation procedures are then provided for a variety of proteins. Applications are then illustrated from a study of how mutation affects the energy landscape for a coiled-coil protein, and for transitions in helix morphology for a DNA duplex. Both systems exhibit an intrinsically multifunnel landscape, with the potential to act as biomolecular switches.

Probing helical transitions in a DNA duplex

The complex conformational change from B-DNA to Z-DNA requires inversion of helix-handedness. Multiple degrees of freedom are intricately coupled during this transition, and formulating an appropriate reaction coordinate that captures the underlying complexity would be problematic. In this contribution, we adopt an alternative approach, based on the potential energy landscape perspective, to construct a kinetic transition network. Microscopic insight into the B → Z transition is provided in terms of geometrically defined discrete paths consisting of local minima and the transition states that connect them. We find that the inversion of handedness can occur via two competing mechanisms, either involving stretched intermediates, or a B-Z junction, in agreement with previous predictions. The organisation of the free energy landscape further suggests that this process is likely to be slow under physiological conditions. Our results represent a key step towards decoding the more intriguing features of the B → Z transition, such as the role of ionic strength and negative supercoiling in reshaping the landscape.

FACT is a sensor of DNA torsional stress in eukaryotic cells

Transitions of B-DNA to alternative DNA structures (ADS) can be triggered by negative torsional strain, which occurs during replication and transcription, and may lead to genomic instability. However, how ADS are recognized in cells is unclear. We found that the binding of candidate anticancer drug, curaxin, to cellular DNA results in uncoiling of nucleosomal DNA, accumulation of negative supercoiling and conversion of multiple regions of genomic DNA into left-handed Z-form. Histone chaperone FACT binds rapidly to the same regions via the SSRP1 subunit in curaxin-treated cells. In vitro binding of purified SSRP1 or its isolated CID domain to a methylated DNA fragment containing alternating purine/pyrimidines, which is prone to Z-DNA transition, is much stronger than to other types of DNA. We propose that FACT can recognize and bind Z-DNA or DNA in transition from a B to Z form. Binding of FACT to these genomic regions triggers a p53 response. Furthermore, FACT has been shown to bind to other types of ADS through a different structural domain, which also leads to p53 activation. Thus, we propose that FACT acts as a sensor of ADS formation in cells. Recognition of ADS by FACT followed by a p53 response may explain the role of FACT in DNA damage prevention.

Nbs1 ChIP-Seq Identifies Off-Target DNA Double-Strand Breaks Induced by AID in Activated Splenic B Cells

Activation-induced cytidine deaminase (AID) is required for initiation of Ig class switch recombination (CSR) and somatic hypermutation (SHM) of antibody genes during immune responses. AID has also been shown to induce chromosomal translocations, mutations, and DNA double-strand breaks (DSBs) involving non-Ig genes in activated B cells. To determine what makes a DNA site a target for AID-induced DSBs, we identify off-target DSBs induced by AID by performing chromatin immunoprecipitation (ChIP) for Nbs1, a protein that binds DSBs, followed by deep sequencing (ChIP-Seq). We detect and characterize hundreds of off-target AID-dependent DSBs. Two types of tandem repeats are highly enriched within the Nbs1-binding sites: long CA repeats, which can form Z-DNA, and tandem pentamers containing the AID target hotspot WGCW. These tandem repeats are not nearly as enriched at AID-independent DSBs, which we also identified. Msh2, a component of the mismatch repair pathway and important for genome stability, increases off-target DSBs, similar to its effect on Ig switch region DSBs, which are required intermediates during CSR. Most of the off-target DSBs are two-ended, consistent with generation during G1 phase, similar to DSBs in Ig switch regions. However, a minority are one-ended, presumably due to conversion of single-strand breaks to DSBs during replication. One-ended DSBs are repaired by processes involving homologous recombination, including break-induced replication repair, which can lead to genome instability. Off-target DSBs, especially those present during S phase, can lead to chromosomal translocations, deletions and gene amplifications, resulting in the high frequency of B cell lymphomas derived from cells that express or have expressed AID.

Activation-induced cytidine deaminase (AID) is required for diversifying antibodies during immune responses, and it does this by introducing mutations and DNA breaks into antibody genes. How AID is targeted is not understood, and it induces chromosomal translocations, mutations, and double-strand breaks (DSBs) at sites other than antibody genes in activated B cells. To determine what makes an off-target DNA site a target for AID-induced DSBs, we identify and characterize hundreds of genome-wide DSBs induced by AID during B cell activation. Interestingly, many of the DSBs are within or adjacent to two types of tandemly repeated simple sequences, which have characteristics that might explain why they are targeted. We find that most of the DSBs are two-ended, consistent with their generation during G1 phase of the cell cycle, which is when AID induces DNA breaks in antibody genes. However, a minority is one-ended, consistent with replication encountering an AID-induced single-strand break, thereby creating a DSB. Both types of off-target DSBs, but especially those present during S phase of the cell cycle, lead to chromosomal translocations, deletions and gene amplifications that can promote B cell lymphomagenesis.

Ion distributions around left- and right-handed DNA and RNA duplexes: a comparative study

The ion atmosphere around nucleic acids is an integral part of their solvated structure. However, detailed aspects of the ionic distribution are difficult to probe experimentally, and comparative studies for different structures of the same sequence are almost non-existent. Here, we have used large-scale molecular dynamics simulations to perform a comparative study of the ion distribution around (5'-CGCGCGCGCGCG-3')2 dodecamers in solution in B-DNA, A-RNA, Z-DNA and Z-RNA forms. The CG sequence is very sensitive to ionic strength and it allows the comparison with the rare but important left-handed forms. The ions investigated include Na(+), K(+) and Mg(2 +), with various concentrations of their chloride salts. Our results quantitatively describe the characteristics of the ionic distributions for different structures at varying ionic strengths, tracing these differences to nucleic acid structure and ion type. Several binding pockets with rather long ion residence times are described, both for the monovalent ions and for the hexahydrated Mg[(H2O)6](2+) ion. The conformations of these binding pockets include direct binding through desolvated ion bridges in the GpC steps in B-DNA and A-RNA; direct binding to backbone oxygens; binding of Mg[(H2O)6](2+) to distant phosphates, resulting in acute bending of A-RNA; tight 'ion traps' in Z-RNA between C-O2 and the C-O2' atoms in GpC steps; and others.

Activation-induced cytidine deaminase induces reproducible DNA breaks at many non-Ig Loci in activated B cells

After immunization or infection, activation-induced cytidine deaminase (AID) initiates diversification of immunoglobulin (Ig) genes in B cells, introducing mutations within the antigen-binding V regions (somatic hypermutation, SHM) and double-strand DNA breaks (DSBs) into switch (S) regions, leading to antibody class switch recombination (CSR). We asked if, during B cell activation, AID also induces DNA breaks at genes other than IgH genes. Using a nonbiased genome-wide approach, we have identified hundreds of reproducible, AID-dependent DSBs in mouse splenic B cells shortly after induction of CSR in culture. Most interestingly, AID induces DSBs at sites syntenic with sites of translocations, deletions, and amplifications found in human B cell lymphomas, including within the oncogene B cell lymphoma11a (bcl11a)/evi9. Unlike AID-induced DSBs in Ig genes, genome-wide AID-dependent DSBs are not restricted to transcribed regions and frequently occur within repeated sequence elements, including CA repeats, non-CA tandem repeats, and SINEs.

Interaction of isoquinoline alkaloids with polymorphic DNA structures

The interaction of berberine, palmatine, and coralyne with the B, Z, and H(L) form of poly[d(G-C)] was studied. Berberine and palmatine showed moderate binding to the B form, while coralyne showed higher binding, as revealed from spectroscopic and thermodynamic data. Berberine and coralyne binding to the B form was exothermic and enthalpy-driven, while palmatine showed exothermic binding which was favored by both negative enthalpy and negative entropy changes. Berberine and palmatine neither bind nor converted the Z-form structure to B form. Coralyne, on the other hand, exhibited a strong binding affinity to Z DNA structure that was enthalpy-driven. Berberine binding to the H(L) form was cooperative, exothermic, and favored by both negative enthalpy and negative entropy changes with the formation of an induced CD band. Palmatine showed weak binding, while coralyne showed a strong binding with the H(L) form. The structural differences in the isoquinoline alkaloids appear to influence the affinity and mode of interactions with these polymorphic DNA structures.

The effect of a promoter polymorphism on the transcription of nitric oxide synthase 1 and its relevance to Parkinson's disease

Transcriptional changes of the enzyme nitric oxide synthase I (NOS1) are believed to play a role in the development of many diseases. The gene for NOS1 has 12 alternative first exons (1A-1L). The 1F exon is one of the most highly utilized first exons in the brain and has a polymorphism ((TG)(m)TA(TG)(n)) located in its promoter region. The polymorphism's length has been suggested to affect NOS1 transcription and play a role in Parkinson's disease (PD); however, the actual influence of the polymorphism on NOS1 transcription has not been studied. To better characterize the links of the polymorphism with PD, a genotyping study was done comparing polymorphism length among 170 PD patients and 150 age-matched controls. The pattern of changes between the two group's allele frequencies shows statistical significance (P = 0.0359). The smallest polymorphism sizes are more predominant among PD patients than controls. To study the effects of this polymorphism on NOS1 gene transcription, reporter gene constructs were made by cloning the NOS1 1F promoter with polymorphism lengths of either 42, 54, or 62 bp in front of the luciferase gene and transfecting them into HeLa or Sk-N-MC cells. NOS1-directed reporter gene constructs with the 62-bp polymorphism increased transcription of luciferase 2.2-fold in HeLa and 1.8-fold in Sk-N-MC cells compared with reporter gene constructs with the 42-bp polymorphism. These data suggest that if smaller polymorphism size contributes to the higher NOS1 levels in PD patients, an as yet unknown transcriptional mechanism is required.

Exploring the relationship between polymorphic (TG/CA)n repeats in intron 1 regions and gene expression

The putative role of (TG/CA)n repeats in the regulation of transcription has recently been reported for several cancer- and disease-related genes, including the genes encoding the epidermal growth factor receptor (EGFR), hydroxysteroid (11-beta) dehydrogenase 2 (HSD11B2) and interferon-gamma (IFNG). These studies indicated a correlation between gene expression levels and the presence or length of (TG/CA)n repeats in their intron 1 regions. A genome-wide search for genes with similar features may provide evidence of whether these dinucleotide repeats represent a class of universal regulators of gene expression, which has recently begun to be investigated as a quantitative complex phenotype. Using a public database of simple repeats, we identified 330 genes containing potentially polymorphic long (TG/CA)n repeats (n >or= 12) in their intron 1 regions. One known physiological pathway, the calcium signalling pathway, was found to be enriched among the genes containing long repeats. In addition, certain biological processes, such as cation transport, signal transduction and ion transport, were found to be enriched in these genes. Genotyping of the long repeats showed that the majority of these dinucleotide repeats were polymorphic in the HapMap CEU (Caucasians from Utah, USA) samples of northern and western European ancestry. Evidence for a significant association between these repeats and gene expression was not observed in the genes selected based on their expression profiles in the HapMap CEU samples. Our current findings, therefore, do not support a role for these repeats as a class of universal gene expression regulators. A more comprehensive evaluation of the relationship between these repeats and gene expression, potentially in other tissues, may be necessary to illustrate their roles in gene regulation in the future.

Methods to determine DNA structural alterations and genetic instability

Chromosomal DNA is a dynamic structure that can adopt a variety of non-canonical (i.e., non-B) conformations. In this regard, at least 10 different forms of non-B DNA conformations have been identified; many of them have been found to be mutagenic, and associated with human disease development. Despite the importance of non-B DNA structures in genetic instability and DNA metabolic processes, mechanisms by which instability occurs remain largely undefined. The purpose of this review is to summarize current methodologies that are used to address questions in the field of non-B DNA structure-induced genetic instability. Advantages and disadvantages of each method will be discussed. A focused effort to further elucidate the mechanisms of non-B DNA-induced genetic instability will lead to a better understanding of how these structure-forming sequences contribute to the development of human disease.

Abundance of dinucleotide repeats and gene expression are inversely correlated: a role for gene function in addition to intron length

High and broad transcription of eukaryotic genes is facilitated by cost minimization, clustered localization in the genome, elevated G+C content, and low nucleosome formation potential. In this scenario, illumination of correlation between abundance of (TG/CA)(n>or=12) repeats, which are negative cis modulators of transcription, and transcriptional levels and other commonly occurring dinucleotide repeats, is required. Three independent microarray datasets were used to examine the correlation of (TG/CA)(n>or=12) and other dinucleotide repeats with gene expression. Compared with the expected equi-distribution pattern under neutral model, highly transcribed genes were poor in repeats, and conversely, weakly transcribed genes were rich in repeats. Furthermore, the inverse correlation between repeat abundance and transcriptional levels appears to be a global phenomenon encompassing all genes regardless of their breadth of transcription. This selective pattern of exclusion of (TG/CA)(n>or=12) and (AT)(n>or=12) repeats in highly transcribed genes is an additional factor along with cost minimization and elevated GC, and therefore, multiple factors govern high transcription of genes. We observed that even after controlling for the effects of GC and average intron lengths, the effect of repeats albeit somewhat weaker was persistent and definite. In the ribosomal protein coding genes, sequence analysis of orthologs suggests that negative selection for repeats perhaps occurred early in evolution. These observations suggest that negative selection of (TG/CA)(n>or=12) microsatellites in the evolution of the highly expressed genes was also controlled by gene function in addition to intron length.

Simple sequence repeats in prokaryotic genomes

Simple sequence repeats (SSRs) in DNA sequences are composed of tandem iterations of short oligonucleotides and may have functional and/or structural properties that distinguish them from general DNA sequences. They are variable in length because of slip-strand mutations and may also affect local structure of the DNA molecule or the encoded proteins. Long SSRs (LSSRs) are common in eukaryotes but rare in most prokaryotes. In pathogens, SSRs can enhance antigenic variance of the pathogen population in a strategy that counteracts the host immune response. We analyze representations of SSRs in >300 prokaryotic genomes and report significant differences among different prokaryotes as well as among different types of SSRs. LSSRs composed of short oligonucleotides (1-4 bp length, designated LSSR(1-4)) are often found in host-adapted pathogens with reduced genomes that are not known to readily survive in a natural environment outside the host. In contrast, LSSRs composed of longer oligonucleotides (5-11 bp length, designated LSSR(5-11)) are found mostly in nonpathogens and opportunistic pathogens with large genomes. Comparisons among SSRs of different lengths suggest that LSSR(1-4) are likely maintained by selection. This is consistent with the established role of some LSSR(1-4) in enhancing antigenic variance. By contrast, abundance of LSSR(5-11) in some genomes may reflect the SSRs' general tendency to expand rather than their specific role in the organisms' physiology. Differences among genomes in terms of SSR representations and their possible interpretations are discussed.

The rare crystallographic structure of d(CGCGCG)(2): the natural spermidine molecule bound to the minor groove of left-handed Z-DNA d(CGCGCG)(2) at 10 degrees C

Several crystal structure analyses of complexes of synthetic polyamine compounds, including N(1)-(2-(2-aminoethylamino))ethyl)ethane-1,2-diamine PA(222) and N(1)-(2-(2-(2-aminoethylamino)ethylamino)ethyl)ethane-1,2-diamine PA(2222), and left-handed Z-DNA d(CGCGCG)(2) have been reported. However, until now, there have been no examples of naturally occurring polyamines bound to the minor groove of the left-handed Z-DNA of d(CGCGCG)(2) molecule. We have found that spermidine, a natural polyamine, is connected to the minor groove of left-handed Z-DNA of d(CGCGCG)(2) molecule in a crystalline complex grown at 10 degrees C. The electron density of the DNA molecule was clear enough to determine that the spermidine was connected in the minor groove of two symmetry related molecules of left-handed Z-DNA d(CGCGCG)(2). This is the first example that a spermidine molecule can form a bridge conformation between two symmetry related molecules of left-handed Z-DNA d(CGCGCG)(2) in the minor groove.

Expoldb: expression linked polymorphism database with inbuilt tools for analysis of expression and simple repeats

BACKGROUND

Quantitative variation in gene expression has been proposed to underlie phenotypic variation among human individuals. A facilitating step towards understanding the basis for gene expression variability is associating genome wide transcription patterns with potential cis modifiers of gene expression.

DESCRIPTION

EXPOLDB, a novel Database, is a new effort addressing this need by providing information on gene expression levels variability across individuals, as well as the presence and features of potentially polymorphic (TG/CA)n repeats. EXPOLDB thus enables associating transcription levels with the presence and length of (TG/CA)n repeats. One of the unique features of this database is the display of expression data for 5 pairs of monozygotic twins, which allows identification of genes whose variability in expression, are influenced by non-genetic factors including environment. In addition to queries by gene name, EXPOLDB allows for queries by a pathway name. Users can also upload their list of HGNC (HUGO (The Human Genome Organisation) Gene Nomenclature Committee) symbols for interrogating expression patterns. The online application 'SimRep' can be used to find simple repeats in a given nucleotide sequence. To help illustrate primary applications, case examples of Housekeeping genes and the RUNX gene family, as well as one example of glycolytic pathway genes are provided.

CONCLUSION

The uniqueness of EXPOLDB is in facilitating the association of genome wide transcription variations with the presence and type of polymorphic repeats while offering the feature for identifying genes whose expression variability are influenced by non genetic factors including environment. In addition, the database allows comprehensive querying including functional information on biochemical pathways of the human genes. EXPOLDB can be accessed at http://expoldb.igib.res.in/expol.

Specific Interactions of Divalent Metal Ions with a DNA Duplex Containing the d(CA)n/(GT)nTandem Repeat

Divalent metal ions are essential for maintaining functional states of the DNA molecule. Their participation in DNA structure is modulated by the base sequence and varies depending on the nature of the ion. The present investigation addresses the interaction of Ca2+ ions with a tandem repeat of two CA dinucleotides, (CA)2/(TG)2. The binding of Ca2+ to the repeat is monitored by nuclear magnetic resonance (NMR) spectroscopy using chemical shift mapping. Parallel experiments monitor binding of Mg2+ ions to the repeat as well as binding of each ion to a DNA duplex in which the (CA)2/(TG)2 repeat is eliminated. The results reveal that the direction and the magnitude of chemical shift changes induced by Ca2+ ions in the NMR spectra of the repeat are different from those induced by Mg2+ ions. The differences between the two cations are significantly diminished by the elimination of the (CA)2/(TG)2 repeat. These findings suggest a specific interaction of Ca2+ ions with the (CA)2/(TG)2 motif. The specificity of the interaction resides in the two A-T base pairs of the repeat, and it involves the major groove of the first A-T base pair and both grooves of the second A-T base pair.

Analysis of distribution indicates diverse functions of simple sequence repeats in Mycoplasma genomes

Simple sequence repeats (SSRs) composed of extensive tandem iterations of a single nucleotide or a short oligonucleotide are rare in most bacterial genomes, but they are common among Mycoplasma. Some of these repeats act as contingency loci in association with families of surface antigens. By contraction or expansion during replication, these SSRs increase genetic variance of the population and facilitate avoidance of the immune response of the host. Occurrence and distribution of SSRs are analyzed in complete genomes of 11 Mycoplasma and 3 related Mollicutes in order to gain insights into functional and evolutionary diversity of the SSRs in Mycoplasma. The results revealed an unexpected variety of SSRs with respect to their distribution and composition and suggest that it is unlikely that all SSRs function as contingency loci or recombination hot spots. Various types of SSRs are most abundant in Mycoplasma hyopneumoniae, whereas Mycoplasma penetrans, Mycoplasma mobile, and Mycoplasma synoviae do not contain unusually long SSRs. Mycoplasma hyopneumoniae and Mycoplasma pulmonis feature abundant short adenine and thymine runs periodically spaced at 11 and 12 bp, respectively, which likely affect the supercoiling propensities of the DNA molecule. Physiological roles of long adenine and thymine runs in M. hyopneumoniae appear independent of location upstream or downstream of genes, unlike contingency loci that are typically located in protein-coding regions or upstream regulatory regions. Comparisons among 3 M. hyopneumoniae strains suggest that the adenine and thymine runs are rarely involved in genome rearrangements. The results indicate that the SSRs in the Mycoplasma genomes play diverse roles, including modulating gene expression as contingency loci, facilitating genome rearrangements via recombination, affecting protein structure and possibly protein-protein interactions, and contributing to the organization of the DNA molecule in the cell.

Determinants of specific binding of HMGB1 protein to hemicatenated DNA loops

Protein HMGB1 has long been known as one of the most abundant non-histone proteins in the nucleus of mammalian cells, and has regained interest recently for its function as an extracellular cytokine. As a DNA-binding protein, HMGB1 facilitates DNA-protein interactions by increasing the flexibility of the double helix, and binds specifically to distorted DNA structures. We have previously observed that HMGB1 binds with extremely high affinity to a novel DNA structure, hemicatenated DNA loops (hcDNA), in which double-stranded DNA fragments containing a tract of poly(CA).poly(TG) form a loop maintained at its base by a hemicatenane. Here, we show that the single HMGB1 domains A and B, the HMG-box domain of sex determination factor SRY, as well as the prokaryotic HMGB1-like protein HU, specifically interact with hcDNA (Kd approximately 0.5 nM). However, the affinity of full-length HMGB1 for hcDNA is three orders of magnitude higher (Kd<0.5 pM) and requires the simultaneous presence of both HMG-box domains A and B plus the acidic C-terminal tail on the molecule. Interestingly, the high affinity of the full-length protein for hcDNA does not decrease in the presence of magnesium. Experiments including a comparison of HMGB1 binding to hcDNA and to minicircles containing the CA/TG sequence, binding studies with HMGB1 mutated at intercalating amino acid residues (involved in recognition of distorted DNA structures), and exonuclease III footprinting, strongly suggest that the hemicatenane, not the DNA loop, is the main determinant of the affinity of HMGB1 for hcDNA. Experiments with supercoiled CA/TG-minicircles did not reveal any involvement of left-handed Z-DNA in HMGB1 binding. Our results point to a tight structural fit between HMGB1 and DNA hemicatenanes under physiological conditions, and suggest that one of the nuclear functions of HMGB1 could be linked to the possible presence of hemicatenanes in the cell.

(TG/CA)n repeats in human gene families: abundance and selective patterns of distribution according to function and gene length

BACKGROUND

Creation of human gene families was facilitated significantly by gene duplication and diversification. The (TG/CA)n repeats exhibit length variability, display genome-wide distribution, and are abundant in the human genome. Accumulation of evidences for their multiple functional roles including regulation of transcription and stimulation of recombination and splicing elect them as functional elements. Here, we report analysis of the distribution of (TG/CA)n repeats in human gene families.

RESULTS

The 1,317 human gene families were classified into six functional classes. Distribution of (TG/CA)n repeats were analyzed both from a global perspective and from a stratified perspective based on their biological properties. The number of genes with repeats decreased with increasing repeat length and several genes (53%) had repeats of multiple types in various combinations. Repeats were positively associated with the class of Signaling and communication whereas, they were negatively associated with the classes of Immune and related functions and of Information. The proportion of genes with (TG/CA)n repeats in each class was proportional to the corresponding average gene length. The repeat distribution pattern in large gene families generally mirrored the global distribution pattern but differed particularly for Collagen gene family, which was rich in repeats. The position and flanking sequences of the repeats of Collagen genes showed high conservation in the Chimpanzee genome. However the majority of these repeats displayed length polymorphism.

CONCLUSION

Positive association of repeats with genes of Signaling and communication points to their role in modulation of transcription. Negative association of repeats in genes of Information relates to the smaller gene length, higher expression and fundamental role in cellular physiology. In genes of Immune and related functions negative association of repeats perhaps relates to the smaller gene length and the directional nature of the recombinogenic processes to generate immune diversity. Thus, multiple factors including gene length, function and directionality of recombinogenic processes steered the observed distribution of (TG/CA)n repeats. Furthermore, the distribution of repeat patterns is consistent with the current model that long repeats tend to contract more than expand whereas, the reverse dynamics operates in short repeats.

Longer dinucleotide repeat polymorphism in matrix metalloproteinase-9 (MMP-9) gene promoter which correlates with higher HTLV-I Tax mediated transcriptional activity influences the risk of HTLV-I associated myelopathy/tropical spastic paraparesis (HAM/TSP)

Matrix metalloproteinase-9 (MMP-9) has been reported to be expressed in various inflammatory disorders including human T cell lymphotropic virus type I (HTLV-I) associated myelopathy/tropical spastic paraparesis (HAM/TSP). HTLV-I-infected T-cells expressed high levels of MMP-9 via viral transactivator Tax mediated activation of the MMP-9 promoter. To investigate whether the d(CA) repeat polymorphism in MMP-9 promoter affects the risk of developing HAM/TSP, we compared the allele frequencies between 200 HAM/TSP patients and 200 HTLV-I seropositive asymptomatic carriers (HCs). The longer d(CA) repeat alleles of MMP-9 promoter, which was associated with higher Tax-mediated transcriptional activity, was more frequently observed in HAM/TSP patients than HCs (p<0.01 by Mann-Whitney U-test). The length alteration of this d(CA) repeat in the MMP-9 promoter may cause phenotypic differences among HTLV-I infected infiltrating cells and may thereby be in part responsible for the development of HAM/TSP.

Interaction of the Zalpha domain of human ADAR1 with a negatively supercoiled plasmid visualized by atomic force microscopy

Interest to the left-handed DNA conformation has been recently boosted by the findings that a number of proteins contain the Zalpha domain, which has been shown to specifically recognize Z-DNA. The biological function of Zalpha is presently unknown, but it has been suggested that it may specifically direct protein regions of Z-DNA induced by negative supercoiling in actively transcribing genes. Many studies, including a crystal structure in complex with Z-DNA, have focused on the human ADAR1 Zalpha domain in isolation. We have hypothesized that the recognition of a Z-DNA sequence by the Zalpha(ADAR1) domain is context specific, occurring under energetic conditions, which favor Z-DNA formation. To test this hypothesis, we have applied atomic force microscopy to image Zalpha(ADAR1) complexed with supercoiled plasmid DNAs. We have demonstrated that the Zalpha(ADAR1) binds specifically to Z-DNA and preferentially to d(CG)(n) inserts, which require less energy for Z-DNA induction compared to other sequences. A notable finding is that site-specific Zalpha binding to d(GC)(13) or d(GC)(2)C(GC)(10) inserts is observed when DNA supercoiling is insufficient to induce Z-DNA formation. These results indicate that Zalpha(ADAR1) binding facilities the B-to-Z transition and provides additional support to the model that Z-DNA binding proteins may regulate biological processes through structure-specific recognition.

The Excitement of Discovery

I had the good luck to start research at the dawn of molecular biology when it was possible to ask fundamental questions about the nature of the nucleic acids and how information is transferred in living systems. The search for answers led me into many different areas, often with the question of how molecular structure leads to biological function. Early work in this period provided some of the roots supporting the current explosive developments in life sciences. Here I give a brief account of my development, describe some contributions, and provide a hint of the exhilaration in discovering new things. Most of all, I had the good fortune to have inspiring teachers, stimulating colleagues, and excellent students.

(TG:CA)(n) repeats in human housekeeping genes

The unravelling of human genome sequence gives a new opportunity to investigate the role of repetitive sequences in gene regulation. Among the various types of repetitive sequences, the dinucleotide (TG:CA)(n) repeats are one of the most abundant in human genome and exhibit polymorphism. Early on, it was observed that the (TG:CA)(n) repeats could modulate gene expression and has the propensity to undergo conformational transitions in in vivo conditions. Recent reports describe the role of polymorphic (TG:CA)(n) repeats in gene regulation in several genes. In this work, we have analysed the distribution of (TG:CA)(n) (n >or= 6) repeats in human 'housekeeping genes' on which recently released Gene Chip data is available. Our results indicate that (i). The number of short intragenic (TG:CA)(n) repeats is significantly higher than the number of long repeats (ii). the proportion of genes with (TG:CA)(n) repeats (n >or= 12 units) had lower mean expression levels compared to those without these repeats, (iii). the genes belonging to the functional class of 'signalling and communication' had a positive association with repeats in contrast to the genes belonging to the 'information' class that were negatively associated with repeats.

Timeline: Z-DNA: the long road to biological function

Biologists were puzzled by the discovery of left-handed Z-DNA because it seemed unnecessary. Z-DNA was stabilized by the negative supercoiling generated by transcription, which indicated a transient localized conformational change. Few laboratories worked on the biology of Z-DNA. However, the discovery that certain classes of proteins bound to Z-DNA with high affinity and great specificity indicated a biological role. The most recent data show that some of these proteins participate in the pathology of poxviruses.

Genetic polymorphisms of matrix metalloproteinases: functional importance in the development of chronic obstructive pulmonary disease?

Chronic obstructive pulmonary disease (COPD) is a major cause of global morbidity and mortality. COPD usually arises from an interaction between both host and environmental risk factors. Cigarette smoking is the major known environmental risk factor for the development of COPD, however, only a minority of smokers (approximately 15 to 20%) develop symptoms. COPD is known to cluster in families, which suggests that there is a genetic predisposition to airflow obstruction. Many candidate genes have been assessed, but the data are often unclear. Here we review evidence that genetic polymorphisms in matrix metalloproteinase genes MMP1, MMP9 and MMP12 may be important in the development of COPD. In a Caucasian population, polymorphisms in the MMP1 and MMP12 genes, but not MMP9, have been suggested to be either causative factors in smoking-related lung injury or are in linkage disequilibrium with other causative polymorphisms. Another study found an association between an MMP9 polymorphism and the development of smoking-induced pulmonary emphysema in Japanese smokers. Understanding the role of genetic polymorphisms in MMP1, MMP9 and MMP12 may help in the discovery of new and more effective therapies.

Structural analysis of hemicatenated DNA loops

BACKGROUND

We have previously isolated a stable alternative DNA structure, which was formed in vitro by reassociation of the strands of DNA fragments containing a 62 bp tract of the CA-microsatellite poly(CA).poly(TG). In the model which was proposed for this structure the double helix is folded into a loop, the base of the loop consists of a DNA junction in which one of the strands of one duplex passes between the two strands of the other duplex, forming a DNA hemicatenane in a hemiknot structure. The hemiknot DNA structures obtained with long CA/TG inserts have been imaged by AFM allowing us to directly visualize the loops.

RESULTS

Here we have analyzed this structure with several different techniques: high-resolution gel electrophoresis, probing by digestion with single stranded DNA-specific nucleases or with DNase I, modification with chemicals specific for unpaired bases, and atomic force microscopy. The data show a change in DNA structure localized to the CA/TG sequence and allow us to better understand the structure of this alternative conformation and the mechanism of its formation.

CONCLUSIONS

The present work is in good agreement with the model of hemicatenated DNA loop proposed previously. In the presence of protein HMGB1, shifted reassociation of the strands of DNA fragments containing a tract of the poly(CA).poly(TG) microsatellite leads to the formation of DNA loops maintained at their base by a hemicatenated junction located within the repetitive sequence. No mobility of the junction along the DNA molecule could be detected under the conditions used. The novel possibility to prepare DNA hemicatenanes should be useful to further study this alternative DNA structure and its involvement in replication or recombination.

Frequent oligonucleotide motifs in genomes of three streptococci

Complete genomes of three closely related Gram-positive bacteria Streptococcus pyogenes, Streptococcus pneumoniae and Lactococcus lactis are analyzed for abundances of short DNA sequence motifs (frequent words). The character and extent of frequent words are strikingly different among these genomes. The frequent words of S.pneumoniae split into three categories: parts of the previously characterized RUP and BOX repetitive elements and a 24 bp tandem repeat in the gene SP1772. The most abundant frequent words of L.lactis are all related to the 13 bp motif, WWNTTACTGACRR or its inverted complement YYGTCAGTAANWW. Distributional analysis of this motif, which we called highly repetitive motif (HRM), indicates its possible dual role. Frequent occurrences immediately downstream of genes suggest a possible role in transcription termination whereas spacings of consecutive HRMs consistent with the DNA helical period are indicative of a protein-binding site. Two regions of the L.lactis genome feature an intriguing pattern of several periodically occurring HRMs separated by precisely 59 bp. In a striking contrast to S.pneumoniae and L.lactis, S.pyogenes contains hardly any frequent words.

Microsatellite variation associated with prolactin expression and growth of salt-challenged tilapia

Biologists have long argued that runs of alternating purines and pyrimidines could form alternative DNA structures, which might regulate transcription. Here, we report that simple sequence repeat polymorphisms in the tilapia prolactin 1 (prl 1) promoter are associated with differences in prl 1 gene expression and the growth response of salt-challenged fishes. Individuals homozygous for long microsatellite alleles express less prl 1 in fresh water but more prl 1 in half-seawater than fishes with other genotypes. Our work provides the first in vivo evidence that differences in microsatellite length among individuals may indeed affect gene expression and that variance in expression has concomitant physiological consequences. These results suggest that dinucleotide microsatellites represent an under-appreciated source of genetic variation for regulatory evolution.

DNA loops and semicatenated DNA junctions

BACKGROUND

Alternative DNA conformations are of particular interest as potential signals to mark important sites on the genome. The structural variability of CA microsatellites is particularly pronounced; these are repetitive poly(CA). poly(TG) DNA sequences spread in all eukaryotic genomes as tracts of up to 60 base pairs long. Many in vitro studies have shown that the structure of poly(CA). poly(TG) can vary markedly from the classical right handed DNA double helix and adopt diverse alternative conformations. Here we have studied the mechanism of formation and the structure of an alternative DNA structure, named Form X, which was observed previously by polyacrylamide gel electrophoresis of DNA fragments containing a tract of the CA microsatellite poly(CA). poly(TG) but had not yet been characterized.

RESULTS

Formation of Form X was found to occur upon reassociation of the strands of a DNA fragment containing a tract of poly(CA). poly(TG), in a process strongly stimulated by the nuclear proteins HMG1 and HMG2. By inserting Form X into DNA minicircles, we show that the DNA strands do not run fully side by side but instead form a DNA knot. When present in a closed DNA molecule, Form X becomes resistant to heating to 100 degrees C and to alkaline pH.

CONCLUSIONS

Our data strongly support a model of Form X consisting in a DNA loop at the base of which the two DNA duplexes cross, with one of the strands of one duplex passing between the strands of the other duplex, and reciprocally, to form a semicatenated DNA junction also called a DNA hemicatenane.

Conservation within artiodactyls of an AATA interrupt in the IGF-I microsatellite for 19-35 million years

Occurrence of an AATA interrupt in the IGF-I microsatellite was investigated in a number of Artiodactyl species, namely pigs, camels, deer, cattle, goats, and sheep. Comparison of DNA sequences in the 5' flank of the microsatellite in these species revealed that the interrupt within the microsatellite is conserved in deer, cattle, sheep, and goats but is absent from camels and pigs. The interrupt was introduced into the Artiodactyl phylogeny after the divergence of the Camelidae but before the divergence of the Cervidae, and thus its time of origin can be estimated to be 19-35 MYA. In contrast to the repeat units which are hypermutable, the interrupt has been conserved for a very long time and may even have suppressed microsatellite variation by inhibiting replication slippage. A 12-bp deletion in the 5' flank of the microsatellite in camels corresponds to a consensus reversed repeat in deer, cattle, sheep, and goats with unknown functional significance. Apart from this deletion, the 5' flank of the microsatellite is highly conserved in Artiodactyl species.

Functional polymorphism in the matrix metalloproteinase-9 promoter as a potential risk factor for intracranial aneurysm

BACKGROUND AND PURPOSE

There is convincing evidence that susceptibility to intracranial aneurysms (ICAs) has a genetic component. However, few studies have sought to identify functional variation in specific candidate genes that may predispose individuals to develop an ICA.

METHODS

ICA cases and controls were genotyped for a simple length polymorphism in the promoter of matrix metalloproteinase-9 (MMP-9) to test for association between variation in the promoter and the occurrence of ICA. Alternative alleles were cloned into an in vitro reporter vector, transfected into human HT1080 fibroblasts, and assayed for promoter activity by beta-gal and luciferase assays. Electrophoretic gel shift assays were used to assess nuclear factor binding.

RESULTS

A length polymorphism in the promoter of MMP-9 was nonrandomly associated with the occurrence of ICA in a case-control study. This polymorphism was shown, by direct sequencing of 36 individuals, to be the only sequence variation within a 736-base pair region proximal to the transcriptional start site of the gene. Variation in the length of this repetitive element was shown to modulate promoter activity in an in vitro reporter assay, with the highest promoter activity being observed in constructs bearing the longest [(CA)23] element. Electrophoretic mobility shift assays were used to show that the (CA) element is bound by a sequence-specific DNA-binding protein.

CONCLUSIONS

Genetic variation in the promoter of the MMP-9 gene results in variation in its expression at the level of transcription. This may result in subtle differences in MMP-9 activity within the circle of Willis, leading to increased susceptibility to ICA formation.

Transcriptional activation of vesicular monoamine transporter 2 in the pre-B cell line Ea3.123

Uptake and storage of monoamines in secretory granules is accomplished by vesicular monoamine transporters, and it is likely that vesicular monoamine transporter 2 (VMAT2) is important for histamine transport in vivo. In the present study we have used the pre-B-cell line Ea3.123 to investigate the mechanisms involved in the transcriptional activation of the VMAT2 gene. In Ea3.123 cells, VMAT2 mRNA abundance was increased following mobilization of intracellular calcium, and this increased mRNA expression was paralleled by changes in l-histidine decarboxylase mRNA, suggesting that VMAT2 may be responsible for sequestration of histamine into secretory vesicles in this cell line. We cloned the 5'-flanking region of the VMAT2 gene and determined its transcriptional start site by primer extension of rat VMAT2 mRNA. There was no TATA or TATA-like sequence upstream of this region; instead there were GC-rich elements, Ca2+/cAMP-response-element- and SP1-binding motifs. Approx. 900 bp upstream of the transcriptional start site was a purine-pyrimidine repeat sequence that may form a Z-DNA structure. A series of 5'-deletional VMAT2-promoter segments cloned upstream of a luciferase reporter were capable of driving transcription and indicated the presence of multiple regulatory elements, while stimulation with ionomycin or PMA resulted in an increased level of the transcriptional activity of the 5'-promoter segments studied.

Conversions of the left-handed form and the protonated form of DNA back to the bound right-handed form by sanguinarine and ethidium: A comparative study

The interaction of sanguinarine and ethidium with right-handed (B-form), left-handed (Z-form) and left-handed protonated (designated as H(L)-form) structures of poly(dG-dC).poly(dG-dC) and poly(dG-me5dC).poly(dG-me5dC) was investigated by measuring the circular dichroism and UV absorption spectral analysis. Both sanguinarine and ethidium bind strongly to the B-form DNA and convert the Z-form and the H(L)-form back to the bound right-handed form. Circular dichroic data also show that the conformation at the binding site is right-handed, even though adjacent regions of the polymer have a left-handed conformation either in Z-form or in H(L)-form. Both the rate and extent of B-form to Z-form transition were decreased by sanguinarine and ethidium under ionic conditions that otherwise favour the left-handed conformation of the polynucleotides. The rate of decrease is faster in the case of ethidium as compared to that of sanguinarine. Scatchard analysis of the spectrophotometric data shows that sanguinarine binds strongly to both the polynucleotides in a non-cooperative manner under B-form conditions, in sharp contrast to the highly-cooperative binding under Z-form and H(L)-form conditions. Correlation of binding isotherms with circular dichroism data indicates that the cooperative binding of sanguinarine under the Z-form and the H(L)-form conditions is associated with a sequential conversion of the polymer from a left-handed to a bound right-handed conformation. Determination of bound alkaloid concentration by spectroscopic titration technique and the measurement of circular dichroic spectra have enabled us to calculate the number of base pairs of Z-form and H(L)-form that adopt a right-handed conformation for each bound alkaloid. Analysis reveals that 2-3 base pairs (bp) of Z-form of poly(dG-dC).poly(dG-dC) and poly(dG-me5dC).poly(dG-me5dC) switch to the right-handed form for each bound sanguinarine, while approximately same number of base pairs switch to the bound right-handed form in complexes with H(L)-form of these polynucleotides. Comparative binding analysis shows that ethidium also converts approximately 2 bp of Z-form or H(L)-form to bound right-handed form under same experimental conditions. Since sanguinarine binds preferentially to alternating GC sequences, which are capable of undergoing the B to Z or B to H(L) transition, these effects may be an important part in understanding its extensive biological activities.

Transcriptional activation of vesicular monoamine transporter 2 in the pre-B cell line Ea3.123

Uptake and storage of monoamines in secretory granules is accomplished by vesicular monoamine transporters, and it is likely that vesicular monoamine transporter 2 (VMAT2) is important for histamine transport in vivo. In the present study we have used the pre-B-cell line Ea3.123 to investigate the mechanisms involved in the transcriptional activation of the VMAT2 gene. In Ea3.123 cells, VMAT2 mRNA abundance was increased following mobilization of intracellular calcium, and this increased mRNA expression was paralleled by changes in l-histidine decarboxylase mRNA, suggesting that VMAT2 may be responsible for sequestration of histamine into secretory vesicles in this cell line. We cloned the 5'-flanking region of the VMAT2 gene and determined its transcriptional start site by primer extension of rat VMAT2 mRNA. There was no TATA or TATA-like sequence upstream of this region; instead there were GC-rich elements, Ca2+/cAMP-response-element- and SP1-binding motifs. Approx. 900 bp upstream of the transcriptional start site was a purine-pyrimidine repeat sequence that may form a Z-DNA structure. A series of 5'-deletional VMAT2-promoter segments cloned upstream of a luciferase reporter were capable of driving transcription and indicated the presence of multiple regulatory elements, while stimulation with ionomycin or PMA resulted in an increased level of the transcriptional activity of the 5'-promoter segments studied.

DNA unwinding in the CYC1 and DED1 yeast promoters

The capacity of promoter DNA of two yeast genes to be unwound was studied. Both promoters, those of the CYC1 and DED1 genes, contain long oligopurine.oligopyrimidine (R.Y) tracts. The two promoters were cloned into negatively supercoiled plasmids, and their sensitivity to single-strand specific nuclease P1 was examined. Extensive P1 cleavage was located within the R.Y tracts, and cleavage sites were mapped. The extent of cleavage was only slightly dependent on P1 concentration, indicating a slow conversion of an intermediate form of DNA into the P1 reactive state. The cleavage required negative supercoiling and was suppressed by NaCl, MgCl2 and spermine. Two-dimensional topoisomer analysis showed that six superhelical turns were opened in the plasmids examined. The results indicate that at sufficient torsional stress, the R.Y tracts can intermittently undergo a transition into an unwound, ready-to-separate state. The oligopurine.oligopyrimidine tracts may thus serve as DNA unwinding centers in the gene promoters where they reside.

A genome-derived (gaa.ttc)24 trinucleotide block binds nuclear protein(s) specifically and forms triple helices

The properties of simple trinucleotide repeats generate increased interest as expansions of certain trinucleotide blocks cause human diseases. Here, we studied protein binding and structural features of a perfect (gaa.ttc)24 tract in its original genomic environment. Electrophoretic mobility shift assays revealed that HeLa nuclear proteins bind to the DNA fragment containing the (gaa.ttc)24 block. Competition experiments using simple (gt.ac)n repeats differing in length and flanking regions showed no cross-reactivity with the major retarded band. For the specific (gaa. ttc)n/protein complex, a binding constant of 9.3x10-9 mol/l was determined. DNase I footprinting revealed protein binding sites located exclusively within the repeat with a preference for the (gaa)24 strand. OsO4 and DEPC modifications followed by electrophoretic and electron microscopical analyses showed that the (gaa.ttc)24 block forms different types of intramolecular triple helices: Under superhelical stress, different H-DNA isomers are evident, whereas exclusively H-Y forms were detected in the relaxed state. Together, these data have functional implications for genomic (gaa.ttc)n tracts.

The identification and characterization of microsatellites in the compact genome of the Japanese pufferfish, Fugu rubripes: perspectives in functional and comparative genomic analyses

Fugu rubripes (Fugu) has one of the smallest recorded vertebrate genomes and is an economic tool for comparative DNA sequence analysis. Initial characterization of 128 kb of Fugu DNA attributed the compactness of this genome, in part, to a sparseness of repetitive DNA sequence compared with mammalian genomic sequences. This paper describes a new and comprehensive analysis in which 501 theoretically possible microsatellites with a repeat unit of one to six bases were used to query two orders of magnitude more Fugu DNA (i.e. 11.338 Mb). A total of 6042 microsatellites were identified and categorized. In decreasing order, the 20 most frequently occurring microsatellites are AC, A, C, AGG, AG, AGC, AAT, AAAT, ACAG, ACGC, ATCC, AAC, ATC, AGGG, AAAG, AAG, AAAC, AT, CCG and TTAGGG. The 20 most frequently occurring microsatellites represent 81.79% of all microsatellites identified. Our results indicate that one microsatellite occurs every 1.876 kb of DNA in Fugu, 11.55% of the microsatellites are detected in open reading frames that are predicted protein coding regions. With respect to the proportion of microsatellites present in open reading frames and the total abundance (bp) of all microsatellites, the genome of Fugu is similar to the genome of many other vertebrate species. Previous estimates performed indicate that approximately 1% of many vertebrate genomes are comprized of microsatellite sequences. However, many differences prevail in the abundance and frequency of the individual microsatellite classes. Many of the frequently occurring microsatellites in Fugu are known to code in other species for regions in proteins such as transcription factors, whilst others are associated with known functions, such as transcription factor binding sites and form part of promoter regions in DNA sequences of genes. Therefore, it is likely that such repeats in genomes have a role in the evolution of genes, regulation of gene expression and consequently the evolution of species.

Homology-directed repair is a major double-strand break repair pathway in mammalian cells

Mammalian cells have been presumed to repair potentially lethal chromosomal double-strand breaks (DSBs) in large part by processes that do not require homology to the break site. This contrasts with Saccharomyces cerevisiae where the major DSB repair pathway is homologous recombination. Recently, it has been determined that DSBs in genomic DNA in mammalian cells can stimulate homologous recombination as much as 3 or 4 orders of magnitude, suggesting that homology-directed repair may play an important role in the repair of chromosomal breaks. To determine whether mammalian cells use recombinational repair at a significant level, we have analyzed the spectrum of repair events at a defined chromosomal break by using direct physical analysis of repair products. When an endonuclease-generated DSB is introduced into one of two direct repeats, homologous repair is found to account for 30-50% of observed repair events. Both noncrossover and deletional homologous repair products are detected, at approximately a 1:3 ratio. These results demonstrate the importance of homologous recombination in the repair of DSBs in mammalian cells. In the remaining observed repair events, DSBs are repaired by nonhomologous processes. The nonhomologous repair events generally result in small deletions or insertions at the break site, although a small fraction of events result in larger chromosomal rearrangements. Interestingly, in two insertions, GT repeats were integrated at one of the broken chromosome ends, suggesting that DSB repair can contribute to the spread of microsatellite sequences in mammalian genomes.

Construction of a Z-DNA-specific restriction endonuclease

Novel restriction enzymes can be created by fusing the nuclease domain of FokI endonuclease with defined DNA binding domains. Recently, we have characterized a domain (Z alpha) from the N-terminal region of human double-stranded RNA adenosine deaminase (hADAR1), which binds the Z-conformation with high specificity. Here we report creation of a conformation-specific endonuclease, Z alpha nuclease, which is a chimera of Z alpha and FokI nuclease. Purified Z alpha nuclease cleaves negatively supercoiled plasmids only when they contain a Z-DNA forming insert, such as (dC-dG)13. The precise location of the cleavage sites was determined by primer extension. Cutting has been mapped to the edge of the B-Z junction, suggesting that Z alpha nuclease binds within the Z-DNA insert, but cleaves in the nearby B-DNA, by using a mechanism similar to type IIs restriction enzymes. These data show that Z alpha binds Z-DNA in an environment similar to that in a cell. Z alpha nuclease, a structure-specific restriction enzyme, may be a useful tool for further study of the biological role of Z-DNA.

Characterization of S1 nuclease sensitive site at transcription initiation region of Attacus ricini rDNA

A single-stranded S1 nuclease hypersensitive site which contains a d(AT)(18) sequence structure located in the 5'-non transcription spacer of silkworm A. ricini ribosomal RNA gene has been reported([1]). Using starved-refed silkworms, another S1 nuclease sensitive site was found existing in the rDNA chromatin, while under merely starving, this S1 sensitive site disappeared([2]). Recently this inducible S1 sensitive site has been further determined. It consists of a d(GT)(10)-d(AT)(10) special DNA sequence at the transcription initiation region, and shows a behavior of ease in DNA-unwinding, indicating that S1 nuclease sensitive sites may have an important function in the regulation of rDNA transcription and replication.

Analysis of left-handed Z-DNA formation in short d(CG)n sequences in Escherichia coli and Halobacterium halobium plasmids. Stabilization by increasing repeat length and DNA supercoiling but not salinity

To evaluate the relative importance of alternating d(CG) sequence length, DNA supercoiling, and salt in left-handed Z-DNA formation, plasmids containing short d(CG)n sequences (n = 3-17) with the capability of replicating in either Escherichia coli or the halophilic archaeum Halobacterium halobium were constructed. Z-DNA conformation in the d(CG)n sequences was assayed by (i) a band shift assay using the Z-DNA-specific Z22 monoclonal antibody (ZIBS assay); (ii) an S1 nuclease cleavage-primer extension assay to map B-Z junctions; and (iii) a BssHII restriction inhibition assay. Using the ZIBS assay on plasmids purified from E. coli, the transition from B-DNA to Z-DNA occurred from d(CG)4, to d(CG)5, with 20% of d(CG)4, and 90% of d(CG)5 in Z-DNA conformation. These findings were consistent with the results of S1 nuclease cleavage observed at B-Z junctions flanking d(CG)4 and d(CG)5 sequences. Resistance to BssHII restriction endonuclease digestion was observed only in supercoiled plasmids containing d(CG)8 or longer sequences, indicating that shorter d(CG)n sequences are in dynamic equilibrium between B- and Z-DNA conformations. When a plasmid containing d(CG)4, was isolated from a topA mutant of E. coli, it contained 25% greater linking deficiency and 40% greater Z-DNA conformation in the alternating d(CG) region. In plasmids purified from H. halobium, which showed 30% greater linking deficiency than from E. coli, 20-40% greater Z-DNA formation was found in d(CG)4-6 sequences. Surprisingly, no significant difference in Z-DNA formation could be detected in d(CG)3-17 sequences in plasmids from either E. coli or H. halobium in the NaCl concentration range of 0.1-4 M.