Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

A glycine-arginine domain in control of the human MRE11 DNA repair protein.

Human MRE11 is a key enzyme in DNA double-strand break repair and genome stability. Human MRE11 bears a glycine-arginine-rich (GAR) motif that is conserved among multicellular eukaryotic species. We investigated how this motif influences MRE11 function. Human MRE11 alone or a complex of MRE11, RAD50, and NBS1 (MRN) was methylated in insect cells, suggesting that this modification is conserved during evolution. We demonstrate that PRMT1 interacts with MRE11 but not with the MRN complex, suggesting that MRE11 arginine methylation occurs prior to the binding of NBS1 and RAD50. Moreover, the first six methylated arginines are essential for the regulation of MRE11 DNA binding and nuclease activity. The inhibition of arginine methylation leads to a reduction in MRE11 and RAD51 focus formation on a unique double-strand break in vivo. Furthermore, the MRE11-methylated GAR domain is sufficient for its targeting to DNA damage foci and colocalization with gamma-H2AX. These studies highlight an important role for the GAR domain in regulating MRE11 function at the biochemical and cellular levels during DNA double-strand break repair.

Results of a collaborative study of the EDNAP group regarding mitochondrial DNA heteroplasmy and segregation in hair shafts.

A collaborative exercise was carried out by the European DNA Profiling Group (EDNAP) in order to evaluate the distribution of mitochondrial DNA (mtDNA) heteroplasmy amongst the hairs of an individual who displays point heteroplasmy in blood and buccal cells. A second aim of the exercise was to study reproducibility of mtDNA sequencing of hairs between laboratories using differing chemistries, further to the first mtDNA reproducibility study carried out by the EDNAP group. Laboratories were asked to type 2 sections from each of 10 hairs, such that each hair was typed by at least two laboratories. Ten laboratories participated in the study, and a total of 55 hairs were typed. The results showed that the C/T point heteroplasmy observed in blood and buccal cells at position 16234 segregated differentially between hairs, such that some hairs showed only C, others only T and the remainder, C/T heteroplasmy at varying ratios. Additionally, differential segregation of heteroplasmic variants was confirmed in independent extracts at positions 16093 and the poly(C) tract at 302-309, whilst a complete A-G transition was confirmed at position 16129 in one hair. Heteroplasmy was observed at position 16195 on both strands of a single extract from one hair segment, but was not observed in the extracts from any other segment of the same hair. Similarly, heteroplasmy at position 16304 was observed on both strands of a single extract from one hair. Additional variants at positions 73, 249 and the HVII poly(C) region were reported by one laboratory; as these were not confirmed in independent extracts, the possibility of contamination cannot be excluded. Additionally, the electrophoresis and detection equipment used by this laboratory was different to those of the other laboratories, and the discrepancies at position 249 and the HVII poly(C) region appear to be due to reading errors that may be associated with this technology. The results, and their implications for forensic mtDNA typing, are discussed in the light of the biology of hair formation.

Use of DNA profiles for investigation using a simulated national DNA database: Part I. Partial SGM Plus profiles.

In traditional criminal investigation, uncertainties are often dealt with using a combination of common sense, practical considerations and experience, but rarely with tailored statistical models. For example, in some countries, in order to search for a given profile in the national DNA database, it must have allelic information for six or more of the ten SGM Plus loci for a simple trace. If the profile does not have this amount of information then it cannot be searched in the national DNA database (NDNAD). This requirement (of a result at six or more loci) is not based on a statistical approach, but rather on the feeling that six or more would be sufficient. A statistical approach, however, could be more rigorous and objective and would take into consideration factors such as the probability of adventitious matches relative to the actual database size and/or investigator's requirements in a sensible way. Therefore, this research was undertaken to establish scientific foundations pertaining to the use of partial SGM Plus loci profiles (or similar) for investigation.

A New Large-DNA-Fragment Delivery System Based on Integrase Activity from an Integrative and Conjugative Element.

