Design and analysis of ChIP-Seq experiments for nuclear factor I DNA-binding proteins

SUMMARY : Eukaryotic DNA interacts with the nuclear proteins using non-covalent ionic interactions. Proteins can recognize specific nucleotide sequences based on the sterical interactions with the DNA and these specific protein-DNA interactions are the basis for many nuclear processes, e.g. gene transcription, chromosomal replication, and recombination. New technology termed ChIP-Seq has been recently developed for the analysis of protein-DNA interactions on a whole genome scale and it is based on immunoprecipitation of chromatin and high-throughput DNA sequencing procedure. ChIP-Seq is a novel technique with a great potential to replace older techniques for mapping of protein-DNA interactions. In this thesis, we bring some new insights into the ChIP-Seq data analysis. First, we point out to some common and so far unknown artifacts of the method. Sequence tag distribution in the genome does not follow uniform distribution and we have found extreme hot-spots of tag accumulation over specific loci in the human and mouse genomes. These artifactual sequence tags accumulations will create false peaks in every ChIP-Seq dataset and we propose different filtering methods to reduce the number of false positives. Next, we propose random sampling as a powerful analytical tool in the ChIP-Seq data analysis that could be used to infer biological knowledge from the massive ChIP-Seq datasets. We created unbiased random sampling algorithm and we used this methodology to reveal some of the important biological properties of Nuclear Factor I DNA binding proteins. Finally, by analyzing the ChIP-Seq data in detail, we revealed that Nuclear Factor I transcription factors mainly act as activators of transcription, and that they are associated with specific chromatin modifications that are markers of open chromatin. We speculate that NFI factors only interact with the DNA wrapped around the nucleosome. We also found multiple loci that indicate possible chromatin barrier activity of NFI proteins, which could suggest the use of NFI binding sequences as chromatin insulators in biotechnology applications. RESUME : L'ADN des eucaryotes interagit avec les protéines nucléaires par des interactions noncovalentes ioniques. Les protéines peuvent reconnaître les séquences nucléotidiques spécifiques basées sur l'interaction stérique avec l'ADN, et des interactions spécifiques contrôlent de nombreux processus nucléaire, p.ex. transcription du gène, la réplication chromosomique, et la recombinaison. Une nouvelle technologie appelée ChIP-Seq a été récemment développée pour l'analyse des interactions protéine-ADN à l'échelle du génome entier et cette approche est basée sur l'immuno-précipitation de la chromatine et sur la procédure de séquençage de l'ADN à haut débit. La nouvelle approche ChIP-Seq a donc un fort potentiel pour remplacer les anciennes techniques de cartographie des interactions protéine-ADN. Dans cette thèse, nous apportons de nouvelles perspectives dans l'analyse des données ChIP-Seq. Tout d'abord, nous avons identifié des artefacts très communs associés à cette méthode qui étaient jusqu'à présent insoupçonnés. La distribution des séquences dans le génome ne suit pas une distribution uniforme et nous avons constaté des positions extrêmes d'accumulation de séquence à des régions spécifiques, des génomes humains et de la souris. Ces accumulations des séquences artéfactuelles créera de faux pics dans toutes les données ChIP-Seq, et nous proposons différentes méthodes de filtrage pour réduire le nombre de faux positifs. Ensuite, nous proposons un nouvel échantillonnage aléatoire comme un outil puissant d'analyse des données ChIP-Seq, ce qui pourraient augmenter l'acquisition de connaissances biologiques à partir des données ChIP-Seq. Nous avons créé un algorithme d'échantillonnage aléatoire et nous avons utilisé cette méthode pour révéler certaines des propriétés biologiques importantes de protéines liant à l'ADN nommés Facteur Nucléaire I (NFI). Enfin, en analysant en détail les données de ChIP-Seq pour la famille de facteurs de transcription nommés Facteur Nucléaire I, nous avons révélé que ces protéines agissent principalement comme des activateurs de transcription, et qu'elles sont associées à des modifications de la chromatine spécifiques qui sont des marqueurs de la chromatine ouverte. Nous pensons que lés facteurs NFI interagir uniquement avec l'ADN enroulé autour du nucléosome. Nous avons également constaté plusieurs régions génomiques qui indiquent une éventuelle activité de barrière chromatinienne des protéines NFI, ce qui pourrait suggérer l'utilisation de séquences de liaison NFI comme séquences isolatrices dans des applications de la biotechnologie.

