MIMAS 3.0 is a Multiomics Information Management and Annotation System.

BACKGROUND: DNA sequence integrity, mRNA concentrations and protein-DNA interactions have been subject to genome-wide analyses based on microarrays with ever increasing efficiency and reliability over the past fifteen years. However, very recently novel technologies for Ultra High-Throughput DNA Sequencing (UHTS) have been harnessed to study these phenomena with unprecedented precision. As a consequence, the extensive bioinformatics environment available for array data management, analysis, interpretation and publication must be extended to include these novel sequencing data types. DESCRIPTION: MIMAS was originally conceived as a simple, convenient and local Microarray Information Management and Annotation System focused on GeneChips for expression profiling studies. MIMAS 3.0 enables users to manage data from high-density oligonucleotide SNP Chips, expression arrays (both 3'UTR and tiling) and promoter arrays, BeadArrays as well as UHTS data using MIAME-compliant standardized vocabulary. Importantly, researchers can export data in MAGE-TAB format and upload them to the EBI's ArrayExpress certified data repository using a one-step procedure. CONCLUSION: We have vastly extended the capability of the system such that it processes the data output of six types of GeneChips (Affymetrix), two different BeadArrays for mRNA and miRNA (Illumina) and the Genome Analyzer (a popular Ultra-High Throughput DNA Sequencer, Illumina), without compromising on its flexibility and user-friendliness. MIMAS, appropriately renamed into Multiomics Information Management and Annotation System, is currently used by scientists working in approximately 50 academic laboratories and genomics platforms in Switzerland and France. MIMAS 3.0 is freely available via http://multiomics.sourceforge.net/.

Vertebrate conserved non coding DNA regions have a high persistence length and a short persistence time.

BACKGROUND: The comparison of complete genomes has revealed surprisingly large numbers of conserved non-protein-coding (CNC) DNA regions. However, the biological function of CNC remains elusive. CNC differ in two aspects from conserved protein-coding regions. They are not conserved across phylum boundaries, and they do not contain readily detectable sub-domains. Here we characterize the persistence length and time of CNC and conserved protein-coding regions in the vertebrate and insect lineages. RESULTS: The persistence length is the length of a genome region over which a certain level of sequence identity is consistently maintained. The persistence time is the evolutionary period during which a conserved region evolves under the same selective constraints.Our main findings are: (i) Insect genomes contain 1.60 times less conserved information than vertebrates; (ii) Vertebrate CNC have a higher persistence length than conserved coding regions or insect CNC; (iii) CNC have shorter persistence times as compared to conserved coding regions in both lineages. CONCLUSION: Higher persistence length of vertebrate CNC indicates that the conserved information in vertebrates and insects is organized in functional elements of different lengths. These findings might be related to the higher morphological complexity of vertebrates and give clues about the structure of active CNC elements.Shorter persistence time might explain the previously puzzling observations of highly conserved CNC within each phylum, and of a lack of conservation between phyla. It suggests that CNC divergence might be a key factor in vertebrate evolution. Further evolutionary studies will help to relate individual CNC to specific developmental processes.

SAMHD1 is mutated recurrently in chronic lymphocytic leukemia and is involved in response to DNA damage.

SAMHD1 is a deoxynucleoside triphosphate triphosphohydrolase and a nuclease that restricts HIV-1 in noncycling cells. Germ-line mutations in SAMHD1 have been described in patients with Aicardi-Goutières syndrome (AGS), a congenital autoimmune disease. In a previous longitudinal whole genome sequencing study of chronic lymphocytic leukemia (CLL), we revealed a SAMHD1 mutation as a potential founding event. Here, we describe an AGS patient carrying a pathogenic germ-line SAMHD1 mutation who developed CLL at 24 years of age. Using clinical trial samples, we show that acquired SAMHD1 mutations are associated with high variant allele frequency and reduced SAMHD1 expression and occur in 11% of relapsed/refractory CLL patients. We provide evidence that SAMHD1 regulates cell proliferation and survival and engages in specific protein interactions in response to DNA damage. We propose that SAMHD1 may have a function in DNA repair and that the presence of SAMHD1 mutations in CLL promotes leukemia development.

Regulation of the DNA-binding and transcriptional activities of Xenopus laevis NFI-X by a novel C-terminal domain.

The nuclear factor I (NFI) family consists of sequence-specific DNA-binding proteins that activate both transcription and adenovirus DNA replication. We have characterized three new members of the NFI family that belong to the Xenopus laevis NFI-X subtype and differ in their C-termini. We show that these polypeptides can activate transcription in HeLa and Drosophila Schneider line 2 cells, using an activation domain that is subdivided into adjacent variable and subtype-specific domains each having independent activation properties in chimeric proteins. Together, these two domains constitute the full NFI-X transactivation potential. In addition, we find that the X. laevis NFI-X proteins are capable of activating adenovirus DNA replication through their conserved N-terminal DNA-binding domains. Surprisingly, their in vitro DNA-binding activities are specifically inhibited by a novel repressor domain contained within the C-terminal part, while the dimerization and replication functions per se are not affected. However, inhibition of DNA-binding activity in vitro is relieved within the cell, as transcriptional activation occurs irrespective of the presence of the repressor domain. Moreover, the region comprising the repressor domain participates in transactivation. Mechanisms that may allow the relief of DNA-binding inhibition in vivo and trigger transcriptional activation are discussed.

DNA's golden jubilee.

Review of : 50 Years of DNA Foreword by Phil Campbell, edited by Julie Clayton and Carina Dennis Palgrave Macmillan

Lesion bypass by the Escherichia coli DNA polymerase V requires assembly of a RecA nucleoprotein filament.

Translesion replication is carried out in Escherichia coli by the SOS-inducible DNA polymerase V (UmuC), an error-prone polymerase, which is specialized for replicating through lesions in DNA, leading to the formation of mutations. Lesion bypass by pol V requires the SOS-regulated proteins UmuD' and RecA and the single-strand DNA-binding protein (SSB). Using an in vitro assay system for translesion replication based on a gapped plasmid carrying a site-specific synthetic abasic site, we show that the assembly of a RecA nucleoprotein filament is required for lesion bypass by pol V. This is based on the reaction requirements for stoichiometric amounts of RecA and for single-stranded gaps longer than 100 nucleotides and on direct visualization of RecA-DNA filaments by electron microscopy. SSB is likely to facilitate the assembly of the RecA nucleoprotein filament; however, it has at least one additional role in lesion bypass. ATPgammaS, which is known to strongly increase binding of RecA to DNA, caused a drastic inhibition of pol V activity. Lesion bypass does not require stoichiometric binding of UmuD' along RecA filaments. In summary, the RecA nucleoprotein filament, previously known to be required for SOS induction and homologous recombination, is also a critical intermediate in translesion replication.

The E.coli RuvAB proteins branch migrate Holliday junctions through heterologous DNA sequences in a reaction facilitated by SSB.

During genetic recombination a heteroduplex joint is formed between two homologous DNA molecules. The heteroduplex joint plays an important role in recombination since it accommodates sequence heterogeneities (mismatches, insertions or deletions) that lead to genetic variation. Two Escherichia coli proteins, RuvA and RuvB, promote the formation of heteroduplex DNA by catalysing the branch migration of crossovers, or Holliday junctions, which link recombining chromosomes. We show that RuvA and RuvB can promote branch migration through 1800 bp of heterologous DNA, in a reaction facilitated by the presence of E.coli single-stranded DNA binding (SSB) protein. Reaction intermediates, containing unpaired heteroduplex regions bound by SSB, were directly visualized by electron microscopy. In the absence of SSB, or when SSB was replaced by a single-strand binding protein from bacteriophage T4 (gene 32 protein), only limited heterologous branch migration was observed. These results show that the RuvAB proteins, which are induced as part of the SOS response to DNA damage, allow genetic recombination and the recombinational repair of DNA to occur in the presence of extensive lengths of heterology.

Selecting control genes for RT-QPCR using public microarray data.

Background: Gene expression analysis has emerged as a major biological research area, with real-time quantitative reverse transcription PCR (RT-QPCR) being one of the most accurate and widely used techniques for expression profiling of selected genes. In order to obtain results that are comparable across assays, a stable normalization strategy is required. In general, the normalization of PCR measurements between different samples uses one to several control genes (e. g. housekeeping genes), from which a baseline reference level is constructed. Thus, the choice of the control genes is of utmost importance, yet there is not a generally accepted standard technique for screening a large number of candidates and identifying the best ones. Results: We propose a novel approach for scoring and ranking candidate genes for their suitability as control genes. Our approach relies on publicly available microarray data and allows the combination of multiple data sets originating from different platforms and/or representing different pathologies. The use of microarray data allows the screening of tens of thousands of genes, producing very comprehensive lists of candidates. We also provide two lists of candidate control genes: one which is breast cancer-specific and one with more general applicability. Two genes from the breast cancer list which had not been previously used as control genes are identified and validated by RT-QPCR. Open source R functions are available at http://www.isrec.isb-sib.ch/similar to vpopovic/research/ Conclusion: We proposed a new method for identifying candidate control genes for RT-QPCR which was able to rank thousands of genes according to some predefined suitability criteria and we applied it to the case of breast cancer. We also empirically showed that translating the results from microarray to PCR platform was achievable.

Differential gene expression in response to mechanical wounding and insect feeding in Arabidopsis.

Wounding in multicellular eukaryotes results in marked changes in gene expression that contribute to tissue defense and repair. Using a cDNA microarray technique, we analyzed the timing, dynamics, and regulation of the expression of 150 genes in mechanically wounded leaves of Arabidopsis. Temporal accumulation of a group of transcripts was correlated with the appearance of oxylipin signals of the jasmonate family. Analysis of the coronatine-insensitive coi1-1 Arabidopsis mutant that is also insensitive to jasmonate allowed us to identify a large number of COI1-dependent and COI1-independent wound-inducible genes. Water stress was found to contribute to the regulation of an unexpectedly large fraction of these genes. Comparing the results of mechanical wounding with damage by feeding larvae of the cabbage butterfly (Pieris rapae) resulted in very different transcript profiles. One gene was specifically induced by insect feeding but not by wounding; moreover, there was a relative lack of water stress-induced gene expression during insect feeding. These results help reveal a feeding strategy of P. rapae that may minimize the activation of a subset of water stress-inducible, defense-related genes.

Geometry and physics of catenanes applied to the study of DNA replication.

The concept of ideal geometric configurations was recently applied to the classification and characterization of various knots. Different knots in their ideal form (i.e., the one requiring the shortest length of a constant-diameter tube to form a given knot) were shown to have an overall compactness proportional to the time-averaged compactness of thermally agitated knotted polymers forming corresponding knots. This was useful for predicting the relative speed of electrophoretic migration of different DNA knots. Here we characterize the ideal geometric configurations of catenanes (called links by mathematicians), i.e., closed curves in space that are topologically linked to each other. We demonstrate that the ideal configurations of different catenanes show interrelations very similar to those observed in the ideal configurations of knots. By analyzing literature data on electrophoretic separations of the torus-type of DNA catenanes with increasing complexity, we observed that their electrophoretic migration is roughly proportional to the overall compactness of ideal representations of the corresponding catenanes. This correlation does not apply, however, to electrophoretic migration of certain replication intermediates, believed up to now to represent the simplest torus-type catenanes. We propose, therefore, that freshly replicated circular DNA molecules, in addition to forming regular catenanes, may also form hemicatenanes.

Implication of DNA methylation and PAX5 factor in the transcriptional regulation of hTERT, the human telomerase reverse transcriptase gene

RESUME La télomérase est une enzyme dite "d'immortalité" qui permet aux cellules de maintenir la longueur de leurs télomères, ce qui confère une capacité de réplication illimitée aux cellules reproductrices et cancéreuses. A l'inverse, les cellules somatiques normales, qui n'expriment pas la télomérase, ont une capacité de réplication limitée. La sous-unité catalytique de la télomérase, hTERT, est définie comme le facteur limitant l'activité télomérasique. Entre activateurs et répresseurs, le rôle de la méthylation de l'ADN et de l'acétylation des histones, de nombreux modèles ont été suggérés. La découverte de l'implication de CTCF dans la régulation transcriptionnelle de hTERT explique en partie le mécanisme de répression de la télomérase dans la plupart des cellules somatiques et sa réactivation dans les cellules tumorales. Dans les cellules télomérase-positives, l'activité inhibitrice de CTCF est bloquée par un mécanisme dépendent ou non de la méthylation. Dans la plupart des carcinomes, une hyperméthylation de la région 5' de hTERT bloque l'effet inhibiteur de CTCF, alors qu'une petite région hypométhylée permet un faible niveau de transcription du gène. Nous avons démontré que la protéine MBD2 se lie spécifiquement sur la région 5' méthylée de hTERT dans différentes lignées cellulaires et qu'elle est impliquée dans la répression partielle de la transcription de hTERT dans les cellules tumorales méthylées. Par contre, nous avons montré que dans les lymphocytes B normaux et néoplasiques, la régulation de hTERT est indépendante de la méthylation. Dans ces cellules, le facteur PAX5 se lie sur la région 5' de hTERT en aval du site d'initiation de la traduction (ATG). L'expression exogène de PAX5 dans les cellules télomérase-négatives active la transcription de hTERT, alors que la répression de PAX5 dans les cellules lymphomateuses inhibe la transcription du gène. PAX5 est donc directement impliqué dans l'activation de l'expression de hTERT dans les lymphocytes B exprimant la télomérase. Ces résultats révèlent des différences entre les niveaux de méthylation de hTERT dans les cellules de carcinomes et les lymphocytes B exprimant la télomérase. La méthylation de hTERT en tant que biomarqueur de cancer a été évaluée, puis appliquée à la détection de métastases. Nous avons ainsi montré que la méthylation de hTERT est positivement corrélée au diagnostic cytologique dans les liquides céphalorachidiens. Nos résultats conduisent à un modèle de régulation de hTERT, qui aide à comprendre comment la transcription de ce gène est régulée par CTCF, avec un mécanisme lié ou non à la méthylation du gène hTERT. La méthylation de hTERT s'est aussi révélée être un nouveau et prometteur biomarqueur de cancer. SUMMARY Human telomerase is an "immortalizing" enzyme that enables cells to maintain telomere length, allowing unlimited replicative capacity to reproductive and cancer cells. Conversely, normal somatic cells that do not express telomerase have a finite replicative capacity. The catalytic subunit of telomerase, hTERT, is defined as the limiting factor for telomerase activity. Between activators and repressors, and the role of DNA methylation and histone acetylation, an abundance of hTERT regulatory models have been suggested. The discovery of the implication of CTCF in the transcriptional regulation of hTERT in part explained the mechanism of silencing of telomerase in most somatic cells and its reactivation in neoplastic cells. In telomerase-positive cells, the inhibitory activity of CTCF is blocked by methylation-dependent and -independent mechanisms. In most carcinoma cells, hypermethylation of the hTERT 5' region has been shown to block the inhibitory effect of CTCF, while a short hypomethylated region allows a low transcription level of the gene. We have demonstrated that MBD2 protein specifically binds the methylated 5' region of hTERT in different cell lines and is therefore involved in the partial repression of hTERT transcription in methylated tumor cells. In contrast, we have shown that in normal and neoplastic B cells, hTERT regulation is methylation-independent. The PAX5 factor has been shown to bind to the hTERT 5'region downstream of the ATG translational start site. Ectopic expression of PAX5 in telomerase-negative cells or repression of PAX5 expression in B lymphoma cells respectively activated and repressed hTERT transcription. Thus, PAX5 is strongly implicated in hTERT expression activation in telomerase-positive B cells. These results reveal differences between the hTERT methylation patterns in telomerase-positive carcinoma cells and telomerase-positive normal B cells. The potential of hTERT methylation as a cancer biomarker was evaluated and applied to the detection of metastasis. We have shown that hTERT methylation correlates with the cytological diagnosis in cerebrospinal fluids. Our results suggest a model of hTERT gene regulation, which helps us to better understand how hTERT transcription is regulated by CTCF in methylation-dependant and independent mechanisms. Our data also indicate that hTERT methylation is a promising new cancer biomarker.

Sleep loss reduces the DNA-binding of BMAL1, CLOCK, and NPAS2 to specific clock genes in the mouse cerebral cortex.

We have previously demonstrated that clock genes contribute to the homeostatic aspect of sleep regulation. Indeed, mutations in some clock genes modify the markers of sleep homeostasis and an increase in homeostatic sleep drive alters clock gene expression in the forebrain. Here, we investigate a possible mechanism by which sleep deprivation (SD) could alter clock gene expression by quantifying DNA-binding of the core-clock transcription factors CLOCK, NPAS2, and BMAL1 to the cis-regulatory sequences of target clock genes in mice. Using chromatin immunoprecipitation (ChIP), we first showed that, as reported for the liver, DNA-binding of CLOCK and BMAL1 to target clock genes changes in function of time-of-day in the cerebral cortex. Tissue extracts were collected at ZT0 (light onset), -6, -12, and -18, and DNA enrichment of E-box or E'-box containing sequences was measured by qPCR. CLOCK and BMAL1 binding to Cry1, Dbp, Per1, and Per2 depended on time-of-day, with maximum values reached at around ZT6. We then observed that SD, performed between ZT0 and -6, significantly decreased DNA-binding of CLOCK and BMAL1 to Dbp, consistent with the observed decrease in Dbp mRNA levels after SD. The DNA-binding of NPAS2 and BMAL1 to Per2 was also decreased by SD, although SD is known to increase Per2 expression in the cortex. DNA-binding to Per1 and Cry1 was not affected by SD. Our results show that the sleep-wake history can affect the clock molecular machinery directly at the level of chromatin binding thereby altering the cortical expression of Dbp and Per2 and likely other targets. Although the precise dynamics of the relationship between DNA-binding and mRNA expression, especially for Per2, remains elusive, the results also suggest that part of the reported circadian changes in DNA-binding of core clock components in tissues peripheral to the suprachiasmatic nuclei could, in fact, be sleep-wake driven.

Impact of normalization on miRNA microarray expression profiling.

Profiling miRNA levels in cells with miRNA microarrays is becoming a widely used technique. Although normalization methods for mRNA gene expression arrays are well established, miRNA array normalization has so far not been investigated in detail. In this study we investigate the impact of normalization on data generated with the Agilent miRNA array platform. We have developed a method to select nonchanging miRNAs (invariants) and use them to compute linear regression normalization coefficients or variance stabilizing normalization (VSN) parameters. We compared the invariants normalization to normalization by scaling, quantile, and VSN with default parameters as well as to no normalization using samples with strong differential expression of miRNAs (heart-brain comparison) and samples where only a few miRNAs are affected (by p53 overexpression in squamous carcinoma cells versus control). All normalization methods performed better than no normalization. Normalization procedures based on the set of invariants and quantile were the most robust over all experimental conditions tested. Our method of invariant selection and normalization is not limited to Agilent miRNA arrays and can be applied to other data sets including those from one color miRNA microarray platforms, focused gene expression arrays, and gene expression analysis using quantitative PCR.

Identification of estrogen-responsive DNA sequences by transient expression experiments in a human breast cancer cell line.

The expression of a hybrid gene formed by the promoter region of the Xenopus laevis vitellogenin gene B1 and the CAT coding region is regulated by estrogen when the gene is transfected into hormone-responsive MCF-7 cells. Furthermore, the 5' flanking region of the gene B1 alone can confer inducibility to heterologous promoters, although to a varying extent depending on the promoter used. Deletion mapping of he vitellogenin hormone-responsive sequences revealed that a 13 bp element 5'-AGTCACTGTGACC-3' at position -334 is essential for estrogen inducibility. We have shown previously that this 13 bp element is present upstream of several liver-specific estrogen-inducible genes.