Single-cell division tracing and transcriptomics reveal cell types and differentiation paths in the regenerating lung.

Understanding the molecular and cellular processes involved in lung epithelial regeneration may fuel the development of therapeutic approaches for lung diseases. We combine mouse models allowing diphtheria toxin-mediated damage of specific epithelial cell types and parallel GFP-labeling of functionally dividing cells with single-cell transcriptomics to characterize the regeneration of the distal lung. We uncover cell types, including Krt13+ basal and Krt15+ club cells, detect an intermediate cell state between basal and goblet cells, reveal goblet cells as actively dividing progenitor cells, and provide evidence that adventitial fibroblasts act as supporting cells in epithelial regeneration. We also show that diphtheria toxin-expressing cells can persist in the lung, express specific inflammatory factors, and transcriptionally resemble a previously undescribed population in the lungs of COVID-19 patients. Our study provides a comprehensive single-cell atlas of the distal lung that characterizes early transcriptional and cellular responses to concise epithelial injury, encompassing proliferation, differentiation, and cell-to-cell interactions.

Multi-omic and functional analysis for classification and treatment of sarcomas with FUS-TFCP2 or EWSR1-TFCP2 fusions.

Linking clinical multi-omics with mechanistic studies may improve the understanding of rare cancers. We leverage two precision oncology programs to investigate rhabdomyosarcoma with FUS/EWSR1-TFCP2 fusions, an orphan malignancy without effective therapies. All tumors exhibit outlier ALK expression, partly accompanied by intragenic deletions and aberrant splicing resulting in ALK variants that are oncogenic and sensitive to ALK inhibitors. Additionally, recurrent CKDN2A/MTAP co-deletions provide a rationale for PRMT5-targeted therapies. Functional studies show that FUS-TFCP2 blocks myogenic differentiation, induces transcription of ALK and truncated TERT, and inhibits DNA repair. Unlike other fusion-driven sarcomas, TFCP2-rearranged tumors exhibit genomic instability and signs of defective homologous recombination. DNA methylation profiling demonstrates a close relationship with undifferentiated sarcomas. In two patients, sarcoma was preceded by benign lesions carrying FUS-TFCP2, indicating stepwise sarcomagenesis. This study illustrates the potential of linking precision oncology with preclinical research to gain insight into the classification, pathogenesis, and therapeutic vulnerabilities of rare cancers.

Single-cell RNA sequencing of mouse lower respiratory tract epithelial cells: A meta-analysis.

The respiratory system is a vital component of our body, essential for both oxygen uptake and immune defense. Knowledge of cellular composition and function in different parts of the respiratory tract provides the basis for a better understanding of the pathological processes involved in various diseases such as chronic respiratory diseases and cancer. Single-cell RNA sequencing (scRNA-seq) is a proficient approach for the identification and transcriptional characterization of cellular phenotypes. Although the mouse is an essential tool for the study of lung development, regeneration, and disease, a scRNA-seq mouse atlas of the lung in which all epithelial cell types are included and annotated systematically is lacking. Here, we established a single-cell transcriptome landscape of the mouse lower respiratory tract by performing a meta-analysis of seven different studies in which mouse lungs and trachea were analyzed by droplet and/or plate-based scRNA-seq technologies. We provide information on the best markers for each epithelial cell type, propose surface markers for the isolation of viable cells, harmonized the annotation of cell types, and compare the mouse single-cell transcriptomes with human scRNA-seq data of the lung.

Transcriptome profiling for precision cancer medicine using shallow nanopore cDNA sequencing.

