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.

Transportin regulates nuclear import of CD44.

CD44 is a facultative cell surface proteoglycan that serves as the principal cell surface receptor for hyaluronan (HA). Studies have shown that in addition to participating in numerous signaling pathways, CD44 becomes internalized upon engagement by ligand and that a portion of its intracellular domain can translocate to the nucleus where it is believed to play a functional role in cell proliferation and survival. However, the mechanisms whereby fragments of CD44 enter the nucleus have not been elucidated. Here we show that CD44 interacts with two import receptors of the importin β superfamily, importin β itself and transportin. Inhibition of importin β-dependent transport failed to block CD44 accumulation in the nucleus. By contrast, inhibition of the transportin-dependent pathway abrogated CD44 import. Mutagenesis of the intracellular domain of CD44 revealed that the 20 membrane-proximal residues contain sequences required for transportin-mediated nuclear transport. Our observations provide evidence that CD44 interacts with importin family members and identify the transportin-dependent pathway as the mechanism whereby full-length CD44 enters the nucleus.

Selection, Recombination and History in a Parasitic Flatworm (Echinococcus) Inferred from Nucleotide Sequences

Three species of flatworms from the genus Echinococcus (E. granulosus, E. multilocularis and E. vogeli) and four strains of E. granulosus (cattle, horse, pig and sheep strains) were analysed by the PCR-SSCP method followed by sequencing, using as targets two non-coding and two coding (one nuclear and one mitochondrial) genomic regions. The sequencing data was used to evaluate hypothesis about the parasite breeding system and the causes of genetic diversification. The calculated recombination parameters suggested that cross-fertilisation was rare in the history of the group. However, the relative rates of substitution in the coding sequences showed that positive selection (instead of purifying selection) drove the evolution of an elastase and neutrophil chemotaxis inhibitor gene (AgB/1). The phylogenetic analyses revealed several ambiguities, indicating that the taxonomic status of the E. granulosus horse strain should be revised

Nuclear rDNA-based molecular clock of the evolution of triatominae (Hemiptera: Reduviidae), vectors of Chagas disease

The evolutionary history and times of divergence of triatomine bug lineages are estimated from molecular clocks inferred from nucleotide sequences of the small subunit SSU (18S) and the second internal transcribed spacer (ITS-2) of the nuclear ribosomal DNA of these reduviids. The 18S rDNA molecular clock rate in Triatominae, and Prosorrhynchan Hemiptera in general, appears to be of 1.8% per 100 million years (my). The ITS-2 molecular clock rate in Triatominae is estimated to be around 0.4-1% per 1 my, indicating that ITS-2 evolves 23-55 times faster than 18S rDNA. Inferred chronological data about the evolution of Triatominae fit well with current hypotheses on their evolutionary histories, but suggest reconsideration of the current taxonomy of North American species complexes.

The second internal transcribed spacer of nuclear ribosomal DNA as a tool for Latin American anopheline taxonomy: a critical review

Among the molecular markers commonly used for mosquito taxonomy, the internal transcribed spacer 2 (ITS2) of the ribosomal DNA is useful for distinguishing among closely-related species. Here we review 178 GenBank accession numbers matching ITS2 sequences of Latin American anophelines. Among those, we found 105 unique sequences corresponding to 35 species. Overall the ITS2 sequences distinguish anopheline species, however, information on intraspecific and geographic variations is scarce. Intraspecific variations ranged from 0.2% to 19% and our analysis indicates that misidentification and/or sequencing errors could be responsible for some of the high values of divergence. Research in Latin American malaria vector taxonomy profited from molecular data provided by single or few field capture mosquitoes. However we propose that caution should be taken and minimum requirements considered in the design of additional studies. Future studies in this field should consider that: (1) voucher specimens, assigned to the DNA sequences, need to be deposited in collections, (2) intraspecific variations should be thoroughly evaluated, (3) ITS2 and other molecular markers, considered as a group, will provide more reliable information, (4) biological data about vector populations are missing and should be prioritized, (5) the molecular markers are most powerful when coupled with traditional taxonomic tools.

Cloning and molecular characterization of Trypanosoma cruzi U2, U4, U5, and U6 small nuclear RNAs