Multiple adenosine 3',5'-cyclic [corrected] monophosphate response element DNA-binding proteins generated by gene diversification and alternative exon splicing.

The protein sequence deduced from the open reading frame of a human placental cDNA encoding a cAMP-responsive enhancer (CRE)-binding protein (CREB-327) has structural features characteristic of several other transcriptional transactivator proteins including jun, fos, C/EBP, myc, and CRE-BP1. Results of Southwestern analysis of nuclear extracts from several different cell lines show that there are multiple CRE-binding proteins, which vary in size in cell lines derived from different tissues and animal species. To examine the molecular diversity of CREB-327 and related proteins at the nucleic acid level, we used labeled cDNAs from human placenta that encode two different CRE-binding proteins (CREB-327 and CRE-BP1) to probe Northern and Southern blots. Both probes hybridized to multiple fragments on Southern blots of genomic DNA from various species. Alternatively, when a human placental c-jun probe was hybridized to the same blot, a single fragment was detected in most cases, consistent with the intronless nature of the human c-jun gene. The CREB-327 probe hybridized to multiple mRNAs, derived from human placenta, ranging in size from 2-9 kilobases. In contrast, the CRE-BP1 probe identified a single 4-kilobase mRNA. Sequence analyses of several overlapping human genomic cosmid clones containing CREB-327 sequences in conjunction with polymerase chain reaction indicates that the CREB-327/341 cDNAs are composed of at least eight or nine exons, and analyses of human placental cDNAs provide direct evidence for at least one alternatively spliced exon. Analyses of mouse/hamster-human hybridoma DNAs by Southern blotting and polymerase chain reaction localizes the CREB-327/341 gene to human chromosome 2. The results indicate that there is a dichotomy of CREB-like proteins, those that are related by overall structure and DNA-binding specificity as well as those that are related by close similarities of primary sequences.

Functional analysis of human POT1 and RPA : insights into the role of single stranded DNA binding proteins at telomeres

