The function of many eukaryotic proteins is regulated by highly dynamic changes in their nucleocytoplasmic distribution. The ability to precisely and reversibly control nuclear translocation would, therefore, allow dissecting and engineering cellular networks. Here we develop a genetically encoded, light-inducible nuclear localization signal (LINuS) based on the LOV2 domain of Avena sativa phototropin 1. LINuS is a small, versatile tag, customizable for different proteins and cell types. LINuS-mediated nuclear import is fast and reversible, and can be tuned at different levels, for instance, by introducing mutations that alter AsLOV2 domain photo-caging properties or by selecting nuclear localization signals (NLSs) of various strengths. We demonstrate the utility of LINuS in mammalian cells by controlling gene expression and entry into mitosis with blue light.

Integrative analysis of multi-omics layers at single cell level is critical for accurate dissection of cell-to-cell variation within certain cell populations. Here we report scCAT-seq, a technique for simultaneously assaying chromatin accessibility and the transcriptome within the same single cell. We show that the combined single cell signatures enable accurate construction of regulatory relationships between cis-regulatory elements and the target genes at single-cell resolution, providing a new dimension of features that helps direct discovery of regulatory patterns specific to distinct cell identities. Moreover, we generate the first single cell integrated map of chromatin accessibility and transcriptome in early embryos and demonstrate the robustness of scCAT-seq in the precise dissection of master transcription factors in cells of distinct states. The ability to obtain these two layers of omics data will help provide more accurate definitions of "single cell state" and enable the deconvolution of regulatory heterogeneity from complex cell populations.

Glioblastoma frequently exhibits therapy-associated subtype transitions to mesenchymal phenotypes with adverse prognosis. Here, we perform multi-omic profiling of 60 glioblastoma primary tumours and use orthogonal analysis of chromatin and RNA-derived gene regulatory networks to identify 38 subtype master regulators, whose cell population-specific activities we further map in published single-cell RNA sequencing data. These analyses identify the oligodendrocyte precursor marker and chromatin modifier SOX10 as a master regulator in RTK I-subtype tumours. In vitro functional studies demonstrate that SOX10 loss causes a subtype switch analogous to the proneural-mesenchymal transition observed in patients at the transcriptomic, epigenetic and phenotypic levels. SOX10 repression in an in vivo syngeneic graft glioblastoma mouse model results in increased tumour invasion, immune cell infiltration and significantly reduced survival, reminiscent of progressive human glioblastoma. These results identify SOX10 as a bona fide master regulator of the RTK I subtype, with both tumour cell-intrinsic and microenvironmental effects.

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.

Multiplexed fluorescence in situ hybridization techniques have enabled cell-type identification, linking transcriptional heterogeneity with spatial heterogeneity of cells. However, inaccurate cell segmentation reduces the efficacy of cell-type identification and tissue characterization. Here, we present a method called Spot-based Spatial cell-type Analysis by Multidimensional mRNA density estimation (SSAM), a robust cell segmentation-free computational framework for identifying cell-types and tissue domains in 2D and 3D. SSAM is applicable to a variety of in situ transcriptomics techniques and capable of integrating prior knowledge of cell types. We apply SSAM to three mouse brain tissue images: the somatosensory cortex imaged by osmFISH, the hypothalamic preoptic region by MERFISH, and the visual cortex by multiplexed smFISH. Here, we show that SSAM detects regions occupied by known cell types that were previously missed and discovers new cell types.

Bringing together cancer genomes from different projects increases power and allows the investigation of pan-cancer, molecular mechanisms. However, working with whole genomes sequenced over several years in different sequencing centres requires a framework to compare the quality of these sequences. We used the Pan-Cancer Analysis of Whole Genomes cohort as a test case to construct such a framework. This cohort contains whole cancer genomes of 2832 donors from 18 sequencing centres. We developed a non-redundant set of five quality control (QC) measurements to establish a star rating system. These QC measures reflect known differences in sequencing protocol and provide a guide to downstream analyses and allow for exclusion of samples of poor quality. We have found that this is an effective framework of quality measures. The implementation of the framework is available at: https://dockstore.org/containers/quay.io/jwerner_dkfz/pancanqc:1.2.2 .

