High Nucleotide Skew Palindromic DNA Sequences Function as Potential Replication Origins due to their Unzipping Propensity

Locations of DNA replication initiation in prokaryotes, called "origins of replication", are well-characterized. However, a mechanistic understanding of the sequence dependence of the local unzipping of double-stranded DNA, the first step towards replication initiation, is lacking. Here, utilizing a Markov chain model that was created to address the directional nature of DNA unzipping and replication, we model the sequence dependence of local melting of double-stranded linear DNA segments. We show that generalized palindromic sequences with high nucleotide skews have a low kinetic barrier for local melting near melting temperatures. This allows for such sequences to function as potential replication origins. We support our claim with evidence for high-skew palindromic sequences within the replication origins of mitochondrial DNA, bacteria, archaea and plasmids.

Unexplored power of CRISPR-Cas9 in neuroscience, a multi-OMICs review

The neuroscience and neurobiology of gene editing to enhance learning and memory is of paramount interest to the scientific community. The advancements of CRISPR system have created avenues to treat neurological disorders by means of versatile modalities varying from expression to suppression of genes and proteins. Neurodegenerative disorders have also been attributed to non-canonical DNA secondary structures by affecting neuron activity through controlling gene expression, nucleosome shape, transcription, translation, replication, and recombination. Changing DNA regulatory elements which could contribute to the fate and function of neurons are thoroughly discussed in this review. This study presents the ability of CRISPR system to boost learning power and memory, treat or cure genetically-based neurological disorders, and alleviate psychiatric diseases by altering the activity and the irritability of the neurons at the synaptic cleft through DNA manipulation, and also, epigenetic modifications using Cas9. We explore and examine how each different OMIC techniques can come useful when altering DNA sequences. Such insight into the underlying relationship between OMICs and cellular behaviors leads us to better neurological and psychiatric therapeutics by intelligently designing and utilizing the CRISPR/Cas9 technology.

Non-canonical DNA structures in the human ribosomal DNA

Non-canonical structures (NCS) refer to the various forms of DNA that differ from the B-conformation described by Watson and Crick. It has been found that these structures are usual components of the genome, actively participating in its essential functions. The present review is focused on the nine kinds of NCS appearing or likely to appear in human ribosomal DNA (rDNA): supercoiling structures, R-loops, G-quadruplexes, i-motifs, DNA triplexes, cruciform structures, DNA bubbles, and A and Z DNA conformations. We discuss the conditions of their generation, including their sequence specificity, distribution within the locus, dynamics, and beneficial and detrimental role in the cell.

Cellular Responses to Widespread DNA Replication Stress

Replicative DNA polymerases are blocked by nearly all types of DNA damage. The resulting DNA replication stress threatens genome stability. DNA replication stress is also caused by depletion of nucleotide pools, DNA polymerase inhibitors, and DNA sequences or structures that are difficult to replicate. Replication stress triggers complex cellular responses that include cell cycle arrest, replication fork collapse to one-ended DNA double-strand breaks, induction of DNA repair, and programmed cell death after excessive damage. Replication stress caused by specific structures (e.g., G-rich sequences that form G-quadruplexes) is localized but occurs during the S phase of every cell division. This review focuses on cellular responses to widespread stress such as that caused by random DNA damage, DNA polymerase inhibition/nucleotide pool depletion, and R-loops. Another form of global replication stress is seen in cancer cells and is termed oncogenic stress, reflecting dysregulated replication origin firing and/or replication fork progression. Replication stress responses are often dysregulated in cancer cells, and this too contributes to ongoing genome instability that can drive cancer progression. Nucleases play critical roles in replication stress responses, including MUS81, EEPD1, Metnase, CtIP, MRE11, EXO1, DNA2-BLM, SLX1-SLX4, XPF-ERCC1-SLX4, Artemis, XPG, FEN1, and TATDN2. Several of these nucleases cleave branched DNA structures at stressed replication forks to promote repair and restart of these forks. We recently defined roles for EEPD1 in restarting stressed replication forks after oxidative DNA damage, and for TATDN2 in mitigating replication stress caused by R-loop accumulation in BRCA1-defective cells. We also discuss how insights into biological responses to genome-wide replication stress can inform novel cancer treatment strategies that exploit synthetic lethal relationships among replication stress response factors.

