Inverted repeats, stem-loops, and cruciforms: significance for initiation of DNA replication

Decoding the Structural Dynamics and Conformational Alternations of DNA Secondary Structures by Single-Molecule FRET Microspectroscopy

In addition to the canonical double helix form, DNA is known to be extrapolated into several other secondary structural patterns involving themselves in inter- and intramolecular type hydrogen bonding. The secondary structures of nucleic acids go through several stages of multiple, complex, and interconvertible heterogeneous conformations. The journey of DNA through these conformers has significant importance and has been monitored thoroughly to establish qualitative and quantitative information about the transition between the unfolded, folded, misfolded, and partially folded states. During this structural interconversion, there always exist specific populations of intermediates, which are short-lived or sometimes even do not accumulate within a heterogeneous population and are challenging to characterize using conventional ensemble techniques. The single-molecule FRET(sm-FRET) microspectroscopic method has the advantages to overcome these limitations and monitors biological phenomena transpiring at a measurable high rate and balanced stochastically over time. Thus, tracing the time trajectory of a particular molecule enables direct measurement of the rate constant of each transition step, including the intermediates that are hidden in the ensemble level due to their low concentrations. This review is focused on the advantages of the employment of single-molecule Forster's resonance energy transfer (sm-FRET), which is worthwhile to access the dynamic architecture and structural transition of various secondary structures that DNA adopts, without letting the donor of one molecule to cross-talk with the acceptor of any other. We have emphasized the studies performed to explore the states of folding and unfolding of several nucleic acid secondary structures, for example, the DNA hairpin, Holliday junction, G-quadruplex, and i-motif.

A Cytotoxic Bis(1,2,3-triazol-5-ylidene)carbazolide Gold(III) Complex Targets DNA by Partial Intercalation

The syntheses of bis(triazolium)carbazole precursors and their corresponding coinage metal (Au, Ag) complexes are reported. For alkylated triazolium salts, di- or tetranuclear complexes with bridging ligands were isolated, while the bis(aryl) analogue afforded a bis(carbene) AuI -CNC pincer complex suitable for oxidation to the redox-stable [AuIII (CNC)Cl]+ cation. Although the ligand salt and the [AuIII (CNC)Cl]+ complex were both notably cytotoxic toward the breast cancer cell line MDA-MB-231, the AuIII complex was somewhat more selective. Electrophoresis, viscometry, UV-vis, CD and LD spectroscopy suggest the cytotoxic [AuIII (CNC)Cl]+ complex behaves as a partial DNA intercalator. In silico screening indicated that the [AuIII (CNC)Cl]+ complex can target DNA three-way junctions with good specificity, several other regular B-DNA forms, and Z-DNA. Multiple hydrophobic π-type interactions involving T and A bases appear to be important for B-form DNA binding, while phosphate O⋅⋅⋅Au interactions evidently underpin Z-DNA binding. The CNC ligand effectively stabilizes the AuIII ion, preventing reduction in the presence of glutathione. Both the redox stability and DNA affinity of the hit compound might be key factors underpinning its cytotoxicity in vitro.

Cruciform Formable Sequences within Pou5f1 Enhancer Are Indispensable for Mouse ES Cell Integrity

DNA can adopt various structures besides the B-form. Among them, cruciform structures are formed on inverted repeat (IR) sequences. While cruciform formable IRs (CFIRs) are sometimes found in regulatory regions of transcription, their function in transcription remains elusive, especially in eukaryotes. We found a cluster of CFIRs within the mouse Pou5f1 enhancer. Here, we demonstrate that this cluster or some member(s) plays an active role in the transcriptional regulation of not only Pou5f1, but also Sox2, Nanog, Klf4 and Esrrb. To clarify in vivo function of the cluster, we performed genome editing using mouse ES cells, in which each of the CFIRs was altered to the corresponding mirror repeat sequence. The alterations reduced the level of the Pou5f1 transcript in the genome-edited cell lines, and elevated those of Sox2, Nanog, Klf4 and Esrrb. Furthermore, transcription of non-coding RNAs (ncRNAs) within the enhancer was also upregulated in the genome-edited cell lines, in a similar manner to Sox2, Nanog, Klf4 and Esrrb. These ncRNAs are hypothesized to control the expression of these four pluripotency genes. The CFIRs present in the Pou5f1 enhancer seem to be important to maintain the integrity of ES cells.

IUPACpal: efficient identification of inverted repeats in IUPAC-encoded DNA sequences

BACKGROUND

An inverted repeat is a DNA sequence followed downstream by its reverse complement, potentially with a gap in the centre. Inverted repeats are found in both prokaryotic and eukaryotic genomes and they have been linked with countless possible functions. Many international consortia provide a comprehensive description of common genetic variation making alternative sequence representations, such as IUPAC encoding, necessary for leveraging the full potential of such broad variation datasets.

RESULTS

We present IUPACPAL, an exact tool for efficient identification of inverted repeats in IUPAC-encoded DNA sequences allowing also for potential mismatches and gaps in the inverted repeats.

CONCLUSION

Within the parameters that were tested, our experimental results show that IUPACPAL compares favourably to a similar application packaged with EMBOSS. We show that IUPACPAL identifies many previously unidentified inverted repeats when compared with EMBOSS, and that this is also performed with orders of magnitude improved speed.

Integrating genome sequence and structural data for statistical learning to predict transcription factor binding sites

We report an approach to predict DNA specificity of the tetracycline repressor (TetR) family transcription regulators (TFRs). First, a genome sequence-based method was streamlined with quantitative P-values defined to filter out reliable predictions. Then, a framework was introduced to incorporate structural data and to train a statistical energy function to score the pairing between TFR and TFR binding site (TFBS) based on sequences. The predictions benchmarked against experiments, TFBSs for 29 out of 30 TFRs were correctly predicted by either the genome sequence-based or the statistical energy-based method. Using P-values or Z-scores as indicators, we estimate that 59.6% of TFRs are covered with relatively reliable predictions by at least one of the two methods, while only 28.7% are covered by the genome sequence-based method alone. Our approach predicts a large number of new TFBs which cannot be correctly retrieved from public databases such as FootprintDB. High-throughput experimental assays suggest that the statistical energy can model the TFBSs of a significant number of TFRs reliably. Thus the energy function may be applied to explore for new TFBSs in respective genomes. It is possible to extend our approach to other transcriptional factor families with sufficient structural information.

A reference catalog of DNA palindromes in the human genome and their variations in 1000 Genomes

A palindrome in DNA is like a palindrome in language, but when read backwards, it is a complement of the forward sequence; effectively, the two halves of a sequence complement each other from its midpoint like in a double strand of DNA. Palindromes are distributed throughout the human genome and play significant roles in gene expression and regulation. Palindromic mutations are linked to many human diseases, such as neuronal disorders, mental retardation, and various cancers. In this work, we computed and analyzed the palindromic sequences in the human genome and studied their conservation in personal genomes using 1000 Genomes data. We found that ~30% of the palindromes exhibit variation, some of which are caused by rare variants. The analysis of disease/trait-associated single-nucleotide polymorphisms in palindromic regions showed that disease-associated risk variants are 14 times more likely to be present in palindromic regions than in other regions. The catalog of palindromes in the reference genome and 1000 Genomes is being made available here with details on their variations in each individual genome to serve as a resource for future and retrospective whole-genome studies identifying statistically significant palindrome variations associated with diseases or traits and their roles in disease mechanisms.