Abstract Telomeres, the natural ends of chromosomes, need to be protected from chromosome end fusions, aberrant homologous recombination and degradation. In humans, chromosome ends are specified through arrays of tandemly repeated 5'-TTAGGG-3' hexamers, ending in a 3' overhang. A complex formed by the six proteins TRF1, TRF2, hRap1, TIN2, TPP1 and POT1 specifically assocìates with and protects telomeres. Telomeres are maintained by semiconservative DNA replication and by a specialized reverse transcriptase, telomerase, that carries an RNA subunit which templates new telomeric repeat synthesis. The telomeric single stranded (ss) DNA binding protein POT1 protects the telomeric 3' overhang and modulates telomerase-mediated telomere elongation. It is possible that POT1 also influences DNA synthesis during semiconservative DNA replication, which is initiated by the DNA polymerase alpha-primase complex. The heterotrimeric ss DNA-binding protein RPA plays essential roles during DNA replication. RPA binds to ss DNA with high affinity in order to stabilize ss DNA and facilitate nascent strand synthesis at the replication fork. Here we investigate how the two proteins RPA and POT1 contribute to telomere maintenance by regulating semi-conservative DNA replication and telomerase. Using chromatin immunoprecipitation experiments, we show that RPA associates with telomeres during S-phase. Analysis of telomere structure in cells shRNA-depleted for RPA and POT1 reveals that loss of RPA and POT1 causes exposure of single-stranded DNA at telomeres, suggestive of incomplete DNA replication. Biochemical experiments using purified recombinant POT1 and RPA show that saturating telomeric oligonucleotides with POT1 or RPA reduces the primase activity of the DNA polymerase alpha-primase complex and the overall activity of telomerase. POT1 and RPA also increase the primer extension by DNA polymerase alpha-primase complex and the processivity of telomerase under certain conditions, although POT1 increases the activities to a greater extent than RPA. We propose that POT1 is required for proper replication of the lagging strand of telomeres and that some phenotypes observed in POT1-depleted cells may stern from incomplete DNA replication rather than de-protection of the single-stranded overhang. Résumé Les télomères, les extrémités normales des chromosomes linéaires, doivent être protégés des fusions chromosomiques, d'événements de recombinaison homologue aberrants et de phénomènes de dégradation. Chez l'Homme, les extrémités des chromosomes sont constitués d'ADN double brin répétitif de séquence 5'-TTAGGG-3', d'une extension simple brin 3' sortante et d'un complexe protéique formé des six facteurs TRF1, TRF2, hRap1, TIN2, TPP1 et POT1 qui, s'associant à cette séquence, protègent l'ADN télomèrique. Les télomères sont maintenus par la télomérase, une transcriptase inverse capable d'allonger l'extension 3' sortante télomérique. POT1 lie l'ADN simple brin télomérique et module l'élongation des télomères par la télomérase. POT1 pourrait en théorie également influencer la réplication semi-conservative de l'ADN. L'ADN-polymérase Pal alpha-primase amorce et initie la synthèse d'ADN. Pendant la réplication, l'ADN simple brin est stabilisé par RPA, un complexe hétérotrimèrique qui lie l'ADN simple brin. RPA facilite la synthèse du brin naissant à la fourche de réplication. Ici nous avons étudié comment ces deux protéines qui lient l'ADN simple brin, RPA et POT1, régulent la réplication des télomères par la télomérase et la machinerie classique de réplication de l'ADN. Par immunoprécipitation de chromatine (ChIP), nous montrons que RPA est localisé aux télomères lors de la phase S du cycle cellulaire. De plus, l'analyse de la structure des télomeres indique que !a perte de RPA ou de POT1 conduit à l'apparition d'ADN simple brin télomérique, suggérant une réplication incomplète de l'ADN télomérique in vivo. Par une approche complémentaire biochimique utilisant les protéines POT1 et RPA recombinantes purifiées, nous montrons également que la liaison de POT1 ou de RPA à des oligonucléotides télomériques bloque l'activité primase du complexe polymérase alpha/primase et réduit l'activité télomérase sur ces substrats. En revanche, leur liaison augmente l'activité ADN-polymérase du complexe polymérase alpha/primase, ainsi que fa processivité de la télomérase dans certaines conditions, POT1 étant le plus efficace des deux facteurs. Nous proposons que POT1 est nécessaire à la réplication du brin retardé au niveau des télomères, ce qui suggère que certains phénotypes des cellules déplétés en POT1 puissent résulter d'une réplication incomplète de l'ADN télémétrique plutôt que d'une déprotection de l'extrémité sortante des télomères.

DNA binding properties of peroxisome proliferator-activated receptor subtypes on various natural peroxisome proliferator response elements. Importance of the 5'-flanking region.

The three subtypes of the peroxisome proliferator-activated receptors (PPARalpha, beta/delta, and gamma) form heterodimers with the 9-cis-retinoic acid receptor (RXR) and bind to a common consensus response element, which consists of a direct repeat of two hexanucleotides spaced by one nucleotide (DR1). As a first step toward understanding the molecular mechanisms determining PPAR subtype specificity, we evaluated by electrophoretic mobility shift assays the binding properties of the three PPAR subtypes, in association with either RXRalpha or RXRgamma, on 16 natural PPAR response elements (PPREs). The main results are as follows. (i) PPARgamma in combination with either RXRalpha or RXRgamma binds more strongly than PPARalpha or PPARbeta to all natural PPREs tested. (ii) The binding of PPAR to strong elements is reinforced if the heterodimerization partner is RXRgamma. In contrast, weak elements favor RXRalpha as heterodimerization partner. (iii) The ordering of the 16 natural PPREs from strong to weak elements does not depend on the core DR1 sequence, which has a relatively uniform degree of conservation, but correlates with the number of identities of the 5'-flanking nucleotides with respect to a consensus element. This 5'-flanking sequence is essential for PPARalpha binding and thus contributes to subtype specificity. As a demonstration of this, the PPARgamma-specific element ARE6 PPRE is able to bind PPARalpha only if its 5'-flanking region is exchanged with that of the more promiscuous HMG PPRE.