Antagonistic roles of canonical and Alternative-RPA in disease-associated tandem CAG repeat instability

Expansions of repeat DNA tracts cause >70 diseases, and ongoing expansions in brains exacerbate disease. During expansion mutations, single-stranded DNAs (ssDNAs) form slipped-DNAs. We find the ssDNA-binding complexes canonical replication protein A (RPA1, RPA2, and RPA3) and Alternative-RPA (RPA1, RPA3, and primate-specific RPA4) are upregulated in Huntington disease and spinocerebellar ataxia type 1 (SCA1) patient brains. Protein interactomes of RPA and Alt-RPA reveal unique and shared partners, including modifiers of CAG instability and disease presentation. RPA enhances in vitro melting, FAN1 excision, and repair of slipped-CAGs and protects against CAG expansions in human cells. RPA overexpression in SCA1 mouse brains ablates expansions, coincident with decreased ATXN1 aggregation, reduced brain DNA damage, improved neuron morphology, and rescued motor phenotypes. In contrast, Alt-RPA inhibits melting, FAN1 excision, and repair of slipped-CAGs and promotes CAG expansions. These findings suggest a functional interplay between the two RPAs where Alt-RPA may antagonistically offset RPA's suppression of disease-associated repeat expansions, which may extend to other DNA processes.

Characterization and sequencing analysis of pLP2.5-11 and pLP3.0-4 novel cryptic plasmids from Lactiplantibacillus plantarum WP72/27

We sequenced and described two cryptic plasmids from Lactiplantibacillus plantarum strain WP72/27, termed pLP2.5-11 (OP831909) and pLP3.0-4 (OP831910). Nucleotide sequencing gave the sizes of pLP2.5-11 and pLP3.0-4 as 2754 and 3197 base pairs, with G + C contents 38.89% and 40.88% and predicted two and eight putative open reading frames, respectively. The RepA protein of pLP2.5-11 shared a 99% identity with pC30il, pLP1 and pC30il, whereas the RepB protein of pLP3.0-4 shared a 98% identity with pXY3, a member of the rolling-circle replication (RCR) pC194 family. The origin of plasmid replication was predicted to consist of inverted and directed repeat sequences upstream of the Rep genes. Sequence analysis predicted that both pLP2.5-11 and pLP3.0-4 plasmids replicate via a rolling-circle process.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03684-y.

Attenuation of a DNA cruciform by a conserved regulator directs T3SS1 mediated virulence in Vibrio parahaemolyticus

Pathogenic Vibrio species account for 3-5 million annual life-threatening human infections. Virulence is driven by bacterial hemolysin and toxin gene expression often positively regulated by the winged helix-turn-helix (wHTH) HlyU transcriptional regulator family and silenced by histone-like nucleoid structural protein (H-NS). In the case of Vibrio parahaemolyticus, HlyU is required for virulence gene expression associated with type 3 Secretion System-1 (T3SS1) although its mechanism of action is not understood. Here, we provide evidence for DNA cruciform attenuation mediated by HlyU binding to support concomitant virulence gene expression. Genetic and biochemical experiments revealed that upon HlyU mediated DNA cruciform attenuation, an intergenic cryptic promoter became accessible allowing for exsA mRNA expression and initiation of an ExsA autoactivation feedback loop at a separate ExsA-dependent promoter. Using a heterologous E. coli expression system, we reconstituted the dual promoter elements which revealed that HlyU binding and DNA cruciform attenuation were strictly required to initiate the ExsA autoactivation loop. The data indicate that HlyU acts to attenuate a transcriptional repressive DNA cruciform to support T3SS1 virulence gene expression and reveals a non-canonical extricating gene regulation mechanism in pathogenic Vibrio species.

Multiple rearrangements and low inter- and intra-species mitogenome sequence variation in the Heterobasidion annosum s.l. species complex

Introduction

Mitochondria are essential organelles in the eukaryotic cells and responsible for the energy production but are also involved in many other functions including virulence of some fungal species. Although the evolution of fungal mitogenomes have been studied at some taxonomic levels there are still many things to be learned from studies of closely related species.

Methods