Transcriptome profiling is a mainstay of translational cancer research and is increasingly finding its way into precision oncology. While bulk RNA sequencing (RNA-seq) is widely available, high investment costs and long data return time are limiting factors for clinical applications. We investigated a portable nanopore long-read sequencing device (MinION, Oxford Nanopore Technologies) for transcriptome profiling of tumors. In particular, we investigated the impact of lower coverage than that of larger sequencing devices by comparing shallow nanopore RNA-seq data with short-read RNA-seq data generated using reversible dye terminator technology (Illumina) for ten samples representing four cancer types. Coupled with ShaNTi (Shallow Nanopore sequencing for Transcriptomics), a newly developed data processing pipeline, a turnaround time of five days was achieved. The correlation of normalized gene-level counts between nanopore and Illumina RNA-seq was high for MinION but not for very low-throughput Flongle flow cells (r = 0.89 and r = 0.24, respectively). A cost-saving approach based on multiplexing of four samples per MinION flow cell maintained a high correlation with Illumina data (r = 0.56-0.86). In addition, we compared the utility of nanopore and Illumina RNA-seq data for analysis tools commonly applied in translational oncology: (1) Shallow nanopore and Illumina RNA-seq were equally useful for inferring signaling pathway activities with PROGENy. (2) Highly expressed genes encoding kinases targeted by clinically approved small-molecule inhibitors were reliably identified by shallow nanopore RNA-seq. (3) In tumor microenvironment composition analysis, quanTIseq performed better than CIBERSORT, likely due to higher average expression of the gene set used for deconvolution. (4) Shallow nanopore RNA-seq was successfully applied to detect fusion genes using the JAFFAL pipeline. These findings suggest that shallow nanopore RNA-seq enables rapid and biologically meaningful transcriptome profiling of tumors, and warrants further exploration in precision cancer medicine studies.

BTBBCL6 dimers as building blocks for reversible drug-induced protein oligomerization.

Here, we characterize the BTB domain of the transcription factor BCL6 (BTBBCL6) as a small-molecule-controlled, reversible oligomerization switch, which oligomerizes upon BI-3802 treatment and de-oligomerizes upon addition of BI-3812. We show that the magnitude of oligomerization can be controlled in vitro by BI-3802 concentration and exposure time. In cellular models, exposure to BI-3802/BI-3812 can drive multiple cycles of foci formation consisting of BTBBCL6 fused to EGFP, which are not degraded due to the lack of a degron. We generated an epidermal growth factor receptor (EGFR)-BTBBCL6 fusion. Treatment with BI-3802, as an ON switch, induced EGFR-BTBBCL6 phosphorylation and activation of downstream effectors, which could in part be reversed by the addition of BI-3812, as an OFF switch. Finally, BI-3802-induced oligomerization of the EGFR-BTBBCL6 fusion enhanced proliferation of an EGF-dependent cell line in absence of EGF. These results demonstrate the successful application of small-molecule-induced, reversible oligomerization as a switch for synthetic biology.

Fusion protein-driven IGF-IR/PI3K/AKT signals deregulate Hippo pathway promoting oncogenic cooperation of YAP1 and FUS-DDIT3 in myxoid liposarcoma.

Myxoid liposarcoma (MLS) represents a common subtype of liposarcoma molecularly characterized by a recurrent chromosomal translocation that generates a chimeric FUS-DDIT3 fusion gene. The FUS-DDIT3 oncoprotein has been shown to be crucial in MLS pathogenesis. Acting as a transcriptional dysregulator, FUS-DDIT3 stimulates proliferation and interferes with adipogenic differentiation. As the fusion protein represents a therapeutically challenging target, a profound understanding of MLS biology is elementary to uncover FUS-DDIT3-dependent molecular vulnerabilities. Recently, a specific reliance on the Hippo pathway effector and transcriptional co-regulator YAP1 was detected in MLS; however, details on the molecular mechanism of FUS-DDIT3-dependent YAP1 activation, and YAP1´s precise mode of action remain unclear. In elaborate in vitro studies, employing RNA interference-based approaches, small-molecule inhibitors, and stimulation experiments with IGF-II, we show that FUS-DDIT3-driven IGF-IR/PI3K/AKT signaling promotes stability and nuclear accumulation of YAP1 via deregulation of the Hippo pathway. Co-immunoprecipitation and proximity ligation assays revealed nuclear co-localization of FUS-DDIT3 and YAP1/TEAD in FUS-DDIT3-expressing mesenchymal stem cells and MLS cell lines. Transcriptome sequencing of MLS cells demonstrated that FUS-DDIT3 and YAP1 co-regulate oncogenic gene signatures related to proliferation, cell cycle progression, apoptosis, and adipogenesis. In adipogenic differentiation assays, we show that YAP1 critically contributes to FUS-DDIT3-mediated adipogenic differentiation arrest. Taken together, our study provides mechanistic insights into a complex FUS-DDIT3-driven network involving IGF-IR/PI3K/AKT signals acting on Hippo/YAP1, and uncovers substantial cooperative effects of YAP1 and FUS-DDIT3 in the pathogenesis of MLS.