Large inverted repeats identified by intra-specific comparison of mitochondrial genomes provide insights into the evolution of Agrocybe aegerita

Genomic structure and content of Agrocybe aegerita mitochondrial DNA contain essential information regarding the evolution of this gourmet mushroom. In this study, eight isolates of A. aegerita were sequenced and assembled into complete mitochondrial genomes. The mtDNA of the isolate Ag0067 contained two genotypes, both of which were quadripartite architecture consisting of two identical inverted repeats, separated by a small single-copy region and a large single-copy region. The only difference was opposite directions of the small single-copy region. The mtDNAs ranged from 116,329 bp to 134,035 bp, harboring two large identical inverted repeats. Genes of plasmid-origin were present in regions flanked by inverted repeat ID2. Most of the core genes evolved at a relatively low rate, whereas five tRNA genes located in corresponding regions of Ag0002:1-14000 and Ag0002:50001-61000 showed higher diversity. A long fragment inversion (10 Kb) was suggested to have occurred during the differentiation of two main clades, leading to two different gene orders. The number and distribution of the introns varied greatly among the A. aegerita mtDNAs. Fast invasion of short insertions likely resulted in the diversity of introns as well as other non-coding regions, increasing the variation of the mtDNAs. We raised a model about the evolution of the large repeats to explain the unusual features of A. aegerita mtDNAs. This study constructed quadripartite architecture of A. aegerita mtDNAs analogous to chloroplast DNA, proposed an interconversion model of the divergent mitochondrial genotypes with large inverted repeats. The findings could increase our knowledge of fungal evolution.

Evolutionary advantage of anti-parallel strand orientation of duplex DNA

DNA in all living systems shares common properties that are remarkably well suited to its function, suggesting refinement by evolution. However, DNA also shares some counter-intuitive properties which confer no obvious benefit, such as strand directionality and anti-parallel strand orientation, which together result in the complicated lagging strand replication. The evolutionary dynamics that led to these properties of DNA remain unknown but their universality suggests that they confer as yet unknown selective advantage to DNA. In this article, we identify an evolutionary advantage of anti-parallel strand orientation of duplex DNA, within a given set of plausible premises. The advantage stems from the increased rate of replication, achieved by dividing the DNA into predictable, independently and simultaneously replicating segments, as opposed to sequentially replicating the entire DNA, thereby parallelizing the replication process. We show that anti-parallel strand orientation is essential for such a replicative organization of DNA, given our premises, the most important of which is the assumption of the presence of sequence-dependent asymmetric cooperativity in DNA.

Polymorphic centromere locations in the pathogenic yeast Candida parapsilosis

Centromeres pose an evolutionary paradox: strongly conserved in function but rapidly changing in sequence and structure. However, in the absence of damage, centromere locations are usually conserved within a species. We report here that isolates of the pathogenic yeast species Candida parapsilosis show within-species polymorphism for the location of centromeres on two of its eight chromosomes. Its old centromeres have an inverted-repeat (IR) structure, whereas its new centromeres have no obvious structural features but are located within 30 kb of the old site. Centromeres can therefore move naturally from one chromosomal site to another, apparently spontaneously and in the absence of any significant changes in DNA sequence. Our observations are consistent with a model in which all centromeres are genetically determined, such as by the presence of short or long IRs or by the ability to form cruciforms. We also find that centromeres have been hotspots for genomic rearrangements in the C. parapsilosis clade.

DoriC 10.0: an updated database of replication origins in prokaryotic genomes including chromosomes and plasmids

DoriC, a database of replication origins, was initially created to present the bacterial oriCs predicted by Ori-Finder or determined by experiments in 2007. DoriC 5.0, an updated database of oriC regions in both bacterial and archaeal genomes, was published in the 2013 Nucleic Acids Research database issue. Now, the latest release DoriC 10, a large-scale update of replication origins in prokaryotic genomes including chromosomes and plasmids, has been presented with a completely redesigned user interface, which is freely available at http://tubic.org/doric/ and http://tubic.tju.edu.cn/doric/. In the current release, the database of DoriC has made significant improvements compared with version 5.0 as follows: (i) inclusion of oriCs on more bacterial chromosomes increased from 1633 to 7580; (ii) inclusion of oriCs on more archaeal chromosomes increased from 86 to 226; (iii) inclusion of 1209 plasmid replication origins retrieved from NCBI annotations or predicted by in silico analysis; (iv) inclusion of more replication origin elements on bacterial chromosomes including DnaA-trio motifs. Now, DoriC becomes the most complete and scalable database of replication origins in prokaryotic genomes, and facilitates the studies in large-scale oriC data mining, strand-biased analyses and replication origin predictions.

Intraspecific Variation within the Utricularia amethystina Species Morphotypes Based on Chloroplast Genomes

Utricularia amethystina Salzm. ex A.St.-Hil. & Girard (Lentibulariaceae) is a highly polymorphic carnivorous plant taxonomically rearranged many times throughout history. Herein, the complete chloroplast genomes (cpDNA) of three U. amethystina morphotypes: purple-, white-, and yellow-flowered, were sequenced, compared, and putative markers for systematic, populations, and evolutionary studies were uncovered. In addition, RNA-Seq and RNA-editing analysis were employed for functional cpDNA evaluation. The cpDNA of three U. amethystina morphotypes exhibits typical quadripartite structure. Fine-grained sequence comparison revealed a high degree of intraspecific genetic variability in all morphotypes, including an exclusive inversion in the psbM and petN genes in U. amethystina yellow. Phylogenetic analyses indicate that U. amethystina morphotypes are monophyletic. Furthermore, in contrast to the terrestrial Utricularia reniformis cpDNA, the U. amethystina morphotypes retain all the plastid NAD(P)H-dehydrogenase (ndh) complex genes. This observation supports the hypothesis that the ndhs in terrestrial Utricularia were independently lost and regained, also suggesting that different habitats (aquatic and terrestrial) are not related to the absence of Utricularia ndhs gene repertoire as previously assumed. Moreover, RNA-Seq analyses recovered similar patterns, including nonsynonymous RNA-editing sites (e.g., rps14 and petB). Collectively, our results bring new insights into the chloroplast genome architecture and evolution of the photosynthesis machinery in the Lentibulariaceae.

The Rich World of p53 DNA Binding Targets: The Role of DNA Structure

The tumor suppressor functions of p53 and its roles in regulating the cell cycle, apoptosis, senescence, and metabolism are accomplished mainly by its interactions with DNA. p53 works as a transcription factor for a significant number of genes. Most p53 target genes contain so-called p53 response elements in their promoters, consisting of 20 bp long canonical consensus sequences. Compared to other transcription factors, which usually bind to one concrete and clearly defined DNA target, the p53 consensus sequence is not strict, but contains two repeats of a 5′RRRCWWGYYY3′ sequence; therefore it varies remarkably among target genes. Moreover, p53 binds also to DNA fragments that at least partially and often completely lack this consensus sequence. p53 also binds with high affinity to a variety of non-B DNA structures including Holliday junctions, cruciform structures, quadruplex DNA, triplex DNA, DNA loops, bulged DNA, and hemicatenane DNA. In this review, we summarize information of the interactions of p53 with various DNA targets and discuss the functional consequences of the rich world of p53 DNA binding targets for its complex regulatory functions.

Non-B-Form DNA Is Enriched at Centromeres

