The human genome-wide distribution of DNA palindromes

Computational analysis on the dissemination of non-B DNA structural motifs in promoter regions of 1180 cellular genomes

The promoter regions of gene regulation are under evolutionary constraints and earlier studies uncovered that they are characterized by enrichment of functional non-B DNA structural signatures like curved DNA, cruciform DNA, G-quadruplex, triple-helical DNA, slipped DNA structures, and Z-DNA. However, these studies are restricted to a few model organisms, single non-B DNA motif types, or whole genomic sequences, and their comparative accumulation in promoter regions of different domains of life has not been reported comprehensively. In this study, for the first time, we investigated the preponderance of non-B DNA-prone motifs in promoter regions in 1180 genomes belonging to 28 taxonomic groups using the non-B DNA Motif Search Tool (nBMST). The trends suggest that they are predominant in promoters compared to the upstream and downstream regions of all three domains of life and variably linked to taxonomic groups. Cruciform DNA motif is the most abundant form of non-B DNA, spanning from archaea to lower eukaryotes. Curved DNA motifs are prominent in host-associated bacteria, and suppressed in mammals. Triplex-DNA and slipped DNA structure repeats are discretely dispersed in all lineages. G-quadruplex motifs are significantly enriched in mammals. We also observed that the unique enrichment of non-B DNA in promoters is strongly linked to genome GC, size, evolutionary time divergence, and ecological adaptations. Overall, our work systematically reports the unique non-B DNA structural landscape of cellular organisms from the perspective of the cis-regulatory code of genomes.

Palindrome-Embedded Hairpin Structure and Its Target-Catalyzed Padlock Cyclization for Label-Free MicroRNA-Initiated Rolling Circle Amplification

Highly sensitive detection of microRNAs (miRNAs) is of great significance in early diagnosis of cancers. Here, we develop a palindrome-embedded hairpin structure and its target-catalyzed padlock cyclization for rolling circle amplification, named PHP-RCA for simplicity, which can be applied in label-free ultrasensitive detection of miRNA. PHP-RCA is a facile system that consists of only an oligonucleotide probe with a palindrome-embedded hairpin structure (PHP). The two ends of PHP were extended as overhangs and designed with the complementary sequences of the target. Hence, the phosphorylated PHP can be cyclized by T4 DNA ligase in the presence of the target that serves as the ligation template. This ligation has formed a palindrome-embedded dumbbell-shaped probe (PDP) that allows phi29 polymerase to perform a typical target-primed RCA on PDP by taking miRNA as a primer, resulting in the production of a lengthy tandem repeat. Benefits from the palindromic sequences and hairpin-shaped structure in padlock double-stranded structures can be infinitely produced during the RCA reaction and provide numerous binding sites for SYBR Green I, a double-stranded dye, achieving a sharp response signal for label-free target detection. We have demonstrated that the proposed system exhibits a good linear range from 0.1 fM to 5 nM with a low detection limit of 0.1 fM, and the non-target miRNA can be clearly distinguished. The advantages of high efficiency, label-free signaling, and the use of only one oligonucleotide component make the PHP-RCA suitable for ultrasensitive, economic, and convenient detection of target miRNAs. This simple and powerful system is expected to provide a promising platform for tumor diagnosis, prognosis, and therapy.

An updated overview of experimental and computational approaches to identify non-canonical DNA/RNA structures with emphasis on G-quadruplexes and R-loops

Multiple types of non-canonical nucleic acid structures play essential roles in DNA recombination and replication, transcription, and genomic instability and have been associated with several human diseases. Thus, an increasing number of experimental and bioinformatics methods have been developed to identify these structures. To date, most reviews have focused on the features of non-canonical DNA/RNA structure formation, experimental approaches to mapping these structures, and the association of these structures with diseases. In addition, two reviews of computational algorithms for the prediction of non-canonical nucleic acid structures have been published. One of these reviews focused only on computational approaches for G4 detection until 2020. The other mainly summarized the computational tools for predicting cruciform, H-DNA and Z-DNA, in which the algorithms discussed were published before 2012. Since then, several experimental and computational methods have been developed. However, a systematic review including the conformation, sequencing mapping methods and computational prediction strategies for these structures has not yet been published. The purpose of this review is to provide an updated overview of conformation, current sequencing technologies and computational identification methods for non-canonical nucleic acid structures, as well as their strengths and weaknesses. We expect that this review will aid in understanding how these structures are characterised and how they contribute to related biological processes and diseases.

