Functional enrichment analysis or gene set enrichment analysis is a basic bioinformatics method that evaluates the biological importance of a list of genes of interest. However, it may produce a long ...list of significant terms with highly redundant information that is difficult to summarize. Current tools to simplify enrichment results by clustering them into groups either still produce redundancy between clusters or do not retain consistent term similarities within clusters. We propose a new method named binary cut for clustering similarity matrices of functional terms. Through comprehensive benchmarks on both simulated and real-world datasets, we demonstrated that binary cut could efficiently cluster functional terms into groups where terms showed consistent similarities within groups and were mutually exclusive between groups. We compared binary cut clustering on the similarity matrices obtained from different similarity measures and found that semantic similarity worked well with binary cut, while similarity matrices based on gene overlap showed less consistent patterns. We implemented the binary cut algorithm in the R package simplifyEnrichment, which additionally provides functionalities for visualizing, summarizing, and comparing the clustering. The simplifyEnrichment package and the documentation are available at https://bioconductor.org/packages/simplifyEnrichment/.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Abstract
Consensus partitioning is an unsupervised method widely used in high-throughput data analysis for revealing subgroups and assigning stability for the classification. However, standard ...consensus partitioning procedures are weak for identifying large numbers of stable subgroups. There are two major issues. First, subgroups with small differences are difficult to be separated if they are simultaneously detected with subgroups with large differences. Second, stability of classification generally decreases as the number of subgroups increases. In this work, we proposed a new strategy to solve these two issues by applying consensus partitioning in a hierarchical procedure. We demonstrated hierarchical consensus partitioning can be efficient to reveal more meaningful subgroups. We also tested the performance of hierarchical consensus partitioning on revealing a great number of subgroups with a large deoxyribonucleic acid methylation dataset. The hierarchical consensus partitioning is implemented in the R package cola with comprehensive functionalities for analysis and visualization. It can also automate the analysis only with a minimum of two lines of code, which generates a detailed HTML report containing the complete analysis. The cola package is available at https://bioconductor.org/packages/cola/.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Cerebral organoids recapitulate the structure and function of the developing human brain in vitro, offering a large potential for personalized therapeutic strategies. The enormous growth of this ...research area over the past decade with its capability for clinical translation makes a non-invasive, automated analysis pipeline of organoids highly desirable. This work presents a novel non-invasive approach to monitor and analyze cerebral organoids over time using high-field magnetic resonance imaging and state-of-the-art tools for automated image analysis. Three specific objectives are addressed, (I) organoid segmentation to investigate organoid development over time, (II) global cysticity classification and (III) local cyst segmentation for organoid quality assessment. We show that organoid growth can be monitored reliably over time and cystic and non-cystic organoids can be separated with high accuracy, with on par or better performance compared to state-of-the-art tools applied to brightfield imaging. Local cyst segmentation is feasible but could be further improved in the future. Overall, these results highlight the potential of the pipeline for clinical application to larger-scale comparative organoid analysis.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK
We studied how intratumoral genetic heterogeneity shapes tumor growth and therapy response for isocitrate dehydrogenase (IDH)-wild-type glioblastoma, a rapidly regrowing tumor. We inferred the ...evolutionary trajectories of matched pairs of primary and relapsed tumors based on deep whole-genome-sequencing data. This analysis suggests both a distant origin of de novo glioblastoma, up to 7 years before diagnosis, and a common path of early tumorigenesis, with one or more of chromosome 7 gain, 9p loss, or 10 loss, at tumor initiation. TERT promoter mutations often occurred later as a prerequisite for rapid growth. In contrast to this common early path, relapsed tumors acquired no stereotypical pattern of mutations and typically regrew from oligoclonal origins, suggesting sparse selective pressure by therapeutic measures.
Display omitted
•We inferred evolutionary trajectories of pairs of primary/relapsed glioblastomas•Chromosome 7 gain, 9p loss, or 10 loss commonly occurred at tumor initiation•TERT promoter mutations often occurred later as a prerequisite for rapid growth•Relapsed tumors typically regrew from oligoclonal origins
By analyzing 21 paired primary and locally relapsed IDH-wild-type glioblastomas (GBM), Körber et al. show that most GBM initiate by gains and losses of specific chromosomes; TERT promoter mutations often occur later as a prerequisite for rapid growth, and relapsed GBM acquire few stereotypical mutations.
Full text
Available for:
GEOZS, IJS, IMTLJ, KILJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBJE, UILJ, UL, UM, UPCLJ, UPUK, ZAGLJ, ZRSKP
Leiomyosarcoma (LMS) is an aggressive mesenchymal malignancy with few therapeutic options. The mechanisms underlying LMS development, including clinically actionable genetic vulnerabilities, are ...largely unknown. Here we show, using whole-exome and transcriptome sequencing, that LMS tumors are characterized by substantial mutational heterogeneity, near-universal inactivation of TP53 and RB1, widespread DNA copy number alterations including chromothripsis, and frequent whole-genome duplication. Furthermore, we detect alternative telomere lengthening in 78% of cases and identify recurrent alterations in telomere maintenance genes such as ATRX, RBL2, and SP100, providing insight into the genetic basis of this mechanism. Finally, most tumors display hallmarks of "BRCAness", including alterations in homologous recombination DNA repair genes, multiple structural rearrangements, and enrichment of specific mutational signatures, and cultured LMS cells are sensitive towards olaparib and cisplatin. This comprehensive study of LMS genomics has uncovered key biological features that may inform future experimental research and enable the design of novel therapies.
Abstract
Inherited disorders of neurotransmitter metabolism are rare neurodevelopmental diseases presenting with movement disorders and global developmental delay. This study presents the results of ...the first standardized deep phenotyping approach and describes the clinical and biochemical presentation at disease onset as well as diagnostic approaches of 275 patients from the registry of the International Working Group on Neurotransmitter related Disorders. The results reveal an increased rate of prematurity, a high risk for being small for gestational age and for congenital microcephaly in some disorders. Age at diagnosis and the diagnostic delay are influenced by the diagnostic methods applied and by disease-specific symptoms. The timepoint of investigation was also a significant factor: delay to diagnosis has decreased in recent years, possibly due to novel diagnostic approaches or raised awareness. Although each disorder has a specific biochemical pattern, we observed confounding exceptions to the rule. The data provide comprehensive insights into the phenotypic spectrum of neurotransmitter disorders.
Chromothripsis is a recently identified mutational phenomenon, by which a presumably single catastrophic event generates extensive genomic rearrangements of one or a few chromosome(s). Considered as ...an early event in tumour development, this form of genome instability plays a prominent role in tumour onset. Chromothripsis prevalence might have been underestimated when using low-resolution methods, and pan-cancer studies based on sequencing are rare. Here we analyse chromothripsis in 28 tumour types covering all major adult cancers (634 tumours, 316 whole-genome and 318 whole-exome sequences). We show that chromothripsis affects a substantial proportion of human cancers, with a prevalence of 49% across all cases. Chromothripsis generates entity-specific genomic alterations driving tumour development, including clinically relevant druggable fusions. Chromothripsis is linked with specific telomere patterns and univocal mutational signatures in distinct tumour entities. Longitudinal analysis of chromothriptic patterns in 24 matched tumour pairs reveals insights in the clonal evolution of tumours with chromothripsis.
DNA methylation patterns delineate clinically relevant subgroups of meningioma. We previously established the six meningioma methylation classes (MC) benign 1–3, intermediate A and B, and malignant. ...Here, we set out to identify subgroup-specific mutational patterns and gene regulation. Whole genome sequencing was performed on 62 samples across all MCs and WHO grades from 62 patients with matched blood control, including 40 sporadic meningiomas and 22 meningiomas arising after radiation (Mrad). RNA sequencing was added for 18 of these cases and chromatin-immunoprecipitation for histone H3 lysine 27 acetylation (H3K27ac) followed by sequencing (ChIP-seq) for 16 samples. Besides the known mutations in meningioma, structural variants were found as the mechanism of
NF2
inactivation in a small subset (5%) of sporadic meningiomas, similar to previous reports for Mrad. Aberrations of
DMD
were found to be enriched in MCs with
NF2
mutations, and
DMD
was among the most differentially upregulated genes in
NF2
mutant compared to
NF2
wild-type cases. The mutational signature AC3, which has been associated with defects in homologous recombination repair (HRR), was detected in both sporadic meningioma and Mrad, but widely distributed across the genome in sporadic cases and enriched near genomic breakpoints in Mrad. Compared to the other MCs, the number of single nucleotide variants matching the AC3 pattern was significantly higher in the malignant MC, which also exhibited higher genomic instability, determined by the numbers of both large segments affected by copy number alterations and breakpoints between large segments. ChIP-seq analysis for H3K27ac revealed a specific activation of genes regulated by the transcription factor FOXM1 in the malignant MC. This analysis also revealed a super enhancer near the
HOXD
gene cluster in this MC, which, together with general upregulation of
HOX
genes in the malignant MC, indicates a role of
HOX
genes in meningioma aggressiveness. This data elucidates the biological mechanisms rendering different epigenetic subgroups of meningiomas, and suggests leveraging HRR as a novel therapeutic target.
Full text
Available for:
EMUNI, FIS, FZAB, GEOZS, GIS, IJS, IMTLJ, KILJ, KISLJ, MFDPS, NLZOH, NUK, OBVAL, OILJ, PNG, SAZU, SBCE, SBJE, SBMB, SBNM, UKNU, UL, UM, UPUK, VKSCE, ZAGLJ
Single-cell genomics technology has transformed our understanding of complex cellular systems. However, excessive cost and a lack of strategies for the purification of newly identified cell types ...impede their functional characterization and large-scale profiling. Here, we have generated high-content single-cell proteo-genomic reference maps of human blood and bone marrow that quantitatively link the expression of up to 197 surface markers to cellular identities and biological processes across all main hematopoietic cell types in healthy aging and leukemia. These reference maps enable the automatic design of cost-effective high-throughput cytometry schemes that outperform state-of-the-art approaches, accurately reflect complex topologies of cellular systems and permit the purification of precisely defined cell states. The systematic integration of cytometry and proteo-genomic data enables the functional capacities of precisely mapped cell states to be measured at the single-cell level. Our study serves as an accessible resource and paves the way for a data-driven era in cytometry.
Full text
Available for:
GEOZS, IJS, IMTLJ, KISLJ, NLZOH, NUK, OILJ, PNG, SAZU, SBCE, SBMB, UL, UM, UPUK, ZAGLJ
Chordomas are rare bone tumors with few therapeutic options. Here we show, using whole-exome and genome sequencing within a precision oncology program, that advanced chordomas (n = 11) may be ...characterized by genomic patterns indicative of defective homologous recombination (HR) DNA repair and alterations affecting HR-related genes, including, for example, deletions and pathogenic germline variants of BRCA2, NBN, and CHEK2. A mutational signature associated with HR deficiency was significantly enriched in 72.7% of samples and co-occurred with genomic instability. The poly(ADP-ribose) polymerase (PARP) inhibitor olaparib, which is preferentially toxic to HR-incompetent cells, led to prolonged clinical benefit in a patient with refractory chordoma, and whole-genome analysis at progression revealed a PARP1 p.T910A mutation predicted to disrupt the autoinhibitory PARP1 helical domain. These findings uncover a therapeutic opportunity in chordoma that warrants further exploration, and provide insight into the mechanisms underlying PARP inhibitor resistance.