Small nuclear RNAs (snRNAs) are important factors in the functioning of eukaryotic cells that form several small complexes with proteins; these ribonucleoprotein particles (U snRNPs) have an essential role in the pre-mRNA processing, particularly in splicing, catalyzed by spliceosomes, large RNA-protein complexes composed of various snRNPs. Even though they are well defined in mammals, snRNPs are still not totally characterized in certain trypanosomatids as Trypanosoma cruzi. For this reason we subjected snRNAs (U2, U4, U5, and U6) from T. cruzi epimastigotes to molecular characterization by polymerase chain reaction (PCR) and reverse transcription-PCR. These amplified sequences were cloned, sequenced, and compared with those other of trypanosomatids. Among these snRNAs, U5 was less conserved and U6 the most conserved. Their respective secondary structures were predicted and compared with known T. brucei structures. In addition, the copy number of each snRNA in the T. cruzi genome was characterized by Southern blotting.

CA88, a nuclear repetitive DNA sequence identified in Schistosoma mansoni, aids in the genotyping of nine Schistosoma species of medical and veterinary importance

CA88 is the first long nuclear repetitive DNA sequence identified in the blood fluke, Schistosoma mansoni. The assembled S. mansoni sequence, which contains the CA88 repeat, has 8,887 nucleotides and at least three repeat units of approximately 360 bp. In addition, CA88 also possesses an internal CA microsatellite, identified as SmBr18. Both PCR and BLAST analysis have been used to analyse and confirm the CA88 sequence in other S. mansoni sequences in the public database. PCR-acquired nuclear repetitive DNA sequence profiles from nine Schistosoma species were used to classify this organism into four genotypes. Included among the nine species analysed were five sequences of both African and Asian lineages that are known to infect humans. Within these genotypes, three of them refer to recognised species groups. A panel of four microsatellite loci, including SmBr18 and three previously published loci, has been used to characterise the nine Schistosoma species. Each species has been identified and classified based on its CA88 DNA fingerprint profile. Furthermore, microsatellite sequences and intra-specific variation have also been observed within the nine Schistosoma species sequences. Taken together, these results support the use of these markers in studying the population dynamics of Schistosoma isolates from endemic areas and also provide new methods for investigating the relationships between different populations of parasites. In addition, these data also indicate that Schistosoma magrebowiei is not a sister taxon to Schistosoma mattheei, prompting a new designation to a basal clade.

Expression and localization of PPARs in the rat ovary during follicular development and the periovulatory period.

PPARs are a family of nuclear hormone receptors involved in various processes that could influence ovarian function. We investigated the cellular localization and expression of PPARs during follicular development in ovarian tissue collected from rats 0, 6, 12, 24, and 48 h post-PMSG. A second group of animals received human CG (hCG) 48 h post-PMSG. Their ovaries were removed 0, 4, 8, 12, and 24 h post-hCG to study the periovulatory period. mRNAs corresponding to the PPAR isotypes (alpha, delta, and gamma) were localized by in situ hybridization. Changes in the levels of mRNA for the PPARs were determined by ribonuclease protection assays. PPAR gamma mRNA was localized primarily to granulosa cells, and levels of expression did not change during follicular development. Four hours post-hCG, levels of mRNA for PPAR gamma decreased (P < 0.05) but not uniformly in all follicles. At 24 h post-hCG, levels of PPAR gamma mRNA were reduced 64%, but some follicles maintained high expression. In contrast, mRNAs for PPAR alpha and delta were located primarily in theca and stroma, and their levels did not change during the intervals studied. To investigate the physiologic significance of PPAR gamma in the ovary, granulosa cells from PMSG-primed rats were cultured for 48 h with prostaglandin J(2) (PGJ(2)) and ciglitazone, PPAR gamma activators. Both compounds increased progesterone and E2 secretion (P < 0.05). These data suggest that PPAR gamma is involved in follicular development, has a negative influence on the luteinization of granulosa cells, and/or regulates the periovulatory shift in steroid production. The more general and steady expression of PPARs alpha and delta indicate that they may play a role in basal ovarian function.

Localization of peroxisome proliferator-activated receptor alpha (PPARα) and N-acyl phosphatidylethanolamine phospholipase D (NAPE-PLD) in cells expressing the Ca(2+)-binding proteins calbindin, calretinin, and parvalbumin in the adult rat hippocampus.

