Discriminative motif discovery in DNA and protein sequences using the DEME algorithm

<p>Abstract</p> <p>Background</p> <p>Motif discovery aims to detect short, highly conserved patterns in a collection of unaligned DNA or protein sequences. Discriminative motif finding algorithms aim to increase the sensitivity and selectivity of motif discovery by utilizing a second set of sequences, and searching only for patterns that can differentiate the two sets of sequences. Potential applications of discriminative motif discovery include discovering transcription factor binding site motifs in ChIP-chip data and finding protein motifs involved in thermal stability using sets of orthologous proteins from thermophilic and mesophilic organisms.</p> <p>Results</p> <p>We describe DEME, a discriminative motif discovery algorithm for use with protein and DNA sequences. Input to DEME is two sets of sequences; a "positive" set and a "negative" set. DEME represents motifs using a probabilistic model, and uses a novel combination of global and local search to find the motif that optimally discriminates between the two sets of sequences. DEME is unique among discriminative motif finders in that it uses an informative Bayesian prior on protein motif columns, allowing it to incorporate prior knowledge of residue characteristics. We also introduce four, synthetic, discriminative motif discovery problems that are designed for evaluating discriminative motif finders in various biologically motivated contexts. We test DEME using these synthetic problems and on two biological problems: finding yeast transcription factor binding motifs in ChIP-chip data, and finding motifs that discriminate between groups of thermophilic and mesophilic orthologous proteins.</p> <p>Conclusion</p> <p>Using artificial data, we show that DEME is more effective than a non-discriminative approach when there are "decoy" motifs or when a variant of the motif is present in the "negative" sequences. With real data, we show that DEME is as good, but not better than non-discriminative algorithms at discovering yeast transcription factor binding motifs. We also show that DEME can find highly informative thermal-stability protein motifs. Binaries for the stand-alone program DEME is free for academic use and is available at <url>http://bioinformatics.org.au/deme/</url></p

Genomic Structure of and Genome-Wide Recombination in the Saccharomyces cerevisiae S288C Progenitor Isolate EM93

The diploid isolate EM93 is the main ancestor to the widely used Saccharomyces cerevisiae haploid laboratory strain, S288C. In this study, we generate a high-resolution overview of the genetic differences between EM93 and S288C. We show that EM93 is heterozygous for >45,000 polymorphisms, including large sequence polymorphisms, such as deletions and a Saccharomyces paradoxus introgression. We also find that many large sequence polymorphisms (LSPs) are associated with Ty-elements and sub-telomeric regions. We identified 2,965 genetic markers, which we then used to genotype 120 EM93 tetrads. In addition to deducing the structures of all EM93 chromosomes, we estimate that the average EM93 meiosis produces 144 detectable recombination events, consisting of 87 crossover and 31 non-crossover gene conversion events. Of the 50 polymorphisms showing the highest levels of non-crossover gene conversions, only three deviated from parity, all of which were near heterozygous LSPs. We find that non-telomeric heterozygous LSPs significantly reduce meiotic recombination in adjacent intervals, while sub-telomeric LSPs have no discernable effect on recombination. We identified 203 recombination hotspots, relatively few of which are hot for both non-crossover gene conversions and crossovers. Strikingly, we find that recombination hotspots show limited conservation. Some novel hotspots are found adjacent to heterozygous LSPs that eliminate other hotspots, suggesting that hotspots may appear and disappear relatively rapidly

Structural Insights into the Quinolone Resistance Mechanism of Mycobacterium tuberculosis DNA Gyrase

