HELLS regulates transcription in T-cell lymphomas by reducing unscheduled R-loops and by facilitating RNAPII progression.

Chromatin modifiers are emerging as major determinants of many types of cancers, including Anaplastic Large Cell Lymphomas (ALCL), a family of highly heterogeneous T-cell lymphomas for which therapeutic options are still limited. HELLS is a multifunctional chromatin remodeling protein that affects genomic instability by participating in the DNA damage response. Although the transcriptional function of HELLS has been suggested, no clues on how HELLS controls transcription are currently available. In this study, by integrating different multi-omics and functional approaches, we characterized the transcriptional landscape of HELLS in ALCL. We explored the clinical impact of its transcriptional program in a large cohort of 44 patients with ALCL. We demonstrated that HELLS, loaded at the level of intronic regions of target promoters, facilitates RNA Polymerase II (RNAPII) progression along the gene bodies by reducing the persistence of co-transcriptional R-loops and promoting DNA damage resolution. Importantly, selective knockdown of HELLS sensitizes ALCL cells to different chemotherapeutic agents, showing a synergistic effect. Collectively, our work unveils the role of HELLS in acting as a gatekeeper of ALCL genome stability providing a rationale for drug design.

Biochemical properties of chromatin domains define genome compartmentalization.

Chromatin three-dimensional (3D) organization inside the cell nucleus determines the separation of euchromatin and heterochromatin domains. Their segregation results in the definition of active and inactive chromatin compartments, whereby the local concentration of associated proteins, RNA and DNA results in the formation of distinct subnuclear structures. Thus, chromatin domains spatially confined in a specific 3D nuclear compartment are expected to share similar epigenetic features and biochemical properties, in terms of accessibility and solubility. Based on this rationale, we developed the 4f-SAMMY-seq to map euchromatin and heterochromatin based on their accessibility and solubility, starting from as little as 10 000 cells. Adopting a tailored bioinformatic data analysis approach we reconstruct also their 3D segregation in active and inactive chromatin compartments and sub-compartments, thus recapitulating the characteristic properties of distinct chromatin states. A key novelty of the new method is the capability to map both the linear segmentation of open and closed chromatin domains, as well as their compartmentalization in one single experiment.

CHD8 suppression impacts on histone H3 lysine 36 trimethylation and alters RNA alternative splicing.

Disruptive mutations in the chromodomain helicase DNA-binding protein 8 gene (CHD8) have been recurrently associated with autism spectrum disorders (ASDs). Here we investigated how chromatin reacts to CHD8 suppression by analyzing a panel of histone modifications in induced pluripotent stem cell-derived neural progenitors. CHD8 suppression led to significant reduction (47.82%) in histone H3K36me3 peaks at gene bodies, particularly impacting on transcriptional elongation chromatin states. H3K36me3 reduction specifically affects highly expressed, CHD8-bound genes and correlates with altered alternative splicing patterns of 462 genes implicated in 'regulation of RNA splicing' and 'mRNA catabolic process'. Mass spectrometry analysis uncovered a novel interaction between CHD8 and the splicing regulator heterogeneous nuclear ribonucleoprotein L (hnRNPL), providing the first mechanistic insights to explain the CHD8 suppression-derived splicing phenotype, partly implicating SETD2, a H3K36me3 methyltransferase. In summary, our results point toward broad molecular consequences of CHD8 suppression, entailing altered histone deposition/maintenance and RNA processing regulation as important regulatory processes in ASD.

Leveraging three-dimensional chromatin architecture for effective reconstruction of enhancer-target gene regulatory interactions.

A growing amount of evidence in literature suggests that germline sequence variants and somatic mutations in non-coding distal regulatory elements may be crucial for defining disease risk and prognostic stratification of patients, in genetic disorders as well as in cancer. Their functional interpretation is challenging because genome-wide enhancer-target gene (ETG) pairing is an open problem in genomics. The solutions proposed so far do not account for the hierarchy of structural domains which define chromatin three-dimensional (3D) architecture. Here we introduce a change of perspective based on the definition of multi-scale structural chromatin domains, integrated in a statistical framework to define ETG pairs. In this work (i) we develop a computational and statistical framework to reconstruct a comprehensive map of ETG pairs leveraging functional genomics data; (ii) we demonstrate that the incorporation of chromatin 3D architecture information improves ETG pairing accuracy and (iii) we use multiple experimental datasets to extensively benchmark our method against previous solutions for the genome-wide reconstruction of ETG pairs. This solution will facilitate the annotation and interpretation of sequence variants in distal non-coding regulatory elements. We expect this to be especially helpful in clinically oriented applications of whole genome sequencing in cancer and undiagnosed genetic diseases research.

