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.

The genomic and transcriptional landscape of primary central nervous system lymphoma.

Primary lymphomas of the central nervous system (PCNSL) are mainly diffuse large B-cell lymphomas (DLBCLs) confined to the central nervous system (CNS). Molecular drivers of PCNSL have not been fully elucidated. Here, we profile and compare the whole-genome and transcriptome landscape of 51 CNS lymphomas (CNSL) to 39 follicular lymphoma and 36 DLBCL cases outside the CNS. We find recurrent mutations in JAK-STAT, NFkB, and B-cell receptor signaling pathways, including hallmark mutations in MYD88 L265P (67%) and CD79B (63%), and CDKN2A deletions (83%). PCNSLs exhibit significantly more focal deletions of HLA-D (6p21) locus as a potential mechanism of immune evasion. Mutational signatures correlating with DNA replication and mitosis are significantly enriched in PCNSL. TERT gene expression is significantly higher in PCNSL compared to activated B-cell (ABC)-DLBCL. Transcriptome analysis clearly distinguishes PCNSL and systemic DLBCL into distinct molecular subtypes. Epstein-Barr virus (EBV)+ CNSL cases lack recurrent mutational hotspots apart from IG and HLA-DRB loci. We show that PCNSL can be clearly distinguished from DLBCL, having distinct expression profiles, IG expression and translocation patterns, as well as specific combinations of genetic alterations.

Alternative lengthening of telomeres in childhood neuroblastoma from genome to proteome.

Telomere maintenance by telomerase activation or alternative lengthening of telomeres (ALT) is a major determinant of poor outcome in neuroblastoma. Here, we screen for ALT in primary and relapsed neuroblastomas (n = 760) and characterize its features using multi-omics profiling. ALT-positive tumors are molecularly distinct from other neuroblastoma subtypes and enriched in a population-based clinical sequencing study cohort for relapsed cases. They display reduced ATRX/DAXX complex abundance, due to either ATRX mutations (55%) or low protein expression. The heterochromatic histone mark H3K9me3 recognized by ATRX is enriched at the telomeres of ALT-positive tumors. Notably, we find a high frequency of telomeric repeat loci with a neuroblastoma ALT-specific hotspot on chr1q42.2 and loss of the adjacent chromosomal segment forming a neo-telomere. ALT-positive neuroblastomas proliferate slowly, which is reflected by a protracted clinical course of disease. Nevertheless, children with an ALT-positive neuroblastoma have dismal outcome.

The landscape of chromothripsis across adult cancer types.

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.

Genomic footprints of activated telomere maintenance mechanisms in cancer.

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.

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 .

Genomic and transcriptomic changes complement each other in the pathogenesis of sporadic Burkitt lymphoma.

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.

Integrative genomic and transcriptomic analysis of leiomyosarcoma.

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.

Epigenomic Profiling of Human CD4+ T Cells Supports a Linear Differentiation Model and Highlights Molecular Regulators of Memory Development.

The impact of epigenetics on the differentiation of memory T (Tmem) cells is poorly defined. We generated deep epigenomes comprising genome-wide profiles of DNA methylation, histone modifications, DNA accessibility, and coding and non-coding RNA expression in naive, central-, effector-, and terminally differentiated CD45RA+ CD4+ Tmem cells from blood and CD69+ Tmem cells from bone marrow (BM-Tmem). We observed a progressive and proliferation-associated global loss of DNA methylation in heterochromatic parts of the genome during Tmem cell differentiation. Furthermore, distinct gradually changing signatures in the epigenome and the transcriptome supported a linear model of memory development in circulating T cells, while tissue-resident BM-Tmem branched off with a unique epigenetic profile. Integrative analyses identified candidate master regulators of Tmem cell differentiation, including the transcription factor FOXP1. This study highlights the importance of epigenomic changes for Tmem cell biology and demonstrates the value of epigenetic data for the identification of lineage regulators.

Consequences of point mutations in melanoma-associated antigen 4 (MAGE-A4) protein: Insights from structural and biophysical studies.

The Melanoma-Associated Antigen A4 (MAGE-A4) protein is a target for cancer therapy. The function of this protein is not well understood. We report the first comprehensive study on key cancer-associated MAGE-A4 mutations and provide analysis on the consequences of these mutations on the structure, folding and stability of the protein. Based on Nuclear Magnetic Resonance and Circular Dichroism, these mutations had no significant effects on the structure and the folding of the protein. Some mutations affected the thermal stability of the protein remarkably. Native mass spectrometry of wild-type MAGE-A4 showed a broad charge state distribution suggestive of a structurally dynamic protein. Significant intensity was found in relatively low charge states, indicative of a predominantly globular form and some population in more extended states. The latter is supported by Ion Mobility measurements. The MAGE-A4 mutants exhibited similar features. These novel molecular insights shed further light on better understanding of these proteins, which are implicated in a wide range of human cancers.

The structure of a transient complex of a nonribosomal peptide synthetase and a cytochrome P450 monooxygenase.

Studying the interplay between nonribosomal peptide synthetases (NRPS), a major source of secondary metabolites, and crucial external modifying enzymes is a challenging task since the interactions involved are often transient in nature. By applying a range of synthetic inhibitor-type compounds, a stabilized complex appropriate for structural analysis was generated for such a tailoring enzyme and an NRPS domain. The complex studied comprises an NRPS peptidyl carrier protein (PCP) domain bound to the Cytochrome P450 enzyme that is crucial for the provision of ß-hydroxylated amino acid precursors in the biosynthesis of the cyclic depsipeptide skyllamycin. The structure reveals that complex formation is governed by hydrophobic interactions, the presence of which can be controlled through minor alterations in PCP structure that enable selectivity amongst multiple highly similar PCP domains.

Blind spots of quantitative RNA-seq: the limits for assessing abundance, differential expression, and isoform switching.

BACKGROUND: RNA-seq is now widely used to quantitatively assess gene expression, expression differences and isoform switching, and promises to deliver results for the entire transcriptome. However, whether the transcriptional state of a gene can be captured accurately depends critically on library preparation, read alignment, expression estimation and the tests for differential expression and isoform switching. There are comparisons available for the individual steps but there is not yet a systematic investigation which specific genes are impacted by biases throughout the entire analysis workflow. It is especially unclear whether for a given gene, with current methods and protocols, expression changes and isoform switches can be detected. RESULTS: For the human genes, we report their detectability under various conditions using different approaches. Overall, we find that the input material has the biggest influence and may, depending on the protocol and RNA degradation, exhibit already strong length-dependent over- and underrepresentation of transcripts. The alignment step aligns for 50% of the isoforms up to 99% of the reads correctly; only in the presence of transcript modifications mainly short isoforms will have a low alignment rate. In our dataset, we found that, depending on the aligner and the input material used, the expression estimation of up to 93% of the genes being accurate within a factor of two; with the deviations being due to ambiguous alignments. Detection of differential expression using a negative-binomial count model works reliably for our simulated data but is dependent on the count accuracy. Interestingly, using the fold-change instead of the p-value as a score for differential expression yields the same performance in the situation of three replicates and the true change being two-fold. Isoform switching is harder to detect and for at least 109 genes the isoform differences evade detection independent of the method used. CONCLUSIONS: RNA-seq is a reliable tool but the repetitive nature of the human genome makes the origin of the reads ambiguous and limits the detectability for certain genes. RNA-seq does not equally well represent isoforms independent of their size which may range from ~200nt to ~100'000nt. Researchers are advised to verify that their target genes do not have extreme properties with respect to repeated regions, GC content, and isoform length and complexity.