During the past few decades, numerous plasmid vectors have been developed for cloning, gene expression analysis, and genetic engineering. Cloning procedures typically rely on PCR amplification, DNA fragment restriction digestion, recovery, and ligation, but increasingly, procedures are being developed to assemble large synthetic DNAs. In this study, we developed a new gene delivery system using the integrase activity of an integrative and conjugative element (ICE). The advantage of the integrase-based delivery is that it can stably introduce a large DNA fragment (at least 75 kb) into one or more specific sites (the gene for glycine-accepting tRNA) on a target chromosome. Integrase recombination activity in Escherichia coli is kept low by using a synthetic hybrid promoter, which, however, is unleashed in the final target host, forcing the integration of the construct. Upon integration, the system is again silenced. Two variants with different genetic features were produced, one in the form of a cloning vector in E. coli and the other as a mini-transposable element by which large DNA constructs assembled in E. coli can be tagged with the integrase gene. We confirmed that the system could successfully introduce cosmid and bacterial artificial chromosome (BAC) DNAs from E. coli into the chromosome of Pseudomonas putida in a site-specific manner. The integrase delivery system works in concert with existing vector systems and could thus be a powerful tool for synthetic constructions of new metabolic pathways in a variety of host bacteria.

Are there two species of Pygmy Shrews (Sorex)? Revisiting the question using DNA sequence data

Pygmy Shrews in North America have variously been considered to be one species (Sorex hoyi) or two species (S. hoyi and S. thompsoni). Currently, only S. hoyi is recognized. In this study, we examine mitochondrial DNA sequence data for the cytochrome b gene to evaluate the level of differentiation and phylogeographic relationships among eleven samples of Pygmy Shrews from across Canada. Pygmy Shrews from eastern Canada (i.e., Ontario, Quebec, New Brunswick, Nova Scotia, and Prince Edward Island) are distinct from Pygmy Shrews from western Canada (Alberta, Yukon) and Alaska. The average level of sequence divergence between these clades (3.3%) falls within the range of values for other recognized pairs of sister species of shrews. A molecular clock based on third position transversion substitutions suggests that these two lineages diverged between 0.44 and 1.67 million years ago. These molecular phylogenetic data. combined with a reinterpretation of previously published morphological data, are suggestive of separate species status for S. hoyi and S. thompsoni as has been previously argued by others. Further analysis of specimens from geographically intermediate areas (e.g., Manitoba. northern Ontario) is required to determine if there is secondary contact and/or introgression between these two putative species.

Tightening of DNA knots by supercoiling facilitates their unknotting by type II DNA topoisomerases.

Using numerical simulations, we compare properties of knotted DNA molecules that are either torsionally relaxed or supercoiled. We observe that DNA supercoiling tightens knotted portions of DNA molecules and accentuates the difference in curvature between knotted and unknotted regions. The increased curvature of knotted regions is expected to make them preferential substrates of type IIA topoisomerases because various earlier experiments have concluded that type IIA DNA topoisomerases preferentially interact with highly curved DNA regions. The supercoiling-induced tightening of DNA knots observed here shows that torsional tension in DNA may serve to expose DNA knots to the unknotting action of type IIA topoisomerases, and thus explains how these topoisomerases could maintain a low knotting equilibrium in vivo, even for long DNA molecules.

Visualisation of human rad52 protein and its complexes with hRad51 and DNA.

The human Rad52 protein stimulates joint molecule formation by hRad51, a homologue of Escherichia coli RecA protein. Electron microscopic analysis of hRad52 shows that it self-associates to form ring structures with a diameter of approximately 10 nm. Each ring contains a hole at its centre. hRad52 binds to single and double-stranded DNA. In the ssDNA-hRad52 complexes, hRad52 was distributed along the length of the DNA, which exhibited a characteristic "beads on a string" appearance. At higher concentrations of hRad52, "super-rings" (approximately 30 nm) were observed and the ssDNA was collapsed upon itself. In contrast, in dsDNA-hRad52 complexes, some regions of the DNA remained protein-free while others, containing hRad52, interacted to form large protein-DNA networks. Saturating concentrations of hRad51 displaced hRad52 from ssDNA, whereas dsDNA-Rad52 complexes (networks) were more resistant to hRad51 invasion and nucleoprotein filament formation. When Rad52-Rad51-DNA complexes were probed with gold-conjugated hRad52 antibodies, the presence of globular hRad52 structures within the Rad51 nucleoprotein filament was observed. These data provide the first direct visualisation of protein-DNA complexes formed by the human Rad51 and Rad52 recombination/repair proteins.