SourceSet: A graphical model approach to identify primary genes in perturbed biological pathways.

Topological gene-set analysis has emerged as a powerful means for omic data interpretation. Although numerous methods for identifying dysregulated genes have been proposed, few of them aim to distinguish genes that are the real source of perturbation from those that merely respond to the signal dysregulation. Here, we propose a new method, called SourceSet, able to distinguish between the primary and the secondary dysregulation within a Gaussian graphical model context. The proposed method compares gene expression profiles in the control and in the perturbed condition and detects the differences in both the mean and the covariance parameters with a series of likelihood ratio tests. The resulting evidence is used to infer the primary and the secondary set, i.e. the genes responsible for the primary dysregulation, and the genes affected by the perturbation through network propagation. The proposed method demonstrates high specificity and sensitivity in different simulated scenarios and on several real biological case studies. In order to fit into the more traditional pathway analysis framework, SourceSet R package also extends the analysis from a single to multiple pathways and provides several graphical outputs, including Cytoscape visualization to browse the results.

simPATHy: a new method for simulating data from perturbed biological PATHways.

Summary: In the omic era, one of the main aims is to discover groups of functionally related genes that drive the difference between different conditions. To this end, a plethora of potentially useful multivariate statistical approaches has been proposed, but their evaluation is hindered by the absence of a gold standard. Here, we propose a method for simulating biological data ­ gene expression, RPKM/FPKM or protein abundances ­ from two conditions, namely, a reference condition and a perturbation of it. Our approach is built upon probabilistic graphical models and is thus especially suited for testing topological approaches. Availability and Implementation: The simPATHy is an R package, it is open source and freely available on CRAN. Contacts: elisa.salviato.2@studenti.unipd.it or chiara.romualdi@unipd.it Supplementary Information: Supplementary data are available at Bioinformatics online.

Exploring the transcriptional cooperation between RUNX2 and its associated elncRNA RAIN.

Recent insights into the mechanisms controlling gene expression identified enhancer-associated long non-coding RNAs (elncRNAs) as master players of transcription in cancers. RUNX2, a mammalian RUNT-related transcription factor, is increasingly recognized in cancer biology for its role in supporting survival and progression also in thyroid cancer (TC). We recently identified, within the RUNX2 locus, a novel elncRNA that we named RAIN (RUNX2 associated intergenic lncRNA). We showed that RAIN and RUNX2 expression correlate in TC, both in vitro and in vivo, and that RAIN promotes RUNX2 expression by interacting with and affecting the activity of the RUNX2 P2 promoter through two distinct mechanisms. Here, we took forward these observations to explore the genome-wide transcriptional function of RAIN and its contribution to the RUNX2-dependent gene expression program in TC. By combining multiple omics data, we demonstrated that RAIN functionally cooperates with RUNX2 to the regulation of a subset of functionally related genes involved in promoting matrix remodeling, migration, and loss of differentiation. We showed that RAIN interacts with RUNX2 and its expression is required for the efficient recruitment of this TF to its target regulatory regions. In addition, our data revealed that besides RUNX2, RAIN governs a hierarchically organized complex transcriptional program by controlling a core of cancer-associated TFs that, in turn, orchestrate the expression of downstream genes. This evidence indicates that the functional cooperation observed between RAIN and RUNX2 can be a diffuse work mechanism for this elncRNA.

Leveraging Tissue-Specific Enhancer-Target Gene Regulatory Networks Identifies Enhancer Somatic Mutations That Functionally Impact Lung Cancer.