Accurate and efficient detection of gene fusions from RNA sequencing data.

The identification of gene fusions from RNA sequencing data is a routine task in cancer research and precision oncology. However, despite the availability of many computational tools, fusion detection remains challenging. Existing methods suffer from poor prediction accuracy and are computationally demanding. We developed Arriba, a novel fusion detection algorithm with high sensitivity and short runtime. When applied to a large collection of published pancreatic cancer samples (n = 803), Arriba identified a variety of driver fusions, many of which affected druggable proteins, including ALK, BRAF, FGFR2, NRG1, NTRK1, NTRK3, RET, and ROS1. The fusions were significantly associated with KRAS wild-type tumors and involved proteins stimulating the MAPK signaling pathway, suggesting that they substitute for activating mutations in KRAS In addition, we confirmed the transforming potential of two novel fusions, RRBP1-RAF1 and RASGRP1-ATP1A1, in cellular assays. These results show Arriba's utility in both basic cancer research and clinical translation.

Small-molecule-induced polymerization triggers degradation of BCL6.

Effective and sustained inhibition of non-enzymatic oncogenic driver proteins is a major pharmacological challenge. The clinical success of thalidomide analogues demonstrates the therapeutic efficacy of drug-induced degradation of transcription factors and other cancer targets1-3, but a substantial subset of proteins are resistant to targeted degradation using existing approaches4,5. Here we report an alternative mechanism of targeted protein degradation, in which a small molecule induces the highly specific, reversible polymerization of a target protein, followed by its sequestration into cellular foci and subsequent degradation. BI-3802 is a small molecule that binds to the Broad-complex, Tramtrack and Bric-à-brac (BTB) domain of the oncogenic transcription factor B cell lymphoma 6 (BCL6) and leads to the proteasomal degradation of BCL66. We use cryo-electron microscopy to reveal how the solvent-exposed moiety of a BCL6-binding molecule contributes to a composite ligand-protein surface that engages BCL6 homodimers to form a supramolecular structure. Drug-induced formation of BCL6 filaments facilitates ubiquitination by the SIAH1 E3 ubiquitin ligase. Our findings demonstrate that a small molecule such as BI-3802 can induce polymerization coupled to highly specific protein degradation, which in the case of BCL6 leads to increased pharmacological activity compared to the effects induced by other BCL6 inhibitors. These findings open new avenues for the development of therapeutic agents and synthetic biology.

Hematopoietic stem and progenitor cell-restricted Cdx2 expression induces transformation to myelodysplasia and acute leukemia.

The caudal-related homeobox transcription factor CDX2 is expressed in leukemic cells but not during normal blood formation. Retroviral overexpression of Cdx2 induces AML in mice, however the developmental stage at which CDX2 exerts its effect is unknown. We developed a conditionally inducible Cdx2 mouse model to determine the effects of in vivo, inducible Cdx2 expression in hematopoietic stem and progenitor cells (HSPCs). Cdx2-transgenic mice develop myelodysplastic syndrome with progression to acute leukemia associated with acquisition of additional driver mutations. Cdx2-expressing HSPCs demonstrate enrichment of hematopoietic-specific enhancers associated with pro-differentiation transcription factors. Furthermore, treatment of Cdx2 AML with azacitidine decreases leukemic burden. Extended scheduling of low-dose azacitidine shows greater efficacy in comparison to intermittent higher-dose azacitidine, linked to more specific epigenetic modulation. Conditional Cdx2 expression in HSPCs is an inducible model of de novo leukemic transformation and can be used to optimize treatment in high-risk AML.

Requirement for LIM kinases in acute myeloid leukemia.