In this study, we have analyzed 60 mitogenomes in the five species of the Heterobasidion annosum sensu lato complex that all are necrotrophic pathogens on conifers.

Results and Discussion

Compared to other fungal genera the genomic and genetic variation between and within species in the complex was low except for multiple rearrangements. Several translocations of large blocks with core genes have occurred between the five species and rearrangements were frequent in intergenic areas. Mitogenome lengths ranged between 108 878 to 116 176 bp, mostly as a result of intron variation. There was a high degree of homology of introns, homing endonuclease genes, and intergenic ORFs among the five Heterobasidion species. Three intergenic ORFs with unknown function (uORF6, uORF8 and uORF9) were found in all five species and was located in conserved synteny blocks. A 13 bp long GC-containing self-complementary palindrome was discovered in many places in the five species that were optional in presence/absence. The within species variation is very low, among 48 H. parviporum mitogenomes, there was only one single intron exchange, and SNP frequency was 0.28% and indel frequency 0.043%. The overall low variation in the Heterobasidion annosum sensu lato complex suggests a slow evolution of the mitogenome.

Dynamic alternative DNA structures in biology and disease

Repetitive elements in the human genome, once considered 'junk DNA', are now known to adopt more than a dozen alternative (that is, non-B) DNA structures, such as self-annealed hairpins, left-handed Z-DNA, three-stranded triplexes (H-DNA) or four-stranded guanine quadruplex structures (G4 DNA). These dynamic conformations can act as functional genomic elements involved in DNA replication and transcription, chromatin organization and genome stability. In addition, recent studies have revealed a role for these alternative structures in triggering error-generating DNA repair processes, thereby actively enabling genome plasticity. As a driving force for genetic variation, non-B DNA structures thus contribute to both disease aetiology and evolution.

Genomic Analysis of Non-B Nucleic Acids Structures in SARS-CoV-2: Potential Key Roles for These Structures in Mutability, Translation, and Replication?

Non-B nucleic acids structures have arisen as key contributors to genetic variation in SARS-CoV-2. Herein, we investigated the presence of defining spike protein mutations falling within inverted repeats (IRs) for 18 SARS-CoV-2 variants, discussed the potential roles of G-quadruplexes (G4s) in SARS-CoV-2 biology, and identified potential pseudoknots within the SARS-CoV-2 genome. Surprisingly, there was a large variation in the number of defining spike protein mutations arising within IRs between variants and these were more likely to occur in the stem region of the predicted hairpin stem-loop secondary structure. Notably, mutations implicated in ACE2 binding and propagation (e.g., ΔH69/V70, N501Y, and D614G) were likely to occur within IRs, whilst mutations involved in antibody neutralization and reduced vaccine efficacy (e.g., T19R, ΔE156, ΔF157, R158G, and G446S) were rarely found within IRs. We also predicted that RNA pseudoknots could predominantly be found within, or next to, 29 mutations found in the SARS-CoV-2 spike protein. Finally, the Omicron variants BA.2, BA.4, BA.5, BA.2.12.1, and BA.2.75 appear to have lost two of the predicted G4-forming sequences found in other variants. These were found in nsp2 and the sequence complementary to the conserved stem-loop II-like motif (S2M) in the 3' untranslated region (UTR). Taken together, non-B nucleic acids structures likely play an integral role in SARS-CoV-2 evolution and genetic diversity.

Nucleic acid folding simulations using a physics-based atomistic free energy model

Performing full-resolution atomistic simulations of nucleic acid folding has remained a challenge for biomolecular modeling. Understanding how nucleic acids fold and how they transition between different folded structures as they unfold and refold has important implications for biology. This paper reports a theoretical model and computer simulation of the ab initio folding of DNA inverted repeat sequences. The formulation is based on an all-atom conformational model of the sugar-phosphate backbone via chain closure, and it incorporates three major molecular-level driving forces-base stacking, counterion-induced backbone self-interactions, and base pairing-via separate analytical theories designed to capture and reproduce the effects of the solvent without requiring explicit water and ions in the simulation. To accelerate computational throughput, a mixed numerical/analytical algorithm for the calculation of the backbone conformational volume is incorporated into the Monte Carlo simulation, and special stochastic sampling techniques were employed to achieve the computational efficiency needed to fold nucleic acids from scratch. This paper describes implementation details, benchmark results, and the advantages and technical challenges with this approach.