Enhancers are noncoding regulatory DNA regions that modulate the transcription of target genes, often over large distances along with the genomic sequence. Enhancer alterations have been associated with various pathological conditions, including cancer. However, the identification and characterization of somatic mutations in noncoding regulatory regions with a functional effect on tumorigenesis and prognosis remain a major challenge. Here, we present a strategy for detecting and characterizing enhancer mutations in a genome-wide analysis of patient cohorts, across three lung cancer subtypes. Lung tissue-specific enhancers were defined by integrating experimental data and public epigenomic profiles, and the genome-wide enhancer-target gene regulatory network of lung cells was constructed by integrating chromatin three-dimensional architecture data. Lung cancers possessed a similar mutation burden at tissue-specific enhancers and exons but with differences in their mutation signatures. Functionally relevant alterations were prioritized on the basis of the pathway-level integration of the effect of a mutation and the frequency of mutations on individual enhancers. The genes enriched for mutated enhancers converged on the regulation of key biological processes and pathways relevant to tumor biology. Recurrent mutations in individual enhancers also affected the expression of target genes, with potential relevance for patient prognosis. Together, these findings show that noncoding regulatory mutations have a potential relevance for cancer pathogenesis and can be exploited for patient classification. SIGNIFICANCE: Mapping enhancer-target gene regulatory interactions and analyzing enhancer mutations at the level of their target genes and pathways reveal convergence of recurrent enhancer mutations on biological processes involved in tumorigenesis and prognosis.

Genome-wide study of salivary miRNAs identifies miR-423-5p as promising diagnostic and prognostic biomarker in oral squamous cell carcinoma.

Survival rates of oral squamous cell carcinoma (OSCC) remained substantially unchanged over the last decades; thus, additional prognostic tools are strongly needed. Salivary miRNAs have emerged as excellent non-invasive cancer biomarker candidates, but their association with OSCC prognosis has not been investigated yet. In this study, we analyzed global salivary miRNA expression in OSCC patients and healthy controls, with the aim to define its diagnostic and prognostic potential. Methods: Saliva was collected from patients with newly diagnosed untreated primary OSCC and healthy controls. Global profiling of salivary miRNAs was carried out through a microarray approach, while signature validation was performed by quantitative real-time PCR (RT-qPCR). A stringent statistical approach for microarray and RT-qPCR data normalization was applied. The diagnostic performance of miRNAs and their correlation with OSCC prognosis were comprehensively analyzed. Results: In total, 25 miRNAs emerged as differentially expressed between OSCC patients and healthy controls and, among them, seven were significantly associated with disease-free survival (DFS). miR-106b-5p, miR-423-5p and miR-193b-3p were expressed at high levels in saliva of OSCC patients and their combination displays the best diagnostic performance (ROC - AUC = 0.98). Moreover, high expression of miR-423-5p was an independent predictor of poor DFS, when included in multivariate survival analysis with the number of positive lymph nodes - the only significant clinical prognosticator. Finally, we observed a significant decrease in miR-423-5p expression in matched post-operative saliva samples, suggesting its potential cancer-specific origin. Conclusion: Salivary miRNAs identified in our cohort of patients show to be accurate in OSCC detection and to effectively stratify patients according to their likelihood of relapse. These results, if validated in an independent set of patients, could be particularly promising for screening/follow-up of high-risk populations and useful for preoperative prognostic assessment.

Comprehensive Profiling of Hypoxia-Related miRNAs Identifies miR-23a-3p Overexpression as a Marker of Platinum Resistance and Poor Prognosis in High-Grade Serous Ovarian Cancer.

The onset of chemo-resistant recurrence represents the principal cause of high-grade serous ovarian carcinoma (HGSOC) death. HGSOC masses are characterized by a hypoxic microenvironment, which contributes to the development of this chemo-resistant phenotype. Hypoxia regulated-miRNAs (HRMs) represent a molecular response of cancer cells to hypoxia and are involved in tumor progression. We investigated the expression of HRMs using miRNA expression data from a total of 273 advanced-stage HGSOC samples. The miRNAs associated with chemoresistance and survival were validated by RT-qPCR and target prediction, and comparative pathway analysis was conducted for target gene identification. Analysis of miRNA expression profiles indicated miR-23a-3p and miR-181c-5p over-expression as associated with chemoresistance and poor PFS. RT-qPCR data confirmed upregulation of miR-23a-3p in tumors from chemoresistant HGSOC patients and its significant association with shorter PFS. In silico miR-23a-3p target prediction and comparative pathway analysis identified platinum drug resistance as the pathway with the highest number of miR-23a-3p target genes. Among them, APAF-1 emerged as the most promising, being downregulated in platinum-resistant patients and in HGSOC chemo-resistant cells. These results highlight miR-23a-3p as a potential biomarker for HGSOC platinum response and prognosis and the miR23a-3p/APAF1 axis as a possible target to overcome platinum-resistance.

