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.

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.

PersonALL: a genetic scoring guide for personalized risk assessment in pediatric B-cell precursor acute lymphoblastic leukemia.

BACKGROUND: Recurrent genetic lesions provide basis for risk assessment in pediatric acute lymphoblastic leukemia (ALL). However, current prognostic classifiers rely on a limited number of predefined sets of alterations. METHODS: Disease-relevant copy number aberrations (CNAs) were screened genome-wide in 260 children with B-cell precursor ALL. Results were integrated with cytogenetic data to improve risk assessment. RESULTS: CNAs were detected in 93.8% (n = 244) of the patients. First, cytogenetic profiles were combined with IKZF1 status (IKZF1normal, IKZF1del and IKZF1plus) and three prognostic subgroups were distinguished with significantly different 5-year event-free survival (EFS) rates, IKAROS-low (n = 215): 86.3%, IKAROS-medium (n = 27): 57.4% and IKAROS-high (n = 18): 37.5%. Second, contribution of genetic aberrations to the clinical outcome was assessed and an aberration-specific score was assigned to each prognostically relevant alteration. By aggregating the scores of aberrations emerging in individual patients, personalized cumulative values were calculated and used for defining four prognostic subgroups with distinct clinical outcomes. Two favorable subgroups included 60% of patients (n = 157) with a 5-year EFS of 96.3% (excellent risk, n = 105) and 87.2% (good risk, n = 52), respectively; while 40% of patients (n = 103) showed high (n = 74) or ultra-poor (n = 29) risk profile (5-year EFS: 67.4% and 39.0%, respectively). CONCLUSIONS: PersonALL, our conceptually novel prognostic classifier considers all combinations of co-segregating genetic alterations, providing a highly personalized patient stratification.

PhosY-secretome profiling combined with kinase-substrate interaction screening defines active c-Src-driven extracellular signaling.

c-Src tyrosine kinase is a renowned key intracellular signaling molecule and a potential target for cancer therapy. Secreted c-Src is a recent observation, but how it contributes to extracellular phosphorylation remains elusive. Using a series of domain deletion mutants, we show that the N-proximal region of c-Src is essential for its secretion. The tissue inhibitor of metalloproteinases 2 (TIMP2) is an extracellular substrate of c-Src. Limited proteolysis-coupled mass spectrometry and mutagenesis studies verify that the Src homology 3 (SH3) domain of c-Src and the P31VHP34 motif of TIMP2 are critical for their interaction. Comparative phosphoproteomic analyses identify an enrichment of PxxP motifs in phosY-containing secretomes from c-Src-expressing cells with cancer-promoting roles. Inhibition of extracellular c-Src using custom SH3-targeting antibodies disrupt kinase-substrate complexes and inhibit cancer cell proliferation. These findings point toward an intricate role for c-Src in generating phosphosecretomes, which will likely influence cell-cell communication, particularly in c-Src-overexpressing cancers.

Next-Generation Sequencing-Based Genomic Profiling of Children with Acute Myeloid Leukemia.

Pediatric acute myeloid leukemia (AML) represents a major cause of childhood leukemic mortality, with only a limited number of studies investigating the molecular landscape of the disease. Here, we present an integrative analysis of cytogenetic and molecular profiles of 75 patients with pediatric AML from a multicentric, real-world patient cohort treated according to AML Berlin-Frankfurt-Münster protocols. Targeted next-generation sequencing of 54 genes revealed 17 genes that were recurrently mutated in >5% of patients. Considerable differences were observed in the mutational profiles compared with previous studies, as BCORL1, CUX1, KDM6A, PHF6, and STAG2 mutations were detected at a higher frequency than previously reported, whereas KIT, NRAS, and KRAS were less frequently mutated. Our study identified novel recurrent mutations at diagnosis in the BCORL1 gene in 9% of the patients. Tumor suppressor gene (PHF6, TP53, and WT1) mutations were found to be associated with induction failure and shorter event-free survival, suggesting important roles of these alterations in resistance to therapy and disease progression. Comparison of the mutational landscape at diagnosis and relapse revealed an enrichment of mutations in tumor suppressor genes (16.2% versus 44.4%) and transcription factors (35.1% versus 55.6%) at relapse. Our findings shed further light on the heterogeneity of pediatric AML and identify previously unappreciated alterations that may lead to improved molecular characterization and risk stratification of pediatric AML.

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.