DNA binding specificity of different STAT proteins. Comparison of in vitro specificity with natural target sites.

STAT transcription factors are expressed in many cell types and bind to similar sequences. However, different STAT gene knock-outs show very distinct phenotypes. To determine whether differences between the binding specificities of STAT proteins account for these effects, we compared the sequences bound by STAT1, STAT5A, STAT5B, and STAT6. One sequence set was selected from random oligonucleotides by recombinant STAT1, STAT5A, or STAT6. For another set including many weak binding sites, we quantified the relative affinities to STAT1, STAT5A, STAT5B, and STAT6. We compared the results to the binding sites in natural STAT target genes identified by others. The experiments confirmed the similar specificity of different STAT proteins. Detailed analysis indicated that STAT5A specificity is more similar to that of STAT6 than that of STAT1, as expected from the evolutionary relationships. The preference of STAT6 for sites in which the half-palindromes (TTC) are separated by four nucleotides (N(4)) was confirmed, but analysis of weak binding sites showed that STAT6 binds fairly well to N(3) sites. As previously reported, STAT1 and STAT5 prefer N(3) sites; however, STAT5A, but not STAT1, weakly binds N(4) sites. None of the STATs bound to half-palindromes. There were no specificity differences between STAT5A and STAT5B.

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.

The N-terminal DNA-binding 'zinc finger' of the oestrogen and glucocorticoid receptors determines target gene specificity.

Steroid hormone receptors activate specific gene transcription by binding as hormone-receptor complexes to short DNA enhancer-like elements termed hormone response elements (HREs). We have shown previously that a highly conserved 66 amino acid region of the oestrogen (ER) and glucocorticoid (GR) receptors, which corresponds to part of the receptor DNA binding domain (region C) is responsible for determining the specificity of target gene activation. This region contains two sub-regions (CI and CII) analogous to the 'zinc-fingers' of the transcription factor TFIIIA. We show here that CI and CII appear to be separate domains both involved in DNA binding. Furthermore, using chimaeric ERs in which either the first (N-terminal) (CI) or second (CII) 'zinc finger' region has been exchanged with that of the GR, indicates that it is the first 'zinc finger' which largely determines target gene specificity. We suggest that receptor recognition of the HRE is analogous to that of the helix-turn-helix DNA binding motif in that the receptor binds to DNA as a dimer with the first 'zinc finger' lying in the major groove recognizing one half of the palindromic HRE, and that protein-DNA interaction is stabilized through non-specific DNA binding and dimer interactions contributed by the second 'zinc finger'.

The loss of GLUT2 expression in the pancreatic beta-cells of diabetic db/db mice is associated with an impaired DNA-binding activity of islet-specific trans-acting factors.

GLUT2 expression is reduced in the pancreatic beta-cells of several diabetic animals. The transcriptional control of the gene in beta-cells involves at least two islet-specific DNA-binding proteins, GTIIa and PDX-1, which also transactivates the insulin, somatostatin and glucokinase genes. In this report, we assessed the DNA-binding activities of GTIIa and PDX-1 to their respective cis-elements of the GLUT2 promoter using nuclear extracts prepared from pancreatic islets of 12 week old db/db diabetic mice. We show that the decreased GLUT2 mRNA expression correlates with a decrease of the GTIIa DNA-binding activity, whereas the PDX-1 binding activity is increased. In these diabetic animals, insulin mRNA expression remains normal. The adjunction of dexamethasone to isolated pancreatic islets, a treatment previously shown to decrease PDX-1 expression in the insulin-secreting HIT-T15 cells, has no effect on the GTIIa and PDX-1 DNA-binding activities. These data suggest that the decreased activity of GTIIa, in contrast to PDX-1, may be a major initial step in the development of the beta-cell dysfunction in this model of diabetes.