Cruciform DNA Structures Act as Legible Templates for Accelerating Homologous Recombination in Transgenic Animals

Inverted repeat (IR) DNA sequences compose cruciform structures. Some genetic disorders are the result of genome inversion or translocation by cruciform DNA structures. The present study examined whether exogenous DNA integration into the chromosomes of transgenic animals was related to cruciform DNA structures. Large imperfect cruciform structures were frequently predicted around predestinated transgene integration sites in host genomes of microinjection-based transgenic (Tg) animals (αLA-LPH Tg goat, Akr1A1eGFP/eGFP Tg mouse, and NFκB-Luc Tg mouse) or CRISPR/Cas9 gene-editing (GE) animals (αLA-AP1 GE mouse). Transgene cassettes were imperfectly matched with their predestinated sequences. According to the analyzed data, we proposed a putative model in which the flexible cruciform DNA structures acted as a legible template for DNA integration into linear DNAs or double-strand break (DSB) alleles. To demonstrate this model, artificial inverted repeat knock-in (KI) reporter plasmids were created to analyze the KI rate using the CRISPR/Cas9 system in NIH3T3 cells. Notably, the KI rate of the 5′ homologous arm inverted repeat donor plasmid (5′IR) with the ROSA gRNA group (31.5%) was significantly higher than the knock-in reporter donor plasmid (KIR) with the ROSA gRNA group (21.3%, p < 0.05). However, the KI rate of the 3′ inverted terminal repeat/inverted repeat donor plasmid (3′ITRIR) group was not different from the KIR group (23.0% vs. 22.0%). These results demonstrated that the legibility of the sequence with the cruciform DNA existing in the transgene promoted homologous recombination (HR) with a higher KI rate. Our findings suggest that flexible cruciform DNAs folded by IR sequences improve the legibility and accelerate DNA 3′-overhang integration into the host genome via homologous recombination machinery.

Structural switching/polymorphism by sequential base substitution at quasi-palindromic SNP site (G → A) in LCR of human β-globin gene cluster

The human β-globin gene Locus Control Region (LCR), a dominant regulator of globin gene expression contains five tissue-specific DNase I-hypersensitive sites (HSs). A single nucleotide polymorphism (SNP) (A → G) present in HS4 region of locus control region (LCR), have shown a notable association between the G allele and the occurrence of β-thalassemia. This SNP site exhibiting a hairpin - duplex equilibrium manifested in A → B like DNA transition has previously been reported from this laboratory. Since, DNA is a dynamic and adaptable molecule, so any change of a single base within a primary DNA sequence can produce major biological consequences commonly manifested in genetic disorders such as sickle cell anemia and β-thalassemia. Herein, the differential behavior of sequential single base substitutions G → A on the quasi-palindromic sequence (d-TGGGGGCCCCA; HPG11) has been explored. A combination of native gel electrophoresis, circular dichroism (CD), and UV-thermal denaturation (Tm) techniques have been used to investigate the structural polymorphism associated with various variants of HPG11 i.e. HPG11A2 to HPG11A5. The CD spectra confirmed that all the HPG11 variants exhibit a hairpin - duplex equilibrium. Oligomer concentration dependence on CD spectra has been correlated with A → B DNA conformational transition. However, as revealed in gel electrophoresis, HPG11A2 → A5 exhibit the formation of a tetramolecular structure (four-way junction) at higher oligomer concentration. UV-melting studies also supported the melting of hairpin, duplex and four-way junction structure. This polymorphism pattern may possibly be significant for DNA-protein recognition, in the process of regulation of LCR in the β-globin gene.

Inverted repeats in coronavirus SARS-CoV-2 genome manifest the evolution events