Metabolic Adaptation as Potential Target in Papillary Renal Cell Carcinomas Based on Their In Situ Metabolic Characteristics.

Metabolic characteristics of kidney cancers have mainly been obtained from the most frequent clear cell renal cell carcinoma (CCRCC) studies. Moreover, the bioenergetic perturbances that affect metabolic adaptation possibilities of papillary renal cell carcinoma (PRCC) have not yet been detailed. Therefore, our study aimed to analyze the in situ metabolic features of PRCC vs. CCRCC tissues and compared the metabolic characteristics of PRCC, CCRCC, and normal tubular epithelial cell lines. The protein and mRNA expressions of the molecular elements in mammalian target of rapamycin (mTOR) and additional metabolic pathways were analyzed in human PRCC cases compared to CCRCC. The metabolic protein expression pattern, metabolite content, mTOR, and metabolic inhibitor sensitivity of renal carcinoma cell lines were also studied and compared with tubular epithelial cells, as "normal" control. We observed higher protein expressions of the "alternative bioenergetic pathway" elements, in correlation with the possible higher glutamine and acetate consumption in PRCC cells instead of higher glycolytic and mTOR activity in CCRCCs. Increased expression of certain metabolic pathway markers correlates with the detected differences in metabolite ratios, as well. The lower lactate/pyruvate, lactate/malate, and higher pyruvate/citrate intracellular metabolite ratios in PRCC compared to CCRCC cell lines suggest that ACHN (PRCC) have lower Warburg glycolytic capacity, less pronounced pyruvate to lactate producing activity and shifted OXPHOS phenotype. However, both studied renal carcinoma cell lines showed higher mTOR activity than tubular epithelial cells cultured in vitro, the metabolite ratio, the enzyme expression profiles, and the higher mitochondrial content also suggest increased importance of mitochondrial functions, including mitochondrial OXPHOS in PRCCs. Additionally, PRCC cells showed significant mTOR inhibitor sensitivity and the used metabolic inhibitors increased the effect of rapamycin in combined treatments. Our study revealed in situ metabolic differences in mTOR and metabolic protein expression patterns of human PRCC and CCRCC tissues as well as in cell lines. These underline the importance in the development of specific new treatment strategies, new mTOR inhibitors, and other anti-metabolic drug combinations in PRCC therapy.

Transcriptional Analysis-Based Alterations Affecting Neuritogenesis of the Peripheral Nervous System in Psoriasis.

An increasing amount of evidence indicates the critical role of the cutaneous nervous system in the initiation and maintenance of psoriatic skin lesions by neurogenic inflammation. However, molecular mechanisms affecting cutaneous neurons are largely uncharacterized. Therefore, we reanalyzed a psoriatic RNA sequencing dataset from published transcriptome experiments of nearly 300 individuals. Using the Ingenuity Pathway Analysis software, we associated several hundreds of differentially expressed transcripts (DETs) to nervous system development and functions. Since neuronal projections were previously reported to be affected in psoriasis, we performed an in-depth analysis of neurite formation-related process. Our in silico analysis suggests that SEMA-PLXN and ROBO-DCC-UNC5 regulating axonal growth and repulsion are differentially affected in non-lesional and lesional skin samples. We identified opposing expressional alterations in secreted ligands for axonal guidance signaling (RTN4/NOGOA, NTNs, SEMAs, SLITs) and non-conventional axon guidance regulating ligands, including WNT5A and their receptors, modulating axon formation. These differences in neuritogenesis may explain the abnormal cutaneous nerve filament formation described in psoriatic skin. The processes also influence T-cell activation and infiltration, thus highlighting an additional angle of the crosstalk between the cutaneous nervous system and the immune responses in psoriasis pathogenesis, in addition to the known neurogenic pro-inflammatory mediators.

SplicingFactory-splicing diversity analysis for transcriptome data.