Tumor-targeting [2]catenane-based grid-patterned periodic DNA monolayer array for in vivo theranostic application

DNA nanotechnology is often used to build various nano-structures for signaling and/or drug delivery, but it essentially suffers from several major limitations, such as a large number of DNA strands and limited targeting ligands. Moreover, there is no report on in vivo two-dimensional DNA arrays because of various technical challenges. By cross-catenating two palindromic DNA rings, herein, we demonstrate a catenane-based grid-patterned periodic DNA monolayer array ([2]GDA) capable of preferentially accumulating in tumor tissues without any targeting ligands, with a thickness equal to the double-helical DNA monolayer (nearly 2 nm). The structural flexibility of [2]GDA enabled it to fold into a spherical object in solution, favoring cellular uptake. Thus, its cellular internalization activity was comparable with that of the commercial lipofectamine 3000. Moreover, [2]GDA retained the structural integrity over 24 h incubation in biological solutions, achieving a 360-fold improvement in in vivo stability. Significantly, anticancer drug-loaded [2]GDA exhibits desirable therapeutic efficacy in tumor-bearing animals without detectable side effects.

Palindromes in DNA-A Risk for Genome Stability and Implications in Cancer

A palindrome in DNA consists of two closely spaced or adjacent inverted repeats. Certain palindromes have important biological functions as parts of various cis-acting elements and protein binding sites. However, many palindromes are known as fragile sites in the genome, sites prone to chromosome breakage which can lead to various genetic rearrangements or even cell death. The ability of certain palindromes to initiate genetic recombination lies in their ability to form secondary structures in DNA which can cause replication stalling and double-strand breaks. Given their recombinogenic nature, it is not surprising that palindromes in the human genome are involved in genetic rearrangements in cancer cells as well as other known recurrent translocations and deletions associated with certain syndromes in humans. Here, we bring an overview of current understanding and knowledge on molecular mechanisms of palindrome recombinogenicity and discuss possible implications of DNA palindromes in carcinogenesis. Furthermore, we overview the data on known palindromic sequences in the human genome and efforts to estimate their number and distribution, as well as underlying mechanisms of genetic rearrangements specific palindromic sequences cause.

Conventional PCR assisted single-component assembly of spherical nucleic acids for simple colorimetric detection of SARS-CoV-2

Continuous identification of suspected infectious cases is crucial to control the recent pandemic caused by the novel human coronavirus SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2). Real-time polymerase chain reaction (real-time PCR) technology cannot be implemented easily and in large scale in some communities due to lack of resources and infrastructures. Here, we report a simple colorimetric strategy derived from linker-based single-component assembly of gold nanoparticle-core spherical nucleic acids (AuNP-core SNAs) for visual detection of PCR products of SARS-CoV-2 ribonucleic acid (RNA) template. A palindromic linker is designed based on SARS-CoV-2 specific E gene to program the identical colloidal SNAs into large assemblies along with a distinct red-to-purple color change. The linker acts as a probe of SARS-CoV-2 RNA in conventional PCR reaction. In the presence of the correct template the palindromic linker, which is complementary to a short region within the target amplicon, is cleaved by 5'-exonuclease activity of deoxyribonucleic acid (DNA) polymerase. Cleavage of the palindromic linker during the amplification process inhibits the single-component assembly formation of SNAs. So, positive and negative viral samples produce simply red and purple colors in the post PCR colorimetric test, respectively. Evaluation of the samples obtained from cases with laboratory-confirmed SARS-CoV-2 infection revealed that our assay can rival with real-time PCR method in sensitivity.

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.