The world faces a great unforeseen challenge through the COVID-19 pandemic caused by coronavirus SARS-CoV-2. The virus genome structure and evolution are positioned front and center for further understanding insights on vaccine development, monitoring of transmission trajectories, and prevention of zoonotic infections of new coronaviruses. Of particular interest are genomic elements Inverse Repeats (IRs), which maintain genome stability, regulate gene expressions, and are the targets of mutations. However, little research attention is given to the IR content analysis in the SARS-CoV-2 genome. In this study, we propose a geometric analysis method and using the method to investigate the distributions of IRs in SARS-CoV-2 and its related coronavirus genomes. The method represents each genomic IR sequence pair as a single point and constructs the geometric shape of the genome using the IRs. Thus, the IR shape can be considered as the signature of the genome. The genomes of different coronaviruses are then compared using the constructed IR shapes. The results demonstrate that SARS-CoV-2 genome, specifically, has an abundance of IRs, and the IRs in coronavirus genomes show an increase during evolution events.

The minicircular and extremely heteroplasmic mitogenome of the holoparasitic plant Rhopalocnemis phalloides

The plastid and nuclear genomes of parasitic plants exhibit deeply altered architectures,^1-13 whereas the few examined mitogenomes range from deeply altered to conventional.^14-20 To provide further insight on mitogenome evolution in parasitic plants, we report the highly modified mitogenome of Rhopalocnemis phalloides, a holoparasite in Balanophoraceae. Its mitogenome is uniquely arranged in 21 minicircular chromosomes that vary in size from 4,949 to 7,861 bp, with a total length of only 130,713 bp. All chromosomes share an identical 896 bp conserved region, with a large stem-loop that acts as the origin of replication, flanked on each side by hypervariable and semi-conserved regions. Similar minicircular structures with shared and unique regions have been observed in parasitic animals and free-living protists,^21-24 suggesting convergent structural evolution. Southern blots confirm both the minicircular structure and the replication origin of the mitochondrial chromosomes. PacBio reads provide evidence for chromosome recombination and rolling-circle replication for the R. phalloides mitogenome. Despite its small size, the mitogenome harbors a typical set of genes and introns within the unique regions of each chromosome, yet introns are the smallest among seed plants and ferns. The mitogenome also exhibits extreme heteroplasmy, predominantly involving short indels and more complex variants, many of which cause potential loss-of-function mutations for some gene copies. All heteroplasmic variants are transcribed, and functional and nonfunctional protein-coding variants are spliced and RNA edited. Our findings offer a unique perspective into how mitogenomes of parasitic plants can be deeply altered and shed light on plant mitogenome 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.

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

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

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.

SARS-CoV-2 hot-spot mutations are significantly enriched within inverted repeats and CpG island loci

SARS-CoV-2 is an intensively investigated virus from the order Nidovirales (Coronaviridae family) that causes COVID-19 disease in humans. Through enormous scientific effort, thousands of viral strains have been sequenced to date, thereby creating a strong background for deep bioinformatics studies of the SARS-CoV-2 genome. In this study, we inspected high-frequency mutations of SARS-CoV-2 and carried out systematic analyses of their overlay with inverted repeat (IR) loci and CpG islands. The main conclusion of our study is that SARS-CoV-2 hot-spot mutations are significantly enriched within both IRs and CpG island loci. This points to their role in genomic instability and may predict further mutational drive of the SARS-CoV-2 genome. Moreover, CpG islands are strongly enriched upstream from viral ORFs and thus could play important roles in transcription and the viral life cycle. We hypothesize that hypermethylation of these loci will decrease the transcription of viral ORFs and could therefore limit the progression of the disease.

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.

Alternative DNA Structures In Vivo: Molecular Evidence and Remaining Questions

Duplex DNA naturally folds into a right-handed double helix in physiological conditions. Some sequences of unusual base composition may nevertheless form alternative structures, as was shown for many repeated sequences in vitro However, evidence for the formation of noncanonical structures in living cells is difficult to gather. It mainly relies on genetic assays demonstrating their function in vivo or through genetic instability reflecting particular properties of such structures. Efforts were made to reveal their existence directly in a living cell, mainly by generating antibodies specific to secondary structures or using chemical ligands selected for their affinity to these structures. Among secondary structure-forming DNAs are G-quadruplexes, human fragile sites containing minisatellites, AT-rich regions, inverted repeats able to form cruciform structures, hairpin-forming CAG/CTG triplet repeats, and triple helices formed by homopurine-homopyrimidine GAA/TTC trinucleotide repeats. Many of these alternative structures are involved in human pathologies, such as neurological or developmental disorders, as in the case of trinucleotide repeats, or cancers triggered by translocations linked to fragile sites. This review will discuss and highlight evidence supporting the formation of alternative DNA structures in vivo and will emphasize the role of the mismatch repair machinery in binding mispaired DNA duplexes, triggering genetic instability.

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.