Proof of concept for microarray-based detection of DNA-binding oncogenes in cell extracts.

The function of DNA-binding proteins is controlled not just by their abundance, but mainly at the level of their activity in terms of their interactions with DNA and protein targets. Moreover, the affinity of such transcription factors to their target sequences is often controlled by co-factors and/or modifications that are not easily assessed from biological samples. Here, we describe a scalable method for monitoring protein-DNA interactions on a microarray surface. This approach was designed to determine the DNA-binding activity of proteins in crude cell extracts, complementing conventional expression profiling arrays. Enzymatic labeling of DNA enables direct normalization of the protein binding to the microarray, allowing the estimation of relative binding affinities. Using DNA sequences covering a range of affinities, we show that the new microarray-based method yields binding strength estimates similar to low-throughput gel mobility-shift assays. The microarray is also of high sensitivity, as it allows the detection of a rare DNA-binding protein from breast cancer cells, the human tumor suppressor AP-2. This approach thus mediates precise and robust assessment of the activity of DNA-binding proteins and takes present DNA-binding assays to a high throughput level.

Biochemical properties of RuvBD113N: a mutation in helicase motif II of the RuvB hexamer affects DNA binding and ATPase activities.

Many DNA helicases utilise the energy derived from nucleoside triphosphate hydrolysis to fuel their actions as molecular motors in a variety of biological processes. In association with RuvA, the E. coli RuvB protein (a hexameric ring helicase), promotes the branch migration of Holliday junctions during genetic recombination and DNA repair. To analyse the relationship between ATP-dependent DNA helicase activity and branch migration, a site-directed mutation was introduced into the helicase II motif of RuvB. Over-expression of RuvBD113N in wild-type E. coli resulted in a dominant negative UVs phenotype. The biochemical properties of RuvBD113N were examined and compared with wild-type RuvB in vitro. The single amino acid substitution resulted in major alterations to the biochemical activities of RuvB, such that RuvBD113N was defective in DNA binding and ATP hydrolysis, while retaining the ability to form hexameric rings and interact with RuvA. RuvBD113N formed heterohexamers with wild-type RuvB, and could inhibit RuvB function by affecting its ability to bind DNA. However, heterohexamers exhibited an ability to promote branch migration in vitro indicating that not all subunits of the ring need to be catalytically competent.

Experimental analysis and computer prediction of CTF/NFI transcription factor DNA binding sites.

Accurate prediction of transcription factor binding sites is needed to unravel the function and regulation of genes discovered in genome sequencing projects. To evaluate current computer prediction tools, we have begun a systematic study of the sequence-specific DNA-binding of a transcription factor belonging to the CTF/NFI family. Using a systematic collection of rationally designed oligonucleotides combined with an in vitro DNA binding assay, we found that the sequence specificity of this protein cannot be represented by a simple consensus sequence or weight matrix. For instance, CTF/NFI uses a flexible DNA binding mode that allows for variations of the binding site length. From the experimental data, we derived a novel prediction method using a generalised profile as a binding site predictor. Experimental evaluation of the generalised profile indicated that it accurately predicts the binding affinity of the transcription factor to natural or synthetic DNA sequences. Furthermore, the in vitro measured binding affinities of a subset of oligonucleotides were found to correlate with their transcriptional activities in transfected cells. The combined computational-experimental approach exemplified in this work thus resulted in an accurate prediction method for CTF/NFI binding sites potentially functioning as regulatory regions in vivo.

Characterization of HbpR binding by site-directed mutagenesis of its DNA-binding site and by deletion of the effector domain.

