RESUMO
Gene Regulatory Networks (GRNs) allow the study of regulation of gene expression of whole genomes. Among the most relevant advantages of using networks to depict this key process, there is the visual representation of large amounts of information and the application of graph theory to generate new knowledge. Nonetheless, despite the many uses of GRNs, it is still difficult and expensive to assign Transcription Factors (TFs) to the regulation of specific genes. ChIP-Seq allows the determination of TF Binding Sites (TFBSs) over whole genomes, but it is still an expensive technique that can only be applied one TF at a time and requires replicates to reduce its noise. Once TFBSs are determined, the assignment of each TF and its binding sites to the regulation of specific genes is not trivial, and it is often performed by carrying out site-specific experiments that are unfeasible to perform in all possible binding sites. Here, we addressed these relevant issues with a two-step methodology using Drosophila melanogaster as a case study. First, our protocol starts by gathering all transcription factor binding sites (TFBSs) determined with ChIP-Seq experiments available at ENCODE and FlyBase. Then each TFBS is used to assign TFs to the regulation of likely target genes based on the TFBS proximity to the transcription start site of all genes. In the final step, to try to select the most likely regulatory TF from those previously assigned to each gene, we employ GENIE3, a random forest-based method, and more than 9,000 RNA-seq experiments from D. melanogaster. Following, we employed known TF protein-protein interactions to estimate the feasibility of regulatory events in our filtered networks. Finally, we show how known interactions between co-regulatory TFs of each gene increase after the second step of our approach, and thus, the consistency of the TF-gene assignment. Also, we employed our methodology to create a network centered on the Drosophila melanogaster gene Hr96 to demonstrate the role of this transcription factor on mitochondrial gene regulation.
RESUMO
Neurofibromatosis type 1 (NF1) patients are more likely to have vitamin D deficiency when compared to the general population. This study aimed to determine the levels of 25-OH-vitamin D [25(OH)D] in individuals with NF1 and disease-unaffected controls and analyze FokI and BsmI VDR gene polymorphisms in a case and in a control group. Vitamin D levels were compared between a group of 45 NF1 patients from Southern Brazil and 45 healthy controls matched by sex, skin type, and age. Genotypic and allelic frequencies of VDR gene polymorphisms were obtained from the same NF1 patients and 150 healthy controls. 25(OH)D deficiency or insufficiency was not more frequent in NF1 patients than in controls (p = 0.074). We also did not observe an association between FokI and BsmI VDR gene polymorphisms and vitamin D levels in NF1 patients, suggesting that their deficient or insufficient biochemical phenotypes are not associated with these genetic variants. The differences between the groups in genotypic and allelic frequencies for FokI and BsmI VDR gene polymorphisms were small and did not reach statistical significance. These polymorphisms are in partial linkage disequilibrium and the haplotype frequencies also did not differ in a significant way between the two groups (p = 0.613).