Mycobacterium tuberculosis DNA gyrase, an indispensable nanomachine involved in the regulation of DNA topology, is the only type II topoisomerase present in this organism and is hence the sole target for quinolone action, a crucial drug active against multidrug-resistant tuberculosis. To understand at an atomic level the quinolone resistance mechanism, which emerges in extensively drug resistant tuberculosis, we performed combined functional, biophysical and structural studies of the two individual domains constituting the catalytic DNA gyrase reaction core, namely the Toprim and the breakage-reunion domains. This allowed us to produce a model of the catalytic reaction core in complex with DNA and a quinolone molecule, identifying original mechanistic properties of quinolone binding and clarifying the relationships between amino acid mutations and resistance phenotype of M. tuberculosis DNA gyrase. These results are compatible with our previous studies on quinolone resistance. Interestingly, the structure of the entire breakage-reunion domain revealed a new interaction, in which the Quinolone-Binding Pocket (QBP) is blocked by the N-terminal helix of a symmetry-related molecule. This interaction provides useful starting points for designing peptide based inhibitors that target DNA gyrase to prevent its binding to DNA

Where Do Introns Come From?

The systems for mRNA surveillance, capping, and cleavage/polyadenylation are proposed to play pivotal roles in the physical establishment and distribution of spliceosomal introns along a transcript

Role of Matrix Metalloproteinases and Therapeutic Benefits of Their Inhibition in Spinal Cord Injury

This review will focus on matrix metalloproteinases (MMPs) and their inhibitors in the context of spinal cord injury (SCI). MMPs have a specific cellular and temporal pattern of expression in the injured spinal cord. Here we consider their diverse functions in the acutely injured cord and during wound healing. Excessive activity of MMPs, and in particular gelatinase B (MMP-9), in the acutely injured cord contributes to disruption of the blood-spinal cord barrier, and the influx of leukocytes into the injured cord, as well as apoptosis. MMP-9 and MMP-2 regulate inflammation and neuropathic pain after peripheral nerve injury and may contribute to SCI-induced pain. Early pharmacologic inhibition of MMPs or the gelatinases (MMP-2 and MMP-9) results in an improvement in long-term neurological recovery and is associated with reduced glial scarring and neuropathic pain. During wound healing, gelatinase A (MMP-2) plays a critical role in limiting the formation of an inhibitory glial scar, and mice that are genetically deficient in this protease showed impaired recovery. Together, these findings illustrate complex, temporally distinct roles of MMPs in SCIs. As early gelatinase activity is detrimental, there is an emerging interest in developing gelatinase-targeted therapeutics that would be specifically tailored to the acute injured spinal cord. Thus, we focus this review on the development of selective gelatinase inhibitors

Genome-Wide Analysis of Heteroduplex DNA in Mismatch Repair–Deficient Yeast Cells Reveals Novel Properties of Meiotic Recombination Pathways

Meiotic DNA double-strand breaks (DSBs) initiate crossover (CO) recombination, which is necessary for accurate chromosome segregation, but DSBs may also repair as non-crossovers (NCOs). Multiple recombination pathways with specific intermediates are expected to lead to COs and NCOs. We revisited the mechanisms of meiotic DSB repair and the regulation of CO formation, by conducting a genome-wide analysis of strand-transfer intermediates associated with recombination events. We performed this analysis in a SK1 × S288C Saccharomyces cerevisiae hybrid lacking the mismatch repair (MMR) protein Msh2, to allow efficient detection of heteroduplex DNAs (hDNAs). First, we observed that the anti-recombinogenic activity of MMR is responsible for a 20% drop in CO number, suggesting that in MMR–proficient cells some DSBs are repaired using the sister chromatid as a template when polymorphisms are present. Second, we observed that a large fraction of NCOs were associated with trans–hDNA tracts constrained to a single chromatid. This unexpected finding is compatible with dissolution of double Holliday junctions (dHJs) during repair, and it suggests the existence of a novel control point for CO formation at the level of the dHJ intermediate, in addition to the previously described control point before the dHJ formation step. Finally, we observed that COs are associated with complex hDNA patterns, confirming that the canonical double-strand break repair model is not sufficient to explain the formation of most COs. We propose that multiple factors contribute to the complexity of recombination intermediates. These factors include repair of nicks and double-stranded gaps, template switches between non-sister and sister chromatids, and HJ branch migration. Finally, the good correlation between the strand transfer properties observed in the absence of and in the presence of Msh2 suggests that the intermediates detected in the absence of Msh2 reflect normal intermediates

Rad3ATR Decorates Critical Chromosomal Domains with γH2A to Protect Genome Integrity during S-Phase in Fission Yeast