Acute myeloid leukemia (AML) is an aggressive disease for which only few targeted therapies are available. Using high-throughput RNA interference (RNAi) screening in AML cell lines, we identified LIM kinase 1 (LIMK1) as a potential novel target for AML treatment. High LIMK1 expression was significantly correlated with shorter survival of AML patients and coincided with FLT3 mutations, KMT2A rearrangements, and elevated HOX gene expression. RNAi- and CRISPR-Cas9-mediated suppression as well as pharmacologic inhibition of LIMK1 and its close homolog LIMK2 reduced colony formation and decreased proliferation due to slowed cell-cycle progression of KMT2A-rearranged AML cell lines and patient-derived xenograft (PDX) samples. This was accompanied by morphologic changes indicative of myeloid differentiation. Transcriptome analysis showed upregulation of several tumor suppressor genes as well as downregulation of HOXA9 targets and mitosis-associated genes in response to LIMK1 suppression, providing a potential mechanistic basis for the anti-leukemic phenotype. Finally, we observed a reciprocal regulation between LIM kinases (LIMK) and CDK6, a kinase known to be involved in the differentiation block of KMT2A-rearranged AML, and addition of the CDK6 inhibitor palbociclib further enhanced the anti-proliferative effect of LIMK inhibition. Together, these data suggest that LIMK are promising targets for AML therapy.

The CDK inhibitor CR8 acts as a molecular glue degrader that depletes cyclin K.

Molecular glue compounds induce protein-protein interactions that, in the context of a ubiquitin ligase, lead to protein degradation1. Unlike traditional enzyme inhibitors, these molecular glue degraders act substoichiometrically to catalyse the rapid depletion of previously inaccessible targets2. They are clinically effective and highly sought-after, but have thus far only been discovered serendipitously. Here, through systematically mining databases for correlations between the cytotoxicity of 4,518 clinical and preclinical small molecules and the expression levels of E3 ligase components across hundreds of human cancer cell lines3-5, we identify CR8-a cyclin-dependent kinase (CDK) inhibitor6-as a compound that acts as a molecular glue degrader. The CDK-bound form of CR8 has a solvent-exposed pyridyl moiety that induces the formation of a complex between CDK12-cyclin K and the CUL4 adaptor protein DDB1, bypassing the requirement for a substrate receptor and presenting cyclin K for ubiquitination and degradation. Our studies demonstrate that chemical alteration of surface-exposed moieties can confer gain-of-function glue properties to an inhibitor, and we propose this as a broader strategy through which target-binding molecules could be converted into molecular glues.

TelomereHunter - in silico estimation of telomere content and composition from cancer genomes.

BACKGROUND: Establishment of telomere maintenance mechanisms is a universal step in tumor development to achieve replicative immortality. These processes leave molecular footprints in cancer genomes in the form of altered telomere content and aberrations in telomere composition. To retrieve these telomere characteristics from high-throughput sequencing data the available computational approaches need to be extended and optimized to fully exploit the information provided by large scale cancer genome data sets. RESULTS: We here present TelomereHunter, a software for the detailed characterization of telomere maintenance mechanism footprints in the genome. The tool is implemented for the analysis of large cancer genome cohorts and provides a variety of diagnostic diagrams as well as machine-readable output for subsequent analysis. A novel key feature is the extraction of singleton telomere variant repeats, which improves the identification and subclassification of the alternative lengthening of telomeres phenotype. We find that whole genome sequencing-derived telomere content estimates strongly correlate with telomere qPCR measurements (r = 0.94). For the first time, we determine the correlation of in silico telomere content quantification from whole genome sequencing and whole genome bisulfite sequencing data derived from the same tumor sample (r = 0.78). An analogous comparison of whole exome sequencing data and whole genome sequencing data measured slightly lower correlation (r = 0.79). However, this is considerably improved by normalization with matched controls (r = 0.91). CONCLUSIONS: TelomereHunter provides new functionality for the analysis of the footprints of telomere maintenance mechanisms in cancer genomes. Besides whole genome sequencing, whole exome sequencing and whole genome bisulfite sequencing are suited for in silico telomere content quantification, especially if matched control samples are available. The software runs under a GPL license and is available at https://www.dkfz.de/en/applied-bioinformatics/telomerehunter/telomerehunter.html .

Defective homologous recombination DNA repair as therapeutic target in advanced chordoma.

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.

Requirement for YAP1 signaling in myxoid liposarcoma.