Pancreatic ductal adenocarcinoma (PDAC) is projected to be the second leading cause of cancer mortality by 2030. Bulk transcriptomic analyses have distinguished 'classical' from 'basal-like' tumors with more aggressive clinical behavior. We derive PDAC organoids from 18 primary tumors and two matched liver metastases, and show that 'classical' and 'basal-like' cells coexist in individual organoids. By single-cell transcriptome analysis of PDAC organoids and primary PDAC, we identify distinct tumor cell states shared across patients, including a cycling progenitor cell state and a differentiated secretory state. Cell states are connected by a differentiation hierarchy, with 'classical' cells concentrated at the endpoint. In an imaging-based drug screen, expression of 'classical' subtype genes correlates with better drug response. Our results thus uncover a functional hierarchy of PDAC cell states linked to transcriptional tumor subtypes, and support the use of PDAC organoids as a clinically relevant model for in vitro studies of tumor heterogeneity.

Intratumoural heterogeneity (ITH) is a major cause of cancer-associated lethality. Extensive genomic ITH has previously been reported in clear cell renal cell carcinoma (ccRCC). Here we address the question whether ITH increases with malignant progression and can hence be exploited as a prognostic marker. Unexpectedly, precision quantitative image analysis reveals that the degree of functional ITH is virtually identical between primary ccRCCs of the lowest stage and advanced, metastatic tumours. Functional ITH was found to show a stage-independent topological pattern with peak proliferative and signalling activities almost exclusively in the tumour periphery. Exome sequencing of matching peripheral and central primary tumour specimens reveals various region-specific mutations. However, these mutations cannot directly explain the zonal pattern suggesting a role of microenvironmental factors in shaping functional ITH. In conclusion, our results indicate that ITH is an early and general characteristic of malignant growth rather than a consequence of malignant progression.

The circadian clock and the cell cycle are major cellular systems that organize global physiology in temporal fashion. It seems conceivable that the potentially conflicting programs are coordinated. We show here that overexpression of MYC in U2OS cells attenuates the clock and conversely promotes cell proliferation while downregulation of MYC strengthens the clock and reduces proliferation. Inhibition of the circadian clock is crucially dependent on the formation of repressive complexes of MYC with MIZ1 and subsequent downregulation of the core clock genes BMAL1 (ARNTL), CLOCK and NPAS2. We show furthermore that BMAL1 expression levels correlate inversely with MYC levels in 102 human lymphomas. Our data suggest that MYC acts as a master coordinator that inversely modulates the impact of cell cycle and circadian clock on gene expression.

Active nucleocytoplasmic transport is a key mechanism underlying protein regulation in eukaryotes. While nuclear protein import can be controlled in space and time with a portfolio of optogenetic tools, protein export has not been tackled so far. Here we present a light-inducible nuclear export system (LEXY) based on a single, genetically encoded tag, which enables precise spatiotemporal control over the export of tagged proteins. A constitutively nuclear, chromatin-anchored LEXY variant expands the method towards light inhibition of endogenous protein export by sequestering cellular CRM1 receptors. We showcase the utility of LEXY for cell biology applications by regulating a synthetic repressor as well as human p53 transcriptional activity with light. LEXY is a powerful addition to the optogenetic toolbox, allowing various novel applications in synthetic and cell biology.

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.

Burkitt lymphoma (BL) is the most common B-cell lymphoma in children. Within the International Cancer Genome Consortium (ICGC), we performed whole genome and transcriptome sequencing of 39 sporadic BL. Here, we unravel interaction of structural, mutational, and transcriptional changes, which contribute to MYC oncogene dysregulation together with the pathognomonic IG-MYC translocation. Moreover, by mapping IGH translocation breakpoints, we provide evidence that the precursor of at least a subset of BL is a B-cell poised to express IGHA. We describe the landscape of mutations, structural variants, and mutational processes, and identified a series of driver genes in the pathogenesis of BL, which can be targeted by various mechanisms, including IG-non MYC translocations, germline and somatic mutations, fusion transcripts, and alternative splicing.

Incomplete understanding of the metastatic process hinders personalized therapy. Here we report the most comprehensive whole-genome study of colorectal metastases vs. matched primary tumors. 65% of somatic mutations originate from a common progenitor, with 15% being tumor- and 19% metastasis-specific, implicating a higher mutation rate in metastases. Tumor- and metastasis-specific mutations harbor elevated levels of BRCAness. We confirm multistage progression with new components ARHGEF7/ARHGEF33. Recurrently mutated non-coding elements include ncRNAs RP11-594N15.3, AC010091, SNHG14, 3' UTRs of FOXP2, DACH2, TRPM3, XKR4, ANO5, CBL, CBLB, the latter four potentially dual protagonists in metastasis and efferocytosis-/PD-L1 mediated immunosuppression. Actionable metastasis-specific lesions include FAT1, FGF1, BRCA2, KDR, and AKT2-, AKT3-, and PDGFRA-3' UTRs. Metastasis specific mutations are enriched in PI3K-Akt signaling, cell adhesion, ECM and hepatic stellate activation genes, suggesting genetic programs for site-specific colonization. Our results put forward hypotheses on tumor and metastasis evolution, and evidence for metastasis-specific events relevant for personalized therapy.