Oncolytic adenovirus and gene therapy with EphA2-BiTE for the treatment of pediatric high-grade gliomas.

BACKGROUND: Pediatric high-grade gliomas (pHGGs) are among the most common and incurable malignant neoplasms of childhood. Despite aggressive, multimodal treatment, the outcome of children with high-grade gliomas has not significantly improved over the past decades, prompting the development of innovative approaches. METHODS: To develop an effective treatment, we aimed at improving the suboptimal antitumor efficacy of oncolytic adenoviruses (OAs) by testing the combination with a gene-therapy approach using a bispecific T-cell engager (BiTE) directed towards the erythropoietin-producing human hepatocellular carcinoma A2 receptor (EphA2), conveyed by a replication-incompetent adenoviral vector (EphA2 adenovirus (EAd)). The combinatorial approach was tested in vitro, in vivo and thoroughly characterized at a molecular level. RESULTS: After confirming the relevance of EphA2 as target in pHGGs, documenting a significant correlation with worse clinical outcome of the patients, we showed that the proposed strategy provides significant EphA2-BiTE amplification and enhanced tumor cell apoptosis, on coculture with T cells. Moreover, T-cell activation through an agonistic anti-CD28 antibody further increased the activation/proliferation profiles and functional response against infected tumor cells, inducing eradication of highly resistant, primary pHGG cells. The gene-expression analysis of tumor cells and T cells, after coculture, revealed the importance of both EphA2-BiTE and costimulation in the proposed system. These in vitro observations translated into significant tumor control in vivo, in both subcutaneous and a more challenging orthotopic model. CONCLUSIONS: The combination of OA and EphA2-BiTE gene therapy strongly enhances the antitumor activity of OA, inducing the eradication of highly resistant tumor cells, thus supporting the clinical translation of the approach.

Global chromatin conformation differences in the Drosophila dosage compensated chromosome X.

In Drosophila melanogaster the single male chromosome X undergoes an average twofold transcriptional upregulation for balancing the transcriptional output between sexes. Previous literature hypothesised that a global change in chromosome structure may accompany this process. However, recent studies based on Hi-C failed to detect these differences. Here we show that global conformational differences are specifically present in the male chromosome X and detectable using Hi-C data on sex-sorted embryos, as well as male and female cell lines, by leveraging custom data analysis solutions. We find the male chromosome X has more mid-/long-range interactions. We also identify differences at structural domain boundaries containing BEAF-32 in conjunction with CP190 or Chromator. Weakening of these domain boundaries in male chromosome X co-localizes with the binding of the dosage compensation complex and its co-factor CLAMP, reported to enhance chromatin accessibility. Together, our data strongly indicate that chromosome X dosage compensation affects global chromosome structure.

Circulating miRNA landscape identifies miR-1246 as promising diagnostic biomarker in high-grade serous ovarian carcinoma: A validation across two independent cohorts.

High-grade serous ovarian carcinoma (HGSOC) is the most lethal gynecologic neoplasm, with five-year survival rate below 30%. Early disease detection is of utmost importance to improve HGSOC cure rate. Sera from 168 HGSOC patients and 65 healthy controls were gathered together from two independent collections and stratified into a training set, for miRNA marker identification, and a validation set, for data validation. An innovative statistical approach for microarray data normalization was developed to identify differentially expressed miRNAs. Signature validation in both the training and validation sets was performed by quantitative Real Time PCR (RT-qPCR). In both the training and validation sets, miR-1246, miR-595 and miR-2278 emerged significantly over expressed in the sera of HGSOC patients compared to healthy controls. Receiver Operating Characteristic curve analysis revealed miR-1246 as the best diagnostic biomarker, with a sensitivity of 87%, a specificity of 77% and an accuracy of 84%. This study is the first step in the identification of circulating miRNAs with diagnostic relevance for HGSOC. According to its specificity and sensitivity, circulating miR-1246 levels are worthy to be further investigated as potential diagnostic biomarker for HGSOC.