Regulation of gene expression by the action of a fungal lncRNA on a transactivator

Long non-coding RNAs (lncRNAs) are crucial factors acting on regulatory processes in eukaryotes. Recently, for the first time in a filamentous fungus, the lncRNA HAX1 was characterized in the ascomycete Trichoderma reesei. In industry, this fungus is widely applied for the high-yield production of cellulases. The lncRNA HAX1 was reported to influence the expression of cellulase-encoding genes; interestingly, this effect is dependent on the presence of its most abundant length. Clearly, HAX1 acts in association with a set of well-described transcription factors to regulate gene expression. In this study, we attempted to elucidate the regulatory strategy of HAX1 and its interactions with the major transcriptional activator Xylanase regulator 1 (Xyr1). We demonstrated that HAX1 interferes with the negative feedback regulatory loop of Xyr1 in a sophisticated manner and thus ultimately has a positive effect on gene expression.

Palindromic Molecular Beacon Based Z-Scheme BiOCl-Au-CdS Photoelectrochemical Biodetection

This work presented an innovative and rationally engineered palindromic molecular beacon (PMB) based "Z-scheme" photoelectrochemical (PEC) biosensing protocol for the selective screening of kanamycin (Kana) through DNA hybridization-induced conformational conversion. Interestingly, the ingeniously designed PMB integrated the multifunctional elements including recognition region, primer-like palindromic fragment, and polymerization-nicking template. The cosensitized structures consisted of CdS quantum dot functionalized hairpin DNA2 (QD-HP2) and region-selectively deposited gold nanoparticles onto {001} facets of bismuth oxychloride (BiOCl-Au). Compared with BiOCl-Au alone, the attachment of CdS QDs onto BiOCl-Au (i.e., BiOCl-Au-CdS QDs) exhibited evidently enhanced photocurrent intensity thanks to the synergistic effect of Z-scheme BiOCl-Au-CdS QDs. After incubation with target Kana, Kana-aptamer binding could induce the exposure of PMB region for hairpin DNA1 (HP1). The exposed palindromic tails hybridized with each other (like a molecular machine) to consume the substrates (dNTPs) and fuels (enzyme) for the releasing of numerous nick fragments (Nick). The as-generated nick fragments could specifically hybridize with the complementary region of QD-HP2, thus resulting in decreasing photocurrent because of the increasing spatial distance for electron transfer between two-type photosensitizers. Under optimum conditions, the PMB-based sensing system exhibited satisfying photocurrent responses toward target Kana within the working range from 50 to 5000 fM at a low detection limit of 29 fM. Impressively, the concept of a palindromic fragment-mediated primer-free biosensing strategy offers a new avenue for advanced development of efficient and convenient biodetection systems.

Symmetric exponential amplification reaction-based DNA nanomachine for the fluorescent detection of nucleic acids

By introducing palindromic sequences into the classical exponential amplification reaction (EXPAR), we constructed a new palindromic fragment-incorporated multifunctional hairpin probe (P-HP)-mediated symmetric exponential amplification reaction (S-EXPAR), to significantly reduce the background signal caused by inherent nonspecific amplification. A G-triplex/ThT complex was used as the signal reporter for the proposed label-free DNA nanomachine. The P-HP consists of five functional regions: a C-rich region (C), a target DNA recognition region (T′), two nicking sites (X′) and a palindromic fragment (P). When target DNA (T) hybridizes with P-HP, the palindromic fragment at the 3′ end of P-HP is fully exposed. Then, the P-HP/T duplexes hybridize with each other through the exposed P, and EXPAR occurs automatically and continuously on both sides of P under the synergistic effect of polymerase and nicking endonuclease. This is called the S-EXPAR assay. In this system, one T converts to a large number of G-triplex fragments, which can combine with ThT within a short time. The G-triplex/ThT complexes formed act as the signal reporter in a label-free and environmentally friendly format. In this way, the limit of detection of this method is as low as 10 pM with a dynamic response range of 10 pM to 300 nM. In addition, this method can detect other nucleic acids by simply changing the T′ region of the P-HP. Thus, the proposed DNA nanomachine is a potential alternative method for nucleic acid detection.