Interaction of a photosensitizer methylene blue with various structural forms (cruciform, bulge duplex and hairpin) of designed DNA sequences

Functionally important, local structural transitions in DNA generate various alternative conformations. Cruciform is one of such alternative DNA structures, usually targeted in genomes by various proteins. Symmetry elements in sequence as inverted repeats are the key factor for cruciform formation, facilitated by the presence of the AT-rich regions. Here, we used biophysical and biochemical techniques such as Gel electrophoresis, Circular dichroism (CD), and UV-thermal melting analysis to explore the structural status of the designed DNA sequences, which had potential to form cruciform structures under physiological conditions. The gel electrophoresis analysis revealed that the designed 53-mer DNA oligonucleotide sequence CR forms an intermolecular bulge duplex with flanking ends, while another sequence CRC adopts an intramolecular hairpin structure with flanking ends. Their equimolar complex (CRCRC) bestowed much-retarded migration due to the formation of a quite intriguing cruciform structure. CD studies confirmed that all the alternative structures (cruciform, bulge duplex, and hairpin with flanking ends) exhibit characteristics of B-DNA type conformation. A triphasic UV-thermal melting curve displayed by the complex formed by the equimolar ratio (CRCRC) is also suggestive of the formation of the cruciform structure. The interaction studies of CR, CRC, and their equimolar complex (1:1) with a photosensitizer methylene blue (MB) indicated that MB could not stabilize the discrete structures formed by CR and CRC sequences, however, the cruciform structure showed a quite significant increment in the melting temperature. Such studies facilitate our understanding of various secondary structures possibly present inside the cell and their interactions with drug/dye molecules.

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.

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.

Novel diversity in mitochondrial genomes of deep-sea Pennatulacea (Cnidaria: Anthozoa: Octocorallia)

We present the first documented complete mitogenomes of deep-sea Pennatulacea, representing nine genera and eight families. These include one species each of the deep-sea genera Funiculina, Halipteris, Protoptilum and Distichoptilum, four species each of Umbellula and Pennatula, three species of Kophobelemnon and two species of Anthoptilum, as well as one species of the epi- and mesobenthic genus Virgularia. Seventeen circular genomes ranged from 18,513 bp (Halipteris cf. finmarchica) to 19,171 bp (Distichoptilum gracile) and contained all genes standard to octocoral mitochondrial genomes (14 protein-coding genes, two ribosomal RNA genes and one transfer RNA). We found at least three different gene orders in Pennatulacea: the ancestral gene order, the gene order found in bamboo corals (Family Isididae), and a novel gene order. The mitogenome of one species of Umbellula has a bipartite genome (∼13 kbp and ∼5 kbp), with good evidence that both parts are circular.

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.

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.

Effect of GCAA stabilizing loops on three- and four-way intramolecular junctions

Tetraloops are a common way of changing the melting behavior of a DNA or RNA structure without changing the sequence of the stem. Because of the ubiquitous nature of tetraloops, our goal is to understand the effect a GCAA tetraloop, which belongs to the GNRA family of tetraloops, has on the unfolding thermodynamics of intramolecular junctions. Specifically, we have described the melting behavior of intramolecular three-way and four-way junctions where a T₅ loop has been replaced with a GCAA tetraloops in different positions. Their thermodynamic profiles, including Δn_Na⁺ and Δn_W, were analyzed based on the position of the tetraloop. We obtained between -16.7 and -27.5 kcal mol^-1 for all junctions studied. The experimental data indicates the influence of the GCAA tetraloop is primarily dictated by the native unfolding of the junction; if the tetraloop is placed on a stem that unfolds as a single domain when the tetraloop is not present, it will unfold as a single domain when the tetraloop is present but with a higher thermal stability. Conversely, if the tetraloop is placed on a stem which unfolds cooperatively with other stems when the tetraloop is not present, the tetraloop will increase the thermal stability of all the stems in the melting domain. The oligonucleotide structure and not the tetraloop itself affects ion uptake; three-way junctions do not gain an increase in ion uptake, but four-way junctions do. This is not the case for water immobilization, where the position of the tetraloop dictates the amount of water immobilized.