Animal and plant centromeres are embedded in repetitive "satellite" DNA, but are thought to be epigenetically specified. To define genetic characteristics of centromeres, we surveyed satellite DNA from diverse eukaryotes and identified variation in <10-bp dyad symmetries predicted to adopt non-B-form conformations. Organisms lacking centromeric dyad symmetries had binding sites for sequence-specific DNA-binding proteins with DNA-bending activity. For example, human and mouse centromeres are depleted for dyad symmetries, but are enriched for non-B-form DNA and are associated with binding sites for the conserved DNA-binding protein CENP-B, which is required for artificial centromere function but is paradoxically nonessential. We also detected dyad symmetries and predicted non-B-form DNA structures at neocentromeres, which form at ectopic loci. We propose that centromeres form at non-B-form DNA because of dyad symmetries or are strengthened by sequence-specific DNA binding proteins. This may resolve the CENP-B paradox and provide a general basis for centromere specification.

A strong structural correlation between short inverted repeat sequences and the polyadenylation signal in yeast and nucleosome exclusion by these inverted repeats

DNA sequences that read the same from 5′ to 3′ in either strand are called inverted repeat sequences or simply IRs. They are found throughout a wide variety of genomes, from prokaryotes to eukaryotes. Despite extensive research, their in vivo functions, if any, remain unclear. Using Saccharomyces cerevisiae, we performed genome-wide analyses for the distribution, occurrence frequency, sequence characteristics and relevance to chromatin structure, for the IRs that reportedly have a cruciform-forming potential. Here, we provide the first comprehensive map of these IRs in the S. cerevisiae genome. The statistically significant enrichment of the IRs was found in the close vicinity of the DNA positions corresponding to polyadenylation [poly(A)] sites and ~ 30 to ~ 60 bp downstream of start codon-coding sites (referred to as ‘start codons’). In the former, ApT- or TpA-rich IRs and A-tract- or T-tract-rich IRs are enriched, while in the latter, different IRs are enriched. Furthermore, we found a strong structural correlation between the former IRs and the poly(A) signal. In the chromatin formed on the gene end regions, the majority of the IRs causes low nucleosome occupancy. The IRs in the region ~ 30 to ~ 60 bp downstream of start codons are located in the + 1 nucleosomes. In contrast, fewer IRs are present in the adjacent region downstream of start codons. The current study suggests that the IRs play similar roles in Escherichia coli and S. cerevisiae to regulate or complete transcription at the RNA level.

Complex Analyses of Short Inverted Repeats in All Sequenced Chloroplast DNAs

Chloroplasts are key organelles in the management of oxygen in algae and plants and are therefore crucial for all living beings that consume oxygen. Chloroplasts typically contain a circular DNA molecule with nucleus-independent replication and heredity. Using "palindrome analyser" we performed complete analyses of short inverted repeats (S-IRs) in all chloroplast DNAs (cpDNAs) available from the NCBI genome database. Our results provide basic parameters of cpDNAs including comparative information on localization, frequency, and differences in S-IR presence. In a total of 2,565 cpDNA sequences available, the average frequency of S-IRs in cpDNA genomes is 45 S-IRs/per kbp, significantly higher than that found in mitochondrial DNA sequences. The frequency of S-IRs in cpDNAs generally decreased with S-IR length, but not for S-IRs 15, 22, 24, or 27 bp long, which are significantly more abundant than S-IRs with other lengths. These results point to the importance of specific S-IRs in cpDNA genomes. Moreover, comparison by Levenshtein distance of S-IR similarities showed that a limited number of S-IR sequences are shared in the majority of cpDNAs. S-IRs are not located randomly in cpDNAs, but are length-dependently enriched in specific locations, including the repeat region, stem, introns, and tRNA regions. The highest enrichment was found for 12 bp and longer S-IRs in the stem-loop region followed by 12 bp and longer S-IRs located before the repeat region. On the other hand, S-IRs are relatively rare in rRNA sequences and around introns. These data show nonrandom and conserved arrangements of S-IRs in chloroplast genomes.

Increased Flexibility between Stems of Intramolecular Three-Way Junctions by the Insertion of Bulges

Intramolecular junctions are a ubiquitous structure within DNA and RNA; three-way junctions in particular have high strain around the junction because of the lack of flexibility, preventing the junctions from adopting conformations that would allow for optimal folding. In this work, we used a combination of calorimetric and spectroscopic techniques to study the unfolding of four intramolecular three-way junctions. The control three-way junction, 3H, has the sequence d(GAAATTGCGCT5GCGCGTGCT5GCACAATTTC), which has three arms of different sequences. We studied three other three-way junctions in which one (2HS1H), two (HS12HS1), and three (HS1HS1HS1) cytosine bulges were placed at the junction to allow the arms to adopt a wider range of conformations that may potentially relieve strain. Through calorimetric studies, it was concluded that bulges produce only minor effects on the enthalpic and thermal stability at physiological salt concentrations for 2HS1H and HS1HS1HS1. HS12HS1 displays the strongest effect, with the GTGC stem lacking a defined transition. In addition to unfolding thermodynamics, the differential binding of counterions, water, and protons was determined. It was found that with each bulge, there was a large increase in the binding of counterions; this correlated with a decrease in the immobilization of structural water molecules. The increase in counterion uptake upon folding likely displaces binding of structural water, which is measured by the osmotic stress method, in favor of electrostricted waters. The cytosine bulges do not affect the binding of protons; this finding indicates that the bulges are not forming base-triplet stacks. These results indicate that bulges in junctions do not affect the unfolding profile or the enthalpy of oligonucleotides but do affect the number and amount of molecules immobilized by the junction.

The Prevalence and Evolutionary Conservation of Inverted Repeats in Proteobacteria

Perfect short inverted repeats (IRs) are known to be enriched in a variety of bacterial and eukaryotic genomes. Currently, it is unclear whether perfect IRs are conserved over evolutionary time scales. In this study, we aimed to characterize the prevalence and evolutionary conservation of IRs across 20 proteobacterial strains. We first identified IRs in Escherichia coli K-12 substr MG1655 and showed that they are overabundant. We next aimed to test whether this overabundance is reflected in the conservation of IRs over evolutionary time scales. To this end, for each perfect IR identified in E. coli MG1655, we collected orthologous sequences from related proteobacterial genomes. We next quantified the evolutionary conservation of these IRs, that is, the presence of the exact same IR across orthologous regions. We observed high conservation of perfect IRs: out of the 234 examined orthologous regions, 145 were more conserved than expected, which is statistically significant even after correcting for multiple testing. Our results together with previous experimental findings support a model in which imperfect IRs are corrected to perfect IRs in a preferential manner via a template switching mechanism.

Requirement or exclusion of inverted repeat sequences with cruciform-forming potential in Escherichia coli revealed by genome-wide analyses