This label-free and ultra-low background signal DNA nanomachine was based on P-HP mediated S-EXPAR and the G-triplex/ThT complex.

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.

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.

Promoter analysis reveals globally differential regulation of human long non-coding RNA and protein-coding genes

Transcriptional regulation of protein-coding genes is increasingly well-understood on a global scale, yet no comparable information exists for long non-coding RNA (lncRNA) genes, which were recently recognized to be as numerous as protein-coding genes in mammalian genomes. We performed a genome-wide comparative analysis of the promoters of human lncRNA and protein-coding genes, finding global differences in specific genetic and epigenetic features relevant to transcriptional regulation. These two groups of genes are hence subject to separate transcriptional regulatory programs, including distinct transcription factor (TF) proteins that significantly favor lncRNA, rather than coding-gene, promoters. We report a specific signature of promoter-proximal transcriptional regulation of lncRNA genes, including several distinct transcription factor binding sites (TFBS). Experimental DNase I hypersensitive site profiles are consistent with active configurations of these lncRNA TFBS sets in diverse human cell types. TFBS ChIP-seq datasets confirm the binding events that we predicted using computational approaches for a subset of factors. For several TFs known to be directly regulated by lncRNAs, we find that their putative TFBSs are enriched at lncRNA promoters, suggesting that the TFs and the lncRNAs may participate in a bidirectional feedback loop regulatory network. Accordingly, cells may be able to modulate lncRNA expression levels independently of mRNA levels via distinct regulatory pathways. Our results also raise the possibility that, given the historical reliance on protein-coding gene catalogs to define the chromatin states of active promoters, a revision of these chromatin signature profiles to incorporate expressed lncRNA genes is warranted in the future.

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.

Competitive superhelical transitions involving cruciform extrusion

A DNA molecule under negative superhelical stress becomes susceptible to transitions to alternate structures. The accessible alternate conformations depend on base sequence and compete for occupancy. We have developed a method to calculate equilibrium distributions among the states available to such systems, as well as their average thermodynamic properties. Here we extend this approach to include superhelical cruciform extrusion at both perfect and imperfect inverted repeat (IR) sequences. We find that short IRs do not extrude cruciforms, even in the absence of competition. But as the length of an IR increases, its extrusion can come to dominate both strand separation and B-Z transitions. Although many IRs are present in human genomic DNA, we find that extrusion-susceptible ones occur infrequently. Moreover, their avoidance of transcription start sites in eukaryotes suggests that cruciform formation is rarely involved in mechanisms of gene regulation. We examine a set of clinically important chromosomal translocation breakpoints that occur at long IRs, whose rearrangement has been proposed to be driven by cruciform extrusion. Our results show that the susceptibilities of these IRs to cruciform formation correspond closely with their observed translocation frequencies.

A method to find palindromes in nucleic acid sequences

Various types of sequences in the human genome are known to play important roles in different aspects of genomic functioning. Among these sequences, palindromic nucleic acid sequences are one such type that have been studied in detail and found to influence a wide variety of genomic characteristics. For a nucleotide sequence to be considered as a palindrome, its complementary strand must read the same in the opposite direction. For example, both the strands i.e the strand going from 5' to 3' and its complementary strand from 3' to 5' must be complementary. A typical nucleotide palindromic sequence would be TATA (5' to 3') and its complimentary sequence from 3' to 5' would be ATAT. Thus, a new method has been developed using dynamic programming to fetch the palindromic nucleic acid sequences. The new method uses less memory and thereby it increases the overall speed and efficiency. The proposed method has been tested using the bacterial (3891 KB bases) and human chromosomal sequences (Chr-18: 74366 kb and Chr-Y: 25554 kb) and the computation time for finding the palindromic sequences is in milli seconds.