MOTIVATION: Alternative splicing contributes to the diversity of RNA found in biological samples. Current tools investigating patterns of alternative splicing check for coordinated changes in the expression or relative ratio of RNA isoforms where specific isoforms are up- or down-regulated in a condition. However, the molecular process of splicing is stochastic and changes in RNA isoform diversity for a gene might arise between samples or conditions. A specific condition can be dominated by a single isoform, while multiple isoforms with similar expression levels can be present in a different condition. These changes might be the result of mutations, drug treatments or differences in the cellular or tissue environment. Here, we present a tool for the characterization and analysis of RNA isoform diversity using isoform level expression measurements. RESULTS: We developed an R package called SplicingFactory, to calculate various RNA isoform diversity metrics, and compare them across conditions. Using the package, we tested the effect of RNA-seq quantification tools, quantification uncertainty, gene expression levels and isoform numbers on the isoform diversity calculation. We analyzed a set of CD34+ hematopoietic stem cells and myelodysplastic syndrome samples and found a set of genes whose isoform diversity change is associated with SF3B1 mutations. AVAILABILITY AND IMPLEMENTATION: The SplicingFactory package is freely available under the GPL-3.0 license from Bioconductor for the Windows, MacOS and Linux operating systems (https://www.bioconductor.org/packages/release/bioc/html/SplicingFactory.html). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Distinct miRNA Expression Signatures of Primary and Secondary Central Nervous System Lymphomas.

Central nervous system (CNS) lymphoma is a rare and aggressive non-Hodgkin lymphoma that might arise in the CNS (primary CNS lymphoma) or disseminates from a systemic lymphoma to the CNS (secondary CNS lymphoma). Dysregulated expression of miRNAs is associated with various pathologic processes, and miRNA expression patterns may have diagnostic, prognostic, and therapeutic implications. However, miRNA expression is understudied in CNS lymphomas. We performed expression analysis of 798 miRNAs in 73 CNS lymphoma samples using the NanoString platform, followed by an analysis to identify potential diagnostic biomarkers characterizing subgroups and to examine differences based on their primary and secondary nature, molecular subtype, mutational patterns, and survival. Thirty-one differentially expressed miRNAs were identified between primary and secondary groups. In addition, 7 more miRNAs were identified associated with a molecular subtype and 25 associated with mutation status. Using unsupervised clustering methods, a small but distinct primary CNS lymphoma subgroup, with characteristically different expression patterns compared with the rest of the cases was defined. Finally, differentially regulated pathways were identified in the above comparisons and the utility of miRNA expression patterns in predicting survival was assessed. Our study identifies a novel CNS lymphoma subgroup defined by distinct miRNAs, proves the importance of specific miRNAs and pathways in the pathogenesis of CNS lymphomas, and provides the basis for future research in defining potential biomarkers.

Screening and monitoring of the BTKC481S mutation in a real-world cohort of patients with relapsed/refractory chronic lymphocytic leukaemia during ibrutinib therapy.

The Bruton's tyrosine kinase (BTK) inhibitor ibrutinib has revolutionised the therapeutic landscape of chronic lymphocytic leukaemia (CLL). Acquired mutations emerging at position C481 in the BTK tyrosine kinase domain are the predominant genetic alterations associated with secondary ibrutinib resistance. To assess the correlation between disease progression, and the emergence and temporal dynamics of the most common resistance mutation BTKC481S , sensitive (10-4 ) time-resolved screening was performed in 83 relapsed/refractory CLL patients during single-agent ibrutinib treatment. With a median follow-up time of 40 months, BTKC481S was detected in 48·2% (40/83) of the patients, with 80·0% (32/40) of them showing disease progression during the examined period. In these 32 cases, representing 72·7% (32/44) of all patients experiencing relapse, emergence of the BTKC481S mutation preceded the symptoms of clinical relapse with a median of nine months. Subsequent Bcl-2 inhibition therapy applied in 28/32 patients harbouring BTKC481S and progressing on ibrutinib conferred clinical and molecular remission across the patients. Our study demonstrates the clinical value of sensitive BTKC481S monitoring with the largest longitudinally analysed real-world patient cohort reported to date and validates the feasibility of an early prediction of relapse in the majority of ibrutinib-treated relapsed/refractory CLL patients experiencing disease progression.

SAMMY-seq reveals early alteration of heterochromatin and deregulation of bivalent genes in Hutchinson-Gilford Progeria Syndrome.

Hutchinson-Gilford progeria syndrome is a genetic disease caused by an aberrant form of Lamin A resulting in chromatin structure disruption, in particular by interfering with lamina associated domains. Early molecular alterations involved in chromatin remodeling have not been identified thus far. Here, we present SAMMY-seq, a high-throughput sequencing-based method for genome-wide characterization of heterochromatin dynamics. Using SAMMY-seq, we detect early stage alterations of heterochromatin structure in progeria primary fibroblasts. These structural changes do not disrupt the distribution of H3K9me3 in early passage cells, thus suggesting that chromatin rearrangements precede H3K9me3 alterations described at later passages. On the other hand, we observe an interplay between changes in chromatin accessibility and Polycomb regulation, with site-specific H3K27me3 variations and transcriptional dysregulation of bivalent genes. We conclude that the correct assembly of lamina associated domains is functionally connected to the Polycomb repression and rapidly lost in early molecular events of progeria pathogenesis.

Postnatal, ontogenic liver growth accomplished by biliary/oval cell proliferation and differentiation.

INTRODUCTION: The liver is well known for its enormous regenerative capacity. If the hepatocytes are compromised the reserve stem cells can regrow the lost tissue by means of oval cells differentiating into hepatocytes. We were curious whether this standby system was able to compensate for ontogenic liver growth arrested by 2-acetylaminofluorene (AAF) treatment or if it can be influenced by cholic acid, known to promote liver growth in several reactions. METHODS: (i) Four weeks-old (60-70g) male F344 rats were kept on standard chow and treated with solvent only, (ii) others were kept on 0,2% cholic acid (CA) enriched diet, (iii) treated with AAF, or (iiii) given a combination of CA diet and AAF treatment (AAF/CA). The proliferative response of epithelial cells was characterized by pulse bromodeoxyuridine labelling. The relative gene expression levels of senescence-related factors and bile acid receptors were determined by quantitative real-time polymerase chain reaction analysis. RESULTS: AAF administration efficiently inhibited the physiological proliferation of hepatocytes in young, male F344 rats after weaning. The activation of stem cells was indicated by the increased proliferation of periportal biliary/oval cells (B/OC). If the rats were fed additionally by cholic acid enriched diet, typical oval cell reaction emerged, subsequently the oval cells differentiated into hepatocytes restituting liver growth. This reaction was mediated by increased production of HGF, IL-6 and SCF by the damaged liver. Moreover, upregulation of FXR expression on B/OC made them competent for bile acids. Our results indicate that endogenous, autocrine mechanisms involved in liver ontogeny are also able to activate the backup regenerative machinery of stem cells.

ATR expands embryonic stem cell fate potential in response to replication stress.

Unrepaired DNA damage during embryonic development can be potentially inherited by a large population of cells. However, the quality control mechanisms that minimize the contribution of damaged cells to developing embryos remain poorly understood. Here, we uncovered an ATR- and CHK1-mediated transcriptional response to replication stress (RS) in mouse embryonic stem cells (ESCs) that induces genes expressed in totipotent two-cell (2C) stage embryos and 2C-like cells. This response is mediated by Dux, a multicopy retrogene defining the cleavage-specific transcriptional program in placental mammals. In response to RS, DUX triggers the transcription of 2C-like markers such as murine endogenous retrovirus-like elements (MERVL) and Zscan4. This response can also be elicited by ETAA1-mediated ATR activation in the absence of RS. ATR-mediated activation of DUX requires GRSF1-dependent post-transcriptional regulation of Dux mRNA. Strikingly, activation of ATR expands ESCs fate potential by extending their contribution to both embryonic and extra-embryonic tissues. These findings define a novel ATR dependent pathway involved in maintaining genome stability in developing embryos by controlling ESCs fate in response to RS.

Large-scale analysis of genome and transcriptome alterations in multiple tumors unveils novel cancer-relevant splicing networks.

Alternative splicing is regulated by multiple RNA-binding proteins and influences the expression of most eukaryotic genes. However, the role of this process in human disease, and particularly in cancer, is only starting to be unveiled. We systematically analyzed mutation, copy number, and gene expression patterns of 1348 RNA-binding protein (RBP) genes in 11 solid tumor types, together with alternative splicing changes in these tumors and the enrichment of binding motifs in the alternatively spliced sequences. Our comprehensive study reveals widespread alterations in the expression of RBP genes, as well as novel mutations and copy number variations in association with multiple alternative splicing changes in cancer drivers and oncogenic pathways. Remarkably, the altered splicing patterns in several tumor types recapitulate those of undifferentiated cells. These patterns are predicted to be mainly controlled by MBNL1 and involve multiple cancer drivers, including the mitotic gene NUMA1 We show that NUMA1 alternative splicing induces enhanced cell proliferation and centrosome amplification in nontumorigenic mammary epithelial cells. Our study uncovers novel splicing networks that potentially contribute to cancer development and progression.

Detection of recurrent alternative splicing switches in tumor samples reveals novel signatures of cancer.

The determination of the alternative splicing isoforms expressed in cancer is fundamental for the development of tumor-specific molecular targets for prognosis and therapy, but it is hindered by the heterogeneity of tumors and the variability across patients. We developed a new computational method, robust to biological and technical variability, which identifies significant transcript isoform changes across multiple samples. We applied this method to more than 4000 samples from the The Cancer Genome Atlas project to obtain novel splicing signatures that are predictive for nine different cancer types, and find a specific signature for basal-like breast tumors involving the tumor-driver CTNND1. Additionally, our method identifies 244 isoform switches, for which the change occurs in the most abundant transcript. Some of these switches occur in known tumor drivers, including PPARG, CCND3, RALGDS, MITF, PRDM1, ABI1 and MYH11, for which the switch implies a change in the protein product. Moreover, some of the switches cannot be described with simple splicing events. Surprisingly, isoform switches are independent of somatic mutations, except for the tumor-suppressor FBLN2 and the oncogene MYH11. Our method reveals novel signatures of cancer in terms of transcript isoforms specifically expressed in tumors, providing novel potential molecular targets for prognosis and therapy. Data and software are available at: http://dx.doi.org/10.6084/m9.figshare.1061917 and https://bitbucket.org/regulatorygenomicsupf/iso-ktsp.

Recombination analysis of Maize dwarf mosaic virus (MDMV) in the Sugarcane mosaic virus (SCMV) subgroup of potyviruses.

Recombination among RNA viruses is a natural phenomenon that appears to have played a significant role in the species development and the evolution of many strains. It also has particular significance for the risk assessment of plants which have been genetically modified for disease resistance by incorporating viral sequences into their genomes. However, the exact recombination events taking place in viral genomes are not investigated in detail for many virus groups. In this analysis, different single-stranded positive-sense RNA potyviruses were compared using various in silico recombination detection methods and new recombination events in the Sugarcane mosaic virus (SCMV) subgroup were detected. For an extended in silico recombination analysis, two of the analyzed Maize dwarf mosaic virus full-length genomes were sequenced additionally during this work. These results strengthen the evidence that recombination is a major driving force in virus evolution, and the emergence of new virus variants in the SCMV subgroup, paired with mutations, could generate viruses with altered biological properties. The intra- and interspecific homolog recombinations seem to be a general trait in this virus group, causing little or no changes to the amino acid of the progenies. However, we found a few breakpoints between the members of SCMV subgroup and the weed-infecting distant relatives, but only a few methods of the RDP3 package predicted these events with low significance level.

Regulation of FAS exon definition and apoptosis by the Ewing sarcoma protein.

The Ewing sarcoma protein EWS is an RNA and DNA binding protein implicated in transcription, pre-mRNA splicing, and DNA damage response. Using CLIP-seq, we identified EWS RNA binding sites in exonic regions near 5' splice sites. A prominent target was exon 6 of the FAS/CD95 receptor, which is alternatively spliced to generate isoforms with opposing activities in programmed cell death. Depletion and overexpression experiments showed that EWS promotes exon 6 inclusion and consequently the synthesis of the proapoptotic FAS/CD95 isoform, whereas an EWS-FLI1 fusion protein characteristic of Ewing sarcomas shows decreased activity. Biochemical analyses revealed that EWS binding promotes the recruitment of U1snRNP and U2AF65 to the splice sites flanking exon 6 and therefore exon definition. Consistent with a role for EWS in the regulation of programmed cell death, cells depleted of EWS show decreased sensitivity to FAS-induced apoptosis, and elevated EWS expression enhances apoptosis in EWS-haploinsufficient Ewing sarcoma cells.