Inverted repeat (IR) sequences are DNA sequences that read the same from 5' to 3' in each strand. Some IRs can form cruciforms under the stress of negative supercoiling, and these IRs are widely found in genomes. However, their biological significance remains unclear. The aim of the current study is to explore this issue further. We constructed the first Escherichia coli genome-wide comprehensive map of IRs with cruciform-forming potential. Based on the map, we performed detailed and quantitative analyses. Here, we report that IRs with cruciform-forming potential are statistically enriched in the following five regions: the adjacent regions downstream of the stop codon-coding sites (referred to as the stop codons), on and around the positions corresponding to mRNA ends (referred to as the gene ends), ~ 20 to ~45 bp upstream of the start codon-coding sites (referred to as the start codons) within the 5'-UTR (untranslated region), ~ 25 to ~ 60 bp downstream of the start codons, and promoter regions. For the adjacent regions downstream of the stop codons and on and around the gene ends, most of the IRs with a repeat unit length of ≥ 8 bp and a spacer size of ≤ 8 bp were parts of the intrinsic terminators, regardless of the location, and presumably used for Rho-independent transcription termination. In contrast, fewer IRs were present in the small region preceding the start codons. In E. coli, IRs with cruciform-forming potential are actively placed or excluded in the regulatory regions for the initiation and termination of transcription and translation, indicating their deep involvement or influence in these processes.

JNSViewer-A JavaScript-based Nucleotide Sequence Viewer for DNA/RNA secondary structures

Many tools are available for visualizing RNA or DNA secondary structures, but there is scarce implementation in JavaScript that provides seamless integration with the increasingly popular web computational platforms. We have developed JNSViewer, a highly interactive web service, which is bundled with several popular tools for DNA/RNA secondary structure prediction and can provide precise and interactive correspondence among nucleotides, dot-bracket data, secondary structure graphs, and genic annotations. In JNSViewer, users can perform RNA secondary structure predictions with different programs and settings, add customized genic annotations in GFF format to structure graphs, search for specific linear motifs, and extract relevant structure graphs of sub-sequences. JNSViewer also allows users to choose a transcript or specific segment of Arabidopsis thaliana genome sequences and predict the corresponding secondary structure. Popular genome browsers (i.e., JBrowse and BrowserGenome) were integrated into JNSViewer to provide powerful visualizations of chromosomal locations, genic annotations, and secondary structures. In addition, we used StructureFold with default settings to predict some RNA structures for Arabidopsis by incorporating in vivo high-throughput RNA structure profiling data and stored the results in our web server, which might be a useful resource for RNA secondary structure studies in plants. JNSViewer is available at http://bioinfolab.miamioh.edu/jnsviewer/index.html.

Amplification of a transgene within a long array of replication origins favors higher gene expression in animal cells

Plasmids with both a mammalian replication initiation region (IR) and a matrix attachment region (MAR) are spontaneously amplified in transfected cells, and generate extrachromosomal double minute (DM) or chromosomal homogeneously staining region (HSR). We previously isolated the shortest core IR (G5) required for gene amplification. In this study, we ligated the G5 DNA to create direct or inverted repeats, mixed the repeats with an expression plasmid, and transfected the mixture into human COLO 320DM or hamster CHO DG44 cells. Consequently, we found that the transfected sequence generated DMs or HSR where, surprisingly, the plasmid sequence was embedded within a long stretch of G5 sequences. The amplified structure from the direct G5 repeats was stable, whereas that from the inverted repeats was not. The amplification might be explained by the efficient replication/multimerization of the G5 repeat and recombination with the co-transfected plasmid in an extrachromosomal context. The product might then be integrated into a chromosome arm to generate a HSR. The expression from the plasmid within the long G5 array was much higher than that from a simple plasmid repeat. Because G5 is a core IR that favors gene expression, a long array of G5 provides an excellent environment for gene expression from the embedded plasmid.

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.

The unique genomic landscape surrounding the EPSPS gene in glyphosate resistant Amaranthus palmeri: a repetitive path to resistance

BACKGROUND

The expanding number and global distributions of herbicide resistant weedy species threaten food, fuel, fiber and bioproduct sustainability and agroecosystem longevity. Amongst the most competitive weeds, Amaranthus palmeri S. Wats has rapidly evolved resistance to glyphosate primarily through massive amplification and insertion of the 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) gene across the genome. Increased EPSPS gene copy numbers results in higher titers of the EPSPS enzyme, the target of glyphosate, and confers resistance to glyphosate treatment. To understand the genomic unit and mechanism of EPSPS gene copy number proliferation, we developed and used a bacterial artificial chromosome (BAC) library from a highly resistant biotype to sequence the local genomic landscape flanking the EPSPS gene.

RESULTS

By sequencing overlapping BACs, a 297 kb sequence was generated, hereafter referred to as the "EPSPS cassette." This region included several putative genes, dense clusters of tandem and inverted repeats, putative helitron and autonomous replication sequences, and regulatory elements. Whole genome shotgun sequencing (WGS) of two biotypes exhibiting high and no resistance to glyphosate was performed to compare genomic representation across the EPSPS cassette. Mapping of sequences for both biotypes to the reference EPSPS cassette revealed significant differences in upstream and downstream sequences relative to EPSPS with regard to both repetitive units and coding content between these biotypes. The differences in sequence may have resulted from a compounded-building mechanism such as repetitive transpositional events. The association of putative helitron sequences with the cassette suggests a possible amplification and distribution mechanism. Flow cytometry revealed that the EPSPS cassette added measurable genomic content.

CONCLUSIONS

The adoption of glyphosate resistant cropping systems in major crops such as corn, soybean, cotton and canola coupled with excessive use of glyphosate herbicide has led to evolved glyphosate resistance in several important weeds. In Amaranthus palmeri, the amplification of the EPSPS cassette, characterized by a complex array of repetitive elements and putative helitron sequences, suggests an adaptive structural genomic mechanism that drives amplification and distribution around the genome. The added genomic content not found in glyphosate sensitive plants may be driving evolution through genome expansion.

A pilot study on the prevalence of DNA palindromes in breast cancer genomes

Background

DNA palindromes are a unique pattern of repeat sequences that are present in the human genome. It consists of a sequence of nucleotides in which the second half is the complement of the first half but appearing in reverse order. These palindromic sequences may have a significant role in DNA replication, transcription and gene regulation processes. They occur frequently in human cancers by clustering at specific locations of the genome that undergo gene amplification and tumorigenesis. Moreover, some studies showed that palindromes are clustered in amplified regions of breast cancer genomes especially in chromosomes (chr) 8 and 11. With the large number of personal genomes and cancer genomes becoming available, it is now possible to study their association to diseases using computational methods. Here, we conducted a pilot study on chromosomes 8 and 11 of cancer genomes to identify computationally the differentially occurring palindromes.

Methods

We processed 69 breast cancer genomes from The Cancer Genome Atlas including serum-normal and tumor genomes, and 1000 Genomes to serve as control group. The Biological Language Modelling Toolkit (BLMT) computes palindromes in whole genomes. We developed a computational pipeline integrating BLMT to compute and compare prevalence of palindromes in personal genomes.

Results

We carried out a pilot study on chr 8 and chr 11 taking into account single nucleotide polymorphisms, insertions and deletions. Of all the palindromes that showed any variation in cancer genomes, 38% of what were near breast cancer genes happened to be the most differentiated palindromes in tumor (i.e. they ranked among the top 25% by our heuristic measure).

Conclusions

These results will shed light on the prevalence of palindromes in oncogenes and the mutations that are present in the palindromic regions that could contribute to genomic rearrangements, and breast cancer progression.

Landscape of somatic mutations in 560 breast cancer whole-genome sequences