Analysis of compensatory substitution and gene evolution on the MAGEA/CSAG-palindrome of the primate X chromosomes

The human X chromosome contains a large number of inverted repeat DNA palindromes. Although arbitrary substitutions destroyed the inverted repeat structure of MAGEA/CSAG-palindrome during the evolutionary process of the primates, most of the substitutions are compensatory. Using maximum parsimony, it is demonstrated that the compensatory substitutions are prone to occur between bases with similar structures on the human, chimpanzee and orangutan MAGEA/CSAG-palindromes. Furthermore, it is found that MAGEA/CSAG genes also exist in orangutan and rhesus monkey palindromes by homologous searching. This suggests that the MAGEA/CSAG-palindrome might predate the divergence of human and other primate lineages. Comparative sequence analysis of the arms and genes on the primate MAGEA/CSAG-palindromes provides possible evidence of subsequently arm to arm gene conversion. These compensatory substitutions on the MAGEA/CSAG-palindrome of the primate X chromosomes play an important role in maintaining their structural symmetry during the process of formation.

Compositional bias is a major determinant of the distribution pattern and abundance of palindromes in Drosophila melanogaster

Palindromic sequences are important DNA motifs related to gene regulation, DNA replication and recombination, and thus, investigating the evolutionary forces shaping the distribution pattern and abundance of palindromes in the genome is substantially important. In this article, we analyzed the abundance of palindromes in the genome, and then explored the possible effects of several genomic factors on the palindrome distribution and abundance in Drosophila melanogaster. Our results show that the palindrome abundance in D. melanogaster deviates from random expectation and the uneven distribution of palindromes across the genome is associated with local GC content, recombination rate, and coding exon density. Our data suggest that base composition is the major determinant of the distribution pattern and abundance of palindromes and the correlation between palindrome density and recombination is a side-product of the effect of compositional bias on the palindrome abundance.

Chromosomal translocations and palindromic AT-rich repeats

Repetitive DNA sequences constitute 30% of the human genome, and are often sites of genomic rearrangement. Recently, it has been found that several constitutional translocations, especially those that involve chromosome 22, take place utilizing palindromic sequences on 22q11 and on the partner chromosome. Analysis of translocation junction fragments shows that the breakpoints of such palindrome-mediated translocations are localized at the center of palindromic AT-rich repeats (PATRRs). The presence of PATRRs at the breakpoints indicates a palindrome-mediated mechanism involved in the generation of these constitutional translocations. Identification of these PATRR-mediated translocations suggests a universal pathway for gross chromosomal rearrangement in the human genome. De novo occurrences of PATRR-mediated translocations can be detected by PCR in normal sperm samples but not somatic cells. Polymorphisms of various PATRRs influence their propensity for adopting a secondary structure, which in turn affects de novo translocation frequency. We propose that the PATRRs form an unstable secondary structure, which leads to double-strand breaks at the center of the PATRR. The double-strand breaks appear to be followed by a non-homologous end-joining repair pathway, ultimately leading to the translocations. This review considers recent findings concerning the mechanism of meiosis-specific, PATRR-mediated translocations.

Mechanisms of unphosphorylated STAT3 transcription factor binding to DNA