The twist, writhe and overall shape of supercoiled DNA change during counterion-induced transition from a loosely to a tightly interwound superhelix. Possible implications for DNA structure in vivo.

A cryo-electron microscopy study of supercoiled DNA molecules freely suspended in cryo-vitrified buffer was combined with Monte Carlo simulations and gel electrophoretic analysis to investigate the role of intersegmental electrostatic repulsion in determining the shape of supercoiled DNA molecules. It is demonstrated here that a decrease of DNA-DNA repulsion by increasing concentrations of counterions causes a higher fraction of the linking number deficit to be partitioned into writhe. When counterions reach concentrations likely to be present under in vivo conditions, naturally supercoiled plasmids adopt a tightly interwound conformation. In these tightly supercoiled DNA molecules the opposing segments of interwound superhelix seem to directly contact each other. This form of supercoiling, where two DNA helices interact laterally, may represent an important functional state of DNA. In the particular case of supercoiled minicircles (178 bp) the delta Lk = -2 topoisomers undergo a sharp structural transition from almost planar circles in low salt buffers to strongly writhed "figure-eight" conformations in buffers containing neutralizing concentrations of counterions. Possible implications of this observed structural transition in DNA are discussed.

Structure of a multisubunit complex that promotes DNA branch migration.

The RuvA and RuvB proteins of Escherichia coli, which are induced in response to DNA damage, are important in the formation of heteroduplex DNA during genetic recombination and related recombinational repair processes. In vitro studies show that RuvA binds Holiday junctions and acts as a specificity factor that targets the RuvB ATPase, a hexameric ring protein, to the junction. Together, RuvA and RuvB promote branch migration, an ATP-dependent reaction that increases the length of the heteroduplex DNA. Electron microscopic visualization of RuvAB now provides a new insight into the mechanism of this process. We observe the formation of a tripartite protein complex in which RuvA binds the crossover and is sandwiched between two hexameric rings of RuvB. The Holliday junction within this complex adopts a square-planar structure. We propose a molecular model for branch migration, a unique feature of which is the role played by the two oppositely oriented RuvB ring motors.

Species identification in mammals from mixed biological samples based on mitochondrial DNA control region length polymorphism.

The ability to identify the species origin of an unknown biological sample is relevant in the fields of human and wildlife forensics. However, the detection of several species mixed in the same sample still remains a challenge. We developed and tested a new approach for mammal DNA identification in mixtures of two or three species, based on the analysis of mitochondrial DNA control region interspecific length polymorphism followed by direct sequencing. Contrary to other published methods dealing with species mixtures, our protocol requires a single universal primer pair and is not based on a pre-defined panel of species. Amplicons can be separated either on agarose gels or using CE. The advantages and limitations of the assay are discussed under different conditions, such as variable template concentration, amplicon sizes and size difference among the amplicons present in the mixture. For the first time, this protocol provides a simple, reliable and flexible method for simultaneous identification of multiple mammalian species from mixtures, without any prior knowledge of the species involved.

DNA supercoiling and its role in DNA decatenation and unknotting.

Chromosomal and plasmid DNA molecules in bacterial cells are maintained under torsional tension and are therefore supercoiled. With the exception of extreme thermophiles, supercoiling has a negative sign, which means that the torsional tension diminishes the DNA helicity and facilitates strand separation. In consequence, negative supercoiling aids such processes as DNA replication or transcription that require global- or local-strand separation. In extreme thermophiles, DNA is positively supercoiled which protects it from thermal denaturation. While the role of DNA supercoiling connected to the control of DNA stability, is thoroughly researched and subject of many reviews, a less known role of DNA supercoiling emerges and consists of aiding DNA topoisomerases in DNA decatenation and unknotting. Although DNA catenanes are natural intermediates in the process of DNA replication of circular DNA molecules, it is necessary that they become very efficiently decatenated, as otherwise the segregation of freshly replicated DNA molecules would be blocked. DNA knots arise as by-products of topoisomerase-mediated intramolecular passages that are needed to facilitate general DNA metabolism, including DNA replication, transcription or recombination. The formed knots are, however, very harmful for cells if not removed efficiently. Here, we overview the role of DNA supercoiling in DNA unknotting and decatenation.

CTCF-T, a paralogue of CTCF, directs sequence specific DNA methylation of the imprinting control region of Igf2/H19