We analysed whole-genome sequences of 560 breast cancers to advance understanding of the driver mutations conferring clonal advantage and the mutational processes generating somatic mutations. We found that 93 protein-coding cancer genes carried probable driver mutations. Some non-coding regions exhibited high mutation frequencies, but most have distinctive structural features probably causing elevated mutation rates and do not contain driver mutations. Mutational signature analysis was extended to genome rearrangements and revealed twelve base substitution and six rearrangement signatures. Three rearrangement signatures, characterized by tandem duplications or deletions, appear associated with defective homologous-recombination-based DNA repair: one with deficient BRCA1 function, another with deficient BRCA1 or BRCA2 function, the cause of the third is unknown. This analysis of all classes of somatic mutation across exons, introns and intergenic regions highlights the repertoire of cancer genes and mutational processes operating, and progresses towards a comprehensive account of the somatic genetic basis of breast cancer.

A bouquet of DNA structures: Emerging diversity

Structural polymorphism of DNA has constantly been evolving from the time of illustration of the double helical model of DNA by Watson and Crick. A variety of non-canonical DNA structures have constantly been documented across the globe. DNA attracted worldwide attention as a carrier of genetic information. In addition to the classical Watson–Crick duplex, DNA can actually adopt diverse structures during its active participation in cellular processes like replication, transcription, recombination and repair. Structures like hairpin, cruciform, triplex, G-triplex, quadruplex, i-motif and other alternative non-canonical DNA structures have been studied at length and have also shown their in vivo occurrence. This review mainly focuses on non-canonical structures adopted by DNA oligonucleotides which have certain prerequisites for their formation in terms of sequence, its length, number and orientation of strands along with varied solution conditions. This conformational polymorphism of DNA might be the basis of different functional properties of a specific set of DNA sequences, further giving some insights for various extremely complicated biological phenomena. Many of these structures have already shown their linkages with diseases like cancer and genetic disorders, hence making them an extremely striking target for structure-specific drug designing and therapeutic applications.

•

DNA can adopt diverse range of structures other than classical Watson–Crick duplex.

•

Discussion of alternate structures like hairpin, cruciform, triplex, quadruplex etc.

•

This review gives some insights for the biological relevance of DNA structures.

Imaging of DNA and Protein-DNA Complexes with Atomic Force Microscopy

This article reviews atomic force microscopy (AFM) studies of DNA structure and dynamics and protein-DNA complexes, including recent advances in the visualization of protein-DNA complexes with the use of cutting-edge, high-speed AFM. Special emphasis is given to direct nanoscale visualization of dynamics of protein-DNA complexes. In the area of DNA structure and dynamics, structural studies of local non-B conformations of DNA and the interplay of local and global DNA conformations are reviewed. The application of time-lapse AFM nanoscale imaging of DNA dynamics is illustrated by studies of Holliday junction branch migration. Structure and dynamics of protein-DNA interactions include problems related to site-specific DNA recombination, DNA replication, and DNA mismatch repair. Studies involving the structure and dynamics of chromatin are also described.

The linear plastid chromosomes of maize: terminal sequences, structures, and implications for DNA replication

The structure of a chromosomal DNA molecule may influence the way in which it is replicated and inherited. For decades plastid DNA (ptDNA) was believed to be circular, with breakage invoked to explain linear forms found upon extraction from the cell. Recent evidence indicates that ptDNA in vivo consists of linear molecules with discrete termini, although these ends were not characterized. We report the sequences of two terminal regions, End1 and End2, for maize (Zea mays L.) ptDNA. We describe structural features of these terminal regions and similarities found in other plant ptDNAs. The terminal sequences are within inverted repeat regions (leading to four genomic isomers) and adjacent to origins of replication. Conceptually, stem-loop structures may be formed following melting of the double-stranded DNA ends. Exonuclease digestion indicates that the ends in maize are unobstructed, but tobacco (Nicotiana tabacum L.) ends may have a 5'-protein. If the terminal structure of ptDNA molecules influences the retention of ptDNA, the unprotected molecular ends in mature leaves of maize may be more susceptible to degradation in vivo than the protected ends in tobacco. The terminal sequences and cumulative GC skew profiles are nearly identical for maize, wheat (Triticum aestivum L.) and rice (Oryza sativa L.), with less similarity among other plants. The linear structure is now confirmed for maize ptDNA and inferred for other plants and suggests a virus-like recombination-dependent replication mechanism for ptDNA. Plastid transformation vectors containing the terminal sequences may increase the chances of success in generating transplastomic cereals.

In Silico/In Vivo Insights into the Functional and Evolutionary Pathway of Pseudomonas aeruginosa Oleate-Diol Synthase. Discovery of a New Bacterial Di-Heme Cytochrome C Peroxidase Subfamily

As previously reported, P. aeruginosa genes PA2077 and PA2078 code for 10S-DOX (10S-Dioxygenase) and 7,10-DS (7,10-Diol Synthase) enzymes involved in long-chain fatty acid oxygenation through the recently described oleate-diol synthase pathway. Analysis of the amino acid sequence of both enzymes revealed the presence of two heme-binding motifs (CXXCH) on each protein. Phylogenetic analysis showed the relation of both proteins to bacterial di-heme cytochrome c peroxidases (Ccps), similar to Xanthomonas sp. 35Y rubber oxidase RoxA. Structural homology modelling of PA2077 and PA2078 was achieved using RoxA (pdb 4b2n) as a template. From the 3D model obtained, presence of significant amino acid variations in the predicted heme-environment was found. Moreover, the presence of palindromic repeats located in enzyme-coding regions, acting as protein evolution elements, is reported here for the first time in P. aeruginosa genome. These observations and the constructed phylogenetic tree of the two proteins, allow the proposal of an evolutionary pathway for P. aeruginosa oleate-diol synthase operon. Taking together the in silico and in vivo results obtained we conclude that enzymes PA2077 and PA2078 are the first described members of a new subfamily of bacterial peroxidases, designated as Fatty acid-di-heme Cytochrome cperoxidases (FadCcp).

Short Inverted Repeats Are Hotspots for Genetic Instability: Relevance to Cancer Genomes

Analyses of chromosomal aberrations in human genetic disorders have revealed that inverted repeat sequences (IRs) often co-localize with endogenous chromosomal instability and breakage hotspots. Approximately 80% of all IRs in the human genome are short (<100 bp), yet the mutagenic potential of such short cruciform-forming sequences has not been characterized. Here, we find that short IRs are enriched at translocation breakpoints in human cancer and stimulate the formation of DNA double-strand breaks (DSBs) and deletions in mammalian and yeast cells. We provide evidence for replication-related mechanisms of IR-induced genetic instability and a novel XPF cleavage-based mechanism independent of DNA replication. These discoveries implicate short IRs as endogenous sources of DNA breakage involved in disease etiology and suggest that these repeats represent a feature of genome plasticity that may contribute to the evolution of the human genome by providing a means for diversity within the population.

detectIR: a novel program for detecting perfect and imperfect inverted repeats using complex numbers and vector calculation

Inverted repeats are present in abundance in both prokaryotic and eukaryotic genomes and can form DNA secondary structures--hairpins and cruciforms that are involved in many important biological processes. Bioinformatics tools for efficient and accurate detection of inverted repeats are desirable, because existing tools are often less accurate and time consuming, sometimes incapable of dealing with genome-scale input data. Here, we present a MATLAB-based program called detectIR for the perfect and imperfect inverted repeat detection that utilizes complex numbers and vector calculation and allows genome-scale data inputs. A novel algorithm is adopted in detectIR to convert the conventional sequence string comparison in inverted repeat detection into vector calculation of complex numbers, allowing non-complementary pairs (mismatches) in the pairing stem and a non-palindromic spacer (loop or gaps) in the middle of inverted repeats. Compared with existing popular tools, our program performs with significantly higher accuracy and efficiency. Using genome sequence data from HIV-1, Arabidopsis thaliana, Homo sapiens and Zea mays for comparison, detectIR can find lots of inverted repeats missed by existing tools whose outputs often contain many invalid cases. detectIR is open source and its source code is freely available at: https://sourceforge.net/projects/detectir.