Phosphorylation of signal transducer and activator of transcription 3 (STAT3) on a single tyrosine residue in response to growth factors, cytokines, interferons, and oncogenes activates its dimerization, translocation to the nucleus, binding to the interferon γ (gamma)-activated sequence (GAS) DNA-binding site and activation of transcription of target genes. STAT3 is constitutively phosphorylated in various cancers and drives gene expression from GAS-containing promoters to promote tumorigenesis. Recently, roles for unphosphorylated STAT3 (U-STAT3) have been described in response to cytokine stimulation, in cancers, and in maintenance of heterochromatin stability. However, the mechanisms underlying U-STAT3 binding to DNA has not been fully investigated. Here, we explore STAT3-DNA interactions by atomic force microscopy (AFM) imaging. We observed that U-STAT3 molecules bind to the GAS DNA-binding site as dimers and monomers. In addition, we observed that U-STAT3 binds to AT-rich DNA sequence sites and recognizes specific DNA structures, such as 4-way junctions and DNA nodes, within negatively supercoiled plasmid DNA. These structures are important for chromatin organization and our data suggest a role for U-STAT3 as a chromatin/genome organizer. Unexpectedly, we found that a C-terminal truncated 67.5-kDa STAT3 isoform recognizes single-stranded spacers within cruciform structures that also have a role in chromatin organization and gene expression. This isoform appears to be abundant in the nuclei of cancer cells and, therefore, may have a role in regulation of gene expression. Taken together, our data highlight novel mechanisms by which U-STAT3 binds to DNA and supports U-STAT3 function as a transcriptional activator and a chromatin/genomic organizer.

Comparative genome analysis suggests characteristics of yeast inverted repeats that are important for transcriptional activity

Inverted repeats are sequences of DNA that, when read in the 5' to 3' direction, have the same sequence on both strands (palindromic portion), with the exception of a small number of nucleotides in the exact center (nonpalindromic spacer). They have been implicated in various DNA-mediated processes including replication, transcription, and genomic instability. At least some of these sequences are capable of forming an alternative DNA structure, called a cruciform, that may be important for mediating these functions. We generated a list of inverted repeats in the Saccharomyces cerevisiae genome and determined which of them are conserved in three related yeasts. We have identified characterisitics of inverted repeats that make them more likely to be conserved than the surrounding DNA and characteristics, such as position and base composition, that make the genes they are associated with likely to be more actively transcribed. This is an important step in determining the functions of this group of genomic elements.

Species-specific typing of DNA based on palindrome frequency patterns

DNA in its natural, double-stranded form may contain palindromes, sequences which read the same from either side because they are identical to their reverse complement on the sister strand. Short palindromes are underrepresented in all kinds of genomes. The frequency distribution of short palindromes exhibits more than twice the inter-species variance of non-palindromic sequences, which renders palindromes optimally suited for the typing of DNA. Here, we show that based on palindrome frequency, DNA sequences can be discriminated to the level of species of origin. By plotting the ratios of actual occurrence to expectancy, we generate palindrome frequency patterns that allow to cluster different sequences of the same genome and to assign plasmids, and in some cases even viruses to their respective host genomes. This finding will be of use in the growing field of metagenomics.

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.

Inverted duplications on acentric markers: mechanism of formation

Acentric inverted duplication (inv dup) markers, the largest group of chromosomal abnormalities with neocentromere formation, are found in patients both with idiopathic mental retardation and with cancer. The mechanism of their formation has been investigated by analyzing the breakpoints and the genotypes of 12 inv dup marker cases (three trisomic, six tetrasomic, two polysomic and one X chromosome derived marker) using a combination of fluorescence in situ hybridization, quantitative SNP array and microsatellite analysis. Inv dup markers were found to form either symmetrically with one breakpoint or asymmetrically with two distinct breakpoints. Genotype analyses revealed that all inv dup markers formed from one single chromatid end. This observation is incompatible with the previously suggested model by which the acentric inv dup markers form through inter-chromosomal U-type exchange. On the basis of the identification of DNA sequence motifs with inverted homologies within all observed breakpoint regions, a new general mechanism is proposed for the acentric inv dup marker formation: following a double-strand break an acentric fragment forms, during either meiosis or mitosis. The open DNA end of the acentric fragment is stabilized by the formation of an intra-chromosomal loop promoted by the presence of sequences with inverted homologies. Likely coinciding with the neocentromere formation, this stabilized fragment is duplicated during an early mitotic event, insuring the marker's survival during cell division and its presence in all cells.

A study on genomic distribution and sequence features of human long inverted repeats reveals species-specific intronic inverted repeats