The N-acylethanolamines (NAEs), oleoylethanolamide (OEA) and palmithylethanolamide (PEA) are known to be endogenous ligands of PPARα receptors, and their presence requires the activation of a specific phospholipase D (NAPE-PLD) associated with intracellular Ca(2+) fluxes. Thus, the identification of a specific population of NAPE-PLD/PPARα-containing neurons that express selective Ca(2+)-binding proteins (CaBPs) may provide a neuroanatomical basis to better understand the PPARα system in the brain. For this purpose, we used double-label immunofluorescence and confocal laser scanning microscopy for the characterization of the co-existence of NAPE-PLD/PPARα and the CaBPs calbindin D28k, calretinin and parvalbumin in the rat hippocampus. PPARα expression was specifically localized in the cell nucleus and, occasionally, in the cytoplasm of the principal cells (dentate granular and CA pyramidal cells) and some non-principal cells of the hippocampus. PPARα was expressed in the calbindin-containing cells of the granular cell layer of the dentate gyrus (DG) and the SP of CA1. These principal PPARα(+)/calbindin(+) cells were closely surrounded by NAPE-PLD(+) fiber varicosities. No pyramidal PPARα(+)/calbindin(+) cells were detected in CA3. Most cells containing parvalbumin expressed both NAPE-PLD and PPARα in the principal layers of the DG and CA1/3. A small number of cells containing PPARα and calretinin was found along the hippocampus. Scattered NAPE-PLD(+)/calretinin(+) cells were specifically detected in CA3. NAPE-PLD(+) puncta surrounded the calretinin(+) cells localized in the principal cells of the DG and CA1. The identification of the hippocampal subpopulations of NAPE-PLD/PPARα-containing neurons that express selective CaBPs should be considered when analyzing the role of NAEs/PPARα-signaling system in the regulation of hippocampal functions.

Genetic consequences of the ice ages on nurseries of the bat Myotis myotis: a mitochondrial and nuclear survey.

Analyses of mitochondrial DNA (mtDNA) control region polymorphism and of variation at 10 nuclear microsatellite loci were used to investigate the mechanisms and genetic consequences of postglacial expansion of Myotis myotis in Europe. Initial sampling consisted of 480 bats genotyped in 24 nursery colonies arranged along a transect of approximately 3000 km. The phylogeographical survey based on mtDNA sequences revealed the existence of major genetic subdivisions across this area, with several suture zones between haplogroups. Such zones of secondary contact were found in the Alps and Rhodopes, whereas other potential barriers to gene flow, like the Pyrenees, did not coincide with genetic discontinuities. Areas of population admixture increased locally the genetic diversity of colonies, which confounded the northward decrease in nucleotide diversity predicted using classical models of postglacial range expansion. However, when analyses were restricted to a subset of 15 nurseries originating from a single presumed glacial refugium, mtDNA polymorphism did indeed support a northwards decrease in diversity. Populations were also highly structured (PhiST = 0.384). Conversely, the same subset of colonies showed no significant latitudinal decrease in microsatellite diversity and much less population structure (FST = 0.010), but pairwise genetic differentiation at these nuclear markers was strongly correlated with increasing geographical distance. Together, this evidence suggests that alleles carried via male bats have maintained enough nuclear gene flow to counteract the effects of recurrent bottlenecks generally associated with recolonization processes. As females are highly philopatric, we argue that the maternally transmitted mtDNA marker better reflects the situation of past, historical gene flow, whereas current levels of gene flow are better reflected by microsatellite markers.

Mitotic functions for SNAP45, a subunit of the small nuclear RNA-activating protein complex SNAPc.

The small nuclear RNA-activating protein complex SNAP(c) is required for transcription of small nuclear RNA genes and binds to a proximal sequence element in their promoters. SNAP(c) contains five types of subunits stably associated with each other. Here we show that one of these polypeptides, SNAP45, also known as PTF delta, localizes to centrosomes during parts of mitosis, as well as to the spindle midzone during anaphase and the mid-body during telophase. Consistent with localization to these mitotic structures, both down- and up-regulation of SNAP45 lead to a G(2)/M arrest with cells displaying abnormal mitotic structures. In contrast, down-regulation of SNAP190, another SNAP(c) subunit, leads to an accumulation of cells with a G(0)/G(1) DNA content. These results are consistent with the proposal that SNAP45 plays two roles in the cell, one as a subunit of the transcription factor SNAP(c) and another as a factor required for proper mitotic progression.

Fusion of the human gene for the polyubiquitination co-effector UEV-1 with Kua, a newly identified gene