In the presence of 2-hydroxybiphenyl, the enhancer binding protein, HbpR, activates the sigma54-dependent P(hbpC) promoter and controls the initial steps of 2-hydroxybiphenyl degradation in Pseudomonas azelaica. In the activation process, an oligomeric HbpR complex of unknown subunit composition binds to an operator region containing two imperfect palindromic sequences. Here, the HbpR-DNA binding interactions were investigated by site-directed mutagenesis of the operator region and by DNA-binding assays using purified HbpR. Mutations that disrupted the twofold symmetry in the palindromes did not affect the binding affinity of HbpR, but various mutations along a 60 bp region, and also outside the direct palindromic sequences, decreased the binding affinity. Footprints of HbpR on mutant operator fragments showed that a partial loss of binding contacts occurs, suggesting that the binding of one HbpR 'protomer' in the oligomeric complex is impaired whilst leaving the other contacts intact. An HbpR variant, devoid of its N-terminal sensing A-domain, was unable to activate transcription from the hbpC promoter while maintaining protection of the operator DNA in footprints. Wild-type HbpR was unable to activate transcription from the hbpC promoter when delta A-HbpR was expressed in the same cell, suggesting the formation of (repressing) hetero-oligomers. This model implies that HbpR can self-associate on its operator DNA without effector recognition or ATP binding. Furthermore, our findings suggest that the N-terminal sensing domain of HbpR is needed to activate the central ATPase domain rather than to repress a constitutively active C domain, as is the case for the related regulatory protein XylR.

Postnatal expression pattern of calcium-binding proteins in organotypic thalamic cultures and in the dorsal thalamus in vivo.

The present study describes the postnatal expression of calbindin, calretinin and parvalbumin and glutamic acid decarboxylase (GAD) and microtubule-associated protein 2 (MAP2) in organotypic monocultures of rat dorsal thalamus compared to the thalamus in vivo. Cultures were maintained for up to 7 weeks. Cortex-conditioned medium improved the survival of thalamic cultures. MAP2-immunoreactive material was present in somata and dendrites of small and large-sized neurons throughout the cultures. Parvalbumin immunoreactivity was present in larger multipolar or bitufted neurons along the edge of a culture. These neurons also displayed strong parvalbumin mRNA and GAD mRNA expression, and GABA immunoreactivity. They likely corresponded to cells of the nucleus reticularis thalami. Parvalbumin mRNA, but neither parvalbumin protein nor GAD mRNA, was expressed in neurons with large somata within the explant. They likely represented relay cells. GAD mRNA, but not parvalbumin mRNA, was expressed in small neurons within the explants. Small neurons also displayed calbindin- and calretinin-immunoreactivity. The small neurons likely represented local circuit neurons. The time course of expression of the calcium-binding proteins revealed that all were present at birth with the predicted molecular weights. A low, but constant parvalbumin expression was observed in vitro without the developmental increase seen in vivo, which most likely represented parvalbumin from afferent sources. In contrast, the explantation transiently downregulated the calretinin and calbindin expression, but the neurons recovered the expression after 14 and 21 days, respectively. In conclusion, thalamic monocultures older than three weeks represent a stable neuronal network containing well differentiated neurons of the nucleus reticularis thalami, relay cells and local circuit neurons.

EZH2 is essential for glioblastoma cancer stem cell maintenance.

Overexpression of the polycomb group protein enhancer of zeste homologue 2 (EZH2) occurs in diverse malignancies, including prostate cancer, breast cancer, and glioblastoma multiforme (GBM). Based on its ability to modulate transcription of key genes implicated in cell cycle control, DNA repair, and cell differentiation, EZH2 is believed to play a crucial role in tissue-specific stem cell maintenance and tumor development. Here, we show that targeted pharmacologic disruption of EZH2 by the S-adenosylhomocysteine hydrolase inhibitor 3-deazaneplanocin A (DZNep), or its specific downregulation by short hairpin RNA (shRNA), strongly impairs GBM cancer stem cell (CSC) self-renewal in vitro and tumor-initiating capacity in vivo. Using genome-wide expression analysis of DZNep-treated GBM CSCs, we found the expression of c-myc, recently reported to be essential for GBM CSCs, to be strongly repressed upon EZH2 depletion. Specific shRNA-mediated downregulation of EZH2 in combination with chromatin immunoprecipitation experiments revealed that c-myc is a direct target of EZH2 in GBM CSCs. Taken together, our observations provide evidence that direct transcriptional regulation of c-myc by EZH2 may constitute a novel mechanism underlying GBM CSC maintenance and suggest that EZH2 may be a valuable new therapeutic target for GBM management.