Structural and biophysical characterization of murine rif1 C terminus reveals high specificity for DNA cruciform structures

Mammalian Rif1 is a key regulator of DNA replication timing, double-stranded DNA break repair, and replication fork restart. Dissecting the molecular functions of Rif1 is essential to understand how it regulates such diverse processes. However, Rif1 is a large protein that lacks well defined functional domains and is predicted to be largely intrinsically disordered; these features have hampered recombinant expression of Rif1 and subsequent functional characterization. Here we applied ESPRIT (expression of soluble proteins by random incremental truncation), an in vitro evolution-like approach, to identify high yielding soluble fragments encompassing conserved regions I and II (CRI and CRII) at the C-terminal region of murine Rif1. NMR analysis showed CRI to be intrinsically disordered, whereas CRII is partially folded. CRII binds cruciform DNA with high selectivity and micromolar affinity and thus represents a functional DNA binding domain. Mutational analysis revealed an α-helical region of CRII to be important for cruciform DNA binding and identified critical residues. Thus, we present the first structural study of the mammalian Rif1, identifying a domain that directly links its function to DNA binding. The high specificity of Rif1 for cruciform structures is significant given the role of this key protein in regulating origin firing and DNA repair.

Impact of bulge loop size on DNA triplet repeat domains: Implications for DNA repair and expansion

Repetitive DNA sequences exhibit complex structural and energy landscapes, populated by metastable, noncanonical states, that favor expansion and deletion events correlated with disease phenotypes. To probe the origins of such genotype-phenotype linkages, we report the impact of sequence and repeat number on properties of (CNG) repeat bulge loops. We find the stability of duplexes with a repeat bulge loop is controlled by two opposing effects; a loop junction-dependent destabilization of the underlying double helix, and a self-structure dependent stabilization of the repeat bulge loop. For small bulge loops, destabilization of the underlying double helix overwhelms any favorable contribution from loop self-structure. As bulge loop size increases, the stabilizing loop structure contribution dominates. The role of sequence on repeat loop stability can be understood in terms of its impact on the opposing influences of junction formation and loop structure. The nature of the bulge loop affects the thermodynamics of these two contributions differently, resulting in unique differences in repeat size-dependent minima in the overall enthalpy, entropy, and free energy changes. Our results define factors that control repeat bulge loop formation; knowledge required to understand how this helix imperfection is linked to DNA expansion, deletion, and disease phenotypes.

A MATLAB-based tool for accurate detection of perfect overlapping and nested inverted repeats in DNA sequences

SUMMARY

Palindromic sequences, or inverted repeats (IRs), in DNA sequences involve important biological processes such as DNA-protein binding, DNA replication and DNA transposition. Development of bioinformatics tools that are capable of accurately detecting perfect IRs can enable genome-wide studies of IR patterns in both prokaryotes and eukaryotes. Different from conventional string-comparison approaches, we propose a novel algorithm that uses a cumulative score system based on a prime number representation of nucleotide bases. We then implemented this algorithm as a MATLAB-based program for perfect IR detection. In comparison with other existing tools, our program demonstrates a high accuracy in detecting nested and overlapping IRs.

AVAILABILITY AND IMPLEMENTATION