Isogenic mutations in the Moraxella catarrhalis CydDC system display pleiotropic phenotypes and reveal the role of a palindrome sequence in its transcriptional regulation

Moraxella catarrhalis is becoming an important human respiratory tract pathogen affecting significant proportions from the population. However, still little is known about its physiology and molecular regulation. To this end, the CydDC, which is a heterodimeric ATP binding cassette transporter that has been shown to contribute to the maintenance of the redox homeostasis across the periplasm in other Gram-negative bacteria, is studied here. Amino acids multiple sequence alignments indicated that M. catarrhalis CydC is different from the CydC proteins of the bacterial species in which this system has been previously studied. These findings prompted further interest in studying this system in M. catarrhalis. Isogenic mutant in the CydDC system showed suppression in growth rate, hypersensitivity to oxidative and reductive stress and increased accumulation of intracellular cysteine levels. In addition, the growth of cydC^- mutant exhibited hypersensitivity to exogenous cysteine; however, it did not display a significant difference from its wild-type counterpart in the murine pulmonary clearance model. Moreover, a palindrome was detected 94bp upstream of the cydD ORF suggesting it might act as a potential regulatory element. Real-time reverse transcription-PCR analysis showed that deletion/change in the palindrome resulted into alterations in the transcription levels of cydC. A better understanding of such system and its regulation helps in developing better ways to combat M. catarrhalis infections.

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.

Palindrome analyser - A new web-based server for predicting and evaluating inverted repeats in nucleotide sequences

DNA cruciform structures play an important role in the regulation of natural processes including gene replication and expression, as well as nucleosome structure and recombination. They have also been implicated in the evolution and development of diseases such as cancer and neurodegenerative disorders. Cruciform structures are formed by inverted repeats, and their stability is enhanced by DNA supercoiling and protein binding. They have received broad attention because of their important roles in biology. Computational approaches to study inverted repeats have allowed detailed analysis of genomes. However, currently there are no easily accessible and user-friendly tools that can analyse inverted repeats, especially among long nucleotide sequences. We have developed a web-based server, Palindrome analyser, which is a user-friendly application for analysing inverted repeats in various DNA (or RNA) sequences including genome sequences and oligonucleotides. It allows users to search and retrieve desired gene/nucleotide sequence entries from the NCBI databases, and provides data on length, sequence, locations and energy required for cruciform formation. Palindrome analyser also features an interactive graphical data representation of the distribution of the inverted repeats, with options for sorting according to the length of inverted repeat, length of loop, and number of mismatches. Palindrome analyser can be accessed at http://bioinformatics.ibp.cz.

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.

The Integron: Adaptation On Demand

The integron is a powerful system which, by capturing, stockpiling, and rearranging new functions carried by gene encoding cassettes, confers upon bacteria a rapid adaptation capability in changing environments. Chromosomally located integrons (CI) have been identified in a large number of environmental Gram-negative bacteria. Integron evolutionary history suggests that these sedentary CIs acquired mobility among bacterial species through their association with transposable elements and conjugative plasmids. As a result of massive antibiotic use, these so-called mobile integrons are now widespread in clinically relevant bacteria and are considered to be the principal agent in the emergence and rise of antibiotic multiresistance in Gram-negative bacteria. Cassette rearrangements are catalyzed by the integron integrase, a site-specific tyrosine recombinase. Central to these reactions is the single-stranded DNA nature of one of the recombination partners, the attC site. This makes the integron a unique recombination system. This review describes the current knowledge on this atypical recombination mechanism, its implications in the reactions involving the different types of sites, attC and attI, and focuses on the tight regulation exerted by the host on integron activity through the control of attC site folding. Furthermore, cassette and integrase expression are also highly controlled by host regulatory networks and the bacterial stress (SOS) response. These intimate connections to the host make the integron a genetically stable and efficient system, granting the bacteria a low cost, highly adaptive evolution potential "on demand".

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.

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.