Myxoid liposarcomas (MLS), malignant tumors of adipocyte origin, are driven by the FUS-DDIT3 fusion gene encoding an aberrant transcription factor. The mechanisms whereby FUS-DDIT3 mediates sarcomagenesis are incompletely understood, and strategies to selectively target MLS cells remain elusive. Here we show, using an unbiased functional genomic approach, that FUS-DDIT3-expressing mesenchymal stem cells and MLS cell lines are dependent on YAP1, a transcriptional co-activator and central effector of the Hippo pathway involved in tissue growth and tumorigenesis, and that increased YAP1 activity is a hallmark of human MLS Mechanistically, FUS-DDIT3 promotes YAP1 expression, nuclear localization, and transcriptional activity and physically associates with YAP1 in the nucleus of MLS cells. Pharmacologic inhibition of YAP1 activity impairs the growth of MLS cells in vitro and in vivo These findings identify overactive YAP1 signaling as unifying feature of MLS development that could represent a novel target for therapeutic intervention.

Exploiting rare driver mutations for precision cancer medicine.

Catalyzed by the ability to develop precision therapies targeting the unique genetic changes that drive individual tumors, sequencing patients' tumor genomes is an increasingly common practice in oncology. In most cancer types, however, a limited number of common mutations are accompanied by a plethora of low-frequency variants whose functional consequences and clinical actionability are often unknown. We here illustrate that this 'long tail' of infrequent molecular alterations includes oncogenic drivers of biological significance that can be the genetic basis of extraordinary responses to systemic cancer therapies. Furthermore, we review current strategies to identify, prioritize, and experimentally validate novel long-tail driver mutations, efforts that will likely provide new insights into the clinically actionable genome and improve outcomes for patients.

NCAM1 (CD56) promotes leukemogenesis and confers drug resistance in AML.

Neural cell adhesion molecule 1 (NCAM1; CD56) is expressed in up to 20% of acute myeloid leukemia (AML) patients. NCAM1 is widely used as a marker of minimal residual disease; however, the biological function of NCAM1 in AML remains elusive. In this study, we investigated the impact of NCAM1 expression on leukemogenesis, drug resistance, and its role as a biomarker to guide therapy. Beside t(8;21) leukemia, NCAM1 expression was found in most molecular AML subgroups at highly heterogeneous expression levels. Using complementary genetic strategies, we demonstrated an essential role of NCAM1 in the regulation of cell survival and stress resistance. Perturbation of NCAM1 induced cell death or differentiation and sensitized leukemic blasts toward genotoxic agents in vitro and in vivo. Furthermore, Ncam1 was highly expressed in leukemic progenitor cells in a murine leukemia model, and genetic depletion of Ncam1 prolonged disease latency and significantly reduced leukemia-initiating cells upon serial transplantation. To further analyze the mechanism of the NCAM1-associated phenotype, we performed phosphoproteomics and transcriptomics in different AML cell lines. NCAM1 expression strongly associated with constitutive activation of the MAPK-signaling pathway, regulation of apoptosis, or glycolysis. Pharmacological inhibition of MEK1/2 specifically inhibited proliferation and sensitized NCAM1+ AML cells to chemotherapy. In summary, our data demonstrate that aberrant expression of NCAM1 is involved in the maintenance of leukemic stem cells and confers stress resistance, likely due to activation of the MAPK pathway. Targeting MEK1/2 sensitizes AML blasts to genotoxic agents, indicating a role for NCAM1 as a biomarker to guide AML treatment.

Cullin 5 is a novel candidate tumor suppressor in renal cell carcinoma involved in the maintenance of genome stability.

Clear cell renal cell carcinoma (ccRCC) is intimately associated with defects in ubiquitin-mediated protein degradation. Herein, we report that deficiency in the E3 ligase subunit cullin 5 (CUL5) promotes chromosomal instability and is an independent negative prognostic factor in ccRCC. CUL5 was initially identified in an RNA interference screen as a novel regulator of centrosome duplication control. We found that depletion of CUL5 rapidly promotes centriole overduplication and mitotic errors. Downregulation of CUL5 also caused an increase of DNA damage that was found to involve impaired DNA double-strand break repair. Using immunohistochemistry, CUL5 protein expression was found to be below detection level in the majority of RCCs. A re-analysis of the TCGA ccRCC cohort showed that a reduced CUL5 gene expression or CUL5 deletion were associated with a significantly worse overall patient survival. In conclusion, our results indicate that CUL5 functions as a novel tumor suppressor with prognostic relevance in ccRCC and is critically involved in the maintenance of genome stability.