SUMMARY Genomic imprinting is an epigenetic mechanism of transcriptional regulation that ensures restriction of expression of a subset of mammalian genes to a single parental allele. The best studied example of imprinted gene regulation is the Igf2/H19 locus, which is also the most commonly altered by loss of imprinting (LOT) in cancer. LOT is associated with numerous hereditary diseases and several childhood, and adult cancers. Differential expression of reciprocal H19 and 1gf2 alleles in somatic cells depends on the methylation status of the imprinting control region (ICR) which regulates binding of CTCF, an ubiquitously expressed 11-zinc finger protein that binds specifically to non-methylated maternal ICR and thereby attenuates expression of Igf2, while it does not bind to methylated paternal ICR, which enables Igf2 expression. Initial ICR methylation occurs during gametogenesis by an as yet unknown mechanism. The accepted hypothesis is that the event of differential maternal and paternal DNA methylation depends on germ-line specific proteins. Our Laboratory identified a novel 11-zinc-finger protein CTCF-T (also known as CTCFL and BORIS) that is uniquely expressed in the male germ-line and is highly homologous within its zinc-finger region with CTCF. The amino-acid sequences flanking the zinc-finger regions of CTCF and CTCF-T have widely diverged, suggesting that though they could bind to the same DNA targets (ICRs) they are likely to have different functions. Interestingly, expression of CTCF-T and CTCF is mutually exclusive; CTCF-T-positive (CTCF-negative) cells occur in the stage of spermatogenesis that coincides with epigenetic reprogramming, including de novo DNA methylation. In our study we demonstrate the role that CTCF-T plays in genomic imprinting. Here we show that CTCF-T binds in vivo to the ICRs of Igf2/H19 and Dlk/Gt12 imprinted genes. In addition, we identified two novel proteins interacting with CTCF-T: a protein arginine methyltransferase PRMT7 and an arginine-rich histone H2A variant that we named trH2A. These interactions were confirmed and show that the two proteins interact with the amino-teiminal region of CTCF-T. Additionally, we show interaction of the amino- terminal region of CTCF-T with histones H1, H2A and H3. These results suggest that CTCF-T is a sequence-specific DNA (ICR) binding protein that associates with histones and recruits PRMT7. Interestingly, PRMT7 has a histone-methyltransferase activity. It has been shown that histone methylation can mark chromatin regions thereby directing DNA-methylation; thus, our hypothesis is that the CTCF-T protein-scaffold directs PRMT7 to methylate histone(s) assembled on ICRs, which marks chromatin for the recruitment of the de novo DNA methyltransferases to methylate DNA. To test this hypothesis, we developed an in vivo DNA-methylation assay using Xenopus laevis' oocytes, where H19 ICR and different expression cDNAs, including CTCF-T, PRMT7 and the de novo DNA methyltransferases (Dnmt3a, Dnmt3b and Dnmt3L) are microinjected into the nucleus. The methylation status of CpGs within the H19 ICR was analysed 48 or 72 hours after injection. Here we demonstrate that CpGs in the ICR are methylated in the presence of both CTCF-T and PRMT7, while control oocytes injected only with ICR did not show any methylation. Additionally, we showed for the first time that Dnmt3L is crucial for the establishment of the imprinting marks on H19 ICR. Moreover, we confirmed that Dnmt3a and Dnmt3b activities are complementary. Our data indicate that all three Dnmt3s are important for efficient de novo DNA methylation. In conclusion, we propose a mechanism for the establishment of de novo imprinting marks during spermatogenesis: the CTCF-T/PRMT7 protein complex directs histone methylation leading to sequence-specific de novo DNA methylation of H19 ICR. RESUME L'empreinte génomique parentale est un mécanisme épigénétique de régulation transcriptionelle qui se traduit par une expression différentielle des deux allèles de certains gènes, en fonction de leur origine parentale. L'exemple le mieux caractérisé de gènes soumis à l'empreinte génomique parentale est le locus Igf2/H19, qui est aussi le plus fréquemment altéré par relaxation d'empreinte (en anglais: loss of imprinting, LOI) dans les cancers. Cette relaxation d'empreinte est aussi associée à de nombreuses maladies héréditaires, ainsi qu'à de nombreux cancers chez l'enfant et l'adulte. Dans les cellules somatiques, les différences d'expression des allèles réciproques H19 et Ig12 est sous le contrôle d'une région ICR (Imprinting Control Region). La méthylation de cette région ICR régule l'ancrage de la protéine à douze doigts de zinc CTCF, qui se lie spécifiquement à l'ICR maternel non-méthylé, atténuant ainsi l'expression de Igf2, alors qu'elle ne s'ancre pas à l'ICR paternel méthyle. Le mécanisme qui accompagne la méthylation initiale de la région ICR durant la gamétogenèse n'a toujours pas été élucidé. L'hypothèse actuelle propose que la différence de méthylation entre l'ADN maternel et paternel résulte de l'expression de protéines propres aux zones germinales. Notre laboratoire a récemment identifié une nouvelle protéine à douze doigts de zinc, CTCF-T (aussi dénommée CTCFL et BORRIS), qui est exprimée uniquement dans les cellules germinales mâles, dont la partie à douze doigts de zinc est fortement homologue à la protéine CTCF. La séquence d'acides aminés de part et d'autre de cette région est quant à elle très divergente, ce qui implique que CTCF-T se lie sans doute au même ADN cible que CTCF, mais possède des fonctions différentes. De plus, l'expression de CTCF-T et de CTCF s'oppose mutuellement; l'expression de la protéine CTCF-T (cellules CTCF-T positives, CTCF negatives) qui a lieu pendant la spermatogenèse coïncide avec la reprogrammation épigénétique, notamment la méthylation de novo de l'ADN. La présente étude démontre le rôle essentiel joué par la protéine CTCF-T dans l'acquisition de l'empreinte génomique parentale. Nous montrons ici que CTCF-T s'associe in vivo avec les régions ICR des loci Igf2/H19 et Dlk/Gt12. Nous avons également identifié deux nouvelles protéines qui interagissent avec CTCF-T : une protéine arginine méthyl transférase PRMT7, et un variant de l'histone H2A, riche en arginine, que nous avons dénommé trH2A. Ces interactions ont été analysées plus en détail, et confinnent que ces deux protéines s'associent avec la région N-terminale de CTCF-T. Aussi, nous présentons une interaction de la région N-terminale de CTCF-T avec les histones H1, H2, et H3. Ces résultats suggèrent que CTCF-T est une protéine qui se lie spécifiquement aux régions ICR, qui s'associe avec différents histones et qui recrute PRMT7. PRMT7 possède une activité méthyl-tansférase envers les histones. Il a été montré que la méthylation des histones marque certains endroits de la chromatine, dirigeant ainsi la méthylation de l'ADN. Notre hypothèse est donc la suivante : la protéine CTCF-T sert de base qui dirige la méthylation des histones par PRMT7 dans les régions ICR, ce qui contribue à marquer la chromatine pour le recrutement de nouvelles méthyl transférases pour méthyler l'ADN. Afin de valider cette hypothèse, nous avons développé un système de méthylation de l'ADN in vivo, dans des oeufs de Xenopus laevis, dans le noyau desquels nous avons mico-injecté la région ICR du locus H19, ainsi que différents vecteurs d'expression pour CTCF-T, PRMT7, et les de novo méthyl transférases (Dnmt3a, Dnmt3b et Dnmt3L). Les CpGs méthyles de la région ICR du locus H19 ont été analysé 48 et 72 heures après l'injection. Cette technique nous a permis de démontrer que les CpGs de la région ICR sont méthyles en présence de CTCF-T et de PRMT7, tandis que les contrôles injectés seulement avec la région ICR ne présentent aucun signe de méthylation. De plus, nous démontrons pour la première fois que la protéine méthyl transférase Dnmt3L est déterminant pour l'établissement de l'empreinte génomique parentale au niveau de la région ICR du locus H19. Aussi, nous confirmons que les activités méthyl transférases de Dnmt3a et Dnmt3b sont complémentaires. Nos données indiquent que les trois protéines Dnmt3 sont impliquées dans la méthylation de l'ADN. En conclusion, nous proposons un mécanisme responsable de la mise en place de nouvelles empreintes génomiques pendant la spermatogenèse : le complexe protéique CTCF-T/PRMT7 dirige la méthylation des histones aboutissant à la méthylation de novo de l'ADN au locus H19.

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.