Cancers require telomere maintenance mechanisms for unlimited replicative potential. They achieve this through TERT activation or alternative telomere lengthening associated with ATRX or DAXX loss. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, we dissect whole-genome sequencing data of over 2500 matched tumor-control samples from 36 different tumor types aggregated within the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium to characterize the genomic footprints of these mechanisms. While the telomere content of tumors with ATRX or DAXX mutations (ATRX/DAXXtrunc) is increased, tumors with TERT modifications show a moderate decrease of telomere content. One quarter of all tumor samples contain somatic integrations of telomeric sequences into non-telomeric DNA. This fraction is increased to 80% prevalence in ATRX/DAXXtrunc tumors, which carry an aberrant telomere variant repeat (TVR) distribution as another genomic marker. The latter feature includes enrichment or depletion of the previously undescribed singleton TVRs TTCGGG and TTTGGG, respectively. Our systematic analysis provides new insight into the recurrent genomic alterations associated with telomere maintenance mechanisms in cancer.

Viruses manipulate cellular metabolism and macromolecule recycling processes like autophagy. Dysregulated metabolism might lead to excessive inflammatory and autoimmune responses as observed in severe and long COVID-19 patients. Here we show that SARS-CoV-2 modulates cellular metabolism and reduces autophagy. Accordingly, compound-driven induction of autophagy limits SARS-CoV-2 propagation. In detail, SARS-CoV-2-infected cells show accumulation of key metabolites, activation of autophagy inhibitors (AKT1, SKP2) and reduction of proteins responsible for autophagy initiation (AMPK, TSC2, ULK1), membrane nucleation, and phagophore formation (BECN1, VPS34, ATG14), as well as autophagosome-lysosome fusion (BECN1, ATG14 oligomers). Consequently, phagophore-incorporated autophagy markers LC3B-II and P62 accumulate, which we confirm in a hamster model and lung samples of COVID-19 patients. Single-nucleus and single-cell sequencing of patient-derived lung and mucosal samples show differential transcriptional regulation of autophagy and immune genes depending on cell type, disease duration, and SARS-CoV-2 replication levels. Targeting of autophagic pathways by exogenous administration of the polyamines spermidine and spermine, the selective AKT1 inhibitor MK-2206, and the BECN1-stabilizing anthelmintic drug niclosamide inhibit SARS-CoV-2 propagation in vitro with IC50 values of 136.7, 7.67, 0.11, and 0.13 μM, respectively. Autophagy-inducing compounds reduce SARS-CoV-2 propagation in primary human lung cells and intestinal organoids emphasizing their potential as treatment options against COVID-19.

The molecular pathogenesis of salivary gland acinic cell carcinoma (AciCC) is poorly understood. The secretory Ca-binding phosphoprotein (SCPP) gene cluster at 4q13 encodes structurally related phosphoproteins of which some are specifically expressed at high levels in the salivary glands and constitute major components of saliva. Here we report on recurrent rearrangements [t(4;9)(q13;q31)] in AciCC that translocate active enhancer regions from the SCPP gene cluster to the region upstream of Nuclear Receptor Subfamily 4 Group A Member 3 (NR4A3) at 9q31. We show that NR4A3 is specifically upregulated in AciCCs, and that active chromatin regions and gene expression signatures in AciCCs are highly correlated with the NR4A3 transcription factor binding motif. Overexpression of NR4A3 in mouse salivary gland cells increases expression of known NR4A3 target genes and has a stimulatory functional effect on cell proliferation. We conclude that NR4A3 is upregulated through enhancer hijacking and has important oncogenic functions in AciCC.

As whole-genome sequencing for cancer genome analysis becomes a clinical tool, a full understanding of the variables affecting sequencing analysis output is required. Here using tumour-normal sample pairs from two different types of cancer, chronic lymphocytic leukaemia and medulloblastoma, we conduct a benchmarking exercise within the context of the International Cancer Genome Consortium. We compare sequencing methods, analysis pipelines and validation methods. We show that using PCR-free methods and increasing sequencing depth to ∼ 100 × shows benefits, as long as the tumour:control coverage ratio remains balanced. We observe widely varying mutation call rates and low concordance among analysis pipelines, reflecting the artefact-prone nature of the raw data and lack of standards for dealing with the artefacts. However, we show that, using the benchmark mutation set we have created, many issues are in fact easy to remedy and have an immediate positive impact on mutation detection accuracy.