The inverted repeats present in a genome play dual roles. They can induce genomic instability and, on the other hand, regulate gene expression. In the present study, we report the distribution and sequence features of recombinogenic long inverted repeats (LIRs) that are capable of forming stable stem-loops or palindromes within the human genome. A total of 2551 LIRs were identified, and 37% of them were located in long introns (largely > 10 kb) of genes. Their distribution appears to be random in introns and is not restrictive, even for regions near intron-exon boundaries. Almost half of them comprise TG/CA-rich repeats, inversely arranged Alu repeats and MADE1 mariners. The remaining LIRs are mostly unique in their sequence features. Comparative studies of human, chimpanzee, rhesus monkey and mouse orthologous genes reveal that human genes have more recombinogenic LIRs than other orthologs, and over 80% are human-specific. The human genes associated with the human-specific LIRs are involved in the pathways of cell communication, development and the nervous system, as based on significantly over-represented Gene Ontology terms. The functional pathways related to the development and functions of the nervous system are not enriched in chimpanzee and mouse orthologs. The findings of the present study provide insight into the role of intronic LIRs in gene regulation and primate speciation.

Size-dependent palindrome-induced intrachromosomal recombination in yeast

Palindromic and quasi-palindromic sequences are important DNA motifs found in various cis-acting genetic elements, but are also known to provoke different types of genetic alterations. The instability of such motifs is clearly size-related and depends on their potential to adopt secondary structures known as hairpins and cruciforms. Here we studied the influence of palindrome size on recombination between two directly repeated copies of the yeast CYC1 gene leading to the loss of the intervening sequence ("pop-out" recombination). We show that palindromes inserted either within one copy or between the two copies of the CYC1 gene become recombinogenic only when they attain a certain critical size and we estimate this critical size to be about 70 bp. With the longest palindrome used in this study (150 bp) we observed a more than 20-fold increase in the pop-out recombination. In the sae2/com1 mutant the palindrome-stimulated recombination was completely abolished. Suppression of palindrome recombinogenicity may be crucial for the maintenance of genetic stability in organisms containing a significant number of large palindromes in their genomes, like humans.

Mus81-dependent double-strand DNA breaks at in vivo-generated cruciform structures in S. cerevisiae

Long DNA palindromes are implicated in chromosomal rearrangement, but their roles in the underlying molecular events remain a matter of conjecture. One notion is that palindromes induce DNA breaks after assuming a cruciform structure, the four-way DNA junction providing a target for cleavage by Holliday junction (HJ)-specific enzymes. Though compelling, few components of the "cruciform resolution" proposal are established. Here we address fundamental properties and genetic dependencies of palindromic DNA metabolism in eukaryotes. Plasmid-borne palindromes introduced into S. cerevisiae are site-specifically broken in vivo, and the breaks exhibit unique hallmarks of an HJ resolvase mechanism. In vivo resolution requires Mus81, for which the bacterial HJ resolvase RusA will substitute. These results provide confirmation of cruciform extrusion and resolution in the context of eukaryotic chromatin. Related observations are that, unchecked by a nuclease function provided by Mre11, episomal palindromes launch a self-perpetuating breakage-fusion-bridge-independent copy number increase termed "escape."

Chromosome evolution with naked eye: palindromic context of the life origin

Based on the representation of the DNA sequence as a two-dimensional (2D) plane walk, we consider the problem of identification and comparison of functional and structural organizations of chromosomes of different organisms. According to the characteristic design of 2D walks we identify telomere sites, palindromes of various sizes and complexity, areas of ribosomal RNA, transposons, as well as diverse satellite sequences. As an interesting result of the application of the 2D walk method, a new duplicated gigantic palindrome in the X human chromosome is detected. A schematic mechanism leading to the formation of such a duplicated palindrome is proposed. Analysis of a large number of the different genomes shows that some chromosomes (or their fragments) of various species appear as imperfect gigantic palindromes, which are disintegrated by many inversions and the mutation drift on different scales. A spread occurrence of these types of sequences in the numerous chromosomes allows us to develop a new insight of some accepted points of the genome evolution in the prebiotic phase.