UEV proteins are enzymatically inactive variants of the E2 ubiquitin-conjugating enzymes that regulate noncanonical elongation of ubiquitin chains. In Saccharomyces cerevisiae, UEV is part of the RAD6-mediated error-free DNA repair pathway. In mammalian cells, UEV proteins can modulate c-FOS transcription and the G2-M transition of the cell cycle. Here we show that the UEV genes from phylogenetically distant organisms present a remarkable conservation in their exon–intron structure. We also show that the human UEV1 gene is fused with the previously unknown gene Kua. In Caenorhabditis elegans and Drosophila melanogaster, Kua and UEV are in separated loci, and are expressed as independent transcripts and proteins. In humans, Kua and UEV1 are adjacent genes, expressed either as separate transcripts encoding independent Kua and UEV1 proteins, or as a hybrid Kua–UEV transcript, encoding a two-domain protein. Kua proteins represent a novel class of conserved proteins with juxtamembrane histidine-rich motifs. Experiments with epitope-tagged proteins show that UEV1A is a nuclear protein, whereas both Kua and Kua–UEV localize to cytoplasmic structures, indicating that the Kua domain determines the cytoplasmic localization of Kua–UEV. Therefore, the addition of a Kua domain to UEV in the fused Kua–UEV protein confers new biological properties to this regulator of variant polyubiquitination.[Kua cDNAs isolated by RT-PCR and described in this paper have been deposited in the GenBank data library under accession nos. AF1155120 (H. sapiens) and AF152361 (D. melanogaster). Genomic clones containing UEV genes: S. cerevisiae, YGL087c (accession no. Z72609); S. pombe, c338 (accession no. AL023781); P. falciparum, MAL3P2 (accession no. AL034558); A. thaliana, F26F24 (accession no. AC005292); C. elegans, F39B2 (accession no. Z92834); D. melanogaster, AC014908; and H. sapiens, 1185N5 (accession no. AL034423). Accession numbers for Kua cDNAs in GenBank dbEST: M. musculus, AA7853; T. cruzi, AI612534. Other Kua-containing sequences: A. thaliana genomic clones F10M23 (accession no. AL035440), F19K23 (accession no. AC000375), and T20K9 (accession no. AC004786).

Forensic identification of urine samples: a comparison between nuclear and mitochondrial DNA markers.

Urine samples from 20 male volunteers of European Caucasian origin were stored at 4 degrees C over a 4-month period in order to compare the identification potential of nuclear DNA (nDNA) and mitochondrial DNA (mtDNA) markers. The amount of nDNA recovered from urines dramatically declined over time. Consequently, nDNA likelihood ratios (LRs) greater than 1,000 were obtained for 100, 70 and 55% of the urines analysed after 6, 60 and 120 days, respectively. For the mtDNA, HVI and HVII sequences were obtained for all samples tested, whatever the period considered. Nevertheless, the highest mtDNA LR of 435 was relatively low compared to its nDNA equivalent. Indeed, LRs obtained with only three nDNA loci could easily exceed this value and are quite easier to obtain. Overall, the joint use of nDNA and mtDNA markers enabled the 20 urine samples to be identified, even after the 4-month period.

Echo-time independent signal modulations for strongly coupled systems in triple echo localization schemes: an extension of S-PRESS editing.

The double spin-echo point resolved spectroscopy sequence (PRESS) is a widely used method and standard in clinical MR spectroscopy. Existence of important J-modulations at constant echo times, depending on the temporal delays between the rf-pulses, have been demonstrated recently for strongly coupled spin systems and were exploited for difference editing, removing singlets from the spectrum (strong-coupling PRESS, S-PRESS). A drawback of this method for in vivo applications is that large signal modulations needed for difference editing occur only at relatively long echo times. In this work we demonstrate that, by simply adding a third refocusing pulse (3S-PRESS), difference editing becomes possible at substantially shorter echo times while, as applied to citrate, more favorable lineshapes can be obtained. For the example of an AB system an analytical description of the MR signal, obtained with this triple refocusing sequence (3S-PRESS), is provided.

Immunoscintigraphic localization of inflammatory lesions: concept, radiolabelling and in vitro testing of a granulocyte specific antibody.

Current nuclear medicine techniques for the localization of inflammatory processes are based on injection of 111In labelled autologous granulocytes which need to be isolated and radiolabelled in vitro before reinjection. A new technique is presented here that obviates the need for cell isolation by the direct intravenous injection of a granulocyte specific 123I labelled monoclonal antibody. In this publication the basic parameters of the antibody granulocyte interaction are described. Antibody binding does not inhibit vital functions of the granulocytes, such as chemotaxis and superoxide generation. Scatchard analysis of binding data reveals an apparent affinity of the antibody for granulocytes of 6.8 X 10(9) l/mol and approximately 7.1 X 10(4) binding sites per cell. Due to the high specificity of the antibody, the only expected interference is from CEA producing tumors.