The source code is freely available on (http://bioinfolab.miamioh.edu/bioinfolab/palindrome.php)

CONTACT

liangc@miamioh.edu or karroje@miamioh.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

The yin and yang of repair mechanisms in DNA structure-induced genetic instability

DNA can adopt a variety of secondary structures that deviate from the canonical Watson-Crick B-DNA form. More than 10 types of non-canonical or non-B DNA secondary structures have been characterized, and the sequences that have the capacity to adopt such structures are very abundant in the human genome. Non-B DNA structures have been implicated in many important biological processes and can serve as sources of genetic instability, implicating them in disease and evolution. Non-B DNA conformations interact with a wide variety of proteins involved in replication, transcription, DNA repair, and chromatin architectural regulation. In this review, we will focus on the interactions of DNA repair proteins with non-B DNA and their roles in genetic instability, as the proteins and DNA involved in such interactions may represent plausible targets for selective therapeutic intervention.

DNA secondary structures: stability and function of G-quadruplex structures

In addition to the canonical double helix, DNA can fold into various other inter- and intramolecular secondary structures. Although many such structures were long thought to be in vitro artefacts, bioinformatics demonstrates that DNA sequences capable of forming these structures are conserved throughout evolution, suggesting the existence of non-B-form DNA in vivo. In addition, genes whose products promote formation or resolution of these structures are found in diverse organisms, and a growing body of work suggests that the resolution of DNA secondary structures is critical for genome integrity. This Review focuses on emerging evidence relating to the characteristics of G-quadruplex structures and the possible influence of such structures on genomic stability and cellular processes, such as transcription.

Conformation effects of CpG methylation on single-stranded DNA oligonucleotides: analysis of the opioid peptide dynorphin-coding sequences

Single-stranded DNA (ssDNA) is characterized by high conformational flexibility that allows these molecules to adopt a variety of conformations. Here we used native polyacrylamide gel electrophoresis (PAGE), circular dichroism (CD) spectroscopy and nuclear magnetic resonance (NMR) spectroscopy to show that cytosine methylation at CpG sites affects the conformational flexibility of short ssDNA molecules. The CpG containing 37-nucleotide PDYN (prodynorphin) fragments were used as model molecules. The presence of secondary DNA structures was evident from differences in oligonucleotide mobilities on PAGE, from CD spectra, and from formation of A-T, G-C, and non-canonical G-T base pairs observed by NMR spectroscopy. The oligonucleotides displayed secondary structures at 4°C, and some also at 37°C. Methylation at CpG sites prompted sequence-dependent formation of novel conformations, or shifted the equilibrium between different existing ssDNA conformations. The effects of methylation on gel mobility and base pairing were comparable in strength to the effects induced by point mutations in the DNA sequences. The conformational effects of methylation may be relevant for epigenetic regulatory events in a chromatin context, including DNA-protein or DNA-DNA recognition in the course of gene transcription, and DNA replication and recombination when double-stranded DNA is unwinded to ssDNA.

Direct measurements of torque during Holliday junction migration

DNA experiences torsional stress resulting from the activities of motor enzymes and bound proteins. The mechanisms by which this torsional stress is dissipated to maintain DNA structural integrity are not fully known. Here, we show that a Holliday junction can limit torsion by coupling rotation to translocation and torque to force. The torque required to mechanically migrate through individual junctions was found to be an order of magnitude smaller than that required to melt DNA. We also directly show that substantially more torque was required to migrate through even a single-base sequence heterology, which has important implications for the activity of junction-migrating enzymes.

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.

Folded DNA in action: hairpin formation and biological functions in prokaryotes

Structured forms of DNA with intrastrand pairing are generated in several cellular processes and are involved in biological functions. These structures may arise on single-stranded DNA (ssDNA) produced during replication, bacterial conjugation, natural transformation, or viral infections. Furthermore, negatively supercoiled DNA can extrude inverted repeats as hairpins in structures called cruciforms. Whether they are on ssDNA or as cruciforms, hairpins can modify the access of proteins to DNA, and in some cases, they can be directly recognized by proteins. Folded DNAs have been found to play an important role in replication, transcription regulation, and recognition of the origins of transfer in conjugative elements. More recently, they were shown to be used as recombination sites. Many of these functions are found on mobile genetic elements likely to be single stranded, including viruses, plasmids, transposons, and integrons, thus giving some clues as to the manner in which they might have evolved. We review here, with special focus on prokaryotes, the functions in which DNA secondary structures play a role and the cellular processes giving rise to them. Finally, we attempt to shed light on the selective pressures leading to the acquisition of functions for DNA secondary structures.

Large inverted repeats within Xp11.2 are present at the breakpoints of isodicentric X chromosomes in Turner syndrome

Turner syndrome (TS) results from whole or partial monosomy X and is mediated by haploinsufficiency of genes that normally escape X-inactivation. Although a 45,X karyotype is observed in half of all TS cases, the most frequent variant TS karyotype includes the isodicentric X chromosome alone [46,X,idic(X)(p11)] or as a mosaic [46,X,idic(X)(p11)/45,X]. Given the mechanism of idic(X)(p11) rearrangement is poorly understood and breakpoint sequence information is unknown, this study sought to investigate the molecular mechanism of idic(X)(p11) formation by determining their precise breakpoint intervals. Karyotype analysis and fluorescence in situ hybridization mapping of eight idic(X)(p11) cell lines and three unbalanced Xp11.2 translocation lines identified the majority of breakpoints within a 5 Mb region, from approximately 53 to 58 Mb, in Xp11.1-p11.22, clustering into four regions. To further refine the breakpoints, a high-resolution oligonucleotide microarray (average of approximately 350 bp) was designed and array-based comparative genomic hybridization (aCGH) was performed on all 11 idic(X)(p11) and Xp11.2 translocation lines. aCGH analyses identified all breakpoint regions, including an idic(X)(p11) line with two potential breakpoints, one breakpoint shared between two idic(X)(p11) lines and two Xp translocations that shared breakpoints with idic(X)(p11) lines. Four of the breakpoint regions included large inverted repeats composed of repetitive gene clusters and segmental duplications, which corresponded to regions of copy-number variation. These data indicate that the rearrangement sites on Xp11.2 that lead to isodicentric chromosome formation and translocations are probably not random and suggest that the complex repetitive architecture of this region predisposes it to rearrangements, some of which are recurrent.

The distribution of inverted repeat sequences in the Saccharomyces cerevisiae genome

Although a variety of possible functions have been proposed for inverted repeat sequences (IRs), it is not known which of them might occur in vivo. We investigate this question by assessing the distributions and properties of IRs in the Saccharomyces cerevisiae (SC) genome. Using the IRFinder algorithm we detect 100,514 IRs having copy length greater than 6 bp and spacer length less than 77 bp. To assess statistical significance we also determine the IR distributions in two types of randomization of the S. cerevisiae genome. We find that the S. cerevisiae genome is significantly enriched in IRs relative to random. The S. cerevisiae IRs are significantly longer and contain fewer imperfections than those from the randomized genomes, suggesting that processes to lengthen and/or correct errors in IRs may be operative in vivo. The S. cerevisiae IRs are highly clustered in intergenic regions, while their occurrence in coding sequences is consistent with random. Clustering is stronger in the 3' flanks of genes than in their 5' flanks. However, the S. cerevisiae genome is not enriched in those IRs that would extrude cruciforms, suggesting that this is not a common event. Various explanations for these results are considered.

Cellular pathways controlling integron cassette site folding

By mobilizing small DNA units, integrons have a major function in the dissemination of antibiotic resistance among bacteria. The acquisition of gene cassettes occurs by recombination between the attI and attC sites catalysed by the IntI1 integron integrase. These recombination reactions use an unconventional mechanism involving a folded single-stranded attC site. We show that cellular bacterial processes delivering ssDNA, such as conjugation and replication, favour proper folding of the attC site. By developing a very sensitive in vivo assay, we also provide evidence that attC sites can recombine as cruciform structures by extrusion from double-stranded DNA. Moreover, we show an influence of DNA superhelicity on attC site extrusion in vitro and in vivo. We show that the proper folding of the attC site depends on both the propensity to form non-recombinogenic structures and the length of their variable terminal structures. These results draw the network of cell processes that regulate integron recombination.

Evidence for a large double-cruciform DNA structure on the X chromosome of human and chimpanzee

The human X chromosome consists of a high number of large inverted repeat (IR) DNA sequences which fulfill all requirements for formation of cruciform DNA structures. Such alternative DNA structures are suggested to have a great impact in altering the chromatin architecture and function. Our comprehensive analysis of the corresponding orthologous nucleotide sequences of an IR sequence from Homo sapiens and Pan troglodytes revealed that most of the nucleotide differences between the two species are symmetrical to the apex of the IR, and that the spacer region of the orthologous IRs are in reverse orientation. We provide evidence that this IR forms a large non-B DNA structure containing two Holliday junctions, allowing intrastrand nucleotide pairing of the arms and interstrand pairing of the spacer region of the IR. This structure would extrude into a large double-cruciform DNA structure providing the molecular basis of translocation events and regulation of gene expression.

Deletion of the cruciform binding domain in CBP/14-3-3 displays reduced origin binding and initiation of DNA replication in budding yeast

Background

Initiation of eukaryotic DNA replication involves many protein-protein and protein-DNA interactions. We have previously shown that 14-3-3 proteins bind cruciform DNA and associate with mammalian and yeast replication origins in a cell cycle dependent manner.

Results

By expressing the human 14-3-3ε, as the sole member of 14-3-3 proteins family in Saccharomyces cerevisiae, we show that 14-3-3ε complements the S. cerevisiae Bmh1/Bmh2 double knockout, conserves its cruciform binding activity, and associates in vivo with the yeast replication origins ARS307. Deletion of the α5-helix, the potential cruciform binding domain of 14-3-3, decreased the cruciform binding activity of the protein as well as its association with the yeast replication origins ARS307 and ARS1. Furthermore, the mutant cells had a reduced ability to stably maintain plasmids bearing one or multiple origins.

Conclusion

14-3-3, a cruciform DNA binding protein, associates with yeast origins of replication and functions as an initiator of DNA replication, presumably through binding to cruciform DNA forming at yeast replicators.

Unusually long palindromes are abundant in mitochondrial control regions of insects and nematodes

Background

Palindromes are known to be involved in a variety of biological processes. In the present investigation we carried out a comprehensive analysis of palindromes in the mitochondrial control regions (CRs) of several animal groups to study their frequency, distribution and architecture to gain insights into the origin of replication of mtDNA.

Methodology/Principal Findings

Many species of Arthropoda, Nematoda, Mollusca and Annelida harbor palindromes and inverted repeats (IRs) in their CRs. Lower animals like cnidarians and higher animal groups like chordates are almost devoid of palindromes and IRs. The study revealed that palindrome occurrence is positively correlated with the AT content of CRs, and that IRs are likely to give rise to longer palindromes.

Conclusions/Significance

The present study attempts to explain possible reasons and gives in silico evidence for absence of palindromes and IRs from CR of vertebrate mtDNA and acquisition and retention of the same in insects. Study of CRs of different animal phyla uncovered unique architecture of this locus, be it high abundance of long palindromes and IRs in CRs of Insecta and Nematoda, or short IRs of 10–20 nucleotides with a spacer region of 12–14 bases in subphylum Chelicerata, or nearly complete of absence of any long palindromes and IRs in Vertebrata, Cnidaria and Echinodermata.

Influence of cationic molecules on the hairpin to duplex equilibria of self-complementary DNA and RNA oligonucleotides

A self-complementary nucleotide sequence can form both a unimolecular hairpin and a bimolecular duplex. In this study, the secondary structures of the self-complementary DNA and RNA oligonucleotides with different sequences and lengths were investigated under various solution conditions by gel electrophoresis, circular dichroism (CD) and electron paramagnetic resonance (EPR) spectroscopy and a ultraviolet (UV) melting analysis. The DNA sequences tended to adopt a hairpin conformation at low cation concentrations, but a bimolecular duplex was preferentially formed at an elevated cationic strength. On the other hand, fully matched RNA sequences adopted a bimolecular duplex regardless of the cation concentration. The thermal melting experiments indicated a greater change in the melting temperature of the bimolecular duplexes (by approximately 20 degrees C) than that of the hairpin (by approximately 10 degrees C) by increasing the NaCl concentration from 10 mM to 1 M. Hairpin formations were also observed for the palindrome DNA sequences derived from Escherichia coli, but association of the complementary palindrome sequences was observed when spermine, one of the major cationic molecules in a cell, existed at the physiological concentration. The results indicate the role of cations for shifting the structural equilibrium toward a nucleotide assembly and implicate nucleotide structures in cells.

Zygote arrest 1 gene in pig, cattle and human: evidence of different transcript variants in male and female germ cells

BACKGROUND

Zygote arrest 1 (ZAR1) is one of the few known oocyte-specific maternal-effect genes essential for the beginning of embryo development discovered in mice. This gene is evolutionary conserved in vertebrates and ZAR1 protein is characterized by the presence of atypical plant homeobox zing finger domain, suggesting its role in transcription regulation. This work was aimed at the study of this gene, which could be one of the key regulators of successful preimplantation development of domestic animals, in pig and cattle, as compared with human.

METHODS

Screenings of somatic cell hybrid panels and in silico research were performed to characterize ZAR1 chromosome localization and sequences. Rapid amplification of cDNA ends was used to obtain full-length cDNAs. Spatio-temporal mRNA expression patterns were studied using Northern blot, reverse transcription coupled to polymerase chain reaction and in situ hybridization.

RESULTS

We demonstrated that ZAR1 is a single copy gene, positioned on chromosome 8 in pig and 6 in cattle, and several variants of correspondent cDNA were cloned from oocytes. Sequence analysis of ZAR1 cDNAs evidenced numerous short inverted repeats within the coding sequences and putative Pumilio-binding and embryo-deadenylation elements within the 3'-untranslated regions, indicating the potential regulation ways. We showed that ZAR1 expressed exclusively in oocytes in pig ovary, persisted during first cleavages in embryos developed in vivo and declined sharply in morulae and blastocysts. ZAR1 mRNA was also detected in testis, and, at lower level, in hypothalamus and pituitary in both species. For the first time, ZAR1 was localized in testicular germ cells, notably in round spermatids. In addition, in pig, cattle and human only shorter ZAR1 transcript variants resulting from alternative splicing were found in testis as compared to oocyte.

CONCLUSION

Our data suggest that in addition to its role in early embryo development highlighted by expression pattern of full-length transcript in oocytes and early embryos, ZAR1 could also be implicated in the regulation of meiosis and post meiotic differentiation of male and female germ cells through expression of shorter splicing variants. Species conservation of ZAR1 expression and regulation underlines the central role of this gene in early reproductive processes.

Inverted repeat structure of the human genome: the X-chromosome contains a preponderance of large, highly homologous inverted repeats that contain testes genes

We have performed the first genome-wide analysis of the Inverted Repeat (IR) structure in the human genome, using a novel and efficient software package called Inverted Repeats Finder (IRF). After masking of known repetitive elements, IRF detected 22,624 human IRs characterized by arm size from 25 bp to >100 kb with at least 75% identity, and spacer length up to 100 kb. This analysis required 6 h on a desktop PC. In all, 166 IRs had arm lengths >8 kb. From this set, IRs were excluded if they were in unfinished/unassembled regions of the genome, or clustered with other closely related IRs, yielding a set of 96 large IRs. Of these, 24 (25%) occurred on the X-chromosome, although it represents only approximately 5% of the genome. Of the X-chromosome IRs, 83.3% were >/=99% identical, compared with 28.8% of autosomal IRs. Eleven IRs from Chromosome X, one from Chromosome 11, and seven already described from Chromosome Y contain genes predominantly expressed in testis. PCR analysis of eight of these IRs correctly amplified the corresponding region in the human genome, and six were also confirmed in gorilla or chimpanzee genomes. Similarity dot-plots revealed that 22 IRs contained further secondary homologous structures partially categorized into three distinct patterns. The prevalence of large highly homologous IRs containing testes genes on the X- and Y-chromosomes suggests a possible role in male germ-line gene expression and/or maintaining sequence integrity by gene conversion.

The human beta-globin replication initiation region consists of two modular independent replicators

Previous studies have shown that mammalian cells contain replicator sequences, which can determine where DNA replication initiates. However, the specific sequences that confer replicator activity were not identified. Here we report a detailed analysis of replicator sequences that dictate initiation of DNA replication from the human beta-globin locus. This analysis suggests that the beta-globin replication initiation region contains two adjacent, redundant replicators. Each replicator was capable of initiating DNA replication independently at ectopic sites. Within each of these two replicators, we identified short, discrete, nonredundant sequences, which cooperatively determine replicator activity. Experiments with somatic cell hybrids further demonstrated that the requirements for initiation at ectopic sites were similar to the requirements for initiation within native human chromosomes. The replicator clustering and redundancy exemplified in the human beta-globin locus may account for the extreme difficulty in identifying replicator sequences in mammalian cells and suggest that mammalian replication initiation sites may be determined by cooperative sequence modules.

Biopolymer Chain Elasticity: A novel concept and a least deformation energy principle predicts backbone and overall folding of DNA TTT hairpins in agreement with NMR distances

A new molecular modelling methodology is presented and shown to apply to all published solution structures of DNA hairpins with TTT in the loop. It is based on the theory of elasticity of thin rods and on the assumption that single-stranded B-DNA behaves as a continuous, unshearable, unstretchable and flexible thin rod. It requires four construction steps: (i) computation of the tri-dimensional trajectory of the elastic line, (ii) global deformation of single-stranded helical DNA onto the elastic line, (iii) optimisation of the nucleoside rotations about the elastic line, (iv) energy minimisation to restore backbone bond lengths and bond angles. This theoretical approach called 'Biopolymer Chain Elasticity' (BCE) is capable of reproducing the tri-dimensional course of the sugar-phosphate chain and, using NMR-derived distances, of reproducing models close to published solution structures. This is shown by computing three different types of distance criteria. The natural description provided by the elastic line and by the new parameter, Omega, which corresponds to the rotation angles of nucleosides about the elastic line, offers a considerable simplification of molecular modelling of hairpin loops. They can be varied independently from each other, since the global shape of the hairpin loop is preserved in all cases.