Schizosaccharomyces pombe Rad3 checkpoint kinase and its human ortholog ATR are essential for maintaining genome integrity in cells treated with genotoxins that damage DNA or arrest replication forks. Rad3 and ATR also function during unperturbed growth, although the events triggering their activation and their critical functions are largely unknown. Here, we use ChIP-on-chip analysis to map genomic loci decorated by phosphorylated histone H2A (γH2A), a Rad3 substrate that establishes a chromatin-based recruitment platform for Crb2 and Brc1 DNA repair/checkpoint proteins. Unexpectedly, γH2A marks a diverse array of genomic features during S-phase, including natural replication fork barriers and a fork breakage site, retrotransposons, heterochromatin in the centromeres and telomeres, and ribosomal RNA (rDNA) repeats. γH2A formation at the centromeres and telomeres is associated with heterochromatin establishment by Clr4 histone methyltransferase. We show that γH2A domains recruit Brc1, a factor involved in repair of damaged replication forks. Brc1 C-terminal BRCT domain binding to γH2A is crucial in the absence of Rqh1Sgs1, a RecQ DNA helicase required for rDNA maintenance whose human homologs are mutated in patients with Werner, Bloom, and Rothmund–Thomson syndromes that are characterized by cancer-predisposition or accelerated aging. We conclude that Rad3 phosphorylates histone H2A to mobilize Brc1 to critical genomic domains during S-phase, and this pathway functions in parallel with Rqh1 DNA helicase in maintaining genome integrity

The Genetic Signatures of Noncoding RNAs

The majority of the genome in animals and plants is transcribed in a developmentally regulated manner to produce large numbers of non–protein-coding RNAs (ncRNAs), whose incidence increases with developmental complexity. There is growing evidence that these transcripts are functional, particularly in the regulation of epigenetic processes, leading to the suggestion that they compose a hitherto hidden layer of genomic programming in humans and other complex organisms. However, to date, very few have been identified in genetic screens. Here I show that this is explicable by an historic emphasis, both phenotypically and technically, on mutations in protein-coding sequences, and by presumptions about the nature of regulatory mutations. Most variations in regulatory sequences produce relatively subtle phenotypic changes, in contrast to mutations in protein-coding sequences that frequently cause catastrophic component failure. Until recently, most mapping projects have focused on protein-coding sequences, and the limited number of identified regulatory mutations have been interpreted as affecting conventional cis-acting promoter and enhancer elements, although these regions are often themselves transcribed. Moreover, ncRNA-directed regulatory circuits underpin most, if not all, complex genetic phenomena in eukaryotes, including RNA interference-related processes such as transcriptional and post-transcriptional gene silencing, position effect variegation, hybrid dysgenesis, chromosome dosage compensation, parental imprinting and allelic exclusion, paramutation, and possibly transvection and transinduction. The next frontier is the identification and functional characterization of the myriad sequence variations that influence quantitative traits, disease susceptibility, and other complex characteristics, which are being shown by genome-wide association studies to lie mostly in noncoding, presumably regulatory, regions. There is every possibility that many of these variations will alter the interactions between regulatory RNAs and their targets, a prospect that should be borne in mind in future functional analyses

Stable transmission of reversible modifications: maintenance of epigenetic information through the cell cycle

Even though every cell in a multicellular organism contains the same genes, the differing spatiotemporal expression of these genes determines the eventual phenotype of a cell. This means that each cell type contains a specific epigenetic program that needs to be replicated through cell divisions, along with the genome, in order to maintain cell identity. The stable inheritance of these programs throughout the cell cycle relies on several epigenetic mechanisms. In this review, DNA methylation and histone methylation by specific histone lysine methyltransferases (KMT) and the Polycomb/Trithorax proteins are considered as the primary mediators of epigenetic inheritance. In addition, non-coding RNAs and nuclear organization are implicated in the stable transfer of epigenetic information. Although most epigenetic modifications are reversible in nature, they can be stably maintained by self-recruitment of modifying protein complexes or maintenance of these complexes or structures through the cell cycle