Benefiting from the intrinsic role of epigenetics to predict patterns of CTCF binding.

Motivation: One of the most relevant mechanisms involved in the determination of chromatin structure is the formation of structural loops that are also related with the conservation of chromatin states. Many of these loops are stabilized by CCCTC-binding factor (CTCF) proteins at their base. Despite the relevance of chromatin structure and the key role of CTCF, the role of the epigenetic factors that are involved in the regulation of CTCF binding, and thus, in the formation of structural loops in the chromatin, is not thoroughly understood. Results: Here we describe a CTCF binding predictor based on Random Forest that employs different epigenetic data and genomic features. Importantly, given the ability of Random Forests to determine the relevance of features for the prediction, our approach also shows how the different types of descriptors impact the binding of CTCF, confirming previous knowledge on the relevance of chromatin accessibility and DNA methylation, but demonstrating the effect of epigenetic modifications on the activity of CTCF. We compared our approach against other predictors and found improved performance in terms of areas under PR and ROC curves (PRAUC-ROCAUC), outperforming current state-of-the-art methods.

Validation of an NGS Panel Designed for Detection of Actionable Mutations in Tumors Common in Latin America.

Next-generation sequencing (NGS) is progressively being used in clinical practice. However, several barriers preclude using this technology for precision oncology in most Latin American countries. To overcome some of these barriers, we have designed a 25-gene panel that contains predictive biomarkers for most current and near-future available therapies in Chile and Latin America. Library preparation was optimized to account for low DNA integrity observed in formalin-fixed paraffin-embedded tissue. The workflow includes an automated bioinformatic pipeline that accounts for the underrepresentation of Latin Americans in genome databases. The panel detected small insertions, deletions, and single nucleotide variants down to allelic frequencies of 0.05 with high sensitivity, specificity, and reproducibility. The workflow was validated in 272 clinical samples from several solid tumor types, including gallbladder (GBC). More than 50 biomarkers were detected in these samples, mainly in BRCA1/2, KRAS, and PIK3CA genes. In GBC, biomarkers for PARP, EGFR, PIK3CA, mTOR, and Hedgehog signaling inhibitors were found. Thus, this small NGS panel is an accurate and sensitive method that may constitute a more cost-efficient alternative to multiple non-NGS assays and costly, large NGS panels. This kind of streamlined assay with automated bioinformatics analysis may facilitate the implementation of precision medicine in Latin America.

Automated generation of context-specific gene regulatory networks with a weighted approach in Drosophila melanogaster.

The regulation of gene expression is a key factor in the development and maintenance of life in all organisms. Even so, little is known at whole genome scale for most genes and contexts. We propose a method, Tool for Weighted Epigenomic Networks in Drosophila melanogaster (Fly T-WEoN), to generate context-specific gene regulatory networks starting from a reference network that contains all known gene regulations in the fly. Unlikely regulations are removed by applying a series of knowledge-based filters. Each of these filters is implemented as an independent module that considers a type of experimental evidence, including DNA methylation, chromatin accessibility, histone modifications and gene expression. Fly T-WEoN is based on heuristic rules that reflect current knowledge on gene regulation in D. melanogaster obtained from the literature. Experimental data files can be generated with several standard procedures and used solely when and if available. Fly T-WEoN is available as a Cytoscape application that permits integration with other tools and facilitates downstream network analysis. In this work, we first demonstrate the reliability of our method to then provide a relevant application case of our tool: early development of D. melanogaster. Fly T-WEoN together with its step-by-step guide is available at https://weon.readthedocs.io.

Early Transcriptional Changes Induced by Wnt/ß-Catenin Signaling in Hippocampal Neurons.

Wnt/ß-catenin signaling modulates brain development and function and its deregulation underlies pathological changes occurring in neurodegenerative and neurodevelopmental disorders. Since one of the main effects of Wnt/ß-catenin signaling is the modulation of target genes, in the present work we examined global transcriptional changes induced by short-term Wnt3a treatment (4 h) in primary cultures of rat hippocampal neurons. RNAseq experiments allowed the identification of 170 differentially expressed genes, including known Wnt/ß-catenin target genes such as Notum, Axin2, and Lef1, as well as novel potential candidates Fam84a, Stk32a, and Itga9. Main biological processes enriched with differentially expressed genes included neural precursor (GO:0061364, p-adjusted = 2.5 × 10-7), forebrain development (GO:0030900, p-adjusted = 7.3 × 10-7), and stem cell differentiation (GO:0048863 p-adjusted = 7.3 × 10-7). Likewise, following activation of the signaling cascade, the expression of a significant number of genes with transcription factor activity (GO:0043565, p-adjusted = 4.1 × 10-6) was induced. We also studied molecular networks enriched upon Wnt3a activation and detected three highly significant expression modules involved in glycerolipid metabolic process (GO:0046486, p-adjusted = 4.5 × 10-19), learning or memory (GO:0007611, p-adjusted = 4.0 × 10-5), and neurotransmitter secretion (GO:0007269, p-adjusted = 5.3 × 10-12). Our results indicate that Wnt/ß-catenin mediated transcription controls multiple biological processes related to neuronal structure and activity that are affected in synaptic dysfunction disorders.

Overrepresentation of glutamate signaling in Alzheimer's disease: network-based pathway enrichment using meta-analysis of genome-wide association studies.

Genome-wide association studies (GWAS) have successfully identified several risk loci for Alzheimer's disease (AD). Nonetheless, these loci do not explain the entire susceptibility of the disease, suggesting that other genetic contributions remain to be identified. Here, we performed a meta-analysis combining data of 4,569 individuals (2,540 cases and 2,029 healthy controls) derived from three publicly available GWAS in AD and replicated a broad genomic region (>248,000 bp) associated with the disease near the APOE/TOMM40 locus in chromosome 19. To detect minor effect size contributions that could help to explain the remaining genetic risk, we conducted network-based pathway analyses either by extracting gene-wise p-values (GW), defined as the single strongest association signal within a gene, or calculated a more stringent gene-based association p-value using the extended Simes (GATES) procedure. Comparison of these strategies revealed that ontological sub-networks (SNs) involved in glutamate signaling were significantly overrepresented in AD (p<2.7×10(-11), p<1.9×10(-11); GW and GATES, respectively). Notably, glutamate signaling SNs were also found to be significantly overrepresented (p<5.1×10(-8)) in the Alzheimer's disease Neuroimaging Initiative (ADNI) study, which was used as a targeted replication sample. Interestingly, components of the glutamate signaling SNs are coordinately expressed in disease-related tissues, which are tightly related to known pathological hallmarks of AD. Our findings suggest that genetic variation within glutamate signaling contributes to the remaining genetic risk of AD and support the notion that functional biological networks should be targeted in future therapies aimed to prevent or treat this devastating neurological disorder.