Changes in cis- and trans-regulatory elements are among the prime sources of genetic and phenotypical variation at species level. The introduction of cis- and trans- regulatory variation has played important roles in driving diversity, phenotypical differentiation, and evolution of humans. Therefore, variation that occurs on cis- and trans- regulatory elements becomes imperative to better understanding of human genetic diversity and its evolution. In this research, around 3360 gene regulatory factors (GRF) from the human genome were catalogued. This catalog includes genes that code for proteins that perform gene regulatory activities such DNA-depending transcription, RNA polymerase II transcription cofactor and co-repressor activity, chromatin binding and remodeling, among other 218 regulatory functions. This GRF catalog allowed us to initially explore how some GRF genes have evolved in humans, archaic humans (Neandertal and Denisovan) and non-human primate species. We discussed the likely phenotypical and medical effects that evolutionary changes in GRF genes may have introduced into the human genome; for instance, traits associated to speech and language capabilities, genomic recombination hotspots, diseases, among others. By using genome-wide datasets, we additionally looked for GRFs likely to be candidates for positive selection in three human populations: Utah Residents with Northern and Western Ancestry (CEU), Han Chinese in Beijing (CHB), and Yoruba in Ibadan (YRI). As result, we produced a set of candidates that gathers genes that may have contributed in shaping the phenotypical diversity currently observed in these populations; for instance, by introducing regulatory diversity at population-specific level. We additionally identified six GRF classes enriched for genes located in regions that are likely candidates for positive selection at population specific level. We found that out of the 41 DNA-binding GRF classes classified so far, six groups exhibited enrichment for genes located on regions that may have been under positive selection: C2H2 zinc finger, KRAB-ZNF zinc finger, Homeo domain, Tryptophan cluster, Fork head/winged helix and, and High-mobility HMG domain. We additionally identified three KRAB-ZNF gene clusters, in the chromosomes one, three, and 16, for the Asian population that exhibit regions with extended haplotype homozygosity EHH (larger than 100 kb). This EHH suggests that these regions have undergone positive selection in CHB population. Finally, considering that a representative fraction of the phenotypic diversity observed between humans and its closely related species are likely explained by changes in cis-regulatory elements (CREs), we investigated putative binding sites for the transcription factor GABPa. Using ChIP-Seq data generated from a human cell line (HEK293T), 11,619 putative GABPa CREs were found, Out of which 224 are putative human-specific. To experimentally validate the transcriptional activity of these human-specific CREs, reporter gene essays and knock-down experiments were performed. Our results supported the functionality of these human-specific GABPa CREs and suggest that at least 1,215 genes are primary targets of GABPa. Finally, further analyses depict scenarios that put together transcriptional regulation by GABPa and the evolution of particular human traits; for instance, cognitive abilities, breast morphology, lipids and glucose metabolic pathways, among others.

Changes in cis- and trans-regulatory elements are among the prime sources of genetic and phenotypical variation at species level. The introduction of cis- and trans- regulatory variation has played important roles in driving diversity, phenotypical differentiation, and evolution of humans. Therefore, variation that occurs on cis- and trans- regulatory elements becomes imperative to better understanding of human genetic diversity and its evolution. In this research, around 3360 gene regulatory factors (GRF) from the human genome were catalogued. This catalog includes genes that code for proteins that perform gene regulatory activities such DNA-depending transcription, RNA polymerase II transcription cofactor and co-repressor activity, chromatin binding and remodeling, among other 218 regulatory functions. This GRF catalog allowed us to initially explore how some GRF genes have evolved in humans, archaic humans (Neandertal and Denisovan) and non-human primate species. We discussed the likely phenotypical and medical effects that evolutionary changes in GRF genes may have introduced into the human genome; for instance, traits associated to speech and language capabilities, genomic recombination hotspots, diseases, among others. By using genome-wide datasets, we additionally looked for GRFs likely to be candidates for positive selection in three human populations: Utah Residents with Northern and Western Ancestry (CEU), Han Chinese in Beijing (CHB), and Yoruba in Ibadan (YRI). As result, we produced a set of candidates that gathers genes that may have contributed in shaping the phenotypical diversity currently observed in these populations; for instance, by introducing regulatory diversity at population-specific level. We additionally identified six GRF classes enriched for genes located in regions that are likely candidates for positive selection at population specific level. We found that out of the 41 DNA-binding GRF classes classified so far, six groups exhibited enrichment for genes located on regions that may have been under positive selection: C2H2 zinc finger, KRAB-ZNF zinc finger, Homeo domain, Tryptophan cluster, Fork head/winged helix and, and High-mobility HMG domain. We additionally identified three KRAB-ZNF gene clusters, in the chromosomes one, three, and 16, for the Asian population that exhibit regions with extended haplotype homozygosity EHH (larger than 100 kb). This EHH suggests that these regions have undergone positive selection in CHB population. Finally, considering that a representative fraction of the phenotypic diversity observed between humans and its closely related species are likely explained by changes in cis-regulatory elements (CREs), we investigated putative binding sites for the transcription factor GABPa. Using ChIP-Seq data generated from a human cell line (HEK293T), 11,619 putative GABPa CREs were found, Out of which 224 are putative human-specific. To experimentally validate the transcriptional activity of these human-specific CREs, reporter gene essays and knock-down experiments were performed. Our results supported the functionality of these human-specific GABPa CREs and suggest that at least 1,215 genes are primary targets of GABPa. Finally, further analyses depict scenarios that put together transcriptional regulation by GABPa and the evolution of particular human traits; for instance, cognitive abilities, breast morphology, lipids and glucose metabolic pathways, among others.

Cis-regulatorische Elemente wie Promotoren und Enhancer, die die Regulation der Transkription von Genen steuern, befinden sich in Regionen des dekondensierten Chromatins. DNase-seq und ATAC-seq sind weit verbreitete Verfahren, um solche offenen Chromatinregionen genomweit zu untersuchen. Die einzel-Nukleotid-Auflösung von DNase-seq wurde des Weiteren genutzt, um Transkriptionsfaktor-Bindungsstellen (TFBS) in regulatorischen Regionen durch TF-Footprinting zu bestimmen. Kürzlich durchgeführte Studien haben jedoch gezeigt, dass DNase I einen Sequenzbias aufweist, welcher nachteilige Auswirkungen auf die Footprinting-Effizienz hat. Auch wurden das Footprinting und die Auswirkungen des Sequenzbias auf ATAC-seq noch nicht umfassend untersucht. In dieser Arbeit nehme ich einen systematischen Vergleich der beiden Methoden vor und zeige, dass die beiden Methoden unterschiedliche Sequenzbiases haben und korrigiere diese protokollspezifischen Biases beim Footprinting. Der Einfluss von Bias-Korrekturen der Footprinting Ergebnisse ist für DNase-seq größer als für ATAC-seq, und Footprinting mit DNase-seq führt zu besseren Ergebnissen in unserer Datensätze. Trotz dieser Unterschiede zeige ich, dass die Integration replizierter Experimente die Ableitung von qualitativ hochwertigen Footprints ermöglicht, wobei die beiden Techniken weitgehend übereinstimmen. Diese Techniken werden ferner eingesetzt, um die cis-regulatorischen Elemente zu charakterisieren, die die Embryogenese der Fruchtfliege Drosophila melanogaster bestimmen. Durch die Verwendung von Embryonen die sich im richtigen Entwicklungsstadium befinden, sowie gewebespezifischer Kernsortierung mit offenem Chromatin-Profiling können zeitlich und gewebespezifisch aufgelöste vermeintliche cis-regulatorische Elemente definiert werden. Zusammengenommen demonstrieren diese Analysen die Fähigkeit der offenen Chromatin-Profilierung und der Computeranalyse zur Aufklärung der Mechanismen der Genregulation. / Cis-regulatory elements such as promoters and enhancers, that govern transcriptional gene regulation, reside in regions of open chromatin. DNase-seq and ATAC-seq are broadly used methods to assay open chromatin regions genome-wide. The single nucleotide resolution of DNase-seq has been further exploited to infer transcription factor binding sites (TFBS) in regulatory regions through TF footprinting. However, recent studies have demonstrated the sequence bias of DNase I and its adverse effects on footprinting efficiency. Furthermore, footprinting and the impact of sequence bias have not been extensively studied for ATAC-seq. In this thesis, I undertake a systematic comparison of the two methods and demonstrate that the two methods have distinct sequence biases and correct for these protocol-specific biases when performing footprinting. The impact of bias correction on footprinting performance is greater for DNase-seq than for ATAC-seq, and footprinting with DNase-seq leads to better performance in our datasets. Despite these differences, I show that integrating replicate experiments allows the inference of high-quality footprints, with substantial agreement between the two techniques. These techniques are further employed to characterize the cis-regulatory elements governing the embryogenesis of a complex organism, the fruit fly Drosophila melanogaster. Combining tight staging of embryos and tissue-specific nuclear sorting with open chromatin profiling, enables the definition of temporally and tissue-specifically resolved putative cis-regulatory elements. Taken together, these analyses demonstrate the power of open chromatin profiling and computational analysis in elucidating the mechanisms of transcriptional gene regulation.